Quantum subspace diagonalization (QSD) methods are quantum-classical hybrid methods, commonly used to find ground and excited state energies by projecting the Hamiltonian to a smaller subspace. In applying these, the choice of subspace basis is critical from the perspectives of basis completeness and efficiency of implementation on quantum computers. In this work, we present Eigenvector Continuation (EC) as a QSD method, where low-energy states of the Hamiltonian at different points in parameter space are chosen as the subspace basis. This unique choice enables rapid evaluation of low-energy spectra, including ground and nearby excited states, with minimal hardware effort. As a particular advantage, EC is able to capture the spectrum across ground state crossovers corresponding to different symmetry sectors of the problem. We demonstrate this method for interacting spin models and molecules.
Quantum subspace diagonalization (QSD) methods are quantum-classical hybrid methods, commonly used to find ground and excited state energies by projecting the Hamiltonian to a smaller subspace. In applying these, the choice of subspace basis is critical from the perspectives of basis completeness and efficiency of implementation on quantum computers. In this work, we present Eigenvector Continuation (EC) as a QSD method, where low-energy states of the Hamiltonian at different points in parameter space are chosen as the subspace basis. This unique choice enables rapid evaluation of low-energy spectra, including ground and nearby excited states, with minimal hardware effort. As a particular advantage, EC is able to capture the spectrum across ground state crossovers corresponding to different symmetry sectors of the problem. We demonstrate this method for interacting spin models and molecules.
We propose quantum-selected configuration interaction, a class of hybrid quantum-classical algorithms for calculating the ground- and excited-state energies of many-electron Hamiltonians on noisy quantum devices. Suppose that an approximate ground state can be prepared on a quantum computer either by variational quantum eigensolver or by some other method. Then, by sampling the state in the computational basis, which is hard for classical computation in general, one can identify the electron configurations that are important for reproducing the ground state. The Hamiltonian in the subspace spanned by those important configurations is diagonalized on classical computers to output the ground-state energy and the corresponding eigenvector. The excited-state energies can be obtained similarly. The result is robust against statistical and physical errors because the noisy quantum devices are used only to define the subspace, and the resulting ground-state energy strictly satisfies the variational principle even in the presence of such errors. The expectation values of various other operators can also be estimated for obtained eigenstates with no additional quantum cost, since the explicit eigenvectors in the subspaces are known. We verified our proposal by numerical simulations and demonstrated it on a quantum device for an eight-qubit molecular Hamiltonian. The proposed algorithms are potentially feasible to tackle some challenging molecules by exploiting quantum devices with several tens of qubits, assisted by high-performance classical computing resources for diagonalization.
Abstract A typical task for classical and quantum computing in chemistry is finding a potential energy surface (PES) along a reaction coordinate, which involves solving the quantum chemistry problem for many points along the reaction path. Developing algorithms to accomplish this task on quantum computers has been an active area of development, yet finding all the relevant eigenstates along the reaction coordinate remains a difficult problem, and determining PESs is thus a costly proposal. In this paper, we demonstrate the use of a eigenvector continuation—a subspace expansion that uses a few eigenstates as a basis—as a tool for rapidly exploring PESs. We apply this to determining the binding PES or torsion PES for several molecules of varying complexity. In all cases, we show that the PES can be captured using relatively few basis states; suggesting that a significant amount of (quantum) computational effort can be saved by making use of already calculated ground states in this manner.
Eigenvector continuation (EC) has recently attracted a lot attention in nuclear structure and reactions as a variational resummation tool for many-body expansions. While previous applications focused on ground-state energies, excited states can be accessed on equal footing. This work is dedicated to a detailed understanding of the emergence of excited states from the eigenvector continuation approach. For numerical applications the one-dimensional quartic anharmonic oscillator is investigated, which represents a strongly non-perturbative quantum system where the use of standard perturbation techniques break down. We discuss how different choices for the construction of the EC manifold affect the quality of the EC resummation and investigate in detail the results from EC for excited states compared to results from a full diagonalization as a function of the basis-space size.
Subspace methods are powerful, noise-resilient methods that can effectively prepare ground states on quantum computers. The challenge is to get a subspace with a small condition number that spans the states of interest using minimal quantum resources. In this work, we will use eigenvector continuation to build a subspace from the low-lying states of a set of Hamiltonians. The basis vectors are prepared using truncated versions of standard state preparation methods such as imaginary-time evolution (ITE), adiabatic state preparation (ASP), and variational quantum eigensolver. By using these truncated methods combined with eigenvector continuation, we can directly improve upon them, obtaining more accurate ground-state energies at a reduced cost. We use several spin systems to demonstrate convergence even when methods like ITE and ASP fail, such as ASP in the presence of level crossings and ITE with vanishing energy gaps. We also showcase the noise resilience of this approach beyond the gains already made by having shallower quantum circuits. Our findings suggest that eigenvector continuation can be used to improve existing state preparation methods in the near term.
Subspace eigenvalue problems arise ubiquitously in quantum chemistry and condensed-matter physics, where the relevant object is often a low-energy manifold rather than a single ground-state wavefunction. In this work, we propose a fault-tolerant quantum algorithm for this subspace-level task based on the Feshbach effective-Hamiltonian formalism. Given block-encoding access to the full Hamiltonian and a chosen d-dimensional reference subspace, the algorithm estimates eigenvalues of states with nonzero overlap with the reference subspace through a local secant fixed-point search. It then implements the associated wave operator and prepares an orthonormal basis whose span approximates the target invariant subspace. The construction combines projected block encodings with quantum singular value transformation (QSVT), which approximates the complementary-space resolvent and thereby provides both the self-energy used for eigenvalue estimation and the wave operator used for eigenstate reconstruction. For target accuracy ε, a single evaluation of the effective Hamiltonian has query complexity O(d3/(g2ε)), up to block-encoding normalization factors, where g is the distance between the target eigenvalue and the nearest pole of the effective Hamiltonian. Under the stated local regularity conditions, the secant search requires only O(loglog(1/ε)) effective-Hamiltonian evaluations to reach the working precision. Classical numerical emulations for an open 4×2 Fermi--Hubbard cluster, all-electron LiH bond stretching, and [Ru(bpy)3]2+ demonstrate the resolution and reconstruction of low-energy states and manifolds across spin-sector crossings, near-degeneracies, and dense excited-state spectra.
Abstract Quantum subspace methods (QSMs) are a class of quantum computing algorithms where the time-independent Schrödinger equation for a quantum system is projected onto a subspace of the underlying Hilbert space. This projection transforms the Schrödinger equation into an eigenvalue problem determined by measurements carried out on a quantum device. The eigenvalue problem is then solved on a classical computer, yielding approximations to ground- and excited-state energies and wavefunctions. QSMs are examples of hybrid quantum–classical methods, where a quantum device supported by classical computational resources is employed to tackle a problem. QSMs are rapidly gaining traction as a strategy to simulate electronic wavefunctions on quantum computers, and thus their design, development, and application is a key research field at the interface between quantum computation and electronic structure (ES). In this review, we provide a self-contained introduction to QSMs, with emphasis on their application to the ES of molecules. We present the theoretical foundations and applications of QSMs, and we discuss their implementation on quantum hardware, illustrating the impact of noise on their performance.