[論文レビュー] Uniform error bounds of time-splitting methods for the nonlinear Dirac equation in the nonrelativistic limit regime.

Weizhu Bao, Yongyong Cai+1 more|arXiv (Cornell University)|Jun 25, 2019

Electromagnetic Simulation and Numerical Methods参考文献 42被引用数 4

ひとこと要約

本稿は、非相対論的極限領域における非線形ディラック方程式に対して、Lie-Trotter（$S_1$）およびStrang（$S_2$）分裂法の解析を行う。ここで解は $O(\varepsilon^2)$ の時間波長で振動する。驚くべきことに、両手法とも、時間ステップ $\tau$ が $O(\varepsilon^2)$ 波長を上回る場合でも、解を正確に捉える超解像性を示し、$\varepsilon$ に対して一様に $1/2$ 階の収束を示す。非共鳴時刻ステップでは、それぞれ $O(\tau)$ および $O(\tau^{3/2})$ の改善された収束率を示す。

ABSTRACT

Super-resolution of the Lie-Trotter splitting ($S_1$) and Strang splitting ($S_2$) is rigorously analyzed for the nonlinear Dirac equation without external magnetic potentials in the nonrelativistic limit regime with a small parameter $0<\varepsilon\leq 1$ inversely proportional to the speed of light. In this limit regime, the solution highly oscillates in time with wavelength at $O(\varepsilon^2)$ in time. The splitting methods surprisingly show super-resolution, in the sense of breaking the resolution constraint under the Shannon's sampling theorem, i.e. the methods can capture the solution accurately even if the time step size $ au$ is much larger than the sampled wavelength at $O(\varepsilon^2)$. Similar to the linear case, $S_1$ and $S_2$ both exhibit $1/2$ order convergence uniformly with respect to $\varepsilon$. Moreover, if $ au$ is non-resonant, i.e. $ au$ is away from certain region determined by $\varepsilon$, $S_1$ would yield an improved uniform first order $O( au)$ error bound, while $S_2$ would give improved uniform $3/2$ order convergence. Numerical results are reported to confirm these rigorous results. Furthermore, we note that super-resolution is still valid for higher order splitting methods.

研究の動機と目的

• 非線形ディラック方程式の非相対論的極限領域における時間分裂法の均一誤差挙動を分析すること。
• 分裂法が、シャノンの定理に従う古典的サンプリング限界を打ち破る超解像性を達成できるかどうかを調査すること。
• $\varepsilon$ に依存しない均一な誤差境界を確立すること。ここで $\varepsilon$ は光速に逆比例する小さなパラメータである。
• 特に非共鳴時刻ステップの下で、改善された収束率が達成される条件を特定すること。

提案手法

• 分析は、外部磁場ポテンシャルのない非線形ディラック方程式に対するLie-Trotter分裂（$S_1$）およびStrang分裂（$S_2$）に焦点を当てる。
• 手法は、各時間ステップごとに線形および非線形部分問題に分解し、逐次的または対称的に解く。
• 時間に $O(\varepsilon^2)$ の振動を示す解を伴う非相対論的極限領域（$0 < \varepsilon \leq 1$）において、厳密な誤差解析を実施する。
• 時間ステップ $\tau$ が $\varepsilon$ に対して共鳴的か非共鳴的かに応じた収束への影響を検討する。
• 漸近展開およびエネルギー推定を用いて理論的境界を導出し、$\varepsilon$ に対して一様な収束性を確立する。
• 理論的結果の妥当性を検証するため、数値実験を実施し、収束率および超解像性挙動を評価する。

実験結果

リサーチクエスチョン

• RQ1 Lie-TrotterおよびStrang分裂法は、非線形ディラック方程式の非相対論的極限領域において、$\varepsilon$ に対して一様な誤差境界を達成できるか？
• RQ2 これらの分裂法は、$\tau \gg \varepsilon^2$ であっても、$O(\varepsilon^2)$-振動解を正確に解像できる超解像性を示すか？
• RQ3 非共鳴時刻ステップの下で、$S_1$ および $S_2$ の収束率は何か？また、共鳴時のケースと比較するとどうなるか？
• RQ4 $S_1$ および $S_2$ を超える高次分裂法に対しても、超解像性は保持されるか？
• RQ5 誤差境界は小さなパラメータ $\varepsilon$ に依存するか？また、$\varepsilon \in (0,1]$ 全体にわたって一様であるか？

主な発見

• Lie-Trotter（$S_1$）およびStrang（$S_2$）分裂法は、すべての $\varepsilon$ に対して一様に $1/2$ 階の収束を示し、時間ステップ $\tau$ に対して一様に成立する。
• 非共鳴時刻ステップの下では、$S_1$ は改善された一様1次収束率 $O(\tau)$ を達成する。
• 非共鳴時刻ステップの下では、$S_2$ は改善された一様 $3/2$ 階収束を示し、すなわち $O(\tau^{3/2})$ である。
• これらの手法は超解像性を示し、$\tau \gg \varepsilon^2$ であっても、$O(\varepsilon^2)$-振動解を正確に解像できる。これは古典的シャノンのサンプリング限界に反する。
• 超解像現象は、$S_1$ および $S_2$ を超える高次分裂法に対しても継続することが示された。
• 数値結果により、理論的誤差境界および超解像性の有効性が、さまざまな $\varepsilon$ および $\tau$ の領域で確認された。