Distorted Fourier transforms 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: How does a shell observation become a spectral transform? At one energy the observation is an amplitude on a surface; over an interval it is a function of momentum. The coarea factor connects these descriptions. In the free square example there are two momenta at each positive energy, whereas the bound-state example requires removing a discrete spectral component before claiming a norm identity.

Resolvent boundary values determine the continuous spectral measure. Their free Fourier traces can be assembled into measurable functions of momentum, giving two distorted Fourier transforms. These maps preserve the norm of the continuous part and turn the perturbed operator into multiplication by the free polynomial.

Retain the hypotheses of Limiting absorption and point spectrum: a real simply characteristic pp without invariant directions, a symmetric short-range differential VV, and the self-adjoint closure HH of p(D)+Vp(D)+V on S\mathcal S. Write

Σ=Z(p)∪A,Ω=R∖Σ,Mλ={ξ:p(ξ)=λ},g(ξ)=∣∇p(ξ)∣.(1) \Sigma=Z(p)\cup\mathcal A,\quad \Omega=\mathbb R\setminus\Sigma,\quad M_\lambda=\{\xi:p(\xi)=\lambda\},\quad g(\xi)=|\nabla p(\xi)|. \tag{1}

The exceptional-set proof shows that Σ\Sigma is closed and countable. Since n≥1n\geq1 and pp has no invariant directions, pp is nonconstant: a constant polynomial would be invariant under every translation. For λ∈Ω\lambda\in\Omega define the strongly continuous BB-valued forcing

a±(λ,f)=(I+VR0,±(λ))−1f,f∈B.(2) a_\pm(\lambda,f)=(I+VR_{0,\pm}(\lambda))^{-1}f, \qquad f\in B. \tag{2}

Use the unitary Fourier transform f^(ξ)=(2π)−n/2∫e−ix⋅ξf(x) dx\widehat f(\xi)=(2\pi)^{-n/2}\int e^{-ix\cdot\xi}f(x)\,dx. The canonical trace Tλ:B→L2(Mλ,dS)T_\lambda:B\to L^2(M_\lambda,dS) is from Global radiation and flux. An L2L^2 Fourier representative by itself need not have a meaningful restriction to MλM_\lambda.

The functional-analysis prerequisite is the full local proof of Self-adjoint spectral calculus with the original domain. Its Cayley-transform construction gives the unique projection-valued measure, bounded Borel calculus and every maximal unbounded multiplier domain for an arbitrary self-adjoint operator, with D(H)={f:∫t2 d(E(t)f,f)<∞}D(H)=\{f:\int t^2\,d(E(t)f,f)<\infty\} equal to the original operator domain. Neither a lower bound nor separability is assumed. Write EE for the spectral measure of HH and

μf(A)=(E(A)f,f),Ed=E(Σ),Ec=E(Ω).(3) \mu_f(A)=(E(A)f,f),\qquad E^d=E(\Sigma),\quad E^c=E(\Omega). \tag{3}

For further reading on spectral calculus and inversion, see Teschl [T], Sections 3.1 and 3.4; for perturbed spectral representations, see Kuroda [K], Section 4.2.

1. Continuous test functions recover the spectral measure

This is the continuous-test counterpart of the interval formula in Resolvents, domains and spectral density, Section 4. Here the explicit general spectral prerequisite permits a free polynomial and its perturbation to be unbounded in both spectral directions.

Lemma 1.1. For a self-adjoint HH, f∈L2f\in L^2 and real χ∈Cc(R)\chi\in C_c(\mathbb R),

∫χ dμf=lim⁡ε↓0±1π∫Rχ(λ)Im⁡((H−λ∓iε)−1f,f) dλ.(4) \int\chi\,d\mu_f =\lim_{\varepsilon\downarrow0} \frac{\pm1}{\pi}\int_{\mathbb R} \chi(\lambda)\operatorname{Im} ((H-\lambda\mp i\varepsilon)^{-1}f,f)\,d\lambda. \tag{4}

Proof. By the spectral theorem the signed imaginary part of the resolvent pairing is

±Im⁡((H−λ∓iε)−1f,f)=∫Rε(t−λ)2+ε2 dμf(t).(5) \pm\operatorname{Im}((H-\lambda\mp i\varepsilon)^{-1}f,f) =\int_{\mathbb R} \frac{\varepsilon}{(t-\lambda)^2+\varepsilon^2}\,d\mu_f(t). \tag{5}

The arctangent normalization gives integral π\pi for this kernel in λ\lambda, while μf(R)=∥f∥22\mu_f(\mathbb R)=\|f\|_2^2. The general product theorem, applied to the finite spectral measure and Lebesgue measure, therefore justifies Fubini even for a signed bounded χ\chi. The right side of (4) before its limit is ∫(Pε∗χ)(t) dμf(t)\int(P_\varepsilon*\chi)(t)\,d\mu_f(t), where Pε(s)=επ(s2+ε2). P_\varepsilon(s)=\frac{\varepsilon}{\pi(s^2+\varepsilon^2)}. This convolution converges uniformly to χ\chi. Indeed, given δ>0\delta>0, the integral on ∣s∣≤δ|s|\leq\delta is bounded by the uniform modulus of continuity ωχ(δ)\omega_\chi(\delta). On its complement the error is bounded by 2∥χ∥∞2\|\chi\|_\infty times a kernel mass at most 2ε/(πδ)2\varepsilon/(\pi\delta). Let ε→0\varepsilon\to0, then δ→0\delta\to0. Integration against the finite measure proves (4). □\square

Theorem 1.2. For f∈Bf\in B, the restriction of μf\mu_f to Ω\Omega has a nonnegative continuous density. With either sign it is

qf(λ)=∫Mλ∣Tλa±(λ,f)(ξ)∣2dS(ξ)g(ξ)=±1πIm⁡(RH,±(λ)f,f).(6) q_f(\lambda) =\int_{M_\lambda} |T_\lambda a_\pm(\lambda,f)(\xi)|^2 \frac{dS(\xi)}{g(\xi)} =\frac{\pm1}{\pi}\operatorname{Im}(R_{H,\pm}(\lambda)f,f). \tag{6}

Proof. Let a=a±(λ,f)a=a_\pm(\lambda,f) and u=R0,±(λ)a=RH,±(λ)fu=R_{0,\pm}(\lambda)a=R_{H,\pm}(\lambda)f. The factorization says f=a+Vuf=a+Vu. The extended endpoint symmetry makes (u,Vu)(u,Vu) real. Hence Im⁡(u,f)=Im⁡(u,a). \operatorname{Im}(u,f)=\operatorname{Im}(u,a). The free forcing-flux formula gives ±2Im⁡(u,a)=2π∫Mλ∣Tλa∣2 dS/g, \pm2\operatorname{Im}(u,a) =2\pi\int_{M_\lambda}|T_\lambda a|^2\,dS/g, proving the second equality in (6) and its nonnegativity. All pairings converge because u∈Xp⊂B∗u\in X_p\subset B^* and f,a,Vu∈Bf,a,Vu\in B. The threshold estimate bounds 1/g1/g uniformly on compact regular energy sets; consequently the integral is finite.

For real χ∈Cc(Ω)\chi\in C_c(\Omega), limiting absorption bounds RH(λ±iε):B→B∗R_H(\lambda\pm i\varepsilon):B\to B^* uniformly on a compact neighborhood of its support. It also makes the pairing on a fixed ff continuous up to the boundary. Dominated convergence in (4) therefore gives

∫χ dμf=∫Ωχ(λ)qf(λ) dλ.(7) \int\chi\,d\mu_f=\int_\Omega\chi(\lambda)q_f(\lambda)\,d\lambda. \tag{7}

The last expression in (6) is continuous, by weak-star continuity of the resolvent against f∈Bf\in B. Thus qfq_f is continuous. Equality for compact continuous tests determines the measures on Ω\Omega. Indeed, if (a,b)(a,b) has compact closure in Ω\Omega, the continuous tents min⁡{1,kdist⁡(λ,R∖(a,b))}\min\{1,k\operatorname{dist}(\lambda,\mathbb R\setminus(a,b))\} increase to its indicator. Monotone convergence gives equality on that interval. Fix such an interval UU; both measures are finite there, and the same argument gives equality on every relatively open subinterval of UU, including UU itself. These intervals are closed under finite intersections and generate the Borel sets of UU. The finite-measure uniqueness proof therefore gives equality on all those Borel sets. Countably many such intervals cover Ω\Omega; disjointifying that cover and using countable additivity gives the claim on all of Ω\Omega. Both signs in (7) give the same measure and continuous density, so their two expressions in (6) agree at every λ∈Ω\lambda\in\Omega. □\square

The factor 1/π1/\pi in (6) and the factor 2π2\pi in free flux cancel. The unitary Fourier convention leaves no further factor in the surface density.

Absolute continuity alone already follows from the locally uniform BB-to-B∗B^* bound in the limiting-absorption lesson. Indeed, on a compact good-energy interval it bounds the positive scalar imaginary part by CI∥f∥B2C_I\|f\|_B^2. Apply the positive-kernel proof of Theorem 5.1 in Resolvents, domains and spectral density, using the general spectral measure supplied here. It gives a density bounded locally by CI∥f∥B2/πC_I\|f\|_B^2/\pi; dense BB tests then exclude singular spectral mass on Ω\Omega. The boundary convergence and free flux calculation above give the stronger continuous density and its exact surface formula (6), which are needed for the spectral transform.

2. Choosing measurable functions with the prescribed surface traces

The change-of-variables and surface-measure facts used here are proved in Coordinate inverses and integration, CI1–CI7. We also give the parameter partition construction, including the extra time parameter used in the next lesson. For any open O⊂RdO\subset\mathbb R^d, put Km={x∈O:∣x∣≤m, dist⁡(x,Rd∖O)≥1/m},m≥1, K_m=\{x\in O:|x|\leq m,\ \operatorname{dist}(x,\mathbb R^d\setminus O)\geq1/m\}, \qquad m\geq1, with distance to the empty set taken as infinity and K0=K−1=∅K_0=K_{-1}=\varnothing. These compact sets exhaust OO and Km⊂int⁡Km+1K_m\subset\operatorname{int}K_{m+1}. Given an open cover of OO, cover each compact band Km∖int⁡Km−1K_m\setminus\operatorname{int}K_{m-1} by finitely many balls whose slightly larger closed balls lie in a member of that cover and in int⁡Km+1∖Km−2\operatorname{int}K_{m+1}\setminus K_{m-2}. This is possible because the band is disjoint from Km−2K_{m-2}. Choose a smooth nonnegative bump supported in each larger ball and positive on the corresponding smaller ball. The collection is locally finite: a neighborhood of any point lies in some int⁡KM\operatorname{int}K_M, and every sufficiently late ball avoids KMK_M. The sum of the bumps is smooth and strictly positive on OO. Dividing each bump by their sum gives the required locally finite smooth partition with supports subordinate to the cover. This proof works for d=1d=1 here and d=2d=2 for the energy-time parameter.

Lemma 2.1. For each f∈Bf\in B, there are measurable functions J±fJ_\pm f on Rn\mathbb R^n such that, for every λ∈Ω\lambda\in\Omega,

(J±f)∣Mλ=Tλa±(λ,f)for almost every point of Mλ.(8) (J_\pm f)|_{M_\lambda}=T_\lambda a_\pm(\lambda,f) \quad\text{for almost every point of }M_\lambda. \tag{8}

They are unique up to Lebesgue-null sets. They may be chosen to vanish on p−1(Σ)p^{-1}(\Sigma).

Proof. Fix a sign and write a(λ)=a±(λ,f)a(\lambda)=a_\pm(\lambda,f), a continuous BB-valued function on Ω\Omega. For each integer j≥1j\geq1, construct a continuous BB-valued approximation hj(λ)h_j(\lambda), each of whose values has smooth compact Fourier support, with ∥hj(λ)−a(λ)∥B≤2−j(λ∈Ω).(9) \|h_j(\lambda)-a(\lambda)\|_B\leq2^{-j} \quad(\lambda\in\Omega). \tag{9} Here are the details. Schwartz functions are dense in BB by the shell approximation in Endpoint spaces and flat energy shells. For a Schwartz function, multiplying its Fourier transform by a smooth cutoff equal to one on a ball of radius tending to infinity converges in every Schwartz seminorm, by Leibniz and its rapidly decreasing derivatives. Fourier inversion on Schwartz space preserves this convergence. For its implication in BB, let v∈Sv\in\mathcal S and choose an integer N>(n+1)/2N>(n+1)/2. The shell of radius Rℓ=2ℓR_\ell=2^\ell has volume at most CnRℓnC_nR_\ell^n. Thus its contribution to ∥v∥B\|v\|_B is at most Cn,NRℓ(n+1)/2−Nsup⁡x(1+∣x∣)N∣v(x)∣C_{n,N}R_\ell^{(n+1)/2-N}\sup_x(1+|x|)^N|v(x)|, with the inner shell bounded by the same seminorm. The geometric series converges. Consequently convergence in Schwartz space implies convergence in BB. Thus F−1Cc∞\mathcal F^{-1}C_c^\infty is dense in BB. At each parameter choose a member of this subspace approximating aa within 2−j−22^{-j-2}; continuity gives a neighborhood on which the same member approximates aa within 2−j−12^{-j-1}. Apply the partition construction above to this cover and take the corresponding convex combinations of the chosen members. Only finitely many terms occur near each parameter, so hjh_j is continuous in BB and its Fourier values are jointly measurable in (λ,ξ)(\lambda,\xi). Its value at every fixed parameter is a finite sum of smooth compact Fourier functions. Convexity proves (9).

Define on p−1(Ω)p^{-1}(\Omega) Gj(ξ)=hj(p(ξ))^(ξ).(10) G_j(\xi)=\widehat{h_j(p(\xi))}(\xi). \tag{10} Each GjG_j is continuous on p−1(Ω)p^{-1}(\Omega): locally in the energy parameter the defining partition has only finitely many terms, each a smooth parameter coefficient times a smooth Fourier function. In particular it is Borel measurable. On MλM_\lambda it is the classical Fourier restriction of hj(λ)h_j(\lambda). Uniform trace bounds on compact regular energy intervals give ∥Gj∣Mλ−Tλa(λ)∥L2(dS)≤CI2−j,λ∈I⋐Ω.(11) \|G_j|_{M_\lambda}-T_\lambda a(\lambda)\|_{L^2(dS)} \leq C_I2^{-j},\qquad \lambda\in I\Subset\Omega. \tag{11} In particular, for every fixed λ\lambda, the sum of the L2(dS)L^2(dS) norms of the successive differences is finite. The nonnegative sums ∑j∣Gj+1−Gj∣\sum_j|G_{j+1}-G_j| have uniformly bounded L2L^2 norm on that surface by the triangle inequality; monotone convergence makes the full sum finite almost everywhere. Thus GjG_j converges pointwise almost everywhere on each fixed MλM_\lambda. Its L2L^2 limit is the trace in (11).

Define J±f(ξ)J_\pm f(\xi) to be this pointwise limit wherever it exists in p−1(Ω)p^{-1}(\Omega), and zero elsewhere. The convergence set is Borel by the countable Cauchy criterion; on it the real and imaginary limits are Borel functions. Extending by zero is therefore Borel measurable. This proves (8) for every regular λ∈Ω\lambda\in\Omega, with a possibly different exceptional surface-null set for each λ\lambda.

For the assertion in ambient measure, use the local energy coordinate p(ξ)p(\xi) and the coarea formula ∫p−1(Ω)F(ξ) dξ=∫Ω∫MλF(ξ) dS(ξ)g(ξ) dλ,F≥0.(12) \int_{p^{-1}(\Omega)}F(\xi)\,d\xi =\int_\Omega\int_{M_\lambda}F(\xi)\,\frac{dS(\xi)}{g(\xi)}\,d\lambda, \qquad F\geq0. \tag{12} On compact subsets of p−1(Ω)p^{-1}(\Omega), gg is bounded below and (11) supplies the required locally integrable estimates. For the Jacobian explicitly, where ∂1p≠0\partial_1p\ne0 write ξ=(h(λ,η),η)\xi=(h(\lambda,\eta),\eta). Implicit differentiation gives ∂λh=(∂1p)−1\partial_\lambda h=(\partial_1p)^{-1} and ∇ηh=−∇ηp/∂1p\nabla_\eta h=-\nabla_\eta p/\partial_1p. Hence the energy-coordinate volume Jacobian is 1/∣∂1p∣1/|\partial_1p|, while dSg=1+∣∇ηh∣2∣∇p∣ dη=dη∣∂1p∣. \frac{dS}{g} =\frac{\sqrt{1+|\nabla_\eta h|^2}}{|\nabla p|}\,d\eta =\frac{d\eta}{|\partial_1p|}. The usual change-of-variables theorem in these coordinate patches and a partition of unity yield (12); increasing compact cutoffs gives its full nonnegative form. Surface-null sets at every λ\lambda therefore give an ambient null set. Also each level set p−1(λ)p^{-1}(\lambda) has Lebesgue measure zero, since p−λp-\lambda is a nonzero polynomial. Here is the full elementary argument. A nonzero one-variable polynomial of degree dd has at most dd distinct roots: if aa is a root, the identities tk−ak=(t−a)∑r=0k−1tk−1−rart^k-a^k=(t-a)\sum_{r=0}^{k-1}t^{k-1-r}a^r factor out t−at-a, and induction on the degree applies to the other roots. A nonzero constant has none. In dimension n>1n>1, write a nonzero polynomial as Q(x′,t)=∑r=0dcr(x′)trQ(x',t)=\sum_{r=0}^d c_r(x')t^r and select one coefficient which is not the zero polynomial. By dimension induction its zero set in x′x' is null. Outside this null set the section in tt has only finitely many roots. Inside any bounded box the zero-set indicator is Borel, its one-dimensional section integral is zero outside the exceptional set, and is bounded by the box length on that set. Tonelli therefore gives zero volume in the box. A countable union of boxes proves the claim in all of Rn\mathbb R^n. Since Σ\Sigma is countable, p−1(Σ)p^{-1}(\Sigma) is null. These facts prove uniqueness almost everywhere and justify the prescribed zero values. □\square

The construction uses canonical BB traces before taking a diagonal in energy and frequency. It preserves the statement on every good surface, as required for subsequent energy-by-energy formulas.

3. The norm identity and the two spectral parts

Theorem 3.1. For f∈Bf\in B, J±f∈L2(dξ)J_\pm f\in L^2(d\xi), and

∥J±f∥L2(dξ)2=∥Ecf∥22.(13) \|J_\pm f\|_{L^2(d\xi)}^2=\|E^cf\|_2^2. \tag{13}

The maps extend uniquely to bounded linear operators J±:L2(dx)⟶L2(dξ). J_\pm:L^2(dx)\longrightarrow L^2(d\xi). They vanish on EdL2E^dL^2 and are isometries on EcL2E^cL^2. The subspace EdL2E^dL^2 is the closed span of all L2L^2 eigenfunctions of HH, and EcL2=Hac(H)E^cL^2=\mathcal H_{\mathrm{ac}}(H). Thus HH has no singular continuous spectrum.

Proof. Exhaust Ω\Omega by compact continuous tests 0≤χk↑10\leq\chi_k\uparrow1. Such a sequence exists because Ω\Omega is open: use an increasing cutoff in ∣λ∣|\lambda| and an increasing cutoff which vanishes where dist⁡(λ,Σ)≤1/k\operatorname{dist}(\lambda,\Sigma)\leq1/k. If Σ\Sigma is empty, only the first cutoff is needed. Monotone convergence in (7), followed by (8) and (12), gives ∥Ecf∥22=∫Ωqf(λ) dλ=∫p−1(Ω)∣J±f(ξ)∣2 dξ. \|E^cf\|_2^2 =\int_\Omega q_f(\lambda)\,d\lambda =\int_{p^{-1}(\Omega)}|J_\pm f(\xi)|^2\,d\xi. The prescribed values off p−1(Ω)p^{-1}(\Omega) prove (13). Linearity holds as an L2L^2 identity: the canonical traces in (8) are linear in ff, so coarea and uniqueness identify the representative for a linear combination with that linear combination of representatives. Formula (13) bounds these operators by one on the dense subspace B⊂L2B\subset L^2, giving their unique bounded extensions. Passing to L2L^2 limits preserves (13), which shows that they vanish on EdL2E^dL^2 and preserve norm on its orthogonal complement.

For an atom λ\lambda, the spectral theorem gives E({λ})L2=ker⁡(H−λ).(14) E(\{\lambda\})L^2=\ker(H-\lambda). \tag{14} Indeed a vector in the range has spectral measure supported at λ\lambda, hence finite second moment and (H−λ)(H-\lambda) norm zero. Conversely an eigenvector has ∫∣t−λ∣2 dμ=0\int|t-\lambda|^2\,d\mu=0, so its measure is supported at λ\lambda and it belongs to that range. Countable strong additivity over Σ\Sigma now makes EdL2E^dL^2 the closed span of the corresponding eigenspaces. The previous lesson identifies every eigenvalue outside Z(p)Z(p) with A\mathcal A, so there are no eigenvectors in EcL2E^cL^2, and all eigenvectors belong to EdL2E^dL^2. No multiplicity or decay assertion at a threshold is needed here.

Let N⊂ΩN\subset\Omega be a Lebesgue-null Borel set. Theorem 1.2 gives ∥E(N)f∥22=μf(N)=0\|E(N)f\|_2^2=\mu_f(N)=0 for f∈Bf\in B. Density and boundedness of E(N)E(N) give E(N)=0E(N)=0 on all L2L^2. For u∈EcL2u\in E^cL^2, its spectral measure also vanishes on Σ\Sigma, so it is absolutely continuous on all of R\mathbb R. Conversely, for an absolutely continuous vector all singleton projections in the countable set Σ\Sigma vanish; strong additivity gives Edu=0E^du=0, so u∈EcL2u\in E^cL^2. On EdL2E^dL^2 all measures are countable sums of atoms. These orthogonal reducing parts sum to L2L^2, excluding a singular continuous component. □\square

Surjectivity of J±:EcL2→L2(dξ)J_\pm:E^cL^2\to L^2(d\xi) will follow from comparison with the time-dependent wave operators. The present theorem proves its isometric range without assuming that comparison.

4. Intertwining the closed operator and its group

Let MpM_p be the maximal self-adjoint multiplication operator by p(ξ)p(\xi) on L2(dξ)L^2(d\xi).

Theorem 4.1. For every f∈D(H)f\in\mathcal D(H),

J±f∈D(Mp),J±Hf=MpJ±f.(15) J_\pm f\in\mathcal D(M_p),\qquad J_\pm Hf=M_pJ_\pm f. \tag{15}

For every f∈L2f\in L^2 and t∈Rt\in\mathbb R,

J±eitHf=eitMpJ±f.(16) J_\pm e^{itH}f=e^{itM_p}J_\pm f. \tag{16}

Proof. First take ϕ∈S\phi\in\mathcal S. Its initial operator image Hϕ=p(D)ϕ+VϕH\phi=p(D)\phi+V\phi belongs to BB: the polynomial part is Schwartz, while ϕ∈Xp\phi\in X_p and the short-range compact extension sends it into BB. At a regular energy, R0,±(λ)(P0−λ)ϕ=ϕ.(17) R_{0,\pm}(\lambda)(P_0-\lambda)\phi=\phi. \tag{17} For example, the nonreal identity is R0(z)(P0−λ)ϕ=ϕ+(z−λ)R0(z)ϕR_0(z)(P_0-\lambda)\phi=\phi+(z-\lambda)R_0(z)\phi; free endpoint bounds make its last term tend to zero distributionally as z→λz\to\lambda on either side. Therefore (H−λ)ϕ=(I+VR0,±(λ))(P0−λ)ϕ. (H-\lambda)\phi =(I+VR_{0,\pm}(\lambda))(P_0-\lambda)\phi. Invert this equality at λ∈Ω\lambda\in\Omega. Taking the canonical Fourier trace of its right forcing gives Tλa±(λ,(H−λ)ϕ)=Tλ(P0−λ)ϕ=0. T_\lambda a_\pm(\lambda,(H-\lambda)\phi) =T_\lambda(P_0-\lambda)\phi=0. Linearity in the forcing, (8) and coarea imply J±Hϕ(ξ)=p(ξ)J±ϕ(ξ)almost everywhere.(18) J_\pm H\phi(\xi)=p(\xi)J_\pm\phi(\xi) \quad\text{almost everywhere}. \tag{18} In particular the right side is L2L^2, proving domain membership for these core vectors.

Since HH is the closure of its Schwartz restriction, any f∈D(H)f\in\mathcal D(H) has ϕj∈S\phi_j\in\mathcal S with ϕj→f\phi_j\to f and Hϕj→HfH\phi_j\to Hf in L2L^2. Boundedness of J±J_\pm gives J±ϕj→J±fJ_\pm\phi_j\to J_\pm f and MpJ±ϕj=J±Hϕj→J±HfM_pJ_\pm\phi_j=J_\pm H\phi_j\to J_\pm Hf. Multiplication by pp is closed. To see this directly, from vj→vv_j\to v and pvj→wpv_j\to w in L2L^2 choose a common subsequence with ∑k(∥vjk−v∥22+∥pvjk−w∥22)<∞\sum_k(\|v_{j_k}-v\|_2^2+\|pv_{j_k}-w\|_2^2)<\infty. Tonelli makes both error sums finite almost everywhere, so both errors tend to zero there. Since pp is finite everywhere, w=pvw=pv almost everywhere; hence vv belongs to the maximal domain and has the required image. This proves (15).

For the group statement, the bounded spectral calculus gives domain invariance and differentiation on D(H)\mathcal D(H). They also follow directly from the second moment: multiplication by eitλe^{it\lambda} preserves ∫λ2 dμf\int\lambda^2\,d\mu_f, and ∣(eihλ−1)/h∣≤∣λ∣|(e^{ih\lambda}-1)/h|\leq|\lambda| permits dominated convergence of the difference quotient. The same statements hold for MpM_p. For f∈D(H)f\in\mathcal D(H) differentiate the L2(dξ)L^2(d\xi)-valued function b(t)=e−itMpJ±eitHf. b(t)=e^{-itM_p}J_\pm e^{itH}f. To justify the product rule, write y(t)=J±eitHfy(t)=J_\pm e^{itH}f. It lies in D(Mp)D(M_p) by (15), has derivative iJ±HeitHfiJ_\pm He^{itH}f, and Mpy(t)=J±eitHHfM_py(t)=J_\pm e^{itH}Hf is continuous. Split the difference quotient of e−itMpy(t)e^{-itM_p}y(t) into the change of yy multiplied by the nearby unitary and the change of the unitary on the fixed domain vector y(t)y(t). Strong continuity treats the first term and the domain derivative treats the second. Its derivative is b′(t)=ie−itMp(J±H−MpJ±)eitHf=0. b'(t)=i e^{-itM_p}(J_\pm H-M_pJ_\pm)e^{itH}f=0. Thus b(t)=b(0)b(t)=b(0), proving (16) on the domain. This domain is dense, and both sides are bounded on L2L^2, so (16) holds for every ff. □\square

Example 4.2 (the free square). For V=0V=0 and p(ξ)=ξ2p(\xi)=\xi^2 in one dimension, J±=FJ_\pm=\mathcal F. Indeed a±(λ,f)=fa_\pm(\lambda,f)=f; for a Schwartz ff, (8) gives its ordinary Fourier values on every nonzero level. The excluded preimage of Σ\Sigma is null by Lemma 2.1, so the two functions agree in ambient measure. Since both operators are bounded and Schwartz functions are dense in L2L^2, the equality extends to every L2L^2 vector. The free measure is absolutely continuous, and for λ>0\lambda>0,

qf(λ)=∣f^(λ)∣2+∣f^(−λ)∣22λ.(19) q_f(\lambda)= \frac{|\widehat f(\sqrt\lambda)|^2+ |\widehat f(-\sqrt\lambda)|^2}{2\sqrt\lambda}. \tag{19}

For λ<0\lambda<0 it is zero. Each one-point surface branch has zero-dimensional measure one, and g(±λ)=2λg(\pm\sqrt\lambda)=2\sqrt\lambda, which explains both terms. The threshold zero is excluded from the continuity assertion.

Example 4.3 (removing a bound state). For H=−∂x2−2sech⁡2xH=-\partial_x^2-2\operatorname{sech}^2x, the bound-state example proves that u=sech⁡xu=\operatorname{sech}x is an eigenfunction at −1-1. Its exponential bound ∣u(x)∣≤2e−∣x∣|u(x)|\leq2e^{-|x|} also puts it in BB: in one dimension the outer shell contribution is at most CRje−Rj/2C R_j e^{-R_j/2}, a summable sequence. For example et≥t3/6e^t\geq t^3/6 bounds this by CRj−2C R_j^{-2}; the inner shell is finite. Thus Ecu=0E^cu=0 and J±u=0J_\pm u=0 in L2(dξ)L^2(d\xi). More precisely, (6) is continuous and nonnegative on Ω\Omega, and its integral for this uu is zero. It vanishes at every λ∈Ω\lambda\in\Omega, so Tλa±(λ,u)=0T_\lambda a_\pm(\lambda,u)=0 on every such surface. The distorted transform removes the bound state by its energy-dependent forcing correction.

Use the conclusion

Check the sign of the Poisson formula, then the exact coarea density and the operator domain after transformation. The transform's isometry on continuous states is a separate step from the later onto theorem.

5. Exercises

Exercise 5.1 (foundation). Prove the estimate ∥Pε∗χ−χ∥∞≤ωχ(δ)+4ε∥χ∥∞πδ. \|P_\varepsilon*\chi-\chi\|_\infty \leq\omega_\chi(\delta) +\frac{4\varepsilon\|\chi\|_\infty}{\pi\delta}. Explain why it applies uniformly even when the centre tt lies outside the support of χ\chi.

Exercise 5.2 (foundation). Derive (19) by changing variables separately on ξ>0\xi>0 and ξ<0\xi<0. Verify that its integral over λ>0\lambda>0 equals ∥f∥22\|f\|_2^2.

Exercise 5.3 (intermediate). Let p(ξ)=ξp(\xi)=\xi on R\mathbb R, and let hλ(ξ)=1{λ}(ξ)h_\lambda(\xi)=\mathbf1_{\{\lambda\}}(\xi). Each hλh_\lambda represents zero in ambient L2(dξ)L^2(d\xi). Compute hp(ξ)(ξ)h_{p(\xi)}(\xi), and explain why this choice cannot represent the canonical trace of the zero BB function.

Exercise 5.4 (intermediate). Suppose f∈Bf\in B and Ecf=0E^cf=0. Prove that both Tλa±(λ,f)T_\lambda a_\pm(\lambda,f) vanish for every λ∈Ω\lambda\in\Omega, not merely almost every energy.

Exercise 5.5 (advanced). Let A,BA,B be self-adjoint operators on Hilbert spaces and JJ a bounded map. Suppose a core C\mathcal C of AA satisfies JC⊂D(B)J\mathcal C\subset\mathcal D(B) and JAϕ=BJϕJA\phi=BJ\phi. Prove domain intertwining and group intertwining on the full spaces, carefully justifying the product derivative.

6. Complete solutions

Solution 5.1. Write the difference as ∫Pε(s)[χ(t−s)−χ(t)] ds\int P_\varepsilon(s)[\chi(t-s)-\chi(t)]\,ds. The integrand difference is at most ωχ(δ)\omega_\chi(\delta) when ∣s∣≤δ|s|\leq\delta. On the complement it is at most 2∥χ∥∞2\|\chi\|_\infty, and ∫∣s∣>δPε(s) ds≤2επ∫δ∞s−2 ds=2επδ. \int_{|s|>\delta}P_\varepsilon(s)\,ds \leq\frac{2\varepsilon}{\pi}\int_\delta^\infty s^{-2}\,ds =\frac{2\varepsilon}{\pi\delta}. The kernel's total mass is one. Adding the two bounds proves the estimate. Neither bound depends on tt; compact continuous χ\chi, extended by zero, is uniformly continuous on the whole line, so the estimate also applies outside its support.

Solution 5.2. On the positive branch λ=ξ2\lambda=\xi^2 has dξ=dλ/(2λ)d\xi=d\lambda/(2\sqrt\lambda); on the negative branch the same absolute Jacobian appears. Therefore, for a nonnegative test χ\chi, ∫χ(ξ2)∣f^(ξ)∣2 dξ=∫0∞χ(λ)∣f^(λ)∣2+∣f^(−λ)∣22λ dλ. \int\chi(\xi^2)|\widehat f(\xi)|^2\,d\xi =\int_0^\infty\chi(\lambda) \frac{|\widehat f(\sqrt\lambda)|^2+ |\widehat f(-\sqrt\lambda)|^2} {2\sqrt\lambda}\,d\lambda. The change of variables is first valid away from zero and then on the full half-lines by monotone convergence. This gives (19). Taking tests increasing to one gives the full Fourier L2L^2 norm, equal to ∥f∥22\|f\|_2^2 by unitary Plancherel.

Solution 5.3. At every ξ\xi, hp(ξ)(ξ)=hξ(ξ)=1h_{p(\xi)}(\xi)=h_\xi(\xi)=1. Thus an arbitrary selection of ambient null-set representatives can create the constant function one on the energy-frequency diagonal, despite representing zero for every fixed parameter in ambient L2L^2. Here Mλ={λ}M_\lambda=\{\lambda\} has zero-dimensional surface measure one. The canonical trace of the zero BB function is zero on that point, while the selected representative has value one. Lemma 2.1 first approximates in BB and controls the actual trace on each surface; an arbitrary ambient representative has no such control.

Solution 5.4. By (13), ∫Ωqf(λ) dλ=0\int_\Omega q_f(\lambda)\,d\lambda=0. The nonnegative continuous function qfq_f must vanish everywhere in the open set Ω\Omega: a positive value would remain bounded below by a positive number on a small interval and give a positive integral. Formula (6) now gives a zero integral of ∣Tλa±∣2/g|T_\lambda a_\pm|^2/g for every λ\lambda. On each compact surface patch gg is positive and finite, so a zero weighted integral forces the trace to vanish almost everywhere there. A countable patch exhaustion gives its vanishing on the whole surface for each fixed energy and both signs.

Solution 5.5. Given f∈D(A)f\in\mathcal D(A), choose ϕj∈C\phi_j\in\mathcal C with ϕj→f\phi_j\to f and Aϕj→AfA\phi_j\to Af. Boundedness gives Jϕj→JfJ\phi_j\to Jf and BJϕj=JAϕj→JAfBJ\phi_j=JA\phi_j\to JAf. Closedness of BB proves Jf∈D(B)Jf\in\mathcal D(B) and BJf=JAfBJf=JAf.

For f∈D(A)f\in\mathcal D(A), the curve y(t)=JeitAfy(t)=Je^{itA}f takes values in D(B)\mathcal D(B), is differentiable as a Hilbert-space curve with y′(t)=iJAeitAfy'(t)=iJAe^{itA}f, and satisfies By(t)=JAeitAfBy(t)=JAe^{itA}f. The latter is continuous in tt, since AeitAf=eitAAfAe^{itA}f=e^{itA}Af. To differentiate e−itBy(t)e^{-itB}y(t), split its difference quotient into the change in yy, multiplied by the nearby unitary, and the change in that unitary on the fixed vector y(t)∈D(B)y(t)\in\mathcal D(B). Strong continuity handles the first term and the domain derivative handles the second. The result is e−itB[y′(t)−iBy(t)]=0e^{-itB}[y'(t)-iBy(t)]=0. Hence e−itBJeitAf=Jfe^{-itB}Je^{itA}f=Jf. Dense-domain approximation, boundedness of JJ, and unitarity extend the identity to all ff.

References