| 查看: 922 | 回复: 3 | |||||
| 当前主题已经存档。 | |||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 gamma337 的 4 个金币 | |||||
[交流]
【求助】分子动力学时间步长
|
|||||
|
下图是我在模拟前对势函数的时间步估计的一个能量输出图，系综为NVE，所示能量为总能（包括时能和动能），开始设定的时间步长为2fs，以后每隔30步增加1fs，共进行了300步，最后是11fs。 在最开始的时候，我们发现体系的能量是有个震荡的，然后才是平稳的。 请问： 1. 最初的微小的震荡是不是正常的？ 2. 也就是在第150步，也就是7fs的时候开始能量缓慢的上升，如果在模拟中，我2-7的哪个时间步长比较合适？ 以上两个问题，请大虾们给点建议，谢谢！！ |
回复此楼» 收录本帖的淘帖专辑推荐
 MD分子动力学 | 感悟人生 |
» 猜你喜欢
期待科研合作，共同发表论文
已经有216人回复
-
如何掌握一门外语？
已经有25人回复
-
物理化学论文润色/翻译怎么收费?
已经有266人回复
-
国家级人才团体课题组招收2026届博士
已经有1人回复
-
气相色谱点不着火
已经有4人回复
-
化工安全有关的复习提纲
已经有0人回复
-
国家级人才课题组招收2026年入学博士
已经有2人回复
-
散金
已经有85人回复
-
浙江师范大学国家杰青杨启华教授团队招收2026年博士研究生
已经有33人回复
» 本主题相关商家推荐: (我也要在这里推广)
leongxf
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 464.5
- 散金: 1
- 帖子: 30
- 在线: 6.2小时
- 虫号: 490352
- 注册: 2008-01-03
- 专业: 金属材料的合金相、相变及
2楼2008-05-28 20:34:17
3楼2008-05-29 09:13:51
leongxf
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 464.5
- 散金: 1
- 帖子: 30
- 在线: 6.2小时
- 虫号: 490352
- 注册: 2008-01-03
- 专业: 金属材料的合金相、相变及
★ ★ ★
lei0736(金币+2,VIP+0):谢谢一直关注 鼓励新人
gamma337(金币+1,VIP+0):这个问题是我在找出来时间步取多大值的时候系统开始不稳定。谢谢你
lei0736(金币+2,VIP+0):谢谢一直关注 鼓励新人
gamma337(金币+1,VIP+0):这个问题是我在找出来时间步取多大值的时候系统开始不稳定。谢谢你
|
昨晚我看错了，如果总能量的图，前面就不应该有震荡（如果你不人为干扰系统），你应该把势能和动能都画出来，看一下，总能量是哪一部分引起的，是因为动能，还是势能，如果是动能，那就是速度越来越大了，体系温度不断升高。如果是势能引起的，就有可能是势函数的问题。 你想模拟尽可能久的时间，是否可以分别作几个体系，分别选取不同的步长，然后比较各体系的情况。 你这样一个体系，变步长，能量的增加都不知什么因素引起的。还有就是选取步长为1fs，由于势函数的错误，也可能导致10000fs后体系量能剧增，更别说几百分秒了。 |
4楼2008-05-29 09:31:25