I'm specifically interested in estimating the growth of the coefficients $(a_n)_{n\in\mathbb{N}}$ of the Taylor series of $$\exp(1/(1-z))$$ centered in $0$: knowing that $\limsup_{n\rightarrow+\infty}|a_n|^{1/n}=1$, what I'm trying to figure out is whether the growth of $(a_n)_{n\in\mathbb{N}}$ is at most polynomial or super-polynomial... I tried to calculate explicitly the coefficients in order to find some regularity in the sequence of the coefficients, but I quickly got lost in the process.
Thanks in advance for any answer or suggestion.
Here's a quick proof of the leading asymptotic term. Use the associated Laguerre polynomial generating function (Gradshteyn 8.975.1) to determine that $$e^{-1}e^{1/(1-z)} = e^{z/(1-z)} = \sum_{n=0}^\infty L_n^{-1}(-1)z^n.$$ Use a contour integral representation of the Laguerre polynomials to show that $$L_n^{-1}(-1)=\frac{1}{n\,2\pi\,i} \oint (1+1/t)^n\,e^t\,dt.$$ The contour will be deformed through the saddle point which occurs at $t_0=\tfrac{1}{2}(-1+\sqrt{4n+1}).$ Saddle point analysis yields the approximation $$L_n^{-1}(-1) \sim \exp{(2\sqrt{n} -1/2)}\, \frac{1}{n\,2\pi\,i}\lim_{a \to \infty}\int_{-ia}^{ia} \exp{(c(t-t_o)^2)}\,dt.$$ In this formula the leading exponential is the asymptotic expansion for $n \to \infty$ of $(1+1/t_0)^n\,e^{t_0}.$ The $c$ in the equation, again keeping only the largest term, is $c=1/\sqrt{n}.$A rigorous application of saddle point method will show that value of all but that on a vertical line segment can be neglected, and that segment can be extended to $\infty$ with negligible affect, and the integral becomes the well-known Gaussian. Collecting results we find that $$L_n^{-1}(-1) \sim \frac{\exp{(2\sqrt{n} -1/2)}}{ 2\sqrt{\pi} \, n^{3/4} } $$ It is interesting that a naive first attempt argued as follows: $$L_n^{-1}(-1) = \frac{1}{n} \sum_{m=1}^n \binom{n}{m} \frac{m}{m!} \sim \sum_{m=1}^n \frac{n^m}{m!} \frac{m}{m!} \sim \sum_{m=1}^\infty \frac{n^m}{m!} \frac{m}{m!}=\frac{1}{\sqrt{n}}\,I_1(2\sqrt{n}).$$ One can use the known asymptotics of the Bessel function and get an expansion that agrees with the above $\textit{except for a factor of }$ $\exp{(-1/2)}.$ I thought that $\binom{n}{m} \sim n^m/m!$ would certainly work since the weight of $1/(m-1)!$ heavily emphasizes the beginning terms of the sum. Numerical evidence shows that it wasn't. I should've remembered that the approximation assumes fixed $m$, but there is a summation so $m$ can get quite large. It's interesting that a mistake only ended up with a different constant before the asymptotic expression.