V. L. Kuzmin

Simulation of optical coherence tomography images by Monte Carlo modeling based on polarization vector approach.

Monte Carlo method is applied for simulation of 2D optical coherence tomography (OCT) images of skin-like model. Layer boundaries in skin model feature curved shape which agrees with physiological structure of human skin. The effect of coherence properties of probing radiation on OCT image formation and speckles in the detected OCT signal is considered. The developed model is employed for image simulation both for conventional and polarization dependent time-domain OCT modalities. Simulation of polarized OCT signal is performed using vector approach developed previously for modeling of electromagnetic field transfer in turbid media.

Simulation of polarization-sensitive optical coherence tomography images by a Monte Carlo method

We introduce a new Monte Carlo (MC) method for simulating optical coherence tomography (OCT) images of complex multilayered turbid scattering media. We demonstrate, for the first time of our knowledge, the use of a MC technique to imitate two-dimensional polarization-sensitive OCT images with nonplanar boundaries of layers in the medium like a human skin. The simulation of polarized low-coherent optical radiation is based on the vector approach generalized from the iterative procedure of the solution of Bethe-Saltpeter equation. The performances of the developed method are demonstrated both for conventional and polarization-sensitive OCT modalities.

Propagation and scattering of light in fluctuating media

Abstract The monograph deals with the problems of the propagation and scattering of light in molecular media. The explicit statistical mechanical averaging procedure for the equations of electrodynamics is developed. It permits to transform the molecular level description into the macroscopic one for the electrodynamics of the fluctuating media. In the framework of such an approach, the problems of the molecular correlation contribution into the dielectric permeability, of the calculation of the reflection coefficients with an account of surface layers and of the multiple light scattering are considered. The developed theory is applied to the description of the critical opalescence, the coherent backscattering enhancement, the light scattering depolarization phenomena and the propagation and scattering of light in anisotropic media, including the case of liquid crystals.

Monte Carlo simulation of coherent effects in multiple scattering

Using a combination of the stochastic Monte Carlo technique and the iteration procedure of the solution to the Bethe–Salpeter equation, it has been shown that the simulation of the optical path of a photon packet undergoing an nth scattering event directly corresponds to the nth–order ladder diagram contribution. In this paper, the Monte Carlo technique is generalized for the simulation of the coherent back–scattering and temporal correlation function of optical radiation scattered within the randomly inhomogeneous turbid medium. The results of simulation demonstrate a good agreement with the diffusing wave theory and experimental results.

Coherent multiple scattering effects and Monte Carlo method

Based on the comparison of the iteration procedure of solving the Bethe-Salpeter equation and the Monte Carlo method, we developed a method for simulating coherent multiple-scattering effects within the framework of a unified stochastic approach. The time correlation function and the interference component were calculated for the coherent backscattering from a multiply scattering medium.

Diffuse photon density wave measurements and Monte Carlo simulations.

Abstract. Diffuse photon density wave (DPDW) methodology is widely used in a number of biomedical applications. Here, we present results of Monte Carlo simulations that employ an effective numerical procedure based upon a description of radiative transfer in terms of the Bethe–Salpeter equation. A multifrequency noncontact DPDW system was used to measure aqueous solutions of intralipid at a wide range of source–detector separation distances, at which the diffusion approximation of the radiative transfer equation is generally considered to be invalid. We find that the signal–noise ratio is larger for the considered algorithm in comparison with the conventional Monte Carlo approach. Experimental data are compared to the Monte Carlo simulations using several values of scattering anisotropy and to the diffusion approximation. Both the Monte Carlo simulations and diffusion approximation were in very good agreement with the experimental data for a wide range of source–detector separations. In addition, measurements with different wavelengths were performed to estimate the size and scattering anisotropy of scatterers.

Modeling of photon density waves in the frequency domain

We have described the transfer of modulated radiation in a random medium in terms of the Bethe-Salpeter equation. Based on the obtained expression for the scattering intensity, we have developed an original technique of modeling the photon density waves in terms of the Monte Carlo method. Expressions for measurable parameters in the frequency domain have been derived, and, based on them, the amplitude and phase of the photon density waves have been calculated. We have studied how the parameters of the photon density waves depend on the scattering anisotropy for model states with the Henyey-Greenstein phase function. The range of applicability of the diffusion approximation for the interpretation of signals of photon density waves has been investigated.

Coherent Backscattering of Circularly Polarized Light from a Disperse Random Medium

To describe propagation of polarized electromagnetic wave within a disperse random medium a new Monte Carlo based technique with an adopted vector formalism has been developed. The technique has been applied for simulation of coherent backscattering of circularly polarized optical radiation from a random scattering medium. It has been found that the sign of helicity of circular polarized light does not change for a medium of point-like scatterers and can change signiflcantly for the scatterers with the higher anisotropy. We conclude that the helicity ∞ip of the circular polarized light can be observed in the tissue-like media. We flnd that this phenomenon manifests itself in case of limited number of scattering events and, apparently, can be attributed to the pulse character of incident radiation rather than to the speciflc form of scattering particles.

Numerical simulation of photon density waves in a biophantom with foreign objects

Based on the expression obtained for the scattering intensity in terms of the Bethe–Salpeter equation, the simulation of diffuse photon density wave (DPDW) is implemented in the context of the Monte Carlo method in random media imitating biological tissues in the presence of foreign objects. The DPDW amplitude and phase are calculated as a function of the source–detector distance with explicit allowance for internal reflection effects; the allowance for the Fresnel reflection is necessary for quantitative agreement with measurements. DPDW parameters are calculated for the first time for a biomodel with foreign objects, and optical parameters of the biomodel, as well as spatial parameters of the problem, are estimated for values at which one can visualize an foreign object embedded into the medium.

Development of a multi-frequency diffuse photon density wave device for the characterization of tissue damage at multiple depths

The ability to determine the depth and degree of cutaneous and subcutaneous tissue damage is critical for medical applications such as burns and pressure ulcers. The Diffuse Photon Density Wave (DPDW) methodology at near infrared wavelengths can be used to non-invasively measure the optical absorption and reduced scattering coefficients of tissue at depths of several millimeters. A multi-frequency DPDW system with one light source and one detector was constructed so that light is focused onto the tissue surface using an optical fiber and lens mounted to a digitally-controlled actuator which changes the distance between light source and detector. A variable RF generator enables the modulation frequency to be selected between 50 to 400MHz. The ability to digitally control both source-detector separation distance and modulation frequency allows for virtually unlimited number of data points, enabling precise selection of the volume and depth of tissue that will be characterized. Suspensions of Intralipid and india ink with known absorption and reduced scattering coefficients were used as optical phantoms to assess device accuracy. Solid silicon phantoms were formulated for stability testing. Standard deviations for amplitude and phase shift readings were found to be 0.9% and 0.2 degrees respectively, over a one hour period. The ability of the system to quantify tissue damage in vivo at multiple depths was tested in a porcine burn model.