[논문 리뷰] Convergence of a Normalized Gradient Algorithm for Computing Ground States

Erwan Faou, Tiphaine Jézéquel|arXiv (Cornell University)|2016. 03. 08.

Advanced Mathematical Physics Problems참고 문헌 14인용 수 19

한 줄 요약

이 논문은 1차원 입자형 비선형 슈뢰딩거 방정식(NLS)의 기본 상태를 계산하기 위한 정규화된 기울기 알고리즘(가상 시간 방법으로 알려진)의 수렴성을 증명한다. 유한 차분 공간 이산화와 선형적 암시적 시간 적분기를 사용하여, 이 방법은 공간 메쉬 크기 $h$와 도메인 절단 $K$에 따라 변하는 오차 한계를 갖는 이산 솔리톤으로 지수적으로 수렴한다. 주요 기여는 $H^1$ 노름에서 명시적인 오차 추정을 포함한 엄밀한 수렴 분석이다.

ABSTRACT

We consider the approximation of the ground state of the one-dimensional cubic nonlinear Schr{\\"o}dinger equation by a normalized gradient algorithm combined with linearly implicit time integrator, and finite difference space approximation. We show that this method, also called imaginary time evolution method in the physics literature, is con-vergent, and we provide error estimates: the algorithm converges exponentially towards a modified solitons that is a space discretization of the exact soliton, with error estimates depending on the discretization parameters.

연구 동기 및 목표

• 1차원 입자형 비선형 슈뢰딩거 방정식의 기본 상태를 계산하기 위한 정규화된 기울기(가상 시간) 방법의 수렴성을 엄밀히 확립하는 것.
• 공간 및 시간 이산화가 존재하는 상황에서 계산된 해와 정확한 솔리톤 사이의 오차를 분석하는 것.
• 선형적 암시적 시간 스킴이 이산 기본 상태를 정확히 유지함으로써 안정성과 정확성을 보장하는지 확인하는 것.
• 메쉬 크기 $h$, 시간 단계 $\tau$, 도메인 절단 $K$에 대한 의존성을 포함한 완전히 이산화된 알고리즘에 대한 정량적 오차 추정을 제공하는 것.

제안 방법

• 알고리즘은 에너지 함수의 $L^2$-기울기의 음수에 해당하는 편미분 방정식 $\partial_t \psi = \frac{1}{2}\Delta\psi + |\psi|^2\psi를 해결하기 위해 선형적 암시적 시간 적분기를 사용한다.
• 스킴은 $\psi_n^* = \psi_n - \tau \widehat{\nabla H}(\psi_n, \psi_n^*)$로 정의되며, 여기서 $-\widehat{\nabla H}(\psi_n, \psi_n^*) = \frac{1}{2}\Delta_h \psi_n^* + |\psi_n|^2 \psi_n^*$로 정의되어, 무조건적 안정성과 기본 상태의 정확한 유지 보장을 보장한다.
• 각 시간 단계 이후, $L^2$ 노름에 따라 $\psi_{n+1} = \psi_n^* / \|\psi_n^*\|_{L^2}$로 정규화하여 단위 질량을 유지한다.
• 공간 이산화에서는 $|x| \leq Kh$ 영역에서 중심 유한 차분을 사용하며, 경계 조건으로 0을 적용하고, 수렴을 측정하기 위해 이산 $H^1$ 노름을 사용한다.
• 이산 기본 상태 $\eta_{h,K}$는 단위 $L^2$ 노름 하에서 에너지 함수의 유일한 최소화자로 정의되며, 이산 고유값 방정식을 만족한다.
• 안정성과 볼록성 분석을 위해 $\eta_{h,K}$ 주변에 국소 좌표계를 구성하여 지수 수렴의 증명을 가능하게 한다.

실험 결과

연구 질문

• RQ1 선형적 암시적 시간 적분기를 사용하는 정규화된 기울기 알고리즘이 공간 및 시간 이산화 하에 1차원 입자형 NLS 방정식의 진정한 기본 상태로 수렴하는가?
• RQ2 이 알고리즘의 수렴 속도는 무엇이며, 메쉬 크기 $h$, 시간 단계 $\tau$, 도메인 절단 $K$와의 의존성은 어떻게 되는가?
• RQ3 선형적 암시적 스킴은 이산 기본 상태를 정확히 유지하는가? 이 성질은 수렴성과 안정성에 어떤 영향을 미치는가?
• RQ4 이산 솔리톤 $\eta_{h,K}$와 정확한 연속 솔리톤 $\eta$ 사이의 거리가 $H^1$ 노름에서 얼마나 가까운가?
• RQ5 반면, 반명시적 또는 완전 암시적 스킴은 수렴하는가? 이산 설정에서의 오차 추정은 무엇인가?

주요 결과

• 알고리즘은 시간에 따라 지수적으로 이산 기본 상태 $\eta_{h,K}$로 수렴하며, $\|\psi_n - \eta_{h,K}\|_h \leq C e^{-cn\tau}$를 만족한다. 여기서 $C, c > 0$는 $h$, $K$, $\tau$와 무관한 상수이다.
• 이산 기본 상태 $\eta_{h,K}$는 이산 고유값 방정식 $\frac{1}{2}(\Delta_h \eta_{h,K})^\ell + |\eta_{h,K}^\ell|^2 \eta_{h,K}^\ell = \lambda_h \eta_{h,K}^\ell$을 만족하여, 스킴이 이를 정확히 유지함을 보장한다.
• 연속 솔리톤 $\eta$와 이산 솔리톤 $\eta_{h,K}$ 사이의 오차는 $\|i_h \eta_{h,K} - \eta\|_{H^1} \leq C(h + \frac{1}{h^2} e^{-C_1 h K})$로 유계이며, 이는 $K$에 대해 스펙트럴 수렴, $h$에 대해 대수적 수렴임을 보여준다.
• 완전히 이산화된 해 $\psi_n$는 $\|i_h \psi_n - \eta\|_{H^1} \leq C(e^{-cn\tau} + h + \frac{1}{h^2} e^{-C_1 h K})$를 만족하며, 시간적 지수 감쇠와 공간 이산화 오차를 통합한다.
• 선형적 암시적 스킴은 기본 상태를 정확히 유지한다. 이 성질은 반명시적 또는 완전 암시적 대안 스킴에서는 공유되지 않으며, 이는 수치적 안정성과 정확성을 향상시킨다.
• 반명시적 또는 완전 암시적 스킴의 경우, 해는 수정된 이산 기본 상태 $\eta_{\tau,h,K}$로 수렴하며, 오차는 $\|i_h \eta_{\tau,h,K} - \eta\|_{H^1} \leq C(\tau + h + \frac{1}{h^2} e^{-C_1 h K})$로 유계이다. 이는 시간 이산화가 추가로 $\tau$-의존 오차를 도입함을 보여준다.