Experimental and Monte Carlo studies of fluence corrections for graphite calorimetry in low- and high-energy clinical proton beams

Purpose: The aim of this study was to determine fluence corrections necessary to convert absorbed dose to graphite, measured by graphite calorimetry, to absorbed dose to water. Fluence corrections were obtained from experiments and Monte Carlo simulations in low- and high-energy proton beams. Methods: Fluence corrections were calculated to account for the difference in fluence between water and graphite at equivalent depths. Measurements were performed with narrow proton beams. Plane-parallel-plate ionization chambers with a large collecting area compared to the beam diameter were used to intercept the whole beam. High- and low-energy proton beams were provided by a scanning and double scattering delivery system, respectively. A mathematical formalism was established to relate fluence corrections derived from Monte Carlo simulations, using the fluka code [A. Ferrari et al., “fluka: A multi-particle transport code,” in CERN 2005-10, INFN/TC 05/11, SLAC-R-773 (2005) and T. T. Böhlen et al., “The fluka Code: Developments and challenges for high energy and medical applications,” Nucl. Data Sheets 120, 211–214 (2014)], to partial fluence corrections measured experimentally. Results: A good agreement was found between the partial fluence corrections derived by Monte Carlo simulations and those determined experimentally. For a high-energy beam of 180 MeV, the fluence corrections from Monte Carlo simulations were found to increase from 0.99 to 1.04 with depth. In the case of a low-energy beam of 60 MeV, the magnitude of fluence corrections was approximately 0.99 at all depths when calculated in the sensitive area of the chamber used in the experiments. Fluence correction calculations were also performed for a larger area and found to increase from 0.99 at the surface to 1.01 at greater depths. Conclusions: Fluence corrections obtained experimentally are partial fluence corrections because they account for differences in the primary and part of the secondary particle fluence. A correction factor, F(d), has been established to relate fluence corrections defined theoretically to partial fluence corrections derived experimentally. The findings presented here are also relevant to water and tissue-equivalent-plastic materials given their carbon content

Time-resolved dosimetry for validation of 4D dose calculation in PBS proton therapy

Four-dimensional dose calculation (4D-DC) is crucial for predicting the dosimetric outcome in the presence of intra-fractional organ motion. Time-resolved dosimetry can provide significant insights into 4D pencil beam scanning dose accumulation and is therefore irreplaceable for benchmarking 4D-DC. In this study a novel approach of time-resolved dosimetry using five PinPoint ionization chambers (ICs) embedded in an anthropomorphic dynamic phantom was employed and validated against beam delivery details. Beam intensity variations as well as the beam delivery time structure were well reflected with an accuracy comparable to the temporal resolution of the IC measurements. The 4D dosimetry approach was further applied for benchmarking the 4D-DC implemented in the RayStation 6.99 treatment planning system. Agreement between computed values and measurements was investigated for (i) partial doses based on individual breathing phases, and (ii) temporally distributed cumulative doses. For varied beam delivery and patient-related parameters the average unsigned dose difference for (i) was 0.04 +/- 0.03 Gy over all considered IC measurement values, while the prescribed physical dose was 2 Gy. By implementing (ii), a strong effect of the dose gradient on measurement accuracy was observed. The gradient originated from scanned beam energy modulation and target motion transversal to the beam. Excluding measurements in the high gradient the relative dose difference between measurements and 4D-DCs for a given treatment plan at the end of delivery was 3.5% on average and 6.6% at maximum over measurement points inside the target. Overall, the agreement between 4D dose measurements in the moving phantom and retrospective 4D-DC was found to be comparable to the static dose differences for all delivery scenarios. The presented 4D-DC has been proven to be suitable for simulating treatment deliveries with various beam- as well as patient-specific parameters and can therefore be employed for dosimetric validation of different motion mitigation techniques

Determination of beam quality correction factors for the Roos plane-parallel ionisation chamber exposed to very high energy electron (VHEE) beams using Geant4

Detailed characterisation of the Roos secondary standard plane-parallel ionisation chamber has been conducted in a novel 200 MeV Very High Energy Electron (VHEE) beam with reference to the standard 12 MeV electron calibration beam used in our experimental work. Stopping-power-ratios and perturbation factors have been determined for both beams and used to calculated the beam quality correction factor using the Geant4 general purpose MC code. These factors have been calculated for a variety of charged particle transport parameters available in Geant4 which were found to pass the Fano cavity test. Stopping-power-ratios for the 12 MeV electron calibration beam quality were found to agree within uncertainties to that quoted by current dosimetry protocols. Perturbation factors were found to vary by up-to 4% for the calibration beam depending on the parameter configuration, compared with only 0.8% for the VHEE beam. Beam quality correction factors were found to describe an approximately 10% lower dose than would be originally calculated if a beam quality correction were not accounted for. Moreover, results presented here largely resolve unphysical chamber measurements, such as collection efficiencies greater than 100%, and assist in the accurate determination of absorbed dose and ion recombination in secondary standard ionisation chambers

Development of optimised tissue-equivalent materials for proton therapy

OBJECTIVE: In proton therapy there is a need for proton optimised tissue-equivalent materials as existing phantom materials can produce large uncertainties in the determination of absorbed dose and range measurements. The aim of this work is to develop and characterise optimised tissue-equivalent materials for proton therapy. APPROACH: A mathematical model was developed to enable the formulation of epoxy-resin based tissue-equivalent materials that are optimised for all relevant interactions of protons with matter, as well as photon interactions, which play a role in the acquisition of CT numbers. This model developed formulations for vertebra bone- and skeletal muscle-equivalent plastic materials. The tissue equivalence of these new materials and commercial bone- and muscle-equivalent plastic materials were theoretical compared against biological tissue compositions. The new materials were manufactured and characterised by their mass density, relative stopping power (RSP) measurements, and CT scans to evaluate their tissue-equivalence. MAIN RESULTS: Results showed that existing tissue-equivalent materials can produce large uncertainties in proton therapy dosimetry. In particular commercial bone materials showed to have a relative difference up to 8 % for range. On the contrary, the best optimised formulations were shown to mimic their target human tissues within 1-2 % for the mass density and RSP. Furthermore, their CT-predicted RSP agreed within 1-2 % of the experimental RSP, confirming their suitability as clinical phantom materials. SIGNIFICANCE: We have developed a tool for the formulation of tissue-equivalent materials optimised for proton dosimetry. Our model has enabled the development of proton optimised tissue-equivalent materials which perform better than existing tissue-equivalent materials. These new materials will enable the advancement of clinical proton phantoms for accurate proton dosimetry

A systematic Monte Carlo study of secondary electron fluence perturbation in clinical proton beams (70-250 MeV) for cylindrical and spherical ion chambers.

A systematic Monte Carlo study of secondary electron fluence perturbation in clinical proton beams (70–250 MeV) for cylindrical and spherical ion chambers. Current dosimetry protocols for clinical protons do not take into account any secondary electron fluence perturbation in ion chambers. In this work we performed a systematic study of secondary electron fluence perturbation factors for spherical and cylindrical ion chambers in proton beams (70–250 MeV). The electron fluence perturbation factor pe was calculated using Monte Carlo transport of protons and secondary electrons. The influence of proton energy cavity wall material (graphite water A150 PMMA polystyrene) cavity radius cavity wall thickness and positioning depth in water is studied. The influence of inelastic nuclear proton interactions is briefly discussed. It was found that pe depends on wall material; the largest values for pe were obtained for ion chambers with A150 walls (pe = 1.009) the smallest values for graphite walls. The perturbation factor was found to be largely independent of proton energy. A slight decrease of pe with cavity radius was obtained especially for low energy protons. The wall thickness was found to have no effect on pe in the range studied (0.025–0.1 cm). The depth of the cavity in a water phantom was also found to have an insignificant effect on pe. Based on the results in the paper for spherical and cylindrical ion chambers a method to calculate pe for a thimble ion chamber is presented. The results presented in this paper for cylindrical and spherical ion chambers are in contradiction to the calculated electron fluence perturbation factors for planar ion chambers in the paper by Casnati et al