24小时热门版块排行榜    

查看: 1057  |  回复: 3
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 gamma337 的 4 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gamma337

银虫 (小有名气)

[交流] 【求助】分子动力学时间步长

下图是我在模拟前对势函数的时间步估计的一个能量输出图,系综为NVE,所示能量为总能(包括时能和动能),开始设定的时间步长为2fs,以后每隔30步增加1fs,共进行了300步,最后是11fs。
在最开始的时候,我们发现体系的能量是有个震荡的,然后才是平稳的。
请问:
1. 最初的微小的震荡是不是正常的?
2. 也就是在第150步,也就是7fs的时候开始能量缓慢的上升,如果在模拟中,我2-7的哪个时间步长比较合适?

以上两个问题,请大虾们给点建议,谢谢!!
回复此楼
美元一起赚: http://usearn.5d6d.com/index.php
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gamma337

银虫 (小有名气)

我们的模拟尽量要做到在有限的时间步内,模拟更长的时间,因此,我们需要一个能确定体系稳定且较大的时间步长。

上个图我们可以看出在5fs的时候体系的能量开始上升,我在模拟的时候实际取的3fs,这个度我很难说的准确,也许取的更小些体系也比较稳定。
美元一起赚: http://usearn.5d6d.com/index.php
3楼2008-05-29 09:13:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

leongxf

银虫 (初入文坛)


lei0736(金币+1,VIP+0):谢谢
不了解你为什么要改变步长,一般分子动力学的步长选1fs,刚开始能量震荡是没错的,搞不懂你的能量为什么不守恒。
2楼2008-05-28 20:34:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leongxf

银虫 (初入文坛)

★ ★ ★
lei0736(金币+2,VIP+0):谢谢一直关注 鼓励新人
gamma337(金币+1,VIP+0):这个问题是我在找出来时间步取多大值的时候系统开始不稳定。谢谢你
昨晚我看错了,如果总能量的图,前面就不应该有震荡(如果你不人为干扰系统),你应该把势能和动能都画出来,看一下,总能量是哪一部分引起的,是因为动能,还是势能,如果是动能,那就是速度越来越大了,体系温度不断升高。如果是势能引起的,就有可能是势函数的问题。
你想模拟尽可能久的时间,是否可以分别作几个体系,分别选取不同的步长,然后比较各体系的情况。
你这样一个体系,变步长,能量的增加都不知什么因素引起的。还有就是选取步长为1fs,由于势函数的错误,也可能导致10000fs后体系量能剧增,更别说几百分秒了。
4楼2008-05-29 09:31:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +57 AI8RaGaPaSCR 2026-08-07 60/3000 2026-08-14 04:45 by bxqFb8Dih2p1
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +11 医学老男孩 2026-08-13 29/1450 2026-08-14 00:04 by 大快人心啦
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
请填处理意见