[논문 리뷰] Ionization Modeling Astrophysical Gaseous Structures. I. The Optically Thin Regime

Christopher W. Churchill, E. S. Klimek외 2명|arXiv (Cornell University)|2014. 09. 02.

Galaxies: Formation, Evolution, Phenomena참고 문헌 12인용 수 4

한 줄 요약

이 논문은 광학적으로 얇은 천체물리학적 기체에서 광이온화, 충돌 이온화, 재결합 과정을 사용하여 평형 이온화 분율을 계산하는 빠르고 정확한 코드인 hartrate를 제시한다. 다양한 밀도와 온도 범위에서 핵심 이온들(H, He, C IV, O VI, Mg II)에 대해 Cloudy 13.03과 뛰어난 일치를 보이며, 이온화 보정 값의 차이가 ≤0.05 dex 이내로 유지되어 우주론적 시뮬레이션에서 효율적인 흡수선 분석을 가능하게 한다.

ABSTRACT

We present a code for modelling the ionization conditions of optically thin astrophysical gas structures. Given the gas hydrogen density, equilibrium temperature, elemental abundances, and the ionizing spectrum, the code solves the equilibrium ionization fractions and number densities for all ions from hydrogen to zinc. The included processes are photoionization, Auger ionization, direct collisional ionization, excitation auto-ionization, charge exchange ionization, two-body radiative recombination, dielectronic recombination, and charge exchange recombination. The ionizing spectrum can be generalized to include the ultraviolet background (UVB) and/or Starburst99 stellar populations of various masses, ages, metallicities, and distances. The ultimate goal with the code is to provide fast computation of the ionization conditions of gas in N-body + hydrodynamics cosmological simulations, in particular adaptive mesh refinement codes, in order to facilitate absorption line analysis of the simulated gas for comparison with observations. To this end, we developed a method to linearize the rate equations and efficiently solve the rate matrix with a minimum number of iterations. Comparing the code to Cloudy 13.03 (Ferland), we find that the hydrogen and helium ionization fractions and the metal species ionization corrections are in excellent agreement. We discuss the science drivers and plans for further development of the ionization code to a full radiative hydrodynamic routine that can be employed for processing the simulations. A stand-alone version of the code has been made publicly available (http://astronomy.nmsu.edu/cwc/Software/Ioncode/).

연구 동기 및 목표

• N-체 + 유체역학적 우주론적 시뮬레이션, 특히 적응 격자 구조(AMR) 코드에 사용할 수 있는 빠르고 정확한 이온화 코드를 개발하는 것.
• 대규모 시뮬레이션된 기체를 효율적으로 후처리하여 관측된 쿼라서 흡수선 데이터와 비교할 수 있는 합성 흡수선을 생성하는 것.
• 실제 이온화 스펙트럼(예: UVB 및 Starburst99 항성 집단 포함)을 사용하여 은하주위 매질(CGM)의 이온화 조건을 모델링하는 것.
• 최소한의 입력 요구 사항으로 관측 이온화 모델링을 위한 독립형 도구를 제공하는 것.
• 자기 차폐 및 주파수 의존성 복사 전달을 고려한 전체 복사 유체역학 시뮬레이션에의 향후 통합을 위한 기반을 마련하는 것.

제안 방법

• 계산 효율성을 확보하기 위해 최소한의 반복 수를 사용하여 이온화 평형 상태의 선형화된 비율 방정식 시스템을 해결하는 코드.
• 8개의 핵심 물리 과정을 포함: 광이온화, 아우거 이온화, 직접 충돌 이온화, 자극화 자가이온화, 전하 교환 이온화, 이체 방사 재결합, 이전자 재결합, 전하 교환 재결합.
• 이온화 스펙트럼은 적색편이(광학적 Haardt & Madau UVB용)와/또는 Starburst99 모델(지정된 질량, 연령, 금속성, 거리 포함)로 정의된다.
• 이온화 평형과 광학적으로 얇은 조건을 가정하여 H에서 Zn까지 모든 이온의 이온화 분율과 수밀도를 계산한다.
• 광이온화, 재결합, 충돌 이온화 단면적에 대해 Verner, AMDPP, ORNL 데이터베이스에서 공개된 원자 데이터를 사용한다.
• 계산 반복 수를 최소화하면서도 높은 정확도를 유지하기 위해 비율 행렬 방정식에 선형화 기법을 적용한다.

실험 결과

연구 질문

• RQ1 hartrate는 광학적으로 얇은 영역에서 기준 코드인 Cloudy 13.03와 비교해 이온화 분율과 보정 값을 얼마나 정확히 재현할 수 있는가?
• RQ2 밀도와 온도 조건이 다양할 때, 핵심 흡수선 이온들(Mg II, C IV, O VI)의 이온화 보정 값이 hartrate와 Cloudy 사이에서 어느 정도 일치하는가?
• RQ3 hartrate는 대규모 우주론적 시뮬레이션의 후처리에 적합한 방식으로 이온화 상태를 효율적으로 계산할 수 있는가?
• RQ4 온도와 밀도가 CGM 내 수소, 헬륨 및 금속 종의 이온화 균형에 어떤 영향을 미치는가?
• RQ5 UVB와 Starburst99 항성 집단과 같은 복합 이온화 원천을 코드가 얼마나 잘 처리하는가?

주요 결과

• hartrate는 수소와 헬륨 이온화 분율에서 Cloudy 13.03와 뛰어난 일치를 보이며, 관측 불확도 수준 이내의 차이를 보인다.
• 밀도 및 온도 영역 −7 ≤ log n_H ≤ 0 및 3 ≤ log T ≤ 6 내에서 Mg II, C IV, O VI의 이온화 보정 값은 ≤0.05 dex 이내로 일치하여 높은 일관성을 보인다.
• log T ≥ 5 조건에서 중성 헬륨 분율에 대해 약 2–3배의 불일치가 관측되었으나, 이는 재결합 데이터의 알려진 불확도로 기인한 것으로 보인다.
• 비율 행렬의 선형화 덕분에 코드는 계산적으로 효율적이며, 대규모 시뮬레이션 출력에서 이온화 상태를 신속하게 계산할 수 있다.
• hartrate의 독립형 버전은 공개되어 있으며, 쿼라서 흡수선 연구를 포함한 관측 데이터에 성공적으로 적용된 바 있다.
• 향후 개선 사항으로 광학적으로 두꺼운 구름에 대한 자기 차폐 및 주파수 의존성 복사 전달 기능이 추가될 예정이며, 업데이트 사항은 온라인으로 제공될 것이다.