Feasibility of the use of a pixelated transmission chamber to measure scattered-radiation in projection radiography: I. Monte Carlo analysis of x-ray scatter propagation and detector modeling.

In this paper, we propose the use of a pixelated transmission chamber, placed between the patient and the imaging detector, to measure the scatter component of a radiation beam impinging on said imaging detector. Using Monte Carlo simulation, a three-parameter model for the propagation of the scatter component in the transmitted beam is first developed. The use of the transmission chamber to determine the model parameters is then modeled, again using Monte Carlo simulation, and the feasibility of this approach is determined. The amount of radiation backscattered from the imaging detector into the transmission chamber was also calculated, for several separation distances between the two. It is shown that at a separation of 10cm, the amount of backscatter radiation is independent of the imaged object and that therefore it can be determined as part of a calibration procedure for the transmission chamber.

Stereotactic breast irradiation with kilovoltage x-ray beams.

The purpose of this work is to determine, using Monte Carlo simulation and a realistic patient model, the characteristics of the resultant absorbed dose distributions when breast tumors are irradiated using small-field stereotactic body radiation therapy (SBRT) with kilovoltage x-ray beams instead of the standard megavoltage energies currently in use. The Rensselaer Polytechnic Institute (RPI) female phantom was used to model a pair of small-field SBRT breast treatments: in one treatment the tumor at depth and another one with the tumor located close to the breast surface. Each treatment consisted of 300 circular beams aimed at the tumor from a plurality of positions. The PENELOPE Monte Carlo code was used to determine the absorbed dose distribution for each beam and subsequently an optimization algorithm determined each beam weight according to a set of prescription goals. Both kilo- and megavoltage beam treatments were modeled, the latter to be used as a reference. Cumulative dose-volume histograms for eleven structures were used to compare the kilovoltage and reference treatments. Integral dose values are also reported. Absorbed dose distributions for the target volumes as well as the organs at risk were within the parameters reported in a clinical trial for both treatments. While for the ipsilateral healthy breast tissue the megavoltage treatment does offer an advantage in terms of less volume irradiated to intermediate doses, for the contralateral structures, breast and lung, the low penetration ability of the kilovoltage treatment results in a lower maximum dose. Skin dose is higher for the kilovoltage treatment but still well within the tolerance limits reported in the clinical trial.

Feasibility of kilovoltage x-ray energy modulation by gaseous media and its application in contrast-enhanced radiotherapy.

PURPOSE: To present a method to modulate the energy contents of a kilovoltage x-ray beam that makes use of a gas as the modulating medium. The method is capable of producing arbitrary x-ray spectra by varying the pressure of the modulating gas and the peak kilovoltage (kVp) of the x-ray beams whose energy is being modulated. METHODS: An aluminum chamber was machined with a 0.5 cm wall thickness, designed to withstand pressures of more than 80 atm. A pressure sensor and electrovalves were used to monitor and regulate the gas pressure. Argon was used as the modulating gas. A CdTe spectrometer was used to measure x-ray spectra for different combinations of kVp and gas pressure, thus obtaining a set of basis x-ray functions. An arbitrary x-ray spectrum can then be formed by the linear combination of such basis functions. In order to show one possible application of the modulation method, a contrast-enhanced radiotherapy prostate treatment was optimized with respect to the x-ray beam energy, without restrictions on the possible shape of the resultant x-ray spectra. RESULTS: The x-ray spectra basis functions obtained display a smooth and gradual variation of their average energy as a function of the gas pressure for a given kVp, sometimes in the order of 1 or 2 keV. This gradual variation would be difficult to obtain with a conventional aluminum or copper filters, as the change in thickness necessary to reproduce the data presented would be in the order of micrometers, making necessary the use of a large number of such filters. Using the modulation method presented here, the authors were able to reconstruct the optimized x-ray spectra from the measured basis functions, for different optimization objectives. CONCLUSIONS: A method has been developed that allows for the controlled modulation of the energy contents of kilovoltage x-ray spectra. The method has been shown to be able to reproduce spectra of arbitrary shape, such as those obtained from the optimization of contrast-enhanced radiotherapy. The method may have other applications as well, such as in the precise matching of diagnostic x-ray catalog spectra.

Monte Carlo modeling of converging small-field contrast-enhanced radiotherapy of prostate.

Radiation therapy using a kilovoltage X-ray source to irradiate a target previously loaded with a radiological contrast agent, contrast-enhanced radiotherapy (CERT), has been shown both theoretically and in a preliminary experimental study to represent a potential alternative to high-energy treatments. It has also been shown, however, to produce an integral dose that can be up to twice that resulting from a conventional megavoltage treatment. In this work, using a realistic patient model and Monte Carlo simulation, a CERT prostate treatment plan is designed that makes use of a plurality of small circular beams aimed at the target in such a way as to minimize the radiological trajectory to the target volume. Gold nanoparticles are assumed to be the contrast agent. Two cases are examined, one with a concentration level in the target of 10 mg-Au per gram of tissue and the second with a concentration of 3 mg-Au per gram of tissue in the target. A background concentration of 1 mg of contrast agent per gram of tissue was assumed everywhere else in both cases. The Cimmino feasibility algorithm was then used to find each beam weight in order to obtain the prescribed target dose, set at 72 Gy to 100% of the tumor volume. It is shown that the approach using the small circular fields, a radiosurgery treatment, produces treatment plans with excellent absorbed dose distributions while at the same time it reduces by up to 60% the non-tumor integral dose imparted to the irradiated subject. A brief discussion on the technology necessary to clinically implement this treatment modality is also presented.

Microdosimetry of X-ray-irradiated gold nanoparticles.

The use of contrast agents, particularly those made of high atomic number elements like gold nanoparticles, to enhance the X-ray absorption properties of tissue has recently gained attention in the context of radiotherapy treatments. Because these contrast agents alter the secondary electron field in the irradiated medium by adding an Auger electron component, it is necessary to determine the change in the microdosimetric spectra brought about by the incorporation of such agents. Using Monte Carlo simulation, it is shown that the linear energy transfer and the beam quality factor in the vicinity of a gold nanoparticle irradiated with kilovoltage X-ray beams increase substantially when compared with irradiation without the gold nanoparticles present.

Monte Carlo modeling and optimization of contrast-enhanced radiotherapy of brain tumors.

Contrast-enhanced radiotherapy involves the use of a kilovoltage x-ray beam to impart a tumoricidal dose to a target into which a radiological contrast agent has previously been loaded in order to increase the x-ray absorption efficiency. In this treatment modality the selection of the proper x-ray spectrum is important since at the energy range of interest the penetration ability of the x-ray beam is limited. For the treatment of brain tumors, the situation is further complicated by the presence of the skull, which also absorbs kilovoltage x-ray in a very efficient manner. In this work, using Monte Carlo simulation, a realistic patient model and the Cimmino algorithm, several irradiation techniques and x-ray spectra are evaluated for two possible clinical scenarios with respect to the location of the target, these being a tumor located at the center of the head and at a position close to the surface of the head. It will be shown that x-ray spectra, such as those produced by a conventional x-ray generator, are capable of producing absorbed dose distributions with excellent uniformity in the target as well as dose differential of at least 20% of the prescribed tumor dose between this and the surrounding brain tissue, when the tumor is located at the center of the head. However, for tumors with a lateral displacement from the center and close to the skull, while the absorbed dose distribution in the target is also quite uniform and the dose to the surrounding brain tissue is within an acceptable range, hot spots in the skull arise which are above what is considered a safe limit. A comparison with previously reported results using mono-energetic x-ray beams such as those produced by a radiation synchrotron is also presented and it is shown that the absorbed dose distributions rendered by this type of beam are very similar to those obtained with a conventional x-ray beam.

Treatment planning considerations in contrast-enhanced radiotherapy: energy and beam aperture optimization.

It has been shown that the use of kilovoltage x-rays in conjunction with a contrast agent incorporated into the tumor can lead to acceptable treatment plans with regard to the absorbed dose distribution produced in the target as well as in the tissue and organs at risk surrounding it. In this work, several key aspects related to the technology and irradiation techniques necessary to clinically implement this treatment modality are addressed by means of Monte Carlo simulation. The Zubal phantom was used to model a prostate radiotherapy treatment, a challenging site due to the depth of the prostate and the presence of bony structures that must be traversed by the x-ray beam on its way to the target. It is assumed that the concentration levels of the enhancing agent present in the tumor are at or below 10 mg per 1 g of tissue. The Monte Carlo code PENELOPE was used to model a commercial x-ray tube having a tungsten target. X-ray energy spectra for several combinations of peak electron energy and added filtration were obtained. For each energy spectrum, a treatment plan was calculated, with the PENELOPE Monte Carlo code, by modeling the irradiation of the patient as 72 independent conformal beams distributed at intervals of 5° around the phantom in order to model a full x-ray source rotation. The Cimmino optimization algorithm was then used to find the optimum beam weight and energy for different treatment strategies. It is shown that for a target dose prescription of 72 Gy covering the whole tumor, the maximum rectal wall and bladder doses are kept below 52 Gy for the largest concentration of contrast agent of 10 mg per 1 g of tissue. It is also shown that concentrations of as little as 5 mg per 1 g of tissue also render dose distributions with excellent sparing of the organs at risk. A treatment strategy to address the presence of non-uniform distributions of the contrast agent in the target is also modeled and discussed.

Contrast-enhanced radiotherapy: feasibility and characteristics of the physical absorbed dose distribution for deep-seated tumors.

Radiotherapy using kilovoltage x-rays in conjunction with contrast agents incorporated into the tumor, gold nanoparticles in particular, could represent a potential alternative to current techniques based on high-energy linear accelerators. In this paper, using the voxelized Zubal phantom in conjunction with the Monte Carlo code PENELOPE to model a prostate cancer treatment, it is shown that in combination with a 360 degrees arc delivery technique, tumoricidal doses of radiation can be delivered to deep-seated tumors while still providing acceptable doses to the skin and other organs at risk for gold concentrations in the tumor within the range of 7-10 mg-Au per gram of tissue. Under these conditions and using a x-ray beam with 90% of the fluence within the range of 80-200 keV, a 72 Gy physical absorbed dose to the prostate can be delivered, while keeping the rectal wall, bladder, skin and femoral heads below 65 Gy, 55 Gy, 40 Gy and 30 Gy, respectively. However, it is also shown that non-uniformities in the contrast agent concentration lead to a severe degradation of the dose distribution and that, therefore, techniques to locally quantify the presence of the contrast agent would be necessary in order to determine the incident x-ray fluence that best reproduces the dosimetry obtained under conditions of uniform contrast agent distribution.

Monte Carlo-derived TLD cross-calibration factors for treatment verification and measurement of skin dose in accelerated partial breast irradiation.

Monte Carlo simulation was employed to calculate the response of TLD-100 chips under irradiation conditions such as those found during accelerated partial breast irradiation with the MammoSite radiation therapy system. The absorbed dose versus radius in the last 0.5 cm of the treated volume was also calculated, employing a resolution of 20 microm, and a function that fits the observed data was determined. Several clinically relevant irradiation conditions were simulated for different combinations of balloon size, balloon-to-surface distance and contents of the contrast solution used to fill the balloon. The thermoluminescent dosemeter (TLD) cross-calibration factors were derived assuming that the calibration of the dosemeters was carried out using a Cobalt 60 beam, and in such a way that they provide a set of parameters that reproduce the function that describes the behavior of the absorbed dose versus radius curve. Such factors may also prove to be useful for those standardized laboratories that provide postal dosimetry services.

Determination of radiotherapy X-ray spectra using a screen-film system.

A method to determine the X-ray spectrum delivered by a medical linear accelerator is presented. This method consists of an analytical calculation of the primary spectrum using the Schiff bremsstrahlung cross-section formula. A correction factor that accounts for the scatter component of the spectrum is estimated by comparing the signal in two screen-film systems to a theoretical prediction using a model of energy deposition in such detectors. The model makes use of the quantum absorption efficiency and the average energy deposited per interacting photon concepts. These two quantities are calculated by means of Monte Carlo simulations of the screen-film systems used. This method is capable of determining the spectrum as a function of the spatial position across a plane perpendicular to the beam central axis. It does not, however, render information about the direction cosines of the X-ray fluence crossing such a plane, a requirement in order to produce a full phase-space file that can be used in conjunction with a Monte Carlo dose calculation engine.

Monte Carlo microdosimetry of 188Re- and 131I-labelled anti-CD20.

The radiolabelled monoclonal antibody anti-CD20 has the property of binding to the CD20 antigen expressed on the cell surface of B-lymphocytes, thus making it a useful tool in the treatment of non-Hodgkin's lymphoma. In this work, the event-by-event Monte Carlo code NOREC is used to calculate the single-event distribution function f(1)(z) in the cell nucleus using the beta spectra of the (188)Re and (131)I radionuclides. The simulated geometry consists of two concentric spheres representing the nucleus and the cell surface embedded in a semi-infinite water medium. An isotropic point source was placed on the cell surface to simulate the binding of the anti-CD20 labelled with either (188)Re or (131)I. The simulations were carried out for two combinations of cell surface and nucleus radii. A method was devised that allows one to calculate the contribution of betas of energy greater than 1 MeV, which cannot be simulated by the NOREC code, to the single-event distribution function. It is shown that disregarding this contribution leads to an overestimation of the frequency-mean specific energy of the order of 9-12%. In general, the antibody radiolabelled with (131)I produces single-event distribution functions that yield higher frequency-mean specific energies.

Characteristics of the photoneutron contamination present in a high-energy radiotherapy treatment room.

The photoneutron contamination arising from a high-energy medical lineal accelerator is calculated using Monte Carlo simulation, as a function of the radiation field size. The information is used to model the neutron propagation in a radiotherapy treatment room and the transmission across concrete mazes. The Monte Carlo code MCNP4C is used to model the main components of a medical lineal accelerator. Simulations were performed to calculate the photoneutron yields and spectra as a function of the radiation field size. The yield of contaminant photoneutrons is observed to increase with the size of the radiation beam, but the energy spectra remain the same, suggesting that the contamination arises from above the movable collimator. The transport of the photoneutrons across a treatment room corroborates the validity of empirical models, but the transmission across a concrete maze produces a dose-equivalent tenth-value layer that differs from previous data.