No Arabic abstract
In a cyclotron-based proton therapy facility, the energy changes are performed by means of a degrader of variable thickness. The interaction of the proton beam with the degrader creates energy tails and increases the beam emittance. A precise model of the degraded beam properties is important not only to better understand the performance of a facility already in operation, but also to support the development of new proton therapy concepts. The exact knowledge of the degraded beam properties, in terms of energy spectrum and transverse phase space, depends on the model used to describe the proton interaction with the degrader material. In this work the model of a graphite degrader has been developed with four Monte Carlo codes: three conventional Monte Carlo codes (FLUKA, GEANT4 and MCNPX) and the multi-purpose particle tracking code OPAL equipped with a simplified Monte Carlo routine. From the comparison between the different codes, we can deduce how the accuracy of the degrader model influences the precision of the beam dynamics model of a possible transport line downstream of the degrader.
Despite the fact that the first-order beam dynamics models allow an approximated evaluation of the beam properties, their contribution is essential during the conceptual design of an accelerator or beamline. However, during the commissioning some of their limitations appear in the comparison against measurements. The extension of the linear model to higher order effects is, therefore, demanded. In this paper, the effects of particle-matter interaction have been included in the model of the transport lines in the proton therapy facility at the Paul Scherrer Institut (PSI) in Switzerland. To improve the performance of the facility, a more precise model was required and has been developed with the multi-particle open source beam dynamics code called OPAL (Object oriented Particle Accelerator Library). In OPAL, the Monte Carlo simulations of Coulomb scattering and energy loss are performed seamless with the particle tracking. Beside the linear optics, the influence of the passive elements (e.g. degrader, collimators, scattering foils and air gaps) on the beam emittance and energy spread can be analysed in the new model. This allows for a significantly improved precision in the prediction of beam transmission and beam properties. The accuracy of the OPAL model has been confirmed by numerous measurements.
We study the propagation of nucleons and nuclei in tissue-like media within a Monte Carlo Model for Heavy-ion Therapy (MCHIT) based on the GEANT4 toolkit (version 8.2). The model takes into account fragmentation of projectile nuclei and secondary interactions of produced nuclear fragments. Model predictions are validated with available experimental data obtained for water and PMMA phantoms irradiated by monoenergetic carbon-ion beams. The MCHIT model describes well (1) the depth-dose distributions in water and PMMA, (2) the doses measured for fragments of certain charge, (3) the distributions of positron emitting nuclear fragments produced by carbon-ion beams, and (4) the energy spectra of secondary neutrons measured at different angles to the beam direction. Radial dose profiles for primary nuclei and for different projectile fragments are calculated and discussed as possible input for evaluation of biological dose distributions. It is shown that at the periphery of the transverse dose profile close to the Bragg peak the dose from secondary nuclear fragments is comparable to the dose from primary nuclei.
Cone beam CT (CBCT) has been widely used for patient setup in image guided radiation therapy (IGRT). Radiation dose from CBCT scans has become a clinical concern. The purposes of this study are 1) to commission a GPU-based Monte Carlo (MC) dose calculation package gCTD for Varian On-Board Imaging (OBI) system and test the calculation accuracy, and 2) to quantitatively evaluate CBCT dose from the OBI system in typical IGRT scan protocols. We first conducted dose measurements in a water phantom. X-ray source model parameters used in gCTD are obtained through a commissioning process. gCTD accuracy is demonstrated by comparing calculations with measurements in water and in CTDI phantoms. 25 brain cancer patients are used to study dose in a standard-dose head protocol, and 25 prostate cancer patients are used to study dose in pelvis protocol and pelvis spotlight protocol. Mean dose to each organ is calculated. Mean dose to 2% voxels that have the highest dose is also computed to quantify the maximum dose. It is found that the mean dose value to an organ varies largely among patients. Moreover, dose distribution is highly non-homogeneous inside an organ. The maximum dose is found to be 1~3 times higher than the mean dose depending on the organ, and is up to 8 times higher for the entire body due to the very high dose region in bony structures. High computational efficiency has also been observed in our studies, such that MC dose calculation time is less than 5 min for a typical case.
During the proton-anti proton collider run several experiments were carried out in order to understand the effect of the beam-beam interaction on backgrounds and lifetimes. In this talk a selection of these experiments will be presented. From these experiments, the importance of relative beam sizes and tune ripple could be demonstrated.
Purpose: To assess the effects of brain movements induced by heartbeat on dose distributions in synchrotron micro- and mini-beam radiaton therapy and to develop a model to help guide decisions and planning for future clinical trials. Methods: The Monte Carlo code PENELOPE was used to simulate the irradiation of a human head phantom with a variety of micro- and mini-beam arrays, with beams narrower than 100mum and above 500mum, respectively, and with radiation fields of 1cmx2cm and 2cmx2cm. The dose in the phantom due to these beams was calculated by superposing the dose profiles obtained for a single beam of 1mumx2cm. A parameter delta, accounting for the total displacement of the brain during the irradiation and due to the cardio-synchronous pulsation, was used to quantify the impact on peak-to-valley dose ratios and the full-width at half-maximum. Results: The difference between the maximum (at the phantom entrance) and the minimum (at the phantom exit) values of the peak-to-valley dose ratio reduces when the parameter $delta$ increases. The full-width at half-maximum remains almost constant with depth for any $delta$ value. Sudden changes in the two quantities are observed at the interfaces between the various tissues (brain, skull and skin) present in the head phantom. The peak-to-valley dose ratio at the center of the head phantom reduces when delta increases, remaining above 70% of the static value only for mini-beams and delta smaller than ~200mum. Conclusions: Optimal setups for brain treatments with synchrotron radiation micro- and mini-beam combs depend on the brain displacement due to cardio-synchronous pulsation. Peak-to-valley dose ratios larger than 90% of the maximum values obtained in the static case occur only for mini-beams and relatively large dose rates.