Orbital-optimized pair-correlated electron simulations on trapped-ion quantum computers

Variational quantum eigensolvers (VQE) are among the most promising approaches for solving electronic structure problems on near-term quantum computers. A critical challenge for VQE in practice is that one needs to strike a balance between the expressivity of the VQE ansatz versus the number of quantum gates required to implement the ansatz, given the reality of noisy quantum operations on near-term quantum computers. In this work, we consider an orbital-optimized pair-correlated approximation to the unitary coupled cluster with singles and doubles (uCCSD) ansatz and report a highly efficient quantum circuit implementation for trapped-ion architectures. We show that orbital optimization can recover significant additional electron correlation energy without sacrificing efficiency through measurements of low-order reduced density matrices (RDMs). In the dissociation of small molecules, the method gives qualitatively accurate predictions in the strongly-correlated regime when running on noise-free quantum simulators. On IonQ's Harmony and Aria trapped-ion quantum computers, we run end-to-end VQE algorithms with up to 12 qubits and 72 variational parameters - the largest full VQE simulation with a correlated wave function on quantum hardware. We find that even without error mitigation techniques, the predicted relative energies across different molecular geometries are in excellent agreement with noise-free simulators.


I. INTRODUCTION
Finding accurate solutions to the electronic structure problem is of great importance to various industries, from modeling pharmaceutical drug docking [1], to designing new materials for light harvesting and CO 2 reduction [2], to elucidating reaction mechanisms in novel battery materials [3].However, the classical computational resources needed to solve the electronic structure problem exactly scales exponentially with the size of systems, which limits routine or practical application to systems with less than 20 electrons.To make the problem tractable on classical computers, various approximate approaches have been developed, each with different trade-offs between cost and accuracy.These approaches include density functional theory (DFT) [4], coupled cluster (CC) [5] methods, density matrix renormalization group methods (DMRG) [6], and quantum Monte Carlo methods (QMC) [7].These methods are routinely applied toward computational chemistry calculations both in academia and in industry.
Despite the abundance of different classical approxima-tions, the electronic structure problem is far from being solved.For example, systems with strongly correlated electronic structure are notoriously difficult to solve.These systems are commonly encountered during bond breaking and formation, as well as when studying systems such as transition-metal-containing catalysts, large π-conjugated systems, and high-temperature superconductors.In these cases, approximate approaches may either fail completely (such as in single-reference methods like DFT or CC), or will be prohibitively costly (such as in multi-reference methods like DMRG or QMC).It is possible that approximate classical approaches will never reliably solve the strong correlation problem.
In contrast, quantum computation [8] has attracted significant attention for its potential to solve certain computational problems more efficiently than with classical computers, especially since IBM launched the first cloud accessible quantum computer and Google demonstrated quantum advantage [9].One of its most promising applications is to solve electronic structure problems efficiently [10]: to illustrate, consider that for a problem containing N spin orbitals, the number of classical bits required to represent the wave function scales combinatorially with N, while on a quantum computer only N qubits are needed.The exponential advantage offered by quantum computers has motivated a great deal of research in developing quantum algorithms to solve the electronic structure problem.
Of these, the variational quantum eigensolver (VQE) algorithm [11][12][13][14][15] is designed specifically for current near-term intermediate scale quantum (NISQ) computers.VQE estimates the ground state of a system by implementing a shallow parameterized circuit, which is classically optimized to variationally minimize the energy expectation value.The VQE algorithm allows the user to select the form of the parameterized circuit.This flexibility allows one to adjust the circuit depth based on the quantum gate fidelity, number of qubits, and desired accuracy.This makes VQE especially suitable for the NISQ era.
There is, however, no free lunch and the ability to run shallower circuits within the VQE comes with two costs.First, the predicted energy in most cases remains approximate, because the accuracy depends on the expressivity of the circuit form.Second, one needs to perform a large number of measurements for VQE.This makes the choice of the ansatz perhaps the most important building block in the VQE algorithm.So what does one choose?Early demonstrations of VQE on quantum hardware utilized the physicallymotivated unitary coupled cluster with singles and doubles (uCCSD) ansatz [16,17].uCCSD is well-known in the quantum chemistry community to be able to treat strongly correlated systems, while remaining classically intractable.As such, A. Peruzzo et al. [11] used the uCCSD ansatz in the first VQE demonstration on a photonic quantum computer to solve for the H 2 molecule in a minimal basis.O'Malley and coworkers [12] performed the same simulations on a superconducting quantum computer with two qubits.In 2019, McCaskey [18] and co-workers simulated metal hydrides in a 2-electron, 2-orbital active space using the uCCSD ansatz on IBM's superconducting quantum computers with four qubits.However, going beyond a minimal active space poses difficulties due to the rapid increase in the number of entangling gates for the uCCSD ansatz.The number of entangling gates in uCCSD (e.g.CNOT ) scales as O(N 4 ), where N is the number of qubits.Even the most efficient implementation of uCCSD circuits contain thousands of entangling gates for small systems [19], which makes it impractical to run on NISQ quantum computers.
Due to the impracticality of the uCCSD ansatz on NISQ quantum computers, hardware efficient ansatzes (HEA) [13,[20][21][22][23] have attracted significant attention.Compared to the uCCSD ansatz, HEAs need significantly shallower circuits.In 2017, researchers from IBM published the first study [13] of using HEAs on superconducting quantum computers to simulate H 2 , LiH, and BeH 2 with 2, 4, and 6 qubits.However, noise in the quantum processing unit (QPU) led to unphysical behavior in the predicted dissociation curve.In 2019, the same researchers [24] use HEAs to demonstrate quantum error mitigation using the zero-noise extrapolation technique with 4 qubits.In 2021, researchers [25] used HEAs to study thermally activated delayed fluorescence (TADF) with 2 qubits on superconducting quantum computers.They found that without using an unscalable error mitigation approach, even a 2-qubit circuit yields qualitatively inaccurate predictions to the relative energy.
The largest VQE simulation performed on quantum hard-ware so far is the Hartree-Fock (HF) study by Google [14], in which they used a superconducting quantum computer to simulate the HF wave function for hydrogen chains up to 12 qubits and 72 entangling gates.However, the calculations faced a considerable amount of hardware noise, necessitating the use of Hartree-Fock specific error mitigation techniques to achieve sufficiently accurate results.A more recent study by Google [26] simulated a cyclobutene ring on a superconducting quantum computer with up to 10 qubits using paircorrelated wave functions.Here, the ansatz was classically pre-optimized on a simulator, leaving the final energy evaluation to be executed on the quantum device.Despite this, this calculation still required classical error mitigation techniques to achieve reasonable results for the final quantum energy evaluation.Trapped-ion quantum computers have several unique advantages over other currently available quantum computing architectures.First, the gate fidelity for trapped-ion qubits is typically higher than for superconducting qubits, which enables users to run deeper circuits.Second, trapped-ion qubits are all-to-all connected.This means one is able to entangle arbitrary pairs of qubits without performing expensive SWAP operations to entangle non-adjacent qubits, which is usually required on systems with interactions that do not form a complete graph.Although both of these advantages should lead to higher fidelity when running VQE circuits, implementations of VQE on trapped-ion quantum computers are rare, and this is mainly due to comparatively limited availability of the trapped ion quantum computing hardware versus superconducting quantum computers.In this work, we fill this gap by performing VQE simulations on two generations of trapped-ion quantum computers constructed by IonQ, Inc.
We consider an approximate ansatz derived from the uCCSD ansatz: the unitary pair CCD (upCCD) ansatz, in which only paired double excitations are retained.This allows us to map the fermionic representation to electron pairs, known as the hard-core boson representation.From this, we show that the unitary pCCD ansatz then requires half the number of qubits to encode the state vector as compared to the uCCSD ansatz.We then introduce the optimal circuit for implementing an arbitrary electron pair excitation in terms of the number of CX gates.The energy expression for the up-CCD ansatz is derived, and we find that at most 3 circuits are needed to compute the energy expectation value, regardless of the size of the system.The shallow circuit structure, along with a constant low number of measurements required, make the upCCD ansatz a perfect candidate on NISQ quantum computers.
The accuracy of the upCCD ansatz depends on the choice of the underlying orbitals.Previous studies on similar wave functions have found that it is necessary to optimize the orbitals along with the cluster amplitudes, especially for strongly correlated systems.In this work, we find that the orbital optimization effects can be incorporated through classical post-processing, and only requires the measurements of one-and two-body reduced density matrices (RDM) of the upCCD ansatz on the quantum device.In our experiments, we observe that failure to use orbital optimization results in highly non-physical energy predictions in the bond-dissociation regime, but physical behavior can be fully recovered by optimizing orbitals together with parameters in the upCCD ansatz.Due to the symmetry of the upCCD ansatz, the energy measurements automatically yield the required measurements for RDMs.This allows us to improve the expressivity of the ansatz, especially for strongly correlated systems without increasing the circuit depth on the quantum computer.
In Table I, we list a collection of techniques used in the study, with inventions in this work marked in italic.Our result (see Table II) is the largest full VQE demonstration on a QPU using a correlated wave function without error mitigation.
The paper is structured as follows.We begin by introducing the upCCD ansatz, then discuss the mapping from electron pairs to Pauli matrices, and an efficient quantum circuit implementation of the ansatz.We then derive the energy expression for the upCCD ansatz.Having laid out the general formalism, we then introduce the orbital optimization of up-CCD using RDMs and the Newton-Raphson algorithm.Results are presented on quantum simulators and IonQ's Harmony and Aria quantum computers for potential energy surface predictions of LiH, H 2 O, and Li 2 O molecule.All the VQE experiments on simulator and quantum computers are end-to-end VQE runs, which means we perform both parameter optimizations and final energy evaluations.We conclude with a summary of our findings and comments on future directions.
Readers are strongly encouraged to read "Methods" section before "Results" section.In "Methods" section, we describe the details of the quantum computer hardware used to run VQE and the specifics of the molecular models used to generate the quantum simulation circuits.We heavily use the notations defined in "Methods" section throughout "Results" section.

A. The VQE Algorithm and Circuit
The unitary pair Coupled Cluster double (upCCD) ansatz is in which T is the pair-double cluster operator, defined as in which i and a are indices for occupied and unoccupied orbitals in the HF state.a † pα (a † pβ ) and a pα (a pβ ) are the Fermionic creation and annihilation operators in the pth spin up (down) orbital.
Each exponentiation of the pair-excitation operator can be efficiently implemented with the following circuit, Once the circuit is defined, one needs to measure the energy expectation value Ψ upCCD |H|Ψ upCCD for the secondquantized Hamiltonian H. Originally, there are O(N 4 ) terms in H, in which N is the number of qubits.However, a majority of them do not contribute to the energy since they break pair symmetry.After eliminating these terms, one finds that only 3 measurements are needed in the X, Y , and Z basis respectively to compute the energy, regardless of the system size.
The upCCD ansatz defined in Equation 1is not invariant to the choice of underlying orbitals.Previous studies [39][40][41][42] on similar wave functions have found that it is necessary to optimize the orbitals along with the cluster amplitudes, especially for strongly correlated systems.The orbital optimized upCCD ansatz is in which there are two different sets of parameters: 1) circuit parameters in the cluster operator T ; 2) orbital rotation parameters in the the orbital rotation operator K, which is defined as where K is an anti-Hermitian matrix and σ indexes the spin.
As shown in the Methods section, we find the optimal set of orbital rotation parameters K pq with the Newton-Raphson (NR) algorithm, in which the energy gradient and Hessian are measured on the quantum computer.Then, the effect of orbital rotation operators can be fully absorbed into one-and two-electron integrals through integral transformation, which is done efficiently on classical computers.In this way doing orbital optimization does not increase the circuit depth or the number of measurements.

B. Experimental Example
All the calculations and experiments were performed using IonQ's in-house quantum chemistry library, which facilitates the preparation and execution of quantum variational algorithms on IonQ's cloud simulators and QPUs.We used the PySCF [43] software suite to compute the molecular integrals necessary to define the second-quantized Hamiltonian, as well as compute the classical FCI reference energies.

LiH
We begin with our bond dissociation results with a simple example, the LiH molecule.The system has only two valence electrons.We freeze the Li 1s orbital, and also exclude the molecular orbitals formed with Li's 2px and 2py orbitals  since they do not contribute to the correlation energy due to symmetry.By doing so we only need 3 qubits and the VQE circuits consists of only 4 CX gates.In Figure 1, we compare the energy predicted by FCI, upCCD, and oo-upCCD.
As shown in the plot, both upCCD and oo-upCCD result in accurate energy predictions when the molecule is in the equilibrium geometry.However, as we enter the stretched region, upCCD gives highly nonphysical predictions.Not only does the energy error increases to tens of millihartrees, but it also exhibits a "hump" in the potential energy surface (PES).Such a nonphysical behavior is primarily due to the break down of the mean-field picture in stretched geometries.To the contrary, the oo-upCCD energy prediction matches FCI in both equilibrium and the stretched regions, which demonstrates the importance of orbital optimization.
We then move from simulators to quantum hardware.In Figure 2, we show the results obtained from IonQ's Harmony quantum computer.The system has 11 all-to-all connected qubits, and the averaged single and two-qubit gate fidelities are 99% and 98%.It has been used in numerous applications, including quantum chemistry [15,34], quantum machine learning [44,45], and finance [46,47].Due to the limited machine time, instead of scanning the entire PES, we selected a few points from squeezed, equilibrium, and stretched geometries.As shown in the plot, the energy measured on noisy quantum hardware is much higher than the simulation results.However, we also find that the amount of error is consistent along the PES.Based on such an observation, we shifted all the measured energies so that the energy at R = 1.6 Angstrom matches the simulation energy.By do-ing so, the shifted energies matches well with the exact energy.This is notable especially with the stretched geometry R = 3.0 Angstrom, in which the shifted energy falls on the simulated PES of oo-upCCD, demonstrating that the orbital optimization effects are successfully captured by the quantum hardware.
Lastly, we ran the same simulation on IonQ Aria: IonQ's latest generation quantum computer and the results are shown in Figure 2. IonQ Aria offers both more qubits and improved gate fidelities over IonQ Harmony.We find that the improved gate fidelities reduces the amount of error in energy by 38%.Once shifted, the relative energy also matches the exact energy within statistical uncertainty.The improvements from Harmony to Aria are not very large in this case due to the simplicity of the circuit, which contains only 4 CX gates.

H 2 O
Our next example is the symmetric double dissociation of H 2 O, as shown in Figure 3 for results obtained on the simulator.We only freeze the O 1s core orbital and keep all other orbitals in the active space.The total number of qubits required is 6 and there are 16 CX gates in the circuit.Again, oo-upCCD produces highly accurate energies compared with FCI.However, unlike LiH, in which oo-upCCD matches FCI exactly, in H 2 O we find the predicted energy error for oo-upCCD is about 20 millihartrees, especially when we are in the stretched geometry.The error is due to the omission of the un-paired excitations in the oo-upCCD ansatz.However, we also note that without orbital optimization, the upCCD ansatz using the HF orbitals yields more than 200 millihartrees of error in energy, again emphasizing the importance of orbital optimization.
Before running the circuit on quantum hardware, we first remove redundant parameters from the ansatz.The redundant parameters are the circuit parameters that do not contribute 0.75 1.00 1. 25  to the energy, and their amplitudes stay zero during the optimization process.For the H 2 O molecule, an example of redundant parameters are the amplitudes that correspond to pair excitations from the non-bonding orbital.In this study, we identify redundant parameters by tracking the evolutions of parameter amplitudes on a noise-free simulator, with all parameters started from zero.Parameters whose amplitudes stay at zero during the entire optimization process are identified as redundant parameters.It is worth noting that such an approach does not scale as the system size, and the running time on simulator becomes prohibitively expensive.Fortunately, there exist scalable approaches for identifying and simulating only non-redundant parameters, such as the gradient based selection used for the ADAPT-VQE [17] method.Upon removal of redundant parameters, we are able to reduce the circuit to contain 4 circuit parameters and 8 CX gates.We then performed the oo-upCCD simulation on IonQ's Aria quantum computer, and the results are shown in Figure 4.The simulation is done on two geometry points: one at the equilibrium geometry and the other one at the stretched geometry.We find that the Aria quantum computer successfully finds the optimal parameters and captures the orbital optimization effects.Similar to LiH, the noise on the hardware introduce a systematic, positive bias to the measured absolute energy, but such a bias is constant at different geometry points.Once we shift the absolute energies by a constant, the energies match the ones measured on a noise-free simulator, which demonstrates that the hardware noise is consistent enough so that the predicted relative energies are accurate.

Li 2 O
Our final example is the symmetric dissociation of the Li 2 O molecule.Li 2 O is one of the secondary reaction products in lithium-air batteries, which is believed to be a candidate for next-generation lithium battery due to its high energy den-0.751.00 1. 25  We then performed the oo-upCCD simulation on the Aria quantum computer.Analogous to H 2 O, we first identify and remove redundant parameters.In this example we find that 1.00 1. 25

III. DISCUSSION
Quantum computers are expected to be able to efficiently solve the electronic structure problem. in principle, the electronic energy can be exactly computed in polynomial time using quantum phase estimation (QPE) [48] or its iterative variant [49].In contrast, the best equivalent classical algorithm (full configuration interaction, or FCI) scales exponentially.In QPE, one implements the time propagator U = exp(−iHt) on the quantum computer and operates it on an efficiently prepared trial state.Assuming the trial state has sufficient overlap with the exact eigenstate |Ψ i , the exact eigenstate's energy is encoded in the phase of the wave function since U |Ψ i = exp(−iE i t) |Ψ i .The phase can be extracted using the quantum Fourier transform (QFT).
While the QPE algorithm can compute the energy levels of molecules exactly, it is impractical on current NISQ computers.In the NISQ era, quantum gates are noisy, and entangling gates are typically an order of magnitude lower in fidelity compared to single qubit gates.This means that one can only perform a limited number of quantum operations to ensure that the results are distinguishable from noise.This poses a significant difficulty for the QPE algorithm, as the implementation of the time propagator is very expensive and yields deep quantum circuits.QPE algorithms without using the time propagator, such as qubitization [50][51][52] have also been developed with improved scaling, but the fact remains that neither algorithm results in circuits that are shallow enough to run on quantum computers without error-corrected qubits.
We therefore focus on VQE, an algorithm expressly designed for NISQ computers.Here, we have developed an efficient VQE algorithm that is able to run on near-term quantum computers with high accuracy.The algorithm employs a chemically-inspired ansatz based on the unitary pair coupled cluster doubles (upCCD) wave function.The upCCD ansatz is obtained from the general unitary CCSD ansatz by retaining only paired double excitations.This allows us to condense electron pairs to the hard-core boson representation and develop an efficient quantum circuit implementation that only requires 2 CX gates to implement one excitation.Since the accuracy of the upCCD ansatz depends on the underlying orbital choice, we developed an orbital optimization algorithm that finds the variationally optimal set of orbitals automatically.We find that orbital optimization can be implemented efficiently by measuring one-and two-body RDMs on a quantum computer and computing integral transformations on a classical computer.
We tested the oo-upCCD VQE approach on the bond dissociation pathways for LiH, H 2 O, and Li 2 O molecule on both quantum simulators and IonQ's Harmony and Aria quantum computers.We find that on quantum simulators, oo-upCCD gives qualitatively accurate predictions to energy both in the weakly correlated and strongly correlated regime.However, upCCD without orbital optimization produces unphysical behavior in the strongly correlated regime.On quantum hardware, we observed that noisy quantum gate operations yield a consistent positive bias for the energy.Such a consistent bias has also been observed before [26].In order to understand this, we have performed simulations with both coherent and incoherent noise models on quantum simulators.We find that if the error rate is low enough (below 1%), both error models produce a constant additive error for different molecular geometries, which aligns with the error rate of the Harmony and Aria quantum computers.The simulation results can be found in the supplementary information.Therefore, although the measured absolute energies can be higher than simulator results by hundreds of millihartrees, the relative energies measured are accurate due to the consistency of errors across the PES (i.e., low non-parallelity error).
As with other seniority zero approaches, oo-upCCD proves effective for describing some strong electron correlations but is unable to deliver quantitative accuracy, a difficulty that may in future be addressed in two different ways.First, one may consider implementing the full unitary-CCSD ansatz with quantum circuits, and pay the price of ending up with very deep circuits that are not practical to run on NISQ devices, even with the most efficient compilation techniques.A more practical way is to trade-off circuit depths with measurements and employ approaches such as the quantum subspace expansion (QSE) [37,53].QSE will be able to solve two problems at the same time: 1) account for correlations contributed from non-bosonic excitations.and 2) account for correlations contributed from orbitals that are outside of the active space.QSE is able to achieve these two goals without increasing circuit depth, by just performing more measurements to com-pute higher order RDMs.
In order to achieve quantitative accuracy on a noisy quantum computer using VQE, one would inevitably perform some form of error mitigation.Over the past few years, various error mitigation methods have been developed, such as noise extrapolation [24], density matrix purification [14,18], symmetry verification [54], randomized compiling [55], and noise-estimation.[56] We believe that an efficient VQE approach combined with measurement based post-processing and noise estimation is a very promising route that harvests the most performance out of near-term quantum computers.Together with continued improvements in quantum hardware, both in terms of qubit number and qubit quality, we will soon see quantum simulation of molecules and materials that surpasses the best classical supercomputers.The experimental demonstration was performed on two generations of quantum processing units (QPU) from IonQ: Harmony and Aria.Both QPUs utilize trapped Ytterbium ions where two states in the ground hyperfine manifold are used as qubit states.These states are manipulated by illuminating individual ions with pulses of 355 nm light that drive Raman transitions between the ground states defining the qubit.By configuring these pulses, arbitrary single qubit gates and Mølmer-Sørenson type two-qubit gates can both be realized.The IonQ Aria QPU (schematic in Figure 7) features not only an order of magnitude better performance in terms of fidelity but also is considerably more robust compared to the IonQ Harmony QPU [57].

B. upCCD Circuit Design
From the electron pair excitation operators, we can define the pair creation and annihilation operators in which a i and a † i are the fermionic annihilation and creation operators on orbital i. α and β indicate spin up and spin down.
They follow bosonic symmetries By performing the Jordan-Wigner Transformation (JWT) to map molecular orbitals to qubits, the pair creation and annihilation operators becomes The pair-excitation operator becomes As one can see from the above equation, after JWT, the pair excitation operator does not have the Pauli-Z strings that occur in general double excitations, due to its bosonic nature.
The exponential of the pair-excitation operator, subtracted by its complex conjugate, becomes One can then show that this is the Givens rotation matrix As has been shown before [27], the Givens rotation matrix belongs to the SO(4) group, which can be implemented in 12 elementary (i.e.R y , R z ) gates and 2 CX gates using the magic gate basis

S S H
The efficient Givens rotation implementation with angle θ is the following circuit.

S
Ry(θ ) in which only two CX gates are required.

C. Hamiltonian and Energy Measurements
Since the upCCD ansatz conserves electron pairs, the terms in the ab initio Hamiltonian that break electron pairs do not contribute to energy.After removing these terms, the Hamiltonian can be written as in which the first term only depends on the number operator and so it can be measured in the computational basis.
The second term only depends on the pair excitation operator defined in Equation 8. Furthermore, we note that the two middle terms in it are associated with purely imaginary coefficients, which do not contribute to energy, so that this term can be measured with all qubits in either the X or Y basis.In summary, only 3 circuits are needed to be run in order to measure the energy expectation value for the upCCD ansatz, compared with a number of measurements that scales as O(N 4 ) (where N is the number of orbitals) if no symmetry is exploited, independent of the size of the system.This makes the upCCD ansatz extremely efficient in terms of number of measurements.

D. Orbital Optimization based on Measurements
The orbital optimization effects can be performed classically with integral transformation.Consider the energy expectation value of the oo-upCCD ansatz.
where K is an anti-Hermitian operator defined in Equation 4. We first organize the elements of the lower triangle of K into the length-d vector κ, where Starting from initial orbitals ( κ = 0), we expand the energy out to second order in κ to obtain where the length-d energy gradient ω and the d × d energy Hessian Q are given by which in turn are functions of the spinless one-and twoelectron reduced density matrices (RDM) Since the spinless RDMs are in the form of expectation values, they can be measured on the quantum computer, and since we only need 1-and 2-RDMs, the cost for measuring them is the same as measuring the energy.Using ω and Q, we can choose a κ that reduces the energy using the Newton-Raphson (NR) method, At this point, if we continue with more NR steps until the energy is minimized, we need to implement exp(K) using quantum circuits, which increases the circuit depth.However, we can instead reset κ to zero by absorbing its effects into oneand two-electron integrals through standard molecular orbital transformation.At this point, another NR step can be taken, and the method can be iterated to convergence.In this way, since κ is always zero, we do not need to implement it with quantum circuits.The VQE algorithm is shown in Algorithm 1, and the detailed expressions for the orbital gradients and Hessians in terms of 1-and 2-RDMs, as well as an example of VQE convergence with respect to optimization iterations, can be found in the supplementary information.
Algorithm 1: VQE Algorithm for oo-upCCD while Energy not converged do Optimize T using SPSA Update circuit parameters based on the optimal T * Optimize K using NR Update orbital parameters based on the optimal K * Rotate one-and two-electron integrals to the new orbital basis Set K to be 0 Compute energy to check convergence end

V. DATA AVAILABILITY
The data presented in this manuscript are available from the corresponding author upon reasonable request.and an additional term (S6) For the first one, one can measure the pth qubit in the Z basis, and multiply the measured results by 2. For the second term, one can measure both the pth and the qth qubit in the Z basis, and if both measurements give 1, this term is 2, otherwise this term is 0. For the third term, we also measure the pth and the qth qubit, and if both measurements give 1, this term is 1, and 0 otherwise.
Then the simple part of the Hamiltonian becomes (pq|pq) (X p X q + iX p Y q − iY p X q +Y p Y q ) = 1 2 ∑ p>q (pq|pq) (X p X q +Y p Y q ) (S9) which can be measured by measuring qubit p and q and X or Y basis.As one can see, only 3 circuits are needed to be run in order to measure the energy expectation value for the upCCD ansatz, and it does not depend on the size of the system.This makes the upCCD ansatz extremely efficient in terms of number of measurements.

C. Analytical Expressions for Orbital Gradient and Hessian
For completeness, we include here expressions for the oo-upCCD orbital rotation gradient and hessian.These should provide everything needed for the Newton-Raphson algorithm we use for orbital optimization.
The energy can be written as

FIG. 2 :
FIG. 2: Dissociation of LiH in STO-3G basis set comparing upCCD, oo-upCCD, and FCI.VQE results are obtained from the IonQ Harmony and the IonQ Aria quantum computer.

FIG. 7 :
FIG. 7: Schematic of the IonQ Aria quantum computer, where two fans of individually addressable beams illuminate a chain of individually imaged ions.

H 2 =
∑ p =q (pq|pq)a † pα a qα a † pβ a qβ (S5) For the first part of the Hamiltonian, we need to compute the following expectation values z p = n pα + n pβ Γ pq = n pα n qα + n pβ n qβ ∆ pq = n pα n qβ .

TABLE I :
Techniques and their corresponding effects used in the paper.

TABLE II :
Comparison between this work and previous publicly-reported VQE chemistry simulations on QPUs.
We freeze the 1s orbital for Li and O, resulting in a circuit with 12 qubits and 64 CX gates.The results on an ideal simulator are shown in Figure5.The difference in energy between oo-upCCD and FCI becomes more noticeable, which is expected since Li 2 O has twice (four times) the number of orbitals as H 2 O (LiH).As a result, there are many more unpaired excited configurations in Li 2 O than H 2 O and LiH, and ignoring them, as is being done in the upCCD ansatz, introduces a more drastic approximation.Again, we find that orbital optimization does not make any noticeable amount of difference in equilibrium geometry, but becomes crucial in stretched geometries.
CX gates) on the quantum hardware, with an additional 66 orbital rotation parameters, for 72 variational parameters total.The results are shown in Figure6at two geometry points: one at the equilibrium geometry and another one at the stretched geometry.Like in H 2 O, we find that despite hardware noise, the predicted relative energy matches the simulator's prediction, and the orbital optimization effects are successfully captured by Aria.