Numerical Calculations. All numerical simulations were done using COMSOL Multiphysics version 5.6. The model is composed of a rectangular volume with a square footprint with sides Λ = 250 nm and periodic boundary conditions were employed for both u and v directions. The DWP unit cell is divided into two parts of air and glass, with one gold nanobrick placed against the glass region. The corners of the nanobrick are rounded with a 5 nm radius. The refractive index of air is set as 1 and that of glass as 1.46 for all wavelengths, while the gold permittivity was interpolated as a function of wavelength from experimental tabulated values28.
Using this model, the complex reflection and transmission coefficients for the glass/OMS/air interface are calculated for both propagation directions (i.e., normal incident from glass or air) and for light linearly polarized along both u and v separately. These are used to calculate the total reflection coefficient \({r}_{\text{F}\text{P}}\) by including the gold substrate with the FP equation19−21
$${r}_{\text{F}\text{P}}={r}_{12}+\frac{{t}_{12}{t}_{21}{r}_{23}{e}^{i2k{n}_{2}{T}_{a}}}{1-{r}_{21}{r}_{23}{e}^{i2k{n}_{2}{T}_{a}}}$$
1
Here, rmn (tmn) denotes the reflection (transmission) coefficients for light incident on material n from material m and the materials are numbered 1, 2, 3 for respectively glass substrate, air and gold substrate, and k is the wavenumber in vacuum. Note that the effect of the OMS on the interface between the air and glass substrate is included in r12, r21, t12 and t21, and that even for normal incidence rFP are polarization-dependent due to the anisotropy of the OMS layer. The reflection coefficient r23 from the air/gold interface is calculated directly using Fresnel equation. An illustration of the DWP geometry is shown in Supplementary Fig. S2 together with plots explaining the behavior of the total reflection coefficient rFP with varying Ta.
For air gaps much smaller than the wavelengths, there is near-field coupling between the nanobricks and the gold substrate in addition to the FP resonances1,20, requiring numerical simulations including also the gold substrate in COMSOL. The results obtained using FP equation were confirmed to give the same results as COMSOL simulations with the whole glass/nanobrick/air/gold substrate model when Ta>100 nm for a wavelength of λ = 800 nm.
As a final comment, to compare with the measurements, the imaginary part of the gold permittivity is increased by three times in the simulations of Figs. 2, 4 and Supplementary Figs. S7, S9-S13, accounting for the surface roughness and grain boundary effects of the fabricated gold nanobricks as well as the increased damping associated with the titanium (Ti) adhesion layers between gold/glass interface.
Fabrication. The OMS for developing MEMS-OMS DWP were fabricated using electron beam lithography (EBL), thin-film deposition, and lift-off techniques. First, a 100-nm-thick poly(methyl methacrylate) (PMMA A2, MicroChem) layer and a 40-nm thick conductive polymer layer (AR-PC 5090, Allresist) were successively spin-coated on the glass substrate (Borofloat 33 wafer, Wafer Universe). Note that the glass substrate was preprocessed to have a 10-µm-high circular pedestal using optical lithography and wet etching, and the OMS pattern was defined on the pedestal using EBL (JEOL JSM-6500F field-emission SEM with a Raith Elphy Quantum lithography system). After development, the OMS layer was formulated by depositing a 1-nm Ti adhesion layer and a 50-nm gold layer (Tornado 400, Cryofox) followed by lift-off in acetone (Supplementary Fig. S4). The pedestal on the glass substrate is very effective for reducing the possible contaminants between the MEMS mirror and OMS surface, thus improving the stability and repeatability of the DWP components. The MEMS mirror is fabricated using standard semiconductor manufacturing processes (Supplementary Fig. S3), in which thin-film lead zirconate titanate (PZT) is incorporated for long-stroke, low-voltage electrical actuation. For use in the MEMS-OMS component, the ultra-flat MEMS mirror was sputtered with a 100 nm gold layer. After the gold deposition, the MEMS mirror surface is inspected with a white light interferometry (Zygo NewView 6300), showing overall good flatness and roughness all over the whole MEMS mirror (i.e., ~ 3 mm diameter) (Supplementary Fig. S3).
The MEM-OMS-based DWP component (Supplementary Fig. S4) was assembled by gluing the MEMS mirror with the glass substrate upon which OMS is structured, and then glued to a printed circuit board (PCB), followed by gold wire bonding process between the MEMS mirror and PCB for enabling simple electrical connection to a voltage controller used to actuate the MEMS mirror.
Characterization. The experimental setup is shown in Supplementary Fig. S8. A collimated fiber-coupled supercontinuum laser (SuperK Extreme, NKT) was directed through a HWP (AHWP10M-980, Thorlabs), a mirror, a linear polarizer (Pol1; LPNIR050-MP2, Thorlabs), two beam splitters (BS1,2; CCM1-BS014, Thorlabs) successively, and then focused onto the DWP samples by an objective (Obj; M Plan Apo, ×20/0.42NA, Mitutoyo). The combination of HWP and Pol1 is used for altering the input LP states as well as the intensity. The reflected light was collected by the same objective and passing through two beam splitters (BS2,3; CCM1-BS014, Thorlabs) and a tube lens (TL; TTL200-S8, Thorlabs), generating the first direct image plane where an iris is placed for filtering out the reflection light within the DWP area. The first direct image is then transformed by a relay lens (RL; AC254-200-B-ML, f = 200 mm, Thorlabs) to the corresponding Fourier image and captured by a CCD camera (CCD; DCC1545M, Thorlabs), according to a 2f configuration. Note that a flip lens (FL; AC254-100-B-ML, f = 100 mm, Thorlabs) is used for switching between the direct and Fourier images, and Stokes analyzer composed of a QWP (AQWP10M-980, Thorlabs) and a linear polarizer (Pol2; LPNIR050-MP2, Thorlabs) is implemented before the CCD camera for performing full Stokes polarimetry23. Two beam splitters are configured for cross-compensating the polarization-dependent phase-shifts in the beam splitters for both incidence and reflection routes.
To obtain the wavelength-resolved full Stokes parameters, we replaced the CCD camera with a fiber-coupled spectrometer (QE Pro, Ocean Optics) and conducted measurements at the Fourier image plane. By rotating the QWP and Pol2, we recorded polarization-resolved spectra of Ix(λ), Iy(λ), Ia(λ), Ib(λ), Ir(λ), Il(λ), and the stokes parameters (s1, s2, s3) are calculated as s1 = (Ix(λ)–Iy(λ))/(Ix(λ)+Iy(λ)), s2 = (Ia(λ)–Ib(λ))/(Ix(λ)+Iy(λ)), s3 = (Ir(λ)–Il(λ))/(Ix(λ)+Iy(λ)). For a reasonable comparison with the simulations, the stokes parameters (s1, s2, s3) are normalized to the polarized proportion of the reflected light beam: \({S}_{\text{1,2},3}=\frac{{s}_{\text{1,2},3}}{\text{D}\text{O}\text{P}}\), and the degree of polarization (DOP) is defined as \(\text{D}\text{O}\text{P}=\sqrt{{s}_{1}^{2}+{s}_{2}^{2}+{s}_{3}^{2}}\).
The coordinates system used is indicated in the lower-left inset of Fig. 1a, with z being the optical axis, x and y are transverse axes in the laboratory frame of reference, while u and v are transverse axes oriented along the long and short sides of the rectangular nanobricks. The angle between u and x is denoted θDWP, while the angle between the x axis and polarization direction of LP incident light is denoted θLP.
To estimate the switching speeds between different orthogonal LP and CP bases (i.e., different DWP status), the setup described above is modified by replacing the input laser and CCD camera with a cw Ti:sapphire laser (Spectra-Physics 3900S, wavelength range: 700 to 1000 nm), and a photodetector (PD; PDA20CS-EC, Thorlabs), respectively. The signals from the PD are acquired with an oscilloscope (DSOX2024A, Keysight). In the measurement, the MEMS-OMS-based DWP is actuated with periodically alternating voltages and different polarizations (i.e., |x>, |y>, |r> and |l>) can be filtered by the Stokes analyzer.