English

$TimeEvolver$: A Program for Time Evolution With Improved Error Bound

Quantum Physics 2022-06-01 v1 Numerical Analysis High Energy Physics - Theory Numerical Analysis Computational Physics

Abstract

We present TimeEvolverTimeEvolver, a program for computing time evolution in a generic quantum system. It relies on well-known Krylov subspace techniques to tackle the problem of multiplying the exponential of a large sparse matrix iHi H, where HH is the Hamiltonian, with an initial vector vv. The fact that HH is Hermitian makes it possible to provide an easily computable bound on the accuracy of the Krylov approximation. Apart from effects of numerical roundoff, the resulting a posteriori error bound is rigorous, which represents a crucial novelty as compared to existing software packages such as ExpokitExpokit (R. Sidje, ACM Trans. Math. Softw. 24 (1) 1998). On a standard notebook, TimeEvolverTimeEvolver allows to compute time evolution with adjustable precision in Hilbert spaces of dimension greater than 10610^6. Additionally, we provide routines for deriving the matrix HH from a more abstract representation of the Hamiltonian operator.

Cite

@article{arxiv.2205.15346,
  title  = {$TimeEvolver$: A Program for Time Evolution With Improved Error Bound},
  author = {Marco Michel and Sebastian Zell},
  journal= {arXiv preprint arXiv:2205.15346},
  year   = {2022}
}

Comments

18 pages, 4 figures

R2 v1 2026-06-24T11:33:37.459Z