comienza con una ecuación homogénea
$$
\frac{\partial{}G(x,t)}{\partial{}t} - \frac{\partial{}^2G(x,t)}{\partial{}x^2} =0
$$
la solución es
$$ G(x,t) = \frac{\mathrm{exp}\left( -\frac{{{x}^{2}}}{ 4 \, a \, t}\right) }{\sqrt{t \; 2 \; \pi }}$$
El área debajo de la gráfica de la función es
$$
A_t = \frac{\int_{-\infty }^{\infty }{{e}^{-\frac{{{x}^{2}}}{4\cdot a\cdot t}}}dx}{\sqrt{2 \pi t}} = \sqrt{2 \pi}
$$
Por lo tanto,
$$
\lim_{t \rightarrow 0^+} \frac{G(x,t)}{A_t} =\delta(t)
$$
Observe que no hay ningún producto de delta.
A continuación, busca una solución en una forma de convolución
$$
T(x,t) = \int_{ -\infty }^\infty \int_{-\infty}^\infty f(x' -x,\tau'-t) G(x',t') dx' \, dt'
$$
y requiere que ambas satisfacen la misma ecuación.
Que es
$$ \frac{\int_{-\infty }^{\infty }\frac{\left( \frac{d}{d\,t} \mathrm{f}\left( x-z,t-u\right) \right) {{e}^{-\frac{{{z z z}^{2}}}{4\cdot a\cdot u}}}}{\sqrt{u}}du}{\sqrt{2 \, \pi }}-\frac{a \int_{-\infty }^{\infty }\frac{\left( \frac{{{d}^{2}}}{d\,{{x}^{2}}} \mathrm{f}\left( x-z,t-u\right) \right) \cdot {{e}^{-\frac{{{z z z}^{2}}}{4\cdot a\cdot u}}}}{\sqrt{u}}du}{\sqrt{2\,\pi }}= \sqrt{2 \pi} \delta(x - z)
$$