We study the spherical pendulum system with an arbitrary potential function $V = V (z)$, which is an integrable system with a first integral whose Hamiltonian flow is periodic. We give an explicit solution to this integrable system and then we compute its action-angle coordinates. In the special case where the potential function is symmetric quadratic like $V = z^2$, we represent its action-angle coordinates in terms of elliptic integrals, and calculate the monodromy.
論文ID : 2509.04207タイトル : Action-angle coordinates of spherical pendulums with symmetric quadratic potentials著者 : Chengle Peng, Xiudi Tang分類 : math.SG(シンプレクティック幾何学)発表日 : 2025年11月7日論文リンク : https://arxiv.org/abs/2509.04207 本論文は、任意のポテンシャル関数 V = V ( z ) V = V(z) V = V ( z ) を持つ球面振子系を研究する。これは第一積分を持つ可積分系であり、そのハミルトン流は周期的である。本論文は、この可積分系の明示解を与え、その作用-角度座標を計算する。ポテンシャル関数が対称二次関数 V = z 2 V = z^2 V = z 2 である特殊な場合、著者は楕円積分を用いて作用-角度座標を表現し、単値性を計算した。
可積分系理論 : ハミルトン動力学において、Liouville可積分系は最大数の関数独立な保存量を持つ重要な動力学系である。その中でも半環面系(semitoric systems)は特殊な部分類であり、一つの成分に回転対称性を示す。球面振子の重要性 : 球面振子は、可積分でありながら焦点-焦点特異点を持つ自然なハミルトン系の数少ない例の一つである。線形ポテンシャル場(例えば地球表面の重力)では、唯一の焦点-焦点点が存在し、これはDuistermaat が作用-角度座標が大域的に存在しないことを発見した古典的な例である。一般的なポテンシャル関数の課題 : より一般的なポテンシャル関数を許容する場合、複数の焦点-焦点点および他の特異点の振る舞いが存在する可能性がある。V ( z ) = z 2 V(z) = z^2 V ( z ) = z 2 の場合、運動量写像の同一ファイバー上に正確に2つの焦点-焦点点が存在する。理論の完成 : 球面振子の位相的性質は研究されているが、その作用-角度座標の明示計算はいまだ不完全である。技術的課題 : 既存の研究の大部分はDirac-Poisson構造と位相的性質に焦点を当てており、一般的な球面振子の完全な分析が欠けている。実用的応用 : 作用-角度座標は量子化と動力学分析において重要な意義を持つ。一般的な球面振子の明示解 : ポテンシャル関数が一般関数 V ( z ) V(z) V ( z ) である球面振子に対して、ハミルトン流 J J J と H H H の明示解(式3.10)を与え、作用座標を定積分として表現した(式5.8)。二次ポテンシャルの完全な分析 : V ( z ) = z 2 V(z) = z^2 V ( z ) = z 2 の場合について深く研究し、楕円積分を用いて結合流の明示公式(式4.3)と作用-角度座標(式5.6)を与えた。単値性の計算 : 焦点-焦点ファイバーの単値性を計算し(式5.7)、可積分系における焦点-焦点ファイバーの理論と一致することを示した。シンプレクティック不変量の明示計算 : 実際にシンプレクティック不変量を計算した数少ない半環面系の例の一つを提供した。球面振子系 ( T ∗ S 2 , ω , F ) (T^*S^2, \omega, F) ( T ∗ S 2 , ω , F ) を研究する。ここで:
T ∗ S 2 T^*S^2 T ∗ S 2 は2次元球面の余接束ω \omega ω は標準的なシンプレクティック構造F = ( J , H ) F = (J, H) F = ( J , H ) は運動量写像であり、J J J は角運動量、H H H は機械的エネルギー目標は系の作用-角度座標を計算することであり、特に対称二次ポテンシャル V ( z ) = z 2 V(z) = z^2 V ( z ) = z 2 の場合を対象とする。
T ∗ S 2 T^*S^2 T ∗ S 2 を T ∗ R 3 T^*\mathbb{R}^3 T ∗ R 3 に埋め込み、位置座標 ( x , y , z ) (x,y,z) ( x , y , z ) と運動量座標 ( u , v , w ) (u,v,w) ( u , v , w ) を使用し、制約条件は以下の通り:
x 2 + y 2 + z 2 = 1 , x u + y v + z w = 0 x^2 + y^2 + z^2 = 1, \quad xu + yv + zw = 0 x 2 + y 2 + z 2 = 1 , xu + y v + z w = 0
3つの座標チャートを定義:
北極座標 U N U_N U N : ( ρ , η , θ , ϕ ) (\rho, \eta, \theta, \phi) ( ρ , η , θ , ϕ ) を使用南極座標 U S U_S U S : ( ρ , η , θ , ϕ ) (\rho, \eta, \theta, \phi) ( ρ , η , θ , ϕ ) を使用赤道座標 U E U_E U E : ( z , w , θ , ϕ ) (z, w, \theta, \phi) ( z , w , θ , ϕ ) を使用J = x v − y u , H = 1 2 ( u 2 + v 2 + w 2 ) + V ( z ) J = xv - yu, \quad H = \frac{1}{2}(u^2 + v^2 + w^2) + V(z) J = xv − y u , H = 2 1 ( u 2 + v 2 + w 2 ) + V ( z )
ここで { J , H } = 0 \{J, H\} = 0 { J , H } = 0 であり、系の可積分性を保証する。
F ( T ∗ S 2 ) = { ( j , h ) ∈ R 2 ∣ j ∈ R , h ≥ j 2 2 } F(T^*S^2) = \left\{(j,h) \in \mathbb{R}^2 \mid j \in \mathbb{R}, h \geq \frac{j^2}{2}\right\} F ( T ∗ S 2 ) = { ( j , h ) ∈ R 2 ∣ j ∈ R , h ≥ 2 j 2 }
北極座標系では、ハミルトン方程式は以下の通り:
d ρ d t = η cos δ , d η d t = − ( ρ η ( 1 − ρ 2 sin 2 δ ) 1 − ρ 2 + ( 1 − ρ 2 ) V ~ N ′ ( ρ ) η ) η cos δ \frac{d\rho}{dt} = \eta \cos \delta, \quad \frac{d\eta}{dt} = -\left(\frac{\rho\eta(1-\rho^2\sin^2\delta)}{1-\rho^2} + \frac{(1-\rho^2)\tilde{V}_N'(\rho)}{\eta}\right)\eta\cos\delta d t d ρ = η cos δ , d t d η = − ( 1 − ρ 2 ρ η ( 1 − ρ 2 s i n 2 δ ) + η ( 1 − ρ 2 ) V ~ N ′ ( ρ ) ) η cos δ
ここで δ = ϕ − θ \delta = \phi - \theta δ = ϕ − θ 、V ~ N ( ρ ) = V ( 1 − ρ 2 ) \tilde{V}_N(\rho) = V(\sqrt{1-\rho^2}) V ~ N ( ρ ) = V ( 1 − ρ 2 ) 。
保存量 ρ η sin δ = j \rho\eta\sin\delta = j ρ η sin δ = j を利用して、方程式を以下のように分離:
d ρ d t = ϵ ( 1 − ρ 2 ) ( 2 ρ 2 ( h − V ~ N ( ρ ) ) − j 2 ) ρ \frac{d\rho}{dt} = \epsilon\frac{\sqrt{(1-\rho^2)(2\rho^2(h-\tilde{V}_N(\rho))-j^2)}}{\rho} d t d ρ = ϵ ρ ( 1 − ρ 2 ) ( 2 ρ 2 ( h − V ~ N ( ρ )) − j 2 )
これは ρ \rho ρ と t t t に関する可分離微分方程式である。
V ( z ) = z 2 V(z) = z^2 V ( z ) = z 2 の場合、パラメータを導入:
k = 1 + h − ( 1 − h ) 2 + 2 j 2 1 + h + ( 1 − h ) 2 + 2 j 2 , n = 1 + h − ( 1 − h ) 2 + 2 j 2 2 k = \sqrt{\frac{1+h-\sqrt{(1-h)^2+2j^2}}{1+h+\sqrt{(1-h)^2+2j^2}}}, \quad n = \frac{1+h-\sqrt{(1-h)^2+2j^2}}{2} k = 1 + h + ( 1 − h ) 2 + 2 j 2 1 + h − ( 1 − h ) 2 + 2 j 2 , n = 2 1 + h − ( 1 − h ) 2 + 2 j 2
ρ ( s , t ) = 1 − n sin 2 γ ( t ) \rho(s,t) = \sqrt{1-n\sin^2\gamma(t)} ρ ( s , t ) = 1 − n sin 2 γ ( t ) θ ( s , t ) = θ 0 + ϵ 2 n j k Π ( γ 0 , n , k ) − ϵ 2 n j k Π ( γ ( t ) , n , k ) + s \theta(s,t) = \theta_0 + \epsilon\sqrt{2n}jk\Pi(\gamma_0,n,k) - \epsilon\sqrt{2n}jk\Pi(\gamma(t),n,k) + s θ ( s , t ) = θ 0 + ϵ 2 n jk Π ( γ 0 , n , k ) − ϵ 2 n jk Π ( γ ( t ) , n , k ) + s
ここで γ ( t ) = am ( F ( γ 0 , k ) − ϵ 2 n t / k , k ) \gamma(t) = \text{am}(F(\gamma_0,k) - \epsilon\sqrt{2n}t/k, k) γ ( t ) = am ( F ( γ 0 , k ) − ϵ 2 n t / k , k ) はヤコビ振幅関数。
これは純粋な数学理論研究であるため、主に以下の方法で検証:
対称性検証 : 系の対称性を利用して周期計算の正確性を検証極限ケース : 特殊なパラメータ値での退化ケースを検証単値性検証 : 計算された単値性が既知の理論と一致することを確認楕円モジュラス: 0 < k ≤ n < 1 0 < k \leq \sqrt{n} < 1 0 < k ≤ n < 1 運動量写像値: ( j , h ) ∈ B r (j,h) \in B_r ( j , h ) ∈ B r (正則値領域) 正則値 ( j , h ) ∈ B r (j,h) \in B_r ( j , h ) ∈ B r に対して、周期格子は:
Λ j , h = ( 1 , 0 ) Z + ( − 2 2 k j n Π ( n , k ) , 2 2 k 1 n K ( k ) ) Z \Lambda_{j,h} = (1,0)\mathbb{Z} + \left(-2\sqrt{2}k\sqrt{\frac{j}{\sqrt{n}}}\Pi(n,k), 2\sqrt{2}k\frac{1}{\sqrt{n}}K(k)\right)\mathbb{Z} Λ j , h = ( 1 , 0 ) Z + ( − 2 2 k n j Π ( n , k ) , 2 2 k n 1 K ( k ) ) Z
A 1 = J A_1 = J A 1 = J A 2 = 2 2 π ( n − 1 n k K ( k ) + n k E ( k ) − ( n − 1 ) ( k 2 − n ) n k Π ( n , k ) ) A_2 = \frac{2\sqrt{2}}{\pi}\left(\frac{n-1}{\sqrt{n}}kK(k) + \frac{\sqrt{n}}{k}E(k) - \frac{(n-1)(k^2-n)}{\sqrt{n}k}\Pi(n,k)\right) A 2 = π 2 2 ( n n − 1 k K ( k ) + k n E ( k ) − n k ( n − 1 ) ( k 2 − n ) Π ( n , k ) )
焦点-焦点値 ( 0 , 1 ) (0,1) ( 0 , 1 ) の周りの単値性は:
dA + ∘ ( dA − ) − 1 = ( 1 0 2 1 ) \text{dA}^+ \circ (\text{dA}^-)^{-1} = \begin{pmatrix} 1 & 0 \\ 2 & 1 \end{pmatrix} dA + ∘ ( dA − ) − 1 = ( 1 2 0 1 )
これは焦点-焦点ファイバーの理論的予測と完全に一致する。
j = 0 , h < 1 j = 0, h < 1 j = 0 , h < 1 : Λ j , h = ( 1 , 0 ) Z + ( 0 , 2 2 K ( h ) ) Z \Lambda_{j,h} = (1,0)\mathbb{Z} + (0, 2\sqrt{2}K(\sqrt{h}))\mathbb{Z} Λ j , h = ( 1 , 0 ) Z + ( 0 , 2 2 K ( h )) Z j = 0 , h > 1 j = 0, h > 1 j = 0 , h > 1 : Λ j , h = ( 1 , 0 ) Z + ( 0 , 2 2 h K ( 1 h ) ) Z \Lambda_{j,h} = (1,0)\mathbb{Z} + (0, \frac{2\sqrt{2}}{\sqrt{h}}K(\frac{1}{\sqrt{h}}))\mathbb{Z} Λ j , h = ( 1 , 0 ) Z + ( 0 , h 2 2 K ( h 1 )) Z Álvaro Pelayo と San Vũ Ngọc による単純半環面系の分類 Joseph Palmer らによる全ての半環面可積分系への分類の拡張 Duistermaat による線形ポテンシャル場での作用-角度座標の大域的非存在性の発見 Efstathiou による V ( z ) = z 2 V(z) = z^2 V ( z ) = z 2 ケースのDirac-Poisson構造の研究 Dullin による球面振子のTaylor級数不変量の計算 Le Floch-Pelayo: 結合角運動量 Alonso-Dullin-Hohloch: 結合スピン振子 完全な解決 : 対称二次ポテンシャル球面振子の完全な作用-角度座標の明示表現を初めて与えた楕円積分表現 : すべての重要な量が楕円積分で表現され、数値計算と理論分析が容易である単値性の確認 : 計算結果が焦点-焦点特異点の一般理論と完全に一致する座標チャートの制限 : 完全な計算には複数の座標チャート間の変換が必要であり、複雑性が増す特殊なポテンシャル関数 : 詳細な分析は二次ポテンシャルに限定され、一般的なポテンシャル関数は積分表現のみ非厳密な半環面系 : J J J が真の写像ではないため、系は厳密には半環面系ではないTaylor級数不変量 : 一般的な球面振子のTaylor級数およびその他の不変量の計算その他のポテンシャル関数 : より多くの特殊なポテンシャル関数への明示計算の拡張量子化への応用 : 作用-角度座標を利用した量子化研究理論的完全性 : ハミルトン方程式から作用-角度座標までの完全な導出チェーンを提供計算の正確性 : すべての公式が明示的であり、検証と応用が容易方法論の革新性 : 楕円積分理論とシンプレクティック幾何学の技巧を巧みに組み合わせた結果の重要性 : シンプレクティック不変量を完全に計算した数少ない半環面系の一つ適用範囲 : 完全な結果は特定の二次ポテンシャル関数にのみ適用可能技術的複雑性 : 楕円関数とシンプレクティック幾何学の深い背景知識が必要実用性の制限 : 純粋数学理論研究として、直接的な応用場面が限定的理論的貢献 : 可積分系理論に重要な明示例を提供方法論的価値 : シンプレクティック幾何学における楕円積分の強力な応用を実証啓発的意義 : 他の可積分系の分析に方法論的指針を提供理論物理学 : 古典力学系の精密な分析数学物理学 : 可積分系とシンプレクティック幾何学の研究数値計算 : 球面振子系の精密な数値計算基盤を提供本論文は可積分系理論、半環面系の分類、楕円積分理論など主要分野をカバーする18篇の重要文献を引用しており、研究に堅実な理論的基礎を提供している。