Resolvents, domains and spectral density

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

Working question: What does a measured resolvent actually determine? For a scalar spectral value tt, the imaginary part of (t−λ−iε)−1(t-\lambda-i\varepsilon)^{-1} is ε/((t−λ)2+ε2)\varepsilon/((t-\lambda)^2+\varepsilon^2). Its height diverges while its integral stays bounded. This separates an operator-norm question from a scalar-measure question before either is used in scattering.

A stationary equation has two parts: a differential expression and the space on which it acts. A resolvent solves that equation away from the spectrum. Scattering begins when the spectral parameter approaches a real energy and the solution ceases to be square integrable. This lesson establishes the operator identities that survive that passage and explains what additional estimate is needed to recover a spectral density.

The prerequisites are Hilbert spaces, Fourier inversion and elementary measure theory. We use the complete bundled proof of Self-adjoint spectral calculus with the original domain. Its Cayley construction supplies the spectral measure and exact maximal multiplier domains for every self-adjoint operator, without a lower-bound or separability assumption. Section 3 supplies the compact-operator Fredholm argument directly, including its Banach-space scope; Compact Fredholm operators and their families is a companion treatment. The freely readable second edition of Gerald Teschl's Mathematical Methods in Quantum Mechanics [T], Theorem 4.3, contains the scalar-kernel proof corresponding to Section 4. Dmitri Yafaev [Y] gives scattering context.

The exact Hilbert-space prerequisites are proved in the bundled elementary Hilbert tools: orthogonal projection, bounded Riesz representation and adjoints. The measure and function-space foundations prove convergence of integrals and completeness of L2L^2. Next read products of sigma-finite measures, Euclidean products, coordinate integration and surface measure, and Fourier normalization. These supply product integration for the possibly singular spectral measures used below, measure uniqueness, Fourier inversion and Plancherel, as well as all substitutions in the examples. Section 5 constructs the bounded density it needs directly.

1. An equation with a specified domain

Our inner products are linear in the first variable. Put Dj=−i∂xjD_j=-i\partial_{x_j}, and use the unitary Fourier transform

Ff(ξ)=(2π)−n/2∫e−ix⋅ξf(x) dx. \mathcal Ff(\xi)=(2\pi)^{-n/2}\int e^{-ix\cdot\xi}f(x)\,dx.

Let p0p_0 be a real polynomial. Its Fourier realization is

H0=F−1Mp0F,D(H0)={u∈L2:p0Fu∈L2}. H_0=\mathcal F^{-1}M_{p_0}\mathcal F,\qquad D(H_0)=\{u\in L^2:p_0\mathcal Fu\in L^2\}.

This is a self-adjoint operator. The maximal multiplication domain is dense, since cutting an arbitrary L2L^2 function to ∣p0∣≤N|p_0|\leq N gives domain vectors converging to it. Multiplication by the real function is symmetric there. If vv is in the adjoint domain with value ww, testing against all L2L^2 functions supported where ∣p0∣≤N|p_0|\leq N gives 1∣p0∣≤Nw=p01∣p0∣≤Nv1_{|p_0|\leq N}w=p_0 1_{|p_0|\leq N}v. Monotone convergence then puts p0vp_0v in L2L^2, and exhaustion gives w=p0vw=p_0v. The unitary Fourier transform carries this assertion back to H0H_0.

Its spectral projections are explicitly EH0(S)=F−1M1p0−1(S)FE_{H_0}(S)=\mathcal F^{-1}M_{1_{p_0^{-1}(S)}}\mathcal F. Indicator products give orthogonal projections, and L2L^2 dominated convergence gives their strong countable additivity. The moment domain and action follow by bounded truncations of p0p_0. The uniqueness proved in the self-adjoint prerequisite identifies these with the original spectral measure. Thus all free spectral measures in the examples below have a proved realization.

If p0p_0 is elliptic of positive order mm, meaning that its homogeneous part of degree mm has no zero on the unit sphere, then

C−1⟨ξ⟩m≤1+∣p0(ξ)∣≤C⟨ξ⟩m. C^{-1}\langle\xi\rangle^m\leq 1+|p_0(\xi)| \leq C\langle\xi\rangle^m.

The upper bound is a polynomial estimate. For the lower bound, ellipticity and compactness of the sphere give ∣p0(ξ)∣≥c∣ξ∣m|p_0(\xi)|\geq c|\xi|^m outside a sufficiently large ball, and the added constant controls the ball. Consequently D(H0)=Hm(Rn)D(H_0)=H^m(\mathbb R^n), with equivalent graph and Sobolev norms. This identification concerns the whole space. A boundary would require another domain.

Theorem 1.1. Let AA be self-adjoint and let VV be bounded and self-adjoint on its Hilbert space. Then

H=A+V,D(H)=D(A), H=A+V,\qquad D(H)=D(A),

is self-adjoint, and its graph norm is equivalent to that of AA.

Proof. Symmetry is immediate on the common domain. For t>∥V∥t>\|V\|, self-adjointness of AA gives RA(±it)=(A∓it)−1R_A(\pm it)=(A\mp it)^{-1} and ∥RA(±it)∥≤t−1\|R_A(\pm it)\|\leq t^{-1}. Thus

H∓it=(I+VRA(±it))(A∓it) H\mp it=(I+VR_A(\pm it))(A\mp it)

maps D(A)D(A) bijectively onto the Hilbert space: the first factor is inverted by its norm-convergent geometric series. A densely defined symmetric operator whose ranges at both itit and −it-it are the whole space is self-adjoint. To verify this last criterion, for u∈D(H∗)u\in D(H^*) solve (H−it)v=(H∗−it)u(H-it)v=(H^*-it)u. Then w=u−vw=u-v satisfies (H∗−it)w=0(H^*-it)w=0. This kernel is the orthogonal complement of ran⁡(H+it)\operatorname{ran}(H+it), hence is zero. Therefore u=v∈D(H)u=v\in D(H). Finally

∥Hu∥≤∥Au∥+∥V∥∥u∥,∥Au∥≤∥Hu∥+∥V∥∥u∥ \|Hu\|\leq\|Au\|+\|V\|\|u\|, \qquad \|Au\|\leq\|Hu\|+\|V\|\|u\|

give both graph-norm comparisons. □\square

For example, −Δ+V-\Delta+V with real V∈L∞(Rn)V\in L^\infty(\mathbb R^n) has domain H2H^2. The theorem does not handle an unbounded potential or a perturbation of the highest derivatives. Such perturbations require estimates on their actual operator domains.

On a bounded smooth region Ω\Omega, the Dirichlet Laplacian has domain H2(Ω)∩H01(Ω)H^2(\Omega)\cap H_0^1(\Omega). The complete bundled proof of The smooth Dirichlet domain and its compact inverse, Theorems 3.1 and 4.1, supplies the full second-order boundary estimate, both domain inclusions, the compact H01H_0^1-to-L2L^2 embedding, positive self-adjoint inverse and compact resolvents. See also Hunter's PDE notes, Theorems 4.27 and 4.30; no higher-order regularity result is imported. Applying Theorem 1.1 to that realization preserves its zero-trace condition. An outgoing solution on Rn\mathbb R^n has a condition at infinity instead; it does not acquire a Dirichlet boundary condition.

2. Which resolvent identities are legitimate?

For a self-adjoint AA, write RA(z)=(A−z)−1R_A(z)=(A-z)^{-1} for Im⁡z≠0\operatorname{Im}z\ne0. This maps the Hilbert space into D(A)D(A). The bundled spectral-measure result gives, for every such AA,

RA(z)=∫(t−z)−1 dEA(t),D(A)={u:∫t2 d(EA(t)u,u)<∞}. R_A(z)=\int (t-z)^{-1}\,dE_A(t),\qquad D(A)=\left\{u:\int t^2\,d(E_A(t)u,u)<\infty\right\}.

In particular ∥RA(z)∥≤∣Im⁡z∣−1\|R_A(z)\|\leq|\operatorname{Im}z|^{-1}. For a general self-adjoint operator the same bound follows from ∥(A−z)u∥≥∣Im⁡z∣∥u∥\|(A-z)u\|\geq|\operatorname{Im}z|\|u\|, applied also to its adjoint to prove surjectivity.

Proposition 2.1. For nonreal z,wz,w,

RA(z)−RA(w)=(z−w)RA(z)RA(w),RA(z)∗=RA(z‾). R_A(z)-R_A(w)=(z-w)R_A(z)R_A(w),\qquad R_A(z)^*=R_A(\overline z).

If H=A+VH=A+V is as in Theorem 1.1, then

RH(z)−RA(z)=−RH(z)VRA(z)=−RA(z)VRH(z), R_H(z)-R_A(z)=-R_H(z)VR_A(z)=-R_A(z)VR_H(z),

and

RH(z)=RA(z)(I+VRA(z))−1. R_H(z)=R_A(z)(I+VR_A(z))^{-1}.

Proof. Apply both sides of the first identity to an arbitrary vector, and insert (A−w)−(A−z)=(z−w)I(A-w)-(A-z)=(z-w)I between the two inverses. Their ranges lie in D(A)D(A), so every application of AA is legitimate. The adjoint identity follows from (A−z)∗=A−z‾(A-z)^*=A-\overline z. For the perturbation formula insert (A−z)−(H−z)=−V(A-z)-(H-z)=-V between the inverses, in either order. Finally I+VRA(z)I+VR_A(z) is invertible because

(I+VRA(z))−1=I−VRH(z). (I+VR_A(z))^{-1}=I-VR_H(z).

Multiplying these two bounded factors on either side, using the preceding identities, gives the identity operator. This also proves the final formula without a smallness hypothesis on VV. □\square

The placement of the factors matters. VRA(z)VR_A(z) and RA(z)VR_A(z)V have different domains and target spaces once weighted estimates replace Hilbert-space bounds. A real-axis version of the formula must first identify the Banach space on which I+VRA(λ±i0)I+VR_A(\lambda\pm i0) acts.

3. Decaying potentials and compact errors

We now make one exact use of Fredholm theory. On a Banach space, a compact perturbation of the identity is Fredholm with index zero. Its kernel and cokernel are finite-dimensional; it is invertible exactly when its kernel is zero. Here is the complete compactness argument. The Banach-space foundation used is Hahn–Banach, Corollary 2.3(1), for extending the finitely many coordinate functionals of a finite-dimensional subspace; its full real and complex proofs are given in Theorems 2.1–2.2 there. The companion reading's Elementary finite-dimensional tools proves bounded coordinates, closed complements and completeness of a quotient by a closed subspace; read that paragraph before the argument below.

Proof of the compact-perturbation input. Let XX be a Banach space, K:X→XK:X\to X compact, and T=I+KT=I+K. Its kernel NN is finite-dimensional: on NN, K=−IK=-I, so its unit ball is relatively compact. An infinite-dimensional normed space has a sequence of unit vectors separated by more than 1/21/2, contradicting that compactness. This follows from the elementary separation lemma for a proper closed subspace Y⊂ZY\subset Z: take z∉Yz\notin Y, set d=dist⁡(z,Y)>0d=\operatorname{dist}(z,Y)>0, choose y∈Yy\in Y with ∥z−y∥<2d\|z-y\|<2d, and normalize z−yz-y to obtain a unit vector at distance more than 1/21/2 from YY. Apply it successively to the spans of the preceding vectors.

There is a constant CC such that

dist⁡(x,N)≤C∥Tx∥(x∈X). \operatorname{dist}(x,N)\leq C\|Tx\|\qquad(x\in X).

Otherwise choose representatives xjx_j with distance one from NN, norm at most two, and Txj→0Tx_j\to0. Subtracting an element of NN gives the norm bound without changing either the distance or image. A subsequence of KxjKx_j converges, so xj=Txj−Kxjx_j=Tx_j-Kx_j converges to a vector in NN, contradicting distance one. This also proves that R=T(X)R=T(X) is closed: if Txj→yTx_j\to y, the displayed bound lets us replace xjx_j modulo NN by bounded representatives. Compactness of their KK-images then gives a convergent subsequence of the representatives and an inverse image of yy.

The quotient X/RX/R is finite-dimensional. For its quotient map qq, the identity q(x)=−q(Kx)q(x)=-q(Kx) holds because Tx∈RTx\in R. Every quotient vector of norm at most one has a representative of norm at most two. Its unit ball is therefore relatively compact, being contained in the image under −qK-qK of that bounded ball. The separation lemma again excludes infinite dimension.

Compactness is preserved by finite sums and by composition on either side with bounded operators. For sums, successively extract subsequences for the finitely many compact images of a bounded sequence; for composition, use continuity on the image side and boundedness on the input side. The metric compactness criterion was proved in the companion reading. Hence every positive power Tm=(I+K)mT^m=(I+K)^m is identity plus compact, by the finite product expansion. We next prove that injectivity and surjectivity are equivalent. If TT is injective but not surjective, the closed ranges Rm=Tm(X)R_m=T^m(X) decrease strictly. They are closed because Tm=I+KmT^m=I+K_m with compact KmK_m, and the preceding closed-range proof applies. Equality of two successive ranges would, by injectivity, imply equality of the preceding pair, eventually contradicting X≠T(X)X\ne T(X). Choose unit xm∈Rm−1x_m\in R_{m-1} at distance more than 1/21/2 from RmR_m. For k>mk>m, the vector Kxm−KxkKx_m-Kx_k equals −xm-x_m plus a vector in RmR_m; hence its norm exceeds 1/21/2. This contradicts compactness. If TT is surjective but has a nonzero kernel, the closed spaces Nm=ker⁡TmN_m=\ker T^m increase strictly: successive preimages of a nonzero kernel vector exhibit the strict inclusions. Choose unit xm∈Nmx_m\in N_m at distance more than 1/21/2 from Nm−1N_{m-1}. For k<mk<m, Kxm−KxkKx_m-Kx_k equals −xm-x_m plus a vector in Nm−1N_{m-1}, giving the same contradiction.

Finally put n=dim⁡Nn=\dim N, c=dim⁡(X/R)c=\dim(X/R), and s=min⁡(n,c)s=\min(n,c). Choose complements X=N⊕X0=R⊕Y0X=N\oplus X_0=R\oplus Y_0, with dim⁡Y0=c\dim Y_0=c. The first complement is the intersection of the kernels of bounded coordinate extensions supplied by Hahn–Banach; the second is obtained by lifting a quotient basis. Define a bounded finite-rank FF to map ss basis vectors of NN to ss independent basis vectors of Y0Y_0, and to vanish on the other kernel basis vectors and on X0X_0. Since T(X0)=RT(X_0)=R, the operator T+FT+F has kernel dimension n−sn-s and cokernel dimension c−sc-s. It is again identity plus compact. If n≠cn\ne c, it would be injective without being surjective, or surjective without being injective, contradicting the preceding paragraph. Thus n=cn=c, which is exactly index zero. When n=0n=0, surjectivity follows, and the displayed distance estimate becomes ∥x∥≤C∥Tx∥\|x\|\leq C\|Tx\|, proving boundedness of the inverse. □\square

We apply this result to a bounded operator on L2L^2; the domain of the differential realization is retained separately below.

Theorem 3.1. Suppose p0p_0 is real and elliptic of order m>0m>0, and VV is a bounded multiplication operator satisfying

ess sup∣x∣>R∣V(x)∣⟶0. \mathop{\mathrm{ess\,sup}}_{|x|>R}|V(x)|\longrightarrow0.

For every nonreal zz, VRH0(z)VR_{H_0}(z) is compact on L2L^2. If VV is real, then RH0+V(z)−RH0(z)R_{H_0+V}(z)-R_{H_0}(z) is compact.

Proof. The multiplier rz(ξ)=(p0(ξ)−z)−1r_z(\xi)=(p_0(\xi)-z)^{-1} is bounded and tends to zero as ∣ξ∣→∞|\xi|\to\infty. Choose a smooth cutoff χN\chi_N equal to one on ∣ξ∣≤N|\xi|\leq N and zero outside ∣ξ∣≤2N|\xi|\leq2N. Then rzχN∈L2r_z\chi_N\in L^2, and its inverse unitary Fourier transform kNk_N is in L2L^2. For VR=V1{∣x∣≤R}V_R=V1_{\{|x|\leq R\}}, the kernel of VRF−1MrzχNFV_R\mathcal F^{-1}M_{r_z\chi_N}\mathcal F is

(2π)−n/2VR(x)kN(x−y). (2\pi)^{-n/2}V_R(x)k_N(x-y).

Its square integral is (2π)−n∥VR∥22∥kN∥22(2\pi)^{-n}\|V_R\|_2^2\|k_N\|_2^2, which is finite. Such an integral operator is compact: finite sums of product functions are dense in the product L2L^2 space, their operators have finite rank, and the kernel L2L^2 norm bounds the operator norm by Cauchy–Schwarz. Here rzχNr_z\chi_N is smooth with compact support, so kNk_N is Schwartz by Fourier inversion. Absolute Fubini first identifies the displayed kernel on compact smooth inputs; boundedness and density extend the identity to all L2L^2. For the density assertion, first truncate an arbitrary product-space L2L^2 function to a bounded box, then use finite simple approximation and the proved Euclidean rectangle approximation. A product-space rectangle is a product of two rectangles. This is the same finite-rank kernel argument proved explicitly in the Dirichlet compactness reading, now with either factor allowed to exhaust the full Euclidean space.

The error from removing the Fourier cutoff is bounded by ∥V∥∞sup⁡∣ξ∣>N∣rz(ξ)∣\|V\|_\infty\sup_{|\xi|>N}|r_z(\xi)|. The error from replacing VV by VRV_R is bounded by ∥V−VR∥∞∥rz∥∞\|V-V_R\|_\infty\|r_z\|_\infty. Let NN and then RR tend to infinity. A norm limit of compact operators is compact, since a finite approximation to the image of a unit ball can be enlarged by the small norm error. The resolvent difference is a bounded operator times VRH0(z)VR_{H_0}(z), by Proposition 2.1, so it is compact too. □\square

If λ\lambda lies outside the spectrum of H0H_0, the same proof gives compactness of VRH0(λ)VR_{H_0}(\lambda). Indeed, p0(ξ0)=λp_0(\xi_0)=\lambda would make normalized indicators of shrinking frequency balls around ξ0\xi_0 into unit vectors whose (H0−λ)(H_0-\lambda)-images tend to zero, contradicting a bounded inverse. Thus p0−λp_0-\lambda never vanishes. Its absolute value has a positive minimum on each compact ball and grows at infinity by ellipticity, so its reciprocal is bounded and tends to zero. The factorization

H0+V−λ=(I+VRH0(λ))(H0−λ) H_0+V-\lambda=(I+VR_{H_0}(\lambda))(H_0-\lambda)

then shows that H0+V−λH_0+V-\lambda, regarded as a bounded map from D(H0)D(H_0) with its graph norm to L2L^2, is Fredholm of index zero. The map H0−λH_0-\lambda is a bounded isomorphism from that complete graph space to L2L^2: its inverse is bounded in graph norm because H0RH0(λ)=I+λRH0(λ)H_0R_{H_0}(\lambda)=I+\lambda R_{H_0}(\lambda). Composition with an isomorphism preserves the kernel dimension, range closedness and cokernel dimension. This conclusion follows from the compactness argument just proved; no index formula for a symbol is being asserted. When its kernel is nonzero, the kernel is a finite-dimensional eigenspace. Inside the free spectrum this factorization has no bounded Hilbert-space inverse and needs different spaces.

4. A spectral measure seen through a resolvent

Let AA be lower bounded and self-adjoint, and put μf(S)=(EA(S)f,f)\mu_f(S)=(E_A(S)f,f). The spectral theorem from Section 2 gives

(RA(λ+iε)f,f)=∫1t−λ−iε dμf(t), (R_A(\lambda+i\varepsilon)f,f) =\int\frac{1}{t-\lambda-i\varepsilon}\,d\mu_f(t),

and therefore

Im⁡(RA(λ+iε)f,f)=∫ε(t−λ)2+ε2 dμf(t). \operatorname{Im}(R_A(\lambda+i\varepsilon)f,f) =\int\frac{\varepsilon}{(t-\lambda)^2+\varepsilon^2}\,d\mu_f(t).

The right side is nonnegative. The sign reverses in the lower half-plane.

Theorem 4.1 (Stone's formula, scalar form). For finite a<ba<b,

lim⁡ε↓01π∫abIm⁡(RA(λ+iε)f,f) dλ=μf((a,b))+12μf({a})+12μf({b}). \begin{gathered} \lim_{\varepsilon\downarrow0}\frac1\pi \int_a^b\operatorname{Im}(R_A(\lambda+i\varepsilon)f,f)\,d\lambda\\ =\mu_f((a,b))+\tfrac12\mu_f(\{a\})+\tfrac12\mu_f(\{b\}). \end{gathered}

Proof. Apply the proved general Tonelli theorem to the finite Borel measure μf\mu_f, of mass ∥f∥2\|f\|^2, and Lebesgue measure on [a,b][a,b]. For each fixed ε>0\varepsilon>0, the kernel ε/((t−λ)2+ε2)\varepsilon/((t-\lambda)^2+\varepsilon^2) is nonnegative and jointly Borel measurable: it is continuous, and open subsets of the real plane are countable unions of rational open rectangles. Both measures are sigma-finite, so the theorem permits exchanging these two integrals even when μf\mu_f is singular. Integrating first in λ\lambda, the inner integral divided by π\pi is

1π[arctan⁡b−tε−arctan⁡a−tε]. \frac1\pi\left[ \arctan\frac{b-t}{\varepsilon}-\arctan\frac{a-t}{\varepsilon} \right].

It is between zero and one. Its pointwise limit is one for a<t<ba<t<b, one half at either endpoint, and zero outside [a,b][a,b]. The measure μf\mu_f has finite mass ∥f∥2\|f\|^2, so dominated convergence proves the formula. □\square

An operator-valued version follows by polarization. For example, with our convention the form B(f,g)=(EA(S)f,g)B(f,g)=(E_A(S)f,g) is recovered from Q(h)=B(h,h)Q(h)=B(h,h) as

B(f,g)=14(Q(f+g)−Q(f−g)+iQ(f+ig)−iQ(f−ig)). \begin{aligned} B(f,g)&=\tfrac14\bigl(Q(f+g)-Q(f-g)\\ &\qquad+iQ(f+ig)-iQ(f-ig)\bigr). \end{aligned}

The operator formula in fact converges strongly. Let hε(t)h_\varepsilon(t) be the normalized arctangent difference in the proof. It is bounded by one and tends pointwise to 1(a,b)+121{a,b}1_{(a,b)}+\tfrac12 1_{\{a,b\}}. The bounded Borel convergence theorem therefore gives strong convergence of hε(A)h_\varepsilon(A) to the corresponding projection sum. This operator is the integral, divided by π\pi, of (RA(λ+iε)−RA(λ−iε))/(2i)(R_A(\lambda+i\varepsilon)-R_A(\lambda-i\varepsilon))/(2i): the resolvent identity makes the integrand norm continuous on the finite interval, the earlier Banach-valued integral exists, and scalar Tonelli followed by polarization identifies all its pairings with hε(A)h_\varepsilon(A). This proves the operator assertion on the whole Hilbert space.

Theorem 4.2. Let II be an open real interval and let ff be a Hilbert-space vector. Suppose

Im⁡(RA(λ+iε)f,f)⟶gf(λ) \operatorname{Im}(R_A(\lambda+i\varepsilon)f,f) \longrightarrow g_f(\lambda)

locally uniformly in λ∈I\lambda\in I, where gfg_f is continuous. Then μf\mu_f on II is absolutely continuous and

dμf(λ)=π−1gf(λ) dλ. d\mu_f(\lambda)=\pi^{-1}g_f(\lambda)\,d\lambda.

Proof. The imaginary parts are locally bounded. At an atom t∈It\in I, the contribution μf({t})/ε\mu_f(\{t\})/\varepsilon to the value at λ=t\lambda=t forces μf({t})=0\mu_f(\{t\})=0. Stone's formula and uniform convergence on [a,b]⊂I[a,b]\subset I now give μf((a,b))=π−1∫abgf(λ) dλ\mu_f((a,b))=\pi^{-1}\int_a^bg_f(\lambda)\,d\lambda. Open intervals determine finite Borel measures locally: atomlessness gives equality on half-open intervals as well, and the proved finite-measure uniqueness lemma gives equality on all Borel sets in any fixed compact subinterval. Exhaustion proves the assertion on II. □\square

If this hypothesis holds for every ff in a dense linear subspace D\mathcal D, then every spectral measure is absolutely continuous on II. Indeed, for a Lebesgue-null Borel set N⊂IN\subset I, ∥EA(N)f∥2=μf(N)=0\|E_A(N)f\|^2=\mu_f(N)=0 for f∈Df\in\mathcal D. Continuity of the orthogonal projection implies EA(N)=0E_A(N)=0. This proves absence of singular spectrum there. Theorem 5.1 below shows that the local bound alone suffices for this last conclusion; convergence supplies the stronger continuous-density formula in Theorem 4.2. Existence of a limit at almost every energy without local control is a different hypothesis and does not give either argument.

Converse for a continuous density. If a finite positive Borel measure μ\mu has continuous density hh on an open interval II, then its Poisson integrals converge locally uniformly to πh\pi h there. Here is the proof needed for the free example. Given compact K⊂IK\subset I, choose a smooth cutoff 0≤χ≤10\leq\chi\leq1, supported in II and equal to one on a neighbourhood of KK. Continuity and positivity of the measure imply h≥0h\geq0: a negative value would give a negative integral over a small interval. Extend g=χhg=\chi h by zero; it is bounded and uniformly continuous on the line. The positive finite measure ν=μ−g dt\nu=\mu-g\,dt is supported at some positive distance dd from KK. Its Poisson integral on KK is at most εμ(R)/d2\varepsilon\mu(\mathbb R)/d^2.

The scalar kernel Pε(r)=ε/(r2+ε2)P_\varepsilon(r)=\varepsilon/(r^2+\varepsilon^2) has integral π\pi, by the proved arctangent primitive. Splitting its convolution with g(λ−r)−g(λ)g(\lambda-r)-g(\lambda) into ∣r∣<δ|r|<\delta and its complement gives, uniformly in λ\lambda, ∣∫Pε(r)g(λ−r) dr−πg(λ)∣≤πωg(δ)+4ε∥g∥∞δ. \begin{gathered} \left|\int P_\varepsilon(r)g(\lambda-r)\,dr-\pi g(\lambda)\right|\\ \leq\pi\omega_g(\delta)+\frac{4\varepsilon\|g\|_\infty}{\delta}. \end{gathered} Indeed the first part is bounded by the modulus of continuity times π\pi, while ∫∣r∣≥δPε(r)dr≤2ε/δ\int_{|r|\geq\delta}P_\varepsilon(r)dr\leq2\varepsilon/\delta. First choose δ\delta small, then ε\varepsilon small. On KK, g=hg=h, and the vanishing contribution of ν\nu proves the claimed uniform limit. No pointwise boundary theorem is assumed.

5. The topology at a real energy

For s∈Rs\in\mathbb R, put

Ls2={f:⟨x⟩sf∈L2},∥f∥Ls2=∥⟨x⟩sf∥2. L^2_s=\{f:\langle x\rangle^sf\in L^2\},\qquad \|f\|_{L^2_s}=\|\langle x\rangle^sf\|_2.

The integral pairing identifies L−s2L^2_{-s} with the continuous dual of Ls2L^2_s, with conjugation according to our inner-product convention. Indeed f↦⟨x⟩sff\mapsto\langle x\rangle^sf is an isometry onto L2L^2, with inverse multiplication by ⟨x⟩−s\langle x\rangle^{-s}. The proved Hilbert representation applied after this isometry gives exactly the pairing with ⟨x⟩sg\langle x\rangle^sg, g∈L2g\in L^2, which lies in L−s2L^2_{-s}. Cauchy–Schwarz gives its norm and the reverse equality follows by taking the corresponding unit vector. When s>0s>0, Ls2⊂L2⊂L−s2L^2_s\subset L^2\subset L^2_{-s}. A limiting absorption estimate has the form

sup⁡λ∈K, 0<ε<1∥RA(λ±iε)∥Ls2→L−s2<∞, \sup_{\lambda\in K,\ 0<\varepsilon<1} \|R_A(\lambda\pm i\varepsilon)\|_{L^2_s\to L^2_{-s}}<\infty,

for specified compact energy sets KK. This bound already has a spectral consequence, before a boundary family is constructed.

Theorem 5.1 (a bound excludes singular spectral mass). Let AA have the spectral-measure representation of Section 4, and let II be an open real interval. Suppose that, for a vector ff and every compact interval J⊂IJ\subset I, there are finite Cf,JC_{f,J} and positive εJ\varepsilon_J such that

0≤Fε,f(λ):=Im⁡(RA(λ+iε)f,f)≤Cf,J,λ∈J,0<ε<εJ. 0\leq F_{\varepsilon,f}(\lambda) :=\operatorname{Im}(R_A(\lambda+i\varepsilon)f,f) \leq C_{f,J},\qquad \lambda\in J,\quad 0<\varepsilon<\varepsilon_J.

Then μf\mu_f is absolutely continuous on II. On the interior of JJ its density satisfies

0≤dμfdλ≤Cf,Jπalmost everywhere. 0\leq\frac{d\mu_f}{d\lambda}\leq\frac{C_{f,J}}{\pi} \quad\hbox{almost everywhere}.

If these bounds hold for every vector in a dense subspace, AA has no singular spectrum in II. In particular, the displayed Ls2L^2_s-to-L−s2L^2_{-s} estimate, for every compact interval in II and s>0s>0, implies this conclusion. Only its upper-half-plane bound is needed.

Proof. Integrate the nonnegative inequality over [a,b]⊂J∘[a,b]\subset J^\circ, then use Theorem 4.1, retaining its endpoint terms:

μf((a,b))+12μf({a})+12μf({b})≤Cf,Jπ(b−a). \mu_f((a,b))+\tfrac12\mu_f(\{a\})+\tfrac12\mu_f(\{b\}) \leq\frac{C_{f,J}}{\pi}(b-a).

At any t∈J∘t\in J^\circ, positivity of the original kernel gives μf({t})/ε≤Fε,f(t)≤Cf,J\mu_f(\{t\})/\varepsilon\leq F_{\varepsilon,f}(t)\leq C_{f,J}. Letting ε↓0\varepsilon\downarrow0 kills every such atom. Thus each open interval in J∘J^\circ has measure at most Cf,J/πC_{f,J}/\pi times its length; intervals meeting an endpoint of J∘J^\circ follow by increasing exhaustion. Every open subset of J∘J^\circ is a countable disjoint union of intervals. To see this, join two points when the interval between them lies in the open set. The equivalence classes are disjoint open intervals; each contains a rational, so there are at most countably many. Countable additivity therefore gives the same bound for that open set.

For a Borel set SS in a smaller compact interval K⊂J∘K\subset J^\circ, cover SS by open sets inside J∘J^\circ whose lengths decrease to its Lebesgue measure. Such covers exist by Lebesgue outer regularity, since KK has positive distance from the complement of J∘J^\circ. Monotonicity then gives

μf(S)≤Cf,Jπ∣S∣. \mu_f(S)\leq\frac{C_{f,J}}{\pi}|S|.

Here is the density construction, including its existence. On a compact interval K⊂J∘K\subset J^\circ, write ν=μf∣K\nu=\mu_f|_K and C=Cf,J/πC=C_{f,J}/\pi. For a simple Borel function φ\varphi on KK, domination gives

∣∫Kφ dν∣≤C∫K∣φ∣ dλ≤C∣K∣1/2∥φ∥L2(K). \left|\int_K\varphi\,d\nu\right| \le C\int_K|\varphi|\,d\lambda \le C|K|^{1/2}\|\varphi\|_{L^2(K)}.

Null sets have zero ν\nu-measure, so this is well defined on L2L^2 classes. Simple-function density extends it to a bounded linear functional on L2(K)L^2(K). The Hilbert representation proved in the prerequisite gives hK∈L2(K)h_K\in L^2(K) with ∫φ dν=∫φhK dλ\int\varphi\,d\nu=\int\varphi h_K\,d\lambda; with linear-first inner products, hKh_K is the conjugate of the representing vector. Testing indicators of all measurable sets shows that hKh_K is real and 0≤hK≤C0\le h_K\le C almost everywhere. For detail, a set on which its imaginary part is at least 1/m1/m, or at most −1/m-1/m, would give a nonreal integral. A set where its real part is below −1/m-1/m, or above C+1/mC+1/m, contradicts respectively positivity or domination. Taking their countable unions proves the assertions.

Thus ν(S)=∫ShK dλ\nu(S)=\int_S h_K\,d\lambda for every Borel S⊂KS\subset K. Densities obtained on overlapping intervals agree almost everywhere, because their difference integrates to zero on every measurable subset of the overlap; the same level-set argument applies. A countable compact exhaustion of J∘J^\circ, and then of II, patches these densities into the asserted local density with the stated bound. This proves the needed bounded-density result without assuming a general Radon–Nikodym theorem.

For the dense-subspace assertion, a null Borel N⊂IN\subset I satisfies ∥EA(N)f∥2=μf(N)=0\|E_A(N)f\|^2=\mu_f(N)=0 on that subspace. The projection has norm at most one, so approximation extends its zero action to the entire Hilbert space. Constants may depend on the test vector.

Finally, if the weighted resolvent norm on JJ is bounded by CJC_J, weighted duality gives

0≤Fε,f(λ)≤∣(RA(λ+iε)f,f)∣≤CJ∥f∥Ls2 2. 0\leq F_{\varepsilon,f}(\lambda) \leq |(R_A(\lambda+i\varepsilon)f,f)| \leq C_J\|f\|_{L^2_s}^{\,2}.

Compactly supported L2L^2 functions belong to Ls2L^2_s and are dense in L2L^2. Apply the preceding argument, with Cf,J=CJ∥f∥Ls22C_{f,J}=C_J\|f\|_{L^2_s}^2. □\square

The proof uses only the positive kernel and the spectral projections. It therefore also applies to a self-adjoint operator unbounded in both directions when its general spectral measure has been supplied. Related smoothness criteria are discussed in Yafaev [Y], Section 2; the scalar measure argument above supplies the consequence used here directly.

To define boundary values one still needs convergence in a specified topology, for example weak convergence against every vector in Ls2L^2_s. A locally uniform operator-norm limit is norm continuous and satisfies Theorem 4.2 on that dense subspace, giving a continuous density. For completeness, at fixed ε>0\varepsilon>0 the resolvent identity bounds the difference at two real energies by ∣λ−μ∣ε−2|\lambda-\mu|\varepsilon^{-2} in L2L^2 operator norm. The continuous inclusions Ls2⊂L2⊂L−s2L^2_s\subset L^2\subset L^2_{-s}, each of norm at most one, give the same continuity in the weighted operator norm. A uniform limit on each compact energy interval is continuous by the triangle inequality. Theorem 5.1 gives a locally bounded density and does not assert continuity or locally uniform convergence of the imaginary parts.

The endpoint model makes this distinction quantitative. Read the earlier prerequisite Flat transport traces and the norm boundary limit for the definitions of B,B∗,B0∗B,B^*,B^*_0, the flat trace and all the following proofs. Its Theorem 5.2 proves for DtD_t that the real-energy boundary solution has exact B∗B^*-distance π∥Tλf∥2\sqrt\pi\|T_\lambda f\|_2 from B0∗B^*_0. Every nonreal resolvent lies in B0∗B^*_0, so a nonzero trace prevents norm convergence although the weak-star boundary exists. Zero trace is also sufficient for norm convergence along any approach in the corresponding half-plane, by the complete zero-integral approximation proof there. Proposition 3.3 proves that the trace operators are nowhere norm continuous, while Theorem 3.2 proves continuity on each fixed forcing term. These arguments require only the earlier measure, Fourier and integral foundations.

Example 5.1. For A=−ΔA=-\Delta and λ>0\lambda>0, the L2L^2 resolvent norm is exactly ε−1\varepsilon^{-1}. The Fourier multiplier has that supremum because every neighbourhood of the sphere ∣ξ∣2=λ|\xi|^2=\lambda has positive measure. It approaches its maximum there. Thus a weighted limit may exist while the L2L^2-operator norm diverges. This is compatible with every preceding theorem.

Example 5.2. For ff with Ff∈Cc∞(Rn)\mathcal Ff\in C_c^\infty(\mathbb R^n), the spectral measure of the free Laplacian on positive energies has density

dμfdλ=12λ∫∣ξ∣=λ∣Ff(ξ)∣2 dS(ξ). \frac{d\mu_f}{d\lambda} =\frac1{2\sqrt\lambda} \int_{|\xi|=\sqrt\lambda}|\mathcal Ff(\xi)|^2\,dS(\xi).

In polar coordinates, dξ=rn−1dr dωd\xi=r^{n-1}dr\,d\omega and dλ=2r drd\lambda=2r\,dr, which proves the formula. On compact positive-energy intervals this density is smooth. The continuous-density Poisson argument proves locally uniform convergence of the upper resolvent's imaginary part to π\pi times this density, in agreement with Stone's formula. The energy zero requires a separate analysis; no regular-energy claim here includes it.

Example 5.3 (a bound with a discontinuous limit). On L2([−2,2],dt)L^2([-2,2],dt), let A=MtA=M_t and f(t)=1f(t)=1. The operator is bounded and self-adjoint on the full space, ∥f∥2=4\|f\|^2=4, and μf\mu_f is Lebesgue measure restricted to [−2,2][-2,2]. Direct integration gives

Fε,f(λ)=arctan⁡2−λε−arctan⁡−2−λε. F_{\varepsilon,f}(\lambda) =\arctan\frac{2-\lambda}{\varepsilon} -\arctan\frac{-2-\lambda}{\varepsilon}.

For every λ\lambda and ε>0\varepsilon>0, this lies between zero and π\pi, the total mass of the positive kernel. Its limit is π\pi for ∣λ∣<2|\lambda|<2, π/2\pi/2 at λ=±2\lambda=\pm2, and zero for ∣λ∣>2|\lambda|>2. Each approximating function is continuous; their discontinuous limit prevents locally uniform convergence on any interval containing either endpoint. Nevertheless Theorem 5.1 gives the sharp density bound dμf/dλ≤1d\mu_f/d\lambda\leq1. This example concerns the scalar bound for this vector, not an operator-norm bound on all of L2([−2,2])L^2([-2,2]).

Bounded imaginary resolvent pairings approach a discontinuous spectral density at the two endpoints.

The curves sample the exact arctangent formula in Example 5.3 at ε=1/2,1/5,1/20\varepsilon=1/2,1/5,1/20. The endpoint dots are exactly (−2,π/2)(-2,\pi/2) and (2,π/2)(2,\pi/2). The horizontal bound π\pi and the density conclusion follow from the proof, rather than from the plotted samples.

Use the conclusion

Compare Example 5.1 with Example 5.3: the first has a diverging Hilbert-space resolvent norm; the second has a bounded scalar observation and a discontinuous boundary density. State which estimate each theorem needs.

6. Exercises

Exercise 6.1 (foundation). Let Aϕ=αϕA\phi=\alpha\phi, with ∥ϕ∥=1\|\phi\|=1. Compute (RA(λ+iε)ϕ,ϕ)(R_A(\lambda+i\varepsilon)\phi,\phi), and evaluate Stone's formula on intervals with α\alpha in the interior and at an endpoint.

Exercise 6.2 (foundation). Take V(x)=2/(1+∣x∣2)V(x)=2/(1+|x|^2). State the domain of −Δ+V-\Delta+V on Rn\mathbb R^n, give an explicit shift making it positive, and decide whether Theorem 3.1 applies. Explain why it does not establish a real-axis limiting absorption estimate.

Exercise 6.3 (intermediate). A formula is proposed as RH(z)=RA(z)(I−VRA(z))−1R_H(z)=R_A(z)(I-VR_A(z))^{-1} for H=A+VH=A+V. Test it when the Hilbert space is C\mathbb C, A=3A=3 and V=2V=2. Derive the correct factorization without commuting any factors.

Exercise 6.4 (intermediate). Suppose RA(λ+iε)R_A(\lambda+i\varepsilon) is locally uniformly bounded from Ls2L^2_s to L−s2L^2_{-s}, with s>0s>0. Prove absence of singular spectrum and give a local bound for dμf/dλd\mu_f/d\lambda, f∈Ls2f\in L^2_s. If there is also a locally uniform operator-norm limit, show that this density is continuous. Identify where the dense test space is used, and explain why Example 5.3 separates the two scalar hypotheses.

Exercise 6.5 (advanced). Let p0(ξ)=ξ12+4ξ22+ξ32p_0(\xi)=\xi_1^2+4\xi_2^2+\xi_3^2. For Ff∈Cc∞(R3)\mathcal Ff\in C_c^\infty(\mathbb R^3), compute the free spectral density by the substitution η=(ξ1,2ξ2,ξ3)\eta=(\xi_1,2\xi_2,\xi_3). Compare it with the surface formula ∫p0=λ∣Ff∣2∣∇p0∣−1 dS\int_{p_0=\lambda}|\mathcal Ff|^2|\nabla p_0|^{-1}\,dS.

7. Complete solutions

Solution 6.1. The scalar value is (α−λ−iε)−1(\alpha-\lambda-i\varepsilon)^{-1}, with imaginary part ε/((α−λ)2+ε2)\varepsilon/((\alpha-\lambda)^2+\varepsilon^2). Its integral divided by π\pi tends to one if a<α<ba<\alpha<b, one half if α=a\alpha=a or α=b\alpha=b, and zero if α∉[a,b]\alpha\notin[a,b]. At λ=α\lambda=\alpha the imaginary part is ε−1\varepsilon^{-1}. Therefore the locally bounded hypothesis in Theorem 4.2 excludes this atom exactly as required.

Solution 6.2. The domain is H2(Rn)H^2(\mathbb R^n), since VV is real and bounded. The operator is already nonnegative; adding II makes it at least II. Moreover ess sup∣x∣>R∣V(x)∣=2/(1+R2)→0\mathop{\mathrm{ess\,sup}}_{|x|>R}|V(x)|=2/(1+R^2)\to0, so the resolvent difference at any nonreal zz is compact. The compactness proof uses the multiplier (∣ξ∣2−z)−1(|\xi|^2-z)^{-1} in L∞L^\infty; its supremum diverges as zz approaches a positive real number. Thus the proof provides no uniform weighted bound at those energies.

Solution 6.3. The actual inverse is (5−z)−1(5-z)^{-1}. The proposed expression gives (1/(3−z))/(1−2/(3−z))=(1−z)−1(1/(3-z))/(1-2/(3-z))=(1-z)^{-1}. The correct identity is

H−z=(I+V(A−z)−1)(A−z) H-z=(I+V(A-z)^{-1})(A-z)

on D(A)D(A). Inverting reverses the factors and gives (A−z)−1(I+V(A−z)−1)−1(A-z)^{-1}(I+V(A-z)^{-1})^{-1}. No interchange of VV with AA is needed.

Solution 6.4. Write CJC_J for the weighted norm bound on a compact energy interval JJ. For f∈Ls2f\in L^2_s, weighted duality bounds the positive scalar imaginary part by CJ∥f∥Ls22C_J\|f\|_{L^2_s}^2. Stone's formula gives interval domination with coefficient CJ∥f∥Ls22/πC_J\|f\|_{L^2_s}^2/\pi, and the value at an atom bounds its mass divided by ε\varepsilon. The atom is therefore zero. Open interval decomposition and outer regularity, as in Theorem 5.1, give the same domination on every Borel set locally. Hence the density is bounded by that coefficient almost everywhere. Spectral projections annihilate the dense Ls2L^2_s test space on null sets and extend by continuity to all L2L^2; this proves absence of singular spectrum without a boundary limit.

With the additional operator-norm limit, the integral pairing bounds the scalar difference by ∥f∥Ls22\|f\|_{L^2_s}^2 times the operator-norm difference. The scalar imaginary part therefore converges locally uniformly to a continuous function, and Theorem 4.2 identifies the density with that function divided by π\pi. Compactly supported L2L^2 functions supply the dense test space; continuity of the projections extends their zero action to arbitrary vectors. The pairing itself does not require an L2L^2 boundary solution. In Example 5.3 the scalar bound is π\pi, but the limit jumps at ±2\pm2, excluding locally uniform convergence there. It demonstrates the difference between the scalar hypotheses without asserting a weighted operator estimate for that example.

Solution 6.5. The Jacobian is dξ=12dηd\xi=\tfrac12d\eta. If G(η)=Ff(η1,η2/2,η3)G(\eta)=\mathcal Ff(\eta_1,\eta_2/2,\eta_3), then

μf((0,λ])=12∫∣η∣2≤λ∣G(η)∣2 dη,dμfdλ=14λ∫∣η∣=λ∣G(η)∣2 dSη. \mu_f((0,\lambda])=\tfrac12\int_{|\eta|^2\leq\lambda}|G(\eta)|^2\,d\eta, \qquad \frac{d\mu_f}{d\lambda} =\frac1{4\sqrt\lambda}\int_{|\eta|=\sqrt\lambda}|G(\eta)|^2\,dS_\eta.

For the coarea comparison, parametrize the ellipsoid by ξ=Tη\xi=T\eta, T=diag⁡(1,1/2,1)T=\operatorname{diag}(1,1/2,1). Surface area transforms by dSξ=∣det⁡T∣ ∣T−Tω∣ dSηdS_\xi=|\det T|\,|T^{-T}\omega|\,dS_\eta, where ω=η/∣η∣\omega=\eta/|\eta|. Also ∇ξp0=2λ T−Tω\nabla_\xi p_0=2\sqrt\lambda\,T^{-T}\omega. Their quotient is dSη/(4λ)dS_\eta/(4\sqrt\lambda), agreeing with the calculation. Both the Jacobian and the gradient are necessary.

To verify the area transformation, choose orthonormal tangent vectors e1,e2e_1,e_2 so that (e1,e2,ω)(e_1,e_2,\omega) is an orthonormal frame. The normal to their transformed tangent plane is nT=T−Tω/∣T−Tω∣n_T=T^{-T}\omega/|T^{-T}\omega|. The transformed normal column has height ∣Tω⋅nT∣=1/∣T−Tω∣|T\omega\cdot n_T|=1/|T^{-T}\omega| above that plane. The determinant is tangent area times this height, as follows by expressing the three columns in an orthonormal frame with last vector nTn_T. Its absolute value is ∣det⁡T∣|\det T|. Hence the tangent area factor is ∣det⁡T∣ ∣T−Tω∣|\det T|\,|T^{-T}\omega|, as used above. The surface-coordinate prerequisite then applies it to each chart.

References