Abstract When simulating the time evolution of quantum many-body systems on a digital quantum computer, one faces the challenges of quantum noise and of the Trotter error due to time discretization. For certain spin chains, it is possible to discretize the time evolution preserving integrability, so that an extensive set of conserved charges are exactly conserved after discretization. In this work we implement, on real quantum computers and on classical simulators, the integrable Trotterization of the spin- 1 / 2 Heisenberg XXX spin chain. We study how quantum noise affects the time evolution of several conserved charges, and observe the decay of the expectation values. We in addition study the early time behaviors of time evolution, which can potentially be used to benchmark quantum devices and algorithms in the future. We also provide an efficient method to generate the conserved charges at higher orders.
Abstract When simulating the time evolution of quantum many-body systems on a digital quantum computer, one faces the challenges of quantum noise and of the Trotter error due to time discretization. For certain spin chains, it is possible to discretize the time evolution preserving integrability, so that an extensive set of conserved charges are exactly conserved after discretization. In this work we implement, on real quantum computers and on classical simulators, the integrable Trotterization of the spin- 1 / 2 Heisenberg XXX spin chain. We study how quantum noise affects the time evolution of several conserved charges, and observe the decay of the expectation values. We in addition study the early time behaviors of time evolution, which can potentially be used to benchmark quantum devices and algorithms in the future. We also provide an efficient method to generate the conserved charges at higher orders.
Digital quantum simulation relies on Trotterization to discretize time evolution into elementary quantum gates. On current quantum processors with notable gate imperfections, there is a critical trade-off between improved accuracy for finer time steps, and increased error rate on account of the larger circuit depth. We present an adaptive Trotterization algorithm to cope with time dependent Hamiltonians, where we propose a concept of piecewise "conserved" quantities to estimate errors in the time evolution between two (nearby) points in time; these allow us to bound the errors accumulated over the full simulation period. They reduce to standard conservation laws in the case of time independent Hamiltonians, for which we first developed an adaptive Trotterization scheme [H. Zhao et al., Making Trotterization adaptive and energy-self-correcting for NISQ devices and beyond, PRX Quantum 4, 030319 (2023).2691-339910.1103/PRXQuantum.4.030319]. We validate the algorithm for a time dependent quantum spin chain, demonstrating that it can outperform the conventional Trotter algorithm with a fixed step size at a controlled error.
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.
We propose a noise-mitigation quantum simulation strategy for near-term quantum devices based on Quantum Circuit Learning (QCL), which is in particular effective for integrable quantum spin chains. The method trains a shallow variational circuit to approximate a deeper time-evolution circuit by learning the conserved charges and only a small amount of dynamical information in the system. Under realistic noise models, the learned circuit maintains both conserved quantities and dynamical observables significantly closer to their true values than the noisy simulation of the original circuit. We demonstrate, on small-scale prototypes, that QCL can act as an effective, physics-informed error mitigation strategy, producing shorter, more robust circuits without exponential sampling overhead.
We discuss a general procedure to construct an integrable real-time Trotterization of interacting lattice models. As an illustrative example, we consider a spin-1/2 chain, with continuous time dynamics described by the isotropic (XXX) Heisenberg Hamiltonian. For periodic boundary conditions, local conservation laws are derived from an inhomogeneous transfer matrix, and a boost operator is constructed. In the continuous time limit, these local charges reduce to the known integrals of motion of the Heisenberg chain. In a simple Kraus representation, we also examine the nonequilibrium setting, where our integrable cellular automaton is driven by stochastic processes at the boundaries. We show explicitly how an exact nonequilibrium steady-state density matrix can be written in terms of a staggered matrix product ansatz, and we propose quasilocal conservation laws for the model with periodic boundary conditions. This simple Trotterization scheme, in particular in the open system framework, could prove to be a useful tool for experimental simulations of the lattice models in terms of trapped ion and atom optics setups.
We develop a new implementation of the real-time full configuration interaction quantum Monte Carlo method, and apply it to the dynamics of the Heisenberg spin chain and the spectrum of the one-dimensional Fröhlich polaron. This is a general method for the numerical integration of the Schrödinger equation to evolve quantum systems in real time, introduced by Guther et al. as an extension of the imaginary-time version of the algorithm that finds ground-state properties. We introduce a new integrator for improved efficiency, and a scaling procedure for better control of simulations. We study the accuracy of the method, including the effect of discretising time and of the stochastic errors introduced by the Monte Carlo elements of the algorithm. Comparing the method to a hardware quantum simulator for the quench dynamics of a Heisenberg spin chain, we find that for this model the discretisation and stochastic errors grow quickly, limiting the effectiveness of the algorithm for finding accurate evolution over long times. We find that spectral information can be generated successfully using an energy-dependent damping, with either a tilted complex contour time evolution or a strong second-order damping. Specifically, we reproduce the low-energy spectrum of the Fr ̈ohlich polaron model and confirm the existence of a bound excited state.
Quantum algorithms for quantum chemistry and other many-body fermionic systems work by expressing the Hamiltonian in a basis of qubits and fragmenting the Hamiltonian into a sum of products of Pauli operators whose exponentials are easily encoded on a quantum device. Applying the product of exponentials, known as Trotterization, leads to an error associated with the non-commutativity of operators. This error can lead to breaking the symmetries of the Hamiltonian because the fragments are not symmetry conserving in general. Nonetheless, many algorithms for time evolution rely on Trotterization, including time evolution and quantum phase estimation. We show that we can express the Hamiltonian in terms of Hermitian excitation operators which map to sums of commuting Pauli strings for any encoding and conserve symmetries corresponding to Abelian groups of symmetry operators. Symmetries corresponding to non-Abelian groups, on the other hand, are not fully conserved by Trotterized Hermitian excitation operators, so we developed ``operator kirigami'' to cut the sum of non-commuting operators by orthogonal projection and to fold terms together using unitary rotations. We tested pools of operators for small molecules and basis sets, and found that electron number and spin symmetry conserving pools led to greater errors that decreased for larger molecules and were negated with second-order Trotterization. Our work shows the potential for testing quantum computing algorithms on classical computers by adapting tools used in electronic structure theory with conserved symmetries.