Monte Carlo modelling of external radiotherapy photon beams

Ant colony algorithm implementation in electron and photon Monte Carlo transport: application to the commissioning of radiosurgery photon beams

PURPOSE

In this work, the authors describe an approach which has been developed to drive the application of different variance-reduction techniques to the Monte Carlo simulation of photon and electron transport in clinical accelerators.

METHODS

The new approach considers the following techniques: Russian roulette, splitting, a modified version of the directional bremsstrahlung splitting, and the azimuthal particle redistribution. Their application is controlled by an ant colony algorithm based on an importance map.

RESULTS

The procedure has been applied to radiosurgery beams. Specifically, the authors have calculated depth-dose profiles, off-axis ratios, and output factors, quantities usually considered in the commissioning of these beams. The agreement between Monte Carlo results and the corresponding measurements is within approximately 3%/0.3 mm for the central axis percentage depth dose and the dose profiles. The importance map generated in the calculation can be used to discuss simulation details in the different parts of the geometry in a simple way. The simulation CPU times are comparable to those needed within other approaches common in this field.

CONCLUSIONS

The new approach is competitive with those previously used in this kind of problems (PSF generation or source models) and has some practical advantages that make it to be a good tool to simulate the radiation transport in problems where the quantities of interest are difficult to obtain because of low statistics.

Four-dimensional dosimetry validation and study in lung radiotherapy using deformable image registration and Monte Carlo techniques

Thoracic cancer treatment presents dosimetric difficulties due to respiratory motion and lung inhomogeneity. Monte Carlo and deformable image registration techniques have been proposed to be used in four-dimensional (4D) dose calculations to overcome the difficulties. This study validates the 4D Monte Carlo dosimetry with measurement, compares 4D dosimetry of different tumor sizes and tumor motion ranges, and demonstrates differences of dose-volume histograms (DVH) with the number of respiratory phases that are included in 4D dosimetry. BEAMnrc was used in dose calculations while an optical flow algorithm was used in deformable image registration and dose mapping. Calculated and measured doses of a moving phantom agreed within 3% at the center of the moving gross tumor volumes (GTV). 4D CT image sets of lung cancer cases were used in the analysis of 4D dosimetry. For a small tumor (12.5 cm³) with motion range of 1.5 cm, reduced tumor volume coverage was observed in the 4D dose with a beam margin of 1 cm. For large tumors and tumors with small motion range (around 1 cm), the 4D dosimetry did not differ appreciably from the static plans. The dose-volume histogram (DVH) analysis shows that the inclusion of only extreme respiratory phases in 4D dosimetry is a reasonable approximation of all-phase inclusion for lung cancer cases similar to the ones studied, which reduces the calculation in 4D dosimetry.

Measurement of depth-dose of linear accelerator and simulation by use of Geant4 computer code

Radiation therapy is an established method of cancer treatment. New technologies in cancer radiotherapy need a more accurate computation of the dose delivered in the radiotherapy treatment plan. This study presents some results of a Geant4-based application for simulation of the absorbed dose distribution given by a medical linear accelerator (LINAC). The LINAC geometry is accurately described in the Monte Carlo code with use of the accelerator manufacturer's specifications. The capability of the software for evaluating the dose distribution has been verified by comparisons with measurements in a water phantom; the comparisons were performed for percentage depth dose (PDD) and profiles for various field sizes and depths, for a 6-MV electron beam. Experimental and calculated dose values were in good agreement both in PDD and in transverse sections of the water phantom.

Implementation of angular response function modeling in SPECT simulations with GATE

Among Monte Carlo simulation codes in medical imaging, the GATE simulation platform is widely used today given its flexibility and accuracy, despite long run times, which in SPECT simulations are mostly spent in tracking photons through the collimators. In this work, a tabulated model of the collimator/detector response was implemented within the GATE framework to significantly reduce the simulation times in SPECT. This implementation uses the angular response function (ARF) model. The performance of the implemented ARF approach has been compared to standard SPECT GATE simulations in terms of the ARF tables' accuracy, overall SPECT system performance and run times. Considering the simulation of the Siemens Symbia T SPECT system using high-energy collimators, differences of less than 1% were measured between the ARF-based and the standard GATE-based simulations, while considering the same noise level in the projections, acceleration factors of up to 180 were obtained when simulating a planar 364 keV source seen with the same SPECT system. The ARF-based and the standard GATE simulation results also agreed very well when considering a four-head SPECT simulation of a realistic Jaszczak phantom filled with iodine-131, with a resulting acceleration factor of 100. In conclusion, the implementation of an ARF-based model of collimator/detector response for SPECT simulations within GATE significantly reduces the simulation run times without compromising accuracy.

Accurate accumulation of dose for improved understanding of radiation effects in normal tissue

The actual distribution of radiation dose accumulated in normal tissues over the complete course of radiation therapy is, in general, poorly quantified. Differences in the patient anatomy between planning and treatment can occur gradually (e.g., tumor regression, resolution of edema) or relatively rapidly (e.g., bladder filling, breathing motion) and these undermine the accuracy of the planned dose distribution. Current efforts to maximize the therapeutic ratio require models that relate the true accumulated dose to clinical outcome. The needed accuracy can only be achieved through the development of robust methods that track the accumulation of dose within the various tissues in the body. Specific needs include the development of segmentation methods, tissue-mapping algorithms, uncertainty estimation, optimal schedules for image-based monitoring, and the development of informatics tools to support subsequent analysis. These developments will not only improve radiation outcomes modeling but will address the technical demands of the adaptive radiotherapy paradigm. The next 5 years need to see academia and industry bring these tools into the hands of the clinician and the clinical scientist.

Fast, accurate photon beam accelerator modeling using BEAMnrc: a systematic investigation of efficiency enhancing methods and cross-section data

In this work, an investigation of efficiency enhancing methods and cross-section data in the BEAMnrc Monte Carlo (MC) code system is presented. Additionally, BEAMnrc was compared with VMC++, another special-purpose MC code system that has recently been enhanced for the simulation of the entire treatment head. BEAMnrc and VMC++ were used to simulate a 6 MV photon beam from a Siemens Primus linear accelerator (linac) and phase space (PHSP) files were generated at 100 cm source-to-surface distance for the 10 x 10 and 40 x 40 cm2 field sizes. The BEAMnrc parameters/techniques under investigation were grouped by (i) photon and bremsstrahlung cross sections, (ii) approximate efficiency improving techniques (AEITs), (iii) variance reduction techniques (VRTs), and (iv) a VRT (bremsstrahlung photon splitting) in combination with an AEIT (charged particle range rejection). The BEAMnrc PHSP file obtained without the efficiency enhancing techniques under study or, when not possible, with their default values (e.g., EXACT algorithm for the boundary crossing algorithm) and with the default cross-section data (PEGS4 and Bethe-Heitler) was used as the "base line" for accuracy verification of the PHSP files generated from the different groups described previously. Subsequently, a selection of the PHSP files was used as input for DOSXYZnrc-based water phantom dose calculations, which were verified against measurements. The performance of the different VRTs and AEITs available in BEAMnrc and of VMC++ was specified by the relative efficiency, i.e., by the efficiency of the MC simulation relative to that of the BEAMnrc base-line calculation. The highest relative efficiencies were approximately 935 (approximately 111 min on a single 2.6 GHz processor) and approximately 200 (approximately 45 min on a single processor) for the 10 x 10 field size with 50 million histories and 40 x 40 cm2 field size with 100 million histories, respectively, using the VRT directional bremsstrahlung splitting (DBS) with no electron splitting. When DBS was used with electron splitting and combined with augmented charged particle range rejection, a technique recently introduced in BEAMnrc, relative efficiencies were approximately 420 (approximately 253 min on a single processor) and approximately 175 (approximately 58 min on a single processor) for the 10 x 10 and 40 x 40 cm2 field sizes, respectively. Calculations of the Siemens Primus treatment head with VMC++ produced relative efficiencies of approximately 1400 (approximately 6 min on a single processor) and approximately 60 (approximately 4 min on a single processor) for the 10 x 10 and 40 x 40 cm2 field sizes, respectively. BEAMnrc PHSP calculations with DBS alone or DBS in combination with charged particle range rejection were more efficient than the other efficiency enhancing techniques used. Using VMC++, accurate simulations of the entire linac treatment head were performed within minutes on a single processor. Noteworthy differences (+/- 1%-3%) in the mean energy, planar fluence, and angular and spectral distributions were observed with the NIST bremsstrahlung cross sections compared with those of Bethe-Heitler (BEAMnrc default bremsstrahlung cross section). However, MC calculated dose distributions in water phantoms (using combinations of VRTs/AEITs and cross-section data) agreed within 2% of measurements. Furthermore, MC calculated dose distributions in a simulated water/air/water phantom, using NIST cross sections, were within 2% agreement with the BEAMnrc Bethe-Heitler default case.

Simulation of large x-ray fields using independently measured source and geometry details

PURPOSE

Obtain an accurate simulation of the dose from the 6 and 18 MV x-ray beams from a Siemens Oncor linear accelerator by comparing simulation to measurement. Constrain the simulation by independently determining parameters of the treatment head and incident beam, in particular, the energy and spot size.

METHODS

Measurements were done with the treatment head in three different configurations: (1) The clinical configuration, (2) the flattening filter removed, and (3) the target and flattening filter removed. Parameters of the incident beam and treatment head were measured directly. Incident beam energy and spectral width were determined from the percent-depth ionization of the raw beam (as described previously), spot size was determined using a spot camera, and the densities of the flattening filters were determined by weighing them. Simulations were done with EGSnrc/BEAMnrc code. An asymmetric simulation was used, including offsets of the spot, primary collimator, and flattening filter from the collimator rotation axis.

RESULTS

Agreement between measurement and simulation was obtained to the least restrictive of 1% or 1 mm at 6 MV, both with and without the flattening filter in place, except for the buildup region. At 18 MV, the agreement was 1.5%/1.5 mm with the flattening filter in place and 1%/1 mm with it removed, except for in the buildup region. In the buildup region, the discrepancy was 2%/2 mm at 18 MV and 1.5%/1.5 mm at 6 MV with the flattening filter either removed or in place. The methodology for measuring the source and geometry parameters for the treatment head simulation is described. Except to determine the density of the flattening filter, no physical modification of the treatment head is necessary to obtain those parameters. In particular, the flattening filter does not need to be removed as was done in this work.

CONCLUSIONS

Good agreement between measured and simulated dose distributions was obtained, even in the buildup region. The simulation was tightly constrained by independent measurements of parameters of the incident beam and treatment head. The method of obtaining the input parameters is described, and can be carried out on a clinical linear accelerator.

A virtual source model of electron contamination of a therapeutic photon beam

The most efficient way of generating particles for Monte Carlo (MC) dose calculation is through a virtual source model (VSM) of the linear accelerator head. We have previously developed a VSM based on three sources: a primary photon source, a secondary photon source and an electron contamination source (Sikora et al 2007). In this work, we present an improvement of the electron contamination source. The VSM of contamination electrons (eVSM) is derived from a full MC simulation of the accelerator head with the BEAMnrc MC system. It comprises a Gaussian source located at the base of the flattening filter. The eVSM models two effects: an energy-dependent source diameter and an angular dependence of the particle fluence. The air scatter of the contamination electrons is approximated by energetic properties of the eVSM so that explicit in-air transport is not required during MC simulation of the dose distributions in the patient. The calculations of electron dose distributions were compared between the eVSM and the full MC simulation. Good agreement was achieved for various rectangular field sizes as well as for complex conformal segment shapes for the contamination electrons of 6 and 15 MV beams. The 3D dose evaluation of the surface dose in a CT-based patient geometry shows high accuracy (2%/2 mm) of the eVSM for both energies. The model has one tunable parameter, the mean energy of the spectrum at the patient surface. High accuracy and efficiency of particle generation make the eVSM a valuable virtual source of contamination electrons for MC treatment planning systems.

The accuracy of EGSnrc, Geant4 and PENELOPE Monte Carlo systems for the simulation of electron scatter in external beam radiotherapy

Three widely used Monte Carlo systems were benchmarked against recently published measurements of the angular distribution of 13 MeV and 20 MeV electrons scattered from foils of different atomic numbers and thicknesses. Source and geometry were simulated in detail to calculate electron fluence profiles 118.2 cm from the exit window. Results were compared to the measured fluence profiles and the characteristic angle where the fluence drops to 1/e of its maximum value. EGSnrc and PENELOPE results, on average, agreed with measurement within 1 standard deviation experimental uncertainty, with EGSnrc estimating slightly lower scatter than measurement and PENELOPE slightly higher scatter. Geant4.9.2 overestimated the characteristic angle for the lower atomic number foils by as much as 10%. Retuning of the scatter distributions in Geant4 led to a much better agreement with measurement, close to that achieved with the other codes. The 3% differences from measurement seen with all codes for at least some of the foils would result in clinically significant errors in the fluence profiles (2%/4 mm), given accurate knowledge of the electron source and treatment head geometry used in radiotherapy. Further improvement in simulation accuracy is needed to achieve 1%/1 mm agreement with measurement for the full range of beam energies, foil atomic number and thickness used in radiotherapy. EGSnrc would achieve this accuracy with an increase in thickness of the mylar sheets in the monitor chamber, PENELOPE with a decrease in thickness.

A comprehensive system for dosimetric commissioning and Monte Carlo validation for the small animal radiation research platform

Our group has constructed the small animal radiation research platform (SARRP) for delivering focal, kilo-voltage radiation to targets in small animals under robotic control using cone-beam CT guidance. The present work was undertaken to support the SARRP's treatment planning capabilities. We have devised a comprehensive system for characterizing the radiation dosimetry in water for the SARRP and have developed a Monte Carlo dose engine with the intent of reproducing these measured results. We find that the SARRP provides sufficient therapeutic dose rates ranging from 102 to 228 cGy min(-1) at 1 cm depth for the available set of high-precision beams ranging from 0.5 to 5 mm in size. In terms of depth-dose, the mean of the absolute percentage differences between the Monte Carlo calculations and measurement is 3.4% over the full range of sampled depths spanning 0.5-7.2 cm for the 3 and 5 mm beams. The measured and computed profiles for these beams agree well overall; of note, good agreement is observed in the profile tails. Especially for the smallest 0.5 and 1 mm beams, including a more realistic description of the effective x-ray source into the Monte Carlo model may be important.

Inference of the optimal pretarget electron beam parameters in a Monte Carlo virtual linac model through simulated annealing

The purpose of this study was to develop an efficient method to determine the optimal intensity distribution of the pretarget electron beam in a Monte Carlo (MC) accelerator model able to most accurately reproduce a set of measured photon field profiles for a given accelerator geometry and nominal photon beam energy. The method has the ability to reduce the number of simulations required to commission a MC accelerator model and has achieved better agreement with measurement than other methods described in literature. The method begins from a cylindrically symmetric pretarget electron beam (radius of 0.5 cm) of uniform intensity. This beam is subdivided into annular regions of fluence for which each region is individually transported through the accelerator head and into a water phantom. A simulated annealing search is then performed to determine the optimal combination of weights of the annular fluences that provide a best match between the measured dose distributions and the weighted sum of annular dose distributions for particular pretarget electron energy. When restricted to Gaussian intensity distributions, the optimization determined an optimal FWHM=1.34 mm for 18.0 MeV electrons, with a RMSE=0.49% on 40 x 40 cm2 lateral profiles. When allowed to deviate from Gaussian intensities a further reduction in RMSE was achieved. For our Clinac 21 EX accelerator MC model (based on the 1996 Varian Oncology Systems, Monte Carlo Project package), the optimal unrestricted intensity distribution was found to be a Gaussian-like solution (18.0 MeV, FWHM= 1.10 mm, 40 x 40 cm2 profile, and RMSE=0.15%) with the presence of an extra focal halo contribution on the order of 10% of the maximum Gaussian intensity. Using the optimally derived intensity, 10 x 10 and 4 x 4 cm2 profiles were found to be in agreement with measurement with a maximum RMSE=0.49%. The optimized Gaussian and unrestricted values of the electron beam FWHM were both within the range of those inferred by focal spot image measurements performed by Jaffray et al. ["X-ray sources of medical linear accelerators: Focal and extra-focal radiation," Med. Phys. 20, 1417-1427 (1993)]. The inference of an extra focal pretarget electron component may be an indicator of a deficiency in the MC model and needs further investigation.

Simulação computacional de um feixe de fótons de 6 MV em diferentes meios heterogêneos utilizando o código PENELOPE

OBJETIVO: Utilizar o código PENELOPE e desenvolver geometrias onde estão presentes heterogeneidades para simular o comportamento do feixe de fótons nessas condições. MATERIAIS E MÉTODOS: Foram feitas simulações do comportamento da radiação ionizante para o caso homogêneo, apenas água, e para os casos heterogêneos, com diferentes materiais. Consideraram-se geometrias cúbicas para os fantomas e geometrias em forma de paralelepípedos para as heterogeneidades com a seguinte composição: tecido simulador de osso e pulmão, seguindo recomendações da International Commission on Radiological Protection, e titânio, alumínio e prata. Definiram-se, como parâmetros de entrada: a energia e o tipo de partícula da fonte, 6 MV de fótons; a distância fonte-superfície de 100 cm; e o campo de radiação de 10x 10 cm². RESULTADOS: Obtiveram-se curvas de percentual de dose em profundidade para todos os casos. Observou-se que em materiais com densidade eletrônica alta, como a prata, a dose absorvida é maior em relação à dose absorvida no fantoma homogêneo, enquanto no tecido simulador de pulmão a dose é menor. CONCLUSÃO: Os resultados obtidos demonstram a importância de se considerar heterogeneidades nos algoritmos dos sistemas de planejamento usados no cálculo da distribuição de dose nos pacientes, evitando-se sub ou superdosagem dos tecidos próximos às heterogeneidades.

Determining superficial dosimetry for the internal canthus from the Monte Carlo simulation of kV photon and MeV electron beams

This paper presents the findings of an investigation into the Monte Carlo simulation of superficial cancer treatments of an internal canthus site using both kilovoltage photons and megavoltage electrons. The EGSnrc system of codes for the Monte Carlo simulation of the transport of electrons and photons through a phantom representative of either a water phantom or treatment site in a patient is utilised. Two clinical treatment units are simulated: the Varian Medical Systems Clinac 2100C accelerator for 6 MeV electron fields and the Pantak Therapax SXT 150 X-ray unit for 100 kVp photon fields. Depth dose, profile and isodose curves for these simulated units are compared against those measured by ion chamber in a PTW Freiburg MP3 water phantom. Good agreement was achieved away from the surface of the phantom between simulated and measured data. Dose distributions are determined for both kV photon and MeV electron fields in the internal canthus site containing lead and tungsten shielding, rapidly sloping surfaces and different density interfaces. There is a relatively high level of deposition of dose in tissue-bone and tissue-cartilage interfaces in the kV photon fields in contrast to the MeV electron fields. This is reflected in the maximum doses in the PTV of the internal canthus field being 12 Gy for kV photons and 4.8 Gy for MeV electrons. From the dose distributions, DVH and dose comparators are used to assess the simulated treatment fields. Any indication as to which modality is preferable to treat the internal canthus requires careful consideration of many different factors, this investigation provides further perspective in being able to assess which modality is appropriate.

Physical considerations on discrepancies in target volume delineation

BACKGROUND AND PURPOSE

To compare the delineations and interpretations of target volumes by physicians in different radio-oncology centers.

MATERIALS AND METHODS

Eleven Swiss radio-oncology centers delineated volumes according to ICRU 50 recommendations for one prostate and one head and neck case. In order to evaluate the consistency of the volume delineations, the following parameters were determined: 1) the target volumes (GTV, CTV and manually expanded PTV) and their extensions in the three main axes and 2) the correlation of the volume delineated by each pair of centers using the ratio of the intersection to the union (called proximity index).

RESULTS

The delineated prostate volume was 105+/-55 cm(3) for the CTV and 218+/-44 cm(3) for the PTV. The delineated head and neck volume was 46+/-15 cm(3) for the GTV, 327+/-154 cm(3) for the CTV and 528+/-106 cm(3) for the PTV. The mean proximity index for the prostate case was 0.50+/-0.13 for the CTV and 0.57+/-0.11 for the PTV. The proximity index for the head and neck case was 0.45+/-0.09 for the GTV, 0.42+/-0.13 for the CTV and 0.59+/-0.06 for the PTV.

CONCLUSIONS

Large discrepancies between all the delineated target volumes were observed. There was an inverse relationship between the CTV volume and the margin between CTV and PTV, leading to less discrepancies in the PTV than is the CTV delineations. There was more spread in the sagittal and frontal planes due to CT pixel anisotropy, which suggests that radiation oncologists should delineate the target volumes not only in the transverse plane, but also in the sagittal and frontal planes to improve the delineation by allowing a consistency check.

High-performance computing for Monte Carlo radiotherapy calculations

We report on the RTGrid project, which investigates approaches for using high-performance computing infrastructures, such as the grid, in order to reduce the turnaround time of Monte Carlo (MC) simulation-based radiotherapy treatment planning. The main aim of this project is to render accurate dose calculations using MC simulations clinically feasible. To this end, we have successfully implemented and deployed the RTGrid distributed simulation framework for MC dose calculations. In this paper, we present the main experimental findings.

Determination of electron energy, spectral width, and beam divergence at the exit window for clinical megavoltage x-ray beams

Monte Carlo simulations of x-ray beams typically take parameters of the electron beam in the accelerating waveguide to be free parameters. In this paper, a methodology is proposed and implemented to determine the energy, spectral width, and beam divergence of the electron source. All treatment head components were removed from the beam path, leaving only the exit window. With the x-ray target and flattener out of the beam, uncertainties in physical characteristics and relative position of the target and flattening filter, and in spot size, did not contribute to uncertainty in the energy. Beam current was lowered to reduce recombination effects. The measured dose distributions were compared with Monte Carlo simulation of the electron beam through the treatment head to extract the electron source characteristics. For the nominal 6 and 18 MV x-ray beams, the energies were 6.51 +/- 0.15 and 13.9 +/- 0.2 MeV, respectively, with the uncertainties resulting from uncertainties in the detector position in the measurement and in the stopping power in the simulations. Gaussian spectral distributions were used, with full widths at half maximum ranging from 20 +/- 4% at 6 MV to 13 +/- 4% at 18 MV required to match the fall-off portion of the percent-depth ionization curve. Profiles at the depth of maximum dose from simulations that used the manufacturer-specified exit window geometry and no beam divergence were 2-3 cm narrower than measured profiles. Two simulation configurations yielding the measured profile width were the manufacturer-specified exit window thickness with electron source divergences of 3.3 degrees at 6 MV and 1.8 degrees at 18 MV and an exit window 40% thicker than the manufacturer's specification with no beam divergence. With the x-ray target in place (and no flattener), comparison of measured to simulated profiles sets upper limits on the electron source divergences of 0.2 degrees at 6 MV and 0.1 degrees at 18 MV. A method of determining source characteristics without mechanical modification of the treatment head, and therefore feasible in clinics, is presented. The energies and spectral widths determined using this method agree with those determined with only the exit window in the beam path.

Monte Carlo modelling of the treatment line of the Proton Therapy Center in Orsay

This paper presents the main results of a Monte Carlo simulation describing the Orsay Proton Therapy Center (CPO) beam line. The project aimed to obtain a prediction of the dose distribution in a water phantom within 2% accuracy in the dose value and a 2 mm of range. The simulation tool used was MCNPX, version 2.5.0, and included all the elements of the CPO beam line. A new algorithm of multiple Coulomb scattering has been incorporated in MCNPX, resulting in a better prediction of the spatial dose distribution and absolute values of the deposited energy. The simulations of 3D dose profiles in water show a very good agreement with measured data to within 2%. We first performed a comparative analysis of the dosimetry in heterogeneous phantoms between the pencil beam algorithm and MCNPX. The simulations give a better agreement with experimental data compared to the pencil beam approach. In a second phase, we simulated the patient-dependent fields along with the spatial dose distributions in a water phantom. The simulated response of a Pixel chamber located 2 m upstream of the water phantom revealed a good agreement with the measured data to within 1%. The results presented herein support the applicability of Monte Carlo models for absolute dosimetry and for design purposes regarding existing and new beam lines at CPO. This work completes a series of publications reporting the progress in the development of a Monte Carlo simulation tool for the CPO beam line dedicated for the treatment of head and neck tumours.

Low dose megavoltage cone beam computed tomography with an unflattened 4 MV beam from a carbon target

Megavoltage cone beam computed tomography (MVCBCT) is routinely used for visualizing anatomical structures and implanted fiducials for patient positioning in radiotherapy. MVCBCT using a 6 MV treatment beam with high atomic number (Z) target and flattening filter in the beamline, as done conventionally, has lower image quality than can be achieved with a MV beam due to heavy filtration of the low-energy bremsstrahlung. The unflattened beam of a low Z target has an abundance of diagnostic energy photons, detected with modern flat panel detectors with much higher efficiency given the same dose to the patient. This principle guided the development of a new megavoltage imaging beamline (IBL) for a commercial radiotherapy linear accelerator. A carbon target was placed in one of the electron primary scattering foil slots on the target-foil slide. A PROM on a function controller board was programed to put the carbon target in place for MVCBCT. A low accelerating potential of 4.2 MV was used for the IBL to restrict leakage of primary electrons through the target such that dose from x rays dominated the signal in the monitor chamber and the patient surface dose. Results from phantom and cadaver images demonstrated that the IBL had much improved image quality over the treatment beam. For similar imaging dose, the IBL improved the contrast-to-noise ratio by as much as a factor of 3 in soft tissue over that of the treatment beam. The IBL increased the spatial resolution by about a factor of 2, allowing the visualization of finer anatomical details. Images of the cadaver contained useful information with doses as low as 1 cGy. The IBL may be installed on certain models of linear accelerators without mechanical modification and results in significant improvement in the image quality with the same dose, or images of the same quality with less than one-third of the dose.

Monte Carlo modeling of a 6 and 18 MV Varian Clinac medical accelerator for in-field and out-of-field dose calculations: development and validation

There is a serious and growing concern about the increased risk of radiation-induced second cancers and late tissue injuries associated with radiation treatment. To better understand and to more accurately quantify non-target organ doses due to scatter and leakage radiation from medical accelerators, a detailed Monte Carlo model of the medical linear accelerator is needed. This paper describes the development and validation of a detailed accelerator model of the Varian Clinac operating at 6 and 18 MV beam energies. Over 100 accelerator components have been defined and integrated using the Monte Carlo code MCNPX. A series of in-field and out-of-field dose validation studies were performed. In-field dose distributions calculated using the accelerator models were tuned to match measurement data that are considered the de facto 'gold standard' for the Varian Clinac accelerator provided by the manufacturer. Field sizes of 4 cm x 4 cm, 10 cm x 10 cm, 20 cm x 20 cm and 40 cm x 40 cm were considered. The local difference between calculated and measured dose on the percent depth dose curve was less than 2% for all locations. The local difference between calculated and measured dose on the dose profile curve was less than 2% in the plateau region and less than 2 mm in the penumbra region for all locations. Out-of-field dose profiles were calculated and compared to measurement data for both beam energies for field sizes of 4 cm x 4 cm, 10 cm x 10 cm and 20 cm x 20 cm. For all field sizes considered in this study, the average local difference between calculated and measured dose for the 6 and 18 MV beams was 14 and 16%, respectively. In addition, a method for determining neutron contamination in the 18 MV operating model was validated by comparing calculated in-air neutron fluence with reported calculations and measurements. The average difference between calculated and measured neutron fluence was 20%. As one of the most detailed accelerator models for both in-field and out-of-field dose calculations, the model will be combined with anatomically realistic computational patient phantoms into a computational framework to calculate non-target organ doses to patients from various radiation treatment plans.

Benchmarking of Monte Carlo simulation of bremsstrahlung from thick targets at radiotherapy energies

Several Monte Carlo systems were benchmarked against published measurements of bremsstrahlung yield from thick targets for 10-30 MV beams. The quantity measured was photon fluence at 1 m per unit energy per incident electron (spectra), and total photon fluence, integrated over energy, per incident electron (photon yield). Results were reported at 10-30 MV on the beam axis for Al and Pb targets and at 15 MV at angles out to 90 degrees for Be, Al, and Pb targets. Beam energy was revised with improved accuracy of 0.5% using an improved energy calibration of the accelerator. Recently released versions of the Monte Carlo systems EGSNRC, GEANT4, and PENELOPE were benchmarked against the published measurements using the revised beam energies. Monte Carlo simulation was capable of calculation of photon yield in the experimental geometry to 5% out to 30 degrees, 10% at wider angles, and photon spectra to 10% at intermediate photon energies, 15% at lower energies. Accuracy of measured photon yield from 0 to 30 degrees was 5%, 1 s.d., increasing to 7% for the larger angles. EGSNRC and PENELOPE results were within 2 s.d. of the measured photon yield at all beam energies and angles, GEANT4 within 3 s.d. Photon yield at nonzero angles for angles covering conventional field sizes used in radiotherapy (out to 10 degrees), measured with an accuracy of 3%, was calculated within 1 s.d. of measurement for EGSNRC, 2 s.d. for PENELOPE and GEANT4. Calculated spectra closely matched measurement at photon energies over 5 MeV. Photon spectra near 5 MeV were underestimated by as much as 10% by all three codes. The photon spectra below 2-3 MeV for the Be and Al targets and small angles were overestimated by up to 15% when using EGSNRC and PENELOPE, 20% with GEANT4. EGSNRC results with the NIST option for the bremsstrahlung cross section were preferred over the alternative cross section available in EGSNRC and over EGS4. GEANT4 results calculated with the "low energy" physics list were more accurate than those calculated with the "standard" physics list.

An overview of Monte Carlo treatment planning for radiotherapy

The implementation of Monte Carlo dose calculation algorithms in clinical radiotherapy treatment planning systems has been anticipated for many years. Despite a continuous increase of interest in Monte Carlo Treatment Planning (MCTP), its introduction into clinical practice has been delayed by the extent of calculation time required. The development of newer and faster MC codes is behind the commercialisation of the first MC-based treatment planning systems. The intended scope of this article is to provide the reader with a compact 'primer' on different approaches to MCTP with particular attention to the latest developments in the field.

VMC++ versus BEAMnrc: a comparison of simulated linear accelerator heads for photon beams

BEAMnrc, a code for simulating medical linear accelerators based on EGSnrc, has been bench-marked and used extensively in the scientific literature and is therefore often considered to be the gold standard for Monte Carlo simulations for radiotherapy applications. However, its long computation times make it too slow for the clinical routine and often even for research purposes without a large investment in computing resources. VMC++ is a much faster code thanks to the intensive use of variance reduction techniques and a much faster implementation of the condensed history technique for charged particle transport. A research version of this code is also capable of simulating the full head of linear accelerators operated in photon mode (excluding multileaf collimators, hard and dynamic wedges). In this work, a validation of the full head simulation at 6 and 18 MV is performed, simulating with VMC++ and BEAMnrc the addition of one head component at a time and comparing the resulting phase space files. For the comparison, photon and electron fluence, photon energy fluence, mean energy, and photon spectra are considered. The largest absolute differences are found in the energy fluences. For all the simulations of the different head components, a very good agreement (differences in energy fluences between VMC++ and BEAMnrc <1%) is obtained. Only a particular case at 6 MV shows a somewhat larger energy fluence difference of 1.4%. Dosimetrically, these phase space differences imply an agreement between both codes at the <1% level, making VMC++ head module suitable for full head simulations with considerable gain in efficiency and without loss of accuracy.

Anatomical imaging for radiotherapy

The goal of radiation therapy is to achieve maximal therapeutic benefit expressed in terms of a high probability of local control of disease with minimal side effects. Physically this often equates to the delivery of a high dose of radiation to the tumour or target region whilst maintaining an acceptably low dose to other tissues, particularly those adjacent to the target. Techniques such as intensity modulated radiotherapy (IMRT), stereotactic radiosurgery and computer planned brachytherapy provide the means to calculate the radiation dose delivery to achieve the desired dose distribution. Imaging is an essential tool in all state of the art planning and delivery techniques: (i) to enable planning of the desired treatment, (ii) to verify the treatment is delivered as planned and (iii) to follow-up treatment outcome to monitor that the treatment has had the desired effect. Clinical imaging techniques can be loosely classified into anatomic methods which measure the basic physical characteristics of tissue such as their density and biological imaging techniques which measure functional characteristics such as metabolism. In this review we consider anatomical imaging techniques. Biological imaging is considered in another article. Anatomical imaging is generally used for goals (i) and (ii) above. Computed tomography (CT) has been the mainstay of anatomical treatment planning for many years, enabling some delineation of soft tissue as well as radiation attenuation estimation for dose prediction. Magnetic resonance imaging is fast becoming widespread alongside CT, enabling superior soft-tissue visualization. Traditionally scanning for treatment planning has relied on the use of a single snapshot scan. Recent years have seen the development of techniques such as 4D CT and adaptive radiotherapy (ART). In 4D CT raw data are encoded with phase information and reconstructed to yield a set of scans detailing motion through the breathing, or cardiac, cycle. In ART a set of scans is taken on different days. Both allow planning to account for variability intrinsic to the patient. Treatment verification has been carried out using a variety of technologies including: MV portal imaging, kV portal/fluoroscopy, MVCT, conebeam kVCT, ultrasound and optical surface imaging. The various methods have their pros and cons. The four x-ray methods involve an extra radiation dose to normal tissue. The portal methods may not generally be used to visualize soft tissue, consequently they are often used in conjunction with implanted fiducial markers. The two CT-based methods allow measurement of inter-fraction variation only. Ultrasound allows soft-tissue measurement with zero dose but requires skilled interpretation, and there is evidence of systematic differences between ultrasound and other data sources, perhaps due to the effects of the probe pressure. Optical imaging also involves zero dose but requires good correlation between the target and the external measurement and thus is often used in conjunction with an x-ray method. The use of anatomical imaging in radiotherapy allows treatment uncertainties to be determined. These include errors between the mean position at treatment and that at planning (the systematic error) and the day-to-day variation in treatment set-up (the random error). Positional variations may also be categorized in terms of inter- and intra-fraction errors. Various empirical treatment margin formulae and intervention approaches exist to determine the optimum strategies for treatment in the presence of these known errors. Other methods exist to try to minimize error margins drastically including the currently available breath-hold techniques and the tracking methods which are largely in development. This paper will review anatomical imaging techniques in radiotherapy and how they are used to boost the therapeutic benefit of the treatment.

Monte Carlo simulation of large electron fields

Two Monte Carlo systems, EGSnrc and Geant4, the latter with two different 'physics lists,' were used to calculate dose distributions in large electron fields used in radiotherapy. Source and geometry parameters were adjusted to match calculated results to measurement. Both codes were capable of accurately reproducing the measured dose distributions of the six electron beams available on the accelerator. Depth penetration matched the average measured with a diode and parallel-plate chamber to 0.04 cm or better. Calculated depth dose curves agreed to 2% with diode measurements in the build-up region, although for the lower beam energies there was a discrepancy of up to 5% in this region when calculated results are compared to parallel-plate measurements. Dose profiles at the depth of maximum dose matched to 2-3% in the central 25 cm of the field, corresponding to the field size of the largest applicator. A 4% match was obtained outside the central region. The discrepancy observed in the bremsstrahlung tail in published results that used EGS4 is no longer evident. Simulations with the different codes and physics lists used different source energies, incident beam angles, thicknesses of the primary foils, and distance between the primary and secondary foil. The true source and geometry parameters were not known with sufficient accuracy to determine which parameter set, including the energy of the source, was closest to the truth. These results underscore the requirement for experimental benchmarks of depth penetration and electron scatter for beam energies and foils relevant to radiotherapy.

A "HOWFARLESS" option to increase efficiency of homogeneous phantom calculations with DOSXYZnrc

This paper describes a "HOWFARLESS" transport option, which has been added to DOSXYZnrc to increase the efficiency of beam commissioning calculations in homogeneous phantoms. The algorithm speeds up charged particle transport by only considering the distance to the extreme outer boundaries of the phantom, thus eliminating the need to stop at voxel boundaries. Dose is deposited by approximating the total curved charged particle steps by two straight-line steps joined at a hinge point. Good agreement with normal simulations is achieved at all beam energies and for all practical maximum step lengths with a 1:1 mixture of approximations based on the initial position/ direction of the particle and on its final position/direction. Use of the "HOWFARLESS" option in phantom calculations for 6 and 18 MV photon beams (10 x 10 cm2 and 40 x 40 cm2 fields) from BEAMnrc-simulated accelerators increases the efficiency at the optimum photon splitting number by a factor of 2.9-5.4 when the exact EGSnrc boundary crossing algorithm (BCA) is used and by 51%-89% when the faster PRESTA-I BCA is employed. The efficiency gain due to the "HOWFARLESS" transport option increases with increasing beam energy and decreases with increasing field/dose voxel size. Efficiency improvement is greater when the efficiency of the particle source itself is not a factor, and in such cases the "HOWFARLESS" option improves the DOSXYZnrc efficiency by up to a factor of 13.1 (exact BCA) or 3.5 (PRESTA-I BCA) for photon beams, and up to a factor of 17.2 (exact BCA) or 5.2 (PRESTA-I BCA) for electron beams.

Dosimetric characteristics of unflattened 6MV photon beams of a clinical linear accelerator: A Monte Carlo study

The purpose of this study is to evaluate the dosimetric properties of a flattening filter free 6 MV photon beam. The 6 MV photon beam of a Varian Clinac 21 EX linac was modeled using the MCNP4C Monte Carlo (MC) code. Dosimetric features including central axis absorbed doses, beam profiles and photon energy spectra were calculated for flattened and unflattened 6 MV photon beams. A substantial increase in the dose rate was seen for the unflattened beam, which was decreased with field size and depth. The penumbra width was decreased less than 0.2 mm (about 5%) and a 25% decrease in out-of-field dose was observed for the unflattened beam. The photon energy spectra were softer for the unflattened beam and the mean energies of spectra were higher for smaller field size. Our study showed that increase in the dose rate and lower out-of-field dose could be considered as practical advantages for unflattened 6 MV beams.

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.

An investigation of dose changes for therapeutic kilovoltage x-ray beams with underlying lead shielding

Kilovoltage x-ray beams are used to treat cancer on or close to the skin surface. Many clinical cases use high atomic number materials as shielding to reduce dose to underlying healthy tissues. In this work, we have investigated the effect on both the surface dose and depth doses in a water phantom with lead shielding at depth in the phantom. The EGSnrc Monte Carlo code was used to simulate the water phantom and to calculate the surface doses and depth doses using primary x-ray beam spectra derived from an analytical model. The x-ray beams were in the energy range of 75-135 kVp with field sizes of 2, 5 and 8 cm diameter. The lead sheet was located beneath the water surface at depths ranging from 0.5-7.5 cm. The surface dose decreased as the lead was positioned closer to the water surface and as the field size was increased. The variation in surface dose as a function of x-ray beam energy was only small but the maximum reduction occurred for the 100 kVp x-ray beam. For the 8 cm diameter field with the lead at 1 cm depth and using the 100 kVp x-ray beam, the surface dose was reduced to 0.898 of the surface dose in the water phantom only. Measured surface dose changes, using a Farmer-type ionization chamber, agreed with the Monte Carlo calculated doses. Calculated depth doses in water with a lead sheet positioned below the surface showed that the dose fall-off increased as the lead was positioned closer to the water surface as compared to the depth dose in the water phantom only. Monte Carlo calculations of the total x-ray beam spectrum at the water surface showed that the total fluence decreased due to a reduction in backscatter from within the water and very little backscatter from the lead. The mean energy of the x-ray spectrum varied less than 1 keV, with the lead at 1 cm beneath the water phantom surface. As the Monte Carlo calculations showed good agreement with the measured results, this method can be used to verify surface dose changes in clinical situations where measurements are difficult. The clinical impact of the use of lead must be considered in the dose prescription for patients being treated with kilovoltage x-ray beams.

Tenth value layers for 60Co gamma rays and for 4, 6, 10, 15, and 18 MV x rays in concrete for beams of cone angles between 0 degrees and 14 degrees calculated by Monte Carlo simulation

The calculation of shielding barrier thicknesses for radiation therapy facilities according to the NCRP formalism is based on the use of broad beams (that is, the maximum possible field sizes). However, in practice, treatment fields used in radiation therapy are, on average, less than half the maximum size. Indeed, many contemporary treatment techniques call for reduced field sizes to reduce co-morbidity and the risk of second cancers. Therefore, published tenth value layers (TVLs) for shielding materials do not apply to these very small fields. There is, hence, a need to determine the TVLs for various beam modalities as a function of field size. The attenuation of (60)Co gamma rays and photons of 4, 6, 10, 15, and 18 MV bremsstrahlung x ray beams by concrete has been studied using the Monte Carlo technique (MCNP version 4C2) for beams of half-opening angles of 0 degrees , 3 degrees , 6 degrees , 9 degrees , 12 degrees , and 14 degrees . The distance between the x-ray source and the distal surface of the shielding wall was fixed at 600 cm, a distance that is typical for modern radiation therapy rooms. The maximum concrete thickness varied between 76.5 cm and 151.5 cm for (60)Co and 18 MV x rays, respectively. Detectors were placed at 630 cm, 700 cm, and 800 cm from the source. TVLs have been determined down to the third TVL. Energy spectra for 4, 6, 10, 15, and 18 MV x rays for 10 x 10 cm(2) and 40 x 40 cm(2) field sizes were used to generate depth dose curves in water that were compared with experimentally measured values.

Patient-specific dosimetry of conventional and intensity modulated radiation therapy using a novel full Monte Carlo phase space reconstruction method from electronic portal images

Electronic portal imagers have promising dosimetric applications in external beam radiation therapy. In this study a patient dose computation algorithm based on Monte Carlo (MC) simulations and on portal images is developed and validated. The patient exit fluence from primary photons is obtained from the portal image after correction for scattered radiation. The scattered radiation at the portal imager and the spectral energy distribution of the primary photons are estimated from MC simulations at the treatment planning stage. The patient exit fluence and the spectral energy distribution of the primary photons are then used to ray-trace the photons from the portal image towards the source through the CT geometry of the patient. Photon weights which reflect the probability of a photon being transmitted are computed during this step. A dedicated MC code is used to transport back these photons from the source through the patient CT geometry to obtain patient dose. Only Compton interactions are considered. This code also produces a reconstructed portal image which is used as a verification tool to ensure that the dose reconstruction is reliable. The dose reconstruction algorithm is compared against MC dose calculation (MCDC) predictions and against measurements in phantom. The reconstructed absolute absorbed doses and the MCDC predictions in homogeneous and heterogeneous phantoms agree within 3% for simple open fields. Comparison with film-measured relative dose distributions for IMRT fields yields agreement within 3 mm, 5%. This novel dose reconstruction algorithm allows for daily patient-specific dosimetry and verification of patient movement.