Algebraic Structure of Three-Flavor Neutrino Oscillations in Constant-Density Matter: Cayley--Hamilton Evolution, DMP Resummation, and Closed-Form Uncertainty Propagation

技术摘要：恒定密度物质中三味中微子振荡的代数结构

问题陈述
虽然中微子振荡物理学已进入精密测量时代，但在恒定密度物质中计算振荡概率仍是下一代实验（如 DUNE、Hyper-Kamiokande (HK) 和 JUNO）的计算瓶颈。这些实验在进行拟合以解析质量排序、$\theta_{23}$ 八重态以及轻子 CP 破坏时，需要进行 $O(10^6)$–$O(10^7)$ 次评估。现有的工具（例如 Denton–Parke NuFast-LBL 算法）虽然具有高速度和高精度，但依赖于缺乏解析透明度的迭代求解器。具体而言，现有文献通常未明确阐述演化算符的代数结构、Denton–Minakata–Parke (DMP) 近似避免在 MSW 共振（$\hat{A} \to 1$）附近发生摄动发散的具体机制，以及用于 NuFIT 6.0 参数基准下严谨不确定性传播的闭式雅可比矩阵（Jacobians）。本文旨在将这些要素综合为一个统一的代数框架。

方法论
作者对恒定密度物质中的三味中微子振荡进行了全面的解析处理，其结构围绕以下支柱展开：

1. 基于 Cayley–Hamilton 的代数演化： 本文明确应用了 Cayley–Hamilton 定理处理味基哈密顿量（$H_f$）。论文证明了演化算符 $S(L) = \exp(-iH_f L)$ 严格是 $H_f$ 的二次多项式。该多项式的系数通过 $H_f$ 的三个实特征值通过一个显式的 Vandermonde 系统确定，作者指出，此前的中微子文献并未明确写出这一过程。
2. 谱解与 DMP 重求和： 特征值通过 Cardano 三角公式导出（恢复了 Zaglauer–Schwarzer 表达式）。本文分析了 Denton–Minakata–Parke (DMP) 近似，将其重新定义为一种重求和（resummation），而非仅仅是摄动展开。通过在 1–3 扇区进行关于物质混合角 $\theta^M_{13}$ 的精确旋转，DMP 方法消除了导致朴素展开在 $\hat{A} \to 1$ 时发散的近简并性。这用一个有效展开参数 $\epsilon_0 \lesssim 0.015$ 取代了无界的项 $(1-\hat{A})^{-1}$，该参数在所有能量下均保持一致有界。
3. 密度矩阵与退相干： 该形式化方法通过 Lindblad 主方程和 $su(3)$中的布洛赫矢量表示扩展到开系统动力学。这使得计算波包退相干修正成为可能，并导出了物质修饰的相干长度（$L_{coh}$）用于陆地基线。
4. 不确定性传播： 作者推导了关于六个 NuFIT 6.0 参数（$\theta_{12}, \theta_{13}, \theta_{23}, \delta_{CP}, \Delta m^2_{21}, \Delta m^2_{31}$）的闭式雅可比矩阵。这些结果通过 Feldman–Cousins 轮廓似然映射进行了验证，以确保在物理边界（例如 $\theta_{13} \to 0$）和循环变量（$\delta_{CP}$）附近的有效性。

主要结果

• 演化算符结构： 本文提供了显式的 Vandermonde 系统（式 15），将三个特征值与二次演化算符的系数联系起来，提供了一条从特征值到概率的直接代数路径，无需进行迭代矩阵指数运算。
• DMP 精度与机制： 研究表明，DMP 近似在陆地能量范围内可达到优于 $10^{-4}$ 的分数精度。本文澄清了这种精度源于对发散级数的重求和，从而产生了一个有界的展开参数 $\epsilon_0$。表 6 将此与其它方案（Madrid, AJLOS, MP）进行了对比，在这些方案中，展开参数随 $O((1-\hat{A})^{-1})$ 缩放。
• 退相干量级： 对于陆地长基线实验（T2K/HK, DUNE），计算得出基线与相干长度之比（$L/L_{coh}$）在 $10^{-3}$ 到 $10^{-2}$ 之间。虽然对于目前的精度而言退相干可以忽略不计，但本文指出，退相干并不像文献中假设的那样完全消失。
• CP 不对称性分解： 总 CP 不对称性被解析地分解为一个“真实”项（依赖于 Jarlско尔德不变量 $J_{CP}$）和一个“虚假”项（物质诱导，关于 $\hat{A}$ 为奇函数）。
• 雅可比基础设施： 提供了关于 $\partial P_{\mu e}/\partial p_k$ 的闭式表达式（式 51–56）。这些表达式允许通过单次矩阵乘法进行线性化不确定性传播，避免了在全局拟合中使用有限差分近似。

意义与主张
本文明确指出，没有任何单个结果（例如特征值解或 DMP 近似本身）是全新的。其主要贡献在于综合与代数透明度：

• 统一语言： 它将 Cayley–Hamilton 结构、DMP 的重求和解释以及密度矩阵退相干联系在一起，形成了一个单一且连贯的框架。
• 解析连续性： 不同于优化速度的迭代求解器（如 NuFast-LBL），该解析框架提供了参数连续性、透明的极限情况（真空、高能、两味情形）以及闭式雅可比矩阵。
• 分析作用： 作者将这项工作定位为不是作为生产拟合中（其中速度至关重要）的 NuFast-LBL 算法的替代品，而是作为解释、验证和不确定性传播的工具。它提供了理解极限并以闭式形式传播相关不确定性的解析基础设施，而迭代求解器由于其构造特性无法提供这些内容。
• 可扩展性： 文中展示了该代数逻辑可以自然地扩展到非标准相互作用（NSI）和惰性中微子情景（需要针对 $4\times4$ 矩阵的三次 Cayley–Hamilton 多项式）。

综上所述，本文提供了使演化算符及其导数变得显式的“簿记”和代数基础，作为精密中微子振荡物理时代的理论锚点。