Electron beam quality correction factors for plane-parallel ionization chambers: Monte Carlo calculations using the PENELOPE system

On the wall perturbation correction for a parallel-plate NACP-02 chamber in clinical electron beams

PURPOSE

In recent years, several Monte Carlo studies have been published concerning the perturbation corrections of a parallel-plate chamber in clinical electron beams. In these studies, a strong depth dependence of the relevant correction factors (p(wall) and P(cav)) for depth beyond the reference depth is recognized and it has been shown that the variation with depth is sensitive to the choice of the chamber's effective point of measurement. Recommendations concerning the positioning of parallel-plate ionization chambers in clinical electron beams are not the same for all current dosimetry protocols. The IAEA TRS-398 as well as the IPEM protocol and the German protocol DIN 6800-2 interpret the depth of measurement within the phantom as the water equivalent depth, i.e., the nonwater equivalence of the entrance window has to be accounted for by shifting the chamber by an amount deltaz. This positioning should ensure that the primary electrons traveling from the surface of the water phantom through the entrance window to the chamber's reference point sustain the same energy loss as the primary electrons in the undisturbed phantom. The objective of the present study is the determination of the shift deltaz for a NACP-02 chamber and the calculation of the resulting wall perturbation correction as a function of depth. Moreover, the contributions of the different chamber walls to the wall perturbation correction are identified.

METHODS

The dose and fluence within the NACP-02 chamber and a wall-less air cavity is calculated using the Monte Carlo code EGSnrc in a water phantom at different depths for different clinical electron beams. In order to determine the necessary shift to account for the nonwater equivalence of the entrance window, the chamber is shifted in steps deltaz around the depth of measurement. The optimal shift deltaz is determined from a comparison of the spectral fluence within the chamber and the bare cavity. The wall perturbation correction is calculated as the ratio between doses for the complete chamber and a wall-less air cavity.

RESULTS

The high energy part of the fluence spectra within the chamber strongly varies even with small chamber shifts, allowing the determination of deltaz within micrometers. For the NACP-02 chamber a shift deltaz = -0.058 cm results. This value is independent of the energy of the primary electrons as well as of the depth within the phantom and it is in good agreement with the value recommended in the German dosimetry protocol. Applying this shift, the calculated wall perturbation correction as a function of depth is varying less than 1% from zero up to the half value depth R50 for electron energies in the range of 6-21 MeV. The remaining depth dependence can mainly be attributed to the scatter properties of the entrance window. When neglecting the nonwater equivalence of the entrance window, the variation of p(wall) with depth is up to 10% and more, especially for low electron energies.

CONCLUSIONS

The variation of the wall perturbation correction for the NACP-02 chamber in clinical electron beams strongly depends on the positioning of the chamber. Applying a shift deltaz = -0.058 cm toward the focus ensures that the primary electron spectrum within the chamber bears the largest resemblance to the fluence of a wall-less cavity. Hence, the influence of the chamber walls on the perturbation correction can be separated out and the residual variation of p(wall) with depth is minimized.

Extraction of depth-dependent perturbation factors for parallel-plate chambers in electron beams using a plastic scintillation detector

PURPOSE

This work presents the experimental extraction of the overall perturbation factor PQ in megavoltage electron beams for NACP-02 and Roos parallel-plate ionization chambers using a plastic scintillation detector (PSD).

METHODS

The authors used a single scanning PSD mounted on a high-precision scanning tank to measure depth-dose curves in 6, 12, and 18 MeV clinical electron beams. The authors also measured depth-dose curves using the NACP-02 and PTW Roos chambers.

RESULTS

The authors found that the perturbation factors for the NACP-02 and Roos chambers increased substantially with depth, especially for low-energy electron beams. The experimental results were in good agreement with the results of Monte Carlo simulations reported by other investigators. The authors also found that using an effective point of measurement (EPOM) placed inside the air cavity reduced the variation of perturbation factors with depth and that the optimal EPOM appears to be energy dependent.

CONCLUSIONS

A PSD can be used to experimentally extract perturbation factors for ionization chambers. The dosimetry protocol recommendations indicating that the point of measurement be placed on the inside face of the front window appear to be incorrect for parallel-plate chambers and result in errors in the R50 of approximately 0.4 mm at 6 MeV, 1.0 mm at 12 MeV, and 1.2 mm at 18 MeV.

On the p(dis) correction factor for cylindrical chambers

The authors of a recent paper (Wang and Rogers 2009 Phys. Med. Biol. 54 1609) have used the Monte Carlo method to simulate the 'classical' experiment made more than 30 years ago by Johansson et al (1978 National and International Standardization of Radiation Dosimetry (Atlanta 1977) vol 2 (Vienna: IAEA) pp 243-70) on the displacement (or replacement) perturbation correction factor p(dis) for cylindrical chambers in 60Co and high-energy photon beams. They conclude that an 'unreasonable normalization at dmax' of the ionization chambers response led to incorrect results, and for the IAEA TRS-398 Code of Practice, which uses ratios of those results, 'the difference in the correction factors can lead to a beam calibration deviation of more than 0.5% for Farmer-like chambers'. The present work critically examines and questions some of the claims and generalized conclusions of the paper. It is demonstrated that for real, commercial Farmer-like chambers, the possible deviations in absorbed dose would be much smaller (typically 0.13%) than those stated by Wang and Rogers, making the impact of their proposed values negligible on practical high-energy photon dosimetry. Differences of the order of 0.4% would only appear at the upper extreme of the energies potentially available for clinical use (around 25 MV) and, because lower energies are more frequently used, the number of radiotherapy photon beams for which the deviations would be larger than say 0.2% is extremely small. This work also raises concerns on the proposed value of pdis for Farmer chambers at the reference quality of 60Co in relation to their impact on electron beam dosimetry, both for direct dose determination using these chambers and for the cross-calibration of plane-parallel chambers. The proposed increase of about 1% in p(dis) (compared with TRS-398) would lower the kQ factors and therefore Dw in electron beams by the same amount. This would yield a severe discrepancy with the current good agreement between electron dosimetry based on an electron cross-calibrated plane-parallel chamber (against a Farmer) or on a directly 60Co calibrated plane-parallel chamber, which is not likely to be in error by 1%. It is suggested that the influence of the 60Co source spectrum used in the simulations may not be negligible for calculations aimed at an uncertainty level of 0.1%.

penMesh--Monte Carlo radiation transport simulation in a triangle mesh geometry

We have developed a general-purpose Monte Carlo simulation code, called penMesh, that combines the accuracy of the radiation transport physics subroutines from PENELOPE and the flexibility of a geometry based on triangle meshes. While the geometric models implemented in most general-purpose codes--such as PENELOPE's quadric geometry--impose some limitations in the shape of the objects that can be simulated, triangle meshes can be used to describe any free-form (arbitrary) object. Triangle meshes are extensively used in computer-aided design and computer graphics. We took advantage of the sophisticated tools already developed in these fields, such as an octree structure and an efficient ray-triangle intersection algorithm, to significantly accelerate the triangle mesh ray-tracing. A detailed description of the new simulation code and its ray-tracing algorithm is provided in this paper. Furthermore, we show how it can be readily used in medical imaging applications thanks to the detailed anatomical phantoms already available. In particular, we present a whole body radiography simulation using a triangulated version of the anthropomorphic NCAT phantom. An example simulation of scatter fraction measurements using a standardized abdomen and lumbar spine phantom, and a benchmark of the triangle mesh and quadric geometries in the ray-tracing of a mathematical breast model, are also presented to show some of the capabilities of penMesh.

Effect of oblique X-ray incidence in flat-panel computed tomography of the breast

We quantify the variation in resolution due to anisotropy caused by oblique X-ray incidence in indirect flat-panel detectors for computed tomography breast imaging systems. We consider a geometry and detector type utilized in breast computed tomography (CT) systems currently being developed. Our methods rely on mantis, a combined X-ray, electron, and optical Monte Carlo transport open source code. The physics models are the most accurate available in general-purpose Monte Carlo packages in the diagnostic energy range. We consider maximum-obliquity angles of 10 ( degrees ) and 13 ( degrees ) at the centers of the 30 and 40 cm detector edges, respectively, and 16 ( degrees ) at the corner of the detector. Our results indicate that blur is asymmetric and that the resolution properties vary significantly with the angle (or location) of incidence. Our results suggest that the asymmetry can be as high as a factor of 2.6 between orthogonal directions. Anisotropy maps predicted by mantis provide an understanding of the effect that such variations have on the imaging system and allow more accurate modeling and optimization of breast CT systems. These maps of anisotropy across the detector could lead to improved reconstruction and help motivate physics-based strategies for computer detection of breast lesions.

Anisotropic imaging performance in breast tomosynthesis

We describe the anisotropy in imaging performance caused by oblique x-ray incidence in indirect detectors for breast tomosynthesis based on columnar scintillator screens. We use MANTIS, a freely available combined x-ray, electron, and optical Monte Carlo transport package which models the indirect detection processes in columnar screens, interaction by interaction. The code has been previously validated against published optical distributions. In this article, initial validation results are provided concerning the blur for particular designs of phosphor screens for which some details with respect to the columnar geometry are available from scanning electron microscopy. The polyenergetic x-ray spectrum utilized comes from a database of experimental data for three different anode/filter/kVp combinations: Mo/Mo at 28 kVp, Rh/Rh at 28 kVp, and W/Al at 42 kVp. The x-ray spectra were then filtered with breast tissue (3, 4, and 6 cm thickness), compression paddle, and support base, according to the oblique paths determined by the incidence angle. The composition of the breast tissue was 50%/50% adipose/glandular tissue mass ratio. Results are reported on the pulse-height statistics of the light output and on spatial blur, expressed as the response of the detector to a pencil beam with a certain incidence angle. Results suggest that the response is nonsymmetrical and that the resolution properties of a tomosynthesis system vary significantly with the angle of x-ray incidence. In contrast, it is found that the noise due to the variability in the number of light photons detected per primary x-ray interaction changes only a few percent. The anisotropy in the response is not less in screens with absorptive backings while the noise introduced by variations in the depth-dependent light output and optical transport is larger. The results suggest that anisotropic imaging performance across the detector area can be incorporated into reconstruction algorithms for improving the image quality of breast tomosynthesis. This study also demonstrates that the assessment of image quality of breast tomosynthesis systems requires a more complete description of the detector response beyond local, center measurements of resolution and noise that assume some degree of symmetry in the detector performance.

SBRT of lung tumours: Monte Carlo simulation with PENELOPE of dose distributions including respiratory motion and comparison with different treatment planning systems

The purpose of this work was to simulate with the Monte Carlo (MC) code PENELOPE the dose distribution in lung tumours including breathing motion in stereotactic body radiation therapy (SBRT). Two phantoms were modelled to simulate a pentagonal cross section with chestwall (unit density), lung (density 0.3 g cm(-3)) and two spherical tumours (unit density) of diameters respectively of 2 cm and 5 cm. The phase-space files (PSF) of four different SBRT field sizes of 6 MV from a Varian accelerator were calculated and used as beam sources to obtain both dose profiles and dose-volume histograms (DVHs) in different volumes of interest. Dose distributions were simulated for five beams impinging on the phantom. The simulations were conducted both for the static case and including the influence of respiratory motion. To reproduce the effect of breathing motion different simulations were performed keeping the beam fixed and displacing the phantom geometry in chosen positions in the cranial and caudal and left-right directions. The final result was obtained by combining the different position with two motion patterns. The MC results were compared with those obtained with three commercial treatment planning systems (TPSs), two based on the pencil beam (PB) algorithm, the TMS-HELAX (Nucletron, Sweden) and Eclipse (Varian Medical System, Palo Alto, CA), and one based on the collapsed cone algorithm (CC), Pinnacle(3) (Philips). Some calculations were also carried out with the analytical anisotropic algorithm (AAA) in the Eclipse system. All calculations with the TPSs were performed without simulated breathing motion, according to clinical practice. In order to compare all the TPSs and MC an absolute dose calibration in Gy/MU was performed. The analysis shows that the dose (Gy/MU) in the central part of the gross tumour volume (GTV) is calculated for both tumour sizes with an accuracy of 2-3% with PB and CC algorithms, compared to MC. At the periphery of the GTV the TPSs overestimate the dose up to 10%, while in the lung tissue close to the GTV PB algorithms overestimate the dose and the CC underestimates it. When clinically relevant breathing motions are included in the MC simulations, the static calculations with the TPSs still give a relatively accurate estimate of the dose in the GTV. On the other hand, the dose at the periphery of the GTV is overestimated, compared to the static case.

Uncertainty estimation in intensity-modulated radiotherapy absolute dosimetry verification

PURPOSE

Intensity-modulated radiotherapy (IMRT) represents an important method for improving RT. The IMRT relative dosimetry checks are well established; however, open questions remain in reference dosimetry with ionization chambers (ICs). The main problem is the departure of the measurement conditions from the reference ones; thus, additional uncertainty is introduced into the dose determination. The goal of this study was to assess this effect systematically.

METHODS AND MATERIALS

Monte Carlo calculations and dosimetric measurements with five different detectors were performed for a number of representative IMRT cases, covering both step-and-shoot and dynamic delivery.

RESULTS

Using ICs with volumes of about 0.125 cm(3) or less, good agreement was observed among the detectors in most of the situations studied. These results also agreed well with the Monte Carlo-calculated nonreference correction factors (c factors). Additionally, we found a general correlation between the IC position relative to a segment and the derived correction factor c, which can be used to estimate the expected overall uncertainty of the treatment.

CONCLUSION

The increase of the reference dose relative standard uncertainty measured with ICs introduced by nonreference conditions when verifying an entire IMRT plan is about 1-1.5%, provided that appropriate small-volume chambers are used. The overall standard uncertainty of the measured IMRT dose amounts to about 2.3%, including the 0.5% of reproducibility and 1.5% of uncertainty associated with the beam calibration factor. Solid state detectors and large-volume chambers are not well suited to IMRT verification dosimetry because of the greater uncertainties. An action level of 5% is appropriate for IMRT verification. Greater discrepancies should lead to a review of the dosimetric procedure, including visual inspection of treatment segments and energy fluence.

Dose distributions in SBRT of lung tumors: Comparison between two different treatment planning algorithms and Monte-Carlo simulation including breathing motions

In stereotactic body radiotherapy (SBRT) of lung tumors, dosimetric problems arise from: 1) the limited accuracy in the dose calculation algorithms in treatment planning systems, and 2) the motions with the respiration of the tumor during treatment. Longitudinal dose distributions have been calculated with Monte Carlo simulation (MC), a pencil beam algorithm (PB) and a collapsed cone algorithm (CC) for two spherical lung tumors (2 cm and 5 cm diameter) in lung tissue, in a phantom situation. Respiratory motions were included by a convolution method, which was validated. In the static situation, the PB significantly overestimates the dose, relative to MC, while the CC gives a relatively accurate estimate. Four different respiratory motion patterns were included in the dose calculation with the MC. A "narrowing" of the longitudinal dose profile of up to 20 mm (at about 90% dose level) is seen relative the static dose profile calculated with the PB.

Monte Carlo simulation of MOSFET detectors for high-energy photon beams using the PENELOPE code

The aim of this work was the Monte Carlo (MC) simulation of the response of commercially available dosimeters based on metal oxide semiconductor field effect transistors (MOSFETs) for radiotherapeutic photon beams using the PENELOPE code. The studied Thomson&Nielsen TN-502-RD MOSFETs have a very small sensitive area of 0.04 mm(2) and a thickness of 0.5 microm which is placed on a flat kapton base and covered by a rounded layer of black epoxy resin. The influence of different metallic and Plastic water build-up caps, together with the orientation of the detector have been investigated for the specific application of MOSFET detectors for entrance in vivo dosimetry. Additionally, the energy dependence of MOSFET detectors for different high-energy photon beams (with energy >1.25 MeV) has been calculated. Calculations were carried out for simulated 6 MV and 18 MV x-ray beams generated by a Varian Clinac 1800 linear accelerator, a Co-60 photon beam from a Theratron 780 unit, and monoenergetic photon beams ranging from 2 MeV to 10 MeV. The results of the validation of the simulated photon beams show that the average difference between MC results and reference data is negligible, within 0.3%. MC simulated results of the effect of the build-up caps on the MOSFET response are in good agreement with experimental measurements, within the uncertainties. In particular, for the 18 MV photon beam the response of the detectors under a tungsten cap is 48% higher than for a 2 cm Plastic water cap and approximately 26% higher when a brass cap is used. This effect is demonstrated to be caused by positron production in the build-up caps of higher atomic number. This work also shows that the MOSFET detectors produce a higher signal when their rounded side is facing the beam (up to 6%) and that there is a significant variation (up to 50%) in the response of the MOSFET for photon energies in the studied energy range. All the results have shown that the PENELOPE code system can successfully reproduce the response of a detector with such a small active area.

The response of prototype plane-parallel ionization chambers in small megavoltage x-ray fields

Accurate small-field dosimetry has become important with the use of multiple small fields in modern radiotherapy treatments such as IMRT and stereotactic radiosurgery. In this study, we investigate the response of a set of prototype plane-parallel ionization chambers, based upon the Exradin T11 chamber, with active volume diameters of 2, 4, 10, and 20 mm, exposed to 6 MV stereotactic radiotherapy x-ray fields. Our goal was to assess their usefulness for accurate small x-ray field dose measurements. The relative ionization response was measured in circular fields (0.5 to 4 cm diameter) as compared to a 10 x 10 cm2 reference field. A large discrepancy (approximately 40%) was found between the relative response in the smallest plane-parallel chamber and other small volume dosimeters (radiochromic film, micro-metal-oxide-semiconductor field-effect transistor and diode) used for comparison. Monte Carlo BEAMnrc simulations were used to simulate the experimental setup in order to investigate the cause of the under-response and to calculate appropriate correction factors that could be applied to experimental measurements. It was found that in small fields, the air cavity of these custom-made research chambers perturbed the secondary electron fluence profile significantly, resulting in decreased fluence within the active volume, which in turn produces a chamber under-response. It is demonstrated that a large correction to the p(fl) correction factor would be required to improve dosimetric accuracy in small fields, and that these factors could be derived using Monte Carlo simulations.

An empirical method for the determination of wall perturbation factors for parallel-plate chambers in high-energy electron beams

The calibration of ion chambers in high-energy electron beams in terms of absorbed dose to water at the National Physical Laboratory requires knowledge of the ratio of perturbation factors in graphite and water phantoms. During a review of data required for the NPL calibration procedure an empirical model was developed to calculate the perturbation due to the rear wall, pwall, of a well-guarded ion chamber in a high-energy electron beam. The overall uncertainty in this method is estimated to be 0.4%, which is the lowest value reported to date. The model reproduces measured data at the 0.1% level or better and indicates that the NACP ion chamber has a nonzero perturbation factor in electron beams due to backscatter from the rear wall. The effect is small (<0.5%) at high energies (R50>4 cm, E0>10 MeV) but becomes large at low energies-up to 1.4% at E0=4 MeV (R50=1.2 cm). The model indicates that there is a nonzero correction for the NACP chamber in both a graphite and water phantom and that material adjacent to the air cavity has a significant effect on the measured ionization. These values are consistent with previous measurements and recent Monte Carlo calculations. The model could be used in the design of ion chambers and in the estimation of corrections for non-homogeneous systems, especially in the absence of accurate Monte Carlo simulations.

Anisotropic imaging performance in indirect x-ray imaging detectors

We report on the variability in imaging system performance due to oblique x-ray incidence, and the associated transport of quanta (both x rays and optical photons) through the phosphor, in columnar indirect digital detectors. The analysis uses MANTIS, a combined x-ray, electron, and optical Monte Carlo transport code freely available. We describe the main features of the simulation method and provide some validation of the phosphor screen models considered in this work. We report x-ray and electron three-dimensional energy deposition distributions and point-response functions (PRFs), including optical spread in columnar phosphor screens of thickness 100 and 500 microm, for 19, 39, 59, and 79 keV monoenergetic x-ray beams incident at 0 degrees, 10 degrees, and 15 degrees. In addition, we present pulse-height spectra for the same phosphor thickness, x-ray energies, and angles of incidence. Our results suggest that the PRF due to the phosphor blur is highly nonsymmetrical, and that the resolution properties of a columnar screen in a tomographic, or tomosynthetic imaging system varies significantly with the angle of x-ray incidence. Moreover, we find that the noise due to the variability in the number of light photons detected per primary x-ray interaction, summarized in the information or Swank factor, is somewhat independent of thickness and incidence angle of the x-ray beam. Our results also suggest that the anisotropy in the PRF is not less in screens with absorptive backings, while the noise introduced by variations in the gain and optical transport is larger. Predictions from MANTIS, after additional validation, can provide the needed understanding of the extent of such variations, and eventually, lead to the incorporation of the changes in imaging performance with incidence angle into the reconstruction algorithms for volumetric x-ray imaging systems.

Configuration of the electron transport algorithm of PENELOPE to simulate ion chambers

The stability of the electron transport algorithm implemented in the Monte Carlo code PENELOPE with respect to variations of its step length is analysed in the context of the simulation of ion chambers used in photon and electron dosimetry. More precisely, the degree of violation of the Fano theorem is quantified (to the 0.1% level) as a function of the simulation parameters that determine the step size. To meet the premises of the theorem, we define an infinite graphite phantom with a cavity delimited by two parallel planes (i.e., a slab) and filled with a 'gas' that has the same composition as graphite but a mass density a thousand-fold smaller. The cavity walls and the gas have identical cross sections, including the density effect associated with inelastic collisions. Electrons with initial kinetic energies equal to 0.01, 0.1, 1, 10 or 20 MeV are generated in the wall and in the gas with a uniform intensity per unit mass. Two configurations, motivated by the design of pancake- and thimble-type chambers, are considered, namely, with the initial direction of emission perpendicular or parallel to the gas-wall interface. This version of the Fano test avoids the need of photon regeneration and the calculation of photon energy absorption coefficients, two ingredients that are common to some alternative definitions of equivalent tests. In order to reduce the number of variables in the analysis, a global new simulation parameter, called the speedup parameter (a), is introduced. It is shown that setting a = 0.2, corresponding to values of the usual PENELOPE parameters of C1 = C2 = 0.02 and values of WCC and WCR that depend on the initial and absorption energies, is appropriate for maximum tolerances of the order of 0.2% with respect to an analogue, i.e., interaction-by-interaction, simulation of the same problem. The precise values of WCC and WCR do not seem to be critical to achieve this level of accuracy. The step-size dependence of the absorbed dose is explained in the light of the properties of PENELOPE's transport mechanics. This work is intended to help users to adopt an optimal configuration that guarantees both a high-accuracy calculation of the absorbed dose and a reasonably short computing time.

Electron stopping power and inelastic mean free path in amino acids and protein over the energy range of 20-20,000 eV

Systematic calculations of stopping power (SPs) and inelastic mean free path (IMFP) values for 20-20,000 eV electrons in a group of 15 amino acids and a simple protein have been performed. The calculations are based on the dielectric response model and take into account the exchange effect between the incident electron and target electrons. The optical energy-loss functions for the 15 investigated amino acids and the protein are evaluated by using an empirical approach, because of the lack of experimental optical data. For all the considered materials, the calculated mean ionization potentials are in good agreement with those given by Bragg's rule, and the evaluated SP values at 20 keV converge well to the Bethe-Bloch predictions. The data shown represent the first results of SP and IMFP, for these 15 amino acids and the protein in the energy range below 20 keV, and might be useful for studies of various radiation effects in these materials. In addition, the average energy deposited by inelastic scattering of the electrons on this group of 15 amino acids, on the protein, on Formvar and on DNA, respectively, has been estimated for energies below 20 keV. The dependences of the average energy deposition on the electron energy are given. These results are important for any detailed studies of radiation-induced inactivation of proteins and the DNA.

Fifty years of Monte Carlo simulations for medical physics

Monte Carlo techniques have become ubiquitous in medical physics over the last 50 years with a doubling of papers on the subject every 5 years between the first PMB paper in 1967 and 2000 when the numbers levelled off. While recognizing the many other roles that Monte Carlo techniques have played in medical physics, this review emphasizes techniques for electron-photon transport simulations. The broad range of codes available is mentioned but there is special emphasis on the EGS4/EGSnrc code system which the author has helped develop for 25 years. The importance of the 1987 Erice Summer School on Monte Carlo techniques is highlighted. As an illustrative example of the role Monte Carlo techniques have played, the history of the correction for wall attenuation and scatter in an ion chamber is presented as it demonstrates the interplay between a specific problem and the development of tools to solve the problem which in turn leads to applications in other areas.

Wall correction factors, Pwall, for parallel-plate ionization chambers

The EGSnrc Monte Carlo user-code CSnrc is used to calculate wall correction factors, Pwall,, for parallel-plate ionization chambers in photon and electron beams. A set of Pwall values, computed at the reference depth in water, is presented for several commonly used parallel-plate chambers. These values differ from the standard assumption of unity used by dosimetry protocols by up to 1.7% for clinical electron beams. Calculations also show that Pwall is strongly dependent on the depth of measurement and can vary by as much as 6% for a 6 MeV beam in moving from a depth of dref to a depth of R50. In photon beams, where there is limited information available regarding Pwall for parallel-plate chambers, CSnrc calculations show Pwall values of up to 2.4% at the reference depth over a range of photon energies. The Pwall values for photon beams are in good agreement with previous estimates of the wall correction but have much lower statistical uncertainties and cover a wider range of photon beam energies.

Monte Carlo correction factors for a Farmer 0.6 cm3ion chamber dose measurement in the build-up region of the 6 MV clinical beam

Reference dosimetry of photon fields is a well-established subject and currently available protocols (such as the IAEA TRS-398 and AAPM TG-51) provide methods for converting the ionization chamber (IC) reading into dose to water, provided reference conditions of charged particle equilibrium (CPE) are fulfilled. But these protocols cannot deal with the build-up region, where the lack of CPE limits the applicability of the cavity theorems and so the chamber correction factors become depth dependent. By explicitly including the IC geometry in the Monte Carlo simulations, depth-dependent dose correction factors are calculated for a PTW 30001 0.6 cm(3) ion chamber in the build-up region of the 6 MV photon beam. The corrected percentage depth dose (PDD) agrees within 2% with that measured using the NACP 02 plane-parallel ion chamber in the build-up region at depths greater than 0.4 cm, where the Farmer chamber wall reaches the phantom surface.

MANTIS: combined x-ray, electron and optical Monte Carlo simulations of indirect radiation imaging systems

We describe MANTIS (Monte carlo x-rAy electroN opTical Imaging Simulation), a tool for simulating imaging systems that tracks x-rays, electrons and optical photons in arbitrary materials and complex geometries. The x-ray and electron transport and involved physics models are from the PENELOPE package, and the optical transport and corresponding physics models are from DETECT-II and include Fresnel refraction and reflection at material boundaries, bulk absorption and scattering. Complex geometries can be handled with the aid of the geometry routines included in PENELOPE. When x-rays or electrons interact and deposit energy in the scintillator, the code generates a number of optical quanta according to a user-selected model for the conversion process. The optical photons are then tracked until they reach an absorption event, which in some cases contributes to the output signal, or escape from the geometry. We demonstrate the capabilities of this new tool with respect to the statistics of the optical signal detected and to the three-dimensional point-response functions corresponding to columnar phosphor screens.

Perturbation correction factors for the NACP-02 plane-parallel ionization chamber in water in high-energy electron beams

Recent dosimetry protocols for clinical high-energy electron beams recommend measurements of absorbed dose-to-water with a plane-parallel or cylindrical ionization chamber. For well-guarded plane-parallel ionization chambers, the ionization chamber perturbation factor in water, p(Q), has a recommended value of unity in all protocols. This assumption was investigated in detail in this study for one of the recommended ionization chambers in the protocols: the Scanditronix NACP-02 plane-parallel ionization chamber. Monte Carlo (MC) simulations of the NACP-02 ionization chamber with the EGSnrc code were validated against backscatter experiments. MC simulations were then used to calculate p(wall), p(cav) and p(Q) perturbation factors and water-to-air Spencer-Attix stopping powers in 4-19 MeV electron beams of a calibration laboratory (NPL), and in 6-22 MeV clinical electron beams from a Varian CL2300 accelerator. Differences between calculated and the currently recommended (Burns et al 1996 Med. Phys. 23 383-8) stopping powers, water-to-air, were found to be limited to 0.9% at depths between the reference depth z(ref) and the depth where the dose has decreased to 50% of the maximum dose, R50. p(wall) was found to exceed unity by 2.3% in the 4 MeV NPL calibration beam at z(ref). For higher energy electron beams p(wall) decreased to a value of about 1%. Combined with a p(cav) about 1% below unity for all energies at z(ref), this was found to cause p(Q) to exceed unity significantly for all energies. In clinical electron beams all three perturbation factors were found to increase with depth. Our findings indicate that the perturbation factors have to be taken into account in calibration procedures and for clinical depth dose measurements with the NACP-02 ionization chamber.