Polarised X-rays in XRF-analysis for improved in vivo detectability of cadmium in man

A technique is described for the in vivo XRF-analysis of cadmium in the kidney cortex of man using plane polarised photons for excitation. The polarised photons are produced by scattering the radiation from an X-ray tube (W anode, 150 kV, 15 mA) in a polymethylmethacrylate disc at a 90 deg angle. The beam paths (X-ray tube to scatterer, scatterer to sample, sample to detector) must represent three mutually orthogonal directions. The minimum detectable concentration for a counting time of 1800 s and a skin-kidney distance of 30 mm is 8 micrograms g-1. This is a factor of 2.5 lower than our earlier method with direct excitation using the 59.5 keV photons from 241Am. The energy imparted has also been lowered from 0.4 to 0.2 mJ. The cadmium concentration in the kidney cortex of six occupationally exposed persons varied between 15 and 170 micrograms g-1.     

Characteristics of low-angle x-ray scattering from some biological samples

Design of medical imaging devices based on the detection of low-angle coherent scattering is a subject of increasing interest. The technique is based on the differences in the distribution of photons coherently scattered from different body tissues. Coherent scattering is also useful in monitoring changes that may occur in a healthy tissue (e.g. carcinoma). In this work, low angle scattering properties of some tissues and tissue-equivalent materials are studied. Special care is given to the possibility of distinguishing between tissues of similar water content (e.g. muscle and blood). For this purpose, a Monte Carlo simulation is updated, introducing molecular form factor data, which include molecular interference effects. This program is used to simulate the angular distribution of scattered photons from two tissue-equivalent materials (lucite and water) and three biological samples (muscle, fat and blood). Simulation results agree well with previously measured angular distributions of scattered photons at 59.54 keV. Scattering from water and lucite is also measured at 8.047 keV. The effects of scattering geometry, sample thickness, incident photon energy and tissue type on the angular distribution of scattered photons are investigated. Results reveal the potential of measuring the full width at half maximum (FWHM) of the scattered photon distribution for tissue characterization. Energies up to 13 keV and sample thickness of 0.3 cm reported maximum differences between investigated samples. These conditions are expected to maximize the potential of using coherent scattering set-ups to monitor changes in biological samples even if their water contents are similar. Present results may act as a guide for the optimization of coherent scattering imaging systems.     

Compton spectroscopy in the diagnostic x-ray energy range. I. Spectrometer design

The optimal design of a Compton spectrometer for measuring photon energy spectra from x-ray tubes in a clinical laboratory is analysed. The demands are: (i) coherent and multiple scattering distort the measurements and must be avoided; (ii) the measuring time should be as short as possible to avoid unnecessary wear on the x-ray tube; and (iii) the impairment in energy resolution due to the scattering geometry should be kept minimal. A scattering angle of 90 degrees is advocated. Scatterers (of low-atomic-number material) in the shape of long circular rods (0.5-4 mm diameter, 20-40 mm long) are preferable to scattering foils. Use of a short focus-scatterer distance (approximately 200 mm) is to be preferred compared to using a large detector area (greater than or equal to 4 mm diameter) in order to establish a sufficiently high count rate in the detector. Short focal distances and a 90 degrees scattering angle are advantages in measuring energy spectra in the gantry of CT machines where the available space is limited. To limit the geometrical energy broadening to less than 1 keV, the spread in scattering angles of registered photons must not exceed 1-2 degrees for incident photon energies of 100-150 keV.     

Monte Carlo simulation of diagnostic x-ray scatter

A Monte Carlo method was developed and implemented to simulate x-ray photon transport. Simulations consisted of a pencil beam of monoenergetic photons with energies from 50 to 110 keV incident on water and aluminum slabs. The dependence of scatter fraction and multiple scattering on x-ray energy, scatterer thickness, and material is reported in both number and energy fluence. The average energy of scattered photons reaching the detector plane is also reported. Comparisons are made to previous x-ray scatter computations.     

Photon backscattering tissue characterization by energy dispersive spectroscopy evaluations

Techniques for in vivo tissue characterization based on scattered photons have usually been confined to evaluating coherent and Compton peaks. However, information can also be obtained from the energy analysis of the Compton scattered distribution. This paper looks at the extension of a technique validated by the authors for characterizing tissues composed of low-atomic-number elements. To this end, an EDXRS (energy dispersive x-ray spectrometry) computer simulation procedure was performed and applied to test the validity of a figure of merit able to characterize binary compounds. This figure of merit is based on the photon fluence values in a restricted energy interval of the measured distribution of incoherently scattered photons. After careful experimental tests with 59.54 keV incident photons at scattering angles down to 60degrees, the simulation procedure was applied to quasi-monochromatic and polychromatic high-radiance sources. The results show that the characterization by the figure of merit, which operates satisfactorily with monochromatic sources, is unsatisfactory in the latter cases, which seem to favour a different parameter for compound characterization.     

The significance of electron binding corrections in Monte Carlo photon transport calculations

Many Monte Carlo simulations ignore coherent scattering events and utilise the Klein-Nishina free electron distribution, rather than the incoherent differential cross-section, for choosing the trajectories of incoherently scattered photons. We assess the accuracy of this model by comparing its results with those of the complete bound electron model (form factor approach), which simulates coherent scattering events, and uses the appropriate bound electron angular scattering distributions. Both analytic and Monte Carlo calculations demonstrate that use of the free electron scattering distributions significantly underestimates the angular distribution of scattered photon energy resulting from low and medium energy photons incident upon carbon, iron, and platinum barriers. In using the free electron approximations to calculate barrier transmission, significant errors occur only for primary photon energies below 100 keV. Implementation of the complete bound electron model reduces the computational efficiency of our Monte Carlo code by only 10-25%.     

Dose measurements in inhomogeneous bone/tissue and lung/tissue phantoms for angiography using synchrotron radiation

For angiography using synchrotron radiation we measured the absorbed dose distribution in inhomogeneous phantoms with thin LiF:Mg, Cu, P, LiF:Mg, Ti thermoluminescent dosimeters (TLDs) in tissue and lung substitutes, and with Mg2SiO4:Tb TLDs in bone substitute for 33.32 keV monoenergetic photons from synchrotron radiation. The energy responses of the TLDs were measured in air for 10-40 keV monoenergetic photons. The values at 30 keV became smaller by 30% for LiF:Mg, Cu, P and larger by 22% for Mg2SiO4:Tb than the ratio of the mass energy absorption coefficients of the TLDs to that of air. These values were used to modify the calculated response of the TLDs in each phantom material. The absorbed dose distribution obtained was compared with that calculated using the Monte Carlo transport code EGS4 expanded to a low-energy region, and their agreement was confirmed taking linear polarization into account. In the bone substitute the dose increased by a factor of 3.9, while behind the bone the dose decreased drastically because of photon attenuation. In the lung substitute a slight dose difference from that in soft tissue was observed because of its different density. The LiF:Mg, Cu, P TLDs exhibited a better energy response, higher sensitivity and wider linear regions than did the other tissue-equivalent TLDs in the low-energy region.     

Beam collimation with polycapillary x-ray optics for high contrast high resolution monochromatic imaging

Monochromatic imaging can provide better contrast and resolution than conventional broadband radiography. In broadband systems, low energy photons do not contribute to the image, but are merely absorbed, while high energy photons produce scattering that degrades the image. By tuning to the optimal energy, one can eliminate undesirable lower and higher energies. Monochromatization is achieved by diffraction from a single crystal. A crystal oriented to diffract at a particular energy, in this case the characteristic line energy, diffracts only those photons within a narrow range of angles. The resultant beam from a divergent source is nearly parallel, but not very intense. To increase the intensity, collimation was performed with polycapillary x-ray optics, which can collect radiation from a divergent source and redirect it into a quasi parallel beam. Contrast and resolution measurements were performed with diffracting crystals with both high and low angular acceptance. Testing was first done at 8 keV with an intense copper rotating anode x-ray source, then 17.5 keV measurements were made with a low power molybdenum source. At 8 keV, subject contrast was a factor of five higher than for the polychromatic case. At 17.5 keV, monochromatic contrast was two times greater than the conventional polychromatic contrast. The subject contrasts measured at both energies were in good agreement with theory. An additional factor of two increase in contrast, for a total gain of four, is expected at 17.5 keV from the removal of scatter. Scatter might be simply removed using an air gap, which does not degrade resolution with a parallel beam.     

Simplified model of low energy x-ray backscattering

The reflection of x-rays from a half space is studied within the framework of a model that assumes multiple isotropic scattering of photons without energy loss. An exactly solvable analytical expression for the angular distribution of reflected photons is derived. The range of validity of the model was determined by the Monte Carlo simulation thereby incorporating energy loss and angular dependences. For water as a scatterer, in the energy range from 10 to 60 keV, which is often used in x-ray diagnostics, the two approaches differ by at most 5%. The analytic results, confirmed by the Monte Carlo simulation, show that the angular distribution of reflected photons for energies greater than 30 keV--where multiple scattering events dominate--may be represented by a cosine law, within a few per cent of accuracy.     

A generalised algorithm for spectral reconstruction in Compton spectroscopy with corrections for coherent scattering

Reconstruction of primary-photon energy spectra from pulse-height distributions obtained in a Compton spectrometer has earlier been performed under the assumption that coherent scattering in the scatterer is negligible. This holds for most clinical x-ray units operated in the range 40-150 kV. In mammography, and to some extent in dental radiography, the relatively high frequency of low-energy photons (less than 30 keV) in the primary beam makes it necessary to extend the algorithms to allow for significant contribution of coherent scattering. This extension is performed as a perturbation calculation to the algorithms developed earlier in which a modified Klein-Nishina scattering cross section was taken as the total scattering cross section. Comparison with energy spectra measured in the primary beam indicates that the Compton spectrometer with the extended algorithm is an excellent instrument for measuring energy spectra with energies down to a few keV.     

Bone composition measured by x-ray scattering.

Ten composite samples consisting of cortical bone and adipose tissue, in known proportions, were made. The intensity of monochromatic x-rays (energy 8 keV) scattered by these samples was determined as a function of the modulus of the scattering vector, K. The ratio of the heights of peaks at K values of around 134 and 22 nm-1 provided a measure of the ratio of adipose tissue to bone mineral in these samples. This method was then used to determine the ratio of adipose tissue to mineral in samples of trabecular bone from 16 vertebral bodies. The results were correlated with measurements of the bone composition determined by ashing (r = 0.66) and histomorphometry (r = 0.66). Furthermore, the ashing and histomorphometry results were correlated with each other (r = 0.68). The feasibility of using higher energy x-rays (35-80 keV) for obtaining the same information from bone within the body is briefly discussed.     

Quantitative assessment of bone mineral by photon scattering: accuracy and precision considerations

A method to determine the bone mineral density of the calcaneum has been reported earlier by our laboratory. In this method, the calcaneum is irradiated by a 60-keV photon beam from 241Am source and both the coherent and Compton scattered photons are detected by a high-purity Ge detector. The bone mineral density is determined by measuring the ratio of coherent-to-Compton scattered photons. The accuracy and the precision (in vitro) of the method are reported in this paper. The accuracy was determined to be 5%. This was obtained by comparing the bone mineral density values of cadaver calcanea measured directly by Archimedes' volume displacement method with the values measured by the scattering method. The precision was determined to be 3% by measuring the bone mineral density of a calibration phantom intermittently over a ten-month period by the scattering method.     

Rapid determination of cancellous bone mineral loss in ovariectomized rats by a subtraction technique

We present a rapid technique for determining cancellous bone mineral changes in small experimental animals. We used the distal centimeter of the right femur from ovariectomized (OX) (N = 30) and shamovariectomized (ShOX) (N = 28) rats, aged 90 days at surgery and killed at times from 125-540 days postsurgery. We used dual photon absorptiometry to scan the segment three times: intact, after parasagittal splitting, and after removing all cancellous bone. We equated the difference between the second and third scans to cancellous bone mineral content (Cn.BMC). To validate this, we compared it with histomorphometrically determined bone volume (BV/TV) of the proximal tibial metaphysis of the same rat. Parasagittally splitting the segment removed no detectable mineral. OX rats had 40% less Cn.BMC than ShOX rats. However, OX rats had 80% lower BV/TV than ShOX rats. The subtraction technique not only makes a rapid, reasonable assessment of cancellous bone loss in OX rats but permits a smaller sample size than histomorphometry. The histomorphometric technique finds a greater difference between OX and ShOX rats because it examines a region where cancellous bone loss is more marked than does the scanning technique. The current technique measures bone of not only the central secondary spongiosa but also the juxtacortical region and the primary spongiosa, where OX-related differences are less prominent. The principles of this subtraction technique proved workable. However, for the future, we recommend a two-scan technique using a dual energy X-ray scanner. It is likely to take only 20-30 min per specimen to assess cancellous bone mineral.     

The determination of electron density by the dual-energy Compton scatter method

The intensity of Compton scattered photons is directly proportional to the electron density of a substance; therefore it is possible to detect changes in electron density by this noninvasive method. The source-detector geometry has been varied to study the effects of geometry on the calculated electron density values, using a 90 degrees scattering angle. The dual-energy Compton scatter method provides a determination of electron density in known samples, with an uncertainty of less than 1% for the optimal geometry. The effective linear attenuation coefficient for Compton scatter has been determined with an uncertainty of 2.1%-3.6% for the optimal geometry used. The applications of the dual-energy Compton scattering technique in radiotherapy dosimetry and skeletal bone densitometry are being investigated.     

X-ray scattering from human breast tissues and breast-equivalent materials

The angular distributions of photons scattered by human breast tissues (adipose and glandular) and by eight breast-equivalent materials (water, polymethylmethacrylate, nylon, polyethylene and four commercial breast-equivalent materials simulating different glandular-adipose proportions) have been measured at a photon energy of 17.44 keV (Kalpha-radiation of Mo). Transmission target geometry has been used with an acceptance of +/- 0.6 degrees and an uncertainty of approximately 7%. Experimental molecular form factors were extracted from diffraction patterns normalizing the number of scattered photons with theoretical data in regions where no structure is expected. Linear attenuation coefficients have been measured for all samples at this energy. The results for water, polymethylmethacrylate, nylon and adipose tissue agree with former reported data. The results for human breast tissues at low and medium scattering angle (1-25 degrees, corresponding to the momentum transfer region between 0.2 and 3 nm(-1)) differ from the breast-equivalent materials. The results for adipose tissue are similar to the corresponding values from commercial breast-equivalent materials while the results for glandular tissue are similar to those for water.     

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.     

An accelerator based system for in vivo neutron activation analysis measurements of manganese in human hand bones

Manganese (Mn) is an essential nutrient for growth and development. Unfortunately, overexposure can lead to neurological damage, which is manifested as a movement disorder marked by tremors. Preclinical symptoms have been found in populations occupationally exposed to the element, and it is suggested that in late stages of the disorder, removing the Mn exposure will not prevent symptoms from progressing. Hence, it is desirable to have a means of monitoring Mn body burden. In vivo neutron activation analysis (IVNAA) is a technique which allows the concentration of some elements to be determined within sites of the body without invasive procedures. Data in the literature suggests that the Mn concentration in bone is greater than other tissues, and that it may be a long term storage site following exposure. Therefore, using the McMaster KN-accelerator to produce neutrons through the 7Li(p,n)7Be reaction, the feasibility of IVNAA for measuring Mn levels in the human hand bone was investigated. Mn is activated through the 55Mn(n,gamma)56Mn reaction, and the 847 keV gamma rays emitted when 56Mn decays are measured outside the body using NaI(Tl) detectors. An optimal incident proton energy of 2.00 MeV was determined from indium foil and microdosimetry measurements. Hand phantom data suggest a minimum detectable limit of approximately 1.8 ppm could be achieved with a reasonably low dose of 50 mSv to the hand (normal manganese levels in the human hand are approximately 1 ppm). It is recommended the technique be developed further to make human in vivo measurements.     

Compton spectroscopy in the diagnostic x-ray energy range. II. Effects of scattering material and energy resolution

The overall performance of a Compton spectrometer and, in particular, its energy resolution are investigated both experimentally and theoretically for different scattering materials. Using low-Z (less than or equal to 8) scatterers of moderate sizes (scatterer diameter d less than or equal to 5 mm), there are negligible disturbances due to coherent and/or multiple scattering at 90 degrees scattering angle and photon energies above 20 keV. Two factors contribute to decreasing the energy resolution compared with that in direct measurements: (i) the velocity distribution of the electrons in the scatterer and (ii) the scattering geometry. Of these, (i) is dominant for photon energies less than or equal to 100 keV. The optimal scattering material is a metal of as low Z as possible, i.e. beryllium. However, polyethylene and lucite are normally sufficiently good scatterers. The scattering geometry may become the dominating factor decreasing energy resolution at high photon energies hv greater than or equal to 150 keV.     

A particle track-repeating algorithm for proton beam dose calculation

A particle track-repeating algorithm has been developed for proton beam dose calculation for radiotherapy. Monoenergetic protons with 250 MeV kinetic energy were simulated in an infinite water phantom using the GEANT3 Monte Carlo code. The changes in location, angle and energy for every transport step and the energy deposition along the track were recorded for the primary protons and all secondary particles. When calculating dose for a patient with a realistic proton beam, the pre-generated particle tracks were repeated in the patient geometry consisting of air, soft tissue and bone. The medium and density for each dose scoring voxel in the patient geometry were derived from patient CT data. The starting point, at which a proton track was repeated, was determined according to the incident proton energy. Thus, any protons with kinetic energy less than 250 MeV can be simulated. Based on the direction of the incident proton, the tracks were first rotated and for the subsequent steps, the scattering angles were simply repeated for air and soft tissue but adjusted properly based on the scattering power for bone. The particle step lengths were adjusted based on the density for air and soft tissue and also on the stopping powers for bone while keeping the energy deposition unchanged in each step. The difference in nuclear interactions and secondary particle generation between water and these materials was ignored. The algorithm has been validated by comparing the dose distributions in uniform water and layered heterogeneous phantoms with those calculated using the GEANT3 code for 120, 150, 180 and 250 MeV proton beams. The differences between them were within 2%. The new algorithm was about 13 times faster than the GEANT3 Monte Carlo code for a uniform phantom geometry and over 700 times faster for a heterogeneous phantom geometry.     

Derivation of linear attenuation coefficients from CT numbers for low-energy photons.

One can estimate photon attenuation properties from the CT number. In a standard method one assumes that the linear attenuation coefficient is proportional to electron density and ignores its nonlinear dependence on atomic number. When the photon energy is lower than about 50 keV, such as for brachytherapy applications, however, photoelectric absorption and Rayleigh scattering become important. Hence the atomic number must be explicitly considered in estimating the linear attenuation coefficient. In this study we propose a method to more accurately estimate the linear attenuation coefficient of low-energy photons from CT numbers. We formulate an equation that relates the CT number to the electron density and the effective atomic number. We use a CT calibration phantom to determine unknown coefficients in the equation. The equation with a given CT number is then solved for the effective atomic number, which in turn is used to calculate the linear attenuation coefficient for low-energy photons. We use the CT phantom to test the new method. The method significantly improves the standard method in estimating the attenuation coefficient at low photon energies (20 keV < or = E < or = 40 keV) for materials with high atomic numbers.