arxiv: v2 [math.ap] 3 Mar 2017

Transcription

1 ON MITTAG-LEFFLER MOMENTS FOR THE BOLTZMANN EQUATION FOR HARD POTENTIALS WITHOUT CUTOFF arxiv: v [math.ap] 3 Mar 17 MAJA TASKOVIĆ, RICARDO J. ALONSO, IRENE M. GAMBA, AND NATAŠA PAVLOVIĆ Abstract. We establish the L 1 weighted propagation properties for solutions of the Boltzmann equation with hard potentials and non-integrable angular components in the collision ernel. Our method identifies null forms by angular averaging and deploys moment estimates of solutions to the Boltzmann equation whose summability is achieved by introducing the new concept of Mittag-Leffler moments - extensions of L 1 exponentially weighted norms. Such L 1 weighted norms of solutions to the Boltzmann equation are, both, generated and propagated in time and the characterization of their corresponding Mittag-Leffler weights depends on the angular singularity and potential rates in the collision ernel. These estimates are a fundamental step in order to obtain L exponentially weighted estimates for solutions of the Boltzmann equation being developed in a follow up wor. 1. Introduction We study generation and propagation in time of L 1 exponentially weighted norms, referred to as exponential moments, associated to probability density functions that solve the Boltzmann equation [1, 11] modeling the evolution of monoatomic rarefied gases. Binary interactions of gas particles are described by transition rates from before and after such interactions, usually referred to as collision ernels. Such ernels are modeled as a product of potential functions of local relative speed and functions of the scattering angle between the pre and post relative velocities. This angular function may or may not be integrable. When integrable, the collision ernel is said to satisfy an angular cutoff condition. The particular case when the angular part of the ernel is bounded, is nown as the Grad s cutoff condition [3]. Otherwise, its non-integrability, referred to as an angular non-cutoff, satisfies specific conditions for details see Section. The concept of exponential moments is associated to the notion of large energy decay rates for tails. A time dependent probability distribution function ft, v is said to have L 1 exponential moment tail behavior of order s and rate rt if, for any fixed t >, R d ft,ve rt v s dv is positive and finite. 1.1 This concept was introduced by Bobylev in [7, 8] and Gamba, Panferov and Villani in [], where they show uniform in time propagation of L 1 Maxwellian tails i.e. Gaussian in v-space, that is s = for several type of collision ernels ranging from 1

2 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ Maxwell-type to hard sphere interactions with angular cutoff conditions, and by Bobylev, Gamba, Panferov in [9] for different values of s,] in the study of inelastic interaction with internal heating sources. These ground breaing wors conceived the idea of controlling exponential moments by proving the summability of power series expansions on a parameter rt. Such formulation was motivated by formally commuting integration in v-space and the infinite sum derived from the power series of the exponential function in 1.1, upon which one obtains R d ft,v r q t v sq Γq +1 = r q t m sq t. 1. Γq +1 The terms m sq t, called polynomial moments, are v sq -weighted L 1 norms of the distribution function fv, t that solves the Boltzmann equation. Representation 1. replaces the quest of L 1 exponential integrability with a given order and rate, with study of summability of infinite sums time series forms. A fundamental technique for accomplishing this tas see[8, 9,, 3, 5, 8] consists of controlling the wea form of the collision operator by the means of angular averaging. These estimates are used to derive a sequence of ordinary differential inequalities for the polynomial moments of the collisional form. These differential inequalities are an algebraic sum of a negative term of moments of highest order and a positive term of bilinear sums of moments of lower orders. Recently Alonso, Cañizo, Gamba and Mohout[1] introduced a new techniquebased on analyzing partial sums corresponding to the infinite sum appearing in 1., to prove the generation of exponential moments with orders up to the potential rate and the propagation of exponential moments with orders up to s =, under an angular integrability condition. It is interesting to note that these results do not rely on the rate of Povzner estimates for angular averaging, and so the resulting order rt may not be optimal. All results mentioned above were developed for the case of an integrable angular collision ernel. This brings us to the setting of this manuscript, the non-cutoff regime. This manuscript focuses on the study of both generation and propagation in time of exponential moments for solutions to the initial value problem for the d-dimensional Boltzmann equation for elastic collisions, in the space homogeneous case, for hard potentials without the angular cutoff assumption. In this direction, Lu and Mouhot [5] showed generation of exponential moments of order up to the potential rate in the collision ernel. In this wor, we considerably extend their result by showing that rates and orders of exponential moments depend on the initial data, as well as potential and angular singularity rates in collision ernels. In order to treat the non-cutoff regime, we develop angular averaged estimates that account for the cancellation of non-integrable angular singularities by means of null forms averaging. The other important component is summability of moments, which is achieved by introducing Mittag-Leffler moments. Indeed, the most significant point of this paper is the introduction of Mittag-Leffler moments, as L 1 Mittag-Leffler weighted norms. They enabled us to extend the range of orders of exponential moments that can be propagated uniformly in time

3 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 3 for the non-cutoff case. To obtain our result, we encounter the need to study 1. whereγq+1isreplacedwithγaq+1, foranonintegera > 1whichisreminiscent ofsomeofthe toolsused in [9], althoughnosumming ofsuch renormalizedmoments was performed there R d ft,v r q t v sq Γaq +1 = r q t m sq t Γaq We observed that the sum appearing on the left-hand side of 1.3 is exactly the well-nown Mittag-Leffler function E a rt v s, where E a is defined as x q E a x := Γaq In analogy to 1.1, this led us to introduce a concept of Mittag-Leffler moments ft,v E a α a v m q tα aq dv = R Γaq +1, 1.5 d which are a natural generalization of exponential moments. Another important aspect of our main result is that the highest order of exponential moment which can be propagated in time, depends continuously on the singularity rate of the angular cross-section. The less singular the angular ernel is, the higher order exponential moment can be propagated. See details in Remar 8. Let us mention one application of L 1 weighted estimates. In [], Gamba, Panferov and Villani, gave a proof to close the open problem of propagation of L - Maxwellian weighted bounds, uniformly in time, to solutions of the Boltzmann equation with hard potential with a cutoff in the angular ernel. Their result follows from an application of a maximum principle of parabolic type, due to the dissipative nature of the collisional integral, and estimates on the Carleman representation of the gain positive part of the collision operator that depend on the L 1 -Maxwellian weighted bounds uniformly propagated in time. We mention here that the extension of such result on propagation of L R d -exponential weights is currently being wored out for the non-cutoff and hard potential case in a forthcoming manuscript [1] using the L 1 weighted estimates obtained in this manuscript. Organization of the paper. Section presents the Boltzmann equation without the angular cutoff condition, exponential and Mittag-Leffer moments and the statements of the two main results of the manuscript - the angular averaged Povzner inequalities with angular singularity cancellation in Lemma.3, and the generation and propagation of Mittag-Leffler moments in Theorem.4. Section 3 contains the proof of the angular averaged Povzner inequalities for non-integrable angular singularity, i.e. Lemma.3. This lemma is the main tool for the formation of ordinary differential inequalities for polynomial moments of all orders, which are covered in Section 4. Section 5 provides details of the proof of the propagation of Mittag-Leffler moments, while in Section 6 we give a new proof of the generation of exponential moments of order up to the rate of potentials. The final section, Appendix, gathers nown and technical yet fundamental results used throughout this manuscript.

4 4 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ. Preliminaries and Main Results.1. The Boltzmann equation. We consider the Cauchy problem for the spatially homogeneous i.e. x-space independent Boltzmann equation { t ft,v = Qf,ft,v, t R +,v R d, d.1 f,v = f v. The function ft,v models the particle density at time t and velocityv of a rarefied gas in which particle collisions are elastic and predominantly binary. The collisional operator Qf, f is a quadratic integral operator defined via Qf,ft,v = f f ff B u,û σdσdv,. R d S d 1 whereweusethe abbreviatednotationf = ft,v, f = ft,v, andf = ft,v. Vectorsv,v denote pre-collisionalvelocitiesand v,v aretheir correspondingpostcollisional velocities. Relative velocity is denoted by u = v v, and its normalization by û = u/ u. Being an elastic interaction of reversible character that conserves momentum v + v = v + v and energy v + v = v + v, pre and postcollisional velocities are related by formulas represented in center of mass V = v +v / and relative velocity u = v v coordinates as follows v = V + u σ, v = V u σ, σ Sd 1..3 The unit vector σ S d 1, referred to as the scattering direction, has the direction of the pre-collisional relative velocity u = v v. We bring to the reader s attention that the pre to post collisional exchange of coordinates satisfy v v = 1 u σ u, v σ û u v u v v v = 1 u σ u. This representation embodies the relation of the exchange of velocity directions as just functions of the relative velocity u and the scattering direction σ. v O The collisional ernel B u,û σ is assumed to tae the form Figure 1. Pre-post collisional velocities B u,û σ = u γ bcosθ,.4 where θ [,π] is the angle between the pre and post collisional relative velocities, and thus it satisfies cosθ = û σ. In this manuscript we wor in the variable hard potentials case, that is < γ 1..5

5 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 5 We assume that the angular ernel is given by a positive measure bû σ over the sphere S d 1. In many models, this function is non-integrable over the sphere, while its weighted integral is finite. In this manuscript we assume that for some β,] the following weighted integral is finite with V d = πd / Γd 1/ being the volume of the d dimensional unit sphere A β := bû σsin β θ dσ S d 1 π = V d bcosθ sin β θ sin d θ dθ <..6 When β = a casethat we donot consider, this condition is nownas the angular cutoff assumption, under which the collisional operator can be split into the gain and loss terms Qf,f = Q + f,f Q f,f,.7 where Q + f,ft,v = f f B u,û σdσdv, R d S d 1 Q f,ft,v = fv f B u,û σdσdv. R d S d 1 In 1963 Grad [3] proposed considering a bounded angular ernel bcosθ and pointed out that different cutoff conditions could be implemented too. Since then the cutoff theory developed extensively, with the belief that removing the singularity of the angular ernel should not affect properties of the equation. Recently, however, it has been observed see for example [4], [14], [15], [16] that the singularity of bcos θ carries regularizing properties. This, in addition to the analytical challenge, motivated further study of the non-cutoff regime. The typical non-cutoff assumption in the literature is the condition.6 with β =. However,weworin the non-cutoffregimewherethe parameterβ,]is allowed to vary and we will see how the strength of the singularity of b influences our main result. In this setting, the splitting.7 is not valid, which is one of the technical challenges that non-cutoff setting brings. In order to address this obstacle we exploit angular cancellation properties for details please see Section 3. Remar 1. In the physically relevant case corresponding to the dimension d = 3, when forces between particles are governed by an inverse power law long range interaction potential φx = Cx p 1, C >, p >, the angular ernel bcosθ has been derived by H. Grad [3] see also [1] and is shown to have the following form bcosθ sinθ C θ 1 ν, θ +, ν = p 1, Note that this model satisfies.6 with any r > ν. γ = p 5, p >..8 p 1

6 6 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ Wea formulation of the collision operator Qf, f. Thans to the symmetries associated to the collisional form Qf,f, defined in the strong form., the collisional operator has a wea formulation that is very important for the analytical manipulation of the equation. Indeed, for any test function φv, v R d, one has see for example [1] 1 Qf,ft,vφvdv = fvfv G φ v,v dv dv,.9 R d R d G φ v,v = φv +φv φv φv B u,û σ dσ. S d 1 The ey aspect of the equation in the wea formulation is expressed in the weight G φ as it carriesall the information about collisions through the collisional ernel B, which is averaged over the unit sphere against test functions φ = φv +φv φv φv. Crucial estimates on the function G φ referred to in the Boltzmann equation literature as Povzner estimates are described below. In the angular cutoff case, positive and negative contributions are treated separately and such estimates are used to estimate the positive part of G φ. A sharp form of angular averaged Povzner estimates from [8, 9, ] is obtained for general test functions φv which are positive and convex. They are crucial for the study of moments summability, the main point of this manuscript. When φv = 1+ v / = v, these estimates, originally developed by Povzner [3], yield ordinary differential inequalities for moment estimates that lead to an existence theory and generation and propagation of moments as developed in Elmroth [18], Desvillettes [13] Wennberg [33] and Mischler, Wennberg [6]. These estimates were also obtained in the non-cutoff case by Wennberg [3] for hard potentials. Uniqueness theory to solutions of the Boltzmann equation for hard potentials was first developed by Di Blassio in [17]. When the angular part of the collision ernel is not integrable, i.e. the non-cutoff case, one needs to expand φ in terms of v v and v v, since both are a multiples of u sin θ/. For this strategy to succeed, the spherical integration variable σ S d 1 must be decomposed as σ = ûcosθ +ωsinθ, corresponding to the polar direction of the relative velocity u, and the azimuthal direction ω S d 1 satisfying u ω =. This decomposition also plays a fundamental role in our derivation of the angular averaged Povzner with singularity cancellation in the proof of Lemma.3. Remar. We note that the identity.9 can also be expressed in a double mixing weighted convolutional form [,, 4] 1 Qf,ft,vφvdv = fvfv u G φ v,ududv R d R d G φ v,u = φv +φv u φv φv ub u,û σ dσ S d 1 since both v and v can be written as functions of v,u and σ from.3, and so the weight function G φ v,u is an average over σ S d 1.

7 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 7.. Moments of solutions to the Boltzmann equation. From the probabilistic viewpoint, moments of a probability distribution density ft, v with respect to the variable v are integrals of such density weighted by functions φv. These are important objects to study as they express average quantities that have significant meaning for the model under consideration. They are the so called observables. In this sense polynomial moments correspond to such integrals for polynomial weights, and exponential moments are for exponential weights. We now recall definitions of polynomial and exponential moments and we here introduce the Mittag-Leffler moments, which are a natural generalization of the exponential moments. Definition.1 Polynomial and exponential moments. Polynomial moment of order q and exponential moment of order s and rate α are respectively defined by: m q t := ft,v v q dv,.1 R d M α,s t := ft,v e α v s dv..11 R d Remar 3. Using the Taylor series expansion, the exponential moment of order s and rate α can also be written as the following sum M α,s t = m qs t α q..1 Remar 4. Polynomial moments can be expressed in terms of the norm of a natural Banach space in the context of the Boltzmann equation. Namely, if we denote L 1 = {f L 1 R d : f v dv = f 1+ v / dv < }, R d R d then Also, note that m q t = f L 1 q t..13 f L 1 q f L 1 q, for any q q..14 Note that this expression is associated to the notion of L 1 exponential tail behavior described in 1.1 and 1.. Consequently, finiteness of exponential moments can be understood as implying that the function ft,v has an exponential tail in v. In this paper, we study whether this property can be generated or propagated in time for the case of variable hard potentials in the non-cutoff case. Because our summability estimates lead to expressions similar to that of.1, yet having Γaq + 1 as a generalization of factorials with non-integer a > 1, we are motivated to use Mittag-Leffler functions, as they are conceived as a generalization

8 8 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ of the Taylor expansion of the exponential function. More precisely, for a parameter a >, Mittag-Leffler function is defined via E a x := x q Γaq Note that for a = 1, the Mittag-Leffler function coincides with the Taylor expansion of the classical exponential function e x. It is also well nown see e.g. [19], page 8. that for any a >, the Mittag-Leffler function asymptotically behaves lie an exponential function of order 1/a, that is E a x e x1/a, asx. Since v is the building bloc for our calculations, we prefer to have x as the argument of Mittag-Leffler function when generalizing e αxs E /s α /s x e αxs, forx..16 Hence, they satisfy the following, with some positive constants c, C c e αxs E /s α /s x C e αxs..17 This motivates our definition of Mittag-Leffler moments. Definition. Mittag-Leffler moment. Mittag-Leffler moment of order s and rate α > of a function f is introduced via R d ft,v E /s α /s v dv..18 Remar 5. In the rest of the paper we will use the fact that Mittag-Leffler moments can be represented as the following sum a time series form, which follows from.15 R d ft,v E /s α /s v dv = Remar 6. Formally, by taing = q s, the above sum becomes m st α Γ +1, s Z m q tα q/s Γ..19 s q +1 that we show it relates to the time series in 1. with the difference being that the summation here goes over the fractions.

9 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 9.3. The main results. There are two important results in this manuscript. The first one relates to the angular averaged Povzner estimate with cancellation. It gives an estimate of the weight function G φ in the wea formulation.9 when the test function is a monomial φv = v rq. We denote this weight function by G rq := G v rq := v rq + v rq v rq v rq B u,û σ dσ. S d 1 Lemma.3. Suppose that the angular ernel bcos θ satisfies the non-cutoff condition.6 with β =. [ Let r,q >. Then the weight function satisfies G rq v,v v v γ A v rq + v rq + A v rq v + v v rq +ε qr/ A qr qr 1 v v v + v qr ],.1 where A = S d π bcosθ sind θdθ is finite by.6. The sequence ε qr/ =: ε q, defined as ε q := A S N π 1 t q 1 sin θ t dt bcosθ sin N θ dθ,. has the following decay properties. If bcos θ satisfies the non-cutoff assumption.6 with β,], then < ε q q 1 β, as q..3 The sequence ε q is the same as in [5]. Its decay properties.3 are also proved in [5], after invoing angular averaging and the dominated convergence theorem. Condition.3 is crucial for finding the highest order s of Mittag-Leffler moment that can be propagated in time. Remar 7. This lemma relies on the polynomial inequality presented in Lemma 3.1. The decay rate of ε q is fundamental for the success of summability arguments, yet is not relevant for the generation and propagation of polynomial moments. In the angular cutoff case when term-by-term techniques were used, the corresponding constant had a rate ε q q r, with r depending on the integrability of b, see [8, 9, ]. When the partial sum technique was employed in [1], the precise rate was not needed any longer. Here however, in the non-cutoff case, the nowledge of the precise decay rate of ε q becomes important again because of extra power of q in the last term of the right-hand side of.3. The second main result, presented as an a priori estimate, consists of two parts. First, under the non-cutoff assumption.6 with β =, we provide a new proof of the generation of exponential moments of order s,γ]. Second, we show the propagation in time of the Mittag-Leffler moments of order s γ,. When s γ,1], β = in the non-cutoff.6 is assumed. When s 1,, the angular ernel is assumed to be less singular. Before we state the theorem, we remind the reader of the following notation L 1 = {f L1 R d : f v dv < }. R d This is the natural Banach norm to solve the Boltzmann equation.

10 1 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ Theorem.4 Generation and Propagation of Exponential-lie moments. Suppose f is a solution to the Boltzmann equation.1 with the collision ernel of the form.4 for hard potentials.5, and with initial data f L 1. a Generation of exponential moments If the angular ernel satisfies the noncutoff condition.6 with β =, then the exponential moment of order γ is generated with a rate rt = α min{t,1}. More precisely, there are positive constants C,α, depending only on b, γ and initial mass and energy, such that ft,v e α min{t,1} v γ dv C, for t..4 R d b Propagation of Mittag-Leffler moments Let s, and suppose that the Mittag-Leffler moment of order s of the initial data f is finite with a rate r = α, that is, R d f v E /s α /s v dv < M..5 Suppose also that the angular cross-section satisfies assumption.6 with β =, if s,1] with β = 4, if s 1,..6 s Then, there exist positive constants C,α, depending only on M, α, b, γ and initial mass and energy such that the Mittag-Leffler moment of order s and rate rt = α remains uniformly bounded in time, that is R d ft,v E /s α /s v dv < C, for t..7 Remar 8. The angular singularity condition β = 4 s in the case of Mittag-Leffler moments of order s 1,, continuously changes from β = for s = 1 to β = for s =. Hence condition β = 4 s continuously interpolates between the most singular ernel typically considered in the literature, which is.6 with β =, and an angular cutoff condition, which corresponds to.6 with β =. This also tells us that in the most singular case one can propagate exponential moments of order s 1, while in angular cutoff cases one can propagate exponential moments of order s to be completely rigorous, Theorem.4 goes up to β > i.e. s <, but [1] already established the case β = i.e. s =. In other words, the less singular the angular ernel is, the higher order exponential moment propagate in time. Remar 9. The propagation result of the theorem can be interpreted in two ways. First, for a Mittag-Leffler or exponential moment of order s to be propagated, the singularity of b should be such that it satisfies.6 with β = 4 s. On the other hand, given an angular ernel b that satisfies condition.6 with a parameter β, ], one can propagate Mittag-Leffler and exponential moments of order s 4 β+.

11 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 11 Remar 1. We note two types of solutions that can be used in the previous theorem. One example are wea solutions, whose existence was proven by Areryd [6] and later extended by Villani [31], under the assumption that initial data has finite mass, energy, entropy and a moment of order +δ, for any δ >. Another type of solutions that could be used are measure wea solutions constructed by Lu and Mouhot [5] see also the result of Morimoto, Wang and Yang [7]. These solutions exist if initial mass and energy are finite, provided that the angular ernel satisfies the following condition π bcosθsind θ1+ logsinθ <, which automatically holds for ernels that satisfy condition.6 with β <. Remar 11. Thans to the fact.17 that Mittag-Leffler functions asymptotically behave lie exponential functions, finiteness of exponential moment of order s is equivalent to finiteness of the corresponding Mittag-Leffler moment. This, in turn implies, as a corollary of Theorem.4 b, the propagation of classical exponential moments. Remar 1. In the case of inverse-power law model described via.8, in which hard potentials correspond to p > 5, the non-cutoff condition.6 is satisfied for β > ν. Hence, Mittag-Leffler moments of orders s < p can be propagated in time. In the graph below y-axis represents the order of exponential tails. The dashed red line mars the highest order of exponential moments that can be generated, while the blue line mars the highest order of Mittag-Leffler moments that can be propagated in time. This graph visually confirms that our propagation result indeed goes beyond the rate of potentials γ. s = p γ = p 5 p 1 p.4. A strategy for proving Theorem.4. Details are provided in Section 5 and Section 6. The proof is inspired by the recent wor [1], where propagation and generation of tail behavior 1. is obtained for angular cutoff regimes. Our goal is to prove that solutions ft,v of the Boltzmann equation for hard potentialsand angularnon-cutoffconditions admit L 1 -Mittag-Lefflermoments with parameters a = s and αt = rt to be found. Because of the asymptotic behavior.16, that would imply that asymptotic limit for large values of v is, indeed, and exponential tail in v-space, with order s and rate rt = αt. Thus, our proof is based on studying partial sums of Mittag-Leffler functions E a α a x, with parameter a = s and with rate αt.

12 1 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ To this end, we wor with n-th partial sums associated to Mittag-Leffler functions, defined as n Ea n α,t = m q t α aq Γaq We need to prove that there exits a positive rate αt and a positive parameter a, both uniform in n, such the sequence of finite sums converges as n. In particular, we need show that Ea n α,t is bounded by a constant independent of time and independent of n. The values for a,α and the bound of the partial sums are found and shown to depend on data parameters given by the collisional ernel characterization and properties of the initial data. In order to achieve all of this, we derive a differential inequality for E n a = En a α,t. The first step in this direction is to obtain differential inequalities for moments m q t, by studying the balance m qt = Qf,ft,v v q dv..9 R d that is a consequence of the Boltzmann equation. The right hand side is estimated by bounding the polynomial moments of the collision operator by non-linear forms of moments m t of order up = q +γ, with < γ 1. This requires finding the estimates of the wea formulation.9 with test functions φv = v. Consequently, we need to estimate the angular integration within the weight function G v qv,v S d 1 v q + v q v q v q bcosθdσ..3 These estimates will lead, thans to.9 and.9, to the following differential inequality for polynomial moments m q K 1 m q+γ + K m q q q +K 3 ε q qq 1 m+γ m q +m m q +γ. 1 =1.31 where K 1 = A C γ, where A was defined in.1 and C γ just depends on the rate of potentials γ. Similarly K and K 3 depend on these data parameters as well. The ey property of this inequality is that the highest order moment of the righthand side comes with a negative sign which is crucial for moments propagation and generation. Another important aspect of this differential inequality is the presence of the factor qq 1 in the last term, which was absent in angular cutoff cases. Because of it, it will be of great importance to now the decay rate for ε q. The second step Section 4 consists in the derivation of a differential inequality for partial sums E n a = E n aα,t obtained by adding n inequalities corresponding to.31 for renormalized polynomial moments m q tα aq /Γaq +1. This will yield d dt En a c q + K 1 Ia,γ n + K 1 c q + K Ea n + ε q q a K 3 CEa n Ia,γ n..3

13 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 13 In particular we obtain an ordinary differential inequality for the partial sum E n a that depends on a shifted partial sum I n a,γ, defined by I n a,γα,t = n m q+γ t α aq..33 Γaq +1 The derivation of the last term in the right hand side of.3 requires a decay property of combinatoric sums of Beta functions. These estimates are very delicate and are presented in detail in Lemma A.4 and Lemma A.5 in the Appendix. The constants K 1,K and K 3 only depend on the singularity conditions.6, and so they are independent of n and on any moment q. The constant c q depends only on a finite number q of moments of the initial data. The choice of q is crucial to control the long time behavior of solutions to inequality.3, and it is done such that ε q q a K 3 < K 1 /, after using condition.3 in Lemma.3. Finally, after showingthat Ia,γ n α,t is bounded below by sum oftwo terms depending linearly on Ea nα,t and on mass m, and nonlinearly on the rate α, we obtain the following differential inequality for partial sums in the case of propagation of initial Mittag-Leffler moments d dt En at K 1 Eat n + K 1m e α1 a α γ α γ +K Propagation estimate. The constant K depends on data parameters characterizing q,c q and K i, i = 1,,3.. In addition, for the generation case, we obtain d dt En γ 1 K1 Eγ n m C q +K Generation estimate. t α Thus, the differential inequalities.3 are reduced to linear ones. Both inequalities have corresponding solutions for choices on parameters a and α that are independent on n and time t, and will depend on q, which depends only on data parameters. 3. Angular averaging lemma This section is about the proof of the angular averaging with cancellation, i.e. Lemma.3, a crucial step for controlling moments and summability of their renormalization by the Gamma function. One of the tools used in the proof is the following estimate on symmetrized convex binomial expansions. Lemma 3.1. [Symmetrized convex binomial expansions estimate] Let a, b, t [,1] and p,1] [,. Then p p ta+1 tb + 1 ta+tb a p b p t1 t a p +b p + t1 t ab p 1 +a p 1 b. 3.1

14 14 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ Proof: Suppose p. The case p,1] can be done analogously. Due to the symmetry of the inequality 3.1, we may without the loss of generality assume that a b. Since all the terms have homogeneity p, the inequality 3.1 is equivalent to showing where Fz is defined by Fz := Fz, z 1, z 1 t1 t p +1 + t1 t z+z p 1 p p. tz+1 t 1 tz+t It is easy to chec that [ F z = p 1 p 1 t1 t z p + t1 tp z p 3 pt p tz +1 t ] p p1 t 1 tz +t. As tz +1 t and 1 tz +t are two convex combinations of z and 1, and since z 1, we have that tz +1 t z and 1 tz +t z. Since p, this implies tz +1 t p z p and 1 tz +t p z p. Therefore, F z p 1 p 1 t1 t z p + t1 tp z p 3 pt z p p1 t z p = t1 tp z p 3. Thus, F z for z 1. So, F z is increasing. Since F 1 =, we have that F z for z 1. Finally using the fact that F1 =, we conclude Fz for z 1. We are now ready to prove the new form of the angular averagingwith cancellation type lemma. For another version, see [5]. Proof of Lemma.3 Recall the definition of the weight G rq G rq v,v := v v γ S d 1 bcosθ sin d θ v rq dσ, 3. where v rq = v rq + v rq v rq v rq. This integral is rigorous even in cases when S d 1 B u,cosθ dσ is unbounded, by an angular cancellation. A natural way of handling the cancellation is to decompose σ S d 1 into θ [,π] and its corresponding azimuthal variable ω S d, i.e. σ = cosθ û+sinθ ω, where S d û = {ω S d 1 : ω û = }. See Figure.

15 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 15 This decomposition allows handling the lac of integrability concentrated at the origin of the polar di- û ˆV rection θ =. However, it requires a specific way of σ decomposing v and v that separates the part that depends on ω. More precisely, v and v θ are decomposed into a sum of a convex combination of the local energies proportional to a function S d j ω of the polar angle θ, and another term depending S d 1 on both the polar angle and ω see the Appendix for details Figure v = E v,v θ + Pθ,ω,.Decomposition3.3 of σ. v = E v,v π θ Pθ,ω. Here Pθ,ω = v v sinθj ω is a null form in ω by averaging, i.e. S d Pθ,ωdω =, and E v,v θ is a convex combination of v and v given by E v,v θ = t v + 1 t v, where t = sin θ. These two fundamental properties mae the weight function G rq v,v well defined for every v and v for sufficiently smooth test functions φ C R d even under the non-cutoff assumption.6 with β =. In fact, Taylor expansions associated to v rq are a sum of a power of E v,v θ, plus a null form in the azimuthal direction, plus a residue proportional to sin θ that will secure the integrability of the angular cross section with respect to the scattering angle θ. Indeed, Taylor expand v rq around Eθ up to the second order to obtain rq v rq = E v,v θ + hsinθj ω 3.4 = E v,v θ rq/ rq + Ev,v θ rq 1 h sinθj ω + rq rq 1 1 h sin θj ω 1 t [ Eθ+th sinθj ω ] rq dt. Similar identity can be obtained for v rq. Since the collisional cross section is independent of the azimuthal integration we will mae use of the following property. Any vector j laying in the plane orthogonal to the direction of u is nullified by multiplication and averaging with respect to the azimuthal direction with respect to u, that is S d j ω dω =. Therefore, we can write G rq v,v as the sum of two integrals on the S d 1 sphere, whose first integrand contains the zero-order order term of the Taylor expansion of,

16 16 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ both, v rq and v rq subtracted by their corresponding un-primed forms, while the second integrand is just the second order term of the Taylor expansion 3.5 G rq v,v = I 1 + I 3.5 π = E v,v θ rq/ +E v,v π θ rq/ v rq v rq bcosθ sin d θ dω dθ S d + rq rq π 1 1 h sin d θ bcosθ j ω 1 t S d [Ev,v θ + th sinθj ω ] rq + [ E v,v π θ th sinθ ] rq dtdωdθ. At this point we use inequality 3.1 to estimate the first integral I 1. We use it with a = v, b = v and t = cos θ, which yields I 1 π S d sin θ v rq + v rq bcosθ sin d θ dθ π sin θ + v rq v + v v rq bcosθ sin d θ dθ = A v rq + v rq + A v rq v + v v. rq The constant A was defined after For the second order term I, we use that j ω 1 and h = v v v v, and that see [5] E v,v θ+thsinθ j ω v + v 1 t 4 sin θ, 3.7 to conclude I rq rq 1 v v S d 1 π sin d θ bcosθ 1 t v + v rq 1 1 t sin θ 4 rq dt dθ. After a simple change of variablest 1 t and recalling the definition of constant ε rq/ in., we see that I ε rq/ A rq rq 1 v v v + v rq. 3.8 Putting together the estimate for I 1 and for I, we obtain the desired estimate on the weight G rq v,v. 4. Ordinary differential inequalities for moments In this section we present two differential inequalities for polynomial moments Proposition 4.1 which will be essential for the proof of Theorem.4. We also

17 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 17 state and prove a result about generation of polynomial moments in the non-cutoff case Proposition 4.. Before we state the proposition, we recall the floor function of a real number, which in the case of a positive real number x R + coincides with the integer part of x x := integer part of x. 4.1 Proposition 4.1. Suppose all the assumptions of Theorem.4 are satisfied. Let q N, and define p = p+1 for any p R to be the integer part of p + 1/. Then for some constants K 1,K,K 3 > depending on γ, bcosθ, dimension d we have the following two ordinary differential inequalities for polynomial moments of the solution f to the Boltzmann equation a The m γ version needed for the generation of exponential moments m γq t K 1m γq+γ + K m γq 4. +K 3 ε qγ/ qγ qγ 1+q γ 1 =1 q γ m γ+γ m γq γ +m γ m γq γ+γ. 1 b The m version needed for propagation of Mittag-Leffler moments m q K 1 m q+γ + K m q 4.3 q q + K 3 ε q qq 1 m+γ m q +m m q +γ. 1 =1 In both cases, the constant K 1 = A C γ, where A was defined in.1 and C γ, to be defined in the proof below, only depends on the γ rate of the hard potentials. Similarly K and K 3, also depend on data, through the dependence on A and C γ. Proof: We start the proof by analyzing m rq with a general polynomial weight v rq. Then by setting r = γ we shall derive a and by setting r = we shall obtain b. Recall that after multiplying the Boltzmann equation.1 by v rq, the wea formulation.9 yields m rq t = 1 ff G rqv,v dv dv. 4.4 R d The weight function G rq can be estimated as in Proposition.3, which yields m rq t A f f v v γ v rq + v rq dvdv R d R d + A f f v v γ v rq v + v v rq dvdv R d R d + A ε rq rq rq/ 1 f f v v γ v v v + v rq dvdv R d R d 4.5

18 18 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ We estimate v v γ via elementary inequalities v v γ Cγ 1 v γ + v γ and v v γ C γ v γ v γ, 4.6 where C γ = min{1, 1 γ } see for example [1]. As an immediate consequence and v v γ v rq + v rq C γ v γ v γ v rq + C γ v γ v γ v rq = C γ v rq+γ + v rq+γ v rq v γ + v γ v rq, 4.7 v v γ v rq v + v v rq Cγ v 1 γ + v γ v rq v + v v rq Cγ 1 v rq v γ + v γ v rq, 4.8 where the last inequality uses Lemma A.1. Combining 4.5 with 4.7 and 4.8 we obtain m rq t A C γ ff v rq+γ + v rq+γ dvdv R d R d + A 1+C 1 γ f f v rq v γ + v γ v rq dvdv R d R d + A ε rq/ C γ rq rq 1 R d ff v γ + v γ v v v + v rq dvdv A C γm tm rq+γ t+ A 1+C 1 γ m γ tm rq t + A ε rq/ rq rq Rd C γ 1 ff v γ + v γ v v v + v rq dvdv. Therefore, since < γ 1, by conservation of mass and energy, m t = m and m γ t m, m rq t K 1m rq+γ t+k m rq t+ K 3 ε rq rq rq/ 1 ff v γ + v γ v v v + v rq dvdv, R d 4.9 where K 1 = A C γ m, K = A 1 + C 1 γ m, and K 3 = A C γ, so these threeconstantsonlydepend onthe initial massandenergy, onthe rateofthepotential γ and on the angular singularity condition.6 that determines the constant A. From here, we proceed to prove a and b separately.

19 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 19 a Setting r = γ in 4.9, applying the following elementary polynomial inequality which is valid for γ,1] v + v γq v γ + v γ q γ, 4.1 and using the polynomial Lemma A. yields m γq t K 1 m γq+γ +K m γq + K 3 ε γq γq γq/ 1 f f v γ + v γ v v v γ + v γ q γ dvdv R d K 1 m γq+γ +K m γq + K 3 ε γq/ q γ = γq γq 1 R d f f v γ + v γ q γ v γ+ v γq γ + v γq γ v γ+ dvdv γq K 1 m γq+γ +K m γq + K 3 ε γq/ q γ = γq 1 q γ m γ++γ m γq γ +m γq γ +γ m γ+ dvdv. Finally, re-indexing to 1 and applying Lemma A.1 yields m γq t K γq γq 1 m γq+γ +K m γq + K 3 ε γq/ 1 1+q γ =1 which completes proof of a. q γ m γ+γ m γq γ +m γq γ+γ m γ dvdv. 1 b Now, we set r = in 4.9 and apply Lemma A. to obtain m q t K 1 m q+γ +K m q + K 3 ε q qq 1 f f v γ + v γ R d q v v = = K 1 m q+γ +K m q + K 3 ε q qq 1 q = q v v q + v q v dvdv R d f f v γ + v γ q v + v q + v q v + dvdv q q = K 1 m q+γ +K m q + K 3 ε q qq 1 m+γ m q +m m q +γ. 1 =1

20 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ The last equality is obtained by re-indexing to 1 and using that 1+ q = q. This completes proof of b. Proposition 4. Polynomial moment bounds for the non-cutoff case. Suppose all the assumptions of Theorem.4 are satisfied. Let f be solution to the homogeneous Boltzmann equation.1 associated to the initial data f. 1 Let the initial mass and energy be finite, i.e. m bounded, then for every p > there exists a constant B rp, depending on rp, γ, m and A from condition.6, such that m rp t B rp max{1,t rp/γ }, for all r R + and t Furthermore, if m rp is finite, then the control can be improved to m rp t B rp, for all r R + and t. 4.1 Proof: These statements can be shown by studying comparison theorems for initial value problems associated with ordinary differential inequalities of the type y t+ay 1+c t Byt, and comparing them to classical Bernoulli s differential equations for the same given initial y. In our context, these inequalities are a result of estimating moments for variable hard potentials, i.e. γ > as indicated in.5. Comparison with Bernoulli type differential equations was classically used in angular cutoff cases in [3, 33, 6, 1]. Also it was used in the proof of propagation of L 1 exponential tails for the derivatives of the solution of the Boltzmann equation by means of geometric series methods in [9,, 3]. In fact, the extension to the non-cutoff case follows in a straightforward way from the moments estimates in Proposition 4.1. Indeed, the moment estimates, from either 4. or 4.3, show that the only negative contribution is on the highest order moment, being either m rq+γ with γ > for r = γ or, respectively. Then, due to the fact that γ >, an application of classical Jensen s inequality with the convex function ϕx = x 1+γ/rp yields m rp+γ t m γ/rp m 1+γ/rp rp t for all t >. Applying this estimate to the negative term in either 4. or 4.3, results in the following estimate m rp B rp m rp K 1 m rp+γ B rp m rp K 1 m 1+γ/rp rp, 4.13 with r either γ in 4., or in 4.3. The constants are K 1 = K 1 γ,a with < γ 1, and A from the angular integrability condition.6; and B rp = B rp K, rp K 3, after using that ε p 1, where K and K 3 also depend on the initial data and collision ernel through γ and A. Therefore, as in [3], we set yt := m rp t, A := K 1, B := B rp and c = γ/rp.

21 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 1 The bound 4.1 then follows by finding an upper solution that solves the associated Bernoulli ODE y t = Byt Ay 1+c t with finite initial polynomial moment y = m rp. This yields that for any t > m rp t [ mrp γ/rp e tbγ/rp + A ] rp/γ B 1 e tbγ/rp [ ] rp/γ A B 1 e tbγ/rp rp/γ A B rp Bγ ebγ/rp rp/γ t rp/γ, t < 1, 1 e Bγ/rp rp/γ, t 1. B rp max{1,t rp/γ }, 4.14 where B rp := K1 B rp rp/γ max { rp γb rp e γbrp/rp rp/γ, 1 e γb rp/rp rp/γ }. Now, since m rp t is a continuous function of time, if m rp is finite for any rp 1, then the bound for strictly positive times we just obtained in 4.14 implies m rp t B rp for possibly different constants B rp. We finally stress that constants B rp depend on rp,γ,m and A from condition Proof of Mittag-Leffler moments propagation Proof of Theorem.4 b. Let us recall representation.19 of the Mittag-Leffler moment of order s and rate α in terms of infinite sums R d ft,v E /s α /s v dv = m q tα q/s Γ. 5.1 s q +1 We introduce abbreviated notation a = s, and note that since s,, we have 1 < a := s <. 5. We consider the n-th partial sum, denoted by Ea n, and the corresponding sum, denoted by Ia,γ n, in which polynomial moments are shifted by γ. In other words, we consider E n a α,t = n m q t α aq Γaq +1, In a,γ α,t = n m q+γ t α aq. Γaq +1

22 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ For each n N, define T n := sup{t E n aα,τ < 4M, for all τ [,t}. 5.3 where the constant M is the one from the initial condition.5. This parameter T n is well-defined and positive. Indeed, since α will be chosen to be, at least, smaller than α, then at time t = we have E n a = n m q α aq Γaq +1 < m q α aq Γaq +1 = f ve /s α /s v dv < 4M, uniformly in n. Therefore, since partial sums are continuous functions of time they are finite sums and each m q t is also continuous function in time t, we conclude that E n a α,t < 4M holds for t on some positive time interval denoted [,t n with t n > and hence T n >. Next, we loo for an ordinary differential inequality that the partial sum E n aα,t satisfies, following the steps presented in Subsection.4. We start by splitting d dt En a α,t into the following two sums, where index q will be fixed later, and then apply the moment differential inequality 4.3 q 1 d dt En a α,t = q 1 m qtα aq Γaq +1 m q tαaq Γaq +1 K 1 + n q=q m n m q+γtαaq Γaq +1 q=q qt α aq Γaq +1 n m + K qtαaq Γaq +1 q=q n ε q qq + K 1αaq q q 3 m +γ m q +m m q +γ Γaq +1 1 q=q =1 =: S K 1 S 1 + K S + K 3 S We estimate each of the four sums S,S 1,S and S 3 separately, with the goal of comparing each of them to the functions Ea nα,t and In a,γ α,t. We remar that the most involving term is S 3. It resembles the corresponding sum in the anglar cutoff case [1], with a crucial difference that our sum S 3 has two extra powers of q, namely qq 1. Therefore, a very sharp calculations is required to control the growth of S 3 as a function of the number q of moments. This is achieved by an appropriaterenormalizationof polynomial moments within S 3 and also by invoing the decay rate of associated combinatoric sums of Beta functions developed in the Appendix A. The term S can be bounded by terms that depends on the initial data and the parameters of the collision cross section. Indeed, from Lemma 4., the propagated polynomial moments can be estimated as follows m p B p and m p B p B p, for anyp >, 5.5 where the constant B p defined in 4.14 depends on γ, the initial p-polynomial moment m p and A from condition.6.

23 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 3 In particular, for < γ < 1, we can fix q, to be chosen later, such that the constant c q := max p I q {B p, B p B p }, with I q = {,...,q +1} 5.6 depends only on q, γ, A from condition.6, and the initial polynomial moments m q, for q I q. Thus, due to the monotoncity of L 1 norms with respect to as presented in.14, both the q-moments and its derivatives, as well as the shifted moments of order q +γ, are controlled by c q as follows m q t, m q+γ t, m q t c q, for all q {,1,,...q }, 5.7 Therefore, for q fixed, to be chosen later, S is estimated by S := q 1 m q αaq Γaq +1 c q q 1 α aq Γaq +1 q 1 α a q c q Γq +1 c q e αa c q, 5.8 for the parameter α small enough to satisfy α < ln 1/a, or equivalently, e αa. 5.9 The second term S 1 is crucial, as it brings the negative contribution that will yield uniform in n and global in time control to an ordinary differential inequality for Ea nα,t. In fact, S 1 is controlled from below by Iaγ n α,t as follows. S 1 := n m q+γ αaq Γaq +1 q=q = I n a,γ q 1 m q+γ α aq Γaq +1, so using 5.7 and the estimate just obtained for S in 5.8, yields the bound from below q S 1 Ia,γ n c 1 α aq q Γaq +1 In a,γ c q. 5.1 The sum S is a part of the partial sum E n a, so S E n a While this term is positiveit will need to be lowerorderthan the onein the negative part of the right hand side. Finally, we estimate S 3 and show that it can be bounded by the product of E n aα,t and I n a,γ α,t. We worout the detailsofthe first term in the sum S 3 := S 3,1 +S 3,, that is the one with m +γ m q. The other sum with m m q +γ can be bounded by following a similar strategy. In order to generate both the partial sum

24 4 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ Ea nα,t and the shifted one In a,γ α,t, we mae use of the following well nown relations between Gamma and Beta functions see also Appendix A. Ba +1,aq +1 = = Γa +1Γaq +1 Γa +1+aq +1 Γa +1Γaq +1 Γaq Therefore, multiplying and dividing products of moments m +γ m q in S 3,1, by Γa +1Γaq +1 yields n ε q qq S 3,1 := 1αaq q q m +γ m q Γaq +1 1 q=q =1 =1 n q q m+γ α a = ε q qq 1 1 Γa +1 q=q Ba +1,aq +1 m q α aq Γaq +1 Γaq + Γaq +1. Notethatthefactors m +γα a Γa+1 and m q α aq Γaq +1 arethebuildingblocsofia,γα,t n and Ea n α,t, respectively. Next,sinceΓaq+/Γaq+1 = aq+1,usingtheinequality a b a b, it follows that n q S 3,1 ε q aq +1qq 1 q=q =1 q =1 q 1 m +γ α a Γa +1 Ba +1,aq +1 m q α aq. Γaq Next we show that the factor q aq +1qq 1 q Ba +1,aq +1 1 =1 on the right hand side of 5.13 grows at most as q a. Indeed, using Lemma A.4, the sum of the Beta functions is bounded by C a aq 1+a. Therefore, S 3,1 is estimated by n q S 3,1 C a ε q q a m +γ α a m q α aq, 5.14 Γa +1 Γaq +1 q=q =1

25 MITTAG-LEFFLER MOMENTS TO SOLUTIONS OF BTE 5 where C a is a possibly different constant that depends on a. Now, by Lemma.3, the factor ε q q a decreases monotonically to zero as q if the angular ernel bcosθ satisfies.6 with β = a. Hence, ε q q a ε q q a, for any q q, 5.15 and thus the term S 3,1 is further estimated by S 3,1 C a ε q q a n q q=q =1 m +γ α a Γa +1 m q α aq Γaq +1. Finally, inspired by [1], we bound this double sum by the product of partial sums E n a I n a,γ. To achieve that, change the order of summation to obtain S 3,1 C a ε q q a C a ε q q a n n = max{q, 1} n = m +γ α a Γa +1 m +γ α a Γa +1 n max{q, 1} m q α aq Γaq +1 m q α aq Γaq C a ε q q a I n a,γ E n a, obtaining the expected control of S 3,1. As mentioned above the estimate of the companion sum S 3, follows in a similar way, so we can assert S 3 C a ε q q a Ea n tin a,γ t Next we obtain an ordinary differential inequality for Eat n depending only on data parameters and Ia,γ n t. Indeed, combining 5.8, 5.1, 5.11 and 5.16 with 5.4 yields d dt En a K 1I n a,γ + c q 1+K 1 + K E n a + ε q q a C a K 3 I n a,γ En a Since, by the definition of time T n, the partial sum E n a is bounded by the constant 4M on the time interval [,T n ], we can estimate, uniformly in n, the following two terms in 5.18 c q 1+K 1 +K E n a c q 1+K 1 +4K M =: K, 5.19 where K depends only on the initial data and q still to be determined. Thus, factoring out Ia,γ n from the remaining two terms in 5.18 yields d dt En a Ia,γ n K 1 ε q q a C a K 3 Ea n + K Ia,γ n K 1 4ε q q a C a K 3 M + K, 5.

26 6 M. TASKOVIĆ, R. J. ALONSO, I. M. GAMBA, AND N. PAVLOVIĆ where in the last inequality we again used that, by the definition of T n, we have Ea n 4M on the closed interval [,T n ]. Now, since ε q q a converges to zero as q tends to infinity by Lemma.3 as bcosθ satisfies.6 with β = a, we can choose large enough q so that For such choice of q we then have K 1 4ε q q a C a K 3 M > K d dt En a K 1 In a,γ + K. 5. The final step consists in finding a lower bound for I n a,γ in terms of En a. The following calculation follows from a revised form of the lower bound given in [1], I n a,γ t := n 1 α γ/ = 1 α γ/ 1 α γ/ 1 α γ/ m q+γ α aq Γaq +1 n n n v q α aq v 1 Γaq +1 α v q+γ α aq v 1 Γaq +1 α ft,v dv Rd v q α aq n Γaq +1 ft,v dv n Ea nt Rd α q α aq ft,v dv Γaq +1 E n at m α qa 1 Γaq +1 ft,v dv v q α aq v < 1 Γaq +1 α ft,v dv > 1 α γ E n at 1 α γ m e αa Therefore, applying inequality 5.3 to 5. yields the following linear differential inequality for the partial sum E n a d dt En a t K 1 E α γ a n t + K 1m e α1 a α γ +K. Then, by the maximum principle for ordinary differential inequalities, K E/s n t = En at M + αγ/ 1 m e α1 a + K K 1 α γ. = M + m e α1 a + αγ/ K 1 K 4M,