Arxiv:2002.05124V5 [Cond-Mat.Mes-Hall] 4 Jun 2020 [45]
Total Page:16
File Type:pdf, Size:1020Kb
Effective Floquet Hamiltonians for periodically-driven twisted bilayer graphene Michael Vogl∗ Department of Physics, The University of Texas at Austin, Austin, TX 78712, USA Martin Rodriguez-Vega∗ Department of Physics, The University of Texas at Austin, Austin, TX 78712, USA and Department of Physics, Northeastern University, Boston, MA 02115, USA Gregory A. Fiete Department of Physics, Northeastern University, Boston, MA 02115, USA and Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA (Dated: June 5, 2020) We derive effective Floquet Hamiltonians for twisted bilayer graphene driven by circularly polar- ized light in two different regimes beyond the weak-drive, high frequency regime. First, we consider a driving protocol relevant for experiments with frequencies smaller than the bandwidth and weak amplitudes and derive an effective Hamiltonian, which through a symmetry analysis, provides an- alytical insight into the rich effects of the drive. We find that circularly polarized light at low frequencies can selectively decrease the strength of AA-type interlayer hopping while leaving the AB-type unaffected. Then, we consider the intermediate frequency, and intermediate-strength drive regime. We provide a compact and accurate effective Hamiltonian which we compare with the Van Vleck expansion and demonstrate that it provides a significantly improved representation of the exact quasienergies. Finally, we discuss the effect of the drive on the symmetries, Fermi velocity and the gap of the Floquet flat bands. I. INTRODUCTION an exponentially-long pre-thermal time regime [59{64] in driven interacting quantum systems allows one to intro- The recent discovery of strong-correlation effects in duce the notion of effective time-independent theories. twisted bilayer graphene (TBG) generated great inter- The development of several techniques to derive effec- est in moir´eheterostructures [1{31] and ways to sim- tive Hamiltonians in different drive regimes led to rapid ulate them [23]. Similar to the behavior in cuprates evolution of the Floquet engineering field [61, 65{86]. [32, 33] at different filling factors superconductivity, For instance, the prediction of an anomalous Hall effect Mott-insulating [4, 16, 34{36] and ferromagnetic be- in single-layer graphene driven by circularly polarized haviour [37, 38] has been observed in TBG. The ex- light [87] has been recently confirmed in experiments [88]. perimental observations were followed by several theo- More generally, there has been an increased interest in retical proposals to explain the observations based on the study of topological transitions induced by periodic the existence of flat bands which appear at special twist drives [87, 89{107]. angles [1, 9, 39]. These flat bands play an essential More recently, the fields of twistronics and Floquet en- role for the emergence of strong correlations because gineering crossed paths in twisted bilayer graphene driven the interaction terms become relatively dominant [40] by circularly polarized light in free space [108{110]. In- over the kinetic energy contributions of the dispersive teresting effects like topological transitions at large twist bands [12, 24, 36, 39{41]. angles using high-frequency drives [108] and the induc- In TBG, the flat bands depend strongly on the twist tion of flat bands using near-infrared light in a wide range angle between the graphene layers, which is experimen- of twist angles [110] were found. These studies are mainly tally difficult to set to a precise values. This challenge has numerical, and only provide analytical descriptions in the lead to several studies proposing different mechanisms to high drive frequency regime, which we will define rigor- correct for deviations from the magic angle. For example, ously in the next section. via pressure [5, 42{44] or light confined in a waveguide The aim of this work is to derive analytical effective arXiv:2002.05124v5 [cond-mat.mes-hall] 4 Jun 2020 [45]. Floquet Hamiltonians that allow us to gain insight into In parallel to the developments on moir´elattices, there twisted bilayer graphene subjected to circularly polarized has been a rapid progress in our understanding of non- light away from the conventional weak drive, high fre- equilibrium systems, both experimentally and theoreti- quency regime of Van Vleck [111, 112], Floquet-Magnus cally, particularly for the case of periodic drives, which or Brillouin-Wigner approximations [74, 113]. Our effec- may be induced by a laser. [46{58]. The existence of tive Floquet Hamiltonians allow us to elucidate the ef- fects of the interplay of moir´elattices and Floquet drives. Particularly, we consider two complementary regimes: i) a regime characterized by weak drive and low frequencies; ∗ These two authors contributed equally. and ii) a regime characterized by intermediate frequen- 2 cies and strong drives. The remainder of the manuscript (a) (b) is organized as follows: in Sec. II we describe the sys- tem we consider; in Sec.III we examine the low-frequency, weak drive limit; and in Sec.IV we address the interme- diate frequency and intermediate strength drive regime. Finally, in Sec.V we present our conclusions and outlook. II. SYSTEM DESCRIPTION A. Static Hamiltonian (c) The starting point of our discussion is the effec- tive Hamiltonian that describes twisted bilayer graphene [1, 25, 39, 114{116] h(−θ=2; k − κ−) T (x) Hk(x) = y ; (1) T (x) h(θ=2; k − κ+) which describes two stacked graphene layers that are ro- tated with respect to each other by an angle θ, as shown in the sketch of figure 1(a). Here, 0 f(R(θ)k) h(θ; k) = γ ; (2) f ∗(R(θ)k) 0 FIG. 1. (Color online) a) Sketch of twisted bilayer graphene is the single-layer graphene Hamiltonian, f(k) = irradiated by circularly polarized light. (b) moir´eBrilloiun ia k − 2 ia k 0 y a k π e 3 0 y +2e 3 sin p0 x − describes the intralayer zone. (c) Band structure for twisted bilayer graphene for w0 = 3 6 ◦ w1 = 110 meV, and θ = 1:05 . The low-energy flat bands are hopping amplitude between nearest-neighbor sites, and highlighted in red. γ = vF =a0, where we use natural units ~ = c = e = 1. The inclusion of the full structure of f(k) means that this Hamiltonian is valid in the full Brillouin zone and not just near a K point. The interlayer hopping matrix would incorporate the two graphene valleys. However, 1 we only consider perturbations induced by light, which X −ibix T (x) = e Ti; (3) cannot induce processes that mix the two valleys. The i=−1 Hamiltonian in the other valley is connected by a C2 ro- 2πn 2πn tation [119]. The symmetries of the continuum model Eq. Ti = w012 + w1 cos σ1 + sin σ2 ; (4) (1) include C3 rotational symmetry about the center of a 3 3 AA region, C2T symmetry (taking into account both val- describes tunneling between the two graphene layers and leys, the TBG presents time-reversal symmetry T ), and encodes a hexagonal pattern that has its origin in that My : y; ky ! −y; −ky mirror symmetry [16, 119, 120]. the two superimposed graphene lattices which develop a In the small-rotation limit, the angle dependence of the moir´epattern (see Fig. 1(b)), where b0 = (0; 0), and graphene sectors can be neglected, leading to an approx- p imate particle-hole symmetry C [120]. b±1 = kθ ± 3=2; 3=2 are the reciprocal lattice vectors. Following Refs. 45, 109, and 110 we introduced an addi- tional parameter w1 into the tunneling term to model re- laxation effects, since AB/BA stacking configurations are B. Driven twisted bilayer graphene energetically favoured over AA configurations [115, 117]. Furthermore, there are indications that AA and AB For the driven system, we assume that circularly po- regions have different interlayer-lattice constants [118]. larized light is applied in a direction normal to the TBG Throughout this work, we fixed γ = vF =a0 = 2:36 eV, plane as sketched in Fig. 1(a). Then, the light enters and a = 2:46 A.˚ For a detailed description of the band 0 via minimal substitution as k ! k~ (t) = k −A cos(Ωt), structure numerical implementation, see the appendix x x x k ! k~ k − A t of [45]. In figure 1(c) we show the band structure for and y y = y sin(Ω ) leaving the tunneling sec- w = w = 110 meV, and θ = 1:05◦, value near the tor almost unaltered. The reason for this is simple. The 0 1 inclusion of light in a tight binding model can be done magic angle. R Rj i R Adr The Hamiltonian in Eq. (1) describes only one val- via a Peirls substution for hoppings tij ! e i tij. ley degree of freedom. A full description of the system The interlayer hopping is dominated by hopping between 3 atoms that are almost exactly on top of each other - as j (t)i = eitjφ(t)i, where jφ(t + 2π=Ω)i = jφ(t)i afterall other atoms are further away and the overlap and is the quasienergy. Replacing j (t)i into the between orbitals is smaller. Therefore for interlayer cou- Schr¨odinger equation leads to [H(x; t) − i@t]jφ(t)i = plings mostly longitudinal components of A contribute in jφ(t)i, which governs the dynamics of the periodic sys- R the line integral R j Adr. Circularly polarized light only tem. The exact solution can be generically obtained Ri has transverse components and therefore has little effect either by constructing the Floquet evolution operator 2π=Ω R −iHF T on interlayer couplings. The time-dependent Hamilto- UF = T exp{−i 0 H(s)dsg = e or by employ- nian is ing the extended-state picture. In the extended-state pic- P inΩt ture, we use the Fourier series jφ(t)i = n e jφni, which leads to P H(n−m) + δ Ωm jφ i = jφ i, m n;m m n h(−θ=2; k~(t) − κ−) T (x) defined in the infinite-dimensional Floquet-Hilbert space H(x; t) = y ; T (x) h(θ=2; k~(t) − κ+) spanned by the direct product of the Hilbert space of the (5) static system and the space spanned by a complete set with H(x; t + 2π=Ω) = H(x; t).