F. Sánchez-Doblado

Ionization chamber dosimetry of small photon fields: a Monte Carlo study on stopping-power ratios for radiosurgery and IMRT beams

Absolute dosimetry with ionization chambers of the narrow photon fields used in stereotactic techniques and IMRT beamlets is constrained by lack of electron equilibrium in the radiation field. It is questionable that stopping-power ratio in dosimetry protocols, obtained for broad photon beams and quasi-electron equilibrium conditions, can be used in the dosimetry of narrow fields while keeping the uncertainty at the same level as for the broad beams used in accelerator calibrations. Monte Carlo simulations have been performed for two 6 MV clinical accelerators (Elekta SL-18 and Siemens Mevatron Primus), equipped with radiosurgery applicators and MLC. Narrow circular and Z-shaped on-axis and off-axis fields, as well as broad IMRT configured beams, have been simulated together with reference 10 x 10 cm2 beams. Phase-space data have been used to generate 3D dose distributions which have been compared satisfactorily with experimental profiles (ion chamber, diodes and film). Photon and electron spectra at various depths in water have been calculated, followed by Spencer-Attix (delta = 10 keV) stopping-power ratio calculations which have been compared to those used in the IAEA TRS-398 code of practice. For water/air and PMMA/air stopping-power ratios, agreements within 0.1% have been obtained for the 10 x 10 cm2 fields. For radiosurgery applicators and narrow MLC beams, the calculated s(w,air) values agree with the reference within +/-0.3%, well within the estimated standard uncertainty of the reference stopping-power ratios (0.5%). Ionization chamber dosimetry of narrow beams at the photon qualities used in this work (6 MV) can therefore be based on stopping-power ratios data in dosimetry protocols. For a modulated 6 MV broad beam used in clinical IMRT, s(w,air) agrees within 0.1% with the value for 10 x 10 cm2, confirming that at low energies IMRT absolute dosimetry can also be based on data for open reference fields. At higher energies (24 MV) the difference in s(w,air) was up to 1.1%, indicating that the use of protocol data for narrow beams in such cases is less accurate than at low energies, and detailed calculations of the dosimetry parameters involved should be performed if similar accuracy to that of 6 MV is sought.

Routine IMRT verification by means of an automated Monte Carlo simulation system

PURPOSE A tool to simulate complete intensity-modulated radiation therapy (IMRT) treatments with the Monte Carlo (MC) method has been developed. This application is based on a distribution model to employ as short processing times as possible for an operative verification. MATERIALS AND METHODS The Clinical Primus-Siemens Linac beam was simulated with MC, using the EGS4 OMEGA-BEAM code package. An additional home-made program prepares the appropriate parameters for the code, using as input the file sent from the planning system to the linac. These parameters are adapted to the simulation code, making physical and clinical subdivisions of the global simulation of the treatment. Each resultant partition is ordered to a client personal computer in a cluster with 47 machines under a Linux environment. The verification procedure starts delivering the treatment on a plastic phantom containing an ionization chamber. If differences are less than 2%, films are inserted at selected planes in the phantom and the treatment is delivered again to evaluate the relative doses. When matching between treatment planning system (TPS), film, and MC is acceptable, a new evaluation of the patient is then performed between TPS and MC. Three different cases are shown to prove the applicability of the verification model. RESULTS Acceptable agreement between the three methods used was obtained. The results are presented using different analysis tools. The actual time employed to simulate the total treatment in each case was no more than 5 h, depending on the number of segments. CONCLUSIONS The MC model presented is fully automated, and results can be achieved within the operative time limits. The procedure is a reliable tool to verify any IMRT treatment.

Automatic determination of primary electron beam parameters in Monte Carlo simulation

In order to obtain realistic and reliable Monte Carlo simulations of medical linac photon beams, an accurate determination of the parameters that define the primary electron beam that hits the target is a fundamental step. In this work we propose a new methodology to commission photon beams in Monte Carlo simulations that ensures the reproducibility of a wide range of clinically useful fields. For such purpose accelerated Monte Carlo simulations of 2 x 2, 10 x 10, and 20 x 20 cm2 fields at SSD = 100 cm are carried out for several combinations of the primary electron beam mean energy and radial FWHM. Then, by performing a simultaneous comparison with the correspondent measurements for these same fields, the best combination is selected. This methodology has been employed to determine the characteristics of the primary electron beams that best reproduce a Siemens PRIMUS and a Varian 2100 CD machine in the Monte Carlo simulations. Excellent agreements were obtained between simulations and measurements for a wide range of field sizes. Because precalculated profiles are stored in databases, the whole commissioning process can be fully automated, avoiding manual fine-tunings. These databases can also be used to characterize any accelerators of the same model from different sites.

Estimation of neutron-equivalent dose in organs of patients undergoing radiotherapy by the use of a novel online digital detector

Neutron peripheral contamination in patients undergoing high-energy photon radiotherapy is considered as a risk factor for secondary cancer induction. Organ-specific neutron-equivalent dose estimation is therefore essential for a reasonable assessment of these associated risks. This work aimed to develop a method to estimate neutron-equivalent doses in multiple organs of radiotherapy patients. The method involved the convolution, at 16 reference points in an anthropomorphic phantom, of the normalized Monte Carlo neutron fluence energy spectra with the kerma and energy-dependent radiation weighting factor. This was then scaled with the total neutron fluence measured with passive detectors, at the same reference points, in order to obtain the equivalent doses in organs. The latter were correlated with the readings of a neutron digital detector located inside the treatment room during phantom irradiation. This digital detector, designed and developed by our group, integrates the thermal neutron fluence. The correlation model, applied to the digital detector readings during patient irradiation, enables the online estimation of neutron-equivalent doses in organs. The model takes into account the specific irradiation site, the field parameters (energy, field size, angle incidence, etc) and the installation (linac and bunker geometry). This method, which is suitable for routine clinical use, will help to systematically generate the dosimetric data essential for the improvement of current risk-estimation models.

A CT-aided PC-based physical treatment planning of TBI: a method for dose calculation.

BACKGROUND AND PURPOSE As for conventional radiotherapy, one of the basic requirements in Total Body Irradiation (TBI) is to know accurately the dose delivered to the entire body. Both the dosimetry and the treatment planning need to be improved. Physical, technical and dosimetrical aspects of TBI have been widely discussed in the literature. However, to our knowledge, no planning systems specifically designed for TBI are commercially available. This article describes a CT-aided PC-based planning system (TBI-Plansys) and its dose calculation algorithm, which applies scatter and inhomogeneity corrections, developed for the TBI technique currently in use at our centre (AP/PA irradiation with patient positioned on his side). MATERIAL AND METHOD A description of the material and method followed in the dosimetrical procedure is included as it constitutes the basis of the proposed dose calculation algorithm (more than 2D). A Windows programming environment has been used to develop the software. RESULTS TBI-Plansys uses patient CT data and indicates absolute and relative dose distributions along midline (at reference points), the transversal axis at the specification point and on transverse sections. The system also calculates the appropriate thicknesses of bolus and shielding to modify undesired dose distributions. TBI-Plansys has been checked against two other well-established systems (beam-zone method and our in vivo semiconductor probe-based system). The checks showed good accuracy with dose differences less than 1% and 3% for homogeneous and inhomogeneous tissues, respectively. CONCLUSIONS CT calculations by TBI-Plansys allow us to detect undesired distributions which may go unnoticed by calculations at only some specific points. The system has shown clear advantages for routine clinical use as it generates more detailed and accurate information than manual calculations and diminishes the time requirements.

The determination of beam quality correction factors: Monte Carlo simulations and measurements

Modern dosimetry protocols are based on the use of ionization chambers provided with a calibration factor in terms of absorbed dose to water. The basic formula to determine the absorbed dose at a users beam contains the well-known beam quality correction factor that is required whenever the quality of radiation used at calibration differs from that of the users radiation. The dosimetry protocols describe the whole ionization chamber calibration procedure and include tabulated beam quality correction factors which refer to 60Co gamma radiation used as calibration quality. They have been calculated for a series of ionization chambers and radiation qualities based on formulae, which are also described in the protocols. In the case of high-energy photon beams, the relative standard uncertainty of the beam quality correction factor is estimated to amount to 1%. In the present work, two alternative methods to determine beam quality correction factors are prescribed-Monte Carlo simulation using the EGSnrc system and an experimental method based on a comparison with a reference chamber. Both Monte Carlo calculations and ratio measurements were carried out for nine chambers at several radiation beams. Four chamber types are not included in the current dosimetry protocols. Beam quality corrections for the reference chamber at two beam qualities were also measured using a calorimeter at a PTB Primary Standards Dosimetry Laboratory. Good agreement between the Monte Carlo calculated (1% uncertainty) and measured (0.5% uncertainty) beam quality correction factors was obtained. Based on these results we propose that beam quality correction factors can be generated both by measurements and by the Monte Carlo simulations with an uncertainty at least comparable to that given in current dosimetry protocols.

Monte Carlo correction factors for a Farmer 0.6 cm3 ion chamber dose measurement in the build-up region of the 6 MV clinical beam

Reference dosimetry of photon fields is a well-established subject and currently available protocols (such as the IAEA TRS-398 and AAPM TG-51) provide methods for converting the ionization chamber (IC) reading into dose to water, provided reference conditions of charged particle equilibrium (CPE) are fulfilled. But these protocols cannot deal with the build-up region, where the lack of CPE limits the applicability of the cavity theorems and so the chamber correction factors become depth dependent. By explicitly including the IC geometry in the Monte Carlo simulations, depth-dependent dose correction factors are calculated for a PTW 30001 0.6 cm(3) ion chamber in the build-up region of the 6 MV photon beam. The corrected percentage depth dose (PDD) agrees within 2% with that measured using the NACP 02 plane-parallel ion chamber in the build-up region at depths greater than 0.4 cm, where the Farmer chamber wall reaches the phantom surface.

Microionization chamber for reference dosimetry in IMRT verification: clinical implications on OAR dosimetric errors

Intensity modulated radiotherapy (IMRT) has become a treatment of choice in many oncological institutions. Small fields or beamlets with sizes of 1 to 5 cm2 are now routinely used in IMRT delivery. Therefore small ionization chambers (IC) with sensitive volumes 0.1 cm3 are generally used for dose verification of an IMRT treatment. The measurement conditions during verification may be quite different from reference conditions normally encountered in clinical beam calibration, so dosimetry of these narrow photon beams pertains to the so-called non-reference conditions for beam calibration. This work aims at estimating the error made when measuring the organ at risks (OAR) absolute dose by a micro ion chamber (microIC) in a typical IMRT treatment. The dose error comes from the assumption that the dosimetric parameters determining the absolute dose are the same as for the reference conditions. We have selected two clinical cases, treated by IMRT, for our dose error evaluations. Detailed geometrical simulation of the microIC and the dose verification set-up was performed. The Monte Carlo (MC) simulation allows us to calculate the dose measured by the chamber as a dose averaged over the air cavity within the ion-chamber active volume (D(air)). The absorbed dose to water (D(water)) is derived as the dose deposited inside the same volume, in the same geometrical position, filled and surrounded by water in the absence of the ion chamber. Therefore, the D(water)/D(air) dose ratio is the MC estimator of the total correction factor needed to convert the absorbed dose in air into the absorbed dose in water. The dose ratio was calculated for the microIC located at the isocentre within the OARs for both clinical cases. The clinical impact of the calculated dose error was found to be negligible for the studied IMRT treatments.

Neutron contamination in radiotherapy: Estimation of second cancers based on measurements in 1377 patients

PURPOSE Second cancer, as a consequence of a curative intent radiotherapy (RT), represents a growing concern nowadays. The unwanted neutron exposure is an important contributor to this risk in patients irradiated with high energy photon beams. The design and development by our group of a neutron digital detector, together with the methodology to estimate, from the detector readings, the neutron equivalent dose in organs, made possible the unprecedented clinical implementation of an online and systematic neutron dosimetry system. The aim of this study was to systematically estimate neutron equivalent dose in organs of a large patient group treated in different installations. PATIENTS AND METHODS Neutron dosimetry was carried out in 1377 adult patients at more than 30 different institutions using the new neutron digital detector located inside the RT room. Second cancer risk estimates were performed applying ICRP risk coefficients. RESULTS Averaged equivalent dose in organs ranges between 0.5 mSv and 129 mSv depending on the type of treatment (dose and beam-on time), the distance to isocenter and the linac model. The mean value of the second cancer risk for our patient group is 1.2%. Reference values are proposed for an overall estimation of the risks in 15 linac models (from 2.8 × 10(-5) to 62.7 × 10(-5)%/MU). CONCLUSIONS The therapeutic benefit of RT must outweigh the second cancer risk. Thus, these results should be taken into account when taking clinical decisions regarding treatment strategy choice during RT planning.

Monte Carlo simulation of complex radiotherapy treatments

Monte Carlo simulation is an accurate way of assessing radiotherapy dose distribution in nonhomogeneous volumes, but it requires long processing times. A new distribution model simulates radiotherapy treatments and runs on a PC network, which reduces the processing time and makes for a powerful treatment-verification tool.