| 5 | 1/1 | 返回列表 |
| 查看: 4154 | 回復(fù): 17 | ||||||
| 【獎(jiǎng)勵(lì)】 本帖被評價(jià)12次,作者nkleof增加金幣 10.8 個(gè) | ||||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||||
[資源]
【分享】態(tài)密度處理程序
|
||||||
|
之前發(fā)過一個(gè)簡單版本的態(tài)密度處理程序,這次更新一下。 split_dos.ksh是根據(jù)JMS NRL的腳本改的,其它程序是我自己寫的。 用法:只能在linux系統(tǒng)下,鍵入各個(gè)程序名稱運(yùn)行就可以,如果遇到平臺問題,sumdos程序無法執(zhí)行,可以用 gcc -lm sum.c 編譯源文件產(chǎn)生可執(zhí)行程序。 運(yùn)行split_dos.ksh會讀取DOSCAR,并生成總的態(tài)密度DOS0和各個(gè)原子的各軌道態(tài)密度DOS1~DOSn。(這里注意在計(jì)算DOS時(shí)INCAR文件里的NEDOS不要設(shè)置為>9999,五位數(shù)的NEDOS值會造成錯(cuò)誤) 然后運(yùn)行sumdos可以產(chǎn)生按照特定原子特定軌道模式生成的態(tài)密度和(注意:適用于ISPIN=2的計(jì)算結(jié)果,即DOSCAR里面共有有20列數(shù)據(jù))。 例: sumdos name -t 1 2 3 sumdos是程序名稱,name是一個(gè)自定義字符串,會加在生成的文件名里,可以作為說明(但是長度只有8個(gè)字符),-t表示生成各個(gè)軌道的態(tài)密度和,(還有兩種模式: -n 表示不加和各軌道態(tài)密度,仍舊生成從s到d各個(gè)軌道的態(tài)密度; -d 表示把p軌道,eg軌道,t2g軌道各自加和。!!仍然注意,和上次一樣,這里的eg指的是第14、15和18、19列數(shù)據(jù),t2g是剩下的3列數(shù)據(jù),如果有各個(gè)d軌道不是按照這樣排列的情形則生成的就不是eg和t2g),后面的數(shù)字表示要加和的原子標(biāo)號。上面的例子就是把第1、2、3號原子的所有軌道的態(tài)密度全部加起來。程序支持采用 n~m 形式表示要處理的第n到m個(gè)標(biāo)號連續(xù)的原子,例如要處理13到26號原子時(shí),不需要一個(gè)一個(gè)寫出原子標(biāo)號,只需要寫 13~26 就可以。 附送一個(gè)將CONTCAR文件轉(zhuǎn)換為可導(dǎo)入Materials Studio的cif文件的程序,使用很簡單,只需要在linux下,vasp任務(wù)的計(jì)算文件夾下,運(yùn)行exe_Ctocif,程序會自動尋找CONTCAR和POTCAR,并生成cif文件。另外也支持在只有POSCAR文件時(shí)轉(zhuǎn)換為cif,但這時(shí)需要手動輸入元素名稱。(!注意:如果是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 ] |
材料計(jì)算模擬實(shí)用技巧 | 材料計(jì)算模擬軟件 | 理論計(jì)算及軟件的使用 | Materials-studio |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +6 | 彼岸t 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +4 | 小小機(jī)靈蟲 2026-03-29 | 10/500 |
|
|
[考研] 289求調(diào)劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +6 | 101次希望 2026-03-28 | 6/300 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +9 | 辰沐5211314 2026-03-26 | 9/450 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
[考研] 壓國家一區(qū)線,求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研]
材料學(xué)碩,求調(diào)劑
6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|