Radiatively Cooled Binary Mass Transfer: Flow Structure, Luminosities, and L2 Outflows Across Mass Transfer Rates

技術的概要：放射冷却を伴う連星質量移動

問題提起
連星系における質量移動（MT）は、重力波源や中心核崩壊型超新星の形成に先行する恒星進化上の重要な過程である。これらの系のモデル化における主要な不確実性は、MT がどの程度保存的か、すなわち、降着天体が移動した物質をすべて保持するのか、それとも角運動量（AM）を運び去りながら系から質量が逃出するのかという点にある。先行研究（例：Scherbak et al. 2025、「Paper I」）は、断熱仮定の下での急速かつ安定した MT の流体力学的シミュレーションを用い、物質が外ラグランジュ点（$L_2$）を通じて流出し、連星外流出（CBO）を形成することを明らかにした。しかし、Paper I は放射冷却を無視しており、これは中程度の MT 率（$\lesssim 10^{-2} M_\odot \text{yr}^{-1}$）において重要となるはずである。先行モデルにおける冷却の欠如は、保存的および非保存的 MT 領域間の遷移、逃出質量が運ぶ特定の AM、およびこれらの過渡現象の観測的特性（光度と温度）に関する未解決の問いを残した。

手法
著者らは、PLUTO コードを用いた 3 次元ニュートン流体力学シミュレーションを実行し、固定された円軌道上にある $10 M_\odot$ の供給星と $5 M_\odot$ の点質量降着天体間の安定した MT をモデル化した。Paper I と異なり、本研究は近似放射冷却モデルを組み込んでいる。

• 冷却処方: 冷却パワー（単位体積あたり）は、光子が軌道面に垂直に逃出すると仮定して推定される。このモデルは MESA から太陽金属量に対する不透明度テーブルを補間し、ガス圧と放射圧のバランスによってガス温度を求解する。冷却項は内部エネルギーを減少させるが、非物理的な値を防ぐための安全下限（例えば、高密度領域での 2000 K の温度下限、および低密度領域での圧力下限）に従う。
• シミュレーショングリッド: 著者らは軌道半長径 $a$ を変化させ、$\sim 10^{-5}$ から $\sim 10^{-1} M_\odot \text{yr}^{-1}$ の物理的 MT 率（$\dot{M}_{\text{phys}}$）の範囲を調査する。4 つの主要なシミュレーションが実施された：lowest_mdot（$2.0 \times 10^{-5}$）、mid_mdot（$2.2 \times 10^{-3}$）、high_mdot（$2.3 \times 10^{-2}$）、および highest_mdot（$1.5 \times 10^{-1}$）。
• 制約: シミュレーションには明示的な粘性、粘性加熱、または詳細な放射輸送は含まれていない。降着天体はプランマー軟化ポテンシャルで表現される。供給星はケルビン・ヘルムホルツ時間スケールでの熱注入によってロシュ限界を越えるように駆動される。

主要な結果

1. 流れの構造と保存性:

   • 質量は供給星から降着天体周囲の降着円盤へ流入する。高い MT 率（$\gtrsim 10^{-3} M_\odot \text{yr}^{-1}$）の場合、赤道方向に集中した顕著な流出が $L_2$ を通じて発生する。
   • 降着天体が保持する質量の割合 $\beta$ は、MT 率の増加に伴い急激に減少する。$\dot{M} \gtrsim 10^{-2} M_\odot \text{yr}^{-1}$ の場合、MT は極めて非保存的である（$\beta \approx 0$）。より低い率（$\dot{M} \lesssim 10^{-3} M_\odot \text{yr}^{-1}$）では、MT は主に保存的である（$\beta \gtrsim 0.9$）が、少量の物質が依然として $L_2$ を通じて逃出する。
   • 極めて非保存的な MT への遷移は、降着天体の円盤外縁における降着光度がエディントン限界を超えるときに概ね発生する。著者らは、非保存的挙動の閾値が Lu et al. (2023) の半解析的予測よりも高い（約 2〜5 倍）ことを指摘している。

2. 角運動量とエネルギー:

   • いかなる場合においても、流出ガスの特定の角運動量（$h_{\text{loss}}$）は、MT 率に関わらず $L_2$ 点の値（$h_{L2}$）にほぼ等しい値に漸近する。これは、少量の流出であっても連星から効率的に AM を取り除くことを意味する。
   • 流出の漸近的速度は、MT 率が低下するにつれて増加する（高い率では軌道速度の約 10% から、低い率では約 50% へ）。ただし、物理単位では、ほとんどのシミュレーションにおいて速度は $30-60 \text{ km s}^{-1}$ の範囲にとどまる。
   • ベルヌーイパラメータは、高い MT 率における流出は束縛状態であるのに対し、低い率における流出はわずかに非束縛状態であることを示している。

3. 光度と観測的外観:

   • ガスの冷却光度（$L$）と温度（$T$）は MT 率とともに増加する。
   • 低い MT 率では、冷却は降着天体の円盤によって支配される。高い率（$\sim 10^{-1} M_\odot \text{yr}^{-1}$）では、CBO が顕著に寄与し、総光度を倍増させる。
   • 最大冷却光度は、最も高い MT 率で $\sim 10^{39} \text{ erg s}^{-1}$、高い MT 率で $\sim 10^{38} \text{ erg s}^{-1}$ に達する。
   • 実効温度（$T_{\text{eff}}$）は、高い率における降着天体の円盤では $\sim 5 \times 10^4 \text{ K}$（紫外線）から、低い率における CBO では $\sim 10^2 \text{ K}$（中赤外線）の範囲にわたる。著者らは、塵の形成が可能であるが、凝縮に必要な温度がグリッドの範囲を超える半径でのみ達成される可能性があるため、シミュレーション領域内では確定的に予測することが困難であると警告している。

意義と主張
本論文は、放射冷却が断熱モデルと比較して安定した MT の力学を著しく変化させる領域を画定すると主張している。MT 率を体系的に変化させることで、著者らは、外側円盤の超エディントン性によって駆動され、系が保存的から極めて非保存的へと遷移する臨界閾値（$\sim 10^{-3}$ から $10^{-2} M_\odot \text{yr}^{-1}$）を特定した。

著者らは、その結果が MT の安定性と合体前の連星コンパクト天体の軌道進化を調査する集団合成計算に必要な入力データを提供すると主張している。彼らは、$L_2$ 流出が運ぶ特定の AM が一貫して高く、これが連星軌道の長期的進化に影響を与えることを強調している。さらに、計算された光度と温度は、SS 433、W Serpentis、または V1309 Sco のような系で観測される可能性がある、急速な MT に伴う過渡現象の観測的外観に関する予測を提供する。著者らは謙虚に、彼らの結果は単一の質量比と太陽金属量に限定されており、より大きな半径での塵の形成をモデル化し、粘性降着物理を組み込むためには将来の研究が必要であると述べている。