| 5 | 1/1 | 返回列表 |
| 查看: 4181 | 回復(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 ] |
材料計算模擬實用技巧 | 材料計算模擬軟件 | 理論計算及軟件的使用 | Materials-studio |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑!生物與醫(yī)藥專碩 +3 | 逆轉(zhuǎn)陸先生 2026-04-01 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +14 | 睿08 2026-03-27 | 15/750 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 5/250 |
|
|
[考研] 348求調(diào)劑 +9 | zzzzyk123 2026-04-01 | 9/450 |
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +14 | 秋二十二 2026-03-30 | 14/700 |
|
|
[考研] 材料調(diào)劑 +10 | Eujd1 2026-03-31 | 11/550 |
|
|
[考研] 086000生物與醫(yī)藥298調(diào)劑求助 +4 | 元元青青 2026-03-31 | 6/300 |
|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +10 | @taotao 2026-03-31 | 11/550 |
|
|
[考研] 初試301,代碼085701環(huán)境工程,本碩一致,四六級已過,有二區(qū)一作,共發(fā)表5篇論文 +3 | axibli 2026-04-01 | 3/150 |
|
|
[考研] 0856求調(diào)劑 +9 | 楒桉 2026-03-28 | 9/450 |
|
|
[考研] 080500-315分復(fù)試調(diào)劑 +9 | 上岸3821 2026-03-31 | 9/450 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +12 | Li李魚 2026-03-26 | 12/600 |
|
|
[考研] 293分求調(diào)劑,外語為俄語 +5 | 加一一九 2026-03-31 | 5/250 |
|
|
[考研] 福建理工大學(xué)材料學(xué)院先進合金團隊招收考研調(diào)劑學(xué)生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|