Abstract
The algorithmic error of digital quantum simulations is usually explored in terms of the spectral norm distance between the actual and ideal evolution operators. In practice, this worst-case error analysis may be unnecessarily pessimistic. To address this, we develop a theory of average-case performance of Hamiltonian simulation with random initial states. We relate the average-case error to the Frobenius norm of the multiplicative error and give upper bounds for the product formula (PF) and truncated Taylor series methods. As applications, we estimate average-case error for the digital Hamiltonian simulation of general lattice Hamiltonians and -local Hamiltonians. In particular, for the nearest-neighbor Heisenberg chain with spins, the error is quadratically reduced from in the worst case to on average for both the PF method and the Taylor series method. Numerical evidence suggests that this theory accurately characterizes the average error for concrete models. We also apply our results to error analysis in the simulation of quantum scrambling.
- Received 16 November 2021
- Revised 1 November 2022
- Accepted 21 November 2022
DOI:https://doi.org/10.1103/PhysRevLett.129.270502
© 2022 American Physical Society