[논문 리뷰] WKB-method for the 1D Schrödinger equation in the semi-classical limit: enhanced phase treatment

Anton Arnold, C. Klein. B. Ujvari|arXiv (Cornell University)|2018. 08. 06.

Quantum Mechanics and Non-Hermitian Physics참고 문헌 25인용 수 3

한 줄 요약

이 논문은 1차원 슈뢰딩거 방정식의 준고전적 극한에서 두 번째 차수 WKB 방법을 향상시켜 오차 분석을 정교화하고, 위상 계산에 수치적 적분 대신 스펙트럼 방법을 도입함으로써, $1/\varepsilon$에 의존하는 위상 통합의 수치 오차를 제거함으로써, 굵은 격자에서도 균일하게 정확하고 고차수 수렴을 달성한다. 이는 작은 $\varepsilon$에서 정확도와 강건성을 크게 향상시킨다.

ABSTRACT

This paper is concerned with the efficient numerical computation of solutions to the 1D stationary Schrödinger equation in the semiclassical limit in the highly oscillatory regime. A previous approach to this problem based on explicitly incorporating the leading terms of the WKB approximation is enhanced in two ways: first a refined error analysis for the method is presented for a not explicitly known WKB phase, and secondly the phase and its derivatives will be computed with spectral methods. The efficiency of the approach is illustrated for several examples.

연구 동기 및 목표

• 표준 수치적 방법이 $\varepsilon \ll 1$일 때 매우 고주기적인 1차원 슈뢰딩거 방정식의 해를 해결하는 데에 비효율적임을 다루기.
• 이전의 WKB 기반 방법에서 정확도를 제한하는 $1/\varepsilon$에 의존하는 수치적 위상 통합 오차를 극복하기.
• 위상과 변환 단계에서의 수치 오차를 포함한 오차 분석을 정교화하여 두 번째 차수 WKB 근사의 점점 더 정확한 점근적 정확도를 향상시키기.
• 저차수 적분(예: 심슨의 법칙) 대신 스펙트럼 방법을 사용하여 위상 계산을 수행함으로써 지수 수렴과 조건 수 문제 감소를 달성하기.
• 스펙트럼 위상 계산이 다양한 $\varepsilon$ 값에 걸쳐 균일하게 정확한 해를 제공할 수 있음을 입증하기, 특히 기계 정밀도 한계에 가까운 경우에도 성립함.

제안 방법

• 원래의 매우 고주기적인 슈뢰딩거 방정식을 해석적 위상과 진폭 분해를 통해 더 매끄러운 시스템으로 변환하기 위해 두 번째 차수 WKB 근사를 적용하기.
• WKB 추측식 $\varphi(x) \sim \exp\left(\frac{1}{\varepsilon}\sum_{p=0}^{\infty} \varepsilon^p \phi_p(x)\right)$을 사용하여 주요 위상과 진폭 보정을 유도하기.
• 고정밀도 위상 함수 $\tilde{\phi}(x) = \int_0^x \left(\sqrt{a(\tau)} - \varepsilon^2 \beta(\tau)\right) d\tau$의 적분을 위해 스펙트럼 콜로케이션(클렌쇼-커티스와 바리센트 보간법을 포함)을 구현하기.
• WKB 근사 오차와 위상에서의 수치적 적분 오차를 모두 고려한 정교한 오차 분석을 도입하여, 스펙트럼 방법을 사용할 경우 후자의 오차가 더 이상 $\mathcal{O}(1/\varepsilon)$가 아니라는 것을 보여주기.
• 위상 제거에 의해 매끄러워진 감소된 변수 $Z(x)$에 대해 표준 유한차분법을 사용하여 굵은 메esh에서 해를 구하기, 위상 제거로 인한 매끄러움을 활용하기.
• 스펙트럼 위상 계산과 저차수 적분(Simpson의 법칙)을 비교하여, 특히 작은 $\varepsilon$에서 조건 수와 수치 오차 측면에서 훨씬 우수한 성능을 보임을 입증하기.

실험 결과

연구 질문

• RQ1 1차원 슈뢰딩거 방정식의 WKB 기반 방법에서 위상 계산의 $1/\varepsilon$에 의존하는 수치 오차를 제거하거나 크게 감소시킬 수 있는가?
• RQ2 수치적 적분 방법의 선택(스펙트럼 vs. 적분)이 준고전적 극한에서 WKB 방법의 수렴 속도와 조건 수에 어떤 영향을 미치는가?
• RQ3 위상과 변환 단계에서의 수치 오차를 포함한 정교한 오차 분석이 더 정확하고 신뢰할 수 있는 수렴 추정치를 도출하는 데 기여하는가?
• RQ4 스펙트럼 방법을 사용한 위상 계산이 다양한 $\varepsilon$ 값에 걸쳐 균일하게 정확한 해를 제공할 수 있는가, 특히 기계 정밀도 한계에 가까운 경우에도 성립하는가?
• RQ5 변형된 변수 $Z(x)$는 고주기적이지만 진폭이 낮은 성질을 가지므로, 추가로 스펙트럼 처리가 가능한가?

주요 결과

• 스펙트럼 위상 계산은 이전의 WKB 방법에서 문제시되었던 $\mathcal{O}(1/\varepsilon)$ 오차 항을 제거하여 균일하게 정확한 해를 가능하게 한다.
• $\varepsilon = 10^{-5}$일 때, 스펙트럼 방법의 오차는 $h \in [10^{-5}, 10^{-2}]$ 범위에서 기계 정밀도 이하로 내려가, 이는 이중 정밀도에서 최적의 정확도를 달성했음을 시사한다.
• WKB 오차의 포화 수준 이상에서 $h$에 대해 두 번째 차수의 수렴을 보이며, 오차 그림의 기울기 분석을 통해 이를 확인하였다.
• 그림 5의 오른쪽 그림(스펙트럼 위상)은 매우 작은 $h$에서 인위적인 오차 증가가 없이, 조건 수가 더 좋고 반올림 오차가 적어진 덕분이다. 반면 심슨의 법칙의 경우 그러한 현상이 나타났다.
• 변형된 변수 $Z$의 오차는 문제적인 $E/\varepsilon$ 항을 포함하지 않는 경계 조건 (3.51)를 만족함으로써, 스펙트럼 방법의 강건성을 확인한다.
• 1차 수준 WKB 방법보다 정확도와 수렴 속도에서 뛰어나며, 두 번째 차수 WKB 오차는 $h \leq 10^{-4} \text{--} 10^{-3}$에서 포화에 도달하는 반면, 스펙트럼 방법은 더 작은 $h$에서도 높은 정확도를 유지한다.