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文件

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

 

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM