Geant4 implementation of inter-atomic interference effect in small-angle coherent X-ray scattering for materials of medical interest

Simulation of gamma-ray shielding properties for materials of medical interest

Radiation that is emitted from unstable nuclei during radioactive decay is an important phenomenon to be used in large fields, and thus, radiation shielding properties are important especially for gamma rays. Thus, in the present work, the radiation shielding properties in terms of linear attenuation coefficients and some other quantities for medical interest materials of water, fat and bone have been obtained. The results simulated by Phy-X/PSD online code the energy range of 10−3–105 MeV.

Evaluation of simulators for x-ray speckle-based phase contrast imaging

X-ray phase contrast imaging (PCI) denotes a group of highly sensitive imaging techniques that permits imaging at scales ranging from nanoscopic to the medical. Recently introduced, speckle-based imaging has seen a rapid development because of its experimental simplicity and its capability to retrieve the refraction, the scattering and the absorption of a sample using a conventional x-ray set-up. Precise simulation would permit to optimise the imaging setups for different applications, but until now works on simulation of x-ray speckle-based PCI have been very few. In this work we evaluate different simulation codes, based on Monte-Carlo, analytical ray-tracing and wave-optics Fresnel propagation. The simulation results are compared to both synchrotron and conventional imaging experiments to permits their validation. We obtain a strong similarity between simulated and experimental data. We discuss the validity and applicability of each approach.

Accuracy of the independent atom approximation in digital tomosynthesis Monte Carlo simulations

Monte Carlo (MC) codes serve as the gold standard simulation tool during design and optimisation of x-ray imaging systems. Such simulations often model Rayleigh scattering based on the Independent Atom Approximation Model (IAM). This model neglects the low range molecular interference (MI) effects of non-crystalline materials such as human tissues. Previous work has found discrepancies in the simulated images of planar x-ray images between IAM and MI models. However, insignificant differences were found for computed tomography (CT) reconstructions. In this work we present Geant4 MC simulations of a flat panel source digital tomosynthesis (DT) system for human extremities. Results show that with a 1:9 scatter to primary ratio (SPR) in the x-ray projections, the DT reconstructions are insensitive to the differences of the IAM and MI models. Therefore, MC codes that use the IAM model are sufficient for the study of DT systems. That is because DT algorithms have a larger effect on image quality than the few percent change in the noise due to a physical model and noise suppression methods make this change even less important. Dependency of this conclusion on SPR must be considered in other DT modalities where SPR might be larger.

Impact of bowtie filter and detector collimation on multislice CT scatter profiles: A simulation study

PURPOSE

To investigate via Monte Carlo simulations, the impact of scan subject size, antiscatter grid (ASG), collimator size, and bowtie filter on the distribution of scatter radiation in a typical realistically modeled third generation 16 slice diagnostic computed tomography (CT) scanner.

METHODS

Full radiation transport was simulated with Geant4 in a realistic CT scanner geometric model, including the imaging phantom, bowtie filter (BTF), collimators and detector assembly, except for the ASGs. An analytical method was employed to quantify the probable transmission through the ASG of each photon intersecting the detector array. Normalized scatter profiles (NSP) and scatter-to-primary-ratio (SPR) profiles were simulated for 90 and 140 kVp beams for different size phantoms and slice thicknesses. The impact of CT scatter on the reconstructed attenuation coefficient factor was also studied as were the modulating effects of phantom- and patient-tissue heterogeneities on scatter profiles. A method to characterize the relative spatial frequency content of sinogram signals was developed to assess the latter.

RESULTS

For the 21.4-cm diameter phantom, NSP and SPR increase linearly with collimator opening for both tube potentials, with the 90 kVp scan exhibiting slightly larger NSP and SPR. The BTF modestly modulates scatter under the phantom center, reducing the prominent off-axis lobes by factors of 1.1-1.3. The ASG reduces scatter on the central axis NSP threefold, and reduces scatter at the detectors outside the phantom shadow by factors of 25 to 500. For the phantoms with diameters of 27 and 32 cm, the scatter increases roughly three- and fourfold, respectively, demonstrating that scatter monotonically increases with phantom size, despite deployment of the ASG and BTF. In the absence of a scan subject, the ASG reduces the signal profile arising photons scattered by the BTF. Without ASG, the in-air scatter profile is relatively flat compared to the scatter profile when the ASG is present. For both 90 and 140 kVp photon spectra, the calculated attenuation coefficient decreases linearly with increasing collimation size. For both homogeneous and heterogeneous objects, NSPs are dominated by low spatial frequency content compared to the primary signal. However, the SPR, which quantifies the local magnitude of nonlinear detector response and is dominated by the high frequency content of the primary profile, can contribute strongly to high-spatial frequency streaking artifacts near high-density structures in reconstructed image artifacts.

CONCLUSION

Public-domain Monte Carlo codes, Geant-4 in particular, is a feasible method for characterizing CT detector response to scattered- and off-focal radiation. Our study demonstrates that the ASG substantially reduces the scatter radiation and reshapes scatter-radiation profiles and affects the accuracy with which the detector array can measure narrow-beam attenuation due its inability to distinguish between true uncollided primary and narrow-angle coherently scattered photons. Hence, incorporating the impact of detector array collimation into the forward-projection signal formation models used by iterative reconstruction algorithms is necessary to use CT for accurately characterizing material properties. While tissue heterogeneities exercise a modest influence on local NPS shape and magnitude, they do not add significant high spatial frequency content.

X-ray Fokker-Planck equation for paraxial imaging

The Fokker-Planck equation can be used in a partially-coherent imaging context to model the evolution of the intensity of a paraxial x-ray wave field with propagation. This forms a natural generalisation of the transport-of-intensity equation. The x-ray Fokker-Planck equation can simultaneously account for both propagation-based phase contrast, and the diffusive effects of sample-induced small-angle x-ray scattering, when forming an x-ray image of a thin sample. Two derivations are given for the Fokker-Planck equation associated with x-ray imaging, together with a Kramers-Moyal generalisation thereof. Both equations are underpinned by the concept of unresolved speckle due to unresolved sample micro-structure. These equations may be applied to the forward problem of modelling image formation in the presence of both coherent and diffusive energy transport. They may also be used to formulate associated inverse problems of retrieving the phase shifts due to a sample placed in an x-ray beam, together with the diffusive properties of the sample. The domain of applicability for the Fokker-Planck and Kramers-Moyal equations for paraxial imaging is at least as broad as that of the transport-of-intensity equation which they generalise, hence the technique is also expected to be useful for paraxial imaging using visible light, electrons and neutrons.