HD 176986 다중행성계에 새롭게 발견된 초지구 HD 176986 d

초록

HARPS와 HARPS‑N 장기 관측을 YARARA와 다차원 가우시안 프로세스 모델링으로 정밀 분석한 결과, 기존에 알려진 두 초지구 b와 c 외에 61.4일 주기의 새로운 초지구 d(최소 질량 ≈ 6.8 M⊕)를 발견하였다. 또한 별의 자기 사이클(≈ 6.7 yr)과 회전 주기(≈ 36 d)를 정확히 측정하였다.

상세 분석

본 연구는 K2.5 V형 별 HD 176986에 대한 757개의 고해상도 스펙트럼(149 H03, 246 H15, 362 HN)을 활용하였다. HARPS는 2015년 광섬유 교체 전후로 두 개의 독립 데이터셋(H03, H15)으로 구분했으며, HARPS‑N은 별도의 HN 데이터셋을 제공한다. 기존 DRS 기반 CCF RV 외에도 YARARA 파이프라인을 적용해 시스템atics(우주선, 텔레리크, 검출기 결함 등)를 정밀 보정하고, 라인‑바이‑라인(LBL) 기법으로 각 스펙트럼 라인의 RV를 개별 추출하였다. 활동 보정 전후의 RMS 차이는 H03 → 3.82 → 3.68 m s⁻¹, H15 → 3.67 → 3.14 m s⁻¹, HN → 3.51 → 3.23 m s⁻¹로, 평균 오차도 0.6 → 0.4 m s⁻¹ 수준으로 크게 향상되었다.

활동 지표(FWHM, BIS, Ca II H&K, Hα 등)와 TESS·ASAS‑SN 광도 데이터를 동시 분석해 별의 회전 주기 36.05 ± 0.7 d와 장기 자기 사이클 2432 ± 60 d를 도출하였다. RV와 활동 지표를 다차원 가우시안 프로세스(GP) 모델에 공동 적용함으로써 활동 신호와 행성 신호를 효과적으로 분리했다. 모델 비교는 Nested Sampling(Dynesty) 기반 베이지안 증거를 사용했으며, 3행성 모델이 2행성 모델 대비 Δln Z ≈ +12로 강력히 선호됨을 확인했다.

새로운 행성 d는 P = 61.376 ± 0.05 d, K ≈ 1.9 m s⁻¹, 최소 질량 M sin i = 6.76 ± 0.9 M⊕를 보이며, 기존 b(6.49 d, 5.36 M⊕)와 c(16.81 d, 9.75 M⊕)와 동적 안정성을 유지한다. d의 공전 반지름은 약 0.25 AU(≈ 0.23 AU)로, 보수적·낙관적 정의의 거주 가능 영역(HZ) 바깥에 위치하지만, 향후 장기 관측으로 HZ 내 저질량 행성을 탐색할 가능성을 열어준다.

활동 보정 전후의 RV 시계열을 비교하면, YARARA‑활동 보정이 없는 경우에도 3행성 신호가 검출되지만, 보정 후에는 신호의 신뢰도와 파라미터 불확도가 현저히 감소한다. 이는 라인‑레벨 보정이 활동 잡음을 효과적으로 억제함을 시사한다. 또한, 별의 금속도