AstroLineage

在线性修正引力 $f(R,T) = R + 2\chi T$ 框架下，流体声速限制 $0 < c_s^2 < 1$（其中 $c_s^2 = \frac{1}{3} + \frac{8\pi}{3\chi}$）约束耦合参数 $\chi \in [-8\pi, 4\pi]$。推导出修正的 Buchdahl 紧致度极限 $\beta_{\rm B} = \frac{128\pi^2 + 8\pi\chi - \chi^2}{288\pi^2}$ 以及主能量条件极限 $\beta_{\rm DEC} = \frac{48\pi^2 + 8\pi\chi - \chi^2}{128\pi^2}$。其轴向微扰有效势 $V(x) = \frac{e^{\nu}}{x^3} \left\{ l(l+1)x+ 4\pi x^3 \left[\rho\left(1+\frac{\chi}{4\pi}\right)-p\left(1+\frac{3\chi}{4\pi}\right)\right] -6 m(x)\right\}$ 具有更深的势阱；采用 NICER 观测半径 $R = 11.4\,\mathrm{km}$ 时，在 $\chi=0.96\pi$、紧致度 $\beta=0.4662702$ 下计算得到引力波回声时延 $\tau_{\rm echo} = 2.39\times 10^2\,\mathrm{ms}$，对应特征回声频率 $f_{\rm echo} = 72.11\,\mathrm{Hz}$，与 GW170817 报道的 $\sim 72\,\mathrm{Hz}$ 信号一致（广义相对论预言为 $5.19\times 10^2\,\mathrm{Hz}$）。

研究发现平滑跨越导致物态方程相较纯强子相变软，降低了中子星最大质量；但当有效袋常数满足 $B_{\rm eff} \ge 180\;{\rm MeV}$ 时，最大质量仍可支撑超过 $2\,M_\odot$ 的天文观测下限。在过渡区，声速平方 $c_s^2$ 产生峰值并随 $B_{\rm eff}$ 增大而更尖锐，高密度下渐近趋于共形极限 $1/3$。对于给定的夸克模型，相变交叉化学势 $\mu_{\rm cros}$ 随着强子物态刚度参数 $\eta$ 的增加而指数减小，且得到的无量纲潮汐形变度 $\Lambda$ 与 GW170817 观测限制高度吻合。

通过放宽最大质量星中夸克核大小先验（仅要求 $P_{\rm max} - P_{\rm trans} > 5\,\mathrm{MeV\,fm^{-3}}$），推断得到的相变密度移向更高密度区间 $\rho_{\rm trans} \in [0.484, 0.741]\,\mathrm{fm^{-3}}$，相变临界质量阈值提升至 $M_{\rm trans} \ge 2.012\,M_\odot$。混合星强子相 EoS 整体偏软，得到更小的半径 $R_{1.4} = 12.964^{+0.334}_{-0.239}\,\mathrm{km}$ 与最大质量 $M_{\rm max} = 2.168^{+0.114}_{-0.121}\,M_\odot$，能更好地吻合 NICER 与引力波潮汐形变观测约束。

在单流体近似下，Higgs 门户模型约束暗物质粒子质量为 $\log_{10}(m_\chi/\mathrm{MeV}) = 3.38_{-0.87}^{+1.05}$（中位值约为 $2.4\,\mathrm{GeV}$），其有效相互作用耦合 $y = 0.104_{-0.063}^{+0.062}$；ALP 模型中核物质偏度系数被显著限制在 $Q_0 = -927.57_{-55.54}^{+69.23}\,\mathrm{MeV}$。各模型预言的最大质量相近（$M_{\rm max} \approx 2.09 - 2.13\,M_\odot$），但在 $1.4\,M_\odot$ 处 Higgs 门户模型呈现更紧致结构（$R_{1.4} = 12.24\,\mathrm{km}$，$\Lambda_{1.4} \approx 361$）及更高的 $f$ 模振荡频率（$\nu_{f,1.4} = 1.83\,\mathrm{kHz}$）。贝叶斯证据量排序显示中子衰变与 ALP 方案优于 Higgs 门户模型。

1. 在零温非旋转（$T=0$）条件下，纯核子（N）和常数耦合超子（Y）模型最大引力质量分别为 $M_G = 2.484\ M_\odot$（重子质量 $M_B = 2.996\ M_\odot, R = 11.997\ \mathrm{km}$）和 $M_G = 2.086\ M_\odot$（$M_B = 2.417\ M_\odot, R = 12.421\ \mathrm{km}$），而密度依赖超子耦合（Y DD）最大质量仅 $1.958\ M_\odot$；$1.4\ M_\odot$ 中子星半径分别为 $13.13\ \mathrm{km}$ 和 $13.12\ \mathrm{km}$，潮汐形变 $\Lambda_{1.4}$ 为 $653$ 和 $651$。2. 中微子囚禁效应会降低给定重子数密度下的奇异数（超子）丰度。3. 应用于引力波事件表明，GW190425（$M_G \sim 3.50\ M_\odot$）在所有模型下均发生直接坍缩（prompt collapse），而 GW170817（$M_G \sim 2.74\ M_\odot$）可形成长寿命残骸。

联合 GW170817、NICER 及大质量脉冲星观测后验，中子星典型参量收敛至 $R(1.4) = 11.7 \pm 0.4\,\mathrm{km}$、$\Lambda(1.4) = 309 \pm 90$、$M_{\rm max} = 2.4 \pm 0.2\,M_{\odot}$。在包含 1000 次双中子星并合的大星表模拟中，经 3 轮自适应提纯后 $\Lambda(1.4)$ 不确定度压窄至 $\pm 16$（相对精度 $\sim 4\%$），成功避免了固定先验下有效样本数从 $N_{\rm eff} \approx 470$ 崩溃至 $22$ 的问题；但向微观核对称能映射受本征退化限制，$\sigma_{L_{\rm sym}}$ 仅从先验的 $29\,\mathrm{MeV}$ 缩小至 $21\,\mathrm{MeV}$，$\sigma_{K_{\rm sym}}$ 几乎保持在 $114\,\mathrm{MeV}$ 未变。

正文总结称，NICER 观测后，约两倍核饱和密度处的压强不确定性以及典型 1.4 个太阳质量中子星的半径不确定性，均约为观测前的一半。作者认为已研究的若干系统误差通常不会造成大的半径偏差，但个别脉冲星的模型、背景和采样问题仍须谨慎处理。

研究发现对于中等相对论性膨胀（$\beta \lesssim 0.3$），光行时延迟效应（使边缘温度升高）与多普勒增亮效应（使视线方向 $\mu \approx 1$ 增亮）作用方向相反并发生近似抵消，从而使观测温度呈现近各向同性，合理解释了 AT2017gfo 的准黑体特征；若忽略该相对论几何动力学效应，推断的动力学温度将被高估最高 $7.6\%$、速度 $\beta$ 将被高估约 $11.9\%$。对 AT2017gfo 拟合给出有效膨胀速度随时间自 $\beta \simeq 0.199 \pm 0.009$ 降至 $0.122 \pm 0.006$，有效观测温度从 $5428\ \mathrm{K}$ 降至 $3061\ \mathrm{K}$，简并组合量 $\kappa M_{\mathrm{ej}}$ 自 $(3.8^{+4.0}_{-1.5})\times 10^{-2}$ 微增至 $(6.5^{+3.9}_{-2.4})\times 10^{-2}\ \mathrm{cm^2\,g^{-1}\,M_\odot}$。

在 Rubin/LSST 巡天节奏下，仅靠千新星光变反演并合时刻 $t_0$ 的 $90\%$ 置信区间宽度中位数为 $5.49\,\mathrm{days}$；若在并合前 $3\,\mathrm{h}$ 引入两个相邻滤光片的测光上限，后验宽度可收缩约 $60\%$；红移先验（固定 $D_L$）带来的改善 $<5\%$。在 AT2017gfo 测光结合单点合成上升期观测的最优场景下，$t_0$ 估计不确定度可压缩至 $0.03\,\mathrm{days}$（$0.72\,\mathrm{hours}$）。对于余辉扣除后的 GRB-KN 样本（如 GRB 211211A），$1\sigma$ 精度可达 $8.7\,\mathrm{hours}$；若首个探测晚于并合 $\gtrsim 4\,\mathrm{days}$（如 GRB 060614），后验弥散则宽达 $>120\,\mathrm{hours}$。

正文报告：偶偶 α 发射核的最优指数为 a=0.66584±0.00014；固定 a=2/3 后，log10(T1/2)=6.653 X−48.975，调整后 R²=0.99061，残差标准差为 0.50 dex。六种理论描述的最优指数平均为 0.660，模型间标准差为 0.028。对36个已测核采用等权加热求和，在0.03–3天内得到局部指数−1.24±0.01。替换222Rn和224Ra寿命后，作者报告加热及光变与实测寿命输入大体一致；但正文关于光变一致性的量化表述存在不一致，不能据此确认其已准确再现 AT2017gfo。

在 Rubin/LSST 巡天节奏下，仅靠千新星光变反演并合时刻 $t_0$ 的 $90\%$ 置信区间宽度中位数为 $5.49\,\mathrm{days}$；若在并合前 $3\,\mathrm{h}$ 引入两个相邻滤光片的测光上限，后验宽度可收缩约 $60\%$；红移先验（固定 $D_L$）带来的改善 $<5\%$。在 AT2017gfo 测光结合单点合成上升期观测的最优场景下，$t_0$ 估计不确定度可压缩至 $0.03\,\mathrm{days}$（$0.72\,\mathrm{hours}$）。对于余辉扣除后的 GRB-KN 样本（如 GRB 211211A），$1\sigma$ 精度可达 $8.7\,\mathrm{hours}$；若首个探测晚于并合 $\gtrsim 4\,\mathrm{days}$（如 GRB 060614），后验弥散则宽达 $>120\,\mathrm{hours}$。

在 Rubin/LSST 巡天节奏下，仅靠千新星光变反演并合时刻 $t_0$ 的 $90\%$ 置信区间宽度中位数为 $5.49\,\mathrm{days}$；若在并合前 $3\,\mathrm{h}$ 引入两个相邻滤光片的测光上限，后验宽度可收缩约 $60\%$；红移先验（固定 $D_L$）带来的改善 $<5\%$。在 AT2017gfo 测光结合单点合成上升期观测的最优场景下，$t_0$ 估计不确定度可压缩至 $0.03\,\mathrm{days}$（$0.72\,\mathrm{hours}$）。对于余辉扣除后的 GRB-KN 样本（如 GRB 211211A），$1\sigma$ 精度可达 $8.7\,\mathrm{hours}$；若首个探测晚于并合 $\gtrsim 4\,\mathrm{days}$（如 GRB 060614），后验弥散则宽达 $>120\,\mathrm{hours}$。