Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Band structure and topological phase transition of photonic time crystals

Open Access Open Access

Abstract

We investigate the band structure and topological phase transition of photonic time crystals (PTC)-systems in which the physical parameter varies periodically in time. We find that the topological phase transition of the PTC system can be revealed by the wave vector gap (k-gap) size, which was induced by the temporal refraction and reflection. Interestingly, a temporal zero-averaged refractive index k-gap is obtained when the PTC system includes a dispersive medium. This special k-gap is invariant with modulation time scaling at a given modulation frequency.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Photonic time crystals (PTC), the physical parameters of which are periodically modulated in time, are temporal analogues of conventional photonic crystals [1–3]. A PTC system is seen as a promising platform for producing amplification and non-reciprocal behavior [4,5]. The breaking of time (space) transition symmetry corresponds to energy (momentum) conservation. That is, temporal (spatial) forward and backward propagating waves will strongly interfere with each other, and inducing a gap in wave vector (frequency) space when the system parameters have temporal (spatial) periodicity [6,7]. Thus, the PTC system can provide a new path for controlling electromagnetic (EM) waves.

In theory, PTC systems can realize all potential practical applications correspond to conventional photonic crystals. We know that the PTC system is required to occur on the time scale of the oscillation of the probe signal in order to observe its band structure features. Thus, setting up such a PTC platform is not easy, especially for high-frequency waves. In the past few years, many strategies are adopted to realize the PTC platforms in experiments [8–10]. It is exciting that a thin-film material with a permittivity that pulsates (uniformly in space) at optical frequencies is realized based on the high nonlinearity effect in epsilon-near-zero materials [10].

Topology in a new degree of freedom for describing an object and has drawn much research attention recently [11,12]. It is actually a mathematical concept with respect to quantities that are preserved under continuous deformations. Topology was extended to various physical systems, such as condensed matter, optics and acoustics, and many others [13–16]. Researchers utilize the topological invariant to characterize the topology of a system. Systems with the same topological invariants are topological equivalent. A topological phase transition occurs when the topological invariant of the system changes, which plays a key role in understanding topological edge states [17].

The topology of an object is robust against the presence of disorder or defects, thus inducing a sea of potential applications in optics, such as unidirectional lasing and photonic circuitry [18,19]. It is well-known that the topological invariant is described by the bulk band geometric phase in a photonic system, i.e., the first Chern number in a two-dimensional periodic system [20–23] and the Zak phase in a one-dimensional periodic system [24–26], which reveals the quantized collective behavior of the EM field on the bulk band. A rigorous relation between the surface states of one-dimensional photonic crystals and the Zak phases of the bulk bands was reported [27]. Inspired by that work, we introduce the topological phase into a PTC system.

The topological phases of a PTC system were proposed very recently [3], the authors calculated the topological invariant of bulk bands and obtained the localized interface state with two PTC of different topological invariants. In this work, we focus on the topological invariants of bulk bands, the k-gap size and their relations. We then reveal the conditions of the topological phase transition of the PTC system. In addition, the dispersive medium for the first time, to the best of our knowledge, is included in the PTC system, which yields a temporal zero-averaged refractive index k-gap that is insensitive to the modulation time scaling at a given modulation frequency.

2. Theoretical method

We consider EM wave propagation inside the PTC with time-dependent permittivity and constant permeability. The constitutive relations of the PTC can be written as D=ε0ϵ(t)E, B=μ0μrH(μr=1), where ε0 and μ0 are the permittivity and permeability in vacuum, respectively, and ϵ(t) is the modulation function of the time-dependent permittivity. We restrict ourselves to the case where ϵ(t)is a periodic step-like function:

ϵ(t)={ϵ1(nTtt1+nT)ϵ2(t1+nTt(n+1)T),
where T=t1+t2 is the modulation period. t1(t2) is the duration correspond to the permittivity ϵ1(ϵ2), and n is an integer. In practice, experimenters typically use pump light to drive the materials for realizing the time-dependent permittivity (PTC system) and use a probe signal to study the properties of the PTC system [10]. The step-like permittivity is easier to realize in experiments than continuously varying one.

Using the fundamental solution to the wave equation and Eq. (1), the electric field of EM wave takes the form

Ei(t)=A0ei(ωit+kr)+B0ei(ωitkr),
where A0 and B0 are amplitudes, and k and ωi=kc/ϵi(i=1,2) denote the wave vector and angular frequency of the EM wave, respectively. One should note that the angular frequency of the EM wave varies with time while the wave vector remains constant. In other words, the EM wave inside the PTC system obeys conservation of momentum at the temporal boundary. Thus, the distinction between the PTC and conventional photonic crystals is the independent variable of the optical parameters, which is responsible for a wide range of exotic phenomena.

We assume the variation of permittivity occurs at the temporal boundary t0. Integrating the Maxwell equations from t0 to t0+ yields

t0t0+(×E)dt=t0t0+Btdt=[B(t0+)B(t0)]
t0t0+(×H)dt=t0t0+Dtdt=D(t0+)D(t0)
The left-hand side of Eq. (3) will vanish in the limit, t0+t0=t0t0=Δt0. Thus, the continuity conditions of the EM fields at the temporal boundary should be
B|t0=B|t0+
and
D|t0=D|t0+.
Combining the Bloch theorem in time, one can obtain the following dispersion relation for the PTC:
cos(ΩT)=cosω1t1cosω2t212(z1z2+z2z1)sinω1t1sinω2t2,
where Ω is the “Bloch frequency”, akin to the “Bloch wave vector” in the case of spatial periodicity, and zi=μi/ϵi(i=1,2). When the absolute value of the right-hand side of Eq. (7) is smaller than or equal to unity, the corresponding Bloch frequencies are real numbers, which means the PTC supports these modes. In contrast, a value larger than unity, i.e., the Bloch frequency is nonzero and imaginary, corresponds to an evanescent mode in PTC.

Generally, the Zak phase is used to characterize the topological invariant of a one-dimensional photonic system. Consulting the definition of the Zak phase for conventional photonic crystals [27], the Zak phase of the bulk band in a PTC system can be written as

θmZak=π/Tπ/T[iunitcelldtε(t)um,Ω*(t)Ωum,Ω(t)]dΩ,
where iunitcelldtε(t)um,Ω*(t)Ωum,Ω(t) is the Berry connection and um,Ω(t) is the unit cell Bloch eigenstate on the mth bulk band with Bloch frequency Ω, which can be obtained using the transfer matrix method.

3. Results and discussion

The PTC system we are about to consider is shown in Fig. 1, where T=t1+t2 is the unit cell size. Here, A and B correspond to the two states of the PTC at time intervals t1 and t2 with refractive indices n1=ϵ1μ1 and n2=ϵ2μ2, respectively. The unit cell with inversion symmetry is marked by the dashed line in Fig. 1. We assume the probe signal is a plane wave, which will be refracted and reflected in the time domain. In the following calculations, the wave vector of the probe signal k is normalized by k0=ω0/c(ω0=1GHz). Moreover, we employ two ancillary parameters g=n1/n2 and f=t1/T, which correspond to the modulation strength and modulation frequency, respectively. According to the definition, zero (infinite) modulation strength means n1<<n2(n1>>n2), thus the bandgap is infinite. In addition, when the modulation strength is unity, the system medium is uniform rather than time dependent. It is foreseeable that the properties of the band gap exhibit a certain symmetry with respect to g=1. Here, we focus on modulation strength ranging from 0.1 to 0.9. Similarly, a zero (unity) modulation frequency means t1<<t2(t1>>t2), and the bandgap will disappear.

 figure: Fig. 1

Fig. 1 Schematic diagram of the PTC system.

Download Full Size | PDF

3.1 Band structure of the PTC

In the limit of weak dispersion, we calculate the band structure of the PTC using Eq. (7), and the result is shown in Fig. 2(a). The basic physical parameters are ϵ1=4, ϵ2=9 and μ1=μ2=1, and the modulation frequency is f=0.3. The calculated result indicates that the band structure exhibits bandgaps in wave vector space (k-gap), which corresponds to a Bloch frequency with nonzero imaginary part. This results from scattering of the probe signal. Compared with previous research that treats permittivity as a continuous function of time [2,7], the size of the higher-order k-gap here is no longer close to zero.

 figure: Fig. 2

Fig. 2 Band structure of the PTC with different modulation frequencies: (a) f=0.3, (b) f=0.4, and (c) f=0.5. (d) The displacement field distribution of the compound PTC system constructed by (PTC1)10 and (PTC2)10. The subsystems PTC1 and PTC2 correspond to the systems (a) and (c), respectively.

Download Full Size | PDF

For convenience, we marked the bulk band 2, k-gap 2, and bulk band 3 with blue letters in Fig. 2. At the same time, we calculate the Zak phases of bulk bands 2 and 3 using Eq. (8), which is shown in red letters. The Zak phase reveals the topological properties of one bulk band and can be used to predict the existence of a localized state [28].

In Figs. 2(b) and 2(c), we examine the band structure when the modulation frequencies of the PTC system are f=0.4 and f=0.5. The remaining physical parameters are the same as those used in Fig. 2(a). One can clearly see that the bulk bands 2 and 3 cross when the modulation frequency increases from 0.3 to 0.4, i.e., k-gap 2 is closed. When the modulation frequency is increases further, k-gap 2 opens again and is accompanied by a change of the Zak phase in bulk bands 2 and 3, which represents a topological phase transition in the PTC system. Therefore, the band structure and bulk band topological phase can be effectively tuned by choosing the appropriate modulation frequency. This makes PTC particularly interesting. Moreover, we also calculate the displacement field distribution of the compound PTC system by transfer matrix method, which is shown in Fig. 2(d). The compound PTC system constructed by two subsystems (PTC1 and PTC2) that undergo the topological phase transition. We find that a localized state forms at the interface between PTC1 (10 periods) and PTC2 (10 periods), and the localized interface state exponential decay on both sides of the interface. In fact, this result can be seen as an analogue of Su-Schrieffer-Heeger model [29] in PTC system.

3.2 The condition of topological phase transition in PTC system

One of the important features of a band structure is their gap size. This can be used to optimize the structure for specific applications, such as improving the properties of microwave antenna, vertical-cavity surface-emitting lasers, and others [30,31]. We find that it also reveals information on the topological phase transition of the PTC system. In the following, we focus our attention on the second k-gap size and investigate the effects of changes in modulation frequency f and modulation strength g in detail.

Figure 3(a) shows a contour map of the second k-gap size Δk in unit k0. The horizontal and vertical axes represent the modulation frequency f and modulation strength g, respectively. We take the basic parameter ϵ2=4, followed by ϵ1=gϵ2. One can clearly find that the second k-gap size varies obviously as the modulation parameters change. Thus, k-gap size can be designed in a predictable manner. In addition, when the modulation strength and modulation frequency satisfy g=f/(1f) (the dashed line in Fig. 3), the second k-gap size is equal to zero, i.e., a topological phase transition occurs. Note that this constraint condition g=f/(1f) corresponds to n1/t1=n2/t2, where ni/ti has the physical meaning of optical thickness in PTC system. In other words, when the states A and B have the same optical thickness, the second k-gap is closed. This is consistent with the conclusion of the conventional photonic crystal system. Therefore, bulk bands 2 and 3 have different topological phases on each side of the dashed line. According to this discovery, when the modulation frequency f=0.4, the occurrence of topological phase transition requires modulation strength g=0.67, which corresponds to the case shown in Fig. 2(b) exactly.

 figure: Fig. 3

Fig. 3 Variation of the second k-gap size Δk with modulation frequency f and modulation strength g: (a) ϵ2=4 and (b) ϵ2=9.

Download Full Size | PDF

Moreover, we also calculate the second k-gap size for ϵ2=9, and the result is shown in Fig. 3(b). In this case, the second k-gap size increases given the same modulation parameters used in Fig. 3(a). Interestingly, the relations between the second k-gap size and the modulation parameters (f and g) remain unchanged, which demonstrates the universality of the results obtained in Fig. 3(a).

We find the topological phase transition will be more complicated for higher order bandgaps. This begs the question: what is the relation between them? In what follows, we move our attention to explaining this phenomenon.

We know that the absolute value of the right-hand side of Eq. (7) must be smaller than or equal to unity when sin(kct2/ϵ2)=0. That is to say the wave vector k=m2πϵ2/ct2(m2N+) must be in the bulk band. Then, we have cos(kct2/ϵ2)=(1)m2, cos(kct1/ϵ1)=(1)m1(m1N+). Thus, cos(ΩT)=(1)m1+m2. Namely, ΩT=0 when m1+m2 is even, and ΩT=±π when m1+m2 is odd. The topological phase transition is equivalent to the degeneracy of band edges. If m=m1+m2 is the order of the k-gap, its degeneracy is just the times of topological phase transitions in the full parameter space. Since m1 and m2 are positive integer numbers, the topological phase transition time of the mth k-gap is m1.

The results above reveal the condition of topological phase transition in PTC system. However, how to observe this topological phase transition in experiment? In conventional photonic crystals system, experimenters measure the topological phase transitions through the reflection phase changes [14]. By making analogies, we believe that the topological phase transition of the PTC system can be observed by measuring its reflection phase in the experiment too.

3.3 Band structure of the PTC with a dispersive medium

Finally, we examine the band structure of the PTC system when region B is a dispersive medium. We assume that the dispersive permittivity and permeability are described using ϵ2=ϵ20α/k2 and μ2=μ20β/k2. We subsequently take ϵ20=μ20=1, α=β=100 and keeping the other physical parameters unchanged (ϵ1=4,μ1=1). This yields the band structure for different unit cell sizes with the same modulation frequency, which are shown in Figs. 4(a) and 4(b). We can find a special k-gap whose size and edges are invariant to modulation time scaling for a given modulation frequency. One thing to note is that this special k-gap is the only gap exhibiting this characteristic, which is different from the other k-gaps that behave as Bragg gaps do. In effective medium theory, this special k-gap is a zero-averaged refractive index gap; the wave vectors of the band edges are corresponding to the two solutions of the zero-averaged refractive index, which have no connection with the modulation time scaling. In contrast, other Bragg k-gaps do not satisfy the zero-averaged refractive index condition and are therefore sensitive to the modulation time scaling.

 figure: Fig. 4

Fig. 4 Band structure of the PTC with a dispersive medium for different unit cell sizes and the same modulation frequency: (a) t1=4ns,t2=8ns, (b) t1=3ns,t2=6ns.

Download Full Size | PDF

4. Conclusions

In summary, we calculated the band structure of PTC, the topological invariants of bulk bands and the variation of the k-gap size. We also obtained the correspondence between the topological phase transition and the modulation parameters of the PTC system. Moreover, a temporal zero-averaged refractive index k-gap that is invariant to modulation time scaling for a given modulation frequency is obtained when the PTC system includes a dispersive medium. Our results provide a new path for controlling EM waves and provide topological insight into PTC systems.

Funding

National Key R&D Program of China (2017YFA0303400) and National Natural Science Foundation of China (NSFC) (91630313).

References

1. N. Wang, Z. Q. Zhang, and C. T. Chan, “Photonic Floquet media with a complex time-periodic permittivity,” Phys. Rev. B 98(8), 085142 (2018). [CrossRef]  

2. J. S. Martínez-Romero and P. Halevi, “Parametric resonances in a temporal photonic crystal slab,” Phys. Rev. A (Coll. Park) 98(5), 053852 (2018). [CrossRef]  

3. E. Lustig, Y. Sharabi, and M. Segev, “Topological aspects of photonic time crystals,” Optica 5(11), 1390–1395 (2018). [CrossRef]  

4. W. B. Cairncross, D. N. Gresh, M. Grau, K. C. Cossel, T. S. Roussy, Y. Ni, Y. Zhou, J. Ye, and E. A. Cornell, “Precision measurement of the electron’s electric dipole moment using trapped molecular ions,” Phys. Rev. Lett. 119(15), 153001 (2017). [CrossRef]   [PubMed]  

5. A. Y. Petrov, “Non-reciprocal modulation via acousto-optics,” Nat. Photonics 12(10), 570 (2018). [CrossRef]  

6. L. Zeng, J. Xu, C. Wang, J. Zhang, Y. Zhao, J. Zeng, and R. Song, “Photonic time crystals,” Sci. Rep. 7(1), 17165 (2017). [CrossRef]   [PubMed]  

7. J. R. Zurita-Sánchez, P. Halevi, and J. C. C. González, “Reflection and transmission of a wave incident on a slab with a time-periodic dielectric function ϵ(t),” Phys. Rev. A 79(5), 053821 (2009). [CrossRef]  

8. K. Narayanan, A. W. Elshaari, and S. F. Preble, “Broadband all-optical modulation in hydrogenated-amorphous silicon waveguides,” Opt. Express 18(10), 9809–9814 (2010). [CrossRef]   [PubMed]  

9. Y. Akihama, Y. Kanamori, and K. Hane, “Ultra-small silicon waveguide coupler switch using gap-variable mechanism,” Opt. Express 19(24), 23658–23663 (2011). [CrossRef]   [PubMed]  

10. S. Vezzoli, V. Bruno, C. DeVault, T. Roger, V. M. Shalaev, A. Boltasseva, M. Ferrera, M. Clerici, A. Dubietis, and D. Faccio, “Optical Time Reversal from Time-Dependent Epsilon-Near-Zero Media,” Phys. Rev. Lett. 120(4), 043902 (2018). [CrossRef]   [PubMed]  

11. L. Lu, J. D. Joannopoulos, and M. Soljačić, “Topological photonics,” Nat. Photonics 8(11), 821–829 (2014). [CrossRef]  

12. A. B. Khanikaev and G. Shvets, “Two-dimensional topological photonics,” Nat. Photonics 11(12), 763–773 (2017). [CrossRef]  

13. X. Zhu, H. X. Wang, C. Xu, Y. Lai, J. H. Jiang, and S. John, “Topological transitions in continuously deformed photonic crystals,” Phys. Rev. B 97(8), 085148 (2018). [CrossRef]  

14. Q. Wang, M. Xiao, H. Liu, S. Zhu, and C. T. Chan, “Measurement of the Zak phase of photonic bands through the interface states of a metasurface/photonic crystal,” Phys. Rev. B 93(4), 041415 (2016). [CrossRef]  

15. W. Zhu, X. Fang, D. Li, Y. Sun, Y. Li, Y. Jing, and H. Chen, “Simultaneous observation of a topological edge state and exceptional point in an open and non-Hermitian acoustic system,” Phys. Rev. Lett. 121(12), 124501 (2018). [CrossRef]   [PubMed]  

16. A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, “The electronic properties of graphene,” Rev. Mod. Phys. 81(1), 109–162 (2009). [CrossRef]  

17. R. M. Lutchyn, J. D. Sau, and S. Das Sarma, “Majorana fermions and a topological phase transition in semiconductor-superconductor heterostructures,” Phys. Rev. Lett. 105(7), 077001 (2010). [CrossRef]   [PubMed]  

18. B. Bahari, A. Ndao, F. Vallini, A. El Amili, Y. Fainman, and B. Kanté, “Nonreciprocal lasing in topological cavities of arbitrary geometries,” Science 358(6363), 636–640 (2017). [CrossRef]   [PubMed]  

19. C. L. Kane and E. J. Mele, “Z2 topological order and the quantum spin Hall effect,” Phys. Rev. Lett. 95(14), 146802 (2005). [CrossRef]   [PubMed]  

20. Q. Wang, M. Xiao, H. Liu, S. Zhu, and C. T. Chan, “Optical interface states protected by synthetic Weyl points,” Phys. Rev. X 7(3), 031032 (2017). [CrossRef]  

21. D. T. Tran, A. Dauphin, N. Goldman, and P. Gaspard, “Topological Hofstadter insulators in a two-dimensional quasicrystal,” Phys. Rev. B Condens. Matter Mater. Phys. 91(8), 085125 (2015). [CrossRef]  

22. B. Y. Xie, H. F. Wang, X. Y. Zhu, M. H. Lu, Z. D. Wang, and Y. F. Chen, “Photonics meets topology,” Opt. Express 26(19), 24531–24550 (2018). [CrossRef]   [PubMed]  

23. X. Y. Guo, Z. Wang, Z. Chi, Q. Zheng, and X. Ren, “Topological polaritons with squeezed photons,” J. Opt. Soc. Am. B 33(6), 1128–1133 (2016). [CrossRef]  

24. M. Xiao, G. Ma, Z. Yang, P. Sheng, Z. Q. Zhang, and C. T. Chan, “Geometric phase and band inversion in periodic acoustic systems,” Nat. Phys. 11(3), 240–244 (2015). [CrossRef]  

25. P. Delplace, D. Ullmo, and G. Montambaux, “Zak phase and the existence of edge states in graphene,” Phys. Rev. B Condens. Matter Mater. Phys. 84(19), 195452 (2011). [CrossRef]  

26. W. Gao, M. Xiao, B. Chen, E. Y. B. Pun, C. T. Chan, and W. Y. Tam, “Controlling interface states in 1D photonic crystals by tuning bulk geometric phases,” Opt. Lett. 42(8), 1500–1503 (2017). [CrossRef]   [PubMed]  

27. M. Xiao, Z. Q. Zhang, and C. T. Chan, “Surface impedance and bulk band geometric phases in one-dimensional systems,” Phys. Rev. X 4(2), 021017 (2014). [CrossRef]  

28. K. H. Choi, C. W. Ling, K. F. Lee, Y. H. Tsang, and K. H. Fung, “Simultaneous multi-frequency topological edge modes between one-dimensional photonic crystals,” Opt. Lett. 41(7), 1644–1647 (2016). [CrossRef]   [PubMed]  

29. W. P. Su, J. R. Schrieffer, and A. J. Heeger, “Solitons in polyacetylene,” Phys. Rev. Lett. 42(25), 1698–1701 (1979). [CrossRef]  

30. M. Deopura, C. K. Ullal, B. Temelkuran, and Y. Fink, “Dielectric omnidirectional visible reflector,” Opt. Lett. 26(15), 1197–1199 (2001). [CrossRef]   [PubMed]  

31. S. Ji, Y. Hong, P. S. Spencer, J. Benedikt, and I. Davies, “Broad tunable photonic microwave generation based on period-one dynamics of optical injection vertical-cavity surface-emitting lasers,” Opt. Express 25(17), 19863–19871 (2017). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1
Fig. 1 Schematic diagram of the PTC system.
Fig. 2
Fig. 2 Band structure of the PTC with different modulation frequencies: (a) f = 0.3 , (b) f = 0.4 , and (c) f = 0.5. (d) The displacement field distribution of the compound PTC system constructed by ( P T C 1 ) 10 and ( P T C 2 ) 10 . The subsystems PTC1 and PTC2 correspond to the systems (a) and (c), respectively.
Fig. 3
Fig. 3 Variation of the second k -gap size Δ k with modulation frequency f and modulation strength g : (a) ϵ 2 = 4 and (b) ϵ 2 = 9.
Fig. 4
Fig. 4 Band structure of the PTC with a dispersive medium for different unit cell sizes and the same modulation frequency: (a) t 1 = 4 n s , t 2 = 8 n s , (b) t 1 = 3 n s , t 2 = 6 n s .

Equations (8)

Equations on this page are rendered with MathJax. Learn more.

ϵ ( t ) = { ϵ 1 ( n T t t 1 + n T ) ϵ 2 ( t 1 + n T t ( n + 1 ) T ) ,
E i ( t ) = A 0 e i ( ω i t + k r ) + B 0 e i ( ω i t k r ) ,
t 0 t 0 + ( × E ) d t = t 0 t 0 + B t d t = [ B ( t 0 + ) B ( t 0 ) ]
t 0 t 0 + ( × H ) d t = t 0 t 0 + D t d t = D ( t 0 + ) D ( t 0 )
B | t 0 = B | t 0 +
D | t 0 = D | t 0 + .
cos ( Ω T ) = cos ω 1 t 1 cos ω 2 t 2 1 2 ( z 1 z 2 + z 2 z 1 ) sin ω 1 t 1 sin ω 2 t 2 ,
θ m Z a k = π / T π / T [ i u n i t c e l l d t ε ( t ) u m , Ω * ( t ) Ω u m , Ω ( t ) ] d Ω ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.