False under the hypotheses as printed. Nonzero coefficients do not ensure the two-column dominant-subspace condition.
For a counterexample, take , eigenvalues
and an orthonormal eigenbasis for which
The two displayed coefficient vectors are orthonormal and can be completed to an orthogonal matrix, so such a real symmetric matrix exists. Every and is nonzero, but
Indeed, with ,
is a nonzero multiple of , while the remaining direction in
is
Thus
By the simultaneous iteration interpretation of the QR algorithm, the leading two columns of the accumulated factor span this same space. The leading block is therefore an orthogonal-coordinate representation of the compression of to this space. Since
we have
for every relevant , rather than . Hence the requested Hausdorff convergence does not hold.
For completeness, the intended statement becomes true if one adds
Let consist of the first two columns of . The Accumulated QR factorization identity gives
After division by , every component along
tends to zero because of the strict spectral gap, while makes the two surviving dominant components independent. Therefore these two-dimensional subspaces converge to the dominant invariant subspace
There are two-by-two orthogonal matrices such that
Because
up to the harmless index convention at , we obtain
in the matrix 2-norm. Orthogonal similarity preserves the spectrum, and the stated spectral perturbation bound now gives