Locally accurate matrix product approximation to thermal states

In one-dimensional quantum systems with short-range interactions, we prove that a thermal state at inverse temperature $\beta=O(1)$ has a matrix product representation with bond dimension $e^{\tilde O(\sqrt{\beta\log(1/\epsilon)})}$ such that all local properties are approximated to accuracy $\epsilon$.