Modeling quasielastic lepton-nucleus interactions with ab initio spectral functions from infinite nuclear matter

技术摘要：利用从头算谱函数对准弹性轻子-原子核相互作用进行建模

问题陈述
下一代中微子振荡实验（如 DUNE、T2HK）需要精确的中微子-原子核相互作用理论模型，以控制系统误差。一个核心挑战是在单一框架内一致地处理不同的相互作用通道和末态相互作用（FSI）。虽然谱函数（SF）形式化方法已被广泛使用，但 FSI 往往被忽略或仅进行近似处理（例如通过光学势）。此外，现有的微观电弱截面计算多局限于使用格林函数蒙特卡洛或耦合簇理论的包容性可观测物理量，往往缺乏对粒子部分 FSI 的一致处理。因此，需要一种能够整合用于初始（空穴）态和末态（粒子）态的“从头算”（ab initio）谱函数的框架，以改进准弹性（QE）散射的描述，特别是在中微子实验相关的低到中动量转移能区。

方法论
作者提出了一个基于局部密度近似（LDA）和脉冲近似（IA）的框架，用于模拟包容性轻子-原子核散射。

1. 理论基础：

   • 谱函数 (SFs)： 本研究利用了基于自洽格林函数（SCGF）理论（采用 Gorkov 框架）推导出的无限对称核物质（NM）的“从头算”谱函数。这些计算采用了来自手征有效场论的现实核哈密顿量（具体为 $\Delta$N2LOGO(394)、$\Delta$N2LOGO(450) 和 N2LOsat(450)），并使用了“ADC(3)-D”截断层级的代数图解构造（ADC）方法。
   • LDA 实现： 有限原子核（$^{12}$C、$^{16}$O、$^{40}$Ca）通过将无限核物质谱函数对核密度分布 $\rho(r)$ 进行积分来建模。该框架同时包含了描述初始束缚核子的空穴谱函数（$S_h$）和考虑 FSI 的粒子谱函数（$S_p$）。
   • 两种模型： 考虑了两种运动学机制：
     • IA+FSI： 对空穴部分和粒子部分均使用“从头算”核物质谱函数。
     • PWIA（平面波脉冲近似）： 使用“从头算”空穴谱函数，但用相对论费米气体模型替换粒子谱函数（仅考虑泡利阻挡），这适用于非相对论“从头算”谱函数不再可靠的高动量转移区域。

2. 神经网络插值：

   • 由于 SCGF 计算仅在能量、动量和密度处提供离散的谱函数点，作者开发了人工神经网络（ANN）来连续插值这些数据。
   • 共训练了六个独立的网络（对应三种哈密顿量与两种同位旋（质子/中子）的每种组合）。
   • 输入特征为能量 ($E$)、动量 ($|p|$)、密度 ($\rho$) 和平滑因子 ($\Gamma$)。输出为谱函数值。
   • 采用了物理启发式的损失函数，通过引入惩罚项来强制执行归一化求和规则（$\int dE S(E) = 1$）。

3. 截面计算：

   • 通过对核体积内空穴和粒子谱函数的乘积进行积分来构建强子张量。
   • 该框架计算了电磁（EM）响应和带电流（CC）响应，包括用于中微子散射的通量平均总截面。

主要贡献

• 一致的 FSI 处理： 引入了一种一致的处理 FSI 的方法，即使用从“头算”得到的粒子谱函数，而非依赖于唯象的光学势。
• 神经网络插值： 开发了用于插值高维“从头算”谱函数数据的 ANN，从而使其能够用于任意运动学变量的截面计算。
• 基准测试： 将模型与针对 $^{12}$C、$^{16}$O 和 $^{40}$Ca 的电子散射实验数据，以及针对 $^{16}$O 的“从头算”LIT-CC（洛伦兹积分变换结合耦合簇）计算进行了基准测试。
• 中微子预测： 该框架被应用于预测 $^{16}$O 的带电流响应以及 $\nu_\mu$-$^{12}$C 的通量平均总截面，并与 LSND 实验数据进行了对比。

结果

• 电子散射：
  • 对于 $^{12}$C，IA+FSI 模型在不同动量转移下，在 $\omega \lesssim \omega_{QE}$ 时与 Rosenbluth 分离的实验纵向响应（$R_L$）数据表现出极佳的一致性。
  • FSI 的引入修正了 QE 峰的位置和强度，而 PWIA 模型往往会将 QE 峰向高能方向移动并高估强度。
  • 横向响应（$R_T$）被系统性地低估了，这与计算中排除了二体电流的情况一致。
  • 对于 $^{16}$O，IA+FSI 结果在 QE 峰附近与 LIT-CC 计算吻合良好，尽管 LDA 方法缺乏 LIT-CC 中观察到的详细壳层结构阈值效应。
  • 对于 $^{40}$Ca，该模型很好地重现了实验数据，在纵向和横向通道中均表现出良好的符合度。
• 中微子散射：
  • 使用 IA+FSI 模型计算的 $\nu_\mu$-$^{12}$C 通量平均总截面与 LSND 实验数据一致。
  • PWIA 模型给出的截面比 IA+FSI 低约 20%，这是由于在较低动量下，QE 峰的偏移导致其超出了积分窗口。
  • 理论不确定性（源于哈密顿量选择和平滑因子 $\Gamma$）在动量转移 $|q| \approx 0.2\text{--}0.6$ GeV 的 QE 峰处评估约为 7–15%。

意义与主张
作者声称，他们的框架提供了“在我们的近似有效范围内与实验数据的一致良好结果”（特别是低到中动量转移区域）。其意义在于：

1. 提高一致性： 在 LDA 内为初始态和末态提供了一致的微观“从头算”描述，这对于中微子振荡实验尤为重要，因为宽能谱通量会抹除详细的壳层结构。
2. 可扩展性： 该框架旨在易于扩展，以包含额外的动力学机制（如介质中 $\pi$ 产生）和其他核哈密顿量。
3. 验证： 通过电子散射数据和 LIT-CC 计算的成功基准测试，验证了该方法在同位旋对称原子核中处理包容性可观物理量的有效性。

论文也谦虚地指出了局限性：LDA 牺牲了详细的壳层结构（使其不太适用于排他性 $(e, e'p)$ 数据），且该模型在极高动量转移下会失效，届于此时相对论修正和超越 QE 的机制（如 $\pi$ 产生）变得占据主导地位。未来的工作计划包括将 SCGF 计算扩展到同位旋非对称物质（与 DUNE 中的 $^{40}$Ar 相关），并将该框架与内核级联模型相结合。