[논문 리뷰] PRISM: A Non-Equilibrium, Multiphase Interstellar Medium Model for Radiation Hydrodynamics Simulations of Galaxies

Harley Katz, Shenghua Liu외 8명|arXiv (Cornell University)|2022. 11. 09.

Astrophysics and Star Formation Studies인용 수 6

한 줄 요약

PRISM는 RAMSES-RTZ 코드에서 실시간 다주파수 복사장 운반을 통합한 115종의 화학 네트워크를 갖춘 비평형 다상성 간성간매체 모델이다. 이 모델은 광범위한 금속성에서 핵심 냉각 및 가열 과정을 정확히 기록하며, 비평형 효과로 인해 전자 분율이 평형 상태에서 3배 이상 벗어나는 것으로 드러나, ISM 열역학 및 관측 가능한 발광선에 상당한 영향을 미친다.

ABSTRACT

We introduce the PRISM interstellar medium (ISM) model for thermochemistry and its implementation in the RAMSES-RTZ code. The model includes a non-equilibrium primordial, metal, and molecular chemistry network for 115 species coupled to on-the-fly multifrequency radiation transport. PRISM accurately accounts for the dominant ISM cooling and heating processes in the low-density regime (i.e. $ρ<10^5\ { m cm^{-3}}$), including photoheating, photoelectric heating, H$_2$ heating/cooling, cosmic-ray heating, H/He cooling, metal-line cooling, CO cooling, and dust cooling (recombination and gas-grain collisions). We validate the model by comparing 1D equilibrium simulations across six dex in metallicity to existing 1D ISM models in the literature. We apply PRISM to high-resolution (4.5 pc) isolated dwarf galaxy simulations that include state-of-the-art models for star formation and stellar feedback to take an inventory of which cooling and heating processes dominate each different gas phase of a galaxy and to understand the importance of non-equilibrium effects. We show that most of the ISM gas is either close to thermal equilibrium or exhibits a slight cooling instability, while from a chemical perspective, the non-equilibrium electron fraction is often more than three times higher or lower than the equilibrium value, which impacts cooling, heating, and observable emission lines. Electron enhancements are attributed to recombination lags while deficits are shown to be due to rapid cosmic-ray heating. The PRISM model and its coupling to RAMSES-RTZ is applicable to a wide variety of astrophysical scenarios, from cosmological simulations to isolated giant molecular clouds, and is particularly useful for understanding how changes to ISM physics impact observable quantities such as metallic emission lines.

연구 동기 및 목표

• 은하 내 화학, 복사, 열역학 간의 복잡한 결합을 정확히 기록하는 물리적으로 정확한 비평형 간성간매체 모델을 개발하는 것.
• 1D 평형 모델의 한계를 해결하기 위해 자가일관된 복사장과 화학을 갖춘 3D 및 시간에 따라 변화하는 ISM 시뮬레이션을 가능하게 하는 것.
• 특히 전자 분율에 미치는 비평형 효과의 영향을 정량화하여 ISM 냉각, 가열 및 관측 가능한 발광선에 미치는 영향을 분석하는 것.
• 금속성 범위 6디엑스에 걸쳐 기존의 1D 평형 ISM 모델과의 정합성을 검증하는 것.
• 고해상도 고립 저금속성 은하 시뮬레이션에 모델을 적용하여 다양한 기체 상에서의 냉각 및 가열 과정의 지배적 요인을 평가하는 것.

제안 방법

• 115종의 종류(예: H2, CO, 금속 이온 포함)를 포함한 비평형 원시성, 금속성 및 분자 화학 네트워크를 구현한다.
• RAMSES-RTZ 코드 내 실시간 다주파수 복사장 운반과 화학 네트워크를 결합하여 복사장과 기체 상태 간의 자가일관된 피드백을 가능하게 한다.
• 주요 ISM 과정 포함: 광가열, 광전자 가열, H2의 가열/냉각, 宇宙선 가열, H/He 냉각, 금속선 냉각, CO 냉각, 재결합 및 기체-먼지 충돌을 통한 먼지 냉각.
• 정확한 냉각 계산을 위해 CLOUDY 캘리브레이션을 거친 미세구조 충돌률을 사용하여 다수의 충돌 파artner(H, H+, H2, e−, He 등)에 적용한다.
• 금속성 범위 6디엑스에 걸쳐 문헌 모델(Koyama & Inutsuka 2000; Wolfire et al. 2003; Bialy & Sternberg 2019; Kim et al. 2022)과의 1D 평형 기준 비교를 수행한다.
• 은하 형성과 항성 피드백을 포함한 3D 고해상도(4.5 pc) 고립 저금속성 은하 시뮬레이션(G8 및 G9)에 모델을 적용하여 상별 냉각/가열 및 비평형 역학을 연구한다.

실험 결과

연구 질문

• RQ1 은하 규모 시뮬레이션에서 전자 분율의 비평형 효과가 ISM 열역학 및 냉각/가열 속도에 어떤 영향을 미치는가?
• RQ2 실제 피드백 조건이 적용된 고립 저금속성 은하에서 다양한 기체 상(예: CNM, WNM, HIM)에서 냉각 및 가열 과정 중 어떤 것이 지배적인가?
• RQ3 시간에 따라 변화하는 복사장과 화학을 포함한 3D 시뮬레이션은 1D 평형 모델과 비교해 원자-분자 전이 및 상 구조 예측에서 얼마나 다를까?
• RQ4 금속성과 복사장의 변화가 전자 분율의 평형 값에서의 벗어남 정도에 어떤 영향을 미치는가?
• RQ5 비평형 전자 분율이 특히 저금속성 환경에서 관측 가능한 발광선에 어떤 영향을 미치는가?

주요 결과

• 비평형 전자 분율은 재결합 지연과 빠른 宇宙선 가열로 인해 평형 값보다 3배 이상 초과하거나 저하될 수 있다.
• 모델은 대부분의 ISM 기체가 열적 평형에 가까운 상태이거나 경미한 냉각 불안정성만을 보이며, 이는 비평형 효과가 상당히 중요하지만 전체 상 구조에서는 지배적이지 않음을 시사한다.
• 3D 시뮬레이션에서 FUV 복사장의 변화, 자기보호, 이온화 피드백를 고려함으로써 원자-분자 전이가 1D 모델보다 빠르게 발생한다.
• 냉각 속도가 재결합 시간스케일을 앞서는 재결합 지연으로 전자 분율 증가가 발생하며, 반면 재결합에 비해 빠른 宇宙선 가열로 전자 분율 감소가 발생한다.
• PRISM 모델은 금속성 범위 6디엑스에 걸쳐 평형 온도-밀도 관계를 정확히 재현하며, 밀도 ≤10³ cm⁻³ 범위에서는 최대 10% 이내의 미미한 차이를 보인다.
• 이 모델은 저금속성 은하에서 금속 농도 측정에 온도 변동의 영향을 설명하는 데 성공적으로 적용되었으며(Cameron et al. 2022), 관측 진단 해석에 유용함을 입증하였다.