A Fano cavity test for Monte Carlo proton transport algorithms.

PURPOSE In the scope of reference dosimetry of radiotherapy beams, Monte Carlo (MC) simulations are widely used to compute ionization chamber dose response accurately. Uncertainties related to the transport algorithm can be verified performing self-consistency tests, i.e., the so-called "Fano cavity test." The Fano cavity test is based on the Fano theorem, which states that under charged particle equilibrium conditions, the charged particle fluence is independent of the mass density of the media as long as the cross-sections are uniform. Such tests have not been performed yet for MC codes simulating proton transport. The objectives of this study are to design a new Fano cavity test for proton MC and to implement the methodology in two MC codes: Geant4 and PENELOPE extended to protons (PENH). METHODS The new Fano test is designed to evaluate the accuracy of proton transport. Virtual particles with an energy of E0 and a mass macroscopic cross section of Σρ are transported, having the ability to generate protons with kinetic energy E0 and to be restored after each interaction, thus providing proton equilibrium. To perform the test, the authors use a simplified simulation model and rigorously demonstrate that the computed cavity dose per incident fluence must equal ΣE0ρ, as expected in classic Fano tests. The implementation of the test is performed in Geant4 and PENH. The geometry used for testing is a 10 × 10 cm(2) parallel virtual field and a cavity (2 × 2 × 0.2 cm(3) size) in a water phantom with dimensions large enough to ensure proton equilibrium. RESULTS For conservative user-defined simulation parameters (leading to small step sizes), both Geant4 and PENH pass the Fano cavity test within 0.1%. However, differences of 0.6% and 0.7% were observed for PENH and Geant4, respectively, using larger step sizes. For PENH, the difference is attributed to the random-hinge method that introduces an artificial energy straggling if step size is not small enough. CONCLUSIONS Using conservative user-defined simulation parameters, both PENH and Geant4 pass the Fano cavity test for proton transport. Our methodology is applicable to any kind of charged particle, provided that the considered MC code is able to track the charged particle considered.

[1]  J. Seuntjens,et al.  Consistency test of the electron transport algorithm in the GEANT4 Monte Carlo code , 2005, Physics in medicine and biology.

[2]  M. S. Yeom,et al.  Monte Carlo calculation of the ionization chamber response to 60Co beam using PENELOPE. , 2006, Medical physics.

[3]  F. Salvat A generic algorithm for Monte Carlo simulation of proton transport , 2013 .

[4]  D. Rogers Fifty years of Monte Carlo simulations for medical physics , 2006, Physics in medicine and biology.

[5]  Iwan Kawrakow,et al.  On the representation of electron multiple elastic-scattering distributions for Monte Carlo calculations , 1998 .

[6]  Iwan Kawrakow,et al.  Ionization chamber gradient effects in nonstandard beam configurations. , 2009, Medical physics.

[7]  I Kawrakow,et al.  Investigation of systematic uncertainties in Monte Carlo-calculated beam quality correction factors , 2010, Physics in medicine and biology.

[8]  V. Smyth Interface effects in the Monte Carlo simulation of electron tracks. , 1986, Medical physics.

[9]  A. Bielajew An analytic theory of the point-source nonuniformity correction factor for thick-walled ionisation chambers in photon beams , 1990 .

[10]  S Stute,et al.  GATE V6: a major enhancement of the GATE simulation platform enabling modelling of CT and radiotherapy , 2011, Physics in medicine and biology.

[11]  E. Sterpin,et al.  Monte Carlo computed machine-specific correction factors for reference dosimetry of TomoTherapy static beam for several ion chambers. , 2012, Medical physics.

[12]  Alex F. Bielajew,et al.  Correction factors for thick-walled ionisation chambers in point-source photon beams , 1990 .

[13]  S. Elles,et al.  Geant4 and Fano cavity test: where are we? , 2008 .

[14]  I. Kawrakow Accurate condensed history Monte Carlo simulation of electron transport. I. EGSnrc, the new EGS4 version. , 2000, Medical physics.

[15]  E Sterpin,et al.  Extension of PENELOPE to protons: simulation of nuclear reactions and benchmark with Geant4. , 2013, Medical physics.

[16]  S Agosteo,et al.  Secondary neutron and photon dose in proton therapy. , 1998, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[17]  U. Fano,et al.  Note on the Bragg-Gray cavity principle for measuring energy dissipation. , 1954, Radiation research.

[18]  H. Bouchard Corrigendum: A theoretical re-examination of Spencer–Attix cavity theory , 2012, Physics in medicine and biology.

[19]  M. Fippel,et al.  A Monte Carlo dose calculation algorithm for proton therapy. , 2004, Medical physics.

[20]  Josep Sempau,et al.  Electron beam quality correction factors for plane-parallel ionization chambers: Monte Carlo calculations using the PENELOPE system. , 2004, Physics in medicine and biology.

[21]  J. Sempau,et al.  PENELOPE-2006: A Code System for Monte Carlo Simulation of Electron and Photon Transport , 2009 .

[22]  D. Rogers,et al.  Monte Carlo calculations of kQ, the beam quality conversion factor. , 2010, Medical physics.

[23]  I. Kawrakow Accurate condensed history Monte Carlo simulation of electron transport. II. Application to ion chamber response simulations. , 2000, Medical physics.

[24]  Steve B. Jiang,et al.  GPU-based fast Monte Carlo dose calculation for proton therapy , 2012, Physics in medicine and biology.

[25]  P Andreo,et al.  Configuration of the electron transport algorithm of PENELOPE to simulate ion chambers , 2006, Physics in medicine and biology.

[26]  J. Seuntjens,et al.  On charged particle equilibrium violation in external photon fields. , 2012, Medical physics.

[27]  F. Verhaegen,et al.  A systematic Monte Carlo study of secondary electron fluence perturbation in clinical proton beams (70-250 MeV) for cylindrical and spherical ion chambers. , 2001, Medical physics.

[28]  I. Kawrakow,et al.  Technical note: overprediction of dose with default PRESTA-I boundary crossing in DOSXYZnrc and BEAMnrc. , 2007, Medical physics.

[29]  H Paganetti,et al.  Nuclear interactions in proton therapy: dose and relative biological effect distributions originating from primary and secondary particles. , 2002, Physics in medicine and biology.

[30]  M R McEwen,et al.  Measured and Monte Carlo calculated kQ factors: Accuracy and comparison. , 2011, Medical physics.