A $z \sim$ 6.2 Quasar on the Local M$_{\rm BH}$-$\sigma_{\rm \ast}$ Relation Quenching Its Host Galaxy from the Aether Survey

기술 요약: 로컬 $M_{BH}$–$\sigma_*$ 관계를 따르는 $z \sim 6.2$ 퀘이사: Aether 서베이를 통한 호스트 은하의 퀜칭(Quenching) 연구

문제 및 동기
초거대질량 블랙홀(SMBH)과 그 호스트 은하의 공동 진화는 국부 우주에서 블랙홀 질량($M_{BH}$)과 호스트 은하의 특성, 구체적으로 별의 속도 분산($\sigma_*$) 및 별의 질량($M_*$) 사이의 긴밀한 상관관계로 입증된다. 그러나 이러한 상관관계가 형성되는 시점과 메커니즘은, 특히 초기 우주($z > 6$)에서 여전히 미해결 과제로 남아 있다. 빅뱅 후 첫 10억 년 이내에 $\sim 10^9 - 10^{10} M_\odot$ 규모의 SMBH를 가진 퀘이사들이 존재하지만, 퀘이사의 압도적인 빛 때문에 희미한 호스트 은하의 $\sigma_*$를 측정하는 것은 매우 어렵다. 이전 연구들은 가스 운동학(예: [C II])이나 저해상도 분광법을 사용하였으며, 이는 종종 블랙홀이 국부 관계에 비해 과하게 무겁다는 결과를 내놓거나 상충되는 결과를 보여왔다. 그러나 최근 JWST 연구들은 일부 고적색편이 천체들이 이미 국부 스케일링 관계와 일치할 수 있음을 시사하고 있다. 나아가, 초기 은하에서 별 형성 과정을 빠르게 멈추게 하는 퀜칭(quenching)의 동력원인 퀘이사 피드백의 메커니즘을 규명하기 위해서는 직접적인 공간 해상도를 가진 증거가 필요하다.

방법론
본 연구는 Aether 서베이(Program ID 5645)의 일환으로 얻은 $z \simeq 6.18$ 퀘이사 J1512+4422에 대한 새로운 JWST/NIRSpec 적분 필드 단위(IFU) 관측 데이터를 제시한다. 관측에는 G395H/F290LP 격자/필터 구성이 사용되었으며, $\sim 2.87–5.14 \, \mu\text{m}$의 파장 범위와 $\lambda/\Delta\lambda \simeq 2700$ ($\sim 130 \, \text{km s}^{-1}$)의 분해능을 제공한다.

분석은 다음과 같은 주요 단계로 진행되었다:

1. 데이터 감소(Data Reduction): IFU 데이터는 노이즈 제거, 플럭스 보존 재투영, 이상치 제거를 위한 커스텀 스크립트를 포함하여 JWST 과학 교정 파이프라인(v1.14.0)을 사용하여 처리되었다.
2. 핵 영역 스펙트럼 피팅: 핵 영역 스펙트럼은 $0.1''$ 구경에서 추출되었으며 BADASS를 사용하여 모델링되었다. 모델에는 퀘이사 멱함수 연속 스펙트럼, Fe II 템플릿, 항성 연속 스펙트럼 및 방출선이 포함되었다. 결정적으로, 별의 연속 스펙트럼과의 퇴행(degeneracy)을 깨기 위해 IFU 데이터의 GALFIT 이미지 분해를 사용하여 퀘이사 연속 스펙트럼 기울기를 독립적으로 제약하였다.
3. 별의 운동학: 별의 속도 분산($\sigma_*$)은 BADASS 내의 pPXF 코드를 사용하여 발머 계열(Balmer series) 흡수선을 피팅함으로써 도출되었다.
4. 유출(Outflow) 분석: 퀘이사 점확산 함수(PSF)를 q3dfit 도구로 차감하여 확장된 선 방출을 분리하였다. 결과로 얻은 공간 해상도 스펙트럼은 속도($v_{50}$)와 선 폭($\sigma$)을 기준으로 가스 유출과 핵 유출을 식별하기 위해 가우시안 성분으로 분해되었다.
5. 에너지 계산: 전자 밀도 범위 $100–1000 \, \text{cm}^{-3}$를 가정하여 핵 영역(미해상) 및 확장 영역(해상) 성분에 대한 질량, 운동량 및 운동 에너지 유출율을 계산하였다.

주요 결과

• 블랙홀 질량 및 별의 속도 분산: 본 연구는 발머 흡수선의 확장을 통해 도출된 $\sigma_* = 288 \pm 60 \, \text{km s}^{-1}$의 별의 속도 분산을 측정하였다. 이 값은 저해상도 고정 슬릿(FS) 분광법으로 도출된 이전의 상한값($<190 \, \text{km s}^{-1}$)보다 유의미하게 높다. 블랙홀 질량은 넓은 H$\alpha$ 선을 기준으로 $M_{BH} \simeq 8.9 \times 10^8 M_\odot$, 넓은 H$\beta$ 선을 기준으로 $M_{BH} \simeq 6.2 \times 10^8 M_\odot$로 추정된다.
• 스케일링 관계에서의 위치: 이러한 측정치를 바탕으로, J1512+4422는 빅뱅 후 불과 $\sim 900$ Myr 만에 국부 $M_{BH}$–$\sigma_*$ 관계(Kormendy & Ho 2013) 위에 놓여 있다. 그러나 블랙홀은 국부 $M_{BH}$–$M_*$ 관계에 비해서는 과하게 무거운 것으로 보이나, 고적색편이에서의 별 질량 측정에 대한 불확실성이 언급되었다.
• 유출 탐지: 핵 영역에서 시작되어 투영 거리 $\sim 3.2$ kpc까지 확장되는 빠른 은하 규모의 유출이 탐지되었다. 핵 유출의 속도는 $\sim 478 \, \text{km s}^{-1}$이며, 확장 유출의 중앙값 속도는 $\sim 352 \, \text{km s}^{-1}$이다. 운동학적으로 정적인 가스의 속도 분산은 단 $\sim 86 \, \text{km s}^{-1}$였으며, 이는 가스 운동학이 이 천체에서 별의 운동학을 대변하는 신뢰할 만한 프록시가 아님을 확인시켜 준다.
• 퀜칭 메커니즘: 유출의 역학적 시간 척도는 $\sim 9$ Myr이며, 이는 최근의 별 형성 억제(지난 10 Myr 동안의 비특이적 별 형성율(sSFR) $< 0.2 \, \text{Gyr}^{-1}$)와 일치한다. 총 질량 유출율($\sim 92.6^{+92.6}_{-74.1} M_\odot \text{yr}^{-1}$)은 현재의 별 형성율($0.9^{+3.8}_{-0.8}$ 또는 $4.3^{+5.8}_{-3.7} M_\odot \text{yr}^{-1}$)을 초과한다. 운동 에너지 유출율은 퀘이사의 볼로미터 광도의 $\sim 0.6\%$이며, 이는 별 형성을 억제하는 데 필요한 음의 피드백(negative feedback)의 이론적 임계값을 충족한다.

의의 및 주장
본 논문은 J1512+4422가 초거대질량 블랙홀과 호스트 은하의 공동 진화, 구체적으로 $M_{BH}$–$\sigma_*$ 관계의 확립이 우주의 첫 10억 년 이내에 일어날 수 있다는 직접적인 증거를 제공한다고 주장한다. 저자들은 퀘이사 구동 피드백이 이러한 천체들을 국부 스케일링 관계에 위치시키고 호스트 은하를 퀜칭하는 핵심 메커니즘이라고 주장한다.

호스트 은하의 특성($\sigma_*$, 별 질량, 크기)은 $z \gtrsim 3$의 정적(quiescent) 및 포스트-스타버스트(post-starburst) 은하(예: GS-9209)의 특성과 유사하게 나타났다. 이는 초기 퀘사르가 활동을 멈춘 후, 높은 적색편이에서 관측되는 거대 정적 은하로 진화한다는 직접적인 진화적 연결 고리를 시사한다. 본 연구는 J1512+4422와 같은 천체의 경우, SMBH의 성장과 호스트 은하의 성장이 효과적인 피드백에 의해 조절되며 초기부터 긴밀하게 결합되어 있음을 강조한다. 또한, 고적색편이 퀘사 true에서 별의 운동학을 정확하게 측정하기 위해서는 고해상도 IFU 데이터가 필수적이며, 저해상도나 가스 기반의 프록시는 오해의 소지가 있는 낮은 속도 분산을 산출할 수 있음을 강조한다.