이중 제트 폭발이 만든 은하 중심 버블의 비밀

초록

본 논문은 은하 중심에서 두 차례에 걸쳐 발사된 AGN 제트가 각각 전방 충격파를 형성해 eROSITA 버블과 Fermi 버블을 만들었다는 시나리오를 3D MHD 시뮬레이션으로 검증한다. 첫 번째 제트는 15 Myr 전, 두 번째는 5 Myr 전에 발생했으며, 각각 약 3.5 × 10⁵⁵ erg와 1.1 × 10⁵⁵ erg의 에너지를 주입한다. 시뮬레이션은 두 버블의 형태, X‑ray 표면 밝기, O VIII/O VII 비율, 라디오 리지, 감마선 방출을 모두 재현하며, 감마선은 두 번째 충격에서 현지 가속된 전자에 의해 발생한다는 결론을 제시한다.

상세 분석

이 연구는 은하 중심의 활동성 핵(AGN) 제트가 두 차례에 걸쳐 발사된다는 가정 하에, 3차원 자기유체역학(MHD) 시뮬레이션을 수행하였다. 첫 번째 제트는 t = 0 Myr에 시작해 1 Myr 동안 지속되며 총 에너지 E₁ = 3.46 × 10⁵⁵ erg를 주입한다. 10 Myr 후 두 번째 제트가 동일한 지속시간으로 E₂ = 1.10 × 10⁵⁵ erg를 투입한다. 두 제트 모두 은하 회전축을 따라 수직으로 발사되며, 초기 플라즈마는 동역학에너지 우세(kinectic‑energy‑dominated) 상태이다. 시뮬레이션 결과, 첫 번째 제트는 외부 전방 충격을 형성해 반경 ≈ 18 kpc(북쪽 ≈ 15 kpc)의 eROSITA 버블을 만들고, 두 번째 제트는 내부 전방 충격을 통해 ≈ 10 kpc 높이의 Fermi 버블을 생성한다. 두 충격 사이의 상호작용은 내부 가스의 재압축과 재가열을 일으켜, 외부 버블 내부에 온도 0.3–0.4 keV, 밀도 10⁻³ cm⁻³ 수준의 얇은 쉘을 만든다. 전자 수밀도와 내부 에너지 밀도는 충격 전후로 크게 증가하며, 자기장 강도는 ∼ 1 μG까지 증폭된다.

X‑ray 합성 맵은 ROSAT 0.11–2.04 keV 대역에서 관측된 고위도 밝기 감소와 저위도 흡수 효과를 재현한다. 특히 0.44–1.21 keV와 0.73–2.04 keV 대역에서는 두 개의 별도 방출 층이 나타나며, 이는 두 전방 충격이 각각 eROSITA와 Fermi 버블의 경계를 정의한다는 것을 시각적으로 확인시킨다. O VIII/O VII 라인 비율 맵은 버블 내부에서 1.5 이하의 마하수(M < 1.5)를 보이며, 약 15 Myr의 동역학적 연령을 지지한다.

감마선 방출은 전자들이 두 번째 충격 전면에서 현지 가속(in‑situ acceleration)된다고 가정하고, 역컴프턴(ICS) 과정을 통해 1–5 GeV 대역의 관측된 표면 밝기를 재현한다. 전자 에너지 밀도는 열에너지 밀도와 비례하도록 설정했으며, 전자 스펙트럼 지수 p ≈ 2.2–2.4를 채택했다. 이 모델은 Fermi 버블의 급격한 가장자리와 거의 균일한 표면 밝기를 자연스럽게 설명한다.

라디오(30 GHz) 편광 시뮬레이션은 충격 압축된 자기장이 토로이달·폴러드 성분을 결합해 관측된 대규모 반대 대칭 편광 구조를 재현한다. 또한, 첫 번째 제트에서 방출된 오래된 전자는 현재 라디오 주파수(≈ 2 GHz)에서 SPASS 리지 형태의 광역 구조를 형성한다는 예측을 제시한다.

전반적으로 이중 제트 시나리오는 단일 폭발 모델이 설명하기 어려운 eROSITA와 Fermi 버블의 크기 차이, 연령 차이, 다중 파장 관측 특성을 일관되게 통합한다.