Investigating the System Configuration of Kepler-451 through Orbital Period Variations: Dynamical and Magnetic Interpretations

技術要約：軌道周期変動を通じたケプラー451（Kepler-451）の系構成の調査

問題提起
ケプラー451は、高温のサブダワーフB型（sdB）星と低質量主系列星からなるポスト共通外層（PCEB）である。このような系における食のタイミング変動（ETV）は、しばしば周連星惑星の証拠として解釈されるが、ケプラー451に関するこれまでの研究では相反する結果が得られている。ある解析では単一の木星質量程度の伴天体が示唆された一方で、別の解析では複雑な3惑星系が提案されている。しかし、PCEBにおける多くの提案された惑星構成は、力学的に不安定であるか、あるいはデータ解析の系統的な誤差によるアーティファクトであることが示されている。本研究の主な目的は、ケプラー451の観測基盤を2024年まで拡張し、観察された周期変動が周連星伴天体に起因するものか、あるいは磁気活動サイクルなどの代替メカニズムに起因するものかを厳密に検証することで、これらの相違を解決することである。

手法
著者らは、2004年から2024年までの包括的なデータセットを用いて、更新された$O-C$（観測値－計算値）図を作成した。このデータセットは以下で構成される：

1. 地上観測： 2014年から2024年の間にトルコの3つの望遠鏡（TUG100、ADYU60、ATA050）を用いて得られた68個の一次食ライトカーブ。
2. 宇宙望遠鏡観測： ケプラー（Kepler）ミッションおよびTESS（トランジット系外惑星探索衛星）のデータ（セクター14, 15, 40, 54, 74, 75, 81, 82を含む）。

以前のケプラー解析における手法の違いに対処するため、著者らは2つの異なるデータセットを定義した：

• DS-A： 脈動に対してプリホワイトニング（白色化）処理を行ったケプラーデータ（J. Krzesinski et al. 2020に従う）。
• DS-B： 明示的なプリホワイトニングを行わずに導出されたケプラー食の中央時刻（E. M. Esmer et al. 2022に従う）。

著者らは、周連星伴天体を検証するために、光行差（LTT）効果を用いてタイミング変動をモデル化した。線形エフェメリス、二次項、および2つまたは3つの伴天体に対するLTT項を含む様々な構成をフィットさせるために、マルコフ連鎖モンテカルロ（MCMC）解析を用いた。モデル化されていないノイズを考慮するため、尤度関数には系統的な不確かさパラメータ（$\sigma_f$）を含めた。

惑星起源と非惑星起源を区別するために、本研究では3つの定式化（標準モデル、2ゾーンモデル、およびスピン・軌道結合モデル）を用いたアップレゲート（Applegate）機構を適用し、磁気活動によって観察された変動を駆動するために必要なエネルギー収支（$\Delta E/E_{sec}$）を算出した。最後に、長期的な力学的安定性を、$10^7$年にわたる$N$体積分（REBOUNDパッケージ）を用い、WHFastインテグレータおよびMEGNOカオス指標を使用して評価した。

主な結果

• データセットへの依存性： 優先されるモデル構成は、使用されたデータセットに依存した。
  • DS-Aの場合： 2伴天体LTTモデル（LTT34）が最良のフィットを示した（RMS = 3.23 s, $\chi^2_\nu$ = 1.23）。第3の伴天体を加えても、フィットの有意な改善は見られなかった。
  • DS-Bの場合： 3伴天体LTTモデル（LTT345）が統計的に最良のフィットを与えた（RMS = 2.31 s, $\chi^2_\nu$ = 1.01）。しかし、最も内側の伴天体の半振幅（1.34 s）は、導出された系統誤差（1.81 s）よりも小さく、この信号は校正または観測の系統的な誤差によるものである可能性を示唆している。
• 磁気的起源のテスト： アップレゲート機構のテストによれば、ほとんどの信号は利用可能な磁気エネルギーを超えるエネルギー収支を必要とし、2ゾーンモデルおよびスピン・軌道結合モデルの下では、磁気的起源を否定した。しかし、標準モデルの下では、両方のデータセットにおける外側のLTT項（DS-Aにおける「第4の天体」およびDS-Bにおける「第5の天体」）のエネルギー比（$\Delta E/E_{sec}$）が1未満であったため、これらの特定の信号については磁気的起源を完全に排除することはできない。
• 力学的安定性：
  • すべての信号を重力的伴天体として扱う構成（例：DS-Bの3伴天体モデル）は、強い相互摂動により短期間（$\sim 10^3$–$10^4$年）で力学的に不安定であることが判明した。
  • 外側の信号（潜在的な磁気活動に起因するもの）を除去した場合、残りの構成（DS-Aでは1つの伴天体を持つ連星、DS-Bでは2つの伴天体を持つ連星）は、少なくとも$10^7$年間、力学的に安定していた。

結論および意義
本論文は、ケプラー451系の最も一貫した解釈として、以下のハイブリッドシナリオを結論付けている：

1. 第二世代惑星： 約3.4 AUに周連星惑星の伴天体が存在するという強い証拠がある（DS-Aモデルでは質量 $\sim 3.75 M_{Jup}$、DS-Bモデルでは $\sim 3.7 M_{Jup}$）。この系のPCEBとしての性質と惑星の質量・軌道半径を考慮すると、この天体は第一世代の生存者ではなく、共通外層フェーズ後の残存円盤から形成された第二世代の惑星である可能性が高い。
2. 磁気変調： より長い周期の信号（多体フィットにおける外側の伴天体に相当するもの）は、惑星体ではなく、M型矮星の二次星における周期的な磁気活動に起因する可能性が高い。
3. 系統的な不確かさ： 約400日の周期信号（DS-Bの3体モデルにおける内側の伴天体）は不確実であり、プリホワイトニングされたデータで解析した場合や、系統誤差の範囲と比較した場合に消失または統計的に有意ではなくなることから、データ解析の系統的な誤差に起因するものと考えられる。

本研究は、約3.4 AUにおける周連星惑星の存在はデータと力学的安定性の制約から強固に支持されている一方で、残りのLTT信号の起源については依然として不透明であることを強調している。著者らは、現在の望遠鏡の回折限界に対する投影角分離（2.3–10.9 mas）が小さいため、直接撮像による確認は現在不可能であると述べている。軌道パラメータを精緻化し、惑星起源と磁気的起源を決定的に識別するためには、今後のタイミング観測および視線速度観測が必要である。