🔭 astrophysics

The DESI DR1 Peculiar Velocity Survey: growth rate measurements from the maximum likelihood fields method

DESI DR1 のデータを用いた最大尤度法による解析により、構造成長率の測定値がΛCDM 宇宙論および一般相対性理論の予測と一致することが確認された。

原著者： Y. Lai, C. Howlett, J. Aguilar, S. Ahlen, A. J. Amsellem, J. Bautista, S. BenZvi, D. Bianchi, C. Blake, D. Brooks, A. Carr, T. Claybaugh, T. M. Davis, A. de la Macorra, P. Doel, K. Douglass, S. Ferrar

公開日 2026-04-13

📖 1 分で読めます☕ さくっと読める

CC BY 4.0

原著者： Y. Lai, C. Howlett, J. Aguilar, S. Ahlen, A. J. Amsellem, J. Bautista, S. BenZvi, D. Bianchi, C. Blake, D. Brooks, A. Carr, T. Claybaugh, T. M. Davis, A. de la Macorra, P. Doel, K. Douglass, S. Ferraro, A. Font-Ribera, J. E. Forero-Romero, E. Gaztañaga, G. Gutierrez, J. Guy, H. K. Herrera-Alcantar, D. Huterer, M. Ishak, R. Joyce, A. Kim, D. Kirkby, T. Kisner, A. Kremin, O. Lahav, C. Lamman, M. Landriau, L. Le Guillou, A. Leauthaud, M. E. Levi, M. Manera, P. Martini, A. Meisner, R. Miquel, J. Moustakas, A. Muñoz-Gutiérrez, S. Nadathur, W. J. Percival, C. Poppett, F. Prada, I. Pérez-Ràfols, F. Qin, C. Ross, G. Rossi, K. Said, E. Sanchez, D. Schlegel, M. Schubnell, H. Seo, J. Silber, D. Sprayberry, G. Tarlé, R. Turner, B. A. Weaver, P. Zarrouk, R. Zhou, H. Zou

原論文は CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) でライセンスされています。 ✨ これは以下の論文のAI生成解説です。著者が執筆または承認したものではありません。技術的な正確性については原論文を参照してください。免責事項の全文を読む

この論文は、**「宇宙の膨張スピードが、目に見えない『重力の力』によってどのように変化しているか」**を、新しい方法で詳しく調べた研究成果です。

専門用語を避け、身近な例え話を使って説明します。

1. 宇宙の「隠れた動き」を測る探検

私たちが夜空を見上げると、星や銀河は遠ざかっているように見えます（宇宙の膨張）。しかし、銀河にはもう一つ、**「見えない風」のような力が働いています。それが「特異速度（Peculiar Velocity）」**です。

例え話：
川（宇宙の膨張）が下流に向かって流れているとします。その川に浮かぶボート（銀河）は、川の流れに乗って下流へ進みます。しかし、ボートの中には、**「川の流れとは関係なく、自分自身で漕いだり、風で流されたりする動き」があります。
この「漕ぐ力」や「風の影響」こそが、銀河同士が重力で引き合い、集まろうとする動き（特異速度）です。この動きを測ることで、「重力が宇宙の構造をどう成長させているか」**がわかります。

2. 今回使った「巨大な網」と「新しい計算機」

この研究では、**DESI（ダークエネルギー分光装置）**という、アメリカの望遠鏡で集めた最新のデータを使いました。

データの規模：
約 8 万個の銀河の「距離と動き」のデータと、約 42 万個の銀河の「位置」のデータを組み合わせて分析しました。これはこれまでのどの研究よりも巨大なデータセットです。
新しい計算方法：
これだけの大量のデータを分析するのは、まるで**「砂漠の砂粒一つ一つを数えながら、その重さを計算する」ようなもので、普通の計算機では何週間もかかってしまいます。
そこで、この研究チームは「JAX（ジャックス）」**という、最新の AI 開発で使われる超高速な計算技術を使い、プログラムを書き直しました。
- 効果： 計算時間が「数週間」から**「1 日」、さらに評価ごとに「1 秒」**という驚異的な速さになりました。これにより、より正確で複雑な計算が可能になりました。

3. 2 つの「ものさし」を合わせて正確に測る

銀河の動きを測るには、2 つの異なる方法（ものさし）を使いました。

** Fundamental Plane（基本平面）：** 楕円銀河（形が丸い銀河）の明るさと大きさの関係から距離を測る方法。
** Tully-Fisher（タリー・フィッシャー）：** 渦巻銀河（形が渦巻きの銀河）の回転速度から距離を測る方法。

これら 2 つの「ものさし」と、銀河の「位置の分布（密度）」を同時に分析することで、「統計的な誤差（ノイズ）」を大幅に減らし、重力の強さをより正確に割り出そうとしました。

4. 結果：アインシュタインの理論は正しいか？

分析の結果、得られた数値は以下の通りです。

発見： 銀河が集まり、宇宙の構造が成長するスピードは、アインシュタインの一般相対性理論（現在の重力の標準モデル）が予測する値と、ほぼ一致していました。
意味： 「重力の法則」は、宇宙の広い範囲でも、これまで通り正しいことが確認されました。もし重力の法則が少し違っていたら（例えば、重力が弱かったり強かったりしたら）、この数値は大きくズレていたはずです。

5. まとめ

この論文は、**「最新の望遠鏡データ」と「AI 並みの超高速計算技術」を組み合わせることで、宇宙の重力の働きをこれまで以上に精密に測定し、「今の宇宙の重力の法則は、アインシュタインの予言通りである」**という結論を導き出した画期的な研究です。

まるで、**「宇宙という巨大なオーケストラの、見えない指揮者（重力）のテンポを、何千もの楽器（銀河）の音から正確に聞き分け、楽譜（物理法則）が正しいかを確認した」**ような作業でした。

DESI DR1 特異速度サーベイ：最大尤度場法による構造成長率の測定に関する技術的概要

本論文は、Dark Energy Spectroscopic Instrument (DESI) の第 1 回データリリース (DR1) における特異速度カタログを用いて、宇宙の大規模構造の成長率を測定した研究です。特に、銀河密度場と特異速度場を同時に最大尤度法 (Maximum Likelihood Fields method) で解析し、線形成長率パラメータ $f\sigma_8$ を制約することに焦点を当てています。

以下に、問題意識、手法、主要な貢献、結果、そして意義について詳細にまとめます。

1. 研究の背景と問題意識

宇宙論モデルの検証: 現在の標準宇宙モデル ( $\Lambda$ CDM) は、宇宙の加速膨張を説明するために宇宙定数を導入していますが、ハッブル定数の不一致（Hubble tension）やダークエネルギーの状態方程式の時間変化など、いくつかの矛盾が指摘されています。
重力理論の検証: 一般相対性理論 (GR) とは異なる重力理論（例：nDGP モデルや $f(R)$ 重力）は、異なる線形成長率を予測します。これらの理論を区別するためには、高精度な構造成長率の測定が不可欠です。
特異速度の重要性: 構造成長率を測定する手法として、赤方偏移空間歪み (RSD) と特異速度 (Peculiar Velocity, PV) があります。特異速度は、銀河の scaling relation（Tully-Fisher 関係や Fundamental Plane）を用いて赤方偏移に依存しない距離を測定することで得られます。特異速度は銀河密度場と相補的であり、両者を組み合わせることでサンプルバリアンスを低減し、統計的不確実性を削減できることが知られています。
DESI DR1 の特徴: 本研究で用いる DESI DR1 特異速度カタログは、約 8 万個の特異速度測定値を含み、過去最大の PV サンプルです。しかし、データ量が膨大であるため、従来の計算手法では解析に莫大な時間がかかるという課題がありました。

2. 手法と理論的枠組み

本研究では、銀河過密度 ( $\delta_g$ ) と対数距離比 ( $\eta$ ) の相関を直接フィットする最大尤度場法を採用しました。

2.1 データセット

特異速度カタログ: DESI DR1 の Fundamental Plane (楕円銀河、約 7.4 万個) と Tully-Fisher (渦巻銀河、約 6,800 個) のカタログを結合。
密度場サンプル: DESI DR1 の Bright Galaxy Sample (BGS) を使用（赤方偏移 $z < 0.1$ の範囲で約 41 万個）。
データ前処理: 外れ値（4 $\sigma$ 以上）の除去、および Tully-Fisher サンプルにおいて $z > 0.05$ の領域を除外（系統誤差の懸念のため）。最終的に 76,616 個の特異速度測定値と 415,523 個の銀河赤方偏移を使用。

2.2 理論的モデリングと計算効率の向上

共分散行列の解析的導出: 銀河過密度と特異速度の自己共分散および相互共分散を解析的に導出。非線形 RSD による減衰（Finger-of-God 効果）をガウスモデルで近似。
グリッド化 (Gridding): 計算コストを削減するため、データを $20 h^{-1}\text{Mpc}$ のグリッドに集約。これにより共分散行列の次元を大幅に削減し、ガウス尤度関数の適用を可能にしました。
ショットノイズ補正: 離散的な銀河分布によるショットノイズを共分散行列の対角成分に追加。
計算効率の革新 (JAX の採用): 従来のコードを Python の JAX ライブラリに書き換え、Just-In-Time (JIT) コンパイルと XLA による線形代数の高速化を実現。これにより、 $11,000 \times 11,000$ 程度の共分散行列に対する尤度計算を 1 回あたり約 1 秒に短縮し、MCMC 解析を数時間〜1 日以内に完了可能にしました。これにより、尤度関数のテイラー展開近似（従来手法）を不要にし、精度を向上させました。

2.3 モックカタログによる検証

AbacusSummit シミュレーションに基づいた 675 個のモックカタログを用いて、バイアスや系統誤差を検証。
結果、設定されたパラメータ（ $\sigma_u = 16 h^{-1}\text{Mpc}$ 、グリッドサイズ $20 h^{-1}\text{Mpc}$ ）において、推定された $f\sigma_8$ が真の値と整合的であり、手法が不偏であることを確認しました。

3. 主要な結果

成長率の測定値:
- 統計的不確実性と系統誤差を考慮した最終的な制約は以下の通りです。
  $f\sigma_8(z_{\text{eff}}=0.07) = 0.483^{+0.080}_{-0.043}(\text{stat}) \pm 0.018(\text{sys})$
- この値は、Planck + $\Lambda$ CDM モデルの予測値 ( $f\sigma_8 \approx 0.449$ ) と 1 $\sigma$ 以内で一致しています。
多 tracer 法の効果:
- 特異速度のみ、密度場のみ、および両者の組み合わせを比較。両者を組み合わせることで統計的不確実性が改善され、SDSS 特異速度カタログの結果よりも約 30% 精度が向上しました。
- 特異速度と密度場を組み合わせることで、 $f\sigma_8$ と非線形速度分散パラメータ ( $\sigma_u$ ) の間のデジェネラシー（重なり）を打破する効果がありました。
成長指数 ( $\gamma$ ) の制約:
- DESI の高赤方偏移データ (ShapeFit) と組み合わせることで、成長指数 $\gamma = 0.58 \pm 0.11$ を導出。これは一般相対性理論の予測 ( $\gamma \approx 0.55$ ) と一致しています。
他の DESI 手法との比較:
- 本研究の「最大尤度場法」は、DESI 内部で開発された「運動量パワースペクトル法」および「速度相関関数法」とも 1 $\sigma$ 以内で一致しており、3 つの異なるアプローチが整合的な結果を与えていることを示しました。

4. 主要な貢献と意義

計算効率の劇的改善: 大規模データ（数十万個の銀河）を用いた最大尤度法解析を、JAX と Cholesky 分解の導入により実用的な時間枠内で実行可能にしました。これにより、尤度近似なしの正確な解析が可能となり、将来のより大規模なデータセット（DR2 以降）への拡張基盤が整いました。
最大の PV サンプルによる制約: 既存の最大の単一特異速度カタログ（SDSS）と比較して、銀河数が約 10 倍、距離比測定値が約 2 倍のデータを用いることで、成長率の制約を約 30% 高精度化しました。
系統誤差の特定と対処: Tully-Fisher サンプルにおいて $z > 0.05$ の領域で Fundamental Plane サンプルと成長率の不一致が見られたため、これを保守的な測定を行うために除外しました。これは、高赤方偏移における Tully-Fisher 関係の未解明の系統誤差を示唆しており、今後の研究課題として明確化されました。
一般相対性理論の検証: 低赤方偏移領域において、GR の予測と整合的な成長率を得ており、現在の標準モデルを支持する結果となりました。

5. 結論

本研究は、DESI DR1 の特異速度データと銀河密度データを最大尤度場法で統合解析し、構造成長率 $f\sigma_8$ を高精度に測定しました。計算手法の革新により、大規模データ解析の精度と効率が飛躍的に向上し、一般相対性理論の予測と矛盾しない結果を得ました。今後は、MOPED などのデータ圧縮技術を用いて共分散行列の次元を削減し、非線形 RSD パラメータ $\sigma_u$ を自由パラメータとして扱うことで、さらなる系統誤差の低減と精度向上が期待されます。