Evaluation of a glassless photographic film scanner for high-gradient radiochromic film dosimetry.

This study evaluates the performance of the Nikon Coolscan 9000 ED film scanner for high-gradient radiochromic film dosimetry. As a reference for comparison, analogue experiments were performed on the Epson Expression 10000XL flatbed scanner. Based on these results, a dosimetric protocol was established for the Nikon scanner and its overall performance for high-gradient dosimetry was evaluated. The Nikon scanner demonstrated a high sensitivity for radiochromic film dosimetry, resulting in more contrast in the digitized image. The scanner's optics also demonstrated excellent stability and did not necessitate warm-up scans prior to data acquisition. Moreover, negative effects of temperature changes of the film inside the scanner were shown to be limited. None of the digitized images showed significant disturbances by moiré-patterns, by virtue of the absence of a glass plate for film positioning. However, scanner response was found to vary considerably across the reading area, requiring an optical density-dependent correction procedure to be incorporated into the scanning protocol. The main limitation of the Nikon Coolscan 9000 ED transmission scanner remains its film size restriction to 6.2 × 20 cm2. Nevertheless, its excellent characteristics render it the preferential tool for high-gradient radiochromic film dosimetry in applications limited to small film sizes, such as dosimetry in the build-up region.

The influence of small field sizes, penumbra, spot size and measurement depth on perturbation factors for microionization chambers.

The purpose of this study was the investigation of perturbation factors for microionization chambers in small field dosimetry and the influence of penumbra for different spot sizes. To this purpose, correlated sampling was implemented in the EGSnrc Monte Carlo (MC) user code cavity: CScavity. CScavity was first benchmarked against results in the literature for an NE2571 chamber. An efficiency increase of 17 was attained for the calculation of a realistic chamber perturbation factor in a water phantom. Calculations have been performed for microionization chambers of type PinPoint 31006 and PinPoint 31016 in full BEAMnrc linac simulations. Investigating the physical backgrounds of the differences for these small field settings, perturbation factors have been split up into (1) central electrode perturbation, (2) wall perturbation, (3) air-to-water perturbation (chamber volume air-to-water) and (4) water volume perturbation (water chamber volume to 1 mm(3) voxel). The influence of different spot sizes, position in penumbra, measuring depth and detector geometry on these perturbation factors has been investigated, in a 0.8 x 0.8 cm(2) field setting. p(cel) for the PP31006 steel electrode shows a variation of up to 1% in the lateral position, but only 0.4% for the PP31016 with an Al electrode. The air-to-water perturbation in the optimal scanning direction for both profiles and depth is most influenced by the radiation field, and only to a small extent the chamber geometry. The PP31016 geometry (shorter, larger radius) requires less total perturbation within the central axis of the field, but results in slightly larger variations off axis in the optimal scanning direction. Smaller spot sizes (0.6 mm FWHM) and sharper penumbras, compared to larger spot sizes (2 mm FWHM), result in larger perturbation starting in the penumbra. The longer geometries of the PP31006/14/15 exhibit in the non-optimal scanning direction large variations in total perturbation (p(tot) 1.201(4) (0.6 mm spot, 3 mm off axis, type A MC uncertainty) to 0.803(4) (5 mm off axis)) mainly due to volume perturbation. Therefore in IMRT settings, when the detector is not always in the optimal scanning direction, the PP31016 geometry requires less extreme perturbation (max p(tot) 1.130(3)) and shows less variation. However, these results suggest that small variations in positioning, spot size or MLC result in large differences in perturbation factors. Therefore even these 0.016 cm(3) ionization chambers are limited in their use for a field setting of 0.8 x 0.8 cm(2), as used in this investigation.

Development of in vitro models for investigating spatially fractionated irradiation: physics and biological results.

We present different in vitro experimental models which allow us to evaluate the effect of spatially fractionated dose distributions on metabolic activity. We irradiated a monolayer of MCF-7/6 human breast cancer cells with a steep and a smooth 6 MV x-ray dose gradient. In the steep gradient model, we irradiated the cells with three separate small fields. We also developed two smooth gradient models. In the first model, the cells are cultured in a T25 flask and irradiated with a smooth dose gradient over the length of the flask, while in the second one, the cells are cultured in a 96-well plate and also irradiated over the length of the plate. In an attempt to correlate the spatially fractionated dose distributions with metabolic activity, the effect of irradiation was evaluated by means of the MTT assay. This assay is used to determine the metabolic activity by measuring the amount of formazan formed after the conversion of MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) by cellular dehydrogenases. The results obtained with our different models suggest a dose-specific effect on metabolic activity, characterized by an increased formazan optical density occurring in the dose range 1.0-4.0 Gy in the steep dose gradient model and in the dose ranges 4.2-6.5 Gy and 2.3-5.1 Gy in the two smooth dose gradient models. The corresponding times for maximal formazan accumulation were 5-7 days in the steep dose gradient model and day 9-13 and day 9-11 in the smooth dose gradient models. Altogether, our results suggest that the MTT assay may be used as a biological dose-response meter to monitor the radiotherapeutic effectiveness.

On the calibration process of film dosimetry: OLS inverse regression versus WLS inverse prediction.

The purpose of this study was both putting forward a statistically correct model for film calibration and the optimization of this process. A reliable calibration is needed in order to perform accurate reference dosimetry with radiographic (Gafchromic) film. Sometimes, an ordinary least squares simple linear (in the parameters) regression is applied to the dose-optical-density (OD) curve with the dose as a function of OD (inverse regression) or sometimes OD as a function of dose (inverse prediction). The application of a simple linear regression fit is an invalid method because heteroscedasticity of the data is not taken into account. This could lead to erroneous results originating from the calibration process itself and thus to a lower accuracy. In this work, we compare the ordinary least squares (OLS) inverse regression method with the correct weighted least squares (WLS) inverse prediction method to create calibration curves. We found that the OLS inverse regression method could lead to a prediction bias of up to 7.3 cGy at 300 cGy and total prediction errors of 3% or more for Gafchromic EBT film. Application of the WLS inverse prediction method resulted in a maximum prediction bias of 1.4 cGy and total prediction errors below 2% in a 0-400 cGy range. We developed a Monte-Carlo-based process to optimize calibrations, depending on the needs of the experiment. This type of thorough analysis can lead to a higher accuracy for film dosimetry.

Monte Carlo modeling of the ModuLeaf miniature MLC for small field dosimetry and quality assurance of the clinical treatment planning system.

The purpose of this investigation was the verification of both the measured data and quality of the implementation of the add-on ModuLeaf miniature multileaf collimator (ML mMLC) into the clinical treatment planning system for conformal stereotactic radiosurgery treatment. To this end the treatment head with ML mMLC was modeled in the BEAMnrc Monte Carlo (MC) code. The 6 MV photon beams used in the setup were first benchmarked with a set of measurements. A total ML mMLC transmission of 1.13% of the 10 x 10 cm2 open field dose was measured and reproduced with the BEAMnrc/DOSXYZnrc code. Correspondence between calculated and measured output factors (OFs) was within 2%. Correspondence between MC and measured profiles was within 2% dose and 2 mm distance, only for the smallest 0.5 x 0.5 cm2 field the results were within 3% dose. In the next step, the MC model was compared with Gafchromic film measurements and Pinnacle(3) 7.4 f (convolution superposition algorithm) calculated dose distributions, using a gamma evaluation comparison, for a multi-beam patient setup delivered to a Lucytrade mark phantom. The gamma evaluation of the MC versus Gafchromic film resulted in 3.4% of points not fulfilling gamma

Precautions and strategies in using a commercial flatbed scanner for radiochromic film dosimetry.

The purpose of this study was to investigate the value of a commercially available flatbed scanner for film dosimetry with radiochromic film for external radiotherapy. The EPSON Pro 1680 Expression scanner was examined as a densitometer for two-dimensional film dosimetry with Gafchromic EBT film. An accurate and efficient scanning procedure was established. Possible drift and warm-up effects of the scanner were studied and the direct physical influence of the scanner light on the radiochromic film was assessed. Next, we investigated the scan field uniformity. Also, we examined if the accuracy of radiochromic film was improved by subtracting the optical density of the unirradiated blank film from the optical density of the irradiated film. To assess the accuracy of Gafchromic EBT film when the EPSON scanner was used as a densitometer, the depth dose of a 2 x 15 cm(2) field and the in-plane and cross-plane profiles of a 15 x 15 cm(2) field were measured and compared with diamond detector measurements. When taking consecutive scans, we found that the optical density taken from the first scan was about 1% higher than the optical density taken from subsequent scans. We attribute this to the warming up of the lamp of the scanner. Longer-term drift of the scanner was found to be absent. We found that the use of a correction matrix was necessary to correct for the non-uniform scanner response over the scan field. Subtracting the optical density of the unirradiated blank film from the irradiated film improves the precision of the Gafchromic EBT film. Depth dose and profile measurements with Gafchromic EBT film and the diamond detector are in agreement within 2.5%. The EPSON Pro 1680 Expression scanner is an excellent tool for accurate two-dimensional film dosimetry with Gafchromic EBT film provided that some precautions and corrections are taken into account.

Decoupling initial electron beam parameters for Monte Carlo photon beam modelling by removing beam-modifying filters from the beam path.

A new method is presented to decouple the parameters of the incident e(-) beam hitting the target of the linear accelerator, which consists essentially in optimizing the agreement between measurements and calculations when the difference filter, which is an additional filter inserted in the linac head to obtain uniform lateral dose-profile curves for the high energy photon beam, and flattening filter are removed from the beam path. This leads to lateral dose-profile curves, which depend only on the mean energy of the incident electron beam, since the effect of the radial intensity distribution of the incident e- beam is negligible when both filters are absent. The location of the primary collimator and the thickness and density of the target are not considered as adjustable parameters, since a satisfactory working Monte Carlo model is obtained for the low energy photon beam (6 MV) of the linac using the same target and primary collimator. This method was applied to conclude that the mean energy of the incident e- beam for the high energy photon beam (18 MV) of our Elekta SLi Plus linac is equal to 14.9 MeV. After optimizing the mean energy, the modelling of the filters, in accordance with the information provided by the manufacturer, can be verified by positioning only one filter in the linac head while the other is removed. It is also demonstrated that the parameter setting for Bremsstrahlung angular sampling in BEAMnrc ('Simple' using the leading term of the Koch and Motz equation or 'KM' using the full equation) leads to different dose-profile curves for the same incident electron energy for the studied 18 MV beam. It is therefore important to perform the calculations in 'KM' mode. Note that both filters are not physically removed from the linac head. All filters remain present in the linac head and are only rotated out of the beam. This makes the described method applicable for practical usage since no recommissioning process is required.

Experimental verification of lung dose with radiochromic film: comparison with Monte Carlo simulations and commercially available treatment planning systems.

The purpose of this study was to assess the absorbed dose in and around lung tissue by performing radiochromic film measurements, Monte Carlo simulations and calculations with superposition convolution algorithms. We considered a layered polystyrene phantom of 12 x 12 x 12 cm3 containing a central cavity of 6 x 6 x 6 cm3 filled with Gammex RMI lung-equivalent material. Two field configurations were investigated, a small 1 x 10 cm2 field and a larger 10 x 10 cm2 field. First, we performed Monte Carlo simulations to investigate the influence of radiochromic film itself on the measured dose distribution when the film intersects a lung-equivalent region and is oriented parallel to the central beam axis. To that end, the film and the lung-equivalent materials were modelled in detail, taking into account their specific composition. Next, measurements were performed with the film oriented both parallel and perpendicular to the central beam axis to verify the results of our Monte Carlo simulations. Finally, we digitized the phantom in two commercially available treatment planning systems, Helax-TMS version 6.1A and Pinnacle version 6.2b, and calculated the absorbed dose in the phantom with their incorporated superposition convolution algorithms to compare with the Monte Carlo simulations. Comparing Monte Carlo simulations with measurements reveals that radiochromic film is a reliable dosimeter in and around lung-equivalent regions when the film is positioned perpendicular to the central beam axis. Radiochromic film is also able to predict the absorbed dose accurately when the film is positioned parallel to the central beam axis through the lung-equivalent region. However, attention must be paid when the film is not positioned along the central beam axis, in which case the film gradually attenuates the beam and decreases the dose measured behind the cavity. This underdosage disappears by offsetting the film a few centimetres. We find deviations of about 3.6% between Monte Carlo and the superposition convolution algorithm of Pinnacle behind the lung region, for both field configurations. Pinnacle is quite accurate in the lung region. Deviations up to 5.6% for the small field are found in the lung region between Monte Carlo and the superposition convolution algorithm of Helax-TMS. Behind the lung region, Helax-TMS is in better agreement with Monte Carlo. Radiochromic film measurements or Monte Carlo simulations are reliable methods to establish the dose in and around lung tissue.

The importance of accurate linear accelerator head modelling for IMRT Monte Carlo calculations.

Two Monte Carlo dose engines for radiotherapy treatment planning, namely a beta release of Peregrine and MCDE (Monte Carlo dose engine), were compared with Helax-TMS (collapsed cone superposition convolution) for a head and neck patient for the Elekta SLi plus linear accelerator. Deviations between the beta release of Peregrine and MCDE up to 10% were obtained in the dose volume histogram of the optical chiasm. It was illustrated that the differences are not caused by the particle transport in the patient, but by the modelling of the Elekta SLi plus accelerator head and more specifically the multileaf collimator (MLC). In MCDE two MLC modules (MLCQ and MLCE) were introduced to study the influence of the tongue-and-groove geometry, leaf bank tilt and leakage on the actual dose volume histograms. Differences in integral dose in the optical chiasm up to 3% between the two modules have been obtained. For single small offset beams though the FWHM of lateral profiles obtained with MLCE can differ by more than 1.5 mm from profiles obtained with MLCQ. Therefore, and because the recent version of MLCE is as fast as MLCQ, we advise to use MLCE for modelling the Elekta MLC. Nevertheless there still remains a large difference (up to 10%) between Peregrine and MCDE. By studying small offset beams we have shown that the profiles obtained with Peregrine are shifted, too wide and too flat compared with MCDE and phantom measurements. The overestimated integral doses for small beam segments explain the deviations observed in the dose volume histograms. The Helax-TMS results are in better agreement with MCDE, although deviations exceeding 5% have been observed in the optical chiasm. Monte Carlo dose deviations of more than 10% as found with Peregrine are unacceptable as an influence on the clinical outcome is possible and as the purpose of Monte Carlo treatment planning is to obtain an accuracy of 2%. We would like to emphasize that only the Elekta MLC has been tested in this work, so it is certainly possible that alpha releases of Peregrine provide more accurate results for other accelerators.

Comparison of build-up dose between Elekta and Varian linear accelerators for high-energy photon beams using radiochromic film and clinical implications for IMRT head and neck treatments.

Skin toxicity has been reported for IMRT of head and neck cancer. The purpose of this study was to investigate the dose in the build-up region delivered by a 6 MV treatment plan for which important skin toxicity was observed. We also investigated if the different designs of the treatment head of an Elekta and a Varian linear accelerator, especially the lower position of the Varian multi-leaf collimator, give rise to different build-up doses. For regular square open beams, the build-up dose along the central beam axis is higher for the Varian machine than for the Elekta machine, both for 6 MV and 18 MV. At the Elekta machine at 18 MV, the superficial dose of a diamond shaped 10 x 10 cm2 field is 3.6% lower than the superficial dose of a regular 10 x 10 cm2 field. This effect is not seen at 6 MV. At the Varian machine, the superficial dose of the diamond shaped field is respectively 3.5 and 14.2% higher than the superficial dose of the regular 10 x 10 cm2 field for 6 MV and 18 MV. Despite the differences measured in build-up dose for single beams between the Elekta and the Varian linear accelerator, there were no measurable differences in superficial dose when a typical IMRT dose plan of 6 MV for a head and neck tumour is executed at the two machines.

MCDE: a new Monte Carlo dose engine for IMRT.

A new accurate Monte Carlo code for IMRT dose computations, MCDE (Monte Carlo dose engine), is introduced. MCDE is based on BEAMnrc/DOSXYZnrc and consequently the accurate EGSnrc electron transport. DOSXYZnrc is reprogrammed as a component module for BEAMnrc. In this way both codes are interconnected elegantly, while maintaining the BEAM structure and only minimal changes to BEAMnrc.mortran are necessary. The treatment head of the Elekta SLiplus linear accelerator is modelled in detail. CT grids consisting of up to 200 slices of 512 x 512 voxels can be introduced and up to 100 beams can be handled simultaneously. The beams and CT data are imported from the treatment planning system GRATIS via a DICOM interface. To enable the handling of up to 50 x 10(6) voxels the system was programmed in Fortran95 to enable dynamic memory management. All region-dependent arrays (dose, statistics, transport arrays) were redefined. A scoring grid was introduced and superimposed on the geometry grid, to be able to limit the number of scoring voxels. The whole system uses approximately 200 MB of RAM and runs on a PC cluster consisting of 38 1.0 GHz processors. A set of in-house made scripts handle the parallellization and the centralization of the Monte Carlo calculations on a server. As an illustration of MCDE, a clinical example is discussed and compared with collapsed cone convolution calculations. At present, the system is still rather slow and is intended to be a tool for reliable verification of IMRT treatment planning in the case of the presence of tissue inhomogeneities such as air cavities.

The value of radiochromic film dosimetry around air cavities: experimental results and Monte Carlo simulations.

In this study we investigate radiochromic film dosimetry around air cavities with particular focus on the perturbation of the dose distribution by the film when the film is parallel to the beam axis. We considered a layered polystyrene phantom containing an air cavity as a model for the air-soft tissue geometry that may occur after surgical resection of a paranasal sinus tumour. A radiochromic film type MD-55 was positioned within the phantom so that it intersected the cavity. Two phantom set-ups were examined. In the first case, the air cavity is at the centre of the phantom, thus the film is lying along the central beam axis. In the second case, the cavity and film are located 2 cm offset from the phantom centre and the central beam axis. In order to examine the influence of the film on the dose distribution and to interpret the film-measured results, Monte Carlo simulations were performed. The film was modelled rigorously to incorporate the composition and structure of the film. Two field configurations, a 1 x 10 cm2 field and a 10 x 10 cm2 field, were examined. The dose behind the air cavity is reduced by 6 to 7% for both field configurations when a film that intersects the cavity contains the central beam axis. This is due to the attenuation exerted by the film when photons cross the cavity. Offsetting the beam to the cavity and the film by 2 cm removes the dose reduction behind the air cavity completely. Another result was that the rebuild-up behind the cavity for the 10 x 10 cm2 field, albeit less significant than for the 1 x 10 cm2 field, could only be measured by the film that was placed offset with respect to the central beam axis. Although radiochromic film is approximately soft-tissue equivalent and energy independent as compared to radiographic films, care should be taken in the case of inhomogeneous phantoms when the film intersects air cavities and contains the beam central axis. Errors in dose measurement can be expected distal to the air cavity due to attenuation in the film itself. This attenuation would not occur in the absence of the film. Both experiments and Monte Carlo computations support this conclusion.