[論文レビュー] A Note on Symplectic Algorithms

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

Cellular Automata and Applications被引用数 5

ひとこと要約

本稿では、非可換微分計算とEuler–Lagrangeコホロジーを用いて、有限次元ハミルトニアン系に対するラグランジュ形式におけるシンプレクティック・アルゴリズムを提示する。解析により、台形積分スキームが解の上にとどまらず関数空間においてもシンプレクティック構造を保存することを示した。保存されるシンプレクティック2形式が閉形式かつ非退化であることを示し、ハミルトニアン枠組みを超えてシンプレクティック構造保存を拡張した。

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.

研究の動機と目的

• 有限次元ハミルトニアン系に対するラグランジュ形式におけるシンプレクティック・アルゴリズムを確立し、既存のハミルトニアンベースのシンプレクティックスキームを補完すること。
• 時間離散化におけるシンプレクティック構造保存の課題に、解空間から関数空間への概念の拡張を用いて対処すること。
• 台形スキームが解軌道上ではなく関数空間上で一般化された意味でシンプレクティックであることを示すために、関数空間における保存されるシンプレクティック2形式を導出すること。
• 離散時間におけるシンプレクティック構造保存の解析のための厳密な枠組みを提供するため、Euler–Lagrangeコホロジー概念を活用すること。

提案手法

• 時間微分がライブニッツ則を満たさない場合に対処するため、離散時間格子上で非可換微分計算（NCDC）を適用する。
• Euler–Lagrange（EL）1形式およびコホロロジー的概念—ゼロEL1形式、余境界EL1形式、およびEL条件—を導入し、シンプレクティック構造保存の解析に用いる。
• 標準的1形式 θ = (∂L/∂q̇)dq から定義されるシンプレクティック2形式 ω = dθ を定義し、EL条件 dE = 0 の下でその時間微分が消えることを示す。
• DEL1形式を構築し、その結果得られる2形式が閉形式かつ非退化であることを示すことで、台形スキームの保存されるシンプレクティック構造を導出する。
• 変数変換 (p̃, q̃) を用いて、保存されるシンプレクティック形式を dṕ ∧ dq̃ として表現し、そのシンプレクティック性を確認する。
• 分離可能および非分離可能なハミルトニアン系の両方に対してこの手法を適用し、それぞれのケースにおける保存されるシンプレクティック2形式の明示的表現を導出する。

実験結果

リサーチクエスチョン

• RQ1 シンプレクティック構造保存を、ハミルトニアン形式にとどまらず、ラグランジュ形式においても厳密に定義・解析することは可能か？
• RQ2 標準的微分計算が非ライブニッツ性のため失敗する状況下で、離散時間積分スキームのシンプレクティシティをどのように解析できるか？
• RQ3 解軌道上ではなく関数空間上で考えた場合、台形スキームは一般化された意味でシンプレクティックであると見なせるか？
• RQ4 Euler–Lagrangeコホロジーは、離散時間におけるシンプレクティック積分スキームの特徴づけにどのような役割を果たすか？
• RQ5 与えられた数値スキームに対して、保存されるシンプレクティック2形式を明示的に構成し、それが閉形式かつ非退化であることを示すことは可能か？

主な発見

• 台形スキームは、シンプレクティック2形式 ω_D^{(n+1)} = ω_D^{(n)} が閉形式かつ非退化であることを保証することにより、解の上にとどまらず関数空間においてもシンプレクティック構造を保存する。
• 分離可能なハミルトニアン H = T(p) + V(q) の場合、保存されるシンプレクティック2形式は ω^{(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̃ と等価である。
• 台形スキームの保存されるシンプレクティック2形式は、明示的に ω_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 は、シンプレクティック構造保存の必要十分条件であり、コホロジーと数値的シンプレクティシティを結びつける。
• 本手法により、解空間から関数空間への解析の拡張を実現し、離散時間格子上でも厳密な微分計算が可能となることから、シンプレクティック構造保存のラグランジュ形式への一般化に成功した。