Analogue of Hamilton-Jacobi theory for the time-evolution operator Michael Vogl,1, 2 Pontus Laurell,1, 3 Aaron D. Barr,1 and Gregory A. Fiete1, 2, 4 1Department of Physics, The University of Texas at Austin, Austin, TX 78712, USA 2Department of Physics, Northeastern University, Boston, MA 02115, USA 3Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 4Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA (Dated: July 16, 2019) In this paper we develop an analogue of Hamilton-Jacobi theory for the time-evolution operator of a quantum many-particle system. The theory offers a useful approach to develop approximations to the time-evolution operator, and also provides a unified framework and starting point for many well-known approximations to the time-evolution operator. In the important special case of peri- odically driven systems at stroboscopic times, we find relatively simple equations for the coupling constants of the Floquet Hamiltonian, where a straightforward truncation of the couplings leads to a powerful class of approximations. Using our theory, we construct a flow chart that illustrates the connection between various common approximations, which also highlights some missing con- nections and associated approximation schemes. These missing connections turn out to imply an analytically accessible approximation that is the \inverse" of a rotating frame approximation and thus has a range of validity complementary to it. We numerically test the various methods on the one-dimensional Ising model to confirm the ranges of validity that one would expect from the approximations used. The theory provides a map of the relations between the growing number of approximations for the time-evolution operator. We describe these relations in a table showing the limitations and advantages of many common approximations, as well as the new approximations introduced in this paper. I. INTRODUCTION operator of finite-size systems. Some approximations we can discuss in a unified manner with the framework intro- duced in this paper include the Dyson-Neumann series, Our understanding of the dynamics of quantum many- [44{46] the Magnus expansion,[47{51] Fer's series,[50{52] particle systems and associated non-equilibrium phenom- Wilcox' series,[50, 51, 53] the rotating frame approxima- ena has seen rapid growth in recent years [1, 2], which tion, [13] and flow equation methods.[54] has resulted from advances in theory|especially Flo- quet systems [3{15]|and experiment, particularly in the The focus of this paper is on the flow equations for preparation of and characterization of non-equilibrium couplings (which we will introduce shortly) in a Hamil- states.[2, 16{22] As the field of non-equilibrium quantum tonian and their relation to the various approximations systems expands, an increasing amount of effort is being mentioned above. We note that important results were devoted to the study of such systems, driven in part by obtained in prior work by flow equation methods in the the wealth of novel phenomenology that the time domain equilibrium case, [55, 56] and the non-equilibrium case permits. For cold atoms trapped in optical lattices, time- [57] making use of the construct of a Sambe space.[58] In dependent driving has enabled [2, 23] coherent control this work, we find an analogue to Hamilton-Jacobi the- of tunneling [24], induction of phase transitions [25, 26], ory for the time-evolution operator that allows one to generation of effective magnetic fields [27], and the mea- make clear connections between various approximation surement of nontrivial topological invariants [28]. Note- schemes. Our formulation takes shape in the form of flow worthy examples from solid state systems include pho- equations for the couplings in the Hamiltonian.[55, 56] toinduced superconductivity, [29, 30] and hidden or oth- Our discussion culminates in a diagram, Fig.1, that erwise inaccessible orders. [31{33] Even more strikingly, interrelates the different approximations and highlights the time domain also allows entirely novel phases, such spots which symmetry suggests can be filled by consid- as time crystals,[34, 35] and non-equilibrium topological ering another limiting flow equation. In this paper, we phases. [9, 36{42] develop this limiting flow equation and find that the ap- A number of numerical and analytical tools have proximations perform as one would have expected from arXiv:1902.07237v3 [cond-mat.stat-mech] 15 Jul 2019 been developed to understand the principal quantity of the order of approximations seen in the diagram in Fig.1. interest|the time evolution operator|which mathemat- We display a table, TableI, that makes clear to the reader ically is a time-ordered exponential. It is not practical to the range of validity for each approximation, along with summarize all known methods to calculate this quantity, its advantages and shortcomings. In particular, we find so we will set our focus on analytically accessible approx- that the flow equations obtained in Ref.[54] are useful imations that can be used for arbitrary forms of time- when truncated like Wegner's, [55, 56] and the approxi- dependence. This restriction will therefore exclude the mation we develop in this paper (by completing our flow vast literature on numerical methods, and a recently de- chart) fills in a gap for an analytically accessible approx- veloped exact method [43] to calculate the time-evolution imation in the intermediate time regime (or equivalently 2 for a time-periodic system in the intermediate frequency If one repeatedly reinserts the left side into the right side regime). Our results provide a comprehensive picture of of the equation one obtains a series in powers of H(t). the approximations to the time-evolution operator that Truncated at second order, this series is given as, can be used for time-periodic quantum many-particle sys- Z t Z t Z t1 0 0 3 tems. The approximation we develop here is particularly U(t) = 1H −i dt H(t )− dt1 dt2H(t1)H(t2)+O(H ): useful when kV k kH0k ∼ !. That is, when the system 0 0 0 is subject to a drive V that is weak compared to the static (4) Neglecting higher order terms in H(t) is valid if H(t) is part of the Hamiltonian H0 but when the static part H0 is not negligible when compared to the drive frequency small compared to ∼ 1=t. In other words, such an expan- !. This regime is relevant to e.g. cavity-QED appli- sion is restricted to sufficiently short times. In addition, cations [59, 60] (where strongly interacting photons are truncating the series destroys its unitarity, which is a se- often subject to weak time-periodic drives), and weakly rious drawback. The loss of unitarity occurs already at driven cold atom quantum ratchets.[61] first order in H(t). Our paper is organized as follows. In Sec.II, we give If one is interested in the evolution of eigenstates of a short introduction to each of the common approxi- a constant Hamiltonian H0, it is advantageous to split −iH0t mations mentioned above. In Sec.III, we establish re- U(t) = UI (t)e because in this case the second factor −iH t iφ lations between the different approximations. In Sec.IV only leads to a phase e 0 j i = e j i. For expecta- we present the new approximation and a diagram of re- tion values of an operator O^ the phase factor disappears ^ y ^ y ^ lations in Fig.1 that were suggested by symmetry. In hOi(t) = h jU (t)OU(t)j i = h jUI (t)OUI (t)j i and UI Sec.V we compare the different approximations for an encodes all relevant information. It now fulfills Ising model. In section VI we calculate the l2 distance Z t 0 0 2 between the exact and the approximate time-evolution UI (t) = 1H − i dt VI (t ) + O(VI ); (5) operators. Finally, in Sec.VII we present our conclusions. 0 iH0t −iH0t A few technical details are relegated to the appendices. where VI = e (H(t) − H0)e . Such a procedure is called the interaction picture (hence the \I" subscript). The expansion for UI (t) is often used in setting up Feyn- II. SUMMARY OF COMMON man diagrams for scattering problems because for small APPROXIMATIONS TO U(t) perturbations V (t) the procedure works well enough to approximate the S-matrix, S = UI (t = 1). The time-evolution operator U(t) fulfills i@ U(t) = H(t)U(t); U(0) = 1 ; (1) t H B. Magnus expansion where 1H is the identity on the Hilbert space of the Hamiltonian, H. We have set Planck's constant, ~ = 1. The broken unitarity in Dyson's approach is a seri- Eq.(1) can be solved by a simple matrix exponential if ous drawback because spurious terms may appear in [H(t1);H(t2)] = 0 for all t1 and t2. However, when calculations|a problem Magnus[47{51] solved. His way [H(t1);H(t2)] 6= 0 Eq.(1)becomes more complicated and of dealing with this issue was by making the ansatz one has to concatenate matrix exponentials for infinites- U(t) = eΩ(t) for the time-evolution operator and search- imal time steps dt, ing for anti-Hermitian Ω instead of U(t). Inserting this U(t) = e−idtH(t−dt)e−idtH(t−2dt) ··· e−idtH(0): (2) ansatz into Eq.(1) and using the general expression for the derivative of the exponential map he found that, Such a procedure, which sometimes is called Trot- terization, reproduces the formal solution: U(t) = dΩ(t) adΩ −i R t dt0H(t0) = i ad H(t); (6) T e 0 in the limit dt ! 0, where T is the time- dt e Ω − 1 order operator and the expression is called a time-ordered where the shorthand adΩ = [Ω;:] was used.