ustbwgs

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 14 (小學(xué)生)
金幣: 1449.6
帖子: 720
在線: 165.5小時(shí)
蟲(chóng)號(hào): 1141993

[交流] 【求助】DOSCAR態(tài)密度里面的數(shù)據(jù)如何分離

我想自己編個(gè)小程序把態(tài)密度里面的數(shù)據(jù)分離出來(lái)，然后再做圖
高手們指點(diǎn)一下，這么多數(shù)據(jù)分別代表什么意思呢？？？
剛開(kāi)始是文件頭
接著有3列，后來(lái)變成10列了，是什么數(shù)據(jù)呢，我是用ICHARG=11數(shù)完之后的結(jié)果》》》

[ Last edited by ustbwgs on 2010-12-5 at 09:59 ]

回復(fù)此樓

» 猜你喜歡

262求調(diào)劑已經(jīng)有3人回復(fù)
303求調(diào)劑已經(jīng)有5人回復(fù)
300求調(diào)劑，材料科學(xué)英一數(shù)二已經(jīng)有17人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有21人回復(fù)
307求調(diào)劑已經(jīng)有15人回復(fù)
286分調(diào)劑已經(jīng)有9人回復(fù)
求調(diào)劑323材料與化工已經(jīng)有8人回復(fù)
322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。已經(jīng)有5人回復(fù)
0703化學(xué)321分求調(diào)劑已經(jīng)有6人回復(fù)
292求調(diào)劑已經(jīng)有10人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

同樣情況下，不同的導(dǎo)出數(shù)據(jù)方法態(tài)密度怎么不一樣已經(jīng)有16人回復(fù)
用originpro 8怎樣畫(huà)態(tài)密度的圖急急急已經(jīng)有19人回復(fù)
DOSCAR畫(huà)態(tài)密度已經(jīng)有16人回復(fù)
態(tài)密度圖跟文獻(xiàn)里面的怎么都不一樣，求救已經(jīng)有36人回復(fù)
【求助】總態(tài)密度與分波態(tài)密度已經(jīng)有9人回復(fù)
【求助】如何用DOSCAR畫(huà)態(tài)密度呢已經(jīng)有12人回復(fù)
【求助】畫(huà)態(tài)密度圖已經(jīng)有7人回復(fù)
【求助】態(tài)密度怎么算？老是出現(xiàn)這種問(wèn)題已經(jīng)有7人回復(fù)
【分享】態(tài)密度處理程序已經(jīng)有17人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2010-12-05 01:15:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

stractor

金蟲(chóng) (著名寫(xiě)手)

1ST強(qiáng)帖: 6
應(yīng)助: 91 (初中生)
金幣: 1487.6
帖子: 2224
在線: 426.8小時(shí)
蟲(chóng)號(hào): 394087

★ ★ ★
ustbwgs(金幣+1):謝謝參與
sunyang1988(金幣+2):辛苦了 2010-12-05 11:22:27

前面的三列是能量，總態(tài)密度和態(tài)密度的積分，即電子數(shù)。
后面的10列是能量和每個(gè)原子各原子軌道的投影態(tài)密度，分別是
s py pz px dxy dyz dz2 dxz dx2
原子順序按POSCAR的順序給出。

[ Last edited by stractor on 2010-12-5 at 20:58 ]

贊一下(2人)

回復(fù)此樓

2樓2010-12-05 11:07:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ustbwgs

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 14 (小學(xué)生)
金幣: 1449.6
帖子: 720
在線: 165.5小時(shí)
蟲(chóng)號(hào): 1141993

引用回帖:

Originally posted by stractor at 2010-12-05 11:07:35:
前面的三列是能量，總態(tài)密度和態(tài)密度的積分，即電子數(shù)。
后面的10列是能量和每個(gè)原子各原子軌道的投影態(tài)密度，分別是
s py pz px dxy dyz dz2 dxz ...

非常感謝你的幫助，后面的10列vasp官方手冊(cè)上說(shuō)是：：：能量，s-DOS(total) S-DSO(mx) S-DOS(my) S-DOS(mz) 然后是P，也是按這個(gè)順序，和你說(shuō)得有點(diǎn)不一樣啊，
不過(guò)按書(shū)上的一直排下去，會(huì)不會(huì)是P完了再D，會(huì)不會(huì)原子序數(shù)大的原子會(huì)有十幾列啊，還是都是10列呢？？？？

贊一下

回復(fù)此樓

3樓2010-12-05 13:52:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunyang1988

至尊木蟲(chóng) (著名寫(xiě)手)

1ST強(qiáng)帖: 1
應(yīng)助: 24 (小學(xué)生)
貴賓: 0.984
金幣: 9397.7
帖子: 1744
在線: 1255.8小時(shí)
蟲(chóng)號(hào): 965316

★ ★ ★
ustbwgs(金幣+1):謝謝參與
ustbwgs(金幣+3): 2010-12-05 15:29:17
youzhizhe(金幣+2):謝謝交流。 2010-12-05 22:49:26

S-DOS沒(méi)有mxmymz，P-DOS才有的，是按照spd這樣排的，沒(méi)有算過(guò)包含f電子的不知道具體情況。
還有要注意的是ISPIN=2時(shí)數(shù)據(jù)會(huì)多一倍，每一個(gè)組分都包含自旋向上與自旋向下的成分

贊一下(2人)

回復(fù)此樓

4樓2010-12-05 14:43:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ustbwgs

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 14 (小學(xué)生)
金幣: 1449.6
帖子: 720
在線: 165.5小時(shí)
蟲(chóng)號(hào): 1141993

引用回帖:

Originally posted by sunyang1988 at 2010-12-05 14:43:49:
S-DOS沒(méi)有mxmymz，P-DOS才有的，是按照spd這樣排的，沒(méi)有算過(guò)包含f電子的不知道具體情況。
還有要注意的是ISPIN=2時(shí)數(shù)據(jù)會(huì)多一倍，每一個(gè)組分都包含自旋向上與自旋向下的成分

謝謝，我也是感覺(jué)S就不可能有啊，還有就是在算態(tài)密度時(shí)，我們一般使LORBIT等于多少呢，取默認(rèn)的不行嗎？有一次我在INCAR中沒(méi)設(shè)置這個(gè)參數(shù)，DOSCAR就沒(méi)輸出了，只有個(gè)文件頭。

贊一下

回復(fù)此樓

5樓2010-12-05 15:27:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

stractor

金蟲(chóng) (著名寫(xiě)手)

1ST強(qiáng)帖: 6
應(yīng)助: 91 (初中生)
金幣: 1487.6
帖子: 2224
在線: 426.8小時(shí)
蟲(chóng)號(hào): 394087

★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
sunyang1988(金幣+1):謝謝交流 2010-12-05 20:58:32

最好設(shè)置LORBIT=11吧，這樣打印的態(tài)密度信息全面一些。
如果不設(shè)置LORBIT=11而采用默認(rèn)值，就只有總的態(tài)密度，也就是沒(méi)有后面各原子的10列數(shù)據(jù)。

[ Last edited by stractor on 2010-12-5 at 20:59 ]

贊一下(3人)

回復(fù)此樓

6樓2010-12-05 20:55:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 262求調(diào)劑 +3	ZZ..000 2026-03-30	3/150	2026-03-30 15:24 by michael2011
[考研] 一志愿鄭大材料工程290求調(diào)劑 +4	Youth_ 2026-03-30	4/200	2026-03-30 13:09 by JourneyLucky
[考研] 化學(xué)308分調(diào)劑 +7	你好明天你好 2026-03-30	8/400	2026-03-30 13:02 by 你好明天你好
[考研] 332求調(diào)劑 +12	032500 2026-03-25	12/600	2026-03-30 12:16 by 學(xué)員8dgXkO
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10	景隨風(fēng) 2026-03-29	15/750	2026-03-30 11:21 by limeifeng
[考研] 327求調(diào)劑 +4	小卡不卡. 2026-03-29	4/200	2026-03-30 06:15 by wxiongid
[考研] 291求調(diào)劑 +7	Y-cap 2026-03-29	8/400	2026-03-29 19:53 by klasasda
[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 求收留 +5	1943443204 2026-03-28	5/250	2026-03-29 12:11 by 無(wú)際的草原
[考研] 356求調(diào)劑 +4	gysy?s?a 2026-03-28	4/200	2026-03-29 10:32 by 唐沐兒
[考研] 315求調(diào)劑 +4	akie... 2026-03-28	5/250	2026-03-28 21:05 by zhq0425
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 325求調(diào)劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚(yú)的貓
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-26 19:57 by nihaoar
[考研] 生物學(xué) 296 求調(diào)劑 +4	朵朵- 2026-03-26	6/300	2026-03-26 19:01 by 不吃魚(yú)的貓
[考研] 281求調(diào)劑 +6	Koxui 2026-03-24	7/350	2026-03-26 15:37 by 無(wú)際的草原
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣(mài)魚(yú) 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 各位老師您好：本人初試372分 +5	jj涌77 2026-03-25	6/300	2026-03-25 14:15 by mapenggao
[考研] 材料考研調(diào)劑生 +3	黃粱一夢(mèng)千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike