1. 程式人生 > 實用技巧 >VASP靜態自洽計算的輸入檔案,執行和輸出檔案

VASP靜態自洽計算的輸入檔案,執行和輸出檔案

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檔案

#電子在導空間上的分佈狀態