应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » MS » MS Discover 动力学求助
111/2返回列表12下一页
查看: 3612  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yida

木虫 (小有名气)

[求助] MS Discover 动力学求助 已有1人参与

我做小分子在水体系中对金属表面的吸附模型，如下图，建好模型结构优化后进行了2ns 的NVT后进行了5ns的NVE，但是NVE任务失败，能量波动阈值采取默认5000，提示能量波动过大。
有以下问题请高手解惑：
体系建的正确的话是不是NVT NVE NPT系综都能运行？也就是说采用以上三种模型优化失败的话就是模型有问题？
为什么结构优化后进行NPT，立马就出错误信息，说能量波动值太大，这是什么原因？

模型.jpg



2ns NVT E-t.jpg



2ns NVT T-t.jpg



5ns NVE E-t.jpg



5ns NVE T-t.jpg
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-04-19 09:49:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

曾一晚

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
yida: 金币+5, ★★★很有帮助 2013-04-19 10:24:21
不一定是模型问题，也有可能是温控、压控算法问题，可以换个算法试试。
赞一下 回复此楼
2楼2013-04-19 10:19:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zmy0714

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
yida: 金币+5, ★★★很有帮助 2013-04-19 10:52:01
先跑NVT,在跑NPT 这样对系统的稳定性要好一些，也不容易报错

能量波动阈值采取默认5000，这个要改一下， 加大十来倍（但估计这个改动对你的模型意义不大）你的温度一下子直接上道了2500K ,这里可以尝试改变温控方法，通过控制温度变化范围，使温度收敛。具体是Velocity scale method、还是Nosé and Andersen methods、Berendsen method，自己琢磨下help文件。希望这些能帮到你。
赞一下 回复此楼
人生第一个五年规划，加油咯
3楼2013-04-19 10:24:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sikasun

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
yida: 金币+2, ★有帮助 2013-04-19 10:52:08
MS的discover模块的NPT就是不好跑的，一般体系要先稳定一些才能跑NPT，不然就出错，你可以先跑NVT让体系充分弛豫以后再NPT
赞一下 回复此楼
4楼2013-04-19 10:25:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yida

木虫 (小有名气)
引用回帖:
3楼: Originally posted by zmy0714 at 2013-04-19 10:24:00
先跑NVT,在跑NPT 这样对系统的稳定性要好一些，也不容易报错

能量波动阈值采取默认5000，这个要改一下， 加大十来倍（但估计这个改动对你的模型意义不大）你的温度一下子直接上道了2500K ,这里可以尝试改变温控方 ...

谢谢，我之前试过先跑NVT，把轨迹文件最后一帧提取出来后再跑NPT，但是收效甚微，你的意思是尝试改变下控温和控压方式是吧，我直接参考文献上的设置，选用的nose，还有个问题想请教下，最终体系平衡的标志是曲线平了，但是是不是温度到达到初始设置温度才能算体系达到平衡啊？
赞一下 回复此楼
5楼2013-04-19 10:51:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yida

木虫 (小有名气)
引用回帖:
4楼: Originally posted by sikasun at 2013-04-19 10:25:36
MS的discover模块的NPT就是不好跑的，一般体系要先稳定一些才能跑NPT，不然就出错，你可以先跑NVT让体系充分弛豫以后再NPT

谢谢
回复此楼
6楼2013-04-19 10:51:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

曾一晚

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
yida: 金币+5, ★★★很有帮助 2013-04-19 18:46:36
引用回帖:
5楼: Originally posted by yida at 2013-04-19 10:51:17
谢谢，我之前试过先跑NVT，把轨迹文件最后一帧提取出来后再跑NPT，但是收效甚微，你的意思是尝试改变下控温和控压方式是吧，我直接参考文献上的设置，选用的nose，还有个问题想请教下，最终体系平衡的标志是曲线平 ...

温度肯定是会一直波动的，只是波动范围大小不同而已，所以说楼主说的那个标准是不一定准确的。下面是别人跟我说的：分子动力学模拟是先确定一个温度，然后根据动能的关系来计算分子的速度，最后潜在的是速度平衡；跑动力学的时候，在很短时间内，温度就平衡了，但是后续还要做很长时间模拟，具体到底平衡没有，还得自己根据结构自己判断。
赞一下 回复此楼
7楼2013-04-19 11:19:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yida

木虫 (小有名气)
引用回帖:
7楼: Originally posted by 曾一晚 at 2013-04-19 11:19:13
温度肯定是会一直波动的，只是波动范围大小不同而已，所以说楼主说的那个标准是不一定准确的。下面是别人跟我说的：分子动力学模拟是先确定一个温度，然后根据动能的关系来计算分子的速度，最后潜在的是速度平衡； ...

非常感谢，具体如何判断能否举个例子？还有就是能不能讲点有关npt的知识？我将感激不尽
赞一下 回复此楼
8楼2013-04-19 18:45:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangjing1083

新虫 (小有名气)
请问楼主有没有画氢键的步骤呢？我是个外行，只是想利用MS模拟一下氢键的强度，麻烦你啦！
赞一下 回复此楼
9楼2015-01-13 18:33:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

njuedu007

禁虫 (著名写手)
★ ★ ★
yida(月只蓝代发): 金币+3, 感谢应助 2015-01-15 20:29:55
本帖内容被屏蔽
10楼2015-01-13 21:15:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yida 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 6/300 2026-11-05 14:06 by My20PsvIExYO
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-05 13:44 by My20PsvIExYO
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-05 13:33 by My20PsvIExYO
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 9/450 2026-11-05 13:22 by My20PsvIExYO
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 13:11 by My20PsvIExYO
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 13:11 by My20PsvIExYO
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 13:11 by My20PsvIExYO
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 ncTVBM6AOJHF 2026-11-04 7/350 2026-11-05 13:00 by My20PsvIExYO
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ncTVBM6AOJHF 2026-11-04 8/400 2026-11-05 13:00 by My20PsvIExYO
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-05 12:51 by My20PsvIExYO
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 6/300 2026-11-05 12:16 by My20PsvIExYO
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 ncTVBM6AOJHF 2026-11-04 7/350 2026-11-05 11:12 by PYZmcHr5kSaq
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 ncTVBM6AOJHF 2026-11-04 8/400 2026-11-05 11:01 by PYZmcHr5kSaq
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ncTVBM6AOJHF 2026-11-04 8/400 2026-11-05 10:59 by PYZmcHr5kSaq
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 10:48 by PYZmcHr5kSaq
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-05 10:48 by PYZmcHr5kSaq
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-05 10:39 by PYZmcHr5kSaq
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 hNpp4Y2CGvSC 2026-11-04 9/450 2026-11-05 09:20 by PYZmcHr5kSaq
[基金申请] 各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
关闭
请填处理意见
关闭