Convex-Gaussianity of fermionic Gibbs states in perturbation theory
We study the structure of Gibbs states in weakly perturbed interacting fermionic systems. First, for a sparse Hamiltonian $H=H_0+V$ with a quadratic term $H_0$ and a non-quadratic perturbation $V$ of scale $\epsilon$, we show that the Gibbs state $\rho_{\beta}$ decomposes into a convex combination of Gaussian states wh...