[论文解读] A Note on Symplectic Algorithms

Hao Yang Guo, Y. Q. Li+1 more|ArXiv.org|Apr 6, 2001

Cellular Automata and Applications被引用 5

一句话总结

本文提出了一种基于非交换微分微积分和欧拉-拉格朗日上同调的拉格朗日形式下的有限维哈密顿系统辛算法。通过证明所保持的辛二维形式在函数空间中是闭且非退化的，该文证明了梯形积分器不仅在解上，而且在函数空间中保持辛结构，从而将辛结构保持的概念从哈密顿框架扩展至更广义的框架。

ABSTRACT

We present the symplectic algorithm in the Lagrangian formalism for the Hamiltonian systems by virtue of the noncommutative differential calculus with respect to the discrete time and the Euler--Lagrange cohomological concepts. We also show that the trapezoidal integrator is symplectic in certain sense.

研究动机与目标

• 建立有限维哈密顿系统在拉格朗日形式下的辛算法，以补充现有的基于哈密顿的辛格式。
• 通过将辛结构保持的概念从解空间扩展至函数空间，解决数值积分中保持辛结构的挑战。
• 通过在函数空间中推导出保持的辛二维形式，证明梯形积分器在广义意义上是辛的。
• 利用欧拉-拉格朗日上同调概念，为离散时间下辛结构保持提供严谨的分析框架。

提出的方法

• 在离散时间格点上应用非交换微分微积分（NCDC），以处理不满足莱布尼茨法则的时间导数。
• 引入欧拉-拉格朗日（EL）一形式及其上同调概念——零EL一形式、上边界EL一形式以及EL条件，以分析辛结构保持性。
• 从标准一形式 θ = (∂L/∂q̇)dq 定义辛二维形式 ω = dθ，并证明在EL条件 dE = 0 下其时间导数为零。
• 通过构造DEL一形式，推导梯形格式的保持辛结构，并证明所得二维形式是闭且非退化的。
• 通过变量变换 (p̃, q̃) 将保持的辛形式表达为 dṕ ∧ dq̃，从而确认其辛性质。
• 将该方法应用于可分与不可分哈密顿系统，分别推导出保持辛二维形式的显式表达式。

实验结果

研究问题

• RQ1 是否可以在拉格朗日形式下而非仅在哈密顿形式下，严格定义并分析辛结构保持？
• RQ2 当标准微分微积分因非莱布尼茨行为而失效时，如何分析离散时间积分器的辛性？
• RQ3 当梯形积分器在函数空间而非仅在解轨迹上考虑时，是否可在广义意义上被认为是辛的？
• RQ4 欧拉-拉格朗日上同调在表征离散时间辛积分器中扮演何种角色？
• RQ5 对于给定的数值格式，是否可以显式构造保持的辛二维形式，并证明其闭且非退化？

主要发现

• 梯形积分器通过确保辛二维形式 ω_D^{(n+1)} = ω_D^{(n)} 是闭且非退化的，在函数空间中保持辛结构，而不仅限于解上。
• 对于可分哈密顿量 H = T(p) + V(q)，保持的辛二维形式为 ω^{(n+1)} = (dp^{(n+1)})^T (1 + h²/4 T_pp^{(n+1)} V_{qq}^{(n+1)}) dq^{(n+1)}。
• 梯形格式的辛结构在变换 p̃^{(n)} = p^{(n)} - h/2 H_q^{(n)}, q̃^{(n)} = q^{(n)} + h/2 H_p^{(n)} 下等价于标准辛形式 dṕ ∧ dq̃。
• 梯形格式保持的辛二维形式显式给出为 ω_D^{(n)} = (dp^{(n)})^T (1 + h²/4 H_pp^{(n)} H_{qq}^{(n)} - h²/4 H_{pq}^{(n)} H_{pq}^{(n)} - H_{pq}^{(n)} H_{pp}^{(n)}) ∧ dq^{(n)} - h²/4 (dq^{(n)})^T H_{qq}^{(n)} H_{pq}^{(n)} ∧ dq^{(n)}。
• EL上同调条件 dE = 0 是辛结构保持的充要条件，将上同调与数值辛性联系起来。
• 该方法通过将分析从解空间提升至函数空间，成功将辛结构保持推广至拉格朗日形式，从而在离散时间格点上实现了严谨的微分微积分。