Automated event generation for S-wave quarkonium and leptonium production in NRQCD and NRQED

基本信息

• 论文ID: 2510.26773
• 标题: Automated event generation for S-wave quarkonium and leptonium production in NRQCD and NRQED
• 作者: Alice Colpani Serri, Chris A. Flett, Jean-Philippe Lansberg, Olivier Mattelaer, Hua-Sheng Shao, Lukas Simon
• 分类: hep-ph (高能物理-唯象学), hep-ex (高能物理-实验), nucl-ex (核物理-实验), nucl-th (核物理-理论)
• 提交时间: 2025年10月30日
• 论文链接: https://arxiv.org/abs/2510.26773

摘要

本文提出了MadGraph5_aMC@NLO框架的扩展，能够在非相对论QCD (NRQCD)和非相对论QED (NRQED)因子化形式下自动计算S波夸克偶素(quarkonium)和轻子偶素(leptonium)产生的领头阶截面。该框架已通过多种基准过程验证，展示了其稳健性和现象学研究的灵活性。该实现的关键优势在于与现有MadGraph5_aMC@NLO功能的无缝集成，不仅支持标准模型计算，还可通过修改的通用费曼输出(UFO)接口处理各种超出标准模型或有效场论场景。此外，该框架与标准蒙特卡洛事件生成器的部分子簇射和强子化功能保持兼容。通过大量实例，作者强调夸克偶素过程的理论研究需要仔细考虑：次领头贡献的影响往往难以仅通过基于耦合层级和速度标度规则的简单计数论证来预测。

研究背景与动机

研究问题

该研究旨在解决夸克偶素和轻子偶素产生过程的自动化事件生成问题，特别是在NRQCD和NRQED理论框架下的实现。

问题重要性

1. 物理意义：夸克偶素(如J/ψ、Υ)是研究QCD微扰和非微扰方面的重要探针。自J/ψ在50年前发现引发"十一月革命"以来，夸克偶素产生机制仍未完全理解。
2. 实验需求：LHC、EIC等对撞机产生大量夸克偶素数据，需要可靠的理论预测工具进行比对。
3. 新物理探索：轻子偶素(如正电子偶素、真μ子偶素、双τ偶素)可用于精确检验QED、寻找新物理和测量基本常数。

现有方法局限性

1. MadOnia：基于MadGraph/MadEvent v4，仅支持单个S波或P波夸克偶素，未移植到现代MG5_aMC框架。
2. HELAC-Onia：支持多个S波和P波夸克偶素，但在处理物理束缚态时需要手动组合短程系数。
3. 通用MC生成器：Pythia、Herwig、Sherpa中的夸克偶素矩阵元和部分子簇射功能尚未完全兼容。
4. 轻子偶素工具缺失：公开可用的轻子偶素产生理论工具极为有限。

研究动机

开发一个全面的、自动化的、用户友好的工具，能够：

• 在NRQCD/NRQED框架下处理任意数量的S波夸克偶素和轻子偶素
• 无缝集成到广泛使用的MG5_aMC框架
• 支持标准模型和超标准模型场景
• 为未来NLO精度计算奠定基础

核心贡献

1. 首次实现：在MG5_aMC框架中首次实现S波夸克偶素和轻子偶素的自动化事件生成
2. 统一框架：提供了处理夸克偶素和轻子偶素的统一NRQCD/NRQED实现
3. 新UFO模型：开发了sm_onia模型，包含专门的boundstates.py文件定义Fock态内容
4. 广泛验证：通过与HELAC-Onia和解析表达式的详细对比验证了实现的正确性
5. 丰富应用：展示了在pp、e⁺e⁻、e⁻p对撞机上的各种产生过程，包括：
   • 单夸克偶素产生及其与电弱玻色子、喷注的关联产生
   • 双/三夸克偶素产生
   • 轻子偶素(正电子偶素、真μ子偶素、双τ偶素)产生
6. 理论洞察：通过大量实例说明次领头贡献的重要性不能简单通过速度标度规则估计

方法详解

理论框架

NRQCD/NRQED因子化

束缚态B的微分截面可表示为： $d\sigma(N_1 + N_2 \to B + X) = \sum_{I_1,I_2,n} \int dx_1 dx_2 f_{I_1/N_1}(x_1) f_{I_2/N_2}(x_2) d\hat{\sigma}(I_1I_2 \to (C_1C_2)[n] + X_p) \langle \mathcal{O}^B_n \rangle$

其中：

• $(C_1C_2)$：重夸克对$(Q\bar{Q}')$或轻子对$(\ell^-\ell'^+)$
• $n = {}^{2S+1}L_J^{[C]}$：量子态(自旋S、轨道角动量L、总角动量J、色态C)
• $d\hat{\sigma}$：短程截面(微扰QCD/QED可计算)
• $\langle \mathcal{O}^B_n \rangle$：长程矩阵元(LDME)

投影方法

通过协变投影方法将振幅投影到特定Fock态：

1. 色投影 (仅对夸克偶素)：
   • 色单态：$P_{C=1} = \delta_{c_3c_4}/\sqrt{N_c}$
   • 色八重态：$P_{C=8} = \sqrt{2} t^{c_{34}}_{c_4c_3}$
2. 自旋投影： $P_S = \frac{\bar{v}_{\lambda_{C_2}}(k_4) \Gamma_S u_{\lambda_{C_1}}(k_3)}{2\sqrt{2m_{C_1}m_{C_2}}}$ 其中$\Gamma_{S=0} = \gamma_5$(单态)，$\Gamma_{S=1} = \not{\varepsilon}^*_{\lambda_s}(K)$(三重态)
3. 轨道角动量投影 (S波：L=0)： $\mathcal{A}^{(n,0)}_{\{[C],S,L=0\}}(r) = \left[\mathcal{A}^{(n,0)}_{\{[C],S\}}(r)\right]_{q=0}$

实现架构

UFO模型扩展

1. 新类Boundstate：类似于Particle类，定义束缚态的Fock态
2. 命名规范：Jpsi(1|3S11)表示J/ψ的主量子数N=1，自旋三重态³S₁，色单态
3. PDG编码：
   • 色单态：使用标准PDG编码(如J/ψ: 443)
   • 色八重态：非标准编码99nqnsnrnLnJ
   • 轻子偶素：新编码±990nℓ₁nℓ₂nrnLnJ

代码生成流程

1. 色投影：在Python层面生成费曼图时应用
2. 自旋投影：通过ALOHA模块生成专用子程序，依赖HELAS库
3. LDME处理：直接整合到振幅平方计算中，通过onia_card.dat输入
4. 相空间积分：采用伪多通道策略，使用VEGAS积分器并行化

运行模式

1. Leading-order模式：完整的LO事件模拟，生成LHEF格式事件
2. Standalone模式：仅生成矩阵元计算代码(Fortran)，用于外部框架集成

技术创新点

1. 自动Fock态组合：用户指定物理束缚态(如Jpsi)，代码自动包含所有相关Fock态贡献
2. 灵活LDME设置：运行时可调整LDME值，无需重新编译
3. BSM兼容性：通过UFO接口支持任意新物理模型(如HEFT示例)
4. 部分子簇射接口：与Pythia等MC生成器无缝对接
5. 轻子偶素支持：首次提供轻子偶素的自动化计算能力

实验设置

计算设置

基本参数

• 夸克质量：$m_c = 1.55$ GeV (3FS), $m_b = 4.7$ GeV, $m_t = 173$ GeV
• 轻子质量：$m_e = 0.511$ MeV, $m_\mu = 105.66$ MeV, $m_\tau = 1.777$ GeV
• 电弱参数：$m_Z = 91.188$ GeV, $m_H = 125$ GeV, $\alpha^{-1}_{G_\mu} = 132.507$
• 重整化/因子化标度：$\mu_R = \mu_F = H_T/2 = \frac{1}{2}\sum_i \sqrt{k^2_{i,T} + m^2_i}$

对撞机设置

1. 质子-质子对撞 (LHC)：$\sqrt{s} = 13$ TeV, PDF4LHC21_40 NNLO PDF
2. 电子-正电子对撞 (Belle/BaBar)：$\sqrt{s} = 10.58$ GeV (Υ(4S)共振)
3. 电子-质子对撞 (EIC)：$\sqrt{s} = 140.7$ GeV ($E_{e^-} = 18$ GeV, $E_p = 275$ GeV)

LDME值

采用表2中的默认值：

• 色单态：通过波函数$R_Q(0)$计算，如$\langle \mathcal{O}^{J/\psi}_{{}^3S_1^{[1]}} \rangle = 1.16$ GeV³
• 色八重态：从实验数据拟合，如$\langle \mathcal{O}^{J/\psi}_{{}^1S_0^{[8]}} \rangle = 0.0146$ GeV³

验证策略

矩阵元层面验证

• 与HELAC-Onia对比42个不同过程的振幅平方
• 相对偏差$\Delta_{\text{rel}} < 10^{-6}$（浮点精度水平）

截面层面验证

• 总截面和微分截面与HELAC-Onia对比
• 拉力(pull)值$< 3$，与统计期望一致

解析表达式验证

• 轻子偶素过程与文献中的解析公式对比
• 在$Z$和$H$玻色子退耦极限下，千分之一精度一致

实验结果

主要结果展示

1. 单夸克偶素产生 (pp对撞，√s = 13 TeV)

2→1过程总截面：

• $\sigma(pp \to \eta_c + X) = 2.9366(5)$ μb
• $\sigma(pp \to J/\psi + X) = 536.14(6)$ nb
• $\sigma(pp \to \eta_b + X) = 5.4935(7)$ μb
• $\sigma(pp \to \Upsilon + X) = 6.0655(4)$ nb

关键发现：

• 赝标量态截面远大于矢量态，因矢量态的2→1 QCD过程因C宇称守恒被禁止
• $\eta_c$截面对低$x$区域PDF极为敏感，不同PDF集可导致截面相对$\eta_b$的大小反转

2. 关联产生 (pp对撞)

与电弱玻色子关联：

• $\sigma(pp \to J/\psi + Z) = 1.3425(4)$ pb (色八重态主导)
• $\sigma(pp \to \Upsilon + Z) = 340.7(1)$ fb (色单态主导)

与喷注/光子关联 (应用运动学切割)：

• $\sigma(pp \to J/\psi + j) = 329.8(2)$ nb ($p_{T,j} > 10$ GeV)
• $\sigma(pp \to J/\psi + \gamma) = 19.13(1)$ nb ($p_{T,\gamma} > 2$ GeV)

3. 双夸克偶素产生

同种夸克偶素对：

• $\sigma(pp \to J/\psi + J/\psi) = 10.756(3)$ nb (色单态主导)
• $\sigma(pp \to \Upsilon + \Upsilon) = 44.63(1)$ pb

混合产生：

• $\sigma(pp \to J/\psi + \Upsilon) = 10.503(3)$ pb (色单态QCD过程在$\mathcal{O}(\alpha^4_s)$消失)

$B_c$对产生：

• $\sigma(pp \to B^+_c + B^-_c) = 472.5(1)$ pb
• $\sigma(pp \to B^{\ast+}_c + B^{\ast-}_c) = 940.9(2)$ pb

4. 三J/ψ产生

$\sigma(pp \to J/\psi + J/\psi + J/\psi) = 1.038(3)$ pb

三色单态贡献仅54.57 fb，因$gg \to 3J/\psi[^3S_1^{[1]}]$因C宇称守恒被禁止，主要来自$q\bar{q}$初态。

轻子偶素产生结果

正电子偶素 (Ps)

pp对撞 ($\sqrt{s} = 13$ TeV)：

• $\sigma(pp \to \text{Ps}_1 + j) = 9.483(2)$ fb

e⁺e⁻对撞 ($\sqrt{s} = 10.58$ GeV)：

• $\sigma(e^+e^- \to \text{Ps}_0 + \gamma) = 3.300(1)$ ab
• $\sigma(e^+e^- \to \text{Ps}_1 + \gamma) = 8.420(2)$ ab

真μ子偶素 (TM)

pp对撞：

• $\sigma(pp \to \text{TM}_1 + j) = 9.460(2)$ fb (与正电子偶素几乎相同)

e⁺e⁻对撞：

• $\sigma(e^+e^- \to \text{TM}_0 + \gamma) = 438.16(1)$ yb (极小，因$4m^2_\mu/s$抑制)
• $\sigma(e^+e^- \to \text{TM}_1 + \gamma) = 8.422(2)$ ab

双τ偶素 (T)

pp对撞：

• $\sigma(pp \to \text{T}_1) = 9.1906(4)$ fb (直接产生，首次可靠计算)
• $\sigma(pp \to \text{T}_1 + j) = 6.055(1)$ fb
• $\sigma(pp \xrightarrow{\gamma\gamma} p\text{T}_0p) = 108.055(7)$ ab (UPC过程)

e⁺e⁻对撞：

• $\sigma(e^+e^- \to \text{T}_0 + \gamma) = 110.001(3)$ zb
• $\sigma(e^+e^- \to \text{T}_1 + \gamma) = 9.576(2)$ ab

特殊案例：J/ψ + H产生 (HEFT框架)

在Higgs有效场论中，通过有效$Hgg$顶点：

• $\sigma(pp \to J/\psi + H) = 1.53^{+0.40}_{-0.29}$ fb (仅色八重态贡献)
• $\sigma(pp \to \Upsilon + H) = 91^{+23}_{-17}$ ab

快度分布特征：

• ${}^1S_0^{[8]}$态在中心快度区出现凹陷($|y| \approx 4$处有峰)
• ${}^3S_1^{[8]}$态在中心快度区有峰
• 凹陷源于有效$Hgg$顶点图，被Yukawa耦合图部分抵消

部分子簇射效应

对$pp \to J/\psi + c\bar{c}$过程，比较固定阶(fLO)和匹配部分子簇射(LO+PS)：

• ${}^3S_1^{[1]}$和${}^1S_0^{[8]}$：PS展宽$p_T$分布
• ${}^3S_1^{[8]}$：PS收窄$p_T$分布
• 总效果：部分抵消，导致适度畸变

消融分析：Fock态贡献

以$pp \to \eta_c$为例，分解为6个部分子通道：

• $gg \to \eta_c[^1S_0^{[1]}]$：2.2995 μb (主导)
• $gg \to \eta_c[^1S_0^{[8]}]$：33.43 nb
• $q\bar{q} \to \eta_c[^3S_1^{[8]}]$：603.7 nb
• 其他通道：$\sim 0$ (量子数守恒禁止)

重要洞察：简单的速度标度估计常失效，因：

1. 量子数守恒导致某些通道消失
2. 传播子增强效应(如$g^* \to (Q\bar{Q})[^3S_1^{[8]}]$)
3. 运动学/动力学抑制或增强

相关工作

通用MC生成器

• Pythia v8.3/Herwig v7.4：引入夸克偶素特定部分子簇射，但与矩阵元功能尚未完全兼容
• 局限性：夸克偶素产生主要通过强子化或衰变，矩阵元实现有限

过程无关的部分子级事件生成器

• MadOnia：基于MadGraph v4，支持单个S/P波夸克偶素，未移植到MG5_aMC
• HELAC-Onia：支持多个S/P波夸克偶素，但需手动组合Fock态，LDME不可运行时调整
• 本工作优势：
  1. 现代MG5_aMC框架，维护活跃
  2. 自动Fock态组合
  3. 运行时可调LDME
  4. 首次支持轻子偶素

过程特定生成器

• SuperChic/STARlight/eSTARlight：专注于排他过程
• EPOS4：重离子碰撞中的夸克偶素，采用Wigner密度矩阵形式
• BCVEGPY：专门用于$B_c$介子强子产生
• 本工作优势：通用性，支持任意S波束缚态组合

截面计算器

• FDCHQHP：S/P波夸克偶素强子产生的NLO计算
• FONLL：非瞬发J/ψ的NLO+NLL计算
• 本工作定位：LO自动化，为NLO扩展奠定基础

符号计算工具包

• FeynOnium/AmpRed：半自动符号计算
• 本工作区别：完全自动化数值计算和事件生成

结论与讨论

主要结论

1. 技术成就：成功在MG5_aMC框架中实现S波夸克偶素和轻子偶素的自动化事件生成
2. 验证充分：通过与HELAC-Onia和解析表达式的详细对比，确认实现的正确性
3. 应用广泛：涵盖pp、e⁺e⁻、e⁻p对撞机的各种产生过程
4. 理论洞察：次领头贡献的重要性不能简单通过速度标度规则估计，需具体过程具体分析
5. 轻子偶素：首次提供系统的轻子偶素产生预测，为未来实验观测提供指导

局限性

1. LO精度：当前仅实现领头阶，NLO修正可能显著（如三J/ψ产生）
2. S波限制：仅支持S波态，P波态实现留待未来工作
3. DPS缺失：不包含双/三部分子散射贡献，这在某些过程中可能主导
4. 前馈效应：未考虑激发态衰变的前馈贡献
5. 相空间积分：采用简化的伪多通道策略，对极复杂过程效率可能不足
6. 低能区域：如$pp \to \text{Ps}_1$等极低能过程，因子化假设可能失效

未来方向

1. NLO实现：基于扩展的FKS减除方案实现NLO精度（已有理论基础67）
2. P波支持：实现轨道和总角动量投影算符
3. 多通道优化：开发针对束缚态的高效相空间参数化
4. 全局拟合：结合实验数据进行LDME全局拟合
5. DPS整合：与DPS框架结合，提供完整预测
6. 重离子碰撞：扩展到核PDF和冷核物质效应

物理影响

1. 实验指导：为LHC、EIC、Belle II等实验的夸克偶素/轻子偶素研究提供可靠预测
2. 机制理解：帮助区分色单态/色八重态机制、SPS/DPS贡献
3. 新物理探索：支持BSM场景下的束缚态产生研究
4. 精密测试：为QCD/QED精密检验提供理论工具

深度评价

优点

技术创新性

1. 架构设计优秀：通过引入Boundstate类和修改UFO格式，优雅地将束缚态整合到MG5_aMC框架
2. 用户友好：命名规范清晰（如Jpsi(1|3S11)），LDME运行时可调，学习曲线平缓
3. 模块化实现：色/自旋/轨道投影分离，便于扩展和维护
4. 多层验证：矩阵元、截面、解析公式三层验证，确保可靠性

科学贡献

1. 首创性：首次在主流MC框架中实现轻子偶素自动化计算
2. 系统性：涵盖3种对撞机类型、40+种过程，展示广泛适用性
3. 洞察深刻：通过大量实例说明速度标度规则的局限性，纠正可能的误解
4. 实用价值：为实验物理学家提供易用工具，降低理论计算门槛

论文质量

1. 结构清晰：理论-实现-验证-应用的逻辑链完整
2. 细节充分：提供完整的参数设置和可复现的命令示例
3. 图表丰富：68页包含5个图和10个表，信息密度高
4. 文献全面：297篇参考文献，充分梳理相关工作

不足与改进建议

方法局限

1. 相空间效率：伪多通道策略对高多重态过程效率有限，建议：
   • 开发束缚态特定的重要性采样
   • 利用机器学习优化相空间映射
2. 低能行为：轻子偶素在极低能区的因子化假设需更仔细论证
3. LDME不确定性：采用固定LDME值，未系统评估不确定性影响

实验设置

1. 切割选择：某些运动学切割（如$p_{T,j} > 2$ GeV）较为任意，缺乏系统性研究
2. PDF依赖性：仅使用单一PDF集，未评估PDF不确定性
3. 标度依赖性：部分过程未给出标度变化不确定性

物理讨论

1. 机制区分：虽展示不同Fock态贡献，但未深入讨论如何实验区分
2. 可观测量：主要关注总截面和单粒子分布，缺少极化、关联等可观测量
3. 背景估计：未讨论实验背景和探测效率

影响力评估

学术影响

• 引用潜力：高，作为首个MG5_aMC束缚态实现，将成为该领域标准工具
• 后续研究：预期催生大量应用研究（LHC数据分析、EIC物理预测等）
• 理论发展：为NLO自动化、DPS整合等后续工作奠定基础

实用价值

• 实验应用：ATLAS、CMS、LHCb、Belle II等实验可直接使用
• 教育价值：降低束缚态物理研究门槛，有利于培养新人
• 开源贡献：代码将通过Launchpad和NLOAccess公开，促进社区发展

可复现性

• 代码可用：承诺公开发布，已提供辅助材料
• 文档详细：论文包含完整的使用示例和参数设置
• 验证充分：多重交叉验证确保结果可信

适用场景

理想应用

1. LHC物理：夸克偶素产生截面预测、机制研究、新物理搜寻
2. EIC物理：深度非弹性散射中的夸克偶素产生、胶子TMD提取
3. B工厂：精密QED检验、轻子偶素寻找
4. 理论研究：LDME拟合、色八重态机制验证、DPS研究

不适用场景

1. 重离子碰撞：缺少热/冷核物质效应，需结合EPOS4等工具
2. 超高精度预测：LO精度不足，需等待NLO实现
3. P波主导过程：如χc产生，需等待P波实现
4. 排他过程：如中心排他产生，需结合SuperChic等专用工具

建议与展望

对作者

1. 尽快实现NLO精度，这是最紧迫的需求
2. 开发P波支持，完善Fock态覆盖
3. 提供LDME全局拟合框架
4. 加强与实验组合作，推动工具应用

对用户

1. 理解LO精度的局限性，谨慎解读绝对预测
2. 进行系统的不确定性评估（标度、PDF、LDME）
3. 结合其他工具（如DPS生成器）获得完整图像
4. 反馈使用体验，促进工具改进

对领域

1. 建立标准化的LDME数据库
2. 发展束缚态物理的NLO自动化技术
3. 推动实验测量覆盖更多过程和可观测量
4. 加强理论-实验协作，实现精密QCD/QED检验

参考文献（精选）

核心理论基础

• 9 Bodwin, Braaten, Lepage (1995): NRQCD因子化形式的奠基性工作
• 64 Caswell, Lepage (1986): NRQED有效拉氏量
• 67 A A H, Shao, Simon (2024): NLO夸克偶素产生的FKS减除

相关工具

• 25 Artoisenet et al. (2008): MadOnia - 前代NRQCD自动化工具
• 27,28 Shao (2013, 2016): HELAC-Onia - 主要对比基准
• 6,7 Alwall et al. (2014, 2018): MadGraph5_aMC@NLO框架

重要应用

• 112 Lansberg, Shao (2013): J/ψ对产生的α⁵s修正重要性
• 194 Shao, Zhang (2019): 三J/ψ产生作为多部分子散射探针
• 4 d'Enterria, Shao (2023): 双τ偶素发现前景

总体评价：这是一篇高质量的方法学论文，技术实现扎实，验证充分，应用广泛。它填补了束缚态物理自动化工具的重要空白，特别是轻子偶素领域。尽管存在LO精度和S波限制等局限，但作为系列工作的第一步，已为未来发展奠定了坚实基础。预期将成为该领域的标准工具，产生广泛影响。建议实验和理论工作者密切关注其后续发展，特别是NLO实现和P波扩展。