[論文レビュー] Statistical Mechanics of the self-gravitating gas: thermodynamic limit, phase diagrams and fractal structures

de Vega, H. J., N. Sánchez|arXiv (Cornell University)|May 27, 2005

Statistical Mechanics and Entropy参考文献 2被引用数 3

ひとこと要約

本稿は、モンテカルロシミュレーション、平均場理論、および低密度展開を用いて、熱平衡状態の自己重力的ガスの包括的な統計力学枠組みを提示する。全系において熱力学的極限が存在することを確立し、η ≈ 1.515 における重力的相転移（ジェイムズ不安定性に起因）を特定し、物理的量が次元なしパラメータ η = Gm²N/(V¹ᐟ³T) のみに依存することを明らかにした。圧縮率および比熱の臨界行動は、収縮ダイナミクスを示唆する。

ABSTRACT

We provide a complete picture to the self-gravitating non-relativistic gas at thermal equilibrium using Monte Carlo simulations, analytic mean field methods (MF) and low density expansions. The system is shown to possess an infinite volume limit in the grand canonical (GCE), canonical (CE) and microcanonical (MCE) ensembles when (N, V) -> infty, keeping N/V^{1/3} fixed. We compute the equation of state (we do not assume it as is customary in hydrodynamics), as well as the energy, free energy, entropy, chemical potential, specific heats, compressibilities and speed of sound; we analyze their properties, signs and singularities. All physical quantities turn out to depend on a single variable eta = G m^2 N /[V^{1/3} T} that is kept fixed in the N -> infty and V -> infty limit. The system is in a gaseous phase for eta < eta_T and collapses into a dense object for eta > eta_T in the CE with the pressure becoming large and negative. At eta \simeq eta_T the isothermal compressibility diverges and the gas collapses. Our Monte Carlo simulations yield eta_T \simeq 1.515. We find that PV/[NT] = f(eta). The function f(eta) has a second Riemann sheet which is only physically realized in the MCE. In the MCE, the collapse phase transition takes place in this second sheet near eta_{MC} = 1.26 and the pressure and temperature are larger in the collapsed phase than in the gaseous phase. Both collapse phase transitions (in the CE and in the MCE) are of zeroth order since the Gibbs free energy has a jump at the transitions.

研究の動機と目的

• 自己重力的ガスが巨大集団、正準、マイクロ正準系において熱力学的極限を有するかどうかを確立すること。
• 事前に形式を仮定せずに、状態方程式および熱力学的性質（エネルギー、エントロピー、比熱、圧縮率）を特定すること。
• ジェイムズ不安定性および収縮に関連する重力的相転移を特定・特徴づけること、特に次元なしパラメータ η との関係を明らかにすること。
• 各系において、ゆらぎおよび関数的デターミナントを計算することで、平均場理論の妥当性を評価すること。
• 流体力学的および局所的状態方程式を分析し、収縮相におけるフラクタル粒子分布の出現を評価すること。

提案手法

• 正準およびマイクロ正準系を数値的に研究するために、立方体幾何におけるモンテカルロシミュレーションを用いた。
• 球対称性を仮定した平均場理論を適用し、状態方程式および熱力学的量の解析的解を導出した。
• 状態方程式を PV/(NT) = f(η) として導出し、ここで η = Gm²N/(V¹ᐟ³T) であり、すべての物理的量がこの1つの変数に依存することを示した。
• 小揺らぎの関数的デターミナントを用いて、平均場理論の妥当性を評価し、各系における臨界点を特定した。
• 局所的熱力学的性質（圧力、密度、音速）を計算し、p(r) = Tρ(r) が局所的に成り立つことを示した。これにより、流体力学的記述が正当化された。
• 重力の長距離性とその平均自由行程への影響を分析し、N→∞ の極限において流体力学的記述が正確に成立することを示した。

実験結果

リサーチクエスチョン

• RQ1 自己重力的ガスは、巨大集団、正準、マイクロ正準系において、明確な熱力学的極限を有するか？
• RQ2 状態方程式は次元なしパラメータ η にどのように依存するか？また、相転移における臨界行動は何か？
• RQ3 重力的相転移の性質は何か？また、ジェイムズ不安定性およびマイクロ正準系における負の比熱とどのように関連するか？
• RQ4 3つの系におけるゆらぎは、平均場理論の妥当性にどのように影響するか？臨界点はどこにあるか？
• RQ5 自己重力的ガスに対して流体力学的記述を正当化できるか？また、熱力学的極限において局所的状態方程式は p(r) = Tρ(r) に簡略化されるか？

主な発見

• 自己重力的ガスは、N→∞ かつ V→∞ の極限で、N/V¹ᐟ³ を固定した条件下で、巨大集団（GCE）、正準（CE）、マイクロ正準（MCE）のすべての系において熱力学的極限を示す。
• 状態方程式は PV/(NT) = f(η) であり、すべての熱力学的量が次元なしパラメータ η = Gm²N/(V¹ᐟ³T) のみに依存する。
• 正準系において、η_T ≈ 1.515 で相転移が発生し、等温圧縮率が発散する。これはジェイムズ不安定性を示唆する。
• マイクロ正準系では、η_MC ≈ 1.26 でゼロ次相転移が発生し、遷移でギブズ自由エネルギーが増加する。これは正準系とは対照的である。
• 平均場状態方程式は1階のアーベル型微分方程式を満たし、正準系およびマイクロ正準系の両方でモンテカルロシミュレーションと極めて良好に一致する。
• 基底状態における粒子分布は、ハウスドルフ次元 D ∈ (1,3) のフラクタル的挙動を示し、密度が増加するにつれて減少する。また、相転移付近で平均粒子間隔が著しく低下する。