lammps跑反应，为什么在两种步长（0.1fs和0.25fs），体系的温度差异如此之大（nve）？

njustwjy · 发表于 Post on 2022-5-24 14:59:36

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

sobereva · 发表于 Post on 2022-5-25 05:44:42

步长越大，动力学积分误差越大，往往导致NVE过程中能量保守性越差

njustwjy · 发表于 Post on 2022-5-25 11:37:13

谢谢sob老师

喵星大佬 · 发表于 Post on 2022-5-25 22:31:15

njustwjy · 发表于 Post on 2022-5-31 10:57:25

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

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

喵星大佬 · 发表于 Post on 2022-5-31 10:59:04

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

很明显是ps

Yaqi · 发表于 Post on 2022-5-31 18:13:40

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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