Calculation of photon path changes due to scatter in Monte Carlo simulations

Abstract

This is the accepted version of the paper.This version of the publication may differ from the final published version. Abstract-Computation using Monte Carlo simulations is widely used for modelling the light-tissue interaction. Despite this, many of the methods used for building such simulations are poorly described in the literature. In particular, a scheme for translating the scatter angles produced from a phase function into updated photon direction vectors is not explicitly reported. To address this, a method for calculating the change in photon direction following a scattering event is described, thus illuminating one of the fundamental 'building blocks' for researchers developing their own Monte Carlo models. The equations derived in this paper may be readily incorporated into applicable Monte Carlo program code. Permanent

Investigating optical path and differential pathlength factor in reflectance photoplethysmography for the assessment of perfusion

“…Detailed calculations of the scattering angles and step sizes are described elsewhere. 29,51 Each of the launched photon packets was introduced with an initial statistical weight w ¼ 1. After each step size, when the photon packet reached an interaction site, a fraction of its weight (w: μ a μ a þμ s ) was modeled to be absorbed within the tissue layer.…”

Differential pathlength factor estimation for brain-like tissue from a single-layer Monte Carlo model

“…It is noteworthy that the Monte Carlo method relies on random sampling of variables from probability distributions [11]. The occurrence of the scattering and absorption events is decided by comparing a random number with the probability of absorption [11,12]. If the photon is not absorbed, scattering takes place.…”

“…If the photon is not absorbed, scattering takes place. The direction of scattering is calculated depending on the Henyey -Greenstein phase function [12]. The position and direction of the scattered photon are updated [12].…”

“…The direction of scattering is calculated depending on the Henyey -Greenstein phase function [12]. The position and direction of the scattered photon are updated [12]. If the scattered photon leaves the tissue medium, a check is made to see if the photon has fulfilled the detection criteria, i.e.…”

Monte Carlo Simulation of the Effect of Human Skin Melanin in Light-Tissue Interactions

“…Moreover, the simulation was limited to the reliable retrieval of data of the anisotropic factors of the epidermis at the selected wavelengths. Instead, the anisotropic factors of skin at similar wavelengths [12,13] were used to estimate g.…”