Peripheral Organ Equivalent Dose Estimation Procedure in Proton Therapy.

The aim of this work is to present a reproducible methodology for the evaluation of total equivalent doses in organs during proton therapy facilities. The methodology is based on measuring the dose equivalent in representative locations inside an anthropomorphic phantom where photon and neutron dosimeters were inserted. The Monte Carlo simulation was needed for obtaining neutron energy distribution inside the phantom. The methodology was implemented for a head irradiation case in the passive proton beam of iThemba Labs (South Africa). Thermoluminescent dosimeter (TLD)-600 and TLD-700 pairs were used as dosimeters inside the phantom and GEANT code for simulations. In addition, Bonner sphere spectrometry was performed inside the treatment room to obtain the neutron spectra, some relevant neutron dosimetric quantities per treatment Gy, and a percentual distribution of neutron fluence and ambient dose equivalent in four energy groups, at two locations. The neutron spectrum at one of those locations was also simulated so that a reasonable agreement between simulation and measurement allowed a validation of the simulation. Results showed that the total out-of-field dose equivalent inside the phantom ranged from 1.4 to 0.28 mSv/Gy, mainly due to the neutron contribution and with a small contribution from photons, 10% on average. The order of magnitude of the equivalent dose in organs was similar, displaying a slow reduction in values as the organ is farther from the target volume. These values were in agreement with those found by other authors in other passive beam facilities under similar irradiation and measurement conditions.

Intensity-modulated radiation therapy and volumetric modulated arc therapy versus conventional conformal techniques at high energy: Dose assessment and impact on second primary cancer in the out-of-field region.

The aim of this work was to estimate peripheral neutron and photon doses associated with the conventional 3D conformal radiotherapy techniques in comparison to modern ones such as Intensity modulated radiation therapy and volumetric modulated arc therapy. Assessment in terms of second cancer incidence ought to peripheral doses was also considered. For that, a dosimetric methodology proposed by the authors has been applied beyond the region where there is no CT information and, thus, treatment planning systems do not calculate and where, nonetheless, about one third of second primary cancers occurs.

Uncomplicated and Cancer-Free Control Probability (UCFCP): A new integral approach to treatment plan optimization in photon radiation therapy.

PURPOSE: Biological treatment plan evaluation does not currently consider second cancer induction from peripheral doses associated to photon radiotherapy. The aim is to propose a methodology to characterize the therapeutic window by means of an integral radiobiological approach, which considers not only Tumour Control Probability (TCP) and Normal Tissue Complication Probability (NTCP) but also Secondary Cancer Probability (SCP). METHODS: Uncomplicated and Cancer-Free Control Probability (UCFCP) function has been proposed assuming a statistically uncorrelated response for tumour and normal tissues. The Poisson's and Lyman's models were chosen for TCP and NTCP calculations, respectively. SCP was modelled as the summation of risks associated to photon and neutron irradiation of radiosensitive organs. For the medium (>4Gy) and low dose regions, mechanistic and linear secondary cancer risks models were used, respectively. Two conformal and intensity-modulated prostate plans at 15MV (same prescription dose) were selected to illustrate the UCFCP features. RESULTS: UCFCP exhibits a bell-shaped behaviour with its maximum inside the therapeutic window. SCP values were not different for the plans analysed (∼2.4%) and agreed with published epidemiological results. Therefore, main differences in UCFCP came from differences in rectal NTCP (18% vs 9% for 3D-CRT and IMRT, respectively). According to UCFCP values, the evaluated IMRT plan ranked first. CONCLUSIONS: The level of SCP was found to be similar to that of NTCP complications which reinforces the importance of considering second cancer risks as part of the possible late sequelae due to treatment. Previous concerns about the effect of peripheral radiation, especially neutrons, in the induction of secondary cancers can be evaluated by quantifying the UCFCP.

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.

Uncertainty estimation in intensity-modulated radiotherapy absolute dosimetry verification.

PURPOSE: Intensity-modulated radiotherapy (IMRT) represents an important method for improving RT. The IMRT relative dosimetry checks are well established; however, open questions remain in reference dosimetry with ionization chambers (ICs). The main problem is the departure of the measurement conditions from the reference ones; thus, additional uncertainty is introduced into the dose determination. The goal of this study was to assess this effect systematically. METHODS AND MATERIALS: Monte Carlo calculations and dosimetric measurements with five different detectors were performed for a number of representative IMRT cases, covering both step-and-shoot and dynamic delivery. RESULTS: Using ICs with volumes of about 0.125 cm(3) or less, good agreement was observed among the detectors in most of the situations studied. These results also agreed well with the Monte Carlo-calculated nonreference correction factors (c factors). Additionally, we found a general correlation between the IC position relative to a segment and the derived correction factor c, which can be used to estimate the expected overall uncertainty of the treatment. CONCLUSION: The increase of the reference dose relative standard uncertainty measured with ICs introduced by nonreference conditions when verifying an entire IMRT plan is about 1-1.5%, provided that appropriate small-volume chambers are used. The overall standard uncertainty of the measured IMRT dose amounts to about 2.3%, including the 0.5% of reproducibility and 1.5% of uncertainty associated with the beam calibration factor. Solid state detectors and large-volume chambers are not well suited to IMRT verification dosimetry because of the greater uncertainties. An action level of 5% is appropriate for IMRT verification. Greater discrepancies should lead to a review of the dosimetric procedure, including visual inspection of treatment segments and energy fluence.

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.

Micro ionization chamber dosimetry in IMRT verification: clinical implications of dosimetric errors in the PTV.

BACKGROUND AND PURPOSE: Absolute dose measurements for Intensity Modulated Radiotherapy (IMRT) beamlets is difficult due to the lack of lateral electron equilibrium. Recently we found that the absolute dosimetry in the penumbra region of the IMRT beamlet, can suffer from significant errors (Capote et al., Med Phys 31 (2004) 2416-2422). This work has the goal to estimate the error made when measuring the Planning Target Volume's (PTV) absolute dose by a micro ion chamber (microIC) in 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. MATERIALS AND METHODS: Two IMRT treatment plans for common prostate carcinoma case, derived by forward and inverse optimisation, were considered. Detailed geometrical simulation of the microIC and the dose verification set-up was performed. The Monte Carlo (MC) simulation allows us to calculate the delivered dose to water and the dose delivered to the active volume of the ion chamber. However, the measured dose in water is usually derived from chamber readings assuming reference conditions. The MC simulation provides needed correction factors for ion chamber dosimetry in non reference conditions. RESULTS: Dose calculations were carried out for some representative beamlets, a combination of segments and for the delivered IMRT treatments. We observe that the largest dose errors (i.e. the largest correction factors) correspond to the smaller contribution of the corresponding IMRT beamlets to the total dose delivered in the ionization chamber within PTV. CONCLUSION: The clinical impact of the calculated dose error in PTV measured dose was found to be negligible for studied IMRT treatments.

Additional dose constraints for analytical beam weighting optimization in IMRT.

This work presents an improvement to an algorithm for analytical beam weighting optimization where a flexible objective function, which considers 'importance factors' for each anatomical region and 'allowed deviations' from the prescribed dose, is defined. This upgrading allows forcing the mean value of the dose distribution to be the desired value, by using Lagrange multipliers. A real case is presented to show the effect of this change.

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 risk's (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.

An EGSnrc Monte Carlo study of the microionization chamber for reference dosimetry of narrow irregular IMRT beamlets.

Intensity modulated radiation therapy (IMRT) has evolved toward the use of many small radiation fields, or "beamlets," to increase the resolution of the intensity map. The size of smaller beamlets can be typically about 1-5 cm2. Therefore small ionization chambers (IC) with sensitive volumes < or = 0.1 cm3 are generally used for dose verification of IMRT treatment. The dosimetry of these narrow photon beams pertains to the so-called nonreference conditions for beam calibration. The use of ion chambers for such narrow beams remains questionable due to the lack of electron equilibrium in most of the field. The present contribution aims to estimate, by the Monte Carlo (MC) method, the total correction needed to convert the IBA-Wellhöfer NAC007 micro IC measured charge in such radiation field to the absolute dose to water. Detailed geometrical simulation of the microionization chamber was performed. The ion chamber was always positioned at a 10 cm depth in water, parallel to the beam axis. The delivered doses to air and water cavity were calculated using the CAVRZ EGSnrc user code. The 6 MV phase-spaces for Primus Clinac (Siemens) used as an input to the CAVRZnrc code were derived by BEAM/EGS4 modeling of the treatment head of the machine along with the multileaf collimator [Sánchez-Doblado et al., Phys. Med. Biol. 48, 2081-2099 (2003)] and contrasted with experimental measurements. Dose calculations were carried out for two irradiation geometries, namely, the reference 10x10 cm2 field and an irregular (approximately 2x2 cm2) IMRT beamlet. The dose measured by the ion chamber is estimated by MC simulation as a dose averaged over the air cavity inside the ion-chamber (Dair). The absorbed dose to water is derived as the dose deposited inside the same volume, in the same geometrical position, filled and surrounded by water (Dwater) in the absence of the ionization chamber. Therefore, the Dwater/Dair dose ratio is a MC direct estimation of the total correction factor needed to convert the absorbed dose in air to absorbed dose to water. The dose ratio was calculated for several chamber positions, starting from the penumbra region around the beamlet along the two diagonals crossing the radiation field. For this quantity from 0 up to a 3% difference is observed between the dose ratio values obtained within the small irregular IMRT beamlet in comparison with the dose ratio derived for the reference 10x10 cm2 field. Greater differences from the reference value up to 9% were obtained in the penumbra region of the small IMRT beamlet.

MLC leaf width impact on the clinical dose distribution: a Monte Carlo approach.

PURPOSE: The influence of the multileaf collimator (MLC) leaf width on the dose distribution in patients treated with conformal radiotherapy and intensity-modulated radiotherapy has been analyzed. This study was based on the Monte Carlo simulation with the beams generated by a linac with the double-focused MLC. MATERIALS AND METHODS: The transmission through the leaves and the exact shape of the penumbra regions are difficult to model by treatment planning system algorithms. An accurate assessment of the dose variations due to the leaf width change can be achieved by means of Monte Carlo simulation. The BEAM/EGS4 code was used at the Hospital of the Virgen Macarena to model a Siemens PRIMUS linac, featuring an MLC with a leaf width projecting 1 cm at the isocenter. Based on this real model, a virtual head was designed while allowing for a variation of the leaf width projection. Both the real linac and the virtual linac, with leaves projecting 0.5 cm, were used to obtain the dose distributions for several treatments. A few disease sites, including the prostate, head and neck, and endometrium, were selected for the design of the conformal and intensity-modulated radiotherapy treatments with a forward planning algorithm sensitive to the different shapes of the volumes of interest. Isodose curves, differential matrix, gamma function, and the dose-volume histograms (DVHs) corresponding to both MLC models were obtained for all cases. The tumor control probability and the normal tissue complication probability were derived for those cases studied featuring the greatest differences between results for both MLCs. RESULTS: The impact on the DVHs of changing leaf width projections at the isocenter from 1.0 cm to 0.5 cm was low. Radiobiologic models showed slightly better tumor control probability/normal tissue complication probability values using the virtual MLC with a leaf width projecting 0.5 cm at isocenter in those cases presenting greater differences in the DVHs. CONCLUSIONS: The impact on the clinical dose distribution due to the MLC leaf width change is low based on the design and conditions used in this study.

Computer optimization of class solutions designed on a beam segmentation basis.

BACKGROUND AND PURPOSE: A method for analytically solving the optimization of beam weighting in radiotherapy treatments using beam segmentation is presented. PATIENTS AND METHODS: A technique has been elaborated permitting the optimization of a flexible objective function which is defined by considering 'importance factors' for each anatomical region and 'allowed deviations' from the prescribed dose. As any change in these importance factors may lead to very different solutions, a statistical tool has been developed which varies the objective function automatically and iteratively to get the best possible results compatible with the chosen class solution. In addition, the spatial symmetry found in many anatomical sites is taken advantage of. Furthermore, freeware code has been written to run this optimization approach. The guidelines to design the beam segmentation used in our institution using organ avoidance criteria, and hence the suitable class solution for different anatomical sites, are given. A treatment-planning study for three anatomical sites is presented and, for two of them, the results obtained with both the suggested and the classical inverse approach are presented. RESULTS: The work presented might be used for beam weighting optimization in any radiotherapy treatment and furthermore, the suggested procedure may successfully confront the intensity-modulated radiation therapy (IMRT) problem and the obtained dose distributions fit the clinical constraints for all anatomical sites studied. When comparing with the classical inverse approach, both results are comparable in terms of dose distribution, but the suggested technique reduces the integral dose as the total number of monitor units is lower. CONCLUSIONS: The developed code performs the optimization with a very low time cost and in addition, this process can be carried out with a conventional treatment-planning system with no need of dedicated IMRT software.

Total skin electron therapy treatment verification: Monte Carlo simulation and beam characteristics of large non-standard electron fields.

Total skin electron therapy (TSET) is a complex technique which requires non-standard measurements and dosimetric procedures. This paper investigates an essential first step towards TSET Monte Carlo (MC) verification. The non-standard 6 MeV 40 x 40 cm2 electron beam at a source to surface distance (SSD) of 100 cm as well as its horizontal projection behind a polymethylmethacrylate (PMMA) screen to SSD = 380 cm were evaluated. The EGS4 OMEGA-BEAM code package running on a Linux home made 47 PCs cluster was used for the MC simulations. Percentage depth-dose curves and profiles were calculated and measured experimentally for the 40 x 40 cm2 field at both SSD = 100 cm and patient surface SSD = 380 cm. The output factor (OF) between the reference 40 x 40 cm2 open field and its horizontal projection as TSET beam at SSD = 380 cm was also measured for comparison with MC results. The accuracy of the simulated beam was validated by the good agreement to within 2% between measured relative dose distributions, including the beam characteristic parameters (R50, R80, R100, Rp, E0) and the MC calculated results. The energy spectrum, fluence and angular distribution at different stages of the beam (at SSD = 100 cm, at SSD = 364.2 cm, behind the PMMA beam spoiler screen and at treatment surface SSD = 380 cm) were derived from MC simulations. Results showed a final decrease in mean energy of almost 56% from the exit window to the treatment surface. A broader angular distribution (FWHM of the angular distribution increased from 13 degrees at SSD = 100 cm to more than 30 degrees at the treatment surface) was fully attributable to the PMMA beam spoiler screen. OF calculations and measurements agreed to less than 1%. The effect of changing the electron energy cut-off from 0.7 MeV to 0.521 MeV and air density fluctuations in the bunker which could affect the MC results were shown to have a negligible impact on the beam fluence distributions. Results proved the applicability of using MC as a treatment verification tool for complex radiotherapy techniques.

Investigation of radiosurgical beam profiles using Monte Carlo method.

An accurate determination of the penumbra of radiosurgery profiles is critical to avoid complications in organs at risk adjacent to the tumor. Conventional detectors may not be accurate enough for small field sizes. The Monte Carlo (MC) method was used to study the behavior of radiosurgical beam profiles at the penumbral region; the BEAM code was also used in this work. Two collimators (2.2- and 0.3-cm diameter) were calculated and compared with empirical measurements obtained with the detectors normally used. The differences found between film dosimetry and MC revealed a systematic error in the reading procedure. In the process, a water phantom was simulated with a layer of the same composition as that of the film. MC calculations with film differed by a small amount from those obtained with the water phantom alone. In conclusion, MC may be used as a verification tool to support dosimetrical procedures with conventional detectors, especially in very small beams such as those used in radiosurgery. Furthermore, it has been proved that the film energy dependence is negligible for fields used in radiosurgery.

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.