🤖 AI Summary
This work investigates the structure of Gibbs states of bounded-degree local fermionic Hamiltonians at high temperatures. We prove that, above a system-size-independent temperature threshold, the Gibbs state admits an exact decomposition as a probabilistic mixture of fermionic Gaussian states—a result constituting the first rigorous demonstration that high-temperature fermionic Gibbs states must possess such Gaussian-mixture structure, thereby revealing their intrinsic classical simulability. Methodologically, our approach integrates operator-norm estimates, analytic continuation of thermal states, and the theory of fermionic quadratic forms; we further design a probabilistic sampling scheme and a Gaussian-state preparation algorithm. Consequently, we construct the first polynomial-time classical algorithm capable of efficiently sampling from and explicitly preparing these high-temperature Gibbs states—overcoming the conventional exponential resource scaling and establishing a new paradigm for classical simulation of thermal equilibrium in strongly correlated fermionic systems.
📝 Abstract
Efficient simulation of a quantum system generally relies on structural properties of the quantum state. Motivated by the recent results by Bakshi et al. on the sudden death of entanglement in high-temperature Gibbs states of quantum spin systems, we study the high-temperature Gibbs states of bounded-degree local fermionic Hamiltonians, which include the special case of geometrically local fermionic systems. We prove that at a sufficiently high temperature that is independent of the system size, the Gibbs state is a probabilistic mixture of fermionic Gaussian states. This forms the basis of an efficient classical algorithm to prepare the Gibbs state by sampling from a distribution of fermionic Gaussian states.