Floquet topological photonic crystals with temporally modulated media

We show that Floquet topological insulating states can exist in two-dimensional photonic crystals made of time-variant optical materials. By arranging the modulating phases, it facilitates effective gauge fields that give rise to topological effects. The band structures demonstrate the existence of topologically non-trivial bandgaps, thereby leading to back-scattering immune unidirectional edge states owing to bulk-edge correspondence. With these first-principle numerical results, we then verify the topological order for every Floquet band via Wilson loop approach. In the final paragraph, the possible experimental implementation for Floquet topological photonics is also discussed. Published by The Optical Society under the terms of the Creative Commons Attribution 4.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal


Introduction
Over the past few decades, the research of temporally modulated systems has been discussed but yet to be fully explored [1][2][3][4][5][6][7][8][9]. Due to the instability of time-varying structures, the progress of wave physics was predominately made in the nature of time-independent systems. However, the recent studies toward this subject have indicated that a time-dependent system could be stable so long as the modulation satisfies time-periodicity. The finding has attracted great attention because the research into this nascent subject of physics is thus initiated. To date, numerous fascinating phenomena have been proposed in classical wave systems, such as non-reciprocal transports [10], effective gauge fields [11][12][13], and wave amplifications [14]. On the other hand, in condensed matter theory the material with periodically dynamic modulation could generate novel propagating modes which travel along the boundary unidirectionally [15] and this unique state originates from the band topology carrying non-zero Chern numbers. It is worth noting that as long as this topological index remains unchanged, the wave propagation is robust against backscattering from defects and impurities. Recently, this temporally modulating scheme has been applied to classical wave systems as well as condensed-matter regime [16]. Particularly in optical systems, photonic Floquet topological insulators have been proposed and experimentally verified in periodically arranged coupled waveguides [17].
Despite the successful realization of waveguide optics, generating Floquet topological insulating states in dynamically modulated optic media remains significantly challenging as such materials do not exist in nature. Therefore, in this article a feasible way to implement time-periodic perturbation and Floquet topological photonic crystals (FPhC) is proposed. Starting from investigating Floquet band structures, topological edge states are verified by numerical calculation of Chern numbers. Then the field profiles showing robustness against defects are demonstrated. Lastly, a Floquet topological metamaterial consisting of metallic wires with temporal modulations is discussed for a viable experimental implementation.

Floquet band structure
In Fig. 1(a), we begin by considering a two-dimensional (2D) FPhC with hexagonal lattice in the microwave range. Each perfect electric conductor (PEC) cylinder is surrounded by three dielectrics given by where ε m , δ, and Ω are dielectric constant, modulating strength, and modulating frequency, respectively. The indices j = 1,2,3 label different modulating phases. As indicated in Ref. [16], this phase arrangement between modulated dielectrics generates an effective magnetic field, which gives rise to topological bandgaps. The radius of PEC cylinder is 0.44a and the thickness of surrounding temporal dielectrics is 0.04a, where a denotes lattice constant. These high filling-ratio PEC cylinders are implemented to form four-band localized modes because these modes interact more efficiently with modulated media. For TE modes, we substitute E z = n E z,n e −i(ω+nΩ)t into the wave equation and obtain [16]  Throughout this article, we assume the permeability of all the materials is unity, whereas the matters with magnetic response should not significantly affect the result. Numerically, as the modulation strength is weak, the harmonic index n in Eq. (2) can be truncated from -2 to 2 that contributes the major composition of the total electric field. By employing FEM-based commercial software COMSOL 5.4, in Fig. 1(b) the blue dash lines represent the band structure in the absence of modulation and a four-band configuration that can be analyzed by a dipolar model is shown [18]. In this unmodulated band structure, there are two quadratic topological degeneracies at Γ point between 1 st and 2 nd bands, and 3 rd and 4 th bands, and one linear (Dirac) degeneracy at K point between 2 nd and 3 rd bands. Due to the spin-like nature of dipolar fields (p x ± ip y ), when a rotation-related field is applied, such as constant magnetic fields [19], Coriolis forces [20], or synthetic gauge fields [21], all the degeneracies are lifted owing to the time-reversal symmetry breaking. This can be seen through the modulated case with perturbation strength δ = 0.15 in Fig. 1(b). The red solid lines in Fig. 1(b) demonstrate a quasi-frequency diagram with four isolated bands. Due to the existence of the effective gauge field, these bands lead to topologically protected edge modes when a boundary is truncated.

Band topology analysis
To confirm the topological order for each band, here we have numerically calculated topological orders to further identify their band topology. Since there exists no coupling between quasifrequency structures, the topological states are therefore characterized by non-zero Chern numbers. The Chern number C n for every band n is defined by the closed line integral of Berry connections around the first Brillouin zone (BZ), [22,23] where k = (k 1 , k 2 ) is Bloch wave-vector, A n,k = −i E z,n ∂ k |E z,n are Berry connections for the n-th band, and θ n,k 1 = ∫ π −π dk 2 A n,k is the evolution of the Berry phase for nth band with fixed k 1 in the BZ. According to Wilson loop approach [23][24][25][26], Chern numbers are equivalent to winding numbers of Berry phase evolution, which can be calculated from the eigenvalues w n of the matrix W. This matrix W(k 1 ), which is defined by a product of the N×N Berry connection matrices along the Wilson loop, reads where |E z,n is Bloch function evaluated via COMSOL. With the matrix W(k 1 ), the corresponding eigenvalues w n,k 1 give rise to Berry phases, namely, We then present the evolution of Berry phases for an effective 1D system with k 1 in Fig. 2 and Chern numbers can be extracted in the following manner: within the region k 1 ∈ [−π, π], θ ∈ [−π, π], one can map the evolution line of Berry phases onto the BZ, thereby determining the topological invariant via the number of times the Berry phases wind around the BZ. Due to time-reversal symmetry breaking, the Berry phases of the four isolated bands in Fig. 2

Robust one-way edge modes
The emergence of topological boundary modes may be one of the most essential properties in topological materials. Based on the preceding analysis, topological one-way edge states are expected between the boundaries with distinct topology in accordance with bulk-edge correspondence. As the Chern numbers from the bottom to top bands in Fig. 1(b) are respectively given by -1,0,0,1, three corresponding bandgaps are topologically non-trivial, thereby supporting the edge states with unidirectional nature. To validate this argument, Fig. 3(a) depicts the projected band structure that presents topological edge states obtained by numerically calculating a 1×10 supercell lattice with truncation in x direction. As expected, one-way propagating modes arise within every bandgap, and their numbers obey bulk-edge correspondence respectively. In the following, we will focus on the middle topological edge modes since the second bandgap is relatively larger than the others. Excited by a line current source with frequency equals 3.0924, the field profile demonstrates that electromagnetic waves travel along the boundary unidirectionally and propagate robustly against the defect, as shown in Fig. 3(b). In addition, when the phase arrangements flip sign, i.e. (φ 1 , φ 2 , φ 3 ) = (0, 2π/3, 4π/3), the wave in Fig. 3(c) propagates along the counter-clockwise direction due to the presence of negative group velocity, but the topological nature remains as it propagates around a sharp corner without reflections.

Possible experimental implementation
In this section, the possible experimental realization is discussed. To implement the proposed temporal dielectrics, a metamaterial scheme in microwave region may be employed as this frequency range suits the use active electrical elements such as time-variant inductors [27]. Consider a ring of conducting thin wires circulating around a PEC cylinder, as shown in Fig. 4(a). Each wire, whose radius equals 0.015a, is connected with a time-variant inductor to introduce temporal modulation. With TE polarization, the current density J z in a lossless wire can be expressed as [27][28] d where L m = L 0 + L 1 [1 + δ cos(Ωt + φ j )]. L 0 and L 1 are the original inductor from wires and the added time-modulated inductor, respectively. α is a coefficient defined by ω 2 p (L 0 + L 1 ), where ω p denotes Drude plasma frequency. Also, the wave equation for electric field E z is Substituting J z = n J z,n e −i(ω+nΩ)t and E z = n E z,n e −i(ω+nΩ)t into Eqs. (6) and (7), two coupled matrix differential equations are then obtained. Again, by exploiting COMSOL, Fig. 4(b) depicts Floquet band structure which exists four isolated bands with Chern numbers equal to -1,0,0,1. Similarly, these topologically non-trivial bands lead to unidirectional edge modes once a domain wall is introduced.

Conclusion
In conclusion, Floquet topological insulating states in a time-modulated photonic crystal have been proposed. The work starts from a FPhC with PEC cylinders surrounded by three time-variant dielectrics. The Floquet band structure is numerically calculated by COMSOL, and the Chern number analysis further proves the distinct topology for each band. We then study the emergence of unidirectional edge states through introducing a 1×10 super-cell lattice. The projected band structure demonstrates that topological edge states arise within all the Floquet bandgaps. At frequency ω/2πc = 3.0924, the propagation exhibits topological protection since no reflection is generated when the waves hit a defect. Lastly, a possible implementation via a combination of time-variant inductors and metallic wires has been proposed. The coupling between the current density and electric field enables us to add dynamic modulations into the system, leading to Floquet topological insulating phases similar to the previous model.

Disclosures
The authors declare that there are no conflicts of interest related to this article.