Beyond binary scission: a generalized three-species cascade breakage model for wormlike micellar solutions

技術要約：虫状ミセル溶液のための一般化された3成分カスケード断裂モデル

問題提起
虫状ミセル（WLM）溶液は、自己組織化されたミセルネットワークの連続的な断裂と再結合によって駆動される、粘弾性、せん断希薄化、およびせん断帯形成（シェアバンディング）を含む複雑なレオロジー挙動を示す。既存の理論的枠組み、特に広く用いられている2成分モデル（Vasquez-Cook-McKinley（VCM）モデルなど）は、ミセル集団を粗い二値表現（典型的には、高度に絡み合った長い鎖と短い断片のみを区別する）に依存している。これらのモデルは特定のマクロな特徴を捉えることには成功しているものの、根本的な限界を抱えている。すなわち、多スケールなカスケード断裂と、それに伴う実際のミセルネットワークに見られる連続的な緩和レジームを明示的に解像できていない。この二値近似は中間的な構造状態をフィルタリングしてしまい、広範な緩和スペクトル（特に実験で観察される動的弾性率の高周波領域での上昇）を正確に記述することを不可能にし、過渡応答やせん断帯形成メカニズムの物理的な解釈を制限している。

手法
これらの限界に対処するため、著者らは中間的なミセル状態を動力学的枠組みに導入した3成分カスケード断裂モデルを開発した。本モデルは、ミセルネットワークが以下の3つの離散的な構造状態で構成されていると仮定している：

1. 種A： ゲルネットワーク（主要な、高度に絡み合った状態）。
2. 種B： 長い鎖（中間状態）。
3. 種C： 短い鎖（絡み合いが解けた断片）。

数学的定式化は、連続的な位相空間分布を用いた動力学理論に基づいている。著者らは、対流による伸長、拡散、および反応項を含む、3つの種の確率密度関数に対する結合したシュモルホフスキー方程式を導出している。主な手法の特徴は以下の通りである：

• 階層的動力学： モデルは、ゲルネットワーク（A）が長い鎖（B）へと断裂し、その後、長い鎖（B）が短い鎖（C）へと断裂するというカスケード経路を定義している。質量保存を保証するために、再結合経路も定義されている。
• 流動依存的な切断： ゲルネットワークおよび長い鎖の断裂速度は、部分伸長鎖対流（PEC）メカニズムを介して局所的なマクロ流に結合している。これにより、変形が切断を加速させる応力依存的な断裂メカニズムが生成される。
• 構成方程式： モデルは、数密度とコンフォメーションテンソルの無次元進化方程式を導出し、これらをマクロなコーシー運動量方程式と結合させている。全応力は、全3種の弾性応力の重ね合わせである。
• 数値実装： 著者らは、均一流条件下（線形粘弾性と定常せん断応答の研究のため）および不均一な円筒型クーエット流（せん断帯形成、過渡的な立ち上がり、および緩和の研究のため）の両方でモデルを解析している。

主な貢献と結果

1. 3モード粘弾性応答：
   均一流の下で、モデルは解析的に離散的な3モード・マックスウェル構造を導き出す。2成分モデルとは異なり、この枠組みは動的貯蔵弾性率（$G'$）および損失弾性率（$G''$）の高周波領域での上昇を捉えることに成功している。著者らは、中間種（B）と短い鎖（C）によって提供される追加の緩和モードが、CPyCl/NaSal溶液の実験的なマスターカーブへの適合性を、特に高周波領域において大幅に改善することを示している。

2. 非単調な流動曲線とせん断帯形成：
   モデルは、構成的不安定性の徴候である、厳密に負の勾配を持つ領域を含む非単調な構成曲線を予測する。不均一なクーエット流シミュレーションにおいて、この不安定性は、共存するせん断帯（高せん断速度領域と低せん断速度領域）の自発的な形成を引き起こす。3成分モデルは、中間種（B）がバンド界面において明確な遷移層を形成することを示しており、これはゲルネットワークが支配的な低せん断バンドと、短い鎖が支配的な高せん断バンドを繋ぐものである。これは、二値モデルよりも解像度の高いバンド界面の微視的構造記述を提供する。

3. 過渡ダイナミクスとヒステリシス：
   モデルは複雑な過渡挙動を捉えている：

• 応力オーバーシュート： 立ち上がり時において、弾性的な負荷と加速されたカスケード断裂の競合により、顕著な応力オーバーシュートが生じる。
• 流動ヒステリシス： モデルはヒステリシスを示し、最終的な定常状態は流動履歴（静止状態からの立ち上がり vs 高せん断からのステップダウン）に依存する。
• 多段階緩和： せん断停止時、モデルは多段階緩和プロセスを予測する。単一の指数関数的減衰ではなく、応力は逐次的な段階を経て緩和される：短い鎖（C）による急速な解放、長い鎖（B）による中間的な調整、そしてゲルネットワーク（A）の緩やかな回復。これにより、実験でしばしば観察されるが簡略化されたモデルでは捉えにくい「長期的な回復タイムスケール」を解明している。

4. 連続体機械的整合性：
   モデルは、大きな振幅のステップひずみにおいて、ロッジ・マイヤーソン関係（$N_1/\sigma_{r\theta} = \gamma_0$）を満たしており、有限変形下でのテンソル的一貫性と物理的妥当性を確認している。

意義と主張
本論文は、提案された3成分フレームワークが、離散的な現象論的ネットワークモデルと連続的な動力学理論（Catesの元のレプテーション・反応理論など）との間の溝を埋めるものであると主張している。中間的な構造状態を制御するカスケード動力学を明示的に解明することで、本モデルは、虫状ミセルのせん断帯形成に関する物理的に解釈可能な枠組みを提供する。

著者らは、このアプローチが以下の間の本質的なつながりを捉えていると断言している：

• カスケード微視的構造の進化（A $\to$ B $\to$ C）。
• 広範な緩和ダイナミクス（3つの異なるタイムスケール）。
• マクロな流動局在化（せん断帯形成）。

その意義は、マクロな計算の扱いやすさを維持しながら、欠落していた中間的な緩和モードを回収できる点にある。これは、複雑なミセル流の豊かな過渡的および定常的なダイナミクスを解釈するための、より包括的な物理的基礎を提供し、将来のマイクロ流体実験や産業応用における解釈を助ける可能性がある。本論文は新しい実験プロトコルを提案するものではなく、既存のレオロジー現象をより良く解釈するための理論的ツールを提供している。