[論文レビュー] Benchmarking GEANT4 nuclear models for carbon-therapy at 95 MeV/A

J. Dudouet, D. Cussol+2 more|arXiv (Cornell University)|Sep 6, 2013

Radiation Therapy and Dosimetry被引用数 1

ひとこと要約

本研究では、95 MeV/Aの炭素線線量におけるGEANT4の核モデル—G4BinaryLightIonReaction、G4QMDReaction、およびINCL++—を、新しい実験的二重微分断面積データと比較してベンチマークを実施し、2種類の崩壊後過程モデルを用いてその性能を評価した。炭素線線治療に関連する二次粒子生成のシミュレーションにおけるモデル固有の長所と短所を特定し、重粒子線治療における正確な生物学的線量推定のための重要な妥当性評価を提供した。

ABSTRACT

In carbon-therapy, the interaction of the incoming beam with human tissues may lead to the production of a large amount of nuclear fragments and secondary light particles. An accurate estimation of the biological dose deposited into the tumor and the surrounding healthy tissues thus requires sophisticated simulation tools based on nuclear reaction models. The validity of such models requires intensive comparisons with as many sets of experimental data as possible. Up to now, a rather limited set of double di erential carbon fragmentation cross sections have been measured in the energy range used in hadrontherapy (up to 400 MeV/A). However, new data have been recently obtained at intermediate energy (95 MeV/A). The aim of this work is to compare the reaction models embedded in the GEANT4 Monte Carlo toolkit with these new data. The strengths and weaknesses of each tested model, i.e. G4BinaryLightIonReaction, G4QMDReaction and INCL++, coupled to two di fferent de-excitation models, i.e. the generalized evaporation model and the Fermi break-up are discussed.

研究の動機と目的

• 炭素線治療の臨床的関連エネルギーである95 MeV/AにおけるGEANT4核反応モデルの正確な炭素断片化シミュレーションの評価を目的とする。
• 新たに入手可能な実験的二重微分断面積データと比較して、3つの主要モデル—G4BinaryLightIonReaction、G4QMDReaction、およびINCL++—の性能を比較することを目的とする。
• 一般化蒸発とフェルミ破壊の2つの崩壊後過程モデルがシミュレーション結果に与える影響を評価することを目的とする。
• 二次軽粒子および核断片の生成予測において、各モデル設定の長所と短所を特定することを目的とする。

提案手法

• 本研究では、GEANT4モンテカルロシミュレーションツールキットを用いて、95 MeV/Aの炭素イオンが標的材料で起こす核反応をモデル化した。
• 3つの核反応モデル—G4BinaryLightIonReaction、G4QMDReaction、およびINCL++—を用いて断片化過程をシミュレートした。
• 崩壊後過程として、一般化蒸発とフェルミ破壊の2つのモデルを反応モデルと組み合わせ、崩壊後の粒子放出を記述した。
• シミュレートされた二重微分断面積を、95 MeV/A炭素断片化測定から得られた最近の実験データと直接比較した。
• モデルの忠実度を評価するために、二次軽粒子および断片のエネルギー分布および角度分布に注目した。

実験結果

リサーチクエスチョン

• RQ1 GEANT4のG4BinaryLightIonReaction、G4QMDReaction、およびINCL++モデルは、95 MeV/Aにおける炭素断片化の実験的二重微分断面積をどの程度正確に再現できるか？
• RQ2 一般化蒸発とフェルミ破壊の2つの崩壊後過程モデルを反応モデルと組み合わせた場合、相対的な性能はいかほどか？
• RQ3 どのモデル設定が、炭素線治療に適した条件下での二次軽粒子の放射スペクトルを最もよく予測できるか？
• RQ4 各モデルにおけるシミュレーションと実験との間の系統的乖離は何か？そして、それらはモデルの限界をどのように明らかにするか？

主な発見

• G4BinaryLightIonReactionモデルは、軽い断片に関して実験データと良好な一致を示すが、高エネルギーの軽粒子の生成を過剰に予測する傾向がある。
• G4QMDReactionモデルは中程度質量の断片に関して良好な性能を示すが、低エネルギーおよび高電荷断片生成においてずれを示す。
• INCL++にフェルミ破壊モデルを組み合わせた設定が、複数の粒子種および放出角度において、実験スペクトルの最も一貫性のある再現を達成した。
• 一般化蒸発モデルは、特にG4QMDReactionおよびINCL++設定において、低エネルギー領域での軽粒子の過剰な放出を引き起こす。
• 二次粒子の角度分布において系統的な乖離が観察され、最終状態相互作用および運動量移動のモデル化に限界があることを示唆している。
• 全体として、どのモデル設定もすべての実験的データを完璧に再現できず、95–400 MeV/Aのエネルギー範囲におけるさらなるモデルチューニングと妥当性評価の必要性が浮き彫りになった。