Ruling out nucleonic direct Urca cooling in low-mass neutron stars using nuclear data

技術要約：低質量中性子星における核子的直接Urca冷却の強い否定

問題提起
高密度物質の状態方程式（EoS）と組成を決定することは、核物理学における中心的な課題である。天体観測、核理論、および実験室での実験は補完的な制約を与えるものであるが、それらが探査する密度やアイソスピン非対称性の領域はそれぞれ異なっている。重大な困難は、核飽和密度付近に限定されている実験室の情報から、中性子星のコアに見られるような超飽和密度への外挿を行うことにある。従来のエネルギー密度汎関数は、飽和密度と超飽和密度の間に人工的な相関を課すことが多く、これがモデル依存の予測を招いている。具体的には、直接Urca（dUrca）冷却過程（$n \to p + e + \bar{\nu}_e$）の低質量中性子星における発現は、対称エネルギーによって支配される陽子分率に極めて敏感であるが、現在の制約では、質量 $\le 1.4 M_\odot$ の星でdUrca冷却が活性化するかどうかについて大きな不確実性が残っている。

手法
著者らは、多様なデータソースを単一の統計的に一貫したモデルである「漸近的因果性を持つメタモデル（asymptotically causal metamodel）」内で統合するために、ベイズ推論フレームワークを採用している。このアプローチは、高密度における因果性を保証しつつ、核物質パラメータ（NMP）の独立した変化を許容することで、硬直的な汎関数形式の限界に対処するものである。

本研究では、以下の3つの主要な制約を統合している：

1. 核構造データ（NSD）： 広範な静的および動的な核観測量（結合エネルギー、半径、双極子分極率など）から導出された、相関を持つNMPの完全な多次元事後分布（[19]における解析に基づく）。これは、低次NMPの完全な相関構造を捉えている。
2. INDRA-FAZIA (IF) データ： 重イオン衝突におけるアイソスピン輸送比の測定値。以前の戦略のように「許容バンド」を定義するのではなく、著者らは各メタモデルの実現における対称エネルギー曲線（$S_\theta(n)$）に沿って、IF尤度を直接評価している。これにより、核飽和密度を中心とした密度依存の制約が得られる。
3. カイラル有効場理論（$\chi$）および天体物理学： 純粋中性子物質に対するカイラルEFTの事前分布、および高質量パルサーの測定値（PSR J0348+0432, PSR J0740+6620）、潮汐変形能（GW170817）、およびNICERの質量-半径事後分布を含む標準的な天体物理学的尤度を組み込んでいる。

グローバルな事後分布は、$\chi$EFT、IF、およびNSDの事前分布と尤度の積を用いて、メタモデルパラメータ $\theta = (\text{NMPs}, X)$ をサンプリングすることによって構築される。得られたEoSは、トルマン・オッペンハイマー・ヴォルコフ（TOV）方程式を解くために使用され、中性子星の観測量を生成する。

主な貢献

• 直接的な尤度評価： 本論文は、密度依存のINDRA-FAZIA尤度を、事前に定義されたバンドに頼るのではなく、対称エネルギー曲線上で直接評価する新しい戦略を実装している。これにより、中間的な近似を回避し、核データの相関構造を保持している。
• 相関の伝播： 相関を持つNMPの多次元事後分布を核構造観測量から中性子星の性質へと伝播させることで、パラメータ間の統計的な依存関係が全過程を通じて維持されることを保証している。
• 柔軟な外挿： 追加の高密度自由度（$X$）を持つ漸近的因果的メタモデルを利用することで、より制限的な汎関数形式において問題となる、実験的に制約された領域（亜飽和密度）と外挿された領域（超飽和密度）の間の人工的な相関を軽減している。

結果
実験室の制約（NSDおよびIF）をカイラルおよび天体物理学的データと組み合わせることで、中性子星の組成、特に低質量星に関する予測が大幅に精緻化された：

• 対称エネルギー： 統合された制約（Full Nuclear, またはFN）は、これら3つのデータソースすべてと整合する対称エネルギーの事後分布を与える。IFデータは、対称エネルギーの値（$E_{sym}$）をその勾配（$L_{sym}$）から区別するのに役立ち、核構造データ単独の場合と比較して $E_{sym}$ に対する制約を強めている。
• 半径の制約： ラボの制約は、低質量中性子星に対して顕著な影響を与える。$2.0 M_\odot$ における半径事後分布の68%信頼区間の幅の減少がわずか13%であるのに対し、IFおよびNSDの制約を加えると、$1.0 M_\odot$ で約40%、$1.4 M_\odot$ で29%縮小する。
• Direct Urcaの発現： 最も重要な知見は、核子的dUrca過程が運動学的に許容されるために必要な最小質量（$M_{dU}$）に関するものである。
  • カイラルおよび天体物理学的制約のみの場合（$\chi$+astro）、$M \le 1.4 M_\odot$ の星でdUrcaが活性である確率は約17.5%である。
  • ラボの情報（IFおよびNSD）を含めた場合（FN+astro）、この確率は 0.26% まで低下する。
  • これは、約67倍の抑制因子を表している。中間的なアンサンブル（IF+$\chi$+astro）では1.2%であり、重イオン衝突からの対称エネルギーの情報がこの抑制の主要な要因であることを示している。
  • したがって、核子的仮説の範囲内では、低質量中性子星におけるdUrca冷却の発現は強く否定される。事後確率は約20%から1%未満へと低下したのである。

意義と限界
本論文は、核飽和付近の物質を探索する実験室実験が、中性子星の組成や高速冷却チャネルの運動学的アクセシビリティに対して定量的な制約を与えられることを示しているが、これらの制約は質量が大きすぎない天体に対して最も効果的である。本研究は、アイソスピン輸送データと核構造観測量の包含が、低質量中性子星における中心部の陽子分率の決定を大幅に改善することを確立した。

著者らはいくつかの注意点を挙げている：

• 結果は、核エネルギーのアイソスピン依存性に関する標準的な放物線近似に基づいている。これからの逸脱は、対称エネルギーと組成の間の接続を変化させる可能性がある。
• dUrca冷却の除外は、あくまで「核子的」仮説に適用されるものである。本論文は、ハイペロンによるdUrca過程による高速冷却を排除しているわけではない。ハイペロンの閾値は依然として制約が乏しく、モデル依存性が高い。
• 高密度への外挿には依然としてモデル依存性が残るが、メタモデルの柔軟性がこれを緩和している。将来の高エネルギー核衝突実験や次世代重力波干渉計からの制約は、特に高質量星に関して、これらの分布をさらに絞り込む可能性がある。

結論として、本研究は多様なデータソースを統合するための堅牢な統計的枠組みを提供し、現在の実験室データが低質量中性子星における核子的直接Urca過程を強く否定していることを示し、それによって彼らの熱進化の理解を精緻化した。