Meh, muy interesante integral! Me puede dar una aproximación heurística, pero creo que alguien va a hacer mejor. Estoy en el autobús y usted sabe, no es fácil.
Me gustaría utilizar la Serie de Taylor para $e^{-a/t}$, por lo tanto
$$\int_0^{+\infty}\sum_{k = 0}^{+\infty} \frac{\left(-a/t\right)^k}{k!}t^{-3/2}\ \text{Erf}(\sqrt{t})\ \text{d}t$$
Y llegamos
$$\sum_{k = 0}^{+\infty}\frac{(-a)^k}{k!}\int_0^{+\infty} t^{-k - 3/2}\ \text{Erf}(\sqrt{t})\ \text{d}t$$
Ahora bien, si llamamos a la simplicidad $b = -k - 3/2$ obtenemos una computable integral (he comprobado en Mathematica), que dice:
$$\int_0^{+\infty} t^{b}\ \text{Erf}(\sqrt{t})\ \text{d}t = -\frac{\Gamma[3/2 + b]}{(1 + b)\sqrt{\pi}} ~~~ \to ~~~ -\frac{\Gamma[-k]}{(-k - 1/2)\sqrt{\pi}}$$
PERO no es la condición a través de este resultado:
$$-\frac{3}{2} < \Re(b) < -1$$
Esto daría entonces
$$-\ \sum_{k = 0}^{+\infty}\frac{(-a)^k}{k!}\frac{\Gamma[-k]}{(-k - 1/2)\sqrt{\pi}}$$
Y aquí me paro porque yo no puedo ir (sobre todo porque no tengo papel y lápices conmigo.. me voy a revisar de nuevo cuando voy a estar en casa).