Research ArticlePHYSICS

Quantum localization bounds Trotter errors in digital quantum simulation

See allHide authors and affiliations

Science Advances  12 Apr 2019:
Vol. 5, no. 4, eaau8342
DOI: 10.1126/sciadv.aau8342


A fundamental challenge in digital quantum simulation (DQS) is the control of an inherent error, which appears when discretizing the time evolution of a quantum many-body system as a sequence of quantum gates, called Trotterization. Here, we show that quantum localization-by constraining the time evolution through quantum interference-strongly bounds these errors for local observables, leading to an error independent of system size and simulation time. DQS is thus intrinsically much more robust than suggested by known error bounds on the global many-body wave function. This robustness is characterized by a sharp threshold as a function of the Trotter step size, which separates a localized region with controllable Trotter errors from a quantum chaotic regime. Our findings show that DQS with comparatively large Trotter steps can retain controlled errors for local observables. It is thus possible to reduce the number of gate operations required to represent the desired time evolution faithfully.


Quantum computers promise to solve certain computational problems exponentially faster than any classical machine (1). A particularly promising application is the solution of quantum many-body problems (2), with large potential impact on quantum chemistry, material science, and fundamental physics. The devices used in this effort can be divided into two major classes: analog quantum simulators, where the Hamiltonian of interest is engineered to mimic the desired quantum many-body physics; and digital quantum simulators (DQSs), where a target time-evolution operator is represented by a sequence of elementary quantum gates. The digital approach is particularly flexible since a universal digital quantum simulator can be freely programmed to simulate the unitary evolution of any many-body Hamiltonian with local interactions (Fig. 1A) (3). Recent experiments have demonstrated remarkable progress in implementing digital quantum simulation (DQS), e.g., by simulating simple molecules in quantum chemistry (46), condensed-matter models (712), and lattice gauge theories (13).

Fig. 1 Trotterized time evolution and resulting error on local observables.

(A) Gate sequence for the digital quantum simulation (DQS) of an Ising model. The desired evolution up to total simulation time t is split into n repeated sequences of length τ = t/n, each decomposed into fundamental quantum gates. The example shows a gate sequence for a four-qubit chain with Ising spin-spin interactions (ZZ) and transverse and longitudinal fields (simulated by single-qubit operations along the X and Z directions on the Bloch sphere). (B) Magnetization dynamics M(t)=N1l=1NSlz(t) in the DQS of the Ising model for N = 20 spins and different Trotter step sizes τ compared to the exact solution. The normalized deviation Δℳ(t)/(hτ)2 with Δℳ(t) = |ℳτ=0(t) − ℳ(t)| from the ideal dynamics shows a collapse of the error dynamics ℳτ=0(t) for sufficiently small τ.

The working principle of DQS is as follows. Suppose that the target Hamiltonian H=l=1MHl can be decomposed into M terms whose time evolution operators Ul(t) = exp(−itHl) can be implemented on the considered quantum computing device. Using the Suzuki-Trotter formula, the full time-evolution operator U(t) = exp(−itH) can be approximated by discretizing it into nN repetitions of the fundamental gates UiU(n)(t)=[U1(tn)U2(tn)UM(tn)]n(1)

This Trotterization comes inherently with an error that can be rigorously upper bounded via the accuracy of the global unitary time-evolution operator (3)U(t)U(n)(t)=t2nl>m=1M[Hl,Hm]+ϵ(2)

Here, ϵ subsumes terms of order t3/n2 and higher. According to Eq. 2, for the lowest-order corrections, the Trotter error on the full time-evolution operator may grow quadratically with total simulation time t and (in generic quantum many-body systems) linearly in the number of simulated degrees of freedom N. It is possible to improve this upper bound, but an error bound that scales less than linear in t is not possible if one is concerned with the entire unitary operator (14). Although the polynomial scaling with both t and N is efficient in a computational complexity sense, it poses a substantial challenge for practical computations (15, 16), seemingly preventing current technology from simulating all but small instances. As we show in this article, these generic bounds on the global many-body wave function overestimate by far the actual error on local observables such as magnetizations or low-order correlation functions. For example, in the DQS of a quantum Ising chain, the deviation of the magnetization dynamics from the ideal evolution can be considerably smaller and remain bounded even at long times (see Fig. 1B). It is the purpose of this article to explain this observation from physical grounds and thus assign a physical interpretation to Trotter errors.

We achieve this by linking Trotter errors to quantum localization. Localization is a ubiquitous phenomenon with many facets. Initially, it has been introduced to understand the absence of transport in systems of free particles with disorder (17). Since then, the concept has been generalized to various contexts such as many-body localization in Hilbert space as absence of quantum ergodicity (18) or energy localization in periodic time-dependent quantum many-body systems as absence of heating in continuously driven systems (19). As we show here, there occurs a related localization in Hilbert space at small Trotter steps that bounds time-discretization errors on local observables occurs. Our aim is to isolate the role of this universal error source, for which we focus in the following mostly on an idealized setting. We discuss the interplay of the Trotter errors with other, platform-dependent error sources in the concluding remarks and in the Supplementary Materials.


Trotter sequences as Floquet systems

In this work, we interpret the Trotterized evolution as a periodically time-dependent quantum many-body system with a period τ = t/n (see Fig. 1). The desired stroboscopic dynamics is therefore governed by an associated Floquet Hamiltonian HF, which we define for later convenience in the following formeiHFτ=U1(τ)U2(τ)UM(τ)(3)

The starting point of our considerations is an analytical expression for HF in the limit of sufficiently small Trotter steps τHF=H+iτ2l>m[Hl,Hm]+O(τ2)(4)

This form, which can be obtained from Eq. 3 via a Magnus expansion, quantifies the Trotterization error on a Hamiltonian level. There remain, however, two fundamental questions that we aim to address in this work: (i) What is the radius of convergence τ* of this expansion? [Mathematically rigorous bounds for the convergence radius of the Magnus expansion do exist, but their applicability to generic quantum many-body systems is not often evident (19, 20).] (ii) What is the influence of corrections to H that appear in HF on the long-time dynamics of local observables? Recent theoretical predictions for heating in generic quantum many-body systems subject to a periodic drive might leave a rather pessimistic impression (2123). We show in this work that the errors on local observables can nevertheless be controlled for all practical purposes. Throughout this work, we adopt the notion of local observables to be operators that include a bounded number of constituents. This definition includes local order parameters and typical correlation functions. These are not only generically the measures to describe the properties of physical systems, they are also the quantities that can be experimentally measured in a scalable way. Moreover, our general argumentation holds for physical Hamiltonians with few-body interactions, which for almost all of our discussion may even be taken to be long-ranged, except if otherwise stated.

Benchmark model: Quantum Ising chain

In the following, we illustrate our discussion with a generic, experimentally relevant model, the quantum Ising chain with Hamiltonian H = HZ + HX, with HZ=Jl=1N1SlzSl+1z+hl=1NSlz and HX=gl=1NSlx. Here, Slγ (γ = x, y, z) denotes spin-1/2 operators at lattice sites l = 1,…, N. These models are paradigmatic workhorses for DQS platforms such as nuclear magnetic resonance (24), trapped ions (7), and superconducting qubits (25). As an initial state, we choose |ψ0〉 = ⊗ l| ↑ 〉l, which can be prepared with high fidelity (7, 25, 26). In the remainder, we use the parameters h/J = g/J = 1 and choose to measure times in the characteristic scale h−1. For details about the simulations including the used gate sequences, see Materials and Methods. Although we focus on this model, our findings also apply to other systems and thus seem generic (see also the Supplementary Materials where we provide a similar analysis for the lattice Schwinger model).

Quantum many-body chaos threshold

As the central result of this work, we connect Trotter errors in DQS with a threshold separating a many-body quantum chaotic region from a localized regime, thus linking the intrinsic accuracy of DQS with a quantum many-body phenomenon. For that purpose, we first investigate the inverse participation ratio (IPR)IPR =vpv2,pv=|ϕv|ψ0|2(5)with |ϕv〉 denoting a full set of eigenstates of the Floquet Hamiltonian HF. The IPR measures the localization properties of the state |ψ0〉 in the eigenbasis |ϕv〉, which is well studied also in the single-particle context (27). In a quantum chaotic delocalized regime, |ψ0〉 is scrambled across the full eigenbasis implying a uniform distribution pvD1, with D as the number of available states in Hilbert space. Since D grows exponentially with the number of degrees of freedom N, we introduce the rate function λD=N1log(D), which exhibits a well-defined thermodynamic limit. Analogously, we define λIPR = −N−1 log(IPR). In Fig. 2A, we show numerical data for the ratio λIPR/λD for the considered benchmark example. For the data in this plot, we take into account the expected leading-order finite-size corrections λD=N1[log(D)log(2)] in the delocalized regime, which can be estimated using random matrix theory (28). As one can see, there appears a sharp threshold separating a quantum chaotic regime at large Trotter steps, where λIPR tends to λD with increasing system size, from a regular region with λIPR/λD<1.

Fig. 2 Localization and quantum chaos in the Trotterized dynamics of the quantum Ising chain.

(A) Rate function λIPR of the IPR, normalized to the maximally achievable value λD describing uniform delocalization over all accessible states. A sharp threshold as a function of the Trotter step size τ separates a localized regime at small τ from a quantum chaotic regime at large τ. (B) The long-time limit ℱ of the OTO correlator also signals a sharp quantum chaos threshold. ℱ is normalized with respect to ℱ0 = 1/8, the theoretical maximum. Full scrambling is only achieved for large Trotter steps.

A strong fingerprint of quantum chaos can also be found in out-of-time–ordered (OTO) correlators, which quantify how fast quantum information scrambles through a many-body system. A typical OTO correlator is of the formF(t)=V(t)WV(t)W(6)where V(t) denotes the time evolution of the operator V in the Heisenberg picture. While quantum chaos via OTO correlators is conventionally diagnosed by considering a late-time exponential growth for operators V and W with finite support in real space (29), here, we consider the asymptotic long-time value of the extensive operator V=W=N1lSlz (30). We estimate the corresponding long-time limit, ℱ = ℱ (t → ∞), via a stroboscopic average F=limnn1l=1nF(lτ).

In Fig. 2B, we present numerical evidence that this quantity detects the many-body quantum chaos threshold that we have seen in the IPR. There is a clear threshold that separates a localized region at small Trotter steps τ, where ℱ > 0, from a quantum chaotic region at large τ, where ℱ → 0. The vanishing OTO correlator in the many-body quantum chaotic regime can be understood directly from the results obtained for the IPR. Consider the spectral decomposition of a local Hermitian operator V = ∑αλα|α〉〈α|, with λα as the eigenvalues and |α〉 as the eigenvectors of V (for the considered magnetization, these are equivalent to the set of spin configurations). The effective Floquet dynamics yields after n periodsV(nτ)=αv,µλαCναCμα*ei(EνEµ)nτ|ϕνϕµ|(7)with the Floquet quasi-energy corresponding to the eigenstate |ϕν〉 and Cνα = 〈ϕv|α〉.The behavior of the IPR suggests that for all spin configurations, pνα|Cνα|2=D1 is uniformly distributed, such that the amplitudes Cνα are almost structureless and contain only a phase information, Cνα=D1/2eiφvα. After sufficiently many Floquet cycles, this phase information is randomized and scrambled by the unitary evolution, except when v = μ, projecting the operator to the so-called diagonal ensemble (31). Thus, for n → ∞, one obtains V(nτ)D1αλα1. Here, D1αλα=D1TrV is equivalent to the infinite-temperature average, which yields a vanishing value for the considered total magnetization. That is, the operator becomes completely scrambled over the full Hilbert space.

Within the localized phase, the amplitudes Cvα contain more structure than only the phase information, which yields a nonzero value for the OTO correlator. For small systems, such as for N = 10 in Fig. 2B, one can observe additional structures in the crossover region, which vanish for larger N. We attribute these to individual quantum many-body resonances, which can be resolved in small systems but which merge for large N.

Robustness of local observables

While the corrections due to time discretization are weak on a Hamiltonian level, as seen in the Magnus expansion in Eq. 4, there is a priori no guarantee that the long-time dynamics is equally well reproduced. It is, e.g., well known for classical chaotic systems that even weak perturbations can grow quickly in time. Here, we provide numerical evidence that in the localized regime, the dynamics of local observables remains constrained and controlled, even in the long-time limit. This asymptotic long-time dynamics is a worst-case scenario for DQS: When Trotter errors on local observables can be controlled in this limit, so can they on shorter times. In the Supplementary Materials, we also discuss the buildup of Trotter errors on short to intermediate times in more detail.

In Fig. 3A, we show the asymptotic long-time value ℳ of the magnetization, M^(t)=N1lSlz(t). One can observe that the many-body quantum chaos threshold identified in the IPR and OTO correlator has a substantial influence on the long-time Trotter error of local observables such as ℳ(t). For large Trotter steps τ, the magnetization acquires its infinite-temperature value, perfectly consistent with the above analysis of the fully delocalized quantum chaotic phase. However, for small Trotter steps, the error Δℳ relative to the targeted dynamics exhibits a quadratic dependence in τ, as we show in Fig. 3B. The origin of these weak Trotter errors can already be identified from the dynamical trajectories of the magnetization shown in the inset of Fig. 1B, where we plot the error Δℳ(t) for different Trotter steps normalized with respect to (hτ)2. We observe a collapse of trajectories corresponding to different τ, with the overall magnitude of the error remaining bounded in time. This finding suggests that in the localized phase, the discretization error on local observables itself shows regular behavior, in the sense that different perturbation strengths as measured by τ do not yield fast diverging expectation values.

Fig. 3 Trotter errors for local observables in the infinite long-time limit for the Ising model.

Both the magnetization ℳ (A) and simulation accuracy QE (C) exhibit a sharp crossover from a regime of controllable Trotter errors for small Trotter steps τ to a regime of strong heating at larger τ. The dashed line in (A) refers to the desired case of the ideal evolution. The Trotter error exhibits a quadratic scaling at small τ for both the deviation of the magnetization, Δℳ = ℳ − ℳτ=0, (B) and QE (D). The solid lines in (B) and (D) represent analytical results obtained perturbatively in the limit of small Trotter steps τ. These results indicate the controlled robustness of digital quantum simulation against Trotter errors, in the long-time limit and largely independent of N.

Simulation accuracy

In the previous sections, we have provided evidence for a sharp threshold between a delocalized and a localized regime with controllable Trotter errors. We now aim to understand the influence of the regular regions onto the dynamics of local observables. We identify as the underlying reason for the weak Trotter errors a dynamical constraint due to an emergent stroboscopic constant of motion in the effective time-periodic problem, which is the Floquet Hamiltonian HF. Although this integral of motion is different from the desired energy conservation of the target Hamiltonian H, the perturbative expansion in Eq. 4 suggests a close connection. It is therefore natural to quantify the accuracy of DQS by measuring how far the system deviates from the desired constant of motion H viaQE(nτ)Eτ(nτ)E0ET=E0(8)

Here, we have introduced Eτ(nτ) = 〈H(nτ)〉τ and E0 = Eτ → 0(nτ) = 〈ψ0|H0〉, where the subindex τ refers to the used Trotter step for the dynamics. In QE(t), we normalize the errors using the system’s energy at infinite temperature, ET==D1TrH. In the idealized limit τ → 0, where the integral of motion HFH, one has QE(t) = 0. In the opposite limit of large Trotter steps, i.e., in the many-body quantum chaotic region, we expect full delocalization over all eigenstates, yielding QE(t) → 1 in the long-time limit. Thus, QE(t) defines a system-independent measure for the simulation accuracy. From an alternative perspective, QE(t) quantifies heating in the effective periodically driven system, as it has been studied previously in the context of energy localization (19).

In Fig. 3C, we show numerical data for the long-time average QE. Again, we find a sharp threshold between the localized and quantum chaotic regimes. For small Trotter steps, QE acquires only a weak quadratic dependence on τ (see Fig. 3D), yieldingQEQE(t)=(τ/τE)α,ττE(9)with α = 2. While τE depends on the microscopic details of the system, we find from our numerics that there is no notable dependence on N even in the asymptotic long-time limit, with potential corrections in the thermodynamic limit N → ∞ that are discussed further below.

To obtain an analytical understanding for the observations of weak Trotter errors on local observables, let us start by considering the Magnus expansion for the Floquet Hamiltonian in Eq. 4, which quantifies the leading-order corrections due to time discretization on a Hamiltonian level. From our numerical results for QE, we anticipate that the target Hamiltonian H is an almost conserved quantity, which motivates us to study the perturbative corrections to strict energy conservation. Using time-dependent perturbation theory up to second order in the Trotter step size τ, we findQE=qE(hτ)2+O[(hτ)3](10)

The explicit derivation and the final formula for qE are given in Materials and Methods. For the considered parameters, we estimate qE = 0.18. As it can be seen in Fig. 3D, this analytical value matches well the numerical results.

To test whether the errors on other local observables are also controlled by the emergent constant of motion in the localized regime, we exemplarily study the corrections to the targeted magnetization dynamics. From time-dependent perturbation theory we obtain ΔM=m(Jτ)2+O[(Jτ3)] with m = 0.05. This theoretical prediction is again very close to the numerical data (see Fig. 3B). As these findings indicate, in the regular region at small Trotter steps, the discretization error on local observables can be captured by time-dependent perturbation theory in the Trotter step size τ—even in the asymptotic long-time limit.

Our observations give a smaller error on local observables than suggested by general considerations on Floquet dynamics in high-frequency regimes (corresponding to small Trotter steps) (32, 33). In these works, it is shown that there exists a static local Hamiltonian H, which approximates, with exponentially small error on local observables, the stroboscopic Floquet long-time dynamics in the regimes relevant also for the present case. In the derivation of these results, however, it turns out that this Hamiltonian is H, in general, different from H. Our results show that the evolution can be approximated by H itself, as is desired within DQS, while the errors in local observables are now polynomial in τ.


As we have shown, intrinsic Trotter errors in DQS are controllable for local observables, with a sharp threshold separating a localized from a many-body quantum chaotic regime. We have achieved this by identifying the Trotterized time evolution on general grounds with an effective time-periodic Floquet problem. As a consequence, the dynamics is constrained by an emergent conserved quantity given by the Floquet Hamiltonian HF in Eq. 3. While the target Hamiltonian H is not conserved in the Trotterized dynamics, in the localized regime, H remains almost conserved up to perturbative corrections for small Trotter steps τ. This finding does not hold anymore in the many-body quantum chaotic regime, where Trotter errors proliferate and become uncontrollable. While we present data here for a specific model and a specific initial state, our arguments remain general; we find similar properties also for other initial conditions and other model systems such as the recently experimentally realized lattice Schwinger model, which we discuss in the Supplementary Materials. Furthermore, analogous behavior is also found in long-ranged spin models (34). Our numerical studies are based on up to N = 20 qubits, which is within realized and expected size ranges of digital quantum simulators (413, 3538).

For experiments, it is of particular interest to assess the precise value of the threshold scale τ*. Theoretically predicting τ* is in general as difficult as solving the desired time evolution. Nevertheless, one can estimate τ* as follows. Before running an experiment, one can numerically calculate QE for small N, yielding a first estimate on τ*. From this starting point, experiments can find an optimal Trotter step at larger N by decreasing τ until sufficient convergence is reached. Once in the perturbative regime, one can use data at nonzero τ to extrapolate to the ideal dynamics in a well-defined way.

For concrete experimental realizations, it is furthermore of relevance how Trotter errors behave on short to intermediate time scales. From Fig. 1B, one can anticipate that the long-time Trotter error emerges already on rather short time scales. In the Supplementary Materials, we study this transient dynamics in more detail. As we show there for the considered Ising chain, the buildup of the sharp threshold between controllable Trotter errors and quantum chaotic behavior can be observed already on experimentally relevant time scales.

While our results appear to be robust upon increasing the number of degrees of freedom, a quantitative extrapolation to N → ∞ would require the numerical study of larger systems. In this context, recent works have argued that generic periodically driven systems will eventually heat up indefinitely in the thermodynamic limit (2123). This might leave a rather pessimistic impression, but, as we explain now, time discretization errors still remain controllable. Even in the worst-case scenario where such an indefinite heating takes place, the energy growth can still be bounded as long as the Hamiltonian has only short-ranged interactions, via |E(t)E0|Ceτ0/τt for τ ≪ τ0 (22, 23, 32, 33). Here, C denotes a constant of dimension energy squared, and τ0 denotes a constant of dimension time, both of which are independent of N. Thus, for a given total simulation time t, one can ensure a maximum allowed error Δ on the simulation accuracy by QE(t), choosing τ according to τ = τ0/log(ct/Δ) with c = C|ET=∞E0|. In this worst-case scenario, the Trotter step size to reach a given accuracy therefore acquires at most a logarithmic dependence on t but remains independent of N. This is still an exponential improvement over the global wave function bounds such as given in Eq. 2. In practice, since it is tunable via τ, this extremely slow intrinsic heating can always be adjusted such that the associated heating rate is smaller than that of other error sources, such that Trotter errors become unimportant.

Therefore, the accuracy of DQS experiments on local observables is limited mainly by extrinsic error sources. While these may in the future be eliminated by error correction (39, 40), for relevant system sizes to solve many-body problems, full error correction is still out of reach with currently available resources. In the Supplementary Materials, we discuss in detail two typical extrinsic error sources, timing errors on individual gates and slow drifts of gate couplings over various shots of the experiment. The slow drifts turn out to be relatively benign, leading only to an effective average over an ensemble of target Hamiltonians. Individual timing errors, however, induce in the limit of small τ, a time scale beyond which the accuracy of DQS is severely affected. In addition, a realistic implementation on a physical device will suffer from other potential imperfections, many of which can be very device specific. Typical error sources include qubit decoherence and faulty pulses such as imperfect swaps between internal levels. Both of these make it highly preferable to use as few gates as possible. In view of these, our results become particularly relevant: As they show, intrinsic errors in DQS remain controlled even with relatively large Trotter steps. This makes it possible to reach a desired simulation time with a reduced number of gates, thus diminishing the influence of extrinsic errors and enhancing the accuracy in DQS for local observables.


Numerical methods and gate sequences

The numerical data shown in this work was obtained for a quantum Ising chain with the HamiltonianH=HZ+HX(11)whereHZ=Jl=1N1SlzSl+1z+hl=1NSlz,HX=gl=1NSlx(12)

Many of the involved contributions in this model Hamiltonian mutually commute. Therefore, only a small set of elementary quantum gates is required to simulate the Trotterized dynamics. We used the following sequence of two gatesU(1)=U1U2,U1=eiτHZ,U2=eiτHX(13)

For the presented simulations of observables, we computed the real-time evolution for 2 × 104 periods, except otherwise noted, using a Lanczos algorithm with full reorthogonalization. Because, for a finite-size system, observables still show remaining temporal fluctuations, we extracted the asymptotic long-time limit of the presented quantities by performing a stroboscopic time average over the last 104 periods. This large number of Trotter steps is far beyond realistic current-day implementations and serves here only as a worst-case scenario. However, when the Trotter errors on local observables can be controlled even in this idealized limit, one can expect the same to hold true on shorter times relevant for current experiments. In the Supplementary Materials, we illustrate in more detail how the Trotter error builds up on short time scales.

The IPR shown in Fig. 2 can, in principle, be obtained either by exact diagonalization or by use of a dynamical evolution. We chose the latter because it allows us to reach larger systems and is, in principle, an experimentally accessible approach. Dynamically, the IPR can be obtained by a stroboscopic meanIPR =limn1nl=1nPl,Pl=|ψ0|eilHFτ|ψ0|2(14)as one can prove by expanding Pl in the eigenbasis of HF, followed by a summation of the resulting geometric series. Note that Pl is the Loschmidt echo, a common indicator for quantum chaotic behavior in single-particle systems (27).

For the computation of the OTO correlator ℱ(t) defined in Eq. 6, we decompose ℱ(nτ) asF(nτ)=ψ1(nτ)|ψ2(nτ)(15)where the two states|ψ1(nτ)=WeiHFnτVeiHFnτ|ψ0(16)|ψ2(nτ)=eiHFnτVeiHFnτW|ψ0(17)can be obtained from forward and backward evolving the quantum many-body state with appropriate insertions of the W and V operators. Since the backward evolution has to be performed for every Trotter step n, the overall runtime of this approach scales proportional to n2. This limits the accessible total simulation time t = nτ. We used n = 103 for the data shown in Fig. 2B and we performed a stroboscopic average over the last 300 periods to obtain an estimate for the asymptotic long-time value.

Trotter errors on local observables from perturbation theory

As mentioned before, the Trotter errors for local observables can be captured using time-dependent perturbation theory in the limit of sufficiently small τ. In the following, we outline how to obtain the analytical expressions for the coefficients qE and m for QE and ℳ, respectively. First, we consider the simulation accuracy QE and, afterward, the Trotter errors on the magnetization ℳ.

For the derivation of the corrections appearing in QE, we utilize the energy of the target Hamiltonian H, and therefore, the simulation accuracy QE exhibits a substantial overlap with the emergent conserved quantity HFHF(nτ)τ=HF= const.(18)

Here, O(nτ)τ=ψ0|eiHFnτOeiHFnτ|ψ0 denotes the full Trotterized time evolution with Trotter step size τ as in the main text. Moreover, we define the expectation values in the initial state via O=ψ0|O|ψ0 and under the ideal time evolution as O(t)=O(t)τ=0.

To obtain all corrections to the desired order, we first express HF using the Magnus expansion up to second order in the Trotter step sizeHF=H+τC1+τ2C2+O(τ3)(19)withC1=i2[HX,HZ],C2=112[HXHZ,[HX,HZ]](20)

For convenience, we restrict the presentation from now on to a sequence of two elementary gates within one period, as we have for the case of the simulated quantum Ising chain. Using the above expansion for HF in combination with the conservation of HF, one obtains for the energy deviationΔE(nτ)=H(nτ)τH,=τΔC1(nτ)+τ2ΔC2(nτ)(21)whereΔCν(nτ)=CνCν(nτ)τ,ν=1,2(22)

As a next step, we use time-dependent perturbation theory to determine the leading order in τ corrections of ΔCv(nτ). For this purpose, we writeeiHFt=eiHtW(t),W(t)=Tei0tdtV(t)(23)with T denoting the time-ordering prescription andV(t)=eiHtVeiHt,V=τC1+τ2C2(24)

For the corrections to ΔE(nτ) quadratic in τ, we need to perform the time-dependent perturbation theory to first order in τ for C1 and can neglect any τ-dependent contributions for C2.

Let us first consider ΔC1(nτ), which givesΔC1(nτ)=C1C1(nτ)iτ0nτdt[C1(t),C1(nτ)](25)

The time integral can be conveniently evaluated by recognizing thatC1=i2[HX,HZ]=i2[H,HZ](26)since H = HX + HZ and thusC1(t)=12ddtHZ(t)(27)

This givesΔC1(nτ)=C1C1(nτ)iτ2[HZ(nτ)HZ,C1(nτ)](28)

In the limit of n → ∞, we can use the general property that expectation values of operators are governed by the so-called diagonal ensemble (31)O(nτ)nλpλλ|O|λ(29)where pλ = |〈λ|ψ0〉|2 and the set of all |λ〉 denotes the eigenstates for the target Hamiltonian H. Using particular properties of the considered protocol, the above result for ΔC1(nτ) can be simplified considerably. We can use, for example, that C1=0 and [HZ,C1(nτ)]=0, because |ψ0〉 is an eigenstate for HZ, which lastly yieldsΔC1(nτ)nτ4λpλλ|[HZ,[HZ,HX]]|λ(30)

For the contributions to ΔE(nτ) that are second order in τ stemming from ΔC2(nτ), we can restrict to the zeroth order in time-dependent perturbation theory for C2(nτ)τ, i.e., we can replace C2(nτ)τC2(nτ). This yieldsΔC2(nτ)nC2λpλλ|C2|λ(31)

Collecting all contributions, we lastly obtainQE=ΔEET=E0=qE(hτ)2+O[(hτ)3](32)withqE=1J2E0[C2λpλλ|C2|λ14λpλλ|[HZ,[HZ,HX]]|λ](33)given that ET=∞ = 0. This expression can be evaluated using full diagonalization, which provides access to all eigenstates |λ〉. For the considered parameters of our simulations, we find qE = 0.18, which is consistent with the full dynamical calculation in the small Trotter step limit (see Fig. 3D).

For estimating the lowest-order corrections in τ for other observables such as the magnetization ℳ, we can not make direct use of the emergent conserved quantity HF as for the energy of the target Hamiltonian. Still, we can perform time-dependent perturbation theory, which we now have to carry out up to second order. Following the same steps as before, we obtain the following expression for the magnetizationΔM(nτ)=M(nτ)τM(nτ)==τ212[{HZ2(nτ),M(nτ)}EZ2M(t)]+iτ26[C1(nτ)C1,M(nτ)]5τ2120nτdtC1(t)HZ(t)M(nτ)+h.c.(34)

Here, {A, B} = AB + BA denotes the anticommutator, and EZ is given by HZ0〉 = EZ0〉. In the limit n → ∞, we can again use that expectation values can be evaluated in the diagonal ensemble. In addition, the expression involving the time integral can be formally solved using the Lehman representation. Last, we obtainΔM(nτ)nm(hτ)2+O[(hτ)3](35)withm=112J2λpλλ|{HZ2,M}EZ2M|λ16J2λpλRe[λ|[HX,M]HZ|λ]+16J2λ,λpλEλEλRe[λ|[HZ,HX]HZ|λλ|M|λ]+16J2λ,λλ|M|λEλEλRe[CλCλ*λ|[HZ,HX]HZ|λ](36)where Cλ = 〈λ|ψ0〉 and Eλ denotes the eigenenergies of the target Hamiltonian H corresponding to the eigenstate |λ〉. Using full diagonalization, we can again evaluate this expression yielding for our model a value of m = 0.05, which we used in Fig. 3B for the asymptotic small τ prediction and which matches well the result from the full dynamics. Notice that for the presented derivation of the perturbative corrections for the magnetization ℳ, we used explicitly that the initial state is an eigenstate of ℳ. Choosing different observables or different initial conditions might yield linearly in τ contributions as the leading-order corrections.


Supplementary material for this article is available at

Section S1. Temporal buildup of the Trotter error threshold

Section S2. Trotter errors in the Ising model for alternative initial condition

Section S3. Trotter errors in the lattice Schwinger model

Section S4. Imperfections

Fig. S1. Temporal buildup of Trotter errors on transient time scales.

Fig. S2. Simulation accuracy QE for initial Neel state.

Fig. S3. IPR for the DQS of the lattice Schwinger model.

Fig. S4. Trotter errors for the DQS of the lattice Schwinger model.

Fig. S5. Timing errors in the dynamics of the simulation accuracy QE(t) for the Ising model.

Reference (41)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We thank A. Polkovnikov for invaluable discussions in the initial stages of this work as well as L. Sieberer, A. Elben, and T. Olsacher for various comments and suggestions on the manuscript. Funding: Work at Innsbruck was supported by ERC Synergy grant UQUAM and the SFB FoQuS (FWF project no. F4016-N23). M.H. acknowledges support by the Deutsche Forschungsgemeinschaft (DFG) via the Gottfried Wilhelm Leibniz Prize program. P.H. acknowledges support by the DFG Collaborative Research Centre SFB 1225 (ISOQUANT) and the ERC Advanced Grant EntangleGen (project ID 694561). Author contributions: P.Z. proposed the research. P.H. and M.H. performed the analytical calculations and analyzed the data. M.H. performed the numerical simulations. All authors discussed the results and wrote the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data available from the first author upon request.
View Abstract

Navigate This Article