Terahertz polarimetric imaging of biological tissues: Monte Carlo modeling of signal contrast mechanisms due to Mie scattering

Many promising biomedical applications have been proposed for terahertz (THz) spectroscopy and diagnostic imaging techniques. Polarimetric imaging systems are generally useful for enhancing imaging contrasts, yet the interplay between THz polarization changes and the random discrete structures in biological samples are not well understood. In this work, we performed Monte Carlo simulations of the propagation of polarized THz waves in skin and adipose tissues based on the Mie scattering from intrinsic structures, such as hair follicles or sweat glands. We show that the polarimetric contrasts are distinctly affected by concentration, size and dielectric properties of the scatterers, as well as the frequency and polarization of the incident THz waves. We describe the experimental requirements for observing and extracting these polarimetric signals due to the low energy and small angular spread of the back-scattered THz radiation. We analyzed the spatially integrated Mueller matrices of samples in the normal-incidence back-scattering geometry. We show that the frequency-dependent degree of polarization (DOP) can be used to infer the concentrations and dielectric contents of the scattering structures. Our modeling approach can be used to inform the design of the imaging modalities and the interpretation of the spectroscopic data in future terahertz biomedical imaging applications.