Documentation

TauCeti.Analysis.ODE.EnergyEstimate.Index

The index of the trajectory operator #

Let A : ℝ → End(E) be a continuous family of operators on a finite-dimensional real inner product space E, with limits A± at ±∞ that have Lyapunov forms (ContinuousLinearMap.HasLyapunovForm), so that the trajectory operator D_A u = u' + A u : W^{1,2}(ℝ, E) → L²(ℝ, E) is Fredholm (TauCeti.WeakDerivStep.isFredholm_trajectoryOperator_of_hasLyapunovForm). Invertible symmetric limits have Lyapunov forms (ContinuousLinearMap.hasLyapunovForm_of_isSymmetric). This file proves:

For the first statement, the affine family (1 - t) A + t B has the same limits as A and B for every real t, so each of its trajectory operators is Fredholm. The trajectory operator depends continuously on its L∞ coefficient (TauCeti.WeakDerivStep.lipschitzWith_trajectoryOperator), so the index is constant along this family. For the second, a constant coefficient A₀ with a spectral gap gives an injective D_{A₀} and an injective formal adjoint D_{-A₀}, by the coercive estimate TauCeti.WeakDerivStep.ker_trajectoryOperator_eq_bot, so D_{A₀} is bijective and has index zero; a family with equal limits A₀ has the index of the constant family.

The first statement is the invariance of the index under deformations of the coefficient family that fix its limits; the second is its value when the limits agree. Together they reduce the computation of the index for arbitrary limits to that of one family for each pair of limits. The third statement is the form that computation takes: the kernel of D_A is the space of square-integrable solutions of v' + A v = 0 (TauCeti.WeakDerivStep.kerTrajectoryOperatorEquiv), and the cokernel is the kernel of the formal adjoint D_{-A}.

For the fourth statement we use the model family ψ(-s) A⁻ + ψ(s) A⁺, where ψ is Mathlib's Real.smoothTransition. On [0, ∞) it is ψ(s) A⁺, a scalar multiple of one symmetric operator, so a solution of v' = -ψ(s) A⁺ v multiplies the coordinate of v along an eigenvector of A⁺ with eigenvalue λ by exp (-λ ∫₀ˢ ψ). It is square integrable on (0, ∞) exactly when v 0 lies in the positive spectral subspace W of A⁺. Likewise it is square integrable on (-∞, 0) exactly when v 0 lies in the negative spectral subspace U of A⁻. So the kernel of the model operator is W ⊓ U, and the same argument for -A identifies its cokernel with Wᗮ ⊓ Uᗮ, because the positive and negative spectral subspaces of an invertible symmetric operator are orthogonal complements. Hence the index is dim (W ⊓ U) - dim (W ⊔ U)ᗮ = dim W + dim U - dim E, which is n⁻(A⁻) - n⁻(A⁺).

For the fifth statement, the index is also constant along (1 - t) A + t B when the limits move, as long as (1 - t) A± + t B± keeps a Lyapunov form for t ∈ [0, 1]. A symmetric Lyapunov form S of A is a Lyapunov form of every (1 - t) A + t S with 0 ≤ t ≤ 1 (ContinuousLinearMap.hasLyapunovForm_one_sub_smul_add_smul), so the limits A± can be deformed to the invertible symmetric operators S±, where the fourth statement applies. For a constant family, computing the index with the forms S and S' at the two ends shows n⁻(S) = n⁻(S').

Main declarations #

References #

theorem TauCeti.WeakDerivStep.index_trajectoryOperator_eq_of_forall_hasLyapunovForm {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {A B : ℝ → E →L[ℝ] E} (hA : Continuous A) (hB : Continuous B) {Aminus Aplus Bminus Bplus : E →L[ℝ] E} (hAbot : Filter.Tendsto A Filter.atBot (nhds Aminus)) (hAtop : Filter.Tendsto A Filter.atTop (nhds Aplus)) (hBbot : Filter.Tendsto B Filter.atBot (nhds Bminus)) (hBtop : Filter.Tendsto B Filter.atTop (nhds Bplus)) (AL BL : ↥(MeasureTheory.Lp (E →L[ℝ] E) ⊤ (MeasureTheory.volume.restrict ↑⊤))) (hAL : ↑↑AL =ᵐ[MeasureTheory.volume.restrict ↑⊤] A) (hBL : ↑↑BL =ᵐ[MeasureTheory.volume.restrict ↑⊤] B) (hminus : ∀ t ∈ Set.Icc 0 1, ((1 - t) • Aminus + t • Bminus).HasLyapunovForm) (hplus : ∀ t ∈ Set.Icc 0 1, ((1 - t) • Aplus + t • Bplus).HasLyapunovForm) :

Homotopy invariance of the index through hyperbolic limits. Let A and B be continuous families of operators on a finite-dimensional real inner product space with limits A± and B± at ±∞, such that (1 - t) A± + t B± has a Lyapunov form for every t ∈ [0, 1]. Then the trajectory operators u ↦ u' + A u and u ↦ u' + B u from W^{1,2}(ℝ, E) to L²(ℝ, E) have the same Fredholm index. Unlike TauCeti.WeakDerivStep.index_trajectoryOperator_eq_of_tendsto, the limits may move.

theorem TauCeti.WeakDerivStep.index_trajectoryOperator_eq_of_tendsto {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {A B : ℝ → E →L[ℝ] E} (hA : Continuous A) (hB : Continuous B) {Aminus Aplus : E →L[ℝ] E} (hminus : Aminus.HasLyapunovForm) (hplus : Aplus.HasLyapunovForm) (hAbot : Filter.Tendsto A Filter.atBot (nhds Aminus)) (hAtop : Filter.Tendsto A Filter.atTop (nhds Aplus)) (hBbot : Filter.Tendsto B Filter.atBot (nhds Aminus)) (hBtop : Filter.Tendsto B Filter.atTop (nhds Aplus)) (AL BL : ↥(MeasureTheory.Lp (E →L[ℝ] E) ⊤ (MeasureTheory.volume.restrict ↑⊤))) (hAL : ↑↑AL =ᵐ[MeasureTheory.volume.restrict ↑⊤] A) (hBL : ↑↑BL =ᵐ[MeasureTheory.volume.restrict ↑⊤] B) :

The index of the trajectory operator depends only on the limits. Let A and B be continuous families of operators on a finite-dimensional real inner product space with the same limits A± at ±∞, which have Lyapunov forms. Then the trajectory operators u ↦ u' + A u and u ↦ u' + B u from W^{1,2}(ℝ, E) to L²(ℝ, E) have the same Fredholm index.

theorem TauCeti.WeakDerivStep.bijective_trajectoryOperator_const {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {A₀ : E →L[ℝ] E} (hsym : (↑A₀).IsSymmetric) {α : ℝ} (hα : 0 < α) (hgap : ∀ (v : E), α * ‖v‖ ≤ ‖A₀ v‖) (AL : ↥(MeasureTheory.Lp (E →L[ℝ] E) ⊤ (MeasureTheory.volume.restrict ↑⊤))) (hAL : ↑↑AL =ᵐ[MeasureTheory.volume.restrict ↑⊤] fun (x : ℝ) => A₀) :

A constant invertible coefficient gives a bijective trajectory operator. For a symmetric operator A₀ on a finite-dimensional real inner product space with α ‖v‖ ≤ ‖A₀ v‖ for some α > 0, the trajectory operator u ↦ u' + A₀ u from W^{1,2}(ℝ, E) to L²(ℝ, E) is bijective. Both it and its formal adjoint -d/ds + A₀ are injective, and its range is closed.

theorem TauCeti.WeakDerivStep.index_trajectoryOperator_eq_zero_of_tendsto {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {A : ℝ → E →L[ℝ] E} (hA : Continuous A) {A₀ : E →L[ℝ] E} (hsym₀ : (↑A₀).IsSymmetric) {α : ℝ} (hα : 0 < α) (hgap : ∀ (v : E), α * ‖v‖ ≤ ‖A₀ v‖) (hbot : Filter.Tendsto A Filter.atBot (nhds A₀)) (htop : Filter.Tendsto A Filter.atTop (nhds A₀)) (AL : ↥(MeasureTheory.Lp (E →L[ℝ] E) ⊤ (MeasureTheory.volume.restrict ↑⊤))) (hAL : ↑↑AL =ᵐ[MeasureTheory.volume.restrict ↑⊤] A) :

Equal limits give index zero. Let A be a continuous family of operators on a finite-dimensional real inner product space with the same invertible symmetric limit A₀ at -∞ and at +∞. Then the trajectory operator u ↦ u' + A u from W^{1,2}(ℝ, E) to L²(ℝ, E) has index zero.

The index of a trajectory operator with closed range counts square-integrable solutions. For a continuous family A of symmetric operators on a finite-dimensional real inner product space E such that the range of D_A u = u' + A u from W^{1,2}(ℝ, E) to L²(ℝ, E) is closed, the index of D_A is the dimension of the space of square-integrable solutions of v' = -A v minus the dimension of the space of square-integrable solutions of v' = A v.

theorem TauCeti.WeakDerivStep.index_trajectoryOperator_eq_finrank_sub_finrank {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {A : ℝ → E →L[ℝ] E} (hA : Continuous A) (hsym : ∀ (s : ℝ), (↑(A s)).IsSymmetric) {Aminus Aplus : E →L[ℝ] E} {α : ℝ} (hα : 0 < α) (hminus : ∀ (v : E), α * ‖v‖ ≤ ‖Aminus v‖) (hplus : ∀ (v : E), α * ‖v‖ ≤ ‖Aplus v‖) (hbot : Filter.Tendsto A Filter.atBot (nhds Aminus)) (htop : Filter.Tendsto A Filter.atTop (nhds Aplus)) (AL : ↥(MeasureTheory.Lp (E →L[ℝ] E) ⊤ (MeasureTheory.volume.restrict ↑⊤))) (hAL : ↑↑AL =ᵐ[MeasureTheory.volume.restrict ↑⊤] A) :

The index of the trajectory operator counts square-integrable solutions. For a continuous family A of symmetric operators on a finite-dimensional real inner product space with invertible limits at both ends, the index of D_A u = u' + A u from W^{1,2}(ℝ, E) to L²(ℝ, E) is the dimension of the space of square-integrable solutions of v' = -A v minus the dimension of the space of square-integrable solutions of v' = A v. This specializes TauCeti.WeakDerivStep.index_trajectoryOperator_eq_finrank_sub_finrank_of_isClosed_range.

The index as a difference of Morse indices #

theorem TauCeti.WeakDerivStep.index_trajectoryOperator_eq_finrank_pointSpectralSubspace_sub_finrank {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {A : ℝ → E →L[ℝ] E} (hA : Continuous A) {Aminus Aplus : E →L[ℝ] E} (hsymminus : (↑Aminus).IsSymmetric) (hsymplus : (↑Aplus).IsSymmetric) {α : ℝ} (hα : 0 < α) (hminus : ∀ (v : E), α * ‖v‖ ≤ ‖Aminus v‖) (hplus : ∀ (v : E), α * ‖v‖ ≤ ‖Aplus v‖) (hbot : Filter.Tendsto A Filter.atBot (nhds Aminus)) (htop : Filter.Tendsto A Filter.atTop (nhds Aplus)) (AL : ↥(MeasureTheory.Lp (E →L[ℝ] E) ⊤ (MeasureTheory.volume.restrict ↑⊤))) (hAL : ↑↑AL =ᵐ[MeasureTheory.volume.restrict ↑⊤] A) :

The index of the trajectory operator is a difference of Morse indices. Let A be a continuous family of operators on a finite-dimensional real inner product space with invertible symmetric limits A± at ±∞. Then the index of the trajectory operator u ↦ u' + A u from W^{1,2}(ℝ, E) to L²(ℝ, E) is the number of negative eigenvalues of A⁻ minus the number of negative eigenvalues of A⁺, counted with multiplicity: the dimensions of the sums of their eigenspaces at negative eigenvalues.

theorem TauCeti.WeakDerivStep.index_trajectoryOperator_eq_finrank_pointSpectralSubspace_sub_finrank_of_forall_le_inner {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {A : ℝ → E →L[ℝ] E} (hA : Continuous A) {Aminus Aplus Sminus Splus : E →L[ℝ] E} (hSminus : (↑Sminus).IsSymmetric) (hSplus : (↑Splus).IsSymmetric) {cminus cplus : ℝ} (hcminus : 0 < cminus) (hcplus : 0 < cplus) (hminus : ∀ (v : E), cminus * ‖v‖ ^ 2 ≤ inner ℝ (Sminus (Aminus v)) v) (hplus : ∀ (v : E), cplus * ‖v‖ ^ 2 ≤ inner ℝ (Splus (Aplus v)) v) (hbot : Filter.Tendsto A Filter.atBot (nhds Aminus)) (htop : Filter.Tendsto A Filter.atTop (nhds Aplus)) (AL : ↥(MeasureTheory.Lp (E →L[ℝ] E) ⊤ (MeasureTheory.volume.restrict ↑⊤))) (hAL : ↑↑AL =ᵐ[MeasureTheory.volume.restrict ↑⊤] A) :

The index of the trajectory operator with hyperbolic limits is a difference of inertia indices of Lyapunov forms. Let A be a continuous family of operators on a finite-dimensional real inner product space whose limits A± at ±∞ have symmetric Lyapunov forms S±: c± ‖v‖² ≤ ⟪S± (A± v), v⟫ with c± > 0. Then the index of the trajectory operator u ↦ u' + A u from W^{1,2}(ℝ, E) to L²(ℝ, E) is the number of negative eigenvalues of S⁻ minus the number of negative eigenvalues of S⁺, counted with multiplicity. Neither A nor its limits need to be symmetric; for invertible symmetric limits one may take S± = A±, which gives TauCeti.WeakDerivStep.index_trajectoryOperator_eq_finrank_pointSpectralSubspace_sub_finrank.

theorem ContinuousLinearMap.finrank_pointSpectralSubspace_Iio_eq_of_forall_le_inner {E : Type u_1} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E] {A S S' : E →L[ℝ] E} (hS : (↑S).IsSymmetric) (hS' : (↑S').IsSymmetric) {c c' : ℝ} (hc : 0 < c) (hc' : 0 < c') (h : ∀ (v : E), c * ‖v‖ ^ 2 ≤ inner ℝ (S (A v)) v) (h' : ∀ (v : E), c' * ‖v‖ ^ 2 ≤ inner ℝ (S' (A v)) v) :

The negative inertia of a Lyapunov form depends only on the operator. If S and S' are symmetric Lyapunov forms of the same operator A on a finite-dimensional real inner product space, c ‖v‖² ≤ ⟪S (A v), v⟫ and c' ‖v‖² ≤ ⟪S' (A v), v⟫ with c, c' > 0, then S and S' have the same number of negative eigenvalues, counted with multiplicity. By the inertia theorem of Ostrowski and Schneider, which is not formalized here, this number is the dimension of the sum of the generalized eigenspaces of A at eigenvalues of negative real part.