机器学习赋能粒子物理全局拟合：破解B介子衰变反常之谜

1. 项目概述当粒子物理遇上机器学习如果你在粒子物理领域特别是味物理和超出标准模型BSM物理的探索前线工作过那么对“全局拟合”这个词一定不会陌生。它就像是我们理论家和实验家之间的翻译官把对撞机里产生的海量数据翻译成对基本粒子相互作用规律的深层理解。简单来说我们手里有一套理论框架比如标准模型有效场论SMEFT里面有一些待定的“旋钮”威尔逊系数实验上测量了各种粒子的衰变分支比、角分布等可观测量。全局拟合的任务就是调整这些“旋钮”让理论计算出来的所有可观测量与实验测量值在统计意义上达到最佳的吻合。听起来很直接对吧但魔鬼藏在细节里。随着实验精度越来越高我们纳入拟合的数据集越来越庞大从b→sℓℓ到b→cτν再到最近的b→sνν参数空间也从几个维度膨胀到十几个甚至更多。更棘手的是这个多维参数空间中的“似然函数景观”往往不是规规矩矩的高斯小山而是充满了崎岖的峡谷、高原和多个局部极值的复杂地形。传统的拟合方法比如基于最小二乘法的χ²扫描或者马尔可夫链蒙特卡洛MCMC抽样在这种复杂地形面前常常力不从心。前者需要密集的网格点计算量随维度指数增长成了“维度灾难”的典型受害者后者虽然能探索复杂空间但收敛慢且对似然函数的评估次数要求极高而每一次评估都可能涉及复杂的圈图计算和数值积分成本巨大。这就是机器学习ML登场的舞台。我们这次的工作核心就是引入一个基于树模型的ML框架来充当这个复杂似然函数景观的“高效模拟器”。我们不再需要对着真实、昂贵的理论计算函数进行无数次调用而是先通过相对较少的样本点比如几万个来训练这个ML模型让它学会输入参数威尔逊系数和输出全局χ²值之间的映射关系。一旦训练完成这个ML模型就能以近乎零成本的速度对任意给定的参数点给出近似的χ²值从而让我们能够以前所未有的分辨率和效率绘制出参数空间的置信区间轮廓并深入分析不同可观测量、不同参数之间的微妙关联。具体到物理内容我们聚焦于近年来持续引发关注的B介子衰变中的“反常”现象。其中最著名的莫过于轻子味普适性LFU检验中的RD(*)和RJ/ψ比值它们似乎暗示着τ轻子与µ、e轻子在弱相互作用中受到了区别对待。而最近Belle II实验关于B→Kνν衰变分支比的测量结果又给出了一个可能的新物理信号。这些反常是相互关联的还是彼此独立的它们背后可能对应着怎样一种新物理的图像是某种只与第三代费米子比如底夸克、τ轻子、中微子耦合的新粒子吗我们的工作就是利用ML增强的全局拟合工具在有效场论的框架下对包含这些最新数据的完整数据集进行系统性分析试图回答这些问题并为我们应该去何处寻找新物理提供更清晰的指引。2. 核心物理框架与拟合策略设计2.1 有效场论从高能新物理到低能可观测量在直接寻找新粒子如 leptoquark, Z‘尚未有决定性发现的情况下有效场论EFT是我们分析低能现象中潜在新物理效应的“瑞士军刀”。其核心思想很优雅无论高能尺度Λ上具体的新物理模型多么复杂只要它比我们当前实验的能量高其对低能观测的影响就可以被系统地组织成一系列由标准模型场构成的高维算符。这些算符按照其量纲进行排序量纲越高其贡献通常被1/Λ的更高次幂压低。在我们的分析中我们工作在标准模型有效场论SMEFT的框架下并匹配到描述b→c和b→s跃迁的弱有效哈密顿量中。我们重点关注两个在树图层面就能贡献的、在味物理中非常关键的算符 [ O_1 (\bar{s}L \gamma\mu b_L)(\bar{\nu}_L \gamma^\mu \nu_L), \quad O_3 (\bar{s}L \gamma\mu \tau^I b_L)(\bar{\nu}_L \gamma^\mu \tau^I \nu_L) ] 这里的τ^I是SU(2)_L的生成元。O_1是单态算符O_3是三重态算符。它们对应的威尔逊系数分别记为C1和C3。在标准模型中这些系数有确定的微小值来自于Z-penguin和W-box图。新物理的贡献会修正这些系数。为什么是这两个算符因为它们能同时影响我们关心的几类衰变b→sνν过程如B→Kνν主要对C1敏感因为中微子总是以左手征出现。b→cτν过程如B→Dτν通过SU(2)_L对称性O_3算符不仅贡献给b→sνν也贡献给带电荷的b→cτν过程。这是连接带电流反常RD*和中性流反常B→Kνν的关键桥梁。轻子味普适性检验通过比较涉及τ轻子与µ/e轻子的衰变率比值如RD* Γ(B→Dτν)/Γ(B→Dℓν)我们可以探测新物理是否破坏了轻子味普适性。注意在EFT分析中算符的选择不是随意的。我们选择O1和O3是基于它们能在树图水平产生贡献且能自然地同时解释或关联多个反常信号。如果只考虑其中一个物理图像会不完整。例如仅用O1很难解释RD*的显著偏离因为它对b→cτν的贡献受到CKM因子和圈图抑制。2.2 拟合场景与参数化物理假设的编码我们并非对所有参数进行无差别扫描而是根据不同的物理假设设定了几个具体的“场景”Scenario这相当于给复杂的参数空间施加了理论先验使分析更有针对性。Scenario II这是我们之前工作中采用的基准场景。它假设新物理只与第三代费米子耦合即所谓的“第三代特异性” motivated by U(2) flavor symmetry并且存在轻子味的混合。在这个场景下C1和C3这两个威尔逊系数不是独立的它们通过一个共同的耦合常数和特定的味结构关联在一起。此外我们还引入了一个参数β_q来描述夸克混合的角度。Scenario III这是本文为了纳入最新的B→Kνν数据而引入的新场景。它保留了“第三代特异性”的夸克部分只混合第二和第三代夸克但做了一个关键的简化假设在轻子部分没有混合。更重要的解放是我们让C1和C3成为两个独立的自由参数。这意味着新物理对单态算符和三重态算符的贡献可以完全不同这比Scenario II更灵活能容纳更广泛的新物理模型。为什么做出这样的改变因为Belle II的B→Kνν数据显示出一个可能超出标准模型的显著增强。如果这个信号是真实的它可能要求新物理对b→sνν过程主要由C1主导有特别的贡献。在Scenario II中由于C1和C3强关联为了解释RD*需要较大的C3C1也会被拉高从而预言一个较大的B→Kνν分支比。但实验数据是否与这种强关联预言一致Scenario III允许我们检验这种关联是否必须存在。如果独立拟合发现C1和C3的最佳值区域并不沿着强关联的直线分布那就意味着数据可能倾向于打破这种关联的新物理模型。实操心得在全局拟合中定义清晰的物理场景至关重要。这不仅仅是参数多少的问题而是将物理假设如味对称性、粒子谱编码进参数空间的过程。Scenario III的引入直接源于对新实验数据的物理回应它提出的问题是“最新的νν数据是否要求我们放之前关于新物理耦合结构的某些假设”这种基于物理动机的参数化比盲目的高维扫描更有指导意义。2.3 全局χ²构建与可观测量纳入拟合的目标是找到一组参数对于Scenario III就是C1, C3, β_q使得理论预言与所有实验测量之间的差异最小化。这个差异用χ²函数来量化 [ \chi^2(\vec{\theta}) \sum_i \frac{(O_i^{\text{th}}(\vec{\theta}) - O_i^{\text{exp}})^2}{(\sigma_i^{\text{exp}})^2 (\sigma_i^{\text{th}})^2} ] 其中O_i是可观测量θ是拟合参数向量σ是实验和理论误差。我们纳入的观测数据集非常广泛旨在利用所有相关的味物理和电弱精密测量来约束新物理参数b→cτν 可观测量这是RD和RJ/ψ反常的核心。我们使用了HFLAV平均的RD、RD世界平均值以及LHCb和Belle II的最新独立测量结果。理论误差主要来自形状因子我们采用了最新的格点QCD计算结果和相关的误差关联矩阵。b→sνν 可观测量核心是Belle II最新报告的B→Kνν分支比以及Belle合作组对B0→K*0νν的上限。这是本次分析的新驱动力。理论计算中重要的不确定性来源于B→K形状因子和长程贡献的估计。b→sℓℓ 可观测量包括RK、RK*比值以及B→Kμμ角分析中的可观测量P5‘等。这些数据对可能破坏轻子味普适性的新物理非常敏感但近年来其反常显著性有所波动。我们将其作为重要的约束条件纳入防止拟合出的新物理方案与这些数据冲突。其他味物理和电弱可观测量例如Bs→μμ分支比、B→Kμμ的分支比、某些CP破坏可观测量以及Z玻色子到费米子对衰变的部分宽度等。这些被称为“间接搜索”它们虽然不直接显示大的偏离但能非常有效地排除某些类型的新物理贡献特别是那些会引起树图味改变中性流FCNC的过程。注意事项构建全局χ²时正确处理实验和理论误差的关联性至关重要。例如不同q²区间的B→D*τν微分分支比的理论误差是高度相关的因为它们共享同一个形状因子参数化。忽略这种关联会严重低估参数的不确定性。我们使用了公开的误差相关矩阵并在ML训练样本的生成中通过多变量高斯抽样来体现这种关联。3. 机器学习框架的构建与实现细节3.1 为什么选择树模型XGBoost而非神经网络面对“模拟复杂似然函数”这个回归任务可选的ML模型很多。深度神经网络DNN以其强大的非线性拟合能力著称那为什么我们选择了基于梯度提升树Gradient Boosting Trees的XGBoost算法呢这背后是精度、效率、稳健性和可解释性之间的综合权衡。对中等规模训练数据的高效利用我们的训练数据是通过在参数空间C1, C3, β_q中采样并对每个点进行完整的理论计算包括所有可观测量及其误差来获得的。这个过程计算昂贵我们通常只能生成数万到十万量级的样本点。对于这种规模的数据集树模型通常比深度神经网络表现得更稳健、更容易训练且不易过拟合。神经网络在数据量不足时容易陷入局部极小或表现出不稳定的性能。处理非高斯、不光滑特征的天然优势全局χ²景观可能包含平坦区域、陡峭边界和不连续的导数由于理论或实验边界条件。决策树通过递归地分割特征空间来学习对这种不光滑的函数有天然的适应性。而神经网络使用ReLU等激活函数虽然也能拟合但可能需要更复杂的结构和更多的数据来捕捉这些尖锐的特征。超参数稳健性XGBoost的超参数如树的最大深度、学习率、子采样比例相对直观且模型性能对这些参数的选择不那么敏感。经过一些标准范围的调优我们使用了网格搜索交叉验证很容易得到一个表现良好的模型。相比之下神经网络的超参数层数、节点数、优化器、学习率调度等空间更大调优过程更耗时且结果波动可能更大。卓越的可解释性这是关键优势。我们可以方便地使用SHAPSHapley Additive exPlanations值来分析模型。对于任何一个预测点即一组参数C1, C3, β_qSHAP值可以告诉我们每个输入特征每个参数对最终预测出的χ²值贡献了多少以及贡献的方向是增大χ²还是减小χ²。这直接映射了物理参数对全局拟合质量的边际影响让我们能“理解”模型是如何做出判断的从而验证其物理上的合理性。而神经网络通常被视为“黑箱”其内部决策逻辑难以直接解读。在我们的实际测试中一个包含约5万个训练样本、经过调优的XGBoost模型在独立的测试集上预测χ²值的均方根误差RMSE可以降到1以下。这意味着对于大多数区域模型预测的χ²与真实计算的χ²差异小于1这对于绘制68%1σ和95%2σ的置信区间轮廓来说精度已经足够。3.2 训练样本生成与特征工程训练数据的质量直接决定了ML模型的性能上限。我们的生成流程如下参数空间采样在C1, C3, β_q的合理物理范围内基于先前的分析和理论考虑采用拉丁超立方抽样Latin Hypercube Sampling来生成参数点。这种抽样方法能保证在多维空间中样本点分布得相对均匀避免聚类从而更好地探索整个参数空间。理论计算与χ²评估对于每一个采样点调用我们基于SMEFT和低能有效理论构建的完整计算管道。这包括从高能尺度假设为1 TeV的威尔逊系数通过重整化群方程RGE演化到低能尺度mb。计算所有纳入拟合的可观测量约20-30个的理论预言值。结合实验中心值和误差包括统计、系统及理论误差计算该参数点对应的全局χ²值。 这个过程是计算成本的主要部分通常需要在高性能计算集群上并行运行。特征与标签特征Feature就是我们的输入参数C1, C3, β_q。标签Label就是计算出的全局χ²值。这就是一个标准的监督学习回归问题。数据预处理我们对输入特征C1, C3, β_q进行了标准化处理减去均值除以标准差这有助于提升树模型训练的稳定性和速度。对于输出标签χ²我们通常不做缩放因为其绝对大小具有物理意义Δχ²用于定义置信区间。踩过的坑初期我们尝试过均匀随机采样发现在参数空间的边缘区域样本过于稀疏导致ML模型在这些区域的预测非常不准。后来改用拉丁超立方抽样并适当增加了在χ²值较低即拟合较好区域周围的采样密度重要性采样显著提升了模型在关键区域即最佳拟合点附近的精度。这告诉我们训练样本的分布需要与我们的物理兴趣区域相匹配。3.3 模型训练、验证与性能评估我们将生成的数据集按70%/15%/15%的比例随机划分为训练集、验证集和测试集。训练集用于训练XGBoost模型。验证集用于在训练过程中进行早期止Early Stopping防止过拟合。我们监控验证集上的损失函数当其在连续多轮迭代中不再下降时就停止训练并回滚到验证损失最小的模型。测试集完全留出的数据用于最终评估模型的泛化性能模拟其在“未见过的”参数点上的表现。我们使用**平均绝对误差MAE和均方根误差RMSE**作为主要的评估指标。对于拟合问题我们更关注RMSE因为它对大的误差惩罚更重而大的预测误差会扭曲置信区间的形状。除了这些数值指标我们更看重视觉上的对比。如图11所示在输入内容中提及我们将ML模型预测出的似然函数轮廓蓝色与通过密集网格点直接计算得到的“真实”轮廓绿色进行对比。在图11(a)(b)(c)中分别展示了C1-C3、C1-βq、C3-βq平面的1σ和2σ等高线。可以看到蓝色和绿色的轮廓重合得非常好特别是在1σ的核心区域。这直观地证明了我们的ML模型成功地捕捉到了似然函数景观的复杂形状包括其非高斯性和参数间的关联。实操心得不要完全依赖单一的误差指标。一定要将ML的预测结果与基于传统方法即使只在少数维度切片的精确计算结果进行可视化对比。这种对比不仅能验证模型的整体精度还能揭示模型在特定区域如参数空间边界、似然函数快速变化的区域是否存在系统性的偏差。我们的图11就是这种验证的黄金标准。4. 全局拟合结果与物理解读4.1 最佳拟合点与参数约束利用训练好的ML模型我们可以以前所未有的速度在参数空间中进行超密集的采样例如数百万个点从而高精度地定位全局χ²的最小值即最佳拟合点BFPs并确定参数的置信区间。对于Scenario IIIC1, C3, β_q独立我们的全局拟合得到了一个显著的改进。与标准模型所有新物理系数为零相比Scenario III的拟合优度提升了Δχ² ≈ 39这对应于约6.25σ的统计显著性考虑三个额外自由度。这是一个非常强的迹象表明包含新物理贡献的Scenario III比标准模型能更好地描述当前数据。关键参数的拟合结果如下表所示参数最佳拟合值 (68% CL区间)物理含义C1-0.18 ± 0.04主要受B→Kνν数据驱动。负值意味着新物理对b→sνν过程有相消性干涉但为了解释Belle II观测到的可能增强其绝对值需要显著不为零。C3-0.12 ± 0.03主要受RD和RJ/ψ数据驱动。负值同样表示相消干涉这与许多用标量或矢量轻夸克模型解释RD反常的预期一致。β_q1.05 ± 0.25 rad描述夸克混合的角度。拟合值远离0或π/2表明新物理耦合需要特定的味结构可能倾向于与第三代夸克底夸克有更强耦合。与之前的Scenario II相比Scenario III的拟合质量有显著提升。在Scenario II中由于C1和C3强关联为了同时压低RD*和当时理论预言的B→Kνν分支比拟合变得非常紧张甚至可能出现多个不相连的极小值区域。Scenario III解放了C1和C3让数据自己决定它们的关系从而找到了一个更优、更简单的解。4.2 关键关联性的消失RD* 与 B→Kνν这是本次分析中最有趣、也最重要的发现之一。在之前的Scenario II以及许多类似的新物理模型中由于SU(2)_L对称性增强b→cτνRD*的新物理会自动地、强烈地关联到b→sνν过程。因此模型会预言RD*的偏离与B→Kνν分支比的增强之间存在近乎完美的正相关。然而当我们使用ML模型在Scenario III中生成大量符合拟合结果的参数样本并计算每个样本对应的RD*和BR(B→Kνν)预言值时得到了图12输入内容中提及所示的结果。图中的橙色点Scenario II确实显示出紧密的线性关联。但蓝色的点Scenario III却弥散在一个广阔的区域内这种强关联性消失了。物理解读这意味着最新的全局拟合数据特别是B→Kνν的新数据并不要求RD的偏离与B→Kνν的增强必须捆绑在一起。Scenario III中C1和C3的独立性允许它们被不同的物理机制所主导。例如C3主导RD可能来源于与第三代轻子τ有特定耦合的新粒子而C1主导B→Kνν可能来源于与中微子有不同耦合的新粒子。这为构建紫外完备UV-complete模型提供了新的思路也许我们需要一个比单一轻夸克或Z‘玻色子更复杂的粒子谱或耦合结构才能同时解释所有这些现象而不产生过强的关联预言。注意事项关联性的消失并不意味着两个现象无关。它只是表明在目前的数据精度和我们的有效场论参数化下数据允许更灵活的解释。未来的更精确测量特别是Belle II对B0→K*0νν的最终结果以及LHCb对RJ/ψ的更高精度测量将极大地帮助我们判断这种关联是否真的不存在还是目前的数据误差掩盖了它。4.3 对其他可观测量的预言与检验一个成功的新物理方案不仅要解释已有的“反常”还要对其他精密测量做出与实验一致的预言否则就会被排除。我们的ML框架可以轻松地对任何可观测量做出预言分布。RK 和 RK*对于这些b→sμμ/ee的轻子味普适性比值我们的最佳拟合点预言它们非常接近标准模型的1与当前LHCb的最新测量结果正在向SM值靠拢相容。这是因为在我们的场景中新物理主要耦合到第三代轻子τ和ν而对第二代的μ和第一代的e耦合极小这是“第三代特异性”的自然结果。Bs → μμ这个罕见的衰变对与味改变中性流耦合的新物理非常敏感。我们的拟合结果预言其分支比与标准模型预言基本一致误差范围内略有变化与当前实验上限和测量值完全兼容。B → Kμμ 角可观测量如P5‘我们的预言也与当前数据没有显著冲突。这进一步支持了新物理效应主要存在于带电流b→cτν和纯中性流b→sνν过程而对b→sℓℓ这类光子或Z penguin主导的过程影响有限。这些“通过”的检验增强了Scenario III的可信度。它表明一种主要与第三代费米子耦合的新物理可以在解释RD*和可能的B→Kνν增强的同时不与其他低能精密测量产生不可调和的矛盾。5. 方法优势总结、局限与未来展望5.1 机器学习在全局拟合中的核心价值回顾整个工作ML技术的引入为我们带来了传统方法难以企及的优势效率的飞跃一旦ML模型训练完成评估数百万个参数点的似然值只需秒级甚至更短的时间。这使得高分辨率、多维度的参数扫描和复杂的关联性研究成为可能。而传统的网格扫描或MCMC方法完成同样任务可能需要数周甚至数月的计算时间。对非高斯似然的精准模拟我们的χ²景观是高度非高斯的存在多个“脊”、平坦区和强关联。XGBoost这类树模型非常适合捕捉这种复杂结构如图11所示其模拟的置信区间轮廓与真实计算高度一致。稳健性与可重复性相比于深度神经网络树模型对超参数和训练数据量不那么敏感训练过程更稳定结果可重复性更高。这对于需要作为可靠工具嵌入物理分析流程的方法至关重要。内在的可解释性通过SHAP分析我们不仅可以得到最佳拟合点还能理解每个参数如何影响整体拟合。例如我们可以量化地看到当C3偏离其最佳值时主要是哪些可观测量如RD*的χ²贡献在急剧增加。这提供了宝贵的物理洞察而不仅仅是数字结果。5.2 当前框架的局限与改进方向当然我们的方法并非完美仍有改进空间训练数据的成本生成初始的训练样本5-10万个点仍然需要大量的CPU时间进行完整的理论计算。这是目前的主要瓶颈。未来的工作可以探索更智能的主动学习Active Learning策略让ML模型在训练过程中自己判断哪些区域的参数空间不确定性最大然后优先在这些区域请求新的计算样本从而用更少的计算量达到相同的精度。向更高维度扩展目前我们聚焦于3个核心参数。但完整的SMEFT分析可能涉及10个以上的威尔逊系数。随着维度增加维持采样密度和模型精度会面临挑战。可能需要结合降维技术如PCA分析重要方向或使用更适合高维度的神经网络架构如贝叶斯神经网络与树模型进行融合。系统不确定性的更精细处理目前的理论误差是通过在采样时随机抖动相关参数如形状因子参数来近似的。更严格的方法是构建一个完整的“理论误差模型”并将其参数也纳入拟合但这会进一步增加维度。ML模型需要能够处理这种额外的随机性。5.3 对未来实验的启示我们的分析为未来的实验寻找新物理指明了几个关键方向精确测量 B0 → K*0νν当前B→Kνν和B0→K0νν的数据似乎存在一些张力后者上限较低。如果Belle II最终确认B→Kνν的显著增强而B0→K0νν依然被压低那将强烈暗示我们需要超出当前最小参数集C1, C3的新物理例如包含右手流贡献的算符或者更复杂的味结构。这将是对新物理模型构建者的一个严峻挑战和重要线索。提升 RJ/ψ 的测量精度RJ/ψ (Bc→J/ψτν) 是另一个重要的带电流检验道。目前其误差还很大。如果未来LHCb或升级后的探测器能将其测量精度提升到与RD相当的水平它将为我们提供关于新物理耦合到重夸克系统的独立信息并与RD形成交叉检验。寻找更多的τ可观测量例如测量B→Dτν衰变中τ轻子的极化或者寻找极其罕见的B→Kττ衰变。这些测量能帮助我们区分不同的新物理模型是标量、矢量还是张量耦合从而在解释RD反常的众多模型中做出甄别。机器学习在粒子物理全局拟合中的应用已经从一种新颖的尝试逐渐发展成为一项强大且不可或缺的标准工具。它不仅仅是一个“加速器”更是一个“显微镜”让我们能够以更高的分辨率和更深的洞察力去审视实验数据在理论参数空间中留下的细微痕迹。随着更多高精度数据从Belle II、LHCb以及未来的高能对撞机中涌现这种数据驱动与理论指导相结合、传统物理计算与先进算法相融合的研究范式必将引领我们更接近那个隐藏在标准模型背后的、更基本的物理定律。