Asymmetric transmission—direction-selective electromagnetic transmission between two ports—is a phenomenon exhibited by two-dimensional chiral systems. The possibility of exploiting this phenomenon in chiral metasurfaces opens exciting possibilities for applications such as optical isolation and routing without external magnetic fields. This work investigates optical asymmetric transmission in chiral plasmonic metasurfaces supporting lattice plasmon modes and unveils its physical origins. We show numerically and experimentally that asymmetric transmission is caused by an unbalanced excitation of such lattice modes by circularly polarized light of opposite handedness. The excitation efficiencies of the lattice modes, and hence, the strength of the asymmetric transmission, are controlled by engineering the in-plane scattering of the individual plasmonic nanoparticles such that the maximum scattering imbalance occurs along one of the in-plane diffraction orders of the metasurface. Furthermore, we show that only the nonzero diffraction orders contribute to this effect. By highlighting the role of the localized plasmon modes supported by the nanoparticle and their radiative coupling to the lattice structure, our study provides a guideline for designing metasurfaces with asymmetric transmission enabled by lattice plasmons.
I. INTRODUCTION
Chirality refers to an intrinsic sense of handedness of a three-dimensional (3D) structure, which remains invariant regardless of the direction of observation. Because of this property, the structure and its mirror image cannot be brought into congruence by a translation and rotation operation.1 Proteins, chemical systems, and many biomolecules are known to be chiral, exhibiting contrasting functions and properties between its mirror-symmetric pairs (enantiomer).2–4 3D chiral systems yield a different optical response when interacting with either handedness of circularly polarized light, leading to phenomena such as circular dichroism and optical activity. Both phenomena are widely used for characterizing the optical response of chiral systems.5 Chirality could also be defined for two-dimensional (2D) structures. However, as planar structures do not possess an intrinsic handedness, chirality is defined by the inability of a 2D structure and its mirror image to be brought into congruence unless lifted from its plane.6 Systems composed of achiral structures are known to exhibit chiral effects, provided that the surface normal and the electric field vector form a chiral triad.6–8 This phenomenon, usually called extrinsic chirality, is indistinguishable from standard (intrinsic) chirality9 and has been reported in a variety of systems, including single nanoparticles.10,11
The sense of handedness of a 2D chiral structure is reversed when it is flipped with respect to an in-plane axis.9 Because of this property, 2D chiral structures can show asymmetric optical transmission, i.e., for a given circular polarization state, the optical transmittance depends on the side of the sample that is being illuminated.12–15 This is in compliance with the Lorentz reciprocity theorem and does not require magnetic fields as for the Faraday effect.12 The phenomenon can also be observed as a difference in reflection and/or a difference in absorption. Occasionally, the term circular dichroism has been used to refer to asymmetric transmission. However, it is essential to note that circular dichroism is exclusive to 3D chiral systems, while asymmetric transmission is exclusive to 2D chiral systems. Although circular dichroism has been reported in 2D arrays through extrinsic chirality, the experimental geometry in such cases tenders the structure 3D chiral.6–8 Another distinguishing feature of asymmetric transmission from circular dichroism is that the former changes the sign of the signal upon reversal of direction of illumination, whereas the latter remains invariant with respect to this operation.
Chiroptical phenomena have been extensively explored in plasmonic systems, such as individual nanostructures,16–19 2D array of plasmonic nanostructures (metasurfaces),20–23 and plasmonic nanostructures in suspension (metafluids).24 Of particular interest is the asymmetric transmission in chiral plasmonic metasurfaces. When the periodic interparticle spacing of such metasurfaces is smaller than the optical wavelength, the asymmetric transmission has been attributed to the simultaneous presence of anisotropy and losses.12,13,25 However, when the interparticle spacing is of the order of the optical wavelength, the metasurface can support lattice plasmon modes enabled by diffractive coupling of localized plasmons,26,27 which, in turn, can enable a different mechanism of asymmetric transmission. The asymmetric transmission of 2D chiral arrays supporting lattice plasmon modes has been reported recently;28,29 however, the origin of this phenomenon remains unclear.
In this work, we unveil the physical origin of the lattice-plasmon-induced asymmetric transmission in diffractive chiral metasurfaces. Through a numerical and experimental analysis, we show that the asymmetric transmission mechanism relies on different lattice plasmon mode excitation efficiencies for left circularly polarized (LCP) and right circularly polarized (RCP) lights. Such excitation efficiencies can be controlled by tailoring the nanoparticle’s in-plane distribution of scattered light for circularly polarized excitation and its alignment with the in-plane diffraction orders of the metasurface. Contrary to previous reports,29 our study shows that the lattice-plasmon-induced asymmetric transmission results solely from the contribution of nonzero diffraction orders.
For simplicity, we analyze this phenomenon for a metasurface composed of an array of achiral nanoparticles—a tetramer consisting of four nanoantennas, as shown in the inset of Fig. 1. We show that near the Rayleigh anomaly condition, where the asymmetric transmission effect is the strongest, the metasurface supports lattice plasmon modes that respond selectively to the polarization handedness, confirming the underlying asymmetric transmission mechanism enabled by the lattice plasmon modes. These results are of general relevance for understanding the role of diffracted waves in 2D chiral systems, particularly for the design of chiral plasmonic metasurfaces.
Schematic illustration of the metasurfaces under investigation. The metasurface is illuminated at normal incidence by circularly polarized light propagating in the positive z direction. The metasurface consists of a square array of gold tetramers embedded in a homogeneous medium. The arrangement and dimensions of the four nanoantennas composing the tetramer are shown in the inset. All dimensions are given in nanometers.
Schematic illustration of the metasurfaces under investigation. The metasurface is illuminated at normal incidence by circularly polarized light propagating in the positive z direction. The metasurface consists of a square array of gold tetramers embedded in a homogeneous medium. The arrangement and dimensions of the four nanoantennas composing the tetramer are shown in the inset. All dimensions are given in nanometers.
II. ASYMMETRIC TRANSMISSION MECHANISM AND NUMERICAL ANALYSIS
Consider the metasurface depicted in Fig. 1, which consists of a planar array of tetramer nanostructures. The tetramers are made of gold, arranged into a squared lattice, and surrounded by a homogeneous medium with the refractive index n = 1.51 (i.e., glass). The array lies on the (x and y) plane, having a lattice spacing Λ = 600 nm. The dimensions of each tetramer are indicated in the inset of Fig. 1. The metasurface is illuminated by a circularly polarized plane wave propagating in the positive z direction and impinging at a normal incidence onto the surface. Although the individual tetramers are achiral, it is possible to create a 2D chiral metasurface by rotating the tetramer at each lattice point around its center.
The metasurface can support lattice plasmon modes when Λ is of the order of the optical wavelength. These modes result from the coupling of light scattered from neighboring particles through grazing diffraction orders or Rayleigh anomalies.26,30 For a normally incident excitation, such Rayleigh anomalies occur at wavelengths given by31
where p and q are the integers indicating the in-plane diffraction orders along the x and y directions, respectively. The Rayleigh anomaly condition for the first diffraction order, , is degenerate and occurs at the wavelength λ1,0 = 906 nm. Similarly, the next Rayleigh anomaly, which corresponds to the diffraction order , occurs at λ1,1 = 640 nm. Lattice plasmon modes propagating along different in-plane diffraction orders are supported at wavelengths close to those described by Eq. (1). The direction of the three main diffraction orders is depicted in Fig. 1.
To understand the asymmetric transmission mechanism exhibited by this metasurface, we first analyze the optical scattering properties of the individual tetramers for LCP and RCP plane-wave illumination. For this, we use a fully vectorial Maxwell equation solver based on the finite-difference time domain (FDTD) method. The permittivity of gold was taken from Ref. 32. Figure 2(a) shows the obtained in-plane angular distribution of the scattered intensity at the wavelength λ = 620 nm, which is close to the particle’s resonance. Here, IL and IR are the in-plane scattered intensities for LCP and RCP excitations, respectively. The difference between these intensities, ΔI = |IR − IL|, will be referred to as the differential in-plane scattering. This quantity is plotted in Fig. 2(b) for the same wavelength as Fig. 2(a) and as a function of wavelength in Fig. 2(c). Observe that regardless of the optical wavelength, the differential in-plane scattering has maxima (ΔImax) and minima (ΔImin) at particular azimuthal angles whose values are odd and even integer multiples of 22.5°, respectively.
Optical response of the isolated nanostructure obtained through FDTD simulations. (a) Polar plot of the far-field scattered intensity as a function of in-plane angle for an individual nanostructure. IR and IL are the normalized scattered intensity for RCP (solid blue) and LCP (solid red) excitations. (b) Normalized plots of differential in-plane scattering, , for a 620 nm excitation wavelength. Superimposed on the plots is the isolated nanostructure for reference. (c) Variation in ΔI as a function of wavelength for the nanostructure in normalized units.
Optical response of the isolated nanostructure obtained through FDTD simulations. (a) Polar plot of the far-field scattered intensity as a function of in-plane angle for an individual nanostructure. IR and IL are the normalized scattered intensity for RCP (solid blue) and LCP (solid red) excitations. (b) Normalized plots of differential in-plane scattering, , for a 620 nm excitation wavelength. Superimposed on the plots is the isolated nanostructure for reference. (c) Variation in ΔI as a function of wavelength for the nanostructure in normalized units.
The individual tetramer exhibits a nonzero differential in-plane scattering in a spectral range overlapping with the Rayleigh anomalies of the array. Consequently, the lattice plasmon excitation efficiency will vary depending on whether the tetramer array is illuminated with RCP or LCP light, leading to different transmittances for each case. In turn, this leads to asymmetric transmission, as reversing the polarization handedness is equivalent to illuminating the metasurface from the opposite side. Moreover, since ΔI varies along the in-plane angle, the differential excitation of lattice plasmon modes by RCP and LCP lights also varies with the relative angle between ΔI and the lattice plasmon mode propagation direction. In particular, when ΔImax (ΔImin) is along the propagation direction of a lattice plasmon mode, the difference in the efficiency of excitation of that lattice plasmon mode will be maximized (minimized).
To analyze this effect, the transmittance of the metasurface as a function of tetramer’s rotation angles is computed and plotted in Fig. 3(a) for RCP illumination; a similar result is obtained for LCP illumination (not shown). These transmittance results obtained for RCP and LCP illuminations will be used for estimating the asymmetric transmission. The relative orientation of the nanostructure for different rotation angles with respect to the lattice is illustrated in Figs. 3(c)–3(f). Figure 3(b) shows the asymmetric transmission calculated from the numerical results as
where TR and TL are the transmittance for RCP and LCP excitations, respectively. The asymmetric transmission vanishes for wavelengths beyond 906 nm because the Rayleigh anomaly condition is not satisfied, and so, the lattice plasmon modes are not supported. On the other hand, the asymmetric transmission is maximized at wavelengths close to the Rayleigh anomaly conditions (λ10 and λ11), where lattice plasmon modes are most efficiently excited.
(a) Numerically calculated transmittance spectra for RCP excitation for the metasurface as a function of rotation of the individual nanostructure. The arrows at the top of the results denote the wavelengths at which the Rayleigh anomalies occur, corresponding to a lattice period of Λ = 600 nm in a homogeneous medium with n = 1.51. (b) Variation in asymmetric transmission calculated using Eq. (2) as a function of the tetramer’s rotation angle. Asymmetric transmission occurs when the maxima of the differential in-plane scattering are aligned along the propagation direction of the grazing diffraction order. (c)–(f) show the schematic illustration of the metasurface when the individual nanostructure at each lattice point is rotated around its center.
(a) Numerically calculated transmittance spectra for RCP excitation for the metasurface as a function of rotation of the individual nanostructure. The arrows at the top of the results denote the wavelengths at which the Rayleigh anomalies occur, corresponding to a lattice period of Λ = 600 nm in a homogeneous medium with n = 1.51. (b) Variation in asymmetric transmission calculated using Eq. (2) as a function of the tetramer’s rotation angle. Asymmetric transmission occurs when the maxima of the differential in-plane scattering are aligned along the propagation direction of the grazing diffraction order. (c)–(f) show the schematic illustration of the metasurface when the individual nanostructure at each lattice point is rotated around its center.
Most importantly, we note that the asymmetric transmission depends on the in-plane rotation of the individual tetramers. Clearly, the largest asymmetric transmission is obtained when the individual nanostructure is rotated by odd integer multiples of 22.5°, which coincides with the angles for ΔImax of the isolated nanostructure [see Fig. 2(b)]. The supplementary material gives a detailed analysis of the asymmetric transmission arising from an array of chiral nanostructures. From these results, we conclude that the value of the differential in-plane scattering of the individual nanostructure, ΔI, along the diffraction orders at the Rayleigh anomaly condition determines the strength of the asymmetric transmission. Furthermore, we note that when lattice plasmon modes are excited, the metasurface can exhibit asymmetric transmission as long as the individual nanostructure exhibits differential in-plane scattering along the lattice plasmon mode propagation direction. Hence, even for a nanostructure having fourfold rotational symmetry, the presence of higher-order diffraction and the associated losses for the higher-orders can result in broken rotational symmetry of the entire process,33 thereby making it irrelevant whether the nanostructure is chiral or achiral and whether it has fourfold rotational symmetry or not.
It is important to note here that isolated nanostructures exhibit differential in-plane scattering only when higher-order electric and magnetic moments are excited (see the supplementary material). Hence, there are three necessary conditions for asymmetric transmission in 2D chiral arrays composed of fourfold rotationally symmetric structures: (i) excitation of higher-order localized plasmonic resonances (resulting in-plane differential scattering), (ii) excitation of lattice plasmon modes, and (iii) an angular overlap between the maxima of in-plane differential scattering and the lattice plasmon mode propagation direction.
Finally, an analysis of the transmission characteristics for RCP and LCP excitations along different diffraction orders, away from the Rayleigh anomaly condition, reveals that the transmission into the zeroth-order is degenerate for both excitations irrespective of the wavelength and nanostructure rotation angles. Thus, the asymmetric transmission observed in Fig. 3(b) is caused solely by the difference in transmission of higher diffraction orders.
III. EXPERIMENTAL RESULTS
Based on the numerical analysis, the fabrication of tetramer metasurfaces with the required orientation of individual structures was carried out. Details about the fabrication process can be found in the supplementary material. A schematic illustration of the experimental setup used for the investigation is shown in Fig. 4(a). A filtered broadband light source is used to carry out spectrally resolved measurements in a range of 630–950 nm. A broadband polarizer followed by a quarter-wave-plate is used to prepare the polarization state of the excitation beam. The beam is then weakly focused to a diameter of ∼50 μm with the help of a convex lens (effective NA ), and the metasurface is precisely placed in the focal spot with the help of a three-axis piezo stage. A detailed description of the experimental setup can be found in the supplementary material.
(a) Sketch of the experimental setup. A broadband light source is filtered to perform spectrally resolved measurements. A broadband polarizer followed by a quarter-wave-plate is used to prepare the polarization state of the excitation beam. The beam is weakly focused (NA = 0.015), producing a focal spot of ∼50 μm diameter on the tetramer metasurface. The transmitted light is collected with a microscope objective of NA = 1.3. The objective’s back focal plane, containing the k-space information, is imaged onto a CCD camera. (b) The k-space image recorded in transmission at an input wavelength of 770 nm. Zeroth and first diffraction orders are visible, which were evaluated for contributions to asymmetric transmission. (c) For direct imaging of the tetramer metasurface, in reflection, at Rayleigh anomaly, we adjust the incident illumination and use a microscope objective with effective NA = 0.1 to illuminate and subsequently collect the reflected light.
(a) Sketch of the experimental setup. A broadband light source is filtered to perform spectrally resolved measurements. A broadband polarizer followed by a quarter-wave-plate is used to prepare the polarization state of the excitation beam. The beam is weakly focused (NA = 0.015), producing a focal spot of ∼50 μm diameter on the tetramer metasurface. The transmitted light is collected with a microscope objective of NA = 1.3. The objective’s back focal plane, containing the k-space information, is imaged onto a CCD camera. (b) The k-space image recorded in transmission at an input wavelength of 770 nm. Zeroth and first diffraction orders are visible, which were evaluated for contributions to asymmetric transmission. (c) For direct imaging of the tetramer metasurface, in reflection, at Rayleigh anomaly, we adjust the incident illumination and use a microscope objective with effective NA = 0.1 to illuminate and subsequently collect the reflected light.
We consider three metasurfaces with in-plane rotation of nanostructures at 22.5°, −22.5°, and 45° to measure the lattice-plasmon-induced asymmetric transmission. For the present metasurface designs, a significant fraction of the transmitted power remains in the zeroth-order, hindering the experimental observations of the effect. Hence, we perform diffraction-resolved measurements to identify the independent power contributions of the transmitted zeroth and first diffraction orders. To this end, the transmitted light is collected with a microscope objective (NA = 1.3), with its back focal plane imaged onto a CCD camera. The back focal plane gives us access to the wave-vector space (k-space) of the transmitted light from which we obtain the diffraction-resolved power contributions (see the supplementary material). The zeroth-order is measured for all wavelengths. The first diffraction orders can be measured only for wavelengths shorter than 780 nm because of the limitation imposed by the maximum collection angle of the microscope objective (60.25°). This limitation prevents us from measuring the power coupled to lattice plasmon modes (near the Rayleigh anomaly wavelengths). For such cases, we employ direct imaging of the metasurface in reflection via leakage radiation microscopy. We adjust the incident illumination and use a microscope objective with effective NA = 0.1 to illuminate and subsequently collect the reflected light. This objective selection in reflection consequently results in a focused spot size smaller than in the transmission case. The use of slightly higher NA for imaging in reflection helps in clearly observing the leaky propagating first diffraction orders, as shown in Fig. 4(c). This technique allows us to obtain qualitative information of the power coupled to lattice plasmon modes by visualizing light leaked as they propagate. Figures 4(b) and 4(c) illustrate, respectively, k-space and leakage radiation microscopy images recorded using the above mentioned techniques.
The experimental measurements of lattice-plasmon-induced asymmetric transmission obtained for the achiral (45° rotation) and chiral (22.5° and −22.5° rotations) metasurfaces are summarized in Fig. 5. The diffraction-resolved asymmetric transmission for the zeroth and first diffraction orders is plotted in Fig. 5(a) as a function of the wavelength; both measurements (dashed lines) and simulations (solid lines) are included in the plots. Diffraction-resolved asymmetric transmission is calculated using Eq. (2), but TR and TL are replaced with and , which are the transmittance in the diffraction order for RCP and LCP excitations, respectively. Leakage radiation microscopy images are used to study the phenomenon at a wavelength close to the Rayleigh anomaly (marked by gray dashed line in spectral plots), by constructing differential measurements obtained from a pixel-by-pixel subtraction of two camera shots taken for excitation with opposite polarization handedness. Such differential leakage radiation microscopy images, as shown in Fig. 5(b), allow us to visualize the asymmetric excitation of lattice plasmon modes responsible for the asymmetric transmission.
Experimental and simulated diffraction-resolved asymmetric transmission results for achiral (45° rotation) and chiral (22.5°–22.5° rotations) metasurfaces. (a) Wavelength vs diffraction-resolved asymmetric transmission for zeroth (top panel) and first diffraction order (bottom panel) for the achiral (green) and chiral (blue and red) metasurfaces. (b) Differential leakage radiation microscopy (LRM) images recorded at 905 nm wavelength, showing asymmetric excitation of lattice plasmon modes for chiral (22.5° and 22.5° rotations) metasurfaces with opposing signs and symmetric excitation of lattice plasmon modes for the achiral metasurface (45°).
Experimental and simulated diffraction-resolved asymmetric transmission results for achiral (45° rotation) and chiral (22.5°–22.5° rotations) metasurfaces. (a) Wavelength vs diffraction-resolved asymmetric transmission for zeroth (top panel) and first diffraction order (bottom panel) for the achiral (green) and chiral (blue and red) metasurfaces. (b) Differential leakage radiation microscopy (LRM) images recorded at 905 nm wavelength, showing asymmetric excitation of lattice plasmon modes for chiral (22.5° and 22.5° rotations) metasurfaces with opposing signs and symmetric excitation of lattice plasmon modes for the achiral metasurface (45°).
First, we discuss the case of the metasurface with tetramers rotated by 45° in-plane [see Fig. 5(a) top-left], which constitutes an achiral metasurface. As expected, this metasurface results in vanishing asymmetric transmission in first diffraction orders, which is confirmed by both experimental and simulation results. The differential leakage radiation microscopy image recorded for this metasurface also confirms vanishing asymmetric propagating modes [see Fig. 5(b), marked in green box]. Next, consider the metasurface constituted by tetramers with in-plane rotations of 22.5° and −22.5° [see Fig. 5(a) top-center and top-right, respectively]. As expected from simulation results shown in Fig. 3(d), the two corresponding cases show opposing signs of asymmetric transmission for the first-order. This is confirmed by the experimental results, as shown in the spectral plot in Fig. 5(a). The differential leakage radiation microscopy image recorded for the chiral arrangements depicts asymmetrically excited lattice plasmon modes with opposing signs [see Fig. 5(b), marked in blue and red boxes].
To demonstrate and validate the sole contribution of the first diffraction orders to the observed asymmetric transmission, we also show the corresponding simulation and experimental results for the zeroth diffraction order not featuring any significant spectral dependence for all three studied metasurfaces [see Fig. 5(a) center]. The absence of asymmetric transmission for all three tetramer metasurfaces for the zeroth-order underscores the essential role of lattice plasmon modes in the asymmetric transmission mechanism.
IV. CONCLUSIONS
In conclusion, we studied the role of lattice plasmon modes on the phenomenon of asymmetric transmission in chiral 2D arrays of plasmonic nanoparticles and explained its physical origin. We showed that the difference in the nanoparticle’s in-plane scattering of LCP and RCP lights can contribute to asymmetric transmission in 2D metasurfaces through the excitation of polarization-handedness selective lattice plasmon modes. The difference in the nanoparticle’s in-plane scattering is attributed to the excitation of higher-order localized plasmonic resonances. The different excitation efficiencies of the lattice plasmon modes, produced by unbalanced in-plane scattering of LCP and RCP lights and its coupling to grazing diffraction orders of the periodic lattice, result in asymmetries in transmission, reflection, and absorption. As the difference in the in-plane scattered intensity varies as a function of the in-plane angle around the nanostructure, rotating the nanostructure in the plane of the metasurface can act as a handle to control the metasurface’s asymmetric transmission. Using diffraction-order-resolved transmission measurements, we demonstrated that the asymmetric transmission results solely from contributions of higher-order diffracted waves, while the zeroth-order transmission is not asymmetric. Furthermore, we showed that, even though the isolated nanostructure has a fourfold rotational symmetry, diffractive metasurfaces can exhibit asymmetric transmission for normal incidence illumination within the spectral range for which lattice plasmon modes are supported. This behavior is contrary to what occurs in a non-diffractive chiral metasurface. Our study sheds new light on asymmetric or chiral phenomena in metasurfaces and the critical role of diffraction effects such as the excitation of lattice plasmon modes.
SUPPLEMENTARY MATERIAL
Details about the fabrication process, experimental setup, and analysis of the asymmetric transmission arising from gammadion metasurface can be found in the supplementary material.
ACKNOWLEDGMENTS
I.D.L. acknowledges the financial support from the Federico Baur Endowed Chair in Nanotechnology. R.W.B. acknowledges support from the Natural Sciences and Engineering Research Council of Canada, the Canada Research Chairs program, the U.S. ARO under Grant No. W911NF-18-1-0337, and U.S. DARPA under Award No. W911NF1810369.
AUTHOR DECLARATIONS
Conflict of Interest
The authors declare no conflicts of interest.
DATA AVAILABILITY
The data that support the findings of this study are available within the article and its supplementary material.