Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Antonios M. Alvertis,Abid Khan,Thomas Iadecola,Peter P. Orth,Norm Tubman

2024-08-02

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. 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.

Quantum Physics,Strongly Correlated Electrons

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to evaluate the accuracy of the Variational Quantum Eigensolver (VQE) when simulating the Hubbard model, with particular attention to the following aspects:

1. **Influence of system size**: Study the accuracy of VQE solutions under different system sizes, especially whether the error of the variational solution will tend to be stable or increase under larger lattice sizes.
2. **Influence of electron - interaction strength**: Explore the influence of different electron - interaction strengths \( U/t \) on the accuracy of VQE simulations, especially the performance in the strong - coupling region \( U/t = 8 \).
3. **Choice of wave - function variational form (Ansatz)**: Compare the influence of different wave - function variational forms (such as Number Preserving (NP) ansatz, Excitation Preserving (EP) ansatz and Unitary Coupled Cluster (UCC) ansatz) on VQE simulation results.
4. **Influence of spatial inhomogeneity and off - site Coulomb interaction**: Study how the accuracy of VQE simulations changes after introducing spatial inhomogeneity parameters and off - site Coulomb interactions.

Through these studies, the author aims to provide references and guidance for future use of quantum computers to solve complex material problems, especially for the accurate simulation of the Hubbard model.