[論文レビュー] A Deflation Technique for Detecting Multiple Liquid Crystal Equilibrium States

James H. Adler, D. B. Emerson+2 more|arXiv (Cornell University)|Jan 27, 2016

Liquid Crystal Research Advancements参考文献 29被引用数 3

ひとこと要約

本稿では、ネマチックおよびキロステリックな液体結晶における複数の異なる平衡状態を、フランク＝フジノ自由エネルギー模型を用いて系ativelyに検出するための、ネストドイテレーションおよびマルチグリッド・ソルバと統合されたデフラーション技法を提案する。この手法により、既に発見済みの解を除外するように非線形系を変更しつつ、有限要素法のスパarsityを維持し、既存のマルチグリッド予め条件付きの再利用を可能にすることで、局所的および大域的最小化子、特に複雑な欠陥構造やヘリカル相を効率的に特定できる。

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.

研究の動機と目的

• 標準的な最小化手法が複数の局所的極値を有するため見過ごされがちな、液体結晶系における複数の異なる平衡状態を系統的に検出するための堅牢な数値的手法の開発。
• 元来偏微分方程式および補完性問題に用いられていたデフラーション技法を、配向子場の点毎に単位長さ制約が課された液体結晶自由エネルギー最小化の制約付き変分枠組みに適応すること。
• デフラーションをネストドイテレーションおよび既存のマルチグリッド・ソルバと統合し、大規模シミュレーションにおける計算効率および収束のロバスト性を向上させること。
• ネマチックおよびキロステリック相において、分岐、ディスクリネーション、およびキラルヘリカル構造を含む複雑な物理現象の解明における本手法の有効性を示すこと。
• 動的解探索を用いて、高エネルギー障壁を有する状態を含む、エネルギー的に有利な配置の正確で体系的な発見を可能とすること。

提案手法

• デフラーション技法は、非線形系に極に類似した補正項を導入することで、既に計算済みの解を抑制し、連続的に異なる平衡状態を発見可能とする。
• この手法は、点毎に ||n|| = 1 を満たす制約を課したフランク＝フジノ自由エネルギー模型の有限要素変分式に埋め込まれる。
• ニュートン反復におけるデフラートド線形系は、元の未デフラート系に特化したマルチグリッド法によって解かれる。これにより、スパarsityと収束速度が維持される。
• ネストドイテレーションをデフラーションと組み合わせ、各メッシュレベルにおける初期推定値を動的に精錬することで収束性を向上させ、計算コストを低減する。
• マルチグリッド・ソルバ内ではブラース＝サラジンの緩和法を用い、グリッド階層全体にわたる解の進行状況を追跡することで、ロバスト性と効率性を向上させる。
• 今後の課題として、分離された最小値の谷間を越えるエネルギー障壁を克服するための一般化されたトンネル戦略の検討がなされる。

実験結果

リサーチクエスチョン

• RQ1 単位長さの配向子制約を有する液体結晶自由エネルギー最小化の制約付き最適化枠組みに、デフラーションが効果的に適応可能か？
• RQ2 デフラーションをネストドイテレーションおよびマルチグリッド・ソルバとどのように統合すれば、複数の解を発見する際の計算効率を維持できるか？
• RQ3 デフラーション法は、ネマチックおよびキロステリックな液体結晶におけるディスクリネーションやヘリカル構造といった複雑な物理現象をどの程度解明できるか？
• RQ4 初期推定値が真の最小値から大きく離れている場合でも、本手法は高エネルギー障壁を有する解を効果的に特定できるか？
• RQ5 挑戦的な配置において収束性を向上させるために、適応的または動的初期推定値戦略をデフラーションプロセスに統合できるか？

主な発見

• デフラーション法は、ディスクリネーションや分岐構造を有する状態を含め、ネマチック液体結晶において複数の異なる平衡状態を効果的に発見した。標準的なソルバが唯一の解に収束する状況でも同様に有効であった。
• キロステリック系では、y軸周りに2π回転する左旋のヘリカル構造が発見され、自由エネルギーは2.984×10⁻⁸に達し、解析的予測と一致した。
• 256×256グリッド上では、最終解の発見に253回のニュートン反復が必要であったが、そのうち63回がデフラーションプロセスに起因しており、解の分離が効果的に行われたことが示された。
• デフラートド系の最細グリッドにおける平均マルチグリッド反復回数は63.0であった。これは、デフラーションに起因する修正にもかかわらず、マルチグリッド収束が依然としてロバストであることを示している。
• キロステリック解のためのワーキングユニット（WU）数は493.0であった。これは最初の解の100.7に比べて顕著に高く、複数状態の発見に伴う計算コストの増加を反映している。
• 本アルゴリズムは高い精度と効率を示し、自由エネルギー値が解析的ベンチマークと10⁻⁸の内側に収束した。これにより、計算された平衡状態の信頼性が確認された。