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

AI summary

Overview

• Research area: Quantum computing for electronic-structure simulation, specifically self-consistent Kohn–Sham density functional theory (KS-DFT).
• Technical level: Advanced.
• Scope: The paper introduces a quantum-circuit eigensolver for KS-DFT that combines plane waves with orthonormalized atomic orbitals and uses a compact subspace to reduce the qubit register needed for repeated calculations.

What This Paper Is About

KS-DFT repeatedly solves for electronic orbitals as the electron density is updated. Representing all-electron orbitals with plane waves alone can require a very large basis, while solving every step in that full basis is costly for quantum hardware. This paper develops a method that combines plane waves and atomic orbitals, then performs most self-consistent steps in a smaller basis while checking when the larger basis must be revisited.

Key Contributions

1. A plane-wave-orthonormalized-orbital (PWOO) basis: It combines a retained set of plane waves with atomic orbitals orthogonalized against those plane waves and each other. This gives a systematic plane-wave component while representing near-nucleus orbital features with localized functions.
2. A first-quantization KS-DFT eigensolver: It encodes one orbital at a time on a qubit register and finds multiple occupied orbitals using a variational circuit, classical angle optimization, deflation, and orthonormalization.
3. A dual-basis strategy: After initial full-basis solves, it solves most SCF steps in a contracted subspace. It checks residuals and updates that subspace when needed, keeping the contracted register size fixed.
4. Tests spanning simulation and hardware readout: The authors test pure plane-wave and PWOO calculations, emulate shot noise for PWOO calculations, and perform a trapped-ion readout of converged states.

Main Findings

• Agreement with classical calculations: For an H₈ chain with the SVWN functional, the state-vector-emulated calculation with 32,768 plane waves agrees with the classical plane-wave result within 0.02 mHa. The PWOO calculations agree with their classical counterparts within 0.44 μHa in the reported table.
• PWOO can use a much smaller basis than pure plane waves: The paper reports that PWOO reaches energies below the STO-3G reference with 16 basis functions, while pure plane waves require more than 4,096 functions to do so. The largest pure plane-wave calculation, with 32,768 functions, is still far from energy convergence.
• Contraction reduces register demands: The quantum calculations use a contracted space with K = 8 functions, which fits on a three-qubit register. The contraction does not eliminate the need for full-basis work: the method uses full-basis bootstrap steps and returns to the full basis for subspace updates.
• The contracted Hamiltonian is cheaper to measure: At the tested density, the full plane-wave Hamiltonian’s Pauli L1 norm grows from 6.97 Ha at n_q = 6 to 233 Ha at n_q = 12, while the corresponding contracted Hamiltonian norms remain around one Hartree. This indicates a much lower per-evaluation shot bound for the contracted problem.
• Shot-based SCF can approach the noise-free solution, but convergence is delicate: The reported STO-3G PWOO shot-SCF runs converge near the noise-free fixed point when using a reduced mixing parameter. The paper attributes persistent fluctuations partly to incomplete convergence of the eigensolver, not just measurement noise.
• The hardware test is a readout, not a hardware SCF calculation: On IonQ Forte-1, three-qubit-register readout reproduces the noise-free PWOO–AO energy difference within one standard deviation. The paper does not report running the full self-consistent loop on the device.
• The work does not claim a speed advantage over classical DFT: The reported timings compare algorithmic variants in the study and are not presented as a quantum-computing advantage.

Methodology in Plain English

The researchers tested the approach on an H₈ chain using the SVWN local-density functional. They first construct an orthonormal basis. In the pure plane-wave version, the basis consists of plane waves; in PWOO, atomic orbitals supplement the plane waves and describe parts of the orbitals that would otherwise require many high-frequency plane waves.

A quantum register represents one KS orbital at a time. A parameterized circuit prepares candidate orbitals, and a classical optimizer adjusts its angles to lower the orbital energy. The method finds occupied orbitals one after another, discourages overlap with already-found orbitals, and orthonormalizes them. Their density is then used to assemble the next KS Hamiltonian, continuing the SCF cycle.

To reduce the register needed for repeated solves, the method first finds a small set of low-energy states in the full basis. It then solves most later steps in the subspace spanned by those states. A residual test checks whether the subspace is still adequate; if not, the method updates it using another full-basis solve.

The study compares these calculations with classical references, simulates shot noise in some PWOO calculations, and tests readout of converged states on a trapped-ion device. The demonstrated molecular calculations use the Γ point; although the formulation is described for periodic systems in principle, the paper says its implementation is restricted to Γ.

Why This Matters

This work explores how to make self-consistent electronic-structure calculations more compatible with quantum registers: the orbital basis can be represented directly, and most repeated eigensolves can use a smaller subspace. Its combination of an all-electron, full-potential formulation with a mixed plane-wave and atomic-orbital basis is a distinctive step toward practical quantum-computing DFT workflows. The calculations remain demonstrations rather than evidence of a quantum advantage.

Potential research and application areas include:

• Molecular electronic structure: Studying molecular orbitals and energies, as demonstrated for the H₈ chain.
• Materials science: Developing quantum-computing methods for electronic-structure problems in materials; periodic use is proposed in principle, but not demonstrated here.
• Condensed-matter physics: Exploring self-consistent electronic calculations relevant to condensed-matter systems, an area where conventional DFT is widely used.
• Methods beyond semilocal DFT: Investigating how the framework might connect to hybrid functionals, DFT+U, or DFT with dynamical mean-field theory; these are discussed as broader DFT contexts, not demonstrated calculations in this paper.

Industry relevance: The work is relevant to quantum-computing hardware and software teams developing algorithms for chemistry and materials simulation. Its hardware result is limited to readout of converged states, and the authors explicitly do not claim a wall-time advantage over classical DFT.

Future Directions

1. Test larger and more varied systems: The paper demonstrates the method on one isolated H₈ chain; broader tests would show how the basis and eigensolver behave for other molecules or periodic systems.
2. Improve shot-based convergence: The observed fluctuations suggest a need to better separate measurement noise from eigensolver error and to improve convergence under noisy density updates.
3. Develop the hardware workflow: A natural next step is to run more of the self-consistent loop on quantum hardware rather than limiting the reported device test to readout.
4. Assess basis and grid convergence: The reported PWOO energies depend on the unit-cell grid and are not yet converged with respect to that grid in the reported correlation-consistent basis calculations.

Target Audience

This paper will be most useful to researchers in quantum algorithms, quantum chemistry, and computational materials science who understand DFT and want to see how orbital-based quantum eigensolvers might be embedded in a self-consistent workflow. It may also interest quantum hardware and software developers evaluating register size, measurement cost, and the gap between emulated calculations and device experiments.

Authors’ 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.