查看: 1390  |  回复: 6
当前主题已经存档。

笑问天

木虫 (著名写手)

[交流] 【求助】分子动力学模拟温度的问题

我用MS的分子动力学模块模拟水分子层在晶面上的吸附,为了消除建模时水分子的初始构象的影响,我在500K下先进行了5ps的平衡。而审稿人却认为在这个温度下水分子都蒸发光了,因此不同意我的做法。但是我也是参考别人的文献才这么做的啊,而且在模拟过程中没发现水分子都跑光光了呀。不过如果采用anneal模块的话,那么水分子确实会分散到整个盒子中去,但是用dynamic的话,就不会跑掉。请问遇到这个问题,要怎么跟审稿人解释呢。

[ Last edited by zdhlover on 2009-12-10 at 15:05 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyx19840628

木虫 (著名写手)

鹰的眼睛,狼的耳朵,豹的速 ...

别人文献上说了吗?
人生的最大遗憾莫过于错误地坚持了不该坚持的,轻易地放弃了不该放弃的……
2楼2009-01-06 20:58:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笑问天

木虫 (著名写手)

文献上没有谈到温度多少度合适,只是简单的说明了下在这个温度点下计算,我也没多想,就照猫画虎的借用了……
3楼2009-01-07 00:12:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuhanhgf2002

金虫 (正式写手)

★ ★ ★ ★ ★
yyx19840628(金币+5,VIP+0):谢谢
因为动力学中采用周期性边界条件来控制连续溶液模型,所以在任何温度下水都不会跑散。这篇文献做的应该是高温动力学模拟,在500k的温度下模拟的目的是想通过提高体系的温度来缩短模拟达到平衡的时间,但是温度不易过高,一般600k以下。你应该把这种高温动力学的思想介绍给审稿人,审稿人不一定是搞计算的,所以回产生误解。还有就是不光是高温动力学模拟,还有人做模拟退火分子动力学,温度还高一些,要到1000K以上,如果照审稿人的思路,材料可能都变性了,模拟毕竟是模拟,跟现实还是有区别的。

[ Last edited by wuhanhgf2002 on 2009-1-9 at 09:02 ]
4楼2009-01-07 08:21:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢分析 受益
你在500K下运行,只是为了消除初始构像的影响,构造一个比较随机和混乱的初始构型而已,这没有任何问题。
为什么500K下没有蒸发呢?因为你并不是在模拟蒸发过程,所谓的500K只是为了给分子一个较大的动能,加快弛豫罢了,你这一步应该是在等容条件下弛豫的,而在等容和体系封闭的条件下,水不可能汽化。你这样跟审稿人讲就是了。
5楼2009-01-07 12:16:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zz422

铜虫 (小有名气)

支持上面的分析
6楼2009-01-07 17:13:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

笑问天

木虫 (著名写手)

谢谢楼上几位高手的分析,平时也有不少人跟我扯这个温度的问题,说实在不明白温度怎么能这样设。我虽然大概知道缘由,可是又说不清道不明。现在要赶紧把高手的分析摘抄下来,铭记在心,随时拿出来显摆,哈哈。
7楼2009-01-07 21:21:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 笑问天 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[硕博家园] 极限拉扯,如何抉择 +5 Isak0320 2026-10-12 7/350 2026-10-12 18:03 by 有人要害朕
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见