On the integral $\int_{e}^{\infty}\frac{t^{1/2}}{\log^{1/2}\left(t\right)}\alpha^{-t/\log\left(t\right)}dt,\,\alpha>1.$ Let $\alpha>1$. I would like to find a closed form or an upper bound of 

$$f\left(\alpha\right)=\int_{e}^{\infty}\frac{t^{1/2}}{\log^{1/2}\left(t\right)}\alpha^{-t/\log\left(t\right)}dt.$$ 

For the closed form I'm very skeptical but I have trouble also for an upper bound. I tried, manipulating a bit, to integrate w.r.t. $\alpha$ since $$\frac{\partial}{\partial\alpha}\alpha^{-t/\log\left(t\right)}=-\frac{t}{\alpha\log\left(t\right)}\alpha^{-t/\log\left(t\right)}$$ but it seems quite useless and at this moment I didn't see a good way to proceed. 
Maybe it is interesting to see, using some trivial substitutions, that $$f\left(\alpha\right)=\int_{e}^{\infty}\frac{\left(e^{3/2}\right)^{-W_{-1}\left(-1/v\right)}}{v\left(-W_{-1}\left(-\frac{1}{v}\right)\right){}^{1/2}}\frac{W_{-1}\left(-\frac{1}{v}\right)}{W_{-1}\left(-\frac{1}{v}\right)+1}\alpha^{-v}dv$$ $$=\int_{e}^{\infty}g\left(w\right)\alpha^{-v}dv$$ where $W_{-1}\left(x\right)$ is the Lambert $W$ function. So it seems that $f(\alpha)$ is somehow connected to the Mellin transform of $g(w).$
Thank you.
 A: A naif but probably efficient approach is to exploit the fact that the logarithm function is approximately constant on short intervals and
$$ \frac{1}{\sqrt{N}}\int_{e^N}^{e^{N+1}}\sqrt{t}\,\alpha^{-t/N}\,dt =\frac{N\sqrt{\pi}}{2\log(\alpha)^{3/2}}\,\text{Erf}\left(\sqrt{\frac{e^N\log\alpha}{N}}\right)$$
can be efficiently approximated through the continued fraction for the error function.
We may also consider this fact: through the Laplace transform
$$ \int_{0}^{+\infty}\sqrt{t}\exp\left(-\frac{t\log\alpha}{N}\right)\,dt = \int_{0}^{+\infty}\mathcal{L}^{-1}\left(\frac{1}{\sqrt{t}}\right)\,\mathcal{L}\left(t \exp\left(-\frac{t\log\alpha}{N}\right)\right)\,ds $$
we get the following integral:
$$ \int_{0}^{+\infty}\frac{N^2}{\sqrt{\pi s}(Ns+\log\alpha)^2}\,ds =\frac{2}{\sqrt{\pi}}\int_{0}^{+\infty}\frac{1}{(s^2+\frac{\log\alpha}{N})^2}\,ds$$
that is simple to estimate in terms of $N$ and $\alpha$. The original integral is a weigthed sum of these integrals, that according to my computations should behave like 
$$\exp\left(-\log(\alpha)^{3/2}\right).$$
But I am probably over-complicating things, and we may recover the same bound by just applying a modified version of Laplace's method to the original integral.
