We present example quantum chemistry programs written with JaqalPaq, a python meta-programming language used to code in Jaqal (Just Another Quantum Assembly Language). These JaqalPaq algorithms are intended to be run on the Quantum Scientific Computing Open User Testbed (QSCOUT) platform at Sandia National Laboratories. Our exemplars use the variational quantum eigensolver (VQE) quantum algorithm to compute the ground state energies of the H2, HeH+, and LiH molecules. Since the exemplars focus on how to program in JaqalPaq, the calculations of the second-quantized Hamiltonians are performed with the PySCF python package, and the mappings of the fermions to qubits are obtained from the OpenFermion python package. Using the emulator functionality of JaqalPaq, we emulate how these exemplars would be executed on an error-free QSCOUT platform and compare the emulated computation of the bond-dissociation curves for these molecules with their exact forms within the relevant basis.
We present example quantum chemistry programs written with JaqalPaq, a python meta-programming language used to code in Jaqal (Just Another Quantum Assembly Language). These JaqalPaq algorithms are intended to be run on the Quantum Scientific Computing Open User Testbed (QSCOUT) platform at Sandia National Laboratories. Our exemplars use the variational quantum eigensolver (VQE) quantum algorithm to compute the ground state energies of the H2, HeH+, and LiH molecules. Since the exemplars focus on how to program in JaqalPaq, the calculations of the second-quantized Hamiltonians are performed with the PySCF python package, and the mappings of the fermions to qubits are obtained from the OpenFermion python package. Using the emulator functionality of JaqalPaq, we emulate how these exemplars would be executed on an error-free QSCOUT platform and compare the emulated computation of the bond-dissociation curves for these molecules with their exact forms within the relevant basis.
Abstract Quantum chemistry (QC) is one of the most promising applications of quantum computing. However, present quantum processing units (QPUs) are still subject to large errors. Therefore, noisy intermediate‐scale quantum (NISQ) hardware is limited in terms of qubit counts/circuit depths. Variational quantum eigensolver (VQE) algorithms can potentially overcome such issues. Here, we introduce the OpenVQE open‐source QC package. It provides tools for using and developing chemically‐inspired adaptive methods derived from unitary coupled cluster (UCC). It facilitates the development and testing of VQE algorithms and is able to use the Atos Quantum Learning Machine (QLM), a general quantum programming framework enabling to write/optimize/simulate quantum computing programs. We present a specific, freely available QLM open‐source module, myQLM‐fermion. We review its key tools for facilitating QC computations (fermionic second quantization, fermion‐spin transforms, etc.). OpenVQE largely extends the QLM's QC capabilities by providing: (i) the functions to generate the different types of excitations beyond the commonly used UCCSD ansatz; (ii) a new Python implementation of the “adaptive derivative assembled pseudo‐Trotter method” (ADAPT‐VQE). Interoperability with other major quantum programming frameworks is ensured thanks to the myQLM‐interop package, which allows users to build their own code and easily execute it on existing QPUs. The combined OpenVQE/myQLM‐fermion libraries facilitate the implementation, testing and development of variational quantum algorithms, while offering access to large molecules as the noiseless Schrödinger‐style dense simulator can reach up to 41 qubits for any circuit. Extensive benchmarks are provided for molecules associated to qubit counts ranging from 4 to 24. We focus on reaching chemical accuracy, reducing the number of circuit gates and optimizing parameters and operators between “fixed‐length” UCC and ADAPT‐VQE ansätze. This article is categorized under: Software > Quantum Chemistry Quantum Computing > Algorithms
Abstract We present a comprehensive quantum chemical study of nine molecules H₂, HeH⁺, LiH, BeH₂, H₂O, NH₃, CH₄, CO₂, and N₂ using the Variational Quantum Eigensolver (VQE) algorithm implemented in Qiskit 2.4.1. Molecular Hamiltonians are constructed in the STO-3G minimal basis set, transformed to qubit operators via the Jordan-Wigner mapping, and optimized using the n-local hardware-efficient ansatz with the SLSQP optimizer. Ground state energies are benchmarked against Full Configuration Interaction (FCI)/STO-3G reference values from authoritative literature sources. VQE achieves accuracies from 96.46% (HeH⁺) to 99.99% (H₂), average 98.68%. Potential energy surfaces are generated for all nine molecules and compared with high level ab initio literature PES references including the Kolos-Wolniewicz H₂ curve, the Polyansky water surface, and the Bytautas-Ruedenberg N₂ FCI surface. Equilibrium bond lengths agree with experiment to within 0.000–0.006 Å. Quantum ansatz circuits are analyzed in detail for all molecules. These results confirm the viability of near-term VQE algorithms for molecular quantum chemistry across diverse bond types and geometries.
This study investigates the Variational Quantum Eigensolver (VQE), a hybrid quantum-classical algorithm that partitions the computational workload between a quantum processing unit and a classical optimiser to determine the ground-state energies of molecular systems. The Qiskit SDK v1.x framework is employed, utilising the PySCFDriver module within the qiskit-nature package to construct second-quantised molecular electronic Hamiltonians for two representative systems: Hydrogen ( ) and Lithium Hydride ( ). For the complete active space is retained without spatial reduction, mapping the full Hamiltonian onto a 4-qubit register. However for the full STO-3G representation demands 12 qubits, rendering direct classical parameter optimisation computationally intractable. This bottleneck is resolved by applying an ActiveSpaceTransformer (num_electrons=2, num_spatial_orbitals=3) to freeze the chemically inactive lithium core electrons, reducing the active register to 6 qubits. Fermionic operators are subsequently converted to Pauli strings via the Jordan-Wigner transformation ( ) enabling the quantum processor to evaluate the molecular Hamiltonian as a circuit while preserving the fermionic anti-commutation relations. Under ideal, noise-free simulation using the StatevectorEstimator primitive paired with the COBYLA optimiser (maxiter=200), all VQE energy points converge precisely onto the Full Configuration Interaction (FCI) classical baseline across the full range of interatomic separations. This confirms that the hardware-efficient TwoLocal ansatz ( rotation layers and entanglement layers with reps=3) possesses sufficient expressive power to capture the molecular ground-state wavefunction. Introducing gate-level noise through AerSimulator backends calibrated from the FakeManilaV2 and FakeNairobiV2 hardware profiles—with optimisation driven by the SPSA algorithm (maxiter=150)—reveals two critical failure modes. A systematic upward energy bias is observed, attributable to depolarising gate errors and measurement readout faults that prevent the optimiser from reaching the true minimum energy eigenvalue. Additionally, severe stochastic oscillations between 1.0 Å and 1.5 Å on the curve induce premature convergence of the SPSA optimiser to spurious local minima. These results demonstrate that unmitigated hardware noise renders VQE incapable of reliably predicting equilibrium bond lengths on current NISQ devices, motivating the future integration of Zero-Noise Extrapolation (ZNE), readout error mitigation, and Virtual Distillation into the computational pipeline. Keywords—Variational Quantum Eigensolver (VQE); Qiskit 1.x; Potential Energy Curve; Molecular Ground-State Energy; Hardware-Efficient Ansatz; Gate-Level Noise Simulation; ActiveSpaceTransformer; NISQ Devices.
We use the Variational Quantum Eigensolver (VQE) as implemented in the Qiskit software package to compute the ground state energy of small molecules derived from water, H2O, and hydrogen cyanide, HCN. The work aims to benchmark algorithms for calculating the electronic structure and energy surfaces of molecules of relevance to prebiotic chemistry, beginning with water and hydrogen cyanide, and to run them on the available simulated and physical quantum hardware. The numerical calculations of the algorithms for small quantum processors allow us to design more efficient protocols to be run in real hardware, as well as to analyze their performance. Future implementations on accessible quantum processing prototypes will benchmark quantum computers and provide tests of quantum advantage with heuristic quantum algorithms.
TenCirChem is an open-source Python library for simulating variational quantum algorithms for quantum computational chemistry. TenCirChem shows high-performance in the simulation of unitary coupled-cluster circuits, using compact representations of quantum states and excitation operators. Additionally, TenCirChem supports noisy circuit simulation and provides algorithms for variational quantum dynamics. TenCirChem ’s capabilities are demonstrated through various examples, such as the calculation of the potential energy curve of H 2 O with a 6-31G(d) basis set using a 34-qubit quantum circuit, the examination of the impact of quantum gate errors on the variational energy of the H 2 molecule, and the exploration of the Marcus inverted region for charge transfer rate based on variational quantum dynamics. Furthermore, TenCirChem is capable of running real quantum hardware experiments, making it a versatile tool for both simulation and experimentation in the field of quantum computational chemistry.
We present computational chemistry data for small molecules (CO, HCl, F2, NH4+, CH4, NH3, H3O+, H2O, BeH2, LiH, OH−, HF, HeH+, H2), obtained by implementing the Unitary Coupled Cluster method with Single and Double excitations (UCCSD) on a quantum computer simulator. We have used the Variational Quantum Eigensolver (VQE) algorithm to extract the ground state energies of these molecules. This energy data represents the expected ground state energy that a quantum computer will produce for the given molecules, on the STO-3G basis. Since there is a lot of interest in the implementation of UCCSD on quantum computers, we hope that our work will serve as a benchmark for future experimental implementations.