Next Article in Journal
High-Speed Spiral-Phase Donut-Modes-Based Hybrid FSO-MMF Communication System by Incorporating OCDMA Scheme
Previous Article in Journal
High-Sensitivity Temperature Sensor Based on the Perfect Metamaterial Absorber in the Terahertz Band
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Extended Representation of Mueller Matrices

1
Instituto Aragonés de Estadística, Gobierno de Aragón, Bernardino Ramazzini 5, 50015 Zaragoza, Spain
2
Departments de Física Aplicada, Universidad de Zaragoza, Pedro Cerbuna 12, 50009 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
Photonics 2023, 10(1), 93; https://doi.org/10.3390/photonics10010093
Submission received: 7 December 2022 / Revised: 9 January 2023 / Accepted: 11 January 2023 / Published: 14 January 2023

Abstract

:
The so-called arbitrary decomposition of a given Mueller matrix into a convex sum of nondepolarizing constituents provides a general framework for parallel decompositions of polarimetric interactions. Even though arbitrary decomposition can be performed through an infinite number of sets of components, the nature of such components is subject to certain restrictions which limit the interpretation of the Mueller matrix in terms of simple configurations. In this communication, a new approach based on the addition of some portion of a perfect depolarizer before the parallel decomposition is introduced, leading to a set of three components which depend, respectively, on the first column, the first row, and the remaining 3 × 3 submatrix of the original Mueller matrix, so that those components inherit, in a decoupled manner, the polarizance vector, the diattenuation vector, and the combined complementary polarimetric information on depolarization and retardance.

1. Introduction

Mueller polarimetry nowadays has an immense variety of applications in science, engineering, remote sensing, etc. Once a Mueller matrix M is measured, the analysis of the information provided by the sixteen real elements of M constitutes a critical stage, in which different procedures are applied, including the identification of descriptors of diattenuation, polarizance, retardance, and depolarization, as well as methods for the serial and parallel decompositions of M into simpler components. In particular, so-called arbitrary and characteristic decompositions [1,2,3] constitute powerful tools that provide sets of peculiar parallel constituents that are susceptible to specific interpretations.
Nevertheless, the scope of potential components with prescribed features is limited by mathematical restrictions that rely on the fact that the coherency matrix C [4], a positive semidefinite Hermitian matrix that is biunivocally associated with the given M (that is, M determines C unambiguously and vice versa, see Equation (4)), does not always satisfy the property rank C = 4 , but its rank can take integer values in the interval 1 rank C 4 [5,6]. Recall that rank C equals the number of independent parallel components of M [1,3].
In this work, a new composition-decomposition approach is presented in which, prior to the application of a specific parallel decomposition, the measured M is combined with a portion of a perfect depolarizer, of which the associated normalized Mueller matrix has the diagonal form M ^ Δ 0 = diag ( 1 , 0 , 0 , 0 ) [6]. Due to the simple structure of M ^ Δ 0 , which corresponds to polarimetric white noise [7], the composed matrix inherits all the anisotropies exhibited by M. Analogously to what occurs in some image diagnostic techniques, where certain contrast agents or colorants are added to the sample in order to improve the images, in the present case, a fully isotropic element is added with the aim of extending the scope of decompositions, which allows certain featured representations to be obtained.
The decompositions dealt with in this communication are focused on identifying sets of three constituents of which the Mueller matrices have the simplest possible structure, regardless of the fact that they are, in general, depolarizing.
The structure of the present communication is organized as follows. The main notions necessary for the development of the new approach are presented in Section 2; Section 3 describes the formulation of the parallel composition of a given Mueller matrix M and a perfect depolarizer; Section 4 deals with the introduction of the homogeneous extended form of M and its decomposition, where the constituents exhibit equal mean intensity coefficients and are analyzed from the mathematical, geometric, and physical points of view; Section 5 describes an alternative extended decomposition of M where the constituents, being of similar nature to those of the homogeneous extended decomposition, exhibit different mean intensity coefficients, and have even simpler forms; Section 6 is devoted to the analysis of the extended decompositions for the particular case in which M lacks either polarizance or diattenuation; and the obtained results are discussed in Section 7.

2. Theoretical Background

The transformation of polarized light by the action of a linear medium (under fixed interaction conditions) can always be represented mathematically as s = M s where s and s are the Stokes vectors that represent the states of polarization of the incident and emerging light beams, respectively, whereas M is the Mueller matrix associated with this kind of interaction, which can always be expressed as [8,9,10]:
M = m 00 M ^ , M ^ ( 1 D T P m ) [ m 1 m 00 ( m 11 m 12 m 13 m 21 m 22 m 23 m 31 m 32 m 33 ) D ( m 01 , m 02 , m 03 ) T m 00 , P ( m 10 , m 20 , m 30 ) T m 00 ]
where m i j   ( i , j = 0 , 1 , 2 , 3 ) denote the elements of M; the superscript T indicates transpose; m 00 is the mean intensity coefficient (MIC), i.e., the ratio between the intensity of the emerging light and the intensity of incident unpolarized light; D and P are the diattenuation and polarizance vectors, with absolute values D (diattenuation) and P (polarizance); and m is the normalized 3 × 3 submatrix associated with M, which provides the complementary information on retardance and depolarization properties.
Regarding the ability of M to preserve the degree of polarization (DOP) of totally polarized incident light, a proper measure is given by the degree of polarimetric purity of M (also called the depolarization index) [11], P Δ , which can be expressed as
P Δ = D 2 + P 2 + 3 P S 2 3 = 2 P P 2 3 + P S 2 [ P P D 2 + P 2 2 ] ,
where P P is the so-called degree of polarizance, or enpolarizance, and P S is the polarimetric dimension index (also called the degree of spherical purity), defined as [6,12,13]
P S m 2 3 [ m 2 1 m 00 k , l = 1 3 m k l 2 ] ,
with m 2 being the Frobenius norm of m.
The maximal degree of polarimetric purity, P Δ = 1 , is exhibited uniquely by nondepolarizing (or pure) media (i.e., media that do not decrease the degree of polarization of totally polarized incident light), whereas P Δ = 0 is characteristic of perfect depolarizers, with the associated Mueller matrix M Δ 0 = m 00 diag ( 1 , 0 , 0 , 0 ) . The maximal value of P S , P S = 1 , implies P Δ = 1 with P P = 0 (pure and nonenpolarizing media), which corresponds uniquely to retarders (regardless of the value of m 00 , i.e., regardless of whether they are transparent or exhibit certain amount of isotropic attenuation), the minimal polarimetric dimension index, P S = 0 , corresponds to media exhibiting m = 0 . The maximal enpolarizance, P P = 1 , implies P Δ = 1 and corresponds to perfect polarizers, whereas the minimal, P P = 0 , is exhibited by nonenpolarizing interactions (either pure or depolarizing) [12,14].
In general, two kinds of decompositions of a Mueller matrix can be performed, namely, serial decompositions (through products of Mueller matrices) and parallel decompositions (through weighted sums of Mueller matrices) [6]. Furthermore, both decompositions can be combined, leading to serial-parallel decompositions [15].
Parallel decompositions, under the scope of which the new approach is developed, consist of representing a Mueller matrix as a convex sum of Mueller matrices. The physical meaning of parallel decompositions is that the incoming electromagnetic wave splits into a set of pencils that interact, without overlapping, with a number of material components that are spatially distributed in the illuminated area, and the emerging pencils are incoherently recombined into the emerging beam.
Thus, the concept of parallel (or additive) composition of Mueller matrices underlies the very concept of the Mueller matrix and obeys certain specific rules; in particular, the coefficients of the Mueller components in the sum should be positive and should add up to one (convex sum) [1,3]. This property is directly linked to the covariance criterion for Mueller matrices, namely, given a Mueller matrix M, its associated Hermitian coherency matrix C ( M ) is positive and semidefinite. The explicit expression of C ( M ) , in terms of the elements m i j of M, is [4]:
C ( M ) = 1 4 ( m 00 + m 11 + m 22 + m 33 m 01 + m 10 i ( m 23 m 32 ) m 02 + m 20 + i ( m 13 m 31 ) m 03 + m 30 i ( m 12 m 21 ) m 01 + m 10 + i ( m 23 m 32 ) m 00 + m 11 m 22 m 33 m 12 + m 21 + i ( m 03 m 30 ) m 13 + m 31 i ( m 02 m 20 ) m 02 + m 20 i ( m 13 m 31 ) m 12 + m 21 i ( m 03 m 30 ) m 00 m 11 + m 22 m 33 m 23 + m 32 + i ( m 01 m 10 ) m 03 + m 30 + i ( m 12 m 21 ) m 13 + m 31 + i ( m 02 m 20 ) m 23 + m 32 i ( m 01 m 10 ) m 00 m 11 m 22 + m 33 )
The passivity constraint (natural linear polarimetric interactions do not amplify the intensity of light) is completely characterized by the inequality m 00 ( 1 + Q ) 1 [16,17], where Q max ( D , P ) . Thus, m 00 1 / ( 1 + Q ) , and therefore media exhibiting nonzero polarizance or diattenuation (called enpolarizing media) necessarily feature m 00 < 1 , whereas the limit m 00 = 1 corresponds to transparent Mueller matrices, which have the general form
M ( 1 0 T 0 m )
By combining the covariance and passivity criteria, physical Mueller matrices are characterized by the ensemble criterion, which means that a given 4 × 4 real matrix X is a Mueller matrix if and only if it can be expressed as convex sum of pure and passive Mueller matrices, which is equivalent to saying that C ( X ) (which, by construction, is a Hermitian matrix) is positive and semidefinite (i.e., the four eigenvalues of C ( X ) are nonnegative), and, in addition, X satisfies the passivity condition x 00 ( 1 + Q ) 1 [16].
Some additional concepts and descriptors that will be useful for further developments and discussions are briefly reviewed below.
Leaving aside systems exhibiting magneto-optic effects, the Mueller matrix that represents the same linear interaction as M, but with the incident and emergent directions of the light probe interchanged, is given by [18,19]:
M r = diag ( 1 , 1 , 1 , 1 ) M T diag ( 1 , 1 , 1 , 1 ) ,
consequently, D ( M r ) = P ( M ) and P ( M r ) = D ( M ) , showing that D and P share a common nature related to the ability of the medium to enpolarize (i.e., to increase the degree of polarization of) unpolarized light incoming in either the forward or reverse directions [1,20]. Note that D, P, and other quantities considered below (when applied to the reverse Mueller matrix), which are defined based on square averages of some Mueller matrix elements, are insensitive to magneto-optic effects (which only affect the signs of certain elements of M).
A fundamental concept related to Mueller matrices is their classification in terms of the auxiliary matrix N G M T G M (with G diag ( 1 , 1 , 1 , 1 ) ). If N is diagonalizable (i.e., there exists an invertible matrix A such that A 1 N A is diagonal), then M can be written in the type-I normal form [21,22,23,24,25,26,27,28]
M = M J 2 M Δ d M J 1 ,
where M J 1 and M J 2 are pure Mueller matrices, and M Δ d (called the type-I canonical Mueller matrix) has the form
M Δ d = diag ( d 0 , d 1 , d 2 , ε d 3 ) [ 0 d 3 d 2 d 1 d 0 ε det M / | det M | ]
M Δ d is a diagonal Mueller matrix representing an intrinsic depolarizer [6,10]. Observe that pure Mueller matrices are always of type-I, in which case M Δ d coincides with the identity matrix [6].
On the other hand, when N is not diagonalizable, M is type-II and it can always be written in the type-II normal form [28]
M = M J 2 M Δ n d M J 1 M Δ n d ( 2 a 0 a 0 0 0 a 0 0 0 0 0 0 a 2 0 0 0 0 a 2 ) [ 0 a 2 a 0 ]
where M J 1 and M J 2 are nonsingular pure Mueller matrices, and M Δ n d is called the type-II canonical Mueller matrix [28].
Leaving aside the MIC, the complete physical information contained in a generic Mueller matrix M can be represented geometrically by means of the pair of ellipsoids E Δ P and E Δ D generated by M and M r , respectively. The canonical depolarizer M Δ (with M Δ representing either M Δ d or M Δ n d , depending on whether M is type-I or type-II) is fully characterized by its associated canonical ellipsoid E Δ . The use of the three characteristic ellipsoids E Δ P , E Δ D , and E Δ leads to a complete and significant geometric view of the properties of M [29].
Consider now the following modified singular value decomposition of the submatrix m of M,
m = m R O m A m R I [ m R i 1 = m R i T , det m R i = + 1 ( i = I , O ) m A diag ( a 1 , a 2 , ε a 3 ) , a 1 a 2 a 3 0 , ε det m / | det m | ]
where the nonnegative parameters ( a 1 , a 2 , a 3 ) are the singular values of m (taken in decreasing order), so that
M R i = ( 1 0 T 0 m R i ) ( i = I , O )
are orthogonal Mueller matrices (representing respective transparent retarders). The arrow form M A ( M ) of M is then defined as [30]
M A ( M ) M R O T M M R I T = m 00 ( 1 D A T P A m A ) [ m A m R O T m m R I T = diag ( a 1 , a 2 , ε a 3 ) a 1 a 2 a 3 0 , ε det m / | det m | D A = m R I D , P A = m R O T P ]
and contains up to ten nonzero elements. The corresponding arrow decomposition of M is defined as
M = M R O M A M R I
Observe that the diattenuation and polarizance vectors of M are recovered from those of M A through the respective transformations D = m R I T D A and P = m R O P A , which preserve the absolute values of the transformed vectors and are determined by the entrance and exit retarders M R I and M R O of M.

3. Parallel Compositions of a Given Mueller Matrix and That of a Perfect Depolarizer

In order to enable decompositions of a Mueller matrix M for which the components exhibit certain essential properties of M in a decoupled and simple manner, the extended form of M is built by adding to M an appropriate proportion of a perfect depolarizer M Δ 0 .
Given a Mueller matrix M (depolarizing or not), it is always possible to build the depolarizing Mueller matrix
M = q m 00 M ^ + ( 1 q ) m 00 M ^ Δ 0 = m 00 ( 1 q D T q P q m ) , ( 0 < q < 1 )
This transformation can also be expressed as follows in terms of the coherency matrices C C ( M ) , C ^ C ( M ^ ) , and C ^ Δ 0 C ( M ^ Δ 0 ) = Ι / 4 (with I being the identity matrix),
C = q m 00 C ^ + ( 1 q ) m 00 ( 1 / 4 ) I , ( 0 < q < 1 )
Obviously, rank C ^ Δ 0 = 4 , and, since q > 0 , together with the fact that the rank of a sum of positive semidefinite Hermitian matrices is greater than or equal to the rank of the addend with largest rank, the coherency matrix C necessarily satisfies rank C = 4 .
Consequently, depending on the value of q, which represents the portion of M with respect to the whole composed matrix, the resulting matrix admits certain parallel decompositions which are not realizable for M itself. In particular, as will be shown in the next section, q = 1 / 3 constitutes a critical value that ensures the realizability of a parallel decomposition of M into three components with very simple structures depending on P, D, and m, respectively, regardless of the value of rank C . Such a decomposition is not possible, in general, when q > 1 / 3 , whereas q < 1 / 3 leads to similar parallel decompositions but with an additional component that is proportional to M ^ Δ 0 with the respective coefficient 1 3 q .

4. Homogeneous Extended Decomposition of a Mueller Matrix

In this section, we consider general Mueller matrices with nonzero polarizance and diattenuation ( ( P > 0 , D > 0 ) . The cases where P = 0 or D = 0 will be dealt with in Section 6.
Taking q = 1 / 3 in Equation (14), the resulting matrix is
M H ( M ) = 1 3 ( m 00 M ^ ) + 2 3 ( m 00 M ^ Δ 0 ) = m 00 ( 1 D T / 3 P / 3 m / 3 )
which will be called the homogeneous extended form of M, where the term homogeneous is used to note that the MICs ( m 00 ) of both components are equal to that of M, which is consistent with the name coined for arbitrary decompositions where all components have equal MICs ( m 00 ) [3].
It is straightforward to prove that M H ( M ) can always be expressed as the following parallel combination of three Mueller matrices
M H ( M ) = 1 3 M m + 1 3 M P + 1 3 M D , M m m 00 ( 1 0 T 0 m ) , M P m 00 ( 1 0 T P 0 ) , M D m 00 ( 1 D T 0 0 ) ,
where m, P, and D appear isolated within respective components. The above decomposition will be called the homogeneous extended decomposition of M, which should be interpreted as a parallel composition of media represented by M m , M P , and M D with equal MICs ( m 00 ) and respective portions (or cross sections) equal to 1/3 (i.e., the intensity I of the incident light probe is shared among the three components with equal intensities I/3).
When applied to the arrow form M A ( M ) , and in accordance with Equation (12), the homogeneous extended decomposition takes the simplified form:
M H ( M A ) = 1 3 M m A + 1 3 M P A + 1 3 M D A [ M m A = M R O T M m M R I T = m 00 diag ( 1 , a 1 , a 2 , ε a 3 ) M P A = m 00 M R O T ( 1 0 T P 0 ) = m 00 ( 1 0 T m R O T P 0 ) = m 00 ( 1 0 T P A 0 ) M D A = m 00 ( 1 D T 0 0 ) M R I T = m 00 ( 1 D T m R I T 0 0 ) = m 00 ( 1 D A T 0 0 ) ]
In this case, the physical information held by M is parameterized through the following sixteen parameters: the MIC, m 00 , of M; the three angular parameters determining the entrance retarder M R I ; the three angular parameters determining the exit retarder M R O ; the polarizance vector P of M; the diattenuation vector D of M; and the three diagonal elements ( a 1 , a 2 , ε a 3 ) of m m A .
The specific properties of each of the components of the homogeneous extended decomposition are analyzed in the following subsections.

4.1. Nonenpolarizing Component

The 3 × 3 submatrix m of the nonenpolarizing component, M m , coincides with that of M. Thus, by considering the procedure used to define the arrow form M A ( M ) [30], m can be written as in Equation (10), m = m R O m A m R I , with m A = diag ( a 1 , a 2 , ε a 3 ) , and therefore M m can be expressed through the following dual retarder transformation [31] (which coincides with the normal form of M m ):
M m = m 00 M R O M ^ Δ m M R I , [ M ^ Δ m = diag ( 1 , a 1 , a 2 , ε a 3 ) ]
where M ^ Δ m is the normalized Mueller matrix of an intrinsic depolarizer, while the entrance and exit equivalent retarders, M R I and M R O , respectively, coincide with those of M.
The eigenvalues (in decreasing order) of the coherency matrix C m associated with M m are
λ m 0 = m 00 ( 1 + a 1 + a 2 + | a 3 | ) , λ m 1 = m 00 ( 1 + a 1 a 2 | a 3 | ) , λ m 2 = m 00 ( 1 a 1 + a 2 | a 3 | ) , λ m 3 = m 00 ( 1 a 1 a 2 + | a 3 | ) .
Note that the above eigenvalues coincide with the diagonal elements of the coherency matrix C associated with M. Furthermore, since M is a Mueller matrix, C is positive and semidefinite, and therefore its diagonal elements are necessarily nonnegative; consequently, the eigenvalues λ m i of the Hermitian matrix C m are nonnegative, showing that C m is a proper coherency matrix and consequently M m is a Mueller matrix.
The integer parameter r m rank C m can take values in the interval 1 r m 4 depending on the nature of the interaction represented by M. Thus, from the point of view of the arbitrary decomposition [3], M m has a number r m of parallel pure components. Moreover, since P P ( M m ) = 0 , the only source of polarimetric purity of M m is P S ( M m ) = P Δ ( M m ) = P S ( M ) .
The passivity of the starting M entails the passivity of M m , that is, m 00 ( 1 + Q ) 1   m 00 1 , where Q max ( D , P ) . Observe also that Equation (19) shows that M m is always a type-I Mueller matrix.
It has been shown that any Mueller matrix admits a meaningful geometric representation by means of the three corresponding characteristic ellipsoids, namely, the forward and reverse ellipsoids, together with the canonical ellipsoid [29]. In the case of M m , its characteristic ellipsoids adopt the simple forms shown in Figure 1.
All three characteristic ellipsoids are centered (the center of the ellipsoid coincides with the origin of the Poincaré sphere) and have the same shape, with semiaxes a 1 , a 2 , a 3   ( a 1 a 2 a 3 ) . The symmetry axes of the canonical ellipsoid E Δ ( M m ) are aligned to the respective axes S 1 , S 2 , S 3 of the Poincaré sphere. Both forward and reverse ellipsoids, E Δ P ( M m ) and E Δ D ( M m ) , are rotated by the respective effects of M R I and M R O .
Degenerate cases occur when a 3 = 0 with a 2 > 0 (the characteristic ellipsoids become ellipses), a 2 = 0 with a 1 > 0 (straight segments), and a 1 = 0 (single point), which corresponds to the trivial limiting case of the zero Mueller matrix.

4.2. Enpolarizing Components

As for the components M P and M D , their associated coherency matrices C P and C D have the following respective sets of eigenvalues:
C P { λ P + = m 00 ( 1 + P ) / 4 ( d o u b l e ) λ P = m 00 ( 1 P ) / 4 ( d o u b l e ) C D { λ D + = m 00 ( 1 + D ) / 4 ( d o u b l e ) λ D = m 00 ( 1 D ) / 4 ( d o u b l e )
which are nonnegative because the conditions P 1 and D 1 are always satisfied by any Mueller matrix M [11].
Thus, C P can have either two or four parallel pure components (two when P = 1 and four when P < 1 ), and the same happens for C D depending on whether D = 1 (two parallel pure components) or D < 1 (four parallel pure components).
As occurs with component M m , the passivity of M implies the passivity of both M P and M D .
Moreover, both M P and M D are type-I singular Mueller matrices, which is evident when they are expressed as the following respective serial decompositions:
M P = M P M R M ^ Δ 0 , M D = M ^ Δ 0 M R M D , M P = m 00 ( 1 P T P m P ) , M D = m 00 ( 1 D T D m D ) , m K 1 K 2 I 3 + 1 1 K 2 K 2 K K T [ K = P , D K = P , D ] ,
where M P and M D represent respective normal diattenuators defined based on vectors P and D [10,32,33,34], M R represents an arbitrary retarder (which plays no role in the definitive forms of M P and M D ), stands for the Kronecker product, and I 3 is the 3 × 3 identity matrix.
Since M P and M D are singular, their characteristic ellipsoids are given by single points [35]. In particular, as shown in Figure 2 and Figure 3, the canonical and reverse ellipsoids of M P , as well as the forward and canonical ellipsoids of M D , are single points located in the center of the Poincaré sphere, whereas the forward ellipsoid of M P and the reverse ellipsoid of M D are given by single points located within or on the surface of the Poincaré sphere depending on whether K < 1 or K = 1   ( K = P , D ) , respectively.

5. Extended Decomposition of a Mueller Matrix

As in Section 4, we consider here general Mueller matrices with nonzero polarizance and diattenuation ( P > 0 , D > 0 ) . The cases where P = 0 or D = 0 will be dealt with in Section 6.
As an alternative to the homogeneous extended form of M, the MIC of the perfect depolarizer added to M can be taken as m 00 ( P + D ) / 2 (instead of m 00 ), so that the extended form of M is defined as
M E ( M ) = 1 3 m 00 M ^ + 2 3 [ m 00 ( P + D ) 2 ] M ^ Δ 0 = 1 3 m 00 ( 1 + P + D D T P m ) [ P > 0 , D > 0 ]
Note that, since P , D 1 , the MIC of the perfect depolarizer, m 00 ( P + D ) / 2 , is smaller than that of the perfect depolarizer in the homogeneous extended form, m 00 .
The extended decomposition of M is then defined as
M E ( M ) = 1 3 M m + 1 3 M P ^ + 1 3 M D ^ , M m m 00 ( 1 0 T 0 m ) , M P ^ m 00 P ( 1 0 T P ^ 0 ) , M D ^ m 00 D ( 1 D ^ T 0 0 ) , [ P ^ = P / P , D ^ = D / D ] .
Here, unlike the homogeneous extended decomposition, the polarizing and diattenuating components exhibit respective polarizance and diattenuation vectors, P ^ and D ^ , having necessarily maximal absolute values | P ^ | = | D ^ | = 1 , and consequently the coherency matrices C P ^ = C ( M P ^ ) and C D ^ = C ( M D ^ ) have the following respective sets of eigenvalues:
C P ^ { λ P + = m 00 / 2 ( d o u b l e ) λ P = 0 ( d o u b l e ) C D ^ { λ D + = m 00 / 2 ( d o u b l e ) λ D = 0 ( d o u b l e )
so that rank C P ^ = rank C D ^ = 2 .
In accordance with this property, both the forward ellipsoid of M P ^ and the reverse ellipsoid of M D ^ are given by single points that are located on the surface of the Poincaré sphere.
The structure of the components of the extended decomposition is even simpler than that of the homogeneous one. In particular, M P ^ and M D ^ can be interpreted by means of the following respective two-component parallel compositions:
M P ^ = 1 2 m 00 P ( 1 P ^ T P ^ P ^ × P ^ T ) + 1 2 m 00 P ( 1 P ^ T P ^ P ^ × P ^ T ) , M D ^ = 1 2 m 00 D ( 1 D ^ T D ^ D ^ × D ^ T ) + 1 2 m 00 D ( 1 D ^ T D ^ D ^ × D ^ T ) ,
where the nondepolarizing components of M P ^ and M D ^ correspond to respective perfect polarizers (i.e., media that fully polarize light entering either in the forward or reverse directions).
Note also that, as with the homogeneous extended decomposition, the passivity of the three components is ensured by the passivity of M itself.
The extended decomposition of the arrow form M A ( M ) has the simple form
M E ( M A ) = 1 3 M m A + 1 3 M P ^ A + 1 3 M D ^ A , [ M m A = M R O T M m M R I T = m 00 diag ( 1 , a 1 , a 2 , ε a 3 ) , M P ^ A = ( m 00 P ) M R O T ( 1 0 T P ^ 0 ) = ( m 00 P ) ( 1 0 T m R O T P ^ 0 ) = ( m 00 P ) ( 1 0 T P ^ A 0 ) , M D ^ A = ( m 00 D ) ( 1 D ^ T 0 0 ) M R I T = ( m 00 D ) ( 1 D ^ T m R I T 0 0 ) = ( m 00 D ) ( 1 D ^ A T 0 0 ) . ]

6. Extended Decompositions of Matrices Lacking Polarizance or Diattenuation

The extended decompositions introduced in Section 4 and Section 5 apply to Mueller matrices with nonzero polarizance and diattenuation ( P > 0 and D > 0 ). When P = 0 or D = 0 , the required ratio between the perfect depolarizer M Δ 0 and the starting M for a well-defined extended form of M, without an excess of M Δ 0 , changes with respect to the case where P > 0 and D > 0 . Appropriate particular forms of extended decompositions are analyzed below, which, together those dealt with in Section 4 and Section 5, provide a complete case analysis.
When P > 0 and D = 0 (nondiattenuating Mueller matrices, denoted as M D = 0 ), the extended decompositions of M take the following forms, where the denominator of the coefficients (polarimetric cross sections) equals the number of components (two):
M H ( M D = 0 ) = 1 2 m 00 M ^ + 1 2 m 00 M ^ Δ 0 = m 00 ( 1 0 T P / 2 m / 2 ) = 1 2 M m + 1 2 M P , M E ( M D = 0 ) = 1 2 m 00 M ^ + 1 2 ( m 00 P ) M ^ Δ 0 = 1 2 m 00 ( 1 + P 0 T P m ) = 1 2 M m + 1 2 M P ^ .
Analogously, when P = 0 and D > 0 (nonpolarizing Mueller matrices, denoted as M P = 0 ),
M H ( M P = 0 ) = 1 2 m 00 M ^ + 1 2 m 00 M ^ Δ 0 = m 00 ( 1 D T / 2 0 m / 2 ) = 1 2 M m + 1 2 M D , M E ( M P = 0 ) = 1 2 m 00 M ^ + 1 2 ( m 00 D ) M ^ Δ 0 = m 00 ( ( 1 + D ) / 2 D T / 2 0 m / 2 ) = 1 2 M m + 1 2 M D ^ .
The Mueller matrices M m , M P , M P ^   M D   M D ^ of the components in the above equations have the forms defined in Equations (17) and (24).
In the case of a nonenpolarizing Mueller matrix M m   ( P = D = 0 ) , it directly has the form of the nonenpolarizing component and therefore its extended form coincides with itself (i.e., the corresponding coefficient for the added perfect depolarizer is zero).

7. Discussion

The extended representation of a given Mueller matrix M involves its convex sum with a perfect depolarizer M Δ 0 , which does not exhibit any anisotropy (or polarimetric preference). Consequently, the anisotropies inherited by the extended representations M H ( M ) and M E ( M ) are precisely those of M. Furthermore, for any given M, both M H ( M ) and M E ( M ) are biunivocally related to M through simple expressions. In fact, the m, P, and D of M H ( M ) are none other than those of M, whereas the structure of M E ( M ) is given by m,  P ^ , D ^ , P, and D (with P = P P ^ and D = D D ^ ).
Furthermore, the extended forms M H ( M ) and M E ( M ) in Equations (16) and (23) involving M ^ Δ 0 , of which the associated coherency matrix is C Δ 0 = I / 4 , ensure that the respective coherency matrices C H ( M ) and C E ( M ) satisfy rank C H ( M ) =   rank C E ( M ) = 4 , which extends substantially the scope of parallel decompositions of the extended forms in comparison to those directly applicable to M, and makes it possible to obtain the extended decompositions considered in Section 4, Section 5 and Section 6.
It is remarkable that all Mueller matrices M m , M P , M D , M P ^ , and M D ^ of the constituents of the extended forms of M are type-I, and are therefore free from the intricate structure exhibited by type-II Mueller matrices [28,36].
Due to their simplicity, the extended forms and decompositions of any given Mueller matrix M can be straightforwardly interpreted regardless of its structural complexity. When the arrow form of M is considered, the analysis becomes even simpler because in that case, once the entrance and exit retarders have been decoupled from M, M m becomes diagonal [37].
An interesting limiting situation is that corresponding to pure Mueller matrices, of which the general structure adopts the form [38]
M J = M R L 2 M R L 0 ( m 00 M ^ D L 0 ) M R L 0 M R L 1 ,
where M R L 1 and M R L 2 are equivalent (entrance and exit) linear retarders (each depending on two angular parameters), M R L 0 represents a horizontal linear retarder (depending on a single parameter), and M ^ D L 0 is the normalized Mueller matrix of a normal [32,33,34] horizontal linear diattenuator which only depends on the polarizance-diattenuation D of M J (recall that pure Mueller matrices necessarily satisfy P = D [39]).
M ^ D L 0 = ( 1 D 0 0 D 1 0 0 0 0 1 D 2 0 0 0 0 1 D 2 )
Therefore, the extended components of the arrow form m 00 M ^ D L 0 of M J are given by
M m = m 00 ( 1 0 0 0 0 1 0 0 0 0 1 D 2 0 0 0 0 1 D 2 ) , M P = m 00 ( 1 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 ) , M D = m 00 ( 1 D 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ) , M P ^ = m 00 ( 1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 ) , M D ^ = m 00 ( 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ) .
Coming back to general Mueller matrices (depolarizing or not), observe that since the extended representation applies to any Mueller matrix without restriction, the following question may arise: given an arbitrary set of normalized matrices of the form M ^ m , M ^ P , and M ^ D , is the 4 × 4 matrix X ^ retrieved as X ^ = 3 M H 2 M ^ Δ 0 from M ^ H = ( 1 / 3 ) M ^ m +   ( 1 / 3 ) M ^ P + ( 1 / 3 ) M D in Equation (17) a (normalized) Mueller matrix? The answer is negative, as evidenced by taking, for instance, the following set of matrices:
M ^ m ( 1 0 0 0 0 x 0 0 0 0 x 0 0 0 0 x ) ( 0 x < 1 ) , M ^ P ( 1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 ) , M ^ D ( 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ) ,
for which X ^ does not satisfy the covariance conditions (i.e., it is not a Mueller matrix) because its associated coherency matrix C ( X ^ ) has, at least, a negative eigenvalue.

8. Conclusions

The Mueller matrices obtained through certain parallel combinations of a given Mueller matrix M and a perfect depolarizer are always susceptible to be submitted to respective kinds of parallel decompositions named the homogeneous extended decomposition and the extended decomposition, the components of which have very simple structures which are directly inherited from the anisotropies exhibited by M.
The parallel composition of M and a perfect depolarizer, with appropriate convex coefficients, only affects the MIC (mean intensity coefficient) of the resulting composed matrix, but it is the key for M to be interpreted in terms of the properties of the components of the corresponding extended decompositions. In particular, two components can be straightforwardly determined from the polarizance and diattenuation vectors of M, respectively, whereas the third component depends exclusively on the 3 × 3 submatrix m of M, which encompasses the remaining polarimetric information. That is to say, once the information on polarizance and diattenuation has been decoupled and allocated to respective parallel components, the structure of the remaining nonenpolarizing component allows for the recovery of the complete polarimetric information (including the depolarization and retardance properties) held by M.
In summary, any Mueller matrix M (depolarizing or nondepolarizing) is susceptible to being represented, uniquely, through the extended representations M H ( M ) and M E ( M ) , and admits respective extended decompositions where the structural properties of M appear decoupled in a very simple manner and encoded into separate components.

Author Contributions

Conceptualization, I.S.J.; and J.J.G.; methodology, I.S.J.; and J.J.G.; writing—original draft preparation, I.S.J.; and J.J.G.; writing—review and editing, I.S.J.; and J.J.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Gil, J.J. Polarimetric characterization of light and media-Physical quantities involved in polarimetric phenomena. Eur. Phys. J. Appl. Phys. 2007, 40, 1–47. [Google Scholar] [CrossRef] [Green Version]
  2. Gil, J.J. Parallel decompositions of Mueller matrices and polarimetric subtraction. EPJ Web Conf. 2010, 5, 04002. [Google Scholar] [CrossRef] [Green Version]
  3. Gil, J.J.; San José, I. Arbitrary decomposition of a Mueller matrix. Opt. Lett. 2019, 44, 5715–5718. [Google Scholar] [CrossRef]
  4. Cloude, S.R. Group theory and polarization algebra. Optik 1986, 75, 26–36. [Google Scholar]
  5. Gil, J.J.; San José, I. Polarimetric subtraction of Mueller matrices. J. Opt. Soc. Am. A 2013, 30, 1078–1088. [Google Scholar] [CrossRef] [PubMed]
  6. Gil, J.J.; Ossikovski, R. Polarized Light and the Mueller Matrix Approach, 2nd Ed; CRC Press: Boca Raton, FL, USA, 2022. [Google Scholar]
  7. Gil, J.J. On optimal filtering of measured Mueller matrices. Appl. Opt. 2016, 55, 5449–5455. [Google Scholar] [CrossRef] [Green Version]
  8. Robson, B.A. The Theory of Polarization Phenomena; Clarendon Press: Oxford, UK, 1975. [Google Scholar]
  9. Xing, Z.-F. On the deterministic and non-deterministic Mueller matrix. J. Mod. Opt. 1992, 39, 461–484. [Google Scholar] [CrossRef]
  10. Lu, S.-Y.; Chipman, R.A. Interpretation of Mueller matrices based on polar decomposition. J. Opt. Soc. Am. A 1996, 13, 1106–1113. [Google Scholar] [CrossRef]
  11. Gil, J.J.; Bernabéu, E. Depolarization and polarization indices of an optical system. Opt. Acta 1986, 33, 185–189. [Google Scholar] [CrossRef]
  12. Gil, J.J. Components of purity of a Mueller matrix. J. Opt. Soc. Am. A 2011, 28, 1578–1585. [Google Scholar] [CrossRef]
  13. Gil, J.J. Thermodynamic reversibility in polarimetry. Photonics 2022, 9, 650. [Google Scholar] [CrossRef]
  14. Gil, J.J. Structure of polarimetric purity of a Mueller matrix and sources of depolarization. Opt. Commun. 2016, 368, 165–173. [Google Scholar] [CrossRef] [Green Version]
  15. Gil, J.J.; San José, I.; Ossikovski, R. Serial–Parallel decompositions of Mueller matrices. J. Opt. Soc. Am. A 2013, 30, 32–50. [Google Scholar] [CrossRef]
  16. Gil, J.J. Characteristic properties of Mueller matrices. J. Opt. Soc. Am. A 2000, 17, 328–334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. San José, I.; Gil, J.J. Characterization of passivity in Mueller matrices. J. Opt. Soc. Am. A 2020, 37, 199–208. [Google Scholar] [CrossRef] [Green Version]
  18. Sekera, Z. Scattering Matrices and Reciprocity Relationships for Various Representations of the State of Polarization. J. Opt. Soc. Am. 1966, 56, 1732–1740. [Google Scholar] [CrossRef]
  19. Schönhofer, A.; Kuball, H.-G. Symmetry properties of the Mueller matrix. Chem. Phys. 1987, 115, 159–167. [Google Scholar] [CrossRef]
  20. Ossikovski, R.; De Martino, A.; Guyot, S. Forward and reverse product decompositions of depolarizing Mueller matrices. Opt. Lett. 2007, 32, 689–691. [Google Scholar] [CrossRef]
  21. Van der Mee, C.V.M. An eigenvalue criterion for matrices transforming Stokes parameters. J. Math. Phys. 1993, 34, 5072–5088. [Google Scholar] [CrossRef]
  22. Sridhar, R.; Simon, R. Normal form for Mueller matrices in polarization optics. J. Mod. Opt. 1994, 41, 1903–1915. [Google Scholar] [CrossRef]
  23. Bolshakov, Y.; van der Mee, C.V.M.; Ran, A.C.M.; Reichstein, B.; Rodman, L. Polar decompositions in finite dimensional indefinite scalar product spaces: Special cases and applications. In Operator Theory: Advances and Applications; Gohberg, I., Ed.; Birkhäuser: Basel, Switzerland, 1996; Volume 87, pp. 61–94. [Google Scholar]
  24. Bolshakov, Y.; van der Mee, C.V.M.; Ran, A.C.M.; Reichstein, B.; Rodman, L. Errata for: Polar decompositions in finite dimensional indefinite scalar product spaces: Special cases and applications. Integral Equ. Oper. Theory 1997, 27, 497–501. [Google Scholar] [CrossRef]
  25. Gopala Rao, A.V.; Mallesh, K.S.; Sudha. On the algebraic characterization of a Mueller matrix in polarization optics. I. Identifying a Mueller matrix from its N matrix. J. Mod. Opt. 1998, 45, 955–987. [Google Scholar]
  26. Gopala Rao, A.V.; Mallesh, K.S.; 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. 1998, 45, 989–999. [Google Scholar]
  27. Ossikovski, R. Analysis of depolarizing Mueller matrices through a symmetric decomposition. J. Opt. Soc. Am. A 2009, 26, 1109–1118. [Google Scholar] [CrossRef] [PubMed]
  28. Ossikovski, R. Canonical forms of depolarizing Mueller matrices. J. Opt. Soc. Am. A 2010, 27, 123–130. [Google Scholar] [CrossRef] [PubMed]
  29. Ossikovski, R.; Gil, J.J.; San José, I. Poincaré sphere mapping by Mueller matrices. J. Opt. Soc. Am. A 2013, 30, 2291–2305. [Google Scholar] [CrossRef]
  30. Gil, J.J. Transmittance constraints in serial decompositions of Mueller matrices. The arrow form of a Mueller matrix. J. Opt. Soc. Am. A 2013, 30, 701–707. [Google Scholar] [CrossRef]
  31. Gil, J.J. Invariant quantities of a Mueller matrix under rotation and retarder transformations. J. Opt. Soc. Am. A 2016, 33, 52–58. [Google Scholar] [CrossRef] [Green Version]
  32. Lu, S.-Y.; Chipman, R.A. Homogeneous and inhomogeneous Jones matrices. J. Opt. Soc. Am. A 1994, 11, 766–773. [Google Scholar] [CrossRef]
  33. Tudor, T.; Gheondea, A. Pauli algebraic forms of normal and non-normal operators. J. Opt. Soc. Am. A 2007, 24, 204–210. [Google Scholar] [CrossRef] [Green Version]
  34. Tudor, T. Interaction of light with the polarization devices: A vectorial Pauli algebraic approach. J. Phys. A-Math. Theor 2008, 41, 41530. [Google Scholar] [CrossRef]
  35. Gil, J.J.; Ossikovski, R.; San José, I. Singular Mueller matrices. J. Opt. Soc. Am. A 2016, 33, 600–609. [Google Scholar] [CrossRef] [Green Version]
  36. San José, I.; Gil, J.J.; Ossikovski, R. Algorithm for the numerical calculation of the serial components of the normal form of depolarizing Mueller matrices. Appl. Opt. 2020, 59, 2291–2297. [Google Scholar] [CrossRef] [PubMed]
  37. Gil, J.J.; San José, I. Universal Synthesizer of Mueller Matrices Based on the Symmetry Properties of the Enpolarizing Ellipsoid. Symmetry 2021, 13, 983. [Google Scholar] [CrossRef]
  38. Gil, J.J.; San José, I. Invariant quantities of a nondepolarizing Mueller matrix. J. Opt. Soc. Am. A 2016, 133, 1307–1312. [Google Scholar] [CrossRef] [PubMed]
  39. Gil, J.J. Determination of Polarization Parameters in Matricial Representation. Theoretical Contribution and Development of an Automatic Measurement Device. Ph.D. Dissertation, University of Zaragoza, Zaragoza, Spain, 1983. Available online: http://zaguan.unizar.es/record/10680/files/TESIS-2013-057.pdf (accessed on 5 December 2022).
Figure 1. Characteristic ellipsoids of the nonenpolarizing component M m of the extended form of a Mueller matrix M: (a) forward ellipsoid E Δ P ( M m ) ; (b) canonical ellipsoid E Δ ( M m ) , of which the semiaxes a 1 , a 2 , a 3 (with a 1 a 2 a 3 ) are aligned with the Poincaré axes S 1 , S 2 , S 3 , respectively; and (c) reverse ellipsoid E Δ D ( M m ) .
Figure 1. Characteristic ellipsoids of the nonenpolarizing component M m of the extended form of a Mueller matrix M: (a) forward ellipsoid E Δ P ( M m ) ; (b) canonical ellipsoid E Δ ( M m ) , of which the semiaxes a 1 , a 2 , a 3 (with a 1 a 2 a 3 ) are aligned with the Poincaré axes S 1 , S 2 , S 3 , respectively; and (c) reverse ellipsoid E Δ D ( M m ) .
Photonics 10 00093 g001
Figure 2. The characteristic ellipsoids of the polarizing component M P of the homogeneous extended form of a Mueller matrix M are given by single points (the three semiaxes are zero-valued): (a) forward ellipsoid, E Δ P ( M P ) , which lies on the surface of the Poincaré sphere if and only if P = 1 ; (b) canonical ellipsoid E Δ ( M P ) ; and (c) reverse ellipsoid E Δ D ( M P ) .
Figure 2. The characteristic ellipsoids of the polarizing component M P of the homogeneous extended form of a Mueller matrix M are given by single points (the three semiaxes are zero-valued): (a) forward ellipsoid, E Δ P ( M P ) , which lies on the surface of the Poincaré sphere if and only if P = 1 ; (b) canonical ellipsoid E Δ ( M P ) ; and (c) reverse ellipsoid E Δ D ( M P ) .
Photonics 10 00093 g002
Figure 3. The characteristic ellipsoids of the polarizing component M D of the homogeneous extended form of a Mueller matrix M are given by single points (the three semiaxes are zero-valued): (a) forward ellipsoid E Δ P ( M D ) ; (b) canonical ellipsoid E Δ ( M D ) ; and (c) reverse ellipsoid, E Δ P ( M D ) , which lies on the surface of the Poincaré sphere if and only if D = 1 .
Figure 3. The characteristic ellipsoids of the polarizing component M D of the homogeneous extended form of a Mueller matrix M are given by single points (the three semiaxes are zero-valued): (a) forward ellipsoid E Δ P ( M D ) ; (b) canonical ellipsoid E Δ ( M D ) ; and (c) reverse ellipsoid, E Δ P ( M D ) , which lies on the surface of the Poincaré sphere if and only if D = 1 .
Photonics 10 00093 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

San José, I.; Gil, J.J. Extended Representation of Mueller Matrices. Photonics 2023, 10, 93. https://doi.org/10.3390/photonics10010093

AMA Style

San José I, Gil JJ. Extended Representation of Mueller Matrices. Photonics. 2023; 10(1):93. https://doi.org/10.3390/photonics10010093

Chicago/Turabian Style

San José, Ignacio, and José J. Gil. 2023. "Extended Representation of Mueller Matrices" Photonics 10, no. 1: 93. https://doi.org/10.3390/photonics10010093

APA Style

San José, I., & Gil, J. J. (2023). Extended Representation of Mueller Matrices. Photonics, 10(1), 93. https://doi.org/10.3390/photonics10010093

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop