Hybrid Quantum-Classical Eigensolver Without Variation or Parametric Gates

The use of near-term quantum devices that lack quantum error correction, for addressing quantum chemistry and physics problems, requires hybrid quantum-classical algorithms and techniques. Here we present a process for obtaining the eigenenergy spectrum of electronic quantum systems. This is achieved by projecting the Hamiltonian of a quantum system onto a limited effective Hilbert space specified by a set of computational bases. From this projection an effective Hamiltonian is obtained. Furthermore, a process for preparing short depth quantum circuits to measure the corresponding diagonal and off-diagonal terms of the effective Hamiltonian is given, whereby quantum entanglement and ancilla qubits are used. The effective Hamiltonian is then diagonalized on a classical computer using numerical algorithms to obtain the eigenvalues. The use case of this approach is demonstrated for ground sate and excited states of BeH$_2$ and LiH molecules, and the density of states, which agrees well with exact solutions. Additionally, hardware demonstration is presented using IBM quantum devices for H$_2$ molecule.

The use of near-term quantum devices that lack quantum error correction, for addressing quantum chemistry and physics problems, requires hybrid quantum-classical algorithms and techniques. Here we present a process for obtaining the eigenenergy spectrum of electronic quantum systems. This is achieved by projecting the Hamiltonian of a quantum system onto a limited effective Hilbert space specified by a set of computational bases. From this projection an effective Hamiltonian is obtained. Furthermore, a process for preparing short depth quantum circuits to measure the corresponding diagonal and off-diagonal terms of the effective Hamiltonian is given, whereby quantum entanglement and ancilla qubits are used. The effective Hamiltonian is then diagonalized on a classical computer using numerical algorithms to obtain the eigenvalues. The use case of this approach is demonstrated for ground sate and excited states of BeH2 and LiH molecules, and the density of states, which agrees well with exact solutions. Additionally, hardware demonstration is presented using IBM quantum devices for H2 molecule.

I. INTRODUCTION
Quantum computers offer the ability to address problems in quantum many-body chemistry and physics by quantum simulation or in a hybrid quantum-classical approach. The latter method is considered the most promising approach for noisy-intermediate scale quantum (NISQ) devices [1]. The prospect and benefits of quantum algorithms, along with suitable hardware, is in overcoming the complexity of the wave-function of a quantum system as it scales exponentially with system size [2]. Therefore developing techniques and algorithms for NISQ era devices that may prove to have some computational advantage themselves, or establish a path towards ideas and foundations that provide advantage for future error-corrected quantum devices, is a worthwhile pursuit.
The leading algorithms intended to be executed on NISQ devices, which aim to determine solutions to an electronic Hamiltonian, are variational in nature [3]. One specific algorithm is the variational quantum eigensolver (VQE), which has been tremendously successful in addressing chemistry and physics problems on quantum hardware and NISQ devices [4][5][6][7][8][9][10][11][12]. However, the restriction or challenge that exist with VQE is the need for prior insight with regard to selecting the trial quantum state, i.e., ansatz circuit. Furthermore, the classical optimization of the ansatz parameters may be a poorly converging problem [13,14] and therefore limiting the applicability of VQE for obtaining results accurate enough for chemical or physical interpretation. Finally, the realization of ansatz circuits that are motivated by domain knowledge, for example the unitary coupled cluster ansatz for chemistry problems [15], may not be directly applicable on NISQ hardware and therefore requires clever modification to obtain hardware efficient ansätze [6,7,16].
In this work, we present a pragmatic hybrid quantum- * Email: jouzdanip@fusion.gat.com; corresponding author classical approach for calculating the eigenenergy spectrum of a quantum system within an effective model. Firstly, an effective Hamiltonian is obtained through measurement of short-depth quantum circuits. The effective Hamiltonian is essentially the projection of the quantum system Hamiltonian onto a limited set of computational bases. The basis set is prepared with the intent of ensuring the dimensions of the corresponding matrix does not grow exponentially with the system size.
In order to evaluate the matrix elements of the effective Hamiltonian, suitable non-parametric quantum circuits are specified. The quantum circuits are designed, executed, and measured. From the result of the measurements, the diagonal and off-diagonal terms of the effective Hamiltonian matrix are obtained. On the classical side, the effective Hamiltonian matrix, with suitable dimensions, is diagonalized numerically using a classical computer. The paper is organized as follows: A short background is presented in Section II. In Section III the steps taken in our hybrid quantum-classical approach are explained in detail. In Section IV we demonstrate the application of this hybrid approach on simple chemical molecules BeH 2 and LiH. In Section V the approach is demonstrated on the IBMQ 5-qubit Valencia quantum processor [17] for H 2 molecule. Finally, we discuss the integration of VQE and the proposed approach in Section VI.

II. BACKGROUND
Consider a quantum many-body system of electrons with the second quantized Hamiltonian: a † i and a i are the creation and annihilation operators, respectively. The anticommutator for the creation and annihilation are given by: a i a † j +a † j a i = δ ij and a i a j +a j a i = a † i a † j + a † j a † i = 0. These rules enforce the non-abelian group statistics for fermions, that is, under exchange of two fermions the wave-function yields a minus sign.
The indices in Eq. (1) refer to single-electron states. The coefficients κ ij and v ijkl are the matrix integrals and v ijkl = ij|V 12 |kl , whereK 1 andV 12 operators correspond to one-and twobody interactions respectively. SinceK 1 andV 12 can depend on other parameters, such as the distance between nuclei, the Hamiltonian in Eq. (1) represents a class of problems. However, this class of problems has the common property that the number of fermions is a conserved value. Strictly speaking, the terms in the Hamiltonian act on fixed-particle-number Hilbert spaces, H N F , that have the correct fermionic antisymmetry, with N F denoting the number of electrons. In this paper, we consider this class of electronic systems where the Hamiltonian is assumed to be in the form of Eq. (1). The coefficients expressed in Eqs. (2)(3) can be obtained using software packages developed for quantum chemistry calculations that perform efficient numerical integration [18].

A. Mapping to qubits & computational basis
The Hamiltonian as written in Eq. (1) can be expressed in the form of qubit operations (i.e., Pauli matrices). This requires a transformation that preserves the anticommutation of the annihilation and creation operators. One transformation that satisfies the criteria, and is based on the physics of spin-lattice models, is the Jordan-Wigner (JW) transformation [19]. The JW-transformed Hamiltonian takes the form in which λ s 's are scalar and a Pauli string operatorĥ s is defined asĥ O s i ∈ {Î,X,Ŷ ,Ẑ} acts on the i-th qubit, {X,Ŷ ,Ẑ} are the three Pauli matrices [20] andÎ is the identity matrix with the number of qubits denoted as N .
If the number ofÎ operators in the tensor product of h s is N − k, we callĥ s a k-local Pauli string operator. Upon the JW transformation of the Hamiltonian, a Fock basis of the second quantization representation is in oneto-one correspondence with a computational basis of the qubits [21]. In other words, a computational basis of |n = |n 0 , n 1 , . . . , n N , with N qubits, where n i ∈ {0, 1}, is equivalent to a an antisymmetric Fock basis. Within the finite, but exponentially large, Hilbert space spanned by 2 N computational basis set, an effective matrix representation of the Hamiltonian may be possible, specifically, if one can efficiently evaluate the matrix elements n |Ĥ |n for an arbitrary computational basis |n and |n . Furthermore, assuming that the dimensions of the resulting effective matrix are relatively small, the matrix can be diagonalized on a classical computer, where its eigenvalues approximate the spectrum of the original Hamiltonian.
In this paper, we show how to evaluate a matrix element n |Ĥ |n for arbitrary computational basis |n and |n , using a quantum circuit that has a circuit depth O(N ). We do so by using ancilla qubits, and thus N + 1 physical resources are needed. In addition, we discuss how to choose an effective subspace for a given electronic Hamiltonian, with a dimension N s , based on physical motivations (see Section III C). The condition N s << 2 N makes it possible to diagonalize the Hamiltonian on a classical computer. In Section IV, we numerically demonstrate this method for simple quantum chemistry systems, focusing on ground state energy and the density of state calculations of the low-energy spectrum. Furthermore, in Section V we use the IBMQ 5-qubit Valencia device to measure the terms for H 2 in the complete computational basis of 4 qubits.

A. Effective Hamiltonian and circuit representation
We first consider a HamiltonianĤ that is expressed in terms of Pauli strings as in Eq. (4). Additionally, a subspace S = {|n }, with N s corresponding computational bases, is considered such that N s << 2 N . Let us define the effective Hamiltonian matrix as the projection ofĤ onto this subspace; that iŝ The next step is to define a simple quantum circuit that utilizes one ancillary qubit to measure n|Ĥ |n matrix element.
The dimension ofĤ ef f depends on the choice of the subspace in S. The choice of the subspace, intuitively, depends on the physics of the problem. However, the focus of this paper is towards quantum chemistry problems, which are a primary application for NISQ devices. For this class of Hamiltonians there is a systematic way to select the appropriate subspace. This is discussed in Section III C.
The evaluation of diagonal terms inĤ ef f , e.g., n|Ĥ |n , is trivially performed by preparing N qubits The off-diagonal matrix elements, e.g., n |Ĥ |n , which are generally complex numbers, can be evaluated using a single ancillary qubit. This can be done by considering the quantum circuits as shown in Fig. 1. The two circuits shown in Fig. 1(a) and (b) are used to calculate the real and imaginary parts of the matrix element, respectively. In both circuits the N + 1 qubits are initially prepared in |ψ int = |0 ⊗ |0 ⊗N state. The qubits are assumed to be enumerated linearly from 1 to N + 1, where the N + 1-th qubit is the control qubit.
After applying the first Hadamard gate on the ancillary qubit, from left to right as shown in the circuits in Fig. 1(a-b), the quantum state of all the qubits is which, after a sequence of controlled-X gates, becomes entangled as Here, control gates (controlled-n and controlled-n ) flip N F qubits (corresponding to the N F occupied electronic states) and prepares the target qubits in the computational basis |n (|n ), conditioned on the state of the control qubit is |1 (|0 ). An example of a controlled-n gate that prepares target qubits in |011 state is shown in Fig. 1(c). In practice, this part of the circuit requires two-qubit gates (e.g., CNOT) and perhaps the need for full connectivity of qubits in order to operate on any two qubits. Full connectivity of qubits could potentially be realized with iontrapped devices [22]. See Section VI for further discussion.
Depending on whether the real part or the imaginary part is calculated, the last gates acting on the control qubit changes. With regard to the real part, after applying the last Hadamard gate on the control qubit in Fig. 1(a), the quantum state is Using this prepared quantum state, one can measurê at the end of the circuit and have Thus, after substituting for the diagonal elements n|Ĥ |n and n |Ĥ |n , using Eq. (8), one obtains the real part of the off-diagonal matrix element n|Ĥ |n . The value m 0 in Eq. (13) is measured on a quantum device using the identity |0 0| = 1 2 (Î +Ẑ), by The imaginary part can be obtained in a similar fashion as done for the real part, but with a slight modification to the circuit as shown in Fig. 1(b). The key addition is a phase-gate, S, before the execution of the last Hadamard gate on the control qubit. This yields the quantum state that now includes a phase factor i before |n . After the last Hadamard gate in Fig. 1(b), and following the same steps, Eqs. (12)- (14), the imaginary part Im n|Ĥ |n ) is obtained. Our approach differs from the typical hybrid quantumclassical paradigm used in ground state chemistry electronic structure calculations in that the quantum hardware is used as a coprocessor to measure these matrix elements. Therefore, no parameterized ansatz or variational optimization is required. In this approach, the depth of the quantum circuit is significantly reduced, however, our method is based on the assumption that the dimensions ofĤ ef f in Eq. (7) are reasonable enough such that it can be diagonalized using classical numerical algorithms.

B. Implementing Measurements
As shown in Eq. (14), the expectation value of the Hamiltonian becomes the weighted sum of expectations for the set of Pauli string operators with respect to the output quantum state of the circuit. Since these operators are not in general commuting, one needs to setup and run a number of different quantum circuits. This number can be up-to all Pauli string operators in the Hamiltonian. Each circuit is then executed many times and every time the qubits are measured, the results are realized in different computational bases. The sampled realization provides a probability distribution and is used to estimate the expectation value of the Pauli string.
Generally, measurements are done either directly or indirectly. In the case of direct measurements, singlequbit rotations are applied to a subset of qubits at the end of the circuit. This subset is identified based on the locations in the tensor-product of the operatorĥ s that are not identityÎ. This set of rotations essentially changes the computational bases in which the given operatorĥ s is diagonal. The direct measurement is commonly used in the experimental demonstration of quantum hardware and VQE [6].
The indirect measurement approach [23] requires a series of controlled gates that are applied to the N target qubits, using one ancillary control qubit. The indirect measurement method is used in iterative quantum phase estimation algorithms [24].
Although the two types of measurement approaches are theoretically equivalent [20], experimentally, there are differences. The benefits and drawbacks of direct and indirect measurements are discussed for example in [25]. The main difference is the number of times the circuit is to be executed to achieve a desired precision , which is O( 1 2 ) and O( 1 ) for direct and indirect measurements respectively. The implementation of a general control-U gate in the indirect measurement is a challenging task [23]. However, in regards to our purposed approach, where the U operator is a single Pauli string, the indirect measurement implementation is straightforward.
The indirect measurement can be adapted as shown in quantum circuit in Fig. 2. An additional control qubit is added to this circuit. The Pauli operators are applied on the k locations ofĥ s when the control is in state |1 . It is straightforward to follow the quantum state of the qubits throughout the circuit. The final state reads as: Upon measuringM = |0 0| ⊗ |0 0| ⊗Î, at the end of this circuit, using Eq. (16), and steps discussed in sec. III A, the real part of the off-diagonal matrix element n |ĥ s |n is obtained. The imaginary part is determined under same steps, while an additional phase gate is applied to the last ancillary qubit, similar to the situation in Fig. 1(b).
In terms of Pauli strings, the value of ψ|M |ψ is obtained from the expectations ψ|Ẑ N +2 |ψ , ψ|Ẑ N +1 |ψ , and ψ| Ẑ N +2 ⊗Ẑ N +1 |ψ , similar to Eq. (14). Note that the approach for evaluating the real and imaginary parts of a matrix element is similar to the interference method introduced in ref. [10], with the difference being that our approach uses an ancillary qubit in order to realize the interference.

C. Preparing the computational basis
The computational basis set S = {|n } with N s << 2 N , needs to be specified in practice. These bases serve as the row and columns ofĤ ef f in Eq. (7). The process and motivation for how to choose this set should be based on the underlying nature and physics of the problem.
In theory, one established approach to approximate the ground states of quantum many-body systems is meanfield theory [26]. In this approach, the true ground state is constructed by perturbing a reference mean-field quantum state. The quantum chemistry field has established theories and techniques for treating such problems. One particularly successful theory and numerical method is coupled cluster (CC), typically referred to as the goldstandard in computational quantum chemistry [27]. In CC one assumes a wave-function ansatz for the ground state. Here |0 is a reference quantum state (e.g., Hartree-Fock) and is considered to be antisymmetric under exchange of two fermions. The oper-atorT =T 1 +T 2 + · · · is a sum over different possible excitation operators with respect to the reference state. Typically, the set of excitation operators inT includes single and double terms (i.e., CCSD), which enables a series representation of the Taylor expansion of eT , but high-order terms can also be added. The coefficients for the excitation operators insideT are determined by variational methods; that is, by minimizing the expectation of the Hamiltonian with respect to the ansatz [27].
Since the exponential operator in the CC ansatz, Eq. (17), is non-unitary, it cannot be directly implemented on gate-based quantum computers, where gates correspond to unitary operators. Thus, a unitary version of the CC ansatz has been introduced and is known as Unitary Coupled Cluster (UCC) [15,28]. Ideally, the implementation of UCC ansatz should be constructed such that the number of gates is minimized so that the circuit depth does not exhaust the current coherence times of NISQ devices. As a result of the this concern, hardware efficient anstaz have been proposed [6] as a substitute.
In the numerical demonstration of this work (see Section IV), we consider a simplified ansatz for the ground state as: where i, j, . . . refer to the unoccupied levels, and ν, β, . . . , refer to occupied levels with respect to single (S ), double (D), and higher order excitation operators. The ansatz in Eq. (18) implies that the true many-body ground state is a superposition of the reference state |0 and all possible S {a † i a ν |0 }, D {a † i a † j a ν a β |0 }, up-toT n |0 excitations, where n inT n is finite and independent of N .
In particular, assuming that it is possible to truncate the series at some low-excitation level such as the D or triples (T ), the number of eigenstates in the expansion of the above ansatz remains a polynomial function in N .
Taking the ansatz in Eq. (18), we specify the set S in the following way. (1) Pick a computational basis as the reference quantum state, that is |0 = |n int . (2) Identify computational bases corresponding to a finite number of excitations such as S and D. These are N s << 2 N bases and polynomial function in N .
In step (1), we identify the initial computational basis |n int by minimizing the n int |Ĥ |n int , in which, for example, a classical Monte Carlo process from spin lattice models [29] can be used. This computational basis is essentially the qubits' configuration that has the lowest energy expectation. Since this step is a classical one, it is performed effectively even for a large number of qubits.
In step (2), once the state |n int is determined, one can rearrange the configuration of the qubits by swapping 1's and 0's within the state |n int . The swapping is done so that the configurations corresponding to S, D, up-toT n are fully realized. The energy expectation corresponding to these configurations (the diagonal elements in the H ef f ) are stored on a classical register. The final result is obtained among the set of configurations whereby N s of are the lowest energies; these are the configurations that are selected. The above steps are demonstrated numerically which is discussed in the next section (see Section IV).

IV. NUMERICAL DEMONSTRATION: LIH AND BEH2
The application of the methodology discussed in Section III is focused on BeH 2 and LiH molecules due to thier relatively small number of electrons and molecular orbital footprint. The number of electrons in BeH 2 and LiH is six and four, respectively. The single-electron molecular spin-orbitals in the second quantized Hamiltonian are constructed using the minimal atomic STO-3G basis set [2,6]. For BeH 2 there are a total of 14 spinorbital, and thus corresponds to 14 qubits; LiH contains 12 spin orbitals and hence 12 qubits. For the purpose of measuring the matrix elements of the effective Hamiltonian, Eq. (7), one additional ancillary qubit is required, as illustrated in the quantum circuits shown in Fig. 1. Therefore, the total number of qubits for the simulation of BeH 2 and LiH is 15 and 13, respectively.
To obtain the coefficients in the Hamiltonian of Eq. (1) -more specifically as defined in Eqs. (2) and (3) we make use of the Psi4 quantum chemistry package [30]. The second quantized Hamiltonian is further transfomed onto a set of Pauli strings and their corresponding weights, Eq. (1), via JW transform using OpenFermion package [31].
In order to construct the potential energy surface for each inter-nuclear distance, R, the steps indicated in the previous paragraph are repeated. The distance R corresponds to the bond length between Be-H or Li-H in a given molecule; both LiH and BeH 2 are linear molecules. We note that these calculations are assuming the total Hamiltonian can be represented using the Born-Oppenheimer approximation, where the dynamics of the core nuclei are neglected. This is standard prac-tice in quantum chemistry calculations [26] and thus at every distance R, the nuclear-nuclear repulsion energy is treated classically and is added to the Hamiltonian as a constant.
For every distance R, a set of computational bases is chosen. This set includes the basis with the lowest energy expectation (the reference configuration) and computational bases that correspond to the low-order excitations (i.e., S, D, etc.) with respect to the reference configuration. In the case of BeH 2 , we construct the effective Hamiltonian matrix by keeping the S, D, and T excitations, and thus the total number of computational bases for BeH 2 is N s = 1588. The low dimension of the effective Hamiltonian matrix makes it possible to diagonalize the matrix on a classical computer using standard numerical techniques [32]. Thus, the lowest energies, including the ground state energy, are obtained at every R. The results for this process are shown in Fig. 3. Fig. 3(a) shows the ground state energy as well as a few low-lying excited states for BeH 2 of the effective Hamiltonian. The exact energies are given by the dashed curves in Fig. 3(a). In Fig. 3(b), the difference between the exact and obtained ground state energy is shown. An error within the hatched area indicates results that are within chemical accuracy. Chemical accuracy is typically identified as ∼ 5 × 10 −3 Hatrees. Figure 3(c) shows the same energy difference for all other excited-state energies. The same process is done for LiH where only S and D configurations are used. The total number of permutations of 4 fermions (495) is reduced to N s = 200, and the results are shown in Fig. 4.
The total number of configurations can be approximated by the highest excitation considered.
For n-electron excitation the number of configurations is Thus, the total number of configurations is less than 1 + · · · + N 2n = O(N 2n ) when N 1. Of course, this argument is valid as long as n can be assumed to be finite and independent of the size of the system N or the number of electrons N F . Under these assumptions, the dimensions of the effective Hamiltonian, H ef f , remains polynomial if one replaces the minimal STO-3G with the extended atomic basis.

A. Density of States
Knowledge of density of states (DOS) can be important in analyzing the thermodynamic behaviour of a system at finite temperature and in the analysis of transition states important in chemical reactions [33]. One advantage of the method proposed in this paper is the insights provided into the low-energy spectrum of the quantum system, more precisely, the degeneracy of the energy levels. To illustrate this, in Fig. 5, the unnormalized DOS of the LiH obtained via the effective Hamiltonian is compared with the exact density, which shows qualitatively good agreement. In this section, we demonstrate the feasibility of the hybrid quantum-classical approach discussed in this paper on a IBMQ hardware device. Specifically, we calculate the eigenenergy spectrum of H 2 molecule using STO-3G minimal basis-set.
Within STO-3G basis set, the total number of spinorbitals for H 2 is four. In contrast to the numerical demonstration, where the computational basis was prepared using the UCC ansatz, here we employ the complete computational basis with conserved number of electrons for H 2 (N F = 2). The number of computational basis set with conserved electron number, is N s = 6, and the bases are: The number of bases is not further reduced; that is, the effective Hamiltonian constructed here is equal to the full Hamiltonian within the N F = 2 electron number. The above set can be interpreted as the collection of all the S and D exciations, as introduced in sub-section III C, above a classical configuration such as |1100 . The Hamiltonian of the problem, stated as a weighted sum of Pauli string operators, as in Eq. (4), has a total of fifteen operators, which can be grouped as: ZIII, IZII, IIZI, IIIZ, ZZII,  ZIZI, ZIIZ, IZZI, IZIZ, IIZZ}  (20) and Here, a Pauli string operator such as ZIII is a shorthand forẐ 1 ⊗Î 2 ⊗Î 3 ⊗Î 4 , as introduced in Eq. (5). Our goal is to use the quantum hardware to evaluate the values for n|ĥ s |n , whereĥ s belongs to G 1 or G 2 , and |n and |n can be any of the bases listed in Eq. (19).
Note that this task is not a difficult one for a classical computer. However, the objective of the paper and hardware demonstration is to establish a hybrid quantumclassical computational scheme in which the quantum processor performs some of the computational tasks. In the simulations of this paper, the task is the evaluation of the phase associated with the exchange of fermions; since n|ĥ s |n can be understood as exchanging fermions on |n via the operatorĥ s , and to arrive at |n . The output can only take discrete values of ±1, ±i, and 0. In each group G 1 and G 2 , Pauli string operators commute and thus, theoretically, can be measured simultaneously. However, in order to perform simultaneous measurement, one needs to adjust the computational bases, by applying further quantum gates on the target qubits. The appended measurement circuit further increases the circuit depth [34]. Thus, practically, the simultaneous measurement may not be advantageous, and, in contrast, induces further noise and decoherence on NISQ devices. Nevertheless, for G 1 , no further circuit depth is required, as they are all measured in the z-basis.
We can further reduce the quantum coprocessor computational load by excluding circuits containing terms that are classically efficient to compute. In eqs. (20)(21), the operators in G 1 do not contribute to off-diagonal matrix element of the Hamiltonian, while operators in G 2 contribute only to the off-diagonal matrix elements. This information is used to reduce the essential number of circuits to be run on the hardware.
The evaluation of an element n|ĥ s |n is performed by assigning an appropriate circuit, as introduced in subsection III A, to a quantum coprocessor. In this work we make use of the IBMQ cloud open devices.
Once a circuit is loaded onto a quantum coprocessor, each circuit is executed a several times. The number of times a circuit is run is referred to as the number of shots. In this work, each circuit was executed for a total of 8, 000 shots. At each run, the qubits are measured, and the collection of the realized configurations is used to evaluate value of n|ĥ s |n .
The preparation of the circuits is performed using the transpile and assemble routines available in the IBM Qiskit API [35], which enables collating individual circuits, corresponding to the diagonal and off-diagonal terms, to prepare each circuit to run as a single job on a target IBMQ device. In this work, we make use of the IBMQ 5-qubit Valencia hardware [17]. The number of circuits -for both real and imaginary parts separately -after the transpile and assemble routines for the diagonal and off-diagonal terms corresponds to 66 and 60 circuits, respectively. The execution timings for the IBMQ QASM simulator range between 0.21-1.5 seconds and for the IBMQ Valencia 5 qubit device 274-298 seconds with 8000 shots per circuit. Due to the limited number of qubits available on this processor, direct measurement is used, that is by measuring all the qubits at the end of the circuit.
As a sanity check, we first use IBMQ QASM simulator [36] to perform all steps involved as outlined in III. The calculations are shown in Fig. 6(a). The results perfectly match the exact values as expected. Furthermore, in Fig. 6(b) the off-diagonal terms in the Hamiltonian are evaluated using IBMQ Valencia device, while the diagonal elements are evaluated using the IBMQ QASM simulator. In both cases, upon diagonalization of the obtained Hamiltonian, the energy spectrum is in good agreement with the exact values. We observe in Fig. 6(b) that there is a slight lifting of the double degeneracy of H 2 , which we associate to the device noise. This is due to the occurance of non-zero values in the measured off-diagonal terms when in fact the exact value should be zero.
Finally, all the matrix elements, that is, diagonal and off-diagonal elements, are evaluated using IBMQ Valencia device. The results are shown in Fig. 7-(a), which was done with and without measurement error mitigation. The use of measurement error mitigation marginally improves the results. The discrepancy between the exact and hardware values, without and with measurement error mitigation, are shown in 7(b-c).
The effect of noise is ubiquitous in current quantum hardware and is inherently complex and difficult to characterize. However, the result of noise due to the measurement process can be mitigated to some degree by using a calibration matrix [35]. The resulting calibration matrix is then used to reduce errors of subsequent circuit measurements. The application of measurement error mitigation is demonstrated in Fig. 7(a) with beneficial impact on error shown in Fig. 7(c).
The calculations in Fig. 7, that is construction of the effective Hamiltonian,Ĥ ef f , are performed through five independent IBMQ job submissions. Each time the energysurface is slightly different, which can be associated to the inherent device errors and perhaps low number of shots used. The error bars in Fig. 7(b-c) indicate the min/max range in errors of the five different calculation.

VI. DISCUSSION AND SUMMARY
Recently an abundance of Variational Quantum Algorithms (VQA) have been introduced to solve electronic quantum many-body problems on NISQ devices [3-5, 8-11, 28]. These VQA proposals mostly rely on optimization of parametric circuits. In this paper we demonstrate that an alternative approach exists, specifically for quantum chemistry problems, which does not require an optimization procedure and parametric ansatz. In this approach, the quantum hardware is used only to prepare a many-body quantum state and to efficiently measure the expectation values of certain target observables with respect to this quantum state. From the output measurements, one can construct what we refer to as an effective Hamiltonian matrix,Ĥ ef f . Upon diagonalization ofĤ ef f using classical eigen-decomposition numerical methods one obtains the ground state and low-lying excited states of the system.
The approach introduced in this paper is particularly suitable for NISQ devices where short quantum circuit depth is essential due to lack of error-correction protocols on these devices. An additional important aspect of this approach is that it provides access to the low-lying energy In the context of VQE and VQE-type algorithms, several attempts have been made to extend the variational approach to excited states [10,11,[37][38][39][40]. Quantum subspace expansion [41,42] for example constructs a set of non-orthogonal bases out of an optimized ansatz, and performs post-processing to obtain excited states. Deflation techniques as described in ref. [11], constructs a pseudo-Hamiltonian in which the ground state is excluded and orthogonality is enforced through regularization. Successful examples are introduced for some lowlying excited states of LiH [8]. In all these previous works, optimization of a parametric ansatz is required therefore necessitates an enormous number of quantum circuit executions and sampling.
The approach introduced in this paper is similar to multistate contracted VQE (MC-VQE) in ref. [10]. The main difference is the application of VQE: MC-VQE is obviously a VQE-type algorithm, our approach is dis-tinct in that it does not require a parametric circuit or variational procedure for optimization. In addition, our method differs since it uses supporting quantum circuit resources, i.e., ancilla qubits, in order to perform interference and measurement that is different from the quantum circuit in ref. [43] and its generalization in ref. [10].
An extension of our approach to VQE type algorithm is possible. This can be done by appending a set of parametric gates that act only on the target qubits to the circuits in Fig. 1. Let us denote this part of the circuit with U (θ), where θ stands for a set of parameters. Then, it can be verified that, given θ, the final matrix element obtained from the circuit after measurement (see Eq. (2) for example) becomes n| U † (θ)ĤU (θ) |n , compared to n|Ĥ |n . The appended parametric circuit U (θ) allows one to project the Hamiltonian onto S(θ) = {U (θ) |n }, for a given θ. This means the N s × N s effective Hamiltonian is now parametric and depends on the value(s) of θ. The optimal parameter(s) are then obtained by minimizing the ground-state energy of the effective Hamiltonian matrix. The reference state {U (θ) |n } can be regarded as the contracted reference states introduced in ref. [10].
One possible limitation of the circuits shown in Fig. 1 is the execution of two-qubit gates corresponding to control operations over n-qubits (i.e., series of CNOT gates). For NISQ devices with hardware-restricted qubit connectivity, this may require a number of SWAP gate operations and therefore can increase the circuit depth and subsequent error rates [44,45]. In essence, the implementation of the circuits described in this paper will depend on the ability to limit circuit depth and associated error rates by NISQ hardware circuit optimization (i.e., scheduling). However, significant improvements in qubit connectivity of various modalities (e.g., ion traps) [22,46] or optimizing quantum circuits against decoherence [47,48] may blunt this concern.
VII. ACKNOWLEDGEMENT The ideas and methodology discussed in Section III have been supported by General Atomics internal R&D funds. Their application towards quantum chemistry problems as discussed in Section IV and Section V is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences, under Award Number DE-SC0020249. We thank Mark Kostuk for his guidance and management as PI under this grant. Additionally, we thank David P. Schissel for feedback and guidance at General Atomics.
Circuit diagrams are rendered using the L A T E X Quantikz package [49], numerical scripts utilize the Python SciPy package [32], and 2D plots are generated with the Python matplotlib package [50]. We acknowledge use of the IBMQ for this work. The views expressed are those of the authors and do not reflect the official policy or position of IBM or the IBMQ team.
Disclaimer: A portion of this report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Refer-ence herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.