将辛算法应用于准经典轨线理论,模拟了Ba+HF反应在扩展的 London-Eyring-Polanyi-Sato(LEPS)势能面上的动力学行为。比较了系统在RK4-AMH4积分,四阶辛积分,六阶辛积分下的总能量守恒情况,结果表明六阶辛算法最好地保持了反应体系的能量守恒,且用时间最短,能够很好地节省计算资源。得到的了六阶辛算法下生成物BaH的振动分布高于以前计算的结果,进一步证明了六阶辛算法保持能量守恒最好的特点。
Symplectic integrators are employed in the classical trajectory calculation of the Ba+HF reaction on the extended London-Eyring-Polanyi-Sato (LEPS) potential energy surface (PES). We compared total energy conservation of the Hamiltonian system among the traditional RK4-ADH4 integrator, and symplectic integrators at fourth, sixth order. The calculations show that the sixth-order symplectic integrator has been the most suitable method for the quasiclassical trajectory calculation, and it use the less computational expense and exhibits better energy conservation. The vibrational population of the sixth-order symplectic integrator is bigger than the former calculations.