Nash estimates and upper bounds for non-homogeneous Kolmogorov equations

We prove a Gaussian upper bound for the fundamental solutions of a class of ultra-parabolic equations in divergence form. The bound is independent on the smoothness of the coefficients and generalizes some classical results by Nash, Aronson and Davies. The class considered has relevant applications in the theory of stochastic processes, in physics and in mathematical finance.