Track the frontier, from condensed matter to cosmology.
Search papers, labs, and topics across Lattice.
This work resolves a fundamental question in quantum many-body physics by determining when weakly interacting fermionic thermal states admit efficient classical simulation. The authors prove that the Gibbs state of any local Hamiltonian $H = H_0 + V$ decomposes into a convex mixture of fermionic Gaussian states at any finite temperature and in arbitrary spatial dimensions, provided the interaction strength falls below a non-zero threshold. Exploiting this convex structure, they construct a polynomial-time classical sampling algorithm using Metropolis updates over an ensemble of tree-structured Gaussian states.
Interacting fermions do not imply quantum hardness: thermal Gibbs states in any dimension are provably just convex mixtures of Gaussian states below a non-zero interaction threshold, admitting efficient classical sampling.
Identifying regimes in which interacting fermionic systems admit efficient classical simulation is a crucial task in quantum many-body physics. In this work, we rigorously establish the Gaussianity and classical simulability of sufficiently weakly interacting fermions. For the Gibbs state of a local interacting fermionic Hamiltonian $H=H_0+V$ at any finite temperature and in any spatial dimension, where $H_0$ is quadratic and $V$ is the interaction term, we show that there exists a positive interaction strength below which the Gibbs state is a convex mixture of fermionic Gaussian states. Moreover, based on this structure, we construct a polynomial-time classical sampler for the Gibbs state via Metropolis updates on a tree whose leaves are Gaussian states.