Propagation of Regularity for Schrödinger Equations with Time Dependent Potentials

Abstract

The dynamics of Schrödinger equation with time dependent potentials of general time dependence is considered. It is shown that for localized in space potentials, there is propagation of regularity which is uniformly bounded in higher Sobolev norms. Unlike the cases where the solution scatter, and then propagation is proved via a standard bootstrap argument, the solutions considered here have a part that does not scatter, as expected in general. For this we introduce propagation estimates that work directly in (e.g.) $H^{2}(\mathbb{R}^{3}).$

We prove Propagation of Regularity for Schrödinger Equations with Time Dependent Potentials, namely, the $H^{2}$ norm of the solution remains uniformly bounded in time. We also prove some general Propagation Estimates that may be of interest.

1 Introduction

The Schrödinger equation with time dependent potential is fundamental; it corresponds to open quantum systems, as well as nonlinear equations. It appears naturally in linearization of non-linear dynamics, in long range N-body scattering, see e.g. Sigal and Soffer (1993) and charge transfer hamiltonians Yajima (1980). See also Chen and Mezher (2026).

Yet, it is under-developed at the theory level, and for a good reason: one cannot use the same type of spectral theory (following linear algebra) as in the time independent case. Hence, the well known methods and results are based on special cases which are mostly perturbative, like small perturbations, adiabatic time dependence, time periodic or vanishing as time goes to infinity.

On the mathematical side, we have some well-known important results on the theory of global existence; see, e.g. Yosida (2012); Reed and Simon (1975); Kitada and Yajima (1982, 1983); Ruiz and Vega (1994); Linares and Ponce (2014). In particular, it is shown that for localized potentials that we consider below, the solutions exist in $H^{2}.$ However, these estimates do not bound the size, since the $H^{2}$ can grow exponentially in time.

In many cases, the relevant dispersive estimates for the solution can be proved, such as the Strichartz estimate or $L^{p}$ decay estimates. This can then be used to prove the uniform boundedness of higher norms by applying Gronwall’s inequalities to derivatives of the equation. This can be done for defocusing Nonlinear wave equations.

For problems with time dependent potentials, it can be done for special types of time dependenceRodnianski and Schlag (2004); Beceanu and Soffer (2019); Soffer and Wu (2025); Rodnianski et al. (2005); Fishman and Soffer (2014); Teufel (2003). Although there are results on scattering theory for time dependent potentials, they do not have the needed dispersive estimates, Kitada and Yajima (1983); Enss and Veselić (1983). For small potentials, the N-body case was treated by Yajima (2025)

In recent years scattering theory was developed to include both time dependent potentials and nonlinear equations with focusing part. Tao (2008); Liu and Soffer (2020, 2025); Soffer and Wu (2022); soffer2025ℒ. They rely on the assumption that the solution is uniformly bounded in $H^{1}.$ They do not prove decay estimates which are useful for propagation of regularity. However, it should be noted that in the works of Tao, and later in the works of Liu-Soffer and Soffer-Wu, the weakly localized part of the solution is shown to be smooth in some cases.

In contrast with these methods, in this work a direct approach is used to control the higher Sobolev norm, avoiding the proofs of decay estimates. Instead, we prove Propagation Estimates(PRES) directly in higher norms. Such estimates are typically proven by the use of Propagation Observables (PROB). PROB is an operator $B(t)$ (or a family of operators depending on time) which is self adjoint and such that its expectation on the state at time $t$, $\langle\psi(t),B(t)\psi(t)\rangle$ has a non-negative time derivative. There are many such known operators which give useful apriori estimates on solutions of both linear and non-linear equations. The first class is composed of operators which correspond to symmetries of the equation, and lead to conservation laws. The constant function implies $L^{2}$ conservation. Functions of the hamiltonian $H$ in the time independent case (energy conservation) etc… A second class of operators are derived from the dispersive nature of the solutions of defocusing equations, and include the Morawetz estimate and the Dilation identity.

Modern estimates are now based on microlocalized objects, where only the leading order term of the derivative is non-negative, and there are higher order corrections, which need to be controlled in $L^{1}(dt).$ We will follow tools from Sigal and Soffer (1986); Graf (1990); Derezinski and Gérard (2013); Hunziker et al. (1999); Sigal and Soffer (1987, 1988); Liu and Soffer (2025); Soffer and Wu (2022); Soffer (2011). See also Soffer (2024) for a review.

So the typical estimate looks like

$$\langle\psi(T),B(T)\psi(T)\rangle-\langle\psi(t),B(t)\psi(t)\rangle=\int_{t}^{T}\langle\psi(s),C^{*}C(s)\psi(s)\rangle ds+$$

$$+\int_{t}^{T}\langle\psi(s),R(s)\psi(s)\rangle ds.\quad\langle R(s)\rangle\in L^{1}(ds).$$

Here we introduce and use the notation $\langle\psi(t),Q(t)\psi(t)\rangle\equiv\langle Q(t)\rangle.$

The proof of the propagation estimates is done by proving on various parts of the phase-space the propagation of regularity using the above argument and generalizations thereof, including summing over time slices, and summing over dyadic (in spectral space) PROBs.

2 Smoothness of the Solution on Compact Domains in Space

We consider the Schrödinger Equation (SE)

$$i\frac{\partial\psi}{\partial t}=-\Delta\psi+V(x,t)\psi\equiv(H_{0}+V)\psi,$$

$$\psi(0)\in H^{s}(\mathbb{R}^{n}),\quad s\geqslant 2,\quad n\geq 3.$$

$V(x,t)$ is assumed to be localized in $x$ and smooth:

$$\left\|\langle x\rangle^{\sigma}V(x,t)\right\|_{L^{\infty}}<c,\quad\sigma>4.$$

$$\left\|\langle x\rangle^{\sigma}\left[\left|\frac{\partial V}{\partial t}\right|+\left|\nabla_{x}V\right|+|\Delta_{x}V|\right]\right\|_{L^{\infty}}<c,\quad\sigma>4.$$

Here, we use the standard definition

$$\langle x\rangle\equiv 1+|x|^{2}.$$

Remark 2.1.

These conditions are not optimal, and apply for $s\leq 2$. In particular, the restrictive $L^{\infty}$ conditions are used to streamline the presentation.

The first step of the proof is to show that the solution is uniformly bounded in $\mathrm{H}^{2}$ (say), locally in space.

Theorem 2.2 (Local Smoothness).

Let $\psi(t)$ be a solution of the (SE) with initial data in $\mathrm{H}^{2}(\mathbb{R}^{n}),n\geq 3.$

Then,

$$\sup_{t}\left\|\langle x\rangle^{-2-\delta}\Delta\psi(t)\right\|_{L^{2}}\lesssim\|\psi(0)\|_{H^{2}}.$$

Proof.

We use the Duhamel representation.

$$\langle x\rangle^{-\sigma}\Delta\psi(t)=\langle x\rangle^{-\sigma}e^{-iH_{0}t}\Delta\psi(0)-i\int_{0}^{t}\langle x\rangle^{-\sigma}\Delta e^{-iH_{0}(t-s)}V\left(x,s\right)\psi(s)ds$$

The first term is bounded in $H^{2}$, since $\psi(0)\in H^{2}$ by assumption. The Duhamel term, we only need to control it in $L^{2}$ for $p\left(\equiv-i\nabla_{x}\right)$ large. Split the integral of time into two parts $\int_{t-\varepsilon}^{t}$ and $\int_{0}^{t-\varepsilon}$. The part $[t-\varepsilon,t]$ is controlled by Bootstrap. The rest is controlled by the smoothing property of the free flow.

$$\langle x\rangle^{-\sigma}\int_{0}^{t-\varepsilon}\Delta e^{i\Delta(t-s)}V(x,s)\psi(s)ds$$ (1)

$$=\int_{0}^{t-\varepsilon}\langle x\rangle^{-\sigma}\Delta\langle x-2pt\rangle^{-2}e^{i\Delta(t-s)}\langle x\rangle^{2}V(x,s)\psi(s)ds$$

$$\tau\equiv t-s$$

Next, we need

Proposition 2.3 (Local Decay and Smoothing).

For $\sigma=2,$

$$\|\left<x\right>^{-\sigma}\langle x-2p\tau\rangle^{-2}\|\leq\mathcal{O}(\left<x\right>^{2}+4p^{2}\tau^{2})^{-1}.$$ (2)

Proof.

The proof follows an algebraic calculation:

$$\begin{gathered}\langle x\rangle^{-2}(-\Delta)\langle x-2p\tau\rangle^{-2}=\langle x\rangle^{-2}\left(-\Delta+x^{2}/4\tau^{2}-x^{2}/4\tau^{2}\right)\langle x-2p\tau\rangle^{-2}\\ =\langle x\rangle^{-2}\left(-\Delta+x^{2}/4\tau^{2}+1\right)\left\langle-\Delta 4\tau^{2}+x^{2}+1\right\rangle^{-1}\left[\left(1+x^{2}-4\Delta\tau^{2}\right)\langle x-2p\tau)^{-2}\right]\\ +\left(O\left(\langle x\rangle^{-2}\right)+O\left(\tau^{-2}\right)\right)\langle x-2p\tau\rangle^{-2}\end{gathered}$$

$$\begin{gathered}\langle x\rangle^{-2}\left(1+x^{2}-4\Delta\tau^{2}\right)\langle x-2p\tau\rangle^{-2}=\langle x\rangle^{-2}\left(1+x^{2}-4\tau^{2}\Delta\right)\left[\frac{-1}{1+x^{2}-4\Delta\tau^{2}}+\frac{1}{\langle x-2p\tau\rangle^{2}}\right]+\\ +\langle x\rangle^{-2}=\langle x\rangle^{-2}\\ +\langle x\rangle^{-2}\left(1+x^{2}-4\tau^{2}\Delta\right)\frac{1}{1+x^{2}-4\Delta\tau^{2}}2A\tau\frac{1}{\langle x-2p\tau\rangle^{2}}\\ =\langle x\rangle^{-2}+\langle x\rangle^{-2}2A\tau\frac{1}{\langle x-2p\tau\rangle^{2}}.\\ A\equiv\frac{1}{2}x\cdot p+\frac{1}{2}p\cdot x=x\cdot p+c\end{gathered}$$

To control the $A$ term, we use (for $\tau<1$)

$$2A\tau =2\tau(x\cdot p+c)=$$

$$=4\tau^{2}\left(c/\tau+\left(\frac{x}{2\tau}\frac{1}{\delta}\right)\cdot(\delta p)\right)$$

$$\leqslant c\tau+\tau^{2}\left(p^{2}\delta^{2}+\frac{1}{\delta^{2}}\frac{x^{2}}{4\tau^{2}}\right)$$

$$=c\tau+4\tau^{2}(-\Delta)\delta^{2}+\frac{1}{\delta^{2}}x^{2}$$

$$=c\tau+4\tau^{2}(-\Delta)\delta^{2}+\delta^{2}x^{2}+\left(\frac{1}{\delta^{2}}-1\right)x^{2}$$

Therefore,

$$\langle x\rangle^{-2}2A\tau\langle x-2p\tau\rangle^{-2}=$$

$$=O(\tau\rangle\langle x\rangle^{-2}+\delta^{2}\langle x\rangle^{-2}\left(-4\tau^{2}\Delta+x^{2}\right)\langle x-2p\tau\rangle^{-2}$$

$$+\left(\frac{1}{\delta^{2}}-1\right)O(1)\langle x-2p\tau\rangle^{-2}$$

By choosing $\delta$ sufficiently small, the $\delta^{2}$ term is absorbed by the LHS. For $\tau>1$, we use the same estimate on

$$\left(\frac{x}{2\tau}\cdot\frac{1}{\delta}\right)\cdot(\delta p)+(\delta p)\cdot\left(\frac{x}{2\tau}\frac{1}{\delta}\right)\text{, which }$$

removes the term $O(\tau)\langle x\rangle^{-2}$. ∎

With the above proposition, we proceed to control the first time interval of (1).

We need to show that the integral over $[0,+t-\varepsilon]$ is uniformly bounded in $L^{2}$. In the control of the $A$ term, we lose a factor of $\langle x\rangle^{-2}$. Therefore, to get the integral bound, we need to use an extra factor of $\langle x\rangle^{-1-\delta}$ : $\left\|\langle x\rangle^{-1-\delta}e^{i\Delta(t-s)}F_{2}(p)\langle x\rangle^{-1-\delta}\right\|_{L^{2}\rightarrow L^{2}}\leq c\langle t-s\rangle^{-1-\delta}.$

The above estimates show that for $|p|\geqslant 1$ the localized solution is uniformly bounded in $H^{2}$. It remains to control the integral part of $[0,t-\varepsilon]$ with low frequency and control the Duhamel term on $[t-\varepsilon,t]$.

Low Frequency

For low frequency, we use the standard estimate $L^{\infty}$ (in three or more dimensions) and $\sigma>n/2:$

$$\left\|\langle x\rangle^{-\sigma}e^{i\Delta(t-s)}f\right\|_{L^{2}}\leqslant\left\|\langle x\rangle^{-\sigma}\right\|_{L^{2}}\left\|e^{i\Delta(t-s)}f\right\|_{L^{\infty}}$$

$$\leqslant\left\|\langle x\rangle^{-\sigma}\right\|_{L^{2}}|t-s|^{-n/2}\|f\|_{L^{1}}.$$

For $\sigma>\frac{n}{2}$ the first term is bounded. For $n\geqslant 3$, the integral over $s$ is bounded.

In our case $f\equiv V(x,s)\psi(s)\in L^{1}$ for $V(x,s)\in L_{s}^{\infty}L_{x}^{2}$.

The interval $[t-\varepsilon,t]$.

We want to control the expression

$$\langle x\rangle^{-\sigma}\int_{t-\varepsilon}^{t}e^{i\Delta(t-s)}\Delta V(x,s)\psi(s)ds.$$

Integration by parts gives:

$$\begin{gathered}\int_{t-\varepsilon}^{t}\Delta e^{i\Delta(t-s)}V(s)\psi(s)ds=\int_{t-\varepsilon}^{t}i\partial_{s}\left(e^{i\Delta(t-s)}V\psi\right)ds\\ -i\int_{t-\varepsilon}^{t}e^{i\Delta(t-s)}\left(\frac{\partial V}{\partial s}\psi(s)+V(-i)i\frac{\partial\psi}{\partial s}\right)ds\\ =iV(t)\psi(t)-ie^{i\Delta\varepsilon}V(t-\varepsilon)\psi(t-\varepsilon)\\ -i\int_{t-\varepsilon}^{t}e^{i\Delta(t-s)}\frac{\partial V}{\partial s}\psi(s)ds\\ -i\int_{t-\varepsilon}^{t}e^{i\Delta(t-s)}V(-i\Delta-iV)\psi(s)ds.\end{gathered}$$

The first two terms are uniformly bounded in $L_{x}^{2}$. The third term times $\langle x\rangle^{-\sigma}$ gives

$$\int_{t-\varepsilon}^{t}\left\|\langle x\rangle^{-\sigma}e^{i\Delta(t-s)}\frac{\partial V}{\partial s}\psi(s)\right\|_{L^{2}}ds$$

$$\leqslant\int_{t-\varepsilon}^{t}\|\langle x\rangle^{-\sigma}\|_{L^{\infty}}|\frac{\partial V}{\partial s}\psi(s)\|_{L_{s}^{\infty}L_{x}^{2}}ds,$$

$$\leq\int_{t-\varepsilon}^{t}\left\|\frac{\partial V}{\partial s}\psi(s)\right\|_{L_{s}^{\infty}L_{x}^{2}}ds<c\epsilon\left\|\frac{\partial V}{\partial s}\right\|_{\infty}\|\psi\|_{2}.$$

The estimate of the last term is similar, where we have $V^{2}\psi(s)$ replacing $\frac{\partial V}{\partial s}$, except that the Laplacian term is bootstrapped:

$$\left\|\langle x\rangle^{-\sigma}\int_{t-\varepsilon}^{t}e^{iH(t-s)}V\Delta\psi(s)ds\right\|_{L^{2}_{x}}$$

$$\quad\leqslant c\varepsilon\|\langle x\rangle^{\sigma}V\|_{L^{\infty}}\left\|\langle x\rangle^{-\sigma}\Delta\psi(s)\right\|_{L_{s}^{\infty}(t-\varepsilon,t)L_{x}^{2}},$$

Hence we need $\left\|V\langle x\rangle^{\sigma}\right\|_{L^{\infty}_{x}}<c$. ∎

3 Control for Large x

In this section, we upgrade control of derivatives for large $x$.

Proposition 3.1.

Under the previous assumptions and Local Smoothness Theorem, we have that

$$\|p\psi(t)\|_{L^{2}}\leqslant t^{a},\quad t\gg 1$$

$1/5=a\quad$ for $|V|\approx\langle x\rangle^{-3-\varepsilon}$ and $|\nabla V|+\left|\partial_{t}V\right|\leqslant\langle x\rangle^{-6-\varepsilon}$.

Proof.

Let $F_{2}(\lambda\geq 1)$ be a smooth characteristic function of the interval $[1,\infty).$ We use the following PROB:

$$\begin{gathered}h(t)=F_{2}\left(\frac{|p|}{t^{\alpha}}\geqslant 1\right)H(t)+H(t)F_{2}\left(\frac{|p|}{t^{\alpha}}\geqslant 1\right)\\ H(t)=p^{2}+V(x,t).\end{gathered}$$

Then,

$$\partial_{t}h(t)= \left\langle i\left[V,F_{2}\right]H+Hi\left[V,F_{2}\right]\right\rangle$$

$$+\left\langle F_{2}\frac{\partial V}{\partial t}+\frac{\partial V}{\partial t}F_{2}\right\rangle$$

$$+\left\langle-\frac{\alpha}{t}\tilde{F}_{2}^{\prime}H+H\left(-\frac{\alpha}{t}\right)\tilde{F}_{2}^{\prime}\right\rangle\equiv I_{1}+I_{2}+I_{3}$$

$I_{2}$ has a localization as a result of the decay in $x$ of $\frac{\partial V}{\partial t}$.

Hence

$$I_{2}= \left\langle F_{2}(\Delta+1)^{-1}(-\Delta+1){\frac{\partial V}{\partial t}}\right\rangle+c\cdot c$$

$$t^{-2\alpha}\left\langle\Delta\psi,\tilde{F}_{2}\frac{\partial V}{\partial t}\psi\right\rangle+c\cdot c\cdot=$$

$$= t^{-2\alpha}\left\langle\Delta\psi,\langle x\rangle^{-\sigma}\langle x\rangle^{\sigma}\tilde{F}_{2}\frac{\partial V}{\partial t}\psi\right\rangle+c\cdot c.$$

$$\leq ct^{-2\alpha}\left\|\langle x\rangle^{-\sigma}\Delta\psi\right\|_{2}\left\|\langle x\rangle^{\sigma}\frac{\partial V}{\partial t}\right\|_{L_{t,x}^{\infty}}\|\psi\|_{2}.$$

$$I_{3}= -\frac{\alpha\cdot 2}{t}\left\langle\tilde{F}_{2}^{\prime}p^{2}\right\rangle+\left\langle-\frac{\alpha}{t}\left(\tilde{F}_{2}^{\prime}V+V\tilde{F}_{2}^{\prime}\right)\right\rangle$$

$$\leq O\left(t^{-1+2\alpha}\right)+O(t^{-1})$$

We used that $\langle x\rangle^{\sigma}\tilde{F}_{2}\frac{\partial V}{\partial t}$ is higher order in $x,t$, since by commuting powers of $x$ through $\tilde{F}_{2}^{\prime}$ we gain a power of $x$ and a power of $t^{-\alpha}$.

It remains to estimate $I_{1}$. First, we note that the leading order of the commutator $i\left[V,F_{2}\right]$ is given by

$$\frac{\partial V}{\partial x}\cdot\frac{\partial F_{2}}{\partial p}=\left(\frac{\partial V}{\partial x}\cdot\tilde{F}_{2}\right)t^{-\alpha};\quad\tilde{F}_{2}\equiv\tilde{F}_{2}(\frac{|p|}{t^{\alpha}}=1).$$

Hence, as before, we can estimate

$$\frac{\partial V}{\partial x}\cdot\frac{\partial F_{2}}{\partial p}p^{2}\sim t^{\alpha}(-\Delta+i)(-\Delta+i)^{-1}\frac{\partial V}{\partial x}\cdot\tilde{F}_{2}(-\Delta+i)^{-1}(-\Delta+i)$$

Assuming we have enough decay in $x$ for $\frac{\partial V}{\partial x}$, we gain from each side $t^{-2\alpha}$. So, all together, such terms decay like

$$t^{-\alpha}\cdot t^{2\alpha}t^{-4\alpha}=t^{-3\alpha}$$

where $t^{-\alpha}$ comes from the commutator, $t^{2\alpha}$ comes from $p^{2}\tilde{F}_{2}^{\prime}$, $t^{-4\alpha}$ comes from each factor of $(-\Delta+c)^{-1}\tilde{F}_{2}^{\prime}$. For this we need that $\frac{\partial V}{\partial x}$ decays at least like $\langle x\rangle^{-4-\varepsilon}$. Using a similar improvement of the estimates of $I_{2}$ and $I_{3}$, we conclude that:

$\partial h(t)\sim t^{-3\alpha}c\left(\left\|\psi_{0}\right\|_{H^{1}}\right)$. Integrating both sides, and remembering that $2\langle H(t)\rangle\sim\langle h(t)\rangle+O\left(t^{+2\alpha}\right)$, we choose $\alpha$ such that

$$t^{2\alpha}\sim\int^{t}s^{-3\alpha}ds\sim t^{-3\alpha+1}$$

This gives $2\alpha=1-3\alpha$ or $\alpha=1/5$.

Consequently, we proved under the assumption of $\langle x\rangle^{-6-\epsilon}$ decay for the derivatives of $V$, that

$$\left\langle p^{2}\right\rangle_{t}\leqslant t^{2/5}$$

∎

4 Incoming and outgoing Waves

We begin with the proof, of general interest, that the incoming part of the solution is smooth (as $\mathrm{H}^{2}$ in our case). Let as before, $A=\frac{1}{2}(x\cdot p+p\cdot x)$ and $M,R\gg 1$.

Then we define by the spectral theorem the following operators Soffer (2011):

$$P^{+}(A)=P_{M,R}^{+}=\frac{1}{2}\left(I+\tanh\left(\frac{A-M}{R}\right)\right)$$

$$P^{-}(A)=\frac{1}{2}\left(I-\tanh\left(\frac{A+M}{R}\right)\right)$$

$$P^{0}(A)=I-P^{-}(A)-P^{+}(A)=\frac{1}{2}\left(\tanh\left(\frac{A+M}{R}\right)-\tanh\left(\frac{A-M}{R}\right)\right)$$

Proposition 4.1 (Local Decay-Incoming Waves).

$$\left\|P^{-}(A)e^{i\Delta\tau}\chi(|p|=K)f(x)\right\|_{L^{2}}\leq\langle K\tau\rangle^{-2}\left\|\chi\langle x\rangle^{\sigma}f\right\|_{H^{2}},\text{with}\,\,\tau>0,K\tau\gg 1.\quad\sigma>2.$$ (3)

Proof.

$$P^{-}(A)e^{i\Delta\tau}\chi f(x)=P^{-}(A)e^{i\Delta\tau}\left(A^{2}+1\right)^{-1}\left(A^{2}+1\right)\chi f$$

$$\left(A^{2}+1\right)\chi f=\tilde{\chi}\left(A^{2}+1\right)\chi f+\left[A^{2},\tilde{\chi}\right]f$$

$$=\tilde{\chi}\left(\left(x\cdot\nabla_{x}\right)^{2}+\left(x\cdot\nabla_{x}+C\right)\right)\chi f+O(K^{-2})\tilde{\chi}^{\prime}(\Delta-1)\chi f.$$

Using the assumption that $f\in H^{2},\langle x\rangle^{2}f\in H^{2}$ the RHS is bounded by

$$c\|\chi f\|_{H^{2}}$$

where $\chi=\chi(|p|\sim k),\quad\tilde{\chi}\chi=\chi$. So, it remains to estimate

$$P^{-}(A)e^{i\Delta\tau}\chi\left(A^{2}+1\right)^{-1}$$

We now move $\left(A^{2}+1\right)^{-1}$ to the Left:

$${\left[\chi,\left(A^{2}+1\right)^{-1}\right]=-\left(A^{2}+1\right)^{-1}\left[\chi,A^{2}+1\right]\left(A^{2}+1\right)^{-1}}$$

$$i\left[A^{2},\chi\right]=Ai[A,\chi]+i[A,\chi]A$$

$$=Ap\chi^{\prime}(p)+p\chi{\prime}(p)A=2p\chi{\prime}(p)A$$

$$+\left(\chi{\prime}(p)+p^{2}\chi^{\prime\prime}\right)A$$

Therefore

$$\begin{gathered}\chi\left(A^{2}+1\right)^{-1}=\left(A^{2}+1\right)^{-1}\chi-\left(A^{2}+1\right)^{-1}G(p)A\left(A^{2}+1\right)^{-1}\\ G(p)=O(1).\end{gathered}$$

So, it is sufficient to estimate

$$P^{+}(A)e^{-i\Delta\tau}\left(A^{2}+1\right)^{-1}\chi(p)=$$

$$=P^{+}(A)\left(\left(A-2p^{2}\tau\right)^{2}+1\right)^{-1}e^{-i\Delta\tau}\chi(p)$$

Formally, since $P^{-}(A)$ projects on $A\leq 0$ up to small corrections,then $A-2p^{2}\tau<0$, so at the "symbol" level the result is clear. However, $P^{-}(A)$ is not a $\psi$DO, and we cannot use symbol calculus. First, we approximate $A-2p^{2}\tau$ by

$$P^{-}(A)A-2p^{2}\tau:$$

$${\left[\left(A-2p^{2}\tau\right)^{2}+1\right]^{-1}-\left[\left(P^{-}(A)A-2p^{2}\tau\right)^{2}+1\right]^{-1}}$$

$$= {\left[\left(P^{-}A-2p^{2}\tau\right)^{2}+1\right]^{-1}\left\{\left(P^{-}A+P^{0}A\right)\left(A-4p^{2}\tau\right)+\left(P^{-}A+P^{0}A\right)A\right\}\times}$$

$${\left[\left(A-2p^{2}\tau\right)^{2}+1\right]^{-1}}$$

The term $P^{-}A\left[\left(P^{-}A-2p^{2}\tau\right)^{2}+1\right]^{-1}$ satisfies the estimate we want, being bounded by

$$\left[\left(1+2p^{2}\tau\right)^{2}+1\right]^{-1}.$$

To control the difference between the resolvents of

$$A_{\tau}^{2}\equiv\left(A-2p^{2}\tau\right)^{2}\text{ and }\left(P^{-}A-2p^{2}\tau\right)^{2}\text{, }$$

we commute the $P^{-}(A)$ through, and use that $P^{+}(A)P^{-}(A)$ is exponentially small $P^{+}(A)P^{-}(A)\sim e^{-|A|/R}$ for $|A|$ large, $P^{+}(A)P^{0}(A)$ is exponentially localized near $A=+M$.

So

$$P^{-}(A)\left\{\left[\langle A_{\tau}\rangle^{2}+1\right]^{-1}-\left[\left\langle P^{-}(A)A_{\tau}\right\rangle^{2}+1\right]^{-1}\right\}$$ (4)

$$=\left(\langle A_{\tau}\rangle^{2}+1\right)^{-1}\left[P^{-}(A)\left(P^{-}+P^{0}\right)4Ap^{2}\tau+\text{ h.o.t. }\right]\left(\left\langle A_{\tau}P^{-}\right\rangle^{2}+1\right)^{-1}$$ (5)

$$+\left[P^{-}(A),\left[\langle A_{\tau}\rangle^{2}+1\right]^{-1}\right]\left[\left\langle P^{-}+P^{0}\right)4Ap^{2}\tau+\text{ h.o.t. }\right]\left(\left\langle A_{\tau}P^{-}\right\rangle^{2}+1\right)^{-1}$$ (6)

$$=\left(\langle A_{\tau}\rangle^{2}+1\right)^{-1}O\left(e^{-|A-M|/R}M\right)4p^{2}\tau\left(\left\langle A_{\tau}P^{-}\right\rangle^{2}+1\right)^{-1}$$ (7)

$$\left.+\left[P^{-}(A),\left[\langle A_{\tau}\rangle^{2}+1\right]^{-1}\right]\left[\left(P^{-}+P^{0}\right)4Ap^{2}\tau+\text{ h.o.t. }\right]\left(\langle A_{\tau}P^{-}\right\rangle^{2}+1\right)^{-1}.$$ (8)

The first term on the RHS is bounded by

$$O\left(e^{-|A-M|/R}M\right)\min\left\{\frac{1}{1+p^{2}\tau}+\frac{1}{1+\langle P^{-}A\rangle}\right\}$$

The second term on the RHS is more difficult to control, since commuting a projection is creating an error term which is not a projection on the other side. This is the reason why one has to use a special form of the projection; see Soffer (2011). To this end, we use the following identity: Let $A$ be the dilation generator, and

$$U_{\theta}=e^{-i\theta A}$$

the group of dilations on $L^{2}\left(\mathbb{R}^{n}\right),\quad\theta$-real. Then for a general operator $B$ on $L^{2}\left(\mathbb{R}^{n}\right),$ we define

$$U_{\theta}BU_{-\theta}\equiv B_{\theta}.$$

We will now extend these identities to the complex $\theta=\theta_{0}+i\theta_{1}\quad\left|\theta_{1}\right|<\delta$, provided that $B$ is a "dilation analytic", in the sense that $B_{\theta}$ makes sense on some common domain $D_{\theta}\subset L^{2}$, for all $\operatorname{such}\theta$.

We then have

Proposition 4.2 (Soffer (2011)).

Let $B$ be an operator such that $B_{\theta}$ is dilation-analytic in $|\Im\theta|<\delta.$ Then,

$$i[P^{-}(A),B]=\cosh^{-1}\frac{(A+M)}{R}[B_{-i/R}-B_{i/R}]\cosh^{-1}\frac{(A+M)}{R}$$ (9)

for all $R>c=\mathcal{O}(1).$

Now we can control the commutator term with $P^{-}(A).$ Clearly, $\langle A-2p^{2}\tau\rangle^{-2}$ is dilation analytic for small $\theta$, since ($\beta$ - real)

$$e^{\beta A}(2p^{2}\tau+1)^{-2}e^{-\beta A}=(2e^{2\beta i}p^{2}\tau+1)^{-2}.$$

$$\tau e^{2\beta i}p^{2}+1=\tau\cos{2\beta}p^{2}+1+i\sin{2\beta}p^{2}\tau,$$

and for $\beta$ sufficiently small, the real part is larger than 1, for $\tau\geq 0.$

Then,

$$1+(A-\tau e^{2\beta i}p^{2})^{2}=1+(A-2p^{2}\tau)^{2}$$

$$-\epsilon(\beta)[Ap^{2}\tau+p^{2}\tau A]+4\cos{4\beta}p^{4}\tau^{2}$$

$$+i[4A\sin{2\beta}p^{2}\tau]+4\sin{4\beta}p^{4}\tau^{2}]\neq 0.$$

Here $\epsilon(\beta)=2-2\cos{2\beta}>0$ for all $\beta\neq 0$ small.

The sign of $Ap^{2}+p^{2}A$ is the same as the sign of $A$, since $Ap^{2}+p^{2}A=2pAp.$

Consequently, the commutator term with $P^{-}(A)$ in (4) using equation (9), is controlled by the norm of

$$\cosh^{-1}{(\beta A)}(P^{-}+P^{0})Ap^{2}\tau\left(\langle A_{\tau}P^{-}\rangle^{2}+1\right)^{-1}$$

$$\leq\frac{M}{\cosh{(\beta A)}}\left(\frac{1}{p^{2}\tau+1}+\frac{1}{\langle AP^{-}\rangle+1}\right).$$

Finally, we note that the resulting estimate of

$$P^{-}(A)[\langle A_{\tau}\rangle^{2}+1]^{-1}\leq\left(\frac{1}{p^{2}\tau+1}+\frac{1}{\langle AP^{-}\rangle+1}\right)$$

can be used to control terms linear in $A$, since they are multiplied from the left by $P^{-}(A)[\langle A_{\tau}\rangle^{2}+1]^{-1}.$ This gives the final statement of the proposition. ∎

We can now use the proposition on Local Decay of Incoming Waves (3) to control the Duhamel term in $H^{2}$, on the range of $P^{-}(A)$. Since the initial data are in $H^{2}$, we only need to control the Duhamel term.

Integrating by parts, we get

$$-i\Delta \int_{0}^{t-\varepsilon}e^{i\Delta(t-s)}V(s)\psi(s)ds$$

$$= e^{i\Delta\varepsilon}V(t-\varepsilon)\psi(t-\varepsilon)-V(t=0)\psi(t=0)$$

$$-\int_{0}^{t-\varepsilon}e^{i\Delta(t-s)}\left(\frac{\partial V}{\partial s}\psi(s)+iV\Delta\psi(s)-iV^{2}\psi(s)\right)ds$$

So, we only need to show that the RHS is uniformly bounded in $L^{2}$. The first and second terms are bounded in $L^{2}$ if $\langle x)^{m}V(x,s)$ is bounded from $H^{2}\rightarrow L^{2}$, since we proved that $\langle x\rangle^{-\sigma}\psi(s)$ is in $H^{2}$. Similarly $\frac{\partial V}{\partial s}\psi-iV\Delta\psi-iV^{2}\psi(s)$ are in $L^{2}$. So, we need to prove the integrability in $s_{\text{. }}$ For this, we use

$$\left\|P^{-}(A)e^{i\Delta(t-s)}g(|p|\geq\delta)\langle x\rangle^{-\sigma}\right\|\leqslant c\langle t-s\rangle^{-2},$$

for $\delta>0$.

We do not have the case $|p|<\delta$, since in this region the solution is in $\mathrm{H}^{s}$ for all $s>0$. The integral in the interval $s\in(t-\varepsilon,t)$ is trivially bounded in $L^{2}$. We conclude:

Theorem 4.3.

Incoming Waves

If the initial data is in $\mathrm{H}^{2}$, then

$$P^{-}(A)\psi(t)\in H^{2}\text{, uniformly in }t\text{. }$$ (10)

5 More On Incoming Waves

As we saw in the previous section, most of the correction to the free flow is outgoing. We will further control the Incoming Waves by PROBs. To control the incoming waves, we need a PRES that registers the fact that solutions become outgoing for large times. The first PROB we use is $F\left(\frac{A+M}{R}\leq-1\right)$.

Proposition 5.1.

The $H^{1}$ norm is bounded for all times, on average, on the negative spectral part of $A$ (incoming).

Proof.

$$\partial_{t}\left\langle F\left(\frac{A+M}{R}\leqslant-1\right)\right\rangle=\frac{1}{R}\left\langle p\psi,\tilde{F}^{\prime}p\psi\right\rangle+\langle i[V,F]\rangle$$ (5,1)

The first term is negative, so we get from (5,1)

$$\langle F\rangle_{T} -\langle F\rangle_{0}+\frac{1}{R}\int_{0}^{T}\langle p\psi,|\tilde{F}\left(\frac{A+M}{R}\sim-1\right)|\,p\psi\rangle ds$$ (11)

$$\lesssim\left|\int_{0}^{T}[V,F]ds\right|$$ (5,2)

We get an estimate of the first derivative in the region $A\sim-M$ (in an interval of size $R$ ), in terms of the potential. Next, we show that the potential term is of higher order in $\left(\frac{1}{R}\right)$ : Write $V$ as $\langle x)^{-\sigma}\langle x\rangle^{\sigma}\vec{V}\langle x\rangle^{\sigma}\langle x\rangle^{-\sigma}\equiv\langle x\rangle^{-\sigma}\tilde{V}\langle x\rangle^{-\sigma}$. Then

$${[V,F]=\left[\langle x\rangle^{-\sigma},F\right]\tilde{V}\langle x\rangle^{-\sigma}+\langle x\rangle^{-\sigma}\tilde{V}\left[\langle x\rangle^{-\sigma},F\right]}$$

$$+\langle x\rangle^{-\sigma}[\tilde{V},F]\langle x\rangle^{-\sigma}$$

$$=\frac{1}{R}\langle x\rangle^{-\sigma}\tilde{F}^{\prime}\tilde{V}\langle x\rangle^{-\sigma}+\langle x\rangle^{-\sigma}\tilde{V}\tilde{F}^{\prime}\langle x\rangle^{-\sigma}\frac{1}{R}$$

$$+\langle x\rangle^{-\sigma}[\tilde{V},F]\langle x\rangle^{-\sigma}$$ (5,3)

Here, we used $[x,F]=\tilde{F}^{\prime}x=x\tilde{F}^{\prime}$ where $\tilde{F}^{\prime}$ stands for a discrete derivative of $F$ in the $\pm i$ direction.

We get another power of $\left(\frac{1}{M}\right)\left(<\frac{1}{R}\right)$ using

$$F(A+i)^{-1}\sim\frac{1}{M}F$$

Then, write $F=F(A+i)^{-1}(A+i)$

$$=\tilde{F}\frac{1}{M}(x\cdot p+c).$$

We can write all terms as

$$\frac{1}{MR}\left\langle p\langle x\rangle^{-\sigma}\psi,O(1)p\langle x\rangle^{-\sigma}\psi\right\rangle.$$

In particular, the term, by the Commutator Expansion Lemma,

$$\langle x\rangle^{-\sigma}[\tilde{V},F]\langle x\rangle^{-\sigma}=$$

$$=\langle x\rangle^{-\sigma}(x\cdot\nabla\tilde{V})\cdot\tilde{F}^{\prime}\frac{1}{R}\langle x\rangle^{-\sigma}$$

$$+\langle x\rangle^{-\sigma}\frac{1}{R^{2}}O_{2}(F)\langle x\rangle^{-\sigma}$$

$$=\frac{1}{R}\langle x\rangle^{-\sigma}(x\cdot\nabla\tilde{V})(A+i)(A+i)^{-1}\tilde{F}^{\prime}\langle x\rangle^{-\sigma}$$

$$+\langle x\rangle^{-\sigma}\frac{1}{R^{2}}O_{2}(F)\langle x\rangle^{-\sigma}$$

$$=\frac{1}{MR}P\langle x\rangle^{-\sigma}(x\cdot\tilde{V})\langle x\rangle^{-\prime}\langle x\rangle^{-\sigma}$$

$$+\langle x\rangle^{-\sigma}\frac{1}{R^{2}}\mathcal{O}_{2}(F)\langle x\rangle^{-\sigma}.$$

We need to know that $\mathcal{O}_{2}(F)$ is bounded.

$$\mathcal{O}_{2}(F)=\int{\hat{F}(\lambda)}e^{-i\lambda A}d\lambda\int_{0}^{\lambda}dse^{isA}\int_{0}^{s}due^{-iuA}(x\cdot\nabla)^{2}\tilde{V}e^{iuA}.$$ (12)

So, an immediate bound in $L^{2}$ is obtained if $\|x^{2}\Delta\tilde{V}\|_{L^{\infty}}<c<\infty.$

But in fact, a weaker condition is sufficient:

$$\int_{0}^{s}due^{-iuA}(x\cdot\nabla)^{2}\tilde{V}e^{iuA}\langle p\rangle^{-1}$$ (13)

be bounded.

The above integral reduces to

$$\int_{0}^{s}dy\frac{f(y)}{y}\langle p\rangle^{-1};\quad f(z)\equiv(z\cdot\nabla_{z})^{2}\tilde{V}(z,t).$$ (14)

Near $x=0$, $\langle p\rangle^{-1}$ bounds $\frac{1}{|x|},$ so the singularity near zero is reduced to $|z|^{2}(\Delta\tilde{V}(z)).$ Up to this point we have an estimate for a fixed $A\sim M.$

We now iterate this estimate to higher Sobolev norm.

We use as PROB

$$\left\langle p\psi(t),AF\left(\frac{A+M}{R}\leq-1\right)p\psi\right\rangle\leq 0.$$ (15)

In this case, the RHS, similar to previous estimate, is given by

$$\int_{0}^{T}\left\langle p^{2}\psi(t),\tilde{F^{\prime}}\left(\frac{A+M}{R}\leqslant-1\right)p^{2}\psi(t)\right\rangle dt$$

$$+\int_{0}^{T}\langle i[V,pAFp]\rangle dt$$

The estimate of the potential is as before. The commutator has two derivatives w.r.t. $x$. Hence, we can pull $(x)^{-\sigma}$ on both sides and one more derivative on each side; then we can use that locally in $x$, the $H^{2}$ norm is uniformly bounded. Therefore, again the potential term is bounded by $T.$ Notice that the LHS is integrated on $\left[0,T^{\prime}\right]$ s.t. the LHS is bounded for such $T^{\prime}$, using the previous estimate.

Upgrading the Estimates We will sum over dyadic intervals covering the negative real line, starting from $-M_{0}.$ Multiplying by $\left(R_{0}2^{n}\right)^{2}(n)^{-1}(n\geqslant 1)$, using for each $n=0,1,2,3...,$  $F_{n}(\frac{A+M_{0}+2^{n}R_{0}}{2^{n}R_{0}}\leq-1).$ After multiplying by $R_{0}^{2}2^{2n}/n$ we take the sum over $n.$

The resulting inequality becomes:

$$\langle\psi(T),A^{2}F(\frac{A+M_{0}}{R_{0}}\leq-1)\psi(T)\rangle$$ (16)

$$-\langle\psi(0),A^{2}F(\frac{A+M_{0}}{R_{0}}\leq-1)\psi(0)\rangle$$ (17)

$$+\int_{0}^{T}\langle p\psi(s),|F(\frac{A+M_{0}}{R_{0}}\leq-1)A\ln^{-1}{\langle A\rangle}|p\psi(s)\rangle ds$$ (18)

$$\approx\int_{0}^{T}\langle p\psi(s)\langle x\rangle^{-\sigma}F(\frac{A+M_{0}}{R_{0}}\leq-1)\ln^{-1}{\langle A\rangle}\langle x\rangle^{-\sigma}p\psi(s)\rangle ds$$ (19)

$$\leq c\|\psi(0)\|^{2}_{H^{2}}T.$$ (20)

Here, we use that

$$\sum_{n}F_{n}\geq F(\frac{A+M_{0}}{R_{0}}\leq-1)$$ (21)

$$\sum_{n}2^{kn}R_{0}^{k}F_{n}\simeq A^{k}F(\frac{A+M_{0}}{R_{0}}\leq-1)\ln{|A|}.$$ (22)

See Liu and Soffer (2025). We therefore conclude that the average over time of

$$\langle p\psi(s),|A|\ln^{-1}{\langle A\rangle}F(\frac{A+M_{0}}{R_{0}}\leq-1)p\psi(s)\rangle$$ (23)

is bounded.

Upgrading to pointwise estimate

To obtain estimates that are better than those bound by $T$, we need to estimate the potential term in terms of the leading term on the RHS of the PRES.

The leading order term is of the form:

$$\int_{0}^{T}\left\langle p\psi,\frac{A}{\ln\langle A\rangle}F\left(\frac{A+M_{0}}{R}\leq-c\right)p\psi\right\rangle_{s}ds$$

(This is derived by multiplying by $R2^{n}/n$ and summing).

The potential term (Assuming first that $V$ is dilation analytic) is:

$$\left.\int_{0}\langle p\psi,x\cdot\tilde{V}|x|^{2}F_{A}\left\langle\frac{A+M_{0}}{R}\leqslant-c\right)p\psi\right\rangle ds$$

Hence $i\left[V,F_{A,M}\right]=c\frac{1/R}{\cosh\left(\frac{A+M}{R}\right)}x\cdot\tilde{V}\frac{1}{\cosh\left(\frac{A-M}{R}\right)}$ with

$$x\cdot\tilde{V}\equiv V_{iR^{-1}}-V_{-iR^{-1}}=V\left(\frac{i|x|}{R}\right)-V\left(\frac{-i|x|}{R}\right)$$

Then, we calculate (using the notation $\operatorname{ch}\equiv\cosh\left(\frac{A+M}{R}\right)$)

$$\left(\frac{1}{\cosh\left(\frac{A+M}{R}\right)}x\cdot\tilde{V}\frac{1}{\cosh\left(\frac{A+M}{R}\right)}\right)=$$

$$= \frac{1}{\operatorname{ch}}\frac{1}{A+i}(A+i)x\cdot\tilde{V}(A+i)\frac{1}{A+i}\frac{1}{ch}$$

$$=$$

$$= c\frac{1}{(i+A)ch}x\cdot\tilde{V}|x|^{2}\frac{1}{(i+A)ch}+\frac{p}{(i+A)ch}x\cdot\tilde{V}|x|^{2}\frac{1}{(i+A)ch}$$

$$+\left\{\frac{1}{\{(i+A)ch}\left[p,A\,ch\right]\frac{1}{(i+A)ch}x\cdot\tilde{v}|x|^{2}\frac{1}{(i+A)ch}+\text{ similar }\right\}$$

Using $\frac{1}{\operatorname{ch}^{2}}=\frac{1}{\cosh^{2}\left(\frac{A+M}{R}\right)}$ localizes $A$ around $-M$ with width $R$, multiplication by $R^{2}2^{n+2}/n$ of the quantity

$$\frac{1}{\cosh^{2}\left(\frac{A+M+2^{n}R}{R}\right)},$$

and summing over $n=0,1,2\ldots$, this sum is essentially

$$(\ln\langle A\rangle)^{-1}(A+i)^{-2}A^{2}F_{A}\left(\frac{A+M}{R}\leq-C\right).$$

So, we project on $\left|\frac{A}{\ln\langle A\rangle}\right|>K\gg\left\||x|^{2}x\cdot\tilde{V}\right\|_{L^{\infty}}$ and we find that the leftover term is

$$\int_{0}^{T}\left\langle p\psi,F_{A}F_{1}E(|A|\leqslant K)\right||x|^{2}x\cdot\tilde{V}|p\psi\rangle_{s}ds$$

But for $K\ll M_{0}$, we have that

$$F_{A}\left(\frac{A+M_{0}}{R}\leq-c\right)F_{A}(|A|\leq K)\approx e^{-M_{0}}.$$

So, the leftover term is bounded by

$$ce^{-M_{0}}T_{0}.$$

We therefore get the following estimate

$$\left\langle\psi(t),{A}^{2}F\left(\frac{A+M_{0}}{R}\leq-C\right)\psi(t)\right\rangle_{T}-\left\langle\psi(0),A^{2}F_{A}\psi(0)\right\rangle_{0}$$ (25)

$$=\int_{0}^{T}\left\langle p\psi(s),(\ln\langle A\rangle)^{-1}(A-\tilde{W}(x))F(|A|>K)F_{A}p\psi(s)\right\rangle ds$$ (26)

$$+\int_{0}^{T}\left\langle p\psi(s),F_{A}\tilde{W}F(|A|\leq K)p\psi(s)\right\rangle ds+h.o.t$$ (27)

The last term on the RHS is bounded by

$$ce^{-M_{0}}T$$

The first term on the RHS is negative since

$$(|A|-\tilde{W})F_{K}(A>K)\sim\tilde{F}_{K}(|A|-\tilde{W})\tilde{F}_{K}\geqslant\theta\tilde{F}_{K}^{2}$$

where $\tilde{F}_{K}^{2}=F_{K}$ and we used that

$$K\gg\||x|^{2}x\cdot\tilde{V}\|_{L^{\infty}}.$$

The h.o.t. come from commutators. Such commutators are of higher order in $(1/R)$ and/or $A,$ as $[p,A]\sim p;[x,F(A/R)]\sim x\tilde{F}^{\prime}/R$.

Up to this point, we used dilation analyticity of the potential to write

$$i[V,F]\sim\frac{1}{ch}x\cdot\tilde{V}\frac{1}{ch}.$$

However, it is not necessary:

$$\langle\psi,i[V,F_{A}]\psi\rangle=\langle\psi,iVF_{A}\psi\rangle-c.c.$$ (28)

$$=\langle\psi,iVF_{A}G_{A}(|A|\leq M_{0}/2)\psi\rangle+\langle\psi,iVF_{A}G_{A}(A\leq-M_{0}/2)\psi\rangle$$ (29)

$$=O(e^{-M_{0}/2})+c\langle p^{2}\psi,V|x|^{4}F_{A}(A+i)^{-4}p^{2}\psi\rangle$$ (30)

$$+c\langle\psi,\tilde{V}F_{A}(A+i)^{-4}\psi\rangle$$ (31)

$$+O(e^{-M_{0}/2})\leq c\|\langle x\rangle^{-\sigma}p^{2}\psi\|\|\langle x\rangle^{-\sigma}F_{A}p^{2}\psi\|M_{0}^{-4}+O(e^{-M_{0}/2})$$ (32)

The integral over time is then controlled by

$$cTe^{-M_{0}/2}$$ (33)

$$+cM_{0}^{-4}\int_{1}^{T}\|\langle x\rangle^{-\sigma}p^{2}\psi\|_{L_{t}^{\infty}L^{2}_{x}}\|\langle x\rangle^{-\sigma}F_{A}p^{2}\psi\|_{L_{t}^{\infty}L^{2}_{x}}ds.$$ (34)

If we now choose $M_{0}\sim t^{\alpha}$, we get the bound

$$\sim\left(\int_{1}^{T}s^{-4\alpha}ds\right)\left\|\langle x\rangle^{-\sigma}F_{A}p^{2}\psi\right\|_{L_{s}^{\infty}L_{x}^{2}}^{2}$$

$$\leq C$$

where $C<\infty$, for

$$\int_{1}^{T}s^{-4\alpha}ds<\infty\quad(\alpha>1/4).$$

We conclude that ((25))

$$\left\langle F_{A}\left(\frac{A+t^{\alpha}}{R}\leq-1\right)A^{2}\right\rangle\approx 1.$$

We conclude that if $|x|\geq 1$, then $\left.\left.\left\langle F_{A}\right|p\right|^{2}\right\rangle\leq 1$. It remains to show that

$$\left\|F_{A}\left(\frac{A+t^{\alpha}}{R}\geqslant-1\right)\Delta\psi(t)\right\|_{2}\sim 0(1)$$

∎

Remark 5.2.

Since we use an analytic PROB, we do not have symmetrization error terms. Hence, all error is of order $M_{0}^{-4}.$ Therefore, multiplying the entire estimate by $M_{n}^{4}$ and summing over $n$, we get a bound on $A^{5-0},$ which for $|x|>1$ also is a bound on the norm $H^{5/2-0}$. We see the expected extra smoothing in the region away from the Propagation Set.

6 Control of the solution on the Propagation Set

High frequency and long-time parts of the solution of the free flow concentrate on the Propagation Set:

$$PS_{E}=\left\{\left.\frac{x}{t}=2p\right\rvert\,\quad p^{2}=E\right\}$$

for the part of the solution localized at energy $E.$ Sharper localization of the solution is constructed by projection on the region $|x-2pt|\leq t^{\alpha}$, $0\leqslant\alpha<1$.

Consider the following $PROB$ :

$$\partial_{t}\left\langle\psi(t),\left(HF_{c}+F_{c}H\right)\psi(t)\right\rangle=\left\langle\psi(t),\left(\frac{\partial V}{\partial t}F_{c}+F_{c}\frac{\partial V}{\partial t}\right)\psi(t)\right\rangle$$ (35)

$$+\left\langle\psi(t),\left(H\frac{\partial F_{c}}{\partial t}+\frac{\partial F_{c}}{\partial t}H\right)\psi(t)\right\rangle$$ (36)

$$+\left\langle\psi(t)\left[Hi\left[V,F_{c}\right]+i\left[V,F_{c}\right]H\right]\psi(t)\right\rangle$$ (37)

$$\text{ Here }H=H_{0}+V(t)=p^{2}+V(t)$$ (38)

$$F_{c}=F_{c}\left(\frac{|x-2pt|}{t^{\alpha}}\leq 1\right)$$ (39)

Proposition 6.1.

The PROB mentioned above ((35)) is uniformly bounded in time.

Proof.

First, we use $F_{c}=F_{c}\left(\frac{|x-2pt|}{s}\leq 1\right.)$

Then, the first term on the right of (35) is of the form

$$\sim\left\langle\langle x\rangle^{-\sigma}F_{c}\right\rangle+c\cdot c.$$

The second term is zero:

$$D_{t}F_{c}=i\left[-\Delta,F_{c}\right]+\frac{\partial F_{c}}{\partial t}=0\quad\textbf{since}\,D_{t}|x-2pt|=0.$$

The third term is of the form

$$\left\langle p^{2}\tilde{V}\tilde{F}_{c}\right\rangle+c.c.$$

So, we need to prove the integrability in time of the above term (depending on $s$ ).

For $\sigma=2$, we get

Lemma 6.2.

$$\langle x\rangle^{-2}\langle x-2pt\rangle^{2}F_{c}\lesssim cF_{c}+O\left(\frac{s}{\langle x\rangle}F_{c}\right)+O\left(\frac{s^{2}}{\langle x\rangle^{2}}F_{c}\right)$$

Hence,

$$\langle x\rangle^{-2}\left(4p^{2}t^{2}\right)F_{c}\leq\langle x\rangle^{-2}\left(4p^{2}t^{2}-x^{2}+x^{2}\right)F_{c}$$

$$\leq F_{c}+\langle x\rangle^{-2}\left((2pt-x)\cdot(2pt+x)+x^{2}-2pt\cdot x\right.$$

$$-x\cdot 2pt)F_{c}$$

$$\leq F_{c}\langle x\rangle^{-2}\left(s(2pt-x)+2|x|s+x^{2}c+|x|s\right)F_{c}$$

$$\leq F_{c}\left[\left(s^{2}/\langle x)^{2}\right)+s/\langle x\rangle+1\right]$$

So,

$$\langle x\rangle^{-2}F_{c}=\langle x\rangle^{-2}\left(1+4p^{2}t^{2}\right)^{-1}\left(1+4p^{2}t^{2}\right)F_{c}$$

$$=\left(1+4p^{2}t^{2}\right)^{-1}\langle x\rangle^{-2}\left(1+4p^{2}t^{2}\right)F_{c}$$

$$\quad-\langle x\rangle^{-2}\left[\langle x\rangle^{2},\left(1+4p^{2}t^{2}\right)^{-1}\right]\langle x\rangle^{-2}\left(1+4p^{2}t^{2}\right)F_{c}$$

The first term on the RHS is bounded by using $A=x\cdot(p-x/(2s))+x^{2}/(2s)$ and using

$$\left(1+4p^{2}t^{2}=1+|x-2pt|^{2}-2x^{2}+8A\right)$$

$$\left(1+4p^{2}t^{2}\right)^{-1}\left[2+\frac{s}{\langle x\rangle}+\frac{s^{2}}{\langle x\rangle^{2}}\right]F_{c}$$

The second term is also bounded by

$$\left(1+4p^{2}t^{2}\right)^{-1}\left[1+\frac{s}{\langle x\rangle}+\frac{s^{2}}{\langle x\rangle^{2}}\right]F_{c}$$

since

$$\begin{gathered}\langle x\rangle^{-2}\left[x^{2},\left(1+4p^{2}t^{2}\right)^{-1}\right]\langle x\rangle^{-2}=\\ \langle x\rangle^{-2}\left(1+4p^{2}t^{2}\right)^{-1}\left[x^{2},4p^{2}t^{2}\right]\left(1+4p^{2}t^{2}\right)^{-1}\langle x\rangle^{-2}\\ -i\left[x,4p^{2}t^{2}\right]=8pt^{2}-\left[x\left[x,4p^{2}t^{2}\right]\right]=8t^{2}\end{gathered}$$

Then

$$\langle x\rangle^{-2}\left(1+4p^{2}t^{2}\right)^{-1}\left(x\cdot pt^{2}+t^{2}\right)\left(1+4p^{2}t^{2}\right)^{-1}\langle x\rangle^{-2}$$

We need to control $pt^{2}$ and $t^{2}$.

$$pt^{2}\left(1+4p^{2}t^{2}\right)^{-1}\langle x\rangle^{-2}=$$

$$\quad pt^{2}\left(1+4p^{2}t^{2}\right)^{-1}|p|\frac{1}{|p|^{2}}\langle x\rangle^{2}=O(1)$$

$$\langle x\rangle^{-2}\left(1+4p^{2}t^{2}\right)^{-1}t^{2}\left(1+4p^{2}t^{2}\right)^{-1}\langle x\rangle^{-2}$$

$$=\langle x\rangle^{-2}|p|^{-1}|p|\left(1+4p^{2}t^{2}\right)^{-1}t^{2}\left(1+4p^{2}t^{2}\right)^{-1}|p|\frac{1}{|p|}\langle x)^{-2}$$

$$=\langle x\rangle^{-1}O(1)\left(1+4p^{2}t^{2}\right)^{-1}\langle x\rangle^{-1}.$$

Next, we observe that for $|p|\geqslant 1$

$$\begin{gathered}|x-2pt|\leq s\\ \text{ with }s\equiv t^{\alpha},\alpha\leq 1\end{gathered}$$

imply, at the symbol level, that for $t\gg t_{0}$, $\langle x\rangle\gtrsim t$. Hence, we expect, as we have used before, that $\quad\langle x\rangle^{-1}sF_{c}$ and $\langle x\rangle^{-2}s^{2}F_{c}$ with $t$ large, and $s\sim t^{\alpha}$, that these expressions are bounded. To conclude

$$\left\langle H(t)F_{c}+F_{c}H(t)\right\rangle_{T}\lesssim\left\langle H(0)F_{c}+F_{c}H(0)\right\rangle$$

$$\quad+\int_{0}^{T}\left(1+t^{-2}O(1)\left\{\left\langle\psi,p^{2}\langle x)^{-\sigma+2}\langle x\rangle^{-2}\tilde{F}_{c}\psi\right\rangle+c.c\right\}dt\right.$$

Finally, we write the $p^{2}$ factor as

$$\left(4p^{2}-\frac{x^{2}}{t^{2}}\right)+\frac{x^{2}}{t^{2}}$$

$$4p^{2}-\frac{x^{2}}{t^{2}}=(2p-x/t)\cdot(2p-x/t+2x/t)=$$

$$=[2pt-x)\cdot(2pt-x)+(2pt-x)\cdot 2x]\frac{1}{t^{2}}$$

we then bound the $(2pt-x)$ factors by $s$, so the overall bound is (on support $\langle x\rangle^{-\sigma+}F_{c}$ )

$$s^{2}/t^{2}+\frac{s\langle x\rangle}{t^{2}}\leqslant\frac{s^{2}}{t^{2}}+\frac{s}{t^{2}}$$

So, by our choice of $s=T^{\alpha}$, we conclude that

$$\left\langle H(T)F_{c}(T)+F_{c}(T)H(T)\right\rangle\stackrel{{\scriptstyle\varepsilon}}{{\sim}}\left\langle H(0)F_{c}+F_{c}H(0)\right\rangle.$$

∎

7 Away From the Propagation Set

Theorem 7.1 (Away from PS).

We have the following bound away from the PS defined as $|x-2pt|\leq t^{\alpha}:$

$$\left\langle\left(\frac{A}{t^{\beta}}\right)^{l+1}\ln\frac{\langle A\rangle}{t^{\beta}}F_{A,C}\left(\frac{\left|2p^{2}t-A\right|}{t^{\alpha}}\geqslant 1\right)\right\rangle_{T}$$ (40)

$$\lesssim\langle\cdots\rangle_{1}+o(1).$$ (41)

Proof.

In this case, we first consider the following PROB:

$$B_{M}\equiv F_{A}\left(\frac{A}{tM}\leq\frac{1}{2}\right)F_{A,C}\left(\frac{2p^{2}t-A}{t^{\alpha}}\geqslant 1\right)+c.c$$ (B.1)

Note that

$$F_{A,C}=e^{+i\Delta t}F_{A}\left(-\frac{A}{t^{\alpha}}\geqslant 1\right)e^{-i\Delta t}.$$ (B.2)

The Heisenberg Derivative w.r.t. $-\Delta$ is negative to leading order:

$$D_{H_{0}}F_{A}=\tilde{F}_{A}^{\prime}\left(\frac{2p^{2}}{Mt}-\frac{A}{Mt\cdot t}\right)$$ (B.3)

This is an identity, since we use an analytic form of $F_{A}$. On the support of $F_{A,C}$ this derivative is

$$\leqslant\tilde{F}_{A}^{\prime}\frac{t^{\alpha}}{Mt}F_{A,C}+c.c.$$

to leading order. The Heisenberg derivative of $F_{A,C}$ is also negative:

$$D_{H_{0}}F_{A,C}=\partial_{t}F_{A,C}=-\frac{\alpha}{t}F_{A,C}^{\prime}.$$ (B.4)

The commutator with the potential is also fast decaying, since on the boundary $\tilde{F}_{A}^{\prime}$, $A\sim Mt$, and on boundary of $\tilde{F}_{A,C}^{\prime}$, we are on the Propagation Set, and $|x|\sim 2pt.$ Now, we note that far away from the PS (prop. set)

$$\frac{2p^{2}}{Mt}\gg\frac{A}{M^{2}t},$$

showing that away from the $PS$, we get a good bound on $p^{2}$. We want to improve this to $p^{4}$, by finding the decay on $M$ of the higher order terms. One such term comes from the potential. With sufficient decay in $|x|$ large, these terms are higher order. The main contribution comes from Symmetrization.

The leading order commutator is negative modulo correction terms coming from symmetrizing. We use, for $A,B$ non-negative operators

$$AB+BA=2A^{1/2}BA^{1/2}+A^{1/2}\left[A^{1/2},B\right]$$

$$+\left[B,A^{1/2}\right]A^{1/2}=2A^{1/2}BA^{1/2}+\left[A^{1/2},\left[A^{1/2},B\right]\right].$$ (B.5)

Using eq. (B.5), we get

$$\tilde{F}_{A}^{\prime}\left(2p^{2}t-A\right)F_{A,C}+F_{A,C}\left(2p^{2}t-A\right)\tilde{F}_{A}^{\prime}$$

$$=2\sqrt{\tilde{F}_{A}^{\prime}}\left[\left(2p^{2}t-A\right)F_{A,C}\right]\sqrt{\tilde{F}_{A}^{\prime}}+$$

$${\left[\sqrt{\tilde{F}_{A}^{\prime}},\left[\sqrt{\tilde{F}_{A}^{\prime}},G_{A,C}\left(2p^{2}t-A\right)\right]\right]\text{ (B.6a) }}$$

$$G_{A,C}=\left(2p^{2}t-A\right)\cdot F_{A,C}\quad(B.6b)$$ (B.6)

By the commutator expansion formula, we have

$${\left[f\left(\frac{A}{t}\right),g\cdot\left(\frac{2p^{2}t-A}{t^{\alpha}}\right)\cdot\left(2p^{2}t-A\right)\right]}$$ (B.7)

$$=\frac{1}{t}f^{\prime}\left(\frac{A}{t}\right)\left[A,g\cdot\left(2p^{2}t-A\right)\right]+\text{ h.o.t. }$$

$$=\frac{1}{t}f^{\prime}\left(\frac{A}{t}\right)\left\{g\cdot 4p^{2}t+g^{\prime}\frac{4p^{2}t}{t^{\alpha}}\cdot\left(2p^{2}t-A\right)\right\}+\text{ h.o.t. }$$

$$\approx\frac{1}{t}f^{\prime}\left(\frac{A}{t}\right)\left\{g\left(4p^{2}t-2A\right)+g\cdot 2A+g^{\prime}\cdot\left\{t^{\alpha}-2A\right\}\right\}+\text{ h.o.t. }$$

The RHS of (B.7) has a leading order term

$$\frac{1}{t}f^{\prime}\left(\frac{A}{t}\right)\left(4p^{2}t-2A\right)g\left(-\frac{A+2p^{2}t}{t^{\alpha}}\geqslant 1\right).$$ (B.8)

The other two terms are of order 1 or lower in $t$, since $\frac{A}{t}f^{\prime}\left(\frac{A}{t}\right)$ is bounded.

Commuting

$${\left[\sqrt{F_{A}^{\prime}},f^{\prime}\left(\frac{A}{t}\right)\left(2p^{2}t-A\right)g\left(\frac{2p^{2}t-A}{t^{\alpha}}\geqslant 1\right)\right]}$$

$$=\frac{1}{t}f^{\prime}\left(\frac{A}{t}\right)^{2}\left[\left(4p^{2}t\right)g+\frac{4p^{2}t}{t^{\alpha}}g^{\prime}\left(2p^{2}t-A\right)\right]+h.o.t.$$ (B.9)

Here, $f\left(\frac{A}{t}\right)\equiv\sqrt{\tilde{F}_{A}^{\prime}}.\quad g\left(\frac{2p^{2}t-A}{t^{\alpha}}\right)\equiv F_{A,C}$. Writing $4p^{2}t$ as $4p^{2}t-2A+2A$, then the $g^{\prime}$ factor in B.9 is bounded by

$$2g^{\prime}+2Ag^{\prime}/t^{\alpha}.$$

So, the leading term is

$$\frac{1}{t}f^{\prime}\left(\frac{A}{t}\right)^{2}4p^{2}tg.$$

The leading order remainder term, from B.7 and B.3 is therefore:

$$\frac{1}{M^{2}t^{2}}f^{\prime}(A/t)^{2}4p^{2}tg\left(\frac{2p^{2}t-A}{t^{\alpha}}\geqslant 1\right)=$$

$$\frac{1}{Mt^{2}}\cdot\frac{1}{M^{2}t^{2}}f^{\prime}\left(\frac{A}{Mt}=1\right)^{2}4p^{2}tg.$$ (B. 10)

The leading order term on the RHS is

$$\frac{1}{tM}\cdot\sqrt{\tilde{F}_{A}^{\prime}\,\left(\frac{A}{t\mu}=1\right)}\left(2p^{2}t-A\right)F_{A,C}\sqrt{\tilde{F}_{A}^{\prime}\,\left(\frac{A}{t\mu}=1\right)}.$$ (B.11)

Since we proved before that $\left\langle p^{2}F\left(|p|>t^{1/5}\right)\right\rangle$ is bounded in time, it follows that the leading remainder term, B.10, is bounded by

$$\left\langle\frac{1}{M^{3}}\frac{1}{t^{3}}f^{2}gt^{2/5}\right\rangle\in L^{1}(dt).$$

In particular, we have the bound

$$\int^{\top}\langle B.10\rangle ds\leqslant\int^{\top}\left\langle\frac{1}{M^{3}}f^{\prime 2}\left(\frac{A}{Ms}\sim 1\right)g\right\rangle_{s}\frac{ds}{s^{13/5}}.$$

Multiply by $M,M^{2}$ or $M^{3}$, where

$$M=M_{0}2^{n}\quad n=0,1,2,\ldots$$

and sum over $n,$ we get the following bounds:

$$\int_{1}^{T}\frac{dt}{t}\left(\frac{A}{t}\right)^{l}\left\langle\tilde{F}\left(\frac{A}{t}\geqslant M_{0}\right)\left(2p^{2}-\frac{A}{t}\right)F_{A,C}\tilde{F}_{A}\left(\frac{A}{t}\geqslant M_{0}\right)\right\rangle_{t}+$$

$$+\left\langle\left(\frac{A}{t}\right)^{l+1}\ln\left\langle\frac{A}{t}\right\rangle F_{A,C}+c.c.\right\rangle_{T}-\left\langle\left(\frac{A}{t}\right)^{l}\ln\langle A\rangle F_{A,C}\right\rangle_{1}$$

$$\leq\left\langle\int_{1}^{T}\left\langle\bar{F}_{A}\left(\frac{A}{M_{0}s}\geqslant 1\right)F_{A,C}\right\rangle\frac{ds}{s^{13/5}}+\text{ h.o.t. }\lesssim O(1).\right.$$

$$l=1,2.$$

Remark 7.2.

The leading error term has decay in time that is more than needed for convergence. Therefore, we can redo the estimate, using

$$B_{A,\beta}\equiv F_{A}\left(\frac{A}{Mt^{\beta}}\leqslant\frac{1}{2}\right)F_{A,c}\left(\frac{2p^{2}-\beta A}{t^{\alpha}}\geqslant 1\right)+c\cdot c.$$

The estimates we obtain show that on support $F_{A,C},\left(\frac{A}{t}\right)^{2}$ and $\left(\frac{A}{t^{d}}\right)^{2}$ are uniformly bounded, respectively.

It also follows that

$$\int_{1}^{T}\left(\frac{A}{t^{\beta}}\right)^{2}\left\langle\tilde{F}_{A}\left(\frac{A}{t^{\beta}}\geqslant M_{0}\right)2p^{2}F_{A,C}\tilde{F}_{A}\right\rangle_{t}\frac{dt}{t^{\beta}}<\infty$$

Next, we note that the region $A-2p^{2}t\geqslant t^{\alpha}$, which is also away from the PS, can be treated similarly: In this case we use

$$B_{M}^{-}\equiv F_{A}\left(\frac{A}{tM}\geqslant\frac{1}{2}\right)F_{A,C}\left(\frac{A-2p^{2}t}{t^{\alpha}}\geqslant 1\right)+\text{ c.c. }$$

and notice that $B_{M}^{-}\geqslant 0+$ h.o.t. and the leading order Heisenberg derivative of $B_{M}$ is negative.

$$D_{H_{0}}F_{A,C}=-\frac{\alpha}{t}F_{A,C}^{\prime}\leqslant 0$$

$$D_{H_{0}}F_{A}=\tilde{F}_{A}^{\prime}\left[p^{2}(Mt)^{-1}-\frac{A}{Mt^{2}}\right],$$

which is also negative on the support of $F_{A,C}$. Similarly with $A/tM$ replaced by $A/t^{\beta}M$, $\beta<1$ not too small. The potential commutator is shown, as before, to have decay in the powers of $M$ and $t$.

The symmetrization works as before, as these operators are the same as before up to the sign change. We conclude in particular that

$$\begin{gathered}\int_{1}^{T}\frac{dt}{t}\left\langle\left(\frac{A}{t^{\beta}}\right)^{l}\cdot\tilde{F}_{A}\left(\frac{A}{t^{\beta}}\geqslant M_{0}\right)p^{2}F_{A,C}\left(\frac{\left|2p^{2}t-A\right|}{t^{\alpha}}\geqslant 1\right)\tilde{F}_{A}\left(\frac{A}{t^{\beta}}\geqslant M_{0}\right)\right\rangle_{t}\\ \leqslant\langle\cdots\rangle_{t=1}+O(1).\quad l=0,1,2.\end{gathered}$$

Furthermore, we have

$$\left\langle\left(\frac{A}{t^{\beta}}\right)^{l+1}\ln\frac{\langle A\rangle}{t^{\beta}}F_{A,C}\left(\frac{\left|2p^{2}t-A\right|}{t^{\alpha}}\geqslant 1\right)\right\rangle_{T}$$ (42)

$$\lesssim\langle\cdots\rangle_{1}+o(1).$$ (43)

∎

8 $H^{2}$ - Regularity on the Outgoing Waves

Theorem 8.1.

Let $\psi(t)$ be a solution of the SE as considered above. Assume the initial condition is well localized in space and in $H^{2}\left(R^{3}\right)$.

Then,
a) If $\sup_{t}\langle\phi(t),H(t)\phi(t)\rangle\lesssim O(1)$ where $\phi(t)=F\left(\frac{A}{t^{\alpha}}\geqslant 1\right)\psi(t)$, then

$$\int_{1}^{T}\left\langle\Delta\psi(t),\left[\left\langle\frac{A}{t^{\alpha}}\right\rangle\ln\left\langle A/t^{\alpha}\right\rangle\right]^{-1}\Delta\psi(t)\right\rangle\frac{dt}{t^{\alpha}}\leq O(1)$$ (C.2)

b) If for suitably localized $G(A)$

$$\sup_{t}\langle\nabla\phi(t),G(A)\langle\frac{A}{t^{\alpha}}\rangle\ln{\langle\frac{A}{t^{\alpha}}\rangle}\nabla\phi(t)\rangle\leq O(1)$$ (C.3a)

then

$$\int_{1}^{T}\frac{dt}{t}\left\langle\Delta\psi(t),G(A)F\left(\frac{A}{t^{\alpha}}\geqslant 1\right)\Delta\psi(t)\right\rangle\leq O(1).$$ (C.3b)

Proof.

$$\partial_{t}\langle\psi(t), \left.F_{A,M}\left(\frac{A}{Mt^{\alpha}}\geq 1\right)H(t)F_{A,M}\left(\frac{A}{Mt^{\alpha}}\geq 1\right)\psi(t)\right\rangle$$

$$= \left\langle\frac{1}{Mt^{\alpha}}\tilde{F}_{A,M}^{\prime}\left[2p^{2}-\frac{\alpha A}{t}\right]H(t)F_{A,M}+c.c.\right\rangle$$

$$+\left\langle\psi(t),F_{A,M}\frac{\partial V}{\partial t}F_{A,M}\psi(t)\right\rangle$$ (C.4)

$$+\left\langle\psi(t),\left(i\left[V,F_{A,M}\right]H(t)F_{A,M}+c.c.\right)\psi(t)\right\rangle\equiv I_{1}+I_{2}+I_{3}$$

We proved in the previous section estimates away from the PS.

So, to deal with the first term, we write

$${\left[2p^{2}-\frac{\alpha A}{t}\right]=\left[2p^{2}-\frac{\alpha A}{t}\right]F_{A,C}^{l}\left(\frac{2p^{2}t-A}{t^{\alpha}}\geqslant 1\right)}$$

$$+ {\left[2p^{2}-\frac{\alpha A}{t}\right]F_{A,C}^{L}\left(\frac{A-2p^{2}t}{t^{\alpha}}\geqslant 1\right)}$$

$$+ {\left[2p^{2}-\frac{\alpha A}{t}\right]F_{A,C}^{0}\left(\frac{\left|A-2p^{2}t\right|}{t^{\alpha}}\leq 1\right)}$$ (C.5)

The first and third terms in equation C.5 are positive (for $2p^{2}t\geq t^{\alpha}$). The second term is supported where $\frac{A}{t}\geqslant 2p^{2}+t^{\alpha-1}$ and multiplied by $\tilde{F}_{A,M}^{{}^{\prime}}$ where $\frac{A}{t^{\alpha}}\sim M$.

Therefore, $Mt^{\alpha-1}\geqslant 2p^{2}+t^{\alpha-1}$.

Consequently,

$$\left\langle\tilde{F}_{A,M}^{\prime}\left[2p^{2}-\frac{\alpha A}{t}\right]F_{A,C}^{L}\left(\frac{A-2p^{2}t}{t^{\alpha}}\geqslant 1\right)H(t)F_{A,M}\right)$$

$$\lesssim \left\langle\tilde{F}_{A,M}^{\prime}M^{2}t^{2\alpha-2}F_{A,C}^{L}F_{A,M}\right\rangle$$ (C.6)

$$\text{ since }F_{2}\left(p^{2}\leq Mt^{\alpha-1}\right)H(t)\lesssim F_{2}Mt^{\alpha-1}+F_{1}V$$

and

$$\tilde{F}_{A,M}^{\prime}\cdot F_{1}V\sim \tilde{F}_{A,M}^{\prime}\langle A\rangle^{-2}F_{1}A^{2}V\sim$$

$$\left\langle Mt^{\alpha}\right\rangle^{-2}Mt^{\alpha-1}\left(|x|^{2}V\right)\sim M^{-1}t^{-1-\alpha}.$$

So,

$$I_{1} \geqslant\langle(Mt^{\alpha})^{-1}\tilde{F}^{\prime}_{A,M}2p^{4}F_{A,C}^{l}F_{A,M}\rangle-\left\langle\left(Mt^{\alpha}\right)^{-1}\tilde{F}_{A,M}^{{}^{\prime}}M^{2}t^{2\alpha-2}F_{A,C}^{l}\left(\frac{A-2p^{2}t}{t^{\alpha}}\geqslant 1\right)F_{A,M}\right\rangle$$

$$I_{2} =\left\langle\psi(t),F_{A,M}\frac{\partial V}{\partial t}F_{A,M}\psi(t)\right\rangle$$ (C.7)

$$\approx\left\langle\psi(t),F_{A,M}\langle A\rangle^{-2}\left(A^{2}+1\right)\langle x\rangle^{-m}\left(A^{2}+1\right)\langle A\rangle^{-2}F_{A,M}\psi(t)\right\rangle$$ (C.8)

$$\approx\left\langle\Delta\psi(t),\langle x\rangle^{-\sigma}\tilde{F}_{A,M}t^{-4\alpha}\langle x\rangle^{-m+4+2\sigma}\tilde{F}_{A,M}\langle x\rangle^{-\sigma}\Delta\psi(t)\right\rangle M^{-4}$$

$$\approx O\left(t^{-4\alpha}\right)\|\Delta\psi(0)\|_{2}^{2}M^{-4},$$

provided $m\geqslant 4+2\sigma;\quad\psi(0)\in H^{2}\left(\mathbb{R}^{3}\right)$.

Remark 8.2.

The above estimate with $m\geqslant 4+2\sigma$ is not necessary, and $\geqslant 2+2\sigma$ is sufficient. In the above we demonstrate the maximal decay possible at the $H^{2}$ level.

The control of $I_{3}$ is similar: In this case, due to the $H(t)$ factor, we can only multiply by $\langle A\rangle^{2}$ and divide. By restricting the domain to the PS or where $A\geqslant 2p^{2}t$, two powers of $\langle A\rangle^{-1}$ give a decay of order $t^{-2}$, so, such regions are bounded by

$$\|\psi(0)\|_{H^{2}}^{2}t^{-2}\left(Mt^{2}\right)^{-1}.$$

In the region $A\leq 2p^{2}t$, the $\langle A\rangle^{-2}$ gives

$$\|\psi(0)\|_{H^{2}}^{2}\left(Mt^{\alpha}\right)^{-2}$$

We also used the Commutator Expansion Lemma for

$$\left[V,F_{A,M}\right]\sim x\cdot\nabla VF_{A,M}^{\prime}\left(Mt^{\alpha}\right)^{-1}+\text{ h.o.t. }$$

So,

$$I_{3}\lesssim\left[\left(Mt^{\alpha}\right)^{-3}+t^{-2}\left(Mt^{\alpha}\right)^{-1}\right]\|\psi(0)\|_{H^{2}}^{2}$$ (C.9)

Combining the equations C.4-C.9 we get

$$\left\langle F_{A,M}H(t)F_{A,M}\right\rangle_{t=T}-\left\langle F_{A,M}H(t)F_{A,M}\right\rangle_{t=1}$$

$$=\int_{1}^{T}\frac{dt}{t^{\alpha}}\left\langle\frac{1}{M}\tilde{F}_{A,M}^{\prime}2p^{2}F_{A,C}^{+}p^{2}\tilde{F}_{A,M}^{\prime}\right\rangle+h.o.t.$$

$$+\int_{1}^{T}t^{2\alpha}\frac{dt}{t^{2}t^{\alpha}}\frac{M^{2}}{M}\left\langle\tilde{F}_{A,M}^{\prime}F_{2}\left(p^{2}\leq Mt^{\alpha-1}\right)F_{A,C}^{-}F_{A,M}+c.c.\right\rangle$$

$$+\int_{1}^{T}dt\left\{O\left(t^{-4\alpha}M^{-4}\right)+\left(Mt^{\alpha}\right)^{-3}+O\left(t^{-2}\right)\left(Mt^{\alpha}\right)^{-1}\right\}^{2}\|\Delta\psi(0)\|_{2}^{2}.$$ (C.10)

Here,

$$F_{A,C}^{+}=F_{A,C}\left(\frac{\left|2p^{2}t-A\right|}{t^{\alpha}}\leq 1\right)+F_{A,C}\left(\frac{2p^{2}t-A}{t^{\alpha}}\geqslant 1\right)$$

We use $\tilde{F}_{A,M}^{\prime}$ to denote generic function of the form $\widetilde{F}_{A,M}\left(\frac{A}{Mt^{\alpha}}=1\right)$.

Next, we will use that

$$\sum_{n=0}^{\infty}\left(M_{0}2^{n}\right)^{l}F_{A,M_{0}2^{n}}\sim\left(\frac{A}{t^{\alpha}}\right)^{l}\ln\left(\langle A\rangle/t^{\alpha}\right)F_{A}\left(\frac{A}{t^{\alpha}}\geqslant M_{0}\right)$$

and

$$\sum_{n=0}^{\infty}\left(M_{0}2^{n}\right)^{l}/nF_{A,M_{0}2^{n}}\sim\left(\frac{A}{t^{\alpha}}\right)^{l}F_{A}\left(\frac{A}{t^{\alpha}}\geq M_{0}\right)$$

We sum over $M=M_{0}2^{n},n=0,1,2,\ldots$ both sides of eq. C.10, and sum over $n$ up to $M\sim cT^{1-\alpha}.$

The LHS becomes

$$\left\langle F_{A}\left(cT^{1-\alpha}\geq\frac{A}{t^{\alpha}}\geq M_{0}\right)H(t)F_{A}\left(cT^{1-\alpha}\geq\frac{A}{t^{\alpha}}\geq M_{0}\right)\ln T\right\rangle_{t}-\left\langle F_{A}\left(cT\geq A\geqslant M_{0}\right)\ln T\,H(1)F_{A}\left(cT\geq A\geq M_{0}\right)\right\rangle.$$

The RHS

The first term

$$\int_{1}^{T}dt\left\langle\langle A\rangle^{-1}F_{A}\left(cT^{1-\alpha}\geq\frac{A}{t^{\alpha}}\geqslant M_{0}\right)p^{2}F_{A,C}^{+}p^{2}F_{A}\left(cT^{1-\alpha}\geq\frac{A}{t^{\alpha}}\geqslant M_{0}\right)\right\rangle\sim$$

$$\frac{1}{cT}\int_{1}^{T}\left\langle F_{A}p^{2}F^{+}_{A,C}p^{2}F_{A}\right\rangle dt.$$

The second term

$$\int_{1}^{T}\frac{dt}{t^{2}}\left\langle AF_{A}\left(cT^{1-\alpha}\geq\frac{A}{t^{\alpha}}\geqslant M_{0}\right)F_{2}F_{A,C}^{-}F_{A}\left(cT^{1-\alpha}\geq\frac{A}{t^{\alpha}}\geqslant 1\right)\right\rangle$$

The third term

$$C+\int_{1}^{T}\frac{dt}{t^{2+\alpha}}\|\psi(0)\|_{H^{2}}^{2}$$

The LHS is uniformly bounded by $\ln T$ since $A\leq cT.$ We conclude that a sequence of times going to infinity exists such that the $H^{2}$ norm remains uniformly bounded on the phase space support of the localization of $A.$ This bound is then extended to all times, using as PROB $p^{2}G(A,t)p^{2}$ with $G$ localizing at the relevant domain in phase space. ∎

High Frequency on the Propagation Set

Theorem 8.3.

We have the following estimate on the PS:

$$\left<H(T)\,F_{c}\,H(T)\right>_{T}\lesssim O(1)\,\|\psi_{0}\|_{L^{2}}\,\|\psi_{0}\|_{H^{2}}.$$ (44)

Proof.

To control the $H^{2}$ norm of the solution on the PS we consider the following PROB:

$$B_{HH}=\left<\psi(t),H(t)\,F_{c}\!\left(\frac{|x-2pt|}{R}\leq 1\right)H(t)\psi(t)\right>.$$ (45)

Here $R$ is a large number that may depend on $t$.

Since

$$D_{H}F_{c}=\partial_{t}F_{c}+i[-\Delta,F_{c}]=-\frac{|x-2pt|}{R^{2}}\,\dot{R}\,F_{c}^{\prime},$$ (46)

it is zero if $\dot{R}=0$.

Hence, for $\dot{R}=0$, we get

$$\partial_{t}B_{HH}=\left<\psi(t),\left(\frac{\partial V}{\partial t}F_{c}\,H(t)+H(t)F_{c}\frac{\partial V}{\partial t}\right)\psi(t)\right>$$ (47)

$$\qquad+\left<\psi(t),H(t)\,i[V,F_{c}]\,H(t)\psi(t)\right>.$$

For arbitrarily large $T$, we want to estimate $B_{HH}(T)$ in the PS; for this, we need to take $R$ at least of order $t^{\alpha}$, for some $\alpha<1$.

We have two types of terms to control:

$$\text{a)}\qquad H(t)\,i[V,F_{c}]\,H(t)=p^{2}\,i[V,F_{c}]\,p^{2}+O(V)\,i[V,F_{c}]$$ (48)

$$\text{b)}\qquad\frac{\partial V}{\partial t}\,F_{c}\,H(t)\sim\langle x\rangle^{-\sigma}F_{c}\,p^{2}+O(V)\,\frac{\partial V}{\partial t}\,F_{c}.$$ (49)

We concentrate on the type (a) that is more difficult. The second part (b) is treated similarly.

The key estimate we use is the gain in regularity and smoothness of the product of $V$ and $\widetilde{F}_{c}$:

$$V(x,t)\,F_{c}=V(x,t)\,\langle x\rangle^{\sigma}\langle x\rangle^{-\sigma}e^{i\Delta t}F_{2}(|p|\geq M)F_{1}\!\left(\frac{|x|}{R}\leq 1\right)e^{-i\Delta t}F_{2}(|p|\geq M)$$ (50)

Therefore, we estimate

$$\langle x\rangle^{-\sigma}e^{i\Delta t}F_{2}F_{1}\!\left(\frac{|x|}{R}\leq 1\right).$$

A first observation is

$$\langle x\rangle^{-\sigma}e^{i\Delta t}F_{2}F_{1}\!\left(\frac{|x|}{R}\leq 1\right)=\langle x\rangle^{-\sigma}\langle x-2pt\rangle^{-2}e^{i\Delta t}\widetilde{F}_{2}R^{2}\widetilde{F}_{1}\!\left(\frac{|x|}{R}\leq 1\right)$$ (51)

where

$$\langle x\rangle^{2}F_{1}\!\left(\frac{|x|}{R}\leq 1\right)=R^{2}\widetilde{F}_{1}\!\left(\frac{|x|}{R}\leq 1\right)$$ (52)

and

$$\widetilde{F}_{2}\approx F_{2}+\mathcal{O}(\frac{1}{M}).$$

Proposition 8.4.

For $\sigma>2$,

$$\langle x\rangle^{-\sigma}e^{i\Delta t}F_{2}F_{1}\!\left(\frac{|x|}{R}\leq 1\right)=O(\langle x\rangle^{-m})\langle pt\rangle^{-2}R^{2}O(1)\,\widetilde{F}_{1},\qquad m=\sigma-2.$$ (53)

Proof.

$$\langle x\rangle^{-2}\langle pt\rangle^{-2}\langle pt\rangle^{2}F_{2}\langle x-2pt\rangle^{-2}F_{2}=$$ (54)

$$\langle x\rangle^{-2}\langle pt\rangle^{-2}\Big(1+p^{2}t^{2}-t\,2x\!\cdot\!p-2pt\!\cdot\!x+x^{2}-x^{2}\Big)F_{2}\langle x-2pt\rangle^{-2}F_{2}=$$

$$\langle x\rangle^{-\sigma}\langle pt\rangle^{-2}\widetilde{F}_{2}+\langle x\rangle^{-\sigma}\langle pt\rangle^{-2}[-x^{2}-t\,2p\!\cdot\!x-t\,2x\!\cdot\!p]\widetilde{F}_{2}\langle x-2pt\rangle^{-2}F_{2}$$

$$\langle x\rangle^{-2}\langle pt\rangle^{-2}(x^{2}-2x\!\cdot\!pt)=$$ (55)

$$=\frac{x^{2}}{2\langle x\rangle^{2}}\langle pt\rangle^{-2}-\langle x\rangle^{-2}\partial_{p}^{2}\langle pt\rangle^{-2}-\langle x\rangle^{-2}2x\langle pt\rangle^{-2}\cdot pt$$

$$\langle x\rangle^{-2}\langle pt\rangle^{-2}2pt\!\cdot\!x=\langle x\rangle^{-2}x\langle pt\rangle^{-2}2pt-\langle x\rangle^{-1}O(1)\langle pt\rangle^{-2}2t$$

Multiplying from the right by $\widetilde{F}_{2}$, we get

$$\langle x\rangle^{-2}\langle pt\rangle^{-2}(x^{2}-2x\!\cdot\!pt+2pt\!\cdot\!x)\widetilde{F}_{2}=\langle pt\rangle^{-2}(1+O(1))\widetilde{F}_{2}$$ (56)

$$\qquad+\langle x\rangle^{-2}\langle pt\rangle^{-2}(1+O(1))\frac{p^{2}t^{2}}{1+p^{2}t^{2}}\frac{1}{M^{2}}\widetilde{F}_{2}$$

$$\qquad+\langle x\rangle^{-2}\langle pt\rangle^{-2}2pt\left(\frac{\widetilde{F}^{\prime}}{M}\right).$$

Hence,

$$\langle x\rangle^{-2}\widetilde{F}_{2}\langle x+2pt\rangle^{-2}\widetilde{F}_{2}=(1+O(1))\langle pt\rangle^{-2}\widetilde{F}_{2}+\langle x\rangle^{-2}\langle pt\rangle^{-1}\left(\frac{\widetilde{F}^{\prime}_{2}}{M}\right)\langle x-2pt\rangle^{-2}F_{2}$$ (57)

Therefore,

$$\langle x\rangle^{-2}(1+O(pt)^{-1})\widetilde{F}_{2}\langle x-2pt\rangle^{-2}\widetilde{F}_{2}=(1+O(1))\langle pt\rangle^{-2}\widetilde{F}_{2}$$ (58)

Therefore,

$$\left(M\gg 1,\quad\langle x\rangle^{-2}\frac{O(pt)^{-1}}{M}\ll 1\right)$$

$$\langle x\rangle^{-2}\widetilde{F}_{2}\langle x-2pt\rangle^{-2}\widetilde{F}_{2}\approx\langle pt\rangle^{-2}\widetilde{F}_{2}$$ (59)

Clearly, it also holds that

$$\langle x\rangle^{-2}\widetilde{F}_{2}\langle x-2pt\rangle^{-2}\widetilde{F}_{2}\lesssim\langle x\rangle^{-2}\widetilde{F}_{2}$$ (60)

We also used that

$$i[x,F_{2}]=-\frac{F_{2}^{\prime}}{M},$$

and therefore

$$[\langle x\rangle^{-\sigma},F_{2}]\sim\langle x\rangle^{-\sigma-1}\left(\frac{F_{2}^{\prime}}{M}+O\!\left(\frac{F_{2}^{\prime}}{M^{2}}\right)\right).$$ (61)

For $|p|\leq\varepsilon$ we have:

$$\langle x\rangle^{-\sigma}\langle pt\rangle^{-2}\langle pt\rangle^{2}\widetilde{F}_{2}\langle x-2pt\rangle^{-2}\widetilde{F}_{2}$$

$$\lesssim\langle pt\rangle^{-2}(1+\varepsilon^{2}t^{2})t^{-3/2}\widetilde{F}_{2}e^{-i\Delta t}$$

$$\lesssim C\,\varepsilon^{2}t^{2}t^{-3/2}\langle pt\rangle^{-2}=C\,\varepsilon^{2}t^{1/2}\langle pt\rangle^{-2}.$$

as operators on

$$L^{2}\to L^{2}(\mathbb{R}^{3})$$

In the above we also used the pointwise time decay estimate $\|\langle x\rangle^{-\sigma}e^{i\Delta t}\langle x\rangle^{-\sigma}\|_{2\rightarrow 2}\leq ct^{-3/2}$ in three dimensions. In higher dimensions, we get similar estimates, also by using $\langle x\rangle^{-\sigma}$ instead of $\langle x\rangle^{-2}$, $\sigma>2$.

In $4$-dimensions, we get with $\sigma=2+\delta$ on the left, the bound

$$C\,\varepsilon^{2}t^{2}t^{-2+\delta}\langle pt\rangle^{-2},$$

since

$$\langle x\rangle^{-2-\delta}\phi\in L^{p}\qquad\text{for all }p>1,\phi\in L^{2}.$$

∎

Corollary 8.5.

$$p^{2}\,i[V,F_{2}F_{c}F_{2}]\,p^{2}=p^{2}\widetilde{V}\langle x\rangle^{-\sigma}\langle x\rangle^{-2}F_{2}F_{c}F_{2}p^{2}+c.c.=$$ (62)

$$=p^{2}\langle x\rangle^{-\sigma}\widetilde{V}\langle pt\rangle^{-2}(1+O(1))\widetilde{F}_{2}p^{2}R^{2}+c.c.$$

$$=p^{2}\langle x\rangle^{-\sigma}\widetilde{V}\langle t\rangle^{-2}O(1)\widetilde{F}_{2}R^{2}+c.c.$$

Consequently,

$$\int_{T_{1}}^{T}\left<p^{2}\,i[V,F_{2}F_{c}\widetilde{F}_{2}]\,p^{2}\right>dt\lesssim\langle T\rangle^{-1}\|\widetilde{V}\|_{\infty}\|\psi(0)\|_{H^{2}}\|\psi(0)\|_{L^{2}}R^{2}.$$ (63)

Since

$$R\sim T^{\alpha},$$

we need

$$\alpha<\frac{1}{2}.$$

This gives control over the high frequency part on the PS:

$$F_{2}\,F_{c}\,\widetilde{F}_{2}.$$

To cover the PS over the time interval $[0,T],$ we break the interval into sub-intervals; the first is $(\delta\ll 1)$

$$[\delta T,T]$$

and we choose

$$R=T^{2/5}.$$

Then, we get

$$\left<H(T)\,F_{c}\,H(T)\right>_{T}-\left<H(T_{1})\,F_{c}\,H(T_{1})\right>_{T_{1}=T^{4/5}}=$$ (64)

$$=\int_{T_{1}}^{T}\partial_{s}\left<\psi(s),H(s)\,F_{c,R}\,H(s)\psi(s)\right>\lesssim T_{1}^{-1}R^{2}\|\psi\|_{H^{2}}\|\psi(0)\|_{L^{2}}$$

$$\lesssim O(1)\,\|\psi_{0}\|_{L^{2}}\,\|\psi_{0}\|_{H^{2}}.$$

Similarly,

$$\left<H(T_{1})\,F_{c}\,H(T_{1})\right>-\left<H(T_{2})\,F_{c}\,H(T_{2})\right>_{T_{2}=T_{1}^{4/5}}$$ (65)

$$\lesssim O(1)\,\|\psi_{0}\|_{L^{2}}\,\|\psi_{0}\|_{H^{2}}.$$

We have

$$T_{n}=T_{n-1}^{4/5}=(T_{n-2})^{\frac{4}{5}\cdot\frac{4}{5}}=\cdots=T_{1}^{(4/5)^{n-1}}=T^{(4/5)^{n}}$$ (66)

Hence, for $T_{n}\sim 1$, we need

$$\left(\frac{4}{5}\right)^{n}\ln T=O(1),$$ (67)

or

$$\ln\ln T\sim n(\ln 5/4).$$ (68)

In order to remove the growth of $\ln\ln T$, it is sufficient to replace the factor $2/5$ by $2/5-0$ in the definition of $F_{c}.$ Then each interval contribution will be suppressed by a factor $T_{n}^{\epsilon},$ so the sum over all contributions is finite, since $T_{n}$ converges to 1 super exponentially fast, and the number of terms is $c\ln\ln T.$ The resulting estimate of the norm $H^{2}$ completes the part of the phase-space left from the previous section. ∎

9 Final Comments

There are problems where higher regularity cannot be achieved without proving Local decay estimates and other $L^{p}$ decay estimates. The approach introduced in this work does not use such methods.

Some simple examples of applications include potentials that correspond to modulated soliton solutions of NLS type equations. The standard example is the soliton part of the dynamics in the modulation equations that couple the changing Soliton to the Radiation. In this case the potential is essentially a smooth localized function with width and amplitude varying in an unknown way in time.

The general case of time periodic, quasi periodic and some almost periodic space-localized potentials is covered.

Also included Saturated Nonlinear Equations, where the nonlinear term is of the general form

$$I(x,t)\frac{|\psi|^{2}}{1+V(x,t)|\psi|^{2}}$$

with smooth localized function $I$, and smooth localized plus a constant function $V.$

More intricate is the case of semilinear equations. First, it should be observed that the Bootstrap step on the time interval $t-\epsilon,t]$ in the proof of Theorem 2.2, can be extended to the nonlinear case. Observe that the a nonlinear bootstrap estimate is sufficient,by simply choosing the size of the time interval, $\epsilon.$ small compared with a power of the size of the $H^{2}$ norm of the Initial Data. This shows that the solution is in $H^{2}$ locally in space and uniform in time. Localization of the interaction term for the next steps can be obtained from the assumption of spherical symmetry. Acknowledgment

A.S. is supported in part by NSF grant DMS-2205931

soffer@math.rutgers.edu