9 May 2016 Retardance of bilayer anisotropic samples consisting of well-aligned cylindrical scatterers and birefringent media
Author Affiliations +
J. of Biomedical Optics, 21(5), 055002 (2016). doi:10.1117/1.JBO.21.5.055002
Abstract
Both cylindrical scatterers and birefringent media may contribute to the anisotropy of tissue, where anisotropy can be characterized using polarization techniques. Our previous studies have shown that a layer of well-aligned cylindrical scatterers displays anisotropic properties similar to those of a piece of birefringent media, whose equivalent extraordinary axis is along the axial direction of the cylinders. We focused on a sample consisting of two layers of anisotropic media, with each layer having a different orientation; the characteristics of this sample were representative of the properties of multilayer fibrous tissues. Using a Mueller matrix decomposition method, we examined in detail how the total retardance and the equivalent extraordinary axis of the bilayered sample varied with changes in the retardance of the two layers and the direction of the extraordinary axis. The results of this study showed that, in such bilayer samples, a layer of well-aligned cylindrical scatterers generated a retardance that behaved exactly like the retardance generated by a piece of birefringent media. The simulated results were also confirmed by the results of experiments using aligned glass fibers.
Guo, Zeng, He, Liu, Du, He, and Ma: Retardance of bilayer anisotropic samples consisting of well-aligned cylindrical scatterers and birefringent media

1.

Introduction

Many biological tissues contain large amounts of fibrous components that exhibit birefringent properties; examples of such components are collagen fibers in tendons, intracellular microtubules, and microfilaments.12.3 In a recent study,4 we showed that scattering by well-aligned cylindrical scatterers displayed anisotropic effects, including a retardance that was similar to that of a piece of birefringent media whose extraordinary axis, and the fast axis, was along the axial direction of the aligned cylinders. We also examined the different contributions to the polarization features provided by the well-aligned cylindrical scatterers and the anisotropic interstitial medium.56.7 Several other groups have tried to quantify the direction of a birefringent media based on polarized light imaging and studying the three-dimensional (3-D) orientation of the optical axis of birefringent turbid media.8,9 It has been shown that for biological tissues containing several layers of different anisotropic properties, the axial variation in the birefringence can result in ambiguity when the measurement is averaged over the sampling volume.1011.12

In this study, we examined in detail how well-aligned cylindrical scatterers and birefringent media contributed to the value and direction of the total retardance of the media. Using Monte Carlo simulations,13 we simulated how the total retardance and the equivalent extraordinary axis of the bilayer sample changed with the change in the value of the retardance and directions of the extraordinary axis of the two anisotropic layers, respectively. The results showed that the retardance of the well-aligned cylindrical scatterers could be regarded as the retardance generated by a piece of birefringent media, whose equivalent extraordinary axis is aligned with the axial direction of the cylinders. We started with the superposition of two layers of uniaxial birefringent media whose axes were in the same plane, perpendicular to the propagation direction of the incident light. Furthermore, we showed that the characteristic features of the bilayer birefringent sample (BBS) were maintained even if one or both of the birefringent layers were replaced by well-aligned cylindrical scatterers.

2.

Theory and Method

2.1.

Mueller Matrix Decomposition Method

Lu and Chipman14 proposed a Mueller matrix containing three physical properties that can be extracted via a polar decomposition process. A Mueller matrix M can be regarded as the product of three constituent “basic” matrices, representing depolarization (MΔ), retardance (MR), and diattenuation (MD), respectively

(1)

M=MΔMRMD.

The depolarization coefficient Δ can be calculated from the elements of MΔ and the linear retardance δ can be calculated from the elements of MR

(2)

Δ=1|trace(mΔ)|3,0<Δ<1,

(3)

δ=cos1[Tr(MR)21],
where MΔ is the submatrix of MΔ. The elements of the retardance vector R=[1r1r2r3] can be calculated as15

(4)

ri=12sinδj,k=13ϵijk(mR)jk,
where ϵijk is the Levi-Civita permutation symbol and MR is the submatrix of MR. θ is the extraordinary axis direction of the linear retardance, which can be calculated as

(5)

θ=12tan1{r2r1}.

2.2.

Monte Carlo Simulations and Scattering Model

In earlier studies, we developed polarized photon scattering models and the corresponding Monte Carlo simulation program for examining the propagation of polarized photons in an anisotropic, turbid sample.13 The simulation program tracks the trajectory and polarization states of the photons; corresponding illuminations of at least four different polarization states of the forward scattering photons are recorded to evaluate the Mueller matrix of the sample

(6)

Sout=M·Sin(IQUV)out=(m11m12m13m14m21m22m23m24m31m32m33m34m41m42m43m44)(IQUV)in.
The scattering models consist of spherical scatterers and infinitely long cylindrical scatterers, as well as interstitial birefringent media.45.6.7,13 The diameter and density of the spherical and cylindrical scatters, the degree of the alignment and spatial orientation of the cylindrical scatters, the optical properties of the scatterers, and the interstitial media were all variables used in the simulations to mimic the properties of turbid media such as tissues.

3.

Results and Discussion

3.1.

Bilayer Birefringent Sample

Typically, the Mueller matrix of a linear retarder can be written as in Ref. 16:

(7)

MR=[10000cos22θ+sin22θcosδsin2θcos  2θ(1cosδ)sin2θsinδ0sin2θcos2θ(1cosδ)sin22θ+cos22θcosδcos2θsinδ0sin2θsinδcos2θsinδcosδ].

We considered a sample consisting of two parallel birefringent layers with retardance values of δ1 and δ2 and extraordinary axes at θ1 and θ2, respectively. Figure 1 shows a schematic illustration of the two layers of birefringent media superposed in the same optical path; here, the planes including both the ordinary and extraordinary axes are parallel to each other and perpendicular to the light propagation direction.

Fig. 1

3-D BBS model.

JBO_21_5_055002_f001.png

The Mueller matrix of such a BBS is

(8)

MR(δ,θ)=MR2(δ2,θ2)·MR1(δ1,θ1).
Based on Eqs. (3), (5), and (7), the δ and θ of MR could be calculated. Assuming that δ1 and θ1 were constant, we then deduced δ and θ as a function of δ2 and θ2, as shown in Fig. 2.

Fig. 2

Keeping δ1 constant, the total retardance δ and orientation of the equivalent extraordinary axis θ changed with changes in δ2; the differently colored lines show data for different angles (0 deg to 90 deg) between the two extraordinary axes of the two layers of birefringent media.

JBO_21_5_055002_f002.png

Figure 2 shows how the total retardance δ and the equivalent extraordinary axis orientation θ varied with changes in the retardance and orientation of the two birefringent layers. In the calculations, the retardance and the orientation of the first layer (δ1,θ1) remained constant, the retardance of the second layer (δ2) changed from 0 to 1 rad, while the angle of intersection between the two extraordinary axes took different values in the range from 0 deg to 90 deg. In Figs. 2(a) and 2(b), δ1=0.25  rad, and in Figs. 2(c) and 2(d), δ1=0.5  rad. The quantitative variation of the retardance and equivalent extraordinary axis from the BBS was independent of the retardance produced by the first layer.

As shown in Figs. 2(a) and 2(c), if the extraordinary axes of the two anisotropic layers were parallel, the total retardance δ was equal to the sum of δ1 and δ2. The contributions from the two layers to the total retardance were additive when the angle of intersection between the two extraordinary axes was smaller than 45 deg, but cancelled each other out when the angle of intersection exceeded 45 deg. If the extraordinary axes of the two layers were perpendicular to each other, or if θ=90  deg, the anisotropy due to the two layers cancelled out exactly at δ2=δ1, resulting in an isotropic bilayer sample. When δ2 was significantly larger than δ1, the slope of the total retardance curve approached 1, indicating that layer 2 provided the dominant contribution to the total retardance δ and to the anisotropic properties of the bilayer sample.

Figures 2(b) and 2(d) show the orientation of the equivalent extraordinary axis of the bilayer sample as a function of δ2. The equivalent extraordinary axis of the bilayer sample rotated gradually from the direction of the first layer axis (when δ2 was much smaller than δ1) to the direction of the second layer axis (when δ2 became much larger than δ1). When the extraordinary axes of the two layers were perpendicular to each other, the equivalent extraordinary axis of the bilayer sample jumped abruptly from the direction of layer 1 (for δ2<δ1) to the direction of layer 2 (for δ2>δ1). The orientation of the extraordinary axis was undefined at δ2=δ1, because the anisotropy of the two layers cancelled each other out completely when the angle of intersection was 90 deg. In other words, under such a condition, the bilayer medium became isotropical.

Following the analysis in the previous section, we further applied Monte Carlo simulations and Mueller matrix decomposition to two types of bilayer anisotropic samples. Figure 3(a) shows a sample with a layer of birefringent media and a layer of well-aligned cylindrical scatterers [Birefringence and cylindrical scatterer sample (BCS)], and Fig. 3(b) shows a sample with two layers of cylindrical scatterers (CCS).

Fig. 3

(a) Cylinder and birefringence medium model and (b) bilayer cylindrical scatterer model.

JBO_21_5_055002_f003.png

It was shown in our previous study that well-aligned cylindrical scatterers contribute to the retardance.4 The extraordinary axis of the scattering-induced birefringence is the extraordinary axis and is along the direction of the well-aligned cylinders. The magnitude of the birefringence is related to the density, size, and order of alignment of the cylindrical scatterers. Next, we investigated whether the characteristic features in Fig. 2 still applied to the total retardance and equivalent extraordinary axis of a bilayer anisotropic scattering sample.

3.2.

Birefringence and Cylindrical Scatterer Sample

For the BCS, we varied the angles between the orientation of the well-aligned cylinders and the extraordinary axis of the birefringent layer from 0 deg to 90 deg. In the simulation, the diameter of the cylinders was 1  μm, the scattering coefficient of the aligned cylinder layer was μs=15  cm1, and the thickness of the sample was 1 cm. In Fig. 4(a), the vertical axis represents the total retardance of the BCS and the horizontal axis is the retardance produced by the birefringent layer. The different shapes in Fig. 4(a) represent different angles between the aligned cylinders and the extraordinary axis of birefringent media, and the solid lines represent the calculated results for the BBS. Here, the nonzero starting value of 0.38 rad corresponded to the retardance produced by the well-aligned cylindrical scatterers. Figure 4 shows that there was good agreement between the simulated results for BCS and the calculated curves for BBS. The equivalent extraordinary axis of BCS also rotated gradually from the direction of aligned cylinders to the direction of the extraordinary axis of the birefringent layer as the retardance of the birefringent layer increased. The very similar anisotropic behaviors observed for the BCS and BBS confirmed that the well-aligned cylindrical scatterers generated birefringence, whose extraordinary axis was along the direction of the aligned cylinders.

Fig. 4

Results from simulations of the total retardance of a cylindrical scatterer and birefringent media system. The different shapes show data for different angles between the axial direction of the cylinders and the extraordinary axis of birefringent media; the same shapes represent the relationship between the retardance of the media and (a) the total retardance and (b) the equivalent extraordinary axis of the system. The solid line represents the equivalent retardance and extraordinary axis of the BBS system when the two extraordinary axis angles are different.

JBO_21_5_055002_f004.png

3.3.

Two Layers of Cylindrical Scatterers Sample

Since many real fibrous tissues have multilayer microstructures with layers aligned in different directions, we considered further samples consisting of two layers of cylindrical scatterers with different orientations and alignments, as shown in Fig. 3(b). It is known from a previous study4 that the birefringence induced by aligned cylindrical scatterers is affected by the diameter, the optical properties, and the alignment of the cylinders, as well as the scattering coefficient and thickness of the layer. In the simulations, we kept the parameters of the first layer constant and changed the scattering coefficient and the other parameters of the second layer to examine in detail how the scattering-induced anisotropy of the two layers affected the total retardance and the equivalent extraordinary axis of the sample. In Figs. 5(a) and 5(b), the cylinders were perfectly aligned without any fluctuations in the orientation. The horizontal axis represents the retardance of the second layer of cylindrical scatterers obtained by decomposing the simulated Mueller matrices. For Figs. 5(c) and 5(d), the orientations of the cylindrical scatterers followed a Gaussian distribution, whose full width at half maximum (FWHM) was η=20  deg. In the simulations, the diameter of the cylinder was 1  μm, the first layer cylinder scattering coefficient was μs=15  cm1, the second layer cylinder scattering coefficient changed from μs=0  cm1 to μs=40  cm1, and the thickness of the two layers was 1 cm, corresponding to a change in retardance from 0 to 1.14 rad. The different shapes show data for different intersection angles between the directions of the two layers of cylindrical scatterers. The solid lines represent the calculated results for BBS. We also investigated a mixed model with two types of cylindrical scatterers and different alignment orientations in the two layers. The simulation results in Fig. 5 were very similar to those shown Figs. 4(a) and 4(b). This good agreement confirmed that the anisotropic optical features from the coupling of the two layers of cylindrical scatterers with different orientations were approximately equivalent to the features observed for the combination of two pieces of birefringent media with different extraordinary axes.

Fig. 5

Results of simulations of the total retardance in a bilayer cylinder scatterers system. The different shapes show data for different angles between the long axis of cylinder and the extraordinary axis of the birefringent media. For (a) and (b), the two layers of cylinder scatterers were well aligned, without fluctuations. For (c) and (d), the cylinder scatterers’ distribution function (FWHM) η was 20 deg.

JBO_21_5_055002_f005.png

We also studied the variation of the depolarization in the bilayer model; the first layer cylinder scattering coefficient was μs=15  cm1, and the second layer cylinder scattering coefficients were μs=20  cm1 and μs=40  cm1. In Fig. 6, the horizontal axis is the angle between the two layers of cylinders, and the vertical axis is the variation of the retardance and depolarization in this model. Similar to our previous study,4 an increased disorder in the orientation of the cylinders resulted in smaller retardance and larger depolarization, which could explain the decreasing retardance and increasing depolarization with increases in the angle between the cylinder scatterers.

Fig. 6

Results of simulations of a bilayer cylinder model. The cylinder parameter was kept constant, and the angle between the two layers of cylinders was changed. (a) The retardance and (b) depolarization values are shown as a function of the intersection angle.

JBO_21_5_055002_f006.png

3.4.

Experiment Results and Analysis of Bilayer Glass Fibers

We made a bilayer fibrous phantom to validate the simulated and calculated results experimentally. In this experiment, we used two slabs of well-aligned glass fibers as the cylindrical scatterers. The two slabs were placed together with four different intersection angles; these angles were 0 deg, 30 deg, 45 deg, and 90 deg, as shown in Fig. 7. The cylindrical scatterers were made from nonbirefringent glass fibers with a refractive index of 1.547 and a radius of 5  μm. The thickness of each glass fiber layer was 1.5  mm, the scattering coefficient was approximately μs=50  cm1, and the absorption coefficient μa was 0. The experimental setup for the Mueller matrix detection was a typical forward measurement configuration. The details of the experimental setup can be found in the literature.4 The mean experimentally measured retardance values corresponding to the intersection angles of 0 deg, 30 deg, 45 deg, and 90 deg are shown in Fig. 8; the results were consistent with the results of the theoretical calculations and the simulation results as shown in Fig. 6(a).

Fig. 7

Two layers of cylinder scatterers with different angles of (a) 0 deg, (b) 30 deg, (c) 45 deg, and (d) 90 deg.

JBO_21_5_055002_f007.png

Fig. 8

Experimental results for two layers of cylindrical scatterers with different angles of 0 deg, 30 deg, 45 deg, and 90 deg.

JBO_21_5_055002_f008.png

4.

Conclusion

Our previous research showed that there are two sources of anisotropy in tissue: one is cylindrical scatterers and the other is birefringent interstitial media. In this study, we investigated the coupling between two-layer anisotropic samples with different extraordinary axis orientations. We applied Monte Carlo simulations and a Mueller matrix decomposition method for three types of bilayer anisotropic samples, and examined how the total retardance and equivalent extraordinary axis varied as a function of the rotation of one of the layers and the consequent variation in its extraordinary axis and retardance values. The calculations for the two anisotropic layers showed that the total retardance was not a linear superposition of the retardance of the two layers. When the birefringence of the first layer was fixed, and the birefringence of the second layer was increased from zero to well above that of the first layer, the equivalent extraordinary axis of the bilayer sample rotated gradually from the first layer to the second layer. During the same process, the total retardance always increased if the intersection angle between the extraordinary axes of the two layers was smaller than 45 deg, but first decreased and then increased when the intersection angle was larger than 45 deg. When the axes of the two layers were orthogonal, the total retardance became zero when the retardance of the two layers was equal. When one or both of the birefringent layers were replaced by well-aligned cylindrical scatters, Monte Carlo simulation results for the bilayer samples displayed the same characteristic features as those shown by the BBS in terms of both the total retardance and the equivalent extraordinary axis. These results confirmed that the well-aligned cylindrical scatterers generated retardance and behaved the same as a layer of birefringent media. A preliminary experiment using glass fibers as cylindrical scatterers was conducted and the results confirm the above conclusions. This study helped to reveal the connection between observed anisotropy and microstructures and provided more evidence for the feasibility of the quantitative characterization of complex tissues.

Acknowledgments

This work has been supported by the National Natural Science Foundation of China, Grants Nos. 11174178, 11374179, 61205199, and 41475125.

References

1. 

B. D. Cameron, Y. Li and A. Nezhuvingal, “Determination of optical scattering properties in turbid media using Mueller matrix imaging,” J. Biomed. Opt. 11(5), 054031 (2006).JBOPFO1083-3668http://dx.doi.org/10.1117/1.2363347Google Scholar

2. 

A. Kienle, F. K. Forster and R. Hibst, “Anisotropy of light propagation in biological tissue,” Opt. Lett. 29(22), 2617–2619 (2004).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.29.002617Google Scholar

3. 

A. Kienle and R. Hibst, “Light guiding in biological tissue due to scattering,” Phys. Rev. Lett. 97(1), 018104 (2006).PRLTAO0031-9007http://dx.doi.org/10.1103/PhysRevLett.97.018104Google Scholar

4. 

Y. Guo et al., “A study on forward scattering Mueller matrix decomposition in anisotropic medium,” Opt. Express 21(15), 18361–18370 (2013).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.21.018361Google Scholar

5. 

E. Du et al., “Two-dimensional backscattering Mueller matrix of sphere-cylinder birefringence media,” J. Biomed. Opt. 17(12), 126016 (2012).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.17.12.126016Google Scholar

6. 

H. He et al., “Application of sphere-cylinder scattering model to skeletal muscle,” Opt. Express 18(14), 15104–15112 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.18.015104Google Scholar

7. 

Y. Guo et al., “Study on retardance due to well-ordered birefringent cylinders in anisotropic scattering media,” J. Biomed. Opt. 19(6), 065001 (2014).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.19.6.065001Google Scholar

8. 

N. Ortega-Quijano, F. Fanjul-Vélez and J. L. Arce-Diego, “Polarimetric study of birefringent turbid media with three-dimensional optic axis orientation,” Biomed. Opt. Express 5(1), 287–292 (2014).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.5.000287Google Scholar

9. 

P. Ghassemi et al., “A new approach for optical assessment of directional anisotropy in turbid media,” J. Biophotonics 9(1–2), 100–108 (2016).http://dx.doi.org/10.1002/jbio.201400124Google Scholar

10. 

T. W. Gilbert et al., “Collagen fiber alignment and biaxial mechanical behavior of porcine urinary bladder derived extracellular matrix,” Biomaterials 29(36), 4775–4782 (2008).BIMADU0142-9612http://dx.doi.org/10.1016/j.biomaterials.2008.08.022Google Scholar

11. 

K. D. Costa, E. J. Lee and J. W. Holmes, “Creating alignment and anisotropy in engineered heart tissue: role of boundary conditions in a model three-dimensional culture system,” Tissue Eng. 9(4), 567–577 (2003).1937-3341http://dx.doi.org/10.1089/107632703768247278Google Scholar

12. 

S. Alali, Y. Wang and I. A. Vitkin, “Detecting axial heterogeneity of birefringence in layered turbid media using polarized light imaging,” Biomed. Opt. Express 3(12), 3251–3263 (2013).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.3.003250Google Scholar

13. 

T. Yun et al., “Monte Carlo simulation of polarized photon scattering in anisotropic media,” Opt. Express 17(19), 16590–16602 (2009).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.17.016590Google Scholar

14. 

S. Y. Lu and R. A. Chipman, “Interpretation of Mueller matrix based on polar decomposition,” J. Opt. Soc. Am. A 13(5), 1106–1113 (1996).JOAOD60740-3232http://dx.doi.org/10.1364/JOSAA.13.001106Google Scholar

15. 

S. Manhas et al., “Mueller matrix approach for determination of optical rotation in chiral turbid media in backscattering geometry,” Opt. Express 14(1), 190–202 (2006).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.14.000190Google Scholar

16. 

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

Biographies for the authors are not available.

Yihong Guo, Nan Zeng, Honghui He, Celong Liu, E. Du, Yonghong He, Hui Ma, "Retardance of bilayer anisotropic samples consisting of well-aligned cylindrical scatterers and birefringent media," Journal of Biomedical Optics 21(5), 055002 (9 May 2016). http://dx.doi.org/10.1117/1.JBO.21.5.055002
JOURNAL ARTICLE
6 PAGES


SHARE
KEYWORDS
Scattering

Monte Carlo methods

Birefringence

Tissues

Anisotropy

Optical fibers

Glasses

Back to Top