An approximate Riemannian heat kernel with the correct delta initial limit and a residual extending with enough regularity to time zero. A cutoff Gaussian times transport coefficients supplies a local construction. Increasing the expansion order, or summing the full asymptotic series smoothly, makes the residual regular enough for a Volterra parametrix correction.
New to topics? Read the docs here!