## Abstract

Decomposition methods have been applied to in-plane Mueller matrix ellipsometric scattering data of the Spectralon reflectance standard. Data were measured at the wavelengths 532 nm and 1500 nm, using an achromatic optimal Mueller matrix scatterometer applying a photomultiplier tube and a high gain InGaAs detector for the two wavelengths. A parametric model with physical significance was deduced through analysis of the product decomposed matrices. It is found that when the data are analyzed as a function of the scattering angle, similar to particle scattering, the matrix elements are largely independent of incidence angle. To the first order, we propose that a Guassian lineshape is appropriate to describe the polarization index, while the decomposed diagonal elements of the retardance matrix have a form resembling Rayleigh single scattering. New models are proposed for the off diagonal elements of the measured Mueller matrix.

© 2013 Optical Society of America

## 1. Introduction

It is important to have a well characterized Mueller matrix of white reflectance standards, such as the Spectralon, in order to obtain correct measurements whenever applying this standard, and in particular, in cases where the polarization of light intervene in the measurement process [1–4]. For un-polarized and linearly polarized illumination, Spectralon white reflectance standards are well documented and often used as calibration standards for optical instruments, but for applications involving other polarization states, more documentation is needed [1]. The Spectralon is the commercial product closest to a Lambertian surface, and is a white diffuse material (based on polytetrafluoroethylene) with excellent reflection properties, produced by Labsphere, USA. It is formed as a thermoplastic resin by heat and pressure treatment. The thermoplastic materials are resin machined into flat, thermally stable, high permittivity surfaces with porosity within the top layer (generally a few tenths of a millimeter thick) that randomizes the phase of the incident light-due to internal multiple reflections-to produce diffuse reflected light [4]. Indeed, the Spectralon has the highest diffuse reflectance values of any known substance (up to 99%) in the ultraviolet, visible, and near-infrared spectral ranges [1, 2], and the reflection properties are spectrally flat between 400 and 1500 nm.

It is further fundamentally interesting to understand how a scattering Mueller matrix from a strongly diffuse medium may be interpreted. It concerns both the interpretation of Mueller matrices resulting from rough surfaces [5–7], powder samples, or extremely diffuse materials such as the Spectralon [2], and the important related problem of Mueller matrices resulting from biological tissue, which has recently regained considerable interest [8–12]. A parametric model of the Mueller matrix elements of the Spectralon diffuser, which relates to basic physical phenomena within the diffuser appears useful in order to develop the understanding of other strongly scattering media, such as biological tissue.

Two recent studies by Svensen et al. [1] and Germer et al. [2] reported the Mueller matrix measurements of the Spectralon using the wavelengths 532 nm [1] and 632.8 nm [2]. In the work of Svensen et al., the Mueller matrix was presented for several incidence angles (0°, 30°, 45° and 60°), and it was shown that the Mueller matrix could to the first order be seen as the sum of an ideal diffuser and a correction matrix [1]. Furthermore, the degree of polarization was observed to be a key parameter in order to develop a parametric model for each of the elements of the Mueller matrix. The origin of the structure in the scattering Mueller matrix elements reported in the work of Svensen et al. was not found, but in a recent study by Germer et al. [2], it was noted that particle scattering is an important process for the Mueller matrix of the Spectralon. Therefore, we have here pursued the idea that the Mueller matrix of the Spectralon appears to belong to a more general class of particle scattering problems.

Considerable effort has been made in order to increase the understanding of the Mueller matrix using various decomposition methods [8,9,13–18]. Decomposition of the Mueller matrix allows to separate different polarization phenomena, and thus to identify the basic building blocks of a complicated Mueller matrix consisting of both dichroism, retardance and diffuse scattering [9]. In the current work, several decomposition methods have been applied to the in-plane Mueller matrix scattering measurements of the Spectralon. As a result, more detailed information about the sample has been extracted, and an improved parametric model for the Mueller matrix has been obtained. In particular, the parametric model was found by studying the data as a function of the Rayleigh scattering angle used in particle scattering experiments, and by using the basis functions revealed from the decomposition.

## 2. Decomposition theory and basis functions

Product or polar decomposition allows the decomposition of a Mueller matrix into a product of elementary Mueller matrices, describing depolarization, retardance and diattenuation [13, 14, 17]. However, a given product decomposition may not always be representing the physical phenomena taking place in the sample, particularly whenever the different polarization altering phenomena are randomly ordered or occur simultaneously. There are several possible product decomposition methods such as the symmetrical decomposition [15], and the more common forward [17] and reverse [13] decompositions which are strongly interrelated [13]. The differential decomposition [12, 16, 19] and the recently described roots decomposition [20, 21] are other possibly more appropriate alternatives. On the other hand, in the work of Svensen et al. [1], a simplified version of the sum (or spectral) decomposition [1, 14] was used to propose a parametrization of the Mueller matrix of the Spectralon.

In the Spectralon sample, many of the scattering phenomena that are altering the polarization, may be envisaged to occur simultaneously, and the imposed ordering of the polarizing phenomena in the product decomposition methods may be incorrect. As such, the forward, reverse and symmetrical polar decomposition are equally (un)valid. Unfortunately, the most natural decomposition method for such a system, the differential decomposition method, is currently strictly only valid for transmission geometries. As a result, the application of a decomposition method is in this work regarded as an empirical mathematical method in order to enhance and reveal an underlying dominating trend in order to parametrize the optical response. We choose therefore to mainly report the application of the product decomposition, due to its numerical stability in the presence of noise, in order to reveal more information about the basis functions describing the scattering Mueller matrix elements.

The measured Mueller matrix (**M**) is, in the forward polar decomposition, decomposed as [17]

**M**

_{Δ}is the depolarization matrix,

**M**

*is the retardance matrix, and*

_{R}**M**

*is the diattenuation matrix. In terms of the forward polar decomposition, the depolarization, retardance and diattenuation matrices are [17]*

_{D}**p**=[

*p*

_{1}

*p*

_{2}

*p*

_{3}]

*and diattenuation vector*

^{T}**d**= [

*d*

_{1}

*d*

_{2}

*d*

_{3}]

*are directly available from the measured Mueller matrix, while the 3 × 3 submatrices*

^{T}**m**

_{Δ},

**m**

*and*

_{R}**m**

*are determined by the polar decomposition formalism [13, 17]. For an ideal diffuser, the*

_{D}**M**

_{Δ}is diagonal and given by

*M*

_{11}·diag(1,

*α*,

*β*,

*γ*), where

*α*,

*β*and

*γ*quantify the degree of polarization of linear, elliptical and circular polarized light, respectively [14, 17]. However, in the work of Lu and Chipman [17] it was proposed that the eigenvalues of

**m**

_{Δ}give the three principal depolarization factors, and hence the average degree of polarization can be defined by [17] The total retardance is calculated from the retardance matrix as [18]

A non-depolarizing matrix can in many cases be made up of the basic building blocks of a diattenuator and a retarder. The block diagonal building block, the so-called isotropic dichroic retarder, commonly encountered in specular reflection from an isotropic surface, or from Mie scattering from small particles is a particularly useful reference

*p̂*,

*ŝ*,

*k̂*} right handed geometry, are given by the normal ellipsometric angles Ψ and Δ [23, 24]

The polarization index *p _{D}* (or degree of purity) [25], which is obtained from the full Mueller matrix was selected as the key figure to model the Spectralon Mueller matrix [1]

*p*≤ 1. When

_{D}*p*= 0, the outgoing light is fully depolarized, and when

_{D}*p*= 1, the outgoing light is fully polarized. This index compares well to the average degree of polarization

_{D}*p*(Eq. (3)).

_{c}## 3. Experimental details

New Mueller matrix data were recorded for the Spectralon in the in-plane scattering geometry shown in Fig. 1. The experimental data are presented according to the right handed coordinate system {*x̂*, *ŷ*, *k̂*} where *x̂* = − *p̂*,*ŷ* =*ŝ* with respect to the right handed system {*ŝ*, *p̂*, *k̂*}, see Fig. 1. All measurements were performed in the incidence plane (see Fig. 1 for the definitions of the angles). The Stokes vector is in this case given as

Mueller matrices for an illumination wavelength of *λ* = 532 nm were measured at every 1° for scattering angles *θ _{S}* (see Fig. 1) relative to the normal of the Spectralon surface from −90° to 90°, with the exception of a dead-zone of ±10° around the back-scattering angle. A rotating diffuser was used prior to the polarization state generator in order to remove speckle effects [1]. Mueller matrices at

*λ*= 1500 nm are reported for every 2° for

*θ*from 10° to 90°, as we will show that the Mueller matrix was found similar to the one at

_{S}*λ*= 532 nm. It is emphasized that the data at

*λ*= 1500 may contain speckle effects, as no rotating diffuser was available for that wavelength. It is noted that all the measured Mueller matrices were found to be physical according to the Cloude filtering criterion [26].

A well conditioned broadband (achromatic) Mueller matrix ellipsometer [27] was used to acquire the experimental data at both wavelengths (see [1] for the experimental details at 532 nm). A 45 mW laser at 1500 nm and an InGaAs detector with variable gain 10^{5} – 10^{9} were used for the near infrared measurement.

## 4. Results and discussion

Figure 2 shows the Mueller matrices recorded at 532 nm for the angles of incidence 0°, 30°, 45°, 60° and 75°, where all elements have been normalized by *M*_{11}. The *M*_{11} element is proportional to the phase function, and is shown in the top left figure. Its form is typical for a near Lambertian diffuser [1,4], although the total intensity has not been calibrated in the current work. It is particularly noted that the Mueller matrix elements have been plotted and analyzed as a function of *θ _{R}*. It is observed from Fig. 2 that the general trend of the normalized Mueller matrix elements appears to be independent of the incidence angle. Deviations from this trend can be observed for light scattered nearly parallel to the sample surface. It is also observed that both block diagonal and off-block diagonal elements are non-zero, particularly for angles close to the sample surface, showing that the Spectralon in these cases is far from an ideal diffuser. The diagonal elements are seen to be strongly dominating for scattering near parallel to the sample surface.

It is interesting to note that the amplitudes of the *m*_{24} ≈ −*m*_{42} elements oscillate with the azimuthal orientation (*χ*) of the sample (i.e. upon rotating the sample around its normal), and that they can be close to zero in a given orientation, see Fig. 3. Such an anisotropic effect was also observed by Germer et al. [2]. To demonstrate smaller variations in the remaining Mueller matrix elements upon sample rotation, we have also plotted Δ*m*_{33} ≡ *m*_{33} − 〈*m*_{33}〉 and Δ(*m*_{34} − *m*_{43}) ≡ (*m*_{34} − 〈*m*_{34}〉) − (*m*_{43} − 〈 *m*_{43}〉), where 〈*m _{ij}*〉 are the averages upon a full 360° rotation, as a function of

*χ*. There may be a weak oscillation in the

*m*

_{34}and

*m*

_{43}elements, while the remaining elements are not observed to oscillate upon azimuthal rotation of the sample, similar to Δ

*m*

_{33}in Fig. 3. It is thus curious that here the linear retardance between the ±45° axis is strongly oscillating, while the linear retardance along the

*x*−

*y*axis is only weakly oscillating, when interpreted using the differential decomposition theory [12, 16]. However, the off-block diagonal elements are non-zero (except possibly

*m*

_{14}and

*m*

_{41}), although only

*m*

_{24}and

*m*

_{42}are observed to be sensitive to the azimuthal orientation of the sample (

*χ*).

Figure 4 shows the recorded Mueller matrix for both wavelengths 1500 nm and 532 nm for the angle of incidence 75°. The Mueller matrix at 1500 nm has a similar form, but with small differences that will be discussed below.

#### 4.1. Model for p_{D}

When the data are plotted against the Rayleigh scattering angle, the depolarization indices *p _{D}* are also observed to become nearly independent of the incidence angles. Previously it was proposed that the following parametric model could represent the data [1]

*p*should be less or equal to one. This is also the case for scattering at a grazing incidence. Therefore, we propose that a peak function, such as the Gaussian function, is a more appropriate choice to fit the data, see Fig. 5. As a result, the following Gaussian function is proposed to model

_{D}*p*

_{D}*p*

_{0}is added. Furthermore, a function

*p*is added in order to express the additional polarization taking place at the backscattering angle

_{BS}*θ*≤

_{i}*θ*and

_{R}*θ*≤ −

_{R}*θ*− 90°. Furthermore, the model does not capture deviations in

_{i}*p*for scattered light nearly parallel to the sample surface, as particularly seen from the data at 532 nm (Fig. 5(a)). Values for the fitted parameters for

_{D}*p*on the entire data set at 532 nm are given in Table 1, and the plot is shown in Fig. 5(a).

_{D}A Gaussian lineshape was also observed for the *p _{D}* data at 1500 nm, as seen in Fig. 5(c) (note that data were here only recorded for positive

*θ*). It shows that

_{R}*p*from 532 nm (dashed line) is less polarized compared to 1500 nm (full line) for scattered light close to the surface plane. For the back scattering direction, 1500 nm is less polarized than 532 nm. Elsewhere, the differences appear small.

_{D}We will now develop an alternative parametric model for the Mueller matrices in Fig. 2 and in Fig. 4, by inspection of the decomposed Mueller matrices using the product decomposition, and looking for a mathematical form of the Mueller elements which may have a common physical interpretation. In particular, the forward polar decomposition allows to inspect the details of the retardance matrix **M*** _{R}*, the diattenuation matrix

**M**

*and the elements of the depolarization matrix*

_{D}**M**

_{Δ}.

#### 4.2. Simplified model for M_{Δ}

The decomposed Mueller matrix was found to result in a close to diagonal **m**_{Δ}. The other off-diagonal elements were therefore for simplicity neglected in the analysis. Let us inspect these diagonal components which we then for simplicity relate directly to *α*, *β* and *γ* of the ”diagonal” depolarization matrix at 532 nm (see Eq. (2)) in Fig. 5(b), where we also neglect the fact that we are not dealing with an ideal diffuser since the polarizance is non-zero [14]. The model in Eq. (11), was used to fit the full datasets for *α*, *β* and *γ* at 532 nm. The resulting parameters from the fits are given in Table 1. Again, it is noted that the parameters describing the back-scattering contribution are only indicative, since we have limited data in this region. The diagonal components *α*, *β* and *γ* of the depolarization matrix at 1500 nm involves another complication, as the measured matrix can be seen to change from a positive to a negative determinant as a function of the Rayleigh scattering angle. We have treated this issue using the recipe described by Ossikovski et al. [28]. However, we observe that at 1500 nm, *γ* has a different shape than at 532 nm, and at 1500 nm we had to use a Lorentzian peak function for a reasonable fit. As a result, the model in Eq. (11), was used to fit the full datasets for *α*, *β* at 1500 nm, while the following Lorentzian function was used to fit *γ*

*α*,

*β*and

*γ*for

*λ*= 1500 nm are given in Table 2.

#### 4.3. Simplified model for M_{R} and M_{D} matrices

Let us now analyze the retardance matrix, **M*** _{R}*. First, inspecting the retardance in Eq. (4), one can observe a similar trend for the scattered light at 532 nm and at 1500 nm, see Fig. 6. In both cases, the retardance increases as a function of the Rayleigh scattering angle

*θ*and crosses 90° after

_{R}*θ*= 90°. There appears to be a definite observable difference between the two wavelengths. The jumps around 95° – 110° in the retardance at 1500 nm, are noise believed to be due to speckles.

_{R}Figures 7(a) and 7(c) shows the normalized *m _{R}*

_{33}and

*m*

_{R}_{44}elements of the retardance matrix

**M**

*for all recorded angles of incidence. Figs. 7(b) and 7(d) shows the retardance data of*

_{R}*m*

_{R}_{34}and

*m*

_{R}_{43}for all recorded angles of incidence. It is interesting to observe that

*m*

_{R}_{33}and

*m*

_{R}_{44}are nearly identical and follow the Rayleigh scattering model for all incidence angles, when analyzed as a function of

*θ*

_{R}*m*

_{R}_{33}and

*m*

_{R}_{44}of approximately

*θ*≈ 7° at 532 nm and

_{B}*θ*≈ 15° at 1500 nm as indicated by the vertical lines in Fig. 7. A similar functional form (without the shift) was used by White [29] and further proposed by Germer et al. [2] for modeling of the Spectralon. The shift

_{B}*θ*, appears to be dispersive (i.e. a funtion of wavelength), although we find that a similar analysis of the Mueller matrix data reported by Germer et al. also gives a shift

_{B}*θ*≈ 7° at 632.8 nm. The shift

_{B}*θ*is suggested related to a pseudo Brewster angle as a result of particle scattering [30]. It is also useful to observe that the total retardance

_{B}*R*in Fig. 6 has a functional form well described by

*R*≈ arccos(

*C*), particularly at 532 nm.

The *m _{R}*

_{34}and

*m*

_{R}_{43}elements, which are normally zero in the Rayleigh scattering model, are proposed to have a related form

*θ*≈ 7° at 532 nm and

_{B}*θ*≈ 15° at 1500 nm. The origin of such a functional form could be the averaging of the Mie scattering from anisotropic shaped particles. On the other hand, one could also imagine that light scattered by particles will travel through a form-birefringent medium causing such a retardance. The retardance appears as the dominant effect, and its effect scales approximately with

_{B}*p*≈

_{D}*α*≈

*β*≈

*γ*in the final Mueller matrix.

Having determined the basic functional forms for the lower right elements of the **m*** _{R}* retardance sub-matrix, we may similarly try to parametrize the remaining elements of

**m**

*, see the experimental data in Fig. 8. Indeed, the function −0.02 + 0.17*

_{R}*R*≈ −

*p*

_{0}/3 +

*R*/6, where

*R*= arccos(

*C*), gives a reasonable fit to the available data sets at 532 nm, for the

*m*

_{R}_{23}element (note that

*m*refers to the retardance matrix

_{Rij}*M*). Similarly, the function 0.02 – 0.23

_{R}*N*≈

*p*

_{0}/3 −

*N*/4 fits the

*m*

_{R}_{24}dataset, while

*m*

_{R}_{22}≈ 1. Adapting the symmetries

*m*

_{R}_{23}≈

*m*

_{R}_{32}and

*m*

_{R}_{24}≈ −

*m*

_{R}_{42}, we obtain a reasonable model for the entire retardance sub-matrix. The experimental decomposed retardance sub-matrix at 532 nm and incidence angles

*θ*= 30°, 45°, 60° and 75° is shown together with the simulations in Fig. 8. This model for the retardance sub-matrix appears reasonable, and a similar model also appears suitable at 1500 nm.

_{i}There are other interesting features in the diattenuation and the polarizance vectors, particularly in the *m*_{12} and *m*_{21} elements, of the original measured Mueller matrix, which appear to scale as
${p}_{D}^{2}\approx {\alpha}^{2}$. Upon searching for the Rayleigh type functional form for the *m*_{12} and *m*_{21} elements, it was observed that

*θ*is found to be similar, although not identical, to the corresponding angle for the retardance matrix. Although both

_{B}*N*,

*S*and

*C*appear as reasonable functional forms to describe the full data sets, there are certain problems with the periodicity of these functions as a function of

*θ*. We therefore believe that a further improved model must properly include other nonlinear effects in both

_{R}*θ*and

_{R}*θ*as a result of refraction, and thus as a function of incidence and the Rayleigh scattering angle. Finally, it is noted that the diattenuation sub-matrix was found to be close to the identity matrix, and any variations therein was neglected in the development of the parametric model.

_{B}#### 4.4. Final parametric model for the Mueller matrix at 532 and 1500 nm

We now return to the original Mueller matrix, and seek a parametric model using the functional forms based on *α*, *β* and *γ*, in addition to *N*, *C* and *S* as determined in the above decomposition analysis. We propose that a simplified parametric model for the Mueller matrix for the Spectralon, can now be obtained using polar product decomposition. First, let
$\mathbf{p}\approx \mathbf{d}\approx {\left[\frac{1}{2}{p}_{0}-4{\alpha}^{2}N,0,0\right]}^{T}$, **m**_{Δ} ≈ diag(*α*, *β*, *γ*), **m*** _{D}* ≈ diag(1, 1, 1) and
${\mathbf{m}}_{R}\approx \left[\begin{array}{ccc}0& 0& 0\\ 0& C& S\\ 0& -S& C\end{array}\right]$, giving the

**M**

_{Δ},

**M**

*and*

_{R}**M**

*matrices in Eq. (2). A direct application of the forward decomposition, i.e. using Eq. (1), will lead to additional terms of the order*

_{D}*α*|

**d**|, |

**d**|

^{2}and smaller, and further

**p**is no longer identical to

**d**as initially assumed. This is an artefact due to the simplifications in the parametrisation of the decomposed matrices. We therefore argue that the observed symmetry of the experimental Mueller matrix must remain, and that the Mueller matrix should also remain physically realizable in the forward scattering direction. As a result, these additional terms can be neglected. The latter approximation gives identical results to simply letting

**p**=

**0**in

**M**

_{Δ}, prior to the matrix multiplication, and a back-substitution of

**p**=

**d**into the final Mueller matrix. The latter approach will result in the following simplified Mueller matrix with a certain physical basis

*C*,

*S*and

*N*are defined as in Eqs. (14), (15) and (16), respectively. We similarly propose a complete parametric model using the above basis functions, that will better fit also the off diagonal elements of the Spectralon Mueller matrices. From trial an error, we parametrize $\mathbf{p}\approx \mathbf{d}\approx {\left[\frac{1}{2}{p}_{0}-4{\alpha}^{2}N,\frac{1}{2}{p}_{0}-2{\alpha}^{2}N,\frac{1}{3}{p}_{0}\right]}^{T}$. We use

**m**

_{Δ}≈ diag(

*α*,

*β*,

*γ*) and

**m**

*≈ diag(1, 1, 1) as above, while the retardance sub-matrix is refined to*

_{D}*R*= arccos(

*C*). Furthermore, the oscillation of the

*m*

_{24}and the −

*m*

_{42}elements have been introduced in the model as a function of

*χ*. Hence, the matrix product of Eq. (1), with

**p**=

**0**, followed with a back-substitution of

**p**=

**d**in the resulting Mueller matrix, result in the following Mueller matrix

*θ*≈ 7°. For 1500 nm, the parameters for the depolarization factors are given in Table 2, and

_{B}*θ*≈ 15°. Note that the

_{B}*m*

_{14}and

*m*

_{41}elements are negligible at 1500 nm, and that a small constant ±

*p*0/4 has been added to the

*m*

_{34}and

*m*

_{43}elements, for a minor adjustment to the experimental data. We have used

*χ*= 30° for the modelling of both wavelengths in Figs. 2 and 4. Similarly, one could propose a small, although more uncertain, oscillation in the

*m*

_{34}and

*m*

_{43}elements.

The Mueller matrix resulting from the parametric model in Eq. (19), is plotted in Fig. 2 together with the experimental data. It is observed, that there is a reasonable correspondence to the range of obtained experimental data. There is a distinct improvement in the correspondence between the model and the experimental results compared to previous work [1], although more parameters have been introduced in the new model. The modelled Mueller matrix for 1500 nm is shown in Fig. 4 together with the data at 532 nm. Only weak differences are observed between the two wavelengths, e.g. *m*_{14} and *m*_{41} have a much smaller offset at 1500 nm, and *m*_{44} is larger at 1500 nm. A self consistency-check was performed by recalculating the polarization index *p _{D}* using Eq. (8) for both wavelengths. This recalculated polarization index was found to be in excellent agreement with both the experimental data and the Gaussian model in Eq. (11), except for the already mentioned deviation for scattering near parallel to the surface.

The parametric model indicate that Mie and Rayleigh scattering appears as an important process for the Spectralon reflectance standard. It is interesting that the complex form of the retardance matrix can be derived from the basis functions *N*, *C* and *S*. However, as with Rayleigh particles, little direct information about the particles themselves are projected in the current model, except for possibly the shift angle *θ _{B}*, which may be regarded as a pseudo-Brewster angle, and thus be related to either the average material properties or the optical properties of the scatterers [30]. The apparent dispersion of

*θ*supports the latter hypothesis. It is striking that the

_{B}*m*

_{12}and the

*m*

_{21}elements are apparently scaling with $\sim {p}_{D}^{2}$, possibly indicating a multiple scattering process. There are experimental uncertainties with respect to the constant offsets in the measured Mueller matrix (particularly

*m*

_{13}and

*m*

_{14}), but we have not been able to relate these to any evident experimental artefacts. As in [1], the

*m*

_{14}and

*m*

_{41}elements have been modelled by constants.

Finally, to reduce the number of free parameters, the result in Eq. (19) can be put into a reduced form

*p*≈

_{D}*p*≈

_{C}*α*≈

*β*≈

*γ*and

*p*

_{0}has been neglected. It is now clearly observed that the lower right 3 × 3 matrix is basically the retardance sub-matrix multiplied with the average depolarisation factor. To complete the discussion, we can now write the matrix as the sum of an ideal diffuser scaling with the average degree of depolarization, and a correction matrix scaling with the average degree of polarization

It appears interesting to study other well determined diffuse samples to investigate whether the above product decomposition method and in particular, the parametric model found for the Spectralon, may be applied to other strongly scattering systems, such as rough surfaces and powder samples. The exact dependency of incidence angle is expected to be revealed with a combination of more data and improved models. As a result, it is envisaged that more information about the scatterers may be extracted from the normalized Mueller matrix.

Preliminary studies of the current dataset, comparing the product decomposition and the current formulation of the differential decomposition [12, 16], give highly similar results for *α*, *β*, *γ* and the linear retardance in the range 0 ≤ *θ _{R}* ≤ 90°, while it is interesting to observe that the current formulation of the differential decomposition method must be revised for

*θ*> 90°. Indeed, this is as expected since the current form of the differential decomposition is only valid for forward scattering. A revised version of the differential decomposition method [12, 16] could be a more appropriate tool for analysis of the scattering Mueller matrix for the Spectralon.

_{R}Finally, it has been proposed to further refine/validate the parametric model for the off-diagonal elements using a few higher order terms in the actual Mie solution. This approach would definitively result in a model based on stronger physical ground, but is out of the scope of the current work.

## 5. Conclusion

A new parametric model for the Mueller matrix of a Spectralon white reflectance standard diffuser has been determined through the use of polar (product) decomposition techniques. New Mueller matrix data at the wavelengths 532 nm and 1500 nm are presented and analyzed. It was found particularly useful to study the Spectralon in terms of the standard ”particle” or Rayleigh scattering angle, which made the depolarization index and the Mueller matrix elements close to independent of the incidence angle. The basic building blocks in the new model is a Gaussian parametric model for the depolarization index, multiplied with Rayleigh like functional forms deduced from the decomposed matrices. A pseudo Brewster angle appeared dispersive and characteristic of the material. The use of decomposition techniques appears promising in order to reveal the basic functional building blocks of a given scattering Mueller matrix, and it is envisaged that the methodology presented here can be used to study other strongly diffuse samples. Finally, a well characterized Mueller matrix of the Spectralon makes it also useful as a calibration standard for polarization sensitive systems.

## Acknowledgment

L.M.S. Aas and J. Maria acknowledge support from The Norwegian Research Center for Solar Cell Technology (project num. 193829).

## References and links

**1. **Ø. Svensen, M. Kildemo, J. Maria, J. J. Stamnes, and Ø. Frette, “Mueller matrix measurements and modeling pertaining to spectralon white reflectance standards,” Opt. Express **20**, 15045–15053 (2012) [CrossRef] [PubMed] .

**2. **T. A. Germer and H. J. Patrick, “Mueller matrix bidirectional reflectance distribution function measurements and modeling of diffuse reflectance standards,” Proc. SPIE **8160**, 81600D (2011) [CrossRef] .

**3. **D. A. Haner, B. T. McGuckin, and C. J. Bruegge, “Polarization characteristics of spectralon illuminated by coherent light,” Appl. Opt. **38**, 6350–6356 (1999) [CrossRef] .

**4. **A. Bhandari, B. Hamre, Ø. Frette, L. Zhao, J. J. Stamnes, and M. Kildemo, “Bidirectional reflectance distribution function of spectralon white reflectance standard illuminated by incoherent unpolarized and plane-polarized light,” Appl. Opt. **50**, 2431–2442 (2011) [CrossRef] [PubMed] .

**5. **T. A. Leskova, P. A. Letnes, A. A. Maradudin, T. Nordam, and I. Simonsen, “The scattering of light from two-dimensional randomly rough surfaces,” Proc. SPIE **8172**, 817208 (2011).

**6. **I. Simonsen, A. A. Maradudin, and T. A. Leskova, “Scattering of electromagnetic waves from two-dimensional randomly rough penetrable surfaces,” Phys. Rev. Lett. **104**, 223904 (2010) [CrossRef] [PubMed] .

**7. **I. Simonsen, A. A. Maradudin, and T. A. Leskova, “Scattering of electromagnetic waves from two-dimensional randomly rough perfectly conducting surfaces: The full angular intensity distribution,” Phys. Rev. A **81**, 013806 (2010) [CrossRef] .

**8. **N. Ghosh, J. Soni, M. Wood, M. Wallenberg, and I. Vitkin, “Mueller matrix polarimetry for the characterization of complex random medium like biological tissues,” Pramana **75**, 1071–1086 (2010) [CrossRef] .

**9. **N. Ghosh and I. A. Vitkin, “Tissue polarimetry: concepts, challenges, applications, and outlook,” J. Biomed. Opt. **16**, 110801 (2011) [CrossRef] [PubMed] .

**10. **A. Pierangelo, S. Manhas, A. Benali, M. R. Antonelli, T. Novikova, P. Validire, B. Gayet, and A. De Martino, “Use of mueller polarimetric imaging for the staging of human colon cancer,” Proc. SPIE **7895**, 78950E (2011) [CrossRef] .

**11. **V. V. Tuchin, L. V. Wang, and D. A. Zimnyakov, *Optical Polarization in Biomedical Applications* (Berlin Heidelberg, 2006).

**12. **S. Kumar, H. Purwar, R. Ossikovski, I. A. Vitkin, and N. Ghosh, “Comparative study of differential matrix and extended polar decomposition formalisms for polarimetric characterization of complex tissue-like turbid media,” Journal of Biomedical Optics **17**, 105006–105006 (2012) [CrossRef] [PubMed] .

**13. **R. Ossikovski, A. D. Martino, and S. Guyot, “Forward and reverse product decompositions of depolarizing mueller matrices,” Opt. Lett. **32**, 689–691 (2007) [CrossRef] [PubMed] .

**14. **R. Ossikovski, M. Anastasiadou, S. Ben Hatit, E. Garcia-Caurel, and A. De Martino, “Depolarizing mueller matrices: how to decompose them?” Phys. Status Solidi (a) **205**, 720–727 (2008) [CrossRef] .

**15. **R. Ossikovski, “Analysis of depolarizing mueller matrices through a symmetric decomposition,” J. Opt. Soc. Am. A **26**, 1109–1118 (2009) [CrossRef] .

**16. **R. Ossikovski, “Differential matrix formalism for depolarizing anisotropic media,” Opt. Lett. **36**, 2330–2332 (2011) [CrossRef] [PubMed] .

**17. **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] .

**18. **S. Manhas, M. K. Swami, P. Buddhiwant, N. Ghosh, P. K. Gupta, and J. Singh, “Mueller matrix approach for determination of optical rotation in chiral turbid media in backscattering geometry,” Opt. Express **14**, 190–202 (2006) [CrossRef] [PubMed] .

**19. **N. Ortega-Quijano and J. L. Arce-Diego, “Mueller matrix differential decomposition.” Opt. lett. **36**, 1942–4 (2011) [CrossRef] [PubMed] .

**20. **H. D. Noble and R. a. Chipman, “Mueller matrix roots algorithm and computational considerations.” Opt. Express **20**, 17–31 (2012) [CrossRef] [PubMed] .

**21. **H. Noble, S. McClain, and R. Chipman, “Mueller matrix roots depolarization parameters,” Appl. Optics **51**, 735–744 (2012) [CrossRef] .

**22. **C. Bohren and D. Huffman, *Absorption and Scattering of Light by Small Particles Interscience* (J. Wiley & Sons, New York, 1983).

**23. **P. Hauge, R. Muller, and C. Smith, “Conventions and formulas for using the mueller-stokes calculus in ellipsometry,” Surf. Sci. **96**, 81–107 (1980) [CrossRef] .

**24. **R. Azzam and N. Bashara, *Ellipsometry and Polarized Light* (North Holland, 1977).

**25. **J. J. Gil and E. Bernabeu, “A depolarization criterion in mueller matrices,” Opt. Acta **32**, 259–261 (1985) [CrossRef] .

**26. **S. R. Cloude, “Conditions for the physical realisability of matrix operators in polarimetry,” Proc. SPIE **1166**, 177–185 (1989) [CrossRef] .

**27. **F. Stabo-Eeg, M. Kildemo, I. S. Nerbø, and M. Lindgren, “Well-conditioned multiple laser mueller matrix ellipsometer,” Opt. Eng. **47**, 073604–073604–9 (2008) [CrossRef] .

**28. **R. Ossikovski, M. Anastasiadou, and A. D. Martino, “Product decompositions of depolarizing mueller matrices with negative determinants,” Opt. Commun. **281**, 2406–2410 (2008) [CrossRef] .

**29. **R. L. White, “Polarization in reflection nebulae. I - Scattering properties of interstellar grains,” Astrophys. J. **229**, 954–961 (1979) [CrossRef] .

**30. **A. Erbe, K. Tauer, and R. Sigel, “Ellipsometric light scattering for the characterization of thin layers on dispersed colloidal particles,” Phys. Rev. E **73**, 031406 (2006) [CrossRef] .