Polynomial Time Quantum Gibbs Sampling for Fermi-Hubbard Model at any Temperature

TL;DR

This paper develops a quantum algorithm for efficiently preparing Gibbs states of fermionic systems at any temperature, including interacting models, with complexity bounds independent of system size.

Contribution

It calculates the spectral gap of the Lindbladian for free fermions, proves a constant gap for interacting fermions, and demonstrates the algorithm's efficiency for large systems.

Findings

Spectral gap of Lindbladian for free fermions computed using third quantisation.

Logarithmic bound on mixing time for free fermions established.

Constant gap proven for perturbed Lindbladian in interacting fermions, enabling efficient Gibbs state preparation.

Abstract

Recently, there have been several advancements in quantum algorithms for Gibbs sampling. These algorithms simulate the dynamics generated by an artificial Lindbladian, which is meticulously constructed to obey a detailed-balance condition with the Gibbs state of interest, ensuring it is a stationary point of the evolution, while simultaneously having efficiently implementable time steps. The overall complexity then depends primarily on the mixing time of the Lindbladian, which can vary drastically, but which has been previously bounded in the regime of high enough temperatures [Rouz\'e et al. arXiv:2403.12691 and arXiv:2411.04885]. In this work, we calculate the spectral gap of the Lindbladian for free fermions using third quantisation, and also prove a logarithmic bound on its mixing time by analysing corresponding covariance matrices. Then we prove a constant gap of the perturbed Lindbladian corresponding to interacting fermions up to some maximal coupling strength. This is achieved by using theorems about stability of the gap for lattice fermions. Our methods apply at any constant temperature and independently of the system size. The gap then provides an upper bound on the mixing time, and hence on the overall complexity of the quantum algorithm, proving that the purified Gibbs state of weakly interacting (quasi-)local fermionic systems of any dimension can be prepared in $\widetilde{\mathcal{O}} (n^3 \operatorname{polylog}(1/\epsilon))$ time on $\mathcal{O}(n)$ qubits, where $n$ denotes the size of the system and $\epsilon$ the desired accuracy. As an application, we explain how to calculate partition functions for the considered systems. We provide exact numerical simulations for small system sizes supporting the theory and also identify different suitable jump operators and filter functions for the sought-after regime of intermediate coupling in the Fermi-Hubbard model.