[論文レビュー] Bertini intra-nuclear cascade implementation in Geant4

A. Heikkinen, N.V. Stepanov+1 more|arXiv (Cornell University)|Jun 2, 2003

Nuclear physics research studies参考文献 16被引用数 68

ひとこと要約

本論文では、Geant4バージョン5.0におけるBertini核内カスケード（INC）モデルの実装を提示し、陽子、中性子、π中間子、光子、および軽い原子核のための中間エネルギー核反応（100 MeV–3 GeV）のシミュレーションを可能にしている。このモデルはBertiniのカスケード力学とエキシトン、準平衡、核分裂、爆発、蒸発モデルを統合しており、鉄および lead の標的における中性子およびπ中間子の二重微分断面積データと良好な一致を示している。

ABSTRACT

We present here a intra-nuclear cascade model implemented in Geant4 5.0. The cascade model is based on re-engineering of INUCL code. Models included are Bertini intra-nuclear cascade model with exitons, pre-equilibrium model, nucleus explosion model, fission model, and evaporation model. Intermediate energy nuclear reactions from 100 MeV to 3 GeV energy are treated for proton, neutron, pions, photon and nuclear isotopes. We represent overview of the models, review results achieved from simulations and make comparisons with experimental data.

研究の動機と目的

• 中間エネルギーのはね出し反応（100 MeV–3 GeV）における陽子、中性子、π中間子、光子、および同位体の核反応を正確にシミュレートできるGeant4シミュレーションフレームワーク内での包括的核内カスケードモデルの開発。
• Bertini INC、エキシトン、準平衡、核分裂、核爆発、蒸発モデルを統合した、Geant4互換の統一フレームワークへの統合。
• 鉄および lead における陽子誘発反応からの中性子およびπ中間子の二重微分断面積実験データとの照合によるモデル性能の検証。
• 将来的な調整および高エネルギー領域への拡張、物理的精度の向上のための基盤の構築。

提案手法

• 相対論的運動論と三領域核密度近似を用いて、INUCLコードに基づいてBertini核内カスケードモデルを再実装。
• 自由核子-核子断面積と領域依存核子密度を用いて、確率的輸送法により粒子相互作用をシミュレート。
• パウリの排他原理を各衝突ステップで適用し、エネルギー閾値（E > 2 MeV）によって二次粒子が脱出可能か、またはさらなる相互作用を引き起こすかを決定。
• カスケード中に動的に更新されるエキシトンモデルにより、核内の残余励起エネルギーの進化を記述。
• カスケード終了後のプロセス（核分裂、核爆発、準平衡放射、統計的蒸発）を、現象論的および統計的手法を用いて逐次的にモデル化。
• 最終状態の粒子放出はDostrovsky法を用いて計算し、励起エネルギーが0.1 MeV未満に低下するまで蒸発を継続。その後、10–15 MeV未満になるまでガンマ線放出を継続。

実験結果

リサーチクエスチョン

• RQ1 256 MeVおよび597 MeVの陽子ビームによる鉄および lead の標的における二重微分中性子放出スペクトルを、Bertini核内カスケードモデルがどの程度正確に再現できるか？
• RQ2 Pb標的における585 MeVの入射陽子エネルギー下でのπ⁺生成断面積を予測する際、実装されたモデルの性能はいかがなものか？
• RQ3 カスケード、準平衡、蒸発を統合したモデルが、100 MeVから3 GeVのエネルギー範囲で実験データをどの程度よく再現できるか？
• RQ4 核分裂および核爆発モデルの組み込みが、高励起状態にある残余核の記述をどの程度改善するか？
• RQ5 多様な核標的および入射粒子に対して、エネルギー保存則と物理的整合性が一貫して保持されるか？

主な発見

• Geant4に実装されたBertini INCモデルは、256 MeVの鉄標的における中性子の二重微分断面積実験データと、7.5°から150°の角度範囲で合理的な一致を示している。
• Pb標的における597 MeVのエネルギーで中性子生成量が、特に60°および120°の角度で良好な一致を示しており、高エネルギー域におけるモデル性能の妥当性が裏付けられている。
• Pb標的における585 MeVの入射陽子エネルギー下でのπ⁺生成断面積予測は、22.5°から135°の角度範囲で中程度から良好な一致を示しており、実験データと整合している。
• モデルは、準平衡および蒸発寄与を含む二次粒子放出のエネルギー的・角度的依存性をうまく再現している。
• カスケード終了後に残存励起エネルギーをチェックすることで、エネルギー保存則が維持されており、シミュレーションの内部整合性が確認された。
• Al、Fe、Pbなどの多様な標的および1–8 GeVの広いエネルギー範囲で性能が確認され、実行時間測定により標準的なコンピュータハードウェアでも実行可能であることが示された。