The FLUKA Code: An Accurate Simulation Tool for Particle Therapy

Monte Carlo (MC) codes are increasingly spreading in the hadrontherapy community due to their detailed description of radiation transport and interaction with matter. The suitability of a MC code for application to hadrontherapy demands accurate and reliable physical models capable of handling all components of the expected radiation field. This becomes extremely important for correctly performing not only physical but also biologically based dose calculations, especially in cases where ions heavier than protons are involved. In addition, accurate prediction of emerging secondary radiation is of utmost importance in innovative areas of research aiming at in vivo treatment verification. This contribution will address the recent developments of the FLUKA MC code and its practical applications in this field. Refinements of the FLUKA nuclear models in the therapeutic energy interval lead to an improved description of the mixed radiation field as shown in the presented benchmarks against experimental data with both 4He and 12C ion beams. Accurate description of ionization energy losses and of particle scattering and interactions lead to the excellent agreement of calculated depth–dose profiles with those measured at leading European hadron therapy centers, both with proton and ion beams. In order to support the application of FLUKA in hospital-based environments, Flair, the FLUKA graphical interface, has been enhanced with the capability of translating CT DICOM images into voxel-based computational phantoms in a fast and well-structured way. The interface is capable of importing also radiotherapy treatment data described in DICOM RT standard. In addition, the interface is equipped with an intuitive PET scanner geometry generator and automatic recording of coincidence events. Clinically, similar cases will be presented both in terms of absorbed dose and biological dose calculations describing the various available features.

[1]  R. Orecchia,et al.  Analysis of pelvis and abdomenal conformal radiotherapy using TCP and NTCP models , 2001 .

[2]  A. J. Cole,et al.  Projectile fragments from 86 MeV/nucleon 12C induced reactions , 1981 .

[3]  G. Molière,et al.  Theorie der Streuung schneller geladener Teilchen , 1955 .

[4]  A Ottolenghi,et al.  The FLUKA code: new developments and application to 1 GeV/n iron beams. , 2005, Advances in space research : the official journal of the Committee on Space Research.

[5]  Julia Bauer,et al.  Clinical implementation and range evaluation of in vivo PET dosimetry for particle irradiation in patients with primary glioma. , 2015, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[6]  Joao Seco,et al.  Simulation of prompt gamma-ray emission during proton radiotherapy , 2012, Physics in medicine and biology.

[7]  Gert Moliere,et al.  Theorie der Streuung schneller geladener Teilchen II Mehrfach-und Vielfachstreuung , 1948 .

[8]  K Parodi,et al.  Monte Carlo simulations to support start-up and treatment planning of scanned proton and carbon ion therapy at a synchrotron-based facility , 2012, Physics in medicine and biology.

[9]  Jong-Won Kim,et al.  Development of dosimetry tools for proton therapy research , 2010 .

[10]  Wousik Kim,et al.  Electron Nonionizing Energy Loss for Device Applications , 2009, IEEE Transactions on Nuclear Science.

[11]  D. Visvikis,et al.  Evaluation of attenuation correction methodology in the allegro PET system , 2004, IEEE Transactions on Nuclear Science.

[12]  K Parodi,et al.  Dose quantification from in-beam positron emission tomography. , 2004, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[13]  Katia Parodi,et al.  Implementation and workflow for PET monitoring of therapeutic ion irradiation: a comparison of in-beam, in-room, and off-line techniques , 2011, Physics in medicine and biology.

[14]  H. Ryde,et al.  The Physics of Heavy Ion Collisions at 100 MeV Per Nucleon , 1983 .

[15]  A. Ferrari,et al.  Electron-Photon Transport in FLUKA: Status , 2001 .

[16]  Katia Parodi,et al.  Clinical implementation of full Monte Carlo dose calculation in proton beam therapy , 2008, Physics in medicine and biology.

[17]  R. Leahy,et al.  Digimouse: a 3D whole body mouse atlas from CT and cryosection data , 2007, Physics in medicine and biology.

[18]  M. Kendall,et al.  Kendall's advanced theory of statistics , 1995 .

[19]  A Ottolenghi,et al.  Estimating mixed field effects: an application supporting the lack of a non-linear component for chromosome aberration induction by neutrons. , 2003, Radiation protection dosimetry.

[20]  M. Pullia,et al.  A new PET prototype for proton therapy: comparison of data and Monte Carlo simulations , 2013 .

[21]  Vlachoudis FLAIR: A POWERFUL BUT USER FRIENDLY GRAPHICAL INTERFACE FOR FLUKA , 2009 .

[22]  Sorge Flavor production in Pb(160A GeV) on Pb collisions: Effect of color ropes and hadronic rescattering. , 1995, Physical review. C, Nuclear physics.

[23]  S. D. Groot,et al.  SOME EXPLICIT FORMULAE FOR THE ANGULAR DISTRIBUTION AND POLARIZATION OF GAMMA RADIATION EMITTED BY ORIENTED NUCLEI , 1953 .

[24]  F. Bloch,et al.  Zur Bremsung rasch bewegter Teilchen beim Durchgang durch Materie , 1933 .

[25]  A. Mairani,et al.  LOW ENERGY NUCLEUS–NUCLEUS REACTIONS: THE BME APPROACH AND ITS INTERFACE WITH FLUKA , 2006 .

[26]  D Robertson,et al.  Optimizing a three-stage Compton camera for measuring prompt gamma rays emitted during proton radiotherapy , 2010, Physics in medicine and biology.

[27]  Klaus Wienhard,et al.  The ECAT HRRT: performance and first clinical application of the new high resolution research tomograph , 2000 .

[28]  E. A. Uehling Penetration of Heavy Charged Particles in Matter , 1954 .

[29]  Reuven Ramaty,et al.  Nuclear Deexcitation Gamma-Ray Lines from Accelerated Particle Interactions , 2001 .

[30]  A. Ferrari,et al.  An improved multiple scattering model for charged particle transport , 1992 .

[31]  C. Case,et al.  MOLIERE'S THEORY OF MULTIPLE SCATTERING. , 1968 .

[32]  M Zaider,et al.  The synergistic effects of different radiations. , 1980, Radiation research.

[33]  Chan Hyeong Kim,et al.  Two-Dimensional Prompt Gamma Measurement Simulation for In Vivo Dose Verification in Proton Therapy: A Monte Carlo Study , 2011 .

[34]  Explore Configuring,et al.  A Simulation Study to , 2004 .

[35]  Byung-Hwi Kang,et al.  Monte Carlo Design Study of a Gamma Detector System to Locate Distal Dose Falloff in Proton Therapy , 2009, IEEE Transactions on Nuclear Science.

[36]  Byung-Hwi Kang,et al.  Monte Carlo Design Study of a Gamma Detector System to Locate Distal Dose Falloff in Proton Therapy , 2009 .

[37]  Lev Davidovich Landau,et al.  On the energy loss of fast particles by ionization , 1944 .

[38]  M Scholz,et al.  The FLUKA Monte Carlo code coupled with the local effect model for biological calculations in carbon ion therapy , 2010, Physics in medicine and biology.

[39]  Katia Parodi,et al.  Accuracy of proton beam range verification using post-treatment positron emission tomography/computed tomography as function of treatment site. , 2011, International journal of radiation oncology, biology, physics.

[40]  R. N. Goble,et al.  Performance evaluation of the microPET P4: a PET system dedicated to animal imaging. , 2001, Physics in medicine and biology.

[41]  D Mason,et al.  SU-E-T-33: Pydicom: An Open Source DICOM Library , 2011 .

[42]  S. M. Seltzer,et al.  Density effect for the ionization loss of charged particles in various substances , 1982 .

[43]  Marc C. Huisman,et al.  Performance evaluation of the Philips MOSAIC small animal PET scanner , 2007, European Journal of Nuclear Medicine and Molecular Imaging.

[44]  Johannes Ranft,et al.  The Monte Carlo Event Generator DPMJET-III , 2000 .

[45]  G Brix,et al.  Performance evaluation of a whole-body PET scanner using the NEMA protocol. National Electrical Manufacturers Association. , 1997, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[46]  R. Bimbot,et al.  Semi-empirical formulae for heavy ion stopping powers in solids in the intermediate energy range , 1989 .

[47]  Katia Parodi,et al.  Monte Carlo-based parametrization of the lateral dose spread for clinical treatment planning of scanned proton and carbon ion beams , 2013, Journal of radiation research.

[48]  Jean Mougey,et al.  Experiments with 1.03 GeV 12C ions , 1982 .

[49]  K Parodi,et al.  In-beam PET monitoring of mono-energetic 16O and 12C beams: experiments and FLUKA simulations for homogeneous targets , 2009, Physics in medicine and biology.

[50]  R. H. Ritchie,et al.  Z 1 3 Effect in the Stopping Power of Matter for Charged Particles , 1972 .

[51]  A. Ferrari,et al.  Physical and biophysical characteristics of a fully modulated 72 MeV therapeutic proton beam: model predictions and experimental data , 1999 .

[52]  I Buvat,et al.  Distributions of secondary particles in proton and carbon-ion therapy: a comparison between GATE/Geant4 and FLUKA Monte Carlo codes , 2013, Physics in medicine and biology.

[53]  Roberta Matheoud,et al.  Performance characteristics obtained for a new 3-dimensional lutetium oxyorthosilicate-based whole-body PET/CT scanner with the National Electrical Manufacturers Association NU 2-2001 standard. , 2005, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[54]  Hou Qing,et al.  Analytic fitting to the mott cross section of electrons , 1995 .

[55]  Yoshitomo Uwamino,et al.  Measurements of double differential neutron production cross sections by 135 MeV/nucleon He, C, Ne and 95 MeV/nucleon Ar ions , 2001 .

[56]  V. Vlachoudis,et al.  The FLUKA Code: Developments and Challenges for High Energy and Medical Applications , 2014 .

[57]  Jong-Won Kim,et al.  Pinhole Camera Measurements of Prompt Gamma-rays for Detection of Beam Range Variation in Proton Therapy , 2009 .

[58]  D Dauvergne,et al.  Design Guidelines for a Double Scattering Compton Camera for Prompt-$\gamma$ Imaging During Ion Beam Therapy: A Monte Carlo Simulation Study , 2011, IEEE Transactions on Nuclear Science.

[59]  Wolfgang Enghardt,et al.  On the effectiveness of ion range determination from in-beam PET data , 2010, Physics in medicine and biology.

[60]  K Parodi,et al.  Integration and evaluation of automated Monte Carlo simulations in the clinical practice of scanned proton and carbon ion beam therapy , 2014, Physics in medicine and biology.

[61]  A Todd-Pokropek,et al.  A file format for the exchange of nuclear medicine image data: a specification of Interfile version 3.3. , 1992, Nuclear medicine communications.

[62]  F Verhaegen,et al.  Monte Carlo calculations of positron emitter yields in proton radiotherapy , 2012, Physics in medicine and biology.

[63]  P. R. Sala,et al.  Cascade particles, nuclear evaporation, and residual nuclei in high energy hadron-nucleus interactions , 1995 .

[64]  I Rinaldi,et al.  An integral test of FLUKA nuclear models with 160 MeV proton beams in multi-layer Faraday cups , 2011, Physics in medicine and biology.

[65]  N. Mott The Scattering of Fast Electrons by Atomic Nuclei , 1929 .

[66]  A Mairani,et al.  Benchmarking nuclear models of FLUKA and GEANT4 for carbon ion therapy , 2010, Physics in medicine and biology.

[67]  Horst Stöcker,et al.  Relativistic quantum molecular dynamics approach to nuclear collisions at ultrarelativistic energies , 1989 .

[68]  W. Lehnert,et al.  Count Rate Performance of the MicroPET Focus 220 Animal Scanner in Singles Transmission Scanning Mode , 2008, IEEE Transactions on Nuclear Science.

[69]  Chan Hyeong Kim,et al.  Determination of the distal dose edge in a human phantom by measuring the prompt gamma distribution: A monte carlo study , 2010 .

[70]  G. Moliere,et al.  Theory of the scattering of fast charged particles. 2. Repeated and multiple scattering , 1948 .

[71]  A Ottolenghi,et al.  The FLUKA code for space applications: recent developments. , 2004, Advances in space research : the official journal of the Committee on Space Research.

[72]  Julia Bauer,et al.  Post-therapeutical β+-activity measurements in comparison to simulations towards in-vivo verification of ion beam therapy , 2011, 2011 IEEE Nuclear Science Symposium Conference Record.

[73]  Fasso New developments in FLUKA modeling hadronic and EM interactions , 1997 .

[74]  H Paganetti,et al.  Adaptation of GEANT4 to Monte Carlo dose calculations based on CT data. , 2004, Medical physics.

[75]  M Scholz,et al.  Rapid calculation of biological effects in ion radiotherapy , 2006, Physics in medicine and biology.

[76]  John N. Dyer,et al.  Resolution of the Σ - -Mass Anomaly , 1963 .

[77]  P. Busca,et al.  Prompt gamma imaging with a slit camera for real-time range control in proton therapy: Experimental validation up to 230 MeV with HICAM and development of a new prototype , 2012 .

[78]  A. Ferrari,et al.  Describing Compton scattering and two-quanta positron annihilation based on Compton profiles: two models suited for the Monte Carlo method , 2012 .

[79]  J. Ziegler The stopping and range of ions in solids vol 1 : The stopping and ranges of ions in matter , 2013 .

[80]  Katia Parodi,et al.  Phase Space Generation for Proton and Carbon Ion Beams for External Users’ Applications at the Heidelberg Ion Therapy Center , 2016, Front. Oncol..

[81]  Sommerer,et al.  Experiments and FLUKA simulations of $^{12}C$ and $^{16}O$ beams for therapy monitoring by means of in-beam Positron Emission Tomography , 2007 .

[82]  James F. Ziegler,et al.  Helium: Stopping Powers and Ranges in All Elemental Matter , 1977 .

[83]  K Parodi,et al.  Dosimetric accuracy assessment of a treatment plan verification system for scanned proton beam radiotherapy: one-year experimental results and Monte Carlo analysis of the involved uncertainties , 2013, Physics in medicine and biology.

[84]  M. Scholz,et al.  Computation of cell survival in heavy ion beams for therapy , 1997, Radiation and environmental biophysics.

[85]  Jinsong Ouyang,et al.  Monitoring proton radiation therapy with in-room PET imaging , 2011, Physics in medicine and biology.

[86]  A Ottolenghi,et al.  Role of shielding in modulating the effects of solar particle events: Monte Carlo calculation of absorbed dose and DNA complex lesions in different organs. , 2004, Advances in space research : the official journal of the Committee on Space Research.

[87]  S. Kailas,et al.  RIPL – Reference Input Parameter Library for Calculation of Nuclear Reactions and Nuclear Data Evaluations , 2009 .

[88]  Chan Hyeong Kim,et al.  Preliminary Study for Determination of Distal Dose Edge by Measuring 90-deg Prompt Gammas with an Array-Type Prompt Gamma Detection System , 2009 .

[89]  M. Bajard,et al.  Monte Carlo Simulations of Prompt-Gamma Emission During Carbon Ion Irradiation , 2009, IEEE Transactions on Nuclear Science.

[90]  P. R. Sala,et al.  FLUKA : Performances and Applications in the Intermediate Energy Range , 2022 .

[91]  P. R. Sala,et al.  The production of residual nuclei in peripheral high energy nucleus-nucleus interactions , 1996 .

[92]  A. Gottlob,et al.  Monte-Carlo study , 1998 .

[93]  H Paganetti,et al.  Clinical CT-based calculations of dose and positron emitter distributions in proton therapy using the FLUKA Monte Carlo code , 2007, Physics in medicine and biology.

[94]  K Parodi,et al.  First experimental-based characterization of oxygen ion beam depth dose distributions at the Heidelberg Ion-Beam Therapy Center , 2012, Physics in medicine and biology.

[95]  Katia Parodi,et al.  PET/CT imaging for treatment verification after proton therapy: a study with plastic phantoms and metallic implants. , 2007, Medical physics.

[96]  A. Ferrari,et al.  FLUKA: Status and Prospects for Hadronic Applications , 2001 .

[97]  Katia Parodi,et al.  Patient study of in vivo verification of beam delivery and range, using positron emission tomography and computed tomography imaging after proton therapy. , 2007, International journal of radiation oncology, biology, physics.

[98]  A. Ferrari,et al.  The Physics of High Energy Reactions , 1997 .

[99]  Zuofeng Li,et al.  Prescribing, Recording, and Reporting Proton-Beam Therapy , 2009 .

[100]  W. Kalender,et al.  Correlation between CT numbers and tissue parameters needed for Monte Carlo simulations of clinical dose distributions , 2000 .

[101]  T. Tessonnier,et al.  Development and application of tools for Monte Carlo based simulations in a particle beam radiotherapy facility. , 2014, Applied radiation and isotopes : including data, instrumentation and methods for use in agriculture, industry and medicine.

[102]  A Ottolenghi,et al.  A Monte Carlo code for a direct estimation of radiation risk. , 2001, Physica medica : PM : an international journal devoted to the applications of physics to medicine and biology : official journal of the Italian Association of Biomedical Physics.

[103]  R. H. Ritchie,et al.  Z13-DEPENDENT STOPPING POWER AND RANGE CONTRIBUTIONS* , 1973 .

[104]  Christopher Kurz,et al.  Implementation and initial clinical experience of offline PET/CT-based verification of scanned carbon ion treatment. , 2013, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[105]  P. W. Chin,et al.  Overview of the FLUKA code , 2014, ICS 2014.

[106]  J. F. Ziegler,et al.  III – Status of Experimental Data , 1977 .

[107]  D. Schardt,et al.  Magnetic scanning system for heavy ion therapy , 1993 .

[108]  Chan Hyeong Kim,et al.  Experimental performance of double-scattering Compton camera with anthropomorphic phantom , 2011 .

[109]  J Fieres,et al.  Experimental characterization of lateral profiles of scanned proton and carbon ion pencil beams for improved beam models in ion therapy treatment planning , 2011, Physics in medicine and biology.

[110]  Shuryak,et al.  Screening of the topological charge in a correlated instanton vacuum. , 1995, Physical review. D, Particles and fields.

[111]  M. J. Berger,et al.  Density effect for the ionization loss of charged particles in various substances , 1984 .

[112]  C. Stearns,et al.  Investigation of the count rate performance of the General Electric Advance positron emission tomograph , 1994, Proceedings of 1994 IEEE Nuclear Science Symposium - NSS'94.

[113]  Space life sciences: radiation risk assessment and radiation measurements in low Earth orbit. , 2004, Advances in space research : the official journal of the Committee on Space Research.

[114]  Akifumi Fukumura,et al.  Measurements of secondary neutrons produced from thick targets bombarded by high-energy helium and carbon ions , 1999 .

[115]  Frederic H Fahey,et al.  Data acquisition in PET imaging. , 2002, Journal of nuclear medicine technology.

[116]  A Mairani,et al.  Dosimetric commissioning and quality assurance of scanned ion beams at the Italian National Center for Oncological Hadrontherapy. , 2015, Medical physics.

[117]  Jong-Won Kim,et al.  Simulation Study on a Prompt Gamma Detection System for Use in Proton Therapy , 2008 .

[118]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

[119]  A. Ferrari,et al.  FLUKA: A Multi-Particle Transport Code , 2005 .

[120]  Chan Hyeong Kim,et al.  Determination of Optimal Energy Window for Measurement of Prompt Gammas from Proton Beam by Monte Carlo Simulations , 2008 .

[121]  J. Ziegler Stopping of energetic light ions in elemental matter , 1999 .

[122]  Albert van der Kogel,et al.  Basic Clinical Radiobiology Fourth Edition , 2009 .

[123]  Vasilis Vlachoudis,et al.  Numerically robust geometry engine for compound solid geometries , 2014, ICS 2014.

[124]  E. Gadioli,et al.  Boltzmann master equation theory of angular distributions in heavy-ion reactions , 1998 .

[125]  P. V. Vavilov IONIZATION LOSSES OF HIGH-ENERGY HEAVY PARTICLES , 1957 .

[126]  Katia Parodi,et al.  Time-of-flight neutron rejection to improve prompt gamma imaging for proton range verification: a simulation study , 2012, Physics in medicine and biology.

[127]  A. Mairani,et al.  The physics of the FLUKA code: Recent developments , 2007 .

[128]  H Paganetti,et al.  Monte Carlo patient study on the comparison of prompt gamma and PET imaging for range verification in proton therapy , 2011, Physics in medicine and biology.