| 2 | 1/1 | 返回列表 |
| 查看: 1137 | 回復(fù): 1 | ||||
[交流]
【求助】siesta中PDOS文件的處理
|
| 請(qǐng)教高手,如題,在siesta中怎么處理PDOS文件! |
siesta問題 |
» 搶金幣啦!回帖就可以得到:
+1/462
+1/84
+1/83
+1/74
+1/63
+1/39
+1/34
+1/18
+1/17
+1/17
+1/16
+1/15
+1/14
+1/9
+1/7
+1/6
+1/5
+1/4
+1/1
+1/1
|
百度來(lái)的哈~~ 一: pdosxm 我主要參考的http://fisica.ehu.es/ag/siesta-extra/issues.html 使得siesta輸出的PDOS文件中沒有*。 另外對(duì)于nspin=2的情況,修改了pdosxml的代碼,使得其能用。 (1) 在arch.make中加入 FPPFLAGS= -DWXML_INIT_FIX $make clean , $make 這樣siesta輸出的PDOS文件中就沒有*。 (2)推薦使用/Util/pdosxml,來(lái)讀取PDOS文件,并轉(zhuǎn)化成gnuplot能處理的形式 使用參見/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ù)分開的PDOS了。 上面修改好的siesta放在 cluster4 /public/soft/siesta/siesta_Util_pdosxml 二 : readpdos 我寫的一個(gè)小代碼,放在 210.45.121.200 ~/share/Seminar_backup/siesta_tutorial/soft_siesta/readpdos 可以對(duì)有*號(hào)的PDOS文件都能處理,但是每次體系,都要重新修改代碼 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 316求調(diào)劑 +5 | 梁茜雯 2026-03-19 | 5/250 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|