[论文解读] Modeling heavy ion ionization loss in the MARS15 code

Igor Rakhno, N. Mokhov+1 more|ArXiv.org|May 5, 2005

Particle Detector Development and Performance参考文献 14被引用 3

一句话总结

本文在MARS15蒙特卡洛代码中提出了一种先进的电离能损模型，引入了Lindhard-Sørensen修正，以考虑有限核尺寸及低/高能效应（如Barkas效应和壳层修正）。该模型与实验数据及NIST表列值高度一致，尤其在质子能量范围内（1–2%以内），以及在从keV/A到GeV/A的宽广能量范围内，对超重离子（如铅和铀）的计算结果误差最高可达10%。

ABSTRACT

The needs of various accelerator and space projects stimulated recent developments to the MARS Monte Carlo code. One of the essential parts of those is heavy ion ionization energy loss. This paper describes an implementation of several corrections to dE/dx in order to take into account the deviations from the Bethe theory at low and high energies as well as the effect of a finite nuclear size at ultra-relativistic energies. Special attention is paid to the transition energy region where the onset of the effect of a finite nuclear size is observed. Comparisons with experimental data and NIST data are presented.

研究动机与目标

• 改进MARS15蒙特卡洛代码中重离子电离能损计算的精度。
• 解决低能和高能区域（尤其是相对论性及超重离子）中与Bethe理论的偏差。
• 在极端相对论区域引入有限核尺寸效应，取代旧有的Mott和Bloch修正。
• 在宽广的能量范围内，通过实验数据和NIST表列值验证该模型。
• 支持加速器、空间和屏蔽应用中涉及重离子的辐射输运模拟。

提出的方法

• 该模型采用三区法：能量高于10 MeV/A时使用Bethe理论，低于1 MeV/A时使用表列质子停止功率，中间区域采用混合插值方法。
• Lindhard-Sørensen形式化方法用于相对论性修正，通过相对论性库仑相移和部分波求和来考虑有限核尺寸效应。
• 当γmₑcR ≫ ħ/2时，采用渐近超相对论极限，得到L_ultra = ln(2c/Rωₚ) - 0.2。
• 低能修正包括通过基于Thomas-Fermi的函数F(V)实现的Barkas效应（z³项），以及对慢速入射粒子的壳层修正（ΔL_shell = -C/Z）。
• 将这些修正整合为统一的电离对数表达式L(β) = L₀(β) + δ/2 + ΣΔLᵢ，其中有效电荷效应根据靶原子序数Z进行调整。
• 采用混合插值程序，确保低能、中间能和高能区域之间的平滑过渡。

实验结果

研究问题

• RQ1 如何在keV/A至GeV/A的全能量范围内准确建模重离子的电离能损？
• RQ2 在超相对论能量下，有限核尺寸对电离能损有何影响？
• RQ3 Barkas效应和壳层修正如何提升低能区与实验数据的一致性？
• RQ4 Lindhard-Sørensen形式化方法与传统Mott和Bloch修正在相对论区域相比有何差异？
• RQ5 MARS15模型在多大程度上复现了NIST表列的质子和α粒子停止功率？

主要发现

• 在全能量范围内，MARS15模型与NIST数据的偏差控制在1.3%以内，优于MCNP5在4 MeV以上时3%的偏差。
• 对于α粒子，模型在400 keV/A以下存在10–15%的偏差，与该能量区域已知的理论-实验偏差一致。
• 与BMA（Bloch-Mott-Ahlen）方法相比，Lindhard-Sørensen修正将电离能损的系统性低估降低了2–3%，尤其在高Z离子（如Xe）中表现显著。
• 对于铅和铀离子，模型在低能区与实验数据一致，偏差在10%以内，并正确捕捉了在超相对论能量下由有限核尺寸引起的电离能损饱和效应。
• 在160 GeV/u能量下，模型与氩气中铅离子的实验数据吻合，证实有限尺寸效应抑制了点状模型中出现的对数上升趋势。
• 模型成功再现了电离能损的密度依赖性，表现为在超相对论能量下，超重铀离子在最低密度靶材中表现出最高的能量损失。