计算化学公社

 找回密码 Forget password
 注册 Register

npt控压出现问题

查看数: 13499 | 评论数: 4 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
mol
发布时间: 2018-1-8 13:55

正文摘要:

各位前辈,小弟目前在模拟原油体系在二氧化硅管内的输运现象,体系如下图: 我想在z方向把油的密度压缩到平衡,使用mdp文件如下: cpp               &n ...

回复 Reply

虔诚的人 发表于 Post on 2018-4-27 19:53:36
sobereva 发表于 2018-4-26 19:39
首先我并不建议用all-bonds约束,引入的人为因素不可忽略
对键进行约束,主要目的是为了消除高频振动, ...

谢谢,sob老师,理解了
sobereva 发表于 Post on 2018-4-26 19:39:05
虔诚的人 发表于 2018-4-26 17:36
sob老师,请问时间不长0.1fs和0.2fs是怎么选择的?为什么time steps=0.2 ps后就没必要设置all-bond了。{: ...

首先我并不建议用all-bonds约束,引入的人为因素不可忽略
对键进行约束,主要目的是为了消除高频振动,从而可以用更大的步长。0.2fs已经算是非常小了,这种情况下没有必要对键做任何约束。一般情况,<=1fs的步长都没必要做任何约束。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
tjuptz + 2

查看全部评分 View all ratings

虔诚的人 发表于 Post on 2018-4-26 17:36:48
sobereva 发表于 2018-1-8 18:20
既然步长都0.2fs了,all-bonds约束也没必要设了。先去了试试
乍看没有明显问题,但一些细节不好说。控压的 ...

sob老师,请问时间不长0.1fs和0.2fs是怎么选择的?为什么time steps=0.2 ps后就没必要设置all-bond了。
sobereva 发表于 Post on 2018-1-8 18:20:38
既然步长都0.2fs了,all-bonds约束也没必要设了。先去了试试
乍看没有明显问题,但一些细节不好说。控压的设定反复再调节调节,再看看top文件有无毛病,反复试试

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

GMT+8, 2025-8-16 18:54 , Processed in 0.160517 second(s), 26 queries , Gzip On.

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