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