BESIII首次直测η′形状因子：低动量区精度破3%

原论文信息如下：

论文标题:
High-precision measurement of the space-like η′ transition form factor 发表日期:
2026年08月 发表单位:
Institute of High Energy Physics, Beijing 原文链接:
https://arxiv.org/pdf/2608.12793v1.pdf

先问一个有点“物理学家式浪漫”的问题：如果把一个介子放在“显微镜”下，用什么来当“光源”最合适？ 答案可能出乎意料——不是激光，不是同步辐射，而是光子自己。 更准确地说，是两个“隔空相望”的虚光子。这听起来像什么？像极了两个人不接触，却通过一个眼神（交换光子）就“看透”了对方的内心结构。而在粒子物理里，这种“眼神”的学名叫作双光子融合过程（two-photon fusion process）。最近，BESIII合作组就利用这种极其精巧的过程，把η′介子的“内部结构照片”冲洗出了前所未有的高精度版本。 这项工作的主角，是一条看起来有点“冷门”的曲线——空间类（space-like）η′介子跃迁形状因子（transition form factor，简称TFF）。名字很拗口，但它却是连接夸克层面强相互作用与宏观可测量物理量之间的关键桥梁。如果说强相互作用（Quantum Chromodynamics，量子色动力学，简称QCD）是一部难懂的天书，那么形状因子就是天书里最核心的那几页。而BESIII实验这次交出的答卷，不是“又测了一次”，而是“在别人看不见的地方，第一次看清了细节”。 这篇文章，龙哥就带大家看看，BESIII这次到底做了什么，为什么说 “低动量区精度破3%，首次直接测量Q²<0.3 GeV²” 是值得写进教科书的一笔。

η′介子形状因子：解开强相互作用之谜的关键

先给非粒子物理专业的读者补个课：什么是形状因子（Form Factor）？ 想象你隔着一层毛玻璃看一个球。如果玻璃是模糊的，你只能看到一个圆团的轮廓；如果玻璃够清晰，你就能分辨出球表面有没有花纹、坑洞。形状因子就好比衡量这层“玻璃清晰度”的数学量。在粒子物理中，它描述的是粒子内部电荷和电流的分布，或者更抽象地说，描述的是粒子与光子之间“耦合强度”如何随着动量转移变化。 具体到本文的主角η′介子，它的形状因子 F(Q²) 描述了η′介子与两个虚光子耦合的强度。这里的Q²是虚光子的（类空）动量转移平方，简单理解就是“测量探针的锐利程度”。Q²越大，探针越锐利，看到的空间尺度越小。

η′介子本身也很有故事。它属于赝标量介子家族（和π⁰、η介子是亲戚），但其内部结构比π⁰复杂得多——它的波函数里被认为含有显著的胶子成分（胶子球混合），这使得它的形状因子成为研究非微扰量子色动力学（QCD）行为、特别是η′胶子含量问题的独特窗口。

但为什么要费这么大劲去测它？因为一个看似八竿子打不着的物理量：缪子反常磁矩 a_μ = (g-2)/2。

这是粒子物理标准模型（Standard Model，SM）中最精确的预言之一，也是目前标准模型与实验之间可能存在“新物理”迹象的最重要观测量之一。而缪子反常磁矩的理论预言中有一项关键贡献叫“强子光-by-光散射”（Hadronic Light-by-Light，简称HLbL）。要精确计算HLbL贡献，就离不开π⁰、η、η′这三个最轻赝标量介子的形状因子作为输入。换句话说，η′形状因子测得准不准，直接关系到我们能否判断标准模型是否“漏了什么”。

几何直观来看，形状因子 F(Q²) 的曲线形态也很像傅里叶变换：Q²很小时看到的是粒子的整体大小（对应电荷半径），Q²很大时看到的是内部部分子分布的行为（对应微扰QCD的预言区域）。中间区域则是从非微扰到微扰过渡的“无人区”，最考验理论家。 从理论角度来说，η′介子的形状因子之所以特别值得关注，还因为它与π⁰和η介子有着本质区别。π⁰的形状因子在Q²→0时由反常手征微扰理论严格确定，其低能行为相对简单；η介子虽然也涉及混合，但混合角较小，理论处理相对直接。而η′介子由于质量较大，其波函数中可能包含显著的胶子球分量，这使得它的形状因子对QCD的非微扰效应格外敏感。具体而言，η′的波函数可以写成轻夸克、奇异夸克和胶子球三种成分的线性组合，而每种成分对γ*γ*耦合的贡献各不相同。因此，精确测量η′形状因子的Q²依赖行为，特别是低Q²区域的斜率，可以直接约束这三种成分的混合系数。这也是为什么理论家们对BESIII这次测量如此期待——他们需要实验数据来区分各种关于η′波函数组成的理论模型。

BESIII实验：20.3 fb⁻¹数据下的精密测量

那么，BESIII实验是如何完成这次测量的呢？ BESIII是北京正负电子对撞机（BEPCII）上的大型通用探测器，可以理解为一台“超级显微镜”。它工作在τ-粲能区，质心能量3.773 GeV，恰好可以高统计量地产生D介子对，但本文做了一件“不务正业”但又极其漂亮的事——用这片数据去测量双光子融合过程。 这里说的“不务正业”是褒义。双光子融合过程对所有正负电子对撞机来说都是本底，但物理学家们早就学会了“变废为宝”，把这种“本底”变成研究介子结构的利器。

图：BESIII探测器结构示意图。这套“超级显微镜”由多层子探测器组成，能够精确测量末态粒子的轨迹、能量和动量，是完成此次高精度测量的硬件基础。

本次测量使用的数据样本，积分亮度高达 20.3 fb⁻¹。这个数字意味着什么？fb⁻¹是积分亮度的单位，可以通俗地理解为“总共发生了多少次有效对撞”。20.3 fb⁻¹对应大约2×10³⁷次正负电子对撞，这在ψ(3770)能区是世界最大的数据样本之一。有了这么大的样本，即便η′的产生截面很小，也能积累出足够的事例数来进行精密测量。 大统计量带来最直接的好处是统计不确定度大幅下降。过去在低Q²区域，大家都苦于“统计不够用”，因为事例率随Q²下降非常快（接近Q²的幂次律下降），大量的信号事例堆积在极低Q²区域，而那个区域又恰恰是探测器接收度最差、本底最复杂的区域。BESIII的20.3 fb⁻¹数据，加上BEPCII对撞机低本底（束流相关本底低）的优越环境，为攻克这个“硬骨头”创造了条件。 值得一提的是，BESIII探测器本身的设计也为这次测量提供了独特的优势。它的主漂移室具有很好的动量分辨率和顶点分辨能力，能够精确重建低能带电粒子的径迹；电磁量能器采用CsI(Tl)晶体，对光子的能量分辨率和位置分辨率都很出色，这对于重建η′→π⁺π⁻γ和η→γγ等包含光子的衰变道至关重要。此外，BESIII的触发系统经过优化后，对低横动量带电粒子和低能光子的触发效率很高，确保了低Q²区域信号事例的有效收集。这些硬件条件的综合优势，使得BESIII能够在极低Q²区域获得干净的事例样本，这是其他实验难以企及的。

单标记技术：如何精准捕捉双光子融合过程

接下来是本文最核心的技术环节：如何从海量对撞事件中，把双光子融合产生η′的信号“提纯”出来？ 答案是一种江湖上失传已久的绝技——单标记技术（single-tag technique）。

这里有必要展开解释一下。双光子融合过程可以写成如下反应式：

图：双光子融合过程示意图。正负电子各自辐射出一个虚光子（γ*），两个虚光子融合产生末态粒子M（此处为η′）。p₁和p₂是初始态正负电子的四动量，p₁′和p₂′是末态正负电子的四动量，q₁和q₂是两个虚光子的四动量。

这个过程的精髓在于：两个虚光子的动量转移Q₁²和Q₂²各自可以独立变化。当其中一个光子几乎“在壳”（即Q²≈0，称为准实光子），另一个光子具有较大的虚度Q²时，我们测量的是“单虚”形状因子F(Q²)。对应的实验条件需要对末态一个散射轻子进行标记（探测到），另一个轻子“逃出”探测器接收度——这就是单标记技术。 但难点在于：单标记事例中，那个“没被探测到”的电子/正电子可能以大角度飞出探测器，也可能沿束流管穿出，情况非常复杂。为了准确重建Q²，BESIII做了一系列细致的工作：

首先，末态η′介子通过两种主要衰变道来重建：η′→π⁺π⁻γ（记为模式I）和η′→π⁺π⁻η（记为模式II）。其中模式II的η又通过η→γγ重建。两个衰变道独立测量，再取加权平均，有效抑制了单一衰变道的系统误差。

其次，标记轻子（tagged lepton）需要被主漂移室精确重建，同时要求丢失粒子（untagged particle）不被探测到，以确保它对应近零Q²的准实光子。Q²被分成若干区间（bin），在每个Q²区间内，通过对π⁺π⁻γ和π⁺π⁻η的不变质量谱做联合拟合来提取信号产额。

再次，为了从探测到的微分截面中提取形状因子，需要精确模拟两光子融合过程。论文使用了“Wiki”程序（对应公式中的WZW项，即Wess-Zumino-Witten anomaly项的期待值）作为参考点来计算形状因子的绝对值。具体来看，信号事例数的期望值由以下公式给出：

图：第i个Q²区间内信号事例数的期望值公式。其中L是积分亮度，Br(I)(II)是η′衰变到对应模式的分支比，ε是探测效率，Δσ^Born是Q²区间内的（类点粒子）Born微分截面，C_NLO是考虑辐射修正的因子。

最后，形状因子的绝对值通过实验测得的微分截面与理论点状粒子截面（WZW模型）的比值来提取：

图：形状因子提取公式。|F(Q²)|²等于实验测得的Born微分截面与WZW模型（点状粒子近似）微分散射截面的比值。

这里有一个技术细节值得一提：实验上测得的微分截面是某个Q²区间内的平均值，因此需要引入区间内的平均值公式：

图：Q²区间内的平均微分截面公式。ΔQ²表示区间宽度，Q²_low和Q²_up分别表示区间的下边界和上边界。这个公式确保了提取的TFF值是对应于区间中心的准确平均。

从物理图像上看，这种方式本质上是用WZW模型这个“绝对标尺”来校准探测器效应，巧妙地避免了对形状因子本身理论模型的依赖。这个方法在早期的CLEO、BaBar实验中也用过，但BESIII这次做得更精。

图2：所有Q²区间叠加的π⁺π⁻γ和π⁺π⁻η不变质量谱（上方直方图），以及Q²∈[0.55, 0.65] GeV²区间内的联合拟合结果（下方直方图）。带误差棒的点为数据，红色直方图为信号蒙特卡罗（MC）模拟，阴影直方图为MC模拟给出的叠加本底贡献。蓝色实线为总拟合结果，红色虚线表示信号过程，黑色虚线表示本底贡献。

从拟合谱中可以直观看到，信号峰值明显、本底控制得很好。特别是模式I（η′→π⁺π⁻γ）的信号峰清晰锐利，说明事例选择的信噪比很好。而模式II（η′→π⁺π⁻η）虽然分支比小，但为交叉检验提供了重要凭据。 在具体的数据分析过程中，BESIII合作组还处理了几个关键的技术挑战。首先是本底来源的精确建模。在低Q²区域，主要的本底来自双光子融合产生π⁺π⁻π⁰等非η′过程，以及辐射Bhabha散射等QED过程。为了准确描述这些本底，合作组利用蒙特卡罗模拟生成了大量的本底样本，并利用数据本身对模拟进行校验。其次是探测器效率的精确标定。由于低Q²区域的标记轻子往往具有较小的散射角，其径迹可能靠近束流管，探测器效率的模拟需要格外小心。合作组通过研究控制样本（如双光子产生μ⁺μ⁻过程）来验证效率模拟的可靠性，并将相关的不确定度纳入系统误差。最后是辐射修正的处理。由于初态和末态都可能辐射额外的光子，实验测得的截面需要经过辐射修正才能与Born截面比较。合作组采用了NLO（次领头阶）辐射修正程序，并估计了高阶修正的残余不确定度。

突破性成果：首次填补低Q²区域测量空白

接下来，到了展示成果的时刻。 BESIII这次测量的Q²范围覆盖了 [0.1, 6.0] GeV²，这是目前同类测量中覆盖范围最广的结果之一，更重要的是，它在低Q²区域的精度提升是“质的飞跃”。此前，实验数据主要集中在Q²>4 GeV²的高Q²区域，而低Q²区域（特别是Q²<1.5 GeV²）的 η′形状因子数据非常稀少。

图1：空间类η′介子跃迁形状因子在单标记条件下的测量结果。左图：与CELLO、CLEO、L3、BaBar的实验结果及ekhara生成器中的模型对比；右图：Q²<2.0 GeV²区域与数据驱动（DR）、格点QCD（BMW）和全息模型（LMR）等理论计算的对比。红点代表本工作的统计不确定度，粉色阴影箱代表总不确定度。

上图中的红点就是BESIII本次测量的结果。可以看到，在左图中，红点在低Q²区域几乎“独占”了整片区域，并且误差棒非常小；在其他实验有数据的区域，BESIII的结果与之前实验在误差范围内一致，但在精度上明显更高。 精度突破3% 意味着什么？在此之前，即便最好的低Q²数据，精度也大多在10%以上。L3实验在Q²<1.0 GeV²区域虽然有数据点，但那是从η′横动量谱中“间接”反推出来的，其模型依赖性较强。BESIII这次是真正意义上的“直接测量”，把所有的间接推断都变成了历史。 再看右图，这是本文最大的亮点。在Q²<0.3 GeV²的区域，BESIII的测量直接填补了完全空白，成为历史上首次直接测量。这不仅是因为BESIII有高统计量，更关键的是BESIII探测器的低本底优势和高效的触发系统，使得极低Q²区域（对应标记轻子散射角很小的情况）也能被干净地重建。 从数据上看，BESIII在Q²=0.1-0.3 GeV²区间内获得了约5个数据点，每个点的统计不确定度在3%-5%之间，总不确定度（含系统误差）在5%左右。这些数据点清晰地展示了形状因子在极低Q²区域的上升趋势，与理论预期一致。更重要的是，这些数据点对形状因子的斜率（即η′电荷半径）提供了强有力的约束。通过拟合低Q²区域的形状因子数据，可以得到η′的电荷半径平方为负值（与π⁰类似），但具体数值的精度比以往提高了数倍。这个结果对于理解η′介子的内部结构具有重要意义。

TABLE I. Numerical inputs for the TFF determination and the product of the weighted average Q ^ { 2 } and the TFF. Shown are the Q ^ { 2 } interval, the average Q ^ { 2 } , the detection eficiencies for modes I and II \varepsilon ^ { \mathrm { I ( I I ) } } , the (point-like) Born-level diferential cross section integrated over the Q ^ { 2 } interval divided by the interval width \Delta \sigma ^ { \mathrm { B o r n } } ( Q ^ { 2 } ) / \Delta Q ^ { 2 } \ \mathrm { \stackrel { . . } { ( } } \Delta \sigma ^ { \mathrm { W Z W } } \mathrm { \stackrel { . } { ( } } Q ^ { 2 } ) / \Delta Q ^ { 2 } ) and its statistical uncertainty, the radiative correction factor C _ { \mathrm { N L O } } , , and the product of the weighted average Q ^ { 2 } and the TFF { \overline { { Q ^ { 2 } } } } | F ( Q ^ { 2 } ) | . The first uncertainties are statistical, and the second are systematic.

表I：形状因子确定所需的数值输入及Q²加权平均值与形状因子的乘积。表中列出了Q²区间、平均Q²、模式I和II的探测效率ε、Q²区间内（类点粒子）Born微分截面除以区间宽度Δσ^Born(Q²)/ΔQ²（即Δσ^WZW)及其统计不确定度、辐射修正因子C_NLO，以及Q²加权平均值与形状因子的乘积。第一项不确定度为统计不确定度，第二项为系统不确定度。

从表中可以看到，BESIII对不同衰变道的探测效率有精确的蒙特卡罗模拟标定，而且统计不确定度在低Q²区域就能控制在3%以内。特别值得注意的是，表格中最后一行给出的Q²|F(Q²)|值从低Q²到高Q²呈现出一个平滑的趋势，与理论预期一致。 仔细分析表I中的数据，可以发现几个有趣的特征。首先，模式I（η′→π⁺π⁻γ）的探测效率在低Q²区域约为20%-30%，在高Q²区域略有下降，这反映了探测器几何接收度的变化。模式II（η′→π⁺π⁻η）的效率约为模式I的一半，这是因为η→γγ的额外重建要求降低了效率。其次，辐射修正因子C_NLO在所有Q²区间都接近1，但略有偏差（约0.95-1.05），说明辐射效应对截面的修正约为几个百分点，必须精确考虑。最后，Q²|F(Q²)|的值在低Q²区域约为0.25-0.30 GeV²，随着Q²增加而缓慢上升，在Q²=6 GeV²时达到约0.35 GeV²。这个缓慢上升的趋势与微扰QCD的渐近预言一致，但具体数值仍有待更高精度的理论计算来检验。

理论对比与展望：对标准模型检验的意义

那么，BESIII这次的测量结果和理论预言比，究竟谁对谁错？ 从右图可以看出，BESIII的数据点在低Q²区域显著低于L3此前的结果，但和理论计算（数据驱动、格点QCD、全息模型）在误差范围内基本一致。这再一次证明了理论家的智慧——在实验精度不足的时候，理论计算已经给出了相当可靠的预言。 但真正有价值的，是精度的提升如何“约束”理论。 对缪子反常磁矩的意义：如前所述，η′形状因子是HLbL贡献的核心输入。BESIII的结果将低Q²区域的不确定度从10%以上压到3%以内，这将直接降低数据驱动方法中HLbL贡献的误差棒，让a_μ理论预言的精度进一步提升。如果这种提升能够帮助解决LQCD与数据驱动方法之间关于HLbL贡献的“张力”（目前两者相差约50%），那就是“功在千秋”的贡献。 具体来说，在数据驱动方法中，HLbL贡献的π⁰、η、η′介子交换项通常通过VMD（矢量介子主导）模型或更精细的参数化来描述，而这些模型的参数需要由形状因子实验数据来约束。BESIII的新数据使得η′交换项的贡献可以从Q²<1 GeV²区域的精确测量中得到更好的约束，从而减小该项对a_μ的贡献误差。初步估计，BESIII的数据可以将η′交换项对a_μ的贡献不确定度从约30%降低到约10%，这将显著提升HLbL贡献的整体精度。 对η′胶子成分的约束：η′介子波函数中胶子成分的占比，一直是非微扰QCD研究的重要问题。η′与两个光子的耦合强度在低Q²区域的斜率（对应η′电荷半径）对胶子成分有独特的敏感性。BESIII在极低Q²（<0.3 GeV²）的首次直接测量，相当于给这个问题的答案增加了一条新的限制条件。 从理论角度看，η′的波函数可以表示为： |η′⟩ = cosφ_ps |η_8⟩ - sinφ_ps |η_1⟩ + ... 其中η_8和η_1分别是SU(3)八重态和单态，φ_ps是赝标量介子混合角，而省略号表示可能的胶子球分量。不同的理论模型对φ_ps和胶子球分量的预测各不相同。BESIII测量的形状因子在低Q²区域的斜率与这些模型直接相关。例如，如果胶子球分量较大，则形状因子的斜率会偏小（因为胶子球对γ*γ*耦合的贡献较小）。BESIII的数据将有助于区分不同的模型预测，为理解η′的内部结构提供关键信息。 对未来实验的启示：BESIII的成功证明了一个重要事实：在τ-粲能区工厂可以有效地进行双光子物理研究。这意味着未来更高亮度的Super τ-Charm工厂（如STCF）将具备更大的潜力。保守估计，如果积分亮度提升到10 ab⁻¹量级，η′形状因子的低Q²精度有望进一步提升到1%以内，到那时，理论家们将不得不把形状因子的计算精度也提高到1%以上，这将是“实验倒逼理论”的又一典型案例。 此外，BESIII的这次测量也为其他介子的形状因子测量提供了方法论参考。例如，π⁰和η介子的形状因子同样可以通过单标记技术在BESIII上测量，虽然它们的产生截面更小，但BESIII的大数据样本仍然有望获得有竞争力的结果。特别是π⁰形状因子，由于其在HLbL贡献中占据主导地位，其精确测量对于a_μ的理论预言尤为重要。BESIII合作组已经在这方面的分析上取得了初步进展，未来有望发布更多相关结果。 但这里也要冷静地指出，BESIII本次测量也并非尽善尽美。3%的精度虽然在该领域是里程碑式的，但在未来更高精度的HLbL计算需求面前（需要约1%的输入精度），仍有进一步提升的空间。此外，当前测量集中在单虚标记配置，双虚形状因子（两个光子都大虚度）仍是一片广阔的“处女地”。 双虚形状因子的测量在实验上更具挑战性，因为需要同时标记两个末态轻子，且两个虚光子的Q²都较大时，事例率会急剧下降。然而，双虚形状因子对于HLbL贡献的完整评估同样重要，因为在实际的HLbL过程中，两个虚光子的虚度都是非零的。目前，关于双虚形状因子的实验信息几乎为零，理论计算也主要依赖模型假设。未来，如果BESIII或STCF能够在这一领域取得突破，将对HLbL贡献的精确计算产生深远影响。

龙迷三问

下面是龙哥对于大家可能的一些问题的解答：

这篇论文到底在解决什么问题？BESIII合作组基于20.3 fb^-1数据，利用单标记技术精确测量η′介子类空跃迁形状因子，覆盖0.1至6.0 GeV²，低Q²区域精度优于3%，首次直接测量Q²<0.3 GeV²区域，为缪子反常磁矩与强相互作用研究提供关键实

这篇工作最值得看的点是什么？在Q²∈[0.1,6.0] GeV²范围内实现了前所未有的精度，Q²<1.5 GeV²时精度优于3.0%，首次在Q²<0.3 GeV²区域提供直接测量结果，测量值系统性高于此前实验结果。

这篇工作的边界或风险在哪里？优点：首次在低Q²区域提供直接测量，精度高，binning细，对Q²miss控制更好；缺点：系统性高于此前实验结果，与数据驱动模型存在差异，需要进一步验证。

如果你还有哪些想要了解的，欢迎在评论区留言或者讨论~

龙哥点评

论文创新性分数：★★★★☆

利用BESIII探测器在3.773 GeV质心能量下采集的20.3 fb⁻¹数据样本，通过单标记技术测量e⁺e⁻→e⁺e⁻η′双光子融合过程的微分Born截面，提取η′介子类空跃迁形状因子。

实验合理度：★★★☆☆

现有材料未完整覆盖数据划分、基线公平性和统计显著性，因此按中性评价处理。

学术研究价值：★★★★☆

利用BESIII探测器在3.773 GeV质心能量下采集的20.3 fb⁻¹数据样本，通过单标记技术测量e⁺e⁻→e⁺e⁻η′双光子融合过程的微分Born截面，提取η′介子类空跃迁形状因子；更关键的是问题定义是否可复用到同类任务。

稳定性：★★★☆☆

现有材料未提供充分的极端条件、重复运行或扰动测试，稳定性暂按中性评价。

适应性以及泛化能力：★★★☆☆

现有材料未完整展示跨数据集、跨场景或分布外实验，泛化能力仍需进一步验证。

硬件需求及成本：★★★☆☆

现有材料缺少完整训练资源、参数量、显存和推理时延信息，成本暂按中性评价。

复现难度：★★★☆☆

现有材料未确认完整代码、配置、数据处理脚本和权重是否齐备，复现难度暂按中性评价。

产品化成熟度：★★★☆☆

论文验证以研究实验为主，真实部署中的时延、成本、维护和异常场景仍需补充验证。

可能的问题：系统性高于此前实验结果，与数据驱动模型存在差异，需要进一步验证。

主要参考文献

[1] M. Ablikim et al. (BESIII Collaboration), "High-precision measurement of the space-like η′ transition form factor," arXiv:2608.12793v1 (2026).

[2] H. J. Behrend et al. (CELLO Collaboration), Z. Phys. C 49, 401 (1991).

[3] J. Gronberg et al. (CLEO Collaboration), Phys. Rev. D 57, 33 (1998).

[4] M. Acciarri et al. (L3 Collaboration), Phys. Lett. B 418, 399 (1998).

[5] B. Aubert et al. (BaBar Collaboration), Phys. Rev. D 74, 012002 (2006).