[논문 리뷰] High order explicit symplectic integrators for the Discrete Non Linear Schrödinger equation

Jehan Boreux, Timotéo Carletti외 1명|arXiv (Cornell University)|2010. 12. 15.

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

한 줄 요약

이 논문은 요시다의 조합 방법을 기반으로 하여 이산 비선형 슈뢰딩거 방정식(DNLS)에 대한 고차수의 명시적 심플렉틱 적분법 가족을 제시한다. 이는 에너지 및 질량 보존이 뛰어나 장기적인 수치 시뮬레이션을 가능하게 한다. 이 방법은 명시적이고 계산적으로 효율적이며, 차수 10까지 도달할 수 있으며, 비심플렉틱 및 암시적 방법보다 장기간에 걸쳐 정확성과 안정성 면에서 뛰어나다.

ABSTRACT

We propose a family of reliable symplectic integrators adapted to the Discrete Non-Linear Schrödinger equation; based on an idea of Yoshida (H. Yoshida, Construction of higher order symplectic integrators, Physics Letters A, 150, 5,6,7, (1990), pp. 262.) we can construct high order numerical schemes, that result to be explicit methods and thus very fast. The performances of the integrators are discussed, studied as functions of the integration time step and compared with some non symplectic methods.

연구 동기 및 목표

• 분離 가능한 해밀토니안 구조를 갖지 않아 표준 심플렉틱 방법에 저항하는 이산 비선형 슈뢰딩거 방정식(DNLS)에 대해 고차수의 명시적 심플렉틱 적분법을 개발하기 위해.
• 뉴턴 반복이 필요한 고비용의 암시적 심플렉틱 방법(예: 가우스-라다우, 생성함수 기반 방법)의 한계를 극복하기 위해, 이는 대규모 DNLS 시뮬레이션에 부적합하다.
• 브리더 및 음의 온도 영역과 같은 비선형 동역학적 행동을 연구하는 데 필수적인 에너지, 질량, 그리고 켤레 관계 $\overline{p} = q$와 같은 핵심 물리적 보존량을 장기간 유지하기 위해.
• 시간 간격과 적분법 차수를 최적화하여 계산 효율성과 정확성을 균형 잡고, 기계 정밀도 수준의 에너지 보존을 유지하면서 CPU 비용을 최소화하기 위해.
• 비준거적 통합 가능한 성질을 지닌 시스템임에도 불구하고, 요시다의 조합 기법을 통해 유도된 명시적 고차수 심플렉틱 방법이 DNLS에 대해 실현 가능하고 우월하다는 것을 입증하기 위해.

제안 방법

• 요시다의 조합 방법을 변형하여 DNLS에 대해 고차수 심플렉틱 적분법을 구성한다. 이는 2차 심플렉틱 사상의 반복적 조합을 통해 차수 2, 4, 6, 8, 10를 달성한다.
• DNLS 시스템의 해밀토니안을 $ H = -i \sum_{l=1}^{N} \left[ \frac{(p_{l+1}-p_l)(q_{l+1}-q_l)}{h^2} - p_l^2 q_l^2 \right] $로 정의한다. 이는 운동에너지와 위치에너지로 분리되지 않아 표준 분할 방법을 적용할 수 없다.
• 선형 부분과 비선형 부분의 해밀토니안 흐름을 각각 $ L_A $ 및 $ L_B $로 나타내며, $ \exp(L_B) $는 벡터-행렬 곱을 사용하여 정의한다.
• 최적화된 계수 $ (c_j, d_j) $를 가진 $ \exp(L_A) $, $ \exp(L_B) $ 및 그 시간 역행 형태의 조합을 통해 적분법을 구현함으로써 심플렉틱성과 고차수 정확성을 확보한다.
• 요시다의 방법을 개선한 버전($ Y^{\text{opt}}_{2m} $)을 사용하여 조합 단계 수를 줄여 CPU 비용을 감소시키면서도 차수를 유지한다.
• FFT 기반 방법을 활용하여 $ \exp(L_B) $의 푸리에 유사한 구조를 이용해 계산을 가속화함으로써, $ \mathcal{O}(N^2) $의 비용을 $ \mathcal{O}(N \log N) $으로 감소시킨다.

실험 결과

연구 질문

• RQ1 분리 가능한 해밀토니안 구조를 갖지 않아 표준 심플렉틱 분할 방법에 저항하는 DNLS 시스템에 대해 고차수의 명시적 심플렉틱 적분법을 구성할 수 있는가?
• RQ2 장기간의 통합 시간 동안 제안된 심플렉틱 적분법의 에너지 및 질량 보존 성능은 비심플렉틱 및 암시적 방법과 비교하여 어떻게 되는가?
• RQ3 목표 정확도와 에너지 보존을 달성하기 위해 적분법 차수, 시간 간격 크기, 계산 비용(CPU 시간) 사이의 상호 상충 관계는 어떠한가?
• RQ4 만약 초기 조건에서 $ \overline{p(t)} = q(t) $ 관계가 성립한다면, 수치적 방법이 이 켤레 관계를 어느 정도 유지하는가?
• RQ5 제안된 방법은 큰 시간 간격을 사용하여도 기계 정밀도 수준의 에너지 보존을 달성할 수 있는가? 이는 DNLS의 장기적 점점 수렴하는 행동을 효율적으로 시뮬레이션할 수 있게 한다.

주요 결과

• 제안된 고차수 명시적 심플렉틱 적분법($ Y_{2m} $)은 $ Y_4 $의 경우 $ \tau \sim 10^{-3} $, $ Y_6 $의 경우 $ \tau \sim 10^{-2} $의 시간 간격을 사용할 때 총 에너지의 상대 오차가 $ 10^{-15} $ 이하로 유지되어 뛰어난 장기적 에너지 보존 성능을 보여준다.
• 질량 적분 $ I = \sum_l p_l q_l $ 는 거의 기계 정밀도 수준에서 유지되어, 방법이 핵심 물리적 보존량을 유지할 수 있음을 확인한다.
• 초기 조건에서 $ \overline{p(t)} = q(t) $ 관계가 성립한다면, 이는 수치적으로 유지되며, 보즈-아인슈타인 응축체와 같은 물리적 시스템을 모델링하는 데 필수적인 특성이다.
• CPU 시간은 고정된 적분법 차수에 대해 약 $ 1/\tau $ 비례하고, 고정된 $ \tau $에 대해 $ 3^m $ 비례하며, 이는 차수 $ 2m $에 따라 지수적 증가를 의미한다. 그러나 최적화된 $ Y^{\text{opt}}_{2m} $ 형태는 표준 $ Y_{2m} $ 대비 CPU 비용을 $ Y_6 $의 경우 1.5배, $ Y_8 $의 경우 2.2배 감소시킨다.
• 유사한 정확도를 확보할 때, $ Y_6(\tau = 10^{-2}) $는 $ Y_4(\tau = 10^{-3}) $보다 CPU 시간을 덜 소비하며, $ Y_8(\tau = 2 \times 10^{-2}) $는 $ Y_6(\tau = 10^{-2}) $보다 더 효율적이다. 이는 더 높은 차수의 방법이 더 큰 시간 간격에서 더 효율적일 수 있음을 보여준다.
• 수정된 룬게-쿠타 4차 방법은 $ Y_4 $보다 빠르지만 $ Y^{\text{opt}}_6 $보다는 느리며, 표준 RK4는 $ Y_4 $보다 훨씬 비싸다. 이는 효율성과 정확성 측면에서 심플렉틱 접근법의 우월성을 강조한다.