論文要約：表面張力効果によるマクロ相への二成分液滴の溶解

タイトル: Dissolution of a two-component drop onto macrophase due to surface tension effect
著者: Alexey Kabalnov (Ink Splat Company)

1. 背景と問題設定

分散系におけるオストワルド熟成（Ostwald ripening）を抑制または遅延させるために、難溶性成分を添加する手法は古くから知られています。しかし、添加物が媒体に対して「完全に不溶」ではなく、「限定的な溶解度（finite solubility）」を持つ場合のメカニズムは十分に解明されていませんでした。

本論文は、この未解決の問題を扱うシリーズの第 1 報として、オストワルド熟成理論に先行する単純なモデル、すなわち「表面張力による過剰ラプラス圧力（excess Laplace pressure）の影響下での、二成分液滴の溶解過程」を解析することを目的としています。このモデルは、成分のモル分率（0 から 1 の全範囲）、モル体積の違い、および活量係数（非理想性）を考慮した一般化された理論の構築を目指しています。

2. 手法と理論的枠組み

2.1 物理モデル

• 系: 媒体中に存在する、2 種類の完全に混和する液体成分（成分 1：易溶性、成分 2：難溶性）からなる球状液滴。
• 駆動力: 液滴内部のラプラス圧力（$2\sigma/r$）による化学ポテンシャルの上昇。
• 境界条件: 液滴は無限大の半径を持つマクロ相（初期組成と平衡状態）に隣接しており、マクロ相の組成は一定と仮定される。
• 基本方程式: 定常状態における各成分の拡散流束（$J_1, J_2$）は、フィックの法則と界面での濃度差に基づいて記述される。

2.2 化学ポテンシャルの解析

液滴の溶解に伴う化学ポテンシャルの変化（$\Delta\mu$）を、以下の 2 つの要因に分解して解析した。

1. ラプラス項: 曲率に起因する圧力上昇による溶解促進。
2. ラウール項: 液滴内部の組成変化（特に難溶性成分の濃縮）による化学ポテンシャルの変化。

2.3 ロック・イン近似（Lock-In Approximation）

難溶性成分（成分 2）の濃度が増加すると、そのラウール項がラプラス項を相殺する状態に達する。この状態を**「ロック・イン状態（lock-in state）」**と呼ぶ。

• 仮定: この状態では、易溶性成分（成分 1）の化学ポテンシャル差がほぼゼロとなり（$\Delta\mu_1 \approx 0$）、成分 1 の溶解は成分 2 の「出血（bleeding）」速度によって制限される。
• ギブス・デュエムの式: 一般溶液（理想・非理想問わず）に対して、化学ポテンシャルの変化を関連付けるためにギブス・デュエムの式を適用し、成分 2 の流束を導出した。

3. 主要な発見と結果

3.1 溶解の 3 つの段階

シミュレーションと理論解析により、溶解過程は以下の 3 つの段階に分類されることを明らかにした。

1. プレ・ロック・イン（Pre-lock-in）:
   • 初期段階。ラプラス圧力により易溶性成分が急速に溶け出し、液滴内の難溶性成分の濃度が急激に上昇する。
2. ロック・イン（Lock-in）:
   • 難溶性成分の濃度上昇によるラウール効果が、ラプラス圧力による溶解駆動力をほぼ完全に相殺する。
   • この状態では、液滴の組成はほぼ一定に保たれ、溶解速度は難溶性成分の拡散速度によって制御される定常状態に入る。
   • この段階では、液滴体積の減少が時間に対して線形（または古典的な立方則に従う）となる。
3. レイト・ロック・イン（Late lock-in）:
   • 液滴が極めて小さくなる終末期。ラプラス圧力が支配的になり、組成変化の仮定が崩れる。
   • 難溶性成分の濃度が急激に上昇し、溶解曲線が非線形（二次項の影響）を示すようになる。

3.2 無次元数による評価

モデルの妥当性を評価するために、2 つの無次元数を導入した。

• 第一ロック・イン数（$L_1$）: ラウール効果とラプラス効果の比率。$L_1 \gg 1$ の場合、ロック・イン状態が長く続き、古典的な立方則（体積の減少が時間に比例）が成立する。
• 第二ロック・イン数（$L_2$）: 難溶性成分の「出血」のしやすさを表す。

数値シミュレーションの結果、$L_1$ が大きい場合（例：界面張力が低い、または難溶性成分の初期濃度が高い）、プレ・ロック・インとレイト・ロック・インの期間は短く、全体の溶解寿命に対してロック・イン状態が支配的であることが確認された。

3.3 新たな溶解速度式（式 31）

従来の経験則（参考文献 [2]）を拡張し、全組成範囲（$x_2=0$ から $1$）および非理想溶液、異なるモル体積をカバーする新しい溶解速度式（式 31）を提案した。
$$w = \frac{4\pi D_2 C_2^* V_m}{\left( \frac{\phi_1}{\phi_1^* V_1} + \frac{\phi_2}{\phi_2^* V_2} \right)^{-1}}$$
（※原文の式 31 の構造に基づく要約）
この式は、純粋な成分 1 または成分 2 としての溶解速度の両極端を滑らかに接続し、モル体積の違いや非理想性を考慮することで、従来の式よりも精度が高いことが数値シミュレーションで検証された（特にモル体積が非対称な場合、従来式は 5 倍もの誤差を生むが、本式は 20% 以内の誤差に収まる）。

4. 意義と結論

• 理論的貢献: 難溶性添加物を持つ二成分系液滴の溶解メカニズムを、化学ポテンシャルと拡散の観点から統一的に記述するモデルを確立した。
• 実用的意義: 提案された式（式 31）は、インクジェット、エマルション、製剤など、多成分分散系の安定性設計において、添加物濃度や物性パラメータを最適化するための強力なツールとなる。
• オストワルド熟成への展開: 本論文で確立された「ロック・イン」メカニズムは、オストワルド熟成の抑制理論（次号論文で詳述予定）の基礎となる。特に $L_1 \gg 1$ の条件下では、粒子径分布が古典的な Lifshits-Slezov-Wagner 理論に従うことが示唆されている。

総じて、本論文は表面張力と溶解度の相互作用による二成分液滴の動的挙動を定量的に予測する枠組みを提供し、分散系の安定化メカニズムの理解を深める重要な一歩である。