H3 이온의 에너지와 분배함수 저밀도 플라스마에서의 전량 양자 통계

초록

이 논문은 경로 적분 몬테카를로(PIMC) 방법을 이용해 저밀도 환경에서 H₃⁺ 이온의 전량 양자 통계량을 160 K부터 15 000 K까지 계산한다. 총 에너지, 분배함수, 자유에너지, 엔트로피, 열용량을 직접 얻고, 4 000 K 이상에서의 해리·재결합 평형과 10 000 K 이상에서의 열이온화 현상을 온도와 밀도에 따라 분석한다. 또한 해리 전 온도 구간에 적용 가능한 새로운 분석적 분배함수 모델을 제시하고, 기존 Neale‑Tennyson 데이터와 비교한다.

상세 분석

본 연구는 H₃⁺ 이온을 구성하는 3개의 양성자와 2개의 전자를 하나의 주기적 초소형 셀에 배치하고, 경로 적분 양자 통계의 기본식 Z = Tr e^{‑βĤ} 를 페어 근사(pair approximation)와 메트로폴리스 알고리즘을 결합한 PIMC 시뮬레이션으로 직접 평가한다. 타임 스텝 τ = 0.03 E_H⁻¹ 로 설정하고, 트롯터 수 M을 온도에 따라 조절해 β = Mτ/k_B 를 만족하도록 하였다. 세 가지 부피(300 a₀³, 100 a₀³, 50 a₀³)에서 각각 평균 질량 밀도 1.3×10⁻⁶, 3.4×10⁻⁵, 2.7×10⁻⁴ g cm⁻³ 를 구현했으며, 이는 저밀도 한계에서 유한 크기 효과를 최소화한다.

시뮬레이션 결과는 온도 4 000 K 이하에서는 H₃⁺ 이온이 거의 완전한 비선형 구조를 유지하며, 총 에너지는 0 K에서의 비adiabatic 기저 상태 에너지(≈‑1.169 E_H)와 일치한다. 4 000 K~10 000 K 구간에서는 온도와 밀도에 따라 H₂ + H⁺, H₂⁺ + H, 2H + H⁺ 등 다양한 파편으로의 해리가 진행되며, 에너지 분포는 히스토그램에서 여러 피크로 명확히 구분된다. 특히 높은 밀도에서는 해리 평형이 억제되고, 낮은 밀도에서는 10 000 K 이상에서 수소 원자와 전자의 열이온화가 관측된다.

분배함수 Z(T)를 직접 얻기 위해 에너지 기대값 ⟨E⟩ = k_B T² ∂f/∂T (f = −F/k_B T) 를 이용, 경험적 형태
⟨E⟩ = k_B T²