| 5 | 1/1 | 返回列表 |
| 查看: 1753 | 回復(fù): 4 | ||
那年的小艾金蟲 (初入文壇)
|
[求助]
關(guān)于態(tài)密度計算的一些問題
|
|
本人先做了晶體的結(jié)構(gòu)優(yōu)化后,將INCAR CONTCAR POTCAR KPOINTS考入新建的DOS文件夾中并將CONTCAR改成了POSCAR,準(zhǔn)備進(jìn)行態(tài)密度的計算。請問INCAR中的設(shè)置哪些需要修改的(體系中有Eu,要考慮自旋) 新手一枚,望各位多加指點。 |
|
VASP中電子態(tài)密度計算的流程 主要分成三步:一、結(jié)構(gòu)優(yōu)化;二、靜態(tài)自洽計算;三、非自洽計算 以Al-FCC為例子 第一步 結(jié)構(gòu)優(yōu)化 輸入文件(INCAR, POTCAR, POSCAR, KPOINT) INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=2 NSW=50 ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在勢庫中拷貝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Direct 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5 KPOINT 文件 Automatic generation 0 Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步 靜態(tài)自洽計算 INCAR: PREC = Medium,ISTART = 0,ICHARG = 2,ISMEAR = -5 輸入文件(INCAR, POTCAR, POSCAR, KPOINT) INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在勢庫中拷貝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 0.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation 0 Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步計算是在結(jié)構(gòu)優(yōu)化的結(jié)果上進(jìn)行的,所以開始第二步的時候,將第一步中的輸入文件INCAR, POTCAR, POSCAR, KPOINT 以及 C* 文件放入靜態(tài)自洽計算中去,并且將CONTCAR 拷貝到 POSCAR中,然后運(yùn)行VASP。計算結(jié)果中的Fermi能是準(zhǔn)確的,需要記下。 第三步 非自洽計算 INCAR: PREC = Medium,ICHARG = 11,ISMEAR = -5,LORBIT = 10 或者 11(這時可不設(shè)RWIGS),ISTART = 1 在第二步自洽計算的基礎(chǔ)上進(jìn)行,修改輸入文件INCAR, POTCAR, POSCAR, KPOINT。 INCAR文件 System=Al ISTART=1 ISMEAR=-5 SIGMA=0.2 ICHARG=11 RWIGS=1.402 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在勢庫中拷貝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 1.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation 0 Mohkorst Pack 21 21 21 0.0 0.0 0.0 VASP中文手冊:在進(jìn)行能帶和DOS計算時,ISMEAR 不能使用N階MP方法。因為MP方法在空軌道上有負(fù)的占據(jù),所以求得的能帶和DOS是不正確的。但是從其它地方看到---“提示:在計算能帶結(jié)構(gòu)時,采ISMEAR = 0或1對結(jié)果的影響非常小,可以認(rèn)為是一樣的。但是不能采用ISMEAR = -5 或-4!盜SMEAR到底多少? 計算能帶:ICHARG = 11 導(dǎo)體的話,用ISMEAR=1; 半導(dǎo)體或絕緣體,用ISMEAR=0 。 計算 DOS: ICHARG = 11 ISMEAR = -5 計算的時候,金屬可用0、1,非金屬不要大過0,體材料可用-4、-5(面的話就用-1、0吧?) 參考http://210.45.66.120/doc/vasp5.0/node224.html 設(shè)置完成后進(jìn)行計算,計算完后,得到包含了態(tài)密度值的DOSCAR文件,采用split_dos對態(tài)密度文件DOSCAR進(jìn)行分割,得到總態(tài)密度DOS0,各個原子的分波態(tài)密度DOS1,DOS2……。另外在運(yùn)行split_dos程序?qū)OSCAR文件分割時,要保證當(dāng)前目錄下有對應(yīng)的OUTCAR和POSCAR文件。 分割后的DOS0,DOS1…等文件的能量值是以費米能級作為能量參考零點。DOS0的第一列數(shù)據(jù)是能量值,單位為eV;第二列數(shù)據(jù)是總態(tài)密度的值,單位 State/eV.unit cell;第三列數(shù)據(jù)是總態(tài)密度的積分值,也就是電子數(shù),單位為electrons。DOS1是第一個原子的分波態(tài)密度值,其中的第一列數(shù)據(jù)是能量值,單位為eV;第二、三、四列數(shù)據(jù)分別對應(yīng)于s、p、d態(tài)的分波態(tài)密度值,單位為State/eV.atom。其他的DOS文件與DOS1類似。 |


| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調(diào)劑 +7 | 梨花珞晚風(fēng) 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 |  ̄^ ̄゜汗 2026-03-19 | 4/200 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|