Scaling up electronic structure calculations on quantum computers: The frozen natural orbital based method of increments

Variational quantum imaginary time evolution for matrix product state Ansatz with tests on transcorrelated Hamiltonians

The matrix product state (MPS) Ansatz offers a promising approach for finding the ground state of molecular Hamiltonians and solving quantum chemistry problems. Building on this concept, the proposed technique of quantum circuit MPS (QCMPS) enables the simulation of chemical systems using a relatively small number of qubits. In this study, we enhance the optimization performance of the QCMPS Ansatz by employing the variational quantum imaginary time evolution (VarQITE) approach. Guided by McLachlan's variational principle, the VarQITE method provides analytical metrics and gradients, resulting in improved convergence efficiency and robustness of the QCMPS. We validate these improvements numerically through simulations of H2, H4, and LiH molecules. In addition, given that VarQITE is applicable to non-Hermitian Hamiltonians, we evaluate its effectiveness in preparing the ground state of transcorrelated Hamiltonians. This approach yields energy estimates comparable to the complete basis set (CBS) limit while using even fewer qubits. In particular, we perform simulations of the beryllium atom and LiH molecule using only three qubits, maintaining high fidelity with the CBS ground state energy of these systems. This qubit reduction is achieved through the combined advantages of both the QCMPS Ansatz and transcorrelation. Our findings demonstrate the potential practicality of this quantum chemistry algorithm on near-term quantum devices.

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.

Linear Scaling Incremental Scheme for Correlation Energies with Embedding Generated Virtuals

A novel incremental scheme is presented including an incremental expansion of the virtual space for the calculation of electron correlation energies, which is compatible with any size-extensive correlation method and scales asymptotically linear for large molecules. The performance is studied for organic molecules, water clusters, and a La(III)-water complex, where the compatibility with pseudopotentials is also examined. The computational requirements are already reduced tremendously for medium-sized water clusters and hydrocarbons with respect to the canonical CCSD as well as the ordinary incremental scheme references. Correlation energies within chemical accuracy have been observed for all studied systems. The novelty of the method is that relatively small virtual spaces are used in combination with tuples of localized occupied spaces. The corresponding orthonormal occupied and virtual orbitals are obtained from QM/QM embedding calculations and can thus be used with standard quantum chemistry codes for correlation calculations. It is presented how relevant virtual spaces are selected and the correlation energies are linked in the new virtual space expansion.

Toward Accurate Post-Born-Oppenheimer Molecular Simulations on Quantum Computers: An Adaptive Variational Eigensolver with Nuclear-Electronic Frozen Natural Orbitals

Nuclear quantum effects such as zero-point energy and hydrogen tunneling play a central role in many biological and chemical processes. The nuclear-electronic orbital (NEO) approach captures these effects by treating selected nuclei quantum mechanically on the same footing as electrons. On classical computers, the resources required for an exact solution of NEO-based models grow exponentially with system size. By contrast, quantum computers offer a means of solving this problem with polynomial scaling. However, due to the limitations of current quantum devices, NEO simulations are confined to the smallest systems described by minimal basis sets, whereas realistic simulations beyond the Born-Oppenheimer approximation require more sophisticated basis sets. For this purpose, we herein extend a hardware-efficient ADAPT-VQE method to the NEO framework in the frozen natural orbital (FNO) basis. We demonstrate on H2 and D2 molecules that the NEO-FNO-ADAPT-VQE method reduces the CNOT count by several orders of magnitude relative to the NEO unitary coupled cluster method with singles and doubles while maintaining the desired accuracy. This extreme reduction in the CNOT gate count is sufficient to permit practical computations employing the NEO method─an important step toward accurate simulations involving nonclassical nuclei and non-Born-Oppenheimer effects on near-term quantum devices. We further show that the method can capture isotope effects, and we demonstrate that inclusion of correlation energy systematically improves the prediction of difference in the zero-point energy (ΔZPE) between isotopes.

Nonadiabatic Nuclear-Electron Dynamics: A Quantum Computing Approach

Coupled quantum electron-nuclear dynamics is often associated with the Born-Huang expansion of the molecular wave function and the appearance of nonadiabatic effects as a perturbation. On the other hand, native multicomponent representations of electrons and nuclei also exist, which do not rely on any a priori approximation. However, their implementation is hampered by prohibitive scaling. Consequently, quantum computers offer a unique opportunity for extending their use to larger systems. Here, we propose a quantum algorithm for simulating the time-evolution of molecular systems and apply it to proton transfer dynamics in malonaldehyde, described as a rigid scaffold. The proposed quantum algorithm can be easily generalized to include the explicit dynamics of the classically described molecular scaffold. We show how entanglement between electronic and nuclear degrees of freedom can persist over long times if electrons do not follow the nuclear displacement adiabatically. The proposed quantum algorithm may become a valid candidate for the study of such phenomena when sufficiently powerful quantum computers become available.

Quantum computing in pharma: A multilayer embedding approach for near future applications

Quantum computers are special purpose machines that are expected to be particularly useful in simulating strongly correlated chemical systems. The quantum computer excels at treating a moderate number of orbitals within an active space in a fully quantum mechanical manner. We present a quantum phase estimation calculation on F₂ in a (2,2) active space on Rigetti's Aspen-11 QPU. While this is a promising start, it also underlines the need for carefully selecting the orbital spaces treated by the quantum computer. In this work, a scheme for selecting such an active space automatically is described and simulated results obtained using both the quantum phase estimation (QPE) and variational quantum eigensolver (VQE) algorithms are presented and combined with a subtractive method to enable accurate description of the environment. The active occupied space is selected from orbitals localized on the chemically relevant fragment of the molecule, while the corresponding virtual space is chosen based on the magnitude of interactions with the occupied space calculated from perturbation theory. This protocol is then applied to two chemical systems of pharmaceutical relevance: the enzyme [Fe] hydrogenase and the photosenzitizer temoporfin. While the sizes of the active spaces currently amenable to a quantum computational treatment are not enough to demonstrate quantum advantage, the procedure outlined here is applicable to any active space size, including those that are outside the reach of classical computation.

Optimized Quantum Phase Estimation for Simulating Electronic States in Various Energy Regimes

While quantum algorithms for simulations exhibit better asymptotic scaling than their classical counterparts, they currently cannot be accurately implemented on real-world devices. Instead, chemists and computer scientists rely on costly classical simulations of these quantum algorithms. In particular, the quantum phase estimation (QPE) algorithm is among several approaches that has attracted much attention in recent years due to its genuine quantum character. However, it is memory-intensive to simulate and intractable for moderate system sizes. This paper discusses the performance and applicability of QPESIM, a new simulation of the QPE algorithm designed to take advantage of modest computational resources. In particular, we demonstrate the versatility of QPESIM in simulating various electronic states by examining the ground and core-level states of H₂O. For these states, we also discuss the effect of the active-space size on the quality of the calculated energies. For the high-energy core-level states, we demonstrate that new QPE simulations for active spaces defined by 15 active orbitals significantly reduce the errors in core-level excitation energies compared to earlier QPE simulations using smaller active spaces.

Quantum Algorithm of the Divide-and-Conquer Unitary Coupled Cluster Method with a Variational Quantum Eigensolver

The variational quantum eigensolver (VQE) with shallow or constant-depth quantum circuits is one of the most pursued approaches in the noisy intermediate-scale quantum (NISQ) devices with incoherent errors. In this study, the divide-and-conquer (DC) linear scaling technique, which divides the entire system into several fragments, is applied to the VQE algorithm based on the unitary coupled cluster (UCC) method, denoted as DC-qUCC/VQE, to reduce the number of required qubits. The unitarity of the UCC ansatz that enables the evaluation of the total energy as well as various molecular properties as expectation values can be easily implemented on quantum devices because the quantum gates are unitary operators themselves. Based on this feature, the present DC-qUCC/VQE algorithm is designed to conserve the total number of electrons in the entire system using the density matrix evaluated on a quantum computer. Numerical assessments clarified that the energy errors of the DC-qUCC/VQE calculations decrease by using the constraint of the total number of electrons. Furthermore, the DC-qUCC/VQE algorithm could reduce the number of quantum gates and shows the possibility of decreasing incoherent errors.

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree-Fock virtual orbitals by natural orbitals (NOs). In this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller-Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties-electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F-Ts), and parity-violating energy differences for H₂Z₂ (Z = O-Se)-we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

Quantum algorithms for electronic structures: basis sets and boundary conditions

The advantages of quantum computers are believed to significantly change the research paradigm of chemical and materials sciences, where computational characterization and theoretical design play an increasingly important role. It is especially desirable to solve the electronic structure problem, a central problem in chemistry and materials science, efficiently and accurately with well-designed quantum algorithms. Various quantum electronic-structure algorithms have been proposed in the literature. In this article, we briefly review recent progress in this direction with a special emphasis on the basis sets and boundary conditions. Compared to classical electronic structure calculations, there are new considerations in choosing a basis set in quantum algorithms. For example, the effect of the basis set on the circuit complexity is very important in quantum algorithm design. Electronic structure calculations should be performed with an appropriate boundary condition. Simply using a wave function ansatz designed for molecular systems in a material system with a periodic boundary condition may lead to significant errors. Artificial boundary conditions can be used to partition a large system into smaller fragments to save quantum resources. The basis sets and boundary conditions are expected to play a crucial role in electronic structure calculations on future quantum computers, especially for realistic systems.

Variational quantum eigensolver simulations with the multireference unitary coupled cluster ansatz: a case study of the C2v quasi-reaction pathway of beryllium insertion to H2 molecule

Variational quantum eigensolver (VQE)-based quantum chemical calculations have been extensively studied as a computational model using noisy intermediate-scale quantum devices. VQE uses a parametrized quantum circuit defined through an “ansatz”...

Cluster many-body expansion: A many-body expansion of the electron correlation energy about a cluster mean field reference

The many-body expansion (MBE) is an efficient tool that has a long history of use for calculating interaction energies, binding energies, lattice energies, and so on. In the past, applications of MBE to correlation energy have been unfeasible for large systems, but recent improvements to computing resources have sparked renewed interest in capturing the correlation energy using the generalized nth order Bethe-Goldstone equation. In this work, we extend this approach, originally proposed for a Slater determinant, to a tensor product state (TPS) based wavefunction. By partitioning the active space into smaller orbital clusters, our approach starts from a cluster mean field reference TPS configuration and includes the correlation contribution of the excited TPSs using the MBE. This method, named cluster MBE (cMBE), improves the convergence of MBE at lower orders compared to directly doing a block-based MBE from a RHF reference. We present numerical results for strongly correlated systems, such as the one- and two-dimensional Hubbard models and the chromium dimer. The performance of the cMBE method is also tested by partitioning the extended π space of several large π-conjugated systems, including a graphene nano-sheet with a very large active space of 114 electrons in 114 orbitals, which would require 10⁶⁶ determinants for the exact FCI solution.