第四讲运动方程的积分算法.pptVIP

第四讲 运动方程的积分算法 黄敏生 积分运动方程的注意点 何为一个好的运动方程积分算法？ 但这并不是很贴切，因为花在积分运动方程上的时间分率（相对于计算相互作用来说）很小，至少对原子或简单分子体系是这样。 计算速度？？ 准确度对较大的时间步长来说更重要，因为所能使用的时间步长越长，单位时间内力的计算量越少。因此，这意味着采用允许用长时间步长的成熟算法是有利的。（然而，保证准确度，不允许发散） 准确度？？ 积分运动方程的注意点 要想获得允许使用长时间步的算法，必须将信息存储在粒子坐标的较高阶导数中。结果是这需要更多的内存。对于一个通常的模拟来说，这并不是一个严重的缺点。因为除非处理很大的体系，存储这些导数所需的内存与即便是在一个通常的工作站上所能获取的总量相比来说仍是很小。 节省内存？？ 能量守恒是一个重要的判据，但实际上需要区分两种能量守恒，即短时间的和长时间的。复杂的高阶算法通常在短时间内（如在几个时间步内）有很好的能量守恒性，然而该方法通常会有所不期望的特征，即时间较长时总能量漂移。（不一定最好） 能量守恒？？ 积分运动方程的注意点 最好有一个算法能同时在短时间和长时间准确地预测所有粒子的轨迹。 不存在 MD模拟所研究的所有体系，体系的轨迹穿过相空间（即对于用由所有粒子坐标和动量所跨过的6N维空间）时敏感地依赖于初始条件。 这意味着两个初始靠得很近的轨迹随时间的演绎将会显著分