Energy-shell factors and outgoing equations

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

Working question: Why does factorizing an energy graph leave an outgoing choice? A first-order factor has homogeneous solutions as well as a particular solution. The outgoing condition supplies the direction of integration, just as in the transport model, but here the solution may initially have only endpoint ball bounds. The proof has to recover genuine Hilbert-space slices rather than assuming them at the start.

Near a regular energy surface, one frequency is a smooth function of the others. A small long-range perturbation moves that graph. Dividing by the graph factor suggests a first-order equation in one spatial coordinate. Such an equation still has many homogeneous solutions. An outgoing condition selects one of them, even when the initial information permits squared mass to grow proportionally to the radius of a ball.

We first prove the outgoing evolution estimate directly from its kernel. We then construct the perturbed energy graph and control every derivative of its polynomial factor and reciprocal. Read Endpoint spaces and flat energy shells for the spatial norms, and Admissible differential perturbations for the differential coefficient class. The measure and Euclidean product proofs give the scalar integral, convergence and slice facts; the same reading proves Fourier inversion and Plancherel. The coordinate inverse and bump proofs supply the regular free graph and smooth cutoffs. Section 3 constructs the vector integrals used here. Freely accessible background references are Teschl's evolution and initial-value texts [T, ODE] and Yafaev's Lectures on scattering theory [Y]. The evolution, rough-slice recovery, root and quotient proofs are given below.

1. What an outgoing condition must control

Write x=(s,z)∈R×Rdx=(s,z)\in\mathbb R\times\mathbb R^d, where d=n−1d=n-1, and put H=L2(Rd)\mathcal H=L^2(\mathbb R^d). For d=0d=0, use H=C\mathcal H=\mathbb C. Our convention is Ds=−i∂sD_s=-i\partial_s. Set

rj=2j,S0={∣x∣<1},Sj={2j−1≤∣x∣<2j}(j≥1),∥f∥B=∑j≥0rj1/2∥f∥L2(Sj),∥v∥B∗=sup⁡j≥0rj−1/2∥v∥L2(Sj).(1) \begin{gathered} r_j=2^j,\qquad S_0=\{|x|<1\},\\ S_j=\{2^{j-1}\le |x|<2^j\}\quad(j\ge1),\\ \|f\|_B=\sum_{j\ge0}r_j^{1/2}\|f\|_{L^2(S_j)},\\ \|v\|_{B^*}=\sup_{j\ge0}r_j^{-1/2}\|v\|_{L^2(S_j)}. \end{gathered} \tag{1}

Shells and balls give the equivalent bound

∫∣x∣<R∣v∣2 dx≤4R∥v∥B∗2,R≥1.(2) \int_{|x|<R}|v|^2\,dx \le 4R\|v\|_{B^*}^2,\qquad R\ge1. \tag{2}

Indeed, choose 2j−1<R≤2j2^{j-1}<R\le2^j and sum the shell estimates through jj. The sum of the radii is at most 2j+1≤4R2^{j+1}\le4R.

The outgoing condition in the ss direction is

lim⁡R→∞1R∫∣x∣<Rs<T∣v(x)∣2 dx=0for every fixed T∈R.(3) \begin{gathered} \lim_{R\to\infty} \frac1R\int_{\substack{|x|<R\\ s<T}}|v(x)|^2\,dx=0 \\ \quad\text{for every fixed }T\in\mathbb R. \end{gathered} \tag{3}

It excludes a wave already present at s=−∞s=-\infty. It permits a nonzero wave escaping toward s=+∞s=+\infty. Exercise 1 shows both distinctions in an exact model.

Two slice estimates will be useful:

∫R∥f(s)∥H ds≤2∥f∥B,∥v∥B∗≤2∥v∥L∞(R;H).(4) \begin{aligned} \int_{\mathbb R}\|f(s)\|_{\mathcal H}\,ds &\le\sqrt2\|f\|_B,\\ \|v\|_{B^*} &\le\sqrt2\|v\|_{L^\infty(\mathbb R;\mathcal H)}. \end{aligned} \tag{4}

For the first, each f1Sjf1_{S_j} has slices supported in [−rj,rj][-r_j,r_j]. Cauchy–Schwarz in ss bounds its slice integral by (2rj)1/2∥f∥L2(Sj)(2r_j)^{1/2}\|f\|_{L^2(S_j)}; sum over jj. For the second, integrate the uniform slice bound on the same interval.

2. A real symbol and its adjoint defect

Let b(s,z,η)b(s,z,\eta) be real and smooth, with support in one fixed compact set of η∈Rd\eta\in\mathbb R^d. Suppose every η\eta derivative is uniformly bounded. Suppose also that there is a nonnegative w∈L1(R)w\in L^1(\mathbb R) such that

∣∂ηα∇zb(s,z,η)∣≤Cαw(s).(5) |\partial_\eta^\alpha\nabla_z b(s,z,\eta)| \le C_\alpha w(s). \tag{5}

Assume its left-quantized operator A(s)=b(s,z,Dz)A(s)=b(s,z,D_z) is continuous in operator norm as a function of ss. Uniform bounds on ∂s∂ηαb\partial_s\partial_\eta^\alpha b, for all finite orders, are a sufficient condition, as the kernel estimate below shows. In dimension zero b(s)b(s) is a real bounded continuous scalar and the adjoint defect is zero.

For d≥1d\ge1 the kernel is

Ks(z,y)=(2π)−d∫ei(z−y)⋅ηb(s,z,η) dη.(6) K_s(z,y)=(2\pi)^{-d}\int e^{i(z-y)\cdot\eta}b(s,z,\eta)\,d\eta. \tag{6}

Integration by parts in η\eta, with compact support, gives for every integer N≥0N\ge0

∣Ks(z,y)∣≤CN(1+∣z−y∣)−N.(7) |K_s(z,y)|\le C_N(1+|z-y|)^{-N}. \tag{7}

For example, applying (1−Δη)k(1-\Delta_\eta)^k to the amplitude yields a factor (1+∣z−y∣2)−k(1+|z-y|^2)^{-k}. Choose 2k≥N2k\ge N.

We will use the following kernel inequality. If

sup⁡z∫∣K(z,y)∣ dy≤M1,sup⁡y∫∣K(z,y)∣ dz≤M2,(8) \begin{gathered} \sup_z\int |K(z,y)|\,dy\le M_1,\\ \qquad \sup_y\int |K(z,y)|\,dz\le M_2, \end{gathered} \tag{8}

then its L2L^2 operator norm is at most (M1M2)1/2(M_1M_2)^{1/2}. To prove it, apply weighted Cauchy–Schwarz to ∫K(z,y)h(y) dy\int K(z,y)h(y)\,dy, then integrate in zz and use Fubini. Thus (7), with N>dN>d, bounds A(s)A(s) uniformly.

Reality of bb identifies the kernel of A(s)−A(s)∗A(s)-A(s)^* as the inverse Fourier integral of b(s,z,η)−b(s,y,η)b(s,z,\eta)-b(s,y,\eta). The segment formula and (5), including its differentiated versions, bound that amplitude's η\eta derivatives by Cαw(s)∣z−y∣C_\alpha w(s)|z-y|. Repeating the integration by parts gives

∣Ks(z,y)−Ks(y,z)‾∣≤CNw(s)∣z−y∣(1+∣z−y∣)−N,∥A(s)−A(s)∗∥H→H≤Cw(s).(9) \begin{aligned} &|K_s(z,y)-\overline{K_s(y,z)}|\\ &\le C_N w(s)|z-y|(1+|z-y|)^{-N},\\ &\|A(s)-A(s)^*\|_{\mathcal H\to\mathcal H}\\ &\le Cw(s). \end{aligned} \tag{9}

Take N>d+1N>d+1 for the second line. A real left symbol need not give a self-adjoint operator. What we have proved is an integrable bound on its failure of self-adjointness.

3. Constructing the outgoing evolution

Theorem 3.1. Under the hypotheses of Section 2, let g∈L1(R;H)g\in L^1(\mathbb R;\mathcal H). There is a solution of

(Ds−A(s))v=g(10) (D_s-A(s))v=g \tag{10}

with a continuous, locally absolutely continuous H\mathcal H-valued representative, and

∥v(s)∥H≤C∫−∞s∥g(t)∥H dt.(11) \|v(s)\|_{\mathcal H} \le C\int_{-\infty}^s\|g(t)\|_{\mathcal H}\,dt. \tag{11}

Every distributional solution v∈B∗v\in B^* of (10) satisfying (3) has this representative and bound. In particular it is unique among the B∗B^* solutions satisfying (3). The constant depends on the symbol bounds and ∫w\int w, and is independent of gg.

Here are the integration facts needed for the full L1L^1 forcing class. The space H=L2(Rd)\mathcal H=L^2(\mathbb R^d) is separable: the Euclidean rectangle-approximation proof makes step functions on finitely many rational boxes with rational real and imaginary values dense. In dimension zero use C\mathbb C. A strongly measurable H\mathcal H-valued function with integrable norm is approximable in integral norm by finite-valued simple functions. Explicitly, first restrict to a bounded time interval and bounded values, losing arbitrarily little norm integral. Cover its values by balls of radius 1/j1/j centered in a countable dense subset, assign each value to the first such ball, and then keep finitely many of the resulting measurable sets. The discarded tail has vanishing measure; the values on the retained interval are bounded. This gives finite simple approximations with integral error tending to zero.

For a simple function F=∑khk1EkF=\sum_k h_k1_{E_k} on disjoint finite-measure sets define ∫F=∑k∣Ek∣hk\int F=\sum_k |E_k|h_k. The triangle inequality gives ∥∫F∥≤∫∥F∥\|\int F\|\le\int\|F\|, so completion defines the integral for the preceding L1L^1 class and makes it independent of the approximation. Pairing with a fixed vector, or applying a bounded operator, commutes with the integral by the same norm estimate. In particular the primitive G(s)=∫asF(t) dtG(s)=\int_a^s F(t)\,dt is locally absolutely continuous: sums of its increments over disjoint intervals are bounded by the norm integral over their union, which tends to zero with that union's measure by truncating the integrable scalar function ∥F∥\|F\|. Testing against a compact smooth scalar function and using scalar Fubini proves G′=FG'=F as a vector distribution. These facts also construct improper integrals whenever the norm has an integrable tail.

For continuous operator-valued functions on compact intervals use Riemann sums in operator norm, as in the continuous-vector integral proof. The same uniform-continuity estimate applies in the complete normed space of bounded operators: an operator-norm Cauchy sequence has a pointwise vector limit which is linear and bounded, and convergence to it is in operator norm. Thus no operator measurability assumption beyond the stated norm continuity is required.

We begin with existence. Since sup⁡s∥A(s)∥≤M\sup_s\|A(s)\|\le M, iteration of U(s,t)=I+i∫tsA(r)U(r,t) dr(12) U(s,t)=I+i\int_t^s A(r)U(r,t)\,dr \tag{12}

converges in operator norm on every bounded time interval. Its term of order kk has norm at most Mk∣s−t∣k/k!M^k|s-t|^k/k!. The same construction works for either orientation of the integral. Iterating the inequality for the difference of two solutions with zero initial value proves uniqueness of (12). Consequently

U(s,t)U(t,r)=U(s,r),U(t,s)=U(s,t)−1.(13) \begin{gathered} U(s,t)U(t,r)=U(s,r),\\ \qquad U(t,s)=U(s,t)^{-1}. \end{gathered} \tag{13}

For u(s)=U(s,t)hu(s)=U(s,t)h, its equation is ∂su=iA(s)u\partial_su=iA(s)u. The inner product is linear in the first variable. Equation (9) gives

∣dds∥u(s)∥2∣≤Cw(s)∥u(s)∥2.(14) \left|\frac{d}{ds}\|u(s)\|^2\right| \le Cw(s)\|u(s)\|^2. \tag{14}

For ε>0\varepsilon>0, divide (14) by ∥u(s)∥2+ε\|u(s)\|^2+\varepsilon. The absolute derivative of log⁡(∥u(s)∥2+ε)\log(\|u(s)\|^2+\varepsilon) is at most Cw(s)Cw(s). Integrate between the two times in either order, exponentiate and let ε\varepsilon decrease to zero. This proves

∥U(s,t)∥≤exp⁡ ⁣(C2∫Rw(r) dr)=:C0(s,t∈R).(15) \begin{gathered} \|U(s,t)\| \\ \le \exp\!\left(\frac C2\int_{\mathbb R}w(r)\,dr\right) \\ =:C_0 \\ \quad(s,t\in\mathbb R). \end{gathered} \tag{15}

Define

v(s)=i∫−∞sU(s,t)g(t) dt,q(s)=∫−∞s∥g(t)∥ dt.(16) \begin{aligned} v(s)&=i\int_{-\infty}^s U(s,t)g(t)\,dt,\\ q(s)&=\int_{-\infty}^s\|g(t)\|\,dt. \end{aligned} \tag{16}

The integral converges in H\mathcal H by (15), and ∥v(s)∥≤C0q(s)\|v(s)\|\le C_0q(s). Its integrand is strongly measurable: approximate gg by simple functions; on every bounded time interval the operator factor is norm continuous, and each resulting vector factor is strongly measurable. Splitting at a finite time aa and using (13) gives v(s)=U(s,a)(v(a)+i∫asU(a,t)g(t) dt). v(s)=U(s,a)\left(v(a)+i\int_a^s U(a,t)g(t)\,dt\right). The bracket is locally absolutely continuous by the primitive construction, and U(s,a)U(s,a) is continuously differentiable in operator norm. The product is therefore locally absolutely continuous. Testing its product rule, or first using simple approximations to gg, gives ∂sv=iAv+ig\partial_sv=iAv+ig in distributions, including the sign in (10). The same formula gives the inhomogeneous integral equation on each finite interval.

The bound and (4) give v∈B∗v\in B^*. To prove (3), fix S<min⁡(T,0)S<\min(T,0). For RR large enough,

1R∫∣x∣<Rs<T∣v∣2 dx≤C02(q(S)2+T−SRq(T)2).(17) \begin{gathered} \frac1R\int_{\substack{|x|<R\\s<T}}|v|^2\,dx \\ \le C_0^2\left(q(S)^2+ \frac{T-S}{R}q(T)^2\right). \end{gathered} \tag{17}

The first term bounds the interval [−R,S][-R,S]; omitting a nonexistent interval only reduces it. Let R→∞R\to\infty, then S→−∞S\to-\infty. Since q(S)→0q(S)\to0, this proves the condition.

4. Recovering slices from a rough outgoing solution

We now prove the remaining assertion of Theorem 3.1. A function in B∗B^* need not have globally square-integrable slices. The proof must produce that conclusion.

The shell estimate implies

∫Rd+1(1+∣x∣)−2p∣v(x)∣2 dx≤Cp∥v∥B∗2,p>12.(18) \begin{gathered} \int_{\mathbb R^{d+1}}(1+|x|)^{-2p}|v(x)|^2\,dx \\ \le C_p\|v\|_{B^*}^2,\\ \qquad p>\tfrac12. \end{gathered} \tag{18}

Indeed, the tail sum is bounded by a constant times ∑j2j(1−2p)\sum_j2^{j(1-2p)}. Fubini therefore gives polynomially weighted L2L^2 slices for almost every ss. Their strong measurability also follows from the rectangle proof: approximate the weighted function by finite product step functions in the space-time L2L^2 norm and take a subsequence whose slice errors tend to zero almost everywhere. On bounded time intervals the weight is equivalent to a weight in zz alone. The rapidly decreasing kernel (7) acts on these slices as a distribution and as a locally square-integrable function. Indeed, for bounded output zz, Cauchy–Schwarz bounds its value by the weighted slice norm times the Ly2L^2_y norm of Ks(z,y)(1+∣y∣)pK_s(z,y)(1+|y|)^p, uniformly finite on bounded output sets by (7). Integration in bounded s,zs,z proves local square integrability.

Choose a real smooth χ\chi, between zero and one, equal to one on the unit ball and zero outside the ball of radius two. Put χR(x)=χ(x/R)\chi_R(x)=\chi(x/R) and vR=χRvv_R=\chi_Rv. The distributional product rule gives

(Ds−A(s))vR=χRg+eR,eR=−iR−1(∂sχ)(x/R)v−[A(s),χR]v.(19) \begin{aligned} (D_s-A(s))v_R&=\chi_Rg+e_R,\\ e_R&=-iR^{-1}(\partial_s\chi)(x/R)v \\ &\quad-[A(s),\chi_R]v. \end{aligned} \tag{19}

Define the finite quantity M(r,T)=∥1{∣x∣<r, s<T}v∥2.(20) M(r,T)=\|1_{\{|x|<r,\ s<T\}}v\|_2. \tag{20}

By (3), M(r,T)=o(r1/2)M(r,T)=o(r^{1/2}) for each fixed TT. The first error in (19), using Cauchy–Schwarz over a time interval of length at most 4R4R, has the bound

∥R−1(∂sχ)(x/R)v∥L1((−∞,T);H)≤CR−1/2M(2R,T)=o(1).(21) \begin{gathered} \|R^{-1}(\partial_s\chi)(x/R)v\| _{L^1((-\infty,T);\mathcal H)} \\ \le CR^{-1/2}M(2R,T)=o(1). \end{gathered} \tag{21}

The commutator kernel is Ks(z,y)(χR(s,y)−χR(s,z)).(22) K_s(z,y)\bigl(\chi_R(s,y)-\chi_R(s,z)\bigr). \tag{22}

It is zero when ∣s∣≥2R|s|\ge2R. The Lipschitz estimate for χ\chi and (7) bound its absolute value by CNR−1∣z−y∣(1+∣z−y∣)−NC_NR^{-1}|z-y|(1+|z-y|)^{-N}. Thus (8) bounds its Lz2L^2_z operator norm by C/RC/R.

First apply it to the near input 1{∣x∣≤4R}v1_{\{|x|\le4R\}}v. Integrating its slice norm for s<Ts<T gives

∥[A,χR](1{∣x∣≤4R}v)∥L1((−∞,T);H)≤CR−1/2M(4R,T)=o(1).(23) \begin{gathered} \|[A,\chi_R](1_{\{|x|\le4R\}}v)\| _{L^1((-\infty,T);\mathcal H)} \\ \le CR^{-1/2}M(4R,T)=o(1). \end{gathered} \tag{23}

The far input requires a separate estimate. Decompose it into annuli rj<∣(s,y)∣≤2rjr_j<|(s,y)|\le2r_j, with rj=4R 2jr_j=4R\,2^j. For a nonzero kernel output, ∣(s,z)∣≤2R|(s,z)|\le2R, whereas the input cutoff is zero. Also ∣s∣≤2R|s|\le2R, so ∣z−y∣≥crj(24) |z-y|\ge c r_j \tag{24}

for a fixed c>0c>0. For example, ∣y∣≥(rj2−4R2)1/2≥3 rj/2|y|\ge(r_j^2-4R^2)^{1/2}\ge\sqrt3\,r_j/2 and ∣z∣≤2R≤rj/2|z|\le2R\le r_j/2.

For each fixed ss, both integrals in (8), restricted to these input and output sets, are at most CNrjd−NC_Nr_j^{d-N}. For the input integral, its volume is at most CrjdCr_j^d; for the output integral, it is at most CRd≤CrjdCR^d\le Cr_j^d. Equation (24) bounds the kernel by CNrj−NC_Nr_j^{-N}. The jj-th far term is consequently bounded in the slice L1L^1 norm by

CR1/2rjd−N∥v1{rj<∣x∣≤2rj}∥2≤CR1/2rjd+1/2−N∥v∥B∗.(25) \begin{aligned} &C R^{1/2}r_j^{d-N} \|v1_{\{r_j<|x|\le2r_j\}}\|_2 \\ &\le C R^{1/2}r_j^{d+1/2-N}\|v\|_{B^*}. \end{aligned} \tag{25}

Here (2) bounds the annular L2L^2 norm by Crj1/2∥v∥B∗C r_j^{1/2}\|v\|_{B^*}. Choose N=d+2N=d+2 and sum the geometric series. Its total is at most CR−1∥v∥B∗CR^{-1}\|v\|_{B^*}. This also proves convergence of the kernel action on the far input. In dimension zero the commutator is identically zero. We have proved, in every dimension,

∥eR∥L1((−∞,T);H)⟶0.(26) \|e_R\|_{L^1((-\infty,T);\mathcal H)}\longrightarrow0. \tag{26}

Now vRv_R is globally L2L^2 and has compact support. Equation (19) implies ∂svR=iAvR+i(χRg+eR)\partial_sv_R=iAv_R+i(\chi_Rg+e_R) in local L1(R;H)L^1(\mathbb R;\mathcal H). It therefore has a locally absolutely continuous representative, as follows directly. Subtract the vector primitive of the right side constructed in Section 3, leaving a locally integrable vector distribution FF with F′=0F'=0. Choose ρ∈Cc∞(R)\rho\in C_c^\infty(\mathbb R), ∫ρ=1\int\rho=1, and set h=∫ρ(s)F(s) ds∈Hh=\int\rho(s)F(s)\,ds\in\mathcal H. Every compact smooth test φ\varphi has ψ=φ−(∫φ)ρ\psi=\varphi-(\int\varphi)\rho of integral zero, so Ψ(s)=∫−∞sψ(t) dt\Psi(s)=\int_{-\infty}^s\psi(t)\,dt is compact and smooth. The identity F′=0F'=0 gives ∫Fψ=0\int F\psi=0, hence ∫Fφ=h∫φ\int F\varphi=h\int\varphi. Thus F=hF=h as a vector distribution and almost everywhere. For the last conclusion pair with a countable dense subset of H\mathcal H; scalar convolution with compact smooth approximate identities and their local L1L^1 convergence, proved in Approximation and convolution, identifies each scalar function almost everywhere. The common null set and density identify the vector. Adding the primitive back gives the claimed representative.

This representative is zero for s<−2Rs<-2R, since it is continuous and agrees there almost everywhere with zero. Variation of constants is valid for this rough forcing as well: multiplying the integral equation by U(a,s)U(a,s), or approximating its L1L^1 right side by simple functions, gives the formula in Section 3. Uniqueness follows by iterating the homogeneous norm inequality, with bound Mk∣s−a∣k/k!M^k|s-a|^k/k! at its kk-th iteration. Thus (15) implies, for s≤Ts\le T,

∥vR(s)∥≤C0q(s)+C0∥eR∥L1((−∞,T);H).(27) \|v_R(s)\|\le C_0q(s)+ C_0\|e_R\|_{L^1((-\infty,T);\mathcal H)}. \tag{27}

Take RR and the upper bound TT through the positive integers. Outside one common null set, vR(s,z)=χR(s,z)v(s,z)v_R(s,z)=\chi_R(s,z)v(s,z). Fatou's lemma in zz, followed by (26), proves ∥v(s)∥H≤C0q(s)\|v(s)\|_{\mathcal H}\le C_0q(s) for almost every ss. The slices now belong to H\mathcal H and are uniformly bounded. Equation (10) gives a locally absolutely continuous representative, and continuity extends the inequality to every ss.

For g=0g=0 this inequality forces v=0v=0. The difference of two B∗B^* solutions satisfying (3) again satisfies (3), by the squared triangle inequality. Thus uniqueness follows. Solutions obeying (11) satisfy (3) by (17), so the existence statement is unique in its asserted class as well. This completes Theorem 3.1. □\square

5. Moving a regular energy graph

Let P(ξ)P(\xi) be a real elliptic polynomial of degree mm. Fix an energy λ\lambda and a regular graph patch

ξ=(ξ1,η),P(E(η),η)=λ,∂ξ1P>0.(28) \begin{gathered} \xi=(\xi_1,\eta),\\ \qquad P(E(\eta),\eta)=\lambda,\\ \qquad \partial_{\xi_1}P>0. \end{gathered} \tag{28}

Work on a compact set of η\eta contained in an open graph patch. Shrink a collar around the graph so that the positive derivative has a fixed lower bound there. Let L(x,ξ)L(x,\xi) be a real polynomial in ξ\xi, of total degree at most mm, with smooth coefficients ℓγ(x)\ell_\gamma(x).

Fix 0<δ<1/30<\delta<1/3, put X=(1+∣x∣2)1/2X=(1+|x|^2)^{1/2}, and define the continuous concave function

μ(t)={δ+t,0≤t≤2,1+(1+δ)t/2,t≥2.(29) \mu(t)= \begin{cases} \delta+t,&0\le t\le2,\\ 1+(1+\delta)t/2,&t\ge2. \end{cases} \tag{29}

Assume the full bounds ∣∂xβℓγ(x)∣≤Cγ,βX−μ(∣β∣)(30) |\partial_x^\beta\ell_\gamma(x)| \le C_{\gamma,\beta}X^{-\mu(|\beta|)} \tag{30}

for every β\beta. Suppose also that L=0L=0 on ∣x∣<R0|x|<R_0, where R0R_0 may be chosen large with the constants in (30) fixed.

Theorem 5.1. For R0R_0 sufficiently large, there is a unique root ξ1=a(x,η)\xi_1=a(x,\eta) near E(η)E(\eta) of P(ξ)+L(x,ξ)=λP(\xi)+L(x,\xi)=\lambda. It is smooth on a neighborhood of the compact patch and satisfies, for every α,β\alpha,\beta,

∣∂ηα∂xβ(a(x,η)−E(η))∣≤Cα,βX−μ(∣β∣).(31) |\partial_\eta^\alpha\partial_x^\beta (a(x,\eta)-E(\eta))| \le C_{\alpha,\beta}X^{-\mu(|\beta|)}. \tag{31}

Proof. Write ξ1=E(η)+h\xi_1=E(\eta)+h. On a fixed small collar, P(E(η)+h,η)−λ=h J(h,η),J(h,η)=∫01∂ξ1P(E(η)+th,η) dt.(32) \begin{gathered} P(E(\eta)+h,\eta)-\lambda=h\,J(h,\eta),\\ J(h,\eta)=\int_0^1 \partial_{\xi_1}P(E(\eta)+th,\eta)\,dt. \end{gathered} \tag{32}

The function JJ has a positive lower bound and bounded derivatives. The root equation becomes h=F(x,h,η),F(x,h,η)=−L(x,E(η)+h,η)J(h,η).(33) \begin{gathered} h=F(x,h,\eta),\\ \qquad F(x,h,\eta)=-\frac{L(x,E(\eta)+h,\eta)}{J(h,\eta)}. \end{gathered} \tag{33}

All mixed derivatives of FF, with physical order ∣β∣|\beta|, are bounded by CX−μ(∣β∣)C X^{-\mu(|\beta|)}. Frequency variables are in a fixed compact collar. Since LL vanishes on the ball and μ(0)=δ>0\mu(0)=\delta>0, increasing R0R_0 makes FF map a fixed small closed hh interval into itself and makes ∣Fh∣≤1/2|F_h|\le1/2, uniformly in x,ηx,\eta. Iteration is a contraction and produces the unique small root.

Here is smooth dependence before differentiating that root. Put p=(x,η)p=(x,\eta). On every compact parameter neighborhood, subtracting two fixed-point equations gives ∣h(p+q)−h(p)∣≤2C∣q∣|h(p+q)-h(p)|\le2C|q|. Taylor's segment formula then yields Δh=h(p+q)−h(p),jp=1−Fh(p,h(p)),jpΔh=DpF(p,h(p))q+o(∣q∣). \begin{gathered} \Delta h=h(p+q)-h(p),\\ j_p=1-F_h(p,h(p)),\\ j_p\Delta h=D_pF(p,h(p))q+o(|q|). \end{gathered} The coefficient has absolute value at least 1/21/2, so this proves differentiability and the continuous derivative formula Dph=(1−Fh)−1DpFD_ph=(1-F_h)^{-1}D_pF. Inductively differentiating that formula proves smoothness at every finite order. Uniqueness of the small root makes these local constructions agree.

For the derivative estimates, concavity gives, whenever b0,…,bq≥0b_0,\ldots,b_q\ge0, ∑j=0qμ(bj)≥μ ⁣(∑j=0qbj)+qδ.(34) \sum_{j=0}^q\mu(b_j) \ge \mu\!\left(\sum_{j=0}^q b_j\right)+q\delta. \tag{34}

To see this, fix their sum B>0B>0. A point in the simplex is the convex combination, with weights bj/Bb_j/B, of its vertices. Concavity of the sum bounds its value below by the common vertex value μ(B)+qμ(0)\mu(B)+q\mu(0). For B=0B=0 the inequality is equality.

At order zero, (33) gives (31). Induct in total mixed derivative order. Differentiating the identity produces the term Fh∂ηα∂xβhF_h\partial_\eta^\alpha\partial_x^\beta h, which we move to the left. Every remaining finite chain term is either a derivative of FF alone, or a derivative of FF multiplied by q≥1q\ge1 lower-order derivatives of hh. Their physical orders b0,…,bqb_0,\ldots,b_q sum to ∣β∣|\beta|. The induction hypothesis and (34) bound the term by CX−μ(∣β∣)C X^{-\mu(|\beta|)}, with additional decay when q>0q>0. Division by 1−Fh1-F_h, whose absolute value is at least 1/21/2, proves (31) at the next order. This includes all frequency derivatives. □\square

6. Dividing the polynomial without losing decay

Ellipticity implies that the homogeneous degree-mm part of PP, evaluated at the first coordinate vector, is nonzero. Thus PP has degree mm in ξ1\xi_1, with a nonzero constant leading coefficient. For R0R_0 large the corresponding coefficient of LL is uniformly small. Polynomial division therefore gives

P(ξ)−λ=(ξ1−E(η))Q0(ξ),P(ξ)+L(x,ξ)−λ=(ξ1−a(x,η))Q(x,ξ).(35) \begin{aligned} P(\xi)-\lambda&=(\xi_1-E(\eta))Q_0(\xi),\\ P(\xi)+L(x,\xi)-\lambda&= (\xi_1-a(x,\eta))Q(x,\xi). \end{aligned} \tag{35}

Both quotients have degree m−1m-1 in ξ1\xi_1. On the chosen collar Q0=J(ξ1−E(η),η)Q_0=J(\xi_1-E(\eta),\eta), so it is bounded away from zero.

Proposition 6.1. On every fixed compact collar as above, every mixed derivative satisfies

∣∂ξα∂xβ(Q−Q0)∣≤Cα,βX−μ(∣β∣),∣∂ξα∂xβ(Q−1−Q0−1)∣≤Cα,β′X−μ(∣β∣).(36) \begin{aligned} |\partial_\xi^\alpha\partial_x^\beta(Q-Q_0)| &\le C_{\alpha,\beta}X^{-\mu(|\beta|)},\\ |\partial_\xi^\alpha\partial_x^\beta (Q^{-1}-Q_0^{-1})| &\le C'_{\alpha,\beta}X^{-\mu(|\beta|)}. \end{aligned} \tag{36}

The second assertion holds after increasing R0R_0 if necessary.

Proof. Subtract the first factorization in (35), using the second root:

L+Q0(a−E)=∑j=0mrj(x,η)ξ1j,Q−Q0=∑j=0m−1qj(x,η)ξ1j.(37) \begin{aligned} L+Q_0(a-E)&=\sum_{j=0}^m r_j(x,\eta)\xi_1^j,\\ Q-Q_0&=\sum_{j=0}^{m-1}q_j(x,\eta)\xi_1^j. \end{aligned} \tag{37}

Every derivative of each rjr_j, with physical order ∣β∣|\beta|, has the first bound in (36), by (30), (31) and the product rule. Coefficient comparison gives

qm−1=rm,qj−1=rj+aqj(1≤j<m).(38) \begin{gathered} q_{m-1}=r_m,\\ \qquad q_{j-1}=r_j+a q_j\quad(1\le j<m). \end{gathered} \tag{38}

Induct downward in jj. Multiplication by E(η)E(\eta) costs no physical decay. For multiplication by a−Ea-E, every product-rule term has two physical orders whose sum is ∣β∣|\beta|. Equations (31) and (34) give at least the required exponent, with an extra δ\delta. This proves all coefficient estimates. The polynomial's powers of ξ1\xi_1 and their derivatives are bounded on the fixed collar, giving the first line of (36).

The zeroth estimate and the large zero region give ∣Q−Q0∣≤∣Q0∣/2|Q-Q_0|\le |Q_0|/2. Hence QQ is uniformly bounded away from zero. Set f=Q−1−Q0−1f=Q^{-1}-Q_0^{-1}. Its zeroth bound follows directly, and Qf=(Q0−Q)Q0−1.(39) Qf=(Q_0-Q)Q_0^{-1}. \tag{39}

The right side has every required mixed bound. Induct on total derivative order in (39), isolating Q∂ξα∂xβfQ\partial_\xi^\alpha\partial_x^\beta f. A differentiated Q0Q_0 has only frequency derivatives and a fixed bound. A differentiated Q−Q0Q-Q_0, multiplied by a lower derivative of ff, is controlled by (34). The induction hypothesis thus bounds every other term by CX−μ(∣β∣)C X^{-\mu(|\beta|)}. Divide by the uniformly nonzero QQ. This proves the second line. □\square

These are exact polynomial identities. In left quantization, a product of operators can also contain terms with frequency and physical derivatives of their symbols. Equation (35) alone does not identify the forcing of a first-order equation obtained from a differential equation.

7. Applying the outgoing theorem to the energy root

Choose real ψ∈Cc∞\psi\in C_c^\infty of the η\eta patch, equal to one on a smaller patch, and define b(s,z,η)=ψ(η)a((s,z),η).(40) b(s,z,\eta)=\psi(\eta)a((s,z),\eta). \tag{40}

Extend it by zero outside the patch. Smoothness follows because the cutoff is supported strictly inside the domain of aa. Theorem 5.1 gives uniform bounds for every η\eta derivative. Physical derivatives of the free part ψE\psi E vanish. Since μ(1)=1+δ\mu(1)=1+\delta, ∣∂ηα∇zb∣+∣∂ηα∂sb∣≤Cα(1+∣s∣)−1−δ.(41) \begin{aligned} |\partial_\eta^\alpha\nabla_z b| +|\partial_\eta^\alpha\partial_s b| &\le C_\alpha(1+|s|)^{-1-\delta}. \end{aligned} \tag{41}

The right side is integrable. The kernel proof gives operator-norm continuity in ss from the same bound on its derivative. All hypotheses of Theorem 3.1 hold. Therefore, if v∈B∗v\in B^*, g∈Bg\in B, (Ds−b(s,z,Dz))v=g,(42) (D_s-b(s,z,D_z))v=g, \tag{42}

and vv satisfies (3), then ∥v(s)∥Lz2≤C∫−∞s∥g(t)∥Lz2 dt.(43) \|v(s)\|_{L^2_z} \le C\int_{-\infty}^s\|g(t)\|_{L^2_z}\,dt. \tag{43}

The theorem also constructs its unique outgoing solution for every g∈L1(R;Lz2)g\in L^1(\mathbb R;L^2_z). Once a frequency-localized differential equation has been reduced to (42), this result supplies its slice estimate and uniqueness. Construction of the long-range outgoing amplitude requires further control of the oscillating energy graph.

Use the conclusion

Follow the outgoing kernel estimate to the slice recovery, then check the perturbed energy root and its reciprocal derivatives. The graph factorization and the outgoing uniqueness principle perform different jobs.

8. Exercises with complete solutions

Exercise 1 — Foundation: a wave may leave a permanent outgoing amplitude. Take A(s)=0A(s)=0, g(s,z)=ρ(s)h(z)g(s,z)=\rho(s)h(z), where ρ∈Cc∞((−1,1))\rho\in C_c^\infty((-1,1)), h∈Cc∞(Rd)h\in C_c^\infty(\mathbb R^d), and c=∫ρc=\int\rho. Find the outgoing solution. Prove (3) and compute its total ball-mass limit divided by RR. Explain why adding a nonzero constant-in-ss vector destroys the outgoing condition.

Solution 1. The formula is v(s,z)=ih(z)∫−∞sρ(t) dtv(s,z)=i h(z)\int_{-\infty}^s\rho(t)\,dt. It is zero for s<−1s<-1 and equals ich(z)ic h(z) for s>1s>1. For fixed TT, only a bounded ss interval contributes to (3), with a bounded slice norm; its quotient tends to zero. The strip ∣s∣≤1|s|\le1 likewise contributes O(1/R)O(1/R) to the total quotient. The part s>1s>1 gives ∣c∣2∫∣z∣<R(1−∣z∣2/R2−R−1)+∣h(z)∣2 dz⟶∣c∣2∥h∥22.(44) \begin{aligned} &|c|^2\int_{|z|<R} \left(\sqrt{1-|z|^2/R^2}-R^{-1}\right)_+ \\ &\quad|h(z)|^2\,dz\\ &\longrightarrow |c|^2\|h\|_2^2. \end{aligned} \tag{44}

Dominated convergence proves the limit. If ch≠0c h\ne0, the solution consequently does not have vanishing total shell mass. Adding k(z)≠0k(z)\ne0 gives the homogeneous solution kk on s<−1s<-1. For any fixed TT, its backward ball quotient has limit ∥k∥22\|k\|_2^2: the length of the interval from −R2−∣z∣2-\sqrt{R^2-|z|^2} to min⁡(T,−1)\min(T,-1), divided by RR, tends to one, and is uniformly bounded for large RR. Dominated convergence gives the stated positive limit. In dimension zero the same computations use ∣h∣2|h|^2 and ∣k∣2|k|^2.

Exercise 2 — Intermediate: the sharp derivative budget. Let δ=1/5\delta=1/5 in (29). Find the least possible value of μ(b0)+μ(b1)+μ(b2)\mu(b_0)+\mu(b_1)+\mu(b_2) for nonnegative bjb_j summing to 33. Compare it with the value for (1,1,1)(1,1,1). Explain how this estimate controls a chain term with two differentiated inner roots.

Solution 2. Equation (34) gives the lower bound μ(3)+2δ=14/5+2/5=16/5\mu(3)+2\delta=14/5+2/5=16/5. The vertex (3,0,0)(3,0,0) attains it, so it is the minimum. For (1,1,1)(1,1,1) the value is 3(6/5)=18/53(6/5)=18/5. A chain term with one outer derivative and two inner root derivatives has physical orders b0,b1,b2b_0,b_1,b_2; its decay is at least X−16/5X^{-16/5}. The target for physical order three is X−14/5X^{-14/5}, so even the least favorable distribution has the extra factor X−2/5X^{-2/5}. Zero physical order does not mean zero derivative: a factor with frequency derivatives only still supplies the exponent μ(0)=1/5\mu(0)=1/5.

Exercise 3 — Intermediate: an anisotropic quadratic shell. Take P(ξ1,η)=ξ12+2η2P(\xi_1,\eta)=\xi_1^2+2\eta^2, λ=1\lambda=1, ∣η∣≤1/4|\eta|\le1/4, and L(x,ξ)=ℓ(x)L(x,\xi)=\ell(x), where ℓ\ell is real, satisfies (30), vanishes on a sufficiently large ball and has ∣ℓ∣≤1/8|\ell|\le1/8. Compute E,a,Q0,QE,a,Q_0,Q, and the two differences in (31) and (36). Identify a collar where the denominators are bounded away from zero.

Solution 3. Put E=(1−2η2)1/2E=(1-2\eta^2)^{1/2} and choose the positive perturbed root a=(1−2η2−ℓ)1/2a=(1-2\eta^2-\ell)^{1/2}. Then E≥7/8E\ge\sqrt{7/8} and a≥3/4a\ge\sqrt{3/4}. Exact division gives a−E=−ℓa+E,Q0=ξ1+E,Q=ξ1+a,Q−Q0=a−E,Q−1−Q0−1=−a−E(ξ1+a)(ξ1+E).(45) \begin{gathered} a-E=-\frac{\ell}{a+E},\\ \qquad Q_0=\xi_1+E,\\ \qquad Q=\xi_1+a,\\ Q-Q_0=a-E,\\ \qquad Q^{-1}-Q_0^{-1} =-\frac{a-E}{(\xi_1+a)(\xi_1+E)}. \end{gathered} \tag{45}

On ∣ξ1−E∣<1/4|\xi_1-E|<1/4, ξ1>1/2\xi_1>1/2, so both quotient denominators have a fixed positive lower bound. Theorem 5.1 and Proposition 6.1 apply directly. Alternatively, differentiate the first formula in (45): the square-root arguments stay in a fixed positive interval, and every physical chain term contains derivatives of ℓ\ell. Equation (34) bounds their products by the required X−μ(∣β∣)X^{-\mu(|\beta|)}. Frequency derivatives have fixed bounded factors on the patch. The remaining formulas give the same estimates by the product and reciprocal rules.

Exercise 4 — Advanced: why reality is not operator symmetry. In one transverse dimension take b(s,z,η)=q(s)a0(z)ψ(η)b(s,z,\eta)=q(s)a_0(z)\psi(\eta), with real smooth compactly supported a0,ψa_0,\psi, and q(s)=(1+s2)−(1+δ)/2q(s)=(1+s^2)^{-(1+\delta)/2}. Show that the adjoint defect has integrable norm. Give choices for which it is nonzero. Derive the C/RC/R bound on its commutator with χR\chi_R.

Solution 4. If k(w)=(2π)−1∫eiwηψ(η) dηk(w)=(2\pi)^{-1}\int e^{iw\eta}\psi(\eta)\,d\eta, the kernel is q(s)a0(z)k(z−y)q(s)a_0(z)k(z-y). Reality of ψ\psi gives k(−w)‾=k(w)\overline{k(-w)}=k(w), so the adjoint-difference kernel is q(s)(a0(z)−a0(y))k(z−y).(46) q(s)\bigl(a_0(z)-a_0(y)\bigr)k(z-y). \tag{46}

Since a0a_0 is Lipschitz and kk is rapidly decreasing, (8) bounds the norm by C∣q(s)∣C|q(s)|. Its integral is finite because q(s)=O(∣s∣−1−δ)q(s)=O(|s|^{-1-\delta}). For a concrete nonzero defect, choose a nonnegative, even, nonzero ψ\psi; then k(0)>0k(0)>0 and kk is nonzero near zero. Choose a0a_0 with a0′(z0)≠0a_0'(z_0)\ne0. Nearby distinct z,yz,y give a nonzero continuous kernel in (46), so the operator is nonzero: tests supported in sufficiently small neighborhoods of these points give a nonzero pairing. Finally the commutator kernel is q(s)a0(z)k(z−y)(χR(s,y)−χR(s,z))q(s)a_0(z)k(z-y)(\chi_R(s,y)-\chi_R(s,z)). The difference is bounded by C∣z−y∣/RC|z-y|/R. Both absolute kernel integrals are bounded by C/RC/R, uniformly in ss, and (8) gives that operator norm.

Exercise 5 — Advanced: the larger forcing class is necessary. Choose h∈Cc∞(Rd)h\in C_c^\infty(\mathbb R^d), supported in ∣z∣<1|z|<1, with ∥h∥2=1\|h\|_2=1, and set ρ(s)=∑j≥12−jj1[3⋅4j, 3⋅4j+1](s),g(s,z)=ρ(s)h(z).(47) \begin{aligned} \rho(s)&=\sum_{j\ge1}\frac{2^{-j}}j 1_{[3\cdot4^j,\ 3\cdot4^j+1]}(s),\\ g(s,z)&=\rho(s)h(z). \end{aligned} \tag{47}

Show that g∈L1(R;Lz2)g\in L^1(\mathbb R;L^2_z) but g∉Bg\notin B. For A=0A=0, verify the outgoing solution's existence, bounded slice norm and condition (3). In dimension zero use h=1h=1.

Solution 5. The slice integral is ∑j≥12−j/j<∞\sum_{j\ge1}2^{-j}/j<\infty. Each pulse is wholly inside the shell with outer radius 22j+2=4⋅4j2^{2j+2}=4\cdot4^j and inner radius 2⋅4j2\cdot4^j; the fixed transverse support and interval of length one preserve both strict inequalities. Distinct pulses occupy distinct shells. That shell's L2L^2 norm is 2−j/j2^{-j}/j, so its contribution to the BB norm is 2j+1(2−j/j)=2/j2^{j+1}(2^{-j}/j)=2/j. The sum diverges.

Nevertheless v(s,z)=ih(z)∫−∞sρ(t) dtv(s,z)=i h(z)\int_{-\infty}^s\rho(t)\,dt is continuous, locally absolutely continuous, and solves Dsv=gD_sv=g. Its slice norm is at most ∑j2−j/j\sum_j2^{-j}/j, so (4) puts it in B∗B^*. For each fixed TT, its nonzero slices with s<Ts<T occupy a bounded interval starting at s=12s=12. Their total squared mass is finite, independently of RR, and hence (3) holds. The uniqueness theorem applies to this forcing even though its BB norm is infinite.

References