A High Performance Partial Wave Analysis Framework for Hadron Spectroscopy

在亚原子世界的深处，物质不仅仅是固体粒子的集合，更是将夸克和胶子结合在一起的力的动态相互作用。这个受量子色动力学规律支配的领域，不仅是质子和中子等普通物质形成的场所，还隐藏着一个更为奇异的景观。物理学家相信，在适当的条件下，这些基本构建模块可以自我排列成奇特的、短寿命的构型，从而挑战粒子通常结合方式的标准规则。其中包括像“胶球”（glueballs）这样完全由力载体组成的“奇异”状态，或是由四个或五个夸克组成的四夸克态和五夸克态。由于这些奇异粒子在衰变为更稳定的碎片之前仅存在于转瞬之间，因此无法被直接观测到。相反，科学家必须扮演“宇宙侦探”的角色，通过分析它们留下的粒子喷射轨迹来重建它们的存在。为此，他们使用一种被称为“部分波分析”（partial wave analysis）的强大统计工具，通过筛选数百万个碰撞事件，来识别出创造这些粒子的共振态所具有的特定量子指纹。

挑战在于海量的数据规模。现代粒子探测器（例如北京的 BESIII 实验）现在记录着数百万个碰撞事件，这一规模已经使用于分析它们的传统方法不堪重负。标准方法要求在分析的每一步中，都要为每一个单独的事件重新计算描述粒子行为的复杂数学模型。随着事件数量增长到数十万甚至数百万，这一过程便成为了瓶颈，消耗大量的计算机内存并耗费数天甚至数周的时间。计算机内存会被中间结果填满，而主处理器与图形卡之间不断的数据传输也使工作进度变得极其缓慢。这种局限性意味着，一些来自现代实验的最具前景的高统计量数据无法得到充分利用，使得发现新奇异物质的可能性变得遥不可及。

为了打破这一障碍，由项本厚（Benhou Xiang）和方双实（Shuangshi Fang）领导的研究团队开发了一种名为 CTPWA 的新软件框架。该系统通过识别出大部分数学工作实际上并不需要重复进行，从而重新构思了分析的执行方式。在粒子衰变过程中，事件的几何结构以及粒子运动的方式是碰撞本身决定的固定事实；它们不会根据正在测试的特定共振态属性而改变。研究人员意识到，他们可以在最开始就将这些基于几何和运动的组件计算出来，并将其存储在计算机的高速内存中。在分析过程中，软件只需更新与共振态属性相关的微小变化部分，而无需为每个事件从头开始重建整个图像。

通过这种方式组织数据，团队能够将整个计算过程转移到图形处理器（GPU）上——这种专门用于渲染视频游戏的芯片，在处理并行数学任务方面表现卓越。由于繁重的几何预计算工作已在前期完成，GPU 可以利用高效的矩阵运算同时处理数百万个事件，而不是逐一处理。这种方法消除了在主计算机和图形卡之间不断移动数据的需求，而这正是以往方法中主要的延迟来源。结果表明，该系统不仅速度更快，而且内存效率更高，能够处理以往会导致软件崩溃的数据集。

团队使用一个特定的粒子衰变模拟（$J/\psi \to K^+K^-\eta$）对这一新框架进行了测试，该模拟涉及一百万个事件。在此测试中，他们模拟了一个场景：允许十种不同的中间共振态属性发生变化，以寻找最佳拟合。使用他们的新方法，包括确定结果不确定性所需的复杂计算在内的整个分析过程，在大约 100 秒内便完成了。相比之下，使用标准图形处理器的现有最佳方法需要花费显著更长的时间，而运行在标准中央处理器上的方法则要慢上数千倍。在某些情况下，传统方法由于内存耗尽而无法完成任务，而新框架却能轻松处理完整的百万级事件数据集。

研究结果表明，这种新方法使得对现代实验产生的大规模数据集进行高精度分析变得可行。通过将这些复杂拟合所需的时间从数天缩短至数分钟，该框架为更详细的强子谱学研究打开了大门。这种效率使物理学家能够更彻底地搜寻奇异粒子（如胶球或混合介子）的微妙信号，这些粒子长期以来一直被理论所预言，但却难以证实。研究人员已将其软件公开发布，为科学界提供了一个探索低能区强相互作用有效自由度的实用工具。凭借这种能力，该领域可以超越过去的限制，将高亮度实验产生的海量数据转化为一张更清晰的亚原子世界图谱。