Bayesian approach study of hybrid neutron stars

技術要約：ハイブリッド中性子星のベイズ的アプローチによる研究

問題提起
中性子星の内部組成、特にハドロン物質から非閉じ込めクォーク物質への相転移の可能性は、天体物理学における未解決の問題である。量子色力学（QCD）は強い相互作用の基本理論であるが、高密度物質のab initio（第一原理）的な記述は現在不可能である。その結果、状態方程式（EOS）を記述するために有効モデルが必要となる。重要な課題は、ハイブリッド星（クォークコアを含む星）が存在し得るのか、もし存在する場合、どのような条件下で存在するのかを決定することである。本研究は、天文学的観測（具体的にはNICERによる質量および半径の測定、およびGW170817の重力波データ）が、クォーク物質モデルのパラメータおよび相転移の開始時期をどのように制約できるかに焦点を当てる。

手法
著者らは、ハイブリッド星モデルのパラメータを制約するためにベイズ推論フレームワークを採用している。この手法は以下の構成要素からなる：

1. 微視的モデル:

   • ハドロン物質: 相対論的平均場（RMF）近似を用いてモデル化されている。以下の3つのパラメータセットが利用された：
     • NL3*ωρ: 核子のみ（硬いEOS）。
     • EL3ωρ: 核子のみ（中間的な硬さ）。
     • EL3ωρY: ハイパーロン（ストレンジバリオン）を含む。
   • クォーク物質: ベクトル型MITバッグモデル（vMIT）を用いてモデル化されている。ラグランジアンには、反発的なベクトル相互作用、バッグ定数（$B$）、およびディラック海への寄与項（$b_4$）が含まれる。
   • 相転移: マクスウェル構成を用いて実装されている。これは、界面において圧力と化学ポテンシャルは連続であるが、バリオン密度は不連続に変化する一次相転移を仮定している。

2. 巨視的計算:

   • 星の構造（質量と半径）は、トルマン・オッペンハイマー・ヴォルコフ（TOV）方程式を解くことによって計算される。
   • 重力波の制約と比較するために、潮汐変形能（$\Lambda$）が計算される。
     物性値のクロスチェックのために、近似的な普遍関係（例：$C-\Lambda$, $\Upsilon-C$, $\Upsilon-\Lambda$）が用いられる。

3. ベイズ推論:

   • パラメータ ($\theta$): バッグ定数 ($B^{1/4}$)、ベクトル結合定数 ($G_v$)、およびディラック海の寄与 ($b_4$)。
   • 事前分布: $B^{1/4} \in [150, 180]$ MeV、$G_v \in [0.2, 0.3]$ fm$^{-2}$、$b_4 \in [0, 2]$ に対して一様分布を仮定した。
   • 尤度関数: 以下から構成される：
     • PSR J0030+0451、PSR J0740+6620、およびGW170817からの質量および半径の制約。
     • 大きなクォークコアを優先するための相転移密度制約（$n_b = 0.275$ fm$^{-3}$ を中心とする）。
     • 最大質量に対する厳格な下限（$M \ge 2.006 M_\odot$）。
   • 証拠（Evidence）の計算: dynesty サンプラーを使用してベイズ証拠を計算し、ベイズ因子を介してモデルを比較した。

主要な貢献および結果

• パラメータ制約: 本研究は、観測的制約を満たす特定のvMITパラメータセットを特定している。最も一貫した結果は、EL3ωρ パラメトリゼーション（核子のみ）において得られ、2つのパラメータ（$G_v, B^{1/4}$）を用い、約 $(0.22, 157.47)$ という値を与えた。

  • この構成は、転移密度 $n_b \approx 1.9 n_0$ および転移質量 $\approx 0.8 M_\odot$ を生成する。
  • また、無次元潮汐変形能 $\Lambda_{1.4} \approx 397$ を与え、これはpQCDの制約（$\Lambda_{1.4} < 400$）と一致している。
  • このモデルの自然対数証拠は $-4.54$ である。

• ハイパーロンおよび第3パラメータの影響:

  • ハイパーロン (EL3ωρY): ハイパーロンを含むモデルは、pQCDの制約に違反する転移密度（$\approx 2.13 n_0$）をもたらし、$\Lambda_{1.4} \approx 488$ という、pQCD限界と矛盾する値を示した。
  • 3パラメータモデル: ディラック海の寄与（$b_4$）を含めると、一般にpQCD制約に違反する結果（例：$\Lambda_{1.4} > 400$）をもたらすか、あるいは無次元慣性モーメントの普遍関係を満たさなくなった。
  • 硬いEOS (NL3*ωρ): より硬い NL3*ωρ パラメトリゼーションを用いたモデルも、転移質量が $1.0-1.2 M_\odot$ 付近となり、$\Lambda_{1.4}$ が 400 を超えるなど、pQCD制約を満たすことができなかった。

• ベイズ因子分析:

  • 2パラメータの EL3ωρ モデルを3パラメータ版と比較したところ、対数ベイズ因子は $0.92$ であり、より単純なモデルを支持する証拠は「決定的なものとは言えない（inconclusive）」ことを示した。
  • しかし、ベイズ因子を物理的制約（pQCDおよび普遍関係）と組み合わせると、2パラメータの EL3ωρ モデルのみがすべての基準に完全に一致する。

• 相転移の特性:

  • 本研究では、EL3ωρ および NL3*ωρ のパラメトリゼーションにおいて、密度 $2.0 n_0$ 未満で相転移が存在し得ることが示された。
  • 推定されたパラメータは、大きなクォークコアを持つハイブリッド星の存在を示唆している。

意義および主張
本論文は、ベイズ推論を高精度な天体物理学的データと組み合わせることで、非閉じ込めクォーク物質の性質を大幅に制約できると主張している。主な結論は、ハドロン物質が EL3ωρ パラメトリゼーションによって記述され、クォーク物質が2パラメータのベクトル型MITバッグモデルによって記述される場合、大きなクォークコアを持つハイブリッド星は、現在の質量、半径、および潮汐変形能のデータと整合する実行可能な解であるということである。

著者らは、ハイパーロンの導入や追加のベクトルチャネルパラメータ（ディラック海）の追加は理論的に動機付けられているものの、これらの拡張は、使用された特定の尤度関数の文脈においては、pQCD制約や普遍関係に違反する結果を招くと強調している。本研究は、「2パラメータ」のアプローチによるベクトル型MITバッグモデルが、EL3ωρ ハドロンEOSと組み合わされたとき、現在の観測状況においてハイブリッド星を最も堅牢に記述できることを示している。この研究は、実行可能なEOSモデルを精緻化し、ハドロン・クォーク相転移を支配するパラメータの不確実性を定量化する役割を果たしている。