Improving the accuracy of the variational quantum eigensolver for molecular systems by the explicitly-correlated perturbative [2]R12correction

Shortcut to chemically accurate quantum computing via density-based basis-set correction

Using GPU-accelerated state-vector emulation, we propose to embed a quantum computing ansatz into density-functional theory via density-based basis-set corrections to obtain quantitative quantum-chemistry results on molecules that would otherwise require brute-force quantum calculations using hundreds of logical qubits. Indeed, accessing a quantitative description of chemical systems while minimizing quantum resources is an essential challenge given the limited qubit capabilities of current quantum processors. We provide a shortcut towards chemically accurate quantum computations by approaching the complete-basis-set limit through coupling the density-based basis-set corrections approach, applied to any given variational ansatz, to an on-the-fly crafting of basis sets specifically adapted to a given system and user-defined qubit budget. The resulting approach self-consistently accelerates the basis-set convergence, improving electronic densities, ground-state energies, and first-order properties (e.g. dipole moments), but can also serve as a classical, a posteriori, energy correction to quantum hardware calculations with expected applications in drug design and materials science.

Towards efficient quantum computing for quantum chemistry: reducing circuit complexity with transcorrelated and adaptive ansatz techniques

The near-term utility of quantum computers is hindered by hardware constraints in the form of noise. One path to achieving noise resilience in hybrid quantum algorithms is to decrease the required circuit depth - the number of applied gates - to solve a given problem. This work demonstrates how to reduce circuit depth by combining the transcorrelated (TC) approach with adaptive quantum ansätze and their implementations in the context of variational quantum imaginary time evolution (AVQITE). The combined TC-AVQITE method is used to calculate ground state energies across the potential energy surfaces of H4, LiH, and H2O. In particular, H4 is a notoriously difficult case where unitary coupled cluster theory, including singles and doubles excitations, fails to provide accurate results. Adding TC yields energies close to the complete basis set (CBS) limit while reducing the number of necessary operators - and thus circuit depth - in the adaptive ansätze. The reduced circuit depth furthermore makes our algorithm more noise-resilient and accelerates convergence. Our study demonstrates that combining the TC method with adaptive ansätze yields compact, noise-resilient, and easy-to-optimize quantum circuits that yield accurate quantum chemistry results close to the CBS limit.

Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method

Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.

Self-Consistent Field Approach for the Variational Quantum Eigensolver: Orbital Optimization Goes Adaptive

We present a self-consistent field (SCF) approach within the adaptive derivative-assembled problem-tailored ansatz variational quantum eigensolver (ADAPT-VQE) framework for efficient quantum simulations of chemical systems on near-term quantum computers. To this end, our ADAPT-VQE-SCF approach combines the idea of generating an ansatz with a small number of parameters, resulting in shallow-depth quantum circuits with a direct minimization of an energy expression that is correct to second order with respect to changes in the molecular orbital basis. Our numerical analysis, including calculations for the transition-metal complex ferrocene [Fe (C5H5)2], indicates that convergence in the self-consistent orbital optimization loop can be reached without a considerable increase in the number of two-qubit gates in the quantum circuit by comparison to a VQE optimization in the initial molecular orbital basis. Moreover, the orbital optimization can be carried out simultaneously within each iteration of the ADAPT-VQE cycle. ADAPT-VQE-SCF thus allows us to implement a routine analogous to the complete active space SCF, a cornerstone of state-of-the-art computational chemistry, in a hardware-efficient manner on near-term quantum computers. Hence, ADAPT-VQE-SCF paves the way toward a paradigm shift for quantitative quantum-chemistry simulations on quantum computers by requiring fewer qubits and opening up for the use of large and flexible atomic orbital basis sets in contrast to earlier methods that are predominantly based on the idea of full active spaces with minimal basis sets.

Applications of noisy quantum computing and quantum error mitigation to "adamantaneland": a benchmarking study for quantum chemistry

The field of quantum computing has the potential to transform quantum chemistry. The variational quantum eigensolver (VQE) algorithm has allowed quantum computing to be applied to chemical problems in the noisy intermediate-scale quantum (NISQ) era. Applications of VQE have generally focused on predicting absolute energies instead of chemical properties that are relative energy differences and that are most interesting to chemists studying a chemical problem. We address this shortcoming by constructing a molecular benchmark data set in this work containing isomers of C10H16 and carbocationic rearrangements of C10H15+, calculated at a high-level of theory. Using the data set, we compared noiseless VQE simulations to conventionally performed density functional and wavefunction theory-based methods to understand the quality of results. We also investigated the effectiveness of a quantum state tomography-based error mitigation technique in applications of VQE under noise (simulated and real). Our findings reveal that the use of quantum error mitigation is crucial in the NISQ era and advantageous to yield almost noiseless quality results.

Chemically Accurate Potential Curves for H2 Molecules Using Explicitly Correlated Qubit-ADAPT

With the recent advances in the development of devices capable of performing quantum computations, a growing interest in finding near-term applications has emerged in many areas of science. In the era of nonfault tolerant quantum devices, algorithms that only require comparably short circuits accompanied by high repetition rates are considered to be a promising approach for assisting classical machines with finding a solution on computationally hard problems. The ADAPT approach previously introduced in Nat. Commun. 10, 3007 (2019) extends the class of variational quantum eigensolver algorithms with dynamically growing ansätze in order to find approximations to the ground and excited state energies of molecules. In this work, the ADAPT algorithm has been combined with a first-quantized formulation for the hydrogen molecule in the Born-Oppenheimer approximation, employing the explicitly correlated basis functions introduced in J. Chem. Phys. 43, 2429 (1965). By the virtue of their explicit electronic correlation properties, it is shown in classically performed simulations that chemical accuracy (<1.6 mHa) can be reached for ground and excited state potential curves using reasonably short circuits.

Basis Set Generation and Optimization in the NISQ Era with Quiqbox.jl

In the noisy intermediate-scale quantum era, ab initio computation of electronic structure problems has become one of the major benchmarks for identifying the boundary between classical and quantum computational power. Basis sets play a key role in the electronic structure methods implemented on both classical and quantum devices. To investigate the consequences of single-particle basis sets, we propose a framework for more customizable basis set generation and optimization. This framework allows composite basis sets to go beyond typical basis set frameworks, such as atomic basis sets, by introducing the concept of mixed-contracted Gaussian-type orbitals. These basis set generations set the stage for more flexible variational optimization of basis set parameters. To realize this framework, we have developed an open-source software package named "Quiqbox" in the Julia programming language. We demonstrate various examples of using Quiqbox for basis set optimization and generation, ranging from optimizing atomic basis sets on the Hartree-Fock level, preparing the initial state for variational quantum eigensolver computation, and constructing basis sets with completely delocalized orbitals. We also include various benchmarks of Quiqbox for basis set optimization and ab initial electronic structure computation.

Nonunitary projective transcorrelation theory inspired by the F12 ansatz

An alternative nonunitary transcorrelation, inspired by the F12 ansatz, is investigated. In contrast to the Jastrow transcorrelation of Boys-Handy, the effective Hamiltonian of this projective transcorrelation features: 1. a series terminating formally at four-body interactions. 2. no spin-contamination within the non-relativistic framework. 3. simultaneous satisfaction of the singlet and triplet first-order cusp conditions. 4. arbitrary choices of pairs for correlation including frozen-core approximations. We discuss the connection between the projective transcorrelation and F12 theory with applications to small molecules, to show that the cusp conditions play an important role to reduce the uncertainty arising from the nonunitary transformation.

Partitioning Quantum Chemistry Simulations with Clifford Circuits

Current quantum computing hardware is restricted by the availability of only few, noisy qubits which limits the investigation of larger, more complex molecules in quantum chemistry calculations on quantum computers in the near term. In this work, we investigate the limits of their classical and near-classical treatment while staying within the framework of quantum circuits and the variational quantum eigensolver. To this end, we consider naive and physically motivated, classically efficient product ansatz for the parametrized wavefunction adapting the separable-pair ansatz form. We combine it with post-treatment to account for interactions between subsystems originating from this ansatz. The classical treatment is given by another quantum circuit that has support between the enforced subsystems and is folded into the Hamiltonian. To avoid an exponential increase in the number of Hamiltonian terms, the entangling operations are constructed from purely Clifford or near-Clifford circuits. While Clifford circuits can be simulated efficiently classically, they are not universal. In order to account for missing expressibility, near-Clifford circuits with only few, selected non-Clifford gates are employed. The exact circuit structure to achieve this objective is molecule-dependent and is constructed using simulated annealing and genetic algorithms. We demonstrate our approach on a set of molecules of interest and investigate the extent of our methodology's reach.

Multiscale quantum algorithms for quantum chemistry

Exploring the potential applications of quantum computers in material design and drug discovery is attracting more and more attention after quantum advantage has been demonstrated using Gaussian boson sampling. However, quantum resource requirements in material and (bio)molecular simulations are far beyond the capacity of near-term quantum devices. In this work, multiscale quantum computing is proposed for quantum simulations of complex systems by integrating multiple computational methods at different scales of resolution. In this framework, most computational methods can be implemented in an efficient way on classical computers, leaving the critical portion of the computation to quantum computers. The simulation scale of quantum computing strongly depends on available quantum resources. As a near-term scheme, we integrate adaptive variational quantum eigensolver algorithms, second-order Møller-Plesset perturbation theory and Hartree-Fock theory within the framework of the many-body expansion fragmentation approach. This new algorithm is applied to model systems consisting of hundreds of orbitals with decent accuracy on the classical simulator. This work should encourage further studies on quantum computing for solving practical material and biochemistry problems.

Performance of a one-parameter correlation factor for transcorrelation: Study on a series of second row atomic and molecular systems

In this work, we investigate the performance of a recently proposed transcorrelated (TC) approach based on a single-parameter correlation factor [E. Giner, J. Chem. Phys. 154, 084119 (2021)] for systems involving more than two electrons. The benefit of such an approach relies on its simplicity as efficient numerical-analytical schemes can be set up to compute the two- and three-body integrals occurring in the effective TC Hamiltonian. To obtain accurate ground state energies within a given basis set, the present TC scheme is coupled to the recently proposed TC-full configuration interaction quantum Monte Carlo method [Cohen et al., J. Chem. Phys. 151, 061101 (2019)]. We report ground state total energies on the Li-Ne series, together with their first cations, computed with increasingly large basis sets and compare to more elaborate correlation factors involving electron-electron-nucleus coordinates. Numerical results on the Li-Ne ionization potentials show that the use of the single-parameter correlation factor brings on average only a slightly lower accuracy (1.2 mH) in a triple-zeta quality basis set with respect to a more sophisticated correlation factor. However, already using a quadruple-zeta quality basis set yields results within chemical accuracy to complete basis set limit results when using this novel single-parameter correlation factor. Calculations on the H₂O, CH₂, and FH molecules show that a similar precision can be obtained within a triple-zeta quality basis set for the atomization energies of molecular systems.