This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Antihelical Edge States in Two-dimensional Photonic Topological Metals

Liangcai Xie    Liang Jin jinliang@nankai.edu.cn    Zhi Song School of Physics, Nankai University, Tianjin 300071, China
Abstract

Topological edge states are the core of topological photonics. Here we introduce the antihelical edge states of time-reversal symmetric topological metals and propose a photonic realization in an anisotropic square lattice of coupled ring resonators, where the clockwise and counterclockwise modes play the role of pseudospins. The antihelical edge states robustly propagate across the corners toward the diagonal of the square lattice: The same (opposite) pseudospins copropagate in the same (opposite) direction on the parallel lattice boundaries; the different pseudospins separate and converge at the opposite corners. The antihelical edge states in the topological metallic phase alter to the helical edge states in the topological insulating phase under a metal-insulator phase transition. The antihelical edge states provide a unique manner of topologically-protected robust light transport applicable for topological purification. Our findings create new opportunities for topological photonics and metamaterials.

Introduction.—The fundamental concepts in condensed matter physics introduced to topological photonics inspire the rapid development of photonic topological states LLu14 ; ABKNPo17 ; TORMP19 . The chiral edge states of topological insulators unidirectionally propagate along boundaries and require the breaking of time-reversal symmetry FHPRL08 ; ZWNa09 ; MHNP11 ; KFNP12 ; GQLPRL13 ; SZhangPRL20 . The degenerate clockwise and counterclockwise modes of ring resonators experience opposite artificial magnetic fields and provide a pseudospin degree of freedom MHNP11 . The helical edge states of time-reversal symmetric topological insulators with different pseudospins unidirectionally propagate in opposite directions. The edge states of topological metals are topologically-protected in the gapless phase. Antichiral edge states have been proposed and implemented on zigzag edges by modifying the next-nearest-neighbor hopping phase of the Haldane model FranzPRL18 ; ZhangBLPRL20 ; ChongYD21 ; LiZYPRL22 . Antichiral edge states propagate along the parallel lattice boundaries in the same direction. Recent progress in antichiral edge states has led to the development of a solution for creating topological metals. Topological metals are difficult to create because of the challenge of separating gapless bands and breaking time-reversal symmetry in photonics. Thus, is it possible to have time-reversal symmetric topological metals with topologically-protected edge states and robust propagation?

Here, we introduce the antihelical edge states [Fig. 1(a)] of time-reversal symmetric topological metals and a photonic realization is proposed in the two-dimensional anisotropic square lattice of coupled ring resonators. The pseudospins are time-reversal symmetric counterparts, and the introduction of pseudospins addresses the difficulty of breaking time-reversal symmetry in photonics. The symmetric component of next-nearest-neighbor couplings creates a nontrivial topology, and the anti-symmetric component of next-nearest-neighbor couplings separates energy bands and supports the antihelical edge states on both horizontal and vertical boundaries. Antihelical edge states robustly copropagate along corners with opposite pseudospins that separate and converge at the opposite corners on the diagonal of a lattice. The antihelical edge states become helical edge states during a metal-insulator phase transition [Fig. 1(b)]. Photonic topological metals provide a new direction for research on topological photonics.

Refer to caption
Figure 1: (a) Robust two-dimensional propagation of chiral, helical, antichiral, and antihelical edge states. Time-reversal symmetry is broken for chiral and antichiral edge states, and time-reversal symmetry is preserved for helical and antihelical edge states. Chiral and helical edge states are in the insulating phases; antichiral and antihelical edge states are in the metallic phases. Different spins of antihelical edge states separate and converge at the opposite corners on the diagonal. (b) Schematic of band and edge state propagation for topological insulator, metal-insulator phase transition, topological metal, and topological phase transition. The direction and length of the arrow indicate the velocity of edge states. At the TPT point, the velocities in the horizontal direction on the top and bottom boundaries are identical; across the TPT point from C=+1C_{\uparrow}=+1 to C=1C_{\uparrow}=-1, the velocity on the bottom boundary exceeds the velocity on the top boundary. (c) Schematic of two-dimensional anisotropic square lattice of coupled ring resonators. (d) Schematic of Hamiltonian for pseudospin-up. (e) Schematic of Hamiltonian for pseudospin-down. For χ=κ/2=1/2\chi=\kappa/2=1/2, J1=1J_{1}=1, J2=1/2J_{2}=1/2, and δ=3/2\delta=3/2, (f) nonzero χ\chi opens the band gap, (g) is the phase diagram for the pseudospin-up, and (h) is the distribution of edge states.

Anisotropic square lattice of ring resonators.—Figure 1(c) presents a two-dimensional anisotropic square lattice of coupled ring resonators (see Supplementary materials A). The clockwise mode (pseudospin-up) and counterclockwise mode (pseudospin-down) are time-reversal counterparts, and they experience the opposite Peierls phases in the horizontal couplings between the nearest-neighbor resonators MHNP11 . The lattice for the pseudospin-up (pseudospin-down) is presented in Fig. 1(d) [Fig. 1(e)]. The horizontal couplings ±iJ1\pm iJ_{1} indicated in green and ±iJ2\pm iJ_{2} indicated in orange are tunneling-direction-dependent, break the time-reversal symmetry of the Hamiltonian for each individual pseudospin, and separately affect the edge states localized on the upper and lower boundaries. The vertical coupling between the nearest-neighbor resonators is κ\kappa indicated in black. The cross couplings between the next-nearest-neighbor resonators are χ1=χ+δ\chi_{1}=\chi+\delta, χ2=χδ\chi_{2}=\chi-\delta LeykamPRL18 . The symmetric component χ\chi opens a band gap [Fig. 1(f)] and creates a nontrivial topology, whereas the anti-symmetric component δ\delta affects the band structure, separates the bands in both the xx and yy directions, and creates antihelical edge states that are localized on the left and right boundaries in the topological metallic phase.

Topological phases.—The Bloch Hamiltonian for the pseudospin-up is h(𝐤)=d0(𝐤)σ0+𝐝(𝐤)σh_{\uparrow}\left(\mathbf{k}\right)=d_{0}\left(\mathbf{k}\right)\sigma_{0}+\mathbf{d\left(\mathbf{k}\right)\cdot\sigma}, where σ0\sigma_{0} is the identical matrix and σ=(σx,σy,σz)\mathbf{\sigma=}\left(\sigma_{x},\sigma_{y},\sigma_{z}\right) is the Pauli matrix. The first term adjusts the band energy and the second term determines the band topology. We have d0(𝐤)=(J1+J2)sinkxd_{0}\mathbf{\left(\mathbf{k}\right)}=(J_{1}+J_{2})\sin k_{x} and 𝐝(𝐤)=𝐫1(kx)𝐫2(ky)\mathbf{d}\left(\mathbf{k}\right)=\mathbf{r}_{1}(k_{x})-\mathbf{r}_{2}(k_{y}) with 𝐫1(kx)=(κ+2χcoskx,2δsinkx,(J1J2)sinkx)\mathbf{r}_{1}(k_{x})=(\kappa+2\chi\cos k_{x},-2\delta\sin k_{x},(J_{1}-J_{2})\sin k_{x}), 𝐫2(ky)=(κcosky,κsinky,0)\mathbf{r}_{2}(k_{y})=(-\kappa\cos k_{y},\kappa\sin k_{y},0). h(𝐤)h_{\uparrow}\left(\mathbf{k}\right) respects the particle-hole symmetry σzhT(𝐤)σz1=h(𝐤)\sigma_{z}h_{\uparrow}^{T}\left(\mathbf{k}\right)\sigma_{z}^{-1}=-h_{\uparrow}\left(-\mathbf{k}\right). The topological phase belongs to the class DD and it is characterized by a \mathbb{Z} topological invariant, that is, the spin-Chern number C=(4π)1|𝐝|3𝐝(kx𝐝×ky𝐝)𝑑kx𝑑kyC_{\uparrow}=-\left(4\pi\right)^{-1}\int\left|\mathbf{d}\right|^{-3}\mathbf{d}\cdot\left(\partial_{k_{x}}\mathbf{d}\times\partial_{k_{y}}\mathbf{d}\right)dk_{x}dk_{y} [Fig. 1(g)]. The band topology is captured by the effective magnetic field 𝐝(𝐤)\mathbf{d}\left(\mathbf{k}\right). After 𝐝(𝐤)\mathbf{d}\left(\mathbf{k}\right) is substituted with 𝐫1(kx)𝐫2(ky)\mathbf{r}_{1}(k_{x})-\mathbf{r}_{2}(k_{y}) in CC_{\uparrow}, the spin-Chern number becomes the definition of a linking number of the two independent periodic vectors 𝐫1(kx)\mathbf{r}_{1}(k_{x}) and 𝐫2(ky)\mathbf{r}_{2}(k_{y}) (see Supplementary materials B). 𝐫1(kx)\mathbf{r}_{1}(k_{x}) is an ellipse that passes through (κ2χ,0,0)\left(\kappa-2\chi,0,0\right) and (κ+2χ,0,0)\left(\kappa+2\chi,0,0\right). 𝐫2(ky)\mathbf{r}_{2}(k_{y}) is a circle centered at the origin with a fixed radius |κ|\left|\kappa\right| in the z=0z=0 plane. Thus, the two closed curves 𝐫1(kx)\mathbf{r}_{1}(k_{x}) and 𝐫2(ky)\mathbf{r}_{2}(k_{y}) are linked when 0<|χ|<|κ|0<\left|\chi\right|<\left|\kappa\right|.

In Fig. 2(a), the red circle 𝐫2(ky)\mathbf{r}_{2}(k_{y}) that is fixed on the z=0z=0 plane is always clockwise; the rotation direction of the blue ellipse 𝐫1(kx)\mathbf{r}_{1}(k_{x}) reverses along the topological phase transition plane J1=J2J_{1}=J_{2}, where the band gap vanishes and 𝐫1(kx)\mathbf{r}_{1}(k_{x}) becomes coplanar to 𝐫2(ky)\mathbf{r}_{2}(k_{y}). The blue ellipse 𝐫1(kx)\mathbf{r}_{1}(k_{x}) circles counterclockwise around the red circle 𝐫2(ky)\mathbf{r}_{2}(k_{y}) once in the region J1>J2J_{1}>J_{2} (C=+1C_{\uparrow}=+1); and the blue ellipse 𝐫1(kx)\mathbf{r}_{1}(k_{x}) circles clockwise around the red circle 𝐫2(ky)\mathbf{r}_{2}(k_{y}) once in the region J1<J2J_{1}<J_{2} (C=1C_{\uparrow}=-1). The right-hand (left-hand) rule identifies C=+1C_{\uparrow}=+1 (C=1C_{\uparrow}=-1): orient the thumb pointing along the arrow of the red circle, curl the rest of the fingers pointing along the arrow of the blue ellipse. Notably, in the phase diagram [Fig. 2(b)], two topological phases with opposite spin-Chern numbers have an identical number of edge states.

Antihelical edge states.—Topological edge states appear in the situation that the next-nearest-neighbor coupling χ\chi is weak than the nearest-neighbor coupling κ\kappa. The nontrivial topology supports helical and antihelical edge states [Fig. 1(h)] that are different from their band structures and their ways of propagation. The antihelical (helical) edge states on the parallel boundaries propagate in the same (opposite) direction, being robust against disorder because of the topological protection and the spatial separation between the edge and bulk states FranzPRL18 .

Refer to caption
Figure 2: (a) Spin-Chern numbers of topological phases indicated by nontrivial topologies of links. Topological phase transition occurs at the red plane. (b) Types of edge states. The entire parameter space J1J_{1}-J2J_{2}-δ\delta is divided into three regions. The cyan regions J1J2<0J_{1}J_{2}<0 are the insulating phase, the red regions δ2<J1J2\delta^{2}<J_{1}J_{2} are the metallic phase with the separable bands being present only in the horizontal direction, and the orange regions 0<J1J2<δ20<J_{1}J_{2}<\delta^{2} are the metallic phase with the separable bands being present in both directions. (c) Propagation of antihelical and helical edge states for pseudospin-up excitation; the black arrows indicate the propagation in the bulk, and they pertain to the topological phase 0<χ<κ0<\chi<\kappa. (d)-(g) Open boundary spectra of the pseudospin-up for the four cases on the left-half of (c). (d) J2=1,δ=1/2J_{2}=-1,\delta=1/2, (e) J2=1,δ=1/2J_{2}=1,\delta=1/2, (f) J2=1,δ=3/2J_{2}=1,\delta=3/2, (g) J2=1,δ=3/2J_{2}=1,\delta=-3/2. The other parameters are κ=1\kappa=1, J1=1J_{1}=1, and χ=1/2\chi=1/2. The simultaneous mirror reflection of the band structures with respect to kx=0k_{x}=0 and ky=0k_{y}=0 in (d)-(g) are the four cases on the right-half of (c). (h)-(k) Robust transports of antihelical edge states performed in a 40×4040\times 40 square lattice for (f) and (g) at the interface between two topological phases C=±1C_{\uparrow}=\pm 1, where the velocities along the top and bottom boundaries in the horizontal direction are the same. Dynamics are robust to random coupling disorder (see Supplementary materials C). (h) and (i) pertain to pseudospin-up excitation; (j) and (k) pertain to pseudospin-down excitation. (h) and (j) are for case (f), (i) and (k) are for case (g).

In Fig. 1(b), the regions J1J2<0J_{1}J_{2}<0 are the topological insulating phase; the helical edge states with different pseudospins propagate clockwise or counterclockwise along the lattice boundaries. The topological phase undergoes a metal-insulator phase transition at J1J2=0J_{1}J_{2}=0. The regions J1J2>0J_{1}J_{2}>0 are the topological metallic phase KamenevPRL18 ; the edge states excited by the pseudospin-up propagate toward the left (or right) on the two parallel boundaries, but the edge states excited by the pseudospin-down propagate toward the right (or left) boundaries for J1,J2>0J_{1},J_{2}>0 (or J1,J2<0J_{1},J_{2}<0). These edge states associated with the two pseudospins constitute the antihelical edge states, and the required counterpropagating modes for the antihelical edge states are the bulk states. The bands touch at the topological phase transition, where topological edge states are degenerate and exhibit identical dispersion and propagation velocity.

The square lattice exhibits nontrivial topology throughout the parameter space J1J_{1}-J2J_{2}-δ\delta [Fig. 2(a)]. The various types of edge states, as distinguished by their propagation, are presented in Fig. 2(b). The helical edge states are in J1J2<0J_{1}J_{2}<0 indicated in cyan and the antihelical edge states are in the other regions J1J2>0J_{1}J_{2}>0. The antihelical edge states differ from the helical edge states in term of their copropagation on the parallel boundaries MHNP11 ; FranzPRL18 . The antihelical edge state excitations with different pseudospins unidirectionally propagate in the opposite directions along the lattice boundaries and enable the separation of the robustly copropagated chiral mode. The helical edge states in the topological insulating phase always appear in both the horizontal and vertical directions; however, the antihelical edge states in the topological metallic phase usually appear only in one direction because of the inseparable band energy in the other direction. Here, the antihelical edge states are simultaneously present in both the horizontal and vertical directions for δ2>J1J2>0\delta^{2}>J_{1}J_{2}>0, where the energy bands are separable in both directions. The antihelical edge states are only present in the horizontal directions for δ2<J1J2\delta^{2}<J_{1}J_{2}. In Fig. 2(b), the surfaces J1J2=0J_{1}J_{2}=0 and δ2=J1J2\delta^{2}=J_{1}J_{2} divide the parameter space J1J_{1}-J2J_{2}-δ\delta into three phases with eight regions. The eight types of edge state propagation for the pseudospin-up excitations are presented in Fig. 2(c). The pseudospin-down edge state excitations propagate in the opposite directions.

The simultaneous sign changes of J1J_{1} and J2J_{2} reverse the propagation directions of all the antihelical edge states. The four cases of edge modes on the left-half and right-half of Fig. 2(c) propagate in opposite directions. The sign change of δ\delta only alters the propagation direction of the antihelical edge states in the vertical direction. The antihelical edge states at the strong χ1\chi_{1} (δ>0\delta>0) pass along the bottom-left and top-right corners and scatter into the bulk at the top-left and bottom-right corners; by contrast, the antihelical edge states at the strong χ2\chi_{2} (δ<0\delta<0) pass along the bottom-right and top-left corners and scatter into the bulk at the top-right and bottom-left corners. In Fig. 2(c), the helical edge states propagate counterclockwise (clockwise) along the lattice boundaries for C=+1C_{\uparrow}=+1 (C=1C_{\uparrow}=-1) in the cyan region J1>0,J2<0J_{1}>0,J_{2}<0 (J1<0,J2>0J_{1}<0,J_{2}>0). The antihelical edge states along the horizontal direction in the red region J1,J2>0J_{1},J_{2}>0 (J1,J2<0J_{1},J_{2}<0) of δ2<J1J2\delta^{2}<J_{1}J_{2}; the edge state excitation propagates leftward (rightward), scatters into the bulk at the corners, and goes backward. The antihelical edge states are present along both the horizontal and vertical directions when 0<J1J2<δ20<J_{1}J_{2}<\delta^{2}; the edge state excitation propagates half a closed-loop along the lattice boundaries, scatters into the bulk at the corners, and goes backward along the diagonal direction with the support of the scattering states. The four cases are respectively distributed in the orange regions (J1,J2,δ>0)\left(J_{1},J_{2},\delta>0\right), (J1,J2<0,δ>0)\left(J_{1},J_{2}<0,\delta>0\right), (J1,J2>0,δ<0)\left(J_{1},J_{2}>0,\delta<0\right), and (J1,J2,δ<0)\left(J_{1},J_{2},\delta<0\right) of 0<J1J2<δ20<J_{1}J_{2}<\delta^{2}.

Figures 2(d)-2(g) present the four cases of the band structures for the pseudospin-up presented on the left-half of Fig. 2(c). The propagation of antihelical edge states in Figs. 2(f) and 2(g) for the pseudospin-up are simulated in Figs. 2(h) and 2(i), and the corresponding propagation for the pseudospin-down are simulated in Figs. 2(j) and 2(k). The separation and convergence of different pseudospins toward the opposite corners of the square lattice are observed.

In conclusion, we have introduced the antihelical edge states in the time-reversal symmetric topological metallic phase, and proposed the photonic realization in coupled ring resonators. The antihelical edge states with different pseudospins propagate in the opposite directions along the corners toward the diagonal of the square lattice. Unconventional copropagation is applicable for the robust spin purification. The antihelical edge states become the helical edge states after undergoing a metal-insulator phase transition. The wide range of reconfigurable robust light propagation enables the flexible control of light flow at the edges. Our findings provide insight into the photonic topological metals and are applicable for the acoustic lattices and other two-dimensional metamaterials. Concepts pertaining to topological metals and antihelical edge states are inspiring in the condensed matters and topological materials.

Acknowledgments

This work was supported by National Natural Science Foundation of China (Grants No. 12222504, No. 11975128, and No. 11874225).

Supplementary Materials

.1 A. Experimental realization of the square lattice

In this section, the experimental realization of the square lattice in the two-dimensional (2D) coupled resonator array is discussed. Figure 3(a) is the schematic of the 2D coupled resonator array. The ring resonators are the primary resonators for the sites of the square lattice. The resonators in green stand for the sites AA and the resonators in orange stand for the sites BB. The linking resonators mediate the photons tunneling among the primary resonators and induce the effective couplings between the nearest-neighbor primary resonators and the next-nearest-neighbor primary resonators. The linking resonators and the primary resonators are coupled through their evanescent fields and the coupling strengths depend on the positions between the linking resonators and the primary resonators. For example, the couplings for the CW mode (pseudospin-up) of the primary resonators are mediated by the CCW/CW mode of the linking resonators.

Refer to caption
Figure 3: Schematic of the 2D square lattice of coupled ring resonators. The ring resonators in green (orange) are the sites AA (BB), and the other stadium resonators are the linking resonators. (a) χ1χ20\chi_{1}\chi_{2}\neq 0. (b) χ2=0\chi_{2}=0. The arrows in the ring resonators indicate the CW mode.

The vertical coupling κ\kappa is reciprocal, the path lengths for the photons tunneling upward and downward between the neighbor ring resonators are equal. The coupling strength is approximately characterized by κ=κl2/Δl\kappa=\kappa_{l}^{2}/\Delta_{l} LYouPRA17 , where κl\kappa_{l} and Δl=ωcωlink\Delta_{l}=\omega_{c}-\omega_{link} are the hopping and detuning between the primary resonators and the linking resonators; and the frequency of the primary (linking) resonators is ωc\omega_{c} (ωlink\omega_{link}).

The horizontal couplings iJ1iJ_{1} and iJ2iJ_{2} are nonreciprocal, carrying the Peierls phase eiπ/2e^{i\pi/2} in the couplings. The Peierls phase is implemented through the optical path length difference for the photons tunneling rightward and leftward between the neighbor ring resonators. The CW mode photons tunneling rightward from the lower half of the linking resonator experience an additional path length l=λ/2l=\lambda/2 than the CW mode photons tunneling leftward from the upper half of the linking resonator MHNP11 , where the wave length is λ\lambda. Therefore, the photons tunneling rightward acquire an extra phase factor eiϕ=eiπl/λ=ie^{i\phi}=e^{i\pi l/\lambda}=i and the photons tunneling leftward acquire an extra phase factor eiϕ=eiπl/λ=ie^{-i\phi}=e^{-i\pi l/\lambda}=-i in the front of the horizontal couplings iJ1iJ_{1} and iJ2iJ_{2}.

The cross coupling χ1\chi_{1} (χ2\chi_{2}) between the primary resonators on the diagonal of the square plaquette is directly (indirectly) mediated by the CCW (CW) mode of the linking resonator along the main diagonal of the square plaquette as indicated by the blue (red) arrows. The main diagonal of the square plaquette refers to the line along the upper-left and the lower-right corners of the square plaquette. The cross couplings χ1\chi_{1} and χ2\chi_{2} are independently mediated through the linking resonators along the diagonals of the square plaquette LeykamPRL18 .

The antihelical edge states can appear in the topological metal phase of the square lattice at the specific case χ1χ2=0\chi_{1}\chi_{2}=0. Thus, the single cross coupling case is adequate for the observation of antihelical edge states in experiments. A concrete system is χ=δ\chi=\delta. In this situation, one of the two cross couplings χ2=0\chi_{2}=0 vanishes as schematically illustrated in Fig. 3(b). This simplifies the setup in the experiment and facilitates the realization of antihelical edge states. The robust propagation with strong χ2\chi_{2} can be observed in the situation χ1=0\chi_{1}=0 through switching the orientation of the linking resonator to alter the connection between the nearest-neighbor resonators on the diagonals of the square plaquette. In this manner, the proposed square lattice can alter between two simple situations χ2=0\chi_{2}=0 and χ1=0\chi_{1}=0.

For example, we consider a simple configuration in experiment that the resonators are neatly arranged in both horizonal and vertical directions of the 2D square lattice. Consequently, the horizontal couplings have equal strengths J1=J2J_{1}=J_{2}. A candidate platform for the possible realization of the photonic topological metal and the antihelical edge states can be chosen as follows. In a situation κ=1\kappa=1, χ=δ=1/2\chi=\delta=1/2, J1=J2=1/4J_{1}=J_{2}=1/4, the two cross couplings are χ1=1\chi_{1}=1, χ2=0\chi_{2}=0 and the cross coupling strength equals to the vertical coupling strength χ1=κ\chi_{1}=\kappa. The round trip length of the resonator is about 7070 μ\mathrm{\mu}m. The resonator supports a single mode transverse electric field at the telecom wavelength 1.551.55 μ\mathrm{\mu}HafeziNP13 . The coupling strengths evanescently decay as the width of the air gap between the neighbor resonators; and the coupling strengths approximately decay from 30\sim 30 GHz\mathrm{GHz} to 5\sim 5 GHz\mathrm{GHz} for the air gap width increasing from 150150 nm\mathrm{nm} to 250250 nm\mathrm{nm} HafeZiNP19 . The topologically robust transport can be experimentally implemented in a 20×2020\times 20 size square lattice of coupled resonators at the vertical coupling strength κ20\kappa\sim 20 GHz\mathrm{GHz}, the cross coupling strength χ120\chi_{1}\sim 20 GHz\mathrm{GHz}, χ2=0\chi_{2}=0, and the horizontal coupling strengths J15J_{1}\sim 5 GHz\mathrm{GHz}, J25J_{2}\sim 5 GHz\mathrm{GHz}. In this situation, the 2D square lattice has a nontrivial topology because of χ=(χ1+χ2)/2<κ\chi=(\chi_{1}+\chi_{2})/2<\kappa; and the system is a topological metal hosting the antihelical edge states in both horizontal and vertical directions because of δ2>J1J2>0\delta^{2}>J_{1}J_{2}>0.

Refer to caption
Figure 4: Berry curvature in the BZ for h(𝐤)h_{\uparrow}\left(\mathbf{k}\right). The spin-Chern number CC_{\uparrow} as a summation of the Berry curvature in the entire BZ yields C=+1C_{\uparrow}=+1 for (a) J1=1,J2=1/2J_{1}=1,J_{2}=1/2, and C=1C_{\uparrow}=-1 for (b) J1=1,J2=3/2J_{1}=1,J_{2}=3/2. Other parameters are κ=1,χ=1/2\kappa=1,\chi=1/2, and δ=1\delta=1.

.2 B. Spin-Chern number characterized by the linking number of the effective magnetic field

Refer to caption
Figure 5: Energy bands for random disorder in all the couplings inside and outside the central one-half area of the square lattice. The upper panels indicate the locations of disorder shaded in gray, and the lower panels are the energy bands for the 1D projection lattices. (a) and (c) are for the disorder in the bulk. (b) and (d) are for the disorder on the boundaries. All the couplings with random disorder are deviated from the set parameters within the range of [20%,20%][-20\%,20\%]. The parameters are χ=κ/2=1/2\chi=\kappa/2=1/2, J1=J2=1J_{1}=J_{2}=1, and δ=3/2\delta=3/2. The square lattice size is 40×4040\times 40.
Refer to caption
Figure 6: Numerical simulations for the robust propagation in the square lattice for random disorder in all the couplings deviated from the set parameters within the range of [10%,10%][-10\%,10\%]. The parameters are χ=κ/2=1/2\chi=\kappa/2=1/2, J1=J2=1J_{1}=J_{2}=1, and δ=3/2\delta=3/2. (a) and (b) are the robust propagations of antihelical edge states for the pseudospin-up in contrast to the robust propagations in the main text Fig. 2(h); (c) and (d) are the robust propagations of antihelical edge states for the pseudospin-down in contrast to the robust propagations in the main text Fig. 2(j). The square lattice size is 40×4040\times 40. The unit of time tt is 1/κ1/\kappa.

The Chern number characterizes the band topology of the 2D topological phase and the Chern number is proportional to the Hall conductance. In a 2D time-reversal invariant system, although the total Chern number is zero, the spin-Chern numbers are quantized. The spin-Chern number of the lower band for the pseudospin-up (CW) mode is defined by

C=12πBZΩ𝑑kx𝑑ky,C_{\uparrow}=\frac{1}{2\pi}\iint_{\mathrm{BZ}}\Omega_{-}dk_{x}dk_{y}, (1)

where Ω=𝐤×𝒜𝐤\Omega_{-}=\nabla_{\mathbf{k}}\times\mathcal{A}_{{}_{\mathbf{k}}} is the Berry curvature, 𝒜𝐤=iΨ(𝐤)||Ψ(𝐤)\mathcal{A}_{{}_{\mathbf{k}}}=-i\left\langle\Psi_{-}(\mathbf{k})\right|\nabla\left|\Psi_{-}(\mathbf{k})\right\rangle is the Berry connection, and |Ψ(𝐤)\left|\Psi_{-}(\mathbf{\ k})\right\rangle is the eigenstate of the lower band WGa15 . The spin-Chern number is an integral of the Berry curvature in the entire Brillouin zone (BZ). Figure 4 provides the numerical results of the Berry curvature in the BZ for the topological phases C=+1C_{\uparrow}=+1 and C=1C_{\uparrow}=-1 TFuku05 .

In the following, we explain the relation between the spin-Chern number and the linking number. The spin-Chern number for the two-band Hamiltonian h(𝐤)=d0σ0+𝐝(𝐤)σh_{\uparrow}\left(\mathbf{k}\right)=d_{0}\sigma_{0}+\mathbf{d(k)\cdot\sigma} for the pseudospin-up is expressed in the form of QWZPRB06

C=14π𝐝(kx𝐝×ky𝐝)|𝐝|3𝑑kx𝑑ky.C_{\uparrow}=-\frac{1}{4\pi}\int\frac{\mathbf{d}\cdot\left(\partial_{k_{x}}\mathbf{d}\times\partial_{k_{y}}\mathbf{d}\right)}{\left|\mathbf{d}\right|^{3}}dk_{x}dk_{y}. (2)

The effective magnetic field in h(𝐤)h_{\uparrow}\left(\mathbf{k}\right) is 𝐝(𝐤)=𝐫1(kx)𝐫2(ky)\mathbf{d(k)}=\mathbf{r}_{1}(k_{x})-\mathbf{r}_{2}(k_{y}) with 𝐫1(kx)=(κ+2χcoskx,2δsinkx,(J1J2)sinkx)\mathbf{r}_{1}(k_{x})=(\kappa+2\chi\cos k_{x},-2\delta\sin k_{x},(J_{1}-J_{2})\sin k_{x}) and 𝐫2(ky)=(κcosky,κsinky,0)\mathbf{r}_{2}(k_{y})=(-\kappa\cos k_{y},\kappa\sin k_{y},0); and d0=(J1+J2)sinkxd_{0}=(J_{1}+J_{2})\sin k_{x}. Substituting 𝐝(𝐤)=𝐫1(kx)𝐫2(ky)\mathbf{d(k)}=\mathbf{r}_{1}(k_{x})-\mathbf{r}_{2}(k_{y}) into equation (2), the spin-Chern number is rewritten in the form of

C=14π02π02π𝐫1𝐫2|𝐫1𝐫2|3(𝐫1kx×𝐫2ky)𝑑kx𝑑ky.C_{\uparrow}=\frac{1}{4\pi}\int\nolimits_{0}^{2\pi}\int\nolimits_{0}^{2\pi}\frac{\mathbf{r}_{1}-\mathbf{r}_{2}}{\left|\mathbf{r}_{1}-\mathbf{r}_{2}\right|^{3}}\cdot\left(\frac{\partial\mathbf{r}_{1}}{\partial k_{x}}\times\frac{\partial\mathbf{r}_{2}}{\partial k_{y}}\right)dk_{x}dk_{y}. (3)

In geometry, equation (3) is the definition of the linking number of two independent closed curves 𝐫1(kx)\mathbf{r}_{1}(k_{x}) and 𝐫2(ky)\mathbf{r}_{2}(k_{y}). The linking number is a topological invariant that characterizes the number of times that 𝐫1(kx)\mathbf{r}_{1}(k_{x}) and 𝐫2(ky)\mathbf{r}_{2}(k_{y}) wrap around each other Hannay05 . Thus, the spin-Chern number is equivalent to the linking number of two closed curves 𝐫1(kx)\mathbf{r}_{1}(k_{x}) and 𝐫2(ky)\mathbf{r}_{2}(k_{y}) of the effective magnetic field 𝐝(𝐤)\mathbf{d}\left(\mathbf{k}\right) for the pseudospin-up Hamiltonian h(𝐤)h_{\uparrow}\left(\mathbf{k}\right). The spin-Chern numbers shown in Fig. 4 are in accords with the representative links shown in Fig. 2(a) of the main text.

For the pseudospin-down mode, the signs of the horizontal couplings J1J_{1} and J2J_{2} change into the opposite J1-J_{1} and J2-J_{2}. The Hamiltonian reads h(𝐤)=d0σ0+(𝐝x,𝐝y,𝐝z)σh_{\downarrow}\left(\mathbf{k}\right)=-d_{0}\sigma_{0}+\left(\mathbf{d}_{x}\mathbf{,d}_{y}\mathbf{,}-\mathbf{d}_{z}\right)\mathbf{\cdot\sigma}. The spin-Chern number obtained from the linking of the effective magnetic field (𝐝x,𝐝y,𝐝z)\left(\mathbf{d}_{x}\mathbf{,d}_{y}\mathbf{,}-\mathbf{d}_{z}\right) yields C=CC_{\downarrow}=-C_{\uparrow}. This is straightforward because that the curve 𝐫2(ky)\mathbf{r}_{2}(k_{y}) and two components of the curve 𝐫1(kx)\mathbf{r}_{1}(k_{x}) in the xx and yy directions are unchanged, and the component of the curve 𝐫1(kx)\mathbf{r}_{1}(k_{x}) in the zz direction changes into the opposite.

.3 C. Robustness of the antihelical edge states

Refer to caption
Figure 7: Numerical simulations for the robust propagation in the square lattice of size 80×8080\times 80. The parameters are χ=κ/2=1/2\chi=\kappa/2=1/2, J1=J2=1J_{1}=J_{2}=1, and δ=3/2\delta=3/2, being identical with the parameters in Fig. 6. The unit of time tt is 1/κ1/\kappa.

In the topological phases of the square lattice, both the helical edge states and the antihelical edge states are topologically-protected and robust to disorder. In comparison, the helical edge states in their propagations are even more stable to the disorder because of the band gap protection. In this section, we provide more details for the robustness of antihelical edge states and their unidirectional propagations in the presence of coupling disorder; notably, the couplings with random disorder satisfy the particle-hole symmetry, which also exists in the square lattice for the pseudospin-up (pseudospin-down). The numerical simulations of the band structures and the robust propagations of antihelical edge states are shown for the imperfect square lattice. The localized edge states and the extended bulk states are spatially separated. Thus, the imperfection in the bulk almost does not affect the edge states and their robust propagations.

The distribution and the localization of the edge states on the boundaries are insensitive to the lattice size; however, the distribution and the extended feature of the bulk states closely depend on the lattice size. The larger lattice size leads to a better spatial separation between the edge states and the bulk states. The influence on the edge states for the disorder in the bulk of the square lattice is slight in comparison with the influence on the edge states for the disorder on the boundaries of the square lattice. This point is elucidated in Fig. 5, where we depict the energy bands for the disorder in the bulk and on the boundaries, respectively; and the random disorder is chosen on all the couplings either inside or outside the central one-half area of the square lattice as schematically illustrated in the upper panels of Fig. 5. The edge states are depicted in Figs. 5(a) and 5(c) for the random disorder in the couplings inside the central one-half area of the square lattice in the bulk; and the edge states are depicted in Figs. 5(b) and 5(d) for the random disorder in the couplings outside the central one-half area of the square lattice on the boundaries. To obtain the 1D projection lattices, the 2D square lattice is set translationally invariant in the xx (yy) direction for Figs. 5(a) and 5(b) [Figs. 5(c) and 5(d)].

In the numerical simulations, the initial excitation has the Gaussian profile in the form of

|Ψ(t0)=Ω1/2kεe(kϵk0)2/(2α2)eiNc(kϵk0)|ψ,\left|\Psi(t_{0})\right\rangle=\Omega^{-1/2}\sum\nolimits_{k_{\varepsilon}}e^{-(k_{\epsilon}-k_{0})^{2}/(2\alpha^{2})}e^{-iN_{c}(k_{\epsilon}-k_{0})}\left|\psi\right\rangle, (4)

where |ψ\left|\psi\right\rangle is the edge mode with the momentum kϵ=kxk_{\epsilon}=k_{x} or kyk_{y}, and k0=πk_{0}=\pi. NcN_{c} is the center of the wave packet and α\alpha controls the width of the Gaussian profile. Ω\Omega normalizes |Ψ(t0)\left|\Psi(t_{0})\right\rangle. The numerical simulations of the robust light propagation for the excitations of antihelical edge states in the presence of random disorder are performed and demonstrated in Fig. 6. All the couplings with random disorder are deviated from the set parameters of the square lattice within the range of [10%,10%][-10\%,10\%]. The edge mode excitations robustly propagate along the boundaries. The dynamics in the numerical simulations are close to the dynamics in the absence of disorder exhibited in the main text Fig. 2; this is a consequence of the limited lattice size 40×4040\times 40 in the numerical simulations. The robust propagation against disorder is more excellent in the numerical simulations performed in the larger size system as shown in Fig. 7.

In the topological metal phase possessing the antihelical edge states in both the horizontal and vertical directions, there are two types of right-angled (9090^{\circ} degrees) corners in the square lattice. One type of corners help the right-angled turning of the edge mode excitations when they propagate along the boundaries; and the other type of corners help the edge mode excitations entering the bulk of the square lattice. Each corner is associated with a cross coupling; and the two types of corners are distinguished from their associated cross couplings χ+δ\chi+\delta and χδ\chi-\delta. The formation of the right-angled corners associated with the weak cross coupling enables the 9090^{\circ} degrees turning across the corners for the edge mode excitations propagating along the lattice boundaries. The formation of the right-angled corners associated with the strong cross coupling enables the edge mode excitations entering the bulk of the square lattice.

References

  • (1) Lu L, Joannopoulos JD, Soljačić M. Topological photonics. Nat Photonics 2014; 8: 821-9.
  • (2) Khanikaev AB, Shvets G. Two-dimensional topological photonics. Nat Photonics 2017; 11: 763-73.
  • (3) Ozawa T, Price HM, Amo A, et al. Topological photonics. Rev Mod Phys 2019; 91: 015006.
  • (4) Haldane FDM, Raghu S. Possible Realization of Directional Optical Waveguides in Photonic Crystals with Broken Time-Reversal Symmetry. Phys Rev Lett 2008; 100: 013904.
  • (5) Wang Z, Chong Y, Joannopoulos JD, et al. Observation of unidirectional backscattering-immune topological electromagnetic states. Nature 2009; 461: 772-5.
  • (6) Hafezi M, Demler EA, Lukin MD, et al. Robust optical delay lines with topological protection. Nat Phys 2011; 7: 907-12.
  • (7) Fang K, Yu Z, Fan S. Realizing effective magnetic field for photons by controlling the phase of dynamic modulation. Nat Photonics 2012; 6: 782-7.
  • (8) Liang GQ, Chong YD. Optical Resonator Analog of a Two-Dimensional Topological Insulator. Phys Rev Lett 2013; 110: 203904.
  • (9) Yang HF, Xu J, Xiong ZF, et al. Optically Reconfigurable Spin-Valley Hall Effect of Light in Coupled Nonlinear Ring Resonator Lattice. Phys Rev Lett 2020; 127: 043904.
  • (10) Colomés E, Franz M. Antichiral Edge States in a Modified Haldane Nanoribbon. Phys Rev Lett 2018; 120: 086603.
  • (11) Zhou P, Liu GG, Yang Y, et al. Observation of Photonic Antichiral Edge States. Phys Rev Lett 2020; 125: 263603.
  • (12) Yang Y, Zhu D, Hang Z, et al. Observation of antichiral edge states in a circuit lattice. Sci China Phys Mech Astron 2021; 64: 257011.
  • (13) Chen J, Li ZY. Prediction and Observation of Robust One-Way Bulk States in a Gyromagnetic Photonic Crystal. Phys Rev Lett 2022; 128: 257401.
  • (14) Leykam D, Mittal S, Hafezi M, et al. Reconfigurable Topological Phases in Next-Nearest-Neighbor Coupled Resonator Lattices. Phys Rev Lett 2018; 121: 023901.
  • (15) Ying X, Kamenev A. Symmetry-Protected Topological Metals. Phys Rev Lett 2018; 121: 086810.
  • (16) Yang F, Liu YC, You L. Anti-𝒫𝒯\mathcal{PT} symmetry in dissipatively coupled optical systems. Phys Rev A 2017; 96: 053845.
  • (17) Hafezi M, Mittal S, Fan J, et al. Imaging topological edge states in silicon photonics. Nat Photonics 2013; 7: 1001-5.
  • (18) Mittal S, Orre VV, Zhu G, et al. Photonic quadrupole topological phases. Nat Photonics 2019; 13: 692-6.
  • (19) Gao W, Lawrence M, Yang B, et al. Topological Photonic Phase in Chiral Hyperbolic Metamaterials. Phys Rev Lett 2015; 114: 037402.
  • (20) Fukui T, Hatsugai Y, Suzuki H. Chern Numbers in Discretized Brillouin Zone: Efficient Method of Computing (Spin) Hall Conductances, J Phys Soc Jpn 2015; 74: 1674-7.
  • (21) Qi XL, Wu YS, Zhang SC. Topological quantization of the spin Hall effect in two-dimensional paramagnetic semiconductors. Phys Rev B 2006; 74: 085308.
  • (22) Dennis MR, Hannay JH. Geometry of Călugăreanu’s theorem, Proc R Soc A 2005; 461: 3245-54.