Abstract:Density-functional theory (DFT) has revolutionized computer simulations in chemistry and material science. A faithful implementation of the theory requires self-consistent calculations. However, this effort involves repeatedly diagonalizing the Hamiltonian, for which a classical algorithm typically requires a computational complexity that scales cubically with respect to the number of electrons. This limits DFT's applicability to large-scale problems with complex chemical environments and microstructures. This article presents a quantum algorithm that has a linear scaling with respect to the number of atoms, which is much smaller than the number of electrons. Our algorithm leverages the quantum singular value transformation (QSVT) to generate a quantum circuit to encode the density-matrix, and an estimation method for computing the output electron density. In addition, we present a randomized block coordinate fixed-point method to accelerate the self-consistent field calculations by reducing the number of components of the electron density that needs to be estimated. The proposed framework is accompanied by a rigorous error analysis that quantifies the function approximation error, the statistical fluctuation, and the iteration complexity. In particular, the analysis of our self-consistent iterations takes into account the measurement noise from the quantum circuit. These advancements offer a promising avenue for tackling large-scale DFT problems, enabling simulations of complex systems that were previously computationally infeasible.

A New History-Driven Algorithm to Calculate the Density of States

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Efficient sampling using Macrocanonical Monte Carlo and density of states mapping

Calculating the Density of States for Large-Scale Molecular Systems

Self-consistent Density Functional Calculations without Diagonalization

The Accuracy and Performance Analysis of the 1/t Wang-Landau Algorithm in the Joint Density of States Estimation

Dynamically Optimized Wang-Landau Sampling with Adaptive Trial Moves and Modification Factors

A new algorithm for Monte Carlo simulation of Ising spin systems

Calculating the many-body density of states on a digital quantum computer

Dynamical computation of the density of states and Bayes factors using nonequilibrium importance sampling

Wide Sampling and Efficient Updating Monte Carlo Algorithms for Dimer Models

Density of states of a two-dimensional XY model from the Wang-Landau algorithm.

Surveying an Energy Landscape

A statistical approach for simulating the density solution of a McKean-Vlasov equation

A portable and flexible implementation of the Wang--Landau algorithm in order to determine the Density of States

A Recipe for Charge Density Prediction

Monte Carlo simulation method for Laughlin-like states in a disk geometry

Reversible algorithm of simulating multivariate densities with multi-hump

Implementation of the Density-functional Theory on Quantum Computers with Linear Scaling with respect to the Number of Atoms

Optimal Ricci curvature Markov chain Monte Carlo methods on finite states

Density of states of a 2D system of soft--sphere fermions by path integral Monte Carlo simulations