| 查看: 2964 | 回复: 18 | |||
gasam金虫 (小有名气)
|
[交流]
分子动力学能量守恒的相关问题 已有5人参与
|
||
|
小弟最近在编分子动力学的程序，在测试NVE能量守恒的时候出现偏移，想咨询下各位专家，这个偏移造成的原因？ 自己考虑会不会是因为potential用了cutoff后，能量有误差？希望大家能指点下小弟。。。。。 用了cutoff后能量是会有点下降？还是上升？还是不变？ 如果不变的话，那个偏移会是什么造成的？我确认了力的式子之类的，都没发现错误。。。。 |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 有谁见过这样的动力学问题？
已经有33人回复
- 扩散系数的计算及体系平衡的判断
已经有4人回复
- 【陈懋章】粘性流体动力学基础
已经有249人回复
- MS模拟中能不能进行一些简单的化学变化？
已经有7人回复
- 请教各位，分子动力学模拟NVE系统，不对速度进行标度可以吗？
已经有7人回复
- 第一原理动力学的系综选择，和各种系综的输入文件Example
已经有13人回复
- 分子动力学模拟长程库仑势的处理问题
已经有8人回复
- 第一性原理如何描述化学反应动力学？
已经有21人回复
- 请问你什么是ab initio md
已经有8人回复
- 动力学与热力学的区别
已经有13人回复
- 【求助】siesta 做NOSE动力学出错，请高手指教
已经有11人回复
- 【求助】直接动力学中能量不守恒问题
已经有3人回复
- 【求助】分子动力学中体系温度不断增加?
已经有16人回复
- 【求助】请问NVE系宗的问题
已经有8人回复
- 【求助】模拟金属定向凝固的软件有哪些？
已经有13人回复
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+3, 感谢指导！ 2013-11-22 15:37:54
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+3, 感谢指导！ 2013-11-22 15:37:54
|
偏差略大，按理说数值误差没有这么大。如果你用的是可变步长，可以尝试增加积分精度。你总的时间是多长？如果增加积分精度后还在这个数量级的话，应该就是截断造成的。我用大概25000项求和展开的势能计算2731858步之后的误差大概是你的十分之一左右。不过因为还要积分在个别点附近趋于无穷大的耦合项，我的步长很小，这么多步也只走了50ns，所以需要和你的时间尺度对比下。你可以尝试先不做截断进行计算，排除有bug的可能。 你是不是哪里不小心用单精度数了？ |
13楼2013-11-22 15:19:36
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
2楼2013-11-22 02:15:52
gasam
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1257.2
- 散金: 1
- 红花: 2
- 帖子: 130
- 在线: 140.2小时
- 虫号: 350788
- 注册: 2007-04-21
- 性别: GG
- 专业: 机械摩擦学与表面技术
3楼2013-11-22 07:53:06
sciencejoy
新虫 (著名写手)
- 应助: 436 (硕士)
- 金币: 11207.1
- 红花: 90
- 帖子: 2974
- 在线: 498小时
- 虫号: 802149
- 注册: 2009-07-02
- 性别: GG
- 专业: 高分子物理与高分子物理化
4楼2013-11-22 09:05:35
yaozhq
金虫 (小有名气)
- 应助: 126 (高中生)
- 金币: 1137.8
- 红花: 7
- 帖子: 260
- 在线: 37.6小时
- 虫号: 1475454
- 注册: 2011-11-04
- 性别: GG
- 专业: 生物大分子结构与功能
5楼2013-11-22 09:34:33
gasam
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1257.2
- 散金: 1
- 红花: 2
- 帖子: 130
- 在线: 140.2小时
- 虫号: 350788
- 注册: 2007-04-21
- 性别: GG
- 专业: 机械摩擦学与表面技术
6楼2013-11-22 10:24:39
gasam
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1257.2
- 散金: 1
- 红花: 2
- 帖子: 130
- 在线: 140.2小时
- 虫号: 350788
- 注册: 2007-04-21
- 性别: GG
- 专业: 机械摩擦学与表面技术
7楼2013-11-22 10:25:50
|
啊 [ 发自小木虫客户端 ] |
8楼2013-11-22 10:41:57
sciencejoy
新虫 (著名写手)
- 应助: 436 (硕士)
- 金币: 11207.1
- 红花: 90
- 帖子: 2974
- 在线: 498小时
- 虫号: 802149
- 注册: 2009-07-02
- 性别: GG
- 专业: 高分子物理与高分子物理化
9楼2013-11-22 12:49:31
gasam
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1257.2
- 散金: 1
- 红花: 2
- 帖子: 130
- 在线: 140.2小时
- 虫号: 350788
- 注册: 2007-04-21
- 性别: GG
- 专业: 机械摩擦学与表面技术
10楼2013-11-22 12:57:55
