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

技術要約：TLCMにおける包括的なロス・シドニー・マクローリン・モデリング・フレームワーク：HD 2685 = TOI-135系への適用

問題提起
ロス・シドニー・マクローリン（RM）効果は、トランジットする系外惑星の空投影スピン・軌道角（$\lambda$）を測定するための主要な手法であり、惑星の移動、動的励起、および潮汐進化に関する重要な知見を提供する。経験的な傾向として、恒星の有効温度（$T_{\text{eff}}$）と傾斜角の間に相関関係があること（高温の恒星はより広い傾斜角分布を示すこと）が示唆されているが、近年の人口統計学的研究は、恒星の多重度、惑星対恒星の質量比、および原始円盤の条件などの要因も重要な役割を果たしていることを示している。これらの集団レベルの診断を精緻化するためには、離心軌道、複数の視線速度（RV）オフセット、および様々な肢端減光則を含む多様な系構成を扱うことができる、堅牢で柔軟かつ検証済みのモデリング・フレームワークが必要である。さらに、HD 2685のような、進化した高温のホスト星を回るホット・ジュピターといった特定の系については、軌道構造と傾斜角を制約するために、拡張された光度および分光データを用いた最新の解析が求められている。

手法
著者らは、Transit and Light Curve Modeller (TLCM) コード内におけるロス・シドニー・マクローリン効果の更新された実装を提示している。主な手法上の進展は以下の通りである：

1. RMアノマリーの数値積分： 特定の肢端減光則（例：二次式）に限定されていた従来の解析的手法とは異なり、新しいフレームワークは、遮蔽された恒星円盤上の部分に対して直接数値積分を行うことで、RMアノマリーを計算する。これにより、モデルは新しい解析的積分を導出することなく、7種類の異なる肢端減光則（線形、2乗、二次式、対数、平方根、Sing-3、およびClaretの4パラメータ）に対応することが可能となる。
2. 対流ブルシフト（CB）モデリング： コードには、$\lambda$ の測定に系統誤差を引き起こす可能性のある現象である、対流ブルシフトの解析的な記述が組み込まれている。モデルは、Liebing et al. (2021) に基づく温度依存の公式を用いて、恒星表面の純粋なブルーシフトを考慮し、高温に対しては指数関数的なカットオフを、低温に対しては一定値を適用する。
3. 結合モデリング機能： RMモジュールは、TLCMの既存の結合光度およびRVフレームワークと結合されている。これは以下の事項を同時にモデリングする：
   • トランジット、掩蔽、および位相曲線（ビーミング、反射、および楕円体効果を含む）。
   • 離心軌道および複数のケプラー信号。
   • 計器のゼロ点オフセットおよびRVジッター。
   • ウェーブレットベースのモデルを用いた相関ノイズ。これは、長時間の高頻度ライトカーブに対して、高密度なガウス過程に代わる計算効率の高い選択肢を提供する。
4. 検証： 更新されたコードは、整列、中間、極、および逆行軌道、ならびに離心軌道や複数のRVオフセットを持つ系を含む、9つの既知の系（WASP-15、HAT-P-1、HAT-P-7、HAT-P-11、HAT-P-32など）を用いて検証された。

HD 2685 = TOI-135への適用
本フレームワークは、単一の進化したF型星（$T_{\text{eff}} \approx 6801$ K）をホストし、約4.1日の軌道を持つ $\sim 1.15 M_{\text{Jup}}$ のホット・ジュピターであるHD 2685系に適用された。

• データ： 本研究では、発見時（セクター1のみ）から大幅に拡張されたデータセットとして、7つのセクター（1, 27, 28, 67, 68, 94, 95）のTESS光度データを活用した。新たに、トランジット中およびトランジット外のHARPS視線速度測定が得られ、これにCORALIE、CHIRON、およびFEROSのアーカイブデータが組み合わされた。
• 解析： ライトカーブとRVデータの結合フィットが、マルコフ連鎖モンテカルロ（DE-MCMC）法を用いて行われた。解析では、投影された恒星回転速度（$V \sin I_*$）の事前分布に対する、導出された傾斜角の感度をテストした。

主な結果

1. 系パラメータ： 結合解析により系パラメータが精緻化され、$a/R_* \approx 7.6$ のスケールされた軌道長半径、$\approx 1.44 R_{\text{Jup}}$ の惑星半径、および $\approx 1.57 R_\odot$ の恒星半径が得られた。離心率はゼロと矛盾しない（$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_*$ の事前分布（$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 であることを見出した。

意義および主張
本論文は、主に2つの領域において意義を主張している：

1. 手法上の進展： 更新されたTLCMは、特定の肢端減光則の解析的な制限を受けない、RMモデリングのための柔軟で数値的に堅牢なフレームワークを提供する。多様な9つの系における成功した検証は、その信頼性を実証している。
2. 集団診断： HD 2685bの測定は、高温の進化した恒星に対応する $T_{\text{eff}} - |\lambda|$ 図の希薄なサンプリング領域におけるデータポイントに寄与する。この結果は、高温の恒星は、放射状エンベロープを持つ恒星と比較して潮汐再整列の効率が低いため、より広い範囲の傾斜角を示すという経験的な傾向と一致している。しかし、著者らは、この系の低い惑星対恒星の質量比（$\sim 0.0008$）が、高質量比の系が整列する傾向にあるのと対照的に、低質量比の系における不整列に関する理論的期待（Rusznak et al. 2025）と一致していることも指摘している。

著者らは、彼らの測定が既存の経験的傾向を支持している一方で、異なる恒星および惑星パラメータ空間におけるスピン・軌道角の分布の散らばりを完全に説明するためには、さらなる観測および理論的研究が必要であると結論付けている。