From suspensions to porous multilayers: microstructure formation and particle packing in drying colloidal films

技術要約：懸濁液から多孔質多層膜へ

問題提起
粒子懸濁液の乾燥は、電気化学デバイス（電解槽、燃料電池、電池など）や印刷エレクトロニクスに使用される機能的な多孔質層を組み立てるための基本的なプロセスである。これらのデバイスの性能は、堆積した薄膜の微細構造特性、具体的には粒子の充填状態と空隙率に決定的に依存している。しかし、蒸発誘起流、粒子間相互作用、および粒子–基板間の相互作用を含む、複雑で結合した物理プロセスにより、これらの微細構造を制御することは依然として大きな課題である。実験的研究では蒸発速度や粒子の特性の影響が探求されてきたが、材料の不透明性と測定上の制約により、微細構造の進化を直接モニタリングすることは困難である。既存の数値モデルは、粒子間の相互作用を簡略化（例：反発ポテンシャルのみを使用）したり、界面における重要な流体力学的および毛管相互作用を無視したりすることが多く、構造形成に対する予測能力に限界がある。

手法
本研究では、厚いコロイド薄膜の乾燥ダイナミクスを数値的に調査するために、格子ボルツマン法–離散要素法（LB-DEM）の結合フレームワークを採用している。

• 流体力学: 流体の流れと蒸発駆動の輸送を解明するために、格子ボルツマン法（LBM）、具体的には擬ポテンシャル多成分モデル（Shan-Chen法）を使用する。D3Q19格子を用い、拡散界面アプローチを採用することで、移動する接触線における応力特異点を自然に回避している。蒸発は、一定の密度境界条件を課すことでモデル化され、これにより拡散支配の質量移動が駆動される。
• 粒子ダイナクション: コロイド粒子は、修正されたバウンスバック境界条件を介して流体と結合した剛体球としてモデル化される。粒子の運動は、流体力（ドラッグ、リフト、および潤滑補正）と粒子間相互作用を考慮した古典的な運動方程式に従う。
• 相互作用: モデルには、ファンデルワールス力（短距離および長距離のカットオフを持つハマー定値によって近似）、重なりを防ぐためのヘルツ・ポテンシャル、および粒子–基板間のレナード＝ジョーンズ・ポテンシャルが含まれる。
• シミュレーション条件: シミュレーションは高速蒸発限界（高ペクレ数、$Pe \gg 1$）において実施され、ブラウン運動および重力沈降は無視される。本研究では、特性的なファンデルワールス引力と毛管力のスケールを比較する無次元粘着パラメータ $\chi = |F^*_{vdW}|/\gamma R$ を系統的に変化させている。

主要な貢献と結果
本研究は、毛管力に対する粒子間相互作用の強さに決定される3つの異なる充填レジームを特定している。

1. 弱い相互作用 ($\chi \approx 0.02$):

   • メカニズム: 毛管力が支配的である。液界面が下降するにつれ、粒子は吸収され、界面において六方晶の配列へと再編成される。
   • 構造: この秩序化されたトップレイヤーがテンプレートとして機能し、層状の組み立てプロセスを誘発する。その結果、軽微な欠陥を伴う、高度に秩序化された多層六方最密充填（HCP）構造が得られる。
   • モデリング: 質量保存と液体拡散に基づく単純な理論モデルが開発された。このモデルは、このレジームにおけるシミュレーションデータと良好な一致を示している。

2. 中間的な相互作用 ($\chi \approx 0.2 - 1.0$):

   • メカニズム: ファンデルワールス力が毛管圧縮と競合し始める。界面に到達する前に、サブサーフェス領域の粒子が小さなクラスターへと凝集する。
   • 構造: これらのクラスターが上層の秩序ある配列を不安定化させ、ランダム最密充填（RCP）レジームへと導く。最終的な構造は周期性を示すものの、弱相互作用の場合と比較して振幅のゆらぎが大きく、無秩序さが増大している。

3. 強い相互作用 ($\chi \approx 10.0 - 14.0$):

   • メカニズム: 強力な引力が乾燥プロセスの初期段階で急速な大規模クラスターへの凝集を引き起こす。これらのクラスターは、毛管圧力による圧縮効果に抵抗するネットワークを形成する。
   • 構造: システムは最密充填層を形成できず、代わりに高度に多孔質なネットワーク状の固体構造（凝集充填、CP）を発達させる。配位数（Coordination number）は大幅に低下し（3～5）、構造は大きな細孔と弱い結合性を特徴とする。

定量的知見

• 空隙率のスケーリング: 最終的な堆積物の空隙率 ($\psi$) は、粘着パラメータに対してべき乗則の関係に従う：$\psi \sim \chi^{2/3}$。低 $\chi$ では空隙率は低く滑らかに変化するが、$\chi$ が増加するにつれて空隙率は急速に上昇し、高粘着時において飽和限界に近づく。
• 配位数: 平均配位数は、粘着パラメータの増加に伴い、約12（HCP）から約6–8（RCP）、そして最終的に約3–5（CP）へと段階的に減少する。これは、秩序だった結晶様構造から、無秩序で結合の弱いネットワークへの転移を反映している。

意義
本論文は、粒子間相互作用が乾燥するコロイド薄膜の最終的な多孔質構造を決定する上で決定的な役割を果たすことを示している。粘着パラメータを系統的に変化させることで、著者らは六方最密充填、ランダム最密充填、および凝集充填の間の連続的な遷移をマッピングしている。これらの知見は、粒子サイズ、材料の選択、またはpHや塩濃度などの媒体特性を通じた粒子相互作用の制御による、機能層の設計のための実用的な指針を提供する。本結果は、輸送特性や電気化学的活性が空隙率と充填密度によって支配されるアプリケーションにおいて、微細構造を最適化するための予測可能な枠組みを提供するものである。