arXiv is now an independent nonprofit! Learn more
License: arXiv.org perpetual non-exclusive license
arXiv:2608.20072v1 [math.AP] 20 Aug 2026

Scattering for the focusing H1/2H^{1/2}-critical nonlinear Schrödinger equation with large data

Qiuye Jia Address: Australian National University Email address: Qiuye.Jia@anu.edu.au
Abstract.

In this article we prove that the solution to the focusing H1/2H^{1/2}-critical nonlinear Schrödinger equation in dimension d5d\geq 5 scatters outside finitely many arbitrarily small spacetime cones. We will discuss the relationship between Theorem 1.5 and the resolution of solitons. There are two key ingredients in the proof: the localized below-ground-state scattering in Theorem 1.1, and the characterization in Theorem 1.2 of the asymptotic behaviour of the residual part (i.e., after subtracting the scattering part).

As a byproduct, we also establish an upgrading mechanism: if the solution is below the ground state (resp. small) along a time sequence in a spacetime cone, then it is also below ground state (resp. small) uniformly for large times in a shrunk spacetime cone. We expect this to be useful in future works using concentration compactness and our spacetime cone approach.

Keywords: 
Focusing nonlinear Schrödinger equation, Compact attractor, Localized scattering, Concentration compactness, Soliton resolution
1991 Mathematics Subject Classification
35Q55, 35B40

1. Introduction

1.1. The setup and main results

In this article, we study the forward-global solutions with uniformly bounded H1H^{1} norm to the focusing nonlinear Schrödinger equation

(1.1) (it+Δ)u=|u|4/(d1)ufor t0,xd,u0=u|t=0H1(d),d5.\displaystyle\begin{split}&(i\partial_{t}+\Delta)u=-|u|^{4/(d-1)}u\qquad\text{for }t\geq 0,\ x\in\mathbb{R}^{d},\\ &u_{0}=u|_{t=0}\in H^{1}(\mathbb{R}^{d}),\quad d\geq 5.\end{split}

The critical scaling level of this equation is H1/2(d)H^{1/2}(\mathbb{R}^{d}) and the actual scaling critical quantity we will use is mass times energy, as defined below. We adopt the convention Δ=i=1dxi2\Delta=\sum_{i=1}^{d}\partial_{x_{i}}^{2}. Throughout the paper, we assume the solution is uniformly bounded in H1H^{1}:

(1.2) supt0u(t)H1<.\displaystyle\sup_{t\geq 0}\|u(t)\|_{H^{1}}<\infty.

For fH1(d)f\in H^{1}(\mathbb{R}^{d}), we define its mass and energy (associated with our particular NLS) by

(1.3) M(f):=d|f|2dx,P(f):=df¯fdx,E(f):=12d|f|2dxd12(d+1)d|f|2(d+1)/(d1)dx.\displaystyle\begin{split}M(f)&:=\int_{\mathbb{R}^{d}}|f|^{2}\,dx,\quad P(f):=\Im\int_{\mathbb{R}^{d}}\overline{f}\nabla f\,dx,\\ E(f)&:=\frac{1}{2}\int_{\mathbb{R}^{d}}|\nabla f|^{2}\,dx-\frac{d-1}{2(d+1)}\int_{\mathbb{R}^{d}}|f|^{2(d+1)/(d-1)}\,dx.\end{split}

The aim of this paper is to show that with arbitrarily large (but fixed) mass or energy, the solution to (1.1) still ‘almost scatters’. By ‘scatter’, as usual, we mean that u(t)u(t) converges to a solution to the free linear Schrödinger equation. More concretely, it is known that the weak limit eitΔue^{-it\Delta}u as tt\to\infty exists and we will denote it by u+u_{+}, which is the ‘radiation data’. See also Theorem 2.7 below for how it lies in a decomposition of uu. Then

(1.4) v(t)=eitΔu+,r(t)=u(t)v(t)\displaystyle v(t)=e^{it\Delta}u_{+},\quad r(t)=u(t)-v(t)

are the free part and the soliton part (or residual part, in the context below) of our solution u(t)u(t). Then ‘u(t)u(t) scatters’ means r(t)0r(t)\to 0.

However, there is a well-known obstruction for this statement to be true: the ‘soliton’ mentioned above in the name of r(t)r(t). Let QQ be the unique positive radial decaying solution of (see e.g. [19, 26])

(1.5) ΔQ+Q|Q|4/(d1)Q=0,\displaystyle-\Delta Q+Q-|Q|^{4/(d-1)}Q=0,

which is called the ‘ground state soliton’. Then eitQ(x)e^{it}Q(x) solves (1.1) as well but it does not scatter. There is a long tradition of avoiding this obstruction and obtaining scattering by imposing a threshold condition stating that the solution is ‘smaller’ than QQ. The threshold governing scattering versus soliton behaviour of solutions to (1.1) is:

(1.6) Θ\displaystyle\Theta :=M(Q)E(Q).\displaystyle:=M(Q)E(Q).

Indeed, this serves as a threshold to give scattering versus soliton dichotomy. See Theorem 4.1 and more discussion on the more complicated situation when it is on or above the threshold in Section 1.2.

The aim of this article is to show that even when the solution is above this threshold in the sense that M(u)E(u)>M(Q)E(Q)M(u)E(u)>M(Q)E(Q), we can still say that uu scatters except on finitely many directions of velocities of soliton components. To this end, it is more convenient to make statements using spacetime cones. For Ωd\Omega\subset\mathbb{R}^{d}, we call a region like

(1.7) C(Ω)={(t,x)×d:tT,x/tΩ},T1C(\Omega)=\{(t,x)\in\mathbb{R}\times\mathbb{R}^{d}:t\geq T,\ x/t\in\Omega\},\quad T\gg 1

a spacetime cone; thus a cutoff χ(x/t)\chi(x/t) localizes to a fixed region in the velocity variable x/tx/t.

Our first result, Theorem 1.1, is the below ground-state scattering localized to a spacetime cone. See Section 1.2 for a discussion of the global version of this result proved previously.

Theorem 1.1.

Let d5d\geq 5, and let uC(0,H1(d))u\in C(\mathbb{R}_{\geq 0};H^{1}(\mathbb{R}^{d})) be a forward-global solution of (1.1) satisfying (1.2). Let u+H1(d)u_{+}\in H^{1}(\mathbb{R}^{d}) be the radiation state of uu. Let χ1,χ2Cc(d)\chi_{1},\chi_{2}\in C_{c}^{\infty}(\mathbb{R}^{d}) be real-valued and satisfy

(1.8) 0χ1,χ21,χ2\displaystyle 0\leq\chi_{1},\chi_{2}\leq 1,\qquad\chi_{2} =1on an open neighborhood of suppχ1.\displaystyle=1\quad\text{on an open neighborhood of }\mathrm{supp}\chi_{1}.

Fix T01T_{0}\geq 1 and δ(0,1)\delta\in(0,1), and assume that χ2(x/t)u(t)\chi_{2}(x/t)u(t) is below threshold in the sense that for every tT0t\geq T_{0}, we have:

(1.9) M(χ2(x/t)u(t))E(χ2(x/t)u(t))\displaystyle M(\chi_{2}(x/t)u(t))E(\chi_{2}(x/t)u(t)) (1δ)M(Q)E(Q),\displaystyle\leq(1-\delta)M(Q)E(Q),
(1.10) χ2(x/t)u(t)L2(χ2(x/t)u(t))L2\displaystyle\|\chi_{2}(x/t)u(t)\|_{L^{2}}\|\nabla(\chi_{2}(x/t)u(t))\|_{L^{2}} <QL2QL2.\displaystyle<\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}}.

Then it scatters in a slightly shrunk cone:

(1.11) limtχ1(x/t)(u(t)eitΔu+)H1\displaystyle\lim_{t\to\infty}\|\chi_{1}(x/t)(u(t)-e^{it\Delta}u_{+})\|_{H^{1}} =0.\displaystyle=0.

With r(t)r(t) as in (1.4), we will use its mass density and current as densities to define measures to identify the asymptotic behaviour of the solution. Concretely, we let X=d{}X=\mathbb{R}^{d}\cup\{\infty\} be the one-point compactification 11 1 Introducing XX is completely for the convenience of writing. One can use d\mathbb{R}^{d} and consider the case with velocity going to infinity by a separate discussion. of d\mathbb{R}^{d}, which one should think of as parametrizing the asymptotic directions. Then we define νt\nu_{t}, JtJ_{t} by:

(1.12) Xhdνt:=dh(x/t)|r(t,x)|2dx,XhdJt:=dh(x/t)(r¯r)(t,x)dx.\displaystyle\int_{X}h\,d\nu_{t}:=\int_{\mathbb{R}^{d}}h(x/t)|r(t,x)|^{2}\,dx,\qquad\int_{X}h\,dJ_{t}:=\int_{\mathbb{R}^{d}}h(x/t)\Im(\overline{r}\nabla r)(t,x)\,dx.

Our Theorem 1.2 below shows that they indeed converge in the weak sense as tt\to\infty and is the technically key step for us to identify those non-scattering directions. The intuition that we should expect this can be seen from ‘assuming’ a soliton resolution type expansion. See the discussion at the beginning of Section 5.

Theorem 1.2.

Let d5d\geq 5, let MM, EE, and PP be as in (1.3), let QQ be as in (1.5), and let Θ\Theta be as in (1.6). Let uC(0,H1(d))u\in C(\mathbb{R}_{\geq 0};H^{1}(\mathbb{R}^{d})) solve

(it+Δ)u=|u|4/(d1)u\displaystyle(i\partial_{t}+\Delta)u=-|u|^{4/(d-1)}u

and satisfy (1.2). Let u+u_{+} be its radiation state, let r(t)r(t) be as in (1.4), put m:=M(u)M(u+)m_{\sharp}:=M(u)-M(u_{+}), and assume m>0m_{\sharp}>0. Define νt\nu_{t} and JtJ_{t} by (1.12). Then there are L1L\geq 1, distinct y1,,yLdy_{1},\ldots,y_{L}\in\mathbb{R}^{d}, and numbers a>0a_{\ell}>0 such that

(1.13) (νt,Jt)(μ,y2μ),μ==1Laδy as t.\displaystyle(\nu_{t},J_{t})\rightharpoonup(\mu_{\infty},\frac{y}{2}\mu_{\infty}),\quad\mu_{\infty}=\sum_{\ell=1}^{L}a_{\ell}\delta_{y_{\ell}}\quad\text{ as }t\to\infty.

Next we discuss the estimates upgrading a condition along a sequence of times to a uniform estimate. The substantive new step is the threshold-free Theorem 1.3.

Theorem 1.3.

Let d5d\geq 5, and let uC(0,H1(d))u\in C(\mathbb{R}_{\geq 0};H^{1}(\mathbb{R}^{d})) solve

(it+Δ)u=|u|4/(d1)u\displaystyle(i\partial_{t}+\Delta)u=-|u|^{4/(d-1)}u

and satisfy (1.2). Let u+u_{+} be its radiation state and let r(t)r(t) be as in (1.4). Suppose that χ0,χ1Cc(d)\chi_{0},\chi_{1}\in C_{c}^{\infty}(\mathbb{R}^{d}) satisfy 0χi10\leq\chi_{i}\leq 1 for i=0,1i=0,1 and

supp(χ0){y:χ1(y)=1}.\displaystyle\mathrm{supp}(\chi_{0})\Subset\{y:\chi_{1}(y)=1\}.

Let TnT_{n}\to\infty, and assume

(1.14) χ1(x/Tn)r(Tn)H1(d)0.\displaystyle\|\chi_{1}(x/T_{n})r(T_{n})\|_{H^{1}(\mathbb{R}^{d})}\to 0.

Then

(1.15) limnsuptTnχ0(x/t)r(t)H1(d)=0.\displaystyle\lim_{n\to\infty}\sup_{t\geq T_{n}}\|\chi_{0}(x/t)r(t)\|_{H^{1}(\mathbb{R}^{d})}=0.

The reason that we want to have such a result upgrading the smallness along a time sequence to the uniform smallness is that the decompositions given by Theorem 2.7 need not carry persistent profile labels. The profiles wjw_{j} and the centers xj,nx_{j,n} in (2.26) are extracted along a subsequence, and Theorem 2.7 does not prevent infinitely many collisions with various possible outcomes. On the other hand, our Theorem 1.3 does not trace an individual profile through those interactions. It instead uses the convergence of the measures (1.12) given by Theorem 1.2.

We will use Theorems 1.2 and 1.3 to prove Theorem 1.4 below. See Section 8.2 for details. It upgrades a threshold bound available only along a sequence of times to one valid for all large times. We will use the upgrading mechanism (more concretely, Theorem 1.4) in the proof of our main result, i.e., Theorem 1.5.

Theorem 1.4.

Let d5d\geq 5, and let uC(0,H1(d))u\in C(\mathbb{R}_{\geq 0};H^{1}(\mathbb{R}^{d})) solve

(it+Δ)u=|u|4/(d1)u\displaystyle(i\partial_{t}+\Delta)u=-|u|^{4/(d-1)}u

satisfying (1.2). Let u+H1(d)u_{+}\in H^{1}(\mathbb{R}^{d}) be the radiation state of uu. Fix χ0,χ1,χ2Cc(d)\chi_{0},\chi_{1},\chi_{2}\in C_{c}^{\infty}(\mathbb{R}^{d}) satisfying 0χi10\leq\chi_{i}\leq 1 for i=0,1,2i=0,1,2 and

(1.16) supp(χ0){y:χ1(y)=1},supp(χ1){y:χ2(y)=1}.\displaystyle\mathrm{supp}(\chi_{0})\Subset\{y:\chi_{1}(y)=1\},\qquad\mathrm{supp}(\chi_{1})\Subset\{y:\chi_{2}(y)=1\}.

Suppose that δ(0,1)\delta\in(0,1), tnt_{n}\to\infty, and

(1.17) M(χ2(x/tn)u(tn))E(χ2(x/tn)u(tn))(1δ)M(Q)E(Q)\displaystyle M(\chi_{2}(x/t_{n})u(t_{n}))E(\chi_{2}(x/t_{n})u(t_{n}))\leq(1-\delta)M(Q)E(Q)

for every nn. Then there is T1T\geq 1 such that, for every tTt\geq T,

(1.18) M(χ0(x/t)u(t))E(χ0(x/t)u(t))(1δ/2)M(Q)E(Q).\displaystyle M(\chi_{0}(x/t)u(t))E(\chi_{0}(x/t)u(t))\leq(1-\delta/2)M(Q)E(Q).

After increasing TT,

(1.19) χ0(x/t)u(t)L2(χ0(x/t)u(t))L2<QL2QL2.\displaystyle\|\chi_{0}(x/t)u(t)\|_{L^{2}}\|\nabla(\chi_{0}(x/t)u(t))\|_{L^{2}}<\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}}.

We now state our main result, which says that the solution to (1.1) scatters except near finitely many directions.

Theorem 1.5.

Let d5d\geq 5, and let uC(0,H1(d))u\in C(\mathbb{R}_{\geq 0};H^{1}(\mathbb{R}^{d})) solve

(it+Δ)u=|u|4/(d1)u\displaystyle(i\partial_{t}+\Delta)u=-|u|^{4/(d-1)}u

and satisfy (1.2). Let u+H1(d)u_{+}\in H^{1}(\mathbb{R}^{d}) be the radiation state of uu. Then there is an integer N0N\geq 0 and distinct v1,,vNdv_{1},\ldots,v_{N}\in\mathbb{R}^{d}. Let V={v1,,vN}V=\{v_{1},\ldots,v_{N}\}. For every ρ>0\rho>0, there is χρC(d)\chi_{\rho}\in C^{\infty}(\mathbb{R}^{d}), with 0χρ10\leq\chi_{\rho}\leq 1 and χρL(d)\nabla\chi_{\rho}\in L^{\infty}(\mathbb{R}^{d}), satisfying

χρ(y)\displaystyle\chi_{\rho}(y) =0\displaystyle=0 if dist(y,V)ρ,\displaystyle\text{if }\mathrm{dist}(y,V)\leq\rho,
χρ(y)\displaystyle\chi_{\rho}(y) =1\displaystyle=1 if dist(y,V)2ρ,\displaystyle\text{if }\mathrm{dist}(y,V)\geq 2\rho,

and

(1.20) limtχρ(x/t)(u(t)eitΔu+)H1(d)=0.\displaystyle\lim_{t\to\infty}\|\chi_{\rho}(x/t)(u(t)-e^{it\Delta}u_{+})\|_{H^{1}(\mathbb{R}^{d})}=0.

Thus χρ1\chi_{\rho}\equiv 1 when N=0N=0, and (1.20) is then precisely forward scattering. More concretely, we have

(1.21) N\displaystyle N u(0)L2supt0u(t)L22M(Q)E(Q)=du(0)L2supt0u(t)L2QL2QL2.\displaystyle\leq\frac{\|u(0)\|_{L^{2}}\sup_{t\geq 0}\|\nabla u(t)\|_{L^{2}}}{\sqrt{2M(Q)E(Q)}}=\frac{\sqrt{d}\,\|u(0)\|_{L^{2}}\sup_{t\geq 0}\|\nabla u(t)\|_{L^{2}}}{\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}}}.

Thus the exceptional spacetime cones about the finitely many velocity directions may be chosen with arbitrarily small aperture.

Remark 1.6.

One can view this as a soliton resolution type result (characterizing all possibilities of non-scattering solitons remains a very difficult step) on the velocity level. This groups all non-scattering solutions with the same asymptotic velocity into a single cluster, and we will need to characterize the sublinear dynamics of these non-scattering solutions in future works.

Our expectation is that the method can be extended to the intercritical range (i.e., between mass-critical and energy-critical). On the other hand, the extension to the mass-critical nonlinearity (like |u|4/du|u|^{4/d}u) and the energy critical case (like |u|4/(d2)u|u|^{4/(d-2)}u) would require some essentially new ingredient, in particular the compact attractor result like Theorem 2.7 in these settings. Also, there is an extra scaling parameter to deal with in these settings.

Proof.

Let r(t)r(t) be as in (1.4), put m=M(u)M(u+)m_{\sharp}=M(u)-M(u_{+}), and let Θ\Theta be as in (1.6). By Proposition 2.10, M(r(t))m0M(r(t))\to m_{\sharp}\geq 0. If m=0m_{\sharp}=0, then, for every R>0R>0,

(1.22) lim suptr(t)H1lim suptP>Rr(t)H1,\displaystyle\limsup_{t\to\infty}\|r(t)\|_{H^{1}}\leq\limsup_{t\to\infty}\|P_{>R}r(t)\|_{H^{1}},

because PRr(t)H1(1+R2)1/2r(t)L20\|P_{\leq R}r(t)\|_{H^{1}}\leq(1+R^{2})^{1/2}\|r(t)\|_{L^{2}}\to 0. Using Proposition 2.8, we know r(t)0r(t)\to 0 in H1H^{1}. Thus the solution itself scatters and we may take N=0N=0 and χρ1\chi_{\rho}\equiv 1.

Now we consider the case m>0m_{\sharp}>0. By Theorem 1.2, we have

(1.23) νtμ=j=1Najδvj,aj>0,vjd,\displaystyle\nu_{t}\rightharpoonup\mu_{\infty}=\sum_{j=1}^{N}a_{j}\delta_{v_{j}},\qquad a_{j}>0,\quad v_{j}\in\mathbb{R}^{d},

where the vjv_{j} are distinct. Fix χ\chi as in Proposition 2.11 and define 𝒜\mathcal{A} by (2.38). We first prove (1.24):

(1.24) 𝒜(vj)Θ(1jN).\displaystyle\mathcal{A}(v_{j})\geq\Theta\qquad(1\leq j\leq N).

If instead 𝒜(vj)<Θ\mathcal{A}(v_{j})<\Theta, the definition (2.38) gives ρ0>0\rho_{0}>0, δ(0,1)\delta\in(0,1), and tnt_{n}\to\infty such that, with

ψ3(y)=χ(yvjρ0),\displaystyle\psi_{3}(y)=\chi\Big(\frac{y-v_{j}}{\rho_{0}}\Big),

one has

M(ψ3(x/tn)u(tn))E(ψ3(x/tn)u(tn))(1δ)Θ.\displaystyle M(\psi_{3}(x/t_{n})u(t_{n}))E(\psi_{3}(x/t_{n})u(t_{n}))\leq(1-\delta)\Theta.

Choose ψ0,ψ1,ψ2Cc(d)\psi_{0},\psi_{1},\psi_{2}\in C_{c}^{\infty}(\mathbb{R}^{d}) satisfying 0ψi10\leq\psi_{i}\leq 1 for i=0,1,2i=0,1,2, ψ0(vj)=1\psi_{0}(v_{j})=1, and

suppψ0{ψ1=1},suppψ1{ψ2=1},suppψ2{ψ3=1}.\displaystyle\mathrm{supp}\psi_{0}\Subset\{\psi_{1}=1\},\qquad\mathrm{supp}\psi_{1}\Subset\{\psi_{2}=1\},\qquad\mathrm{supp}\psi_{2}\Subset\{\psi_{3}=1\}.

Then Theorem 1.4, applied to ψ1,ψ2,ψ3\psi_{1},\psi_{2},\psi_{3} (being χ0,χ1,χ2\chi_{0},\chi_{1},\chi_{2} there respectively), gives the two eventual below-threshold hypotheses of Theorem 1.1 for ψ0,ψ1\psi_{0},\psi_{1}. Hence

ψ0(x/t)r(t)H10.\displaystyle\|\psi_{0}(x/t)r(t)\|_{H^{1}}\to 0.

On the other hand, (1.23) gives

ψ0(x/t)r(t)L22=ψ02dνtψ02dμaj>0,\displaystyle\|\psi_{0}(x/t)r(t)\|_{L^{2}}^{2}=\int\psi_{0}^{2}\,d\nu_{t}\to\int\psi_{0}^{2}\,d\mu_{\infty}\geq a_{j}>0,

a contradiction. This proves (1.24).

Taking c0=Θc_{0}=\Theta in Proposition 2.11, we obtain the first bound in (1.21). Then (3.1) gives the last equality.

Set V={v1,,vN}V=\{v_{1},\ldots,v_{N}\}. For a prescribed ρ>0\rho>0, choose χρ\chi_{\rho} as in Theorem 1.5; then dist(suppχρ,V)ρ\mathrm{dist}(\mathrm{supp}\chi_{\rho},V)\geq\rho. Since the full trajectory in (1.23) converges to a measure supported on VV, Proposition 6.3, applied with An=nA_{n}=n, Bn=n+1B_{n}=n+1, and K=VK=V, yields

supntn+1χρ(x/t)r(t)H10.\displaystyle\sup_{n\leq t\leq n+1}\|\chi_{\rho}(x/t)r(t)\|_{H^{1}}\to 0.

This is (1.20). ∎

1.2. Links to the existing literature

The threshold in (1.9)–(1.10) is the ground-state threshold of the intercritical focusing NLS. See Weinstein’s work [26] on its relationship with the sharp Gagliardo–Nirenberg inequality and the characterization and uniqueness given by Kwong [19].

Such a threshold for scattering arose in the pioneering work treating the energy-critical NLS by Kenig and Merle [17], where the concentration-compactness and rigidity method was introduced to prove below-ground-state scattering in the radial case. The below ground-state scattering in the non-radial case for the energy-critical NLS in d5d\geq 5 is proven by Killip and Visan [18]. In the intercritical range (i.e., mass supercritical and energy subcritical) the sharp form of the threshold is due to Holmer and Roudenko [16], for radial data and for the three-dimensional cubic equation, which is (1.1) with d=3d=3. Afterwards, Duyckaerts, Holmer, and Roudenko [7] removed the radial assumption in that dimension.

Fang, Xie, and Cazenave [10] prove below ground-state scattering in general dimension and for every intercritical power, and Guevara [14] also proved below ground-state scattering in that generality. We use this below ground-state scattering theory in Proposition 3.2. See also [1] for more discussions.

A second route on this question uses the Morawetz estimate. Dodson and Murphy proved below ground-state scattering first for the three-dimensional cubic equation [5] with radial data, and then, through an interaction Morawetz estimate, for non-radial data in every dimension d3d\geq 3 for the H˙1/2\dot{H}^{1/2}-critical equation [6], which has the same type of nonlinearity as (1.1); we use this result as Theorem 4.1. The two-dimensional radial case is treated in [2]. Any of [10, 14, 6] would serve equally well as the input to Section 4.

Our Theorem 1.1 is the analogue of the global below ground-state scattering results, localized to a spacetime cone as in (1.7). On the other hand, in the global results [7, 10, 14, 6], the threshold hypothesis is imposed on the initial data globally and yields scattering of the solution as a whole. Our Theorem 1.1 aims to allow large initial data and include soliton components while deriving scattering in the region that is away from these soliton components. So Theorem 1.1 fits in a more general picture described by the soliton resolution conjecture. In our setting, no global threshold assumption is made: the solution is arbitrary, subject only to a uniform H1H^{1} bound. Finally, the hypothesis is asked for at all large times, whereas concentration compactness produces information only along sequences of times. Theorem 1.4 closes precisely that gap, which is why it is isolated as a separate statement. After obtaining the localized scattering in Theorem 1.1, there can only be finitely many directions on which scattering fails, which is our Theorem 1.5.

On the other hand, starting with the work of Duyckaerts and Merle [8], there have been many successes in characterizing the dynamics of NLS solutions that are on or above the threshold determined by the ground-state soliton. The works that are most closely related to our current setting are the works of Duyckaerts and Roudenko [9] and the work of Nakanishi and Schlag [20].

The limitation d5d\geq 5 that we have in our main results comes from Theorem 2.7, which is [23, Theorem 1.28]. It is used to give a profile decomposition for us. See also [21][24][25].

We should also mention that the spacetime cone approach is inspired by the work of Gell-Redman, Gomes and Hassell [12, 11], and the long tradition of microlocal analysis. Although these works are not directly used on a technical level, they have provided conceptual guidance to the author.

1.3. The structure of the paper

The structure of the paper is as follows. We recall some facts from the existing literature in Sections 23 on linear and nonlinear solutions (mostly the compact attractor and the ground state soliton). Then we prove Theorem 1.1 in Section 4. In Sections 56 we study the asymptotic behaviour of the measure using some parts of our solution as the density. This will be the key step to identify those directions on which we could potentially have solitons. Afterwards, Section 7 proves Theorem 1.2. Finally, we will prove our byproduct: the estimates upgrading the smallness or the below threshold property along a time sequence to a uniform property in Sections 8.1 and 8.2.

Acknowledgements

The project grew out of an ongoing project with Andrew Hassell. The author is very grateful to Andrew Hassell for many helpful conversations and sharing the idea of using spacetime cones. The author is supported by the Australian Research Council through grant FL220100072.

2. Preliminaries: Radiation and compact attractor

2.1. Free radiation

We begin by recording the standard estimates for the free Schrödinger equation that are used throughout the paper. The estimates in Proposition 2.1 are standard; see [22, Theorem 2.3] or [4, Chapter 2].

Proposition 2.1 (Dispersive and Strichartz estimates).

Let d2d\geq 2. For every t0t\neq 0, we have

(2.1) eitΔϕL(d)(4π|t|)d/2ϕL1(d)(ϕL1(d)).\displaystyle\|e^{it\Delta}\phi\|_{L^{\infty}(\mathbb{R}^{d})}\leq(4\pi|t|)^{-d/2}\|\phi\|_{L^{1}(\mathbb{R}^{d})}\qquad(\phi\in L^{1}(\mathbb{R}^{d})).

Then interpolation and the unitarity of eitΔe^{it\Delta} on L2(d)L^{2}(\mathbb{R}^{d}) give:

(2.2) eitΔϕLp(d)Cd,p|t|d(121p)ϕLp(d)(ϕLp(d)),\displaystyle\|e^{it\Delta}\phi\|_{L^{p}(\mathbb{R}^{d})}\leq C_{d,p}|t|^{-d(\frac{1}{2}-\frac{1}{p})}\|\phi\|_{L^{p^{\prime}}(\mathbb{R}^{d})}\qquad(\phi\in L^{p^{\prime}}(\mathbb{R}^{d})),

where 2p2\leq p\leq\infty and 1p+1p=1\frac{1}{p^{\prime}}+\frac{1}{p}=1. We say that a pair (q,r)(q,r) is admissible if

(2.3) 2q+dr=d2,2q,r,\displaystyle\frac{2}{q}+\frac{d}{r}=\frac{d}{2},\quad 2\leq q,r\leq\infty,

excluding (q,r)=(2,)(q,r)=(2,\infty) when d=2d=2. Then, for all admissible (q,r)(q,r) and (q~,r~)(\tilde{q},\tilde{r}), every interval II\subset\mathbb{R}, and every t0It_{0}\in I,

(2.4) eitΔfLtqLxr(I×d)\displaystyle\|e^{it\Delta}f\|_{L_{t}^{q}L_{x}^{r}(I\times\mathbb{R}^{d})} Cd,qfL2(d),\displaystyle\leq C_{d,q}\|f\|_{L^{2}(\mathbb{R}^{d})},
(2.5) t0tei(ts)ΔF(s)𝑑sLtqLxr(I×d)\displaystyle\Big\|\int_{t_{0}}^{t}e^{i(t-s)\Delta}F(s)\,ds\Big\|_{L_{t}^{q}L_{x}^{r}(I\times\mathbb{R}^{d})} Cd,q,q~FLtq~Lxr~(I×d),\displaystyle\leq C_{d,q,\tilde{q}}\|F\|_{L_{t}^{\tilde{q}^{\prime}}L_{x}^{\tilde{r}^{\prime}}(I\times\mathbb{R}^{d})},

with constants independent of II and t0t_{0}.

We next record some other properties of the linear solutions.

Lemma 2.2.

Let d5d\geq 5. There is a constant CdC_{d} such that

(2.6) eitΔfL2(d+1)/(d1)(t×xd)CdfH1(d)\displaystyle\|e^{it\Delta}f\|_{L^{2(d+1)/(d-1)}(\mathbb{R}_{t}\times\mathbb{R}_{x}^{d})}\leq C_{d}\|f\|_{H^{1}(\mathbb{R}^{d})}

for every fH1(d)f\in H^{1}(\mathbb{R}^{d}). In particular,

(2.7) eitΔfL2(d+1)/(d1)2(d+1)/(d1)𝑑t<\displaystyle\int_{\mathbb{R}}\|e^{it\Delta}f\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}\,dt<\infty

and

(2.8) limTTeitΔfL2(d+1)/(d1)2(d+1)/(d1)dtt=0.\displaystyle\lim_{T\to\infty}\int_{T}^{\infty}\|e^{it\Delta}f\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}\frac{dt}{t}=0.
Proof.

Set s=1d+1s_{*}=\frac{1}{d+1}, r=2d(d+1)d2d+2r_{*}=\frac{2d(d+1)}{d^{2}-d+2}. Then (2(d+1)/(d1),r)(2(d+1)/(d-1),r_{*}) is admissible in the sense of (2.3). Apply (2.4) to |D|sf|D|^{s_{*}}f to obtain

|D|seitΔfLt2(d+1)/(d1)LxrCdfH˙s.\displaystyle\||D|^{s_{*}}e^{it\Delta}f\|_{L_{t}^{2(d+1)/(d-1)}L_{x}^{r_{*}}}\leq C_{d}\|f\|_{\dot{H}^{s_{*}}}.

Since

d12(d+1)=1rsd,\displaystyle\frac{d-1}{2(d+1)}=\frac{1}{r_{*}}-\frac{s_{*}}{d},

Sobolev embedding gives

eitΔfLx2(d+1)/(d1)Cd|D|seitΔfLxr.\displaystyle\|e^{it\Delta}f\|_{L_{x}^{2(d+1)/(d-1)}}\leq C_{d}\||D|^{s_{*}}e^{it\Delta}f\|_{L_{x}^{r_{*}}}.

Combining inequalities above, we have

eitΔfLt,x2(d+1)/(d1)CdfH˙sfL21sfL2sfH1,\displaystyle\|e^{it\Delta}f\|_{L_{t,x}^{2(d+1)/(d-1)}}\leq C_{d}\|f\|_{\dot{H}^{s_{*}}}\leq\|f\|_{L^{2}}^{1-s_{*}}\|\nabla f\|_{L^{2}}^{s_{*}}\leq\|f\|_{H^{1}},

which is (2.6). Raising (2.6) to the power 2(d+1)/(d1)2(d+1)/(d-1) proves (2.7), and (2.8) follows from this directly. ∎

We record Lemma 2.3, a modified version of the dispersive estimate. It follows directly from (2.2) if the initial data is in L2(d+1)d+3L^{\frac{2(d+1)}{d+3}}, but needs a bit more justification if we only require the L2L^{2}-level decay of the initial data.

Lemma 2.3.

Let d5d\geq 5. For every fH1(d)f\in H^{1}(\mathbb{R}^{d}),

(2.9) lim|t|eitΔfL2(d+1)/(d1)(d)=0.\displaystyle\lim_{|t|\to\infty}\|e^{it\Delta}f\|_{L^{2(d+1)/(d-1)}(\mathbb{R}^{d})}=0.
Proof.

Since 2<2(d+1)/(d1)<2d/(d2)2<2(d+1)/(d-1)<2d/(d-2), Sobolev embedding gives

(2.10) hL2(d+1)/(d1)CdhH1.\displaystyle\|h\|_{L^{2(d+1)/(d-1)}}\leq C_{d}\|h\|_{H^{1}}.

For a Schwartz function ϕ\phi, the case p=2(d+1)/(d1)p=2(d+1)/(d-1) of (2.2) gives

(2.11) eitΔϕL2(d+1)/(d1)Cd|t|d/(d+1)ϕL2(d+1)/(d+3).\displaystyle\|e^{it\Delta}\phi\|_{L^{2(d+1)/(d-1)}}\leq C_{d}|t|^{-d/(d+1)}\|\phi\|_{L^{2(d+1)/(d+3)}}.

For fH1f\in H^{1}, choose Schwartz ϕ\phi with fϕH1<ϵ\|f-\phi\|_{H^{1}}<\epsilon. The H1H^{1}-unitarity of eitΔe^{it\Delta} and (2.10) give

supteitΔ(fϕ)L2(d+1)/(d1)Cdϵ.\displaystyle\sup_{t}\|e^{it\Delta}(f-\phi)\|_{L^{2(d+1)/(d-1)}}\leq C_{d}\epsilon.

Using (2.11), eitΔϕL2(d+1)/(d1)ϵ\|e^{it\Delta}\phi\|_{L^{2(d+1)/(d-1)}}\leq\epsilon for sufficiently large |t||t|. Therefore

lim sup|t|eitΔfL2(d+1)/(d1)(Cd+1)ϵ.\displaystyle\limsup_{|t|\to\infty}\|e^{it\Delta}f\|_{L^{2(d+1)/(d-1)}}\leq(C_{d}+1)\epsilon.

Sending ϵ0\epsilon\to 0 proves (2.9). ∎

Next we discuss the behavior of the measure given by the mass and current density. The intuition behind the limits below is that if we expand a scattering solution to the Schrödinger equation, then the leading order term (with x/2tx/2t as the variable) is the Fourier transform of the initial data. The convention for the Fourier transform we use in this paper is

f^(ξ)=(2π)ddeixξf(x)dx.\displaystyle\widehat{f}(\xi)=(2\pi)^{-d}\int_{\mathbb{R}^{d}}e^{-ix\cdot\xi}f(x)\,dx.
Proposition 2.4.

Let d1d\geq 1, u0H1(d)u_{0}\in H^{1}(\mathbb{R}^{d}), vu0(t)=eitΔu0v_{u_{0}}(t)=e^{it\Delta}u_{0} for t>0t>0, and let jj be as in (2.29). Define the measure ρu0\rho_{u_{0}} by

dΨ(y)dρu0(y)=(2π)ddΨ(2ξ)|u0^(ξ)|2𝑑ξ\displaystyle\int_{\mathbb{R}^{d}}\Psi(y)\,d\rho_{u_{0}}(y)=(2\pi)^{d}\int_{\mathbb{R}^{d}}\Psi(2\xi)|\widehat{u_{0}}(\xi)|^{2}\,d\xi

for bounded continuous Ψ\Psi. Define μtu0\mu_{t}^{u_{0}} and Jtu0J_{t}^{u_{0}} by

(2.12) Ψdμtu0=dΨ(x/t)|vu0(t,x)|2𝑑x,\displaystyle\int\Psi\,d\mu_{t}^{u_{0}}=\int_{\mathbb{R}^{d}}\Psi(x/t)|v_{u_{0}}(t,x)|^{2}\,dx,
(2.13) ΨdJtu0=dΨ(x/t)j(vu0(t))(x)𝑑x.\displaystyle\int\Psi\,dJ_{t}^{u_{0}}=\int_{\mathbb{R}^{d}}\Psi(x/t)j(v_{u_{0}}(t))(x)\,dx.

Then, as tt\to\infty,

(2.14) μtu0\displaystyle\mu_{t}^{u_{0}} ρu0,Jtu0y2ρu0in total variation.\displaystyle\to\rho_{u_{0}},\;J_{t}^{u_{0}}\to\frac{y}{2}\rho_{u_{0}}\;\text{in total variation}.

Moreover, μtu0\mu_{t}^{u_{0}} and |Jtu0||J_{t}^{u_{0}}| are uniformly tight for all sufficiently large tt: for every ϵ>0\epsilon>0, there are R,T>0R,T>0 such that, for all tTt\geq T,

(2.15) μtu0({yd:|y|>R})+|Jtu0|({yd:|y|>R})<ϵ.\displaystyle\mu_{t}^{u_{0}}(\{y\in\mathbb{R}^{d}:|y|>R\})+|J_{t}^{u_{0}}|(\{y\in\mathbb{R}^{d}:|y|>R\})<\epsilon.
Proof.

Set

gt(x)=ei|x|2/(4t)u0(x).\displaystyle g_{t}(x)=e^{i|x|^{2}/(4t)}u_{0}(x).

Since vu0=eitΔu0v_{u_{0}}=e^{it\Delta}u_{0}, we have

vu0(t,2tξ)=(πt)d/2γdeit|ξ|2g^t(ξ),|γd|=1.\displaystyle v_{u_{0}}(t,2t\xi)=\Big(\frac{\pi}{t}\Big)^{d/2}\gamma_{d}e^{it|\xi|^{2}}\widehat{g}_{t}(\xi),\qquad|\gamma_{d}|=1.

After the change of variables x=2tξx=2t\xi,

(2.16) Ψ(x/t)|vu0(t,x)|2𝑑x=(2π)dΨ(2ξ)|g^t(ξ)|2𝑑ξ.\displaystyle\int\Psi(x/t)|v_{u_{0}}(t,x)|^{2}\,dx=(2\pi)^{d}\int\Psi(2\xi)|\widehat{g}_{t}(\xi)|^{2}\,d\xi.

Multiplication by ei|x|2/(4t)e^{i|x|^{2}/(4t)} converges strongly to the identity on L2L^{2}. Thus gtu0g_{t}\to u_{0} in L2L^{2}, g^tu0^\widehat{g}_{t}\to\widehat{u_{0}} in L2L^{2}, and

(2.17) |g^t|2|u0^|2L1(2π)d(gtL2+u0L2)gtu0L20.\displaystyle\big\||\widehat{g}_{t}|^{2}-|\widehat{u_{0}}|^{2}\big\|_{L^{1}}\leq(2\pi)^{-d}(\|g_{t}\|_{L^{2}}+\|u_{0}\|_{L^{2}})\|g_{t}-u_{0}\|_{L^{2}}\to 0.

Equations (2.16)–(2.17) prove total-variation convergence of the mass measures.

Suppose first that u0u_{0} is Schwartz. A direct computation gives (with ξ=x/(2t)\xi=x/(2t))

j(vu0(t))(x)=(πt)d(ξ|g^t(ξ)|2+12t(g^t(ξ)¯ξg^t(ξ))).\displaystyle j(v_{u_{0}}(t))(x)=\Big(\frac{\pi}{t}\Big)^{d}\Big(\xi|\widehat{g}_{t}(\xi)|^{2}+\frac{1}{2t}\Im(\overline{\widehat{g}_{t}(\xi)}\nabla_{\xi}\widehat{g}_{t}(\xi))\Big).

Consequently, we have

(2.18) Ψ(x/t)dJtu0=(2π)dΨ(2ξ)ξ|g^t(ξ)|2𝑑ξ+Rt,\displaystyle\int\Psi(x/t)\,dJ_{t}^{u_{0}}=(2\pi)^{d}\int\Psi(2\xi)\xi|\widehat{g}_{t}(\xi)|^{2}\,d\xi+R_{t},

where RtR_{t} satisfies |Rt|ΨL2tu0L2xu0L2|R_{t}|\leq\frac{\|\Psi\|_{L^{\infty}}}{2t}\|u_{0}\|_{L^{2}}\|xu_{0}\|_{L^{2}}. For Schwartz u0u_{0}, we have gtu0g_{t}\to u_{0} in H1H^{1}. This is because

gt=ei|x|2/(4t)(u0+ix2tu0).\displaystyle\nabla g_{t}=e^{i|x|^{2}/(4t)}\Big(\nabla u_{0}+\frac{ix}{2t}u_{0}\Big).

This gives gtu0g_{t}\to u_{0} in H1H^{1} for Schwartz u0u_{0} as tt\to\infty by dominated convergence. Taking Fourier transform, we have ξg^tξu0^\xi\widehat{g}_{t}\to\xi\widehat{u_{0}} in L2L^{2}, and consequently

ξ|g^t|2ξ|u0^|2in L1.\displaystyle\xi|\widehat{g}_{t}|^{2}\to\xi|\widehat{u_{0}}|^{2}\quad\text{in }L^{1}.

Therefore the right-hand side of (2.18) (hence the left-hand side as well) converges to

(2π)dΨ(2ξ)ξ|u0^(ξ)|2𝑑ξ=Ψ(y)y2dρu0(y).\displaystyle(2\pi)^{d}\int\Psi(2\xi)\xi|\widehat{u_{0}}(\xi)|^{2}\,d\xi=\int\Psi(y)\frac{y}{2}\,d\rho_{u_{0}}(y).

For general u0H1u_{0}\in H^{1}, choose Schwartz u0,nu0u_{0,n}\to u_{0} in H1H^{1}. Uniformly in tt,

(2.19) |Ψd(μtu0μtu0,n)|ΨL(u0L2+u0,nL2)u0u0,nL2,\displaystyle|\int\Psi\,d(\mu_{t}^{u_{0}}-\mu_{t}^{u_{0,n}})|\leq\|\Psi\|_{L^{\infty}}(\|u_{0}\|_{L^{2}}+\|u_{0,n}\|_{L^{2}})\|u_{0}-u_{0,n}\|_{L^{2}},

and

(2.20) |Ψd(Jtu0Jtu0,n)|ΨL(u0u0,nL2u0L2+u0,nL2(u0u0,n)L2).\displaystyle|\int\Psi\,d(J_{t}^{u_{0}}-J_{t}^{u_{0,n}})|\leq\|\Psi\|_{L^{\infty}}(\|u_{0}-u_{0,n}\|_{L^{2}}\|\nabla u_{0}\|_{L^{2}}+\|u_{0,n}\|_{L^{2}}\|\nabla(u_{0}-u_{0,n})\|_{L^{2}}).

The estimates (2.19) and (2.20) also compare the two limiting measures on the right-hand side of (2.14). First let tt\to\infty and then nn\to\infty.

Now we prove the tightness in (2.15). For every measurable set EE, by Cauchy–Schwarz we have

(2.21) |Jtu0|(E)μtu0(E)1/2u0L2.\displaystyle|J_{t}^{u_{0}}|(E)\leq\mu_{t}^{u_{0}}(E)^{1/2}\|\nabla u_{0}\|_{L^{2}}.

Then (2.14) controls the mass term in (2.15), while (2.21) controls the current term. ∎

The following lemma shows that the majority of the mass of the solution can’t move arbitrarily fast.

Lemma 2.5.

Let d1d\geq 1, p>0p>0, and uu be a solution to

(it+Δ)u=|u|pu.\displaystyle(i\partial_{t}+\Delta)u=-|u|^{p}u.

Suppose uC(0,H1(d))u\in C(\mathbb{R}_{\geq 0};H^{1}(\mathbb{R}^{d})) and B=supt0u(t)H1<B=\sup_{t\geq 0}\|u(t)\|_{H^{1}}<\infty. Let u+H1u_{+}\in H^{1}, and let v(t)v(t) and r(t)r(t) be as in (1.4). If ηC(d)\eta\in C^{\infty}(\mathbb{R}^{d}), 0η10\leq\eta\leq 1, η=0\eta=0 on {|y|1}\{|y|\leq 1\}, and η=1\eta=1 on {|y|2}\{|y|\geq 2\}, then

(2.22) limRsupt1(η(x/(Rt))u(t)L2+η(x/(Rt))v(t)L2+η(x/(Rt))r(t)L2)=0.\displaystyle\lim_{R\to\infty}\sup_{t\geq 1}(\|\eta(x/(Rt))u(t)\|_{L^{2}}+\|\eta(x/(Rt))v(t)\|_{L^{2}}+\|\eta(x/(Rt))r(t)\|_{L^{2}})=0.
Proof.

Fix R1R\geq 1 and t1t\geq 1. For s[0,t]s\in[0,t], define

FR,t(s)=η(x/(Rt))2|u(s,x)|2𝑑x.\displaystyle F_{R,t}(s)=\int\eta(x/(Rt))^{2}|u(s,x)|^{2}\,dx.

Differentiating in ss and using (1.1), we know that ddsFR,t(s)\frac{d}{ds}F_{R,t}(s) is an integral involving uu and u\nabla u, and we have 22 2 Rigorously speaking, we need to approximate by nice solutions first.

|ddsFR,t(s)|CηRtu(s)L2u(s)L2CηB2Rt.\displaystyle\big|\frac{d}{ds}F_{R,t}(s)\big|\leq\frac{C_{\eta}}{Rt}\|u(s)\|_{L^{2}}\|\nabla u(s)\|_{L^{2}}\leq\frac{C_{\eta}B^{2}}{Rt}.

Integrating, we obtain

(2.23) FR,t(t)|x|Rt|u(0,x)|2𝑑x+CηB2R|x|R|u(0,x)|2𝑑x+CηB2R.\displaystyle F_{R,t}(t)\leq\int_{|x|\geq Rt}|u(0,x)|^{2}\,dx+\frac{C_{\eta}B^{2}}{R}\leq\int_{|x|\geq R}|u(0,x)|^{2}\,dx+\frac{C_{\eta}B^{2}}{R}.

The right-hand side of (2.23) is independent of tt and tends to zero as RR\to\infty.

Applying the same argument in the proof of (2.23) to vv gives

supt1η(x/(Rt))v(t)L22|x|R|u+(x)|2𝑑x+Cηu+H12R.\displaystyle\sup_{t\geq 1}\|\eta(x/(Rt))v(t)\|_{L^{2}}^{2}\leq\int_{|x|\geq R}|u_{+}(x)|^{2}\,dx+\frac{C_{\eta}\|u_{+}\|_{H^{1}}^{2}}{R}.

Finally, η(x/(Rt))r=η(x/(Rt))uη(x/(Rt))v\eta(x/(Rt))r=\eta(x/(Rt))u-\eta(x/(Rt))v, so applying the triangle inequality proves (2.22). ∎

2.2. Compact attractor

In this subsection, we recall some results from [23]. For every ff for which the NLS flow is defined, denote the solution at time tt with initial data ff by S(t)fS(t)f. A set KH1(d)K\subset H^{1}(\mathbb{R}^{d}) is precompact modulo at most JJ translations if there is a compact set K0H1(d)K_{0}\subset H^{1}(\mathbb{R}^{d}) such that every element of KK is a sum of at most JJ translates of elements of K0K_{0}. Then we have the following decay estimate that is uniform over KK.

Lemma 2.6.

[23, Corollary B.6] Let KH1(d)K\subset H^{1}(\mathbb{R}^{d}) be precompact modulo at most JJ translations. For every 2<q2d/(d2)2<q\leq 2d/(d-2),

(2.24) limt±supfKeitΔfLq\displaystyle\lim_{t\to\pm\infty}\sup_{f\in K}\|e^{it\Delta}f\|_{L^{q}} =0.\displaystyle=0.

For every R>0R>0,

(2.25) limt±supfKsupx0d|xx0|R(|eitΔf(x)|2+|eitΔf(x)|2)dx\displaystyle\lim_{t\to\pm\infty}\sup_{f\in K}\sup_{x_{0}\in\mathbb{R}^{d}}\int_{|x-x_{0}|\leq R}(|e^{it\Delta}f(x)|^{2}+|\nabla e^{it\Delta}f(x)|^{2})\,dx =0.\displaystyle=0.

Then we come to an important black box for us, the profile decomposition of the solutions to (1.1).

Theorem 2.7.

Let d5d\geq 5, and let uC(0,H1(d))u\in C(\mathbb{R}_{\geq 0};H^{1}(\mathbb{R}^{d})) solve

(it+Δ)u=|u|4/(d1)u\displaystyle(i\partial_{t}+\Delta)u=-|u|^{4/(d-1)}u

and satisfy (1.2). There exist an integer J1J\geq 1, a unique radiation state u+H1(d)u_{+}\in H^{1}(\mathbb{R}^{d}), and a closed translation-invariant set 𝒦\mathcal{K}, precompact modulo at most JJ translations. The forward NLS solution operator S(t)S(t) is defined on 𝒦\mathcal{K} and leaves it forward invariant. Writing J𝒦J\mathcal{K} for the set of sums of JJ elements of 𝒦\mathcal{K}, one has

distH1(u(t)eitΔu+,J𝒦)0(t).\displaystyle\operatorname{dist}_{H^{1}}\big(u(t)-e^{it\Delta}u_{+},J\mathcal{K}\big)\to 0\qquad(t\to\infty).

Moreover, for every sequence tnt_{n}\to\infty there are a subsequence, profiles w1,,wJ𝒦w_{1},\ldots,w_{J}\in\mathcal{K}, and centers x1,n,,xJ,ndx_{1,n},\ldots,x_{J,n}\in\mathbb{R}^{d} satisfying

(2.26) u(tn)=eitnΔu++j=1Jwj(xj,n)+oH1(1),\displaystyle u(t_{n})=e^{it_{n}\Delta}u_{+}+\sum_{j=1}^{J}w_{j}(\,\cdot-x_{j,n})+o_{H^{1}}(1),

and

(2.27) |xj,nxk,n|(jk).\displaystyle|x_{j,n}-x_{k,n}|\to\infty\qquad(j\neq k).

Verifying that our nonlinearity satisfies the conditions in [23] is standard, and in fact this is already listed as an example in [23] immediately after the conditions are stated.

Next we show that the soliton part can’t have arbitrarily large frequency. For ydy\in\mathbb{R}^{d}, we denote the translation action by τyf(x)=f(xy)\tau_{y}f(x)=f(x-y).

Proposition 2.8.

Let r(t)r(t) be as in (1.4). If P>RP_{>R} is the Fourier projection onto {|ξ|>R}\{|\xi|>R\}, then

(2.28) limRlim suptP>Rr(t)H1=0.\displaystyle\lim_{R\to\infty}\limsup_{t\to\infty}\|P_{>R}r(t)\|_{H^{1}}=0.
Proof.

By Theorem 2.7, the residual approaches a set contained in a finite sum of translates of one compact subset KH1K\subset H^{1}. Thus there are MM\in\mathbb{N} and η(t)0\eta(t)\to 0 such that, for all sufficiently large tt,

r(t)m=1Mτxm(t)km(t)H1η(t),km(t)K.\displaystyle\big\|r(t)-\sum_{m=1}^{M}\tau_{x_{m}(t)}k_{m}(t)\big\|_{H^{1}}\leq\eta(t),\qquad k_{m}(t)\in K.

Compactness of KK implies

supkKP>RkH10.\displaystyle\sup_{k\in K}\|P_{>R}k\|_{H^{1}}\to 0.

Indeed, take a finite ϵ\epsilon-net in KK, use dominated convergence in Fourier space for the net points, and use that P>RP_{>R} has norm at most one on H1H^{1}. Since Fourier projection commutes with translations,

P>Rr(t)H1η(t)+MsupkKP>RkH1.\displaystyle\|P_{>R}r(t)\|_{H^{1}}\leq\eta(t)+M\sup_{k\in K}\|P_{>R}k\|_{H^{1}}.

Taking the late-time limsup and then RR\to\infty proves (2.28). ∎

Next we estimate the current of the solution. Define:

(2.29) j(f):=(f¯f).\displaystyle j(f):=\Im(\overline{f}\nabla f).

Asymptotically, the soliton part and the radiation part tend to decouple, and we have the following asymptotic orthogonality.

Proposition 2.9.

Let v(t)v(t) and r(t)r(t) be as in (1.4). Then

(2.30) d(|r||v|+|r||v|+|v||r|)𝑑x0.\displaystyle\int_{\mathbb{R}^{d}}\big(|r||v|+|r||\nabla v|+|v||\nabla r|\big)\,dx\to 0.

Consequently,

(2.31) |u|2|v|2|r|2L10,\displaystyle\big\||u|^{2}-|v|^{2}-|r|^{2}\big\|_{L^{1}}\to 0,
(2.32) j(u)j(v)j(r)L10.\displaystyle\|j(u)-j(v)-j(r)\|_{L^{1}}\to 0.
Proof.

For fL2f\in L^{2} and fixed R>0R>0,

(2.33) supydeitΔfL2(B(y,R))0.\displaystyle\sup_{y\in\mathbb{R}^{d}}\|e^{it\Delta}f\|_{L^{2}(B(y,R))}\to 0.

Approximate ff by a compactly supported smooth function, use (2.1) on the approximation, and use L2L^{2} unitarity for the error. Apply (2.33) to u+u_{+} and each component of u+\nabla u_{+}.

For fixed wH1w\in H^{1}, arbitrary centers xnx_{n}, and tnt_{n}\to\infty, split w(xn)w(\cdot-x_{n}) into a fixed-radius ball and an H1H^{1}-small tail. Cauchy–Schwarz and (2.33) show

|w(xxn)||v(tn,x)|𝑑x0,\displaystyle\int|w(x-x_{n})||v(t_{n},x)|\,dx\to 0,
(2.34) |w(xxn)||v(tn,x)|𝑑x0,|w(xxn)||v(tn,x)|𝑑x0.\displaystyle\int|w(x-x_{n})||\nabla v(t_{n},x)|\,dx\to 0,\qquad\int|\nabla w(x-x_{n})||v(t_{n},x)|\,dx\to 0.

For an arbitrary sequence tnt_{n}\to\infty, apply Theorem 2.7:

(2.35) r(tn)=j=1Jwj(xj,n)+en,en0 in H1.\displaystyle r(t_{n})=\sum_{j=1}^{J}w_{j}(\cdot-x_{j,n})+e_{n},\qquad e_{n}\to 0\text{ in }H^{1}.

Thus (2.30) holds on every extracted subsequence and therefore along the full time variable. Also, we have

(|u|2|v|2)|r|2=2(v¯r),j(u)j(v)j(r)=(v¯r+r¯v).\displaystyle\big(|u|^{2}-|v|^{2}\big)-|r|^{2}=2\Re(\overline{v}r),\quad j(u)-j(v)-j(r)=\Im(\overline{v}\nabla r+\overline{r}\nabla v).

Taking L1L^{1} norms gives (2.31)–(2.32). Pushforward does not increase total variation. ∎

We consider the mass, energy, momentum of the soliton part below. By the decoupling property of the soliton part and the radiation part mentioned above, these quantities are obtained by subtracting the mass, energy, or momentum of the radiation part from that of the entire solution. This decoupling property will be used many times throughout the paper.

Proposition 2.10.

Let MM, EE, and PP be as in (1.3). Define

(2.36) m:=M(u)M(u+),e:=E(u)12u+L22,p:=P(u)P(u+),R:=me|p|22.\displaystyle\begin{split}m_{\sharp}&:=M(u)-M(u_{+}),\qquad e_{\sharp}:=E(u)-\frac{1}{2}\|\nabla u_{+}\|_{L^{2}}^{2},\\ p_{\sharp}&:=P(u)-P(u_{+}),\qquad R_{\sharp}:=m_{\sharp}e_{\sharp}-\frac{|p_{\sharp}|^{2}}{2}.\end{split}

For v(t)v(t) and r(t)r(t) as in (1.4),

(2.37) M(r(t))m,E(r(t))e,P(r(t))p.\displaystyle M(r(t))\to m_{\sharp},\qquad E(r(t))\to e_{\sharp},\qquad P(r(t))\to p_{\sharp}.
Proof.

The limits of M(r(t))M(r(t)) and P(r(t))P(r(t)) follow directly from Proposition 2.9. It remains to prove the energy limit. The argument used in the proof of Proposition 2.9 also gives

v(t),r(t)0.\displaystyle\langle\nabla v(t),\nabla r(t)\rangle\to 0.

More precisely, along an arbitrary sequence tnt_{n}\to\infty, take the profile decomposition (2.35). The argument proving (2.34), now applied to the components of u+\nabla u_{+} and wj\nabla w_{j}, makes the pairing of v(tn)\nabla v(t_{n}) with every translated profile tend to zero, while Cauchy–Schwarz handles the H1H^{1}-small remainder. Since the original sequence was arbitrary, the displayed limit holds for tt\to\infty.

We have v(t)L2(d+1)/(d1)0\|v(t)\|_{L^{2(d+1)/(d-1)}}\to 0 by Lemma 2.3. The pointwise power inequality and the uniform H1H^{1} bounds give

|u(t)L2(d+1)/(d1)2(d+1)/(d1)r(t)L2(d+1)/(d1)2(d+1)/(d1)|v(t)L2(d+1)/(d1)(u(t)L2(d+1)/(d1)(d+3)/(d1)+r(t)L2(d+1)/(d1)(d+3)/(d1))0.\displaystyle\big|\|u(t)\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}-\|r(t)\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}\big|\lesssim\|v(t)\|_{L^{2(d+1)/(d-1)}}\big(\|u(t)\|_{L^{2(d+1)/(d-1)}}^{(d+3)/(d-1)}+\|r(t)\|_{L^{2(d+1)/(d-1)}}^{(d+3)/(d-1)}\big)\to 0.

Using u=v+ru=v+r and the definition of the energy, we therefore have

E(u)12v(t)L22E(r(t))=\displaystyle E(u)-\frac{1}{2}\|\nabla v(t)\|_{L^{2}}^{2}-E(r(t))= v(t),r(t)\displaystyle\Re\langle\nabla v(t),\nabla r(t)\rangle
d12(d+1)(u(t)L2(d+1)/(d1)2(d+1)/(d1)r(t)L2(d+1)/(d1)2(d+1)/(d1))0.\displaystyle-\frac{d-1}{2(d+1)}\big(\|u(t)\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}-\|r(t)\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}\big)\to 0.

Finally, v(t)L2=u+L2\|\nabla v(t)\|_{L^{2}}=\|\nabla u_{+}\|_{L^{2}} by unitarity of the free flow. This proves the energy limit and hence (2.37). ∎

2.3. Finiteness of non-scattering directions

Now we prove that there can only be finitely many directions on which the localized product of mass and energy does not tend to zero in a narrow cone near these directions. To quantify this, for vdv\in\mathbb{R}^{d}, we define

(2.38) 𝒜(v):=lim infρ0lim inftM(χ(x/tvρ)u(t))E(χ(x/tvρ)u(t)).\displaystyle\mathcal{A}(v):=\liminf_{\rho\downarrow 0}\liminf_{t\to\infty}M\Big(\chi\Big(\frac{x/t-v}{\rho}\Big)u(t)\Big)E\Big(\chi\Big(\frac{x/t-v}{\rho}\Big)u(t)\Big).

Then our result is the following, which is a crucial step in proving Theorem 1.5.

Proposition 2.11.

Let uu be a global H1H^{1} solution to (1.1) satisfying (1.2). Fix a bump function χCc(d)\chi\in C_{c}^{\infty}(\mathbb{R}^{d}), with 0χ10\leq\chi\leq 1, supported in B(0,1)B(0,1), equal to one near the origin. For every c0>0c_{0}>0, there are only finitely many vdv\in\mathbb{R}^{d} such that 𝒜(v)c0\mathcal{A}(v)\geq c_{0}. More precisely,

(2.39) #{vd:𝒜(v)c0}u(0)L2supt0u(t)L22c0.\displaystyle\#\{v\in\mathbb{R}^{d}:\mathcal{A}(v)\geq c_{0}\}\leq\frac{\|u(0)\|_{L^{2}}\sup_{t\geq 0}\|\nabla u(t)\|_{L^{2}}}{\sqrt{2c_{0}}}.
Proof.

Take any NN distinct velocities v1,,vNv_{1},\ldots,v_{N} satisfying

(2.40) 𝒜(vj)c0,1jN.\displaystyle\mathcal{A}(v_{j})\geq c_{0},\qquad 1\leq j\leq N.

We prove a bound for NN independent of the chosen velocities. Fix 0<ϵ<c0/20<\epsilon<c_{0}/2. By (2.38) and (2.40), for each jj, whenever ρ>0\rho>0 is sufficiently small,

(2.41) lim inftM(χ(x/tvjρ)u(t))E(χ(x/tvjρ)u(t))c0ϵ.\displaystyle\liminf_{t\to\infty}M\Big(\chi\Big(\frac{x/t-v_{j}}{\rho}\Big)u(t)\Big)E\Big(\chi\Big(\frac{x/t-v_{j}}{\rho}\Big)u(t)\Big)\geq c_{0}-\epsilon.

Since v1,,vNv_{1},\ldots,v_{N} are distinct, we can choose sufficiently small ρ1,,ρN>0\rho_{1},\ldots,\rho_{N}>0 so that the balls B(vj,ρj)B(v_{j},\rho_{j}), 1jN1\leq j\leq N, are pairwise disjoint and

(2.42) lim inftM(χ(x/tvjρj)u(t))E(χ(x/tvjρj)u(t))c0ϵ\displaystyle\liminf_{t\to\infty}M\Big(\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\Big)E\Big(\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\Big)\geq c_{0}-\epsilon

for every 1jN1\leq j\leq N.

Since there are only finitely many jj, there exists T>0T>0 such that, for every tTt\geq T and every 1jN1\leq j\leq N,

(2.43) M(χ(x/tvjρj)u(t))E(χ(x/tvjρj)u(t))c02ϵ.\displaystyle M\Big(\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\Big)E\Big(\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\Big)\geq c_{0}-2\epsilon.

For the energy defined in (1.3), we have OPENE(x/tvjρj)u(t))12|(x/tvjρj)u(t)L22E(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t))\leq\frac{1}{2}\|\nabla\|\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\|_{L^{2}}^{2}. So

(2.44) χ(x/tvjρj)u(t)L2(χ(x/tvjρj)u(t))L22(c02ϵ).\displaystyle\big\|\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\big\|_{L^{2}}\big\|\nabla\Big(\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\Big)\big\|_{L^{2}}\geq\sqrt{2(c_{0}-2\epsilon)}.

Summing (2.44) over jj and applying the Cauchy–Schwarz inequality gives

(2.45) N2(c02ϵ)j=1Nχ(x/tvjρj)u(t)L2(χ(x/tvjρj)u(t))L2(j=1Nχ(x/tvjρj)u(t)L22)1/2(j=1N(χ(x/tvjρj)u(t))L22)1/2.\displaystyle\begin{split}N\sqrt{2(c_{0}-2\epsilon)}&\leq\sum_{j=1}^{N}\big\|\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\big\|_{L^{2}}\big\|\nabla\Big(\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\Big)\big\|_{L^{2}}\\ &\leq\Big(\sum_{j=1}^{N}\big\|\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\big\|_{L^{2}}^{2}\Big)^{1/2}\Big(\sum_{j=1}^{N}\big\|\nabla\Big(\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\Big)\big\|_{L^{2}}^{2}\Big)^{1/2}.\end{split}

By construction, the supports of the cutoffs are pairwise disjoint for every t>0t>0, and 0χ10\leq\chi\leq 1. Therefore, by conservation of mass,

(2.46) j=1Nχ(x/tvjρj)u(t)L22\displaystyle\sum_{j=1}^{N}\big\|\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\big\|_{L^{2}}^{2} =dj=1Nχ(x/tvjρj)2|u(t,x)|2𝑑xu(t)L22=u(0)L22.\displaystyle=\int_{\mathbb{R}^{d}}\sum_{j=1}^{N}\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)^{2}|u(t,x)|^{2}\,dx\leq\|u(t)\|_{L^{2}}^{2}=\|u(0)\|_{L^{2}}^{2}.

Moreover,

(2.47) (χ(x/tvjρj)u(t))=χ(x/tvjρj)u(t)+1tρj(χ)(x/tvjρj)u(t).\displaystyle\nabla\Big(\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\Big)=\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)\nabla u(t)+\frac{1}{t\rho_{j}}(\nabla\chi)\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t).

Using again the pairwise disjointness of the supports and the triangle inequality in the direct sum over jj, we obtain

(2.48) (j=1N(χ(x/tvjρj)u(t))L22)1/2u(t)L2+χLtmax1jN1ρju(0)L2.\displaystyle\Big(\sum_{j=1}^{N}\big\|\nabla\Big(\chi\Big(\frac{x/t-v_{j}}{\rho_{j}}\Big)u(t)\Big)\big\|_{L^{2}}^{2}\Big)^{1/2}\leq\|\nabla u(t)\|_{L^{2}}+\frac{\|\nabla\chi\|_{L^{\infty}}}{t}\max_{1\leq j\leq N}\frac{1}{\rho_{j}}\|u(0)\|_{L^{2}}.

Combining (2.45), (2.46), and (2.48), for every tTt\geq T we have

(2.49) N2(c02ϵ)u(0)L2(sups0u(s)L2+χLtmax1jN1ρju(0)L2).\displaystyle N\sqrt{2(c_{0}-2\epsilon)}\leq\|u(0)\|_{L^{2}}\Big(\sup_{s\geq 0}\|\nabla u(s)\|_{L^{2}}+\frac{\|\nabla\chi\|_{L^{\infty}}}{t}\max_{1\leq j\leq N}\frac{1}{\rho_{j}}\|u(0)\|_{L^{2}}\Big).

Letting tt\to\infty in (2.49) gives

(2.50) N2(c02ϵ)u(0)L2supt0u(t)L2.\displaystyle N\sqrt{2(c_{0}-2\epsilon)}\leq\|u(0)\|_{L^{2}}\sup_{t\geq 0}\|\nabla u(t)\|_{L^{2}}.

Finally, letting ϵ0\epsilon\to 0 in (2.50), we obtain (2.39). ∎

3. Preliminaries about the soliton and the upgrading mechanism

We recall some basic facts about the soliton in this section and prove some preliminary results about the upgrading mechanism that improves an estimate along a time sequence to an estimate that is uniform for all large times.

Proposition 3.1 (Ground-state identities and threshold trapping).

Let d5d\geq 5, let MM and EE be defined by (1.3), and let QQ be the ground state in (1.5). Then

(3.1) QL22+QL22=QL2(d+1)/(d1)2(d+1)/(d1)=(d+1)M(Q),QL22=dM(Q),E(Q)=12M(Q)=12dQL22,M(Q)E(Q)=(QL2QL2)22d>0.\displaystyle\begin{aligned} \|\nabla Q\|_{L^{2}}^{2}+\|Q\|_{L^{2}}^{2}&=\|Q\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}=(d+1)M(Q),\\ \|\nabla Q\|_{L^{2}}^{2}&=dM(Q),\qquad E(Q)=\frac{1}{2}M(Q)=\frac{1}{2d}\|\nabla Q\|_{L^{2}}^{2},\\ M(Q)E(Q)&=\frac{(\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}})^{2}}{2d}>0.\end{aligned}

For every fH1(d)f\in H^{1}(\mathbb{R}^{d}), we have the sharp Gagliardo–Nirenberg inequality

(3.2) fL2(d+1)/(d1)2(d+1)/(d1)\displaystyle\|f\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)} CGNfL22/(d1)fL22d/(d1),\displaystyle\leq C_{\mathrm{GN}}\|f\|_{L^{2}}^{2/(d-1)}\|\nabla f\|_{L^{2}}^{2d/(d-1)},

where

(3.3) CGN\displaystyle C_{\mathrm{GN}} =d+1d(QL2QL2)2/(d1).\displaystyle=\frac{d+1}{d}(\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}})^{-2/(d-1)}.

Set

(3.4) s(f)\displaystyle s(f) =fL2fL2QL2QL2.\displaystyle=\frac{\|f\|_{L^{2}}\|\nabla f\|_{L^{2}}}{\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}}}.

Then, without assuming a mass-energy bound,

(3.5) E(f)\displaystyle E(f) (12d12ds(f)2/(d1))fL22.\displaystyle\geq\Big(\frac{1}{2}-\frac{d-1}{2d}s(f)^{2/(d-1)}\Big)\|\nabla f\|_{L^{2}}^{2}.

In particular, E(f)0E(f)\geq 0 whenever s(f)1s(f)\leq 1. If, in addition, for some δ(0,1)\delta\in(0,1),

(3.6) M(f)E(f)\displaystyle M(f)E(f) (1δ)M(Q)E(Q),s(f)<1,\displaystyle\leq(1-\delta)M(Q)E(Q),\qquad s(f)<1,

then

(3.7) fL2fL2\displaystyle\|f\|_{L^{2}}\|\nabla f\|_{L^{2}} 1δQL2QL2.\displaystyle\leq\sqrt{1-\delta}\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}}.
Proof.

Multiplying (1.5) by QQ and integrating by parts gives

(3.8) QL22+QL22\displaystyle\|\nabla Q\|_{L^{2}}^{2}+\|Q\|_{L^{2}}^{2} =QL2(d+1)/(d1)2(d+1)/(d1).\displaystyle=\|Q\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}.

The Pohozaev identity (in the form given in [3]) gives

(3.9) d22QL22+d2QL22\displaystyle\frac{d-2}{2}\|\nabla Q\|_{L^{2}}^{2}+\frac{d}{2}\|Q\|_{L^{2}}^{2} =d(d1)2(d+1)QL2(d+1)/(d1)2(d+1)/(d1).\displaystyle=\frac{d(d-1)}{2(d+1)}\|Q\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}.

Combining (3.8) and (3.9) gives

QL22=dQL22,QL2(d+1)/(d1)2(d+1)/(d1)=(d+1)QL22,\displaystyle\|\nabla Q\|_{L^{2}}^{2}=d\|Q\|_{L^{2}}^{2},\qquad\|Q\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}=(d+1)\|Q\|_{L^{2}}^{2},

hence E(Q)=12QL22E(Q)=\frac{1}{2}\|Q\|_{L^{2}}^{2}, and all the formulas in (3.1) follow.

Now we turn to the remaining claims. Weinstein’s result [26] shows that the sharp constant in (3.2) is attained by a positive ground state, which solves (1.5); Kwong’s uniqueness theorem then identifies it with QQ [19]. Hence equality holds in (3.2) for QQ. Combining this with (3.1) yields (3.3). Substituting this constant into (3.2) gives (3.5). The coefficient on the right-hand side of (3.5) is nonnegative when s(f)1s(f)\leq 1, proving the E(f)0E(f)\geq 0 part.

Assume now (3.6). Multiplying (3.5) by M(f)M(f) and combining with (3.1) gives

(3.10) M(f)E(f)M(Q)E(Q)\displaystyle\frac{M(f)E(f)}{M(Q)E(Q)} ds(f)2(d1)s(f)2d/(d1)=s(f)2(d(d1)s(f)2/(d1))s(f)2,\displaystyle\geq ds(f)^{2}-(d-1)s(f)^{2d/(d-1)}=s(f)^{2}(d-(d-1)s(f)^{2/(d-1)})\geq s(f)^{2},

where the last inequality used s(f)<1s(f)<1. Combining this with (3.6) gives s(f)1δs(f)\leq\sqrt{1-\delta}, which is (3.7). ∎

As one would expect, if a nonlinear solution has a compact trajectory (say in H1H^{1}), then it should not scatter and it should be ‘above threshold’. This is indeed the case as we shall show now. On the other hand, the bound M(w)E(w)ΘM(w)E(w)\geq\Theta is not Galilean invariant. For fixed M(w)M(w), the quantity E(w)E(w) is minimized in the ‘rest frame’, that is the one in which the momentum is zero. Then being above threshold in this frame gives a stronger inequality in terms of our original ww. It is this form below that enters our proof in Section 5.

Proposition 3.2.

Let MM, EE, PP, and Θ\Theta be as in (1.3) and (1.6). Then every w𝒦{0}w\in\mathcal{K}\setminus\{0\} satisfies

(3.11) E(w)|P(w)|22M(w)>0,M(w)E(w)|P(w)|22Θ.\displaystyle E(w)-\frac{|P(w)|^{2}}{2M(w)}>0,\qquad M(w)E(w)-\frac{|P(w)|^{2}}{2}\geq\Theta.
Proof.

Fix w𝒦{0}w\in\mathcal{K}\setminus\{0\}, and let W(t)=S(t)w𝒦W(t)=S(t)w\in\mathcal{K} be the solution of the NLS in Theorem 2.7 with W(0)=wW(0)=w. Since 𝒦\mathcal{K} is precompact modulo at most JJ translations, WW is uniformly bounded in H1H^{1} for all t0t\geq 0. Let PP be as in (1.3) and w0w\neq 0. Set ξ=P(w)/M(w)\xi=-P(w)/M(w) and define

(3.12) W~(t,x)=ei(xξt|ξ|2)W(t,x2ξt),\displaystyle\widetilde{W}(t,x)=e^{i(x\cdot\xi-t|\xi|^{2})}W(t,x-2\xi t),

which solves (1.1), and is still uniformly bounded in H1H^{1}. In particular, we have

M(W~(0))\displaystyle M(\widetilde{W}(0)) =M(w),P(W~(0))=0,\displaystyle=M(w),\qquad P(\widetilde{W}(0))=0,
E(W~(0))\displaystyle E(\widetilde{W}(0)) =E(w)+ξP(w)+12|ξ|2M(w)=E(w)|P(w)|22M(w).\displaystyle=E(w)+\xi\cdot P(w)+\frac{1}{2}|\xi|^{2}M(w)=E(w)-\frac{|P(w)|^{2}}{2M(w)}.

When E(W~(0))<0E(\widetilde{W}(0))<0, the finite time blow-up is well-known (see [13]).

If E(W~(0))=0E(\widetilde{W}(0))=0, then [14, Lemma 2.16] rules out Part I of [14, Theorem A*] for nonzero data, while Part II contradicts the same global boundedness. Hence E(W~(0))>0E(\widetilde{W}(0))>0, which is the first inequality in (3.11). Multiplying the displayed identity for E(W~(0))E(\widetilde{W}(0)) by M(W~(0))=M(w)M(\widetilde{W}(0))=M(w) gives

M(W~(0))E(W~(0))=M(w)E(w)|P(w)|22.\displaystyle M(\widetilde{W}(0))E(\widetilde{W}(0))=M(w)E(w)-\frac{|P(w)|^{2}}{2}.

If M(w)E(w)|P(w)|22<ΘM(w)E(w)-\frac{|P(w)|^{2}}{2}<\Theta, then M(W~(0))E(W~(0))<ΘM(\widetilde{W}(0))E(\widetilde{W}(0))<\Theta. Combining with the uniform boundedness in H1H^{1}, [14, Theorems A and A*] show that W~\widetilde{W} scatters in the forward direction, and so does WW.

If WW scattered forward, then for some ϕ+H1\phi_{+}\in H^{1}, eitΔW(t)ϕ+e^{-it\Delta}W(t)\to\phi_{+} in H1H^{1}, and hence in LrL^{r} for 2<r<2d/(d2)2<r<2d/(d-2). However, since W(t)𝒦W(t)\in\mathcal{K}, Lemma 2.6 gives

(3.13) eitΔW(t)Lrsupf𝒦eitΔfLr0(t).\displaystyle\|e^{-it\Delta}W(t)\|_{L^{r}}\leq\sup_{f\in\mathcal{K}}\|e^{-it\Delta}f\|_{L^{r}}\to 0\qquad(t\to\infty).

So we have ϕ+=0\phi_{+}=0 and the conservation of mass shows

M(w)=eitΔW(t)L22ϕ+L22=0,\displaystyle M(w)=\|e^{-it\Delta}W(t)\|_{L^{2}}^{2}\to\|\phi_{+}\|_{L^{2}}^{2}=0,

whcih is a contradiction and completes the proof. ∎

Two results below are the starting points of our upgrading mechanism. We will prove our upgrading estimates by contradiction. The goal of this part is to show that if the uniform estimate for all large times fails, then the mass m=M(u)M(u+)m_{\sharp}=M(u)-M(u_{+}) of the soliton part is strictly positive. For why this is useful, see Proposition 5.9 below and the proof of Theorem 1.3 in Section 8.1.

Lemma 3.3.

Let d1d\geq 1, uC(1,H1)u\in C(\mathbb{R}_{\geq 1};H^{1}), u+H1u_{+}\in H^{1}, and vv and rr as in (1.4). Let χ0,χ1Cc(d)\chi_{0},\chi_{1}\in C_{c}^{\infty}(\mathbb{R}^{d}) satisfy χ1=1\chi_{1}=1 on suppχ0\mathrm{supp}\chi_{0}. Suppose TnT_{n}\to\infty and

χ1(x/Tn)r(Tn)H10.\displaystyle\|\chi_{1}(x/T_{n})r(T_{n})\|_{H^{1}}\to 0.

If

(3.14) lim supnsuptTnχ0(x/t)r(t)H1>0,\displaystyle\limsup_{n\to\infty}\sup_{t\geq T_{n}}\|\chi_{0}(x/t)r(t)\|_{H^{1}}>0,

where the limsup may be infinite, then after passage to a subsequence there are ϵ>0\epsilon>0 and finite sn>Tns_{n}>T_{n}, with sns_{n}\to\infty, such that

χ0(x/sn)r(sn)H1=ϵ,\displaystyle\|\chi_{0}(x/s_{n})r(s_{n})\|_{H^{1}}=\epsilon,

and we have χ0(x/t)r(t)H1<ϵ\|\chi_{0}(x/t)r(t)\|_{H^{1}}<\epsilon for t(Tn,sn)t\in(T_{n},s_{n}).

Proof.

Since χ0=χ0χ1\chi_{0}=\chi_{0}\chi_{1}, we have

χ0(x/Tn)r(Tn)H1Cχ0χ1(x/Tn)r(Tn)H10.\displaystyle\|\chi_{0}(x/T_{n})r(T_{n})\|_{H^{1}}\leq C_{\chi_{0}}\|\chi_{1}(x/T_{n})r(T_{n})\|_{H^{1}}\to 0.

Since r(t)r(t) is H1H^{1}-continuous in tt, we know χ0(x/t)r(t)H1\|\chi_{0}(x/t)r(t)\|_{H^{1}} is continuous in tt.

If (3.14) is true, then the quantity will fluctuate between 0 and lim supnsuptTnχ0(x/t)r(t)H1\limsup_{n\to\infty}\sup_{t\geq T_{n}}\|\chi_{0}(x/t)r(t)\|_{H^{1}} infinitely many times. Then, fixing a sufficiently small ϵ\epsilon, we have χ0(x/Tn)r(Tn)H1<ϵ\|\chi_{0}(x/T_{n})r(T_{n})\|_{H^{1}}<\epsilon for large nn. Then we define

sn=inf{tTn:χ0(x/t)r(t)H1ϵ},\displaystyle s_{n}=\inf\{t\geq T_{n}:\|\chi_{0}(x/t)r(t)\|_{H^{1}}\geq\epsilon\},

and sns_{n} will have the claimed property. ∎

The lemma above gives:

Proposition 3.4.

Let the setting be as in Lemma 3.3. Suppose (1.15) fails (or equivalently (3.14) is true). Then m=M(u)M(u+)>0m_{\sharp}=M(u)-M(u_{+})>0.

Proof.

We have M(r(t))mM(r(t))\to m_{\sharp} by Proposition 2.10. If m=0m_{\sharp}=0, arguing as in the discussion after (1.22), we know r(t)H10\|r(t)\|_{H^{1}}\to 0, contradicting the failure of (1.15). Thus m>0m_{\sharp}>0. ∎

4. Below threshold scattering in a spacetime cone

We prove Theorem 1.1 in this section. We first recall the scattering result to which we resort.

Theorem 4.1.

[6, Theorem 1.1], [14, 10] Let d3d\geq 3. For u0H1(d)u_{0}\in H^{1}(\mathbb{R}^{d}) such that

(4.1) M(u0)E(u0)\displaystyle M(u_{0})E(u_{0}) <M(Q)E(Q),\displaystyle<M(Q)E(Q),
(4.2) u0L2u0L2\displaystyle\|u_{0}\|_{L^{2}}\|\nabla u_{0}\|_{L^{2}} <QL2QL2,\displaystyle<\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}},

the solution uu of (1.1) with u(0)=u0u(0)=u_{0} is global and scatters in both time directions. We consider only the forward direction: there exists u+H1(d)u_{+}\in H^{1}(\mathbb{R}^{d}) such that

(4.3) u(t)eitΔu+H1\displaystyle\|u(t)-e^{it\Delta}u_{+}\|_{H^{1}} 0as t.\displaystyle\to 0\quad\text{as }t\to\infty.

Using Theorem 4.1, we now prove Theorem 1.1.

Proof of Theorem 1.1.

We first apply Theorem 2.7 to uu and use notation there. Take a sequence tnt_{n}\to\infty. After passing to a subsequence, Theorem 2.7 gives the decomposition in (2.26) with centers satisfying (2.27). Passing to a subsequence again, each of the following three bounded sequences converges to a finite number:

(4.4) M(χ2(x/tn)u(tn)),(χ2(x/tn)u(tn))L22,E(χ2(x/tn)u(tn)),\displaystyle M(\chi_{2}(x/t_{n})u(t_{n})),\qquad\|\nabla(\chi_{2}(x/t_{n})u(t_{n}))\|_{L^{2}}^{2},\qquad E(\chi_{2}(x/t_{n})u(t_{n})),

and, for every jj, the sequence xj,n/tnx_{j,n}/t_{n} either converges in d\mathbb{R}^{d} or tends to infinity:

(4.5) xj,ntn\displaystyle\frac{x_{j,n}}{t_{n}} vjX=d{}.\displaystyle\to v_{j}\in X=\mathbb{R}^{d}\cup\{\infty\}.

We will need the mass and energy decouplings:

(4.6) M(χ2(x/tn)u(tn))\displaystyle M(\chi_{2}(x/t_{n})u(t_{n})) =M(χ2(x/tn)eitnΔu+)+j=1JM(cjwj)+o(1),\displaystyle=M(\chi_{2}(x/t_{n})e^{it_{n}\Delta}u_{+})+\sum_{j=1}^{J}M(c_{j}w_{j})+o(1),
(4.7) (χ2(x/tn)u(tn))L22\displaystyle\|\nabla(\chi_{2}(x/t_{n})u(t_{n}))\|_{L^{2}}^{2} =(χ2(x/tn)eitnΔu+)L22+j=1J(cjwj)L22+o(1),\displaystyle=\|\nabla(\chi_{2}(x/t_{n})e^{it_{n}\Delta}u_{+})\|_{L^{2}}^{2}+\sum_{j=1}^{J}\|\nabla(c_{j}w_{j})\|_{L^{2}}^{2}+o(1),
(4.8) χ2(x/tn)u(tn)L2(d+1)/(d1)2(d+1)/(d1)\displaystyle\|\chi_{2}(x/t_{n})u(t_{n})\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)} =j=1JcjwjL2(d+1)/(d1)2(d+1)/(d1)+o(1).\displaystyle=\sum_{j=1}^{J}\|c_{j}w_{j}\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}+o(1).
(4.9) E(χ2(x/tn)u(tn))\displaystyle E(\chi_{2}(x/t_{n})u(t_{n})) =12(χ2(x/tn)eitnΔu+)L22+j=1JE(cjwj)+o(1).\displaystyle=\frac{1}{2}\|\nabla(\chi_{2}(x/t_{n})e^{it_{n}\Delta}u_{+})\|_{L^{2}}^{2}+\sum_{j=1}^{J}E(c_{j}w_{j})+o(1).

See e.g. [18, Lemma 3.3][7, Lemma 2.2, Lemma 2.3] for the proof. The nonlinear contribution of the radiation is absent from (4.8) and (4.9) because of the ‘dispersive estimate’ in Lemma 2.3.

We now use the threshold estimates from Proposition 3.1. For all sufficiently large nn, one has tnT0t_{n}\geq T_{0}, so (1.9)–(1.10) apply to f=χ2(x/tn)u(tn)f=\chi_{2}(x/t_{n})u(t_{n}). Applying (3.7) to this ff, and using (4.6) and (4.7), every fixed profile satisfies

(4.10) cjwjL2(cjwj)L2\displaystyle\|c_{j}w_{j}\|_{L^{2}}\|\nabla(c_{j}w_{j})\|_{L^{2}} 1δQL2QL2<QL2QL2.\displaystyle\leq\sqrt{1-\delta}\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}}<\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}}.

Applying (3.5) with f=cjwjf=c_{j}w_{j}, and using (4.10), gives

(4.11) E(cjwj)\displaystyle E(c_{j}w_{j}) 0.\displaystyle\geq 0.

Thus the quadratic radiation term and every profile-energy term displayed explicitly on the right side of (4.9) are nonnegative.

Suppose vjsuppχ1v_{j}\in\mathrm{supp}\chi_{1}. By (1.8), cj=χ2(vj)=1c_{j}=\chi_{2}(v_{j})=1.

(4.12) M(wj)\displaystyle M(w_{j}) limnM(χ2(x/tn)u(tn)),\displaystyle\leq\lim_{n\to\infty}M(\chi_{2}(x/t_{n})u(t_{n})),
(4.13) E(wj)\displaystyle E(w_{j}) limnE(χ2(x/tn)u(tn)).\displaystyle\leq\lim_{n\to\infty}E(\chi_{2}(x/t_{n})u(t_{n})).

The limits on the right-hand sides of (4.12)–(4.13) exist as explained in (4.4). The inequalities (4.12)–(4.13) follow from (4.6), (4.9), and (4.11). All quantities in (4.12)–(4.13) are nonnegative by (4.6), (4.9), and (4.11). Using (1.9), this gives

(4.14) M(wj)E(wj)\displaystyle M(w_{j})E(w_{j}) (1δ)M(Q)E(Q)<M(Q)E(Q).\displaystyle\leq(1-\delta)M(Q)E(Q)<M(Q)E(Q).

If wj0w_{j}\neq 0, then (3.11) gives

M(wj)E(wj)M(wj)E(wj)|P(wj)|22Θ=M(Q)E(Q),\displaystyle M(w_{j})E(w_{j})\geq M(w_{j})E(w_{j})-\frac{|P(w_{j})|^{2}}{2}\geq\Theta=M(Q)E(Q),

contradicting (4.14). Therefore

(4.15) vjsuppχ1\displaystyle v_{j}\in\mathrm{supp}\chi_{1} wj=0.\displaystyle\quad\Longrightarrow\quad w_{j}=0.

For every remaining profile, we know either vj=v_{j}=\infty or vjsuppχ1v_{j}\notin\mathrm{supp}\chi_{1}, then we have

(4.16) χ1(x/tn)τxj,nwjH1\displaystyle\|\chi_{1}(x/t_{n})\tau_{x_{j,n}}w_{j}\|_{H^{1}} 0.\displaystyle\to 0.

The oH1(1)o_{H^{1}}(1) term in (2.26) remains oH1(1)o_{H^{1}}(1) after multiplication by χ1(x/tn)\chi_{1}(x/t_{n}).

Subtracting the radiation term from (2.26), multiplying by χ1(x/tn)\chi_{1}(x/t_{n}), and using (4.15) and (4.16), we obtain

(4.17) χ1(x/tn)(u(tn)eitnΔu+)H1\displaystyle\|\chi_{1}(x/t_{n})(u(t_{n})-e^{it_{n}\Delta}u_{+})\|_{H^{1}} 0.\displaystyle\to 0.

Since the sequence tnt_{n}\to\infty is arbitrary, this proves (1.11). ∎

5. The atomic measures associated with the non-scattering part

In this section, we discuss the property of the measures (1.12) as tt\to\infty. In particular, we show that they can only concentrate at finitely many asymptotic directions. This gives a characterization of the non-scattering part of our solution u(t)u(t). Equivalently, this says that asymptotically the non-scattering part r(t)r(t) can only concentrate on finitely many directions.

To explain the intuition, suppose one indeed has a resolution of solitons. Then one should expect an asymptotic expansion of u(t)u(t) like

(5.1) u(t)Ctd/2u^+(x/2t)+i=1JQi,\displaystyle u(t)\sim Ct^{-d/2}\widehat{u}_{+}(x/2t)+\sum_{i=1}^{J}Q_{i},

where QiQ_{i} is the soliton component with velocity yiy_{i} and u^+\widehat{u}_{+} is the Fourier transform of the ‘free part’ u+u_{+}. See [15] for the discussion about the first term in the defocusing setting. If we consider the space of velocities parametrized by x/tx/t, then as tt\to\infty, the radiation part (i.e., the first term above) tends to a continuous measure with density |u^+|2|\widehat{u}_{+}|^{2}, while the remaining soliton components tend to a sum of Dirac measures. In addition, since, as we have seen, the soliton components and the radiation component decouple in the long-time asymptotics, if we consider the measure with density |r(t)|2=|u(t)eitΔu+|2|r(t)|^{2}=|u(t)-e^{it\Delta}u_{+}|^{2}, then we remove the continuous part directly (that is, the cross terms tend to zero) and are left with a sum of atomic (i.e., Dirac) measures. We will justify this intuition on the weak limit level.

Proposition 5.1.

Let JJ and u+u_{+} be as in Theorem 2.7, and define r(t)r(t) by (1.4). Let X=d{}X=\mathbb{R}^{d}\cup\{\infty\} be the one-point compactification of d\mathbb{R}^{d}, and let νt\nu_{t} be the residual mass measure defined in (1.12).

For every sequence tnt_{n}\to\infty, after passing to a subsequence, not relabeled, there are functions w1,,wJH1w_{1},\ldots,w_{J}\in H^{1}, centers xj,ndx_{j,n}\in\mathbb{R}^{d}, points yjXy_{j}\in X, and errors enH1e_{n}\in H^{1} such that

(5.2) r(tn)\displaystyle r(t_{n}) =j=1Jwj(xj,n)+en,enH10,\displaystyle=\sum_{j=1}^{J}w_{j}(\cdot-x_{j,n})+e_{n},\qquad\|e_{n}\|_{H^{1}}\to 0,
(5.3) xj,ntn\displaystyle\frac{x_{j,n}}{t_{n}} yj(1jJ),\displaystyle\to y_{j}\qquad(1\leq j\leq J),

and the centers satisfy (2.27). Along this subsequence,

(5.4) νtnj=1JM(wj)δyj.\displaystyle\nu_{t_{n}}\rightharpoonup\sum_{j=1}^{J}M(w_{j})\delta_{y_{j}}.

Here the functions wjw_{j} and the points yjy_{j} may depend on the original sequence.

Proof.

The expression of r(t)r(t) in (5.2) follows from (2.26). Since XX is compact33 3 Again, this is just a convenience of writing. One can separately consider sequences for which xj,n/tnx_{j,n}/t_{n} remains bounded and sequences for which xj,n/tnx_{j,n}/t_{n}\to\infty. and there are only finitely many centers, a further subsequence satisfies xj,n/tnyjx_{j,n}/t_{n}\to y_{j} in XX for every jj.

Let ΦC(X)\Phi\in C(X). The error ene_{n} contributes o(1)o(1) because

(5.5) |dΦ(x/tn)(|r(tn,x)|2|j=1Jwj(xxj,n)|2)𝑑x|\displaystyle\big|\int_{\mathbb{R}^{d}}\Phi(x/t_{n})\Big(|r(t_{n},x)|^{2}-\big|\sum_{j=1}^{J}w_{j}(x-x_{j,n})\big|^{2}\Big)\,dx\big| ΦL(X)enL2(2j=1JwjL2+enL2)0.\displaystyle\leq\|\Phi\|_{L^{\infty}(X)}\|e_{n}\|_{L^{2}}\Big(2\sum_{j=1}^{J}\|w_{j}\|_{L^{2}}+\|e_{n}\|_{L^{2}}\Big)\to 0.

For each jj, the change of variables x=xj,n+zx=x_{j,n}+z gives

(5.6) dΦ(x/tn)|wj(xxj,n)|2𝑑x\displaystyle\int_{\mathbb{R}^{d}}\Phi(x/t_{n})|w_{j}(x-x_{j,n})|^{2}\,dx =dΦ(xj,ntn+ztn)|wj(z)|2𝑑zΦ(yj)M(wj).\displaystyle=\int_{\mathbb{R}^{d}}\Phi\Big(\frac{x_{j,n}}{t_{n}}+\frac{z}{t_{n}}\Big)|w_{j}(z)|^{2}\,dz\to\Phi(y_{j})M(w_{j}).

Indeed, the argument of Φ\Phi converges to yjy_{j} in XX for every fixed zz, and dominated convergence applies.

For jkj\neq k, (2.27) gives

(5.7) dΦ(x/tn)wj(xxj,n)wk(xxk,n)¯𝑑x0.\displaystyle\int_{\mathbb{R}^{d}}\Phi(x/t_{n})w_{j}(x-x_{j,n})\overline{w_{k}(x-x_{k,n})}\,dx\to 0.

To see this, approximate wjw_{j} and wkw_{k} in L2L^{2} by compactly supported functions. Their translates have disjoint supports for all sufficiently large nn by (2.27), and Cauchy–Schwarz controls the approximation errors uniformly in nn.

Expanding |j=1Jwj(xxj,n)|2|\sum_{j=1}^{J}w_{j}(x-x_{j,n})|^{2} and using (5.5), (5.6), and (5.7), together with the definition of νt\nu_{t} in (1.12), yields

(5.8) XΦdνtn\displaystyle\int_{X}\Phi\,d\nu_{t_{n}} j=1JΦ(yj)M(wj)=XΦd(j=1JM(wj)δyj),\displaystyle\to\sum_{j=1}^{J}\Phi(y_{j})M(w_{j})=\int_{X}\Phi\,d\Big(\sum_{j=1}^{J}M(w_{j})\delta_{y_{j}}\Big),

which proves (5.4). ∎

Now we introduce the distance and the class of measures describing the limit we will consider. Fix a metric ρX\rho_{X} on XX such that the induced topology is the topology of the one-point compactification of d\mathbb{R}^{d}. For example, we can identify XX with 𝕊d\mathbb{S}^{d} using the stereographic projection and use the round sphere metric. Then we define

dBL(μ,σ)=sup{|XΦd(μσ)|:ΦL1,LipρX(Φ)1},\displaystyle d_{\mathrm{BL}}(\mu,\sigma)=\sup\big\{|\int_{X}\Phi\,d(\mu-\sigma)|:\|\Phi\|_{L^{\infty}}\leq 1,\ \mathrm{Lip}_{\rho_{X}}(\Phi)\leq 1\big\},

where BL stands for ‘bounded Lipschitz’. Let

dJ((μ,J),(σ,H))=dBL(μ,σ)+k=1ddBL(Jk,Hk).\displaystyle d_{J}((\mu,J),(\sigma,H))=d_{\mathrm{BL}}(\mu,\sigma)+\sum_{k=1}^{d}d_{\mathrm{BL}}(J_{k},H_{k}).

We use 𝒟\mathcal{D}_{\sharp} to denote the class of pairs (μ,J)(\mu,J) such that

(5.9) μ==1Laδy,J==1Lbδy,\displaystyle\mu=\sum_{\ell=1}^{L}a_{\ell}\delta_{y_{\ell}},\qquad J=\sum_{\ell=1}^{L}b_{\ell}\delta_{y_{\ell}},

where the yXy_{\ell}\in X are distinct, a>0a_{\ell}>0, bdb_{\ell}\in\mathbb{R}^{d}, and nn_{\ell}\in\mathbb{N}, and where the following conditions hold. Let m,e,p,Rm_{\sharp},e_{\sharp},p_{\sharp},R_{\sharp} be defined in (2.36), and define

(5.10) η=pm.\displaystyle\eta=\frac{p_{\sharp}}{m_{\sharp}}.

We require:

a=m,b=p,\displaystyle\sum_{\ell}a_{\ell}=m_{\sharp},\qquad\sum_{\ell}b_{\ell}=p_{\sharp},
(5.11) m=1Ln2Θa+m2=1La|baη|2R.\displaystyle m_{\sharp}\sum_{\ell=1}^{L}\frac{n_{\ell}^{2}\Theta}{a_{\ell}}+\frac{m_{\sharp}}{2}\sum_{\ell=1}^{L}a_{\ell}\big|\frac{b_{\ell}}{a_{\ell}}-\eta\big|^{2}\leq R_{\sharp}.
Remark 5.2.

In the proofs below, the reader can simply take n=1n_{\ell}=1 and collect all terms with the same yy_{\ell} together. We still include this parameter here because we believe that this is the conceptually correct formulation for future work treating the dynamics of soliton components at scales sublinear in tt, and we want to prove these technical lemmas at that level of generality. See (5.13) below for an example of the natural choice of nn_{\ell}.

Then we have the following result.

Proposition 5.3.

Let 𝒟\mathcal{D}_{\sharp} be the class defined by (5.9)–(5.11). If m>0m_{\sharp}>0, then R>0R_{\sharp}>0. With JtJ_{t} defined in (1.12), we have

(5.12) distdJ((νt,Jt),𝒟)0.\displaystyle\mathrm{dist}_{d_{J}}((\nu_{t},J_{t}),\mathcal{D}_{\sharp})\to 0.

For every profile decomposition (5.2)–(5.3) whose associated measures (1.12) converge to an atomic pair (μ,J)(\mu,J) of the form (5.9), nn_{\ell} may be chosen as

(5.13) n=#{j:wj0,xj,ntny}.\displaystyle n_{\ell}=\#\big\{j:w_{j}\neq 0,\ \frac{x_{j,n}}{t_{n}}\to y_{\ell}\big\}.

If N=nN=\sum_{\ell}n_{\ell}, then

(5.14) amn2ΘR,N2ΘR.\displaystyle a_{\ell}\geq\frac{m_{\sharp}n_{\ell}^{2}\Theta}{R_{\sharp}},\qquad N^{2}\Theta\leq R_{\sharp}.
Proof.

Take an arbitrary tnt_{n}\to\infty and apply Proposition 5.1 with profiles wjw_{j}.

By the same arguement referred for obtaining (4.6)-(4.9), applied to (5.2), without the cut-off, and the same argument applied to j(t)j(t) in place of |u|2|u|^{2}, together with Proposition 2.10, we have

(5.15) jM(wj)=m,jE(wj)=e,jP(wj)=p.\displaystyle\sum_{j}M(w_{j})=m_{\sharp},\qquad\sum_{j}E(w_{j})=e_{\sharp},\qquad\sum_{j}P(w_{j})=p_{\sharp}.

We already have the convergence of νtn\nu_{t_{n}} in (5.4). Now we prove

(5.16) JtnjP(wj)δyj.\displaystyle J_{t_{n}}\rightharpoonup\sum_{j}P(w_{j})\delta_{y_{j}}.

Indeed, the H1H^{1}-small remainder changes the current in L1L^{1} by o(1)o(1). The cross terms also have no contribution as nn\to\infty for the same reason as in the preceding decoupling argument. Finally, each diagonal term converges to give terms on the right-hand side of (5.16).

Define

νj=P(wj)M(wj),Rj=M(wj)E(wj)|P(wj)|22Θ>0,\displaystyle\nu_{j}=\frac{P(w_{j})}{M(w_{j})},\qquad R_{j}=M(w_{j})E(w_{j})-\frac{|P(w_{j})|^{2}}{2}\geq\Theta>0,

where the last inequality used Proposition 3.2. Using (5.15), a direct computation gives

(5.17) Rm=j(RjM(wj)+M(wj)2|νjη|2).\displaystyle\frac{R_{\sharp}}{m_{\sharp}}=\sum_{j}\Big(\frac{R_{j}}{M(w_{j})}+\frac{M(w_{j})}{2}|\nu_{j}-\eta|^{2}\Big).

This proves R>0R_{\sharp}>0.

Partition the indices jj with wj0w_{j}\neq 0 according to the limit of their normalized centers in (5.3), and write GG_{\ell} for the set of indices whose limit is yy_{\ell}. For each GG_{\ell}, set

n=|G|,a=jGM(wj),b=jGP(wj).\displaystyle n_{\ell}=|G_{\ell}|,\quad a_{\ell}=\sum_{j\in G_{\ell}}M(w_{j}),\quad b_{\ell}=\sum_{j\in G_{\ell}}P(w_{j}).

Using RjΘR_{j}\geq\Theta recalled above and Cauchy-Schwarz, we have

(5.18) jGRjM(wj)ΘjG1M(wj)n2Θa.\displaystyle\sum_{j\in G_{\ell}}\frac{R_{j}}{M(w_{j})}\geq\Theta\sum_{j\in G_{\ell}}\frac{1}{M(w_{j})}\geq\frac{n_{\ell}^{2}\Theta}{a_{\ell}}.

Noticing that (M(wj)(νjba),)(\sqrt{M(w_{j})}(\nu_{j}-\frac{b_{\ell}}{a_{\ell}}),...) and (,M(wj)(baη),)(...,\sqrt{M(w_{j})}(\frac{b_{\ell}}{a_{\ell}}-\eta),...), where components run through all jGj\in G_{\ell}, are orthogonal, we have

(5.19) jGM(wj)|νjη|2=jGM(wj)|νjba|2+a|baη|2a|baη|2.\displaystyle\sum_{j\in G_{\ell}}M(w_{j})|\nu_{j}-\eta|^{2}=\sum_{j\in G_{\ell}}M(w_{j})|\nu_{j}-\frac{b_{\ell}}{a_{\ell}}|^{2}+a_{\ell}\big|\frac{b_{\ell}}{a_{\ell}}-\eta\big|^{2}\geq a_{\ell}\big|\frac{b_{\ell}}{a_{\ell}}-\eta\big|^{2}.

Inserting (5.18) and (5.19) into (5.17) gives (5.11).

The mass and current measures have uniformly bounded total variation, so weak convergence on compact XX is convergence in dJd_{J}. Since our sequence tnt_{n} is arbitrary, we have proved (5.12).

Finally, each term in (5.11) is nonnegative, giving the bound in (5.14). Also

n2a(n)2a=N2m.\displaystyle\sum_{\ell}\frac{n_{\ell}^{2}}{a_{\ell}}\geq\frac{(\sum_{\ell}n_{\ell})^{2}}{\sum_{\ell}a_{\ell}}=\frac{N^{2}}{m_{\sharp}}.

This yields N2ΘRN^{2}\Theta\leq R_{\sharp}. ∎

We next record the compactness and connectedness of the set of subsequential limits.

Proposition 5.4.

Let Ω\Omega_{\sharp} be the set of all dJd_{J}-limits of (νtn,Jtn)(\nu_{t_{n}},J_{t_{n}}) along tnt_{n}\to\infty. Then 𝒟\mathcal{D}_{\sharp} is compact and Ω\Omega_{\sharp} is a nonempty compact connected subset of 𝒟\mathcal{D}_{\sharp}.

Remark 5.5.

In particular, this shows that if (νt,Jt)(\nu_{t},J_{t}) does not converge (i.e., if the sequential limits are not all the same), the failure cannot consist of shifts within a discrete set of configurations of the asymptotic atomic measure in (5.9).

Proof.

We prove 𝒟\mathcal{D}_{\sharp} is compact first. Given a sequence (μk,Jk)(\mu_{k},J_{k}) in 𝒟\mathcal{D}_{\sharp}, write

(5.20) μk==1Lka,kδy,k,Jk==1Lkb,kδy,k,\displaystyle\mu_{k}=\sum_{\ell=1}^{L_{k}}a_{\ell,k}\delta_{y_{\ell,k}},\qquad J_{k}=\sum_{\ell=1}^{L_{k}}b_{\ell,k}\delta_{y_{\ell,k}},

where the y,ky_{\ell,k} are distinct and the corresponding multiplicities are n,kn_{\ell,k}. By (5.14), we have

(5.21) Lk=1Lkn,kRΘ,a,kmn,k2ΘRmΘR.\displaystyle L_{k}\leq\sum_{\ell=1}^{L_{k}}n_{\ell,k}\leq\sqrt{\frac{R_{\sharp}}{\Theta}},\qquad a_{\ell,k}\geq\frac{m_{\sharp}n_{\ell,k}^{2}\Theta}{R_{\sharp}}\geq\frac{m_{\sharp}\Theta}{R_{\sharp}}.

Moreover, using (5.11), we have

(5.22) a,k|b,ka,kη|22Rm,\displaystyle a_{\ell,k}\big|\frac{b_{\ell,k}}{a_{\ell,k}}-\eta\big|^{2}\leq\frac{2R_{\sharp}}{m_{\sharp}},

and hence

(5.23) |b,k|a,k|η|+2Ra,kmm|η|+2R.\displaystyle|b_{\ell,k}|\leq a_{\ell,k}|\eta|+\sqrt{\frac{2R_{\sharp}a_{\ell,k}}{m_{\sharp}}}\leq m_{\sharp}|\eta|+\sqrt{2R_{\sharp}}.

So LkL_{k} and n,kn_{\ell,k} are uniformly bounded. After passing to a subsequence and relabeling, we may therefore assume

(5.24) Lk=L,n,k=n.\displaystyle L_{k}=L,\qquad n_{\ell,k}=n_{\ell}.

Similarly, since XX is compact, passing to a subsequence again, we may assume

(5.25) y,ky,a,ka,b,kb\displaystyle y_{\ell,k}\to y_{\ell},\qquad a_{\ell,k}\to a_{\ell},\qquad b_{\ell,k}\to b_{\ell}

for every 1L1\leq\ell\leq L. The lower bound for a,ka_{\ell,k} gives a>0a_{\ell}>0.

The limiting points yy_{\ell} need not be distinct. Suppose, for example,

(5.26) y1==ys.\displaystyle y_{1}=\cdots=y_{s}.

Replace these terms by one term at y1y_{1} with coefficient =1sa\sum_{\ell=1}^{s}a_{\ell} and =1sb\sum_{\ell=1}^{s}b_{\ell}, and multiplicity =1sn\sum_{\ell=1}^{s}n_{\ell}. Cauchy–Schwarz gives

(5.27) (=1sn)2=1sa=1sn2a,(=1sa)|=1sb=1saη|2=|=1sa(baη)|2=1sa=1sa|baη|2.\displaystyle\frac{(\sum_{\ell=1}^{s}n_{\ell})^{2}}{\sum_{\ell=1}^{s}a_{\ell}}\leq\sum_{\ell=1}^{s}\frac{n_{\ell}^{2}}{a_{\ell}},\quad\Big(\sum_{\ell=1}^{s}a_{\ell}\Big)\big|\frac{\sum_{\ell=1}^{s}b_{\ell}}{\sum_{\ell=1}^{s}a_{\ell}}-\eta\big|^{2}=\frac{\Big|\sum_{\ell=1}^{s}a_{\ell}\Big(\frac{b_{\ell}}{a_{\ell}}-\eta\Big)\Big|^{2}}{\sum_{\ell=1}^{s}a_{\ell}}\leq\sum_{\ell=1}^{s}a_{\ell}\big|\frac{b_{\ell}}{a_{\ell}}-\eta\big|^{2}.

Thus collecting the same limiting positions together does not increase either term in (5.11). Repeating this produces an expression that is a sum of delta functions at distinct limiting positions and still belongs to 𝒟\mathcal{D}_{\sharp}. The convergence of the positions, masses and momenta (i.e., those coefficients) gives convergence in dJd_{J}. Hence every sequence in 𝒟\mathcal{D}_{\sharp} has a subsequence converging in dJd_{J} to an element of 𝒟\mathcal{D}_{\sharp}, and we know that 𝒟\mathcal{D}_{\sharp} is compact.

The pair (νt,Jt)(\nu_{t},J_{t}) is continuous in tt with respect to dJd_{J}. This follows from H1H^{1}-continuity of r(t)r(t), dominated convergence for the self-similar multiplier, and

j(r(tn))j(r(t))L10when r(tn)r(t) in H1.\displaystyle\|j(r(t_{n}))-j(r(t))\|_{L^{1}}\to 0\quad\text{when }r(t_{n})\to r(t)\text{ in }H^{1}.

Its range is relatively compact because the mass and current total variations are uniformly bounded. For T1T\geq 1, consider

KT={(νt,Jt):tT}¯.\displaystyle K_{T}=\overline{\{(\nu_{t},J_{t}):t\geq T\}}.

Each KTK_{T} is nonempty, compact, and connected, and the family is nested (decreasing in TT). Its intersection is Ω\Omega_{\sharp}, so the latter is nonempty, compact, and connected. Since 𝒟\mathcal{D}_{\sharp} is compact, (5.12) shows Ω𝒟\Omega_{\sharp}\subset\mathcal{D}_{\sharp}. ∎

Next we consider the family of subsequential limits of our measures (νt,Jt)(\nu_{t},J_{t}) parametrized by a time translation parameter. Since the proof is quite involved, we briefly sketch the idea of the proof first. The goal is to derive the limiting behaviour of (νt,Jt)(\nu_{t},J_{t}). As Proposition 2.9 shows, roughly speaking, one can approximate |r|2|r|^{2} by |u|2|v|2|u|^{2}-|v|^{2}, and we will reduce our problem to investigating the measure with density |u|2|v|2|u|^{2}-|v|^{2}. Similarly, for the measure that originally has density j(r)j(r), we will approximate it by the measure with density j(u)j(v)j(u)-j(v). After integrating the time derivative of such a measure (see (5.31) below), we obtain an identity that is analogous to (5.29) at finite time. Then, roughly speaking, the goal is to ‘take the limit’ in this identity. See also DAD_{A} in (6.3) below for the motivation to consider such measures.

Another ingredient is that, in order to obtain a nice family of limits parametrized by the delayed time τ\tau, we prove the equicontinuity in τ\tau of the auxiliary sequence of measures above and then apply the Arzelà–Ascoli theorem. It is when we compute the time derivative of these measures that our PDEs enter, through the local conservation law in (5.35). Now we state our result.

Theorem 5.6.

Let r(t)r(t) be defined by (1.4), and suppose m:=M(u)M(u+)>0m_{\sharp}:=M(u)-M(u_{+})>0. Define νt\nu_{t} and JtJ_{t} by (1.12), and let dJd_{J} and 𝒟\mathcal{D}_{\sharp} be as in Proposition 5.3. For every sequence sns_{n}\to\infty, after potentially passing to a subsequence, there exists a family (μ(τ),J(τ))(\mu(\tau),J(\tau)) that is continuous in τ\tau with respect to dJd_{J}, such that, for every S>0S>0,

(5.28) sup|τ|SdJ((νesn+τ,Jesn+τ),(μ(τ),J(τ)))0.\displaystyle\sup_{|\tau|\leq S}d_{J}\big((\nu_{e^{s_{n}+\tau}},J_{e^{s_{n}+\tau}}),(\mu(\tau),J(\tau))\big)\to 0.

We call such a map a subsequential limit family. For every τ\tau, the pair belongs to 𝒟\mathcal{D}_{\sharp} and μ(τ)({})=0\mu(\tau)(\{\infty\})=0. Writing Jk(τ)J_{k}(\tau) for the kk-th scalar component of J(τ)J(\tau), we have

suppJk(τ)suppμ(τ)(1kd),\displaystyle\mathrm{supp}J_{k}(\tau)\subseteq\mathrm{supp}\mu(\tau)\qquad(1\leq k\leq d),

and suppμ(τ)\mathrm{supp}\mu(\tau) is a finite subset of d\mathbb{R}^{d}. For every real-valued ϕCc(d)\phi\in C_{c}^{\infty}(\mathbb{R}^{d}) and α<β\alpha<\beta, we have

(5.29) ϕ𝑑μ(β)ϕ𝑑μ(α)=αβ(2dϕ(y)dJ(τ)(y)dyϕ(y)𝑑μ(τ)(y))𝑑τ.\displaystyle\int\phi\,d\mu(\beta)-\int\phi\,d\mu(\alpha)=\int_{\alpha}^{\beta}\Big(2\int_{\mathbb{R}^{d}}\nabla\phi(y)\cdot dJ(\tau)(y)-\int_{\mathbb{R}^{d}}y\cdot\nabla\phi(y)\,d\mu(\tau)(y)\Big)d\tau.
Proof.

Let v(t)=eitΔu+v(t)=e^{it\Delta}u_{+} be as in (1.4), and let j(f)=(f¯f)j(f)=\Im(\overline{f}\nabla f) as in (2.29). Since we have uniform H1H^{1} bounds for uu, vv, and rr, Cauchy–Schwarz gives

(5.30) supt(|u(t)|2|v(t)|2L1+j(u(t))j(v(t))L1)<.\displaystyle\sup_{t}\Big(\||u(t)|^{2}-|v(t)|^{2}\|_{L^{1}}+\|j(u(t))-j(v(t))\|_{L^{1}}\Big)<\infty.

Set tn,τ=esn+τt_{n,\tau}=e^{s_{n}+\tau}. For a real-valued ϕCc(d)\phi\in C_{c}^{\infty}(\mathbb{R}^{d}), define

(5.31) Aϕ(s)=dϕ(x/es)(|u|2|v|2)(es,x)𝑑x.\displaystyle A_{\phi}(s)=\int_{\mathbb{R}^{d}}\phi(x/e^{s})(|u|^{2}-|v|^{2})(e^{s},x)\,dx.

To compute the time derivative of AϕA_{\phi}, we use the mass-current conservation laws for uu and vv (due to (1.1) and the fact that vv satisfies the linear Schrödinger equation):

t|u|2+2j(u)=0,t|v|2+2j(v)=0.\displaystyle\partial_{t}|u|^{2}+2\nabla\cdot j(u)=0,\quad\partial_{t}|v|^{2}+2\nabla\cdot j(v)=0.

This gives t(|u|2|v|2)=2(j(u)j(v))\partial_{t}(|u|^{2}-|v|^{2})=-2\nabla\cdot(j(u)-j(v)), and we have

(5.32) Aϕ(s)=dϕ(x/es)(2(j(u)j(v))(es,x)(x/es)(|u|2|v|2)(es,x))𝑑x.\displaystyle A_{\phi}^{\prime}(s)=\int_{\mathbb{R}^{d}}\nabla\phi(x/e^{s})\cdot\Big(2(j(u)-j(v))(e^{s},x)-(x/e^{s})(|u|^{2}-|v|^{2})(e^{s},x)\Big)\,dx.

In particular, for α<β\alpha<\beta, integrating (5.32) from sn+αs_{n}+\alpha to sn+βs_{n}+\beta gives

(5.33) Aϕ(sn+β)Aϕ(sn+α)=2αβdϕ(x/tn,τ)(j(u)j(v))(tn,τ,x)𝑑x𝑑ταβd(x/tn,τ)ϕ(x/tn,τ)(|u|2|v|2)(tn,τ,x)dxdτ.\displaystyle\begin{split}A_{\phi}(s_{n}+\beta)-A_{\phi}(s_{n}+\alpha)&=2\int_{\alpha}^{\beta}\int_{\mathbb{R}^{d}}\nabla\phi(x/t_{n,\tau})\cdot(j(u)-j(v))(t_{n,\tau},x)\,dx\,d\tau\\ &\quad-\int_{\alpha}^{\beta}\int_{\mathbb{R}^{d}}(x/t_{n,\tau})\cdot\nabla\phi(x/t_{n,\tau})(|u|^{2}-|v|^{2})(t_{n,\tau},x)\,dx\,d\tau.\end{split}

This is the finite-time equality that we will pass to the limit. The right-hand side of (5.32) is uniformly bounded by (5.30).

We also need to consider the measure with j(u)j(v)j(u)-j(v) being its density. This is for both the joint convergence in (5.28) and for the current term in (5.33). For 1kd1\leq k\leq d, set

(5.34) Bϕ,k(s)=dϕ(x/es)(jk(u)jk(v))(es,x)𝑑x.\displaystyle B_{\phi,k}(s)=\int_{\mathbb{R}^{d}}\phi(x/e^{s})\big(j_{k}(u)-j_{k}(v)\big)(e^{s},x)\,dx.

The local momentum laws for uu and vv reads

(5.35) tjk(u)+m=1dm(2(ku¯mu)2d+1δkm|u|2(d+1)/(d1))12kΔ|u|2=0,tjk(v)+m=1dm(2(kv¯mv))12kΔ|v|2=0.\displaystyle\begin{split}\partial_{t}j_{k}(u)&+\sum_{m=1}^{d}\partial_{m}\Big(2\Re(\partial_{k}\overline{u}\,\partial_{m}u)-\frac{2}{d+1}\delta_{km}|u|^{2(d+1)/(d-1)}\Big)-\frac{1}{2}\partial_{k}\Delta|u|^{2}=0,\\ \partial_{t}j_{k}(v)&+\sum_{m=1}^{d}\partial_{m}\Big(2\Re(\partial_{k}\overline{v}\,\partial_{m}v)\Big)-\frac{1}{2}\partial_{k}\Delta|v|^{2}=0.\end{split}

Using (5.35), we have

(5.36) Bϕ,k(s)=d(x/es)ϕ(x/es)(jk(u)jk(v))(es,x)dx+m=1ddmϕ(x/es)(2(ku¯mukv¯mv)2d+1δkm|u|2(d+1)/(d1))(es,x)dx12e2sd(kΔϕ)(x/es)(|u|2|v|2)(es,x)dx.\displaystyle\begin{split}B_{\phi,k}^{\prime}(s)=&-\int_{\mathbb{R}^{d}}(x/e^{s})\cdot\nabla\phi(x/e^{s})\big(j_{k}(u)-j_{k}(v)\big)(e^{s},x)\,dx\\ &+\sum_{m=1}^{d}\int_{\mathbb{R}^{d}}\partial_{m}\phi(x/e^{s})\Big(2\Re(\partial_{k}\overline{u}\,\partial_{m}u-\partial_{k}\overline{v}\,\partial_{m}v)-\frac{2}{d+1}\delta_{km}|u|^{2(d+1)/(d-1)}\Big)(e^{s},x)\,dx\\ &-\frac{1}{2e^{2s}}\int_{\mathbb{R}^{d}}(\partial_{k}\Delta\phi)(x/e^{s})(|u|^{2}-|v|^{2})(e^{s},x)\,dx.\end{split}

The first line of (5.36) is uniformly bounded by (5.30). The quadratic terms in the second line are controlled by the uniform H1H^{1} bounds, and for the nonlinear term we use the embedding H1L2(d+1)/(d1)H^{1}\hookrightarrow L^{2(d+1)/(d-1)} from (2.10). The last line is O(e2s)O(e^{-2s}) by (5.30).

Thus, for every fixed S>0S>0, the translated functions Aϕ(sn+τ)A_{\phi}(s_{n}+\tau) and Bϕ,k(sn+τ)B_{\phi,k}(s_{n}+\tau) are Lipschitz with a uniform Lipschitz constant for |τ|S|\tau|\leq S and all sufficiently large nn.

Now we explain our approximation step. If |τ|S|\tau|\leq S, then tn,τesnSt_{n,\tau}\geq e^{s_{n}-S}\to\infty. Consequently, (2.31) and (2.32) show that, for every ΦC(X,)\Phi\in C(X;\mathbb{R}), we have

(5.37) sup|τ|S|dΦ(x/tn,τ)(|u|2|v|2|r|2)(tn,τ,x)𝑑x|0,\displaystyle\sup_{|\tau|\leq S}\Big|\int_{\mathbb{R}^{d}}\Phi(x/t_{n,\tau})\big(|u|^{2}-|v|^{2}-|r|^{2}\big)(t_{n,\tau},x)\,dx\Big|\to 0,
(5.38) sup|τ|S|dΦ(x/tn,τ)(j(u)j(v)j(r))(tn,τ,x)𝑑x|0.\displaystyle\sup_{|\tau|\leq S}\Big|\int_{\mathbb{R}^{d}}\Phi(x/t_{n,\tau})\big(j(u)-j(v)-j(r)\big)(t_{n,\tau},x)\,dx\Big|\to 0.

We now prove (5.28). Choose a countable dense family {Φm}m1C(X,)\{\Phi_{m}\}_{m\geq 1}\subset C(X;\mathbb{R}) containing the constant function 11, such that every Φm\Phi_{m} restricts to a constant plus a function in Cc(d)C_{c}^{\infty}(\mathbb{R}^{d}), and set

𝒱=span{Φm:m1}.\displaystyle\mathcal{V}=\operatorname{span}_{\mathbb{R}}\{\Phi_{m}:m\geq 1\}.

Write Φm=cm+ϕm\Phi_{m}=c_{m}+\phi_{m}, where cmc_{m}\in\mathbb{R} and ϕmCc(d)\phi_{m}\in C_{c}^{\infty}(\mathbb{R}^{d}). Conservation of mass and momentum gives

d(|u(t)|2|v(t)|2)𝑑x=m,d(j(u(t))j(v(t)))𝑑x=p,\displaystyle\int_{\mathbb{R}^{d}}(|u(t)|^{2}-|v(t)|^{2})\,dx=m_{\sharp},\qquad\int_{\mathbb{R}^{d}}\big(j(u(t))-j(v(t))\big)\,dx=p_{\sharp},

and hence

dΦm(x/tn,τ)(|u|2|v|2)(tn,τ,x)𝑑x=cmm+Aϕm(sn+τ),\displaystyle\int_{\mathbb{R}^{d}}\Phi_{m}(x/t_{n,\tau})(|u|^{2}-|v|^{2})(t_{n,\tau},x)\,dx=c_{m}m_{\sharp}+A_{\phi_{m}}(s_{n}+\tau),
dΦm(x/tn,τ)(jk(u)jk(v))(tn,τ,x)𝑑x=cm(p)k+Bϕm,k(sn+τ).\displaystyle\int_{\mathbb{R}^{d}}\Phi_{m}(x/t_{n,\tau})\big(j_{k}(u)-j_{k}(v)\big)(t_{n,\tau},x)\,dx=c_{m}(p_{\sharp})_{k}+B_{\phi_{m},k}(s_{n}+\tau).

These are uniformly bounded by our uniform H1H^{1} property of u,vu,v, and they are Lipschitz in τ\tau over any fixed compact interval. The Arzelà–Ascoli theorem and a diagonal extraction over mm, the finitely many current components kk, and the bounds |τ|S|\tau|\leq S with SS\in\mathbb{N} therefore give one subsequence on which all these pairings converge locally uniformly in τ\tau. By linearity, the same holds for every Φ𝒱\Phi\in\mathcal{V}. The decoupling estimates (5.37) and (5.38) show that the corresponding pairings of (νtn,τ,Jtn,τ)(\nu_{t_{n,\tau}},J_{t_{n,\tau}}) have the same limits.

Let TV\|\cdot\|_{\mathrm{TV}} denote total variation. By (1.12),

(5.39) νtTV=r(t)L22,JtTVr(t)L2r(t)L2.\displaystyle\|\nu_{t}\|_{\mathrm{TV}}=\|r(t)\|_{L^{2}}^{2},\qquad\|J_{t}\|_{\mathrm{TV}}\leq\|r(t)\|_{L^{2}}\|\nabla r(t)\|_{L^{2}}.

Thus each scalar limit at each τ\tau defines a linear functional on 𝒱\mathcal{V} bounded by CΦLC\|\Phi\|_{L^{\infty}}, with CC independent of τ\tau. Since 𝒱\mathcal{V} is dense in C(X,)C(X;\mathbb{R}), each such functional has a unique extension to C(X,)C(X;\mathbb{R}) and is represented by one of the measures μ(τ),J1(τ),,Jd(τ)\mu(\tau),J_{1}(\tau),\ldots,J_{d}(\tau), whose total variation is at most CC.

Now we upgrade our results from 𝒱\mathcal{V} to general ΦC(X,)\Phi\in C(X;\mathbb{R}). Indeed, for Ψ𝒱\Psi\in\mathcal{V},

(5.40) |XΦd(νtn,τμ(τ))|\displaystyle\big|\int_{X}\Phi\,d(\nu_{t_{n,\tau}}-\mu(\tau))\big|\leq |XΨd(νtn,τμ(τ))|+2CΦΨL.\displaystyle\big|\int_{X}\Psi\,d(\nu_{t_{n,\tau}}-\mu(\tau))\big|+2C\|\Phi-\Psi\|_{L^{\infty}}.

The same estimate holds for each current component. Approximating Φ\Phi uniformly by an element of 𝒱\mathcal{V} proves locally uniform convergence for every fixed continuous test function.

Fix S>0S>0 and let

={ΦC(X):ΦL1,LipρX(Φ)1}.\displaystyle\mathcal{B}=\{\Phi\in C(X):\|\Phi\|_{L^{\infty}}\leq 1,\ \mathrm{Lip}_{\rho_{X}}(\Phi)\leq 1\}.

Since XX is compact, \mathcal{B} is compact in C(X)C(X) by Arzelà–Ascoli theorem. Given ϵ>0\epsilon>0, choose a finite ϵ\epsilon-net Ψ1,,ΨN\Psi_{1},\ldots,\Psi_{N} for \mathcal{B}. For |τ|S|\tau|\leq S, (5.40) gives

(5.41) dJ((νtn,τ,Jtn,τ),(μ(τ),J(τ)))\displaystyle d_{J}\big((\nu_{t_{n,\tau}},J_{t_{n,\tau}}),(\mu(\tau),J(\tau))\big)\leq max1iN|XΨid(νtn,τμ(τ))|\displaystyle\max_{1\leq i\leq N}\big|\int_{X}\Psi_{i}\,d(\nu_{t_{n,\tau}}-\mu(\tau))\big|
+\displaystyle+ k=1dmax1iN|XΨid(Jtn,τ,kJk(τ))|+2(d+1)Cϵ.\displaystyle\sum_{k=1}^{d}\max_{1\leq i\leq N}\big|\int_{X}\Psi_{i}\,d(J_{t_{n,\tau},k}-J_{k}(\tau))\big|+2(d+1)C\epsilon.

The first two terms on the right tend to zero uniformly for |τ|S|\tau|\leq S, because the net is finite. Hence

lim supnsup|τ|SdJ((νtn,τ,Jtn,τ),(μ(τ),J(τ)))2(d+1)Cϵ.\displaystyle\limsup_{n\to\infty}\sup_{|\tau|\leq S}d_{J}\big((\nu_{t_{n,\tau}},J_{t_{n,\tau}}),(\mu(\tau),J(\tau))\big)\leq 2(d+1)C\epsilon.

Letting ϵ0\epsilon\to 0 proves (5.28).

For each nn, (νtn,τ,Jtn,τ)(\nu_{t_{n,\tau}},J_{t_{n,\tau}}) is continuous in τ\tau with respect to dJd_{J}. Its locally uniform limit (μ(τ),J(τ))(\mu(\tau),J(\tau)) is therefore also continuous in τ\tau with respect to dJd_{J}. Now for each fixed τ\tau, tn,τt_{n,\tau}\to\infty. So we have

(μ(τ),J(τ))𝒟\displaystyle(\mu(\tau),J(\tau))\in\mathcal{D}_{\sharp}

by Proposition 5.3 and the compactness of 𝒟\mathcal{D}_{\sharp} from Proposition 5.4. By the definition of 𝒟\mathcal{D}_{\sharp}, the measure μ(τ)\mu(\tau) has finite support in XX and

suppJk(τ)suppμ(τ)(1kd).\displaystyle\mathrm{supp}J_{k}(\tau)\subseteq\mathrm{supp}\mu(\tau)\qquad(1\leq k\leq d).

By Lemma 2.5, we have μ(τ)({})=0\mu(\tau)(\{\infty\})=0, or equivalently suppμ(τ)d\mathrm{supp}\mu(\tau)\subset\mathbb{R}^{d} is compact.

We finally return to the exact finite-time identity (5.33). Fix α<β\alpha<\beta. For γ{α,β}\gamma\in\{\alpha,\beta\}, (5.28) and (5.37) give

Aϕ(sn+γ)Xϕ𝑑μ(γ).\displaystyle A_{\phi}(s_{n}+\gamma)\to\int_{X}\phi\,d\mu(\gamma).

Moreover, uniformly for ατβ\alpha\leq\tau\leq\beta, (5.28), (5.37), and (5.38) give

dϕ(x/tn,τ)(j(u)j(v))(tn,τ,x)𝑑xXϕ(y)dJ(τ)(y),\displaystyle\int_{\mathbb{R}^{d}}\nabla\phi(x/t_{n,\tau})\cdot(j(u)-j(v))(t_{n,\tau},x)\,dx\to\int_{X}\nabla\phi(y)\cdot dJ(\tau)(y),
d(x/tn,τ)ϕ(x/tn,τ)(|u|2|v|2)(tn,τ,x)𝑑xXyϕ(y)𝑑μ(τ)(y).\displaystyle\int_{\mathbb{R}^{d}}(x/t_{n,\tau})\cdot\nabla\phi(x/t_{n,\tau})(|u|^{2}-|v|^{2})(t_{n,\tau},x)\,dx\to\int_{X}y\cdot\nabla\phi(y)\,d\mu(\tau)(y).

The functions ϕ(y)\nabla\phi(y) and yϕ(y)y\cdot\nabla\phi(y) extend continuously to XX, with value zero at \infty. We may therefore pass to the limit in (5.33). Since μ(τ)\mu(\tau) and J(τ)J(\tau) have no support at \infty, the resulting identity is exactly (5.29). ∎

The identity (5.29) in Theorem 5.6 does not control the momentum and this is the issue that we will handle next using a similar argument. Throughout the rest of this section we fix a subsequential limit family (μ(τ),J(τ))(\mu(\tau),J(\tau)) obtained along sns_{n}\to\infty in Theorem 5.6, write tn,τ=esn+τt_{n,\tau}=e^{s_{n}+\tau}, and let v(t)v(t) be as in (1.4). Finally, for 1k,md1\leq k,m\leq d we define a (signed) measure Σn,km(τ)\Sigma_{n,km}(\tau) on XX by

Xϕ(y)dΣn,km(τ)(y)=\displaystyle\int_{X}\phi(y)\,d\Sigma_{n,km}(\tau)(y)= dϕ(x/tn,τ)(2(ku¯mukv¯mv)2d+1δkm|u|2(d+1)/(d1))(tn,τ,x)𝑑x\displaystyle\int_{\mathbb{R}^{d}}\phi(x/t_{n,\tau})\Big(2\Re(\partial_{k}\overline{u}\,\partial_{m}u-\partial_{k}\overline{v}\,\partial_{m}v)-\frac{2}{d+1}\delta_{km}|u|^{2(d+1)/(d-1)}\Big)(t_{n,\tau},x)\,dx

for ϕC(X)\phi\in C(X).

Proposition 5.7.

After passing to a further subsequence, there are signed measures Σkm(τ)\Sigma_{km}(\tau) on XX, uniformly bounded in total variation, such that

(5.42) XΦ(τ,y)dΣn,km(τ)(y)𝑑τ\displaystyle\int_{\mathbb{R}}\int_{X}\Phi(\tau,y)\,d\Sigma_{n,km}(\tau)(y)\,d\tau XΦ(τ,y)dΣkm(τ)(y)𝑑τ\displaystyle\to\int_{\mathbb{R}}\int_{X}\Phi(\tau,y)\,d\Sigma_{km}(\tau)(y)\,d\tau

whenever Φ(τ,)\Phi(\tau,\cdot) is integrable in τ\tau with values in C(X)C(X) and vanishes outside a bounded interval. For every χCc(d)\chi\in C_{c}^{\infty}(\mathbb{R}^{d}), α<β\alpha<\beta, and 1kd1\leq k\leq d,

(5.43) dχdJk(β)dχdJk(α)\displaystyle\int_{\mathbb{R}^{d}}\chi\,dJ_{k}(\beta)-\int_{\mathbb{R}^{d}}\chi\,dJ_{k}(\alpha) =αβ(dyχ(y)dJk(τ)(y)+m=1ddmχ(y)dΣkm(τ)(y))dτ.\displaystyle=\int_{\alpha}^{\beta}\Big(-\int_{\mathbb{R}^{d}}y\cdot\nabla\chi(y)\,dJ_{k}(\tau)(y)+\sum_{m=1}^{d}\int_{\mathbb{R}^{d}}\partial_{m}\chi(y)\,d\Sigma_{km}(\tau)(y)\Big)d\tau.

Let II\subset\mathbb{R} be a compact interval, and let KdK\subset\mathbb{R}^{d} be compact. Suppose that

(5.44) suppμ(τ)K,τI,\displaystyle\mathrm{supp}\mu(\tau)\subset K,\quad\forall\tau\in I,

Then 44 4 Rigorously speaking, Σkm\Sigma_{km} is only defined a.e. for τ\tau. But it will only be involved in integrals and we ignore this issue below.

(5.45) suppΣkm(τ)Kfor almost every τI.\displaystyle\mathrm{supp}\,\Sigma_{km}(\tau)\subset K\qquad\text{for almost every }\tau\in I.
Proof.

The uniform H1H^{1} bound in Theorem 2.7, the H1H^{1}-unitarity of eitΔe^{it\Delta}, and (2.10) give

Σn,km(τ)TVu(esn+τ)L22+v(esn+τ)L22+u(esn+τ)L2(d+1)/(d1)2(d+1)/(d1)1.\displaystyle\|\Sigma_{n,km}(\tau)\|_{\mathrm{TV}}\lesssim\|\nabla u(e^{s_{n}+\tau})\|_{L^{2}}^{2}+\|\nabla v(e^{s_{n}+\tau})\|_{L^{2}}^{2}+\|u(e^{s_{n}+\tau})\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}\lesssim 1.

On each bounded interval I0I_{0}, integration against Σn,km(τ)\Sigma_{n,km}(\tau) defines a uniformly bounded functional on L1(I0,C(X))L^{1}(I_{0};C(X)). So we have a subsequential weak-star limit, and the representation of the dual gives signed measures Σkm(τ)\Sigma_{km}(\tau), defined for almost every τ\tau, whose pairings with functions in C(X)C(X) are measurable and whose total variations have the same uniform bound. A diagonal extraction over the components and over an increasing sequence of intervals gives (5.42).

For χCc(d)\chi\in C_{c}^{\infty}(\mathbb{R}^{d}), use the current pairing Bχ,kB_{\chi,k} defined in (5.34). By the chain rule and (5.36), the function Bχ,k(sn+τ)B_{\chi,k}(s_{n}+\tau) is locally absolutely continuous in τ\tau and has, for almost every τ\tau, the derivative given there with ϕ=χ\phi=\chi and s=sn+τs=s_{n}+\tau. By the definition of Σn,km\Sigma_{n,km}, the stress term in that identity is

m=1dXmχ(y)dΣn,km(τ)(y).\displaystyle\sum_{m=1}^{d}\int_{X}\partial_{m}\chi(y)\,d\Sigma_{n,km}(\tau)(y).

By (5.28) and (5.38), the endpoint values Bχ,k(sn+γ)B_{\chi,k}(s_{n}+\gamma), γ{α,β}\gamma\in\{\alpha,\beta\}, converge to the corresponding pairings with JkJ_{k}. The first term on the right-hand side of (5.36), evaluated with ϕ=χ\phi=\chi and s=sn+τs=s_{n}+\tau, converges uniformly on bounded τ\tau-intervals to

dyχ(y)dJk(τ)(y).\displaystyle-\int_{\mathbb{R}^{d}}y\cdot\nabla\chi(y)\,dJ_{k}(\tau)(y).

The last term in (5.36), with the same substitution, tends uniformly to zero by (5.30), since esn+τe^{s_{n}+\tau}\to\infty uniformly on bounded τ\tau-intervals. Integrating (5.36) from sn+αs_{n}+\alpha to sn+βs_{n}+\beta, changing variables s=sn+τs=s_{n}+\tau, and using (5.42) with Φ(τ,y)=𝟏{ατβ}(τ)mχ(y)\Phi(\tau,y)=\mathbf{1}_{\{\alpha\leq\tau\leq\beta\}}(\tau)\partial_{m}\chi(y), we obtain (5.43).

It remains to prove (5.45). Take χ~C(I×d)\tilde{\chi}\in C^{\infty}(I\times\mathbb{R}^{d}) that has bounded derivatives, and is constant in yy outside a fixed compact set, and vanishes on a neighborhood of KK. Then we will show

(5.46) supτIχ~(τ,x/esn+τ)r(esn+τ)H10.\displaystyle\sup_{\tau\in I}\|\tilde{\chi}(\tau,x/e^{s_{n}+\tau})r(e^{s_{n}+\tau})\|_{H^{1}}\to 0.

If (5.46) is not true, we can choose τnτI\tau_{n}\to\tau_{*}\in I along which the norm stays bounded away from zero. Apply Theorem 2.7 at the times esn+τne^{s_{n}+\tau_{n}} and pass to a subsequence on which all normalized profile centers converge in XX. The convergence (5.28), the dJd_{J}-continuity in Theorem 5.6, and Proposition 5.1 show that the limiting center of every nonzero profile belongs to suppμ(τ)\mathrm{supp}\mu(\tau_{*}). Multiplication by χ~(τn,x/esn+τn)\tilde{\chi}(\tau_{n},x/e^{s_{n}+\tau_{n}}) therefore sends each translated profile to zero in H1H^{1} by dominated convergence. The same is true of the oH1(1)o_{H^{1}}(1) remainder in Theorem 2.7; the term created by differentiating the cutoff has the additional factor e(sn+τn)e^{-(s_{n}+\tau_{n})}. This contradiction proves (5.46).

Let ΦC(I×X)\Phi\in C(I\times X) be supported away from KK. Choose χ~\tilde{\chi} satisfying the hypotheses of (5.46) and equal to one on a neighborhood of suppΦ\mathrm{supp}\Phi. The product rule gives, uniformly on II,

(5.47) χ~(τ,x/esn+τ)r(esn+τ)L2\displaystyle\|\tilde{\chi}(\tau,x/e^{s_{n}+\tau})\nabla r(e^{s_{n}+\tau})\|_{L^{2}} (χ~(τ,x/esn+τ)r(esn+τ))L2+e(sn+τ)yχ~Lr(esn+τ)L20.\displaystyle\leq\|\nabla\big(\tilde{\chi}(\tau,x/e^{s_{n}+\tau})r(e^{s_{n}+\tau})\big)\|_{L^{2}}+e^{-(s_{n}+\tau)}\|\nabla_{y}\tilde{\chi}\|_{L^{\infty}}\|r(e^{s_{n}+\tau})\|_{L^{2}}\to 0.

Since u=v+ru=v+r, (5.47) and the uniform H1H^{1} bounds make the quadratic part of

XΦ(τ,y)dΣn,km(τ)(y)\displaystyle\int_{X}\Phi(\tau,y)\,d\Sigma_{n,km}(\tau)(y)

tend to zero uniformly for τI\tau\in I. Its nonlinear part also tends to zero uniformly, because (2.10) controls the localized rr term and Lemma 2.3 gives

supτIv(esn+τ)L2(d+1)/(d1)0.\displaystyle\sup_{\tau\in I}\|v(e^{s_{n}+\tau})\|_{L^{2(d+1)/(d-1)}}\to 0.

After integration in τ\tau, (5.42) shows that the signed measure Σkm(τ)dτ\Sigma_{km}(\tau)\,d\tau on I×XI\times X annihilates every continuous test supported in (I×X)(I×K)(I\times X)\setminus(I\times K), which proves (5.45). ∎

Now we show that the asymptotic velocities at which (νt,Jt)(\nu_{t},J_{t}) are concentrated cannot be arbitrarily large. Conceptually this is similar to Lemma 2.5.

Proposition 5.8.

There is a constant Ru<R_{u}<\infty, depending only on uu, such that every subsequential limit family satisfies

(5.48) suppμ(τ){yd:|y|Ru}(τ).\displaystyle\mathrm{supp}\mu(\tau)\subseteq\{y\in\mathbb{R}^{d}:|y|\leq R_{u}\}\qquad(\tau\in\mathbb{R}).

In addition, (5.48) also holds with suppμ(τ)\mathrm{supp}\mu(\tau) replaced by suppJk(τ)\mathrm{supp}J_{k}(\tau), 1kd1\leq k\leq d.

Proof.

Fix τ\tau and write the atomic masses of μ(τ)\mu(\tau) as aa_{\ell}. We prove (5.48). Since every (μ(τ),J(τ))(\mu(\tau),J(\tau)) belongs to 𝒟\mathcal{D}_{\sharp} by Theorem 5.6, (5.14) gives

aa:=mΘR>0\displaystyle a_{\ell}\geq a_{*}:=\frac{m_{\sharp}\Theta}{R_{\sharp}}>0

for every atom of every subsequential limit family. Let ϑ\vartheta be the cutoff in Lemma 2.5. Choose R<R<\infty so large that

supt1ϑ(x/(Rt))r(t)L22<a2.\displaystyle\sup_{t\geq 1}\|\vartheta(x/(Rt))r(t)\|_{L^{2}}^{2}<\frac{a_{*}}{2}.

Regard ϑ(y/R)2\vartheta(y/R)^{2} as a function on XX that is continuous in yy, with value one at \infty. For each fixed τ\tau, (5.28) gives

Xϑ(y/R)2𝑑μ(τ)(y)=limnϑ(x/(Rtn,τ))r(tn,τ)L22a2.\displaystyle\int_{X}\vartheta(y/R)^{2}\,d\mu(\tau)(y)=\lim_{n\to\infty}\|\vartheta(x/(Rt_{n,\tau}))r(t_{n,\tau})\|_{L^{2}}^{2}\leq\frac{a_{*}}{2}.

Because ϑ(y/R)=1\vartheta(y/R)=1 when |y|2R|y|\geq 2R, no atom of mass at least aa_{*} can lie there. Thus (5.48) holds with Ru=2RR_{u}=2R. The assertion for Jk(τ)J_{k}(\tau) follows from suppJk(τ)suppμ(τ)\mathrm{supp}J_{k}(\tau)\subseteq\mathrm{supp}\mu(\tau) in Theorem 5.6.

We record the atomic form of the limiting mass measure.

Proposition 5.9.

Let the setting and notation be as in Theorem 5.6. Assume that νtμ\nu_{t}\rightharpoonup\mu_{\infty} weakly on XX as tt\to\infty. Then there are L1L\geq 1, distinct y1,,yLdy_{1},\ldots,y_{L}\in\mathbb{R}^{d}, and numbers a>0a_{\ell}>0 such that

μ==1Laδy.\displaystyle\mu_{\infty}=\sum_{\ell=1}^{L}a_{\ell}\delta_{y_{\ell}}.
Proof.

Let tnt_{n}\to\infty be arbitrary and set sn=logtns_{n}=\log t_{n}. After passing to a subsequence, Theorem 5.6 gives a subsequential limit family (μ(τ),J(τ))(\mu(\tau),J(\tau)). At τ=0\tau=0, the assumed convergence of νt\nu_{t} and (5.28) give μ(0)=μ\mu(0)=\mu_{\infty}. In addition, Theorem 5.6 also gives that suppμ(0)\mathrm{supp}\mu(0) is in a finite subset of d\mathbb{R}^{d}, which gives the stated representation. ∎

6. Uniform in time estimates

We prove some upgrading estimates in this section. The upgrading here is from estimates at a sequence of times to estimates on a sequence of time intervals. We begin with some computations preparing for our ‘propagation estimates’. Here we take uCc(d,)u\in C_{c}^{\infty}(\mathbb{R}^{d};\mathbb{C}); it will eventually be a truncated version of our function.

Lemma 6.1.

Let bCc(d,d)b\in C_{c}^{\infty}(\mathbb{R}^{d};\mathbb{R}^{d}), cCc(d,)c\in C_{c}^{\infty}(\mathbb{R}^{d};\mathbb{R}), and define a(y,ξ)=b(y)ξ+c(y)a(y,\xi)=b(y)\cdot\xi+c(y). For t>0t>0, set At=Opw(a(x/t,ξ))A_{t}=\mathrm{Op}^{w}(a(x/t,\xi)). For uCc(d,)u\in C_{c}^{\infty}(\mathbb{R}^{d};\mathbb{C}), define Muf=|u|4/(d1)fM_{u}f=|u|^{4/(d-1)}f. Then

(6.1) i(AtMuMuAt)u,u=2d+1t1d(b)(x/t)|u|2(d+1)/(d1)dx,\displaystyle i\langle(A_{t}M_{u}-M_{u}A_{t})u,u\rangle=-\frac{2}{d+1}t^{-1}\int_{\mathbb{R}^{d}}(\nabla\cdot b)(x/t)|u|^{2(d+1)/(d-1)}\,dx,

and this functional has a unique continuous extension to H1H^{1}.

Proof.

For fH1(d)f\in H^{1}(\mathbb{R}^{d}), the definition of Weyl quantization gives

(6.2) Atf=ib(x/t)fi2t(b)(x/t)f+c(x/t)f.\displaystyle A_{t}f=-ib(x/t)\cdot\nabla f-\frac{i}{2t}(\nabla\cdot b)(x/t)f+c(x/t)f.

The function Muu=|u|4/(d1)uM_{u}u=|u|^{4/(d-1)}u belongs to H1(d)H^{1}(\mathbb{R}^{d}). Because bb, cc, and |u|4/(d1)|u|^{4/(d-1)} are real-valued, integration by parts shows that AtA_{t} is symmetric as a first-order form on H1H^{1} and that MuM_{u} is symmetric. Consequently,

(AtMuMuAt)u,u=Muu,AtuMuu,Atu¯.\displaystyle\langle(A_{t}M_{u}-M_{u}A_{t})u,u\rangle=\langle M_{u}u,A_{t}u\rangle-\overline{\langle M_{u}u,A_{t}u\rangle}.

Using (6.2) with f=uf=u, we obtain

i(AtMuMuAt)u,u\displaystyle i\langle(A_{t}M_{u}-M_{u}A_{t})u,u\rangle =2iMuu,Atu\displaystyle=2\Re\langle iM_{u}u,A_{t}u\rangle
=db(x/t)|u|4/(d1)(|u|2)dx\displaystyle=-\int_{\mathbb{R}^{d}}b(x/t)\cdot|u|^{4/(d-1)}\nabla(|u|^{2})\,dx
t1d(b)(x/t)|u|2(d+1)/(d1)dx.\displaystyle\quad-t^{-1}\int_{\mathbb{R}^{d}}(\nabla\cdot b)(x/t)|u|^{2(d+1)/(d-1)}\,dx.

Since |u|4/(d1)(|u|2)=d1d+1(|u|2(d+1)/(d1))|u|^{4/(d-1)}\nabla(|u|^{2})=\frac{d-1}{d+1}\nabla(|u|^{2(d+1)/(d-1)}), integration by parts proves (6.1).

Using ||z|2(d+1)/(d1)|w|2(d+1)/(d1)|Cd(|z|(d+3)/(d1)+|w|(d+3)/(d1))|zw|\big||z|^{2(d+1)/(d-1)}-|w|^{2(d+1)/(d-1)}\big|\leq C_{d}\big(|z|^{(d+3)/(d-1)}+|w|^{(d+3)/(d-1)}\big)|z-w| and Hölder’s inequality, we have

|d(b)(x/t)(|u|2(d+1)/(d1)|u~|2(d+1)/(d1))𝑑x|\displaystyle\Big|\int_{\mathbb{R}^{d}}(\nabla\cdot b)(x/t)\big(|u|^{2(d+1)/(d-1)}-|\widetilde{u}|^{2(d+1)/(d-1)}\big)\,dx\Big| CdbL(uL2(d+1)/(d1)(d+3)/(d1)CLOSE\displaystyle\leq C_{d}\|\nabla\cdot b\|_{L^{\infty}}\big(\|u\|_{L^{2(d+1)/(d-1)}}^{(d+3)/(d-1)}
OPEN+u~L2(d+1)/(d1)(d+3)/(d1))uu~L2(d+1)/(d1).\displaystyle\qquad+\|\widetilde{u}\|_{L^{2(d+1)/(d-1)}}^{(d+3)/(d-1)}\big)\|u-\widetilde{u}\|_{L^{2(d+1)/(d-1)}}.

Combining this estimate with the Sobolev embedding in (2.10) shows that the right-hand side of (6.1) is continuous in the H1H^{1} norm, and the conclusion follows. ∎

Now we compute the time derivative of the action DAD_{A} below. This will be used in Proposition 7.4 and eventually this will be used to approximate LqL_{q} in (7.36).

Proposition 6.2.

Let u+H1(d)u_{+}\in H^{1}(\mathbb{R}^{d}), and let v(t)v(t) and r(t)r(t) be as in (1.4). Let

bCc(d,d),cCc(d,),a(y,ξ)=b(y)ξ+c(y),\displaystyle b\in C_{c}^{\infty}(\mathbb{R}^{d};\mathbb{R}^{d}),\qquad c\in C_{c}^{\infty}(\mathbb{R}^{d};\mathbb{R}),\qquad a(y,\xi)=b(y)\cdot\xi+c(y),

and define V=(2ξy)yV=(2\xi-y)\cdot\nabla_{y}. Let At=Opw(a(x/t,ξ))A_{t}=\mathrm{Op}^{w}(a(x/t,\xi)), which is the same as in (6.2). Then

(6.3) DA(t)=Atu(t),u(t)Atv(t),v(t)\displaystyle D_{A}(t)=\langle A_{t}u(t),u(t)\rangle-\langle A_{t}v(t),v(t)\rangle

is locally absolutely continuous for t1t\geq 1 and satisfies

(6.4) DA(t)=t1Opw((Va)(x/t,ξ))r(t),r(t)+2t1Opw((Va)(x/t,ξ))r(t),v(t)2d+1t1d(b)(x/t)|u(t,x)|2(d+1)/(d1)dx\displaystyle\begin{aligned} D_{A}^{\prime}(t)&=t^{-1}\langle\mathrm{Op}^{w}((Va)(x/t,\xi))r(t),r(t)\rangle+2t^{-1}\Re\langle\mathrm{Op}^{w}((Va)(x/t,\xi))r(t),v(t)\rangle\\ &\quad-\frac{2}{d+1}t^{-1}\int_{\mathbb{R}^{d}}(\nabla\cdot b)(x/t)|u(t,x)|^{2(d+1)/(d-1)}\,dx\end{aligned}

for almost every t1t\geq 1. The first two pairings on the right-hand side of (6.4) are H1H^{-1}H1H^{1} pairings.

Proof.

Integration by parts yields

Atf,f=db(x/t)(f¯f)dx+dc(x/t)|f|2dx.\displaystyle\langle A_{t}f,f\rangle=\int_{\mathbb{R}^{d}}b(x/t)\cdot\Im(\overline{f}\nabla f)\,dx+\int_{\mathbb{R}^{d}}c(x/t)|f|^{2}\,dx.

Consequently,

(6.5) |Atf,f|\displaystyle|\langle A_{t}f,f\rangle| bLfL2fL2+cLfL22fH12.\displaystyle\leq\|b\|_{L^{\infty}}\|f\|_{L^{2}}\|\nabla f\|_{L^{2}}+\|c\|_{L^{\infty}}\|f\|_{L^{2}}^{2}\lesssim\|f\|_{H^{1}}^{2}.

We next prove (6.4). Using

tAt=t1Opw((yya)(x/t,ξ)),i((Δ)AtAt(Δ))=2t1Opw((ξya)(x/t,ξ)),\displaystyle\partial_{t}A_{t}=-t^{-1}\mathrm{Op}^{w}((y\cdot\nabla_{y}a)(x/t,\xi)),\quad i\big((-\Delta)A_{t}-A_{t}(-\Delta)\big)=2t^{-1}\mathrm{Op}^{w}((\xi\cdot\nabla_{y}a)(x/t,\xi)),

we have

(6.6) tAt+i((Δ)AtAt(Δ))=t1Bt,\displaystyle\partial_{t}A_{t}+i\big((-\Delta)A_{t}-A_{t}(-\Delta)\big)=t^{-1}B_{t},

where Bt=Opw((Va)(x/t,ξ))B_{t}=\mathrm{Op}^{w}((Va)(x/t,\xi)). Since

(6.7) Va=\displaystyle Va= ξT(Db+DbT)ξ+(2cDby)ξyc,\displaystyle\xi^{T}(Db+Db^{T})\xi+(2\nabla c-Db\,y)\cdot\xi-y\cdot\nabla c,

BtB_{t} is a second-order pseudodifferential operator with symbol bounds uniform in t1t\geq 1. Thus Bt:H1H1B_{t}:H^{1}\to H^{-1} is bounded uniformly in tt and Bt,\langle B_{t}\cdot,\cdot\rangle is symmetric. In particular,

(6.8) |Btf,g|fH1gH1,f,gH1(d),t1.\displaystyle|\langle B_{t}f,g\rangle|\lesssim\|f\|_{H^{1}}\|g\|_{H^{1}},\qquad f,g\in H^{1}(\mathbb{R}^{d}),\quad t\geq 1.

So we can group and ungroup pairings below.

A direct computation shows

(6.9) DA(t)\displaystyle D_{A}^{\prime}(t) =t1(Btu(t),u(t)Btv(t),v(t))2d+1t1d(b)(x/t)|u(t,x)|2(d+1)/(d1)𝑑x.\displaystyle=t^{-1}\big(\langle B_{t}u(t),u(t)\rangle-\langle B_{t}v(t),v(t)\rangle\big)-\frac{2}{d+1}t^{-1}\int_{\mathbb{R}^{d}}(\nabla\cdot b)(x/t)|u(t,x)|^{2(d+1)/(d-1)}\,dx.

Finally, since u=r+vu=r+v, the symmetric property of BtB_{t} gives

Btu,uBtv,v=Btr,r+2Btr,v.\displaystyle\langle B_{t}u,u\rangle-\langle B_{t}v,v\rangle=\langle B_{t}r,r\rangle+2\Re\langle B_{t}r,v\rangle.

Substituting this into (6.9) proves (6.4). ∎

Now we justify the following intuition. Suppose that νt\nu_{t}, which has density |r|2|r|^{2}, converges to a measure μ\mu with suppμK\mathrm{supp}\mu\subset K and that suppχ\mathrm{supp}\chi is separated from KK. Then χ(x/t)r\chi(x/t)r should converge to 00.

Proposition 6.3.

Let rr be defined by (1.4), and let νt\nu_{t} be the residual mass measure defined in (1.12). Let 1AnBn1\leq A_{n}\leq B_{n} and AnA_{n}\to\infty. Suppose there is a fixed closed set KdK\subset\mathbb{R}^{d} such that whenever nkn_{k}\to\infty, AnktkBnkA_{n_{k}}\leq t_{k}\leq B_{n_{k}}, and νtkμ\nu_{t_{k}}\rightharpoonup\mu weakly on X=d{}X=\mathbb{R}^{d}\cup\{\infty\}, one has suppμK\mathrm{supp}\mu\subset K. If χC1(d)\chi\in C^{1}(\mathbb{R}^{d}), χ,χL\chi,\nabla\chi\in L^{\infty}, and dist(suppχ,K)>0\mathrm{dist}(\mathrm{supp}\chi,K)>0, then

(6.10) supAntBnχ(x/t)r(t)H10.\displaystyle\sup_{A_{n}\leq t\leq B_{n}}\|\chi(x/t)r(t)\|_{H^{1}}\to 0.
Proof.

If (6.10) fails, then, after passing to a subsequence, there are ϵ>0\epsilon>0, indices nkn_{k}\to\infty, and AnktkBnkA_{n_{k}}\leq t_{k}\leq B_{n_{k}} such that

χ(x/tk)r(tk)H1ϵ.\displaystyle\|\chi(x/t_{k})r(t_{k})\|_{H^{1}}\geq\epsilon.

Since AnkA_{n_{k}}\to\infty, we have tkt_{k}\to\infty. By Proposition 5.1, after passing to a further subsequence, there are w1,,wJH1w_{1},\ldots,w_{J}\in H^{1}, centers xj,kx_{j,k}, points yjXy_{j}\in X, and ek0e_{k}\to 0 in H1H^{1} such that

(6.11) r(tk)=j=1Jwj(xj,k)+ek,xj,ktk\displaystyle r(t_{k})=\sum_{j=1}^{J}w_{j}(\cdot-x_{j,k})+e_{k},\quad\frac{x_{j,k}}{t_{k}} yj,νtkj=1JM(wj)δyj.\displaystyle\to y_{j},\quad\nu_{t_{k}}\rightharpoonup\sum_{j=1}^{J}M(w_{j})\delta_{y_{j}}.

The hypothesis on weak limits therefore implies wj0yjKw_{j}\neq 0\Longrightarrow y_{j}\in K. For every such jj and fixed zdz\in\mathbb{R}^{d},

xj,ktk+ztkyj,χ(xj,ktk+ztk)0,\displaystyle\frac{x_{j,k}}{t_{k}}+\frac{z}{t_{k}}\to y_{j},\qquad\chi\Big(\frac{x_{j,k}}{t_{k}}+\frac{z}{t_{k}}\Big)\to 0,

where the second limit follows from dist(suppχ,K)>0\mathrm{dist}(\mathrm{supp}\chi,K)>0. The product rule gives

(6.12) (χ(x/t)wj)=χ(x/t)wj+t1(χ)(x/t)wj,\displaystyle\nabla(\chi(x/t)w_{j})=\chi(x/t)\nabla w_{j}+t^{-1}(\nabla\chi)(x/t)w_{j},

The term containing χ\nabla\chi is bounded by tk1χLwjL2t_{k}^{-1}\|\nabla\chi\|_{L^{\infty}}\|w_{j}\|_{L^{2}}. Hence we have

(6.13) χ(x/tk)wj(xj,k)H10.\displaystyle\|\chi(x/t_{k})w_{j}(\cdot-x_{j,k})\|_{H^{1}}\to 0.

Similarly, (6.12) with wjw_{j} replaced by eke_{k} gives

(6.14) χ(x/tk)ekH1C(χL+χL)ekH10.\displaystyle\|\chi(x/t_{k})e_{k}\|_{H^{1}}\leq C(\|\chi\|_{L^{\infty}}+\|\nabla\chi\|_{L^{\infty}})\|e_{k}\|_{H^{1}}\to 0.

Applying (6.13) and (6.14) to (6.11) contradicts the lower bound ϵ\epsilon. ∎

Conversely, if (1.15) fails (or equivalently (3.14) holds), then Lemma 3.3 gives a time sequence along which the H1H^{1} norm of localized rr is bounded below. Now we upgrade this to an L2L^{2} lower bound.

Proposition 6.4.

Let d5d\geq 5, and let uC(0,H1(d))u\in C(\mathbb{R}_{\geq 0};H^{1}(\mathbb{R}^{d})) solve

(it+Δ)u=|u|4/(d1)u\displaystyle(i\partial_{t}+\Delta)u=-|u|^{4/(d-1)}u

and satisfy (1.2). Let u+H1(d)u_{+}\in H^{1}(\mathbb{R}^{d}) be the radiation state as in Theorem 2.7, and let rr be defined by (1.4). Let χ0,χ1Cc(d)\chi_{0},\chi_{1}\in C_{c}^{\infty}(\mathbb{R}^{d}) satisfy

χ1=1on suppχ0.\displaystyle\chi_{1}=1\qquad\text{on }\mathrm{supp}\chi_{0}.

Suppose that TnT_{n}\to\infty,

(6.15) χ1(x/Tn)r(Tn)H10,\displaystyle\|\chi_{1}(x/T_{n})r(T_{n})\|_{H^{1}}\to 0,

and that (3.14) holds. After passing to a subsequence and relabeling, let ϵ>0\epsilon>0 and sn>Tns_{n}>T_{n} be the times given by Lemma 3.3, so that

(6.16) χ0(x/sn)r(sn)H1\displaystyle\|\chi_{0}(x/s_{n})r(s_{n})\|_{H^{1}} =ϵ.\displaystyle=\epsilon.

Then there are m>0m_{*}>0 and n0n_{0} such that

(6.17) χ0(x/sn)r(sn)L2\displaystyle\|\chi_{0}(x/s_{n})r(s_{n})\|_{L^{2}} m,nn0.\displaystyle\geq m_{*},\qquad n\geq n_{0}.
Proof.

We first prove

(6.18) limLlim supnP>2L(χ0(x/sn)r(sn))H1=0.\displaystyle\lim_{L\to\infty}\limsup_{n\to\infty}\|P_{>2L}(\chi_{0}(x/s_{n})r(s_{n}))\|_{H^{1}}=0.

Here P>LP_{>L} and PLP_{\leq L} denote the Fourier projections to {|ξ|>L}\{|\xi|>L\} and {|ξ|L}\{|\xi|\leq L\}, respectively.

For L1L\geq 1, decompose

(6.19) P>2L(χ0(x/sn)r(sn))\displaystyle P_{>2L}(\chi_{0}(x/s_{n})r(s_{n})) =P>2L(χ0(x/sn)P>Lr(sn))+P>2L(χ0(x/sn)PLr(sn)).\displaystyle=P_{>2L}(\chi_{0}(x/s_{n})P_{>L}r(s_{n}))+P_{>2L}(\chi_{0}(x/s_{n})P_{\leq L}r(s_{n})).

Multiplication by χ0(x/sn)\chi_{0}(x/s_{n}) is bounded on H1H^{1} (for sn1s_{n}\geq 1). So we have

(6.20) P>2L(χ0(x/sn)P>Lr(sn))H1Cχ0P>Lr(sn)H1.\displaystyle\|P_{>2L}(\chi_{0}(x/s_{n})P_{>L}r(s_{n}))\|_{H^{1}}\leq C_{\chi_{0}}\|P_{>L}r(s_{n})\|_{H^{1}}.

For the second term on the right-hand side of (6.19), we take Fourier transform. Then we write the Fourier transform of P>2L(χ0(x/sn)PLr(sn))P_{>2L}(\chi_{0}(x/s_{n})P_{\leq L}r(s_{n})) as a convolution. If the total frequency satisfies |ξ|>2L|\xi|>2L and the frequency in PLr(sn)^\widehat{P_{\leq L}r(s_{n})} satisfies |η|L|\eta|\leq L, then |ξη|>L|\xi-\eta|>L. So we have

(6.21) P>2L(χ0(x/sn)PLr(sn))H1\displaystyle\|P_{>2L}(\chi_{0}(x/s_{n})P_{\leq L}r(s_{n}))\|_{H^{1}} C|PLr(sn)||ζ|>snLH1(1+|ζ|)|χ0^(ζ)|𝑑ζ.\displaystyle\leq C\|P_{\leq L}r(s_{n})\|_{H^{1}}\int_{|\zeta|>s_{n}L}(1+|\zeta|)|\widehat{\chi_{0}}(\zeta)|\,d\zeta.

Since χ0Cc(d)\chi_{0}\in C_{c}^{\infty}(\mathbb{R}^{d}), its Fourier transform is rapidly decreasing. Hence, for every fixed LL,

|ζ|>snL(1+|ζ|)|χ0^(ζ)|𝑑ζ0.\displaystyle\int_{|\zeta|>s_{n}L}(1+|\zeta|)|\widehat{\chi_{0}}(\zeta)|\,d\zeta\to 0.

Since we have supt0r(t)H1<\sup_{t\geq 0}\|r(t)\|_{H^{1}}<\infty, combining (6.20) and (6.21) gives

lim supnP>2L(χ0(x/sn)r(sn))H1Cχ0lim supnP>Lr(sn)H1,\displaystyle\limsup_{n\to\infty}\|P_{>2L}(\chi_{0}(x/s_{n})r(s_{n}))\|_{H^{1}}\leq C_{\chi_{0}}\limsup_{n\to\infty}\|P_{>L}r(s_{n})\|_{H^{1}},

which gives (6.18) by Proposition 2.8.

Choose LL sufficiently large and then n0n_{0} sufficiently large so that

P>2L(χ0(x/sn)r(sn))H1ϵ2,nn0.\displaystyle\|P_{>2L}(\chi_{0}(x/s_{n})r(s_{n}))\|_{H^{1}}\leq\frac{\epsilon}{2},\qquad n\geq n_{0}.

Since

χ0(x/sn)r(sn)=P2L(χ0(x/sn)r(sn))+P>2L(χ0(x/sn)r(sn))\displaystyle\chi_{0}(x/s_{n})r(s_{n})=P_{\leq 2L}(\chi_{0}(x/s_{n})r(s_{n}))+P_{>2L}(\chi_{0}(x/s_{n})r(s_{n}))

and χ0(x/sn)r(sn)H1=ϵ\|\chi_{0}(x/s_{n})r(s_{n})\|_{H^{1}}=\epsilon, the triangle inequality gives

P2L(χ0(x/sn)r(sn))H1ϵ2.\displaystyle\|P_{\leq 2L}(\chi_{0}(x/s_{n})r(s_{n}))\|_{H^{1}}\geq\frac{\epsilon}{2}.

Since 1+|ξ|21+4L21+|\xi|^{2}\leq 1+4L^{2} on the Fourier support of P2LP_{\leq 2L}, we have

P2L(χ0(x/sn)r(sn))H1\displaystyle\|P_{\leq 2L}(\chi_{0}(x/s_{n})r(s_{n}))\|_{H^{1}} 1+4L2P2L(χ0(x/sn)r(sn))L21+4L2χ0(x/sn)r(sn)L2.\displaystyle\leq\sqrt{1+4L^{2}}\|P_{\leq 2L}(\chi_{0}(x/s_{n})r(s_{n}))\|_{L^{2}}\leq\sqrt{1+4L^{2}}\|\chi_{0}(x/s_{n})r(s_{n})\|_{L^{2}}.

It follows that

χ0(x/sn)r(sn)L2ϵ21+4L2=:m>0\displaystyle\|\chi_{0}(x/s_{n})r(s_{n})\|_{L^{2}}\geq\frac{\epsilon}{2\sqrt{1+4L^{2}}}=:m_{*}>0

for every nn0n\geq n_{0}, which proves (6.17). ∎

7. The limiting measure and the proof of Theorem 1.2

We prove Theorem 1.2 in this section. We sketch the overall idea first. Let Ω\Omega_{\sharp} be the set of subsequential limits of (νt,Jt)(\nu_{t},J_{t}) as tt\to\infty in Proposition 5.4, which is compact and connected. Now the convergence of (νt,Jt)(\nu_{t},J_{t}) is equivalent to this set being a single point. We first find a convex function Ψ\Psi such that Ψ𝑑μ\int\Psi d\mu (as a function of μ\mu) has a unique maximizer μ\mu_{*}. Then we prove that elements of Ω\Omega_{\sharp} have the following rigidity: if the first component is μ\mu_{*}, then the momentum component has to be J=y2μJ_{*}=\frac{y}{2}\mu_{*}. If the convergence fails, then there will be two time sequences such that (νt,Jt)(\nu_{t},J_{t}) is very close to (μ,J)(\mu_{*},J_{*}) along one sequence but stays a fixed positive distance from (μ,J)(\mu_{*},J_{*}) along the other.

Then we will show that the elements of Ω\Omega_{\sharp} have the following coercive property: if (μ,J)(\mu,J) stays a fixed positive distance from (μ,J)(\mu_{*},J_{*}), then Ψ𝑑μ\int\Psi\,d\mu is below Ψdμ\int\Psi\,d\mu_{*} by a fixed positive amount. Then we show that this is ‘repulsive’ enough to make the deficit from Ψdμ\int\Psi\,d\mu_{*} grow exponentially in the delayed time parameter |τ||\tau| from Theorem 5.6. This will lead to a contradiction. The last step is difficult, so we will spend most of this section computing and estimating the currents involved in this step.

We will consider (νt,Jt)(\nu_{t},J_{t}) reparametrized on an exponential scale:

(7.1) F(s)=(νes,Jes).\displaystyle F(s)=(\nu_{e^{s}},J_{e^{s}}).
Proposition 7.1.

Assume FF in (7.1) does not converge as ss\to\infty. There is a real polynomial Ψ\Psi, a neighborhood UU of a common ball containing the supports of all subsequential limits, and λΨ>0\lambda_{\Psi}>0 such that

zTD2Ψ(y)zλΨ|z|2\displaystyle z^{T}D^{2}\Psi(y)z\geq\lambda_{\Psi}|z|^{2}

for yUy\in U, and

Ψ𝑑μ\displaystyle\int\Psi\,d\mu

has a unique maximizer μ\mu_{*} over the projection of the set Ω\Omega_{\sharp} in Proposition 5.4 to the mass component.

Proof.

By Proposition 5.8 and (5.14), for an element of Ω\Omega_{\sharp}, the mass component (i.e., the μ\mu-component) is supported in one fixed compact ball and has at most

Nmax=R/Θ\displaystyle N_{\max}=\big\lfloor\sqrt{R_{\sharp}/\Theta}\big\rfloor

atoms. The coordinate monomials of total degree at most 2Nmax12N_{\max}-1 separate all such measures. Indeed, the difference of two such measures has at most 2Nmax2N_{\max} support points, and Lagrange-type polynomials of degree at most 2Nmax12N_{\max}-1 recover every signed weight.

Let g1,,gAg_{1},\ldots,g_{A} be those monomials and ϕ0(y)=|y|2\phi_{0}(y)=|y|^{2}. Choose CjC_{j} large enough that

ϕj=Cjϕ0+gj\displaystyle\phi_{j}=C_{j}\phi_{0}+g_{j}

has positive-definite Hessian on a neighborhood of the support ball. The moment map

L(μ)=(ϕ0𝑑μ,,ϕA𝑑μ)\displaystyle L(\mu)=\Big(\int\phi_{0}\,d\mu,\ldots,\int\phi_{A}\,d\mu\Big)

is continuous and injective on the compact set Ωmass\Omega_{\mathrm{mass}} formed by the projection of Ω\Omega_{\sharp} to the mass component.

For a compact set CA+1C\subset\mathbb{R}^{A+1}, define

hC(w)=maxxCwx.\displaystyle h_{C}(w)=\max_{x\in C}w\cdot x.

This function is Lipschitz in ww, hence differentiable almost everywhere. If xCx\in C maximizes wxw\cdot x, then

hC(w+z)hC(w)+zxzA+1.\displaystyle h_{C}(w+z)\geq h_{C}(w)+z\cdot x\;\forall z\in\mathbb{R}^{A+1}.

At a point where hCh_{C} is differentiable, this inequality, applied to both zz and z-z, forces x=hC(w)x=\nabla h_{C}(w). Thus the maximizer xx is unique. Apply this with C=L(Ωmass)C=L(\Omega_{\mathrm{mass}}), choose such a vector w=(w0,,wA)w=(w_{0},\ldots,w_{A}) with every wj>0w_{j}>0, and set

Ψ=j=0Awjϕj.\displaystyle\Psi=\sum_{j=0}^{A}w_{j}\phi_{j}.

Because every wjw_{j} is positive, the Hessian of Ψ\Psi has a uniform positive lower bound on a neighborhood of the support ball. The injectivity of LL then gives a unique measure maximizing Ψ𝑑μ\int\Psi\,d\mu over Ωmass\Omega_{\mathrm{mass}}. ∎

Now we prove the rigidity of the limiting measures in Ω\Omega_{\sharp}: if the first component is μ\mu_{*}, then the momentum component is determined by it as well. That is, yy is indeed the velocity, and the momentum equals one-half of the velocity times the mass.

Proposition 7.2 (Uniqueness at the maximizing measure).

Let Ψ\Psi and

μ=i=1Laiδxi,ai>0,\displaystyle\mu_{*}=\sum_{i=1}^{L}a_{i}^{*}\delta_{x_{i}^{*}},\qquad a_{i}^{*}>0,

be as in Proposition 7.1, with the xix_{i}^{*} distinct, and set

(7.2) Z=(μ,y2μ).\displaystyle Z_{*}=\Big(\mu_{*},\frac{y}{2}\mu_{*}\Big).

If (μ(τ),J(τ))(\mu(\tau),J(\tau)) is a subsequential limit family from Theorem 5.6 and μ(0)=μ\mu(0)=\mu_{*}, then

(7.3) (μ(τ),J(τ))=Z(τ).\displaystyle(\mu(\tau),J(\tau))=Z_{*}\qquad(\tau\in\mathbb{R}).

Consequently, ZZ_{*} is the only subsequential limit of FF in (7.1) whose first component is μ\mu_{*}.

Proof.

Let (μ(τ),J(τ))(\mu(\tau),J(\tau)) be such a subsequential limit family, obtained along sns_{n}\to\infty, and put

H(τ)=Ψ𝑑μ(τ).\displaystyle H(\tau)=\int\Psi\,d\mu(\tau).

For each fixed τ\tau, (5.28) shows that (μ(τ),J(τ))(\mu(\tau),J(\tau)) is a subsequential limit of FF. The maximizing property in Proposition 7.1 therefore gives

H(τ)H(0)(τ).\displaystyle H(\tau)\leq H(0)\qquad(\tau\in\mathbb{R}).

Let UU and λΨ\lambda_{\Psi} be as in Proposition 7.1. Choose pairwise disjoint balls UiU_{i} containing xix_{i}^{*}, with Ui¯U\overline{U_{i}}\subset U. The dJd_{J}-continuity in Theorem 5.6 and the atom-mass lower bound in Proposition 5.3 imply that, for some ϵ>0\epsilon>0,

(7.4) suppμ(τ)i=1LUi(|τ|<ϵ).\displaystyle\mathrm{supp}\mu(\tau)\subset\bigcup_{i=1}^{L}U_{i}\qquad(|\tau|<\epsilon).

Indeed, otherwise an atom of uniformly positive mass would remain outside these balls along a sequence τ0\tau\to 0, contradicting μ(τ)μ\mu(\tau)\rightharpoonup\mu_{*}.

Choose χiCc(d)\chi_{i}\in C_{c}^{\infty}(\mathbb{R}^{d}) equal to one on a neighborhood of Ui¯\overline{U_{i}} and zero near every Uj¯\overline{U_{j}} with jij\neq i, and define

Pi(τ)=χi𝑑J(τ),Bi(τ)=yχi(y)𝑑μ(τ)(y).\displaystyle P_{i}(\tau)=\int\chi_{i}\,dJ(\tau),\qquad B_{i}(\tau)=\int y\chi_{i}(y)\,d\mu(\tau)(y).

The support containment (7.4) and the support relation in Theorem 5.6 show that χi\nabla\chi_{i} vanishes on the supports of both μ(τ)\mu(\tau) and J(τ)J(\tau). Thus (5.29), first with ϕ=χi\phi=\chi_{i}, gives

χi𝑑μ(τ)=ai(|τ|<ϵ).\displaystyle\int\chi_{i}\,d\mu(\tau)=a_{i}^{*}\qquad(|\tau|<\epsilon).

Now we pass to the further subsequence given by Proposition 5.7. This does not change the subsequential limit family (μ(τ),J(τ))(\mu(\tau),J(\tau)) because of (5.28). By (5.45), we have

(7.5) suppΣkm(τ)i=1LUi¯\displaystyle\mathrm{supp}\Sigma_{km}(\tau)\subset\bigcup_{i=1}^{L}\overline{U_{i}}

for almost every |τ|<ϵ|\tau|<\epsilon. Since χi\nabla\chi_{i} vanishes on the union in (7.5), (5.43) shows that Pi(τ)P_{i}(\tau) is constant and we will denote it by PiP_{i} below. Applying (5.29) with ϕ(y)=χi(y)yk\phi(y)=\chi_{i}(y)y_{k} for each coordinate kk gives

Bi=2PiBi.\displaystyle B_{i}^{\prime}=2P_{i}-B_{i}.

Since Bi(0)=aixiB_{i}(0)=a_{i}^{*}x_{i}^{*}, it follows that

(7.6) Bi(τ)aixi=(1eτ)(2Piaixi).\displaystyle B_{i}(\tau)-a_{i}^{*}x_{i}^{*}=(1-e^{-\tau})(2P_{i}-a_{i}^{*}x_{i}^{*}).

The Hessian bound in Proposition 7.1 gives, for yUiy\in U_{i},

Ψ(y)Ψ(xi)+Ψ(xi)(yxi)+λΨ2|yxi|2.\displaystyle\Psi(y)\geq\Psi(x_{i}^{*})+\nabla\Psi(x_{i}^{*})\cdot(y-x_{i}^{*})+\frac{\lambda_{\Psi}}{2}|y-x_{i}^{*}|^{2}.

Integrating over the UiU_{i}, summing, using (7.6), and applying Cauchy–Schwarz yields

(7.7) H(τ)H(0)\displaystyle H(\tau)-H(0) (1eτ)iΨ(xi)(2Piaixi)+λΨ2iUi|yxi|2𝑑μ(τ)(y)\displaystyle\geq(1-e^{-\tau})\sum_{i}\nabla\Psi(x_{i}^{*})\cdot(2P_{i}-a_{i}^{*}x_{i}^{*})+\frac{\lambda_{\Psi}}{2}\sum_{i}\int_{U_{i}}|y-x_{i}^{*}|^{2}\,d\mu(\tau)(y)
(1eτ)iΨ(xi)(2Piaixi)+λΨ2(1eτ)2i|2Piaixi|2ai(|τ|<ϵ).\displaystyle\geq(1-e^{-\tau})\sum_{i}\nabla\Psi(x_{i}^{*})\cdot(2P_{i}-a_{i}^{*}x_{i}^{*})+\frac{\lambda_{\Psi}}{2}(1-e^{-\tau})^{2}\sum_{i}\frac{|2P_{i}-a_{i}^{*}x_{i}^{*}|^{2}}{a_{i}^{*}}\qquad(|\tau|<\epsilon).

Together with H(τ)H(0)H(\tau)\leq H(0) for both signs of τ\tau, this first makes the coefficient of 1eτ1-e^{-\tau} vanish and then gives

2Pi=aixi(1iL).\displaystyle 2P_{i}=a_{i}^{*}x_{i}^{*}\qquad(1\leq i\leq L).

Substituting 2Pi=aixi2P_{i}=a_{i}^{*}x_{i}^{*} into (7.7) and using H(τ)H(0)H(\tau)\leq H(0) give

0H(τ)H(0)λΨ2iUi|yxi|2𝑑μ(τ)(y).\displaystyle 0\geq H(\tau)-H(0)\geq\frac{\lambda_{\Psi}}{2}\sum_{i}\int_{U_{i}}|y-x_{i}^{*}|^{2}\,d\mu(\tau)(y).

Hence μ(τ)=μ\mu(\tau)=\mu_{*} for |τ|<ϵ|\tau|<\epsilon. The support relation in Theorem 5.6, together with Pi=aixi/2P_{i}=a_{i}^{*}x_{i}^{*}/2, shows that PiP_{i} is the vector weight of J(τ)J(\tau) at xix_{i}^{*} and gives

J(τ)=y2μ(|τ|<ϵ).\displaystyle J(\tau)=\frac{y}{2}\mu_{*}\qquad(|\tau|<\epsilon).

The set of τ\tau for which (μ(τ),J(τ))=Z(\mu(\tau),J(\tau))=Z_{*} is closed by the dJd_{J}-continuity in Theorem 5.6. On the other hand, if (μ(τ0),J(τ0))=Z(\mu(\tau_{0}),J(\tau_{0}))=Z_{*}, then (5.28) shows that the family (μ(τ0+σ),J(τ0+σ))(\mu(\tau_{0}+\sigma),J(\tau_{0}+\sigma)) is the subsequential limit family obtained along sn+τ0s_{n}+\tau_{0}, so the preceding local argument applies at τ0\tau_{0}. So this set of τ\tau is also open. Hence (μ(τ),J(τ))Z(\mu(\tau),J(\tau))\equiv Z_{*}.

Finally, let (μ,J)(\mu_{*},J) be any subsequential limit of FF. By the definition in Proposition 5.4, there is a sequence sns_{n}\to\infty such that F(sn)(μ,J)F(s_{n})\to(\mu_{*},J). After passing to a subsequence, Theorem 5.6 and (5.28) give a subsequential limit family (μ(τ),J(τ))(\mu(\tau),J(\tau)) with (μ(0),J(0))=(μ,J)(\mu(0),J(0))=(\mu_{*},J). Equation (7.3), now proved, gives (μ,J)=Z(\mu_{*},J)=Z_{*}. At least one such point exists by the definition of μ\mu_{*} in Proposition 7.1. ∎

Now we approximate Ψ\Psi by a simpler function that is almost piecewise linear. The advantage of using this qq rather than using Ψ\Psi above directly is that its second- and higher-order derivatives are easier to control.

Proposition 7.3.

Let Ψ\Psi and μ\mu_{*} be as in Proposition 7.1, with μ=jajδxj\mu_{*}=\sum_{j}a_{j}^{*}\delta_{x_{j}^{*}}. There are disjoint balls UjU_{j} about xjx_{j}^{*}, a smooth convex function q:dq:\mathbb{R}^{d}\to\mathbb{R}, and nested smooth cutoffs χ0,θ~,χ1\chi_{0},\tilde{\theta},\chi_{1} such that:

  1. (1)

    For every jj and every yUjy\in U_{j},

    (7.8) q(y)=Ψ(xj)+Ψ(xj)(yxj).\displaystyle q(y)=\Psi(x_{j}^{*})+\nabla\Psi(x_{j}^{*})\cdot(y-x_{j}^{*}).
  2. (2)

    D2q0D^{2}q\geq 0, and q\nabla q and yq(y)q(y)y\cdot\nabla q(y)-q(y), together with all their derivatives, are bounded. Moreover, D2q=0D^{2}q=0 outside {χ0=1}\{\chi_{0}=1\}, θ~=1\tilde{\theta}=1 near the supports of χ0,χ0,Δχ0\chi_{0},\nabla\chi_{0},\Delta\chi_{0}, χ1=1\chi_{1}=1 near suppθ~\mathrm{supp}\tilde{\theta}, and yχ0(y)y\cdot\nabla\chi_{0}(y) is bounded.

  3. (3)

    dist(suppχ1,jUj¯)>0\mathrm{dist}(\mathrm{supp}\chi_{1},\bigcup_{j}\overline{U_{j}})>0.

Proof.

Let U,λΨ>0U,\lambda_{\Psi}>0 be as in Proposition 7.1. In addition, we can choose it so that the line segment joining any two atoms xj,xkx_{j}^{*},x_{k}^{*} lies in UU. For each atom xjx_{j}^{*}, define the affine function

j(y)\displaystyle\ell_{j}(y) =Ψ(xj)+Ψ(xj)(yxj).\displaystyle=\Psi(x_{j}^{*})+\nabla\Psi(x_{j}^{*})\cdot(y-x_{j}^{*}).

If jkj\neq k, as Ψ\Psi is strictly convex (see Proposition 7.1), we have

j(xj)k(xj)\displaystyle\ell_{j}(x_{j}^{*})-\ell_{k}(x_{j}^{*}) =Ψ(xj)Ψ(xk)Ψ(xk)(xjxk)\displaystyle=\Psi(x_{j}^{*})-\Psi(x_{k}^{*})-\nabla\Psi(x_{k}^{*})\cdot(x_{j}^{*}-x_{k}^{*})
=01(1s)(xjxk)TD2Ψ(xk+s(xjxk))(xjxk)𝑑s\displaystyle=\int_{0}^{1}(1-s)(x_{j}^{*}-x_{k}^{*})^{T}D^{2}\Psi(x_{k}^{*}+s(x_{j}^{*}-x_{k}^{*}))(x_{j}^{*}-x_{k}^{*})\,ds
λΨ2|xjxk|2>0.\displaystyle\geq\frac{\lambda_{\Psi}}{2}|x_{j}^{*}-x_{k}^{*}|^{2}>0.

Since there are only finitely many atoms, we may choose rj>0r_{j}>0 such that the balls B(xj,2rj)B(x_{j}^{*},2r_{j}) are pairwise disjoint, are contained in UU, and satisfy

j(y)\displaystyle\ell_{j}(y) >k(y)for yB(xj,2rj),kj.\displaystyle>\ell_{k}(y)\qquad\text{for }y\in B(x_{j}^{*},2r_{j}),\ k\neq j.

Define

m(y)\displaystyle m(y) =maxjj(y).\displaystyle=\max_{j}\ell_{j}(y).

Then mm is convex and globally Lipschitz, and m=jm=\ell_{j} on B(xj,2rj)B(x_{j}^{*},2r_{j}). Choose a nonnegative radially symmetric function ρϵCc(d)\rho_{\epsilon}\in C_{c}^{\infty}(\mathbb{R}^{d}) such that

suppρϵ\displaystyle\mathrm{supp}\rho_{\epsilon} B(0,ϵ),dρϵ(z)𝑑z=1,0<ϵ<minjrj,\displaystyle\subset B(0,\epsilon),\qquad\int_{\mathbb{R}^{d}}\rho_{\epsilon}(z)\,dz=1,\qquad 0<\epsilon<\min_{j}r_{j},

and define

(7.9) q=mρϵ.\displaystyle q=m*\rho_{\epsilon}.

Because mm is convex and ρϵ0\rho_{\epsilon}\geq 0, qq is smooth and D2q0D^{2}q\geq 0.

Suppose that yB(xj,rj)y\in B(x_{j}^{*},r_{j}) and zsuppρϵz\in\mathrm{supp}\rho_{\epsilon}. Then yzB(xj,2rj)y-z\in B(x_{j}^{*},2r_{j}), so m(yz)=j(yz)m(y-z)=\ell_{j}(y-z). Since ρϵ\rho_{\epsilon} is even,

q(y)\displaystyle q(y) =dj(yz)ρϵ(z)𝑑z=j(y)Ψ(xj)dzρϵ(z)𝑑z=j(y).\displaystyle=\int_{\mathbb{R}^{d}}\ell_{j}(y-z)\rho_{\epsilon}(z)\,dz=\ell_{j}(y)-\nabla\Psi(x_{j}^{*})\cdot\int_{\mathbb{R}^{d}}z\rho_{\epsilon}(z)\,dz=\ell_{j}(y).

Thus

q(y)\displaystyle q(y) =Ψ(xj)+Ψ(xj)(yxj)for yB(xj,rj).\displaystyle=\Psi(x_{j}^{*})+\nabla\Psi(x_{j}^{*})\cdot(y-x_{j}^{*})\qquad\text{for }y\in B(x_{j}^{*},r_{j}).

Since mm is Lipschitz, its weak gradient is bounded. Therefore

q=(m)ρϵ,αq=(m)αρϵ.\displaystyle\nabla q=(\nabla m)*\rho_{\epsilon},\quad\partial^{\alpha}\nabla q=(\nabla m)*\partial^{\alpha}\rho_{\epsilon}.

Hence q\nabla q and all its derivatives are bounded.

At almost every point ww, one of the affine functions j\ell_{j} realizes the maximum defining mm and

wm(w)m(w)\displaystyle w\cdot\nabla m(w)-m(w) =xjΨ(xj)Ψ(xj).\displaystyle=x_{j}^{*}\cdot\nabla\Psi(x_{j}^{*})-\Psi(x_{j}^{*}).

The right side takes only finitely many values. Hence wm(w)m(w)w\cdot\nabla m(w)-m(w) is bounded almost everywhere. By the definition of qq in (7.9), we have

yq(y)q(y)\displaystyle y\cdot\nabla q(y)-q(y) =d((yz)m(yz)m(yz))ρϵ(z)𝑑z+dzm(yz)ρϵ(z)𝑑z.\displaystyle=\int_{\mathbb{R}^{d}}((y-z)\cdot\nabla m(y-z)-m(y-z))\rho_{\epsilon}(z)\,dz+\int_{\mathbb{R}^{d}}z\cdot\nabla m(y-z)\rho_{\epsilon}(z)\,dz.

The first term is the convolution of a bounded function with ρϵ\rho_{\epsilon}. The second is a finite sum of convolutions of the bounded functions νm\partial_{\nu}m with the compactly supported smooth functions zνρϵ(z)z_{\nu}\rho_{\epsilon}(z). Differentiation in yy only falls on the smooth one, and we know that yq(y)q(y)y\cdot\nabla q(y)-q(y) and all its derivatives are bounded.

We now construct the cutoffs. Set

Uj\displaystyle U_{j} =B(xj,rj/20).\displaystyle=B(x_{j}^{*},r_{j}/20).

Choose smooth cutoffs χ1,θ~,χ0:d\chi_{1},\tilde{\theta},\chi_{0}:\mathbb{R}^{d}\to\mathbb{R} satisfying 0χ1,θ~,χ010\leq\chi_{1},\tilde{\theta},\chi_{0}\leq 1 and the following prescribed values:

  1. (1)

    χ1(y)=0\chi_{1}(y)=0 if |yxj|rj/10|y-x_{j}^{*}|\leq r_{j}/10 for at least one jj, and χ1(y)=1\chi_{1}(y)=1 if |yxj|rj/5|y-x_{j}^{*}|\geq r_{j}/5 for every jj.

  2. (2)

    θ~(y)=0\tilde{\theta}(y)=0 if |yxj|rj/4|y-x_{j}^{*}|\leq r_{j}/4 for at least one jj, and θ~(y)=1\tilde{\theta}(y)=1 if |yxj|rj/2|y-x_{j}^{*}|\geq r_{j}/2 for every jj.

  3. (3)

    χ0(y)=0\chi_{0}(y)=0 if |yxj|3rj/5|y-x_{j}^{*}|\leq 3r_{j}/5 for at least one jj, and χ0(y)=1\chi_{0}(y)=1 if |yxj|4rj/5|y-x_{j}^{*}|\geq 4r_{j}/5 for every jj.

Then one can verify that all desired properties are satisfied. ∎

In the proof of the convergence of (νt,Jt)(\nu_{t},J_{t}) in (1.12), we will again use the approximation strategy. That is, we use j(u)j(v)j(u)-j(v) to approximate j(r)j(r) and use |u|2|v|2|u|^{2}-|v|^{2} to approximate |r|2|r|^{2}. So we consider

(7.10) Aq(t)\displaystyle A_{q}(t) =q(x/t)(j(u)j(v))𝑑x12(xtq(x/t)q(x/t))(|u|2|v|2)𝑑x,\displaystyle=\int\nabla q(x/t)\cdot\big(j(u)-j(v)\big)\,dx-\frac{1}{2}\int\Big(\frac{x}{t}\cdot\nabla q(x/t)-q(x/t)\Big)(|u|^{2}-|v|^{2})\,dx,

which is DAD_{A} in Proposition 6.2 with b=q(y)b=\nabla q(y) , c=12(yq(y)q(y))c=-\frac{1}{2}(y\cdot\nabla q(y)-q(y)). Here qq is the smooth convex function as in Proposition 7.3 above. Then we show that AqA_{q} is ‘non-decreasing’ in the following sense.

Proposition 7.4.

Let u+H1(d)u_{+}\in H^{1}(\mathbb{R}^{d}), and let v(t)v(t) and r(t)r(t) be as in (1.4). Let q,Aqq,A_{q} be as above. Let χ0,θ~,χ1\chi_{0},\tilde{\theta},\chi_{1} be as in Proposition 7.3, with D2q=0D^{2}q=0 outside {χ0=1}\{\chi_{0}=1\}. If AnA_{n}\to\infty, An<BnA_{n}<B_{n}, and

supAntBnχ1(x/t)r(t)H10,\displaystyle\sup_{A_{n}\leq t\leq B_{n}}\|\chi_{1}(x/t)r(t)\|_{H^{1}}\to 0,

then

lim infn(Aq(Bn)Aq(An))0.\displaystyle\liminf_{n\to\infty}\big(A_{q}(B_{n})-A_{q}(A_{n})\big)\geq 0.
Proof.

In the notation of Proposition 6.2, we have

a=q(y)ξyq(y)q(y)2.\displaystyle a=\nabla q(y)\cdot\xi-\frac{y\cdot\nabla q(y)-q(y)}{2}.

Our goal is to prove that the right-hand side of (6.4) is asymptotically nonnegative, up to a decaying negative term, as tt\to\infty. We still use V=(2ξy)yV=(2\xi-y)\cdot\nabla_{y} and we have

Va=V(q(y)ξyq(y)q(y)2)=12(2ξy)TD2q(y)(2ξy)0,\displaystyle Va=V\big(\nabla q(y)\cdot\xi-\frac{y\cdot\nabla q(y)-q(y)}{2}\big)=\frac{1}{2}(2\xi-y)^{T}D^{2}q(y)(2\xi-y)\geq 0,

where we used the convexity D2q0D^{2}q\geq 0.

Now we put things into the context of Proposition 6.2, where coefficients are required to have compact support 55 5 Another way is to relax the condition of Proposition 6.2 to only requiring bounds on b,cb,c for the proof to go through., by truncating objects above. Choose a bump function χCc(d)\chi\in C_{c}^{\infty}(\mathbb{R}^{d}), nonincreasing in |y||y|, such that χ=1\chi=1 on {|y|1}\{|y|\leq 1\} and χ=0\chi=0 on {|y|2}\{|y|\geq 2\}, and set

χR(y)=χ(y/R).\displaystyle\chi_{R}(y)=\chi(y/R).

Applying Proposition 6.2 with b=χR(y)q(y)b=\chi_{R}(y)\nabla q(y) and c=12χR(y)(yq(y)q(y))c=-\frac{1}{2}\chi_{R}(y)(y\cdot\nabla q(y)-q(y)) gives the corresponding truncated AqA_{q}:

(7.11) Aq,R(t)\displaystyle A_{q,R}(t) =χR(x/t)q(x/t)(j(u)j(v))dx12χR(x/t)(xtq(x/t)q(x/t))(|u|2|v|2)dx.\displaystyle=\int\chi_{R}(x/t)\nabla q(x/t)\cdot\big(j(u)-j(v)\big)\,dx-\frac{1}{2}\int\chi_{R}(x/t)(\frac{x}{t}\cdot\nabla q(x/t)-q(x/t))(|u|^{2}-|v|^{2})\,dx.

Then the corresponding truncated ‘VaVa’ in (6.7) is

(7.12) V(χR(y)(q(y)ξyq(y)q(y)2))=12χR(y)(2ξy)TD2q(y)(2ξy)+((2ξy)χR(y))(q(y)ξyq(y)q(y)2).\displaystyle\begin{aligned} &V(\chi_{R}(y)(\nabla q(y)\cdot\xi-\frac{y\cdot\nabla q(y)-q(y)}{2}))\\ &=\frac{1}{2}\chi_{R}(y)(2\xi-y)^{T}D^{2}q(y)(2\xi-y)+((2\xi-y)\cdot\nabla\chi_{R}(y))(\nabla q(y)\cdot\xi-\frac{y\cdot\nabla q(y)-q(y)}{2}).\end{aligned}

Now we consider the error term introduced by χR\chi_{R}. The second term on the right side of (7.12) is a polynomial of degree at most two in ξ\xi whose coefficients are O(R1)O(R^{-1}) and supported in {R|y|2R}\{R\leq|y|\leq 2R\}. As in the derivation of (6.8), its quantization is a pseudodifferential operator with O(R1)O(R^{-1}) symbols. The corresponding quadratic form is bounded by the H1H^{1} norms of uu and vv over {R|x/t|2R}\{R\leq|x/t|\leq 2R\}, with times O(R1)O(R^{-1}) level constant. Similarly, using (7.10) and (7.11), we know Aq,R(t)Aq(t)A_{q,R}(t)\to A_{q}(t) as RR\to\infty for every fixed tt.

The first term on the right side of (7.12) is nonnegative on the symbol level, and now we consider its operator level positivity after quantization. The symbol 2ξy2\xi-y quantizes to be Lt=2Dx/tL_{t}=2D-x/t, where D=iD=-i\nabla. We have

(7.13) Opw(12χR(x/t)(2ξx/t)TD2q(x/t)(2ξx/t))f,f=12χR(x/t)(Ltf)D2q(x/t)(Ltf)dx12t2j,k(yjyk(χR(D2q)jk))(x/t)|f|2dx.\displaystyle\begin{aligned} &\langle\mathrm{Op}^{w}\big(\frac{1}{2}\chi_{R}(x/t)(2\xi-x/t)^{T}D^{2}q(x/t)(2\xi-x/t)\big)f,f\rangle\\ &\quad=\frac{1}{2}\int\chi_{R}(x/t)(L_{t}f)^{*}D^{2}q(x/t)(L_{t}f)\,dx-\frac{1}{2t^{2}}\int\sum_{j,k}\big(\partial_{y_{j}}\partial_{y_{k}}(\chi_{R}(D^{2}q)_{jk})\big)(x/t)|f|^{2}\,dx.\end{aligned}

The first term is positive, and the second term is O(t2)O(t^{-2}), hence integrable.

The first integral on the right-hand side of (7.13) is non-decreasing in RR, and we define (allowed to be infinity):

(7.14) K(f,t)=limRχR(x/t)(Ltf)D2q(x/t)(Ltf)𝑑x.\displaystyle K(f,t)=\lim_{R\to\infty}\int\chi_{R}(x/t)(L_{t}f)^{*}D^{2}q(x/t)(L_{t}f)\,dx.

The strategy for the remaining part is as follows: we first take f=rf=r, which gives the term we currently have. If we add another multiple of K(v,t)K(v,t), then these two terms control the second term in (6.4). So the negative part we are left to control is K(v,t)K(v,t) and we show that it is integrable in time below. In particular, this shows that it has tail tending to zero.

More concretely, let AtA_{t} be as in Proposition 6.2 with χR\chi_{R} inserted, and consider

(7.15) Atv(t),v(t)=χR(x/t)q(x/t)j(v(t))dx12χR(x/t)(xtq(x/t)q(x/t))|v(t)|2dx.\displaystyle\langle A_{t}v(t),v(t)\rangle=\int\chi_{R}(x/t)\nabla q(x/t)\cdot j(v(t))\,dx-\frac{1}{2}\int\chi_{R}(x/t)(\frac{x}{t}\cdot\nabla q(x/t)-q(x/t))|v(t)|^{2}\,dx.

Let Bt=Opw((Va)(x/t,ξ))B_{t}=\mathrm{Op}^{w}((Va)(x/t,\xi)) with aa as above. Since (it+Δ)v=0(i\partial_{t}+\Delta)v=0, using (6.6), we have

(7.16) ddtAtv(t),v(t)=t1Btv(t),v(t),\displaystyle\frac{d}{dt}\langle A_{t}v(t),v(t)\rangle=t^{-1}\langle B_{t}v(t),v(t)\rangle,

which is the second term in the bracket on the right-hand side of (6.9). The right-hand side of (7.16) is evaluated in (7.13).

Now we integrate this from 11 to T1T\gg 1 and let RR\to\infty. Terms containing derivatives of χR\chi_{R} tend to zero as the discussion after (7.12) when RR\to\infty. Then we have

(q(x/t)j(v(t))𝑑x12(xtq(x/t)q(x/t))|v(t)|2𝑑x)|t=1t=T\displaystyle\big(\int\nabla q(x/t)\cdot j(v(t))\,dx-\frac{1}{2}\int(\frac{x}{t}\cdot\nabla q(x/t)-q(x/t))|v(t)|^{2}\,dx\big)\big|_{t=1}^{t=T}
=121TK(v,t)dtt121T(Δ2q)(x/t)|v(t,x)|2𝑑xdtt3.\displaystyle\quad=\frac{1}{2}\int_{1}^{T}K(v,t)\frac{dt}{t}-\frac{1}{2}\int_{1}^{T}\int(\Delta^{2}q)(x/t)|v(t,x)|^{2}\,dx\frac{dt}{t^{3}}.

The left-hand side is bounded uniformly in TT, because q\nabla q and yqqy\cdot\nabla q-q are bounded by the construction in Proposition 7.3 and v(t)H1\|v(t)\|_{H^{1}} is preserved. The last integral is bounded by Δ2qLu+L221t3𝑑t\|\Delta^{2}q\|_{L^{\infty}}\|u_{+}\|_{L^{2}}^{2}\int_{1}^{\infty}t^{-3}\,dt. Therefore

(7.17) 1K(v,t)dtt<.\displaystyle\int_{1}^{\infty}K(v,t)\frac{dt}{t}<\infty.

We now return to apply (6.4) to Aq,RA_{q,R} (playing the role of DAD_{A} there) from (7.11) to obtain

(7.18) Aq,R(t)=t1Opw((Va)(x/t,ξ))r(t),r(t)+2t1Opw((Va)(x/t,ξ))r(t),v(t)2d+1t1(χRΔq+χRq)(x/t)|u(t,x)|2(d+1)/(d1)dx,\displaystyle\begin{aligned} A_{q,R}^{\prime}(t)&=t^{-1}\langle\mathrm{Op}^{w}((Va)(x/t,\xi))r(t),r(t)\rangle+2t^{-1}\Re\langle\mathrm{Op}^{w}((Va)(x/t,\xi))r(t),v(t)\rangle\\ &\quad-\frac{2}{d+1}t^{-1}\int(\chi_{R}\Delta q+\nabla\chi_{R}\cdot\nabla q)(x/t)|u(t,x)|^{2(d+1)/(d-1)}\,dx,\end{aligned}

with VaVa given by (7.12). The second term of (7.12) enters the first two terms of (7.18) only through integrals over the region {R|x/t|2R}\{R\leq|x/t|\leq 2R\} and has the derivative of χR\chi_{R}, which is O(R1)O(R^{-1}). So it has no contribution (for fixed time interval) as RR\to\infty.

The cross term (i.e., the second term on the right-hand side of (7.18)) can be controlled it by quadratic terms in rr and vv. Using (7.13) with f=rf=r and f=vf=v, we have

Opw(12χR(x/t)(2ξx/t)TD2q(x/t)(2ξx/t))r,r+2Opw(12χR(x/t)(2ξx/t)TD2q(x/t)(2ξx/t))r,v\displaystyle\langle\mathrm{Op}^{w}\big(\frac{1}{2}\chi_{R}(x/t)(2\xi-x/t)^{T}D^{2}q(x/t)(2\xi-x/t)\big)r,r\rangle+2\Re\langle\mathrm{Op}^{w}\big(\frac{1}{2}\chi_{R}(x/t)(2\xi-x/t)^{T}D^{2}q(x/t)(2\xi-x/t)\big)r,v\rangle
14χR(x/t)(Ltr)D2q(x/t)(Ltr)𝑑xCχR(x/t)(Ltv)D2q(x/t)(Ltv)𝑑xCt2.\displaystyle\geq\frac{1}{4}\int\chi_{R}(x/t)(L_{t}r)^{*}D^{2}q(x/t)(L_{t}r)\,dx-C\int\chi_{R}(x/t)(L_{t}v)^{*}D^{2}q(x/t)(L_{t}v)\,dx-Ct^{-2}.

Then we substitute this into (7.18) and integrate it in tt from SS to TT and send RR\to\infty. We have 66 6 The first term is nonnegative and can be discarded before taking the limit to avoid the issue of finiteness of that term.

(7.19) Aq(T)Aq(S)\displaystyle A_{q}(T)-A_{q}(S) CSTK(v,t)dttCSTdtt32d+1ST(Δq)(x/t)|u(t,x)|2(d+1)/(d1)dxdtt.\displaystyle\geq-C\int_{S}^{T}K(v,t)\frac{dt}{t}-C\int_{S}^{T}\frac{dt}{t^{3}}-\frac{2}{d+1}\int_{S}^{T}\int(\Delta q)(x/t)|u(t,x)|^{2(d+1)/(d-1)}\,dx\frac{dt}{t}.

Of the terms on the right-hand side, it is straightforward that the tail of the second tends to zero as S,TS,T\to\infty. We also proved above that the first one has tail tending to zero.

Now we consider the last term. We will use that for any sequence An,BnA_{n},B_{n}\to\infty, An<BnA_{n}<B_{n}, we have

(7.20) AnBn(Δq)(x/t)|r(t,x)|2(d+1)/(d1)𝑑xdtt0.\displaystyle\int_{A_{n}}^{B_{n}}\int(\Delta q)(x/t)|r(t,x)|^{2(d+1)/(d-1)}\,dx\frac{dt}{t}\to 0.

We postpone its proof to Proposition 7.5, and the conditions there are indeed satisfied by the construction in Proposition 7.3. On the other hand, Lemma 2.2 gives

(7.21) AnBn(Δq)(x/t)|v(t,x)|2(d+1)/(d1)𝑑xdtt\displaystyle\int_{A_{n}}^{B_{n}}\int(\Delta q)(x/t)|v(t,x)|^{2(d+1)/(d-1)}\,dx\frac{dt}{t} ΔqLAnv(t)L2(d+1)/(d1)2(d+1)/(d1)dtt0.\displaystyle\leq\|\Delta q\|_{L^{\infty}}\int_{A_{n}}^{\infty}\|v(t)\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}\frac{dt}{t}\to 0.

Since u=r+vu=r+v, we have |u|2(d+1)/(d1)Cd(|r|2(d+1)/(d1)+|v|2(d+1)/(d1))|u|^{2(d+1)/(d-1)}\leq C_{d}(|r|^{2(d+1)/(d-1)}+|v|^{2(d+1)/(d-1)}), so (7.20) and (7.21) give

(7.22) AnBn(Δq)(x/t)|u(t,x)|2(d+1)/(d1)𝑑xdtt0.\displaystyle\int_{A_{n}}^{B_{n}}\int(\Delta q)(x/t)|u(t,x)|^{2(d+1)/(d-1)}\,dx\frac{dt}{t}\to 0.

Finally, taking S=AnS=A_{n} and T=BnT=B_{n} in (7.19) and combining this with the preceding tail estimates gives

lim infn(Aq(Bn)Aq(An))0.\displaystyle\liminf_{n\to\infty}(A_{q}(B_{n})-A_{q}(A_{n}))\geq 0.

Now we prove (7.20) used above.

Proposition 7.5.

Let u+H1(d)u_{+}\in H^{1}(\mathbb{R}^{d}), and let v(t)v(t) and r(t)r(t) be as in (1.4). Choose θ,χ0,χ1,θ~C(d)\theta,\chi_{0},\chi_{1},\tilde{\theta}\in C^{\infty}(\mathbb{R}^{d}) with bounded derivatives such that: 0χ0,χ1,θ~10\leq\chi_{0},\chi_{1},\tilde{\theta}\leq 1, 0θ10\leq\theta\lesssim 1, and θ=0\theta=0 outside {χ0=1}\{\chi_{0}=1\}, θ~=1\tilde{\theta}=1 near the supports of χ0,χ0,Δχ0\chi_{0},\nabla\chi_{0},\Delta\chi_{0}, and χ1=1\chi_{1}=1 near suppθ~\mathrm{supp}\tilde{\theta}. Assume yχ0(y)y\cdot\nabla\chi_{0}(y) is bounded. If AnA_{n}\to\infty, An<BnA_{n}<B_{n}, and

ϵn:=supAntBnχ1(x/t)r(t)H10,\displaystyle\epsilon_{n}:=\sup_{A_{n}\leq t\leq B_{n}}\|\chi_{1}(x/t)r(t)\|_{H^{1}}\to 0,

then

(7.23) AnBnθ(x/t)|r(t,x)|2(d+1)/(d1)𝑑xdtt0.\displaystyle\int_{A_{n}}^{B_{n}}\int\theta(x/t)|r(t,x)|^{2(d+1)/(d-1)}\,dx\frac{dt}{t}\to 0.
Proof.

By definition we have

(7.24) (it+Δ+|u|4/(d1))r\displaystyle(i\partial_{t}+\Delta+|u|^{4/(d-1)})r =|u|4/(d1)v.\displaystyle=-|u|^{4/(d-1)}v.

First, we prove the localized smallness of |u|4/(d1)|u|^{4/(d-1)}. Interpolation between L2L^{2} and L2d/(d2)L^{2d/(d-2)}, followed by the Sobolev embedding H1L2d/(d2)H^{1}\hookrightarrow L^{2d/(d-2)}, gives

(7.25) fL2d/(d1)fL21/2fL2d/(d2)1/2CfH1.\displaystyle\|f\|_{L^{2d/(d-1)}}\leq\|f\|_{L^{2}}^{1/2}\|f\|_{L^{2d/(d-2)}}^{1/2}\leq C\|f\|_{H^{1}}.

Because χ1=1\chi_{1}=1 on a neighborhood of suppθ~\mathrm{supp}\tilde{\theta}, for every AntBnA_{n}\leq t\leq B_{n}, using (7.25), we have

θ~(x/t)|r(t)|4/(d1)Ld/2\displaystyle\|\tilde{\theta}(x/t)|r(t)|^{4/(d-1)}\|_{L^{d/2}} χ1(x/t)r(t)L2d/(d1)4/(d1)Cχ1(x/t)r(t)H14/(d1)Cϵn4/(d1).\displaystyle\leq\|\chi_{1}(x/t)r(t)\|_{L^{2d/(d-1)}}^{4/(d-1)}\leq C\|\chi_{1}(x/t)r(t)\|_{H^{1}}^{4/(d-1)}\leq C\epsilon_{n}^{4/(d-1)}.

Using Lemma 2.3, we have

v(t)L2(d+1)/(d1)0as t.\displaystyle\|v(t)\|_{L^{2(d+1)/(d-1)}}\to 0\qquad\text{as }t\to\infty.

The free flow is unitary on L2L^{2}, so v(t)L2=u+L2\|v(t)\|_{L^{2}}=\|u_{+}\|_{L^{2}}. Interpolating these two bounds gives

v(t)L2d/(d1)\displaystyle\|v(t)\|_{L^{2d/(d-1)}} u+L2(d1)/(2d)v(t)L2(d+1)/(d1)(d+1)/(2d)0.\displaystyle\leq\|u_{+}\|_{L^{2}}^{(d-1)/(2d)}\|v(t)\|_{L^{2(d+1)/(d-1)}}^{(d+1)/(2d)}\to 0.

Since |u|4/(d1)=|r+v|4/(d1)|r|4/(d1)+|v|4/(d1)|u|^{4/(d-1)}=|r+v|^{4/(d-1)}\lesssim|r|^{4/(d-1)}+|v|^{4/(d-1)} 77 7 In fact, in our case, 4d11\frac{4}{d-1}\leq 1, and we have |r+v|4/(d1)|r|4/(d1)+|v|4/(d1)|r+v|^{4/(d-1)}\leq|r|^{4/(d-1)}+|v|^{4/(d-1)}., we have

(7.26) δn:=supAntBnθ~(x/t)|u(t)|4/(d1)Ld/2Cϵn4/(d1)+suptAnv(t)L2d/(d1)4/(d1)0.\displaystyle\begin{aligned} \delta_{n}&:=\sup_{A_{n}\leq t\leq B_{n}}\|\tilde{\theta}(x/t)|u(t)|^{4/(d-1)}\|_{L^{d/2}}\leq C\epsilon_{n}^{4/(d-1)}+\sup_{t\geq A_{n}}\|v(t)\|_{L^{2d/(d-1)}}^{4/(d-1)}\to 0.\end{aligned}

We next multiply (7.24) by χ0(x/t)\chi_{0}(x/t). Since θ~=1\tilde{\theta}=1 on a neighborhood of suppχ0\mathrm{supp}\chi_{0},

θ~(x/t)χ0(x/t)=χ0(x/t).\displaystyle\tilde{\theta}(x/t)\chi_{0}(x/t)=\chi_{0}(x/t).

Then we have

(7.27) (it+Δ+θ~(x/t)|u(t,x)|4/(d1))(χ0(x/t)r(t,x))=it1((x/t)χ0(x/t))r(t,x)+2t1χ0(x/t)r(t,x)+t2(Δχ0)(x/t)r(t,x)χ0(x/t)|u(t,x)|4/(d1)v(t,x).\displaystyle\begin{aligned} (i\partial_{t}+\Delta+\tilde{\theta}(x/t)|u(t,x)|^{4/(d-1)})(\chi_{0}(x/t)r(t,x))=&-it^{-1}((x/t)\cdot\nabla\chi_{0}(x/t))r(t,x)+2t^{-1}\nabla\chi_{0}(x/t)\cdot\nabla r(t,x)\\ &\qquad+t^{-2}(\Delta\chi_{0})(x/t)r(t,x)-\chi_{0}(x/t)|u(t,x)|^{4/(d-1)}v(t,x).\end{aligned}

For every nonnegative integer kk such that 2kAn<Bn2^{k}A_{n}<B_{n}, define

In,k={t:2kAntmin(2k+1An,Bn)}.\displaystyle I_{n,k}=\{t\in\mathbb{R}:2^{k}A_{n}\leq t\leq\min(2^{k+1}A_{n},B_{n})\}.

These intervals cover every tt satisfying AntBnA_{n}\leq t\leq B_{n} and have disjoint interiors.

For such an interval, define

χ0rS(In,k)\displaystyle\|\chi_{0}r\|_{S(I_{n,k})} :=χ0(x/t)r(t)LtLx2(In,k)+χ0(x/t)r(t)Lt2Lx2d/(d2)(In,k).\displaystyle:=\|\chi_{0}(x/t)r(t)\|_{L_{t}^{\infty}L_{x}^{2}(I_{n,k})}+\|\chi_{0}(x/t)r(t)\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}.

Here S(In,k)S(I_{n,k}) denotes the Strichartz norm. By the choice of cut-offs, we know

χ1=1 on a neighborhood of suppχ0suppχ0suppΔχ0.\displaystyle\chi_{1}=1\quad\text{ on a neighborhood of }\mathrm{supp}\chi_{0}\cup\mathrm{supp}\nabla\chi_{0}\cup\mathrm{supp}\Delta\chi_{0}.

In particular,

χ0(x/(2kAn))r(2kAn)L2\displaystyle\|\chi_{0}(x/(2^{k}A_{n}))r(2^{k}A_{n})\|_{L^{2}} χ1(x/(2kAn))r(2kAn)L2ϵn.\displaystyle\leq\|\chi_{1}(x/(2^{k}A_{n}))r(2^{k}A_{n})\|_{L^{2}}\leq\epsilon_{n}.

We now estimate the three commutator terms in (7.27). Since yχ0(y)y\cdot\nabla\chi_{0}(y) is bounded by hypothesis and is supported where χ00\nabla\chi_{0}\neq 0,

((x/t)χ0(x/t))r(t)L2\displaystyle\|((x/t)\cdot\nabla\chi_{0}(x/t))r(t)\|_{L^{2}} Cχ1(x/t)r(t)L2Cϵn\displaystyle\leq C\|\chi_{1}(x/t)r(t)\|_{L^{2}}\leq C\epsilon_{n}

for AntBnA_{n}\leq t\leq B_{n}.

On a neighborhood of suppχ0\mathrm{supp}\nabla\chi_{0}, one has both χ1=1\chi_{1}=1 and χ1=0\nabla\chi_{1}=0. Hence

χ0(x/t)r(t)=χ0(x/t)(χ1(x/t)r(t)),\displaystyle\nabla\chi_{0}(x/t)\cdot\nabla r(t)=\nabla\chi_{0}(x/t)\cdot\nabla(\chi_{1}(x/t)r(t)),

and therefore

χ0(x/t)r(t)L2Cϵn.\displaystyle\|\nabla\chi_{0}(x/t)\cdot\nabla r(t)\|_{L^{2}}\leq C\epsilon_{n}.

The residual rr of (1.4) is uniformly bounded in H1H^{1}, because uu is uniformly bounded in H1H^{1} and the free flow preserves the H1H^{1}-norm of u+u_{+}. Thus

(Δχ0)(x/t)r(t)L2C.\displaystyle\|(\Delta\chi_{0})(x/t)r(t)\|_{L^{2}}\leq C.

Since In,kdttlog2\int_{I_{n,k}}\frac{dt}{t}\leq\log 2, In,kdtt212kAn\int_{I_{n,k}}\frac{dt}{t^{2}}\leq\frac{1}{2^{k}A_{n}}, the first three terms on the right-hand side of (7.27) satisfy

it1((x/t)χ0(x/t))r+2t1χ0(x/t)r+t2(Δχ0)(x/t)rLt1Lx2(In,k)\displaystyle\|-it^{-1}((x/t)\cdot\nabla\chi_{0}(x/t))r+2t^{-1}\nabla\chi_{0}(x/t)\cdot\nabla r+t^{-2}(\Delta\chi_{0})(x/t)r\|_{L_{t}^{1}L_{x}^{2}(I_{n,k})} Cϵn+C2kAn.\displaystyle\leq C\epsilon_{n}+\frac{C}{2^{k}A_{n}}.

Using Hölder’s inequality and δn\delta_{n} defined in (7.26), we have

θ~(x/t)|u|4/(d1)χ0(x/t)rLt2Lx2d/(d+2)(In,k)\displaystyle\|\tilde{\theta}(x/t)|u|^{4/(d-1)}\chi_{0}(x/t)r\|_{L_{t}^{2}L_{x}^{2d/(d+2)}(I_{n,k})} θ~(x/t)|u(t)|4/(d1)LLd/2χ0(x/t)rLt2Lx2d/(d2)(In,k)δnχ0rS(In,k).\displaystyle\leq\|\tilde{\theta}(x/t)|u(t)|^{4/(d-1)}\|_{L^{\infty}L^{d/2}}\|\chi_{0}(x/t)r\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}\leq\delta_{n}\|\chi_{0}r\|_{S(I_{n,k})}.

Since χ0θ~\chi_{0}\leq\tilde{\theta} by our choice, we have

χ0(x/t)|u|4/(d1)vLt2Lx2d/(d+2)(In,k)\displaystyle\|\chi_{0}(x/t)|u|^{4/(d-1)}v\|_{L_{t}^{2}L_{x}^{2d/(d+2)}(I_{n,k})} δnvLt2Lx2d/(d2)(In,k).\displaystyle\leq\delta_{n}\|v\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}.

Applying (2.4) and (2.5) to (7.27) gives

χ0rS(In,k)\displaystyle\|\chi_{0}r\|_{S(I_{n,k})} C(ϵn+12kAn+δnχ0rS(In,k)+δnvLt2Lx2d/(d2)(In,k)).\displaystyle\leq C(\epsilon_{n}+\frac{1}{2^{k}A_{n}}+\delta_{n}\|\chi_{0}r\|_{S(I_{n,k})}+\delta_{n}\|v\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}).

Since δn0\delta_{n}\to 0 as nn\to\infty, the term containing χ0rS(In,k)\|\chi_{0}r\|_{S(I_{n,k})} on the right can be absorbed for all sufficiently large nn. Thus

(7.28) χ0rS(In,k)\displaystyle\|\chi_{0}r\|_{S(I_{n,k})} C(ϵn+12kAn+δnvLt2Lx2d/(d2)(In,k)).\displaystyle\leq C(\epsilon_{n}+\frac{1}{2^{k}A_{n}}+\delta_{n}\|v\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}).

Applying (2.4) with (q,r)=(2,2d/(d2))(q,r)=(2,2d/(d-2)) to vv, we have

kvLt2Lx2d/(d2)(In,k)2Cdu+L22.\displaystyle\sum_{k}\|v\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}^{2}\leq C_{d}\|u_{+}\|_{L^{2}}^{2}.

Now we prove (7.23). For every fH1f\in H^{1}, Hölder’s inequality and Sobolev embedding give

fL2(d+1)/(d1)2(d+1)/(d1)fL2d/(d1)4/(d1)fL2d/(d2)2CfH14/(d1)fL2d/(d2)2.\displaystyle\|f\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}\leq\|f\|_{L^{2d/(d-1)}}^{4/(d-1)}\|f\|_{L^{2d/(d-2)}}^{2}\leq C\|f\|_{H^{1}}^{4/(d-1)}\|f\|_{L^{2d/(d-2)}}^{2}.

Since χ0(x/t)r(t)\chi_{0}(x/t)r(t) is uniformly bounded in H1H^{1}, for all sufficiently large nn and AntBnA_{n}\leq t\leq B_{n}, we have

χ0(x/t)r(t)L2(d+1)/(d1)2(d+1)/(d1)Cχ0(x/t)r(t)L2d/(d2)2.\displaystyle\|\chi_{0}(x/t)r(t)\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}\leq C\|\chi_{0}(x/t)r(t)\|_{L^{2d/(d-2)}}^{2}.

Since θ=0\theta=0 outside {χ0=1}\{\chi_{0}=1\} and 0θ10\leq\theta\lesssim 1,

θ(x/t)|r(t,x)|2(d+1)/(d1)|χ0(x/t)r(t,x)|2(d+1)/(d1).\displaystyle\theta(x/t)|r(t,x)|^{2(d+1)/(d-1)}\lesssim|\chi_{0}(x/t)r(t,x)|^{2(d+1)/(d-1)}.

Therefore

In,kdθ(x/t)|r(t,x)|2(d+1)/(d1)𝑑xdtt\displaystyle\int_{I_{n,k}}\int_{\mathbb{R}^{d}}\theta(x/t)|r(t,x)|^{2(d+1)/(d-1)}\,dx\,\frac{dt}{t} C2kAnχ0(x/t)r(t)Lt2Lx2d/(d2)(In,k)2C2kAnχ0rS(In,k)2.\displaystyle\leq\frac{C}{2^{k}A_{n}}\|\chi_{0}(x/t)r(t)\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}^{2}\leq\frac{C}{2^{k}A_{n}}\|\chi_{0}r\|_{S(I_{n,k})}^{2}.

Using (7.28),

In,kdθ(x/t)|r(t,x)|2(d+1)/(d1)𝑑xdtt\displaystyle\int_{I_{n,k}}\int_{\mathbb{R}^{d}}\theta(x/t)|r(t,x)|^{2(d+1)/(d-1)}\,dx\,\frac{dt}{t} Cϵn22kAn+C(2kAn)3+Cδn22kAnvLt2Lx2d/(d2)(In,k)2.\displaystyle\leq\frac{C\epsilon_{n}^{2}}{2^{k}A_{n}}+\frac{C}{(2^{k}A_{n})^{3}}+\frac{C\delta_{n}^{2}}{2^{k}A_{n}}\|v\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}^{2}.

Summing over kk, and using

k12kAn\displaystyle\sum_{k}\frac{1}{2^{k}A_{n}} 2An,k1(2kAn)3CAn3,\displaystyle\leq\frac{2}{A_{n}},\qquad\sum_{k}\frac{1}{(2^{k}A_{n})^{3}}\leq\frac{C}{A_{n}^{3}},

we obtain

AnBndθ(x/t)|r(t,x)|2(d+1)/(d1)𝑑xdtt\displaystyle\int_{A_{n}}^{B_{n}}\int_{\mathbb{R}^{d}}\theta(x/t)|r(t,x)|^{2(d+1)/(d-1)}\,dx\,\frac{dt}{t} CAn1ϵn2+CAn3+Cδn2k12kAnvLt2Lx2d/(d2)(In,k)2.\displaystyle\leq CA_{n}^{-1}\epsilon_{n}^{2}+CA_{n}^{-3}+C\delta_{n}^{2}\sum_{k}\frac{1}{2^{k}A_{n}}\|v\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}^{2}.

For the last term on the right-hand side, we have

k12kAnvLt2Lx2d/(d2)(In,k)2\displaystyle\sum_{k}\frac{1}{2^{k}A_{n}}\|v\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}^{2} An1kvLt2Lx2d/(d2)(In,k)2CdAn1u+L22.\displaystyle\leq A_{n}^{-1}\sum_{k}\|v\|_{L_{t}^{2}L_{x}^{2d/(d-2)}(I_{n,k})}^{2}\leq C_{d}A_{n}^{-1}\|u_{+}\|_{L^{2}}^{2}.

Consequently,

AnBndθ(x/t)|r(t,x)|2(d+1)/(d1)𝑑xdtt\displaystyle\int_{A_{n}}^{B_{n}}\int_{\mathbb{R}^{d}}\theta(x/t)|r(t,x)|^{2(d+1)/(d-1)}\,dx\,\frac{dt}{t} C(An1ϵn2+An3+An1δn2u+L22),\displaystyle\leq C(A_{n}^{-1}\epsilon_{n}^{2}+A_{n}^{-3}+A_{n}^{-1}\delta_{n}^{2}\|u_{+}\|_{L^{2}}^{2}),

where the right-hand side tends to 0 as nn\to\infty since AnA_{n}\to\infty, ϵn0\epsilon_{n}\to 0, and δn0\delta_{n}\to 0. This concludes the proof. ∎

Finally, we prove Theorem 1.2.

Proof of Theorem 1.2.

Recall F(s)=(νes,Jes)F(s)=(\nu_{e^{s}},J_{e^{s}}) defined in (7.1). Assume that F(s)F(s) does not converge as ss\to\infty. Let Ω\Omega_{\sharp} be the set of subsequential limits from Proposition 5.4, which is also the set of all limits of F(sn)F(s_{n}) along sequences sns_{n}\to\infty. By Proposition 5.4, Ω\Omega_{\sharp} is a nonempty compact connected set, and because FF does not converge, it contains more than one point.

Let μ,Ψ\mu_{*},\Psi be as in Proposition 7.1. By Proposition 7.2, the only element of Ω\Omega_{\sharp} whose first component is μ\mu_{*} is Z=(μ,y2μ)Z_{*}=(\mu_{*},\frac{y}{2}\mu_{*}). In addition, if sns_{n}\to\infty, F(sn)ZF(s_{n})\to Z_{*}, and after passing to a subsequence, F(sn+τ)(μ(τ),J(τ))F(s_{n}+\tau)\to(\mu(\tau),J(\tau)) for τ\tau in every bounded interval, then

(μ(τ),J(τ))=Z for every τ.\displaystyle(\mu(\tau),J(\tau))=Z_{*}\quad\text{ for every }\tau\in\mathbb{R}.

Applying Proposition 7.3 and writing

μ=jajδxj,\displaystyle\mu_{*}=\sum_{j}a_{j}^{*}\delta_{x_{j}^{*}},

we obtain disjoint neighborhoods UjU_{j} of the points xjx_{j}^{*}, a smooth convex function qq such that q\nabla q and yq(y)q(y)y\cdot\nabla q(y)-q(y) are bounded, and the cutoffs appearing in Propositions 6.3 and 7.4. On each UjU_{j}, we have q(y)q(y) given by (7.8). Define

Q=qdμ=Ψdμ.\displaystyle Q_{*}=\int q\,d\mu_{*}=\int\Psi\,d\mu_{*}.

By Proposition 5.3, all coefficients of delta functions here has a uniform lower bound. A argument by contradiction shows: for measures (μ,J)(\mu,J) that are sufficiently close to ZZ_{*}, these points has to be very close to at least one of xjx_{j}^{*} as well. This means that we can choose a sufficiently small fixed closed dJd_{J}-ball 𝔅\mathfrak{B} centered at ZZ_{*} such that every element (μ,J)Ω(\mu,J)\in\Omega_{\sharp} lying in 𝔅\mathfrak{B} satisfies

suppμjUj.\displaystyle\mathrm{supp}\mu\subset\bigcup_{j}U_{j}.

We can also choose 𝔅\mathfrak{B} small enough so that it does not contain the entire Ω\Omega_{\sharp}. Consider the elements of Ω\Omega_{\sharp} lying on 𝔅\partial\mathfrak{B}, which is a compact set. Since Proposition 5.8 places the support of all relevant measures in a fixed compact ball, so

Ψ𝑑μ\displaystyle\int\Psi d\mu

is continuous in (μ,J)(\mu,J) on Ω\Omega_{\sharp}.

Now we show that no element (μ,J)Ω𝔅(\mu,J)\in\Omega_{\sharp}\cap\partial\mathfrak{B} can satisfy

Ψ𝑑μ=Ψdμ.\displaystyle\int\Psi\,d\mu=\int\Psi\,d\mu_{*}.

Indeed, the uniqueness statement in Proposition 7.1 would give μ=μ\mu=\mu_{*}, and Proposition 7.2 would then give (μ,J)=Z(\mu,J)=Z_{*}, which is impossible because ZZ_{*} is the center of 𝔅\mathfrak{B}. Therefore compactness and the continuity in μ\mu give a fixed number δ0>0\delta_{0}>0 such that every (μ,J)Ω𝔅(\mu,J)\in\Omega_{\sharp}\cap\partial\mathfrak{B} satisfies

(7.29) Ψ𝑑μΨdμδ0.\displaystyle\int\Psi d\mu\leq\int\Psi d\mu_{*}-\delta_{0}.

Since ZΩZ_{*}\in\Omega_{\sharp}, choose αn\alpha_{n}\to\infty such that

F(αn)Z.\displaystyle F(\alpha_{n})\to Z_{*}.

For sufficiently large nn, F(αn)F(\alpha_{n}) lies inside 𝔅\mathfrak{B}. Because Ω\Omega_{\sharp} contains a point outside 𝔅\mathfrak{B}, there are arbitrarily late times at which FF lies outside 𝔅\mathfrak{B}. The map FF is continuous in dJd_{J} by Proposition 5.4, so let βn>αn\beta_{n}>\alpha_{n} be the first time after αn\alpha_{n} at which FF reaches 𝔅\partial\mathfrak{B}.

We now prove:

(7.30) βnαn.\displaystyle\beta_{n}-\alpha_{n}\to\infty.

Otherwise, after passing to a subsequence, we can find a constant h0h\geq 0 such that

βnαnh.\displaystyle\beta_{n}-\alpha_{n}\to h.

By Theorem 5.6, after passing to another subsequence, there is a subsequential limit family (μα(τ),Jα(τ))(\mu_{\alpha}(\tau),J_{\alpha}(\tau)), defined for every τ\tau\in\mathbb{R}, such that

F(αn+τ)(μα(τ),Jα(τ))\displaystyle F(\alpha_{n}+\tau)\to(\mu_{\alpha}(\tau),J_{\alpha}(\tau))

locally uniformly for τ\tau in bounded intervals. At τ=0\tau=0,

(μα(0),Jα(0))=limnF(αn)=Z.\displaystyle(\mu_{\alpha}(0),J_{\alpha}(0))=\lim_{n\to\infty}F(\alpha_{n})=Z_{*}.

So by Proposition 7.2, we have

(μα(τ),Jα(τ))=Zfor every τ.\displaystyle(\mu_{\alpha}(\tau),J_{\alpha}(\tau))=Z_{*}\qquad\text{for every }\tau\in\mathbb{R}.

In particular,

F(βn)=F(αn+(βnαn))(μα(h),Jα(h))=Z.\displaystyle F(\beta_{n})=F(\alpha_{n}+(\beta_{n}-\alpha_{n}))\to(\mu_{\alpha}(h),J_{\alpha}(h))=Z_{*}.

This is impossible because every F(βn)F(\beta_{n}) lies in 𝔅\partial\mathfrak{B}. This proves (7.30).

Applying Theorem 5.6 to the time sequence βn\beta_{n}, after passing to a subsequence,

(7.31) F(βn+τ)(μβ(τ),Jβ(τ))\displaystyle F(\beta_{n}+\tau)\to(\mu_{\beta}(\tau),J_{\beta}(\tau))

locally uniformly for τ\tau in bounded intervals. For every fixed τ0\tau\leq 0, equation (7.30) implies that, for all sufficiently large nn,

αnβn+τβn.\displaystyle\alpha_{n}\leq\beta_{n}+\tau\leq\beta_{n}.

By the definition of βn\beta_{n}, F(s)F(s) remains in 𝔅\mathfrak{B} for every s[αn,βn]s\in[\alpha_{n},\beta_{n}]. Therefore (μβ(τ),Jβ(τ))(\mu_{\beta}(\tau),J_{\beta}(\tau)) lies in 𝔅\mathfrak{B} for every τ0\tau\leq 0. Moreover, since βn+τ\beta_{n}+\tau\to\infty, each (μβ(τ),Jβ(τ))(\mu_{\beta}(\tau),J_{\beta}(\tau)) belongs to Ω\Omega_{\sharp}. By the choice of 𝔅\mathfrak{B},

suppμβ(τ)jUj(τ0).\displaystyle\mathrm{supp}\mu_{\beta}(\tau)\subset\bigcup_{j}U_{j}\qquad(\tau\leq 0).

At τ=0\tau=0, equation (7.31) and the fact that every F(βn)F(\beta_{n}) lies in 𝔅\partial\mathfrak{B} show that (μβ(0),Jβ(0))(\mu_{\beta}(0),J_{\beta}(0)) also lies in 𝔅\partial\mathfrak{B}.

Let qq be the smooth convex function as in Proposition 7.3. Define

(7.32) Q(τ)=qdμβ(τ).\displaystyle Q(\tau)=\int q\,d\mu_{\beta}(\tau).

Since qq is given by (7.8) on UjU_{j}, and the expression in (7.8) agrees with a tangent plane of the convex function Ψ\Psi, we have q(y)Ψ(y)q(y)\leq\Psi(y) on UjU_{j}. Since (μβ(τ),Jβ(τ))Ω(\mu_{\beta}(\tau),J_{\beta}(\tau))\in\Omega_{\sharp}, the maximizing property of μ\mu_{*} gives

Ψdμβ(τ)Ψdμ.\displaystyle\int\Psi\,d\mu_{\beta}(\tau)\leq\int\Psi\,d\mu_{*}.

Therefore, for every τ0\tau\leq 0,

Q(τ)=qdμβ(τ)Ψdμβ(τ)Ψdμ=Q.\displaystyle Q(\tau)=\int q\,d\mu_{\beta}(\tau)\leq\int\Psi\,d\mu_{\beta}(\tau)\leq\int\Psi\,d\mu_{*}=Q_{*}.

At τ=0\tau=0, (μβ(0),Jβ(0))(\mu_{\beta}(0),J_{\beta}(0)) lies in 𝔅\partial\mathfrak{B}, so (7.29) gives

Q(0)Ψdμβ(0)Ψdμδ0=Qδ0.\displaystyle Q(0)\leq\int\Psi\,d\mu_{\beta}(0)\leq\int\Psi\,d\mu_{*}-\delta_{0}=Q_{*}-\delta_{0}.

Thus

(7.33) Q(0)Qδ0,Q(τ)Qτ0.\displaystyle Q(0)\leq Q_{*}-\delta_{0},\quad Q(\tau)\leq Q_{*}\;\forall\tau\leq 0.

Fix τ0\tau\leq 0 and consider the times satisfying eαnteβn+τe^{\alpha_{n}}\leq t\leq e^{\beta_{n}+\tau}, which exist by (7.30) for sufficiently large nn. Take a sequence of times tnt_{n} satisfying

eαntneβn+τ.\displaystyle e^{\alpha_{n}}\leq t_{n}\leq e^{\beta_{n}+\tau}.

Then F(logtn)F(\log t_{n}) lies in 𝔅\mathfrak{B}. Every limit of a sequence of these states therefore belongs to Ω\Omega_{\sharp}, lies in 𝔅\mathfrak{B}, and has its mass support contained in jUj\bigcup_{j}U_{j}. The cutoff χ1\chi_{1} given by Proposition 7.3 is supported a positive distance away from jUj\bigcup_{j}U_{j}. Applying Proposition 6.3 with K=jUj¯K=\bigcup_{j}\overline{U_{j}} gives

(7.34) supeαnteβn+τχ1(x/t)r(t)H10.\displaystyle\sup_{e^{\alpha_{n}}\leq t\leq e^{\beta_{n}+\tau}}\|\chi_{1}(x/t)r(t)\|_{H^{1}}\to 0.

By (7.34), Proposition 7.4 can now be applied to these intervals and, with AqA_{q} defined in (7.10), gives:

(7.35) lim infn(Aq(eβn+τ)Aq(eαn))0.\displaystyle\liminf_{n\to\infty}\big(A_{q}(e^{\beta_{n}+\tau})-A_{q}(e^{\alpha_{n}})\big)\geq 0.

For a pair (μ,J)(\mu,J) of measures as before, define

(7.36) Lq(μ,J)\displaystyle L_{q}(\mu,J) :=q(y)dJ(y)12(yq(y)q(y))𝑑μ(y).\displaystyle:=\int\nabla q(y)\cdot dJ(y)-\frac{1}{2}\int\big(y\cdot\nabla q(y)-q(y)\big)\,d\mu(y).

By this definition and the definition of (νt,Jt)(\nu_{t},J_{t}) in (1.12),

Lq(νt,Jt)\displaystyle L_{q}(\nu_{t},J_{t}) =dq(x/t)j(r(t,x))𝑑x12d(xtq(x/t)q(x/t))|r(t,x)|2𝑑x.\displaystyle=\int_{\mathbb{R}^{d}}\nabla q(x/t)\cdot j(r(t,x))\,dx-\frac{1}{2}\int_{\mathbb{R}^{d}}\Big(\frac{x}{t}\cdot\nabla q(x/t)-q(x/t)\Big)|r(t,x)|^{2}\,dx.

Since u(t,x)=v(t,x)+r(t,x)u(t,x)=v(t,x)+r(t,x), combining the estimates in Proposition 2.9, we have

|Aq(t)Lq(νt,Jt)|\displaystyle|A_{q}(t)-L_{q}(\nu_{t},J_{t})| |q|dL(|v(t,x)||r(t,x)|+|r(t,x)||v(t,x)|)𝑑x\displaystyle\leq\|\nabla q\|_{L^{\infty}}\int_{\mathbb{R}^{d}}\big(|v(t,x)||\nabla r(t,x)|+|r(t,x)||\nabla v(t,x)|\big)\,dx
(7.37) +|yqq|dL|v(t,x)||r(t,x)|𝑑x0.\displaystyle\quad+\|y\cdot\nabla q-q\|_{L^{\infty}}\int_{\mathbb{R}^{d}}|v(t,x)||r(t,x)|\,dx\to 0.

By the choice of αn\alpha_{n} and (7.2),

(7.38) F(αn)=(νeαn,Jeαn)(μ,y2μ).\displaystyle F(\alpha_{n})=(\nu_{e^{\alpha_{n}}},J_{e^{\alpha_{n}}})\to\Big(\mu_{*},\frac{y}{2}\mu_{*}\Big).

Choose χRCc(d)\chi_{R}\in C_{c}^{\infty}(\mathbb{R}^{d}) satisfying 0χR10\leq\chi_{R}\leq 1, equal to one on BRB_{R}, and supported in B2RB_{2R}. For fixed RR, the functions χRq\chi_{R}\nabla q and χR(yqq)\chi_{R}(y\cdot\nabla q-q) extend continuously to XX, so (7.38) may be tested against them. Using Lemma 2.5, we know

|x|>Rt|r(t,x)|2𝑑x0, as R,\displaystyle\int_{|x|>Rt}|r(t,x)|^{2}\,dx\to 0,\text{ as }R\to\infty,

and this controls any L2L^{2}-level term as well. In addition

|x|>Rt|j(r(t,x))|𝑑x(|x|>Rt|r(t,x)|2𝑑x)1/2r(t)L2\displaystyle\int_{|x|>Rt}|j(r(t,x))|\,dx\leq\Big(\int_{|x|>Rt}|r(t,x)|^{2}\,dx\Big)^{1/2}\|\nabla r(t)\|_{L^{2}}

and the uniform H1H^{1} bound on rr control the current tails. The limit pair is supported in a fixed ball by Proposition 5.8. Thus, first letting nn\to\infty and then RR\to\infty, we obtain

(7.39) Lq(νeαn,Jeαn)Lq(μ,y2μ).\displaystyle L_{q}(\nu_{e^{\alpha_{n}}},J_{e^{\alpha_{n}}})\to L_{q}\Big(\mu_{*},\frac{y}{2}\mu_{*}\Big).

Combining (7.39) with (7.37) gives

(7.40) Aq(eαn)Lq(μ,y2μ).\displaystyle A_{q}(e^{\alpha_{n}})\to L_{q}(\mu_{*},\frac{y}{2}\mu_{*}).

Using (7.36), we have

Lq(μ,y2μ)\displaystyle L_{q}(\mu_{*},\frac{y}{2}\mu_{*}) =q(y)y2dμ(y)12(yq(y)q(y))dμ(y)=12qdμ=Q2.\displaystyle=\int\nabla q(y)\cdot\frac{y}{2}\,d\mu_{*}(y)-\frac{1}{2}\int\big(y\cdot\nabla q(y)-q(y)\big)\,d\mu_{*}(y)=\frac{1}{2}\int q\,d\mu_{*}=\frac{Q_{*}}{2}.

Similarly, (7.31) gives

Aq(eβn+τ)Lq(μβ(τ),Jβ(τ)).\displaystyle A_{q}(e^{\beta_{n}+\tau})\to L_{q}(\mu_{\beta}(\tau),J_{\beta}(\tau)).

Combining (7.35), (7.40), and (7.39), we obtain

(7.41) Lq(μβ(τ),Jβ(τ))Q2(τ0).\displaystyle L_{q}(\mu_{\beta}(\tau),J_{\beta}(\tau))\geq\frac{Q_{*}}{2}\qquad(\tau\leq 0).

The measures in (7.31) satisfy (5.29) from Theorem 5.6. Although (5.29) is stated for compactly supported smooth functions, Proposition 5.8 places the support of every μβ(τ)\mu_{\beta}(\tau) and Jβ(τ)J_{\beta}(\tau) in one fixed compact ball. We may therefore apply (5.29) to a compactly supported smooth function that agrees with qq on that ball. Equation (5.29) shows that QQ is absolutely continuous and gives, for almost every τ\tau,

(7.42) Q(τ)=q(y)(2dJβ(τ)(y)ydμβ(τ)(y)).\displaystyle Q^{\prime}(\tau)=\int\nabla q(y)\cdot\big(2\,dJ_{\beta}(\tau)(y)-y\,d\mu_{\beta}(\tau)(y)\big).

Using (7.42) and the definition of LqL_{q},

(7.43) 2Lq(μβ(τ),Jβ(τ))\displaystyle 2L_{q}(\mu_{\beta}(\tau),J_{\beta}(\tau)) =2q(y)dJβ(τ)(y)(yq(y)q(y))dμβ(τ)(y)=Q(τ)+Q(τ).\displaystyle=2\int\nabla q(y)\cdot dJ_{\beta}(\tau)(y)-\int\big(y\cdot\nabla q(y)-q(y)\big)\,d\mu_{\beta}(\tau)(y)=Q^{\prime}(\tau)+Q(\tau).

Combining (7.41) and (7.43) gives

(7.44) Q(τ)+Q(τ)Q(τ0).\displaystyle Q^{\prime}(\tau)+Q(\tau)\geq Q_{*}\qquad(\tau\leq 0).

Applying a Grönwall type argument, we have

Q(τ)eτ(Q(0)Q)+Q.\displaystyle Q(\tau)\leq e^{-\tau}(Q(0)-Q_{*})+Q_{*}.

So |Q(τ)||Q(\tau)| can be arbitrarily large by taking τ\tau to be very negative. But every μβ(τ)\mu_{\beta}(\tau) has the fixed total mass mm_{\sharp}, and Proposition 5.8 places its support in one fixed compact ball. Since qq is bounded on that ball,

Q(τ)=qdμβ(τ)\displaystyle Q(\tau)=\int q\,d\mu_{\beta}(\tau)

is uniformly bounded for all τ\tau, which is a contradition.

So the assumption that F(s)=(νes,Jes)F(s)=(\nu_{e^{s}},J_{e^{s}}) does not converge is false. Since t=est=e^{s}, this is exactly the joint weak convergence of (νt,Jt)(\nu_{t},J_{t}) as tt\to\infty.

Denote the joint weak limit by (μ,J)(\mu_{\infty},J_{\infty}) and μ\mu_{\infty} takes the form in Proposition 5.9. Apply Theorem 5.6 to a sequence sns_{n}\to\infty. Since esn+τe^{s_{n}+\tau}\to\infty for every fixed τ\tau, the convergence just proved and (5.28) show that the resulting subsequential limit family is constant and equal to (μ,J)(\mu_{\infty},J_{\infty}). The support assertion in Theorem 5.6 therefore gives J==1LbδyJ_{\infty}=\sum_{\ell=1}^{L}b_{\ell}\delta_{y_{\ell}} for some bdb_{\ell}\in\mathbb{R}^{d}. Substituting these representations into (5.29) gives

(7.45) =1L(2bay)ϕ(y)=0\displaystyle\sum_{\ell=1}^{L}(2b_{\ell}-a_{\ell}y_{\ell})\cdot\nabla\phi(y_{\ell})=0

for every real-valued ϕCc(d)\phi\in C_{c}^{\infty}(\mathbb{R}^{d}). So we have b=ay/2b_{\ell}=a_{\ell}y_{\ell}/2, which means J=yμ/2J_{\infty}=y\mu_{\infty}/2 as stated in (1.13). ∎

8. Sequential to uniform upgrading

We now come to our byproduct: the upgrading estimates in Theorem 1.3 and Theorem 1.4 in this section.

8.1. Sequential to uniform upgrading I: smallness

We prove Theorem 1.3 first. Our Theorem 1.2 established in the previous section already gives a very strong characterization of r(t)r(t), and the remaining work is not so difficult.

Proof of Theorem 1.3.

Suppose (1.15) fails. Then Proposition 3.4, since (1.14) now holds, shows that m=M(u)M(u+)>0m_{\sharp}=M(u)-M(u_{+})>0. By Theorem 1.2, we have

νtμ,μ==1Laδy,a>0,yd.\displaystyle\nu_{t}\rightharpoonup\mu_{\infty},\quad\mu_{\infty}=\sum_{\ell=1}^{L}a_{\ell}\delta_{y_{\ell}},\qquad a_{\ell}>0,\quad y_{\ell}\in\mathbb{R}^{d}.

By the definition of νt\nu_{t} in (1.12), the weak convergence gives

(8.1) χ1(x/Tn)r(Tn)L22=1Laχ1(y)2.\displaystyle\|\chi_{1}(x/T_{n})r(T_{n})\|_{L^{2}}^{2}\to\sum_{\ell=1}^{L}a_{\ell}\chi_{1}(y_{\ell})^{2}.

The left side tends to zero by (1.14), and we conclude

(8.2) χ1(y)=0(1L).\displaystyle\chi_{1}(y_{\ell})=0\qquad(1\leq\ell\leq L).

Applying Proposition 6.4, after potentially passing to a subsequence, we have sn>Tns_{n}>T_{n} and a constant m>0m_{*}>0 such that

χ0(x/sn)r(sn)L2m.\displaystyle\|\chi_{0}(x/s_{n})r(s_{n})\|_{L^{2}}\geq m_{*}.

In the same way as (8.1) (i.e., using νtμ\nu_{t}\rightharpoonup\mu_{\infty}), we have

χ0(x/sn)r(sn)L22=1Laχ0(y)2m2.\displaystyle\|\chi_{0}(x/s_{n})r(s_{n})\|_{L^{2}}^{2}\to\sum_{\ell=1}^{L}a_{\ell}\chi_{0}(y_{\ell})^{2}\geq m_{*}^{2}.

Thus χ0(y)0\chi_{0}(y_{\ell})\neq 0 for some \ell. Since suppχ0{y:χ1(y)=1}\mathrm{supp}\chi_{0}\Subset\{y:\chi_{1}(y)=1\}, this implies χ1(y)=1\chi_{1}(y_{\ell})=1, contradicting (8.2). ∎

8.2. Sequential to uniform upgrade II: below threshold

We prove Theorem 1.4 in this section.

Proof of Theorem 1.4.

Let r(t)r(t) and v(t)v(t) be as in (1.4), and let Θ\Theta be as in (1.6). For i=0,2i=0,2, define

mi=(2π)dχi(2ξ)2|u^+(ξ)|2𝑑ξ,ki=(2π)dχi(2ξ)2|ξ|2|u^+(ξ)|2𝑑ξ.\displaystyle m_{i}=(2\pi)^{d}\int\chi_{i}(2\xi)^{2}|\widehat{u}_{+}(\xi)|^{2}\,d\xi,\quad k_{i}=(2\pi)^{d}\int\chi_{i}(2\xi)^{2}|\xi|^{2}|\widehat{u}_{+}(\xi)|^{2}\,d\xi.

Take an arbitrary subsequence of tnt_{n}, still denoted by tnt_{n}. By Theorem 2.7, after passing to a further subsequence,

u(tn)=eitnΔu++j=1Jwj(xj,n)+en,en0in H1,|xj,nxk,n|(jk).\displaystyle u(t_{n})=e^{it_{n}\Delta}u_{+}+\sum_{j=1}^{J}w_{j}(\cdot-x_{j,n})+e_{n},\quad e_{n}\to 0\quad\text{in }H^{1},\quad|x_{j,n}-x_{k,n}|\to\infty\qquad(j\neq k).

After passing to a further subsequence, for every jj,

xj,ntnyjX=d{}.\displaystyle\frac{x_{j,n}}{t_{n}}\to y_{j}\in X=\mathbb{R}^{d}\cup\{\infty\}.

Set

aj={χ2(yj),yjd,0,yj=.\displaystyle a_{j}=\begin{cases}\chi_{2}(y_{j}),&y_{j}\in\mathbb{R}^{d},\\ 0,&y_{j}=\infty.\end{cases}

By (4.4), after passing to a further subsequence, the limits

Mseq:=limnM(χ2(x/tn)u(tn)),Eseq:=limnE(χ2(x/tn)u(tn))\displaystyle M_{\mathrm{seq}}:=\lim_{n\to\infty}M(\chi_{2}(x/t_{n})u(t_{n})),\quad E_{\mathrm{seq}}:=\lim_{n\to\infty}E(\chi_{2}(x/t_{n})u(t_{n}))

exist. Using the results referred after (4.6)-(4.9), we have

χ2(x/tn)wj(xj,n)ajwj(xj,n)H10.\displaystyle\|\chi_{2}(x/t_{n})w_{j}(\cdot-x_{j,n})-a_{j}w_{j}(\cdot-x_{j,n})\|_{H^{1}}\to 0.

Using the same referred results, we also have the analogue of (4.6)

(8.3) Mseq=limnM(χ2(x/tn)eitnΔu+)+j=1Jaj2M(wj).\displaystyle M_{\mathrm{seq}}=\lim_{n\to\infty}M(\chi_{2}(x/t_{n})e^{it_{n}\Delta}u_{+})+\sum_{j=1}^{J}a_{j}^{2}M(w_{j}).

Here we used Lemma 2.3 to see that the contribution of the radiation term to the nonlinear term in the energy tends to zero. Also, the (referred) argument proving (4.9) gives

(8.4) Eseq=limn12(χ2(x/tn)eitnΔu+)L22+j=1JE(ajwj).\displaystyle E_{\mathrm{seq}}=\lim_{n\to\infty}\frac{1}{2}\|\nabla(\chi_{2}(x/t_{n})e^{it_{n}\Delta}u_{+})\|_{L^{2}}^{2}+\sum_{j=1}^{J}E(a_{j}w_{j}).

For 0a10\leq a\leq 1, we have

(8.5) E(aw)=a2E(w)+d12(d+1)(a2a2(d+1)/(d1))wL2(d+1)/(d1)2(d+1)/(d1).\displaystyle E(aw)=a^{2}E(w)+\frac{d-1}{2(d+1)}\big(a^{2}-a^{2(d+1)/(d-1)}\big)\|w\|_{L^{2(d+1)/(d-1)}}^{2(d+1)/(d-1)}.

Since 2(d+1)/(d1)>22(d+1)/(d-1)>2, the second term in (8.5) is nonnegative. For wj0w_{j}\neq 0, Proposition 3.2 gives

E(wj)>0,M(wj)E(wj)Θ.\displaystyle E(w_{j})>0,\qquad M(w_{j})E(w_{j})\geq\Theta.

It follows from (8.5) that

E(ajwj)0\displaystyle E(a_{j}w_{j})\geq 0

for every jj. Thus every term on the right-hand sides of (8.3) and (8.4) is nonnegative.

Then we apply Proposition 2.4, applied to both χ2(x/tn)eitnΔu+\chi_{2}(x/t_{n})e^{it_{n}\Delta}u_{+} and its gradient to obtain

limnM(χ2(x/tn)eitnΔu+)=m2,limn12(χ2(x/tn)eitnΔu+)L22=k22.\displaystyle\lim_{n\to\infty}M(\chi_{2}(x/t_{n})e^{it_{n}\Delta}u_{+})=m_{2},\quad\lim_{n\to\infty}\frac{1}{2}\|\nabla(\chi_{2}(x/t_{n})e^{it_{n}\Delta}u_{+})\|_{L^{2}}^{2}=\frac{k_{2}}{2}.

Terms introduced by differentiating χ(x/tn)\chi(x/t_{n}) will tend to zero since they carry an extra tn1t_{n}^{-1} factor. Passing to the limit in (1.17) gives

MseqEseq(1δ)Θ.\displaystyle M_{\mathrm{seq}}E_{\mathrm{seq}}\leq(1-\delta)\Theta.

The nonnegativity in (8.3) and (8.4) therefore gives

(8.6) m2k22(1δ)Θ.\displaystyle\frac{m_{2}k_{2}}{2}\leq(1-\delta)\Theta.

Suppose that wj0w_{j}\neq 0 and yjsuppχ1y_{j}\in\mathrm{supp}\chi_{1} for some jj. The second inclusion in (1.16) gives aj=χ2(yj)=1a_{j}=\chi_{2}(y_{j})=1. Equations (8.3) and (8.4) imply

M(wj)Mseq,E(wj)Eseq.\displaystyle M(w_{j})\leq M_{\mathrm{seq}},\qquad E(w_{j})\leq E_{\mathrm{seq}}.

Hence

M(wj)E(wj)MseqEseq(1δ)Θ<Θ,\displaystyle M(w_{j})E(w_{j})\leq M_{\mathrm{seq}}E_{\mathrm{seq}}\leq(1-\delta)\Theta<\Theta,

contradicting Proposition 3.2. Therefore

wj0yjsuppχ1.\displaystyle w_{j}\neq 0\quad\implies\quad y_{j}\notin\mathrm{supp}\chi_{1}.

For every nonzero profile, either yjdsuppχ1y_{j}\in\mathbb{R}^{d}\setminus\mathrm{supp}\chi_{1} or yj=y_{j}=\infty. The argument proving (5.6), now with the cutoff χ1\chi_{1}, gives

χ1(x/tn)wj(xj,n)H10\displaystyle\|\chi_{1}(x/t_{n})w_{j}(\cdot-x_{j,n})\|_{H^{1}}\to 0

for every jj. Subtracting the radiation term from (2.26) gives

r(tn)=j=1Jwj(xj,n)+en.\displaystyle r(t_{n})=\sum_{j=1}^{J}w_{j}(\cdot-x_{j,n})+e_{n}.

Multiplication by χ1(x/tn)\chi_{1}(x/t_{n}) is uniformly bounded on H1H^{1}. Therefore

χ1(x/tn)r(tn)H1\displaystyle\|\chi_{1}(x/t_{n})r(t_{n})\|_{H^{1}} j=1Jχ1(x/tn)wj(xj,n)H1\displaystyle\leq\sum_{j=1}^{J}\|\chi_{1}(x/t_{n})w_{j}(\cdot-x_{j,n})\|_{H^{1}}
+CenH10.\displaystyle\quad+C\|e_{n}\|_{H^{1}}\to 0.

Since the initial subsequence was arbitrary, the subsequence criterion gives

(8.7) χ1(x/tn)r(tn)H10.\displaystyle\|\chi_{1}(x/t_{n})r(t_{n})\|_{H^{1}}\to 0.

The first inclusion in (1.16) and (8.7) verify the hypotheses of Theorem 1.3 for the cutoffs χ0,χ1\chi_{0},\chi_{1} and the sequence tnt_{n}. Therefore

(8.8) χ0(x/t)u(t)χ0(x/t)v(t)H10.\displaystyle\|\chi_{0}(x/t)u(t)-\chi_{0}(x/t)v(t)\|_{H^{1}}\to 0.

On the other hand, similar as above, we apply Proposition 2.4 to χ0(x/t)v(t)\chi_{0}(x/t)v(t) and its derivatves 88 8 Rigorously, Proposition 2.4 required H1H^{1}-regularity. But the first claim concerning the measure μtu0\mu_{t}^{u_{0}} only requires L2L^{2} membership of u0u_{0}. to obtain

M(χ0(x/t)v(t))m0,E(χ0(x/t)v(t))k02.\displaystyle M(\chi_{0}(x/t)v(t))\to m_{0},\quad E(\chi_{0}(x/t)v(t))\to\frac{k_{0}}{2}.

The nonlinear term in the energy does not contribute to E(χ0(x/t)v(t))E(\chi_{0}(x/t)v(t)) in the limit by Lemma 2.3. Since the difference of the mass and energy can be bounded by the H1H^{1}-norm, combining the equation above and (8.8) shows

M(χ0(x/t)u(t))m0,E(χ0(x/t)u(t))k02.\displaystyle M(\chi_{0}(x/t)u(t))\to m_{0},\quad E(\chi_{0}(x/t)u(t))\to\frac{k_{0}}{2}.

By (1.16), we know 0χ0χ20\leq\chi_{0}\leq\chi_{2}. The terms involving the derivative of χi\chi_{i} does not contribute to the limit since they carry an extra t1t^{-1} factor. So (8.6) enforces

(8.9) m0k02(1δ)Θ,\displaystyle\frac{m_{0}k_{0}}{2}\leq(1-\delta)\Theta,

which proves (1.18) for all sufficiently large tt. Now we prove (1.19). Recalling the last part of (3.1), using (8.9) above, we have

limtχ0(x/t)u(t)L2(χ0(x/t)u(t))L2=m0k01δdQL2QL2<QL2QL2.\displaystyle\lim_{t\to\infty}\|\chi_{0}(x/t)u(t)\|_{L^{2}}\|\nabla(\chi_{0}(x/t)u(t))\|_{L^{2}}=\sqrt{m_{0}k_{0}}\leq\sqrt{\frac{1-\delta}{d}}\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}}<\|Q\|_{L^{2}}\|\nabla Q\|_{L^{2}}.

References

  • [1] T. Akahori and H. Nawa (2013) Blowup and scattering problems for the nonlinear Schrödinger equations. Kyoto Journal of Mathematics 53 (3), pp. 629–672. Cited by: §1.2.
  • [2] A. K. Arora, B. Dodson, and J. Murphy (2020) Scattering below the ground state for the 2d radial nonlinear Schrödinger equation. Proceedings of the American Mathematical Society 148 (4), pp. 1653–1663. External Links: 1906.00515 Cited by: §1.2.
  • [3] H. Berestycki and P.-L. Lions (1983) Nonlinear scalar field equations. I. Existence of a ground state. Arch. Rational Mech. Anal. 82 (4), pp. 313–345. External Links: ISSN 0003-9527, Document, Link, MathReview (Wei Ming Ni) Cited by: §3.
  • [4] T. Cazenave (2003) Semilinear Schrödinger equations. Courant Lecture Notes in Mathematics, Vol. 10, New York University, Courant Institute of Mathematical Sciences, New York; American Mathematical Society, Providence, RI. External Links: ISBN 0-8218-3399-5, Document, Link, MathReview (Woodford W. Zachary) Cited by: §2.1.
  • [5] B. Dodson and J. Murphy (2017) A new proof of scattering below the ground state for the 3d radial focusing cubic NLS. Proceedings of the American Mathematical Society 145 (11), pp. 4859–4867. External Links: 1611.04195 Cited by: §1.2.
  • [6] B. Dodson and J. Murphy (2018) A new proof of scattering below the ground state for the non-radial focusing NLS. Mathematical Research Letters 25 (6), pp. 1805–1825. External Links: 1712.09962 Cited by: §1.2, §1.2, Theorem 4.1.
  • [7] T. Duyckaerts, J. Holmer, and S. Roudenko (2008) Scattering for the non-radial 3D cubic nonlinear Schrödinger equation. Math. Res. Lett. 15 (6), pp. 1233–1250. External Links: ISSN 1073-2780, Document, Link, MathReview (Yoshihisa Nakamura) Cited by: §1.2, §1.2, §4.
  • [8] T. Duyckaerts and F. Merle (2009) Dynamic of threshold solutions for energy-critical NLS. Geom. Funct. Anal. 18 (6), pp. 1787–1840. External Links: ISSN 1016-443X,1420-8970, Document, Link, MathReview (Olivier J. Goubet) Cited by: §1.2.
  • [9] T. Duyckaerts and S. Roudenko (2010) Threshold solutions for the focusing 3D cubic Schrödinger equation. Rev. Mat. Iberoam. 26 (1), pp. 1–56. External Links: ISSN 0213-2230,2235-0616, Document, Link, MathReview (Yuichiro Kawahara) Cited by: §1.2.
  • [10] D. Fang, J. Xie, and T. Cazenave (2011) Scattering for the focusing energy-subcritical nonlinear Schrödinger equation. Sci. China Math. 54 (10), pp. 2037–2062. External Links: ISSN 1674-7283,1869-1862, Document, Link, MathReview (Yoshihisa Nakamura) Cited by: §1.2, §1.2, §1.2, Theorem 4.1.
  • [11] J. Gell-Redman, S. Gomes, and A. Hassell (2023) Scattering regularity for small data solutions of the nonlinear Schrödinger equation. arXiv preprint arXiv:2305.12429. Cited by: §1.2.
  • [12] J. Gell-Redman, S. Gomes, and A. Hassell (2025) Propagation estimates and Fredholm analysis for the time-dependent Schrödinger equation. Amer. J. Math. 147 (6), pp. 1577–1652. External Links: ISSN 0002-9327,1080-6377, MathReview Entry Cited by: §1.2.
  • [13] R. T. Glassey (1977) On the blowing up of solutions to the Cauchy problem for nonlinear Schrödinger equations. J. Math. Phys. 18 (9), pp. 1794–1797. External Links: ISSN 0022-2488,1089-7658, Document, Link, MathReview (A. A. Arsen\cprimeev) Cited by: §3.
  • [14] C. D. Guevara (2014) Global behavior of finite energy solutions to the dd-dimensional focusing nonlinear Schrödinger equation. Appl. Math. Res. Express. AMRX 2014 (2), pp. 177–243. External Links: ISSN 1687-1200,1687-1197, Document, Link, MathReview (Eduardo Colorado) Cited by: §1.2, §1.2, §1.2, §3, §3, Theorem 4.1.
  • [15] A. Hassell and Q. Jia (2026) The final state problem for the nonlinear Schrödinger equation in dimensions 1, 2 and 3. Pure Appl. Anal. 8 (1), pp. 1–39. External Links: ISSN 2578-5893,2578-5885, Document, Link, MathReview Entry Cited by: §5.
  • [16] J. Holmer and S. Roudenko (2008) A sharp condition for scattering of the radial 3D cubic nonlinear Schrödinger equation. Communications in Mathematical Physics 282 (2), pp. 435–467. External Links: Document Cited by: §1.2.
  • [17] C. E. Kenig and F. Merle (2006) Global well-posedness, scattering and blow-up for the energy-critical, focusing, non-linear Schrödinger equation in the radial case. Inventiones Mathematicae 166 (3), pp. 645–675. Cited by: §1.2.
  • [18] R. Killip and M. Visan (2010) The focusing energy-critical nonlinear Schrödinger equation in dimensions five and higher. Amer. J. Math. 132 (2), pp. 361–424. External Links: ISSN 0002-9327,1080-6377, Document, Link, MathReview (Benedetta Pellacci) Cited by: §1.2, §4.
  • [19] M. K. Kwong (1989) Uniqueness of positive solutions of Δuu+up=0\Delta u-u+u^{p}=0 in n\mathbb{R}^{n}. Archive for Rational Mechanics and Analysis 105 (3), pp. 243–266. External Links: Document Cited by: §1.1, §1.2, §3.
  • [20] K. Nakanishi and W. Schlag (2012) Global dynamics above the ground state energy for the cubic NLS equation in 3D. Calc. Var. Partial Differential Equations 44 (1-2), pp. 1–45. External Links: ISSN 0944-2669,1432-0835, Document, Link, MathReview (Tohru Ozawa) Cited by: §1.2.
  • [21] T. Tao (2004) On the asymptotic behavior of large radial data for a focusing non-linear Schrödinger equation. Dynamics of Partial Differential Equations 1 (1), pp. 1–48. External Links: math/0309428 Cited by: §1.2.
  • [22] T. Tao (2006) Nonlinear dispersive equations. CBMS Regional Conference Series in Mathematics, Vol. 106, Conference Board of the Mathematical Sciences, Washington, DC; by the American Mathematical Society, Providence, RI. Note: Local and global analysis External Links: ISBN 0-8218-4143-2, Document, Link, MathReview (Sebastian Herr) Cited by: §2.1.
  • [23] T. Tao (2007) A (concentration-)compact attractor for high-dimensional non-linear Schrödinger equations. Dyn. Partial Differ. Equ. 4 (1), pp. 1–53. External Links: ISSN 1548-159X,2163-7873, Document, Link, MathReview (W.-H. Steeb) Cited by: §1.2, §2.2, §2.2, Lemma 2.6.
  • [24] T. Tao (2008) A global compact attractor for high-dimensional defocusing non-linear Schrödinger equations with potential. Dynamics of Partial Differential Equations 5 (2), pp. 101–116. External Links: 0805.1544 Cited by: §1.2.
  • [25] T. Tao (2009) Why are solitons stable?. Bulletin of the American Mathematical Society (N.S.) 46 (1), pp. 1–33. External Links: Document Cited by: §1.2.
  • [26] M. I. Weinstein (1983) Nonlinear Schrödinger equations and sharp interpolation estimates. Communications in Mathematical Physics 87 (4), pp. 567–576. External Links: Document Cited by: §1.1, §1.2, §3.