| 查看: 950 | 回复: 3 | |||||
| 当前主题已经存档。 | |||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 gamma337 的 4 个金币 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[交流]
【求助】分子动力学时间步长
|
|||||
|
下图是我在模拟前对势函数的时间步估计的一个能量输出图，系综为NVE，所示能量为总能（包括时能和动能），开始设定的时间步长为2fs，以后每隔30步增加1fs，共进行了300步，最后是11fs。 在最开始的时候，我们发现体系的能量是有个震荡的，然后才是平稳的。 请问： 1. 最初的微小的震荡是不是正常的？ 2. 也就是在第150步，也就是7fs的时候开始能量缓慢的上升，如果在模拟中，我2-7的哪个时间步长比较合适？ 以上两个问题，请大虾们给点建议，谢谢！！ |
回复此楼» 收录本帖的淘帖专辑推荐
 MD分子动力学 | 感悟人生 |
» 猜你喜欢
福州大学梁宇航副教授招收2026年申核制博士研究生/硕士研究生（理论计算方向）
已经有21人回复
-
福州大学梁宇航副教授招收2026年申核制博士研究生/硕士研究生（理论计算方向）
已经有20人回复
-
物理化学论文润色/翻译怎么收费?
已经有135人回复
-
深圳理工大学招聘AI4Science与分子模拟方向科研人员
已经有0人回复
-
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
-
全奖博士 英国利物浦大学 × 台湾清华大学 联合培养
已经有0人回复
-
寻博导
已经有10人回复
-
26申博求博导-第一性原理/计算化学-可提前进组
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
3楼2008-05-29 09:13:51
leongxf
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 464.5
- 散金: 1
- 帖子: 30
- 在线: 6.2小时
- 虫号: 490352
- 注册: 2008-01-03
- 专业: 金属材料的合金相、相变及
2楼2008-05-28 20:34:17
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