## 1 Introduction

Electromagnetic fields which carry chirality – in their simplest form, left- and right-circularly polarized plane waves – deserve huge interest as they interact with matter chirality, enabling for instance to discriminate enantiomers in chemistry, which are ultimately connected with key features of living organisms. Indeed, many biomolecules have a specific handedness (*homochirality*), and it is not yet clear why nature has decided to go in that precise direction [1]. From a more application-oriented point of view, the pharmaceutical industry constantly seeks for effective methods to discriminate stereoisomers, an application where chiral light-matter interaction could prove useful.

To date, the most common technique to prepare and analyze chiral light is to employ birefringent plates and linear polarizers that convert light to and from linear polarization, as the technology of direct sources and detectors of chiral light is still in its infancy [2], [3], [4], [5]. Last advances in nanotechnology are however revolutionizing chiral optical devices [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17]. Far- and near-field chiral electromagnetic responses have been indeed observed in a variety of artificially structured systems, where the shape of the machined elements must have a three-dimensional character if geometric and electromagnetic chirality has to be attained in its most rigorous form, because of the requirement of the absence of any mirror symmetry plane [18], [19], [20], [21]. Several proposed structures hence rely on volumetric fabrication techniques, which however suffer from either scarce throughput or limited flexibility. Thus, less demanding fabrication technologies – i.e. planar technologies – may also be employed, as witnessed by some reports [22], [23], [24]. For instance, a dielectric film with a non-centrosymmetric partially etched planar pattern was proposed as a simple gateway towards strong chiro-optical phenomena [25]. Indeed, there is a wide interest in developing subwavelength-patterned high-index dielectrics, both for applications and for fundamental research: from one side, they enable the synthesis of flat lenses, polarimeters, spectrometers, nonlinear components and computer-generated holograms [26], [27], [28], [29], [30], [31], [32], [33]; from another side, they exhibit a variety of intriguing phenomena such as Fano lineshapes, perfect forward scattering, geometric phase effects, and bound states in the continuum resonances [34], [35], [36], [37], [38], [39], [40].

Most of these effects arise from complex electromagnetic behaviors whose essence can be grasped understanding the interplay between localized and delocalized resonances. In other words, these systems often rely on the co-presence, and on the competition, between guided wave phenomena and antenna-like responses, as highlighted by the prototypical transition between the guided mode filter regime, the photonic crystal regime, and the independent-particle Mie-scattering regime [41]. In this work we have reported on the observation of photonic bands in a chiral metasurface, highlighting that the chiral response shows fingerprints of both guided wave and locally resonant phenomena. The object under investigation is the simplest conceivable 2d-chiral dielectric metasurface: a slab perforated with L-shaped holes [42]. In general, a 2d-chiral object is defined as follows: given the family П of planes perpendicular to a given plane *π_{z}* (the subscript

*z*refers to the fact that

*π*is usually chosen to be a constant-

_{z}*z*plane), an object is 2d-chiral if it cannot be superimposed with itself after a mirror reflection for any plane belonging to П. A 2d-chiral object is also 3d-chiral if it cannot be superimposed with itself after a mirror reflection for any plane parallel to

*π*. An object that is 2d-chiral but not 3d-chiral has a mirror symmetry plane of the type

_{z}*π*. 2d-chiral objects may have vanishing or finite thickness along

_{z}*z*. Our metasurface has a finite thickness; this property allows it to implement interesting chiral electromagnetic features such as chiro-optical far-field response at normal incidence [19]. To our knowledge, our work is the first report of a 2d-chiral, 3d-achiral, dielectric metasurface; moreover, it is the first report of angularly-resolved spectroscopy of 2d-chiral, 3d-achiral metasurfaces. The metasurface under analysis also exhibits superchiral near-fields with incident unpolarized light. The response at normal incidence will be analyzed by means of a singular-value decomposition approach that reveals the operational capabilities of any 2d-chiral photonic device. In addition, we have shown that the extremely simple L-shaped structure is also a

*minimal*one: by tuning a small number of geometrical parameters, it is possible to access a very wide set of the transmission matrices characteristic of a general 2d-chiral metasurface. In this sense, we approached the solution of the inverse design problem of 2d-chiral metasurfaces [43].

## 2 Chiral response at normal incidence

The metasurface under investigation is illustrated in Figure 1A. It consists of a 220 nm thick gallium arsenide membrane, patterned with L-shaped holes arranged over a square lattice. The geometric parameters of the holes are reported in Figure 1B, superimposed to the scanning electron microscope (SEM) image of a fabricated sample. While in forthcoming analysis the parameters *a*, *f*_{1}…*f*_{4} are allowed to vary, the first part of the article deals with a specific choice of parameters, which correspond indeed to the SEM image. In detail, we have *a*=1134 nm, *f*_{1}=0.76, *f*_{2}=0.58, *f*_{3}=*f*_{4}=0.327. The sample fabrication process, whose details are reported in the Supplementary Materials, allows to obtain a frame-supported membrane, freely accessible from both sides to perform optical measurements with a moderately focused (≈50 μm spot) near-infrared (1–2 μm wavelength) beam. The beam originates from a supercontinuum source, filtered by means of an acousto-optic tunable filter yielding a spectral bandwidth of ≈2 nm. The system is controlled through an automated software developed in C++. The polarization state is prepared with a Glan-Taylor polarizer followed by a *λ*/4 superachromatic waveplate; no polarization analysis was performed on the beam transmitted after the sample. The sample was mounted on a rotating support that enables to measure the angularly-resolved transmittance, as illustrated in Figure 1A. A first set of measurements has been collected at normal incidence (α=0). The experimental data are reported in Figure 1C, and are compared with the outcome of a numeric model (rigorous coupled wave analysis, RCWA – see the Supplementary Materials for details). The spectra consist of a series of quite narrow dips, some of them having different shapes and depths depending on the polarization state of the incident light. From the transmittances for left (right) circularly polarized light, respectively *front incidence* and *back incidence* have been collected from the same sample, but illuminated from either the top surface or the bottom surface. The trace *enantiomer* derives instead from a second sample, which differs from the first one solely for the L-shaped hole sense of chirality at fabrication stage: the second sample, when observed from the front side, is described by the same parameters of the first sample exception made for the exchange *f*_{1}↔*f*_{2}. A good matching between the CPDT spectrum of the enantiomer with that of the original sample’s back side, as well as the fact that the front and back side CPDT are opposite to each other, indicates that the objects under investigation are almost perfectly 2-d chiral. Indeed, while in principle a through-hole perforated homogeneous membrane has a mirror plane parallel to its surface, fabrication imperfections may have lead the actual sample to deviate from this ideal characteristic. The data indicate instead that the two samples are essentially the same object, from both geometric and electromagnetic points of view.

The data reported so far indicate that a 2d-chiral metasurface with a very simple design exhibits strong circular polarization differential transmittance in a narrow band close to the optical telecommunication window. It should be emphasized that the observed CPDT is not accompanied by absorption; rather, it originates from a redistribution of the incident energy among the transmitted and reflected beams (no diffraction is present, thanks to the subwavelength dimensions of the pattern). In this sense, the presented object is rather a beam splitter sensitive to the circular polarization, and the observed CPDT is inherently different from the absorption circular dichroism (or circular diattenuation, CD), properly occurring, for instance, in isotropic solutions of chiral molecules. However, from an observer point of view who can only access the transmitted beam, CPDT is indistinguishable from CD. Hence, our metasurface can be employed to mimick the CD of target objects, such as naturally occurring molecules.

## 3 Chiral response at oblique incidence and band structure

In the transmission spectra reported in Figure 1, a rich structure can be noticed, whose physical origin deserves attention *per se* and in view of applications. In order to get further insight into the nature of the resonances leading to CPDT peaks, a powerful method is to measure the metasurface transmittance at different angles of incidence (i.e. for different orientations of the incident wavevector). Angularly-resolved measurements have been recently employed to reveal special features of polarization phenomena like magneto-optic effects in quasi-ordered structures [44] and asymmetric transmission in low-symmetry plasmonic hole arrays [45]. With angularly-resolved measurements, one may distinguish between dispersive and non-dispersive resonances, which originate from different physical phenomena: guided mode resonances and antenna-like resonances, respectively. In Figure 2A and B, we have reported the CPDT mapped in the frequency-wavevector space for our metasurface. Here, the in-plane wavevector is the projection of the incident light wavevector on the metasurface plane. In the experiment only positive angles have been studied, while the model data are available for both positive and negative angles. The good matching between experimental and model data allows to acquire confidence in the model as a whole, which shows interesting features. First, strongly dispersive bands, with positive and negative slope, suggest that the metasurface optical response has important contributions from guided mode resonances (also known as quasi-guided modes); this picture is supported by an empty-lattice band-folding model illustrated in the Supplementary Material. These resonances are responsible for the narrow dips observed in the transmission spectra and for the sharp oscillations in the CPDT spectra.

Second, by comparing the CPDT observed at opposite angles, it can be noticed that in the region at the center of the map (−20°<α<20°, 190 THz<ν<210 THz) the CPDT is not symmetric with respect to the exchange α↔−*α*. The effect can be better noticed looking at the curves in Figure 2C, where the CPDT as a function of the incidence angle has been reported for three fixed values of light frequency (labels i–iii in the figure). We attribute this effect to the presence of a spectrally and angularly broad chiral resonance, which we referred to as L-hole resonance. This resonance has a behavior that contrasts in two ways that of the guided mode resonances: first, the broad angular response of the L-hole resonance suggests that it has a spatially localized nature, as opposed to the delocalized, traveling-wave nature of the guided mode resonances (which instead follow a dispersive energy-wavevector curve). Second, the L-hole resonance induces the property CPDT(α)≠−CPDT(−α) (observed in the central region of the map, around point A in Figure 2B), while the guided-mode resonances display CPDT(α)=−CPDT(−α) (dispersive bands, point B).

Seeking for further confirmation of this view, we have calculated (see the Supplementary Materials for details) the optical near fields of the metasurface unit cell, upon the illumination conditions that excite the two types of resonances. The data are plotted in Figure 2D, where both the energy density (color map) and the Poynting vector (arrows) are reported. While the energy density distributions do not show remarkable behaviors, the Poynting vector field shows interesting properties that are consistent with the picture sketched above. When the metasurface is illuminated with the energy-wavevector pair labeled as “A” in Figure 2B, the Poynting vector in the unit cell has a strongly inhomogeneous distribution: for left polarized excitation, it even “winds up” around the L-shaped hole. For right-polarized excitation the Poynting vector has a different, yet irregular, distribution. We have interpreted these irregular distributions as arising from the L-hole resonance: a localized resonance dominated by multipoles, similar to the reports of [25]. On the contrary, when the metasurface is illuminated with the energy-wavevector pair labeled “B”, the Poynting vector shows much more regular distributions, typical of a guided mode that propagates parallel to the slab. Noticeably, the direction of the Poynting vector is opposite with respect to that of the in-plane wavevector, consistently with the negative dispersion of the photonic band where point B lies on.

## 4 Superchiral near-fields

Besides far-field angularly-dispersive chiro-optic response, which is of interest – for instance – for filters, holograms, and multiplexers, near-field optical chirality plays a fundamental role in the interaction with chiral matter. First relegated to the role of pure mathematical curiosity, electromagnetic (e.m.) chirality is now regarded as a crucial quantity to be taken care of when the optical detection of chiral molecules is under investigation [46], [47], [48], [49], [50]. The ability of a photonic structure to enhance the e.m. chirality is quantified by the *normalized e.m. chirality*, defined as *c* is the speed of light, **E** and **B** are the fields at the point of interest, and **E**_{inc} is the electric field vector of the incident plane wave. In the absence of any photonic structure, the plane wave would freely propagate and one would have *superchiral field*). Figure 3 shows that the 2d-chiral L-shape patterned membrane has a strong near-field e.m. chirality, also at normal incidence, and also when considering the spatial average of *λ*=1.45 μm. We have attributed these anomalies to the influence of the chiral localized resonance, “mode A”, discussed above, whose interplay with the narrower guided mode resonances leads to such features. In the considered spectral region other interesting phenomena take place. First, the metasurface displays

## 5 Singular value decomposition analysis of the metasurface operation

We will now focus back to the far-field metasurface response, in the attempt to understand in a comprehensive way its amplitude, phase, and polarization response. This analysis relies on the properties of the *T*-matrix (transmission matrix, or Jones matrix), where the device response concerning light transmission is completely encoded.^{1} The *T*-matrix can be written over different bases; the most common choices being that of linearly or circularly polarized waves. The following Section identifies the *T*-matrix written in these bases as *T*_{L} and *T*_{C}, respectively. For instance, *T*_{C} operates over Jones vectors whose elements are the right- and left-handed circularly polarized components of incident/transmitted light: **t**=*T*_{C}**i**, with (**t**=*t*_{R},*t*_{L})and **i**=(*i*_{R},*i*_{L}).

The form of *T*, and hence the possible operations that a metasurface can implement, are dictated by geometrical symmetry properties of the pattern lattice and unit cell. For instance, the L-shaped hole structure belongs to the more general category of objects that have a single mirror symmetry plane *π_{z}* (see Section 1); these objects are usually said to belong to the

*M*symmetry class. Under normal incidence, the

_{x,y}*T*

_{L}matrix of such metasurface is symmetric:

*T*

_{C}, in particular to its singular value decomposition (SVD). The SVD is an algebraic operation that reveals the structure of non-unitary linear operators, similarly to what eigenvalue decomposition (ED) does about unitary operators; a short recall of the properties of the SVD are given in the Supplementary Material. In the case under analysis we believe that the SVD is better suited with respect to ED, since the

*T*-matrix is in general non-unitary. Indeed, a dielectric metasurface is lossless in its full response (reflection plus transmission), but it appears lossy as the sole transmission is considered. The SVD of

*T*naturally describes the energy-handling behavior of the metasurface, since the square of its two singular values are the maximum and minimum wave intensity that can be transmitted by the metasurface (assuming unit intensity incident beam) [51]. Thus, our approach goes beyond what was reported in past studies on planar chiral metamaterials which did not rely on the SVD [52], [53], although an application of SVD to Müller-Jones calculus, not specific for chiral materials, is reported in [54]. We thus believe that SVD of chiral metasurfaces may provide a complementary viewpoint with respect to other decomposition methods (pseudopolar and integral decomposition [55], [56]) employed in polarization optics.

Algebraic manipulations (see the Supplementary Material) leveraging on the symmetry of *T*_{L} imply that the unitary matrices *V* and *W* appearing in the SVD of *T*_{C}, i.e. *T*_{C}=VΣW^{†}, must fulfil *σ*_{1,2}, and that the notation *V* and *W* implies that these matrices can be parametrized in the following insightful form:

where *F*, transformed by the basis change identified by the unitary matrix *D*_{1}. Noticeably, *D*_{1} is the matrix that describes, in the circular polarization basis, a rotation of the reference system about the axis perpendicular to the metasurface plane. In other words, the operation of an *M _{x,y}*-symmetric metasurface can be summarized in the operation encoded by the operator

*M*-symmetric metasurface, it is sufficient in essence to specify its

_{x,y}*F*: after that, the full space of transmission matrices can be saturated by simply rotating the metasurface about the axis perpendicular to the metasurface plane. It should be noticed that the above parametrization involves six real parameters, consistently with the fact that we are describing the symmetric Jones matrix

*T*

_{L}.

We now pay attention to understand the meaning of the parameters entering the matrices introduced above. To this end one should explicitly describe the operation of *T*_{C} over a generic incident field vector **i**. By decomposing it over the basis defined by the columns **w**^{(1,2)} of the unitary matrix *W*, one has **w**^{(1,2)}, reverses the handedness of these states (by converting **w**^{(1,2)}→**v**^{(1,2)}), and scales the energy content of each polarization state, according to the singular values *σ*_{1,2}. The location of the polarization states is connected with the parameters *ψ* and *θ* that appear in the matrix decomposition reported above (Figure 4). As expected, *ψ* is connected to a physical rotation (a rotation in the plane of the coordinate axes identifying the *x*- and *y*-related Stokes parameters *S*_{1,2}). Instead, *θ* is directly connected with the “chiral” Stokes parameter *S*_{3}. In the following, we mostly dropped the use of *θ* and rather defined the metasurface functionality referring to the *S*_{3} values of

To complete the picture one should also consider the phase response of the metasurface, which is described by the parameters *R* components of **w**^{(1,2)}, since

The algebraic and parametric analysis performed so far is valid at each individual wavelength. More information can however be obtained by studying the spectral dispersion of the SVD parameters, whose frequency-dependent response can be eventually correlated with quantities of more direct experimental access. In Figure 5A we have plotted the calculated spectral dispersion of the singular values for the L-shaped hole array described in the previous sections. The data follow from numerical SVD of the Jones matrix calculated from RCWA (see Supplementary Material); an experimental measurement of the SVD parameters is possible, in principle, relying on Müller-Jones polarimetry (exception made for the determination of *ϕ* that would require interferometry). In the analysis of Figure 5A singular values are not sorted in decreasing order, rather, they are sorted such as their trend with respect to the wavelength is smooth. The spectrum reveals a structure with narrow dips, mostly occurring in pairs, which stand out of a background where both *σ*’s are close to 1 (for *λ*<1.45 μm) or where one *σ* is close to 1 and the other is close to 0.5 (*λ*>1.45 μm). Noticeably, these dips occur at the same wavelengths where the spectra of transmission, circular polarization differential transmittance, and near-field chirality also show peaks or dips. It is also interesting to notice that the dips reach zero at those wavelengths, the metasurface does not transmit the radiation which is incident with the polarization dictated by the corresponding right singular vector. Figure 5B reports the spectral dependence of the third Stokes parameter of the right singular vectors. Here it can be noticed that the *σ*_{1,2}, are present. Finer features, often occurring close to the points where *σ*_{1}=*σ*_{2}, have an unclear interpretation and may be connected to singular value degeneracy. In the wavelength region around 1.45 μm a reversal of the signs of *σ*_{1,2} and on *ϕ* and 2Δ, which also show a rich behavior, with rapid changes occurring close to the zeros of the singular values. In general, the spectral behavior of the SVD parameter is rather complex. While this fact, from one side, opens further challenges in view of a comprehensive interpretation, it also holds promises in view of metasurface engineering: in essence, there is a wide space where the metasurface functionality can be chosen and designed.

## 6 Inverse problem and metasurface minimality

Provided with the SVD formalism, and motivated by the conclusions of the previous Section, an intriguing question is whether it is possible to solve the inverse problem: given a target metasurface function (i.e. a target *T* matrix, or better a set of target parameters *σ*_{1,2}, *ϕ*, and *ψ*), does it exist a metasurface unit cell shape that induces it? We have determined that appropriately shaped L-holes in a 220 nm thick dielectric slab with a refractive index n=3.38 (i.e. that of GaAs at 1.55 μm), are capable of implementing an arbitrary *T* matrix of the class pertinent to *M _{x,y}* objects. In this analysis the technologically relevant case of normal incidence is assumed. We approached the inverse problem by choosing a large set (120 combinations) of target parameters; the complete list is reported in the Supplementary Material. By running an optimization algorithm we have been able to identify, for all the targets, a specific geometry of the L-shaped hole such that the calculated

*T*matrix matches the target. It is interesting to notice that the L-shaped hole has 5 degrees of freedom (

*a*,

*f*

_{1}…

*f*

_{4}), the same as the number of non-trivial degrees of freedom of the

*T*matrix of

*M*objects:

_{x,y}*σ*

_{1,2},

*ϕ*. In this sense we have defined the metasurface as

*minimal*. We recall that the parameter

*ψ*is trivial in the sense defined in the previous Section, since it can be targeted by simply rotating the metasurface. While the complete solution dataset is available in the Supplementary Material, we have reported here four exemplary cases. Three of them are of immediate and intuitive interpretation; they are illustrated in Figure 6A–C. The target objects are polarizers, which are characterized by

*σ*

_{1}=1,

*σ*

_{2}=0. The first (A) is an ordinary circular polarizer, the second (B) is a circular polarizer that also imprints a phase delay, and the third (C) is an elliptical polarizer. More precisely, the objects under consideration are polarizers

*and*handedness inverters: thanks to the special structure of the singular vectors, a

*M*object with

_{x,y}*σ*

_{1}=1 and

*σ*

_{2}=0 eliminates one of the polarization components of the incident field, and transmits the other with inverted handedness (see Figure 4 and the related discussion). The fourth example, Figure 6D, is an implementation of the most general type of

*T*-matrices that can be targeted by means of

*M*objects. Here, all the target parameters (

_{x,y}*σ*

_{1,2},

*ϕ*) assume non-trivial values. The input and output singular states are elliptic polarization states, there is a phase difference between the input singular states, there is a phase difference between the output singular states, and the singular values are different from each other, and different from 0 to 1. In other words the metasurface behaves as the combination of a diattenuator and a retarder acting over elliptical states. Implementing this operation by means of conventional optical elements would require a bulky stack of layers, which are here replaced by a simple, subwavelength film where the operation is fully encoded in the shape of the lithographically defined holes. Hence, one might envisage ultracompact and ultralight optical components acting over all the key parameters of light – amplitude, phase and polarization. What is more, the possibility to access the full

*T*-matrix parameter space by means of a simple hole-shape tuning of a fixed-thickness membrane allows to implement space-variant metasurfaces capable of performing more advanced operations with respect to what is known to date [30] with possible applications to beam shaping, holography, and cryptography. As a final remark we have highlighted the targeting procedure, which was performed at a wavelength of 1.55 μm that is of direct interest to the telecommunication technology, is fully scalable, thanks to the scale invariance of Maxwell’s equations.

## 7 Conclusions

In conclusion we have reported the observation of various chiro-optical phenomena occurring in a 2d-chiral patterned dielectric slab, i.e. in a patterned slab that exhibits a single mirror symmetry plane parallel to the slab itself. Circular polarization differential transmittance is arranged in dispersive bands, and shows the fingerprints of localized and delocalized photonic resonances. It is in particular the effect of a localized resonance that induces an intrinsic chiral response on the metasurface, with consequences on the near-field chirality that shows superchirality and other anomalous behaviors. Relying on the singular value analysis of the transmission matrix, we identified the key parameters describing in full the operation of a 2d-chiral slab. We also noticed connections between the singular value spectrum and the above cited phenomena. Finally, we showed that the L-shaped hole structure, i.e. the most intuitive 2d-chiral pattern, is also a minimal one, as it allows to implement arbitrary transmission matrices with a minimal number of parameters. This result about inverse design might open the way towards advanced space-variant metasurfaces that exploit in full explore the phase, amplitude, and polarization degrees of freedom of light.

We like to acknowledge Alberto Bordin, who participated in an early stage of the project, Stefano Luin (Scuola Normale Superiore, Pisa) for the support on data fitting, Francesca Bontempi (Scuola Superiore Sant’Anna, Pisa) for ellipsometric measurements of the dielectric thin film, Sara Nocentini and Lorenzo Pattelli (LENS, Firenze) for their useful discussions and precious support for the spectroscopic measurements. This work was in part supported by the European Commission through the project PHENOMEN (H2020-EU-713450).

^{}

**Author contributions:** S.Z. has conceived the research and performed the electromagnetic simulations and the SVD analysis. S.Z. has fabricated the metasurface, from a thin film grown epitaxially by G.B. The spectroscopic measurements have been performed by S.Z. with the support of G.M. and F.R. The data have been interpreted by all the authors. The manuscript has been written by S.Z., with comments by all the authors.

## References

- [1]↑
Lough WJ, Wainer IW. Chirality in natural and applied sciences. Oxford, UK, Blackwell Science, 2002.

- [2]↑
Lobanov SV, Tikhodeev SG, Gippius NA, et al. Controlling circular polarization of light emitted by quantum dots using chiral photonic crystal slabs. Phys Rev B 2015;92:205309.

- [3]↑
Dyakov SA, Semenenko VA, Gippius NA, Tikhodeev SG. Magnetic field free circularly polarized thermal emission from a chiral metasurface. Phys Rev B 2018;98: 235416.

- [4]↑
Konishi K, Nomura M, Kumagai N, Iwamoto S, Arakawa Y, Kuwata-Gonokami M. Circularly polarized light emission from semiconductor planar chiral nanostructures. Phys Rev Lett 2011;106:057402.

- [5]↑
Söllner I, Mahmoodian S, Hansen SL, et al. Deterministic photon–emitter coupling in chiral photonic circuits. Nat Nanotechnol 2015;10:775–8.

- [6]↑
Pfeiffer C, Grbic A. Bianisotropic metasurfaces for optimal polarization control: analysis and synthesis. Phys Rev Appl 2014;2:044011.

- [7]↑
Zhao Y, Belkin M, Alù A. Twisted optical metamaterials for planarized ultrathin broadband circular polarizers. Nat Commun 2012;3:870.

- [8]↑
Zhao Y, Askarpour AN, Sun L, Shi J, Li X, Alù A. Chirality detection of enantiomers using twisted optical metamaterials. Nat Commun 2017;8:14180.

- [9]↑
Hentschel M, Schäferling M, Duan X, Giessen H, Liu N. Chiral plasmonics. Sci Adv 2017;3:e1602735.

- [10]↑
Kong X-T, Besteiro LV, Wang Z, Govorov AO. Plasmonic chirality and circular dichroism in bioassembled and nonbiological systems: theoretical background and recent progress. Adv Mat 2018;1:1801790.

- [11]↑
Asadchy VS, Díaz-Rubio A, Tretyakov SA. Bianisotropic metasurfaces: physics and applications. Nanophotonics 2018;7:1069.

- [12]↑
Tullius R, Karimullah AS, Rodier M, et al. “Superchiral” spectroscopy: detection of protein higher order hierarchical structure with chiral plasmonic nanostructures. J Amer Chem Soc 2015;137:8380–3.

- [13]↑
Schaeferling M, Engheta N, Giessen H, Weiss T. Reducing the complexity: enantioselective chiral near-fields by diagonal slit and mirror configuration. ACS Photonics 2016;3: 1076–84.

- [14]↑
Vázquez-Guardado A, Chanda D. Superchiral light generation on degenerate achiral surfaces. Phys Rev Lett 2018;120:137601.

- [15]↑
Shaltout A, Liu J, Kildishev A, Shalaev V. Photonic spin Hall effect in plasmon metasurfaces for on-chip chiroptical spectroscopy. Optica 2015;2:860–3.

- [16]↑
Liu L, Yang D, Wan W, Yang H, Gong Q, Li Y. Fast fabrication of silver helical metamaterial with single-exposure femtosecond laser photoreduction. Nanophotonics 2019;8:1087.

- [17]↑
Zhang R, Zhao Q, Wang X, Gao W, Li J, Tam WY. Measuring circular phase-dichroism of chiral metasurface. Nanophotonics 2019;8:909.

- [18]↑
Plum E, Liu X-X, Fedotov VA, Chen Y, Tsai DP, Zheludev NI. Metamaterials: optical activity without chirality. Phys Rev Lett 2009;102:113902.

- [19]↑
Plum E, Zheludev NI. Chirality and anisotropy of planar metamaterials. In: Maradudin AA, editor. Structured surfaces as optical metamaterials. Cambridge, UK, Cambridge University Press, 2011:94–157.

- [20]↑
Fernandez-Corbaton I, Fruhnert M, Rockstuhl C. Objects of maximum electromagnetic chirality. Phys Rev X 2016;6:031013.

- [21]↑
Garcia-Santiago X, Burger S, Rockstuhl C, Fernandez-Corbaton I. Measuring the electromagnetic chirality of 2d arrays under normal illumination. Opt Lett 2017;42:4075–8.

- [22]↑
Kuwata-Gonokami M, Saito N, Ino Y, et al. Giant optical activity in quasi-two-dimensional planar nanostructures. Phys Rev Lett 2005;95:227401.

- [23]↑
Wu C, Arju N, Kelp G, et al. Spectrally selective chiral silicon metasurfaces based on infrared fano resonances. Nat Commun 2014;5:3892.

- [24]↑
Ye W, Yuan X, Guo C, Zhang J, Yang B, Zhang S. Large chiroptical effects in planar chiral metamaterials. Phys Rev Appl 2017;7:054003.

- [25]↑
Zhu AY, Chen WT, Zaidi A, et al. Giant intrinsic chiro-optical activity in planar dielectric nanostructures. Light Sci Appl 2018;7:17158.

- [26]↑
Zhao R, Sain B, Wei Q, et al. Multichannel vectorial holographic display and encryption. Light Sci Appl 2018;7:95.

- [27]↑
Arbabi A, Horie Y, Bagheri M, Faraon A. Dielectric metasurfaces for complete control of phase and polarization with subwavelength spatial resolution and high transmission. Nat Nanotech 2015;10:937–43.

- [28]↑
Karakasoglu I, Xiao M, Fan S. Polarization control with dielectric helix metasurfaces and arrays. Opt Exp 2018;26:21664–74.

- [29]↑
Zhu AY, Chen W-T, Khorasaninejad M, et al. Ultra-compact visible chiral spectrometer with meta-lenses. APL Photon 2017;2:036103.

- [30]↑
Balthasar Mueller JP, Rubin NA, Devlin RC, Groever B, Capasso F. Metasurface polarization optics: independent phase control of arbitrary orthogonal states of polarization. Phys Rev Lett 2017;118:113901.

- [33]↑
Kamali SM, Arbabi E, Arbabi A, Faraon A. A review of dielectric optical metasurfaces for wavefront control. Nanophotonics 2018;7:1041.

- [34]↑
Limonov MF, Rybin MV, Poddubny AN, Kivshar YS. Fano resonances in photonics. Nat Photon 2017;11:543–54.

- [35]↑
Fu YH, Kuznetsov AI, Miroshnichenko AE, Yu YF, Luk’yanchuk B. Directional visible light scattering by silicon nanoparticles. Nat Commun 2013;4:1527.

- [36]↑
Khanikaev AB, Wu C, Shvets G. Fano-resonant metamaterials and their applications. Nanophotonics 2013;2:247.

- [37]↑
Person S, Jain M, Lapin Z, Sáenz JJ, Wicks G, Novotny L. Demonstration of zero optical backscattering from single nanoparticles. Nano Lett 2013;13:1806–9.

- [38]↑
Kim J, Li Y, Miskiewicz MN, Oh C, Kudenov MW, Escuti MJ. Fabrication of ideal geometric-phase holograms with arbitrary wavefronts. Optica 2015;2:958–64.

- [39]↑
Hsu CW, Zhen B, Stone AD, Joannopoulos JD, Soljačić M. Bound states in the continuum. Nat Rev Mater 2016;1:16048.

- [40]↑
Koshelev K, Favraud G, Bogdanov A, Kivshar Y, Fratalocchi A. Nonradiating photonics with resonant dielectric nanostructures. Nanophotonics 2019;8:725.

- [41]↑
Collin S. Nanostructure arrays in free-space: optical properties and applications. Rep Prog Phys 2014;77:126402.

- [42]↑
Menzel C, Rockstuhl C, Lederer F. Advanced jones calculus for the classification of periodic metamaterials. Phys Rev A 2010;82:053811.

- [43]↑
Molesky S, Lin Z, Piggott AY, Jin W, Vucković J, Rodriguez AW. Inverse design in nanophotonics. Nat Photon 2018;12: 659–70.

- [44]↑
Kalish AN, Komarov RS, Kozhaev MA, et al. Magnetoplasmonic quasicrystals: an approach for multiband magneto-optical response. Optica 2018;5:617–23.

- [45]↑
Arteaga O, Maoz BM, Nichols S, Markovich G, Kahr B. Complete polarimetry on the asymmetric transmission through subwavelength hole arrays. Opt Express 2014;22:13719–32.

- [46]↑
Vázquez-Lozano JE, Martínez A. Optical chirality in dispersive and lossy media. Phys Rev Lett 2018;121:043901.

- [47]↑
Alpeggiani F, Bliokh KY, Nori F, Kuipers L. Electromagnetic helicity in complex media. Phys Rev Lett 2018;120:243605.

- [48]↑
Tang Y, Cohen AE. Optical chirality and its interaction with matter. Phys Rev Lett 2010;104:163901.

- [49]↑
Tang Y, Cohen AE. Enhanced enantioselectivity in excitation of chiral molecules by superchiral light. Science 2011;332:333–6.

- [50]↑
Mun J, Rho J. Importance of higher-order multipole transitions on chiral nearfield interactions. Nanophotonics 2019;8:941.

- [51]↑
Ge L, Feng L. Contrasting eigenvalue and singular-value spectra for lasing and antilasing in a 𝒫𝒯-symmetric periodic structure. Phys Rev A 2017;95:013813.

- [52]↑
Bai B, Svirko Y, Turunen J, Vallius T. Optical activity in planar chiral metamaterials: theoretical study. Phys Rev A 2007;76:023811.

- [53]↑
Bai B, Konishi K, Meng X, et al. Mechanism of the large polarization rotation effect in the all-dielectric artificially chiral nanogratings. Opt Exp 2009;17:688–96.

- [54]↑
Ossikovski R. Interpretation of nondepolarizing Mueller matrices based on singular-value decomposition. J Opt Soc Am A 2008;25:473–82.

- [55]↑
Arteaga O, Canillas A. Pseudopolar decomposition of the Jones and Mueller-Jones exponential polarization matrices. J Opt Soc Am A 2009;26:783–93.

- [56]↑
Ossikovski R, Kuntman MA, Arteaga O. Anisotropic integral decomposition of depolarizing Mueller matrices. OSA Continuum 2019;2:1900–7.

## Footnotes

## Supplementary Material

The online version of this article offers supplementary material (https://doi.org/10.1515/nanoph-2019-0321).

## Footnotes

^{1}

We recall that we are limiting ourselves to the regime where the wavelength is larger than the pattern periodicity, i.e. to the *metasurface regime*, where no diffraction takes place.