[논문 리뷰] Split-step Fourier methods for the Gross-Pitaevskii equation

Juha Javanainen, Janne Ruostekoski|arXiv (Cornell University)|2004. 11. 05.

Photonic and Optical Devices참고 문헌 2인용 수 28

한 줄 요약

이 논문은 고체-플라즈마 방정식(Gross-Pitaevskii equation, GPE)에 대한 스플릿-스텝 푸리에 방법이, 비선형 항에 가장 최근에 갱신된 파동 함수를 사용할 경우, 선형 슈뢰딩거 방정식에 대한 동일한 방법과 동일한 시간 스텝 정확도를 달성함을 입증한다. 기호 계산을 통해 저자들은 이 간단한 규칙이 테스트된 모든 스플릿-스텝 방법, 고차수 방법을 포함하여 알고리즘 정확도의 순서를 유지함을 증명하며, 임의의 공간 차원에서 최소한의 스플릿-스텝 방법에 대해 이 규칙이 일반적으로 성립할 것이라 추측한다.

ABSTRACT

We perform a systematic study of the accuracy of split-step Fourier transform methods for the time dependent Gross-Pitaevskii equation using symbolic calculation. Provided the most recent approximation for the wave function is always used in the nonlinear atom-atom interaction potential energy, every split-step algorithm we have tried has the same-order time stepping error for the Gross-Pitaevskii equation and the Schroedinger equation.

연구 동기 및 목표

• 시간에 의존하는 고체-플라즈마 방정식(Gross-Pitaevskii equation, GPE)에 대한 스플릿-스텝 푸리에 방법의 정확도를 체계적으로 분석하는 것.
• 초저온 원자 물리학에서 핵심적인 역할을 하는 GPE의 수치적 구현에서 오랫동안 애매하게 여겨진 문제를 해결하는 것, 즉 명확한 수렴 보장 없이 비공식적이고 검증되지 않은 방법을 사용하는 많은 연구 그룹의 문제를 해결하는 것.
• 선형 슈뢰딩거 방정식에 대해 스플릿-스텝 푸리에 방법이 얻는 유리한 정확도 특성이 비선형 GPE로 확장되었을 때 그대로 유지되는지 확인하는 것.
• 특정 알고리즘에 관계없이 고차수 시간 정확도를 보장하는, 계산적으로 비용이 들지 않는 일반적인 규칙을 규명하는 것.

제안 방법

• 시간 스텝 $ h $ 에 대한 기호 멱급수 전개를 사용하여 GPE의 정확한 시간 진동 연산자와 스플릿-스텝 근사 간의 체계적 비교를 수행하였다.
• 비가환적인 운동에너지 및 위치에너지 연산자의 합의 지수를 전개하기 위해 연산자 대수를 사용하였으며, $ h $ 에 대한 순서별로 교환자 항을 추적하였다.
• 오차 항을 특정 순서까지 상쇄시키기 위해 스플릿-스텝 순서의 계수 $ \alpha_i $ 및 $ \beta_i $ 에 대한 다변수 다항식 방정식을 수립하였다.
• 기호 조작을 위해 Mathematica를 사용하여 지수를 멱급수로 전개하고, 파동 함수에 순차적으로 적용함으로써 스플릿-스텝 알고리즘을 구현하였다.
• 위치 단계와 운동량 단계를 각각 시작 단계로 사용한 $ \mathcal{O}(h^3) $, $ \mathcal{O}(h^4) $, $ \mathcal{O}(h^5) $ 방법을 포함한 여러 스플릿-스텝 스킴을 테스트하였다.
• 1차원 두阶도수 미분을 라플라스 연산자로 대체하여 다차원 GPE에 동일한 방법을 적용함으로써 차원 간 일관성을 확인하였다.

실험 결과

연구 질문

• RQ1 스플릿-스텝 푸리에 방법을 비선형 고체-플라즈마 방정식(Gross-Pitaevskii equation)에 적용할 경우 시간 스텝 오차의 순서는 무엇인가?
• RQ2 동일한 알고리즘을 적용할 때 선형 슈뢰딩거 방정식에 대해 스플릿-스텝 방법이 얻는 정확도 특성이 비선형 GPE로 확장되었을 때 그대로 유지되는가?
• RQ3 GPE의 비선형 항 $ g|\psi|^2 $ 를 다룰 때, 선형 경우와 동일한 정확도 순서를 유지하기 위한 특정 규칙이 존재하는가?
• RQ4 정확도가 단계의 순서(예: 위치 또는 운동량에서 시작하는지) 또는 스플릿-스텝 스킴의 계수 선택에 따라 달라지는가?
• RQ5 관찰된 정확도 유지 특성이 임의의 공간 차원 및 고차수 스플릿-스텝 방법으로 일반화될 수 있는가?

주요 결과

• 스플릿-스텝 푸리에 방법이 GPE에 적용될 경우, 비선형 포텐셜 항 $ g|\psi|^2 $ 에 가장 최근에 이용 가능한 파동 함수를 사용할 경우, 선형 슈뢰딩거 방정식의 경우와 동일한 시간 스텝 오차 순서를 달성한다.
• 테스트된 모든 스플릿-스텝 스킴—$ \mathcal{O}(h^3) $, $ \mathcal{O}(h^4) $, $ \mathcal{O}(h^5) $—에 대해 비선형 항에서 $ c_0 = 0 $, $ |c_1| = 1 $ 이고, $ |\psi|^2 $ 에 가장 최근의 $ \psi $ 를 사용할 경우 $ \mathcal{O}(h^3) $ 정확도를 확보한다.
• 가장 최근의 파동 함수를 사용하는 조건은 정확도를 유지하기 위해 뿐만 아니라, 비선형 케이스에서 선형 케이스와 동일한 오차 순서를 유지하기 위해 필수적임이 테스트된 스킴들에 기반해 나타난다.
• 저자들은 이 규칙—$ |\psi|^2 $ 에 가장 최신의 $ \psi $ 를 사용하는 것—이 모든 최소한의 스플릿-스텝 방법에 대해 차원이나 순서에 관계없이 일반적으로 성립할 것이라 추측한다.
• 이 방법의 정확도는 두 차원 및 세 차원에서도 유지되며, 명시적인 $ \mathcal{O}(h^3) $ 세 지수 스킴에 대한 검증을 통해 확인되었다.
• 이러한 강건성의 근본적인 이유는 아직 알 수 없으며, 노름을 보존하는 비선형 미분방정식에 더 깊은 수학적 구조가 존재할 가능성을 시사한다.