Algebraic coefficients of fractional parentage for three particles in the oscillator Jacobi basis

技术摘要：三粒子振荡子雅可比基底下的代数系数（Fractional Parentage）

问题陈述
在平移不变的少体核物理学中，计算通常使用雅可比坐标（Jacobi coordinates）进行，以消除虚假的质心激发。然而，在这一基底中强制执行泡利原理（反对称性）具有挑战性。与单粒子基底不同（在单粒子基底中，粒子交换仅重新标记状态），在雅可比基底中交换核子会混合坐标。因此，寻找轨道分数亲缘系数（CFP）——即在置换群 $S(3)$ 下正确变换的基态线性组合——需要对每一个能量（$E$）和角动量（$L$）块内的置换算符（具体为类和算符 $\hat{\Lambda} = \hat{P}_{13} + \hat{P}_{23}$）进行对角化。

标准方法涉及利用谐振子（HO）括号构建该算符的矩阵，并应用数值特征值求解器。这一过程计算成本高昂（对于块维度 $n$，复杂度为 $O(n^3)$），且必须在每个块内重复进行，并会在本应是精确整数的量中引入数值舍入误差。此外，生成的特征向量通常需要事后分析以识别其对称性类别和配对结构。

方法论
作者提出了一种闭式构造方法，消除了对数值对角化的需求。该方法的核心依赖于一个特定的酉变换 $\hat{W} = \hat{H}_1 \hat{\Phi}$，它对类算符进行了“修饰”（dressing）：

• $\hat{H}_1$：质量参数为 $d=1$ 的谐振子括号算符。
• $\hat{\Phi}$：一个固定的复数化算符，由提升至振荡子空间的对角相位矩阵 $F = \text{diag}(1, i)$ 定义。

推导过程是将一个几何恒等式从雅可比坐标平面提升到振荡子 Fock 空间。在平面上，通过 $F$ 和 Hadamard 镜像 $T_1$（对应于 $\hat{H}_1$）对两个转置反射之和（$P_{13} + P_{23}$）进行共轭运算，可以将该和转化为一个带有实系数的简单标签交换。

导出的核心恒等式为：
$$\hat{W} \hat{\Lambda} \hat{W}^\dagger = \sigma \hat{J} \hat{Z}$$
其中：

• $\sigma = (-1)^{E-L}$ 是块内的常数相位。
• $\hat{J}$ 是标签交换算符（交换两个雅可比量子数对）。
• $\hat{Z}$ 是一个对角算符，其整数特征值为 $z_i = 2\cos[\pi(e_1 - e_2)/3] \in \{\pm 1, \pm 2\}$。

这一变换表明，当被 $\hat{W}$ 修饰后，类算符不再是一个需要对角化的稠密矩阵，而是一个简单的标签置换乘以一个整数。

主要贡献与结果

1. 闭式 CFP：CFP 被显式构造为：
   $$|K_{i,s}\rangle = \sqrt{2} \hat{Q}_s \hat{W}^\dagger |i\rangle$$
   其中 $\hat{Q}_s$ 是选择特定 $e_2$-宇称的投影算符。作者证明，每一个 CFP 仅仅是 $d=1$ HO 括号表的 $\sqrt{2}$ 倍缩放、带符号装饰的列，并受限于一个宇称扇区。
2. 精确的对称性标记：类算符的特征值是确定的整数（$\pm 2, \pm 1$），直接由状态的量子数决定，从而无需进行数值特征值提取。对称性类别（$S(3)$ 不可约表示）可立即分配：
   • $z_i = \pm 2$ 对应于全对称 $[3]$ 或全反对称 $[111]$ 态。
   • $|z_i| = 1$ 对应于混合对称 $[21]$ 对。
3. 算法效率：生成 CFP 的计算成本从 $O(n^3)$（特征值求解器）降低到 $O(n^2)$（矩阵组装和列选择）。特征值求解器实际上从 CFP 步骤中被移除，并完全转移到了为 $\hat{H}_1$ 提供支持的括号引擎中。
4. 数值验证：
   • 该方法已用 Fortran 实现，并在高达 $E=50$ 的 1,301 个块（总计 912,366 个 CFP）上进行了测试。
   • $S(3)$ 类分配在每个块中都是精确的。
   • 残余数值误差被发现仅源于用于计算 $\hat{H}_1$ 的基于 SU(3) 的括号引擎。在双精度下，对于 $E \le 35$，误差保持在 $10^{-8}$ 以下。对于更高能量，使用四倍精度构建使误差恢复到 $10^{-22}$ 或更好。
   • 该构造通过了一个独立的生产级三核子代码（Trys_kunai）进行了验证。在完整的流水线中用代数生成的 CFP 替换数值生成的 CFP，得到了相同的 $^3\text{He}$ 基态能量（差异 $< 10^{-9}$ MeV），证实了该方法的物理正确性。

意义
本文声称，在振荡子雅可比基底中确定轨道 CFP 时，特征值求解器是不必要的。系数已经“表格化”在 $d=1$ 谐振子括号本身之中。通过利用特定的酉算符 $\hat{W}$，宇称、相位和对称性分类的复杂簿记工作坍缩为了一个单一的算符表达式。

其意义在于：

• 统一性：一个公式适用于所有 $E$、$L$ 和对称性类别，包括此前需要单独处理的情况（例如 $e_1=e_2$ 的中心区域）。
• 精确性：该方法产生精确的整数特征值，并避免了迭代对角化固有的舍入误差累积。
• 性能：虽然由于括号引擎本身的成本，整体加速效果有所缓和，但在 $E=50$ 时，CFP 特定的步骤比最快的可用特征值求解器快 15 倍。

作者指出，剩余的数值限制并非源于 CFP 的构造，而是源于括号引擎（SU(3) 塔构造）的具体实现。他们建议，一个闭式括号引擎将消除流水线中的所有数值误差。这项工作也指出了向超球谐波基底和四粒子情形扩展的潜力。