| 2 | 1/1 | 返回列表 |
| 查看: 1136 | 回復(fù): 1 | ||||
maomao-666木蟲(chóng) (小有名氣)
|
[交流]
【求助】siesta中PDOS文件的處理
|
| 請(qǐng)教高手,如題,在siesta中怎么處理PDOS文件! |
siesta問(wèn)題 |
» 搶金幣啦!回帖就可以得到:
+1/185
+1/95
+1/83
+1/79
+1/74
+1/53
+1/36
+1/35
+1/34
+1/34
+1/31
+1/31
+1/18
+1/7
+1/6
+5/5
+1/5
+3/2
+1/1
+1/1
鐵蟲(chóng) (初入文壇)
|
百度來(lái)的哈~~ 一: pdosxm 我主要參考的http://fisica.ehu.es/ag/siesta-extra/issues.html 使得siesta輸出的PDOS文件中沒(méi)有*。 另外對(duì)于nspin=2的情況,修改了pdosxml的代碼,使得其能用。 (1) 在arch.make中加入 FPPFLAGS= -DWXML_INIT_FIX $make clean , $make 這樣siesta輸出的PDOS文件中就沒(méi)有*。 (2)推薦使用/Util/pdosxml,來(lái)讀取PDOS文件,并轉(zhuǎn)化成gnuplot能處理的形式 使用參見(jiàn)/Util/pdosxml/README 主要步驟如下: 1 在m_orbital_chooser.f90 中選好你要的軌道 2 make 3 ./pdos h2o_dos.PDOS >out 采用原來(lái)的代碼,計(jì)算nspin=2時(shí)會(huì)出錯(cuò),http://fisica.ehu.es/ag/siesta-extra/issues.html給出的 m_pdos.f90文件 其實(shí)和原來(lái)的一樣。我們要自己改一下代碼 在 m_pdos.f90 大概150多行的地方: if (spin_polarized) then dos1 = dos1 + data(1:ndata:2) dos2 = dos2 + data(2:ndata:2) else dos1 = dos1 + data endif || 替換成 || || \ / if (spin_polarized) then !shanghui add here to make nspin=2's pdos do i=1,ndata/2 dos1(i) = dos1(i) + data(2*i-1) dos2(i) = dos2(i) - data(2*i) enddo !shanghui end add here else dos1 = dos1 + data endif 這樣你就能做出正負(fù)分開(kāi)的PDOS了。 上面修改好的siesta放在 cluster4 /public/soft/siesta/siesta_Util_pdosxml 二 : readpdos 我寫(xiě)的一個(gè)小代碼,放在 210.45.121.200 ~/share/Seminar_backup/siesta_tutorial/soft_siesta/readpdos 可以對(duì)有*號(hào)的PDOS文件都能處理,但是每次體系,都要重新修改代碼 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | Mqqqqqq 2026-03-19 | 5/250 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 265求調(diào)劑 +8 | 梁梁校校 2026-03-17 | 8/400 |
|
|
[考研] 286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 288求調(diào)劑 +15 | 于海海海海 2026-03-19 | 15/750 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|