## 1.

## INTRODUCTION

In the second-order approach, the polarization effects of a material medium on the electromagnetic wave that interacts with it are characterized by means of the corresponding Mueller matrix. Beyond its mere role as the 4×4 matrix that transforms the Stokes vectors of the incoming polarization states into the Stokes vectors of the outgoing states, the Mueller matrix contains, in a intricate manner, rich information on the physical nature and properties of the medium.

This work is devoted to identify and interpret the main physical quantities involved in a Mueller matrix. After an analysis of each significant quantity, the notions of components of purity, indices of polarimetric purity as well as the sets of quantities that remain invariant under certain transformations are revisited and interpreted on the basis of a common framework.

For the sake of self-consistency of this work, it is worth considering briefly the concept of Mueller matrix. The transformation of the state of polarization of an electromagnetic beam caused by its linear interaction with a material medium can be formulated mathematically as **Ms** = **s**′, where **s** and **s**′ are the respective Stokes vectors of the input and output states of polarization and **M** is the Mueller matrix associated with the medium (for the given particular conditions of interaction). In general, the medium can exhibit dispersive effects or certain degree of heterogeneity over the area illuminated by the incident electromagnetic beam, producing depolarization. Therefore, the emerging electromagnetic wave is composed of a number of incoherent contributions (mutually coherent or incoherent), in such a manner that the polarimetric effect of the medium not always can be represented by means of a Mueller-Jones matrix (deterministic case).

In other words, while some polarimetric interactions have a deterministic nature and can be formulated either by the Jones calculus or by the Stokes-Mueller calculus (in which case the Mueller matrix **M** is *pure*, also called *nondepolarizing* or *Mueller-Jones* matrix, that is, **M** is derived from a unique Jones matrix), in the most general case, **M** is a convex sum, or ensemble average, of pure Mueller matrices (i.e. a linear combination of pure Mueller matrices with positive coefficients that sum to one). This is a consequence of the fact that, essentially, any macroscopic polarimetric effect is due to a sort of scattering by a myriad of elementary electronic excitations, each one with a well-defined Jones matrix, and the whole action of the medium is a statistical composition of such *pure* interactions. This physical requisite is called the *ensemble criterion* [1] (also *covariance criterion* or *Cloude’s criterion* [2]).

Note that the ensemble criterion imposes additional restrictions beyond the fact that a Mueller matrix is a *Stokes matrix* transforming Stokes vectors into Stokes vectors [3]. A simple example of a Stokes matrix that is not a Mueller matrix is the Minkowski metric **G** ≡ diag (1, −1, −1, −1), which plays a key role in the Stokes-Mueller algebra [5-12] but has no associated Jones matrix.

Furthermore, leaving aside certain artificial arrangements [13], the linear polarimetric effects are *passive* [3,5,14,15], in the sense that the *intensity* (power density flux) of the output electromagnetic wave never exceeds the intensity of the input wave.

Thus, by considering jointly the ensemble and passivity criteria, we can state that a 4×4 real matrix **M** is a Mueller matrix if and only if it can be expressed as

where **M*** _{Ji}* are pure and passive Mueller matrices (i.e.

**M**

*are associated with respective passive Jones matrices). For a more detailed analysis of this subject we refer the reader to Ref. [1].*

_{Ji}Before identifying the relevant physical parameters involved in a general Mueller matrix **M**, let us first recall that the structure of the information contained in **M** is closely related to its block expression [16],

where **M**̂ is the normalized version of **M**.

For some analyses it is also useful to parameterize the submatrix **m** as follows [17]

so that **M** is expressed in terms of *m*_{00} and the five *constitutive vectors*

It is also worth to consider briefly the reciprocity properties of Mueller matrices, that is, the relation between **M**and the Mueller matrix **M*** ^{r}* that corresponds to the same medium as

**M**, but interchanging input and output directions of the probe wave interacting with it. It has been proved that (leaving aside systems involving magnetooptic effects, whose reciprocity property

**M**

*=*

^{r}**M**differs from the usual rule) the

*reciprocal Mueller matrix*

**M**

*is given by [18,19]*

^{r}where the superscript *T* indicates transposed matrix. Therefore, the knowledge of a depolarizing Mueller matrix **M** fully determines the reciprocal Mueller matrix **M*** ^{r}*. In fact, it has been proved that the constraining inequalities that characterize the Mueller matrices are invariant with respect to the change of

**M**

*by*

^{r}**M**

*[5,1], and consequently all the physical conditions for a matrix*

^{T}**M**to be a Mueller matrix must also be satisfied by

**M**

*and*

^{r}**M**

*. In other words,*

^{T}**M**is a Mueller matrix if and only if

**M**

*(or*

^{T}**M**

*) is a Mueller matrix.*

^{r}The state of polarization of an electromagnetic beam with a well defined direction of propagation can be represented by means of the corresponding Stokes vector

where *I* is the *intensity*, *P* is the *degree of polarization* (a measure of how close is the state represented by **s** to a totally polarized one), while ϕ and χ are respectively the *azimuth* and *ellipticity angle* of the average polarization ellipse (called the *characteristic polarization ellipse* [20]).

The whole set of polarization states **s**̂ ≡ **s**/ *I* (i.e., states with unit intensity and arbitrary values of *P*, ϕ and χ) constitute the Poincaré sphere. States with *P* = 1 (said to be *pure*, or *totally polarized*) constitute the surface of the sphere, while states satisfying 0 ≤ *P* < 1 lie inside the sphere (*P* = 0 corresponding to the origin).

To inspect some polarimetric properties of **M**, it will be useful to keep in mind the following transformations, through **M** and through **M*** ^{T}*, of an input state

**s**̂ into the respective output states

**s**

*and*

^{f}**s**

^{r}## 2.

## MEAN INTENSITY COEFFICIENT

Let us consider the pair of normalized pure Stokes vectors

which are said to be *mutually orthogonal* (i.e., they satisfy ) and are represented by antipodal points on the surface of the Poincaré sphere.

The effect of a medium represented by a Mueller matrix **M** on the following equiprobable incoherent mixture of the states **s**̂* _{p+}* and

**s**̂

_{p−}is given by the outgoing Stokes vector

whose intensity is *m*_{00}. Since each pure state **s**̂ * _{p+}* has its respective orthogonal state

**s**̂

*, it turns out that an equiprobable incoherent mixture of all pure states (thus covering the entire surface of the Poincaré sphere) is given by a unpolarized state*

_{p−}and therefore, from (11), the intensity of the corresponding output state is *m*_{00}. Furthermore, if we consider a unpolarized input state with intensity *I*, the intensity of the transformed state is *m*_{00} *I*. This fact provides the physical meaning of *m*_{00} and justifies calling it the *mean intensity coefficient* of **M** (also termed *transmittance* -or reflectance- or *gain*).

## 3.

## POLARIZANCE AND DIATTENUATION

A given Mueller matrix transforms an input unpolarized state as shown in Eq. (11), so that the ability of **M** to polarize input unpolarized states is characterized by the *polarizance vector* **P**, whose absolute value *P* is the degree of polarization of the transformed state and is called the *polarizance* of **M** [21,22]

The upper limit *P* = 1 of *P* (which is restricted to the interval 0 ≤ *P* ≤ 1) is achieved for media that fully polarize the input unpolarized states. Thus, a medium satisfying *P* = 1 is called a *polarizer* and his associated Mueller matrix **M*** _{P̂D}* has the form [23,24]

Moreover, the polarizance of a Mueller matrix **M** can be expressed as [17]

where the *degree of linear polarizance P _{L}* is a measure of the ability of

**M**to transform a unpolarized state into a linearly polarized state (0 ≤

*P*≤ 1), while the

_{L}*degree of circular polarizance P*is a measure of the ability of

_{C}**M**to C;transform a unpolarized state into a circularly polarized one (|

*P*| ≤ 1, with

_{C}*P*= −1 for a left-handed circular polarizer and

_{C}*P*= 1 for a right-handed circular polarizer).

_{C}To analyze the physical information involved in the *diattenuation vector* **D**, let us observe that, from Eq. (7), the and minimum intensities for the output state **s*** ^{f}* =

**Mŝ**are given by

which correspond respectively to the following mutually-orthogonal input states

so that the absolute value *D* of **D** can be expressed as [5,25,1]

showing that the *diattenuation D* is a measure of the relative difference between the maximum and minimum intensities of the output states with respect to all possible input polarization states. *D* is limited by 0 ≤ *D* ≤ 1, where *D* = 1 corresponds to media for which there exists an input state (**s**_{−D̂}) whose corresponding output state has zero intensity. When *D* = 1 the medium is called an *analyzer* and has the general form [23,24]

As for polarizance, diattenuation can be expressed as a quadratic average of linear diattenuation *D _{L}* and circular diattenuation

*D*[17]

_{C}Let us now consider the reverse **M*** ^{r}* and transpose

**M**

*Mueller matrices associated with*

^{T}**M**, and observe that, from the mere definition of

*P*,

*D*,

*D*,

_{L}*D*,

_{C}*P*,

_{L}*P*, these quantities satisfy

_{C}In other words, the polarizance of **M** is also the diattenuation (or *reciprocal polarizance*) of **M*** ^{r}* and

**M**

*, while the diattenuation of*

^{T}**M**is also the polarizance (or

*reciprocal diattenuation*) of

**M**

*and*

^{r}**M**

*[4].*

^{T}## 4.

## POLARIZANCE ANGLE

Vectors **P** and **D**, whose respective physical meanings have been analyzed in Section 3, can be represented in the Poincaré sphere, and the *polarizance angle η* formed by them and defined as

is an interesting quantity involved in **M** that remains invariant under certain transformations of **M** (this subject will be considered in further sections on retarder and rotation transformations of **M**). In particular, the role played by *η* in the classification of nondepolarizing Mueller matrices is studied in Ref. [27].

## 5.

## DEGREE OF POLARIZANCE

As we have seen, despite the fact that both polarizance *P* and diattenuation *D* have respective specific physical meanings, they rely on a common dual physical nature, so that, for some purposes, it is useful to define the *degree of polarizance P _{P}* [28]

which provides a measure of the combined contribution of *P* and *D* to polarimetric purity. *P _{P}* is restricted to the range 0 ≤

*P*≤ 1,

_{P}*P*= 1 corresponding to a pure polarizer and

_{P}*P*= 0 corresponding to a

_{P}*nonpolarizing Mueller matrix*

**M**(with zero diattenuation and zero polarizance).

Concerning the notion of the degree of polarizance, it is worth to recall that, contrary to that it would seem at first sight, *P _{P}* is related to the ability of the medium represented by

**M**to reduce the degree of polarization of certain input states of polarization [29]. To analyze this feature, let us consider the example of a

*normal diattenuator*[30] (which is a kind of pure Mueller matrix [1]), with diattenuation 0 <

*D*< 1, whose Mueller matrix has the form [23]

The action of **M***D* on the partially polarized state (1,**D*** ^{T}*)

*, whose degree of polarization*

^{T}*D*equals the diattenuation of

**M**

*D*, produces an output partially polarized state represented by the following Stokes vector

whose degree of polarization is *D* (1 + *D*^{2}) < *D**D* and, consequently, the fact that *D* > 0 implies that **M***D* reduces the degree of polarization of some input partially polarized states.

## 6.

## DEGREE OF SPHERICAL PURITY

The *degree of spherical purity, P _{S}*, is defined as [29,31]

so that it is limited by 0 ≤ *P _{S}* ≤ 1 and constitutes a measure of how close is

**M**to the Mueller matrix of a retarder (in general elliptic, and regardless of the effective value of its retardance) [26]. That is,

*P*= 1 corresponds to a Mueller matrix whose normalized version is the Mueller matrix

_{S}**M**

*of a retarder (i.e. and det*

_{R}**M**

*= +1). The lower limit and*

_{R}*P*= 0 corresponds to a depolarizing medium satisfying

_{S}**m = 0**, which is fully characterized by

*m*

_{00},

**P**and

**D**.

## 7.

## DEGREE OF POLARIMETRIC PURITY

Polarimetric purity refers to the closeness of a Mueller matrix to that of a deterministic nondepolarizing medium (i.e., a *pure* medium, hence characterized by the fact that it preserves the degree of polarization of any input totally polarized state in both forward and reverse directions). An overall measure of the degree of polarimetric purity of **M** is given by the *depolarization index P*_{Δ} defined as [22]

This quantity is limited by 0 ≤ *P*_{Δ} ≤ 1, the upper limit *P*_{Δ} = 1 corresponding to a pure medium (regardless of its particular retardance and diattenuation properties) and the lower limit *P*_{Δ} = 0 corresponding to an ideal depolarizer with associated Mueller matrix **M**_{Δ0} = *m*_{00} diag (1, 0, 0, 0).

In accordance with the physical meaning of *P*_{Δ}, an overall measure of the depolarizing power (or polarimetric randomness) of the medium is given by the *depolarizance* [26]

Pure Mueller matrices are characterized by *P*_{Δ} = 1 (*D*_{Δ} = 0), while Mueller matrices satisfying *P*_{Δ} < 1 (*D*_{Δ} > 0) are called *nonpure* or *depolarizing*.

Despite the fact that *P*_{Δ} constitutes a well-defined overall measure of the polarimetric purity, it has been shown that it does not provide enough information for a complete parameterization of the polarimetric purity of **M** [4,31]. Such parameterization can be obtained from two complementary sets of quantities called the *components of purity* and the *indices of polarimetric purity* which will be described in respective sections.

## 8.

## COMPONENTS OF PURITY

The block expression of a Mueller matrix in Eq. (2) is particularly significant because it reflects relevant properties of **M** concerning the structure of the physical information contained in it. Accordingly, and from the inspection of the right term of Eq. (26) the quantities *P*, *D* and *P _{S}*, taken as a set, are called the

*components of purity*(hereafter CP) of

**M**. The analysis of the particular values of the CP allows for a classification of Mueller matrices (see Refs. [20,26]).

The parameters *P _{P}* and

*P*, taken as a set, are called the

_{S}*sources of purity*of

**M**, which represent complementary contributions to the overall polarimetric purity given by

*P*

_{Δ}[28]

and their representation in the purity figure [26] provides a meaningful view of the different types of Mueller matrices (see Fig. 1 in Section 9).

## 9.

## INDICES OF POLARIMETRIC PURITY

While the components of purity provide information of how the polarimetric purity is shared among polarizance, diattenuation and spherical purity, a complementary view of how the polarimetric purity (or, conversely, its randomness) is quantitatively structured, is achieved from the eigenvalues of the covariance matrix **H** associated with **M** and defined as

where ⊗ stands for the Kronecker product and **σ*** _{k}* are the Pauli matrices (ordered as usual in Optics)

Since **H** is a positive semidefinite Hermitian matrix, its eigenvalues are nonnegative, and, taken in decreasing order, are denoted by *λ _{i}* (

*λ*

_{0}≥

*λ*

_{1}≥

*λ*

_{2}≥

*λ*

_{3}). Thus,

**H**can be expressed as

where **U** is the unitary matrix whose columns are the eigenvectors **u**_{i} and the so-called *spectral decomposition*, or *Cloude’s decomposition* of **H** is defined as [2,3]

where *λ*̂ ≡ *λ* (tr **H**) and **H**_{Ji} are the covariance matrices associated with respective pure Mueller matrices **M***Ji* (the subscript *J* is used to indicate that these matrices are associated with pure Mueller matrices, that is, each **H*** _{Ji}* has only one nonzero eigenvalue

*λ*

*). The*

_{i}*spectral decomposition*of

**M**is formulated as follows in terms of the pure Mueller matrices

**M**

*≡*

_{Ji}*m*

_{00}

**M̂**

*with the same mean intensity coefficients [3] equal to*

_{Ji}*m*

_{00}= tr

**H**

The above decomposition means that, from a polarimetric point of view, the medium represented by **M** is equivalent to a *parallel composition* [32,33] of the pure media represented by **M*** _{Ji}*, that is, the medium behaves like if the input electromagnetic beam was shared among a set of pencils interacting separately with the pure components represented by

**M**

*, in such a manner that the emerging pencils recombine incoherently into a whole output beam [3,4]. Consequently, the normalized eigenvalues*

_{Ji}*λ̂*of

_{i}**H**represent the relative weights of the pure components

**M**

*Ji*in the spectral decomposition.

An objective view of how the polarimetric purity is structured in **M** is obtained by rearranging the spectral decomposition in the form of the following *trivial*, or *characteristic*, decomposition [4]

where the quantities *P*_{1}, *P*_{2} and *P*_{3} are the so-called *indices of polarimetric purity* (IPP) of **M**, defined as follows from the normalized eigenvalues of **H** [34]

while **M̂**_{J0}, **M̂**_{J1}, **M̂**_{J2} and **M̂**_{Δ0} are normalized Mueller matrices associated with the covariance matrices **Ĥ**_{J 0}, **Ĥ**_{J1}, **Ĥ**_{J2}, and **Ĥ**_{Δ0} defined as

where **I** is the 4×4 identity matrix. The above expressions allows for the following interpretation of the characteristic components of **M**:

•

**M**_{J1}≡*m*_{00}**M**̂_{J1}is called the*characteristic pure component*(note that it coincides with the first component of the spectral decomposition) and its relative weight in the characteristic decomposition (34) is precisely the first index of polarimetric purity, or*degree of polarization P*_{1}of**M**.•

**M**_{2}≡*m*_{00}**M̂**_{2}represents a*2D depolarizer*, and its relative weight is given by the difference*P*_{2}−*P*_{1}.•

**M**_{3}≡*m*_{00}**M**̂_{3}represents a*3D depolarizer*, and its relative weight is given by the difference*P*_{3}−*P*_{2}.•

**M**_{Δ0}≡*m*_{00}diag (1, 0, 0, 0) represents an*ideal depolarizer*(or*4D depolarizer*), whose relative weight is given by the difference 1 −*P*_{3}.

The overall degree of polarimetric purity *P*_{Δ} can be calculated through the following weighted quadratic average of the IPP [34]

It is remarkable that, from the role played by the IPP in the characteristic decomposition and from the scaled structure of their values [34]

the IPP constitute a privileged set of parameters providing quantitative information of how the polarimetric purity (or, conversely, the polarimetric randomness) is organized in the medium represented by **M**. Note that the IPP are insensitive to the specific features of **M** relative to its CP.

Although IPP and CP provide independent and complementary information, note that the set (*P*_{1,} *P*_{2,} *P*_{3,} *P*, *D*) of *purity parameters* is sufficient to calculate *P*_{Δ} and *P _{S}*. The purity figure (Fig. 1) provides a powerful geometric tool to inspect the intricate relations between the IPP and the CP. A detailed analysis of the physical interpretation of both sets of parameters, can be found in Ref. [26]

## 10.

## RETARDER TRANSFORMATIONS

The term *retarder* is used for pure media with zero diattenuation-polarizance, *P* = *D* = 0. The Mueller matrix of an ideal retarder **M*** _{R}* has the general form

so that a given Mueller matrix **M** corresponds to a retarder if and only if **M*** ^{T}* =

**M**

^{−1}, that is, if and only if

**M**is orthogonal.

The action of **M*** _{R}* can be represented in the Poincaré sphere as a rotation by the angle Δ about the axis defined by the unit

*Pauli vector*

**u**

*of*

_{R}**M**

*, so that any input state is transformed into an output state obtained through such rotation. The two eigenstates of*

_{R}**M**

*are given by the mutually orthogonal and totally polarized Stokes vectors*

_{R}Thus, a given retarder can also be represented by means of its corresponding *straight retardance vector* **R** defined as

where **u*** _{R}* is the Pauli vector determining the azimuth and ellipticity of the fast eigenstate of

**M**

*and Δ is the*

_{R}*retardance*. Note that

**R**has been defined so as to satisfy 0 ≤

*R*≤ 1 which allows for its representation in the Poincaré sphere.

Unlike Mueller matrices of retarders, Mueller matrices (pure or nonpure) with nonzero diattenuation or polarizance, are necessarily associated with media that exhibit forward or reverse diattenuation, which is a property linked to selective absorption (dichroism) or selective reflection, hence producing power loss, in the sense that such effects cannot be compensated by means of serial combinations with additional passive media. Therefore, unlike the action of diattenuators, the action of a retarder produces lossless effects on the interacting electromagnetic wave. In fact, given **M*** _{R}*, always exists a passive pure Mueller matrix

**M**

_{R}

^{−1}(which also corresponds to a retarder) whose serial combination with the retarder represented by

**M**

*,*

_{R}**M**

_{R}

^{−1}

**M**

_{R}=

**I**produces a neutral effect.

To get a deeper knowledge of the information contained in **M**, it is therefore interesting to analyze its possible lossless transformations through its serial combination with other Mueller matrices. The most general form of a lossless transformation of **M** is that of a *dual retarder transformation* [17]

whose physical interpretation is that the medium represented by **M** is sandwiched by the retarders represented by **M**_{R1} (*entrance retarder*) and **M**_{R2} (*exit retarder*), so that this serial combination is represented by the resulting Mueller matrix **M**′. Dual retarder transformations preserve important quantities like *m*_{00}, *D*, *P*, *P _{S}*,

*P*

_{1},

*P*

_{2},

*P*

_{3}and

*P*

_{Δ}, thus including those that fully characterize the depolarizing properties of

**M**; this is the reason why Mueller matriceslinked to

**M**through dual retarder transformations are said to be

*invariant-equivalent*to

**M**[31].

Among the infinite Mueller matrices that are invariant-equivalent to **M**, it is worth to pay attention to the so-called *arrow form* **M*** _{A}* of

**M**[35], which is built from the singular value decomposition of the submatrix

**m**of

**M**

where **m*** _{RI}* and

_{m}*are proper orthogonal matrices, and*

_{RO}**m**

*is the diagonal matrix whose elements are defined as follows from the singular values (*

_{A}*a*

_{1},

*a*

_{2},

*a*

_{3}) of

**m**,

The arrow form **M*** _{A}* of

**M**is defined as [35]

so that the dual retarder transformation **M** = **M**_{RO}**M**_{A}**M*** _{RI}* is called the

*arrow decomposition*of

**M**.

The main feature of the arrow form **M*** _{A}* is that the six off-diagonal elements of

**m**

*are zero, so that the ten nonzero elements of*

_{A}**M**

*provide all the information of*

_{A}**M**that is invariant under dual retarder transformations. In fact,

**M**

*is built from*

_{A}*m*

_{00}together with (1) the diattenuation vector

**D**

*≡*

_{A}**m**

_{RI}**D**, (2) the polarizance vector

**P**

*≡*

_{A}**m**

^{T}

_{RO}**P**, and (3) the

*spherical vector*

whose absolute value equals the degree of spherical purity

Other interesting class of retarder transformations is the *single retarder transformation*, defined as an orthogonal similarity transformation of **M** [17]

whose physical interpretation is that the medium represented by **M** is sandwiched by two identical retarders whose fast eigenstates are mutually orthogonal.

## 11.

## ROTATION TRANSFORMATIONS

This section is devoted to dual retarder transformations that can be physically performed by means of rotations of the respective laboratory reference frames used for the representation of the Stokes vectors of the input and output electromagnetic beam [17].

The matrix **M*** _{G}* that transforms a Mueller matrix

**M**into the Mueller matrix of the same medium but with respect to a rotated laboratory reference frame, is commonly called a rotator and has the form

where *θ* is the angle rotated. Note that **M*** _{G}* (

*θ*) is a Mueller matrix. In fact, it coincides with the Mueller matrix of circular retarder with retardance 2

*θ*[1].

A *dual-rotation transformation* [17] is a kind of dual retarder transformation where **M**_{R2} and **M**_{R1} have the form of respective rotation matrices **M**_{R2} = **M*** _{G}* (

*θ*

_{2}),

**M**

_{R1}=

**M**

*(*

_{G}*θ*

_{1})

By using the five-vector expression of **M** in terms of the elements of its constitutive vectors **D**, **P**, **k**, **r** and **q**, [17]

the dual rotation transformations can be expressed as [17]

An especially interesting subclass of dual rotation transformations is that constituted by *single rotation transformations* [17], which correspond to the case that the laboratory reference frames for the representation of input and output polarization states are rotated jointly (and, hence, applicable for the case of direct transmission experiments where input and output reference frames coincide)

When expressed in terms of the elements of the five constitutive vectors of **M**, the single rotation transformation adopts the form [17]

## 12.

## PHYSICAL PARAMETERIZATION OF A MUELLER MATRIX

The retarder and rotation transformations provide meaningful ways for describing the physical information contained in **M** in terms of certain sets of quantities.

From the point of view of dual retarder transformations and, in particular, from the arrow decomposition of **M**, we get that a complete set of sixteen physical parameters is provided by

that is, by

• The mean intensity coefficient

*m*_{00}• The

*entrance straight retardance vector***R**, which contains all the information of the entrance equivalent retarder through its Pauli vector_{I}**u**and its retardance Δ_{I}. Thus,_{I}**R**fully determines the entrance retarding properties of_{I}**M**.• The

*exit straight retardance vector***R**, which contains all the information of the exit equivalent retarder through its Pauli vector_{O}**u**and its retardance_{O}*Δ*. Thus,_{O}**R**fully determines the exit retarding properties of_{O}**M**.• The diattenuation vector , which fully determines the diattenuation properties of

**M**.• The polarizance vector

**P**=**m**_{RO}**P**, which fully determines the polarizance properties of_{A}**M**.• The three indices of polarimetric purity (

*P*_{1},*P*_{2},*P*_{3}), which fully determines the quantitative depolarizing properties of**M**.

As we have seen in previous sections, the *quantitative* characterization of the polarimetric purity of **M** is given by the IPP: *P*_{1}, *P*_{2}, *P*_{3} (which, obviously, also determine the quantitative structure of the polarimetric randomness and depolarizing properties of **M**), while the *qualitative* characterization of the polarimetric purity (or of the corresponding sources of depolarization) of **M** is given by the CP: *P*, *D*, *P _{S}*. Moreover, it has been demonstrated [36] that any depolarizing Mueller matrix

**M**has an associated reference pure Mueller matrix

**M**

*(*

_{J}**M**) to which

**M**is reduced when the IPP of

**M**are replaced by

*P*

_{1}=

*P*

_{2}=

*P*

_{3}= 1. That is, the IPP act as regulators of the quantitative structure of the polarimetric purity of

**M**, while the CP provide detailed information on the sources and nature of the polarimetric purity associated with

**M**.

In virtue of the dual retarder transformation of a pure Mueller matrix, which can always be written as [1,37]

where the central horizontal diattenuator has the form

the *entrance equivalent retarder* **M**_{R1} is determined by a sequence of a linear retarder (with retardance Δ* _{I}* and whose fast linear eigenstate has azimuth

*ϕ*) and a horizontal linear retarder with retardanceΔ/2, while the

_{I}*exit equivalent retarder*

**M**

_{R2}is determined by a sequence of a horizontal linear retarder with retardance Δ/2 and a linear retarder (with retardance Δ

*O*and whose fast linear eigenstate has azimuth

*ϕ*).

_{O}Therefore, when the IPP of **M** are replaced by their values *P*_{1} = *P*_{2} = *P*_{3} = 1 corresponding to total purity, the six parameters associated with the set composed of the entrance and exit retarders, **M*** _{RI}* and

**M**

*, of the arrow decomposition of*

_{RO}**M**,

*into the five parameters*

_{collapse}*ϕ*, Δ

_{I}*, Δ,*

_{I}*ϕ*, Δ

_{O}*, while the nine parameters of*

_{O}**M**

*into the two parameters of*

_{A collapse}**M**

_{DL0}(

*m*

_{00},

*D*). Consequently, contrary to the apparent necessity of nine parameters for the complete integral quantitative and qualitative characterization of the depolarization properties of

**M**, they are fully determined by its five purity parameters

*P*

_{1},

*P*

_{2},

*P*

_{3},

*P*,

*D*. Obviously, the complete joint characterization of retarding, polarizing, diattenuating and depolarizing properties of

**M**is given by

**R**

*,*

_{I}**R**

*,*

_{O}**P**,

**D**,

*P*

_{1},

*P*

_{2},

*P*

_{3}(fifteen parameters that, together with the mean intensity coefficient

*m*

_{00}, complete the sixteen physical parameters that fully characterize

**M**).

Thus, the above *physical parameterization* of **M** has important consequences for the physical interpretation of Mueller matrices. In fact, the dual retarder transformation that links **M** and its arrow form **M*** _{A}* is reversible in the sense that it does not involve diattenuation effects (the CP of

**M**

*are equal to those of*

_{A}**M**) or changes in the structure of polarimetric purity (the IPP of

**M**

*are equal to those of*

_{A}**M**). Conversely, transformations involving diattenuators or depolarizers necessarily affect to the CP or to the IPP.

Obviously, there are many quantities that can be derived from the above minimum set of sixteen independent parameters in Eq. (54), like, for instance, *P*_{Δ}, *D*_{Δ} and *P _{S}*.

It has been demonstrated [17] that other alternative parameterizations of **M**, based on dual retarder transformations, are possible, as for instance the one obtained by replacing *P*_{1}, *P*_{2} and *P*_{3} in the set (52) by the scalar quantities *P _{S}*, det

**M**(whose physical meaning is studied in Refs. [24,38]), and

**P**

^{T}**mD**(whose physical interpretation would require further analysis). Other possible choice is that constituted by

## 13.

## OTHER PARAMETERIZATIONS OF A MUELLER MATRIX

Concerning the parameterization of **M** in terms of the single retarder transformation, let us observe that a complete set is constituted by [17]

Note that **P**^{T}**D** = *PD* cos*η*, where *η* is the polarizance angle. Nevertheless, the physical interpretation of **P**^{T}**m D**, **P**^{T}**m**^{T}**D** and tr**M** is not straightforward and require additional study.

From the point of view of dual rotation transformations, different sets of sixteen mutually independent parameters can be chosen, as for example, from Eq. (51), [17]

where, in analogy to the linear and circular components of vectors **P** and **D**, the *linear* and *circular* components of vectors **r**, **q** and **k** have been defined as

And finally, from the single rotation transformation, the following parameterization in terms of sixteen independent parameters is obtained [17]

or, alternatively

## 14.

## CONCLUSION

The relevant physical information contained in a Mueller matrix **M** has been identified, decoupled and interpreted on the basis of the arrow decomposition of **M**:

• The mean intensity coefficient of

**M**(mean transmittance or gain) is given by*m*_{00}(one parameter).• The

*birefringence*properties of**M**are fully determined by the entrance and exit straight retardation vectors,**R**and_{I}**R**(3+3 parameters)._{O}• The

*polarizance*and*diattenuation*properties are fully determined by the respective polarizance and diattenuation vectors**P**and**D**(3+3 parameters).• The quantitative structure of the polarimetric purity (or, conversely, the polarimetric randomness) is fully determined by the indices of polarimetric purity

*P*_{1},*P*_{2},*P*_{3}(3 parameters).

Complete quantitative and qualitative information on the depolarizing properties of **M** is provided by the set of five *purity parameters P*_{1}, *P*_{2}, *P*_{3}, *P*, *D*, from which other relevant physically invariant quantities relative to polarimetric purity, like *P _{S}*,

*P*

_{Δ}or

*D*

_{Δ}can be calculated

**Funding**. This research was supported by Ministerio de Economía y Competitividad, grants FIS2011-22496 and FIS2014-58303-P, and by Gobierno de Aragón (E99).

## REFERENCES

Gil, J. J. and Ossikovski, R., [Polarized light and the Mueller matrix approach] CRC Press. Boca Raton (2016).Google Scholar

Cloude, S. R., “Group theory and polarization algebra,” Optik 75, 26-36 (1986).Google Scholar

Gil, J. J., “Characteristic properties of Mueller matrices,” J. Opt. Soc. Am. A 17, 328-334 (2000). https://doi.org/10.1364/JOSAA.17.000328Google Scholar

Gil, J. J., “Polarimetric characterization of light and media,” Eur. Phys. J. Appl. Phys. 40, 1-47 (2007). https://doi.org/10.1051/epjap:2007153Google Scholar

Gil, J. J., [Determination of polarization parameters in matricial representation. Theoretical contribution and development of an automatic measurement device]. PhD dissertation, University of Zaragoza (1983). Available from http://zaguan.unizar.es/record/10680/files/TESIS-2013-057.pdf (1983)Google Scholar

Arnal, P. M., [Modelo matricial para el estudio de fenómenos de polarización de la luz]. PhD dissertation. University of Zaragoza (1990).Google Scholar

van der Mee, C. V. M., “An eigenvalue criterion for matrices transforming Stokes parameters,” J. Math. Phys. 34, 5072-5088 (1993). https://doi.org/10.1063/1.530343Google Scholar

Bolshakov, Y. van der Mee, C. V. M., Ran, A. C. M., Reichstein, B. and Rodman, L., “Polar decompositions in finite dimensional indefinite scalar product spaces: special cases and applications,” in [Operator Theory: Advances and Applications, I]. Gohberg, ed., Vol. 87, pp. 61-94. Birkhäuser Verlag, (1996)Google Scholar

Bolshakov, Y. van der Mee, C. V. M., Ran, A. C. M., Reichstein, B. and Rodman, L., “Errata for: Polar decompositions in finite dimensional indefinite scalar product spaces: special cases and applications,” Integral Equation Oper. Theory 27, 497–501 (1997). https://doi.org/10.1007/BF01192129Google Scholar

Gopala Rao, A.V. Mallesh, K.S. and Sudha, “On the algebraic characterization of a Mueller matrix in polarization optics. I. Identifying a Mueller matrix from its *N* matrix”, J. Mod. Opt. 45, 955-987 (1998).Google Scholar

Gopala Rao, A.V. Mallesh, K.S. and Sudha, “On the algebraic characterization of a Mueller matrix in polarization optics. II. Necessary and sufficient conditions for Jones derived Mueller matrices,” J. Mod. Opt. 45, 989-999 (1998).Google Scholar

Ossikovski, R.. “Analysis of depolarizing Mueller matrices through a symmetric decomposition”, J. Opt. Soc. Am. A 26, 1109-1118 (2009). https://doi.org/10.1364/JOSAA.26.001109Google Scholar

Opatrný, T. and Perina, J.. “Non-image-forming polarization optical devices and Lorentz transformations - an analogy,” Phys. Lett. A 181, 199-202 (1993). https://doi.org/10.1016/0375-9601(93)90639-HGoogle Scholar

Barakat, R. “Conditions for the physical realizability of polarization matrices characterizing passive systems,” J. Mod. Opt. 34, 1535-1544 (1987). https://doi.org/10.1080/09500348714551471Google Scholar

Kostinski, A. B. and Givens, R. C., “On the gain of a passive linear depolarizing system,” J. Mod. Opt. 39, 1947-1952 (1992). https://doi.org/10.1080/09500349214552001Google Scholar

Xing, Z-F, “On the deterministic and non-deterministic Mueller matrix,” J. Mod. Opt. 39, 461-484 (1992). https://doi.org/10.1080/09500349214550471Google Scholar

Gil, J. J., “Invariant quantities of a Mueller matrix under rotation and retarder transformations,” J. Opt. Soc, Am. A 33, 52-58 (2016). https://doi.org/10.1364/JOSAA.33.000052Google Scholar

Sekera, Z., “Scattering Matrices and Reciprocity Relationships for Various Representations of the State of Polarization,” J. Opt. Soc. Am. 56, 1732-1740 (1966). https://doi.org/10.1364/JOSA.56.001732Google Scholar

Schönhofer, A. and Kuball H. -G., “Symmetry properties of the Mueller matrix,” Chem. Phys. 115, 159-167 (1987). https://doi.org/10.1016/0301-0104(87)80030-7Google Scholar

Gil, J. J., “Components of purity of a 3D polarization state,” J. Opt. Soc. Am. A 33, 40-43 (2016). https://doi.org/10.1364/JOSAA.33.000040Google Scholar

Shurcliff, W. A., [Polarized Light]. Harvard U. Press, Cambridge (1962).Google Scholar

Gil, J. J. and Bernabéu, E., “Depolarization and polarization indices of an optical system,” Opt. Acta 33, 185-189 (1986). https://doi.org/10.1080/713821924Google Scholar

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

Gil, J. J., Ossikovski, R., and San José, I., “Singular Mueller matrices,” J. Opt. Soc Am. A 33, 600-609 (2016). https://doi.org/10.1364/JOSAA.33.000600Google Scholar

Lu, S.-Y. and Chipman, R. A., “Homogeneous and inhomogeneous Jones matrices,” J. Opt. Soc. Am. A 11, 766-773 (1994). https://doi.org/10.1364/JOSAA.11.000766Google Scholar

Gil, J. J., “Structure of polarimetric purity of a Mueller matrix and sources of depolarization,” Opt. Commun. 368, 165-173 (2016). https://doi.org/10.1016/j.optcom.2016.01.092Google Scholar

Gil, J. J. and San José, I., “Two-vector representation of a nondepolarizing Mueller matrix,” arXiv:1602.01794 [physics.optics] (2016). https://doi.org/10.1016/j.optcom.2016.04.043Google Scholar

Gil, J. J., “Components of purity of a Mueller matrix,” J. Opt. Soc. Am. A 28, 1578-1585 (2011). https://doi.org/10.1364/JOSAA.28.001578Google Scholar

Simon, R., “Nondepolarizing systems and degree of polarization,” Opt. Commun. 77, 349-354 (1990). https://doi.org/10.1016/0030-4018(90)90123-BGoogle Scholar

Tudor, T., “Generalized observables in polarization optics,”. J. Phys. A 36, 9577-9590 (2003). https://doi.org/10.1088/0305-4470/36/36/309Google Scholar

Gil, J. J., “Review on Mueller matrix algebra for the analysis of polarimetric measurements,” J. Ap. Remote Sens. 8, 081599-37 (2014). https://doi.org/10.1117/1.JRS.8.081599Google Scholar

Gil, J. J., San José, I. and Ossikovski, R., “Serial-parallel decompositions of Mueller matrices,” J. Opt. Soc. Am. A 30, 32–50 (2013). https://doi.org/10.1364/JOSAA.30.000032Google Scholar

Gil, J. J. and San José, I., “Polarimetric subtraction of Mueller matrices,” J. Opt. Soc. Am. A 30, 1078-1088 (2013). https://doi.org/10.1364/JOSAA.30.001078Google Scholar

San Jose, I. and Gil, J. J., “Invariant indices of polarimetric purity. Generalized indices of purity for nxn covariance matrices,” Opt. Commun. 284, 38-47 (2011). https://doi.org/10.1016/j.optcom.2010.08.077Google Scholar

Gil, J. J., “Transmittance constraints in serial decompositions of depolarizing Mueller matrices. The arrow form of a Mueller matrix,” J. Opt. Soc. Am. A 30, 701-707 (2013). https://doi.org/10.1364/JOSAA.30.000701Google Scholar

Gil, J. J., “From a nondepolarizing Mueller matrix to a depolarizing Mueller matrix,” J. Opt. Soc. Am. A 31, 2736-2743 (2014). https://doi.org/10.1364/JOSAA.31.002736Google Scholar

Ossikovski, R., “Interpretation of nondepolarizing Mueller matrices based on singular-value decomposition,” J. Opt. Soc. Am. A 25, 473-482 (2008). https://doi.org/10.1364/JOSAA.25.000473Google Scholar

Gil, J. J. and San José, I., “Reduced form of a Mueller matrix,” J. Mod Opt. Accepted (2016). https://doi.org/10.1080/09500340.2016.1162337Google Scholar