[論文レビュー] Electromagnetic pulse propagation in passive media by path integral methods

Sergei V. Shabanov|ArXiv.org|Dec 15, 2003

Plasmonic and Surface Plasmon Research参考文献 11被引用数 3

ひとこと要約

本稿では、量子理論における経路積分法を用いて、非可逆的で分散的かつ吸収的媒質におけるマクスウェル方程式の時間領域数値解法を提案する。ハミルトニアン形式と擬スペクトル技法（特に高速フーリエ変換）を活用することで、指数的収束、数値的分散なし、正確なガウスの法則の満足、無条件安定性を達成し、1時間ステップあたりの計算コストはO(N log₂N)である。

ABSTRACT

A novel time domain solver of Maxwell's equations in passive (dispersive and absorbing) media is proposed. The method is based on the path integral formalism of quantum theory and entails the use of ({\it i}) the Hamiltonian formalism and ({\it ii}) pseudospectral methods (the fast Fourier transform, in particular) of solving differential equations. In contrast to finite differencing schemes, the path integral based algorithm has no artificial numerical dispersion (dispersive errors), operates at the Nyquist limit (two grid points per shortest wavelength in the wavepacket) and exhibits an exponential convergence as the grid size increases, which, in turn, should lead to a higher accuracy. The Gauss law holds exactly with no extra computational cost. Each time step requires $O(N\log_2 N)$ elementary operations where $N$ is the grid size. It can also be applied to simulations of electromagnetic waves in passive media whose properties are time dependent when conventional stationary (scattering matrix) methods are inapplicable. The stability and accuracy of the algorithm are investigated in detail.

研究の動機と目的

• 有限差分法と比較して、より優れた精度と安定性を有する、非可逆的（分散的および吸収的）媒質におけるマクスウェル方程式の時間領域数値アルゴリズムの開発。
• 計算量子物理学の原則（特に経路積分形式と擬スペクトル法）を、古典的電磁波伝播に応用すること。
• グリッドの細分化に伴う指数的収束の達成、人工的な数値的分散の排除、分散媒質におけるエネルギーの正確な保存（ユニタリティ）。
• 従来の定常的（散乱行列）手法が適用できない時間依存媒質特性のシミュレーションを可能にすること。
• 複雑な媒質における電磁パルス伝播のための、有限差分法や有限要素法の代替として、安定的かつ効率的かつ高精度な代替手法の提供。

提案手法

• 一般化された位相空間における1次時間微分方程式に変換することで、非可逆的媒質におけるマクスウェル方程式をハミルトニアン形式で再定式化する。
• 時間発展演算子に経路積分形式（ネルソンの構成法）を適用し、量子的枠組み内での正確な時間発展を可能にする。
• 高速フーリエ変換（FFT）をコア計算エンジンとして用い、O(N log₂N)の複雑度とグリッド細分化に伴う指数的収束を達成する。
• スプリットステップまたはループフロッグスキームを用いてアルゴリズムを実装し、ハミルトニアン演算子を運動エネルギー的およびポテンシャル的部類に分割して、効率的な時間ステッピングを実現する。
• 境界および界面を、境界条件を直接課す代わりに、不連続な媒質パラメータ（例：導電率や誘電率プロファイル）によってモデル化する。
• ハミルトニアンおよび位相空間形式の構造のおかげで、追加計算コストなしにガウスの法則を正確に満足する。

実験結果

リサーチクエスチョン

• RQ1 量子理論における経路積分法を、非可逆的で分散的および吸収的媒質における時間領域マクスウェル方程式に効果的に適応できるか？
• RQ2 提案された擬スペクトル経路積分アルゴリズムは、従来の有限差分法とは異なり、指数的収束を達成し、数値的分散を排除できるか？
• RQ3 追加の計算コストなしに、分散媒質において正確なガウスの法則の満足とエネルギー保存（ユニタリティ）を維持できるか？
• RQ4 定常的散乱行列手法が適用できない時間依存媒質のシナリオにおいて、このアルゴリズムはどのように性能を発揮するか？
• RQ5 滑らかでない、または不連続な媒質パラメータを有する問題に、FFTに基づく経路積分法を適用する際の計算的および数値的課題は何か？

主な発見

• グリッドサイズの増大に伴い指数的収束を示す。有限差分法の多項式的収束率と比較して顕著に優れている。
• ネルソンの経路積分構成法と擬スペクトル法の適用により、数値的分散が存在せず、ナイキスト周波数限界（最短波長あたり2グリッド点）での正確な波動伝搬が保証される。
• ハミルトニアンおよび位相空間形式の構造のおかげで、追加計算コストなしにガウスの法則が正確に満足される。これは定理8.2で証明済み。
• 多共鳴的ローレンツモデルでは、アルゴリズムがユニタリであることが確認され、分散媒質内での波パッケージエネルギーが完全に保存され、減衰がない。
• 有限差分法や有限要素法とは異なり、条件に依存しない安定性を示す。
• 導電性層におけるほぼ完全な吸収を達成するための必要条件が導出された：σₗ < ω²L/(8πc) かつ σₗ > (n+1)c/(8πL)。これにより、反射および透過が抑制される。