Forward-modelling Milky Way Cepheids: selection effects and physical priors in the Gaia-HST calibration

この論文は、銀河系セファイド変光星の Gaia 視差を用いた距離梯子の較正において、銀河の構造と HST 観測の選択効果を明示的に組み込んだ完全な前方モデル化ベイズ解析を採用することで、従来の一様体積事前分布に基づく手法がもたらす系統的バイアスを解消し、Hubble 定数の局所測定値と初期宇宙推測値との間の「Hubble 緊張」を支持する結果を得たことを示しています。

要約

1. 背景：宇宙の「ものさし」は狂っている？

私たちが宇宙の大きさを測るには、何段階もの「ものさし」が必要です。

1. 1 段目： 銀河系内の星（セファイド変光星）の距離を測る。
2. 2 段目： その星を使って、近くの銀河の距離を測る。
3. 3 段目： 近くの銀河を使って、遠くの宇宙の膨張速度（ハッブル定数）を計算する。

しかし、ここで大きな問題が起きました。

• 方法 A（宇宙の赤ちゃんの頃の写真から計算）： 宇宙の膨張速度は「約 67」。
• 方法 B（今の宇宙の星を使って計算）： 宇宙の膨張速度は「約 73」。

この 2 つの値が一致しないことを**「ハッブルの緊張（Hubble Tension）」と呼びます。これは、物理学の大きな謎です。「もしかしたら、新しい物理法則があるのか？」と期待されたのですが、実は「ものさし（距離測定）の計算方法に、ちょっとした『見落とし』があったのではないか？」**という疑念がありました。

2. この論文の核心：「見えないフィルター」を忘れない

この研究チームは、**「私たちが観測している星は、宇宙の全人口の代表ではない」**という点に注目しました。

例え話：「街の人口調査」

想像してください。あなたが「ある街の平均身長」を調べるために、街を歩いているとします。

• あなたの調査方法（従来のやり方）： 「見つけた人」の身長を測って平均を出した。
• 問題点： あなたが歩いたのは「公園」だけでした。公園には「背の高い子供」しかいません（背の低い子供は公園に入れない、あるいは見えない）。
• 結果： 「この街の平均身長は 180cm だ！」と間違った結論を出してしまいます。

この論文の著者たちは、**「公園（観測データ）には、特定の条件（明るさや距離）を満たす星しかいない」という「選択のフィルター」**を、計算に組み込まなければいけないと指摘しました。

3. 彼らがやったこと：「未来から過去を予測する」

従来の計算方法は、「観測したデータ」を元に「星の性質」を推測する（後ろ向き）やり方でした。
しかし、この論文では**「Forward-modelling（前方モデリング）」**という新しいアプローチを取りました。

例え話：「料理のレシピ作り」

• 従来の方法： 出来上がった料理（観測データ）を見て、「たぶん塩は小さじ 1 だったな」と推測する。
• この論文の方法： 「もし塩を小さじ 1 入れたら、どんな味がするだろう？」とシミュレーション（料理）を作り、実際に食べた味（観測データ）と比べてみる。
  • もし味が合わなければ、「塩の量（パラメータ）」を調整して、またシミュレーションする。
  • これを何千回も繰り返して、最も確実な「レシピ（正しい距離と性質）」を見つける。

この方法のすごい点は、**「公園（観測データ）に誰が入れるか（選択フィルター）」と「街の地形（銀河の形）」**を、シミュレーションの中に完璧に再現したことです。

4. 重要な発見：「均一な分布」という思い込みの罠

以前の研究（Högås & Mörtsell さんという方々）は、「星は宇宙に均一に広がっているはずだ」という仮定（均一な分布）を使って計算しました。これにより、ハッブル定数の矛盾が小さくなった（2σ 程度まで減った）と報告されました。

しかし、この論文はそれを**「誤り」**だと証明しました。

• なぜ誤りなのか？
  • 「均一に広がっている」と仮定して計算すると、**「遠くの星もたくさんいるはずだ」**という前提になってしまいます。
  • でも、実際には「遠くの星は暗すぎて見えない（フィルター）」のです。
  • この「見えない星」を無視して計算すると、**「星は実際より明るく見える」**という誤った結論になり、結果として「宇宙の膨張速度」の計算が狂ってしまいます。

結論：
「均一な分布」という仮定だけで計算すると、「選択のフィルター（誰が見えるか）」を無視することになり、「星の明るさ（距離の基準）」が約 0.05 マグニチュード（明るさの単位）もずれてしまうことが分かりました。これは、ハッブル定数の矛盾を「解消した」のではなく、**「計算ミスで矛盾を小さく見せただけ」**だったのです。

5. 最終的な結論：矛盾は残ったが、信頼性は上がった

この論文の新しい計算方法（フィルターと銀河の形を正しく考慮した方法）で計算し直した結果：

• 従来の「矛盾（5σ）」は解消されませんでした。
• 計算結果は、従来の「SH0ES」というチームの値と非常に良く一致しました。

これはどういう意味？
「ハッブルの緊張」は、計算のミスではなく、「本当に新しい物理法則（未知の何か）」がある可能性が高いことを示唆しています。もし計算方法に問題があったなら、この新しい方法で直るはずですが、直らなかったからです。

まとめ

• 問題： 宇宙の膨張速度を測ると、2 つの違う答えが出てしまう。
• 原因の疑い： 「距離の測り方（星の選び方）」にミスがあったのではないか？
• この論文の貢献：
  • 「見えない星（フィルター）」を無視して計算すると、結果が狂うことを証明した。
  • 「均一な分布」という安易な仮定を使わず、**「実際に観測された星だけが選ばれた理由」**を完璧に計算に組み込んだ。
• 結果： 計算ミスはなかったことが確認され、**「宇宙の謎（ハッブルの緊張）は、本当に新しい物理学の発見を待っている」**という結論がより強固になった。

つまり、この論文は**「計算のやり方を完璧に整えて、本当に『新しい物理』を探す準備が整った」**と宣言したような、非常に重要な研究なのです。

技術概要

以下は、提示された論文「Forward-modelling Milky Way Cepheids: selection effects and physical priors in the Gaia–HST calibration」の詳細な技術的サマリーです。

1. 研究の背景と問題提起

ハッブル定数（$H_0$）の値を巡る「ハッブル・テンション（Hubble Tension）」は、宇宙マイクロ波背景放射（CMB）から推定される初期宇宙の値と、局所距離梯子（Distance Ladder）を用いた現在の宇宙の値との間に約 5$\sigma$ の不一致が存在する重大な問題です。この局所距離梯子の較正において、銀河系内のセファイド変光星（Cepheids）の距離梯子の第 1 段階（銀河系内較正）は極めて重要です。

近年、Gaia EDR3 による高精度の視差測定が可能になりましたが、従来の SH0ES チームによる解析（Riess et al. 2021, 2022）や、Högås & Mörtsell (2026; 以下 HM26) による新しいアプローチには、統計的な仮定に関する課題が残されています。特に、HM26 は銀河系セファイドに対して「体積一様（uniform-in-volume）」な距離事前分布を採用し、選択効果（selection effects）を明示的にモデル化しませんでした。HM26 はこのアプローチによりハッブル・テンションが緩和されると報告しましたが、本論文はこれが選択効果を無視した結果生じたアーティファクト（人工的な結果）であると主張しています。

核心的な問題:

• 銀河系セファイドのサンプルは、観測キャンペーン（HST Cycle 22 と 27）の選別基準（明るさ、視差、周期、赤化など）によって制限されており、母集団全体を代表する「体積一様」なサンプルではありません。
• 選択効果を無視して物理的に不適切な事前分布（例：体積一様分布）を適用すると、距離梯子の較正パラメータ（周期 - 光度関係のゼロ点や Gaia 視差のゼロ点オフセット）に重大なバイアスが生じます。

2. 手法：ベイズ前方モデル（Forward-Modelling）

本論文では、銀河系セファイドのサンプルを星ごとにモデル化する、完全なベイズ前方モデルを構築しました。このアプローチは、観測データを生成する過程を明示的にシミュレーションし、事後分布を推論します。

主要な構成要素:

1. データセット:

   • C22 (Cycle 22): 周期 $P > 8$ 日、距離 $\sim 6.6$ kpc までの遠方セファイド（44 星）。
   • C27 (Cycle 27): 近傍の高視差セファイド（22 星）。
   • LMC & N4258: 幾何学的距離（大マゼラン雲の連星、N4258 のメーザー）で較正された外部銀河のセファイド（それぞれ 69 星、443 星）。
   • 観測量：ウェーゼンハイト等級（$m_W^H$）、周期（$P$）、金属量（$[O/H]$）、Gaia EDR3 視差（$\varpi_{obs}$）。

2. 物理的事前分布（Physical Priors）:

   • 銀河円盤モデル: 銀河系セファイドの距離事前分布として、銀河の薄い円盤構造を反映したモデルを採用しました（半径方向と垂直方向のスケール長を持つ指数関数分布）。
   • HM26 との対比: HM26 が採用した「体積一様（$d^2$）」事前分布ではなく、銀河の幾何学的構造に即した事前分布を使用します。

3. 選択効果のモデル化（Selection Modelling）:

   • 観測サンプルがどのように選ばれたかを確率的にモデル化します。
   • 検出確率（Detection Probability）: 等級、視差、周期、赤化（$A_H$）のカット条件を、滑らかな閾値関数（誤差関数 $\Phi$）として表現し、積分計算によって検出確率を導出しました。
   • これにより、観測された星の分布と、モデルから生成される「検出された」サンプルの分布を一致させることが可能になります。

4. 推論フレームワーク:

   • パラメータ: 周期 - 光度関係（PL 関係）のパラメータ（ゼロ点 $M_W^H,1$、傾き $b_W$、金属量係数 $Z_W$）、Gaia 視差のゼロ点オフセット（$\delta\varpi$）、各キャンペーン固有の固有散乱（$\sigma_{int}$）、銀河系セファイドの個々の距離（$d_i$）などを同時に推論します。
   • アルゴリズム: JAX 言語と NUTS（No-U-Turn Sampler）を用いた MCMC 法で事後分布をサンプリングしました。
   • 解析的周辺化: 金属量や周期などの潜在変数、および選択関数の積分を解析的に行うことで、計算コストを削減し、高次元の積分を距離と天球位置への積分に還元しました。

3. 主要な成果と結果

1. PL 関係パラメータと Gaia オフセットの推定:

   • 本モデル（円盤事前分布＋選択効果モデル）による推定結果：
     • PL ゼロ点：$M_W^H,1 = -5.909 \pm 0.022$ mag
     • Gaia 視差オフセット：$\delta\varpi = -12.4 \pm 5.3$ $\mu$as
   • これらの値は、従来の SH0ES の最大尤度法による結果（R22）と $<0.7\sigma$ の範囲で一致しており、局所距離梯子の較正が統計的に堅牢であることを示しています。

2. 選択効果と事前分布の影響の定量化:

   • HM26 手法の再現とバイアスの特定: 選択効果を無視し、体積一様事前分布を採用した場合、PL ゼロ点は約 0.05 mag 明るく（$M_W^H,1 \approx -5.95$）、Gaia オフセットは約 -14 $\mu$as ほど負の方向にシフトしました。
   • 原因: このバイアスは、主に「選択効果の無視」によって引き起こされます。選択効果を正しくモデル化すれば、事前分布の形状（円盤型か体積一様か）の影響は微小（ゼロ点で 0.002 mag 未満）になります。
   • HM26 の結論の誤り: HM26 が報告したハッブル・テンションの緩和（$H_0$ の低下）は、不適切なモデル（選択効果の欠落）によるバイアスに起因しており、新しい物理的発見ではありません。

3. モデル検証（Posterior Predictive Checks）:

   • 本モデルから生成された模擬データ（パララックス、等級、周期の分布）は、観測データと非常に良く一致しました（KS 検定 p 値は 0.2 以上）。
   • 一方、選択効果を無視したモデルは、観測されたパララックス分布や等級分布を再現できず、明らかにデータと矛盾することが示されました。

4. 固有散乱（Intrinsic Scatter）:

   • C22 サンプル（遠方）の固有散乱は $\sim 0.07$ mag、C27（近傍）は $\sim 0.04$ mag と推定されました。C22 の大きな散乱は、遠方の視線方向における残存する赤化の不均一性による可能性が高いと結論付けられています。

4. ハッブル定数への影響

• 本論文で推定された PL ゼロ点は、SH0ES のベースライン値と統計的に有意な差がありません。
• このゼロ点を R22 のフル・距離梯子解析に適用しても、$H_0$ の値は $73.12 \pm 1.08$km s$^{-1}$Mpc$^{-1}$となり、従来の$73.04 \pm 1.01$ とほぼ変化しません。
• したがって、本論文の手法はハッブル・テンションを解消するものではなく、局所距離梯子の較正が既存の値と整合的であることを確認し、テンションが $\sim 5\sigma$ で残存することを再確認しました。

5. 意義と結論

• 統計的厳密性の確立: 距離梯子の解析において、選択効果（Selection Effects）を明示的にモデル化することの重要性を強く示しました。特に、Gaia 視差を用いた高精密な較正では、サンプルの不完全性を無視すると重大なバイアスが生じます。
• HM26 に対する反証: HM26 が報告したハッブル・テンションの緩和は、モデルの欠陥（選択効果の無視と不適切な事前分布）によるアーティファクトであることを実証しました。
• 将来への展望: この前方モデルの枠組みは、第 2 段階（セファイド宿主銀河）および第 3 段階（Ia 型超新星）の解析へ拡張可能です。将来、選択基準が明確に定義された専用観測キャンペーンに対して、この手法を適用することで、より信頼性の高い宇宙論パラメータの推定が可能になると期待されます。

要約すれば、本論文は「Gaia と HST データを用いた銀河系セファイドの較正において、選択効果を厳密に扱ったベイズ前方モデルを構築し、それが従来の結果と整合的であることを示した。一方で、選択効果を無視したモデル（HM26）はバイアスを含んでおり、ハッブル・テンションの緩和を主張することは誤りである」という重要な結論を導き出しています。