Changran Geng

Analysis of the track- and dose-averaged LET and LET spectra in proton therapy using the GEANT4 Monte Carlo code

PURPOSE The motivation of this study was to find and eliminate the cause of errors in dose-averaged linear energy transfer (LET) calculations from therapeutic protons in small targets, such as biological cell layers, calculated using the geant 4 Monte Carlo code. Furthermore, the purpose was also to provide a recommendation to select an appropriate LET quantity from geant 4 simulations to correlate with biological effectiveness of therapeutic protons. METHODS The authors developed a particle tracking step based strategy to calculate the average LET quantities (track-averaged LET, LETt and dose-averaged LET, LETd) using geant 4 for different tracking step size limits. A step size limit refers to the maximally allowable tracking step length. The authors investigated how the tracking step size limit influenced the calculated LETt and LETd of protons with six different step limits ranging from 1 to 500 μm in a water phantom irradiated by a 79.7-MeV clinical proton beam. In addition, the authors analyzed the detailed stochastic energy deposition information including fluence spectra and dose spectra of the energy-deposition-per-step of protons. As a reference, the authors also calculated the averaged LET and analyzed the LET spectra combining the Monte Carlo method and the deterministic method. Relative biological effectiveness (RBE) calculations were performed to illustrate the impact of different LET calculation methods on the RBE-weighted dose. RESULTS Simulation results showed that the step limit effect was small for LETt but significant for LETd. This resulted from differences in the energy-deposition-per-step between the fluence spectra and dose spectra at different depths in the phantom. Using the Monte Carlo particle tracking method in geant 4 can result in incorrect LETd calculation results in the dose plateau region for small step limits. The erroneous LETd results can be attributed to the algorithm to determine fluctuations in energy deposition along the tracking step in geant 4. The incorrect LETd values lead to substantial differences in the calculated RBE. CONCLUSIONS When the geant 4 particle tracking method is used to calculate the average LET values within targets with a small step limit, such as smaller than 500 μm, the authors recommend the use of LETt in the dose plateau region and LETd around the Bragg peak. For a large step limit, i.e., 500 μm, LETd is recommended along the whole Bragg curve. The transition point depends on beam parameters and can be found by determining the location where the gradient of the ratio of LETd and LETt becomes positive.

Validation of the radiobiology toolkit TOPAS-nBio in simple DNA geometries

Computational simulations offer a powerful tool for quantitatively investigating radiation interactions with biological tissue and can help bridge the gap between physics, chemistry and biology. The TOPAS collaboration is tackling this challenge by extending the current Monte Carlo tool to allow for sub-cellular in silico simulations in a new extension, TOPAS-nBio. TOPAS wraps and extends the Geant4 Monte Carlo simulation toolkit and the new extension allows the modeling of particles down to vibrational energies (∼2eV) within realistic biological geometries. Here we present a validation of biological geometries available in TOPAS-nBio, by comparing our results to two previously published studies. We compare the prediction of strand breaks in a simple linear DNA strand from TOPAS-nBio to a published Monte Carlo track structure simulation study. While TOPAS-nBio confirms the trend in strand break generation, it predicts a higher frequency of events below an energy of 17.5eV compared to the alternative Monte Carlo track structure study. This is due to differences in the physics models used by each code. We also compare the experimental measurement of strand breaks from incident protons in DNA plasmids to TOPAS-nBio simulations. Our results show good agreement of single and double strand breaks predicting a similar increase in the strand break yield with increasing LET.

Assessing the radiation-induced second cancer risk in proton therapy for pediatric brain tumors: the impact of employing a patient-specific aperture in pencil beam scanning.

The purpose of this study was to compare the radiation-induced second cancer risks for in-field and out-of-field organs and tissues for pencil beam scanning (PBS) and passive scattering proton therapy (PPT) and assess the impact of adding patient-specific apertures to sharpen the penumbra in pencil beam scanning for pediatric brain tumor patients. Five proton therapy plans were created for each of three pediatric patients using PPT as well as PBS with two spot sizes (average sigma of ~17 mm and ~8 mm at isocenter) and choice of patient-specific apertures. The lifetime attributable second malignancy risks for both in-field and out-of-field tissues and organs were compared among five delivery techniques. The risk for in-field tissues was calculated using the organ equivalent dose, which is determined by the dose volume histogram. For out-of-field organs, the organ-specific dose equivalent from secondary neutrons was calculated using Monte Carlo and anthropomorphic pediatric phantoms. We find that either for small spot size PBS or for large spot size PBS, a patient-specific aperture reduces the in-field cancer risk to values lower than that for PPT. The reduction for large spot sizes (on average 43%) is larger than for small spot sizes (on average 21%). For out-of-field organs, the risk varies only marginally by employing a patient-specific aperture (on average from  -2% to 16% with increasing distance from the tumor), but is still one to two orders of magnitude lower than that for PPT. In conclusion, when pencil beam spot sizes are large, the addition of apertures to sharpen the penumbra decreases the in-field radiation-induced secondary cancer risk. There is a slight increase in out-of-field cancer risk as a result of neutron scatter from the aperture, but this risk is by far outweighed by the in-field risk benefit from using an aperture with a large PBS spot size. In general, the risk for developing a second malignancy in out-of-field organs for PBS remains much lower compared to PPT even if apertures are being applied.

GEANT4 calculations of neutron dose in radiation protection using a homogeneous phantom and a Chinese hybrid male phantom

The purpose of this study is to verify the feasibility of applying GEANT4 (version 10.01) in neutron dose calculations in radiation protection by comparing the calculation results with MCNP5. The depth dose distributions are investigated in a homogeneous phantom, and the fluence-to-dose conversion coefficients are calculated for different organs in the Chinese hybrid male phantom for neutrons with energy ranging from 1 × 10(-9) to 10 MeV. By comparing the simulation results between GEANT4 and MCNP5, it is shown that using the high-precision (HP) neutron physics list, GEANT4 produces the closest simulation results to MCNP5. However, differences could be observed when the neutron energy is lower than 1 × 10(-6) MeV. Activating the thermal scattering with an S matrix correction in GEANT4 with HP and MCNP5 in thermal energy range can reduce the difference between these two codes.

Monte Carlo study of the beam shaping assembly optimization for providing high epithermal neutron flux for BNCT based on D–T neutron generator

In this study, a beam shaping assembly with high epithermal neutron flux output was designed based on a D–T neutron generator using Monte Carlo N-particle Transport Code. D2O-54Fe and AlF3-60Ni interlayer moderator, efficient multiplier, filters, and reflector were used to improve the neutron beam quality according to the requirements of boron neutron capture therapy while maintaining high flux of epithermal neutron beam. In addition, the dose performance of the beam from our proposed facility was assessed in the Snyder head phantom. The simulation results proved that the proposed neutron beam was applicable to the treatment of deep-seated brain tumor.

Determination of the relationship between dose deposition and Cerenkov photons in homogeneous and heterogeneous phantoms during radiotherapy using Monte Carlo method

To explore the response relationship between Cerenkov photons and dose deposition, which is the theoretical premise of dose measurement based on Cerenkov effect, Geant4 was used to simulate the process of homogeneous and heterogeneous phantoms irradiated with monoenergetic or clinical beam energy spectrum. Results showed certain response relationship between Cerenkov photons and dose deposition regardless of the phantom used under different irradiation conditions. However, deviation observed in the axial distribution of dose deposition as characterized by Cerenkov photons was larger under electron beam than under photon beam.

Dose assessment for the fetus considering scattered and secondary radiation from photon and proton therapy when treating a brain tumor of the mother.

The goal of this work was to determine the scattered photon dose and secondary neutron dose and resulting risk for the sensitive fetus from photon and proton radiotherapy when treating a brain tumor during pregnancy. Anthropomorphic pregnancy phantoms with three stages (3-, 6-, 9-month) based on ICRP reference parameters were implemented in Monte Carlo platform TOPAS, to evaluate the scattered dose and secondary neutron dose and dose equivalent. To evaluate the dose equivalent, dose averaged quality factors were considered for neutrons. This study compared three treatment modalities: passive scattering and pencil beam scanning proton therapy (PPT and PBS) and 6-MV 3D conformal photon therapy. The results show that, for 3D conformal photon therapy, the scattered photon dose equivalent to the fetal body increases from 0.011 to 0.030 mSv per treatment Gy with increasing stage of gestation. For PBS, the neutron dose equivalent to the fetal body was significantly lower, i.e. increasing from 1.5 × 10(-3) to 2.5 × 10(-3) mSv per treatment Gy with increasing stage of gestation. For PPT, the neutron dose equivalent of the fetus decreases from 0.17 to 0.13 mSv per treatment Gy with the growing fetus. The ratios of dose equivalents to the fetus for a 52.2 Gy(RBE) course of radiation therapy to a typical CT scan of the mothers head ranged from 3.4-4.4 for PBS, 30-41 for 3D conformal photon therapy and 180-500 for PPT, respectively. The attained dose to a fetus from the three modalities is far lower than the thresholds of malformation, severe mental retardation and lethal death. The childhood cancer excessive absolute risk was estimated using a linear no-threshold dose-response relationship. The risk would be 1.0 (95% CI: 0.6, 1.6) and 0.1 (95% CI: -0.01, 0.52) in 10(5) for the 9-month fetus for PBS with a prescribed dose of 52.2 Gy(RBE). The increased risks for PPT and photon therapy are about two and one orders of magnitude larger than that for PBS, respectively. We can conclude that a pregnant woman with a brain tumor could be treated with pencil beam scanning with acceptable risks to the fetus.

A Monte Carlo-based radiation safety assessment for astronauts in an environment with confined magnetic field shielding

The active shielding technique has great potential for radiation protection in space exploration because it has the advantage of a significant mass saving compared with the passive shielding technique. This paper demonstrates a Monte Carlo-based approach to evaluating the shielding effectiveness of the active shielding technique using confined magnetic fields (CMFs). The International Commission on Radiological Protection reference anthropomorphic phantom, as well as the toroidal CMF, was modeled using the Monte Carlo toolkit Geant4. The penetrating primary particle fluence, organ-specific dose equivalent, and male effective dose were calculated for particles in galactic cosmic radiation (GCR) and solar particle events (SPEs). Results show that the SPE protons can be easily shielded against, even almost completely deflected, by the toroidal magnetic field. GCR particles can also be more effectively shielded against by increasing the magnetic field strength. Our results also show that the introduction of a structural Al wall in the CMF did not provide additional shielding for GCR; in fact it can weaken the total shielding effect of the CMF. This study demonstrated the feasibility of accurately determining the radiation field inside the environment and evaluating the organ dose equivalents for astronauts under active shielding using the CMF.

Prediction of Treatment Response for Combined Chemo- and Radiation Therapy for Non-Small Cell Lung Cancer Patients Using a Bio-Mathematical Model

The goal of this work was to develop a mathematical model to predict Kaplan–Meier survival curves for chemotherapy combined with radiation in Non-Small Cell Lung Cancer patients for use in clinical trial design. The Gompertz model was used to describe tumor growth, radiation effect was simulated by the linear-quadratic model with an α/β-ratio of 10, and chemotherapy effect was based on the log-cell kill model. To account for repopulation during treatment, we considered two independent methods: 1) kickoff-repopulation using exponential growth with a decreased volume doubling time, or 2) Gompertz-repopulation using the gradually accelerating growth rate with tumor shrinkage. The input parameters were independently estimated by fitting to the SEER database for untreated tumors, RTOG-8808 for radiation only, and RTOG-9410 for sequential chemo-radiation. Applying the model, the benefit from concurrent chemo-radiation comparing to sequential for stage III patients was predicted to be a 6.6% and 6.2% improvement in overall survival for 3 and 5-years respectively, comparing well to the 5.3% and 4.5% observed in RTOG-9410. In summary, a mathematical model was developed to model tumor growth over extended periods of time, and can be used for the optimization of combined chemo-radiation scheduling and sequencing.

Influence of Neutron Sources and 10B Concentration on Boron Neutron Capture Therapy for Shallow and Deeper Non-small Cell Lung Cancer.

Abstract Boron Neutron Capture Therapy (BNCT) is a radiotherapy that combines biological targeting and high Linear Energy Transfer (LET). It is considered a potential therapeutic approach for non-small cell lung cancer (NSCLC). It could avoid the inaccurate treatment caused by the lung motion during radiotherapy, because the dose deposition mainly depends on the boron localization and neutron source. Thus, 10B concentration and neutron sources are both principal factors of BNCT, and they play significant roles in the curative effect of BNCT for different cases. The purpose was to explore the feasibility of BNCT treatment for NSCLC with either of two neutron sources (the epithermal reactor at the Massachusetts Institute of Technology named “MIT source” and the accelerator neutron source designed in Argentina named “MEC source”) and various boron concentrations. Shallow and deeper lung tumors were defined in the Chinese hybrid radiation phantom, and the Monte Carlo method was used to calculate the dose to tumors and healthy organs. The MEC source was more appropriate to treat the shallow tumor (depth of 6 cm) with a shorter treatment time. However, the MIT source was more suitable for deep lung tumor (depth of 9 cm) treatment, as the MEC source is more likely to exceed the skin dose limit. Thus, a neutron source consisting of more fast neutrons is not necessarily suitable for deep treatment of lung tumors. Theoretical distribution of 10B in tumors and organs at risk (especially skin) was obtained to meet the treatable requirement of BNCT, which may provide the references to identify the feasibility of BNCT for the treatment of lung cancer using these two neutron sources in future clinical applications.