Optimization of the Monte Carlo code for modeling of photon migration in tissue

The Monte Carlo method is frequently used to simulate light transport in turbid media because of its simplicity and flexibility, allowing to analyze complicated geometrical structures. Monte Carlo simulations are, however, time consuming because of the necessity to track the paths of individual photons. The time consuming computation is mainly associated with the calculation of the logarithmic and trigonometric functions as well as the generation of pseudo-random numbers. In this paper, the Monte Carlo algorithm was developed and optimized, by approximation of the logarithmic and trigonometric functions. The approximations were based on polynomial and rational functions, and the errors of these approximations are less than 1% of the values of the original functions. The proposed algorithm was verified by simulations of the time-resolved reflectance at several source-detector separations. The results of the calculation using the approximated algorithm were compared with those of the Monte Carlo simulations obtained with an exact computation of the logarithm and trigonometric functions as well as with the solution of the diffusion equation. The errors of the moments of the simulated distributions of times of flight of photons (total number of photons, mean time of flight and variance) are less than 2% for a range of optical properties, typical of living tissues. The proposed approximated algorithm allows to speed up the Monte Carlo simulations by a factor of 4. The developed code can be used on parallel machines, allowing for further acceleration.

[1]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[2]  B. Wilson,et al.  Monte Carlo modeling of light propagation in highly scattering tissues. I. Model predictions and comparison with diffusion theory , 1989, IEEE Transactions on Biomedical Engineering.

[3]  S Andersson-Engels,et al.  Real-time method for fitting time-resolved reflectance and transmittance measurements with a monte carlo model. , 1998, Applied optics.

[4]  Lihong V. Wang,et al.  Monte Carlo Modeling of Light Transport in Multi-layered Tissues in Standard C , 1992 .

[5]  S. Prahl.  Light transport in tissue , 1988 .

[6]  J. Spanier,et al.  Perturbation Monte Carlo methods to solve inverse photon migration problems in heterogeneous tissues. , 2001, Optics letters.

[7]  S. Arridge,et al.  Optical imaging in medicine: II. Modelling and reconstruction , 1997, Physics in medicine and biology.

[8]  M. H. Koelink,et al.  Condensed Monte Carlo simulations for the description of light transport. , 1993, Applied optics.

[9]  H. Kahn,et al.  Random sampling (Monte Carlo) techniques in neutron attenuation problems--I. , 1950, Nucleonics.

[10]  H. Kahn.  Random sampling (Monte Carlo) techniques in neutron attenuation problems--II. , 1950, Nucleonics.

[11]  Gilbert Jarry,et al.  A simulation method for the study of laser transillumination of biological tissues , 2006, Annals of Biomedical Engineering.

[12]  Shu Tezuka,et al.  On the lattice structure of the add-with-carry and subtract-with-borrow random number generators , 1993, TOMC.

[13]  M S Patterson,et al.  Determination of the optical properties of turbid media from a single Monte Carlo simulation , 1996, Physics in medicine and biology.

[14]  D. Delpy,et al.  Optical Imaging in Medicine , 1998, CLEO/Europe Conference on Lasers and Electro-Optics.

[15]  A. Villringer,et al.  Time-resolved multidistance near-infrared spectroscopy of the adult head: intracerebral and extracerebral absorption changes from moments of distribution of times of flight of photons. , 2004, Applied optics.

[16]  A. Villringer,et al.  Determining changes in NIR absorption using a layered model of the human head , 2001, Physics in medicine and biology.

[17]  Joshua J. Bloch.  Effective Java : programming language guide , 2001 .

[18]  David Thomas,et al.  The Art in Computer Programming , 2001 .

[19]  R. Alcouffe,et al.  Comparison of finite-difference transport and diffusion calculations for photon migration in homogeneous and heterogeneous tissues. , 1998, Physics in medicine and biology.

[20]  Y Phaneendra Kumar,et al.  Reconstruction of optical properties of low-scattering tissue using derivative estimated through perturbation Monte-Carlo method. , 2004, Journal of biomedical optics.

[21]  B. Wilson,et al.  Monte Carlo modeling of light propagation in highly scattering tissues. II. Comparison with measurements in phantoms , 1989, IEEE Transactions on Biomedical Engineering.

[22]  Antonio Pifferi,et al.  Accelerated Monte Carlo models to simulate fluorescence spectra from layered tissues. , 2003, Journal of the Optical Society of America. A, Optics, image science, and vision.

[23]  D. Boas,et al.  Three dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head. , 2002, Optics express.

[24]  L Wang,et al.  MCML--Monte Carlo modeling of light transport in multi-layered tissues. , 1995, Computer methods and programs in biomedicine.

[25]  F. James.  A Review of Pseudorandom Number Generators , 1990 .

[26]  A. Witt,et al.  Multiple scattering in reflection nebulae. I - A Monte Carlo approach. II - Uniform plane-parallel nebulae with foreground stars. III - Nebulae with embedded illuminating stars , 1977 .

[27]  M. Schweiger,et al.  Theoretical and experimental investigation of near-infrared light propagation in a model of the adult head. , 1997, Applied optics.

[28]  George Marsaglia,et al.  Toward a universal random number generator , 1987 .

[29]  W Steenbergen,et al.  Doppler Monte Carlo simulations of light scattering in tissue to support laser-Doppler perfusion measurements. , 1999, Technology and health care : official journal of the European Society for Engineering and Medicine.

[30]  Heidrun Wabnitz,et al.  Evaluation of optical properties of highly scattering media by moments of distributions of times of flight of photons. , 2003, Applied optics.

[31]  I. Vattulainen,et al.  A comparative study of some pseudorandom number generators , 1995 .

[32]  Annamaria Kisslinger,et al.  Multiple processor version of a Monte Carlo code for photon transport in turbid media , 2000 .

[33]  D. Delpy,et al.  Near-infrared light propagation in an adult head model. I. Modeling of low-level scattering in the cerebrospinal fluid layer. , 2003, Applied optics.

[34]  H.J.C.M. Sterenborg,et al.  Skin optics , 1989, IEEE Transactions on Biomedical Engineering.

[35]  Multiple scattering in reflection nebulae. IV. The multiplicity of scattering. , 1977 .