| 查看: 1027 | 回复: 3 | |||||
| 当前主题已经存档。 | |||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 gamma337 的 4 个金币 | |||||
[交流]
【求助】分子动力学时间步长
|
|||||
|
下图是我在模拟前对势函数的时间步估计的一个能量输出图，系综为NVE，所示能量为总能（包括时能和动能），开始设定的时间步长为2fs，以后每隔30步增加1fs，共进行了300步，最后是11fs。 在最开始的时候，我们发现体系的能量是有个震荡的，然后才是平稳的。 请问： 1. 最初的微小的震荡是不是正常的？ 2. 也就是在第150步，也就是7fs的时候开始能量缓慢的上升，如果在模拟中，我2-7的哪个时间步长比较合适？ 以上两个问题，请大虾们给点建议，谢谢！！ |
回复此楼» 收录本帖的淘帖专辑推荐
 MD分子动力学 | 感悟人生 |
» 猜你喜欢
γ-NiOOH的cif文件求助
已经有0人回复
-
锆合金吸氢材料
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有55人回复
-
国家青年基金祈福
已经有54人回复
-
分子筛催化剂
已经有1人回复
-
青岛大学国家优青团队急招一名26年秋季入学博士（报名截止时间为6月17号下午1点）
已经有8人回复
-
福建理工大学资源化学专业正式获批并启动招生
已经有20人回复
-
求助电化学测试原理，电化学测试技术分析等PPT
已经有0人回复
-
Crystals（IF=2.9）- Special Issue invatation- Functional Materials
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
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