Generalized Non-linear Bayesian Pulsar Timing with Enterprise

宇宙の静寂の奥深く、暗闇の中で灯台のように回転しているのは、その物質のわずかスプーン1杯分が地球上で10億トンの重さになるほど高密度な中性子星である。パルサーとして知られるこれらの星の中には、最高の研究所にある原子時計にも匹敵する精度で回転し、電波のビームを時計仕掛けのような規則正しさで地球へと掃引するものがある。天文学者がこれらの信号を捉えるとき、彼らはパルスの到着時刻をマイクロ秒の端数まで正確に測定することができる。この驚異的な精度により、これらの遠方の星は宇宙の実験室へと姿を変え、科学者が極限の環境下における重力の振る舞いを含む、物理学の基本法則を検証することを可能にしている。しかし、これらの星を道具として活用するためには、まずその回転から伴星の重力に至るまで、あらゆる要素を考慮した完璧な運動モデルを構築しなければならない。もしモデルがわずかにずれていたり、データに隠れたノイズが含まれていたりすると、星の特性の測定値がドリフト（漂流）し、物質の本質について誤った結論を導き出すことになりかねない。

研究チームは現在、数十年にわたり使用されてきた伝統的な手法を超え、これらのモデルを構築するための新しい方法を開発した。パルサーの動きを単純化された直線近似に頼るのではなく、タイミングデータとその中のノイズを単一の相互接続されたシステムとして扱う、強力で柔軟なアプローチを用いたのである。この手法を4つの特定のパルサーに適用した結果、彼らの新技術は、従来の研究よりも、星の質量や軌道に対してより厳密で正確な制約を与えることが多いことが判明した。いくつかのケースでは、星が以前考えられていたよりもわずかに異なる回転をしていることを発見し、さらに初めて、あるパルサーにおいてデータの中に隠れていた特定の種類の内部ノイズの明確な証拠を見出した。この研究は、コンピュータに単純な形状を強制するのではなく、データの持つ全容的な複雑さを探索させることで、宇宙で最も重い天体についてより多くのことを学べることを示唆している。

研究者たちは、中性子星が伴天体（多くの場合、白色矮星）の周囲を回っている連星パルサー4つに焦点を当てた。これらの星は重力のダンスに縛り付けられているため、そのタイミング信号には軌道運動の刻印が刻まれている。パルスの到着方法を分析することで、科学者はパルサーと伴天体の質量を測定できる。これは、中性子星の質量が、地球上では再現不可能な密度の状態における物質の状態を教えてくれるため、極めて重要である。もし中性子星が重すぎればブラックホールへと崩壊する可能性があり、もし十分に軽ければ、核力がそのような圧力の下でどのように維持されているかを明らかにする。長年、科学者はこれらの質量を決定するために、一般化最小二乗法と呼ばれる手法を用いてきた。この手法は、観測されたパルス時刻と予測された時刻との差を最小化する最適なパラメータのセットを見つけ出すものである。しかし、この手法は、残差（エラー）がいかなるものであれ、小さく線形であり、主要なモデルに対する単純な補正として扱えるという仮定に基づいている。

アンドリュー・カイザー氏らによるこの新しい研究は、Enterpriseと呼ばれるソフトウェアパッケージを用いたベイズ的アプローチを用いることで、この仮定に挑戦している。この枠組みにおいて、研究者たちは単に「最良の適合」を見つけるだけでなく、与えられたデータに対して各解がどの程度尤もらしいかを計量しながら、起こりうるすべての解の景観（ランドスケープ）を探索する。彼らは3つの異なるタイミングモデルの扱い方をテストした。第一の手法は、伝統的なアプローチを模したもので、タイミングパラメータを、ノイズを分析する前に数学的にデータから取り除かれる線形補正として扱う。第二の手法は、最も計算負荷が高いもので、タイミングモデルを完全な一般非線形信号として扱い、分析の各ステップごとにモデル全体を再計算する。第三の手法は、これら2つを組み合わせたハイブリッドである。これらのアプローチを比較することで、チームは、伝統的な線形仮定が重要な詳細を隠してしまっていないか、あるいは結果を歪めていないかを確認しようとした。

彼らがパルサーPSR J2043+1711にこれらの手法を適用した際、予想外の結果が得られた。これまでの研究では、固有のレッドノイズ（パルサーの自転における緩やかでランダムな変動の一種）の有意な証拠は見つかっていなかった。しかし、彼らの完全な一般ベイズ法を用いることで、チームはこのノイズを初めて検出した。彼らは、パルサーの自転が単に完璧に安定しているだけでなく、微細な低周波のジッター（揺らぎ）を含んでいることを発見した。この発見は、パルサーのタイミングパラメータのイメージを一変させ、その回転速度と伴天体の質量の推定値を変化させた。この研究は、もしこのノイズが無視されたり、簡略化されたモデルで扱われたりした場合、タイミングパラメータがそれを補うようにシフトし、星の物理的特性についてわずかに異なる結論を導いてしまうことを示した。

パルサーPSR J1600–3053については、新しい手法が星の質量に対してより厳密な制約を与えたことをチームは見出した。以前の測定では、大きな誤差範囲を伴って約2.0太陽質量の質量が示唆されていた。今回の新しい分析によって、この範囲は大幅に絞り込まれ、より精密な推定値が得られた。同様に、極めて高い質量を持つことで知られるパルサーPSR J0740+6620についても、チームの結果は従来の知見と一致していたものの、より鋭い制約を与えた。彼らは、この星が太陽の約2.06倍の質量を持つことを、非常に小さな誤差範囲とともに確認した。この精度は、中性子星の内部で物質が何を起こしているかについての特定の理論を排除するのに不可欠である。もし質量が高すぎれば、現在の理解における物理法則を破ることになり、強い核力の理論を修正する必要があることを示唆することになるからである。

最も劇的な結果は、パルサーPSR J1640+2224からもたらされた。以前の研究では、このパルサーは驚くほど高い質量推定値を示しており、時には太陽の3倍を超える値を示していたが、これは中性子星としては物理的に不可能である。研究者たちは、この高い質量は星の真の性質ではなく、モデリング手法によるアーティファクト（偽の産物）ではないかと疑った。これを検証するため、彼らはモデルに物理的な制限を課し、質量を太陽の3倍未満に制限した。これを行ったところ、推定質量は約2.2太陽質量という、より妥当な値へと低下した。この結果は、極めて重要な点を示している。すなわち、数学的モデルがタイミングパラメータとノイズの間の複雑な相互作用を説明できるほど柔軟でない場合、物理的に不可能な結果を生み出す可能性があるということである。モデルに全範囲の可能性を探索させ、その上で物理的制約を適用することで、チームはより現実的な星の姿を取り戻すことができた。

また、この研究は、電波望遠鏡によって収集される膨大なデータをどのように扱うかという問題にも取り組んだ。パルサー・タイミング・アレイは、多くの望遠鏡から長年にわたって収集されたデータを集めるが、このデータはしばしば、星と星の間にある星間物質からのノイズによって汚染されている。研究者たちは、彼らの完全な一般手法が、パルサーの運動による真の信号と、環境からのノイズを区別することに優れていることを見出した。多くの場合、伝統的な線形手法は、一部のノイズをタイミングパラメータの中に吸収してしまい、星の運動を実際とはわずかに異なるものに見せていた。タイミングモデルとノイズを結合されたシステムとして扱うことで、新しい手法はそれらをより効果的に分離することができた。

この研究の主要な教訓の一つは、統計的手法の選択が重要であるということである。研究者たちは、伝統的な線形アプローチは多くのパルサーに対してうまく機能するものの、データが複雑であったり、ノイズが顕著であったりする場合には失敗し得ることを発見した。完全な一般ベイズ法は、計算コストは高いものの、データを分析するためのより堅牢な方法を提供する。それは、科学者が起こりうるすべての解の範囲を見渡し、異なるパラメータが互いにどのように関連しているかを理解することを可能にする。これは、多くのパルサーの精密なタイミングに依存している将来の重力波探索にとって特に重要である。もしタイミングモデルが正確でなければ、重力波によって引き起こされる時空の微細なさざ波を見逃したり、誤解したりする可能性がある。

チームはまた、異なるデータセットが結果にどのように影響するかについても調査した。彼らは、10年以上にわたる観測を含むNANOGravコラボレーションのデータを分析した。彼らは、データが蓄積されるにつれて、パルサー質量の制約はより厳密になる一方で、中心値が時としてシフトすることを発見した。これは、タイミングモデルがいまだ進化過程にあり、初期の研究による最良適合値が最終的な答えではない可能性があることを示唆している。研究者たちは、さらなるデータが集まることで、モデルが安定し、測定値が異なるデータセット間で一貫していくことを期待している。

結局のところ、この論文は、宇宙の声を聞くための私たちの道具を洗練させることについての物語である。単純化された仮定から離れ、データの持つ全容的な複雑さを受け入れることで、研究者たちは、私たちが宇宙で最も極端な天体についてより多くのことを学べることを示した。彼らの研究は、パルサー・タイミングがどのように行われるべきかという新しい標準を提供し、中性子星の質量と軌道の測定が可能な限り正確であることを保証するものである。これは単なる技術的な改善ではない。それは、物質と重力の根本的な性質を理解するためのステップである。私たちがこれらの宇宙の時計に耳を傾せ続けるとき、この研究で開発された手法は、単一の中性子星の回転から、数十億光年離れたブラックホールの衝突に至るまで、宇宙の最も微かな囁きを聞き取る助けとなるだろう。