含力梯度显辛算法拓广于摄动二体问题.pdfVIP

第27卷 第4期 江 西 科 学 V01．27No．4 2009年 8月 ． JIANGXI SCIENCE Aug．2009 文章编号：1001—3679(2009)04一o5o6一o5 含力梯度显辛算法拓广于摄动二体问题 徐 佳 (南昌大学物理系，江西 南昌330031) 摘要：将 + 分解形式的哈密顿系统对应的含力梯度显辛算法推广到具有 ％ +日 分解形式(其中 为动 能，为势能，而 日0和 均可积且前者是主要部分)，并将其应用于求解摄动二体问题。能量相对误差表明 风 +日1分解的含梯度显辛算法明显优于相应的 + 分解方法。 关键词：拓广；含力梯度显辛算法；摄动二体问题 中图分类号： 文献标识码：A ExplicitSymplecticIntegratorswithForceGradient Extend tothePerturbedTwo-bodyProblem XUJia (DepartmentofPhysics，NanchangUniversity，JiangxiNanchang330031PRC) Abstract：ThispaperextendsexplicitsymplecticintegratorswithforcegradientscorrespondingHam · ihoniansystemwhichhasdecomposedformasT+Vtodecomposedofmr asHo+H (whereTandV iskineticandpotentialenergy．Boht HoandHlareintegrable，andhteformerishtemainpartofthe Hamihonina function)．Thenusehtemtosolvetheperturbedtwo—bodyproblemandfindhtatexplicit symplecticintergatorswiht rgadientsdecomposedas +H1arebetterthan correspondingdecom— po sedasT+Vbyanalyzingtherelativeenergyerroroftheperturbedtwo—bodyproblem． Keywords：Extension，Explicitsymplecticintergatorswithforcergadient，Theperturbedwto—body problem 者 一主要沿着Ruth思路给予辛算法迅猛的发 0 引言 展，主要用于研究经典力学和量子力学中的多体 辛算法是一种能够保持能量和相空间辛结构 问题，但二阶以上 +I，显辛算法分解因式中都 的数值方法，长时间积分的数值结果能较好的保 存在负系数 问题 J。这将 导致在研究 Monte 持系统的主要特征，迄今为止是研究哈密顿系统 Carlo模型时，数值计算结果出现极大化。为了克 长期定性演化的最佳数值积分工具。 服这一问题，Suzuki_6提出将交杂的多项式指数 20世纪8O年代初，冯康…和 Ruht 分别建 算子引人常见的高阶辛算法，并可以用来求解量 立了隐式辛算法和求解具有动能 和势能 两 子多体问题。由于其分解因式的系数比较复杂， 可积分解形式的哈密顿系统的显辛算法 (简称 71 并未引起广大学者的注意。只有 Chin7【给予了 + 显辛算法)来克服经典算法的不足(本论文中 关注，并以较简