[論文レビュー] 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モンテカルロコードにおける高度な電離エネルギー損失モデルを提示する。このモデルは、有限核サイズおよび低・高エネルギー効果（Barkas効果、シェル補正など）を考慮したLindhard-Sørensen補正を組み込む。モデルは実験データおよびNISTの表値と非常に良好な一致を示し、特に1–2%の誤差で陽子に対して、10%以内で鉛やウランのような超重イオンに対して、keV/AからGeV/Aの広いエネルギー範囲で妥当である。

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 を得る。
• 低エネルギー補正にはBarkas効果（z³項）が、Thomas-Fermiに基づく関数F(V)を介して取り入れられ、遅い入射粒子に対してシェル補正（Δ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%以内の一致を示し、4 MeV以上でMCNP5の3%のずれを改善した。
• アルファ粒子に関しては、400 keV/A未満で10–15%のずれを示し、この領域における理論と実験の既知の不一致と整合的である。
• Lindhard-Sørensen補正により、BMA（Bloch-Mott-Ahlen）アプローチに比べて2–3%の系統的過小評価が軽減され、特にXeのような高Zイオンで顕著である。
• 鉛およびウランイオンに関しては、低エネルギー領域で10%以内の良好な一致を示し、超相対論的エネルギー領域における有限核サイズ効果による電離損失の飽和を正しく捉えている。
• 160 GeV/uで、アルゴン中の鉛イオンの実験データと一致し、点状モデルが示す対数的上昇が有限核サイズ効果によって抑制されることを確認した。
• モデルは電離損失の密度依存性を成功裏に再現し、最低密度の標的に、超相対論的ウランイオンで最大のエネルギー損失を示した。