In this paper, we first prove that the cubic, defocusing nonlinear Schrödinger equation on the two dimensional hyperbolic space with radial initial data in Hs(H2) is globally well-posed and scatters when s > 3/4. Then, we extend the result to nonlinearities of order p > 3. The result is proved by extending the high–low method of Bourgain in the hyperbolic setting and by using a Morawetz type estimate proved by Staffilani and Ionescu.

In this paper, we first consider the cubic nonlinear Schrödinger (NLS) initial value problem on the hyperbolic plane H2,

itu+ΔH2u=|u|2u,tR,xH2,u(0,x)=ϕ(x),
(1.1)

where u = u(t, x) is a complex-value function in spacetime R×H2 and ϕ is a radial initial datum.

The solution of (1.1) conserves both the mass

M(u(t)):=H2|u(t,x)|2dx=M(u0)
(1.2)

and the energy

E(u(t)):=H212|H2u(t,x)|2+14|u(t,x)|4dx=E(u0).
(1.3)

Conservation laws of mass and energy give the control of the L2 and 1 norms of the solutions, respectively.

Our goal in this paper is to prove the global well-posedness and scattering of (1.1) with the regularity of the initial data below H1.

In order to best frame the problem and to emphasize its challenges, we start by recalling the results in Rd, a setting that has been extensively considered in recent years. Consider the evolution equation in (1.4) with general non-linearities,

itu+Δu=|u|p1u,p>1,
(1.4)

in Rd. Let us first recall that the critical scaling exponent in Rd is

scd22p1.
(1.5)

It is well-known that in the sub-critical and critical regimes (s > sc and s = sc, respectively), the initial value problem (1.4) is locally well-posed;1–4 we recall here that with local well-posedness, we mean local in time existence, uniqueness, and continuous dependence of the data to the solution map. Thanks to the conservation laws of energy and mass, the H1-subcritical initial value problem and the L2-subcritical initial value problem are globally well-posed in the energy space H1 and mass space L2, respectively. The questions about scattering are much more delicate. What we mean with scattering here will be made more precise later (see, for example, Theorem 1.1), but in general terms, with scattering, we intend that the nonlinear solution as time goes to infinity approaches a linear one.

Before we talk about the global results in the more general subcritical case with data with regularity between L2 (mass) and H1 (energy), that is, in Hs, 0 < s < 1, let us denote with gMp the regularity index above which one obtains the global well-posedness for the NLS problem on the manifold M with power nonlinearity p and with sMp the index above which one obtains scattering (with the global well-posedness) again on the manifold M with power nonlinearity p.

The very first global well-posedness result in the subcritical case between the two (mass and energy) conservation laws (0 < s < 1) was given by Bourgain,5 where he developed the high–low method to prove global well-posedness for the cubic (p = 3) NLS in two dimensions for initial data in Hs, s > 3/5. According to the above notation, the regularity index in Ref. 5 is gR23=3/5.

We now describe the high–low method of Bourgain because it is the inspiration for part of our current work. To start, the initial datum is decomposed into a (smoother) low frequency part and a (rougher) high frequency part. The first step is to solve the NLS globally for the smoother part, for which the energy is finite, and then solve a difference equation for the rougher part. The miracle in this argument, which allows one to continue with an iteration, is that, in fact, the Duhamel term in the solution to the difference equation is small and smoother in an interval of time that is inversely proportional to the size of the low frequency part of the initial datum. At the next iteration, one merges this smoother part with the evolution of the low frequency part of the datum and repeats. It is worth mentioning that in order to obtain the miracle step, Bourgain5 used a Fourier transform based space Xs,b that captures particularly well the behavior of solutions with low regularity initial datum. Let us remark that there is no scattering result in the high–low method proposed by Bourgain.

Colliander–Keel–Staffilani–Takaoka–Tao6 improved the global well-posedness index gR23 of the initial data to 4/7 by introducing a different method, now known as the I-method. This is also based on an iterative argument. One first defines a Fourier multiplier that smooths out the initial data into the energy space and proves that the energy of the smoothed solution is almost conserved, that is, at each iteration, the growth of such modified energy is uniformly small. The index gR23 is derived by keeping the accumulation of energy controlled. As a result, Colliander et al.6 obtained a polynomial growth of the sub-energy Sobolev norm of the global solution. The cubic NLS in R3 was also considered in Ref. 6, and the index gR33=5/6.

Later, in Ref. 7, by combining the Morawetz estimate with the I-method and a bootstrapping argument, the same authors were able to lower the index gR33 to 4/5 and proved for the first time that the global solution also scatters, and hence, gR33=sR33=4/5. It is important to recall though that Bourgain8 had already proved global well-posedness for general data with index gR33=11/13 and scattering for radial data with index sR33=5/7.

To prove scattering, one needs to show that a spacetime norm of the solution is uniformly bounded. To this end, an iteration of local well-posedness would not suffice. Instead, one uses a Morawetz estimate that gives a uniform bound of the L4 spacetime norm of the solution, combined with the I-method. More in detail, one splits the timeline into a finite number of intervals Ij, of possible infinite length, on which the LIj4 of the solution is small. The smallness allows for a better spacetime bound of the global solution on each interval Ij, and then, one uses an iteration on the finite number of these intervals, which finally gives the desired spacetime uniform bound for the solution and hence scattering.

More results on the high–low method and the I-method both in Rd and compact manifolds can be found in Refs. 7, 921, and 22.

We now consider the initial value problem,

itu+ΔHdu=|u|p1u,tR,xHd,u(0,x)=ϕ(x),
(1.6)

with p > 1. Compared to what we recalled above, we expect even better results in Hd. In fact, the negative curvature of the ambient manifold allows for more dispersion in Hd than in the Euclidean spaces. Mathematically, we can see this in the Strichartz estimates on Hd, a family of estimates that is broader than the one obtained for the Euclidean space (see Refs. 36 and 37). The fact that the family of Strichartz estimates is larger in Hd reminds us of another case in which this is true. In fact also for the wave equation, the Strichartz estimates form a larger family. In this case though, it is not the curvature of the ambient manifold that generates a larger number of estimates, but instead it is the fact that the wave operator has a strong smoothing effect pointwise in time, a property that is not enjoyed by the Schrödinger operator. As a consequence, when one considers a nonlinear wave equation, the smoother and more plentiful estimates provide more suitable control of the nonlinear terms, and this is the reason why in the nonlinear wave setting, the miracle step in the high–low method in Ref. 23 does not need the Fourier type spaces Xs,b mentioned above. However, in contrast, the larger range of the Strichartz estimates for the Schrödinger operator in the hyperbolic space still is not readily enough to handle the miracle step since although one obtains better spacetime estimates, there is no pointwise smoothing effect; hence, the context we work in is more challenging than the one in Ref. 23. At this point, one may guess that using some hyperbolic version of the space Xs,b may do the trick. While this is indeed the case when the problem is posed in Td, see, for example, Ref. 11, in Hd the space, it is not clear how to define the Fourier transform based Xs,b type spaces in a way that is useful to handle nonlinearities. A naive definition using the Helgason–Fourier transform in Ref. 24 is deficient because of the following two reasons: first, the eigenfunctions of the Laplace–Beltrami operator on Hd lead to a very different Fourier inversion formula and Plancherel theorem. In particular, we cannot claim that the Fourier transform of a product is the convolution of Fourier transforms, which is a fundamental fact used in the estimates of nonlinear terms via the space Xs,b. Second, the frequency localization based on the Helgason–Fourier transform does not behave well in Lp(Hd), which causes difficulties in defining an effective Littlewood–Paley decomposition. We anticipate here that our approach to recover the miracle step, where one has to prove a gain of smoothness for the solution to the nonlinear difference equation, takes advantage of a Kato type smoothing effect. This smoothing is not pointwise in time, like for the wave operator, but average in time, hence much weaker. In order to make up for this weakness, we need to use a maximal function estimate combined with a better Sobolev embedding, which, in turn, forces us to assume radial symmetry for our initial data. We expect though that our global well-posedness and scattering results are true in general, and we believe that the more sophisticated smoothing effect in Ref. 25 may play an important role.

We now move to a summary of results that have been proved in the context of well-posedness and scattering for NLS in Hd. Although the initial value problem (1.6) cannot be properly scaled, we still use the same index sc defined in (1.5) to guide us in gauging the difficulty of proving the global well-posedness and scattering for (1.6). The subcritical initial value problem in the hyperbolic setting was first considered in Ref. 26, where the authors proved scattering for a family of power-type nonlinearity NLS with radial H1 data. Banica and Duyckaerts27 showed global well-posedness, scattering, and blow-up results for energy-subcritical focusing NLS also on the hyperbolic space. Recall that an NLS equation is called focusing when the nonlinearity in (1.4) has a negative sign, that is, itu+ΔH2u=|u|p1u.

In the critical setting Ionescu et al.28 proved global well-posedness and scattering of the energy-critical NLS in H3. This result uses an ad hoc profile decomposition technique to transfer the already available result in R3 (Ref. 29) into the H3 setting. A similar technique was also used in T3 for the same energy critical problem.30 We do not think that this method, which is well suited for critical settings, could work in our subcritical setting, when the initial data are in Hs(H2),0<s<1, but it may work to transfer in H2 the result that Dodson proved for mass critical in R2.31 To the best of the authors’ knowledge, there are no known subcritical global well-posedness and scattering results with initial data not at the conservation law level in hyperbolic spaces.

We now state the main result of this work for the initial value problem (1.1). Later, in Sec. V, we state a similar result for the more general version (1.6) with p > 3.

Theorem 1.1.
The initial value problem (1.1) with radial initial dataϕHs(H2)withs > 3/4 is globally well-posed and scattering holds, that is, there existsu±Hs(H2)such that
limt±u(t)eitΔH2u±Hxs(Hd)=0.
(1.7)

Remark 1.2.

Here, we conduct a discussion on the indices of regularity for the global well-posedness, and we make a comparison with other results.

As discussed above, the equivalent case we consider here but in R2 was treated by Bourgain without redial symmetry using the Xs,b space. Since we cannot use the same approach in H2, we decided first to rework this case using different tools such as Kato smoothing effect, maximal function estimates, and better Sobolev embedding. We did this because in R2, we have a Littlewood–Paley decomposition that works very well. Using these tools in the implementation of the high–low method, we obtained that the cubic radial NLS is globally well posed when s > 4/5, that is, gR23=4/5. Recall that Bourgain’s result gives gR23=3/5, which is better than what we can do in R2, and it is for general data. However, what we achieved in this first step is a blueprint that is generalizable to the H2 space. We did not report our calculations for R2 here, but the curious reader can check a longer version of this work in Ref. 32.

One notes that the index sH23=3/4 that we attained in Theorem 1.1 is smaller than the one we obtained in R2, where we worked out only the global well-posedness, not the scattering. This is because of the better radial Sobolev embedding in H2 and of the help coming from the strong Morawetz estimate used in the local theory.

Now, let us discuss a little bit of history concerning the indices of regularity for the global well-posedness. In Bourgain’s paper, where the high–low method was introduced,5 the global existence index is gR23=3/5. Later in Ref. 6, where the I-method was used, the index gR23 was improved to 4/7, and later in Ref. 33, thanks to a sophisticated treatment of the Fourier multiplier involved in the I-method, the global existence index gR23 was lowered further to 1/2. The global well-posedness of cubic NLS in two dimensions with H1/2 data was proved in Ref. 17. In addition, gR23 was improved to 1/3 in Ref. 9 and to 1/4 in Ref. 13. In Refs. 34, 35, and 31 Dodson proved global well-posedness and scattering for the mass-critical NLS in any dimension. As a consequence, via the persistence of regularity property, mass-critical NLS equations with any subcritical initial data are globally well-posed as well.

In this subsection, we summarize the three main parts of the proof of Theorem 1.1. In general terms, we combine the high–low method with a Morawetz type estimate that gives a bound for the spacetime L4 norm.

The first part of the proof deals with the analysis of the energy increment. Following Bourgain’s high–low method, we first decompose the initial datum into a high and low frequency part. Then, we write the solution u as the sum of the linear evolution of the high frequency part and a remainder ζ that solve a difference equation that evolves from the low frequency part of the original initial datum. In this first step, we assume that in an interval [0, τ], where τ could be infinity, the L4 spacetime norm of the solution is small. We then prove an estimate for the energy increment of ζ. This is the content of Proposition 3.1. To prove this energy increment estimate, we further decompose ζ = ζ1 + ζ2, where ζ1 is the nonlinear solution starting from the low frequency part of the datum and ζ2 solves a difference equation with zero datum. This part is similar to the high–low method of Bourgain, but here the interval of time is not small and the smallness comes from the L4 norm. The miracle step is then to be able to show that ζ2 is smoother and small in the appropriate norms. In the second part of the proof, we assume that the total L4 spacetime norm of the solution is bounded and we subdivide the timeline into finitely many intervals in which this norm is small. Here, we apply the first part described above, and we prove a global energy increment for ζ, this is Proposition 3.2. In the last part, we use a bootstrapping argument to show that indeed, the L4 norms of the solution u are bounded. This part requires a modification of the Morawetz estimate in Ref. 36 (see Proposition 4.1), and it uses the global energy increment proved in Proposition 3.2.

To summarize, the rest of this paper is organized as follows: In Sec. II, we discuss the geometry of the domain H2 and collect the useful analysis tools in H2. Next, in Sec. III, we present the calculation of the energy increment of the smoother part of the solution. In Sec. IV, we run a bootstrapping argument based on the estimates derived from Sec. III and modified Morawetz estimates and complete the Proof of Theorem 1.1. Finally, in Sec. V, we generalize this theorem for p > 3.

We define

fLtqLxr(I×H2):=IH2|f(t,x)|rdxqrdt1q,

where I is a time interval.

We use the Japanese bracket notation in the following sense:

ΩfX=fX+ΩfX,

where X is one of the normed spaces we use below.

We adopt the usual notation that AB or BA to denote an estimate of the form ACB for some constant 0 < C < depending only on the a priori fixed constants of the problem. If ABA, we write AB or AB.

We consider the Minkowski space R2+1 with the standard Minkowski metric,

(dx0)2+(dx1)2+(dx2)2,

and we define the bilinear form on R2+1×R2+1,

x,y=x0y0x1y1x2y2.

The hyperbolic space H2 is defined as

H2={xR2+1:x,x=1 and x0>0}.

An alternative definition for the hyperbolic space is

H2={x=(t,s)R2+1,(t,s)=(coshr,sinhrω),r0,ωS1}.

One has

dt=sinhrdr,ds=coshrωdr+sinhrdω

and the metric induced on H2 is dr2 + sinh2r dω2, where 2 is the metric on the sphere S1.

Then, one can rewrite integrals as

H2f(x)dx=0S1f(r,ω)sinhrdrdω.

Let 0 = {(1, 0, 0)} denote the origin of H2. The distance of a point to 0 is

d((coshr,sinhrω),0)=r.

More generally, the distance between two arbitrary points is

d(x,x)=cosh1([x,x]).

The general definition of the Laplace–Beltrami operator is given by

ΔH2=r2+coshrsinhrr+1sinh2rΔS1.

Remark 2.1.

The form of the Laplace–Beltrami operator implies that there will be no scaling symmetry in H2 as we usually have in the Rd setting.

In this subsection, we recall some important and classical analysis developed for the hyperbolic spaces.

1. Fourier transform on Hd

For θSd1 and λ a real number, the functions of the type

hλ,θ(x)=[x,Λ(θ)]iλd12,

where Λ(θ) denotes the point of Rd+1 given by (1, θ), are generalized eigenfunctions of the Laplace–Beltrami operator. Indeed, we have

ΔHdhλ,θ=λ2+(d1)24hλ,θ.

The Fourier transform on Hd is defined as

f̂(λ,θ)Hdhλ,θ(x)f(x)dx,

and the Fourier inversion formula on Hd takes the form of

f(x)=Sd1h¯λ,θ(x)f̂(λ,θ)dθdλ|c(λ)|2,

where c(λ) is the Harish–Chandra coefficient

1c(λ)2=12(2π)dΓ(iλ+d12)2Γ(iλ)2.

2. Strichartz estimates

In this subsection, we recall the Strichartz estimates proved in the hyperbolic space. We say that a couple (q, r) is admissible if (1/q, 1/r) belong to the triangle

Td=(1q,1r)(0,12]×(0,12)2q+drd2(0,12).

We have the following theorem by Anker–Piefelice.36 

Theorem 2.2
(Strichartz estimates in Ref. 37). Assume thatuis the solution to the inhomogeneous initial value problem,
itu+ΔHdu=F,tR,xHd,u(0,x)=ϕ(x).
Then, for any admissible exponents (q, r) and(q̃,r̃), we have the Strichartz estimates,
uLtqLxr(R×Hd)ϕLx2(Hd)+FLtq̃Lxr̃(R×Hd).

Remark 2.3.

Strichartz estimates are better in Hd in the sense that the set Td of admissible pairs for Hd is much wider than the corresponding set Id for Rd, which is just the lower edge of the triangle (see also Fig. 1).

FIG. 1.

Strichartz admissible pair regions (d = 2 and d ≥ 3) for the hyperbolic space.

FIG. 1.

Strichartz admissible pair regions (d = 2 and d ≥ 3) for the hyperbolic space.

Close modal

Definition 2.4
(Strichartz spaces). We define the Banach space
S0(I)=fC(I:L2(H2)):fS0(I)=sup(q,r) admissible fLtqLxr(I×H2)<.
In addition, we define the Banach space Sσ(I), where σ > 0,
Sσ(I)=fC(I:Hσ(H2)):fSσ(I)=(Δ)σ2fS0(I)<.

3. Local smoothing estimates in the hyperbolic space

Theorem 2.5
(Theorem 1.2 in Ref. 38: Local smoothing estimates in H2). For anyɛ > 0,
x12ε12eitΔfLt,x2(R×H2)fLx2(H2),x12ε0tei(ts)ΔF(s,x)dsLt,x2(R×H2)x12+εFLt,x2(R×H2).

Remark 2.6.

Kaizuka38 considered more general manifolds that they are denoted with X. To obtain the above theorem, one needs to take p(λ)=λ2, p(D)=ΔXρ2, and m = 2.

4. Heat-flow-based Littlewood–Paley projections and functional inequalities on H2

The Littlewood–Paley projections on H2 that we use in this paper are based on the linear heat equation esΔ. It turned out, in fact, that for us, this is a great substitute for the standard Littlewood–Paley decomposition used in Rd, since in Hd, one cannot localize in frequencies efficiently. We report below several results that first appeared in Ref. 25.

Definition 2.7
(Section 2.7.1 in Ref. 25: Heat-flow-based Littlewood–Paley projections). For any s > 0, we define
Psf=esΔf,Psf=s(Δ)esΔf.
By the fundamental theorem of calculus, it is straightforward to verify that
Psf=sPsfdss for s>0.
In particular, we have
f=0Psfdss,
which is the basic identity that relates f to its Littlewood–Paley resolution {Psf}s(0,). We also have
Psf=0sPsfdss.

Remark 2.8.

Intuitively, Psf may be interpreted as a projection of f to frequencies comparable to s−1/2. Ps and Ps can be viewed as the projections into low and high frequencies, respectively.

Lemma 2.9
(Lemma 2.5 in Ref. 25). Let 1 < p < andpq. Letρ0satisfy
0<ρ02<12min{1p,11p}.
ForfLxp(H2)ands > 0, we have
esΔfLxq(H2)+sΔesΔfLxq(H2)s(1p1q)eρ02sfLxp(H2).

Remark 2.10.
In particular, if p = q in Lemma 2.9, using Definition 2.7, we have
PsfLxp(H2)+PsfLxp(H2)fLxp(H2).

Lemma 2.11
(Corollary 2.7 in Ref. 25).Let 0 < α < 1 and 1 < p < . ForfLp(H2), we have
sα(Δ)αesΔfLxp(H2)fLxp(H2).

Lemma 2.12
(Lemma 2.9 in Ref. 25: Boundedness of Riesz transform). LetfC0(H2). Then, for 1 < p < , we have
fLxp(H2)(Δ)12fLxp(H2).

Lemma 2.13
(Lemma 2.10 in Ref. 25: Lp interpolation inequalities). LetfC0(H2). Then, for any 0 ≤ βαand 1 < p < , we have
(Δ)βfLxp(H2)fLxp(H2)1βα(Δ)αfLxp(H2)βα.
Moreover, for 1 < p < ,pq, we have
fLxq(H2)fLxp(H2)1θ(Δ)αfLxp(H2)θ, where 0<θ=1α(1p1q)<1.

Lemma 2.14
(Sobolev embedding).
Ws,p(Hd)Lp(Hd),if 1<pq< and 1p1q=sd.

Lemma 2.15
(Lemma 2.12 in Ref. 25: Gagliardo–Nirenberg inequality). LetfC0(H2). Then, for any 1 < p < ,pq, and 0 < θ < 1, we have
fLxq(H2)fLxp(H2)1θfLxp(H2)θ,where 1q=1pθ2.
In particular, for anys > 0,
fLx(H2)fLx4(H2)12fLx4(H2)12fLx2(H2)14fLx2(H2)12ΔfLx2(H2)14s12sΔfLx2(H2).

Lemma 2.16
(Proposition 2.14 in Ref. 25: Sobolev product rule). Forσ > 0, we have
fgHxσ(H2)fLx(H2)gHxσ(H2)+fHxσ(H2)gLx(H2).

Lemma 2.17
(General Sobolev product rule). Forσ > 0, we have
fgWxσ,r(H2)fWxσ,p1(H2)gLxp2(H2)+fLxq1(H2)gWxσ,q2(H2).

Remark 2.18.

Lemma 2.17 allows more possible Lp norms than Lemma 2.16 in the product rule. The Proof of Lemma 2.17 is using Triebel’s argument in Ref. 39 (see, for example, Secs. 7.2.2 and 7.2.4) and can also be found in Proposition 2.14 of Ref. 25. This proof relies on a localization lemma (see, for example, Lemma 2.16 in Ref. 25) to reduce to the standard Sobolev product rule.

Lemma 2.19
(Bernstein inequalities). For 0 ≤ β < α < β + 1,
(Δ)βPsfLx2(H2)sαβ(Δ)αfLx2(H2),(Δ)αPsfLx2(H2)sβα(Δ)βfLx2(H2).

Proof.
Using Definition 2.7, Lemma 2.11, and the formal property (−Δ)a(−Δ)b = (−Δ)a+b,
(Δ)βPsfLx2(H2)=0s(Δ)βt(Δ)etΔfdttLx2(H2)0stαβ1t1+βα(Δ)1+βαetΔ(Δ)αfLx2(H2)dt0stαβ1(Δ)αfLx2(H2)dt=sαβ(Δ)αfLx2(H2).
Similarly,
(Δ)αPsfLx2(H2)=sβαsαβ(Δ)αβesΔ(Δ)βfLx2(H2)sβα(Δ)βfLx2(H2).

5. Radial Sobolev embeddings

Lemma 2.20
(Lemma 2.13 in Ref. 25: Radial Sobolev embeddings in H2). For anys > 0 and any functionfradial,
sinh12(r)fLx(H2)fLx2(H2)12fLx2(H2)12s14sΔ12fLx2(H2).

Corollary 2.21
(Frequency localized radial Sobolev embeddings in H2). For anys > 0 andfradial,
sinh12(r)PsfLx(H2)s14Ps2fLx2(H2).

Proof.
Taking f = Psf in Lemma 2.20, we have
sinh12(r)PsfLx(H2)s14sΔ12PsfLx2(H2)=s14PsfLx2(H2)+(sΔ)12PsfLx2(H2).
Both terms in parentheses are bounded by Ps2fLx2(H2). In fact, by Definition 2.7, Remark 2.10, and Lemma 2.11,
PsfLx2(H2)=(sΔ)esΔfLx2(H2)=es2Δ(sΔes2Δf)Lx2(H2)Ps2fLx2(H2),(sΔ)12PsfLx2(H2)=(sΔ)12(sΔ)esΔfLx2(H2)=(sΔ)12es2Δ(sΔes2Δf)Lx2(H2)Ps2fLx2(H2).

Corollary 2.22.
For anys > 0, 1/4 < α < 1 andfradial,
sinh12(r)fLx(H2)fLx2(H2)114α(Δ)αfLx2(H2)14α.

Proof.
We write f into its Littlewood–Paley decomposition, and using Corollary 2.21, we have
sinh12(r)f=0sinh12(r)Psfdss0s14Ps2fLx2(H2)dss=0Ts14Ps2fLx2(H2)dss+Ts14Ps2fLx2(H2)dss:=I+II.
Here, T is a constant that will be chosen later. Using Remark 2.10 and Lemma 2.11, we have the following two estimates for PsfLxp(H2). For 0 < θ < 1,
PsfLx2(H2)=sΔesΔfLx2(H2)fLx2(H2),
(2.1)
PsfLx2(H2)=sθs1θ(Δ)1θesΔ(Δ)θfLx2(H2)sθ(Δ)θfLx2(H2).
(2.2)
Now, by (2.2) and (2.1), for 1/4 < α < 1,
I=0Ts14Ps2fLx2(H2)dss0Ts14sα(Δ)αfLx2(H2)dss=Tα14(Δ)αfLx2(H2),II=Ts14Ps2fLx2(H2)dss(Ts14dss)fLx2(H2)=T14fLx2(H2).
Therefore, for any T > 0,
sinh12(r)fLx(H2)Tα14(Δ)αfLx2(H2)+T14fLx2(H2).
Optimizing the choice of T, we obtain for 1/4 < α < 1,
sinh12(r)fLx(H2)(Δ)αfLx2(H2)14αfLx2(H2)114α.

In this section, we analyze a certain energy increment. As mentioned in the Introduction, we present a modified Morawetz type estimate in Subsection IV A and in Subsection IV B, we conclude the global well-posedness and scattering proof by showing that the spacetime L4 norm of the solution is uniformly bounded.

Let us recall schematically below the heat-flow-based Littlewood–Paley projections,

Psf=s(Δ)esΔf a projection to frequencies comparable to s12,Psf=esΔf=sPsfdss a projection to frequencies lower than s12,Psf=0sPsfdss a projection to frequencies higher than s12.

Now, we decompose the initial data ϕ into a low frequency component η0=P>s0ϕ and a high frequency component ψ0=Ps0ϕ, where s01 is a fixed large frequency and will be determined later in the proof. Note that s01/2 plays the same role as N0 in Ref. 5.

Using the above decomposition, we would like to write u into the sum of the following two solutions ψ and ζ, where ψ=eitΔPs0ϕ solves the linear Schrödinger with high frequency data,

itψ+ΔH2ψ=0,ψ(0,x)=ψ0=Ps0ϕ,
(3.1)

and ζ solves the difference equation with low frequency data,

itζ+ΔH2ζ=u2u=G(ζ,ψ),ζ(0,x)=η0=P>s0ϕ,
(3.2)

and here, G(ζ,ψ)=ζ+ψ2(ζ+ψ)=ζ2ζ+O(ζ2ψ)+O(ζψ2)+O(ψ3).

The main results in this section are a local energy increment (Proposition 3.1) and a conditional global energy increment (Proposition 3.2) for the solution ζ.

Proposition 3.1
(Local energy increment). Consideruas in (1.1) defined onI×H2, whereI = [0, τ], such that
uLt,x4(I×H2)4=ε
(3.3)
for some universal constantɛ. Then, fors > 3/4 and sufficiently smalls0, the solutionζ, under the decompositionu = ψ + ζdefined as in (3.1) and (3.2), satisfies the following energy increment:
E(ζ(τ))E(ζ(0))+Cs032s54.

Proposition 3.2
(Conditional global energy increment). Consideruas in (1.1) defined on[0,T]×H2, where
uLt,x4([0,T]×H2)4M
for some constantM. Then, fors > 3/4 and sufficiently smalls0, the energy ofζsatisfies the following energy increment:
E(ζ(T))E(ζ(0))+CMεs032s54,
whereɛis the small constant in Proposition 3.1.

Remark 3.3.

T could be infinity. In fact, the ultimate goal of this paper is to show that the spacetime L4 norm is bounded for all time intervals, which implies scattering.

To analyze the behavior of the solution ζ more carefully, we first make a further decomposition. That is, we would like to separate the differential equation (3.2) into a cubic NLS with low frequency data,

Itζ1+ΔH2ζ1=ζ12ζ1,ζ1(0,x)=η0=P>s0ϕ,
(3.4)

and a difference equation with zero initial value,

itζ2+ΔH2ζ2=u2uζ12ζ1,ζ2(0,x)=0.
(3.5)

Hence, ζ = ζ1 + ζ2, and the full solution u is the sum of these three solutions,

u=ζ1+ζ2+ψ.
(3.6)

It is worth mentioning that the decomposition in Bourgain’s work5 is a cubic NLS with low frequency data,

itu0+Δu0=u02u0,u0(0,x)=ϕ0(x)=P<N0ϕ,
(3.7)

and a difference equation with high frequency data,

itv+Δv=F(u0,v),v(0,x)=ψ0(x)=PN0ϕ,
(3.8)

where F(u0,v)=v2v+2u0v2+ū0v2+u02v¯+2u02v. Then, the full solution u is u = u0 + v. In our work, we need to be more careful.

Note that s01/2 plays the same role as N0 in Ref. 5. When comparing these two decomposition, we can relate them in the following sense: ζ1 is the same as u0 and ψ + ζ2 is, in fact, v, where ψ is the linear solution in (3.8) and ζ2 is the Duhamel term w in (3.8).

ζ1u0,ζ2w,ψeitΔPs0ϕ.

Step 1: Understanding the decomposed initial data. Recall that we decomposed the initial data ϕ = η0 + ψ0, where η0=P>s0ϕ and ψ0=Ps0ϕ. Here, we list several facts of the decomposed initial data η0 and ψ0.

Fact 3.4.
For the low frequency data η0,
(1)η0H1,η0Hx1(H2)s012(s1),(2)η0Hxσ(H2)s0σ2(s1) for 0<σ<1,(3)E(η0)s0s1.

In fact, (1) follows from Bernstein inequality (Lemma 2.19), (2) is by interpolating L2 and H1 norms, and (3) is due to Sobolev embedding (Lemma 2.14) and (2).

Fact 3.5.
For the high frequency data ψ0,
(1)ψ0Lx2(H2)s012s,(2)ψ0Hxs(H2)1.

Here, (1) follows from Bernstein inequality (Lemma 2.19), while (2) is due to the fact that ϕ being in Hs.

Step 2: Estimation on the solutionψof(3.1).

In fact, the solution ψ of the linear equation (3.1) is global, although it lives in a rough space Hs. Moreover, from the linear Strichartz estimates, Lemma 2.19 and (1) in Fact 3.5, one has

ψLt,x4(R×H2)=eitΔψ0Lt,x4(R×H2)ψ0Lx2(H2)s012s.
(3.9)

More generally,

ψS0(R)s012s and ψSσ(R)s012(sσ) for 0σs.
(3.10)

Step 3: Estimation on the solutionζ1of(3.4).

Lemma 3.6.
Due to the low frequency componentη0(x) ofϕbeing inH1,ζ1(t) is a global solution andζ1(t)Hx1(H2)is conserved. More precisely,
(1)ζ1(t) exists globally and E(ζ1)(t)=E(η0)s0s1,(2)ζ1Lt,x4(I×H2)4ε+s02s+ζ2Lt,x4(I×H2)4.

Remark 3.7.

Ultimately, we will show ζ1Lt,x4(I×H2)4ε in Corollary 3.9, and here, (2) is an intermediate step.

Proof of Lemma 3.6.
First, with the conservation of E(ζ1) and (3) in Fact 3.4, it is easy to see that
E(ζ1(t))E(η0)s0s1.
With η0 being in H1, thanks to Ref. 36, ζ1 is globally well-posed, which proves (1).
For (2), recall the decomposition of u in (3.6), then we simply use the triangle inequality, the assumption (3.3), and (3.9) and obtain
ζ1Lt,x4(I×H2)4uLt,x4(I×H2)4+ψLt,x4(I×H2)4+ζ2Lt,x4(I×H2)4ε+s02s+ζ2Lt,x4(I×H2)4.

Step 4: Estimation on the solutionζ2of(3.5)and extra estimates onζ1.

Recall (3.5)

itζ2+ΔH2ζ2=F(ζ1,ζ2,ψ),ζ2(0,x)=0,

where

F(ζ1,ζ2,ψ)=u2uζ12ζ1=O(ζ23)+O(ψ3)+O(ζ22ζ1)+O(ψ2ζ1)+O(ζ2ζ12)+O(ψζ12).

Lemma 3.8.
The solutionζ2satisfies the following estimates onI:
(1)ζ2Lt,x4(I×H2)s012s,(2)ζ2LtLx2(I×H2)s012s,(3)ζ2LtHx1(I×H2)s0s34.

Proof of Lemma 3.8.
Noting that (3.5) has zero initial value, we write out the integral equation using its Duhamel formula,
ζ2(t)=i0tei(ts)ΔH2Fds.
By Strichartz estimates and Hölder inequality, we have
ζ2Lt,x4(I×H2)ζ2Lt,x4(I×H2)3+ψLt,x4(I×H2)3+ζ2Lt,x4(I×H2)ζ1Lt,x4(I×H2)2+ψLt,x4(I×H2)ζ1Lt,x4(I×H2)2.
(3.11)
Note that here there should be two more nonlinear terms in F that contribute to FLt,x4/3 in (3.11), which are O(ζ22ζ1) and O(ψ2ζ1). However, we dropped them, since their contributions are controlled by a multiple of those of the four nonlinear terms that are written in (3.11). We will also drop them in the rest of this paper.
Using (3.9) and (2) in Lemma 3.6, we write (3.11) into
ζ2Lt,x4(I×H2)ζ2Lt,x4(I×H2)3+s032s+ζ2Lt,x4(I×H2)(ε12+s0s+ζ2Lt,x4(I×H2)2)+s012s(ε12+s0s+ζ2Lt,x4(I×H2)2).
Noting that initially ζ2(0) = 0, then by a continuity argument, we obtain
ζ2Lt,x4(I×H2)s012s,
which proves (1). The estimate in (3.11) also works for ζ2LtLx2(I×H2), and hence, (2) holds.

We postpone the Proof of (3) to step 6.□

With enough estimates on ζ2 in hand, as we promised in Lemma 3.6, we will finish the analysis of ζ1.

Corollary 3.9.

As a consequence of (1) in Lemma 3.8,

  1. we improve the bound in Lemma 3.6 byζ1Lt,x4(I×H2)4ε+s02sεand

  2. obtainζ1Sσ(I)s0σ2(s1), where 0 ≤ σ ≤ 1.

Proof of Corollary 3.9.

Combining (2) in Lemma 3.6 and (1) in Lemma 3.8, it is easy to see that (1) holds.

For (2), we use the integral equation corresponding to the initial value problem via the Duhamel principle, Strichartz estimates, and (1), and we obtain
ζ1S0(I)η0Lx2(H2)+ζ12ζ1Lt,x43(I×H2)1+ζ1Lt,x4(I×H2)31+ε34,ζ1S1(I)η0Lx2(H2)+ζ12ζ1Lt,x43(I×H2)s012(s1)+ζ1S1(I)ε12.
Therefore,
ζ1S0(I)1andζ1S1(I)s012(s1).
(3.12)
Then, interpolating (3.12), we have that for 0 < σ < 1,
ζ1Sσ(I)s0σ2(s1).

Step 5: Local energy increment.

Now, we are ready to compute the energy increment from 0 to τ and show such an increment is as described in Proposition 3.1. That is, we will show

E(ζ(τ))=E(ζ1(τ)+ζ2(τ))=E(ζ1(τ))+E(ζ1(τ)+ζ2(τ))E(ζ1(τ))E(η0)+Cs032s54.

In fact, a direct computation of the difference of the energy gives

E(ζ1(τ)+ζ2(τ))E(ζ1(τ))ζ1(τ)Hx1(H2)+ζ2(τ)Hx1(H2)ζ2(τ)Hx1(H2)+ζ1(τ)+ζ2(τ)3ζ2(τ)Lx1(H2):=I+II.

By the energy conservation of ζ1, (1) in Lemma 3.6, and (3) in Lemma 3.8,

I=ζ1(τ)Hx1(H2)+ζ2(τ)Hx1(H2)ζ2(τ)Hx1(H2)E(η0)12ζ2(τ)Hx1(H2)+ζ2(τ)Hx1(H2)2s012(s1)s0s34+s02(s34)s032s54,

for s > 1/2. [ζ1(τ)Hx1ζ2(τ)Hx1 dominates in I.]

Using Sobolev embedding, Lemma 2.13, (2) and (3) in Lemma 3.8, and (1) in Lemma 3.6, we have the following L4 norm estimates for ζ1 and ζ2:

ζ2(τ)Lx4(H2)ζ2(τ)Hx12(H2)ζ2(τ)Lx2(H2)12ζ2(τ)Hx1(H2)12s034s38,ζ1(τ)Lx4(H2)E(ζ1)14s014(s1).

Combining with Hölder inequality, we compute

II=ζ1(τ)+ζ2(τ)3ζ2(τ)Lx1(H2)ζ1(τ)Lx4(H2)+ζ2(τ)Lx4(H2)3ζ2(τ)Lx4(H2)(s014(s1)+s034s38)3s034s38s032s98

for s > 1/4. [ζ1(τ)Lx4(H2)3ζ2(τ)Lx4(H2) dominates in II, and I dominates in E(ζ(τ)) − E(ζ(0)).]

Now, we finish the calculation of the analysis of the energy increment in Proposition 3.1.

Step 6: Proof of (3) in Lemma 3.8.

Before proving (3) in Lemma 3.8, we first state the following lemma:

Lemma 3.10.
FortIdefined in (3.3), we have for 0 ≤ σs,
ζ2Sσ(I)s012(σs+sσ)s012(sσ).
(3.13)

Proof of Lemma 3.10.
For 0 < σ < s, by the integral equation, Bernstein inequality (Lemma 2.19), Lemma 2.17, Strichartz inequalities, and (3.3),
ζ2Sσ(I)Δσ2O(ζ23)Lt,x43(I×H2)+Δσ2O(ψ3)Lt,x43(I×H2)+Δσ2O(ζ2ζ12)Lt,x43(I×H2)+Δσ2O(ψζ12)Lt,x43(I×H2)ζ2Sσ(I)ζ2Lt,x4(I×H2)2+ψSσ(I)ψLt,x4(I×H2)2+ζ1Sσ(I)ζ1Lt,x4(I×H2)ζ2Lt,x4(I×H2)+ζ2Sσ(I)ζ1Lt,x4(I×H2)2+ζ2Sσ(I)ζ1Lt,x4(I×H2)2+ζ1Sσ(I)ζ1Lt,x4(I×H2)ψLt,x4(I×H2).
(3.14)
Using (1) in Lemma 3.8, (3.9), (3.10), and Corollary 3.9, we write (3.14) into
ζ2Sσ(I)ζ2Sσ(I)s0s+s012(sσ)s0s+s0σ2(s1)ε14s012s+ζ2Sσ(I)ε12+ζ2Sσ(I)ε12+s0σ2(s1)ε14s012s.
Then, we have
ζ2Sσ(I)s0σ2(s1)s012s=s012(σs+sσ).

Finally, we arrive at the Proof of (3) in Lemma 3.8.

Proof of (3) in Lemma 3.8.

In this step, we prove the smoothness of the solution ζ2 using the local smoothing estimate and the radial assumption of the initial data. In fact, this is the only place where the radial assumption is used, and all other steps work for all general data.

First, by Strichartz inequalities, we write
ζ2LtLx2(I×H2)O(ζ23)Lt,x43(I×H2)+O(ψ3)Lt,x43(I×H2)+O(ζ2ζ12)Lt,x43(I×H2)+O(ψζ12)Lt,x43(I×H2):=I+II+III+IV.
(3.15)
Before estimating IIV, we compute the following norms that are needed in the rest of this proof:

Claim 3.11.
For r=x,
(1)sinh12(r)ψLt,x(R×H2)s012s14,(2)sinh12(r)ζ1Lt,x(I×H2)s014(s1),(3)sinh12(r)ζ2Lt,x(I×H2)s012s14.

Proof of Claim 3.11.
Using the radial Sobolev embedding (Corollary 2.22), Strichartz estimates, and Fact 3.5,
sinh12(r)ψLt,x(R×H2)ψLtLx2(R×H2)114α(Δ)αψLtLx2(R×H2)14αψ0Lx2(H2)114αψ0Hx2α(H2)14αs012s×(114α)s012(s2α)×14α=s012s14.
By Corollaries 2.22 and 3.9,
sinh12(r)ζ1Lt,x(I×H2)ζ1LtLx2(I×H2)12ζ1LtLx2(I×H2)12s014(s1).
By Corollary 2.22 and Lemma 3.10,
sinh12(r)ζ2Lt,x(I×H2)ζ2LtLx2(I×H2)114α(Δ)αζ2LtLx2(I×H2)14αs012s×(114α)s012(2αs+s2α)×14α=s034s14.

Let us also recall some estimates from Subsection III A [(3.10), (3.9), Corollary 3.9, Lemmas 3.10, and 3.8],

ψSσ(R)s012(sσ) for 0σs,ψLt,x4(R×H2)s012s,ζ1Sσ(I)s0σ2(s1) for 0σ1,ζ1Lt,x4(I×H2)4ε,ζ2Sσ(I)s012(sσ) for 0σs,ζ2Lt,x4(I×H2)s012s.
(3.16)

Now, we continue working on (3.15). By Lemma 3.10,

I=O(ζ23)Lt,x43(I×H2)ζ2LtLx2(I×H2)ζ2Lt83Lx8(I×H2)2s0sζ2LtLx2(I×H2).

This term will be absorbed by the left-hand side of (3.15).

For II, we employ the local smoothing estimate. Since ψ is a linear solution, the linear version should be enough for this term. To implement the local smoothing estimate, we would like to introduce the weight x1/2ε1, where ɛ1 is a small positive number, and split out half derivative from the full gradient. Then, by the chain rule and Hölder inequality, we write

O(ψ3)Lt,x43(I×H2)x12ε112(12ψ)Lt,x2(I×H2)x12+ε1O(ψ2)Lt,x4(I×H2).
(3.17)

Now, we compute the two factors above separately. Using the linear local smoothing estimate (Lemma II C 3) and Lemma 2.19, we write the first factor into

x12ε112(12ψ)Lt,x2(I×H2)12Ps0ϕLx2(H2)s012(s12).
(3.18)

To estimate the second factor in (3.17), by Hölder inequality, Sobolev embedding, (3.10), (3.9), and Claim 3.11, we have

x12+ε1O(ψ2)Lt,x4(I×H2)χ{x1}O(ψ2)Lt,x4(I×H2)+χ{x>1}x12+ε1O(ψ2)Lt,x4(I×H2)ψLt,x8(I×H2)2+sinh12(r)ψLt,x(I×H2)ψLt,x4(I×H2)s0s12+s012s14s012ss0s12
(3.19)

for all s. Then, combining (3.18) and (3.19), we obtain

II(3.17)s012(s12)s0s12=s032s34.

For III and IV, we write them in a similar way as in (3.17),

IIIO(ζ2ζ12)Lt,x43(I×H2)x12ε112(12ζ2)Lt,x2(I×H2)x12+ε1O(ζ12)Lt,x4(I×H2)
(3.20)
+x12ε112(12ζ1)Lt,x2(I×H2)x12+ε1O(ζ1ζ2)Lt,x4(I×H2),IVO(ψζ12)Lt,x43(I×H2)x12ε112(12ψ)Lt,x2(I×H2)x12+ε1O(ζ12)Lt,x4(I×H2)
(3.21)
+x12ε112(12ζ1)Lt,x2(I×H2)x12+ε1O(ζ1ψ)Lt,x4(I×H2).

Noting that the estimates for ζ2 in Claim 3.11 and (3.16) are equally good or better than those for ψ, we will only present the estimation of IV. Our goal here is to prove that IV is bounded by s0s3/4. In fact, III is also bounded by s0s3/4.

We first start with the terms with no derivative and positive weights. By Hölder inequality, Sobolev embedding, Corollary 3.9, and Claim 3.11,

x12+ε1O(ζ12)Lt,x4(I×H2)ζ1Lt,x8(I×H2)2+sinh12(r)ζ1Lt,x(I×H2)ζ1Lt,x4(I×H2)s012(s1)+s014(s1)s012(s1)

for s < 1.

By Hölder inequality, Sobolev embedding, Corollary 3.9, (3.10), and Claim 3.11,

x12+ε1O(ζ1ψ)Lt,x4(I×H2)ζ1Lt,x12(I×H2)ψLt,x6(I×H2)+sinh12(r)ψLt,x(I×H2)ζ1Lt,x4(I×H2)s013(s1)s012(s13)+s012s14ε14s012s14

for s > 3/4.

Then, we focus on the terms with derivatives and negative weights.

Note that we have treated x12ε112(12ψ)Lt,x2(I×H2) in (3.18).

By Hölder inequality, Fact 3.4, Sobolev embedding, Corollary 3.9, and Claim 3.11,

x12ε112(12ζ1)Lt,x2(I×H2)12η0Lx2(H2)+x12+ε1ζ12ζ1Lt,x2(I×H2)s014(s1)+ζ1Lt,x6(I×H2)3+sinh12(r)ζ1Lt,x(I×H2)ζ1Lt,x4(I×H2)2s014(s1)+s016(s1)×3+s014(s1)ε12s012(s1)

for s < 1.

Therefore, continue from (3.20) and (3.21),

IIIIVs012(s12)s012(s1)+s012(s1)s012s14s0s34 for all s.

Combining all the terms IIV, we write (3.15) into

ζ2LtLx2(I×H2)I+II+III+IVs0sζ2LtLx2(I×H2)+s032s34+s0s34+s0s34.

Then, we have

ζ2LtLx2(I×H2)s0s34.

This concludes the Proof of (3) in Lemma 3.8 and completes the analysis of the energy increment.□

We divide the time interval [0, T] into [0, T] = ∪iIi = ∪i[ai, ai+1] such that on each Ii,

uLt,x4(Ii×H2)4=ε.
(3.22)

Hence,

#IiMε.

Let us remark that the length of such small intervals could be very long, and if some of them is an infinite interval, say [ak, ), then we just call ak+1 = .

On the first interval I1 = [0, a1], we can apply Proposition 3.1 and have the local energy increment

E(ζ(a1))E(ζ(0))+Cs032s54.

On the second time interval I2 = [a1, a2], we solve ζ by solving a cubic NLS with smoother data,

itζ1(1)+ΔH2ζ1(1)=ζ1(1)2ζ1(1),ζ1(1)(a1,x)=ζ(a1)=ζ1(a1)+ζ2(a1),
(3.23)

and a difference equation with zero initial value,

itζ2(1)+ΔH2ζ2(1)=u2uζ1(1)2ζ1(1),ζ2(1)(a1,x)=0.
(3.24)

Hence, ζ=ζ1(1)+ζ2(1), and the full solution will be u=ζ+ψ=ζ1(1)+ζ2(1)+ψ. Applying Proposition 3.1 again, we see that

E(ζ(a2))E(ζ(a1))+Cs032s54.

The reason why we are safe to apply Proposition 3.1 on I2 is that new decomposed initial data ζ(a1) and ψ(a1) satisfy all the facts in Facts 3.4 and 3.5, and all the calculations that we did in Proposition 3.1 will apply to the new systems (3.23) and (3.24). In particular, the size of new initial data ζ(a1) in energy is the size of ζ(0) in (3.2) plus a small error from ζ2(a1), which can be seen from Proposition 3.1,

E(ζ(a1))E(ζ(0))+Cs032s54small errorE(ζ(0))s0s1.

In addition, the H1 norm of ζ(a1) can be thought as the H1 norm of ζ(0) plus a small error,

ζ(a1)Hx1(H2)ζ1(a1)Hx1(H2)+ζ2(a1)Hx1(H2)s012(s1)size of H1 norm of ζ(0)+s0s34small errors012(s1).

Then, we can continue this iteration as long as the accumulated energy increment does not surpass the size of the initial energy of ζ(0), which guarantees that the setup for the smoother component remains the same size in the next iteration. That is,

CMεs032s54E(ζ(0))s0s1,

which gives

Ms012s+14.
(3.25)

The total energy increment is

E(ζ(T))E(ζ(0))+CMεs032s54.
(3.26)

For now, we finish the Proof of Proposition 3.2 and give the choice of s0.

Remark 3.12
(Boundedness of Hs norm of u). As a consequence of Proposition 3.2, we conclude that the Hs norm of u has the following bound:
u(T)Hxs(H2)s012(s1)s.
In fact, (3.26) implies the boundedness of the H1 norm of ζ(T),
ζ(T)Hx1(H2)2E(ζ(T))E(ζ(0))+CMεs032s54s0s1.
(3.27)
The triangle inequality and the mass conservation laws of u and ψ with (3.10) give the boundedness of L2 norm of ζ(T),
ζ(T)Lx2(H2)u(T)Lx2(H2)+ψ(T)Lx2(H2)u(0)Lx2(H2)+s012s2u(0)Lx2(H2).
(3.28)
Then, the Hs bound ζ(T) follows from the interpolation (3.28) and (3.27),
ζ(T)Hxs(H2)ζ(T)Lx2(H2)1sζ(T)Hx1(H2)ss012(s1)s.
(3.29)
Therefore, the Hs norm of u(T) is bounded due to (3.29) and the fact ψHs(H2),
u(T)Hxs(H2)ζ(T)Hxs(H2)+ψ(T)Hxs(H2)s012(s1)s+1s012(s1)s.
Consequently, we also have the bound of the Hσ (0 < σ < s) norm by interpolating the Hs with L2 norms,
u(T)Hxσ(H2)u(T)Hxs(H2)σsu(T)Lx2(H2)1σ2s0σ2(s1).

In this section, we use a bootstrapping argument to finally show the global well-posedness and scattering results stated in Theorem 1.1. In fact, we will show that the Morawetz norm is uniformly bounded, and once we have it, the proof of scattering will be standard.

We first recall that the Morawetz estimate of the cubic NLS on H2 in Ref. 36, when u is the solution to the cubic NLS equation itu+ΔH2u=u2u, reads as

uLt,x4([t1,t2]×H2)4uLtLx2([t1,t2]×H2)uLtHx1([t1,t2]×H2).

Proposition 4.1.
If we modify the NLS equation, that is,usolves
itu+ΔH2u=u2u+N,
then the modified Morawetz estimate becomes
uLt,x4([t1,t2]×H2)4uLtLx2([t1,t2]×H2)uLtHx1([t1,t2]×H2)+NūLt,x1([t1,t2]×H2)+Nu¯Lt,x1([t1,t2]×H2).
(4.1)

Remark 4.2.

The Proof of Proposition 4.1 is very similar as the Proof in Ref. 36. The difference is that we consider a more general nonlinear term, which mainly gives two extra terms that account for the two extra terms in (4.1).

Now, we are ready for the bootstrapping argument.

Step 1: Setup of the open–close argument.Define

W:=T:uLt,x4([0,T]×H2)4M,

where M > 0 is a constant. W is closed and non-empty. Now, we want to show that W is open. If T1W, then due to the local well-posedness theory and Remark 3.12, for some T0 > T1 and T0 sufficiently close to T1, we have

uLt,x4([0,T0]×H2)42M.

In fact, the Hs norm of u(T1) is bounded, then using a standard local well-posedness argument, we can continue the solution u from time T1 at least for a short time. Within such a short time period, due to the sub-criticality, the spacetime L4 norm of u will be bounded by twice the Hσ, for σ arbitrarily small, norm at T1, which is of order s0σ/2(s1) (see Remark 3.12). Hence, we want to ensure that

s0σ2(s1)<M14s014(12s+14),

which is achieved for any s > 3/4 by taking σ small enough (say, σ = 1/4). This guarantees the existence of such T0.

Now, we show that T0W, that is,

uLt,x4([0,T0]×H2)4M.
(4.2)

Recalling the decomposition of the solution u in (3.1) and (3.2) and using the smallness of L4 norm of ψ in (3.9), one can reduce (4.2) to

ζLt,x4([0,T0]×H2)412M.
(4.3)

Now, we will prove the improved bound of Lt,x4 in (4.3) in steps 2 and 3.

Step 2: Improving the bound forLt,x4. The modified Morawetz estimate in (4.1) that now gives

ζLt,x4([0,T0]×H2)4ζLtLx2([0,T0]×H2)ζLtHx1([0,T0]×H2)+Nζ¯Lt,x1([0,T0]×H2)+Nζ¯Lt1Lx1([0,T0]×H2).
(4.4)

Recall Eq. (3.2) that ζ satisfies, then in our case, N is given by

N=u2uζ2ζ=ψ+ζ2(ψ+ζ)ζ2ζ=ψ2ψ+O(ψ2ζ)+O(ψζ2).

Now, we estimate the right-hand side terms in (4.4).

For the second term in (4.4), by Hölder inequality, (3.9), and (4.3), we have

Nζ¯Lt,x1([0,T0]×H2)ψLt,x4([0,T0]×H2)3ζLt,x4([0,T0]×H2)+ψLt,x4([0,T0]×H2)ζLt,x4([0,T0]×H2)3s032sζLt,x4([0,T0]×H2)+s012sζLt,x4([0,T0]×H2)3s032sM14+s012sM34.

We write the last term in (4.4) as

Nζ¯Lt,x1([0,T0]×H2)O(ζ2ψ)ζ¯Lt,x1([0,T0]×H2)+O(ψ3)ζ¯Lt,x1([0,T0]×H2)ζLt,x4([0,T0]×H2)2ψLt,x4([0,T0]×H2)ζLt,x4([0,T0]×H2)+ψLt,x4([0,T0]×H2)3ζLt,x4([0,T0]×H2).
(4.5)

Claim 4.3.
We claim that
(1)ζLtLx2([0,T0]×H2)s012(s1),(2)ζLt,x4([0,T0]×H2)4Ms02(s1).

Assuming Claim 4.3, we continue the estimation of (4.5),

(4.5)M12s012s(Ms02(s1))14+s032s(Ms02(s1))14=M34s0s12+M14s02s12.

Using the same calculation as in (3.28), we have

ζ(t)Lx2(H2)u(t)Lx2(H2)+ψ(t)Lx2(H2)u(0)Lx2(H2)+s012s2u(0)Lx2(H2),

then the first term in (4.4) is bounded by

ζLtLx2([0,T0]×H2)ζLtHx1([0,T0]×H2)u(0)Lx2(H2)s012(s1).

Now, (4.4) becomes

ζLt,x4([0,T0]×H2)4s012(s1)+(s032sM14+s012sM34)+(M34s0s12+M14s02s12).

To close the argument, we need the following inequality holds for Ms01/2s+1/4:

s012(s1)+(s032sM14+s012sM34)+(M34s0s12+M14s02s12)<12M.
(4.6)

This requirement of (4.6) can be achieved for s > 3/4. Now, we are left to prove Claim 4.3.

Step 3: Proof of Claim 4.3.

Proof of Claim 4.3.
The total energy increment from Proposition 3.2 implies (1),
ζLtLx2([0,T0]×H2)2suptE(ζ(t))E(ζ(0))+Mεs032s54s0s1.
To estimate ζLt4Lx4([0,T0]×H2), we consider the subintervals Ii’s defined in Proposition 3.2. We claim that
ζLt,x4(Ii×H2)ζ(ai)Lx2(H2).
(4.7)
In fact, by Strichartz estimates,
ζLt,x4(Ii×H2)ζ(ai)Lx2(H2)+(ζ+ψ)3Lt,x43(Ii×H2)ζ(ai)Lx2(H2)+O(ζ3)Lt,x43(Ii×H2)+O(ψ3)Lt,x43(Ii×H2).
(4.8)
The second term in (4.8) will be absorbed by the left-hand side of (4.8),
O(ζ3)Lt,x43(Ii×H2)ζLt,x4(Ii×H2)ζLt,x4(Ii×H2)2ζLt,x4(Ii×H2)ε12.
The last inequality in the above equation is due to (3.9) and the construction of Ij (3.22),
ζLt,x4(Ii×H2)ψLt,x4(Ii×H2)+uLt,x4(Ii×H2)s012s+ε14ε14.
(4.9)
For the last term in (4.8), the same calculation as in (3.17) gives
O(ψ3)Lt,x43(Ii×H2)s032s34.
Then, (4.8) becomes
ζLt,x4(Ii×H2)ζ(ai)Lx2(H2)+ζLt,x4(Ii×H2)ε12+s032s34.
Therefore, the claim (4.7) follows.
Putting all the small intervals together and using (2), we get
ζLt,x4([0,T0]×H2)4MεsupIiζ(ai)Lx2(H2)4MεζLtLx2([0,T0]×H2)4Mεs02(s1).
Now, we finish the Proof of Claim 4.3.□

Our result for the cubic NLS in Theorem 1.1 can be generalized to a larger class of nonlinearities. In fact, we have the following result.

Theorem 5.1.

The initial value problem (1.6) with radial initial dataϕHs(H2)is globally well-posed and scatters inHs(H2)whens>3p63p5.

Note that the scaling of (1.6) is sc=12p1.

We will briefly present how the method presented above for the cubic NLS can be generalized to nonlinearities of order p.

  1. As in Sec. III, we decompose u in ψ and ζ, where ψ=eitΔPs0ϕ solves the linear Schrödinger with high frequency data and ζ solves the difference equation with low frequency data,
    itψ+ΔH2ψ=0,ψ(0,x)=ψ0=Ps0ϕ,itζ+ΔH2ζ=up1u,ζ(0,x)=η0=P>s0ϕ.
    (5.1)
    Then, using similar analysis, we obtain the global energy increment given the boundedness of u in the critical spacetime Lt,x2(p1) (see Propositions 5.2 and 5.3). The analogs of all the estimates that we used in the cubic case can be found in (5.3).
  2. Similarly, a bootstrapping argument on the Lt,x2(p1) norm gives both the global existence and scattering.

Note that the only difference in the general case is that the spacetime Lt,x2(p1) in the local theory (Proposition 5.2) is different from the Morawetz norm. Hence, in the bootstrapping argument, an intermediate step is needed. In fact, in this step, we first obtain and improve the estimates on the Morawetz norm, then we bootstrap the Lt,x2(p1) norm with the better Morawetz bound. Note that Lt,x2(p1) agrees with the Morawetz norm when p = 3, and hence, such a step is not needed in Sec. IV.

We now present the analogs of Propositions 3.1 and 3.2 on the energy increment.

Proposition 5.2
(Local energy increment). Consideruas in (1.6) defined onI×H2, whereI = [0, τ], such that
uLt,x2(p1)(I×H2)2(p1)=ε
for some universal constantɛ. Then, fors>pp+1and sufficiently smalls0, the solutionζ, under the decompositionu = ψ + ζdefined as in (5.1), satisfies the following energy increment:
E(ζ(τ))E(ζ(0))+Cs0p+34sp+24.

Proposition 5.3
(Conditional global energy increment). Consideruas in (1.6) defined on[0,T]×H2, where
uLt,x2(p1)([0,T]×H2)2(p1)M
for some constantM. Then, fors>pp+1and sufficiently smalls0, the energy ofζsatisfies the following energy increment:
E(ζ(T))E(ζ(0))+CMεs0p+34sp+24,
whereɛis the small constant in Proposition 5.2.

Within the proofs of above two propositions, we need a further decomposition for ζ as is (3.4)–(3.6),

itζ1+ΔH2ζ1=ζ1p1ζ1,ζ1(0,x)=η0=P>s0ϕ,itζ2+ΔH2ζ2=up1uζ1p1ζ1,ζ2(0,x)=0.
(5.2)

Hence, the full solution u is the sum of these three solutions u = ζ1 + ζ2 + ψ.

The analog of (3.16) can be computed similarly as follows:

ψSσ(R)s012(sσ) for 0σs,ψLt,x2(p1)(R×H2)s012(ssc),ζ1Sσ(I)ϕHxsc1 for 0σscs0σsc2(1sc)(s1) for scσ1,ζ1Lt,x2(p1)(I×H2)2(p1)ϕHxscε,
(5.3)
ζ2Sσ(I)ϕHxscs012(sσ) for 0σs,ζ2Lt,x2(p1)(I×H2)ϕHxscs012(ssc).

Claim 3.11 will be the same in the general setting. Most importantly, the analog of (3) in Lemma 3.8 is

ζ2LtHx1(I×H2)s0p+14sp4.

The choice of s0 in Proposition 5.3 is given by

Ms012(1s1sc12).
(5.4)

It is also worth mentioning that the hidden constant in (5.3) at the second iteration is bounded by the Hsc norm of ϕ plus a small error,

(ζ1+ζ2)(a1)Hxsc(H2)ζ1(a1)Hxsc(H2)+ζ2(a1)Hxsc(H2)ϕHxsc(H2)+s012ss012sc(12+1s1sc)

for s>pp+1. Then, the accumulated gain of this hidden constant will be dominated by the size of the Hsc norm of ϕ, hence not growing.

We consider

W:=T:uLt,x2(p1)([0,T]×H2)2(p1)M.

We are then reduced to showing that for T0 chosen in the same manner as in step 1 in Sec. IV,

ζLt,x2(p1)([0,T0]×H2)2(p1)12M.

Things are different here. First, interpolating ζLt,x2+([0,T0]×H2)2+M with the bound of the Lt,x2(p1) norm in the assumption gives an estimate on the Morawetz norm,

ζLt,xp+1([0,T0]×H2)p+1M.
(5.5)

Using (5.5) and the modified Morawetz estimate (4.1), we obtain as before

ζLt,xp+1([0,T0]×H2)p+1s012(s1).

If we simply require s01/2(s1)<M, here there will be no room to improve the Lt,x2(p1) norm at all. Hence, to this end, we demand it to be much smaller than M, that is, for α ∈ (0, 1),

s012(s1)MαM,
(5.6)

hence recalling (5.4), we get the first restriction on s

s>1α2+α(p1).

With this better Morawetz bound, we can improve the Lt,x2(p1) norm by making it smaller than M by Hölder inequality, (5.6), and Proposition 5.3 with the choice of M as in (5.4),

ζLt,x2(p1)([0,T0]×H2)2(p1)ζLt,xp+1([0,T0]×H2)p+1ζLt,x([0,T0]×H2)p3+Mα1ζLtLx2+([0,T0]×H2)p3+Mαs012(s1)(p3)M.

This with (5.4) implies the second restriction on s,

s>11α2(p3)+(1α)(p1).

Therefore, combining both conditions on s, we choose s>sH2p to be the best possible scattering index, where

sH2p=minα(0,1)max1α2+α(p1),11α2(p3)+(1α)(p1)=113p5=3p63p5.

G.S. and X.Y. graciously acknowledge the support from the Jarve Seed Fund. G.S. was also supported, in part, by the grant NSF (Grant No. DMS-1764403) and X.Y. by an AMS-Simons Travel Grants. G.S. and X.Y. would also like to thank A. Lawrie and S. Shahshahani for very insightful conversations.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
T.
Cazenave
,
Semilinear Schrödinger Equations
, Courant Lecture Notes in Mathematics Vol. 10 (
New York University, Courant Institute of Mathematical Sciences,
AMS
,
Providence, RI
,
2003
).
2.
T.
Cazenave
and
F. B.
Weissler
, “
The Cauchy problem for the nonlinear Schrödinger equation in H1
,”
Manuscripta Math.
61
,
477
494
(
1988
).
3.
T.
Cazenave
and
F.
Weissler
, “
Some remarks on the nonlinear Schrödinger equation in the subcritical case
,” in
New Methods and Results in Nonlinear Field Equations (Bielefeld, 1987)
, Lecture Notes in Physics Vol. 347 (
Springer
,
Berlin
,
1989
).
4.
T.
Cazenave
and
F. B.
Weissler
, “
The Cauchy problem for the critical nonlinear Schrödinger equation in Hs
,”
Nonlinear Anal.
14
,
807
836
(
1990
).
5.
J.
Bourgain
, “
Refinements of Strichartz’ inequality and applications to 2D-NLS with critical nonlinearity
,”
Int. Math. Res. Not.
1998
(
5
),
253
283
.
6.
J.
Colliander
,
M.
Keel
,
G.
Staffilani
,
H.
Takaoka
, and
T.
Tao
, “
Almost conservation laws and global rough solutions to a nonlinear Schrödinger equation
,”
Math. Res. Lett.
9
,
659
686
(
2002
).
7.
J.
Colliander
,
M.
Keel
,
G.
Staffilani
,
H.
Takaoka
, and
T.
Tao
, “
Global existence and scattering for rough solutions of a nonlinear Schrödinger equation on R3
,”
Commun. Pure Appl. Anal.
57
,
987
1014
(
2004
).
8.
J.
Bourgain
, “
Scattering in the energy space and below for 3D NLS
,”
J. Anal. Math.
75
,
267
297
(
1998
).
9.
J.
Colliander
and
T.
Roy
, “
Bootstrapped Morawetz estimates and resonant decomposition for low regularity global solutions of cubic NLS on R2
,”
Commun. Pure Appl. Anal.
10
(
2
),
397
414
(
2011
).
10.
D.
De Silva
,
N.
Pavlović
,
G.
Staffilani
, and
N.
Tzirakis
, “
Global well-posedness for the L2 critical nonlinear Schrödinger equation in higher dimensions
,”
Commun. Pure Appl. Anal.
6
(
4
),
1023
1041
(
2007
).
11.
D.
De Silva
,
N.
Pavlović
,
G.
Staffilani
, and
N.
Tzirakis
, “
Global well-posedness for a periodic nonlinear Schrödinger equation in 1D and 2D
,”
Discrete Contin. Dyn. Syst.
19
(
1
),
37
65
(
2007
).
12.
D.
De Silva
,
N.
Pavlović
,
G.
Staffilani
, and
N.
Tzirakis
, “
Global well-posedness and polynomial bounds for the defocusing L2-critical nonlinear Schrödinger equation in R
,”
Commun. Partial Differ. Equations
33
(
7-9
),
1395
1429
(
2008
).
13.
B.
Dodson
, “
Improved almost Morawetz estimates for the cubic nonlinear Schrödinger equation
,”
Commun. Pure Appl. Anal.
10
(
1
),
127
140
(
2011
).
14.
B. G.
Dodson
, “
Global well-posedness for the defocusing, quintic nonlinear Schrödinger equation in one dimension for low regularity data
,”
Int. Math. Res. Not.
2012
(
4
),
870
893
.
15.
B.
Dodson
, “
Global well-posedness and scattering for the defocusing, cubic nonlinear Schrödinger equation when n = 3 via a linear-nonlinear decomposition
,”
Discrete Contin. Dyn. Syst.
33
(
5
),
1905
1926
(
2013
).
16.
B.
Dodson
, “
Global well-posedness and scattering for nonlinear Schrödinger equations with algebraic nonlinearity when d = 2, 3 and u0 is radial
,”
Cambridge. J. Math.
7
(
3
),
283
318
(
2019
).
17.
Y. F.
Fang
and
M. G.
Grillakis
, “
On the global existence of rough solutions of the cubic defocusing Schrödinger equation in R2+1
,”
J. Hyperbolic Differ. Equations
04
,
233
257
(
2007
).
18.
C.
Guo
and
S.
Cui
, “
Global existence for 2D nonlinear Schrödinger equations via high-low frequency decomposition method
,”
J. Math. Anal. Appl.
324
(
2
),
882
907
(
2006
).
19.
Z.
Hani
, “
Global well-posedness of the cubic nonlinear Schrödinger equation on closed manifolds
,”
Commun. Partial Differ. Equations
37
(
7
),
1186
1236
(
2012
).
20.
Y.
Li
,
Y.
Wu
, and
G.
Xu
, “
Global well-posedness for the mass-critical nonlinear Schrödinger equation on T
,”
J. Differ. Equations
250
(
6
),
2715
2736
(
2011
).
21.
Q.
Su
, “
Global well posedness and scattering for the defocusing, cubic NLS in R3
,”
Math. Res. Lett.
19
(
2
),
431
451
(
2012
).
22.
N.
Tzirakis
, “
The Cauchy problem for the semilinear quintic Schrödinger equation in one dimension
,”
Differ. Integr. Equations
18
(
8
),
947
960
(
2005
).
23.
C.
Kenig
,
G.
Ponce
, and
L.
Vega
, “
Global well-posedness for semi-linear wave equations
,”
Commun. Partial Differ. Equations
25
(
9-10
),
1741
1752
(
2000
).
25.
A.
Lawrie
,
J.
Lührmann
,
S.-J.
Oh
, and
S.
Shahshahani
, “
Asymptotic stability of harmonic maps on the hyperbolic plane under the Schrödinger maps evolution
,” arXiv:1909.06899.
26.
V.
Banica
,
R.
Carles
, and
G.
Staffilani
, “
Scattering theory for radial nonlinear Schrödinger equations on hyperbolic space
,”
Geom. Funct. Anal.
18
(
2
),
367
399
(
2008
).
27.
V.
Banica
and
T.
Duyckaerts
, “
Global existence, scattering and blow-up for the focusing NLS on the hyperbolic space
,”
Dyn. Partial Differ. Equations
12
(
1
),
53
96
(
2015
).
28.
A.
Ionescu
,
B.
Pausader
, and
G.
Staffilani
, “
On the global well-posedness of energy-critical Schrödinger equations in curved spaces
,”
Anal. PDE
5
(
4
),
705
746
(
2012
).
29.
J.
Colliander
,
M.
Keel
,
G.
Staffilani
,
H.
Takaoka
, and
T.
Tao
, “
Global well-posedness and scattering for the energy-critical nonlinear Schrödinger equation in R3
,”
Ann. Math.
167
(
3
),
767
865
(
2008
).
30.
A. D.
Ionescu
and
B.
Pausader
, “
The energy-critical defocusing NLS on T3
,”
Duke Math. J.
161
(
8
),
1581
1612
(
2012
).
31.
B.
Dodson
, “
Global well-posedness and scattering for the defocusing, L2-critical, nonlinear Schrödinger equation when d = 2
,”
Duke Math. J.
165
(
18
),
3435
3516
(
2016
).
32.
G.
Staffilani
and
X.
Yu
, “
On the high-low method for NLS on the hyperbolic space
,” arXiv:2004.05711.
33.
J.
Colliander
,
M.
Keel
,
G.
Staffilani
,
H.
Takaoka
, and
T.
Tao
, “
Resonant decompositions and the I-method for the cubic nonlinear Schrödinger equation on R2
,”
Discrete Contin. Dyn. Syst.
21
,
665
686
(
2008
).
34.
B.
Dodson
, “
Global well-posedness and scattering for the defocusing, L2-critical nonlinear Schrödinger equation when d ≥ 3
,”
J. Am. Math. Soc.
25
(
2
),
429
463
(
2012
).
35.
B.
Dodson
, “
Global well-posedness and scattering for the defocusing, L2-critical, nonlinear Schrödinger equation when d = 1
,”
Am. J. Math.
138
(
2
),
531
569
(
2016
).
36.
A. D.
Ionescu
and
G.
Staffilani
, “
Semilinear Schrödinger flows on hyperbolic spaces: Scattering H1
,”
Math. Ann.
345
(
1
),
133
158
(
2009
).
37.
J.-P.
Anker
and
V.
Pierfelice
, “
Nonlinear Schrödinger equation on real hyperbolic spaces
,”
Ann. Inst. Henri Poincaré Anal. Non Linéaire
26
(
5
),
1853
1869
(
2009
).
38.
K.
Kaizuka
, “
Resolvent estimates on symmetric spaces of noncompact type
,”
J. Math. Soc. Jpn.
66
(
3
),
895
926
(
2014
).
39.
H.
Triebel
,
Theory of Function Spaces II
, Monographs in Mathematics Vol. 84 (
Birkhäuser Verlag
,
Basel
,
1992
).