Frank Verhaegen

SpekCalc: a program to calculate photon spectra from tungsten anode x-ray tubes

A software program, SpekCalc, is presented for the calculation of x-ray spectra from tungsten anode x-ray tubes. SpekCalc was designed primarily for use in a medical physics context, for both research and education purposes, but may also be of interest to those working with x-ray tubes in industry. Noteworthy is the particularly wide range of tube potentials (40-300 kVp) and anode angles (recommended: 6-30 degrees) that can be modelled: the program is therefore potentially of use to those working in superficial/orthovoltage radiotherapy, as well as diagnostic radiology. The utility is free to download and is based on a deterministic model of x-ray spectrum generation (Poludniowski 2007 Med. Phys. 34 2175). Filtration can be applied for seven materials (air, water, Be, Al, Cu, Sn and W). In this note SpekCalc is described and illustrative examples are shown. Predictions are compared to those of a state-of-the-art Monte Carlo code (BEAMnrc) and, where possible, to an alternative, widely-used, spectrum calculation program (IPEM78).

Monte Carlo modelling of external radiotherapy photon beams

An essential requirement for successful radiation therapy is that the discrepancies between dose distributions calculated at the treatment planning stage and those delivered to the patient are minimized. An important component in the treatment planning process is the accurate calculation of dose distributions. The most accurate way to do this is by Monte Carlo calculation of particle transport, first in the geometry of the external or internal source followed by tracking the transport and energy deposition in the tissues of interest. Additionally, Monte Carlo simulations allow one to investigate the influence of source components on beams of a particular type and their contaminant particles. Since the mid 1990s, there has been an enormous increase in Monte Carlo studies dealing specifically with the subject of the present review, i.e., external photon beam Monte Carlo calculations, aided by the advent of new codes and fast computers. The foundations for this work were laid from the late 1970s until the early 1990s. In this paper we will review the progress made in this field over the last 25 years. The review will be focused mainly on Monte Carlo modelling of linear accelerator treatment heads but sections will also be devoted to kilovoltage x-ray units and 60Co teletherapy sources.

High-Resolution, Small Animal Radiation Research Platform With X-Ray Tomographic Guidance Capabilities

PURPOSE To demonstrate the computed tomography, conformal irradiation, and treatment planning capabilities of a small animal radiation research platform (SARRP). METHODS AND MATERIALS The SARRP uses a dual-focal spot, constant voltage X-ray source mounted on a gantry with a source-to-isocenter distance of 35 cm. Gantry rotation is limited to 120 degrees from vertical. X-rays of 80-100 kVp from the smaller 0.4-mm focal spot are used for imaging. Both 0.4-mm and 3.0-mm focal spots operate at 225 kVp for irradiation. Robotic translate/rotate stages are used to position the animal. Cone-beam computed tomography is achieved by rotating the horizontal animal between the stationary X-ray source and a flat-panel detector. The radiation beams range from 0.5 mm in diameter to 60 x 60 mm(2). Dosimetry is measured with radiochromic films. Monte Carlo dose calculations are used for treatment planning. The combination of gantry and robotic stage motions facilitate conformal irradiation. RESULTS The SARRP spans 3 ft x 4 ft x 6 ft (width x length x height). Depending on the filtration, the isocenter dose outputs at a 1-cm depth in water were 22-375 cGy/min from the smallest to the largest radiation fields. The 20-80% dose falloff spanned 0.16 mm. Cone-beam computed tomography with 0.6 x 0.6 x 0.6 mm(3) voxel resolution was acquired with a dose of <1 cGy. Treatment planning was performed at submillimeter resolution. CONCLUSION The capability of the SARRP to deliver highly focal beams to multiple animal model systems provides new research opportunities that more realistically bridge laboratory research and clinical translation.

Characterization of scattered radiation in kV CBCT images using Monte Carlo simulations

Kilovoltage (kV) cone beam computed tomography (CBCT) images suffer from a substantial scatter contribution. In this study, Monte Carlo (MC) simulations are used to evaluate the scattered radiation present in projection images. These predicted scatter distributions are also used as a scatter correction technique. Images were acquired using a kV CBCT bench top system. The EGSnrc MC code was used to model the flat panel imager, the phantoms, and the x-ray source. The x-ray source model was validated using first and second half-value layers (HVL) and profile measurements. The HVLs and the profile were found to agree within 3% and 6%, respectively. MC simulated and measured projection images for a cylindrical water phantom and for an anthropomorphic head phantom agreed within 8% and 10%. A modified version of the DOSXYZnrc MC code was used to score phase space files with identified scattered and primary particles behind the phantoms. The cone angle, the source-to-detector distance, the phantom geometry, and the energy were varied to determine their effect on the scattered radiation distribution. A scatter correction technique was developed in which the MC predicted scatter distribution is subtracted from the projections prior to reconstruction. Preliminary testing of the procedure was done with an anthropomorphic head phantom and a contrast phantom. Contrast and profile measurements were obtained for the scatter corrected and noncorrected images. An improvement of 3% for contrast between solid water and a liver insert and 11% between solid water and a Teflon insert were obtained and a significant reduction in cupping and streaking artifacts was observed.

Dual-energy CT-based material extraction for tissue segmentation in Monte Carlo dose calculations

Monte Carlo (MC) dose calculations are performed on patient geometries derived from computed tomography (CT) images. For most available MC codes, the Hounsfield units (HU) in each voxel of a CT image have to be converted into mass density (rho) and material type. This is typically done with a (HU; rho) calibration curve which may lead to mis-assignment of media. In this work, an improved material segmentation using dual-energy CT-based material extraction is presented. For this purpose, the differences in extracted effective atomic numbers Z and the relative electron densities rho(e) of each voxel are used. Dual-energy CT material extraction based on parametrization of the linear attenuation coefficient for 17 tissue-equivalent inserts inside a solid water phantom was done. Scans of the phantom were acquired at 100 kVp and 140 kVp from which Z and rho(e) values of each insert were derived. The mean errors on Z and rho(e) extraction were 2.8% and 1.8%, respectively. Phantom dose calculations were performed for 250 kVp and 18 MV photon beams and an 18 MeV electron beam in the EGSnrc/DOSXYZnrc code. Two material assignments were used: the conventional (HU; rho) and the novel (HU; rho, Z) dual-energy CT tissue segmentation. The dose calculation errors using the conventional tissue segmentation were as high as 17% in a mis-assigned soft bone tissue-equivalent material for the 250 kVp photon beam. Similarly, the errors for the 18 MeV electron beam and the 18 MV photon beam were up to 6% and 3% in some mis-assigned media. The assignment of all tissue-equivalent inserts was accurate using the novel dual-energy CT material assignment. As a result, the dose calculation errors were below 1% in all beam arrangements. Comparable improvement in dose calculation accuracy is expected for human tissues. The dual-energy tissue segmentation offers a significantly higher accuracy compared to the conventional single-energy segmentation.

Small animal radiotherapy research platforms

Advances in conformal radiation therapy and advancements in pre-clinical radiotherapy research have recently stimulated the development of precise micro-irradiators for small animals such as mice and rats. These devices are often kilovolt x-ray radiation sources combined with high-resolution CT imaging equipment for image guidance, as the latter allows precise and accurate beam positioning. This is similar to modern human radiotherapy practice. These devices are considered a major step forward compared to the current standard of animal experimentation in cancer radiobiology research. The availability of this novel equipment enables a wide variety of pre-clinical experiments on the synergy of radiation with other therapies, complex radiation schemes, sub-target boost studies, hypofractionated radiotherapy, contrast-enhanced radiotherapy and studies of relative biological effectiveness, to name just a few examples. In this review we discuss the required irradiation and imaging capabilities of small animal radiation research platforms. We describe the need for improved small animal radiotherapy research and highlight pioneering efforts, some of which led recently to commercially available prototypes. From this, it will be clear that much further development is still needed, on both the irradiation side and imaging side. We discuss at length the need for improved treatment planning tools for small animal platforms, and the current lack of a standard therein. Finally, we mention some recent experimental work using the early animal radiation research platforms, and the potential they offer for advancing radiobiology research.

Accuracy of the photon and electron physics in GEANT4 for radiotherapy applications

This work involves a validation of the photon and electron transport of the GEANT4 particle simulation toolkit for radiotherapy physics applications. We examine the cross sections and sampling algorithms of the three electromagnetic physics models in version 4.6.1 of the toolkit: Standard, Low-energy, and Penelope. The depth dose distributions in water for incident monoenergetic and clinical beams are compared to the EGSNRC results. In photon beam simulations, all three models agree with EGSNRC to within 2%, except for the buildup region. Larger deviations are found for incident electron beams, and the differences are affected by user-imposed electron step limitations. Particle distributions through thin layers of clinical target materials, and perturbation effects near high-Z and low-Z interfaces are also investigated. The electron step size artifacts observed in our studies indicate potential problems with the condensed history algorithm. A careful selection of physics processes and transport parameters is needed for optimum efficiency and accuracy.

Spectroscopic characterization of a novel electronic brachytherapy system

The Axxent developed by Xoft Inc. is a novel electronic brachytherapy system capable of generating x-rays up to 50 keV. These low energy photon-emitting sources merit attention not only because of their ability to vary the dosimetric properties of the radiation, but also because of the radiobiological effects of low energy x-rays. The objective of this study is to characterize the x-ray source and to model it using the Geant4 Monte Carlo code. Spectral and attenuation curve measurements are performed at various peak voltages and angles and the source is characterized in terms of spectrum and half-value layers (HVLs). Also, the effects of source variation and source aging are quantified. Bremsstrahlung splitting, phase-space scoring and particle-tagging features are implemented in the Geant4 code, which is bench-marked against BEAMnrc simulations. HVLs from spectral measurements, attenuation curve measurements and Geant4 simulations mostly agree within uncertainty. However, there are discrepancies between measurements and simulations for photons emitted on the source transverse plane (90 degrees).

Correction of CT artifacts and its influence on Monte Carlo dose calculations

Computed tomography (CT) images of patients having metallic implants or dental fillings exhibit severe streaking artifacts. These artifacts may disallow tumor and organ delineation and compromise dose calculation outcomes in radiotherapy. We used a sinogram interpolation metal streaking artifact correction algorithm on several phantoms of exact-known compositions and on a prostate patient with two hip prostheses. We compared original CT images and artifact-corrected images of both. To evaluate the effect of the artifact correction on dose calculations, we performed Monte Carlo dose calculation in the EGSnrc/DOSXYZnrc code. For the phantoms, we performed calculations in the exact geometry, in the original CT geometry and in the artifact-corrected geometry for photon and electron beams. The maximum errors in 6 MV photon beam dose calculation were found to exceed 25% in original CT images when the standard DOSXYZnrc/CTCREATE calibration is used but less than 2% in artifact-corrected images when an extended calibration is used. The extended calibration includes an extra calibration point for a metal. The patient dose volume histograms of a hypothetical target irradiated by five 18 MV photon beams in a hypothetical treatment differ significantly in the original CT geometry and in the artifact-corrected geometry. This was found to be mostly due to miss-assignment of tissue voxels to air due to metal artifacts. We also developed a simple Monte Carlo model for a CT scanner and we simulated the contribution of scatter and beam hardening to metal streaking artifacts. We found that whereas beam hardening has a minor effect on metal artifacts, scatter is an important cause of these artifacts.

Monte Carlo modelling of radiotherapy kV x-ray units

To obtain accurate information for absorbed dose calculations in water for kilovoltage x-rays, the photon spectrum, planar fluence and the angular distribution of the photons at the collimator exit of the x-ray unit have to be known. The only way to obtain this information is by Monte Carlo (MC) simulation. Compared with the situation for high-energy photons and electrons, where in recent years numerous papers have been devoted to MC modelling of complete clinical accelerator units, there is a lack of similar work for kV x-ray units. A reliable MC model for a kV x-ray unit would allow the output information to be used in a treatment planning system for regular and irregular treatment fields. Furthermore, with MC simulation, perturbation factors of dose-measuring devices, such as those specified in codes of practice, can be calculated. In this work, the MC code EGS4/BEAM was used to build realistic models of two complete x-ray units. The tungsten target, exit window, collimator, additional filtration and applicator were taken into account. For some aspects of the work, a comparison was made with the simulations from another MC code, MCNP4B. The contribution to the characteristic radiation from electron impact ionization and from the photoelectric effect of reabsorbed bremsstrahlung photons was studied. Calculated and measured photon fluence spectra in air and half-value layers for a Philips MCN410 tube were compared for several anode voltages and additional filtrations. Results from the two codes agreed well, and the agreement with measured spectra was found to be good for energies above 50 keV but rather less good below that energy. For a Siemens Stabilipan 2 Th300 x-ray tube, HVLs and dose distributions in water were compared with measurements for several clinical x-ray qualities. For most of the combinations of radiation qualities and applicators, good agreement was obtained, although there were also some cases where the agreement was not so good. Electron contamination and photon build-up at the water surface were studied using MC simulation. The influence of depth on the photon spectral distribution was investigated. Both EGS4/BEAM and MCNP4B, in their default versions, handle inadequately the production of characteristic x-rays. This was found to have only a minor influence on the calculated dosimetric quantities. Simulations with MCNP4B required the use of several variance reduction techniques in order to obtain results within reasonable calculation times.