Diffuse light propagation in a turbid medium with varying refractive index: Monte Carlo modeling in a spherically symmetrical geometry

We report on the development of Monte Carlo software that can model media with spatially varying scattering coefficient, absorption, and refractive index. The varying refractive index is implemented by calculating curved photon paths in the medium. The results of the numerical simulations are compared with analytical solutions obtained using the diffusion approximation. The model under investigation is a scattering medium that contains a spherically symmetrical inclusion (inhomogeneity) created by variation in optical properties and having no sharp boundaries. The following steady-state cases are considered: (a) a nonabsorbing medium with a spherically symmetrical varying refractive index, (b) an inclusion with varying absorption and scattering coefficients and constant refractive index, and (c) an inclusion with varying absorption, scattering, and refractive index. In the latter case it is shown that the interplay between the absorption coefficient and the refractive index may create the effect of a hidden inclusion.