[논문 리뷰] High order exponential integrators for nonlinear Schrödinger equations with application to rotating Bose-Einstein condensates

Christophe Besse, Guillaume Dujardin외 1명|arXiv (Cornell University)|2015. 07. 02.

Numerical methods for differential equations참고 문헌 3인용 수 6

한 줄 요약

이 논문은 도는 보즈아인스타인 응축체(BEC)를 모델링하는 비선형 슈뢰딩거 방정식을 해결하기 위해 고차수 지수적 적분법—특히 지수적 룬게쿠타(ERK) 및 로즌 방법—을 제안한다. 방정식을 회전 항을 제거하도록 변수를 바꾸어 변환함으로써 저자들은 고차수 정확도와 뛰어난 수렴 성질을 달성하였으며, 특히 낮은 오차 상수를 갖는 장시간 시뮬레이션에서 전통적인 스플릿스텝 방법보다 정확도와 계산 효율성 면에서 뛰어나다.

ABSTRACT

This article deals with the numerical integration in time of nonlinear Schrödinger equations. The main application is the numerical simulation of rotating Bose-Einstein condensates. The authors perform a change of unknown so that the rotation term disappears and they obtain as a result a nonautonomous nonlinear Schrödinger equation. They consider exponential integrators such as exponential Runge--Kutta methods and Lawson methods. They provide an analysis of the order of convergence and some preservation properties of these methods in a simplified setting and they supplement their results with numerical experiments with realistic physical parameters. Moreover, they compare these methods with the classical split-step methods applied to the same problem.

연구 동기 및 목표

• 도는 보즈아인스타인 응축체(BEC) 역학에서 발생하는 비자율 비선형 슈뢰딩거 방정식을 위한 고차수 시간 적분 방법을 개발하기 위해.
• 단순화된 소볼레프 공간 프레임워크에서 지수적 룬게쿠타 및 로즌 방법의 수렴성과 보존 성질을 분석하기 위해.
• 이러한 지수적 방법들이 전통적인 스플릿스텝 방법과 비교하여 정확도, 보존 불변량 유지, 계산 비용 측면에서 어떻게 성능을 보이는지 비교하기 위해.
• 실제 2차원 BEC 시뮬레이션에서 도는 것과 비대칭 포물형 포물형 퍼텐셜을 갖는 경우에 이러한 방법의 효과성을 입증하기 위해.
• 지수적 적분법이 장시간에 걸쳐 스플릿팅 방법에서 흔히 발생하는 공진 현상을 피할 수 있음을 보여주기 위해.

제안 방법

• 도는 BEC 방정식을 변수를 바꾸어 회전 항을 제거함으로써 비자율 비선형 슈뢰딩거 방정식으로 변환한다.
• 변환된 방정식에 지수적 룬게쿠타(ERK) 및 로즌 방법을 적용하며, 선형 부분의 정확한 적분을 위해 선형 연산자의 지수를 활용한다.
• 경계 조건을 피하고 라플라스 연산자의 스펙트럼 계산을 단순화하기 위해 주기적 토러스 도메인 $\mathbb{T}_\delta^d$ 를 사용한다.
• 낮은 오차 상수와 높은 정확도를 확보하기 위해 단계 수가 높은 6차 수준의 가우스형 ERK 및 로즌 방법을 구현한다.
• 특히 2차원 시뮬레이션에서 높은 공간 정확도를 확보하기 위해 허프만 스펙트럼 방법을 사용한 가짜 스펙트럼 공간 이산화를 구현한다.
• 1차원 및 2차원 시험 케이스 전반에 걸쳐 운동량, 질량, 에너지에 해당하는 오차 노름 $\mathcal{E}_{P,h}$, $\mathcal{E}_{M,h}$, 및 $\mathcal{E}_{E,h}$ 를 사용하여 방법들을 비교한다.

실험 결과

연구 질문

• RQ1 고차수 지수적 적분법은 도는 BEC 시뮬레이션에서 전통적인 스플릿스텝 방법보다 뛰어난 정확도를 달성할 수 있는가?
• RQ2 정확히 보존되지 않는다는 점을 감안할 때, 지수적 룬게쿠타 및 로즌 방법은 스플릿팅 방법과 비슷한 정확도로 질량 및 기타 보존 불변량을 유지하는가?
• RQ3 비자율적이고 물리적으로 의미 있는 설정에서 ERK 및 로즌 방법의 오차 상수와 수렴 속도는 스플릿팅 방법과 어떻게 비교되는가?
• RQ4 지수적 적분법은 스플릿팅 방법에서 흔히 발생하는 장시간에 걸친 공진 현상을 피할 수 있는가?
• RQ5 실제 2차원 BEC 문제에 대해 암시적 지수적 방법의 계산 효율성은 명시적 스플릿팅 방법보다 어떻게 되는가?

주요 결과

• 가우스-ERK 방법은 로즌 및 스플릿팅 방법보다 훨씬 낮은 오차 상수를 달성하여 모든 테스트 설정에서 유한 시간 정확도 면에서 뛰어난 성능을 보였다.
• 시간 간격이 약 $10^{-4}$ 수준일 때, 가우스-ERK 방법의 상대적 질량 오차는 정확히 질량을 보존하는 로즌 및 스플릿팅 방법과 유사한 수준이었으며, 이는 질량을 정확히 보존하지는 않음에도 불구하고 성능이 뛰어나다는 것을 의미한다.
• 암시적 성격에도 불구하고 1차원 시뮬레이션에서 가우스-ERK 방법의 계산 시간은 로즌 및 스플릿팅 방법보다 낮았다.
• 장시간 시뮬레이션에서 가우스-ERK 방법은 공진 현상이 없었고, 반면 스플릿팅 방법은 이러한 불안정성을 보였다.
• 실제 매개변수($\beta=1000$, $\Omega=0.9$, $\gamma_x=1.05$, $\gamma_y=0.95$)를 갖는 2차원 시뮬레이션에서 가우스-ERK 방법은 이전 연구에서 관측된 비틀림 격자(dipole lattice)의 구조를 성공적으로 재현하였다.
• 이 방법들은 비대칭 포물형 및 포물형-포물형 4차 퍼텐셜을 포함한 일반적인 포물형 포탄에 대해서도 효과적이었으며, 이는 가짜 스펙트럼 방법이 실패하는 영역에서도 성능을 유지하였다.