论文摘要
MCTDH(Multiconfiguration time-dependent Hartree)方法是一种有效的含时近似方法,它在研究多原子体系的分子动力学问题时具有很大的优越性。MCTDH方法的精确性和有效性在很大程度上依赖于求解其运动方程所使用的积分算法。因为MCTDH运动方程是非线性一阶耦合微分方程,所以不能使用常规高效的传播子,如分离算符,Chebyshev等,对其进行求解。寻找具有高精度和高稳定性的积分算法对MCTDH方法的发展非常重要。而在经典力学领域广泛应用的辛算法,已做为一种新的传播子引入到精确求解含时薛定谔方程的量子力学方法中,研究证明辛算法是一种稳定而高效的传播子,应用简单而适用面广。我们把辛算法引入到MCTDH方程的求解当中(简称为辛MCTDH方法),并计算了H+H2、H+D2、Cl+H2等反应体系初态选择的反应几率,以测试该方法的计算精度和长时传播的稳定性。结果表明,辛算法适合于对MCTDH方程的求解,辛MCTDH方法比原有MCTDH方法具有更高的精度,在长时传播中也具有稳定性。在本文的最后,我们还简要讨论了如何进一步提高辛MCTDH方法的计算精度和效率。