Abstract
This paper presents a snapshot spectroscopic Mueller matrix polarimetry based on spectral modulation. The polarization state generator consists of a linear polarizer in front of two high-order retarders, and the polarization state analyzer is formed by two non-polarization beam splitters incorporated with three high-order retarder/linear analyzer pairs. It can simultaneously generate three modulated spectra used for reconstructing the 16 spectroscopic Mueller elements of the sample. Since each of the modulated spectra produces seven separate channels equally spaced in the Fourier domain, the channel bandwidth can be enhanced efficiently compared with the conventional spectrally modulated spectroscopic Mueller matrix polarimetry. The feasibility of the proposed spectroscopic Mueller matrix polarimetry is demonstrated by the experimental measurement of an achromatic quarter-wave plate.
© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
Spectroscopic Mueller matrix polarimetry (SMMP) is a well-established technique that can measure all 16 elements of the 4×4 Mueller matrix of a material sample as a function of wavelength, which plays an important role in many scientific and technological areas, such as material characterization, biomedical optics, and nano-structure metrology [1–3]. A typical SMMP system consists of a light source and a polarization state generator (PSG) that generates spectroscopic polarization states to interact with the sample, a polarization state analyzer (PSA) equipped with a spectrometer that analyses the spectroscopic polarization states after the interaction. Up to now, several types of SMMP have been developed [4–12]. Among them, spectrally modulated SMMP (SM-SMMP) is a more attractive approach due to its capacity to acquire the spectra of the 16 Mueller elements of the sample in a single integration period of the detector array without mechanically moving components or electric polarization modulation devices (snapshot mode). As such, SM-SMMP is considered as an appropriate solution for in situ monitoring of complex fast processes.
In conventional SM-SMMP, the PSG and PSA are both composed of two high-order retarders and a linear polarizer, which can produce a modulated spectrum that contains separate channels within the Fourier domain. By extracting the desired channels, adequate equations to reconstruct the 16 spectroscopic Mueller elements can be obtained [10,11]. Narrow channel bandwidth is the key factor that restricts the applications of conventional SM-SMMP. As the modulated spectrum comprises more than 25 channels, the channel bandwidth of the conventional SM-SMMP is less than 1/25 of the total bandwidth. Such a narrow bandwidth is more susceptible to cross-talk between adjacent channels. If cross-talk occurs, the spectroscopic Mueller elements cannot be accurately measured. Thus, conventional SM-SMMP is limited to measure the sample whose Mueller elements are varying slowly with wavelength (information contained in the channel are bandlimited to within the bandwidth of a single channel), it may produce unsatisfactory results when measuring the spectroscopic Mueller elements with sharp spectral lines. Recently, several pioneering works have achieved an increase in channel bandwidth of spectrally modulated Stokes polarimetry [13–18]. More recently, our group proposed a hybrid spectro-temporal modulated SMMP, which can efficiently enhance the channel bandwidth [19]. However, the mechanically movable component in the system makes it unable to operate in snapshot mode. Therefore, much work remains to be done in order to accomplish a snapshot SM-SMMP with large channel bandwidth.
In this paper, we propose a modified snapshot spectrally modulated SMMP (MSM-SMMP), which increases the channel bandwidth by a factor of more than 25/7 compared to the conventional SM-SMMP. The optical configuration and operating principle of MSM-SMMP are described in Section 2. Section 3 is the experimental demonstration, and the conclusion is included in Section 4.
2. Theoretical analysis
2.1 Instrumentation
The proposed MSM-SMMP configuration scheme is illustrated in Fig. 1. It is consists of five parts: 1) a collimating optics with a Tungsten-Halogen lamp as a broadband light source, 2) the polarization state generator (PSG) consists of a linear polarizer P and two high-order retarders R01 and R02, 3) the measured sample, 4) the polarization state analyzer (PSA) consists of two non-polarization beam splitters (NPS1 and NPS2) and three high-order retarders R1, R2, R3 respectively followed by a linear analyzer A1, A2, A3, 5) a triple-spectrum sensing module consists of three parabolic mirrors, three multimode optical fibers with a diameter of 1mm and three array sensor spectrometers. The high-order retarders R01, R02, R1, R2, and R3 are built of quartz crystals with the same thickness of 23mm. The polarizer P, the analyzers A1∼A3 are α-BBO Glan-Taylor prisms with an extinction ratio of less than 5×10−6. The transmission axes of the polarizer P and the analyzers A1∼A3 are oriented at ${0^ \circ }$ to x-axis. The fast axis of high-order retarders R01, R02, R1, R2, and R3 are oriented at ${45^ \circ }$, ${0^ \circ }$, ${45^ \circ }$, $- {68.4^ \circ }$, and $- {21.6^ \circ }$ to x-axis, respectively.
2.2 Measurement method
Using Mueller calculations, the triple modulated spectra ${I_n}(\sigma ),n = 1,2,3$ recorded by the spectrometers can be determined as
Eq. (1) shows that each of the modulated spectra ${I_n}(\sigma )$ includes seven different carrier frequency components 0, ${\pm} \delta (\sigma )$, ${\pm} 2\delta (\sigma )$ and ${\pm} 3\delta (\sigma )$, which respectively carries the information about $1/32{K_n}{S_0}(\sigma ){D_{0n}}$, $1/32{K_n}{S_0}(\sigma ){D_{ {\pm} 1n}}$, $1/32{K_n}{S_0}(\sigma ){D_{ {\pm} 2n}}$, and $1/32{K_n}{S_0}(\sigma ){D_{ {\pm} 3n}}$. Therefore, inverse Fourier transform of every single ${I_n}(\sigma )$ produces seven separated channels equally spaced in the optical path difference (OPD):
Like the conventional SM-SMMP, to reconstruct the 16 Mueller elements, the ${D_{zn}}$ coefficients must be obtained firstly [11]. By filtering the channels ${A_{n0}}(h)$, ${A_{n1}}(h - L)$, ${A_{n2}}(h - 2L)$, ${A_{n3}}(h - 3L)$ and performing Fourier transform, we can obtain the following equations:
Since the phase terms ${e^{i\delta }}$, ${e^{2i\delta }}$, ${e^{3i\delta }}$ and ${K_n}{S_0}(\sigma )$ in Eq. (3a)∼Eq. (3d) are independent of the measured sample, they can be easily calibrated by measuring a linear polarizer with its polarization direction aligned at 45°relative to x axis [21–24]. Then, ${D_{0n}}$, ${D_{1n}}$, ${D_{2n}}$, and ${D_{3n}}$ can be determined. To reconstruct the 16 Mueller elements ${m_{kl}}(\sigma )(k,l \in [0,3])$, a 21-dimension vector F is defined as
From the magnitudes of ${D_{1n}}$, ${D_{2n}}$, ${D_{3n}}$, and ${D_{0n}}$ in Table 1, the $21 \times 16$-dimension transformation Matrix $\textrm{W}$ can be found, and we have
Then, the vector ${\mathbf{V}_\mathbf{m}}$ can be obtained by
where ${\mathbf{W}^\textrm{ + }}$ is the pseudo-inverse of $\mathbf{W}$.In the above measurement method, the specific ${\alpha _n}$ (${45^ \circ }$, $- {68.4^ \circ }$, and $- {21.6^ \circ }$) are obtained by shuffled complex evolution for minimizing the condition number (CN) of $\mathbf{W}$ that ensures the insensitivity of the reconstructed Mueller elements to fluctuations in vector $\mathbf{F}$[13,25]. The matrix $\mathbf{W}$ has a minimum CN of 43.43, which is similar to the CN of measurement matrix of the dual-rotating retarder SMMP in Ref. [26]. In addition, each of ${C_n}(h)$ has the same total bandwidth, which is defined as the twice summation of the OPDs produced by the high-order retarders R01, R02, and Rn, and the seven separate channels in ${C_n}(h)$ are equally spaced in OPD domain. Thus, the channel bandwidth of the filtered channels for the reconstruction is 1/7 of the total bandwidth of every single ${C_n}(h)$.
3. Experimental results
To demonstrate the feasibility of the proposed MSM-SMMP, we firstly performed a measurement on an achromatic quarter wave plate (AQWP) with an azimuth of ${12^ \circ }$. Figure 2 shows the triple modulated spectra ${I_n}(\sigma )$ recorded by the spectrometers. As can be seen, ${I_n}(\sigma )$ are modulated with the wavenumber dependences of the phase retardances of the high-order retarders. The inverse Fourier transforms of ${I_n}(\sigma )$ that produce the separate channels in OPD are shown in Fig. 3. The channels in Fig. 3(a), Fig. 3(b), and Fig. 3(c) have almost the same total bandwidth of 1.24mm. By performing the reconstruction described in Section 2, the 16 spectroscopic Mueller elements of the AQWP are reconstructed, which are shown in Fig. 4. The zero theoretical values of ${m_{01}}(\sigma )$, ${m_{02}}(\sigma )$, ${m_{03}}(\sigma )$, ${m_{10}}(\sigma )$, ${m_{20}}(\sigma )$, ${m_{30}}(\sigma )$, and ${m_{33}}(\sigma )$ can be observed as well as the symmetry in ${m_{12}}(\sigma )\textrm{ = }{m_{21}}(\sigma )$, ${m_{13}}(\sigma )\textrm{ = } - {m_{31}}(\sigma )$, and ${m_{23}}(\sigma )\textrm{ = } - {m_{32}}(\sigma )$[20].
The integration time of the spectrometers is set to be 8ms in the experiment. The measurement time of 50ms is a suggestion in this paper for averaging and signal processing steps. The measurement accuracy of the proposed MSM-SMMP is estimated through comparisons with the results obtained by a commercial SMMP (VASE ELLIPSOMETER, J.A.Woollam, Inc., USA) that performed the same measurement. The experimental results and the Mueller elements measured by the commercial equipment are highly consistent, and the residual errors shown in the inset of Fig. 4 are less than 0.02 over most of the spectral range. The residual errors come mainly due to the calibration errors of the system parameters, including the azimuths of the system components, the retardances of the high-order retarders, and a slight error can be due to the signal processing performed in the Fourier domain.
To further investigate the feasibility of the proposed MSM-SMMP, the AQWP rotated at the azimuths, denoted by $\theta$, in the range of ${0^ \circ }\sim {90^ \circ }$ with an interval of ${22.5^ \circ }$ are measured. Figure 5 illustrates the reconstructed 16 Mueller matrix elements at the azimuths corresponding to the wavenumber of $16340c{m^{ - 1}}$, $13340c{m^{ - 1}}$ and $10340c{m^{ - 1}}$. As can be seen, the elements ${m_{01}}$, ${m_{02}}$, ${m_{03}}$, ${m_{10}}$, ${m_{20}}$, ${m_{30}}$, and ${m_{33}}$ are independent with the azimuth of the measured AQWP, which is consistent with the fact that these elements have zero theoretical values. For a perfect AQWP, the elements varied with azimuth have the function form of ${m_{11}} = {\cos ^2}2\theta$, ${m_{21}} = {m_{12}} = \cos 2\theta \sin 2\theta$, ${m_{31}} ={-} {m_{13}} ={-} \sin 2\theta$, ${m_{32}} ={-} {m_{23}} ={-} \cos 2\theta$, and ${m_{22}} = {\sin ^2}2\theta$[20]. Since the AQWP has a retardance accuracy of $\lambda /100$ in the spectral range of 325nm∼1100nm, the reconstructed elements ${m_{11}}$, ${m_{12}}$, ${m_{13}}$, ${m_{21}}$, ${m_{22}}$, ${m_{23}}$, ${m_{31}}$, and ${m_{32}}$ are mainly dependent on the azimuths, and have good agreement with the theoretical values, which can be also observed in Fig. 5. These results prove the feasibility of the proposed MSM-SMMP in the snapshot measurement of the complete spectroscopic Mueller matrix, and the measurement speed of 50ms can satisfy most applications for in-situ monitoring of fast processes [27–29]. An even faster measurement speed can be achieved by setting the integration time of the spectrometers to a lower value. Furthermore, it can be expected in the case that the conventional SM-SMMP with 25 channels and has a total bandwidth of 1.24mm with equally spaced channel structure [10,12], the proposed MSM-SMMP can improve the channel bandwidth by a factor of 25/7 than that of the conventional SM-SMMP.
4. Conclusion
In this work, a modified spectroscopic Mueller matrix polarimetry based on spectral modulation has been developed. Compared with conventional SM-SMMP, the proposed MSM-SMMP has larger channel bandwidth due to the cancelation of carrier frequencies in the modulated spectrum. It can acquire the complete spectroscopic Mueller matrix in snapshot mode. The reliable measurement with large channel bandwidth makes the proposed snapshot MSM-SMMP ideal for in-situ monitoring of fast processes in various industrial applications.
Funding
International Cooperation and Exchange Programme (51961145305); China Postdoctoral Science Foundation (2020M673448); National Natural Science Foundation of China (61805193); Youth Innovation Promotion Association of the Chinese Academy of Sciences (2018440).
Acknowledgments
The authors thank the anonymous reviewers for their helpful comments and constructive suggestions.
Disclosures
The authors declare no conflicts of interest.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
References
1. H. Song, M. Gu, Y. Fang, X. Ho, H. Chen, S. Jiang, and Liu, “Complex optical conductivity of two-dimensional mos2: A striking layer dependency,” J. Phys. Chem. Lett. 10(20), 6246–6252 (2019). [CrossRef]
2. H. Arwin, “Application of ellipsometry techniques to biological materials,” Thin Solid Films 519(9), 2589–2592 (2011). [CrossRef]
3. X. Chen, S. Liu, C. Zhang, H. Jiang, Z. Ma, T. Sun, and Z. Xu, “Accurate characterization of nanoimprinted resist patterns using mueller matrix ellipsometry,” Opt. Express 22(12), 15165–15177 (2014). [CrossRef]
4. N. Quan, C. Zhang, T. Mu, S. Li, and C. You, “Estimation variance of dual-rotating-retarder Mueller matrix polarimeter in the presence of Gaussian thermal noise and poisson shot noise,” J. Opt. 22(2), 025701 (2020). [CrossRef]
5. J. López-Téllez and N. Bruce, “Mueller-matrix polarimeter using analysis of the nonlinear voltage–retardance relationship for liquid-crystal variable retarders,” Appl. Opt. 53(24), 5359–5366 (2014). [CrossRef]
6. S. Alali, A. Gribble, and I. A. Vitkin, “Rapid wide-field mueller matrix polarimetry imaging based on four photoelastic modulators with no moving parts,” Opt. Lett. 41(5), 1038–1041 (2016). [CrossRef]
7. S. Zhang, H. Jiang, H. Gu, X. Chen, and S. Liu, “High-speed Mueller matrix ellipsometer with microsecond temporal resolution,” Opt. Express 28(8), 10873–10887 (2020). [CrossRef]
8. Liam J. K. Cross and Dennis K. Hore, “Dual-modulator broadband infrared Mueller matrix ellipsometry,” Appl. Opt. 51(21), 5100–5110 (2012). [CrossRef]
9. L. M. Sandvik Aas, P.G. Ellingsen, B. E. Fladmark, P.A. Letnes, and M. Kildemo, “Overdetermined broadband spectroscopic Mueller matrix polarimeter designed by genetic algorithms,” Opt. Express 21(7), 8753–8762 (2013). [CrossRef]
10. M. Dubreuil, S. Rivet, B. Le Jeune, and J. Cariou, “Snapshot mueller matrix polarimeter by wavelength polarization coding,” Opt. Express 15(21), 13660–13668 (2007). [CrossRef]
11. N. Hagen, K. Oka, and E. L. Dereniak, “Snapshot Mueller matrix spectropolarimeter,” Opt. Lett. 32(15), 2100–2102 (2007). [CrossRef]
12. P. Lemaillet, S. Rivet, and B. L. Jeune, “Optimization of a snapshot Mueller matrix polarimeter,” Opt. Lett. 33(2), 144–146 (2008). [CrossRef]
13. J. M. Craven, M.W. Kudenov, and E. L. Dereniak, “False signature reduction in channeled spectropolarimetry,” Opt. Eng. 49(5), 053602 (2010). [CrossRef]
14. Kira A. Hart, M. Kupinski, D. Wu, and R. A. Chipman, “First results from an uncooled LWIR polarimeter forcubesat deployment,” Opt. Eng. 59(7), 075103 (2020). [CrossRef]
15. C. Zhang, N. Quan, and T. Mu, “Stokes imaging spectropolarimeter based on channeled polarimetry with full-resolution spectra and aliasing reduction,” Appl. Opt. 57(21), 6128–6134 (2018). [CrossRef]
16. Q. Li, F. Lu, X. Wang, and C. Zhu, “Low crosstalk polarization-difference channeled imaging spectropolarimeter using double-Wollaston prism,” Opt. Express 27(8), 11734–11747 (2019). [CrossRef]
17. A. Altaqui and M. W. Kudenov, “Phase-shifting interferometrybased Fourier transform channeled spectropolarimeter,” Appl. Opt. 58(7), 1830–1840 (2019). [CrossRef]
18. W. Ren, C. Fu, D. Wu, Y. Xie, and G. R. Arce, “Channeled compressive imaging spectropolarimeter,” Opt. Express 27(3), 2197–2211 (2019). [CrossRef]
19. N. Quan, C. Zhang, T. Mu, and C. You, “Spectroscopic Mueller matrix polarimeter based on spectro-temporal modulation,” Opt. Express 28(25), 37758–37772 (2020). [CrossRef]
20. R. A. Chipman, “Polarimetry,” in Handbook of Optics.McGraw-Hill, Inc. (1999).
21. K. Oka and T. Kato, “Spectroscopic polarimetry with a channeled spectrum,” Opt. Lett. 24(21), 1475–1477 (1999). [CrossRef]
22. F. Snik, T. Karalidi, and C. U. Keller, “Spectral modulation for full linear polarimetry,” Appl. Opt. 48(7), 1337–1346 (2009). [CrossRef]
23. A. S. Alenin and J. S. Tyo, “Generalized channeled polarimetry,” J. Opt. Soc. Am. A 31(5), 1013–1022 (2014). [CrossRef]
24. C. F. LaCasse, R. A. Chipman, and J. S. Tyo, “Band limited data reconstruction in modulated polarimeters,” Opt. Express 19(16), 14976–14989 (2011). [CrossRef]
25. Q. Y. Duan, V. K. Gupta, and S. Sorooshian, “Shuffled complex evolution approach for effective and efficient global minimization,” J Optim Theory Appl 76(3), 501–521 (1993). [CrossRef]
26. M.H. Smith, “Optimization of a dual-rotating-retarder Mueller matrix polarimeter,” Appl. Opt. 41(13), 2488–2493 (2002). [CrossRef]
27. J. Lee, J. Koh, and R. W. Collins, “Multichannel Mueller matrix ellipsometer for real-time spectroscopy of anisotropic surfaces and films,” Opt. Lett. 25(21), 1573–1575 (2000). [CrossRef]
28. A. Mock, T. Carlson, J. Vanderslice, J. Mohrmann, J. A. Woollam, E. Schubert, and M. Schubert, “Multiple-layered effective medium approximation approach to modeling environmental effects on alumina passivated highly porous silicon nanostructured thin films measured by in-situ mueller matrix ellipsometry,” Appl. Surf. Sci 421, 663–666 (2017). [CrossRef]
29. L. L. Deibler and M. H. Smith, “Measurement of the complex refractive index of isotropic materials with Mueller matrix polarimetry,” Appl. Opt. 40(22), 3659–3667 (2001). [CrossRef]