Low-Mass Neutron Stars and Effective Phase Transitions from a Hybrid Van der Waals-Polytropic Equation of State

技術要約：低質量中性子星とハイブリッド・ファン・デル・ワールス–ポリトロープ状態方程式による実効的な相転移

問題提起
中性子星（NS）の内部構造と状態方程式（EoS）は、高エネルギーおよび核天体物理学における中心的な課題である。観測データは、大質量パルサー（$\sim 2 M_\odot$）の存在を裏付けている一方で、中性子星の低質量限界については依然として未解決の問題であり、近年では$\sim 0.77\text{--}0.98 M_\odot$程度の質量を持つ候補も示唆されている。これらの天体をモデル化する際の主な困難は、現実的な微視的EoSモデルが複雑な多体計算や膨大なパラメータセットを伴うことが多く、解析的に扱うことが困難である点にある。この複雑さは、軟化、硬化、あるいは相転移に類する挙動を駆動している特定の成分を系統的に特定することを妨げている。さらに、高密度物質における相転移は理論的に予想されているものの、その兆候は完全な数値モデルの中ではマクロな観測量の中に隠されてしまうことが多い。著者らは、詳細な微視的物理学に依存することなく、簡略化された解析的に扱いやすいEoSが、相転移に関連する主要な現象論を再現し、かつ低質量中性子星の構成を支持できるかどうかという問題に取り組んでいる。

手法
著者らは、修正ファン・デル・ワールス（vdW）型のコアとポリトロープ型の地殻を結合させた、区分的なハイブリッド状態方程式を構築している。このモデルは、非線形な密度依存性を組み込みつつ、解析的に扱いやすく設計されている。

1. 状態方程式の構築:

   • コア ($\varepsilon > \varepsilon_0$): 修正vdW式によってモデル化される：
     $$p(\varepsilon) = \frac{\alpha_1}{(\varepsilon_0/\varepsilon) - 1} \left(\frac{\varepsilon}{\varepsilon_0}\right)^{\tau_1} + \beta_1 \left(\frac{\varepsilon}{\varepsilon_0}\right)^{\sigma_1}$$
     ここで、$\alpha_1$は反発的な熱的寄与を制御し、$\beta_1$は相互作用に起因する高密度項を表し、指数$\tau_1$および$\sigma_1$はこれらの項の密度依存性を支配する。温度は緩やかな密度依存性 $T(\varepsilon) \propto \varepsilon^{\tau_1}$ を持つと仮定される。
   • 地殻 ($\varepsilon \le \varepsilon_0$): 二項ポリトロープ形式によってモデル化される：
     $$p(\varepsilon) = \alpha_2 \left(\frac{\varepsilon}{\varepsilon_0}\right)^{\tau_2} + \beta_2 \left(\frac{\varepsilon}{\varepsilon_0}\right)^{\sigma_2}$$
   • 界面: コアと地殻の間の遷移は、半径 $r_1$ の周囲の「ぼやけた」領域として扱われ、圧力とエネルギー密度は急激な不連続性を強制するのではなく、近似的に一致させられる。

2. 理論的枠組み:

   • 熱力学的制約: モデルは熱的安定性（$dp/d\varepsilon \ge 0$）および因果律（$v_s^2 \le c^2$）を課す。
   • 臨界密度: 著者らは、圧力－エネルギー密度関係 $p(\varepsilon)$ の曲率が変化する臨界エネルギー密度（$\varepsilon_c$）に関する条件を導出する。これらはパラメータ $(\tau_1, \sigma_1)$ によって決定される。
   • 非相対論的極限: 著者らは、コアおよび地殻に対する一般化されたレーン・エムデン方程式を導出する。界面付近では、遷移層を不連続な境界ではなく拡散境界として記述するために、有効結合関数 $A(\eta)$（密度依存の導電率として機能する）を定義する。
   • 相対論的取り扱い: 様々なパラメータ領域（臨界密度解析によって定義される）にわたって、トルマン・オッペンハイマー・ヴォルコフ（TOV）方程式を数値積分する。

3. 診断ツール:

   • 研究では質量－半径（$M-R$）系列を分析する。
   • 相転移の兆候（例：プラトーや傾きの変化）を特定するために、バリオン化学ポテンシャル $\mu(\varepsilon)$ が算出される。

主要な貢献と結果

• パラメータ空間の分類: 著者らは、正の臨界密度（$\varepsilon_c$）の存在およびその数に基づいて、$(\tau_1, \sigma_1)$ パラメータ空間をマッピングした。

  • 正の根を持たない領域は、単調で硬いEoSに対応する。
  • 一つの根を持つ領域は、緩やかな曲率の変化（軟化）を示す。
  • 二つの正の根を持つ領域は、「スピノーダル的」な挙動を示し、二つの変曲点を持ち、真の微視的な相転移を必要とせずに、第一種相転移の熱力学的兆候（メタ安定性と混合相）を模倣する。

• 低質量中性子星の構成: TOV方程式の数値積分により、重力質量が $0.99\text{--}2.05 M_\odot$ の範囲にある安定した恒星構成が得られた。具体的には、本モデルは質量が $1.0\text{--}1.1 M_\odot$、半径が $\approx 11$ km の構成を正常に再現しており、これは近年の低質量中性子星の候補と一致している。

• 「相転移」の性質:

  • 化学ポテンシャル $\mu(\varepsilon)$ の分析により、本モデルは強い非線形性と曲率の変化を示すものの、教科書的な第一種相転移に特徴的な平坦なプラトーは生成しないことが明らかになった。
  • 代わりに、「相転移的な」挙候は、名目上のコア－地殻のマッチング点 $\varepsilon_0$ よりも著しく高い密度において、$\mu(\varepsilon)$ および $dp/d\varepsilon$ の顕著な傾きの変化として現れる。
  • 「物理的な」遷移領域は、vdW型のコア項の競合によって駆動されており、人工的な区分的な接合部ではなく、恒星コアの深部に位置している。

• 界面のダイナミクス: 非相対論的解析は、コア－地殻の界面が有効係数 $A(\eta)$ によって制御される薄い拡散層として振る舞うことを示している。二つの臨界密度が存在する領域では、$A(\eta)$ は界面付近で極小値をとり、これは圧力項の競合が最も激しい場所で局所的な有効導電率の低下が生じることを示唆している。

意義と主張
本論文の主要な意義は、相転移の兆候および低質量構成の出現を探求するための、制御された完全な解析的テストベッドを提供することにあると著者らは主張している。著者らは、このモデルが詳細な微視的物理学や特定の核相互作用を再現することを意図していないことを明示している。むしろ、以下のことを示すためのものである。

1. 簡略化された区分的なEoSにおける少数のパラメータによって、複雑な熱力学的挙動（実効的な相転移など）や低質量中性子星の構成を制御できること。
2. （化学ポテンシャルによって診断される）相転移の兆候の所在は、必ずしも区分的なEoSの構造的接合点と一致するわけではなく、むしろ高密度における内部の非線形性を反映していること。
3. 低質量中性子星（$M \lesssim 1.1 M_\odot$）は簡略化された現象論的モデルと両立可能であり、その存在は粗い構造的議論のみによって否定されることはないが、その実際の形成は超新星爆発のダイナミクスや微視的物理学に依存すること。

本研究は、モデルは現象論的ではあるものの、圧力－エネルギー密度関係の非線形性が恒星の全体的な構造に与える影響を追跡する上で効果的であり、解析的モデルにおいてどのような条件下で相転移的な特徴が生じるかを調査するための透明性の高い枠組みを提供すると結論付けている。