OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 17, Iss. 19 — Sep. 14, 2009
  • pp: 16590–16602
« Show journal navigation

Monte Carlo simulation of polarized photon scattering in anisotropic media

Tianliang Yun, Nan Zeng, Wei Li, Dongzhi Li, Xiaoyu Jiang, and Hui Ma  »View Author Affiliations


Optics Express, Vol. 17, Issue 19, pp. 16590-16602 (2009)
http://dx.doi.org/10.1364/OE.17.016590


View Full Text Article

Acrobat PDF (257 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We present a Monte Carlo simulation program for the propagation of polarized photons in an anisotropic scattering model consisting of poly-dispersed spherical and infinite long cylindrical scatterers. The cylinders are aligned following a Gaussian distribution. Densities and sizes of the spherical and cylindrical scatterers, as well as the orientation of the cylinders are variables for the simulation of different anisotropic media. The good agreement between the simulation and experimental results of polarization imaging confirms the validity of the polarization-dependent Monte Carlo simulation program.

© 2009 OSA

1. Introduction

Monte Carlo simulation has been an effective method for the investigation of light propagation in turbid media, such as biological tissues. The source code of a Monte Carlo program for simulating the behavior of photon transport in multilayered turbid media was made public [1

1. L. H. Wang, S. L. Jacques, and L.-Q. Zheng, “MCML--Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed. 47(2), 131–146 (1995). [PubMed]

]. In the program, the photons are unpolarized and the scattering phase function follows the Henyey-Greenstein function for isotropic media.

In the past few years there are increasingly interests in the light-scattering properties of tissues for therapeutic and diagnostic purposes. These applications motivate efforts in improving numerical simulations for real biological samples, which often contain complex tissue microstructure. Some groups have investigated light scattering from nonspherical particles, such as disc [10

10. J. He, A. Karlsson, J. Swartling, and S. Andersson-Engels, “Light scattering by multiple red blood cells,” J. Opt. Soc. Am. A 21, 1953–1961 (2004).

] and ellipsoid [11

11. J. D. Keener, K. J. Chalut, J. W. Pyhtila, and A. Wax, “Application of Mie theory to determine the structure of spheroidal scatterers in biological materials,” Opt. Lett. 32(10), 1326–1328 (2007). [PubMed]

]. Many tissues contain microstructure of ordered elongated subunits, such as the myofibrils in muscles or the collagen fibers in skin, tendons, or ligaments. They are anisotropic and cannot be described by an isotropic scattering model, such as randomly distributed spherical or non spherical scatterers. Recently, Monte Carlo simulations of unpolarized photon propagation in anisotropic media, such as dentin [12

12. A. Kienle and R. Hibst, “Light guiding in biological tissue due to scattering,” Phys. Rev. Lett. 97(1), 018104 (2006). [PubMed]

] and softwood [13

13. A. Kienle, C. D’Andrea, F. Foschum, P. Taroni, and A. Pifferi, “Light propagation in dry and wet softwood,” Opt. Express 16(13), 9895–9906 (2008). [PubMed]

], have been reported. The anisotropic media were approximated to a mixture of aligned cylinders of infinite length and solid microspheres [14

14. A. Kienle, F. K. Forster, and R. Hibst, “Anisotropy of light propagation in biological tissue,” Opt. Lett. 29(22), 2617–2619 (2004). [PubMed]

].

In this paper, we use a similar sphere-cylinder model for anisotropic scattering media. This model consists of both spherical and infinite long cylindrical scattering particles. A polarization-dependent Monte Carlo program is developed to simulate the behaviors of polarized photons as they propagate through the medium. Densities and sizes of the spherical and cylindrical scatterers, as well as the orientation of the cylinders are variables in the simulation of different anisotropic media. The good agreement between the simulation and the experimental results from polarization imaging confirms the validity of this anisotropic scattering model and of the polarization-dependent Monte Carlo program.

2. Methods: Monte Carlo simulation

2.1 Polarized photon scattering at an infinitely long cylinder

Using the continuity conditions on the boundary of a cylinder, the scalar wave equation 2ψ+k2ψ=0 has been solved analytically in cylindrical polar coordinate for light scattering by an infinitely long cylinder. The scattered fields in matrix form and relevant parameters are given as [15

15. C. F. Bohren, and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley and Sons, New York, 1983).

]:
(EsEs)=ei3π/42πkrsinζeik(rsinζzcosζ)(T1T3T4T2)(EiEi)
(1)
{T1=bnIeinΘ=b0I+2n=1bnIcos(nΘ)T2=anIIeinΘ=a0II+2n=1anIIcos(nΘ)T3=anIeinΘ=2in=1anIsin(nΘ)T4=bnIIeinΘ=2in=1bnIIsin(nΘ)
(2)
{anI=CnVnBnDnWnVn+iDn2,anII=AnVniCnDnWnVn+iDn2bnI=WnBn+iDnCnWnVn+iDn2,bnII=iCnWn+AnDnWnVn+iDn2
(3)
{An=iξ[ξJn(η)Jn(ξ)ηJn(η)Jn(ξ)]Bn=ξ[m2ξJn(η)Jn(ξ)ηJn(η)Jn(ξ)]Cn=ncosζηJn(η)Jn(ξ)(ξ2η21)Dn=ncosζηJn(η)Hn(1)(ξ)(ξ2η21)Wn=iξ[ηJn(η)Hn(1)(ξ)ξJn(η)Hn(1)(ξ)]Vn=ξ[m2ξJn(η)Hn(1)(ξ)ηJn(η)Hn(1)(ξ)]
(4)
{Θ=πθξ=xsinζη=xm2cos2ζx=ka
(5)
where Ei and Ei represent the incident electrical field, r and z are cylindrical polar coordinates, a is the radius of the cylinder, m is the refractive index of the cylinder relative to that of the surrounding medium, k is the wave number, ζ and θ are the incident angle and the angle between the incident and scattering planes as shown in Fig. 1
Fig. 1 Scheme of light scattering characteristics by an infinitely long cylinder
. Jn is the Bessel function and Hn(1)=Jn+iYn is the Hankel function. The Mueller matrix M(ζ,θ) can be derived from Eq. (1) according to the definition of Stokes vector:
M(ζ,θ)=2πkrsinζ(m11m12m13m14m21m22m23m24m31m32m33m34m41m42m43m44)
(6)
where
{m11=(|T1|2+|T2|2+|T3|2+|T4|2)/2m12=(|T1|2|T2|2+|T3|2|T4|2)/2m13=Re{T1T4+T2T3}m14=Im{T1T4T2T3}m21=(|T1|2|T2|2|T3|2+|T4|2)/2m22=(|T1|2+|T2|2|T3|2|T4|2)/2m23=Re{T1T4T2T3}m24=Im{T1T4+T2T3}m31=Re{T1T3+T2T4}m32=Re{T1T3T2T4}m33=Re{T1T2+T3T4}m34=Im{T1T2+T3T4}m41=Im{T1*T3+T2T4}m42=Im{T1*T3T2T4}m43=Im{T1*T2T3*T4}m44=Re{T1T2T3T4}
(7)
The matrix elements of M(ζ,θ)can be calculated using Eq. (2-7). In our simulation, the angle resolutions of ζ and θ of M(ζ,θ) are 1°.

2.1.1 The phase function

The spatial distribution of scattered photons for oblique incidence at an infinitely long cylinder follows a conical shape around the cylinder with the half angle ζ which is the same as the incident angle [see Fig. 1]. The scattering probability around the cone depends on the scattering direction defined by θ. When a scattering occurs, a new reference frame is defined by the incident direction of the photon and the orientation of the cylinder. The Stokes vector must be transformed so that the polarization is properly represented with respect to the new reference frame. Then the new Stokes vector of scattered photon at any scattering direction θ can be calculated as follow:
Ss(Is,Qs,Us,Vs)=M(ζ,θ)RS
(8)
where R is the rotational matrix transforming the initial polarization state S to the new incident reference frame, M(ζ,θ) is the Mueller matrix whose elements are calculated according to Eq. (2-7). We then obtain the scattering phase functions corresponding to different incident angle ζ and polarization state S by normalizing the Stokes vectorIs(θ). For unpolarized incoming photons, the scattering phase functions at different incident angles (Fig. 2(a)
Fig. 2 (a) Phase functions for unpolarized light with different incident angles (b) Phase functions for different polarizations of normally incident light
) are identical to the previously published data (Fig. 2(b) in Ref [16

16. A. Kienle, F. Forster, R. Diebolder, and R. Hibst, “Light propagation in dentin: influence of microstructure on anisotropy,” Phys. Med. Biol. 48(7–N), 14 (2003).

].) which was calculated using the unpolarized photon scattering model. Figure 2(b) shows phase functions corresponding to normal incidence but different polarizations. In both figures, we choose 2 µm as the diameter of cylindrical scatterers. The refractive indices of the cylinder and its surrounding are 1.33 and 1.52 respectively. The wavelength of the incident light is 633 nm.

2.1.2 The scattering coefficient

The scattering coefficient of the microspheres μs,sph is a constant once their size and refractive index are set. However, the scattering coefficient of cylinders, μs,cyl, varies with the angle ζ between the direction of incident photon and the cylinder [16

16. A. Kienle, F. Forster, R. Diebolder, and R. Hibst, “Light propagation in dentin: influence of microstructure on anisotropy,” Phys. Med. Biol. 48(7–N), 14 (2003).

]:
μs,cyl(ζ)=Qsca(ζ)dCA
(9)
where d and CA are the diameter and density of the cylinders. Qsca is the scattering efficiency of a single cylinder, which is a function of the angle ζ and the polarization state of the incident photon. In the program, we first determine ζ according to the spatial orientation of the cylinder and the propagation direction of the photon. If the electric field is parallel to the cylinder, efficiency Qsca,//(ζ) is calculated using the analytic solution of the Maxwell equation [15

15. C. F. Bohren, and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley and Sons, New York, 1983).

]. For arbitrary polarization states S=[I,Q,U,V], Qsca,S(ζ) is calculated as follow:

Qsca,S(ζ)=Ws,SWs,//Qsca,//(ζ)=θ=02πIs,S(θ)θ=02πIs,//(θ)Qsca,//(ζ)             =θ=02πM(ζ,θ)(I,Q,U,V)T(1)θ=02πM(ζ,θ)(1,1,0,0)T(1)Qsca,//(ζ)
(10)

2.2 The framework of program

To simulate anisotropic behavior of the media, the new program introduces infinitely long cylinders into the scattering media. The number of different cylinders as well as their sizes, concentrations and orientation functions are preset parameters in the input data file. New subroutines are added for calculations of cylindrical scatters, such as scattering probability, phase function, rotation of the Stokes vector reference plane etc.

A flow-chart of the program is shown in Fig. 3
Fig. 3 Flow chart of Monte Carlo program of polarized light scattering in medium consisting of spheres and cylinders
. For polarized Monte Carlo, the polarization states of photons are tracked step by step in addition to their spatial positions and direction of motions. The detailed descriptions of the program are as follow:

  • 1) The program first establishes the laboratory coordinate system: XYplane represents the surface of semi-infinite medium and Z is the depth direction. Two unit vectors v and w define the E field reference frame. v=[vx,vy,vz]=[0,1,0]. A photon of arbitrary polarization state is represented by a Stokes vector S=[I,Q,U,V] in the reference frame. A photon is launched at an initial position and direction determined by the geometry and the incident angle of the beam. The propagation direction of the photon is represented by a unit vector u, and w=u×v.
  • 2) Before being scattered, the photon moves to a new position which is calculated based on the extinction coefficient μt and a pseudo random number δ generated in the interval (0, 1).

    Δs=ln(δ)μt
    (11)

    where μt=μa+μs,sph+μs,cyl(ζ), μa is the absorption coefficient, μs,sph and μs,cyl are the scattering coefficient of spherical and cylindrical scatterers respectively. μs,cyl is a function of ζ as mentioned in section 2.1.2. Thus, the new position of the photon is updated as follow:

    {x=x+uxΔsy=y+uyΔsz=z+uzΔs
    (12)

    As in other Monte Carlo simulations [1

    1. L. H. Wang, S. L. Jacques, and L.-Q. Zheng, “MCML--Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed. 47(2), 131–146 (1995). [PubMed]

    ,4

    4. J. Ramella-Roman, S. Prahl, and S. L. Jacques, “Three Monte Carlo programs of polarized light transport into scattering media: part I,” Opt. Express 13(12), 4420–4438 (2005). [PubMed]

    ], we use a weight W to describe the absorption of light. By updating this weight after every step according to μs/(μs+μa), the weight of the photon Wwhich initially equals to 1 becomes

    Wafter=Wbefore(μsμs+μa)
    (13)

    where Wbefore is the weight before this step.

  • 3) The photon may hit either a spherical or cylindrical scatter,, so we need to make the choice before the next scattering event occurs. First, we determine ζ according to the spatial orientation distribution of the cylinders and the propagation direction of the photon, and calculate μs,cyl using Eq. (9). Then considering pcyl=μs,cyl/(μs,sph+μs,cyl) for the probability of scattering at cylinders and psph=1pcyl for spheres, we make a statistical choice on what type of scatterer the photon hits.
  • 4) The program carries out rotation of the Stokes vector reference frame and computes the phase function with the respective scattering matrix M. Details are significantly different for scattering at spheres and cylinders. The procedures for spherical scatterers are similar to a published work [4

    4. J. Ramella-Roman, S. Prahl, and S. L. Jacques, “Three Monte Carlo programs of polarized light transport into scattering media: part I,” Opt. Express 13(12), 4420–4438 (2005). [PubMed]

    ,5

    5. J. C. Ramella-Roman, S. A. Prahl, and S. L. Jacques, “Three Monte Carlo programs of polarized light transport into scattering media: part II,” Opt. Express 13(25), 10392–10405 (2005). [PubMed]

    ]. The procedure for scattering at a cylinder is as follow.
    • (a) Once the spatial orientation of cylinder (l) and the direction of photon propagation (u0) are known, the incident plane and the new reference frame (vi) are determined. vi=l×u0.
    • (b) The Stokes vector is rotated by an angle β for its reference plane to coincide with the incident plane, as shown in Fig. 4. β is the angle between vi and v0. The rotational matrix is

      R(β)=[10000cos(2β)sin(2β)00sin(2β)cos(2β)00001]
      (14)

Fig. 4 Infinite cylinder obliquely illuminated by a plane wave
  • (c) The new polarization state is calculated as follow,

    S(θ)=M(ζ,θ)R(β)S
    (15)

    The phase function and a random number determine the scattering direction us.

  • (d) Update the new reference frame of Stokes vector: vs=l×us.
    • 5) Repeat the above steps. If the weight W falls below a preset threshold, the photon is absorbed and the calculation terminates. If the photon moves beyond the boundary of the medium, the corresponding Stokes vector is calculated as

      Sout=WS
      (16)

      Before the photon is collected by the detector, its polarization state needs to be transformed so that its reference plane coincides with the detector. The first rotation about the direction of photon propagation u makes vz=0. The angle of this rotation ε is given by

      ε=0 when vz=0 and uz=0ε=tan1(vzwz) in all other cases
      (17)

      The second rotation by an angle φ about the Z-axis makes v parallel to the Y-axis. φ is calculated as

      φ=tan1(uyux) for backscattered photonφ=tan1(uyux) for transmitted photon
      (18)

      Thus the final Stokes vector is:

      Sfinal=R(φ)R(ε)Sout
      (19)

    • 6) Repeat above procedures for the next photon until the preset photon number is reached.

2.3 Model of the anisotropic medium

In this work, we approximate the anisotropic medium, such as the fibrous tissues, to a mixture of poly-dispersed microspheres and cylindrical scatterers of infinite length. The sizes and concentrations of the spheres and cylinders can be varied independently. The cylinders are aligned following a Gaussian distribution. The orientation and standard deviation Δη of the distribution function are also variable. Based on such a simple sphere-cylinder model, more complicated structure can be simulated.

3. Performance tests of the Monte Carlo program

3.1 Computation speed

In a typical test run using a personal computer with Intel Core2 Duo E6300 1.86 GHz processor and 2 G memory, it costs about 50 minutes to complete a simulation with the following parameters: thickness 2 cm, scattering coefficient 100 cm1, absorption coefficient 1 cm1, one million photons. With the same number of photons, the thickness and the total scattering coefficient of the medium significantly influence the computation time. The average scattering numbers of photons increase with the thickness or the total scattering coefficient. For polarization imaging simulation, the diffused photons quickly loose their polarization as they scatter around in the medium and may have little contribution in polarization sensitive detection. Thus, one may cut down the computation time by setting up a threshold for the maximum number of scattering a photon may experience.

3.2 Comparison with simulation of unpolarized light scattering in anisotropic medium

To test the validity of the program, we simulate polarized photon back scattering in an anisotropic scattering medium and compared the results with a previously published simulation for unpolarized photon back scattering using the same medium [14

14. A. Kienle, F. K. Forster, and R. Hibst, “Anisotropy of light propagation in biological tissue,” Opt. Lett. 29(22), 2617–2619 (2004). [PubMed]

]. We use linearly polarized photons with random polarization angles to simulate the unpolarized situation. The scattering medium is semi-infinite and contains infinitely long cylindrical scatters only. The cylinders are aligned in the plane and around the x axis following a Gaussian distribution with standard deviation Δη=10°. The diameter of the cylinders is 2 μm and the absorption coefficient is μa=0.1 cm −1. The scattering coefficient is chosen as 1420 cm −1 which was not given explicitly in Ref [14

14. A. Kienle, F. K. Forster, and R. Hibst, “Anisotropy of light propagation in biological tissue,” Opt. Lett. 29(22), 2617–2619 (2004). [PubMed]

].

In the simulation, linearly polarized photons with random polarization angle incident perpendicularly onto the surface of the medium (the X-Y plane) at x=y=0. Stokes vectors of the backscattered photons are calculated to obtain the backward spatial distribution of S=[I,Q,U,V]. The dimensionless spatially resolved reflectance R(x,y) is then calculated by normalizing I(x,y) to its maximum. The newly simulated result is shown in Fig. 5
Fig. 5 R(x, y) pattern obtained from Monte Carlo simulation of a semi-infinite turbid medium containing cylinders along the axis with Δη = 10°
using the same scale and format as Fig. 2(a) in Ref [14

14. A. Kienle, F. K. Forster, and R. Hibst, “Anisotropy of light propagation in biological tissue,” Opt. Lett. 29(22), 2617–2619 (2004). [PubMed]

]. The plot scales with scattering coefficient. For 1420 cm −1, the two figures look similar.

3.3 Comparison with polarization imaging experiments

More tests are carried out to compare the simulated results of reflective polarization imaging against those from experiments.

3.3.1 Experiment setup

The schematic diagram of the experimental setup is shown in Fig. 6
Fig. 6 Scheme of experimental setup
[17

17. X. Y. Jiang, N. Zeng, Y. H. He, and H. Ma, “Investigation of Linear Polarization Difference Imaging Based on Rotation of Incident and Backscattered Polarization Angles,” Pro. Biochem. Biophys. 34, 659 (2007).

]. The 650 nm light from a 1 W LED with proper attenuation is collimated by a lens and passes through a linear polarizer P1 and then illuminates the sample at an angle of 20° to the normal direction of its surface. Backscattered photons from the sample first pass though a linear polarizer P2, then are recorded by a CCD.

3.3.2 Imaging and data analysis

3.3.3 Result of the experiment

In the experiment, the sample is a three-layers medium which consists of a slab of well aligned silk fibers submerged in microsphere solution. The first and third layers are 1 cm thick solutions of 1.5 μm diameter polystyrene microsphere in water (International Laboratory, USA). The refractive index of the microsphere and water are 1.59 and 1.33 respectively. The scattering coefficients are 19.5 cm −1, 13 cm −1, 6.5 cm −1 and 3.25 cm −1 determined by the concentration of the microspheres. The second layer contains only well aligned silk fibers (provided by Guangxi Institute of Supervision and Testing on Product Quality). The thickness of the silk fiber slab is 0.1 cm and the refractive index is 1.56. The diameter of the cylindrical scatters is taken as 1.5 μm, since both literature and electron microscope images show that silk fibers contains micro fibers of such size..

We simulate the experimental results using both the sphere-only and sphere-cylinder models. The sphere-only model cannot give reasonable good match between the experimental and simulated results.

For the sphere-cylinder model, we first chose the experimental data with fixed microspheres scattering coefficient 3.3 cm −1, simulate the LDP curve using different parameters for the silk fiber slab until a good match between simulated and experimental results is reached. Then, we simulate the experimental results corresponding to different microsphere scattering coefficient. As shown in Fig. 7
Fig. 7 Dependence of Linear Differential Polarization to incident polarization angle
, the simulations of different samples using the sphere-cylinder model agree well with the experiments. The scattering coefficient of silks fiber slab is 70 cm −1. The standard deviation of the Gaussian distribution for spatial alignment of the fibers is 5°.

4. Conclusion

In summary, we present a Monte Carlo program that simulates the propagation of polarized light in anisotropic media. The medium for the simulation is composed of spherical and infinitely long cylindrical scatterers. Its anisotropic characters can be adjusted by changing the orientation distribution of the cylindrical scatterers, and the composition and the diameter of the two types of scatterers. The program is capable of simulating the behavior of polarized light propagation in anisotropic media, obtaining the spatial distribution and polarization states of the scattered photons. To validate the program, we demonstrated that a simulation using photons of randomly linear polarization agrees well with a previously published simulation using unpolarized photons. We also carried on more tests using samples which contain spheres and cylinders. The good agreement between the simulated and experimental results for such sphere-cylinder samples presents further evidence to the validity of the program.

Acknowledgement

This work has been supported by National Natural Science Foundation of China (grants 60578003 and 60778044), and Ministry of Science and Technology (grant 2006CB70570).

References and links

1.

L. H. Wang, S. L. Jacques, and L.-Q. Zheng, “MCML--Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed. 47(2), 131–146 (1995). [PubMed]

2.

G. W. Kattawar and G. N. Plass, “Radiance and polarization of multiple scattered light from haze and clouds,” Appl. Opt. 7(8), 1519–1527 (1968). [PubMed]

3.

P. Bruscaglioni, G. Zaccanti, and Q. Wei, “Transmission of a pulsed polarized light beam through thick turbid media: numerical results,” Appl. Opt. 32(30), 6142–6150 (1993). [PubMed]

4.

J. Ramella-Roman, S. Prahl, and S. L. Jacques, “Three Monte Carlo programs of polarized light transport into scattering media: part I,” Opt. Express 13(12), 4420–4438 (2005). [PubMed]

5.

J. C. Ramella-Roman, S. A. Prahl, and S. L. Jacques, “Three Monte Carlo programs of polarized light transport into scattering media: part II,” Opt. Express 13(25), 10392–10405 (2005). [PubMed]

6.

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).

7.

B. D. Cameron, M. J. Rakovic, M. Mehrübeoglu, G. W. Kattawar, S. Rastegar, L. V. Wang, and G. L. Coté, “Measurement and calculation of the two-dimensional backscattering Mueller matrix of a turbid medium,” Opt. Lett. 23(7), 485–487 (1998).

8.

M. J. Rakovi, G. W. Kattawar, M. B. Mehru Beo Lu, B. D. Cameron, L. V. Wang, S. Rastegar, and G. L. Coté, “Light backscattering polarization patterns from turbid media: theory and experiment,” Appl. Opt. 38(15), 3399–3408 (1999).

9.

D. Côté and I. Vitkin, “Robust concentration determination of optically active molecules in turbid media with validated three-dimensional polarization sensitive Monte Carlo calculations,” Opt. Express 13(1), 148–163 (2005). [PubMed]

10.

J. He, A. Karlsson, J. Swartling, and S. Andersson-Engels, “Light scattering by multiple red blood cells,” J. Opt. Soc. Am. A 21, 1953–1961 (2004).

11.

J. D. Keener, K. J. Chalut, J. W. Pyhtila, and A. Wax, “Application of Mie theory to determine the structure of spheroidal scatterers in biological materials,” Opt. Lett. 32(10), 1326–1328 (2007). [PubMed]

12.

A. Kienle and R. Hibst, “Light guiding in biological tissue due to scattering,” Phys. Rev. Lett. 97(1), 018104 (2006). [PubMed]

13.

A. Kienle, C. D’Andrea, F. Foschum, P. Taroni, and A. Pifferi, “Light propagation in dry and wet softwood,” Opt. Express 16(13), 9895–9906 (2008). [PubMed]

14.

A. Kienle, F. K. Forster, and R. Hibst, “Anisotropy of light propagation in biological tissue,” Opt. Lett. 29(22), 2617–2619 (2004). [PubMed]

15.

C. F. Bohren, and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley and Sons, New York, 1983).

16.

A. Kienle, F. Forster, R. Diebolder, and R. Hibst, “Light propagation in dentin: influence of microstructure on anisotropy,” Phys. Med. Biol. 48(7–N), 14 (2003).

17.

X. Y. Jiang, N. Zeng, Y. H. He, and H. Ma, “Investigation of Linear Polarization Difference Imaging Based on Rotation of Incident and Backscattered Polarization Angles,” Pro. Biochem. Biophys. 34, 659 (2007).

OCIS Codes
(170.3660) Medical optics and biotechnology : Light propagation in tissues
(290.5855) Scattering : Scattering, polarization

ToC Category:
Scattering

History
Original Manuscript: June 19, 2009
Revised Manuscript: August 29, 2009
Manuscript Accepted: August 29, 2009
Published: September 2, 2009

Citation
Tianliang Yun, Nan Zeng, Wei Li, Dongzhi Li, Xiaoyu Jiang, and Hui Ma, "Monte Carlo simulation of polarized photon scattering in anisotropic media," Opt. Express 17, 16590-16602 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-19-16590


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. L. H. Wang, S. L. Jacques, and L.-Q. Zheng, “MCML--Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Methods Programs Biomed. 47(2), 131–146 (1995). [PubMed]
  2. G. W. Kattawar and G. N. Plass, “Radiance and polarization of multiple scattered light from haze and clouds,” Appl. Opt. 7(8), 1519–1527 (1968). [PubMed]
  3. P. Bruscaglioni, G. Zaccanti, and Q. Wei, “Transmission of a pulsed polarized light beam through thick turbid media: numerical results,” Appl. Opt. 32(30), 6142–6150 (1993). [PubMed]
  4. J. Ramella-Roman, S. Prahl, and S. L. Jacques, “Three Monte Carlo programs of polarized light transport into scattering media: part I,” Opt. Express 13(12), 4420–4438 (2005). [PubMed]
  5. J. C. Ramella-Roman, S. A. Prahl, and S. L. Jacques, “Three Monte Carlo programs of polarized light transport into scattering media: part II,” Opt. Express 13(25), 10392–10405 (2005). [PubMed]
  6. 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).
  7. B. D. Cameron, M. J. Rakovic, M. Mehrübeoglu, G. W. Kattawar, S. Rastegar, L. V. Wang, and G. L. Coté, “Measurement and calculation of the two-dimensional backscattering Mueller matrix of a turbid medium,” Opt. Lett. 23(7), 485–487 (1998).
  8. M. J. Rakovi, G. W. Kattawar, M. B. Mehru Beo Lu, B. D. Cameron, L. V. Wang, S. Rastegar, and G. L. Coté, “Light backscattering polarization patterns from turbid media: theory and experiment,” Appl. Opt. 38(15), 3399–3408 (1999).
  9. D. Côté and I. Vitkin, “Robust concentration determination of optically active molecules in turbid media with validated three-dimensional polarization sensitive Monte Carlo calculations,” Opt. Express 13(1), 148–163 (2005). [PubMed]
  10. J. He, A. Karlsson, J. Swartling, and S. Andersson-Engels, “Light scattering by multiple red blood cells,” J. Opt. Soc. Am. A 21, 1953–1961 (2004).
  11. J. D. Keener, K. J. Chalut, J. W. Pyhtila, and A. Wax, “Application of Mie theory to determine the structure of spheroidal scatterers in biological materials,” Opt. Lett. 32(10), 1326–1328 (2007). [PubMed]
  12. A. Kienle and R. Hibst, “Light guiding in biological tissue due to scattering,” Phys. Rev. Lett. 97(1), 018104 (2006). [PubMed]
  13. A. Kienle, C. D’Andrea, F. Foschum, P. Taroni, and A. Pifferi, “Light propagation in dry and wet softwood,” Opt. Express 16(13), 9895–9906 (2008). [PubMed]
  14. A. Kienle, F. K. Forster, and R. Hibst, “Anisotropy of light propagation in biological tissue,” Opt. Lett. 29(22), 2617–2619 (2004). [PubMed]
  15. C. F. Bohren, and D. R. Huffman, Absorption and Scattering of Light by Small Particles (John Wiley and Sons, New York, 1983).
  16. A. Kienle, F. Forster, R. Diebolder, and R. Hibst, “Light propagation in dentin: influence of microstructure on anisotropy,” Phys. Med. Biol. 48(7–N), 14 (2003).
  17. X. Y. Jiang, N. Zeng, Y. H. He, and H. Ma, “Investigation of Linear Polarization Difference Imaging Based on Rotation of Incident and Backscattered Polarization Angles,” Pro. Biochem. Biophys. 34, 659 (2007).

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