## Exploring photonic topology

Scattering topological effects are being explored in a variety of electronic and optical materials systems owing to their robustness against defects (see the Perspective by Özdemir). Yang *et al.* designed and fabricated an ideal optical analog of a three-dimensional Weyl system. Angular transmission measurements revealed four Weyl points at the same energy, as well as the signature helicoidal arcs associated with such an exotic topological system. Zhou *et al.* theoretically proposed and experimentally demonstrated the formation of a topologically protected bulk Fermi arc. They attributed the formation of the arc to the topological nature of paired exceptional points (points at which gain and loss in the system are matched). Photonic crystals may provide a powerful platform for studying exotic properties of topological electronic systems and may also be used to develop optical devices that exploit topological properties of light-matter interactions.

## Abstract

Weyl points are the crossings of linearly dispersing energy bands of three-dimensional crystals, providing the opportunity to explore a variety of intriguing phenomena such as topologically protected surface states and chiral anomalies. However, the lack of an ideal Weyl system in which the Weyl points all exist at the same energy and are separated from any other bands poses a serious limitation to the further development of Weyl physics and potential applications. By experimentally characterizing a microwave photonic crystal of saddle-shaped metallic coils, we observed ideal Weyl points that are related to each other through symmetry operations. Topological surface states exhibiting helicoidal structure have also been demonstrated. Our system provides a photonic platform for exploring ideal Weyl systems and developing possible topological devices.

Topology is the mathematics of conserved properties under continuous deformations, and the recent study of band topologies is yielding a suite of fascinating interface transport phenomena that include one-way propagation of energy and previously unknown relativistic behavior. The two-dimensional honeycomb lattice is the most studied in the exploration of topological phenomena. Made famous by graphene (*1*), the energy-momentum dispersion in a honeycomb system is linear, and the crossings of bands in energy-momentum space are known as Dirac points. The transport of the quasiparticles around these points is massless, and this remarkable transport behavior is associated with the “hidden” symmetry associated with its two identical sublattices. Weyl points are the characteristic of an analogous phenomenon when the lattice is extended to three dimensions (*2*–*5*). In electronic systems, materials exhibiting Weyl points are known as Weyl semimetals, and Weyl fermion is the solution to the massless Dirac equation. Each Weyl point can be assigned an integer “charge” based on its chirality, known as the Chern number, and much like magnetic monopoles, Weyl points are only ever found in pairs of opposite charge. Just like Dirac points, these Weyl points also exist in photonic systems, but unlike Dirac points, they can only exist once either (or both) time-reversal or space-inversion symmetry of the crystal is broken. To date, Weyl points of various forms have been proposed and realized in several boson or fermion systems (*2*–*4*, *6*–*13*). Among them, the presence of surface state arcs as one of the fingerprints of Weyl systems has been observed.

However, demonstration of more fundamental topological features of Weyl points—such as the helicoidal dispersion, which yields the open Fermi arcs of topological surface states (*14*)—has been hindered by the complicated configuration of energy bands at the Weyl energy. Moreover, some realistic and innovative device applications critically depend on a simple embodiment of Weyl systems (*5*). Thus, an ideal Weyl system (*15*–*17*) has attracted much attention because in such systems, all Weyl nodes are symmetry-related, residing at the same energy with a large momentum separation and devoid of nontopological bands in a sufficiently large energy interval.

Although Weyl degeneracies can be readily found by breaking either time-reversal or inversion symmetry (*5*), or both, the experimental realization of a truly ideal Weyl system has not yet been reported. Here, we explore the microwave response of a three-dimensional photonic crystal composed of metallic inclusions (termed a “meta-crystal”) in order to realize an ideal Weyl system protected by D_{2d} point symmetry. Our meta-crystal exhibits four Weyl points at the same energy, the minimum number allowed in the presence of time-reversal symmetry. By placing an excitation point-source on one surface of the crystal, and scanning the near fields on the opposite surface, we observed the intriguing helicoidal structure of topological surface states: a physical representation of Riemann surface generated by a multivalued function (*14*).

Our meta-crystal design offers an ideal platform for the investigation of various unconventional physics in Weyl systems. The symmetry of the studied meta-crystal belongs to the simple tetragonal lattice with symmorphic space group (no. 115). The basis comprises a saddle-shaped connective metallic coil (Fig. 1, A and B) that possesses D_{2d} ( in Hermann-Mauguin notation) point group symmetry. The system has no spatial inversion. These metallic elements support localized electromagnetic resonances with current distributions that can be expanded into multipolar modes (*18*). In an effective medium model (supplementary materials section 4 and fig. S1) (*19*), these resonances collectively exhibit a bi-anisotropic effect, leading to a directionally dependent chirality response (*20*). Here, the unavoidable crossings between the longitudinal mode (LM) with negative dispersion and the transverse electric mode (TM) with positive dispersion along Γ–M result in the formation of a type-I Weyl point (Fig. 1C and fig. S2) (*19*, *21*). Analysis via the irreducible representation of the point group shows that these two modes belong to two different classes, with eigenvalues ±1 of C_{2} rotation along Γ–M (Fig. 1D), where level repulsion is forbidden (supplementary materials section 5 and fig. S3) (*19*). The other three Weyl points are obtained after application of the D_{2d} symmetry operation. For instance, three twofold rotation symmetries (C_{2} and 2C_{2}′) combined with time-reversal symmetry guarantee that these four Weyl nodes are located on the Γ–M at the same frequency, at which application of two mirror symmetries (σ* _{x}* and σ

*) reverses the corresponding topological charges. The simulated band structure along high-symmetry lines is shown in Fig. 1E (as defined in Fig. 1D) in the Brillouin zone (BZ), where a pair of Weyl points resides at the same frequency. As such, these Weyl degeneracies appear in a relatively large energy window (~2.1 GHz around the frequency of the Weyl point) (Fig. 1E, blue shaded region) that is also devoid of other bulk bands and hence unequivocally facilitates their experimental identification.*

_{y}The linear band crossings of the bulk states around the Weyl point are confirmed with angle-resolved transmission measurements (*9*). In order to couple energy across the meta-crystal, the momentum of the bulk states must be matched to the in-plane momentum of an incident wave; a sample with special crystal cutting is fabricated (Fig. 2A), in which the crystal orientation forms an angle of 26.57° with one of the cutting boundaries. Compared with the global axis, *xyz*, a local coordinate, *uvw*, is defined. The length (along *u*), width (along *v*), and height (along *w*) of the sample are 300, 100, and 300 mm, respectively. Two angles, θ and ϕ (Fig. 2, A and B), are scanned to obtain the angle-resolved transmission spectrum. With this specific crystal cutting, when ϕ *=* 0°, Weyl points in BZ are projected along the scan wave vector *k*_{p}, which is related to θ, as indicated in Fig. 2C. Obviously, two of the projected Weyl points are located within the light circle (magenta circle) at the Weyl frequency (13.5 GHz). Thus, even a plane wave illuminated directly from air onto the sample can address these two Weyl points. Comparisons between the simulation and experiment results are shown in Fig. 2, D, E, and F. In Fig. 2D, with ϕ *=* 0°, a linear gapless energy dispersion is obtained, and the density of states vanishes at the Weyl frequency because of the absence of other bulk states at the same frequency in an ideal Weyl system. After rotating the sample to ϕ *=* 30° and 60° around the *v* axis, a complete gap is observed as expected.

Another direct manifestation of the topological aspects of a Weyl system is the exotic topological surface states taking the form of arcs connecting the topologically distinct bulk states. Following a closed contour around an end of the arcs, one moves between the lower (valence) and upper (conduction) bands (*14*), which is a direct consequence of the chiral characteristic of Weyl nodes, as schematically shown in Fig. 3A. It is well known that the gapless surface states of Weyl crystals take the form of helicoid Riemann surfaces (*14*), where the bulk Weyl points correspond to the poles and zeros adopting the sign of their respective Chern numbers. Recently, it was shown that topological surface states of double Weyl systems can be analytically expressed, across the entire BZ, as the double-periodic Weierstrass elliptic function (*22*). Because the Weierstrass elliptic function has one second-order pole and one second-order zero, it is not the most fundamental expression of the Weyl surface states. Here, we show that our ideal-Weyl meta-crystal of four Weyl points has surface states whose dispersion is topologically equivalent to the argument of Jacobi elliptic function cn(*z,m*) of two poles and two zeros on the complex plane. cn(*z*,*m*) is a meromorphic function with periods 4*K*(m) and 4*K*(1 – m), where *K* is the complete elliptic integrals of the first kind. For our system, the mapping is given by ω(*k _{x}*,

*k*) ~ arg{cn[(

_{y}*k*–

_{x}*k*)/2 + (

_{y}*k*+

_{x}*k*)

_{y}*i*/2, 1/2]}, as plotted in Fig. 3A.

The helicoidal structure of the surface arcs was probed by using the transmitted near-field scanning configuration, with the excitation source located at the center of the bottom layer of the meta-crystal stack (Fig. 3B, setup “a”), where the detecting probe can raster-scan the top surface so as to map out both the bulk and surface modes. Another configuration (fig. S4B, setup “b”), in which the excitation source is positioned at the edge or corner of the top surface, is also used to identify the surface states. These two setups provide complementary information for the observation of helicoid surface states. In all near-field measurements, we set the scanning step as 1 mm (*a*/3), providing a large surface momentum space in the range of (–3π/*a*, 3π/*a*)^{2} after the Fourier transformation. The helicoid structure of the surface arc was experimentally measured and numerically simulated and is presented as a series of equifrequency contours between 12.6 and 14.0 GHz (Fig. 3, C and E, in experiment, and Fig. 3, D and F, in simulation).

As shown in Fig. 3, C and D, at 13.1 GHz, which is below the Weyl frequency, the Fourier transformation of the experimentally measured field distribution shows the presence of four symmetrically displaced elliptical bulk states with the same size located along the diagonal directions. We clearly observed two surface arcs running across the BZ boundaries and connecting the neighboring bulk states with opposite topological charges. In the vicinity of the air equifrequency contour (circle), there exists a surface ellipse. The surface ellipse joins and reroutes the surface arc at higher frequencies (Fig. 3, E and F). Indeed, the surface ellipse and surface arcs together form the same unified helicoid surface in the dispersion of the surface states.

With increasing frequency, the top surface arc that emerged from the Weyl node with positive and negative topological charge rotates anticlockwise and clockwise, respectively. The observed rotation of the helicoid surface state around a Weyl node can therefore be used to detect the chirality of the Weyl node (*23*). As mentioned above, at lower frequencies each surface arc connects between the bulk states through the BZ boundary, whereas the surface ellipse expands gradually with increasing frequency. Between 13.5 and 13.6 GHz, the surface arc and surface ellipse connect with each other and then transition into a new configuration: a direct surface arc connecting between the bulk states within the BZ, and a surface ellipse centered at its edge. The evolution of the surface arc configuration across the measured frequency range matches topologically with that described by the Jacobi elliptic function shown in Fig. 3A. At the frequency of 14.3 GHz, the surface arcs appear to be linear (fig. S4F) (*19*), leading to nearly diffractionless propagation of the surface wave in the real space (fig. S4C) (*19*). Slightly away from the Weyl frequency, the equi-energy contour of the bulk state consists of four very small spheres enclosing the Weyl points. It is expected that the interference between them results in a chessboard-like interference pattern in real space, which is experimentally confirmed as a spatial frequency filter (fig. S5) (*17*, *19*). In addition, the dimensional reduction from ideal Weyl points to graphene-like dispersion is shown in fig. S6 (supplementary materials section 8) (*19*). We also analyzed that the presence of dielectric loss of our system does not affect the existence of the Weyl points (supplementary materials section 9) (*19*).

The designed ideal Weyl system presented here opens up opportunities for studying intriguing physics and offers a prototype platform for realistic device applications. In photonics, besides the topologically nontrivial surface states supported by Weyl materials, the diverging Berry curvature (*5*) close to Weyl points provides a new degree of freedom in controlling the transport of optical wave packets and may lead to the observation of a gigantic Hall effect for light (*24*). Furthermore, any phenomena related to the conical dispersion of the light cone may be observed around Weyl points, such as diverging and diminishing scattering cross sections (*25*). The vanishing density of states at Weyl frequencies also provides a robust platform for controlling light-matter interaction when emitters are embedded inside photonic Weyl materials.

## Supplementary Materials

www.sciencemag.org/content/359/6379/1013/suppl/DC1

Materials and Methods

Supplementary Text

Figs. S1 to S6

This is an article distributed under the terms of the Science Journals Default License.

## References and Notes

**Acknowledgments:**This work was financially supported by the European Research Council Consolidator Grant (Topological), Horizon 2020 Action Project grant 734578 (D-SPA), and Leverhulme Trust (grant RPG-2012-674). S.Z. acknowledges support from the Royal Society and Wolfson Foundation. B.Y. acknowledges support from China Scholarship Council (grant 201306110041). Q.G. acknowledges the financial support of the National Natural Science Foundation of China (grant 11604216). Y.X. acknowledges support from the National Natural Science Foundation of China (grant 61490713). L.E.B. and A.H. acknowledge financial support from the Engineering and Physical Sciences Research Council of the United Kingdom (grant EP/L015331/1). C.F. was supported by the National Key Research and Development Program of China under grant 2016YFA0302400 and by the National Natural Science Foundation of China (NSFC) under grant 11674370. L.L. was supported by the National Key R&D Program of China under grants 2017YFA0303800 and 2016YFA0302400 and by NSFC under project 11721404. Near-field scanning data were collected by a vector network analyzer controlled with

*xyz*-stage at G31 at the Department of Physics and Astronomy, University of Exeter, UK. All data created during this research are openly available from the University of Exeter’s institutional repository at http://hdl.handle.net/10871/30744.