Theta diffusion stability criterion ($G(s)=[1-(1-a)s]/[1+as],s\ge 0$)

On an orthogonal diffusion eigenmode, the theta method has the displayed amplification with $s=-\Delta t\lambda$. The identity $(1+as{)}^2-[1-(1-a)s{]}^2=2s+(2a-1)s^2$ proves contraction at all nonnegative $s$ exactly when $a\ge 1/2$. For $a<1/2$, contraction holds up to $s\le 2/(1-2a)$, with a positive denominator in that range. Orthogonality is what transfers these modal bounds to a full L2 norm estimate.

At $a=1/2$, arbitrarily stiff heat modes approach amplification $-1$: the Crank-Nicolson method is contractive but does not rapidly damp unresolved rough components. The Backward Euler method instead has limit zero. Thus unconditional linear stability and strong stiff decay are different requirements. High smooth-data consistency order alone does not imply a uniform high-order estimate down to the initial instant for arbitrary L2 initial conditions.