The age of the Universe from a large sample of the oldest Galactic stars

技術要約：最古の銀河系恒星の大規模サンプルによる宇宙の年齢

問題提起
本論文は、銀河系における最古の恒星の潜在的な年齢（$A_\star$）を推定することにより、宇宙の年齢（$A_U$）を決定することを目的としている。この測定値は、$A_U$ の決定的な下限値となる（ここで $A_U = A_\star + t_f$ であり、$t_f \approx 0.2$ Gyr は最初の長寿命星の形成時間である）。この研究は、「ハッブル・テンション」によって動機付けられている。これは、標準的な $\Lambda$CDM モデルに基づく宇宙マイクロ波背景放射（CMB）から導出されたハッブル定数（$H_0^{CMB} \approx 67.4$ km/s/Mpc）と、局所的な距離梯子による測定値（$H_0^{local} \approx 73.5$ km/s/Mpc）との間の統計的に有意な不一致である。

再結合前における新しい物理学を伴うこれらのテンションへの解決策は、高い $H_0$ を受け入れつつ CMB データとの整合性を維持するために、通常、より若い宇宙（$A_U \approx 12.9 \pm 0.2$ Gyr）を予測する。対照的に、標準的な CMB 校正済み $\Lambda$CDM モデルは $A_U \approx 13.8$ Gyr を予測する。著者らは、結果をいずれかのモデルに偏らせるような強い宇宙論的事前分布を課すことなく、高品質な銀河系恒星の大規模サンプルから $A_\star$ を導出することで、これら競合する予測を検証することを目指している。

手法
著者らは、LAMOST DR7 分光と Gaia eDR3 アストロメトリを組み合わせた Xiang & Rix (2022, 以下 X22) カタログに含まれる 247,103 個の亜巨星（subgiant）のサンプルを利用している。恒星の年齢は、当初、Yonsei-Yale (YY) 恒星進化モデル（isochrones）を用いて推定された。宇宙論的推論のための年齢推定の信頼性を確保するため、著者らは、宇宙論的なバイアスを導入することなく、問題のある年齢決定（例：高年齢域での非現実的な低不確かさを持つもの）を持つ恒星を除去するように設計された、一連の厳格な集団レベルの品質管理カットを適用した。

主な品質管理ステップは以下の通りである：

1. アストロメトリ・カット： パララックス（年周視差）の不確かさを $\sigma_\varpi/\varpi < 0.1$、かつ $\varpi > 0.2$ mas と規定し、サンプルを 5 kpc 以内に限定した。
2. 連星汚染： Gaia データの ipd_frac_multi_peak 値が高い恒星を除外し、境界領域にある連星を除去した。
3. 化学的整合性： 金属量（[Fe/H]）および $\alpha$ 要素組成（[$\alpha$/Fe]）に対する年齢依存のカットを適用した。著者らは、古い恒星は金属欠乏であり、かつ $\alpha$ 元素が濃縮されている必要があるという、銀河化学進化モデルと一致する条件を課している。古いように見えるが金属量が多い、あるいは $\alpha$ 元素が少ない恒星は、測定エラーである可能性が高いため拒絶された。
4. カタログ間検証： YY 年齢を、Gaia DR3 データのみから導出された FLAME 年齢と比較した。年齢と質量の関係に対して線形トレンドラインを当てはめて欠損している FLAME 年齢を補完した後、YY 年齢と FLAME 年齢の間で $3\sigma$ の外れ値除去を行った。

これらのカットの後、最終的な公称サンプルには 155,600 個の恒星が含まれた。著者らは、潜在的な年齢分布 $P(A)$ を再構成し $A_\star$ を推論するために、2 つの主要な統計的手法を用いた：

• 反復再構成（Iterative Reconstruction）： 集団分布が個々の恒星の事後分布を洗練させる事前分布として機能する、非パラメトリックな手法。収束するまで反復される。
• MCMC 再構成： 100 個の年齢ビン（0.1–20.1 Gyr）において $P(A)$ を推論するためのマルコフ連鎖モンテカルロ法。著者らは、ノイズを滑らかにするための勾配ペナルティを含め、分布の対数事後分布をモデル化した。$A_\star$ を推論するために、再構成された確率密度とその不確かさの比（$P/\sigma$）を分析し、分布が減少傾向から平坦な「バイアス床（bias floor）」（アルゴリズムによって $P > 0$ が強制される領域）へと遷移する鋭い「クロスオーバー」点を特定した。

主な結果
解析により、サンプル中の最古の恒星の年齢の堅牢な推定値が得られた：
$$A_\star = 13.73^{+0.18}_{-0.15} \text{ Gyr}$$
この結果は、公称サンプルおよび解析手法から導出されたものである。著者らは、品質管理カットを変化させることで、$A_\star$ の堅牢性をテストした：

• 金属量カット： [Fe/H] の上限を緩和する（より金属に富む恒星を許可する）と、$A_\star$ はわずかに上昇した（最大 $14.02^{+0.18}_{-0.15}$ Gyr）。一方で、上限を厳しくすると $A_\star$ は低下した（最小 $13.31^{+0.21}_{-0.18}$ Gyr）。しかし、著者らは、0.2 dex だけカットを厳しくすることは、銀河の年齢-金属量トレンドと一致する恒星を除去し、Gaia-Sausage/Enceladus 合体のような既知の集団を除外してしまうため、非物理的であると主張している。
• $\alpha$ 元素増強およびカタログ間カット： [$\alpha$/Fe] の下限および YY-FLAME の一致の厳密さ（$2.5\sigma$ 対 $3\sigma$）の変化は、$A_\star$ の値を公称結果と誤差範囲内で一致する結果をもたらした。

「初期化学（early chemistry）」サブサンプル（低 [Fe/H] かつ高 [$\alpha$/Fe] を持つ 1,556 個の恒星）に焦点を当てた二次解析では、$P(A)$ に対するパラメトリックな遅延-$\tau$ モデルを用いて、$A_\star = 13.85^{+0.05}_{-0.04}$ Gyr という一貫した結果が得られ、主たる知見をさらに裏付けた。

意義および主張
本論文は、推論された $A_\star$ が標準的な CMB 校正済み $\Lambda$CDM 予測（$t_f = 0.2$ Gyr と仮定した場合の $A_\star^{CMB} \approx 13.6$ Gyr）と一致している一方で、「初期時刻（early-time）」の解決策（$A_\star^{ETS} \approx 12.7$ Gyr）とは有意な（$>2\sigma$）不一致があることを主張している。

著者らは以下の通り主張している：

1. 新しい物理学への制約： 標準的な $\Lambda$CDM の年齢との一致は、再結合前の新しい物理学のみに依存するハッブル・テンションの解決策に対して疑念を投げかける。なぜなら、それらの解決策は通常、$12.9 \pm 0.2$ Gyr という宇宙年齢を必要とするからである。このような低い年齢を著者らの結果と両立させるには、恒星モデルの不確かさが現在理解されているよりも大幅に大きい必要があるが、本サンプルの低金属量や独立したアステロセイズモロジーによる制約を考慮すると、それは起こり得ないとされている。
2. 手法の堅牢性： 本研究は、集団レベルの品質管理カット（化学組成およびカタログ間の整合性に基づく）が、年齢推定を歪ませる可能性のある外れ値を除去するために不可欠であることを示している。MCMC 再構成の使用により、潜在的な年齢分布とその不確かさを厳密に取り扱うことが可能となった。
3. 先行研究との整合性： この結果は、最古の球状星団や個別の恒星（例：HD 140283）のこれまでの推定値と一致しており、宇宙が約 13.8 億年であることを裏付ける見解を強化している。

結論として、恒星モデルの不確かさは存在するものの、それらは、ハッブル・テンションの初期時刻の解決策が予測する低い年齢と、推論された高い最古の恒星の年齢とを両立させるには不十分であり、したがって標準的な宇宙年齢を支持するものであると本論文は結んでいる。