A quantum algorithm for estimating the determinant

Abstract

We present a quantum algorithm for estimating the matrix determinant based on quantum spectral sampling. The algorithm estimates the logarithm of the determinant of an $n \times n$ positive sparse matrix to an accuracy $\epsilon$ in time ${\cal O}(\log n/\epsilon^3)$, exponentially faster than previously existing classical or quantum algorithms that scale linearly in $n$. The quantum spectral sampling algorithm generalizes to estimating any quantity $\sum_j f(\lambda_j)$, where $\lambda_j$ are the matrix eigenvalues. For example, the algorithm allows the efficient estimation of the partition function $Z(\beta) =\sum_j e^{-\beta E_j}$ of a Hamiltonian system with energy eigenvalues $E_j$, and of the entropy $ S =-\sum_j p_j \log p_j$ of a density matrix with eigenvalues $p_j$.