24小时热门版块排行榜    

查看: 3004  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
信息提示
请填处理意见