초고온 목성의 극한 바람

초록

케플러-9b의 일광면을 고해상도 광학 스펙트럼으로 관측해, 위상에 따라 변하는 도플러 이동을 측정하였다. 회전 속도(6.4 km s⁻¹)를 넘어서는 -13.4 km s⁻¹에서 -0.4 km s⁻¹까지의 순풍 이동을 확인했으며, 이는 일광면에 11.7 km s⁻¹에 달하는 초음속 동-야풍이 존재함을 의미한다.

상세 분석

본 연구는 Keck Planet Finder(KPF) 광학 고해상도 분광기를 이용해 초고온 목성(KELT‑9 b)의 위상곡선 스펙트럼을 ϕ = 0.05–0.48 구간에서 연속적으로 수집하였다. 각 노출 150 s로 설정해 행성의 반대편 반사속도(Kp ≈ 245 km s⁻¹)로 인한 스미어링을 최소화했으며, 평균 S/N ≈ 400을 달성했다. 데이터 전처리 단계에서는 레이저 주파수 보정, 블레이즈 보정, molecfit을 통한 수증기·산소 대기 흡수 제거, 그리고 별 스펙트럼을 마스터 템플릿으로 나눠 별 라인을 제거하였다. 이후 고차 다항식으로 연속체를 정규화하고, 강한 대기 흡수 구간을 5σ 이상 클리핑해 전체 파장 1.2 %를 마스크 처리했다.

모델링에는 petitRADTRANS와 Guillot(2010)식 P‑T 프로파일을 사용해 적외선·가시광선 불투과도(κ_IR, γ)와 평형 온도(T_eq)를 자유 파라미터로 설정하였다. 화학 평형은 easychem으로 실시간 계산했으며, Fe, Ti, Ca, Mg, Si 등 금속 원소의 전리와 라인 리스트는 Kurucz(2018)를 기반으로 pyROX를 통해 9000 K까지 확장하였다. 모델 스펙트럼은 R = 10⁶ 원본 라인표를 3배 다운샘플링해 연산 효율을 높인 뒤, KPF 해상도(R ≈ 98 000)와 Gaussian 커널로 컨볼루션했다. 행성‑별 플럭스 비율은 (R_p/R_)² · F_p/F_ 로 계산하고, 최종적으로 150 pixel(≈2–4 Å) Gaussian 스무딩을 빼 연속체를 제거하였다.

교차상관(CCF) 분석은 모델 템플릿을 -300 ~ +300 km s⁻¹ 범위에서 1 km s⁻¹ 간격으로 이동시키며, 관측 불확실성(σ₀)으로 가중된 공분산 행렬을 적용해 수행했다. Kp‑vs‑v_sys 행렬을 구축해 최적 Kp ≈ 245 km s⁻¹와 v_sys ≈ -17.6 km s⁻¹를 도출했으며, 위상에 따라 CCF 피크가 -13.4 km s⁻¹에서 -0.4 km s⁻¹까지 이동하는 것을 확인했다. 이 이동 폭은 행성의 자전 속도(6.4 km s⁻¹)를 크게 초과하므로, 순수 자전이 아닌 동‑야풍이 주된 원인임을 입증한다. 전이된 풍속은 위상별 베이즈 추정으로 ϕ ≈ 0.1에서 11.7 ± 0.6 km s⁻¹, ϕ ≈ 0.45에서는 거의 0 km s⁻¹에 수렴한다.

3D 일반 순환 모델(GCM)과 비교했을 때, 관측된 풍속은 약한 마찰(drag) 파라미터와 일치한다. 이는 전자 전도성 대기와 자기장 결합에 의한 마찰이 상대적으로 약해 열 재분배 효율이 높음을 시사한다. 실제 Spitzer·TESS·CHEOPS 위상곡선에서 추정된 재분배 효율(≈0.3)과도 일관된다. 화학적으로는 Fe I/II, Ti I/II, Ca I/II, Mg I, Si I이 강하게 검출되었으며, 온도 구조는 고도 10⁻³–10⁻⁴ bar에서 온도 역전(>4000 K)이 뚜렷함을 보여준다. 밤면 온도는 제한된 위상 데이터로부터 상부 대기(≤10⁻⁴ bar)에서 2500–3000 K 수준으로 추정되었다.

결론적으로, 고해상도 위상곡선 스펙트로스코피는 초고온 목성의 대기 역학을 직접 측정할 수 있는 강력한 도구이며, KELT‑9 b는 현재까지 관측된 가장 강력한 동‑야풍(≈12 km s⁻¹)을 보유한 행성으로 확인되었다. 이러한 관측은 극한 방사선 환경에서의 대기 순환, 마찰, 그리고 전자기 상호작용을 이해하는 데 중요한 벤치마크를 제공한다.