[논문 리뷰] A Deflation Technique for Detecting Multiple Liquid Crystal Equilibrium States

James H. Adler, D. B. Emerson외 2명|arXiv (Cornell University)|2016. 01. 27.

Liquid Crystal Research Advancements참고 문헌 29인용 수 3

한 줄 요약

이 논문은 프랭크-후스텐 자유에너지 모델을 사용하여 난류 및 콜레스테릭 액정에서 다수의 서로 다른 평형 구성을 체계적으로 탐지하기 위해 내재된 반복법과 다중 격자 해법기법과 통합된 탈락 기법을 소개한다. 이 방법은 이전에 발견된 해를 제거하면서도 유한요소의 희박성과 기존의 다중 격자 예비조건화 기법을 유지함으로써 국소 및 전역 최소화 해를 효율적으로 탐색하며, 복잡한 결함 구조와 나선형 상을 포함한 다양한 해를 발견할 수 있다.

ABSTRACT

Multiple equilibrium states arise in many physical systems, including various types of liquid crystal structures. Having the ability to reliably compute such states enables more accurate physical analysis and understanding of experimental behavior. This paper adapts and extends a deflation technique for the computation of multiple distinct solutions arising in the context of modeling equilibrium configurations of nematic and cholesteric liquid crystals. The deflation method is applied as part of an overall free-energy variational approach and is modified to fit the framework of optimization of a functional with pointwise constraints. It is shown that multigrid methods designed for the undeflated systems may be applied to efficiently solve the linear systems arising in the application of deflation. For the numerical algorithm, the deflation approach is interwoven with nested iteration, creating a dynamic and efficient method that further enables the discovery of distinct solutions. Finally, four numerical experiments are performed demonstrating the efficacy and accuracy of the algorithm in detecting important physical phenomena, including bifurcation and disclination behaviors. The final numerical experiment expands the algorithm to model cholesteric liquid crystals and illustrates the full discovery power of the deflation process.

연구 동기 및 목표

• 표준 최소화 기법이 다수의 국소 극값으로 인해 놓치기 쉬운 액정 시스템에서 다수의 서로 다른 평형 상태를 탐지하기 위한 강력한 수치적 방법을 개발하는 것.
• 원래는 PDE 및 보조 문제에 사용된 탈락 기법을 액정 자유에너지 최소화의 제약 조건이 있는 변분 프레임워크에 적응시키는 것 — 특히 도메인 장의 국소 단위 길이 제약 조건 ||n|| = 1을 포함한다.
• 탈락 기법을 내재된 반복법과 기존의 다중 격자 해법기법과 통합하여 대규모 시뮬레이션에서 계산 효율성과 수렴 안정성을 향상시키는 것.
• 이 방법이 난류 및 콜레스테릭 상에서 분기, 불연속성 결함, 나선형 구조와 같은 복잡한 물리적 현상을 어떻게 해결할 수 있는지 보여주는 것.
• 동적 해 탐색을 통해 높은 에너지 장벽을 가진 해를 포함한 에너지적으로 유리한 구성 상태를 정확하고 체계적으로 탐색할 수 있도록 하는 것.

제안 방법

• 탈락 기법은 비선형 시스템에 다극자 유사 보정 항을 도입하여 이전에 계산된 해를 억제함으로써 순차적으로 서로 다른 평형 상태를 탐지하도록 수정한다.
• 이 방법은 도메인 장의 국소 단위 길이 제약 조건 ||n|| = 1을 만족하는 프랭크-후스텐 자유에너지 모델의 유한요소 변분 수식에 내장된다.
• 뉴턴 반복 과정에서 발생하는 탈락된 선형 시스템은 원래의 탈락되지 않은 시스템에 대해 설계된 다중 격자 방법을 그대로 사용하여 해석되며, 희박성과 수렴 속도를 유지한다.
• 내재된 반복법을 탈락 기법과 융합하여 각 격자 수준에서의 초기 추정치를 동적으로 정밀화함으로써 수렴 성능을 향상시키고 계산 비용을 절감한다.
• 다중 격자 해법기법 내부에서 브라에스-사라진의 이완 기법을 사용하고, 계층적 격자 간의 해 진전을 추적함으로써 수렴 안정성과 효율성을 향상시킨다.
• 미래 작업으로는 분리된 국소 최소값 골짜기 사이의 에너지 장벽을 극복하기 위한 일반화된 터널링 전략을 고려한다.

실험 결과

연구 질문

• RQ1 단위 길이 도메인 장 제약 조건이 있는 액정 자유에너지 최소화의 제약 조건이 있는 최적화 프레임워크에 탈락 기법을 효과적으로 적용할 수 있는가?
• RQ2 탈락 기법을 내재된 반복법과 다중 격자 해법기법과 통합하여 다수의 해를 탐지하면서도 계산 효율성을 유지할 수 있는가?
• RQ3 탈락 기법이 난류 및 콜레스테릭 액정에서 불연속성 결함과 나선형 구조와 같은 복잡한 물리적 현상을 어느 정도 해결할 수 있는가?
• RQ4 초기 추정치가 진정한 최소값에서 멀리 떨어져 있을 경우, 높은 에너지 장벽을 가진 해를 탐지하는 데 이 방법은 어떻게 성능을 발휘하는가?
• RQ5 어려운 구성 상태에서 수렴 성능을 향상시키기 위해 탈락 과정을 적응형 또는 동적 초기 추정치 전략으로 개선할 수 있는가?

주요 결과

• 탈락 기법은 표준 해법기법이 오직 하나의 해로 수렴하는 상황에서도 난류 액정에서 불연속성 결함과 분기된 구조를 포함한 다수의 서로 다른 평형 상태를 성공적으로 탐지하였다.
• 콜레스테릭 시스템의 경우, y축을 중심으로 2π의 회전을 가진 왼쪽 나선형 구조를 발견하였으며, 자유 에너지는 2.984×10⁻⁸로 분석 예측과 일치하였다.
• 256×256 격자에서 알고리즘은 최종 해를 발견하기 위해 총 253회의 뉴턴 반복을 수행하였으며, 이 중 63회의 반복이 탈락 과정에 기인하여 효과적인 해의 고립이 이루어졌음을 시사한다.
• 최종 격자에서 탈락된 시스템의 평균 다중 격자 반복 횟수는 63.0이었으며, 이는 탈락에 의해 유도된 수정에도 불구하고 다중 격자 수렴이 여전히 강력함을 보여준다.
• 콜레스테릭 해의 경우 작업 단위(WU) 수는 493.0으로 첫 번째 해의 100.7보다 크게 증가하였으며, 이는 다수의 상태를 탐지하는 데 소요된 계산 비용을 반영한다.
• 알고리즘이 높은 정확도와 효율성을 보였으며, 자유 에너지 값이 분석 기준치와 10⁻⁸ 이내로 수렴하여 계산된 평형 상태의 신뢰성을 확인하였다.