气体探测器中气体增益波动的显微模拟与统计分析

《Applied Radiation and Isotopes》：Microscopic simulation and statistical analysis of gas gain fluctuations in gaseous detectors

　　Ijaz Ahmed | Muhammad Saqlain | Farzana Ahmad | Jamil Muhammad
物理学系，巴基斯坦伊斯兰堡联邦乌尔都语艺术、科学与技术大学

**摘要**
现代气体探测器的性能受到电荷倍增过程中随机波动的限制。本研究采用高保真微观建模框架（Garfield++和Magboltz），分析了基线三元混合物（97% C2H2F4、2.5% iC4H10、0.5% SF6）中的气体增益波动。我们评估了在降低的电场（E/p）高达450 kV/cm Torr以及间隙宽度从0.01毫米到0.022毫米的条件下，相对方差（f）和第一个Townsend系数（α）。系统不确定性分析表明，非电离非弹性碰撞是导致增益随机性的主要因素。重要的是，在高电场下相对方差趋于稳定；当E/p > 250 kV/cm Torr时，f与间隙宽度无关，稳定在0.3到0.5之间。此外，该研究还量化了电离机制与激发机制之间的竞争关系。电离产额模拟显示，从指数Furry统计分布到峰值Polya分布的转变是由高电场强度下的电离碰撞饱和引起的。结果被综合成一条通用增益曲线，将f与电离效率（α/E）相关联，为探测器分辨率提供了一种可靠的预测模型。这些发现为在高通量辐射环境中（如大型强子对撞机LHC所遇到的环境）转向环保、低GWP气体混合物奠定了重要基础。

**引言**
气体辐射探测器的发展是实验核物理学和粒子物理学史上的一个基石（Stephan, 1993; Sauli, 2014）。自盖革-米勒计数器问世以来，这些仪器已经演变为能够进行高速跟踪和精确能量测量的高度复杂系统（Knoll, 2010; Blum et al., 2008）。这一发展过程中的一个关键节点是1968年Georges Charpak发明的多丝正比室（MWPC），它实现了在复杂多轨迹环境中的电子数据采集（Charpak et al., 1968; Charpak and Sauli, 1984）。在现代高能物理学中，由于气体探测器的辐射耐受性和成本效益，它们仍然不可或缺（Rigoletti, 2000; Veenhof, 1998; Sauli, 1997）。本研究探讨了电荷倍增的随机性质，这是探测器性能的基本限制。

当入射的电离粒子穿过气体体积时，检测过程开始，粒子主要通过库仑相互作用与气体分子相互作用而损失能量（Martin, 2013; Tavernier, 2010）。初级离子对的平均数量n?由沉积的总能量E和产生离子对所需的平均能量W值决定（Saha, 2013; El-Banbi et al., 2018）：n? = E/W。这种电离的统计分布受到Fano因子F的限制，Fano因子解释了电离事件之间的相关性（Fano, 1947; Alkhazov, 1970）：σ2 = F * n? = FEW。对于惰性气体和常见淬火剂，准确的W值对于预测固有能量分辨率至关重要（Saha, 2013; Morrissey, 2009）。初级电子在外加电场E的作用下向阳极漂移，而离子则向阴极移动（Lippmann, 2003）。在靠近细阳极线的地方，电场强度急剧增加。当电场强度足以引起二次电离时，会发生Townsend雪崩（Davydov, 2004; Legler, 1961）。电子数量n的增长由第一个Townsend系数α描述（Davydov, 2004; Kowalski et al., 1992）：n = n? * exp(α(x)dx。对于圆柱形几何结构，气体增益M可以用Diethorn公式表示（Diethorn, 1956）：lnM = Vr * ln(b/a) * ln(2ΔV) / (lnV / paln(b/a) - lnK)，其中V是电压，a和b是电极半径，p是压力（Charpak and Sauli, 1984; Kowalski et al., 1992）。

气体探测器的响应可根据电压范围进行分类，如图1（CERN）所示。在正比区域，信号高度与初始能量沉积成线性关系（Sahoo, 2015; Lamarsh, 1983）。随着电场的增加，探测器进入盖革-米勒区域，在该区域比例性不再成立（Knoll, 2003; Knoll, 2010）。尽管平均增益是可预测的，但雪崩大小的波动限制了能量分辨率（Demir and Tapan, 2004; Schindler, 2012b）。这些波动源于路径长度和电子能量分布的变化（Schindler, 2012b; Biagi, 1999）。Garfield++和Magboltz仿真工具包用于进行电子传输的微观跟踪（Schindler, 2012b; Veenhof, 1998）。对C2H2F4和iC4H10混合物的模拟提供了现代探测器配置中雪崩统计的详细分析（Schindler, 2012b; Rigoletti, 2000）。

**主要目标**
本研究旨在严格分析气体探测器中的气体增益波动，重点是从传统的宏观建模转向高保真微观仿真。通过使用Garfield++和Magboltz，我们量化了Townsend雪崩过程中的随机变化，并确定了这些变化对能量和空间分辨率的影响。分析的很大一部分集中在复杂气体混合物（尤其是四氟乙烷（C2H2F4）、异丁烷（iC4H10）和六氟化硫（SF6）在现代探测器架构（如电阻板室RPCs和气体电子倍增器GEMs）中的性能。此外，研究还试图将模拟的雪崩统计数据与现有的实验数据集进行验证，从而为优化气体增益谱提供一个可靠的预测框架。最终，这项工作提供了关于电场配置和气体组成如何相互作用的更深入的物理理解。

**电子传输和第一个Townsend系数**
在这些研究中，控制电荷倍增的基本过程是Townsend雪崩。当初级电子在气体介质中移动时，它们会经历多次弹性和非弹性碰撞。电子种群的增长由第一个Townsend系数α表征，α表示每个单位路径长度产生的离子对数量。这里，α对降低的电场（E/p）的依赖性是通过微 vervías.
**微观仿真框架**
我们使用Garfield++仿真工具包，这是一个专为气体探测器详细建模设计的面向对象框架（Veenhof, 1998; Schindler, 2012b）。与依赖平均传输系数的传统宏观代码不同，Garfield++可以对单个电子进行微观跟踪。这种方法对于研究气体增益波动至关重要，因为它考虑了每个电子的确切碰撞历史。该框架与HEED集成，用于初级...

**结果与讨论**
以下分析通过一系列特征性仿真，研究了控制电子雪崩统计的物理参数。评估重点关注决定增益波动的微观相互作用以及检测系统的整体稳定性，从而验证了微观跟踪方法与既定物理原理的一致性。图3展示了从电子-分子碰撞中系统不确定性的传播过程。

**结论**
本研究成功地将气体增益波动的分析从传统近似方法过渡到使用Garfield++和Magboltz工具包的高保真微观框架。通过逐事件蒙特卡罗跟踪电子传输，这些研究对控制气体探测器中电荷倍增的随机过程提供了明确的表征。本研究的一个重要成就是对仿真的严格验证。

**作者贡献声明**
Ijaz Ahmed：撰写——原始草稿、方法论、研究、概念化。
Muhammad Saqlain：可视化、软件、资源、方法论、研究。
Farzana Ahmad：撰写——原始草稿、方法论、研究、正式分析、概念化。
Jamil Muhammad：撰写——审稿与编辑、撰写——原始草稿、监督、研究、正式分析。

**未引用参考文献**
Biagi (1988), Buchtela (2003), Davydov (0000), ELGAS Ltd (2026), Fabjan (2020), Francais (2014), Francais (2017), Mandelli et al. (2020), Med Phys (0000), Pizzimento et al. (2022a), Schindler (2012a), Sipaj (2012), Veenhof and Renner (2010)