[论文解读] Monte Carlo Simulations of Secondary Cosmic-Ray Variations in Atmospheric Electric Fields : Implications for Long Duration Electron and Gamma-ray Emissions from Thunderclouds

Harufumi Tsuchiya|arXiv (Cornell University)|Feb 26, 2026

Lightning and Electromagnetic Phenomena被引用 0

一句话总结

该论文使用 PHITS 蒙特卡洛模拟和真实的二级宇宙射线谱来研究雷暴云电场如何改变二级宇宙射线，并识别在产生长持续电子和伽马射线辐射（包括高于 40 MeV 的伽马射线）时的条件。

ABSTRACT

Monte Carlo simulations were conducted using the Particle and Heavy Ion Transport code System (PHITS) to investigate the role of secondary cosmic rays in the generation of long-duration bursts from thunderclouds and to clarify the conditions of the electric field region responsible for particle acceleration. The simulations utilized realistic secondary cosmic-ray spectra, including gamma rays, electrons, positrons, and muons, as input. The simulation results indicate that gamma rays provide the dominant supply of seed electrons for long-duration bursts, regardless of the geometry or strength of the electric field region. They also reveal the structure and strength of the electric field region required to produce gamma rays exceeding several tens of MeV, which have so far been detected only by high-altitude observations. Furthermore, the fluxes of long-duration bursts estimated from the simulation results were compared with observational data to constrain the properties of the electric field region. In particular, the comparison with measurements at Yangbajing, located at an altitude of 4.3~km, helps narrow down the possible range of electric field strengths and configurations.

研究动机与目标

• 阐明并量化稳态二级宇宙射线与局部雷暴云电场（EFs）相互作用的机理与影响。
• 确定在 EF 内最有效提供相对论性雪崩种子电子的二级宇宙射线分量。
• 表征产生超过几十 MeV 的伽马射线的 EF 区几何形状与强度，并与高空观测结果进行比较。
• 利用对杨巴景观观测（Yangbajing）测量的长期 bursts 的模型化通量来约束电场属性。

提出的方法

• 使用 PHITS（ver.3.34）模拟 gamma 射线、电子、正电子和μ子在建模的 EF 区域中的传播与相互作用。
• 基于 EXPACS 输出在海拔 4.3 公里处注入真实的二级宇宙射线谱作为初始输入。
• 将 EF 区域建模为具有可调宽度 W、长度 L、云基高度 H 的矩形立方体，并将电子向地面加速。
• 考虑包括制动辐射、离子化损失、散射、光电效应、康普顿散射和对对产生等在内的传播与倍增过程。
• 分析逃离 EF 区域并到达地面或恰在 EF 区底部之下的电子与伽马射线的能量谱与空间分布。

Figure 1: Secondary cosmic-ray spectra expected at Yangbajing (4.3 km a.s.l.), calculated by EXPACS Sato ( 2015 )

实验结果

研究问题

• RQ1 二级宇宙射线的哪些分量（伽马射线、电子、正电子、μ子）在雷暴云 EF 区域内最有效提供高能电子？
• RQ2 EF 强度与 EF 区几何形状（W、L、H）如何影响长期 bursts 的产生、光谱与地面通量？
• RQ3 哪些 EF 配置能够产生超过 40 MeV 的伽马射线，模型结果与高海拔观测结果如何比较？
• RQ4 伽玛射线分量如何促进和推动雷暴云中的相对论性逃逸电子雪崩（RREA）的起始与发展？

主要发现

• 伽马射线在各 EF 几何与强度下提供了主导的种子电子，用于长期 bursts。
• 当 EF 区域长度较长时（L 较大），在 EF 强度超过 RREA 阈值时会增强电子与伽马射线通量；若 EF 未达到阈值，则效应减弱。
• 增大 W（宽度）通常会因横向散射和制动辐射光子的再相互作用而提高电子与伽马射线通量。
• 电子对累积电子光谱贡献显著，在各种条件下伽马射线贡献大约占约 60%。
• 在强 EF 下的伽马射线光谱类似于 RREA-形式光谱，较长的 L 时电子平均能量约为 9–12 MeV；而非常短的 L 可能产生更高特征光子能量（可达 ~40 MeV）。
• 通过将对杨巴景观观测（4.3 公里 a.s.l.）的模型化通量增强与观测进行比较，本研究对 EF 强度与配置范围进行了约束，表明地面时段 bursts 对 EF 强度与几何高度敏感。

Figure 2: Assumed structure and orientation of the electric field. The assumed EF region is shown in light blue. Dashed arrows indicate the parameters $W$ , $L$ , and $H$ used in the simulations. The EF orientation is represented by a solid arrow. The origin of the simulation coordinate system is de