Relaxation dynamics and long-time tails explain shear-induced diffusion of soft athermal particles near jamming

この論文「Relaxation dynamics and long-time tails explain shear-induced diffusion of soft athermal particles near jamming（ジャミング近傍のソフト・アサーマル粒子のせん断誘起拡散を説明する緩和ダイナミクスと長時間テール）」の技術的サマリーを以下に示します。

1. 研究の背景と問題設定

ソフトなアサーマル粒子（泡、エマルジョン、コロイド懸濁液、粉体など）は、臨界充填率 $\phi_J$ において「ジャミング転移」を起こすことが知られており、その機械的・レオロジー的性質は臨界現象として記述されます。
既往の研究では、せん断誘起拡散係数 $D$ と系の特徴的な長さスケール（クラスターサイズや系サイズ）との関係が主に議論されてきました。しかし、熱平衡系（古典流体やガラス）では拡散が「緩和ダイナミクス」によって説明されるのに対し、アサーマル系におけるせん断誘起拡散のメカニズム、特に緩和時間と拡散係数の関係、および長時間テール（long-time tail）の役割については未解明な点が多く残されていました。

本研究の主な問いは以下の通りです：

• ジャミング近傍におけるせん断誘起拡散係数 $D$ は、Green-Kubo (GK) 公式を通じて、粒子の速度自己相関関数（VACF）の緩和ダイナミクスとどのように関連しているか？
• ジャミング状態（$\phi > \phi_J$）の準静的せん断条件下で、VACF がどのような振る舞いを示し、それが拡散係数の発散にどう寄与するか？

2. 手法

• 数値シミュレーション: 2 次元のソフト・アサーマル粒子系に対する分子動力学（MD）シミュレーションを実施。
• モデル: 直径 $d_0$ と $1.4d_0$ の 2 種類の粒子からなる系（全粒子数 $N=8192$）。粒子間力は重なり量に比例する弾性力、および速度に比例する減衰力（過減衰ダイナミクス）で記述。
• 条件: 充填率 $\phi$ をジャミング点 $\phi_J \approx 0.8433$ 付近（$0.8 \le \phi \le 0.9$）で変化させ、せん断率 $\dot{\gamma}$ を広範囲（$10^{-7} \le \dot{\gamma}t_0 \le 10^{-2}$）で変化させる。
• 解析手法:
  • 横方向の平均二乗変位（MSD）から拡散係数 $D$ を算出。
  • Green-Kubo (GK) 公式を用いて、$D$ を平均二乗速度 $\langle v_y^2 \rangle$ と速度自己相関関数 $C(\tau)$ の時間積分の積として表現。
  • 臨界スケーリング解析を行い、臨界指数を推定。

3. 主要な結果と発見

A. 平均二乗速度の臨界スケーリング

GK 公式に含まれる平均二乗横速度 $\langle v_y^2 \rangle$ は、ジャミング近傍で臨界スケーリングに従うことが確認されました。

• ニュートン流体領域（$\phi < \phi_J$）: $\langle v_y^2 \rangle \sim \dot{\gamma}^2$ に比例し、粘度の発散と整合します。
• 準静的ジャミング領域（$\phi > \phi_J$）: $\langle v_y^2 \rangle \sim \dot{\gamma}$ に比例し、降伏応力のスケーリングと整合します。
• これらの振る舞いは、スケーリング関数 $F(x)$ によって統一的に記述され、臨界指数 $\psi \approx 5.39$ および $\lambda \approx 3.87$ が得られました。

B. 緩和時間と相関関数の振る舞い

速度自己相関関数 $C(\tau)$ の解析により、以下の 2 つの異なる領域が明らかになりました。

1. ジャミング以下、または高いせん断率（$\phi < \phi_J$ または $\dot{\gamma}$ 大）:
   • $C(\tau)$ は引き伸ばされた指数関数（stretched exponential）$C(\tau) = \exp[-(\tau/\tau^*)^\alpha]$ でよく記述されます。
   • 緩和時間 $\tau^*$ も臨界スケーリングに従い、指数 $\epsilon \approx 2.57$ が得られました。
2. ジャミング以上、かつ極めて低いせん断率（$\phi > \phi_J$ かつ $\dot{\gamma} \to 0$）:
   • $C(\tau)$ は引き伸ばされた指数関数では記述できず、長時間テール（power-law decay）$C(\tau) \sim (\dot{\gamma}\tau)^{-1}$ を示します。
   • この領域では、GK 公式の時間積分が 2 次元系において発散し、熱力学極限（$L \to \infty$）で拡散係数 $D$ が無限大になることを意味します。

C. 縮約 GK 公式と臨界指数の関係式

引き伸ばされた指数関数領域において、GK 公式を解析的に積分することで「縮約 GK 公式」$D = \frac{\tau^* \langle v_y^2 \rangle}{\alpha} \Gamma(1/\alpha)$ を導出しました。
これに各物理量の臨界スケーリングを代入することで、異なる物理量間の臨界指数を結びつけるスケーリング関係式を提案しました：
$$\nu + \psi - \epsilon = \lambda$$
ここで、$\nu$ は拡散係数の指数、$\psi$ は速度の指数、$\epsilon$ は緩和時間の指数、$\lambda$ はせん断率の指数です。数値シミュレーションの結果（$\nu \approx 0.939, \psi \approx 5.39, \epsilon \approx 2.57, \lambda \approx 3.87$）はこの関係式を高精度で満たしています。

D. 有限サイズ効果と長時間テールの関連

ジャミング以上の準静的領域で観測される $D$ の発散は、系サイズ $L$ に依存する有限サイズ効果（$D/\dot{\gamma} \sim L^\beta$）として以前から報告されていましたが、本研究ではこれを速度自己相関関数の長時間テールによる GK 公式の積分発散として解釈し、メカニズムを解明しました。

4. 意義と結論

• メカニズムの解明: 従来の研究が「長さスケール」に焦点を当てていたのに対し、本研究は「時間スケール（緩和ダイナミクス）」に基づいてせん断誘起拡散を説明しました。
• 2 次元と 3 次元の違い: 古典流体と同様に、2 次元のアサーマル系では長時間テールにより GK 公式が破綻（積分発散）しますが、3 次元系では発散しない可能性が示唆されています（3 次元での検証は今後の課題）。
• 理論的枠組みの確立: 臨界現象の枠組みを用いて、拡散係数、速度、緩和時間のスケーリングを統一的に記述する関係式を提案し、数値的に検証しました。
• 実用的意義: 食品、医薬品、パーソナルケア製品の製造など、ソフトマターの輸送特性を制御する上で、ジャミング近傍の拡散挙動を理解することは極めて重要です。

この論文は、ジャミング近傍の非平衡系における輸送現象を、緩和ダイナミクスと長時間相関の観点から定量的に理解するための重要な一歩を示しています。