Consider a Toeplitz matrix $T$, indexed by $\mathbb{N}_0 \times \mathbb{N}_0$. given by the sequence $t_k,k \in \mathbb{Z}$ where $t_k \geq 0,\sum_{k=-\infty}^\infty t_k=1$. By this I mean that $T_{i,i+k}=t_k$ for all $i \in \mathbb{N}_0$ and $k \in \{ -i,-(i-1),\dots,0,1,\dots \}$.
It is easy to see that $T$ defines a bounded operator on all the $\ell^p(\mathbb{N}_0)$ spaces. Suppose additionally that there are vectors $b,x^*$ with $Tx^*+b=x^*$.
I would like to understand under what assumptions the iteration $x^{(k+1)}= Tx^{(k)}+b$ will converge to $x^*$.
It is not hard to see that this will not happen if $\lim_{i \to \infty} x^{(0)}_i$ and $\lim_{i \to \infty} x^*_i$ both exist but are different. Thus $\ell^\infty$ perturbations will in general not decay. But it seems that this should be the only obstacle, i.e. $\ell^p$ perturbations should decay.
However, Szego's theorem tells us that the spectrum of $T$ accumulates around $1$, since it traces out the curve $\left \{ \sum_{k=-\infty}^\infty t_k e^{2 \pi i k \lambda} : \lambda \in [0,1] \right \}$. So it seems that there is no hope of using a simple Banach fixed point theorem argument. Is there some workaround? Perhaps some kind of Riemann-Lebesgue lemma type argument based on the idea that the eigenvectors whose eigenvalues are near $1$ cannot contribute too much to any fixed $x^{(0)}$? Or perhaps some kind of Perron-Frobenius type result?