Determination of neutron star radius from pulse profile modeling using profile likelihood

技術要約：プロファイル尤度を用いたパルス・プロファイル・モデリングによる中性子星半径の決定

問題提起
Neutron Star Interior Composition Explorer (NICER) は、中性子星の質量と半径を決定するために、主に X-ray Pulse Simulation and Inference (X-PSI) パッケージを用いた X 線パルス・プロファイル・モデリングを利用している。現在の標準的な手法は、ネステッド・サンプリング・アルゴリズム（例：MultiNest）を用いたベイズ推論に依存しており、事後分布を推定している。しかし、半径が 11 km を超える場合、平坦な尤度曲面においてサンプラーの選択や収束の問題に起因する可能性がある、推定半径の不一致が異なる研究間で観察されている。さらに、広い事前分布範囲に対するベイズ的周辺化は、パラメータ推定を偏らせるボリューム効果を導入する可能性がある。ベイズ手法は天体物理学における標準であるが、頻度論的推論は素粒子物理学において依然として支配的である。本研究は、X-PSI から導出される基礎物理学の制約の堅牢性を確保するために、ベイズ的アプローチと頻度論的アプローチを比較する必要性に対処するものである。

手法
著者らは、基礎となる物理モデルおよび尤度実装を標準的なベイズ的アプローチと同一に保ったまま、既存の X-PSI フレームワーク内にプロファイル尤度技術を用いた頻度論的推論パイプラインを実装した。

• データソース: 解析には、H25 においてネステッド・サンプラーをテストするために使用された合成 NICER + XMM-Newton データセット（syntX1 と指定）を用いた。
• モデル構成: 本研究では、2 つの単一温度かつ単純連結な円形ホットスポットを特徴とする ST-U パルス・プロファイル・モデルを採用した。中性子星の大気は NSX 完全電離水素モデルを用いてモデリングされ、光子の伝播は、扁平なシュヴァルツシルト時空における一般相対論的レイ・トレーシングを通じて計算された。
• 頻度論的実装:
  • パラメータ・ベクトルは、関心のあるパラメータ（赤道半径、$R_{eq}$）と、妨害パラメータ（$\alpha$）に分割された。
  • 妨害パラメータ（質量、距離、傾斜角、ホットスポットの幾何学、温度、および星間水素柱密度を含む）は、固定された $R_{eq}$ に対して $L(R_{eq}, \alpha)$ を $\alpha$ に関して最大化することにより、プロファイル化された。
  • 最適化には、Nelder-Mead 単純形アルゴリズム（scipy.optimize.minimize）を用いた。
  • 信頼区間は、1 自由度の $\chi^2$ 分布を仮定し、$\Delta\chi^2$ 統計量（$\Delta\chi^2 = -2 \ln L_{profile} + 2 \ln L_{max}$）を用いて構築された。
• 比較ベースライン: 結果は、同一の設定および低解像度の X-PSI 設定を用いた MultiNest ネステッド・サンプリング・アルゴリズムによる標準的なベイズ解析と比較された。

主な貢献

1. 原理実証的適用: 本研究は、X-PSI パッケージを用いて中性子星半径を決定するための、頻度論的プロファイル尤度推論の初の実装例を示すものである。
2. 手法の比較: 同一の合成データに対して、ベイズ的推論と頻度論的推論のパイプラインを直接比較し、推論手法のみを変数として分離した。
3. オープンサイエンス: 著者らは、さらなる研究を促進するために、X-PSI 用の頻度論的推論コードの修正版を公開した。

結果
注入された真の半径が 12.176 km である合成データセット syntX1 を用いた結果：

• ベイズ推論: 半径の周辺化された 68% 信用区間は $R_{68,bayes} = 11.970^{+0.180}_{-0.180}$ km であった。この結果は、注入値と $1.1\sigma$ 以内で一致した。
• 頻度論的推論: プロファイル尤度は 12.096 km の最良適合半径を与えた。68% 信頼区間は $r_{68} = 12.096^{+0.064}_{-0.182}$ km であり、95% 信頼区間は $r_{95} = 12.096^{+0.316}_{-0.397}$ km であった。頻度論的推定は、注入値と $0.4\sigma$ 以内で一致した。
• 精度: 頻度論的手法の精度は、ベイズ解析と同等であることが示された。
• 計算効率: 頻度論的解析には約 9 CPU コア時間が必要であったのに対し、ベイズ解析には約 3,840 CPU コア時間（48 コアのノードで 80 時間）を要した。したがって、頻度論的アプローチは数桁高速である。

意義および主張
本論文は、本研究を**原理実証（proof-of-principle）**研究として位置づけている。著者らは、頻度論的プロファイル尤度推論を用いることで、合成データを用いて真の中性子星半径を $1\sigma$ 未満の精度で正常に回収できると主張している。この手法は、ベイズ解析と同等の精度を提供しつつ、計算コストを大幅に削減する。著者らは、この手法が X-PSI におけるパラメータ推定のための堅牢な代替手段となり、ベイズ推論を補完する価値のある手法であると結論付けている。彼らは、今後の課題としてこの技術を NICER 望遠鏡からの実際のパルサー・データに適用することを明言しているが、本稿ではそのような適用は提示されていない。