论文部分内容阅读
Pihajoki提出了一种适用于不可分离哈密顿系统的二阶显式相空间扩充的的对称蛙跳积分方法。在这个基础上,在扩充的相空间中如何组合这些变量的方法上提出了一些见解。数据测试表明,连续的置换坐标和动量可以使类似的蛙跳算法得到更加准确的计算结果和最佳的长时间稳定性,也提出了一个新的方法用来构建许多四阶显式相空间扩充的类辛算法,每种算法都是由六个常用的不经过任何排列的蛙跳算法构成的。这种构建模式由四部分组成:置换坐标,没有置换的三个二阶蛙跳格式的积分,置换动量,没有置换的三个二阶蛙跳格式的积分。同理,在扩充的相空间中六阶,八阶,甚至更高阶显式类辛算法都是可以得到的。用一些常见的不可分离的哈密顿系统为例,就像有后牛顿项的无旋转的致密双星系统来证明提出的四阶算法要优于目前已知的算法,也包括了在 Chin系统中的四阶辛算法和四阶显式和隐式中点法混合的辛算法。当给予适当的初始值时,可以发现显式相空间扩充的类辛算法更适合用于不可分离的哈密顿模型中,这些模型包括太阳系系统中常用的带微小扰动的重力系统和具有后牛顿项的自旋的致密双星系统。 之后将之用于更加复杂的哈密顿系统中,并讨论二阶,四阶,八阶的性质,讨论相对论核-壳系统的测地线方程的一些物理学性质。最稳定的圆形轨道的半径在赤道平面只取决于四极矩。给定的扁圆的四极矩会导致两个最稳定的圆形轨道存在,它们的半径大于史瓦西半径。然而,一个扁长的四极矩只对应一个最稳定的圆轨道,其半径小于史瓦西半径。在一般的测地线轨道中,上述的四阶显式相空间空充的对称算法能有效的运用到这个哈密顿系统中,尽管它们相对论中核与壳不能分离开来。借助这样的积分方法,可以快速的计算这个轨道并且不会伴随能量误差的增长。这些核壳之间极矩对测地线方程的影响是由数值估计的。