Analysing gelation transition through fractional viscoelasticity and Mittag-Leffler-Prabhakar function

技術要約：分数階粘弾性とミッタク・レフラー・プラバカール関数を通じたゲル化転移の解析

問題提起
粘弾性液体から粘弾性固体（またはその逆）への変容プロセスであるゲル化転移は、スケールフリーなべき乗則粘弾性を特徴とする臨界ゲル状態によって記述される。分数階微積分、特にスプリングポット（spring-pot）要素は、このべき乗則挙動を記述するための自然な数学的言語を提供するが、既存の文献には重大な限界が存在する。従来の分数階粘弾性モデルは、主に臨界ゲル点における状態のみに適用されているか、あるいはゲル前（ゾル）とゲル後（ゲル）の各領域に対して独立して適用される傾向があった。決定的なことに、既存の枠組みでは、系が臨界点を通過する際の、架橋度（$p$）に対する動的弾性率およびその一次導関数の連続性を明示的に強制するものがない。この欠落により、対称性の条件（$\kappa_S = \kappa_G = \kappa$）およびハイパースケーリング関係（$n = z/(z+s)$）を、経験的な観察結果としてではなく、理論的な必然性として厳密に導出することが妨げられている。さらに、標準的な分数階機械組立体は、許容される全臨界指数の範囲において、物理的制約（緩和の単調性など）を満たせない場合が多い。

手法
著者らは、ゲル化転移における線形粘弾性特性の進化をモデル化するために、分数階粘好性と3パラメータのミッタク・レフラー・プラバカール（Mittag-Leffler-Prabhakar）関数に基づいた統一的な理論的枠組みを開発した。その手法は以下の通りである：

1. 制約の定式化： モデルは、固有の熱力学的制約（因果律、非負の弾性率、単調性）および臨界ゲル化理論に見られる特定のスケーリング関係を満たすように構築されている。最も重要な点は、本フレームワークが、臨界ゲル点（$p_c$）における複素弾性率 $G^*(\omega)$ およびその $p$ に関する一次導関数の連続性を強制している点である。
2. ゲル前（Pre-Gel）モデリング：
   • モデル1： 2つのスプリングポットと1つのダッシュポットの直列結合。これは2パラメータのミッタク・レフラー関数による緩和弾性率を与える。しかし、このモデルは $n + \kappa < 1$ であるシステムに限定され、臨界点近傍（$\varepsilon \to 0$）においてのみ漸近的に正確である。
   • モデル2： 3パラメータのミッタク・レフラー・プラバカール関数を用いた一般化。このモデルは $n + \kappa < 1$ の制限を取り除き、すべての $n, \kappa \in (0,1)$ に対して物理的に許容可能な記述を可能にする。また、プラバカール指数 $\gamma_S$ が粘度発散指数 $s$ に等しいことを一意に決定する。
3. ゲル後（Post-Gel）モデリング：
   • モデル1： フックの法則に従うスプリング（平衡弾性率を提供）と、スプリングポットの並列配列。
   • モデル2： プラバカール関数の打ち切り級数展開。これは、特定の係数置換の下でモデル1と代数的に等価であることが示されているが、より制約された、物理的に構造化されたパラメータセットを提供する。
4. 連続性の強制： ゲル前モデルとゲル後モデルを同時に解く。動的弾性率の導関数が $p_c$ で一致するという要求により、両側の緩和スケーリング指数が同一（$\kappa_S = \kappa_G = \kappa$）であり、かつ主要項の偏差係数が等しいことが強制される。
5. 実験的検証： 提案されたモデルは、以下の2つの異なる実験データセットを用いて検証された：
   • ゾル－ゲル転移（$n=0.45$）を起こす化学架橋型ポリジメチルシロキサン（PDMS）の時間領域緩和弾性率データ。
   • ゲル－ゾル転移（$n=0.77$）を起こす熱可逆性ポリビニルアルコール（PVA）溶液の周波数領域動的弾性率データ。

主な貢献と結果

• 統一された枠組み： 本論文は、連続性の制約を通じてゲル前およびゲル後の領域を厳密に結びつけ、転移を分離した状態ではなく、単一の連続的なプロセスとして扱う、初の分数階粘弾性フレームワークを提示している。
• ハイパースケーリングの理論的必然性： 動的弾性率の一次導関数の連続性を強制することにより、著者らは、対称性の条件（$\kappa_S = \kappa_G$）およびハイパースケーリング関係（$n = z/(z+s)$）が、経験的な偶然の一致ではなく、物理的一貫性を保つために必要な理論的必然であることを証明した。
• プラバカール関数の優位性： プラバカールに基づく「ゲル前モデル2」は、標準的なスプリングポット／ダッシュポット組立体よりも優れていることが示された。このモデルは、漸近極限で消失するダッシュポット要素を必要とすることなく、パラメータ $\gamma_S = s$ を通じてゼロせん断粘度の発散（$\eta_0 \sim \varepsilon^{-s}$）を本質的に捉えている。これは、プラバカール関数が機械的なアナログよりも、ゲル前の緩和ダイナミクスを表現する上でより自然な表現であることを示唆している。
• モデルに依存しないフィンガープリント： 解析により、臨界ゲル状態における、モデルに依存しない周波数独立なレオロジー的フィンガープリントを特定した：すなわち、制御パラメータに関する損失弾性と貯蔵弾性率の導関数の比、$(\partial G'' / \partial G')_{\varepsilon \to 0} = \tan((n-\kappa)\pi/2)$ である。この量は、臨界指数 $n$ と $\kappa$ によって一意に決定される。
• 実験との一致： 提案されたモデルは、化学架橋系（不可逆）および物理的会合系（可逆）の両方の実験データを再現することに成功した。適合度は臨界点近傍で極めて高く、系がゾルまたはゲル領域へと遠ざかるにつれて偏差が増大するが、これは導出された式の漸近的な性質と一致している。

意義
本論文の主要な意義は、分数階粘弾性モデリングとゲル化の基礎的なスケーリング則との間の断絶を解消した点にある。導関数の連続性を厳格な物理的制約として統合することで、本研究はハイパースケーリング関係を、経験的な観察から理論的な要請へと昇華させた。プラバカール指数 $\gamma$ を、クラスターサイズ分布や緩和スペクトルに関連する物理量（粘度指数 $s$）と等しい、測定可能な物理的パラメータとして特定したことは、数学的なフィッティングパラメータと根底にある物理学との間の新たな繋がりを提供している。最後に、モデルに依存しない周波数独立なフィンガープリント $(\partial G'' / \partial G')$ の特定は、ゲル化の化学的性質や経路に依存しない、臨界ゲル状態を特徴付けるための堅牢かつ普遍的な基準を提供するものである。