The FLUKA Monte Carlo code coupled with the local effect model for biological calculations in carbon ion therapy

Clinical Monte Carlo (MC) calculations for carbon ion therapy have to provide absorbed and RBE-weighted dose. The latter is defined as the product of the dose and the relative biological effectiveness (RBE). At the GSI Helmholtzzentrum für Schwerionenforschung as well as at the Heidelberg Ion Therapy Center (HIT), the RBE values are calculated according to the local effect model (LEM). In this paper, we describe the approach followed for coupling the FLUKA MC code with the LEM and its application to dose and RBE-weighted dose calculations for a superimposition of two opposed (12)C ion fields as applied in therapeutic irradiations. The obtained results are compared with the available experimental data of CHO (Chinese hamster ovary) cell survival and the outcomes of the GSI analytical treatment planning code TRiP98. Some discrepancies have been observed between the analytical and MC calculations of absorbed physical dose profiles, which can be explained by the differences between the laterally integrated depth-dose distributions in water used as input basic data in TRiP98 and the FLUKA recalculated ones. On the other hand, taking into account the differences in the physical beam modeling, the FLUKA-based biological calculations of the CHO cell survival profiles are found in good agreement with the experimental data as well with the TRiP98 predictions. The developed approach that combines the MC transport/interaction capability with the same biological model as in the treatment planning system (TPS) will be used at HIT to support validation/improvement of both dose and RBE-weighted dose calculations performed by the analytical TPS.

[1]  K. Parodi,et al.  Clinical calculations of physical and biological effective dose distributions in proton and carbon ion therapy using the FLUKA Monte Carlo code , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[2]  A. Kellerer,et al.  A generalized formulation of dual radiation action. , 2012, Radiation research.

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

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

[5]  D. Schardt,et al.  Experimental fragmentation studies with 12C therapy beams. , 2006, Radiation protection dosimetry.

[6]  A. Brahme,et al.  Depth absorbed dose and LET distributions of therapeutic 1H, 4He, 7Li, and 12C beams. , 2007, Medical physics.

[7]  Igor Mishustin,et al.  Neutrons from fragmentation of light nuclei in tissue-like media: a study with the GEANT4 toolkit. , 2005, Physics in medicine and biology.

[8]  Harold O. Wyckoff,et al.  International Commission ON Radiation Units and Measurements (ICRU). , 1974, The American journal of roentgenology, radium therapy, and nuclear medicine.

[9]  A Gemmel,et al.  Biological dose optimization with multiple ion fields , 2008, Physics in medicine and biology.

[10]  M. Pelliccioni,et al.  On the Conversion Coefficients from Fluence to Ambient Dose Equivalent , 1994 .

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

[12]  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.

[13]  M Scholz,et al.  Treatment planning for heavy-ion radiotherapy: calculation and optimization of biologically effective dose. , 2000, Physics in medicine and biology.

[14]  Katia Parodi,et al.  Investigating the accuracy of the FLUKA code for transport of therapeutic ion beams in matter , 2006, Physics in medicine and biology.

[15]  K Niita,et al.  Simulations of an accelerator-based shielding experiment using the particle and heavy-ion transport code system PHITS. , 2005, Advances in space research : the official journal of the Committee on Space Research.

[16]  O Jäkel,et al.  Treatment planning for heavy-ion radiotherapy: physical beam model and dose optimization. , 2000, Physics in medicine and biology.

[17]  A. Brahme,et al.  Ion beam transport in tissue-like media using the Monte Carlo code SHIELD-HIT. , 2004, Physics in medicine and biology.

[18]  M Krämer,et al.  Biological dosimetry of complex ion radiation fields. , 2003, Physics in medicine and biology.

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

[20]  U Oelfke,et al.  Measurements of relative biological effectiveness of the 70 MeV proton beam at TRIUMF using Chinese hamster V79 cells and the high-precision cell sorter assay. , 1996, Radiation research.

[21]  M Scholz,et al.  Track structure and the calculation of biological effects of heavy charged particles. , 1996, Advances in space research : the official journal of the Committee on Space Research.

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

[23]  M. Scholz,et al.  The Increased Biological Effectiveness of Heavy Charged Particles: From Radiobiology to Treatment Planning , 2003, Technology in cancer research & treatment.

[24]  Michael Scholz,et al.  Biophysical calculation of cell survival probabilities using amorphous track structure models for heavy-ion irradiation , 2008, Physics in medicine and biology.

[25]  Hiroshi Nakashima,et al.  PHITS: A particle and heavy ion transport code system , 2006 .

[26]  O Jäkel,et al.  Treatment planning for heavy ion radiotherapy: clinical implementation and application. , 2001, Physics in medicine and biology.

[27]  F. Cerutti,et al.  The FLUKA code: Description and benchmarking , 2007 .