Height fluctuations in non-integrable classical dimers

We rigorously establish the asymptotic equivalence between the height function of interacting dimers on the square lattice and the massless Gaussian free field. Our theorem explains the microscopic origin of the sine-Gordon field theory description away from the free fermion point, which has previously been elusive. We use a novel technique, based on the combination of discrete holomorphicity with exact, constructive, renormalization group methods, which has the potential of being applicable to a variety of other non-integrable models at or close to criticality.