The aim of this paper is twofold. On the one hand, the study of gradient Schrödinger operators on manifolds with density φ. We classify the space of solutions when the underlying manifold is φ−parabolic. As an application, we extend the Naber-Yau Liouville Theorem, and we will prove that a complete manifold with density is φ−parabolic if, and only if, it has finite φ−capacity. Moreover, we show that the linear space given by the kernel of a nonnegative gradient Schrödinger operators is one dimensional provided there exists a bounded function on it and the underlying manifold is φ−parabolic.On the other hand, the topological and geometric classification of complete weighted H φ −stable hypersurfaces immersed in a manifold with density (N , g, φ) satisfying a lower bound on its Bakry-Émery-Ricci tensor. Also, we classify weighted stable surfaces in a threemanifold with density whose Perelman scalar curvature, in short, P-scalar curvature, satisfies R ∞ φ + |∇φ| 2 4 ≥ 0. Here, the P-scalar curvature is defined as R ∞ φ = R − 2∆ g φ − |∇ g φ| 2 , being R the scalar curvature of (N , g).Finally, we discuss the relationship of manifolds with density, Mean Curvature Flow (MCF), Ricci Flow and Optimal Transportation Theory. We obtain classification results for stable self-similiar solutions to the MCF, and also for stable translating solitons to the MCF (as far as we know, the first classification result). Moreover, for gradient Ricci solitons, we recover the Hamilton-Ivey-Perelman classification assuming only a L 2 −type bound on the scalar curvature and a inequality between the scalar curvature and the Ricci tensor. We also classify gradient Ricci solitons when the scalar curvature does not change sign. We finish by classifying critical transportation plans for the Boltzman entropy on φ−parabolic manifolds. In particular, we recover the case for the Gaussian measure.