Computation of operator exponentials using the Dunford-Cauchy integral

We consider an n-qubit quantum system with a Hamiltonian, defined by an expansion in the Pauli basis, and propose a new algorithm for classical computing the exponential of the Hamiltonian. The algorithm is based on the representation of the exponential by the Dunford-Cauchy integral, followed by an efficient computation of the resolvent, and is suitable for Hamiltonians that are sparse in the Pauli basis. The practical efficiency of the algorithm is demonstrated by two illustrative examples.