Quantum computing has emerged as a promising platform for simulating strongly correlated systems in chemistry, for which the standard quantum chemistry methods are either qualitatively inaccurate or too expensive. However, due to the hardware limitations of the available noisy near-term quantum devices, their application is currently limited only to small chemical systems. One way for extending the range of applicability can be achieved within the quantum embedding approach. Herein, we employ the projection-based embedding method for combining the variational quantum eigensolver (VQE) algorithm, although not limited to, with density functional theory (DFT). The developed VQE-in-DFT method is then implemented efficiently on a real quantum device and employed for simulating the triple bond breaking process in butyronitrile. The results presented herein show that the developed method is a promising approach for simulating systems with a strongly correlated fragment on a quantum computer.
Quantum computing has emerged as a promising platform for simulating strongly correlated systems in chemistry, for which the standard quantum chemistry methods are either qualitatively inaccurate or too expensive. However, due to the hardware limitations of the available noisy near-term quantum devices, their application is currently limited only to small chemical systems. One way for extending the range of applicability can be achieved within the quantum embedding approach. Herein, we employ the projection-based embedding method for combining the variational quantum eigensolver (VQE) algorithm, although not limited to, with density functional theory (DFT). The developed VQE-in-DFT method is then implemented efficiently on a real quantum device and employed for simulating the triple bond breaking process in butyronitrile. The results presented herein show that the developed method is a promising approach for simulating systems with a strongly correlated fragment on a quantum computer.
Simulating strongly correlated molecular systems on near-term quantum hardware remains challenging because current hardware offers limited quantum volume and moderate-fidelity qubits. One way to address this challenge is through bootstrap embedding (BE). Bootstrap embedding partitions molecules into smaller fragments that are then embedded in the ``bath'' of other fragments iteratively. Bootstrap embedding is appealing for quantum simulation because fragmentation reduces the qubit requirements of each fragment problem. In this work, we develop a quantum bootstrap embedding (QBE) workflow that uses variational quantum eigensolver (VQE) fragment solvers and examines the algorithmic choices that determine the success of the overall VQE-QBE procedure. To improve efficiency, we introduce FastAdaptVQE, a sparse-matrix-accelerated form of the adaptive variational quantum eigensolver (ADAPT-VQE) that replaces symbolic commutator evaluation with direct statevector calculations, and MatrixFreeAdaptVQE, a matrix-free extension that removes the sparse-matrix memory bottleneck that arises for larger fragments. We also explore modifying the ADAPT-VQE operator-selection step by replacing the purely greedy choice with a lookahead strategy. Benchmarks on H4, F2, and selected geometries of H6 reach chemical accuracy, within 1 kcal/mol of BE results obtained with a full configuration interaction (FCI) solver. For n-butane, density-matching residuals converged, but energy errors remained outside chemical accuracy. These results show that combining QBE with VQE can accurately compute molecular energies, while density convergence alone does not ensure energy accuracy. This work lays the foundation for extending these energy calculations to larger molecular systems and quantum materials on near-term quantum hardware.
Quantum computers open up new avenues for modelling the physical properties of materials and molecules. Density Functional Theory (DFT) is the gold standard classical algorithm for predicting these properties, but relies on approximations of the unknown universal functional, limiting its general applicability for many fundamental and technologically relevant systems. In this work we develop a hybrid quantum/classical algorithm called quantum enhanced DFT (QEDFT) that systematically constructs quantum approximations of the universal functional using data obtained from a quantum computer. We benchmark the QEDFT algorithm on the Fermi-Hubbard model, both numerically and on data from experiments on real quantum hardware. We find that QEDFT surpasses the quality of groundstate results obtained from Hartree-Fock DFT, as well as from direct application of conventional quantum algorithms such as VQE. Furthermore, we demonstrate that QEDFT works even when only noisy, low-depth quantum computation is available, by benchmarking the algorithm on data obtained from Google's quantum computer. We further show how QEDFT also captures quintessential properties of strongly correlated Mott physics for large Fermi-Hubbard systems using functionals generated on much smaller system sizes. Our results indicate that QEDFT can be applied to realistic materials and molecular systems, and has the potential to outperform the direct application of either DFT or VQE alone, without the requirement of large scale or fully fault-tolerant quantum computers.
Variational quantum eigensolver (VQE) is a hybrid quantum-classical algorithm designed for noisy intermediate-scale quantum (NISQ) computers. It is promising for quantum chemical calculations (QCC) because it can calculate the ground-state energy of a target molecule. Although VQE has a potential to achieve a higher accuracy than classical approximation methods in QCC, it is challenging to achieve it on current NISQ computers due to the significant impact of noises. Density matrix embedding theory (DMET) is a well-known technique to divide a molecule into multiple fragments, which is available to mitigate the noise impact on VQE. However, our preliminary evaluation shows that the naive combination of DMET and VQE does not outperform a gold standard classical method. In this work, we present three approaches to mitigate the noise impact for the DMET+VQE combination. (1) The size of quantum circuits used by VQE is decreased by reducing the number of bath orbitals which represent interactions between multiple fragments in DMET. (2) Reduced density matrices (RDMs), which are used to calculate a molecular energy in DMET, are calculated accurately based on expectation values obtained by executing quantum circuits using a noise-less quantum computer simulator. (3) The parameters of a quantum circuit optimized by VQE are refined with mathematical post-processing. The evaluation using a noisy quantum computer simulator shows that our approaches significantly improve the accuracy of the DMET+VQE combination. Moreover, we demonstrate that on a real NISQ device, the DMET+VQE combination applying our three approaches achieves a higher accuracy than the gold standard classical method.
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.
High Resolution Image Download MS PowerPoint Slide This study explores the application of quantum computing to metal cluster analysis through the development and implementation of a quantum-DFT embedding workflow. Classical computational methods, while transformative, often face limitations in achieving chemical accuracy and computational efficiency, particularly for nanoscale systems. To address these challenges, we integrate the Variational Quantum Eigensolver (VQE) with density functional theory (DFT), leveraging the capabilities of quantum computing aiming to improve the modeling of electronic structures. Aluminum and gold clusters were used as model systems to test the established workflow. The workflow successfully determined electronic properties for aluminum clusters up to Al 7 – . Although gold clusters were used as a test case to investigate the potential reduction of nitric oxide (NO), memory limitations, the lack of relativistic corrections, and the inability to handle open-shell systems presented challenges that underscore the need for advancements in quantum hardware and algorithms. This proof-of-concept study demonstrates the potential of quantum DFT embedding to advance materials discovery, including applications in catalysis and nanomaterial design, while providing insights into the current limitations of near-term quantum devices.
Quantum mechanics has introduced a new theoretical framework for the study of molecules, enabling the prediction of properties and dynamics through the solution of the Schrödinger equation applied to these systems. However, solving this equation is computationally expensive, which has led to the development of various mathematical frameworks and computational methods designed to balance the available resources with the desired level of accuracy. In particular, quantum computers have emerged as a promising technology with the potential to address these problems more efficiently in the coming decades, whether through reductions in memory, time, and energy consumption, i. e., reductions in computational complexity or by enhancing precision. This research field is known as Quantum Simulation. Given the current technological limitations of quantum computers, Variational Quantum Algorithms (VQAs), especially the Variational Quantum Eigensolver (VQE), stand out as a feasible approach to demonstrating advantages over classical methods in the near term. This feasibility arises from their lower demand for quantum gates and the reduced depth of the circuits required for their implementation. In this work, we present the application of quantum mechanics to the study of molecules, provide an introduction to the fundamentals of quantum computing, and explore the integration of these fields by employing the VQE in molecular simulations. Finally, we discuss the spatial and temporal complexity associated with the algorithm, highlighting its implications and challenges.