[논문 리뷰] A Toolkit For Steady States of Nonlinear Wave Equations: Continuous Time Nesterov and Exponential Time Differencing Schemes

C. B. Ward, N. Whitaker외 2명|arXiv (Cornell University)|2017. 10. 13.

Nonlinear Photonic Systems참고 문헌 14인용 수 5

한 줄 요약

이 논문은 비선형 파동 방정식, 특히 비선형 슈뢰딩거 방정식의 정 steady 상태를 계산하기 위한 새로운 툴킷을 제안한다. 이 툴킷은 연속 시간 Nesterov 가속 기법과 지수 시간 차분(ETD) 및 스펙트럼 재정규화를 조합하여 사용한다. 제안된 방법들—ACTN(가속된 연속 시간 Nesterov)과 ETDV(변동하는 스텝 크기를 가진 지수 시간 차분)—는 기존 표준 방법들에 비해 훨씬 더 빠른 수렴 속도를 보이며, 1D 및 2D 예제에서 반복 수를 최대 50% 감소시킨다. 이는 비대칭 상태 및 자극 상태로의 수렴을 가능하게 한다.

ABSTRACT

Several methods exist for finding ground (as well as excited) states of nonlinear waves equations. In this paper we first introduce two modifications of the so-called accelerated imaginary-time evolution method (AITEM). In our first modification, time integration of the underlying gradient flow is done using exponential time differencing instead of using more standard methods. In the second modification, we present a generalization of the gradient flow model, motivated by the work of Nesterov, as well as that of Candes and collaborators. Additionally, we consider combinations of these methods with the so-called spectral renormalization scheme. Finally, we apply these techniques to the so-called Squared Operator Method, enabling convergence to excited states. Various examples are shown to illustrate the effectiveness of these new schemes, comparing them to standard ones established in the literature. In most cases, we find significant reductions in the number of iterations needed to reach convergence.

연구 동기 및 목표

• 비선형 파동 방정식, 특히 비선형 슈뢰딩거 방정식의 기초 상태 및 자극 상태를 더 빠르고 안정적으로 계산하기 위한 반복적 방법을 개발하는 것.
• 고급 시간 적분 기법을 통합하여 기존의 허구 시간 진화 방법의 수렴 속도와 안정성을 향상시키는 것.
• 제곱 연산자 방법(SOM)을 통해 기존 방법으로는 도달할 수 없는 비대칭 상태 및 자극 상태로의 수렴을 가능하게 하는 것.
• 기존의 AITEM, ITEM 및 스펙트럼 재정규화와 같은 방법들과 비교하여 새로운 방법들의 성능을 다양한 물리적 조건에서 평가하는 것.
• 비선형 광학, 원자 물리학 및 유체역학 분야의 연구자들이 실용적이고 계산 효율성이 높은 툴킷을 활용할 수 있도록 제공하는 것.

제안 방법

• 비선형 파동 방정식의 에너지 기능의 경사 흐름을 해결하기 위해 Nesterov의 가속 경사 방법의 연속 시간 변형인 ACTN을 제안한다.
• 경사 흐름을 통합하기 위해 지수 시간 차분을 사용하는 ETDV를 도입하여, 기존의 명시적 방법에 비해 안정성과 수렴 속도를 향상시킨다.
• ACTN과 ETDV를 스펙트럼 재정규화와 조합하여, 비틀림이 있거나 대칭성이 있는 문제들에 대한 수렴 성능을 향상시킨다.
• 제곱 연산자 방법(SOM)을 적용하여 원래의 정 steady 상태 문제를 비기본 상태, 즉 자극 상태 및 비대칭 상태로 수렴할 수 있는 시스템으로 변환한다.
• ETDV에서 변수 스텝 크기 전략을 사용하여 시간 스텝을 적응적으로 조절하고 수렴 효율성을 향상시킨다.
• 레이일리 몰입도 기반의 수정된 라그랑주 승수 갱신 규칙을 사용하여 진화 과정 중에 화학적 위치를 동적으로 조정한다.

실험 결과

연구 질문

• RQ1 연속 시간 Nesterov 가속 기법이 비선형 슈뢰딩거 방정식의 정 steady 상태에 수렴하기 위해 필요한 반복 수를 크게 줄일 수 있는가?
• RQ2 지수 시간 차분(ETD)이 비선형 파동 정 steady 상태 문제의 수렴 속도와 안정성 측면에서 기존 시간 적분 방법보다 뛰어나게 작용하는가?
• RQ3 ACTN과 ETDV를 스펙트럼 재정규화와 조합하면 다중 Well 잠재함수에서 비대칭 상태 또는 자극 상태로의 수렴 성능이 향상되는가?
• RQ4 AITEM 및 ITEM과 같은 기존 방법들과 비교하여 새로운 방법들은 다양한 물리적 영역에서 반복 수와 안정성 측면에서 어떻게 성능을 발휘하는가?
• RQ5 제곱 연산자 방법(SOM)이 가속 Nesterov 또는 ETD 방법과 결합되었을 때 비기본 상태를 효과적으로 타겟팅할 수 있는가?

주요 결과

• ACTN 방법은 1D 단일 Well 잠재함수에서 기초 상태에 수렴하기 위해 필요한 반복 수를 AITEM의 106회에서 59회로 줄여 44% 감소시켰다.
• 1D 이중 Well 잠재함수에서 유일하게 ETDV 방법만 비대칭 기초 상태에 수렴했으며, 다른 방법들은 실패하거나 매우 느리게 수렴했다.
• 포물선 포획이 있는 2D 비집속 NLS에서 ACTN과 ETDV는 각각 107~114회의 반복으로 수렴했고, AITEM(177회)과 AITEMRe(171회)보다 뛰어났다.
• 주기적 잠재함수를 가진 2D 집속 NLS에서 모든 방법들이 갭 솔리톤에 수렴했지만, ACTN과 ETDV는 AITEM과 AITEMRe보다 더 적은 반복 수를 요구했다.
• ACTN과 제곱 연산자 방법(SOM)의 조합은 1D에서 다중 봉우리 및 비대칭 해를 포함한 자극 상태에 성공적으로 수렴했으며, 이는 기존 방법으로는 도달할 수 없는 상태였다.
• ETDV 방법은 복잡한 잠재함수에서 뛰어난 안정성을 보였고, 이중 Well 케이스에서 비대칭 기초 상태에 수렴한 유일한 방법이었다.