I would like to know how the integral \begin{equation} \int_{-\infty}^{+\infty}\frac{e^{-x^2}}{e^{-x} + 1}dx \end{equation} can be evaluated. Mathematica says the value is $\sqrt{\pi}/2$. I thought I could get this using Residue theorem but I'm having trouble showing that the arc \begin{equation} \int_{0}^\pi\frac{e^{-R^2e^{2i\theta}}}{e^{-Re^{i\theta}} + 1}iRe^{i\theta}d\theta \end{equation} would go to zero as $R\rightarrow 0$, mainly because $e^{-R^2e^{2i\theta}}$ is actually large when $\theta$ is around $\pi/2$.
Ultimately, I'm actually interested in a slightly more general integral which Mathematica won't evaluate at all: \begin{equation} \int_{-\infty}^{+\infty}\frac{e^{-a(x-b)^2}}{e^{-x} + 1}dx. \end{equation} So I'm hoping that understanding the simple case above will help me.