Monte Carlo simulation of source-excited in vivo x-ray fluorescence measurements of heavy metals

Development of a (170)Tm source for mercury monitoring studies in humans using XRF

The goals of the present study were to develop a (170)Tm radioisotope and generate a K XRF spectrum of mercury. Thulium foil and thulium oxide powder were both tested for impurities and the latter was found to be a better prospect for further studies. The (170)Tm radioisotope was developed from thulium oxide powder following the method of disolution and absorption. A suitable source holder and collimator were also designed based on Monte Carlo simulations. Using the radioisotope thus developed, a mercury XRF spectrum was successfully generated.

Uncertainty calculations for the measurement of in vivo bone lead by x-ray fluorescence

In order to quantify the bone lead concentration from an in vivo x-ray fluorescence measurement, typically two estimates of the lead concentration are determined by comparing the normalized x-ray peak amplitudes from the Kalpha(1) and Kbeta(1) features to those of the calibration phantoms. In each case, the normalization consists of taking the ratio of the x-ray peak amplitude to the amplitude of the coherently scattered photon peak in the spectrum. These two Pb concentration estimates are then used to determine the weighted mean lead concentration of that sample. In calculating the uncertainties of these measurements, it is important to include any covariance terms where appropriate. When determining the uncertainty of the lead concentrations from each x-ray peak, the standard approach does not include covariance between the x-ray peaks and the coherently scattered feature. These spectral features originate from two distinct physical processes, and therefore no covariance between these features can exist. Through experimental and simulated data, we confirm that there is no observed covariance between the detected Pb x-ray peaks and the coherently scattered photon signal, as expected. This is in direct contrast to recent work published by Brito (2006 Phys. Med. Biol. 51 6125-39). There is, however, covariance introduced in the calculation of the weighted mean lead concentration due to the common coherent normalization. This must be accounted for in calculating the uncertainty of the weighted mean lead concentration, as is currently the case. We propose here an alternative approach to calculating the weighted mean lead concentration in such a way as to eliminate the covariance introduced by the common coherent normalization. It should be emphasized that this alternative approach will only apply in situations in which the calibration line intercept is not included in the calculation of the Pb concentration from the spectral data: when the source of the intercept is well characterized and known to come from trace contamination by Pb in the plaster of Paris calibration standards. In our approach, the coherent normalization is only applied to one parameter and we no longer take a weighted mean of correlated quantities. Our proposed alternative calculation has essentially no effect on the calculated error of the mean lead concentration, indicating that the existing method of accounting for this covariance is sufficient.

Monte Carlo simulation of an anthropometric phantom used for calibrating in vivo K-XRF spectroscopy measurements of stable lead in bone

An anthropometric surrogate (phantom) of the human leg was defined in the constructs of the Monte Carlo N Particle (MCNP) code to predict the response when used in calibrating K x-ray fluorescence (K-XRF) spectrometry measurements of stable lead in bone. The predicted response compared favorably with measurements using the anthropometric phantom containing a tibia with increasing stable lead content. These benchmark measurements confirmed the validity of a modified MCNP code to accurately simulate K-XRF spectrometry measurements of stable lead in bone. A second, cylindrical leg phantom was simulated to determine whether the shape of the calibration phantom is a significant factor in evaluating K-XRF performance. Simulations of the cylindrical and anthropometric calibration phantoms suggest that a cylindrical calibration standard overestimates lead content of a human leg up to 4%. A two-way analysis of variance determined that phantom shape is a statistically significant factor in predicting the K-XRF response. These results suggest that an anthropometric phantom provides a more accurate calibration standard compared to the conventional cylindrical shape, and that a cylindrical shape introduces a 4% positive bias in measured lead values.

Modification to the Monte Carlo N-particle code for simulating direct, in vivo measurement of stable lead in bone

Monte Carlo N-Particle version 4C (MCNP4C) was used to simulate photon interactions associated with in vivo x-ray fluorescence (XRF) measurement of stable lead in bone. Experimental measurements, performed using a cylindrical anthropometric phantom (i.e., surrogate) of the human leg made from tissue substitutes for muscle and bone, revealed a significant difference between the intensity of the observed and predicted coherent backscatter peak. The observed difference was due to the failure of MCNP4C to simulate photon scatter associated with greater than six inverse angstroms of momentum transfer. The MCNP4C source code, photon directory, and photon library were modified to incorporate atomic form factors up to 7.1 inverse angstroms for the high Z elements defined in the K XRF simulation. The intensity of the predicted coherent photon backscatter peak at 88 keV using the modified code increased from 3.50 x 10(-9) to 8.59 x 10(-7) (roughly two orders of magnitude) and compares favorably with the experimental measurements.

A Monte Carlo (MC) based individual calibration method for in vivo x-ray fluorescence analysis (XRF)

X-ray fluorescence analysis (XRF) is a non-invasive method that can be used for in vivo determination of thyroid iodine content. System calibrations with phantoms resembling the neck may give misleading results in the cases when the measurement situation largely differs from the calibration situation. In such cases, Monte Carlo (MC) simulations offer a possibility of improving the calibration by better accounting for individual features of the measured subjects. This study investigates the prospects of implementing MC simulations in a calibration procedure applicable to in vivo XRF measurements. Simulations were performed with Penelope 2005 to examine a procedure where a parameter, independent of the iodine concentration, was used to get an estimate of the expected detector signal if the thyroid had been measured outside the neck. An attempt to increase the simulation speed and reduce the variance by exclusion of electrons and by implementation of interaction forcing was conducted. Special attention was given to the geometry features: analysed volume, source-sample-detector distances, thyroid lobe size and position in the neck. Implementation of interaction forcing and exclusion of electrons had no obvious adverse effect on the quotients while the simulation time involved in an individual calibration was low enough to be clinically feasible.

Monte Carlo simulations of in vivo K-shell X-ray fluorescence bone lead measurement and implications for radiation dosimetry

In order to improve measurement precision and decrease minimum detectable limit, recent applications of K-shell X-ray fluorescence (KXRF) bone lead measurement have used shorter source-to-sample (S-S) distances (approximately 0.5 cm) than the traditionally standard values ranging between 2.0 and 3.0 cm. This alteration will have an effect on subject radiation dose. This paper reports a comprehensive Monte Carlo study performed to investigate the radiation energy deposition values delivered to the leg of model human subjects of various ages. The simulations were run for models approximating 1-year, 5-year, and adult subjects, assuming lead concentrations of 10 microg/g in bone and tracing 500 million photons in each simulation. Trials were performed over a range of S-S distances, from 0.5 to 6.0 cm. The energy deposition due to photoelectric and Compton processes occurring in bone and soft tissue are presented. For each subject age, the Monte Carlo analysis demonstrates that the amount of energy deposited in the bone is increased as the sample is moved closer to the source (from 3.0 to 0.5 cm). The amount of energy deposited in the bone was found to increase by approximately 91% (1-year old), 66% (5-year old), and 41% (adult). The amount of energy deposited to the leg sample as a whole increased by approximately 43% (1-year old), 32% (5-year old), and 21% (adult). Results are used to estimate the changes in the amount of dose received by subjects of different ages.

Estimation of a method detection limit for an in vivo XRF arsenic detection system

An x-ray fluorescence measurement system has been developed with an 125I source to detect arsenic in superficial layers of phantoms and tissue. Based on in vivo measurements, in conjunction with Monte Carlo simulations, the detection limit for arsenic in skin ranges between 2.6+/-0.5 and 5.7+/-1.1 microg g(-1), depending on skin thickness and assuming that arsenic is uniformly distributed in the skin. The effect of skin arsenic distribution was also examined.

The feasibility of measuring silver concentrations in vivo with x-ray fluorescence

X-ray fluorescence (XRF) has been demonstrated to be an extremely useful technique for measuring trace quantities of heavy metals in various tissues in the body. This study investigates the applicability of XRF to the measurement of silver concentrations in skin. The system chosen employs an 125I source to excite the silver K x-rays, with the source, sample and detector arranged in a 90 degrees geometry. Experiments with silver-doped skin phantoms indicate that a minimum detectable concentration of 3-4 ppm is possible in 10-20 min measurement periods. Based on estimates of silver concentrations in the skin of patients suffering from argyria, the proposed system has sufficient sensitivity to warrant further investigation into its usefulness for non-invasive monitoring of exposed populations. Specifically, such a measurement may well allow for the identification of individuals at risk of subsequently exhibiting argyria, an irreversible skin pigmentation arising from silver exposure.

Monte Carlo simulation of x-ray spectra in diagnostic radiology and mammography using MCNP4C

The general purpose Monte Carlo N-particle radiation transport computer code (MCNP4C) was used for the simulation of x-ray spectra in diagnostic radiology and mammography. The electrons were transported until they slow down and stop in the target. Both bremsstrahlung and characteristic x-ray production were considered in this work. We focus on the simulation of various target/filter combinations to investigate the effect of tube voltage, target material and filter thickness on x-ray spectra in the diagnostic radiology and mammography energy ranges. The simulated x-ray spectra were compared with experimental measurements and spectra calculated by IPEM report number 78. In addition, the anode heel effect and off-axis x-ray spectra were assessed for different anode angles and target materials and the results were compared with EGS4-based Monte Carlo simulations and measured data. Quantitative evaluation of the differences between our Monte Carlo simulated and comparison spectra was performed using student's t-test statistical analysis. Generally, there is a good agreement between the simulated x-ray and comparison spectra, although there are systematic differences between the simulated and reference spectra especially in the K-characteristic x-rays intensity. Nevertheless, no statistically significant differences have been observed between IPEM spectra and the simulated spectra. It has been shown that the difference between MCNP simulated spectra and IPEM spectra in the low energy range is the result of the overestimation of characteristic photons following the normalization procedure. The transmission curves produced by MCNP4C have good agreement with the IPEM report especially for tube voltages of 50 kV and 80 kV. The systematic discrepancy for higher tube voltages is the result of systematic differences between the corresponding spectra.

Monte Carlo investigations of distance-dependent effects on energy deposition in K-shell x-ray fluorescence bone lead measurement

Radiation energy deposition results are presented from a Monte Carlo code simulating the lower part of a leg during an in vivo 109Cd K-shell x-ray fluorescence (KXRF) bone lead measurement. The simulations were run for a leg phantom model representing an adult subject, assuming concentrations of 10 microg Pb per gram bone mineral and tracing 500 million photons in each simulation. Trials were performed over a range (0.5-6.0 cm) of source-to-sample (S-S) distances. Energies deposited due to Compton and photoelectric processes occurring in the bone and the soft tissue were obtained. The data show an increase in the amount of energy deposited in the bone as the sample is moved closer to the source (from 2.0 cm to 0.5 cm). However, there is a decrease in the amount of energy deposited in the soft tissue as the sample is moved closer to the source over the same distance interval. In decreasing the S-S distance from 2.0 cm to 0.5 cm, the amount of energy deposited in the sample as a whole was found to increase by 11%. By calculating the energy deposition in the bone and in the soft tissue as a fraction of the total energy deposited in the sample, the corresponding changes are quantified as a function of S-S distance. Similarly, the proportions of energy deposited via the photoelectric effect and Compton scattering are presented as a function of S-S distance.

A comparison between EGS4 and MCNP computer modeling of an in vivo X-ray fluorescence system

The Monte Carlo computer codes EGS4 and MCNP were used to develop a theoretical model of a 180 degrees geometry in vivo X-ray fluorescence system for the measurement of platinum concentration in head and neck tumors. The model included specification of the photon source, collimators, phantoms and detector. Theoretical results were compared and evaluated against X-ray fluorescence data obtained experimentally from an existing system developed by the Swansea In Vivo Analysis and Cancer Research Group. The EGS4 results agreed well with the MCNP results. However, agreement between the measured spectral shape obtained using the experimental X-ray fluorescence system and the simulated spectral shape obtained using the two Monte Carlo codes was relatively poor. The main reason for the disagreement between the results arises from the basic assumptions which the two codes used in their calculations. Both codes assume a "free" electron model for Compton interactions. This assumption will underestimate the results and invalidates any predicted and experimental spectra when compared with each other.

Normalisation with coherent scatter signal: improvements in the calibration procedure of the 57Co-based in vivo XRF bone-Pb measurement

The feasibility of a normalised calibration variable to account for interpatient variability for in vivo 57Co XRF (X-ray fluorescence) finger bone-lead measurements was assessed. Normalising the lead X-ray intensities to the coherent scatter signal was investigated by experiment and Monte Carlo simulation. The X-ray to coherent ratios for a fixed lead concentration were within 5-10% of the mean, within uncertainty, over a physiologically relevant range of finger bone sizes and overlying tissue thicknesses. This is an acceptable level of variation to introduce, as it is less than the uncertainty of a typical in vivo measurement. Normalisation has several advantages compared with the current method of correcting for interpatient variation, as it eliminates the need for transporting extensive equipment to on-site measurements, reduces the subject dose by a factor of two, and increases the objectivity of the bone-Pb assessment.