Solid adsorption: the missing mechanism for surfactant contact lines -- a phase-field approach

技術要約：固相吸着：界面における界面活性剤接触線の欠落したメカニズム – フェーズフィールド・アプローチによる検討

問題提起
界面活性剤は、界面張力や濡れ性を変化させることで、多相系において極めて重要な役割を果たす。多様な系における堅牢な実験的事実として、界面活性剤は平衡接触角を系統的に減少させ、親水性および疎水性の両方の表面を、より高い濡れ性へとシフトさせる。しかし、可溶性界面活性剤を含む移動接触線の問題に関する既存の数値計算フレームワークには、根本的な限界が存在する：それらは固相壁への界面活性剤の吸着を無視している。流体－流体間の界面張力の減少のみを考慮する現在のモデルは、界面活性剤が単にヤングの方程式を再スケール化するだけであると予測する。これは、実験結果との定性的な不一致を招く。すなわち、これらのモデルは親水性基質をより親水的に、疎水性基質をより疎水的にすると予測してしまうのである。さらに、これらのモデルは「オートフォビング（自発的撥水現象）」（液滴の自発的な収縮）のような現象を説明できない。この現象には、固液界面と固気界面の非対称な修飾が必要となる。著者らは、固相吸着の欠如こそが、この不一致を解消するために必要な「欠落したメカニズム」であると特定している。

手法
著者らは、界面および固相表面への吸着の両方を組み込んだ、可溶性界面活性剤を含む二相流のための熱力学的に整合したフェーズフィールドモデルを開発した。このフレームワークは、熱力学第二法則に準拠した変分原理から導出されている。

• 支配方程式: モデルは、非圧縮ナビエ・ストック方程式を、相場変数（$\phi$）および界面活性剤濃度（$\psi$）に関するカーン・ヒリアード方程式と結合させている。バルクの自由エネルギーには、混合エネルギー、界面活性剤の混合エントロピー、界面吸着エネルギー、および界面活性剤の溶解性に関するバルク・ペナルティが含まれる。
• 固相吸着メカニズム: 本研究の主要な革新は、壁の自由エネルギーの定式化である。著者らは、局所的な界面活性剤濃度として固液表面張力の減少をモデル化するために、ラングミュア・シスズコウスキーの状態で方程式を導入している。これにより、固液界面と固気界面に対して非対称な吸着強度（$\beta_{s1}$ および $\beta_{s2}$）を設定することが可能となる。
• 境界条件: 変分法による導出により、壁における相場および界面活性剤濃度の動力学的境界条件が得られる。決定的なことに、界面活性剤の自由エネルギーに勾配ペナルティ項（$|\nabla \psi|^2$）を含めることで、界面活性剤の式は空間に関して4階の微分方程式となる。この数学的構造は、境界における界面活性剤濃度の非自明な法線勾配を規定するために不可欠であり、これにより固相吸着の効果をバルク流体分布へと伝達する。
• 数値実装: システムは、時間積分に2次のアダムス・バッシュフォース法を用い、ナビエ・ストック方程式には分数ステップ法を用いた、スタッガード・デカルト格子上で解かれる。変則的な界面活性剤移動率から生じる変数係数ラプラシアンを扱うために、方向依存の順方向・逆方向分割戦略が採用されている。

主な結果
モデルは、クリーンな液滴の広がりに関する解析的ベンチマークに対して検証され、実験データおよび先行する数値的研究と比較された。

1. 濡れ性の不一致の解決: 固相吸着を含めることで、モデルはすべての初期接触角に対して、普遍的な親水性へのシフト（接触角の減少）を予測する。これは、界面活性剤の添加によって親水性および疎水性の両方の表面がより濡れやすくなるという実験的観察と一致する。対照的に、壁吸着を含まないシミュレーションでは、基質の既存の濡れ特性を増幅させるという誤った傾向が再現された。
2. 界面曲率とマランゴニ応力: 固相吸着の導入により、接触線付近での流体－流体界面の微細ながらも重要な曲がりが明らかになった。この曲率は結合メカニズムによって生じる：固相吸着が壁付近に界面活性剤濃度の勾配を生み出し、それがマランゴニ応力を発生させてバルク流を駆動する。非圧縮性の制約がこの流れを局所的な圧力場へと変換し、それが法線応力の跳びを修正して、界面の曲率を調整させる。この特徴は、固相吸着を欠いたモデルには存在しない。
3. オートフォビングの予測: 吸着の非対称性を反転させること（固気界面への強い吸着と、固液界面への弱い吸着）、によって、モデルはオートフォビングを正常に予測した。このレジームでは、液滴はより大きな平衡接触角へと自発的に収縮する。シミュレーションでは、文献で報告されている実験範囲と一致する、約$15^\circ$から$19^\circ$の接触角増加が示された。この現象は、接触線の前方にある「乾燥した」基質への界面活性剤の優先的な蓄積によって駆動され、固液間の表面張力に対して固気間の表面張力を劇的に減少させる。
4. スケーリング則: 著者らは、壁における界面活性剤吸着層の厚さに関するスケーリング則を導出した。彼らは、層の厚さが界面活性剤の拡散パラメータに依存するレジーム（低勾配レジーム）と、壁における界面活性剤の強さに依存するレジーム（高勾配レジーム）の2つのレジームを特定した。

意義および主張
本論文は、固相表面への吸着の導入が、界面活性剤を含む接触線ダイナミクスの予測的モデリングに必要な「欠落したメカニズム」を提供すると主張している。著者らは、このメカニズムなしには、数値モデルは濡れ性が向上するという普遍的な実験的傾向を定性的に再現することも、オートフォビングのような複雑な現象を説明することもできないと断言している。

本研究は以下のことを確立している：

• 固相吸着は単なる二次的な効果ではなく、接触線における力の基本バランスを変化させる中心的な駆動力であること。
• 提案されたモデルの熱力学的整合性、特に勾配ペナルティを含む自由エネルギー関数からの境界条件の導出は、界面活性剤分布、マランゴニ応力、および界面曲率の結合を捉えるために不可欠であること。
• 本モデルは、先行する数値予測と実験的事実との間の溝を埋めることに成功しており、壁吸着の非対称性のみに基づいて、広がり（spreading）の促進と自発的な脱濡れ（dewetting）の両方を説明できる統一的なフレームワークを提供していること。

著者らは、本研究が先行する定式化と実験との間の定性的な不一致を解決したことを結論付けており、固相吸着が、界面活性剤による濡れ現象を正確にモデル化するために不可欠な決定要因であることを示している。