Approximating the ground states of strongly interacting electron systems in quantum chemistry and condensed matter physics is expected to be one of the earliest applications of quantum computers. In this paper, we prepare highly accurate ground states of the Fermi-Hubbard model for small grids up to 6 sites (12 qubits) by using an interpretable, adaptive variational quantum eigensolver(VQE) called ADAPT-VQE. In contrast with non-adaptive VQE, this algorithm builds a system-specific ansatz by adding an optimal gate built from one-body or two-body fermionic operators at each step. We show this adaptive method outperforms the non-adaptive counterpart in terms of fewer variational parameters, short gate depth, and scaling with the system size. The fidelity and energy of the prepared state appear to improve asymptotically with ansatz depth. We also demonstrate the application of adaptive variational methods by preparing excited states and Green functions using a proposed ADAPT-SSVQE algorithm. Lower depth, asymptotic convergence, noise tolerance of a variational approach, and a highly controllable, system-specific ansatz make the adaptive variational methods particularly well-suited for NISQ devices.Approximating the ground states of strongly interacting electron systems in quantum chemistry and condensed matter physics is expected to be one of the earliest applications of quantum computers. In this paper, we prepare highly accurate ground states of the Fermi-Hubbard model for small grids up to 6 sites (12 qubits) by using an interpretable, adaptive variational quantum eigensolver(VQE) called ADAPT-VQE. In contrast with non-adaptive VQE, this algorithm builds a system-specific ansatz by adding an optimal gate built from one-body or two-body fermionic operators at each step. We show this adaptive method outperforms the non-adaptive counterpart in terms of fewer variational parameters, short gate depth, and scaling with the system size. The fidelity and energy of the prepared state appear to improve asymptotically with ansatz depth. We also demonstrate the application of adaptive variational methods by preparing excited states and Green functions using a proposed ADAPT-SSVQE algorithm. Lower depth, asymptotic convergence, noise tolerance of a variational approach, and a highly controllable, system-specific ansatz make the adaptive variational methods particularly well-suited for NISQ devices.

Sep 24, 20211 citation

Eigenstate solutions of the Fermi-Hubbard model via symmetry-enhanced variational quantum eigensolver

Yao, Shaohui, Wenyu Wang

The Variational Quantum Eigensolver (VQE), as a hybrid quantum-classical algorithm, is an important tool for effective quantum computing in the current noisy intermediate-scale quantum (NISQ) era. However, the traditional hardware-efficient ansatz without taking into account symmetries requires more computational resources to explore the unnecessary regions in the Hilbert space. The conventional Subspace-Search VQE (SSVQE) algorithm, which can calculate excited states, is also unable to effectively handle degenerate states since the loss function only contains the expectation value of the Hamiltonian. In this study, the energy eigenstates of the one-dimensional Fermi-Hubbard model with two lattice sites and the two-dimensional Hubbard model with four lattice sites are calculated. By incorporating symmetries into the quantum circuits and loss function, we find that both the ground state and excited state calculations are improved greatly compared to the case without symmetries. The enhancement in excited state calculations is particularly significant. This is because quantum circuits that conserve the particle number are used, and appropriate penalty terms are added to the loss function, enabling the optimization process to correctly identify degenerate states. The results are verified through repeated simulations.The Variational Quantum Eigensolver (VQE), as a hybrid quantum-classical algorithm, is an important tool for effective quantum computing in the current noisy intermediate-scale quantum (NISQ) era. However, the traditional hardware-efficient ansatz without taking into account symmetries requires more computational resources to explore the unnecessary regions in the Hilbert space. The conventional Subspace-Search VQE (SSVQE) algorithm, which can calculate excited states, is also unable to effectively handle degenerate states since the loss function only contains the expectation value of the Hamiltonian. In this study, the energy eigenstates of the one-dimensional Fermi-Hubbard model with two lattice sites and the two-dimensional Hubbard model with four lattice sites are calculated. By incorporating symmetries into the quantum circuits and loss function, we find that both the ground state and excited state calculations are improved greatly compared to the case without symmetries. The enhancement in excited state calculations is particularly significant. This is because quantum circuits that conserve the particle number are used, and appropriate penalty terms are added to the loss function, enabling the optimization process to correctly identify degenerate states. The results are verified through repeated simulations.

arXiv:2501.15903Jan 27, 20250 citations

Adaptive variational quantum eigensolvers for highly excited states

Feng Zhang, Niladri Gomes, Yongxin Yao +2 more

Highly excited states of quantum many-body systems are central objects in the study of quantum dynamics and thermalization that challenge classical computational methods due to their volume-law entanglement content. In this work, we explore the potential of variational quantum algorithms to approximate such states. We propose an adaptive variational quantum eigensolver (VQE) for excited states (X) that self-generates a variational ansatz for arbitrary eigenstates of a many-body Hamiltonian $H$ by attempting to minimize the energy variance with respect to $H$ . We benchmark the method by applying it to an Ising spin chain with integrable and nonintegrable regimes, where we calculate various quantities of interest, including the total energy, magnetization density, and entanglement entropy. We also compare the performance of adaptive VQE-X to an adaptive variant of the folded-spectrum method. For both methods, we find a strong dependence of the algorithm's performance on the choice of operator pool used for the adaptive construction of the ansatz. In particular, an operator pool including long-range two-body gates accelerates the convergence of both algorithms in the nonintegrable regime. We also study the scaling of the number of variational parameters with system size, finding that an exponentially large number of parameters may be necessary to approximate individual highly excited states. Nevertheless, we argue that these methods lay a foundation for the use of quantum algorithms to study finite-energy-density properties of many-body systems.Highly excited states of quantum many-body systems are central objects in the study of quantum dynamics and thermalization that challenge classical computational methods due to their volume-law entanglement content. In this work, we explore the potential of variational quantum algorithms to approximate such states. We propose an adaptive variational quantum eigensolver (VQE) for excited states (X) that self-generates a variational ansatz for arbitrary eigenstates of a many-body Hamiltonian $H$ by attempting to minimize the energy variance with respect to $H$. We benchmark the method by applying it to an Ising spin chain with integrable and nonintegrable regimes, where we calculate various quantities of interest, including the total energy, magnetization density, and entanglement entropy. We also compare the performance of adaptive VQE-X to an adaptive variant of the folded-spectrum method. For both methods, we find a strong dependence of the algorithm's performance on the choice of operator pool used for the adaptive construction of the ansatz. In particular, an operator pool including long-range two-body gates accelerates the convergence of both algorithms in the nonintegrable regime. We also study the scaling of the number of variational parameters with system size, finding that an exponentially large number of parameters may be necessary to approximate individual highly excited states. Nevertheless, we argue that these methods lay a foundation for the use of quantum algorithms to study finite-energy-density properties of many-body systems.

arXiv:2104.12636Physical review. B./Physical review. BAug 30, 202150 citations

Low-depth circuit ansatz for preparing correlated fermionic states on a quantum computer

Pierre-Luc Dallaire-Demers, Jonathan Romero, Libor Veis +2 more

Abstract Quantum simulations are bound to be one of the main applications of near-term quantum computers. Quantum chemistry and condensed matter physics are expected to benefit from these technological developments. Several quantum simulation methods are known to prepare a state on a quantum computer and measure the desired observables. The most resource economic procedure is the variational quantum eigensolver (VQE), which has traditionally employed unitary coupled cluster as the ansatz to approximate ground states of many-body fermionic Hamiltonians. A significant caveat of the method is that the initial state of the procedure is a single reference product state from a classical Hartree–Fock calculation with no pairing correlations, hence it cannot represent superconducting states. In this work, we propose to improve the method by initializing the algorithm with a more general fermionic Gaussian state, an idea borrowed from the field of nuclear physics. We show how this Gaussian reference state can be prepared with a linear-depth circuit of quantum matchgates. By augmenting the set of available gates with nearest-neighbor phase coupling, we generate a low-depth circuit ansatz that can accurately prepare the ground state of correlated fermionic systems. This extends the range of applicability of the VQE to systems with strong pairing correlations such as superconductors, atomic nuclei, and topological materials.Abstract Quantum simulations are bound to be one of the main applications of near-term quantum computers. Quantum chemistry and condensed matter physics are expected to benefit from these technological developments. Several quantum simulation methods are known to prepare a state on a quantum computer and measure the desired observables. The most resource economic procedure is the variational quantum eigensolver (VQE), which has traditionally employed unitary coupled cluster as the ansatz to approximate ground states of many-body fermionic Hamiltonians. A significant caveat of the method is that the initial state of the procedure is a single reference product state from a classical Hartree–Fock calculation with no pairing correlations, hence it cannot represent superconducting states. In this work, we propose to improve the method by initializing the algorithm with a more general fermionic Gaussian state, an idea borrowed from the field of nuclear physics. We show how this Gaussian reference state can be prepared with a linear-depth circuit of quantum matchgates. By augmenting the set of available gates with nearest-neighbor phase coupling, we generate a low-depth circuit ansatz that can accurately prepare the ground state of correlated fermionic systems. This extends the range of applicability of the VQE to systems with strong pairing correlations such as superconductors, atomic nuclei, and topological materials.

arXiv:1801.01053Quantum Science and TechnologyAug 7, 2019119 citations

Fermionic adaptive sampling theory for variational quantum eigensolvers

Marco Majland, Patrick Ettenhuber, Nikolaj Thomas Zinner

Quantum chemistry has been identified as one of the most promising areas where quantum computing can have a tremendous impact. For current noisy intermediate-scale quantum (NISQ) devices, one of the best available methods to prepare approximate wave functions on quantum computers is the adaptive derivative-assembled pseudo-Trotter Ansatz variational quantum eigensolver (ADAPT-VQE). However, ADAPT-VQE suffers from a significant measurement overhead when estimating the importance of operators in the wave function. In this work, we propose fermionic adaptive sampling theory VQE (FAST-VQE), a method for selecting operators based on importance metrics solely derived from the populations of Slater determinants in the wave function. Thus, our method mitigates measurement overheads for ADAPT-VQE as it is only dependent on the populations of Slater determinants which can simply be determined by measurements in the computational basis. We introduce two heuristic importance metrics, one based on selected configuration interaction with perturbation theory and one based on approximate gradients. In state vector and finite shot simulations, FAST-VQE using the heuristic metric based on approximate gradients converges at the same rate or faster than ADAPT-VQE and requires dramatically fewer shots.Quantum chemistry has been identified as one of the most promising areas where quantum computing can have a tremendous impact. For current noisy intermediate-scale quantum (NISQ) devices, one of the best available methods to prepare approximate wave functions on quantum computers is the adaptive derivative-assembled pseudo-Trotter Ansatz variational quantum eigensolver (ADAPT-VQE). However, ADAPT-VQE suffers from a significant measurement overhead when estimating the importance of operators in the wave function. In this work, we propose fermionic adaptive sampling theory VQE (FAST-VQE), a method for selecting operators based on importance metrics solely derived from the populations of Slater determinants in the wave function. Thus, our method mitigates measurement overheads for ADAPT-VQE as it is only dependent on the populations of Slater determinants which can simply be determined by measurements in the computational basis. We introduce two heuristic importance metrics, one based on selected configuration interaction with perturbation theory and one based on approximate gradients. In state vector and finite shot simulations, FAST-VQE using the heuristic metric based on approximate gradients converges at the same rate or faster than ADAPT-VQE and requires dramatically fewer shots.

arXiv:2303.07417Physical Review ANov 21, 20237 citations

Qubit-excitation-based adaptive variational quantum eigensolver

Yordan S. Yordanov, V. Armaos, Crispin H. W. Barnes +1 more

Abstract Molecular simulations with the variational quantum eigensolver (VQE) are a promising application for emerging noisy intermediate-scale quantum computers. Constructing accurate molecular ansatze that are easy to optimize and implemented by shallow quantum circuits is crucial for the successful implementation of such simulations. Ansatze are, generally, constructed as series of fermionic-excitation evolutions. Instead, we demonstrate the usefulness of constructing ansatze with ``qubit-excitation evolutions', which, contrary to fermionic excitation evolutions, obey ``qubit commutation relations'. We show that qubit excitation evolutions, despite the lack of some of the physical features of fermionic excitation evolutions, accurately construct ansatze, while requiring asymptotically fewer gates. Utilizing qubit excitation evolutions, we introduce the iterative qubit excitation based VQE (IQEB-VQE) algorithm. The IQEB-VQE performs molecular simulations using a problem-tailored ansatz, grown iteratively by appending evolutions of single and double qubit excitation operators. By performing numerical simulations for small molecules, we benchmark the IQEB-VQE, and compare it against other competitive VQE algorithms. In terms of circuit efficiency and time complexity, we find that the IQEB-VQE systematically outperforms the previously most circuit-efficient, practically scalable VQE algorithms.Abstract Molecular simulations with the variational quantum eigensolver (VQE) are a promising application for emerging noisy intermediate-scale quantum computers. Constructing accurate molecular ansatze that are easy to optimize and implemented by shallow quantum circuits is crucial for the successful implementation of such simulations. Ansatze are, generally, constructed as series of fermionic-excitation evolutions. Instead, we demonstrate the usefulness of constructing ansatze with ``qubit-excitation evolutions', which, contrary to fermionic excitation evolutions, obey ``qubit commutation relations'. We show that qubit excitation evolutions, despite the lack of some of the physical features of fermionic excitation evolutions, accurately construct ansatze, while requiring asymptotically fewer gates. Utilizing qubit excitation evolutions, we introduce the iterative qubit excitation based VQE (IQEB-VQE) algorithm. The IQEB-VQE performs molecular simulations using a problem-tailored ansatz, grown iteratively by appending evolutions of single and double qubit excitation operators. By performing numerical simulations for small molecules, we benchmark the IQEB-VQE, and compare it against other competitive VQE algorithms. In terms of circuit efficiency and time complexity, we find that the IQEB-VQE systematically outperforms the previously most circuit-efficient, practically scalable VQE algorithms.

arXiv:2011.10540Communications PhysicsOct 14, 2021145 citations

H

by attempting to minimize the energy variance with respect to

H

. We benchmark the method by applying it to an Ising spin chain with integrable and nonintegrable regimes, where we calculate various quantities of interest, including the total energy, magnetization density, and entanglement entropy. We also compare the performance of adaptive VQE-X to an adaptive variant of the folded-spectrum method. For both methods, we find a strong dependence of the algorithm's performance on the choice of operator pool used for the adaptive construction of the ansatz. In particular, an operator pool including long-range two-body gates accelerates the convergence of both algorithms in the nonintegrable regime. We also study the scaling of the number of variational parameters with system size, finding that an exponentially large number of parameters may be necessary to approximate individual highly excited states. Nevertheless, we argue that these methods lay a foundation for the use of quantum algorithms to study finite-energy-density properties of many-body systems.