计算化学公社

标题: 求助, VASP进行NVT模拟，温度一直再下降，SMASS这个值如何调整？ [打印本页]

作者
Author: 薛定谔猫之谜 时间: 2019-12-17 20:23
标题: 求助, VASP进行NVT模拟，温度一直再下降，SMASS这个值如何调整？
本帖最后由薛定谔猫之谜于 2019-12-17 20:30 编辑

原计划使用VASP做NVT MD模拟，设定了温度400K, 但是跑完20ps后，发现温度越来越低

(右图)显然已经不是NVT了。VASP教程上说是热浴的频率和系统的声子频率不一致

是SMASS设置不合理。请问这里系统的声子频率怎么得到呢？具体怎么得到的这个频率？

(Fig. 1是温度随着时间步数的变化；Fig. 2 是结构图，其中绿, 碳;白，氢；棕，钯)

附

INCAR：

# Startparameter forthis run:

PREC = N

ISPIN = 1

ENCUT = 400

# Electronic Relaxation

EDIFF = 1E-05

ALGO = Fast

ISTART = 1

ICHARG = 2

NELMDL = -3

NELMIN = 5

NELM = 200

LREAL = Auto

LCHARG = .F.

LWAVE = .F.

#Exchange correlation treatment

GGA = PE

#Molecular Dynamic

IBRION = 0

NSW = 10000

EDIFFG = -1E-02

POTIM = 2

SMASS = 0

TEBEG = 400

TEEND = 400

MDALGO = 1

ISYM = 0

POMASS = 10 12 106 14 # H C Pd N

作者
Author: 薛定谔猫之谜 时间: 2019-12-20 11:08
1. 切换成Nose-hoover恒温器，即MDALGO = 2。SMASS才有意义。
2. 使用Anderson恒温器，没有设置ANDERSEN_PROB = 0≤[real]≤1，相当于跑的NVE系综。

作者
Author: Weldingspock 时间: 2023-6-15 09:31
楼主解决了吗，怎么计算声子谱呢

作者
Author: wwjie 时间: 2024-3-10 15:56
遇到同样的问题

作者
Author: wwjie 时间: 2024-3-10 15:57
你的MDALGO设置的不对吧

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