[论文解读] Monte Carlo simulations of the electron-gas interactions in the KATRIN experiment

J. Kellerer, F. Spanier|arXiv (Cornell University)|Dec 31, 2021

Atomic and Molecular Physics参考文献 28被引用 3

一句话总结

本文提出了一种新的蒙特卡洛模拟代码 KARL，用于模拟 KATRIN 实验中无窗气态氚源（WGTS）内的电子-气体相互作用，通过追踪电子和离子以计算能量谱、粒子密度和电流分布。其主要贡献在于实现了电子与离子动力学的自洽、空间分辨模拟，揭示了显著的径向离子电流以及与位置相关的谱特征，这对于等离子体势建模和中微子质量探测灵敏度至关重要。

ABSTRACT

At the KATRIN experiment, the electron antineutrino mass is inferred from the shape of the $\beta$-decay spectrum of tritium. Important systematic effects in the Windowless Gaseous Tritium Source (WGTS) of the experiment include the energy loss by electron scattering, and the extended starting potential. In the WGTS, primary high-energy electrons from $\beta$-decay produce an extended secondary spectrum of electrons through various atomic and molecular processes including ionization, recombination, cluster formation and scattering. In addition to providing data essential to the simulation of energy loss processes, the electron spectrum also provides information important in the simulation of plasma processes. These simulations will then provide an insight on the starting potential. Here, a Monte Carlo approach is used to model the electron spectrum in the source for a given magnetic and electric field configuration. The spectrum is evaluated at different positions within the WGTS, which allows for a direct analysis of the spectrum close to the rear wall and detector end of the experiment. Alongside electrons, also ions are tracked by the simulation, resulting in a full description of the currents in the source.

研究动机与目标

• 开发一种自洽的蒙特卡洛模拟方法，用于模拟 KATRIN WGTS 中带电粒子的动力学行为，包括电子和离子。
• 解析源内部空间依赖的电子能量谱、粒子密度和电流分布。
• 提供准确的粒子电流估计值，特别是径向离子电流，这对于建模等离子体势效应至关重要。
• 通过包含完整的粒子追踪、次级过程和场耦合，改进以往的模拟方法。
• 支持与总电流和径向电流实验测量结果的直接比较，助力等离子体势重建。

提出的方法

• 采用蒙特卡洛方法模拟在 2.5 T 磁场和电场构型下电子和离子的单个轨迹。
• 使用更新的截面数据建模电子与氚分子的相互作用，包括电离、复合和弹性散射过程。
• 追踪次级粒子，并通过自洽的密度场反馈其对局部等离子体密度的影响。
• 在多个径向和轴向位置设置虚拟屏障，以计算粒子电流密度和能量谱。
• 利用洛伦兹力作用下的运动方程积分粒子运动，并包含弹性散射引起的回旋运动效应。
• 通过已知物理极限和先前结果对模拟进行验证，尤其关注电子能量损失和谱形特征。

实验结果

研究问题

• RQ1 电子能量谱在 WGTS 内部的空间分布如何变化，特别是在后壁和束流管附近？
• RQ2 电子和离子的径向与轴向电流密度是多少？它们如何随磁场强度和源条件变化？
• RQ3 次级电子过程（电离、复合、散射）在多大程度上改变了源中的主电子谱？
• RQ4 粒子密度和电流分布如何依赖于 WGTS 中入口气体密度和温度？
• RQ5 该模拟能否重现可测量的电流（如后壁和 DPS 电流），并预测不可测量的组分（如径向离子通量）？

主要发现

• 模拟的电子谱表现出三个显著区域：β衰变区、热化区和次级电子区，各区域的相对贡献随源内位置而变化。
• 朝向束流管壁的径向离子电流不可忽略，当磁场为 2.5 T 时，在距后壁 5 cm 处最大可达 8 × 10^6 cm⁻²s⁻¹。
• 离子密度分布表现出强烈的径向依赖性，向束流管壁方向增加，当磁场为 2.5 T 且入口密度为 6.3 × 10¹⁴ cm⁻³ 时，峰值密度超过 10^6 cm⁻³。
• 径向电流和密度分布对磁场强度敏感，低磁场（1 T）下电流减小，原因在于回旋半径增大及散射引起的径向扩散增强。
• 模拟揭示，电子能量损失和次级电子产生显著改变了后壁和束流管附近（尤其是低能区）的谱形。
• 该代码成功再现了已知的谱特征，并通过引入完整的离子追踪和空间分辨诊断，扩展了以往模拟的范围，实现了与未来径向电流测量的直接对比。