Arxiv:1707.04503V2 [Physics.Optics] 4 Mar 2018 2

Manipulating Type-I and Type-II Dirac Polaritons in Cavity-Embedded Honeycomb Metasurfaces Charlie-Ray Mann,1 Thomas J. Sturges,1 Guillaume Weick,2 William L. Barnes,1 and Eros Mariani1 1EPSRC Centre for Doctoral Training in Metamaterials (XM2), Department of Physics and Astronomy, University of Exeter, Exeter, EX4 4QL, United Kingdom. 2Universite´ de Strasbourg, CNRS, Institut de Physique et Chimie des Materiaux´ de Strasbourg, UMR 7504, F-67000 Strasbourg, France. Pseudorelativistic Dirac quasiparticles have emerged in a plethora of artificial graphene systems that mimic the underlying honeycomb symmetry of graphene. However, it is notoriously difficult to manipulate their properties without modifying the lattice structure. Here we theoretically investigate polaritons supported by honeycomb metasurfaces and, despite the trivial nature of the resonant elements, we unveil rich Dirac physics stemming from a non-trivial winding in the light-matter interaction. A new kind of type-II Dirac point emerges which simultaneously exists with its conventional type-I counterpart. By modifying only the photonic environment via an enclosing cavity, one can manipulate the location of the type-II Dirac points, leading to distinct polariton phases. This enables one to alter the fundamental properties of the emergent Dirac polaritons while preserving the lattice structure—a unique scenario which has no analog in real or artificial graphene systems so far. Exploiting the photonic environment will thus give rise to unexplored Dirac physics at the subwavelength scale. arXiv:1707.04503v2 [physics.optics] 4 Mar 2018 2 The groundbreaking discovery of monolayer graphene [1] has inspired an extensive quest to emulate massless Dirac quasiparticles in a myriad of distinct artificial graphene systems [2–11], ranging from ultracold atoms in optical lattices [3] to evanescently coupled photonic waveguide arrays [4]. Owing to their honeycomb symmetry, linear band-degeneracies man- ifest in the quasiparticle spectrum which we call conventional Dirac points (CDPs). These belong to the ubiquitous type- I class of two-dimensional (2D) Dirac points that are characterized by Dirac cones with closed isofrequency contours. As a result, the corresponding quasiparticles are described by the rather exotic 2D massless Dirac Hamiltonian [12], in stark contrast to the more conventional Schrodinger¨ equation. Thus, these emergent Dirac quasiparticles offer fundamental insight into pseudorelativistic phenomena such as the iconic Klein Figure 1. Schematic of a cavity-embedded honeycomb metasur- paradox, corresponding to the perfect transmission through ar- face. The honeycomb array of meta-atoms is composed of two in- bitrarily large potential barriers [13]. The latter is responsible equivalent (A andp B) hexagonal sublattices—definedp by lattice vec- for the suppression of backscattering and antilocalization of 3 3 3 3 tors a1 = a(− ; ) and a2 = a( ; )—which are connected 2 2 2 2 p Dirac quasiparticles, which are highly desirable properties for 3 1 by nearest-neighbor vectors e1 = a(0; −1), e2 = a( ; ) and p 2 2 efficient quasiparticle propagation in novel devices. 3 1 e3 = a(− 2 ; 2 ), where a is the subwavelength nearest-neighbor Since the existence of type-I CDPs is intrinsically linked separation. Each meta-atom is modelled as an electric dipole, ori- to the honeycomb structure, the fundamental properties of the ented normal to the plane of the lattice. The honeycomb metasurface massless Dirac quasiparticles are notoriously robust and diffi- is embedded inside a photonic cavity of height L, which is composed cult to manipulate. However, by exploiting meticulous control of two perfectly conducting metallic plates, enabling one to modify over the lattice structure, artificial graphene systems have en- the photonic environment while preserving the lattice structure. This abled the exploration of Dirac quasiparticles in new regimes general model can be realized across the electromagnetic spectrum, that are difficult, if not impossible to achieve in graphene from arrays of plasmonic nanorods to microwave helical resonators. itself [14–19]. An archetypal example which has attracted considerable interest is the paradigm of strain-engineering, where it has been shown that lattice anisotropy can induce the non-local interactions, mediated by the electromagnetic en- merging and annihilation of type-I CDPs [3, 14–16, 20–23], vironment, give rise to rich and unique Dirac physics at the and that aperiodicity can generate large pseudomagnetic fields subwavelength scale. In particular, we unveil the emergence [17, 24]. of a new kind of type-II Dirac point, which simultaneously ex- This concept of Hamiltonian engineering, and the hunt for ists with the type-I CDPs. In contrast to the latter, the former exotic quasiparticles, has recently entered the realm of polari- are characterized by critically tilted Dirac cones with open, tonics [25–30]. The true potential of polaritons lies in their hyperbolic isofrequency contours. hybrid nature, where their light and matter constituents can Recently, the discovery of type-II Dirac/Weyl semimet- be manipulated independently, thereby providing additional als has sparked a growing interest in the prospect of tunable degrees of freedom. As an example, recent work has type-II Dirac/Weyl points since the corresponding ‘Lorentz- shown the tantalizing prospect of engineering novel topologi- violating’ Dirac quasiparticles exhibit markedly different cal polaritons by introducing a winding coupling between or- properties from their type-I counterparts [31–35]. After dinary photons and excitons [27]. Shortly thereafter, these their realization, the pursuit for electromagnetic analogs com- non-trivial interactions were exploited to gap out photonic menced [36–42]. However, just like type-I CDPs, the ex- Dirac points as an alternative route to topological polaritons istence of type-II Dirac/Weyl points is usually tied to judi- [29]. ciously engineered lattice structures, and are thus impossible In this work we exploit the hybrid nature of polaritons in to manipulate without modifying the structure [36–42]. a different setting, namely metamaterials. We explore Dirac In stark contrast, here the existence of the type-II Dirac physics in an entirely new scenario, shifting the focus from the points is intrinsically linked to the hybrid nature of the po- lattice structure and its deformations to the effect of manipu- laritons, emerging from a non-trivial winding in the light- lating the surrounding photonic environment. In particular, matter interaction. Crucially, the combination of this unique we theoretically study the polaritons supported by crystalline physical origin and the truly 2D nature of the metasurface metasurfaces consisting of a honeycomb array of resonant, results in an unprecedented tunability of emergent Dirac po- dipolar meta-atoms. Such a structure is imminently realizable laritons by modifying only the surrounding photonic environ- across the electromagnetic spectrum from arrays of plasmonic ment. Specifically, we show that by embedding the honey- nanoparticles to microwave helical resonators. Despite the el- comb metasurface inside a planar photonic cavity and sim- ementary nature of these metasurfaces, we show that complex ply changing the cavity height, one can induce novel phase 3 transitions, including the multi-merging of type-I and type- Emergence of Type-I Dirac Points II Dirac points and the annihilation of type-II Dirac points. This striking tunability results in distinct polariton phases, de- The matter Hamiltonian within the nearest-neighbor approxi- spite the preserved lattice structure. In particular, we unveil mation reads a morphing between a linear and a parabolic spectrum ac- X y y ~ X y companied by a change in topological Berry phase, and an Hmat = ~!~0 aqaq + bqbq + ~Ω fqbqaq + H.c. ; unprecedented environment-induced inversion of chirality, all q q of which have no analog in graphene or artificial graphene (1) systems studied thus far. Therefore, this unique paradigm of where, for brevity, we have not presented the non-resonant !~ exploiting the photonic environment will give rise to novel terms (see Methods for derivation). In equation (1), 0 is the Ω~ Dirac-related phenomena at the subwavelength scale, such as renormalized resonant frequency and is the renormalized anomalous Klein tunnelling, negative refraction and pseudo- Coulombic interaction strength due to the cavity-induced im- magnetic Landau levels, which can all be tuned via the pho- age dipoles (see Methods for their dependence on the cav- ay by tonic environment alone. ity height). The bosonic operators q and q create quanta of the quasistatic collective-dipole modes that extend across the A and B sublattices, respectively, with wavevector q in the first Brillouin zone (see Fig. 2a). Finally, the function P3 fq = j=1 exp iq · ej encodes the honeycomb geometry RESULTS of the lattice with nearest-neighbor vectors ej (see Fig. 1). We diagonalize the matter Hamiltonian (1) as Hmat = P P !matβy β βy = Hamiltonian Formulation τ=± q ~ qτ qτ qτ where the bosonic operators qτ y qj qτ i create quasistatic collective-dipole normal modes with dispersion !mat =! ~ + τΩ~jf j. Here, τ indexes the While metamaterials have traditionally been described in qτ 0 q upper (τ = +1) and lower (τ = −1) bands and y = terms of macroscopic effective properties [36, 39, 43], the q (ay ; by ) is a spinor creation operator. The spinors j i = importance of crystallinity is becoming

Load more