Kick bimodality of neutron stars and mode dependence of their parameters

问题陈述
观测表明，中子星（NSs）拥有的空间速度显著高于其前身星，这一现象归因于超新星爆发过程中获得的初生踢（natal kicks）。虽然理论建模和群体合成（例如 Igoshev 2020）已支持存在双峰踢速度分布（包含低速和高速两种模式），但观测确认仍具挑战性。由于径向速度分量未知、出生地不确定以及在银河系势场中的运动演化，从观测到的二维横向速度推导初生踢分布并非易事。此外，现有研究尚未深入调查特定无线电脉冲星物理参数是否与这些不同的踢速度模式相关，特别是要区分是由选择偏差还是内在物理差异所致。

方法论
作者分析了选自 ATNF 脉冲星目录（v2.7.0）的 202 颗正常孤立无线电脉冲星样本。筛选标准要求具有可靠的视差测量值或色散量距离 < 10 kpc，且具有已知的自行、自转参数（$P, \dot{P}$），并排除了回收脉冲星、联星系统及奇异源。

为了将这些脉冲星分类为低速或高速模式，作者采用了轨迹积分模型：

1. 轨迹积分： 对于每颗脉冲星，基于视差或色散量采样 100 个距离值。利用银河坐标和自行，生成 100 个现今位置。
2. 速度重建： 由于径向速度（$v_r$）无法观测，作者假设相对于局部标准静止系（LSR）的速度方向是各向同性的。他们从截断均匀分布（Uniform(0.05$\pi$, 0.95$\pi$)）中对视线与三维速度向量之间的夹角 $\theta$ 进行采样，以避免极端异常值，从而为每个距离生成 $10^3$ 个速度向量。这导致每颗脉冲星产生 $10^5$ 个位置-速度组合。
3. 反向积分： 使用银河系势场（McMillan 2017），将轨迹进行反向积分，时间跨度达 $1.5 \times \tau_c$（特征年龄）。出生地通过与银盘的交点或在 $t = -\tau_c$ 时（若未发生交点）来确定。
4. 踢速度计算： 在每个潜在出生地，计算踢速度（$v_{kick} = |\vec{v} - \vec{v}_{LSR}|$）。
5. 分类： 利用 Igoshev (2020) 的双峰踢分布参数（两个麦克斯韦分布的加权和：$\sigma_1 = 56$ km s$^{-1}$, $\sigma_2 = 336$ km s$^{-1}$, 权重 $w_1=0.2, w_2=0.8$），作者计算了每颗脉冲星属于低速或高速模式的概率。如果脉冲星在所有可能的出生地处概率均超过 0.5，则将其归类为“明确”属于某一模式；否则，根据多数出生地的结果进行归类。

关键结果

• 模式分布： 样本中约 30%（202 颗中的 60 颗）属于低速模式，这与此前对年轻脉冲星的估计一致。其余约 70% 属于高速模式。
• 选择偏差效应： 正如预期，高速脉冲星平均而言比低速脉冲星更年轻（自转减慢年龄较短）且距离更远。这些差异通过 Kolmogorov-Smirnov (K-S) 测试得到证实（$p < 0.01$），并被归因于观测选择效应（在较大距离处难以测量微小速度，以及高速天体离开巡天体积）。
• 脉冲宽度： 在两种模式之间未发现脉冲宽度（$W_{10}$ 和 $W_{50}$）存在显著差异（$p \approx 0.5$），这表明该参数不受强选择偏差或内在模式依赖性的影响。
• 磁场相关性： 在表面磁场分布中检测到了统计学上的显著差异（$p \approx 0.02$）。低速模式的脉冲星在低磁场（$B \lesssim 10^{12}$ G）方面表现出过剩。值得注意的是，在高速模式中未发现磁场 $B \lesssim 10^{11}$ G 的脉冲星。即使在考虑了距离和年龄后，这种差异依然存在，尽管作者指出这可能部分受到磁场随时间衰减的影响。

意义与主张
本文声称提供了一种将正常无线电脉冲星详细分类为双峰踢速度组的方法，并分析了由此产生的参数分布。其主要意义在于发现了踢速度模式与磁场强度之间的相关性。作者指出，虽然选择偏差和磁场衰减可能导致观察到的差异，但这种差异的起源尚不完全明确。他们认为，如果双峰踢结构源于超新星爆发物理学的差异，那么磁场差异可能具有类似的起源，即可能将中子星踢机制与外部偶极磁场联系起来。然而，作者保持了谦逊的态度，承认需要详细的群体合成建模来确认这些起源，且不能完全排除选择效应的影响。本研究并非提出新的实验应用，而是通过对现有数据的统计分析来完善对中子星运动学和磁性质的理解。