Hostname: page-component-54dcc4c588-mz6gc Total loading time: 0 Render date: 2025-09-12T10:53:31.338Z Has data issue: false hasContentIssue false

Asymptotic attraction with algebraic rates toward fronts of dispersive–diffusive Burgers equations

Published online by Cambridge University Press:  29 August 2025

Milena Stanislavova
Affiliation:
Department of Mathematics, University of Alabama – Birmingham, Birmingham, AL, USA
Atanas G. Stefanov*
Affiliation:
Department of Mathematics, University of Alabama – Birmingham, Birmingham, AL, USA
*
Corresponding author: Atanas G. Stefanov; Email: stefanov@uab.edu
Rights & Permissions [Opens in a new window]

Abstract

Burgers equation is a classic model, which arises in numerous applications. At its very core, it is a simple conservation law, which serves as a toy model for various dynamics phenomena. In particular, it supports explicit heteroclinic solutions, both fronts and backs. Their stability has been studied in detail. There has been substantial interest in considering dispersive and/or diffusive modifications, which present novel dynamical paradigms in such a simple setting. More specifically, the KdV–Burgers model has been shown to support unique fronts (not all of them monotone!) with fixed values at $\pm \infty$. Many articles, among which [11], [9] and [10], have studied the question of stability of monotone (or close to monotone) fronts. In a breakthrough paper [2], the authors have extended these results in several different directions. They have considered a wider range of models. The fronts do not need to be monotone but are subject of a spectral condition instead. Most importantly, the method allows for large perturbations, as long as the heteroclinic conditions at $\pm \infty$ are met. That is, there is asymptotic attraction to the said fronts or equivalently the limit set consists of one point. The purpose of this paper is to extend the results of [2] by providing explicit algebraic rates of convergence as $t\to \infty$. We bootstrap these results from the results in [2] using additional energy estimates for two important examples, namely KdV–Burgers and the fractional Burgers problem. These rates are likely not optimal, but we conjecture that they are algebraic nonetheless.

Information

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2025. Published by Cambridge University Press.

1. Introduction

The Burger’s equation

(1.1)\begin{equation} u_t - u_{xx} +u u_x=0, u(0,x)=u_0(x) \end{equation}

where $u:\mathbb R_+\times \mathbb R \to \mathbb R$, is a particularly simple example of a conservation law in one spatial dimension. As a simple model, which resembles (both in its form as well as its behaviour) many distinguished equations in mathematical physics, it is an ubiquitous object in the modern dynamical systems theory. The Cauchy problem, along with many other properties, has been explored in the literature. In particular, and specifically for this model, the Cole–Hopf transformation has played an important role as it transforms (1.1) into a linear problem. More specifically, the assignment $u=-2\frac{w_x}{w}$ reduces (1.1) to the linear heat equation for w, namely

\begin{equation*}w_t-w_{xx}=0, w(0,x)=e^{-\frac{1}{2}\int_0^x u_0(y) dy}\end{equation*}

In this way, one can write explicit solutions and derive particular properties of the conservation laws, which are not transparent in the usual general implicit formulations.

Along with the Burgers equation, related models have been analysed. The KdV–Burgers model

(1.2)\begin{equation} u_t-u_{xx}+u u_x= \nu u_{xxx}, \end{equation}

which describes, among other things, bores in viscous fluids [Reference Su and Gardner12], and weak plasma shocks propagating transversely to a magnetic field. Many similar models were analysed recently, among them the following general model [Reference Barker, Bronski, Hur and Yang2],

(1.3)\begin{equation} u_t - u_{xx}+u u_x = {\mathcal L} u. \end{equation}

Here, ${\mathcal L}$ is a multiplier operator, given by $\widehat{{\mathcal L} f}(k) = l(k) \hat{f}(k)$. Now, we must impose assumptions on ${\mathcal L}$. Namely, it is required that ${\mathcal L}$ maps real-valued functions into real-valued functions, and also its quadratic form is non-positive. That is, it satisfies

(1.4)\begin{eqnarray} & & \langle {\mathcal L} f,f \rangle =\Re \int_{-\infty}^\infty l(k) |\hat{f}(k)|^2 dk\leq 0 \end{eqnarray}
(1.5)\begin{eqnarray} & & {\mathcal L} [1]=0 \end{eqnarray}

The condition (1.5) can be alternatively stated as $l[0]=0$ and l is continuous at zero, but the form (1.5) is more suitable in this context. Note that (1.3) incorporates (1.2), but includes much more. In fact, we list several physically relevant cases: the KdV–Burgers, ${\mathcal L}=\nu \partial_{xxx}$, the Benjamin–Ono–Burgers, ${\mathcal L}=\partial_x |\partial_x|$, the Hilbert–Burgers case ${\mathcal L}= \mathcal H=\partial_x|\partial_x|^{-1}$, the fractional Burgers Footnote equation 1 with ${\mathcal L}=-\sum_{j=1}^N a_j (-\partial_{xx})^{\alpha_j}, \ \ \ 0 \lt \alpha_1 \lt \alpha_2 \lt \ldots \lt \alpha_N \lt 1, a_j \gt 0$, etc. All these could be of interest.

Understanding of the long-time dynamics of such models is clearly important, and as it turns out a challenging problem. As is usual in conservation laws, the standard background is starting with heteroclinic data with enough decay towards the limiting values at $\pm \infty$. To this end, assume that the initial data u 0 satisfies

(1.6)\begin{equation} \lim_{x\to \pm \infty} u_0(x)=u_\pm \end{equation}

where $u_- \gt u_+$ are constants.

1.1. Galilean transformation: reductions and existence of fronts

Let us record two elementary properties of Equation (1.3). For the Burgers equation, i.e. the model (1.3) with ${\mathcal L}=0$, we have true invariances with respect to the Galilean transformation

\begin{equation*}u(t,x)=v(t,x- ct)+c,\end{equation*}

as well as the scaling transformation

\begin{equation*}u\to u_\lambda:=\lambda u(\lambda^2 t, \lambda x), \lambda \gt 0\end{equation*}

These two transformations, in general, do not commute with the flow of (1.3), but they transform it into a model of the same type, i.e. it still satisfies (1.4), (1.5). Since our goal is to understand the long-term dynamics of heteroclinic initial data, we can use these transformations to our advantage as follows: given $u_{\pm}, u_- \gt u_+$, introduce $c=\frac{u_- + u_+}{u_--u_+}, \lambda=\frac{u_- - u_+}{2} \gt 0$ and a function v, so that

\begin{equation*} u(t,x)=\lambda U(\lambda^2 t, \lambda x - c\lambda^2 t) + c \lambda \end{equation*}

Then, the function U will satisfy the Footnote equation 2

(1.7)\begin{equation} U_t - U_{xx}+U U_x={\mathcal L}_1 U, \end{equation}

where $\widehat{{\mathcal L}_1 g}(k)=\lambda^{-2} l(\lambda k) \hat{g}(k)$, $\lim_{x\to \mp \infty} U(x) =\pm 1$. As we shall see, ${\mathcal L}_1$ also satisfies the relevant assumptions, (1.4), (1.5), so it will suffice to work with the model, where $U_{\mp}=\pm 1$. Naturally, and by the choice of the transformations, this corresponds to a case, where c = 0 as well. Thus, we restrict our attention to the steady-state solutions of (1.3) with end values $\mp 1$, which in this case are functions $\phi:\mathbb R\to \mathbb R$,

(1.8)\begin{equation} -\phi_{xx}+\phi \phi_x={\mathcal L} \phi; \ \ \lim_{x\to \mp \infty} \phi(x)=\pm 1. \end{equation}

The existence theory for such a problem is, of course, highly non-trivial. In fact, as in previous papers, we shall make assumptions to this effect. We mention in this regard that for the KdV–Burger’s model, the existence of such monotone ϕ was proved in [Reference Bona, Rajopadhye and Schonbek4, Reference Bona and Schonbek5]. It is worth adding that these fronts are monotone if and only if $|\nu|\leq \frac{1}{4}$, [Reference Bona, Rajopadhye and Schonbek4]. In addition, such solutions are unique (up to a translation), and if one specifies the values $u_{\pm}$. In fact, in the case $\nu=-\frac{1}{10}$, there is the explicit formulaFootnote 3

\begin{equation*} u(t,x)=\frac{6}{5}\left[\textrm{sech}^2\left(x+\frac{12}{5}t\right)- 2\tanh\left(x+\frac{12}{5}t\right)-2\right]. \end{equation*}

See [Reference Jeffrey and Mohamad7].

1.2. Asymptotic attraction towards heteroclinic fronts

We now review the literature regarding earlier asymptotic stability results. Pego [Reference Pego11] has shown that the monotonic fronts are orbitally asymptotically stable. More specifically, he showed that any (small) perturbation of the front converges to some translate of the front. In later works, Khodja [Reference Khodja8] and Naumkin–Shishmarev [Reference Naumkin and Shishmarev9, Reference Naumkin and Shishmarev10] showed that the KdV–Burgers fronts are asymptotically stable, even for some “slightly” non-monotonic fronts, i.e. $\nu: 0 \lt \nu-\frac{1}{4} \lt \lt 1$.

Next, we introduce the results of Barker–Bronski–Hur–Yang [Reference Barker, Bronski, Hur and Yang2], which represents a real breakthrough in the study of this problem. Specifically, the authors improved the aforementioned results in at least three important directions. Firstly, they consider all (non-necessarily monotonic) fronts, subject to a natural and specific spectral condition,Footnote 4 see Theorem 1.1. Secondly, they extend the range of applicability to a wide range of examples, driven by possibly a mixture of dissipative and/or dispersive linear operators, subject to (1.4). Thirdly, and most impressively, they allow large data perturbations for asymptotic attraction type of result vis-a-vis asymptotic stability. That is, as long as the initial data approaches the fixed heteroclinics with good decay properties, and the fronts exist (with the prescribed spectral condition), the solution will stay close to translates (defined dynamically) of the front, all in L 2 norm.

We state the precise result below.

Theorem 1.1 (Barker–Bronski–Hur–Yang [Reference Barker, Bronski, Hur and Yang2])

Suppose that ${\mathcal L}$ satisfies the dispersivity property (1.4), while the front ϕ satisfies the profile equation (1.8) and in addition.

(1.9)\begin{equation} \lim_{x\to \mp \infty} \phi(x)=\pm 1, \phi_x, \phi_{xx}\in L^2(\mathbb R), \int_{\mathbb R} (1+|x|) |\phi_x(x)| dx \lt \infty. \end{equation}

Finally, assume that for some $\epsilon\in (0,1)$, the Schrödinger operator

(1.10)\begin{equation} \mathcal H_\epsilon:=-(1-\epsilon) \partial_{x}^2+\frac{1}{2} \phi_x \ \ {has\ exactly\ one\ negative\ eigenvalue}. \end{equation}

Assume that u 0 is a heteroclinic data, so that $ \lim_{x\to \mp \infty} u_0(x)=\pm 1$ and more precisely, $u_0(x)-\phi\in L^2(\mathbb R)$. Then, there exists a continuous function $x_0(t): \mathbb R_+\to \mathbb R$, so that the solution to (1.3) can be written in the form

(1.11)\begin{equation} u(t, x)=\phi(x-x_0(t))+v(t, x-x_0(t)), \end{equation}

so that $t\to \|v(t, \cdot)\|_{L^2(\mathbb R)}$ is a monotonically decreasing function of time and

(1.12)\begin{equation} \lim_{t\to \infty} \|v(t, \cdot)\|_{L^p(\mathbb R)}=0, p\in (2, \infty]. \end{equation}

Let us record an energy estimate, which came up in the proof of Theorem 1.1. More precisely, and under the assumptions of Theorem 1.1, and specifically (1.10), we have

(1.13)\begin{equation} \partial_t \|v\|_{L^2}^2\leq -C \|v_x\|_{L^2}^2. \end{equation}

We shall make extensive use of (1.13) in the sequel, but note right away that the monotonicity of the mapping $t\to \|v(t, \cdot)\|_{L^2(\mathbb R)}$ is a direct consequence of it.

Next, a few remarks are in order.

  1. (1) Note that the results in Theorem 1.1 do not require smallness of the perturbation $u-\phi$. In other words, this is a strong asymptotic attraction result rather than a mere asymptotic stability result.

  2. (2) Related to the previous point, this implies a rather strong uniqueness (up to the usual translation) theorem for such fronts ϕ. Namely, fronts ϕ with the properties (1.9), (1.10), if any, are unique. The uniqueness result was certainly explored in the literature in certain contexts (see for example [Reference Jeffrey and Mohamad7] for the KdV–Burgers case), but not in the generality of the model (1.3).

  3. (3) In the statement of Theorem 1.1, the particular values of the limits $\phi_{\pm}=\lim_{x\to \pm\infty} \phi(x)$, are set to be $\phi_{\pm}=\mp 1$ only out of convenience. Assuming that solutions exist for some other pair , subject to $\phi_{\pm}: \phi_- \gt \phi_+$, one can retrace the proof in [Reference Barker, Bronski, Hur and Yang2] to the same results, as stated in Theorem 1.1.

  4. (4) Note that since $\int_{-\infty}^{+\infty}\phi_x(x) dx=-2$, the Schrödinger operators $\mathcal H_\epsilon$ do have negative spectrum, which consists of eigenvalues. The condition (1.10) requires that the spectrum is minimal, in the sense that it consists of only one eigenvalue.

  5. (5) Related to the previous item, we note that the spectral condition (1.10) has been rigorously verified in [Reference Barker, Bronski, Hur and Yang2], for the (unique) monotone KdV–Burger’s fronts, i.e. for the case $|\nu| \lt \frac{1}{4}$. In the companion paper [Reference Barker, Bronski, Hur and Yang3], the authors further explored the validity of (1.10) beyond the monotone cases. In fact, informal numerics show that (1.10) holds all the way up to $|\nu| \lt 4.1$, while rigorous numerics confirm this up to $|\nu| \lt 3.9$.

1.3. Main result: pointwise decay for the KdV–Burgers equation

Our first result concerns the standard KdV–Burgers model. We bootstrap upon the results of [Reference Barker, Bronski, Hur and Yang2] to obtain pointwise in time bounds for various Lebesgue norms. The precise statement is below.

Theorem 1.2 (Pointwise decay to the fronts for the KdV–Burgers model)

For the KdV–Burgers model, i.e. ${\mathcal L}=\nu \partial_{xxx}$, under the assumptions of Theorem 1.1 and

(1.14)\begin{equation} \int_{\mathbb R} |u_0(x)-\phi(x)|^2 (1+|x|) dx \lt \infty, \end{equation}

the decomposition (1.11) still holds, with

(1.15)\begin{equation} \|v(t, \cdot)\|_{L^2}\leq \frac{C}{\sqrt{t}}. \end{equation}

In addition, for every $q\in (1,2)$, there is a decay $\lim_{t\to \infty} \|v(t, \cdot)\|_{L^q}=0$, and in fact for every δ > 0, there exists $C=C_\delta$, so that

(1.16)\begin{eqnarray} \|v(t, \cdot)\|_{L^q}\leq C_\delta t^{-(1-\frac{1}{q})+\delta}, 1 \lt q \lt 2. \end{eqnarray}

For p > 2, the bound is

(1.17)\begin{equation} \|v(t, \cdot)\|_{L^p}\leq C t^{-\frac{1}{p}}, 2\leq p \lt \infty. \end{equation}

Some important remarks to be made herein.

  1. (1) While we do require a more localized data, see (1.14), we still allow such perturbations to be large.

  2. (2) The decay rate (1.15) is faster, compared to the linear heat equation, where we only have $\|v(t)\|_{L^2(\mathbb R)}\sim t^{-\frac{1}{4}}$. This enhanced decay rate can be thought of as an anomalous dissipation of sorts, caused by the non-trivial background around the front solution and the lack of edge resonance.

  3. (3) The estimate (1.17) is obtained as an interpolation between the L 2 bounds (1.15) and the uniform $L^\infty$ bound (1.12), obtained earlier in [Reference Barker, Bronski, Hur and Yang2].

1.4. Main result: pointwise decay for the fractional Burgers equation

For clarity, we consider the fractional Burgers equation to be the evolution problem

(1.18)\begin{equation} u_t - u_{xx}+\sum_{j=1}^N a_j (-\partial_{xx})^{\alpha_j} u + u u_x=0, \end{equation}

subject to initial data $u(0,x)=u_0(x)$. Clearly, the evolution of (1.18) preserves odd initial data. We have the following, more precise than (1.12) result.

Theorem 1.3. Suppose that ϕ is an odd stationary front of the fractional Burger equation

(1.19)\begin{equation} -\phi_{xx}+\sum_{j=1}^N a_j (-\partial_{xx})^{\alpha_j} \phi +\phi\phi_x=0, \end{equation}

where $0 \lt \alpha_1 \lt \alpha_2 \lt \ldots \lt \alpha_N \lt 1, a_j\geq 0$. In addition, we assume that ϕ satisfies $\phi_x, \phi_{xx} \in L^2(\mathbb R)$ and the spectral assumption (1.10).

Suppose now that $u_0=\phi+v_0$, where $v_0\in L^1\cap L^2(\mathbb R)$ is an odd function. Then, the solution u of (1.3), with initial data u 0 can be written as

\begin{equation*}u(t,x)=\phi(x)+v(t,x),\end{equation*}

where for t > 2, one has the bounds

(1.20)\begin{eqnarray} \|v(t, \cdot)\|_{L^p} &\leq & C \|v_0\|_{L^2\cap L^1} \left(\frac{\ln(t)}{t}\right)^{\frac{1-\frac{1}{p}}{2}}, 1 \lt p \lt 2 \end{eqnarray}
(1.21)\begin{eqnarray} \|v(t, \cdot)\|_{L^q} &\leq & C \left(\frac{\ln(t)}{t}\right)^{\frac{7}{24}-\frac{1}{12 q}}, 2 \lt q\leq \infty \end{eqnarray}
Remarks.
  • Note that as before, the perturbation v 0 may be large, as long as it is well-localized, i.e. $v_0\in L^1\cap L^2(\mathbb R)$.

  • Related to the previous item, it is worth noting that our proof fails under an assumption such as $v_0\in L^q\cap L^2, 1 \lt q$. That is, we cannot establish the algebraic decay rates (for any $p\in (1, \infty]$!) in (1.20), (1.21), without assuming first that $v_0\in L^1(\mathbb R)$. At this point, it is unclear whether this is a technical obstacle or one indeed needs to require L 1 decay of the perturbation.

  • As before, Theorem 1.3 provides a strong uniqueness property of the odd fronts ϕ satisfying (1.19).

  • It is an open question whether or not the algebraic decay rates stated in (1.20) are sharp. While the logarithmic corrections are most certainly not, one might wonder about the other bounds in (1.20), (1.21).

  • The assumption for existence of solution of (1.19), which in addition satisfies $\phi_x, \phi_{xx} \in L^2(\mathbb R)$ and the spectral assumption (1.10) is a non-trivial one. In some earlier works, see for example [Reference Achleitner, Cuesta and Hittmeir1], stationary fronts were constructed for similar fractional KdV–Burger’s problems, but with a non-trivial KdV term. The scheme in [Reference Achleitner, Cuesta and Hittmeir1] can be applied to (1.19). More specifically, assume that we have only one fractional term in (1.19), $(-\partial_{xx})^\alpha=D^{2\alpha}, 0 \lt \alpha \lt \frac{1}{2}$. Then, we look for solution $\phi:\lim_{x\to \pm \infty}\phi(x)=\mp 1$,

    (1.22)\begin{equation} \phi_x =-\partial_{x}[-\partial_{xx}+D^{2\alpha}]^{-1} [\phi_x \phi]=K\ast [\phi_x\phi], \ \ \hat{K}(\xi)= \frac{2\pi i \xi}{4\pi^2 \xi^2 +(2\pi)^{2\alpha} |\xi|^{2\alpha}}. \end{equation}

    It can be shown that $|K(x)|\leq C |x|^{-(2-2\alpha)}, |x| \gt \gt 1$. We conjecture that this equation is solvable and the precise decay rate is inherited by ϕx, due to the integral formulation (1.22). That is $|\phi_x(x)|\leq C |x|^{-(2-2\alpha)}, |x| \gt \gt 1$, and so $\phi_x\in L^1(\mathbb R)$, whereas $|\phi_{xx}(x)|\leq C |x|^{-(3-2\alpha)}, |x| \gt \gt 1$, etc. In such a way, ϕx will never satisfy the assumptionFootnote 5 $(1+|x|)\phi_x(x)\in L^1(\mathbb R)$. On the other hand, the symmetric operator $\mathcal H_\epsilon=-(1-\epsilon)\partial_{xx}+\frac{1}{2}\phi_x$ may still satisfy the assumption (1.10), which is the only one used in the proof of Theorem 1.1.

We plan our paper as follows. In Section 2, we introduce some basic notations, as well as the Sobolev embedding and various related interpolation inequalities. We introduce the notion of a diffusive semigroup, as this allows us for an effective energy estimates for the fractional Burgers evolution in $L^p, p \gt 1$ norms. In Section 3, we prove Theorem 1.2. Specifically, for the case of the KdV–Burgers model, we show enhanced energy estimates for the weighted norms $\left(\int_{\mathbb R} |v(t,x)|^2 |x| dx\right)^{\frac{1}{2}}$, which allows us to bootstrap the results of [Reference Barker, Bronski, Hur and Yang2] to pointwise in time decay estimates for $\|v(t, \cdot)\|_{L^2}$. In Section 4, we prove Theorem 1.3. More specifically, we use the integration by parts techniques of [Reference Chamorro and Lemarié-Rieusset6] to derive energy estimates for $\|v(t, \cdot)\|_{L^p}, p \gt 1$ and ultimately for p = 1. We again use in a critical way the energy estimate (1.13) from [Reference Barker, Bronski, Hur and Yang2], and in conjunction with the a priori estimates for $\|v(t, \cdot)\|_{L^1}$, we conclude pointwise in time decay for the solution norms.

2. Preliminaries

We shall use standard notations for the Lebesgue spaces $L^p, 1\leq p\leq \infty$. More generally, the Sobolev spaces $W^{k,p}(\mathbb R)$ for integer $k, 1 \lt p \lt \infty$, are defined via the norms

\begin{equation*}\|f\|_{W^{k,p}}=\|f\|_{L^p}+ \sum_{l=1}^k \|\partial_x^l f\|_{L^p}\end{equation*}

2.1. Fourier transform, fractional Burger’s and the Bernstein inequality

We take the Fourier transform and its inverse in the form

\begin{equation*} \hat{f}(\xi)=\int_{\mathbb R} f(x) e^{-2\pi i x\xi} dx, \ \ f(x)=\int_{\mathbb R} \hat{f}(\xi) e^{2\pi i x\xi} d\xi. \end{equation*}

In this notation, we have $\widehat{f''}(\xi)=-4\pi^2 \xi^2 \hat{f}(\xi)$. More generally, we introduce the fractional operators $(-\partial_{xx})^\alpha$ as acting on Schwartz functions via the symbol $(2\pi|\xi|)^{2\alpha}$. More precisely,

\begin{equation*} \widehat{(-\partial_{xx})^\alpha f}(\xi)= (2\pi|\xi|)^{2\alpha} \hat{f}(\xi). \end{equation*}

In fact, this allows one to introduce a companion spaces to $W^{s,p}$, where $s \gt 0, 1 \lt p \lt \infty$ is not necessarily an integer. More precisely, these are given by the norm

\begin{equation*} \|f\|_{W^{s,p}}=\|f\|_{L^p}+\| (-\partial_{xx})^{\frac{s}{2}} f\|_{L^p}. \end{equation*}

These spaces are equivalent to the previously defined $W^{k,p}$ for integer k. The Sobolev inequality in one spatial dimension takes the simple form

\begin{equation*} \|u\|_{L^q(\mathbb R)}\leq C \|u\|_{W^{s,p}} \ \ 1 \lt p \lt q \lt \infty, \ \ s\geq \frac{1}{p}-\frac{1}{q}. \end{equation*}

Next, we discuss Bernstein’s inequality, which is a relatively crude, but helpful form of Sobolev embedding. Consider a function $f: {\mathbf R}^n\to {\mathbf R}^n$ and a set of finite Lebesgue measure A. Let $f_A: \widehat{f_A}(\xi)=\hat{f}(\xi) \chi_A(\xi)$, where $\chi_A(x)=\left\{\begin{array}{cc} 1 & x\in A\\ 0 & x\notin A \end{array} \right..$ Then, for all $1\leq p\leq \infty$,

(2.1)\begin{equation} \|f_A\|_{L^q({\mathbf R}^n)}\leq |A|^{\frac{1}{p}-\frac{1}{q}} \|f\|_{L^p({\mathbf R}^n)}. \end{equation}

An easy estimate, which follows immediately from Plancherel’s identity, is

\begin{equation*} \|\partial_x f_A\|_q\leq C \min_{\xi\in A} |\xi| \|f_A\|_q, 1\leq q\leq \infty. \end{equation*}

We will use the rough Fourier cutoff projections

\begin{equation*} f_{ \lt \epsilon}: \widehat{f_{ \lt \epsilon}}(\xi):=\hat{f}(\xi)\chi_{(-\epsilon, \epsilon)}(\xi), \ \ f_{ \gt \epsilon}:=f- f_{ \lt \epsilon}. \end{equation*}

We also need the Gagliardo–Nirenberg estimate (or the log-convexity of the map $p\to \|f\|_p$). More precisely, for every $1\leq p \lt q \lt r\leq \infty: \frac{1}{q}=\frac{\theta}{p}+\frac{1-\theta}{r}$,

\begin{equation*} \|f\|_q\leq \|f\|_p^\theta \|f\|_r^{1-\theta}. \end{equation*}

Finally, there is the estimate, $ \|f\|_{L^\infty}^2 \leq C \|f'\|_{L^2} \|f\|_{L^2}, $ which will also be useful.

2.2. Diffusion semigroups

The notion of diffusive semigroups will play a role in our considerations. Shortly, we say that $\{T(t)\}_{t\geq 0}$ is a C 0 semigroup on a Banach space, if $T(t)\in B(X)$, $T(t+s)=T(t)T(s), T(0)=Id$ and for every $f\in X: \lim_{t\to 0+} \|T(t) f - f\|_X=0$. It is standard that one can always associate with such objects a generator L, a generally unbounded operator, which can be introduced on a domain $D(L)=\{f\in X: \lim_{t\to 0+} \frac{T(t) f - f}{t} \ \ \textrm{exists}\}$ as

\begin{equation*} L f:=\lim_{t\to 0+} \frac{T(t) f - f}{t}, f\in D(L). \end{equation*}

Naturally, $T(t) f\in D(L)$ solves the X valued differential equation $x'=L x, x(0)=f$ for any $f\in X$, so we informally denote it by $T(t)=e^{t L}$.

Definition 2.1. We say that a C 0 semigroup on $L^2({\mathbf R}^n)$, T(t), is a diffusive semigroup, if $T(t), t\geq 0$ is self-adjoint, and given by a probability density kernel, i.e.

\begin{equation*} T(t) f(x) = \int_{{\mathbf R}^n} p_t(x-y) f(y) dy, f\in L^2({\mathbf R}^n), p_t(x)\geq 0, \int_{{\mathbf R}^n} p_t(x) dx = 1. \end{equation*}

We observe that the generator L in such cases is necessarily a non-positive operator, since for all $t\geq 0$ and a test function f, we have

\begin{align*} \langle T(t) f,f \rangle &= \int p_t(x-y) f(y) \bar{f}(x) dx dy \leq \\ & \leq \left(\int p_t(x-y) |f (y)|^2 dx dy\right)^{\frac{1}{2}} \left(\int p_t(x-y) |f (x)|^2 dx dy\right)^{\frac{1}{2}}=\|f\|^2, \end{align*}

whence $\langle L f,f \rangle \leq \limsup_{t\to 0+} \frac{\langle T(t) f - f,f \rangle}{t}\leq 0$.

It is standard that the operator $L=\partial_{xx}$ is a generator of such semigroup. More generally, the fractional Burgers operators $L=-\sum_{j=1}^N a_j (-\partial_{xx})^{\alpha_j}, a_j\geq 0, 0 \lt \alpha_1 \lt \ldots \alpha_N\leq 1$ also generate diffusion semigroups,Footnote 6 see for example [Reference Chamorro and Lemarié-Rieusset6]. An important property of such semigroups is the following.

Proposition 2.2 (Reference Chamorro and Lemarié-Rieusset6, Theorem 3.2)

Assume that L generates a diffusive semigroup. Then, for all $1\leq p \lt \infty$, we have

(2.2)\begin{equation} \int_{{\mathbf R}^n} f|f|^{p-2} L(f) dx\leq 0. \end{equation}

If, in addition, $ 2\leq p \lt \infty$, then

(2.3)\begin{equation} -p \int f|f|^{p-2} L(f) dx \geq \int |\sqrt{-L}(f|f|^{\frac{p}{2}-1})|^2 dx \end{equation}

In particular, we have that for each $0 \lt \alpha\leq 1, 2\leq p \lt \infty$

(2.4)\begin{equation} \int_{{\mathbf R}^n} ((-\Delta)^\alpha f) f|f|^{p-2} dx \geq 0. \end{equation}

3. Attraction to fronts: the Burger-KdV model

We start by revisiting some of the results obtained in the proof of Theorem 1.1. Indeed, in the course of the proof, the authors set $\dot{x}_0(t)=\gamma \int_{-\infty}^\infty \phi_x(x) v(x) dx$, withFootnote 7 γ > 1. This is done to ensure that the rank-one perturbation $\mathcal H_\gamma=-\partial_{xx}+\frac{1}{2}\phi_x+\gamma \langle \phi_x,\cdot \rangle \phi_x$ of the important linearized operatorFootnote 8 $\mathcal H=-\partial_{xx}+\frac{1}{2}\phi_x$ has the property $\mathcal H_\gamma\geq 0$. Two important features of the construction are

(3.1)\begin{eqnarray} & & t\to \|v(t, \cdot)\|_{L^2}^2 \ \ \textrm{is monotonically decreasing} \end{eqnarray}
(3.2)\begin{eqnarray} & & \int_0^t |\dot{x}_0(s)|^2 ds+ \int_0^t \|v_x(s, \cdot)\|_{L^2}^2 ds \leq C \|v_0\|_{L^2}^2 \end{eqnarray}

We make substantial use of these two facts in our subsequent analysis. Next, we write the equation for the residual $v: u(x,t)=\phi(x-x_0(t))+v(t, x-x_0(t))$. This takes the form

(3.3)\begin{equation} v_t - v_{xx} -\dot{x}_0 v_x-\dot{x}_0 \phi_x+\phi_x v+\phi v_x+ v v_x+\nu v_{xxx}=0. \end{equation}

We setup an energy estimate for v in the space $L^2_1(\mathbb R)$. More precisely, we shall show that

\begin{equation*} \int_{-\infty}^\infty |v(t,x)|^2 |x| dx\leq C, t\geq 0. \end{equation*}

Of course, the global solution v produced in [Reference Barker, Bronski, Hur and Yang2] does not a priori live in the space $L^2_1(\mathbb R)$, so one begins by only formally setting the energy estimate for it. We will ignore the details of this justification, as this is standard by considering instead a quantity of the form $\int_{-\infty}^\infty |v(t,x)|^2 |x| \zeta(x/N)dx$ for $N \gt \gt 1$ and some smooth cutoff function ζ. With these caveats, we proceed to bound the $L^2_1(\mathbb R)$ norm. Taking a dot product of (3.3) with $v |x|$, we have

\begin{align*} &\frac{1}{2} \partial_t \int v^2(t,x) |x| dx -\int v_{xx} v |x| dx -\dot{x}_0 \int v_x v |x| dx - \dot{x}_0 \int \phi_x v |x| dx+\int \phi_x v v |x| dx \,+\\ & \quad +\int \phi v_x v |x| dx + \int v v_x v |x| dx +\nu \int v_{xxx} v |x| dx =0. \end{align*}

We now proceed to estimate various terms, by keeping in mind the relation (3.2), which guides us about the acceptable error terms. After integration by parts and algebraic manipulations, we have

\begin{align*} -\int v_{xx} v |x| dx &= \int (v_x)^2 |x| dx + \int_{-\infty}^\infty v_x v sgn(x) dx \geq -v^2(t,0)\geq -\|v(t, \cdot)\|_{L^\infty}^2,\\ \left| \dot{x}_0 \int v_x v |x| dx \right| &= \left| \frac{\dot{x}_0}{2} \int v^2 sgn(x) dx \right| \leq C |\dot{x}_0 (t)| \|v(t, \cdot)\|_{L^2}^2\leq C |\dot{x}_0 (t)| \|v(t, \cdot)\|_{L^2}\\ \left| \dot{x}_0 \int \phi_x v |x| dx\right|&\leq C |\dot{x}_0 (t)| \|v(t, \cdot)\|_{L^2}, \\ \left|\int v v_x v |x| dx\right|&=\frac{1}{3} \left|\int v^3 sgn(x)dx\right| \leq C \|v(t, \cdot)\|_{L^3}^3\leq C \|v\|_{H^{\frac{1}{6}}}^3\leq C \|v_x(t, \cdot)\|_{L^2}^{\frac{1}{2}} \|v(t, \cdot)\|_{L^2}^{\frac{5}{2}}\\ &\leq C \|v_x(t, \cdot)\|_{L^2}^{\frac{1}{2}} \|v(t, \cdot)\|_{L^2}^{\frac{3}{2}} \end{align*}

Also,

\begin{equation*} \int \phi_x v^2 |x| dx\geq -\|v(t, \cdot)\|_{L^\infty}^2 \int |\phi_x(x)| |x| dx =-C \|v(t, \cdot)\|_{L^\infty}^2. \end{equation*}

Here, in addition to Sobolev embedding and Gagliardo–Nirenberg estimates, we used the bound $\|v(t, \cdot)\|_{L^2}\leq C$. Next,

\begin{eqnarray*} & & \int \phi v_x v |x| dx=-\frac{1}{2} \left( \int \phi_x(x) v^2(x) |x| dx + \int \phi(x) v^2(t,x) sgn(x) dx \right)\geq \\ & & \geq - C \|v(t, \cdot)\|_{L^\infty}^2-\frac{1}{2}\int \phi(x) v^2(t,x) sgn(x) dx \end{eqnarray*}

Now, note that $-\phi(x) sgn(x) \gt 0$ for all x and also there exists A, so that $-\phi(x) sgn(x) \gt \frac{1}{2}$ for $|x| \gt A$. This last statement follows from $\lim_{x\to \pm \infty} \phi(x) = \mp 1$. So, we further estimate

\begin{eqnarray*} & & -\frac{1}{2}\int \phi(x) v^2(t,x) sgn(x) dx \gt -\frac{1}{2}\int_{|x| \gt A} \phi(x) v^2(t,x) sgn(x) dx \gt \frac{1}{4} \int_{|x| \gt A} v^2(t,x) dx = \\ & & = \frac{1}{4} \int v^2(t,x) dx -\frac{1}{4} \int_{-A}^A v^2(t,x) dx \geq \frac{1}{4} \int v^2(t,x) dx - C \|v(t, \cdot)\|_{L^\infty}^2, \end{eqnarray*}

Finally, for the term $\int v_x'' v |x| dx$, we have

\begin{equation*} \left|\int v_x'' v |x| dx \right|= \left| \frac{1}{2} \int (v_x)^2 sgn(x) - \int v_x' v sgn(x) dx \right| \leq 2 |v(0)| |v_x(0)| + C \|v_x(t, \cdot)\|_{L^2}^2. \end{equation*}

Putting it all together and estimating $-\|v(t, \cdot)\|_{L^\infty}^2 \geq - C\|v_x(t, \cdot)\|_{L^2} \|v(t, \cdot)\|_{L^2}$, we obtain

\begin{eqnarray*} & & \frac{1}{2} \partial_t \int v^2(t,x) |x| dx + \frac{1}{4} \int v^2(t,x) dx \leq C( |\dot{x}_0 (t)| \|v(t, \cdot)\|_{L^2}+ |v(t,0)| |v_x(t,0)| )+ \\ & & + C( \|v_x(t, \cdot)\|_{L^2}^2+\|v_x(t, \cdot)\|_{L^2} \|v(t, \cdot)\|_{L^2}+ \|v_x(t, \cdot)\|_{L^2}^{\frac{1}{2}} \|v(t, \cdot)\|_{L^2}^{\frac{3}{2}}). \end{eqnarray*}

We are almost ready to derive an estimate on the desired norms, except that it is hard to control $|v_x(t,0)| $. Indeed, the easiest way out would be to use the norm $\|v_x(t, \cdot)\|_{L^\infty}$, for which we unfortunately lack a priori estimates. Instead, we employ a simple averaging argument to reduce to the quantity $\|v_x(t, \cdot)\|_{L^2}$, for which we do have a priori bound, see (3.2).

Indeed, instead of finding the bound for $ \int v^2(t,x) |x| dx$, we might employ an estimate for $ \int v^2(t,x) |x-a| dx$, where $a\in (-1,1)$. Running through identical estimates, we arrive at

\begin{eqnarray*} & & \frac{1}{2} \partial_t \int v^2(t,x) |x-a| dx + \frac{1}{4} \int v^2(t,x) dx \leq C( |\dot{x}_0 (t)| \|v(t, \cdot)\|_{L^2}+|v(t,a)| |v_x(t,a)| )+\\ & & + C(\|v_x(t, \cdot)\|_{L^2}^2+\|v_x(t, \cdot)\|_{L^2} \|v(t, \cdot)\|_{L^2}+ \|v_x(t, \cdot)\|_{L^2}^{\frac{1}{2}} \|v(t, \cdot)\|_{L^2}^{\frac{3}{2}}). \end{eqnarray*}

Integrating on the time interval $(0,t)$ brings about

\begin{eqnarray*} & & \frac{1}{2} \int v^2(t,x) |x-a| dx + \frac{1}{4} \int_0^t \int v^2(s,x) dx ds \leq \int v_0^2(x) |x-a| dx + \int_0^t |v(s,a)| |v_x(s,a)| ds +\\ & & + C( \int_0^t |\dot{x}_0 (s)| \|v(s, \cdot)\|_{L^2} + \|v_x(s, \cdot)\|_{L^2}^2 +\|v_x(s, \cdot)\|_{L^2} \|v(s, \cdot)\|_{L^2}+ \|v_x(s, \cdot)\|_{L^2}^{\frac{1}{2}} \|v(s, \cdot)\|_{L^2}^{\frac{3}{2}} ) ds. \end{eqnarray*}

Taking an average of the previous estimate, i.e. integrating the previous estimate $\int_{-1}^1 \ldots da$, and in doing so, recall the formula

\begin{equation*} \int_{-1}^1 |x-a| da=\left\{ \begin{array}{cc} 2|x| & |x| \gt 1 \\ x^2+1 & |x|\leq 1 \end{array} \right. \in (2|x|, 2|x|+1). \end{equation*}

This yields the estimate

\begin{align*} & \frac{1}{2} \int v^2(t,x) |x| dx + \frac{1}{2} \int_0^t \int v^2(s,x) dx ds \leq \int v_0^2(x) (2|x|+1) dx + \\ & \quad+ \int_{-1}^1 \int_0^t |v(s,a)| |v_x(s,a)| ) ds da+ \\ & \quad+ C \int_0^t [|\dot{x}_0 (s)| \|v(s, \cdot)\|_{L^2} + \|v_x(s, \cdot)\|_{L^2}^2+ \|v_x(s, \cdot)\|_{L^2} \|v(s, \cdot)\|_{L^2}+ \|v_x(s, \cdot)\|_{L^2}^{\frac{1}{2}} \|v(s, \cdot)\|_{L^2}^{\frac{3}{2}}] ds . \end{align*}

Now, based on (3.2) and Hölder’s inequality

\begin{eqnarray*} & & \int_0^t |\dot{x}_0 (s)| \|v(s, \cdot)\|_{L^2} ds \leq \left( \int_0^t |\dot{x}_0 (s)|^2 ds \right)^{\frac{1}{2}} \left(\int_0^t \|v(s, \cdot)\|_{L^2}^2 ds\right)^{\frac{1}{2}} \leq C \|v\|_{L^2_{t,x}}, \\ & & \int_0^t \|v_x(s, \cdot)\|_{L^2}^2 ds \leq C, \\ & & \int_0^t \|v_x(s, \cdot)\|_{L^2} \|v(s, \cdot)\|_{L^2} ds\leq C \left(\int_0^t \|v_x(s, \cdot)\|_{L^2}^2\right)^{\frac{1}{2}} \left(\int_0^t \|v(s, \cdot)\|_{L^2}^2\right)^{\frac{1}{2}} \leq C \|v\|_{L^2_{tx}}, \\ & & \int_0^t \|v_x(s, \cdot)\|_{L^2}^{\frac{1}{2}} \|v(s, \cdot)\|_{L^2}^{\frac{3}{2}} ds \leq \left(\int_0^t \|v_x(s, \cdot)\|_{L^2}^2 ds\right)^{\frac{1}{4}} \left(\int_0^t \|v(s, \cdot)\|_{L^2}^2 ds\right)^{\frac{3}{4}}\leq C \|v\|_{L^2_{t,x}}^{\frac{3}{2}}, \\ & & \int_{-1}^1 \int_0^t |v(s,a)| |v_x(s,a)| da \leq \left(\int_0^t \|v_x(s, \cdot)\|_{L^2}^2 ds\right)^{\frac{1}{2}} \left(\int_0^t \|v(s, \cdot)\|_{L^2}^2 ds\right)^{\frac{1}{2}} \leq C \|v\|_{L^2_{t,x}}. \end{eqnarray*}

All in all, we conclude that

(3.4)\begin{equation} \int v^2(t,x) |x| dx+ \frac{1}{2} \|v\|_{L^2_{t x}}^2 \leq C \|v\|_{L^2_{t,x}}^{\frac{3}{2}}+C. \end{equation}

We can hide $\|v\|_{L^2_{t,x}}^{\frac{3}{2}}$ on the right-hand side behind $\frac{1}{2} \|v\|_{L^2_{t,x}}^2$. More precisely, by Young’s inequality, for any choice of ϵ > 0,

\begin{equation*} \|v\|_{L^2_{t,x}}^{\frac{3}{2}}\leq \frac{3}{4} \epsilon \|v\|_{L^2_{t,x}}^2+ \frac{1}{4} \epsilon^{-3}. \end{equation*}

With an appropriately small choice of ϵ, we enter this on the right-hand side of (3.4), which leads us to the bound

(3.5)\begin{equation} \sup_{0 \lt t \lt T} \int v^2(t,x) |x| dx+ \int_0^T \|v(t, \cdot)\|_{L^2_{x}}^2 dt \leq C, \end{equation}

for each T > 0. This bound has two important consequences. Firstly, using the Hölder inequality, we have that for each $1 \lt p \lt 2$,

\begin{eqnarray*} & & \int |v(x)|^p dx = \int_{|x| \lt 1} |v(x)|^p dx + \int_{|x|\geq 1} |v(x)|^p dx\leq \\ & & \leq C \|v\|_{L^2}^p+ \left(\int_{|x| \gt 1} v^2(x) |x| dx \right)^{\frac{p}{2}} \left(\int_{|x|\geq 1} |x|^{-\frac{p}{2-p}}dx \right)^{1-\frac{p}{2}}\leq C_p \left(\int |v|^2 (1+|x|) dx \right)^{\frac{p}{2}}, \end{eqnarray*}

for which we have good bounds from (3.5). It follows that for any p > 1, $\sup_{0 \lt t \lt T} \int |v(x)|^p\leq C_p$. Note, however, that the constant Cp blows up as $p\to 1+$, due to $\lim_{p\to 1+} \left(\int_{|x|\geq 1} |x|^{-\frac{p}{2-p}}dx \right)=+\infty$.

Secondly, by taking into account (3.1) and (3.5), we have that for each t > 0,

\begin{equation*} t \|v(t, \cdot)\|_{L^2}^2\leq \int_{0}^{t} \|v(s, \cdot)\|_{L^2}^2 ds \leq \int_0^{\infty} \|v(s, \cdot)\|_{L^2}^2 ds\leq C, \end{equation*}

whence it follows that v obeys the bound $\|v(t, \cdot)\|_{L^2}\leq C t^{-\frac{1}{2}}$. Interpolating this with the uniform bounds for $\|v(t, \cdot)\|_{L^p}, 1 \lt p \lt 2$, obtained just before, we have that for every δ > 0, there is $C=C(\delta)$, so that

\begin{equation*} \|v(t, \cdot)\|_{L^p}\leq C t^{-(1-\frac{1}{p})+\delta}, 1 \lt p \lt 2. \end{equation*}

4. Asymptotic attraction to odd fronts for the generalized Burgers model

In this section, we present an extension of the results in [Reference Barker, Bronski, Hur and Yang2], when the fronts and the respective perturbations are odd.

Now that ${\mathcal L}=-\sum_{j=1}^N a_j (-\partial_{xx})^{\alpha_j} $ preserves the parity of the inputs, it is clear that if u 0 is odd, so is $v_0=u_0-\phi$. Clearly, the evolution preserves the odd solutions, whence $x\to v(t,x)$ is odd for all $t\geq 0$. Recall that by the construction in [Reference Barker, Bronski, Hur and Yang2], the translation function x 0 was selected (see formula (20) there) via the formula $ \dot{x}_0(t)=\gamma \langle \phi_x,v \rangle=0, $ since ϕx is even, while v is odd. It follows by Theorem 1.1 that for this particular front and odd perturbations u 0, $u(t,x)=\phi(x)+v(t,x)$, where

(4.1)\begin{equation} \lim_{t\to \infty} \|v(t, \cdot)\|_{L^p}=0, 2 \lt p\leq \infty. \end{equation}

Recall that v also satisfies the energy estimate (1.13).

We now aim to establish uniform estimates for the norm $\|v(t, \cdot)\|_{L^1(\mathbb R)}$. To this end, take arbitrary large T and assume $0\leq t \lt T$. Note that v satisfies

(4.2)\begin{equation} v_t - v_{xx} - {\mathcal L} v +\phi_x v+\phi v_x+ v v_x=0. \end{equation}

For each p > 1, multiply (4.2) by $v |v|^{p-2}$ and integrate in x. We obtain,

\begin{equation*} \int v_t v |v|^{p-2} dx = \frac{1}{p} \partial_t \int |v|^p dx, \end{equation*}

Integration by parts yields for each p > 1,

\begin{eqnarray*} & & \int v v_xv |v|^{p-2} dx = \int v_x |v|^p dx = - \int v (|v|^p)' dx =- p\int v v_xv |v|^{p-2} dx, \end{eqnarray*}

whence

\begin{equation*} \int v v_xv |v|^{p-2} dx=0. \end{equation*}

By (2.4), we have the coercivity estimates

\begin{equation*} -\int v_{xx} v |v|^{p-2} dx \geq 0, -\int ({\mathcal L} v) v |v|^{p-2} dx =\sum_{j=1}^N a_j \int ((-\partial_{xx})^{\alpha_j} v) v |v|^{p-2} dx \geq 0. \end{equation*}

Next,

\begin{equation*} \int \phi v_x v |v|^{p-2} dx = \frac{1}{p} \int \phi (|v|^p)' dx = - \frac{1}{p} \int \phi_x |v|^p dx. \end{equation*}

All in all, we get

\begin{equation*} \partial_t \int |v|^p dx + (p-1) \int \phi_x |v|^{p} dx \leq 0. \end{equation*}

Denoting $M_T:=\sup_{0 \lt t \lt T} \|v(t, \cdot)\|_{L^\infty}$, and integrating the above expression in $(0,t)$, we obtain

(4.3)\begin{equation} \int |v(t, x)|^p dx\leq \|v_0\|_{L^p}^p+ (p-1) \int_0^t \int_{\mathbb R} |\phi_x(x)| |v(t,x)|^{p} dx\leq \|v_0\|_{L^p}^p+(p-1) \|\phi_x\|_{L^1}T M_T^p. \end{equation}

Note that all the constants in (4.3) are written out explicitly. This allows us to take a limit $p\to 1+$. We obtain

(4.4)\begin{equation} \|v(t, \cdot)\|_{L^1}\leq \|v_0\|_{L^1}, \end{equation}

and this is valid for all T > 0, so (4.4) is a global estimate.

This provides another important ingredient to the estimates (4.1), which we have borrowed from Theorem 1.1. Fix ϵ > 0. By Plancherel’s identity

\begin{equation*} I(t):=\|v(t, \cdot)\|_{L^2}^2=\|v_{ \lt \epsilon}(t, \cdot)\|_{L^2}^2+\|v_{ \gt \epsilon}(t, \cdot)\|_{L^2}^2=:I_{ \lt \epsilon}(t)+I_{ \gt \epsilon}(t). \end{equation*}

We now use (1.13) to conclude

\begin{equation*} \partial_t(I_{ \lt \epsilon}(t)+I_{ \gt \epsilon}(t))\leq - C \|v_x(t)\|_{L^2}^2 \lt - C \|\partial_x v_{ \gt \epsilon}(t)\|_{L^2}^2\leq -C \epsilon^2 I_{ \gt \epsilon}(t). \end{equation*}

It follows that

\begin{equation*} I_{ \gt \epsilon}'(t)+C\epsilon^2 I_{ \gt \epsilon}(t)+ I_{ \lt \epsilon}'(t)\leq 0. \end{equation*}

Multiplying by $e^{C \epsilon^2 t}$ and integrating yields

(4.5)\begin{equation} e^{C \epsilon^2 t} I_{ \gt \epsilon}(t) - I_{ \gt \epsilon}(0) +\int_0^t e^{C \epsilon^2 s} I_{ \lt \epsilon}'(s) ds\leq 0. \end{equation}

However, integration by parts yields

\begin{eqnarray*} -\int_0^t e^{C \epsilon^2 s} I_{ \lt \epsilon}'(s) ds=-I_{ \lt \epsilon} (s)e^{C \epsilon^2 s}|_0^t + C\epsilon^2 \int_0^t e^{C \epsilon^2 s} I_{ \lt \epsilon}(s) ds\leq I_{ \lt \epsilon} (0)+C\epsilon^2 \int_0^t e^{C \epsilon^2 s} I_{ \lt \epsilon}(s) ds. \end{eqnarray*}

Plugging this back in (4.5) and multiplying by $e^{-C \epsilon^2 t}$, yields

(4.6)\begin{equation} I_{ \gt \epsilon}(t)\leq I(0) e^{-C \epsilon^2 t}+ C\epsilon^2 \int_0^t e^{C \epsilon^2 (s-t)} I_{ \lt \epsilon}(s) ds\leq I(0) e^{-C \epsilon^2 t}+\max_{0 \lt s \lt t} I_{ \lt \epsilon}(s). \end{equation}

Now, we apply the Bernstein’s inequality for $A=(-\epsilon, \epsilon)$, and (4.4) to obtain

(4.7)\begin{equation} I_{ \lt \epsilon}(s)= \|v_{ \lt \epsilon}(s, \cdot)\|_{L^2(\mathbb R)}^2\leq C \epsilon \|v(s, \cdot)\|_{L^1(\mathbb R)}^2 \leq C\epsilon \|v_0\|_{L^1(\mathbb R)}^2. \end{equation}

Adding $I_{ \lt \epsilon}(t)$ on both sides of (4.6) and using (4.7), we obtain

(4.8)\begin{equation} \|v(t, \cdot)\|_{L^2}^2=I(t)\leq e^{- C \epsilon^2 t} \|v_0\|_{L^2}^2+ C\epsilon \|v_0\|_{L^1}^2\leq C_1 (e^{- C_1 \epsilon^2 t} +\epsilon)\|v_0\|_{L^2\cap L^1}^2, \end{equation}

for some absolute constant C 1, i.e. independent on $\epsilon, t$.

This inequality is true for every fixed ϵ > 0, and every t, but since clearly the best strategy is to optimize (4.8) in ϵ given t, we should proceed with some care. To this end, fix large $t \gt \gt 1$, and then select the unique $\epsilon \lt \lt 1: e^{- C_1 \epsilon^2 t}=\epsilon$. Clearly, such choice of ϵ is possible, let us compute its asymptotic in terms of t. Since we need to solve,

\begin{equation*} \frac{\epsilon^2}{\ln(\epsilon)}=-\frac{1}{C_1 t}, \end{equation*}

its unique solution obeys the asymptotic

(4.9)\begin{equation} \epsilon\sim \sqrt{\frac{\ln(t)}{t}}. \end{equation}

It follows that

(4.10)\begin{equation} \|v(t, \cdot)\|_{L^2}^2 \leq C_3 \frac{\sqrt{\ln(t)}}{\sqrt{t}} \|v_0\|_{L^2\cap L^1}^2 \end{equation}

We conclude that $\|v(t, \cdot)\|_{L^2}\leq C \|v_0\|_{L^2\cap L^1} \left(\frac{\ln(t)}{t}\right)^{\frac{1}{4}}$, as stated. Interpolating between the bounds (4.10) and (4.4), the Gagliardo-Nirenberg’s inequality now yields

\begin{equation*} \|v(t, \cdot)\|_{L^q}\leq C \|v_0\|_{L^2\cap L^1} \left(\frac{\ln(t)}{t}\right)^{\frac{1-\frac{1}{q}}{2}}, 1 \lt q \lt 2, \end{equation*}

for intermediate values of q. While one can similarly interpolate with the bound (4.1), this will not be a particularly good decay, especially near $q=\infty$.

Instead, we use an alternative estimate based again on (1.13). Indeed from (3.2) and (1.13) , it is clear that for every T > 0

\begin{equation*} \int_T^\infty \|v_x(s, \cdot)\|_{L^2}^2 ds \leq \|v(T, \cdot)\|_{L^2}^2\lesssim \left(\frac{\ln(T)}{T}\right)^{\frac{1}{2}}, \end{equation*}

where the last estimate follows from the L 2 bounds (4.10), Also, for every $T \gt \gt 1$, we have that there exists $T_0\in (T, 2T)$, so that

\begin{equation*} T \|v_x(T_0, \cdot)\|_{L^2}^2 \leq \int_T^{2T} \|v_x(s, \cdot)\|_{L^2}^2 ds \lt \int_T^\infty \|v_x(s, \cdot)\|_{L^2}^2 ds\lesssim \left(\frac{\ln(T)}{T}\right)^{\frac{1}{2}}. \end{equation*}

In other words,

(4.11)\begin{equation} \|v_x(T_0, \cdot)\|_{L^2}^2 \leq C \frac{\sqrt{\ln(T_0)}}{T_0^{\frac{3}{2}}}, \int_{T_0}^\infty \|v_x(s, \cdot)\|_{L^2}^2 ds \leq C \left(\frac{\ln(T_0)}{T_0}\right)^{\frac{1}{2}}. \end{equation}

For large T 0, denote $\delta:=\left(\frac{\ln(T_0)}{T_0}\right)^{\frac{1}{2}} \gt \gt C \frac{\sqrt{\ln(T_0)}}{T_0^{\frac{3}{2}}}$. A this point, we refer to the argument in [Reference Barker, Bronski, Hur and Yang2], which establishes that starting with (4.11), we conclude

(4.12)\begin{equation} \sup_{t \gt T_0} \|v_x(t)\|_{L^2}\leq C \delta^{\frac{2}{3}}\sim \left(\frac{\ln(T_0)}{T_0}\right)^{\frac{1}{3}} \end{equation}

Since this can be done for at least one point $T_0\in (2^k, 2^{k+1}), k \gt \gt 1$, we obtain the bound

\begin{equation*} \|v_x(t)\|_{L^2}\leq C \left(\frac{\ln(t)}{t}\right)^{\frac{1}{3}} \end{equation*}

By Sobolev embedding, followed by the Gagliardo-Nirenberg’s inequality we conclude

\begin{equation*} \|v(t, \cdot)\|_{L^q}\leq \|v(t, \cdot)\|_{\dot{H}^{\frac{1}{2}-\frac{1}{q}}}\leq \|v(t, \cdot)\|_{L^2}^{\frac{1}{2}+\frac{1}{q}}\|v_x(t, \cdot)\|_{L^2}^{\frac{1}{2}-\frac{1}{q}}\leq C \left(\frac{\ln(t)}{t}\right)^{\frac{7}{24}-\frac{1}{12 q}}. \end{equation*}

Footnotes

1 For the definition of the fractional derivatives please consult Section 2 below.

2 Here note that the condition (1.5) is instrumental in establishing (1.7), as ${\mathcal L}$ annihilates the constants.

3 Unfortunately, for all other values of ν, it does not seem likely that such an explicit formula exist, but our arguments herein are independent of such formula.

4 Which appears to be necessary.

5 And it is possible that due to the slow decay of the potential ϕx, the operator $\mathcal H_\epsilon$ may have a complicated point spectrum, including accumulation of e-values at zero.

6 Here the assumption $\alpha_j\leq 1$ is crucial, since it is well-known that the kernels of $e^{-t (-\partial_{xx})^\alpha}$ are sign-changing, as soon as α > 1.

7 But if one uses different values of $\phi_\pm$, the value of γ must change according to $\gamma \gt \frac{2}{\phi_- -\phi_+}$.

8 Recall that $\mathcal H$ has a negative e-value.

References

Achleitner, F, Cuesta, CM and Hittmeir, S (2014) Travelling waves for a non-local Korteweg–de Vries–Burgers equation Journal of Differential Equations 257, 720758.10.1016/j.jde.2014.04.015CrossRefGoogle Scholar
Barker, B, Bronski, J, Hur, V and Yang, Z (2022) Asymptotic stability of sharp fronts. I. One bound state implies stability, https://arxiv.org/abs/2112.04700v2.Google Scholar
Barker, B, Bronski, J, Hur, V and Yang, Z Asymptotic stability of sharp fronts: Analysis and rigorous computation, https://arxiv.org/pdf/2112.04700.Google Scholar
Bona, J, Rajopadhye, S and Schonbek, ME (1994) Models for propagation of bores. I. Two-dimensional theory. Differential and Integral Equations 7, 699734.10.57262/die/1370267701CrossRefGoogle Scholar
Bona, J and Schonbek, M (1985) Traveling wave solutions to Korteweg-de Vries-Burgers equation. Proc Roy Soc Edinburgh 101A, 207226.10.1017/S0308210500020783CrossRefGoogle Scholar
Chamorro, D and Lemarié-Rieusset, PG (2012) Quasi-geostrophic equations, nonlinear Bernstein inequalities and α-stable processes. Revista MatemáTica Iberoamericana 28, 11091122.Google Scholar
Jeffrey, A and Mohamad, M (1991) Exact solutions to the KdV-Burgers equation Wave Motion 14, 369375.10.1016/0165-2125(91)90031-ICrossRefGoogle Scholar
Khodja, M (1989) Nonlinear stability of oscillatory traveling waves for some systems of hyperbolic conservation laws Thesis (Ph.D.)–University of Michigan.Google Scholar
Naumkin, PI and Shishmarev, AI (1991) A problem on the decay of step-like data for the KortewegdeVries-Burgers equation. Funktsionalʼnyi Analiz i ego Prilozheniya 25, 2132.Google Scholar
Naumkin, PI and Shishmarev, AI (1992) A problem on the decay of step-like data for the Korteweg-deVries-Burgers equation. Funktsionalʼnyi Analiz i ego Prilozheniya 26, 8893.Google Scholar
Pego, RL (1985) Remarks on the stability of shock profiles for conservation laws with dissipation. Transactions of the American Mathematical Society 291, 353361.10.1090/S0002-9947-1985-0797065-0CrossRefGoogle Scholar
Su, C and Gardner, C (1969) Korteweg-de Vries equation and generalizations. III. Derivation of the Korteweg-de Vries equation and Burgers equation. Journal of Mathematical Physics 10, 536539.10.1063/1.1664873CrossRefGoogle Scholar