24小时热门版块排行榜    

查看: 2967  |  回复: 15

newhoperoad

木虫 (正式写手)


[交流] 碳纳米管MD模拟中积分时间步长选取时存在的问题

在使用MS的discover模拟碳纳米管,随意建立一个碳纳米管模型,然后在setup中选择compass力场,其他选择默认;dynamics中选择NPT系综,步数设置为500000,积分时间步长设定为1fs 或者0.5fs,温度使用Anderrsen控制在室温298K,其他为默认值;然后开始计算,计算中出现的问题是:温度波动很大,无法保持温度恒定在298K,这是为什么呢?
对于强悍的MS,不应该在这么基本的运算上出问题啊?请指教,谢谢!
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )
★ ★ ★ ★ ★
newhoperoad: 金币+5 2014-01-03 22:32:31
如果一定要精确控温
用Velocity scale恒温器
12楼2014-01-03 17:06:48
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

newhoperoad

木虫 (正式写手)


我猜想,是不是因为碳纳米管是自身非常稳定的刚性结构,在自身结构平衡的基础上,任何一点的外界环境变化,都会引起分子结构整体连锁式的运动,从而使动能增加,温度升高。不知道是不是这样的,或者有别的什么原因?
2楼2014-01-03 13:58:11
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sciencejoy

新虫 (著名写手)


跑平衡了没有?
3楼2014-01-03 15:32:22
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hanhan1020

木虫 (正式写手)


★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
newhoperoad: 金币+5, 第一个回复,谢谢啦 2014-01-03 20:00:52
恒定在298K?你所谓的恒定是什么意思?只要在298+—10k内都算可以,如果是后期波动大 考虑噪音项影响
4楼2014-01-03 15:52:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

newhoperoad

木虫 (正式写手)


引用回帖:
4楼: Originally posted by hanhan1020 at 2014-01-03 15:52:14
恒定在298K?你所谓的恒定是什么意思?只要在298+—10k内都算可以,如果是后期波动大 考虑噪音项影响

波动超过了+-10K,请问您说的噪音项指的是精度吗?

[ Last edited by newhoperoad on 2014-1-3 at 16:15 ]
5楼2014-01-03 16:09:24
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

newhoperoad

木虫 (正式写手)


引用回帖:
3楼: Originally posted by sciencejoy at 2014-01-03 15:32:22
跑平衡了没有?

跑够5ns,还是波动比较大,应该是@hanhan1020 说的噪音问题
6楼2014-01-03 16:13:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hanhan1020

木虫 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
你改控温器  Nose 试试
7楼2014-01-03 16:17:11
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hanhan1020

木虫 (正式写手)


噪音项 可以通过控温器的选择减小的,你用nose试试
8楼2014-01-03 16:18:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

newhoperoad

木虫 (正式写手)


引用回帖:
8楼: Originally posted by hanhan1020 at 2014-01-03 16:18:42
噪音项 可以通过控温器的选择减小的,你用nose试试

我用nose试了一下,还是存在温度波动过大的问题,还可能是什么原因呢?
9楼2014-01-03 16:42:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hanhan1020

木虫 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
9楼: Originally posted by newhoperoad at 2014-01-03 16:42:31
我用nose试了一下,还是存在温度波动过大的问题,还可能是什么原因呢?...

温度波动有多大?
10楼2014-01-03 16:46:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
温度上下涨落是正常现象,平衡时围绕一个稳定值即可。
11楼2014-01-03 17:05:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

newhoperoad

木虫 (正式写手)


引用回帖:
10楼: Originally posted by hanhan1020 at 2014-01-03 16:46:56
温度波动有多大?...

温度波动在50~70k
13楼2014-01-03 19:22:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

newhoperoad

木虫 (正式写手)


引用回帖:
12楼: Originally posted by 月只蓝 at 2014-01-03 17:06:48
如果一定要精确控温
用Velocity scale恒温器

Velocity Scale 确实实现了控温,但是动画图片显示碳纳米管振动非常剧烈,有突然的较大的位移,结构整体有跳跃,这是怎么回事啊?用其他方法控温,振动没有那么剧烈,运动是连续的。
14楼2014-01-03 19:58:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hanhan1020

木虫 (正式写手)


多优化一下结构
15楼2014-01-03 20:01:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

newhoperoad

木虫 (正式写手)


引用回帖:
15楼: Originally posted by hanhan1020 at 2014-01-03 20:01:34
多优化一下结构

优化过了,碳纳米管本身就是优化好的结构,结果没什么变化,还是存在这个问题
16楼2014-01-03 22:29:29
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 newhoperoad 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:33 by ZPa0EcMwuECS
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:22 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:18 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 01:03 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见