中文

含小势的 Dirac 方程长时间动力学时间分裂方法的改进一致误差界

数值分析 2021-12-08 v1 数值分析

摘要

我们建立了含小电磁势的 Dirac 方程长时间动力学时间分裂方法的改进一致误差界,该小势由一个无量纲参数 ε(0,1]\varepsilon\in (0, 1] 刻画,表示势的振幅。我们首先用时间分裂方法对 Dirac 方程进行时间半离散,随后用 Fourier 拟谱方法在空间上全离散。利用二阶时间分裂方法对 Dirac 方程的酉流性质,我们证明了半离散和全离散分别在任意时间 t[0,Tε]t\in[0,T_\varepsilon](其中 Tε=T/εT_\varepsilon = T/\varepsilonT>0T > 0)上的一致误差界 C(t)τ2C(t)\tau^2C(t)(hm+τ2)C(t)(h^m+\tau^2),这些界对 ε(0,1]\varepsilon \in (0, 1] 一致成立,其中 τ\tau 为时间步长,hh 为网格尺寸,m2m\geq 2 依赖于解的正则性,且 C(t)=C0+C1εtC0+C1TC(t) = C_0 + C_1\varepsilon t\le C_0+C_1T 关于 tt 至多线性增长,C00C_0\ge0C1>0C_1>0 为独立于 tthhτ\tauε\varepsilon 的两个常数。随后通过采用正则性补偿振荡(RCO)技术,该技术利用解的正则性控制高频模态、利用相位相消和能量方法控制低频模态,我们建立了半离散和全离散直至长时间 TεT_\varepsilon 的改进一致误差界 O(ετ2)O(\varepsilon\tau^2)O(hm+ετ2)O(h^m +\varepsilon\tau^2)。我们给出了数值结果以验证误差界并表明其为尖锐的。同时还提供了 Dirac 方程不同时间离散化精度的比较。

关键词

引用

@article{arxiv.2112.03616,
  title  = {Improved uniform error bounds on time-splitting methods for the long-time dynamics of the Dirac equation with small potentials},
  author = {Weizhu Bao and Yue Feng and Jia Yin},
  journal= {arXiv preprint arXiv:2112.03616},
  year   = {2021}
}

备注

22 pages, 6 figures