Rigorous Maximum Likelihood Estimation for Quantum States

TL;DR

This paper introduces a scalable, low-memory maximum likelihood estimation method for quantum state tomography that efficiently recovers large quantum states without requiring explicit density matrix formation.

Contribution

It reformulates MLE using Burer-Monteiro factorization, eliminating constraints analytically, and develops a low-memory algorithm suitable for large quantum states with competitive performance.

Findings

Able to recover 20-qubit states on a laptop in under 5 hours

Requires O(d log d) complexity per iteration for fixed rank

Demonstrates competitive performance with state-of-the-art methods

Abstract

Existing quantum state tomography methods are limited in scalability due to their high computation and memory demands, making them impractical for recovery of large quantum states. In this work, we address these limitations by reformulating the maximum likelihood estimation (MLE) problem using the Burer-Monteiro factorization, resulting in a non-convex but low-rank parameterization of the density matrix. We derive a fully unconstrained formulation by analytically eliminating the trace-one and positive semidefinite constraints, thereby avoiding the need for projection steps during optimization. Furthermore, we determine the Lagrange multiplier associated with the unit-trace constraint a priori, reducing computational overhead. The resulting formulation is amenable to scalable first-order optimization, and we demonstrate its tractability using limited-memory BFGS (L-BFGS). Importantly, we also propose a low-memory version of the above algorithm to fully recover certain large quantum states with Pauli-based POVM measurements. Our low-memory algorithm avoids explicitly forming any density matrix, and does not require the density matrix to have a matrix product state (MPS) or other tensor structure. For a fixed number of measurements and fixed rank, our algorithm requires just $\mathcal{O}(d \log d)$ complexity per iteration to recover a $d \times d$ density matrix. Additionally, we derive a useful error bound that can be used to give a rigorous termination criterion. We numerically demonstrate that our method is competitive with state-of-the-art algorithms for moderately sized problems, and then demonstrate that our method can solve a 20-qubit problem on a laptop in under 5 hours.