In Appendix A of [this paper](https://arxiv.org/pdf/cond-mat/9904278.pdf), it is claimed that the asymptotic behaviour of

$$\phi_1(y,\lambda)=\frac{1}{\Gamma(\frac{1-\lambda}{2})}\int_0^\infty dt~e^{-t}\cos(2y\sqrt{t})t^{-\frac{1+\lambda}{2}},$$

for $\lambda_2=\text{Im}(\lambda)\to+\infty$ (where $\lambda=\lambda_1+i\lambda_2$), is

$$\phi_1(y,\lambda_1+i\lambda_2)\sim\cosh[y\sqrt{\lambda_2-i\lambda_1}(1+i)]\exp(-y^2/2) .$$

I have been trying various asymptotic methods to arrive at this result, such as the method of stationary phase, without success. How might one derive this result?

<!---->