記錄一下最近計算的不同壓力下Cr3C2的彈性常數

1. vasp直接計算

vasp5.1版本后支持設置IBRION=6來直接計算材料的彈性常數，最直接的做法就是在INCAR中設置以下三個參數

IBRION = 6 ISIF = 3 NFREE = 4

具體的INCAR設置參見

ISTART=0 ICHARG=2 PREC=Accurate GGA = PE ADDGRID =.TRUE.# Electronic relaxation ENCUT=550 EDIFF=1E-8ISMEAR=0 SIGMA=0.05 NELMIN=5# Ionic relaxation NSW=1 EDIFFG=-1e-3 ISIF=3IBRION=6 #計算彈性常數 POTIM=0.015 NFREE = 4# Other Tags PSTRESS=50 #計算在5GPa壓強下的彈性常數#Write flags LWAVE=.FALSE. LCHARG=.FALSE.KSPACING = 0.25 #設置K點密度

計算結束以后就會在OUTCAR中輸出彈性常數（注意單位是kBar，1GPa=10kBar）

Cr3C2屬于正交晶系，總共有9個獨立矩陣元(C11,C12, C13,C22，C23，C33，C44，C55，C66)

2. vaspkit計算

要計算不同壓強下的彈性常數需要現在不同壓強下進行結構優化，再分別計算彈性常數。具體計算細節后續參見博客。

計算結果如下

除了彈性常數外，還可以得到柔順系數，并且楊氏模量，體積模量等都不需要手動去算，都會直接輸出。

總結

以上是生活随笔為你收集整理的VASP计算弹性常数的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。