← 返回 PaperDaily
大模型与智能体
BESIII首次观测Σ⁰Σ̄⁰η衰变,12%规则再遭检验
粒子物理界的“寻宝游戏”又有新斩获!BESIII合作组利用百亿级数据样本,首次观测到J/ψ和ψ(3686)衰变到Σ⁰Σ̄⁰η的稀有过程,直接检验粲偶素衰变中的12%规则。想了解粒子如何“变装”成超子对加介子?这篇带你一探究竟。
龙哥读论文
发布于 2026-08-20 00:20:08
阅读 2
查看原文
原论文信息如下:
在粒子物理的世界里,有一个非常有趣的“家族遗传定律”——12%规则。它预言,同一家族的两个粒子J/ψ和ψ(3686)(粲偶素家族的兄弟俩),衰变到完全相同末态的分支比,应当大致等于12%。这一定律源于微扰量子色动力学(QCD)的理论预言,但实验上却发现了不少“叛逆者”,其中最著名的就是“ρπ之谜”:ψ(3686)衰变到ρπ被强烈抑制,与12%规则的预言差了将近两个数量级!
为什么有的衰变道乖乖听话,有的却“叛逆”?这正是研究粲偶素强子衰变的核心问题之一。BESIII合作组这次把目光锁定在了一个此前从未被观测过的末态——Σ⁰Σ̄⁰η,用百亿级的J/ψ数据和数十亿级的ψ(3686)数据,首次“抓住”了这两个稀有衰变过程,并对12%规则进行了新一轮检验。
引言:12%规则与粲偶素衰变
要理解这项工作的意义,得先从粲偶素说起。粲偶素(charmonium)是由一个粲夸克(c)和一个反粲夸克(c̄)组成的束缚态,J/ψ和ψ(3686)就是粲偶素家族里最出名的两位。它们都是自旋为1的矢量介子,质量分别为3.097 GeV/c²和3.686 GeV/c²,就像是同一对父母生下的哥哥和弟弟。
按照微扰QCD的理论预言,ψ(3686)衰变到某个末态h的分支比与J/ψ衰变到同一末态的分支比之比,应该约等于二者到正负电子对分支比之比,即约12%。这一定律可以用公式简洁表达:Q_h = B(ψ(3686)→h) / B(J/ψ→h) ≈ 12%。在QCD的框架下,这个比值之所以成立,是因为两个粒子衰变到同一强子末态的机制相似,只是相空间和波函数重叠有所差异。
然而实验发现,不少衰变道并不遵守这条规则。1983年,Mark II实验首先发现ψ(3686)→ρπ相对于J/ψ→ρπ被强烈抑制,形成了著名的“ρπ之谜”。后续实验又陆续发现了一些偏离:比如J/ψ和ψ(3686)衰变到ΛΛ̄π⁰这个同位旋破坏过程,与12%规则有明显偏离;而衰变到ΛΛ̄η(同位旋守恒过程)却与规则一致。这种“有的道符合、有的道偏离”的现象,让物理学家意识到12%规则并不是铁律,其偏离机制可能蕴藏着超越微扰QCD的新物理线索。
正是这些未解之谜,推动着BESIII合作组不断寻找新的衰变道来检验12%规则。本论文关注的是J/ψ和ψ(3686)衰变到Σ⁰Σ̄⁰η的强子过程。Σ⁰是含奇异夸克的超子,η是赝标量介子,这个末态属于“重子-反重子-赝标量介子”(BB̄P)家族,是一类能够用来研究激发态重子(如Σ*)和重子-反重子阈值增强现象的理想过程。BESIII拥有世界上最大的J/ψ和ψ(3686)数据样本,正是寻找这类稀有衰变的最佳场所。
方法概述:BESIII探测器与数据样本
BESIII探测器运行在北京正负电子对撞机BEPCII上,这个“粒子工厂”能在1.84到4.95 GeV的质心系能量范围内提供高亮度的正负电子对撞。BESIII探测器的主体是一个圆柱形结构,覆盖了93%的立体角,内部包括多层漂移室(MDC)用于测量带电粒子的动量、飞行时间探测器(TOF)用于粒子鉴别、CsI(Tl)电磁量能器(EMC)用于测量光子和电子能量,整个系统置于1.0 T的超导螺线管磁场中。动量分辨率在1 GeV/c时达到0.5%,光子能量分辨率在桶部为2.5%,这些性能指标为精确测量稀有衰变提供了坚实基础。
本工作使用的数据量非常惊人:(10087±44)×10⁶个J/ψ事件和(2712±14)×10⁶个ψ(3686)事件。这些数据是BESIII合作组多年运行积累的成果,也是世界上最大的粲偶素数据样本。正是依托如此庞大的统计量,研究人员才有能力去寻找分支比低至10⁻⁵量级的稀有衰变过程。
核心设计:事件选择与信号重建
要在海量数据中找到目标信号,无异于大海捞针。研究团队设计了一套完整的事件选择流程,逐层筛选出候选事件。
首先看末态粒子的衰变链。Σ⁰和反Σ⁰超子分别通过γΛ和γΛ̄衰变重建,η介子通过γγ衰变重建,而Λ和Λ̄又通过pπ⁻和p̄π⁺衰变重建。因此最终需要重建的末态是:p、p̄、π⁺、π⁻和四个光子。整个过程就像剥洋葱一样,从最外层的稳定粒子逐层向内还原出短寿命的中间态。
在带电径迹的选择上,要求极角满足|cosθ|<0.93,并且离对撞点的最近距离在横向小于10 cm、纵向小于20 cm。粒子鉴别(PID)联合了漂移室的能量损失(dE/dx)和飞行时间(TOF)信息,为每一条径迹分配最可能的粒子类型。光子候选则要求能量在桶部大于25 MeV、端盖大于50 MeV,并且时间信息与事件起始时间一致。
Λ和Λ̄的寻找是关键的中间步骤。研究者将带正电的质子候选与带负电的π介子候选组合,进行二次顶点拟合——因为Λ是长寿命粒子,会飞行一段距离后才衰变,次级顶点与主对撞点有明显分离。每个事件选择χ²最小的Λ和Λ̄候选。
接下来是运动学拟合。研究者对ΛΛ̄γγγγ假设做四约束(4C)运动学拟合,要求所有末态粒子的总四动量与初始的正负电子对撞系统能量守恒。对于有多余光子组合的事件,保留χ²最小的组合。J/ψ和ψ(3686)分别要求χ²<60和χ²<40,这些阈值通过最大化S/√(S+B)来优化。
Σ⁰、Σ̄⁰和η候选通过质量窗口选择:|M(γΛ)−m_Σ⁰|<15 MeV/c²、|M(γΛ̄)−m_Σ̄⁰|<15 MeV/c²、|M(γγ)−m_η|<25 MeV/c²。为了压低含π⁰的背景,要求|M(γγ)−m_π⁰|>20 MeV/c²。对于ψ(3686)样本,还额外要求|M(Σ⁰Σ̄⁰)−m_J/ψ|>30 MeV/c²,以剔除ψ(3686)→ηJ/ψ、J/ψ→Σ⁰Σ̄⁰的级联衰变背景。对于J/ψ样本,则额外做一次ΛΛ̄γγγ假设的4C拟合,以剔除少一个光子的背景过程。
背景研究利用包含100亿J/ψ事件和27亿ψ(3686)事件的inclusive MC样本。在γγ不变质量谱上没有观察到显著的峰状背景,因此采用Σ⁰Σ̄⁰二维边带法来估计组合本底。边带区域取在[1.1475, 1.1625] GeV/c²和[1.2225, 1.2375] GeV/c²,通过边带数据的比例因子来推算信号区内的本底贡献。
核心原理推导:分支比测量与效率修正
分支比的测量公式是:B(ψ′→Σ⁰Σ̄⁰η) = N_sig / (N_ψ′ × ε_w × B₁² × B₂² × B₃),其中N_sig是信号产量,N_ψ′是J/ψ或ψ(3686)的总事件数,ε_w是加权探测效率,B₁、B₂、B₃分别是Σ⁰→γΛ、Λ→pπ首次观测!J/ψ和ψ(3686)衰变到Σ⁰Σ̄⁰η的新发现
在粒子物理实验的世界里,“首次观测”这四个字的分量极重——它意味着信号显著性必须超过5个标准偏差(5σ),远超“证据”的门槛。BESIII合作组这次的成果,就是同时把J/ψ→Σ⁰Σ̄⁰η和ψ(3686)→Σ⁰Σ̄⁰η两个过程都推过了这条线。
研究人员在(10087±44)×10⁶个J/ψ事件和(2712±14)×10⁶个ψ(3686)事件中寻找信号。这个末态的衰变链相当“绕”:Σ⁰会迅速衰变成γΛ,反Σ⁰衰变成γΛ̄,η衰变成γγ,而Λ/Λ̄再衰变成pπ⁻/p̄π⁺。所以探测器最终看到的是一堆质子、反质子、π介子和光子的组合。要在海量数据中把这种组合精确拼装出来,靠的是严格的事件选择加运动学拟合,每一步都像在玩一场高难度的拼图游戏。
在完成初选之后,研究者对γγ不变质量谱在η信号区附近做拟合,可以清晰地看到η介子产生的信号峰。图2展示了J/ψ样本(上图)和ψ(3686)样本(下图)中M(γγ)的拟合结果,带误差棒的黑点是数据,蓝色实线是总拟合结果,红色虚线是信号分量,蓝色虚线是平滑本底。
B(J/ψ→Σ⁰Σ̄⁰η) = (7.5 ± 0.3 ± 0.8) × 10⁻⁵
B(ψ(3686)→Σ⁰Σ̄⁰η) = (1.3 ± 0.1 ± 0.1) × 10⁻⁵
其中第一个不确定度为统计不确定度,第二个为系统不确定度。从数值上可以直观感受到,J/ψ衰变到Σ⁰Σ̄⁰η的概率大约是十万分之七点五,而ψ(3686)衰变到同一末态的概率大约是十万分之一多一点。这些过程属于典型的稀有衰变,如果没有百亿量级的“家底”,根本不可能找到它们。
12%规则再受检验:Q值17.3%意味着什么?
有了两个分支比,自然就可以计算它们的比值Q。根据论文中的定义:
其中B(ψ(3686)→h)和B(J/ψ→h)分别表示ψ(3686)和J/ψ衰变到同一末态h的分支比。这个比值预期约为12%,因为它大致等于二者到正负电子对分支比的比值。12%规则的物理根源在于微扰QCD:当粲偶素衰变到强子末态时,如果不考虑非微扰效应的细节差异,ψ(3686)和J/ψ衰变到同一末态的部分宽度之比应当与到轻子对的部分宽度之比一致。
Q = B(ψ(3686)→Σ⁰Σ̄⁰η) / B(J/ψ→Σ⁰Σ̄⁰η) = (17.3 ± 1.5 ± 1.7)%
这个17.3%比12%高出了一截,看起来似乎有“超标”的嫌疑。但物理学家判断“符合”还是“偏离”,不能只看中心值,还得看误差。把统计不确定度1.5%和系统不确定度1.7%相加,总不确定度大约2.3%。也就是说,17.3%和12%之间的差距大约是5.3个百分点,换算成标准偏差大约是2倍多。论文原文给出的是“与12%规则在3.0σ内一致”,其实就是说:虽然中心值看着偏高,但误差还挺大,眼下还不能断定规则被打破。
这个结果放在粲偶素衰变的大图景里非常有意思。此前实验发现,J/ψ和ψ(3686)衰变到ΛΛ̄η(同位旋守恒过程)时Q值符合12%规则;而衰变到ΛΛ̄π⁰(同位旋破坏过程)时则明显偏离。如今Σ⁰Σ̄⁰η给出的中心值17.3%又比12%高了一些,虽然统计上还不显著,但方向上似乎暗示着不同超子末态之间的动力学机制存在细微差别。
除了测量分支比,论文还专门检查了Σ⁰η、Σ̄⁰η和Σ⁰Σ̄⁰的不变质量谱,想看看有没有中间共振态或者阈值增强结构的蛛丝马迹。结果如图3所示,数据分布的形状和纯相空间(PHSP)蒙特卡罗模拟基本一致,没有观察到显著的峰状结构。这说明该三体衰变中中间态动力学很可能不占主导,或者说其贡献在当前统计精度下还分辨不出来。
虽然没有看到新粒子“现身”,但这类搜索本身就是对理论模型的重要约束。不同的理论框架(如胶子机制、碎裂机制、重子偶素等)对阈值增强结构的预言不同,多一个通道的“干净”结果,就多一份排除或限制的判据。
BESIII探测器:世界最大粲偶素数据样本的采集者
搞出这么大动静的BESIII探测器,到底什么来头?BESIII全称是北京谱仪III(Beijing Spectrometer III),运行在北京正负电子对撞机BEPCII上。BEPCII是一个质心系能量覆盖2.0到4.95 GeV的“粒子工厂”,专门在粲偶素能区附近提供高亮度的正负电子对撞。所谓“粒子工厂”,就是能像工厂流水线一样大批量生产某种粒子的对撞机——BEPCII在J/ψ能区的峰值亮度达到了1×10³³ cm⁻²s⁻¹量级,这也是为什么它能积累出全球最大粲偶素数据样本的根本原因。
BESIII探测器的主体是一个大型圆柱形装置,从内到外主要包含几个关键子系统。最内层是多层漂移室(MDC),负责测量带电粒子的径迹和动量,在1 GeV/c的动量下分辨率约为0.5%。紧挨着的是飞行时间探测器(TOF),利用粒子飞行时间差异来做粒子鉴别(PID),时间分辨率约为80皮秒(桶部)。再往外是CsI(Tl)晶体电磁量能器(EMC),用来测量光子和电子的能量,在桶部对1 GeV光子的能量分辨率约为2.5%。最外层是μ子探测器,用来鉴别μ子。整个探测器置于一个1.0 T的超导螺线管磁场中,实现对带电粒子动量的精确测量。
这些性能参数意味着什么?简单说,BESIII能同时精确重建带电粒子和中性粒子(光子),而且粒子鉴别能力强、能量分辨率高。对于Σ⁰Σ̄⁰η这种包含多个光子和多个超子的复杂末态来说,探测器的hermeticity(立体角覆盖达93%)和能量分辨率就是生命线——覆盖不够全,光子容易“漏”掉;分辨率不够好,η峰就被背景淹没。
本工作使用的数据量是(10087±44)×10⁶个J/ψ事件和(2712±14)×10⁶个ψ(3686)事件。这组数字是BESIII合作组多年运行积累的成果,也正是因为这组世界最大样本的存在,研究者才有底气去挑战分支比低至10⁻⁵量级的稀有过程。换一个探测器,换一个小一号的样本,这事儿还真干不成。
从数据到结论:事件选择与分支比测量的关键技术
说了半天“事件选择”,那到底是怎么从海量数据里把信号挑出来的?这里面的技术含量,值得展开讲一讲。
首先看衰变链的“剥洋葱”策略。Σ⁰是短寿命粒子,在探测器里根本没法直接看到,只能通过它的衰变产物γΛ来重建。Λ虽然也是短寿命粒子,但它有个特点——飞行距离比较长(cτ约7.9 cm),所以Λ→pπ⁻的衰变顶点可以和主对撞点明显分开。利用这个特征,研究者对pπ⁻组合做“二次顶点拟合”,把Λ候选的飞行信息也用上,能大幅压低组合本底。η介子则通过γ→γγ的两个光子重建。这样一来,整个事件的末态就是p、p̄、π⁺、π⁻和四个光子。
接下来是运动学拟合。研究者把“ΛΛ̄γγγγ”的假设做四约束(4C)运动学拟合,其中“4C”指的是能量守恒加上三维动量守恒这四个约束条件。拟合的目标是让所有末态粒子的总四动量与初始正负电子对撞系统的四动量严格匹配,同时允许对粒子的动量和能量做微小调整。对于那些有多余光子组合的事件,保留拟合χ²最小的组合。J/ψ和ψ(3686)样本分别要求χ²<60和χ²<40,这个阈值不是拍脑袋定的,而是通过最大化S/√(S+B)(信号显著性)来扫描优化的。
在初选完成后,下一个关键步骤是计算探测效率。效率不是简单的一个数,而是逐事件加权得到的。加权效率的定义为:
其中wᵢ是每个蒙特卡罗(MC)产生事件对应的权重。之所以要加权,是因为真实的信号过程可能存在中间态或者特定的角分布,而简单的相空间MC不一定能完全描述这些特征。研究者会根据数据中的分布(比如Σ⁰η不变质量谱)对MC样本重新加权,让MC与数据在运动学变量上尽可能一致,从而减小因MC模型不准确带来的效率偏差。
公式中各项的物理含义是:N_sig是拟合得到的信号事件数;N_ψ′是J/ψ或ψ(3686)的总事件数;ε_w是前面说的加权探测效率;B₁是Σ⁰→γΛ的分支比;B₂是Λ→pπ⁻的分支比;B₃是η→γγ的分支比。由于Σ⁰和Λ在末态中各出现两次(正反各一次),所以B₁和B₂都取平方。把这些因子全部除掉,才能从“看到的事件数”反推回“真实产生的概率”。
背景压制方面,论文采用了两条腿走路的策略。一方面用包含百亿级事件的inclusive MC样本研究可能的物理背景来源;另一方面对γγ不变质量谱进行拟合,用“平滑本底+信号峰”来提取信号产额。论文还特别检查了Σ⁰Σ̄⁰的二维边带区域,确保没有漏掉的峰状背景。对于ψ(3686)样本,还额外剔除了ψ(3686)→ηJ/ψ、J/ψ→Σ⁰Σ̄⁰的级联衰变污染——这类背景在运动学上非常接近真实信号,如果不小心处理,会把分支比测高。
系统误差全面评估:测量结果的可靠性与精度
分支比的测量精度不仅取决于统计量,还取决于系统误差的控制。这篇论文在系统误差评估上做得很细,几乎把每一项可能的偏差来源都列了个遍。表I汇总了所有相对系统不确定度的来源。
从表中可以看到,系统误差的来源主要包括以下几个方面。角分布不确定性 :信号产生过程的角分布可能偏离MC模拟的假设,导致效率计算出现偏差,论文通过改变产生过程的模型来估计这部分不确定度。质量窗口选择 :Σ⁰和η的信号窗口宽度变化会影响最终产额,论文通过改变窗口大小来评估。此外还有背景形状参数化 、中间态分支比(B₁、B₂、B₃)的世界平均值精度 、产额拟合方法 和总事件数N_ψ′的标定误差 等。
具体到数值,J/ψ道的总系统误差约为10.7%(0.8/7.5),ψ(3686)道的总系统误差约为7.7%(0.1/1.3)。J/ψ道系统误差偏大,主要是因为该道的修正因子(角分布等)不确定性更高。相比统计误差(J/ψ道约4%,ψ(3686)道约7.7%),系统误差在J/ψ道占主导,在ψ(3686)道则与统计误差相当。
从整体来看,这个系统误差体系是成熟且可信的。BESIII合作组在粲偶素物理上有多年积累,对这套误差评估方法已经打磨得相当完善。值得一提的是,Q值测量中许多系统误差在比值中会部分抵消(比如探测效率相关的误差),所以Q值的总系统误差(约1.7%)相对于单个分支比的系统误差来说已经压得很低,这也让Q值成为比绝对分支比更“干净”的物理观测量。
龙迷三问
这篇论文到底在解决什么问题? BESIII利用百亿J/ψ和ψ(3686)样本,首次观测到Σ⁰Σ̄⁰η衰变道,测量分支比并检验12%规则,Q值为17.3%,在3.0σ内与理论一致。
这篇工作最值得看的点是什么? 首次观测到J/ψ→Σ⁰Σ̄⁰η和ψ(3686)→Σ⁰Σ̄⁰η两个衰变道,测得分支比分别为(7.5±0.3±0.8)×10⁻⁵和(1.3±0.1±0.1)×10⁻⁵,Q值为(17.3±1.5±1.7)%,与12%规则在3.0σ内一致。
这篇工作的边界或风险在哪里? 优点:基于BESIII探测器采集的世界最大J/ψ和ψ(3686)数据样本,统计量充足;首次观测到两个衰变道,填补了该领域的空白;采用事件逐例加权方法有效处理了数据与MC模拟的差异;系统误差分析全面细致。缺点:MC模型引入的系统误差较大(J/ψ道为7.1%);未发现中间共振态或阈值增强结构,物理信息有限;对12%规则的检验精度受限于统计量。
如果你还有哪些想要了解的,欢迎在评论区留言或者讨论~
龙哥点评
论文创新性分数:★★★★☆
首次观测两个粲偶素衰变道,创新性硬核;分析范式继承自BESIII成熟体系,方法本身的新意有限。
实验合理度:★★★★★
基于世界最大数据样本,事件选择、背景估计、系统误差体系完备,结论可信度高。
学术研究价值:★★★★☆
对检验12%规则、约束BB̄P三体衰变的动力学模型都有直接贡献;但未找到显著中间态,理论判别力有限。
稳定性:★★★★☆
信号显著、拟合稳定,分支比精度合理;但稀有过程的统计量决定了结果对样本依赖较强。
适应性以及泛化能力:★★★☆☆
方法针对粲偶素衰变到超子对加介子的特定过程,可推广到同类BB̄P末态分析,但无法泛化到其他领域。
硬件需求及成本:★★★☆☆
依赖BEPCII大型加速器设施,数据获取周期长、运行成本高;但离线分析阶段仅需常规计算资源。
复现难度:★★☆☆☆
数据非公开,需加入BESIII合作组方可复现;但分析流程论文描述详细,合作组内部可完整复现。
产品化成熟度:★☆☆☆☆
基础研究前沿探索,不涉及产品化。
可能的问题:
Q值中心值偏高但统计不足,无法给出“违反”或“符合”的强结论;未找到中间态使得该道对重子谱研究的增量有限。
[1] BESIII Collaboration. First measurements of the branching fractions of J/ψ and ψ(3686) → Σ⁰Σ̄⁰η. arXiv:2608.16214.
[2] BESIII Collaboration. Design and construction of the BESIII detector. Nucl. Instrum. Methods Phys. Res., Sect. A, 2010, 614: 345.
[3] 原文中关于12%规则、ρπ之谜和粲偶素衰变的详细引用文献列表,参见论文参考文献章节。
*本文仅代表个人理解及观点,不构成任何论文审核或者项目落地推荐意见,具体以相关组织评审结果为准。欢迎就论文内容交流探讨,理性发言哦~ 想了解更多原文细节的小伙伴,可以点击 "阅读原文", 查看更多原论文细节哦!
想和龙哥一起追踪BESIII、LHCb这些前沿实验的最新结果吗?想知道粒子物理如何检验标准模型、探索新物理?
欢迎加入龙哥读论文粉丝群,
扫描下方二维码或者添加龙哥助手微信号加群 :kangjinlonghelper。
一定要备注:研究方向+地点+学校/公司+昵称(如 高能物理+北京+中科院+小王) ,根据格式备注,可更快被通过且邀请进群。
『龙哥读论文』微信群目前包含:图像处理、大模型及智能体、自动驾驶及机器人、AI医疗及AI金融5个群,高能物理同好也欢迎来建新群~
*本文仅代表个人理解及观点,不构成任何论文审核或者项目落地推荐意见,具体以相关组织评审结果为准。欢迎就论文内容交流探讨,理性发言哦~ 想了解更多原文细节的小伙伴,可以点击 "阅读原文", 查看更多原论文细节哦!