求助VASP热力学积分计算自由能

phyang

根据VASP-wiki以及PRL 2018 121, 195701 的介绍，VASP可以直接设置SCALEE来表示热力学积分里面的lambda进行热力学积分计算自由能。我最近想算液态铁的自由能，对其进行了尝试，发现当这个lambda设置较小的时候（接近理想气体），MD跑了一两百步直接崩，因为原子靠太近了（如图所示）。根据文献的建议，当lambda较小的时候，建议用LANGEVIN热浴或者ANDERSEN热浴，我均尝试了，都失败了。有没有大佬用过这个功能啊，提供宝贵的建议，不胜感激！
这是我的 INCAR，kpoints用的一个gamma点

ISTART =  0
LREAL  = A
ENCUT  = 600
PREC   =  N
LWAVE  = F
LCHARG = .FTRUE.        (Write CHGCAR or not)
NPAR = 4

Electronic Relaxation
ISMEAR =  -1
SIGMA  =  0.4308
EDIFF  =  1E-05

LMAXMIX  = 4
#ALGO = Fast

Molecular Dynamics
IBRION =  0            (Activate MD)
NSW    =  50000
POTIM  =  0.5            (Timestep in fs)
#SMASS  =  3
TEBEG  =  5000     (Start temperature K)
TEEND  =  5000
MDALGO =  1
ANDERSEN_PROB = 0.02
ISYM   =  0          (Switch symmetry off)
ISIF = 2

LSCALAPACK = F

SCALEE = 0.000354605899

NBANDS = 500