← 返回 PaperDaily
大模型与智能体
BESIII首次看到粲偶素双奇异末态，分支比仅2.84×10⁻⁶
一个粲偶素粒子“临终”前，罕见地吐出了两代奇异粒子的组合。BESIII在27.12亿个事例里翻箱倒柜，终于以11.5σ的把握抓住了这个概率只有百万分之几的衰变。
龙哥读论文
阅读 1
查看原文
原论文信息如下：
在进入今天的主角之前，不妨先想一个问题：如果给你27.12亿张彩票，要你在里面找出几十张“隐藏款”，每张都必须同时满足四五个苛刻条件，你会不会觉得头大？
这恰恰就是高能物理实验家常年的日常。只不过他们手里的彩票，是一堆正负电子对撞后产生的粒子事例；而他们要寻找的“隐藏款”，则是一些极其稀有的粒子衰变过程。
最近的BESIII合作组就干成了一件这样的事：他们在(2.712±0.014)×10⁹个ψ(3686)粒子衰变事例中，首次清清楚楚地看到了ψ(3686)→pK⁻K_S⁰Ξ̄⁰+c.c.这个衰变过程，统计显著性高达11.5σ。
11.5σ是什么概念？在高能物理里，发现一个新粒子或新过程，通常要5σ才算“官宣”。5σ意味着这件事是随机涨落的概率大约只有三百五十万分之一。而11.5σ，基本上属于“把实锤焊死了”的级别。
更妙的是，这个过程的末态非常有意思：一个质子p、一个负K介子K⁻、一个中性K介子K_S⁰，外加一个双重奇异重子Ξ̄⁰。翻译成大白话就是：粲偶素在“临终”前，一口气吐出了一堆带奇异数的产物。这种“双奇异”加“双重奇异重子”的组合，在粲偶素衰变里是出了名的稀罕货。
那这道题到底难在哪？BESIII又是怎么从二十多亿个事例里把这不到一百个“幸运儿”挑出来的？这篇文章就带大家把整个破案过程捋一遍。
先认识一下“嫌疑人”：ψ(3686)到底是什么？
ψ(3686)这名字听起来像某种暗号，其实它是一个粲偶素态，可以粗略理解成一个粲夸克c和一个反粲夸克c̄抱在一起组成的束缚态。粲夸克是第二代夸克，质量比上夸克、下夸克大得多，所以c-c̄这个“夸克对”表现得像个小小的重原子系统，非相对论量子力学也能把它描述得八九不离十。
这种粲偶素粒子是不稳定的，它随时可能衰变。按照微扰量子色动力学（QCD）的标准剧本，ψ(3686)衰变的主要方式，是里面的c和c̄先湮灭，变成三个胶子或者一个虚光子，然后再“重组”出各种强子。绝大多数已知的ψ(3686)衰变末态里，最多只带一个s-s̄夸克对。
为什么带两个s-s̄对就这么难？
因为奇异夸克s比u、d夸克重，产生一个s-s̄对需要额外能量。在粲偶素衰变的能量预算里，想同时凑出两对奇异夸克，就像让你同时背两袋大米上六楼——不是不可能，但确实累得够呛。因此，这类“多奇异”衰变过程极其稀少，但恰恰因为它们稀少，才对粲偶素的衰变机制格外敏感，是检验理论模型的绝佳探针。
论文盯上的这个ψ(3686)→pK⁻K_S⁰Ξ̄⁰，就是典型代表。你看这个末态：K⁻和K_S⁰都含奇异夸克，Ξ̄⁰是ss双奇异重子，相当于末态里同时出现了“两份”奇异夸克来源。过去BESIII曾观察过同类的ψ(3686)→Λ(1520)(→pK⁻)K⁻Ξ̄⁺衰变，分支比在10⁻⁶量级。而这次这个新通道，跟它是同位旋对称性的“亲戚”，理论上也该能看见，只是谁也没想到，实验上真正抓到它，需要等到积累了27亿个ψ(3686)事例之后。
27亿个事例是怎么攒出来的？
想要看清楚ψ(3686)那点“家底”，首先得有足够多的ψ(3686)。这里的样本来自北京正负电子对撞机II（BEPCII）和BESIII探测器。BEPCII是一个运行在质心能量1.84到4.95 GeV区间的对撞机，峰值亮度达到1.1×10³³ cm⁻²s⁻¹。简单说，它就像一条繁忙的“粒子高速公路”，让电子和正电子迎头相撞，然后产生各种感兴趣的粒子。
BESIII探测器则像一个巨大的“高速相机”，围绕对撞点层层包裹。它的主体是一个基于氦气的多层漂移室，用来精确测量带电粒子的动量和能量损失；外面套着飞行时间探测器，负责测粒子的飞行时间；再外面是CsI(Tl)电磁量能器，用来探测光子和电子；整个系统泡在1.0 T的超导螺线管磁场里，靠粒子在磁场里拐弯的曲率来算动量。这套组合拳打下来，带电粒子动量分辨率在1 GeV/c时能达到0.5%，相当能打。
首次发现ψ(3686)多奇异数衰变
在粒子物理里，能被称作“罕见”的过程很多，可这次BESIII合作组盯上的ψ(3686)→pK⁻KS0Ξ̄0+c.c.，即使在罕见事件里也算得上“有排面”：末态同时带出了两条奇异数线索——一个负K介子K⁻、一个中性K介子KS0，外加一个双重奇异重子Ξ̄0。这是粲偶素衰变中极难出现的多奇异组合。
为什么难？因为奇异夸克s比上夸克u和下夸克d都重。在粲偶素这类由正反粲夸克组成的束缚态衰变中，产生一对s-s̄夸克本来就要消耗额外能量，想同时凑出“两个奇异数来源”更是难上加难。过去人们熟悉的ψ(3686)衰变，末态里大多只含一个s-s̄对；而这次观测的末态，相当于要求粲偶素把一份并不富裕的能量预算，同时花在两条独立的奇异夸克产生路径上。
这类过程之所以让实验物理学家牵肠挂肚，是因为它提供了相当“挑剔”的探针。粲偶素内部的一对正反粲夸克湮灭后，是主要变成三个胶子，还是经过电磁过程变成虚光子？这些不同机制产生的多奇异强子末态比例并不相同。能够测量一个具体的多奇异衰变分支比，就能帮理论家看清楚：微扰量子色动力学（QCD，即描述夸克和胶子强相互作用的基本理论）在处理奇异夸克产生时，到底还有多少账对不上。
事实上，BESIII此前已经观察过与之同类的ψ(3686)→Λ(1520)(→pK⁻)K⁻Ξ̄⁺衰变，分支比在10⁻⁶量级。Λ(1520)是一个超子激发态，它衰变到pK⁻后，会留下一个跟本次通道非常相似的形状。理论上的同位旋对称性暗示，这次测的ψ(3686)→pK⁻KS0Ξ̄0也应该能出现，只是谁也没想到，要等到BESIII攒够27.12亿个ψ(3686)事例之后，它才肯露出真容。
这次的结果相当漂亮：信号产额87.1±11.4个，统计显著性11.5σ，分支比测得(2.84±0.40±0.25)×10⁻⁶。对一个总事例数高达27.12亿的样本来说，最终能被确认的净信号只有81.5个左右——换句话说，这真是一场不折不扣的“大海捞针”。
BESIII探测器与数据分析方法
像这种多奇异衰落道，实验上面临的第一关不是物理分析，而是如何从探测器信号里把粒子一个个认出来。BESIII探测器记录的是正负电子对撞后产生的各种粒子径迹和能量沉积，但实验真正关心的那些粒子往往“活不长”。比如Ξ̄⁰反超子，它一旦产生就会迅速衰变，根本不可能直接飞进探测器末端。实验上能做的，只能是先找到它衰变出的下一级粒子，再逐级向上重建。
这条重建链相当长。KS0要重建到π⁺π⁻；Ξ̄⁰则先衰变成Λ̄π⁰，其中的Λ̄再衰变到p̄π⁺，π⁰再衰变到两个光子γγ。也就是说，一个看似简单的末态pK⁻KS0Ξ̄0，在探测器里实际要识别的是p、K⁻、π⁺、π⁻、p̄以及两个光子。为了让这些粒子组合尽量干净，分析中还要做顶点拟合、飞行长度要求和衰变质量窗筛选。
最有杀伤力的筛选工具是一个五约束（5C）运动学拟合。为什么叫5C？因为其中四个约束来自能量-动量守恒，要把正负电子对撞的初态四动量与候选末态的全体动量精确对准；第五个约束来自π⁰衰变到γγ的质量约束，保证两个光子确实能组成一个π⁰。只有能够通过这套严格配对的候选者，才值得留下来做进一步分析。
光靠筛选还不够，因为探测器响应和效率必须靠蒙特卡洛（Monte Carlo，简称MC）模拟来校准。所谓MC模拟，就是按照已知物理过程在计算机里“虚拟地”重放一遍粒子衰变和探测器响应，看看一个真实信号会长成什么样。本文生成了一大批包容性MC样本，还额外制作了50万个混合信号MC事例。混合比例约4:1，一部分是直接的四体相空间衰变，另一部分则带有Λ(1520)中间共振成分，这个比例是通过检查数据里pK⁻不变质量谱确定的。
把数据和混合信号MC放到一起比较，就能直观检验模拟是否可靠。图2给出了多个两体不变质量谱的对比：带误差棒的黑点是ψ(3686)数据，蓝色实线是混合信号MC，黄色和粉色阴影分别代表来自不同本底源的估计。
图2：e⁺e⁻→pK⁻KS0Ξ̄0候选事例的两体不变质量分布。数据与混合信号MC样本总体符合良好。带误差棒的点表示质心能量√s=3.686 GeV数据，蓝色实线为混合信号MC；黄色阴影为包含式MC样本归一化本底，粉色阴影为KS0和Ξ̄0数据边带归一化本底。
数据点和蓝色信号MC在多个质量谱上都能对齐，本底则被压在很低的水平。这说明整套选择条件和MC模型是可靠的，可以放心进入下一步：把信号数量精确地“数”出来。
信号提取与分支比测量
要把“看起来有信号”变成“可以写进论文的信号数”，核心做法是二维无分箱最大似然拟合。所谓二维，是把每个候选同时放到由M(π⁺π⁻)和M(Λ̄π⁰)组成的平面上，信号事例会集中在两个不变质量的峰区重叠处。所谓无分箱，则是不对观测量做粗粒化分bin，把每个候选携带的细节信息都保留下来，在样本量不大时尤其有优势。
既然是二维拟合，本底也得按二维来分类。真正的信号要求KS0和Ξ̄0同时存在；可很多本底候选只有其中一个是真的，另一个是随机组合凑出来的。于是论文把本底拆成三类：含KS0但不含Ξ̄0的本底；含Ξ̄0但不含KS0的本底；以及KS0和Ξ̄0都不是真的组合本底。通过同时拟合信号与三类本底，就能尽量公平地把信号从“候选人堆”里分离出来。
拟合结果由图3给出。数据用带误差棒的点显示，蓝色实线是整体拟合结果。下方不同颜色的虚线分别代表不同类型的本底和信号成分：红色虚线是含KS0但非Ξ̄0的本底，绿色虚线是含Ξ̄0但非KS0的本底，蓝色虚线是信号，黄色虚线则是组合本底。
图3：对√s=3.686 GeV ψ(3686)数据样本中e⁺e⁻→pK⁻KS0Ξ̄0候选事例，在M(π⁺π⁻)与M(Λ̄π⁰)分布上的拟合投影。数据点为带误差棒的黑点，蓝色实线为拟合结果，各色虚线表示不同信号与本底成分。
最终，ψ(3686)→pK⁻KS0Ξ̄0在3.686 GeV数据中的信号产额被确定为87.1±11.4。统计显著性用最大似然比来判定：把包含信号成分和不包含信号成分的两个拟合的最大似然值做差，得到Δ(lnL)，再换算成正态分布的标准差数。本文算出的结果是11.5σ——这意味着信号来自随机涨落的可能性微乎其微，属于“实锤中的实锤”。
不过，87.1个信号并不是全部来自ψ(3686)衰变，还要小心扣除一种“冒牌货”：连续本底。正负电子对撞不一定要先形成ψ(3686)，也可以直接通过连续的e⁺e⁻湮灭过程产生强子。为了估计这部分贡献，论文使用了质心能量√s=3.773 GeV的离共振数据样本。这个能量点远离ψ(3686)共振峰，产生的几乎全是连续强子本底，因此可以作为“参照组”。
图4展示的就是这个3.773 GeV参照样本的拟合投影。可以看出，该样本中也有少量候选事件，需要通过亮度、能量依赖因子和探测效率换算到3.686 GeV能量点。
图4：对√s=3.773 GeV离共振数据样本，在M(π⁺π⁻)与M(Λ̄π⁰)分布上的拟合投影。带误差棒的点为数据，各色曲线含义与图3一致。该样本主要用于估计连续本底贡献。
扣完连续本底后，净信号产额Nnet=81.5±11.5。有了净信号产额，分支比就可以用下面这个核心公式来计算。式子里的Nψ(3686)是总的ψ(3686)事例数，ε是探测效率，后面四个分支比的乘积来自KS0、Ξ̄0、Λ̄和π⁰各自中间衰变道的已知分支比乘积。
公式（1）：ψ(3686)→pK⁻KS0Ξ̄0的分支比计算公式。Nnet为扣除连续本底后的净信号产额，Nψ(3686)为数据中ψ(3686)总数，∏iBi是四个中间衰变分支比的乘积，ε为信号MC给出的探测效率。
连续本底从3.773 GeV换算到3.686 GeV时，需要引入一个标度因子fc。净信号产额Nnet与标度因子fc的计算方式如下：
公式（2）：净信号产额Nnet等于3.686 GeV数据中拟合出的信号产额N3.686，减去折算后的连续本底产额Ncontinuum,3.686。
公式（3）：标度因子fc由3.686 GeV与3.773 GeV两个能量点的积分亮度比、能量依赖因子以及连续过程探测效率比共同决定。文中取能量依赖幂指数n=1，对应1/s型的截面行为。
把各数值代入后，最终得到的探测效率为2.43%，分支比为(2.84±0.40)×10⁻⁶。表1汇总了信号产额、连续本底估计、探测效率以及分支比等关键量。表中误差只含统计部分。
表1：分支比测量中的关键数值汇总。N3.686=87.1±11.4是ψ(3686)数据中的拟合信号产额；Ncontinuum,3.773=31.2±6.6是3.773 GeV离共振样本中的产额；换算到3.686 GeV后连续本底为5.6±1.2；净信号产额Nnet=81.5±11.5；探测效率ϵ=2.43%；分支比B=(2.84±0.40)×10⁻⁶。
除了统计检验，论文还用50个统计独立的伪实验对整个分析流程做了验证。每个伪实验都按真实数据的样本量抽样，再走一遍和真实数据完全相同的分析流程。结果平均测量值与输入值一致，说明这套信号提取方法本身没有明显的系统偏向。
系统误差的精细评估
分支比测量不能只看统计误差。所谓系统误差，是指由探测器校准、事例选择、模拟模型等环节的不完美带来的偏差，它不会因为样本量增大而自动消失，必须逐项盘查。
第一类误差来自带电粒子重建。质子与K⁻的径迹重建和粒子鉴别（PID）误差，分别通过J/ψ→pp̄π⁺π⁻等对照样本来标定，得到数据与MC之间约1%的差异。由于末态里有一个质子和一个K⁻，径迹和PID两项误差合起来各占2.0%。KS0重建的研究使用J/ψ→K*±K̄∓等控制样本，误差1.5%；Λ重建使用J/ψ→pK⁻Λ等控制样本，误差2.0%；π⁰重建使用D⁰介子衰变控制样本，误差2.0%。
第二类误差来自分析方法和MC模型本身。二维拟合的误差由信号形状变化、本底函数阶数改变以及拟合范围变化三部分合成，共2.0%；MC模型误差则通过比较名义信号样本与纯相空间样本的效率差异来估计，为2.9%；5C运动学拟合的螺旋参数修正带来1.6%的误差；信号MC样本量有限也会贡献0.9%的统计误差。
第三类误差相对隐蔽，是共振振幅与连续本底振幅之间可能的干涉效应。在ψ(3686)共振峰附近，信号产生的振幅和连续强子产生的振幅会像两条波一样叠加：合起来的强度不只是简单求和，还取决于两者的相对相位。如下式所示，总截面由连续振幅ac与共振振幅aR的相干叠加决定。
公式（4）：总截面σtotal的表达式。ac是连续振幅，aR是共振振幅，φ是二者的相对相位。若相位未知，干涉可能增强也可能压低信号，进而影响分支比的估计。
公式（5）：干涉对共振分支比的相对影响rR。A由连续截面与分支比决定，B来自共振参数，φ仍是相对相位。由于现有离共振数据量不足以确定φ，论文假设φ在0到360度范围内均匀分布，再保守地取边界值，最终把这项干涉效应定为6.3%的系统误差。
公式（6）：连续过程截面σc的计算式。Nc是离共振数据中拟合得到的产额，L为积分亮度，ε为连续过程探测效率，再乘上能量依赖因子换算到3.686 GeV。
把所有来源按平方和方式合成，最终的系统误差为8.7%。表2给出了完整的误差清单。可以看到，最大的单项误差正是上面讨论的潜在干涉效应6.3%，其次是MC模型2.9%，其余各项大多在2%左右。这说明该结果已经做到了实验条件允许范围内的相当精细。
表2：分支比测量中各项系统误差来源及相对大小。总计为8.7%。其中“潜在干涉”是最大贡献项，来自共振与连续振幅之间的相位不确定性。
物理意义与未来展望
从物理上看，这次测量有两层价值。第一层，它首次给出了一个“两个奇异介子+一个双奇异重子”的ψ(3686)衰变分支比，实验数值直接为多奇异粲偶素衰变机制提供了参照点。第二层，在pK⁻不变质量谱中可以看到Λ(1520)贡献的迹象。Λ(1520)这个超子激发态的产生，与之前测过的K⁻Λ(1520)Ξ̄⁺通道通过同位旋对称性相连。也就是说，这次的新结果等于在另一条同位旋通道上又钉下了一颗“实验钉子”。
当然，论文并没有声称已经完全理解了衰变机制。当前数据量下，能看到Λ(1520)的迹象，但还不足以对中间共振结构做完整的振幅分析。干涉相位也仍然受限于离共振样本的统计量。这些悬而未决的问题，正是未来实验继续发力的方向。
下一步值得期待的是未来超级tau-粲工厂（Super Tau-Charm Factory）。那类设施能提供比现在高得多的ψ(3686)统计量，到时不仅可以对这种多奇异衰变道做精细的中间共振研究，还有可能发现新的强子态结构。对普通读者来说，这次测量真正迷人的地方在于：即便标准模型已经相当成功，强相互作用的非微扰细节里仍然藏着许多等待实验去锚定的“暗礁”。每多一个稀有衰变道的精确测量，就多一块帮理论模型校准航向的基石。
龙迷三问
这篇论文到底在解决什么问题？BESIII实验利用27.12亿ψ(3686)事例，首次观测到ψ(3686)→pK⁻K_S⁰Ξ̄⁰+c.c.衰变，统计显著性11.5σ，分支比(2.84±0.40±0.25)×10⁻⁶，为多奇异粲偶素衰变研究提供关键实验参考。
这篇工作最值得看的点是什么？首次观测到ψ(3686)→pK⁻Kₛ⁰Ξ̄⁰+c.c.衰变，统计显著性达11.5σ，分支比测量为(2.84±0.40±0.25)×10⁻⁶，系统误差控制良好（总系统不确定度8.7%）。
这篇工作的边界或风险在哪里？优点：实验方法严谨，利用大规模数据样本和精细的事例选择标准，通过二维拟合有效分离信号与背景，系统误差评估全面细致。缺点：连续谱背景贡献的估计依赖于模型假设（n=1），存在一定不确定性；Λ(1520)中间共振态的贡献尚未精确确定。
如果你还有哪些想要了解的，欢迎在评论区留言或者讨论~
龙哥点评
论文创新性分数：★★★★☆
这是一项实验发现型成果，首次观测到ψ(3686)衰变到“两个奇异介子+一个双奇异重子”的末态，并给出首个分支比测量。方法层面沿用BESIII成熟的粒子重建与拟合框架，创新更多体现在物理结果的首创性上。
实验合理度：★★★★★
分析流程完整且扎实。KS0、Λ、π⁰等重建效率都使用独立控制样本标定，连续本底扣除使用离共振数据，还做了50个伪实验验证无偏性，各项系统误差逐一量化，实验设计几乎没有明显漏洞。
学术研究价值：★★★★☆
多奇异粲偶素衰变是检验强相互作用产生机制的重要窗口。这个分支比可以作为理论模型的新约束，也为未来超tau-粲工厂的相关研究提供了实验坐标。不过它属于特定物理问题的积累，不改变粒子物理整体框架。
稳定性：★★★★☆
信号统计显著性高达11.5σ，核心结果非常稳固。但信号产额本身只有87个，连续本底干涉又贡献了6.3%的系统误差，说明结果在统计和系统两个层面上仍有可见的“窄缝”。相对大部分高能物理稀有衰变测量，这已经属于高稳定性结果。
适应性以及泛化能力：★★★☆☆
分析方法可以推广到BESIII数据中其他含KS0、Λ或Ξ的超子稀有衰变测量，在这个范围内泛化能力不错。但它不是通用方法，离开高能对撞实验环境后并不适用。
硬件需求及成本：★★☆☆☆
本工作依赖北京正负电子对撞机BEPCII和BESIII谱仪这类大科学装置，建设和运行成本高昂。好在数据已经采集完成，单次分析本身的算力需求相对可控，主要成本集中在装置维护与大规模MC模拟上。
复现难度：★★☆☆☆
对粒子物理以外的小组来说，几乎无法复现，因为需要专用加速器、探测器、庞大的MC样本和整套离线软件框架。论文分析流程描述相当细致，如果具备BESIII合作组的数据使用条件，可以按相同步骤验证。
产品化成熟度：★☆☆☆☆
纯基础科学研究，不面向产品化。其产出是物理知识和实验基准，短期不会有直接的技术产品转化。
可能的问题：信号产额仅87个，即使显著性很高，统计样本仍偏小；pK⁻质量谱中Λ(1520)的贡献只能作为迹象，尚未做分波或振幅分析；共振与连续本底干涉的相对相位无法由现有离共振数据确定，只能假设均匀先验做保守估计。更细致的机制结论有待未来更高统计量数据。
[1] M. Ablikim et al. (BESIII Collaboration), Observation of ψ(3686) → pK−K0S Ξ¯0 + c.c., arXiv:2609.02051, https://arxiv.org/pdf/2609.02051v1.pdf
微观世界越挖越“奇异”，BESIII这回把双奇异末态都翻出来了！想第一时间看懂高能对撞实验里的新信号？来『龙哥读论文』粉丝群，和大伙一起聊聊“多少个σ才算发现”那些事儿。欢迎加入龙哥读论文粉丝群，
扫描下方二维码或者添加龙哥助手微信号加群：kangjinlonghelper。
一定要备注：研究方向+地点+学校/公司+昵称（如 粒子物理+北京+中科院高能所+小龙人），根据格式备注，可更快被通过且邀请进群。
*本文仅代表个人理解及观点，不构成任何论文审核或者项目落地推荐意见，具体以相关组织评审结果为准。欢迎就论文内容交流探讨，理性发言哦～ 想了解更多原文细节的小伙伴，可以点击"阅读原文"，查看更多原论文细节哦！