Modified waves and the direction of escape

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

Working question: Which spatial cone corresponds to the chosen time direction? A packet of velocity vv follows x≈tvx\approx tv. When t<0t<0, that ray points toward −v-v, so decay imposed only near the positive-velocity cone can miss the packet entirely. The modifier must cancel the phase along the rays selected by its sign before its residual is integrated.

A long-range force changes the phase of a travelling wave even after its coefficient becomes small. A modifier follows that phase. The remaining error can then be integrated in time. We prove this for a real polynomial free operator, a rough short-range differential perturbation, and any self-adjoint realization of their sum. The direction in which a packet escapes determines which spatial cone must have short-range decay.

Read Smooth long-range phases from Hamilton trajectories for the real phase and all its frequency derivatives. The arbitrary-self-adjoint measure and equality with the original second-moment domain are proved in Self-adjoint spectral calculus with the original domain. Its unitary-group proof and both directions of the generator-domain criterion will be used below. The maximal polynomial Fourier realization and its spectral projections are proved in Resolvents, domains and spectral density, Section 1. These statements include operators with no lower bound, on arbitrary Hilbert spaces. Section 2 below proves the uniform nonstationary integral estimate. The complete finite-derivative L2L^2 estimate and its exact dilation convention are proved in A finite-derivative bound for left quantization. Its measure and Euclidean product proofs supply dominated convergence, Fatou and Tonelli. Approximation and convolution, Theorems 2.1 and 3.1, supplies the smooth approximation used below. The required finite seminorms are estimated in Section 4. The elementary real powers, logarithms and trigonometry supply the scalar formulas used in the examples.

The free texts of Oh [O] and Teschl [T] provide background on evolution and characteristics. The comparison argument below integrates the evolution error after inserting the stated phase modifier.

1. A packet has a velocity and a time direction

Use the unitary Fourier transform and Dj=−i∂xjD_j=-i\partial_{x_j}. Let P0P_0 be a real nonconstant polynomial on Rn\mathbb R^n. Its maximal Fourier multiplication operator is

H0=P0(D),D(H0)={u:P0u^∈L2},v(ξ)=∇P0(ξ),Ω={ξ:v(ξ)≠0}.(1) \begin{gathered} H_0=P_0(D),\qquad \mathcal D(H_0)=\{u:P_0\widehat u\in L^2\},\\ v(\xi)=\nabla P_0(\xi),\qquad \Omega=\{\xi:v(\xi)\ne0\}. \end{gathered} \tag{1}

A free packet with frequencies near ξ0\xi^0 travels near x=tv(ξ0)x=tv(\xi^0). Positive time samples the direction v(ξ0)v(\xi^0); negative time samples its opposite. A cone here is open and invariant under positive dilations. It need not contain the negative of any of its vectors.

Write a finite-order differential perturbation as

V(x,D)=∑∣α∣≤pvα(x)Dα=VL(x,D)+VS(x,D),∫Rn∑α∣vα(x)∣2X−2M dx<∞,X=(1+∣x∣2)1/2.(2) \begin{gathered} V(x,D)=\sum_{|\alpha|\le p}v_\alpha(x)D^\alpha\\ =V_L(x,D)+V_S(x,D),\\ \int_{\mathbb R^n}\sum_\alpha |v_\alpha(x)|^2X^{-2M}\,dx<\infty,\\ X=(1+|x|^2)^{1/2}. \end{gathered} \tag{2}

The coefficients in (2) are measurable and can be complex. The operator H0+VH_0+V is assumed symmetric on Schwartz space and to have a self-adjoint extension HH. The weighted condition makes V:S→L2V:\mathcal S\to L^2 continuous: after increasing MM to a nonnegative integer if necessary, ∥vαDαu∥2≤∥vαX−M∥2sup⁡x∣XMDαu(x)∣\|v_\alpha D^\alpha u\|_2\le\|v_\alpha X^{-M}\|_2\sup_x|X^M D^\alpha u(x)|, and the finite sum is bounded by Schwartz seminorms. It does not assert that VV is bounded on L2L^2 or that D(H)=D(H0)\mathcal D(H)=\mathcal D(H_0).

The long coefficients are real and smooth. Fix 0<δ<1/30<\delta<1/3, and assume

∣∂βvα,L(x)∣≤CαβX−m(∣β∣),m(j)=j+δ(0≤j≤2),m(j)=1+(1+δ)j/2(j≥2).(3) \begin{gathered} |\partial^\beta v_{\alpha,L}(x)| \le C_{\alpha\beta}X^{-m(|\beta|)},\\ m(j)=j+\delta\quad(0\le j\le2),\\ m(j)=1+(1+\delta)j/2\quad(j\ge2). \end{gathered} \tag{3}

The phase theorem in the prerequisite gives a real smooth W(ξ,t)W(\xi,t). On every compact K⊂ΩK\subset\Omega, for both sufficiently late time directions, it satisfies

y(ξ,t)=∂ξW(ξ,t),∂tW=P0(ξ)+VL(y(ξ,t),ξ),∣∂ξα(W−tP0)∣≤Cα,K∣t∣1+∣α∣−m(∣α∣).(4) \begin{gathered} y(\xi,t)=\partial_\xi W(\xi,t),\\ \partial_tW=P_0(\xi)+V_L(y(\xi,t),\xi),\\ |\partial_\xi^\alpha(W-tP_0)| \le C_{\alpha,K}|t|^{1+|\alpha|-m(|\alpha|)} . \end{gathered} \tag{4}

Every frequency derivative is included. The threshold depends on KK. In particular y/t=v+O(∣t∣−δ)y/t=v+O(|t|^{-\delta}). The modifier e−iW(D,t)e^{-iW(D,t)} is unitary for every tt.

For a cone Γ\Gamma, measure the short coefficients by their unit-ball L2L^2 envelope:

AΓ(r)=sup⁡r<∣z∣<2r(∫∣y−z∣<1y∈Γ∑α∣vα,S(y)∣2 dy)1/2,∫1∞AΓ(r) dr<∞.(5) \begin{aligned} A_\Gamma(r)&=\sup_{r<|z|<2r}\\ &\quad\left(\int_{\substack{|y-z|<1\\ y\in\Gamma}} \sum_\alpha|v_{\alpha,S}(y)|^2\,dy\right)^{1/2},\\ &\quad\int_1^\infty A_\Gamma(r)\,dr<\infty. \end{aligned} \tag{5}

This permits narrow or irregular peaks. Only their local square integrals in the selected cone are controlled. The envelope is a measurable function of r>0r>0: if AΓ(r0)>aA_\Gamma(r_0)>a, one center zz realizes a value greater than aa, and its strict inequalities r<∣z∣<2rr<|z|<2r persist for rr near r0r_0. Thus {r:AΓ(r)>a}\{r:A_\Gamma(r)>a\} is open. On a compact interval of positive radii every relevant unit ball lies in one bounded ball, so local square integrability bounds the envelope there.

Fix σ∈{+1,−1}\sigma\in\{+1,-1\}. Our short-range hypothesis for this time direction is: for almost every ξ∈Ω\xi\in\Omega, there is a cone Γ\Gamma containing σv(ξ)\sigma v(\xi) for which (5) holds. Sections 2–4 prove the error estimates. Section 5 uses them to construct the limit at t→σ∞t\to\sigma\infty. Section 6 explains why the sign in this hypothesis is necessary.

2. Concentration despite large high derivatives of the phase

Let u^=a∈Cc∞(Ω)\widehat u=a\in C_c^\infty(\Omega), and let ω\omega be a neighborhood of v(supp⁡a)v(\operatorname{supp}a). Put ut=e−iW(D,t)uu_t=e^{-iW(D,t)}u.

Lemma 2.1 (concentration). For every multiindex α\alpha and every NN, for sufficiently large ∣t∣|t|,

∣∂xαut(x)∣≤Cα,N(∣x∣+∣t∣)−N,x/t∉ω.(6) |\partial_x^\alpha u_t(x)| \le C_{\alpha,N}(|x|+|t|)^{-N}, \qquad x/t\notin\omega . \tag{6}

Proof. Set

τ=∣t∣,c=(1−δ)/2,r=1−c,q=τc,L=∣x∣+τ,R=W−tP0.(7) \begin{gathered} \tau=|t|,\qquad c=(1-\delta)/2,\qquad r=1-c,\\ q=\tau^c,\qquad L=|x|+\tau,\qquad R=W-tP_0. \end{gathered} \tag{7}

On the compact frequency support, Rξ=O(τ1−δ)R_\xi=O(\tau^{1-\delta}) and every derivative of RR of order j≥2j\ge2 is O(τcj)O(\tau^{cj}). The first estimate says that the actual velocity is close to the free velocity. The higher estimates can grow with jj; a fixed unscaled phase estimate would not have uniform constants.

Outside the stated velocity neighborhood,

∣x−tv(ξ)∣≥γL,∣x−y(ξ,t)∣≥(γ/2)L(8) |x-tv(\xi)|\ge\gamma L,\qquad |x-y(\xi,t)|\ge(\gamma/2)L \tag{8}

on a slightly enlarged compact frequency set, once τ\tau is large. For bounded x/tx/t, use the positive distance from the compact velocity set to the closed complement of ω\omega. For large ∣x/t∣|x/t|, use the triangle inequality. The error O(τ1−δ)O(\tau^{1-\delta}) is smaller than γL/2\gamma L/2.

Change variables to η=qξ\eta=q\xi. Here is a smooth lattice partition with uniform bounds. Using the explicit smooth bump construction, choose b≥0b\ge0, supported in (−1,1)n(-1,1)^n and strictly positive on [−1/2,1/2]n[-1/2,1/2]^n. Set S(η)=∑g∈Znb(η−g),χ0(η)=b(η)/S(η). S(\eta)=\sum_{g\in\mathbb Z^n}b(\eta-g), \qquad \chi_0(\eta)=b(\eta)/S(\eta). Only a uniformly bounded number of summands occur near any point. Hence SS is smooth and periodic, with all derivatives bounded. Every point is within the closed half-cube of some lattice point; the positive minimum of bb on that cube bounds SS away from zero. Thus χ0\chi_0 has bounded derivatives, fixed compact support, and periodicity of SS gives ∑gχ0(η−g)=1\sum_g\chi_0(\eta-g)=1. At most CqnCq^n terms meet qsupp⁡aq\operatorname{supp}a. After translation η=g+z\eta=g+z, every cell has a fixed support and uniformly bounded amplitude derivatives. Its oscillatory parameter and normalized phase are

λ=L/q,ϕg(z)=x⋅(g+z)/q−W((g+z)/q,t)λ.(9) \begin{gathered} \lambda=L/q,\\ \phi_g(z)= \frac{x\cdot(g+z)/q-W((g+z)/q,t)}{\lambda}. \end{gathered} \tag{9}

Subtract ϕg(0)\phi_g(0); this changes only a scalar of modulus one. By (8), the normalized gradient has norm at least γ/2\gamma/2. All its derivatives are bounded uniformly. For order j≥2j\ge2, the free part contributes at most Cτq−j/λ≤Cq1−jC\tau q^{-j}/\lambda\le Cq^{1-j}. The correction contributes at most Cτcjq−j/λ=C/λC\tau^{cj}q^{-j}/\lambda=C/\lambda. The first-order upper bound follows directly from ∣y∣≤Cτ|y|\le C\tau. After the constant is subtracted, the zeroth-order bound follows on the fixed cell from the first-order bound. Enlarging the frequency support slightly supplies the fixed neighborhood needed for these estimates.

Here is the nonstationary bound in the exact uniform form needed. On a cell set Bg=∇ϕg/∣∇ϕg∣2B_g=\nabla\phi_g/|\nabla\phi_g|^2. The lower gradient bound and the finite derivative bounds just established bound every fixed finite derivative of BgB_g, uniformly in the cell and in x,tx,t. For the compactly supported cell amplitude aga_g, integration by parts gives

∫eiλϕgag=−1iλ∫eiλϕgdiv⁡(Bgag). \int e^{i\lambda\phi_g}a_g =-\frac1{i\lambda}\int e^{i\lambda\phi_g} \operatorname{div}(B_g a_g).

Iterate this identity AA times. Each resulting amplitude is a finite sum of products of derivatives of BgB_g and derivatives of aga_g of orders at most AA; their supports lie in the same fixed compact cell. Their L1L^1 norms are uniformly bounded. Hence the cell integral is bounded by CAλ−AC_A\lambda^{-A}, with a constant independent of the lattice cell and of x,tx,t.

The Jacobian q−nq^{-n} cancels the number O(qn)O(q^n) of cells. Since q≤Lcq\le L^c, we have λ≥L1−c\lambda\ge L^{1-c}. Increase AA to obtain any prescribed power in (6). Spatial differentiation replaces aa by a fixed polynomial times aa, so the argument applies unchanged. The variable x/tx/t and the estimates in τ\tau include negative time. ∎

We will also need a time-dependent amplitude. The same proof works when it has a fixed compact frequency support and

∣∂ξαat∣≤Cαq∣α∣.(10) |\partial_\xi^\alpha a_t|\le C_\alpha q^{|\alpha|}. \tag{10}

Its rescaled cell derivatives are bounded, which is all the cancellation proof used.

3. Rough coefficients are integrable along the selected rays

Suppose the packet support KK has one cone Γ\Gamma as in (5), containing σv(K)\sigma v(K). Shrink a bounded velocity neighborhood ω\omega so its closure avoids zero and σω‾⊂Γ\sigma\overline\omega\subset\Gamma. For t=στt=\sigma\tau, the set tωt\omega lies inside Γ\Gamma, between radii aτa\tau and bτb\tau, with 0<a<b0<a<b.

Lemma 3.1. On this time half-line,

∫T∞∥VS(x,D)uστ∥2 dτ<∞.(11) \int_T^\infty \|V_S(x,D)u_{\sigma\tau}\|_2\,d\tau<\infty . \tag{11}

Proof. A unit ball meeting tωt\omega has its center in a slightly wider annulus, with radii still comparable to τ\tau. Cover that annulus by finitely many overlapping intervals (djτ,2djτ)(d_j\tau,2d_j\tau). Their constants dj>0d_j>0 are fixed. Every coefficient square integral on such a ball, restricted to tωt\omega, is at most B(τ)2B(\tau)^2, where

B(τ)=∑j=1JAΓ(djτ),∫T∞B(τ) dτ<∞.(12) B(\tau)=\sum_{j=1}^J A_\Gamma(d_j\tau),\qquad \int_T^\infty B(\tau)\,d\tau<\infty . \tag{12}

The last assertion follows by substituting r=djτr=d_j\tau. Overlapping intervals avoid any problem at strict annulus endpoints.

For an integer s>n/2s>n/2, a translated fixed cutoff, Fourier inversion, Cauchy–Schwarz and Plancherel give the local bound

sup⁡B(z,1)∣∂αw∣2≤C∑∣β∣≤s∫B(z,2)∣∂α+βw(y)∣2 dy.(13) \begin{aligned} \sup_{B(z,1)}|\partial^\alpha w|^2 &\le C\sum_{|\beta|\le s} \int_{B(z,2)} |\partial^{\alpha+\beta}w(y)|^2\,dy . \end{aligned} \tag{13}

Indeed, apply the Fourier bound for the supremum to a cutoff times ∂αw\partial^\alpha w. The weight (1+∣ξ∣2)−s(1+|\xi|^2)^{-s} is integrable. Integer-order Plancherel and Leibniz bound its HsH^s norm by the displayed local derivatives. All constants are independent of the center zz.

Insert 1=∣B(0,1)∣−1∫1{∣y−z∣<1} dz1=|B(0,1)|^{-1}\int 1_{\{|y-z|<1\}}\,dz in the square norm on tωt\omega. Cauchy–Schwarz in the finite coefficient sum, (12), (13) and Tonelli yield

∥1tωVSut∥22≤CB(τ)2∑∣ν∣≤p+s∥∂νut∥22≤CuB(τ)2.(14) \begin{aligned} \|1_{t\omega}V_Su_t\|_2^2 &\le C B(\tau)^2 \sum_{|\nu|\le p+s}\|\partial^\nu u_t\|_2^2\\ &\le C_u B(\tau)^2 . \end{aligned} \tag{14}

The norms in the middle line are constant in time, because differentiation commutes with the unitary Fourier modifier.

The short coefficients also satisfy a weighted square-integrability condition as in (2): subtract the decaying long coefficients and, if needed, increase MM. On the complement of tωt\omega, (6) therefore bounds the square norm by CNsup⁡xX2M(∣x∣+τ)−2NC_N\sup_x X^{2M}(|x|+\tau)^{-2N}. Choose N≥M+2N\ge M+2. This is at most Cτ−4C\tau^{-4}. Consequently

∥VSut∥2≤Cu(B(τ)+τ−2).(15) \|V_Su_t\|_2 \le C_u\bigl(B(\tau)+\tau^{-2}\bigr). \tag{15}

Both terms are integrable. This proves (11) for the chosen sign. ∎

No derivatives of the rough coefficients were taken. The proof uses their local square integrals, rather than replacing them by a pointwise bound.

4. The long-range difference gains an integrable power

We next compare a smooth coefficient at the physical position xx with its value at the phase position y=Wξy=W_\xi. Select a small velocity ball centered at v0≠0v^0\ne0, of radius less than ∣v0∣/2|v^0|/2. Choose a compact frequency cutoff χ0\chi_0 and a velocity cutoff χ\chi, so χ=1\chi=1 near all velocities on the packet support. On their supports, for large τ\tau, both x/tx/t and y/ty/t lie in this ball. The joining segment then stays at distance comparable to τ\tau from zero, for either time sign.

Let ff be one long coefficient, satisfying (3). Set

fj(x,y)=∫01∂jf(y+s(x−y)) ds,f(x)−f(y)=∑j(xj−yj)fj(x,y),Gj,t(x,ξ)=χ(x/t)χ0(ξ)fj(x,y(ξ,t)).(16) \begin{gathered} f_j(x,y)= \int_0^1\partial_jf(y+s(x-y))\,ds,\\ f(x)-f(y)=\sum_j(x_j-y_j)f_j(x,y),\\ G_{j,t}(x,\xi)= \chi(x/t)\chi_0(\xi)f_j(x,y(\xi,t)). \end{gathered} \tag{16}

Joint derivatives of fjf_j of total order ll have size at most Clτ−m(l+1)C_l\tau^{-m(l+1)} in this region.

Lemma 4.1 (all symbol derivatives). For every α,β\alpha,\beta,

∣∂ξα∂xβGj,t∣≤Cαβτ∣α∣−m(∣α∣+∣β∣+1).(17) \begin{aligned} |\partial_\xi^\alpha\partial_x^\beta G_{j,t}| &\le C_{\alpha\beta} \tau^{|\alpha|-m(|\alpha|+|\beta|+1)} . \end{aligned} \tag{17}

Proof. The positive-order derivatives of yy have bounds Ckτa(k)C_k\tau^{a(k)}, where

a(1)=1,a(k)=c(k+1)(k≥2).(18) a(1)=1,\qquad a(k)=c(k+1)\quad(k\ge2). \tag{18}

This follows from (4), including the free part of size O(τ)O(\tau). Put k=∣α∣k=|\alpha|, b=∣β∣b=|\beta|. For k=0k=0, the joint derivative bound on fjf_j already gives (17). For k>0k>0, a chain-rule term has q≥1q\ge1 outer yy derivatives and positive inner orders kik_i, with ∑iki=k\sum_i k_i=k. Let ℓ\ell count the inner orders at least two. The others equal one, so

∑ia(ki)=q−ℓ+c(k−q+2ℓ)=ck+rq−δℓ.(19) \begin{aligned} \sum_i a(k_i) &=q-\ell+c(k-q+2\ell)\\ &=ck+rq-\delta\ell . \end{aligned} \tag{19}

Here r=(1+δ)/2r=(1+\delta)/2. Since q+b+1≥2q+b+1\ge2, the linear part of mm gives the exponent

−m(q+b+1)+∑ia(ki)=ck−1−r(b+1)−δℓ≤k−m(k+b+1).(20) \begin{aligned} &-m(q+b+1)+\sum_i a(k_i)\\ &=ck-1-r(b+1)-\delta\ell\\ &\le k-m(k+b+1). \end{aligned} \tag{20}

This controls every finite chain partition. Frequency derivatives falling on χ0\chi_0 are harmless because k−m(k+b+1)k-m(k+b+1) is nondecreasing in kk. Each position derivative falling on χ(x/t)\chi(x/t) contributes τ−1\tau^{-1}; the increments of mm are at most one, so the same target at the total position order follows. This proves (17). ∎

We need only a fixed finite number of these bounds to control the operators. Both Gj,tG_{j,t} and its first frequency derivative obey the weaker convenient estimate

∣∂ξα∂xβat∣≤Cαβτ−1−δτc∣α∣−r∣β∣,at=Gj,t or ∂ξjGj,t.(21) \begin{gathered} |\partial_\xi^\alpha\partial_x^\beta a_t|\\ \le C_{\alpha\beta}\tau^{-1-\delta} \tau^{c|\alpha|-r|\beta|},\\ a_t=G_{j,t}\ \text{or }\partial_{\xi_j}G_{j,t}. \end{gathered} \tag{21}

For GG, the zeroth estimate is exact and the positive-order exponent in (17) is smaller by cc. For ∂ξjG\partial_{\xi_j}G, all relevant arguments of mm are at least two, and (21) is exactly its exponent.

Let at(x,D)a_t(x,D) mean left, or Kohn–Nirenberg, quantization. The unitary dilation Uτw(z)=τnc/2w(τcz)U_\tau w(z)=\tau^{nc/2}w(\tau^cz) conjugates its symbol to at(τcz,τ−cη)a_t(\tau^cz,\tau^{-c}\eta). Divide this symbol by τ−1−δ\tau^{-1-\delta}. Each derivative is bounded by a constant times τ(c−r)∣β∣=τ−δ∣β∣≤1\tau^{(c-r)|\beta|}=\tau^{-\delta|\beta|}\le1. Choose one integer N>n/2N>n/2. The finite-derivative L2L^2 theorem requires only ∣α∣≤2N|\alpha|\le2N, ∣β∣≤2N|\beta|\le2N for each conjugated symbol. For both symbols in (21), (17) therefore uses at most 2N+12N+1 frequency and 2N2N position derivatives of GG.

The segment formula uses at most 4N+24N+2 derivatives of ff and 2N+22N+2 derivatives of WW, all supplied by (3)–(4). Thus the theorem applies with constants independent of τ\tau, and

∥Gj,t(x,D)∥2→2+∥∂ξjGj,t(x,D)∥2→2≤Cτ−1−δ.(22) \begin{gathered} \|G_{j,t}(x,D)\|_{2\to2}\\ +\|\partial_{\xi_j}G_{j,t}(x,D)\|_{2\to2}\\ \le C\tau^{-1-\delta}. \end{gathered} \tag{22}

This reduction needs no additional symbol-class theorem at a moving scale.

The order of the factors matters. Left quantization gives

∑jGj,t(x,D)(xj−yj(D,t))=χ(x/t)(f(x)−f(y(D,t)))χ0(D)−i∑j∂ξjGj,t(x,D).(23) \begin{aligned} &\sum_jG_{j,t}(x,D)(x_j-y_j(D,t))\\ &=\chi(x/t)\bigl(f(x)-f(y(D,t))\bigr)\\ &\quad\chi_0(D) -i\sum_j\partial_{\xi_j}G_{j,t}(x,D). \end{aligned} \tag{23}

To verify the sign, use G(x,D)xj=xjG(x,D)−i(∂ξjG)(x,D)G(x,D)x_j=x_jG(x,D)-i(\partial_{\xi_j}G)(x,D), then the scalar identity in (16). On a packet,

(xj−yj(D,t))ut=e−iW(D,t)(xju).(24) (x_j-y_j(D,t))u_t=e^{-iW(D,t)}(x_j u). \tag{24}

Fourier differentiation proves this identity directly. Choose χ0=1\chi_0=1 on a neighborhood of the Fourier support of uu, which also contains the supports of its frequency derivatives. Equations (22)–(24) imply

∥χ(x/t)(f(x)−f(y(D,t)))ut∥2≤Cτ−1−δ(∑j∥xju∥2+∥u∥2).(25) \begin{aligned} &\|\chi(x/t)(f(x)-f(y(D,t)))u_t\|_2\\ &\le C\tau^{-1-\delta}\\ &\quad\left(\sum_j\|x_j u\|_2+\|u\|_2\right). \end{aligned} \tag{25}

Both terms outside χ\chi are rapidly small. The f(x)utf(x)u_t term follows from (6) and boundedness of ff. For the multiplier term, the finite chain rule gives

∣∂ξαf(y(ξ,t))∣≤Cατc∣α∣−δ(26) |\partial_\xi^\alpha f(y(\xi,t))| \le C_\alpha\tau^{c|\alpha|-\delta} \tag{26}

on the fixed frequency support. For order zero this is the coefficient decay. At positive order kk, if the outer order is one, the exponent is −δ-\delta for k=1k=1, and −1−δ+c(k+1)≤ck−δ-1-\delta+c(k+1)\le ck-\delta for k≥2k\ge2. If the outer order is q≥2q\ge2, (19) and m(q)=1+rqm(q)=1+rq give ck−1−δℓ≤ck−δck-1-\delta\ell\le ck-\delta. This proves (26), so the time-dependent-amplitude version (10) of concentration applies to f(y(ξ,t))u^f(y(\xi,t))\widehat u.

Apply (25) to every packet DαuD^\alpha u in the finite differential sum. The rapid tails have arbitrarily integrable powers. We have proved

∫T∞∥(VL(x,D)−VL(y(D,στ),D))uστ∥2 dτ<∞.(27) \begin{gathered} \int_T^\infty \|(V_L(x,D)\\ -V_L(y(D,\sigma\tau),D)) u_{\sigma\tau}\|_2\,d\tau<\infty . \end{gathered} \tag{27}

To make the rapid-tail step explicit, a pointwise bound CN(∣x∣+τ)−NC_N(|x|+\tau)^{-N} has L2L^2 norm at most CNτn/2−NC_N\tau^{n/2-N}, by substituting x=τzx=\tau z in its square integral when 2N>n2N>n. Choose N>n/2+1N>n/2+1 to obtain a time-integrable bound. The amplitude estimate (26) and the fixed compact Fourier supports allow this choice for every differential term. This part is valid for both time directions; the short-range cone condition is where the sign enters.

5. Integrate the error, then extend by density

Theorem 5.1 (modified wave operators). Assume (2)–(4) and the self-adjoint extension described in Section 1. Fix a sign σ\sigma. Suppose (5) holds in a cone containing σv(ξ)\sigma v(\xi) for almost every regular frequency ξ\xi. Then the strong limit

Wσu=lim⁡t→σ∞eitHe−iW(D,t)u,u∈L2,(28) \mathcal W_\sigma u =\lim_{t\to\sigma\infty} e^{itH}e^{-iW(D,t)}u,\qquad u\in L^2, \tag{28}

exists and is an isometry. It satisfies, for every real ss,

eisHWσ=WσeisH0,WσD(H0)⊂D(H),HWσu=WσH0u.(29) \begin{gathered} e^{isH}\mathcal W_\sigma =\mathcal W_\sigma e^{isH_0},\\ \mathcal W_\sigma\mathcal D(H_0)\subset\mathcal D(H), \qquad H\mathcal W_\sigma u=\mathcal W_\sigma H_0u. \end{gathered} \tag{29}

If the hypothesis holds for both signs, both operators exist. The theorem imposes no ellipticity, no lower bound, and no equality of the full free and perturbed domains.

Proof. Let Ωσ\Omega_\sigma be the set of regular frequencies whose signed velocity belongs to at least one cone satisfying (5). This set is open: membership of a continuous velocity in a fixed open cone persists nearby. Its complement inside Ω\Omega is null by hypothesis. The critical set is null as well. Indeed, one partial derivative of the nonconstant polynomial is a nonzero polynomial. In one variable a nonzero polynomial has finitely many roots: factoring out one root reduces the degree, so induction bounds their number by the degree. In dimension nn, write a nonzero polynomial as ∑jaj(ξ′)ξnj\sum_j a_j(\xi')\xi_n^j, with at least one nonzero coefficient polynomial aja_j. Outside that coefficient's null zero set, supplied by induction in dimension, every one-dimensional fiber has finitely many roots. Tonelli on each bounded box makes the full zero set null; the exceptional base set contributes zero because the box has bounded height. A countable union of boxes covers the space. Applying this fact to the nonzero partial derivative proves the critical-set assertion and hence that Ωσ\Omega_\sigma has full measure.

For completeness, put Kj={ξ:∣ξ∣≤j,dist⁡(ξ,Ωσc)≥1/j}. K_j=\{\xi:|\xi|\le j, \operatorname{dist}(\xi,\Omega_\sigma^c)\ge1/j\}. When the complement is empty take its distance to be infinity. These compact subsets of Ωσ\Omega_\sigma increase to that open set. For any a∈L2a\in L^2, dominated convergence gives 1Kja→a1_{K_j}a\to a. For fixed jj, convolution with a smooth unit-integral bump supported in a ball of radius ε<1/(2j)\varepsilon<1/(2j) gives a compact smooth function supported inside Ωσ\Omega_\sigma. The programme convolution proof makes its L2L^2 error tend to zero with ε\varepsilon. Plancherel therefore proves that the selected packets are dense in L2L^2.

Cover the compact support of such a packet by finitely many frequency neighborhoods, each mapped by the signed velocity into one good cone and by the velocity into one ball separated from zero. The finite bump partition proof splits the packet into finitely many compact smooth pieces in those neighborhoods. Lemma 3.1 and (27) give an integrable residual on the chosen time end for each piece.

For bounded time intervals, ut=e−iW(D,t)uu_t=e^{-iW(D,t)}u is a C1C^1 Schwartz path. Its Fourier transform has fixed compact support and a smooth amplitude. The weighted condition makes VutVu_t continuous in L2L^2. The extension HH agrees with H0+VH_0+V on this initial Schwartz domain, so ut∈D(H)u_t\in\mathcal D(H) and HutHu_t is continuous. The unitary-group product rule yields

ddt(eitHut)=ieitH(H0+V−Wt(D,t))ut=ieitH(VS+VL(x,D)−VL(y(D,t),D))ut(30) \begin{aligned} \frac{d}{dt}(e^{itH}u_t)&\\ &=i e^{itH}(H_0+V-W_t(D,t))u_t\\ &=i e^{itH} \bigl(V_S+V_L(x,D)\\ &\quad -V_L(y(D,t),D)\bigr)u_t \end{aligned} \tag{30}

on the sufficiently late half-line. The last equality uses the exact equation in (4). The product rule follows by splitting a difference quotient into the change in the vector and the group derivative on a fixed domain vector. Graph continuity verifies the domain requirements.

The integral of the norm of (30) is finite. On a finite interval its continuous vector integral exists by the Riemann-sum construction in Approximation and convolution, Section 5. Pair with an arbitrary fixed vector and apply the scalar fundamental theorem; commuting the pairing with the Riemann limit proves that this integral equals the difference of the endpoint vectors. The integral norm is bounded by the integral of the norm. Its vanishing tail therefore makes eitHute^{itH}u_t Cauchy at the selected end. Every approximating operator is unitary, so the limit preserves the packet's norm. For a general uu, choose a dense packet hh. If At=eitHe−iW(D,t)A_t=e^{itH}e^{-iW(D,t)}, then ∥At(u−h)−Ar(u−h)∥≤2∥u−h∥\|A_t(u-h)-A_r(u-h)\|\le2\|u-h\|. This proves convergence and isometry on all of L2L^2.

To prove intertwining, fix s∈Rs\in\mathbb R. On every compact K⊂ΩK\subset\Omega,

W(ξ,t)−W(ξ,t−s)−sP0(ξ)=∫t−stVL(y(ξ,r),ξ) dr=OK,s(∣t∣−δ).(31) \begin{aligned} &W(\xi,t)-W(\xi,t-s)-sP_0(\xi)\\ &=\int_{t-s}^{t}V_L(y(\xi,r),\xi)\,dr\\ &=O_{K,s}(|t|^{-\delta}). \end{aligned} \tag{31}

The oriented integral includes negative ss. Its exponential tends uniformly to one on KK. The multiplier has modulus one, and Ω\Omega has full measure, so dominated convergence proves its strong convergence to the identity on all of L2L^2. Reindexing the limit by t+st+s gives

eisHWσe−isH0u=lim⁡t→σ∞eitHe−iW(D,t)ei(W(D,t)−W(D,t−s)−sH0)u=Wσu.(32) \begin{gathered} e^{isH}\mathcal W_\sigma e^{-isH_0}u\\ =\lim_{t\to\sigma\infty} e^{itH}e^{-iW(D,t)}\\ \quad e^{i(W(D,t)-W(D,t-s)-sH_0)}u\\ =\mathcal W_\sigma u . \end{gathered} \tag{32}

For u∈D(H0)u\in\mathcal D(H_0), the group difference quotient in (29) therefore converges on Wσu\mathcal W_\sigma u to iWσH0ui\mathcal W_\sigma H_0u. The general self-adjoint spectral domain condition identifies this limit with iHWσuiH\mathcal W_\sigma u. More explicitly, boundedness of the quotients and Fatou in their spectral square integrals imply the finite second moment; dominated convergence then identifies the derivative. This proves the domain inclusion and operator identity. ∎

Spectral measures of the range. The range is closed: a convergent sequence of image vectors comes from a Cauchy sequence, since the map is an isometry. Every vector in the range has an absolutely continuous spectral measure, as follows.

First consider the free measure of uu, which is the pushforward of ∣u^(ξ)∣2dξ|\widehat u(\xi)|^2d\xi by P0P_0. If N⊂RN\subset\mathbb R is a null Borel set, then P0−1(N)P_0^{-1}(N) is null. The critical set was proved null above. Near each remaining point, some partial derivative is nonzero and the coordinate inverse theorem gives the coordinates (P0(ξ),ξ′)(P_0(\xi),\xi'). Tonelli makes N×Rn−1N\times\mathbb R^{n-1} null, and the proved change-of-variables formula gives null measure to its preimage in that chart. A countable chart cover suffices: choose from the countable rational-ball base balls lying in these chart neighborhoods and covering the regular set. Thus every free scalar spectral measure is absolutely continuous.

Next, let JJ be any isometry satisfying eisHJ=JeisH0e^{isH}J=Je^{isH_0}, in particular J=WσJ=\mathcal W_\sigma. For g=Jug=Ju, isometry and the group identity give (e−itHg,g)=(e−itH0u,u)(e^{-itH}g,g)=(e^{-itH_0}u,u). For any self-adjoint AA and Im⁡z>0\operatorname{Im}z>0, scalar spectral integration yields

(RA(z)f,f)=i∫0∞eitz(e−itAf,f) dt. (R_A(z)f,f) =i\int_0^\infty e^{itz}(e^{-itA}f,f)\,dt.

Indeed, the absolute double integral is bounded by ∥f∥2/Im⁡z\|f\|^2/\operatorname{Im}z. Fubini therefore applies, and the inner scalar integral is i∫0∞eit(z−λ)dt=(λ−z)−1i\int_0^\infty e^{it(z-\lambda)}dt=(\lambda-z)^{-1}. The two scalar resolvent functions for gg and uu are consequently equal. They determine equal finite spectral measures. To check this last step including atoms, multiply the imaginary part at a+iεa+i\varepsilon by ε\varepsilon; its kernel ε2/((λ−a)2+ε2)\varepsilon^2/((\lambda-a)^2+\varepsilon^2) is bounded by one and tends to the indicator of {a}\{a\}. Dominated convergence gives equality of atom masses. The proved scalar Stone formula then gives equality on every bounded open interval, including its half-endpoint terms. Equality of total masses follows from isometry. The finite-measure uniqueness lemma, applied to intervals and the whole line, gives equality on every Borel set. Hence the measure of JuJu is absolutely continuous.

Equality of the range with the whole absolutely continuous subspace is a further assertion. Here is the isometry calculation used in Exercise 7.5. If J+J_+ and J−J_- are isometries with the same range K\mathcal K, then J±∗J±=IJ_\pm^*J_\pm=I and J±J±∗=PKJ_\pm J_\pm^*=P_{\mathcal K}. The latter identity follows directly: it fixes each vector J±fJ_\pm f and annihilates its orthogonal complement by the adjoint pairing. Therefore, for S=J+∗J−S=J_+^*J_-, both S∗SS^*S and SS∗SS^* are the identity. Without equality of ranges, the adjoint bound gives only ∥S∥≤1\|S\|\le1. These adjoints and orthogonal projections have their elementary proof in Hilbert-space tools.

6. A cone for positive velocity can miss every negative-time ray

On the line, choose a real smooth 0≤χ−≤10\le\chi_-\le1, equal to one for x≤−1x\le-1 and zero for x≥0x\ge0. For b>0b>0, put

H0=D3,H=D3+bχ−(x),VL=0,W(ξ,t)=tξ3.(33) \begin{gathered} H_0=D^3,\qquad H=D^3+b\chi_-(x),\\ V_L=0,\qquad W(\xi,t)=t\xi^3 . \end{gathered} \tag{33}

The operator HH is self-adjoint on H3(R)H^3(\mathbb R). Fourier multiplication gives self-adjointness of D3D^3, and adding the bounded real multiplier preserves self-adjointness on that domain. To verify this directly, for self-adjoint AA and bounded symmetric BB, the adjoint equation for A+BA+B says that the adjoint equation for AA holds with the vector w−Bzw-Bz. Thus D((A+B)∗)=D(A∗)\mathcal D((A+B)^*)=\mathcal D(A^*), and (A+B)∗=A+B(A+B)^*=A+B.

The same elementary domain calculation also justifies unitary changes of variables used in Exercise 7.4. For a unitary UU, testing adjoint pairings shows (UAU∗)∗=UA∗U∗(UAU^*)^*=UA^*U^* with domain UD(A∗)UD(A^*). Its spectral projections are UEA(B)U∗UE_A(B)U^*, by bounded truncation of the moment formula and uniqueness of the spectral measure. In particular its group is UeitAU∗Ue^{itA}U^*. If UU multiplies by a smooth scalar of modulus one with bounded first derivative, the weak product rule gives U:H1→H1U:H^1\to H^1; the same holds for U∗U^*. One can verify this rule first on compact smooth functions and pass to the limit using the proved H1H^1 density and the two bounded multiplier norms. Thus the differential conjugation in that exercise holds on the whole stated domain.

The weighted coefficient condition holds with M=1M=1. Every nonzero frequency has velocity 3ξ2>03\xi^2>0. In the cone (0,∞)(0,\infty), the short coefficient vanishes, so (5) is exactly zero. The positive wave operator exists by Theorem 5.1.

The negative limit with this modifier fails. Take nonzero uu with u^∈Cc∞((a,d))\widehat u\in C_c^\infty((a,d)), 0<a<d0<a<d. For negative time and x≥−1x\ge-1, the free phase derivative is x+3∣t∣ξ2x+3|t|\xi^2, bounded below by a positive multiple of ∣x∣+∣t∣|x|+|t| for sufficiently large ∣t∣|t|. All normalized phase derivatives are bounded. The nonstationary estimate and integration in the half-line give

∥1{x≥−1}e−itH0u∥2≤CN∣t∣1/2−N.(34) \|1_{\{x\ge-1\}}e^{-itH_0}u\|_2 \le C_N|t|^{1/2-N}. \tag{34}

Compare first with the correct constant background on the left: Bt=eitHe−it(H0+b)B_t=e^{itH}e^{-it(H_0+b)}. Its derivative on this packet is ieitHb(χ−−1)e−it(H0+b)ui e^{itH}b(\chi_--1)e^{-it(H_0+b)}u. The coefficient vanishes for x≤−1x\le-1. Equation (34) with N=2N=2 is integrable, so Cook's argument gives Btu→ZuB_tu\to Zu as t→−∞t\to-\infty. Since each BtB_t is unitary, ∥Zu∥=∥u∥>0\|Zu\|=\|u\|>0.

But the proposed comparison satisfies

eitHe−itH0u=eibtBtu,tn=−2πn/b,sn=−(2n+1)π/b.(35) \begin{gathered} e^{itH}e^{-itH_0}u=e^{ibt}B_tu,\\ t_n=-2\pi n/b,\qquad s_n=-(2n+1)\pi/b. \end{gathered} \tag{35}

Its limits along these two sequences are ZuZu and −Zu-Zu. They differ. The missing negative cone permits a nonzero background phase on the left; the positive cone contains no information about it.

Example (33) shows that a cone condition at the positive velocity alone does not control both time directions. Theorem 5.1 uses the signed velocity required by the general nonelliptic assertion.

When the short coefficients satisfy (5) in every direction, both signs are available. The example addresses the directional hypothesis when coefficients are allowed to behave differently in different cones.

Use the conclusion

Apply the signed-velocity counterexample before reading the general existence proof. Keep the frequency derivatives of the phase, rough short-range coefficient bounds and arbitrary self-adjoint realization as distinct inputs.

7. Five exercises with complete solutions

Exercise 7.1 (first check: frequency is not escape direction). On the line let P0(ξ)=ξ3+ξP_0(\xi)=\xi^3+\xi. Find the regular set and the velocities of packets near ξ=1\xi=1 and ξ=−1\xi=-1. For each time sign, identify the relevant half-line. Does changing the sign of frequency exchange those half-lines?

Solution 7.1. The derivative is 3ξ2+13\xi^2+1, which is positive everywhere. Thus Ω=R\Omega=\mathbb R, and both selected packets have velocity 44. Their positive-time centers are near 4t>04t>0; their negative-time centers are near 4t<04t<0. The good cone for the positive limit is the positive half-line, and for the negative limit it is the negative half-line. The two frequency signs have the same velocity. Frequency sign alone does not determine the direction sampled by this polynomial.

Exercise 7.2 (a decay condition without a spare power). In Rn\mathbb R^n, let fβ(x)=1(e+X) [log⁡(e+X)]β,β≥0.(36) f_\beta(x)= \frac{1}{(e+X)\,[\log(e+X)]^\beta}, \qquad \beta\ge0 . \tag{36} Use this real multiplication potential as VSV_S, with VL=0V_L=0. For the cone Γ=Rn∖{0}\Gamma=\mathbb R^n\setminus\{0\}, determine exactly when (5) holds. Explain why the result gives wave existence for any real nonconstant P0P_0, although fβf_\beta is not bounded by a fixed CX−1−εCX^{-1-\varepsilon} when β>1\beta>1.

Solution 7.2. For large rr, a unit ball with center r<∣z∣<2rr<|z|<2r has XX comparable to rr, and log⁡(e+X)\log(e+X) comparable to log⁡r\log r. Its volume is fixed; removal of the single point zero changes no integral. Both upper and lower bounds therefore give

AΓ(r)≍1r(log⁡r)β.(37) A_\Gamma(r)\asymp \frac1{r(\log r)^\beta}. \tag{37}

The substitution s=log⁡rs=\log r makes its tail integral ∫s−β ds\int s^{-\beta}\,ds, which is finite exactly for β>1\beta>1. The bounded initial radial interval is harmless. The potential is bounded and real, so H0+fβH_0+f_\beta is self-adjoint on D(H0)\mathcal D(H_0), by the bounded-perturbation argument in Section 6. The weighted condition holds after choosing MM large enough. The cone contains both signed velocities, so Theorem 5.1 supplies both limits when β>1\beta>1. For every ε>0\varepsilon>0, the ratio of fβf_\beta to X−1−εX^{-1-\varepsilon} is comparable to Xε/(log⁡X)βX^\varepsilon/(\log X)^\beta, which tends to infinity. The integral criterion still applies. For β≤1\beta\le1, this calculation says that the stated sufficient hypothesis fails; it makes no claim about whether a wave operator can exist by another argument.

Exercise 7.3 (the derivative count and the operator sign). Set δ=1/5\delta=1/5. Compute the exponent in (17) for ∣α∣=3|\alpha|=3, ∣β∣=2|\beta|=2, and for the same derivatives of ∂ξjGj,t\partial_{\xi_j}G_{j,t}. Compare with (21). Check the sign of the last term in (23), and bound the integral from TT to infinity of an operator error bounded by Cτ−1−δC\tau^{-1-\delta}.

Solution 7.3. Here c=2/5c=2/5, r=3/5r=3/5, and m(j)=1+3j/5m(j)=1+3j/5 for j≥2j\ge2. The two exact exponents are

3−m(6)=3−23/5=−8/5,4−m(7)=4−26/5=−6/5.(38) \begin{aligned} 3-m(6)&=3-23/5=-8/5,\\ 4-m(7)&=4-26/5=-6/5. \end{aligned} \tag{38}

The convenient bound (21) has exponent −6/5+(2/5)3−(3/5)2=−6/5-6/5+(2/5)3-(3/5)2=-6/5. It is weaker than the first exact bound by 2/5=c2/5=c, and equals the second. After the dilation and division by τ−6/5\tau^{-6/5}, these derivatives have exponent at most (c−r)2=−2/5(c-r)2=-2/5, so they are bounded for τ≥1\tau\ge1.

Fourier differentiation gives [xj,G(x,D)]=i(∂ξjG)(x,D)[x_j,G(x,D)]=i(\partial_{\xi_j}G)(x,D). Therefore G(x,D)xj=xjG(x,D)−i(∂ξjG)(x,D)G(x,D)x_j=x_jG(x,D)-i(\partial_{\xi_j}G)(x,D), proving the minus sign in (23). Finally

∫T∞Cτ−6/5 dτ=5CT−1/5.(39) \int_T^\infty C\tau^{-6/5}\,d\tau =5C T^{-1/5}. \tag{39}

The frequency-derivative commutator is integrable with exactly the same operator bound as the main coefficient difference.

Exercise 7.4 (solve a complete long-range comparison). Let a>0a>0, b∈Rb\in\mathbb R, 0<δ<1/30<\delta<1/3, and H0=aDH_0=aD, H=aD+b(1+x2)−δ/2H=aD+b(1+x^2)^{-\delta/2}. Define F(0)=0F(0)=0, F′(x)=b(1+x2)−δ/2/aF'(x)=b(1+x^2)^{-\delta/2}/a. Find a phase satisfying (4), calculate both modified wave operators, and determine whether their ranges are complete.

Solution 7.4. The real multiplication potential is bounded. Its jj-th derivative is O(X−j−δ)O(X^{-j-\delta}), which implies (3): the exponents agree through order two, and for j≥2j\ge2 the difference j+δ−m(j)=c(j−2)j+\delta-m(j)=c(j-2) is nonnegative. All velocities equal aa. Put

W(ξ,t)=atξ+F(at).(40) W(\xi,t)=at\xi+F(at). \tag{40}

Then Wξ=atW_\xi=at and Wt=aξ+b(1+a2t2)−δ/2W_t=a\xi+b(1+a^2t^2)^{-\delta/2}, giving the exact equation. The correction is O(∣t∣1−δ)O(|t|^{1-\delta}); all its positive frequency derivatives vanish. Thus every bound in (4) holds.

The unitary multiplier Uu=e−iF(x)uUu=e^{-iF(x)}u preserves H1H^1, because F′F' is bounded. Differentiation gives HU=UH0HU=UH_0. Hence H=UH0U∗H=UH_0U^* on H1H^1, and its evolution is UeitH0U∗Ue^{itH_0}U^*. Since eitH0u(x)=u(x+at)e^{itH_0}u(x)=u(x+at), the comparison operator is

eitHe−iW(D,t)u(x)=e−iF(x)ei(F(x+at)−F(at))u(x).(41) \begin{gathered} e^{itH}e^{-iW(D,t)}u(x)\\ =e^{-iF(x)} e^{i(F(x+at)-F(at))}u(x). \end{gathered} \tag{41}

For each fixed xx, F(x+at)−F(at)F(x+at)-F(at) is the integral of F′F' over an interval of fixed length. It tends to zero at both time ends because F′F' tends to zero. The exponential has modulus one, so dominated convergence in the square norm gives W+=W−=U\mathcal W_+=\mathcal W_-=U. These operators are unitary onto L2L^2. The transport operator aDaD has entirely absolutely continuous spectrum by Fourier multiplication and the linear change of variable λ=aξ\lambda=a\xi. Its unitary conjugate HH has the same spectral type. Thus both ranges equal Hac(H)=L2\mathcal H_{\mathrm{ac}}(H)=L^2.

Exercise 7.5 (the freedom in a fixed phase). Assume both limits in (28) exist. At each time end, replace W(ξ,t)W(\xi,t) by W(ξ,t)+Cσ(ξ)W(\xi,t)+C_\sigma(\xi), where CσC_\sigma is a real measurable finite function. Find the new wave operators, their domain intertwining, and the change in S=W+∗W−S=\mathcal W_+^*\mathcal W_-. This replacement need not solve the same Hamilton–Jacobi equation.

Solution 7.5. Write Mσ=e−iCσ(D)M_\sigma=e^{-iC_\sigma(D)}, a unitary Fourier multiplier. The new approximant is exactly AtMσA_tM_\sigma, so

W~σ=WσMσ,S~=eiC+(D)Se−iC−(D).(42) \begin{gathered} \widetilde{\mathcal W}_\sigma=\mathcal W_\sigma M_\sigma,\\ \widetilde S=e^{iC_+(D)}S e^{-iC_-(D)} . \end{gathered} \tag{42}

Multiplication by a function of modulus one preserves the Fourier condition P0u^∈L2P_0\widehat u\in L^2, so MσM_\sigma preserves D(H0)\mathcal D(H_0) and commutes there with H0H_0. It also commutes with the free unitary group. Equation (29) therefore proves both the group and domain identities for the new wave operators. Their ranges and isometry properties are unchanged. If the original two ranges coincide, SS is unitary by the isometry calculation in Wave operators and modified phases, and the displayed formula keeps S~\widetilde S unitary. If their ranges have not been shown to coincide, both SS and S~\widetilde S are only known from these assumptions to be contractions. A fixed phase changes the comparison's coordinates; it cannot supply a missing completeness proof.

References

[O] Sung-Jin Oh, Lecture Notes for Math 222A, free evolving lecture notes, University of California, Berkeley, Fall 2023, §2.4.1, pp. 23–24.

[T] Gerald Teschl, Mathematical Methods in Quantum Mechanics: With Applications to Schrödinger Operators, free author's online second edition, 2014, Theorem 5.1, pp. 145–146, Theorem 12.2 and Lemma 12.3, pp. 284–285.

[HW] Lars Hörmander, The existence of wave operators in scattering theory, freely readable journal scan, 1976, Theorem 3.9 and its proof, pp. 83–86, and Lemma A.1, p. 88.