3.1 数值积分算法 3.1 数值积分算法:分子动力学模拟的“时间之舵”与“力之刻度” 倘若将分子动力学(Molecular Dynamics, MD)比作一场在原子尺度上演的精密交响——那么牛顿第二定律 $ \mathbf{F}i = mi \mathbf{a}i = mi \ddot{\mathbf{r}}i $ 就是它的总谱;系综(Ensemble)理论是舞台的物理法则——规定这场演出必须在恒温、恒压或孤立条件下进行;而数值积分算法,则是指挥家手中那支不可见却无处不在的指挥棒:它不创造力,却决定力如何被“听见”;它不定义运动,却裁定运动如何被“看见”。没有这支指挥棒,再优美的谱子也只是一叠静默的纸;