Variational algorithms for strongly correlated chemical and materials systems are one of the most promising applications of near-term quantum computers. We present an extension to the variational quantum eigensolver that approximates the ground state of a system by solving a generalized eigenvalue problem in a subspace spanned by a collection of parametrized quantum states. This allows for the systematic improvement of a logical wavefunction ansatz without a significant increase in circuit complexity. To minimize the circuit complexity of this approach, we propose a strategy for efficiently measuring the Hamiltonian and overlap matrix elements between states parametrized by circuits that commute with the total particle number operator. This strategy doubles the size of the state preparation circuits but not their depth, while adding a small number of additional two-qubit gates relative to standard variational quantum eigensolver. We also propose a classical Monte Carlo scheme to estimate the uncertainty in the ground state energy caused by a finite number of measurements of the matrix elements. We explain how this Monte Carlo procedure can be extended to adaptively schedule the required measurements, reducing the number of circuit executions necessary for a given accuracy. We apply these ideas to two model strongly correlated systems, a square configuration of H4 and the π-system of hexatriene (C6H8).
Variational algorithms for strongly correlated chemical and materials systems are one of the most promising applications of near-term quantum computers. We present an extension to the variational quantum eigensolver that approximates the ground state of a system by solving a generalized eigenvalue problem in a subspace spanned by a collection of parametrized quantum states. This allows for the systematic improvement of a logical wavefunction ansatz without a significant increase in circuit complexity. To minimize the circuit complexity of this approach, we propose a strategy for efficiently measuring the Hamiltonian and overlap matrix elements between states parametrized by circuits that commute with the total particle number operator. This strategy doubles the size of the state preparation circuits but not their depth, while adding a small number of additional two-qubit gates relative to standard variational quantum eigensolver. We also propose a classical Monte Carlo scheme to estimate the uncertainty in the ground state energy caused by a finite number of measurements of the matrix elements. We explain how this Monte Carlo procedure can be extended to adaptively schedule the required measurements, reducing the number of circuit executions necessary for a given accuracy. We apply these ideas to two model strongly correlated systems, a square configuration of H4 and the π-system of hexatriene (C6H8).
Variational algorithms for strongly correlated chemical and materials systems are one of the most promising applications of near-term quantum computers. We present an extension to the variational quantum eigensolver that approximates the ground state of a system by solving a generalized eigenvalue problem in a subspace spanned by a collection of parametrized quantum states. This allows for the systematic improvement of a logical wavefunction ansatz without a significant increase in circuit complexity. To minimize the circuit complexity of this approach, we propose a strategy for efficiently measuring the Hamiltonian and overlap matrix elements between states parametrized by circuits that commute with the total particle number operator. This strategy doubles the size of the state preparation circuits but not their depth, while adding a small number of additional two-qubit gates relative to standard variational quantum eigensolver. We also propose a classical Monte Carlo scheme to estimate the uncertainty in the ground state energy caused by a finite number of measurements of the matrix elements. We explain how this Monte Carlo procedure can be extended to adaptively schedule the required measurements, reducing the number of circuit executions necessary for a given accuracy. We apply these ideas to two model strongly correlated systems, a square configuration of H 4 and the π -system of hexatriene (C 6 H 8 ).
Quantum simulation of chemical systems is one of the most promising near-term applications of quantum computers. The variational quantum eigensolver, a leading algorithm for molecular simulations on quantum hardware, has a serious limitation in that it typically relies on a pre-selected wavefunction ansatz that results in approximate wavefunctions and energies. Here we present an arbitrarily accurate variational algorithm that, instead of fixing an ansatz upfront, grows it systematically one operator at a time in a way dictated by the molecule being simulated. This generates an ansatz with a small number of parameters, leading to shallow-depth circuits. We present numerical simulations, including for a prototypical strongly correlated molecule, which show that our algorithm performs much better than a unitary coupled cluster approach, in terms of both circuit depth and chemical accuracy. Our results highlight the potential of our adaptive algorithm for exact simulations with present-day and near-term quantum hardware.
We propose a modification of the Variational Quantum Eigensolver algorithm for electronic structure optimization using quantum computers, named nonunitary Variational Quantum Eigensolver (nu-VQE), in which a nonunitary operator is combined with the original system Hamiltonian leading to a new variational problem with a simplified wave function ansatz. In the present work, as nonunitary operator, we use the Jastrow factor, inspired from classical Quantum Monte Carlo techniques for simulation of strongly correlated electrons. The method is applied to prototypical molecular Hamiltonians for which we obtain accurate ground-state energies with shallower circuits, at the cost of an increased number of measurements. Finally, we also show that this method achieves an important error mitigation effect that drastically improves the quality of the results for VQE optimizations on today's noisy quantum computers. The absolute error in the calculated energy within our scheme is 1 order of magnitude smaller than the corresponding result using traditional VQE methods, with the same circuit depth.
We describe the contextual subspace variational quantum eigensolver (CS-VQE), a hybrid quantum-classical algorithm for approximating the ground state energy of a Hamiltonian. The approximation to the ground state energy is obtained as the sum of two contributions. The first contribution comes from a noncontextual approximation to the Hamiltonian, and is computed classically. The second contribution is obtained by using the variational quantum eigensolver (VQE) technique to compute a contextual correction on a quantum processor. In general the VQE computation of the contextual correction uses fewer qubits and measurements than the VQE computation of the original problem. Varying the number of qubits used for the contextual correction adjusts the quality of the approximation. We simulate CS-VQE on tapered Hamiltonians for small molecules, and find that the number of qubits required to reach chemical accuracy can be reduced by more than a factor of two. The number of terms required to compute the contextual correction can be reduced by more than a factor of ten, without the use of other measurement reduction schemes. This indicates that CS-VQE is a promising approach for eigenvalue computations on noisy intermediate-scale quantum devices.
The Non-Orthogonal Quantum Eigensolver (NOQE) provides an accurate framework for electronic-structure calculations, but the estimation of its Hamiltonian and overlap matrix elements relies on sampling and requires O(1/ε2) circuit repetitions to achieve additive precision ε. Here, we reformulate this matrix-element estimation step as a collection of amplitude-estimation tasks and integrate iterative quantum amplitude estimation into the NOQE workflow. The resulting protocol achieves near-Heisenberg query complexity O(1/ε) for these estimation tasks, by replacing incoherent statistical averaging with coherent amplitude amplification. We present explicit circuit constructions and the corresponding implementation procedure. Numerical simulations for the electronic states of the hydrogen molecule show that the proposed method reaches chemical accuracy with substantially fewer total queries than the original sampling-based protocol. Overall, this work provides a measurement-efficient route to high-precision energy estimation and illustrates how sampling-limited quantum algorithms can be systematically reformulated to leverage quantum coherence and achieve lower measurement costs.
The variational quantum eigensolver (VQE) is a method of choice to solve the electronic structure problem for molecules on near-term gate-based quantum computers.However, the circuit depth is expected to grow significantly with the problem size.Increased depth can both degrade the accuracy of the results and reduce trainability.In this work, we propose an approach to reduce ansatz circuit depth.Our approach, called "PermVQE," adds an additional optimization loop to the VQE that permutes qubits in order to solve for the qubit Hamiltonian that maximally localizes correlations in the ground state.The choice of permutations is based on mutual information, which is a measure of interaction between electrons and/or holes in spin-orbitals.Encoding strongly entangled spin-orbitals into proximal qubits on a quantum chip naturally reduces the circuit depth needed to prepare the ground state.For representative molecular systems, LiH, H 2 , (H 2 ) 2 , H = 4 , H + 3 , and N 2 , we demonstrate that placing entangled qubits in close proximity leads to shallower depth circuits required to reach a given eigenvalue-eigenvector accuracy.The approach is designed for hardware-efficient ansatz of any qubit connectivity, and examples are demonstrated for linear and twodimensional grid architectures.The main ideas can also be applied to simulate molecules with other ansatz as well as variational quantum algorithms beyond the VQE.In particular, we demonstrate the beneficial effect of qubit permutations to build fermionic-adaptive derivative assembled pseudo-Trotter ansatz on a linear qubit connectivity architecture with nearly a twofold reduction of the number of controlled NOT gates.