Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Abstract

Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to $32$ qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.

References in corpus (50)

• Variational Quantum Algorithms
• The theory of variational hybrid quantum-classical algorithms
• Complexity in Strongly Correlated Electronic Systems
• Quantum Cluster Theories
• Theory of Intertwined Orders in High Temperature Superconductors
• The Variational Quantum Eigensolver: a review of methods and best practices
• The ITensor Software Library for Tensor Network Calculations
• An adaptive variational algorithm for exact molecular simulations on a quantum computer
• Frequency-dependent local interactions and low-energy effective models from electronic structure calculations
• Towards Practical Quantum Variational Algorithms
• The Hubbard Model
• Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms
• qubit-ADAPT-VQE: An adaptive algorithm for constructing hardware-efficient ansatze on a quantum processor
• Order and quantum phase transitions in the cuprate superconductors
• Superconductivity and the Pseudogap in the two-dimensional Hubbard model
• Solving strongly correlated electron models on a quantum computer
• Strongly correlated excitonic insulator in atomic double layers
• Strongly Correlated Superconductivity
• Barren Plateaus in Variational Quantum Computing
• Superconductivity in the Hubbard model and its interplay with charge stripes and next-nearest hopping t'
• Strategies for solving the Fermi-Hubbard model on near-term quantum computers
• Antiferromagnetic correlations in the metallic strongly correlated transition metal oxide LaNiO${}_3$
• Quantum Monte Carlo study of the two-dimensional fermion Hubbard Model
• Variational Quantum State Eigensolver
• RESPACK: An ab initio tool for derivation of effective low-energy model of material
• Two-band model for magnetism and superconductivity in nickelates
• Observing ground-state properties of the Fermi-Hubbard model using a scalable algorithm on a quantum computer
• Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method
• Strong Coupling Nature of the Excitonic Insulator State in Ta${}_2$NiSe${}_5$
• Band insulator to Mott insulator transition in a bilayer Hubbard model
• Adaptive Variational Quantum Imaginary Time Evolution Approach for Ground State Preparation
• Excitonic Instability at the Spin-State Transition in the Two-Band Hubbard Model
• Ab Initio Full Cell GW+DMFT for Correlated Materials
• Resource Estimation for Quantum Variational Simulations of the Hubbard Model
• Benchmarking adaptive variational quantum eigensolvers
• Finding the ground state of the Hubbard model by variational methods on a quantum computer with gate errors
• Gutzwiller Hybrid Quantum-Classical Computing Approach for Correlated Materials
• Ab initio downfolding study of the iron-based ladder superconductor BaFe${}_2$S${}_3$
• From real materials to model Hamiltonians with density matrix downfolding
• Real-Time Krylov Theory for Quantum Computing Algorithms
• Dynamical Mean-Field Theory Simulations with the Adaptive Sampling Configuration Interaction Method
• A method to efficiently simulate the thermodynamical properties of the Fermi-Hubbard model on a quantum computer
• Metal-insulator transitions, superconductivity, and magnetism in the two-band Hubbard model
• Entanglement of localized states
• Comparative study of adaptive variational quantum eigensolvers for multi-orbital impurity models
• Competition between band and Mott insulator in the bilayer Hubbard model: a dynamical cluster approximation study
• Simulating strongly interacting Hubbard chains with the Variational Hamiltonian Ansatz on a quantum computer
• Quantum subspace expansion in the presence of hardware noise
• Compressing Hamiltonians with ab initio downfolding for simulating strongly-correlated materials on quantum computers
• From Heisenberg to Hubbard: An initial state for the shallow quantum simulation of correlated electrons

Cited by in corpus (2)

• Efficient Fermi-Hubbard model ground-state preparation by coupling to a classical reservoir in the instantaneous-response limit
• Variational Quantum Subspace Construction via Symmetry-Preserving Cost Functions