JWST MIRI가 밝힌 은하핵 실리케이트 방출과 토러스 구조

초록

JWST/MIRI MRS로 관측한 근거리 QSO2 J1010의 핵심 0.5–1.1 kpc 영역에서 9.7, 18, 23 µm 실리케이트 특징이 강하게 방출됨을 확인했다. 베이지안 툴 BayesClumpy2를 이용해 CLUMPY와 CAT3D‑WIND 토러스 모델을 적합했으며, CAT3D‑WIND가 marginal likelihood와 잔차 면에서 우수했다. 두 모델 모두 중간 정도의 토러스 커버링 팩터(C_T≈0.5–0.7)와 낮은 시야각(i≈13–50°)을 요구한다. 네 개 QSO2가 Eddington 비율‑칼럼밀도 도표의 “blow‑out” 영역에 위치해 방사압에 의한 가스·먼지 제거가 진행 중임을 시사한다. 이는 JWST로 관측된 Seyfert 2와 대비돼, AGN 광도와 Eddington 비율이 높을수록 토러스 커버링 팩터가 감소한다는 시나리오를 뒷받침한다.

상세 분석

본 연구는 JWST/MIRI MRS가 제공하는 0.3–0.8″(≈0.5–1.1 kpc) 고해상도 중적외선 스펙트럼을 활용해, 전통적으로 광학적으로 type‑2 로 분류된 QSO2 J1010의 핵심 영역에서 실리케이트 9.7 µm, 18 µm, 23 µm 방출을 명확히 검출한 점이 가장 큰 혁신이다. 이전 Spitzer/IRS는 수십 kpc 규모의 빔을 사용해 호스트 은하의 먼지와 혼합된 신호를 제공했으나, JWST는 핵심 구역만을 분리함으로써 직접적으로 가열된 뜨거운 먼지(수백 K~1000 K)의 방출을 관측했다.

스펙트럼 전처리 단계에서는 PAH와 강한 라인들을 BayesClumpy2에 내장된 Donnan et al. (2024) 모델로 피팅 후 제거했으며, 남은 연속 스펙트럼이 토러스 모델링에 사용되었다. 두 토러스 모델은 구조적 차이를 가진다. CLUMPY는 구형 구름이 토러스 형태로 분포하는 순수 클러디 모델이며, 주요 파라미터는 구름 수(N₀), 구름 광학 깊이(τ_V), 토러스 각도(σ), 반경 분포(q), 그리고 시야각(i)이다. 반면 CAT3D‑WIND은 평면 디스크와 방사압에 의해 가속된 먼지 바람(wind) 성분을 결합한 하이브리드 모델로, 바람의 개시 각도(θ_w), 바람 내 구름 수(N_w), 그리고 바람과 디스크의 상대적 커버링 팩터를 추가한다.

베이지안 적합 결과, 두 모델 모두 C_T≈0.45–0.66(불확실성 포함)와 i≈13°–50°를 요구한다. 특히 CAT3D‑WIND은 marginal likelihood가 CLUMPY보다 약 1.8배 높으며, 9.7 µm 실리케이트 피크와 23 µm 결정성 실리케이트 밴드의 미세 구조를 더 정확히 재현한다. 이는 바람 성분이 직접적인 고온 구름을 노출시켜 실리케이트 방출을 강화시키는 물리적 메커니즘과 일치한다.

또한, J1010을 포함한 4개의 QSO2가 Eddington 비율(λ_Edd≈0.1–1)과 CO 기반 열량밀도(N_CO≈10^22.3 cm⁻²) 조합으로 “blow‑out” 혹은 “forbidden” 영역에 위치한다는 점은 방사압 피드백이 토러스 구조를 얇게 만들고, 커버링 팩터를 감소시킨다는 가설을 실증한다. 반대로 JWST로 관측된 Seyfert 2(λ_Edd≈10⁻³–10⁻², N_H≈10^22.5–10^24 cm⁻²)는 전형적인 실리케이트 흡수를 보이며, 토러스가 보다 두껍고 시야각이 큰 상태임을 시사한다.

이러한 차이는 AGN 광도와 Eddington 비율이 토러스 물리와 구조를 조절한다는 이론적 프레임워크(예: Ricci et al. 2017, Hönig & Kishimoto 2017)와 일맥상통한다. 고광도 QSO2는 강한 방사압이 먼지 구름을 휘발시키거나 바람으로 끌어올려, 토러스의 평균 시야각을 감소시키고, 실리케이트 방출을 가능하게 만든다. 반면 저광도 Seyfert 2는 방사압이 충분히 강하지 않아 토러스가 유지되고, 실리케이트는 주로 그림자 구름에서 흡수된다.

연구의 한계로는 모델 파라미터 간 강한 상관관계(예: C_T와 i)와 제한된 파장 범위(5–25 µm)에서만 적합을 수행했기 때문에, 장파장(>30 µm) FIR 데이터와 X‑ray 흡수 컬럼을 동시에 고려한 다중 파장 베이지안 분석이 필요하다. 또한, 바람 성분의 물리적 파라미터(속도, 밀도 분포)는 현재 모델에 고정돼 있어, 실제 관측된 분자/이온화 가스 흐름과의 정량적 비교가 부족하다. 향후 ALMA와 JWST NIRSpec을 결합한 전방위적 관측이 이러한 불확실성을 해소할 것이다.