A comprehensive Rossiter-Mclaughlin Modelling Framework in TLCM: Application to HD 2685 $=$ TOI-135 system

기술 요약: TLCM 내의 포괄적인 로스피터-맥클러린 모델링 프레임워크: HD 2685 = TOI-135 시스템에 대한 적용

문제 제기
로스피터-맥클러린(Rossiter-McLaughlin, RM) 효과는 통과하는 외계행성의 하늘 투영 스핀-궤도 각도($\lambda$)를 측정하는 주요 방법으로, 행성 이동, 역학적 흥분 및 조석 진화에 대한 중요한 통찰을 제공한다. 경험적 추세는 항성의 유효 온도($T_{\text{eff}}$)와 이각(obliquity) 사이의 상관관계(뜨거운 별이 더 넓은 이각 분포를 보임)를 시사하지만, 최근의 인구 통계 연구는 항성 다중성, 행성 대 항성 질량비, 그리고 원시 원반 조건과 같은 요인들이 중요한 역할을 한다는 것을 보여준다. 이러한 인구 수준의 진단 도구를 정교화하기 위해서는 eccentricity가 있는 궤도, 다양한 radial velocity(RV) 오프셋, 그리고 다양한 림 암흑(limb-darkening) 법칙을 처리할 수 있는 견고하고 유연하며 검증된 모델링 프레임워크가 필요하다. 나아가, 진화한 뜨거운 모항성을 도는 핫 주피터인 HD 2685와 같은 특정 시스템은 궤도 구조와 이각을 제약하기 위해 확장된 광도 및 분광 데이터셋을 통한 업데이트된 분석을 요구한다.

방법론
저자들은 TLCM(Transit and Light Curve Modeller) 코드 내에서 로스피터-맥클러린 효과의 업데이트된 구현을 제시한다. 주요 방법론적 발전 사항은 다음과 같다:

1. RM 이상 현상의 수치 적분: 특정 림 암흑 법칙(예: 이차식)에 국한되었던 기존의 해석적 접근 방식과 달리, 새로운 프레임워크는 가려진 항성 디스크 영역에 대한 직접적인 수치 적분을 통해 RM 이상 현상을 계산한다. 이를 통해 모델은 새로운 해석적 적분을 유도하지 않고도 7가지의 서로 다른 림 암흑 법칙(선형, power-2, 이차식, 로그, 제곱근, Sing-3, Claret 4-매개변수)을 수용할 수 있다.
2. 대류 블루시프트(Convective Blueshift, CB) 모델링: 코드는 $\lambda$ 측정에서 계통 오차를 유발할 수 있는 현상인 대류 블루시프트에 대한 해석적 설명을 포함한다. 모델은 Liebing 등(2021)으로부터 유도된 온도 의존 공식을 사용하여 항성 표면의 순 블루시프트를 설명하며, 고온에서는 지수적 컷오프를 적용하고 저온에서는 상수 값을 적용한다.
3. 결합 모델링 능력: RM 모듈은 TLCM의 기존 결합 광도 및 RV 프레임워크와 결합된다. 이는 다음을 동시에 모델링한다:
   • 통과(transit), 식(occultation) 및 위상 곡선(beaming, reflection, ellipsoidal 효과 포함).
   • 이심 궤도(eccentric orbits) 및 다중 케플러 신호.
   • 기기 제로 포인트 오프셋 및 RV 지터(jitter).
   • 긴 고해상도 광도 곡선을 위해 밀집된 가우시안 프로세스의 효율적인 대안을 제공하는 웨이브릿 기반 모델을 사용한 상관 잡음(correlated noise).
4. 검증: 업데이트된 코드는 정렬된 궤도, 중간 이각, 극성 및 역행 궤도, 그리고 이심 궤도 및 다중 RV 오프셋을 가진 시스템을 포함한 9개의 알려진 시스템(WASP-15, HAT-P-1, HAT-P-7, HAT-P-11, HAT-P-32 등 포함)을 통해 검증되었다.

HD 2685 = TOI-135에 대한 적용
본 프레임워크는 약 $1.15 M_{\text{Jup}}$의 핫 주피터를 약 4.1일 궤도로 거느린 단일한 진화된 F형 항성($T_{\text{eff}} \approx 6801$ K)인 HD 2685 시스템에 적용되었다.

• 데이터: 본 연구는 TESS 광도 데이터를 7개 섹터(1, 27, 28, 67, 68, 94, 95)로부터 활용하여, 원래의 발견(섹터 1만 사용) 때보다 데이터셋을 크게 확장하였다. 통과 중 및 통과 외 구간의 새로운 HARPS 방사 속도 측정이 수행되었으며, 이는 CORALIE, CHIRON, FEROS의 아카이브 데이터와 결합되었다.
• 분석: 광도 곡선과 RV 데이터를 이용한 결합 피팅이 DE-MCMC 방법을 사용하여 수행되었다. 분석은 투영된 항성 회전 속도($V \sin I_*$)에 대한 사전 확률(prior)이 유도된 이각($\lambda$)의 민감도를 테스트하였다.

주요 결과

1. 시스템 파라미터: 결합 분석을 통해 시스템 파라미터를 정밀하게 도출하였으며, $a/R_* \approx 7.6$의 척도 반경, 약 $1.44 R_{\text{Jup}}$의 행성 반지름, 약 $1.57 R_\odot$의 항성 반지름을 얻었다. 이심률은 0에 가까운 것으로 나타났으며($e \approx 0.004$), 이는 이전의 추정치와 약간 다르지만 원형 궤도임을 시사한다.
2. 로스피터-맥클러린 탐지: 명확한 RM 신호가 검출되었다. 고정된 분광 사전 확률($V \sin I_* = 15.4 \pm 0.2$ km/s)에 대해, 하늘 투영 이각은 $\lambda = 55.6^{+10.9}_{-11.9}$ 도로 측정되었다.
3. 강건성: $V \sin I_*$에 대한 더 넓은 사전 확률(uniform prior $U(5, 25)$ km/s 포함)을 사용한 테스트 결과, $\lambda$의 공식 불확실성은 증가했지만($+13^\circ, -24^\circ$까지), 중심값은 $56^\circ$ 부근에서 안정적으로 유지되었다. 이는 이 시스템이 중간 정도의 투영 이각을 가지고 있음을 확인시켜 준다.
4. 지터 보정: 사용된 서로 다른 기기들에 대해 RV 지터 값을 재도출하였으며, HARPS의 경우 46 m/s, CORALIE는 48 m/s, CHIRON은 69 m/s, FEROS는 23.5 m/s로 나타났다.

의의 및 주장
본 논문은 두 가지 주요 영역에서 의의를 갖는다:

1. 방법론적 발전: 업데이트된 TLCM은 특정 림 암흑 법칙의 해석적 제한에 구애받지 않는, RM 모델링을 위한 유연하고 수치적으로 강건한 프레임워크를 제공한다. 다양한 시스템에 걸친 성공적인 검증은 그 신뢰성을 입증한다.
2. 인구 통계 진단: HD 2685b의 측정값은 뜨겁고 진화한 항성들이 중간 이각을 보이는 $T_{\text{eff}} - |\lambda|$ 다이어그램의 희소한 샘플 영역에 데이터 포인트를 추가한다. 이 결과는 뜨거운 별들이 차가운 별들에 비해 더 넓은 이각 분포를 보인다는 경험적 추세와 일치하며, 이는 복사 외포를 가진 별들에서 조석 재정렬의 효율성이 낮기 때문인 것으로 보인다. 그러나 저자들은 이 시스템의 낮은 행성 대 항성 질량비($\sim 0.0008$)가 정렬 상태를 유지하는 경향이 있는 고질량비 시스템과 대조적으로, 불일치(misalignment)가 발생하는 저질량비 시스템에 대한 이론적 기대(Rusznak et al. 2025)와 일치한다는 점을 언급한다.

저자들은 본 연구가 기존의 경험적 추세를 뒷받침하고 있지만, 서로 다른 항성 및 행성 파라미터 공간 전반에 걸친 스핀-궤도 각도 분포의 산포를 완전히 설명하기 위해서는 추가적인 관측 및 이론적 연구가 필요하다고 결론짓는다.