Tiziano Binzoni

The exact homogenized probability density function, for a photon making a step of length s has been analytically derived for a binary (isotropic-Poisson) statistical mixture with unmatched refractive indexes. The companions, exact, homogenized probability density functions for a photon to change direction (“scatter”), with polar ϑ and azimuthal φ angles, and the homogenized albedo, have also been obtained analytically. These functions also apply to negative refractive indexes and can reduce the number of Monte Carlo simulations needed for photon propagation in complex binary (isotropic-Poisson) statistical mixtures from hundreds to just one, for an equivalent homogeneous medium. Note, that this is not an approximate approach, but a mathematically equivalent and exact result. Additionally, some tutorial examples of homogenized Monte Carlo simulations are also given. Published by the American Physical Society 2024

Monte Carlo (MC) simulations allowing to describe photons propagation in statistical mixtures represent an interest that goes way beyond the domain of optics, and can cover, e.g., nuclear reactor physics, image analysis or life science just to name a few. MC simulations are considered a “gold standard” because they give exact solutions (in the statistical sense), however, in the case of statistical mixtures their implementation is often extremely complex. For this reason, the aim of the present contribution is to propose a new approach that should allow us in the future to simplify the MC approach. This is done through an explanatory example, i.e.; by deriving the ‘exact’ analytical expression for the probability density function of photons’ random steps (single step function, SSF) propagating in a medium represented as a binary (isotropic-Poisson) statistical mixture. The use of the SSF reduces the problem to an ‘equivalent’ homogeneous medium behaving exactly as the original binary statistical mixture. This will reduce hundreds MC simulations, allowing to obtain one set of wanted parameters, to only one equivalent simple MC simulation. To the best of our knowledge the analytically ‘exact’ SSF for a binary (isotropic-Poisson) statistical mixture has never been derived before.

Monte Carlo (MC) simulations allowing to describe photons propagation in statistical mixtures represent an interest that goes way beyond the domain of optics, and can cover, e.g., nuclear reactor physics, image analysis or life science just to name a few. MC simulations are considered a ``gold standard'' because they give exact solutions (in the statistical sense), however, in the case of statistical mixtures they are enormously time consuming and their implementation is often extremely complex. For this reason, the aim of the present contribution is to propose a new approach that should allow us in the future to simplify the MC approach. This is done through an explanatory example, i.e.; by deriving the `exact' analytical expression for the probability density function of photons' random steps (single step function, SSF) propagating in a medium represented as a binary (isotropic-Poisson) statistical mixture. The use of the SSF reduces the problem to an `equivalent' homogeneous medium behaving exactly as the original binary statistical mixture. This will reduce hundreds time-consuming MC simulations to only one equivalent simple MC simulation. To the best of our knowledge the analytically `exact' SSF for a binary (isotropic-Poisson) statistical mixture has never been derived before.

Significance: Code verification is an unavoidable step prior to using a Monte Carlo (MC) code. Indeed, in biomedical optics, a widespread verification procedure for MC codes is still missing. Analytical benchmarks that can be easily used for the verification of different MC routines offer an important resource. Aim: We aim to provide a two-step verification procedure for MC codes enabling the two main tasks of an MC simulator: (1) the generation of photons’ trajectories and (2) the intersections of trajectories with boundaries separating the regions with different optical properties. The proposed method is purely based on elementary analytical benchmarks, therefore, the correctness of an MC code can be assessed with a one-sample t-test. Approach: The two-step verification is based on the following two analytical benchmarks: (1) the exact analytical formulas for the statistical moments of the spatial coordinates where the scattering events occur in an infinite medium and (2) the exact invariant solutions of the radiative transfer equation for radiance, fluence rate, and mean path length in media subjected to a Lambertian illumination. Results: We carried out a wide set of comparisons between MC results and the two analytical benchmarks for a wide range of optical properties (from non-scattering to highly scattering media, with different types of scattering functions) in an infinite non-absorbing medium (step 1) and in a non-absorbing slab (step 2). The deviations between MC results and exact analytical values are usually within two standard errors (i.e., t-tests not rejected at a 5% level of significance). The comparisons show that the accuracy of the verification increases with the number of simulated trajectories so that, in principle, an arbitrary accuracy can be obtained. Conclusions: Given the simplicity of the verification method proposed, we envision that it can be widely used in the field of biomedical optics.

Anomalous radiative transfer (ART) theory is a generalization of classical radiative transfer theory. The present tutorial wants to show how Monte Carlo (MC) codes describing photons transport in anomalous media can be implemented. It is shown that the heart of the method consists in suitably describing, in a "non-classical" manner, photons steps starting from fixed light sources or from boundaries separating regions of the medium with different optical properties. To give a better intuition of the importance of these particular photon step lengths, it is also shown numerically that the described approach is essential to preserve the invariance property for light propagation. An interesting byproduct of the MC method for ART, is that it allow us to simplify the structure of "classical" MC codes, utilized e.g. in biomedical optics.

In the present contribution, three means of measuring the geometrical and topological complexity of photons’ paths in random media are proposed. This is realized by investigating the behavior of the average crossing number, the mean writhe, and the minimal crossing number of photons’ paths generated by Monte Carlo (MC) simulations, for different sets of optical parameters. It is observed that the complexity of the photons’ paths increases for increasing light source/detector spacing, and that highly “knotted” paths are formed. Due to the particular rules utilized to generate the MC photons’ paths, the present results may have an interest not only for the biomedical optics community, but also from a pure mathematical point of view.

In the present contribution three means of measuring the geometrical and topological complexity of photons' paths in random media are proposed. This is realized by investigating the behavior of the average crossing number, the mean writhe, and the minimal crossing number of photons' paths generated by Monte Carlo (MC) simulations, for different sets of optical parameters. It is observed that the complexity of the photons' paths increases for increasing light source/detector spacing, and that highly "knotted" paths are formed. Due to the particular rules utilized to generate the MC photons' paths, the present results may have an interest not only for the biomedical optics community, but also from a pure mathematical point of view.

A generalized time-independent correlation transport equation (GCTE) is proposed for the field autocorrelation function. The GCTE generalizes various models for anomalous transport of photons and takes into account the possible presence of a static background. In a tutorial example, the GCTE is solved for a homogeneous semi-infinite medium in reflectance configuration through Monte Carlo simulations. The chosen anomalous photon transport model also includes the classic and the “generalized” Lambert–Beer’s law (depending on the choice of parameters). A numerical algorithm allowing generation of the related anomalous random photon steps is also given. The clear influence of anomalous transport on the field autocorrelation function is shown and discussed for the proposed specific examples by comparing the general results with the classical case (Lambert–Beer’s law).

In the present contribution, the theoretical relationship existing between the blood flow index measured by diffuse correlation spectroscopy and single-photon laser-Doppler flowmetry (SP-LDF) is investigated. A specific mathematical description that accounts for the properties of single-photon detectors for SP-LDP was developed. Static background has also been considered and, to the best of our knowledge, this has never been included before in SP-LDF analytical theories. The comparisons were realized for two SP-LDF implementations: for “classical” and “fast” algorithms. “Classical” SP-LDF is not sensitive to small changes on the number of detected speckles and coherence length of the laser, usually described by a unique parameter “beta.” This is a strong point when assessing blood flow index, e.g., in humans, where “beta” is particularly difficult to be determined in real time. The proposed theory may be utilized, e.g., to investigate other SP-LDF setups and optical/physiological parameter ranges or, generally, to optimize real SP-LDF instrumentation.