计算化学公社

标题: lammps跑反应，为什么在两种步长（0.1fs和0.25fs），体系的温度差异如此之大（nve）？ [打印本页]

作者
Author: njustwjy 时间: 2022-5-24 14:59
标题: lammps跑反应，为什么在两种步长（0.1fs和0.25fs），体系的温度差异如此之大（nve）？
用Lammps跑反应，前70ps是nvt，体系升到1000K，然后用nve系综，发现在两种步长下（分别是0.1fs和0.25fs），最终的体系温差差异很大，相差了将近2000k，这是为什么呢？我一直以为步长差距不大，跑出来的反应过程应该是差不多的，还希望各位老师帮忙解答一下，学生在此表示感谢。

作者
Author: sobereva 时间: 2022-5-25 05:44
步长越大，动力学积分误差越大，往往导致NVE过程中能量保守性越差

作者
Author: njustwjy 时间: 2022-5-25 11:37
谢谢sob老师

作者
Author: 喵星大佬 时间: 2022-5-25 22:31
https://forum.taichi.graphics/t/ ... ecular-dynamics/992
看看这里面的就知道了

作者
Author: njustwjy 时间: 2022-5-31 10:57

喵星大佬发表于 2022-5-25 22:31
https://forum.taichi.graphics/t/homework-1-bechmarking-time-integrators-for-molecular-dynamics/992
...

感谢，看了以后很有收获。文中提到的time的unit是fs吗？0.004fs也太小了吧

作者
Author: 喵星大佬 时间: 2022-5-31 10:59

njustwjy 发表于 2022-5-31 10:57
感谢，看了以后很有收获。文中提到的time的unit是fs吗？0.004fs也太小了吧

很明显是ps

作者
Author: Yaqi 时间: 2022-5-31 18:13

喵星大佬发表于 2022-5-25 22:31
https://forum.taichi.graphics/t/homework-1-bechmarking-time-integrators-for-molecular-dynamics/992
...

这个网页里面的可视化很好看，您知道是怎么实现的吗？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)