Wolf-Rayet stellar evolution models with improved treatment of the atmosphere

技術的サマリー：大気処理を改善したウォルフ・ライエ星の進化モデル

問題提起
大質量星の進化モデルは、歴史的に外境界条件を定義するためにエディントン灰色大気近似に依存してきた。この近似は、非局所熱平衡（non-LTE）効果、非灰色放射輸送、恒星風、および線遮蔽といった、大質量星大気に固有の複雑な物理現象を捉えきれない。その結果、標準的なモデルは、詳細な大気コードを用いたスペクトル分析から導き出されたものよりも著しく高い有効温度（$T_{\text{eff}}$）をウォルフ・ライエ（WR）星に対して予測する。この不一致は、進化済みの大質量星をヘルツシュプルング・ラッセル（HR）図上に正確に配置するのを妨げ、銀河進化におけるそれらの役割の定量的理解を制限している。

手法
著者らは、最先端の大気モデルを大質量星用の恒星進化コードSTAREVOL（バージョン 3.40）に直接統合する方法を開発した。中核的な戦略は、進化コードと放射輸送コードの時間のかかるリアルタイム結合を回避し、事前に計算された大気のグリッドを使用することである。

1. 大気グリッドの構築： CMFGENを用いて、非 LTE、球対称の大気モデルのグリッドを計算した。グリッドパラメータには、光度（$L$）、有効温度（$T_{\text{eff}}$）、表面重力（$g$）、質量損失率（$\dot{M}$）、および表面化学組成が含まれた。
2. 補間スキーム： 進化中に境界条件を提供するため、著者らは Delaunay 三角分割に基づく補間手法（delaunaysparse アルゴリズムを使用）を採用した。補間は、整合光学深さ $\tau_{\text{match}} = 25$ における温度 $T(\tau_{\text{match}})$、$\log(g)$、および $\log(\dot{M})$ で定義された正規化されたパラメータ空間で行われた。光度は補間変数として使用されず、進化モデルの光度に一致させるために補間後に修正された。
3. 進化への結合：
   • 進化モデルの数値表面は、拡散近似の有効性を確保するのに十分な深さである $\tau_0 = 10$ に設定された。
   • 各時間ステップにおいて、コードは現在の恒星パラメータに基づいて CMFGEN グリッドから温度、密度、および半径のプロファイルを補間した。
   • これらの補間されたプロファイルは、内部構造方程式（特に $\tau_0$ における温度と密度）の境界条件として機能した。
   • 補間された半径プロファイルは、光球半径（$R_{\text{eff}}$）を定義し、次の反復のための有効温度（$T_{\text{eff}}$）と表面重力を計算するために使用された。
4. モデルパラメータ： この手法は、太陽金属量における初期質量 50 $M_\odot$ および 60 $M_\odot$ の非回転モデルに適用された。質量損失 prescriptions は、OB 相、LBV 相、WR 相の各段階に合わせて調整され、特に WR 相については、エディントン極限付近で流体力学的に一貫した大気構造を計算できるように特定の修正が加えられた。

主要な結果

• 有効温度の低下： 詳細な大気の取り込みは、WR 相中の予測される有効温度を著しく低下させる。古典的なエディントン灰色モデルと比較して、CMFGEN ベースのモデルは、WR 進化時間の約半分において 20 kK を超える $T_{\text{eff}}$ の低下を示す。
• 光球の膨張： $T_{\text{eff}}$ の低下は、主に恒星風によって引き起こされる光球半径の膨張（特定の事例では最大で約 80%）に起因しており、これは静的な平面平行灰色大気には存在しない特徴である。
• 内部構造への影響なし： 表面特性に著しい変化があったにもかかわらず、内部構造（高光学深さにおける温度、密度、および半径のプロファイル）およびモデルの化学進化は、古典的な計算と実質的に同一のままである。大気から深部内部へのフィードバックは無視できる。
• 観測との一致： 修正された軌跡は、銀河系および大マゼラン雲（LMC）の WR 星（WN 型および WC 型の両方）の観測と非常に良い一致を示す。具体的には、モデルは WNh 星の位置を成功裡に再現し、観測データと整合する有効温度（例えば、WC 星が約 85,000 K 付近に現れること）で水素を含まない WN 星および WC 星への遷移を再現する。
• ポストプロセッシングの同等性： 著者らは、古典的モデルに対して同じ補間手法を用いて、古典的モデルの内部パラメータに基づいて $T_{\text{eff}}$ と $R_{\text{eff}}$ を再計算するポストプロセッシング補正を適用した場合、直接大気を組み込んで計算された軌跡とほぼ区別がつかない進化軌跡が得られることを示した。これは、主な影響が内部進化ではなく表面の定義にあることを確認するものである。

意義と主張
本論文は、進化済みの大質量星を HR 図上で正しく再現し、理論的な $T_{\text{eff}}$ と観測的なスペクトル分析の間の長年の不一致を解決するために、詳細な大気モデルの取り込みが必要であると主張している。ただし、著者らはこの改善の範囲については控えめである：

• この手法はモデルの外観（表面特性）を修正するが、内部進化や化学進化を変化させるものではない。
• この研究は、太陽金属量を持つ 50 および 60 $M_\odot$ の非回転モデルに限定されている。
• 使用された質量損失率は、大気計算における数値的安定性を確保するために 3 分の 1 に削減された。したがって、報告された温度低下は、水素を含まない WR 相における期待される効果の下限を表している。
• 著者らは、亜表面対流層（例えば、鉄の不透明度の突起）の処理や、WNh 風の起動メカニズムの特定の物理に関する限界を認めており、すべての進化段階にわたって完全に自己一貫した処理を達成するためにはさらなる研究が必要であると指摘している。

結論として、本研究は、詳細な大気が WR 星の観測可能な特性を定義する上で重要である一方で、これらの星の内部進化は大気境界条件の選択に対して頑健であることを確立し、多くの文脈において直接結合に代わる効率的な代替手段としてポストプロセッシング補正を使用することを正当化している。