胡晓凡(00:00:08):
咱们开始吧！30点下开直播到我就这么说！二点半开始报告。整了一个 indico，谢谢老师！Indico cdex 网页。本次会议的。

胡晓凡(00:01:06):
二点半。12点半吧。

胡晓凡(00:01:24):
行了。

胡晓凡(00:01:47):
是领导。你现在是很高的？然后滴滴也对我会将 PPT 传上去。因为这个东西是有摄像头的奖学金什么的，发个有的技术有的会议记录有。

胡晓凡(00:04:42):
我们开始吧，各位老师，各位同学大家下午好，今天是我们本学年第一次的这个 thu seminar 就是以前的这个大组会，可以看到相比以前的大组会，它有两个区别，一个，它就是由以前的一周一次改成了两周一次，再它就是这个组会的形式就更正式，我们今天的两个报告分别是对撞机组的，从 ds02317到 bs05700近阈值效应与重轻强子谱学而中微子的报告是。

胡晓凡(00:05:27):
反应堆中微子能谱解谱方法简介每个报告安排是20分钟，再加后面的10分钟提问，现在有请卢思源同学来为我们带来对撞机组的报告。你可以用键盘，我会及时。【卢思源】大家好，今天我想从 ds0222317出发。讨论新出现的 bs0star 5700候选态以及它背后的精神动力学，这两个它有一个比较值得比较的特点，它们都位于两个位于相关的两个强子阈值以下，而且宽度很窄。接下来我会先介绍实验结结果与对称性分类，在讨论阈值效应怎样影响质量与宽度，最后说明需要哪些观测量。才能进一步区分不同的理论解释。

胡晓凡(00:06:38):
首先看这次的分析结果，lhcb 利用今年8月份 lhcb 利用 run one 和 run two 总计。有个 in inverse 被罚的数据，在 bs0派零的质量分布中观察到了一个新的窄结构，整体显著度超过七倍。西格玛。它的质量约为5699mev 宽度目前给出的上限在90%的置信水平下小于9.8个 m0，这里不能把焦点8mb 当成此的宽度，因为实验还没有显著分辨出它的本征宽度。

胡晓凡(00:07:20):
对于今天的讨论，最关键的是峰与 bk 阈值的位置关系最低的 bk 阈值大约在5773MV。这个结构比预期低了大约74开米，这意味着它不能通过产生真实的 B 和 K 介子发生衰变。如果把它解释，我们可以看到这个 bsbs0和派零他们两个都是零负的，如果他们从最低阶的 S 波衰变的话。那么 ds05700它就是一个它也它是 jp 等于零正的认真的太，但是如果把它解释为这个 cp 等于零正的 ds 态。就会出现两个相互关联的问题，那为什么坐在阈值以下这个位置又怎样解释它的窄宽度，后面我会以 gp 等于零正的假设作为讨论前提。不过量子数仍然需要独立的实验检验。为了回答刚才的问题，我按这个顺序展开，我先用 hqss，也就是重夸克自旋对称性。说明这个态在重氢介子谱中应该处于什么位置，然后回顾 ds40422317，它在 csb 体系中已经提出类似的问题，可以作为理解 bsr 体系的参照。接下来讨论耦合道与阈值效应，尤其是为什么一个已经关闭的睡眠道仍然能够改变质量。为什么 S 波的阈值效应特别突出？最后再把这些动力学效应与守正伴随机制联系起来。讨论怎样利用实验结果的约束约来约束不同的理论图景？

胡晓凡(00:09:13):
先介绍后面需要用到的 H two ss 一个重轻介子包含一个重夸克以及围绕它运动的轻自由度。当重夸克质量远大于强相互作用的典型尺度时。重夸克自旋方向对轻自由度运动的影响会比较小，因此在重夸克质量趋于无穷大的极限下。可以把重夸克自旋与轻自由度分开处理。从从从这个展开式可以看到最低阶的最低阶不显示依赖中华和自旋到了 mq 分之一阶。动能项仍然与自旋无关，而四磁项开始显示一代自旋，从而破坏 H two ss 对称这个对称性，同一个二重态内部的质量劈裂就与这类修正有关。由于修正受到重夸克质量倒数的抑制，hqss 在 bsr 体系中通常比在 CSR 体系中更加准确，它对谱学分类的作用是我们可以先确定轻自由度。

胡晓凡(00:10:24):
对总角动量 gl，再与重夸克的自旋1/2耦合得到整个介子的总角动量。对于今天关心的 P 波激发，这 L 可以是1/2，也可以是3/2前一种形成。零正和一正的二重态后一种形成一正二重态。虽然两组里都有一正，但它们的轻自由度结构不同，舞蹈衰变方式也不同。

胡晓凡(00:10:51):
下一页就看这个区别，对于 jl 等于1/2这一组零正和一正可以分别通过 S 波耦合到 bk 和 B star k。如果通道开放相空间足够，这类衰变通常会带来较大的宽度。对于 gl 等于3/2这一组相关衰变主要通过地波进行。它在阈值附近受到更高次动量因子的抑制，因此通常比较窄。这里需要区分两个角动量概念，介子内部的 P 波激发与衰变末态之间的相对分波是两回事。内部属于 P 波的介子，完全可以通过 S 波衰变，所以对于一个零正态，不能只根据量子数判断它一定很宽，还必须看它处于。

胡晓凡(00:11:41):
Bk 阈值的上方还是下方，如果低于阈值，原本可能很强的 S 波衰变。就会被运动学条件关闭。这个情况在 cs8体系中有一个典型例子，就是 ds042317它的质量约为2317个 MV 位于 bk 阈值以下大约40个 MV，这里是一个它的示意图。所以由于它的位置位于这个 dk 阈值之下，所以说通常一些很重要的 dk 强势频道是关闭的，让你观察到的 ds 派零衰变涉及的。同位旋破缺会受到抑制，这为它的窄宽度提供了自然解释，我们还要进一步问它的质量为什么会落到阈值以下？要不再卖。

胡晓凡(00:12:53):
没电了有。对。但我们还要进一步问它的质量为什么会落在阈值以下，这也是态发现后受到广泛关注的原因，它的质量低于一些传统夸克模型的预期。说明附近 dk 道的作用可能不能忽略这里的异常低是指相对于具体模型而言的，在耦合道的描述中。一个紧致的 cs 阀核心可以与 dk 连续谱发生耦合，最终的物理质量会包含这种耦合带来的修正，而不是直接等于没有计入这部分作用时的质量参数。因此，我们希望在同一个框架中理解它的低质量阈值下位置以及窄宽度新的 bs 结构为我们提供了一个机会。

胡晓凡(00:13:54):
可以检验类似的动力学机制，在 bs8体系中是否仍然成立。回到 bs 过去已经确立的 bs 义务585830和 bs2star5840通常归入 jl 等于3/2的二重态。它们都是比较窄的态，与前面讨论的第一步衰变特征一致。相比之下，gl 等于1/2的零正一正态，二重态长期缺少明确的实验确认。这次出现的 ds0star 5700为其中的零正态提供了一个自然候选，如果这个解释成立，下一步很重要的任务就是寻找同一二重态中的。一生伴随态这里的一证要并不是已知的 dse555830它们的总总的 gp 是一样的。自夸克及线下，它们属于不同的轻自由度分类，找到这个大水菜之后，我们才能对单个候选态的解释，推进到对整个二重态的检验。

胡晓凡(00:15:03):
现在可以很具体的回答这个零正候选态为什么能够很窄，如果它的量子数确实是零正，那么它的 bk 的衰变。本可以通过 S 波发生这个耦合本来有可能带来很大的强衰变宽度，但是它的质量低于 bk，阈值大约74个 mev 没有足够的能量产生真实的 B 和 k。所以这个设备到实质上会被关闭。与此同时，实验看到的 bs 派零衰变，对于通常对于我们假设的零正的这批假假设来说。涉及同位旋破缺，因此受到了抑制。这就说明载宽度并不能直接推出 bq 和很弱，即使耦合很强，只要没有相空间，这个道也不能贡献真实的衰变宽度。因此质量，位置和射量宽度必须结合形态，不过这里我还有一个问题，bk 到了不能产生真实衰变，是否意味着我们可以在理论中忽略它？下一页说明答案是否定的。先看这页上面的传播子 m0表示模型中没有计入这部分耦合修正时的质量参数。

胡晓凡(00:16:26):
西格玛 bk 表示 bk 到带来的自能修正，物理态对应的破对应的极点由传播子分母为零的位置确定这个极点能量可以写成。表达式，E 等于 M 减2分之 I 伽 I 乘伽马，因此极点的实部对应质量 M 虚部对应负的半宽度。再把智能分成实部和虚部两部分就能理解他们为什么分别与质量修正和宽度有关。智能会改变负极点的位置，从而影响它的。实部和虚部，当能量低于 bk 阈值的时候，bb 和 K 不能同时成为真值的在壳中间态。因此，在这里采用的近似下 PK 到在阈值下的10能量处的智能吸附为零，贡献对应的衰变宽度。但是虚的 bk 中间过程仍然存在，智能10步一般不会消失，质量仍然会受到修正，同时其他开放到仍然可以使这个态就。非零的总宽度，因此一个很窄的态仍然可能受到显著的耦合到质量修正，具体移动多少向哪个方向移动。需要结合相互作用和重整化方案进行计算。如果修正很重要，紧致核心和 dk 连续谱都需要纳入描述。这些只比较不同分波在阈值附近的行为，首先看下面的动量关系，在量体阈值附近相对动量的大小正比于能量与阈值之差绝对值的平方根。

胡晓凡(00:18:17):
这个平方根意味着阈值处存在特殊的非解析的结构，也就是这里说这里标出的 branch point，它这个非解析的意思是它不能在0点处。做泰勒展开，再看下面的分布规律，阈值相关的非解析，其贡献对于角动量为 L 的分波带有 P 的22l 加一次方的动量因子，因此 S 波对应 p。

胡晓凡(00:18:50):
P 波对应10的立方 D 波对应10的5次方当动量趋近于零时，高斯幂会更快的变小。这里的大小关系是在。同一典型尺度下比较小动量展开后的抑制程度，具体贡献还要乘上各自的耦合系数。所以 S 波保留的最明显的平方根型阈值效应更容易使线型发生明显变形。或者形成一直一直尖点，这并不保证每个 S 波阈值附近都会产生心态，向他说明这类阈值附近的动力学不能轻易忽略。考虑耦合到之后，对物理态组成的描述也会更丰富，我们可以把这个态理解成紧致的 bs8核心与 bk 连续谱成分的叠加。随着相互作用参数改变两部分的相对重要性也会发生变化，这些从左到右表示的就是从紧致成分较重要，经过明显混合的再到 bk 成分较重要的不同情况。下面列出的 ds0star2317的组分比例来自于某一个耦合道分析中的例子，用用以说明一个物理，它也可以同时包含多种成分。

胡晓凡(00:20:12):
不过这些比例很依赖具体模型，自由度选择定义不能直接当成实验测出的概率也不能直接搬到 esr 体系。因此说态存在紧致核心。同时又具有显著的两强子成分并不矛盾，真正需要进一步追问的是，这种组成判断对于相关道的耦合散射行为和衰变过程能给出哪些可以检验的结果。除了阈值动力学，还可以从守正对称性的角度理解这些态之间的关系，前面的 hqss 把相同与自旋。

胡晓凡(00:20:53):
先相同语圈自旋等于零和自旋等于一的态联系起来，例如零零-1负构成一组零正一正构成另一组手中伴随态伴随机制的进一步讨论这两组。相反，宇称的二重态之间是否存在由轻自由度决定的。决定的关系，在这一机制中正负宇称之间的质量差与轻自由度的自发守恒对称破缺有关，如果这个基是占主的，那么宇称劈裂对中文。对中位夸克种类的依赖可能相对较弱。也就是说，从 cs8体系换到 ds8体系以后，这些质量差仍然处在相近的尺度上。不过。不过这需要，这是需要检验的动力学判断，实际的物理质量还会受到有限重夸克质量修正以及附近耦合道的影响。这些这样可以具体比较一下信息化体系和运算化体系的质量差。先看 cs8体系，ds042317和 ds 的质量约为。

胡晓凡(00:22:06):
质量差约为349mev，而 ds12460ds2的质量差约为348mb2者非常接近。如果把新的 bs04点5700作为零零正态候选。那么它与 bs bs 的质量差约为332mev 也处在相近的尺度上。这个规律值得从守正伴随机制的角度研究，但仅凭这些质量差，还不能确定唯一确定动力学机制。原因是 bk 和 B star K 两个偶道之间也受到 hqss 的约束，相应的阈值效应可能产生相关的数量移动，从而给出类似的序列规律。因此这张表的最后一行不重要 esr 体系中的一正半学态还没有得到确认，只有找到它并精确测量质量才能比较两组宇称奇点以及正宇称二重态内部的自旋奇点对理论提出完整的约束。前面介绍了几种理论图景关注的物理问题有所不同，所以需用来需要用来检验它们的观测量也不完全相同。这里在这里由时间关系我就略过了，到时候也可以把今天的报告归纳一下。最主要的结论是，域下窄态仍可以具有显著的耦合到效应。新的 bs 结构与 dsn star2317。他们具有值得研究的对应关系，二者都位于相关的 S 波阈值以下主导强衰变通道关闭，因此可以表现为窄态，但关闭衰变道并不会消除虚过程带来的。

胡晓凡(00:23:55):
质量修正，理解它们的质量位置时需要认真考虑连续辅助作用 hvss 和守征矩阵态之间的关系提供了分类依据。而具体的动力学解释还要通过更多观测量来检验，其中寻找一正伴随态，并联合测量质量，宽度和耦合是进一步区分不同解释的重要方向。

胡晓凡(00:24:22):
还有行，那我接着说，最最后做一个方法上的延伸讨论这类强子谱学分析和中微子实验之间有哪些可以交流的经验。首先两类分析都需要从明确的物理假设出发。在 lhcb 谱学分析中，我们通过比较不同的振幅和共振结构，在中微子实验中则可能比较不同的质量参数相互作用模型。或者新的物理假设都需要把这些假设转化为探测器，实际能够观测的分布对 ltcd 需要从振幅和得到。得到质量和角角分布，同时考虑干涉接收度，分辨率以及背景对于中微子实验需要把束流截面和探测器响应结合起来。得到重建后的事件分布，接下来通过联合似然进行统计推断，把系统误差纳入分析，同时检查参数相关性和退化方向。最后还需要利用替代参数法控制样本和模拟实验检验结论是否稳健，当然两类分析的具体难点不同，强子振幅分析特别关注相位干涉以及。基点的解析，延拓中微子分析则经常面对核效应能量重建和震荡参数之间的检定，那他们有一个共同需求模型，不仅要描述一个局部分布。能够解释相关样本，并经得起系统误差和分析方案变化的检验。我的报告就到这里，谢谢大家！

胡晓凡(00:26:03):
正好20分钟！这边今天由于王亚朋同学他在出差，他在阅读摘要之后委托我提问一个问题，就是在这个第二页。这个 bs0star5700这个态在分析过程中，该峰的真伪是怎么辨别的，还有它的显著度和系统误差是怎么估计？这个我没有太了解布氏通风，这个我相信需要在如果说这一页有一个引用文献的话，我们去就会引用文献对我们去引用，去点击它引用的地方就能知道更多信息。请问大家还有什么问题？如果没有的话，我提一个就是第5页，它主导的衰变是 S 波衰变和 D 波衰变。它能有 P 波衰变吗？它允许吗？博是可以报说它存在的是不主要？我不太清楚，我觉得应该找到这个东西的源，我们再去看，谢谢！

胡晓凡(00:27:42):
还有你说要，那你怎么表达到底？24层以上。用用用什么道具扫对你你不是要你说，关键是找到他们，从从那个理论模型的角度找到一阵会就很有帮助，但是怎么找？不说有找动机很强，但是具体怎么讲我不是很清楚。我问你不，这个是那个，这不是我做的是个公开结果吗？这是已公开结果？从后分对，不是我做的分析是讲的你调研好。好好好，确实不要看没关系，还有刚才你我刚才。好家伙关于 B wps70一些细节看一下做了？好像他好像没有论文，他有一个论文好像没有公布吧，就是他有一个新闻，还是随便那个官方信不信，现在已经使用有，或许我们可以找到，如果有新闻的话，它一定是在一个什么会议上。如果能找到会议 PPT 的行知道行，我太可怜了，我之后看。

胡晓凡(00:29:37):
大家还有好的，谢谢四院长，我们下一个报告中微子组的报告人是孙浩哲同学，我们需要特别注意，我们这次的一个改变就是每一个报告都。都需要有一个提问人，请问刘继安同学在场吗？好的，请您认真听报告，待会儿您将作为这个报告的提问人好的用用对用用上下翻页就行。OK 好了。好【孙昊哲】大家好，我给大家简单介绍一下，前几天看了一篇这个大湾解谱的工作，然后它当然里面包括了一些这个大湾在2021年发表的一些独立解谱的工作，以及在2025年发表的这个联合解谱的工作。这是今天的大概的这个大纲，首先是 motivation 然后介绍用到几种解谱的方法，简单介绍一下在大家完成上实现的过程，以及最后一个总结。

胡晓凡(00:30:47):
首先是对大家班做一个简单的介绍，大亚湾是一个中微子反应堆中微子实验，短基线反应堆中微子实验，主要目标是测量中微子混合角 cd3。它的反应堆有6个反应堆堆芯，然后是商用的压水堆，然后探测器有四个近点探测器和四个远点探测器。近点探测器的位置，大概离这个堆芯的距离是560m 或者到600m，这些是加权过的这个距离。然后所以近点探测器能够在这个大的震荡发生之前，精确比较精确测量原本的反应堆的中微子能谱反中微子能谱。

胡晓凡(00:31:25):
然后就是 motivation 就是我们为什么要做一个解谱工作，因为就反应堆中微子而言，反应堆中微子，目前看它的测量的结果与理论的模型之间有存在不符，然后比如通量不符以及谱型的不相符。所以我们实验中只能测到每个探测器的探测谱，如果我们想在不同实验之间做个对比的话，可能是需要一个反推回去一个中子能谱。可以方便对比，然后另外其实这个解谱工作不在于这个反应堆中微子，或者说甚至不只在于物理中存在，它是一个很常见的场景，就是一个原始的信号，经过一次的响应之后，它会变一个变成另外一个形状。除此之外，再加上一些噪声，最后成成为我们真正的测量结果，然后我们需要从这个真实从测量结果反推回真实的谱。这是我们的一个动机好，然后简单介绍一下几种。

胡晓凡(00:32:20):
使用到的 unfolding 的方法。首先看我们的问题，还是由一个真实的信号 S 经过探测器响应 R 以及加上一个噪声之后，它就变成了一个测量谱 p。然后我们如果从贝叶斯的这个概率的角度来考虑，其实这个响应矩阵就相当于是一个这个条件概率，然后我们把这个公式写出来，就发现我们主要有一个前沿的这个公式前沿的。先验的这个模型，然后我们就可以通过这个响应探测器响应，得到一个后验的结果，然后我们再把这个后验结果反反过来，再带回到相应的模型里就可以不断迭代。然后最终迭代的结果就是会使我们会预我们预期它会与真实的结果相近，这是这就是贝叶斯迭代的这个思想。

胡晓凡(00:33:11):
但是，这个理念是一样的，但是具体的这个计算方法有很多，比如说在这里具体的一种这个 rl 的迭代方法，它公式具体就这么写，然后不再具体解释了，其中比较重要的一点是它是有一个比较特殊的参数就是只有一个参数迭代的次数 k。当这个 K 趋于零的话，其实就是把先验的概率就其实我们输入的先验概先验的这个模型，而当这个迭代次数。

胡晓凡(00:33:43):
或者较少的时候，它肯定是受先验的输入影响会很大，但是当迭代次数多了之后，在这个数学上可以证明它的迭代结果会收敛到最大似然估计。那我们的估计的参数是这个一个信号的每个并的大小，那如果它收敛到最大似然估计的话就是最终的这个迭代的无限次结果就是把我们的谱原本拿出来了，所以这也是我们不需要的，所以我们要在这个迭代次数之间做一个取舍，找到一个合适的迭代次数，然后迭代次数选择之后我们会再简单介绍。这是从贝叶斯的角度考虑，那我们就从把它写成矩阵的语言，从最简单矩阵方向来想，其实这个问题最直直白的想法就是我既然知道一个。响应的矩阵 r，然后那么我们就希望直接把 R 对它求逆，然后就直接从这个我们测量结果反推回这个真实的结果，比如说我们把这个 R 逆写出来，大概就是一个。在这么一个形式，但这里面这个 R 转置乘 R 的逆它转置乘 r，它并不是一个性，不一定是一个性质很好的矩阵，它里面可能有很多非常小的特征值，如果我们直接取逆的话，这些它的就会把某些模式的噪声，经过放大，就不会得到真实的谱形，所以我们的想法就是忽略掉这些小的模式，而保留那些真正与物理相关的大的模式。

胡晓凡(00:35:21):
这个想法就是通过奇异值分解，然后忽略掉小的奇异值就是 svd 方法。这 svd 分解任意的实数矩阵都可以拆成两个正交矩阵，这个 U 和 v，然后以及乘一个对角矩阵这个 sigma 的乘积。其中的 UV 是正交矩阵，然后西格玛是对角矩阵，然后这个对角元，这个西格玛一到西格玛 N，我们习惯上把最大值放在前面，然后我们把这个 V 移到左边来，在这个 S 空间就是真实的物理信号的空间里的任意一组。一组模式 V 它就会对应到我们探测到探测空间的一组模式 U 那么我们的现在目标就是想要把真与有真实物理意义的这个 V 模式挑选出来，然后而把没有。把噪声对应这种无意义的模式给它去除掉，就是我们如何选取正则化的过程，很简单，想法就是说我们把每一个。其值给它加上一个权重，然后保留那些重要的，然后忽略那些比较小的那正则化的过程，当然有很多种，比如说第一种直接做一个截断，我们只选取前 N 个大的奇异值，然后保留到，然后去除掉其他的所有奇异值。或者说用一种比较平滑的方法，把这个权重写成一个，其值与人为的惩罚因子，这个 tau 的1个×1，一个相除的结果。

胡晓凡(00:36:50):
当然还有很多其他的正则化的方法，还有就是我们的 winner svd 它其实严格意义上并不算一种正则化的方法，它是让它不再是人为去给选取一个正则化的参数，而是把用用于基于这个信噪比的这个想法，然后把这个噪声在某一个模式上的。这个方差给提取出来，然后比较奇异值与这个方差大小，然后来这个保来看是否保留这个结果就简单介绍，简单对比一下，其实就是贝叶斯方法，它的机制是通过提前终止迭代过程，然后得到我们想说想要的结果，然后正则化方式，这个 svd 方法都是通过，选取特定的模式保留某些奇异值，而。去除某些模式，他们都需要一些参数来决定最终的结果，比如说迭代方法，它的参数就是迭代次数 n，然后这种截断就是选取截几个，然后其他这种惩罚项的方式就是惩罚项的大小是为人为来确定的，而至于这个 yod svd 方法，它不是人为确定的，而是，去计算噪声的方差。

胡晓凡(00:38:04):
不同方法它的参数选择参数的优化也会有不同，比如说前三种的话，就需要人为去优优化这个参数就一般需要用这个做蒙卡产生大量的模拟数据，然后来从模拟中研究。哪一种参数会比较适合，当然也对于是否有输入的依赖也不同，比如说贝叶斯的话，方法就会依赖于这个输入，而其他方法不依赖于这个。先的输入好，然后那我们要优化这个参数的话，我们就想首先知道用什么来衡量一个解谱它的结果是好还是坏，我们用的方法就是这个均方差。它的这个它偏差的平方，加上这个统计方差，所以对于这个贝叶斯方法，我们就需要选取这个与实验相同条件产生大量的蒙卡数据，然后来然后选取不同迭代次数来对比每次不同迭代次数。产生的这个均方差的大小，但做的时候有一些注意的要点就是我们输入的模型与我们产生就是我们作为输入的先验的模型，与我们产生蒙卡的真实谱是不能相同的。这个选取谱的时候也要选取输入模型的时候，也要根据这个我们对于输入模型的了解程度，根据它对它的方差来进行选择，而对于其他的这些 svd 的截断方法，惩罚项的方法，也跟这个相同就是产生大量的蒙卡数据，然后不断的扫描参数找到最优势的结果。

胡晓凡(00:39:40):
而对于 wind svd 方法，刚才也介绍了它的原理的是计算某一个模式，i 项上这个方差的大小。就方差投影在这个模式 I 上的大小是多少，然后用这个权重因子来压低这部分，当然实际的做的过程中，一般来讲都是把这个方差投重新投影到这个协方差矩阵上。这些都是数学上的处理就是不太我要追溯那上面就简单讲了几种所谓的定的方法。其中就是贝叶斯迭代的方法，然后 svd 正则化的方法以及 W svd 的方法，但实际上正则化的方法不止这样几种。刚才说这些情况下，其实可以人为自己选择，可以选择去抹平到一些这个，比如说可以针对性的抹平到二阶导或者抹平到一抹平一阶导等等。

胡晓凡(00:40:36):
根据我们需要的模型去做选择，然后简单介绍一下在大湾上，我们实际工作在大亚湾这种反应堆中微子实验中需要。需要考虑的一些过程，大湾是一个商用的一家水堆，然后它主要是产生的反应堆中微子是由四种核素衰变。就是铀235铀238钚239以及钚241。主要的来源是铀235和钚239，它们大概占了百分之大概九十左右。然后这个衰裂变产额，由于这个钚241是钚239产生的，所以。

胡晓凡(00:41:15):
就是这个8241与8239的身身份份额基本上近似于正比，所以我们没有办法把它们解耦开，所以我们只能解解一个谱是铀235的谱，以及这个钚二钚241和钚2391个线性组合的谱，就是这个钚的 combine 谱。对于这些小的谱小的影响，比如说铀238和钚二三。我波224是波241的谱，我们就把它作为约束项放到拟合里面。然后我们在大湾反应堆，我们能记录到每每一周的能谱以及每一周的瞬变份的裂变份额，所以我们可以对时间上的所有的这些，每每周的数据做一个整体的拟合来分解得到。235241的谱。OK 然后我们把我们得到了这个铀235的谱以及钚的谱之后，我们还会同时得到它们的协方差矩阵。OK 这是我们通过这种时间上的双裂变份额得到的分离的这个元素的分离的核素的慢中子能谱，然后刚才说到一个问题就是。比如说像贝叶斯迭代方法的话，我们需要一个先验的模型输入。对先验模型输入的话，就需要一个一个对于比如说铀二铀238和钚241的这个先验的谱或者铀235以及钚的这个约束的谱，这些谱是怎么计算的，需要通过这个计算来得到。那这个计算有两种方法，第一种就是这个转换反 conversion model method，它的方法是我们通过磁谱仪的实验，可以直接测量每一种。

胡晓凡(00:43:10):
核素热中子裂变之后，照射热中子裂变之后，所有的裂变碎片累加产生的这个贝塔能谱，然后我们在这个累加的谱里，引入大量这个贝塔谱的分支来拟合掉这样贝塔谱，然后利用这个贝塔衰变的动能关系，反推回这个中微子能谱，然后再把这所有这些拟合出来的结果相加得到这个反的中微子能谱，这就是转化法。那另一种方法就是 summation model 就是加和法，我们就可以把这个核数据里每一条分支比它的每一个衰变的分支比都考虑清楚，然后把每一个衰变的这个中微子能谱都加在一起，产生。

胡晓凡(00:43:48):
最后的总和的谱是两种方法，然后现在我们有了这个先验的输入，有了这个约束的谱，然后我们也通过前面的这个分周能谱裂变份额计算出来的分解谱，那我们怎么如何同时解谱？恒信谱，它的思想就是把这三张谱拼接成一张谱。比如说我如果对于每个谱，比如说1231235的谱或者布的谱或者总总的谱分成25个并的话，我们可以把这三个25拼成一张总共的75个并的大谱。然后同时我们之前得到的拟合过程中得到的这个。协方差矩阵，所以我们现在有一个75的75个并的合成的一张谱，以及75×75的协方差矩阵谱，那么我们就可以在此基础上利用这个贝叶斯迭代的方法，或者说这个 svd 的方法来反解出这来直接反解出这三张谱就是两张铀和钚的谱以及总的。

胡晓凡(00:44:45):
总总谱。Ok. 然后对于选取的方法，对于不同的要求，我们的选取方法也不一样。在单独在大湾2021年单独解谱的时候，发现这个 vana svd 的方法它的均方根是最小的，所以它的效果是最好的，当时也是按这种方法来算的，但实际上解联合谱的时候，这 svd 方，这个 winner svd 方法就会有一些它的缺点。就是 nsvd 分解的时候，它的某一个模式可能不再代表真实的物理结构了。比如说我们解谱的时候，铀235的谱与钚的谱肯定会存在负相关，因为总谱是确定的，所以它们会存在负相关。很有可能在 svd 分解的时候，将某一个模式分解到了这个铀235加布二，加上布的谱对应一个结构，此时他！此时它的权重因子，不但能代表真实的信噪比，所以我们在为了避免这种现象，最后这个大班在最终的。

胡晓凡(00:45:45):
这个解谱过程中还是采用了这个 svd 正则化的方法来解谱好，然后就是最后就是简单总结一下，一方面是这个报告简单介绍一下这个贝叶斯正则化方贝叶斯迭代方法。sv 方法怎么进行这个解谱，然后以及这个正则化方法的选取参数是如何选取的效果是如何评估的？然后。接着就是大班这个文章中，大湾用的这个正则化 svd 的方法，然后利用反应堆堆芯演化与同做同位素分离，然后加联合的解谱，得到的这个最高的。最高精度的反对中子总统，然后发现这个总谱的流量就总谱的通量与这个 ki 模型和 sm2023模型比较吻合，但是它的通量是低于这个 huber mer 的模型。得到的这个谱型与所有的模型均不一致，在5m1附近会发现有一个 B 存在，然后也会验证了这个铀235是通量亏损主要来源好，那今天我的这个报告就这些。

胡晓凡(00:46:54):
他这个时间有点多，我们请刘纪安同学提到两个问题。【刘骥安】我有几个问题，第一个是如果是结果写出来的东西和模型不一致的话，我怎么看？我们这个信号是不一致的地方，到底是这个？你希望模板还是它的规模不一样，就是像我刚才说的就是可以用蒙卡方式来评判你这个解谱工作本身它是否合理，就是可以用蒙卡产生这个假数据，然后通过蒙卡结果看它的均方根是否是比较小的，如果发现这个谱本身就这个解谱方法本身没有问题的话，其他方向，比如说我们对探测器不了解，或者说我们对稀盐概型的了解是错误的等等。这是军工科技？对。不能与生产 C 是相比，就是这样。贝叶斯的迭代方法是需要一个先验的模型作为输入，然后它才能开始迭代。我们在这个这是输入，这是先验的输入模型，另外就是我们做模拟实验的时候，其实我们也不知道真实的这个模型是什么样的，所以我们还有一个模拟用到的输入的模型。他的这两个模型，它是必须要是不能是相同的。

胡晓凡(00:48:51):
大家还有问题？如果没有的话，我提一个别说就是这一页说到它使用的是蒙卡来验证这个根据测量谱反解中微子能谱这个做的对不对，那么有没有一个真实的一个实验的一个？一种能不能在实验上有一种非常能否非常简单的过程，我们可以用它来做一个 calibration 的事情。我想就是做一个你的用一个已知的谱来证明我们这个来检验我们的解谱方法是有效的，对这个我还没考虑过，不过我觉得是可行吧，我再去调一下吧，这个从从这个原理上来讲，它是行得通的，但是。由于统计网络等问题，评判的时候，反正需要考虑一下评判的标准，谢谢我看看桌子。

胡晓凡(00:50:04):
【徐闯】对。迭代是第一次是我的，比如说我们的大大弯的工作，就是它的第一次的输入谱是通过核数据库里的所有的这个反应过程加和作为一个结果。对因为就当然可以对，其实就是它会就是当然肯定是越精确越好，但是由于我们因为做实验之前也不知道具体的结果是什么样嘛，所以。

胡晓凡(00:50:42):
对也到最后收敛对，但是这个迭代方法它的最终判断，不是以收敛为主的，因为对它是以这个均方根的这个标准为为主。哪个。刚才前边都做了有到就是大概就是这公式其实就是 P 等于 R 乘 S 加加加一，它就是把所有的涨落就是统计涨落以及这个系统下涨落都放到这个最后这个一里了，就让这个一代表所有的统计涨落。你在你迭代的过程中，你使用它是不是也有一个？你等你说在迭代过程中，你现在大概算一个平方多少？

胡晓凡(00:51:59):
这个我确实不太清楚，我最近再去调研他对因为因为就是我理解他的想法就是这个一次和就是一个均值是零的一个过程，就是它的期望是零，然后它有一个统计涨落的过程，然后所以他就没在这里去把一次和加到这个。这个迭代过程里面去，我是这么理解的。大家没有问题，那我们谢谢涛哥，我们这次的。这个30秒我们就结束了。下一次的报告人我们暂定是刘继安和王亚红。他这一场。

胡晓凡(00:52:59):
预期？