hamilon函数分解与hamilon时的混合辛积分器嵌入法.docxVIP

hamilon函数分解与hamilon时的混合辛积分器嵌入法

1.其他经典应用

辛积分器具有保能量和辛结结构的优点。这是一种解决方案，用于研究哈mile系统的数值积分。主要用于研究经典力学、分子动力学、力学、声波、离子和太阳能动态等领域的哈mile定性发展。自30年代以来，已经在这些领域发表了一些关于这种方法的理论和应用研究。

尽管如此,辛积分器在广义相对论中的讨论和运用却非常少见.根源在于广义相对论引力系统大都用度规张量表示,似乎很难找到一个通常意义的Hamilton函数与之对应.更重要的原因在于相对论引力系统一般不能分解成两个可积部分,从而导致显辛方法无法应用,只能采用隐辛方法,这样会降低积分效率.事实上,度规张量对应了一个超Hamilton,这从理论上解决了相对论引力系统应用辛积分器的困难.当然,这种辛格式只能是隐式的,如二阶隐式中点法作为一种辛方法就适合使用.当将度规张量做后牛顿近似展开并采用正则坐标和动量时,可以得到后牛顿Hamilton构型.它由牛顿项和后牛顿项两部分构成,前者属于Kepler主要部分,而后者为摄动部分,这就是所谓的Hamilton摄动分解方法.文献把描述银河系中心的Kerr度规正是如此处理的.然后,该文用分析方法给出牛顿项的解析解,又借助隐中点法求出后牛顿项的数值解,再将它们对称组合构成一个二阶混合辛积分器,后面称为隐中点嵌入法.就作者所知,这种嵌套法早在10年前已被Liao讨论,只不过它当时是为经典Hamilton问题的求解而设计.另一方面,也可以不采用嵌入方式而直接利用隐中点法求解整个后牛顿Hamilton系统,但近期研究表明嵌入法与直接法计算效率相当,但嵌入法拥有数值精度优势.这主要是由于采用Hamilton摄动分解形式导致两部分之间存有后牛顿阶小量级的差别,小量级会出现在截断误差里而引起误差减小.这尤其在Wisdom-Holman二阶leapfrog辛积分器计算太阳系n体问题时变得更为明显.

除参量化后牛顿(PPN)方法外,由中子星或黑洞构成的相对论自旋致密双星二体动力学也是近10多年来后牛顿天体力学和相对论天体物理受到特别关注的几个问题之一.因为自旋致密双星是LIGO和VIRGO等引力波探测器最为理想的波源选择.数据分析处理需要用到高精度引力波形的理论模板,而这种高精度主要由高阶后牛顿近似提供.基于这一点,自旋致密双星动力学从两方面开展研究.一方面,尽可能建立高阶后牛顿运动方程与引力波形,如谐和坐标下的后牛顿Lagrangian构型[19—21]以及ADM坐标下的后牛顿Hamilton构型[21—24].轨道部分现已取得3.5阶后牛顿精度,而自旋部分精确到4阶.另一方面,致力于自旋致密双星是否混沌的研究,因为混沌存在不利于数据分析时剔除信号中的噪声.应当着重指出的是尽管在耗散系统等经典动力学问题中取得了有关混沌理论研究与应用研究的丰硕成果[28—33],但关于自旋致密双星的混沌研究还很不成熟,并且所得结果曾经在文献中一度非常混乱,甚至相互矛盾.例如,Levin利用分形方法发现一体自旋且质量相近或二体自旋且自旋效果仅限于自旋与轨道耦合的2阶后牛顿Lagrangian系统存在混沌,但文献却指出一体自旋且质量相近的2阶后牛顿Hamiltonian是不可能出现混沌的.如果无视两种近似描述构型的差异来看待这些结果,那么自然觉得结论相悖.Levin认为Lagrangian系统与Hamiltonian系统只是近似但不完全等价,因此两者有不同的动力学现象都是正确的.这是由动力学方程不同而引起混沌结果的差异.不同的混沌识别方法的使用也曾带来混沌结果的极大争议.起初,Levin利用分形方法发现2阶后牛顿Lagrangian系统存在混沌;可是,文献用Lyapunov指数方法没有找到混沌;紧接着文献又重新计算Lyapunov指数并且获得正Lyapunov指数作为有力证据来驳斥前文,认为前文没有获得正确Lyapunov指数是由重整化不当引起的.Wu等不持这种观点,而认为他们所采用的Lyapunov指数计算方法是有差异的,前者采用了原始定义求极限方法以便尽可能取得稳定的Lyapunov指数,需要相当长的积分时间才能完成,但后者是借助拟合方法来计算Lyapunov指数,只要较短的积分时间即可.进一步指出自旋致密双星大都会出现并合(merger)且混沌即使存在也会很弱,因而Lyapunov指数不是识别混沌的最佳工具,至于其他混沌指标如0-1指标(0-1test)是否适合还有待探讨,但两邻近轨道的快速Lyapunov指标是适合的.Wu等还澄清了文献关于初始条件或参数变化而引起动力学跃迁的争议,认为对于其他固定的初始条件和参数的不同选择来说各种跃迁结果都是可能的,并建议用初始条件和参数的全面综合来看待这种跃迁关系.总