[论文解读] 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

一句话总结

本文分析了非相对论极限 regime 中非线性狄拉克方程的李-特罗特ter（$S_1$）与斯特朗格（$S_2$）分裂方法，其中解在时间上呈现 $O(\varepsilon^2)$ 的振荡波长。令人惊讶的是，这两种方法均表现出超分辨率特性，即使时间步长 $\tau$ 超过 $O(\varepsilon^2)$ 的波长，也能准确捕捉解，且在 $\varepsilon$ 上保持一致的 $1/2$ 阶收敛；在非共振 $\tau$ 条件下，收敛率分别提升至 $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.

研究动机与目标

• 分析非线性狄拉克方程在非相对论极限 regime 中时间分裂方法的均匀误差行为。
• 研究分裂方法是否能够实现超分辨率，从而突破香农定理下的经典采样极限。
• 建立与小参数 $\\varepsilon$ 无关的均匀误差界，其中 $\\varepsilon$ 与光速成反比。
• 确定在何种条件下可实现更高的收敛速率，特别是针对非共振时间步长。

提出的方法

• 分析聚焦于非线性狄拉克方程的李-特罗特ter分裂（$S_1$）与斯特朗格分裂（$S_2$），且不考虑外部磁势。
• 方法将方程分解为线性和非线性子问题，分别在每个时间步长内顺序或对称求解。
• 在非相对论极限 regime 中进行严格误差分析，其中 $0 < \varepsilon \leq 1$，解表现出时间上的 $O(\\varepsilon^2)$ 振荡行为。
• 研究时间步长 $\tau$ 对收敛性的影响，特别区分 $\tau$ 相对于 $\varepsilon$ 的共振与非共振情形。
• 通过渐近展开与能量估计推导理论界，建立关于 $\varepsilon$ 的一致收敛性。
• 进行数值实验以验证理论发现的收敛速率与超分辨率行为。

实验结果

研究问题

• RQ1 李-特罗特ter与斯特朗格分裂方法是否能在非线性狄拉克方程的非相对论极限 regime 中实现与 $\varepsilon$ 无关的均匀误差界？
• RQ2 这些分裂方法是否表现出超分辨率，使得即使时间步长 $\tau \gg \varepsilon^2$，也能准确解析 $O(\varepsilon^2)$ 振荡解？
• RQ3 当 $\tau$ 为非共振时，$S_1$ 与 $S_2$ 的收敛速率如何？与共振情形相比有何差异？
• RQ4 超分辨率现象是否在超越 $S_1$ 与 $S_2$ 的高阶分裂方法中依然成立？
• RQ5 误差界如何依赖于小参数 $\varepsilon$？在 $\varepsilon \in (0,1]$ 范围内是否保持一致？

主要发现

• 李-特罗特ter（$S_1$）与斯特朗格（$S_2$）分裂方法在时间步长 $\tau$ 上对小参数 $\varepsilon$ 均实现一致的 $1/2$ 阶收敛。
• 当时间步长 $\tau$ 为非共振时，$S_1$ 实现改进的均匀一阶收敛率 $O(\tau)$。
• 在非共振 $\tau$ 条件下，$S_2$ 展现出改进的均匀 $3/2$ 阶收敛，即 $O(\tau^{3/2})$。
• 这些方法表现出超分辨率，即使 $\tau \gg \varepsilon^2$，也能准确解析 $O(\varepsilon^2)$ 振荡解，违反经典香农采样极限。
• 超分辨率现象被证明在超越 $S_1$ 与 $S_2$ 的高阶分裂方法中依然成立。
• 数值结果验证了理论误差界，并确认了在各种 $\varepsilon$ 与 $\tau$ 参数范围内的超分辨率有效性。