🔭 astrophysics

Revisiting GW170817 at milliarcsecond scale: high-precision constraints on jet geometry and $H_0$

本論文は、GW170817 に関する利用可能なすべての VLBI データに対して、新しい流体力学モデルを用いた堅牢なベイズ解析を提示し、ジェット観測角度を精密に制限するとともに、ハッブル定数を $65.5\pm4.4$ km s $^{-1}$ Mpc $^{-1}$ と導出しており、この値は初期宇宙のプランク観測結果と後期宇宙の SH0ES 結果のいずれよりも前者に一致する。

原著者： Kelly Gourdji, Adam T. Deller, Chris Flynn, Taya Govreen-Segal, Cullan Howlett, Kunal P. Mooley, Ehud Nakar

公開日 2026-05-13

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

CC BY 4.0

原著者： Kelly Gourdji, Adam T. Deller, Chris Flynn, Taya Govreen-Segal, Cullan Howlett, Kunal P. Mooley, Ehud Nakar

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

宇宙を巨大な膨張する風船だと想像してみてください。何十年もの間、科学者たちは、その風船がどのくらいの速さで膨張しているのかについて議論を続けてきました。この膨張率を「ハッブル定数（ $H_0$ ）」と呼びます。問題は、この値を「初期宇宙」（ビッグバンの残光など）を用いて測定するとある数値が得られるのに対し、「後期宇宙」（近くの爆発する星など）を用いて測定すると、それとは異なりわずかに速い数値が得られることです。この不一致は「ハッブルのひずみ」として知られています。

2017 年、天文学者たちは画期的な発見を成し遂げました。彼らは、2 つの中性子星が衝突する「宇宙の衝突」を検知したのです。この「GW170817」と名付けられた事象は、重力波（時空のさざなみ）と光（明るい閃光）の両方を放出したため、特別でした。これにより、宇宙の膨張率を測定するのに役立つ「標準サイレン（宇宙の灯台）」となったのです。

しかし、課題がありました。その灯台までの距離を測定するには、灯台が私たちに対してどのように傾いているかを正確に知る必要があります。コマを横から見るのと、上から見るのとでは、見え方が異なります。傾きを誤って推測すれば、距離も誤って推測することになり、宇宙の膨張率の計算もずれてしまいます。

「盲点」の問題
GW170817 の場合、重力波検出器は傾きをあまり正確に特定できませんでした。なぜなら、この事象はほぼ正面から（銃の銃口をまっすぐ覗き込むように）観測されたからです。「左右の回転」と「上下の回転」を区別することが困難でした。

これを解決するために、科学者たちは衝突後に残る減衰する光、「残光」を観測する必要がありました。衝突の破片が宇宙空間へと急激に飛び出すと、粒子のジェットが形成されます。このジェットが天球上を時間とともに移動する様子を監視できれば、角度を特定することができます。

この論文の新たなアプローチ：高解像度マップ
ケリー・ガーディと共同研究者によって書かれたこの論文は、ぼやけた道路地図から高解像度の衛星画像へとアップグレードするようなものです。

従来の方法： 過去の研究では、時間の経過に伴う光の明るさ（「光度曲線」）を観察することで、ジェットの形状を推測しようとしました。しかし、これは車のエンジン音だけを聞いて車の形状を推測しようとするようなものです。エンジン音は車の形状だけでなく、その走行速度にも依存するため、これは厄介な問題です。
新しい方法： 著者たちは「VLBI（超長基線電波干渉法）」という技術を用いました。これは、地球全体にわたって電波望遠鏡を連結し、地球サイズの巨大な望遠鏡を1つ作り出すようなものです。これにより、彼らは「ミリ秒角」（地球から月にある硬貨を見るようなもの）の細部まで見えるほど鮮明な視界を得ることができました。

彼らは単に光の塊の「中心」を見るだけでなく、 sophisticated なコンピュータモデルを用いて、ガスジェットが宇宙空間に衝突した際にどのように振る舞うかをシミュレーションし、電波の形状と動き全体を直接分析しました。彼らは4 つの異なる時間点からのデータを組み合わせました。

0 日目： 衝突の正確な瞬間（ハッブル宇宙望遠鏡によって発見）。
75 日目および 230 日目： 望遠鏡ネットワークからの電波観測。
207 日目： 重要な「チェックポイント」として機能した高解像度の電波観測。

結果：より鮮明な角度
複雑な3次元モデルをこの超高精細なデータに適合させることで、彼らは以前よりもはるかに高い精度で観測角度を決定しました。

角度： ジェットは、私たちの視線から約18〜20 度傾いていることが判明しました（以前の推定値は 14〜28 度と、やや曖昧な範囲でした）。
距離： 角度をより正確に知ることができたため、距離もより正確に計算できました。衝突は約4,400 万光年の距離で起こったことがわかりました。
膨張率（ $H_0$ ）： この新しい正確な距離を用いて、ハッブル定数は65.5 km/s/Mpcと計算されました。

なぜこれが重要なのか
この結果は、宇宙論的な議論における「決着をつけるもの」です。

「初期宇宙」（プランク衛星）は、その率が67.4であると述べています。
「後期宇宙」（SH0ES チーム）は、その率が73.0であると述べています。
この新しい測定値は 65.5 です。

著者たちは、彼らの結果が「初期宇宙」の数値に非常に近い（ごく小さな誤差の範囲内）が、「後期宇宙」の数値よりも有意に低いことを指摘しています。これは、おそらく「後期宇宙」の測定値がわずかにずれているか、あるいは宇宙の膨張に関する私たちの理解の再考が必要であることを示唆しています。

留保事項
著者たちは、彼らの結果がいくつかの要素に依存していることを慎重に指摘しています。

ジェットモデル： 彼らは、ジェットがコンピュータシミュレーションと完全に同じように振る舞うと仮定しています。もし実際のジェットが奇妙であったり、衝突の残骸が予想よりも遅く輝くような隠れたエネルギー源を持っていたりする場合、角度の計算は変動する可能性があります。
局所運動： 銀河は宇宙の膨張に伴って動くだけでなく、独自の「固有運動」（局所的な重力による漂流）も持っています。著者たちはこれらの局所的な漂流を推定して平均化する必要があり、これにより不確実性が加わります。

まとめ
この論文は、回転するコマのぼやけた写真を、それがどのように傾いているかを正確に見るためのスーパーレンズで撮影するようなものです。それによって、彼らは衝突までの距離をより正確に測定し、宇宙が「ビッグバン」の測定値に「近くの星」の測定値よりも密接に一致する速度で膨張していることを発見しました。これはハッブルのひずみの謎を完全に解決するものではありませんが、「より遅い」膨張率に傾く、非常に強力で独立した証拠を提供するものです。

技術的サマリー：GW170817 のミリ秒角スケールでの再検討

問題提起
連星中性子星合体 GW170817 の検出は、ハッブル定数（ $H_0$ ）の「標準サイレン」測定における最初の機会を提供した。しかし、この測定の精度は、光度距離（ $D_L$ ）と連星系の傾斜角（ $\iota$ ）の間の強い共分散によって本質的に制限されている。重力波（GW）データのみでは、GW170817 のような正面を向いた事象の傾斜角を精密に制約することはできないが、合体によって放出された相対論的ジェットは、この幾何学に対する電磁気的な手がかりを提供する。残光の光度曲線を用いて観測角度（ $\theta_v$ ）を制約しようとした以前の試みは、観測角度とジェット開口角（ $\theta_{cp}$ ）の間の幾何学的な縮退によって制限されてきた。一方、ジェットを超光速固有運動として捉えた超長基線干渉計（VLBI）観測は、直交する制約を提供するが、歴史的には離散的な半解析的ジェットモデルの適合や、パラメータ空間の限られた部分集合への依存に留まり、完全な空間情報（例えば源の大きさ）や系統的な不確実性を厳密に扱うことが欠けていた。

手法
本研究は、GW170817 の観測幾何学、およびそれに続く $H_0$ に対する堅牢な制約を抽出するために設計された包括的なベイズ的視平面モデル適合フレームワークを提示する。手法は以下の構成要素を統合している：

データ統合: 分析には、利用可能な 4 つのミリ秒角スケールの測位データセットがすべて組み込まれている：
- Mooley ら (2018) による、2 つの高精度アレイ（HSA）VLBI エポック（合体後 75 日目と 230 日目）。
- Ghirlanda ら (2019) による、1 つのグローバル VLBI エポック（207 日目）。これは HSA 較正との整合性を確保するため、著者らが独立して再処理を行った。
- 合体のアンカー（0 日目）として機能する、キロノバの位置を測定したハッブル宇宙望遠鏡（HST）のデータ（8 日目）。
流体力学モデル化: 離散的な事前計算モデルに依存するのではなく、著者らは 2 次元相対論的シミュレーションに基づく新しい流体力学的ジェット残光モデルの連続体に対して適用されるスケーリング形式を採用している（Govreen-Segal & Nakar 2023）。これらのモデルは、4 つのジェットエネルギー分布（トップハット型および指数 $b=3, 4, 6$ のべき乗則）を網羅している。このスケーリングにより、モデルの特定のパラメータ比に対する不変性を利用し、任意の観測角度（ $\theta_v$ ）、開口角（ $\theta_{cp}$ ）、および観測者時間に対してモデルを調整することが可能となる。
ベイズ推論: フレームワークは、画像平面のデコンボリューションに伴う情報損失を回避し、GPU 加速ソフトウェア GALARIO とネストドサンプラー Dynesty（Bilby 経由）を用いて、観測された電波干渉計の可視度（visibilities）に対して直接動作する。モデルには、各エポックに対する測位オフセット（ $\delta_{ra}, \delta_{dec}$ ）およびフラックススケール不確実性（ $\delta_F$ ）を含む系統的な不確実性を厳密に処理するためのノイズパラメータが含まれている。
宇宙論的制約: $H_0$ を測定するために、フレームワークは GW 尤度（傾斜角と距離のカーネル密度推定量を介して）および固有速度補正を組み込んでいる。著者らは、銀河団カタログと固有速度再構成マップの 28 通りの組み合わせを検討し、これらの系統的な不確実性に対して周辺化するためにベイズモデル平均化（BMA）を採用している。

主要な貢献

連続パラメータ空間: 離散モデルの適合を超えて、ジェット幾何学（ $\theta_v$ および $\theta_{cp}$ ）の連続的な範囲を探索することで、特定のジェット構造の選択によって導入されるバイアスを回避している。
視平面適合: 画像平面の重心ではなく、可視度に対して直接適合することにより、縮退を解くために決定的な源の大きさの制約を含む、データの完全な空間情報を活用している。
系統的な不確実性の処理: 測位およびフラックス較正の不確実性に対するノイズパラメータの導入と、固有速度補正に対する BMA アプローチの併用は、固定された補正を仮定することが多かった以前の研究よりも堅牢な誤差予算を提供する。
独立したグローバル VLBI 再処理: HSA データとのフレーム整合性を確保するため、著者らはグローバル VLBI データを独立して再処理し、以前の再処理に対する測位オフセットを特定・修正した。

結果

観測角度の制約:
- 表面輝度変動の測定に基づき、光度距離を $D_L = 40.7$ Mpc に固定すると仮定した場合、ジェット幾何学に関する光度曲線に基づく事前分布を組み込むことで、観測角度は $\theta_v = 18.3^\circ - 20.3^\circ$ （68% 信頼区間）と得られる。
- $D_L$ および固有速度補正に対して周辺化した場合、好まれる観測角度はわずかに $\theta_v = 16.8^\circ - 19.2^\circ$ へシフトする。
光度距離: 距離の直接適合は $D_L = 44.0 \pm 1.6$ Mpc を与える。これは表面輝度変動法から導出された $40.7 \pm 2.4$ Mpc の値よりも約 $1\sigma$ 大きく、観測角度の制約が仮定された距離に敏感であることを浮き彫りにしている。
ハッブル定数（ $H_0$ ）: 28 の固有速度補正シナリオにわたるベイズモデル平均化された測定値は、 $H_0 = 65.5 \pm 4.4$ $H_{0} = 65.5 \pm 4.4$ km s $^{-1}$ $^{- 1}$ Mpc $^{-1}$ $^{- 1}$ である。
- この値は、初期宇宙のプランク協力団の値（ $67.4 \pm 0.5$ km s $^{-1}$ Mpc $^{-1}$ ）と $0.5\sigma$ の範囲内にある。
- 後期宇宙の SH0ES 測定値（ $73.04 \pm 1.04$ km s $^{-1}$ Mpc $^{-1}$ ）よりも $1.7\sigma$ 低い。

意義と主張
本論文は、この独立した後期宇宙の $H_0$ 測定が、SH0ES の結果よりも初期宇宙のプランク結果とより一致していると主張している。著者らは、これを表面的に受け取れば、ハッブル定数問題（Hubble tension）を解決するための非標準宇宙論モデルの必要性に疑問を投げかける可能性があると示唆している。しかし、著者らはこの主張に対して控えめであり、以下のいくつかの留保事項を明示的に指摘している：

モデルの誤指定: 残光光度曲線の後期進化（ $\gtrsim 350$ 日）は、標準的なシンクロトロンモデルではよく適合せず、過剰放射の源は不明である。著者らは、分析（より早期のエポックに焦点を当てたもの）がモデルの不正確さよりも幾何学によって支配されていると論じているが、未モデル化の物理が結果をバイアスする可能性を認めている。
ジェットのアライメント: 分析では、ジェットが軌道角運動量軸に沿って放出されると仮定している。もし有意なミスマッチが存在すれば、 $H_0$ の推定値に非線形的なバイアスを生じさせる可能性がある。
固有速度の支配性: 著者らは、観測角度の制約そのものではなく、固有速度補正が $H_0$ 測定における不確実性の支配的な源（約 5% の分数的不確実性に寄与）であると特定している。

本研究は、ジェット残光の高精度 VLBI 観測が近いうちに傾斜角を制約するための独自の強力なツールであり続けるが、ハッブル定数問題の決定的な解決を達成するには、より多くの明るいサイレンのサンプルと、ジェット物理および固有速度のモデル化の改善が必要であると結論付けている。

関連論文