[论文解读] Low Regularity Exponential-Type Integrators for Semilinear Schrödinger Equations

Alexander Ostermann, Katharina Schratz|Repository KITopen (Karlsruhe Institute of Technology)|Mar 24, 2016

Advanced Mathematical Physics Problems参考文献 12被引用 4

一句话总结

本文提出了一种用于半线性薛定谔方程的新颖的一阶指数型积分器，可在解属于 $H^{r+1}$ 时实现 $H^r$ 范数下的收敛性，仅需额外一个导数的正则性——远低于经典方法的要求。该格式采用扭曲变量公式，并通过 $\varphi_1$-函数精确积分刚性项，即使在低正则性解下也能实现最优收敛。

ABSTRACT

We introduce low regularity exponential-type integrators for nonlinear Schrödinger equations for which first-order convergence only requires the boundedness of one additional derivative of the solution. More precisely, we will prove first-order convergence in H$^r$ for solutions in H$^r$$^+$$^1$ (r > d/2) of the derived schemes. This allows us lower regularity assumptions on the data than for instance required for classical splitting or exponential integration schemes. For one dimensional quadratic Schrödinger equations we can even prove first-order convergence without any loss of regularity. Numerical experiments underline the favorable error behavior of the newly introduced exponential-type integrators for low regularity solutions compared to classical splitting and exponential integration schemes.

研究动机与目标

• 解决经典分裂法与指数积分器在半线性薛定谔方程中对高正则性要求的问题。
• 开发一种时间积分器，使得在解的正则性假设最小化时仍能保持一阶收敛性。
• 将正则性要求从经典格式所需的两个额外导数降低至仅一个。
• 通过数值实验展示该方法在非光滑与低正则性解下的优越收敛行为。
• 将指数积分器的适用范围扩展至解光滑性有限的问题，特别是在低正则性区域。

提出的方法

• 该方法基于对扭曲变量 $v(t) = e^{-it\Delta}u(t)$ 应用杜哈梅尔公式，将刚性线性演化与非线性项解耦。
• 通过定义为 $\varphi_1(z) = (e^z - 1)/z$ 的 $\varphi_1$-函数精确处理涉及 $\overline{u}^{2p}$ 的主导非线性项，以保持结构与精度。
• 积分器构造为 $u^{n+1} = e^{i\tau\Delta}\left[u^n - i\mu\tau (u^n)^{p+1} \varphi_1(-2i\tau\Delta)(\overline{u^n})^p\right]$，确保一阶收敛。
• 该方法通过直接精确求解刚性部分，避免对杜哈梅尔公式中的积分进行近似，从而降低局部截断误差。
• 在 $r > d/2$ 条件下，于 $H^r$-范数中进行分析，证明了解具有 $H^{r+1}$ 正则性时的收敛性。
• 通过光滑与非光滑初值的数值验证，与李分裂法、斯特兰格分裂法及经典指数积分器进行比较。

Figure 1: Initial values ( 42 ) normalized in $L^{2}$ for two different values of $\vartheta$ . Left: $\vartheta=\frac{3}{2}$ . Right: $\vartheta=5$ .

实验结果

研究问题

• RQ1 能否为半线性薛定谔方程构造一种一阶指数积分器，其正则性要求仅比解空间 $H^r$ 额外多一个导数？
• RQ2 所提出的积分器是否能在解正则性较低时仍保持一阶收敛性，特别是在经典格式因阶数降低而失效的情况下？
• RQ3 对于非光滑或低正则性初值，该新积分器在精度与收敛速率上相较于经典分裂法与指数积分器表现如何？
• RQ4 该方法能否在非线性项为非整数幂 $p$ 时（如非线性项 $|u|^{2p}u$）实现最优收敛，尤其当非线性项不光滑时？
• RQ5 在小非线性情形（如二次薛定谔方程）下，该方法是否仍保持稳定与收敛？

主要发现

• 所提出的指数型积分器在解属于 $H^{r+1}$ 时，于 $H^r$ 范数下实现一阶收敛，且仅需额外一个导数的正则性，而经典方法需两个导数。
• 在一维二次薛定谔方程（$p=1$）中，该方法在无正则性损失下实现一阶收敛。
• 在非光滑初值的数值实验中，经典指数积分器无法收敛，而李/斯特兰格分裂法出现严重的阶数降低至 1/2 阶。
• 新积分器在非光滑解下仍保持一阶收敛，与理论预测一致，并得到数值结果验证。
• 对于非整数指数 $p$，新积分器在 $p = 1/4$ 时仍保持一阶收敛，而经典分裂法显著恶化。
• 在小非线性情形（$\mu = 0.01$）下，新积分器对光滑与非光滑数据均保持高度准确且一阶收敛，优于经典方法。

Figure 2: Error in the energy space $H^{1}$ of the Lie splitting ( 38 ) (blue, dotted circle), Strang splitting ( 39 ) (magenta, star), classical exponential integrator ( 40 ) (green, dashed) and exponential-type integration scheme ( 25 ) (red, dash dotted). Left picture: $H^{3/2}$ solutions. Right