The variational quantum eigensolver (VQE) and its variants, which is a method for finding eigenstates and eigenenergies of a given Hamiltonian, are appealing applications of near-term quantum computers. Although the eigenenergies are certainly important quantities which determines properties of a given system, their derivatives with respect to parameters of the system, such as positions of nuclei if we target a quantum chemistry problem, are also crucial to analyze the system. Here, we describe methods to evaluate analytical derivatives of the eigenenergy of a given Hamiltonian, including the excited state energy as well as the ground state energy, with respect to the system parameters in the framework of the VQE. We give explicit, low-depth quantum circuits which can measure essential quantities to evaluate energy derivatives, incorporating with proof-of-principle numerical simulations. This work extends the theory of the variational quantum eigensolver, by enabling it to measure more physical properties of a quantum system than before and to perform the geometry optimization of a molecule.

The variational quantum eigensolver (VQE), a variational algorithm to obtain an approximated ground state of a given Hamiltonian, is an appealing application of near-term quantum computers. The original work [Peruzzo et al.; Nat. Commun.; 5, 4213 (2014)] focused only on finding a ground state, whereas the excited states can also induce interesting phenomena in molecules and materials. Calculating excited states is, in general, a more difficult task than finding ground states for classical computers. To extend the framework to excited states, we here propose an algorithm, the subspace-search variational quantum eigensolver (SSVQE). This algorithm searches a low energy subspace by supplying orthogonal input states to the variational ansatz and relies on the unitarity of transformations to ensure the orthogonality of output states. The k-th excited state is obtained as the highest energy state in the low energy subspace. The proposed algorithm consists only of two parameter optimization procedures and does not employ any ancilla qubits. The disuse of the ancilla qubits is a great improvement from the existing proposals for excited states, which have utilized the swap test, making our proposal a truly near-term quantum algorithm. We further generalize the SSVQE to obtain all excited states up to the k-th by only single optimization procedure. From numerical simulations, we verify the proposed algorithms. This work greatly extends the applicable domain of the VQE to excited states and their related properties like a transition amplitude without sacrificing any feasibility of it.

The variational quantum eigensolver (VQE) is an attracting possible application of near-term quantum computers. Originally, the aim of the VQE is to find a ground state for a given specific Hamiltonian. It is achieved by minimizing the expectation value of the Hamiltonian with respect to an ansatz state by tuning parameters θ on a quantum circuit which constructs the ansatz. Here we consider an extended problem of the VQE, namely, our objective in this work is to “generalize” the optimized output of the VQE just like machine learning. We aim to find ground states for a given set of Hamiltonians {H(x)}, where x is a parameter which specifies the quantum system under consideration, such as geometries of atoms of a molecule. Our approach is to train the circuit on the small number of x’s. Specifically, we employ the interpolation of the optimal circuit parameter determined at different x’s, assuming that the circuit parameter θ has simple dependency on a hidden parameter x as θ(x). We show by numerical simulations that, using an ansatz which we call the Hamiltonian-alternating ansatz, the optimal circuit parameters can be interpolated to give near-optimal ground states in between the trained x’s. The proposed method can greatly reduce, on a rough estimation by a few orders of magnitude, the time required to obtain ground states for different Hamiltonians by the VQE. Once generalized, the ansatz circuit can predict the ground state without optimizing the circuit parameter θ in a certain range of x.