Recent studies have shown that Trotter errors are highly initial-state dependent and that standard upper bounds often substantially overestimate them. However, the mechanism underlying anomalously small Trotter errors and a systematic route to identifying error-resilient states remain unclear. Using interaction-picture perturbation theory, we derive an analytical expression for the leading-order Trotter error in the eigenbasis of the Hamiltonian. Our analysis shows that initial states supported on spectrally commensurate energy ladders exhibit strongly suppressed error growth together with persistent Loschmidt revivals. We refer to such states as Trotter scars. To identify such states in practice, we further introduce a general variational framework for finding error-minimizing initial states for a given Hamiltonian. Applying this framework to several spin models, we find optimized states whose spectral support and dynamical behavior agree with the perturbative prediction. Our results reveal the spectral origin of Trotter-error resilience and provide a practical strategy for discovering error-resilient states in digital quantum simulation.
Recent studies have shown that Trotter errors are highly initial-state dependent and that standard upper bounds often substantially overestimate them. However, the mechanism underlying anomalously small Trotter errors and a systematic route to identifying error-resilient states remain unclear. Using interaction-picture perturbation theory, we derive an analytical expression for the leading-order Trotter error in the eigenbasis of the Hamiltonian. Our analysis shows that initial states supported on spectrally commensurate energy ladders exhibit strongly suppressed error growth together with persistent Loschmidt revivals. We refer to such states as Trotter scars. To identify such states in practice, we further introduce a general variational framework for finding error-minimizing initial states for a given Hamiltonian. Applying this framework to several spin models, we find optimized states whose spectral support and dynamical behavior agree with the perturbative prediction. Our results reveal the spectral origin of Trotter-error resilience and provide a practical strategy for discovering error-resilient states in digital quantum simulation.
Trotter errors in digitized quantum dynamics arise from approximating time-ordered evolution under noncommuting Hamiltonian terms with a product formula. In the adiabatic regime, such errors are known to exhibit long-time self-healing [Phys. Rev. Lett. \textbf{131}, 060602 (2023)], where discretization effects are effectively suppressed. Here we show that self-healing persists at finite evolution times once nonadiabatic errors induced by finite-speed ramps are compensated. Using counterdiabatic driving to cancel diabatic transitions and isolate discretization effects, we study both noninteracting and interacting spin models and characterize the finite-time scaling with the Trotter steps and the total evolution time. In the instantaneous eigenbasis of the driven Hamiltonian, the leading digital error maps to an effective harmonic perturbation whose dominant Fourier component yields an analytic upper bound on the finite-time Trotter error and reveals the phase-cancellation mechanism underlying self-healing. Our results establish finite-time self-healing as a generic feature of digitized counterdiabatic protocols, clarify its mechanism beyond the long-time adiabatic limit, and provide practical guidance for high-fidelity state preparation on gate-based quantum processors.
Trotter product formulas are a leading approach for Hamiltonian simulation on quantum computers, yet their practical performance has remained difficult to assess due to the challenge of accurately estimating the Trotter error. In this work, we develop new theoretical results, algorithms, and software tools that advance the state-of-the-art in Trotter error estimation by orders of magnitude in both scale and accuracy. On the theoretical side, we prove that in the asymptotic limit the error of a product formula depends on the diagonal elements of the Baker-Campbell-Hausdorff (BCH) error operator in the eigenbasis of the Hamiltonian, rather than its full spectral norm -- yielding an improved scaling for Hamiltonian simulation using product formulas. On the algorithmic side, we introduce a compact representation of the BCH expansion that reduces the number of commutators from O(n3) to O(n) for second-order, and from O(n5) to O(n2) for fourth-order formulas on n fragments, complemented by an importance sampling scheme to further reduce the computational cost. We provide implementations of these techniques in software and demonstrate their power on two applications: (i) X-ray absorption spectroscopy of an electronic Hamiltonian (Li4Mn2O) at up to 56 qubits using tensor networks; and (ii) vibronic dynamics of naphthalene at over 100 qubits using ML-MCTDH, where we find that naive analytical bounds overestimate the required number of Trotter steps by nearly five orders of magnitude. Our framework enables, for the first time, the accurate estimation of Trotter error at practically relevant system sizes, providing a foundation for fair algorithmic comparisons and rational design of product formulas.
Quantum many-body scars are an intriguing dynamical regime in which quantum systems exhibit coherent dynamics and long-range correlations when prepared in certain initial states. We use this combination of coherence and many-body correlations to benchmark the performance of present-day quantum computing devices by using them to simulate the dynamics of an antiferromagnetic initial state in mixed-field Ising chains of up to 19 sites. In addition to calculating the dynamics of local observables, we also calculate the Loschmidt echo and a nontrivial unequal-time connected correlation function that witnesses long-range many-body correlations in the scarred dynamics. We find coherent dynamics to persist over up to 40 Trotter steps even in the presence of various sources of error. To obtain these results, we leverage a variety of error-mitigation techniques including noise tailoring, zero-noise extrapolation, dynamical decoupling, and physically motivated postselection of measurement results. Crucially, we also find that using pulse-level control to implement the Ising interaction yields a substantial improvement over the standard controlled-not-based compilation of this interaction. Our results demonstrate the power of error-mitigation techniques and pulse-level control to probe many-body coherence and correlation effects on present-day quantum hardware.
Simulating quantum dynamics beyond the reach of classical computers is one of the main envisioned applications of quantum computers. The most promising quantum algorithms to this end in the near term are the simplest, which use the Trotter formula and its higher-order variants to approximate the dynamics of interest. The approximation error of these algorithms is often poorly understood, even in the most basic and topical cases where the target Hamiltonian decomposes into two realizable terms: H=H_{1}+H_{2}. Recent studies have reported anomalously low approximation error with unexpected scaling in such cases, which they attribute to quantum interference between the errors from different steps of the algorithm. Here, we provide a simpler picture of these effects by relating the Trotter formula to its second-order variant for such H=H_{1}+H_{2} cases. Our method generalizes state-of-the-art error bounds without the technical caveats of prior studies, and elucidates how each part of the total error arises from the underlying quantum circuit. We compare our bound to the true error numerically, and find a close match over many orders of magnitude in the simulation parameters. Our findings further reduce the required circuit depth for the least experimentally demanding quantum simulation algorithms, and illustrate a useful method for bounding simulation error more broadly.
Accurate digital quantum simulation at long times is limited by the accumulation of errors inherent to approximate simulation. Here we introduce RL-Trotter, a reinforcement-learning framework that treats unavoidable approximation errors as resources for error correction rather than merely imperfections to suppress. We show that low-dimensional information from conservation laws, such as the energy and energy variance, provides a sufficient learning signal to guide the agent, which learns to adapt a single scalar---the next Trotter step size---without access to the target wave function. By optimizing the entire long-time evolution rather than individual steps, RL-Trotter discovers self-correcting sequences in which later errors compensate for those accumulated earlier, increasing the accuracy of the long-time dynamics. The learned policies are intrinsically robust to measurement noise, substantially reducing measurement overhead. They also generalize to previously unseen, physically similar initial states and transfer from small, classically simulable systems to systems an order of magnitude larger. This enables a practical protocol based on classical pretraining followed by direct deployment or limited fine-tuning on quantum hardware. Our results establish a broader perspective for quantum algorithms: errors in approximate evolution can be orchestrated into resources for accurate and resource-efficient quantum dynamics.
Abstract We investigated the stability of quantum many-body scars under perturbations using a PXP model. Using numerical simulations, we analyzed the fidelity and correlation functions of the system to track its evolution and assess the persistence of revivals. We explored two distinct types of perturbations: imperfections in the initial state, which naturally arise in experimental settings, and disorder in the Hamiltonian, which was introduced as a computational tool to probe the system’s response to external disturbances. We found that while certain initial state modifications can partially stabilize oscillations by restricting the system’s accessible dynamics, others accelerate thermalization by increasing overlap with the full eigenstate spectrum. Similarly, disorder in the Hamiltonian gradually suppresses revivals and enhances thermalization, but scars demonstrate robustness for small perturbations. Our results provide insights into the conditions under which quantum many-body scars persist, offering valuable implications for experimental realizations and controlled quantum dynamics.