[論文レビュー] Analytic Expressions for Exponentials of Specific Hamiltonian Matrices

C. Baumgarten|arXiv (Cornell University)|Mar 3, 2017

Laser-Matter Interactions and Applications参考文献 4被引用数 3

ひとこと要約

本稿では、ケイリー・ハミルトンの定理とファデーヴ・レヴリエ法を用いて、4×4および6×6ハミルトニアン行列の行列指数関数を解析的に計算する手法を提示する。この手法により、時間に依存する係数を用いた行列のべき乗の線形結合として指数関数を表現することで、線形力学系の時間発展演算子を高速かつシンプレクティック性を保ったまま計算可能となる。特にイオンビーム光学分野において、複数の時間点やビームライン位置での繰り返し評価を著しく高速化できる。

ABSTRACT

Hamiltonian matrices appear in a variety or problems in physics and engineering, mostly related to the time evolution of linear dynamical systems as for instance in ion beam optics. The time evolution is given by symplectic transfer matrices which are the exponentials of the corresponding Hamiltonian matrices. We describe a method to compute analytic formulas for the matrix exponentials of Hamiltonian matrices of dimensions $4 imes 4$ and $6 imes 6$. The method is based on the Cayley-Hamilton theorem and the Faddeev-LeVerrier method to compute the coefficients of the characteristic polynomial. The presented method is extended to the solutions of $2\,n imes 2\,n$-matrices when the roots of the characteristic polynomials are computed numerically. The main advantage of this method is a speedup for cases in which the exponential has to be computed for a number of different points in time or positions along the beamline.

研究の動機と目的

• ハミルトニアン行列、特に4×4および6×6の場合に、高速かつ解析的な行列指数関数の計算手法を開発すること。
• 得られる転送行列のシンプレクティック性を保証し、エネルギーおよびエミッタンスなどの物理的保存量を維持すること。
• 線形ビームダイナミクスにおいて、複数の時間点やビームライン位置での行列指数関数の繰り返し計算を効率的に行えるようにすること。
• 解析的解が困難な場合に備え、2n×2nハミルトニアン行列への一般化を図ること。
• 完全な固有値分解が計算コストが高いため、級数展開やシンプレクティックな対角化に代わる数値的に安定した代替手法を提供すること。

提案手法

• ケイリー・ハミルトンの定理を用いて、行列指数関数を次数2n−1までのハミルトニアン行列のべき乗の線形結合として表現する。
• ファデーヴ・レヴリエ法を用いて、行列の偶数次べき乗のトレースから特徴多項式の係数を計算する。
• ハミルトニアン行列の固有値（±λのペアとして現れる）を用いて、三角関数および多項式項を用いた係数関数x_k(τ)の解析的表現を導出する。
• ゼロ固有値ペアが存在する場合には、最初の2m項に対してx_k(τ) = τ^k / k! の単項式を用いる。ここでmはゼロ固有値ペアの数である。
• 特徴多項式の係数から導かれる同伴行列構造を用いて、係数関数x_k(τ)の微分方程式系を再帰的に解く。
• 4×4および6×6行列に対しては解析的に適用し、2n×2n行列への一般化では必要に応じて数値的固有値求解を用いる。

実験結果

リサーチクエスチョン

• RQ1 4×4または6×6ハミルトニアン行列の行列指数関数を、シンプレクティック性を保証しつつ解析的に計算する方法は何か？
• RQ2 線形ビームダイナミクスにおいて、複数の時間点やビームライン位置で繰り返し行列指数関数を計算する際の最も効率的な方法は何か？
• RQ3 ファデーヴ・レヴリエ法をケイリー・ハミルトンの定理と効果的に組み合わせることで、行列指数関数の係数の閉形式表現を導出できるか？
• RQ4 解析的固有値解が得られない場合に、任意の2n×2nハミルトニアン行列へこの手法をどのように一般化できるか？
• RQ5 固有値が純虚数またはゼロである場合、exp(Fτ)の展開における係数関数x_k(τ)の構造はどのようなものか？

主な発見

• 本手法は、ケイリー・ハミルトンの定理とファデーヴ・レヴリエ法を用いて、4×4および6×6ハミルトニアン行列の行列指数関数の解析的表現を提供する。
• 1組の虚数固有値ωをもつ4×4行列の場合、係数関数にはx₄(τ) = (cos(ωτ)−1)/ω⁴ + τ²/(2ω²) および x₅(τ) = (sin(ωτ)−ωτ)/ω⁵ + τ³/(6ω²) といった項が含まれる。
• 本手法により得られる行列指数関数はシンプレクティック性を保証し、エネルギーおよびエミッタンスなどの物理的保存量が維持される。
• 係数関数x_k(τ)を事前に計算して再利用できるため、複数の時間点やビームライン位置での行列指数関数計算において顕著な高速化が達成できる。
• 2n×2n行列の場合、固有値を数値的に計算し、係数関数を決定するための再帰的方程式系を用いることで一般化が可能となる。
• 係数関数は、非ゼロ固有値からの三角関数項と、ゼロ固有値ペアからの奇数次多項式項の和として構成され、同伴行列系を用いて全解が導出される。