The initial-value problem for the cubic-quintic NLS with non-vanishing boundary conditions
Rowan Killip, Jason Murphy, Monica Visan

TL;DR
This paper studies the initial-value problem for a cubic-quintic nonlinear Schrödinger equation with non-vanishing boundary conditions, proving scattering for small initial data in a weighted Sobolev space.
Contribution
It establishes a scattering result for the cubic-quintic NLS with non-zero boundary conditions, extending understanding of long-term behavior of solutions.
Findings
Proves scattering for small initial data in weighted Sobolev spaces.
Analyzes stability of constant solutions under the cubic-quintic NLS.
Extends scattering theory to non-vanishing boundary conditions.
Abstract
We consider the initial-value problem for the cubic-quintic NLS \[ (i\partial_t+\Delta)\psi=\alpha_1 \psi-\alpha_{3}\vert \psi\vert^2 \psi+\alpha_5\vert \psi\vert^4 \psi \] in three spatial dimensions in the class of solutions with as . Here , , and are such that is an energetically stable equilibrium solution to this equation. Normalizing the boundary condition to as , we study the associated initial-value problem for and prove a scattering result for small initial data in a weighted Sobolev space.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsAdvanced Mathematical Physics Problems
The initial-value problem for the cubic-quintic NLS with non-vanishing boundary conditions
Rowan Killip
Department of Mathematics, UCLA, Los Angeles, USA
,
Jason Murphy
Department of Mathematics and Statistics, Missouri University of Science and Technology
and
Monica Visan
Department of Mathematics, UCLA, Los Angeles, USA
Abstract.
We consider the initial-value problem for the cubic-quintic NLS
[TABLE]
in three spatial dimensions in the class of solutions with as . Here , , and are such that is an energetically stable equilibrium solution to this equation. Normalizing the boundary condition to as , we study the associated initial-value problem for and prove a scattering result for small initial data in a weighted Sobolev space.
1. Introduction
We study the initial-value problem for the cubic-quintic nonlinear Schrödinger equation (NLS) with non-vanishing boundary conditions in three space dimensions:
[TABLE]
We restrict attention to parameters for which the polynomial
[TABLE]
has a positive root with . This guarantees the energetic stability of the constant solution and we are interested in the dynamics of perturbations to this equilibrium. Correspondingly, we impose the boundary condition
[TABLE]
on our solutions to (1.1).
Rescaling both space-time and the values of the solution, we may reduce (1.1) and (1.2) to
[TABLE]
for some . This is the Hamiltonian evolution associated to the (conserved) energy
[TABLE]
In particular, we see that is always a local minimum of the energy. When , the energy is unbounded below which makes the system susceptible to wave collapse for large initial perturbations; however, we shall focus here on small perturbations. When , the potential energy has a local maximum at , while for , this energy has a second local minimum at .
The model (1.3) describes the behavior of a localized disturbance in an infinite expanse of quantum fluid that is otherwise quiescent. The particular case is known as the Gross–Pitaevskii equation and has received much attention. In particular, the asymptotic stability of the equilibrium solution to the Gross–Pitaevskii equation was proved in works of Gustafson, Nakanishi, and Tsai [10, 11, 12]. We will extend this result to the more general model (1.3), while also permitting a wider class of initial data.
The inclusion of an additional parameter in (1.3) has allowed researchers to better fit the behavior of several real physical systems and (1.1) has been used as a model in superfluidity [7, 8], descriptions of bosons [1] and of defectons [20], the theory of ferromagnetic and molecular chains [21, 22], and in nuclear hydrodynamics [14]. For comparison, the initial-value problem (1.1) with data decaying at infinity describes the dynamics of a finite body of fluid and it was studied in [18].
The boundary condition (1.2) may be further simplified to the following:
[TABLE]
Indeed, finite energy functions obeying (1.2) have a limiting phase as , which we can normalize to zero; furthermore, the dynamics of (1.1) preserve the value of this phase, so that the boundary condition is independent of time. See [5] for these observations in the case of the Gross–Pitaevskii equation. The fact that the boundary condition is independent of time breaks the gauge invariance of (1.1); in particular, we cannot use a time-dependent phase factor to remove the linear term in this equation. The linear term ultimately leads to weaker wave-like dispersion at low frequencies, which presents a key challenge in understanding the long-time behavior of solutions.
As we study perturbations of the constant solution , it is natural to introduce the function defined via . Using (1.3), we have the following equation for :
[TABLE]
where the nonlinearity is given by
[TABLE]
We may also write , where represents the polynomial of degree in and appearing in .
To put (1.6) into the standard framework of dispersive equations, it is convenient to diagonalize the equation. To do this, we employ operators and defined by
[TABLE]
The function then satisfies the following equation
[TABLE]
In our previous work [16] we considered the final-state problem for (1.6) (with ) and constructed solutions scattering to prescribed asymptotic states. Here we consider the initial-value problem for small localized data (and ).
Theorem 1.1**.**
Let with and . Suppose
[TABLE]
is sufficiently small. Then there exists a unique global solution to (1.7) with that scatters in both time directions, that is, there exist unique so that
[TABLE]
Furthermore, the convergence in (1.9) holds at a rate of for some .
We will outline the proof of this result in Section 1.1. Before turning to that subject, we will first make a series of remarks on the nature of this theorem and its relation to prior work. In particular, we will discuss the meaning of our theorem in the variables . Correspondingly, we define
[TABLE]
Remarks**.**
1. The norm appearing here occurs as the specialization to initial data of the more general -norm, defined in (1.12).
2. It was already known that our assumptions on imply global existence and uniqueness of solutions to (1.6). This follows from the well-posedness theory for (1.6) that we developed in our previous work [16]. Specifically, our arguments in [16] yield large-data global well-posedness for (1.6) in the energy space
[TABLE]
when and small-data global well-posedness for the remaining values of . Note that is not a vector space, but can naturally be interpreted as a complete metric space; see [5, 16]. By using (3.1) below, it is easy to see that under the hypotheses of Theorem 1.1, the initial data and so the previous analysis applies. This reasoning also shows that (1.9) implies that
[TABLE]
3. To properly recast the hypotheses on in terms of , it is natural to separate the high and low frequencies using Littewood–Paley operators; see Section 2 for the precise definitions. Self evidently, we have
[TABLE]
while setting in the proof of Lemma 3.3 allows one to deduce
[TABLE]
Analogously, the scattering statement can be recast as follows:
[TABLE]
as . Here .
As our hypotheses do not guarantee that nor are square-integrable, one should not expect that one can remove the operator from the last term above. On the other hand, due to the quantitative rate of convergence we obtain in Theorem 1.1, one can show that
[TABLE]
for some .
4. As the previous remark showed, our hypotheses on the real and imaginary parts of the physical field are different. While this may seem peculiar in a Schrödinger-like setting, it is perfectly normal in studies on the wave equation, where one invariably adopts differing norms for the displacement and velocity components. As the dispersion relation of the linear equation underlying (1.1) is wave-like at low frequencies (it has a conical singularity at the origin), it is natural that our hypotheses on , which is analogous to displacement, and those on , which is analogous to velocity, differ by exactly one derivative at low frequencies.
5. When , (1.3) becomes the Gross–Pitaevskii equation:
[TABLE]
The question of the long-time behavior of solutions to this equation was the central topic of a recent series of papers [10, 11, 12]. The last of these shows scattering in the three-dimensional setting (as considered here) and much of what we do follows closely in the footsteps of that paper. In actuality, the paper [12] has had a profound impact far beyond its original scope by codifying and popularizing (contemporaneously with [6]) what has become known as the space-time resonance method.
More concretely, the paper [12] proves scattering for solutions to (1.10) under the hypothesis that and belong to . These are stronger hypotheses than those employed in this paper in the sense that we have no weighted hypotheses on the derivatives beyond the angular regularity requirement . The nature of this improvement is most dramatic in the case of spherically symmetric data, for which .
6. In the special case of radial data, however, the contribution of this paper is eclipsed by the more recent work [9] on (1.10), which appeared while this paper was being finalized. Nevertheless, for general data, the hypotheses of the two works are incomparable. Specifically, [9] has weaker hypotheses on the low-frequency portion of the solution, since they do not require weighted decay in ; on the other hand, they do place more stringent hypotheses on the high frequencies, requiring that belongs to and not merely .
1.1. Strategy of the proof
The proof of Theorem 1.1 relies on a bootstrap argument. For this, we introduce a few specific norms. Letting denote any time interval, we first consider the Strichartz-type norm
[TABLE]
where the space-time norms are taken over .
Next, we introduce the vector field and the ‘energy norm’
[TABLE]
We refer to the component of this norm as the ‘weighted norm’.
For convenience, we also introduce the notation
[TABLE]
We will prove a bootstrap estimate of the form
[TABLE]
which, for small initial data, closes to give control in the -norm for all time. Using global -norm bounds, we can then prove scattering by standard arguments (see Section 11).
To explain our choice of spaces, we compare (1.7) with related scale-invariant nonlinear Schrödinger equations of the form
[TABLE]
The scaling symmetry of (1.13) identifies a scaling-critical space of initial data, namely , where . The regularity associated to a quintic nonlinearity is , while the regularity associated to a quadratic nonlinearity is . As (1.7) has a quintic nonlinear term, it is natural to seek control in a Strichartz space at regularity, namely, the -norm.
For Schrödinger equations with negative critical regularity, one seeks control of the solution in weighted spaces rather than Sobolev spaces of negative order, in which NLS is known to be illposed. From the point of view of scaling, prescribing is like prescribing . However, it is not expected that solutions will remain bounded in a weighted -norm. Indeed, this is not even case for the linear equation. Instead, in the case of (1.13), one endeavors to prove -bounds for the Galilean quantity , which in turn implies decay for the solution via Klainerman–Sobolev-type inequalities. As one can check that , we see that the operator defined above is the analogue of for the problem (1.7).
For the case of (1.13), the operator essentially obeys a chain rule when applied to the nonlinearity, which greatly aids the analysis. This fact relies crucially on the gauge-invariance of the nonlinearity, i.e. the symmetry . In the case of (1.7), the nonlinearity is not gauge-invariant. Because of this, estimating in becomes significantly more challenging. In fact, establishing a suitable estimate for this term is the principal difficulty in this paper.
Broadly speaking, we follow an approach known as the method of space-time resonances (cf. [6, 12]), which we now briefly describe in our setting. Introducing the interaction variable , we have
[TABLE]
To estimate the latter term, we first use the Duhamel formula for (cf. (4.1) below) to write an integral formula for . Because it is the quadratic terms that will ultimately be the most difficult to estimate, let us describe the technique for a single quadratic nonlinear term, say . For such a term, we would be led to estimate in the following term:
[TABLE]
where the phase . If the derivative lands on , then after applying Plancherel, we recover a copy of and have a term that is amenable to Strichartz estimates. We will see that such terms are not too difficult to estimate (cf. Section 5).
However, if the derivative lands on , then we are faced with estimating the following ‘phase derivative’ term:
[TABLE]
Such terms are significantly more difficult to estimate. In particular, we need to exhibit additional decay to overcome the factor . The idea is to exploit oscillation in the phase, which as usual is achieved via integration by parts.
On ‘time non-resonant’ regions (i.e. regions on which ), one can use the identity
[TABLE]
to integrate by parts with respect to , which either cancels the in (1.14) or leads to additional copies of the solution (which decays, due to the bootstrap assumption) via (1.7). On ‘space non-resonant’ regions (i.e. regions on which ), one can use the identity
[TABLE]
to integrate by parts with respect to , which cancels the in (1.14); if the derivative then lands on a copy of , we recover another copy of .
The strategy is then to decompose frequency space into regions of non-resonance and estimate the contribution of each region separately. Note that each type of quadratic term (namely, , , and ) will lead to a different phase in (1.14), and hence requires its own decomposition. For us, the most difficult case will come from the nonlinearity. In fact, in this case we will need to rely on one other notion of non-resonance, which we have called ‘angular non-resonance’; this refers to the non-vanishing of
[TABLE]
and uses a different integration by parts identity (see Section 6.3). It is only in treating the angular non-resonant region that we will need control over the angular momentum type quantity .
After integrating by parts on a non-resonant region, we can again write the resulting expressions back in terms of and . An example of a resulting term for the nonlinearity (in the time non-resonant case, say) is
[TABLE]
where is the bilinear multiplier with symbol with denoting a cutoff to the time non-resonant region. In general, the terms we need to estimate will now involve a bilinear multiplier applied to , , , or ; the particular multipliers that appear depend on whether the region is space, time, or angular non-resonant (see Sections 6.1–6.3). The multipliers that appear will not typically be amenable to standard bilinear estimates such as the Coifman–Meyer theorem. Indeed, even many of the cutoff functions used to partition frequency space into appropriate non-resonant regions already suffer from this problem. As in [12], it is essential to curtail the associated losses by exploiting the fact that the final result is to be estimated in . In this paper, this key bilinear estimate appears as Proposition 2.9, which extends an analogous result from [12].
There is one additional twist to what we described above, namely, the fact that a direct implementation of the space-time resonance approach with the quadratic nonlinearity
[TABLE]
in (1.7) does not work. Problems arise in quadratic frequency interactions with a small output frequency. For the first two terms above, the factor of provides cancellation at zero output frequency; however, the third term is missing this factor. For this reason, we employ a normal form transformation to remove the quadratic term and carry out the space-time resonance method with the normal form of the equation. See Section 4 for further discussion.
Having completed the outline of the proof, it now seems pertinent to discuss more fully its relation to the works [9, 12] on (1.10). First and foremost, our improvement on [12] is based on the discovery and subsequent exploitation of additional non-resonance phenomena in the case of high-high frequency interactions. Most important of these is the angular non-resonance discussed in Section 9.4; however, the use of time non-resonance in Section 9.3 is also novel and essential. On top of this, we believe that our paper introduces a number of individually minor, but collectively significant, improvements to the analysis in [12]. Two examples of this are the introduction and systematic use of the bound
[TABLE]
(see Lemma 3.2) and the extension of their bilinear estimate to include the end-point case (see Proposition 2.9). One end-result of these improvements is that we are able to present our result in a paper of comparable length while including a much more thorough exposition of the details. In truth, the brevity of [12] rendered us unable to reconstruct their arguments in several places.
There is less of a connection between what we do here and the paper [9], since the key point of that paper is to exploit additional Strichartz estimates that become available for radial data or data with higher angular regularity. The failure of such estimates in general manifests, for example, in the fact that while the traditional NLS is ill-posed in spaces of negative regularity, well-posedness can be restored by passing to radial data, or data with additional angular regularity.
1.2. Organization of the paper
In Section 2, we set up notation and collect some useful lemmas, including the bilinear estimate Proposition 2.9. In Section 3, we collect some consequences of the boundedness of the energy norm, including some decay estimates and control over the Strichartz norm. In Section 4, we discuss a normal form transformation used to ameliorate the effect of the quadratic terms in the nonlinearity. In Section 5, we begin to estimate the weighted norm, dealing with the cubic and higher terms, along with the quadratic terms without phase derivatives. Sections 6–10 are dedicated to estimating the quadratic terms containing the phase derivatives (see Proposition 6.1). These sections comprise the heart of the paper. Finally, in Section 11, we collect all of the estimates and complete the proof of Theorem 1.1.
Acknowledgements
We are grateful to an anonymous referee for their suggestion of the parameterization of the nonlinearity appearing in (1.3); the original version of this paper employed a different parameterization, which resulted in us treating only what are now the cases .
This work was partially supported by a grant from the Simons Foundation (#342360 to Rowan Killip). R. K. was further supported by NSF grants DMS-1265868 and DMS-1600942. J. M. was supported by the NSF Postdoctoral Fellowship DMS-1400706 at the University of California, Berkeley. M. V. was supported by NSF grant DMS-1500707. Part of the work on this project was supported by the NSF grant DMS-1440140, while the authors were in residence at the Mathematical Sciences Research Institute in Berkeley, California, during the Fall 2015 semester.
2. Notation and useful lemmas
For non-negative quantities and , we write to denote for some . We write to denote for some small . We write if and . We write and . We write to denote for some small .
We write a complex-valued function as . We write to indicate that either or may appear. When is a monomial, we write to denote a finite linear combination of the factors of , where Mikhlin multipliers (for example, Littlewood–Paley projections or the operator defined below) and/or complex conjugation may be additionally applied in each factor. We extend Ø to polynomials via .
In what follows, denotes derivatives with respect to with fixed and . On the other hand, indicates derivatives with respect to with fixed and . Similarly, indicates derivatives with respect to with fixed and .
For a time interval , we write for the Banach space of functions equipped with the norm
[TABLE]
with the usual adjustments when . If , we write . We will often abbreviate
[TABLE]
We write for the Hölder dual of , i.e. the solution to .
At times, we will make use of the Lorentz spaces defined via the quasi-norm
[TABLE]
where and . We have and for . We also have the following generalized Hölder and Hardy–Littlewood–Sobolev inequalities:
Lemma 2.1** (Hölder and Hardy–Littlewood–Sobolev in Lorentz spaces, [13, 19]).**
*
The following estimates hold:*
- (i)
[TABLE]
whenever and satisfy and .
- (ii)
[TABLE]
whenever and satisfy and .
Lemma 2.1(ii) implies the following Sobolev embedding on :
[TABLE]
We define the Fourier transform on via
[TABLE]
We define via . We make use of the following Fourier multiplier operators:
[TABLE]
We also use the notation .
We employ the standard Littlewood–Paley theory. Let be a radial bump function supported in and equal to one on the unit ball. Let . For , we define the Littlewood–Paley projections
[TABLE]
We also write , and similarly for the other operators. These operators commute with all other Fourier multiplier operators. They are self-adjoint and bounded on every and space for and . They obey the following standard estimates.
Lemma 2.2** (Bernstein).**
Let and . Then
[TABLE]
The operators are not true projections, in the sense that . As a substitute, we introduce the ‘fattened’ Littlewood–Paley multiplier and define the operators analogously to the above. We then have .
2.1. Linear estimates
In this section, we record some estimates for the linear propagator .
Proposition 2.3** (Dispersive estimates).**
For and ,
[TABLE]
For and ,
[TABLE]
The proof of this proposition is a standard application of stationary phase methods. The key information is the eigenvalues of the Hessian of ; these are
[TABLE]
and the latter has multiplicity .
From these dispersive estimates, one can deduce the following Strichartz estimates (see [15], for example). In the following, we call an admissible pair if and . We call a dual admissible pair if is an admissible pair.
Proposition 2.4** (Strichartz estimates).**
Let be a time interval and . Let be an admissible pair and let be a dual admissible pair. Then
[TABLE]
2.2. Bilinear estimates
We will frequently encounter bilinear Fourier multiplier operators of the following type.
Definition 2.5** (Bilinear operators).**
Given a function , we define
[TABLE]
For functions we write for the multiplier with symbol given by the pointwise product .
Remark 2.6**.**
The symbol corresponds to the pointwise product. The condition characterizes those multipliers that map real-valued functions to real-valued functions.
We will often work with bilinear multipliers on the Fourier side. Note that
[TABLE]
We will consistently use the notation
[TABLE]
and write
[TABLE]
We will need several estimates for bilinear operators. We begin with some standard results concerning sufficiently smooth multipliers.
Definition 2.7** (Coifman–Meyer–Mikhlin).**
We call a Fourier multiplier a Mikhlin multiplier if
[TABLE]
for all multiindices up to sufficiently high order.
We call a bilinear multiplier a Coifman–Meyer multiplier if
[TABLE]
for all multiindices up to sufficiently high order.
We call a bilinear multiplier a Coifman–Meyer–Mikhlin multiplier if it may be written in the form
[TABLE]
where are Mikhlin and is Coifman–Meyer.
Combining the Coifman–Meyer theorem [3, 4] together with the standard Mikhlin multiplier theorem, we obtain the following result:
Proposition 2.8** (Coifman–Meyer–Mikhlin estimate [3, 4]).**
If is Coifman–Meyer–Mikhlin, then
[TABLE]
for all and satisfying .
The next bilinear estimates allow for more singular multipliers. These estimates are a slight strengthening of similar estimates appearing in [12] (cf. [12, Lemma 9.1, Corollary 9.3]). Pointwise symbol estimates are replaced by conditions involving -based Sobolev regularity. We need both a time-independent bilinear estimate and a Strichartz-type bilinear estimate. Given a symbol , we will consider the following ‘norm’:
[TABLE]
where for fixed we consider as a function of via or as a function of via .
Proposition 2.9** (Bilinear estimates).**
Let be the symbols of bilinear operators .
- (i)
For any ,
[TABLE]
- (ii)
For any dual admissible pair and any satisfying and ,
[TABLE]
Remark 2.10**.**
No assumptions are needed for the operator beyond boundedness. In Lemma 4.1(ii) below we will see an example of a bilinear multiplier that is readily seen to be bounded (by relatively elementary means) but which is not of Coifman–Meyer–Mikhlin type. A larger class of symbols obeying bounds is admitted in [12]; this can be immediately recovered from the above by dyadic splitting and the triangle inequality. We prefer to formulate our result in terms of the norm (2.4) since in practical instances, this seems easiest to bound in an effective manner.
Proof.
By symmetry, it suffices to prove the result with the first factor on the right-hand side of (2.4). Define so that
[TABLE]
As
[TABLE]
uniformly for (by Plancherel), we may choose so that
[TABLE]
uniformly for .
(i) Fix and . By Plancherel,
[TABLE]
Using Plancherel again and the fact that translation is an isometry on , we deduce
[TABLE]
By Cauchy–Schwarz and (2.5),
[TABLE]
The estimate in (i) follows.
(ii) In this case, we write
[TABLE]
Thus, applying Strichartz (Proposition 2.4) and estimating as above,
[TABLE]
The estimate in (ii) follows. ∎
3. Decay via the energy norm
In this section we prove that control over the energy norm
[TABLE]
where implies decay for the solution , as well as for . In fact, one does not need to use the term to prove decay. The estimates we prove are in the spirit of Klainerman–Sobolev inequalities. As a consequence of the decay estimates, we can then prove an a priori estimate for the Strichartz-norm of . At the end of this section, we also record some bounds for norms with weights that will play a role in controlling the -norm of .
We begin with an elementary lemma.
Lemma 3.1** (Control of low frequencies).**
[TABLE]
Proof.
As the -norm controls the -norm, it suffices to consider the low frequencies. For this, we use the Sobolev embedding (2.1) and Hölder’s inequality (cf. Lemma 2.1) to estimate
[TABLE]
The result follows.∎
Next, we prove decay.
Lemma 3.2** (Decay).**
Fix .
- (i)
**
- (ii)
For ,
[TABLE]
Consequently,
[TABLE]
- (iii)
For any ,
[TABLE]
Proof.
First note that by Sobolev embedding,
[TABLE]
Thus, to prove (i) and (ii), it suffices to show that control over implies decay for times .
To this end, we note that the dispersive estimate (Proposition 2.3) and Hölder’s inequality imply
[TABLE]
giving (i).
Using instead the frequency-localized dispersive estimate (2.2), we easily obtain (3.2). Thus, for ,
[TABLE]
which together with (i) settles (ii).
The estimates in (iii) follow by interpolating between (i) and (ii).∎
We will next prove control over the Strichartz norm of . As the nonlinearity is described in terms of , it is useful to record the following lemma.
Lemma 3.3** (Control of ).**
Let and let be a time interval. Then with all space-time norms over ,
[TABLE]
Proof.
As is bounded on and and commutes with any Fourier multiplier operator whose symbol is radial, it is clear that
[TABLE]
For the low frequencies, we use Bernstein and (3.1) to see that
[TABLE]
Next, we use Lemma 3.2 to estimate
[TABLE]
Finally, using , unitarity of , and boundedness of Riesz potentials,
[TABLE]
The result follows. ∎
We can now prove control of the Strichartz norm (cf. (1.11)).
Proposition 3.4**.**
Let and suppose solves (1.7). Write . There exists such that
[TABLE]
Proof.
Recalling that commutes with , the estimate (3.4) follows immediately from Proposition 2.4.
We turn to (3.5). An application of Proposition 2.4 shows that to prove (3.5), we need to estimate
[TABLE]
To do this, we will rely on the decay for provided by Lemma 3.2, as well as Lemma 3.3. We also note that obeys Leibniz and chain rules.
For the quadratic and cubic terms, we may estimate
[TABLE]
for .
We turn to the quartic and quintic terms. Writing and distributing derivatives, we see that to estimate the terms involving , it suffices to consider the contribution of the following terms:
[TABLE]
Recalling boundedness of (and the fact that this operator commutes with ), we can first bound
[TABLE]
Finally, we note that by Bernstein and Lemma 3.3 (and the fact that commutes with frequency projections), we have
[TABLE]
Thus, estimating with the same spaces as above, we have
[TABLE]
The result follows. ∎
Remark 3.5**.**
Note that the proof above gives
[TABLE]
which we will use in Section 10.
The bulk of the remainder of the paper is devoted to controlling the weighted norm. Taking the Fourier transform, we compute
[TABLE]
To estimate the contribution of the cubic and higher order terms in the nonlinearity to the weighted norm, we will write as in (3.7) and estimate each resulting piece separately. This requires bounds on , for which we rely on the following lemma.
Lemma 3.6** (Bounds for weights).**
Let and . Then
[TABLE]
Proof.
As the commutator is bounded on and for and , it suffices to prove the estimates (3.8) and (3.9) for .
To this end, we write
[TABLE]
As the commutator , we obtain
[TABLE]
Next, using Lemmas 3.1 and 3.2 together with Sobolev embedding, we estimate
[TABLE]
The result follows. ∎
4. Normal form transformation
In this section, we discuss normal form transformations for (1.7). The use of normal forms in PDE originated in the work of Shatah [23] and has since become a widely-used technique in the setting of nonlinear dispersive equations. Briefly, the idea is to look for a change of variables with the effect of removing the most problematic terms from the nonlinearity.
In our setting it is the quadratic terms, namely,
[TABLE]
that are the most difficult to estimate. Recall also that and . In our previous work concerning the final-state problem [16], we used a normal form transformation to eliminate the term from the nonlinearity. In [16], we worked with energy-space solutions (i.e. no weighted assumption) and the poor spatial decay of was the biggest difficulty in handling the quadratic terms; in particular, we had only , compared to . Due to the nature of our arguments (which entailed testing against a dense subclass of test functions, namely, those supported away from the origin in frequency space), the interaction output frequency did not play an important role. In particular, the presence of in the real part of the nonlinearity offered no real advantage.
In this work, to treat the quadratic terms will require a careful decomposition of frequency space into ‘non-resonant’ regions, which requires an accounting of all of the various quadratic frequency interactions; see Sections 6–9. In this analysis, the operator will be crucial in our treatment of interactions with small output frequency. Accordingly, the absence of in the imaginary part of the nonlinearity is problematic. Note also that the weighted assumption gives better control over (cf. Section 3), which already makes the term less intimidating in this setting. Thus, in this work, we will employ a normal form transformation to remove the term . The authors of [12] also elected to remove this term in order to treat the initial-value problem for Gross–Pitaevskii, for the same reasons. As we will see, the normal form transformation that removes the term also cancels the two copies of appearing in the term, which actually puts the and terms on essentially equal footing. It is worth noting that the symmetry conditions required of the normal form transformation make it impossible to remove all of the quadratic terms.
The standard approach to normal form transformations is to introduce a new variable of the form
[TABLE]
where is a bilinear operator that is chosen so that
[TABLE]
cancels the problematic term(s) in the nonlinearity. Note that this will add some new terms to the nonlinearity, as well. One then proves estimates for , which (hopefully) solves a better equation than , and then inverts the normal form to deduce information about .
We take the following closely-related approach. First, we rewrite (1.7) in integral form:
[TABLE]
Given a bilinear operator as above, we can write
[TABLE]
As before, we wish to choose to cancel the problematic quadratic terms. In place of inverting the transformation, we will instead need to estimate the quadratic boundary term above.
We turn now to the details. We take
[TABLE]
where are symmetric, bilinear Fourier multiplier operators to be determined below. Using (1.7), we have
[TABLE]
In order to cancel the term , we need to choose and so that
[TABLE]
Imposing symmetry on and then leads to the unique choice
[TABLE]
where the symbol of is given by
[TABLE]
With this choice, we write
[TABLE]
where the are defined as follows:
[TABLE]
The structure of for will not be too important. The quadratic terms, however, have some nice structure that will be important in the analysis below:
[TABLE]
where
[TABLE]
One can check that
[TABLE]
In particular, is Coifman–Meyer, but is not. In fact, is not even Coifman–Meyer–Mikhlin; nonetheless, it is bounded, as the following lemma shows.
Lemma 4.1**.**
- (i)
The operator is Coifman–Meyer.
- (ii)
For any Coifman–Meyer–Mikhlin multiplier , the operator with symbol
[TABLE]
maps boundedly whenever and In particular, boundedly.
Proof.
Item (i) follows from (4.6). We focus on (ii) and argue by interpolation: Writing
[TABLE]
for , it suffices to show that the operator corresponding to maps boundedly. To this end, we note that and are Coifman–Meyer–Mikhlin with norms bounded by for some . Indeed,
[TABLE]
while for one exchanges and in the first three factors. Boundedness of can then be deduced by the three lines theorem (see, for example, [24, Lemma 4.2]). ∎
5. The weighted norm, part I
In this section, we begin estimating the -norm of , where is a solution to (1.7). To do this, we write the Duhamel formula for using the normal form (4.3).
We first consider the contribution of everything except the quadratic terms appearing in the nonlinearity.
Proposition 5.1**.**
Let and let if and otherwise. Let be a solution to (1.7). There exists such that for all ,
[TABLE]
Here is as in (4.2) and the are as in (4.4).
Proof.
To begin, recall that , so that
[TABLE]
This, together with the fact that is unitary on , immediately implies (5.1).
We turn to (5.2) and fix . By Plancherel, it is equivalent to estimate
[TABLE]
where and .
Distributing the derivative, we are led to estimate the following terms:
[TABLE]
where the multipliers and have symbols
[TABLE]
Note that while the operator with symbol is not Coifman–Meyer, it is Coifman–Meyer–Mikhlin. Indeed, this is apparent from
[TABLE]
Thus, using Proposition 2.8 and Lemma 3.2, we estimate
[TABLE]
Similarly,
[TABLE]
is Coifman–Meyer–Mikhlin. Thus, we can use Lemma 3.1 and Lemma 3.2 to estimate
[TABLE]
For the third term in (5.5), we split
[TABLE]
We estimate the contribution of the first summand by
[TABLE]
For the second summand, we note that is still a Coifman–Meyer multiplier. Thus
[TABLE]
Finally, we consider (5.3). For this, we again use (5.4). The nonlinearity consists of two types of terms, namely, the original terms coming from and the terms coming from the normal form.
For the original terms, we recall from (3.7) that we may write
[TABLE]
and we note that the commutator is given by Applying Proposition 2.4, we get
[TABLE]
where and denotes any dual Strichartz norm.
For the terms coming from the normal form, it suffices to estimate
[TABLE]
Once again, we use (5.6) and write
[TABLE]
where has the symbol . Note that and are Coifman–Meyer–Mikhlin operators.
Applying Strichartz, it remains to estimate the terms on the right-hand sides of (5.7) and (5.8) in any dual Strichartz space. Note that by (3.6), the contribution of with is acceptable.
We first consider the cubic contributions in the remaining terms. By Lemma 3.6 and Lemma 3.2, we may estimate as follows:
[TABLE]
Next, by Proposition 2.8 and Lemma 3.3,
[TABLE]
Let us consider now the quartic contributions. Arguing as above,
[TABLE]
and
[TABLE]
We now turn to the quintic contributions. Arguing as above,
[TABLE]
Next,
[TABLE]
Finally, we consider the sextic contribution, which comes only from the normal form. We rely on Sobolev embedding and the fact that is Coifman–Meyer. We get
[TABLE]
Using that and are Coifman–Meyer–Mikhlin, we estimate
[TABLE]
This completes the proof of Proposition 5.1. ∎
We turn to estimating the contribution of the quadratic terms arising from the nonlinearity in (4.3), that is,
[TABLE]
where we recall
[TABLE]
Thus, by Plancherel we are led to estimate terms of the form
[TABLE]
where , , and (cf. (4.5)). In the remainder of this section, we consider the easier terms in which the derivative misses the phase.
Proposition 5.2**.**
Let and let if and otherwise. Let be a solution to (1.7) and . For all ,
[TABLE]
where as in (4.5).
Proof.
If the derivative lands on , then we apply Proposition 2.4, Lemma 4.1, and Lemma 3.2 to estimate
[TABLE]
We next consider the operators
[TABLE]
Distributing the derivative, we first consider the term in which hits the factor in the multiplier . One can verify that this term is of the form
[TABLE]
where is Coifman–Meyer–Mikhlin. Thus, by Lemmas 4.1, 3.1, and 3.2,
[TABLE]
All remaining terms are either Coifman–Meyer–Mikhlin symbols or products of Coifman–Meyer–Mikhlin symbols with the multiplier
[TABLE]
Thus, by Proposition 2.8, Lemma 4.1, and Lemma 3.2, the contribution of these terms is bounded by
[TABLE]
This completes the proof of the proposition.∎
6. The weighted norm, part II
In this section we begin to estimate the contribution of the quadratic nonlinearity when the derivative hits the phase. The goal of the next several sections is to prove the following:
Proposition 6.1**.**
Let and let if and otherwise. There exists such that for all ,
[TABLE]
where as in (4.5) and .
We briefly recall the general strategy. We consider separately the phases arising from the , , and terms. For each phase, we decompose into regions on which we have ‘non-resonance’, which refers to the non-vanishing of or some particular derivative of . We carry out these decompositions in Section 7.1, Section 8, and Section 9. On each such region, we perform an integration by parts in the term
[TABLE]
using an identity that capitalizes on the particular from of non-resonance; see (6.2), (6.6), and (6.11) below. The terms arising from integration by parts using a given identity are all of a similar form; we estimate these in Sections 10.1–10.3, relying heavily on Proposition 2.9 and Lemma 3.2.
We will always perform our decompositions and integrations by parts at fixed frequencies . As are not true projections, we will also use the fattened Littlewood–Paley projections , which ensure that .
As noted in the introduction, many of the bilinear multipliers used to partition frequency space will not not obey uniform bounds; indeed in some cases they constitute the worst term.
6.1. Time non-resonance
Time non-resonance refers to the non-vanishing of . Suppose is a cutoff to a region on which . We wish to estimate the contribution of
[TABLE]
to (6.1). As on the support of , we may use the identity
[TABLE]
to integrate by parts in the term above. Defining the multiplier
[TABLE]
we see that it suffices to estimate the following terms:
[TABLE]
Note that because we are working on dyadic time intervals, we can estimate the two terms in (6.4) in the same way. Note also that to arrive at (6.5), we have used the fact that
[TABLE]
6.2. Space non-resonance
Space non-resonance refers to the non-vanishing of . Suppose is a cutoff to a region on which . We wish to estimate the contribution of
[TABLE]
to (6.1). As on the support of , we may integrate by parts using the identity
[TABLE]
Defining the multipliers
[TABLE]
we see that it suffices to estimate the following terms:
[TABLE]
6.3. Angular non-resonance
For the phase corresponding to the nonlinearity, namely,
[TABLE]
we will use a different identity to integrate by parts on one region of frequency space. We adopt the following notation for the projections of into the directions orthogonal to and :
[TABLE]
We will refer to the non-vanishing of as angular non-resonance. Supposing is a cutoff to a region on which , we wish to estimate the contribution of
[TABLE]
to (6.1). As , one has
[TABLE]
Thus, as on the support of , we may integrate by parts in the term above using the identity
[TABLE]
Defining
[TABLE]
we are first led to estimate the following term, which arises when the divergence misses both copies of :
[TABLE]
The terms in which the divergence hits a copy of take the form
[TABLE]
For these terms, we first rely on the vector identity
[TABLE]
to write
[TABLE]
Exploiting and the radiality of , this further reduces to
[TABLE]
We next write
[TABLE]
Thus, defining
[TABLE]
we are led to estimate the following terms:
[TABLE]
7. Non-resonant decompositions I: Preliminaries and terms
The next three sections are devoted to the decompositions of frequency space into non-resonant regions for the phases corresponding to each quadratic nonlinear term.
The following lemma illustrates how we estimate many of the multipliers appearing below.
Lemma 7.1** (Typical estimate).**
Consider a multiplier of the form
[TABLE]
Suppose , and that for fixed the multiplier vanishes except for those on a set of volume . Suppose further that
[TABLE]
Then
[TABLE]
where we recall the notation from (2.4).
Proof.
Note that for , we have
[TABLE]
By the Leibniz rule, (7.1), and Hölder’s inequality, we arrive at the estimate
[TABLE]
which gives the result. ∎
In general, we will use cutoffs to decompose frequency space into regions of time, space, or angular non-resonance, as quantified by the bounds we can prove for and its derivatives. We will also need to keep track of the bounds satisfied by the derivatives of the cutoff functions.
We will consistently use the following notation.
Definition 7.2**.**
For we let be the angle defined via . We define
[TABLE]
where . We also denote . Note that .
For reference, we collect some bounds on the derivatives of , all of which follow from explicit computation.
Lemma 7.3** (Derivative bounds for ).**
Write where
[TABLE]
Then
[TABLE]
Consequently,
[TABLE]
The next lemma will be of frequent use in the following sections.
Lemma 7.4**.**
For any ,
[TABLE]
By the triangle inequality, one can replace \bigl{|}|x|-|y|\bigr{|} by .
Proof.
Without loss of generality, assume . Write
[TABLE]
By the fundamental theorem of calculus and Lemma 7.3,
[TABLE]
As and
[TABLE]
the result follows. ∎
7.1. Non-resonant decomposition for terms
The ‘decomposition’ for terms is particularly simple, due to the fact that there is always time non-resonance.
Proposition 7.5** (Non-resonant decomposition for terms).**
Fix . Define
[TABLE]
Then
[TABLE]
Remark 7.6**.**
The multiplier does depend on and . This will be the case with all the multipliers in the following sections. We have elected to omit any explicit reference to this dependence in our notation to reduce clutter.
Proof.
Using Lemma 7.3, one can check that
[TABLE]
Thus,
[TABLE]
Using this together with (7.1) and (4.6), we deduce (7.2). ∎
8. Non-resonant decompositions II: terms
In this section we carry out the decomposition into non-resonant regions for the terms in (6.1). This corresponds to the phase . Recall that we always work at fixed frequencies . By symmetry, it suffices to consider the case
[TABLE]
Proposition 8.1** (Non-resonant decomposition for terms).**
Given , there exists a decomposition such that the following holds: the multipliers
[TABLE]
satisfy
[TABLE]
8.1. Region 1. Time non-resonance
Define
[TABLE]
This is a Coifman–Meyer multiplier that restricts to the region
[TABLE]
In particular,
[TABLE]
As in (6.3), we define the multiplier
[TABLE]
Lemma 8.2** (Region 1 bounds).**
The following bound holds:
[TABLE]
Proof.
Using (8.2) and Lemma 7.3, one can check
[TABLE]
Thus,
[TABLE]
Using this together with the cutoff bounds ((7.1), (8.3)) and (4.6), we deduce
[TABLE]
which implies (8.5). ∎
Remark 8.3**.**
On the complement of Region 1 (i.e. on the support of ), we have
[TABLE]
8.2. Region 2. Time non-resonance
Let be a partition of unity adapted to a maximal -separated set on . Define
[TABLE]
and let
[TABLE]
Note that
[TABLE]
We define
[TABLE]
and
[TABLE]
Lemma 8.4** (Region 2 bounds).**
The following bound holds:
[TABLE]
Proof.
First note that on the support of we have . Thus,
[TABLE]
This implies , which in turn implies
[TABLE]
Using this together with Lemma 7.3 (and recalling ), one also finds
[TABLE]
It follows that
[TABLE]
Using this together with the cutoff bounds ((7.1), (8.3), (8.7)) and (4.6), we deduce
[TABLE]
which gives (8.9).∎
Remark 8.5**.**
On the complement of Region 2 (i.e. on the support of ), we have
[TABLE]
8.3. Region 3. Time non-resonance
Let be a partition of unity adapted to a maximal -separated set on . Define
[TABLE]
We let
[TABLE]
One can check directly that
[TABLE]
We define
[TABLE]
and let
[TABLE]
Lemma 8.6** (Region 3 bounds).**
The following bound holds:
[TABLE]
Proof.
We begin by collecting some facts related to the support of .
First note that by the definition of , we have
[TABLE]
so that
[TABLE]
Thus, using the law of sines, (8.10), and (8.6), we deduce
[TABLE]
In particular, , so that
[TABLE]
Next, note that for fixed , the multiplier restricts to a set of volume
[TABLE]
We now claim the lower bound
[TABLE]
To justify this, we treat the two alternatives in (8.13) separately, beginning with the first. Writing
[TABLE]
a direct computation gives
[TABLE]
On the other hand, recalling (8.6), we have
[TABLE]
As direct computation also gives
[TABLE]
we can use (8.14) and continue from above to estimate
[TABLE]
Thus (8.17) holds when .
We now turn to the justification of (8.17) under the second alternative in (8.13). In fact, we will show something stronger, namely,
[TABLE]
Indeed, under this condition, . Thus
[TABLE]
This implies and therefore and so (8.18) follows.
We next claim the upper bound
[TABLE]
Indeed, this follows from Lemma 7.4 and (8.14):
[TABLE]
To proceed, we break into three cases.
Case 1. We first suppose
[TABLE]
For which the bounds established in (8.17) and (8.19) will be sufficient. Indeed from these and Lemma 7.3, we have
[TABLE]
Thus,
[TABLE]
Combining this with the cutoff bounds ((7.1), (8.3), (8.7), (8.11)), (4.6), and the volume bound (8.16), we deduce
[TABLE]
which is acceptable.
Case 2. Next, we suppose
[TABLE]
In this case, we will again use the bounds in and (8.19); however, we will need the following additional estimate:
[TABLE]
To prove this, we first note that by (8.15) and (8.21), we have
[TABLE]
Thus, we may use Lemma 7.4, (8.1), (8.6), and (8.21) to bound
[TABLE]
Noting that Lemma 7.3 also gives
[TABLE]
we use (8.17), (8.19), and (8.22) to deduce
[TABLE]
Combining this with the cutoff bounds ((7.1), (8.3), (8.7), (8.11)), (4.6) and the volume bound (8.16), and recalling that in this regime, we deduce
[TABLE]
which is acceptable.
Case 3. Finally, we suppose
[TABLE]
Noting that Lemma 7.3 implies
[TABLE]
we use (8.18) (which holds under the assumption (8.23)) and (8.19) to deduce:
[TABLE]
Combining this with the cutoff bounds (7.1) ((8.3), (8.7), (8.11)), (4.6), and the volume bound (8.16), and recalling that in this regime, we deduce
[TABLE]
which is acceptable.
This completes the proof of Lemma 8.6. ∎
Remark 8.7**.**
On the complement of Region 3 (i.e. on the support of ), we have
[TABLE]
Thus
[TABLE]
By the law of sines and (8.6), the same is true of . In particular,
[TABLE]
8.4. Region 4. Space non-resonance
We finally define
[TABLE]
and let
[TABLE]
We first claim that
[TABLE]
We begin by using Lemma 7.4, (8.10), and the law of sines to estimate
[TABLE]
Next, using (8.26),
[TABLE]
Using (8.29) together with the first estimate in (8.28), we deduce the first bound in (8.27). On the other hand,
[TABLE]
This bound, together with the second estimate in (8.28), gives the second bound in (8.27).
Note that using (8.29) together with (8.6) and the law of sines, we also get
[TABLE]
We are now in a position to prove bounds for the multipliers and .
Lemma 8.8** (Region 4 bounds).**
The following bounds hold:
[TABLE]
Proof.
To begin, recall that For the terms in which the divergence misses the factor , we claim that we get an upper bound of , which (using (8.32)) one can check is acceptable. Indeed, using the cutoff bounds (7.1), (8.3), (8.7), together with (4.6), we see that if the divergence hits the product
[TABLE]
then we will have an additional factor of , but otherwise we can argue exactly the same as for . If the divergence hits , then we get an additional factor of (cf. (8.11)). However, in this case, we can also use the better volume bound (8.16). Thus, once again we face an additional factor of , which is acceptable.
Thus, to estimate , it suffices to treat and the multiplier
[TABLE]
We will first prove the general bounds for and , and then give the improvements when and .
Using Lemma 7.3, we find
[TABLE]
Using (8.27) and (8.31) as well (along with (8.25)), we find:
[TABLE]
We now claim that
[TABLE]
which implies the first bound in (8.32).
Proof of (8.35).
If no derivatives hit the product
[TABLE]
then we use the cutoff bounds ((7.1), (8.3), (8.7)), together with (8.34) (with ), (4.6), and the upper bound to estimate the contribution to the -norms by
[TABLE]
for , which is acceptable.
If any derivative hits , then we are in a position to use the volume bound (8.16). In particular, in the case that all derivatives land on , we may use (8.11), (8.16), and (8.34) (with ) to estimate the contribution to the -norms by
[TABLE]
for , which is acceptable. By using the second bound in (8.34) (with ) and the volume bound (8.16), we get the same estimate for when one derivative lands on each factor in (8.36), which is again an acceptable contribution.
It remains to consider the case when all derivatives land on for . For this, we use the first bound in (8.34), and for fixed we use spherical coordinates in the variable (with as the north pole) to compute the -norm. In particular, recalling (8.24), we may estimate the contribution to the -norms by
[TABLE]
which is acceptable. This completes the proof of (8.35). (Note that we do not need the case here; however, we will use it below.) ∎
We next claim
[TABLE]
which implies the first bound in (8.33).
Proof of (8.38).
We argue similarly to the case of (8.35). If no derivatives hit
[TABLE]
then we use (7.1), (8.3), (8.7), (4.6), the second bound in (8.34) (with ), and the upper bound to estimate the contribution to the -norms by
[TABLE]
for , which is acceptable.
If any derivative hits , then we are in a position to use the volume bound (8.16). In particular, in the case that all derivatives land on , we may use (8.11), (8.16), and the second bound in (8.34) (with ) to estimate the contribution to the -norms by
[TABLE]
for , which is also acceptable. By using the second bound in (8.34) (with ) and the volume bound (8.16), we get the same estimate for when one derivative lands on each factor in (8.39), which is again an acceptable contribution.
It remains to consider the case when all derivatives land on for . For this, we estimate as before; that is we use the first bound in (8.34), and for fixed we use spherical coordinates in the variable (with as the north pole) to compute the -norm. By (8.37), we estimate the contribution to the -norms by
[TABLE]
for , which is acceptable.∎
So far, we have established the first estimates in (8.32) and (8.33). To complete the proof of Lemma 8.8, we need to consider the case
[TABLE]
Using Lemma 7.3 and (8.40), we find
[TABLE]
Thus, recalling (8.27) and (8.30), we find:
[TABLE]
Using this together with (7.1), (8.3), (8.7), (8.11), and (4.6), it is not hard to verify that
[TABLE]
which gives the second bound in (8.32). Similarly,
[TABLE]
giving the second bound in (8.33). This completes the proof of Lemma 8.8. ∎
Proof of Proposition 8.1.
The proof of Proposition 8.1 follows immediately from Lemmas 8.2, 8.4, 8.6, and 8.8. One may also check that the cutoffs sum to . Note that we always choose the worst multiplier bound in the statement of Proposition 8.1. For example, the bounds for and when are better than those stated; however, the stated bound is attained by . ∎
9. Non-resonant decompositions III: terms
In this section we carry out the decomposition into non-resonant regions for the terms in (6.1). By symmetry, it suffices to consider the phase
[TABLE]
As before, we will work at fixed frequencies . To simplify notation, we will suppress the subscripts on our multipliers.
The result of this section is the following.
Proposition 9.1** (Non-resonant decomposition for terms).**
Given and , there exists a decomposition on the support of such that the following holds:
- •
For , define
[TABLE]
Then
[TABLE]
- •
For , define
[TABLE]
Then
[TABLE]
and
[TABLE]
with having the same support properties as .
- •
Finally, define
[TABLE]
Then
[TABLE]
with in the support of these multipliers.
To begin, we record a lemma that allows us to exploit some cancellation in derivatives of .
Lemma 9.2**.**
Let . For ,
[TABLE]
Proof.
Using and , we deduce
[TABLE]
Thus, the result is a consequence of the fundamental theorem of calculus and Lemma 7.3.∎
9.1. Region 1. Time non-resonance
Define
[TABLE]
and vanishing otherwise. We then define
[TABLE]
Where this is non-zero, we have
[TABLE]
Lemma 9.3** (Region 1 bounds).**
The following bound holds:
[TABLE]
Moreover, in the support of .
Proof.
Using (9.1) and Lemma 7.3, we deduce for :
[TABLE]
Thus,
[TABLE]
Recalling also (7.1), we deduce
[TABLE]
which gives the lemma. ∎
Remark 9.4**.**
In the remaining cutoffs, a factor will enforce the constraint
[TABLE]
9.2. Region 2. Time non-resonance
Let be a partition of unity adapted to a maximal -separated set on . Define
[TABLE]
and let
[TABLE]
One can check that
[TABLE]
We define
[TABLE]
and
[TABLE]
Lemma 9.5** (Region 2 bounds).**
The following bound holds:
[TABLE]
Moreover, in the support of .
Proof.
On the support of we have . Thus and so
[TABLE]
In particular, using (9.2), we deduce that
[TABLE]
which justifies the last assertion in the lemma, as well as
[TABLE]
These relations in turn yield
[TABLE]
We note that by Lemma 9.2, we have
[TABLE]
while by Lemma 7.3,
[TABLE]
Using this together with (9.6), we find:
[TABLE]
Recalling the cutoff bounds (7.1) and (9.3) (and relying on the factor ), we deduce
[TABLE]
which completes the proof of Lemma 9.5. ∎
Remark 9.6**.**
On the complement of Region 2, we have
[TABLE]
9.3. Region 3. Time non-resonance
Let and be as in Region 2. Define
[TABLE]
and let
[TABLE]
One verifies that
[TABLE]
We define
[TABLE]
and vanishing otherwise (see also (9.2)). Note that will be chosen below.
We define
[TABLE]
Lemma 9.7** (Region 3 bounds).**
The following bound holds,
[TABLE]
Moreover, on the support of .
Proof.
From (9.2) and definition of , we have that
[TABLE]
in the support of . This justifies the last assertion in the lemma.
Turning to bounding , we first claim that
[TABLE]
To see this, we introduce the function
[TABLE]
A direct computation gives
[TABLE]
Note that
[TABLE]
so that
[TABLE]
Noting also that in the support of , we have
[TABLE]
provided we choose the constant appropriately in (9.11).
We next claim
[TABLE]
[TABLE]
which gives the case of (9.14). The cases follow from Lemma 7.3, recalling that (9.12) holds.
Finally, using Lemma 9.2 and (9.12), we have
[TABLE]
Using (9.13), (9.14), and (9.15), we deduce:
[TABLE]
Using this together with the cutoff bounds ((7.1), (9.3), (9.10)), (4.6), and relying on the factor , we deduce
[TABLE]
which completes the proof of the Lemma 9.7. ∎
Remark 9.8**.**
On the support of , we have
[TABLE]
9.4. Region 4. Angular non-resonance
Define
[TABLE]
and vanishing otherwise (see also (9.2)). Note that
[TABLE]
On the support of , we will use the integration by parts described in Section 6.3.
Define
[TABLE]
where we recall the notation .
Recalling the discussion in Section 6.3, we need to prove multiplier bounds for
[TABLE]
Lemma 9.9** (Region 4 bounds).**
The following bounds hold:
[TABLE]
Moreover, in the support of the multipliers above.
Proof.
First note that in the support of , we have
[TABLE]
Using (9.16) and (9.19) one readily checks
[TABLE]
For higher derivatives, one can check that
[TABLE]
From Lemma 7.3, we also have
[TABLE]
Using the last three bounds, one deduces
[TABLE]
Next, note that
[TABLE]
Using this together with the cutoff bounds ((7.1), (9.3), (9.10), (9.17)), (4.6), and (9.21) (and recalling the presence of ), we deduce
[TABLE]
giving the first estimate in (9.18).
Next, we recall (9.19) and note that
[TABLE]
Combining this with the bounds above (i.e. all of the cutoff bounds and the estimate (9.21)), one can deduce
[TABLE]
giving the second estimate in (9.18).
Finally, we compute directly
[TABLE]
Thus one can check
[TABLE]
and we can again estimate as above to deduce
[TABLE]
giving the final estimate in (9.18). Note that in contrast to the first two estimates in (9.18), we use the fact that . ∎
9.5. Region 5. Space non-resonance
Let
[TABLE]
and vanishing otherwise. One can readily check that
[TABLE]
and also that on the support of ,
[TABLE]
We define the multipliers
[TABLE]
Lemma 9.10** (Region 5 bounds).**
The following bounds hold:
[TABLE]
with and on the support of these multipliers.
Proof.
By Lemma 7.3, (9.23), and the fundamental theorem of calculus,
[TABLE]
Similarly,
[TABLE]
For higher derivatives, we use Lemma 7.3 and (9.23) to estimate
[TABLE]
Using this together with (9.24) and (9.25), we deduce
[TABLE]
This, together with the cutoff bounds ((7.1), (9.3), (9.10), (9.17), and (9.22)) and (4.6), implies
[TABLE]
giving the first bound. Distributing the derivative in and using the same bounds above yields
[TABLE]
giving the second bound. This completes the proof of Lemma 9.10. ∎
Remark 9.11**.**
Combining the cutoff functions defined thus far, we have
[TABLE]
9.6. Region 6. Time non-resonance
Let be a partition of unity adapted to a maximal -separated set on . Define
[TABLE]
and let
[TABLE]
Note that
[TABLE]
We define
[TABLE]
and vanishing otherwise. Note that
[TABLE]
We define the multiplier
[TABLE]
Lemma 9.12** (Region 6 bounds).**
The following bound holds:
[TABLE]
Moreover, on the support of .
Proof.
First note that for fixed , the multiplier restricts to a set of volume
[TABLE]
By construction, on the support of ; thus
[TABLE]
We first claim the lower bound
[TABLE]
We begin by writing
[TABLE]
By direct computation, one verifies
[TABLE]
For the remaining piece, we use the fundamental theorem of calculus, Lemma 7.3, and (9.29) to get
[TABLE]
Thus (9.30) holds.
Next, we claim
[TABLE]
Indeed, using Lemma 7.4, we first get
[TABLE]
By (9.8), the law of sines, and (9.29),
[TABLE]
Thus (recalling ), (9.31) holds.
For higher derivatives, Lemma 7.3 gives
[TABLE]
Next, we will show
[TABLE]
In fact, for this follows directly from Lemma 9.2. For , we first have by Lemma 7.4,
[TABLE]
As , the law of sines and (9.29) yield
[TABLE]
as needed.
Using (9.30), (9.31), (9.32), and (9.33), we deduce
[TABLE]
Combining this with the cutoff bounds ((7.1), (9.3), (9.10), (9.17), (9.22), and (9.27)), (4.6), and the volume bound (9.28), we deduce
[TABLE]
giving the desired bound. This completes the proof of Lemma 9.12. ∎
Remark 9.13**.**
On the support of , one has
[TABLE]
so that
[TABLE]
Before proceeding to Regions 7 and 8, let us pause to derive a lower bound that we will be using in both: Assume first that . Then
[TABLE]
As \bigl{(}\tfrac{\xi_{1}}{|\xi_{1}|}+\tfrac{\xi_{2}}{|\xi_{2}|}\bigr{)}\cdot\tfrac{\xi_{1}}{|\xi_{1}|}=1+\cos(\theta_{12})\geq 0, it follows that
[TABLE]
Repeating the preceding argument with roles reversed, shows that (9.36) actually holds whenever .
9.7. Region 7. Space non-resonance
Let be a partition of unity adapted to a maximal -separated set . Define
[TABLE]
and set
[TABLE]
We now define
[TABLE]
and vanishing otherwise. Note that
[TABLE]
We also define
[TABLE]
Lemma 9.14** (Region 7 bounds).**
The following bounds hold:
[TABLE]
Moreover, on the support of these multipliers.
Proof.
Note that the cutoff guarantees
[TABLE]
Furthermore, for fixed , this restricts to a set of volume
[TABLE]
We first establish
[TABLE]
Indeed, using Lemma 7.4 with (9.8), the law of sines, and (9.38), we find
[TABLE]
which gives the case . The other cases follow from Lemma 7.3, recalling .
We next claim
[TABLE]
Recall from (9.8) that ; furthermore, by the law of sines, and (9.38), we have
[TABLE]
We conclude . Recalling (9.38) and the fact that , we deduce that
[TABLE]
for some . Thus , and so (9.36) gives (9.41).
We next note that Lemma 7.3 immediately gives
[TABLE]
Using (9.40), (9.41), and (9.42), one verifies
[TABLE]
We are now in a position to estimate the multipliers. Comparing (9.43) and (4.6) with all of the cutoff bounds ((7.1), (9.3), (9.27), and (9.37)), and exploiting the volume bound (9.39), one finds that
[TABLE]
This gives the first multiplier bound.
We turn to estimating . If the divergence lands on any term other than
[TABLE]
we will face an additional , but otherwise we can estimate exactly as above. Thus, we get an acceptable contribution for these terms (cf. ), and we are left to estimate
[TABLE]
Recalling the cutoff bounds (9.27) and (9.37), we see that we will lose an additional when estimating . In particular, using (9.43), the cutoff bounds, and the volume bound (9.39), we deduce
[TABLE]
giving the second multiplier bound. This completes the proof of Lemma 9.14.∎
Remark 9.15**.**
On the support of , we have
[TABLE]
9.8. Region 8. Space non-resonance
We define
[TABLE]
and vanishing otherwise. Corresponding to this we define multipliers
[TABLE]
Lemma 9.16** (Region 8 bounds).**
The following bounds hold:
[TABLE]
with on the support of these multipliers.
Proof.
To begin, recall that . For the terms in which the divergence misses the factor , we claim that we get an upper bound of , which is an acceptable bound. Indeed, using the cutoff bounds (7.1), (9.3), as well as (4.6), we first see that if the divergence hits the product
[TABLE]
then we get an additional factor of , but we will otherwise be able to estimate exactly as we do for below. If the divergence hits the product
[TABLE]
then we will get an additional factor of (cf. (9.27) and (9.37)); however, in this case, we will always be able to use the improved volume bounds in (9.28) or (9.39), which gains an additional in the arguments presented below. In particular, we again face only an additional factor of .
Thus, to treat , it suffices to get suitable estimates for and the multiplier
[TABLE]
We begin by applying Lemma 7.4 and using (9.8) and the law of sines to deduce
[TABLE]
Case 1. Recalling that , we first consider the case
[TABLE]
Note that in this case, \bigl{|}|\xi_{1}|-|\xi_{2}|\bigr{|}\gtrsim N_{1}, so that (9.46) and (9.36) give
[TABLE]
Thus by (9.36), the law of sines, (9.47), and (9.44), we have
[TABLE]
For higher derivatives, we have by Lemma 7.3,
[TABLE]
Using these bounds, one verifies
[TABLE]
Combining the second estimate in (9.49) with the cutoff bounds ((7.1), (9.3), (9.27), (9.37)), (4.6), and noting that in this regime, we deduce
[TABLE]
which is acceptable. (The worst terms arise when all of the derivatives hit .)
We turn to estimating in Case 1 (cf. (9.45)). We wish to show
[TABLE]
For this, it will suffice to prove
[TABLE]
First, if no derivatives hit
[TABLE]
then we use the cutoff bounds ((7.1), and (9.3)), (4.6), and the second estimate in (9.49) (with ) to estimate the contribution to the -norms by
[TABLE]
for , which is acceptable.
Next, if any derivative hits , then we may use the volume bound in (9.28) or (9.39). In particular, when all derivatives land on these cutoffs, we may use (9.27), (9.37), and the second estimate in (9.49) (with ) to estimate the contribution to the -norms by
[TABLE]
for , which is acceptable.
It remains to consider the following two situations:
- (i)
, with one derivative hitting and the other hitting the cutoffs ,
- (ii)
and all derivatives land on .
For these terms, for fixed we use spherical coordinates in (with as the north pole) to compute the -norms. We also recall (9.44).
Using (9.27), (9.37), and the first bound in (9.49) (with ), we estimate the contribution of terms in (i) to the -norm by
[TABLE]
which is acceptable in light of (9.50).
Similarly, we estimate the contribution of the terms in (ii) to the -norms by
[TABLE]
which is acceptable. This completes the proof of (9.50) and hence of Lemma 9.16 in Case 1.
Case 2. Suppose now that
[TABLE]
By (9.36), the law of sines, and (9.44), we have
[TABLE]
Combining this with the upper bound (9.46), we find
[TABLE]
For higher derivatives of , we appeal to Lemma 7.3 to find
[TABLE]
For higher derivatives of , Lemma 9.2 and (9.51) give
[TABLE]
Combining (9.51), (9.53), (9.52), (9.54), (9.55) and recalling the bounds (9.44), we deduce
[TABLE]
We are now in a position to prove the multiplier bounds, beginning with . We use the second estimate in (9.56), the cutoff bounds ((7.1), (9.3), (9.27), (9.37)), (4.6), the fact that in this regime, and (9.51) to arrive at the estimate
[TABLE]
which is acceptable.
We turn to estimating in Case 2 (cf. (9.45)). In this case, we need to show
[TABLE]
for which it will suffice to prove
[TABLE]
First, if no derivatives hit
[TABLE]
then we use the cutoff bounds ((7.1), and (9.3)), (4.6), and the second estimate in (9.56) (with ) to estimate the contribution to the -norms by
[TABLE]
for , which is acceptable.
If any derivative hits , then we may use the volume bound in (9.28) or (9.39). In particular, when all derivatives land on these cutoffs, we may use (9.27), (9.37), and the second estimate in (9.56) (with ) to estimate the contribution to the -norms by
[TABLE]
for , which is acceptable.
It remains to consider the following two situations:
- (i)
, with one derivative hitting and the other hitting the cutoffs ,
- (ii)
and all derivatives land on .
For these terms, for fixed we again use spherical coordinates in (with as the north pole) to compute the -norms. Recall also (9.44).
Using (9.27), (9.37), and the first bound in (9.56) (with ), we estimate the contribution of terms in (i) to the -norm by
[TABLE]
which is acceptable.
Similarly, we estimate the contribution of the terms in (ii) to the -norms by
[TABLE]
which is acceptable. This completes the proof of Lemma 9.16 in Case 2. ∎
Proof of Proposition 9.1.
Proposition 9.1 follows from Lemmas 9.3, 9.5, 9.7, 9.9, 9.10, 9.12, 9.14, and 9.16. One may also check that the multipliers sum to on the support of . Note that in the statement of Proposition 9.1, we always choose the worst bounds; for example, the bounds obtained for in Lemma 9.16 are strictly worse than those for in Lemma 9.14. ∎
10. Estimation of non-resonant terms
In this section, we complete the proof of Proposition 6.1 by estimating the contribution of the time, space, and angular non-resonant regions of frequency space to (6.1).
Recall that in Sections 6.1–6.3, we collected the terms that we need to estimate in order to deal with each type of non-resonant region, namely (6.4) and (6.5) for time non-resonant regions, (6.9) and (6.10) for space non-resonant regions, and (6.14), (6.17)–(6.19) for angular non-resonant regions. These involve certain bilinear operators, namely, for time non-resonant terms (cf. (6.3)), and for space non-resonant terms (cf. (6.7), (6.8)), and , , for angular non-resonant terms (cf. (6.13), (6.15), (6.16)). In Sections 7–9, we decomposed frequency space into non-resonant regions for each type of quadratic nonlinearity and proved bounds for the resulting bilinear operators (cf. Propositions 7.5, 8.1, and 9.1).
In this section, we will show that the bounds established in Propositions 7.5, 8.1, and 9.1 suffice to prove the desired estimate in Proposition 6.1.
Our main tools will be the bilinear estimate Proposition 2.9 and the decay estimates in Lemma 3.2. For convenience, we record some particular consequences of Lemma 3.2 here. We let and be as in Proposition 6.1, that is, if and otherwise. Then for , we have
[TABLE]
Indeed, (10.1) follows from interpolation between and , Lemma 3.2, and the -control of . We apply (10.1) most often with , i.e.
[TABLE]
We also use the following consequence of Bernstein and Lemma 3.2:
[TABLE]
Remark 10.1**.**
As in Proposition 6.1 we work on dyadic intervals, when we apply Proposition 2.9 with a dual Strichartz norm , we always have the option to bound the contribution of all time integrals by .
10.1. Estimation of time non-resonant terms
To estimate the contribution of time non-resonant regions to (6.1), we need to bound terms of the form (6.4) and (6.5), which involve the multiplier (cf. (6.3)). To estimate these terms, we will use the bilinear estimate Proposition 2.9.
By virtue of symmetries inherent in the terms that we need to estimate and in light of the multiplier bounds in Propositions 7.5, 8.1, and 9.1, we see that we may reduce to the case ; moreover, we see that it is natural to cover this with three regimes: (i) , (ii) , and (iii) . Of course, the second regime is contained in the third, but we have to deal with a different multiplier bound in this particular case (cf. Regions 2 and 3 in Proposition 9.1). We record here the bounds we must contend with in each regime:
[TABLE]
1. We begin by estimating terms of the form (6.4). We can treat the cases and together, using the worse multiplier bound of in (10.2). Using the dual Strichartz norm with Proposition 2.9, we estimate these terms via Bernstein, Lemma 3.2, and (10.1):
[TABLE]
which is acceptable.
2. We estimate terms of the form (6.4) in the regime using Proposition 2.9 with the dual Strichartz norm and Lemma 3.2:
[TABLE]
which is acceptable.
3. We next consider terms of the form (6.5) in the regime . Note that since , it suffices to treat terms in which falls on the nonlinearity.
Using Proposition 2.9 with the dual Strichartz norm , (10.2), (10.1), and Lemma 3.2, we may bound the contribution of the quadratic part of by
[TABLE]
We estimate the contribution of cubic and higher terms via (10.2), Bernstein, Lemma 3.2, and (3.6): for ,
[TABLE]
which is acceptable.
4. We next estimate terms of the form (6.5) in the regime .
We first consider the case when falls on the nonlinearity. To estimate the contribution of the quadratic terms, we use Proposition 2.9 with the dual Strichartz norm , (10.2), Bernstein, and Lemma 3.2:
[TABLE]
which is acceptable. We next consider the contribution of the cubic and higher terms. We use Proposition 2.9 with the dual Strichartz norm , (10.2), Bernstein, Lemma 3.2, and (3.6): for , we bound these terms by
[TABLE]
which is acceptable.
We turn to the case when lands on the nonlinearity. For the quadratic terms, we use Proposition 2.9 with the dual Strichartz norm . We get the contribution
[TABLE]
which is acceptable. For the cubic and higher terms, we instead use and rely on (3.6) and Bernstein: for we get the contribution
[TABLE]
which is acceptable.
5. Next, we estimate terms of the form (6.5) in the regime . To begin, we consider only the contribution of the cubic and higher terms.
We first consider the case when falls on the nonlinearity. We use Proposition 2.9 with the dual Strichartz norm , (3.6), Bernstein, and Lemma 3.2 to obtain for the contribution
[TABLE]
which is acceptable.
Next, we consider the case when falls on the nonlinearity. We again use Proposition 2.9 and dual Strichartz norm , (3.6), Bernstein, and Lemma 3.2 to obtain for the contribution
[TABLE]
which is acceptable.
6. Finally, we treat the contribution of the quadratic terms in the nonlinearity to (6.5) when . We will use Proposition 2.9 with the dual Strichartz norm .
Recalling that and , we note that the quadratic terms in (1.7) are of the form
[TABLE]
If lands on the term, we estimate via Lemma 3.2:
[TABLE]
which is acceptable. If lands on the term, we get the contribution
[TABLE]
which is acceptable.
If lands on the term, then we estimate by Bernstein and Lemma 3.2:
[TABLE]
which is acceptable.
Finally, if lands on the term, we get the contribution
[TABLE]
which is acceptable. This completes the estimation of the contribution of the time non-resonant regions.
10.2. Estimation of space non-resonant terms
To estimate the contribution of the space non-resonant regions to (6.1), we need to bound terms of the form (6.9), involving the multiplier, and (6.10), involving the multiplier. Again, we will rely on Proposition 2.9.
Note that we only had space non-resonant regions for the and nonlinearities; the multiplier bounds we established appear in Propositions 8.1 and 9.1. In all cases of space non-resonance, we had . Examining the multiplier bounds, it is natural to split into two cases: and . We record here the worst bounds for each multiplier in these two regimes:
[TABLE]
We first consider terms of the form (6.9) in the case . We will use Proposition 2.9 with the dual Strichartz norm . We estimate such terms via (10.3), Bernstein, and Lemma 3.2:
[TABLE]
which is acceptable.
We next consider terms of the form (6.9) in the case . We will use Proposition 2.9 with the dual Strichartz norm . We estimate such terms via (10.3) and (10.1):
[TABLE]
which is acceptable.
We turn to terms of the form (6.10) in the case . We will use Proposition 2.9 with the dual Strichartz norm . We estimate such terms via (10.3), using (3.2) for and Lemma 3.2(iii) for :
[TABLE]
which is acceptable.
Finally, we consider terms of the form (6.10) in the case . We will use Proposition 2.9 with the dual Strichartz norm . We estimate such terms via (10.3), Lemma 3.2, and Bernstein:
[TABLE]
which is acceptable. This completes the estimation of the contribution of space non-resonant regions to (6.1).
10.3. Estimation of angular non-resonant terms
To estimate the contribution of angular non-resonant regions to (6.1), we need to bound the terms (6.14), (6.17), (6.18), and (6.19).
There is only one angular non-resonant region in our decompositions. The relevant multiplier bounds appear in Proposition 9.1; in particular, the bounds are all . Note also that in the angular non-resonant region, we have the condition .
To estimate the term (6.14), we use Proposition 2.9 with the dual Strichartz norm and (10.1):
[TABLE]
which is acceptable.
We next consider (6.17) and (6.18). As , it suffices to treat (6.17). Using Proposition 2.9 with the dual Strichartz norm and (10.1), we get the contribution
[TABLE]
which is acceptable.
Finally, we estimate (6.19). Arguing as above, we estimate this term by
[TABLE]
which is acceptable. This completes the estimation of the contribution of the angular non-resonant regions to (6.1), and hence the proof of Proposition 6.1.
11. Proof of the main result
In this section, we finally prove the main result, Theorem 1.1. As discussed in the Introduction, [16] already guarantees that we have a unique global solution to (1.7) for initial data as in Theorem 1.1. Thus, it remains to show that the solution scatters. We will prove scattering forward in time.
Collecting the results of Proposition 3.4, Proposition 5.1, Proposition 5.2, and Proposition 6.1, we arrive at the following a priori estimate:
[TABLE]
for some , where . Using this, one can show that if is sufficiently small, then the solution obeys
[TABLE]
To prove scattering, it suffices to show that is Cauchy in , is Cauchy in , and is Cauchy in as .
We first use the Duhamel formula (4.1), Proposition 3.4, and (11.1) to estimate
[TABLE]
for . Thus is Cauchy in as . Arguing similarly (using Proposition 3.4 and recalling (1.11)) also shows that is Cauchy in as .
To show that is Cauchy in , we use the normal form of the equation; that is, we use the Duhamel formula with (4.3). Then, for , we use Proposition 5.1, Proposition 5.2, Proposition 6.1, and (11.1) to estimate
[TABLE]
Thus is Cauchy in as .
Finally, note that the argument above gives a rate of convergence of to the scattering state. This completes the proof of Theorem 1.1.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] I. V. Barashenkov, A. D. Gocheva, V. G. Makhankov, and I. V. Puzynin, Stability of the soliton-like “bubbles”. Phys. D 34 (1989), no. 1–2, 240–254. MR 0982390
- 2[2] M. Christ and M. Weinstein, Dispersion of small amplitude solutions of the generalized Korteweg–de Vries equation . J. Funct. Anal. 100 (1991), no. 1, 87–109. MR 1124294
- 3[3] R. Coifman and Y. Meyer, Commutateurs d’intégrales singulières et opérateurs multilinéaires . Ann. Inst. Fourier (Grenoble) 28 (1978), no. 3, xi, 177–202. MR 0511821
- 4[4] R. Coifman and Y. Meyer, Ondelettes et opérateurs. III. Opérateurs multilinéaires. Actualités Mathématiques, Hermann, Paris, 1991. pp. i–xii and 383–538. MR 1160989
- 5[5] P. Gérard, The Cauchy problem for the Gross–Pitaevskii equation. Ann. Inst. H. Poincaré Anal. Non Linéaire 23 (2006), no. 5, 765–779. MR 2259616
- 6[6] P. Germain, N. Masmoudi, and J. Shatah, Global solutions for 3D quadratic Schrödinger equations. Int. Math. Res. Not. IMRN 2009, no. 3, 414–432. MR 2482120
- 7[7] V. L. Ginzburg and L. P. Pitaevskii, On the theory of superfluidity. Sov. Phys. JETP 34 (1958), no. 7, 858–861; translation from Zh. Eksp. Teor. Fiz., 32 (1958), 1240–1245. MR 0105929
- 8[8] V. L. Ginzburg and A. A. Sobyanin, Superfluidity of helium II near the λ 𝜆 \lambda -point. Sov. Phys. Usp. 19 (1976), no. 10: 773.
