辛算法在N+O2反应准经典轨线计算中的应用

基于Sayós等通过精确的abinitio计算得到的基态势能面,采用四阶显式辛格式计算了N+O2反应体系的准经典轨线,并与四阶Runge-Kutta法的结果进行了比较.结果显示,四阶Runge-Kutta法不能保持反应体系的能量守恒,计算结果没有真实地描述反应体系的运动;然而,四阶显式辛格式保持了哈密顿系统的辛结构,计算中能很好地保持反应体系的能量守恒,在N+O2反应体系的准经典轨线计算中明显优于四阶Runge-Kutta法.