A comprehensive study of $δ$ Scuti-type pulsators in eclipsing binaries: oscillating eclipsing Algols

技術要約：食連星における$\delta$ Scuti型脈動星に関する包括的研究

問題提起
食連星（EB）であり、かつ二重線分光連星（SB2）でもある系は、恒星進化モデルを検証するために不可欠な質量や半径といった絶対的な恒星パラメータを決定するための重要な実験場である。しかし、基本パラメータだけでは、恒星の内部構造や化学組成を明らかにできないことが多い。$\delta$ Scuti型の脈動星は、アステロセイズモロジーを通じて恒星内部を探る経路を提供するが、その特性評価は高品質な視線速度（RV）データの不足によって阻まれている。RV測定がなければ絶対質量を決定できず、これらの系の進化史が不完全なまま残されてしまう。さらに、これらの系の中でも特定のサブセットである振動食アルゴル型（oEA）系は、極端な質量比を持つ半離隔配置を含んでおり、これは質量移動によって恒星の寿命やヘルツシュプルング・ラッセル（HR）図上の位置を変化させた歴史を示唆している。質量移動が中間質量星の脈動特性や進化にどのように影響するかを包括的に理解するには、測光、分光、および進化モデルの同時解析が必要である。

手法
著者らは、CRÉME（Comprehensive Research with Échelles on the Most Interesting Eclipsing Binaries）サーベイによって特定された5つのoEA系（TIC 64437380 (HD 139774)、TIC 82474821 (SW Pup)、TIC 126945917 (HD 202042)、TIC 10756751 (GP Cet)、TIC 37817410 (TZ Eri)）の多角的な解析を提示する。

1. データ取得:

   • 測光: 高精度かつ短周期（2分間隔）のライトカーブ（LC）を、Transiting Exoplanet Survey Satellite（TESS）から取得した。データはlightkurveおよびwotanパッケージを用いて処理され、フラックスの正規化と長期的なトレンドの除去が行われた。
   • 分光: 多時期の高分解能スペクトルが、主にFEROS (MPG-2.2m) および CHIRON (SMARTS 1.5m) を用いて取得され、GP CetについてはHIDESおよびIRCSからの追加データも使用された。

2. 軌道および恒星パラメータの導出:

   • 視線速度: RVは二次元相互相関法（toddar）を用いて計算された。軌道解は、周期、系統速度、および速度半振幅をフィッティングするv2fitを用いて導出された。
   • 分光解析: 主星の大気パラメータ（$T_{\text{eff}}$, $\log g$, $v \sin i$）は、iSpecフレームワークとGSSP（Grid Search in Stellar Parameters）モジュールを使用して決定された。個々の成分スペクトルを分離するために、分光解離（spectral disentangling）はdisentangling_shift_and_addコードを用いて行われた。
   • ライトカーブ・モデリング: phoebe2コードを用いて、ロッシュ幾何学、リムダークニング、および重力ダークニングを組み込んだLCのモデリングを行った。モデリングでは、二次星がロッシュ・ローブを満たす半離隔の制約を利用した。傾斜角、半径比、温度比などのパラメータは、Nelder-Mead法および差分進化アルゴリズムを用いて最適化され、不確実性はマルコフ連鎖モンテカルロ（MCMC）サンプリングを介して推定された。

3. 進化モデリング:

   • MESA（Modules for Experiments in Stellar Astrophysics）コードを用い、5,000個の連星進化モデルからなるグリッドを構築した。
   • グリッドは、初期軌道周期（1.2–5.0日）、初期質量（0.5–2.5 $M_\odot$）、金属量、対流オーバーシュート、および質量移動効率（$\beta$）を変化させた。
   • モデルは、マハラノビス距離を用いて観測された質量および半径と比較され、Case A（中心部水素燃焼中の質量移動）とCase B（中心部水素枯渇後の質量移動）のシナリオを区別した。

4. 脈動解析:

   • 周波数解析は、連星モデルのフィット残差に対して行われた。
   • 残差からは、軌道の高調波および長期的なトレンドが除去された。
   • Period04ソフトウェアを用い、プレホワイトニング（pre-whitening）手順を用いて、振動周波数（0–80 d$^{-1}$）を探索した。信号対雑音比（SNR） $\ge$ 5の周波数を有意とみなした。

主な結果

• 系パラメータ: すべての系において、精密な軌道および恒星パラメータが導出された。これらの系は低い質量比（$q \approx 0.1$–0.3）を示しており、その半離隔アルゴル型の性質を裏付けている。主成分は$\delta$ Scuti型の脈動星（質量 $\sim$1.7–2.4 $M_\odot$）であり、二次成分はより低温で、光度が低く、膨張しており、ロッシュ・ローブを満たしている。
• 進化史: MESAシミュレーションは、すべての系が重大な質量移動を経験したことを示唆している。
  • TIC 82474821: 最良の適合は、非保存的なCase Bの質量移動の歴史を示唆しており、系の年齢は0.5年から3 Gyrの間と推定される。
  • TIC 10756751: モデルは、ほぼ保存的なCase Aの質量移動を示しており、年齢範囲は1.86–8.13 Gyrである。
  • TIC 37817410: フィットは、過剰な太陽金属量とドナーのエンベロープにおける不安定なCNO燃焼サイクルを伴うCase Bのシナリオを指し示している。
  • TIC 64437380 & TIC 126945917: フィットの精度は低かったものの（最大50%の偏差）、結果は依然としてCase Bの分類を支持している。ただし、著者らは、初期パラメータの範囲（特に初期軌道周期と二次星の質量）を精緻化することで、フィットの質を改善できると述べている。
• 脈動: すべての系の残差において有意な振動周波数が検出され、それらが脈動連星であることが確認された。解析により、これらの周波数を連星の軌道信号から分離することに成功した。

意義と主張
本論文は、よく特性化されたoEA系のカタログを増やすことに貢献すると主張しており、恒星の構造および進化の理論モデルをテストするための不可欠なベンチマークを提供している。食連星からの精密な幾何学的制約とアステロセイズミック・データとを組み合わせることで、本研究は、恒星の幾何学と内部構造を同時に制約することの実現可能性を示している。

著者らは、進化モデルが概ね観測された系を再現しているものの、不一致が残っていると控えめに結論づけている。彼らはこれを、対流オーバーシュート、質量移動効率、角運動量喪失、および恒星回転といった物理的扱いの不確実性に起因するものとしている。したがって、本論文は決定的な解決策というよりも、基礎的な一歩として位置づけており、形成シナリオをさらに制約し、質量移動によって誘発される内部構造の変化を理解するためには、モード同定のための多色測光とセイスミック・モデリングを伴う今後の研究が必要であると述べている。本研究は、中間質量星の進化経路を形作る上での連星相互作用の決定的な役割を強調している。