Nous revisitons la méthode classique des paramétrices pour en déduire une minoration et une majoration gaussiennes, pour la solution fondamentale d'un opérateur parabolique général sous forme non divergentielle. Nous utilisons ensuite le fait que la fonction de Neumann Green, d'un opérateur parabolique général sur un ouvert borné régulier, peut être construite comme somme de la solution fondamentale et d'une intégrale de type simple couche parabolique pour établir une minoration gaussienne pour cette fonction de Neumann Green. Le point clef de la preuve réside dans l'effet régularisant, en temps, de l'intégrale de type simple couche. Nous démontrons aussi que cette approche peut être adaptée pour démontrer une minoration gaussienne pour la fonction de Green-Neumann correspondante à l'opérateur de Laplace-Beltrami sur un ouvert régulier d'une variété riemannienne compacte sans bord. Nous démontrons ensuite une nouvelle majoration gaussienne pour la fonction de Neumann Green correspondante à l'opérateur de Laplace-Beltrami sur un ouvert Lipschitz d'une variété riemannienne complète. L'intérêt de cette nouvelle majoration est qu'elle ne contient pas le terme habituel d'une exponentielle en temps. Finalement, comme application des estimations gaussiennes, nous donnons un résultat de compacité des potentiels isospectraux en relation avec une formule asymptotique pour les noyaux de la chaleur