Underdosage of the upper-airway mucosa for small fields as used in intensity-modulated radiation therapy: a comparison between radiochromic film measurements, Monte Carlo simulations, and collapsed cone convolution calculations

Synergistic Radiosensitization Mediated by Chemodynamic Therapy via a Novel Biodegradable Peroxidases Mimicking Nanohybrid

Purpose

Reactive oxygen species (ROS) are practically essential in radiotherapy to damage cancer cells; however, they are always inadequate for some malignant entities. Here, we designed a biodegradable mesoporous silica decorated with hemin and glucose oxidase (GOD@Hemin-MSN) to generate a chemodynamic therapy in order to enhance the killing capacity of radiotherapy.

Methods

Mesoporous silica, as an outstanding drug carrier, can deliver hemin and glucose oxidase to the tumor site. With high level of metabolism activity, cancer cells are abundant in glucose, which can be oxidized into H₂O₂ by glucose oxidase (GOD) on site. The generated H₂O₂ is subsequently converted into intracellular ROS, especially hydroxyl radical within the tumor microenvironment, by hemin, which has mimetic peroxidase properties. By this means, the ROS can be supplemented or enriched to facilitate the killing of tumor cells.

Results

The chemodynamic therapy induced by GOD@Hemin-MSN produced quantities of ROS, which compensated for their inadequacy as a result of radiotherapy, and exhibited remarkable antitumor efficacy, with a tumor inhibition rate of 91.5% in A549 tumor-bearing mice.

Conclusion

This work has validated GOD@Hemin-MSN as a radiosensitizer in chemodynamic therapy, which showed biocompatibility and potential for translational application.

Accuracy Evaluation of Collapsed Cone Convolution Superposition Algorithms for the Nasopharynx Interface in the Early Stage of Nasopharyngeal Carcinoma

This study combined the use of radiation dosimeteric measurements and a custom-made anthropomorphic phantom in order to evaluate the accuracy of therapeutic dose calculations at the nasopharyngeal air-tissue interface. The doses at the nasopharyngeal air-tissue interface obtained utilizing the Pinnacle and TomoTherapy TPS, which are based on collapsed cone convolution superposition (CCCS) algorithms, were evaluated and measured under single $10\times 10\text{c}{\text{m}}^2$ , $2\times 2\text{c}{\text{m}}^2$ , two parallel opposed $2\times 2\text{c}{\text{m}}^2$ and clinical fields for early stage of nasopharyngeal carcinoma by using EBT3, GR-200F, and TLD 100. At the air-tissue interface under a $10\times 10\text{c}{\text{m}}^2$ field, the TPS dose calculation values were in good agreement with the dosimeter measurement with all differences within 3.5%. When measured the single field $2\times 2\text{c}{\text{m}}^2$ , the differences between the average dose were measured at the distal interface for EBT3, GR-200F, and TLD-100 and the calculation values were -15.8%, -16.4%, and -4.9%, respectively. When using the clinical techniques such as IMRT, VMAT, and tomotherapy, the measurement results at the interface for all three techniques did not imply under dose. Small-field sizes will lead to dose overestimation at the nasopharyngeal air-tissue interface due to electronic disequilibrium when using CCCS algorithms. However, under clinical applications of multiangle irradiation, the dose errors caused by this effect were not significant.

Maintaining dosimetric quality when switching to a Monte Carlo dose engine for head and neck volumetric-modulated arc therapy planning

Head and neck cancers present challenges in radiation treatment planning due to the large number of critical structures near the target(s) and highly heterogeneous tissue composition. While Monte Carlo (MC) dose calculations currently offer the most accurate approximation of dose deposition in tissue, the switch to MC presents challenges in preserving the parameters of care. The differences in dose‐to‐tissue were widely discussed in the literature, but mostly in the context of recalculating the existing plans rather than reoptimizing with the MC dose engine. Also, the target dose homogeneity received less attention. We adhere to strict dose homogeneity objectives in clinical practice. In this study, we started with 21 clinical volumetric‐modulated arc therapy (VMAT) plans previously developed in Pinnacle treatment planning system. Those plans were recalculated “as is” with RayStation (RS) MC algorithm and then reoptimized in RS with both collapsed cone (CC) and MC algorithms. MC statistical uncertainty (0.3%) was selected carefully to balance the dose computation time (1–2 min) with the planning target volume (PTV) dose‐volume histogram (DVH) shape approaching that of a “noise‐free” calculation. When the hot spot in head and neck MC‐based treatment planning is defined as dose to 0.03 cc, it is exceedingly difficult to limit it to 105% of the prescription dose, as we were used to with the CC algorithm. The average hot spot after optimization and calculation with RS MC was statistically significantly higher compared to Pinnacle and RS CC algorithms by 1.2 and 1.0 %, respectively. The 95% confidence interval (CI) observed in this study suggests that in most cases a hot spot of ≤107% is achievable. Compared to the 95% CI for the previous clinical plans recalculated with RS MC “as is” (upper limit 108%), in real terms this result is at least as good or better than the historic plans.

Comparative study of the collapsed cone convolution and Monte Carlo algorithms for radiation therapy planning of canine sinonasal tumors reveals significant dosimetric differences

Computer-based radiation therapy requires high targeting and dosimetric precision. Analytical dosimetric algorithms typically are fast and clinically viable but can have increasing errors near air-bone interfaces. These are commonly found within dogs undergoing radiation planning for sinonasal cancer. This retrospective methods comparison study is designed to compare the dosimetry of both tumor volumes and organs at risk and quantify the differences between collapsed cone convolution (CCC) and Monte Carlo (MC) algorithms. Canine sinonasal tumor plans were optimized with CCC and then recalculated by MC with identical control points and monitor units. Planning target volume (PTV)_air , PTV_{soft tissue} , and PTV_bone were created to analyze the dose discrepancy within the PTV. Thirty imaging sets of dogs were included. Monte Carlo served as the gold standard calculation for the dosimetric comparison. Collapsed cone convolution overestimated the mean dose (D_mean ) to PTV and PTV_{soft tissue} by 0.9% and 0.5%, respectively (both P < 0.001). Collapsed cone convolution overestimated D_mean to PTV_bone by 3% (P < 0.001). Collapsed cone convolution underestimated the near-maximum dose (D₂ ) to PTV_air by 1.1% (P < 0.001), and underestimated conformity index and homogeneity index in PTV (both P < 0.001). Mean doses of contralateral and ipsilateral eyes were overestimated by CCC by 1.6% and 1.7%, respectively (both P < 0.001). Near-maximum doses of skin and brain were overestimated by CCC by 2.2% and 0.7%, respectively (both P < 0.001). As clinical accessibility of Monte Carlo becomes more widespread, dose constraints may need to be re-evaluated with appropriate plan evaluation and follow-up.

Report of AAPM Task Group 155: Megavoltage photon beam dosimetry in small fields and non-equilibrium conditions

Small-field dosimetry used in advance treatment technologies poses challenges due to loss of lateral charged particle equilibrium (LCPE), occlusion of the primary photon source, and the limited choice of suitable radiation detectors. These challenges greatly influence dosimetric accuracy. Many high-profile radiation incidents have demonstrated a poor understanding of appropriate methodology for small-field dosimetry. These incidents are a cause for concern because the use of small fields in various specialized radiation treatment techniques continues to grow rapidly. Reference and relative dosimetry in small and composite fields are the subject of the International Atomic Energy Agency (IAEA) dosimetry code of practice that has been published as TRS-483 and an AAPM summary publication (IAEA TRS 483; Dosimetry of small static fields used in external beam radiotherapy: An IAEA/AAPM International Code of Practice for reference and relative dose determination, Technical Report Series No. 483; Palmans et al., Med Phys 45(11):e1123, 2018). The charge of AAPM task group 155 (TG-155) is to summarize current knowledge on small-field dosimetry and to provide recommendations of best practices for relative dose determination in small megavoltage photon beams. An overview of the issue of LCPE and the changes in photon beam perturbations with decreasing field size is provided. Recommendations are included on appropriate detector systems and measurement methodologies. Existing published data on dosimetric parameters in small photon fields (e.g., percentage depth dose, tissue phantom ratio/tissue maximum ratio, off-axis ratios, and field output factors) together with the necessary perturbation corrections for various detectors are reviewed. A discussion on errors and an uncertainty analysis in measurements is provided. The design of beam models in treatment planning systems to simulate small fields necessitates special attention on the influence of the primary beam source and collimating devices in the computation of energy fluence and dose. The general requirements for fluence and dose calculation engines suitable for modeling dose in small fields are reviewed. Implementations in commercial treatment planning systems vary widely, and the aims of this report are to provide insight for the medical physicist and guidance to developers of beams models for radiotherapy treatment planning systems.

Monte Carlo as quality control tool of stereotactic body radiation therapy treatment plans

PURPOSE/OBJECTIVE

The objective of this study was to verify the accuracy of treatment plans of stereotactic body radiation therapy (SBRT) and to verify the feasibility of the use of Monte Carlo (MC) as quality control (QC) on a daily basis.

MATERIAL/METHODS

Using EGSnrc, a MC model of Agility™ linear accelerator was created. Various measurements (Percentage depth dose (PDD), Profiles and Output factors) were done for different fields sizes from 1x1 up to 40x40 (cm²). An iterative model optimization was performed to achieve adequate parameters of MC simulation. 40 SBRT patient's dosimetry plans were calculated by Monaco™ 3.1.1. CT images, RT-STRUCT and RT-PLAN files from Monaco™ being used as input for Moderato MC code. Finally, dose volume histogram (DVH) and paired t-tests for each contour were used for dosimetry comparison of the Monaco™ and MC.

RESULTS

Validation of MC model was successful, as <2% difference comparing to measurements for all field's sizes. The main energy of electron source incident on the target was 5.8 MeV, and the full width at half maximum (FWHM) of Gaussian electron source were 0.09 and 0.2 (cm) in X and Y directions, respectively. For 40 treatment plan comparisons, the minimum absolute difference of mean dose of planning treatment planning (PTV) was 0.1% while the maximum was 6.3%. The minimum absolute difference of Max dose of PTV was 0.2% while the maximum was 8.1%.

CONCLUSION

SBRT treatment plans of Monaco agreed with MC results. It possible to use MC for treatment plans verifications as independent QC tool.

Assessing localized dosimetric effects due to unplanned gas cavities during pelvic MR-guided radiotherapy using Monte Carlo simulations

PURPOSE

It has been proposed that beam modulation and opposing beam configurations can cancel effects of the Electron Return Effect (ERE) during MR-guided radiotherapy (MRgRT). However, this may not always be the case for unplanned gas cavities outside of the target in the pelvic region. We evaluate dosimetric effects, including effects in the rectal wall, due to unplanned spherical air cavities during MRgRT.

METHODS

Nine virtual cuboid water phantoms containing spherical air cavities (0.5-7.5 cm diameter) and a reference phantom without air were created. Monte Carlo dose calculations of 7 MV photons under the influence of a 1.5 T transverse magnetic field were produced using Monaco 5.19.02 Treatment Planning System (TPS) (Elekta AB, Stockholm, Sweden). Cavities in the path of a single and multiple beam plans were considered. Dose distributions of phantoms with and without air cavities were compared (ΔD% ) using a spherical coordinate system originating in the center of the cavity. Effects in the rectal wall were quantified by comparing dose volume histogram (DVH) parameters for solid and gaseous filling from simulated rectal wall structures.

RESULTS

Max(ΔD% ) of ~70% and 20% were observed around large cavities in the path of a single and multiple beam plans, respectively. Approximately 45 cm3 of phantom surrounding the largest cavity in a single beam received dose changes of >10%. Dmean in the rectal wall was unchanged when comparing gaseous and solid filling in the path of a single beam; however, D1cc and Dmax increased by up to ~45% and ~63%, respectively.

CONCLUSIONS

Unplanned gas cavities in the path of a single beam during pelvic MRgRT with a 1.5 T transverse magnetic field cause dose changes which may impact toxicity in the rectal wall, depending on local dose and fractionation. Effects are reduced but not eliminated with a five-beam plan.

Skin dose calculation during radiotherapy of head and neck cancer using deformable image registration of planning and mega-voltage computed tomography scans

Background and Purpose

Head-Neck (HN) patients may experience severe acute skin complications that can cause treatment interruption and increase the risk of late fibrosis. This study assessed a method for accurately monitoring skin dose changes during helical tomotherapy for HN cancer based on deformable image registration of planning computed tomography (CT) and mega-voltage CT (MVCT).

Materials and Methods

Planning CTs of nine patients were deformably registered to mid-treatment MVCT (MV15) images resulting in CTdef images. The original plans were recalculated on both CTdef and mid-treatment kilo-voltage CT (CT15) taken as ground truth. Superficial layers (SL) of the body with thicknesses of 2, 3 and 5 mm (SL2, SL3, SL5) were considered as derma surrogates. SL V95%, V97%, V98%, V100%, V102%, V105% and V107% of the prescribed PTV dose were extracted for CT15/CTdef and compared (considering patients with skin dose > 95%). For comparison, doses were calculated directly on the calibrated MVCT and analyzed in the same way.

Results

Differences between SL2/SL3/SL5 V95%-V107% in CT15/CTdef were very small: for eight of nine patients the difference between the considered SL2 Vd% computed on CTdef and CT15 was less than 1.4 cm³ for all d%. A larger value was found when using MVCT for skin dose calculation (4.8 cm³ for SL2), although CTdef body contour matched CT15 body with accuracy similar to that of MV15.

Conclusions

Deforming the planning CT-to-MVCT was shown to be accurate considering external body contours and skin DVHs. The method was able to accurately identify superficial overdosing.

Experimental Validation of Monte Carlo Simulations Based on a Virtual Source Model for TomoTherapy in a RANDO Phantom

A virtual source model for Monte Carlo simulations of helical TomoTherapy has been developed previously by the authors. The purpose of this work is to perform experiments in an anthropomorphic (RANDO) phantom with the same order of complexity as in clinical treatments to validate the virtual source model to be used for quality assurance secondary check on TomoTherapy patient planning dose. Helical TomoTherapy involves complex delivery pattern with irregular beam apertures and couch movement during irradiation. Monte Carlo simulation, as the most accurate dose algorithm, is desirable in radiation dosimetry. Current Monte Carlo simulations for helical TomoTherapy adopt the full Monte Carlo model, which includes detailed modeling of individual machine component, and thus, large phase space files are required at different scoring planes. As an alternative approach, we developed a virtual source model without using the large phase space files for the patient dose calculations previously. In this work, we apply the simulation system to recompute the patient doses, which were generated by the treatment planning system in an anthropomorphic phantom to mimic the real patient treatments. We performed thermoluminescence dosimeter point dose and film measurements to compare with Monte Carlo results. Thermoluminescence dosimeter measurements show that the relative difference in both Monte Carlo and treatment planning system is within 3%, with the largest difference less than 5% for both the test plans. The film measurements demonstrated 85.7% and 98.4% passing rate using the 3 mm/3% acceptance criterion for the head and neck and lung cases, respectively. Over 95% passing rate is achieved if 4 mm/4% criterion is applied. For the dose–volume histograms, very good agreement is obtained between the Monte Carlo and treatment planning system method for both cases. The experimental results demonstrate that the virtual source model Monte Carlo system can be a viable option for the accurate dose calculation of helical TomoTherapy.

Clinical and dosimetric implications of intensity-modulated radiotherapy for early-stage glottic carcinoma

Conventional parallel-opposed radiotherapy (PORT) is the established standard technique for early-stage glottic carcinoma. However, case reports have reported the utility of intensity-modulated radiotherapy (IMRT) and volumetric-modulated arc therapy (VMAT) with or without image guidance (image-guided radiotherapy, IGRT) in select patients. The proposed advantages of IMRT/VMAT include sparing of the carotid artery, thyroid gland, and the remaining functional larynx, although these benefits remain unclear. The following case study presents a patient with multiple vascular comorbidities treated with VMAT for early-stage glottic carcinoma. A detailed explanation of the corresponding treatment details, dose-volume histogram (DVH) analysis, and a review of the relevant literature are provided. Conventional PORT remains the standard of care for early-stage glottic carcinoma. IMRT or VMAT may be beneficial for select patients, although great care is necessary to avoid a geographical miss. Clinical data supporting the benefit of CRT are lacking. Therefore, these techniques should be used with caution and only in selected patients.

A virtual source model for Monte Carlo simulation of helical tomotherapy

The purpose of this study was to present a Monte Carlo (MC) simulation method based on a virtual source, jaw, and MLC model to calculate dose in patient for helical tomotherapy without the need of calculating phase‐space files (PSFs). Current studies on the tomotherapy MC simulation adopt a full MC model, which includes extensive modeling of radiation source, primary and secondary jaws, and multileaf collimator (MLC). In the full MC model, PSFs need to be created at different scoring planes to facilitate the patient dose calculations. In the present work, the virtual source model (VSM) we established was based on the gold standard beam data of a tomotherapy unit, which can be exported from the treatment planning station (TPS). The TPS‐generated sinograms were extracted from the archived patient XML (eXtensible Markup Language) files. The fluence map for the MC sampling was created by incorporating the percentage leaf open time (LOT) with leaf filter, jaw penumbra, and leaf latency contained from sinogram files. The VSM was validated for various geometry setups and clinical situations involving heterogeneous media and delivery quality assurance (DQA) cases. An agreement of <1% was obtained between the measured and simulated results for percent depth doses (PDDs) and open beam profiles for all three jaw settings in the VSM commissioning. The accuracy of the VSM leaf filter model was verified in comparing the measured and simulated results for a Picket Fence pattern. An agreement of <2% was achieved between the presented VSM and a published full MC model for heterogeneous phantoms. For complex clinical head and neck (HN) cases, the VSM‐based MC simulation of DQA plans agreed with the film measurement with 98% of planar dose pixels passing on the 2%/2 mm gamma criteria. For patient treatment plans, results showed comparable dose‐volume histograms (DVHs) for planning target volumes (PTVs) and organs at risk (OARs). Deviations observed in this study were consistent with literature. The VSM‐based MC simulation approach can be feasibly built from the gold standard beam model of a tomotherapy unit. The accuracy of the VSM was validated against measurements in homogeneous media, as well as published full MC model in heterogeneous media.

PACS numbers: 87.53.‐j, 87.55.K‐

Evaluation of a novel secondary check tool for intensity-modulated radiotherapy treatment planning

The purpose of this study was to assess the accuracy and efficacy of an automated treatment plan verification, or “secondary check”, tool (Mobius3D), which uses a reference dataset to perform an independent three‐dimensional dose verification of the treatment planning system (TPS) dose calculation and assesses plan quality by comparing dose‐volume histograms to reference benchmarks. The accuracy of the Mobius3D (M3D) system was evaluated by comparing dose calculations from IMRT and VMAT plans with measurements in phantom geometries and with TPS calculated dose distributions in prostate, lung, and head and neck patients (ten each). For the patient cases, instances of DVH limits exceeding reference values were also recorded. M3D showed agreement with measured point and planar doses that was comparable to the TPS in phantom geometries. No statistically significant differences (p<0.05) were noted. M3D dose distributions from VMAT plans in patient cases were in good agreement with the TPS, with an average of 99.5% of dose points showing γ5%,3mm<1. The M3D system also identified several plans that had exceeded dose‐volume limits specified by RTOG protocols for those sites. The M3D system showed dosimetric accuracy comparable with the TPS, and identified several plans that exceeded dosimetric benchmarks. The M3D system possesses the potential to enhance the current treatment plan verification paradigm and improve safety in the clinical treatment planning and review process.

PACS number: 87.55.D‐, 87.55.km, 87.55.Qr,

Evaluation of the accuracy of dose delivery for IMRT based on transit dosimetry

The authors have developed and evaluated a new method that uses transit dose in the treatment room in verification of dose delivery to patients. Five intensity modulated radiotherapy (IMRT) plans were selected from actual cancer patients, and transit dose was measured using MapCHECK 2 and an anthropomorphic phantom. The criteria used as a verification tool for the dose delivery to the patient were gamma-index-based dose comparison between the computed dose and measured dose. When the doses were delivered to an anthropomorphic phantom normally, the average passing rate was 95.2% based on a gamma index analysis. This feasibility study suggested that transit dose-based quality assurance can provide information about the accuracy of an inhomogeneity correction algorithm and patient positioning during treatment, allowing its use as a verification tool for actual dose delivery to patients in the treatment room.

A review on the use of grid-based Boltzmann equation solvers for dose calculation in external photon beam treatment planning

Deterministic linear Boltzmann transport equation (D-LBTE) solvers have recently been developed, and one of the latest available software codes, Acuros XB, has been implemented in a commercial treatment planning system for radiotherapy photon beam dose calculation. One of the major limitations of most commercially available model-based algorithms for photon dose calculation is the ability to account for the effect of electron transport. This induces some errors in patient dose calculations, especially near heterogeneous interfaces between low and high density media such as tissue/lung interfaces. D-LBTE solvers have a high potential of producing accurate dose distributions in and near heterogeneous media in the human body. Extensive previous investigations have proved that D-LBTE solvers were able to produce comparable dose calculation accuracy as Monte Carlo methods with a reasonable speed good enough for clinical use. The current paper reviews the dosimetric evaluations of D-LBTE solvers for external beam photon radiotherapy. This content summarizes and discusses dosimetric validations for D-LBTE solvers in both homogeneous and heterogeneous media under different circumstances and also the clinical impact on various diseases due to the conversion of dose calculation from a conventional convolution/superposition algorithm to a recently released D-LBTE solver.

Radiochromic film based transit dosimetry for verification of dose delivery with intensity modulated radiotherapy

PURPOSE

To evaluate the transit dose based patient specific quality assurance (QA) of intensity modulated radiation therapy (IMRT) for verification of the accuracy of dose delivered to the patient.

METHODS

Five IMRT plans were selected and utilized to irradiate a homogeneous plastic water phantom and an inhomogeneous anthropomorphic phantom. The transit dose distribution was measured with radiochromic film and was compared with the computed dose map on the same plane using a gamma index with a 3% dose and a 3 mm distance-to-dose agreement tolerance limit.

RESULTS

While the average gamma index for comparisons of dose distributions was less than one for 98.9% of all pixels from the transit dose with the homogeneous phantom, the passing rate was reduced to 95.0% for the transit dose with the inhomogeneous phantom. Transit doses due to a 5 mm setup error may cause up to a 50% failure rate of the gamma index.

CONCLUSIONS

Transit dose based IMRT QA may be superior to the traditional QA method since the former can show whether the inhomogeneity correction algorithm from TPS is accurate. In addition, transit dose based IMRT QA can be used to verify the accuracy of the dose delivered to the patient during treatment by revealing significant increases in the failure rate of the gamma index resulting from errors in patient positioning during treatment.

Bone and mucosal dosimetry in skin radiation therapy: a Monte Carlo study using kilovoltage photon and megavoltage electron beams

This study examines variations of bone and mucosal doses with variable soft tissue and bone thicknesses, mimicking the oral or nasal cavity in skin radiation therapy. Monte Carlo simulations (EGSnrc-based codes) using the clinical kilovoltage (kVp) photon and megavoltage (MeV) electron beams, and the pencil-beam algorithm (Pinnacle(3) treatment planning system) using the MeV electron beams were performed in dose calculations. Phase-space files for the 105 and 220 kVp beams (Gulmay D3225 x-ray machine), and the 4 and 6 MeV electron beams (Varian 21 EX linear accelerator) with a field size of 5 cm diameter were generated using the BEAMnrc code, and verified using measurements. Inhomogeneous phantoms containing uniform water, bone and air layers were irradiated by the kVp photon and MeV electron beams. Relative depth, bone and mucosal doses were calculated for the uniform water and bone layers which were varied in thickness in the ranges of 0.5-2 cm and 0.2-1 cm. A uniform water layer of bolus with thickness equal to the depth of maximum dose (d(max)) of the electron beams (0.7 cm for 4 MeV and 1.5 cm for 6 MeV) was added on top of the phantom to ensure that the maximum dose was at the phantom surface. From our Monte Carlo results, the 4 and 6 MeV electron beams were found to produce insignificant bone and mucosal dose (<1%), when the uniform water layer at the phantom surface was thicker than 1.5 cm. When considering the 0.5 cm thin uniform water and bone layers, the 4 MeV electron beam deposited less bone and mucosal dose than the 6 MeV beam. Moreover, it was found that the 105 kVp beam produced more than twice the dose to bone than the 220 kVp beam when the uniform water thickness at the phantom surface was small (0.5 cm). However, the difference in bone dose enhancement between the 105 and 220 kVp beams became smaller when the thicknesses of the uniform water and bone layers in the phantom increased. Dose in the second bone layer interfacing with air was found to be higher for the 220 kVp beam than that of the 105 kVp beam, when the bone thickness was 1 cm. In this study, dose deviations of bone and mucosal layers of 18% and 17% were found between our results from Monte Carlo simulation and the pencil-beam algorithm, which overestimated the doses. Relative depth, bone and mucosal doses were studied by varying the beam nature, beam energy and thicknesses of the bone and uniform water using an inhomogeneous phantom to model the oral or nasal cavity. While the dose distribution in the pharynx region is unavailable due to the lack of a commercial treatment planning system commissioned for kVp beam planning in skin radiation therapy, our study provided an essential insight into the radiation staff to justify and estimate bone and mucosal dose.

The accuracy of dose calculations by anisotropic analytical algorithms for stereotactic radiotherapy in nasopharyngeal carcinoma

Nasopharyngeal tumors are commonly treated with intensity-modulated radiotherapy techniques. For photon dose calculations, problems related to loss of lateral electronic equilibrium exist when small fields are used. The anisotropic analytical algorithm (AAA) implemented in Varian Eclipse was developed to replace the pencil beam convolution (PBC) algorithm for more accurate dose prediction in an inhomogeneous medium. The purpose of this study was to investigate the accuracy of the AAA for predicting interface doses for intensity-modulated stereotactic radiotherapy boost of nasopharyngeal tumors. The central axis depth dose data and dose profiles of phantoms with rectangular air cavities for small fields were measured using a 6 MV beam. In addition, the air-tissue interface doses from six different intensity-modulated stereotactic radiotherapy plans were measured in an anthropomorphic phantom. The nasopharyngeal region of the phantom was especially modified to simulate the air cavities of a typical patient. The measured data were compared to the data calculated by both the AAA and the PBC algorithm. When using single small fields in rectangular air cavity phantoms, both AAA and PBC overestimated the central axis dose at and beyond the first few millimeters of the air-water interface. Although the AAA performs better than the PBC algorithm, its calculated interface dose could still be more than three times that of the measured dose when a 2 × 2 cm(2) field was used. Testing of the algorithms using the anthropomorphic phantom showed that the maximum overestimation by the PBC algorithm was 20.7%, while that by the AAA was 8.3%. When multiple fields were used in a patient geometry, the dose prediction errors of the AAA would be substantially reduced compared with those from a single field. However, overestimation of more than 3% could still be found at some points at the air-tissue interface.

In vivo real-time rectal wall dosimetry for prostate radiotherapy

Rectal balloons are used in external beam prostate radiotherapy to provide reproducible anatomy and rectal dose reductions. This is an investigation into the combination of a MOSFET radiation detector with a rectal balloon for realtime in vivo rectal wall dosimetry. The MOSFET used in the study is a radiation detector that provides a water equivalent depth of measurement of 70 microm. Two MOSFETs were combined in a face-to-face orientation. The reproducibility, sensitivity and angular dependence were measured for the dual MOSFET in a 6 MV photon beam. The dual MOSFET was combined with a rectal balloon and irradiated with hypothetical prostate treatments in a phantom. The anterior rectal wall dose was measured in real time and compared with the planning system calculated dose. The dual MOSFET showed angular dependence within +/-2.5% in the azimuth and +2.5%/-4% in the polar axes. When compared with an ion chamber measurement in a phantom, the dual MOSFET agreed within 2.5% for a range of radiation path lengths and incident angles. The dual MOSFET had reproducible sensitivity for fraction sizes of 2-10 Gy. For the hypothetical prostate treatments the measured anterior rectal wall dose was 2.6 and 3.2% lower than the calculated dose for 3DCRT and IMRT plans. This was expected due to limitations of the dose calculation method used at the balloon cavity interface. A dual MOSFET combined with a commercial rectal balloon was shown to provide reproducible measurements of the anterior rectal wall dose in real time. The measured anterior rectal wall dose agreed with the expected dose from the treatment plan for 3DCRT and IMRT plans. The dual MOSFET could be read out in real time during the irradiation, providing the capability for real-time dose monitoring of the rectal wall dose during treatment.

Stereotactic body radiation therapy: the report of AAPM Task Group 101

Task Group 101 of the AAPM has prepared this report for medical physicists, clinicians, and therapists in order to outline the best practice guidelines for the external-beam radiation therapy technique referred to as stereotactic body radiation therapy (SBRT). The task group report includes a review of the literature to identify reported clinical findings and expected outcomes for this treatment modality. Information is provided for establishing a SBRT program, including protocols, equipment, resources, and QA procedures. Additionally, suggestions for developing consistent documentation for prescribing, reporting, and recording SBRT treatment delivery is provided.

Monte Carlo vs. pencil beam based optimization of stereotactic lung IMRT

Background

The purpose of the present study is to compare finite size pencil beam (fsPB) and Monte Carlo (MC) based optimization of lung intensity-modulated stereotactic radiotherapy (lung IMSRT).

Materials and methods

A fsPB and a MC algorithm as implemented in a biological IMRT planning system were validated by film measurements in a static lung phantom. Then, they were applied for static lung IMSRT planning based on three different geometrical patient models (one phase static CT, density overwrite one phase static CT, average CT) of the same patient. Both 6 and 15 MV beam energies were used. The resulting treatment plans were compared by how well they fulfilled the prescribed optimization constraints both for the dose distributions calculated on the static patient models and for the accumulated dose, recalculated with MC on each of 8 CTs of a 4DCT set.

Results

In the phantom measurements, the MC dose engine showed discrepancies < 2%, while the fsPB dose engine showed discrepancies of up to 8% in the presence of lateral electron disequilibrium in the target. In the patient plan optimization, this translates into violations of organ at risk constraints and unpredictable target doses for the fsPB optimized plans. For the 4D MC recalculated dose distribution, MC optimized plans always underestimate the target doses, but the organ at risk doses were comparable. The results depend on the static patient model, and the smallest discrepancy was found for the MC optimized plan on the density overwrite one phase static CT model.

Conclusions

It is feasible to employ the MC dose engine for optimization of lung IMSRT and the plans are superior to fsPB. Use of static patient models introduces a bias in the MC dose distribution compared to the 4D MC recalculated dose, but this bias is predictable and therefore MC based optimization on static patient models is considered safe.

Breast irradiation with three conformal photon fields for patients with high lung involvement

Since 2001, 50 breast cancer patients, for whom extensive lung/heart involvement was expected from the use of conventional tangential 2-field technique (2F) owing to complex anatomies, were irradiated using a 3-field conformal technique (3F). Dose plans were designed for both 3F and 2F and a dose volume histogram analysis on ipsilateral lung, heart, and target was conducted. The 3F technique allowed a significant reduction in ipsilateral lung irradiation: mean dose dropped from 16.0+/-3.8 (2F) to 12.0+/-2.7 Gy (3F) and V(45Gy) from 20.7+/-6.8 (2F) to 3.2+/-2.9% (3F). Similarly, in patients irradiated to the left breast, heart mean dose was reduced from 8.1 Gy (2F) to 6.8 Gy (3F) and D(15%) from 19.0 Gy to 14.0 Gy. All differences reached a high degree of significance. The target coverage was not clinically compromised since the slight reduction using 3F compared with 2F is limited to V(95%) while V(90%) difference, even if statistically significant, is small: 98.2+/-1.8% (3F) and 98.8+/-1.6 (2F). A preliminary report on clinical follow up is also included; with a mean follow up of 15.8 months, no pulmonary or cardiac complications were observed.

Endo-rectal balloon cavity dosimetry in a phantom: performance under IMRT and helical tomotherapy beams

BACKGROUND AND PURPOSE

The use of endo-rectal balloons as immobilisation devices in external beam radiotherapy for prostate cancer has led to improved target position reproducibility and a decrease in rectal toxicity. The air cavity created by an endo-rectal balloon in photon radiotherapy perturbs the dose distribution. In this study, the effect of the balloon cavity on the dose distribution and the accuracy to which two treatment planning systems calculate the dose distribution were investigated.

MATERIALS AND METHODS

Single beams as well as 3D conformal, conventional IMRT and helical tomotherapy treatment plans were investigated using a specifically constructed phantom. Radiochromic film was used to measure the cavity wall doses and cavity wall DVHs.

RESULTS

For a 70 Gy prescription dose both the Pinnacle and TomoTherapy TPSs over-predicted the anterior cavity wall dose by 1.43 Gy, 3.92 Gy and 2.67 Gy for 3D conformal, conventional IMRT and helical tomotherapy, respectively. The posterior cavity wall dose was under-predicted by 2.62 Gy, 2.01 Gy and 4.79 Gy for 3D conformal, conventional IMRT and helical tomotherapy, respectively. An over-prediction by the Pinnacle RTPS of the V50, V60, V65 and V70 values for the cavity wall DVH was measured for the 3D conformal and conventional IMRT cases. These reductions may lead to a less than expected rectal toxicity. The TomoTherapy RTPS under-predicted the V50, V60, V65 and V70 values which may lead to higher rectal toxicity than predicted.

CONCLUSION

Calculation of dose around an air cavity created by an endo-rectal balloon provides a challenge for radiotherapy planning systems. Various electronic disequilibrium situations exist due to the cavity, which can lead to a lower anterior rectal wall and higher posterior rectal wall dose than that calculated by planning systems. This has consequences for comparisons of dose volume constraints between different modalities.

Feasibility of a multigroup deterministic solution method for three-dimensional radiotherapy dose calculations

PURPOSE

To investigate the potential of a novel deterministic solver, Attila, for external photon beam radiotherapy dose calculations.

METHODS AND MATERIALS

Two hypothetical cases for prostate and head-and-neck cancer photon beam treatment plans were calculated using Attila and EGSnrc Monte Carlo simulations. Open beams were modeled as isotropic photon point sources collimated to specified field sizes. The sources had a realistic energy spectrum calculated by Monte Carlo for a Varian Clinac 2100 operated in a 6-MV photon mode. The Attila computational grids consisted of 106,000 elements, or 424,000 spatial degrees of freedom, for the prostate case, and 123,000 tetrahedral elements, or 492,000 spatial degrees of freedom, for the head-and-neck cases.

RESULTS

For both cases, results demonstrate excellent agreement between Attila and EGSnrc in all areas, including the build-up regions, near heterogeneities, and at the beam penumbra. Dose agreement for 99% of the voxels was within the 3% (relative point-wise difference) or 3-mm distance-to-agreement criterion. Localized differences between the Attila and EGSnrc results were observed at bone and soft-tissue interfaces and are attributable to the effect of voxel material homogenization in calculating dose-to-medium in EGSnrc. For both cases, Attila calculation times were <20 central processing unit minutes on a single 2.2-GHz AMD Opteron processor.

CONCLUSIONS

The methods in Attila have the potential to be the basis for an efficient dose engine for patient-specific treatment planning, providing accuracy similar to that obtained by Monte Carlo.

Computational Challenges for Image-Guided Radiation Therapy: Framework and Current Research

It is arguable that the imaging and delivery hardware necessary for delivering real-time adaptive image-guided radiotherapy is available on high-end linear accelerators. Robust and computationally efficient software is the limiting factor in achieving highly accurate and precise radiotherapy to the constantly changing anatomy of a cancer patient. The limitations are not caused by the availability of algorithms but rather issues of reliability, integration, and calculation time. However, each of the software components is an active area of research and development at academic and commercial centers. This article describes the software solutions in 4 broad areas: deformable image registration, adaptive replanning, real-time image guidance, and dose calculation and accumulation. Given the pace of technological advancement, the integration of these software solutions to develop real-time adaptive image-guided radiotherapy and the associated challenges they bring will be implemented to varying degrees by all major manufacturers over the coming years.

Monte carlo evaluation of the AAA treatment planning algorithm in a heterogeneous multilayer phantom and IMRT clinical treatments for an Elekta SL25 linear accelerator

The Anisotropic Analytical Algorithm (AAA) is a new pencil beam convolution/superposition algorithm proposed by Varian for photon dose calculations. The configuration of AAA depends on linear accelerator design and specifications. The purpose of this study was to investigate the accuracy of AAA for an Elekta SL25 linear accelerator for small fields and intensity modulated radiation therapy (IMRT) treatments in inhomogeneous media. The accuracy of AAA was evaluated in two studies. First, AAA was compared both with Monte Carlo (MC) and the measurements in an inhomogeneous phantom simulating lung equivalent tissues and bone ribs. The algorithm was tested under lateral electronic disequilibrium conditions, using small fields (2 x 2 cm(2)). Good agreement was generally achieved for depth dose and profiles, with deviations generally below 3% in lung inhomogeneities and below 5% at interfaces. However, the effects of attenuation and scattering close to the bone ribs were not fully taken into account by AAA, and small inhomogeneities may lead to planning errors. Second, AAA and MC were compared for IMRT plans in clinical conditions, i.e., dose calculations in a computed tomography scan of a patient. One ethmoid tumor, one orophaxynx and two lung tumors are presented in this paper. Small differences were found between the dose volume histograms. For instance, a 1.7% difference for the mean planning target volume dose was obtained for the ethmoid case. Since better agreement was achieved for the same plans but in homogeneous conditions, these differences must be attributed to the handling of inhomogeneities by AAA. Therefore, inherent assumptions of the algorithm, principally the assumption of independent depth and lateral directions in the scaling of the kernels, were slightly influencing AAA's validity in inhomogeneities. However, AAA showed a good accuracy overall and a great ability to handle small fields in inhomogeneous media compared to other pencil beam convolution algorithms.

Accuracy of inhomogeneity correction algorithm in intensity-modulated radiotherapy of head-and-neck tumors

We examined the degree of calculated-to-measured dose difference for nasopharyngeal target volume in intensity-modulated radiotherapy (IMRT) based on the observed/expected ratio using patient anatomy with humanoid head-and-neck phantom. The plans were designed with a clinical treatment planning system that uses a measurement-based pencil beam dose-calculation algorithm. Two kinds of IMRT plans, which give a direct indication of the error introduced in routine treatment planning, were categorized and evaluated. The experimental results show that when the beams pass through the oral cavity in anthropomorphic head-and-neck phantom, the average dose difference becomes significant, revealing about 10% dose difference to prescribed dose at isocenter. To investigate both the physical reasons of the dose discrepancy and the inhomogeneity effect, we performed the 10 cases of IMRT quality assurance (QA) with plastic and humanoid phantoms. Our result suggests that the transient electronic disequilibrium with the increased lateral electron range may cause the inaccuracy of dose calculation algorithm, and the effectiveness of the inhomogeneity corrections used in IMRT plans should be evaluated to ensure meaningful quality assurance and delivery.

A finite size pencil beam algorithm for IMRT dose optimization: density corrections

For beamlet-based IMRT optimization, fast and less accurate dose computation algorithms are frequently used, while more accurate algorithms are needed to recompute the final dose for verification. In order to speed up the optimization process and ensure close proximity between dose in optimization and verification, proper consideration of dose gradients and tissue inhomogeneity effects should be ensured at every stage of the optimization. Due to their speed, pencil beam algorithms are often used for precalculation of beamlet dose distributions in IMRT treatment planning systems. However, accounting for tissue heterogeneities with these models requires the use of approximate rescaling methods. Recently, a finite size pencil beam (fsPB) algorithm, based on a simple and small set of data, was proposed which was specifically designed for the purpose of dose pre-computation in beamlet-based IMRT. The present work describes the incorporation of 3D density corrections, based on Monte Carlo simulations in heterogeneous phantoms, into this method improving the algorithm accuracy in inhomogeneous geometries while keeping its original speed and simplicity of commissioning. The algorithm affords the full accuracy of 3D density corrections at every stage of the optimization, hence providing the means for density related fluence modulation like penumbra shaping at field edges.

Comparison of 6MV and 18MV photons for IMRT treatment of lung cancer

BACKGROUND AND PURPOSE

To compare 6 MV and 18 MV photon intensity modulated radiotherapy (IMRT) for non-small cell lung cancer.

MATERIALS AND METHODS

Doses for a cohort of 10 patients, typical for our department, were computed with a commercially available convolution/superposition (CS) algorithm. Final dose computation was also performed with a dedicated IMRT Monte Carlo dose engine (MCDE).

RESULTS

CS plans showed higher D(95%) (Gy) for the GTV (68.13 vs 67.36, p=0.004) and CTV (67.23 vs 66.87, p=0.028) with 18 than with 6 MV photons. MCDE computations demonstrated higher doses with 6 MV than 18 MV in D(95%) for the PTV (64.62 vs 63.64, p=0.009), PTV(optim) (65.48 vs 64.83, p=0.014) and CTV (66.22 vs 65.64, p=0.027). Dose inhomogeneity was lower with 18 than with 6 MV photons for GTV (0.08 vs 0.09, p=0.007) and CTV (0.10 vs 0.11, p=0.045) in CS but not MCDE plans. 6 MV photons significantly (D(33%); p=0.045) spared the esophagus in MCDE plans. Observed dose differences between lower and higher energy IMRT plans were dependent on the individual patient.

CONCLUSIONS

Selection of photon energy depends on priority ranking of endpoints and individual patients. In the absence of highly accurate dose computation algorithms such as CS and MCDE, 6 MV photons may be the prudent choice.

Optimization of extracranial stereotactic radiation therapy of small lung lesions using accurate dose calculation algorithms

BACKGROUND

The aim of this study was to compare and to validate different dose calculation algorithms for the use in radiation therapy of small lung lesions and to optimize the treatment planning using accurate dose calculation algorithms.

METHODS

A 9-field conformal treatment plan was generated on an inhomogeneous phantom with lung mimics and a soft tissue equivalent insert, mimicking a lung tumor. The dose distribution was calculated with the Pencil Beam and Collapsed Cone algorithms implemented in Masterplan (Nucletron) and the Monte Carlo system XVMC and validated using Gafchromic EBT films. Differences in dose distribution were evaluated. The plans were then optimized by adding segments to the outer shell of the target in order to increase the dose near the interface to the lung.

RESULTS

The Pencil Beam algorithm overestimated the dose by up to 15% compared to the measurements. Collapsed Cone and Monte Carlo predicted the dose more accurately with a maximum difference of -8% and -3% respectively compared to the film. Plan optimization by adding small segments to the peripheral parts of the target, creating a 2-step fluence modulation, allowed to increase target coverage and homogeneity as compared to the uncorrected 9 field plan.

CONCLUSION

The use of forward 2-step fluence modulation in radiotherapy of small lung lesions allows the improvement of tumor coverage and dose homogeneity as compared to non-modulated treatment plans and may thus help to increase the local tumor control probability. While the Collapsed Cone algorithm is closer to measurements than the Pencil Beam algorithm, both algorithms are limited at tissue/lung interfaces, leaving Monte-Carlo the most accurate algorithm for dose prediction.

EPR dosimetry of radiotherapy photon beams in inhomogeneous media using alanine films

In the current work, EPR (electron paramagnetic resonance) dosimetry using alanine films (134 microm thick) was utilized for dose measurements in inhomogeneous phantoms irradiated with radiotherapy photon beams. The main phantom material was PMMA, while either Styrofoam or aluminium was introduced as an inhomogeneity. The phantoms were irradiated to a maximum dose of about 30 Gy with 6 or 15 MV photons. The performance of the alanine film dosimeters was investigated and compared to results from ion chamber dosimetry, Monte Carlo simulations and radiotherapy treatment planning calculations. It was found that the alanine film dosimeters had a linear dose response above approximately 5 Gy, while a background signal obscured the response at lower dose levels. For doses between 5 and 60 Gy, the standard deviation of single alanine film dose estimates was about 2%. The alanine film dose estimates yielded results comparable to those from the Monte Carlo simulations and the ion chamber measurements, with absolute differences between estimates in the order of 1-15%. The treatment planning calculations exhibited limited applicability. The current work shows that alanine film dosimetry is a method suitable for estimating radiotherapeutical doses and for dose measurements in inhomogeneous media.

Comparison of dose-volume histograms of IMRT treatment plans for ethmoid sinus cancer computed by advanced treatment planning systems including Monte Carlo

BACKGROUND AND PURPOSE

To recompute clinical intensity-modulated treatment plans for ethmoid sinus cancer and to compare quantitatively the dose-volume histograms (DVHs) of the planning target volume (PTV) and the optic organs at risk.

MATERIAL AND METHODS

Ten step-and-shoot intensity-modulated treatment plans were enrolled in this study. Large natural and surgical air cavities challenged the calculation systems. Each optimized treatment plan was recalculated by two superposition convolution (TMS and Pinnacle) and a Monte Carlo system (MCDE). To compare the resulting DVHs, a one-way ANOVA for repeated measurements was performed and multiple pairwise comparisons were made.

RESULTS

The tails of the PTV-DVHs were significantly higher for the Monte Carlo system. The DVHs of the critical organs displayed some statistically but not always clinically significant differences. For the individual patients, the three planning systems sometimes reproduced clinically discrepant DVHs that were not significantly different when averaged over all patients.

CONCLUSIONS

Dose to air cavities contains computational uncertainty. As this dose is clinically irrelevant and optimizing it is meaningless, we recommended extracting the air from the PTV when constructing the PTV-DVH. The planning systems considered reproduce DVHs that are significantly different, especially in the tail region of PTV-DVHs.

Accuracy of patient dose calculation for lung IMRT: A comparison of Monte Carlo, convolution/superposition, and pencil beam computations

The accuracy of dose computation within the lungs depends strongly on the performance of the calculation algorithm in regions of electronic disequilibrium that arise near tissue inhomogeneities with large density variations. There is a lack of data evaluating the performance of highly developed analytical dose calculation algorithms compared to Monte Carlo computations in a clinical setting. We compared full Monte Carlo calculations (performed by our Monte Carlo dose engine MCDE) with two different commercial convolution/superposition (CS) implementations (Pinnacle-CS and Helax-TMS's collapsed cone model Helax-CC) and one pencil beam algorithm (Helax-TMS's pencil beam model Helax-PB) for 10 intensity modulated radiation therapy (IMRT) lung cancer patients. Treatment plans were created for two photon beam qualities (6 and 18 MV). For each dose calculation algorithm, patient, and beam quality, the following set of clinically relevant dose-volume values was reported: (i) minimal, median, and maximal dose (Dmin, D50, and Dmax) for the gross tumor and planning target volumes (GTV and PTV); (ii) the volume of the lungs (excluding the GTV) receiving at least 20 and 30 Gy (V20 and V30) and the mean lung dose; (iii) the 33rd percentile dose (D33) and Dmax delivered to the heart and the expanded esophagus; and (iv) Dmax for the expanded spinal cord. Statistical analysis was performed by means of one-way analysis of variance for repeated measurements and Tukey pairwise comparison of means. Pinnacle-CS showed an excellent agreement with MCDE within the target structures, whereas the best correspondence for the organs at risk (OARs) was found between Helax-CC and MCDE. Results from Helax-PB were unsatisfying for both targets and OARs. Additionally, individual patient results were analyzed. Within the target structures, deviations above 5% were found in one patient for the comparison of MCDE and Helax-CC, while all differences between MCDE and Pinnacle-CS were below 5%. For both Pinnacle-CS and Helax-CC, deviations from MCDE above 5% were found within the OARs: within the lungs for two (6 MV) and six (18 MV) patients for Pinnacle-CS, and within other OARs for two patients for Helax-CC (for Dmax of the heart and D33 of the expanded esophagus) but only for 6 MV. For one patient, all four algorithms were used to recompute the dose after replacing all computed tomography voxels within the patient's skin contour by water. This made all differences above 5% between MCDE and the other dose calculation algorithms disappear. Thus, the observed deviations mainly arose from differences in particle transport modeling within the lungs, and the commissioning of the algorithms was adequately performed (or the commissioning was less important for this type of treatment). In conclusion, not one pair of the dose calculation algorithms we investigated could provide results that were consistent within 5% for all 10 patients for the set of clinically relevant dose-volume indices studied. As the results from both CS algorithms differed significantly, care should be taken when evaluating treatment plans as the choice of dose calculation algorithm may influence clinical results. Full Monte Carlo provides a great benchmarking tool for evaluating the performance of other algorithms for patient dose computations.

Dose build-up behind air cavities for Co-60, 4, 6 and 8 MV. Measurements and Monte Carlo simulations

It has been shown in several studies that the build-up in photon beams behind air cavities (such as in the head and neck) increases with energy. In this study this effect is investigated over a broad range of energies that have been used for treating head and neck tumours. The study addresses the question of whether an energy lower than 6 MV is desirable and is based on measurements and Monte Carlo (MC) simulations. In a PMMA phantom containing an air cavity (3 x 16 x 3 cm3 at 3 cm depth) an ionization chamber (Capintec PS-033) was used to measure the dose build-up behind the cavity for 4, 6 and 8 MV beam qualities for different field sizes (from 3 x 6 cm2 to 8 x 8 cm2). MC simulations were made using the EGSnrc code for the same geometry and energies as well as for Co-60. Measurements and MC simulations agree well when the fixed-separation plane-parallel chamber measurements have been corrected for the expected over-response in the build-up region. This work demonstrates that the build-up effect of 6 MV is 'closer' to the build-up effect of 8 MV than to that of 4 MV. This suggests that if the build-up effect is of concern when the target volume is in the vicinity of air cavities, 4 MV should be preferred over both 6 MV and 8 MV. This work also shows that the build-up effect for Co-60 is significantly smaller than that of 4 MV. Moreover, the build-up effect increases as the field size decreases. With the increasing use of IMRT (and radiosurgery), small fields are used more frequently making these issues even more relevant. This should be taken into consideration when choosing the accelerator energies for a radiotherapy department.

Comparison of PENELOPE Monte Carlo dose calculations with Fricke dosimeter and ionization chamber measurements in heterogeneous phantoms (18 MeV electron and 12 MV photon beams)

Different measurements of depth-dose curves and dose profiles were performed in heterogeneous phantoms and compared to dose distributions calculated by a Monte Carlo code. These heterogeneous phantoms consisted of lung and/or bone heterogeneities. Irradiations and simulations were carried out for an 18 MeV electron beam and a 12 MV photon beam. Depth-dose curves were measured with Fricke dosimeters and with plane and cylindrical ionization chambers. Dose profiles were measured with a small cylindrical ionization chamber at different depths. The LINAC was modelled using the PENELOPE code and phase space files were used as input data for the calculations of the dose distributions in every simulation. The detectors (Fricke dosimeters and ionization chambers) were not modelled in the geometry. There is generally a good agreement between the measurements and PENELOPE. Some discrepancies exist, near interfaces, between the ionization chamber and PENELOPE due to the attenuation of the lower energy electrons by the wall of the ionization chamber.