Densely packed particle raft at vertically vibrated air-water interface

技術要約：垂直振動する気液界面における高密度粒子ラフト

問題提起
流体界面に形成される粒子ラフトは、固体的弾性から液体的流動に至るまで、豊かな機械的挙動を示す。準静的な変形（引張歪み、せん断、インデンテーション）に対するラフトの応答については十分に文書化されているが、強固で均質な動的励起下でのダイナミクスについては、未だ理解が進んでいない。特に、垂直振動環境における毛管相互作用、流体力学的力、および粒子の充填の相互作用は複雑である。従来の研究は、低い充填率や局所的な波の生成に焦点を当ててきたが、高密度の粒子的ラフトがグローバルな垂直振動に対してどのように応答するか、特に集団的なパターンの出現、ガラス的ダイナミクス、およびキャビテーションに関する理解には空白が残されている。

手法
著者らは、脱イオン水（$h = 5$ mm）を満たした正方形の容器（$W = 10.9$ cm）を垂直に振動させる実験装置を構築した。容器は、駆動周波数 $f \in [20, 95]$ Hz、振幅 $A \in [0.02, 0.2]$ mm の正弦波振動 $Z(t) = A \sin(2\pi ft)$ を受ける。

• 粒子: 密度 $\rho_p \approx 15$ kg m$^{-3}$、平均直径 $d = 1.4 \pm 0.2$ mm および $d = 2.1 \pm 0.2$ mm の閉細胞構造を持つ発泡スチロール球を使用した。接触角は $87^\circ \pm 6^\circ$ と測定された。
• 充填: 単層の厳密な構造を形成し、充填率 $\phi$ は容器面積に対する投影粒子面積の比として定義され、$\phi \approx 0.90$ までテストされた。
• 撮像: パターン形成は、カメラを用いて 73 Hz ($2448 \times 2048$ ピクセル) でモニタリングされた。粒子軌跡とサイクル内のダイナミクスは、高速カメラを用いて 460 Hz ($2448 \times 1024$ ピクセル) および特定の軌跡分析のために 186 Hz で捕捉された。
• 解析: 粒子位置は円形ハフ変換を介して抽出された。分散関係は、時間平均画像の2次元フーリエ変換によって決定された。粒子ダイナミクスは、平均二乗変位（MSD）、変位確率分布、および動的不均一性を測定するための4点感受率（$\chi_4$）を用いて定量化された。

主要な貢献と結果
$\phi, f, A$ を系統的に変化させることにより、本研究は以下の5つの異なるダイナミクス・レジームを特定し、相図としてまとめた。

1. 規則的な定常波: 中程度の周波数と振幅において、ラフトは古典的なファラデー波パターン（正方形、ストライプ）を、駆動周波数の半分（$f/2$）で示しながら振動する。

   • 分散関係: 著者らは波の分散を測定し、高い充填率（$\phi \approx 0.90$）では古典的な分散関係が成立しないことを発見した。有効表面張力（$\gamma_e$）と面外曲げ弾性率（$B$）を考慮した修正された関係式が、データをより良く適合させた。
   • 創発的弾性: $\phi$ が増加するにつれて、$\gamma_e$ は単調に減少し、一方で $B$ は急速に増加する。これは、高密度に充填された粒子が持続的な接触力ネットワークを発達させ、界面に顕著な曲げ剛性を付与していることを示している。

2. カオス的波: 周波数が小〜中程度で振幅が増加すると、系は秩序だった定常波から時空間カオス的パターンへと遷移する。粒子の存在は定常波を歪ませ、純粋な液体システムと比較して、欠陥形成やカオスへの追加の経路を提供することがある。

3. 空隙（ボイド）形成: 低い充填率（$\phi = 0.77$）および低い振幅において、粒子被覆が不十分なため、定常波と共存する空隙が形成される。

4. ガラス的レジーム: 高周波数（$f \ge 65$ Hz）かつ低振幅において、定常波は消失する。粒子は隣接粒子によって形成されたケージの中で「熱的」な運動を示す。

   • 劣拡散: MSD解析は、低い $\phi$ での正常拡散（$\alpha \approx 1$）から、高い $\phi$ での劣拡散挙動（$\alpha \approx 0.5$）への遷移を明らかにした。これは、ガラス転移に近づく過冷却液体に特徴的な挙動である。
   • 不均一性: 系は空間的および時間的な動的不均一性を示す。4点感受率 $\chi_4$ は遅延時間 $\tau^*$ においてピークに達し、この $\tau^*$ は $\phi$ に対して指数関数的に増加する。$\chi_4$ のピーク値は非単調に変化し、これは波によるデコレーション（無相関化）と密度駆動のケージングとの間の競合を示唆している。

5. キャビティ（空洞）形成: 高周波数（$f \ge 75$ Hz）および臨界振幅（$A^*$）において、大きなキャビティがラフト内に核生成し、拡大する。

   • メカニズム: キャビティ内には自由表面波が存在する一方、周囲のラフトはガラス状態を維持している。キャビティのサイズは、高密度ラフトの閉鎖圧力と、波と粒子の相互作用（波の腹から粒子が遠ざかるドリフト）によって駆動される外向きのドリフトとのバランスによって安定化される。
   • スケーリング: キャビティ形成に必要な臨界振幅 $A^*$ は、周波数の増加とともに減少する。安定化したキャビティサイズ $L^*$ は、特徴的な注入エネルギー $E^* = A^{*2}f^2/(gd)$ に対して線形にスケールする。
   • 形状遷移: $f \approx 85$ Hz 以下では、キャビティは不規則である。これより高い周波数では、キャビティは円形かつ安定し、周囲のラフトからの抵抗が液体的なものから弾性的なものへと遷移していることを示している。

意義と主張
本論文は、垂直振動する粒子ラフトが、非平衡物理学を調査するための多目的なモデルシステムとして機能することを主張している。

• 界面力学: 本研究は、高密度粒子ラフトが、純粋な流体界面とは異なる、調整可能な有効表面張力および創発的弾性（曲げ弾性率）を示すことを実証している。
• ガラス的ダイナミクス: 本システムは、非熱的な系において、ガラス転移、劣拡散輸送、および動的不均一性を研究するための、マクロで観察可能なプラットフォームを提供する。
• 複雑な相互作用: ここでの相互作用は、理想化されたポテンシャル（例：レナード・ジョーンズ・ポテンシャル）とは異なり、複雑な多体的な毛管および流体力学的効果から生じている。振動パラメータを通じてこれらの相互作用を制御できる能力は、単純な弾塑性挙動を超えた、現実世界の高密度粒子集合体の理解への道を開くものである。

著者らは、カオス的レジームにおける過渡的な挙動や、粒子と流体の相互作用の全容としての複雑さは現在の研究の範囲外であることを認めつつ、特定されたレジームの定常状態の挙動を特徴付けることに焦点を当て、研究の範囲を控えめに設定している。