Question.
Let $(z_n)$ be a complex sequence such that $$ z_0 = 0,\quad z_1=1, \quad z_{n}=\frac13(z_{n-1}+z_{n-2}+z_{n-1}z_{n-2})\quad (n>1)\tag{1} $$ Prove that the limit of the sequence exists.
Observations.
This is a complex sequence; one cannot expect to apply the monotone convergence theorem directly.(Added notes: Thanks to paw88789's comment below, this is actually a real sequence. The $z_n$'s are originally coefficients of a complex power series (in the linked question mentioned later) which could be complex. But I overlooked the assumption that it is immediate from the definition that they are actually all real. One could generalize this problem by replacing $z_1$ with some non-real complex number though. But that's not the current focus. Since this sequence is actually real and bounded from below, one can directly apply the monotone convergence theorem if one can show that the sequence is "eventually" decreasing.)
One straightforward idea is to show that the sequence is Cauchy. However, it seems difficult to directly estimate $|z_m-z_n|$ using (1).
The recurrence in (1) is very close to linear, where one can use the characteristic function techniques. But it is not very clear how to pass from "linear" to "nonlinear".
Remarks.
This is a follow-up question to another one of mine. The question above is an essential step in one approach to establish the calculation in the linked question. I find this isolated step itself is interesting. Since I figured out a detailed proof enlighted by a sketchy comment under the linked post, I will post one proof as an answer below.
I would also like to see other approaches to the question.