The asymptotic decay rate and other qualitative features of solutions to parabolic equations set in noncompact domains of RN, or Riemannian manifolds, depend in general on a suitable notion of the geometry of the domain at infinity, as well as on the structure of the equation. We report on some recent and on-going work on the interplay of these two factors (joint project with prof. A.Tedeev).