OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 20, Iss. 15 — Jul. 16, 2012
  • pp: 16618–16630

GPU accelerated electric field Monte Carlo simulation of light propagation in turbid media using a finite-size beam model

Yaru Wang, Pengcheng Li, Chao Jiang, Jia Wang, and Qingming Luo  »View Author Affiliations

Optics Express, Vol. 20, Issue 15, pp. 16618-16630 (2012)

View Full Text Article

Enhanced HTML    Acrobat PDF (5068 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



An electric field Monte Carlo (EMC) simulation directly traces the complex electric field vectors in multiple scattering and estimates the electric field in a preferred direction. The full vectorial nature of EMC makes it a powerful and flexible tool to simulate the coherence and polarization phenomena of light. As a numerical method, EMC needs to launch a large number of photons to achieve an accurate result, making it time-consuming. Previously, EMC did not account for the beam size. Because of the stochastic character of the instantaneous electric field in the simulation, the convolution method alone is unsuitable for the Monte Carlo simulation of photon energy for a beam with a finite size. It is necessary to launch photons from all possible locations to simulate a finite-size beam, which results in a significant increase in the computational burden. In order to accelerate the simulation, a parallel implementation of the electric field Monte Carlo simulation based on the compute unified device architecture (CUDA) running on a graphics processing unit (GPU) is presented in this paper. Our program, which is optimized for Fermi architecture, is able to simulate the coherence phenomenon of a finite-size beam normally incident on turbid media. A maximum speedup of over 370x is achieved with a GTX480 GPU, compared with that obtained using an Intel i3-2120 CPU.

© 2012 OSA

OCIS Codes
(030.6140) Coherence and statistical optics : Speckle
(170.5280) Medical optics and biotechnology : Photon migration
(290.1350) Scattering : Backscattering
(290.4210) Scattering : Multiple scattering
(290.7050) Scattering : Turbid media

ToC Category:

Original Manuscript: April 19, 2012
Revised Manuscript: June 29, 2012
Manuscript Accepted: July 1, 2012
Published: July 9, 2012

Virtual Issues
Vol. 7, Iss. 9 Virtual Journal for Biomedical Optics

Yaru Wang, Pengcheng Li, Chao Jiang, Jia Wang, and Qingming Luo, "GPU accelerated electric field Monte Carlo simulation of light propagation in turbid media using a finite-size beam model," Opt. Express 20, 16618-16630 (2012)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. A. Ishimaru, Wave Propagation and Scattering in Random Media (Academic, New York, 1999), Vol.1.
  2. J. W. Goodman, “Statistical properties of laser speckle patterns,” in Laser speckle and related phenomena, 2nd ed., J. C. Dainty, ed. (Springer-Verlag, New York, 1984), pp. 9–75.
  3. J. W. Goodman, Speckle Phenomena in Optics: Theory and Applications (Roberts and Company, Englewood, Colorado, 2007).
  4. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and et, “Optical coherence tomography,” Science254(5035), 1178–1181 (1991). [CrossRef] [PubMed]
  5. R. Drezek, A. Dunn, and R. Richards-Kortum, “Light scattering from cells: finite-difference time-domain simulations and goniometric measurements,” Appl. Opt.38(16), 3651–3661 (1999). [CrossRef] [PubMed]
  6. L. Wang, S. L. Jacques, and L. Zheng, “MCML–Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Meth. Prog. Biomed.47(2), 131–146 (1995). [CrossRef]
  7. B. C. Wilson and G. Adam, “A Monte Carlo model for the absorption and flux distributions of light in tissue,” Med. Phys.10(6), 824–830 (1983). [CrossRef] [PubMed]
  8. S. Bartel and A. H. Hielscher, “Monte Carlo simulations of the diffuse backscattering Mueller matrix for highly scattering media,” Appl. Opt.39(10), 1580–1588 (2000). [CrossRef] [PubMed]
  9. M. Xu, “Electric field Monte Carlo simulation of polarized light propagation in turbid media,” Opt. Express12(26), 6530–6539 (2004). [CrossRef] [PubMed]
  10. S. Moon, D. Kim, and E. Sim, “Monte Carlo study of coherent diffuse photon transport in a homogeneous turbid medium: a degree-of-coherence based approach,” Appl. Opt.47(3), 336–345 (2008). [CrossRef] [PubMed]
  11. J. Sawicki, N. Kastor, and M. Xu, “Electric field Monte Carlo simulation of coherent backscattering of polarized light by a turbid medium containing Mie scatterers,” Opt. Express16(8), 5728–5738 (2008). [CrossRef] [PubMed]
  12. M. Xu and R. R. Alfano, “Random walk of polarized light in turbid media,” Phys. Rev. Lett.95(21), 213901 (2005). [CrossRef] [PubMed]
  13. M. Xu and R. R. Alfano, “Circular polarization memory of light,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys.72(6), 065601 (2005). [CrossRef] [PubMed]
  14. K. G. Phillips, M. Xu, S. Gayen, and R. Alfano, “Time-resolved ring structure of circularly polarized beams backscattered from forward scattering media,” Opt. Express13(20), 7954–7969 (2005). [CrossRef] [PubMed]
  15. R. Liao, H. Zhu, Y. Huang, and J. Lv, “Monte Carlo modelling of OCT with finite-size-spot photon beam,” Chin. Opt. Lett.3, S346–S347 (2005).
  16. C. K. Hayakawa, V. Venugopalan, V. V. Krishnamachari, and E. O. Potma, “Amplitude and phase of tightly focused laser beams in turbid media,” Phys. Rev. Lett.103(4), 043903 (2009). [CrossRef] [PubMed]
  17. M. Šormaz and P. Jenny, “Contrast improvement by selecting ballistic-photons using polarization gating,” Opt. Express18(23), 23746–23755 (2010). [CrossRef] [PubMed]
  18. C. K. Hayakawa, E. O. Potma, and V. Venugopalan, “Electric field Monte Carlo simulations of focal field distributions produced by tightly focused laser beams in tissues,” Biomed. Opt. Express2(2), 278–290 (2011). [CrossRef] [PubMed]
  19. J. Von Neumann, “Various techniques used in connection with random digits,” J. Res. Natl. Bur. Stand.5, 36–38 (1951).
  20. E. Alerstam, W. C. Y. Lo, T. D. Han, J. Rose, S. Andersson-Engels, and L. Lilge, “Next-generation acceleration and code optimization for light transport in turbid media using GPUs,” Biomed. Opt. Express1(2), 658–675 (2010). [CrossRef] [PubMed]
  21. G. Marsaglia and A. Zaman, “A new class of random number generators,” Ann. Appl. Probab.1(3), 462–480 (1991). [CrossRef]
  22. Nvidia Corporation, “CUDA programming guide 4.2,” (2012).
  23. G. Sendra, H. Rabal, M. Trivi, and R. Arizaga, “Numerical model for simulation of dynamic speckle reference patterns,” Opt. Commun.282(18), 3693–3700 (2009). [CrossRef]
  24. D. D. Duncan and S. J. Kirkpatrick, “The copula: a tool for simulating speckle dynamics,” J. Opt. Soc. Am. A25(1), 231–237 (2008). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited