The role of fundamental solution in Potential and Regularity Theory for subelliptic PDE

In this survey we consider a general Hormander type operator, represented as a sum of squares of vector fields plus a drift and we outline the central role of the fundamental solution in developing Potential and Regularity Theory for solutions of related PDEs. After recalling the Gaussian behavior at infinity of the kernel, we show some mean value formulas on the level sets of the fundamental solution, which are the starting point to obtain a comprehensive parallel of the classical Potential Theory. Then we show that a precise knowledge of the fundamental solution leads to global regularity results, namely estimates at the boundary or on the whole space. Finally in the problem of regularity of non linear differential equations we need an ad hoc modification of the parametrix method, based on the properties of the fundamental solution of an approximating problem.