First, define a sequence $F_0,F_1,\dots$ of functions by $$F_0(x,z) = z,$$ $$F_k(x,z)=x\exp\left(F_{k-1}(x,z)\right) \quad \text{for }k\geq1.$$ So, for example, $$F_1(x,z) = x e^z, \quad F_2(x,z)=xe^{xe^z},\dots$$ etc. Also, set $F_{-1}(x,z)=0$. Now, let $$G(x,z) = \sum_{k=0}^\infty \left(F_k(x,z) - F_{k-1}(x,x)\right).$$ That is, $$G(x,z) = z + \left(xe^z - x\right)+\left(xe^{xe^z} - xe^x\right) + \dots$$

What I would like to do is to get some information (it doesn't have to be amazingly strong information...) about the asymptotics of the coefficient of the $x^{n-j}z^j$ term in the power series for $G(x,z)$.

Question: Does anyone know whether I have any hope in extracting any information from this generating function? If so, any ideas about what I should do/try? Even a pointer to something in the literature which might help me would be great!

By the way, the function $G(x,z)$ is closely linked to the Lambert $W$ Function. In particular, (I think) it is not hard to see that $$G(x,x)=\sum_{n=1}^\infty \frac{n^{n-1} x^n}{n!}$$ and it is well known that this function is the solution to the functional equation $$G(x,x) = x\exp(G(x,x)).$$ The thing that makes this question tricky therefore seems to be the presence of the second variable, $z$.