 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4167 | 回復(fù): 17 | ||||||
| 【獎勵】 本帖被評價12次，作者nkleof增加金幣 10.8 個 | ||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||||
[資源]
【分享】態(tài)密度處理程序
|
||||||
|
之前發(fā)過一個簡單版本的態(tài)密度處理程序，這次更新一下。 split_dos.ksh是根據(jù)JMS NRL的腳本改的，其它程序是我自己寫的。 用法：只能在linux系統(tǒng)下，鍵入各個程序名稱運行就可以，如果遇到平臺問題，sumdos程序無法執(zhí)行，可以用 gcc -lm sum.c 編譯源文件產(chǎn)生可執(zhí)行程序。 運行split_dos.ksh會讀取DOSCAR，并生成總的態(tài)密度DOS0和各個原子的各軌道態(tài)密度DOS1~DOSn。（這里注意在計算DOS時INCAR文件里的NEDOS不要設(shè)置為>9999，五位數(shù)的NEDOS值會造成錯誤） 然后運行sumdos可以產(chǎn)生按照特定原子特定軌道模式生成的態(tài)密度和（注意：適用于ISPIN=2的計算結(jié)果，即DOSCAR里面共有有20列數(shù)據(jù)）。 例： sumdos name -t 1 2 3 sumdos是程序名稱，name是一個自定義字符串，會加在生成的文件名里，可以作為說明（但是長度只有8個字符），-t表示生成各個軌道的態(tài)密度和，（還有兩種模式: -n 表示不加和各軌道態(tài)密度，仍舊生成從s到d各個軌道的態(tài)密度； -d 表示把p軌道，eg軌道，t2g軌道各自加和。��！仍然注意，和上次一樣，這里的eg指的是第14、15和18、19列數(shù)據(jù)，t2g是剩下的3列數(shù)據(jù)，如果有各個d軌道不是按照這樣排列的情形則生成的就不是eg和t2g），后面的數(shù)字表示要加和的原子標(biāo)號。上面的例子就是把第1、2、3號原子的所有軌道的態(tài)密度全部加起來。程序支持采用 n~m 形式表示要處理的第n到m個標(biāo)號連續(xù)的原子，例如要處理13到26號原子時，不需要一個一個寫出原子標(biāo)號，只需要寫 13~26 就可以。 附送一個將CONTCAR文件轉(zhuǎn)換為可導(dǎo)入Materials Studio的cif文件的程序，使用很簡單，只需要在linux下，vasp任務(wù)的計算文件夾下，運行exe_Ctocif，程序會自動尋找CONTCAR和POTCAR，并生成cif文件。另外也支持在只有POSCAR文件時轉(zhuǎn)換為cif，但這時需要手動輸入元素名稱。（��！注意：如果是vasp5版本的CONTCAR，請先刪除其中的元素符號那一行。） 例如： exe_Ctocif exe_Ctocif POSCAR Si O 這兩種方式都可以，前一種自動從贗勢里找元素名稱，后一種手動輸入元素名。 http://ishare.iask.sina.com.cn/f/9316742.html [ Last edited by nkleof on 2010-8-16 at 11:10 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 材料計算模擬實用技巧 | 材料計算模擬軟件 | 理論計算及軟件的使用 | Materials-studio |
» 猜你喜歡
295材料工程專碩求調(diào)劑
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有8人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
279求調(diào)劑
已經(jīng)有8人回復(fù)
-
328求調(diào)劑
已經(jīng)有8人回復(fù)
-
322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級。
已經(jīng)有9人回復(fù)
-
307求調(diào)劑
已經(jīng)有17人回復(fù)
-
286求調(diào)劑
已經(jīng)有5人回復(fù)
-
085600 295分求調(diào)劑
已經(jīng)有6人回復(fù)
-
287求調(diào)劑
已經(jīng)有13人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- linux shell循環(huán)的問題
已經(jīng)有3人回復(fù)
- VASP計算中，分波態(tài)密度問題
已經(jīng)有6人回復(fù)
- vasp計算含有磷原子結(jié)構(gòu)的態(tài)密度時磷原子為什么還有d態(tài)
已經(jīng)有14人回復(fù)
- VASP 計算DOS，態(tài)密度中費米能的作用
已經(jīng)有3人回復(fù)
- 用什么軟件把三張態(tài)密度圖畫在一個圖中？
已經(jīng)有11人回復(fù)
- 總的態(tài)密度為什么和分波態(tài)密度不相等？
已經(jīng)有19人回復(fù)
- 求處理vasp能帶的小程序
已經(jīng)有14人回復(fù)
- ENCUT 的腳本程序怎么用
已經(jīng)有15人回復(fù)
- vasp計算結(jié)果處理態(tài)密度要用到的程序
已經(jīng)有18人回復(fù)
- VASP計算態(tài)密度出現(xiàn)問題
已經(jīng)有12人回復(fù)
- 【求助】總態(tài)密度與分波態(tài)密度
已經(jīng)有9人回復(fù)
- 【求助】CASTEP中分析態(tài)密度
已經(jīng)有7人回復(fù)
- 【求助】求VASP計算ZnO能帶和態(tài)密度的詳細(xì)輸入設(shè)置
已經(jīng)有15人回復(fù)
- 【求助】態(tài)密度中6d貢獻(xiàn)部分丟失
已經(jīng)有3人回復(fù)
- 【求助】VASP計算電子態(tài)密度
已經(jīng)有34人回復(fù)
- 【求助】關(guān)于鋰原子vasp計算后的1s2s2p分波態(tài)密度處理
已經(jīng)有12人回復(fù)
3樓2010-08-16 11:09:21
2樓2010-08-16 10:57:56
5樓2010-08-16 12:02:14
cenwanglai
榮譽版主 (知名作家)
- 1ST強帖: 5
- 應(yīng)助: 46 (小學(xué)生)
- 貴賓: 8.842
- 金幣: 7440.4
- 帖子: 5306
- 在線: 1961.4小時
- 蟲號: 537452
★★★★★ 五星級,優(yōu)秀推薦
|
支持分享和原創(chuàng)。 split_dos源代碼中的邏輯關(guān)系理清楚沒？form=0/1， 及 # Total is first DOS if [ $form -eq 0 ]; then i=0 else i=0 fi 這個沒理解為什么要這樣寫。 還有后面的 while [ $i -le 0 ]; do #echo $i $start $end if [ $i -gt 0 ]; then sed -n ''$i' p' tmp.dat2 >| DOS$i fi if [ $i -gt 0 ]; then 這個循環(huán)從i<=0開始，但是上面無論怎樣，看起來i的初值是0呢。為什么這樣做循環(huán)呢？奇怪，好像可以大量簡化，即刪掉while [ $i -le 0 ]這段循環(huán)。 你說呢？ |
6樓2010-08-16 14:38:36
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +8 | HanBeiNingZC 2026-03-24 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +5 | 小小機(jī)靈蟲 2026-03-29 | 11/550 |
|
|
|
[考研] 297求調(diào)劑 +13 | 田洪有 2026-03-26 | 13/650 |
|
|
|
[考研] 317求調(diào)劑 +10 | 蛋黃咸肉粽 2026-03-26 | 10/500 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +10 | 麻辣魷魚 2026-03-27 | 11/550 |
|
|
|
[考研] 289求調(diào)劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運的三木 2026-03-29 | 5/250 |
|
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研]
 085404電子信息284分求調(diào)劑
+4
|
13659058978 2026-03-24 | 4/200 |
|