OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 19, Iss. 22 — Oct. 24, 2011
  • pp: 21643–21651
« Show journal navigation

Scattering-phase theorem: anomalous diffraction by forward-peaked scattering media

Min Xu  »View Author Affiliations


Optics Express, Vol. 19, Issue 22, pp. 21643-21651 (2011)
http://dx.doi.org/10.1364/OE.19.021643


View Full Text Article

Acrobat PDF (1186 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

The scattering-phase theorem states that the values of scattering and reduced scattering coefficients of the bulk random media are proportional to the variance of the phase and the variance of the phase gradient, respectively, of the phase map of light passing through one thin slice of the medium. We report a new derivation of the scattering phase theorem and provide the correct form of the relation between the variance of phase gradient and the reduced scattering coefficient. We show the scattering-phase theorem is the consequence of anomalous diffraction by a thin slice of forward-peaked scattering media. A new set of scattering-phase relations with relaxed requirement on the thickness of the slice are provided. The condition for the scattering-phase theorem to be valid is discussed and illustrated with simulated data. The scattering-phase theorem is then applied to determine the scattering coefficient μs, the reduced scattering coefficient μs, and the anisotropy factor g for polystyrene sphere and Intralipid-20% suspensions with excellent accuracy from quantitative phase imaging of respective thin slices. The spatially-resolved μs, μs and g maps obtained via such a scattering-phase relationship may find general applications in the characterization of the optical property of homogeneous and heterogeneous random media.

© 2011 OSA

1. Introduction

Multiple scattering of light by random media is ubiquitous in nature. Multiple scattering by tissue, cloud, and other random media withholds direct image of such systems. Indirect characterizing and imaging these systems with multiple scattering light has attracted immense interest due to its practical importance and noninvasiveness nature. In the limit after light being scattered a sufficient number of times, light diffuses in the random medium and light diffusion is characterized by the reduced scattering coefficient μs The transport mean free path, given by the inverse of μs, can be significantly larger than the distance that light travels between consecutive scattering events, μs1, the inverse of the scattering coefficient. Their ratio, μs1/μs1, is typically 10 – 100 in strongly forward-peaked scattering media such as biological tissue probed by visible or near infrared light [1

1. W. F. Cheong, S. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron. 26, 2166–2185 (1990). [CrossRef]

]. In the other extreme, light transmitted through a thin slice of forward scattering media of thickness Lμs1 suffers minimal scattering with its unscattered intensity decreasing according to the Beer’s law. The phase map ϕ of the transmitted light wave can be measured using quantitative phase imaging. The two extreme cases of light propagation in random media–diffusion of multiply scattered light and transmission of minimally scattered light–has been recently suggested inherently connected first by Wang et al. [2

2. Z. Wang, H. Ding, K. Tangella, and G. Popescu, “A scattering-phase theorem,” in Biomedical Optics (BIOMED) Topical Meeting and Tabletop Exhibit, p. BTuD111p (OSA, 2010).

, 3

3. Z. Wang, H. Ding, and G. Popescu, “Scattering-phase theorem.” Opt. Lett. 36(7), 1215–1217 (2011). [CrossRef] [PubMed]

] and later by Iftikhar et al. [4

4. M. Iftikhar, B. DeAngelo, G. Arzumanov, P. Shanley, Z. Xu, and M. Xu, “Characterizing scattering property of random media from phase map of a thin slice: the scattering-phase theorem and the intensity propagation equation approach,” in Optical Tomography and Spectroscopy of Tissue IX, B. J. Tromberg, A. G. Yodh, M. Tamura, E. M. Sevick-Muraca, and R. R. Alfano, eds., vol. 7896 of Proceedings of SPIE, p. 78961O (SPIE, Bellingham, WA, 2011).

]. The values of μs and μs of the bulk media are found to be proportional to the variance of the phase and the variance of the phase gradient, respectively, of the phase map of light passing through one thin slice of the medium. This is so called “scattering-phase theorem.”

In this paper, we report first a new derivation of the scattering phase theorem and provide the correct relation between the variance of phase gradient and μs. The anisotropy factor, g ≡ 1 – μs/μs, an important parameter linked to the morphology of the scatterers in the medium, can then be derived directly from the phase map. More importantly, we show the scattering-phase theorem is the consequence of anomalous diffraction by a thin slice of forward-peaked scattering media. A set of μsϕ, μsϕ, and gϕ relations are provided, for the first time, with relaxed requirement on the thickness of the slice. The condition for the scattering-phase theorem to be valid is discussed and illustrated with simulated data. The scattering-phase theorem is then applied to determine the scattering coefficient, the reduced scattering coefficient and the anisotropy factor for polystyrene sphere and Intralipid-20% suspensions with excellent accuracy from their quantitative phase maps measured by differential interference contrast microscopy. The paper ends with a discussion of the significance and applications of this scattering-phase relationship.

2. Theory

Let’s consider a thin slice of random medium of thickness L illuminated by a plane wave of unit intensity. The spatially resolved phase map ϕ(ρ) for wave transmission is expressed as ϕ(ρ)=k0Ldzm(ρ,z) where k ≡ 2πn0/λ is the wave number, n0 is the background refractive index, λ is the wavelength of light in vacuum, and m is the relative refractive index at position (ρ,z) with ρ and z the lateral and axial coordinates, respectively. The fluctuation in relative refractive index δmm – 1 satisfies 〈δm〉 = 0 where 〈〉 means the spatial average. The phase map ϕ(ρ) can be readily measured with quantitative phase imaging approaches [5

5. A. Barty, K. A. Nugent, D. Paganin, and A. Roberts, “Quantitative optical phase microscopy.” Opt. Lett. 23(11), 817–819 (1998). [CrossRef]

9

9. S. S. Kou, L. Waller, G. Barbastathis, and C. J. R. Sheppard, “Transport-of-intensity approach to differential interference contrast (TI-DIC) microscopy for quantitative phase imaging,” Opt. Lett. 35(3), 447–449 (2010). [CrossRef] [PubMed]

].

The relation between the scattering coefficient μs of the bulk medium and the variance of the phase has been obtained based on the decomposition of the transmitted statistically homogeneous wave field U into its spatial average and a spatially varying component U(ρ) = U0(ρ) +U1(ρ) and the fact that U0 = 〈U〉 corresponds to the unscattered wave and U1 is the scattered component [2

2. Z. Wang, H. Ding, K. Tangella, and G. Popescu, “A scattering-phase theorem,” in Biomedical Optics (BIOMED) Topical Meeting and Tabletop Exhibit, p. BTuD111p (OSA, 2010).

, 10

10. H. Ding, Z. Wang, F. Nguyen, S. A. Boppart, and G. Popescu, “Fourier Transform Light Scattering of Inhomogeneous and Dynamic Structures,” Phys. Rev. Lett. 101(23), 238102 (2008). [CrossRef] [PubMed]

]. When the thickness of the thin slice Lμs1, the intensity of the unscattered wave is expressed as |U0|2 = |〈e(ρ)〉|2 = exp(−μsL) by the Beer’s law. Hence μsL = −2ln|〈eiΔϕ(ρ)〉| where Δϕϕϕ=k0Ldzδm(ρ,z). Since |Δϕ| ≪ 1 as implied by Lμs1, this reduces to
μsL=(Δϕ)2
(1)
if we apply the well-known cumulant expansion theorem [11

11. S. H. Ma, Statistical Mechanics (World Scientific, 1985).

] and write eiΔϕ(ρ)=exp(iΔϕ12(Δϕ)2). The distribution of the phase needs not to follow a Gaussian distribution for Eq. (1) to be valid.

Both relations between the scattering coefficient μs and the variance of the phase, and the reduced scattering coefficient μs and the variance of the phase gradient are the consequence of anomalous diffraction by a thin slice of forward-peaked scattering media and the requirement of Lμs1 can be relaxed. Following the treatment of anomalous diffraction by van de Hulst [12

12. H. C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).

], the scattering amplitude of light into direction θ due to the thin slice is given by
S(θ)=ik22π1+cosθ2(1eiΔϕ(ρ))exp(iksρ)dρ
(2)
where ks is the propagation direction of the scattered light and s is a unit direction vector using the Huygens’ principle [12

12. H. C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).

, 13

13. The factor i appears in Eq. (2) and did not appear in anomalous diffraction by optically soft particles described by Hulst in Ref 12. This difference originates from the fact that the scattering wave into direction θ is proportional to −iS(θ) in the Hulst convention and S(θ) in the contemporary convention adopted here.

]. The presence of the thin slice alters the field on the z = L plane to eiΔϕ(ρ) from 1 and hence the scattered wave is eiΔϕ(ρ) – 1 whereas eiΔϕ(ρ) is the total wave on that plane [12

12. H. C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).

]. We could replace cosθ in Eq. (2) by 1 as scattering is forward-peaked. The scattering cross section Csca = 4πk−2S(0) by the optical extinction theorem is then found to be
Csca=2(1cosΔϕ)dρ.
(3)
The reduced scattering cross section Csca = k−2 ∫(1 – cosθ)|S(θ)|2 dΩ can be simplified by first writing 1cosθ=s2/2 where s is the projection of s on the lateral plane and rewriting Csca as
Csca=18π2dsdρ(1eiΔϕ(ρ))ddρexp(iksρ)×dρ(1eiΔϕ(ρ))ddρexp(iksρ).
(4)
By performing partial integration on d/ and d/′ in Eq. (4) and then integrating over s, we have
Csca=12k2dρdρ[ddρ(1eiΔϕ(ρ))][ddρ(1eiΔϕ(ρ))]δ(ρρ)=12k2|ddρ(1eiΔϕ(ρ))|2dρ,
(5)
which reduces to
Csca=12k2|ϕ|2dρ.
(6)
As the scattering and the reduced scattering cross sections are given by μsAL and μsAL, respectively, by definition for the thin slice of area A, we obtain
μsL=2A(1cosΔϕ)dρ=21cosΔϕ
(7)
and
μsL=12k2A|ϕ|2dρ=12k2|ϕ|2
(8)
from Eqs. (3) and (6). In addition, the anisotropy factor g ≡ 1 – μs/μs, representing the mean cosine of the scattering angle, is given by
g=1|ϕ|24k2(1cosΔϕ)1|ϕ|22k2(Δϕ)2.
(9)

Equations (7), (8) and (9) constitute the main result for the scattering phase theorem. The μsϕ relation (7) reduces to the known expression (1) under the condition μsL ≪ 1, or equivalently, |Δϕ| ≪ 1. These relations share the same origin as the anomalous diffraction by optically soft particles introduced by Hulst [12

12. H. C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).

] which has found wide applications in light scattering [14

14. S. A. Ackerman and G. L. Stephens, “The absorption of solar radiation by cloud droplets: an application of anomalous diffraction theory,” J. Atmos. Sci. 44(12), 1574–1588 (1987). [CrossRef]

18

18. M. Xu and A. Katz, Light Scattering Reviews, vol. III, chap. Statistical Interpretation of Light Anomalous Diffraction by Small Particles and its Applications in Bio-agent Detection and Monitoring, pp. 27–68 (Springer, 2008). [CrossRef]

]. Equations (7), (8) and (9) are valid for forward-peaked scattering media as long as the ray does not deviate from the forward direction. The scattering-phase theorem is applicable to a slice of homogeneous or inhomogeneous medium. In the latter case, a map of μs, μs and g can be computed from the phase map using spatial averaging over local regions rather than the whole slice.

3. Simulations and experiments

We performed simulations to validate the scattering-phase theorem for a random medium. In simulation, the fluctuation of the refractive index of the medium Rn(r) = 〈δm(r′)δm(r′ + r)〉 is assumed to be the Whittle-Matern correlation function [19

19. P. Guttorp and T. Gneiting, “On the Whittle-Matrn correlation family,” Tech. Rep. NRCSE-TRS No. 080, NRCSE, University of Washington (2005).

] given by:
Rn(r)=(δm)2γ(rl)
(10)
with
γ(rl)=21ν|Γ(ν)|1(rl)νKν(rl)
(11)
where Kν (·) is the modified Bessel function of the second kind. The Whittle-Matern correlation function has been used extensively to model turbulence and refractive index fluctuation in biological tissue [20

20. J. M. Schmitt and G. Kumar, “Turbulent nature of refractive-index variations in biological tissue,” Opt. Lett. 21, 1310–1312 (1996). [CrossRef] [PubMed]

, 21

21. V. Turzhitsky, A. Radosevich, J. D. Rogers, A. Taflove, and V. Backman, “A predictive model of backscattering at subdiffusion length scales,” Biomed. Opt. Express 1(3), 1034–1046 (2010), URL http://www.opticsinfobase.org/boe/abstract.cfm?URI=boe-1-3-1034.

]. The typical values are 〈(δm)〉2 = 0.012, l ∼ 0.5μm, and n0 = 1.367 for biological tissue [22

22. M. Xu, T. T. Wu, and J. Y. Qu, “Unified Mie and fractal scattering by cells and experimental study on application in optical characterization of cellular and subcellular structures,” J. Biomed. Opt. 13, 038802 (2008).

, 23

23. T. T. Wu, J. Y. Qu, and M. Xu, “Unified Mie and fractal scattering by biological cells and subcellular structures,” Opt. Lett. 32, 2324–2326 (2007). [CrossRef] [PubMed]

]. The Fourier transform of the correlation function is given by
R^n(q)=(δm)2Γ(ν+3/2)π3/2|Γ(ν)|l3(1+q2l2)ν3/2
(12)
when ν > −3/2. Light scattering by the random medium is fully described by the power spectrum of the fluctuation of the refractive index. Following [22

22. M. Xu, T. T. Wu, and J. Y. Qu, “Unified Mie and fractal scattering by cells and experimental study on application in optical characterization of cellular and subcellular structures,” J. Biomed. Opt. 13, 038802 (2008).

,24

24. M. Xu and R. R. Alfano, “Fractal mechanisms of light scattering in biological tissue and cells,” Opt. Lett. 30, 3051–3053 (2005). [CrossRef] [PubMed]

], the scattering coefficient and the reduced scattering coefficient are given by
μs=2π1/2k2l(δm)2Γ(ν+1/2)|Γ(ν)|[1(1+4X2)ν1/2],
(13)
and
μs=π1/2l1(δm)2Γ(ν+1/2)|Γ(ν)|1ν1/2×[1(1+4X2)ν1/2[1+4X2(ν+1/2)]],
(14)
respectively, in this model where Xkl is the size parameter.

We set the strength of refractive index fluctuation 〈(δm)2〉 = 0.012, the correlation length l = 0.5μm, the background refractive index of the sample n0 = 1.367, and the wavelength of the incident beam λ = 0.5μm in the simulation. The random field inside a box of size 10l × 10l × L with varying thickness L = l, 5l, 20l, and 100l was simulated using RandomFields [25

25. M. Schlather, “An introduction to positive-definite functions and to unconditional simulation of random fields,” Tech. Rep. ST-99-10, Lancaster University (1999).

] with a specified spacial resolution. The phase map was generated by line integration. The gradient of the phase was computed from the phase map using the finite difference between neighboring phases. Total 15 simulations were performed for each set of parameters with their mean and standard deviation being reported hereafter.

Figure 1, from left- to right-hand direction, displays the normalized phase map ( ΔOPL/δm2l) where the optical path length fluctuation is given by ΔOPL=0Ldzδm(ρ,z), the ratio of 2 〈1 – cosΔϕ〉 over μsL, and the ratio of (2k2)−1 〈|∇ϕ|2〉 over μsL for various ν. The normalized phase map is shown for thin slices of thickness L = 20l. The scattering coefficient and the anisotropy factor for the bulk random medium are, μsl = 0.023 and g = 0.988 in the case of ν = 1.0, μsl = 0.015 and g = 0.968 in the case of ν = 0.5, and μsl = 0.0040 and g = 0.915 in the case of ν = 0.1, respectively. The thickness of the samples covers the range starting from μsL ≪ 1 to μsL > 1. The two ratios 2〈1 – cosΔϕ〉/μsL and (2k2)−1 〈|∇ϕ|2〉/μsL are expected to be unity according to Eqs. (7) and (8). Figure 3 shows the former ratio approaches unity when the thickness of the medium is at least 5l. The value of μs can be computed from the phase map at all levels of resolution. On the other hand, the resolution matters for probing μs. The latter ratio approaches unity and the best estimation for μs is obtained only when the resolution is 0.1l − 0.2l and the thickness L ≥ 5l. Insufficient resolution results in an underestimation of μs.

We then examined the light scattering properties of polystyrene sphere and Intralipid-20% suspensions by applying the scattering-phase theorem to the quantitative phase map of respective thin slice measured with a differential interference contrast (DIC) microscope (Axiovert 40CFL, Zeiss). The light source was a Halogen 35W lamp filtered by a 550nm narrow-band filter. The numerical aperture for the condenser and objective (APlan 40×) were 0.2 and 0.5, respectively. The pixel size for the recorded images was 0.064μm using Canon 5D Mark II. The quantitative phase map for a monolayer of polystyrene sphere suspension (size: 8.31μm) in water and a thin film (thickness: 4μm) of Intralipid-20% suspension on a glass microscope slide were computed from in-focus and out-of-focus (δz = 1μm) DIC images under Köhler illumination using the transport-of-intensity approach [9

9. S. S. Kou, L. Waller, G. Barbastathis, and C. J. R. Sheppard, “Transport-of-intensity approach to differential interference contrast (TI-DIC) microscopy for quantitative phase imaging,” Opt. Lett. 35(3), 447–449 (2010). [CrossRef] [PubMed]

]. Figure 2 shows the computed optical path length maps ΔOPL for the two samples. The scattering property for each individual spheres can be analyzed by applying the scattering-phase theorem to the region in the phase map being occupied by the sphere. For example, the region highlighted by white dash lines for the central sphere yields μs = 0.234μm−1, μs = 0.0202μm−1 and g = 0.91 with an area 61.0μm2. The scattering and reduced scattering cross sections are 118μm2 and 10.2μm2. The mean scattering and reduced scattering cross sections for all the spheres contained in the displayed section are 116μm2 and 9.8μm2, respectively. These values are in excellent agreement with the theoretical prediction for a polystyrene sphere of the specified size (Csca = 125μm2, Csca = 9.5μm2 and g = 0.92) computed with a Mie code [26

26. W. J. Wiscombe, “Improved Mie scattering algorithms,” Appl. Opt. 19, 1505–1509 (1980). [CrossRef] [PubMed]

].

Fig. 2 The optical path length map ΔOPL for a monolayer of polystyrene sphere suspension (size: 8.31μm) in water (left) and a thin film (thickness: 4μm) of Intralipid-20% suspension (right).

The scattering and the reduced scattering coefficients for Intralipid-20% suspension are found to be 0.136μm−1 and 0.001μm−1, respectively, from the whole section displayed in Fig. 2. The former agrees with the known μs value (0.139μm−1) whereas the latter dramatically underestimates μs (0.031μm−1) at 550nm [27

27. H. J. van Staveren, C. J. M. Moes, J. van Marle, S. A. Prahl, and M. J. C. van GemertJ, “Light scattering in Intralipid-10% in the wavelength range of 400–1100nm,” Appl. Opt. 30(31), 4507–4514 (1991). [CrossRef] [PubMed]

, 28

28. A. Giusto, R. Saija, M. A. Iat, P. Denti, F. Borghese, and O. I. Sindoni, “Optical properties of high-density dispersions of particles: application to intralipid solutions,” Appl. Opt. 42(21), 4375–4380 (2003). [CrossRef] [PubMed]

]. This behavior is expected as the characteristic correlation length for the Intralipid suspension is sub-wavelength [29

29. M. Xu, M. Alrubaiee, and R. R. Alfano, “Fractal mechanism of light scattering for tissue optical biopsy,” in Optical Biopsy VI, R. R. Alfano and A. Katz, eds., vol. 6091 of Proceedings of SPIE, p. 60910E (2006).

] and the resolution of the phase map is insufficient to provide an accurate estimation of μs directly (see Fig. 1). The quality of μs estimation, however, can be significantly improved by properly taking into account light diffraction in the microscope and sharpening the phase map accordingly. This procedure yields the new value of μs to be 0.022μm−1, agreeing reasonably well with the real value. The detail will be published elsewhere.

Fig. 1 The normalized phase map ( ΔOPL/δm2l), the ratio of 2〈1 – cosΔϕ〉 over μsL, and the ratio of (2k2)–1 〈|∇ϕ|2〉 over μsL are displayed, from left- to right-hand direction, for a thin slice of random medium of varying thickness L with the refractive index fluctuation following the Whittle-Matern correlation function of ν = 1.0 (top row), ν = 0.5 (middle row) and ν = 0.1 (bottom row). The normalized phase map is shown for L = 20l.

4. Discussion

The μsϕ and μsϕ relations can be justified intuitively as the following. Light scattering (μs) depends on the fluctuation of the refractive index which emerges as the variance in the phase map for light transmission through a thin slice. Light reduced scattering (μs) reflects the deviation of the equal-phase wave front away from the forward direction which is described by the local tilt (gradient) in the phase for light transmission through a thin slice. Assuming the thin slice of sample of thickness L is uniformly divided into N = L/l layers with l the correlation length of the random medium, Δϕ (and ∇ϕ) is the summation of N independent random numbers from the N layers. Hence the spatial average 〈(Δϕ)2〉 (and 〈|∇ϕ|2〉) scales with N rather than N2. These considerations lead to μsL ∝ (Δϕ)2 and μsL ∝ 〈|∇ϕ|2〉. In cases such as a monolayer of scatterers of size much larger than the wavelength, the condition that |Δϕ| ≪ 1 is not satisfied, the more general μsϕ relation (7) should be used whereas the μsϕ relation remains the same provided the rays do not deviate from the forward direction (the scatterers are optically soft).

The above argument also explains that the thickness of the sample should be at least multiple l (with a sufficient large N) to obtain the values of μs and μs correctly from the phase map as observed in the simulation. To properly compute the local tilt in the phase to obtain μs with finite difference, the separation between the two points must be smaller than the size of the scattering structure. The separation at the order of 0.1l – 0.2l may be optimal as suggested by the simulation.

Finally, we would like to point out that the limiting form of the scattering-phase theorem when |Δϕ| ≪ 1 has also been obtained previously by us using another approach [4

4. M. Iftikhar, B. DeAngelo, G. Arzumanov, P. Shanley, Z. Xu, and M. Xu, “Characterizing scattering property of random media from phase map of a thin slice: the scattering-phase theorem and the intensity propagation equation approach,” in Optical Tomography and Spectroscopy of Tissue IX, B. J. Tromberg, A. G. Yodh, M. Tamura, E. M. Sevick-Muraca, and R. R. Alfano, eds., vol. 7896 of Proceedings of SPIE, p. 78961O (SPIE, Bellingham, WA, 2011).

] through analyzing the cross correlation 〈Δϕ(ρϕ(ρ′)〉 between two points ρ and ρ′ on the phase map for light transmission through a thin slice of a weakly scattering random medium. The scattering-phase theorem in this limit is equivalent to
Δϕ(ρ)Δϕ(ρ)=μsL12(kΔρ)2μsL
(15)
where Δρ ≡ |ρρ′| is the distance between the two points [4

4. M. Iftikhar, B. DeAngelo, G. Arzumanov, P. Shanley, Z. Xu, and M. Xu, “Characterizing scattering property of random media from phase map of a thin slice: the scattering-phase theorem and the intensity propagation equation approach,” in Optical Tomography and Spectroscopy of Tissue IX, B. J. Tromberg, A. G. Yodh, M. Tamura, E. M. Sevick-Muraca, and R. R. Alfano, eds., vol. 7896 of Proceedings of SPIE, p. 78961O (SPIE, Bellingham, WA, 2011).

]. Ref [2

2. Z. Wang, H. Ding, K. Tangella, and G. Popescu, “A scattering-phase theorem,” in Biomedical Optics (BIOMED) Topical Meeting and Tabletop Exhibit, p. BTuD111p (OSA, 2010).

, 3

3. Z. Wang, H. Ding, and G. Popescu, “Scattering-phase theorem.” Opt. Lett. 36(7), 1215–1217 (2011). [CrossRef] [PubMed]

] presented a different expression for the gϕ relation. The difference originates from the scattered wave was assumed to be eiΔϕ(ρ) in our notation in Ref [2

2. Z. Wang, H. Ding, K. Tangella, and G. Popescu, “A scattering-phase theorem,” in Biomedical Optics (BIOMED) Topical Meeting and Tabletop Exhibit, p. BTuD111p (OSA, 2010).

, 3

3. Z. Wang, H. Ding, and G. Popescu, “Scattering-phase theorem.” Opt. Lett. 36(7), 1215–1217 (2011). [CrossRef] [PubMed]

]. Since the presence of the thin slice alters the field on the z = L plane to eiΔϕ(ρ) from 1, the scattered wave is [eiΔϕ(ρ) – 1] whereas eiΔϕ(ρ) is the total wave on that plane. The probability density for light scattering into direction q = ks, hence, is given by
P(q)=|U˜(q)|2|U˜(q)|2dq
(16)
where
U˜(q)=1(2π)2[eiΔϕ(ρ)1]exp(iρq)dρ.
(17)
One could follow the procedure outlined in Ref [2

2. Z. Wang, H. Ding, K. Tangella, and G. Popescu, “A scattering-phase theorem,” in Biomedical Optics (BIOMED) Topical Meeting and Tabletop Exhibit, p. BTuD111p (OSA, 2010).

,3

3. Z. Wang, H. Ding, and G. Popescu, “Scattering-phase theorem.” Opt. Lett. 36(7), 1215–1217 (2011). [CrossRef] [PubMed]

] and reach the gϕ relation (9) if the correct probability density Eq. (16) for light scattering into direction q is used.

5. Conclusion

In summary, we have derived the scattering phase theorem and provided the correct relation between the variance of phase gradient and the reduced scattering coefficient. More importantly, the scattering-phase theorem is shown to be the consequence of anomalous diffraction by a thin slice of forward-peaked scattering media. A set of μsϕ, μsϕ, and gϕ relations have been provided, for the first time, with relaxed requirement on the thickness of the slice. The condition for the scattering-phase theorem to be valid has been discussed and illustrated with simulated data. The scattering-phase theorem has been applied to determine successfully the scattering coefficient, the reduced scattering coefficient and the anisotropy factor for polystyrene sphere and Intralipid-20% suspensions from their respective quantitative phase map of a thin slice.

The characterization of the scattering properties (μs, μs, and g) of biological tissue and cells has been a challenging and important problem in biomedical optics [30

30. S. Menon, Q. Su, and R. Grobe, “Determination of g and μ using multiply scattered light in turbid media,” Phys. Rev. Lett. 94, 153904 (2005). [CrossRef] [PubMed]

]. This scattering-phase relationship establishes a new means to characterize the scattering properties of these samples. The spatially-resolved μs, μs and g maps obtained via such a scattering-phase relationship will provide detailed local maps for scattering structures which may be of important diagnosis value, and may find applications in the characterization of the optical property of homogeneous and heterogeneous random media in general.

Acknowledgments

MX acknowledges Research Corporation, NIH ( 1R15EB009224) and DOD ( W81XWH-10-1-0526) for their support.

References and links

1.

W. F. Cheong, S. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron. 26, 2166–2185 (1990). [CrossRef]

2.

Z. Wang, H. Ding, K. Tangella, and G. Popescu, “A scattering-phase theorem,” in Biomedical Optics (BIOMED) Topical Meeting and Tabletop Exhibit, p. BTuD111p (OSA, 2010).

3.

Z. Wang, H. Ding, and G. Popescu, “Scattering-phase theorem.” Opt. Lett. 36(7), 1215–1217 (2011). [CrossRef] [PubMed]

4.

M. Iftikhar, B. DeAngelo, G. Arzumanov, P. Shanley, Z. Xu, and M. Xu, “Characterizing scattering property of random media from phase map of a thin slice: the scattering-phase theorem and the intensity propagation equation approach,” in Optical Tomography and Spectroscopy of Tissue IX, B. J. Tromberg, A. G. Yodh, M. Tamura, E. M. Sevick-Muraca, and R. R. Alfano, eds., vol. 7896 of Proceedings of SPIE, p. 78961O (SPIE, Bellingham, WA, 2011).

5.

A. Barty, K. A. Nugent, D. Paganin, and A. Roberts, “Quantitative optical phase microscopy.” Opt. Lett. 23(11), 817–819 (1998). [CrossRef]

6.

M. R. Arnison, K. G. Larkin, C. J. R. Sheppard, N. I. Smith, and C. J. Cogswell, “Linear phase imaging using differential interference contrast microscopy.” J. Microsc. 214(1), 7–12 (2004). [CrossRef] [PubMed]

7.

G. Popescu, T. Ikeda, R. R. Dasari, and M. S. Feld, “Diffraction phase microscopy for quantifying cell structure and dynamics,” Opt. Lett. 31(6), 775–777 (2006). [CrossRef] [PubMed]

8.

W. S. Rockward, A. L. Thomas, B. Zhao, and C. A. DiMarzio, “Quantitative phase measurements using optical quadrature microscopy,” Appl. Opt. 47(10), 1684–1696 (2008). [CrossRef] [PubMed]

9.

S. S. Kou, L. Waller, G. Barbastathis, and C. J. R. Sheppard, “Transport-of-intensity approach to differential interference contrast (TI-DIC) microscopy for quantitative phase imaging,” Opt. Lett. 35(3), 447–449 (2010). [CrossRef] [PubMed]

10.

H. Ding, Z. Wang, F. Nguyen, S. A. Boppart, and G. Popescu, “Fourier Transform Light Scattering of Inhomogeneous and Dynamic Structures,” Phys. Rev. Lett. 101(23), 238102 (2008). [CrossRef] [PubMed]

11.

S. H. Ma, Statistical Mechanics (World Scientific, 1985).

12.

H. C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).

13.

The factor i appears in Eq. (2) and did not appear in anomalous diffraction by optically soft particles described by Hulst in Ref 12. This difference originates from the fact that the scattering wave into direction θ is proportional to −iS(θ) in the Hulst convention and S(θ) in the contemporary convention adopted here.

14.

S. A. Ackerman and G. L. Stephens, “The absorption of solar radiation by cloud droplets: an application of anomalous diffraction theory,” J. Atmos. Sci. 44(12), 1574–1588 (1987). [CrossRef]

15.

P. Chýlek and J. D. Klett, “Extinction cross sections of nonspherical particles in the anomalous diffraction approximation,” J. Opt. Soc. Am. A 8, 274–281 (1991). [CrossRef]

16.

M. Xu, M. Lax, and R. R. Alfano, “Light anomalous diffraction using geometrical path statistics of rays and Gaussian ray approximation,” Opt. Lett. 28, 179–181 (2003). [CrossRef] [PubMed]

17.

P. Yang, Z. Zhang, B. Baum, H.-L. Huang, and Y. Hu, “A new look at anomalous diffraction theory (ADT): Algorithm in cumulative projected-area distribution domain and modified ADT,” J. Quant. Spectrosc. Radiat. Transf. 89, 421–442 (2004). [CrossRef]

18.

M. Xu and A. Katz, Light Scattering Reviews, vol. III, chap. Statistical Interpretation of Light Anomalous Diffraction by Small Particles and its Applications in Bio-agent Detection and Monitoring, pp. 27–68 (Springer, 2008). [CrossRef]

19.

P. Guttorp and T. Gneiting, “On the Whittle-Matrn correlation family,” Tech. Rep. NRCSE-TRS No. 080, NRCSE, University of Washington (2005).

20.

J. M. Schmitt and G. Kumar, “Turbulent nature of refractive-index variations in biological tissue,” Opt. Lett. 21, 1310–1312 (1996). [CrossRef] [PubMed]

21.

V. Turzhitsky, A. Radosevich, J. D. Rogers, A. Taflove, and V. Backman, “A predictive model of backscattering at subdiffusion length scales,” Biomed. Opt. Express 1(3), 1034–1046 (2010), URL http://www.opticsinfobase.org/boe/abstract.cfm?URI=boe-1-3-1034.

22.

M. Xu, T. T. Wu, and J. Y. Qu, “Unified Mie and fractal scattering by cells and experimental study on application in optical characterization of cellular and subcellular structures,” J. Biomed. Opt. 13, 038802 (2008).

23.

T. T. Wu, J. Y. Qu, and M. Xu, “Unified Mie and fractal scattering by biological cells and subcellular structures,” Opt. Lett. 32, 2324–2326 (2007). [CrossRef] [PubMed]

24.

M. Xu and R. R. Alfano, “Fractal mechanisms of light scattering in biological tissue and cells,” Opt. Lett. 30, 3051–3053 (2005). [CrossRef] [PubMed]

25.

M. Schlather, “An introduction to positive-definite functions and to unconditional simulation of random fields,” Tech. Rep. ST-99-10, Lancaster University (1999).

26.

W. J. Wiscombe, “Improved Mie scattering algorithms,” Appl. Opt. 19, 1505–1509 (1980). [CrossRef] [PubMed]

27.

H. J. van Staveren, C. J. M. Moes, J. van Marle, S. A. Prahl, and M. J. C. van GemertJ, “Light scattering in Intralipid-10% in the wavelength range of 400–1100nm,” Appl. Opt. 30(31), 4507–4514 (1991). [CrossRef] [PubMed]

28.

A. Giusto, R. Saija, M. A. Iat, P. Denti, F. Borghese, and O. I. Sindoni, “Optical properties of high-density dispersions of particles: application to intralipid solutions,” Appl. Opt. 42(21), 4375–4380 (2003). [CrossRef] [PubMed]

29.

M. Xu, M. Alrubaiee, and R. R. Alfano, “Fractal mechanism of light scattering for tissue optical biopsy,” in Optical Biopsy VI, R. R. Alfano and A. Katz, eds., vol. 6091 of Proceedings of SPIE, p. 60910E (2006).

30.

S. Menon, Q. Su, and R. Grobe, “Determination of g and μ using multiply scattered light in turbid media,” Phys. Rev. Lett. 94, 153904 (2005). [CrossRef] [PubMed]

OCIS Codes
(180.3170) Microscopy : Interference microscopy
(290.5820) Scattering : Scattering measurements
(290.7050) Scattering : Turbid media
(290.5825) Scattering : Scattering theory

ToC Category:
Scattering

History
Original Manuscript: July 28, 2011
Revised Manuscript: September 16, 2011
Manuscript Accepted: September 21, 2011
Published: October 19, 2011

Virtual Issues
Vol. 6, Iss. 11 Virtual Journal for Biomedical Optics

Citation
Min Xu, "Scattering-phase theorem: anomalous diffraction by forward-peaked scattering media," Opt. Express 19, 21643-21651 (2011)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-22-21643


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. W. F. Cheong, S. Prahl, and A. J. Welch, “A review of the optical properties of biological tissues,” IEEE J. Quantum Electron.26, 2166–2185 (1990). [CrossRef]
  2. Z. Wang, H. Ding, K. Tangella, and G. Popescu, “A scattering-phase theorem,” in Biomedical Optics (BIOMED) Topical Meeting and Tabletop Exhibit, p. BTuD111p (OSA, 2010).
  3. Z. Wang, H. Ding, and G. Popescu, “Scattering-phase theorem.” Opt. Lett.36(7), 1215–1217 (2011). [CrossRef] [PubMed]
  4. M. Iftikhar, B. DeAngelo, G. Arzumanov, P. Shanley, Z. Xu, and M. Xu, “Characterizing scattering property of random media from phase map of a thin slice: the scattering-phase theorem and the intensity propagation equation approach,” in Optical Tomography and Spectroscopy of Tissue IX, B. J. Tromberg, A. G. Yodh, M. Tamura, E. M. Sevick-Muraca, and R. R. Alfano, eds., vol. 7896 of Proceedings of SPIE, p. 78961O (SPIE, Bellingham, WA, 2011).
  5. A. Barty, K. A. Nugent, D. Paganin, and A. Roberts, “Quantitative optical phase microscopy.” Opt. Lett.23(11), 817–819 (1998). [CrossRef]
  6. M. R. Arnison, K. G. Larkin, C. J. R. Sheppard, N. I. Smith, and C. J. Cogswell, “Linear phase imaging using differential interference contrast microscopy.” J. Microsc.214(1), 7–12 (2004). [CrossRef] [PubMed]
  7. G. Popescu, T. Ikeda, R. R. Dasari, and M. S. Feld, “Diffraction phase microscopy for quantifying cell structure and dynamics,” Opt. Lett.31(6), 775–777 (2006). [CrossRef] [PubMed]
  8. W. S. Rockward, A. L. Thomas, B. Zhao, and C. A. DiMarzio, “Quantitative phase measurements using optical quadrature microscopy,” Appl. Opt.47(10), 1684–1696 (2008). [CrossRef] [PubMed]
  9. S. S. Kou, L. Waller, G. Barbastathis, and C. J. R. Sheppard, “Transport-of-intensity approach to differential interference contrast (TI-DIC) microscopy for quantitative phase imaging,” Opt. Lett.35(3), 447–449 (2010). [CrossRef] [PubMed]
  10. H. Ding, Z. Wang, F. Nguyen, S. A. Boppart, and G. Popescu, “Fourier Transform Light Scattering of Inhomogeneous and Dynamic Structures,” Phys. Rev. Lett.101(23), 238102 (2008). [CrossRef] [PubMed]
  11. S. H. Ma, Statistical Mechanics (World Scientific, 1985).
  12. H. C. van de Hulst, Light Scattering by Small Particles (Dover, 1981).
  13. The factor i appears in Eq. (2) and did not appear in anomalous diffraction by optically soft particles described by Hulst in Ref 12. This difference originates from the fact that the scattering wave into direction θ is proportional to −iS(θ) in the Hulst convention and S(θ) in the contemporary convention adopted here.
  14. S. A. Ackerman and G. L. Stephens, “The absorption of solar radiation by cloud droplets: an application of anomalous diffraction theory,” J. Atmos. Sci.44(12), 1574–1588 (1987). [CrossRef]
  15. P. Chýlek and J. D. Klett, “Extinction cross sections of nonspherical particles in the anomalous diffraction approximation,” J. Opt. Soc. Am. A8, 274–281 (1991). [CrossRef]
  16. M. Xu, M. Lax, and R. R. Alfano, “Light anomalous diffraction using geometrical path statistics of rays and Gaussian ray approximation,” Opt. Lett.28, 179–181 (2003). [CrossRef] [PubMed]
  17. P. Yang, Z. Zhang, B. Baum, H.-L. Huang, and Y. Hu, “A new look at anomalous diffraction theory (ADT): Algorithm in cumulative projected-area distribution domain and modified ADT,” J. Quant. Spectrosc. Radiat. Transf.89, 421–442 (2004). [CrossRef]
  18. M. Xu and A. Katz, Light Scattering Reviews, vol. III, chap. Statistical Interpretation of Light Anomalous Diffraction by Small Particles and its Applications in Bio-agent Detection and Monitoring, pp. 27–68 (Springer, 2008). [CrossRef]
  19. P. Guttorp and T. Gneiting, “On the Whittle-Matrn correlation family,” Tech. Rep. NRCSE-TRS No. 080, NRCSE, University of Washington (2005).
  20. J. M. Schmitt and G. Kumar, “Turbulent nature of refractive-index variations in biological tissue,” Opt. Lett.21, 1310–1312 (1996). [CrossRef] [PubMed]
  21. V. Turzhitsky, A. Radosevich, J. D. Rogers, A. Taflove, and V. Backman, “A predictive model of backscattering at subdiffusion length scales,” Biomed. Opt. Express1(3), 1034–1046 (2010), URL http://www.opticsinfobase.org/boe/abstract.cfm?URI=boe-1-3-1034 .
  22. M. Xu, T. T. Wu, and J. Y. Qu, “Unified Mie and fractal scattering by cells and experimental study on application in optical characterization of cellular and subcellular structures,” J. Biomed. Opt.13, 038802 (2008).
  23. T. T. Wu, J. Y. Qu, and M. Xu, “Unified Mie and fractal scattering by biological cells and subcellular structures,” Opt. Lett.32, 2324–2326 (2007). [CrossRef] [PubMed]
  24. M. Xu and R. R. Alfano, “Fractal mechanisms of light scattering in biological tissue and cells,” Opt. Lett.30, 3051–3053 (2005). [CrossRef] [PubMed]
  25. M. Schlather, “An introduction to positive-definite functions and to unconditional simulation of random fields,” Tech. Rep. ST-99-10, Lancaster University (1999).
  26. W. J. Wiscombe, “Improved Mie scattering algorithms,” Appl. Opt.19, 1505–1509 (1980). [CrossRef] [PubMed]
  27. H. J. van Staveren, C. J. M. Moes, J. van Marle, S. A. Prahl, and M. J. C. van GemertJ, “Light scattering in Intralipid-10% in the wavelength range of 400–1100nm,” Appl. Opt.30(31), 4507–4514 (1991). [CrossRef] [PubMed]
  28. A. Giusto, R. Saija, M. A. Iat, P. Denti, F. Borghese, and O. I. Sindoni, “Optical properties of high-density dispersions of particles: application to intralipid solutions,” Appl. Opt.42(21), 4375–4380 (2003). [CrossRef] [PubMed]
  29. M. Xu, M. Alrubaiee, and R. R. Alfano, “Fractal mechanism of light scattering for tissue optical biopsy,” in Optical Biopsy VI, R. R. Alfano and A. Katz, eds., vol. 6091 of Proceedings of SPIE, p. 60910E (2006).
  30. S. Menon, Q. Su, and R. Grobe, “Determination of g and μ using multiply scattered light in turbid media,” Phys. Rev. Lett.94, 153904 (2005). [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.

Figures

Fig. 2 Fig. 1
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited