Morrey's Derivation of Hydrodynamics from Statistical Mechanics: A Modern Exposition

技術要約：統計力学からの流体力学の導出におけるモリの理論：現代的解説

問題提起
取り組まれている中心的な問題は、相互作用する$N$個の粒子のミクロなニュートン力学から、マクロな流体の運動方程式（具体的には圧縮性オイラー方程式）を数学的に導出することである。速度論的な経路（ボルツマン・グラッド極限）は、希薄ガスの場合は流体方程式の導出に成功しているが、これは理想気体の状態方程式に限定される。モリによる1955年のプログラムは、異なるスケーリング・レジーム（$N\epsilon^3 = O(1)$）を採用し、制約付き位相空間分布の列を構成することで、このギャップを埋めようと試みている。本論文は、モリの議論を現代的に、かつ透明性の高い形で解説し、形式的な導出と未証明の仮説を厳密に分離し、残された正確な解析学的欠陥を特定するものである。

手法および設定
本論文では、$m_N$の質量を持つ$N$個の同一の点質量粒子からなる系を考える。粒子間の相互作用ポテンシャルは$\Phi_N(r) = \epsilon^{-s} \Phi(r/\epsilon)$であり、$\epsilon \to 0$（かつ$N \to \infty$）とする。スケーリング則は$m_N = D\epsilon^3$および$N m_N = O(1)$であり、これにより全質量を有限に保ちつつ、相互作用の範囲を微視的なスケールへと縮小させている。

手法は、以下の3つの明確なフェーズで進行する：

1. 形式的なバランス則（第2節）：
   不変な$N$粒子位相空間分布が存在すると仮定することで、簡約周辺分布に対するBBGKY階層を導出する。局所的な変数（粒子を中心とし、相対座標を$\epsilon$でスケーリングしたもの）を導入することにより、この階層を局所密度$\phi_{Nl}$を用いた形式に書き換える。この方程式を積分することで、質量、運動量、エネルギーに関する形式的なバランス則が得られる。極めて重要な点として、これらの方程式は閉じておらず、高次の周辺分布や相互作用項に依存しており、それらは極限分布によって決定される必要がある。

2. 制約付き測度の構成（第3–5節）：
   あらかじめ不変測度を仮定するのではなく、モリは制約多様体$M_N$上に支持される一連の測度$\mu_{M_N}$を構成する。この多様体は、セル平均の質量、運動量、およびエネルギーが、規定されたマクロな場$(\rho, u, e)$と一致するすべての微視的状態からなる。

• 測度は、制約のヤコビアンを通じて、リウヴィル測度を$M_N$へ投影することによって幾何学的に定義される。
• 本論文は、周辺密度$\phi_{Nl}$を、組合せ論的な因子と条件付き体積比へと明示的に分解する。
• 重要なステップは、相互作用エネルギーを含む「条件付き体積比」を用いてセルごとの制約を解くことである。これにより、残りのエネルギー予算に依存するカットオフ関数を持つギブス分布に類似した、有限$N$の公式（式5.22）が得られる。

3. 流体力学的極限とエルゴード仮定（第3, 4, 6–8節）：
   議論の核心は、以下の2つの未証明の柱に依存している：

• エルゴード型仮定（仮定3）： マクロな場の時間微分を、制約多様体$M_N$上の空間平均で置き換えられると仮定している。これは、粒子の軌道がマクロな進化に対して、制約多様体を十分に速く探索することを前提としている。
• ギブス形式への収束（定理3）： 有限$N$の周辺分布$\phi_{Nl}$が、相関関数$g_l$と正規化関数$C(\rho, b)$によって特徴付けられる極限のギブス型分布に収束すると仮定している。

主な貢献と結果

• スケーリング・レジームの明確化： 本論文は、モルのレジーム（$N\epsilon^3 = O(1)$）が、ボルツマン・グラッド極限（$N\epsilon^{d-1} = O(1)$）とは根本的に異なることを強調している。このレジームは液相の状態に適しており、圧力は多角形式に固定されず、相互作用ポテンシャルから導かれるため、より柔軟な状態方程式を許容する。
• 周辺分布の明示的な構成： 本論文は、有限$N$の周辺分布の導出（式5.22）を詳細に示し、制約付き体積比がいかにして関数$C(\rho, b)$とその微分によって支配される構造を導くかを明らかにしている。
• 状態方程式の導出： ギブス形式の解析を通じて、圧力$p$と内部エネルギー$\varepsilon$を、正規化関数$C(\rho, b)$および相関関数$g_2$を用いて明示的に導出している。
  • 圧力は $p = \frac{\rho}{2A} - \frac{3}{4A} \frac{\rho^2 C_\rho}{C}$ で与えられる（ここで$A$は逆温度パラメータ）。
  • この結果は、状態方程式が、理想気体の法則を強制されるのではなく、相互作用ポテンシャル$\Phi$を通じて$C$によって決定されることを示している。
• エントロピー構造： 本論文は、関数$C(\rho, b)$がエントロピー$S$を生成することを示している。具体的には、$A(d\varepsilon - p d(1/\rho))$の組み合わせが完全微分$dS$になることが示され、輸送則$S_t + u \cdot \nabla S = 0$ が導かれる。これにより、導出された系が圧縮性オイラー方程式のエントロピー保存則を満たすことが確認される。
• 相関関数の形式解： 第7節では、マイヤーのクラスター展開を用いて相関関数$g_l$の形式解を提供し、ギブス仮説の下で空間相関が再帰的に決定できることを示している。

特定された欠陥と限界
本論文の主張は控えめであり、モリの導出が現代的な意味で完全には厳密ではないことを明示している。本解説は、以下の決定的な欠陥を孤立させている：

1. 測度の不変性： 構成された測度$\mu_{M_N}$は、制約多様体上に投影されたフローの下では不変であるが、真のニュートン力学のフローの下では不変ではない。本論文は、その偏差が流体力学的極限において無視可能（$O(\epsilon)$）であると仮定する予想2に依拠している。
2. エルゴード仮説： 仮定3（時間微分を空間平均で置き換えること）は証明されていない。低次のモーメントについてはコンパクト性を用いた議論（定理2）により部分的な正当化が示されているが、高次のモーメント（特にエネルギー）の有界性は依然として未解決の問題である（予想1）。
3. 微分の収束： 相関関数$g_{Nkl}$とその微分（特に空間微分）の一様収束は、厳密には確立されていない。予想4は、極限方程式が成立するために、これらの微分が一様に有界であることを仮定している。
4. 単一セルの制限： 明示的な有限$N$の公式の導出は、規定された粒子が単一のセル内に存在するという単純化された仮定（予想3）に依存しており、セル間の相互作用を無視している。本論文は、この誤差は極限において消失すると主張しているが、それを証明してはいない。

意義
本研究の意義は、モリの1955年の議論を構造的に明確にし、アクセス可能にした「現代的解説」としての役割にある。本論文はヒルベルトの第6問題を解決すると主張するのではなく、むしろ困難さがどこにあるのかを正確に描き出すものである。制約付き測数の厳密な幾何学的構成と、ヒューリスティックなエルゴード仮定を分離することで、本論文は、もしこれらが解決されれば、一般的な状態方程式を持つ高密度流体に対する流体力学の導出を完結させることになる、明確に定義された一連の未解決問題（予想1–4）を特定している。これにより、このレジームにおける状態方程式の柔軟性が、相互作用ポテンシャルの直接的な帰結であることを示しており、これは速度論による硬直した理想気体の結果とは対照的である。