Incorporating neutron star physics into gravitational wave inference with physics-informed priors using normalizing flows

中性子星やブラックホールのような高密度の恒星残骸が互いに螺旋を描きながら接近し、衝突すると、時空の織物に「重力波」と呼ばれるさざ波を送ります。地球上の検出器はこれらの微かな信号を捉え、科学者が、わずか数十年前までは不可能であった方法で宇宙の「音」を聴くことを可能にしています。衝突の間に何が起こったのかを理解するために、研究者たちは、波の音からそれらを生み出した天体の特性へと、統計的な手法を用いて逆算します。このプロセスには、信号が届く前から、それらの天体がどのような性質を持つ可能性があるかについて、教育的な推測を行うことが求められます。伝統的に、科学者は「アグノスティック（無知）」な推測を使用してきました。これは、天体の質量や内部構造に関する事前知識を持たないことを前提とし、あらゆる可能性を等しく起こり得るものとして扱うものです。しかし、このアプローチは、極限の圧力下における物質の振る舞いや、銀河内にどのように中性子星が分布しているかといった、既存の膨大な知識を無視してしまいます。

新しい研究は、こうした推測を扱うためのよりスマートな方法を提案しています。研究者たちは、現実世界の物理学を重力波信号の解析に直接組み込む手法を開発しました。中性子星の内部組成や、既知の天体の典型的な質量に関する既存の知識を無視する代わりに、彼らは最初からこの情報を組み込んだシステムを構築しました。核物理学の法則や既知の恒星の観測に基づいた人工知能の一種を用いることで、彼らは「物理学に基づいた（physics-informed）」期待値を作り出しました。この手法を3つの特定の宇宙衝突から得られた実際のデータに適用したところ、衝突した天体の性質を裏付けるだけでなく、従来のメソッドよりもはるかに鮮明な特性の姿を描き出すことに成功しました。

この研究の核心は、これら衝突から得られるデータをどのように解釈するかという点にあります。重力波信号が検出されるとき、それは天体の質量、自転、そして接近するにつれて内部構造がどのように変形するかといった情報の複雑な混合物です。この情報を抽出するために、研究者はベイズ推論、すなわち、新しい証拠が得られるたびに仮説の確率を更新していく数学的枠組みに依拠しています。この更新の出発点は「事前分布（prior）」であり、これは新しいデータを見る前に、科学者が真実であると信じている事柄を表します。長年、標準的な慣習は、バイアスを避けるためにはそれが最も安全であると考え、あらゆる質量や内部構造が可能であると仮定する「無情報事前分布」を使用することでした。しかし、私たちは中性子星についてすでに多くのことを知っています。彼らが極めて高密度な物質で構成されており、特定の物理法則に従っていることも知っていますし、私たちの銀河系内にある多くの中性子星の質量も測定しています。新しい研究は、この知識を捨て去ることは機会損失であると主張しています。

これを解決するために、チームは中性子星の物理学に関する現在の理解を、解析に直接エンコードする柔軟なシステムを作成しました。彼らはまず、「状態方程式」に対する制約を集めました。これは、本質的に、中性子星の内部に見られる極限の密度下で物質がどのように振る舞うかを記述するルールブックです。このルールブックは、核理論の実験や、高速回転するパルサー（中性子星の一種）の観測に基づいています。また、彼らは私たちの銀河における放射観測から、高質量の中性子星が通常どの程度の質量であるかというデータも収集しました。これらの情報を用いて、彼らは、もし衝突がこれらの物理法則に従う天体によって引き起こされたとしたら、どのような衝突が見えるかというシナリオを数百万通りシミュレーションしました。

これらのシミュレーションは、「ノーマライジング・フロー（正規化流）」として知られる機械学習モデルの訓練に使用されました。このモデルを、最も可能性の高い物理的シナリオの形状を学習する、非常に効率的な「地図製作者」と考えてください。一度訓練されると、このモデルは洗練されたフィルターとして機能します。解析が質量や内部構造のあらゆる理論的な組み合わせの間を彷脱することを許すのではなく、物理学が天体が最も存在しやすいと示唆する領域へと探索を導きます。これにより、解析は最も妥当な説明に注意を集中させることができ、結果をより精密にし、結論をより強固なものにします。

研究者たちは、アメリカ、イタリア、日本にある観測所によって検出された3つの実際の重力波イベントに対して、この新しいアプローチをテストしました。最初のイベント、GW170817は、光による観測でも捉えられた二つの中性子星の衝突でした。第二のイベント、GW190425は、別の重い成分を持つ中性子星の衝突でした。第三のイベント、GW230529は、中性子星とブラックホールの衝突であり、以前は確実に分類することが困難であったイベントです。あらゆるケースにおいて、物理学に基づいた手法は、従来のガイドなしのアプローチよりも明確な宇宙の姿を提供しました。

最初のイベントであるGW170817について、新手法は、衝突した天体が確かに二つの中性子星であることを確認しました。さらに重要なことに、この手法は、これらの中の「内部の硬さ」に関する異なる理論を区別することを可能にしました。データは「ソフト（軟らかい）」な状態方程式、つまり、星の内部の物質が比較的圧縮しやすいことを強く支持していました。この結論は、重力波データと光の観測を組み合わせた他の研究の結果とも一致していましたが、新しい手法は、重力波信号のみを用いてこの結論に到達しました。また、衝突がどの程度の距離で発生したかについても、より精密な測定を行い、推定距離を以前の解析よりもわずかに遠くに押し上げました。

第二のイベント、GW190425は、衝突した星々が私たちの銀河で見られる典型的な中性子星よりもはるかに重かったため、難題を提示しました。伝統的な解析では、彼らの正確な性質を特定することに苦慮しました。しかし、物理学に基づいたアプローチは、その異常な質量にもかかわらず、彼らを連星中性子星システムであるとはっきりと特定しました。また、この手法は、二つの星が同じ重さではない可能性が高いことも明らかにしましたが、これは従来のガイドなしの手法が見逃していた詳細です。解析は、このイベントに対して、信号が比較的微弱であった事実とより整合する、以前考えられていたものよりも大幅に大きい距離を示唆しました。

第三のイベント、GW230529は、中性子星とブラックホールによるものでした。信号が微弱であったため、天体が正確に何であるか、あるいはどのように回転しているかを判断することが困難でした。新しい手法は、可能性を絞り込むことで、ここで優れた能力を発揮しました。この手法は、このイベントを中性子星とブラックホールの衝突であると断定し、他の可能性を高い信頼度で排除しました。さらに、物理的な制約を取り入れることで、二つの天体の間の質量比について、より厳密な推定値を提供しました。これは結果として、以前は不確実であったブラックホールのスピンのより明確な姿を与えました。また、この手法は、信号の微弱さに矛盾しないよう、このイベントの推定距離をさらに遠くに押し上げました。

これら3つのイベントすべてに共通する重要な発見は、物理学に基づいた事前分布が、一貫して衝突までの距離の推定値を大きくしたことです。これは、新しい手法が天体の質量をより厳密に制約したために起こりました。信号の大きさは質量と距離の両方に依存するため、質量をより正確に知ることは、距離の計算を調整させることになります。これらのケースでは、その調整は、天体が以前考えられていたよりも遠方にあることを意味していました。この変化は、これらの衝突に伴う光を探している天文学者にとって現実的な影響を及ぼします。なぜなら、より遠方の天体はより暗く、発見が困難になるからです。

この研究はまた、このアプローチが曖昧なイベントを整理するための強力なツールになり得ることを示しました。GW230529の場合、この手法は、衝突の性質について決定的な証拠を提供し、これは以前は微弱な信号では達成が困難であったことです。研究者たちは、物理学の法則を用いて探索を導くことで、より大きな信号やより優れた検出器を待つことなく、同じデータからより多くの情報を抽出できることを示しました。

将来に向けて、研究者たちは、彼らの手法が柔軟に設計されていることを強調しています。中性子星に関する理解が進んだり、新しいタイプの観測が可能になったりすれば、システムを更新することができます。これは解析プロセス全体の完全な刷新を必要とするものではなく、単にコンピュータがデータをナビゲートするために使用する「地図」を更新するだけです。つまり、私たちが高密度な物質の内部についてより多くを学ぶにつれて、宇宙の音を解釈する私たちの能力も自動的に向上していくのです。この研究は、データを白紙の状態として扱うのではなく、既知の自然界の法則との「対話」として扱うことへの転換を意味しており、これにより、私たちはより鮮明に宇宙の声を聴くことができるようになるのです。