Abstract:The Lie-Trotter formula, together with its higher-order generalizations, provides a direct approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of prior approaches based on truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, including simulations of second-quantized plane-wave electronic structure, $k$-local Hamiltonians, rapidly decaying power-law interactions, clustered Hamiltonians, the transverse field Ising model, and quantum ferromagnets, nearly matching or even outperforming the best previous results. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound as a byproduct. Our analysis reproduces known tight bounds for first- and second-order formulas. Our higher-order bound overestimates the complexity of simulating a one-dimensional Heisenberg model with an even-odd ordering of terms by only a factor of $5$, and is close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both asymptotic scaling and constant prefactor.

Reducing the truncation error in Taylor model multiplication

A Matlab Toolbox for Positive Fractional Time Derivative Modeling of Arbitrarily Frequency-Dependent Viscosity

Aggressively Truncated Taylor Series Method For Accurate Computation Of Exponentials Of Essentially Nonnegative Matrices

Model Reduction of Parametric Differential-Algebraic Systems by Balanced Truncation

Error bounds for frequency weighted balanced truncation and relative error model reduction

Automatically Bounding the Taylor Remainder Series: Tighter Bounds and New Applications

Explicit and implicit error inhibiting schemes with post-processing

Enhancing Interpolation and Approximation Error Estimates Using a Novel Taylor-like Formula

Quality of Approximate Balanced Truncation

Measurable Taylor's Theorem: An Elementary Proof

Model Reduction Using the Balanced Truncation Method and the Padé Approximation Method

A new second order Taylor-like theorem with an optimized reduced remainder

Balanced Truncation Reduced Models of Quadratic-Bilinear Systems in Time Interval

Truncation Error Calculation Based on Richardson Extrapolation for Variable-Step Collaborative Simulation

Taylor-Series Expansion Based Numerical Methods: A Primer, Performance Benchmarking and New Approaches for Problems with Non-smooth Solutions

Modeling still matters: a surprising instance of catastrophic floating point errors in mathematical biology and numerical methods for ODEs

Improving the numerical stability of fast matrix multiplication

Exponentially Reduced Circuit Depths Using Trotter Error Mitigation

A Theory of Trotter Error

Truncation Error of the Network Simulation Method: Chaotic Dynamical Systems in Mechanical Engineering

Approximation of the derivatives beyond Taylor expansion