Quantum chemistry applications on quantum computers currently rely heavily on the variational quantum eigensolver (VQE) algorithm. This hybrid quantum-classical algorithm aims at finding ground-state solutions of molecular systems based on the variational principle. VQE calculations can be systematically implemented for perturbations to each molecular degree of freedom, generating a Born-Oppenheimer potential-energy surface (PES) for the molecule. The PES can then be used to derive thermodynamic properties, which are often desirable for applications in chemical engineering and materials design. It is clear from this process that quantum chemistry applications contain a substantial classical computing component in addition to steps that can be performed using a quantum computer. In order to design efficient work flows that take full advantage of each hardware type, it is critical to consider the entire process so that the high-accuracy electronic energies possible from quantum computing are not squandered in the process of calculating thermodynamic properties. We present a summary of the hybrid quantum-classical work flow to compute thermodynamic properties. This work flow contains many options that can significantly affect the efficiency and the accuracy of the results, including classical optimizer attributes, number of ansatz repetitions, and how the vibrational Schr\"odinger equation is solved to determine vibrational modes. We also analyze the effects of these options by employing robust statistics along with simulations and experiments on actual quantum hardware. We show that, through careful selection of work flow options, nearly order-of-magnitude increases in accuracy are possible at equivalent computing time.
Quantum chemistry applications on quantum computers currently rely heavily on the variational quantum eigensolver (VQE) algorithm. This hybrid quantum-classical algorithm aims at finding ground-state solutions of molecular systems based on the variational principle. VQE calculations can be systematically implemented for perturbations to each molecular degree of freedom, generating a Born-Oppenheimer potential-energy surface (PES) for the molecule. The PES can then be used to derive thermodynamic properties, which are often desirable for applications in chemical engineering and materials design. It is clear from this process that quantum chemistry applications contain a substantial classical computing component in addition to steps that can be performed using a quantum computer. In order to design efficient work flows that take full advantage of each hardware type, it is critical to consider the entire process so that the high-accuracy electronic energies possible from quantum computing are not squandered in the process of calculating thermodynamic properties. We present a summary of the hybrid quantum-classical work flow to compute thermodynamic properties. This work flow contains many options that can significantly affect the efficiency and the accuracy of the results, including classical optimizer attributes, number of ansatz repetitions, and how the vibrational Schr"odinger equation is solved to determine vibrational modes. We also analyze the effects of these options by employing robust statistics along with simulations and experiments on actual quantum hardware. We show that, through careful selection of work flow options, nearly order-of-magnitude increases in accuracy are possible at equivalent computing time.
New approaches into computational quantum chemistry can be developed through the use of quantum computing. While universal, fault-tolerant quantum computers are still not available, and we want to utilize today's noisy quantum processors. One of their flagship applications is the variational quantum eigensolver (VQE)-an algorithm for calculating the minimum energy of a physical Hamiltonian. In this study, we investigate how various types of errors affect the VQE and how to efficiently use the available resources to produce precise computational results. We utilize a simulator of a noisy quantum device, an exact statevector simulator, and physical quantum hardware to study the VQE algorithm for molecular hydrogen. We find that the optimal method of running the hybrid classical-quantum optimization is to: (i) allow some noise in intermediate energy evaluations, using fewer shots per step and fewer optimization iterations, but ensure a high final readout precision; (ii) emphasize efficient problem encoding and ansatz parametrization; and (iii) run all experiments within a short time-frame, avoiding parameter drift with time. Nevertheless, current publicly available quantum resources are still very noisy and scarce/expensive, and even when using them efficiently, it is quite difficult to perform trustworthy calculations of molecular energies.
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.
Quantum computing promises to efficiently and accurately solve many important problems in quantum chemistry which elude classical solvers, such as the electronic structure problem of highly correlated materials. Two leading methods in solving the ground state problem are the Variational Quantum Eigensolver (VQE) and Adiabatic Quantum Computing (AQC) algorithms. VQE often struggles with convergence due to the energy landscape being highly non-convex and the existence of barren plateaux, and implementing AQC is beyond the capabilities of current quantum devices as it requires deep circuits. Adiabatically-inspired algorithms aim to fill this gap. In this paper, we first present a unifying framework for these algorithms and then benchmark the following methods: the Adiabatically Assisted VQE (AAVQE) (Garcia-Saez and Latorre (2018)), the Variational Adiabatic Quantum Computing (VAQC) (Harwood et al (2022)), and the Adiabatic Quantum Computing with Parametrized Quantum Circuits (AQC-PQC) (Kolotouros et al (2025)) algorithms. Second, we introduce a novel hybrid approach termed G-AQC-PQC, which generalizes the AQC-PQC method, and combines adiabatic-inspired initialization with the low-memory BFGS optimizer, reducing the quantum computational cost of the method. Third, we compare the accuracy of the methods for chemistry applications using the beryllium hydride molecule (BeH2). We compare the approaches across a number of different choices (ansätze types, depth, discretization steps, initial Hamiltonian, adiabatic schedules and method used). Our results show that the G-AQC-PQC outperforms conventional VQE. We further discuss limitations such as the zero-gradient problem and identify regimes where adiabatically-inspired methods offer a tangible advantage for near-term quantum chemistry applications.
Exploring many-body chemical systems on classical computers often involves solving the Schrödinger equation. However, this approach is frequently limited by the exponential increase in the dimensionality of the Hamiltonian as the number of degrees of freedom increases. In contrast, quantum computing, specifically through the variational quantum eigensolver (VQE) framework, shows promise in overcoming this exponential cost. VQE can utilize the collective properties of quantum states to model the wavefunction in polynomial time. Despite the current limitations of quantum hardware, significant advances have been made in the development of VQE-based algorithms. In this review, we provide an overview of emerging protocols, focusing on their applications in simulating the ground state, excited state, and vibrational properties of chemical systems. By examining notable algorithmic advancements and applications, this review aims to shed light on the challenges and potential of VQE-based algorithms in addressing relevant chemical problems.
In this work, we present a comprehensive comparative study between two different Ansatz quantum circuits , specifically the Unitary Vibrational Coupled Cluster (UVCC) and the Compact Heuristic Circuit (CHC) with use in computational chemistry . By incorporating these circuits into the Variational Quantum Eigensolver (VQE), we aim to find the vibrational ground state energy of small molecules. Our results, obtained using both the UVCC and CHC ansatzes with and without quantum hardware noise, are benchmarked against classical computational methods for energy state determination. We emphasize the efficiency and accuracy of the CHC Ansatz , which significantly reduces circuit complexity without sacrificing the fidelity of the results. This characteristic demonstrates the potential of the CHC Ansatz for scalable quantum chemistry applications, particularly in the Noisy Intermediate-Scale Quantum (NISQ) era. Furthermore, we employ the CHC Ansatz in conjunction with the Variational Quantum Deflation (VQD) algorithm to determine the excited vibrational state energies, which is compared against the reliable method of quantum Equation of Motion(qEOM). Through this comparative study, we aim to provide valuable insights into the practical applications of these Ansatz in quantum simulations of molecular systems, highlighting their respective advantages and potential for future quantum chemistry research.
Variational quantum eigensolvers (VQEs) are among the most promising quantum algorithms for solving electronic structure problems in quantum chemistry, particularly during the noisy intermediate-scale quantum (NISQ) era. In this study, we investigate the capabilities and limitations of VQE algorithms implemented on current quantum hardware for determining molecular ground-state energies, focusing on the adaptive derivative-assembled pseudo-Trotter ansatz VQE (ADAPT-VQE). To address the significant computational challenges posed by molecular Hamiltonians, we explore various well known strategies to simplify the Hamiltonian, optimize the ansatz, and improve classical parameter optimization through modifications of the COBYLA optimizer. These enhancements are integrated into a tailored quantum computing implementation designed to minimize the circuit depth and computational cost. Using benzene as a benchmark system, we demonstrate the application of these optimizations on an IBM quantum computer. Despite these improvements, our results highlight the limitations imposed by current quantum hardware, particularly the impact of quantum noise on state preparation and energy measurement. The noise levels in today's devices prevent meaningful evaluations of molecular Hamiltonians with sufficient accuracy to produce reliable quantum chemical insights. Finally, we extrapolate the requirements for future quantum hardware to enable practical and scalable quantum chemistry calculations using VQE algorithms. This work provides an assessment of current quantum algorithms for molecular modeling on real quantum hardware, highlighting the impact of noise and hardware limitations on the achievable accuracy.