以下は、提出された論文「The YREC Stellar Evolution Code: Public Data Release」の技術的な要約です。

論文概要

タイトル: The YREC Stellar Evolution Code: Public Data Release（YREC 恒星進化コード：公開データリリース）
著者: Marc H. Pinsonneault 他（オハイオ州立大学、ハワイ大学、ヴァンダービルト大学など）
日付: 2026 年 3 月 30 日（ドラフト版）

1. 背景と課題 (Problem)

恒星の構造と進化の理論は、質量保存、エネルギー保存、運動量保存、熱力学第二法則などの基礎物理法則に基づいていますが、その解析解を得ることは極めて困難です。1950 年代以降、数値計算コードが開発され、恒星の生涯を解明する上で不可欠なツールとなりました。しかし、現在利用可能な主要な恒星進化コードは限られており（MESA, CESAM, Eggleton STARS, ESTER など）、多くのコードは非公開または特定のユーザーグループに限定されています。

また、既存のコードの中には、モジュール化されたユーザー入力や教育用途への配慮が不足しているもの、あるいは特定の物理過程（回転、磁場、星斑など）の扱いに特化しすぎているものがあります。YREC（Yale Rotating Evolution Code）は 60 年以上の歴史を持つコードですが、その最新バージョンの公開、ドキュメントの整備、および教育・研究コミュニティへのアクセスの容易化が求められていました。

2. 手法とコードの構造 (Methodology)

YREC は、予備主系列（pre-MS）からヘリウム燃焼期までの広範な質量範囲（褐色矮星から高質量星まで）の恒星進化をシミュレーションするための「緩和法（relaxation code）」です。

• 数値手法: ヘニイ（Henyey）法を採用し、恒星構造の 4 つの連立微分方程式を有限差分方程式に変換して解きます。時間ステップごとに構造変数を反復計算し、収束させます。
• 境界条件: 大気モデル（Gray 大気、Kurucz、Allard など）を上部境界条件として使用し、内部解と外層（エンベロープ）の積分をマッチングさせる「フィッティングポイント」手法を採用しています。MESA と異なり、光球より深い位置をフィッティングポイントに設定することで計算の安定性を高めています。
• 物理過程:
  • 状態方程式 (EOS): OPAL、SCV、Yale EoS から選択可能。
  • 不透明度: 原子（OPAL, OP）、分子（Ferguson et al.）、伝導（Cassisi et al.）のテーブルを使用。
  • 核反応率: 水素燃焼には Bahcall & Pinsonneault (1992) 法、ヘリウム燃焼には Caughlan & Fowler (1988) 法などを使用。
  • 対流: 混合長理論（MLT）とシュワルツシルト基準を使用。オーバーシュート（step overshoot）や半対流（semi-convection）のオプションも備えています。
  • 拡散: 重力沈降と熱拡散（微小拡散）を考慮可能。
  • 回転と混合: 擬似 2 次元近似を用いて回転による構造変化を計算。角運動量の輸送（子午面循環、せん断不安定、ゴールドシュタイン・シュバルツ・フリッケ不安定など）と、磁気風による角運動量損失（ブレーキング）をモデル化できます。
  • 星斑: 星斑によるフラックスの遮蔽と表面境界条件の変更を考慮し、半径の膨張やリチウム枯渇への影響をシミュレート可能です。
• 実行モード: 通常の進化計算（時間依存）と、既知の解を基準としたスケーリング計算（ZAMS, ZAHB 等での再スケーリング）の 2 モードがあります。

3. 主な貢献 (Key Contributions)

1. YREC コードの公開リリース: ソースコード、入力ファイルライブラリ、ドキュメント、テストスイート、スクリプトを GitHub および Zenodo で公開しました。
2. 教育ツールの強化: 軽量で高速な設計を活かし、WebAssembly 化されたインタラクティブな Web 版（yreclab）を提供。学生がリアルタイムで恒星進化をシミュレート・操作できる環境を整備しました。
3. 包括的なテストスイート: 太陽モデル、低質量星、褐色矮星、高質量星、回転モデル、星斑モデルなど、多様な物理シナリオに対するテストケースを提供し、コードの回帰テストとユーザーの検証を支援します。
4. 補助ツールの提供:
   • rotevol: 非回転モデルから角運動量進化を高速に追跡するトレーサーコード（YREC 本編の 100 分の 1 の速度）。
   • kiauhoku / ISTARY: 恒星モデルグリッドから等齢線（isochrones）を生成するツールセット。
   • スクリプト: 太陽モデルの較正、パラメータグリッドの自動実行、出力ファイルの解析（Python/Pandas 連携）などを支援するツール群。

4. 結果と検証 (Results)

• 精度と性能:
  • 太陽モデルにおいて、音速プロファイルやニュートリノフラックスが観測値と良好に一致することを確認しました。
  • 低質量星から高質量星までの進化軌道（HR 図）を、MIST や PARSEC などの他の主要コードと比較し、主系列帯や赤色巨星分枝（RGB）での一致を確認しました。
  • 計算速度の比較では、MESA と同等のメッシュ解像度において、YREC は ZAMS までの計算で約 8 倍、TAMS（主系列終点）から RGB 先端までは約 6 倍高速であることを示しました。簡易的な状態方程式（Yale EoS）を使用すればさらに高速化可能です。
• 物理的発見の再現:
  • 回転: 回転による混合と角運動量損失のモデルが、星団の自転周期 - 年齢関係（Gyrochronology）を再現できることを示しました。
  • 星斑: 星斑の存在が低質量星の半径膨張やリチウム枯渇パターンに与える影響を定量化しました。
  • ヘリウム燃焼: 低質量星におけるヘリウムフラッシュ後の水平分枝（HB）進化において、オーバーシュートや半対流の扱いがコアの成長と光度に与える影響を詳細に議論しました。
• 制限事項:
  • ヘリウムフラッシュ（非対称な核反応）を直接通過させる機能は現在のバージョンでは削除されており、2 ステップ手順（RGB 先端から ZAHB への再スケーリング）が必要です。
  • 高質量星（16 太陽質量以上）では放射圧駆動の星風による大気不安定が発生し、計算が失敗する場合があります。
  • 褐色矮星の低温・高密度領域では、状態方程式や境界条件のテーブルの限界により、特定の質量・年齢範囲で注意が必要です。

5. 意義と将来展望 (Significance)

この論文は、YREC コードを現代のオープンソースコミュニティに再導入し、研究と教育の両面でその利用を促進する重要なステップです。

• 研究面: 回転、磁場、星斑などの複雑な物理過程を柔軟に組み合わせて恒星進化を研究できる強力なツールを提供します。特に、星斑が若い星の進化に与える影響や、回転混合が化学進化に及ぼす効果の解明に寄与します。
• 教育面: Web 版（yreclab）の導入により、学生が高度な数値シミュレーションを直感的に体験できるようになり、恒星構造と進化の教育の民主化に貢献します。
• コミュニティ: 公開されたスクリプトやツール（rotevol, kiauhoku など）は、大規模なモデルグリッドの生成や等齢線の作成を効率化し、天文学研究のワークフローを改善します。

今後は、対流境界の扱いや質量損失の物理モデルの更新、ハードコーディングされた部分の一般化など、コードのさらなる改良と拡張が予定されています。