AGN에서 방출되는 다양한 우주선 피드백이 대질량 은하의 특성에 미치는 영향

초록

FIRE‑3 고해상도 시뮬레이션을 이용해 AGN가 방출하는 우주선(CR)의 주입 효율과 전파 방식을 다양하게 바꿔 보았다. 모든 모델이 별 질량‑암흑질량 관계와 SMBH‑σ 관계를 만족하며 은하를 효과적으로 퀜칭했지만, CGM의 온도·밀도·CR 에너지 밀도는 1 dex 이상 차이를 보였다. 따라서 CGM 관측이 CR‑AGN 피드백의 미세 물리학을 제한하는 핵심 도구가 될 수 있다.

상세 분석

본 연구는 FIRE‑3 프레임워크에 포함된 최신 CR‑MHD 모듈을 활용해 네 개의 질량이 10¹³ M⊙ 수준인 ‘m13’ 은하를 대상으로 네 가지 CR 피드백 설정을 비교하였다. 첫 번째는 고정 파워‑법칙 형태의 확산 계수(‘CD’)이며, 두 번째는 은하 내부 환경에 따라 확산 계수가 변하는 ‘VD’ 모델이다. VD 모델은 Milky Way ISM 관측(VOYAGER, AMS‑02)을 재현하도록 캘리브레이션된 외부 구동(External‑Driving) 스키마를 사용한다. 각 VD 모델에 대해 ϵ_BH,cr=10⁻⁴, 10⁻³, 3×10⁻³ 세 가지 주입 효율을 적용해 총 7개의 시뮬레이션을 수행하였다.

시뮬레이션 결과는 크게 두 축으로 해석할 수 있다. (1) 은하 전체 규모의 규제 – 모든 모델이 z=0에서 M★‑Mhalo, M_BH‑σ, 그리고 300 Myr 평균 SFR 등 관측 스케일링 관계에 들어맞는다. 특히 ϵ_BH,cr이 10⁻³ 수준이면 별 형성률이 0.1 dex 이하로 억제되어 ‘quenched’ 영역에 위치한다. VD 모델은 동일 효율에서도 CD 모델보다 약간 더 강력한 규제를 보이며, 이는 변동 확산 계수가 고밀도 중심부에서 CR 압력을 효율적으로 전달하기 때문이다. (2) CGM 물리량의 민감도 – 반면, 은하 주변 0.1–1 R_200 영역의 가스 온도·밀도·CR 에너지 밀도는 모델마다 1–2 dex 차이를 보인다. 높은 ϵ_BH,cr(VDHiCR)에서는 CR 압력이 외부 가스를 팽창시켜 온도가 10⁶ K 이하로 낮아지고, 중성 원소 이온화 상태가 변한다. 반대로 낮은 효율(VDLoCR)에서는 CR가 거의 무시되므로 전통적인 열·동역학 피드백(기계적·복사)만으로 CGM가 유지된다. 이러한 차이는 CR 확산 계수의 반경 의존성, 특히 외부 구동 모델이 저밀도 외부 가스에서 확산을 억제하고 흐름을 강화하는 메커니즘에 기인한다.

핵심적인 물리적 함의는 다음과 같다. 첫째, AGN‑CR 피드백은 은하 퀜칭에 충분히 강력하지만, CGM 구조를 크게 재구성한다는 점이다. 둘째, CR 전파 파라미터(확산/스트리밍)와 주입 효율이 서로 보완적으로 작용한다. 예를 들어, 낮은 효율이지만 높은 확산 계수(VD)도 충분히 퀜칭을 달성한다. 셋째, 관측 가능한 CGM 지표(OVI, NeVIII, X‑ray 표면 밝기 등)가 CR‑AGN 모델을 구별할 수 있는 강력한 도구가 된다. 따라서 향후 다파장(UV, X‑ray, SZ) 합성 관측을 통해 CGM 전자 밀도·온도·비열을 정밀 측정하면, 현재 시뮬레이션에 남아 있는 ‘미세 물리학’(CR 스캐터링, 파동‑입자 상호작용) 파라미터를 제한할 수 있다.