Disk-Regulated Mass Transfer Between Rotating Non-Degenerate Stars: Insights from Be and sdOB Binaries

技術要約：回転する非縮退星間における円盤制御型の質量移動

問題提起
非縮退星間における質量移動は、バイナリ進化において極めて重要でありながら、未だ十分に理解されていないプロセスである。中心的な不確実性は、降着効率（転送された質量のうち、受容星が保持する割合）にある。標準的な連星個体群合成（BPS）および詳細な恒星進化シミュレーション（例：MESA）では、通常、「回転制限型」の降着処方箋が採用されている。この枠組みでは、受容星が臨界破断回転速度に近づくと降着が制限されるため、多くの場合、非常に低い降着効率を招くことになる。

しかし、Be星連星、特にBe星＋準白色矮星型O/B型（sdOB）系からの観測的制約は、著しく高い降着効率（しばしば $\gtrsim 30\%$ から $>50\%$）を示唆している。既存の回転制限モデルは、これらの系におけるBe星の質量を系統的に過小評価しており、Be+sdOB連星の観測されたパラメータ空間を再現できていない。先行研究によれば、「回転制限型」および「熱制限型」の処方箋のいずれも、単独では観測されたBe+sdOB個体群から導き出される効率の要件を満たすことができないことも示されている。

手法
著者らは、これらの不一致に対処するため、物理的根拠に基づいた「円盤―恒星結合（disk-star coupling）」質量降着処方箋を提案する。その核心となる物理メカニズムは、降着円盤または放出円盤が受容星へと輸送される角運動量を制御し、受容星が臨界回転限界付近で回転している場合でも質量降着を継続させるという仮説である。

1. 角運動量の制御： Popham & Narayan (1991) の定常状態薄円盤解に基づき、著者らは受容される恒星の比角運動量 ($j_{acc}$) を、恒星の回転率 ($\Omega/\Omega_{crit}$) の関数としてモデル化している。

   • 低回転率において、$j_{acc}$ は恒星表面におけるケプラー回転値に近い。
   • 恒星が閾値回転率（採用された値は $\Omega/\Omega_{crit} \approx 0.9$）に近づくと、円盤のせん断応力により、角運動量が恒星から円盤へと外向きに輸送される。その結果、$j_{acc}$ は急激に減少し、負の値にさえなり得る。これにより、恒星が臨界回転を超えないようにしながら、質量降着を継続することが可能となる。
   • 著者らは、1次元シミュレーションにおける数値的安定性を確保するため、厳密な閾値によって生じる不連続性を回避するように、平滑化された関数形式で $j_{acc}$ を実装している。

2. 保存則と効率の計算： このモデルは、恒星から抽出された過剰な角運動量が円盤を通じて輸送され、円盤の外縁部 ($R_{out}$) での質量放出を介して除去されると仮定している。質量保存および角運動量保存を適用することで、降着効率 ($\eta = \dot{M}_{acc}/\dot{M}_{tr}$) が自己整合的に計算される。この効率は、円盤端における比角運動量と降着半径における比角運動量の比に依存する。

3. シミュレーション・フレームワーク： 著者らは、この処方箋をPOSYDONフレームワーク内に実装し、MESA (Modules for Experiments in Stellar Astrophysics) を用いて連星モデルを進化させた。

   • 超臨界回転を防止するための暗黙的な風（wind）スキームを無効にし、代わりに、超臨界領域における回転強化風と減少した降着率に頼ることで、亜臨界回転を回復させる手法を採用した。
   • 太陽組成を用い、主星質量 ($2-9 M_\odot$)、質量比 ($0.4-0.95$)、および軌道周期を変化させた連星モデルのグリッドをシミュレートした。
   • また、コア・オーバーシューティングの取り扱いが最終的な成分質量に与える影響についても調査した（標準的なPOSYDONの設定と、低減されたオーバーシューティング処方箋を比較）。

主な貢献と結果

• 解析的推定： 著者らは、降着効率の解析的な式を導出した。円盤降着レジームにおいては、効率は恒星の回転、質量比、および恒星半径とロッシュ・ローブ半径の比に依存することを示している。直接衝突（direct-impact）レジームにおいては、回転と質量比に依存する。全体的な質量加重降着効率は、比角運動量が急激に減少する閾値回転率付近で予想される値に近いことが判明した。

• 連星進化の比較：

  • 回転制限モデル： この標準的なモデルでは、受容星は転送された質量のわずかな割合（テストケースでは $\sim 0.2 M_\odot$）を降着した後に臨界回転に達するため、効率が低くなる（ほとんどの系で $<10\%$）。これにより、結果として得られるBe星は、観測値よりも系統的に質量が小さくなる。
  • 円盤―恒星結合モデル： このモデルでは、受容星はより多くの質量を降着することが可能であり（例：テストケースの急速相において $\sim 1.6 M_\odot$）、効率は $\sim 40-60\%$ に達する。円盤が角運動量を抽出することで、恒星は臨界回転閾値付近に留まりながら、膨張し、かつ降着を続けることができる。

• 観測の再現：

  • Be+sdOB系に適用した場合、回転制限モデルは、特に広軌道（Case B）の系において、観測と一致する質量を持つBe星を生成できない。
  • 円盤―恒星結合モデルは、観測範囲の大部分をカバーする質量を持つBe星を生成する。
  • 著者らは、最終的な成分質量が想定されるオーバーシューティングに敏感であることを発見した。この新しい降着処方箋と低減されたオーバーシューティングを組み合わせたモデルは、標準的な高オーバーシューティングモデルよりも、Be+sdOB系の成分質量に関して観測結果に近い結果をもたらす。

• 軌道周期： 円盤―恒星結合モデルは質量保持を増加させる一方で、得られる軌道周期は回転制限モデルと同等か、あるいはわずかに長くなる。これは、より保守的な質量移動によってシステムからの総質量損失が減少するため、より極端な質量比が生じ、それが軌道膨張を加速させ、円盤端からの角運動量損失を部分的に相殺するためである。

意義と主張
本論文は、観測されたBe+sdOB連星の特性と現在のシミュレーション結果との間の不一致は、円盤制御型の質量移動モデルを採用することで解決できると主張している。著者らの主張は以下の通りである：

1. 回転制限型の降着処方箋は、Be+sdOB系におけるBe星の質量を系統的に過小評価している。
2. 角運動量輸送の制御を可能にする円盤―恒星結合メカニズムは、観測的制約（多くの系で効率が $>50\%$ であるという要件）と一致する高い降着効率を可能にする。
3. この新しい降着処方箋と低減されたコア・オーバーシューティングを組み合わせることは、Be+sdOB連星の質量と軌道周期に対して、より物理的に一貫した説明を提供する。

著者らは、実際の効率は円盤構造、磁場、および円盤形成と質量放出（例：L2点やホットスポット経由）の具体的な条件といった不確実な要因に依存することを認めている。しかし、円盤―恒星結合モデルは、潮汐によるスピンダウンが非効率的な広軌道連星において、大質量Be星が存在することを説明するための、必要な物理的メカニズムを提供すると結論付けている。