[논문 리뷰] Convex Splitting Method for the Calculation of Transition States of Energy Functional

Shuting Gu, Xiang Zhou|arXiv (Cornell University)|2016. 10. 23.

Block Copolymer Self-Assembly참고 문헌 26인용 수 3

한 줄 요약

이 논문은 에너지 함수의 색인-1 안장점(전이 상태)을 효율적으로 계산하기 위해 반복 최소화 형식(IMF) 내에서 볼록 분할 방법을 도입한다. 이는 Allen-Cahn 및 Cahn-Hilliard 시스템과 같은 에너지 함수의 전이 상태를 계산하는 데 유용하다. 각 IMF 반복에서 보조 함수를 볼록 분할함으로써, 이 방법은 무조건적인 에너지 안정성을 보장하고 큰 시간 간격을 허용하여 1D 및 2D 모델, 특히 Landau-Brazovskii 함수에 대해 계산 효율을 크게 향상시킨다.

ABSTRACT

Among numerical methods for partial differential equations arising from steepest descent dynamics of energy functionals (e.g., Allen-Cahn and Cahn-Hilliard equations), the convex splitting method is well-known to maintain unconditional energy stability for a large time step size. In this work, we show how to use the convex splitting idea to find transition states, i.e., index-1 saddle points of the same energy functionals. Based on the iterative minimization formulation (IMF) for saddle points (SIAM J. Numer. Anal., vol. 53, p1786, 2015), we introduce the convex splitting method to minimize the auxiliary functional at each cycle of the IMF. We present a general principle of constructing convex splitting forms for these auxiliary functionals and show how to avoid solving nonlinear equations. The new numerical scheme based on the convex splitting method allows for large time step sizes. The new methods are tested for the one dimensional Ginzburg-Landau energy functional in the search of the Allen-Cahn or Cahn-Hilliard types of transition states. We provide the numerical results of transition states for the two dimensional Landau-Brazovskii energy functional for diblock copolymers.

연구 동기 및 목표

• 공간적으로 확장된 시스템에서 에너지 함수의 색인-1 안장점(전이 상태)을 계산하기 위한 효율적인 수치적 방법을 개발하는 것.
• 기울기 흐름 시뮬레이션에서 에너지 안정성으로 잘 알려진 볼록 분할 방법을 반복 최소화 형식(IMF)에 적응시켜 안장점 계산에 활용하는 것.
• IMF 내의 보조 최소화 하위 문제의 해법에서 큰 시간 간격을 허용하여 기존 방법의 안정성 제약를 극복하는 것.
• H^{-1} 및 L^2 계량에서 1D Ginzburg-Landau 및 2D Landau-Brazovskii 에너지 함수에 대해 이 방법의 강건성과 효율성을 입증하는 것.
• 상태의 기울기와 수직인 조건을 부여하여 이동 모드로 인한 허위 수렴 문제와 같은 실용적 과제를 해결하는 것.

제안 방법

• 이 방법은 각 IMF 반복에서 보조 함수의 볼록 분할을 통해 하위 문제를 해결하는 방식으로 반복 최소화 형식(IMF)에 볼록 분할 접근법을 통합한다.
• 보조 함수 $L$를 볼록 및 오목 부분으로 분해할 수 있는 일반 원칙을 제안하여 큰 시간 간격에서도 무조건적인 에너지 안정성을 확보한다.
• 반복 단계에서 비선형 방정식을 풀지 않기 때문에 반복적 선형화된 형식을 사용함으로써 계산 효율성을 향상시킨다.
• 보조 함수에서 사용하는 최소 모드가 공간 기울기 $\partial_x \phi$와 수직임을 보장함으로써, 자명한 이동 모드로의 수렴을 방지한다.
• 정확도와 계산 비용의 균형을 유지하기 위해 하위 문제에 적응형 정지 기준을 적용하여 IMF의 이차 수렴 속도를 유지한다.
• 1D Ginzburg-Landau 및 2D Landau-Brazovskii 모델에 대해 이 방법을 검증하였으며, 수치 결과는 큰 시간 간격에서도 안정적인 수렴을 보여준다.

실험 결과

연구 질문

• RQ1 볼록 분할 방법은 반복 최소화 형식(IMF)에 효과적으로 적응되어 색인-1 안장점을 계산하는 데 사용될 수 있는가?
• RQ2 IMF 내의 보조 함수에 대한 볼록 분할 접근법이 무조건적인 에너지 안정성을 보장하고 큰 시간 간격을 허용하는가?
• RQ3 이 방법은 시스템의 이동 불변성으로 인한 허위 해로의 수렴을 방지하는 데 얼마나 효과적인가?
• RQ4 기존 방법에 비해 IMF 내 볼록 분할을 사용할 경우 전이 상태 계산의 계산 효율성이 얼마나 향상되는가?
• RQ5 이 방법은 다블 블록 고분자에 대한 Landau-Brazovskii 모델과 같은 다른 에너지 함수로 일반화될 수 있는가?

주요 결과

• IMF 프레임워크 내에서 적용된 볼록 분할 방법은 무조건적인 에너지 안정성 시간 적분을 가능하게 하여 안정성 손실 없이 큰 시간 간격을 허용한다.
• 1D Ginzburg-Landau 에너지에 대해, 이 방법은 각각 $L^2$ 및 $H^{-1}$ 계량에서 Allen-Cahn 및 Cahn-Hilliard 유형의 전이 상태를 성공적으로 계산한다.
• 2D Landau-Brazovskii 모델에서는 큰 시간 간격에서도 안정적인 전이 상태로 수렴함을 보여, 고차원에서의 강건성을 입증한다.
• 최소 모드가 상태의 공간 기울기와 수직임을 보장함으로써 자명한 이동 모드로의 수렴을 방지하여 물리적으로 의미 있는 안장점 탐색을 보장한다.
• 광범위한 수치 실험을 통해 이 방법이 이차 수렴 속도를 유지하고 기존 방법에 비해 계산 효율을 크게 향상시킴을 확인하였다.
• 이 접근법은 기존의 기울기 흐름을 위한 에너지 안정적 방법을 안장점 계산에 효율적인 도구로 재사용할 수 있는 일반적인 프레임워크를 제공한다.