[論文レビュー] Nuclear halo of a 177 MeV proton beam in water: theory, measurement and parameterization

B. Gottschalk, Ethan W. Cascio+2 more|arXiv (Cornell University)|Sep 5, 2014

Radiation Therapy and Dosimetry参考文献 37被引用数 9

ひとこと要約

本研究は、水中における177 MeVプロトンペンチルビームの核殻を測定・モデル化し、モデル依存およびモデル非依存のフィット手法を用いてコア、ハロー、オーラの各線量成分に分離する。ハロー半径はビームの進行距離のおよそ1/3であることが判明し、積分深さ線量（T(w)）をパrameterizationに用いることは疑問視され、ブラッグピークチャネルが正確なペンチルビームモデル化に不要であることが示された。モデル非依存のフィット手法は、9%の全般的なRMS誤差を示し、より優れた移植性と正確性を示した。

ABSTRACT

The dose distribution of a monoenergetic pencil beam in water consists of an electromagnetic "core", a "halo" from charged nuclear secondaries, and a much larger "aura" from neutral secondaries. These regions overlap, but each has distinct spatial characteristics. We have measured the core/halo using a 177MeV test beam offset in a water tank. The beam monitor was a fluence calibrated plane parallel ionization chamber (IC) and the field chamber, a dose calibrated Exradin T1, so the dose measurements are absolute (MeV/g/p). We performed depth-dose scans at ten displacements from the beam axis ranging from 0 to 10cm. The dose spans five orders of magnitude, and the transition from halo to aura is clearly visible. We have performed model-dependent (MD) and model-independent (MI) fits to the data. The MD fit separates the dose into core, elastic/inelastic nuclear, nonelastic nuclear and aura terms, and achieves a global rms measurement/fit ratio of 15%. The MI fit uses cubic splines and the same ratio is 9%. We review the literature, in particular the use of Pedroni's parametrization of the core/halo. Several papers improve on his Gaussian transverse distribution of the halo, but all retain his T(w), the radial integral of the depth-dose multiplying both the core and halo terms and motivating measurements with large "Bragg peak chambers" (BPCs). We argue that this use of T(w), which by its definition includes energy deposition by nuclear secondaries, is incorrect. T(w) should be replaced in the core term, and in at least part of the halo, by a purely electromagnetic mass stopping power. BPC measurements are unnecessary, and irrelevant to parameterizing the pencil beam.

研究の動機と目的

• 177 MeVプロトンペンチルビームの水中における空間線量分布を測定・モデル化し、コア、ハロー、オーラの各成分を区別する。
• コアおよびハローのパrameterizationに一般的に用いられるT(w)、すなわち深さ線量の半径方向積分の使用を疑問視し、核二次粒子によるエネルギー損失が誤って含まれていると主張する。
• 25パラメータを用いたモデル依存フィットと比較して、立方スプラインを用いたモデル非依存フィットが、コア／ハロー線量分布の特徴づけにおいてより優れた移植性と正確性を示すことを示す。
• ブラッグピークチャネルが、特に小領域線量測定において、概念的に誤りであり、正確なペンチルビーム線量測定に不要であることを示す。
• モンテカルロシミュレーションのベンチマークを提供するため、ハロー領域における電磁的、弾性的、非弾性的核反応プロセスの寄与を正確に定量する。

提案手法

• 10個の半径方向オフセット（0–10 cm）で深さ線量スキャンを実施。線量測定にはフィルランス補正済みのイオン化線量計と、線量補正済みのExradin T1線量計（MeV/g/p）を用いた。
• 25パラメータを用いたモデル依存（MD）フィットを実施し、電磁的（コア）、弾性核反応、非弾性核反応、オーラ成分に線量を分解する。
• 深さおよび半径方向に立方スプラインを用いたモデル非依存（MI）フィットを実施し、二次微分密度に基づく最適化された初期スプライン点を用いた双線形補間を適用した。
• 運動学的解析を用いて最大ハロー半径を推定し、運動量保存則および重心系散乱角に基づき、ビーム進行距離のおよそ1/3であることが示された。
• コアおよびハローのパラメータ化においてT(w)を、核二次粒子の寄与を含まない純粋な電磁的質量停止力S_em/ρ_waterに置き換えることで、核二次粒子による汚染を排除した。
• 文献データ、モンテカルロシミュレーション（Gate/Geant4）、既存のパラメータ化（Pedroni, Zhang, Zhu）との比較を通じて結果を妥当性確認し、ブラッグピークチャネル測定値とも比較した。

実験結果

リサーチクエスチョン

• RQ1 177 MeVプロトンビームの水中における核ハローの空間的拡がりと半径依存性は何か？また、ビーム進行距離とどのように関係しているか？
• RQ2 深さ線量の半径方向積分T(w)を用いたコアおよびハローのモデリングがなぜ問題であるのか？より物理的に正確な代替手法は何か？
• RQ3 コア／ハロー線量分布の特徴づけにおいて、モデル依存フィットとモデル非依存フィットは、正確性と移植性の観点でどのように比較されるか？
• RQ4 ブラッグピークチャネルは、ペンチルビーム線量測定においてどれほど誤差や冗長性を生じるのか？その使用は排除可能か？
• RQ5 電磁的、弾性的、非弾性的核プロセスが、ハロー領域の総線量に果たす定量的寄与は何か？

主な発見

• 最大ハロー半径は、酸素および水素における弾性および準弾性散乱の運動学的制約に一致し、ビーム進行距離のおよそ1/3である。
• モデル非依存フィットは、測定値とフィット値の間で9%の全般的なRMS比を達成し、モデル依存フィット（15% RMS）を上回る性能を示した。また、概念的に単純で、より移植性に優れている。
• コアおよびハローのパラメータ化にT(w)を用いることは誤りである。核二次粒子によるエネルギー損失が含まれるためであり、S_em/ρ_water（純粋な電磁的停止力）に置き換えることで、物理的正確性が向上する。
• ブラッグピークチャネルは、ペンチルビームパラメータ化に不要である。その測定値は冗長であり、非コア成分を含むため、誤解を招く可能性がある。
• ハローは特に大半径領域で総線量に顕著な寄与を示し、フィールドサイズ依存の誤差を避けるために、小領域線量測定において考慮すべきである。
• ビームサイズを補正したPedroniの元々のパラメータ化は、コア領域で良好な一致を示したが、ハロー領域ではほとんど変化せず、T(w)に基づくモデルの限界を浮き彫りにした。