Effective Hamiltonians, prethermalization and slow energy absorption in periodically driven many-body systems

Effective Hamiltonians, prethermalization and slow energy absorption in periodically driven many-body systems

Dmitry A. Abanin Department of Theoretical Physics, University of Geneva    Wojciech De Roeck Instituut voor Theoretische Fysica, KU Leuven, Belgium    Wen Wei Ho Department of Theoretical Physics, University of Geneva    François Huveneers CEREMADE, Universite Paris-Dauphine, France
July 12, 2019
Abstract

We establish some general dynamical properties of quantum many-body systems that are subject to a high-frequency periodic driving. We prove that such systems have a quasi-conserved extensive quantity , which plays the role of an effective static Hamiltonian. The dynamics of the system (e.g., evolution of any local observable) is well-approximated by the evolution with the Hamiltonian up to time , which is exponentially large in the driving frequency. We further show that the energy absorption rate is exponentially small in the driving frequency. In cases where is ergodic, the driven system prethermalizes to a thermal state described by at intermediate times , eventually heating up to an infinite-temperature state after times . Our results indicate that rapidly driven many-body systems generically exhibit prethermalization and very slow heating. We briefly discuss implications for experiments which realize topological states by periodic driving.

pacs:
73.43.Cd, 05.30.Jp, 37.10.Jk, 71.10.Fd

I Introduction

Recent advances in laser cooling techniques have resulted in experimental realizations of well-isolated, highly tunable quantum many-body systems of cold atoms Bloch et al. (2008). A rich experimental toolbox of available quantum optics, combined with the systems’ slow intrinsic time scales, allow for a preparation of non-equilibrium many-body states and also a precise characterization of their quantum evolution. This has made the study of different dynamical regimes in many-body systems one of the forefront directions in modern condensed matter physics (for a review, see Ref. Polkovnikov et al., 2011).

Conventional wisdom suggests that in a majority of many-body systems, the Hamiltonian time evolution starting from a non-equilibrium state should lead to thermalization at sufficiently long times: that is, physical observables reach thermal values, given by the microcanonical ensemble. Thermalization in such ergodic systems is understood in terms of the properties of individual eigenstates themselves – observables measured in these eigenstates are already thermal, as encapsulated by the eigenstate thermalization hypothesis (ETH) Deutsch (1991); Srednicki (1994); Rigol et al. (2008). However, while ETH implies eventual thermalization, it does not make predictions regarding the intermediate-time dynamics of the system. Therefore, much work has been dedicated to studying how thermal equilibrium emerges in different many-body systems.

In particular, there is a class of systems, which exhibit the phenomenon of prethermalization Berges et al. (2004); Eckstein et al. (2009); Moeckel and Kehrein (2010); Mathey and Polkovnikov (2010). Such systems have a set of approximate conservation laws, in addition to energy; therefore, at intermediate time scales they equilibrate to a state given by the generalized Gibbs ensemble, which is restricted by those conservation laws. Full thermal equilibrium is reached at much longer time scales, set by the relaxation times of the approximate integrals of motion. Prethermalization has been experimentally observed in a nearly integrable one-dimensional Bose gas Gring et al. (2012).

In this paper, we establish some general properties of dynamics of periodically driven many-body systems (Floquet systems). Periodic driving in quantum systems has recently attracted much theoretical and experimental attention, because, amongst many applications, it provides a tool for inducing effective magnetic fields, and for modifying topological properties of Bloch bands Oka and Aoki (2009); Kitagawa et al. (2010); Lindner et al. (2011); Rudner et al. (2013). Indeed, since periodic driving is naturally realized in cold atomic systems by applying electromagnetic fields, topologically non-trivial Bloch bands (Floquet topological insulators) in non-interacting systems have been observed experimentallyAidelsburger et al. (2013); Jotzu et al. (2014); Aidelsburger et al. (2015). However, since periodic driving breaks energy conservation, driven ergodic (many-body) systems are expected to heat up, eventually evolving into a featureless, infinite-temperature state Ponte et al. (2015a); D’Alessio and Rigol (2014); Lazarides et al. (2014); foo (). Thus, many-body effects are expected to generally make such Floquet systems unstable. Below, we derive general bounds for energy absorption rates in periodically driven many-body systems, which can be applied for instance to understand the lifetimes of Floquet topological insulators.

As the main result of the paper, we show that rapidly driven many-body systems with local interactions generally have a local, quasi-conserved extensive quantity, , which plays the role of an effective Hamiltonian. At times , the time evolution of any local observable is well-approximated by the Hamiltonian evolution with the time independent Hamiltonian . Thus, assuming that the Hamiltonian is ergodic, the system exhibits prethermalization to a thermal state described by the Hamiltonian , with an effective temperature set by the initial “energy” . The quasi-conservation of is destroyed at timescale , when energy absorption occurs and an infinite-temperature state is formed. We show that the heating timescale is exponentially large in the driving frequency :

(1)

where is a numerical constant of order 1, and has the meaning of a maximum energy per particle or spin, precisely defined below. Thus, rapidly driven many-body systems generically have a very long prethermalization regime, and absorb energy exponentially slowly in the driving frequency. We emphasize that these results are non-perturbative; they generalize and complement our previous work, Ref. Abanin et al., 2015, where bounds on linear-response heating rates were proven. As an implication of our result, we show that the measurement of a local operator time evolved with the effective Hamiltonian is close to the measurement of the same operator but exactly time evolved, up to exponentially long times.

The structure of the rest of the paper is as follows. In Sec. II, we define the set-up and present the central idea of the transformation used to obtain our results. Then, in Sec. III, we work out in detail, using the method presented, the optimal order of the transformation to obtain the effective Hamiltonian and also the heating time scale for which this effective Hamiltonian is valid. Next, in Sec. IV, we present the implications of our result for the observation of a local operator. Lastly, we end with a discussion in Sec. V.

Ii Set-up and outline of method

We consider a quantum many-body system subject to a drive with a period , described by a time dependent Hamiltonian:

(2)

where is time independent, and, without loss of generality, the time average of the driving term is chosen to be zero, . We focus on the case of a lattice system with locally bounded Hilbert space. In other words, the Hilbert space of site is finite-dimensional, as is the case for fermions, spins, as well as hard-core bosons. We also restrict to one-dimensional systems, but this is not crucial to the method, see also Ref. Abanin et al., 2015. Both and are assumed to be local many-body operators, that is, they can be written as a sum of local terms:

(3)

where runs over all lattice sites, . The locality of the interactions means that each term acts non-trivially on at most adjacent sites . (e.g., for the nearest-neighbor Heisenberg model, ); we refer to as the range of the operator. Each term is bounded by a constant interaction strength :

(4)

We will focus on the case when the driving frequency is large (or equivalently, the driving period is small) compared to these local energy scales, that is, .

Now, the unitary dynamics of the system is described by the time evolution operator , which obeys the equation:

(5)

where is the identity operator.

Floquet theory (for a review, see Ref. Bukov et al., 2015) predicts that the solution of Eq.(5) can be written in the following form:

(6)

where is a time periodic unitary such that , and is a time independent Floquet Hamiltonian. In particular, the evolution operator over one period is given by:

(7)

Thus, the evolution of the system at stroboscopic times is governed by the time independent Hamiltonian . Note that the choice of is not unique: given a particular and projectors onto its eigenstates with eigenvalues , the Hamiltonian is also a valid Floquet Hamiltonian for any .

Typically, there is no closed-form solution of Eq.(6), and one relies on iterative schemes such as the Magnus expansion to obtain for high-frequency drives (for a recent review, see Refs. Bukov et al., 2015; Blanes et al., 2009; Eckardt and Anisimovas, 2015). In this approach, is expanded in terms of powers of (equivalently, of inverse frequency ), , where . The formal solution of Eqs.(5, 6) then gives expressed in terms of nested commutators of at different times. However, the Magnus expansion is only known to converge for bounded Hamiltonians, such that , with ,  Blanes et al. (2009). Since many-body systems have extensive energies and do not satisfy this condition, the Magnus expansion is expected not to converge in this case. Indeed, the existence of a quasi-local Floquet Hamiltonian would imply that the system does not heat up to an infinite-temperature state at long times, contrary to the general arguments based on the ETH Ponte et al. (2015a).

Therefore, we propose an alternative approach. The central idea is as follows: we unitarily transform the Hamiltonian, systematically removing time dependent terms at increasing order in . Truncating the procedure at some optimal order (defined below), we obtain a quasi-conserved time independent Hamiltonian operator .

More concretely, we transform the system’s wavefunction by a time periodic unitary , such that :

(8)

Importantly, the wavefunction coincides with the original wavefunction at stroboscopic times . Its evolution is described by the Schroedinger equation

(9)

with a modified Hamiltonian:

(10)

Thus, the transformation defines a new periodic Hamiltonian , which gives the same stroboscopic evolution as the original Hamiltonian .

For our purposes, it is convenient to write the operator as an exponential of a periodic operator , which is anti-Hermitian, , and to represent as an -degree polynomial in the driving period :

(11)

Here, the order of the polynomial should be treated as a parameter to be optimized in a manner described below. Using Duhamel’s formula, for , Eq. (10) can be rewritten as follows:

(12)

where , which gives an expansion of naturally in powers of .

We will show below that the operators can be chosen to get rid of the time dependence of of order for , leaving behind a time-dependent piece of order . Furthermore, we will show that for a given and , there exists an optimal , for which this driving term’s norm (suitably defined) becomes minimal. For a many-body system with local interactions, we find that the optimal , and for this , the driving term’s norm is exponentially reduced by a factor of . The time-independent part of the corresponding Hamiltonian then represents a quasi-conserved energy, valid for an exponentially long time .

Iii Method, optimal order, and heating time scale

We now utilize the transformation outlined in the previous section to transform the original Hamiltonian. We derive the optimal order at which the remaining driving term becomes minimal, which gives us both the effective Hamiltonian and the heating time scale .

iii.1 Simple example: single rotating frame transformation

To get some familiarity regarding the use of our approach before going into full generality, it is instructive to first consider the simple example of a transformation for , i.e., a single rotating frame transformation, so that , and is chosen such that the driving term of order is eliminated in Eq. (12).

Since the zeroth-order contribution in (12) is given by , we define by:

(13)

With this choice of , of Eq. (12) can be expanded in “powers of ”,

(14)

where is the term of order :

(15)

To the first order in , the rotated Hamiltonian is given by:

where is the time-independent part of , and is the new driving term (with zero time-average) at this order. A straightforward calculation shows that coincides with the second order of the Magnus expansion.

We see that is the order piece of and is time independent, while the remaining piece that appears at orders and higher is still time dependent, and represents the new driving term. Thus the rotated Hamiltonian can be written as

(16)

Contrasted to the original Hamiltonian , it appears that the new driving term’s norm has been reduced by a factor of .

However, there is an important distinction to be made between and the original Hamiltonian . In a many-body system, the rotated Hamiltonian in Eq.(14) is now quasi-local instead of being strictly local. This is because involves nested commutators of and , and the norm of each term decreases exponentially with for sufficiently rapid driving. To establish this, we note that each term is extensive and can be written as . We denote the maximum local(l) norm of as , and use the following fact: for any two extensive operators , of range , respectively, such that , , has a range of at most , and , with norm

(17)

This is because each operator can commute non-trivially with at most operators . Repeatedly applying this estimate to the operators , that enter Eq.(15), and using the fact that has range , and (which follows from Eq.(13)), we obtain:

(18)

and the range of equals . Thus, this establishes the quasi-locality of the Hamiltonian .

Further, by using an appropriately weighted local norm (see appendix (A.3) and also Ref. Abanin et al., 2015), the size of is . Therefore, we see that the transformation reduces the amplitude of the time dependent term by a factor of order , while at the same time making the Hamiltonian quasi-local, and renormalizing its time independent part.

iii.2 General case

Next, we proceed to the general case of . Then, , where as mentioned, is chosen such that the only time dependent terms in the Hamiltonian are of order . This condition gives us a set of recursive relations for : we use them to ‘absorb’ the time dependent pieces of order in for . To derive these relations, we first note that the term of the order in has the following form:

(19)

where is expressed in terms of :

(20)

(For , we set .) We can separate into a time independent part, , and a time dependent part with zero average over one period:

(21)

We eliminate the time dependent term of the order in (see Eq.(19)) by choosing as follows:

(22)

for . In particular, for , is given by Eq.(13).

Relations (20,21,22) define the transformation which makes the time dependent terms in the Hamiltonian of the order :

(23)

In a manner similar to the simple example of considered before, the full Hamiltonian and time dependent term can be shown to be quasi-local (see appendix (A.3)).

Now, let us now estimate the norm of . We argue that there is an optimal which we call , for which the procedure we have outlined before approximatively minimizes the local norm of . This has physical consequences for both the heating time scale and the observation of a local observable, for example. Thus, should be chosen as .

To this end, we prove a number of inequalities for the norms of various operators, and (refer to the appendix and to Ref. Abanin et al., 2015 for generalizations). In the following, there will appear constants , etc., which depend on the microscopic details of the system such as and , but importantly not on the driving period . It is to be understood that these constants can be different for different objects in question that are being bounded. Now, for , we have

(24)

with a combinatorial constant of order . The other operators have then derived bounds since and . For , we have

(25)

The factor in the above bounds arises because of the many-body nature of the system: involves nested commutators of , . Eq. (24) shows that there are two competing effects which control the behavior of : suppression of by a factor of , and its growth due to . Eventually, the factorial dominates and therefore for , the local norm of stops decreasing with .

The optimal that we have to choose for is roughly the same as the one to choose to minimize the norm of or (see appendix). From the right hand side of Eq. (24) or Eq. (25), we obtain

(26)

with (independent of ) defined in appendix (A.3). This gives us the following bound on ,

(27)

for , which in turn gives us an estimate on the remainder:

(28)

As already indicated, the truly useful version of this bound also expresses that local terms in with large range are additionally damped, and this is indeed captured by the use of a stronger norm in the appendix.

Furthermore, at this optimal order, the time independent part of the transformed Hamiltonian is a physical, local many-body Hamiltonian

(29)

and differs from the original Hamiltonian by a sum of small local terms, more precisely

(30)

.

Eq. (28) together with (30) imply that the energy absorption rate (per volume) is exponentially small, giving us a characteristic heating time scale that scales like

(31)

The operator is therefore a quasi-conserved extensive quantity, playing the role of an effective static Hamiltonian, and it can be used to accurately describe stroboscopic dynamics up to times .

Iv Implications: evolution of a local observable

Next, we spell out the consequences of the existence of this effective Hamiltonian for the time evolution of a local observable , with . Let us consider the difference between evolved in time using the exact time evolution operator and the time evolution generated by the effective Hamiltonian . The difference

(32)

can be recast, using the frame transformation and Duhamel’s formula, as

(33)

where is the evolution from time to generated by the time dependent Hamiltonian . The norm of the difference can be bounded using the unitarity of as

(34)

which can be controlled by the Lieb-Robinson bound, see Refs. Lieb and Robinson, 1972; Nachtergaele and Sims, 2010. Indeed, let us first pretend that the range of local terms in is maximally , then the Lieb-Robinson bound yields

(35)

where is the Lieb-Robinson velocity of , which can be chosen to be . Here is a numerical constant of order .

The bound (35) expresses that only those terms in that have support within distance of the support of , contribute to the commutator, see Ref. Abanin et al., 2015 for a more detailed derivation of such bounds. Since in our case the support of local terms in can grow arbitrarily large (because it is quasi-local), we however need to use the exponential decay in support of the norm of each local term in to derive Eq. (35), in which case depends on the decay constant. We omit this straightforward calculation and refer to Ref. Abanin et al., 2015.

Using , we conclude that the difference Eq. (32) grows as with and hence it remains small up to an exponentially long time . Thus, a measurement of that is time evolved by the effective Hamiltonian will be close to the measurement of that is time evolved by the exact Hamiltonian , for an exponentially long time.

V Discussion and conclusion

In this paper, we considered many-body systems subject to a high-frequency periodic driving. We have shown that there is a broad time window, , in which stroboscopic dynamics of such systems is controlled by an effective time independent Hamiltonian . We have used a series of “gauge”, time periodic unitary transformations to effectively reduce the strength of the driving term and to establish the existence of . The advantage of our approach compared to the standard Magnus expansion Bukov et al. (2015); Blanes et al. (2009); Eckardt and Anisimovas (2015) is that it allows us to control the magnitude of the driving terms after the transformations.

We note that recently Canovi et al. Canovi et al. (2016) and Bukov et al. Bukov et al. (2015) discussed prethermalization in weakly interacting driven systems. Our results complement these works: we have shown that (rapidly) driven interacting systems generically exhibit a broad prethermalization regime, which can be observed in a quench experiment as follows. Let us initially prepare the system in some non-equilibrium state , and subject it to a rapid periodic drive. At times the system will reach a steady state, in which physical observables have thermal values, , where the density matrix , with being the effective temperature set by the energy density of the initial state. Thus, at times the system appears as if it is not heating up. The system will absorb energy and relax to a featureless, infinite-temperature state beyond times . We expect this phenomenon to be observable in driven system of cold atoms and spins (assuming relaxation of spins due to phonons is slow).

Finally, we briefly discuss the implications of our results for the current efforts to realize topologically non-trivial strongly correlated states (e.g., fractional Chern insulators) in periodically driven systems. Experimentally, one tries to design a drive for which the ground state of an effective time independent Hamiltonian (usually calculated within low-order Magnus expansion) is topologically non-trivial. A central challenge is to prepare the system in a ground state of the effective Hamiltonian. Since we have shown that the dynamics of the system is controlled by up to exponentially long times, one can envision that the “Floquet fractional Chern insulators” can be prepared as follows. Let us assume that the system can be initially prepared in a (topologically trivial) ground state of the Hamiltonian . Then, the driving is switched on adiabatically to the value which corresponds to the desired effective Hamiltonian . However, the switching should also be done quickly compared to to avoid energy absorption. Since and describe different phases, the system will necessarily go through a quantum critical point (QCP), and excitations will be created via a Kibble-Zurek mechanism. The number of excitations can be minimized by designing a non-linear passage through the QCP Barankov and Polkovnikov (2008). We leave a detailed exploration of these ideas for future workHo and Abanin (2016).

Note added. Recently, a related result, Ref. Mori et al., 2016 appeared, building on Ref. Kuwahara et al., 2016 (local driving). Ref. Mori et al., 2016 proves a similar bound for the absorption rate in driven systems using a different approach (namely, studying evolution over one driving period).

Vi Acknowledgements

D.A. acknowledges support by Alfred Sloan Foundation. W.D.R. also thanks the DFG (German Research Fund), the Belgian Interuniversity Attraction Pole (P07/18 Dygest), as well as ANR grant JCJC for financial support.

References

Appendix A Technical estimates and proofs

Here, we provide a proof of the bounds on the terms of the renormalized Hamiltonian and the bound on the remainder .

a.1 Setup

Let us first recall the norm that we are using. We write operators , periodic in time, as a sum of local terms

where runs over the sites of the (finite but large) volume and is an operator that acts nontrivially on the sites where , with independent of and called the “range of ”. The local norm is then defined as

In what follows, we mostly drop the dependence on from the notation. Let us now list the important bounds on local norms that we claim: The operators have range and their local norms are bounded as

(S1)
(S2)
(S3)

We write , so that (S1) is consistent with the fact that the range of the original operator is and its local norm is . The bounds in (S2) follow immediately because time averaging of local term does not increase its norm (here we use that the norm was defined as the supremum over time). The bound (S3) follows from (S1) for a given because the integral over one period yields an additional factor , i.e. using . Hence we have now in particular established the above bounds for and we have shown that the bound (S1) implies the others, for a given . Therefore, to complete an inductive proof, it suffices to prove (S1) for while assuming the other bounds for all . To achieve this, we use (20):

(S4)

The right hand side of (S4) is a sum of local operators, all of which have support not greater than adjacent sites. We estimate the norm of each of these local operators by using repeatedly (if have overlapping support) and the above bounds for . Then we sum the bounds on all terms that have site as the leftmost site of their support, to get a bound for . The result is, separately for the first (S5) and second (S6) term of (S4):

(S5)
(S6)

where if statement holds true, and otherwise, and we abbreviated

The sum is over all sequences of discrete intervals (sets of adjacent sites) such that we have the following conditions.

  1. All interval lengths are multiples of : .

  2. For , is nonempty.

  3. .

Intersection condition stems from the structure of nested commutators. Condition 4 says that we consider terms the support of which starts at site .

To conclude the proof of the bounds (S1), we have to show that, for some (-independent) choice of the constant , the sum of (S5) and (S6) is bounded by . It is sufficient to prove a bound on (S6), as (S5) reduces to that case upon increasing , with such that . For the same reason, we can replace by . Hence, we show Lemma 1.

Lemma 1.

For some independent of ,

(S7)

a.2 Proof of Lemma 1

For simplicity, we set . The case follows analogously. We set and we let be such that . Note that because there are at least two overlapping intervals and the sum of their lengths is . We write for the sequence of lengths. First, we dominate the sum on the left-hand side of (S7) as

(S8)

where the factor accounts for the choice of position of the intervals in the stretch . There are intervals but at least one of them has to be placed such that its leftmost point is at , therefore we have instead of .

We use the upper bound in Stirling’s formula

to get

(S9)

Here and below we use for numerical constants that do not depend on , their value can change from line to line. To deal with the product over such factors we define with

where is the number of ’large’ naturals in the sequence :

Plugging (S9) into (S8) and using the above notation, we get

(S8) (S10)
(S11)

where we have also used the lower bound in Stirling’s formula. For , we find the bound

(S12)

In words, short intervals yield small factors , but constrains .

Proof of (S12). Note that implies . There are ways to choose large naturals from . This yields the bound

The sums over the dummy variable are estimated as

where the constant in the first inequality of course depends on and diverges when . ∎

Plugging (S12) into (S11) and using , we get

(S13)

The first sum over is trivially bounded by . For the second sum over , we use to get the upper bound

provided and . It follows that (S13) can be made smaller than by choosing large enough.

a.3 Bound on remainder

We start immediately from the expression (10) that we repeat here

(S14)

We will plug in with and estimate the local terms. Since there will be terms of any range, it is beneficial to introduce here a weighted norm. If with supported on a finite interval , then we put, for some ,

where the sup ranges over sites