Anders Ahnesjö

Small fields: nonequilibrium radiation dosimetry.

Advances in radiation treatment with beamlet-based intensity modulation, image-guided radiation therapy, and stereotactic radiosurgery (including specialized equipments like CyberKnife, Gamma Knife, tomotherapy, and high-resolution multileaf collimating systems) have resulted in the use of reduced treatment fields to a subcentimeter scale. Compared to the traditional radiotherapy with fields > or =4 x 4 cm2, this can result in significant uncertainty in the accuracy of clinical dosimetry. The dosimetry of small fields is challenging due to nonequilibrium conditions created as a consequence of the secondary electron track lengths and the source size projected through the collimating system that are comparable to the treatment field size. It is further complicated by the prolonged electron tracks in the presence of low-density inhomogeneities. Also, radiation detectors introduced into such fields usually perturb the level of disequilibrium. Hence, the dosimetric accuracy previously achieved for standard radiotherapy applications is at risk for both absolute and relative dose determination. This article summarizes the present knowledge and gives an insight into the future procedures to handle the nonequilibrium radiation dosimetry problems. It is anticipated that new miniature detectors with controlled perturbations and corrections will be available to meet the demand for accurate measurements. It is also expected that the Monte Carlo techniques will increasingly be used in assessing the accuracy, verification, and calculation of dose, and will aid perturbation calculations of detectors used in small and highly conformal radiation beams. rican Association of Physicists in Medicine.

Dose calculations for external photon beams in radiotherapy

Dose calculation methods for photon beams are reviewed in the context of radiation therapy treatment planning. Following introductory summaries on photon beam characteristics and clinical requirements on dose calculations, calculation methods are described in order of increasing explicitness of particle transport. The simplest are dose ratio factorizations limited to point dose estimates useful for checking other more general, but also more complex, approaches. Some methods incorporate detailed modelling of scatter dose through differentiation of measured data combined with various integration techniques. State-of-the-art methods based on point or pencil kernels, which are derived through Monte Carlo simulations, to characterize secondary particle transport are presented in some detail. Explicit particle transport methods, such as Monte Carlo, are briefly summarized. The extensive literature on beam characterization and handling of treatment head scatter is reviewed in the context of providing phase space data for kernel based and/or direct Monte Carlo dose calculations. Finally, a brief overview of inverse methods for optimization and dose reconstruction is provided.

Accelerator beam data commissioning equipment and procedures: Report of the TG-106 of the Therapy Physics Committee of the AAPM

For commissioning a linear accelerator for clinical use, medical physicists are faced with many challenges including the need for precision, a variety of testing methods, data validation, the lack of standards, and time constraints. Since commissioning beam data are treated as a reference and ultimately used by treatment planning systems, it is vitally important that the collected data are of the highest quality to avoid dosimetric and patient treatment errors that may subsequently lead to a poor radiation outcome. Beam data commissioning should be performed with appropriate knowledge and proper tools and should be independent of the person collecting the data. To achieve this goal, Task Group 106 (TG-106) of the Therapy Physics Committee of the American Association of Physicists in Medicine was formed to review the practical aspects as well as the physics of linear accelerator commissioning. The report provides guidelines and recommendations on the proper selection of phantoms and detectors, setting up of a phantom for data acquisition (both scanning and no-scanning data), procedures for acquiring specific photon and electron beam parameters and methods to reduce measurement errors (<1%), beam data processing and detector size convolution for accurate profiles. The TG-106 also provides a brief discussion on the emerging trend in Monte Carlo simulation techniques in photon and electron beam commissioning. The procedures described in this report should assist a qualified medical physicist in either measuring a complete set of beam data, or in verifying a subset of data before initial use or for periodic quality assurance measurements. By combining practical experience with theoretical discussion, this document sets a new standard for beam data commissioning.

Implementation of pencil kernel and depth penetration algorithms for treatment planning of proton beams

The implementation of two algorithms for calculating dose distributions for radiation therapy treatment planning of intermediate energy proton beams is described. A pencil kernel algorithm and a depth penetration algorithm have been incorporated into a commercial three dimensional treatment planning system (Helax-TMS, Helax AB, Sweden) to allow conformal planning techniques using irregularly shaped fields, proton range modulation, range modification and dose calculation for non-coplanar beams. The pencil kernel algorithm is developed from the Fermi Eyges formalism and Molière multiple-scattering theory with range straggling corrections applied. The depth penetration algorithm is based on the energy loss in the continuous slowing down approximation with simple correction factors applied to the beam penumbra region and has been implemented for fast, interactive treatment planning. Modelling of the effects of air gaps and range modifying device thickness and position are implicit to both algorithms. Measured and calculated dose values are compared for a therapeutic proton beam in both homogeneous and heterogeneous phantoms of varying complexity. Both algorithms model the beam penumbra as a function of depth in a homogeneous phantom with acceptable accuracy. Results show that the pencil kernel algorithm is required for modelling the dose perturbation effects from scattering in heterogeneous media.

Report of AAPM Therapy Physics Committee Task Group 74: in-air output ratio, Sc, for megavoltage photon beams.

The concept of in-air output ratio (Sc) was introduced to characterize how the incident photon fluence per monitor unit (or unit time for a Co-60 unit) varies with collimator settings. However, there has been much confusion regarding the measurement technique to be used that has prevented the accurate and consistent determination of Sc. The main thrust of the report is to devise a theoretical and measurement formalism that ensures interinstitutional consistency of Sc. The in-air output ratio, Sc, is defined as the ratio of primary collision water kerma in free-space, Kp, per monitor unit between an arbitrary collimator setting and the reference collimator setting at the same location. Miniphantoms with sufficient lateral and longitudinal thicknesses to eliminate electron contamination and maintain transient electron equilibrium are recommended for the measurement of Sc. The authors present a correction formalism to extrapolate the correct Sc from the measured values using high-Z miniphantom. Miniphantoms made of high-Z material are used to measure Sc for small fields (e.g., IMRT or stereotactic radiosurgery). This report presents a review of the components of Sc, including headscatter, source-obscuring, and monitor-backscattering effects. A review of calculation methods (Monte Carlo and empirical) used to calculate Sc for arbitrary shaped fields is presented. The authors discussed the use of Sc in photon dose calculation algorithms, in particular, monitor unit calculation. Finally, a summary of Sc data (from RPC and other institutions) is included for QA purposes.

Dose calculation in brachytherapy for a source using a primary and scatter dose separation technique

A dose calculation algorithm for brachytherapy is presented that reduces errors in absolute dose calculation and facilitates new techniques for modelling heterogeneity effects from tissues, internal shields and superficially positioned sources. The algorithm is based on Monte Carlo simulations for specific source and applicator combinations. The dose is scored separately, in absolute units, for the primary and different categories of scatter according to the photon scatter generation. Radial dose distributions for the primary dose and the total scatter dose are parametrized using functions based on simple one-dimensional transport theory. The fitted radial parameters are functions of the angle to the long axis of the source to account for the anisotropy of the dose distribution. The kerma in air at the reference point 1 m from the source is also simulated using Monte Carlo techniques and both the dose and kerma are normalized per source emitted radiant energy. The calculated kerma per radiant energy is used together with the measured reference air kerma rate and the ratio of the dose to the kerma to calibrate the calculated absolute dose rate. Data are presented for an 192Ir cylindrical source, in combination with water, nylon and stainless steel applicators. Values of the radial dose profiles, specific dose rate constants and corrections to the air kerma for attenuation and scatter in air are calculated. Anisotropy functions for the 192Ir source and a water-equivalent applicator are compared to published values. The effects of the applicator wall material on the radial dose distribution are also discussed.

Modeling transmission and scatter for photon beam attenuators.

The development of treatment planning methods in radiation therapy requires dose calculation methods that are both accurate and general enough to provide a dose per unit monitor setting for a broad variety of fields and beam modifiers. The purpose of this work was to develop models for calculation of scatter and transmission for photon beam attenuators such as compensating filters, wedges, and block trays. The attenuation of the beam is calculated using a spectrum of the beam, and a correction factor based on attenuation measurements. Small angle coherent scatter and electron binding effects on scattering cross sections are considered by use of a correction factor. Quality changes in beam penetrability and energy fluence to dose conversion are modeled by use of the calculated primary beam spectrum after passage through the attenuator. The beam spectra are derived by the depth dose effective method, i.e., by minimizing the difference between measured and calculated depth dose distributions, where the calculated distributions are derived by superposing data from a database for monoenergetic photons. The attenuator scatter is integrated over the area viewed from the calculation point of view using first scatter theory. Calculations are simplified by replacing the energy and angular-dependent cross-section formulas with the forward scatter constant r2(0) and a set of parametrized correction functions. The set of corrections include functions for the Compton energy loss, scatter attenuation, and secondary bremsstrahlung production. The effect of charged particle contamination is bypassed by avoiding use of dmax for absolute dose calibrations. The results of the model are compared with scatter measurements in air for copper and lead filters and with dose to a water phantom for lead filters for 4 and 18 MV. For attenuated beams, downstream of the buildup region, the calculated results agree with measurements on the 1.5% level. The accuracy was slightly less in situations where the scatter component is very large, as for very large fields with very short filter to detector distances. The implementation of the model into treatment planning systems is discussed.

The collapsed cone superposition algorithm applied to scatter dose calculations in brachytherapy

Methods for scatter dose calculations in brachytherapy have been developed based on the collapsed cone superposition algorithm. The methods account for effects on the scatter dose caused by the three-dimensional distribution of heterogeneities in the irradiated volume and are considerably faster than methods based on straightforward superposition of kernels or direct Monte Carlo simulations. Use of a successive-scattering approach, in which the dose contribution from once- and multiply scattered photons are calculated separately, was found superior to conventional superposition using a single point kernel for all scatter generations. Use of the successive-scattering approach significantly reduces artifacts stemming from steep fluence gradients, typical of the brachytherapy geometry and critical for the collapsed cone approximation. The algorithm is tested versus Monte Carlo simulations for point sources of energies 28.4, 100, 350, and 662 keV. Results agree well for both a homogeneous water phantom and an air-water half-phantom.

The IMRT information process—mastering the degrees of freedom in external beam therapy

The techniques and procedures for intensity-modulated radiation therapy (IMRT) are reviewed in the context of the information process central to treatment planning and delivery of IMRT. A presentation is given of the evolution of the information based radiotherapy workflow and dose delivery techniques, as well as the volume and planning concepts for relating the dose information to image based patient representations. The formulation of the dose shaping process as an optimization problem is described. The different steps in the calculation flow for determination of machine parameters for dose delivery are described starting from the formulation of optimization objectives over dose calculation to optimization procedures. Finally, the main elements of the quality assurance procedure necessary for implementing IMRT clinically are reviewed.

Accounting for high Z shields in brachytherapy using collapsed cone superposition for scatter dose calculation.

Common clinical brachytherapy treatment planning algorithms perform at best one-dimensional corrections for high Z heterogeneities that will be inaccurate for intermediate energies (60-100 keV). The development of fast methods for a three-dimensional dose calculation to account for high Z materials in this energy range is important, e.g., to fully utilize the potential of patient individualized shields using isotopes such as 241Am and 169Yb. In this work we use the collapsed cone superposition algorithm to calculate the scatter dose contribution around partly lead-shielded point sources at 60, 100, and 350 keV. Methods to scale point kernels for water into kernels for high Z materials are derived. The scaling accounts for scattered photon spectral differences between materials and thus goes beyond the common density scaling approach. Compared to Monte Carlo simulations, the results of our algorithm yield agreements on the unshielded side to within 3% at 350 and 60 keV and to within 7% at 100 keV out to distances of 8 cm from the source. The effect of the shield in the center of the unshielded region is small at 350 keV but significant and occurs at short distances at 100 and 60 keV. At 60 keV, the shield causes a dose reduction of around 10%, 1 cm from the source on the unshielded side. At 100 keV, the reverse effect is seen, the insertion of shields leading to the total dose being increased by about 10% at 1 cm. That one-dimensional algorithms are incapable of predicting these changes shows the importance of accounting for the full three-dimensional geometry in correctly determining the scatter dose contribution.