[論文レビュー] Pendulum Integration and Elliptic Functions

Pedro L. Garrido, Giovanni Gallavotti|arXiv (Cornell University)|Dec 12, 2008

Advanced Differential Equations and Dynamical Systems参考文献 11被引用数 3

ひとこと要約

本稿は、ヤコビの楕円関数を用いて、不安定平衡点近くの古典的 pendulum の正規化超双曲的座標を厳密に導出する。先行研究の誤りを是正し、超双曲的固定点と楕円的固定点の直接的な対応関係を確立する。主な貢献は、楕円積分およびモジュラーパラメータで表されるハミルトニアンの閉形式表現であり、$ g_0(x') $ を用いたpendulum 動作の正確な統合を可能にする。エネルギーおよび作用変数の明示的べき級数および積表示が得られる。

ABSTRACT

Revisiting canonical integration of the classical pendulum around its unstable equilibrium, normal hyperbolic canonical coordinates are constructed

研究の動機と目的

• 不安定固定点近くのpendulumの正規化正準形を、楕円関数論を用いて再導出すること。
• 特に [6, Appendix 9] における正規化超双曲的座標の構成に関する誤りを是正すること。
• 不安定平衡点近くの超双曲的力学と安定平衡点近くの楕円的力学との間の明確な対応関係を確立すること。
• モジュラーパラメータで表したエネルギーおよび作用変数の明示的べき級数および積表示を提供すること。

提案手法

• 時間発展を引数 $ u = gt/k $ の関数としてパラメータ表示するために、ヤコビの楕円関数 $ \text{sn}, \text{cn}, \text{dn} $ を用いる。
• エネルギーおよび作用変数を表現するために、モジュラー・パラメータ $ x' = e^{-\frac{\tau}{\tau'}} $ を導入し、ここで $ \tau = \frac{\text{K}(h')}{\text{K}(h)} $ である。
• ハミルトニアン $ G(pq) = \frac{d\tilde{\text{U}}}{d(pq)}(0) = g $ を導出し、積 $ pq $ に依存することを示し、正規化超双曲的性質を確認する。
• 生成関数 $ S_s' $ および $ R_s' $ を、$ p' = \sqrt{x_s'} \cos(g_0^{(s)}t) $、$ q' = \sqrt{x_s'} \sin(g_0^{(s)}t) $ を含む無限級数を用いて構成する。
• $ (p,q) = (a_s(x_s')p', a_s(x_s')q') $ という変換を適用し、$ a_s^2(z) = -16I \frac{d}{dz}g_0^{(s)}(z) $ とすることで正準構造を保証する。
• 積表示 $ g_0(x') = g \prod_{n=1}^\infty \left( \frac{1 + x'^n}{1 - x'^n} \right)^2 $ を用いて、$ g_0 $ 及びその対数微分をテータ関数で表現する。

実験結果

リサーチクエスチョン

• RQ1 不安定平衡点近くのpendulumの正規化正準形は、どのように楕円関数論を用いて厳密に導出可能か？
• RQ2 先行研究の[6, Appendix 9] における正規化超双曲的座標の構成に、どのような是正が必要か？
• RQ3 エネルギー $ U $ はモジュラー・パラメータ $ x' $ とどのように関係し、そのべき級数または積表示の構造は何か？
• RQ4 解析接続によって、超双曲的（不安定）と楕円的（安定）の間の正確な変換関係は何か？
• RQ5 楕円関数パラメータ表示から、作用角変数をどのように明示的に構成できるか？

主な発見

• エネルギー $ U = 2g^2I / k^2 $ はモジュラー・パラメータ $ x' $ で表され、$ g_0(x') = \frac{\pi}{2} \frac{g}{h' \text{K}(h)} $ であり、積表示 $ g_0(x') = g \prod_{n=1}^\infty \left( \frac{1 + x'^n}{1 - x'^n} \right)^2 $ を持つ。
• エネルギー $ U(x') $ は $ \frac{d}{dx'}U(x') = g_0(x') \frac{d}{dx'}(x' a(x')^2) $ を満たし、ここで $ a(x')^2 = 8I \frac{d}{dx'}g_0(x') $ である。これにより、エネルギーの変動がモジュラー微分と関連づけられる。
• 安定系ではエネルギーが $ U_s(x_s') = 32Ig_s^2 x_s' \prod_{n=1}^\infty \left( \frac{1 + x_s'^{2n}}{1 + x_s'^{2n-1}} \right)^8 $ で与えられ、$ x_s' = -x' $ であることから、超双曲的系と双対的構造を持つことが示される。
• 安定系の正規形ハミルトニアンは $ \mathcal{U}_s(x) = 32Ig_s^2 W\left( \frac{x}{64Ig_s} \right) $ であり、ここで $ W(z) = z(1 - 2z - 4z^2 - 20z^3 - \cdots) $ であり、既知のモジュラー形式と一致する。
• $ (B,\beta) $ から正準 $ (p,q) $ への変換は、$ p' = \sqrt{x_s'} \cos(g_0^{(s)}t) $、$ q' = \sqrt{x_s'} \sin(g_0^{(s)}t) $ を経て、$ a_s(x_s') $ によるスケーリングを施すことで構成され、ヤコビアンが 1 となるように保証される。
• 関数 $ g_0(x') $ の対数微分は $ 4 \sum_{n=1}^\infty \frac{n x'^{n-1}}{1 - x'^{2n}} $ であり、これは $ \frac{1}{2} \frac{d^2}{dz^2} \log \theta_4(z, x') \big|_{z=0} $ に等しく、テータ関数と関連づけられる。