Self-consistent field theory of semiflexible nematics: Density-nematic coupling, anisotropic elasticity, and defect core sizes

技術要約：半柔軟性ネマティックにおける自己整合場理論

問題提起
小分子熱的液晶の線形応答はよく理解されている一方で、高分子、ミセル集合体、生物学的フィラメントからなるリオトロピック（液晶）液晶の挙動は、理論的および実験的に大きな課題を提示している。中心的な問題は、リオトロピック系で観察される大きな弾性異方性の起源と、鎖の密度とネマティック秩序との間の結合である。Onsagerモデルに大きく依拠する既存の微視的理論は、空間的に変化する半柔軟性ネマティックにおける弾性、トポロジカル欠陥、および密度・秩序結合を同時に扱うことに苦慮している。さらに、マクロな連続体理論（Oseen–FrankおよびLandau–de Gennes）は、通常、基礎となる分子統計との明示的な関連を持たない現象論的な係数を導入しており、弾性異方性および欠陥コア構造の微視的な起源は未探索のまま残されている。

手法
著者らは、半柔軟性ネマティックをワーム状鎖としてモデル化する自己整合場理論（SCFT）の枠組みを採用している。系は以下の要素を含むハミルトニアンによって定義される：

1. 曲げエネルギー： 永続長（$l_p$）および経路長（$L$）によって支配される。
2. Maier–Saupe (MS) 相互作用： ネマティック秩序を駆動する異方的な配向相互作用（$u_2$）。
3. 等方的な排除体積： リオトロピック条件を表す溶媒媒介相互作用（$u_0$）（溶媒効果をFlory–Hugginsパラメータを介して暗黙的に扱う）。

研究は3つの段階で進行する：

• 一様状態解析： 著者らは、鎖セグメント密度（$\rho$）および一軸秩序パラメータ（$S$）に依存する有効自由エネルギー汎関数を導出する。これにより、等方相－ネマティック相共存（バイノーダル）および安定限界（スピンノーダル）の境界の決定が可能となる。
• 線形応答解析： 密度、一軸秩序、およびディレクター配向の摂動に対する系の応答を調べるために、弱い外部場を導入する。これにより、逆感受率行列が得られ、そこからFrank弾性定数（$K_1, K_2, K_3$）および相関長が抽出される。
• 欠陥構造シミュレーション： 完全なSCFT方程式を$+1/2$ディスクリネーションに対して数値的に解き、ディレクターのプロファイルとコアサイズを解析し、剛直鎖の極限と半柔軟性レジームを比較する。

主要な貢献と結果

• 密度－ネマティック結合と安定性：
  導出された有効自由エネルギーは、密度（$\rho$）と秩序度（$S$）の間の強い結合を明らかにする。著者らは、強い密度－ネマティック結合がネマティック状態を不安定化させることを示している。スピンノーダル不安定性の位置は、排除体積とネマティック相互作用の強さの比（$u_0/u_2$）に決定的に依存する。対照的に、ディレクターの歪みは高次においてのみ密度や秩序の変化と結合するため、線形応答レジームでは実質的にデカップル（非結合）したままである。

• 弾性異方性と柔軟性：
  本研究は、鎖の柔軟性（$L/l_p$）がFrank弾性定数をどのように変化させるかを定量化している。

  • 剛直鎖： $K_3$（曲げ）は配向強度とともに急速に増加するが、$K_1$（スプレイル）および$K_2$（ツイスト）は同程度のままである。
  • 半柔軟性鎖： 柔軟性の増加は、ツイスト（$K_2$）および曲げ（$K_3$）の弾性を強く抑制するが、スプレイル弾性（$K_1$）は比較的弱く影響を受ける。
  • クロスオーバー： この微分的な抑制により、柔軟性が増すにつれて、曲げ支配からスプレイル支配の弾性へのクロスオーバーが生じる。これは、$K_1/K_3$の比の反転（$<1$から$>1$へ）によって裏付けられる。

• 相関長と異方性：
  線形応答解析は、逆感受率の対角成分の勾配係数が、ディレクターに平行な相関長が短い（$K_{\parallel} < K_{\perp}$）ことを示唆している一方で、密度と秩序の間の結合が有効係数を繰り込み（renormalize）を行うことを明らかにしている。これにより、ディレクターに平行な相関長が実際にはより大きい（$\xi_{\parallel} > \xi_{\perp}$）という異方性が生じる。この知見は、等方相－ネマティック界面および欠陥コア構造に関する実験的観察と一致している。

• 欠陥コアサイズ：
  $+1/2$ディスクリネーションの数値シミュレーションは、コアサイズが異方的であり、密度勾配がディレクターに平行な方向において半径が大きくなることを示している。本研究は、特徴的な欠陥コアサイズが平衡相関長（$\xi_Q$）に比例することを確立しており、これはLandau–de Gennesのスケーリングと一致している。著者らは、リオトロピック・クロモニクスがしばしば異常に大きなコアを示す（おそらく集合体の切断による）ことに触れつつ、多分散性を除外している現在のモデルが、密度ゆらぎに起因するコアサイズと相関長の間のスケーリング関係を捉えていると述べている。

• スピンノーダル挙動：
  系がネマティック・スピンノーダル（$u_0$の変動によって制御される）に近づくと、相関長は $\xi \sim (u_0 - u_{0,i})^{-1/2}$ のスケーリングに従って発散する。しかし、欠陥構成においては、コアにおける局所的な等方性ドメイン（タクトイド）の形成が、スピンノーダルに達する前にネマティック秩序を抑制することで、この発散を阻止する。

意義
本論文は、半柔軟性ネマティックにおける分子の柔軟性、密度－ネマティック秩序結合、および弾性異方性を結びつける統一された微視的枠組みを提供する。Frank弾性定数および相関長をSCFT形式内の鎖の統計から直接導出することで、本研究は、リオトロピック系で見られる大きな弾性異方性の微視的な起源を解明している。得られた結果は、欠陥構造を記述するための連続体理論（Landau–de Gennesなど）の使用を正当化すると同時に、弾性応答と安定限界を正確に予測するためには、密度ゆらぎと鎖の柔軟性を組み込む必要があることを強調している。著者らは、本フレームワークが自己組織化された多分散鎖や可逆的凝集への将来的な拡張の基礎を提供すると結論づけているが、現在の理論的分析の範囲を超えて、具体的な新しい実験プロトコルや応用を提案してはいない。