Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2609.33934

Quantum-Computing Self-Consistent Kohn-Sham DFT: Plane-Wave-Orthonormalized Orbitals with a Dual-Basis Quantum Eigensolver

Lazaro Calderin

quant-phphysics.chem-phphysics.comp-phSubmitted Sep 27, 2026• Announced Sep 29, 2026• 0 citations
arXivHTMLPDF

Abstract

We present a first-quantization, all-electron, full-potential Kohn-Sham density functional theory method for quantum computers, based on a classically optimized quantum-circuit eigensolver for the self-consistent field (SCF) eigenproblems in any orthonormal basis. The plane-wave-orthonormalized-orbital (PWOO) basis combines plane waves with atomic orbitals (AOs) orthonormalized against them and one another. The eigensolver has a dual-basis form: after full-basis bootstrap steps, it solves each SCF step in a contracted basis of $K$ functions ($\lceil\log_2 K\rceil$ qubits), updated in the full basis when the residual outside it exceeds a threshold. For an H$_8$ chain (SVWN, $K=8$), the state-vector-emulated SCF agrees with classical calculations within 0.02 mHa for 32768 plane waves and PWOO bases up to cc-pVTZ; contraction changes these PWOO energies by at most 5.7 $μ$Ha. With emulated shot noise, the STO-3G PWOO SCF converges near the noise-free fixed point. On IonQ Forte-1, three-qubit-register readout reproduces the noise-free PWOO-AO energy difference within one standard deviation.

Related papers

Self-Consistent Field Approach for the Variational Quantum Eigensolver: Orbital Optimization Goes Adaptive

Aaron Fitzpatrick, Anton Nykänen, N. Walter Talarico +4 more

We present a self-consistent field (SCF) approach within the adaptive derivative-assembled problem-tailored ansatz variational quantum eigensolver (ADAPT-VQE) framework for efficient quantum simulations of chemical systems on near-term quantum computers. To this end, our ADAPT-VQE-SCF approach combines the idea of generating an ansatz with a small number of parameters, resulting in shallow-depth quantum circuits with a direct minimization of an energy expression that is correct to second order with respect to changes in the molecular orbital basis. Our numerical analysis, including calculations for the transition-metal complex ferrocene [Fe (C 5 H 5 ) 2 ], indicates that convergence in the self-consistent orbital optimization loop can be reached without a considerable increase in the number of two-qubit gates in the quantum circuit by comparison to a VQE optimization in the initial molecular orbital basis. Moreover, the orbital optimization can be carried out simultaneously within each iteration of the ADAPT-VQE cycle. ADAPT-VQE-SCF thus allows us to implement a routine analogous to the complete active space SCF, a cornerstone of state-of-the-art computational chemistry, in a hardware-efficient manner on near-term quantum computers. Hence, ADAPT-VQE-SCF paves the way toward a paradigm shift for quantitative quantum-chemistry simulations on quantum computers by requiring fewer qubits and opening up for the use of large and flexible atomic orbital basis sets in contrast to earlier methods that are predominantly based on the idea of full active spaces with minimal basis sets.

arXiv:2212.11405•The Journal of Physical Chemistry A•Mar 28, 2024•22 citations

Quantum algorithms for density functional theory with minimal readout

Yuansheng Zhao, Hirofumi Nishi, Taichi Kosugi +5 more

While quantum computers have shown significant promise for electronic structure calculations, their potential to accelerate density functional theory (DFT) calculations remains unclear. In this work, we present a qubit-efficient encoding scheme for wavefunctions in Kohn--Sham (KS) DFT, together with a quantum algorithm that computes all occupied orbitals simultaneously. We further show that our algorithm is particularly well suited to the Harris functional, enabling the total energy to be evaluated with a potential exponential speedup over classical approaches by entirely avoiding the costly readout of the electronic density. In addition, we propose a second method for achieving self-consistent DFT calculations using multiple copies of the wavefunction, which likewise circumvents density readout. The applicability of our algorithms is demonstrated through several numerical examples, and their efficiency is compared with that of existing approaches.

quant-pharXiv:2605.29774•May 29, 2026•0 citations

A Tensor-product approach for large scale electronic structure calculations using Kohn–Sham density functional theory

Phani Motamarri, Vikram Gavini, Thomas Blesgen

Quantum-mechanical calculations based on Kohn–Sham density functional theory (DFT) played a significant role in accurately predicting various aspects of materials behavior over the past decade. The Kohn–Sham approach to DFT reduces the many-body Schrodinger (eigen value) problem of interacting electrons into an equivalent problem of noninteracting electrons in an effective mean field that is governed by electron-density. Despite the reduced computational complexity of Kohn–Sham DFT, large-scale DFT calculations are still computationally very demanding with the resulting computational complexity scaling cubically with number of atoms in a given materials system. Numerical algorithms with reduced scaling behavior which are robust, computationally efficient and scalable on parallel computing architectures are always desirable to enable simulations at larger scales and on more complex systems. Following this line of thought, this study explores the use of tensor structured methods for ab-initio numerical solution of Kohn–Sham equations arising in DFT calculations. Earlier studies on tensor-structured methods have been quite successful in the accurate calculation of Hartree and the nonlocal exchange operators arising in the Hartree–Fock equations. A recent investigation of low-rank Tucker-type decomposition of the electron-density of large aluminum clusters (obtained from the finite-element discretization of orbital free DFT) shows the exponential decay of approximation error with respect to Tucker rank (number of tensor-basis functions in Tucker type representation). The results also indicate a smaller Tucker rank for the accurate representation of the electron density and is only weakly dependent on the system sizes studied. The promising success of tensor-structured techniques in resolving the electronic structure of material systems has enabled us to take a step further. In this study, we propose a systematic way of computing a globally adapted Tucker-type basis for solving the Kohn–Sham DFT problem by using a separable approximation of the Kohn–Sham Hamiltonian. Further, the resulting Kohn–Sham eigenvalue problem is projected into the aforementioned Tucker basis and is solved for ground-state energy using a self-consistent field iteration. The rank of the resulting Tucker representation and the computational complexity of the calculation are examined on representative benchmark examples involving metallic and insulating systems.

•Jan 1, 2014•0 citations

An algorithm for atom-centered lossy compression of the atomic orbital basis in density functional theory calculations

Anthony O. Lara, Justin J. Talbot, Zhe Wang +1 more

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-zeta (QZ) size are required to adequately converge Kohn-Sham density functional theory (DFT) calculations towards the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. In this work, we develop and test a natural atomic orbital (NAO) scheme in which the NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. The NAO representation enables one-center compression of the AO basis in a manner that is optimal for a given threshold, by discarding NAOs with occupation numbers below that threshold. Extensive tests using the Hartree-Fock functional suggest that a threshold of $10^{-5}$ can yield a compression factor (ratio of AO to compressed NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold $10^{-7}$, while the compression factor is typically between 2 and 2.5.

physics.chem-phCross-listedarXiv:2512.00118•Dec 1, 2025

Quantum Solvers for Plane-Wave Hamiltonians: Abridging Virtual Spaces Through the Optimization of Pairwise Correlations

Eric J. Bylaska, Duo Song, Nicholas P. Bauman +3 more

For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.

arXiv:2009.00080•Frontiers in Chemistry•Mar 18, 2021•19 citations

Generalized Kohn-Sham system in one-matrix functional theory

Ryan Requist, Oleg Pankratov

A system of electrons in a local or nonlocal external potential can be described with one-matrix functional theory (1MFT), which is similar to density-functional theory (DFT) but takes the one-particle reduced density matrix (one-matrix) instead of the density as its basic variable. Within 1MFT, Gilbert derived [Phys. Rev. B 12, 2111 (1975)] effective single-particle equations analogous to the Kohn-Sham (KS) equations in DFT. The self-consistent solution of these 1MFT-KS equations reproduces not only the density of the original electron system but also its one-matrix. While in DFT it is usually possible to reproduce the density using KS orbitals with integer (0 or 1) occupancy, in 1MFT reproducing the one-matrix requires in general fractional occupancies. The variational principle implies that the KS eigenvalues of all fractionally occupied orbitals must collapse at self-consistency to a single level. We show that as a consequence of the degeneracy, the iteration of the KS equations is intrinsically divergent. Fortunately, the level-shifting method, commonly introduced in Hartree-Fock calculations, is always able to force convergence. We introduce an alternative derivation of the 1MFT-KS equations that allows control of the eigenvalue collapse by constraining the occupancies. As an explicit example, we apply the 1MFT-KS scheme to calculate the ground state one-matrix of an exactly solvable two-site Hubbard model.

arXiv:0711.2996•Physical Review B•Jun 25, 2008•37 citations
We present a first-quantization, all-electron, full-potential Kohn-Sham density functional theory method for quantum computers, based on a classically optimized quantum-circuit eigensolver for the self-consistent field (SCF) eigenproblems in any orthonormal basis. The plane-wave-orthonormalized-orbital (PWOO) basis combines plane waves with atomic orbitals (AOs) orthonormalized against them and one another. The eigensolver has a dual-basis form: after full-basis bootstrap steps, it solves each SCF step in a contracted basis of KKK functions (⌈log⁡2K⌉\lceil\log_2 K\rceil⌈log2​K⌉ qubits), updated in the full basis when the residual outside it exceeds a threshold. For an H8_88​ chain (SVWN, K=8K=8K=8), the state-vector-emulated SCF agrees with classical calculations within 0.02 mHa for 32768 plane waves and PWOO bases up to cc-pVTZ; contraction changes these PWOO energies by at most 5.7 μμμHa. With emulated shot noise, the STO-3G PWOO SCF converges near the noise-free fixed point. On IonQ Forte-1, three-qubit-register readout reproduces the noise-free PWOO-AO energy difference within one standard deviation.
We present a self-consistent field (SCF) approach within the adaptive derivative-assembled problem-tailored ansatz variational quantum eigensolver (ADAPT-VQE) framework for efficient quantum simulations of chemical systems on near-term quantum computers. To this end, our ADAPT-VQE-SCF approach combines the idea of generating an ansatz with a small number of parameters, resulting in shallow-depth quantum circuits with a direct minimization of an energy expression that is correct to second order with respect to changes in the molecular orbital basis. Our numerical analysis, including calculations for the transition-metal complex ferrocene [Fe (C 5 H 5 ) 2 ], indicates that convergence in the self-consistent orbital optimization loop can be reached without a considerable increase in the number of two-qubit gates in the quantum circuit by comparison to a VQE optimization in the initial molecular orbital basis. Moreover, the orbital optimization can be carried out simultaneously within each iteration of the ADAPT-VQE cycle. ADAPT-VQE-SCF thus allows us to implement a routine analogous to the complete active space SCF, a cornerstone of state-of-the-art computational chemistry, in a hardware-efficient manner on near-term quantum computers. Hence, ADAPT-VQE-SCF paves the way toward a paradigm shift for quantitative quantum-chemistry simulations on quantum computers by requiring fewer qubits and opening up for the use of large and flexible atomic orbital basis sets in contrast to earlier methods that are predominantly based on the idea of full active spaces with minimal basis sets.
While quantum computers have shown significant promise for electronic structure calculations, their potential to accelerate density functional theory (DFT) calculations remains unclear. In this work, we present a qubit-efficient encoding scheme for wavefunctions in Kohn--Sham (KS) DFT, together with a quantum algorithm that computes all occupied orbitals simultaneously. We further show that our algorithm is particularly well suited to the Harris functional, enabling the total energy to be evaluated with a potential exponential speedup over classical approaches by entirely avoiding the costly readout of the electronic density. In addition, we propose a second method for achieving self-consistent DFT calculations using multiple copies of the wavefunction, which likewise circumvents density readout. The applicability of our algorithms is demonstrated through several numerical examples, and their efficiency is compared with that of existing approaches.
Quantum-mechanical calculations based on Kohn–Sham density functional theory (DFT) played a significant role in accurately predicting various aspects of materials behavior over the past decade. The Kohn–Sham approach to DFT reduces the many-body Schrodinger (eigen value) problem of interacting electrons into an equivalent problem of noninteracting electrons in an effective mean field that is governed by electron-density. Despite the reduced computational complexity of Kohn–Sham DFT, large-scale DFT calculations are still computationally very demanding with the resulting computational complexity scaling cubically with number of atoms in a given materials system. Numerical algorithms with reduced scaling behavior which are robust, computationally efficient and scalable on parallel computing architectures are always desirable to enable simulations at larger scales and on more complex systems. Following this line of thought, this study explores the use of tensor structured methods for ab-initio numerical solution of Kohn–Sham equations arising in DFT calculations. Earlier studies on tensor-structured methods have been quite successful in the accurate calculation of Hartree and the nonlocal exchange operators arising in the Hartree–Fock equations. A recent investigation of low-rank Tucker-type decomposition of the electron-density of large aluminum clusters (obtained from the finite-element discretization of orbital free DFT) shows the exponential decay of approximation error with respect to Tucker rank (number of tensor-basis functions in Tucker type representation). The results also indicate a smaller Tucker rank for the accurate representation of the electron density and is only weakly dependent on the system sizes studied. The promising success of tensor-structured techniques in resolving the electronic structure of material systems has enabled us to take a step further. In this study, we propose a systematic way of computing a globally adapted Tucker-type basis for solving the Kohn–Sham DFT problem by using a separable approximation of the Kohn–Sham Hamiltonian. Further, the resulting Kohn–Sham eigenvalue problem is projected into the aforementioned Tucker basis and is solved for ground-state energy using a self-consistent field iteration. The rank of the resulting Tucker representation and the computational complexity of the calculation are examined on representative benchmark examples involving metallic and insulating systems.
Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-zeta (QZ) size are required to adequately converge Kohn-Sham density functional theory (DFT) calculations towards the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. In this work, we develop and test a natural atomic orbital (NAO) scheme in which the NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. The NAO representation enables one-center compression of the AO basis in a manner that is optimal for a given threshold, by discarding NAOs with occupation numbers below that threshold. Extensive tests using the Hartree-Fock functional suggest that a threshold of 10−510^{-5}10−5 can yield a compression factor (ratio of AO to compressed NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10−710^{-7}10−7, while the compression factor is typically between 2 and 2.5.
For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.
A system of electrons in a local or nonlocal external potential can be described with one-matrix functional theory (1MFT), which is similar to density-functional theory (DFT) but takes the one-particle reduced density matrix (one-matrix) instead of the density as its basic variable. Within 1MFT, Gilbert derived [Phys. Rev. B 12, 2111 (1975)] effective single-particle equations analogous to the Kohn-Sham (KS) equations in DFT. The self-consistent solution of these 1MFT-KS equations reproduces not only the density of the original electron system but also its one-matrix. While in DFT it is usually possible to reproduce the density using KS orbitals with integer (0 or 1) occupancy, in 1MFT reproducing the one-matrix requires in general fractional occupancies. The variational principle implies that the KS eigenvalues of all fractionally occupied orbitals must collapse at self-consistency to a single level. We show that as a consequence of the degeneracy, the iteration of the KS equations is intrinsically divergent. Fortunately, the level-shifting method, commonly introduced in Hartree-Fock calculations, is always able to force convergence. We introduce an alternative derivation of the 1MFT-KS equations that allows control of the eigenvalue collapse by constraining the occupancies. As an explicit example, we apply the 1MFT-KS scheme to calculate the ground state one-matrix of an exactly solvable two-site Hubbard model.