Per-eigenvector Davis sin 2θ and off-diagonal tan 2θ bounds #
Let T and H be symmetric operators on a real or complex inner-product space. Suppose
U is invariant under T, and T has two separated quadratic-form bounds: at least
b * ‖u‖² on U and at most a * ‖z‖² on Uᗮ. For a unit eigenvector x of T + H,
write x = y + z with y ∈ U and z ∈ Uᗮ.
The sharp per-eigenvector double-angle product estimate is
(b-a) * ‖y‖ * ‖z‖ ≤ ε when ‖H v‖ ≤ ε * ‖v‖.
With an orthogonal projection onto U, θ = arccos ‖P_U x‖ gives the angle form
(b-a) * sin (2θ) ≤ 2ε. If b > a, these imply the usual divided
bound sin (2θ) ≤ 2ε / (b-a); positivity of the gap is required for division.
If the two diagonal blocks of H with respect to U ⊕ Uᗮ vanish, the sharper,
homogeneous off-diagonal product estimate is
(b-a) * ‖y‖ * ‖z‖ ≤ |‖y‖² - ‖z‖²| * ε.
For a unit vector, b > a, and cos (2θ) ≠ 0, this implies the absolute tangent
bound |tan (2θ)| ≤ 2ε / (b-a). The product estimate itself requires neither
unit normalization nor a tangent pole restriction.
These are eigenvector bounds, not the finite-dimensional unitarily invariant
subspace sin 2Θ theorem or the general sin Θ theorems.
They do not require finite-dimensionality. The proof uses a phase-free block-energy
identity and real two-vector quadratic-form identities.
Main declarations #
LinearMap.IsSymmetric.sin_two_theta_le_of_memandLinearMap.IsSymmetric.sin_two_theta_le.LinearMap.IsSymmetric.sin_two_arccos_le.LinearMap.IsSymmetric.tan_two_theta_le_of_memandLinearMap.IsSymmetric.tan_two_theta_le.
References and provenance #
- C. Davis, The rotation of eigenvectors by a perturbation, J. Math. Anal. Appl. 6 (1963), 159–173.
- C. Davis and W. M. Kahan, The rotation of eigenvectors by a perturbation. III, SIAM J. Numer. Anal. 7 (1970), 1–46.
- Adapted from
ForTauCeti/Analysis/InnerProductSpace/DoubleAngle/Vector.leanin the AIQ-Kitware DKPS formalization (Apache-2.0); attribution retained in the source header.
Davis's sin 2θ theorem, orthogonal-decomposition form. Let T, H
be symmetric, U a T-invariant subspace with the quadratic form of T at
least b * ‖·‖ ^ 2 on U and at most a * ‖·‖ ^ 2 on Uᗮ, and let
y + z (y ∈ U, z ∈ Uᗮ) be a unit eigenvector of T + H with real
eigenvalue μ — no location assumption on μ. If ‖H v‖ ≤ ε * ‖v‖ for
all v, then
(b - a) * (‖y‖ * ‖z‖) ≤ ε.
For the angle θ of the unit eigenvector with U, 2 * ‖y‖ * ‖z‖ = sin (2θ).
Only when b > a can the product inequality be divided by the spectral gap to
obtain sin (2θ) ≤ 2ε / (b-a). The product statement also holds for b ≤ a and
does not require an orthogonal projection onto U.
Davis's sin 2θ theorem, per-eigenvector product form. Under the
hypotheses of sin_two_theta_le_of_mem, for a unit eigenvector x of
T + H and P = U.starProjection,
(b - a) * (‖P x‖ * ‖x - P x‖) ≤ ε,
When b > a, this implies the divided sin (2θ) ≤ 2ε / (b-a) bound for
the angle θ between x and U.
Davis's sin 2θ theorem, angle form. Under the hypotheses of
sin_two_theta_le, with θ = arccos ‖P x‖ the angle between the unit
eigenvector x and the invariant subspace U,
(b - a) * sin (2θ) ≤ 2 * ε, with no division by b-a.
Off-diagonal double-angle product bound. Suppose the diagonal
blocks of H relative to U ⊕ Uᗮ vanish. For any (possibly nonunit) eigenvector
y+z of T+H, with y ∈ U and z ∈ Uᗮ,
(b-a) * ‖y‖ * ‖z‖ ≤ |‖y‖² - ‖z‖²| * ε.
The result is homogeneous and has no positivity-of-gap or transversality side conditions.
For unit eigenvectors, with b>a and cos (2θ) ≠ 0, it yields
|tan (2θ)| ≤ 2ε/(b-a); the absolute value and the pole exclusion are essential.
The proof uses the same block-energy identity as the sine estimate.
Off-diagonal double-angle product bound, projection form.
Under the hypotheses of tan_two_theta_le_of_mem, for any eigenvector x of
T+H and P=U.starProjection,
(b-a) * ‖P x‖ * ‖x-P x‖ ≤ |‖P x‖² - ‖x-P x‖²| * ε.
The statement is homogeneous: ‖x‖=1 is not required. For unit x, b>a,
and away from cos (2θ)=0, this gives the absolute tan (2θ) bound.