The famous, yet unsolved, Fermi-Hubbard model for strongly-correlated electronic systems is a prominent target for quantum computers. However, accurately representing the Fermi-Hubbard ground state for large instances may be beyond the reach of near-term quantum hardware. Here we show experimentally that an efficient, low-depth variational quantum algorithm with few parameters can reproduce important qualitative features of medium-size instances of the Fermi-Hubbard model. We address 1 × 8 and 2 × 4 instances on 16 qubits on a superconducting quantum processor, substantially larger than previous work based on less scalable compression techniques, and going beyond the family of 1D Fermi-Hubbard instances, which are solvable classically. Consistent with predictions for the ground state, we observe the onset of the metal-insulator transition and Friedel oscillations in 1D, and antiferromagnetic order in both 1D and 2D. We use a variety of error-mitigation techniques, including symmetries of the Fermi-Hubbard model and a recently developed technique tailored to simulating fermionic systems. We also introduce a new variational optimisation algorithm based on iterative Bayesian updates of a local surrogate model.
The Fermi-Hubbard model is a plausible target to be solved by a quantum computer using the variational quantum eigensolver algorithm. However, problem sizes beyond the reach of classical exact diagonalisation are also beyond the reach of current quantum computing hardware. Here we use a simple method which compresses the first nontrivial subcase of the Hubbard model -- with one spin-up and one spin-down fermion -- enabling larger instances to be addressed using current quantum computing hardware. We implement this method on a superconducting quantum hardware platform for the case of the
2×1 Hubbard model, including error-mitigation techniques, and show that the ground state is found with relatively high accuracy.
It is expected that the simulation of correlated fermions in chemistry and material science will be one of the first practical applications of quantum processors. Given the rapid evolution of quantum hardware, it is increasingly important to develop robust benchmarking techniques to gauge the capacity of quantum hardware specifically for the purpose of fermionic simulation. Here we propose using the one-dimensional Fermi-Hubbard model as an application benchmark for variational quantum simulations on near-term quantum devices. Since the one-dimensional Hubbard model is both strongly correlated and exactly solvable with the Bethe ansatz, it provides a reference ground state energy that a given device with limited coherence will be able to approximate up to a maximal size. The length of the largest chain that can be simulated provides an effective fermionic length. We use variational quantum eigensolver to approximate the ground state energy values of Fermi-Hubbard instances and show how the fermionic length benchmark can be used in practice to assess the performance of bounded-depth devices in a scalable fashion.
We report experimental digital quantum simulation of the one-dimensional Fermi-Hubbard model on a superconducting quantum processor at a scale beyond the reach of exact statevector simulation and challenging for state-of-the-art tensor-network methods. We encode this problem using up to 120 qubits through an efficient mapping that reduces circuit complexity, and we improve accuracy through error suppression to simulate dynamical evolution using up to 90 Trotter steps. From a vacancy defect introduced in the middle of an
L=31-site (62-qubit) Néel initial state, we directly observe spin-charge separation to
t=9 in natural units and quantitatively extract velocities
vc,
vs, which match classical simulations across model parameters, and analytics in appropriate regimes. We then extend experiments to
L=60 fermions (120 qubits) and evolution times to
t=6 using 30 Trotter steps. Quantum-processor outputs agree quantitatively with approximate classical simulations performed using a time-dependent variational principle (TDVP) solver and we observe that increasing the TDVP bond dimension through
χ=4096 expands the range of evolution times within which agreement has RMSE
∼1% before the approaches diverge. Owing to the large scale of the simulation and the use of overhead-free error-suppression techniques, for simulated evolution times at the limit of quantum/classical agreement (
t≳5 in natural hopping units), the wall-clock runtime of the quantum processor is nearly three orders of magnitude faster than TDVP simulations with
χ=4096 and implemented using the ITensor and TeNPy packages. These results establish contemporary digital quantum processors as a versatile, quantitatively accurate, and competitive platform for the study of fermionic many-body dynamics in regimes where leading classical methods can become prohibitively expensive.
Quantum simulation and near-term quantum algorithms have gained significant momentum over the past few years, with researchers exploring quantum-based approaches to solve classically intractable problems. Out of many research problems, certain problems are well suited for quantum algorithms. Fermi Hubbard model based on lattice structures is one of the quantum-mechanical models used in many-body and condensed matter physics to capture the essence of electron-correlated systems. In this work, we show one of the first kinds of experiments on simulating the ground state properties of 1D and 2D Fermi Hubbard-based lattice models compared with the reference ground state energy values using Qiskit Nature. We utilize the Variational Quantum Eigensolver (VQE) and different gradient-based optimizers to achieve results that closely converge towards the exact values, showcasing the potential of variational quantum algorithms.
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.
As the advances in quantum hardware bring us into the noisy intermediate-scale quantum (NISQ) era, one possible task we can perform without quantum error correction using NISQ machines is the variational quantum eigensolver (VQE) due to its shallow depth. A specific problem that we can tackle is the strongly interacting Fermi-Hubbard model, which is classically intractable and has practical implications in areas like superconductivity. In this paper, we outline the details about the gate sequence, the measurement scheme, and the relevant error-mitigation techniques for the implementation of the Hubbard VQE on a NISQ platform. We perform resource estimation for both silicon spin qubits and superconducting qubits for a 50-qubit simulation, which cannot be solved exactly via classical means, and find similar results. The number of two-qubit gates required is on the order of 20 000. Hence, to suppress the mean circuit-error count to a level such that we can obtain meaningful results with the aid of error mitigation, we need to achieve a two-qubit gate error rate of approximately ${10}^{\ensuremath{-}4}$. When searching for the ground state, we need a few days for one gradient-descent iteration, which is impractical. This can be reduced to around 10 min if we distribute our task among hundreds of quantum-processing units. Hence, implementing a 50-qubit Hubbard model VQE on a NISQ machine can be on the brink of being feasible in near term, but further optimization of our simulation scheme, improvements in the gate fidelity, improvements in the optimization scheme and advances in the error-mitigation techniques are needed to overcome the remaining obstacles. The scalability of the hardware platform is also essential to overcome the runtime issue via parallelization, which can be done on one single silicon multicore processor or across multiple superconducting processors.