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)$ 决定。
   • 非相对论极限： 作者为核心和壳层推导了广义朗道–埃姆登（Lane–Emden）方程。在界面附近，他们定义了一个有效耦合函数 $A(\eta)$（充当密度依赖的电导率），将过渡层描述为一个扩散边界而非不连续面。
   • 相对论处理： 在各种由临界密度分析定义的参数区间内，对托尔曼–奥本海默–沃尔科夫（TOV）方程进行数值积分。

3. 诊断工具：

   • 研究分析了质量–半径（$M-R$）序列。
   • 计算了重子化学势 $\mu(\varepsilon)$，以识别相变的特征（例如，平台期或斜率变化）。

主要贡献与结果

• 参数空间分类： 作者绘制了 $(\tau_1, \sigma_1)$ 参数空间图，以识别基于正临界密度（$\varepsilon_c$）存在及其数量的不同机制。

  • 无正根的机制对应于单调、硬化的 EoS。
  • 具有一个根的机制表现出轻微的曲率变化（软化）。
  • 具有两个正根的机制表现出“类旋节线（spinodal-like）”行为，具有两个拐点，这在不需要真实微观物理相变的情况下，模拟了热力学特征（亚稳态和混合相）。

• 低质量中子星构型： TOV 方程的数值积分得到了引力质量范围在 $0.99\text{--}2.05 M_\odot$ 之间的稳定恒星构型。具体而言，该模型成功重现了质量约为 $1.0\text{--}1.1 M_\odot$ 且半径 $\approx 11$ km 的构型，这与近期的低质量中子星候选天体相一致。

• “相变”的本质：

  • 对化学势 $\mu(\varepsilon)$ 的分析表明，尽管模型表现出强烈的非线性和曲率变化，但它并未产生教科书式一级相变所特有的平坦平台。
  • 相反，“类相变”行为表现为 $\mu(\varepsilon)$ 和 $dp/d\varepsilon$ 在显著高于名义核心-壳层匹配点 ($\varepsilon_0$) 的密度处发生的显著斜率变化。
  • “物理”转换区是由 vdW 型核心项之间的竞争驱动的，且位于恒星核心深处，而非人工设置的分段连接处。

• 界面动力学： 非相对论分析表明，核心-壳层界面表现为一个受有效系数 $A(\eta)$ 控制的薄扩散层。在具有两个临界密度的机制中，$A(\eta)$ 在界面附近出现极小值，表明在压力项竞争最剧烈的地方存在局部的有效电导率降低。

意义与主张
本文声称其主要意义在于提供了一个受控且完全解析的测试平台，用于探索相变特征和低质量构型的涌现。作者明确指出，该模型并非旨在重现详细的微观物理或特定的核相互作用。相反，它旨在证明：

1. 在简化的分段式 EoS 中，少量的参数可以控制复杂热力学行为（如有效相变）和低质量中子星构型的出现。
2. 通过化学势诊断的相变特征位置并不一定与分段式 EoS 的结构连接点重合；相反，它反映了高密度下的内部非线性。
3. 低质量中子星（$M \lesssim 1.1 M_\odot$）与简化的现象学模型是兼容的，这表明仅凭粗略的结构论据无法排除其存在的可能性，尽管它们的实际形成取决于超新星动力学和微观物理。

研究结论认为，尽管该模型是现象学的，但它有效地追踪了压力-能量密度关系中的非线性如何从 EoS 影响到全局恒星结构，为研究在解析模型中何时会出现类相变特征提供了一个透明的框架。