Ernesto Mainegra-Hing

Efficient x-ray tube simulations

This article describes an efficiency study of directional bremsstrahlung splitting (DBS) for x-ray tube modeling. DBS is shown to be up to five or six orders of magnitude more efficient at 50 or 135 kV tube potential than a simulation without splitting, and 60 times more efficient compared to uniform bremsstrahlung splitting. A methodology is presented to determine the optimum splitting number for a given situation using a second degree polynomial expression derived from theoretical considerations. Very large optimum splitting numbers are found for small fields (5 mm radius) at 1 m from the x-ray source, which are relevant for half-value layer (HVL) calculations and for simulations related to primary air kerma standards. Two approaches for the calculation of kerma at a plane and inside a volume using track-length estimation are implemented in BEAMnrc, a user-code from the EGSnrc Monte Carlo simulation system for photon and electron transport. A practical application of DBS to HVL calculations for a Comet MXR-320 x-ray tube is reported. The agreement with measured HVLs at different constant tube potentials is found to be better than 2.3%.

Calculation of photon energy deposition kernels and electron dose point kernels in water

Effects of changes in the physics of EGSnrc compared to EGS4/PRESTA on energy deposition kernels for monoenergetic photons and on dose point kernels for beta sources in water are investigated. In the diagnostic energy range, Compton binding corrections were found to increase the primary energy fraction up to 4.5% at 30 keV with a corresponding reduction of the scatter component of the kernels. Rayleigh scattered photons significantly increase the scatter component of the kernels and reduce the primary energy fraction with a maximum 12% reduction also at 30 keV where the Rayleigh cross section in water has its maximum value. Sampling the photo-electron angular distribution produces a redistribution of the energy deposited by primaries around the interaction site causing differences of up to 2.7 times in the backscattered energy fraction at 20 keV. Above the pair production threshold, the dose distribution versus angle of the primary dose component is significantly different from the EGS4 results. This is related to the more accurate angular sampling of the electron-positron pair direction in EGSnrc as opposed to using a fixed angle approximation in default EGS4. Total energy fractions for photon beams obtained with EGSnrc and EGS4 are almost the same within 0.2%. This fact suggests that the estimate of the total dose at a given point inside an infinite homogeneous water phantom irradiated by broad beams of photons will be very similar for kernels calculated with both codes. However, at interfaces or near boundaries results can be very different especially in the diagnostic energy range. EGSnrc calculated kernels for monoenergetic electrons (50 keV, 100 keV, and 1 MeV) and beta spectra (32P and 90Y) are in excellent agreement with reported EGS4 values except at 1 MeV where inclusion of spin effects in EGSnrc produces an increase of the effective range of electrons. Comparison at 1 MeV with an ETRAN calculation of the electron dose point kernel shows excellent agreement.

Calculations for plane-parallel ion chambers in 60Co beams using the EGSnrc Monte Carlo code.

The EGSnrc Monte Carlo simulation system is used to obtain, for 10 plane-parallel ionization chambers in 60Co beams, the correction factors Kcomp and Pwall that account for the nonequivalence of the chamber wall material to the buildup cap and the phantom material, respectively. A more robust calculation method has been used compared to that used in previous works. A minor conceptual error related to the axial nonuniformity correction factor, Kan, has been identified and shown to have an effect of about 0.2%. The assumption that Pwall in-phantom is numerically equal to Kcomp calculated for a water buildup cap is shown to be accurate to better than 0.06%, thereby justifying the use of Kcomp calculations which are much more efficient. The effect on the calculated dose to the air in the cavity of the particle production threshold and transport energies used in the simulations is studied. Uncertainties in the calculated correction factors due to uncertainties in the photon and electron cross-section data are studied. They are 0.14% and 0.24%, respectively (1 standard deviation), for Kcomp factors. The uncertainties on Kwall factors are 0.03% from photon cross-section uncertainties and negligible from electron cross-section uncertainties. A comparison with previous EGS4/PRESTA calculations shows that present results are systematically higher by an average of 0.8%, ranging from 0.4% up to 1.4%. The present results are in better agreement with reported experimental values.

Optimizing non‐Pb radiation shielding materials using bilayers

PURPOSE The objective of this study was to demonstrate that the weight of non-Pb radiation shielding materials can be minimized by structuring the material as a bilayer composed of different metal-powder-embedded elastomer layers. METHODS Measurements and Monte Carlo (MC) calculations were performed to study the attenuation properties of several non-Pb metal bilayers over the x-ray energy range 30-150 keV. Metals for the layers were chosen on the basis of low cost, nontoxicity, and complementary photoelectric absorption characteristics. The EGSnrc user code cavity.cpp was used to calculate the resultant x-ray fluence spectra after attenuation by these metal layers. Air kerma attenuation was measured using commercially manufactured metal/elastomer test layers. These layers were irradiated using the primary standard calibration beams at the Institute for National Measurement Standards in Ottawa, Canada utilizing the six x-ray beam qualities recommended in the German Standard DIN 6857. Both the measurements and the calculations were designed to approximate surface irradiation as well as penetrating radiation at 10 mm depth in soft tissue. The MC modeling point and the position of the measurement detector for surface irradiation were both directly against the downstream face of the attenuating material, as recommended in DIN 6857. RESULTS The low-Z upstream/high-Z downstream ordering of the metal bilayers provided substantially more attenuation than the reverse order. Optimal percentages of each metal in each bilayer were determined for each x-ray radiation beam quality. CONCLUSIONS Depending on the x-ray quality, appropriate choices of two complementary metal-embedded elastomer layers can decrease the weight of radiation shielding garments by up to 25% compared to Pb-based elastomer garments while providing equivalent attenuation.

Variance reduction techniques for fast Monte Carlo CBCT scatter correction calculations

Several variance reduction techniques improving the efficiency of the Monte Carlo estimation of the scatter contribution to a cone beam computed tomography (CBCT) scan were implemented in egs_ctct, an EGSnrc-based application for CBCT-related calculations. The largest impact on the efficiency comes from the splitting + Russian Roulette techniques which are described in detail. The fixed splitting technique is outperformed by both the position-dependent importance splitting (PDIS) and the region-dependent importance splitting (RDIS). The superiority of PDIS over RDIS observed for a water phantom with bone inserts is not observed when applying these techniques to a more realistic human chest phantom. A maximum efficiency improvement of several orders of magnitude over an analog calculation is obtained. A scatter calculation combining the reported efficiency gain with a smoothing algorithm is already in the proximity of being of practical use if a medium size computer cluster is available.

On the accuracy of techniques for obtaining the calibration coefficient N{sub K} of {sup 192}Ir HDR brachytherapy sources

The accuracy of interpolation or averaging procedures for obtaining the calibration coefficient NK for Ir192 high-dose-rate brachytherapy sources has been investigated using the EGSnrc Monte Carlo simulation system. It is shown that the widely used two-point averaging procedure of Goetsch et al. [Med. Phys. 18, 462 (1991)] has some conceptual problems. Most importantly, they recommended, as did the IAEA, averaging AwallNK values whereas one should average 1∕NK values. In practice this and other issues are shown to have little effect except for Goetsch et al.s methods for determining Awall values. Their method of generalizing the Awall values measured in one geometry to other geometries is incorrect by up to 2%. However, these errors in Awall values cause systematic errors of only 0.3% in Ir192 calibration coefficients. It is shown that Awall values need not be included in the averaging technique at all, thereby simplifying the technique considerably. It is demonstrated that as long as ion chambers with a flat response are used and/or very heavily filtered 250kV (or higher) beams of x rays are used in the averaging, then almost all techniques can provide adequate accuracy.

Fast Monte Carlo calculation of scatter corrections for CBCT images

A very fast Monte Carlo algorithm for the calculation of the scatter contribution in cone beam computed tomography, implemented within the EGSnrc framework, is presented. Based on the combination of several variance reduction techniques, an efficiency improvement of three orders of magnitude over an analog simulation is obtained. A denoising algorithm applied to the computed scatter distribution is shown to further increase the efficiency of the calculation by about a factor of 10. An iterative scatter correction algorithm is proposed and its feasibility is demonstrated on three different phantoms.

Monte Carlo reference data sets for imaging research: Executive summary of the report of AAPM Research Committee Task Group 195

The use of Monte Carlo simulations in diagnostic medical imaging research is widespread due to its flexibility and ability to estimate quantities that are challenging to measure empirically. However, any new Monte Carlo simulation code needs to be validated before it can be used reliably. The type and degree of validation required depends on the goals of the research project, but, typically, such validation involves either comparison of simulation results to physical measurements or to previously published results obtained with established Monte Carlo codes. The former is complicated due to nuances of experimental conditions and uncertainty, while the latter is challenging due to typical graphical presentation and lack of simulation details in previous publications. In addition, entering the field of Monte Carlo simulations in general involves a steep learning curve. It is not a simple task to learn how to program and interpret a Monte Carlo simulation, even when using one of the publicly available code packages. This Task Group report provides a common reference for benchmarking Monte Carlo simulations across a range of Monte Carlo codes and simulation scenarios. In the report, all simulation conditions are provided for six different Monte Carlo simulation cases that involve common x-ray based imaging research areas. The results obtained for the six cases using four publicly available Monte Carlo software packages are included in tabular form. In addition to a full description of all simulation conditions and results, a discussion and comparison of results among the Monte Carlo packages and the lessons learned during the compilation of these results are included. This abridged version of the report includes only an introductory description of the six cases and a brief example of the results of one of the cases. This work provides an investigator the necessary information to benchmark his/her Monte Carlo simulation software against the reference cases included here before performing his/her own novel research. In addition, an investigator entering the field of Monte Carlo simulations can use these descriptions and results as a self-teaching tool to ensure that he/she is able to perform a specific simulation correctly. Finally, educators can assign these cases as learning projects as part of course objectives or training programs.

Patient-specific scatter correction in clinical cone beam computed tomography imaging made possible by the combination of Monte Carlo simulations and a ray tracing algorithm

Abstract Purpose. Cone beam computed tomography (CBCT) image quality is limited by scattered photons. Monte Carlo (MC) simulations provide the ability of predicting the patient-specific scatter contamination in clinical CBCT imaging. Lengthy simulations prevent MC-based scatter correction from being fully implemented in a clinical setting. This study investigates the combination of using fast MC simulations to predict scatter distributions with a ray tracing algorithm to allow calibration between simulated and clinical CBCT images. Material and methods. An EGSnrc-based user code (egs_cbct), was used to perform MC simulations of an Elekta XVI CBCT imaging system. A 60keV x-ray source was used, and air kerma scored at the detector plane. Several variance reduction techniques (VRTs) were used to increase the scatter calculation efficiency. Three patient phantoms based on CT scans were simulated, namely a brain, a thorax and a pelvis scan. A ray tracing algorithm was used to calculate the detector signal due to primary photons. A total of 288 projections were simulated, one for each thread on the computer cluster used for the investigation. Results. Scatter distributions for the brain, thorax and pelvis scan were simulated within 2% statistical uncertainty in two hours per scan. Within the same time, the ray tracing algorithm provided the primary signal for each of the projections. Thus, all the data needed for MC-based scatter correction in clinical CBCT imaging was obtained within two hours per patient, using a full simulation of the clinical CBCT geometry. Conclusions. This study shows that use of MC-based scatter corrections in CBCT imaging has a great potential to improve CBCT image quality. By use of powerful VRTs to predict scatter distributions and a ray tracing algorithm to calculate the primary signal, it is possible to obtain the necessary data for patient specific MC scatter correction within two hours per patient.

Hounsfield unit recovery in clinical cone beam CT images of the thorax acquired for image guided radiation therapy

A comprehensive artefact correction method for clinical cone beam CT (CBCT) images acquired for image guided radiation therapy (IGRT) on a commercial system is presented. The method is demonstrated to reduce artefacts and recover CT-like Hounsfield units (HU) in reconstructed CBCT images of five lung cancer patients. Projection image based artefact corrections of image lag, detector scatter, body scatter and beam hardening are described and applied to CBCT images of five lung cancer patients. Image quality is evaluated through visual appearance of the reconstructed images, HU-correspondence with the planning CT images, and total volume HU error. Artefacts are reduced and CT-like HUs are recovered in the artefact corrected CBCT images. Visual inspection confirms that artefacts are indeed suppressed by the proposed method, and the HU root mean square difference between reconstructed CBCTs and the reference CT images are reduced by 31% when using the artefact corrections compared to the standard clinical CBCT reconstruction. A versatile artefact correction method for clinical CBCT images acquired for IGRT has been developed. HU values are recovered in the corrected CBCT images. The proposed method relies on post processing of clinical projection images, and does not require patient specific optimisation. It is thus a powerful tool for image quality improvement of large numbers of CBCT images.