计算化学公社

标题: 求助VASP热力学积分计算自由能 [打印本页]

作者
Author: phyang 时间: 2024-11-17 09:56
标题: 求助VASP热力学积分计算自由能
根据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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)