[논문 리뷰] Statistical Mechanics of the self-gravitating gas: thermodynamic limit, phase diagrams and fractal structures

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

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)로, 모든 물리적 양이 이 단일 변수에 의존함을 보였다.
• 작은 변동의 함수행렬식을 사용하여 평균장 이론의 타당성을 평가하고 각 계열에서 임계점을 확인하였다.
• 국소 열역학적 성질(압력, 밀도, 음속)을 계산하고, p(r) = Tρ(r)가 국소적으로 성립함을 입증하여 유체역학 기술의 타당성을 확인하였다.
• 중력의 장거리 성질과 평균 자유로움에 미치는 영향을 분석하여, N→∞ 근처에서 유체역학 기술이 정확해짐을 보였다.

실험 결과

연구 질문

• RQ1 자중가스는 광역계, 정열계, 미세계계에서 잘 정의된 열역학적 한계를 가지는가?
• RQ2 상태방정식은 차원이 없는 매개변수 η에 어떻게 의존하며, 상전이에서의 임계 행동은 어떠한가?
• RQ3 중력 상전이의 성격은 무엇이며, 잼 불안정성과 정열계에서의 음의 비열과의 관계는 어떠한가?
• RQ4 세 계열에서의 변동은 평균장 이론의 타당성에 어떻게 영향을 미치며, 임계점은 어디에 위치하는가?
• 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에서 제0차 상전이가 발생하며, 이 전이에서 구리에르지가 증가한다. 이는 정열계와 대조된다.
• 평균장 상태방정식은 제1종 아벨 유형의 미분방정식을 만족하며, 정열계와 미세계계 양쪽에서 몬테카를로 시뮬레이션과 매우 잘 일치한다.
• 기본 상태에서 입자 분포는 허스트 소수 D ∈ (1,3)를 가지며, 밀도가 증가함에 따라 감소한다. 또한 상전이에서 평균 입자 간격이 급격히 감소한다.