| 查看: 1017 | 回复: 3 | |||||
| 当前主题已经存档。 | |||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 gamma337 的 4 个金币 | |||||
[交流]
【求助】分子动力学时间步长
|
|||||
|
下图是我在模拟前对势函数的时间步估计的一个能量输出图，系综为NVE，所示能量为总能（包括时能和动能），开始设定的时间步长为2fs，以后每隔30步增加1fs，共进行了300步，最后是11fs。 在最开始的时候，我们发现体系的能量是有个震荡的，然后才是平稳的。 请问： 1. 最初的微小的震荡是不是正常的？ 2. 也就是在第150步，也就是7fs的时候开始能量缓慢的上升，如果在模拟中，我2-7的哪个时间步长比较合适？ 以上两个问题，请大虾们给点建议，谢谢！！ |
回复此楼» 收录本帖的淘帖专辑推荐
 MD分子动力学 | 感悟人生 |
» 猜你喜欢
求有机合成大神指点三硫酸乙烯酯（CAS：2793408-99-6）的合成路线
已经有10人回复
-
科研工作者面临的几个矛盾，你会怎么选
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有300人回复
-
钯炭、铂炭、钌炭等加氢催化剂
已经有23人回复
-
AgI的FT-IR谱图
已经有1人回复
-
阴离子交换膜电解槽电压衰减原因
已经有0人回复
-
铜的植物缓蚀剂研究
已经有0人回复
-
26年电池方向博士申请
已经有2人回复
-
【急招】 北京工业大学 第二批博士申请-杰青课题组1-2名额，25日截止！！
已经有10人回复
-
1，6-己二胺甲基化合成四甲基己二胺
已经有2人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
10