Abstract
Topological phase transitions in condensed matter systems give rise to exotic states of matter such as topological insulators, superconductors, and superfluids. Photonic topological systems open a whole new realm of research and technological opportunities, exhibiting a number of important distinctions from their condensed matter counterparts. Photonic modes can leak into free space, which makes it possible to probe topological photonic phases by spectroscopic means via Fano resonances. Based on this idea, we develop a technique to retrieve the topological properties of alldielectric metasurfaces from the measured farfield scattering characteristics. Collected angleresolved spectra provide the momentumdependent frequencies and lifetimes of the photonic modes that enable the retrieval of the effective Hamiltonian and extraction of the topological invariant. Our results demonstrate how the topological states of open nonHermitian systems can be explored via farfield measurements, thus paving a way to the design of metasurfaces with unique scattering characteristics controlled via topological effects.
Introduction
Topological phase transitions in twodimensional (2D) condensed matter systems have attracted an enormous interest^{1,2,3,4} crowned by the Nobel Prize in Physics in 2016. With the advent of photonic topological insulators^{5,6,7,8,9,10,11,12} the research domain expanded to include topological transitions for light^{13,14}. However, most of the topological photonic systems considered until today have largely ignored the fact that photonic modes can couple to the continuum of freespace modes, and the exploration of the topological modes has been limited to the nearfield properties. Only recently have topological phenomena in photonic systems associated with leaky states attracted attention in the context of nonHermitian systems^{15,16}, including Floquet systems^{17,18} and nonradiative modes in the continuum^{19}.
While the leakage of photonic modes to the freespace continuum can be considered as a challenge, as it makes the system nonHermitian, it may also offer an alternative route to explore the topology of photonic bands with farfield measurements, provided that the radiative decay of the topological modes is controllable and nondestructive. This is especially relevant in the context of metasurfaces—ultrathin photonic structures enabling a variety of novel optical devices, including flat lenses and waveplates to control polarization and angular momentum of optical beams^{20,21,22,23,24,25,26}. Endowing metasurfaces with topological properties may open a whole new realm of opportunities in the field of light scattering and wavefront control, as the synthetic gauge field for light have been shown to controllably affect all degrees of freedom associated with electromagnetic radiation^{24,27}. Provided the metasurface is engineered to sustain topological photonic modes above the light line, they may couple to the radiative continuum, giving rise to unique scattering features, such as Fano resonances discriminating light by its polarization, angular momentum, or any other synthetic degree of freedom.
In recent years, it has become a traditional approach to look at the topological properties of photonic systems through the prism of the bulkboundary correspondence principle^{28}, focusing on investigation of their edge characteristics, including chiral and helical edge states. Although this strategy to explore topological properties^{6,9,29,30} is quite appealing due to the distinctive features such as backscatteringimmune propagation, it can divert us from new and potentially beneficial properties of bulk modes stemming from their topological nature.
To explore these ideas, we design and fabricate a metasurface that supports radiative photonic modes in the nearinfrared spectral range and exhibits a topological transition. As the parameters of the metasurface are tuned, the coupling of topological modes to the radiation continuum enables the observation of a topological transition, which is accompanied by the inversion of bright and dark modes, referred to as band crossing. Angleresolved spectroscopy allows for the direct extraction of spectral positions, intensity, and bandwidth of the peaks corresponding to the excitation of bulk states of the metasurface from the transmission and reflection spectra. As demonstrated below, these parameters provide valuable information about the nature of the eigenmodes supported by the structure, enabling the univocal retrieval of topological invariants from farfield measurements.
Results
Topological metasurface
The metasurface consists of cylindrical Si pillars arranged into hexagon clusters with edge length R, placed at the sites of a triangular lattice with period a (Fig. 1a, b). The topological properties of its infinite 2D analog possessing C_{6v} symmetry have been the subject of several recent studies^{31,32,33,34,35,36,37}. For a/R = 3, this system represents a conventional honeycomb lattice (with the unit cell containing two cylinders), which exhibits Dirac cones centered at K and K′ points of the Brillouin zone. However, if the lattice symmetry is reduced by clustering six neighboring pillars so that \(a{\mathrm{/}}R\not = 3\) (distorted lattice), the size of the unit cell increases leading to the reshaping of the first Brillouin zone as shown in Fig. 1c. As a consequence, the Dirac points appear in the vicinity of the Γ point due to the bands folding in the distorted lattice as shown in Fig. 1d. Additionally, the interaction between the valleys of the former honeycomb lattice caused by such symmetry reduction leads to the opening of photonic bandgaps (Fig. 1e, g). This interaction can be viewed as synthetic spin–orbit coupling between two pseudospin states corresponding to the former valley degree of freedom of the unperturbed honeycomb lattice. Previous studies suggested that the shrunken structure with a/R > 3 is topologically trivial, whereas the expanded system with a/R < 3 is topologically nontrivial^{31} (see also Supplementary Note 1), which, as shown below, remains true for the case of the open system. The important distinction, however, is that in the case of the open metasurface, the folded modes appear within the light cone, and, therefore, are radiatively coupled to the continuum of free space, which enables their farfield characterization^{38} by spectroscopic methods (Fig. 1f, h).
To perform spectroscopy measurements of metasurfaces with two topologically distinct geometries, we fabricated two sets of samples (Fig. 1a, b) with the same lattice constant a = 750 nm, radius of the pillars r = 75 nm, and height of the pillars h = 1000 nm. The sizes of the clusters R shown in Fig. 1a, b were chosen to be a/R = 3.15 and a/R = 2.85 for shrunken and expanded structures, respectively.
Effective Hamiltonian
Although the topological properties of the perturbed honeycomb lattice have been previously explored under the tight binding approximation^{31}, here we take a different approach, fully considering the electromagnetic nature of the system. Correspondingly, the effective Hamiltonian near the Γ point of the Brillouin zone is obtained directly from Maxwell’s equations by the plane wave expansion method^{8,12} (Supplementary Note 2). To describe the properties of photonic bands of bulk modes in the 2D photonic crystal in the vicinity of the Γ point, we construct an effective 4 × 4 Hamiltonian. Focusing on vertical dipolar modes with E_{ z } component of the electric field directed along the axis of Si rods, we find the results in compliance with the tight binding model (Supplementary Note 1 and refs.^{31,32}). The derived effective 4 × 4 Hamiltonian assumes the following form^{31}:
where the 2 × 2 matrices \(\hat H_ \pm\) and \(\hat K\) are
μ(k) = μ + βk^{2}, μ and β are the mass term and band parabolicity, respectively. This form of the Hamiltonian corresponds to the basis choice
where, \(\left {p_ \pm } \right\rangle = p_x \pm i{\kern 1pt} p_y\) are the circularly polarized dipolar modes and \(\left {d_ \pm } \right\rangle = d_{x^2  y^2} \pm i{\kern 1pt} d_{xy}\) are the circularly polarized quadrupolar modes, which both originate from modes supported by an isolated metamolecule (the hexamer of 6 dielectric rods).
The four eigenstates of the Hamiltonian Eq. (1) exhibit a pairwise degeneracy at the Γ point k = 0 (Supplementary Note 1 and 2), with eigenvalues equal to μ (dipolar modes) and −μ (quadrupolar modes). Numerically computed field distribution of the modes (shown in Fig. 1i, j) demonstrate inplane and outofplane variation, respectively, and confirm their dipolar and quadrupolar form and strong confinement to the pillars. Transition from shrunken to expanded design is accompanied by the band inversion for dipolar and quadrupolar eigenstates, and by the reversal of the sign of the mass term μ. In addition, the degeneracy between left and righthanded circularly polarized modes is removed for nonzero values of k due to the term being proportional to α. However, for the topological invariant calculation, this term is inessential and can be dropped^{31}. As a result of this approximation, the Hamiltonian splits into two independent 2 × 2 blocks, and the effective pseudospin can be introduced. A straightforward calculation^{39} of spin Chern number (Supplementary Note 1) yields
thus, to evaluate the topological invariant one has to extract the mass term μ and the parabolicity of the bands β from the experimentally measured spectra.
Retrieval of topological state via farfield spectroscopy
The modes of interest are located above the light cone, and, for the case of a metasurface of a finite thickness, they couple to the radiation continuum of an open background medium^{38}. As a result, the modes acquire a finite lifetime and the corresponding eigenfrequencies become complexvalued. Thus, the silicon pillars constituting the metasurface effectively act as laterally coupled cavities, and the leaky modes of the structure can be excited by the polarization currents induced at the interface of the metasurface by the incident fields. However, taking into consideration the symmetry of the modes, normally incident light can only couple to the dipolar modes, which thus represent superradiant modes, whereas coupling to the quadrupolar modes is suppressed because of the symmetry mismatch with the incident field, and they represent subradiant modes. The quadrupolar modes can still be excited indirectly through coupling to the dipolar modes at oblique incidence, because of the hybridization of dipolar and quadrupolar modes away from the Γ point. This mechanism is referred to as an extrinsic coupling to dark modes in Fanoresonant systems and is caused by the symmetry reduction due to the finite value of the inplane wave vector components.
Of special interest are the effects of the topological transition on the radiative coupling of the modes to the continuum, which should follow the band inversion described above, giving rise to switching of the bright (dipolar) and dark (quadrupolar) modes, and, thus, enabling control over the farfield scattering properties of the metasurface.
To support our claims, the complex photonic band structures for both the shrunken and expanded metasurfaces were calculated in Comsol Multiphysics with radiative decay fully considered. The results are shown in Fig. 1e, g. Indeed, the dipolar and quadrupolar modes at the Γ point appear to be bright and dark, respectively, as reflected by the radiative quality factors (essentially, the lifetimes of the modes). Even at oblique incidence, when modes are hybridized, the respective bands may have radiative quality factors that differ by orders of magnitude. The numerically calculated extinction spectra at oblique incidence are shown in Fig. 1f, h, and confirm excitation of the two eigenmodes, which manifest as two peaks. The predominantly dipolar and quadrupolar modes can be easily discriminated by their distinct bandwidth and the amplitude of the corresponding peaks. This allows one to establish the symmetry of the mode without the need to look into the nearfield pattern and to extract the mass term μ needed for evaluation of the topological invariant.
To formalize our description in the context of the topological properties of leaky modes, we develop a unified approach utilizing the temporal coupled mode theory (CMT) along with the electromagnetic effective Hamiltonian description. We write the CMT equations for the amplitudes of the righthanded and lefthanded circularly polarized leaky modes coupled to the external source in the following blockdiagonal form:
where \(\left {\psi _ \pm } \right\rangle = \left( {\left {p_ \pm } \right\rangle ,\left {d_ \pm } \right\rangle } \right)^T\) is the mode composed of p (dipole) and d (quadrupole) eigenmodes of the system. The first term of Eq. (6) describes the evolution of the coupled modes in a closed system, \(\chi\) describes the coupling strength of the system to the external field \(E_{{\mathrm{in}}} = E_0{\mathrm{/}}\sqrt 2\) (E_{0} is the field amplitude), γ_{0} represents nonradiative losses in the structure, and \(\gamma _{\mathrm{r}} = \chi ^2{\mathrm{/}}2\) describes the radiative losses. The use of this form of dynamic equation can be also justified by applying electromagnetic perturbation theory for an open system as shown in Supplementary Note 4.
Once Eq. (6) is solved, one can obtain expressions for the transmission (and reflection) coefficients:
where ψ_{±}(p) is the first (dipolar) component of the twocomponent wave function \(\left {\psi _ \pm } \right\rangle\). Since the results are the same for lefthanded and righthanded circular polarizations, from now on we omit the ± subscript. Plugging the Hamiltonian into Eq. (7), we derive the expression for the extinction in the form
Equation (8) suggests in particular that two peaks in the \(\tilde R\) spectrum can be observed due to excitation of dipolar and quadrupolar eigenmodes. Thus, the frequency of the most pronounced peak at incidence close to normal is around f_{0} + μ, whereas the frequency of the less intense peak is found around f_{0} − μ, where f_{0} is the midgap frequency.
Overall, there are six parameters in the effective model that control the entire dispersion and topological properties of the farfield response: μ, \(\gamma _{\mathrm{r}} = \chi ^2{\mathrm{/}}2\), γ_{0}, v, and β (see Eqs. (2) and (6)). The numerical values of parameters comprising the effective Hamiltonian can be obtained by fitting the experimental data. Next, by using the developed technique, we analyze the experimental data for both shrunken (a/R = 3.15) and expanded (a/R = 2.85) structures and extract their topological characteristics.
To confirm our theoretical predictions, a set of shrunk and expanded silicon samples on top of sapphire (Al_{2}O_{3}) substrates was fabricated and tested experimentally. The structures were illuminated using the halogen source Ocean Optics HL2000LL. The transmission was measured for both structures by the spectrometer Ocean Optics NIR Quest NQ 512–2.2 in the range of the wavelengths from 896 to 2142 nm corresponding to the frequency range (1.40 ÷ 3.35) × 10^{14} Hz. The measured transmittance was normalized to the one of a sapphire substrate. Figure 2 shows the color map of the quantity \(\tilde R = 1  T{\mathrm{/}}T_0\) as a function of the wavelength and the incidence angle for both the shrunken and expanded structures.
In the wavelength range from 1600 to 1900 nm, two peaks are clearly observed. The detailed spectroscopy results for the whole studied spectral range are provided in Supplementary Note 3. In agreement with our numerical calculations, the experimental data clearly show that the most intense and broad peak for the shrunken structure is the lowfrequency one, whereas for the expanded structure the highfrequency peak is more pronounced. Thus, the shrunken structure is characterized by a negative effective mass μ, while the expanded structure is described by a positive effective mass μ.
Next, we applied the above analytical model to fit the measured data with Eqs. (6) and (8), as detailed in the Supplementary Note 3. The results of the simplest fitting algorithm for the two particular angles of incidence, 0° and 16°, are illustrated in Fig. 3 for the cases of both shrunken and expanded structures. The extracted parameters of the effective Hamiltonian are listed in Table 1. Importantly, these numbers appear to be nearly independent of the fitting algorithm (see details in Supplementary Note 3). From the measured transmittance spectra, we recover again that the shrunken structure is topologically trivial, since the effective mass μ and parabolicity parameter β have the same sign, yielding zero spin Chern number. On the contrary, the expanded structure appears to be topologically nontrivial, due to the opposite signs of μ and β, yielding a spin Chern number of 1. It is important to mention that for nonHermitian systems (including the one considered here) the definition of topological invariant has been a subject of active discussion^{15,16,17,18,39}. In Supplementary Note 7, following ref. ^{39}, we show that the topological invariant remains unaffected by the presence of loss for the case of the metasurface studied here.
Interestingly, the topological properties have been directly extracted previously from the bulk spectrum for 1D systems only^{40}. A different approach based on the edge properties was taken by Mittal et al.^{41}, who introduced a unit quantum flux at the edge and observed a shift in the position of the edge spectrum. Our results prove for the first time that it is possible to retrieve the topological properties of an open 2D system via farfield measurements in a robust and direct way.
Farfield and nearfield observation of edge states
In order to confirm that the structure under study indeed supports topological edge modes, we have performed additional experimental studies of the domain wall formed between the two topologically distinct structures (shrunken and expanded) and further visualized edge modes with the farfield and nearfield experimental techniques. The farfield experimental studies have been performed on a sample consisting of an array of domain walls formed by repetitive stitching of 12 unit cells of shrunken and expanded domains along a zigzag cut of the lattice. The schematic of the experimental setup is shown in Fig. 4a with an SEM image of the sample near the domain wall in Fig. 4b. The results of numerical calculations performed within the use of the CMT are presented in Fig. 4c and reveal the presence of the edge states within the bandgap of bulk modes. The examination of numerically computed field profiles of the edge modes overlaid with the SEM image confirms their localization at the domain wall as shown in falsecolor plot in Fig. 4b. The experimental farfield measurements of sample’s reflectance are presented in Fig. 4b and confirm the presence of edge states within the topological bandgap. Interestingly, despite the presence of the additional scattering channels due to diffraction, the modes have relatively narrow bandwidth and they can be easily distinguished from the bulk modes by their distinctive linear (1D Diraclike) dispersion at small tangential wavenumbers. Thus, the leakage of the edge states appears to be strong enough so that they can be excited and measured in farfield region, yet, the nonHermiticity it introduces is not strong enough to destroy the topological state of the system. Thus, the presence of the edge states in our system serves as another piece of evidence of the topological properties of the structure despite its open nature.
Additional proofofconcept experimental nearfield studies of the edge modes were conducted in the microwave spectral range. Note that, due to the dielectric nature of the structure and scalability of Maxwell’s equations, the obtained results can be directly transferred to the structures operating in the optical domain.
We considered a metasurface consisting of dielectric rods arranged into the lattice with the permittivity of the rods ε = 10, lattice period a = 37.5 mm, radius and height of dielectric rods are r = 1.5 mm and h = 10 mm, respectively. The cluster size R (Fig. 1a, b) was chosen such that a/R = 3.33 and a/R = 2.73 for shrunken and expanded structures, respectively. For these parameters, the bandgap is characterized by the midgap frequency of 6.79 GHz with a total bandwidth of 310 MHz.
The experimental sample was constructed by embedding commercially available dielectric powder (Eccostock HiK with the permittivity ε = 10 and loss tangent tan δ = 0.0007) (Ceramic powder datasheet https://www.eccosorb.com/productseccostockhikpowder.htm) into the machineprocessed substrate made from Styrofoam material (permittivity of 1.15 and negligible losses) in which the holes were drilled (Fig. 5a). All the measurements were performed in an anechoic chamber using an automatic mechanical nearfield scanning setup and an electric field probe connected to the receiving port of the vector network analyzer (Agilent E8362C). The probe was oriented normally with respect to the interface of the structure. The nearfield was scanned at 2 mm distance from the back interface of the metasurface to avoid a direct contact between the probe and the sample. To excite the edge state, the subwavelength dipole source whose polarization could be adjusted from circular to linear was placed next to the domain wall in the center of the expanded unit cell. Figure 5b shows the results of full wave simulations of the field distribution of the edge state. The corresponding experimental map presented in Fig. 5c clearly shows that the mode excited inside the bandgap at frequency f = 6.84 GHz is indeed localized to the domain wall and represents the edge state. To confirm the topological properties of the edge states, we performed nearfield excitation both by linearly and circularly polarized dipole fields. As can be seen from Fig. 5b, c, the linearly polarized dipole leads to bidirectional excitation of the edge state. In contrast, the circularly polarization excitation leads to oneway propagation of the edge state evidencing the locking of the polarization and the propagation direction, thus, confirming the topological nature of the corresponding electromagnetic mode^{42,43}.
These results prove independently that the metasurface considered here exhibits a topological transition when its design is changed from shrunken to expanded structure, thus validating the introduced farfield measurement technique for the detection of a topological transition in open photonic systems such as metasurfaces.
Discussion
In this paper, we have introduced the concept of topological metasurfaces and demonstrated that their scattering characteristics, radiative quality factors of the modes in particular, can be controlled by synthetic gauge fields. The developed formalism allows us to relate topological properties and the farfield scattering, thus enabling the extraction of the effective Hamiltonian and the topological invariant (pseudospin Chern number) of the metasurface from the measured transmission/reflection spectra. Coupling of photonic modes of the structure to the freespace modes allows for the direct probing and visualization of the topological phase transitions in the open photonic system. We believe that our results open a new avenue for designing metasurfaces with desirable scattering properties controllable via their topological characteristics. Provided the properties of metasurfaces are designed to stem from edge states, their response may exhibit topological robustness to continuous modifications and disorder, which can be of immense value for practical applications.
Methods
Modeling and calculations
To derive the effective Hamiltonian and topological properties of the photonic crystal, we use the tight binding approach (Supplementary Note 1), as well as the plane wave expansion method (Supplementary Note 2). To calculate the farfield scattering properties of the structure, we employ the CMT (Supplementary Note 3). Perturbative analysis of the radiative losses for the metasurface is performed by the guided mode expansion method (Supplementary Note 4). The calculation of scattering from the topological metagrating based on the CMT is detailed in Supplementary Note 6. The role of nonHermiticity and topological invariants in the presence of radiative losses are discussed in Supplementary Note 7.
Sample fabrication
A scheme of the fabrication procedure is outlined in Supplementary Note 5 and Supplementary Fig. 6.
Data availability
The data that support the findings of this study are available from the corresponding author upon request.
References
 1.
Kosterlitz, J. M. & Thouless, D. J. Ordering, metastability and phase transitions in twodimensional systems. J. Phys. C. 6, 1181–1203 (1973).
 2.
Bernevig, B. A., Hughes, T. L. & Zhang, S.C. Quantum spin Hall effect and topological phase transition in HgTe quantum wells. Science 314, 1757–1761 (2006).
 3.
Lutchyn, R. M., Sau, J. D. & Sarma, S. D. Majorana fermions and a topological phase transition in semiconductorsuperconductor heterostructures. Phys. Rev. Lett. 105, 077001 (2010).
 4.
Dagvadorj, G. et al. Nonequilibrium phase transition in a twodimensional driven open quantum system. Phys. Rev. X 5, 041028 (2015).
 5.
Raghu, S. & Haldane, F. D. M. Analogs of quantumHalleffect edge states in photonic crystals. Phys. Rev. A. 78, 033834 (2008).
 6.
Wang, Z., Chong, Y., Joannopoulos, J. D. & Soljačić, M. Observation of unidirectional backscatteringimmune topological electromagnetic states. Nature 461, 772–775 (2009).
 7.
Hafezi, M., Demler, E. A., Lukin, M. D. & Taylor, J. M. Robust optical delay lines with topological protection. Nat. Phys. 7, 907–912 (2011).
 8.
Khanikaev, A. B. et al. Photonic topological insulators. Nat. Mater. 12, 233 (2013).
 9.
Hafezi, M., Mittal, S., Fan, J., Migdall, A. & Taylor, J. M. Imaging topological edge states in silicon photonics. Nat. Photonics 7, 1001–1005 (2013).
 10.
Rechtsman, M. C. et al. Photonic Floquet topological insulators. Nature 496, 196–200 (2013).
 11.
Ozawa, T., Price, H. M., Goldman, N., Zilberberg, O. & Carusotto, I. Synthetic dimensions in integrated photonics: from optical isolation to fourdimensional quantum Hall physics. Phys. Rev. A. 93, 043827 (2016).
 12.
Slobozhanyuk, A. et al. Threedimensional alldielectric photonic topological insulator. Nat. Photonics 11, 130–136 (2016).
 13.
Lu, L., Joannopoulos, J. D. & Soljačić, M. Topological photonics. Nat. Photonics 8, 821–829 (2014).
 14.
Lu, L., Joannopoulos, J. D. & Soljačić, M. Topological states in photonic systems. Nat. Phys. 12, 626 (2016).
 15.
Esaki, K., Sato, M., Hasebe, K. & Kohmoto, M. Edge states and topological phases in nonHermitian systems. Phys. Rev. B 84, 205128 (2011).
 16.
Leykam, D., Bliokh, K. Y., Huang, C., Chong, Y. D. & Nori, F. Edge modes, degeneracies, and topological numbers in nonhermitian systems. Phys. Rev. Lett. 118, 040401 (2017).
 17.
Zeuner, J. M. et al. Observation of a topological transition in the bulk of a nonhermitian system. Phys. Rev. Lett. 115, 040402 (2015).
 18.
Weimann, S. et al. Topologically protected bound states in photonic parity–timesymmetric crystals. Nat. Mater. 16, 433–438 (2016).
 19.
Hsu, C. W., Zhen, B., Stone, A. D., Joannopoulos, J. D. & Soljačić, M. Bound states in the continuum. Nat. Rev. Mater. 1, 16048 (2016).
 20.
Capasso, N. Y. F. Flat optics with designer metasurfaces. Nat. Mater. 13, 139–150 (2014).
 21.
Kildishev, A. V., Boltasseva, A. & Shalaev, V. M. Planar photonics with metasurfaces. Science 339, 1232009 (2013).
 22.
Minovich, A. E. et al. Functional and nonlinear optical metasurfaces. Laser & Photon. Rev. 9, 195–213 (2015).
 23.
Bliokh, K. Y., RodrguezFortuño, F. J., Nori, F. & Zayats, A. V. Spin–orbit interactions of light. Nat. Photonics 9, 796–808 (2015).
 24.
Bliokh, K. Y. & Nori, F. Transverse and longitudinal angular momenta of light. Phys. Rep. 592, 1–38 (2015).
 25.
Glybovski, S. B., Tretyakov, S. A., Belov, P. A., Kivshar, Y. S. & Simovski, C. R. Metasurfaces: from microwaves to visible. Phys. Rep. 634, 1–72 (2016).
 26.
Xiao, S. et al. Spindependent optics with metasurfaces. Nanophotonics 6, 215–234 (2017). https://doi.org/10.1515/nanoph20160121
 27.
Bliokh, K. Y., Smirnova, D. & Nori, F. Quantum spin Hall effect of light. Science 348, 1448–1451 (2015).
 28.
Bernevig, A. B. & Hughes, T. L. Topological insulators and topological superconductors (Princeton University Press, New Jersey, 2013).
 29.
Ma, T., Khanikaev, A. B., Mousavi, S. H. & Shvets, G. Guiding electromagnetic waves around sharp corners: topologically protected photonic transport in metawaveguides. Phys. Rev. Lett. 114, 127401 (2015).
 30.
Slobozhanyuk, A. P. et al. Experimental demonstration of topological effects in bianisotropic metamaterials. Sci. Rep. 6, 22270 (2016).
 31.
Wu, L.H. & Hu, X. Scheme for achieving a topological photonic crystal by using dielectric material. Phys. Rev. Lett. 114, 23901 (2015).
 32.
Barik, S., Miyake, H., DeGottardi, W., Waks, E. & Hafezi, M. Twodimensionally confined topological edge states in photonic crystals. New. J. Phys. 18, 113013 (2016).
 33.
Dong, J.W., Chen, X.D., Zhu, H., Wang, Y. & Zhang, X. Valley photonic crystals for control of spin and topology. Nat. Mater. 16, 298–302 (2016).
 34.
Ma, T. & Shvets, G. AllSi valleyHall photonic topological insulator. New. J. Phys. 18, 025012 (2016).
 35.
Noh, J. et al. Topological protection of photonic midgap cavity modes. Preprint at https://arxiv.org/abs/1611.02373 (2016).
 36.
Yves, S. et al. Crystalline metamaterials for topological properties at subwavelength scales. Nat. Commun. 8, 16023 (2017).
 37.
Yang, Y. et al. Visualization of unidirectional optical waveguide using topological photonic crystals made of dielectric material. Preprint at https://arxiv.org/abs/1610.07780 (2016).
 38.
Reagan, E. C. et al. Direct imaging of isofrequency contours in photonic structures. Sci. Adv. 2, e1601591 (2016).
 39.
Shen, S.Q. Topological Insulators. Dirac Equation in Condensed Matters (Springer, Berlin, 2012).
 40.
Atala, M. et al. Direct measurement of the zak phase in topological bloch bands. Nat. Phys. 9, 795 (2013).
 41.
Mittal, S., Ganeshan, S., Fan, J., Vaezi, A. & Hafezi, M. Measurement of topological invariants in a 2d photonic system. Nat. Photonics 10, 180–183 (2016).
 42.
Yang, Y. et al. Visualization of unidirectional optical waveguide using topological photonic crystals made of dielectric material. Preprint at https://arxiv.org/abs/1610.07780 (2017).
 43.
Barik, S. et al. A topological quantum optics interface. Preprint at https://arxiv.org/abs/1711.00478 (2017).
Acknowledgements
The work was supported by the National Science Foundation grants CMMI1537294 and EFRI1641069. Research carried out in part at the Center for Functional Nanomaterials, Brookhaven National Laboratory, which is supported by the U.S. Department of Energy, Office of Basic Energy Sciences, under Contract No. DESC0012704. The theoretical models were developed with the support of the Russian Science Foundation (Grant No. 161910538).
Author information
Affiliations
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Gorlach, M.A., Ni, X., Smirnova, D.A. et al. Farfield probing of leaky topological states in alldielectric metasurfaces. Nat Commun 9, 909 (2018). https://doi.org/10.1038/s41467018033309
Received:
Accepted:
Published:
Further reading

Tailoring topological edge states with photonic crystal nanobeam cavities
Scientific Reports (2021)

Experimental observation of topological Z2 excitonpolaritons in transition metal dichalcogenide monolayers
Nature Communications (2021)

Probing bulk topological invariants using leaky photonic lattices
Nature Physics (2021)

A highperformance topological bulk laser based on bandinversioninduced reflection
Nature Nanotechnology (2020)

Roomtemperature lasing from nanophotonic topological cavities
Light: Science & Applications (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.