Quadrupole topological photonic crystals

Download PDF Article Open access Published: 19 June 2020 Quadrupole topological photonic crystals Li He1, Zachariah Addison1, Eugene J. Mele1 & …Bo Zhen  ORCID:

orcid.org/0000-0001-8815-38931 Show authors Nature Communications volume 11, Article number: 3119 (2020) Cite this article

11k Accesses

Metrics details

Quadrupole topological phases, exhibiting protected boundary states that are themselves topological insulators of lower dimensions, have recently been of great interest. Extensions of these

ideas from current tight binding models to continuum theories for realistic materials require the identification of quantized invariants describing the bulk quadrupole order. Here we

identify the analog of quadrupole order in Maxwell’s equations for a gyromagnetic photonic crystal (PhC) through a double-band-inversion process. The quadrupole moment is quantized by the

simultaneous presence of crystalline symmetry and broken time-reversal symmetry, which is confirmed using three independent methods: analysis of symmetry eigenvalues, numerical calculations

of the nested Wannier bands and the expectation value of the quadrupole operator. Furthermore, we reveal the boundary manifestations of quadrupole phases as quantized edge polarizations and

fractional corner charges. The latter are the consequence of a filling anomaly of energy bands as first predicted in electronic systems.

quadrupole topological insulator using orbital-induced synthetic flux Article Open access 03 November 2022 Real higher-order Weyl photonic crystal Article Open access 20 October 2023

Bulk-local-density-of-state correspondence in topological insulators Article Open access 14 November 2023 Introduction

Symmetries play a pivotal role in understanding and classifying various topological phases of matter1,2,3,4. In periodic media, systems with different symmetries can admit different

classifications characterized by quantized topological invariants. Furthermore, interface states in symmetry-protected topological (SPT) phases can only robustly exist, between two

topologically distinct regions, when required bulk symmetries are preserved at the boundaries. In the simplest example of SPT phases, a one-dimensional (1D) Su-Schrieffer-Heeger (SSH) model,

the topological invariant—Zak phase5—is only quantized when inversion symmetry is preserved, leading to edge states that are protected by non-zero bulk polarizations. Recently, the concept

of polarization, or bulk dipole moment in crystals, has been generalized to include multipole moments, such as quadrupole and octupole moments, leading to the discovery of higher-order

topological insulators (HOTIs)6,7,8,9,10. In particular, these HOTIs have vanishing dipole densities but non-zero higher-order multiple moments, which can be quantized by certain crystalline

symmetries, such as reflection and rotation. In contrast to conventional topological insulators (TIs) that support gapless boundary states, HOTIs exhibit protected boundaries that are,

themselves, TIs in lower dimensions.

Recently, quadrupole TIs have been demonstrated in a number of classical systems, ranging from microwave11 and optics12 to acoustics13 and circuits14,15. Most experiments are analyzed as

lattice models, following the tight-binding approximation. In these lattice models, the adopted symmetry required to quantize the quadrupole moment is reflection or fourfold rotation (C4)

with threaded π-flux, and hence these systems are time-reversal invariant. Unfortunately, most realistic systems with subwavelength features, such as photonic crystals, cannot be modeled

discretely but require a different continuum approach. For example, the practical implementations of many lattice models no long preserve quantized bulk quadrupole moments in the continuum

theory. In addition, most experimental works solely used the existence of in-gap corner states as the measure of bulk quadrupole topology. The validity relies on the presence of additional

chiral symmetry in the lattice model, which is often not preserved in the continuum theory.

Here, we find solutions to continuous Maxwell's equations in gyromagnetic photonic crystals that are the electrodynamic analogs to quadrupole topological phases. The proposed topological

PhCs have quantized bulk quadrupole moments, which are protected by the simultaneous presence of crystalline symmetries and broken time-reversal symmetry (T)16. In particular, we show it is

essential to break time-reversal symmetry—to open the energy gap—while preserving crystalline symmetries—to quantize bulk dipole and quadrupole moments. In addition, we note the existence of

corner states is neither sufficient nor necessary condition for quadrupole phases. Instead, we validate the bulk quadrupole nature through analyzing the Wannier band polarization and its

manifestations at boundaries as quantized edge polarizations and fractional corner charges. All calculations are based on realistic parameters that are readily available in the microwave

regime17,18.

We start by presenting the topological phase transition between a trivial (Fig. 1b) and a quadrupole two-dimensional (2D) PhC (Fig. 1c). The 2D PhC consists of gyromagnetic rods in air and

is homogeneous along the out-of-plane (z) direction. Experimentally, this boundary condition can also be realized using metals17,18. The gapless transition point is achieved in a 2 × 2

super-cell structure with four square rods (Fig. 1a). All rods are identical in shape and are of the same gyromagnetic material of Yttrium Iron Garnet (YIG) with isotropic dielectric

permittivity of ϵ = 15ϵ0 and inplane permeability μ = 14μ0. To break time-reversal symmetry, an external magnetic field is applied along the out-of-plane direction (z), which induces

complex-valued off-diagonal terms in the permeability tensor of YIG19:

}_{0}\end{array}\right]$$ (1)

This gyromagnetic response, \({\mu }_{xy}={\mu }_{yx}^{* }=i\kappa\), breaks T but preserves C4 and Mx(y)T. Here Mx(y) is mirror reflection that transforms x(y) to −x(y). At the phase

transition, all rods are placed at a/2 away from their neighbors; the corresponding band structure for TM modes (Ez, Hx, Hy) has twofold (fourfold) degeneracies at the center (corner) of the

folded Brilluion zone Γ (M). Both degeneracies are lifted when the four rods are simultaneously displaced inward (d < 0, Fig. 1b) or outward (d > 0, Fig. 1c) along the diagonal lines of the

unit cell. On either side of the transition point, the band structure is fully gapped owing to T-breaking (shaded in yellow) and supports two topologically distinct phases determined by the

displacement d: for the choice of cell in Fig. 1, inward displacements with negative d give rise to trivial phases, whereas outward displacements with positive d correspond to quadrupole

phases.

a Band degeneracies are created at the center (Γ) and corner (M) of the folded Brilluion zone for a 2 × 2 super-cell PhC (inset). An external magnetic field (B) induces gyromagnetic

responses in YIG rods (gray squares), which breaks T but preserves C4. The gyromagnetic response κ = 12.4μ0 is used in the calculations. b, c As the rods simultaneously move inward or

outward, the second gap is opened (shaded in yellow), but with different quadrupole phases: inward (outward) displacements with negative (positive) d correspond to trivial (non-trivial)

quadrupole phases. d, e Ez mode profiles for the second and third modes at M feature a band inversion between the modes with C4 index of +i (green) and −i (red) in the trivial and

non-trivial phases. f Complete topological phase diagram of the PhCs, including quadrupole phases (green), trivial phases (orange), and Chern insulating phases (purple). Systems are gapless

along the axes (red). NI normal insulator, QTI quadrupole insulator, CI Chern insulator.

Next, we present our calculations of the quadrupole topological invariant qxy, using three different approaches, and demonstrate the phase transition between qxy = 0 and qxy = 1/2 by

displacing the dielectric rods. We start by evaluating qxy using the C4 eigenvalues at high-symmetry k points of all bands below the gap6,16:

){r}_{4}^{+* }({\rm{M}})={r}_{4}^{-}(\Gamma ){r}_{4}^{-* }({\rm{M}}).$$ (2)

Here, \({r}_{4}^{+}\) (\({r}_{4}^{-}\)) is the C4 eigenvalue of a mode with C2 eigenvalue r2 = +1(−1); naturally, \({r}_{4}^{+}={\pm}\!{1}\) and \({r}_{4}^{-}={\pm}\!{i}\). Accordingly, a

quadrupole topological phase transition happens when two pairs of bands switch their C4 eigenvalues at the same time—a process we call “double-band-inversion”. Specifically, the double-band

inversion happens in our system when d changes from negative to positive: through this process, the TM mode at M with a phase winding of +2π in the Ez mode profile (r4 = i, labeled as green)

switches position with the one with the winding of −2π (r4 = −i, label in red); meanwhile, the two modes with r4 = ±1 at M also switch positions. On the other hand, all C4 eigenvalues at Γ

remain unchanged. Using Eq. (2), we identify PhCs with d  > 0 to be topologically non-trivial, with bulk quadrupole invariant qxy = 1/2, and PhCs with d < 0 to be topologically trivial with

qxy = 0.

A more-comprehensive topological phase diagram of our system is shown in Fig. 1f, determined by displacement d and strength of gyromagnetic response κ. Both quadrupole topological insulators

(green) and trivial insulators (orange) are identified, with the super-cell structure (d = 0) being the transition in between the two phases (y axis). Interestingly, Chern insulating phases

(purple) are observed at large displacements, consistent with the broken time-reversal symmetry19,20.

To confirm the quadrupole topology of our PhC, we explicit show that its dipole moments are zero (px = py = 0), but its quadrupole moment is non-zero (qxy = 1/2). To this end, we perform two

separate sets of calculations using two different methods: (1) the nested Wilson loop formulation6; (2) and the expectation value of the exponentiated quadrupole operator21,22, and show

they reach the same conclusions. Here, we present the nested Wilson loop calculations, leaving the second method in Supplementary Note 3. We start by computing the band structure and mode

profiles (Ez) of the TM modes for a particular displacement of dielectric rods d = a/4 − w/2 (Fig. 2a, same as in Fig. 1c). Using these as input, we compute the Wannier bands νx,y for the

two lowest-energy bands (Supplementary Note 1).

a Wannier bands, νx(ky) and νy(kx), are calculated for the first and second bands of the quadrupole PhC. Results show the bulk dipole momentum is zero (px = py = 0). b Calculations of the

nested Wilson loops, \({p}_{y}^{{\nu }_{x}^{-}}={p}_{x}^{{\nu }_{y}^{-}}=0.5\), show the bulk quadrupole momentum qxy is non-trivial and quantized to 0.5. c, d Similar calculations repeated

for a trivial PhC, showing zero bulk dipole moments (px = py = 0) and zero quadrupole moments (\({p}_{y}^{{\nu }_{x}^{-}}={p}_{x}^{{\nu }_{y}^{-}}=0\)).

The results (Fig. 2a) show that the Wannier bands are related as \(\{{\nu }_{x}^{j}({k}_{y})\}\to \{-{\nu }_{x}^{j}(-{k}_{y})\}\) mod 1 owing to the presence of C2 symmetry (Supplementary

Note 2). Therefore, the bulk polarization px, which is simply the integral of the two Wannier bands over the 1D B.Z., vanishes: \({p}_{x}={\sum }_{n = {\pm}

}{\int}_{{\rm{B}}.{\rm{Z}}.}d{k}_{y}\,{\nu }_{x}^{n}({k}_{y})=0\). A similar argument can also be applied to bulk polarization py, proving the bulk dipole moments are zero: px = py = 0.

Meanwhile, owing to the existence of a gap in the Wannier spectrum, the upper and lower Wannier bands can be separated into two sectors, labeled as \({\nu }_{x}^{\pm }\) (\({\nu }_{y}^{\pm

}\)), respectively. To confirm the quadrupole topology, we compute the polarizations of the Wannier bands within one sector, \({p}_{y}^{{\nu }_{x}^{-}}\) and \({p}_{x}^{{\nu }_{y}^{-}}\),

using the nested Wilson loop formulation6,23. As shown in Fig. 2b, both Wannier band polarizations are quantized to be 1/2—due to the simultaneous preservation of MxT and MyT symmetries

(Supplementary Note 2)—which further confirms our PhC has a non-zero bulk quadrupole moment: \({q}_{xy}=2{p}_{x}^{{\nu }_{y}^{-}}{p}_{y}^{{\nu }_{x}^{-}}=1/2\).

To compare, we repeat the same set of calculations for a different unit cell with a negative displacement of d = −(a/4 − w/2) (Fig. 2c). As shown, the Wannier bands are also gapped and sum

to zero, meaning the bulk dipoles px,y remain zero. However, the nested Wilson loops, \({p}_{y}^{{\nu }_{x}^{-}}\) and \({p}_{x}^{{\nu }_{y}^{-}}\), are also zero, leading to a trivial

quadrupole moment qxy = 0. This correspondence between positive (negative) displacements leading to non-trivial (trivial) quadrupole phases is consistent with our preceding conclusions based

on C4 symmetry eigenvalues. Importantly, the non-trivial and trivial PhCs in Fig. 2 can be simply related to each other, by shifting the choice of the unit cell center by a/2. This

observation leads to an intuitive understanding of the difference in quadrupole moment between the two phases as discussed in Supplementary Note 5.

bands

Next, we present the physical consequences of quadrupole topological PhCs at interfaces, originating from the bulk-edge correspondence of the Wannier bands. Following classical

electromagnetism, a non-zero bulk quadrupole moment in a finite system is manifested as edge polarizations at its boundaries. Here, we study the 1D interfaces between quadrupole (trivial)

PhCs and perfect electric conductors (PECs). We find that the non-trivial quadrupole topology indeed leads to a quantized edge polarization along the interface, which is absent for a trivial

PhC. Specifically, we consider a strip of 20 unit cells of the quadrupole (trivial) PhC design, which satisfies periodic boundary condition in the x—direction and closed boundary condition

in the y—direction, owing to the two PECs on top and bottom as shown in Fig. 3a, b. The energy dispersions ω(kx) of the quadrupole and trivial strips (Fig. 3c, d) share the same bulk energy

spectra (gray areas) but have different edge dispersions (gray solid lines), owing to their different edge terminations. The Wannier centers νx are calculated for the two different strips

based on their energy dispersions and Bloch mode profiles (Fig. 3e, f). For the topological strip, two additional Wannier states (red circles) are found to emerge outside the Wannier bands

(blue) at the middle of the gap (±0.5), which is protected by the non-trivial bulk quadrupole moment. This can be understood in a similar way as the in-gap edge states in the 1D SSH model

with non-trivial bulk polarization. The quantization of the mid-gap Wannier states at ±0.5 is due to the additional MxT symmetry in our system, which is retained even at the boundary of a

finite strip as shown in Fig. 3a. (Supplementary Note 2). In comparison, no additional Wannier states are found inside the Wannier gaps for the trivial strip, consistent with the lack of

bulk quadrupole moment.

a, b Two 1D interfaces are created between a strip of quadrupole (trivial) PhC (same as Fig. 1b,c) and perfect electric conductors (PECs). c, d Energy dispersions of the two strip setups. e

For the non-trivial setup, Wannier centers of the eigenmodes show two Wannier states at ±0.5, which are outside the bulk Wannier bands owing to the non-trivial bulk quadrupole moment. f In

comparison, no in-gap Wannier states are observed in the trivial setup. g Calculation of the spatial density of polarization px(y) shows the in-gap Wannier states are localized at the top (y

 = 20a) and bottom interfaces (y = 0). h Similar calculations repeated for the trivial PhC, showing no edge polarization, which is consistent with the trivial quadrupole moment.

image

Furthermore, the two mid-gap Wannier states in the topological strip are spatially localized at the top and bottom edges. To demonstrate this, we further study the spatial distribution of

Wannier states by calculating the polarization density px as a function of position along y (Supplementary Note 6). In order to choose a definite sign of the polarization, we introduce an

infinitesimal perturbation to break the C2 symmetry of the semi-infinite strip. For the quadrupole PhC, as shown in Fig. 3g, there are non-zero polarization densities developed near the two

edges at y = 0 and y = 20a, and the edge polarizations are quantized to ±0.5. In comparison, neither edge nor bulk polarization are observed in the trivial PhC (Fig. 3h).

corner charges and filling anomaly

Finally, we show the physical consequence of quadrupole PhCs as localized 0D corner states, which are the photonic analogs of states responsible for filling anomalies and fractional charges

in an electronic setting24. To this end, we solve the eigenstates in a finite 2D quadrupole PhC enclosed by PECs, with a thin air gap in between (Fig. 4a). The eigenstates are labeled

according to their energies, with the lowest-energy state labeled as 1. Aside from delocalized bulk states, four degenerate states—199 through 202—are found, with their energies inside the

bulk energy gap (Fig. 4b). Their mode profiles (Ez) confirm that these states are spatially localized at the four corners, with one example shown in Fig. 4c. Owing to the lack of chiral

symmetry in our system, which is generic in the continuum theory expanded around non-zero frequencies, the energy of the four corner states is not pinned at the center of the bulk energy

gap25; instead they can be shifted, even immersed into the bulk continuum, by modifying the refractive index at the corners. This renders the appearance of corner states less of an essential

signature of quadrupole phases10,11. In fact, corner states are also found in other higher-order topological phases with vanishing bulk quadrupole moments23,24,25,26,27,28,29,30.

Topologically protected corner states and filling anomaly in a 2D system.

a For a finite PhC with N×N unit cells enclosed by PECs, the counting of the number of eigenstates below and above an energy gap are distinct between quadrupole and trivial phases. b

Eigenstates of a 10 × 10 quadrupole PhC, showing four degenerate corner states inside the bulk gap, shaded in yellow. c Mode profile (∣Ez∣) of one of the corner states. d Spatial

distribution of the accumulative electromagnetic energy density for the lowest 200 energy eigenstates in the quadrupole PhC, showing fractional occupations (2 ± 0.5) at the corners.

size image

Instead, here we illustrate the non-trivial quadrupole nature of our PhC using the filling anomaly, by counting the number of energy eigenstates below and above a given bulk energy gap10

(Fig. 4a). Specifically, even though trivial samples may support corner states, they originate from either the top or bottom band alone, leaving 2N2−4n states in the bulk continuum (n is an

integer). On the other hand, for non-trivial PhCs, the number of states below and above the energy gap are both 2N2 − 2 = 198 (N = 10 in our case). This indicates the four degenerate corner

states we observed are “contributed” by both the top and bottom bands together, and thus proves the non-trivial quadrupole topology of our design. As a consequence, for a quadrupole PhC,

quantized fractional charges appear at four corners of a finite sized system (Fig. 4d) when calculate the spatial distribution of the lowest 200 energy states, in a similar vein as the

calculations of charge density in electronic systems at “half-filling” (here, we have introduced an infinitesimal perturbation to break C4 symmetry in order to split the fourfold degenerate

corner states). Remarkably, these corner charges are shared by two convergent dipoles on the two perpendicular edges, as the magnitude of the corner charges is equal to the edge

polarizations. This further confirms that these corner charges originate from a non-zero bulk quadrupole moment. We point out that the observed fractional corner charges arise from the

fundamental difference in the counting of bulk states10, and was recently proposed in electronic systems by Benalcazar et al.24 as a filling anomaly: a mismatch between number of states in

an energy band and the number of electrons required for charge neutrality.

In summary, we present quadrupole topological photonic crystals with truly quantized invariants and the physical consequences at material’s edges and corners. The proposed gyromagnetic PhCs

can be readily realized in the microwave regime. Meanwhile, the coexistence of multiple topological phases in our system, both quadrupole TIs and Chern insulators, provides a versatile

platform to further demonstrate topological photonic circuits with protected elements immune to disorders in various dimensions. Finally, our findings of inducing quadrupole phase

transitions and quantizing quadrupole moments—through crystalline symmetries in conjunction with broken time-reversal symmetry—can also be applied to other wave systems, including electrons,

phonons, and polaritons.

The band structures and mode profiles are calculated using the Finite Element Method in COMSOL Multiphysics 5.4. Specifically, we compute the band structures and mode profiles in a 2D

geometry with periodic boundary conditions along all directions. The corresponding Wannier bands are calculated by using the Bloch mode profiles as input (Supplementary Note 1).

availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

The code that supports the plots within this paper and other findings of this study is available from the corresponding authors upon reasonable request.

References Chiu, C.-K., Teo, J. C. Y., Schnyder, A. P. & Shinsei, R. Classification of topological quantum matter with symmetries. Rev. Mod. Phys. 88, 035005 (2016).

Article  ADS  Google Scholar 

Kitaev, A. Periodic table for topological insulators and superconductors. AIP conference proceedings 1134, 22–30 (2009).

Ryu, S., Schnyder, A. P., Furusaki, A. & Ludwig, A. W. W. Topological insulators and superconductors: tenfold way and dimensional hierarchy. N. J. Phys. 12, 065010 (2010).

Article  Google Scholar 

Schnyder, A. P., Ryu, S., Furusaki, A. & Ludwig, A. W. W. Classification of topological insulators and superconductors in three spatial dimensions. Phys. Rev. B 78, 195125 (2008).

Article  ADS  Google Scholar 

Zak, J. Berry’s phase for energy bands in solids. Phys. Rev. Lett. 62, 2747 (1989).

Article  ADS  CAS  Google Scholar 

Benalcazar, W. A., Bernevig, B. A. & Hughes, T. L. Electric multipole moments, topological multipole moment pumping, and chiral hinge states in crystalline insulators. Phys. Rev. B 96,

245115 (2017).

Article  ADS  Google Scholar 

Benalcazar, W. A., Bernevig, B. A. & Hughes, T. L. Quantized electric multipole insulators. Science 357, 61–66 (2017).

Article  ADS  MathSciNet  CAS  Google Scholar 

Langbehn, J., Peng, Y., Trifunovic, L., von Oppen, F. & Brouwer, P. W. Reflection-symmetric second-order topological insulators and superconductors. Phys. Rev. Lett. 119, 246401 (2017).

Article  ADS  Google Scholar 

Schindler, F. et al. Higher-order topological insulators. Sci. Adv. 4, eaat0346 (2018).

Article  ADS  Google Scholar 

Song, Z., Fang, Z. & Fang, C. (d- 2)-dimensional edge states of rotation symmetry protected topological states. Phys. Rev. Lett. 119, 246402 (2017).

Article  ADS  Google Scholar 

Peterson, C. W., Benalcazar, W. A., Hughes, T. L. & Bahl, G. A. quantized microwave quadrupole insulator with topologically protected corner states. Nature 555, 346–350 (2018).

Article  ADS  CAS  Google Scholar 

Mittal, S. et al. Photonic quadrupole topological phases. Nat. Photonics 13, 692–696 (2019).

Article  ADS  CAS  Google Scholar 

Serra-Garcia, M. et al. Observation of a phononic quadrupole topological insulator. Nature 555, 342–345 (2018).

Article  ADS  CAS  Google Scholar 

Imhof, S. et al. Topolectrical-circuit realization of topological corner modes. Nat. Phys. 14, 925–929 (2018).

Article  CAS  Google Scholar 

Serra-Garcia, M., Süsstrunk, R. & Huber, S. D. Observation of quadrupole transitions and edge mode topology in an lc circuit network. Phys. Rev. B 99, 020304 (2019).

Article  ADS  CAS  Google Scholar 

Wieder, B. J. et al. Strong and fragile topological dirac semimetals with higher-order fermi arcs. Nat. Commun. 11, 1–13 (2020).

Article  Google Scholar 

Skirlo, S. A. et al. Experimental observation of large chern numbers in photonic crystals. Phys. Rev. Lett. 115, 253901 (2015).

Article  ADS  Google Scholar 

Wang, Z., Chong, Y., Joannopoulos, J. D. & Soljačić, M. Observation of unidirectional backscattering-immune topological electromagnetic states. Nature 461, 772–775 (2009).

Article  ADS  CAS  Google Scholar 

Wang, Z., Chong, Y. D., Joannopoulos, J. D. & Soljačić, M. Reflection-free one-way edge modes in a gyromagnetic photonic crystal. Phys. Rev. Lett. 100, 013905 (2008).

Article  ADS  Google Scholar 

Skirlo, S. A., Lu, L. & Soljačić, M. Multimode one-way waveguides of large chern numbers. Phys. Rev. Lett. 113, 113904 (2014).

Article  ADS  Google Scholar 

Kang, B., Shiozaki, K. & Cho, G. Y. Many-body order parameters for multipoles in solids. Phys. Rev. B 100, 245134 (2019).

Article  ADS  CAS  Google Scholar 

Wheeler, W. A., Wagner, L. K. & Hughes, T. L. Many-body electric multipole operators in extended systems. Phys. Rev. B 100, 245135 (2019).

Article  ADS  CAS  Google Scholar 

Xie, B.-Y. et al. Second-order photonic topological insulator with corner states. Phys. Rev. B 98, 205147 (2018).

Article  ADS  CAS  Google Scholar 

Benalcazar, W. A., Li, T. & Hughes, T. L. Quantization of fractional corner charge in c n-symmetric higher-order topological crystalline insulators. Phys. Rev. B 99, 245151 (2019).

Article  ADS  CAS  Google Scholar 

Noh, J. et al. Topological protection of photonic mid-gap defect modes. Nat. Photonics 12, 408–415 (2018).

Article  ADS  CAS  Google Scholar 

Chen, X.-D. et al. Direct observation of corner states in second-order topological photonic crystal slabs. Phys. Rev. Lett. 122, 233902 (2019).

Article  ADS  CAS  Google Scholar 

Hassan, A. E. L. et al. Corner states of light in photonic waveguides. Nat. Photonics 13, 697–700 (2019).

Article  ADS  Google Scholar 

Ni, X., Weiner, M., Alù, A. & Khanikaev, A. B. Observation of higher-order topological acoustic states protected by generalized chiral symmetry. Nat. Mater. 18, 113–120 (2019).

Article  ADS  CAS  Google Scholar 

Xue, H., Yang, Y., Gao, F., Chong, Y. & Zhang, B. Acoustic higher-order topological insulator on a kagome lattice. Nat. Mater. 18, 108–112 (2019).

Article  ADS  CAS  Google Scholar 

Zhang, X. et al. Second-order topology and multidimensional topological transitions in sonic crystals. Nat. Phys. 15, 582–588 (2019).

Article  CAS  Google Scholar 

Download references

This work was partly supported by the National Science Foundation through the University of Pennsylvania Materials Research Science and Engineering Center DMR-1720530 and the US Office of

Naval Research (ONR) Multidisciplinary University Research Initiative (MURI) grant N00014-20-1-2325 on Robust Photonic Materials with High-Order Topological Protection. L.H. was supported by

the Air Force Office of Scientific Research under award number FA9550-18-1-0133. Work by Z.A. and E.M. on the symmetry analysis of the quadrupole phase was supported by the Department of

Energy, Office of Basic Energy Sciences under grant DE FG02 84ER45118. B.Z. was supported by the Army Research Office under award contract W911NF-19-1-0087.

Author informationAuthors and Affiliations Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA, 19104, USA

Li He, Zachariah Addison, Eugene J. Mele & Bo Zhen

AuthorsLi HeView author publications You can also search for this author inPubMed Google Scholar

Zachariah AddisonView author publications You can also search for this author inPubMed Google Scholar

Eugene J. MeleView author publications You can also search for this author inPubMed Google Scholar

Bo ZhenView author publications You can also search for this author inPubMed Google Scholar

L.H. and B.Z. conceived the idea. L.H. carried out numerical simulations. L.H., Z.A., E.M., and B.Z. discussed and interpreted the results. L.H and B.Z. wrote the paper with contribution

from all authors. B.Z. supervised the project.

Corresponding author Correspondence to Bo Zhen.

The authors declare no competing interests.

Peer review information Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

informationPeer Review fileRights 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/.

Reprints and permissions

About this articleCite this article He, L., Addison, Z., Mele, E.J. et al. Quadrupole topological photonic crystals. Nat Commun 11, 3119 (2020). https://doi.org/10.1038/s41467-020-16916-z

Download citation

Received: 14 November 2019

Accepted: 29 May 2020

Published: 19 June 2020

DOI: https://doi.org/10.1038/s41467-020-16916-z

Share this article Anyone you share the following link with will be able to read this content:

Get shareable link Sorry, a shareable link is not currently available for this article.

Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative