辛算法在准经典轨迹计算中的应用：H＋H2体系  

作　　者：吴韬[1] 沈长圣[1] 居冠之 居宁 

机构地区：南京大学配位化学国家重点实验室,化学系,亚微观固态化学研究所,南京,210093[1] 南京大学配位化学国家重点实验室,化学系,亚微观固态化学研究所,南京,210093 山东大学晶体材料国家重点实验室, 济南,250100[2] Department,of,Mathematics,Duke,University,Durham,NC,27708-0320,USA[3]

出　　处：《中国科学（B辑）》2002年第1期46-50,共5页Science in China(Series B)

基　　金：国家自然科学基金

摘　　要：采用辛算法模拟了H + H2在BKMP势能面上的长时间动力学行为, 并和Runge-Kutta方法计算的结果进行了比较. 结果表明, 采用Runge-Kutta方法模拟长时动力学行为时体系中出现动力学耗散行为, 导致无法真实的模拟体系的物理实际; 而采用辛算法模拟精确度更高, 四阶辛算法误差只有Runge-Kutta法的约百分之一, 六阶辛算法的误差则更小, 体系中也未见动力学耗散行为. 这些结果说明, 在进行长时动力学行为模拟的时候, 在保持体系能量守恒性上, 辛算法应当是一种更好的积分方法.

关 键 词：辛算法 准经典轨迹法 H+H2体系 BKMP势能面 长时动力学行为 动力耗散 

分 类 号：O643.1[理学—物理化学]