Finite element analysis of a nematic liquid crystal Landau-de Gennes model with quartic elastic terms

技術要約：4次弾性項を含むネマチック液晶ランダウ・ド・ジェンヌ・モデルの有限要素解析

問題提起
本論文は、4次の弾性エネルギー項を組み込んだ、ネマチック液晶のための一般化されたランダウ・ド・ジェンヌ・モデルの数値解析に取り組んでいる。標準的なランダウ・ド・ジェンヌ・モデルは2次の弾性エネルギー密度を用いるが、近年の研究において、Golovatyら[15]は、弾性定数（$K_1, K_2, K_3$）が大きく異なる領域において、相転移をより正確に捉え、Oseen-Frankダイレクター場モデルと厳密に接続するために、高次の一般化を提案した。著者らは、これらの4次項を含む特定のエネルギー汎関数 $F(Q)$ を調査している。具体的には、Golovatyらが提案したモデルに対し、$|Q|^2|\nabla Q|^2$ に比例する項（$L_5$ で表される）を付加したものである。対処すべき主な課題は、このようなモデルに対する証明可能な安定性を備えた数値スキームが既存しないこと、および、メッシュサイズがゼロに収束する際の離散最小化子の厳密な収束証明である。

手法
著者らは、エネルギー汎関数 $F(Q)$ の勾配流に対する完全離散有限要素スキームを開発している。

• モデルの定式化: エネルギー $F(Q)$ は、テンソル係数 $S_1 = \frac{s_0}{3}I + Q$ および $S_2 = \frac{2s_0}{3}I - Q$ を伴う $\text{div } Q$ と $\text{curl } Q$ に関する5つの弾性項、およびバルクエネルギー $W(Q)$ と追加の4次項 $L_5|Q|^2|\nabla Q|^2$ で構成される。
• 数値スキーム: 半黙示的な時間ステップ手法を用いて、非線形黙示的な有限要素スキーム（スキーム3.1）を構築している。このスキームは、離散的なエネルギー散逸則を満たすように設計されており、メッシュサイズ $h$ や時間ステップ $\Delta t$ に依存しない無条件の安定性を保証する。
• 存在と一意性: 各時間ステップにおける非線形系の可解性は、Leray-Schauderの不動点定理を用いて確立される。一意性は、時間ステップサイズに関する制約（$\Delta t \leq Ch^{d+2}$）の下で証明される。
• 収束解析: 著者らは、離散エネルギー $F_h$ が $h \to 0$ のときに連続エネルギー $F$ へ $\Gamma$-収束することを厳密に証明している。この証明の極めて重要な要素は、回復列（recovery sequence）の構築であり、これは $L_5 > 0$ 項によって提供されるエネルギーの強凸性（coercivity）に依存している。著者らは、$Q$ テンソル場の潜在的な非有界性や滑らかさの欠如に対処するために、切断演算子（truncation operators）とモルフィケーション（mollification）を利用している。

主要な貢献

1. エネルギー安定スキーム: 本論文は、4次の弾性項を持つランダウ・ド・ジェンヌ・モデルの勾配流に対して、証明可能であり、かつエネルギー的に安定で収束する初の数値スキームを提示している。
2. $\Gamma$-収束の証明: 著者らは、離散エネルギーの全域的最小化子が、連続エネルギーの全域的最小化子に収束するという厳密な証明を提供している。彼らは、$L_5$ 項の包含が、この証明に必要な回復列を構築するために数学的に必要であることを示しているが、数値実験の結果は、$L_5$ が小さい場合やゼロの場合でもダイナミクスは同様であることを示唆している。
3. 理論の拡張: 本研究は、Golovatyら[15]によって確立された理論的枠組みを拡張し、相関長さが消失する極限におけるOseen-Frankモデルとの接続および一軸状態の回復に関する彼らの結果が、単純な適応によってこの修正モデルにおいても保持されることを示している。

数値結果
著者らは、スキームを検証し、4次モデルの物理的含意を探求するために、2次元における数値実験を行っている。

• 収束テスト: スキームは、空間要素幅および時間ステップサイズの双方に対して2次の収束性を示す。
• 相転移: 等方相からネマチック相への相転移のシミュレーションは、物理的な観察[14]と一致する、非球面対称な境界を持つ欠陥やタクトイド（ネマチック液滴）の形成を正常に再現している。
• 弾性定数の格差: 弾性定数（特に $L_1$）を変化させることで、シミュレーションは異なる欠陥ダイナミクスを捉える。例えば、次数1および次数-1のタクトイドのシミュレーションでは、4次モデルは標準的な2次のランダウ・ド・ジェンヌ・モデルとは著しく異なる欠陥構造（V字型のテールや特定の渦の配置など）を生み出す。
• 2次モデルとの比較: 4次モデルは、標準的な2次の1定数近似よりも豊かなダイナミクスを捉えられることが示されている。具体的には、4次モデルは、物理実験で観察されているものの標準的な2次モデルでは捉えられない特徴である、キンク（折れ曲がり）を持つ境界を持つタクトイドの形成を可能にする。

意義と主張
本論文は、提案された数値手法が堅牢であり、標準的な2次モデルでは到達不可能な複雑な液晶挙動をシミュレートできると主張している。著者らは、$L_5$ 項が$\Gamma$-収束の数学的証明には不可欠である一方で、数値結果は $L_5$ が小さい、あるいはゼロの場合でもモデルは同様に振る舞うことを示しており、これはOseen-Frankエネルギーから導出された4次項の物理的な妥当性を示唆していると強調している。本研究は、高度に格差のある弾性定数や複雑な欠陥ダイナミクスを持つ材料を研究するためのツールとして、高次のランダウ・ド・ジェンヌ・モデルの厳密な数値解析の基礎を提供するものである。また、著者らは、$\Gamma$-収束のために開発された手法が、4次のエネルギー項を持つ他の関連モデルにも適用可能である可能性についても言及している。