Elasticity from Extended Structural Dynamics: Stress, Strain, and Elastic Moduli as Emergent Collective Properties

技術要約：拡張構造力学からの弾性

問題提起
ナビエ・ストークス・フーリエ（NSF）方程式によって定式化される古典的な流体力学は、物質を構造を持たない点粒子の集合体とする理想化に基づいている。この仮定は根本的な限界をもたらす。すなわち、放物型輸送方程式は信号伝播速度を無限大に設定してしまい、衝撃波面には数値的安定性のために人工的な粘性が要求され、さらに、この理論には異方性輸送、非フーリエ熱現象、およびメンペラ効果のようなマルチタイムスケールの緩和効果を扱うメカニズムが欠如している。これまでの拡張手法（例：マイクロポーラー流体や拡張不可逆熱力学）は、これらの問題に対して現象論的に対処してきたが、それらの多くは第一原理からの厳密な導出を欠いているか、あるいはマクロな輸送係数と微視的な粒子の幾何学的形状および動力学を直接結びつけることに失敗している。

手法
本論文では、構成要素を配向、角運動量、および内部変形モードを有する空間的に広がった物体として扱うキネティック・フレームワークである**拡張構造力学（Extended Structural Dynamics: ESD）**を導入する。その手法は以下の通りである：

1. 拡張相空間： 本フレームワークは、位置（$\mathbf{r}$）、運動量（$\mathbf{p}$）、配向（$R \in SO(3)$）、物体固定系角運動量（$\mathbf{L}$）、および内部変形座標（$\{\xi_n, \pi_n\}$）を含む拡張された相空間 $T^*(\mathbb{R}^3 \times SO(3) \times \Xi)$ を仮定する。
2. キネティック方程式： 1粒子分布関数 $f(\mathbf{r}, t)$ の進化は、ハミルトニアン・ストリーミングのためのポアソン括弧と、二体衝突および構造的交換を考慮した衝突演算子 $C[f]$ を含む、拡張されたボルツマン方程式によって支配される。
3. チャップマン・エンスコグ展開： 局所平衡の周りで、クヌーセン数のべき乗を用いた系統的な展開が行われる。
4. 閉鎖（クロージャ）： 導出には、モード依存のBGK（Bhatnagar–Gross–Krook）緩和モデルを利用し、並進、回転、および内部自由度に対してそれぞれ異なる緩和時間（$\tau_{trans}, \tau_{rot}, \tau_{int}$）を割り当てる。
5. フレーム不変性： 導出過程において、$SO(3)$ ストリーミング項から自然に現れるコローショナル（ジャウマン）微分を用いることで、構成則が客観的（フレーム不変）であることを明示的に保証する。

主要な貢献と結果

• 一般化された流体力学方程式： 展開により、主に2つの双曲型・放物型輸送法則が得られる：

  1. スピン動力学： 配向緩和を流体の渦度へと結合させる、スピン密度 $\mathbf{s}$（粗視化された角運動量）の動力学方程式。これには、構造拡散項と、渦度と平均粒子角速度との差に比例する緩和項が含まれる。
  2. 熱流束緩和： 熱流束 $\mathbf{q}$ に対するカテトー型の方程式を導入し、有限の熱緩和時間 $\tau_q$ を与える。これにより、双曲型熱方程式が導かれ、有限の熱波速度が示唆される。

• 創発的輸送係数：

  • 構造粘性（$\eta_{struct}$）： 有効粘性は、並進成分と構造成分の和であることが示される。$\eta_{struct}$ は粒子の幾何学的形状（球体では消失する）、熱エネルギー、および配向緩和時間 $\tau_{rot}$ に依存してスケールする。
  • 回転熱伝導率（$\kappa_{rot}$）： 回転への寄与に対する同様のスケーリングが導出される。
  • 有限の伝播速度： 本フレームワークは、粒子の有限の慣性モーメント（$I > 0$）から直接導かれる、運動量および熱の両方における有限の信号速度を予測する。

• 衝撃波の正則化： 本論文は、有限の配向慣性が圧縮中の瞬時な再配向を防ぐことを実証している。これにより、衝撃波面の固有の正則化がもたらされ、有限の衝撃波幅 $\delta \sim \sqrt{D_s \tau_{rot}}$ が予測される。非対称分子の場合、この幅は古典的な予測よりも大幅に広く（例：SF$_6$ で $\sim 7$ nm）、広幅化した衝撃波プロファイルの実験的観察結果を説明できる可能性がある。

• メンペラ効果： 内部エネルギーを並進温度（$T_{trans}$）と回転温度（$T_{rot}$）に分解することにより、本論文は二温度緩和モデルを導出する。これにより、メンペラ効果（より熱い系がより冷たい系よりも速く冷却される現象）の基準を確立する：システムは内部平衡から遠い状態（$T_{trans} \gg T_{rot}$）から開始し、急速な内部平衡化が追加の冷却チャネルを提供する必要がある。

  • 定量的推定： コロイダル楕円体について、本モデルはクロスオーバー時間 $t_{cross} \approx 12$ ms を予測しており、これは現在の光ピンセット実験の分解能の範囲内にある。

意義と主張
本論文は、ESDを、マイクロポーラー流体理論や拡張不可逆熱力学（EIT）のような構成的なアプローチとは対照的な、生成的なフレームワークとして位置づけている。

• 導出 vs 仮定： スピン場を構成則として仮定するマイクロポーラー理論とは異なり、ESDはスピン場を分布関数のモーメントとして導出し、輸送係数を粒子の幾何学的形状と動力学から導出する。
• 存在論的必然性： 有限の緩和時間と双曲型の性質は、単なるアドホックな修正ではなく、実在する粒子が有限の慣性モーメント（$I > 0$）を持つことから生じる存在論的な必然性として提示されている。
• 主張の謙抑性： 著者らは、自らの結果を以下の3つのレベルに明確に区別している：
  1. 導出された結果： 定性的な特徴（スピン・バランスの存在、双曲型・放物型構造、フレーム不変性）は堅牢であり、BGKによる閉鎖に依存しない。
  2. スケーリング推定： 関数的な依存関係（例：$\eta_{struct} \sim a^2 \tau_{rot}$）は物理的に意味があるが、BGK近似のため厳密な数値的前係数は欠いている。
  3. ヒューリスティックな推定： 特定の数値的予測（例：$t_{cross} \approx 12$ ms, $\delta \approx 7$ nm）は、リアプノフ不安定性を伴う $\tau_{rot}$ の推測に基づく加法的な公式に依存している。これらは、パラメータフリーの予測ではなく、検証可能な帰結を示唆するためのオーダー・オブ・マグニチュード（桁違い）の例示として提示されている。

結論として、ESDは、構造的効果を自然に組み込んだ非平衡熱力学の統一的な第一原理フレームワークを提供し、異方性輸送、有限の信号速度、および古典的な流体力学では捉えられないマルチスケール緩和現象を説明する経路を提示している。