[논문 리뷰] Extended Hot Gas Halos Around Starburst Galaxies

Kohji Tomisaka, Joel N. Bregman|arXiv (Cornell University)|1992. 09. 12.

Astrophysics and Star Formation Studies인용 수 6

한 줄 요약

이 논문은 항성 폭발 은하 주변의 확장된 뜨거운 가스 허브를 유체역학 시뮬레이션을 사용해 모델링하며, 초신성에 의해 유도되는 바람이 충격파를 형성하고 세 단계로 진화하는 것을 보여준다: 디스크 확장, 디스크 탈출, 그리고 허브 확산. 시뮬레이션은 GINGA 대역에서 3×10³⁹에서 10⁴⁰ erg/s의 X선 복사원도와 Einstein/ROSAT 대역에서 약 10⁴¹ erg/s의 X선 복사원도를 예측하며, 허브 밀도가 10⁻³ cm⁻³일 경우 50 Myr 내로 기체가 40–50 kpc까지 도달한다.

ABSTRACT

Reanalysis of Einstein IPC data and new observations from the GINGA LAC indicate the presence of extended X-ray emission (10-50 kpc) around the starburst galaxy M82. Here we model this emission by calculating numerical hydrodynamic simulations of the starburst event to much later times and larger scales than previously considered. For our models, we adopt a supernova rate of 0.1 ${ m yr}^{-1}$, and an extended low-density static halo that is bound to the galaxy. There are three stages to the evolution of the wind-blown bubble and the propagation of the shock front: the bubble expands in an almost uniform density disk gas, with a deceleration of the shock front ($t \alt $ 3.6 Myr); breakout from the disk and the upward acceleration of the shock front (3.6 Myr $\alt t \alt$ 18 Myr); propagation into the halo, leading to a more spherical system and shock deceleration (18 Myr $\alt t$). For a halo density of $10^{-3} { m cm}^{-3}$, the outflow reaches a distance of 40-50 kpc from the center of the starburst galaxy in 50 Myr. We calculate the time evolution of the X-ray luminosity and find that the extended starburst emits $3 imes 10^{39}\lcgs$ to $10^{40}\lcgs$ in the GINGA LAC band and $\sim 10^{41}\lcgs$ in the Einstein or ROSAT HRI band. The degree of the ionization equilibrium in the outflow and its effect on the iron K$α$ line emission are discussed.

연구 동기 및 목표

• 관측된 M82 은하 주변 10–50 kpc 범위에서의 확장된 X선 방사원의 기원을 설명하는 것.
• 은하 허브 환경에서 초신성에 의해 유도되는 바람의 장기적 진화를 조사하는 것.
• 气泡 형성에서 허브 확산에 이르는 충격파의 유체역학적 진화를 모델링하는 것.
• 뜨거운 기체 분출의 X선 복사원도와 이온화 상태를 정량화하는 것.
• 이온화 평형이 철 Kα 선 방출에 미치는 영향을 평가하는 것.

제안 방법

• 초신성 발생률이 0.1 yr⁻¹인 50 Myr 동안의 항성 폭발에 의해 유도되는 바람의 진화를 추적하는 수치적 유체역학 시뮬레이션을 수행한다.
• 은하에 결속된 정적이고 낮은 밀도의 허브(10⁻³ cm⁻³)를 가정하며, 초기 조건으로 균일한 밀도의 디스크를 포함한다.
• 시뮬레이션은 세 단계를 포착한다: 충격파가 디스크 내에서 둔화되는 단계(t < 3.6 Myr), 디스크에서 탈출하는 단계(3.6 Myr < t < 18 Myr), 그리고 허브로 구형 확산되는 단계(t > 18 Myr).
• X선 복사원도는 GINGA LAC 및 Einstein/ROSAT HRI의 스펙트럼 대역을 사용하여 시간에 따라 계산된다.
• 이온화 평형과 철 Kα 선 방출을 모델링하여 기체 분출의 스펙트럼 특징을 평가한다.
• 시뮬레이션은 대규모 공간 스케일(최대 50 kpc)과 늦은 시간(50 Myr)까지 수행되어 이전 연구를 초월한다.

실험 결과

연구 질문

• RQ1 M82 주변 10–50 kpc에서 관측된 확장된 X선 방사원의 원인은 무엇인가?
• RQ2 충격파는 시간이 지남에 따라 디스크에서 은하 허브로 어떻게 진화하는가?
• RQ3 다양한 에너지 대역에서 뜨거운 기체 분출의 예측된 X선 복사원도는 무엇인가?
• RQ4 이온화 평형은 철 Kα 선 방출 강도에 어떤 영향을 미치는가?
• RQ5 실제 허브 밀도 조건 하에서 50 Myr 시점에 뜨거운 기체 허브의 최대 확장 범위는 얼마인가?

주요 결과

• 유체역학 시뮬레이션은 초신성에 의해 유도되는 바람이 50 kpc까지 확장되는 뜨거운 가스 허브를 생성할 수 있음을 보여주며, 관측된 확장된 X선 방사원을 재현한다.
• 충격파는 세 가지 명백한 단계로 진화한다: 디스크 내 둔화, 디스크 탈출, 그리고 허브로의 구형 확산으로, 전이 시점은 약 3.6 Myr와 18 Myr이다.
• 허브 밀도가 10⁻³ cm⁻³일 경우, 기체 분출은 50 Myr 이내에 은하 중심에서 40–50 kpc까지 도달한다.
• 예측된 X선 복사원도는 GINGA LAC 대역에서 3×10³⁹에서 10⁴⁰ erg/s이며, Einstein 또는 ROSAT HRI 대역에서는 약 10⁴¹ erg/s이다.
• 기체 분출에서는 이온화 평형이 완전히 확립되지 않아, 예측된 철 Kα 선 방출에 영향을 미친다.
• 이 모델은 M82의 확장된 허브의 형태와 X선 특성을 일관되게 설명한다.