Degree of polarization surfaces and maps for analysis of depolarization

doi:10.1364/OPEX.12.004941

« Show journal navigation

Degree of polarization surfaces and maps for analysis of depolarization

Brian DeBoo, J. Sasian, and R. Chipman »View Author Affiliations

Optics Express, Vol. 12, Issue 20, pp. 4941-4958 (2004)
http://dx.doi.org/10.1364/OPEX.12.004941

View Full Text Article

Acrobat PDF (1938 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (30)
Cited By
Figures (16)
Multimedia (4)
Metrics

Abstract

The concept of degree of polarization surfaces is introduced as an aid to classifying the depolarization properties of Mueller matrices. Degree of polarization surfaces provide a visualization of the dependence of depolarization on incident polarization state. The surfaces result from a non-uniform contraction of the Poincaré sphere corresponding to the depolarization properties encoded in a Mueller matrix. For a given Mueller matrix, the degree of polarization surface is defined by moving each point on the unit Poincaré sphere radially inward until its distance from the origin equals the output state degree of polarization for the corresponding input state. Of the sixteen elements in a Mueller matrix, twelve contribute to the shape of the degree of polarization surface, yielding a complex family of surfaces. The surface shapes associated with the numerator and denominator of the degree of polarization function are analyzed separately. Protrusion of the numerator surface through the denominator surface at any point indicates non-physical Mueller matrices. Degree of polarization maps are plots of the degree of polarization on flat projections of the sphere. These maps reveal depolarization patterns in a manner well suited for quantifying the degree of polarization variations, making degree of polarization surfaces and maps valuable tools for categorizing and classifying the depolarization properties of Mueller matrices.

1. Inroduction

The transformation from an incident polarization state into an exiting polarization state which occurs during a linear interaction between light and matter (e.g. transmission through a retarder or polarizer, or reflection from a thin film) is commonly described using one of two polarization calculi, the Jones calculus or the Mueller calculus. The Jones calculus, which contains absolute amplitude and phase information, is more useful for describing fully polarized light and non-depolarizing optical devices [1–3

R.C. Jones, “A new calculus for the treatment of optical systems: I. Description and discussion of the calculus,” J. Opt. Soc. Am. 31, 488–493 (1941). [CrossRef]

]. All Jones matrices are also expressible as Mueller matrices and these non-depolarizing Mueller matrices are well understood [4

K.D. Abhyankar and A.L. Fymat, “Relations between the elements of the phase matrix for scattering,” J. Math Phys. 10, 1935–1938 (1969). [CrossRef]

, 5

E.S. Fry and G.W. Kattawar, “Relationships between elements of the Stokes matrix,” Appl. Opt. 20, 2811–2814 (1981). [CrossRef] [PubMed]

]. The Mueller calculus is a more general representation for polarization interactions, which applies to incoherent states and describes polarized, partially polarized, or unpolarized light and can quantify depolarization. Among its other properties, a Mueller matrix describes how incident polarized light is depolarized for any incident state and describes the variations of depolarization with polarization state which occur. A single Jones matrix cannot describe depolarization. Sets of Jones matrices can describe depolarization (e.g., Jones matrices as a function of wavelength or position) but such Jones matrices cannot be measured in a straightforward and convenient manner to characterize the depolarization of an optical element, liquid crystal cell, scattering surface, or other depolarizing sample [6

W. Swindell, ed., Benchmark Papers in Optics: Polarized Light , (Hutchinson and Ross, Stroudsburg, 1975).

, 7

R.A. Chipman, “Polarimetry,” in Handbook of Optics , Vol. II, (McGraw Hill, New York, 1995).

]. Thus Mueller matrices are almost always used to characterize depolarization

Depolarization is the reduction of the degree of polarization of light. In the Mueller calculus depolarization can be pictured as a coupling of polarized into unpolarized light, where polarized light is incident and the exiting Stokes vector can be mathematically separated into a fully polarized and an unpolarized Stokes vector. Depolarization of some optical devices has been described (e.g. liquid crystals [8

J.L. Pezzaniti, S.C. McClain, R.A. Chipman, and S.Y. Lu, “Depolarization in liquid crystal TV’s,” Opt. Lett. 18, 2071–2073 (1993). [CrossRef] [PubMed]

]), and commercial depolarizers are available. However, reports of depolarization measurements in the literature (see, e.g., [9–20

S.-M.F. Nee and T.-W. Nee, “Principal Mueller matrix of reflection and scattering for a one-dimensional rough surface,” Opt. Eng. 41, 994–1001 (2002). [CrossRef]

]) have been relatively limited. This limited number of published depolarization measurements may in part be due to the only recent commercialization of Mueller matrix polarimeters, which offer the most straightforward way to measure depolarization. Lenses, mirrors, filters, and other typical optical elements exhibit very small amounts of depolarization, typically less than a few tenths of a percent. In contrast, the depolarization of most diffusely reflecting objects such as paints, metal and wood surfaces, natural materials, etc, is significant, varying from a few percent to 100% (i.e. complete depolarization).

A single-valued depolarization metric, the depolarization index, has been introduced to describe the degree to which a Mueller matrix depolarizes incident states [7

R.A. Chipman, “Polarimetry,” in Handbook of Optics , Vol. II, (McGraw Hill, New York, 1995).

, 21

J.J. Gil and E. Bernabeau, “Depolarization and polarization indices of an optical system,” Opt. Acta 33, 185–189 (1986). [CrossRef]

]. However, such a single number metric cannot describe the full complexity of depolarization associated with a Mueller matrix.

This capability of Mueller matrices to define depolarization for all incident states is useful, particularly when strong depolarization occurs [22

B. DeBoo, “Investigation of polarization scatter properties using active imaging polarimetry,” PhD dissertation, Optical Sciences Center, University of Arizona (2004).

]. Depolarization is associated with a reduction in the degree of polarization ( DoP ) of incident states. Here, DoP surfaces and maps are introduced to quantify and visualize the complete depolarizing properties of a Mueller matrix for all incident states.

Practically, DoP is a measure of the randomness of polarization in a light beam, a property characterized by how much of this beam may be blocked by a polarizer. Mathematically, on the Poincaré sphere, the DoP represents the distance of a normalized Stokes vector’s last three components from the origin. The surface of the unit Poincaré sphere has DoP =1 and represents all fully polarized states [23

D. Goldstein, Polarized Light, Second Ed., Revised and Expanded , (Marcel Dekker, New York, 2003). [CrossRef]

]. A depolarizing interaction causes fully polarized Stokes states on the surface of the Poincaré sphere to emerge with DoP ≤1. The DoP surface plots these exiting DoP values along the same radial vectors from the origin that define the corresponding input states. A similar plot was discussed in a different context by Williams [10

M.W. Williams, “Depolarization and cross polarization in ellipsometry of rough surfaces.” Appl. Opt. 25, 3616–3621 (1986). [CrossRef] [PubMed]

]. A DoP map is a two-dimensional contour plot of the DoP surface versus two coordinates defining input polarization states (i.e. polarization ellipse orientation and ellipticity).

This paper investigates the way that a general depolarizing Mueller matrix affects DoP surface and map shapes, thus providing new methods for classifying Mueller matrices and their behaviors. The DoP surfaces are explored by examining separately the changes imparted to the surface shape by the diattenuation, retardance, polarizance, and depolarization properties encoded in a Mueller matrix. DoP maps and surfaces are introduced in Section 2. The denominator of the DoP function, which contains diattenuation properties of the Mueller matrix, is treated in Section 3. The retardance and part of the depolarizing properties of the Mueller matrix are covered in Section 4, via a singular valued decomposition, while Section 5 addresses the additional depolarization effect of polarizance on the DoP surface. Section 6 presents an analysis of complete DoP surfaces and maps. Section 7 addresses the physical realizability of a Mueller matrix using a criterion related to the DoP numerator and denominator. Section 8 explores a family of Mueller matrices and the evolution of their DoP surfaces and maps to better understand how features of each plot relate to specific portions of the Mueller matrix.

Describing depolarization using the DoP of polarized states offers significant advantages over single-valued methods. DoP surfaces and maps provide an insightful visualization of how a Mueller matrix depolarizes all incident polarized states. The information content of DoP surfaces and maps is greater than other depolarization metrics. For example, maxima and minima in the output DoP are readily observed with DoP surfaces, and the values of these extremes easily quantified with DoP maps. DoP surfaces and maps are thorough in describing the behaviors of depolarizing Mueller matrices and are examined in detail here.

2. The degree of polarization

The degree of polarization ( DoP ) characterizes the randomness of a polarization state. The DoP of Stokes vector S=(S₀,S₁,S₂,S₃) is

DoP (S) = \frac{\sqrt{S_{1}^{2} + S_{2}^{2} + S_{3}^{2}}}{S_{0}} .

(1)

When DoP =0, the light is unpolarized and all ideal polarizers block half the beam. When DoP =1, the beam is completely polarized and some ideal polarizer, either linear, circular, or elliptical, will completely block the beam. Thus, (1+ DoP )/2 is the fraction of a beam that can be blocked by an ideal polarizer.

Normalized Stokes vectors with DoP = 1 are fully polarized states which lie on the surface of the unit Poincaré sphere. The Stokes vectors on the surface of the Poincaré sphere can be parameterized as

S (θ, φ) = (\begin{matrix} 1 \\ Cos (2 θ) Cos (φ) \\ Sin (2 θ) Cos (φ) \\ Sin (φ) \end{matrix}),

(2)

where θ is the polarization orientation (one half the longitude on a globe), and φ is the latitude on the Poincaré sphere. For example, (θ, φ)=(0°, 0°) is horizontal linearly polarized light, (45°, 0°) is 45° polarized light, and (θ, 90°) represents right circularly polarized light for all θ. The Degree of Circular Polarization ( DoCP ) of a Stokes vector is defined as DoCP = S₃/S₀ = Sin(φ). A circular polarizer will block (1+| DoCP |)/2 of a beam; the polarizer’s helicity is right or left depending on whether the state is in the upper or lower hemisphere of the Poincaré sphere, respectively.

3. Degree of polarization surfaces and maps

The DoP surface for a Mueller matrix, M, is formed by moving normalized Stokes vectors, S, on the surface of the Poincaré sphere radially inward to a distance DoP (S′=M-S) from the origin, plotted for all incident S on the surface of the Poincaré sphere, given in Eq. (2). The DoP surface results from the product of a scalar, the DoP , and a vector, (S₁, S₂, S₃), formed from the last three elements of the normalized Stokes vector,

DoPSurface (M, S) = \frac{\sqrt{{S_{1}}^{'} {(M, S)}^{2} + {S_{2}}^{'} {(M, S)}^{2} + {S_{3}}^{'} {(M, S)}^{2}}}{{S_{0}}^{'} (M, S)} (S_{1}, S_{2}, S_{3}),

(3)

for all (

S_{1}^{2}

S_{2}^{2}

S_{3}^{2}

)^1/2=1.

The DoP map for a Mueller matrix is a contour plot of the DoP of exiting light as a function of the incident polarized state and represents a “flattened” DoP surface. In this paper, the DoP map is plotted with axes θ (polarization ellipse major axis orientation) and DoCP, but there is some flexibility in the choice of parameterization of the polarized Stokes vectors. In general the DoP map provides easier visualization of maxima, minima, saddles, and other features of the depolarization variation than the DoP surface.

Figure 1 shows a DoP surface and its corresponding DoP map for an example Mueller matrix with depolarization,

M_{1} = (\begin{matrix} 1.0 & - 0.226 & 0.069 & 0.196 \\ - 0.030 & 0.052 & 0.357 & - 0.336 \\ 0.069 & - 0.454 & - 0.266 & - 0.194 \\ 0.196 & - 0.336 & 0.194 & 0.584 \end{matrix}) .

(4)

Fig. 1. Degree of polarization surface (a) and degree of polarization map (b) of a typical depolarizing Mueller matrix showing two local maxima, two local minima, and two saddles.

The distance from the origin to the DoP surface is the output DoP for the corresponding incident state on the Poincaré sphere. Remember that the output polarization state is in general different from the input polarization state, and this output state information is not contained in the DoP Surface. Where the surface is pinched toward the origin, those incident Stokes states are more depolarized by the Mueller matrix. In this example, (160°, 0.1) is depolarized the most and (175°, -0.7) is depolarized the least.

The DoP surface is analyzed here by decomposing the surface into terms relating to specific polarization properties of the Mueller matrix. The output Stokes vectors for M as a function of the input Stokes vector, S, are

(\begin{matrix} {S_{0}}^{'} \\ {S_{1}}^{'} \\ {S_{2}}^{'} \\ {S_{3}}^{'} \end{matrix}) = (\begin{matrix} m_{00} & m_{01} & m_{02} & m_{03} \\ m_{10} & m_{11} & m_{12} & m_{23} \\ m_{20} & m_{21} & m_{22} & m_{23} \\ m_{30} & m_{31} & m_{32} & m_{33} \end{matrix}) \cdot (\begin{matrix} S_{0} \\ S_{1} \\ S_{2} \\ S_{3} \end{matrix}) = (\begin{matrix} m_{00} S_{0} + m_{01} S_{1} + m_{02} S_{2} + m_{03} S_{3} \\ m_{10} S_{0} + m_{11} S_{1} + m_{12} S_{2} + m_{13} S_{3} \\ m_{20} S_{0} + m_{21} S_{1} + m_{22} S_{2} + m_{23} S_{3} \\ m_{30} S_{0} + m_{31} S_{1} + m_{32} S_{2} + m_{33} S_{3} \end{matrix}) .

(5)

Substituting Eq. (5) into Eq. (3), the scalar part of the DoP surface takes the rather cumbersome form

DoP Surface (M, S) = \frac{DoP Numerator}{DoP Deno \min ator} = \frac{{(\sum_{k = 1}^{3} {(m_{k 0} S_{0} + m_{k 1} S_{1} + m_{k 2} S_{2} + m_{k 3} S_{3})}^{2})}^{\frac{1}{2}}}{m_{00} S_{0} + m_{01} S_{1} + m_{02} S_{2} + m_{03} S_{3}} .

(6)

4. The denominator of the DoP surface

Comparing Eqs. (3) and (6) shows that the denominator of the DoP surface function is

S_{0}^{'} (M, S) = m_{00} S_{0} + m_{01} S_{1} + m_{02} S_{2} + m_{03} S_{3},

(7)

the output flux, S ₀′, for input S. Note that S ₀′ depends only on the first row of M. The m₀₀ element is the intensity throughput of the Mueller matrix for unpolarized light and describes losses associated with absorption, reflection or scattering. The elements m₀₁, m₀₂, and m₀₃ characterize the diattenuation, the tendency of M to act as a partial polarizer [23

D. Goldstein, Polarized Light, Second Ed., Revised and Expanded , (Marcel Dekker, New York, 2003). [CrossRef]

The shape of S ₀′(M,S), the denominator surface, is a limaÇon of revolution [24

R.C. Yates, A Handbook on Curves and Their Properties , (J. W. Edwards, Ann Arbor, 1952).

DoPDenominator (M, S) = m_{00} + m_{01} Cos (2 θ) Cos (φ) + m_{02} Sin (2 θ) Cos (φ) + m_{03} Sin (φ) .

(8)

The simplest polar form of a limaÇon is

r = 1 + ξ Cos (α), (0 \leq α < 2 π),

(9)

plotted in Fig. 2 as ξ varies. LimaÇons belong to a more general family of curves, the botanic curves, with polar form

r = 1 + ξ Cos (cα), (0 \leq α < 2 π),

(10)

where “c” describes the number of “petals” of the botanic curve. Limaçons are botanic curves of order c=1 with only a single petal.

Fig. 2. Family of limaçons. DoP denominator surfaces are limaçons rotated about their axis of symmetry.

Consider the depolarizing Mueller matrix M ₀ constructed from the sum of horizontal linear polarizer and 45° linear polarizer Mueller matrices,

M_{0} = \frac{1}{2} (\begin{matrix} 1 & 1 & 0 & 0 \\ 1 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{matrix}) + \frac{1}{2} (\begin{matrix} 1 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 \\ 1 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 \end{matrix}) = (\begin{matrix} 1 & 0.5 & 0.5 & 0 \\ 0.5 & 0.5 & 0 & 0 \\ 0.5 & 0 & 0.5 & 0 \\ 0 & 0 & 0 & 0 \end{matrix}),

(11)

which physically corresponds to covering an aperture with equal fractions of the two polarizers. The denominator of the DoP surface for M₀ is

DoP Denominator (M_{0}, S) = 1 + 0.5 Cos (2 θ) Cos (φ) + 0.5 Sin (2 θ) Cos (φ),

(12)

which represents the intensity transmittance for any input Stokes vector, and is plotted in Fig. 3. Note the dimple in the limaçon shape corresponding to minimum transmittance, which occurs along the axis of rotational symmetry of the denominator surface. This axis of symmetry defines the diattenuation axis through the Poincaré sphere.

Fig. 3. DoP denominator surface for Mueller matrix M₀ of Eq. (11). The red line denotes the diattenuation axis, the axis of rotational symmetry of the surface.

The Mueller matrix elements m₀₁, m₀₂, and m₀₃ define three degrees of freedom, two of which determine the diattenuation axis orientation and latitude as

({(2 θ)}_{axis}, ϕ_{axis}) = (ArcCos [\frac{m_{01}}{\sqrt{{m_{01}}^{2} + {m_{02}}^{2}}}], ArcSin [\frac{m_{03}}{\sqrt{{m_{01}}^{2} + {m_{02}}^{2} + {m_{03}}^{2}}}]) .

(13)

The diattenuation axis passes through the maximum and the minimum (dimple) of the DoP denominator surface, as well as the Poincaré sphere points yielding these transmittance extema, T_max=((2θ)_axis, ϕ_axis) and T_min= (π+(2θ)_axis, -ϕ_axis), respectively. The third degree of freedom defined by Mueller elements m₀₁, m₀₂, and m₀₃ describes the diattenuation (also “diattenuation magnitude”), d, as

d = \frac{\sqrt{{m_{01}}^{2} + {m_{02}}^{2} + {m_{03}}^{2}}}{m_{00}} = \frac{T_{max} - T_{min}}{T_{max} + T_{min}},

(14)

which describes the degree to which M is a partial polarizer, given by the maximum and minimum transmittance [23

D. Goldstein, Polarized Light, Second Ed., Revised and Expanded , (Marcel Dekker, New York, 2003). [CrossRef]

]. For ideal polarizers (d =1), the dimple in the DoP denominator surface becomes a cusp with its point at the origin.

5. The degree of polarization numerator

The DoP numerator surface, described in Eq. (6), is more complex than the denominator surface because it depends on S₁′,S₂′, and S₃′ and includes a square root of the sum of three terms. The DoP numerator depends on the twelve elements in the bottom three rows of the Mueller matrix, which encode the polarizance, retardance, and depolarization properties of the matrix [7

R.A. Chipman, “Polarimetry,” in Handbook of Optics , Vol. II, (McGraw Hill, New York, 1995).

, 23

D. Goldstein, Polarized Light, Second Ed., Revised and Expanded , (Marcel Dekker, New York, 2003). [CrossRef]

Consider the last three Mueller matrix rows separated into the sum of two matrices:

(\begin{matrix} \dots & \dots & \dots & \dots \\ m_{10} & m_{11} & m_{12} & m_{13} \\ m_{20} & m_{21} & m_{22} & m_{23} \\ m_{30} & m_{31} & m_{32} & m_{33} \end{matrix}) = (\begin{matrix} \dots & \dots & \dots & \dots \\ m_{10} & 0 & 0 & 0 \\ m_{20} & 0 & 0 & 0 \\ m_{30} & 0 & 0 & 0 \end{matrix}) + (\begin{matrix} \dots & \dots & \dots & \dots \\ 0 & m_{11} & m_{12} & m_{13} \\ 0 & m_{21} & m_{22} & m_{23} \\ 0 & m_{31} & m_{32} & m_{33} \end{matrix}) .

(15)

The first column of M is the polarizance vector, P, the output Stokes vector when unpolarized light is incident,

P (M) = M \cdot (\begin{matrix} 1 \\ 0 \\ 0 \\ 0 \end{matrix}) = (\begin{matrix} m_{00} \\ m_{10} \\ m_{20} \\ m_{30} \end{matrix}) .

(16)

The lower right 3×3 submatrix of M,

Q = (\begin{matrix} m_{11} & m_{12} & m_{13} \\ m_{21} & m_{22} & m_{23} \\ m_{31} & m_{32} & m_{33} \end{matrix}),

(17)

is an orthogonal matrix (real, unitary [25

G.B. Arfken and H.J. Weber, Mathematical Methods for Physicists , (Academic Press, San Diego, 2001).

]) for retarders. Diattenuators, as well as depolarizers (uniform and non-uniform), generate Hermitian [25

G.B. Arfken and H.J. Weber, Mathematical Methods for Physicists , (Academic Press, San Diego, 2001).

] Q.

Submatrix Q can be decomposed into the product of a unitary matrix and a Hermitian matrix in two ways [26

S.-Y. Lu, “An interpretation of polarization matrices,” PhD dissertation, Dept. of Physics, University of Alabama at Huntsville (1995).

, 27

S.-Y. Lu and R.A. Chipman, “Interpretation of Mueller matrices based on polar decomposition,” J. Opt. Soc. Am. A 13, 1106–1113 (1996). [CrossRef]

Q = HU = U' H' .

(18)

Unitary matrices represent rotations and preserve magnitude (i.e. impart no contraction or expansion to a surface). Retarders rotate the Poincaré sphere about the fast and slow axis of the retarder, so Mueller matrices and Q matrices for retarders are unitary.

The Hermitian component, H or H’, of submatrix Q (see Eq. (18)) is diagonalizable into the form of a similarity transformation [25

G.B. Arfken and H.J. Weber, Mathematical Methods for Physicists , (Academic Press, San Diego, 2001).

]

H = U_{1} D {U_{1}}^{- 1} or H' = {U_{2}}^{- 1} D U_{2},

(19)

where U₁ and U₂ are unitary and D is a diagonal matrix of H’s eigenvalues. The three diagonal elements of D operate on a sphere to contract the sphere along three orthogonal coordinate axes. The amount of contraction, 1- DoP (S′), is related to the amount of depolarization. Consider the three cases D₁ , with all diagonal elements equal, D₂ , with two diagonal elements equal, and D₃ , with all three diagonal elements different,

D_{1} = (\begin{matrix} a & 0 & 0 \\ 0 & a & 0 \\ 0 & 0 & a \end{matrix}); D_{2} = (\begin{matrix} a & 0 & 0 \\ 0 & a & 0 \\ 0 & 0 & b \end{matrix}); D_{3} = (\begin{matrix} a & 0 & 0 \\ 0 & b & 0 \\ 0 & 0 & c \end{matrix}) .

(20)

D₁ shrinks the Poincaré sphere uniformly into a smaller sphere. The corresponding Mueller matrix, UD, given by

UD = (\begin{matrix} 1 & 0 & 0 & 0 \\ 0 & a & 0 & 0 \\ 0 & 0 & a & 0 \\ 0 & 0 & 0 & a \end{matrix})

(21)

is termed a uniform depolarizer Mueller matrix since UD depolarizes all incident states equally. D₂ and D₃ indicate nonuniform contractions along orthogonal axes associated with nonuniform depolarization. Figures 4–6 show the depolarization-contracted surfaces and cross-sections associated with D₁ , D₂ , and D₃ , respectively. The cross sections all belong to the family of botanic curves of order c=2 (see Eq. (10)), having two petals, where the circle in Fig. 4 is the special case of ξ=0.

Fig. 4. Surface shape and associated cross section after contraction of unit Poincaré sphere by uniform depolarization matrix of type D₁ with {a,a,a}={0.7,0.7,0.7}.

Fig. 5. Surface shape (a) and associated cross sections (b), (c), (d), after contraction of unit Poincaré sphere by nonuniform depolarization matrix of type D₂ with {a,b,a}={0.6,0.2,0.6}.

Fig. 6. Surface shape (a) and associated cross sections (b), (c), (d), after contraction of unit Poincaré sphere by nonuniform depolarization matrix of type D₃ with {a,b,c}={0.4,0.2,0.8}.

The singular value decomposition (SVD) of Q is obtained by substituting Eq. (19) into Eq. (18),

Q = (U_{1} D {U_{1}}^{- 1}) U_{3} = U_{0} ({U_{2}}^{- 1} D U_{2}) = U_{1} D U_{2} .

(22)

The product of two unitary matrices (

{U_{1}}^{- 1}

U₃ or U₀

{U_{2}}^{- 1}

) is unitary just as successive rotations can be replaced by a single rotation. Note that the retardance components U₁ and U₂ of the decomposition are the same regardless of the decomposition order chosen in Eq. (18). The form U₁DU₂ relates the retardance and depolarization to the shape of the DoP numerator surface. Matrix U₂ rotates the Poincaré sphere prior to D contracting it along three axes. This is followed by a final rotation U₁ that rotates polarization states into their final output states. Since U₁ (being unitary) preserves vector magnitude, and the DoP surface is plotted in terms of incident Stokes states, U₁ bears no effect (rotation or otherwise) on the DoP numerator surface. U₁ does however affect a related DoP surface representation plotted in terms of output Stokes states.

In the SVD of M₁ of Eq. (4),

Q_{1} = (\begin{matrix} 0.4224 & 0.2741 & 0.8640 \\ - 0.0474 & - 0.9452 & 0.3230 \\ - 0.9052 & 0.1774 & 0.3863 \end{matrix}) • (\begin{matrix} 0.7477 & 0 & 0 \\ 0 & 0.5769 & 0 \\ 0 & 0 & 0.3981 \end{matrix}) • (\begin{matrix} 0.4652 & - 0.0166 & - 0.8850 \\ 0.6665 & 0.6645 & 0.3379 \\ - 0.5825 & 0.7471 & - 0.3202 \end{matrix})

= U_{1} D U_{2},

(23)

neither U₁ nor U₂ affects the DoP numerator shape, and only U₂ affects the orientation. The contracting effect of D₁ (along axes first rotated by U₂ ) is shown in Fig. 7. The Poincaré sphere is contracted along three orthogonal axes to radii of 0.7477, 0.5769 and 0.3981. The first contraction occurs along an axis specified by the top row of U₂ , and so forth. The resulting surface is not always convex.

Fig. 7. Final shape of surface from interaction of Q₁ with incident fully polarized Stokes vectors.

The SVD handles singular matrices, such as ideal polarizers, without complication (e.g. Q₀ of matrix M₀ in Eq. (11)).

6. Effect of polarizance on DoP numerator

The DoP numerator surface depends on the sum of Q and the polarizance as shown in Eqs. (6) and (15). This combination drags and distorts the Q-surface in the direction of the polarizance vector (Eq. (16)), moving the origin a distance equal to the magnitude of the polarizance vector, P. Before the effect of polarizance, the DoP numerator surface coincides with the Q-surface, and the square root in the DoP numerator represents the magnitudes of the radial vectors to this surface. Because the summation of P with Q occurs inside the DoP numerator square root, the induced effect on the DoP numerator surface is not merely a pure translation, but rather a translation and distortion. This transformation is animated for M₀ (Eq. (11)) in Fig. 8, and for M₁ (Eq. (4)) in Fig. 9, where the final shapes in Figs. 8 and 9 are the DoP numerator surfaces for M₀ and M₁ , respectively.

Fig. 8. (1.33 MB) Animation showing the effect of polarizance on the numerator surface for Mueller matrix M₀.

Fig. 9. (1.33MB) Animation showing the effect of polarizance on the numerator surface for Mueller matrix M₁.

Formation of the DoP numerator surface is visualized in three steps, a rotation by U₂ , the subsequent contraction of the Poincaré sphere in three orthogonal Stokes dimensions by D, and the shift and distortion (associated with the DoP numerator square root) in the direction of the three-component polarizance vector.

7. DoP surface and DoP map

The DoP surface is the quotient of the DoP numerator and DoP denominator surfaces. For example, Fig. 1(a) is the DoP surface for M₁ (Eq. (4)), and Fig. 10 is the DoP surface for M₀ (Eq. (11)).

Fig. 10. DoP surface for Mueller matrix M₀.

Figure 11 contains the DoP map for Mueller matrix M₀ (the DoP map for M₁ is Fig. 1(b)).

Fig. 11. DoP map for Mueller matrix M₀.

DoP maps and surfaces may exhibit one or two maxima and one or two minima. These maxima or minima may be degenerate for entire circles of incident states around the Poincaré sphere, as in the minimum of Fig. 11 for M₀ . Two local maxima also appear in the DoP map for M₀ . M₁ has two minima and two maxima. After an extensive search, cases where the DoP maps contain three or more maxima or minima have not been found.

8. Physical realizability of a Mueller matrix

When the output DoP for any incident state is less than zero or greater than one, that matrix is not a physically realizable Mueller matrix. Many relationships among Mueller matrix elements to ensure physical realizability have been published [28–30

C. Brosseau, Fundamentals of Polarized Light: A Statistical Optics Approach , (John Wiley and Sons, Inc., New York, 1998).

]. Such relationships can be understood geometrically in terms of the numerator and denominator surfaces. The numerator surface must not protrude from the denominator surface at any point or the DoP at that point is greater than one. The two surfaces may be tangent and the DoP is one at points of tangency. Similarly, the origin of the Poincaré sphere must lie within or on both surfaces, or the DoP will be negative. When the numerator surface passes through the origin, the corresponding state is completely depolarized.

Plotting the DoP numerator and denominator surfaces together for a given Mueller matrix aids the visualization of these geometrical relationships and is useful for establishing physical realizability. For example, the numerator and denominator surfaces for M₁ (Eq. (4)) are plotted in Fig. 12. The denominator surface surrounds the numerator surface, so M₁ is physical.

Fig. 12. DoP numerator surface (colored) plotted within DoP denominator surface (wireframe). The numerator surface must lie entirely within the denominator surface for physically realizable Mueller matrices, to maintain a DoP less than or equal to one.

Consider a family of matrices connecting the ideal depolarizer with M₀ parameterized by coefficient i ,

M' = (\begin{matrix} 1 & i & i & 0 \\ i & i & 0 & 0 \\ i & 0 & i & 0 \\ 0 & 0 & 0 & 0 \end{matrix}) .

(24)

As i increases the numerator surface uniformly grows starting as a single point at the origin. Figure 13 animates the growth of the numerator surface within the denominator surface for i ranging from 0.0 to 0.7. When i = 0.5 the numerator surface is tangent to the denominator surface at the two maxima of Fig. 11. When i exceeds 0.5, the numerator surface protrudes from the denominator surface and the Mueller matrices are non-physical.

Fig. 13. (967 KB) DoP numerator surface plotted within DoP denominator surface, as Mueller matrix parameter “i” of Eq. (24) is increased in increments of 0.1 from zero to 0.7.

9. A family of DoP maps

A family of Mueller matrices will be analyzed to demonstrate how certain matrix properties relate to properties of DoP surfaces and maps. Observing the metamorphosis of DoP surfaces and maps clarifies the depolarizing effects of the matrix family. Consider the family of Mueller matrices

M_{j} = LP (0) + LP (j \frac{π}{12}), j ∊ {0,1, \dots 11}

(25)

whose members are formed by summation of a horizontal linear polarizer and another linear polarizer, LP(α), with its polarization axis at an angle α to horizontal. Figure 14 shows the DoP surfaces for each member of this family plotted within the Poincaré sphere, and Fig. 15 the corresponding DoP maps. Figure 16 animates the evolution of DoP surfaces.

This family always has two maxima with DoP = 1 which occur for Stokes vectors orthogonal to the two polarizer axes. An incident state orthogonal to one of the polarizers is completely blocked by that polarizer and light only transmits through the other polarizer, thus emerging with DoP = 1. States aligned with either polarizer also partially transmit through the second polarizer (unless the two polarizers are orthogonal) and the contribution of these two beams reduces the DoP . Since one of the polarizers forming every member of Eq. (25) is a horizontal linear polarizer, every DoP surface and map in this family shows a maximum DoP = 1 for vertically polarized light. The other DoP maximum rotates around the Poincaré sphere as the second polarizer rotates. Half way between these maxima is a valley of increasing width and depth (see Fig. 14). When the second polarizer axis is at 90° to the first, the output DoP maxima occur for input vertically and horizontally polarized light and the valley reaches a DoP = 0. The DoP surface in this case separates into two tangent spheres. The DoP = 0 output states occur for input states on the great circle of the Poincaré sphere through 45°, right circular, 135°, and left circular states.

Fig. 14. Family of DoP surfaces for LP(0)+LP(jπ/12). The DoP surfaces are plotted within the normalized Poincaré sphere (unit radius) to ensure that the DoP remains physical (radius less than one) throughout the family.

Fig. 15. DoP maps for Mueller matrix family of Eq. (25). Note that the map always has two maxima and one mimimum (except for the degenerate case of a sphere). Also note that the scale changes between plots to better convey information using the same number of contours.

Fig. 16. (1.46 MB) Evolution of DoP surfaces through members of Mueller matrix family in Eq. (25).

10. Conclusion

DoP maps and surfaces represent the variation of depolarization of a Mueller matrix for all fully polarized incident states. For non-depolarizing Mueller matrices, the DoP surface is the unit sphere and the DoP map is unity everywhere. For depolarizing Mueller matrices the DoP surface contracts toward the origin by an amount equal to the depolarization for that incident state. Thus, the maps and surfaces indicate the variation of depolarization with incident state.

Twelve of sixteen Mueller matrix degrees of freedom affect the shape of DoP surfaces, with three related to the DoP denominator arising from the diattenuation properties encoded in the first row of the Mueller matrix. Two of three denominator degrees of freedom define a diattenuation axis and the third the diattenuation magnitude which relates to the size of the dimple in the limaçon cross-section of the denominator surface.

The DoP numerator surface function has nine degrees of freedom relating to depolarization, retardance, and polarizance. The unitary matrix U₁ in the singular value decomposition of Mueller submatrix Q is a rotation of the output polarization state due to retardance which does not change the DoP surface. Three depolarization degrees of freedom in matrix D contract the Poincaré sphere along three axes rotated from the Stokes basis states by the retardance effects of unitary matrix U₂ . The polarizance vector acts via an additive term that drags and distorts the surface in the direction of the polarizance vector. DoP surfaces may also be plotted in terms of the output Stokes states rather than input Stokes states (making visible the rotation effects of U₁ ) but this representation is not explored here.

For a physically realizable Mueller matrix, the DoP numerator surface must lie entirely within the denominator surface and must incorporate the Poincaré sphere origin (either on or within the surface). Physically realizable DoP maps have only been observed with one or two maxima and one or two minima. It is postulated that the DoP maps cannot have more than two maxima or minima, but such a conjecture is not proven here.

DoP surfaces and DoP maps yield a detailed picture of depolarization and are of useful for understanding and classifying depolarizing Mueller matrices.

Acknowledgments

The authors wish to thank Sandia National Laboratories for its support of this work. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000. This paper is a section of a dissertation submitted in partial fulfillment of the requirements for a PhD degree in Optical Sciences at the University of Arizona. The authors acknowledge the help of Bridget Ford and Eric Shields. Special thanks go to Justin Wolfe and Neil Beaudry of the Optical Sciences Center, University of Arizona, for their assistance.

References and Links

1.	R.C. Jones, “A new calculus for the treatment of optical systems: I. Description and discussion of the calculus,” J. Opt. Soc. Am. 31, 488–493 (1941). [CrossRef]
2.	R.C. Jones, “New calculus for the treatment of optical systems: VIII. Electromagnetic theory,” J. Opt. Soc. Am. 46, 126–131 (1956). [CrossRef]
3.	M. Izdebski, W. Kucharczyk, and R.E. Raab, “Application of the Jones calculus for a modulated double-refracted light beam propagating in a homogeneous and nondepolarizing electro-optic uniaxial crystal,” J. Opt. Soc. Am. A 21, 132–139 (2004). [CrossRef]
4.	K.D. Abhyankar and A.L. Fymat, “Relations between the elements of the phase matrix for scattering,” J. Math Phys. 10, 1935–1938 (1969). [CrossRef]
5.	E.S. Fry and G.W. Kattawar, “Relationships between elements of the Stokes matrix,” Appl. Opt. 20, 2811–2814 (1981). [CrossRef] [PubMed]
6.	W. Swindell, ed., Benchmark Papers in Optics: Polarized Light , (Hutchinson and Ross, Stroudsburg, 1975).
7.	R.A. Chipman, “Polarimetry,” in Handbook of Optics , Vol. II, (McGraw Hill, New York, 1995).
8.	J.L. Pezzaniti, S.C. McClain, R.A. Chipman, and S.Y. Lu, “Depolarization in liquid crystal TV’s,” Opt. Lett. 18, 2071–2073 (1993). [CrossRef] [PubMed]
9.	S.-M.F. Nee and T.-W. Nee, “Principal Mueller matrix of reflection and scattering for a one-dimensional rough surface,” Opt. Eng. 41, 994–1001 (2002). [CrossRef]
10.	M.W. Williams, “Depolarization and cross polarization in ellipsometry of rough surfaces.” Appl. Opt. 25, 3616–3621 (1986). [CrossRef] [PubMed]
11.	B. Laude-Boulesteix, A. DeMartino, B. Drévilon, and L. Schwartz, “Mueller polarimetric imaging system with liquid crystals,” Appl. Opt. 43, 2824–2832 (2004). [CrossRef] [PubMed]
12.	N. Ghosh, H.S. Patel, and P.K. Gupta, “Depolarization of light in tissue phantoms—effect of a distribution in the size of scatterers,” Opt. Express 11, 2198–2205 (2003). [CrossRef] [PubMed]
13.	K. Sassen, “Cirrus cloud iridescence: a rare case study,” Appl. Opt. 42, 486–491 (2003). [CrossRef] [PubMed]
14.	G.E. Jellison, C.O. Griffiths, D.E. Holcomb, and C.M. Rouleau, “Transmission two-modulator generalized ellipsometry measurements,” Appl. Opt. 41, 6555–6566 (2002). [CrossRef] [PubMed]
15.	W.-N. Chen, C.-W. Chiang, and J.-B. Nee, “Lidar ratio and depolarization ratio for cirrus clouds,” Appl. Opt. 41, 6470–6476 (2002). [CrossRef] [PubMed]
16.	V.A. Ruiz-Cortés and J.C. Dainty, “Experimental light-scattering measurements from large-scale composite randomly rough surfaces,” J. Opt. Soc. Am. A 19, 2043–2052 (2002). [CrossRef]
17.	I. Shoji, Y. Sato, S. Kurimura, V. Lupei, T. Taira, A. Ikesue, and K. Yoshida, “Thermal-birefringence-induced depolarization in NdYAG ceramics,” Opt. Lett. 27, 234–236 (2002). [CrossRef]
18.	S.-M. F. Nee, “Depolarization and principal Mueller matrix measured by null ellipsometry,” Appl. Opt. 40, 4933–4939 (2001). [CrossRef]
19.	M. Moscoso, J.B. Keller, and G. Papanicolaou, “Depolarization and blurring of optical images by biological tissue,” J. Opt. Soc. Am. A 18, 948–960 (2001). [CrossRef]
20.	S.-M.F. Nee, “Depolarization and retardation of a birefringent slab,” J. Opt. Soc. Am. A 17, 2067–2073 (2000). [CrossRef]
21.	J.J. Gil and E. Bernabeau, “Depolarization and polarization indices of an optical system,” Opt. Acta 33, 185–189 (1986). [CrossRef]
22.	B. DeBoo, “Investigation of polarization scatter properties using active imaging polarimetry,” PhD dissertation, Optical Sciences Center, University of Arizona (2004).
23.	D. Goldstein, Polarized Light, Second Ed., Revised and Expanded , (Marcel Dekker, New York, 2003). [CrossRef]
24.	R.C. Yates, A Handbook on Curves and Their Properties , (J. W. Edwards, Ann Arbor, 1952).
25.	G.B. Arfken and H.J. Weber, Mathematical Methods for Physicists , (Academic Press, San Diego, 2001).
26.	S.-Y. Lu, “An interpretation of polarization matrices,” PhD dissertation, Dept. of Physics, University of Alabama at Huntsville (1995).
27.	S.-Y. Lu and R.A. Chipman, “Interpretation of Mueller matrices based on polar decomposition,” J. Opt. Soc. Am. A 13, 1106–1113 (1996). [CrossRef]
28.	C. Brosseau, Fundamentals of Polarized Light: A Statistical Optics Approach , (John Wiley and Sons, Inc., New York, 1998).
29.	C.R. Givens and A.B. Kostinski, “A simple necessary and sufficient condition on physically realizable Mueller matrices,” J. Mod. Opt. 41, 471–481 (1994).
30.	J. Cariou, B. Le Jeune, J. Lotrian, and Y. Guern, “Polarization effects of seawater and underwater targets,” Appl. Opt. 29, 1689- (1990). [CrossRef] [PubMed]

OCIS Codes
(120.5410) Instrumentation, measurement, and metrology : Polarimetry
(230.5440) Optical devices : Polarization-selective devices
(260.5430) Physical optics : Polarization

ToC Category:
Research Papers

History
Original Manuscript: August 20, 2004
Revised Manuscript: September 24, 2004
Published: October 4, 2004

Citation
Brian DeBoo, J. Sasian, and R. Chipman, "Degree of polarization surfaces and maps for analysis of depolarization," Opt. Express 12, 4941-4958 (2004)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-12-20-4941

Sort: Journal | Reset

References

R.C. Jones, �??A new calculus for the treatment of optical systems: I. Description and discussion of the calculus,�?? J. Opt. Soc. Am. 31, 488-493 (1941). [CrossRef]
R.C. Jones, �??New calculus for the treatment of optical systems: VIII. Electromagnetic theory,�?? J. Opt. Soc. Am. 46, 126-131 (1956). [CrossRef]
M. Izdebski, W. Kucharczyk, and R.E. Raab, �??Application of the Jones calculus for a modulated doublerefracted light beam propagating in a homogeneous and nondepolarizing electro-optic uniaxial crystal,�?? J. Opt. Soc. Am. A 21, 132-139 (2004). [CrossRef]
K.D. Abhyankar and A.L. Fymat, �??Relations between the elements of the phase matrix for scattering,�?? J. Math Phys. 10, 1935-1938 (1969). [CrossRef]
E.S. Fry and G.W. Kattawar, �??Relationships between elements of the Stokes matrix,�?? Appl. Opt. 20, 2811-2814 (1981). [CrossRef] [PubMed]
W. Swindell, ed., Benchmark Papers in Optics: Polarized Light, (Hutchinson and Ross, Stroudsburg, 1975).
R.A. Chipman, �??Polarimetry,�?? in Handbook of Optics, Vol. II, (McGraw Hill, New York, 1995).
J.L. Pezzaniti, S.C. McClain, R.A. Chipman, and S.Y. Lu, �??Depolarization in liquid crystal TV�??s,�?? Opt. Lett. 18, 2071-2073 (1993). [CrossRef] [PubMed]
S.-M.F. Nee and T.-W. Nee, �??Principal Mueller matrix of reflection and scattering for a one-dimensional rough surface,�?? Opt. Eng. 41, 994-1001 (2002). [CrossRef]
M.W. Williams, �??Depolarization and cross polarization in ellipsometry of rough surfaces.�?? Appl. Opt. 25, 3616-3621 (1986). [CrossRef] [PubMed]
B. Laude-Boulesteix, A. DeMartino, B. Drévilon, and L. Schwartz, �??Mueller polarimetric imaging system with liquid crystals,�?? Appl. Opt. 43, 2824-2832 (2004). [CrossRef] [PubMed]
N. Ghosh, H.S. Patel, and P.K. Gupta, �??Depolarization of light in tissue phantoms�??effect of a distribution in the size of scatterers,�?? Opt. Express 11, 2198-2205 (2003). [CrossRef] [PubMed]
K. Sassen, �??Cirrus cloud iridescence: a rare case study,�?? Appl. Opt. 42, 486-491 (2003). [CrossRef] [PubMed]
G.E. Jellison, C.O. Griffiths, D.E. Holcomb, and C.M. Rouleau, �??Transmission two-modulator generalized ellipsometry measurements,�?? Appl. Opt. 41, 6555-6566 (2002). [CrossRef] [PubMed]
W.-N. Chen, C.-W. Chiang, J.-B. Nee, �??Lidar ratio and depolarization ratio for cirrus clouds,�?? Appl. Opt. 41, 6470-6476 (2002). [CrossRef] [PubMed]
V.A. Ruiz-Cortés, and J.C. Dainty, �??Experimental light-scattering measurements from large-scale composite randomly rough surfaces,�?? J. Opt. Soc. Am. A 19, 2043-2052 (2002). [CrossRef]
I. Shoji, Y. Sato, S. Kurimura, V. Lupei, T. Taira, A. Ikesue, and K. Yoshida, �??Thermal-birefringenceinduced depolarization in NdYAG ceramics,�?? Opt. Lett. 27, 234-236 (2002). [CrossRef]
S.-M. F. Nee, �??Depolarization and principal Mueller matrix measured by null ellipsometry,�?? Appl. Opt. 40, 4933-4939 (2001). [CrossRef]
M. Moscoso, J.B. Keller, and G. Papanicolaou, �??Depolarization and blurring of optical images by biological tissue,�?? J. Opt. Soc. Am. A 18, 948-960 (2001). [CrossRef]
S.-M.F. Nee, �??Depolarization and retardation of a birefringent slab,�?? J. Opt. Soc. Am. A 17, 2067-2073 (2000). [CrossRef]
J.J. Gil and E. Bernabeau, �??Depolarization and polarization indices of an optical system,�?? Opt. Acta 33, 185-189 (1986). [CrossRef]
B. DeBoo, �??Investigation of polarization scatter properties using active imaging polarimetry,�?? PhD dissertation, Optical Sciences Center, University of Arizona (2004).
D. Goldstein, Polarized Light, Second Ed., Revised and Expanded, (Marcel Dekker, New York, 2003). [CrossRef]
R.C. Yates, A Handbook on Curves and Their Properties, (J. W. Edwards, Ann Arbor, 1952).
G.B. Arfken, and H.J. Weber, Mathematical Methods for Physicists, (Academic Press, San Diego, 2001).
S.-Y. Lu, �??An interpretation of polarization matrices,�?? PhD dissertation, Dept. of Physics, University of Alabama at Huntsville (1995).
S.-Y. Lu and R.A. Chipman, �??Interpretation of Mueller matrices based on polar decomposition,�?? J. Opt. Soc. Am. A 13, 1106-1113 (1996). [CrossRef]
C. Brosseau, Fundamentals of Polarized Light: A Statistical Optics Approach, (John Wiley and Sons, Inc., New York, 1998).
C.R. Givens and A.B. Kostinski, �??A simple necessary and sufficient condition on physically realizable Mueller matrices,�?? J. Mod. Opt. 41, 471-481 (1994).
J. Cariou, B. Le Jeune, J. Lotrian, and Y. Guern, �??Polarization effects of seawater and underwater targets,�?? Appl. Opt. 29, 1689- (1990) [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.

Click here to see a list of articles that cite this paper

Figures


Fig. 1.	Fig. 2.	Fig. 3.


Fig. 4.	Fig. 5.	Fig. 6.


Fig. 7.	Fig. 8.	Fig. 9.


Fig. 10.	Fig. 11.	Fig. 12.


Fig. 13.	Fig. 14.	Fig. 15.


Fig. 16.