[논문 리뷰] Uniform error bounds of time-splitting methods for the nonlinear Dirac equation in the nonrelativistic limit regime.

Weizhu Bao, Yongyong Cai외 1명|arXiv (Cornell University)|2019. 06. 25.

Electromagnetic Simulation and Numerical Methods참고 문헌 42인용 수 4

한 줄 요약

이 논문은 비상대론적 극한 영역에서 비선형 디랙 방정식에 대한 리-트로터 ($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.

연구 동기 및 목표

• 비상대론적 극한 영역에서 비선형 디랙 방정식에 대한 시간 분할 방법의 균일한 오차 행동을 분석하는 것.
• 분할 방법이 샤논 정리에 따른 고전적 샘플링 한계를 뛰어넘어 초해상도를 달성할 수 있는지 조사하는 것.
• 작은 매개수 $\varepsilon$에 대해 독립적인 균일한 오차 한계를 확립하는 것. 이는 빛의 속도의 역수와 비례한다.
• 특히 비공진성 시간 스텝 조건에서 향상된 수렴 속도를 달성할 수 있는 조건을 규명하는 것.

제안 방법

• 분석은 외부 자기장 퍼텐셜이 없는 비선형 디랙 방정식에 대한 리-트로터 분할 ($S_1$) 및 스트랑 분할 ($S_2$)에 집중한다.
• 이 방법들은 각 시간 스텝 동안 선형 및 비선형 부분문제로 방정식을 분해하여 순차적 또는 대칭적으로 풀이한다.
• 시간에 따른 $O(\varepsilon^2)$-진동 행동를 보이는 비상대론적 극한 영역에서, $0 < \varepsilon \leq 1$ 조건에서 엄밀한 오차 분석을 수행한다.
• 시간 스텝 $\tau$ 가 수렴에 미치는 영향을 분석하며, 특히 $\varepsilon$에 대한 공진성과 비공진성 $\tau$ 조건을 구분한다.
• 점근 전개와 에너지 추정을 사용하여 $\varepsilon$에 대해 균일한 수렴을 확립하는 이론적 한계를 유도한다.
• 수렴 속도와 초해상도 행동의 이론적 결과를 검증하기 위해 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 리-트로터 및 스트랑 분할 방법은 비선형 디랙 방정식의 비상대론적 극한 영역에서 균일한 오차 한계를 달성할 수 있는가?
• 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]$ 전역에서 균일한가?

주요 결과

• 리-트로터 ($S_1$) 및 스트랑 ($S_2$) 분할 방법은 작은 매개수 $\varepsilon$ 에 대해 시간 스텝 $\tau$ 에 대해 균일한 $1/2$차 수렴을 달성한다.
• 시간 스텝 $\tau$ 가 비공진성일 경우, $S_1$ 은 향상된 균일한 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$ 영역에서 확인한다.