应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 其他 » 分子动力学能量守恒的相关问题
51/1返回列表
查看: 2927  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gasam

金虫 (小有名气)

[交流] 分子动力学能量守恒的相关问题 已有5人参与

小弟最近在编分子动力学的程序，在测试NVE能量守恒的时候出现偏移，想咨询下各位专家，这个偏移造成的原因？
自己考虑会不会是因为potential用了cutoff后，能量有误差？希望大家能指点下小弟。。。。。  用了cutoff后能量是会有点下降？还是上升？还是不变？
如果不变的话，那个偏移会是什么造成的？我确认了力的式子之类的，都没发现错误。。。。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-21 20:54:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gasam

金虫 (小有名气)
引用回帖:
2楼: Originally posted by virtualzx at 2013-11-22 03:15:52
cutoff肯定会有误差。在相同量纲下力应该等于能量导数的相反数，你用比较高阶的方法求一下数值导数，看是不是和你的力吻合得足够好。这个偏差就会导致能量不守恒

谢谢啊，但是cutoff引起的误差有没有什么规律啊，比如说比真实的能量大或者是小？

[ 发自小木虫客户端 ]
赞一下 回复此楼
3楼2013-11-22 07:53:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

virtualzx

木虫 (著名写手)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+2, 感谢指导！ 2013-11-22 15:33:59
cutoff肯定会有误差。在相同量纲下力应该等于能量导数的相反数，你用比较高阶的方法求一下数值导数，看是不是和你的力吻合得足够好。这个偏差就会导致能量不守恒
赞一下 回复此楼
2楼2013-11-22 02:15:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sciencejoy

新虫 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+1, 鼓励交流！ 2013-11-22 15:34:05
干嘛自己编分子动力学程序？那么多现成的
赞一下 回复此楼
4楼2013-11-22 09:05:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaozhq

金虫 (小有名气)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+1, 鼓励交流！ 2013-11-22 15:34:13
计算本来就有误差 不考虑算法 两个浮点数的运算直接就会丢失很多有效数字 这个要看最后偏差有多大了
赞一下 回复此楼
5楼2013-11-22 09:34:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭