I am currently interested in the following sequence: $$\begin{cases}u_0 & = & \alpha\\u_n & = & u_{n-1}^2-n\end{cases}$$ where $\alpha > C \approx 1.75793275... $ with $C$ being the Nested Radical Constant.

I'm interested to get the behaviour of this sequence as $n$ grows to infinity.

Furthermore, even though I would love to have an answer for any $\alpha$, you can for the sake of simplicity simply consider the case where $\alpha = 2$. In what follows, I consider only this case.

The OEIS know little about this sequence (http://oeis.org/A198959), and is not really helpful.

Now, I know (meaning I proved) very few about this sequence. It is obviously non-negative, and non-decreasing. Plus it goes to infinity.

Then, I feel like we should have $u_n \sim \lambda^{2^n}$ for some $1 < \lambda \leqslant \alpha$. The good news is that, when I compute the first hundred of terms, this equivalent seems pretty correct, with (in the case $\alpha=2$) $\lambda \approx 1.613590596957970...$.

Unfortunately, I am stuck here. I have been unable to prove this equivalent, or to find the following terms in the expansion of this sequence. And Plouffe's inverter (now Inverse Symbolic Calculator) does not recognize this number.

Do any of you have any ideas, clues, proofs, or links to papers related to this sequence? As it is just a matter of mathematical curiosity, I would take any relevant advance (a proof of the equivalent, a way to get the following terms, a close value for $\lambda$ (or an implicit formula, or involving $C$ itself, or an algorithm, etc.), a solution for the general problem with $\alpha \neq 2$, etc.).

Of course, as I did not proved the equivalent, it may be just false, in that case a counterargument would be welcomed too :-).

Thanks all.