This paper proves that solutions to a perturbed sine-Gordon equation stay close to a virtual solitary manifold constructed via iterative deformation, with stability and accuracy improving with the perturbation's differentiability.
Contribution
It introduces a method to construct a virtual solitary manifold for the perturbed sine-Gordon equation and proves stability with high accuracy depending on the smoothness of the perturbation.
Findings
01
Solutions remain close to the virtual solitary manifold over long times.
02
The virtual manifold is constructed through iterative deformation.
03
Stability and accuracy improve with the differentiability of the perturbation.
Abstract
We study the perturbed sine-Gordon equation θtt−θxx+sinθ=F(ε,x), where F is of differentiability class Cn in ε and the first k derivatives vanish at 0, i.e., ∂εlF(0,⋅)=0 for 0≤l≤k. We construct implicitly a virtual solitary manifold by deformation of the classical solitary manifold in n iteration steps. Our main result establishes that the initial value problem with an appropriate initial state εn-close to the virtual solitary manifold has a unique solution which follows up to time 1/(C~ε2k+1) and errors of order εn a trajectory on the virtual solitary manifold. The trajectory on the virtual solitary manifold is described by two parameters which satisfy a system of ODEs. In contrast to previous works our stability result yields arbitrarily…
{F~(ε,ξ,x):=F(ε,x)χ(ξ),where χ∈C∞(R),χ(ξ)=1 for ∣ξ∣≤∣ξs∣+3 and χ(ξ)=0 for ∣ξ∣≥∣ξs∣+4.
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.
Full text
\newpagestyle
main
\headrule\sethead[0][\chaptertitle][] \sectiontitle0
We study the perturbed sine-Gordon equation
θtt−θxx+sinθ=F(ε,x),
where F is of differentiability class Cn in ε and the first k derivatives vanish at [math], i.e., ∂εlF(0,⋅)=0 for 0≤l≤k.
We construct implicitly a virtual solitary manifold by deformation of the classical solitary manifold in n iteration steps.
Our main result establishes that the initial value problem with an appropriate initial state εn-close to the virtual solitary manifold has a unique solution which follows up to time 1/(C~ε2k+1) and errors of order εn a trajectory on the virtual solitary manifold. The trajectory on the virtual solitary manifold is described by two parameters which satisfy a system of ODEs.
In contrast to previous works our stability result yields arbitrarily high accuracy as long as the perturbation F is sufficiently often differentiable.
1 Introduction
The perturbed sine-Gordon equation
[TABLE]
is a Hamiltonian evolution equation with Hamiltonian given by
[TABLE]
and the symplectic form given by
[TABLE]
where
[TABLE]
In first order formulation (1) can be written as a system:
[TABLE]
The unperturbed sine-Gordon equation (F(ε,x)=0),
admits soliton solutions
(θ0(ξ(t),u(t),x)ψ0(ξ(t),u(t),x)),
where
[TABLE]
Here the functions (θ0,ψ0) are defined by
[TABLE]
where
[TABLE]
and θK satisfies θK′′(x)=sinθK(x) with boundary conditions θK(x)→(2π0) as x→±∞.
The states (θ0(a,v,⋅)ψ0(a,v,⋅)) form the classical two-dimensional solitary manifold
[TABLE]
In the present paper, we assume that the perturbation term F in (1) is of differentiability class Cn in ε and that its first k derivatives vanish at 0, i.e.,
[TABLE]
We construct
a virtual solitary manifold Snε, which is adjusted to the perturbation F.
The construction can be thought of as a successive distortion of the classical solitary manifold S0. It is based on an iteration scheme composed
of n steps, where in each iteration step a specific PDE will be solved implicitly.
In the last iteration step we obtain an implicit solution (θnε(ξ,u,x),ψnε(ξ,u,x)) which defines
the virtual solitary manifold
[TABLE]
We consider for ξs∈R and ε≪1 the Cauchy problem
[TABLE]
with initial data εn-close to the virtual solitary manifold Snε, i.e.,
[TABLE]
Further, we suppose that
(v(0,⋅),w(0,⋅)) is symplectic orthogonal to the tangent space of Snε at
(θnε(ξs,us,⋅),ψnε(ξs,us,⋅))
and that the smallness assumption
[TABLE]
is satisfied.
Our main theorem shows that, under the mentioned assumptions on the initial data, the Cauchy problem (7) has a unique solution (θ,ψ) for times
[TABLE]
which may be written in the form
[TABLE]
Furthermore, the solution remains εn-close to the manifold Snε, i.e.,
[TABLE]
and the parameters (ξˉ(t),uˉ(t)) satisfy the ODEs
[TABLE]
with initial data
ξˉ(0)=ξs,uˉ(0)=us, where λnε is defined implicitly. The time scale (8)
is nontrivial and the
parameters ξˉ,uˉ describe
a fixed nontrivial perturbation of the uniform linear motion
as ε→0 if the perturbation F satisfies condition (28) mentioned below.
This result yields a fairly accurate description of the solution (θ,ψ) to the Cauchy problem (7), since we are able to control the dynamics on the virtual solitary manifold Snε by the ODEs (10) and the dynamics of the transversal component (v(t,⋅),w(t,⋅)) by the upper bound on its norm (9).
The higher the differentiability class Cn
of the perturbation F the higher is the accuracy of our stability statement.
The time scale of the result is the larger
the more first derivatives of F vanish at 0.
A precise statement is found in section 2.
Let us mention related works and give some background to our paper.
Orbital stability of soliton solutions under perturbations of the initial data has been proven for the (unperturbed) sine-Gordon equation (see [HPW82], [Stu12, Section 4]).
In [Stu92] D. M. Stuart investigated the perturbed sine-Gordon equation,
[TABLE]
where the perturbation g=g(θ) is a smooth function such that g0(Z)=g(θK(Z))∈L2(dZ) and ε≪1. He proved that there exists T∗=O(ε1) such that the corresponding initial value problem with initial data
[TABLE]
has a unique solution of the form
[TABLE]
where θ~∈C([0,T∗],H1),θT∈C([0,T∗],L2) and
[TABLE]
Here p~,u~,C~,dTdu~,dTdC~,∣θ~∣H1(R) are bounded independent of ε and u0,C0 are solutions of certain explicitly given modulation equations.
The proof is based on an orthogonal
decomposition of the solution into an oscillatory part and a one-dimensional
"zero-mode" term.
The sine-Gordon equation arises in various physical phenomena such as dynamics of long Josephson junctions [ZHQ95, KM89], dislocations in crystals
[FK39], waves in ferromagnetic materials [Mik78], etc. In [Sky61] T. H. R. Skyrme proposed the equation to model elementary particles. Dynamics of solitons under constant electric field were examined numerically in [IC79].
In the present paper, we investigate virtual solitons in the presence of a time independent electric field F(ε,x), which is a physically relevant problem.
There are also many other long (but finite)-time results for different equations with external potentials such as [FGJS04, JFGS06, HZ07, Hol11].
For instance, in [HZ08] J. Holmer and M. Zworski considered the Gross-Pitaevskii equation
[TABLE]
with a slowly varying smooth potential V(x)=W(εx), where W∈C3(R,R).
They proved that up to time εlog(1/ε) and errors of size ε2 in H1, the solution is a soliton evolving according to the classical dynamics of a natural effective Hamiltonian.
A common method for investigating the interaction of solitons with external potentials is
to decompose
the dynamics
in a neighbourhood of the classical solitary manifold and to
apply Lyapunov-type arguments afterwards.
Beside other main ingredients this approach has been chosen in [FGJS04, JFGS06, HZ07, HZ08, Hol11].
In the present paper we extend this method for the perturbed sine-Gordon equation by
introducing, the adjusted to the perturbation F, virtual solitary manifold.
We decompose the dynamics in a part on the virtual (rather than classical) solitary manifold and a transversal part, before we proceed with a customized Lyapunov method.
Utilizing the virtual solitary manifold is crucial for the high accuracy in our stability result.
In the broadest sense, a similar technique has been used in [HL12, Section 4, Section 5] for the NLS equation.
In the present paper,
in the iteration scheme for the construction of the virtual solitary manifold, each iteration corresponds to a correction of the classical solitary manifold.
These corrections are represented by approximate solutions of a specific PDE (24) mentioned below,
where the accuracy of the approximate solutions increases with each iteration step by order 1 in ε.
So in [HL12, Section 4, Section 5],
the solitary manifold has been corrected once, which corresponds to the correction executed in the first iteration in the construction of the virtual solitary manifold in the present paper.
In our approach, the construction of the virtual solitary manifold relies on the implicit function theorem.
However, the correction in [HL12, Section 4, Section 5] was done in the form
of a direct asymptotic expansion
and not by employing the implicit function theorem.
In that sense our paper presents a new point of view.
A further remarkable point is that unlike [HL12, Section 4, Section 5] we are able to carry out arbitrarily many corrections in the case of the sine-Gordon equation
as long as the perturbation F is sufficiently often differentiable in ε.
The accuracy of our stability result is the higher the more corrections are possible, wheras the number of possible corrections is determined by the differentiability class of the perturbation F.
We abstained from considering perturbations of type εW(εx) for the following reason. In our approach we do need the assumption that the perturbation F is differentiable with respect to ε, but
there does not exist a function W=0,W∈L2(R)
such that the mapping
ε↦εW(ε⋅)
is differentiable in L2(R).
Further results on long time soliton asymptotics and orbital stability for different equations can be found in
[Wei86, Ben76, Bon75, MP12, SW90, BP92, IKV12, KMMn17, CMnPS16].
This paper is based on [Mas16, Part IV], where many of the computations are presented in greater detail.
Now let us comment on our techniques. The virtual solitary manifold is constructed by the following iteration scheme.
Let firstly ε↦F~(ε) be a general function of differentiability class Cn
mapping into a specific
Sobolev space
such that F~(ε) depends on (ξ,x)
and
F~(0)=0.
F~ will be specified later.
The function (θ0,ψ0), given by (4), is a solution of
[TABLE]
which is the equation characterizing the classical solitons.
In the first iteration step we modify
G0(θ,ψ)=0
by introducing an
additional unknown variable λ and
adding some terms involving (θ0,ψ0) and F~.
The modified equation is of the form
[TABLE]
Here and in the following iterations the functions θ,ψ depend on (ξ,u,x) and
λ depends on (ξ,u).
We solve G1ε(θ,ψ,λ)=0 implicitly for (θ,ψ,λ) in terms of ε and denote the solution by (θ1ε,ψ1ε,λ1ε).
In the next iteration step we modify G1ε(θ,ψ,λ)=0
by adding some terms involving (θ1ε,ψ1ε) and solve the modified equation of the form
[TABLE]
implicitly for (θ,ψ,λ) in terms of ε.
Due to the assumption that ε↦F~(ε)
is of differentiability class Cn,
it is possible to iterate this
procedure until we obtain in the nth step
the equation
[TABLE]
We solve Gnε(θ,ψ,λ)=0 implicitly for (θ,ψ,λ) in terms of ε and
denote the solution by (θnε,ψnε,λnε).
The
existence of the implicit solutions ε↦(θjε,ψjε,λjε) for 1≤j≤n
is ensured by the implicit function theorem.
In the actual proof,
we consider for functional analytic reasons the translated maps and solve the equations
[TABLE]
for (θ^,ψ^,λ) in terms of ε. This is caused, among others, by the fact that θ0(ξ,u,x)→0 as ∣x∣→∞ for fixed ξ and u.
We denote the solutions to the equations G~jε(θ^,ψ^,λ)=0 by (θ^jε,ψ^jε,λjε), where (θjε,ψjε,λjε)=(θ0+θ^jε,ψ0+ψ^jε,λjε).
The application of the implicit function theorem
relies on the fact that
(0,0,0,0)
solves all equations in a particular point, i.e.,
G~j0(0,0,0)=0.
As a consequence of the construction, the
solution obtained in the jth iteration
ε↦(θjε,ψjε,λjε) solves the equation
[TABLE]
up to errors of order εj+1 for 1≤j≤n.
In order to define the virtual solitary manifold we apply this iteration scheme on a specific F~, which is
a truncated version of the perturbation term F in (7), given by
[TABLE]
From now on we denote by (θnε,ψnε,λnε) the solution obtained in the nth iteration by application of the iteration scheme on the specific F~ given by (25).
The first two components of
(θnε,ψnε,λnε) define the virtual solitary manifold (6).
A further consequence of the construction and assumption (5) is that the functions
(θnε(ξ(t),u(t),x),ψnε(ξ(t),u(t),x))
solve the perturbed sine-Gordon equation (3)
up to an error of order εn+k+1 as long as (ξ(t),u(t)) satisfy the ODE system
ξ˙(t)=u(t), u˙(t)=λnε(ξ(t),u(t)).
We call these approximate solutions virtual solitons. In the further proof they play a role which is comparable to that of classical solitons, for instance, in the proof of orbital stability (F=0) of classical solitons (see [Stu12, Section 4]).
The idea of deforming the classical solitary manifold and
utilizing thereby implicitly defined functions appears in [Stu12] with the purpose of rewriting the Hamiltonian in a neighbourhood of the manifold of virtual solitons (see [Stu12, Section 3]).
The virtual solitons in our paper
and the corresponding manifold (6) are defined by
equations and an iteration scheme that were not considered in [Stu12].
The existence of a local solution of the Cauchy problem (7) follows from the contraction mapping theorem.
In the following approach we
derive some bounds which
imply that the local solution is continuable and that estimate (9) is satisfied on the relevant time scale.
We decompose the solution of (7)
into a point on the virtual solitary manifold Snε and a transversal component, i.e.,
[TABLE]
where the parameters (ξ(t),u(t)) are chosen in such a way that the transversal component (v(t,⋅),w(t,⋅)) is symplectic orthogonal to the tangent space of Snε at the corresponding point.
This symplectic decomposition is possible in a small uniform distance to the virtual solitary manifold due to the implicit function theorem.
The energy
[TABLE]
and the momentum
[TABLE]
are conserved quantities of the unperturbed sine-Gordon equation. We make use of this fact and achieve control over the transversal component of the solution (v,w) by utilizing an almost conserved
Lyapunov function, given by
[TABLE]
where (v,w) and (ξ,u) are such as in (26). Lε is the quadratic part of
[TABLE]
The Lyapunov function is bounded from below in terms of ∣v(t,⋅)∣H1(R)2+∣w(t,⋅)∣L2(R)2,
which is a consequence of symplectic orthogonality in decomposition (26) and of
spectral properties of the operator
−∂Z2+cosθK(Z).
The parameters (ξ,u)
satisfy ODEs (10) up to errors of order εn+k+1, which goes
back (among others) to the
construction of the virtual solitary manifold, especially to the fact that (θnε,ψnε,λnε) solves (22) with F~ given by (25).
This property of
(ξ,u) and once again the mentioned fact about (θnε,ψnε,λnε) allow us to
control the Lyapunov function from above.
Therefore we are able to estimate the norm of the transversal component (v,w) and obtain ultimately bound (9).
Using Gronwall’s lemma we pass from the approximate equations for the parameters (ξ,u) to the exact ODEs (10).
Finally let us explain
under which conditions
the result provides a nontrivial dynamics on the virtual solitary manifold
as ε→0.
The linearization of
(θ^,ψ^,λ)↦G~nε(θ^,ψ^,λ)
carried out at
(θ^,ψ^,λ)=(0,0,0), ε=0
is invertible
and we denote the linearization by
[TABLE]
Thus there exist functions (θˉ,ψˉ,λˉ) such that the (k+1)th derivative of F~ with respect to ε, evaluated at ε=0, can be written in the form
[TABLE]
Here the functions θˉ,ψˉ depend on (ξ,u,x) and λˉ depends on (ξ,u).
The ODEs (10) can be rescaled in time by introducing s=εβ(k)t with β(k)=2k+1,
ξ^(s)=ξˉ(s/εβ(k)),
and
u^(s)=εβ(k)1uˉ(s/εβ(k))
such that the
corresponding transformed ODEs have the form
[TABLE]
As ε→0, the transformed ODEs converge to ODEs that describe a fixed nontrivial perturbation of the uniform linear motion if
the next condition is satisfied:
[TABLE]
This
is for the following reason.
Due to (5) ∂εlF~(0)=0 for 1≤l≤k and differentiation of Gnε(θnε,ψnε,λnε)=0
with respect to ε yields (see proof of theorem 3.4):
[TABLE]
Using invertibility of M21, condition (28) and the fact that λn0=0
it follows that 0=λnε=O(εk+1), which implies the claim.
The paper is organized as follows. In section 2, we formulate the main result. In section 3, we construct the virtual solitary manifold.
We prove in section 4 that in a uniform distance to the virtual solitary manifold a decomposition into symplectically orthogonal components is possible. The existence of a local solution (θ,ψ) with initial state close to the virtual solitary manifold is established in section 5. In section 6, we derive modulation equations for the parameters that describe the position on the manifold. We introduce a Lyapunov function and compute its time derivative in section 7. A lower bound on the Lyapunov function is proved in section 8. In section 9, we prove our main result, theorem 2.2.
Some preliminary decompositions are showed in Appendix A. These decompositions are used in Appendix B, where we prove that the linearizations considered in section 3 are invertible.
Notation and Conventions
For a Hilbert space H its inner product is denoted by ⟨⋅,⋅⟩H, the orthogonal complement of a closed subspace M
in H by M⊥,H, the orthogonal projection on M by (⋅)M and the span of v1,…,vp∈H by ⟨v1,…,vp⟩ .
For functions λ depending on (ξ,u) and functions θ depending on (ξ,u,x) the notation λ(ξ,u)=λ(u)(ξ), θ(ξ,u,x)=θ(u)(ξ,x) is used. γ without an argument denotes always γ(u). Occasionally we drop the dependence of functions on certain variables. We also denote occasionally by ∥⋅∥ the norm of an operator and drop the spaces in the notation. We write Lx2(R),Hξ,xk(R2) and so on for the Lebesgue and Sobolev spaces when we wish to emphasize the variables of integration.
2 Main Result
To formulate our result precisely, we need some definitions.
Definition 2.1.
Let α,n∈N and u∗>0. Let us denote by I(u∗):=[−u∗,u∗].
(a)
Hk,α(R)* denotes the weighted Sobolev space of functions with finite norm*
[TABLE]
(b)
Hk,α(R2)* denotes the weighted Sobolev space of functions with finite norm*
[TABLE]
(c)
\underaccent{\bar}{ Y}^{\alpha}* is the space H3,α(R2)⊕H2,α(R2)⊕H2,α(R)
with the finite norm*
*For l∈N and 0<U<u∗ we introduce the parameter area
*
[TABLE]
where V(l,U,u∗):=lu∗−U.
The weighted Sobolev spaces in definition 2.1 (a), (b) are defined as in [Kop15]. We are now ready to state our main result.
Theorem 2.2.
Let n,k∈N,
n≥1, k+1≤n.
Assume that ξs∈R, F∈Cn((−1,1),H1,1(R)) and ∂εlF(0,⋅)=0 for 0≤l≤k.
Then there exist ε0,u∗,C~>0 and a map
[TABLE]
of class Cn such that the following holds.
Let ε∈(0,ε0) and 0<U<u∗. Consider the Cauchy problem
[TABLE]
*where (θnε,ψnε,λnε)=(θ0+θ^nε,ψ0+ψ^nε,λnε) with (θ0,ψ0) given by (4)
such that the following assumptions are satisfied:
*
(a)
∣us∣≤C~ε2k+1;
(b)
Nε(θ(0,⋅),ψ(0,⋅),ξs,us)=0,
where
Nε=(N1ε,N2ε):L∞(R)×L2(R)×Σ(2,U,u∗)→R2 is given by
[TABLE]
and the symplectic form Ω is given by (2);
(c)
∣v(0,⋅)∣H1(R)2+∣w(0,⋅)∣L2(R)2≤ε2n, where (v(0,⋅),w(0,⋅)) is given by
(30).
Then the Cauchy problem defined by (30)
has a unique solution on the time interval
[TABLE]
The solution may be written in the form
[TABLE]
where v,w,
have regularity
(v(t),w(t))∈C1([0,T],H1(R)⊕L2(R))
and ξˉ,uˉ solve the ODEs
[TABLE]
with initial data
ξˉ(0)=ξs,uˉ(0)=us
such that
[TABLE]
The constant C~ depends on F and ξs.
The parameters ξˉ,uˉ describe
a fixed nontrivial perturbation of the uniform linear motion
as ε→0 if condition (28) is satisfied.
3 Construction of the Virtual Solitary Manifold
3.1 Iteration Scheme
Let α,n∈N.
In this subsection,
we establish
the iteration scheme presented in the introduction.
We implement the scheme for a general function F~:(−1,1)→H1,α(R2),ε↦F~(ε) of class Cn that satisfies F~(0)=0.
For being able to apply the implicit function theorem in the proof of existence of iterative solutions, we need to show that the corresponding linearizations of
(θ^,ψ^,λ)↦G~jε(θ^,ψ^,λ) carried out at
(θ^,ψ^,λ)=(0,0,0), ε=0
are invertible (G~j given by (16)-(23)).
This is done in the following proposition, which is a main ingredient in the construction of the virtual solitary manifold.
We start with a definition.
Definition 3.1.
(a)
\underaccent{\bar}{ Z}^{\alpha}* is the space H2,α(R2)⊕H1,α(R2)
with the finite norm*
[TABLE]
(b)
Znα=Znα(u∗)* is the space \bigg{\{}z=(v,w)\in C^{n}(I(u_{*}),\underaccent{\bar}{ Z}^{\alpha}):\|z\|_{Z_{n}^{\alpha}(u_{*})}<\infty\bigg{\}}\,
with the finite norm*
[TABLE]
(c)
We denote by
t1(ξ,u,x):=(∂ξθ0(ξ,u,x)∂ξψ0(ξ,u,x))
and by
t2(ξ,u,x):=(∂uθ0(ξ,u,x)∂uψ0(ξ,u,x)), where u∈(−1,1),ξ,x∈R.
Recall that the spaces \underaccent{\bar}{ Y}^{\alpha}, Ynα(u∗) were defined in section 2. We set Ynα:=Ynα(u∗).
Proposition 3.2.
There exists uα>0 such that the operator
Mnα:Ynα(u∗)→Znα(u∗),(θ,ψ,λ)↦Mnα(θ,ψ,λ),
given by
[TABLE]
is invertible if u∗<uα.
For proof see appendix A and appendix B.
We formalize the iteration scheme in the following theorem. The maps G~j
are defined on spaces of different regularity in u such that the regularity of the spaces decreases with increasing j, which ensures well-definedness of G~j.
Theorem 3.3.
*Let J=(−1,1), u∗<uα
and let F~:J→H1,α(R2),ε↦F~(ε) be a Cn function such that F~(0)=0.
Let G~1 be given by
*
[TABLE]
where G1 is defined by (16). Then there exists ε∗>0 and
a map
[TABLE]
of class Cn such that
G~1ε(θ^1ε,ψ^1ε,λ1ε)=0.
Let G~2 be given by
[TABLE]
where G2 is defined by (19) with (θ1ε,ψ1ε,λ1ε)=(θ0+θ^1ε,ψ0+ψ^1ε,λ1ε). Then there exists
a map
[TABLE]
of class Cn such that
G~2ε(θ^2ε,ψ^2ε,λ2ε)=0.
This process can be continued successively to arrive at G~n be given by
[TABLE]
where Gn is defined by (22) with (θn−1ε,ψn−1ε,λn−1ε)=(θ0+θ^n−1ε,ψ0+ψ^n−1ε,λn−1ε).
Ultimately there exists
a map
[TABLE]
of class Cn such that
G~nε(θ^nε,ψ^nε,λnε)=0
and we set (θnε,ψnε,λnε)=(θ0+θ^nε,ψ0+ψ^nε,λnε).
Proof 0.
We skip u∗ in the notation.
Notice that
G~10(0,0,0)=G10(θ0,ψ0,0)=0.
The derivative of
G~1:J×Yn+1α→Zn+1α
with respect to (θ^,ψ^,λ) evaluated at (ε,θ^,ψ^,λ)=(0,0,0,0)
is
Mn+1α, which is invertible due to proposition 3.2. By the implicit function theorem there exists a
ε1∗>0 and a map
[TABLE]
of class Cn such that
G~1ε(θ^1ε,ψ^1ε,λ1ε)=0.
We continue successively this process until we
obtain
that the derivative of
G~n:J×Y2α→Z2α
is
M2α.
This yields by using
the same argument combined with G~n0(0,0,0)=0 that there exists
εn∗>0 and a map
[TABLE]
of class Cn such that
G~nε(θ^nε,ψ^nε,λnε)=0.
We set ε∗=min{ε1∗,ε2∗,…,εn∗}.
□
In the following theorem we state the properties of the nth iterative solution from theorem 3.3.
Theorem 3.4.
Let the assumptions of theorem 3.3 hold. Then the following relations are satisfied.
(a)
(∂εjθn−10,∂εjψn−10,∂εjλn−10)=(∂εjθn0,∂εjψn0,∂εjλn0)* for j=0,…,n−1, where n≥2.*
The derivatives of the iterative solutions coincide at [math] in the following way:
(∂εjθ10,∂εjψ10,∂εjλ10)=(∂εjθ20,∂εjψ20,∂εjλ20) for j=0,1;(∂εjθ20,∂εjψ20,∂εjλ20)=(∂εjθ30,∂εjψ30,∂εjλ30) for j=0,1,2 and so on up to the identities for the nth iterative solution stated in theorem 3.4 (a).
(a) The claim can be proved by induction on n. We show the induction step.
Assume that the claim is true for all integers less than or equal to n−2.
Let 0≤j≤n−1. The fact that
the solutions
from theorem 3.3 satisfy
[TABLE]
and that the injections
H1(R)⊂L∞(R),H2(R2)⊂L∞(R2)
are continuous [Bre11, Corollary 9.13] yields the justification for using the Leibniz’s formula and Faà di Bruno’s formula.
Thus we obtain for the
j-th derivatives with respect to ε, evaluated at
ε=0:
[TABLE]
[TABLE]
*where l=l1+l2+…+lj and the sum is taken over all l1,l2,…,lj for which l1+2l2+…+jlj=j.
Subtracting (34) from (33) yields the claim due to proposition 3.2.
(b) follows from (a), theorem 3.3 and Taylor’s formula.
□*
3.2 Virtual Solitary Manifold
From now on we set α:=1.
In this subsection, we apply theorem 3.3 on a specific F~ and define the virtual solitary manifold by the solution obtained in the nth iteration.
Definition 3.6.
Let F,ξs be from theorem 2.2 and Ξ:=Ξ(ξs):=∣ξs∣+3.
We set
F~(ε,ξ,x):=F(ε,x)χ(ξ),
where χ is a smooth cutoff function
with χ(ξ)=1 for ∣ξ∣≤Ξ and χ(ξ)=0 for ∣ξ∣≥Ξ+1.
The next lemma follows immediately from theorem 3.4 and from the assumptions on F in theorem 2.2.
F~∈Cn((−1,1),H1,1(R2))* and
∂εlF~(0,⋅,⋅)=0 for 0≤l≤k.*
(c)
00λnεY21(u∗)=O(εk+1).**
Lemma 3.8.
Let v∈H1(R2). Then there exists b>0 such that
[TABLE]
Proof 0.
This follows from applying Morrey’s embedding Theorem to the variable ξ.
□
We solve iteratively the equations in
theorem 3.3 with the specific F~(ε,ξ,x):=F(ε,x)χ(k,ξ) from definition 3.6
and define by the nth solution (θnε,ψnε,λnε) the virtual solitary manifold Snε.
We utilize the truncated version of F rather than F itself in order to make sure that the maps G~j in theorem 3.3 are well defined.
theorem 3.3 is applicable to F~
due to lemma 3.7.
Definition 3.9.
Let u1 be from proposition 3.2 (case α=1). We fix a specific u∗ such that 0<u∗<u1.
Let F~ be from definition 3.6.
Let ε∗>0 be the constant and let (θnε,ψnε,λnε) be the nth solution obtained from application of theorem 3.3 to F~.
We set
[TABLE]
and call Snε the virtual solitary manifold.
Remark 3.10.
**
(a)
*From now on we
denote by (θnε,ψnε,λnε)=(θ0+θ^nε,ψ0+ψ^nε,λnε) always the *nth solution utilized in definition 3.9.
(b)
The vectors
[TABLE]
are tangent vectors of the manifold Snε at the point (θnε(ξ,u,⋅),ψnε(ξ,u,⋅)) and form a basis of the tangent space at this point.
4 Symplectic Orthogonal Decomposition
Let from now on u∗,ε∗ be always from definition 3.9 and let U be fixed such that 0<U<u∗.
We consider V(l,U,u∗),Σ(l,U,u∗) introduced in definition 2.1 (e)
and the function Nε:L∞(R)×L2(R)×Σ(2,U,u∗)→R2 defined as in (31) for these specific u∗,U. For simplicity of further notation we set Σ(l,U,u∗)=Σ(l) and V(l,U,u∗)=V(l).
In this chapter we will choose ε0 sufficiently small and consider ε∈(0,ε0].
We show that if (θ,ψ)∈L∞(R)⊕L2(R) is close enough (in the L∞(R)⊕L2(R) norm) to the region
[TABLE]
of the virtual solitary manifold
Snε, then there exists a unique (ξ,u)∈Σ(2) such that we are able to decompose the solution
[TABLE]
in a point on the virtual solitary manifold (θnε(ξ,u,⋅),ψnε(ξ,u,⋅)) and a transversal component (v(⋅),w(⋅)), which
is symplectic orthogonal to the tangent vectors
t1,nε(ξ,u,⋅)
and
t2,nε(ξ,u,⋅)
at the corresponding point of the manifold Snε,
i.e., the orthogonality condition
[TABLE]
is satisfied.
We prove that the symplectic decomposition is possible in a small uniform distance to the manifold Snε, where the distance might depend on ε0 but does not depend on ε.
Remark 4.1.
*In theorem 3.3 we have solved the equations defining (θnε,ψnε,λnε) in weighted spaces. One of the reasons for working in weighted Sobolev spaces was to make sure that Nε:L∞(R)×L2(R)×Σ(2)→R2 is
well defined.
*
We start with a definition and some elementary lemmas which will be used later.
A straight forward computation yields the following lemma.
Lemma 4.3.
Let ε∈(0,ε∗). It holds that
[TABLE]
Lemma 4.4.
Let ε0>0 be sufficiently small.
There exist constants
c=c(U)>0,C=C(U)>0,
such that ∀ε∈(0,ε0],(ξ,u)∈R×[−U−V(2),U+V(2)]:
[TABLE]
Proof 0.
Using lemma 3.8 and
continuity of
ε↦(θ^nε,ψ^nε,λnε)
(see theorem 3.3) we obtain for sufficiently small ε0:
∀ε∈(0,ε0],(ξ,u)∈R×[−U−V(2),U+V(2)]:∣knε(ξ,u)∣<m/2,
which implies the claim.
□
The next lemma provides that the symplectic decomposition described above is possible. In the proof we will take derivatives of (θnε,ψnε) up to second order with respect to ξ and u. This was the reason for solving, in section 3, the equations defining (θnε,ψnε,λnε) in spaces of higher regularity in ξ and u.
Lemma 4.5.
Let ε0>0 be sufficiently small.
Let
[TABLE]
There exists r>0 such that if ε∈(0,ε0]
and p≤r then for any (θ,ψ)∈OU,pε there exists a unique
(ξ,u)∈Σ(2) such that
[TABLE]
and the map
(θ,ψ)↦(ξ(θ,ψ),u(θ,ψ))
is in C1(OU,pε,Σ(2)).
Proof 0.
*Let ε0∈(0,ε∗).
We will specify ε0 later in this proof.
Let ε∈(0,ε0]. Notice that the map
ε↦(θ^nε,ψ^nε,λnε)
from theorem 3.3
is continuous and it holds
Σ(4)⊂Σ(3)⊂Σ(2).
Consider (ξ0,u0)∈Σ(3).
lemma 4.3 yields that
*
[TABLE]
Using lemma 3.8 we obtain for sufficiently small ε0 for all
ε∈(0,ε0]:
∣knε(ξ0,u0)∣≤2m
and thus
[TABLE]
*We prove that there exist r>0,δˉ>0,ε0>0
such that
∀ε∈(0,ε0],(ξ0,u0)∈Σ(3) there exist balls
Br(θnε(ξ0,u0,⋅),ψnε(ξ0,u0,⋅))⊂L∞(R)⊕L2(R),
Bδˉ(ξ0,u0)⊂Σ(2),
and a map
*
[TABLE]
such that
Nε(θ,ψ,Tξ0,u0ε(θ,ψ))=0
on
Br(θnε(ξ0,u0,⋅),ψnε(ξ0,u0,⋅)). Therefore we refer to [Dei85, Theorem 15.1] and check their proof of the implicit function theorem, whereas we show that r and δˉ do not depend on ε and on (ξ0,u0).
We introduce
[TABLE]
Notice that
Nˉξ0,u0ε(0,0,0,0)=(0,0).
We set
Kξ0,u0ε:=D(ξ,u)Nˉξ0,u0ε(0,0,0,0)
and
[TABLE]
which is well defined due to (36).
Due to lemma 3.8 it holds for a sufficiently small ε0 that
[TABLE]
In this proof we denote by ∥⋅∥ the maximum row sum norm of a 2×2 matrix induced by the maximum norm ∣⋅∣∞ in R2.
We claim that ∃k∈(0,1),δˉ>0,ε0>0∀ε∈(0,ε0],(ξ0,u0)∈Σ(3)∀((θ,ψ),(ξ,u))∈Bδˉ(0)×Bδˉ(0):∥D(ξ,u)Sξ0,u0ε(θ,ψ,ξ,u)∥≤k<1.
Due to (35) it holds that
[TABLE]
*The claim follows by using lemma 3.8, lemma 4.4 and estimating each entry of D(ξ,u)Sξ0,u0ε(θ,ψ,ξ,u), for instance:
*
[TABLE]
Similarly as above one shows that ∃r≤δˉ,ε0>0∀ε∈(0,ε0],(ξ0,u0)∈Σ(3)∀(θ,ψ)∈Br(0):∣Sξ0,u0ε(θ,ψ,0,0)∣∞<δˉ(1−k),
which completes the proof.
□
5 Existence of Dynamics and the Orthogonal Component
We argue similar to [Stu98, Proof of theorem 2.1]. Let ε0 be from lemma 4.5 and ε∈(0,ε0]. In order to make use of existence theory we
consider the problem
[TABLE]
By [Mar76, Theorem VIII 2.1, Theorem VIII 3.2
] there exists a local solution (see also [Stu98, Proof of theorem 2.1], [Stu92, p.434
]), where
[TABLE]
(θ,ψ) given by θ(t,x)=vˉ(t,x)+θnε(ξs,us,x) and ψ(t,x)=wˉ(t,x)+ψnε(ξs,us,x) solves obviously locally the Cauchy problem (30) and (θ,ψ)∈C1([0,Tloc],L∞(R)⊕L2(R)) due to Morrey’s embedding theorem.
We are going to obtain a bound in section 9 which will imply that the local solutions are indeed continuable.
So from now we assume that (vˉ,wˉ)∈C1([0,T],H1(R)⊕L2(R)) is a solution of (37)-(38) and (θ,ψ) is a solution of (30) such that (θ,ψ)∈C1([0,T],L∞(R)⊕L2(R)), where T>0.
In the following we define, similar to Σ(l,U,u∗), a new parameter area, where the parameter ξ is bounded.
Given (θ,ψ) we choose the parameters (ξ(t),u(t)) according to lemma 4.5 and define (v,w) as follows:
[TABLE]
(v(t,x),w(t,x)) is well defined for t≥0 so small that
∣v(t)∣L∞(R)+∣w(t)∣L2(R)≤r
and
(ξ(t),u(t))∈Σ(4,Ξ),
where r is from lemma 4.5.
We formalize this in the following definition.
Notice that (ξs,us)=(ξ(0),u(0))∈Σ(4,Ξ).
The transversal component (v(t,x),w(t,x)) is well defined for 0≤t≤t∗
since, among others, we choose ε such that
ε∈(0,ε0] with ε0 from lemma 4.5 and the initial data such that ∣v(0)∣L∞(R)+∣w(0)∣L2(R)≤2r with (v(0),w(0)) given by (30).
Lemma 5.3.
Let T=min{t∗,T} and let (v,w) be defined by (39)-(40). Then
(v,w)∈C1([0,T],H1(R)⊕L2(R)).
Proof 0.
This follows by using (39)-(40)
and the fact that (vˉ,wˉ)∈C1([0,T],H1(R)⊕L2(R)), since the difference (θK(γ(u0)(⋅−ξ0))−θK(γ(uˉ)(⋅−ξˉ))) is in L2(R) for all (ξ0,u0),(ξˉ,uˉ)∈R×(−1,1).
□
In the following lemma we point out the relation between F and (θnε,ψnε,λnε). Notice that there appears F instead of F~ in the equation above.
Moreover, we state the rates of convergence of Rnε(ξ,u,⋅) and λnε(ξ,u) which will be needed in the proof of the modulation equations for the parameters (ξ(t),u(t)) in the next section and in the proof of the main result in section 9.
Lemma 5.4.
It holds that for a.e. (ξ,u,x)∈Σ(4,Ξ)×R
[TABLE]
and
∣[Rnε(ξ,u,⋅)]1∣L2(R)=O(εn+k+1),∣[Rnε(ξ,u,⋅)]2∣L2(R)=O(εn+k+1),∣λnε(ξ,u)∣=O(εk+1),∣∂1λnε(ξ,u)∣=O(εk+1),∣∂2λnε(ξ,u)∣=O(εk+1)
uniformly in (ξ,u)∈Σ(4,Ξ).
Proof 0.
The first identity follows due to LABEL:ITnquantitativ and lemma 3.7.
Using LABEL:ITnquantitativ, lemma 3.8 and Morrey’s embedding theorem we obtain for all (ξ,u)∈Σ(4,Ξ):
[TABLE]
and
[TABLE]
The other cases can be treated analogously.
□
We compute the time derivatives of v and w, which will be needed in the following sections.
for times t∈[0,t∗], where R~(v)=O(∣v∣Hx1(R)3) and Rnε(ξ,u,x) is from theorem 3.4 (b).
Proof 0.
By taking the time derivatives of (v,w) and using lemma 5.4, (30) we obtain
[TABLE]
and
[TABLE]
where we have expanded the term sin(θnε(ξ,u,x)+v(x)).
□
6 Modulation Equations
In the following lemma we derive modulation equations for the parameters (ξ(t),u(t)).
Lemma 6.1.
There exists an ε0>0 such that the following statement holds. Let ε∈(0,ε0] and let (v,w) be given by (39)-(40)
with (ξ,u) obtained from lemma 4.5. Let
[TABLE]
Then it holds for t∈[0,t∗] that
[TABLE]
where C depends on F and ξs.
Proof 0.
The technique we use is similar to that in the proof of [IKV12, Lemma 6.2].
Let ε0∈(0,ε∗) with ε∗ from definition 3.9 and let ε∈(0,ε0). Further in the proof we will make some more assumptions on ε0.
We start with some definitions and set
[TABLE]
Now we consider for any (ξˉ,uˉ)∈R×[−U−V(2),U+V(2)], (vˉ,wˉ)∈H1(R)×L2(R) the matrix:
[TABLE]
We use lemma 4.4, lemma 3.8 and Hölder’s inequality similar to the proof of lemma 4.5
and obtain for all (ξˉ,uˉ)∈R×[−U−V(2),U+V(2)], (vˉ,wˉ)∈H1(R)×L2(R):
[TABLE]
where we denote by ∥⋅∥ a matrix norm. Let I=I2 be the identity matrix of dimension 2. Due to (41) we are able to choose ε0>0 such that if ∣vˉ∣H1(R),∣wˉ∣L2(R)≤ε0 then the matrix
[TABLE]
*is invertible by von Neumann’s theorem.
Using
(39)-(40) we express the orthogonality condition Nε(θ,ψ,ξ,u)=0 from lemma 4.5 in terms of (v,w,ξ,u)
and take its derivative with respect to t.
For simplicity of notation, we drop (θ,ψ,ξ,u) and obtain
in matrix form:
*
[TABLE]
where
M=Mnε(ξ,u,v,w), Ω=Ωnε(ξ,u), P1=P1,nε(ξ,u,v,w), P2=P2,nε(ξ,u,v,w),
[TABLE]
and
[TABLE]
If ∣v∣H1(R),∣w∣L2(R)≤ε0 then we obtain
as mentioned above
by von
Neumann’s theorem that
[TABLE]
We make a further assumption on ε0, namely that ε0 should be so small that the convergence rates in lemma 5.4
are satisfied.
Now we consider P1 and P2.
The zeroth-order Taylor’s approximations (in ε) of expressions (42)-(43) and (44)-(45) respectively are
[TABLE]
where Lξ,u
is given in definition A.4.
Integration by parts and symplectic orthogonality yield that these approximations
vanish, which can also be deduced from lemma A.5.
Thus we obtain from lemma 5.4,
and similar arguments as above
[TABLE]
□
7 Lyapunov Function
In this section we introduce the Lyapunov function and calculate its time derivative.
Definition 7.1.
Let (v,w) be given by (39)-(40), with (ξ,u) obtained from lemma 4.5.
We define the Lyapunov function Lε by
[TABLE]
and the auxiliary function L by
[TABLE]
Lemma 7.2.
It holds that
[TABLE]
Proof 0.
We use a similar technique as in the proof of
[KSK97, Lemma 2.1].
We can assume that the initial data of our problem have compact support. This allows us to do the following computations (integration by parts etc.). The claim for non-compactly supported initial data follows by density arguments. We obtain the stated formula by taking the time derivative of (46), where we use
[TABLE]
and
∫∂xv(x)∂x2v(x)+w(x)∂xw(x)dx=0.
□
8 Lower Bound
Here we introduce a functional E and prove a lower bound on E by using symplectic orthogonality combined with functional analytic arguments. This will imply a lower bound on the Lyapunov function Lε,
which will play a key role in the proof of the main result.
Definition 8.1.
For (v,w)∈H1(R)×L2(R), (ξ,u)∈R×(−1,1) we set
[TABLE]
where Z=γ(x−ξ) and vZ(x)=∂Zv(γZ+ξ)=γ1∂xv(x).
A straightforward computation yields the following lemma.
Lemma 8.2.
For (v,w)∈H1(R)×L2(R), (ξ,u)∈R×(−1,1) it holds that
[TABLE]
Recalling the relations (39)-(40) we introduce a notation in order to be able to express the orthogonality conditions in terms of the variables (v,w,ξ,u) instead of the variables
(θ,ψ,ξ,u).
Definition 8.3.
For (v,w)∈H1(R)×L2(R), (ξ,u)∈R×(−u∗,u∗) we set
[TABLE]
Now we prove a lower bound on the functional E.
Lemma 8.4.
Let ε0>0 be sufficiently small.
There exists c>0 such that if ε∈(0,ε0),
(ξ,u)∈[−Ξ,Ξ]×[−U−V(2),U+V(2)]⊂R×(−1,1) and (v,w)∈H1(R)×L2(R)
satisfy
[TABLE]
then
[TABLE]
Proof 0.
We follow closely [Stu12] and [Stu98]. This proof is a slight modification of the proof of [Stu12, Lemma 4.3].
First of all we choose ε0 such that ε0∈(0,ε∗) with ε∗ from definition 3.9. We will specify ε0 later.
Notice that the operator −∂Z2+cosθK(Z) is nonnegative. It has (see [Stu92]) an one
dimensional null space spanned by θK′(⋅) and the essential spectrum [1,∞).
We argue by contradiction and assume that the result claimed is false: ∀j∈N∃εj∈(0,ε0],(ξj,uj)∈[−Ξ,Ξ]×[−U−V(2),U+V(2)],(vˉj,wˉj)∈H1(R)×L2(R):
[TABLE]
This statement is also true for the sequences
vj:=vˉj(∣vˉj∣H1(R)2+∣wˉj∣L2(R)2)−21
and
wj:=wˉj(∣vˉj∣H1(R)2+∣wˉj∣L2(R)2)−21.
Assuming that ∣vj∣L2(R)→j→∞0 we obtain
∣(vj)x∣L2(R)→j→∞0 and
∣wj∣L2(R)→j→∞0. This is a contradiction to the fact that ∣vj∣H1(R)2+∣wj∣L2(R)2=1∀j∈N. By passing to a subsequence we may assume (without loss of generality) that there exists δˉ>0 such that
[TABLE]
Since (vj,wj) is bounded in H1(R)×L2(R) we may assume that vj⇀H1(R)v and wj⇀L2(R)w by taking subsequences.
Due to Rellich’s theorem we may assume again by passing to subsequences that vj→L2(Ω)v, where Ω⊂R is bounded and open. Passing to a further subsequence we assume almost everywhere convergence.
Due to the fact that
[TABLE]
and that −∂Z2+cosθK(Z) is a nonnegative operator we obtain the estimate
[TABLE]
where we used integration by parts and substitution.
We may extract a subsequence such that uj→Ru, ξj→Rξ and εj→Rε^.
It follows from (47) and from the previous estimate that
[TABLE]
for a sufficiently large r~. As a consequence, (48) and the strong convergence on bounded domains yield
∫{x∈R:∣x∣≤r~}v2(x)dx≥δˉ, from which it follows that v≡0. Weak convergence and the continuity of
ε↦(θ^nε,ψ^nε,λnε)
(see theorem 3.3)
imply using the triangle inequality that
[TABLE]
and
[TABLE]
Due to (49) we are able to apply Fatou’s lemma for a sufficiently large r~ and obtain
[TABLE]
where we have used that (vj) converges almost everywhere.
The dominated convergence theorem yields:
[TABLE]
(47) together with (51)-(54) imply that
E(v,w,ξ,u)=0.
This yields (v(x),w(x))=η(θK′(γ(u)(x−ξ)),−uγ(u)θK′′(γ(u)(x−ξ))) for some η=0, since v≡0.
Using lemma 3.8 and the notation from definition 4.2 we choose
ε0 sufficiently small so that for all (ξ,u)
[TABLE]
which implies that Nˇ2ε^(v,w,ξ,u)=0. This yields a contradiction to (50).
□
Remark 8.5.
Let (v,w) be given by (39)-(40), with (ξ,u) obtained from lemma 4.5. It holds that
L(t)=E(v(t),w(t),ξ(t),u(t)).
We prove first a version of theorem 2.2 with approximate equations for the parameters (ξ,u).
Theorem 9.1.
Assume that the assumptions
of theorem 2.2
on n,k,ξs and F
are satisfied.
There exist ε0,u∗,C~>0
and a map
[TABLE]
*of class Cn such that the following holds. Let ε∈(0,ε0).
Consider the Cauchy problem
*
[TABLE]
where (θnε,ψnε,λnε)=(θ0+θ^nε,ψ0+ψ^nε,λnε)
and (ξs,us)=(ξ(0),u(0))∈R×(−1,1) such that the following assumptions are satisfied:
(a)
∣us∣≤C~ε2k+1.
(b)
Nε(θ(0,⋅),ψ(0,⋅),ξs,us)=0.
(c)
∣v(0,⋅)∣H1(R)2+∣w(0,⋅)∣L2(R)2≤ε2n, where (v(0,⋅),w(0,⋅)) is given by (56).
Then the Cauchy problem defined by (56)
has a unique solution on the time interval
[TABLE]
The solution may be written in the form
[TABLE]
where v,w,ξ,u have regularity
(v(t),w(t))∈C1([0,T],H1(R)⊕L2(R))
and
(ξ(t),u(t))∈C1([0,T],R×(−1,1))
such that the symplectic orthogonality condition
[TABLE]
is satisfied. There exist positive constants c,C such that
[TABLE]
and
[TABLE]
The constants c,C depend on F and ξs.
Notice that the previous theorem describes the dynamics less precisely than theorem 2.2. However, in the previous theorem
the orthogonality condition is satisfied which does not have to hold in theorem 2.2.
The proof of theorem 9.1 needs some preparation. The existence of ε0>0, u∗>0 and the map (55) is ensured by theorem 3.3.
Now we suppose that (56)
has a solution and we make some assumptions on (v,w) given by (39)-(40) and on (ξ,u) obtained from lemma 4.5. Then the following lemma yields us more accurate information on (v,w) and (ξ,u).
Lemma 9.2.
Assume that the assumptions
of theorem 2.2
on n,k,ξs,F
are satisfied and let 0<δ<1/32.
There exist ε0,C~>0 such that the following statement holds. Let ε∈(0,ε0).
Assume that (56)
has a solution (θ,ψ) on [0,T] such that
[TABLE]
Suppose that
0≤T≤t∗≤T.
Suppose that (v,w) is given by (39)-(40), with (ξ,u) obtained from lemma 4.5 such that
[TABLE]
Then, provided
[TABLE]
there exist c,C>0 such that
(1)
∀t∈[0,T](ξ(t),u(t))∈Σ(5,Ξ),**
(2)
\left|{v}\right|_{L^{\infty}({[{0},{T}]},H^{1}(\mathbb{R}))}^{2}+\left|{w}\right|_{L^{\infty}({[{0},{T}]},L^{2}(\mathbb{R}))}^{2}\leq\frac{1}{\color[rgb]{0,0,0}c}(L(0)+C\varepsilon^{2n}),*
where c is from lemma 8.4 and C depends on F,ξs.*
Remark 9.3.
*Notice that the assumption T≤t∗ yields us the information:
∀t∈[0,T](ξ(t),u(t))∈Σ(4,Ξ).*
Proof 0.
*Choose ε0 sufficiently small, in particular such that the lemmas used below can be applied and such that for any ε∈(0,ε0) the following statement holds: if (v,w)∈H1(R)×L2(R) satisfies ∣v∣H1(R)2+∣w∣L2(R)2≤ε2n−δ then it holds that ∣v∣L∞(R)+∣w∣L2(R)≤2r, where r is from lemma 4.5. This can be ensured by Morrey’s embedding theorem.
Then, using lemma 5.4, it follows that there exists C~>1 such that ∀t∈[0,T]:
[TABLE]
This implies (1) by choosing ε0 small enough and utilizing ∣us∣≤C~ε2k+1.
By using lemma 8.4, lemma 8.2, lemma 7.2 and lemma 5.4 we obtain for times
0≤t≤T≤C~εβ(k)1,
estimate (2):
[TABLE]
□
Theorem 9.4.
Assume that the assumptions
of theorem 2.2
on n,k,ξs and F
are satisfied.
There exists ε0,C~>0 such that the following statement holds. Let ε∈(0,ε0).
Assume that
(56)
has a solution (θ,ψ) on [0,T] such that
[TABLE]
Suppose that
0≤T≤T
and that the assumptions (a),(b),(c) of theorem 9.1 are satisfied.
Then, provided
[TABLE]
it holds that (v,w) given by (39)-(40) is well defined for times [0,T] and
there exists c^>0 such that
The previous theorem implies that the local solution of (56)
discussed in section 5 is indeed continuable up to times 1/(C~εβ(k)) for ε∈(0,ε0). theorem 9.4 and lemma 6.1 yield the approximate equations for the parameters (ξ,u), which conclude the proof of theorem 9.1.
9.2 ODE Analysis
In this subsection we lay the groundwork for passing from the approximate equations for the parameters (ξ,u)
in theorem 9.1 to the ODEs given by (32).
We start with a preparing lemma.
Lemma 9.5.
There exists ε0>0 such that the following statement holds. Let ε∈(0,ε0).
Let β(k)=2k+1.
Let ξ~=ξ~(s), u~=u~(s),
ϵ1=ϵ1(s), ϵ2=ϵ2(s) be C1 real-valued
functions.
Suppose that
[TABLE]
on [0,T] for j=1,2.
Assume that on [0,T],
[TABLE]
Let ξ^=ξ^(s) and u^=u^(s) be C1 real-valued
functions which satisfy the exact equations
[TABLE]
Then
there exists c>0 such that the estimates
[TABLE]
hold on [0,T].
Proof 0.
We follow very closely [HZ08, Lemma 6.1].
We choose ε0 so small that the convergence rates in lemma 5.4 are satisfied for all ε∈(0,ε0). Let ε∈(0,ε0).
Let x=x(s) and y=y(s) be C1
real-valued functions, C≥1, and let (x,y) satisfy the differential
inequalities:
[TABLE]
For z(s)=x2+y2 the following estimate holds
[TABLE]
It follows from Gronwall’s lemma that z(s)≤z(0)e4Cs.
Thus
[TABLE]
Now we recall Duhamel’s formula.
Let X(s):R→R2 be a two-vector function,
X0∈R2 a two-vector, and
A(s):R→(2×2 matrices) a 2×2 matrix function. We consider the ODE system
[TABLE]
and denote its solution by X(s)=S(s,s′)X0 such that
[TABLE]
Let F(s):R→R2 be a 2-vector function. We can express the
solution of the inhomogeneous ODE system
[TABLE]
with initial condition X(0)=0 by Duhamel’s formula
[TABLE]
Let U=u^−u~ and Ξ=ξ^−ξ~. These functions satisfy
[TABLE]
Let
[TABLE]
[TABLE]
We set
[TABLE]
and obtain by Duhamel’s formula:
[TABLE]
We use lemma 5.4 and apply (57) with
[TABLE]
It follows that
[TABLE]
Using (58) we obtain that on [0,T]
[TABLE]
*which yields the claim.
*□
In the following we show the relation between the parameters (ξ,u) selected by the implicit function theorem according to lemma 4.5 and the solutions (ξ^,u^) of the exact ODEs from the previous lemma.
Lemma 9.6.
Assume that the assumptions
of theorem 9.1 are satisfied.
There exists ε0>0 such that the following statement holds. Let ε∈(0,ε0),
β(k)=2k+1 and
s=εβ(k)t,
where
[TABLE]
Let (ξ,u) be the parameters selected according to lemma 4.5 and
(ξ^,u^) from lemma 9.5.
Then there exists c>0 such that
[TABLE]
Proof 0.
We choose ε0 as in theorem 9.1.
Let ε∈(0,ε0) and
theorem 9.1 yields the dynamics with the parameters (ξ,u) selected by the implicit function theorem according to lemma 4.5 on the time interval
0≤t≤(C~εβ(k))−1.
Using lemma 9.6 and the the triangle inequality we can replace (ξ(t),u(t))
with
(ξˉ(t),uˉ(t)):=(ξ^(εβ(k)t),εβ(k)u^(εβ(k)t)). The claim follows after possibly increasing the constant C~ in the proof of theorem 9.1.
□
Appendix A Preliminary Decompositions
Let α,n∈N. Here we prove some decompositions for certain Sobolev spaces on R and on R2.
We start with the spaces on R and prove an orthogonal decomposition of H1(R)⊕L2(R).
We prove the case α=0 that implies the claim. We consider the case (ξ,u)=(0,0). The proof for a general (ξ,u)∈R×(−1,1) works in the same way.
L0,00 is self-adjoint and [math] is an isolated point of σ(L0,00).
l:=L0,00∣H2(R)∩⟨θK′⟩⊥ is self-adjoint and has a bounded inverse (see [HS96, Proposition 6.6]). Notice that ranL0,00=ranl. Let yn=Mxn→L2y. Boundness yields xn=l−1yn→L2l−1y, where l−1 denotes the bounded extension of l−1 on the closure ranl. Since l∗=l is a closed operator (see [HS96, Proposition 4.9]), we obtain l(l−1y)=y.
□
lemma A.2 and the inverse mapping theorem imply the following lemma.
M^ξ,uα∈L(Hξ,u,⊥2,α(R)⊕R,L2,α(R))* and M^ξ,uα is one-to-one and onto.*
Definition A.4.
We define the following operators.
(a)
Lξ,u:H2(R)⊕H1(R)→H1(R)⊕L2(R)* given by*
[TABLE]
(b)
{\hat{\cal L}}_{\xi,u}:\Big{[}H^{2}(\mathbb{R})\oplus H^{1}(\mathbb{R})\Big{]}\cap({\rm ker}\,{\cal L}_{\xi,u})^{\perp,L^{2}\oplus L^{2}}\to H^{1}(\mathbb{R})\oplus L^{2}(\mathbb{R})* given by*
[TABLE]
Lemma A.5 (orthogonal sum).
[TABLE]
Proof 0.
"⊃": clear.
"⊂": Let (vˉ,wˉ)∈H1(R)⊕L2(R).
Orthogonal decomposition of L2(R)⊕L2(R) yields that there exists μ(ξ,u)∈R and (θn,ψn)∈H2(R)⊕H1(R) such that
[TABLE]
since
kerLξ,u∗=⟨Jt1(ξ,u,⋅)⟩ due to lemma A.3.
Hence \begin{pmatrix}v\\
w\end{pmatrix}:={\scriptscriptstyle L^{2}\oplus L^{2}-}\lim_{n\to\infty}{\cal L}_{\xi,u}\begin{pmatrix}\theta_{n}\\
\psi_{n}\end{pmatrix}\in\langle\mathbb{J}t_{1}(\xi,u,\cdot)\rangle^{\perp,L^{2}\oplus L^{2}}\cap\Big{[}H^{1}(\mathbb{R})\oplus L^{2}(\mathbb{R})\Big{]} and thus
(−uγv′(⋅)+γw(⋅))∈⟨θK′(γ(⋅−ξ)⟩⊥,L2=ranL^ξ,u
due to lemma A.3 (α=0).
By setting
{\tilde{\theta}}(x):={[\hat{L}_{\xi,u}]^{-1}}\Big{(}-u\gamma\partial_{x}v(x)+\gamma w(x)\Big{)},
ψ~(x):=−u∂xθ~(x)−v(x) we obtain that (θ~,ψ~)∈H2(R)⊕H1(R) due to lemma A.3 and that
L^ξ,u(θ~ψ~)(kerLξ,u)⊥,L2⊕L2=(vw).
□
We turn to decompositions of Sobolev spaces on R2. Here
lemma A.5 will be needed later in the proofs. We start with a definition and some preparing lemmas.
Unlike the one-dimensional case we consider ξ not as
a fixed parameter anymore, but as a new variable.
The claim follows from lemma A.3 combined with the fact that the operator −∂Z2+cosθK(Z) is nonnegative with essential spectrum [1,∞) and with one
dimensional null space spanned by θK′(⋅) (see [Stu92]).
□
Using that
Lu is self-adjoint and [math] is an isolated point of σ(Lu) one proves
analogously to lemma A.2
the following lemma.
Lemma A.8.
ranLuα* is closed with respect to L2,α(R2).*
Analogously to
lemma A.3 we obtain the next lemma.
Lemma A.9.
(a)
L2,α(R2)=ranLuα⊕L2,αkerLuα.
(b)
L2,α(R2)=ranL^uα⊕L2,αkerLuα.
(c)
L^uα∈L(Hu,⊥2,α(R2),Lu,⊥2,α(R2)).
(d)
[L^uα]−1∈L(Lu,⊥2,α(R2),Hu,⊥2,α(R2)).
(e)
M^uα,M^α* are one-to-one, onto, bounded and the inverse mappings are also bounded.*
Lemma A.10.
Let mα:Yα→Zα,(θ,λ)↦mα(θ,λ) be the linear operator, given by
[TABLE]
Then mα is one-to-one, onto and bounded, i.e., [mα]−1 is bounded.
Proof 0.
mα* is well defined and mα
one-to-one due to
lemma A.9. In order to see that
mα is onto let v∈Zα. Due to
lemma A.9 for each u∈I there exists (θ(u),λ(u))∈Hu,⊥2,α(R2)⊕H2,α(R) such that*
[TABLE]
where θˉ(u)(ξ,Z)=θ(u)(ξ,γZ+ξ).
It holds for h∈H2(R2) the inequality
[TABLE]
Using (60) for h(ξ,x)=(1+∣ξ∣2+∣x∣2)2α[θ(u)−θ(uˉ)](ξ,x), (59) and lemma A.9 we obtain:
[TABLE]
This implies that
(θ,λ)∈Yα,
since v∈Zα.
The inverse mapping theorem yields that [mα]−1 is bounded, since mα is bounded.
□
In the following we introduce the operator L^uα that will be used for the main decomposition of this appendix in corollary A.13.
Definition A.11.
We define the following operators.
(a)
Luα:H2,α(R2)⊕H1,α(R2)→H1,α(R2)⊕L2,α(R2)* given by*
It follows from lemma A.7 that ker[Luα]∗={λJt1(u),λ∈H2,α(R)}. One proves first the case α=0 analogously to the proof of lemma A.5 which can be used to deduce the case α=0.
□
Corollary A.13 (direct sum).
[TABLE]
Proof 0.
"⊃": clear.
"⊂": Let (v,w)∈H1,α(R2)⊕L2,α(R2) then
there exists due to lemma A.12(θ,ψ)=(θ(u),ψ(u))∈H2,α(R2)⊕H1,α(R2)∩(kerLuα)⊥,L2,α⊕L2,α and λ=λ(u)∈H2,α(R) such that
[TABLE]
Assume without loss of generality
∣λ∣H2,α(R)=0, then
⟨λt1(u),λJt2(u)⟩L2,a(R2)⊕L2,a(R2)=0.
Thus
due to lemma A.12 there exist
(θˉ,ψˉ)=(θˉ(u),ψˉ(u))∈H2,α(R2)⊕H1,α(R2)∩(kerLuα)⊥,L2,α⊕L2,α and
0=λˉ=λˉ(u)∈H2,α(R) such that
[TABLE]
This is an identity in H1,α(R2)⊕L2,α(R2).
Fixing ξ and pairing this identity with Jt1(ξ,u,⋅) in
Lx2(R)⊕Lx2(R) yields due to lemma A.5 for a.e. ξ∈R
the identity
λ(ξ,u)=η(u)λˉ(ξ,u),
where
η(u):=γ(u)−3m−1(u2γ3∣θK′′∣L2(R)2+γ∣θK′∣L2(R)2)∈R.
Thus using (61) we obtain
Let α,n∈N. We want to show that the linear operator Mnα:Ynα(u∗)→Znα(u∗)
is invertible if u∗ is small. The operator Mnα contains derivatives with respect to ξ and x which makes it difficult to analyze it. Therefore we consider first an operator M~α:Yα→Zα,
which contains only derivatives with respect to x. Using corollary A.13 we prove that Mnα is invertible.
Definition B.1.
We define the following operators.
(a)
Muα:H2,α(R2)⊕H1,α(R2)⊕H2,α(R)→H1,α(R2)⊕L2,α(R2)* given by*
Zα=Zα(u∗)* is the space \bigg{\{}z=(v,w)\in C(I(u_{*}),\bar{Z}^{\alpha}):\|z\|_{Z^{\alpha}(u_{*})}<\infty\bigg{\}}\,
with the finite norm*
[TABLE]
Lemma B.2.
The linear operator
M~α:Yα(u∗)→Zα(u∗),(θ,ψ,λ)↦M~α(θ,ψ,λ),
given by
[TABLE]
is invertible.
Proof 0.
M~α* is onto: Let (v,w)∈Zα. Due to corollary A.13 for all u∈I there exist
(\theta(u),\psi(u))\in[H^{2,\alpha}(\mathbb{R}^{2})\oplus H^{1,\alpha}(\mathbb{R}^{2})\Big{]}\cap\left({\rm ker}\,{\cal L}_{u}^{\alpha}\right)^{\perp,L^{2,\alpha}\oplus L^{2,\alpha}}
and λ(u)∈H2,α(R) such that*
[TABLE]
This is an identity in H1,α(R2)⊕L2,α(R2).
By fixing ξ, using lemma A.5 and pairing (62) with Jtξ(ξ,u,x)
in Lx2,α(R)⊕Lx2,α(R) we obtain
λ∈C(I,L2,α(R)).
Further (62) yields
[TABLE]
where Z=γ(x−ξ).
Hence (((1+u2)−1γ−3θ)kerLuα,λ)∈Yα due to lemma A.10, since (1+u2)−1γ−3(w−u∂xv)∈Zα and λ∈C(I,L2,α(R)).
Thus
(θ,ψ,λ)∈Yα and M~α(θ,ψ,λ)=(v,w).
M~α is one-to-one due to corollary A.13. The inverse mapping theorem yields the claim.
□
Next, we want to show that the operator norm of
M^u−1
is bounded by a function, which is continuous in u. We start with a preparing lemma.
Lemma B.3 (Norm of [M^uα]−1).
There exists a constant cα>0 such that
[TABLE]
Proof 0.
Let
∣v∣L2,α(R2)≤1.
Due to lemma A.9 and lemma A.7 there exists
(θ,λ)∈Hu,⊥2,α(R2)⊕H2,α(R),
such that
[TABLE]
where θˉ(u)(ξ,Z)=θ(u)(ξ,γZ+ξ).
Using (60)
for h(ξ,x)=(1+∣ξ∣2+∣x∣2)2αθ(ξ,x)
we obtain
[TABLE]
□
Lemma B.4 (Norm of [M^uα]−1).
There exists a continuous function Cα:(−1,1)→R such that
[TABLE]
Proof 0.
Let
∣(v,w)∣H1,α(R2)⊕L2,α(R2)≤1.
Due to corollary A.13 there exists
(\theta,\psi,\lambda)\in\Big{[}H^{2,\alpha}(\mathbb{R}^{2})\oplus H^{1,\alpha}(\mathbb{R}^{2})\Big{]}\cap{\rm ker}\,{\cal L}_{u}^{\perp,L^{2,\alpha}\oplus L^{2,\alpha}}\oplus H^{2,\alpha}(\mathbb{R})\,
such that
(vw)=M^uα(θ,ψ,λ).
lemma B.3 and the expression for
(θ,ψ,λ) from the proof of lemma B.2 imply the claim.
□
Definition B.5.
We define the following operators.
(a)
Kuα:H2,α(R2)⊕H1,α(R2)→H1,α(R2)⊕L2,α(R2)* given by *
[TABLE]
(b)
Nuα:H2,α(R2)⊕H1,α(R2)⊕H2,α(R)→H1,α(R2)⊕L2,α(R2)* given by*
Using von Neumann’s theorem we are able to prove now that for small u∗ an extension of the operator Mnα:Ynα(u∗)→Znα(u∗) is invertible. Before proceeding to the proof we introduce the following definition in order to specify u∗.
Definition B.6.
Let Cα be a specific fixed function from lemma B.4. Set
[TABLE]
Corollary B.7.
The linear operator
Mα:Yα(u∗)→Zα(u∗),(θ,ψ,λ)↦Mα(θ,ψ,λ),
given by
[TABLE]
is invertible if u∗<u~α.
Proof 0.
The operator M~α
is invertible by lemma B.2
and it holds for its operator norm that
[M~α]−1≤sup∣u∣≤u∗Cα(u).
Let Pα:Yα(u∗)→Zα(u∗) be given by
[TABLE]
It holds that
∥Pα∥≤sup∣u∣≤u∗∣u∣
and thus
[TABLE]
due to (63), since u∗<u~α.
Hence
Pα+M~α=Mα is invertible by von Neumann’s theorem.
□
Analogously to corollary B.7
one shows first that the corresponding operator on spaces of higher regularity in (ξ,x) is invertible.
The claim of proposition 3.2 for the operator on spaces of higher regularity in u and in (ξ,x) follows by using difference quotients, orthogonal projection and the inverse mapping theorem.
□
Acknowledgements
My sincere gratitude goes to my PhD advisor Markus Kunze for the continuous support, his patience and motivation.
I would like to express my deep appreciation to Justin Holmer for numerous helpful and fruitful discussions.
Bibliography32
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[Ben 76] T. Brooke Benjamin. Applications of Leray-Schauder degree theory to problems of hydrodynamic stability. Math. Proc. Cambridge Philos. Soc. , 79(2):373–392, 1976.
2[Bon 75] J. Bona. On the stability theory of solitary waves. Proc. Roy. Soc. London Ser. A , 344(1638):363–374, 1975.
3[BP 92] V. S. Buslaev and G. S. Perel’man. On nonlinear scattering of states which are close to a soliton. Astérisque , (210):6, 49–63, 1992. Méthodes semi-classiques, Vol. 2 (Nantes, 1991).
4[Bre 11] Haim Brezis. Functional analysis, Sobolev spaces and partial differential equations . Universitext. Springer, New York, 2011.
7[FK 39] Y. I. Frenkel, T. Kontorova. J. Phys. Acad. Sci. USSR 1, 137, 1939.
8[FGJS 04] J. Fröhlich, S. Gustafson, B. L. G. Jonsson, and I. M. Sigal. Solitary wave dynamics in an external potential. Comm. Math. Phys. , 250(3):613–642, 2004.