An Analytical Two Incompressible Fluid Star with a Mixed Ordinary Dark Matter Core and an Ordinary Matter Envelope

技術要約：混合型通常物質・ダークマターコアと通常物質エンベロープを持つ、解析的な二成分非圧縮流体星

問題提起
本論文は、コンパクト星における（ダークマターと解釈される）独立して保存される第二の物質成分の重力効果を調査するための、解析的に扱いやすいモデルの必要性に取り組んでいる。定密度の一成分流体星であるシュヴァルツシルト内部解（SIS）は、一般相対性理論における基本的なベンチマークであるが、ダークマターが混入した現実的なシナリオでは、複雑な数値計算が必要となるか、あるいは解析的な制御が困難になることが多い。具体的には、著者らは、コア内に共存しつつも星の表面には到達しない、中心に集中したダークマター成分の影響を分離することを目的としている。この構成は、異なる流体が個別のコア領域とシェル領域を占める従来の二層モデルとは異なり、ここでは流体がコア内で共存しており、一方の流体のみが表面まで延びている。

手法
著者らは、以下の二つの領域からなる静的かつ球対称なモデルを構築している。

1. 領域 I (コア、$0 \le r \le R_i$): 二つの独立して保存される非圧縮完全流体（通常物質 $\rho_o$ およびダークマター $\rho_D$）を含む。これらは共通の時空幾何学を介してのみ相互作用し、物質やエネルギーの直接的な交換はない。全エネルギー密度は $\rho_T = \rho_o + \rho_D$ である。
2. 領域 II (エンベロープ、$R_i \le r \le R_e$): 定密度 $\rho_o$ の通常物質のみを含む。
3. 領域 III (外部、$r > R_e$): シュヴァルツシルト計量によって記述される真空。

系はアインシュタイン場の方程式および各流体の保存則 ($\nabla_\nu T^{\mu\nu}_o = 0$ および $\nabla_\nu T^{\mu\nu}_D = 0$) に従う。著者らは、結合されたトルマン・オッペンハイマー・ヴォルコフ（TOV）方程式を解析的に解いている。主な手順は以下の通りである：

• コア内における全圧力 $p_T(r)$ および個別の圧力 $p_o(r), p_D(r)$ の明示的な表現の導出。
• 両領域における計量関数 $\Phi(r)$ および $\Psi(r)$ の解出。
• 境界条件の適用：界面 $R_i$ における計量成分 ($g_{tt}, g_{rr}$) および圧力の連続性、ならびに表面 $R_e$ における圧力の消失 ($p_o(R_e)=0$) およびダークマターのコア境界における圧力の消失 ($p_D(R_i)=0$)。
• 無次元パラメータの定義：密度比 $f = \rho_D/\rho_o$、半径比 $\alpha = R_e/R_i$、およびコンパクトネスパラメータ $C_i = 2M_i/R_i$ と $C_e = 2M_e/R_e$。

主な貢献と結果

1. 解析解: 主要な貢献は、この混合コア・単一エンベロープ構成に対する圧力プロファイルと計量関数の明示的な解析表現の導出である。中心圧力 $p_{TC}$ は自由パラメータではなく、$f, \alpha,$ およびコンパクトネスに依存する境界条件および整合条件によって決定される。
2. ブッダハル的な臨界コンパクトネス: 著者らは、中心圧力の発散 ($p_{TC} \to \infty$) に関連する臨界条件を導出した。これは、二成分流体系におけるブッダハル的な限界を定義するものである。
   • 臨界コンパクトネスは、ダークマターの割合 $f$ と相対的なコアサイズ $\alpha$ に明示的に依存する。
   • ダークマターが消失する極限 ($f \to 0$) かつコアが星を満たす極限 ($\alpha \to 1$) では、標準的なシュヴァルツシルトの結果 $2M/R = 8/9$ が回収される。
   • 数値解析により、固定された $\alpha$ に対して、コアのコンパクトネス $C_i$ の臨界値は、ダークマターの割合 $f$ が増加するにつれて減少することが明らかになった。
3. パラメータ空間と許容性: 計量関数が実数であり、圧力が有限かつ正である物理的に許容されるパラメータ空間を特定している。図1は、様々な $\alpha$ 値に対する $(C_i, f)$ 平面における臨界境界を示しており、小さな $f$ においては約リニアな関係があることを示している。
4. 質量-半径関係と縮退: 質量-半径系列の分析により、ダークマターの存在は与えられた半径に対して全質量を増加させる（$1 + f\alpha^{-3}$ に比例してスケールする）一方で、異なる内部構成が同一のグローバルなコンパクトネス $C = 2M/R$ を持ち得ることが示された。具体的には、$f$ と $\alpha$ の異なる組み合わせを持つ構成が、同一のグローバルなコンパクトネスを持ちながら、異なる内部物質分布を持つことが可能である。
5. 計量の振る舞い: 時間成分の計量 $g_{tt}$ は、コア・エンベロープ界面において連続であることが示されているが、その径方向のプロファイルは、閉じ込められたダークマター成分によって修正されている。

意義と主張
本論文は、この研究を、現実的なダークマター混入中性子星の微視的なモデルとしてではなく、「制御された解析的実験室」として位置づけている。著者らは、この構成が以下のための有用なベンチマークを提供すると主張している：

• 中心に集中した第二成分の重力的影響を分離すること。
• 詳細な状態方程式を含むより現実的なモデルでテスト可能な、定性的な傾向（$f$ や $\alpha$ に対する臨界コンパクトネスの依存性など）を特定すること。
• 二成分混合コアから単一成分エンベロープへの遷移を研究するための、扱いやすい枠組みを提供すること。

著者らは、導出された臨界コンパクトネスは中心圧力の発散（静的な許容条件）に関連するものであり、動的な安定性の境界として解釈すべきではないことを明記している。動的な安定性を確立するには、結合された径方向の摂動に関する別途の解析が必要であると述べている。本モデルは、複雑な微視物理学から重力的効果を切り離すためのベースラインを提供するためのものである。このモデルは、重力的効果を複雑な微視物理学から分離し、将来的なコンパクト天体へのダークマター混入に関する研究のための基礎を提供するものである。