EDITAR : Originalmente hice la integral sobre $(-\infty,\infty)$ y obtuvo una solución elemental, en lugar de sobre el intervalo $[0,\infty)$ donde obtenemos funciones especiales. Esto se ha modificado.
Definir $I(a,b) = \int_0^{+\infty} \exp(-a t^2-bt)dt$ para $a\in \mathbb{R}_+$ y $b\in\mathbb{C}$ . Si completamos el cuadrado, obtenemos
$$I(a,b) = \int_0^{+\infty} \exp\left[ -a\left(t+\frac{b}{2a}\right)^2+\frac{b^2}{4a} \right]dt$$
$$=\exp\left(\frac{b^2}{4a}\right) \int_0^{+\infty} \exp\left[ -a\left(t+\frac{b}{2a}\right)^2\right]dt$$
$$=\exp\left(\frac{b^2}{4a}\right)\int_{b/2a}^\infty e^{-au^2}du$$
$$=\frac{1}{2} \sqrt{\frac{\pi}{a}}\exp\left(\frac{b^2}{4a}\right) \mathrm{erfc}\left(\frac{b}{2\sqrt{a}}\right).$$
En lo anterior, $\mathrm{erfc}(\cdot)$ denota el complemento función de error (que se define para los números complejos). Ahora observa que el seno de la integral original se puede dividir en exponenciales complejas, dando
$$(*)=\frac{I(a^2,b-ci)-I(a^2,b+ci)}{2i}.$$