Potential analysis for a class of diffusion equations: a Gaussian bounds approach