No Arabic abstract
Background: Treatment verification with PET imaging in charged particle therapy is conventionally done by comparing measurements of spatial distributions with Monte Carlo (MC) predictions. However, decay curves can provide additional independent information about the treatment and the irradiated tissue. Most studies performed so far focus on long time intervals. Here we investigate the reliability of MC predictions of space and time (decay rate) profiles shortly after irradiation, and we show how the decay rates can give an indication about the elements of which the phantom is made up. Methods and Materials: Various phantoms were irradiated in clinical and near-clinical conditions at the Cyclotron Centre of the Bronowice proton therapy centre. PET data were acquired with a planar 16x16 cm$^2$ PET system. MC simulations of particle interactions and photon propagation in the phantoms were performed using the FLUKA code. The analysis included a comparison between experimental data and MC simulations of space and time profiles, as well as a fitting procedure to obtain the various isotope contributions in the phantoms. Results and conclusions: There was a good agreement between data and MC predictions in 1-dimensional space and decay rate distributions. The fractions of $^{11}$C, $^{15}$O and $^{10}$C that were obtained by fitting the decay rates with multiple simple exponentials generally agreed well with the MC expectations. We found a small excess of $^{10}$C in data compared to what was predicted in MC, which was clear especially in the PE phantom.
PET imaging is a non-invasive technique for particle range verification in proton therapy. It is based on measuring the beta+ annihilations caused by nuclear interactions of the protons in the patient. In this work we present measurements for proton range verification in phantoms, performed at the CNAO particle therapy treatment center in Pavia, Italy, with our 10 x 10 cm^2 planar PET prototype DoPET. PMMA phantoms were irradiated with mono-energetic proton beams and clinical treatment plans, and PET data were acquired during and shortly after proton irradiation. We created 1-D profiles of the beta+ activity along the proton beam-axis, and evaluated the difference between the proximal rise and the distal fall-off position of the activity distribution. A good agreement with FLUKA Monte Carlo predictions was obtained. We also assessed the system response when the PMMA phantom contained an air cavity. The system was able to detect these cavities quickly after irradiation.
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.
A broad-beam-delivery system for heavy-charged-particle radiotherapy often employs multiple collimators and a range-compensating filter, which potentially offer complex beam customization. In treatment planning, it is however difficult for a conventional pencil-beam algorithm to deal with these structures due to beam-size growth during transport. This study aims to resolve the problem with a novel computational model. The pencil beams are initially defined at the range compensating filter with angular-acceptance correction for the upstream collimators followed by the range compensation effects. They are individually transported with possible splitting near the downstream collimator edges to deal with its fine structure. The dose distribution for a carbon-ion beam was calculated and compared with existing experimental data. The penumbra sizes of various collimator edges agreed between them to a submillimeter level. This beam-customization model will complete an accurate and efficient dose-calculation algorithm for treatment planning with heavy charged particles.
A model for beam customization with collimators and a range-compensating filter based on the phase-space theory for beam transport is presented for dose distribution calculation in treatment planning of radiotherapy with protons and heavier ions. Independent handling of pencil beams in conventional pencil-beam algorithms causes unphysical collimator-height dependence in the middle of large fields, which is resolved by the framework comprised of generation, transport, collimation, regeneration, range-compensation, and edge-sharpening processes with a matrix of pencil beams. The model was verified to be consistent with measurement and analytic estimation at a submillimeter level in penumbra of individual collimators with a combinational-collimated carbon-ion beam. The model computation is fast, accurate, and readily applicable to pencil-beam algorithms in treatment planning with capability of combinational collimation to make best use of the beam-customization devices.
Measurements performed with the purpose of characterizing the charged secondary radiation for dose release monitoring in particle therapy are reported. Charged secondary yields, energy spectra and emission profiles produced in poly-methyl methacrylate (PMMA) target by 4He and 12C beams of different therapeutic energies were measured at 60 and 90 degree with respect to the primary beam direction. The secondary yields of protons produced along the primary beam path in PMMA target were obtained. The energy spectra of charged secondaries were obtained from time-of-flight information, whereas the emission profiles were reconstructed exploiting tracking detector information. The measured charged secondary yields and emission profiles are in agreement with the results reported in literature and confirm the feasibility of ion beam therapy range monitoring using 12C ion beam. The feasibility of range monitoring using charged secondary particles is also suggested for 4He ion beam.