Hugo Bouchard

GPUMCD: A new GPU‐oriented Monte Carlo dose calculation platform

PURPOSE Monte Carlo methods are considered as the gold standard for dosimetric computations in radiotherapy. Their execution time is, however, still an obstacle to the routine use of Monte Carlo packages in a clinical setting. To address this problem, a completely new, and designed from the ground up for the GPU, Monte Carlo dose calculation package for voxelized geometries is proposed: GPUMCD. METHOD GPUMCD implements a coupled photon-electron Monte Carlo simulation for energies in the range of 0.01-20 MeV. An analog simulation of photon interactions is used and a class II condensed history method has been implemented for the simulation of electrons. A new GPU random number generator, some divergence reduction methods, as well as other optimization strategies are also described. GPUMCD was run on a NVIDIA GTX480, while single threaded implementations of EGSnrc and DPM were run on an Intel Core i7 860. RESULTS Dosimetric results obtained with GPUMCD were compared to EGSnrc. In all but one test case, 98% or more of all significant voxels passed the gamma criteria of 2%-2 mm. In terms of execution speed and efficiency, GPUMCD is more than 900 times faster than EGSnrc and more than 200 times faster than DPM, a Monte Carlo package aiming fast executions. Absolute execution times of less than 0.3 s are found for the simulation of 1M electrons and 4M photons in water for monoenergetic beams of 15 MeV, including GPU-CPU memory transfers. CONCLUSION GPUMCD, a new GPU-oriented Monte Carlo dose calculation platform, has been compared to EGSnrc and DPM in terms of dosimetric results and execution speed. Its accuracy and speed make it an interesting solution for full Monte Carlo dose calculation in radiation oncology.

On the characterization and uncertainty analysis of radiochromic film dosimetry.

Radiochromic film is a dosimeter of choice in applications requiring high spatial resolution, two dimensional measurements, or minimum perturbation of the beam fluence. Since the measurement uncertainty in Gafchromic film dosimetry is thought to be significant compared to that of ionization chambers, a rigorous method to evaluate measurement uncertainties is desired. This article provides a method that takes into account the correlation between fit parameters as well as single dose values in order to obtain accurate uncertainties in absolute and relative measurements. A complete portrait of all sources of uncertainty in Gafchromic film dosimetry is given. The parametrization of variance as a function of the number of averaged pixels is obtained in order to accurately predict the uncertainty as a function of the size of the region of interest. The choice of functional form for the sensitometric curve is based on four criteria and a convergence of global net optical density uncertainty to 0.0013 is demonstrated. A minimum number of 12 points is recommended to characterize the sensitometric curve to a sufficient precision on the uncertainty estimation. Uncertainty levels of 0.9% on absolute dose measurements and 0.45% on relative measurements are achieved using a 12-point calibration curve with 220 cGy and repeating measurements five times. Uncertainties of 0.8% and 0.4% are achievable when using 35 points during film characterization. Ignoring covariance terms is shown to lead to errors in the estimation of uncertainty.

Detector dose response in megavoltage small photon beams. I. Theoretical concepts

PURPOSE To explain the reasons for significant quality correction factors in megavoltage small photon fields and clarify the underlying concepts relevant to dosimetry under such conditions. METHODS The validity of cavity theory and the requirement of charged particle equilibrium (CPE) are addressed from a theoretical point of view in the context of nonstandard beams. Perturbation effects are described into four main subeffects, explaining their nature and pointing out their relative importance in small photon fields. RESULTS It is demonstrated that the failure to meet classical cavity theory requirements, such as CPE, is not the reason for significant quality correction factors. On the contrary, it is shown that the lack of CPE alone cannot explain these corrections and that what matters most, apart from volume averaging effects, is the relationship between the lack of CPE in the small field itself and the density of the detector cavity. The density perturbation effect is explained based on Fanos theorem, describing the compensating effect of two main contributions to cavity absorbed dose. Using the same approach, perturbation effects arising from the difference in atomic properties of the cavity medium and the presence of extracameral components are explained. Volume averaging effects are also discussed in detail. CONCLUSIONS Quality correction factors of small megavoltage photon fields are mainly due to differences in electron density between water and the detector medium and to volume averaging over the detector cavity. Other effects, such as the presence of extracameral components and differences in atomic properties of the detection medium with respect to water, can also play an accentuated role in small photon fields compared to standard beams.

Reference dosimetry using radiochromic film

The objectives of this study are to identify and quantify factors that influence radiochromic film dose response and to determine whether such films are suitable for reference dosimetry. The influence of several parameters that may introduce systematic dose errors when performing reference dose measurements were investigated. The effect of the film storage temperature was determined by comparing the performance of three lots of GAFCHROMIC EBT2 films stored at either 4°C or room temperature. The effect of high (>80%) or low (<20%) relative humidity was also determined. Doses measured in optimal conditions with EBT and EBT2 films were then compared with an A12 ionization chamber measurement. Intensity‐modulated radiation therapy quality controls using EBT2 films were also performed in reference dose. The results obtained using reference dose measurements were compared with those obtained using relative dose measurements. Storing the film at 4°C improves the stability of the film over time, but does not eliminate the noncatalytic film development, seen as a rise in optical density over time in the absence of radiation. Relative humidity variations ranging from 80% to 20% have a strong impact on the optical density and could introduce dose errors of up to 15% if the humidity were not controlled during the film storage period. During the scanning procedure, the film temperature influences the optical density that is measured. When controlling for these three parameters, the dose differences between EBT or EBT2 and the A12 chamber are found to be within ±4% (2σ level) over a dose range of 20–350 cGy. Our results also demonstrate the limitation of the Anisotropic Analytical Algorithm for dose calculation of highly modulated treatment plans. PACS numbers: 87.55.Qr; 87.56.Fc

Correction-less dosimetry of nonstandard photon fields: a new criterion to determine the usability of radiation detectors.

In the IAEA-AAPM dosimetry formalism, detector measurements in general nonstandard conditions are corrected using the factor k(f(clin),f(msr))(Q(clin),Q(msr)). This factor needs to be evaluated on a case-by-case basis which is difficult to accomplish in practice. The present paper aims to provide a method that allows neglecting correction factors for small and composite IMRT fields by first determining a radiation detectors usability in these fields. Detailed models of nine radiation detectors are built: four ionization chambers (NE2571, A12, A1SL, A14), three small field detectors (PTW31018 microLion, PTW60003 natural diamond, PTW60012 unshielded diode) and two near water-equivalent detectors (alanine, W1 scintillating fiber). Using the egs_chamber Monte Carlo code, dose response functions at 6 MV and 25 MV are sampled for each detector and their corresponding volume of water. These functions are then used with a newly derived criterion to evaluate an upper bound ξ(f(ns),f(msr))(Q(ns),Q(msr)) on the variable ε(f(ns),f(msr))(Q(ns),Q(msr)) if no field collimation/modulation occurs over a given perturbation zone. The variable ε(f(ns),f(msr))(Q(ns),Q(msr)) is defined as the absolute value of the relative deviation from unity of a nonstandard field quality correction factor k(f(ns),f(msr))(Q(ns),Q(msr)). Using the same criterion, perturbation zones are evaluated by finding the smallest field size allowed for correction-less dosimetry with a given tolerance ξ(f(ns),f(msr))(Q(ns),Q(msr)). For composite fields, the sensitivity of detectors to the non-uniformity of virtual symmetric collapsed beams over regions of interest specified by the criterion is studied to estimate an upper bound ξ ̃(f(ns),f(ref))(Q(ns),Q) on ε(f(ns),f(ref))(Q(ns),Q) for a given beam flatness. Finally, a newly defined perturbation function is used to minimize the perturbations of the microLion chamber through density compensation. The theoretical criterion shows good agreement with full Monte Carlo simulations of ε(f(ns),f(msr))(Q(ns) Q(msr)). Perturbation zones are shown to be sensitive to both the energy of the beam and the orientation of the detector. The density-compensated microLion shows significantly improved response in both axial and radial orientations in small and composite IMRT fields. Finally, the new Exradin W1 scintillator is shown to have ξ(f(ns),f(msr))(Q(ns),Q(msr)) values under 1% in small fields. The methods presented in this work theoretically show that correction-less dosimetry of nonstandard field can be accomplished by knowing the limit of usability of radiation detectors in these conditions. Potential applications include small field output factor measurements and absolute absorbed dose to water verification in the QA of clinical IMRT fields.

The potential of dual‐energy CT to reduce proton beam range uncertainties

Purpose Dual‐energy CT (DECT) promises improvements in estimating stopping power ratios (SPRs) for proton therapy treatment planning. Although several comparable mathematical formalisms have been proposed in literature, the optimal techniques to characterize human tissue SPRs with DECT in a clinical environment are not fully established. The aim of this work is to compare the most robust DECT methods against conventional single‐energy CT (SECT) in conditions reproducing a clinical environment, where CT artifacts and noise play a major role on the accuracy of these techniques. Methods Available DECT tissue characterization methods are investigated and their ability to predict SPRs is compared in three contexts: (a) a theoretical environment using the XCOM cross section database; (b) experimental data using a dual‐source CT scanner on a calibration phantom; (c) simulations of a virtual humanoid phantom with the ImaSim software. The latter comparison accounts for uncertainties caused by CT artifacts and noise, but leaves aside other sources of uncertainties such as CT grid size and the I‐values. To evaluate the clinical impact, a beam range calculation model is used to predict errors from the probability distribution functions determined with ImaSim simulations. Range errors caused by SPR errors in soft tissues and bones are investigated. Results Range error estimations demonstrate that DECT has the potential of reducing proton beam range uncertainties by 0.4% in soft tissues using low noise levels of 12 and 8 HU in DECT, corresponding to 7 HU in SECT. For range uncertainties caused by the transport of protons through bones, the reduction in range uncertainties for the same levels of noise is found to be up to 0.6 to 1.1 mm for bone thicknesses ranging from 1 to 5 cm, respectively. We also show that for double the amount noise, i.e., 14 HU in SECT and 24 and 16 HU for DECT, the advantages of DECT in soft tissues are lost over SECT. In bones however, the reduction in range uncertainties is found to be between 0.5 and 0.9 mm for bone thicknesses ranging from 1 to 5 cm, respectively. Conclusion DECT has a clear potential to improve proton beam range predictions over SECT in proton therapy. However, in the current state high levels of noise remain problematic for DECT characterization methods and do not allow getting the full benefits of this technology. Future work should focus on adapting DECT methods to noise and investigate methods based on raw‐data to reduce CT artifacts.

Investigation of three radiation detectors for accurate measurement of absorbed dose in nonstandard fields

PURPOSE To establish accurate experimental dosimetry techniques for reference dose measurements in nonstandard composite fields. METHODS A cylindrical PMMA phantom filled with water was constructed, at the center of which reference absorbed dose to water for a head and neck IMRT delivery was measured. Based on the proposed new formalism for reference dosimetry of nonstandard fields [Alfonso et al., Med. Phys. 35, 5179-5186 (2008)], a candidate plan-class specific reference (pcsr) field for a typical head and neck IMRT delivery was created on the CT images of the phantom. The absorbed dose to water in the pcsr field normalized to that in a reference 10×10cm2 field was measured using three radiation detectors: Gafchromic® EBT films, a diamond detector, and a guarded liquid-filled ionization chamber developed in-house (GLIC-03). Pcsr correction factors kQpcsr,Qfpcsr,fref were determined for five different types of air-filled ionization chambers (Exradin A12, NE2571, Exradin A1SL, Exradin A14, and PinPoint® 31006) in a fully rotated delivery and in a delivery with the same MLC settings and weights but from a single gantry angle (a collapsed delivery). RESULTS The combined standard uncertainty in measuring the correction factor kQpcsr,Qfpcsr,fref using the three dosimetry techniques was 0.3%. For all the air-filled ionization chambers and the pcsr field tested, the correction factor was not different from unity by more than ±0.8%. For the fully rotated delivery, the correction factors were in a narrow range of 0.9955-0.9986, while in the collapsed delivery, they were in a slightly broader range of 0.9922-1.0048. In the collapsed delivery, the Farmer-type chambers (Exradin A12 and NE2571) had very similar correction factors (0.9922 and 0.9931, respectively), whereas the correction factors for the smaller chambers showed more distinct chamber-type dependence. CONCLUSIONS The authors have established three experimental dosimetry techniques that allow reference measurements of nonstandard field correction factors kQpcsr,Qfpcsr,fref for air-filled ionization chambers at the 0.3% 1σ uncertainty level. These techniques can be used to determine criteria for the selection of plan-class specific reference fields and ultimately improve clinical reference dosimetry of nonstandard fields.

Detector dose response in megavoltage small photon beams. II. Pencil beam perturbation effects

PURPOSE To quantify detector perturbation effects in megavoltage small photon fields and support the theoretical explanation on the nature of quality correction factors in these conditions. METHODS In this second paper, a modern approach to radiation dosimetry is defined for any detector and applied to small photon fields. Fanos theorem is adapted in the form of a cavity theory and applied in the context of nonstandard beams to express four main effects in the form of perturbation factors. The pencil-beam decomposition method is detailed and adapted to the calculation of perturbation factors and quality correction factors. The approach defines a perturbation function which, for a given field size or beam modulation, entirely determines these dosimetric factors. Monte Carlo calculations are performed in different cavity sizes for different detection materials, electron densities, and extracameral components. RESULTS Perturbation effects are detailed with calculated perturbation functions, showing the relative magnitude of the effects as well as the geometrical extent to which collimating or modulating the beam impacts the dosimetric factors. The existence of a perturbation zone around the detector cavity is demonstrated and the approach is discussed and linked to previous approaches in the literature to determine critical field sizes. CONCLUSIONS Monte Carlo simulations are valuable to describe pencil beam perturbation effects and detail the nature of dosimetric factors in megavoltage small photon fields. In practice, it is shown that dosimetric factors could be avoided if the field size remains larger than the detector perturbation zone. However, given a detector and beam quality, a full account for the detector geometry is necessary to determine critical field sizes.

Backscatter factors and mass energy-absorption coefficient ratios for diagnostic radiology dosimetry

Backscatter factors, B, and mass energy-absorption coefficient ratios, (μ(en)/ρ)(w, air), for the determination of the surface dose in diagnostic radiology were calculated using Monte Carlo simulations. The main purpose was to extend the range of available data to qualities used in modern x-ray techniques, particularly for interventional radiology. A comprehensive database for mono-energetic photons between 4 and 150 keV and different field sizes was created for a 15 cm thick water phantom. Backscattered spectra were calculated with the PENELOPE Monte Carlo system, scoring track-length fluence differential in energy with negligible statistical uncertainty; using the Monte Carlo computed spectra, B factors and (μ(en)/ρ)(w, air) were then calculated numerically for each energy. Weighted averaging procedures were subsequently used to convolve incident clinical spectra with mono-energetic data. The method was benchmarked against full Monte Carlo calculations of incident clinical spectra obtaining differences within 0.3-0.6%. The technique used enables the calculation of B and (μ(en)/ρ)(w, air) for any incident spectrum without further time-consuming Monte Carlo simulations. The adequacy of the extended dosimetry data to a broader range of clinical qualities than those currently available, while keeping consistency with existing data, was confirmed through detailed comparisons. Mono-energetic and spectra-averaged values were compared with published data, including those in ICRU Report 74 and IAEA TRS-457, finding average differences of 0.6%. Results are provided in comprehensive tables appropriated for clinical use. Additional qualities can easily be calculated using a designed GUI interface in conjunction with software to generate incident photon spectra.

A Monte Carlo method to evaluate the impact of positioning errors on detector response and quality correction factors in nonstandard beams

During experimental procedures, an adequate evaluation of all sources of uncertainty is necessary to obtain an overall uncertainty budget. In specific radiation dosimetry applications where a single detector is used, common methods to evaluate uncertainties caused by setup positioning errors are not applicable when the dose gradient is not known a priori. This study describes a method to compute these uncertainties using the Monte Carlo method. A mathematical formalism is developed to calculate unbiased estimates of the uncertainties. The method is implemented in egs_chamber, an EGSnrc-based code that allows for the efficient calculation of detector doses and dose ratios. The correct implementation of the method into the egs_chamber code is validated with an extensive series of tests. The accuracy of the developed mathematical formalism is verified by comparing egs_chamber simulation results to the theoretical expectation in an ideal situation where the uncertainty can be computed analytically. Three examples of uncertainties are considered for realistic models of an Exradin A12 ionization chamber and a PTW 60012 diode, and results are computed for parameters representing nearly realistic positioning error probability distributions. Results of practical examples show that uncertainties caused by positioning errors can be significant during IMRT reference dosimetry as well as small field output factor measurements. The method described in this paper is of interest in the study of single-detector response uncertainties during nonstandard beam measurements, both in the scope of daily routine as well as when developing new dosimetry protocols. It is pointed out that such uncertainties should be considered in new protocols devoted to single-detector measurements in regions with unpredictable dose gradients. The method is available within the egs_chamber code in the latest official release of the EGSnrc system.