From the answer of this question, it has been shown that the infinite integral $$ f_v(x,y) = \int_0^\infty e^{-xu-y \sqrt{u^2+v^2}} \, \mathrm{d} u $$ can conveniently be expressed in terms of infinite series involving modified Bessel functions of the second kind of the form $$ f_v(x,y) = \sum_{n = 0}^\infty v \left( \frac{x^2 v}{y} \right)^n \left( \frac{2^{-n}}{n!}\, K_{n+1}(vy) - \frac{x}{2n+1} \frac{2^n n!}{(2n)!} \left( \frac{2v}{\pi y} \right)^\frac{1}{2} \, K_{n+\frac{3}{2}} (vy) \right) , $$ where $x \in \mathbb{C}$ and $y,v \in \mathbb{R}_+^*$. In particular, for $x=0$, it follows readily that $$ f_v (0,y) = \int_0^\infty e^{-y \sqrt{u^2+v^2}} \, \mathrm{d}u = v K_1(vy) \, . $$
I would like to know under which circumstances the series converges to the above function. For instance, for $x=3/2$, $y=4/5$, and $v=1/2$, it can be checked that the series is diverging.
What I know is the asymptotic expansion of the factorial $n! \sim \sqrt{2\pi n} (n/e)^n$ as $n \to \infty$. However, modified Bessel functions of the second kind do not have asymptotic expansions as $n \to \infty$ apparently.
Any thoughts on how to proceed with this?