Phased dipole arrays for the polarization-maintaining near-field directionality
We begin with the analysis of the essential role of the metasurface-waveguide coupler in the polarization-maintaining near-field directionality. Here we consider a two-dimensional case, where the guided waves propagate only along the \(\pm \widehat{x}\)direction [Fig. 1a-d]. At the frequency of interest, the meta-waveguide supports both p-polarized and s-polarized guided waves with a same in-plane wavevector, namely \({k}_{\text{p},\text{g}\text{w}}={k}_{\text{s},\text{g}\text{w}}={k}_{\text{g}\text{w}}\); see the dispersion of guided waves in Fig. S4. The metasurface is set to be anisotropic and have a certain phase gradient along the \(\widehat{x}\) direction. That is, under the normal incidence of p-polarized (s-polarized) propagating waves, the metasurface is designed to be equivalent to a phased \({p}_{\text{x}} \left({p}_{\text{y}}\right)\) dipole array with \(N\) elements, pointing in the \(\widehat{x}\) (\(\widehat{y}\)) direction and placed along the \(x\) axis with equal spacing \(a=\frac{2\pi }{{N}_{0}\bullet {k}_{\text{g}\text{w}}}\), where \({N}_{0}\) is an integer. Each dipolar element has a progressive phase shift \({\alpha }_{\text{p}}=\frac{2\pi }{{N}_{0}}\) (\({\alpha }_{\text{s}}=-\frac{2\pi }{{N}_{0}}\)) relative to its adjacent element under the incidence of p-polarized (s-polarized) propagating waves.
By rigorously solving Maxwell’s equations, we have \({E}_{\text{z}}\left(\text{x},\text{z}\right)=\int {E}_{\text{z}}\left({k}_{\text{x}}\right){e}^{i{k}_{x}x}{e}^{i{k}_{z}|z-{z}_{0}|}d{k}_{x}\) and \({H}_{\text{z}}\left(\text{x},\text{z}\right)=0\) for the phased \({p}_{\text{x}}\) dipole array in free space [40], where
$${E}_{\text{z}}\left({k}_{\text{x}}\right)=\frac{i}{8{\pi }^{2}{\epsilon }_{0}}{k}_{\text{x}}{p}_{\text{x}0}\bullet \sum _{j=1}^{N}{e}^{\left(j-1\right)\bullet i(\frac{2\pi }{{N}_{0}}-{k}_{\text{x}}a)}$$
1
\({k}_{z}=\sqrt{{k}_{0}^{2}-{k}_{\text{x}}^{2}}\) , \({p}_{\text{x}0}\) is the moment of each dipole, \({z}_{0}\) is the vertical position of all dipoles in the \(\widehat{z}\) direction, and \({\epsilon }_{0}\) is the permittivity in free space. By contrast, we have \({E}_{\text{z}}\left(\text{x},\text{z}\right)=0\) and \({H}_{\text{z}}\left(\text{x},\text{z}\right)=\int {H}_{\text{z}}\left({k}_{\text{x}}\right){e}^{i{k}_{x}x}{e}^{i{k}_{z}|z-{z}_{0}|}d{k}_{x}\) for the phased \({p}_{\text{y}}\) dipole array in free space, where
$${H}_{\text{z}}\left({k}_{\text{x}}\right)=\frac{i\omega }{8{\pi }^{2}}\frac{{k}_{\text{x}}{p}_{\text{y}0}}{{k}_{\text{z}}}\bullet \sum _{j=1}^{N}{e}^{(j-1)\bullet i(-\frac{2\pi }{{N}_{0}}-{k}_{\text{x}}a)}$$
2
and \({p}_{\text{y}0}\) is the moment of each dipole. Detailed derivation of \({E}_{\text{z}}\left({k}_{\text{x}},z\right)\) and \({H}_{\text{z}}\left({k}_{\text{x}},z\right)\) is given in supplementary Section S1. Since \({E}_{\text{z}}\left(\text{x},\text{z}\right)\) (\({H}_{\text{z}}\left(\text{x},\text{z}\right)\)) is a representative field component for p-polarized (s-polarized) waves, we can conclude that only p-polarized (s-polarized) guided waves would be excited by the phased \({p}_{\text{x}}\) (\({p}_{\text{y}}\)) dipole array under the illumination of p-polarized (s-polarized) propagating waves. Therefore, such a combined metasurface-waveguide coupler could well preserve the polarization state during the coupling between guided waves and propagating waves.
The spatial frequency spectra of the phased \({p}_{\text{x}}\) and \({p}_{\text{y}}\) dipole arrays are shown in Fig. 1b&d, according to equations (1–2). For illustration here and below, \({N}_{0}=5\) is used, the frequency is \(\omega /2\pi =9.5\) GHz, and the in-plane wavevector of guided waves is \({k}_{\text{g}\text{w}}=1.05{k}_{0}\), where \({k}_{0}=\omega /c\) and \(c\) is the speed of light in free space. The ultrahigh spectral asymmetry at \({k}_{\text{x}}=\pm {k}_{\text{g}\text{w}}\) is achievable if \(N\) is large enough [e.g.\(N=25\) in Fig. 1b&d]. To be specific, we have \(\frac{\left|{E}_{\text{z}}\left({k}_{\text{g}\text{w}}\right)\right|}{\left|{E}_{\text{z}}\left(-{k}_{\text{g}\text{w}}\right)\right|}\gg 1\) for the \({p}_{\text{x}}\) dipole array in Fig. 1b but \(\frac{\left|{H}_{\text{z}}\left({k}_{\text{g}\text{w}}\right)\right|}{\left|{H}_{\text{z}}\left(-{k}_{\text{g}\text{w}}\right)\right|}\ll 1\) for the \({p}_{\text{y}}\) dipole array in Fig. 1d. Such an ultrahigh spectral asymmetry in the \(k\) space implies the possible realization of not only polarization-maintaining but also polarization-dependent near-field directionality. That is, the directional excitation of guided waves with their polarization same as the incident propagating waves can occur; moreover, the excited guided waves with different polarizations would flow to opposite directions; see the schematic in Fig. 1a&c.
Design Of The Reflection-free, Anisotropic, And Gradient Metasurface
We now proceed to discuss a feasible design methodology for the desired metasurface in Fig. 2. Without loss of generality, each supercell of metasurface has \({N}_{0}\) unit cells in the \(\widehat{x}\) direction; under the normal incidence of light, the polarization of the transmitted fields beneath each unit cell is the same as that of the incident fields [Fig. 2a&d], by exploiting the unique capability that anisotropic metasurface can control the transmitted fields for orthogonal polarizations independently. This way, the transmitted electric field \({E}_{m,n}\)beneath each unit cell is simply proportional to the corresponding transmission coefficient \({t}_{m,n}\), namely \({E}_{m,n}:{E}_{m,n+1}={t}_{m,n}:{t}_{m,n+1}\), where the subscript \(m=p\) or \(s\) represents the incident p-polarized or s-polarized waves, and the subscript \(n\) indicates the \({n}^{\text{t}\text{h}}\) unit cell within the supercell. Due to the deep-subwavelength size of unit cells, the transmitted field beneath each unit cell can be reasonably treated as a secondary point source with a dipole moment of \({p}_{m,n}\) [Fig. 2b&e], whose orientation is dictated by the incident electric fields. In other words, we have \({p}_{m,n}:{p}_{m,n+1}={E}_{m,n}:{E}_{m,n+1}\).
If the metasurface is reflection-free under the normal incidence of light, the corresponding transmission coefficient of each unit cell should ideally have a magnitude of unity, namely \(\left|{t}_{m,n}\right|=1\). Under this condition, we have \(\left|{p}_{m,n}\right|\bullet {e}^{iArg\left({p}_{m,n}\right)}:\left|{p}_{m,n+1}\right|\bullet {e}^{iArg\left({p}_{m,n+1}\right)}={e}^{iArg\left({t}_{m,n}\right)}:{e}^{iArg\left({t}_{m,n+1}\right)}\). As such, \(\left|{p}_{m,n+1}\right|=\left|{p}_{m,n}\right|\) and \(Arg\left({p}_{m,n+1}\right)-Arg\left({p}_{m,n}\right)=Arg\left({t}_{m,n+1}\right)-Arg\left({t}_{m,n}\right)\) can be obtained. Moreover, if \(Arg\left({t}_{\text{p},n+1}\right)-Arg\left({t}_{\text{p},n}\right)={k}_{\text{g}\text{w}}\bullet a={\alpha }_{\text{p}}\) and \(Arg\left({t}_{\text{s},n+1}\right)-Arg\left({t}_{\text{s},n}\right)=-{k}_{\text{g}\text{w}}\bullet a={\alpha }_{\text{s}}\), the effective dipole arrays \({p}_{m,n}\) exactly correspond to the phased dipole arrays proposed in Fig. 1 with the designated progressive phase shifts. Therefore, \(\left|{t}_{m,n}\right|=1\) and \(Arg\left({t}_{m,n+1}\right)-Arg\left({t}_{m,n}\right)={\alpha }_{m}\) for each unit cell in the supercell of metasurface are actually the key conditions to enable the polarization-maintaining near-field directionality.
Due to the recent advancement in metasurfaces, these key conditions can be realized. Inspired by the isotropic ABA meta-particles [41–45], the anisotropic ABBA meta-particles can be adopted for the design of each unit cell in the metasurface, because they can enable us to achieve arbitrary transmission phases, along with a high transmission amplitude. For illustration, we show one example in Fig. 2c&f; see the fabricated metasurface in Fig. S3. The phase condition of \(Arg\left({t}_{m,n+1}\right)-Arg\left({t}_{m,n}\right)={\alpha }_{m}\) is achieved, as can be seen from the perfect match between the theoretical and numerical results in Fig. 2c&f. The magnitude condition of \(\left|{t}_{m,n}\right|=1\) is approximately realized, since we always have \(\left|{t}_{m,n}\right|>0.9\). The slight discrepancy in amplitude is mainly caused by the material loss (i.e. the material loss in the dielectric).
Near-field measurement of the polarization-maintaining near-field directionality
The fabricated metasurface can now be combined with the judiciously designed meta-waveguide. In order to achieve the high conversion efficiency between propagating waves and guided waves, the vertical distance between the metasurface and the meta-waveguide should be optimized. Here the optimized distance is 10 mm [Fig. S5]. Based on this metasurface-waveguide coupler, we carry out the microwave measurement in Figs. 3–4. In order to clearly demonstrate the polarization-maintaining near-field directionality with high coupling efficiency, both the near-field and far-field scattering measurements are implemented.
Figure 3 shows the near-field measurement. Under the normal incidence of p-polarized propagating waves [Fig. 3a-c], only p-polarized guided waves are efficiently excited. Moreover, 97.3% of the excited p-polarized guided waves propagates to the right side (namely \(+\widehat{x}\) direction). Similarly, under the normal incidence of s-polarized propagating waves [Fig. 3d-f], only s-polarized guided waves are excited. By contrast, 95.2% of the excited s-polarized guided waves propagates to the left side (\(-\widehat{x}\) direction). These measured results in Fig. 3c&f exhibit good agreements with the numerical ones in Fig. 3a-b&d-e, as carried out by the FDTD simulation.
Observation of the polarization-maintaining near-field directionality with high coupling efficiency for arbitrary incident polarization states
To quantitatively characterize the polarization-maintaining and polarization-dependent near-field directionality, Fig. 4 shows the coupling efficiency, namely the ratio between the excited power of guided waves propagating to the designated direction (\({P}_{\text{p},\text{g}\text{w},\text{r}\text{i}\text{g}\text{h}\text{t}}+{P}_{\text{s},\text{g}\text{w},\text{l}\text{e}\text{f}\text{t}}\)) and the total incident power of propagating waves (\({P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}={P}_{\text{p},\text{p}\text{w}}+{P}_{\text{s},\text{p}\text{w}}\)). Here \({P}_{\text{p},\text{p}\text{w}}\) and \({P}_{\text{s},\text{p}\text{w}}\) stand for the incident powers of p-polarized and s-polarized propagating waves, respectively; \({P}_{\text{p},\text{g}\text{w},\text{r}\text{i}\text{g}\text{h}\text{t}}\) (\({P}_{\text{s},\text{g}\text{w},\text{l}\text{e}\text{f}\text{t}}\)) corresponds to the excited power of the desired p-polarized (s-polarized) guided waves propagating to the designated right (left) side. For illustration, we show \({P}_{\text{p},\text{g}\text{w},\text{r}\text{i}\text{g}\text{h}\text{t}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}\) and \({P}_{\text{s},\text{g}\text{w},\text{l}\text{e}\text{f}\text{t}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}\) in Fig. 4a&b, respectively. Under the normal incidence of p-polarized propagating waves, we have \({P}_{\text{p},\text{g}\text{w},\text{r}\text{i}\text{g}\text{h}\text{t}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}=0.88\) [Fig. 4a]. Similarly, we have \({P}_{\text{s},\text{g}\text{w},\text{l}\text{e}\text{f}\text{t}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}=0.86\) [Fig. 4b] if the incident waves are s-polarized. These measured results clearly indicate the high coupling efficiency in the observed polarization-maintaining near-field directionality.
Actually, the incident light can be arbitrarily polarized and have a random mixture of p-polarized and s-polarized propagating waves. We thus also investigate the dependence of the coupling efficiency on the incident polarization state in Fig. 4a&b. Under the normal incidence, the incident polarization state is closely related to the angle \(\varphi\) between the direction of incident electric fields and the \(\widehat{x}\) direction; see the inset in Fig. 4b. By this definition, the incident light is the p-polarized (s-polarized) propagating wave if \(\varphi ={0}^{o}\) (\(\varphi ={90}^{o}\)), and the values of \({P}_{\text{p},\text{p}\text{w}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}\) and \({P}_{\text{s},\text{p}\text{w}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}\) would be highly dependent on \(\varphi\) as shown in Fig. 4a&b. Remarkably, the measured \({P}_{\text{p},\text{g}\text{w},\text{r}\text{i}\text{g}\text{h}\text{t}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}\) as a function of \(\varphi\) almost has a same variation tendency with that of \({P}_{\text{p},\text{p}\text{w}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}\), while the variation tendency of \({P}_{\text{s},\text{g}\text{w},\text{l}\text{e}\text{f}\text{t}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}\) agrees well with that of \({P}_{\text{s},\text{p}\text{w}}/{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}\). Moreover, we always have \(\frac{{P}_{\text{p},\text{g}\text{w},\text{r}\text{i}\text{g}\text{h}\text{t}}+{P}_{\text{s},\text{g}\text{w},\text{l}\text{e}\text{f}\text{t}}}{{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}}\ge 85\%\) [Fig. 4c] when \(\varphi\) changes from \({0}^{o}\) to \(9{0}^{o}\). This further indicates that the polarization-maintaining near-field directionality with high coupling efficiency can occur for arbitrary incident polarization states. This feature is also verified by the far-field measurement. From Fig. 4c, we always have \(\frac{{P}_{\text{s}\text{c}\text{a}\text{t}\text{t}\text{e}\text{r}\text{e}\text{d}}}{{P}_{\text{i}\text{n}\text{c}\text{i}\text{d}\text{e}\text{n}\text{t}}}\le 11.2\%\) for arbitrary value of \(\varphi\), where \({P}_{\text{s}\text{c}\text{a}\text{t}\text{t}\text{e}\text{r}\text{e}\text{d}}\) stands for all scattered power into propagating waves; see the calculation in Fig. S6. In addition, a minor proportion of power (around 2%) would be dissipated during the coupling, due to the material loss in the designed metasurface.