各位老师，本人最近在使用vasp计算petn晶体的弹性常数。步骤如下，首先优化晶体，优化的incar和kpoints如下：
incar
system=petn
ismear=0
sigma=0.05
encut=1000
ediff=0.0001
ediffg=0.001
prec=accurate
potim=0.0100
lcharg=false
lwave=false
ibrion=2
isif=3
nfree=4
nsw=100
lreal=auto

kpoints
k
0
gamma
   3   3   4
0.0  0.0  0.0

优化结束后（56步优化结束），将优化得到的contcar文件重新命名为poscar文件（计算弹性常数用）。使用下列incar作为计算，kpoints保持不变：
system=petn
ismear=0
sigma=0.05
encut=1000
ediff=0.0001
ediffg=0.001
prec=accurate
potim=0.01
lcharg=false
lwave=false
ibrion=6
isif=4
nfree=4
nsw=1
lreal=auto
为了测得更多结果，更改encut的值为：400,500,600,700,800,1000，其他的所有都保持不变，为什么计算的弹性常数差距如此大：


encut=400.png



encut=500.png



encut=600.png



encut=700.png



encut=800.png



encut=1000.png 返回小木虫查看更多