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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 其他 » 分子动力学能量守恒的相关问题
北京石油化工学院2026年研究生招生接收调剂公告
192/2返回列表上一页12
查看: 2985  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gasam

金虫 (小有名气)
引用回帖:
9楼: Originally posted by sciencejoy at 2013-11-22 13:49:31
开源的程序可以自由修改，没有版权问题。难道你不知道lammps, dl_poly, gromacs, ……？...

最关键是老师逼我自己写，说不准用那些开源的

[ 发自小木虫客户端 ]
赞一下 回复此楼
11楼2013-11-22 13:09:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

moyunfei8899

铜虫 (小有名气)
直接用lammps
回复此楼
12楼2013-11-22 13:13:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+3, 感谢指导！ 2013-11-22 15:37:54
引用回帖:
6楼: Originally posted by gasam at 2013-11-21 21:24:39
100000step
温度大概2200K，每个原子大概偏0.18Kcal/mol
...

偏差略大，按理说数值误差没有这么大。如果你用的是可变步长，可以尝试增加积分精度。你总的时间是多长？如果增加积分精度后还在这个数量级的话，应该就是截断造成的。我用大概25000项求和展开的势能计算2731858步之后的误差大概是你的十分之一左右。不过因为还要积分在个别点附近趋于无穷大的耦合项，我的步长很小，这么多步也只走了50ns，所以需要和你的时间尺度对比下。你可以尝试先不做截断进行计算，排除有bug的可能。

你是不是哪里不小心用单精度数了？
赞一下(1人) 回复此楼
13楼2013-11-22 15:19:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gasam

金虫 (小有名气)
引用回帖:
13楼: Originally posted by virtualzx at 2013-11-22 16:19:36
偏差略大，按理说数值误差没有这么大。如果你用的是可变步长，可以尝试增加积分精度。你总的时间是多长？如果增加积分精度后还在这个数量级的话，应该就是截断造成的。我用大概25000项求和展开的势能计算2731858步 ...

我写的是基于bond order，积分步长比较小，0.05fs,10万step,我也觉得误差有点大，我推了下理论值要是再小十分之一就在合理范围内。但是我感觉用了cutoff后nve本身能量就应该不是水平的，应该是有点略微下降，不知道我自己的理解对不
赞一下 回复此楼
14楼2013-11-22 20:32:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+1, 感谢指导！ 2013-11-23 09:51:29
引用回帖:
14楼: Originally posted by gasam at 2013-11-22 07:32:25
我写的是基于bond order，积分步长比较小，0.05fs,10万step,我也觉得误差有点大，我推了下理论值要是再小十分之一就在合理范围内。但是我感觉用了cutoff后nve本身能量就应该不是水平的，应该是有点略微下降，不知道 ...

上升或下降取决于你cutoff了什么项以及分子怎么运动。如果实在需要更高精度只好不做cutoff或者提升cutoff threshold再算一次试试了。
赞一下 回复此楼
15楼2013-11-22 22:06:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sciencejoy

新虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
11楼: Originally posted by gasam at 2013-11-22 13:09:52
最关键是老师逼我自己写，说不准用那些开源的
...

这样的老板很讨厌，自以为是锻炼学生，其实白白浪费时间。能修改开源包里的程序才是真正的锻炼学生。
赞一下 回复此楼
16楼2013-11-23 09:46:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gasam

金虫 (小有名气)
引用回帖:
15楼: Originally posted by virtualzx at 2013-11-22 23:06:44
上升或下降取决于你cutoff了什么项以及分子怎么运动。如果实在需要更高精度只好不做cutoff或者提升cutoff threshold再算一次试试了。...

好的，谢了啊

[ 发自小木虫客户端 ]
回复此楼
17楼2013-11-23 12:18:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaozhq

金虫 (小有名气)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+1, 鼓励交流！ 2013-11-25 18:07:31
月只蓝: 金币+1, 鼓励交流！ 2013-11-26 11:26:58
引用回帖:
17楼: Originally posted by gasam at 2013-11-23 12:18:20
好的，谢了啊
...

想到一种方法 就是先用GNU MP 做一个高精度的极短时间模拟 然后通过比较就能观察出来偏差主要从什么地方产生的了 不过想起来就挺麻烦 而且必须很慢……
赞一下 回复此楼
18楼2013-11-25 17:14:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gasam

金虫 (小有名气)
引用回帖:
18楼: Originally posted by yaozhq at 2013-11-25 18:14:33
想到一种方法 就是先用GNU MP 做一个高精度的极短时间模拟 然后通过比较就能观察出来偏差主要从什么地方产生的了 不过想起来就挺麻烦 而且必须很慢……...

这是什么方法？我相关资料吗

[ 发自小木虫客户端 ]
赞一下 回复此楼
19楼2013-11-26 09:17:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gasam 的主题更新
192/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物与医药求调剂 +5 heguanhua 2026-04-05 6/300 2026-04-05 22:58 by Hdyxbekcb
[考研] 复试调剂 +13 呼呼？~+123456 2026-04-05 13/650 2026-04-05 22:07 by 醉翁wl
[考研] 085602调剂 初试总分335 +8 19123253302 2026-04-05 8/400 2026-04-05 20:44 by 南航~万老师
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 22408 总分320，一篇论文二作，两个国三，求调剂 +3 Leomulufu 2026-04-04 5/250 2026-04-05 19:04 by chongya
[考研] 275求调剂 +16 waltzh 2026-04-01 16/800 2026-04-05 17:14 by Hdyxbekcb
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 求调剂 +7 xzghyuj 2026-04-04 7/350 2026-04-04 22:25 by oooqiao
[考研] 348分环境工程·调剂 +10 吴彦祖24k 2026-04-03 11/550 2026-04-04 14:19 by 无际的草原
[考研] 0856调剂 +8 曲听筠 2026-03-30 8/400 2026-04-04 08:46 by tianyyysss
[考研] 求调剂 +3 wos666 2026-04-03 3/150 2026-04-04 05:16 by gswylq
[考研] 350一志愿北京航空航天大学08500材料科学与工程求调剂 +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by 无际的草原
[考研] 一志愿华东理工大学，080500学硕，317分，求调剂 +13 s1145 2026-03-31 15/750 2026-04-03 11:44 by msi123
[考研] 312 化工或制药调剂 +8 小小墨123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[考研] 一志愿郑大材料工程290求调剂 +20 Youth_ 2026-03-30 20/1000 2026-04-02 14:48 by 5896
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[考研] 348求调剂 +6 吴彦祖24k 2026-04-02 6/300 2026-04-02 14:07 by 给你你注意休息
[考研] 324求调剂 +5 想上学求调 2026-04-01 6/300 2026-04-02 10:16 by sanrepian
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
信息提示
关闭
请填处理意见
关闭