The quest for chemical accuracy in predicting molecular properties and reactivity remains a central challenge in computational chemistry. While Density Functional Theory (DFT) offers a favorable cost-to-accuracy ratio, its inherent approximations limit its ability to consistently achieve sub-kJ/mol accuracy, particularly for complex systems. Conversely, high-level wave function theory (WFT) methods, often termed "post-Hartree-Fock" (post-HF) methods, deliver unparalleled accuracy but are prohibitively expensive for systems beyond a few dozen atoms. This paper explores the synergistic potential of hybrid DFT/post-HF schemes in the context of emerging exascale computing architectures. We delineate how combining DFT's efficiency for electron correlation and the rigorous treatment of dynamic and static correlation by post-HF methods can overcome individual limitations. Exascale computing, with its unprecedented power and massive parallelism, provides a unique opportunity to accelerate these sophisticated hybrid methods. We discuss algorithmic challenges, software engineering requirements, and novel parallelization strategies to harness exascale capabilities effectively. By detailing specific hybrid schemes, such as double-hybrid functionals and explicitly correlated methods, we demonstrate their promise for pushing the boundaries of chemical accuracy for previously inaccessible systems. The ultimate goal is to enable predictive simulations that drive discovery in materials science, drug design, and catalysis, making chemical accuracy a routine harvest rather than an exceptional yield.
The quest for chemical accuracy in predicting molecular properties and reactivity remains a central challenge in computational chemistry. While Density Functional Theory (DFT) offers a favorable cost-to-accuracy ratio, its inherent approximations limit its ability to consistently achieve sub-kJ/mol accuracy, particularly for complex systems. Conversely, high-level wave function theory (WFT) methods, often termed "post-Hartree-Fock" (post-HF) methods, deliver unparalleled accuracy but are prohibitively expensive for systems beyond a few dozen atoms. This paper explores the synergistic potential of hybrid DFT/post-HF schemes in the context of emerging exascale computing architectures. We delineate how combining DFT's efficiency for electron correlation and the rigorous treatment of dynamic and static correlation by post-HF methods can overcome individual limitations. Exascale computing, with its unprecedented power and massive parallelism, provides a unique opportunity to accelerate these sophisticated hybrid methods. We discuss algorithmic challenges, software engineering requirements, and novel parallelization strategies to harness exascale capabilities effectively. By detailing specific hybrid schemes, such as double-hybrid functionals and explicitly correlated methods, we demonstrate their promise for pushing the boundaries of chemical accuracy for previously inaccessible systems. The ultimate goal is to enable predictive simulations that drive discovery in materials science, drug design, and catalysis, making chemical accuracy a routine harvest rather than an exceptional yield.
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 chemistry plays a crucial role in advancing fundamental understanding of chemical processes and driving innovation in energy, medicine, and materials science. Recent progress in quantum computing has opened new possibilities for molecular simulations that are beyond the practical limits of classical approaches. Academic studies have demonstrated hybrid Density Functional Theory (DFT) and Variational Quantum Eigensolver (VQE) benchmarks on small transition-metal systems, such as iron porphyrin and heme analogues, using current noisy intermediate-scale quantum (NISQ) hardware. Although full-scale simulations of complex biological systems such as the complete Cytochrome P450 (CYP450) active site remain a long-term goal due to limitations in qubit numbers, coherence times, and error correction, the underlying methodologies are now well established. At present, quantum chemistry applications in drug discovery remain largely experimental, and widespread practical medical impact is expected to require further technological advances over the next 5–10 years. Nevertheless, quantum software and algorithm development for chemical applications is progressing rapidly. This perspective summarizes recent advances in quantum computing algorithms, hardware, and software relevant to chemistry, and critically discusses the remaining challenges and opportunities for applying quantum computing to chemical problems, with particular emphasis on drug discovery and development.
Hybrid artificial intelligence–quantum (AI–Q) models arise at the intersection of two rapidly advancing computational paradigms, promising unprecedented capabilities for real-time molecular structure prediction. In classical computational chemistry, the accurate determination of ground-state energies and optimized geometries relies on ab initio methods such as density functional theory (DFT) and coupled-cluster approaches, which, while highly precise, exhibit steep polynomial or exponential scaling with system size. Conversely, deep neural networks (DNNs) have demonstrated remarkable success in learning potential energy surfaces (PES) from large datasets, achieving rapid inference but sometimes lacking the ultimate chemical accuracy required for predictive applications in drug discovery or materials design. On the other hand, variational quantum eigensolvers (VQEs) leverage the intrinsic parallelism of quantum hardware to approximate electronic structure problems, yet they remain constrained by qubit coherence times, gate fidelities, and circuit depth limitations inherent to noisy intermediate-scale quantum (NISQ) devices.
Accurate prediction of molecular structures and properties is vital for chemistry; materials science, and drug discovery, yet classical electronic-structure methods often fail for strongly correlated systems, large basis sets, and complex potential-energy landscapes. Quantum technology encompassing quantum computing, quantum machine learning, and hybrid quantum classical strategies offers a fundamentally new paradigm by encoding many-electron wave functions directly on qubits and exploiting superposition and entanglement to explore exponentially large Hilbert spaces. This review synthesizes recent algorithmic and hardware advances relevant to molecular modelling, including the Variational Quantum Eigensolver (VQE), Quantum Phase Estimation (QPE), quantum unitary coupled-cluster (q-UCC) families, equation-of-motion and subspace methods for excited states, tensor-network hybrids, and quantum kernel and variational QML approaches. We examine noise-aware hybrid workflows, error-mitigation techniques, symmetry-preserving ansätze, and operator-factorization methods that reduce measurement and gate overhead. Representative applications are discussed: ground- and excited-state energy prediction, potential-energy surface mapping, geometry and transition-state optimization, spectroscopic property estimation (IR, UV–Vis, NMR, EPR), and reaction-dynamics scenarios where non-adiabatic effects and conical intersections dominate. Resource estimates and scaling analyses clarify current NISQ limitations qubit counts, circuit depth, shot complexity and delineate the roadmap to fault-tolerant QPE for chemical accuracy. We compare quantum approaches with classical baselines (DFT, CCSD (T), multireference methods), identifying domains where quantum methods already show promise (strong correlation, multi-reference dissociation, spin-state ordering) and where classical methods remain competitive. Finally, we highlight near-term industrial opportunities in drug design, catalysis, CO2 capture, and energy materials, and outline critical research directions: algorithmic reductions in measurement/precision cost, hardware improvements in fidelity and connectivity, scalable ansatz design, and integrated software stacks for reproducible hybrid simulations. Together, these developments indicate that while practical, large-scale quantum advantage for general chemistry remains future work, quantum technologies are rapidly maturing into powerful tools for targeted molecular problems that are intractable with existing classical techniques.
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.
Electronic structure calculations have the potential to predict key matter transformations for applications of strategic technological importance, from drug discovery to material science and catalysis. However, a predictive physicochemical characterization of these processes often requires accurate quantum chemical modeling of complex molecular systems with hundreds to thousands of atoms. Due to the computationally demanding nature of electronic structure calculations and the complexity of modern high-performance computing hardware, quantum chemistry software has historically failed to operate at such large molecular scales with accuracy and speed that are useful in practice. In this paper, novel algorithms and software are presented that enable extreme-scale quantum chemistry capabilities with particular emphasis on exascale calculations. This includes the development and application of the multi-Graphics Processing Unit (GPU) library LibCChem 2.0 as part of the General Atomic and Molecular Electronic Structure System package and of the standalone Extreme-scale Electronic Structure System (EXESS), designed from the ground up for scaling on thousands of GPUs to perform high-performance accurate quantum chemistry calculations at unprecedented speed and molecular scales. Among various results, we report that the EXESS implementation enables Hartree-Fock/cc-pVDZ plus RI-MP2/cc-pVDZ/cc-pVDZ-RIFIT calculations on an ionic liquid system with 623 016 electrons and 146 592 atoms in less than 45 min using 27 600 GPUs on the Summit supercomputer with a 94.6% parallel efficiency.