Carole K. Hayakawa

Perturbation Monte Carlo methods to solve inverse photon migration problems in heterogeneous tissues

We introduce a novel and efficient method to provide solutions to inverse photon migration problems in heterogeneous turbid media. The method extracts derivative information from a single Monte Carlo simulation to permit the rapid determination of rates of change in the detected photon signal with respect to perturbations in background tissue optical properties. We then feed this derivative information to a nonlinear optimization algorithm to determine the optical properties of the tissue heterogeneity under examination. We demonstrate the use of this approach to solve rapidly a two-region inverse problem of photon migration in the transport regime, for which diffusion-approximation-based approaches are not applicable.

Perturbation and differential Monte Carlo methods for measurement of optical properties in a layered epithelial tissue model.

The use of perturbation and differential Monte Carlo (pMC/dMC) methods in conjunction with nonlinear optimization algorithms were proposed recently as a means to solve inverse photon migration problems in regionwise heterogeneous turbid media. We demonstrate the application of pMC/dMC methods for the recovery of optical properties in a two-layer extended epithelial tissue model from experimental measurements of spatially resolved diffuse reflectance. The results demonstrate that pMC/dMC methods provide a rapid and accurate approach to solve two-region inverse photon migration problems in the transport regime, that is, on spatial scales smaller than a transport mean free path and in media where optical scattering need not dominate absorption. The pMC/dMC approach is found to be effective over a broad range of absorption (50 to 400%) and scattering (70 to 130%) perturbations. The recovery of optical properties from spatially resolved diffuse reflectance measurements is examined for different sets of source-detector separation. These results provide some guidance for the design of compact fiber-based probes to determine and isolate optical properties from both epithelial and stromal layers of superficial tissues.

Use of the δ-P 1 approximation for recovery of optical absorption, scattering, and asymmetry coefficients in turbid media

We introduce a robust method to recover optical absorption, reduced scattering, and single-scattering asymmetry coefficients (microa, micros, g1) of infinite turbid media over a range of (micros/microa) spanning 3 orders of magnitude. This is accomplished through the spatially resolved measurement of irradiance at source-detector separations spanning 0.25-8 transport mean free paths (l*). These measurements are rapidly processed by a multistaged nonlinear optimization algorithm in which the measured irradiances are compared with predictions given by the delta-P1 variant of the diffusion approximation to the Boltzmann transport equation. The ability of the delta-P1 model to accurately describe radiative transport within media of arbitrary albedo and on spatial scales comparable to l* is the key element enabling the separation of g1 from micros.

Analysis of single Monte Carlo methods for prediction of reflectance from turbid media

Starting from the radiative transport equation we derive the scaling relationships that enable a single Monte Carlo (MC) simulation to predict the spatially- and temporally-resolved reflectance from homogeneous semi-infinite media with arbitrary scattering and absorption coefficients. This derivation shows that a rigorous application of this single Monte Carlo (sMC) approach requires the rescaling to be done individually for each photon biography. We examine the accuracy of the sMC method when processing simulations on an individual photon basis and also demonstrate the use of adaptive binning and interpolation using non-uniform rational B-splines (NURBS) to achieve order of magnitude reductions in the relative error as compared to the use of uniform binning and linear interpolation. This improved implementation for sMC simulation serves as a fast and accurate solver to address both forward and inverse problems and is available for use at http://www.virtualphotonics.org/.

Quantitative spectroscopy of superficial turbid media

We report a novel diffuse optical spectroscopy probe design for determining optical properties of superficial volumes of turbid samples. The fiber-based probe employs a highly scattering layer placed in contact with the sample of interest. This layer diffuses photons from a collimated light source before they enter the sample and provides a basis for describing light transported in superficial media by the diffusion approximation. We compare the performance of this modified two-layer diffusion model with Monte Carlo simulations. A set of experiments that demonstrate the feasibility of this method in turbid tissue phantoms is also presented. Optical properties deduced by this approach are in good agreement with those derived by use of a benchmark method for determining optical properties. The average interrogation depth of the probe design investigated here is estimated to be less than 1 mm.

Electric field Monte Carlo simulations of focal field distributions produced by tightly focused laser beams in tissues

The focal field distribution of tightly focused laser beams in turbid media is sensitive to optical scattering and therefore of direct relevance to image quality in confocal and nonlinear microscopy. A model that considers both the influence of scattering and diffraction on the amplitude and phase of the electric field in focused beam geometries is required to describe these distorted focal fields. We combine an electric field Monte Carlo approach that simulates the electric field propagation in turbid media with an angular-spectrum representation of diffraction theory to analyze the effect of tissue scattering properties on the focal field. In particular, we examine the impact of variations in the scattering coefficient (µs), single-scattering anisotropy (g), of the turbid medium and the numerical aperture of the focusing lens on the focal volume at various depths. The model predicts a scattering-induced broadening, amplitude loss, and depolarization of the focal field that corroborates experimental results. We find that both the width and the amplitude of the focal field are dictated primarily by µs with little influence from g. In addition, our model confirms that the depolarization rate is small compared to the amplitude loss of the tightly focused field.

Investigation of a probe design for facilitating the uses of the standard photon diffusion equation at short source-detector separations: Monte Carlo simulations

We design a special diffusing probe to investigate the optical properties of human skin in vivo. The special geometry of the probe enables a modified two-layer (MTL) diffusion model to precisely describe the photon transport even when the source-detector separation is shorter than 3 mean free paths. We provide a frequency domain comparison between the Monte Carlo model and the diffusion model in both the MTL geometry and conventional semiinfinite geometry. We show that using the Monte Carlo model as a benchmark method, the MTL diffusion theory performs better than the diffusion theory in the semiinfinite geometry. In addition, we carry out Monte Carlo simulations with the goal of investigating the dependence of the interrogation depth of this probe on several parameters including source-detector separation, sample optical properties, and properties of the diffusing high-scattering layer. From the simulations, we find that the optical properties of samples modulate the interrogation volume greatly, and the source-detector separation and the thickness of the diffusing layer are the two dominant probe parameters that impact the interrogation volume. Our simulation results provide design guidelines for a MTL geometry probe.

Comparative analysis of discrete and continuous absorption weighting estimators used in Monte Carlo simulations of radiative transport in turbid media

We examine the relative error of Monte Carlo simulations of radiative transport that employ two commonly used estimators that account for absorption differently, either discretely, at interaction points, or continuously, between interaction points. We provide a rigorous derivation of these discrete and continuous absorption weighting estimators within a stochastic model that we show to be equivalent to an analytic model, based on the radiative transport equation (RTE). We establish that both absorption weighting estimators are unbiased and, therefore, converge to the solution of the RTE. An analysis of spatially resolved reflectance predictions provided by these two estimators reveals no advantage to either in cases of highly scattering and highly anisotropic media. However, for moderate to highly absorbing media or isotropically scattering media, the discrete estimator provides smaller errors at proximal source locations while the continuous estimator provides smaller errors at distal locations. The origin of these differing variance characteristics can be understood through examination of the distribution of exiting photon weights.

Estimating optical properties in layered tissues by use of the Born approximation of the radiative transport equation

We use the Born approximation of the radiative transport equation to recover simultaneously the absorption and scattering coefficients in a single layer of a two-layer tissue sample from reflectance data. This method reduces the estimation of both optical properties to a single linear, least-squares problem. It is valid over length scales smaller than a transport mean free path and hence is useful for epithelial tissue layers. We demonstrate the accuracy of this method by using spatially resolved reflectance data computed with Monte Carlo simulations.

Determination of Optical Properties of Superficial Volumes of Layered Tissue Phantoms

Previously, we reported the design of a new diffusing probe that employs a standard two-layer diffusion model to recover the optical properties of turbid samples. This particular probe had a source-detector separation of 2.5 mm and performance was validated with Monte Carlo simulations and homogeneous phantom experiments. The goal of the current study is to characterize the performance of this new method in the context of two-layer phantoms that mimic the optical properties of human skin. We analyze the accuracy of the recovered top layer optical properties and their dependences on the thickness of the top layer of two-layer phantoms. Our results demonstrate that the optical properties of the top layer can be accurately determined with a 1.6 mm source-detector separation diffusing probe when this layer thickness is as thin as 1 mm. Monte Carlo simulations illustrate that the interrogation depth can be further decreased by shortening the source-detector separation.