Chronos: Towards a self-consistent and absolute stellar age scale. I. A Bayesian hierarchical lithium-age model: Validation on the Pleiades cluster

技術的サマリー：Chronos：I. ベイズ階層リチウム年齢モデル：プレアデス星団における検証

問題提起
自己整合性があり絶対的な恒星年齢スケールを確立することは、天体物理学における中心的な課題であり続けている。特に、異なる年代測定手法間の系統的な相互較正が欠如しており、これらがしばしば最大で 2 倍の不一致をもたらすことがその主な要因である。リチウム枯渇境界（LDB）は、低質量恒星集団に対して堅牢で物理的に根拠のあるクロノメーターを提供するが、その適用は不均一な温度スケールと天体物理学的な分散によって現在制限されている。具体的には、完全対流型の超低温矮星（UCD）と部分的に放射型の FGK 星との間の遷移は、しばしば鋭い温度カットとして扱われ、領域遷移全体にわたって不確実性を一貫して伝播させることに失敗している。さらに、FGK 集団は恒星自転によって駆動される顕著なリチウム分散を示し、高速自転星は低速自転星と比較して増大したリチウム量と増加したばらつきを示す。標準的な手法は、しばしば単一値の処方箋や経験的フィットに依存しており、この内在的な多様性やリチウム枯渇の非線形物理を完全に捉えていない。

手法
著者らは、有効温度（$T_{\text{eff}}$）、リチウム量（$A_{\text{Li}}$）、および大域的な星団年齢を共同モデル化することで恒星集団の年齢を推論するように設計された、ベイズ階層モデル（BHM）であるChronosを提示する。この枠組みは、以下の 3 つの中核的構成要素を統合している：

1. 階層構造：このモデルは、観測データ（$T_{\text{eff,obs}}, A_{\text{Li,obs}}$）を、潜在する真の恒星パラメータのノイズを含む実現として扱う。恒星領域間の硬い境界を避ける代わりに、完全対流型と部分的に放射型の恒星間の遷移をモデル化するために、滑らかなロジスティック関数を使用する。
2. ニューラルネットワーク代理モデル：BT-Settl 進化モデルで訓練された多層パーセプトロン（MLP）が、理論的リチウム枯渇曲線の代理モデルとして機能する。このネットワークは、年齢と$T_{\text{eff}}$の関数として平均理論リチウム量（$\mu_{A_{\text{Li}}}$）を予測し、物理モデルの高速かつ連続的な近似を提供する。
3. 自転を考慮した混合モデル：FGK 星における自転誘起のリチウム増強を考慮するため、このモデルは真のリチウム量に対して 2 成分混合分布を採用する。これは、「低速自転星」成分（理論予測を中心とした正規分布としてモデル化）と「高速自転星」成分（振幅オフセット$\Delta A_{\text{Li,rot}}$と追加分散を伴う歪正規分布としてモデル化）から構成される。ある恒星が高速自転集団に属する確率は、そのスペクトル領域と大域的な自転割合パラメータ（$p_{\text{rot}}$）の関数として扱われる。

推論は、PyMC に実装されたノ・Uターン・サンプラー（NUTS）を用いたマルコフ連鎖モンテカルロ（MCMC）によって行われる。このモデルは、同じ階層仮定から生成された合成データセットを用いて検証され、122 個のプレアデス星団メンバーからなる均質なデータセットに適用された。

主要な貢献

• 統合枠組み：Chronos は、リチウムに基づく恒星年代学を、任意の温度閾値を課すことなく、星団年齢、内在的分散、および自転パラメータの同時推論を可能にする、一貫した階層推論問題として再構成する。
• 確率的領域遷移：このモデルは、UCD と FGK 星間の遷移を連続的な確率関数を通じて明示的に処理し、LDB 全体にわたって一貫した不確実性の伝播を確保する。
• 自転モデリング：自転効果を組み込む混合モデルを採用することで、この枠組みは FGK 星で観測される非対称なリチウム分散を考慮し、内在的な天体物理学的ばらつきと自転誘起の増強を区別する。
• 代理モデル化：進化モデルで訓練されたニューラルネットワークの使用により、複雑な理論グリッドと統計的モデリングの間のギャップを埋めながら、ベイズ推論ループ内でリチウム枯渇物理を効率的に評価することを可能にする。

結果
プレアデス星団に適用された Chronos は、Age = 124.53$^{+3.34}_{-2.70}$ Myrを中心とする事後年齢分布をもたらす。この結果は、古典的な LDB 推定値および最近の等年齢線フィット結果と一致するが、より tight な不確実性制約を特徴とする。このモデルは、20–400 Myr の範囲で合成テストにおいて入力年齢を成功裡に回復し、LDB 感応領域（$T_{\text{eff}} < 3700$ K）における対象数が限られていても堅牢であることを示した。

収束統計（$\hat{R} \approx 1$）、残差分析、事後予測チェック（PPC）を含む診断チェックは、モデルの安定性を確認した。分析により、推定された年齢は主に急峻な LDB 遷移付近の恒星（$T_{\text{eff}} \sim 2800$ K）によって駆動されており、自転パラメータとはほとんど相関していないことが明らかになった。これは、年齢決定が自転誘起のリチウム増強の特定のパラメータ化に対して堅牢であることを示している。また、このモデルは高速自転星の割合（$p_{\text{rot}} \approx 0.32$）およびリチウム増強の大きさ（$\Delta A_{\text{Li,rot}} \approx 1.78$ dex）を成功裡に制約した。

意義
本論文は、Chronos を、自己整合性があり絶対的な恒星年齢スケールに向けた方法論的に堅牢な一歩として位置づけている。リチウムに基づく年代測定を階層推論問題として定式化できることを実証することで、この研究は、若い恒星集団（1–600 Myr）における年齢決定のための柔軟かつ統計的に厳密な枠組みを提供する。著者らは、このアプローチが複数の年齢診断の相対的な制約能力を自己整合的に評価することを可能にすると強調している。現在の検証はプレアデスと太陽金属量モデルに限定されているが、この枠組みは、金属量を次元として含めること、および将来の反復において追加の年齢感応観測量を取り入れることを拡張可能に設計されている。この研究は、すべての年代測定不一致を解決すると主張するものではないが、リチウム診断をより広範で統合された年齢スケールに統合するための、物理的に解釈可能で統計的に健全な基盤を確立するものである。