计算化学公社

 找回密码 Forget password
 注册 Register
Views: 17335|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 求助, VASP进行NVT模拟,温度一直再下降,SMASS这个值如何调整?

[复制链接 Copy URL]

6

帖子

0

威望

161

eV
积分
167

Level 3 能力者

本帖最后由 薛定谔猫之谜 于 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

fig2.gif (12.24 KB, 下载次数 Times of downloads: 71)

Fig. 2 20ps左右的动态结构变化

Fig. 2 20ps左右的动态结构变化

图片1.jpg (36.02 KB, 下载次数 Times of downloads: 53)

Fig. 1 温度随时间步数的变化

Fig. 1 温度随时间步数的变化

6

帖子

0

威望

161

eV
积分
167

Level 3 能力者

2#
 楼主 Author| 发表于 Post on 2019-12-20 11:08:27 | 只看该作者 Only view this author
1. 切换成Nose-hoover恒温器,即MDALGO =  2。SMASS才有意义。
2. 使用Anderson恒温器,没有设置ANDERSEN_PROB = 0≤[real]≤1,相当于跑的NVE系综。

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
weilu + 4

查看全部评分 View all ratings

207

帖子

1

威望

1433

eV
积分
1660

Level 5 (御坂)

3#
发表于 Post on 2023-6-15 09:31:49 | 只看该作者 Only view this author
楼主解决了吗,怎么计算声子谱呢

21

帖子

0

威望

450

eV
积分
471

Level 3 能力者

4#
发表于 Post on 2024-3-10 15:56:57 | 只看该作者 Only view this author
遇到同样的问题

21

帖子

0

威望

450

eV
积分
471

Level 3 能力者

5#
发表于 Post on 2024-3-10 15:57:48 | 只看该作者 Only view this author
你的MDALGO设置的不对吧

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-18 02:35 , Processed in 1.243567 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list