The Monte Carlo modelling of in vivo X-ray fluorescence measurement of lead in tissue

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.

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.

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

This paper reports on the Monte Carlo simulation of in vivo x-ray fluorescence (XRF) measurements. Our model is an improvement on previously reported simulations in that it relies on a theoretical basis for modelling Compton momentum broadening as well as detector efficiency. Furthermore, this model is an accurate simulation of experimentally detected spectra when comparisons are made in absolute counts; preceding models have generally only achieved agreement with spectra normalized to unit area. Our code is sufficiently flexible to be applied to the investigation of numerous source-excited in vivo XRF systems. Thus far the simulation has been applied to the modelling of two different systems. The first application was the investigation of various aspects of a new in vivo XRF system, the measurement of uranium in bone with 57Co in a backscatter (approximately 180 degrees) geometry. The Monte Carlo simulation was critical in assessing the potential of applying XRF to the measurement of uranium in bone. Currently the Monte Carlo code is being used to evaluate a potential means of simplifying an established in vivo XRF system, the measurement of lead in bone with 57Co in a 90 degrees geometry. The results from these simulations may demonstrate that calibration procedures can be significantly simplified and subject dose may be reduced. As well as providing an excellent tool for optimizing designs of new systems and improving existing techniques, this model can be used in the investigation of the dosimetry of various XRF systems. Our simulation allows a detailed understanding of the numerous processes involved when heavy metal concentrations are measured in vivo with XRF.

Development of the specific purpose Monte Carlo code CEARXRF for the design and use of in vivo X-ray fluorescence analysis systems for lead in bone

X-ray fluorescence (XRF) systems have been increasingly used for in vivo toxic trace-element analysis in the human body, such as lead in the tibia. Monte Carlo simulation can provide an efficient and flexible method for designing and using in vivo XRF systems. The Monte Carlo code CEARXRF has been developed specifically to simulate the complete pulse height spectrum of energy-dispersive XRF systems. This code is capable of tracking photons in a general geometry and modelling all of the physics of photon interactions in the energy range 1-150 keV for elements Z = 1-94, including primary and higher degree excitations of K and L XRF, the Doppler broadening of Compton-scattered photon energies, and the polarization effects in low-energy photon scatterings. The scattering background for minimum detectable concentration (MDC) analysis may be simulated more accurately by taking into account Doppler broadening in the distribution of the Compton-scattered photon energy due to electron-binding effects. The use of polarized excitation photons has been shown to be important in producing a low scattering background and good measurement sensitivity. The code has two very unique and important features: (1) complete composition and density correlated sampling that is extremely useful for studying measurement sensitivity to small changes in sample composition and density; and (2) Monte Carlo library spectra calculation for the determination of elemental amounts by the Monte Carlo-Library Least-Squares (MCLLS) method. The capability of CEARXRF to aid the design and optimization of in vivo XRF analysis has been verified by modelling hypothesized lead K and L XRF measurement systems.

Optimization of in vivo X-ray fluorescence analysis methods for bone lead by simulation with the Monte Carlo code CEARXRF

In the design of X-ray fluorescence (XRF) systems for in vivo measurements of lead in human bone, the most important considerations are the minimum detectable concentration (MDC), and accuracy and precision. Possible design optimizations can be investigated much more easily and economically by Monte Carlo simulation than by experiment. The specific purpose Monte Carlo code CEARXRF has been used in the present study for: (1) improving the MDC of a hypothesized in vivo 109Cd source-based KXRF system and a 109Cd source or X-ray tube source-based LXRF system by investigating the effects of source polarization and source-bone-detector geometry modification on reducing the scattering background, and (2) investigating the effects of sample variables, such as overlying skin thickness on the MDC and the lead XRF intensity precision. In addition, the feasibility of the Monte Carlo-Library Least-Squares (MCLLS) approach has been investigated in a preliminary fashion for 109Cd-based KXRF spectroscopy analysis.