Propagation of stretched exponential moments for the Kac equation and Boltzmann equation with Maxwell molecules

TL;DR

This paper proves the propagation of stretched exponential moments over time for solutions to the Boltzmann and Kac equations with Maxwell molecules, using Mittag-Leffler moments to handle angular singularities.

Contribution

It establishes the propagation of stretched exponential moments for Maxwell molecules' Boltzmann and Kac equations, including non-cutoff cases, using Mittag-Leffler moments.

Findings

Propagation of moments depends on angular kernel singularity.

Results apply to both cutoff and non-cutoff cases.

Uses Mittag-Leffler moments to handle singularities.

Abstract

We study the spatially homogeneous Boltzmann equation for Maxwell molecules, and its $1$-dimensional model, the Kac equation. We prove propagation in time of stretched exponential moments of their weak solutions, both for the angular cutoff and the angular non-cutoff case. The order of the stretched exponential moments in question depends on the singularity rate of the angular kernel of the Boltzmann and the Kac equation. One of the main tools we use are Mittag-Leffler moments, which generalize the exponential ones.

Equations

$${\partial }_{t}f(t,v)={\int }_{\mathbb{R}}{\int }_{-\pi }^{\pi }\left(f^{\prime }{f}_{\ast }^{\prime }-f{f}_{\ast }\right)b_K(|\theta |)\mathrm{d}\theta \mathrm{d}{v}_{\ast }.$$

$${\partial }_{t}f(t,v)={\int }_{{\mathbb{R}}^d}{\int }_{S^{d-1}}\left(f^{\prime }{f}_{\ast }^{\prime }-f{f}_{\ast }\right)|v-v_{\ast }{|}^{\gamma }{b}_B\left(\frac{v-v_{\ast }}{|v-v_{\ast }|}\cdot \sigma \right)\mathrm{d}\sigma \mathrm{d}{v}_{\ast },$$

$${\partial }_{t}f(t,v)={\int }_{{\mathbb{R}}^d}{\int }_{S^{d-1}}\left(f^{\prime }{f}_{\ast }^{\prime }-f{f}_{\ast }\right)b_B\left(\frac{v-v_{\ast }}{|v-v_{\ast }|}\cdot \sigma \right)\mathrm{d}\sigma \mathrm{d}{v}_{\ast }.$$

$$m_q(t):={\int }_{{\mathbb{R}}^d}f(t,v)⟨v{⟩}^q\mathrm{d}v.$$

$$M_{\alpha ,s}(t):={\int }_{{\mathbb{R}}^d}f(t,v)e^{\alpha ⟨v{⟩}^s}\mathrm{d}v,\alpha >0.$$

$$M_{\alpha ,s}(t)={\int }_{{\mathbb{R}}^d}f(t,v)\sum _{q=0}^{\infty }\frac{{\alpha }^q⟨v{⟩}^{sq}}{q!}=\sum _{q=0}^{\infty }\frac{{\alpha }^q{m}_{sq}(t)}{\Gamma (q+1)}.$$

$${\partial }_{t}f(t,v)=K(f,f)(t,v),$$

$$K(f,f)(t,v)={\int }_{\mathbb{R}}{\int }_{-\pi }^{\pi }\left(f^{\prime }{f}_{\ast }^{\prime }-f{f}_{\ast }\right)b_K(|\theta |)\mathrm{d}\theta \mathrm{d}{v}_{\ast },$$

$$v^{\prime 2}+v_{\ast }^{\prime 2}=v^2+v_{\ast }^2,$$

$$v^{\prime }$$

$$v_{\ast }^{\prime }$$

$${\int }_{-\pi }^{\pi }{b}_K(|\theta |){\sin }^{\kappa }\theta \mathrm{d}\theta <\infty ,\text{for some}\kappa \in [0,2].$$

$$\begin{gathered} {\int }_{\mathbb{R}}K(f,f)\varphi (v)\mathrm{d}v \\ =\frac{1}{2}{\int }_{\mathbb{R}}{\int }_{\mathbb{R}}{\int }_{-\pi }^{\pi }f{f}_{\ast }\left(\varphi (v^{\prime })+\varphi (v_{\ast }^{\prime })-\varphi (v)-\varphi (v_{\ast })\right)b_K(|\theta |)\mathrm{d}\theta \mathrm{d}{v}_{\ast }\mathrm{d}v. \end{gathered}$$

$$\{\begin{array}{rl}{\partial }_{t}f(t,v)& =K(f,f)(t,v)t\in {\mathbb{R}}_{+},v\in \mathbb{R},\\ f(0,v)& =f_0(v),\end{array}$$

$${\int }_{\mathbb{R}}{f}_0(v)\left(⟨v{⟩}^2+\left|\log f_0(v)\right|\right)\mathrm{d}v<\infty .$$

$${\partial }_t{\int }_{\mathbb{R}}f\varphi (v)\mathrm{d}v={\int }_{\mathbb{R}}{\int }_{\mathbb{R}}{K}^{\varphi }(v,v_{\ast })f{f}_{\ast }\mathrm{d}v\mathrm{d}{v}_{\ast },$$

$$K^{\varphi }(v,v_{\ast })={\int }_{-\pi }^{\pi }\left(\varphi (v^{\prime })-\varphi (v)\right)b_K(|\theta |)\mathrm{d}\theta .$$

$${\int }_{\mathbb{R}}f(t,v)\mathrm{d}v={\int }_{\mathbb{R}}{f}_0(v)\mathrm{d}v.$$

$${\int }_{\mathbb{R}}f(t,v)v^2\mathrm{d}v={\int }_{\mathbb{R}}{f}_0(v)v^2\mathrm{d}v,$$

$${\int }_{\mathbb{R}}{f}_0(v)\left(1+|v{|}^{2p}\right)\mathrm{d}v<C,$$

$${\partial }_{t}f(t,v)=Q(f,f)(t,v).$$

$$Q(f,f)(t,v)={\int }_{{\mathbb{R}}^d}{\int }_{S^d-1}\left(f^{\prime }{f}_{\ast }^{\prime }-f{f}_{\ast }\right)|v-v_{\ast }{|}^{\gamma }{b}_B\left(\hat{u}\cdot \sigma \right)\mathrm{d}\sigma \mathrm{d}{v}_{\ast },$$

$$v^{\prime }+v_{\ast }^{\prime }$$

$$|v^{\prime }{|}^2+|v_{\ast }^{\prime }{|}^2$$

$$v^{\prime }$$

$$v_{\ast }^{\prime }$$

$$\gamma =0.$$

$${\int }_0^{\pi }{b}_B(\cos \theta ){\sin }^{\beta }\theta {\sin }^{d-2}\theta d\theta <\infty ,\text{for some}\beta \in [0,2].$$

$$b_B(\cos \theta )\sin \theta \sim {\theta }^{-\frac{3}{2}},\theta \to 0.$$

$$\begin{gathered} {\int }_{{\mathbb{R}}^d}Q(f,f)\varphi (v)\mathrm{d}v \\ =\frac{1}{2}{\int }_{{\mathbb{R}}^{2d}}f{f}_{\ast }|v-v_{\ast }{|}^{\gamma }{\int }_{S^{d-1}}\left(\varphi (v^{\prime })+\varphi (v_{\ast }^{\prime })-\varphi (v)-\varphi (v_{\ast })\right)b_B(\hat{u}\cdot \sigma )\mathrm{d}\sigma \mathrm{d}{v}_{\ast }\mathrm{d}v. \end{gathered}$$

$$\{\begin{array}{rl}{\partial }_{t}f(t,v)& =Q(f,f)(t,v)t\in {\mathbb{R}}_{+},v\in {\mathbb{R}}^d,\\ f(0,v)& =f_0(v),\end{array}$$

$${\int }_{{\mathbb{R}}^d}{f}_0(v)\left(1+|v{|}^2+\log (1+f_0(v))\right)dv<+\infty .$$

Full text

Milana Pavić-Čolić

Department of Mathematics and Informatics

Faculty of Sciences, University of Novi Sad

Trg Dositeja Obradovića 4, 21000 Novi Sad, Serbia

Maja Tasković

Department of Mathematics

University of Pennsylvania

David Rittenhouse Lab.

209 South 33rd Street, Philadelphia, PA 19104

1. Introduction

In this paper we study exponential tails (exponentially weighted $L^{1}$ norms) of weak solutions to the Kac equation and the spatially homogeneous Boltzmann equation for Maxwell molecules. We show propagation in time of such tails, both in the so-called cutoff and the non-cutoff case.

Both the Kac equation and the Boltzmann equation model the evolution of a probability distribution of particles inside a gas interacting via binary collisions. The models we consider are spatially homogeneous, which means that the probability distribution $f(t,v)$ depends only on time $t$, velocity $v$, but not of the spatial variable $x$. The Kac equation is a model for a 1-dimensional spatially homogeneous gas in which collisions conserve the mass and the energy, but not the momentum. On the other hand, the spatially homogeneous Boltzmann equation describes a gas in a $d$-dimensional space, with $d\geq 2$, in which particles collisions are elastic, meaning they conserve the mass, momentum and energy.

The probability distribution function $f(t,v)$, for time $t\in\mathbb{R}^{+}$ and velocity $v\in\mathbb{R}^{d}$ (with $d=1$ for the Kac equation and $d\geq 2$ for the Boltzmann equation), changes due to the free transport and collisions. In the case of the spatially homogeneous Kac equation its evolution is modeled by the following equation

[TABLE]

The spatially homogeneous Boltzmann equation on the other hand reads

[TABLE]

which in the case of Maxwell molecules ($\gamma=0$) reduces to

[TABLE]

Details about the notation employed are contained in Section 2 for the Kac equation and in Section 3 for the Boltzmann equation. For now we only remark that for both equations we consider angular kernels $b_{K}$ and $b_{B}$ that may of may not be integrable. When the angular singularity is non-integrable, our results depend on the singularity rate of the kernels.

The Kac equation (1) and the corresponding Boltzmann equation for Maxwell molecules (3) share many properties (one notable difference is that the Kac equation does not conserve the momentum). In particular, both equations propagate polynomial and exponential moments, whose definitions we now recall.

Definition 1.1**.**

The polynomial moment of order $q$ of the distribution function $f$ is defined by

[TABLE]

Definition 1.2**.**

The stretched exponential moment of order $s$ and rate $\alpha$ of the distribution function $f$ is defined by

[TABLE]

In this paper, the special case when $s=2$ is referred to as the Maxwellian moment.

When $f$ solves the Kac equation, the dimension $d$ in these formulas is one. We also remark that we use the following notation $\langle x\rangle:=\sqrt{1+x_{1}^{2}+\dots x_{d}^{2}}$, for any $x=(x_{1},\dots,x_{d})\in\mathbb{R}^{d}$, $d\geq 1$. The results presented in this paper are also valid when the moments are defined with absolute values $|v|$ in place of $\langle v\rangle$.

In the case of the Kac equation, the study of stretched exponential moments goes back to [8]. There, the constant angular kernel is considered, and the propagation of stretched exponential moments of orders $s=1$ and $s=2$ is proved.

For the Boltzmann equation, propagation of Maxwellian moments was proved in the case of Maxwell molecules $\gamma=0$ in [3, 4] and recently in [6] via Fourier transform techniques. The theory was later extended to hard potentials $\gamma\in(0,1]$ in the context of Maxwellian moments in [5, 7, 10, 6], and in the context of stretched exponential moments in [12, 1]. Finally, the propagation of stretched exponential moments for the non-integrable angular kernels are studied in [11, 13].

In this paper, we generalize results of [8] to include more general orders of stretched exponentials, namely $s\in(0,2]$. The angular kernels that we study are more general and may or may not be integrable. In the case of non-integrable angular kernels, the singularity rate affects the order of moments that propagate in time. In addition, we apply same technique to prove propagation of stretched exponential moments for the Boltzmann equation with $\gamma=0$, thus extending the result of [13].

We point out that the method we employ in this paper differs from the approach in [8]. Elegant calculations for exponential moments (5) of order $s=1$ and $s=2$ in [8] are done directly at the level of exponential moments. In this manuscript, we take a different route and express exponential moments as infinite sums of polynomial moments and then strive to show that such infinite sums are finite. Such approach has been first developed in the context of the Boltzmann equation in [5], where the following fundamental relation was noted

[TABLE]

Finiteness of such sums can be studied by proving term-by-term geometric decay, or by showing that partial sums are uniformly bounded.

Our proof is inspired by the works [1, 13], where the partial sum approach is developed. Moreover, motivated by [13], we exploit the notion of Mittag-Leffler moments, which serve as a generalization of stretched exponential moments and which are very flexible for the calculations at hand. We recall the definition and the motivation for Mittag-Leffler moments in Section 5.

The paper is organized as follows. A brief review of the Kac equation in provided in Section 2, while the review of the Boltzmann equation is contained in Section 3. In Section 4 we state our main result. Section 5 recalls the notion of Mittag-Leffler functions and moments, one of the main tools in the proof of the main theorem. Section 6 contains another key tool - an angular averaging lemma with cancellation. In Section 7, the angular averaging lemma is used to derive differential inequalities satisfied by polynomial moments of the solution to the Cauchy problem under the consideration. Finally, in Section 8 we provide the proof of the main theorem. The Appendix lists auxiliary Lemmas.

2. The spatially homogeneous Kac equation

The Kac model statistically describes the state of the gas in one dimension. The main object is the distribution function $f(t,x,v)\geq 0$ which depends on time $t\geq 0$, space position $x\in\mathbb{R}$ and velocity $v\in\mathbb{R}$, and which changes in time due to the free transport and collisions between gas particles. Assuming that collisions are binary and that they conserve mass and energy, but not momentum, the evolution of the distribution function is determined by the Kac equation.

In this paper we assume that the distribution function does not depend on the space position $x$, i.e. $f:=f(t,v)$. In that case, $f$ satisfies the spatially homogeneous Kac equation

[TABLE]

where the collision operator $K(f,f)$ is defined by

[TABLE]

with the standard abbreviations $f_{*}:=f(t,v_{*})$, $f^{\prime}:=f(t,v^{\prime})$, $f^{\prime}_{*}:=f(t,v^{\prime}_{*})$.

The velocities $v^{\prime},v^{\prime}_{*}$ and $v,v_{*}$ denote the pre and post-collisional velocities for the pair of colliding particles, respectively. A collision conserves the energy of the two particles

[TABLE]

so by introducing a parameter $\theta\in[-\pi,\pi]$, the collision rules read

[TABLE]

Note that a 2-dimensional vector $(v^{\prime},v^{\prime}_{*})$ can be viewed as a rotation of the 2-dimensional vector $(v,v_{*})$ by the angle $\theta$.

In this paper we assume that the angular kernel $b_{K}(|\theta|)\geq 0$ satisfies the following assumption

[TABLE]

The case $\kappa=0$ corresponds to the so called Grad’s cutoff case, when the angular kernel is integrable on $[-\pi,\pi]$. Otherwise, when $\kappa$ is strictly positive, i.e. the non-cutoff case, $b_{K}(|\theta|)$ is allowed to have $\kappa$ more degrees of singularity at $\theta=0$.

2.1. Weak formulation of the collision operator

Since the Jacobian of the transformation (9) is unit, for a test function $\phi(v)$, the weak formulation of the collision operator $K(f,f)$ reads

[TABLE]

2.2. Weak solutions of the Kac equation

We recall the definition of a weak solution to the Cauchy problem for the Kac equation

[TABLE]

whose existence was proved in [8] for the cutoff case, i.e. $\kappa=0$, and in [9] for the non-cutoff case, i.e. $\kappa\in(0,2]$.

Definition 2.1**.**

Let $f_{0}\geq 0$ be a function defined on $\mathbb{R}^{d}$ with finite mass, energy and entropy, i.e.

[TABLE]

Then we say $f\geq 0$ is a weak solution to the Cauchy problem (12) with $K(f,f)$ given by (8) if $f(t,v)\in L^{\infty}\left(\left[0,+\infty\right);L^{1}_{2}\right)$, and for all test functions $\phi\in W^{2,\infty}(\mathbb{R}_{v})$ we have

[TABLE]

where

[TABLE]

For these solutions conservation of mass holds

[TABLE]

while the energy decreases in time. However, the energy is conserved, that is

[TABLE]

if there exists $C>0$ such that

[TABLE]

for some $p\geq 2$. For details, see [9].

3. The spatially homogeneous Boltzmann equation

The state of gas particles which at a time $t\in\mathbb{R}^{+}$ have a position $x\in\mathbb{R}^{d}$ and velocity $v\in\mathbb{R}^{d}$, $d\geq 2$, is statistically described by the distribution function $f=f(t,x,v)\geq 0$. The evolution of such distribution function is modeled by the Boltzmann equation, which takes into account the effects of the free transport and collisions on $f$. The collisions are assumed to be binary and elastic, that is, they conserve mass, momentum and energy for any pair of colliding particles.

When the distribution function is independent of the spatial variable $x$, that is $f:=f(t,v)\geq 0$, which is the so called spatially homogeneous case, the Boltzmann equation reads

[TABLE]

The collision operator $Q(f,f)$ is defined by

[TABLE]

with the standard abbreviations $f_{*}:=f(t,v_{*})$, $f^{\prime}:=f(t,v^{\prime})$, $f^{\prime}_{*}:=f(t,v^{\prime}_{*})$. For a pair of particles, vectors $v^{\prime},v^{\prime}_{*}$ denote pre-collisional velocities, while vectors $v,v_{*}$ denote their post-collisional velocities. Local momentum and energy are conserved, i.e.

[TABLE]

Thus by introducing a parameter $\sigma\in S^{d-1}$, the collision laws can be expressed as

[TABLE]

The unit vector $\sigma\in S^{d-1}$ has the direction of the relative velocity $u^{\prime}=v^{\prime}-v^{\prime}_{*}$, while the normalization of the relative velocity $u=v-v_{*}$ is denoted by $\hat{u}:=\frac{u}{|u|}$. The angle between these two directions, denoted by $\theta$, is called the scattering angle and it satisfies $\hat{u}\cdot\sigma=\cos\theta$.

Due to physical considerations, the parameter $\gamma$ is a number in the range $(-d,1]$. In this paper we consider the Maxwell molecules model, which corresponds to

[TABLE]

The angular kernel $b_{B}(\hat{u}\cdot\sigma)=b_{B}(\cos\theta)$ is a non-negative function that encodes the likelihood of collisions between particles. It has a singularity for $\sigma$ that satisfies $\hat{u}\cdot\sigma=1$, i.e. $\theta=0$, which may or may not be integrable in $\sigma\in S^{d-1}$. Its integrability is often referred to as the angular cutoff, while its non-integrability is referred to as the non-cutoff case. In this paper we assume that

[TABLE]

The case $\beta=0$ corresponds to $b_{B}(\hat{u}\cdot\sigma)$ being integrable in $\sigma\in S^{d-1}$, i.e. it corresponds to the cutoff case. When $\beta>0$, then the angular kernel $b_{B}$ is allowed to have $\beta$ more degrees of singularity compared to the cutoff case.

In particular, in the case of inverse power-law potentials for the Maxwell molecules, the interaction potential in 3 dimensions is of the form $V(r)=r^{-4}$. Then a nonintegrable singularity of the function $b_{B}$ is known

[TABLE]

Therefore, $\beta$ should satisfy $\beta>\frac{1}{2}$.

3.1. Weak formulation of the collision operator

Since the Jacobian of the pre to post collision transformation is unit and due to the symmetries of the kernel, for any sufficiently smooth test function $\phi(v)$, the weak formulation of the collision operator $Q(f,f)$ reads

[TABLE]

3.2. Weak solutions to the Boltzmann equation

We recall the definition of a weak solution to the Cauchy problem for the Boltzmann equation

[TABLE]

whose existence in three dimensions and for the angular kernel (19) with $\beta\in[0,2]$ is proved in [2, 14].

Definition 3.1**.**

Let $f_{0}\geq 0$ be a function defined in $\mathbb{R}^{d}$ with finite mass, energy and entropy

[TABLE]

Then we say $f$ is a weak solution to the Cauchy problem (21) if it satisfies the following conditions

• •

$f\geq 0,\;f\in C(\mathbb{R}^{+};\mathcal{D}^{\prime}(\mathbb{R}^{d}))\,\cap\,f\in L^{1}([0,T];L^{1}_{2+\gamma})$

• •

$f(0,v)=f_{0}(v)$

• •

$\forall t\geq 0$:   $\int f(t,v)\psi(v)dv=\int f_{0}(v)\psi(v)dv$, for $\psi(v)=1,v_{1},...,v_{d},|v|^{2}$

• •

$f(t,\cdot)\in L\log L$ and $\forall t\geq 0:\,\int f(t,v)\log f(t,v)dv\leq\int f_{0}(v)\log f_{0}dv$

• •

$\forall\phi(t,v)\in C^{1}(\mathbb{R}^{+},C^{\infty}_{0}(\mathbb{R}^{3}))$, $\forall t\geq 0$ we have that

[TABLE]

4. The main results

Our main result establishes propagation of stretched exponential moments for the Kac equation and for the Boltzmann equation corresponding to Maxwell molecules.

Theorem 4.1**.**

Suppose initial datum $f_{0}\geq 0$ has finite mass, energy and entropy, i.e. (13) in case of the Kac equation and (22) in the case of the Boltzmann equation.

• (a)

Kac equation: Let $f(t,v)$ be an associated weak solution to the Cauchy problem (12), with (8) and with the angular kernel satisfying (10) with $\kappa\in[0,2]$. If

[TABLE]

then for every $\alpha_{0}>0$ there exists $0<\alpha\leq\alpha_{0}$ and a constant $C>0$ (depending only on the initial data and $\kappa$) so that

[TABLE]

• (b)

Boltzmann equation for Maxwell molecules: Let $f(t,v)$ be an associated weak solution to the Cauchy problem (21) with the angular kernel satisfying (19) with $\beta\in[0,2]$. If

[TABLE]

then for every $\alpha_{0}>0$ there exists $0<\alpha\leq\alpha_{0}$ and a constant $C>0$ (depending only on the initial data and $\beta$) so that

[TABLE]

Remark 1**.**

We make several remarks about this result.

1. (i)

The order $s$ of the stretched exponential moment that propagates in time depends on the singularity rate of the angular kernel. According to (23) and (25), the more singular the kernel is, the smaller the $s$ is. 2. (ii)

The Maxwellian moment $s=2$ can be reached only in the cutoff case i.e. $\kappa=0$ for the Kac equation or $\beta=0$ for the Boltzmann equation. 3. (iii)

The cutoff Kac equation was studied in [8], where propagation of moments of order $s=1$ and $s=2$ was proved. We extend this result by allowing $s\in[0,2]$, and by considering the non-cutoff kernels too. 4. (iv)

The cutoff Boltzmann equation for Maxwell molecules was studied in [3, 4], where propagation of Maxwellian moments ($s=2$) is proved. We extend this result by allowing $s\in[0,2]$. 5. (v)

The non-cutoff Boltzmann equation for hard potentials $\gamma>0$ was studied in [11] where generation of stretched exponential moments of order $s=\gamma$ was proved, and [13] where propagation of stretched exponential moments was proved depending on the singularity rate of the angular kernel. We extend the work of [13] to include the case $\gamma=0$.

5. Mittag-Leffler moments

In this section, we recall the definition of Mittag-Leffler moments, first introduced in [13]. They are a generalization of stretched exponential moments, and they are convenient for the study of exponential decay properties of a function $f$. Namely, these moments are the $L^{1}$ norms weighted with Mittag-Leffler functions which asymptotically behave like exponentials. More precisely, a Mittag-Leffler function with parameter $a>0$ is defined by

[TABLE]

Note that $\mathcal{E}_{1}(x)$ is simply the Maclaurin series of $e^{x}$, while it is well-known that for $a>0$

[TABLE]

Therefore,

[TABLE]

This motivated the definition of Mittag-Leffler moment [13]

Definition 5.1**.**

The Mittag-Leffler moment of a rate $\alpha>0$ and an order $s>0$ is defined via

[TABLE]

for any $t\geq 0$.

Remark 2**.**

Due to the asymptotic behavior of Mittag-Leffler functions, the finiteness of the stretched exponential moment $M_{\alpha,s}(t)$ at any time $t>0$ is equivalent to the finiteness of the corresponding Mittag-Leffler moment $\mathcal{M}_{\alpha,s}(t)$.

6. Angular averaging lemmas with cancellation

Before proving Theorem 4.1, we provide an estimate of the angular part of the weak formulation (11) and (20) when the test function is a monomial $\phi(v)=\langle v\rangle^{2q}$. These bounds will be later used to derive a differential inequality for polynomial moment in Lemma 7.1.

Lemma 6.1**.**

Let $q\geq 2$.

• (a)

Kac equation: Suppose that the angular kernel of the Kac equation $b_{K}$ satisfies the assumption (10). Then

[TABLE]

where

[TABLE]

and

[TABLE]

• (b)

Boltzmann equation for Maxwell molecules: Suppose that the angular kernel $b_{B}$ satisfies the assumption (19).

[TABLE]

where

[TABLE]

and

[TABLE]

Remark 3**.**

The sequences $\{\varepsilon_{\kappa,q}\}_{q}$ and $\{\varepsilon_{\beta,q}\}_{q}$ are decreasing to zero with a certain decay rate depending on the angular singularity rate $\kappa\in[0,2]$ in the case of the Kac equation and $\beta\in[0,2]$ in the case of the Boltzmann equation, [11],

[TABLE]

Proof of Lemma 6.1.

The proof of part (b) can be found in [13, Lemma 2.3]. Thus, here we provide only the proof of part (a).

If $E(\theta)$ denotes the following convex combination of particle energies

[TABLE]

then, using the collision rules (9), we obtain

[TABLE]

Taylor expansion of $\langle v^{\prime}\rangle^{2q}$ around $E(\theta)$ up to the second order yields

[TABLE]

Analogous expression can be written for $\langle v^{\prime}_{*}\rangle$ as well.

The first order term in the above expression is an odd function in $\theta$, which nullifies by integration over the even domain $[-\pi,\pi]$. Therefore, we can write

[TABLE]

where

[TABLE]

We now proceed to estimate the terms $I_{1}$ and $I_{2}$ separately.

Term $I_{1}$

The term $I_{1}$ is estimated by an application of Lemma A.1. Indeed, for $t=\cos^{2}\theta$, so that $1-t=\sin^{2}\theta$, and $a=\langle v\rangle^{2}$, $b=\langle v_{*}\rangle^{2}$, recalling (33) we obtain:

[TABLE]

Therefore, recalling (27), we have

[TABLE]

Term $I_{2}$

By Cauchy-Schwartz inequality, $-2\,t\,vv_{*}\sin\theta\cos\theta\leq t\,E(\pi-\theta)$. Thus,

[TABLE]

where the last inequality follows from

[TABLE]

The same estimate holds for $E(\pi-\theta)+2tvv_{*}\sin\theta\cos\theta$. Therefore, recalling (28) the definition of $\varepsilon_{\kappa,q}$, we have

[TABLE]

Adding estimates for terms $I_{1}$ and $I_{2}$ completes the proof of lemma. ∎

7. Bounds on polynomial moments

Lemma 7.1**.**

With assumptions and notations of Lemma 6.1, we have the following differential inequalities for a polynomial moment $m_{2q}$

• (a)

In the case of the Kac equation,

[TABLE]

• (b)

In the case of the Boltzmann equation for Maxwell molecules,

[TABLE]

Proof.

Multiplying the Kac equation (7) with $\langle v\rangle^{2q}$ and integrating with respect to $v$, we obtain an equation for the polynomial moment $m_{2q}$

[TABLE]

Using the weak formulation (11) one has

[TABLE]

Applying Lemma 6.1 and Lemma A.2 yields

[TABLE]

It remains to change index $k$ in the sum and the part (a) is proven. The proof of the part (b) can be done in an analogous way. ∎

Lemma 7.2** (Propagation of polynomial moments ).**

Suppose $f$ is a weak solution to the Cauchy problem of either Kac equation (12) or the Boltzmann equation for Maxwell molecules (21). Then for every $q\geq 0$, we have

[TABLE]

where the constant $C^{*}_{q}>0$, is uniform in time, and depends on $q$ and the first $q$ moments of the initial data.

Proof.

Applying the inequality (57) to the differential inequality (35) yields

[TABLE]

where $C_{q}=C_{1}m_{2}(0)+C_{1}q(q-1)\,2^{q-2}$. Therefore,

[TABLE]

Applying this inequality inductively, for an integer $q\in\mathbb{N}$, we have

[TABLE]

Therefore, every even moment is bounded uniformly in time.

Moments whose order is not an even integer can be interpolated by even moments. For example, if $0<2q-2\leq p\leq 2q$, then

[TABLE]

Hence, polynomial moments of non-even orders are bounded uniformly in time as well.

∎

Remark 4**.**

We also note that derivatives of polynomial moments are uniformly bounded in time. Namely, applying the inequality (57) to the differential inequality (35) yields

[TABLE]

where $C_{q}=C_{1}m_{2}(0)+4q(q-1)\,\varepsilon_{\kappa,q}\,2^{q-2}$. Lemma 7.2 implies that $m_{2q-2}(t)$ is bounded uniformly in time by a constant $C^{*}_{q-1}$. Thus,

[TABLE]

8. Proof of Theorem 4.1

Proof of Theorem 4.1 (a).

Recall from Remark 2 that finiteness of the stretched exponential moment $M_{\alpha,s}(t)$ is equivalent to finiteness of the Mittag-Leffler moment of the same rate and order. Therefore, we set out to prove finiteness of Mittag-Leffler moment of order $s$ and rate $\alpha$ that will be determined later:

[TABLE]

where

[TABLE]

The case $a=1$ corresponds to $s=2$ i.e. Maxwellian moments. Propagation of such moments can be esablished according to (23) only in the cutoff case $\kappa=0$. This result (propagation of Maxwellian moments in the cutoff case) was alrady established in [8]. Thus, we here focus on the case when $a>1$.

The goal is to prove that partial sums of (43)

[TABLE]

are bounded uniformly in time and $n$.

From the differential inequality for polynomial moments (35), and by denoting

[TABLE]

for any $q_{0}\geq 3$, we obtain the following differential inequality for the partial sum $E^{n}$:

[TABLE]

We proceed to estimate each $S_{i}$, $i=0,1,2,3$ separately.

For later purposes, we introduce a constant which will be an upper bound for the first $q_{0}-1$ polynomial moments and their derivatives. Let,

[TABLE]

where $C_{q}^{*}$ is the constant from Lemma 7.2, and $C_{q}C_{q-1}^{*}$ is from Remark 4. Then

[TABLE]

Term $S_{0}$

Since the mass is conserved (14), i.e. $m_{0}^{\prime}=0$, the first term in the sum $S_{0}$ is equal to zero. Hence, by (46) we have

[TABLE]

Reindexing the sum and using the monotonicity of the Gamma function $\Gamma(aq+a+1)\geq\Gamma(q+1)$, for $a>1$, $q=0,1,2,\dots$, we obtain

[TABLE]

for $\alpha$ small enough so that

[TABLE]

Term $S_{1}$

Using the bound (46) and parameter $\alpha$ chosen so that (48) holds, we have

[TABLE]

Term $S_{2}$

Using again the monotonicity of the Gamma function, we obtain

[TABLE]

Term $S_{3}$

Using the property of the Beta function $B(x,y)=\frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}$, the term $S_{3}$ can be rearranged

[TABLE]

where the last estimate follows by an application of Lemma A.4.

Since by (23) we have

[TABLE]

Remark 3 implies that $\left\{q^{2-a}\varepsilon_{\kappa,q}\right\}_{q}$ is a decreasing sequence and

[TABLE]

If we denote $c_{a}=aC_{a}$, then the monotonicity of $\left\{q^{2-a}\varepsilon_{\kappa,q}\right\}_{q}$ yields

[TABLE]

Going back to (44) and applying the bounds (8), (8), (8), (8) we obtain a differential inequality for the partial sum $E_{n}$

[TABLE]

Due to the conservation of mass, i.e. $m_{0}(t)=m_{0}(0)$, and the dissipation of energy, i.e. $m_{2}(t)\leq m_{2}(0)$ for the weak solution $f$, we have

[TABLE]

To show that such $E^{n}(t)$ is uniformly bounded in time and $n$, we define

[TABLE]

where $M_{0}$ is the bound on the initial data in (4.1), with the goal of proving that $T_{n}=\infty$ for all $n\in\mathbb{N}$.

The number $T_{n}$ is well-defined and positive. Indeed, since $\alpha<\alpha_{0}$, at time $t=0$ we have

[TABLE]

uniformly in $n$, by (4.1). Since $E^{n}(t)$ are continuous functions of $t$, $E^{n}(t)<4M_{0}$ for $t$ on some positive time interval $[0,t_{n})$, $t_{n}>0$. Therefore, $T_{n}>0$.

Also, since $E^{n}(t)\leq M_{0}$ holds on the time interval $[0,T_{n}]$, from (53) we obtain the following differential inequality

[TABLE]

We conclude that

[TABLE]

First, since $q_{0}^{2-a}\varepsilon_{\kappa,q_{0}}$ converges to zero as $q_{0}$ tends to infinity, we can choose $q_{0}$ large enough so that

[TABLE]

Then, we choose $\alpha$ sufficiently small so that

[TABLE]

Therefore, applying estimates (55), (56) and $m_{0}(0)\leq M_{0}$ to the differential inequality (8) yields

[TABLE]

for any $t\in[0,T_{n}]$. Therefore, the strict inequality $E^{n}<4M_{0}$ holds on the closed interval $[0,T_{n}]$ for each $n$. But, since $E^{n}(t)$ is continuous function in $t$, the inequality $E^{n}(t)<4M_{0}$ holds on a slightly larger interval $[0,T_{n}+\mu)$, $\mu>0$. This contradicts definition of $T_{n}$ unless $T_{n}=+\infty$ for all $n$. Therefore,

[TABLE]

Hence, letting $n\rightarrow\infty$, we conclude that $\mathcal{M}_{\alpha,s}(t)<4M_{0}$ for all $t\geq 0$.

Part (b) of Theorem (4.1) can be proved completely analogously to the proof of part (a). This is due to the similarity of the differential inequalities for polynomial moments (35) and (36). In addition, according to (31) the decay rate of the sequences $\varepsilon_{k},q$ and $\varepsilon_{\beta},q$ depends in the exact same way on the the singularity rate of the angular kernel ($\kappa$ for the Kac equation and $\beta$ for the Boltzmann equation). ∎

Appendix A Auxiliary results

Lemma A.1**.**

Let $a,b\geq 0$, $t\in[0,1]$ and $p\in(0,1]\cup[2,\infty)$. Then

[TABLE]

Lemma A.2**.**

Assume $p>1$ and let $k_{p}=\lfloor\frac{p+1}{2}\rfloor$. Then, for all $x,y>0$ the following inequality holds:

[TABLE]

Lemma A.3**.**

Let $b\leq a\leq\frac{p}{2}$. Then for any $x,y\geq 0$

[TABLE]

Remark 5**.**

The above lemma is useful for comparing products of moments whose total homogeneity is the same. Namely,

[TABLE]

Lemma A.4**.**

Let $p\geq 3$ and $k_{p}=\lfloor\frac{p+1}{2}\rfloor$. Then for any $a>1$ we have

[TABLE]

where the constant $C_{a}$ depends only on $a$.

Acknowledgments

Authors are grateful to Laurent Desvillettes for suggesting us to study the problem considered in this paper. We also thank Ricardo J. Alonso for discussions on the subject. We would like to thank Irene M. Gamba and Nataša Pavlović for their valuable remarks. The work of M. P.-Č. was partially supported by the NSF grant NSF-DMS-RNMS-1107465 and by Project No. ON174016 of the Serbian Ministry of Education, Science and Technological Development. The work of M.T. has been supported by NSF grants DMS-1413064, NSF-DMS-RNMS-1107465 and DMS-1516228. Support from the Institute of Computational Engineering and Sciences (ICES) at the University of Texas Austin is gratefully acknowledged.