Vincent Duchêne

Dag Nilsson

Erik Wahlén

Solitary wave solutions to a class of modified Green-Naghdi systems

We provide the existence and asymptotic description of solitary wave solutions to a class of modified Green-Naghdi systems, modeling the propagation of long surface or internal waves. This class was recently proposed by Duchêne, Israwi and Talhouk [23] in order to improve the frequency dispersion of the original Green-Naghdi system while maintaining the same precision. The solitary waves are constructed from the solutions of a constrained minimization problem. The main difficulties stem from the fact that the functional at stake involves low order non-local operators, intertwining multiplications and convolutions through Fourier multipliers.

Introduction

Motivation

In this work, we study solitary traveling waves for a class of long-wave models for the propagation of surface and internal waves. Starting with the serendipitous discovery and experimental investigation by John Scott Russell, the study of solitary waves at the surface of a thin layer of water in a canal has a rich history [START_REF] Darrigol | The spirited horse, the engineer, and the mathematician: water waves in nineteenth-century hydrodynamics[END_REF]. In particular, it is well-known that the most widely used nonlinear and dispersive models for the propagation of surface gravity waves, such as the Korteweg-de Vries equation or the Boussinesq and Green-Naghdi systems, admit explicit families of solitary waves [START_REF] Boussinesq | Théorie des ondes et des remous qui se propagent le long d'un canal rectangulaire horizontal, en communiquant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond[END_REF][START_REF] Rayleigh | On waves[END_REF][START_REF] Korteweg | On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves[END_REF][START_REF] Serre | Contribution à l'étude des écoulements permanents et variables dans les canaux[END_REF][START_REF] Chen | Exact solutions of various Boussinesq systems[END_REF]. These equations can be derived as asymptotic models for the so-called water waves system, describing the motion of a two-dimensional layer of ideal, incompressible, homogeneous, irrotational fluid with a free surface and a flat impermeable bottom; we let the reader refer to [START_REF] Lannes | The water waves problem[END_REF] and references therein for a detailed account of the rigorous justification of these models. Among them, the Green-Naghdi model is the most precise, in the sense that it does not assume that the surface deformation is small. However, the validity of all these models relies on the hypothesis that the depth of the layer is thin compared with the horizontal wavelength of the flow and, as expected, the models do not describe the system accurately (for instance the dispersion relation of infinitesimally small waves) in a deep water situation. In order to tackle this issue, one of the authors has recently proposed in [START_REF] Duchêne | A new class of two-layer Green-Naghdi systems with improved frequency dispersion[END_REF] a new family of models:

     ∂ t ζ + ∂ x w = 0, ∂ t h -1 w + Q F [h](h -1 w) + g∂ x ζ + 1 2 ∂ x (h -1 w) 2 = ∂ x R F [h, h -1 w] , (1.1)
where

Q F [h]u def = - 1 3 h -1 ∂ x F h 3 ∂ x F{u} , R F [h, u] def = 1 3 uh -1 ∂ x F h 3 ∂ x F{u} + 1 2 h∂ x F{u} 2 .
Here, ζ is the surface deformation, h = d + ζ the total depth (where d is the depth of the layer at rest), u the layer-averaged horizontal velocity, w = hu the horizontal momentum and g the gravitational acceleration; see Figure 1. Finally, F def = F(D) is a Fourier multiplier, i.e.

F{ϕ}(k) = F(k) ϕ(k).
The original Green-Naghdi model is recovered when setting F(k) ≡ 1. Any other choice satisfying F(k) = 1 + O(k 2 ) enjoys the same precision (in the sense of consistency) in the shallow-water regime and the specific choice of F(k) = 3 d|k| tanh(d|k|) -3 d 2 |k| 2 allows to obtain a model whose linearization around constant states fits exactly with the one of the water waves system. Compared with other strategies such as the Benjamin-Bona-Mahony trick [START_REF] Benjamin | Model equations for long waves in nonlinear dispersive systems[END_REF][START_REF] Madsen | A new form of the Boussinesq equations with improved linear dispersion characteristics[END_REF] or using different choices of velocity unknowns [START_REF] Bona | A model for the two-way propagation of water waves in a channel[END_REF][START_REF] Nwogu | Alternative form of boussinesq equations for nearshore wave propagation[END_REF], an additional advantage of (1.1) is that it preserves the Hamiltonian structure of the model, which turns out to play a key role since the existence of solitary waves will be deduced from a variational principle. The study of [START_REF] Duchêne | A new class of two-layer Green-Naghdi systems with improved frequency dispersion[END_REF] is not restricted to surface propagation, but is rather dedicated to the propagation of internal waves at the interface between two immiscible fluids, confined above and below by rigid, impermeable and flat boundaries. Such a configuration appears naturally as a model for the ocean, as salinity and temperature may induce sharp density stratification, so that internal solitary waves are observed in many places [START_REF] Ostrovsky | Do internal solitons exist in the ocean?[END_REF][START_REF] Jackson | An atlas of internal solitary-like waves and their properties[END_REF][START_REF] Helfrich | Long nonlinear internal waves[END_REF]. Due to the weak density contrast, the observed solitary waves typically have much larger amplitude than their surface counterpart, hence the bilayer extension of the Green-Naghdi system introduced by [START_REF] Miyata | Long internal waves of large amplitude[END_REF][START_REF] Tseva | Unsteady long waves in a two-layer fluid[END_REF][START_REF] Choi | Fully nonlinear internal waves in a two-fluid system[END_REF], often called Miyata-Choi-Camassa model, is a very natural choice. It however suffers from strong Kelvin-Helmholtz instabilities -in fact stronger than the ones of the water waves system for large frequencies -and the work in [START_REF] Duchêne | A new class of two-layer Green-Naghdi systems with improved frequency dispersion[END_REF] was motivated by taming these instabilities. The modified bilayer system reads

       ∂ t ζ + ∂ x w = 0, ∂ t h 1 + γh 2 h 1 h 2 w + Q F γ,δ [ζ]w + (γ + δ)∂ x ζ + 1 2 ∂ x h 2 1 -γh 2 2 (h 1 h 2 ) 2 w 2 = ∂ x R F γ,δ [ζ, w] (1.2)
where we denote

h 1 = 1 -ζ, h 2 = δ -1 + ζ, Q F γ,δ [ζ]w def = Q F 2 [h 2 ](h -1 2 w) + γQ F 1 [h 1 ](h -1 1 w) and R F γ,δ [ζ, w] def = R F 2 [h 2 , h -1 2 w] -γR F 1 [h 1 , h -1 1 w],
with

Q F i [h i ]u i def = - 1 3 h -1 i ∂ x F i h 3 i ∂ x F i {u i } , R F i [h i , u i ] def = 1 3 u i h -1 i ∂ x F i h 3 i ∂ x F i {u i } + 1 2 h i ∂ x F i {u i } 2 .
layer and w = h 1 h 2 (u 2 -γu 1 )/(h 1 + γh 2 ) is the shear momentum. In this formulation we have used dimensionless variables, so that the depth at rest of the upper layer is scaled to 1, whereas the one of the lower layer is δ -1 , in which δ is the ratio of the depth at rest of the upper layer to the depth at rest of the lower layer (see Figure 1). Similarly, γ is the ratio of the upper layer over the lower layer densities. As a consequence of our scaling, the celerity of infinitesimally small and long waves is c 0 = 1. Once again, F i (i = 1, 2) are Fourier multipliers. The choice F id i (k) ≡ 1 yields the Miyata-Choi-Camassa model while

F imp i (k) = 3 δ -1 i |k| tanh(δ -1 i |k|) - 3 δ -2 i |k| 2 ,
with convention δ 1 = 1, δ 2 = δ, fits the behavior of the full bilayer Euler system around constant states, and thus gives hope for an improved precision when weak nonlinearities are involved. Note that compared to equations ( 7)-( 9) in [START_REF] Duchêne | A new class of two-layer Green-Naghdi systems with improved frequency dispersion[END_REF] we have scaled the variables so that the shallowness parameter µ and amplitude parameter do not appear in the equations. This is for notational convenience since the parameters do not play a direct role in our results. On the other hand, we only expect the above model to be relevant for describing water waves in the regime µ 1 and the solutions that we construct in the end are found in the long-wave regime , µ 1. In the following, we study solitary waves for the bilayer system (1.2), noting that setting γ = 0 immediately yields the corresponding result for the one-layer situation, namely system (1.1). Our results are valid for a large class of parameters γ, δ and Fourier multipliers F 1 , F 2 , described hereafter. Our results are twofold: i. We prove the existence of a family of solitary wave solutions for system (1.2);

ii. We provide an asymptotic description for this family in the long-wave regime.

These solitary waves are constructed from the Euler-Lagrange equation associated with a constrained minimization problem, as made possible by the Hamiltonian structure of system (1.2). There are however several difficulties compared with standard works in the literature following a similar strategy (see e.g. [START_REF] Pava | Nonlinear dispersive equations. Existence and stability of solitary and periodic travelling wave solutions[END_REF] and references therein). Our functional cannot be written as the sum of the linear dispersive contribution and the nonlinear pointwise contribution: Fourier multipliers and nonlinearities are entangled. What is more, the operators involved are typically of low order (F is a smoothing operator). In order to deal with this situation, we follow a strategy based on penalization and concentration-compactness used in a number of recent papers on the water waves problem (see e.g. [START_REF] Buffoni | Existence and conditional energetic stability of capillary-gravity solitary water waves by minimisation[END_REF][START_REF] Buffoni | Existence and conditional energetic stability of three-dimensional fully localised solitary gravity-capillary water waves[END_REF][START_REF] Groves | Existence and conditional energetic stability of solitary gravity-capillary water waves with constant vorticity[END_REF] and references therein) and in particular, in a recent work by one of the authors on nonlocal model equations with weak dispersion, [START_REF] Ehrnström | On the existence and stability of solitarywave solutions to a class of evolution equations of Whitham type[END_REF]. Thus we show that the strategy therein may be favorably applied to bidirectional systems of equations in addition to unidirectional scalar equations such as the Whitham equation. Roughly speaking, the strategy is the following. The minimization problem is first solved in periodic domains using a penalization argument do deal with the fact that the energy functional is not coercive. This allows to construct a special minimizing sequence for the real line problem by letting the period tend to infinity, which is essential to rule-out the dichotomy scenario in Lions' concentration-compactness principle. The long-wave description follows from precise asymptotic estimates and standard properties of the limiting (Korteweg-de Vries) model. When the Fourier multipliers F i have sufficiently high order, we can in fact avoid the penalization argument and consider the minimization problem on the real line directly, since any minimizing sequence is then also a special minimizing sequence. In particular, this is the case for the original Miyata-Choi-Camassa model (and of course also the Green-Naghdi system).

Our existence proof unfortunately gives no information about stability, since our variational formulation does not involve conserved functionals; see the discussion in Section 1.2. If sufficiently strong surface tension is included in the model, we expect that a different variational formulation could be used which also yields a conditional stability result (see [START_REF] Buffoni | Existence and conditional energetic stability of capillary-gravity solitary water waves by minimisation[END_REF][START_REF] Buffoni | Existence and conditional energetic stability of three-dimensional fully localised solitary gravity-capillary water waves[END_REF][START_REF] Groves | Existence and conditional energetic stability of solitary gravity-capillary water waves with constant vorticity[END_REF]). A similar situation appears e.g. in the study of Boussinesq systems [START_REF] Chen | Solitary-wave solutions to Boussinesq systems with large surface tension[END_REF][START_REF]Existence of traveling-wave solutions to Boussinesq systems[END_REF].

The minimization problem

We now set up the minimization problem which allows to obtain solitary waves of system (1.2). We seek traveling waves of (1.2), namely solutions of the form (abusing notation) ζ(t, x) = ζ(x -ct), w(t, x) = w(x -ct); from which we deduce

-c∂ x ζ + ∂ x w = 0 ; -c∂ x A F γ,δ [ζ]w + (γ + δ)∂ x ζ + 1 2 ∂ x h 2 1 -γh 2 2 h 2 1 h 2 2 w 2 -∂ x R F γ,δ [ζ, w] = 0,
where we denote

A F γ,δ [ζ]w def = A F 2 [h 2 ](h -1 2 w) + γA F 1 [h 1 ](h -1 1 w), A F i [h i ]u i def = u i + Q F i [h i ]u i .
Integrating these equations and using the assumption (since we restrict ourselves to solitary waves)

lim |x|→∞ ζ(x) = lim |x|→∞ w(x) = 0 yields the system of equations        -cζ + w = 0, -cA F γ,δ [ζ]w + (γ + δ)ζ + 1 2 h 2 1 -γh 2 2 (h 1 h 2 ) 2 w 2 -R F γ,δ [ζ, w] = 0. (1.3)
We now observe that system (1.2) enjoys a Hamiltonian structure. Indeed, define the functional

H(ζ, w) def = 1 2 ∞ -∞ (γ + δ)ζ 2 + wA F γ,δ [ζ]w dx.
Under reasonable assumptions on F 1 , F 2 , and for sufficiently regular ζ, A F γ,δ [ζ] defines a well-defined, symmetric, positive definite operator [START_REF] Duchêne | A new class of two-layer Green-Naghdi systems with improved frequency dispersion[END_REF]. We may thus introduce the variable

v def = A F γ,δ [ζ]w, (1.4)
and write

H(ζ, v) def = 1 2 ∞ -∞ (γ + δ)ζ 2 + v(A F γ,δ [ζ]) -1 v dx.
It is now straightforward to check that (1.2) can be written in terms of functional derivatives of H:

∂ t ζ = -∂ x (d v H) ; ∂ t v = -∂ x (d ζ H) . (1.5)
It is therefore tempting to seek traveling waves through the time-independent quantities

H(ζ, v) = 1 2 ∞ -∞ (γ + δ)ζ 2 + vA F γ,δ [ζ] -1 v dx and I(ζ, v) def = ∞ -∞ ζv dx.
Indeed any solution to (1.3) is a critical point of the functional H -cI:

dH(ζ, v) -cdI(ζ, v) = 0 ⇐⇒ d v H -cζ = 0 and d ζ H -cv = 0,
which, by (1.5), is the desired system of equations. However, notice that from Weyl's essential spectrum theorem, one has

spec ess d 2 H(ζ, v) -cd 2 I(ζ, v) = spec ess d 2 H ∞ (ζ, v) -cd 2 I ∞ (ζ, v)
where H ∞ , I ∞ are the asymptotic operators as |x| → ∞, i.e.

d 2 H ∞ (ζ, v)-cd 2 I ∞ (ζ, v) = d 2 H(0, 0)-cd 2 I(0, 0) = γ + δ -c -c (γ + δ -γ 3 (∂ x F 1 ) 2 -1 3δ (∂ x F 2 ) 2 ) -1 .
Since the spectrum of the above operator has both negative and positive components, the desired critical point is neither a minimizer nor a maximizer, as noticed (for the Green-Naghdi system) in [START_REF] Li | Hamiltonian structure and linear stability of solitary waves of the Green-Naghdi equations[END_REF].

We will obtain solutions to (1.3) from a constrained minimization problem depending solely on the variable ζ. Notice that for each fixed c and ζ,

the functional v → H(ζ, v) -cI(ζ, v) has a unique critical point, v c,ζ = cA F γ,δ [ζ]ζ . Substituting v c,ζ into H(ζ, v) -cI(ζ, v), we obtain H(ζ, v c,ζ ) -cI(ζ, v c,ζ ) = 1 2 ∞ -∞ (γ + δ)ζ 2 -c 2 ζA F γ,δ [ζ]ζ dx = γ + δ 2 ζ 2 L 2 - c 2 2 I(ζ, A F γ,δ [ζ]ζ). Observe now that (ζ, v) is a critical point of H(ζ, v) -cI(ζ, v) if and only if ζ is a critical point of H(ζ, v c,ζ ) -cI(ζ, v c,ζ ) and v = v c,ζ . We thus define E(ζ) def = I(ζ, A F γ,δ [ζ]ζ) = ∞ -∞ ζA F γ,δ [ζ]ζ dx = ∞ -∞ h 1 + γh 2 h 1 h 2 ζ 2 + γ 3 (1 -ζ) 3 ∂ x F 1 { ζ 1 -ζ } 2 + 1 3 (δ -1 + ζ) 3 ∂ x F 2 { ζ δ -1 + ζ } 2 dx = γE(ζ) + E(ζ), (1.6)
where

E(ζ) = ∞ -∞ ζ 2 1 -ζ + 1 3 (1 -ζ) 3 ∂ x F 1 { ζ 1 -ζ } 2 dx, E(ζ) = ∞ -∞ ζ 2 δ -1 + ζ + 1 3 (δ -1 + ζ) 3 ∂ x F 2 { ζ δ -1 + ζ } 2 dx
and look for critical points of

H(ζ, v c,ζ ) -cI(ζ, v c,ζ
) by considering the minimization problem

arg min E(ζ), (γ + δ) ζ 2 L 2 = q , (1.7)
with c -2 acting as a Lagrange multiplier. Another way of thinking of this reduction is to solve the first equation in (1.3) for w and substitute the solution w = cζ into the second equation, yielding

c 2 A F γ,δ [ζ]ζ = (γ + δ)ζ + c 2 2 h 2 1 -γh 2 2 h 2 1 h 2 2 ζ 2 -c 2 R F γ,δ [ζ, ζ]. (1.8) Computing dE(ζ) = 2 h 1 + γh 2 h 1 h 2 ζ- h 2 1 -γh 2 2 h 2 1 h 2 2 ζ 2 - 2 3 δ -1 h -2 2 ∂ x F 2 h 3 2 ∂ x F 2 {h -1 2 ζ} - 2γ 3 h -2 1 ∂ x F 1 h 3 1 ∂ x F 1 {h -1 1 ζ} + h 2 ∂ x F 2 {h -1 2 ζ} 2 -γ h 1 ∂ x F 1 {h -1 1 ζ} 2 , (1.9)
we find that

1 2 dE(ζ) = A F γ,δ [ζ]ζ - 1 2 h 2 1 -γh 2 2 h 2 1 h 2 2 ζ 2 + R F γ,δ [ζ, ζ],
which allows us to recognize (1.8) as the Euler-Lagrange equation for (1.7).

Statement of the results

For the sake of readability, we postpone to Section 2 the definition and (standard) notations of the functional spaces used herein. The class of Fourier multipliers for which our main result is valid is the following.

Definition 1.1 (Admissible class of Fourier multipliers).

i.

F(k) = F(|k|) and 0 < F ≤ 1; ii. F ∈ C 2 (R), F(0) = 1 and F (0) = 0;
iii. There exists an integer j ≥ 2 such that

∂ j k (kF(k)) ∈ L 2 (R);
iv. There exists θ ∈ [0, 1) and C F ± > 0 such that

C F -(1 + |k|) -θ ≤ F(k) ≤ C F + (1 + |k|) -θ .
We also introduce a second class of strongly admissible Fourier multipliers which is used in our second result. Assumption 1.4 (Admissible parameters). In the following, we fix γ ≥ 0, δ ∈ (0, ∞) such that δ 2 -γ = 0. We also fix a positive number ν such that ν ≥ 1 -θ and ν > 1/2 (the second condition is automatically satisfied if θ < 1/2). Finally, fix R an arbitrary positive constant.

Remark 1.5. Our results hold for any values of the parameters (γ, δ) ∈ [0, ∞) × (0, ∞) such that δ 2 = γ, although admissible values for q 0 depend on the choice of the parameters. However, not all parameters are physically relevant in the oceanographic context. When γ > 1, the upper fluid is heavier than the lower fluid, and the system suffers from strong Rayleigh-Taylor instabilities [START_REF] Chandrasekhar | Hydrodynamic and hydromagnetic stability[END_REF]. In the bilayer setting, the use of the rigid-lid assumption is well-grounded only when the density contrast, 1 -γ, is small. In this situation, one may use the Boussinesq approximation, that is set γ = 1; see [START_REF] Duchêne | On the rigid-lid approximation for two shallow layers of immiscible fluids with small density contrast[END_REF] in the dispersionless setting. Notice however that system (1.2) exhibits unstable modes that are reminiscent of Kelvin-Helmholtz instabilities when the Fourier multipliers F i satisfy Definition 1.1,iv with θ ∈ [0, 1); see [START_REF] Duchêne | A new class of two-layer Green-Naghdi systems with improved frequency dispersion[END_REF]. It is therefore noteworthy that internal solitary waves in the ocean and in laboratory experiments are remarkably stable and fit very well with the Miyata-Choi-Camassa predictions [START_REF] Helfrich | Long nonlinear internal waves[END_REF]. The sign of the parameter δ 2 -γ is known to determine whether long solitary waves are of elevation or depression type, as corroborated by Theorem 1.7. At the critical value δ 2 = γ, the first-order model would be the modified (cubic) KdV equation, predicting that no solitary wave exists [START_REF] Djordjevic | The fission and disintegration of internal solitary waves moving over two-dimensional topography[END_REF].

We study the constrained minimization problem arg min

ζ∈V q,R E(ζ), (1.10) with V q,R = {ζ ∈ H ν (R) : ζ H ν (R) < R, (γ + δ) ζ 2 L 2 (R) = q},
and q ∈ (0, q 0 ), with q 0 sufficiently small. Notice in particular that as soon as q is sufficiently small ζ L ∞ < min(1, δ -1 ) (by Lemma 2.1 thereafter and since ν > 1/2) and E(ζ) is well-defined (by Lemmata 2.3 and 2.4 and since ν ≥ 1 -θ) for any ζ ∈ V q,R . Any solution will satisfy the Euler-Lagrange equation

dE(ζ) + 2α(γ + δ)ζ = 0, (1.11)
where α is a Lagrange multiplier. Equation (1.11) is exactly (1.8) with (-α) -1 = c 2 , and therefore provides a traveling-wave solution to (1.2).

Our goal is to prove the following theorems.

Theorem 1.6. Let γ, δ, ν, R satisfying Assumption 1.4 and F i , i = 1, 2 be admissible in the sense of Definition 1.1. Let D q,R be the set of minimizers of E over V q,R . Then there exists q 0 > 0 such that for all q ∈ (0, q 0 ), the following statements hold:

• The set D q,R is nonempty and each element in D q,R solves the traveling wave equation (1.8), with c 2 = (-α) -1 > 1. Thus for any ζ ∈ D q,R , ζ(x ± ct), w ± = ±cζ(x ± ct) is a supercritical solitary wave solution to (1.2).

• For any minimizing sequence

{ζ n } n∈N for E in V q,R such that sup n∈N ζ n H ν (R) < R, there exists a sequence {x n } n∈N of real numbers such that a subsequence of {ζ n (•+x n )} n∈N converges (strongly in H ν (R) if ν = 1 -θ > 1/2; weakly in H ν (R)
and strongly in H s (R) for s ∈ [0, ν) otherwise) to an element in D q,R .

• There exist two constants m, M > 0 such that

ζ 2 H ν (R) ≤ M q and c -2 = -α ≤ 1 -mq 2 3 ,
uniformly over q ∈ (0, q 0 ) and ζ ∈ D q,R .

Theorem 1.7. In addition to the hypotheses of Theorem 1.6, assume that F i , i = 1, 2, are strongly admissible in the sense of Definition 1.2. Then there exists q 0 > 0 such that for any q ∈ (0, q 0 ), each ζ ∈ D q,R belongs to H s (R) for any s ≥ 0 and

sup ζ∈D q,R inf x0∈R q -2 3 ζ(q -1/3 •) -ξ KdV (• -x 0 ) H 1 (R) = O(q 1 6 )
where

ξ KdV (x) = α 0 (γ + δ) δ 2 -γ sech 2 1 2 3α 0 (γ + δ) γ + δ -1 x
is the unique (up to translation) solution of the KdV equation (5.2) and

α 0 = 3 4 (δ 2 -γ) 4 (γ + δ) 4 (γ + δ -1 ) 1 3
.

In addition, the number α, defined in Theorem 1.6, satisfies

α + 1 = q 2 3 α 0 + O(q 5 6 ),
uniformly over q ∈ (0, q 0 ) and ζ ∈ D q,R .

Technical results

In the following, we denote C(λ 1 , λ 2 , . . . ) a positive constant depending non-decreasingly on the parameters λ 1 , λ 2 , . . . . We write A B when A ≤ CB with C a nonnegative constant whose value is of no importance. We do not display the dependence with respect to the parameters γ, δ, C Fi ± and regularity indexes.

Functional setting on the real line Here and thereafter, we denote L 2 (R) the standard Lebesgue space of square-integrable functions, endowed with the norm

f L 2 = ( ∞ -∞ |f (x)| 2 dx) 1/2 .
The real inner product of

f 1 , f 2 ∈ L 2 (R) is denoted by f 1 , f 2 = R f 1 (x)f 2 (x)dx.
We use the same notation for duality pairings which are clear from the context. The space L ∞ (R) consists of all essentially bounded, Lebesgue-measurable functions f , endowed with the norm f L ∞ = ess sup x∈R |f (x)|. For any real constant s ∈ R, H s (R) denotes the Sobolev space of all tempered distributions f with finite norm f

H s = Λ s f L 2 < ∞, where Λ is the pseudo-differential operator Λ = (1 -∂ 2 x ) 1 2 . For n ∈ N, C n (R)
is the space of functions having continuous derivatives up to order n, and C ∞ (R) = n∈N C n (R). The Schwartz space is denoted S(R) and the tempered distributions S (R). We use the following convention for the Fourier transform:

F f (k) = f (k) def = 1 √ 2π R f (x)e -ixk dx.
We start with standard estimates in Sobolev spaces. The following interpolation estimates are standard and used without reference in our proofs.

Lemma 2.1 (Interpolation estimates). Let f ∈ H µ (R), with µ > 1/2. i. One has f ∈ L ∞ (R) and f L ∞ f 1-1 2µ L 2 f 1 2µ H µ .
ii. For any δ ∈ (0, µ), one has f ∈ H µ-δ (R) and

f H µ-δ ≤ f δ µ L 2 f 1-δ µ H µ .
The following Lemma is given for instance in [START_REF] Benzoni-Gavage | Multidimensional hyperbolic partial differential equations. First-order systems and applications[END_REF]Theorem C.12].

Lemma 2.2 (Composition estimate). Let G be a smooth function vanishing at 0, and

f ∈ H µ (R) with µ > 1/2. Then G • f ∈ H µ (R) and we have G • f H µ ≤ C( f L ∞ ) f H µ . Lemma 2.3 (Product estimates). i. For any f, g ∈ L ∞ (R) ∩ H s (R) with s ≥ 0, one has f g ∈ H s (R) and f g H s f H s g L ∞ + g H s f L ∞ .
ii. For any f ∈ H s (R), g ∈ H t (R) with s + t ≥ 0, and let r such that min(s, t) ≥ r and r < s + t -1/2. Then one has f g ∈ H r (R) and

f g H r f H s g H t .
iii. For any

ζ ∈ L ∞ (R) such that ζ L ∞ ≤ 1 -h 0 with h 0 > 0 and any f ∈ L ∞ (R), one has f 1 + ζ L ∞ ≤ C(h -1 0 ) f L ∞ . iv. For any ζ ∈ H µ (R) with µ > 1/2 such that ζ L ∞ ≤ 1 -h 0 with h 0 > 0 and any f ∈ H s (R) with s ∈ [-µ, µ], one has f 1 + ζ H s ≤ C(h -1 0 , ζ H µ ) f H s .
Proof.

(R), f g H s g H µ f H s . Hence f 1 + ζ H s ≤ f H s + f ζ 1 + ζ H s ≤ f H s + ζ 1 + ζ H µ f H s We conclude by ζ 1 + ζ H µ C(h -1 0 ) ζ H µ ,
where we have used Lemma 2.2, and the estimate is proved.

The following Lemma justifies the assumptions of admissible Fourier multipliers in Definition 1.1.

Lemma 2.4 (Properties of admissible Fourier multipliers). Any admissible Fourier multipler (in the sense of Definition 1.1), F i , satisfies the following.

i. The linear operator

∂ x F i (D) is bounded from H s (R) to H s-1+θ (R), for any s ∈ R, and
∂ x F i H s →H s-1+θ C Fi + .
Moreover, for any ζ ∈ H s+1-θ , one has

ζ 2 H s + (C Fi + ) -2 ∂ x F i {ζ} 2 H s ζ 2 H s+1-θ ζ 2 H s + (C Fi -) -2 ∂ x F i {ζ} 2 H s . ii. Let ϕ ∈ C ∞ (R) with compact support and [∂ x F i , ϕ]ζ = ∂ x F i {ϕζ} -ϕ∂ x F i {ζ}. Then [∂ x F i , ϕ]ζ L 2 ϕ L 1 ζ H 1-θ .
iii. There exists j ≥ 2 and C j such that for any ζ ∈ L 2 (R) with compact support

|∂ x F i {ζ}|(x) ≤ C j dist(x, supp(ζ)) j ζ L 2 , for a.a. x ∈ R \ supp(ζ).
Proof. The first result is obvious from Definition 1.1,i and the definition of Sobolev spaces. For the second, we shall first prove that the function

G i : k → kF i (k) satisfies |G i (k)| k 1-θ . (2.1)
To this aim, let us first consider G ∈ S(R) and χ a smooth cut-off function, such that χ(k) = 1 for |k| ≤ 1/2 and χ(k) = 0 for |k| ≥ 1. We decompose

|G |(k) ≤ |χ(D)G |(k) + |(1 -χ(D))G |(k).
For the first contribution, one has

|χ(D)G |(k) = 1 √ 2π R χ(ξ)G (k + ξ)dξ = 1 √ 2π R ( χ) (ξ)G(k + ξ)dξ sup ξ∈R |G(k + ξ)| k + ξ 1-θ k 1-θ • 1-θ χ L 1 ,
and the second contribution satisfies for any j ≥ 2,

|(1 -χ(D))G |(k) (1 -χ(ξ))|ξ| Ĝ(ξ) L 1 ξ -(j-1) |ξ| j Ĝ(ξ) L 1 G (j) L 2
, by the Cauchy-Schwarz inequality and Parseval's theorem. Thus we find, for any j ≥ 2,

|G |(k) • θ-1 G L ∞ k 1-θ + G (j) L 2 .
The same estimate applies to G(k) = kF i (k) by smooth approximation, and (2.1) follows from Definition 1.1. Now, we note that the Fourier transform of

[∂ x F i , ϕ]ζ is given by 1 √ 2π R (ikF i (k) -isF i (s)) φ(k -s) ζ(s)ds.
By (2.1) and the mean value theorem, we find that

|kF i (k) -sF i (s)| (1 + |s|) 1-θ |k -s|.
Due to Young's inequality and Parseval's theorem, we find

[∂ x F i , ϕ]ζ L 2 ϕ L 1 (1 + | • |) 1-θ ζ L 2 ϕ L 1 ζ H 1-θ .
For the third result, let us assume at first that the kernel

K i def = F -1 (ikF i (k)) ∈ L 2 (R). Then one has |∂ x F i {ζ}|(x) = 1 √ 2π R K i (x -y)ζ(y) dy = 1 √ 2π supp(ζ) (x -y) j K i (x -y)ζ(y) (x -y) j dy ≤ |K i,j | * |ζ| (x) √ 2π dist(x, supp(ζ)) j ζ L 2 dist(x, supp(ζ)) j ,
where we denote K i,j (x) = x j K i (x), remark that K i,j ∈ L 2 (R) by Definition 1.1,iii and Plancherel's theorem, and apply the Cauchy-Schwarz inequality to the convolution. If K i / ∈ L 2 (R), we obtain the result by regularizing K i (i.e. smoothly truncating F i ) and passing to the limit. Lemma 2.5. Let γ ≥ 0, δ > 0, µ > 1/2 and F i be admissible Fourier multipliers. Assume that

ζ ∈ H µ (R) is such that 1 -ζ L ∞ ≥ h 0 , δ -1 -ζ L ∞ ≥ h 0 , with h 0 > 0. Then there exist a constant C 0 = C(h -1 0 , ζ H µ ) such that C -1 0 ζ 2 H 1-θ ≤ E(ζ) ≤ C 0 ζ 2 H 1-θ .
Proof. We first deal with the contribution of E(ζ) defined in (1.6). By Lemma 2.4,i we get that

E(ζ) ≤ C( ζ L ∞ ) ζ 1 -ζ 2 H 1-θ and ζ 1 -ζ 2 H 1-θ ≤ C(h -1 0 )E(ζ).
By Lemma 2.3,iv, one has

ζ 1 -ζ H 1-θ ≤ C(h -1 0 , ζ H µ ) ζ H 1-θ ,
and the triangle inequality together with Lemma 2.3,ii yields

ζ H 1-θ ζ 1 -ζ H 1-θ + ζ 2 1 -ζ H 1-θ ζ 1 -ζ H 1-θ + ζ H µ ζ 1 -ζ H 1-θ .
Collecting the above information, we find that

C -1 0 ζ 2 H 1-θ ≤ E(ζ) ≤ C 0 ζ 2 H 1-θ , with C 0 = C(h -1 0 , ζ H µ )
. Similar estimates hold for E(ζ), and thus for

E(ζ) = γE(ζ) + E(ζ).
Lemma 2.6. Let γ ≥ 0, δ > 0, µ > 1/2 and F i be admissible Fourier multipliers. Assume that, for j ∈ {1, 2},

ζ j ∈ H µ (R) is such that 1 -ζ j L ∞ ≥ h 0 and δ -1 -ζ j L ∞ ≥ h 0 , with h 0 > 0. Then one has E(ζ 1 ) -E(ζ 2 ) ≤ C(h -1 0 , ζ 1 H µ , ζ 2 H µ ) ζ 1 -ζ 2 H µ .
Proof. As previously, we detail the result for E(ζ), as the similar estimate for E(ζ) is obtained in the same way. One has

E(ζ 1 ) -E(ζ 2 ) = R ζ 2 1 1 -ζ 1 - ζ 2 2 1 -ζ 2 + 1 3 (1 -ζ 1 ) 3 -(1 -ζ 2 ) 3 ∂ x F 1 { ζ 1 1 -ζ 1 } 2 + 1 3 (1 -ζ 2 ) 3 ∂ x F 1 { ζ 1 1 -ζ 1 } 2 -∂ x F 1 { ζ 2 1 -ζ 2 } 2 dx,
By Lemma 2.3, iii, and the Cauchy-Schwarz inequality, we immediately have

R ζ 2 1 1 -ζ 1 - ζ 2 2 1 -ζ 2 dx ≤ C(h -1 0 , ζ 1 L ∞ , ζ 2 L ∞ )( ζ 1 L 2 + ζ 2 L 2 ) ζ 1 -ζ 2 L 2 .
Similarly, we find by Lemma 2.4,i

R (1 -ζ 1 ) 3 -(1 -ζ 2 ) 3 ∂ x F 1 { ζ 1 1 -ζ 1 } 2 dx ≤ C( ζ 1 L ∞ , ζ 2 L ∞ ) ζ 1 -ζ 2 L ∞ ζ 1 1 -ζ 1 2 H 1-θ ,
and by Lemma 2.3,iv,

ζ 1 1 -ζ 1 2 H 1-θ ≤ C(h -1 0 , ζ 1 H µ ). Finally, R (1 -ζ 2 ) 3 ∂ x F 1 { ζ 1 1 -ζ 1 } 2 -∂ x F 1 { ζ 2 1 -ζ 2 } 2 dx ≤ C( ζ 2 L ∞ ) ζ 1 1 -ζ 1 - ζ 2 1 -ζ 2 H 1-θ ζ 1 1 -ζ 1 + ζ 2 1 -ζ 2 H 1-θ ,
and we conclude by Lemma 2.3, iv

ζ 1 1 -ζ 1 - ζ 2 1 -ζ 2 H 1-θ ≤ C(h -1 0 , ζ 1 H µ , ζ 2 H µ ) ζ 1 -ζ 2 H 1-θ , and ζ 1 1 -ζ 1 + ζ 2 1 -ζ 2 H 1-θ ≤ C(h -1 0 , ζ 1 H µ , ζ 2 H µ ).
The result is proved.

Lemma 2.7. Let γ ≥ 0, δ > 0, and F i be admissible Fourier multipliers. Let l ∈ {1, 2, 3} and

ζ ∈ H l (R) such that 1 -ζ L ∞ ≥ h 0 and δ -1 -ζ L ∞ ≥ h 0 , with h 0 > 0. Then one can decompose E(ζ) = R (γ + δ)ζ 2 + (γ -δ 2 )ζ 3 + γ + δ -1 3 (∂ x ζ) 2 dx + E rem (ζ), and
dE(ζ), ζ = R 2(γ + δ)ζ 2 + 3(γ -δ 2 )ζ 3 + 2 γ + δ -1 3 (∂ x ζ) 2 dx + dE rem (ζ), ζ ,
where

|E rem | + | dE rem (ζ), ζ | ≤ C(h -1 0 , ζ H 1 ) ζ 2 L ∞ ζ 2 L 2 + ζ L ∞ ∂ x ζ 2 L 2 + ∂ l x ζ L 2 ∂ x ζ L 2 .
Proof. We consider E(ζ); the corresponding expansion for E(ζ) is obtained similarly. We write

E(ζ) = R ζ 2 + ζ 3 + 1 3 (∂ x ζ) 2 dx + E rem (ζ),
where

E rem (ζ) = R ζ 4 1 -ζ dx + 1 3 R (1 -ζ) 3 ∂ x { ζ 1 -ζ } 2 -(∂ x ζ) 2 dx + R (1 -ζ) 3 ∂ x F 1 { ζ 1 -ζ } 2 -∂ x { ζ 1 -ζ } 2 dx.
Note that by the Cauchy-Schwarz inequality

R ζ 4 1 -ζ dx ≤ ζ 2 L ∞ ζ 2 L 2 h 0 and R (1 -ζ) 3 ∂ x ζ 1 -ζ 2 -(∂ x ζ) 2 dx = R ζ(∂ x ζ) 2 1 -ζ dx ≤ ζ L ∞ ∂ x ζ 2 L 2 h 0 . Moreover R (1 -ζ) 3 ∂ x F 1 ζ 1 -ζ 2 -∂ x ζ 1 -ζ 2 dx ≤ R (1 + |ζ|) 3 (∂ x F 1 -∂ x ) ζ 1 -ζ (∂ x F 1 + ∂ x ) ζ 1 -ζ dx.
Applying the Cauchy-Schwarz inequality, Plancherel's theorem and the estimates

|F 1 (k) -1| |k| l-1 , |F 1 (k) + 1| 1,
(by Definition 1.1,i and ii), we deduce

R (1 -ζ) 3 ∂ x F 1 ζ 1 -ζ 2 -∂ x ζ 1 -ζ 2 dx ≤ (1 + ζ L ∞ ) 3 ∂ l x ζ 1 -ζ L 2 ∂ x ζ 1 -ζ L 2 ≤ C( ζ H µ ) ∂ l x ζ L 2 ∂ x ζ L 2 ,
where the last inequality follows from Leibniz's rule and standard bilinear estimates

dE(ζ), ζ = R 2 h 1 + γh 2 h 1 h 2 ζ 2 - h 2 1 -γh 2 2 h 2 1 h 2 2 ζ 3 + 2 3 δ -1 h 3 2 ∂ x F 2 {h -1 2 ζ} ∂ x F 2 {h -2 2 ζ} + 2γ 3 h 3 1 ∂ x F 1 {h -1 1 ζ} ∂ x F 1 {h -2 1 ζ} + ζ h 2 ∂ x F 2 {h -1 2 ζ} 2 -γζ h 1 ∂ x F 1 {h -1 1 ζ} 2 dx, (2.2)
and we do not detail for the sake of conciseness.

Lemma 2.8. Let γ ≥ 0, δ > 0, µ > 1/2 and F i be admissible Fourier multipliers such that µ ≥ 1-θ.

Let ζ ∈ H µ (R) such that 1 -ζ L ∞ ≥ h 0 and δ -1 -ζ L ∞ ≥ h 0 , with h 0 > 0. Then one can decompose E(ζ) = E 2 (ζ) + E 3 (ζ) + E (1) rem (ζ) and dE(ζ), ζ = 2E 2 (ζ) + 3E 3 (ζ) + E (2) rem (ζ)
, where

E 2 (ζ) = R (γ + δ)ζ 2 + γ 1 3 (∂ x F 1 {ζ}) 2 + δ -1 1 3 (∂ x F 2 {ζ}) 2 dx, E 3 (ζ) = R (γ -δ 2 )ζ 3 -γζ(∂ x F 1 {ζ}) 2 + ζ(∂ x F 2 {ζ}) 2 + γ 2 3 (∂ x F 1 {ζ})(∂ x F 1 {ζ 2 }) - 2 3 (∂ x F 2 {ζ})(∂ x F 2 {ζ 2 }) dx.
Moreover, one has

E 2 (ζ) ≥ (γ + δ) ζ 2 L 2 and |E 3 (ζ)| ≤ C(h -1 0 , ζ H µ ) ζ L ∞ ζ 2 H 1-θ , ∀j ∈ {1, 2}, |E (j) rem (ζ)| ≤ C(h -1 0 , ζ H µ ) ζ 2 L ∞ ζ 2 H 1-θ , Proof.
The estimate on E 3 is straightforward by the Cauchy-Schwarz inequality and applying Lemma 2.4,i and Lemma 2.3,i. We detail the estimate for

|E (j) rem (ζ)|.
As above, we focus on the terms involving E, the terms involving E being obtained identically. One has

E (1) rem (ζ) = R ζ 4 1 -ζ + ζ 2 - ζ 3 3 ∂ x F 1 ζ 1 -ζ 2 dx + 1 3 R (1 -3ζ) ∂ x F 1 ζ 1 -ζ 2 -(∂ x F 1 {ζ}) 2 -2∂ x F 1 {ζ}∂ x F 1 {ζ 2 } + 3ζ(∂ x F 1 {ζ}) 2 dx.
We estimate each bracket separately. First note that, by the Cauchy-Schwarz inequality, Lemma 2.4,i and Lemma 2.3,iii,iv one has

R ζ 4 1 -ζ dx ≤ ζ 2 L ∞ ζ 2 L 2 h 0 , R ζ 2 - ζ 3 3 ∂ x F 1 ζ 1 -ζ 2 dx ≤ C(h -1 0 , ζ H µ ) ζ 2 L ∞ ζ 2 H 1-θ .
Next consider

1 3 R (1 -3ζ) ∂ x F 1 ζ 1 -ζ 2 -(∂ x F 1 {ζ}) 2 -2∂ x F 1 {ζ}∂ x F 1 {ζ 2 } + 3ζ(∂ x F 1 {ζ}) 2 dx = 1 3 R (1 -3ζ) ∂ x F 1 ζ 2 1 -ζ 2 + 2(∂ x F 1 {ζ})(∂ x F 1 ζ 3 1 -ζ ) -6ζ(∂ x F 1 {ζ}) ∂ x F 1 ζ 2 1 -ζ dx.
where we used the identity

∂ x F 1 ζ 1-ζ = ∂ x F 1 {ζ} + ∂ x F 1 ζ 2 1-ζ .
It follows from the Cauchy-Schwarz inequality, Lemma 2.4,i and Lemma 2.3,i,iv that

1 3 R ∂ x F 1 ζ 1 -ζ 2 -(∂ x F 1 {ζ}) 2 -2∂ x F 1 {ζ}∂ x F 1 {ζ 2 } dx ζ 2 1 -ζ 2 H 1-θ + ζ H 1-θ ζ 3 1 -ζ H 1-θ ≤ C(h -1 0 , ζ H µ ) ζ 2 L ∞ ζ 2 H 1-θ and R ζ ∂ x F 1 ζ 2 1 -ζ 2 + 2(∂ x F 1 {ζ}) ∂ x F 1 ζ 2 1 -ζ dx ζ L ∞ ζ 2 1 -ζ 2 H 1-θ + ζ L ∞ ζ H 1-θ ζ 2 1 -ζ H 1-θ ≤ C(h -1 0 , ζ H µ ) ζ 2 L ∞ ζ 2 H 1-θ .
This concludes the estimate for |E We do not detail for the sake of conciseness.

Periodic functional setting Given P > 0, we denote L 2 P the space of P -periodic, locally squareintegrable functions, endowed with the norm

u L 2 P = u L 2 (-P/2,P/2) def = P/2 -P/2 |u(x)| 2 dx 1 2
.

The Fourier coefficients of u ∈ L 2 P are defined by

ûk def = 1 √ P P/2 -P/2 u(x)e -2iπkx P dx, u(x) = 1 √ P k∈Z ûk e 2iπkx P .
We define, for s ≥ 0,

H s P def = u ∈ L 2 P , u 2 H s P def = k∈Z 1 + 4π 2 k 2 P 2 s |û k | 2 < ∞ .
The Fourier multiplier operator Λ :

S (R) → S (R) is defined as usual by Λ = (1 -∂ 2 x ) 1 2
. It maps periodic distributions to periodic distributions and we have

Λu k = 1 + 4π 2 k 2 P 2 1 2
ûk .

Thus

∂ x F i u k = 2πik P F i (2πk/P )û k .
For any s > 1/2, the continuous embedding where

u L ∞ ≤ 1 √ P k∈Z |û k | ≤ u H s P × 1 √ P k∈Z 1 (1 + 4π 2 k 2 P 2 ) s 1
E P (ζ) = P/2 -P/2 ζ 2 1 -ζ + 1 3 (1 -ζ) 3 ∂ x F 1 { ζ 1 -ζ } 2 dx, E P (ζ) = P/2 -P/2 ζ 2 δ -1 + ζ + 1 3 (δ -1 + ζ) 3 ∂ x F 2 { ζ δ -1 + ζ } 2 dx.

The periodic problem

Our first task is to construct periodic traveling-wave solutions with large periods by considering the periodic minimization problem corresponding to (1.10). We will use this in the next section to construct a special minimizing sequence for (1.10), which is useful when ν > 1 -θ. When θ < 1/2 and ν = 1 -θ, any minimizing sequence has the special property and therefore it is strictly speaking unnecessary to first consider the periodic minimization problem. Nevertheless, we consider here all possible parameters in order to highlight some interesting differences between the cases ν = 1 -θ and ν > 1 -θ.

We ensure that the hypotheses of Section 2, namely

ζ ∈ H ν P and ζ L ∞ < min(1, δ -1 )
will be satisfied through a penalization argument. To this aim, we fix R > 0 and restrict ourselves to values q ∈ (0, q 0 ) sufficiently small so that ζ

H ν P (R) ≤ 2R and (γ + δ) ζ 2 L 2 P = q ensures (by Lemma 2.1,i in the periodic setting) that ζ L ∞ < min(1, δ -1
) -h 0 with some h 0 > 0, uniformly with respect to P ≥ P 0 sufficiently large (and likewise in the real line setting). We then define : [0, (2R) 2 ) → [0, ∞) a smooth, non-decreasing penalization function, satisfying i. (t) = 0 for 0 ≤ t ≤ R 2 ;

ii. (t) → ∞ as t (2R) 2 ;

iii. For any a 1 ∈ (0, 1), there exists M 1 , M 2 > 0 and a 2 > 1 such that

(t) ≤ M 1 (t) a1 + M 2 (t) a2 ; (3.1) for instance : (R 2 , (2R) 2 ) t → ((2R) 2 -t) -1 exp( 1 R 2 -t ). Now consider the functional E P, (ζ) def = ( ζ 2 H ν P ) + E P (ζ)
and the constraint set

V P,q,2R def = ζ ∈ H ν P , (γ + δ) ζ 2 L 2 P = q and ζ H ν P < 2R .
In the following we solve for q sufficiently small and P sufficiently large arg min

ζ∈V P,q,2R E P, (ζ). (3.2)
Lemma 3.1. There exists q 0 > 0 such that for any q ∈ (0, q 0 ), the functional E P, : V P,q,2R → R is weakly lower semi-continuous, bounded from below and

E P, (ζ) → ∞ as ζ H ν P 2R.
In particular, it has a minimizer ζ P ∈ V P,q,2R , which satisfies the Euler-Lagrange equation

2 ( ζ P 2 H ν P )Λ 2ν ζ P + dE P (ζ P ) + 2α P (γ + δ)ζ P = 0 (3.3)
for some Lagrange multiplier α P (ζ P ) ∈ R.

Proof. As explained above, we restrict ourselves to q ∈ (0, q 0 ) so that E P, is well-defined, and in particular

sup ζ∈V P,q,2R ζ L ∞ < min(1, δ -1
( ζ 2 H ν P ) ≤ lim inf n→∞ ( ζ n 2 H ν P ). Now, by the Rellich-Kondrachov theorem, ζ n ζ in H ν P implies that ζ n → ζ in H s P , for s ∈ (1/2, ν); and in particular, ζ n → ζ in L ∞ . Moreover, one has sup n ζ n L ∞ < min(1, δ -1 ) since ζ n ∈ V P,q,2R , and therefore ζn 1-ζn → ζ 1-ζ in L ∞ ,
∂ x F 1 ζn 1-ζn ∂ x F 1 ζ 1-ζ in L 2
P by Lemma 2.4,i, and finally

(1 -ζ n ) 3 2 ∂ x F 1 ζ n 1 -ζ n (1 -ζ) 3 2 ∂ x F 1 ζ 1 -ζ in L 2 P .
In the same way, we find

(δ -1 + ζ n ) 3 2 ∂ x F 2 ζ n δ -1 + ζ n (δ -1 + ζ) 3 2 ∂ x F 2 ζ δ -1 + ζ in L 2 P .
and

(1 -ζ n ) -1 2 ζ n → (1 -ζ) -1 2 ζ, (δ -1 + ζ n ) -1 2 ζ n → (δ -1 + ζ) -1 2 ζ in L 2 P .
The result follows from the weak lower semi-continuity of • L 2 P . Now we wish to prove that ζ P ∈ V P,q,R , and in particular satisfies the Euler-Lagrange equation

dE P (ζ P ) + 2α P (γ + δ)ζ P = 0.
From this point on, we heavily make use of the property (see Assumption 1.4)

γ -δ 2 = 0.
without explicit references in the statements. Lemma 3.2. There exists m > 0 and q 0 > 0 such that for any q ∈ (0, q 0 ),

I q def = inf{E(ζ), ζ ∈ V q,R } < q(1 -mq 2 3 )
and there exists P q > 0 such that

I P, ,q def = inf{E P, (ζ), ζ ∈ V P,q,2R } < q(1 -mq 2 3 )
for any P ≥ P q .

Proof. Let us first consider the case of the real line. Consider ψ ∈ C ∞ (R) with compact support, such that (γ + δ) ψ

2 L 2 = 1; and denote ψ λ : x → λ 1 2 ψ(λx). One has R ψ 3 λ dx = λ 1 2 R ψ 3 dx and ∂ x ψ λ L 2 = λ ∂ x ψ L 2 .
It follows that, for the case when γ -δ 2 < 0, one can choose ψ ≥ 0 and λ small enough so that

R (γ -δ 2 )ψ 3 λ + γ + δ -1 3 (∂ x ψ λ ) 2 dx def = -2m < 0.
If γ -δ 2 > 0, we instead let ψ ≤ 0 and again choose λ small enough so that the above holds. Now, consider φ q : x → q 2 3 ψ λ (q 1 3 x). One has

φ q L ∞ ≤ q 2 3 ψ λ L ∞ , R φ 3 q dx = q 5 3 R ψ 3 λ dx and ∂ n x φ q L 2 = q 1 2 + n 3 ∂ n x ψ λ L 2 (n ∈ N).
In particular, for q sufficiently small, φ q H ν < R; and by Lemma 2.7 with l = 3,

E(φ q ) = (γ + δ) R φ 2 q dx + R (γ -δ 2 )φ 3 q + γ + δ -1 3 (∂ x φ q ) 3 dx + O(q 7 3 ) = q -2mq 5 3 + O(q 7
3 ).

The result follows in the real-line setting.

We now deduce the result in the periodic setting. By taking P ≥ P q sufficiently large, we may ensure supp φ q ⊂ (-P/2, P/2), and define φ P,q = j∈Z φ q (x -jP ) ∈ H n P (n ∈ N).

One has (γ + δ) φ P,q

2 L 2 P = (γ + δ) φ q 2 L 2 = q
and φ P,q H ν P < R for q sufficiently small, so that as well as E P, (φ P,q ) = E P (φ P,q ) = (γ + δ)

P/2 -P/2 φ 2 P,q dx + P/2 -P/2 (γ -δ 2 )φ 3 P,q + γ + δ -1 3 (∂ x φ P,q ) 2 dx + O(q 7 3 ) = q -2mq 5 3 + O(q 7
3 ).

The result is proved.

Lemma 3.3. There exists q 0 > 0 such that for any q ∈ (0, q 0 ), one has

∀P ≥ P q , |α P + 1| < 1 2 ,
where α P is defined in Lemma 3.1 and P q in Lemma 3.2.

Proof. We use the Euler-Lagrange equation satisfied by ζ P , namely

Let us now use Lemma 3.2, which asserts

( ζ P 2 H ν P ) + E P (ζ P ) < q(1 -mq 2 3 ) ≤ q. (3.5) Remark that, since 1 -ζ, δ -1 + ζ ≥ h 0 > 0, one has E P (ζ P ) ≥ P/2 -P/2 γ ζ 2 P 1 -ζ P + ζ 2 P δ -1 + ζ P dx = (γ + δ) P/2 -P/2 ζ 2 P dx + γ P 2 -P 2 ζ 3 P 1 -ζ P dx -δ P 2 -P 2 ζ 3 P δ -1 + ζ P dx = q + O(q 1+ 2ν ),
where = ν -1/2 > 0 and we use in the last estimate that

ζ P 2 L ∞ ζ P 2-1 ν L 2 P ζ P 1 ν H ν P = O(q 2ν-1 2ν
), by the interpolation estimate Lemma 2.1,i in the periodic-setting.

Combining with (3.5) yields

E P (ζ P ) = q + O(q 1+ 2ν )
and

( ζ P 2 H ν P ) = O(q 1+ 2ν ).
Hence, by (3.1)

( ζ P 2 H ν P ) = O(q 1+ 4ν ).
Now, by Lemma 2.5 and using once again (3.5), one has

ζ P 2 H 1-θ P E P (ζ P ) = O(q).
Thus by Lemma 2.8, one has

|E 3,P (ζ P )| ζ P L ∞ ζ P 2 H 1-θ P = O(q 1+ 2ν )
and

|E (2) rem,P (ζ P )| + |E (1) rem,P (ζ P )| ζ P 2 L ∞ ζ P 2 H 1-θ P = O(q 1+ ν ).
Plugging the above estimates into (3.4) yields

-α P q = q + O(q 1+ 4ν ),
and the proof is complete.

Lemma 3.4. Let q 0 , q ∈ (0, q 0 ) and P q be as in Lemma 3.3. There exists M > 0 such that one has

ζ 2 H ν P ≤ M q
uniformly over q ∈ (0, q 0 ), P ≥ P q and ζ in the set of minimizers of E P, over V P,q,2R .

Proof. Recall the Euler-Lagrange equation:

2 ( ζ P 2 H ν P )Λ 2ν ζ P + dE P (ζ P ) + 2α P ζ P = 0. (3.6)
It follows from the proof of Lemma 3.3 that for q 0 sufficiently small

ζ P 2 H 1-θ P q with 0 ≤ θ < 1.
Thus the result is proved if ν = 1-θ, and we focus below on the situation ν > 1-θ.

In this case we obtain the desired estimate in a similar fashion after finite induction. Indeed, define

r n = min(ν -(1 -θ), n(1 -θ))
, n ∈ N, and assume that ζ P 2 H rn P q. Note that this is satisfied for n = 0 by assumption. We will show below that

ζ P 2 H 1-θ+rn P ζ P 2 H rn P q.
(3.7)

Since 1 -θ > 0, the desired result follows by finite induction.

Let us now prove (3.7

where we define dE P and dE P from (1.9) as in (1.6). In particular,

dE(ζ P ), Λ 2rn ζ P = 2ζ P 1 -ζ P , Λ 2rn ζ P + ζ 2 P (1 -ζ P ) 2 , Λ 2rn ζ P + (1 -ζ P ) 2 ∂ x F ζ P 1 -ζ P 2 , Λ 2rn ζ P + 2 3 (1 -ζ P ) 3 ∂ x F 1 ζ P 1 -ζ P , ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 . (3.10)
We estimate each term of (3.10), using that (γ+δ) It is clear that the second term can be estimated in the same way. Next we see that

ζ P 2 L 2 P = q, ζ P H ν P < 2R and ζ P L ∞ < 1-h 0 ,
(1 -ζ P ) 2 ∂ x F ζ P 1 -ζ P 2 , Λ 2rn ζ P ≤ (1 -ζ P ) 2 ∂ x F 1 ζ P 1 -ζ P 2 H θ-1+rn P ζ P H 1-θ+rn P ∂ x F 1 ζ P 1 -ζ P 2 H θ-1+rn P ζ P H 1-θ+rn P ,
and Lemma 2.3,ii, Lemma 2.4,i and Lemma 2.1,ii yield for any 0 < < min(ν + θ -1, 1 -θ, ν -1/2),

∂ x F 1 ζ P 1 -ζ P 2 H θ-1+rn P ∂ x F 1 ζ P 1 -ζ P H rn P ∂ x F 1 ζ P 1 -ζ P H ν+θ-1- P ζ P H 1-θ+rn P ζ P H ν- P ζ P H 1-θ+rn P q 2ν .
We therefore have that

(1 -ζ P ) 2 ∂ x F ζ P 1 -ζ P 2 , Λ 2rn ζ P q 2ν ζ P 2 H 1-θ+rn P .
(3.12)

We next consider the remaining term in (3.10) and note that

(1 -ζ P ) 3 ∂ x F 1 ζ P 1 -ζ P , ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 = ∂ x F 1 {ζ P }, ∂ x F 1 {Λ 2rn ζ P } + (-3ζ P + 3ζ 2 P -ζ 3 P )∂ x F 1 ζ P 1 -ζ P , ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 I + ∂ x F 1 ζ P 1 -ζ P -ζ P , ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 II + ∂ x F 1 {ζ P }, ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 -Λ 2rn ζ P III . (3.13)
First we see that, by Lemma 2.4,i,

ζ 2 H rn P + ∂ x F 1 {ζ P }, ∂ x F 1 {Λ 2rn ζ P } ζ 2 H 1-θ+rn P . (3.14)
We estimate I in (3.13) proceeding as previously:

ζ P ∂ x F 1 ζ P 1 -ζ P , ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 ∂ x F 1 ζ P 1 -ζ P H rn P ζ P ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 H -rn P ζ P H 1-θ+rn P ζ P H ν- P ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 H -rn P ζ P H ν- P ζ P 2 H 1-θ+rn P (3.15) q 2ν ζ P 2 H 1-θ+rn P , (3.16)
where we choose 0 < < min{ν -1/2, 1 -θ}. The remaining terms in I are of higher order and can be estimated in the same way. Next we estimate II:

∂ x F 1 ζ 2 P 1 -ζ P , ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 ζ 2 P H 1-θ+rn P Λ 2rn ζ P H 1-θ-rn P ζ P L ∞ ζ P 2 H 1-θ+rn P (3.17) q ν-1/2 2ν ζ P 2 H 1-θ+rn P , (3.18)
where we used Lemma 2.3,i and Lemma 2.1,i. Finally consider III: proceeding as above,

∂ x F 1 {ζ P }, ∂ x F 1 Λ 2rn ζ P (1 -ζ P ) 2 -Λ 2rn ζ P ζ P H 1-θ+rn P 2ζ P -ζ 2 P (1 -ζ P ) 2 Λ 2rn ζ P H 1-θ-rn P ≤ ζ P 2 H 1-θ+rn P 2ζ P -ζ 2 P (1 -ζ P ) 2 H ν- P (3.19)
q 2ν ζ P 2 H 1-θ+rn P , (3.20
+ ζ P 2 H 1-θ+rn P q 2ν .
Thus one may choose q sufficiently small, so that (3.7) holds. This concludes the proof.

We now collect the preceding results and deduce solutions of the non-penalized periodic problem.

Theorem 3.5 (Existence of periodic minimizers).

There exists q 0 > 0 such that for any q ∈ (0, q 0 ), one can define P q > 0 and the following holds. For each P ≥ P q , there exists ζ P ∈ V P,q,R such that

E P (ζ P ) = inf ζ∈V P,q,R E P (ζ) def = I P,q .
and the Euler-Lagrange equation holds with α P ∈ (-3/2, -1/2):

dE P (ζ P ) + 2α P (γ + δ)ζ P = 0. (3.21)
Furthermore, there exists M > 0, independent of q, such that

ζ P 2 H ν P ≤ M q
uniformly with respect to P ≥ P q .

Proof. From Lemma 3.4, any minimizer of E P, over V P,q,2R satisfies, for q 0 sufficiently small and P ≥ P q sufficiently large,

ζ P 2 H ν P ≤ M q < R 2 .
Thus the Euler-Lagrange equation (3.3) becomes (3.21), and the control on α P is stated in Lemma 3.3. Moreover, since E P, = E P over V P,q,R , ζ P minimizes E P over V P,q,R . The theorem is proved.

Remark 3.6. If θ ∈ [0, 1/2) and ν = 1 -θ, then the functional E P is coercive on V P,q,R by Lemma 2.5, and it isn't necessary to consider the penalized functional E P, to construct periodic minimizers. Indeed, one can minimize E P over V P,q,R directly, noting that any minimizing sequence satisfies (up to subsequences) sup n ζ P,n 2 H ν P ≤ M q < R 2 if q 0 is sufficiently small.

The real line problem

The construction of a minimizer for the real line problem, (1.10), will follow from Lions' concentration compactness principle. The main difficulty consists in excluding the "dichotomy" scenario. To this aim, we shall use a special minimizing sequence (satisfying the additional estimate ζ n 2 H ν q) to show that the function q → I q is strictly subhomogeneous (see Proposition 4.2), which implies that it is also strictly subadditive (Corollary 4.3). This special subsequence is constructed from the solutions of the periodic problem, obtained in Theorem 3.5, with period P n → ∞.

A special minimizing sequence

Theorem 4.1 (Special minimizing sequence for E). There exists q 0 > 0 such that for any q ∈ (0, q 0 ), one can define constants m, M > 0 and a sequence {ζ n } n∈N satisfying

(γ + δ) ζ n 2 L 2 = q, ζ n 2 H ν ≤ M q and lim n→∞ E(ζ n ) = I q def = inf ζ∈V q,R E(ζ) < q(1 -mq 2
3 ).

Proof. The estimate on I q was proved in Lemma 3.2; thus we only need to construct a minimizing sequence satisfying ζ n 2 H ν ≤ M q. If ν = 1-θ, then any minimizing sequence satisfies this property as a consequence of Lemma 2.5, so we assume in the following that ν > 1 -θ. Let q 0 be sufficiently small so that Theorem 3.5 holds. By the construction of [24, p. 2918 and proof of Theorem 3.8], one obtains, for any P n sufficiently large,

x n ∈ R, ζ Pn ∈ H ν Pn and ζ n ∈ H ν (R) such that ζ Pn -ζ Pn (• -x n ) L 2 Pn → 0 (P n → ∞) (4.1)
where ζ Pn is defined by Theorem 3. H ν ≤ M q < R 2 provided that P n is sufficiently large and q 0 is sufficiently small; and ζ n ∈ V q,R by (4.3). Thus there only remains to prove that ζ n is a minimizing sequence.

Notice that by (4.1) and (4.4) and Lemma 2.1,ii,

ζ Pn -ζ Pn (• -x n ) H ν Pn → 0 (4.5)
for any ν ∈ [0, ν).

One has

E Pn ( ζ Pn ) -E(ζ n ) = 1 3 Pn/2 -Pn/2 (1 -ζ n ) 3 ∂ x F 1 { ζ Pn 1 -ζ Pn } 2 -∂ x F 1 { ζ n 1 -ζ n } 2 dx - 1 3 R\(-Pn/2,Pn/2) (1 -ζ n ) 3 ∂ x F 1 { ζ n 1 -ζ n } 2 dx.
Using Lemma 2.4,i, (4.4) and the Cauchy-Schwarz inequality, we deduce

Pn/2 -Pn/2 (1 -ζ n ) 3 ∂ x F 1 { ζ Pn 1 -ζ Pn } 2 -∂ x F 1 { ζ n 1 -ζ n } 2 dx C( ζ n H ν ) ∂ x F 1 { ζ Pn 1 -ζ Pn } -∂ x F 1 { ζ n 1 -ζ n } L 2 (-Pn/2,Pn/2) .
Notice now that, by uniqueness of the Fourier decomposition in L 2 Pn , one has the identity

∂ x F 1 { ζ Pn 1 -ζ Pn }(x) = l∈Z ∂ x F 1 { ζ n 1 -ζ n } (x + lP n ),
and therefore, by Lemma 2.4,iii and (4.2), one has

∂ x F 1 { ζ Pn 1 -ζ Pn } -∂ x F 1 { ζ n 1 -ζ n } 2 L 2 (-Pn/2,Pn/2) = Pn/2 -Pn/2   |l|≥1 |∂ x F 1 { ζ n 1 -ζ n }(y + lP n )|   2 dy Pn/2 -Pn/2   |l|≥1 1 (P 1 2 n + (l -1)P n ) j   2 dy → 0 (P n → ∞).
since j ≥ 2. Similarly, Lemma 2.4,iii and (4.2) yield R\(-Pn/2,Pn/2)

(1 -ζ n ) 3 ∂ x F 1 { ζ n 1 -ζ n } 2 dx ≤ C( ζ H ν ) ∞ Pn/2 1 (x -P n /2 + P 1 2 n ) 2j dx → 0 (P n → ∞).
The component E satisfies the same bounds, thus we proved

E Pn ( ζ Pn ) -E(ζ n ) → 0 (P n → ∞).
Now by Lemma 2.6 (which holds in the periodic setting and uniformly with respect to P > 0) with ν replaced by some ν ∈ (1/2, ν) and (4.5), one has

E Pn ( ζ Pn ) -I Pn,q = E Pn ( ζ Pn ) -E Pn (ζ Pn (• -x n )) → 0 (P n → ∞).
Thus we found that

I q ≤ E(ζ n ) = I Pn,q + o(1) (P n → ∞).
There remains to prove the converse inequality. For any > 0, there exists ζ ∈ V q,R such that

E(ζ) ≤ I q + 3 .
By the same argument as above, we construct by smoothly truncating and rescaling, ζ ∈ V q,R such that supp ζ ∈ (-P , P ), and

E( ζ) ≤ E(ζ) + 3 .
Then for P n ≥ 2P , one has ζPn = j∈Z ζ(• + jP n ) ∈ V P,q,R and, as above,

E Pn ( ζPn ) -E( ζ) → 0 (P n → ∞).
Combining the above yields, for P n sufficiently large,

I Pn,q ≤ E Pn ( ζPn ) ≤ E( ζ) + 3 ≤ I q + .
Thus we proved that

E(ζ n ) → I q (P n → ∞).
This concludes the proof.

The following proposition is essential to rule out the "dichotomy" scenario in Lions' concentrationcompactness principle (see below). Proposition 4.2. There exists q 0 > 0 such that the map q → I q is strictly subhomogeneous for q ∈ (0, q 0 ): I aq < aI q whenever 0 < q < aq < q 0 .

Proof. Let us consider ζ n the special minimizing sequence defined in Theorem 4.1. We first fix a 0 > 1, and restrict q 0 > 0 if necessary, so that for any a ∈ (1, a 0 ] and q ∈ (0, q 0 ) such that aq < q 0 , one has a

1 2 ζ n 2 H ν ≤ M aq ≤ M q 0 < R 2 .
Thus we have, by definition of I aq and Lemma 2.8,

I aq ≤ E(a 1 2 ζ n ) = aE 2 (ζ n ) + a 3 2 E 3 (ζ n ) + E (1) rem (a 1 2 ζ n ) = aE(ζ n ) + (a 3 2 -a)E 3 (ζ n ) + E (1) rem (a 1 2 ζ n ) -aE (1) rem (ζ n ).
Moreover, by Theorem 4.1, one has

lim n→∞ E(ζ n ) = I q < q(1 -mq 2
3 ), and Lemma 2.8 yields

-E 2 (ζ n ) ≤ -q and E (1) rem (ζ n ) ζ n 4 H ν a 2 0 q 2 .
It follows that one has for q ∈ (0, q 0 ) with q 0 sufficiently small and n sufficiently large,

E 3 (ζ n ) = E(ζ n ) -E 2 (ζ n ) -E (1) rem (ζ n ) ≤ - 1 2 mq 5 3 .
Thus we find for n sufficiently large,

I aq ≤ aI q -(a 3 2 -a)(m/2)q 5 3 + lim sup n→∞ E (1) rem (a 1 2 ζ n ) -aE (1) rem (ζ n ) . (4.6)
We now estimate the last contribution, treating separately E

(1)

rem and E (1) rem in the same spirit as in the proof of Lemma 2.8. Consider E (1) rem for instance. We develop each contribution in

E (1) rem (a 1 2 ζ n ) using Neumann series in powers of a 1 2 ζ n : E (1) rem (a 1 2 ζ n ) = R k≥4 (a 1 2 ζ n ) k dx + R k1+k2+k3≥4 c k1,k2,k3 (a 1 2 ζ n ) k1 (∂ x F 1 {(a 1 2 ζ n ) k2 })(∂ x F 1 {(a 1 2 ζ n ) k3 }) dx.
The series are absolutely convergent provided q 0 is sufficiently small, and start at index k = 4, as pointed out in the proof of Lemma 2.8. We now subtract the contributions of aE

rem (a

1 2 ζ n ) -aE (1) rem (ζ n )| ≤ k≥4 (a k 2 -a) ζ n k-2 L ∞ ζ n 2 L 2 + k1+k2+k3≥4 |c k1,k2,k3 |(a k 1 +k 2 +k 3 2 -a) ζ n k1 L ∞ ∂ x F 1 {ζ k2 n } L 2 ∂ x F 1 {ζ k3 n } L 2 .
Using that

|a k 2 -a| ≤ (a 3 2 -a)(k -2)a k-3
2 , Lemma 2.4,i, that H ν is a Banach algebra as well as the continuous embedding H ν ⊂ L ∞ , we find that one can restrict q 0 > 0 such that the above series is convergent and yields

|E (1) rem (a 1 2 ζ n ) -aE (1) rem (ζ n )| ≤ C(a 0 )(a 3 2 -a)q 2 ,
uniformly over q ∈ (0, q 0 ) and a ∈ (1, a 0 ] such that aq < q 0 . Plugging this estimate in (4.6) and restricting q 0 if necessary, we deduce I aq < aI q for 0 < q < aq < q 0 , a ∈ (1, a 0 ].

Consider now the case when a ∈ (1, a p 0 ] for an integer p ≥ 2. Then a 1 p ∈ (1, a 0 ] and so

I aq = I a 1 p a p-1 p q < a 1 p I a p-1 p q = a 1 p I a 1 p a p-2 p q < a 2 p I a p-2 p q < . . . < aI q .
The result is proved.

By a standard argument, Proposition 4.2 induces the subadditivity of the map q → I q .

Corollary 4.3. There exists q 0 > 0 such that the map q → I q is strictly subadditive for q ∈ (0, q 0 ):

I q1+q2 < I q1 + I q2 whenever 0 < q 1 < q 1 + q 2 < q 0 .

Concentration-compactness; proof of Theorem 1.6

We now prove Theorem 1.6. Let us first recall Lions' concentration compactness principle [START_REF] Lions | The concentration-compactness principle in the calculus of variations. The locally compact case. I[END_REF]. We shall apply Theorem 4.4 to

e n = γ ζ 2 n 1 -ζ n + 1 3 (1 -ζ n ) 3 ∂ x F 1 { ζ n 1 -ζ n } 2 + ζ 2 n δ -1 + ζ n + 1 3 (δ -1 + ζ n ) 3 ∂ x F 2 { ζ n δ -1 + ζ n } 2 ,
where

ζ n is a minimizing sequence of E over V q,R with sup n ζ n 2 H ν < R 2 .
Such a sequence is known to exist provided that q ∈ (0, q 0 ) is sufficiently small, by Theorem 4.1 (and any minimizing sequence is valid when ν = 1 -θ, by Lemma 2.5; see Remark 3.6). The choice of density is inspired by the recent paper [START_REF] Arnesen | Existence of solitary-wave solutions to nonlocal equations[END_REF], and allows (contrarily to the more evident choice e n = ζ 2 n ) to show, when ν = 1-θ, that the constructed limit satisfies E(η) = I q and is therefore a solution to the constrained minimization problem (1.10). Notice that

R e n dx = E(ζ n ) → I q (n → ∞)
and that there exists a constant C such that

ζ n 2 L 2 (J) = J |ζ n | 2 dx ≤ C J e n dx (4.7)
for any interval J ⊆ R.

We exclude the two first scenarii in Lemmata 4.5 and 4.6, below. Thus the concentration scenario holds and, using (4.7), we find that there exists {x n } n∈N ⊂ R such that for any > 0, there exists r > 0 with η n L 2 (|x|>r) < , Now, consider χ a smooth cut-off function, such that χ(x) = 1 for |x| ≤ r and χ(x) = 0 for |x| ≥ 2r. One has χη n χη weakly in H ν (R), and by compact embedding [3, Corollary 2.96], we find that one can extract a subsequence, still denoted η n , such that χ(η n -η) L 2 (R) ≤ for n sufficiently large. Combining the above estimates, we find that the subsequence satisfies η n -η L 2 (R) < 3 for n sufficiently large. By Cantor's diagonal extraction process, we construct a subsequence satisfying η n -η L 2 → 0; and by interpolation,

where {η n } n∈N def = {ζ n (• + x n )} n∈N . Since sup n∈N η n H ν (R) < R, there exists η ∈ H ν (R) satisfying η H ν (R) < R
η n -η H s → 0 for any s ∈ [0, ν). In particular (γ + δ) η 2 L 2 = q, and recall η H ν ≤ sup n ζ n H ν < R, thus η ∈ V q,R . If ν > 1 -θ, we deduce E(η n ) → E(η)
as n → ∞ by Lemma 2.6. If on the other hand ν = 1-θ we use the weak lower semi-continuity argument in the proof of Lemma 3.1 to deduce that I q ≤ E(η) ≤ lim n→∞ E(η n ) = I q . In either case we have that E(η) = I q .

The constructed function η ∈ H ν (R) is therefore a solution to the constrained minimization problem (1.10). In particular, it solves the Euler-Lagrange equation (1.11) with α < 0 provided that q ∈ (0, q 0 ) is sufficiently small (proceeding as in Lemma 3.3), and therefore satisfies (1.8) with

c 2 = (-α) -1 > 0.
This proves the first item of Theorem 1.6, as well as the second item -except for the strong convergence in H ν (R) when ν = 1-θ > 1/2. This result follows from the fact that weak convergence together with convergence of the norm implies strong convergence in a Hilbert space (applied to

(γ 1/2 (1 -ζ n ) 3/2 ∂ x F 1 { ζn 1-ζn } , (δ -1 + ζ n ) 3/2 ∂ x F 2 { ζn δ -1 +ζn } ) ∈ (L 2 (R)) 2
). There remains to prove the estimates of the third item. Proceeding as in Lemma 3.4, we find

ζ 2 H ν ≤ M q,
uniformly over the minimizers of E over V q,R . Moreover, by Lemma 2.8, one has

-α q = -α(γ + δ) ζ 2 L 2 = 1 2 dE(ζ), ζ = E 2 (ζ) + 3 2 E 3 (ζ) + 1 2 E (2) rem (ζ) = 3 2 E(ζ) - 1 2 E 2 (ζ) - 3 2 E (1) rem (ζ) + 1 2 E (2) rem (ζ).
where

E 2 (ζ) ≥ q and |E (1) rem (ζ)| + |E (2) rem (ζ)| = O(q 2 ). Altogether, using that E(ζ) < q(1 -mq 2 3
) by Lemma 3.2, we find

-α q < 3 2 q(1 -mq 2 3 ) - 1 2 q + O(q 2 ) = q 1 - 3 2 mq 2 3 + O(q 2 ),
and the result follows. Theorem 1.6 is proved.

Lemma 4.5 (Excluding "vanishing"). No subsequence of {e n } n∈N has the "vanishing" property.

Proof. By Lemmata 2.8 and 3.2, one has for n sufficiently large

q(1 -mq 2 3 ) > E(ζ n ) = E 2 (ζ n ) + E 3 (ζ n ) + E (1) rem (ζ n ) ≥ q + E 3 (ζ n ) + E (1) rem (ζ n )
and hence mq

5 3 ≤ |E 3 (ζ n )| + |E (1) rem (ζ n )| ζ n L ∞ .
On the other hand, one has

ζ n L ∞ ((x-1 2 ,x+ 1 2 )) ≤ ϕ x ζ n L ∞ (R) ≤ ϕ x ζ n 1-1 2ν L 2 (R) ϕ x ζ n 1 2ν H ν (R) ≤ C ζ n 1-1 2ν L 2 ((x-1,x+1) ζ n 1 2ν H ν (R) ,
where ϕ x = ϕ(• -x) with ϕ a smooth function such that ϕ = 1 for |x| ≤ 1/2, ϕ = 0 for |x| ≥ 1, and 0 ≤ ϕ ≤ 1 otherwise; and using Lemma 2.1,i and Lemma 2.

3,ii. Since C is independent of x ∈ R, this shows that ζ n L ∞ ≤ CR 1 2ν sup x∈R ζ n 1-1 2ν L 2 ((x-1,x+1)) .
Hence one has for n sufficiently large

q 5 3 sup x∈R ζ n 1-1 2ν L 2 ((x-1,x+1)) ,
from which, using (4.7), it follows that "vanishing" cannot occur.

Lemma 4.6 (Excluding "dichotomy"). No subsequence of {e n } n∈N has the "dichotomy" property.

Proof. We denote by χ ∈ C ∞ (R + ) a non-increasing function with

χ(r) = 1 if 0 ≤ r ≤ 1 and χ(r) = 0 if r ≥ 2, (4.8)
and such that

χ = χ 2 1 , 1 -χ = χ 2 2
where χ 1 and χ 2 are smooth. For instance, set χ(r) = 1 -(1 -χ 2 (r)) 2 with χ ∈ C ∞ (R + ) nonincreasing and satisfying (4.8). Define

η n = ζ n (• + x n ), and
η (1) n (x) = η n (x)χ(|x|/M n ) and η (2) n (x) = η n (x) 1 -χ(2|x|/N n ) , noting that supp(η (1) n ) ⊂ [-2M n , 2M n ] and supp(η (2) n ) ⊂ R \ [-N n /2, N n /2].
After possibly extracting a subsequence, we can assume that

η (1) n 2 L 2 → q * γ + δ (4.9)
with q * ∈ [0, q]. For n sufficiently large, one has

N n > N n /2 > 2M n > M n and therefore η (2) n 2 L 2 = η n 2 L 2 (|x|>Nn) + η (2) n 2 L 2 (Mn<|x|<Nn) = q γ + δ -η (1) n 2 L 2 -η n 2 L 2 (Mn<|x|<Nn) + η (1) n 2 L 2 (Mn<|x|<Nn) + η (2) n 2 L 2 (Mn<|x|<Nn) → q -q * γ + δ (4.10) since η (1) n 2
L 2 (Mn<|x|<Nn) + η (2) n 2 L 2 (Mn<|x|<Nn) ≤ 2 η n 2 L 2 (Mn<|x|<Nn) and η n 2 L 2 (Mn<|x|<Nn) ≤ C Mn<|x-xn|<Nn
e n dx → 0 by (4.7) and the assumption of the dichotomy scenario.

We claim that E(η

n ) → I * . To show this, note that

E(η (1) n ) = R (η (1)
n ) 2 1 -η (1) n + 1 3 (1 -η (1) n ) 3 ∂ x F 1 { η (1) n 1 -η (1) n } 2 dx, where R (η (1)
n ) 2 1 -η (1) n dx - |x|≤Mn η 2 n 1 -η n dx Mn≤|x-xn|≤Nn η 2 n dx → 0.
We next find that

R (1 -η (1) n ) 3 ∂ x F 1 { η (1) n 1 -η (1) n } 2 -(1 -η (1) n ) 3 ∂ x F 1 { η (1) n 1 -η n } 2 dx = R (1 -η (1) n ) 3 ∂ x F 1 { η (1) n 1 -η (1) n } + ∂ x F 1 { η (1) n 1 -η n } ∂ x F 1 η (1) n (η (1)
n -η n ) (1 -η (1) n )(1 -η n ) dx.
Noting that η (1) n (η (1) n -

η n ) = -χ 2 1 (|x|/M n )χ 2 2 (|x|/M n )η 2 n so that (γ + δ) η (1) n 2 L 2 = q and (γ + δ) η (2) n 2
L 2 = q -q * for any n ∈ N. One easily checks that

lim sup n→∞ η (1) n H ν < R, lim sup n→∞ η (2) n H ν < R and that lim n→∞ (E( η (1) n ) -E(η (1) n )) = lim n→∞ (E( η (2) n ) -E(η (2) n )) = 0.
Thus we arrive at the following contradiction to Corollary 4.3:

I q < I q * + I q-q * ≤ lim n→∞ (E( η (1) n ) + E( η (2) n )) = I * + I q -I * = I q .
This concludes the proof of Lemma 4.6.

Long-wave asymptotics

In this section we prove that the solutions of (1.8) obtained in Theorem 1.6 are approximated by solutions of the corresponding KdV equation in the long-wave regime, i.e. letting q → 0 in the constrained minimization problem (1.10). Indeed, if we introduce the scaling

ζ(x) = S KdV (ξ)(x) def = q 2 3 ξ(q 1 3 x) (5.1)
in (1.11) and denote α + 1 = α 0 q 2 3 , then we find that the leading order part of the equation as q → 0 is

α 0 (γ + δ)ξ + 3(γ -δ 2 )ξ 2 2 - (γ + δ -1 ) 3 ∂ 2 x ξ = 0. (5.2)
Recall (see e.g. [START_REF] Pava | Nonlinear dispersive equations. Existence and stability of solitary and periodic travelling wave solutions[END_REF]) that ξ ∈ L 2 (R) satisfying (5.2) uniquely defines (up to spatial translation) a solitary-wave solution of the KdV equation, with explicit formula

ξ KdV (x) = α 0 (γ + δ) δ 2 -γ sech 2 1 2 3α 0 (γ + δ) γ + δ -1 x .
Equation (5.2) can also be obtained as the Euler-Lagrange equation associated with the minimizer of the scalar functional E KdV (consistently with Lemma 2.7)

E KdV (ξ) = R (γ -δ 2 )ξ 3 + (γ + δ -1 ) 3 (∂ x ξ) 2 dx, over the set U 1 def = {ξ ∈ H 1 (R) : (γ + δ) ξ 2 L 2 = 1}
. Indeed, any minimizer satisfies the Euler-Lagrange equation

dE KdV (ξ) + 2(γ + δ)α 0 ξ = 0, (5.3)
which is (5.2) with α 0 the Lagrange multiplier. Testing the constraint (γ + δ) ξ 2 L 2 = 1 with the above explicit formula, we find that

(γ + δ)α 0 = 3 4 (δ 2 -γ) 4 (γ + δ)(γ + δ -1 ) 1 3
.

(5.4)

Additional computations show that

I KdV = inf{E KdV (ξ) : ξ ∈ U 1 } = E KdV (ξ KdV ) = - 3 5 α 0 .
We aim at proving that the variational characterization of (5.2), and therefore its explicit solutions, approximate (after suitable rescaling) the corresponding one of (1.8), namely (1.10), in the limit q → 0.

Refined estimates

We start by establishing estimates on ζ ∈ D q,R the set of minimizers of E over V q,R , as provided by Theorem 1.6. Here and below, we rely on extra assumptions on the Fourier multipliers, which are assumed to be strongly admissible, in the sense of Definition 1.2.

Lemma 5.1. There exists q 0 > 0 such that ζ ∈ H s for any s ≥ 0, and there exists M s > 0 such that ζ 2 H s ≤ M s q uniformly for q ∈ (0, q 0 ) and ζ ∈ D q,R .

Proof. Once the regularity property ζ ∈ H s has been established, the corresponding estimate is obtained as in the proof of Lemma 3.4, thus we focus only on the regularity issue. This follows from the Euler-Lagrange equation (1.11) and elliptic estimates. However, the ellipticity property is not straightforward to ascertain when γ = 0, and we will make use of paradifferential calculus. These tools are recalled in Appendix A.

By assumption, one has ζ ∈ H ν with ν > 1/2 and ν ≥ 1 -θ > 0. We fix ∈ (0, ν -1/2) and r = min(1 -θ, ν -1/2 -) > 0. We show below that ζ ∈ H ν satisfying (1.11) yields ζ ∈ H ν+r , and the argument can be bootstrapped to obtain arbitrarily high regularity, ζ ∈ H s , s ≥ 0.

First we write (1.11) as the equality, valid in H -ν ,

2 3 h -2 2 ∂ x F 2 h 3 2 ∂ x F 2 {h -1 2 ζ} + 2γ 3 h -2 1 ∂ x F 1 h 3 1 ∂ x F 1 {h -1 1 ζ} = 2α(γ + δ)ζ + 2 h 1 + γh 2 h 1 h 2 ζ - h 2 1 -γh 2 2 h 2 1 h 2 2 |ζ| 2 + h 2 ∂ x F 2 {h -1 2 ζ} 2 -γ h 1 ∂ x F 1 {h -1 1 ζ} 2 def = R(ζ) (5.5) denoting h 1 = 1 -ζ, h 2 = δ -1 + ζ,
(ζ) ∈ H 2(ν-(1-θ))-1/2-in the case 1/2 < ν ≤ 1/2 + (1 -θ), and R(ζ) ∈ H ν-(1-θ) if ν > 1/2 + (1 -θ).
In other words, we find

R(ζ) ∈ H ν-2(1-θ)+r .
(5.6)

Above, we used that ζ L ∞ < min(1, δ -1 ) and therefore h 1 (x) n -1 ∈ H ν and h 2 (x) n -(δ -1 ) n ∈ H ν for any n ∈ Z. This holds as well in the Hölder space W r,∞ since r ∈ (0, ν -1/2). In particular, we have

∀n ∈ Z, h 1 (x) n ∈ Γ 0 r and ∂ x F i ∈ Γ 1-θ r , recalling Definition A.1. By Lemma A.6, we find ζh -1 1 -T h -2 1
ζ ∈ H ν+r , and Lemma 2.4 and Lemma 2.3 yield

h -2 1 ∂ x F 1 h 3 1 ∂ x F 1 {ζh -1 1 -T h -2 1 ζ} ∈ H ν-2(1-θ)+r .
(5.7) -θ) . We deduce by Lemma A.4 that

By Lemma A.3, one has ∂ x F 1 T h -2 1 ζ = T ikF1(k) T h -2 1 ζ ∈ H ν-(1
ikF 1 (k)h -2 1 ∈ Γ 1-θ r and ∂ x F 1 T h -2 1 ζ -T ikF1(k)h -2 1 ζ ∈ H ν-(1-θ)+r ,
from which we deduce as above -θ) and Lemma A.5, one obtains

h -2 1 ∂ x F 1 h 3 1 ∂ x F 1 {T h -2 1 ζ} -T ikF1(k)h -2 1 ζ ∈ H ν-2(1-θ)+r . (5.8) Using that T ikF1(k)h -2 1 ζ ∈ H ν-(1
(h 3 1 -T h 3 1 )T ikF1(k)h -2 1 ζ ∈ H ν-(1-θ)+r .
As above, it follows by Lemma 2.4 and Lemma 2.3 that

h -2 1 ∂ x F 1 (h 3 1 -T h 3 1 )T ikF1(k)h -2 1 ζ ∈ H ν-2(1-θ)+r .
(5.9)

We use again Lemma A.3 and Lemma A.4 to deduce that -(kF

1 (k)) 2 h 1 ∈ Γ 2(1-θ) r and h -2 1 ∂ x F 1 T h 3 1 T ikF1(k)h -2 1 ζ -T -(kF1(k)) 2 h1 ζ ∈ H ν-2(1-θ)+r .
(5.10)

Finally, Lemma A.5 yields

(h -2 1 -T h -2 1 )T -(kF1(k)) 2 h1 ζ ∈ H ν-2(1-θ)+r (5.11)
and Lemma A.4 yields

(T h -2 1 T -(kF1(k)) 2 h1 -T -(kF1(k)) 2 h -1 1 )ζ ∈ H ν-2(1-θ)+r . (5.12)
Collecting (5.7)-(5.12), we proved

h -2 1 ∂ x F 1 h 3 1 ∂ x F 1 {h -1 1 ζ} -T -(kF1(k)) 2 h -1 1 ζ ∈ H ν-2(1-θ)+r .
(5.13) By (5.6), (5.13) and the corresponding estimate for the second contribution in the left-hand side of (5.5), one finds

T 2 3 h -1 2 (ikF2(k)) 2 + 2γ 3 h -1 1 (ikF1(k)) 2 ζ ∈ H ν-2(1-θ)+r . Moreover, since ζ ∈ H ν , one has 2 3 h -1 2 (x) + 2γ 3 h -1 1 (x) ∈ Γ 0 r and therefore T 2 3 h -1 2 (x)+ 2γ 3 h -1 1 (x) ζ ∈ H ν ⊂ H ν-2(1-θ)+r . Adding the two terms yields T a(x,k) ζ ∈ H ν-2(1-θ)+r (5.14) with a(x, k) def = 2 3 h -1 2 (x) 1 + (kF 2 (k)) 2 + 2γ 3 h -1 1 (x) 1 + (kF 1 (k)) 2 . Notice that a(x, k) ∈ Γ 2(1-θ) r and a(x, k) -1 ∈ Γ -2(1-θ) r .
In particular, Lemma A.3 and (5.14) yield

T a(x,k) -1 T a(x,k) ζ ∈ H ν+r .
Additionally, by Lemma A.4, we have

ζ -T a(x,k) -1 T a(x,k) ζ = T a(x,k) -1 a(x,k) ζ -T a(x,k) -1 T a(x,k) ζ ∈ H ν+r .
Adding the two terms shows that ζ ∈ H ν+r , which concludes the proof.

Remark 5.2. In the one-layer situation, namely γ = 0, the use of paradifferential calculus is not necessary, and Lemma 5.1 can be obtained through a direct use of Lemmata 2.2 and 2.3. In particular, Lemma 5.1 and subsequent results hold for (non-necessarily strongly) admissible Fourier multipliers, in the sense of Definition 1.1.

The following lemma shows that the minimizers of E over V q,R , as provided by Theorem 1.6, scale as (5.1).

Lemma 5.3. There exists q 0 > 0 and C > 0 such that the estimates

ζ L ∞ ≤ Cq 2 3 ,
(5.15)

∂ x ζ 2 L 2 ≤ Cq 5 3 , (5.16
)

∂ 2 x ζ 2 L 2 ≤ Cq 7 3
(5.17)

hold uniformly for q ∈ (0, q 0 ) and ζ ∈ D q,R , the set of minimizers of E over V q,R .

Proof. Let ζ be minimizer over V q,R . Since 2(γ + δ)αζ + dE(ζ) = 0, we get from Lemma 2.8 that

2α(γ + δ)ζ + dE 2 (ζ) = 2α(γ + δ)ζ + dE(ζ) -dE 3 (ζ) -dE (1) rem (ζ) = -dE 3 (ζ) -dE (1) rem (ζ), (5.18)
where

dE 3 (ζ) = γdE 3 (ζ) + dE 3 (ζ),
and

dE 3 (ζ) = 3ζ 2 -∂ x F 1 {ζ} 2 + 2∂ x F 1 {ζ∂ x F 1 {ζ}} - 2 3 ∂ x F 1 {∂ x F 1 {ζ 2 }} - 4 3 ζ∂ x F 1 {∂ x F 1 {ζ}}, dE 3 (ζ) = -3δ 2 ζ 2 + ∂ x F 2 {ζ} 2 -2∂ x F 2 {ζ∂F 2 {ζ}} + 2 3 ∂ x F 2 {∂ x F 2 {ζ 2 }} + 4 3 ∂ x F 2 {∂ x F 2 {ζ}}.
We also have that

dE 2 (ζ) = γdE 2 (ζ) + dE 2 (ζ) = 2(γ + δ)ζ - 2 3 γ∂ x F 1 {∂ x F 1 {ζ}} + δ -1 ∂ x F 2 {∂ x F 2 {ζ}} .
In frequency space equation (5.18) becomes

2 (γ + δ)α + γ + δ + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 ) ζ(k) = -F dE 3 (ζ) + dE (1) rem (ζ) (k).
By using the estimate for α in Theorem 1.6, we deduce

| ζ(k)| ≤ 1 2 F dE 3 (ζ) + dE (1) rem (ζ) (k) mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )
.

(5.19)

The estimates follow from (5.19) and a suitable decomposition into high-and low-frequency components. In order to estimate the right-hand-side, we heavily make use of Lemma 5.1:

ζ 2 H n
q for all n ∈ N. This will be used again throughout the proof without reference.

We first deduce from Lemma 2.3 that

F dE 3 (ζ) L 1 dE 3 (ζ) H 1 q and F dE 3 (ζ) L ∞ ζ 2 L 2 + ∂ x ζ 2 L 2 + ζ L 2 ∂ 2 x ζ L 2 q and, similarly, F dE (1) rem (ζ) L 1 + F dE (1) rem (ζ) L ∞ q
By the definition of admissible Fourier multipliers in (1.1), there exists c 0 , k 0 > 0 such that

∀k ∈ R \ [-k 0 , k 0 ], |k|F(k) ≥ c 0 .
We also assume that F(k) > 0, and therefore there exists c 0 > 0 such that

∀k ∈ [-k 0 , k 0 ], F(k) ≥ c 0 .
As a consequence, we have

sup k∈R\[-k0,k0] 1 mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )
1 and k0 -k0

1 mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 ) dk q -1 3 .

Now, we decompose

ζ L ∞ ≤ 1 √ 2π ζ L 1 ≤ 1 2 R 1 mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 ) F dE 3 (ζ) + dE (1) rem (ζ) (k) dk.
(5.20) into high-and low-frequency components and estimate each part. For the low frequency part we have k0 -k0

1 mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 ) F dE 3 (ζ) + dE (1) rem (ζ) (k) dk ≤ F dE 3 (ζ) + dE (1) rem (ζ) (k) L ∞ × k0 -k0 1 mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 ) dk q 2 3 .
For the high-frequency part

R\[-k0,k0] 1 mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 ) F dE 3 (ζ) + dE rem (ζ) (k) dk ≤ F dE 3 (ζ) + dE rem (ζ) (k) L 1 sup k∈R\[-k0,k0] 1 mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )
q.

Combining the above estimates in (5.20) gives us the inequality (5.15).

Let us now turn to (5.16). By (5.19) we have

∂ x ζ 2 L 2 ≤ 1 4 R k 2 |F(dE 3 (ζ) + dE (1) rem (ζ))| 2 (mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 dk (5.21)
We estimate the low-frequency part as above:

k0 -k0 k 2 |F(dE 3 (ζ) + dE (1) rem (ζ))(k)| 2 (mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 dk ≤ F(dE 3 (ζ) + dE (1) rem (ζ))(k) 2 L ∞ × k0 -k0 k 2 (mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 dk q 5 3 .
As for the high frequency part, we notice that kF(dE (1) rem (ζ)) L 2 q and kF(dE

3 (ζ)) L 2 ζ L ∞ ∂ x ζ L 2 + ∂ x ζ L 2 ∂ 2 x ζ H 1 + ζ L ∞ ∂ 3 x ζ L 2 q 7 6 + q 1 2 ∂ x ζ L 2 ,
where we used (5.15). It follows that

R\[-k0,k0] k 2 |F(dE 3 (ζ) + dE (1) rem (ζ))(k)| 2 (cq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 dk ≤ kF(dE 3 (ζ)+dE (1) rem (ζ)) 2 L 2 × sup k∈R\[-k0,k0] 1 (mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 q 7 3 +q ∂ x ζ 2 L 2 .
Combining the high and low frequency estimates into (5.21) gives us

∂ x ζ 2 L 2 q 5 3 + q ∂ x ζ 2
L 2 , hence, for q 0 sufficiently small we get (5.16).

We conclude with the proof of (5.17). By (5.19) we have

∂ 2 x ζ 2 L 2 ≤ 1 4 R k 4 |F(dE 3 (ζ) + dE (1) rem
(ζ))| 2 (mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 dk. (5.22)
Now we remark that by the Cauchy-Schwarz inequality and integration by parts, one has

∀j ≥ 2, ∂ j x ζ 2 L 2 ≤ ∂ x ζ L 2 ∂ 2j-1 x ζ L 2 q 4 3 .
Since ζ L ∞ and ∂ x ζ L 2 satisfy similar estimates by (5.15) and (5.16), we have

k 2 F(dE 3 (ζ)) L 2 ζ L ∞ ∂ 2 x ζ L 2 + ∂ x ζ H 1 ∂ x ζ L 2 + ζ L ∞ ∂ 4 x ζ L 2 + ∂ x ζ H 1 ∂ 3 x ζ L 2 + ∂ 2 x ζ H 1 ∂ 2 x ζ L 2 q 4 3 , and
kF(dE 3 (ζ)) L ∞ ζ L 2 ∂ x ζ L 2 + ∂ x ζ L 2 ∂ 2 x ζ L 2 + ζ L 2 ∂ 3 x ζ L 2 q 8 6 + q 1 2 ∂ 3 x ζ L 2 q 8 6 + q 1 2 ∂ 2 x ζ 1 2 L 2 ∂ 4 x ζ 1 2 L 2 q 8 6 + q 5 6 ∂ 2 x ζ 1 2 L 2 q 8 6 + q 1 3 ∂ 2 x ζ L 2 , in addition to kF(E (1) rem (ζ)) L ∞ + k 2 F(E (1) rem (ζ)) L 2 q 3 2 .
Thus, proceeding as above, we find k0 -k0

k 4 |F(dE 3 (ζ) + dE (1) rem (ζ))(k)| 2 (mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 dk ≤ kF(dE 3 (ζ)+dE (1) rem (ζ))(k) 2 L ∞ × k0 -k0 k 2 (mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 dk q 7 3 +q 1 3 ∂ 2 x ζ 2 L 2 and R\[-k0,k0] k 4 |F(dE 3 (ζ) + dE (1) rem (ζ))(k)| 2 (mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 dk ≤ k 2 F(dE 3 (ζ) + dE (1) rem (ζ)) 2 L 2 × sup k∈R\[-k0,k0] 1 (mq 2 3 + 1 3 (γ(kF 1 (k)) 2 + δ -1 (kF 2 (k)) 2 )) 2 q 8 3 .
Plugging these estimates into (5.22) and restricting q ∈ (0, q 0 ) if necessary yields (5.17), and the proof is complete.

Convergence results; proof of Theorem 1.7

We are now in position to relate the minimizers of E in D q,R with the corresponding solution of the KdV equation. We first compare

I KdV = inf{E KdV (ξ) : ξ ∈ U 1 } and I q = inf{E(ζ) : ζ ∈ V q,R }.
Lemma 5.4. There exists q 0 > 0 such that the quantities I q and I KdV satisfy

I q = q + E KdV (ζ) + O(q 2 ), uniformly over minimizers of E in V q,R , (5.23)
I q = q + q 5 3 I KdV + O(q 2 ) = q - 3 5 α 0 q 5 3 + O(q 2 ). (5.24)
uniformly over q ∈ (0, q 0 ).

Proof. Recall that, from Lemma 2.7, one has for any ζ ∈ H 2 ,

E(ζ) = (γ + δ) ζ 2 L 2 + E KdV (ζ) + E rem (ζ), (5.25)
with

|E rem (ζ)| ≤ C(h -1 0 , ζ H 1 )( ζ 2 L ∞ ζ 2 L 2 + ζ L ∞ ∂ x ζ 2 L 2 + ∂ 2 x ζ L 2 ∂ x ζ L 2 ). Let ζ be a minimizer of E in V q,R
and note that ζ ∈ H 2 by Lemma 5.1. Using Lemma 5.3 and (5.25), we obtain

I q = E(ζ) = q + E KdV (ζ) + O(q 2 ).
Introducing ξ = S -1 KdV (ζ), we find that ξ ∈ U 1 and

E KdV (ζ) = q 5 3 E KdV (ξ) ≥ q 5 3 I KdV .
Thus we found I q ≥ q + q 5 3 I KdV + O(q 2 ).

Similarly, notice that ζ = S KdV (ξ KdV ) satisfies ζ ∈ V q,R (for q sufficiently small) and, by (5.25)

I q ≤ E( ζ) = q + E KdV ( ζ) + O(q 2 ). Since E KdV ( ζ) = q 5 3 E KdV (ξ KdV ) = q 5
3 I KdV , we deduce

I q ≤ q + q 5 3 I KdV + O(q 2 ).
We have thus proved (5.24).

This next result is the first part of Theorem 1.7, which relates the minimizers of E in V q,R with the minimizers of E KdV in U 1 .

Theorem 5.5. Let q 0 > 0 be such that Theorem 1.6 and Lemma 5.4 hold. Then for any q ∈ (0, q 0 ) and ζ ∈ D q,R , there exists x ζ ∈ R such that

q -2 3 ζ(q -1 3 •) -ξ KdV (• -x ζ ) H 1 q 1 6 ,
uniformly with respect to q ∈ (0, q 0 ) and ζ ∈ D q,R .

Proof. Assume that there exists > 0 and a sequence ζ n ∈ D qn,R with q n 0 such that

∀n ∈ N, inf x0∈R q -2 3 n ζ n (q -1 3 n •) -ξ KdV (• -x 0 ) H 1 ≥ . (5.26) Denote for simplicity ξ n (x) = q -2 3 n ζ n (q -1 3 n x)
. From (5.23) in Lemma 5.4, we have

I qn = E(ζ n ) = q n + E KdV (ζ n ) + O(q 2 n ) = q n + q 5 3 n E KdV (ξ n ) + O(q 2 n ).
By (5.24) in Lemma 5.4, we deduce that

E KdV (ξ n ) -I KdV = O(q 1 3 n ).
In particular {ξ n } n∈N is a minimizing sequence for E KdV satisfying the constraint (γ +δ) ξ n 2 L 2 = 1. It follows [START_REF] Pava | Nonlinear dispersive equations. Existence and stability of solitary and periodic travelling wave solutions[END_REF] that there exists a sequence {x n } n∈N such that

ξ n (• -x n ) -ξ KdV H 1 → 0,
which contradicts (5.26).

The quantitative estimate follows from the argument in [START_REF] Bona | Stability and instability of solitary waves of Korteweg-de Vries type[END_REF]. From the above, we may apply Lemma 4.1 therein and define uniquely x ζ such that ξ, ξ KdV (• -x ζ ) = 0, where we denote ξ = q -23 ζ(q -1 3 •). Following the above estimates and [7, Lemma 5.2], we find

ξ -ξ KdV (• -x ζ ) 2 H 1 E KdV (ξ) -E KdV (ξ KdV ) q 1 3 .
This concludes the proof.

Next we prove the second part of Theorem 1.7, which relates the Lagrange multipliers α with the one of the KdV equation, α 0 .

Theorem 5.6. The number α, defined in Theorem 1.6, satisfies 6 ) as q 0. This implies that 6 ) as q 0 and therefore dE(ζ), ζ = 2q + q 6 Numerical study

α + 1 = q 2 3 α 0 + O(q 5
(ζ) -ξ KdV (• -x ζ ) H 1 = O(q 1
dE KdV (S -1 KdV (ζ)), S -1 KdV ζ -dE KdV (ξ KdV ), ξ KdV = O(q 1
In this section, we provide numerical illustrations of our results as well as some numerical experiments for situations which are not covered by our results. We first describe our numerical scheme, before discussing the outcome of these simulations.

Description of the numerical scheme Our numerical scheme computes solutions for (1.8) for given value of c (and hence does not not follow the minimization strategy developed in this work).

Because we seek smooth localized solutions and our operators involve Fourier multipliers, it is very natural to discretize the problem through spectral methods [START_REF] Trefethen | Software, Environments, and Tools[END_REF]. We are thus left with the problem of finding a root for a nonlinear function defined in a finite (but large) dimensional space. To this aim, we employ the Matlab script fsolve which implements the so-called trust-region dogleg algorithm [START_REF] Conn | Trust-region methods[END_REF] based on Newton's method. For an efficient and successful outcome of the method, it is important to have a fairly precise initial guess. To this aim, we use the exact solution of the Green-Naghdi model, which is either explicit (in the one-layer situation [START_REF] Serre | Contribution à l'étude des écoulements permanents et variables dans les canaux[END_REF]) or obtained as the solution of an ordinary differential equation (in the bi-layer situation [START_REF] Miyata | Long internal waves of large amplitude[END_REF][START_REF] Choi | Fully nonlinear internal waves in a two-fluid system[END_REF]) that we solve numerically. Our solutions are compared with the corresponding ones of the full Euler system. To compute the latter, we use the Matlab script developed by Per-Olav Rusås and documented in [START_REF] Grue | Properties of large-amplitude internal waves[END_REF] in the bilayer configuration while in the one-layer case, the Matlab script of Clamond and Dutykh [START_REF] Clamond | Fast accurate computation of the fully nonlinear solitary surface gravity waves[END_REF] offer faster and more accurate results (although limited to relatively small velocities).

Two-layer setting The solitary-wave solutions of the Miyata-Choi-Camassa system have been studied in the original papers of [START_REF] Miyata | Long internal waves of large amplitude[END_REF][START_REF] Choi | Fully nonlinear internal waves in a two-fluid system[END_REF]. In particular we know that for a given amplitude, or a given velocity, there exists at most one solitary wave (up to spatial translations). The solitary waves are of elevation if δ 2 -γ > 0, of depression if δ 2 -γ < 0, and do not exist if δ 2 = γ. Contrarily to the one-layer situation, the bilayer Green-Naghdi model admits solitary waves only for a finite range of velocities (resp. amplitudes), c ∈ (1, c max (γ, δ)) (resp. |a| ∈ (0, a max (γ, δ))). With our choice of parameters (namely γ = 1, δ = 1/2), one has c max = 1 + 1/8 ≈ 1.06066 and |a max | = 1/2.

As the velocity approaches c max , the solitary waves broadens and its mass keeps increasing. These type of profiles or often refered to as "table-top" profiles, and lead to bore profiles in the limit c → c max . When the velocity is small the numerically computed solitary wave solutions of the bilayer original (F i = 1) and full dispersion (F i = F imp i ) Green-Naghdi systems and the one of the water waves systems (and to a lesser extent the KdV model) agree, so that the curves corresponding to the three former models are indistinguishable in see Figure 2(a). For larger velocities, as in Figure 2(b), the numerically computed solitary wave solutions of the Green-Naghdi and water waves systems is very different from the sech 2 profile of the solitary wave solution to the Korteweg-de Vries equation. It is interesting to see that both the original and full dispersion Green-Naghdi models offer good approximations, even in this "large velocity" limit (the normalized l 2 difference of the computed solutions is ≈ 2.10 -3 in both cases). This means that the internal solitary wave keeps a long-wave feature even for large velocities. These observations were already documented and corroborated by laboratory experiments in [START_REF] Michallet | Experimental study of interfacial solitary waves[END_REF][START_REF] Grue | Properties of large-amplitude internal waves[END_REF][START_REF] Camassa | On the realm of validity of strongly nonlinear asymptotic approximations for internal waves[END_REF].

One-layer setting In the one-layer setting, the script by Clamond and Dutykh [START_REF] Clamond | Fast accurate computation of the fully nonlinear solitary surface gravity waves[END_REF] allows to have a very precise numerical computation of the water waves solitary solution, from which the numerical solutions of the Green-Naghdi models can be compared. In this setting, namely γ = 0 and δ = 1, we have an explicit solution for the Green-Naghdi model [START_REF] Serre | Contribution à l'étude des écoulements permanents et variables dans les canaux[END_REF]:

ζ GN (x) = (c 2 -1) sech 2 ( 1 2 3 c 2 -1 c 2 x) = c 2 ζ KdV (x).
In Figure 3, we compute the solitary waves for our models with different (small) values of the velocity, rescaled by S -1 KdV . One clearly sees, as predicted by Theorem 1.7 and the above formula, that the solitary waves converge towards ξ KdV after rescaling, as c 1. One also sees that the water waves solution is closer to the one predicted by the model with full dispersion than the original Green-Naghdi model. Figure 4 shows that the convergence rate is indeed quadratic for the full dispersion model whereas it is only linear for the original Green-Naghdi model (and therefore only qualitatively better than the KdV model).

A Paradifferential calculus

The definitions and properties below are collected from [START_REF] Métivier | Para-differential calculus and applications to the Cauchy problem for nonlinear systems[END_REF]; see also [START_REF] Chemin | Oxford Lecture Series in Mathematics and its Applications[END_REF][START_REF] Benzoni-Gavage | Multidimensional hyperbolic partial differential equations. First-order systems and applications[END_REF] for relevant references. where ψ(D, ξ) is the Fourier multiplier associated with ψ(η, ξ) (here, ξ is a parameter). The operator is defined for u ∈ H s (R) by density and continuous linear extension.

Figure 1 :

 1 Figure 1: Sketch of the domain and notations in the one-layer and bilayer situations.

Definition 1 . 2 (Proposition 1 . 3 .

 1213 Strongly admissible class of Fourier multipliers). An admissible Fourier multipler F in the sense of Definition 1.1 is strongly admissible if F ∈ C ∞ (R) and for each j ∈ N there exists a constant C j such that |∂ j k F(k)| ≤ C j (1 + |k|) -θ-j . Notice the following. The two aforementioned examples, namely F id i and F imp i are strongly admissible, and satisfy Definition 1.1,iv with (respectively) θ = 0, 1/2.

[ 5 ,

 5 Prop. C.12]. Combining the above estimates together with similar calculations for E yields the desired estimate for |E rem |. The estimate for | dE rem (ζ), ζ | follows in the same way when decomposing

( 1 )

 1 rem (ζ)|. The estimate for |E(1) rem (ζ)| is obtained identically, and the one for |E

( 2 )

 2 rem (ζ)| are obtained using similar estimates when decomposing dE(ζ), ζ given in (2.2).

( 1 )

 1 rem (ζ n ) and by the triangle and Cauchy-Schwarz inequalities, |E

Theorem 4 . 4 (•

 44 Concentration-compactness). Any sequence {e n } n∈N ⊂ L 1 (R) of non-negative functions such that lim n→∞ R e n dx = I > 0 admits a subsequence, denoted again {e n } n∈N , for which one of the following phenomena occurs. • (Vanishing) For each r > 0, one has lim (Dichotomy) There are real sequences {x n } n∈N , {M n } n∈N , {N n } n∈N ⊂ R and I * ∈ (0, I) such that M n , N n → ∞, M n /N n → 0, and xn+Mn xn-Mn e n dx → I * and xn+Nn xn-Nn e n dx → I * as n → ∞. • (Concentration) There exists a sequence {x n } n∈N ⊂ R with the property that for each > 0, there exists r > 0 with xn+r xn-r e n dx ≥ Ifor all n ∈ N.

3 , 5 3

 35 uniformly over D q,R .Proof. By Lemma 2.7, we havedE(ζ), ζ = 2(γ + δ) R ζ 2 dx + dE KdV (ζ), ζ + dE rem (ζ), ζwhere, using Lemma 5.3, one has dE rem (ζ), ζ q 7 uniformly for minimizers of E in V q,R , and sodE(ζ), ζ = 2q + q dE KdV (S -1 KdV (ζ)), S -1 KdV (ζ) + O(qBy Theorem 5.5 there exists x ζ such that S -1 KdV

5 3 dE

 3 KdV (ξ KdV ), ξ KdV + O(q 116 ). Now recall the Euler-Lagrange equations (1.11) and (5.3), which yield immediately2α(ζ) q = -dE(ζ), ζ , 2α 0 = -dE KdV (ξ KdV ), ξ KdV ,and the result follows.

  Small velocity, c = 1.005.

  Large velocity, c = 1.06065.

Figure 2 :

 2 Figure 2: Comparison of the bilayer Green-Naghdi models and the water waves system (γ = 1, δ = 1/2).

  Rescaled solitary waves for c = 1.025, 1.01, 1.002.

Figure 3 :

 3 Figure 3: Comparison of the solutions of the KdV and Green-Naghdi models and the water waves system in the one-layer setting (γ = 0, δ = 1).

10

 10

Figure 4 :

 4 Figure 4: Convergence rate. Log-log plot of the normalized l 2 norm of the error as a function of c -1.

Definition A. 1 (

 1 Symbols). Given m ∈ R and r ≥ 0, we denote Γ m r the space of distributions a(x, ξ) on R 2 such that for almost anyx ∈ R, ξ → a(x, ξ) ∈ C ∞ (R), and ∀α ∈ N, ∃C α > 0 such that ∀ξ ∈ R, ∂ α ξ a(•, ξ) W r,∞ ≤ C α (1 + |ξ|) m-α ,where W r,∞ denote the Hölder space (Lipschitz for integer values).Below, we use an admissible cut-off function ψ in the sense of[START_REF] Métivier | Para-differential calculus and applications to the Cauchy problem for nonlinear systems[END_REF] Definition 5.1.4] and define paradifferential operators as follows (the constant factor depends on the choice of convention for the Fourier transform).

Definition A. 2 (

 2 Paradifferential operators). For a ∈ Γ m 0 and u ∈ S(R), we defineT a u(x) ), e ix• ψ(D, •)a(x, •) (S(R),S (R)) ,

  Similarly, the operators ∂ x F i extend to operators from S (R) to itself, and maps smoothly H s

	u	2 H s P	=	P/2 -P/2	uΛ 2s u dx
	and Λ m is an isomorphism from H s P to H s-m	P into H s-1+θ P	, acting on the
	Fourier coefficients by pointwise multiplication:

P

for any s, m ∈ R.

  ). The argument is standard; see e.g. [42, §I.1,I.2]. Since E P, (ζ) → ∞ as ζ H ν P 2R, any minimizing sequence is bounded, and therefore weakly convergent (up to a subsequence) in the (reflexive) Hilbert space H ν P . We only need to show that for any ζ n ∈ V P,q,2R such that ζ n ζ weakly in H ν P , one has 0 ≤ E P, (ζ) ≤ lim inf Notice first that by the weak lower semi-continuity of • H ν

n→∞ E P, (ζ n ). P , and since is non-decreasing, one has

  Using the decompositions in Lemma 2.8 (changing the domain of integration to [-P/2, P/2]) yieldsdE P (ζ P ), ζ P = 2E2,P (ζ P ) + 3E 3,P (ζ P ) + E (ζ P ) + E 3,P (ζ P ) + E

	This equation is well-defined in (H ν P ) and testing against ζ P ∈ H ν P yields
	2 ( ζ P	2 H ν P	) ζ P	2 H ν P	+ 2α P (γ + δ) ζ P	2 L 2 P	+ dE P (ζ P ), ζ P = 0.
											(2) rem,P (ζ P )
				= 2E P (2) rem,P (ζ P ) -2E	(1) rem,P (ζ P ),
	so that one obtains the identity
	-α P q = E P (ζ P ) +	1 2	E 3,P (ζ P ) +	1 2	E rem,P (ζ P ) -E (2)	(1) rem,P (ζ P ) + ( ζ P	2 H ν P	) ζ P	2 H ν P	.
			2 ( ζ P	2 H ν

P

)Λ 2ν ζ P + dE P (ζ P ) + 2α P (γ + δ)ζ P = 0.

  and recalling α ∈ (-3/2, -1/2). Using Lemma 2.3, ii and Lemma 2.4,i, one easily checks that R

Vincent Duchêne, Dag Nilsson and Erik WahlénJune 27, 2017

).

Acknowledgements. V. Duchêne was partially supported by the Agence Nationale de la Recherche (project ANR-13-BS01-0003-01 DYFICOLTI). D. Nilsson and E. Wahlén were supported by the Swedish Research Council (grant no. 621-2012-3753).

we can estimate R (1 -η (1) n

n (η

→ 0 by Lemma 2.3 i and Lemma 2.1,i, and using η

as the remainder term is bounded by a constant times Mn≤|x-xn|≤Nn e n dx. An analogous argument for E reveals that

and by similar reasoning one finds that

We next claim that q * > 0. Indeed, if q * = 0, we set

.

By (4.10) and since q * = 0, one has c n → 1. Thus we note

n ) → I q -I * < I q as n → ∞. We obtain a similar contradiction involving η [START_REF] Pava | Nonlinear dispersive equations. Existence and stability of solitary and periodic travelling wave solutions[END_REF] n and (4.9) if we assume that q * = q. Hence, 0 < q * < q. In view of the above, we can rescale η (1)

(γ + δ)

n and η (2) n def = (q -q * ) 1 2

(γ + δ)

The following lemma is a direct application of the above definitions [START_REF] Métivier | Para-differential calculus and applications to the Cauchy problem for nonlinear systems[END_REF]Theorem 5.1.15].

Lemma A.3. For any r ≥ 0 and a ∈ Γ m r ⊂ Γ m 0 , and for all s ∈ R, the operator T a extends in a unique way to a bounded operator from H s+m to H s .

If a(ξ) is a symbol independent of x, then T a = a(D), the corresponding Fourier multiplier.

The main tool we use is the following composition property [35, Theorem 6. and T a T b -T ab is a bounded operator from H s+m+m -r to H s , for any s ∈ R.

Of particular interest is the case when the symbol a(x) ∈ L ∞ is independent of ξ. The admissible cut-off function can be constructed so that the paraproduct T a u corresponds to a standard Littlewood-Paley decomposition of the product au. This allows to show that au-T a u is a smoothing operator provided that a is sufficiently regular.

Lemma A.5. Let v ∈ H s and u ∈ H t , and r ≥ 0. Then uv-T v u ∈ H r provided that s+t ≥ 0, s ≥ r and s + t > r + 1/2.

The definitions of the paraproduct in [START_REF] Chemin | Oxford Lecture Series in Mathematics and its Applications[END_REF] and [START_REF] Métivier | Para-differential calculus and applications to the Cauchy problem for nonlinear systems[END_REF] differ slightly but it is not hard to show that [13, Theorem 2.4.1] still holds for the paraproduct as it is defined in [START_REF] Métivier | Para-differential calculus and applications to the Cauchy problem for nonlinear systems[END_REF], and Lemma A.5 follows directly from this theorem.

We conclude with the following lemma, displayed in [START_REF] Métivier | Para-differential calculus and applications to the Cauchy problem for nonlinear systems[END_REF]Theorem 5.2.4] Lemma A.6. Let G ∈ C ∞ (R) be such that G(0) = 0. If u ∈ H s with s > 1/2, then G(u)-T G (u) u ∈ H 2s-1/2 .