VASP輸入文件
1.INCAR文件
# 計算的體系的名稱
SYSTEM = Hybrid
# I/O設置(讀入、讀出)
ISTART = 0 # 0代表一個全新的計算
ICHARG = 2
LWAVE = .FALSE.
LCHARG = .TRUE. #輸出電子密度文件
LVOT = .TRUE. #loacl potential文件
LELF = .TRUE. #輸出電子局域化函數
LORBIT = 11 #輸出材料的分波態密度
#以上四個文件一般在需要計算准確電子結構(自洽運算)時才會出現
# Electronic Relaxation(電子步)
ENCUT = 600 #平面階段能,由贗勢決定(POTCAR中的ENMAX*1.25~1.50)(數值越大,精度越高)
NELM = 100 #最大電子步數量,默認60步(難收斂體系,如過渡金屬,可設置為200、300~)
ALGO = Fast #自洽循環算法(Fast,Normal,VaryFast)
PREC = Accurate #精度設置
ISMEAR = 0 # 金屬0+,非金屬0-(金屬不可以設置為負數)
SIGMA = 0.05 # 划分布里淵區,展寬
EDIFF = 1E-5 #電子步收斂精度
AMIX = 0.1 #
BMIX = 0.01 #加快收斂
NEDOT = 2000
EMIN = -10.0 #兩個
EMAX = 10.0 #細化能級
# Ionic Relaxation (離子步)
IBRION = 2 #離子弛豫的算法(2較穩定)
NSW = 300 #最大離子步的個數
EDIFFG = -1E-2 #離子步收斂標准(默認是EDIFF*10)
ISIF = 3 # 對晶胞的弛豫方法(3:全弛豫 2:固定體積的弛豫 4:固定體積但允許形狀改變)
# Polarization (含磁性體系)
ISPIN = 2 # 自旋極化
LOSRBIT = .False. #考慮非線性磁性(默認不考慮)
# Parallization(並行計算的相關設定)
LREAL = Auto #計算進行在實空間或是導空間
NPAR = 1 #取節點數量或核心數量的根號值
LPLANE = .FALSE. #針對平面波的設置
-
‘#’ 表示注釋符,后面的內容不會被執行
2.KPOINTS文件
# 划分布里淵區網格
Automatic mesh
0
Gamma #划分網格
8 8 8
0.0 0.0 0.0
# 另一種KPOINTS文件形式
Automatic mesh
0
Auto
25.0
3.POSCAR文件
# 通過VESTA建模得到
# 包括體系名稱,晶胞基矢信息,原子元素種類、數量和具體坐標
4.POTCAR文件
# 從贗勢庫中得到的贗勢信息,一般不需要進行操作
VASP程序的執行
mpirun -np x vasp_std < INCAR
# x是指CPU的線程數(核心數)
VASP輸出文件
1.OSZICAR
#第一行內容
F #體系的自由能
E0 #最后體系的總能量
dE #體系能量的改變量
# 使用linux命令得到文件中我們需要的數據
grep E0 OSZICAR #打印出OSZICAR文件中所有含有'E0'的行
2.OUTCAR文件
# 包含VASP全部輸入信息的總結(日志文件)
# 最后一行:計算成本,磁性信息,電荷信息
grep 'reached required accurary' OUTCAR
# 若OUTCAR文件中存在上述字符串,則表明計算成功收斂;若沒有,則未收斂
3.DOSCAR文件
# 材料的態密度信息(材料在不同能級上各自的占據情況/概率)(5列、3列)
# 若在電子步設置了LORBIT=11,則會輸出分波態密度(19列)更大
4.CHGCAR文件
#電子在空間位置上的分布,可直接通過VESTA文件畫圖
5.LOCALCAR文件
#電子在空間中的能級分布
6.ELFCAR文件
#定義的電子局域化函數
7.EIGENVAL文件
#電子在導空間上的分布狀態