| 5 | 1/1 | 返回列表 |
| 查看: 1774 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
那年的小艾金蟲 (初入文壇)
|
[求助]
關(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中,然后運行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……。另外在運行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類似。 |


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26年申博 +3 | staryer 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 安全工程 285 求調(diào)劑 +3 | Xinyu56 2026-04-01 | 4/200 |
|
|
[考研] 286求調(diào)劑 +16 | PolarBear11 2026-03-26 | 16/800 |
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +22 | 材料學(xué)257求調(diào)劑 2026-03-28 | 23/1150 |
|
|
[考研] 353求調(diào)劑 +4 | 拉鉤不許變 2026-04-01 | 4/200 |
|
|
[考研] 英一數(shù)二生物信息學(xué)287分,本科生物科學(xué),求調(diào)劑 +3 | 碧水xyz 2026-03-29 | 4/200 |
|
|
[考研] 324求調(diào)劑 +4 | 想上學(xué)求調(diào) 2026-04-01 | 5/250 |
|
|
[考研] 085600,321分求調(diào)劑 +13 | 大饞小子 2026-03-31 | 13/650 |
|
|
[考研] 一志愿中農(nóng)0710生物學(xué),微生物方向總分338求調(diào)劑 +3 | 柒xxxx. 2026-03-26 | 3/150 |
|
|
[考研] 262求調(diào)劑 +9 | 勵志一定發(fā)文章 2026-03-31 | 10/500 |
|
|
[考研] 318求調(diào)劑 +8 | 七憶77 2026-04-01 | 8/400 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +6 | 篤行致遠(yuǎn). 2026-03-28 | 6/300 |
|
|
[考研] 【調(diào)劑】一志愿廈大生物與醫(yī)藥調(diào)劑 +3 | Echo蝦米 2026-03-31 | 3/150 |
|
|
[考研] 085600,專業(yè)課化工原理,320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|