Monte Carlo simulation of light fluence in tissue in a cylindrical diffusing fibre geometry

The propagation of light emitted by a linear light diffuser in a cylindrical hollow organ was investigated by means of the Monte Carlo (MC) method. The height and radius of the cavity, scattering (mu(s)) (or reduced scattering, mu(s)') and absorption (mu(a)) coefficients, anisotropy (g), and refractive indices of the media involved (nl, nr) are required as input data by the MC code, as are characteristics of the light diffuser (length, delivered power and emission profile). Results of our MC model were tested by measuring the light fluence rate in a tissue-simulating phantom (mu(a) = 0.5 cm(-1), mu(s) = 23 cm(-1) and g = 0.75) irradiated at 633 nm with a cylindrical diffuser. Since geometric and optical parameters determine the behaviour of light propagation in tissue, MC simulations with different sets of input parameters were carried out to provide qualitative as well as quantitative data useful in planning photodynamic therapy. Data are reported on light penetration and fluence rate build-up at mu(a) and mu(s)' values ranging between 0.1 and 5 cm(-1) and 2.5 and 50 cm(-1), respectively. Furthermore, results suggest that a shift and spread could occur in the isofluence curves along the symmetry axis, which depend on the diameter of the treated lumen as well as on the emission profile of the light diffuser. Using our data it is possible to estimate how inaccuracy in knowledge of the optical coefficients can affect (i.e. usually by increasing) the light dose scheduled at a certain depth into tissue.