A self-consistent single-fluid framework for neutron stars admixed with mirror dark matter

技術要約：ミラーダークマターが混入した中性子星のための自己整合的な単一流体フレームワーク

問題提起
ダークマター（DM）の性質は、その質量、結合強度、および相互作用モードが直接検出実験によって未決定であるため、依然として開かれた問いである。中性子星（NS）のようなコンパクト天体は、間接的なDM検出のための潜在的な天体物理学的実験室として機能するが、既存のモデリング手法には重大な限界がある。

1. 二流体形式論： これらはDMとバリオン物質（BM）を、重力を介してのみ相互作用する別々の流体として扱う。これにより、DMコアやハローの形成は可能だが、状態方程式（EOS）を変化させ得る非重力的な相互作用を捉えることができない。
2. 既存の単一流体アプローチ： これらは通常、DMフェルミ運動量（$k_F^\chi$）を外部パラメータとして固定する。その結果、星全体でDM数密度（$n_\chi$）がほぼ一様になり、高密度コアに向かって相対的なDM割合（$n_\chi/n_B$）が減少してしまう。このプロファイルは、散乱と熱化がDM含有量を局所的なバリオン環境に結合させるはずであるという、現実的な捕獲シナリオと矛盾する。
3. 中性子崩壊平衡モデル： 一部の手法は、DM集団を決定するために $\mu_\chi = \mu_n$ を課す。これは極限的なケースとしては有用であるが、有限の捕獲率や散乱断面積を無視しているため、DM存在量を過大評価する可能性があり、強い反発的な自己相互作用を人工的に導入しない限り、非物理的なDM集団を招く恐れがある。

したがって、現実的なレベルでDM含有量を制御し、局所的なDM密度がバリオン密度を追跡することを可能にし、かつ相互作用の強さを粒子物理学的な制約に直接結びつける、自己整合的なフレームワークが必要とされている。

手法
著者らは、ダークセクターが可視バリオンセクターの正確なコピー（同一の粒子内容、質量、および相互作用を持つ）であるミラーダークマターのシナリオに基づいた、自己整合的な単一流体フレームワークを開発した。主要な手法構成要素は以下の通りである：

• 相互作用メカニズム： 両セクターは、接触ベクトル電流–電流相互作用（$L_{DB} = \alpha (\bar{\psi}_D \gamma^\mu \psi_D)(\bar{\psi}_B \gamma^\mu \psi_B)$）を介して結合されている。これは、両セクターの有効化学ポテンシャルに対して、それぞれ $-\alpha n_B$ および $-\alpha n_D$ という相互平均場シフトを生成する。
• グローバルな制約： フェルミ運動量を固定する代わりに、本フレームワークではグローバルなDM比率 $F_D = N_D/N_B$（全ダーク粒子数と全バリオン粒数の比）を固定する。両セクターは同じ恒星体積を占有するため、これは局所的な関係 $n_D = F_D n_B$ を課し、DM密度が星全体を通じてバリオン密度と自己整合的にスケールすることを保証する。
• 状態方程式（EOS）： 本研究では、対称エネルギーの勾配が異なる4つの相対論的平均場（RMF）モデル、NL3$\omega\rho$、FSU2R、NL3、およびDDME2（熱力学的整合性のために再配置自己エネルギー項を必要とする密度依存モデル）を用いている。
• パラメータ制約： 結合定数 $\alpha$ は、地上での直接検出実験（具体的には $\sim 1$ GeV のDMに対する CRESST-III の制限）によって制約されたスピン非依存のDM–核子散乱断面積（$\sigma_{DN}$）に直接関連付けられている。
• 数値解法： 結合された平均場方程式は、バリオン数保存および $\beta$ 平衡（および電荷中性）が各セクターに対して個別に課された状態で同時に解かれる。得られたEOSは、恒星構造を決定するためのトルマン・オッペンハイマー・ヴォルコフ（TOV）方程式を解くために使用される。

主な貢献

1. 熱力学的整合性： 本フレームワークは、固定密度処置に代わる物理的に動機付けられた選択肢を提供し、局所的なDM密度が自然にバリオン物質密度に従うことを可能にする。これは、固定された $k_F$ モデルによる人工的な密度プロファイルや、平衡のみを考慮したモデルにおけるDMの過大評価を回避する。
2. ミクロ・マクロの架け橋： 結合強度 $\alpha$ を $\sigma_{DN}$ に関連付けることで、本研究は、コンパクト天体の観測量と地上実験からの粒子物理学的制約との間に直接的な定量的リンクを確立している。
3. 統一的扱い： DMを単なる重力成分としてではなく、非重力的相互作用を通じてEOS、粒子組成、およびグローバルな恒星特性を変化させる、動的に結合した流体として扱っている。

結果
このフレームワークをミラーDMが混入した中性子物質に適用したところ、テストされたすべてのEOSモデルにおいて以下の系統的な知見が得られた：

• EOSの軟化： DM–BM相互作用は、高密度物質の結合を弱め、その非圧縮性（$K_0$）を減少させる。$F_D$ が 0% から 10% へ増加するにつれ、EOSは一定のエネルギー密度においてより低い圧力へとシフトする。
• マクロな変化：
  • 最大質量： 最大支持質量（$M_{max}$）は約 7–9% 減少する（例：NL3$\omega\rho$ の場合、2.76 から 2.56 $M_\odot$ へ）。
  • コンパクト性： 中性子星はよりコンパクトになり、中心密度（$n_c$）が増加し、1.4 $M_\odot$ の星の半径（$R$）は約 0.8–1.0 km 減少する。
  • 結合エネルギー： 単位核子あたりの結合エネルギーは負の値が小さくなり、システム全体の安定性の低下を示す。
• 微視的な組成： 相互作用は物質をより中性子に富んだ組成へと駆動する。DM誘起の化学ポテンシャルシフトの存在下で $\beta$ 平衡と電荷中性を満たすために、陽子分率（$Y_p$）は特定の密度において系統的に抑制される。
• ダイレクト・ウリカ（DU）過程： 陽子分率の抑制により、ダイレクト・ウリカ過程の開始密度（$n_{DU}$）はより高い恒星密度へとシフトする。
  • 硬い対称エネルギーを持つモデル（例：NL3）では、このシフトにより急速冷却の開始がより質量の大きい星へと移動する。
  • 軟い対称エネルギーを持つモデル（例：NL3$\omega\rho$, FSU2R）では、コンパクト性の増加により、開始がより質量の小さい星へと移動するか、あるいは DDME2 の場合、支持される質量範囲内で DU 過程が全く発生しない可能性がある。
• 観測との適合性： DM混入を持つ DDME2 モデルは、マルチメッセンジャー制約（GW170817、HESS J1731-347、および大質量パルサー）と最も良好な一致を示している。これは $\sim 2 M_\odot$ の質量限界を維持しつつ、半径の制約にも近づいているが、HESS J1731-347 の領域を満たすには比較的大きな $F_D$ を必要とする。

意義と主張
本論文は、ミラーDMの混入が中性子星の微視的な組成（中性子への富化の増大）とマクロな構造（コンパクト性の増大、最大質量の減少）の両方を根本的に変容させることを主張している。

本研究の意義は、熱力学的に整合し、観測によって制約されたプラットフォームを提供することにある。著者らは、本フレームワークを用いることで、中性子星の測定値が直接検出実験によって探索されているのと同じパラメータ空間内で解釈可能になると主張している。予測されるシグネチャー、具体的には最大質量の減少、コンパクト性の増大、およびダイレクト・ウリカ冷却開始のシフトは、マルチメッセンジャー・シグネチャーとなり得るものである。これは、コンパクト星の極限環境におけるダークマター物理学をテストするための有望な手段を提供し、この単一流体アプローチを従来の二流体モデルや固定密度モデルから区別するものである。