No Arabic abstract
We study energy deposition by light nuclei in tissue-like media taking into account nuclear fragmentation reactions, in particular, production of secondary neutrons. The calculations are carried out within a Monte Carlo model for Heavy-Ion Therapy (MCHIT) based on the GEANT4 toolkit. Experimental data on depth-dose distributions for 135A-400A MeV C-12 and O-18 beams are described very well without any adjustment of the model parameters. This gives confidence in successful use of the GEANT4 toolkit for MC simulations of cancer therapy with beams of light nuclei. The energy deposition due to secondary neutrons produced by C-12 and Ne-20 beams in a (40-50 cm)^3 water phantom is estimated to 1-2% of the total dose, that is only slightly above the neutron contribution (~1%) induced by a 200 MeV proton beam.
It is well known from numerous experiments that nuclear multifragmentation is a dominating mechanism for production of intermediate-mass fragments in nucleus-nucleus collisions at energies above 100 A MeV. In this paper we investigate the validity and performance of the Fermi break-up model and the statistical multifragmentation model implemented as parts of the Geant4 toolkit. We study the impact of violent nuclear disintegration reactions on the depth-dose profiles and yields of secondary fragments for beams of light and medium-weight nuclei propagating in extended media. Implications for ion-beam cancer therapy and shielding from cosmic radiation are discussed.
We study the energy deposition by light and heavy nuclei in tissue-like media as used for cancer therapy. The depth-dose distributions for protons, $^{3}$He, $^{12}$C, $^{20}$Ne, and $^{58}$Ni nuclei are calculated within a Monte Carlo model based on the GEANT4 toolkit. These distributions are compared with each other and with available experimental data. It is demonstrated that nuclear fragmentation reactions essentially reduce the peak-to-plateau ratio of the dose profiles for deeply penetrating energetic ions heavier than $^{3}$He. On the other hand, all projectiles up to $^{20}$Ne were found equally suitable for therapeutic use at low penetration depths.
We model the responses of Tissue-Equivalent Proportional Counters (TEPC) to radiation fields of therapeutic C-12 beams in a water phantom and to quasi-monoenergetic neutrons in a PMMA phantom. Simulations are performed with the Monte Carlo model for Heavy Ion Therapy (MCHIT) based on the Geant4 toolkit. The shapes of the calculated lineal energy spectra agree well with measurements in both cases. The influence of fragmentation reactions on the TEPC response to a narrow pencil-like beam with its width smaller than the TEPC diameter is investigated by Monte Carlo modeling. It is found that total lineal energy spectra are not very sensitive to the choice of the nuclear fragmentation model used. The calculated frequency-mean lineal energy y_f differs from the data on the axis of a therapeutic beam by less than 10% and by 10-20% at other TEPC positions. The validation of MCHIT with neutron beams gives us confidence in estimating the contributions to lineal energy spectra due to secondary neutrons produced in water by C-12 nuclei. As found, the neutron contribution at 10 cm distance from the beam axis amounts to ~ 50% close the entrance to the phantom and decreases to ~ 25% at the depth of the Bragg peak and beyond it. The presented results can help in evaluating biological out-of-field doses in carbon-ion therapy.
We study the spatial distributions of $beta^+$-activity produced by therapeutic beams of $^3$He and $^{12}$C ions in various tissue-like materials. The calculations were performed within a Monte Carlo model for Heavy-Ion Therapy (MCHIT) based on the GEANT4 toolkit. The contributions from $^{10,11}$C, $^{13}$N, $^{14,15}$O, $^{17,18}$F and $^{30}$P positron-emitting nuclei were calculated and compared with experimental data obtained during and after irradiation. Positron emitting nuclei are created by $^{12}$C beam in fragmentation reactions of projectile and target nuclei. This leads to a $beta^+$-activity profile characterised by a noticeable peak located close to the Bragg peak in the corresponding depth-dose distribution. On the contrary, as the most of positron-emitting nuclei are produced by $^3$He beam in target fragmentation reactions, the calculated total $beta^+$-activity during or soon after the irradiation period is evenly distributed within the projectile range. However, we predict also the presence of $^{13}$N, $^{14}$O, $^{17,18}$F created in charge-transfer reactions by low-energy $^3$He ions close to the end of their range in several tissue-like media. The time evolution of $beta^+$-activity profiles was investigated for both kinds of beams. Due to the production of $^{18}$F nuclide the $beta^+$-activity profile measured 2 or 3 hours after irradiation with $^{3}$He ions will have a distinct peak correlated with the maximum of depth-dose distribution. We found certain advantages of low-energy $^{3}$He beams over low-energy proton beams for reliable PET monitoring during particle therapy of shallow located tumours. In this case the distal edge of $beta^+$-activity distribution from $^{17}$F nuclei clearly marks the range of $^{3}$He in tissues.
Over the past few decades, researchers have developed several approaches such as the Reference Phantom Method (RPM) to estimate ultrasound attenuation coefficient (AC) and backscatter coefficient (BSC). AC and BSC can help to discriminate pathology from normal tissue during in-vivo imaging. In this paper, we propose a new RPM model to simultaneously compute AC and BSC for harmonic imaging and a normalized score that combines the two parameters as a measure of disease progression. The model utilizes the spectral difference between two regions of interest, the first, a proximal, close to the probe and second, a distal, away from the probe. We have implemented an algorithm based on the model and shown that it provides accurate and stable estimates to within 5% of AC and BSC for simulated received echo from post-focal depths of a homogeneous liver-like medium. For practical applications with time gain and time frequency compensated in-phase and quadrature (IQ) data from ultrasound scanner, the method has been approximated and generalized to estimate AC and BSC for tissue layer underlying a more attenuative subcutaneous layer. The angular spectrum approach for ultrasound propagation in biological tissue is employed as a virtual Reference Phantom (VRP). The VRP is calibrated with a fixed probe and scanning protocol for application to liver tissue. In a feasibility study with 16 subjects, the method is able to separate 9/11 cases of progressive non-alcoholic fatty liver disease from 5 normal. In particular, it is able to separate 4/5 cases of non-alcoholic steato-hepatitis and early fibrosis (F<=2) from normal tissue. More extensive clinical studies are needed to assess the full capability of this model for screening and monitoring disease progression in liver and other tissues.