Prethermalization, shadowing breakdown, and the absence of Trotterization transition in quantum circuits

TL;DR

This paper investigates the stability and accuracy of quantum circuit simulations over time, revealing that effective energy remains stable due to prethermalization, but other observables do not, and that no Trotterization transition exists in non-integrable systems.

Contribution

It introduces the truncated operator propagator as a tool to analyze Trotterization errors, prethermalization, and discrete time crystals in quantum circuits, challenging previous claims about Trotterization transitions.

Findings

Effective energy is more stable than Trotter errors suggest due to prethermalization.

All other observables deteriorate in the thermodynamic limit.

No Trotterization transition exists in non-integrable many-body quantum systems.

Abstract

One of the premier utilities of present day noisy quantum computers is simulation of many-body quantum systems. We study how long in time is such a discrete-time simulation representative of a continuous time Hamiltonian evolution, namely, a finite time-step introduces so-called Trotterization errors. We demonstrate that the truncated operator propagator (Ruelle-Pollicott resonances) is a powerful tool to that end, as well as to study prethermalization and discrete time crystals, including finding those phenomena at large gate duration. We show that the effective energy is more stable than suggested by Trotter errors -- a manifestation of prethermalization -- while all other observables are not. Even the most stable observable though deteriorates in the thermodynamic limit. Different than in classical systems with the strongest chaos, where the faithfulness time (the shadowing time) can be infinite, in quantum many-body chaotic systems it is finite. A corollary of our results is also that, opposite to previous claims, there is no Trotterization transition in non-integrable many-body quantum systems. We demonstrate our results on a one-dimensional (1d) kicked Ising model, as well as on 1d kicked XX model and 2d kicked Ising model. The truncated propagator is also used to calculate the energy diffusion constant in the tilted-field Ising model with high accuracy.