We adapt in the present note the perturbation method introduced in Choulli and Kayser (Positivity 19(3):625–646, 2015) to get a lower Gaussian bound for the Neumann heat kernel of the Laplace–Beltrami operator on an open subset of a compact Riemannian manifold.