Self-adjoint short-range operators

Written by GPT-6.1 Sol (OpenAI) and GPT-6 Astra (OpenAI). Self-checked by the writing AI. Original exposition: CC0.

Working question: Which closed realization is being scattered? A formal expression p(D)+Vp(D)+V is not yet a self-adjoint operator. A bounded potential permits the elementary domain argument in the first lesson, but a general differential perturbation needs its own closure and factorization. Keeping that distinction prevents an unproved relative-bound assumption from entering the wave construction.

A symmetric differential perturbation can be unbounded on L2L^2. Compactness on the endpoint graph space still gives a self-adjoint closure and an exact resolvent factorization. The operator domain will be the closure of its Schwartz graph; every resolvent composition below is justified on that domain.

Let pp be real, simply characteristic, with no invariant direction, and let P0=p(D)P_0=p(D) on its maximal L2L^2 domain. Let VV be a short-range finite-order local differential expression with Lloc2L^2_{\mathrm{loc}} coefficients, in the sense of Short-range compactness and local tests. Use its graph space XpX_p, unique compact extension V:Xp→BV:X_p\to B, and distribution-to-norm convergence theorem. Assume VV is symmetric on compact smooth functions:

(Vϕ,ψ)=(ϕ,Vψ),ϕ,ψ∈Cc∞.(1) (V\phi,\psi)=(\phi,V\psi),\qquad \phi,\psi\in C_c^\infty. \tag{1}

The pairing is linear in its first argument. The compact alternative needed here and in the later boundary analysis is proved next on an arbitrary Banach space. Its finite-dimensional facts—bounded coordinates, compact unit spheres, closed finite-dimensional subspaces and complete quotients—are proved in Compact Fredholm operators, elementary tools. We use the complete shell duality, density and embeddings of Endpoint spaces and flat energy shells. The Fourier transform on Schwartz functions and tempered distributions is proved in Measure and Fourier foundations.

Compact alternative used below. If CC is compact on a Banach space YY, then A=I+CA=I+C has closed range, finite-dimensional kernel, and dim⁡ker⁡A=dim⁡(Y/ran⁡A). \dim\ker A=\dim(Y/\operatorname{ran}A). Consequently injectivity implies a bounded inverse on all of YY.

Proof. The unit ball in ker⁡A\ker A is precompact, because x=−Cxx=-Cx there. A normed space with precompact unit ball is finite dimensional: in an infinite-dimensional space, successively choose unit vectors at distance greater than 1/21/2 from the spans of the preceding vectors. The elementary selection works for any proper closed subspace MM: take v∉Mv\notin M, put d=dist⁡(v,M)>0d=\operatorname{dist}(v,M)>0, choose w∈Mw\in M with ∥v−w∥<2d\|v-w\|<2d, and normalize v−wv-w. Its distance from MM exceeds 1/21/2. Applying this to the closed finite-dimensional spans gives a sequence with no convergent subsequence.

There is a constant c>0c>0 such that ∥Ax∥≥cdist⁡(x,ker⁡A)\|Ax\|\geq c\operatorname{dist}(x,\ker A). Otherwise choose quotient classes of norm one with images tending to zero and representatives xjx_j of norm at most two. A subsequence of CxjCx_j converges by compactness; xj=Axj−Cxjx_j=Ax_j-Cx_j then converges to a vector of ker⁡A\ker A, contradicting its quotient norm one. If AxjAx_j converges, subtract kernel vectors so that these representatives are bounded, using this estimate. Compactness again gives a subsequence for which CxjCx_j, and hence xjx_j, converge. Its limit maps to the prescribed output limit. Thus the range is closed. The same arguments apply to every Ak=I+CkA^k=I+C_k, since its finite polynomial remainder CkC_k is compact. Indeed the image of a precompact set under a bounded map is precompact, and the sum of finitely many precompact sets is precompact: choose convergent subsequences successively for each summand. Thus every positive power of CC, and then the finite binomial remainder of (I+C)k(I+C)^k, is compact.

Put Nk=ker⁡AkN_k=\ker A^k and Rk=ran⁡AkR_k=\operatorname{ran}A^k, with N0=0N_0=0, R0=YR_0=Y. The increasing chain NkN_k stabilizes. If it did not, choose unit xk∈Nkx_k\in N_k at distance greater than 1/21/2 from Nk−1N_{k-1}. For l<kl<k, both AxkAx_k and CxlCx_l belong to Nk−1N_{k-1}, so ∥Cxk−Cxl∥>1/2\|Cx_k-Cx_l\|>1/2, contradicting compactness. Likewise the decreasing chain RkR_k stabilizes: if all its inclusions were strict, choose unit xk∈Rk−1x_k\in R_{k-1} at distance greater than 1/21/2 from RkR_k. For l>kl>k, both AxkAx_k and CxlCx_l belong to RkR_k, giving the same contradiction. These choices use the closedness just proved and the elementary subspace selection above. Equality in the range chain propagates by applying AA; equality in the kernel chain propagates by taking the inverse image under AA.

Choose mm beyond both stabilization indices. Then Y=Nm⊕RmY=N_m\oplus R_m. Indeed, if x=Amy∈Nmx=A^my\in N_m, then y∈N2m=Nmy\in N_{2m}=N_m and x=0x=0. For arbitrary xx, Amx∈Rm=R2mA^mx\in R_m=R_{2m}; write Amx=A2myA^mx=A^{2m}y to get x−Amy∈Nmx-A^my\in N_m. Both summands are closed and NmN_m is finite dimensional. The projection onto NmN_m is bounded because the quotient map Y→Y/RmY\to Y/R_m is one-to-one on this finite-dimensional summand: its norm on that summand is bounded below on the compact unit sphere. Its complementary projection Q:Y→RmQ:Y\to R_m is bounded too. On RmR_m, AA is injective by the zero intersection and onto by stabilization. Since ker⁡A⊂Nm\ker A\subset N_m, for x∈Rmx\in R_m and k∈ker⁡Ak\in\ker A we have x=Q(x−k)x=Q(x-k). Hence ∥x∥≤∥Q∥dist⁡(x,ker⁡A)≤∥Q∥c−1∥Ax∥. \|x\|\leq\|Q\|\operatorname{dist}(x,\ker A) \leq \|Q\|c^{-1}\|Ax\|. This proves boundedness of the inverse on RmR_m, which complements the generalized kernel NmN_m. On NmN_m, ordinary finite-dimensional rank–nullity gives equal kernel dimension and range codimension. Since AA preserves both summands and is onto RmR_m, this proves the displayed equality on YY. When ker⁡A=0\ker A=0, the quotient lower bound becomes ∥Ax∥≥c∥x∥\|Ax\|\geq c\|x\|; the equality gives surjectivity and this bound gives boundedness of the inverse. □\square

1. Extending symmetry to endpoint waves

Lemma 1.1. For all u,v∈Xpu,v\in X_p,

(Vu,v)=(u,Vv).(2) (Vu,v)=(u,Vv). \tag{2}

Both pairings are absolutely convergent.

Proof. The graph norm controls u,v∈B∗u,v\in B^*, while Vu,Vv∈BVu,Vv\in B. Shell Cauchy–Schwarz proves absolute integrability and continuity of the pairings in the graph norms. Smooth density therefore reduces the assertion to smooth u,v∈Xpu,v\in X_p.

Choose a real compact smooth χ\chi, equal to one near zero, and set uR=χ(x/R)u, vR=χ(x/R)vu_R=\chi(x/R)u,\ v_R=\chi(x/R)v. These functions are compact and smooth. Polynomial Leibniz and the uniform bounds on the derivatives of χ(x/R)\chi(x/R), R≥1R\geq1, show that uR,vRu_R,v_R are uniformly bounded in XpX_p. They converge to u,vu,v in distributions. Theorem 4.1 of Short-range compactness and local tests gives VuR→Vu, VvR→VvVu_R\to Vu,\ Vv_R\to Vv in BB.

Apply (1) to uR,vRu_R,v_R. For example,

∣(VuR−Vu,vR)∣≤∥VuR−Vu∥B∥vR∥B∗⟶0. |(Vu_R-Vu,v_R)| \leq\|Vu_R-Vu\|_B\|v_R\|_{B^*}\longrightarrow0.

Also (Vu,vR)→(Vu,v)(Vu,v_R)\to(Vu,v) by dominated convergence, since ∣Vu v∣|Vu\,v| is integrable and the cutoffs are uniformly bounded. The other side has the same two steps. This proves (2) for smooth inputs, then for all inputs by graph-norm density. □\square

2. Essential self-adjointness

Initially let

P=P0+V,D(P)=S(Rn). P=P_0+V,\qquad \mathcal D(P)=\mathcal S(\mathbb R^n).

Schwartz functions belong to XpX_p, so VV maps this domain to B⊂L2B\subset L^2. The maximal multiplier P0P_0 is self-adjoint by Resolvents, domains and spectral density.

Theorem 2.1. PP is essentially self-adjoint. Write H=P‾H=\overline P for its self-adjoint closure.

Proof. Lemma 1.1 and the real Fourier symbol make PP symmetric on its dense domain. It is closable: if uj→0u_j\to0 and Puj→hPu_j\to h in L2L^2, then (h,v)=lim⁡(uj,Pv)=0(h,v)=\lim(u_j,Pv)=0 for all vv in the dense domain, hence h=0h=0. Its closure is closed and symmetric.

Corollary 4.3 of Short-range compactness and local tests gives

∥VRp(it)∥B→B⟶0(∣t∣→∞).(3) \|VR_p(it)\|_{B\to B}\longrightarrow0 \quad(|t|\to\infty). \tag{3}

Choose one T>0T>0 sufficiently large that this norm is less than one at both t=Tt=T and t=−Tt=-T. The convergent geometric series makes I+VRp(±iT)I+VR_p(\pm iT) invertible on BB. Explicitly, for a bounded CC with ∥C∥<1\|C\|<1, completeness of BB defines Sf=∑k≥0(−C)kfS f=\sum_{k\ge0}(-C)^kf, with ∥S∥≤(1−∥C∥)−1\|S\|\le(1-\|C\|)^{-1}. The geometric tail bound gives convergence in operator norm; multiplying its finite sums on either side by I+CI+C leaves a remainder of norm at most ∥C∥N+1\|C\|^{N+1}. Passing to the limit proves both inverse identities.

For nonreal zz, P0−zP_0-z is a bijection on S\mathcal S. Indeed its inverse is Fourier multiplication by 1/(p−z)1/(p-z); the denominator is bounded below by ∣Im⁡z∣|\operatorname{Im}z|, and every derivative of the reciprocal has at most polynomial growth. It therefore maps Schwartz functions to Schwartz functions. On that domain the exact factorization is

(P−z)u=(I+VRp(z))(P0−z)u.(4) (P-z)u=(I+VR_p(z))(P_0-z)u. \tag{4}

Consequently, for the chosen z=±iTz=\pm iT, Ran⁡(P∓iT)=(I+VRp(±iT))S\operatorname{Ran}(P\mp iT)=(I+VR_p(\pm iT))\mathcal S. Schwartz functions are dense in BB, and an invertible bounded map preserves density. Both these ranges are therefore dense in BB, hence in L2L^2.

To complete the operator argument explicitly, let A=P‾A=\overline P. Symmetry gives

∥(A−it)u∥22=∥Au∥22+t2∥u∥22,t∈R.(5) \|(A-it)u\|_2^2=\|Au\|_2^2+t^2\|u\|_2^2,\qquad t\in\mathbb R. \tag{5}

For t≠0t\ne0, this estimate and closedness show that Ran⁡(A−it)\operatorname{Ran}(A-it) is closed: a Cauchy output sequence makes both its inputs and their AA images Cauchy. Since it contains the dense initial range, it equals L2L^2. Thus A−itA-it is onto for both chosen signs.

The orthogonal-complement identity Ran⁡(A−it)⊥=ker⁡(A∗+it)\operatorname{Ran}(A-it)^\perp=\ker(A^*+it) shows that the two corresponding adjoint defect kernels vanish. If v∈D(A∗)v\in\mathcal D(A^*), use the onto range at itit to choose w∈D(A)w\in\mathcal D(A) with (A−it)w=(A∗−it)v(A-it)w=(A^*-it)v. Then v−w∈ker⁡(A∗−it)=0v-w\in\ker(A^*-it)=0, using the onto range at the opposite sign. Hence v=w∈D(A)v=w\in\mathcal D(A). The reverse inclusion follows from symmetry, so A=A∗A=A^*. □\square

The theorem identifies the operator as the closure of its Schwartz restriction. For a general differential VV, an equality of its domain with the free maximal domain would need an additional argument.

3. Justifying the resolvent factorization on BB

For Im⁡z≠0\operatorname{Im}z\ne0, write RH(z)=(H−z)−1R_H(z)=(H-z)^{-1}. This inverse exists on all of L2L^2: if z=a+ibz=a+ib, symmetry and self-adjointness give

∥(H−z)u∥22=∥(H−a)u∥22+b2∥u∥22. \|(H-z)u\|_2^2=\|(H-a)u\|_2^2+b^2\|u\|_2^2.

The same closed-range argument as above applies to H−aH-a. The orthogonal complement of the range is ker⁡(H−z‾)=0\ker(H-\overline z)=0, by the same estimate, so the range is also dense and hence equals L2L^2. Consequently ∥RH(z)∥L2→L2≤∣b∣−1\|R_H(z)\|_{L^2\to L^2}\le |b|^{-1}, and its range lies in D(H)\mathcal D(H).

Lemma 3.1. If f∈Bf\in B, then Rp(z)f∈D(H)R_p(z)f\in\mathcal D(H) and

(H−z)Rp(z)f=f+VRp(z)f.(6) (H-z)R_p(z)f=f+VR_p(z)f. \tag{6}

Proof. Choose fj∈Sf_j\in\mathcal S converging to ff in BB, hence in L2L^2. Put uj=Rp(z)fj∈Su_j=R_p(z)f_j\in\mathcal S. The ordinary resolvent bound gives uj→Rp(z)fu_j\to R_p(z)f in L2L^2. Theorem 1.1 of Global polynomial resolvent estimates, applied to each derivative of the polynomial, gives uj→Rp(z)fu_j\to R_p(z)f in XpX_p, and therefore Vuj→VRp(z)fVu_j\to VR_p(z)f in B⊂L2B\subset L^2. Also P0uj=fj+zuj⟶f+zRp(z)f P_0u_j=f_j+zu_j\longrightarrow f+zR_p(z)f in L2L^2. Thus PujPu_j converges to the sum of these two limits. Closedness of HH, as the closure of this restriction, proves membership and (6). □\square

Theorem 3.2. For every nonreal zz, I+VRp(z)I+VR_p(z) is invertible on BB, and

RH(z)f=Rp(z)(I+VRp(z))−1f∈Xp,f∈B.(7) R_H(z)f=R_p(z)(I+VR_p(z))^{-1}f\in X_p,\qquad f\in B. \tag{7}

Both resolvent identities are valid on BB:

Rp(z)=RH(z)+RH(z)VRp(z),(8) R_p(z)=R_H(z)+R_H(z)VR_p(z), \tag{8} Rp(z)=RH(z)+Rp(z)VRH(z).(9) R_p(z)=R_H(z)+R_p(z)VR_H(z). \tag{9}

Proof. The free resolvent maps BB boundedly to XpX_p, and V:Xp→BV:X_p\to B is compact, so VRp(z)VR_p(z) is compact on BB. If (I+VRp(z))f=0(I+VR_p(z))f=0, Lemma 3.1 gives (H−z)Rp(z)f=0(H-z)R_p(z)f=0. The nonreal resolvent of the self-adjoint HH is injective, so Rp(z)f=0R_p(z)f=0, then f=0f=0. Apply the compact alternative proved above to get the bounded inverse on BB.

Take g=(I+VRp(z))−1fg=(I+VR_p(z))^{-1}f. Equation (6) gives (H−z)Rp(z)g=f(H-z)R_p(z)g=f, which proves (7). In particular every occurrence of VRH(z)fVR_H(z)f in (9) is defined in BB.

Applying RH(z)R_H(z) to (6) gives (8) directly. Alternatively, (7) gives f=g+VRp(z)gf=g+VR_p(z)g and RH(z)f=Rp(z)gR_H(z)f=R_p(z)g; applying Rp(z)R_p(z) to the first equality gives (9). All compositions have now been assigned their domains and spaces. □\square

No real-boundary invertibility is asserted here. At real energies a nonzero kernel may occur; the limiting-absorption and point-spectrum analysis must handle it separately.

4. Local continuity and a concrete potential

Corollary 4.1. On each open half-plane, RH(z)R_H(z) is locally bounded as an operator B→XpB\to X_p and is weak-star continuous in every graph component.

Proof. The free map Rp(z):B→XpR_p(z):B\to X_p is locally bounded and has weak-star continuous components. The short-range convergence theorem upgrades this to norm continuity of VRp(z)VR_p(z) on BB in operator norm. To check the uniform assertion, if it failed at zz, choose zj→zz_j\to z and ∥fj∥B≤1\|f_j\|_B\leq1 for which ∥V(Rp(zj)−Rp(z))fj∥B\|V(R_p(z_j)-R_p(z))f_j\|_B stays positive. The difference is uniformly bounded in XpX_p by the compact energy bound. Its L2L^2 norm tends to zero, because the Hilbert resolvent identity gives ∥(Rp(zj)−Rp(z))fj∥2≤∣zj−z∣∣Im⁡zj∣ ∣Im⁡z∣∥fj∥2. \|(R_p(z_j)-R_p(z))f_j\|_2 \leq\frac{|z_j-z|}{|\operatorname{Im}z_j|\,|\operatorname{Im}z|} \|f_j\|_2. The convergence theorem contradicts that positive lower bound.

The identity T(z′)−T(z)=T(z′) [VRp(z)−VRp(z′)] T(z),T(z)=(I+VRp(z))−1, T(z')-T(z)=T(z')\,[VR_p(z)-VR_p(z')]\,T(z), \qquad T(z)=(I+VR_p(z))^{-1}, and a geometric series near each zz give local norm continuity and local boundedness of TT. For a fixed f∈Bf\in B, subtract (7) at the two parameters:

RH(z′)f−RH(z)f=Rp(z′)[T(z′)−T(z)]f+[Rp(z′)−Rp(z)]T(z)f. \begin{aligned} R_H(z')f-R_H(z)f &=R_p(z')[T(z')-T(z)]f\\ &\quad+[R_p(z')-R_p(z)]T(z)f. \end{aligned}

The first term tends to zero in XpX_p, by the local bound for the free resolvent and norm continuity of TT. The second tends weak-star to zero in each graph component, by the free endpoint theorem with its fixed forcing T(z)f∈BT(z)f\in B. Local boundedness follows directly from the product of the two local operator bounds. This proves the claim. □\square

Example 4.2. Let p(ξ)=ξ1−ξ22p(\xi)=\xi_1-\xi_2^2, and let VV be a real measurable function satisfying ∣V(x)∣≤C(1+∣x∣)−1−δ,δ>0. |V(x)|\leq C(1+|x|)^{-1-\delta},\qquad\delta>0. The hypotheses on the free polynomial can be checked directly:

p~(ξ)2=p(ξ)2+4ξ22+5,∣∇p(ξ)∣2=1+4ξ22. \widetilde p(\xi)^2=p(\xi)^2+4\xi_2^2+5, \qquad |\nabla p(\xi)|^2=1+4\xi_2^2.

Thus p~≤5 (∣p∣+∣∇p∣)\widetilde p\le\sqrt5\,(|p|+|\nabla p|), so pp is simply characteristic. If a vector ww were invariant, comparing p(ξ+tw)=p(ξ)p(\xi+tw)=p(\xi) for all ξ,t\xi,t gives w2=0w_2=0 from the coefficient of tξ2t\xi_2, and then w1=0w_1=0. Its invariant space is therefore zero. Strength properness and the decaying-coefficient criterion make this multiplication short range. Its real values give (1). The theorem supplies a self-adjoint closure and (7), even though the free operator is nonelliptic and is not lower bounded. In this bounded-potential example, the bounded-perturbation theorem makes A=P0+VA=P_0+V self-adjoint on the free maximal domain. It extends the initial Schwartz operator, so its closed graph contains that of HH, giving H⊂AH\subset A. Taking adjoints reverses this inclusion: directly from the adjoint definition, A∗⊂H∗A^*\subset H^*. Since both operators are self-adjoint, A⊂HA\subset H as well. They are equal, proving the claimed domain identification for this example.

Use the conclusion

Verify the extended symmetry pairing before using the deficiency argument. In the resolvent factorization, identify the space on which the compact factor is inverted and the domain to which the resulting solution belongs.

5. Exercises

Exercise 5.1 (foundation). In the scalar model P0=aP_0=a, V=bV=b, verify (7)-(9) with a,b∈Ra,b\in\mathbb R. Use the calculation to check the plus sign in I+VRp(z)I+VR_p(z).

Exercise 5.2 (foundation). Prove that a closed symmetric operator AA has closed range for A−itA-it, t≠0t\ne0, using (5). Identify where closedness, rather than symmetry alone, is used.

Exercise 5.3 (intermediate). Verify that multiplication by 1/(p−z)1/(p-z) maps S\mathcal S to itself for a real polynomial and nonreal zz. Explain why the same global Schwartz assertion cannot simply be used at a real regular energy.

Exercise 5.4 (intermediate). In Lemma 1.1's cutoff limit, prove separately the norm-error estimate and the dominated-convergence step for (VuR,vR)(Vu_R,v_R). Explain why a distributional limit alone would not justify their pairing limit.

Exercise 5.5 (advanced). Let C(z)=VRp(z)C(z)=VR_p(z) and suppose I+C(z0)I+C(z_0) is invertible on BB. Derive a sufficient norm-neighborhood condition for I+C(z)I+C(z) to remain invertible. Give an explicit bound for its inverse and derive the inverse-difference formula used in Corollary 4.1.

6. Complete solutions

Solution 5.1. Here Rp(z)=(a−z)−1R_p(z)=(a-z)^{-1}, so Rp(z)(1+bRp(z))−1=1a−za−za+b−z=1a+b−z. R_p(z)\left(1+bR_p(z)\right)^{-1} =\frac1{a-z}\frac{a-z}{a+b-z} =\frac1{a+b-z}. Also Rp−RH=b/[(a−z)(a+b−z)]R_p-R_H=b/[(a-z)(a+b-z)], which equals both products on the right of (8) and (9). A minus sign in the factor would instead produce the denominator a−b−za-b-z.

Solution 5.2. If (A−it)uj(A-it)u_j is Cauchy, apply (5) to uj−uku_j-u_k. Both uju_j and AujAu_j are Cauchy in L2L^2, with limits u,hu,h. Closedness says u∈D(A)u\in\mathcal D(A) and Au=hAu=h. The original output limit is therefore (A−it)u(A-it)u, proving it lies in the range. Symmetry supplies (5); closedness supplies membership of the limiting input in the domain.

Solution 5.3. Every derivative of 1/(p−z)1/(p-z) is a finite sum of products of polynomial derivatives divided by a positive power of p−zp-z. The denominator's absolute value is at least ∣Im⁡z∣>0|\operatorname{Im}z|>0, so all these derivatives grow at most polynomially. Leibniz's rule then shows that multiplying any Schwartz Fourier function by the reciprocal preserves every Schwartz seminorm. At a real energy with a nonempty regular shell, the denominator vanishes on that shell. The upper and lower boundary reciprocals then contain a principal-value pole and a delta surface term, so they are distributions rather than smooth multipliers of polynomial growth. A real energy with an empty shell need not have this obstruction: for p(ξ)=∣ξ∣2p(\xi)=|\xi|^2 and λ=−1\lambda=-1, the reciprocal is (1+∣ξ∣2)−1(1+|\xi|^2)^{-1}. It is smooth, and all its derivatives have at most polynomial growth, so it preserves Schwartz space by the same Leibniz argument. Its boundary value has no singular surface contribution.

Solution 5.4. Write (VuR,vR)−(Vu,v)=(VuR−Vu,vR)+(Vu,vR−v). (Vu_R,v_R)-(Vu,v) =(Vu_R-Vu,v_R)+(Vu,v_R-v). The first term is bounded by ∥VuR−Vu∥B∥vR∥B∗\|Vu_R-Vu\|_B\|v_R\|_{B^*}, which tends to zero by strong compactness convergence and the uniform cutoff norm bound. The second is the integral of Vu v‾ (χ(x/R)−1)Vu\,\overline v\,(\chi(x/R)-1), bounded in absolute value by a constant times the integrable function ∣Vu v∣|Vu\,v|; its pointwise factor tends to zero. Dominated convergence applies. A bare distributional limit does not control the first norm error or an integral against a varying, potentially noncompact endpoint wave.

Solution 5.5. Put T0=(I+C(z0))−1T_0=(I+C(z_0))^{-1}. Factor I+C(z)=[I+(C(z)−C(z0))T0](I+C(z0)). I+C(z)= \left[I+(C(z)-C(z_0))T_0\right](I+C(z_0)). If ∥(C(z)−C(z0))T0∥≤q<1\|(C(z)-C(z_0))T_0\|\leq q<1, the geometric series in the bracket converges. Hence ∥(I+C(z))−1∥≤∥T0∥/(1−q). \|(I+C(z))^{-1}\|\leq\|T_0\|/(1-q). For any two invertible factors A′,AA',A, multiplying their difference by the inverses gives A′−1−A−1=A′−1(A−A′)A−1A'^{-1}-A^{-1}=A'^{-1}(A-A')A^{-1}. Taking A′=I+C(z′)A'=I+C(z'), A=I+C(z)A=I+C(z) yields exactly the displayed formula in Corollary 4.1.

References