Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Degree of polarization of light scattered from correlated surface and bulk disorders

Open Access Open Access

Abstract

Using a single scattering theory, we derive the expression of the degree of polarization of the light scattered from a layer exhibiting both surface and volume scattering. The expression puts forward the intimate connection between the degree of polarization and the statistical correlation between surface and volume disorders. It also permits a quantitative analysis of depolarization for uncorrelated, partially correlated and perfectly correlated disorders. We show that measuring the degree of polarization could allow one to assess the surface-volume correlation function, and that, reciprocally, the degree of polarization could be engineered by an appropriate design of the correlation function.

© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Polarimetric measurements are key elements in the toolbox for the characterization of complex photonic structures, including thin films, metamaterials, photonic crystals, plasmonic gratings [13], or disordered materials such as colloidal suspensions [46] and rough surfaces [79]. Polarization analysis is also of great interest for systems displaying both surface and volume disorder [10,11]. In this context, depolarization measurements have shown their ability to discriminate between surface and bulk scattering. The approach has been implemented on highly scattering samples [1214], where multiple scattering from the bulk is the main source of depolarization. Interestingly, depolarization can also reveal information on weakly scattering systems, where the interaction with light occurs chiefly through single scattering, and in which volume and surface disorders may contribute with similar weights. It is often assumed that single scattering does not produce depolarization, which is actually not true for systems exhibiting (at least) two types of disorders with different polarization responses [15]. Examples of such systems are clouds of particles of different species [16], media with depolarizing dielectric heterogeneities [17], dielectric films with rough interfaces [18], pseudo-random overlaid gratings [19,20], or samples with a rough surface and volume dielectric fluctuations [21,22]. Recently, perfect depolarization has even been predicted in the single scattering regime, for a system combining uncorrelated surface and volume disorders [15].

An open question is whether depolarization of the light scattered by a system with surface and volume disorders can provide information on the existence of statistical correlations between the two types of disorder. The purpose of this paper is to examine this question in the case of weakly disordered samples, in which surface and volume disorders contribute through single scattering. To proceed, we establish a general relation between the degree of polarization of the scattered light and the cross-correlation function between the surface roughness and the dielectric fluctuations in the volume. Based on this relation, we address several issues, such as the conditions to get full depolarization of the incident light, or the possibility to engineer the surface-volume correlation function to produce a prescribed value of the degree of polarization of the scattered light.

The paper is organized as follows. In section 2, we introduce the geometry and the statistical model, focusing on the description of the cross-correlation function between the surface and volume disorders. In section 3, we summarize the scattering theory that was described initially in Ref. [23], and derive the expression of the degree of polarization. Based on this expression, we examine in section 4 the general conditions to get depolarization of the scattered light. In section 5, we analyze the behavior of the degree of polarization for correlated surface-volume scattering. In particular, we discuss the possibility of maximizing depolarization, and of designing the surface-volume cross-correlation function to reach a prescribed form of the degree of polarization of the scattered light. Finally, we summarize the main results in section 6.

2. Scattering geometry and statistical model

We consider a scattering layer with average thickness $L$ separating two semi-infinite media, and exhibiting both surface and volume disorders [Fig. 1(a)]. We take direction $x_3$ to be normal to the layer which is assumed to be of infinite extent along directions $x_1$ and $x_2$. The layer has a rough upper surface, described by a profile $x_3 = \zeta (\mathop {\mathbf {x}_{\parallel }^{}})$, with $\mathop {\mathbf {x}_{\parallel }^{}}=(x_1,x_2)$. Its lower interface is flat, and coincides with the plane $x_3=-L$. The external upper and lower media, corresponding to the regions $x_3 > \zeta (\mathop {\mathbf {x}_{\parallel }^{}})$ and $x_3 < -L$, have real dielectric functions $\varepsilon _1$ and $\varepsilon _2$, respectively. The layer also exhibits volume disorder, described by a dielectric function $\varepsilon (\mathop {\mathbf {x}_{}^{}})=\varepsilon _2+\Delta \varepsilon (\mathop {\mathbf {x}_{}^{}})$ fluctuating around the average value $\varepsilon _2$. The geometry is depicted in Fig. 1.

 figure: Fig. 1.

Fig. 1. (a) Cross-section of the scattering layer in the plane $(x_1,x_3)$, showing both a rough surface and volume dielectric fluctuations. (b) Schematics defining the incident and scattering wave vectors. The incident wave vector lies in the plane $(x_1,x_3)$, with an in-plane component $\mathop {\mathbf {p}_{0}^{}}$ and a direction defined by the polar angle of incidence $\theta _0$. The scattered wave vector has an arbitrary in-plane component $\mathop {\mathbf {p}_{}^{}}$, and a direction defined by the scattering angles $(\theta,\phi )$.

Download Full Size | PDF

In this study we will focus on the role of statistical correlations between the rough surface and the bulk dielectric fluctuations on depolarization. Depending on the dependence of $\Delta \varepsilon$ with respect to the longitudinal direction $x_3$, different types of layers can be defined. Here we will consider dielectric fluctuations taking constant values across the layer, with $\Delta \varepsilon$ depending only on $\mathop {\mathbf {x}_{\parallel }^{}}$. This type of disorder corresponds to the picture in Fig. 1(a), and was referred to as surface-like configuration in Ref. [23]. A real life example of dielectric layers with a rough top free interface and quasi-two-dimensional dielectric fluctuations occurs naturally through the phase segregation of barium borosilicate glass thin films. Upon annealing different phases appear and deform the top surface due to their largely different viscosity and the confining geometry for the film [24,25].

In order to define the statistical model, we start by writing the dielectric function of the whole system in the form

$$\varepsilon (\mathop{\mathbf{x}_{}^{}}) = \varepsilon_1 + \mathrm{H} \Big( \zeta(\mathop{\mathbf{x}_{{\parallel}}^{}}) - x_3 \Big) \Big( \varepsilon_2 -\varepsilon_1 + \Delta \varepsilon (\mathop{\mathbf{x}_{{\parallel}}^{}}) \, \mathrm{H} (x_3 + L) \Big) \: ,$$
where $\mathrm {H}$ is the Heaviside step function. The surface profile $\zeta$ and the dielectric fluctuation $\Delta \varepsilon$ are assumed to be realizations of correlated, zero mean and stationary Gaussian stochastic processes. In these conditions, the stochastic process defining the dielectric function $\varepsilon$ is fully characterized by $\langle \zeta (\mathop {\mathbf {x}_{\parallel }^{}}) \rangle = 0$, $\langle \Delta \varepsilon (\mathop {\mathbf {x}_{\parallel }^{}}) \rangle = 0$ and
$$\begin{aligned} \;\;\;\;\left\langle \zeta (\mathop{\mathbf{x}_{{\parallel}}^{}}) \zeta (\mathop{\mathbf{x}_{{\parallel}}^{\prime}}) \right\rangle = \sigma_\zeta^2 \, W_\zeta (\mathop{\mathbf{x}_{{\parallel}}^{}}-\mathop{\mathbf{x}_{{\parallel}}^{\prime}}) \: , \end{aligned}$$
$$\begin{aligned} \left\langle \Delta \varepsilon (\mathop{\mathbf{x}_{{\parallel}}^{}}) \Delta \varepsilon (\mathop{\mathbf{x}_{{\parallel}}^{\prime}}) \right\rangle = \sigma_\varepsilon^2 \, W_\varepsilon (\mathop{\mathbf{x}_{{\parallel}}^{}}-\mathop{\mathbf{x}_{{\parallel}}^{\prime}}) \: , \end{aligned}$$
$$\begin{aligned} \;\;\;\;\;\;\;\left\langle \zeta (\mathop{\mathbf{x}_{{\parallel}}^{}}) \Delta \varepsilon (\mathop{\mathbf{x}_{{\parallel}}^{\prime}}) \right\rangle = \sigma_\zeta \sigma_\varepsilon \, W_{\zeta \varepsilon} (\mathop{\mathbf{x}_{{\parallel}}^{}}-\mathop{\mathbf{x}_{{\parallel}}^{\prime}}) \: , \end{aligned}$$
where the angle brackets denote an ensemble average. Equations (2a) and (2b) define the surface and dielectric autocorrelation functions $W_\zeta$ and $W_\varepsilon$, and standard deviations $\sigma _\zeta \geq 0$ and $\sigma _\varepsilon \geq 0$. Equation (2c) defines the cross-correlation function of the processes $\zeta$ and $\Delta \varepsilon$. The full definition of the stochastic processes requires explicit expressions for $W_\zeta$, $W_\varepsilon$ and $W_{\zeta \varepsilon }$. A convenient model, introduced in Ref. [23], assumes Gaussian autocorrelation functions given by
$$\begin{aligned}W_\zeta (\mathop{\mathbf{x}_{{\parallel}}^{}}-\mathop{\mathbf{x}_{{\parallel}}^{\prime}}) = \exp \Bigg( - \frac{\big| \mathop{\mathbf{x}_{{\parallel}}^{}} -\mathop{\mathbf{x}_{{\parallel}}^{\prime}} \big|^2}{\ell_\zeta^2} \Bigg) \: , \end{aligned}$$
$$\begin{aligned}W_\varepsilon (\mathop{\mathbf{x}_{{\parallel}}^{}}-\mathop{\mathbf{x}_{{\parallel}}^{\prime}}) = \exp \Bigg( - \frac{\big| \mathop{\mathbf{x}_{{\parallel}}^{}} -\mathop{\mathbf{x}_{{\parallel}}^{\prime}} \big|^2}{\ell_\varepsilon^2} \Bigg) \: , \end{aligned}$$
where $\ell _\zeta$ and $\ell _\varepsilon$ are the correlation lengths of the surface roughness and the dielectric volume fluctuations, respectively. The cross-correlation function can be modeled via a power spectral density of the form
$$\widetilde{W}_{\zeta \varepsilon} (\mathop{\mathbf{p}_{}^{}}) = \gamma(\mathop{\mathbf{p}_{}^{}}) \, \widetilde{W}_\zeta^{1/2}(\mathop{\mathbf{p}_{}^{}}) \widetilde{W}_\varepsilon^{1/2}(\mathop{\mathbf{p}_{}^{}}) \: ,$$
where $\widetilde {f}(\mathop {\mathbf {p}_{}^{}})$ denotes the two-dimensional Fourier transform of a function $f(\mathop {\mathbf {x}_{\parallel }^{}})$. This specific form of the cross-spectral power density is consistent with the constraints imposed by the nature of the covariance matrix, that has to be real, symmetric and positive definite [23]. The factor $\gamma (\mathop {\mathbf {p}_{}^{}})$, which will be denoted by spectral correlation modulator, has to satisfy $|\gamma | \leq 1$ and $\gamma (-\mathop {\mathbf {p}_{}^{}}) = \gamma ^*(\mathop {\mathbf {p}_{}^{}})$ [26].

3. Degree of polarization in the single scattering regime

Our purpose is to connect the degree of polarization of the light backscattered (reflected) from the scattering layer, upon illumination by a monochromatic plane wave with angular frequency $\omega$ incident from medium 1. The complex amplitude of the incident plane wave is taken of the form

$$\mathop{\mathbf{E}_{0}^{}} (\mathop{\mathbf{x}_{}^{}}) = \sum_{\nu = p,s } E_{0,\nu} \, \mathop{\mathbf{\hat{e}}_{1,\nu}^{-}} (\mathop{\mathbf{p}_{0}^{}}) \, \exp [ i \mathop{\mathbf{p}_{0}^{}} \cdot \mathop{\mathbf{x}_{{\parallel}}^{}} - i \alpha_1(\mathop{\mathbf{p}_{0}^{}}) x_3 ] \: ,$$
where $\alpha _1(\mathop {\mathbf {p}_{0}^{}})$ is the normal component of the wave vector, and $\mathop {\mathbf {\hat {e}}_{1,\nu }^{-}} (\mathop {\mathbf {p}_{0}^{}})$ are unit vectors defining the $s$ and $p$ polarizations. These vectors are defined in medium $j=1$ and medium $j=2$ by the following relations
$$\begin{aligned}\alpha_j(\mathop{\mathbf{p}_{}^{}}) = \left( \varepsilon_j k_0^2 - \mathop{\mathbf{p}_{}^{2}} \right)^{1/2} \:, \: \mathrm{Re} (\alpha_j) \geq 0, \: \mathrm{Im} (\alpha_j) \geq 0 \: , \end{aligned}$$
$$\begin{aligned} \mathop{\mathbf{\hat{e}}_{j,s}^{{\pm}}} (\mathop{\mathbf{p}_{}^{}}) = \mathop{\mathbf{\hat{e}}_{3}^{}} \times \mathop{\mathbf{\hat{p}}_{}^{}} \: , \end{aligned}$$
$$\begin{aligned} \mathop{\mathbf{\hat{e}}_{j,p}^{{\pm}}} (\mathop{\mathbf{p}_{}^{}}) = \frac{\pm \alpha_j(\mathop{\mathbf{p}_{}^{}}) \mathop{\mathbf{\hat{p}}_{}^{}} - |\mathop{\mathbf{p}_{}^{}}| \, \mathop{\mathbf{\hat{e}}_{3}^{}}}{\sqrt{\varepsilon_j} k_0} \: .\end{aligned}$$

In these relations $\mathop {\mathbf {\hat {p}}_{}^{}}=\mathop {\mathbf {p}_{}^{}}/|\mathop {\mathbf {p}_{}^{}}|$, $\mathop {\mathbf {\hat {e}}_{3}^{}}$ is the unit vector along the positive $x_3$ axis and $k_0=\omega /c=2\pi /\lambda$ with $c$ the speed of light in vacuum. The meaning of the different wave vectors and polarization vectors is illustrated in Fig. 2.

 figure: Fig. 2.

Fig. 2. Definition of the wave vectors for the incident and scattered fields, and of the unit vectors defining the $s$ and $p$ polarization components.

Download Full Size | PDF

The purpose of this work is to characterize the degree of polarization of the scattered field $\mathop {\mathbf {E}_{s}^{}}(\mathop {\mathbf {x}_{\parallel }^{}},x_3)$ for an observation point in reflection, i.e. for $x_3 > \zeta (\mathop {\mathbf {x}_{\parallel }^{}})$. The Fourier transform of the field with respect to $\mathop {\mathbf {x}_{\parallel }^{}}$ can be written in the form

$$\mathop{\mathbf{\widetilde{E}}_{s}^{}} (\mathop{\mathbf{p}_{}^{}},x_3) = \sum_{\mu = p,s } \mathop{\mathbf{\hat{e}}_{1,\mu}^+} (\mathop{\mathbf{p}_{}^{}}) \sum_{\nu = p,s} R_{\mu \nu}^{}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \, E_{0,\nu} \: \exp \Big( i \alpha_1(\mathop{\mathbf{p}_{}^{}}) \, x_3 \Big) \: ,$$
where the reflection amplitude $R_{\mu \nu }(\mathop {\mathbf {p}_{}^{}},\mathop {\mathbf {p}_{0}^{}})$ connects a scattered wave in state $(\mathop {\mathbf {p}_{}^{}},\mu )$ to an incident wave in state $(\mathop {\mathbf {p}_{0}^{}},\nu )$. In the single scattering regime, the scattered field can be written as the sum of a contribution from the rough surface and a contribution from the volume dielectric fluctuations [23,27]. In terms of the reflection amplitude, this means that
$$R_{\mu \nu} = R_{\zeta, \mu \nu} + R_{\varepsilon, \mu \nu} \: ,$$
where $R_{\zeta, \mu \nu }$ and $R_{\varepsilon, \mu \nu }$ are the surface and volume reflection amplitudes, respectively. For a weakly scattering layer, such that the conditions of small surface amplitude ($\sqrt {\varepsilon _j} k_0 \sigma \ll 1$) and small thickness ($\sqrt {\varepsilon _j} k_0 L \ll 1$) are satisfied, the reflection amplitudes have analytical expressions. They can be written as [23]
$$\begin{aligned} R_{\zeta, \mu \nu}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = s(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \, \rho_{\zeta, \mu \nu}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \: , \end{aligned}$$
$$\begin{aligned} R_{\varepsilon, \mu \nu} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = v(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \: \rho_{\varepsilon, \mu \nu}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \: . \end{aligned}$$

Each reflection amplitude is the product of a random contribution from the surface or the volume and of a deterministic polarization coupling factor. The contributions from surface and volume disorders take the following forms

$$\begin{aligned} \;\;\;\;\;\;\;\;\; s(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = \frac{i k_0^2}{2 \alpha_2(\mathop{\mathbf{p}_{}^{}})} \, (\varepsilon_2 -\varepsilon_1) \widetilde{\zeta}(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \;,\end{aligned}$$
$$\begin{aligned}v(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = \frac{i k_0^2}{2\alpha_2(\mathop{\mathbf{p}_{}^{}})} \Delta \widetilde{\varepsilon} (\mathop{\mathbf{p}_{}^{}} - \mathop{\mathbf{p}_{0}^{}}) \, L \: , \end{aligned}$$
where $\widetilde {\zeta }$ and $\Delta \widetilde {\varepsilon }$ are the Fourier transforms of the surface profile function and of the dielectric fluctuation. The polarization coupling factors are given by
$$\begin{aligned}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; \rho_{\zeta, \mu \nu} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = t_{12}^{(\mu)}(\mathop{\mathbf{p}_{}^{}}) \, \mathop{\mathbf{\hat{e}}_{2,\mu}^+}(\mathop{\mathbf{p}_{}^{}}) \cdot \Big[ \mathop{\mathbf{\hat{e}}_{1,\nu}^{-}}(\mathop{\mathbf{p}_{0}^{}}) + r_{21}^{(\nu)}(\mathop{\mathbf{p}_{0}^{}}) \mathop{\mathbf{\hat{e}}_{1,\nu}^+}(\mathop{\mathbf{p}_{0}^{}}) \Big] \: , \end{aligned}$$
$$\begin{aligned}\rho_{\varepsilon, \mu \nu} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = t_{12}^{(\mu)}(\mathop{\mathbf{p}_{}^{}}) \, \mathop{\mathbf{\hat{e}}_{2,\mu}^+}(\mathop{\mathbf{p}_{}^{}}) \cdot \mathop{\mathbf{\hat{e}}_{2,\nu}^{-}}(\mathop{\mathbf{p}_{0}^{}}) \, t_{21}^{(\nu)}(\mathop{\mathbf{p}_{0}^{}}) \: , \end{aligned}$$
where $r_{ji}^{(\nu )}$ and $t_{ji}^{(\nu )}$ are the Fresnel reflection and transmission amplitudes for a $\nu$-polarized plane wave incident on a planar surface from medium $i$ to medium $j$ [see for example Ref. [23], Eq. (A4)]. The polarization coupling factors depend only on the geometry of the reference system, namely, a planar interface between two homogeneous media with dielectric functions $\varepsilon _1$ and $\varepsilon _2$. Physically, they describe the polarization response of an electric dipole source radiating in the reference medium [23,28].

The polarization coupling factors have interesting properties, that will be useful in the following. First, it can be verified that the surface and volume factors are different only for $\mu = \nu = p$. Second, it is also interesting to note that for normal incidence ($\mathop {\mathbf {p}_{0}^{}} = \mathop {\mathbf {0}_{}^{}}$) the two polarization coupling factors in Eq. (11) are equal for all $\mathop {\mathbf {p}_{}^{}}$ and any pair of polarization states ($\mu$, $\nu$). Finally, they are real-valued functions in the radiative region ($|\mathbf {p}| < \sqrt {\varepsilon _1} k_0$). In summary, the polarization coupling factors satisfy

$$\begin{aligned} \rho_{\zeta, \mu s} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = \rho_{\varepsilon, \mu s}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \equiv \rho_{\mu s} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \: , \end{aligned}$$
$$\begin{aligned} \rho_{\zeta, s p}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = \rho_{\varepsilon, s p} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \equiv \rho_{s p} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \: , \end{aligned}$$
$$\begin{aligned} \rho_{\zeta, \mu \nu}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{0}_{}^{}}) = \rho_{\varepsilon, \mu \nu} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{0}_{}^{}})\: , \end{aligned}$$
$$\begin{aligned} \rho_{\zeta, \mu \nu}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \in \mathbb{R} \ \ \mathrm{and} \ \ \rho_{\varepsilon, \mu \nu}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \in \mathbb{R} \ \ \mathrm{for} \ \ |\mathop{\mathbf{p}_{}^{}}|, |\mathop{\mathbf{p}_{0}^{}}| \leq \sqrt{\varepsilon}_1 k_0 \: . \end{aligned}$$

We now turn to the expression for the degree of polarization of the backscattered light. For an incident plane wave in state $(\mathop {\mathbf {p}_{0}^{}},\nu )$, the degree of polarization of a wave scattered in direction $\mathop {\mathbf {p}_{}^{}}$ is defined as [29]

$$\mathcal{P}^{(\nu)}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = \left(1 - 4 \frac{\det \mathop{\mathbf{J}_{}^{(\nu)}}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}})}{\big[ \mathrm{Tr} \mathop{\mathbf{J}_{}^{(\nu)}}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \big]^2} \right)^{1/2} \: ,$$
where $\mathop {\mathbf {J}_{}^{(\nu )}}$ is the Jones coherence matrix with matrix elements
$$J^{(\nu)}_{\mu \mu^\prime} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = \left\langle R_{\mu \nu}^{} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) R_{\mu^\prime \nu}^{*} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \right\rangle \: .$$

We see that the degree of polarization is directly obtained from the reflection amplitude $R_{\mu \nu }$. It characterizes the statistical correlation between different vector components of the scattered field, given a state of polarization of the incident field (different incident states can lead to different degrees of polarization). Explicit expressions for the determinant and the trace of the coherence matrix can be obtained by inserting Eq. (9) into Eq. (14). For uncorrelated surface and volume disorders, we would simply have

$$\begin{aligned} \det \mathop{\mathbf{J}_{\mathrm{unco}}^{(\nu)}} = \left\langle |s|^2 \right\rangle \left\langle |v|^2 \right\rangle \left| \rho_{\zeta, p \nu} \rho_{\varepsilon, s \nu} - \rho_{\varepsilon, p \nu} \rho_{\zeta, s \nu} \right|^2 \:, \end{aligned}$$
$$\begin{aligned} \mathrm{Tr} \mathop{\mathbf{J}_{\mathrm{unco}}^{(\nu)}} = \left\langle |s|^2 \right\rangle \left( |\rho_{\zeta, p \nu}|^2 + |\rho_{\zeta, s \nu}|^2 \right) + \left\langle |v|^2 \right\rangle \left( |\rho_{\varepsilon, p \nu}|^2 + |\rho_{\varepsilon, s \nu}|^2 \right) \: . \end{aligned}$$

In the presence of surface-volume correlations, additional contributions have to be taken into account, and we find that

$$\begin{aligned} \det \mathop{\mathbf{J}_{}^{(\nu)}} = &\det \mathop{\mathbf{J}_{\mathrm{unco}}^{(\nu)}} + 2 \mathrm{Re} \left( \left\langle s v^* \right\rangle \rho_{\zeta, s\nu} \rho_{\varepsilon, s\nu}^* \right) \Big[ \left\langle |s|^2\right\rangle |\rho_{\zeta, p \nu} |^2 + \left\langle |v|^2\right\rangle |\rho_{\varepsilon, p \nu} |^2 \Big]\\ &+ 2 \mathrm{Re} \left( \left\langle s v^* \right\rangle \rho_{\zeta, p\nu} \rho_{\varepsilon, p \nu}^* \right) \Big[ \left\langle |s|^2\right\rangle |\rho_{\zeta, s \nu} |^2 + \left\langle |v|^2\right\rangle |\rho_{\varepsilon, s \nu} |^2 \Big]\\ &+ 4 \mathrm{Re} \left( \left\langle s v^* \right\rangle \rho_{\zeta, p\nu} \rho_{\varepsilon, p \nu}^* \right) \mathrm{Re} \left( \left\langle s v^* \right\rangle \rho_{\zeta, s\nu} \rho_{\varepsilon, s\nu}^* \right)\\ &- \Big[ \left\langle |s|^2\right\rangle \rho_{\zeta, p \nu} \rho_{\zeta, s \nu}^* + \left\langle |v|^2\right\rangle \rho_{\varepsilon, p \nu} \rho_{\varepsilon, s \nu}^* \Big] \Big[ \left\langle s v^* \right\rangle \rho_{\zeta, s \nu} \rho_{\varepsilon, p \nu}^* + \left\langle v s^*\right\rangle \rho_{\varepsilon, s \nu} \rho_{\zeta, p \nu}^* \Big]\\ &- \Big[ \left\langle |s|^2\right\rangle \rho_{\zeta, s \nu} \rho_{\zeta, p \nu}^* + \left\langle |v|^2\right\rangle \rho_{\varepsilon, s \nu} \rho_{\varepsilon, p \nu}^* \Big] \Big[ \left\langle s v^* \right\rangle \rho_{\zeta, p \nu} \rho_{\varepsilon, s \nu}^* + \left\langle v s^* \right\rangle \rho_{\varepsilon, p \nu} \rho_{\zeta, s \nu}^* \Big]\\ &- \Big[ \left\langle sv^* \right\rangle \rho_{\zeta, p \nu} \rho_{\varepsilon, s \nu}^* + \left\langle v s^*\right\rangle \rho_{\varepsilon, p \nu} \rho_{\zeta, s \nu}^* \Big] \Big[ \left\langle s v^* \right\rangle \rho_{\zeta, s \nu} \rho_{\varepsilon, p \nu}^* + \left\langle v s^* \right\rangle \rho_{\varepsilon, s \nu} \rho_{\zeta, p \nu}^* \Big] \: , \end{aligned}$$
$$\begin{aligned} \mathrm{Tr} \mathop{\mathbf{J}_{}^{(\nu)}} = \mathrm{Tr} \mathop{\mathbf{J}_{\mathrm{unco}}^{(\nu)}} + 2 \mathrm{Re} \Big[ \left\langle sv^* \right\rangle \big( \rho_{\zeta, p \nu} \rho_{\varepsilon, p \nu}^* + \rho_{\zeta, s \nu} \rho_{\varepsilon, s \nu}^* \big) \Big] \: .\end{aligned}$$

Note that the quantities $\left \langle |s|^2\right \rangle$, $\left \langle |v|^2\right \rangle$, and $\left \langle s v^*\right \rangle$ in Eqs. (15) and (16) are proportional to the surface power spectrum $\widetilde {W}_\zeta$, the power spectrum of the dielectric fluctuation $\widetilde {W}_{\varepsilon }$, and the cross-power spectrum $\widetilde {W}_{\zeta \varepsilon }$, respectively (full expressions are given in Appendix D of Ref. [23]).

4. Conditions for depolarization

Equations (13) and (16) provide a general expression of the degree of polarization for a weakly disordered layer in the single scattering regime. This expression allows us to analyze the conditions for depolarization of the scattered light, given an incident polarized plane wave. From the properties (12) of the polarization coupling factors, we easily find that for normal incidence ($\mathop {\mathbf {p}_{0}^{}}=\mathop {\mathbf {0}_{}^{}}$), and independently of the incident polarization, one has $\mathop {\mathcal {P}_{}^{(\nu )}}(\mathop {\mathbf {p}_{}^{}},\mathop {\mathbf {0}_{}^{}}) = 1$, meaning that the scattered waves remain perfectly polarized. This can be seen by noticing that the polarization coupling factors in Eq. (16) are equal in this case, thus canceling $\mathrm {det} \mathop {\mathbf {J}_{}^{(\nu )}}$. We also find that $\mathop {\mathcal {P}_{}^{(s)}}(\mathop {\mathbf {p}_{}^{}},\mathop {\mathbf {p}_{0}^{}}) = 1$, meaning that no depolarization occurs for an $s$-polarized incident wave. Indeed, for an incident $s$-polarized wave, the field scattered by the surface and the field scattered by the volume are produced in the same polarization state, for all realizations of the scattering medium. Thus, independently of the scattering amplitudes $s(\mathop {\mathbf {p}_{}^{}},\mathop {\mathbf {p}_{0}^{}})$ and $v(\mathop {\mathbf {p}_{}^{}},\mathop {\mathbf {p}_{0}^{}})$, the total scattered field is always perfectly polarized. These two results lead to the conclusion that depolarization in the single scattering regime can only occur for a $p$-polarized incident wave at oblique incidence.

For such a wave, it is also interesting to note that substantial depolarization in the single scattering regime can only be observed for two scattering processes (surface and volume) with similar strengths. Indeed, if one of the processes dominates over the other, then the degree of polarization tends to unity. Consider, for example, the extreme case $\Delta \varepsilon = 0$ and $\zeta \neq 0$ where surface scattering dominates. In this case the scattering amplitude $v(\mathop {\mathbf {p}_{}^{}},\mathop {\mathbf {p}_{0}^{}})$ vanishes, so that $\mathrm {det} \mathop {\mathbf {J}_{}^{(p)}} = 0$ and $\mathop {\mathcal {P}_{}^{(p)}}= 1$, as can be seen from Eqs. (13) and (16). The same analysis holds for $\zeta = 0$ and $\Delta \varepsilon \neq 0$ where volume scattering dominates. This analysis is consistent with the well-known fact that for a either surface or volume scattering, there is no depolarization in the single scattering regime. Conversely, when surface and volume scattering occur simultaneously, the scattered field is the sum of a field scattered by the surface and a field scattered by the volume weighted by random factors (the scattering amplitudes $s$ and $v$). The resulting polarization state is stochastic, and the degree of polarization can decrease. Qualitatively, to observe substantial depolarization, we can deduce that the two scattering processes (surface and volume) must have different polarization responses and comparable strengths.

5. Connecting the degree of polarization to surface-volume correlations

Having these considerations in mind, we will focus on the case of a $p$-polarized wave at oblique incidence interacting with a layer with surface and volume disorders with equal strengths, meaning that $|\varepsilon _2 - \varepsilon _1| \sigma _\zeta = \sigma _\varepsilon L$, and equal correlation lengths $\ell _\zeta = \ell _\varepsilon$. Under theses conditions, the correlation functions of surface and volume disorders are identical, $W_\zeta = W_\varepsilon$, and the scattering amplitudes have equal average intensities $\langle |s|^2 \rangle = \langle |v|^2 \rangle$. To analyze the influence of surface-volume correlations on depolarization, it is useful to recast the degree of polarization in the form (see Appendix A)

$$\mathop{\mathcal{P}_{}^{(p)}}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = \left[ 1 - \frac{4 \rho_{sp}^2 (\rho_{\zeta,pp} - \rho_{\varepsilon,pp})^2 \big(1 - |\gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}})|^2 \big)}{\Big[ \rho_{\zeta, pp}^2 + \rho_{\varepsilon,pp}^2 + 2 \rho_{sp}^2 + 2 \mathrm{Re}\big( \gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \big) \big( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \big) \Big]^2} \right]^{1/2} \! ,$$
where $\gamma (\mathop {\mathbf {p}_{}^{}})$ is the spectral correlation modulator defined in Eq. (4).

5.1 Vanishing or perfect correlation

The particular case of uncorrelated surface and volume disorders, corresponding to $\gamma (\mathop {\mathbf {p}_{}^{}})=0$, has been examined in detail in Ref. [15]. It was shown that perfect depolarization can be achieved in specific scattering directions. From Eq. (17), one immediately finds that for $\gamma (\mathop {\mathbf {p}_{}^{}})=0$ the degree of polarization vanishes when the equality

$$4 \, \rho_{sp}^2 \, \Big( \rho_{\zeta, pp} - \rho_{\varepsilon, pp} \Big)^2 = \Big[ \rho_{\zeta,pp}^2 + \rho_{\varepsilon,pp}^2 + 2 \rho_{sp}^2 \Big]^2$$
is satisfied. The remarkable fact is that observation directions such that $\rho _{sp} = \rho _{\zeta,pp}= -\rho _{\varepsilon,pp}$ exist, for which condition (18) holds [15]. An example of the angular distribution of the degree of polarization for uncorrelated disorders is presented in Fig. 3(a). We observe perfect depolarization ($\mathop {\mathcal {P}_{}^{(p)}}=0$) for two scattering directions, symmetrically positioned with respect to the plane of incidence, for which Eq. (18) is satisfied. In these directions, the field scattered by the surface and the field scattered by the volume are orthogonal and weighted by uncorrelated amplitudes $s$ and $v$ with equal average intensities, leading to perfect depolarization [15]. Experimentally, partial depolarization was observed by Germer et al. on steel samples exhibiting uncorrelated surface roughness and volume heterogeneity [21].

Another extreme situation is that of perfect surface-volume correlation, corresponding to $|\gamma | = 1$. In this case, we find that $\mathcal {P}^{(p)}=1$, independently of the behavior of the polarization coupling factors. This means that single scattering from two perfectly correlated random processes does not induce any depolarization. Indeed, when the two scattering processes are perfectly correlated, the scattering amplitudes $s$ and $v$ are connected by a simple (complex-valued) multiplicative constant. Consequently, even though the polarization states for surface and volume scattering are expected to be different, the resulting scattered field possesses a deterministic polarization state, hence a unit degree of polarization. The phenomenon of absence of depolarization for perfectly correlated disorders is in fact more general. For instance, this was demonstrated theoretically and experimentally for dielectric films of SiO$_2$ on silicon obtained by a thermally grown oxide from an initial rough surface of a silicon substrate. The degree of polarization was shown to be unity for a wide range of scattering angles, and the recovered correlation factor close to one over a broad spectral range [18].

 figure: Fig. 3.

Fig. 3. Degree of polarization $\mathcal {P}^{(p)}$ versus the observation direction $\mathop {\mathbf {\hat {p}}_{}^{}} = (\hat {p}_1,\hat {p}_2)$. (a) Uncorrelated surface and volume disorders. (b) Uniform correlation with $\gamma (\mathop {\mathbf {p}_{}^{}} ) = 1/2$. (c) Shift correlation with $\gamma (\mathop {\mathbf {p}_{}^{}}) = \exp \big ( i 5 \hat {p}_1 \big ) / 2$. For all cases the angles of incidence are $\theta _0 = 75^\circ$ and $\phi _0 = 0^\circ$. Layer parameters: $\varepsilon _1 = 1$, $\varepsilon _2 = 2.25$, $\sigma _\varepsilon = 0.36$, $L =\lambda / 20$, $\ell _{\varepsilon } = \ell _{\zeta } = \lambda / 2$, $\sigma _\zeta = 1.4 \times 10^{-2} \lambda$. The parameters are chosen such that $\left \langle |s|^2 \right \rangle = \left \langle |v|^2 \right \rangle$.

Download Full Size | PDF

5.2 Partial correlation

In the presence of partial correlation between surface and volume disorders, we expect partial or perfect depolarization of the scattered light but with different features compared with uncorrelated disorders. The direct connection between the degree of polarization and the spectral correlation modulator given by Eq. (17) allows us to study the process quantitatively. In the following we examine a few situations of particular interest.

Uniform correlation — A uniform partial correlation is characterized by $\gamma (\mathop {\mathbf {p}_{}^{}}) = \gamma _0$, with $|\gamma _0| < 1$. Figure 3(b) shows the angular distribution of the degree of polarization for $\gamma _0 = 1/2$. By comparison to Fig. 3(a), we observe that the directions of perfect depolarization are shifted compared to those of the uncorrelated case, with a shift towards larger or smaller azimuthal angles, depending on the sign of $\gamma _0$.

Shift correlation — Another particular case is the wave vector dependent correlation modulator $\gamma (\mathop {\mathbf {p}_{}^{}}) = \gamma _0 \exp (i \mathop {\mathbf {a}_{}^{}} \cdot \mathop {\mathbf {p}_{}^{}})$ with $\mathop {\mathbf {a}_{}^{}}$ a constant (spatial) vector. In real space, this form of correlation corresponds to a surface profile and dielectric fluctuations that are scaled and shifted copies of each other, such that $\sigma _\varepsilon \zeta (\mathop {\mathbf {x}_{\parallel }^{}}-\mathop {\mathbf {a}_{}^{}}) = \pm \sigma _\zeta \Delta \varepsilon (\mathop {\mathbf {x}_{\parallel }^{}})$ (for $\gamma _0=\pm 1$). We show in Fig. 3(c) the angular distribution of the degree of polarization for the shift correlation with $\mathop {\mathbf {a}_{}^{}} = 5 \lambda \mathop {\mathbf {\hat {e}}_{1}^{}}$ and $\gamma _0 = 1/2$. Since $|\gamma (\mathop {\mathbf {p}_{}^{}})| = 1/2$ as for the uniform correlation examined previously, the role played by the phase term in $\gamma (\mathop {\mathbf {p}_{}^{}})$ is directly revealed by comparison with Fig. 3(b). We observe in Fig. 3(c) partial depolarization fringes whose positions are controlled by the real part of $\gamma (\mathop {\mathbf {p}_{}^{}})$. These fringes in the degree of polarization are reminiscent of similar fringes observed in the angular distribution of the diffuse intensity [23]. In the context of pseudo-random grating layers, similar fringes were observed in the degree of circular polarization when the top surface profile of the overlay was shifted compared to that of the substrate and were used to deduce the spatial shift between the two surface profiles in Refs. [19,20].

Qualitative picture for partial depolarization — We have seen that for perfectly correlated disorders ($|\gamma |=1$), the scattered field possesses a well-defined polarization state, leading to $\mathcal {P}^{(p)}=1$. Conversely, for uncorrelated surface and volume disorders, the fields scattered by the surface and the volume are non-colinear and weighted by random uncorrelated amplitudes $s$ and $v$. The resulting field is partially polarized, as illustrated schematically in Fig. 4(a), and even fully depolarized when the conditions illustrated in Fig. 4(b) are met. Starting from a vanishing degree of polarization for uncorrelated processes, increasing the surface-volume correlation can be seen as a repolarization mechanism. Indeed, even a partial correlation links the weighting amplitudes $s$ and $v$. Consequently, even for orthogonal surface and volume polarization states, the distribution of the resulting field becomes anisotropic as illustrated in Fig. 4(c), leading to partial repolarization for $|\gamma |<1$ and even total repolarization for $|\gamma |=1$. Reversely, starting in the uncorrelated case from a direction for which the polarization states are not orthogonal, correlation of the amplitudes can correct for this lack of orthogonality to make the resulting field distribution isotropic. This is the case for the directions of perfect depolarization in Fig. 3(b) for instance.

Engineering the degree of polarization — It is also instructive to examine the possibility of shaping the degree of polarization $\mathop {\mathcal {P}_{}^{(p)}}$ by an appropriate design of the surface-volume statistical correlation. For instance, one could seek to cancel the degree of polarization over a range of observation directions, or to set it to some prescribed value. We first note that $\mathop {\mathcal {P}_{}^{(p)}}=1$ in the plane of incidence (since $\rho _{sp} = 0$) independently of the surface-volume correlation, so that shaping is meaningful only for observation directions outside the plane of incidence. For a given direction of incidence $\mathop {\mathbf {p}_{0}^{}}$, consider the problem of minimizing the degree of polarization given by Eq. (17) in an observation direction $\mathbf {p}$, with $\gamma$ as the free parameter. We note that the minimizer is necessarily real and satisfies $|\gamma |<1$. Indeed, by writing $\gamma = |\gamma | \exp (i \phi )$ and assuming a fixed modulus $|\gamma |$, minimizing Eq. (17) is equivalent to minimizing

$$\mathrm{Re}\big( \gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \big) \big( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \big) = |\gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) | \cos \phi \big( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \big) \: .$$

 figure: Fig. 4.

Fig. 4. Schematic illustration of the conditions for (a) partial depolarization or (b) perfect depolarization of the light scattered from uncorrelated surface and volume disorders. Partial correlation can be seen as a repolarization mechanism, as illustrated in (c). The red and blue arrows represent the deterministic polarization states in the $(\mathop {\mathbf {\hat {e}}_{s}^{}}, \mathop {\mathbf {\hat {e}}_{p}^{}})$ basis. Summing these polarization states with random weights $s$ and $v$ produces the total field whose distribution is represented by the colored area.

Download Full Size | PDF

Depending on the sign of the factor $\rho _{\zeta,pp} \rho _{\varepsilon,pp} + \rho _{sp}^2$, we find that the minimum is reached for $\phi =0$ or $\phi =\pi$, forcing $\gamma$ to be real. The problem of minimizing $\mathop {\mathcal {P}_{}^{(p)}}$ is thus reduced to a one-dimensional problem with variable $\gamma$ and is analyzed in detail in Appendix B. The minimum is found for an optimal correlation modulator $\gamma _\star$ given by

$$\gamma_\star (\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) ={-} \frac{2 \left( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \right)}{\rho_{\zeta, pp}^2 + \rho_{\varepsilon,pp}^2 + 2 \rho_{sp}^2} \: ,$$
and the corresponding minimum value of the degree of polarization is
$$\mathop{\mathcal{P}_{{\star}}^{(p)}}(\mathop{\mathbf{p}_{}^{}}, \mathop{\mathbf{p}_{0}^{}}) = \frac{\left| \rho_{\zeta, pp} + \rho_{\varepsilon,pp} \right|}{ \sqrt{\left( \rho_{\zeta, pp} + \rho_{\varepsilon,pp} \right)^2 + 4 \rho_{sp}^2 } } \: .$$

Note that these expressions are consistent with the existence of directions exhibiting perfect depolarization for uncorrelated disorders. Indeed, in the absence of surface-volume correlation, the directions $\mathop {\mathbf {p}_{}^{}}$ of perfect depolarization are characterized by the equation $\rho _{\zeta,pp} \rho _{\varepsilon,pp} + \rho _{sp}^2 = 0$ [15]. When this condition is satisfied, we immediately find from Eq. (20) that $\gamma _\star = 0$. Since in these directions we also have $\rho _{\zeta,pp}= -\rho _{\varepsilon,pp}$ (see section 5.1), we also find that the corresponding degree of polarization vanishes.

 figure: Fig. 5.

Fig. 5. (a) Optimal correlation modulator $\gamma _\star$ versus the observation direction $\mathop {\mathbf {\hat {p}}_{}^{}} = (\hat {p}_1,\hat {p}_2)$ and (b) the corresponding degree of polarization $\mathop {\mathcal {P}_{\star }^{(p)}}(\mathop {\mathbf {p}_{}^{}})$. The black or white dots indicate the directions of perfect depolarization for $\gamma =0$ taken from Fig. 3(a). The remaining parameters are the same as in Fig. 3.

Download Full Size | PDF

To get a more general picture in the presence of surface-volume correlations, we show the optimal correlation modulator $\gamma _\star$ versus the observation direction in Fig. 5(a), and the corresponding degree of polarization in Fig. 5(b). The two directions corresponding to perfect depolarization for uncorrelated disorders are indicated by the black or white dots in Fig. 5, both lying on the contours $\gamma _\star = 0$ and $\mathop {\mathcal {P}_{\star }^{(p)}} = 0$ as expected. We also see in Fig. 5(b) that the degree of polarization vanishes on a contour indicated by the dark red color, and defining a range of observation angles over which the scattered light is fully depolarized. To define this contour, we can set the right-hand side of Eq. (21) to zero, which leads to the condition $\rho _{\zeta, pp} + \rho _{\varepsilon,pp} = 0$. Recalling the definition of the polarization coupling factors in Eq. (11), we find that perfect depolarization is obtained for observation directions $\mathop {\mathbf {p}_{}^{}}$ satisfying

$$t_{12}^{(p)}(\mathop{\mathbf{p}_{}^{}}) \mathop{\mathbf{\hat{e}}_{2,p}^+}(\mathop{\mathbf{p}_{}^{}}) \cdot \Big[ \mathop{\mathbf{\hat{e}}_{1,p}^{-}}(\mathop{\mathbf{p}_{0}^{}}) + r_{21}^{(p)}(\mathop{\mathbf{p}_{0}^{}}) \mathop{\mathbf{\hat{e}}_{1,p}^+}(\mathop{\mathbf{p}_{0}^{}}) + \mathop{\mathbf{\hat{e}}_{2,p}^{-}}(\mathop{\mathbf{p}_{0}^{}}) \, t_{21}^{(p)}(\mathop{\mathbf{p}_{0}^{}})\Big] = 0 \: .$$

This implicit equation defines the dark red contour in Fig. 5(b).

Remarks — At this point we can make a few remarks. The discussion and equations presented previously in this section were obtained under the assumption that the surface and volume disorders had equal strengths and equal correlation lengths. It is difficult to establish general results if these assumptions are not met. Nevertheless, a general trend can be identified as follows. For a fixed direction of incidence, the expression for the degree of polarization depends on the scattering direction. For different forms of the power spectra and equal global scattering strength (i.e., integrated over all directions) one may have a region of scattering directions for which $\left \langle |s|^2\right \rangle > \left \langle |v|^2\right \rangle$, a region for which $\left \langle |s|^2\right \rangle < \left \langle |v|^2\right \rangle$, and a contour on which $\left \langle |s|^2\right \rangle = \left \langle |v|^2\right \rangle$. We may expect to observe identical depolarization on this contour as that observed in the case studied in this section, and a degree of polarization going progressively towards unity away from the contour. The typical width around this contour within which significant depolarization will be observed should depend on the rate at which the two power spectra differ from one another away from the contour $\left \langle |s|^2\right \rangle = \left \langle |v|^2\right \rangle$; the faster they differ the narrower this characteristic width. Again, since depolarization will strongly depend on the forms of the power spectra $\left \langle |s|^2\right \rangle$ and $\left \langle |v|^2\right \rangle$ and the scattering strengths, it is difficult to determine the general depolarization behavior. For the sake of clarity in communicating the main results of this work, we assumed a particular class of scattering systems for which the power spectra and the strengths are the same. We leave a discussion of the more general case for future work.

6. Conclusion

In summary, we have derived a general expression for the degree of polarization of the light scattered from a weakly scattering layer exhibiting both surface and volume scattering. This expression puts forward the direct connection between the degree of polarization and the cross-correlation function of the surface and volume disorders. We have analyzed depolarization of the backscattered light for uncorrelated, perfectly correlated and partially correlated disorders. The analysis shows that measuring the degree of polarization could be used, in principle, to assess the statistical correlation between the surface roughness and the bulk dielectric fluctuations. In addition, an appropriate shaping of the correlation function could be used to shape the degree of polarization over a range of scattering angles. Manufacturing such samples could be done by a combination of techniques such as direct laser writing and lithography for example, or via self-organization processes such as the recently demonstrated phase separation of thin barium borosilicate layers [24,25]. As an alternative, engineering the cross-correlation between the surfaces of a randomly rough film (for which the general analysis developed in this work also holds) could potentially be simpler to achieve, either by lithography techniques [19,20] or by successive and controlled exposure of speckle patterns onto photosensitive coatings.

A. Expression for the degree of polarization for a p-polarized incident wave

In this appendix we derive Eq. (17). For correlated surface and volume disorders, making use of properties (12), it is easy to show that Eq. (16a) can be rewritten as

$$\begin{aligned}\mathrm{det} \mathop{\mathbf{J}_{}^{(p)}} = &\left\langle |s|^2 \right\rangle \left\langle |v|^2 \right\rangle \left| \rho_{\zeta, p p} - \rho_{\varepsilon, p p} \right|^2 \rho_{s p}^2 + 4 \mathrm{Re} \big( \left\langle s v^* \right\rangle \big)^2 \rho_{\zeta,pp} \, \rho_{\varepsilon, pp} \, \rho_{sp}^2\\ &- \Big| \rho_{\zeta,pp} \left\langle s v^* \right\rangle + \rho_{\varepsilon, pp} \left\langle v s^* \right\rangle \Big|^2 \rho_{sp}^2 \: . \end{aligned}$$

From the expressions for the various covariances derived in Ref. [23] (see Eqs. (D4-D6) in [23]) we obtain

$$\begin{aligned}\mathrm{det} \mathop{\mathbf{J}_{}^{(p)}} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = &\left[ \frac{k_0^4}{4 |\alpha_2(\mathop{\mathbf{p}_{}^{}})|^2} \right]^2 (\varepsilon_2 -\varepsilon_1)^2 \sigma_\zeta^2 \sigma_\varepsilon^2 L^2 \widetilde{W}_\zeta (\mathop{\mathbf{p}_{}^{}} - \mathop{\mathbf{p}_{0}^{}}) \widetilde{W}_\varepsilon (\mathop{\mathbf{p}_{}^{}} - \mathop{\mathbf{p}_{0}^{}}) \, \rho_{s p}^2(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}})\\ &\times \Bigg[ \left| \rho_{\zeta, p p}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) - \rho_{\varepsilon, p p}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \right|^2 + 4 \mathrm{Re} ( \gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) )^2 \rho_{\zeta,pp}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \, \rho_{\varepsilon, pp}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}})\\ &- \Big| \rho_{\zeta,pp}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) + \rho_{\varepsilon, pp}(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \gamma^*(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \Big|^2 \Bigg]\\ = &\left[ \frac{k_0^4}{4 |\alpha_2(\mathop{\mathbf{p}_{}^{}})|^2} \right]^2 (\varepsilon_2 -\varepsilon_1)^2 \sigma_\zeta^2 \sigma_\varepsilon^2 L^2 \widetilde{W}_\zeta (\mathop{\mathbf{p}_{}^{}} - \mathop{\mathbf{p}_{0}^{}}) \widetilde{W}_\varepsilon (\mathop{\mathbf{p}_{}^{}} - \mathop{\mathbf{p}_{0}^{}}) \, \rho_{s p}^2(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}})\\ &\times \Bigg[ \left| \rho_{\zeta, p p}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) - \rho_{\varepsilon, p p}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \right|^2\\ &+ \rho_{\zeta, p p}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \rho_{\varepsilon, p p}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \Big( 4 \mathrm{Re} ( \gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) )^2 - \gamma^2(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) -\gamma^{*2}(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \Big)\\ &- |\gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}})|^2 \Big( \rho_{\zeta,pp}^2(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) + \rho_{\varepsilon, pp}^2(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \Big) \Bigg] \: . \end{aligned}$$

Making use of the identity $4 (\mathrm {Re} \, z)^2 - z^2 -z^{*2} = 2 |z|^2$ valid for any complex number $z$, we finally obtain

$$\begin{aligned} \mathrm{det} \mathop{\mathbf{J}_{}^{(p)}} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = \: &\left[ \frac{k_0^4}{4 |\alpha_2(\mathop{\mathbf{p}_{}^{}})|^2} \right]^2 (\varepsilon_2 -\varepsilon_1)^2 \sigma_\zeta^2 \sigma_\varepsilon^2 L^2 \widetilde{W}_\zeta (\mathop{\mathbf{p}_{}^{}} - \mathop{\mathbf{p}_{0}^{}}) \widetilde{W}_\varepsilon (\mathop{\mathbf{p}_{}^{}} - \mathop{\mathbf{p}_{0}^{}}) \, \rho_{s p}^2(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}})\\ &\times \left| \rho_{\zeta, p p}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) - \rho_{\varepsilon, p p}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \right|^2 \Big( 1 - |\gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}})|^2 \Big)\\ = \: &\mathrm{det} \mathop{\mathbf{J}_{\mathrm{uncor}}^{(p)}} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) \, \Big( 1 - |\gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}})|^2 \Big) \: . \end{aligned}$$

To complete the derivation of the degree of polarization, we need to compute the trace of the coherence matrix, and we obtain

$$\begin{aligned}\mathrm{Tr} \mathop{\mathbf{J}_{}^{(p)}} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = &\frac{k_0^4}{4 |\alpha_2(\mathop{\mathbf{p}_{}^{}})|^2} \Bigg[ (\varepsilon_2{\,-\,}\varepsilon_1)^2 \sigma_\zeta^2 \widetilde{W}_\zeta (\mathop{\mathbf{p}_{}^{}}{\,-\,}\mathop{\mathbf{p}_{0}^{}}) \big( \rho_{\zeta,pp}^2{\,+\,}\rho_{sp}^2 \big){\,+\,}\sigma_\varepsilon^2 L^2 \widetilde{W}_\varepsilon (\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \big( \rho_{\varepsilon,pp}^2{\,+\,}\rho_{sp}^2 \big)\\ &+ 2 \mathrm{Re} \big( \gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \big) (\varepsilon_2 -\varepsilon_1) \sigma_\zeta \widetilde{W}_\zeta^{1/2} (\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \sigma_\varepsilon L \widetilde{W}_\varepsilon^{1/2} (\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \Big( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \Big) \Bigg] \\ = & \mathrm{Tr} \mathop{\mathbf{J}_{\mathrm{unco}}^{(p)}} (\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}})\\ &+ 2 \mathrm{Re} \big( \gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \big) (\varepsilon_2 -\varepsilon_1) \sigma_\zeta \widetilde{W}_\zeta^{1/2} (\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \sigma_\varepsilon L \widetilde{W}_\varepsilon^{1/2} (\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \Big( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \Big) \: . \end{aligned}$$

In the conditions $|\varepsilon _2 - \varepsilon _1| \sigma _\zeta = \sigma _\varepsilon L$ and $\ell _\zeta = \ell _\varepsilon$, that are assumed in the main text, we obtain after some simplifications

$$\mathop{\mathcal{P}_{}^{(p)}}(\mathop{\mathbf{p}_{}^{}},\mathop{\mathbf{p}_{0}^{}}) = \left[ 1 - \frac{4 \rho_{sp}^2 (\rho_{\zeta,pp} - \rho_{\varepsilon,pp})^2 \big(1 - |\gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}})|^2 \big)}{\Big[ \rho_{\zeta, pp}^2 + \rho_{\varepsilon,pp}^2 + 2 \rho_{sp}^2 + 2 \mathrm{Re}\big[ \gamma(\mathop{\mathbf{p}_{}^{}}-\mathop{\mathbf{p}_{0}^{}}) \big] \big( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \big) \Big]^2} \right]^{1/2} \: ,$$
which is Eq. (17) in the main text.

B. Minimization of the degree of polarization

In this appendix we derive Eqs. (20) and (21). We have seen in section 5 that the degree of polarization may be minimized with $\gamma$ as a free parameter, and that the minimizer is real valued. This means that we can search the point $|\gamma _\star |$ such that $\partial \mathcal {P}^{(p)} / \partial |\gamma | = 0$, or equivalently, $\partial \left (\mathcal {P}^{(p)} \right )^2/ \partial |\gamma | = 0$. Using the notations $A = 4 \rho _{sp}^2 (\rho _{\zeta,pp} - \rho _{\varepsilon,pp})^2$, $B = \rho _{\zeta, pp}^2 + \rho _{\varepsilon,pp}^2 + 2 \rho _{sp}^2$, and $C = 2 \big | \rho _{\zeta,pp} \rho _{\varepsilon,pp} + \rho _{sp}^2 \big |$, we can write

$$\left[ \mathop{\mathcal{P}_{}^{(p)}} \right]^2 = 1 - \frac{A\big(1 - |\gamma|^2 \big)}{\Big[ B - C |\gamma| \Big]^2} \: ,$$
from which we find that
$$\frac{\partial \left[ \mathop{\mathcal{P}_{}^{(p)}} \right]^2}{\partial |\gamma| } = 2A \frac{|\gamma| (B - C |\gamma|) - C(1 - |\gamma|^2) }{\Big[ B - C |\gamma| \Big]^3} \: .$$

The minimizer $|\gamma _\star |$ is the solution to the equation

$$|\gamma^\star| (B - C |\gamma_\star|) - C(1 - |\gamma_\star|^2) = 0 \: ,$$
which immediately leads to
$$|\gamma_\star| = \frac{C}{B} = \frac{2 \big| \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \big|}{\rho_{\zeta, pp}^2 + \rho_{\varepsilon,pp}^2 + 2 \rho_{sp}^2} \: .$$

Since $\mathrm {sign}(\gamma _\star ) = - \mathrm {sign} ( \rho _{\zeta,pp} \rho _{\varepsilon,pp} + \rho _{sp}^2 )$, we end up with

$$\gamma_\star{=}{-} \frac{2 \left( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \right)}{\rho_{\zeta, pp}^2 + \rho_{\varepsilon,pp}^2 + 2 \rho_{sp}^2} \: .$$

By inserting Eq. (32) into Eq. (28) we also find that

$$\begin{aligned}\left[ \mathop{\mathcal{P}_{{\star}}^{(p)}} \right]^2 &= 1 - \frac{A\big(1 - \frac{C^2}{B^2} \big)}{\Big[ B - \frac{C^2}{B} \Big]^2} = 1 - \frac{A}{B^2 - C^2} \: ,\\ &= 1 - \frac{4 \rho_{sp}^2 (\rho_{\zeta,pp} - \rho_{\varepsilon,pp})^2}{ \left( \rho_{\zeta, pp}^2 + \rho_{\varepsilon,pp}^2 + 2 \rho_{sp}^2 \right)^2 - 4 \left( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \right)^2 } \: . \end{aligned}$$

The denominator of the second term on the right-hand side can be recast as

$$\begin{aligned}\left( \rho_{\zeta, pp}^2 + \rho_{\varepsilon,pp}^2 + 2 \rho_{sp}^2 \right)^2 - 4 \left( \rho_{\zeta,pp} \rho_{\varepsilon,pp} + \rho_{sp}^2 \right)^2 =\\ \left( \rho_{\zeta, pp} - \rho_{\varepsilon,pp}\right)^2 \left[ \left( \rho_{\zeta, pp} + \rho_{\varepsilon,pp} \right)^2 + 4 \rho_{sp}^2 \right] \: , \end{aligned}$$
which finally leads to
$$\left[ \mathop{\mathcal{P}_{{\star}}^{(p)}} \right]^2 = \frac{\left( \rho_{\zeta, pp} + \rho_{\varepsilon,pp} \right)^2}{ \left( \rho_{\zeta, pp} + \rho_{\varepsilon,pp} \right)^2 + 4 \rho_{sp}^2 } \: .$$

This completes the derivation of Eqs. (20) and (21) in the main text.

Funding

Agence Nationale de la Recherche (ANR-10-IDEX-0001-02 PSL*, ANR-10-LABX-24, ANR-15-CHIN-0003).

Acknowledgments

We are grateful to Eugenio Méndez for fruitful discussions.

Disclosures

The authors declare no conflicts of interest.

Data Availability

Numerical data underlying the results presented in this paper were produced using a homemade computer program. They are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. T. Oates, H. Wormeester, and H. Arwin, “Characterization of plasmonic effects in thin films and metamaterials using spectroscopic ellipsometry,” Prog. Surf. Sci. 86(11-12), 328–376 (2011). [CrossRef]  

2. T. Brakstad, M. Kildemo, Z. Ghadyani, and I. Simonsen, “Dispersion of polarization coupling, localized and collective plasmon modes in a metallic photonic crystal mapped by Mueller Matrix Ellipsometry,” Opt. Express 23(17), 22800–22815 (2015). [CrossRef]  

3. M. Wang, A. Löhle, B. Gompf, M. Dressel, and A. Berrier, “Physical interpretation of Mueller matrix spectra: a versatile method applied to gold gratings,” Opt. Express 25(6), 6983–6996 (2017). [CrossRef]  

4. C. M. Lam and A. Ishimaru, “Calculation of Mueller matrices and polarization signatures for a slab of random medium using vector radiative transfer,” IEEE Trans. Antennas Propagat. 41(7), 851–862 (1993). [CrossRef]  

5. A. H. Hielscher, A. A. Eick, J. R. Mourant, D. Shen, J. P. Freyer, and I. J. Bigio, “Diffuse backscattering Mueller matrices of highly scattering media,” Opt. Express 1(13), 441–453 (1997). [CrossRef]  

6. K. Drozdowicz-Tomsia, F. Xie, N. Calander, I. Gryczynski, K. Gryczynski, and E. M. Goldys, “Depolarized light scattering from colloidal gold nanoparticles,” Chem. Phys. Lett. 468(1-3), 69–74 (2009). [CrossRef]  

7. J. Ellis, P. Caillard, and A. Dogariu, “Off-diagonal Mueller matrix elements in backscattering from highly diffusive media,” J. Opt. Soc. Am. A 19(1), 43–48 (2002). [CrossRef]  

8. P. A. Letnes, A. A. Maradudin, T. Nordam, and I. Simonsen, “Calculation of the Mueller matrix for scattering of light from two-dimensional rough surfaces,” Phys. Rev. A 86(3), 031803 (2012). [CrossRef]  

9. I. Simonsen, “Optics of surface disordered systems: a random walk through rough surface scattering phenomena,” Eur. Phys. J. Spec. Top. 181(1), 1–103 (2010). [CrossRef]  

10. C. M. Lam and A. Ishimaru, “Mueller matrix calculation for a slab of random medium with both random rough surfaces and discrete particles,” IEEE Trans. Antennas Propagat. 42(2), 145–156 (1994). [CrossRef]  

11. T. A. Germer, “Angular dependence and polarization of out-of-plane optical scattering from particulate contamination, subsurface defects, and surface microroughness,” Appl. Opt. 36(33), 8798–8805 (1997). [CrossRef]  

12. J. Sorrentini, M. Zerrad, and C. Amra, “Statistical signatures of random media and their correlation to polarization properties,” Opt. Lett. 34(16), 2429–2431 (2009). [CrossRef]  

13. J. Dupont, X. Orlik, A. Ghabbach, M. Zerrad, G. Soriano, and C. Amra, “Polarization analysis of speckle field below its transverse correlation width : application to surface and bulk scattering,” Opt. Express 22(20), 24133–24141 (2014). [CrossRef]  

14. A. Ghabbach, M. Zerrad, G. Soriano, S. Liukaityte, and C. Amra, “Depolarization and enpolarization DOP histograms measured for surface and bulk speckle patterns,” Opt. Express 22(18), 21427–21440 (2014). [CrossRef]  

15. J.-P. Banon, I. Simonsen, and R. Carminati, “Perfect depolarization in single scattering of light from uncorrelated surface and volume disorder,” Opt. Lett. 45(23), 6354–6357 (2020). [CrossRef]  

16. A. D. Buckingham and M. J. Stephen, “A theory of the depolarization of light scattered by a dense medium,” Trans. Faraday Soc. 53, 884–893 (1957). [CrossRef]  

17. R. Ossikovski and O. Arteaga, “Statistical meaning of the differential Mueller matrix of depolarizing homogeneous media,” Opt. Lett. 39(15), 4470–4473 (2014). [CrossRef]  

18. T. A. Germer, “Measurement of roughness of two interfaces of a dielectric film by scattering ellipsometry,” Phys. Rev. Lett. 85(2), 349–352 (2000). [CrossRef]  

19. T. A. Germer, “Measurement of lithographic overlay by light-scattering ellipsometry,” Proc. SPIE 4780, 72–79 (2002). [CrossRef]  

20. T. A. Germer and M. J. Fasolka, “Characterizing surface roughness of thin films by polarized light scattering,” Proc. SPIE 5188, 264–275 (2003). [CrossRef]  

21. T. A. Germer, T. Rinder, and H. Rothe, “Polarized light-scattering measurements of polished and etched steel surfaces,” Proc. SPIE 4100, 148–155 (2000). [CrossRef]  

22. T. A. Germer and M. E. Nadal, “Modeling the appearance of special effect pigment coatings,” Proc. SPIE 4447, 77–86 (2001). [CrossRef]  

23. J.-P. Banon, I. Simonsen, and R. Carminati, “Single scattering of polarized light by correlated surface and volume disorder,” Phys. Rev. A 101(5), 053847 (2020). [CrossRef]  

24. B. Bouteille, J. T. Fonné, E. Burov, E. Gouillart, H. Henry, H. Montigaud, P. Jop, and D. Vandembroucq, “Slow coarsening of ultra-confined phase-separated glass thin films,” Appl. Phys. Lett. 120(5), 051602 (2022). [CrossRef]  

25. B. Bouteille, “Séparation de phase dans les couches minces de verre pour la nanostructuration de surface,” PhD Theses, Sorbonne Université (2020).

26. L. Mandel and E. Wolf, Optical coherence and quantum optics (Cambridge University, 1995).

27. J. M. Elson, “Theory of light scattering from a rough surface with an inhomogeneous dielectric permittivity,” Phys. Rev. B 30(10), 5460–5480 (1984). [CrossRef]  

28. J.-P. Banon, Ø. S. Hetland, and I. Simonsen, “Physics of polarized light scattering from weakly rough dielectric surfaces: Yoneda and Brewster scattering phenomena,” Phys. Rev. A 99(2), 023834 (2019). [CrossRef]  

29. M. Born and E. Wolf, Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light (Cambridge University, 1999), 7th ed.

Data Availability

Numerical data underlying the results presented in this paper were produced using a homemade computer program. They are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. (a) Cross-section of the scattering layer in the plane $(x_1,x_3)$, showing both a rough surface and volume dielectric fluctuations. (b) Schematics defining the incident and scattering wave vectors. The incident wave vector lies in the plane $(x_1,x_3)$, with an in-plane component $\mathop {\mathbf {p}_{0}^{}}$ and a direction defined by the polar angle of incidence $\theta _0$. The scattered wave vector has an arbitrary in-plane component $\mathop {\mathbf {p}_{}^{}}$, and a direction defined by the scattering angles $(\theta,\phi )$.
Fig. 2.
Fig. 2. Definition of the wave vectors for the incident and scattered fields, and of the unit vectors defining the $s$ and $p$ polarization components.
Fig. 3.
Fig. 3. Degree of polarization $\mathcal {P}^{(p)}$ versus the observation direction $\mathop {\mathbf {\hat {p}}_{}^{}} = (\hat {p}_1,\hat {p}_2)$. (a) Uncorrelated surface and volume disorders. (b) Uniform correlation with $\gamma (\mathop {\mathbf {p}_{}^{}} ) = 1/2$. (c) Shift correlation with $\gamma (\mathop {\mathbf {p}_{}^{}}) = \exp \big ( i 5 \hat {p}_1 \big ) / 2$. For all cases the angles of incidence are $\theta _0 = 75^\circ$ and $\phi _0 = 0^\circ$. Layer parameters: $\varepsilon _1 = 1$, $\varepsilon _2 = 2.25$, $\sigma _\varepsilon = 0.36$, $L =\lambda / 20$, $\ell _{\varepsilon } = \ell _{\zeta } = \lambda / 2$, $\sigma _\zeta = 1.4 \times 10^{-2} \lambda$. The parameters are chosen such that $\left \langle |s|^2 \right \rangle = \left \langle |v|^2 \right \rangle$.
Fig. 4.
Fig. 4. Schematic illustration of the conditions for (a) partial depolarization or (b) perfect depolarization of the light scattered from uncorrelated surface and volume disorders. Partial correlation can be seen as a repolarization mechanism, as illustrated in (c). The red and blue arrows represent the deterministic polarization states in the $(\mathop {\mathbf {\hat {e}}_{s}^{}}, \mathop {\mathbf {\hat {e}}_{p}^{}})$ basis. Summing these polarization states with random weights $s$ and $v$ produces the total field whose distribution is represented by the colored area.
Fig. 5.
Fig. 5. (a) Optimal correlation modulator $\gamma _\star$ versus the observation direction $\mathop {\mathbf {\hat {p}}_{}^{}} = (\hat {p}_1,\hat {p}_2)$ and (b) the corresponding degree of polarization $\mathop {\mathcal {P}_{\star }^{(p)}}(\mathop {\mathbf {p}_{}^{}})$. The black or white dots indicate the directions of perfect depolarization for $\gamma =0$ taken from Fig. 3(a). The remaining parameters are the same as in Fig. 3.

Equations (48)

Equations on this page are rendered with MathJax. Learn more.

ε ( x ) = ε 1 + H ( ζ ( x ) x 3 ) ( ε 2 ε 1 + Δ ε ( x ) H ( x 3 + L ) ) ,
ζ ( x ) ζ ( x ) = σ ζ 2 W ζ ( x x ) ,
Δ ε ( x ) Δ ε ( x ) = σ ε 2 W ε ( x x ) ,
ζ ( x ) Δ ε ( x ) = σ ζ σ ε W ζ ε ( x x ) ,
W ζ ( x x ) = exp ( | x x | 2 ζ 2 ) ,
W ε ( x x ) = exp ( | x x | 2 ε 2 ) ,
W ~ ζ ε ( p ) = γ ( p ) W ~ ζ 1 / 2 ( p ) W ~ ε 1 / 2 ( p ) ,
E 0 ( x ) = ν = p , s E 0 , ν e ^ 1 , ν ( p 0 ) exp [ i p 0 x i α 1 ( p 0 ) x 3 ] ,
α j ( p ) = ( ε j k 0 2 p 2 ) 1 / 2 , R e ( α j ) 0 , I m ( α j ) 0 ,
e ^ j , s ± ( p ) = e ^ 3 × p ^ ,
e ^ j , p ± ( p ) = ± α j ( p ) p ^ | p | e ^ 3 ε j k 0 .
E ~ s ( p , x 3 ) = μ = p , s e ^ 1 , μ + ( p ) ν = p , s R μ ν ( p , p 0 ) E 0 , ν exp ( i α 1 ( p ) x 3 ) ,
R μ ν = R ζ , μ ν + R ε , μ ν ,
R ζ , μ ν ( p , p 0 ) = s ( p , p 0 ) ρ ζ , μ ν ( p , p 0 ) ,
R ε , μ ν ( p , p 0 ) = v ( p , p 0 ) ρ ε , μ ν ( p , p 0 ) .
s ( p , p 0 ) = i k 0 2 2 α 2 ( p ) ( ε 2 ε 1 ) ζ ~ ( p p 0 ) ,
v ( p , p 0 ) = i k 0 2 2 α 2 ( p ) Δ ε ~ ( p p 0 ) L ,
ρ ζ , μ ν ( p , p 0 ) = t 12 ( μ ) ( p ) e ^ 2 , μ + ( p ) [ e ^ 1 , ν ( p 0 ) + r 21 ( ν ) ( p 0 ) e ^ 1 , ν + ( p 0 ) ] ,
ρ ε , μ ν ( p , p 0 ) = t 12 ( μ ) ( p ) e ^ 2 , μ + ( p ) e ^ 2 , ν ( p 0 ) t 21 ( ν ) ( p 0 ) ,
ρ ζ , μ s ( p , p 0 ) = ρ ε , μ s ( p , p 0 ) ρ μ s ( p , p 0 ) ,
ρ ζ , s p ( p , p 0 ) = ρ ε , s p ( p , p 0 ) ρ s p ( p , p 0 ) ,
ρ ζ , μ ν ( p , 0 ) = ρ ε , μ ν ( p , 0 ) ,
ρ ζ , μ ν ( p , p 0 ) R     a n d     ρ ε , μ ν ( p , p 0 ) R     f o r     | p | , | p 0 | ε 1 k 0 .
P ( ν ) ( p , p 0 ) = ( 1 4 det J ( ν ) ( p , p 0 ) [ T r J ( ν ) ( p , p 0 ) ] 2 ) 1 / 2 ,
J μ μ ( ν ) ( p , p 0 ) = R μ ν ( p , p 0 ) R μ ν ( p , p 0 ) .
det J u n c o ( ν ) = | s | 2 | v | 2 | ρ ζ , p ν ρ ε , s ν ρ ε , p ν ρ ζ , s ν | 2 ,
T r J u n c o ( ν ) = | s | 2 ( | ρ ζ , p ν | 2 + | ρ ζ , s ν | 2 ) + | v | 2 ( | ρ ε , p ν | 2 + | ρ ε , s ν | 2 ) .
det J ( ν ) = det J u n c o ( ν ) + 2 R e ( s v ρ ζ , s ν ρ ε , s ν ) [ | s | 2 | ρ ζ , p ν | 2 + | v | 2 | ρ ε , p ν | 2 ] + 2 R e ( s v ρ ζ , p ν ρ ε , p ν ) [ | s | 2 | ρ ζ , s ν | 2 + | v | 2 | ρ ε , s ν | 2 ] + 4 R e ( s v ρ ζ , p ν ρ ε , p ν ) R e ( s v ρ ζ , s ν ρ ε , s ν ) [ | s | 2 ρ ζ , p ν ρ ζ , s ν + | v | 2 ρ ε , p ν ρ ε , s ν ] [ s v ρ ζ , s ν ρ ε , p ν + v s ρ ε , s ν ρ ζ , p ν ] [ | s | 2 ρ ζ , s ν ρ ζ , p ν + | v | 2 ρ ε , s ν ρ ε , p ν ] [ s v ρ ζ , p ν ρ ε , s ν + v s ρ ε , p ν ρ ζ , s ν ] [ s v ρ ζ , p ν ρ ε , s ν + v s ρ ε , p ν ρ ζ , s ν ] [ s v ρ ζ , s ν ρ ε , p ν + v s ρ ε , s ν ρ ζ , p ν ] ,
T r J ( ν ) = T r J u n c o ( ν ) + 2 R e [ s v ( ρ ζ , p ν ρ ε , p ν + ρ ζ , s ν ρ ε , s ν ) ] .
P ( p ) ( p , p 0 ) = [ 1 4 ρ s p 2 ( ρ ζ , p p ρ ε , p p ) 2 ( 1 | γ ( p p 0 ) | 2 ) [ ρ ζ , p p 2 + ρ ε , p p 2 + 2 ρ s p 2 + 2 R e ( γ ( p p 0 ) ) ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) ] 2 ] 1 / 2 ,
4 ρ s p 2 ( ρ ζ , p p ρ ε , p p ) 2 = [ ρ ζ , p p 2 + ρ ε , p p 2 + 2 ρ s p 2 ] 2
R e ( γ ( p p 0 ) ) ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) = | γ ( p p 0 ) | cos ϕ ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) .
γ ( p p 0 ) = 2 ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) ρ ζ , p p 2 + ρ ε , p p 2 + 2 ρ s p 2 ,
P ( p ) ( p , p 0 ) = | ρ ζ , p p + ρ ε , p p | ( ρ ζ , p p + ρ ε , p p ) 2 + 4 ρ s p 2 .
t 12 ( p ) ( p ) e ^ 2 , p + ( p ) [ e ^ 1 , p ( p 0 ) + r 21 ( p ) ( p 0 ) e ^ 1 , p + ( p 0 ) + e ^ 2 , p ( p 0 ) t 21 ( p ) ( p 0 ) ] = 0 .
d e t J ( p ) = | s | 2 | v | 2 | ρ ζ , p p ρ ε , p p | 2 ρ s p 2 + 4 R e ( s v ) 2 ρ ζ , p p ρ ε , p p ρ s p 2 | ρ ζ , p p s v + ρ ε , p p v s | 2 ρ s p 2 .
d e t J ( p ) ( p , p 0 ) = [ k 0 4 4 | α 2 ( p ) | 2 ] 2 ( ε 2 ε 1 ) 2 σ ζ 2 σ ε 2 L 2 W ~ ζ ( p p 0 ) W ~ ε ( p p 0 ) ρ s p 2 ( p , p 0 ) × [ | ρ ζ , p p ( p , p 0 ) ρ ε , p p ( p , p 0 ) | 2 + 4 R e ( γ ( p p 0 ) ) 2 ρ ζ , p p ( p , p 0 ) ρ ε , p p ( p , p 0 ) | ρ ζ , p p ( p , p 0 ) γ ( p p 0 ) + ρ ε , p p ( p p 0 ) γ ( p p 0 ) | 2 ] = [ k 0 4 4 | α 2 ( p ) | 2 ] 2 ( ε 2 ε 1 ) 2 σ ζ 2 σ ε 2 L 2 W ~ ζ ( p p 0 ) W ~ ε ( p p 0 ) ρ s p 2 ( p , p 0 ) × [ | ρ ζ , p p ( p , p 0 ) ρ ε , p p ( p , p 0 ) | 2 + ρ ζ , p p ( p , p 0 ) ρ ε , p p ( p , p 0 ) ( 4 R e ( γ ( p p 0 ) ) 2 γ 2 ( p p 0 ) γ 2 ( p p 0 ) ) | γ ( p p 0 ) | 2 ( ρ ζ , p p 2 ( p , p 0 ) + ρ ε , p p 2 ( p , p 0 ) ) ] .
d e t J ( p ) ( p , p 0 ) = [ k 0 4 4 | α 2 ( p ) | 2 ] 2 ( ε 2 ε 1 ) 2 σ ζ 2 σ ε 2 L 2 W ~ ζ ( p p 0 ) W ~ ε ( p p 0 ) ρ s p 2 ( p , p 0 ) × | ρ ζ , p p ( p , p 0 ) ρ ε , p p ( p , p 0 ) | 2 ( 1 | γ ( p p 0 ) | 2 ) = d e t J u n c o r ( p ) ( p , p 0 ) ( 1 | γ ( p p 0 ) | 2 ) .
T r J ( p ) ( p , p 0 ) = k 0 4 4 | α 2 ( p ) | 2 [ ( ε 2 ε 1 ) 2 σ ζ 2 W ~ ζ ( p p 0 ) ( ρ ζ , p p 2 + ρ s p 2 ) + σ ε 2 L 2 W ~ ε ( p p 0 ) ( ρ ε , p p 2 + ρ s p 2 ) + 2 R e ( γ ( p p 0 ) ) ( ε 2 ε 1 ) σ ζ W ~ ζ 1 / 2 ( p p 0 ) σ ε L W ~ ε 1 / 2 ( p p 0 ) ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) ] = T r J u n c o ( p ) ( p , p 0 ) + 2 R e ( γ ( p p 0 ) ) ( ε 2 ε 1 ) σ ζ W ~ ζ 1 / 2 ( p p 0 ) σ ε L W ~ ε 1 / 2 ( p p 0 ) ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) .
P ( p ) ( p , p 0 ) = [ 1 4 ρ s p 2 ( ρ ζ , p p ρ ε , p p ) 2 ( 1 | γ ( p p 0 ) | 2 ) [ ρ ζ , p p 2 + ρ ε , p p 2 + 2 ρ s p 2 + 2 R e [ γ ( p p 0 ) ] ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) ] 2 ] 1 / 2 ,
[ P ( p ) ] 2 = 1 A ( 1 | γ | 2 ) [ B C | γ | ] 2 ,
[ P ( p ) ] 2 | γ | = 2 A | γ | ( B C | γ | ) C ( 1 | γ | 2 ) [ B C | γ | ] 3 .
| γ | ( B C | γ | ) C ( 1 | γ | 2 ) = 0 ,
| γ | = C B = 2 | ρ ζ , p p ρ ε , p p + ρ s p 2 | ρ ζ , p p 2 + ρ ε , p p 2 + 2 ρ s p 2 .
γ = 2 ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) ρ ζ , p p 2 + ρ ε , p p 2 + 2 ρ s p 2 .
[ P ( p ) ] 2 = 1 A ( 1 C 2 B 2 ) [ B C 2 B ] 2 = 1 A B 2 C 2 , = 1 4 ρ s p 2 ( ρ ζ , p p ρ ε , p p ) 2 ( ρ ζ , p p 2 + ρ ε , p p 2 + 2 ρ s p 2 ) 2 4 ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) 2 .
( ρ ζ , p p 2 + ρ ε , p p 2 + 2 ρ s p 2 ) 2 4 ( ρ ζ , p p ρ ε , p p + ρ s p 2 ) 2 = ( ρ ζ , p p ρ ε , p p ) 2 [ ( ρ ζ , p p + ρ ε , p p ) 2 + 4 ρ s p 2 ] ,
[ P ( p ) ] 2 = ( ρ ζ , p p + ρ ε , p p ) 2 ( ρ ζ , p p + ρ ε , p p ) 2 + 4 ρ s p 2 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.