10-5-5总结

今天基本上全天在研究分子动力学模拟，问题出现在系统能量守恒的问题上，理论上能量（动能+势能）应该是守恒的，但是模拟出来能量却在波动，请教了一下专业人士，解答是不应该波动，如果波动也应该是由于计算机的舍入误差因此的，按照这种说法，能量的有效数位起码要十几位- -！，而我的模拟结果只有4位保持不变，一直在这个范围波动。我的理解的原因可能是离散化求速度，位置，动能，势能等引起的，而且系统的边界条件是镜像的，这样会不会对理论值有影响呢？总之，我觉得离散化求速度，位置，动能，势能，肯定会引起问题的。明天问问专业人士看看吧。