VASP计算小体系NVT分子动力学时温度波动大的问题如何解决

zhk9393

如题所示，我在用VASP计算大约20个原子的分子动力学NVT模拟时，遇到了最终温度波动大约在20-30K的问题（最终温度为300K）。想问大佬们可以通过哪些方法来减少最终温度的波动。
以下为我进行NVT模拟的参数：
PREC = Med
ISPIN = 1
ISTART = 0
ICHARG = 1
GGA    = PE
ENCUT  = 500
LREAL  = Auto
NCORE  = 12
ISIF = 2
EDIFFG = -0.02
#NPAR = 4
ISMEAR = 0
SIGMA = 0.05
ALGO = Fast
NELMIN = 5
NELM = 200
EDIFF = 1E-5
ISYM = 0
IBRION = 0
NSW = 2000
POTIM = 1
SMASS = -1
NBLOCK = 4
TEBEG = 300
TEEND = 300
NWRITE = 1
LWAVE = .F.
LCHARG = .F.
IVDW=11

MIOSAMA

感觉20个原子体系跑MD太小了，可以扩胞算一下

ZHJiao

关注一下 ，遇到了同样的问题