We derive upper bounds on the fluctuations of a class of random surfaces of the ∇φ-type with convex interaction potentials. The Brascamp-Lieb concentration inequality provides an upper bound on these fluctuations for uniformly convex potentials. We extend these results to twice continuously differentiable convex potentials whose second derivative grows asymptotically like a polynomial and may vanish on an (arbitrarily large) interval. Specifically, we prove that, when the underlying graph is the d-dimensional torus of side length L, the variance of the height is smaller than C ln L in two dimensions and remains bounded in dimension d ≥ 3.The proof makes use of the Helffer-Sjöstrand representation formula (originally introduced by Helffer and Sjöstrand (1994) and used by Naddaf andSpencer (1997) andGiacomin, Olla Spohn (2001) to identify the scaling limit of the model), the anchored Nash inequality (and the corresponding on-diagonal heat kernel upper bound) established by Mourrat and Otto (2016) and Efron's monotonicity theorem for log-concave measures (Efron (1965)). Contents 1 Introduction 1 2 Notation and preliminary results 8 3 Fluctuation estimates for the Langevin dynamic 13 4 On diagonal upper bound for the heat kernel 21 References 36