We study universal aspects of thermalization induced by Trotterization, a procedure routinely used in gate-based quantum computation. We use the reduced-Bardeen-Cooper-Schrieffer model-quantum integrable with a classically integrable mean-field limit-where the effects of Trotter chaos are expected to be particularly stark. The resulting Trotterized chaotic dynamics is characterized by its Lyapunov spectrum and rescaled Kolmogorov-Sinai entropy. The chaos quantifiers depend on the Trotterization time step τ. We observe a Trotter transition at a finite step value τ_{c}≈sqrt[N]. While the dynamics is weakly chaotic for time steps τ≪τ_{c}, the regime of large Trotterization steps is characterized by short temporal correlations. We derive two different scaling laws for the two different regimes by numerically fitting the maximum Lyapunov exponent data. The scaling law of the large τ limit agrees well with the one derived from the kicked top map. Beyond its relevance to current quantum computers, our work opens other directions-such as probing observables like the Loschmidt echo, which lie beyond standard mean-field description--across the Trotter transition we uncover.
We study universal aspects of thermalization induced by Trotterization, a procedure routinely used in gate-based quantum computation. We use the reduced-Bardeen-Cooper-Schrieffer model-quantum integrable with a classically integrable mean-field limit-where the effects of Trotter chaos are expected to be particularly stark. The resulting Trotterized chaotic dynamics is characterized by its Lyapunov spectrum and rescaled Kolmogorov-Sinai entropy. The chaos quantifiers depend on the Trotterization time step τ. We observe a Trotter transition at a finite step value τ_{c}≈sqrt[N]. While the dynamics is weakly chaotic for time steps τ≪τ_{c}, the regime of large Trotterization steps is characterized by short temporal correlations. We derive two different scaling laws for the two different regimes by numerically fitting the maximum Lyapunov exponent data. The scaling law of the large τ limit agrees well with the one derived from the kicked top map. Beyond its relevance to current quantum computers, our work opens other directions-such as probing observables like the Loschmidt echo, which lie beyond standard mean-field description--across the Trotter transition we uncover.
We investigate many-body dynamics where the evolution is governed by unitary circuits through the lens of “Krylov complexity,” a recently proposed measure of complexity and quantum chaos. We extend the formalism of Krylov complexity to unitary circuit dynamics and focus on Floquet circuits arising as the Trotter decomposition of Hamiltonian dynamics. For short Trotter steps the results from Hamiltonian dynamics are recovered, whereas a large Trotter step results in different universal behavior characterized by the existence of local : operators with vanishing autocorrelation functions, as exemplified in dual-unitary circuits. These operators exhibit maximal complexity growth, act as a memoryless bath for the dynamics, and can be directly probed in current quantum computing setups. These two regimes are separated by a crossover in chaotic systems. Conversely, we find that free integrable systems exhibit a nonanalytic transition between these different regimes, where maximally ergodic operators appear at a critical Trotter step.
Digital quantum simulation (DQS) is one of the most promising applications of quantum computing. Trotterised DQS models quantum dynamics by decomposing the evolution into a sequence of gate-based Trotter steps. Using larger Trotter steps reduces the required gate count, but it is now understood that the Trotterisation step size cannot be increased beyond a certain threshold value, beyond which there is a sudden breakdown in simulation accuracy, due to the onset of quantum chaos. To date, this behaviour has only been studied for closed systems not subjected to the environmental noise and decoherence to which real-world quantum computers are highly sensitive. In this work, we study a range of experimentally relevant open-system quantum simulator models and show that Trotterisation thresholds persist even when open-system effects are included, even in the long-time, steady-state limit, once all coherence has been lost. We show that open DQS can accurately reproduce the dynamics of open target systems prior to the threshold, but still exhibits a sudden performance breakdown beyond it. We introduce new, statistical techniques to rigorously characterise the onset of quantum chaos beyond the Trotterisation threshold for both approximately unitary and genuine dissipative quantum chaos, based on random matrix theory (RMT), building on techniques developed for closed systems. Finally, we demonstrate that dissipative quantum chaos emerges when the dissipation is in balance with the presence of underlying unitary quantum chaos; the dissipation must be strong enough to induce genuinely dissipative behaviour, yet not so strong as to suppress the underlying chaotic structure.
The Trotter-Suzuki decomposition is a promising avenue for digital quantum simulation (DQS), approximating continuous-time dynamics by discrete Trotter steps of duration τ. Recent work suggested that DQS is typically characterized by a sharp Trotter transition: when τ is increased beyond a threshold value, approximation errors become uncontrolled at large times due to the onset of quantum chaos. Here, we contrast this picture with the case of integrable DQS. We focus on a simple quench from a spin-wave state in the prototypical XXZ Heisenberg spin chain, and study its integrable Trotterized evolution as a function of τ. Because of its exact local conservation laws, the system does not heat up to infinite temperature and the late-time properties of the dynamics are captured by a discrete generalized Gibbs ensemble (dGGE). By means of exact calculations we find that, for small τ, the dGGE depends analytically on the Trotter step, implying that discretization errors remain bounded even at infinite times. Conversely, the dGGE changes abruptly at a threshold value τ_{th}, signaling a novel type of Trotter transition. We show that the latter can be detected locally, as it is associated with the appearance of a nonzero staggered magnetization with a subtle dependence on τ. We highlight the differences between continuous and discrete GGEs, suggesting the latter as novel interesting nonequilibrium states exclusive to digital platforms.
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 of previous claims, there is no Trotterization transition in nonintegrable many-body quantum systems. We demonstrate our results on a one-dimensional (1D) kicked Ising model, as well as on a 1D kicked XX model and a 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.
Digital quantum simulation (DQS) is one of the most promising paths for achieving first useful real-world applications for quantum processors. Yet even assuming rapid progress in device engineering and development of fault-tolerant quantum processors, algorithmic resource optimisation will long remain crucial to exploit their full power. Currently, Trotterisation provides state-of-the-art resource scaling. And recent theoretical studies of Trotterised Ising models suggest that even better performance than expected may be possible up to a distinct breakdown threshold in empirical performance. Here, we study multiple paradigmatic DQS models with experimentally realisable Trotterisations, and evidence the universality of a range of Trotterisation performance behaviours, including not only the threshold, but also new features in the pre-threshold regime that is most important for practical applications. In each model, we observe a distinct Trotterisation threshold shared across widely varying performance signatures; we further show that an onset of quantum chaotic dynamics causes the performance breakdown and is directly induced by digitisation errors. In the important pre-threshold regime, we are able to identify new distinct regimes displaying qualitatively different quasiperiodic performance behaviours, and show analytic behaviour for properly defined operational Trotter errors. Our results rely crucially on diverse new analytical tools, and provide a previously missing unified picture of Trotterisation behaviour across local observables, the global quantum state, and the full Trotterised unitary. This work provides new insights and tools for addressing important questions about the algorithm performance and underlying theoretical principles of sufficiently complex Trotterisation-based DQS, that will help in extracting maximum simulation power from future quantum processors.
We propose a novel framework to characterize the thermalization of many-body dynamical systems close to integrable limits using the scaling properties of the full Lyapunov spectrum. We use a classical unitary map model to investigate macroscopic weakly nonintegrable dynamics beyond the limits set by the KAM regime. We perform our analysis in two fundamentally distinct long-range and short-range integrable limits which stem from the type of nonintegrable perturbations. Long-range limits result in a single parameter scaling of the Lyapunov spectrum, with the inverse largest Lyapunov exponent being the only diverging control parameter and the rescaled spectrum approaching an analytical function. Short-range limits result in a dramatic slowing down of thermalization which manifests through the rescaled Lyapunov spectrum approaching a non-analytic function. An additional diverging length scale controls the exponential suppression of all Lyapunov exponents relative to the largest one.