Radiative decays of dynamically generated pentaquarks in the chiral unitary approach: the $P_c(4457)\to P_c(4312)\,\gamma$ transition

技术摘要：手征幺正性方法中动态生成五夸克态的辐射衰变

问题陈述
LHCb 实验观测到了三个窄共振态五夸克态（$P_c(4312)$、$P_c(4440)$ 和 $P_c(4457)$），它们出现在 $J/\psi p$ 能谱中。虽然分子模型表明这些是由于耦合道相互作用产生的 $S$ 波强子分子，但其内部结构仍有待最终确认。检验分子假设的一个关键测试是辐射跃迁 $P_c(4457) \to P_c(4312)\gamma$。在分子图像中，光子耦合到组成部分的介子-重子组分；而在紧凑夸克模型（双夸克-双夸克-反夸克或强子粲偶素）中，光子则耦合到夸克核。以往的研究在处理此类辐射跃迁时存在局限性，通常仅考虑单一主导的介子辐射圈，或者忽略了完整的耦合道内容及复杂的残数相位。本文旨在解决系统性、规范不变地处理动态生成五夸克态辐射衰变的需求，并特别将其应用于 $P_c(4457) \to P_c(4312)\gamma$ 跃迁。

方法论
作者采用了结合重夸克自旋对称性和局部隐函数相互作用的手征幺正性方法。$P_c(4312)$ 和 $P_c(4457)$ 被视为由粲介子与粲重子的散射而动态生成的极点。

• 形式体系： 计算利用耦合道振幅 $T = [1 - V G]^{-1} V$，其中 $V$ 是相互作用核，$G$ 是介子-重子圈函数。这些态由复能量平面中的极点及其残数（$g_i$）来定义。
• 辐射振幅： 光子耦合到两个极点的完整介子-重子组分。该跃迁通过 19 个不同的 $M1$ 三角圈进行建模。这些圈被分为两类：Class M（光子从介子线上发射）和 Class B（光子从重子线上发射）。
• 规范不变性： 通过包含所有必要的图来确保计算的规范不变性。由于 $S$ 波顶点中接触项（Kroll-Ruderman 项）消失，因此 19 个三角圈之和构成了完整的规范不变振幅。
• 顶点与对称性：
  • 强相互作用顶点由耦合道解的残数确定。
  • 重介子的电磁顶点（$D^* \to D\gamma$）归一化至实验数据（$\bar{D}^{*0} \to \bar{D}^0\gamma$）。
  • 轻夸克磁矩和重子跃迁矩使用朴素夸克模型和重夸克自旋对称性进行估算。
  • 计算明确保留了残数的复相位，这对于干涉效应至关重要。
• 计算： 19 个圈积分被简化为具有闭合解析形式的单数值求积。作者还测试了结果对圈动量上软高斯形式因子的敏感性，以便与之前的唯象研究进行比较。

主要贡献

1. 完整的耦合道处理： 不同于以往侧重于单一主导通道（如 $\bar{D}^*\Sigma_c$）的工作，本研究包含了通过电磁顶点连接两个极点的所有开放和封闭通道，从而形成了 19 个三角圈之和。
2. 相位敏感性： 本文证明了耦合道解中复数残数的相对相位显著影响干涉模式和最终的衰变宽度。
3. 系统性框架： 它提供了一个用于计算动态生成态辐射衰变的规范不变模板，可应用于其他跃迁（例如 $P_c(4440) \to P_c(4312)\gamma$）以及奇数扇区（$P_{cs}$ 态），且无需引入新的参数。

结果

• 衰变宽度： 辐射宽度 $\Gamma(P_c(4457) \to P_c(4312)\gamma)$ 的中心值为 6.7 keV，保守范围约为 2 至 9 keV。
• 主导贡献：
  • 最主要的贡献来自 $\bar{D}^{*0} \to \bar{D}^0\gamma$ 介子辐射圈（圈 a）。
  • 第二大贡献来自于近阈值的 $\bar{D}^*\Lambda_c$ 圈（圈 e），涉及 $\Sigma_c \to \Lambda_c \gamma$ 跃迁。
  • 所有圈的相干和大于非相干和（3.82 keV），表明存在相长干涉。
• 相位敏感性： 宽度对 $\bar{D}^*\Lambda_c$ 残数与 $\Sigma_c \to \Lambda_c$ 磁矩之间的相对符号高度敏感。破坏性干涉情景使宽度降至 ~1.9 keV，而相长干涉情景则产生 ~6.7 keV。
• 形式因子效应： 应用与手征幺正性方法重整化标度一致的软高斯形式因子，会将宽度降低至约 4.2 keV。由于主导圈具有长程性质，其受形式因子的影响较小，但总宽度下降了约三分之一。
• 级联率： 假设中心宽度值，级联衰变 $\Lambda_b^0 \to J/\psi p K^- \gamma$ 的分支比约为 $1.3 \times 10^{-9}$。
• 文献对比： 该结果大于早期的分子计算（参考文献 [34, 36]），后者给出的结果为 1.1–2.2 keV。这种差异归因于包含了完整的通道内容、开放的 $\bar{D}^*\Lambda_c$ 通道以及复数残数相位的相长干涉。当对主导图应用相同的软形式因子时，结果与早期的发现一致，这证实了差异在于额外的圈和相位。

意义与主张
论文声称提供了第一个关于动态生成态五夸克辐射衰变的完整、规范不变的计算。作者提出了三个具体的观测物理量来测试这些态的分子性质：

1. 多极矩成分： 预测该跃迁是纯 $M1$ 模式。显著的 $E2$ 成分将与最小的 $S$ 波分子图像相矛盾。
2. 辐射线比例： 提出的比率 $R_\gamma = \Gamma(P_c(4457) \to P_c(4312)\gamma) / \Gamma(P_c(4440) \to P_c(4312)\gamma)$ 是一个干净的观测量，对自旋分配和残数结构敏感，且能很大程度上抵消模型依赖的不确定性。
3. 结合能依赖性： 宽度对结合能（特别是 $P_c(4457)$ 极点靠近 $\bar{D}^*\Sigma_c$ 阈值的程度）的敏感性，为分子波函数提供了一个不同于线形分析的约束。

作者得出结论，测量该辐射宽度不仅可以测试分子假设，还可以确定无法仅通过线形数据获得的残数相对相位。该方法被呈现为一个稳健的框架，可扩展到其他奇异强子。