Quantum simulations of molecular electronic structure hold the potential for computational advantage but are inherently limited by the linear scaling of qubit requirements with respect to system size. Quantum Density Functional Theory (Q-DFT) addresses this constraint by reformulating the interacting many-body problem in terms of a non-interacting Kohn--Sham reference system. Here, we present and experimentally validate a Q-DFT algorithm that requires only $\mathcal{O}(\log_2(N))$ qubits to simulate a system of $N$ orbitals, thereby reducing resource overhead. More precisely, we solve the Kohn--Sham problem for the lithium hydride (LiH) molecule using the ensemble extension of the variational quantum eigensolver paradigm that minimizes the energy of the system in the subspace spanned by the ground and first excited states. To execute this task, we encode the problem with qubits on a quantum photonic platform, utilizing an architecture that combines a deterministic quantum dot single-photon source with a programmable integrated photonic circuit. By evaluating the molecular ground- and first excited-state properties, this work establishes a pathway for simulating larger molecular systems, demonstrating the viability of integrating qubit-efficient algorithms with advanced reconfigurable photonic hardware for near-term quantum chemistry applications.
AI Summary
The paper investigates whether qubit-efficient quantum density functional theory (Q-DFT) can be implemented on programmable photonic hardware to solve molecular Kohn–Sham problems. It applies an existing binary encoding requiring $\lceil\log_2 N\rceil$ qubits for $N$ basis orbitals to LiH in a reduced STO-3G basis containing four spatial orbitals, using two qubits. An equal-weight ensemble VQE optimizes the subspace of the two occupied Kohn–Sham orbitals; density-matrix tomography supplies the electronic density.
The implementation combines a quantum-dot single-photon source with an eight-mode programmable interferometer. Two dual-rail qubits undergo a four-parameter rotation ansatz containing one post-selected CNOT, optimized using COBYLA. State-vector simulations complete six self-consistent iterations, reaching an energy error below $10^{-6}$ Hartree and a density error around $10^{-3}$ relative to PySCF. However, finite-sampling simulations do not fully converge within six iterations even with $10^9$ shots, plateauing near $10^{-3}$ Hartree energy error and $10^{-1}$ density error.
Consequently, the experiment optimizes six Hamiltonians generated offline rather than closing the self-consistent loop on hardware. The measured photon indistinguishability is approximately $89\%$, and the reconstructed ensemble density matrix has $90.5\%$ fidelity with the ideal reference. Noise simulations incorporating distinguishability and counting statistics broadly reproduce the optimization behavior, but experimental densities do not quantitatively recover the ionic bond. The stated limitations include additional unmodeled control errors, optical losses, high sampling requirements, and probabilistic gates. The logarithmic qubit count does not establish a computational advantage, and accuracy remains dependent on the approximate exchange–correlation functional.
Little sign of AI-generated textOverclaimsIncremental
Central claim may not hold: The abstract presents experimental validation of Q-DFT, but the experiment only optimizes Kohn–Sham Hamiltonians generated offline and never feeds measured densities back into a self-consistent calculation.
The equal-weight ensemble VQE determines the two-dimensional occupied Kohn–Sham subspace rather than individual eigenstates, so it does not directly establish the abstract's claimed ground- and first-excited-state properties or molecular excitation properties.
The experimental results do not quantitatively reproduce the target electronic structure, with the paper explicitly stating that they are insufficient to accurately describe LiH's ionic bond.
Written by gpt-6.1-sol from the paper's full text. Automated: it can be wrong.
Quantum simulations of molecular electronic structure hold the potential for computational advantage but are inherently limited by the linear scaling of qubit requirements with respect to system size. Quantum Density Functional Theory (Q-DFT) addresses this constraint by reformulating the interacting many-body problem in terms of a non-interacting Kohn--Sham reference system. Here, we present and experimentally validate a Q-DFT algorithm that requires only O(log2(N)) qubits to simulate a system of N orbitals, thereby reducing resource overhead. More precisely, we solve the Kohn--Sham problem for the lithium hydride (LiH) molecule using the ensemble extension of the variational quantum eigensolver paradigm that minimizes the energy of the system in the subspace spanned by the ground and first excited states. To execute this task, we encode the problem with qubits on a quantum photonic platform, utilizing an architecture that combines a deterministic quantum dot single-photon source with a programmable integrated photonic circuit. By evaluating the molecular ground- and first excited-state properties, this work establishes a pathway for simulating larger molecular systems, demonstrating the viability of integrating qubit-efficient algorithms with advanced reconfigurable photonic hardware for near-term quantum chemistry applications.
The paper investigates whether qubit-efficient quantum density functional theory (Q-DFT) can be implemented on programmable photonic hardware to solve molecular Kohn–Sham problems. It applies an existing binary encoding requiring ⌈log2N⌉ qubits for N basis orbitals to LiH in a reduced STO-3G basis containing four spatial orbitals, using two qubits. An equal-weight ensemble VQE optimizes the subspace of the two occupied Kohn–Sham orbitals; density-matrix tomography supplies the electronic density.
The implementation combines a quantum-dot single-photon source with an eight-mode programmable interferometer. Two dual-rail qubits undergo a four-parameter rotation ansatz containing one post-selected CNOT, optimized using COBYLA. State-vector simulations complete six self-consistent iterations, reaching an energy error below 10−6 Hartree and a density error around 10−3 relative to PySCF. However, finite-sampling simulations do not fully converge within six iterations even with 109 shots, plateauing near 10−3 Hartree energy error and 10−1 density error.
Consequently, the experiment optimizes six Hamiltonians generated offline rather than closing the self-consistent loop on hardware. The measured photon indistinguishability is approximately 89%, and the reconstructed ensemble density matrix has 90.5% fidelity with the ideal reference. Noise simulations incorporating distinguishability and counting statistics broadly reproduce the optimization behavior, but experimental densities do not quantitatively recover the ionic bond. The stated limitations include additional unmodeled control errors, optical losses, high sampling requirements, and probabilistic gates. The logarithmic qubit count does not establish a computational advantage, and accuracy remains dependent on the approximate exchange–correlation functional.
Solving the electronic structure problem is a notorious challenge in quantum chemistry and material science. Variational quantum eigensolver (VQE) is a promising hybrid classical-quantum algorithm for finding the lowest-energy configuration of a molecular system. However, it typically requires many qubits and quantum gates with substantial quantum circuit depth to accurately represent the electronic wave function of complex structures. Here, we propose an alternative approach to solve the electronic structure problem using VQE with a single qudit. Our approach exploits a high-dimensional orbital angular momentum state of a heralded single photon and notably reduces the required quantum resources compared to conventional multi-qubit–based VQE. We experimentally demonstrate that our single-qudit–based VQE can efficiently estimate the ground state energy of hydrogen (H 2 ) and lithium hydride (LiH) molecular systems corresponding to two- and four-qubit systems, respectively. We believe that our scheme opens a pathway to perform a large-scale quantum simulation for solving more complex problems in quantum chemistry and material science.
Quantum chemistry models atoms and molecules using quantum mechanics, with the electronic structure problem at its core. This problem, which determines the energy and properties of electrons in a molecule, is exponentially hard to solve exactly on classical computers. To approximate solutions, effective models introduce layers of simplifications. However, these approximations can fail for strongly correlated systems, necessitating alternative approaches. Quantum computers, capable of processing quantum information natively, offer a promising path to simulate molecular systems beyond classical reach, with potential applications in catalysis, materials science, and photochemistry. Despite rapid advancements in quantum hardware, practical quantum computations remain challenging due to limited resources. This thesis optimizes quantum resource usage in computational chemistry, introducing novel embedding techniques and exploring applications beyond ground-state energies. Chapter 1 establishes the foundation by reviewing quantum chemistry methods, from Hartree-Fock and Density Functional Theory (DFT) to active space approaches. It also provides an introduction to quantum computing for chemistry, focusing on near-term algorithms. The following chapters present novel methods aimed at reducing quantum resource requirements and tackling problems classical methods struggle with. Chapter 2 introduces a technique to reduce the λ-norm of the electronic structure Hamiltonian, a key measure of quantum resource demand. By transforming molecular orbitals, this method enhances sampling efficiency in Noisy Intermediate-Scale Quantum (NISQ) algorithms and reduces runtime in Fault-Tolerant Quantum Computing (FTQC). Localized orbitals outperform canonical ones, and an orbital-optimization scheme is proposed to minimize the λ-norm. Benchmarks on hydrogen chains, FeMoco, and Ruthenium complexes demonstrate its effectiveness. Chapter 3 presents the FragPT2 framework for simulating large systems with multiple interacting fragments. This method embeds each fragment in the mean field of others and applies multi-reference perturbation theory to capture inter-fragment correlations such as dispersion and charge transfer efficiently. The approach reduces computational overhead and enables scalable simulations, as shown in applications to \ce{N2} dimers, aromatic dimers, and butadiene. Chapter 4 focuses on photochemical systems, where multiple electronic states must be treated equally. The state-averaged orbital-optimized variational quantum eigensolver (SAOOVQE) is extended with analytical energy gradients and non-adiabatic couplings, enabling geometry optimizations and simulations of non-adiabatic dynamics on quantum hardware. This is demonstrated by locating the conical intersection in formaldimine, highlighting SAOOVQE’s potential for studying photoisomerization in biomolecules. Building on this, Chapter 5 develops a quantum algorithm to detect conical intersections, where the Born-Oppenheimer approximation breaks down. These intersections are characterized by a discrete Berry phase (0 or π), representing their topological nature. The algorithm extracts this phase as a single bit of information, ensuring resilience to noise. Analytical proof and numerical validation using a toy molecular model illustrate its effectiveness for studying light-perception biochemistry and related photochemical processes. Finally, Chapter 6 explores machine learning density functionals trained on noisy quantum data. While DFT is computationally efficient, its accuracy depends on approximated exchange-correlation functionals. This work proposes a neural network-based approach to learn the exact functional from quantum data, showing that training on noisy outputs from NISQ devices can mitigate sampling errors. Applications to the Fermi-Hubbard model and other systems demonstrate its potential to generalize and provide accurate results under noise. This thesis advances quantum computational chemistry by optimizing quantum resource usage, developing new embedding strategies, and expanding the scope of quantum simulations beyond ground-state properties. These contributions pave the way for practical quantum algorithms that could revolutionize molecular simulations in the near future.
Scientific simulations demand methods combining scalability with predictive accuracy. Density Functional Theory (DFT) on High-Performance Computing (HPC) enables large-scale electronic-structure simulations but is limited by approximations affecting strongly correlated systems and band-gap predictions. Quantum computing offers a pathway to address this, though current Noisy Intermediate-Scale Quantum (NISQ) hardware remains constrained by qubit resources, noise, and execution cost.
This work presents a hybrid DFT-Quantum Embedding (QDFT) framework integrating classical HPC-based DFT with a quantum electronic-structure solver. Large systems are partitioned to isolate a chemically relevant active space, treated via the Variational Quantum Eigensolver (VQE), while the remaining degrees of freedom are described by DFT. The framework incorporates active-space selection, embedded Hamiltonian construction, symmetry preservation, operator mapping, self-consistent density updating, and modular classical-quantum coupling.
We focus on noiseless quantum simulation to systematically evaluate accuracy, convergence, active-space dependence, computational cost, and HPC scalability without hardware noise. Detailed profiling identifies computational bottlenecks and highlights limitations of CPU-based quantum simulation. A QPU runtime-estimation methodology is additionally developed to assess execution requirements on actual quantum hardware.
Results demonstrate quantum embedding's potential to improve selected electronic-structure properties while retaining classical HPC's scalability. Noisy quantum simulation and QPU execution remain key future directions, providing a pathway toward practical, scalable HPC-quantum hybrid simulations as hardware matures.
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.
In the near future, material and drug design may be aided by quantum computer assisted simulations. These have the potential to target chemical systems intractable by the most powerful classical computers. However, the resources offered by contemporary quantum computers are still limited, restricting the simulations to very simple molecules. In order to rapidly scale up to more interesting molecular systems, we propose the embedding of the quantum electronic structure calculation into a classically computed environment obtained at the Hartree-Fock (HF) or density functional theory (DFT) level of theory. This result is achieved by constructing an effective Hamiltonian that incorporates a mean field potential describing the action of the inactive electrons on a selected Active Space (AS). The ground state of the AS Hamiltonian is then determined by means of the variational quantum eigensolver algorithm. We show that with the proposed HF and DFT embedding schemes, we can obtain significant energy corrections to the reference HF and DFT calculations for a number of simple molecules in their strongly correlated limit (the dissociation regime) as well as for systems of the size of the oxirane molecule.
Quantum computers have the potential to deliver speed-ups for solving certain important problems that are intractable for classical counterparts, making them a promising avenue for advancing modern computation. However, many quantum algorithms require deep quantum circuits, which are challenging to implement on current noisy devices. To address this limitation, variational quantum algorithms (VQAs) have been actively developed, enabling practical quantum computing in the noisy intermediate-scale quantum (NISQ) era. Among them, the variational quantum eigensolver (VQE) stands out as a leading approach for solving problems in quantum chemistry, many-body physics, and even integer factorization. The VQE algorithm can be implemented on various quantum hardware platforms, including photonic systems, quantum dots, trapped ions, neutral atoms, and superconducting circuits. In particular, photonic platforms offer several advantages: they operate at room temperature, exhibit low decoherence, and support multiple degrees of freedom, making them suitable for scalable, high-dimensional quantum computation. Here we present methodologies for realizing VQE on photonic systems, highlighting their potential for practical quantum computing. We first provide a theoretical overview of the VQE framework, focusing on the procedure for variationally estimating ground state energies. We then explore how photonic systems can implement these processes, showing that a wide variety of problems can be addressed using either multiple qubit states or a single qudit state.
Bruno Senjean