Abstract:Significance: Perturbation and differential Monte Carlo (pMC/dMC) methods, used in conjunction with nonlinear optimization methods, have been successfully applied to solve inverse problems in diffuse optics. Application of pMC to systems over a large range of optical properties requires optimal "placement" of baseline conventional Monte Carlo (cMC) simulations to minimize the pMC variance. The inability to predict the growth in pMC solution uncertainty with perturbation size limits the application of pMC, especially for multispectral datasets where the variation of optical properties can be substantial. Aim: We aim to predict the variation of pMC variance with perturbation size without explicit computation of perturbed photon weights. Our proposed method can be used to determine the range of optical properties over which pMC predictions provide sufficient accuracy. This method can be used to specify the optical properties for the reference cMC simulations that pMC utilizes to provide accurate predictions over a desired optical property range. Approach: We utilize a conventional error propagation methodology to calculate changes in pMC relative error for Monte Carlo simulations. We demonstrate this methodology for spatially resolved diffuse reflectance measurements with ±20% scattering perturbations. We examine the performance of our method for reference simulations spanning a broad range of optical properties relevant for diffuse optical imaging of biological tissues. Our predictions are computed using the variance, covariance, and skewness of the photon weight, path length, and collision distributions generated by the reference simulation. Results: We find that our methodology performs best when used in conjunction with reference cMC simulations that utilize Russian Roulette (RR) method. Specifically, we demonstrate that for a proximal detector placed immediately adjacent to the source, we can estimate the pMC relative error within 5% of the true value for scattering perturbations in the range of [-15%,+20%]. For a distal detector placed at ∼3 transport mean free paths relative to the source, our method provides relative error estimates within 20% for scattering perturbations in the range of [-8%,+15%]. Moreover, reference simulations performed at lower (μs'/μa) values showed better performance for both proximal and distal detectors. Conclusions: These findings indicate that reference simulations utilizing continuous absorption weighting (CAW) with the Russian Roulette method and executed using optical properties with a low (μs'/μa) ratio spanning the desired range of μs values, are highly advantageous for the deployment of pMC to obtain radiative transport estimates over a wide range of optical properties.

Probing the position-dependent optical energy fluence rate in three-dimensional scattering samples

Finite element modeling of light propagation in turbid media under illumination of a continuous-wave beam

Comparison Between the Pseudo-Spectral Time Domain Method and the Discrete Dipole Approximation for Light Scattering Simulations

A Simulation Study on the Flash X-Ray Spectra Spatial Distribution

Simulation of Light Scattering in Biological Tissues: the Coherent Component

3D spatially-resolved optical energy density enhanced by wavefront shaping

Adaptive handling of Rayleigh and Raman scatter of fluorescence data based on evaluation of the degree of spectral overlap

Analysis of relative error in perturbation Monte Carlo simulations of radiative transport

Poynting-Stokes tensor and radiative transfer in discrete random media: the microphysical paradigm

A Sequential Method for Measuring the Optical Properties of Two-Layer Media with Spatially-Resolved Diffuse Reflectance: Simulation Study

Estimation of Photon Distribution within Biological Tissue Using Monte Carlo Simulation

Universal Theory of Light Scattering of Randomly Oriented Particles: A Fluctuational-Electrodynamics Approach for Modeling of Light Transport in Disordered Nanostructures

Dependent scattering and plasmon coupling in concentrated suspensions of optically hard nanoparticles

Scattering of Electromagnetic Waves from 3d Multilayer Random Rough Surfaces Based on the Second-Order Small Perturbation Method: Energy Conservation, Reflectivity, and Emissivity

Experimental determination of effective light transport properties in fully anisotropic media

Spectral Determination of Μa, Μs and G of One Thick Turbid Sample from Three Scattered Light Signals

Multiple and dependent scattering by densely packed discrete spheres: Comparison of radiative transfer and Maxwell theory

A finite element approach for modeling photon transport in tissue

Characterization of the Optical Properties of Photoluminescent Turbid Media Using an Integrating Sphere and Monte Carlo Simulations

Fast finite difference solver for optical microscopy in deep biological tissue

A low-time complexity semi-analytic Monte Carlo radiative transfer model: Application to optical characteristics of complex spatial targets