Primero, tengamos en cuenta que, debido a la rareza de la función seno, tenemos $(-x)\sin(-x)=(-x)\cdot$ $(-\sin x)=x\sin x\iff I=2\int_0^\infty xe^{-\frac{x^2}2}\sin kx\,dx$. En segundo lugar, desde la $\sin kx=\Im(e^{ikx})$, $I=2\cdot\Im\left[\int_0^\infty xe^{-\left(\frac{x^2}2-ikx\right)}dx\right]$ donde $\Im(z)=\Im(a+bi)=b$. Ahora, vamos a pagar un vistazo más de cerca en el exponente : $\frac{x^2}2-ikx=\frac12(x^2-2ikx)=\frac12\left[(x-ik)^2+k^2\right]=\frac{t^2+k^2}2$ donde $t=x-ik$, y $dt$ $=dx$. Entonces:
$$I=2\,\Im\left[\int_{0-ik}^{\infty-ik}(t+ik)e^{-\frac{t^2+k^2}2}dt\right]=2\,\Im\left[\int_{0-ik}^{\infty-ik}te^{-\frac{t^2+k^2}2}dt+ik\int_{0-ik}^{\infty-ik}e^{-\frac{t^2+k^2}2}dt\right]=$$
$$=2\,\Im\left[\int_0^\infty e^{-u}du+ike^{-\frac{k^2}2}\int_{0-ik}^{\infty-ik}e^{-\frac{t^2}2}dt\right]=2\,\Im\left[1+ike^{-\frac{k^2}2}\sqrt\frac\pi2\left(1+\text{Erf}\left(\tfrac{ki}{\sqrt2}\right)\right)\right]=$$
$=ke^{-\frac{k^2}2}\sqrt{2\pi}$ , ya que la función de error puramente imaginario argumento es puramente imaginaria así, lo que significa que yo veces sí no posee ninguna parte imaginaria, y por lo tanto no influyen en el resultado final.