The effect of relative humidity on the polarization Mueller matrix under the oil smoke environment

For the variation of the polarized Mueller matrix of oil smoke particles under different relative humidity levels, the polarized single scattering characteristics of oil smoke particles are studied by using the Mie scattering theory, and the multiple scattering simulation is implemented with the Monte Carlo method. Variation in relative humidity is achieved by changing in mixing of the oil smoke and the water fog particles during the same dry particle filling time. Using the 36 sets of polarized Mueller matrices method, the Mueller matrix patterns of oil smoke were calculated for four conditions of 0%, 10%, 40%, and 95% relative humidity, respectively. We can verify the simulation’s correctness from the simulation and the experimental results. Specifically, as the relative humidity increases, the size of the Mueller matrix pattern increases, and the patterns of m₂₂, m₃₃, and m₄₄ related to the depolarization characteristic change significantly. Furthermore, the scattering depolarization coefficients of the Mueller matrix polar decomposition increase with the increasing relative humidity, with the ability of depolarization being continuously enhanced. This method determines differences in relative humidity using intuitive measurements of stabilized scattering patterns, which can present a theoretical basis for the impact of environmental variation on polarization detection.

1 Introduction

When polarized light is irradiated on a strong scattering medium, multiple scattering will occur within this medium, and the polarization information of the transmission light will change, so polarization is widely used in the detection of the properties of the medium [1–3]. However, as the current detection methods are primarily focused on the study of the polarization states obtained from the detection terminal, more research on the properties of the transmission medium is needed. To get more properties of the medium, the polarization information of the scattering medium is expressed in the form of a matrix. The information during the transmission can be carried and displayed using the matrix property [4–6]. It can complement polarization detection and has a vital research value in the field of polarization detection.

For the optical system, a parametric model was constructed to describe the Mueller matrix of optical elements. Liu et al. proposed a method for the comprehensive characterization of an arbitrary composite waveplate based on spectroscopic Mueller matrix polarimetry and described its polarization properties by using four polarization parameters, namely, retardance, fast axis azimuth, rotation angle, and diattenuation angle [7]. For the complex media samples, to avoid the direct study of the samples and make verification easy and accurate, the Mueller matrix patterns are applied to analyze the polarization properties of samples in forward transmission by a third-party mapping method. Currently, many researchers have applied Mueller matrices to various fields. Wang et al. compared the polarization parameters extracted from the 3 × 3 and 4 × 4 Mueller matrix decomposition methods at directions of transmission and reflection, and the parameters derived from the two types of matrices display similar qualitative relations to changes in the microstructure of the biological tissues [8]. Khaliq et al. compared linear depolarization and linear retardance for both Mueller matrix polar decomposition and Mueller matrix transformation methods, and the study indicated that the optical variables of Mueller matrix polar decomposition and Mueller matrix transformation methods are sensitively correlated to the microstructure of the tissue samples [9]. Zhang et al. presented both experimental measurements and Monte Carlo-based simulations of transmission polarization to analyze the effects of concentration variation on the water fog particles. They demonstrated better persistence of circular polarization than linear polarization in more serious fog [10]. Liu et al. derived the depolarization index to characterize scattering media by estimating the transmittance with Monte Carlo simulation and enhanced the underwater vision under diverse water turbidity [11]. Sanz et al. measured the Mueller matrix of microstructured surfaces by using the polar decomposition method to detect the defects on the scattering flat surface [12]. Liu et al. introduced Mueller matrix polarimetry (MMP) to characterize nanoimprinted grating structures, and noticeable depolarization effects from measured data were observed [13]. For further research, they investigated the non-uniform depolarization properties of SiO₂ thin films, two-dimensional Si grating, and three-dimensional Si cylinder grating by Mueller matrix polar decomposition [14]. Han et al. introduced a rapid Mueller matrix imaging system to research the structural characteristics of starch granules by testing different incoming and outgoing polarizations [15].

In this paper, we extend the research to consider the effect of the variation in relative humidity of the oil smoke environment on the polarization Mueller matrix. We present the Mie scattering theory of spherical particles combined with the Monte Carlo algorithm for calculating the polarization multiple scattering. We also analyze the depolarization characteristics under the different relative humidity levels by using the Mueller matrix polar decomposition method. To validate our practice, we compare the simulation method with measurements of oil smoke produced by glycerol atomization. This paper is organized as follows. Section 2 focuses on the experimental procedure for oil smoke particles under variable relative humidity levels, such as experimental setup, particle size measurement of dry particles, and close humidity monitoring in testing. In Section 3, we mainly describe the effect of relative humidity on the size and refractive index of oil smoke particles based on Hanel’s empirical growth equation and present the theoretical analysis of the Mueller matrix by using Monte Carlo simulation and the Mueller matrix polar decomposition method. Section 4 describes the simulation and experimental results and shows the consistency between simulation and measurement results. Finally, in Section 5, we summarize the results and further work.

2 Mueller matrix testing methods

2.1 Experimental setup

Figure 1 shows the Mueller matrix testing setup. A light beam emitted by an 83.06-mW laser at a wavelength of 532 nm passes through an attenuator (ATT) to adjust the illumination power. A rotatable linear polarizer (LP1, LPVISC050-MP2, Thorlabs, US) and a quarter wave-plate (W1, AQWP10M-580, Thorlabs, US) composed of the polarization state generator (PSG) are placed in front of the ATT to control the incident polarization state as one of horizontal (H), vertical (V), ${45}^{\circ }$ (P), ${135}^{\circ }$ (M) linear polarization states, and right-hand (R) and left-hand (L) circular polarization states. After beam expansion by the expander (BE, BE02-05-A, Thorlabs, US), the photons forward-scattered from the particles in the environment modulation system (EMS) pass through a polarization state receiver (PSR). The PSR consists of a quarter-wave plate (W2, AQWP10M-580, Thorlabs, US) and a linear polarizer (P2, LPVISC050-MP2, Thorlabs, US). Finally, the patterns can be recorded by using a CCD camera (MER-502-79U3C, Daheng Optics, China) with a pixel size of 3.45 $\mu \mathrm{m}$ ×3.45 $\mu \mathrm{m}$ and 2448 × 2048 pixels. In the environment modulation system (EMS), the oil smoke environment was generated by glycerol ultrasonic atomization with a refractive index of 1.47 (sample data provided by Antari Company, Taiwan). Before relative humidity testing, we measured the smoke particle size by using Malvern Spraytec (STP5311, made by Malvern Instruments Limited, Britain), shown in Section 2.2. The monitoring method of the relative humidity and its effect on the optical thickness is described in Section 2.3.

FIGURE 1

FIGURE 1. Schematic of the Mueller matrix measurement setup. ATT: attenuator; LP1 and LP2: linear polarizer; W1 and W2: quarter wave-plate; BE: beam expander; EMS: environment modulation system; CCD: charge coupled device; PSG: polarization state generator; PSR: polarization state receiver.

2.2 Particle size measurement of dry particles

During testing, to ensure the consistency of the concentration of dry particles at each relative humidity level, we control the amount of dry particles in the EMS by using the particle filling time. In this paper, we choose a dry particle filling time of 5 s, and then the volume distribution of the particle size can be measured directly by Spraytec, which is shown in Figure 2A. To ensure suitability of the testing data for theory simulation analysis, we converted the cumulative volume distribution to the number distribution in our previous study [16]. Figure 2B shows the cumulative number distribution of oil smoke, and the oil smoke particle has a mean radius value of 0.77 $\mu \mathrm{m}$.

FIGURE 2

FIGURE 2. Size distribution of oil smoke particles at the smoke filling time of 5 s. (A) Cumulative volume distribution; (B) cumulative number distribution.

2.3 Relative humidity monitoring in the EMS

For monitoring of the relative humidity in the EMS, this paper uses ultrasonic atomization of water (ZS-10Z, Hangzhou Zhengdao Electrical Equipment Co., Ltd, China) for humidity adjustment and the relative humidity calibrator (RHCL-2, OMEGA Engineering Inc, China) for real-time measurement in the EMS. The water fog filling will stop once the relative humidity meets the testing requirement. During the experiments, the relative humidity levels of 0%, 10%, 40%, and 95% were chosen as the low, medium, and high typical relative humidity. We characterize the medium concentration by the optical thickness according to Beer’s law to learn the concentration variation at different levels of relative humidity. The optical thickness measurement schematic diagram is shown in Figure 3, in which the light power meter (LPM, PM100D, made by Thorlabs, US) is applied for outgoing light intensity testing. The optical thickness $\tau$ can be obtained from the incident light intensity $I_1$ and emergent light intensity $I_2$, which can be expressed as [17]

$I_2=I_1\exp \left(-\tau \right)=I_1\exp \left(-\rho k_{\mathrm{e}}l\right),(1)$

where $l$ is the thickness of the medium in the EMS (units of cm), $\rho$ is the medium concentration (units of #/cm³), and $k_{\mathrm{e}}=k_s+k_a$ is the extinction cross-section (units of cm²).

FIGURE 3

FIGURE 3. Schematic diagram of optical thickness measurement.

Figure 4 shows the stability of the emergent light intensity as a function of sampling time for one of the groups at 40% relative humidity. The result shows that the smoke environment is stable 3–10 min after the smoke filling, so we can conduct testing during this time. When the relative humidity levels are 0%, 10%, 40%, and 95%, the corresponding optical thicknesses are 2.02, 2.08, 2.45, and 3.26, respectively. Increase in the optical thickness is due to the larger relative humidity resulting in the increase of the extinction cross-section $k_{\mathrm{e}}$ in Eq. 1.

FIGURE 4

FIGURE 4. Intensity of the emergent light under the relative humidity of 40%.

3 Monte Carlo calculations under different levels of relative humidity

3.1 Effect of relative humidity on oil smoke particles

Due to the hygroscopicity of the oil smoke particles, the radius and refractive index of the oil smoke particles after water absorption can be obtained from the Hanel’s empirical growth equation [18]. When the relative humidity $Hr$ varies from 0 to 1, the particle size $R$ and refractive index (real part $M_{re}$ and imaginary part $M_{ie}$) can be calculated by the empirical growth equation and the volume weighted average method, respectively, which can be expressed as

$\frac{R}{R_0}={\left(1-Hr\right)}^{-\frac{1}{u}},(2)$

$M_{re}=M_{rw}+\left(M_{ro}-M_{rw}\right){\left(R/R_o\right)}^{-3},(3)$

$M_{ie}=\frac{M_{iw}}{M_{rw}^2}+\left[\frac{M_{io}}{\left(M_{ro}^2+2\right)}-\frac{M_{iw}}{M_{rw}^2}\right]{\left(\frac{R}{R_o}\right)}^{-3}\left(M_{re}^2+2\right),(4)$

where $M_{rw}$ and $M_{iw}$ are the real and the imaginary parts of the refractive index of pure water particles, respectively; $M_{ro}$ and $M_{io}$ denote the real and imaginary parts of the refractive index of the dry oil smoke particles, respectively; $\mu$ is the quality growth factor, which takes the value of 4 [19]; $R_0$ is the radius of the dry oil smoke particle. Figure 5 shows the variation of particle radius and the refractive index with the relative humidity. The imaginary part of the particle refractive index after moisture absorption is tiny, so it can be ignored.

FIGURE 5

FIGURE 5. Effect of relative humidity on oil smoke particle parameters. (A) Particle radius; (B) real part of the refractive index.

3.2 Monte Carlo simulation

A Monte Carlo method based on Ramella-Roman et al. [3] is used to perform the calculations for the multiple scattering of oil smoke particles. It involves the change of the transmission energy and the direction of the photons, which is finally characterized by polarization.

The Stokes vector uses four parameters $I$, $Q$, $U$, and $V$ to describe the polarized light [20], in which $I$, $Q$, $U$, and $V$ correspond to the total intensity of the light, linear polarization in the horizontal or vertical direction, linear polarization in $+{45}^{\mathrm{o}}$ or $-{45}^{\mathrm{o}}$ direction, and circular polarization, respectively. The Stokes vector of the light can be expressed as [21]

$\begin{array}{l}I=E_1^2+E_2^2\\ Q=E_1^2-E_2^2\\ U=2E_1{E}_2\cos \theta \\ V=2E_1{E}_2\sin \theta ,\end{array}(5)$

where $E_1$ and $E_2$ denote parallel and perpendicular components of plane wave amplitudes, respectively, and $\theta$ is the phase difference between $E_1$ and $E_2$.

When the polarized light passes through the medium, the Mueller matrix $\mathit{M}\left(\theta ,\varphi \right)$ of the medium connects the incident Stokes vector ${\mathit{S}}_{\text{in}}$ with the outgoing vector ${\mathit{S}}_{\text{out}}$ of the transmission light, which is

${\mathit{S}}_{\text{out}}=\mathit{M}\left(\theta ,\varphi \right)\cdot {\mathit{S}}_{\text{in}},(6)$

$\left(\begin{array}{c}{I}_{out}\\ Q_{out}\\ U_{out}\\ V_{out}\end{array}\right)=\left[\begin{array}{cccc}{m}_{11}& m_{12}& m_{13}& m_{14}\\ m_{21}& m_{22}& m_{23}& m_{24}\\ m_{31}& m_{32}& m_{33}& m_{34}\\ m_{41}& m_{42}& m_{43}& m_{44}\end{array}\right]\left(\begin{array}{c}{I}_{in}\\ Q_{in}\\ U_{in}\\ V_{in}\end{array}\right),(7)$

where $\mathit{M}\left(\theta ,\varphi \right)=\mathit{R}\left(-\varphi \right)\mathit{S}\left(\theta \right)\mathit{R}\left(-\varphi \right)$, $\mathit{R}\left(-\varphi \right)$ is the rotation matrix for the Stokes vector of light transformation between the reference plane and the scattering plane.

For the spherical, isotropic, and homogeneous particles, the Mueller matrix is

$\mathit{S}\left(\theta \right)=\left[\begin{array}{cccc}{S}_{11}& S_{12}& 0& 0\\ S_{12}& S_{11}& 0& 0\\ 0& 0& S_{33}& S_{34}\\ 0& 0& -S_{34}& S_{33}\end{array}\right],(8)$

and the rotation matrix can be written as

$\mathit{R}\left(-\varphi \right)=\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& \cos \left(2\varphi \right)& -\sin \left(2\varphi \right)& 0\\ 0& \sin \left(2\varphi \right)& \cos \left(2\varphi \right)& 0\\ 0& 0& 0& 1\end{array}\right],(9)$

so $\mathit{M}\left(\theta ,\varphi \right)$ can be expressed as [22]

$\begin{array}{c}\mathit{M}\left(\theta \right)=\left[\begin{array}{cccc}{m}_{11}& m_{12}& m_{13}& m_{14}\\ m_{21}& m_{22}& m_{23}& m_{24}\\ m_{31}& m_{32}& m_{33}& m_{34}\\ m_{41}& m_{42}& m_{43}& m_{44}\end{array}\right]=\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& \cos \left(2\varphi \right)& -\sin \left(2\varphi \right)& 0\\ 0& \sin \left(2\varphi \right)& \cos \left(2\varphi \right)& 0\\ 0& 0& 0& 1\end{array}\right]\left[\begin{array}{cccc}{S}_{11}& S_{12}& 0& 0\\ S_{12}& S_{11}& 0& 0\\ 0& 0& S_{33}& S_{34}\\ 0& 0& -S_{34}& S_{33}\end{array}\right]\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& \cos \left(2\varphi \right)& -\sin \left(2\varphi \right)& 0\\ 0& \sin \left(2\varphi \right)& \cos \left(2\varphi \right)& 0\\ 0& 0& 0& 1\end{array}\right]\\ =\left[\begin{array}{cccc}{S}_{11}& S_{12}& 0& 0\\ S_{12}\cos \left(2\varphi \right)& S_{11}\cos \left(2\varphi \right)& -S_{33}\sin \left(2\varphi \right)& -S_{34}\sin \left(2\varphi \right)\\ S_{12}\sin \left(2\varphi \right)& S_{11}\sin \left(2\varphi \right)& S_{33}\cos \left(2\varphi \right)& S_{34}\cos \left(2\varphi \right)\\ 0& 0& -S_{34}& S_{33}\end{array}\right]\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& \cos \left(2\varphi \right)& -\sin \left(2\varphi \right)& 0\\ 0& \sin \left(2\varphi \right)& \cos \left(2\varphi \right)& 0\\ 0& 0& 0& 1\end{array}\right]\\ =\left[\begin{array}{cccc}{S}_{11}& S_{12}\cos \left(2\varphi \right)& -S_{12}\sin \left(2\varphi \right)& 0\\ S_{12}\cos \left(2\varphi \right)& S_{11}{\cos }^2\left(2\varphi \right)-S_{33}{\sin }^2\left(2\varphi \right)& -S_{11}\cos \left(2\varphi \right)\sin \left(2\varphi \right)-S_{33}\sin \left(2\varphi \right)\cos \left(2\varphi \right)& -S_{34}\sin \left(2\varphi \right)\\ S_{12}\sin \left(2\varphi \right)& S_{11}\cos \left(2\varphi \right)\sin \left(2\varphi \right)+S_{33}\sin \left(2\varphi \right)\cos \left(2\varphi \right)& -S_{11}{\cos }^2\left(2\varphi \right)+S_{33}{\sin }^2\left(2\varphi \right)& S_{34}\cos \left(2\varphi \right)\\ 0& -S_{34}\sin \left(2\varphi \right)& -S_{34}\cos \left(2\varphi \right)& S_{33}\end{array}\right]\end{array}.(10)$

Photon energy loss can occur in the process of the transmission; once the photon energy weight decreases below a certain threshold or the photon flies off the boundary, the photon transmission is terminated. After the photon has been scattered $n$ times, the energy weight $W_n$ is

$W_n=W_{n-1}\cdot \frac{k_s}{\left(k_s+k_a\right)},(11)$

where $k_s$ and $k_a$ are the scattering and absorption cross-sections of the polydispersed system, respectively, which can be obtained from

$k_s=\pi {\displaystyle \underset{r_1}{\overset{r_2}{\int }}}{r}^2{Q}_{sc}n\left(r\right)dr,(12)$

$k_a=\pi {\displaystyle \underset{r_1}{\overset{r_2}{\int }}}{r}^2{Q}_{ab}n\left(r\right)dr,(13)$

in which $r_1$ and $r_2$ are the lower and upper limits of the particle radius, respectively; $Q_{sc}$ and $Q_{ab}$ are the scattering and absorption factors, respectively; and $n\left(r\right)$ is the number distribution of the particle size.

3.3 Mueller matrix analysis

The Mueller matrix can be acquired in various ways; in this paper, we use a combination of 36 groups to adjust the polarization state generator (PSG) and the polarization state receiver (PSR) parts to $\mathrm{H}$, $\mathrm{V}$, $\mathrm{P}$, $\mathrm{M}$, $\mathrm{L}$, and $\mathrm{R}$, respectively. The detailed combination of the 36 groups is shown in Table 1.

TABLE 1

TABLE 1. Combination of 36 groups of the Mueller matrix.

According to the aforementioned combinations, the Mueller matrix can be derived as [23]

$\begin{array}{l}\mathit{M}=\left[\begin{array}{cccc}{m}_{11}& m_{12}& m_{13}& m_{14}\\ m_{21}& m_{22}& m_{23}& m_{24}\\ m_{31}& m_{32}& m_{33}& m_{34}\\ m_{41}& m_{42}& m_{43}& m_{44}\end{array}\right]=\\ \frac{1}{2}\left[\begin{array}{cccc}{S}_{\text{HH}}+S_{\text{HV}}+S_{\text{VH}}+S_{\text{VV}}& S_{\text{HH}}+S_{\text{HV}}-S_{\text{VH}}-S_{\text{VV}}& S_{\text{PH}}+S_{\text{PV}}-S_{\text{MP}}-S_{\text{MM}}& S_{\text{RH}}+S_{\text{RV}}-S_{\text{LH}}-S_{\text{LV}}\\ S_{\text{HH}}+S_{\text{VH}}-S_{\text{HV}}-S_{\text{VV}}& S_{\text{HH}}+S_{\text{VV}}-S_{\text{HV}}-S_{\text{VH}}& S_{\text{PH}}+S_{\text{MV}}-S_{\text{MH}}-S_{\text{PV}}& S_{\text{RH}}+S_{\text{LV}}-S_{\text{RV}}-S_{\text{LH}}\\ S_{\text{HP}}+S_{\text{VP}}-S_{\text{HM}}-S_{\text{VM}}& S_{\text{HP}}+S_{\text{VM}}-S_{\text{VP}}-S_{\text{HM}}& S_{\text{PP}}+S_{\text{MM}}-S_{\text{PM}}-S_{\text{MP}}& S_{\text{RP}}+S_{\text{LM}}-S_{\text{LP}}-S_{\text{RM}}\\ S_{\text{HR}}+S_{\text{VR}}-S_{\text{RL}}-S_{\text{LL}}& S_{\text{HR}}+S_{\text{VL}}-S_{\text{VR}}-S_{\text{HL}}& S_{\mathrm{PR}}+S_{\text{ML}}-S_{\text{MR}}-S_{\text{PL}}& S_{\text{RR}}+S_{\text{LL}}-S_{\text{RL}}-S_{\text{LR}}\end{array}\right],\end{array}(14)$

where $S$ indicates the light intensity and the subscripts $\mathrm{H}$, $\mathrm{V}$, $\mathrm{P}$, $\mathrm{M}$, $\mathrm{L}$, and $\mathrm{R}$ are the corresponding polarization states described in Table 1.

3.4 Mueller matrix polar decomposition

With the polar decomposition method, the Mueller matrix $\mathit{M}$ can be decomposed into the product of the depolarization matrix ${\mathit{M}}_{\Delta }$, retardance matrix ${\mathit{M}}_{\mathrm{R}}$, and diattenuation matrix ${\mathit{M}}_{\mathrm{D}}$ [24], which is

$\mathit{M}={\mathit{M}}_{\Delta }{\mathit{M}}_{\mathrm{R}}{\mathit{M}}_{\mathrm{D}},(15)$

in which ${\mathit{M}}_{\mathrm{R}}$ and ${\mathit{M}}_{\mathrm{D}}$ represent the transformation process of polarizations from one type of fully polarized light to another one, as well as ${\mathit{M}}_{\Delta }$ indicates the conversion of fully polarized light to partially polarized light [25,26]. ${\mathit{M}}^{\prime }$ contains only the retarder and depolarizer, which can be written as

${\mathit{M}}^{\prime }={\mathit{M}\mathit{M}}_{\mathrm{D}}^{\left(-1\right)}={\mathit{M}}_{\Delta }{\mathit{M}}_{\mathrm{R}}=\left[\begin{array}{cc}1& {\overrightarrow{\mathbf{0}}}^{\left(\mathrm{T}\right)}\\ {\overrightarrow{\mathit{P}}}_{\Delta }& {\mathit{m}}_{\Delta }\end{array}\right]\left[\begin{array}{cc}1& {\overrightarrow{\mathbf{0}}}^{\left(\mathrm{T}\right)}\\ {\overrightarrow{\mathbf{0}}}^{\left(\mathrm{T}\right)}& {\mathit{m}}_{\mathrm{R}}\end{array}\right]=\left[\begin{array}{cc}1& {\overrightarrow{\mathbf{0}}}^{\left(\mathrm{T}\right)}\\ {\overrightarrow{\mathit{P}}}_{\Delta }& {\mathit{m}}_{\Delta }{\mathit{m}}_{\mathrm{R}}\end{array}\right]=\left[\begin{array}{cc}1& {\overrightarrow{\mathbf{0}}}^{\left(\mathrm{T}\right)}\\ {\overrightarrow{\mathit{P}}}_{\Delta }& {\mathit{m}}^{\prime }\end{array}\right],(16)$

where ${\overrightarrow{\mathit{P}}}_{\Delta }$ is the polarizance of the depolarizer and ${\mathit{m}}^{\prime }={\mathit{m}}_{\Delta }{\mathit{m}}_{\mathrm{R}}$, ${\mathit{m}}_{\Delta }$, and ${\mathit{m}}_{\mathrm{R}}$ are the $3\times 3$ submatrix of the depolarization matrix and retardance matrix.

We set $\sqrt{{\lambda }_1}$, $\sqrt{{\lambda }_2}$, and $\sqrt{{\lambda }_3}$ as the eigenvalues of $m_{\Delta }$, which can be calculated as

$\begin{array}{l}{\mathit{m}}_{\Delta }=\pm {\left[{\mathit{m}}^{\prime }{\left({\mathit{m}}^{\prime }\right)}^{\mathrm{T}}+\left(\sqrt{{\lambda }_1{\lambda }_2}+\sqrt{{\lambda }_2{\lambda }_3}+\sqrt{{\lambda }_3{\lambda }_1}\right)\mathit{I}\right]}^{-1}\\ \begin{array}{cc}& \times \left[\left(\sqrt{{\lambda }_1}+\sqrt{{\lambda }_2}+\sqrt{{\lambda }_3}\right){\mathit{m}}^{\prime }{\left({\mathit{m}}^{\prime }\right)}^{\mathrm{T}}+\sqrt{{\lambda }_1{\lambda }_2{\lambda }_3}\mathit{I}\right].\end{array}\end{array}(17)$

The total depolarizer coefficient can be written as

$\Delta =1-\frac{\left|tr\left({\mathit{m}}_{\Delta }\right)\right|}{3}=1-\frac{\left|tr\left({\mathit{M}}_{\Delta }\right)-1\right|}{3},0\le \Delta \le 1,(18)$

where ${\mathit{M}}_{\Delta }$ can be obtained from Eqs 16, 17.

4 Simulation and experimental result and analysis

4.1 Simulation results

Our goal in this work is to investigate the effect of relative humidity on the polarization Mueller matrix. Since the empty EMS without particles has little impact on the Mueller matrix of the oil smoke environment, we chose dry particles and low, medium, and high typical humidity for the study. The Mueller matrix patterns of oil smoke particles at different relative humidity are simulated by the Monte Carlo simulation method in Section 3.2 with the initial parameters of the wavelength of 532 $\text{nm}$ and photon number of 10⁷. The results are shown in Figure 6, in which Figures 6A, B, C, D show the Mueller matrix simulation with each element pattern pixels of 100 × 100 at the relative humidity of 0%, 10%, 40%, and 95% respectively. We can see that the 16 elements of the Mueller matrix exhibit symmetry, which conforms to the cosine distribution according to Eq. 10. With the increasing relative humidity, the scattering pattern size of the Mueller matrix is found to increase. This is due to different relative humidity leading to variations in particle size and refractive index, which leads to a rise in the scattering and absorption coefficients of the medium and a more significant deflection of the photons scattered in the medium. Moreover, the number of photons in the pattern that deviate from the incident center point increases after transmission, which increases pattern sizes.

FIGURE 6

FIGURE 6. Mueller matrix simulation pattern under different relative humidity. (A) 0%; (B) 10%; (C) 40%; (D) 95%.

To further quantify the effect of different relative humidity on each element of the Mueller matrix, we first normalize the pixel intensity values of each element at different relative humidity by m₁₁ to eliminate the effects of light intensity and keep the intensity value of each pixel in the range of [−1 1] [27], and then according to the appearance frequency of the intensity values, the interval segmentation statistics of intensity values are performed for each normalized matrix element m₁₁–m₄₄ to get the number of intensity values in a small interval of different central values. Finally, the normalized intensity frequency distribution of each element is established, shown in Figure 7. In Figure 7, the horizontal coordinate is the central value of the intervals and the vertical coordinate is the number of intensity values occurring in this division interval. As can be seen from the figures, the three diagonal elements m₂₂, m₃₃, and m₄₄ change more significantly as the relative humidity increases. The depolarization phenomenon is more significant as observed from the Mueller matrix polarization decomposition in Section 3.4. This phenomenon is also consistent with the conclusion that the scattering depolarization of spherical, isotropic, and homogeneous particles is more significant than the retardance and diattenuation [28].

FIGURE 7

FIGURE 7. Normalized intensity frequency distribution curves of Mueller matrix elements under different relative humidity.

The depolarization coefficient pattern of Figure 6 under different relative humidity is solved by Eq. 15. Figure 8 shows the Mueller matrix depolarization coefficient patterns' three-dimensional simulation results for 0%, 10%, 40%, and 95% relative humidity. The figures show that the minimum value of the oil smoke depolarization coefficient pattern increases gradually with increasing relative humidity, which are 0.426, 0.511, 0.565, and 0.816, respectively. This indicates that the depolarization phenomenon becomes apparent as the relative humidity of the oil smoke increases.

FIGURE 8

FIGURE 8. Three-dimensional diagram of Mueller matrix simulation depolarization under different relative humidity. (A) 0%; (B) 10%; (C) 40%; (D) 95%.

4.2 Experimental results

The PSG and PSA parts are adjusted in the measurement according to the data in Table 1, and 16 sets of Mueller matrix element patterns can be obtained by calculating Eq. 14. To correspond to the image size of the Mueller matrix elements in the simulations, the raw patterns captured by the CCD with the size of 2448 $\times$ 2048 are cropped to the 500 × 500 pixel images with the center position as the original point. Based on the simulation analysis and our previous research for spherical, isotropic, and homogeneous particles [10], in this paper, we focus on the pattern variations of the three matrix elements m₂₂, m₃₃, and m₄₄ related to depolarization. Figures 9–12 show the testing results under the relative humidity of 0%, 10%, 40%, and 95%, respectively. From the results, we can see that with the increasing relative humidity, the patterns calculated by Eq. 14 keep getting more prominent, which is consistent with the trend of simulation results. Figure 13 shows the schematic diagram of the quantitative analysis process, the edge points of the m₂₂, m₃₃, and m₄₄ patterns in the horizontal leftward and rightward and vertical upward and downward directions are used as reference points, and the absolute values of the pixel intensities in the four directions are summed for the m₂₂, m₃₃, and m₄₄ patterns under the relative humidity of 0%, 10%, 40%, and 95%. According to calculation, at the four relative humidity levels, the absolute sum of the pixel intensities of m₂₂ in four directions are 1.09, 1.12, 2.46, and 2.81, for m₃₃ are 0.86, 0.92, 0.99, and 1.34, and for m₄₄ are 0.44, 0.61, 1.11, and 1.71, respectively. It is further verified that as the relative humidity increases, the photon motion becomes more deflected with the rise of the medium scattering coefficient, leading to the more significant deviation of the pattern depolarization terms m₂₂, m₃₃, and m₄₄ from the central point of propagation.

FIGURE 9

FIGURE 9. Experimental diagram of testing results under the relative humidity of 0%.

FIGURE 10

FIGURE 10. Experimental diagram of testing results under the relative humidity of 10%.

FIGURE 11

FIGURE 11. Experimental diagram of testing results under the relative humidity of 40%.

FIGURE 12

FIGURE 12. Experimental diagram of testing results under the relative humidity of 95%.

FIGURE 13

FIGURE 13. Schematic diagram of the quantitative analysis process.

4.3 Agreement between simulation and experiment

To verify the accuracy of the simulations, we compared the simulation result of Mueller matrix elements with the experimental images, and the mean square error (MSE) was calculated [29]. Each Mueller matrix element’s simulated and measured patterns under the same relative humidity are subjected to a point-taking operation. During the calculation, the center pixel is used as the center of the circle, and the points are taken at 45° intervals in the 360° direction of the circle. To make the values consistent, the edge points of the pattern are taken in these eight directions to complete the values compatible with the simulation and experiment. Finally, the normalized intensity values of the taken points are calculated for the MSE, which can be expressed as

$MSE=\frac{1}{n}{\displaystyle \sum _1^n}{\left(R-R_m\right)}^2,(19)$

where $R$ and $R_m$ are the simulated and measured normalized intensity values, respectively, and $n=8$ is the number of samples at different directions. Table 2 shows the MSE calculation of simulation and measurement for matrix elements m₂₂, m₃₃, and m₄₄ at the relative humidity of 0%, 10%, 40%, and 95%.

TABLE 2

TABLE 2. MSE between simulation and experiment.

It is apparent that the average value of the MSE with different relative humidity levels is 0.00538. This is due to the fact that the relative humidity controlled in the experiment cannot reach the ideal state, resulting in the simulation and experimental patterns not being fully consistent. However, the simulation and the observed images have the exact change trend, which can indicate the correctness of the simulation.

5 Conclusion

In this paper, we present Mie scattering and the Monte Carlo method as the scattering model to study the effect of relative humidity variations on Mueller matrix patterns and depolarization of smoke particles. On this basis, an experimental setup is established to test Mueller matrix with controllable environmental parameters. The experimental and simulation results show that (1) the Mueller matrix pattern increases with the increasing relative humidity. (2) The frequency distribution curves of m₂₂, m₃₃, and m₄₄ patterns change most obviously with the increasing relative humidity, which can reflect the more significant depolarization characteristics of the oil smoke particles. (3) The depolarization of the oil smoke particles gradually increases as relative humidity increases, which is useful to analyze the decrease of degree of polarization (DOP) for polarized incident light in a higher-relative humidity environment. (4) Comparing the simulations with the experiments, the average value of the MSE reaches 0.00538. Therefore, the study of the Mueller matrix can reflect the characteristics of the oil smoke environment under different relative humidity and extend the application range of the polarization Mueller matrix.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Materials; further inquiries can be directed to the corresponding author.

Author contributions

CS: writing–original draft and writing–review and editing. SZ: writing–original draft and writing–review and editing. QF: writing–review and editing. JZ: writing–review and editing. JD: writing–review and editing. YL: writing–review and editing.

Funding

The authors declare financial support was received for the research, authorship, and/or publication of this article. The work was supported by the National Natural Science Foundation of China (62127813; 61890963), Jilin Scientific and Technological Development Program (20210201093GX), and Education Department Program of Jilin Province (JJKH20220737KJ; JJKH20220738KJ).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Han Y, Ding L, Wang YP, Zheng HY, Fang L. Polarized light illuminated scattering characteristics of single airborne particle. J Quant Spectrosc Ra (2021) 266:107568. doi:10.1016/j.jqsrt.2021.107568

CrossRef Full Text | Google Scholar

2. Bruce NC, López-Téllez JM, Rodríguez-Núñez O, Rodríguez-Herrera OG. Permitted experimental errors for optimized variable-retarder Mueller-matrix polarimeters. Opt Express (2018) 26:13693–704. doi:10.1364/oe.26.013693

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Wang K, Fu Q, Wang B, Shi HD, Liu J, Wang Q, et al. Simplified analysis and suppression of polarization aberration in planar symmetric optical systems. Opt Express (2023) 31, 30750–30766. doi:10.1364/OE.498455

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Alali S, Kim A, Vurgun N, Wood MFG, Vitkin IA, Ahmad M. Quantitative correlation between light depolarization and transport albedo of various porcine tissues. J Biomed Opt (2012) 17, 045004. doi:10.1117/1.JBO.17.4.045004

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Zhao YZ, Li YH, He WJ, Liu Y, Fu YJ. Polarization scattering imaging experiment based on Mueller matrix. Opt Commun (2021) 490:126892. doi:10.1016/j.optcom.2021.126892

CrossRef Full Text | Google Scholar

6. Wang HY, Hu HF, Jiang JF, Li JQ, Li XB, Zhang WH, et al. Polarization differential imaging in turbid water via Mueller matrix and illumination modulation. Opt Commun (2021) 499:127274. doi:10.1016/j.optcom.2021.127274

CrossRef Full Text | Google Scholar

7. Gu HG, Chen XG, Shi YT, Jiang H, Zhang CW, Gong P, et al. Comprehensive characterization of a general composite waveplate by spectroscopic Mueller matrix polarimetry. Opt Express (2018) 26:25408. doi:10.1364/OE.26.025408

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wang YF, Guo YH, Zeng N, Cheng DS, He HH, Ma H. Study on the validity of 3×3 Mueller matrix decomposition. J Biomed Opt (2015) 20:065003. doi:10.1117/1.jbo.20.6.065003

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Khaliq A, Ashraf S, Gul B, Ahmad I. Comparative study of 3 x 3 Mueller matrix transformation and polar decomposition. Opt Commun (2021) 485:126756. doi:10.1016/j.optcom.2021.126756

CrossRef Full Text | Google Scholar

10. Zhang S, Shen CB, Feng XQ, Lu ST, Zhan JT, Duan J, et al. Light transmission polarization patterns and depolarization from water fog concentration variation. Opt Eng (2023) 62:014104. doi:10.1117/1.OE.62.1.014104

CrossRef Full Text | Google Scholar

11. Liu T, Sun T, He HH, Liu SX, Dong Y, Wu J, et al. Comparative study of the imaging contrasts of Mueller matrix derived parameters between transmission and backscattering polarimetry. Biomed Opt Express (2018) 9:4413–28. doi:10.1364/BOE.9.004413

PubMed Abstract | CrossRef Full Text | Google Scholar

12. SanzSaiz JMJM, González F, Moreno F. Polar decomposition of the Mueller matrix: a polarimetric rule of thumb for square-profile surface structure recognition. Appl Opt (2011) 50:3781–8. doi:10.1364/AO.50.003781

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chen XG, Zhang CW, Liu SY. Depolarization effects from nanoimprinted grating structures as measured by Mueller matrix polarimetry. APPL PHYS LETT (2013) 103(15):151605. doi:10.1063/1.4824760

CrossRef Full Text | Google Scholar

14. Dong ZQ, Gu HG, Zhu JL, Shi YT, Nie L, Lyu JY, et al. Nonuniform depolarization properties of typical nanostructures and potential applications. Opt Lett (2020) 45(7):1910–3. doi:10.1364/OL.389732

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Han CY, Du CY, Chen DF. Evaluation of structural and molecular variation of starch granules during the gelatinization process by using the rapid Mueller matrix imaging polarimetry system. Opt Express (2018) 26:15851. doi:10.1364/OE.26.015851

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhang S, Zhan JT, Fu Q, Duan J, Li YC, Jiang HL. Effects of environment variation of glycerol smoke particles on the persistence of linear and circular polarization. Opt Express (2020) 28:20236–48. doi:10.1364/OE.395428

PubMed Abstract | CrossRef Full Text | Google Scholar

17. van der Laan JD, Wright JB, Kemme SA, Scrymgeour DA. Superior signal persistence of circularly polarized light in polydisperse, real-world fog environments. Appl Opt (2018) 57:5464. doi:10.1364/AO.57.005464

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Hnel G. The properties of atmospheric aerosol particles as functions of the relative humidity at thermodynamic equilibrium with the surrounding moist air. Adv Geophys (1976) 19:73–188. doi:10.1016/s0065-2687(08)60142-9

CrossRef Full Text | Google Scholar

19. Fitzgerald JW. Effect of relative humidity on the aerosol backscattering coefficient at 0694- and 106-μm wavelengths. Appl Opt (1984) 23:411. doi:10.1364/AO.23.000411

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Bartel S, Hielscher AH. Monte Carlo simulations of the diffuse backscattering mueller matrix for highly scattering media. Appl. Opt. (2000) 39:1580–1588. doi:10.1364/AO.39.001580

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Nee JB. Multichannel polarization lidar measurements of aerosols and cirrus clouds. Eur Phys J Plus (2020) 237:07020. doi:10.1051/epjconf/202023707020

CrossRef Full Text | Google Scholar

22. Milun JR, George WK, Mehru˝ beo M, Brent DC, Lihong VW, Sohi R, et al. Light backscattering polarization patterns from turbid media: theory and experiment. Appl Opt (1999) 38:3399. doi:10.1364/AO.38.003399

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Du E, He HH, Zeng N, Guo YH, Liao R, He YH, et al. Two-dimensional backscattering Mueller matrix of sphere–cylinder birefringence media. J Biomed Opt (2012) 17:126016–1260167. doi:10.1117/1.JBO.17.12.126016

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Lu SY, Chipman RA. Interpretation of Mueller matrices based on polar decomposition. J OPT SOC AM A (1996) 13:1106–13. doi:10.1364/JOSAA.13.001106

CrossRef Full Text | Google Scholar

25. He C, Chang J, Salter PS, Shen YX, Dai B, Li PC, et al. Revealing complex optical phenomena through vectorial metrics. Adv Photon (2022) 4. doi:10.1117/1.AP.4.2.026001

CrossRef Full Text | Google Scholar

26. Chen ZH, Yao Y, Zhu YH, Ma H. Removing the dichroism and retardance artifacts in a collinear backscattering Mueller matrix imaging system. Opt Express (2018) 26:28288–301. doi:10.1364/OE.26.028288

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Andreas HH, Judith RM, Irving JB. Influence of particle size and concentration on the diffuse backscattering of polarized light from tissue phantoms and biological cell suspensions. Appl Opt (1997) 36:125–35. doi:10.1364/AO.36.000125

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Sheng W, Li WP, Qi J, Liu T, He HH, Dong Y, et al. Quantitative analysis of 4 × 4 mueller matrix transformation parameters for biomedical imaging. Photonics (2019) 6:34. doi:10.3390/photonics6010034

CrossRef Full Text | Google Scholar

29. Raz A, Mendlovic D. Sequential filtering for color image acquisition. Opt Express (2014) 22:26878–83. doi:10.1364/OE.22.026878

PubMed Abstract | CrossRef Full Text | Google Scholar