The CRIMSON survey I: super-stellar SiO in the directly imaged companion TWA 5 B from high-resolution M-band spectroscopy

기술 요약: CRIMSON 서베이 I: 직접 영상화된 동반성 TWA 5 B의 초거성(Super-stellar) SiO

문제 및 동기
거대 행성의 형성사를 이해하려면 현재의 대기 조성과 원시 행성계 원반의 화학적 분할(chemical partitioning) 사이의 연결 고리를 찾아야 한다. 휘발성 종(예: CO, H$_2$O)은 C/O 비를 통해 형성 위치를 추적하는 데 널리 사용되어 왔으나, 이 측정값들은 원반 내의 복잡한 열적 및 화학적 처리 과정으로 인한 퇴보(degeneracy) 문제를 겪는다. 이를 보완하는 접근법은 난성분(refractory elements, 예: Si, Fe)을 측정하는 것인데, 이들은 원반의 대부분의 과정 동안 고체 상태로 남아 있다. 난성분과 휘발성 원소의 비율은 고체와 가스의 강착(accretion)을 규명할 수 있다. 그러나 직접 영상화된 행성에서 난성분 종을 검출하는 것은 매우 어렵다. 초고온 목성(ultra-hot Jupiters)의 경우, 가시광선 영역에서 원자 또는 이온 형태로 난성분이 자주 검출되지만, 이러한 종들은 매우 높은 온도($T_{eq} \gtrsim 2000$ K)에서만 제한적으로 나타난다. 더 차가운 직접 영상화된 동반성들의 경우, 난성분은 구름(예: 마그네슘 규산염)으로 응축되어 기상(gas phase)을 결핍시킬 것으로 예상된다. 초기 난성분 함량과 구름의 특성을 추적하기 위해서는 고해상도 분광법을 통해 분자 형태의 난성분에 접근해야 한다. M-밴드(3.5–5.2 $\mu$m)는 실리케이트 구름이 없는 뜨거운 대기에서 실리콘의 주요 운반체인 기체 상태의 일산화규소(SiO)를 검출할 수 있는 창을 제공하며, 이는 대기의 난성분 인벤토리와 구름 응축 상태를 직접적으로 조사할 수 있는 도구가 된다.

방법론
본 연구는 VLT의 CRIRES+ 기기를 사용하여 직접 영상화된 동반성 TWA 5 B에 대한 고해상도($R \sim 92,000$) M-밴드 분광 측정을 수행한 CRIMSON 서베이의 첫 번째 결과를 제시한다. TWA 5 B는 젊고($\sim 10$ Myr) 뜨거운($T_{eq} \approx 2400$ K) 아성간(sub-stellar) 동반성으로, 유의미한 구름 응축이 일어나기 전 기상 SiO를 검출하기에 이상적인 대상이다.

• 관측: 2025년 3월 19일에 총 2.2시간의 노출 시간을 갖는 두 개의 관측 블록이 수행되었다. 관측에는 0.2" 슬릿과 M-밴드 관측에 내재된 높은 열적 배경을 제거하기 위한 ABBA 노딩 패턴이 사용되었다.
• 데이터 감소 및 추출: 본 연구의 핵심적인 방법론적 기여는 커스텀 스펙트럼 추출 루틴의 개발이다. 저자들은 곡선 형태의 스펙트럼 트레이스에 대한 검출기의 이산적 공간 샘플링(0.056"/pixel)으로 인해 발생하는 CRIRES+ 2D 스펙트럼의 주기적인 계통 오차를 확인하였다. 이를 완화하기 위해, 스펙트럼을 추출하기 전 PSF를 재구성하기 위해 스플라인 보간법을 사용하여 100배 오버샘플링하는 오버샘플링 추출 루틴을 구현하였다. 이를 통해 계통 오차를 플럭스 수준의 $\le 10^{-3}$ 이하로 줄였다.
• 후처리: 데이터는 molecfit을 사용하여 지구 대기 흡수를 보정하였고, 파장 영역에서 적용된 주성분 분석(PCA)을 통해 항성 오염을 제거하였다. 결정적으로, 동반성의 스펙트럼 특징이 침식되는 것을 방지하기 위해 동반성 신호가 포함된 공간 행을 제외하고 PCA 고유 벡터를 계산하였다.
• 분석: 분자 검출은 대기 모델 템플릿과의 교차 상관(cross-correlation)을 통해 이루어졌다. 대기 특성(화학적 풍부도, 역학, 온도-압력(TP) 프로파일)은 베이지안 중첩 샘플링(Bayesian nested sampling) 접근법을 사용하는 HyDRA 프레임워크를 통해 역산(retrieval)되었다. 역산에는 H$_2$O, CO, SiO에 대한 "자유 화학(free chemistry)" 방식이 사용되었으며, 매개변수화된 구름 덱 모델이 포함되었다.

주요 결과

1. 분자 검출: 본 연구는 직접 영상화된 동반성에서 기체 상태의 SiO를 최초로 확정 검출하였으며, 신호 대 잡음비(S/N)는 7.5이다. 휘발성 종의 강력한 검출 또한 확인되었다: CO (S/N = 9.1) 및 H$_2$O (S/N = 18.8). 결합 모델은 S/N = 21.6로 검출되었다.
2. 대기 풍부도: 역산된 SiO의 부피 혼합비(VMR)는 $\log(\text{SiO}) = -3.56^{+0.42}_{-0.32}$이다. 이 풍부도는 매우 높으며, 이는 대기의 실리콘이 거의 전적으로 기상 SiO 내에 존재하며 응축물로의 결핍이 미미함을 의미한다.
3. 구름 특성: 높은 SiO 풍부도와 함께 역산된 낮은 구름 불투명도($\log(\kappa) \approx -7.0$) 및 높은 압력 쪽으로 밀려난 불확실한 구름 덱 압력은, 조사된 압력 범위($-2 \lesssim \log(P [\text{bar}]) \lesssim 0$)에서 상당한 마그네슘-실리케이트 구름 응축이 존재하지 않음을 나타낸다. 이는 해당 천체의 온도가 이러한 구름들의 응축 임계값보다 높다는 점과 일치한다.
4. 원소 비율 및 형성:
   • Si/H: 대기는 항성보다 높은 실리콘 과잉($[\text{Si/H}]_\star = 1.41^{+0.42}_{-0.32}$)을 보인다. 이는 항성 대비 약 25배의 농축을 나타낸다.
   • 휘발성-대-난성분 비율: C/O 비는 항성과 유사하며($\sim 0.26$), O/Si 및 C/Si는 약간 항성보다 낮다.
   • 형성 시나리오: 이러한 비율은 두 가지 형성 경로, 즉 CO 설선(snowline) 내부 또는 외부에서의 핵 성장(core-accretion) 모델, 또는 중력 불안정성(gravitational instability)에 의한 붕괴 후 상당한 고체 농축 모델과 일치한다. 그러나 TWA 5 B의 높은 질량($25 \pm 5 M_J$)과 넓은 궤도($\sim 86$ AU)는 핵 성장을 해당 시스템의 추정된 원반 질량과 조화시키기 어렵게 만든다. 따라서 붕괴 후 약 $20 M_\oplus$의 실리케이트가 풍부한 고체를 강착하는 중력 불안정성 모델이 실리콘 농축을 설명하기 위한 유효한 대안으로 제시된다.
5. 역학: 투영 회전 속도는 $v \sin(i) = 14.19^{+0.88}_{-0.79}$ km/s로 측정되었다.

의의 및 주장
본 논문은 TWA 5 B에서 기체 상태의 SiO를 검출한 것이 뜨거운 거대 행성의 구름 특성을 위한 독특한 진단 도구를 제공한다고 주장한다. SiO는 뜨거운 대기에서 실리콘의 주요 운반체이며 오직 응축에 의해서만 결핍되기 때문에, 그 풍부도는 마그네슘-실리케이트 구름의 존재 여부를 판단하는 직접적인 "리트머스 시험지" 역할을 한다. 저자들은 이 측정이 TWA 5 B에서 그러한 구름의 부재를 확인해주며, 이를 통해 대기의 난성분 함량을 순수하게 측정할 수 있게 한다고 논한다.

나아가, 본 연구는 고해상도 M-밴드 분광법이 분자 형태의 난성분 화학에 접근함으로써 휘발성 단독 연구에서 발생하는 퇴보 문제를 해결할 수 있음을 입증한다. 저자들은 직접 영상화된 행성 및 갈색 왜성 집단에 대한 SiO 풍부도 조사가 형성되는 주요 구름 종(예: 포스테라이트 vs 엔스타타이트)과 그 응축 온도를 결정할 수 있으며, 이는 구름 물리학을 규명하기 위한 우주 기반 중적외선 관측의 지상 기반 대안을 제공할 것이라고 상정한다. 마지막으로, 본 연구는 미래의 ELT용 METIS와 같은 기기에서 직면할 수 있는 검출기 샘플링 계통 오차를 극복하기 위해 커스텀 추출 기법이 필수적임을 강조한다.