Envelope Inflation and outflow Driven by Energy Deposition in Massive Stars

技術要約：大質量星におけるエネルギー沈着に起因する外層の膨張と流出

問題提起
進化した大質量星は、しばしば高質量放出や流出を示すことがあり、これらはルミナス・ブルー・バリアブル（LBV）や超新星前駆星（pre-SN）の進化過程で観測される現象である。超エディントン風や連星相互作用がその駆動メカニズムとして提案されているが、衝撃的なエネルギー沈着に対する恒星外層の具体的な流体ダイナミクス的応答については、まだ十分に理解されていない。先行研究では、連続的なエネルギー沈着やパラメトリックなエネルギー注入が探求されてきたが、局所的かつ衝撃的なエネルギー入力によって、どのようにして恒星の外層が静水圧的な膨張から動的な流出へと移行するかという詳細な調査が欠けていることが多い。本研究は、大質量星の外層が突然の局所的なエネルギー沈着に対してどのように反応するかという空白を埋めるものであり、特に、このエネルギーが外層を単に膨張させるのか、それとも物質を解離させて流出を形成させるのかという条件を定量化する。

手法
著者らは、1次元流体恒星進化コードであるMESAを用い、ポストメインシーケンス相にある$70 M_\odot$の星をシミュレートしている。本研究では、理想化された「サーマル・ボム（熱爆弾）」アプローチを採用しており、これは、外層の結合エネルギーが仮説上の伴星の重力エネルギーと等しくなる臨界半径（$R_{crit} \approx 23.21 R_\odot$）に位置する、狭いラグランジュ質量範囲（半径幅 $\sim 0.5 R_\odot$ に相当）内に、エネルギーを瞬時に沈着させるものである。

シミュレーションでは、以下の2つの異なる物理的レジームを比較している：

1. ケース1（静水圧的）: mlt++対流処理を含み、電子散乱の増強を無視している。この設定は流出の生成を抑制し、準静水圧的な外層膨張の研究を可能にする。
2. ケース2（流体的）: mlt++処理を無効化し、電子散乱不透明度を人工的に100倍に高めている。この構成は、放射圧支配の外層を強制することで、物質がエネルギー的に解離することを可能にし、動的な応答を検証する。

本研究では、3つのパラメータを異なるケース間で系統的に変化させている：

• ケース2A: 加熱領域の幅が $0.5 R_\odot$ であるベースラインとなる流体的実行。
• ケース2B: 加熱領域の幅を $0.8 R_\odot$ に広げた流体的実行。
• ケース2C: ケース2Aと比較して、加熱領域の幅を一定に保ったまま、より高い沈着エネルギー（$E_{dep}$）を与えた流体的実行。

シミュレーションでは、結合された物質と解離した物質への遷移に焦点を当てながら、1.2年間の期間にわたって全エネルギー、密度、質量放出率、および運動エネルギーの進化を追跡している。

主な結果

• 静水圧的応答 vs 流体的応答: 静水圧的なケース（ケース1）では、エネルギー沈着によって外層は大幅に膨張する（半径は約1.09倍増加する）が、重力的に結合されたままである。流出は生成されず、外層の全エネルギーは負のままとなる。対照的に、流体的なケース（ケース2A, 2B, 2C）では、沈着されたエネルギーが外層を解離させるのに十分であり、持続的な流出の形成につながる。
• 流出の発現: 流体モデルにおいて、流出は全エネルギーが正となる半径 $R_{out} \approx 58.1 R_\odot$ で開始される。この半径はエネルギー沈着部位（$23\text{--}23.5 R_\odot$）とは異なっており、エネルギーが外側へ伝播し、より大きな半径で物質を解離させることを示している。
• 加熱領域の幅の影響（ケース2B）: 加熱領域の幅を $0.5 R_\odot$ から $0.8 R_\odot$ に広げると、総解離質量と平均質量放出率（$\dot{M}$）が高くなる。しかし、より広い加熱領域は、より狭い局所的な沈着と比較して、恒星の膨張をわずかに抑制する。広範な領域はより多くの質量範囲にエネルギーを与え、より持続的な質量負荷型の流出を作り出す一方で、狭い領域は急峻な圧力勾配を生み出し、より衝撃的ではあるが質量の少ない放出を駆動する。
• 沈着エネルギーの影響（ケース2C）: 沈着エネルギーを増やすと、より激しく急速な流出が引き起こされる。直感に反して、これは、より低いエネルギーのケース2Aよりも、最終的な恒星の構成がやや高温で、光度が低く、膨張も小さいという結果をもたらす。著者らは、これを、沈着されたエネルギーが運動エネルギーへと急速に変換され、外層の長期的な膨張を駆動する前にシステムからエネルギーが奪われるためであると説明している。
• 二峰性挙動: 高エネルギーのケース2Cでは、流出プロファイルが二峰性の構造を示しており、これはバースト中に複数の放出エピソードが発生したことを示唆している。この特徴は、低エネルギーまたは広幅のケースには見られない。
• 質量放出と速度: シミュレーションによる質量放出率は最大 $\sim 10^{-2} M_\odot \text{yr}^{-1}$ に達し、流出速度は $\sim 4.5 \times 10^3 \text{ km s}^{-1}$ である。沈着されたエネルギーのうち、放出物によって持ち去られる割合は時間の経過とともに減少する（例：ケース2Aでは、$t=0.4$ 年時点の $\sim 24\%$ から、$t=1.2$ 年時点では無視できるレベルまで減少）。

意義および主張
本論文は、衝撃的なエネルギー入力に対する大質量星の外層の一般的な流体ダイナミクス的応答の「ベースケース」を確立したと主張している。主要な貢献は、応答の性質（膨張か流出か）が、対流の扱いと不透明度、およびエネルギー沈着の空間分布と大きさに極めて敏感であることを示した点にある。

著者らは、これらのモデルが理想化されたものであり、特定の観測事象（$\eta$ Carinaの大噴出など）を再現したり、特定の物理的トリガー（連星力学や波の輸送など）を自己整合的に再現することを目的としていないことを強調している。代わりに、本研究は外層の基本的な流体ダイナミクス的挙動を分離して抽出している。結果は、エネルギー沈着が強い流出を駆動し得る一方で、質量放出の効率と外層膨張の間の関係は、沈着パラメータに対して非線形に依存することを示唆している。具体的には、より高いエネルギー入力は、より急速なエネルギー損失を放出物を通じて引き起こし、その後の恒星の膨張を制限する場合がある。著者らは、これらの知見が、定常的な光度駆動型の質量放出に関する既存の研究を補完するものであると述べている。