In this report a new EGS4 version, called EGSnrc to reflect the substantial changes made to the original code is reported, which incorporates a new any-angle multiple elastic scattering theory, an improved electron-step algorithm, a correct implementation of the fictitious cross section method for sampling distances between discrete interactions, a more accurate evaluation of energy loss, as well as an exact boundary crossing algorithm. It is demonstrated that EGSnrc allows for an artifact free Monte Carlo simulation of ion chamber response and backscattering, situations that have been considered in the past as the two of the most stringent tests of condensed history Monte Carlo codes. A detailed discussion of the effect of the various components of the condensed history simulation of electron transport on the simulated ion chamber response is given in the accompanying paper. ©

The Monte Carlo (MC) method has been shown through many research studies to calculate accurate dose distributions for clinical radiotherapy, particularly in heterogeneous patient tissues where the effects of electron transport cannot be accurately handled with conventional, deterministic dose algorithms. Despite its proven accuracy and the potential for improved dose distributions to influence treatment outcomes, the long calculation times previously associated with MC simulation rendered this method impractical for routine clinical treatment planning. However, the development of faster codes optimized for radiotherapy calculations and improvements in computer processor technology have substantially reduced calculation times to, in some instances, within minutes on a single processor. These advances have motivated several major treatment planning system vendors to embark upon the path of MC techniques. Several commercial vendors have already released or are currently in the process of releasing MC algorithms for photon and/or electron beam treatment planning. Consequently, the accessibility and use of MC treatment planning algorithms may well become widespread in the radiotherapy community. With MC simulation, dose is computed stochastically using first principles; this method is therefore quite different from conventional dose algorithms. Issues such as statistical uncertainties, the use of variance reduction techniques, the ability to account for geometric details in the accelerator treatment head simulation, and other features, are all unique components of a MC treatment planning algorithm. Successful implementation by the clinical physicist of such a system will require an understanding of the basic principles of MC techniques. The purpose of this report, while providing education and review on the use of MC simulation in radiotherapy planning, is to set out, for both users and developers, the salient issues associated with clinical implementation and experimental verification of MC dose algorithms. As the MC method is an emerging technology, this report is not meant to be prescriptive. Rather, it is intended as a preliminary report to review the tenets of the MC method and to provide the framework upon which to build a comprehensive program for commissioning and routine quality assurance of MC-based treatment planning systems.

No abstract

This article presents the implementation of several variance reduction techniques that dramatically improve the simulation efficiency of ion chamber dose and perturbation factor calculations. The cavity user code for the EGSnrc Monte Carlo code system is extended by photon cross-section enhancement (XCSE), an intermediate phase-space storage (IPSS) technique, and a correlated sampling (CS) scheme. XCSE increases the density of photon interaction sites inside and in the vicinity of the chamber and results-in combination with a Russian Roulette game for electrons that cannot reach the cavity volume-in an increased efficiency of up to a factor of 350 for calculating dose in a Farmer type chamber placed at 10 cm depth in a water phantom. In combination with the IPSS and CS techniques, the efficiency for the calculation of the central electrode perturbation factor Pcel can be increased by up to three orders of magnitude for a single chamber location and by nearly four orders of magnitude when considering the Pcel variation with depth or with distance from the central axis in a large field photon beam. The intermediate storage of the phase-space properties of particles entering a volume that contains many possible chamber locations leads to efficiency improvements by a factor larger than 500 when computing a profile of chamber doses in the field of a linear accelerator photon beam. All techniques are combined in a new EGSnrc user code egs_chamber. Optimum settings for the variance reduction parameters are investigated and are reported for a Farmer type ion chamber. A few example calculations illustrating the capabilities of the egs_chamber code are presented.

A new model for calculating electron beam dose has been developed. The algorithm is based on a two- or three-dimensional geometry defined by computerized tomography (CT) images. The Monte Carlo technique was used to solve the electron transport equation. However, in contrast to conventional Monte Carlo models (EGS4) several approximations and simplifications in the description of elementary electron processes were introduced reducing in this manner the computational time by a factor of about 35 without significant loss in accuracy. The Monte Carlo computer program does not need any precalculated data. The random access memory required is about 16 Mbytes for a 128(2) X 50 matrix, depending on the resolution of the CT cube. The Voxel Monte Carlo model (VMC) was tested in comparison to calculations by EGS4 and the "Hogstrom algorithm" (MDAH) using several fictive phantoms. In all cases a good coincidence has been found between EGS4 and VMC, especially near tissue inhomogeneities, whereas the MDAH algorithm has produced dose underestimations of up to 40%.