Improving the runtime of quantum phase estimation for chemistry through basis set optimization

Quantum phase estimation (QPE) is a promising quantum algorithm for obtaining molecular ground-state energies with chemical accuracy. However, its computational cost, dominated by the Hamiltonian 1-norm $λ$ and the cost of the block encoding, scales at least quadratically with the number of molecular orbitals, making it challenging to incorporate dynamic correlation by enlarging the active space. In this work, we investigate two strategies to mitigate this cost through the optimization of the basis set. First, we investigate whether adjusting the coefficients of Gaussian basis functions can minimize the 1-norm while preserving the accuracy of the ground state energy. Although this method leads to a reduction in the 1-norm up to 10%, this reduction is system-dependent and diminishes with increasing molecular size. Second, we demonstrate that employing a large-basis-set frozen natural orbital (FNO) strategy results in a substantial reduction in QPE resources without compromising accuracy. We study a dataset of 58 small organic molecules and the dissociation curve of N2, and demonstrate that an active space constructed from orbitals derived from larger basis sets captures correlation effects more effectively. This approach yields up to an 80% reduction in the 1-norm $λ$ and also leads to a 55% reduction in the number of orbitals. Our results highlight that improving the quality, not just the size, of the orbital basis is a viable strategy for extending QPE to include dynamical correlation, making progress toward scalable and chemically accurate quantum simulations with tractable resource requirements.