[论文解读] MARS Code Developments

N. Mokhov, S. Striganov+4 more|ArXiv.org|Dec 15, 1998

Particle physics theoretical and experimental studies参考文献 14被引用 32

一句话总结

该论文针对复杂介质中强子与电磁级联过程的模拟，显著提升了MARS蒙特卡罗代码的性能，包括引入新的核截面库、改进了π介子产生与μ子吸收的模型，以及优化了电磁相互作用的处理方式。更新后的mars13(98)代码在1 MeV至100 TeV的宽广能量范围内与实验数据高度一致，验证了其在加速器屏蔽与探测器设计应用中的可靠性。

ABSTRACT

Recent developments in the physical model of 1 MeV to 100 TeV hadron and lepton interactions with nuclei and atoms are described. These include a new nuclear cross section library, a model for soft pion production, the cascade-exciton model, the dual parton model, deuteron-nucleus and neutrino-nucleus interaction models, detailed description of muon, pion, and antiproton absorption and a unified treatment of muon and charged hadron electromagnetic interactions with matter. New algorithms are implemented into the MARS13(98)} Monte Carlo code and benchmarked against experimental data. The code capabilities to simulate cascades and generate a variety of results in complex media have been also enhanced.

研究动机与目标

• 在宽广能量范围（1 MeV–100 TeV）内提升强子与轻子与物质相互作用的物理建模精度。
• 利用新编纂的数据与理论框架，提升强子-核与光子-核相互作用的核截面精度。
• 开发并实现高能强子相互作用中软π介子产生、级联-激发模型及双部分子动力学的先进模型。
• 统一并优化μ子与带电强子电磁相互作用的处理方式，包括能量损失与轫致辐射。
• 扩展代码在评分、几何处理及与MCNP和ANSYS等其他模拟工具集成方面的能力，以支持复杂工程应用。

提出的方法

• 采用基于级联-激发模型（CEM）改进的参数化方法与新编纂的数据，构建了适用于强子（p, n, π±, K±, p̄）在1 MeV至100 TeV能量范围内的新核截面库。
• 应用Glauber多重散射理论，结合核子密度分布（费米型分布），计算高能区强子-核的总截面、非弹性截面与弹性截面。
• 提出一种新的现象学模型用于pA相互作用中的π介子产生，采用依赖于横动量弱、束流与π介子动量中等的比值函数R^{pA→π±X}。
• 用简化算法替代复杂的μ子轫致辐射形式化处理，该算法在总截面与dE/dx上与精确公式一致，误差小于1%。
• 通过Bhabha公式实现δ电子能量采样与受限能量损失计算，提升电磁相互作用建模精度。
• 集成自动MCNP接口以实现低能中子与光子输运的精确模拟，并引入面向对象的几何建模与基于网格的直方图评分方法，提升评分灵活性。

实验结果

研究问题

• RQ1 如何在1 MeV–100 TeV能量范围内提升强子-核总截面与微分截面的计算精度？
• RQ2 在pA相互作用中，新的现象学π介子产生模型在前向与后向半球是否能较好地再现实验数据？
• RQ3 更新后的μ子与带电强子电磁相互作用模型与实验测量的能量损失及角分布数据相比表现如何？
• RQ4 MARS13(98)代码能否可靠模拟复杂非均质介质中的级联过程，且具备改进的评分与材料组成处理能力？
• RQ5 在中子衰减与剂量计算等基准测试场景中，MARS13(98)与GEANT、MCNP、FLUKA等成熟代码相比性能表现如何？

主要发现

• pA相互作用中π介子产生的新模型在全动量范围内与实验数据吻合良好，尤其在π介子动量为0.5–5 GeV/c的困难区域表现优异。
• MARS13(98)对8 GeV/c质子束轰击厚铅靶产生的π介子能谱模拟结果与实验数据基本一致，在低动量区域优于GEANT。
• 使用MARS13(98)计算的经过锂吸收体（50 cm与150 cm）后的μ子能量谱与角分布，与GEANT和ICool的结果高度一致，验证了新电磁相互作用模型的可靠性。
• 简化后的μ子轫致辐射算法在总截面与dE/dx上与精确公式一致，误差小于1%，在保持高精度的同时显著降低计算成本。
• 与SATIF-4数据的基准测试表明，MARS13(98)在混凝土中子衰减（100–400 MeV）与组织等效体模中吸收剂量（0.1–10 GeV）方面，与HETC、ANISN、MCNP等代码结果高度一致。
• 与MCNP的自动接口实现了低能中子与光子输运的精确模拟，新材料处理方法无需材料均质化，显著提升了模拟保真度。