Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions and Power Like Tails

Similar documents
SOLUTIONS OF THE LINEAR BOLTZMANN EQUATION AND SOME DIRICHLET SERIES

ON THE SELF-SIMILAR ASYMPTOTICS FOR GENERALIZED NON-LINEAR KINETIC MAXWELL MODELS

Modelling and numerical methods for the diffusion of impurities in a gas

Exponential tail behavior for solutions to the homogeneous Boltzmann equation

Moment Inequalities and High-Energy Tails for Boltzmann Equations with Inelastic Interactions

Exponential moments for the homogeneous Kac equation

Existence of global solutions to the Cauchy problem for the inelastic Boltzmann equation with near-vacuum data

REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS

Classical inequalities for the Boltzmann collision operator.

On Weak Solutions to the Linear Boltzmann Equation with Inelastic Coulomb Collisions

Anomalous transport of particles in Plasma physics

Kinetic models of Maxwell type. A brief history Part I

Spectral - Lagrangian methods for Collisional Models of Non - Equilibrium Statistical States

On high energy tails in inelastic gases arxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Oct 2005

Review of Fundamental Equations Supplementary notes on Section 1.2 and 1.3

Exact Solutions of the Einstein Equations

A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Global Weak Solution to the Boltzmann-Enskog equation

A Solution of the Two-dimensional Boltzmann Transport Equation

Fractional diffusion limit for collisional kinetic equations: A moments method

Velocity Distribution in a Viscous Granular Gas

Entropic structure of the Landau equation. Coulomb interaction

The Equivalence of Ergodicity and Weak Mixing for Infinitely Divisible Processes1

The first order quasi-linear PDEs

Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones

Chapter 2 Non-linear Equations

Anomalous velocity distributions in inelastic Maxwell gases

Tests of a convergent numerical scheme for nonlinear Boltzmann-like models with chemical reactions: two, three and four species

Random Averaging. Eli Ben-Naim Los Alamos National Laboratory. Paul Krapivsky (Boston University) John Machta (University of Massachusetts)

Stochastic Particle Methods for Rarefied Gases

On the Boltzmann equation: global solutions in one spatial dimension

at time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))

On the Optimal Choice of Coefficients in a Truncated Wild Sum and Approximate Solutions for the Kac Equation

Copyright by Ricardo Jose Alonso 2008

The existence of Burnett coefficients in the periodic Lorentz gas

Fluid Equations for Rarefied Gases

Asymptotic stability of an evolutionary nonlinear Boltzmann-type equation

A conservative spectral method for the Boltzmann equation with anisotropic scattering and the grazing collisions limit

Kinetic theory of gases

Isometric elastic deformations

A fast spectral method for the inelastic Boltzmann collision operator and application to heated granular gases

On the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions

QUANTUM MODELS FOR SEMICONDUCTORS AND NONLINEAR DIFFUSION EQUATIONS OF FOURTH ORDER

The Boltzmann Equation and Its Applications

Tanaka Theorem for Inelastic Maxwell Models

1 Review of di erential calculus

z x = f x (x, y, a, b), z y = f y (x, y, a, b). F(x, y, z, z x, z y ) = 0. This is a PDE for the unknown function of two independent variables.

Boltzmann Equation and Hydrodynamics beyond Navier-Stokes

Bulk equations and Knudsen layers for the regularized 13 moment equations

Riemann Zeta Function and Prime Number Distribution

Complex Analysis MATH 6300 Fall 2013 Homework 4

Topological properties of Z p and Q p and Euclidean models

Nonclassical Particle Transport in Heterogeneous Materials

Author(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)

LARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION

Monte Carlo methods for kinetic equations

Physical Modeling of Multiphase flow. Boltzmann method

In this chapter we study several functions that are useful in calculus and other areas of mathematics.

Entropy-dissipation methods I: Fokker-Planck equations

Multiple Random Variables

HW2 Solutions. MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22]

Landau s Fermi Liquid Theory

HW 3 Solution Key. Classical Mechanics I HW # 3 Solution Key

hal , version 1-22 Nov 2009

Lecture 9. = 1+z + 2! + z3. 1 = 0, it follows that the radius of convergence of (1) is.

The propagation of chaos for a rarefied gas of hard spheres

Laplace s Equation. Chapter Mean Value Formulas

From Newton s law to the linear Boltzmann equation without cut-off

Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones

arxiv:math/ v1 [math.ap] 28 Oct 2005

First-Passage Statistics of Extreme Values

From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction

Asymptotic of Enumerative Invariants in CP 2

Differential equations, comprehensive exam topics and sample questions

Hylomorphic solitons and their dynamics

Kinetic Models for Granular Flow

A Generalization of Bernoulli's Inequality

Is Weak Pseudo-Hermiticity Weaker than Pseudo-Hermiticity?

Type II hidden symmetries of the Laplace equation

Thermodynamic limit for a system of interacting fermions in a random medium. Pieces one-dimensional model

arxiv: v1 [math.ap] 11 Jan 2014

A comparative study of no-time-counter and majorant collision frequency numerical schemes in DSMC

On the Validity of the Assumption of Local Equilibrium in Non-Equilibrium Thermodynamics

Math 266, Midterm Exam 1

Derivation of quantum hydrodynamic equations with Fermi-Dirac and Bose-Einstein statistics

Equations of linear stellar oscillations

Generalized Burgers equations and Miura Map in nonabelian ring. nonabelian rings as integrable systems.

Long time behavior of non autonomous Fokker Planck equations and the cooling of granular gases

Relevant sections from AMATH 351 Course Notes (Wainwright): 1.3 Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): 1.1.

Complex Homework Summer 2014

Properties of the S-matrix

Mathematical modelling of collective behavior

Pure Quantum States Are Fundamental, Mixtures (Composite States) Are Mathematical Constructions: An Argument Using Algorithmic Information Theory

Mathematics 104 Fall Term 2006 Solutions to Final Exam. sin(ln t) dt = e x sin(x) dx.

Fluid dynamics for a vapor-gas mixture derived from kinetic theory

Spiral Wave Chimeras

EXPONENTIAL AND ALGEBRAIC RELAXATION IN KINETIC MODELS FOR WEALTH DISTRIBUTION

Monte Carlo method projective estimators for angular and temporal characteristics evaluation of polarized radiation

Non-real zeroes of real entire functions and their derivatives

Transcription:

Journal of Statistical Physics ( C 6 ) DOI: 1.17/s1955-6-944-8 Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions and Power Like Tails A. V. Bobylev 1 and I. M. Gamba Received July 1, 5; accepted: February 15, 6 We consider the Boltzmann equations for mixtures of Maxwell gases. It is shown that in certain limiting case the equations admit self-similar solutions that can be constructed in explicit form. More precisely, the solutions have simple explicit integral representations. The most interesting solutions have finite energy and power like tails. This shows that power like tails can appear not just for granular particles (Maxwell models are far from reality in this case), but also in the system of particles interacting in accordance with laws of classical mechanics. In addition, non-existence of positive self-similar solutions with finite moments of any order is proven for a wide class of Maxwell models. KEY WORDS: Boltzmann equation gas mixture Maxwell models exact solutions high energy tails To Carlo Cercignani, on his 65th Birthday 1. INTRODUCTION In this paper we continue the study of self-similar solutions for various physical systems described by the Boltzmann equations with Maxwell collision kernels. (1 4) It was understood long ago (5) that the key mathematical tool for such equations is the Fourier transform in the velocity space. A detailed analytical theory of the classical spatially homogeneous Boltzmann equation for Maxwell molecules was mainly completed in the 8 s (see (6) for a review). There were almost no new essential results in this field during the 9 s, except for some interesting pure mathematical publications, in particular. (7,8) 1 Department of Mathematics, Karlstad University, Karlstad, SE-65188 Sweden; e-mail: alexander. bobylev@kau.se Department of Mathematics, The University of Texas at Austin, Austin, TX 7871-18 U.S.A.; e-mail: gamba@math.utexas.edu -4715 C 6 Springer Science+Business Media, Inc.

Bobylev and Gamba Quite unexpectedly, Maxwell models again became the subject of many publications in the beginning of the s. The starting point was an idea to use such models for descriptions of inelastic (granular) gases. Inelastic Maxwell models were introduced in (1) (see also (9) for the one dimensional case). It was clear that all the analytical techniques previously developed for classical Maxwell models can be used in the inelastical case almost without changes. Many references to papers published in s by physicists can be found in the book. (1) One interesting result (absent in the elastic case) was the appearance of selfsimilar solutions with power like tails. It was conjectured in (11) and later proved in, (3,4) that such solutions represent asymptotic states for a wide class of initial data. On the other hand, inelastic Maxwell models are just a rough approximation for the inelastic hard sphere model and they give usually wrong answers to the question of large velocity asymptotics (a mathematically rigorous study of such asymptotics for a hard sphere model can be found in (1) ). Some new results in the theory of classical (elastic) Boltzmann equation for Maxwell molecules were also recently published in: (,3) Self-similar solutions (two of which were found in explicit form) and the proof that such solutions represent a large time asymptotics for initial data with infinite energy, clarification of the old Krook-Wu conjecture, (13) etc. It is clear that both elastic and inelastic Maxwell models must be studied from a unified point of view. An interesting question arises in connection with power-like tails for high velocities: Is it possible to observe a similar effect (an appearance of power-like tails from initial data with exponential tails) in the system of particles interacting according to laws of classical mechanics (i.e. without inelasticity assumption)? This is the main question for the present paper. We shall see below that the answer to this question is probably affirmative. The key idea is to consider a mixture of classical Maxwell gases and to find a corresponding limiting case for which such behavior can be in principle observed. The paper is organized as follows. First we consider the Boltzmann equation for Maxwell mixtures and pass to the Fourier representation (Sec ). Then we study a binary mixture and show that corresponding equations formally admit a class of self-similar solutions (Sec 3). In order to simplify the problem we pass to the limit that corresponds to a one component gas in the presence of the thermostat with fixed temperature T (Sec 4). The general problem can be reduced to the case T = (cold thermostat). Then we consider the case of the model cross section (pseudo-maxwell molecules with isotropic scattering) and construct a family of exact self-similar solutions (Sec 5). These solutions are studied in detail in Sections 6 and 7. The new solutions have a lot in common with exact solutions from. () They, however, have finite energy and therefore are more interesting for applications. We did not try to prove neither existence of such solutions in the more general case nor to show that they are large time asymptotic states for a wide

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions class of initial data. This is done, as a particular case, in our paper (14) jointly with C. Cercignani. Instead, we prove in Section 8 a general statement (applied, in particular, to inelastic Maxwell models and elastic models in the cold thermostat) concerning non-existence of positive self-similar solutions with finite moments of any order. Thus it is sufficient in many cases to prove that such positive solution does exist, then it always has just a finite number of even integer moments for all values of parameters of the equation.. MAXWELL MIXTURES We consider a spatially homogeneous mixture of N Maxwell gases. Each component of the mixture is characterized by the molecular mass m i and the distribution function f i = f i (v, t), i = 1,...,N, v R 3 and t R + denote velocity and time respectively. The distribution functions are normalized in such a way that R 3 dv f i (v, t) = ρ i ρ i is the number density of the i th component of the mixture. Note that the quantities ρ i, i = 1,...,N are preserved in time. The interaction between particles is described by the matrix of Maxwell type differential cross-sections σ ij ( u,θ) = 1 u g ij(cos θ), i, j = 1,...,N. u is the relative speed of colliding particles, θ [,π] is the scattering angle. In case of true Maxwell molecules, i.e. particles interacting with potentials U ij = α ij r, α 4 ij >, r > denotes a distance between interacting particles, the following formulas are valid (15) ( ) 1/ αij g ij (cos θ) = g(cos θ), m ij = m im j, i, j = 1,...,N. m ij m i + m j We shall assume below the same kind of formulas for g ij (cos θ) with an arbitrary function g(cos θ) (pseudo-maxwell particles). The corresponding system of Boltzmann equations reads f i t = N j=1 R 3 S dv dω g ij ( u ω u ) [ f i (v ) f j (v ) f i(v) f j (v )] (.1)

Bobylev and Gamba the pair (v,v ) are pre-collisional velocities v = (m i v + m j v + m j u ω)(m i + m j ) 1, v = (m iv + m j v m i u ω)(m i + m j ) 1 ; i, j = 1,...,N (.) with respect to the post-collisional velocities (v, v ). The Fourier transform of Eqs. (.1) (.), for ϕ j (k, t) = dv f j (v, t) e ik v ; k R 3 R 3 leads to equations (6) ϕ i t = N S(ϕ i,ϕ j ) (.3) J=1 ( ) k ω S(ϕ i,ϕ j ) = dω g ij [ϕ i (k + ij S k, t) ϕ j (k ij, t) ϕ i(k, t) ϕ j (, t)] (.4) k + ij = m ik + m j k ω, k ij = m j (k k ω) m i + m j m i + m j so, k = k + ij + k ij ; i, j = 1,...,N. We note that ( ) m (k j ij ) = 4 k s, k + ij m i + m = k m i m j (1 4 j (m i + m j ) s s = 1 ( 1 k ω ), s 1. k Then, we consider isotropic solutions of Eqs. (.3) (.4) and obtain ϕ i (k, t) = ϕ i ( k m i, t ( k ) t ϕ i, t = m i ), ), N S( ϕ i, ϕ j ), (.5) j=1

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions with and S( ϕ i, ϕ j ) = [ ( k ) dsg ij (s) ϕ i (1 β ij s), t m i ϕ j ( k m j β ij m j m i s, t ) ϕ i ( k m i, t ) ] ϕ(, t) β ij = 4m im j (m i + m j ), G ij = 4πg ij (1 s) for i, j = 1,...,N. We omit the tildes and denote x = k m i in each of Eq. (.5). Then the resulting set of equations becomes t ϕ i(x, t) = N j=1 dsg ij (s)[ϕ i (x(1 β ij s))ϕ j (xβ ij s) ϕ i (x)ϕ j ()] m i m j β ij = 4 (m i + m j ) 1, β ij = β ji, β ii = 1, i, j = 1,...,N. Therefore, the most general system of isotropic Fourier transformed Boltzmann equations for Maxwell mixtures reads ϕ i t = N γ ij ϕ i ((1 β ij s)x)ϕ j (β ij sx) ϕ i (x)ϕ j (), (.6) j=1, for any function A(s), s [, 1], and A(s) = A(s)G(s)ds, G(s) = 4πg(1 s), γ ij = αij m ij, β ij = 4m ij m i m j, m ij = m im j m i + m j ; i, j = 1,...,N. (.7) 3. BINARY MIXTURE We consider a special case N = in Eqs. (.6), (.7) and denote m 1 = M, m = m, ϕ 1 = ρ 1 ϕ(x, t), ϕ = ρ ψ(x, t),

Bobylev and Gamba such that ϕ(, t) = ψ(, t) = 1. Then, we obtain ( ) 1/ ( ) 1/ α11 α1 ϕ t = ρ 1 ϕ,ϕ +ρ ϕ,ψ β M m 1 ( ) 1/ ( ) 1/ α α1 ψ t = ρ ψ, ψ +ρ 1 ψ, ϕ β (3.1) m m 1 β = 4mM (m + M), m 1 = mm m + M, ϕ,ψ β = following the notation of (.7), and thus, ds G(s){ϕ((1 βs)x, t)ψ(βsx, t) ϕ(x, t)ψ(, t)} ϕ,ϕ = ϕ,ϕ 1, ψ, ψ = ψ, ψ 1. We recall the connection of functions ϕ(x, t) and ψ(x, t) with corresponding solutions f 1, ( v, t) of the Boltzmann Eq. (.1) ( k ) ρ 1 ϕ M, t = dv f 1 ( v, t)e ik v, R 3 ( k ) ρ ψ m, t = dv f ( v, t)e ik v R 3 The usual definition of kinetic temperatures is given by equalities T i = m i dv v f i ( v, t), i = 1,. 3ρ i R 3 Then, one can easily verify that T 1 (t) = ϕ (, t), T (t) = ψ (, t). the primes denote derivatives on x. The equilibrium temperatures T eq of the binary mixture reads T eq = ρ 1T 1 + ρ T = const. ρ 1 + ρ The relaxation process in the binary mixture described by Eq. (3.1) leads to usual Maxwell asymptotics states ϕ t exp( T eq x), ψ t exp( T eq x).

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions By using Eq. (3.1), one can easily verify (at the formal level) that Therefore, dt 1 dt λ = = λρ (T 1 T ), ( α1 m 1 dt dt ) 1/ β s, s = = λρ 1 (T T 1 ), ds G(s)s. T 1 (t) = T eq + ρ 1 e t, T (t) = T eq + ρ e t, = ρ 1ρ (T 1 () T ()), = λ(ρ 1 + ρ ). (3.) ρ 1 + ρ It is easy to see that Eq. (3.1) formally admit the following class of self-similar solutions ϕ(x, t) = (x e t ) e T eq x, ψ(x, t) = (x e t )e T eq x It is, however, difficult to investigate such solutions (in particular, to prove that corresponding distribution functions are positive) in the most general case. Therefore we consider a simplified problem. 4. WEAKLY COUPLED BINARY MIXTURE If the masses M and m are fixed, then Eq. (3.1) contain five positive parameters p i,α ij, j = 1,. We shall consider below a special limiting case of Eq. (3.1) (weakly interacting gases) such that α 1, ρ, ρ α1 = const. (4.1) We assume that the other parameters ρ 1,α 11 and α remain constant and denote ϕ(x, t) = ϕ(x, t),ψ(x, t) = ψ(x, t), ( ) 1/ α11 t = ρ 1 t, (4.) M θ = ρ ( ) α1 M 1/ = const. ρ 1 α 11 m 1 Then we formally obtain (tildes are omitted below) ϕ t = ϕ,ϕ +θ ϕ,ψ β, ψ, ψ =, (4.3) it is assumed that the functions ϕ(x, t), ψ(x, t) and their time derivatives remain finite in the limit (4.1). The limiting temperatures (see Eq. (3.)) are given

by the identities Bobylev and Gamba T (t) = T () = T eq = const., and T 1 (t) = T () + (T 1 () T ())e θβ s t in the notation (4.) (tildes are omitted). Thus <ψ,ψ>=, ψ(, t) = 1, ψ (, t) = T () ψ(x, t) = e T ()x, and we reduce Eq. (4.3) to the unique equation for ϕ(x, t) ϕ = ϕ,ϕ +θ ϕ,e T()x β. (4.4) t The general case T () > can be reduced to the case T () = bysubstitution ϕ(x, t) = ϕ(x, t)exp( T ()x). Eq. (4.4) shows that, in the limiting case (4.1), the second component of the mixture plays a role of a thermostat with the fixed temperature T (), moreover, it is enough to consider the case T () = (cold thermostat). Then Eq. (4.4), in explicit form, reads ϕ t = dsg(s){ϕ(sx)ϕ[(1 s)x] + θϕ[(1 βs)x] ϕ(x)[ϕ() + θψ()]}, ϕ() = ψ() = 1, (4.5) the argument t of ϕ(x, t) is omitted. We consider below Eq. (4.5) assuming that ϕ( k, ) is a characteristic function (Fourier transform of a probability measure in R 3 ). Then Eq. (4.5) describes a homogeneous cooling process in the system of particles that interact between themselves and with the cold thermostat. Though all interactions are elastic (at the microlevel), such system has much in common with the gas of inelastic particles considered in. (3,4) It was proved in these papers that the general inelastic Maxwell model has the self-similar asymptotics in a certain precise sense. We conjecture the same asymptotic property for solutions of Eq. (4.5) (see (14) for its proof). In this paper we construct some explicit examples of self-similar solutions and show that such solutions have power-like tails for large velocities. 5. EXACT SOLUTIONS IN THE FOURIER-LAPLACE REPRESENTATION We consider Eq. (4.5) with β = 1 (equal masses m = M = 1) and G(s) = 1 (isotropic scattering). Then Eq. (4.5) reads ϕ t = dsϕ[(1 s)x][ϕ(sx) + θ] (1 + θ)ϕ(x), ϕ() = 1. (5.1)

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions From the physical point of view, this is a model for a mixture of two weakly interacting gases, which consist of particles with identical masses. One can assume that the sorts of particles differ, say, by color. Our goal is to describe a family of self-similar solutions of Eq. (5.1) such that ϕ(x, t) = ψ(xe µt ), ψ(x) = 1 ax p, x, (5.) The parameters θ, p > and µ R will be determined later. From now on x denotes the self-similar variable. (The notation ψ(xe µt ) should not be confused with one for the function ψ(x, t) from Sections 3 and 4.) Substituting Eq. (5.) into Eq. (5.1) we obtain to. () µxψ (x) (1 + θ)ψ(x) + 1 ψ (ψ + θ) =, x ψ 1 ψ = x dyψ 1 (y)ψ (x y). (5.3) This equation can be simplified by the use of the Laplace transform similarly w(z) = L(ψ)(z) = provided ψ(x) < A exp(z x), with some positive A and z. First, we recall properties of the Laplace transform and ψ(x)e zx, Re z > z, (5.4) L(xψ) = w (z), L(x ψ) = w (z) and L(ψ ) = zw(z) ψ(), L(x ψ ) = d dz (zw ψ()) = (zw(z)). Then we obtain the following equation for w(z): ( µ(zw) + (1 + θ)w + w w + θ ) =. z Next, we denote u(z) = zw(z) = so that the above equation is transformed to dx e x ψ ( ) x, (5.5) z µz u + (1 + θ)zu + u(u 1) =. (5.6) The next step is to simplify this equation by standard substitutions. We denote z = z q, u(z) = ũ( z), (5.7)

Bobylev and Gamba and obtain the equation for ũ( z) (tildes are omitted) µq z u + q[1 + θ + µ(q 1)]zu + u(u 1) =. Then, setting u(z) = z y(z) + B, B = const., (5.8) we obtain the following equation for y(z): µq z 4 y + z 4 y + αz 3 y + βz y + B(B 1) =, (5.9) the parameters α and β are given by the relations α = q(5µq + 1 + θ µ), β = B 1 + 4µq + q(1 + θ µ) Now, the parameters q and B can be chosen in such a way that α = β =. Thus, we take q = 1 + θ µ 5µ and obtain that Eq. (5.9) is reduced to, B = 6µq + 1, µq y + y B(B 1) + =. (5.1) z 4 This is the simplest standard form, for which Eq. (5.6) can be reduced in the general case. Finally, Eq. (5.1) is of the Painlevé type if and only if B = or B = 1. Otherwise, it has moving logarithmic singularities and does not have any simple nontrivial analytic solutions (this is just a repetition of arguments on a similar equation considered in () ). Therefore we consider the two special cases. In both cases Case 1: B =, which implies Case : B = 1, which implies θ = µ 1 5µq. (5.11) 6µq = 1, so that y = 6y. (5.1) 6µq = 1, so that y = 6y. (5.13) The general solution of the non-linear ODE y = by is expressed in terms of Weierstrass elliptic function. Similarly to, () we can show that just the simplest

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions exact solutions, namely Case 1: y(z) = (c 1 + z) and Case : y(z) = (c + z), with constant c 1 and c, are solutions to Eqs. (5.1) and (5.13) respectively, that lead to solutions of Eq. (5.6) satisfying appropriate boundary conditions at infinity. Coming back to the original notation (see the transformations (5.7), (5.8)), we obtain, therefore, two different exact solutions of Eq. (5.6) Case 1: u(z) = (1 + c 1 z q ) for 6µq = 1; Case : u(z) = 1 (1 + c z q ) for 6µq = 1. (5.14) We have that, in both cases, the coupling constant θ must satisfy (5.11). The constants c 1 and c are determined by the boundary conditions as follows. The given asymptotics in (5.) for ψ(x), x, leads to the asymptotics for u(z), as defined in (5.5), at infinity u(z) 1 + b z, z, p b is a non-zero constant. Recalling that ψ(x) 1for positive solutions of the Boltzmann equation, we can assume, without loss of generality, that u(z) 1 + 1, z. (5.15) z p Then we obtain, for the above two cases in (5.14), the following formulas satisfying the boundary condition (5.15) at infinity: Case 1: q = p; u(z) = θ = (3p 1)(1 p) p. ( 1 + 1 z p ), µ = 1 6p ; Case :q = p ; u(z) = 1 (1 + z p/ ), µ = 3p ; (3p + 1)( 1p) θ =. (5.16) 3p The result can be formulated in the following way. Proposition 5.1. condition Eq. (5.1) has exact self-similar solutions (5.) satisfying the ψ(x) 1 x p, Ɣ(p + 1) x, p >,

Bobylev and Gamba for the following values of the parameters θ(p) and µ(p): 1:µ(p) = 1 (3p 1)(1 p), θ(p) = ; 6p 6p :µ(p) = (3p + 1)( p), θ(p) = ; (5.17) 3p 3p The solutions of Eq. (5.1) are given by equalities [ ] ψ i (x) = L 1 ui (z), i = 1,, (5.18) z with u 1, (z) from Eq. (5.16), for cases 1 and respectively. ( ) The solutions have a physical meaning if, both, θ and ψ k is the Fourier transform of a positive function (measure). The first condition leads to inequalities 1 3 p 1 in the case 1, and to p in the case. The second condition will be discussed in the next section. 6. DISTRIBUTION FUNCTIONS First we evaluate the inverse Laplace transforms (5.18). In the case 1 we obtain [ ) ] ψ i (x) = L 1 1 (1 + z p 1, z 3 p 1. The general formula from (17) leads to [ ] 1 ψ 1 (x) = ( 1/p x), (x) = L 1 z (1 + z p ) sin pπ = pπ This case corresponds to θ = (1 + s cos pπ) ds e xs 1/p (1 + s + s cos pπ). (6.1) (3p 1)(1 p) 6p (6.) in Eq. (5.1). We note that θ = forp = 1 3, 1. Eq. (5.1) in such cases is the Fourier transformed Boltzmann equation for one-component gas. The exact solutions (6.1) with p = 1, 1 3 were already obtained in.() Eq. (6.1) therefore yields a generalization of these solutions to the case of binary mixture (with equal masses) provided the parameter θ is given (for given p [1/3, 1/]) in Eq. (6.). The corresponding

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions solutions of the Boltzmann equations are positive solutions with infinite energy. All their properties can be studied in the same way as in. () The case is more interesting since it includes also solutions with finite energy. The inverse Laplace transform ψ (x) = L 1 { 1 z [ 1 ]} 1 (1 + z p/ ) can be evaluated in the following way. We denote [ ] (x) = L 1 1, (1 + z p/ ) then x ψ (x) = 1 dy (y) = dy (y), dy (y) = 1. (6.3) x The function (x) can be expressed through the integral (x) = 1 πi C dz e xz (1 + z p/ ), < p the contour C lies around the negative half of the real axis (see (17) for details). Then we obtain (x) = 1 π dr e rx A(r), A(x) = Im [1 + (re iπ ) p/ ] = r p/ (1 + r p/ cos pπ )sin pπ (1 + r p + r p/ cos pπ ) Coming back to the function ψ (x) (6.3), we obtain ψ (x) = 1 dr e rx A(r). (6.4) π r The final result ψ (x) = 4sin pπ ds(1 + s cos pπ ) e xs/p pπ (1 + s + s cos pπ, < p, ) is obtained by substitution s = r p/ in the integral (6.4). We remind to the reader that the corresponding distribution function f ( v, t), that solves the Boltzmann equation, reads f ( v, t) = e 3µt/ F( v e µt/ ), µ = 3p, (6.5)

Bobylev and Gamba Noting that ( k F[F] = dvf( v )e ik v ) = ψ. R 3 v k T e T e = F[M T ( v )], M T ( v ) = (πt ), 3/ we obtain the integral representation of F( v ): F( v ) = 4sin pπ (1 + s cos pπ ds ) pπ (1 + s + s cos pπ M ) s/p( v ). (6.6) This function is obviously positive for < p 1. We note that the solution (6.5) (6.6) corresponds to the value in Eq. (5.1). θ = (3p + 1)( p) 3p 7. SOLUTIONS WITH FINITE ENERGY We consider in more detail the most interesting (for applications) case p = 1 in Eq. (6.6). Then We denote µ = 3, θ = 4 3, f ( v, t) = et F( v e t/3 ), (7.1) F( v ) = 4 π exp( v /s ) ds (πs ) 3/ (1 + s ). (7.) y = v, F( v ) = (π 5 ) 1/ (y), (7.3) r (y) = dr (1 + r) e ry = 1 r e r dr y 3 (1 + r. (7.4) ) y Asymptotic expansion of (y) for large positive y follows from integration of the formal series ( 1 + r ) ( ) r n = ( 1) n (n + 1). y y n=

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions Thus we obtain from Eq. (7.4) (y) 1 y 3 n= n (n + 1)(n + )! ( 1), y (7.5) y n In order to describe a behavior of (y) for small positive y we transform the first integral in Eq. (7.4) in the following way Noting that (y) = E m (y) = 1 1 dr (r 1) e (r 1)y = e y [E r (y) E 1 (y) + E (y)], dr e ry r m. we obtain E (y) = e y y, E (y) = e y ye 1 (y), E 1 (y) = We note that 1 y ds e s s (y) = 1 + 1 y ( + y)ey E 1 (y), ds e s s + = 1 ds e s s ds (e s 1) s = + y γ, 577 is the Euler constant. Therefore y ds (e s 1) s ln y. ds e s ln s = γ, E 1 (y) = (γ + ln y) + ds (1 e s ) s and (y) = 1 + 1 y ] [e y + ( + y) y (γ + ln y) ds (es 1). (7.6) s Thus, the asymptotic equality (7.5) and the formula (7.6) describe the behavior of the distribution function F(v) in (7.3), for large and small values of v.we obtain F( v ) = ( ) 5/ [ ( )] 1 1 1 + O, v, π v 6 v

Bobylev and Gamba F( v ) = 1/ 1 π 5/ v [1 + v ln v +O( v )], v. (7.7) All the exact solutions can be generalized to the case of the thermostat with finite temperature T (see Eq. (4.4) with T () = T ). Then Eq. (5.1) is replaced by the following equation ϕ t = + θ ϕ() = 1. ds{ϕ(sx)ϕ[(1 s)x] ϕ(x)ϕ()} ds{ϕ[(1 s)x]e Tsx ϕ(x)}, This equation can be reduced to Eq. (5.1) by substitution ϕ(x, t) = ˆϕ(x, t)e Tx. The corresponding self-similar solutions read ψ(x) satisfies Eq. (5.3). ϕ(x, t) = ψ(xe µt )e Tx, 8. SELF-SIMILAR SOLUTIONS AND POWER LIKE TAILS We consider in this section a more general class of equations for the function ϕ(x, t): ϕ t = + θ dsg(s){ϕ(a(s)x)ϕ[b(s)x] ϕ(x)ϕ()} dsh(s){ϕ[c(s)x] ϕ(x)}, ϕ() = 1, (8.1) with non-negative functions G(s), H(s), a(s), b(s) and c(s) with s [, 1]. We also assume that G(s), H(s) are integrable on [, 1], and a(s) 1, b(s) 1, c(s) 1; s 1. The function ϕ(x, t) is understood as the Fourier transform ϕ ( k, t) = dv f ( v, t) e ik v ; k R d, d = 1,,..., (8.) R d of a time dependent probability density f ( v, t) R d.

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions Eq. (8.1) allows to consider from a unified point of view two different kind of Maxwell models: I Inelastic Maxwell models, (1,3) H(s) =, a(s) = sz, b(s) = 1 sz( z), z = 1 + e, e 1. (8.3) II Maxwell mixtures described by Eq. (4.5), H(s) = θg(s), a(s) = s, b(s) = 1 s, c(s) = 1 βs. (8.4) The condition of integrability of G(s) and H(s) can be easily weakened. We do not do it here in order to simplify proofs. Our main goal in this section is to prove, roughly speaking, that self-similar solutions (distribution functions f ( v, t) from Eq. (8.)) have power-like tails. More precisely, we are going to show that such distribution functions can not have finite moments of any order. Eq. (8.1) admits (formally) a class of self-similar solutions ϕ(x, t) = ψ(xe µt ), µ >, (8.5) ϕ(x, t) satisfies the equation µxψ = + dsg(s)[ψ(a(s)x)ψ((b(s)x) ψ(x)] dsh(s)[ψ(c(s)x) ψ(x)], ψ() = 1, µ >. (8.6) The corresponding function f ( v, t) R d (see Eq. (8.)) reads f ( v, t) = e d µt µt F( v e ), ψ( k ) = dvf( v )e ik v. R d (8.7) Our goal is to prove the following general fact: if such a function F( v ) (generalized density of a probability measure in R d ) does exists, then it can not have finite moments of all orders. We assume the opposite and represent the integral in Eq. (8.1) as a formal series ψ( k ( 1) n ) = (n)! α n(d)m n k n, n= m n = dv F( v ) v n, R d n =, 1,...;

Bobylev and Gamba α n (1) = 1, α n (d) = 1 dω(ω ω) n, if d. S d 1 S d 1 S d 1 is the usual measure of the unit sphere S d 1 in R d,ω S d 1 is an arbitrary unit vector. Hence we obtain ψ( k ) = n= ( 1) n ψ n k n n!, ψ = 1, ψ n = n! (n!) α n(d) m n, n = 1,,... (8.8) The convergence of the Taylor series (8.8) is irrelevant for our goals. The only important point is that ψ(x) is infinitely differentiable for all x < (see any textbook in probability theory, for example (16) ) and ψ (n) () = ψ n >, n =, 1,... (8.9) On the other hand, the equations for ψ n can be easily obtained by substitution of the series (8.8) into Eq. (8.6). Then ψ = 1, ψ 1 [µ λ(1)] =, ψ n [µn λ(n)] + n 1 G(k, n k)ψ k ψ n k =, n =, 3,..., λ(n) = k=1 dsg(s)[1 a n (s) b n (s)] + dsh(s)[1 c n (s)], ( ) k + l 1 G(k, l) = dsg(s)a k (s)b l (s), k, l = 1,,... k Now we can use conditions (8.9). First we obtain µ = λ(1) and recall that µ> by assumption (8.5). Then we note that G(k, l) for all k, l = 1,,... Therefore ψ n [ µn + λ(n)] = λ(n) µ>, n =, 3,... n On the other hand, λ(n) ds[g(s) + H(s)] <, and, therefore, we get a contradiction. Thus, the following statement has been proven.

Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions Proposition 8.1. Eq. (8.6), µ>, G(s), H(s) L + [, 1], a(s) 1, b(s) 1, c(s) 1; s [, 1], does not have infinitely differentiable at x = solutions satisfying conditions (8.9). Corollary 8.. The corresponding probability density F( v ) from (8.7) cannot have finite moments of all orders. We can now apply the result to the inelastic Maxwell models (8.3) and conclude that the similar statement proved in our first paper (1) on that subject (see, (1) Section 5, Theorem 5.1) is valid not just for almost all, but for all values of the restitution coefficient e [, 1]. Consequently we can revise now a statement from (3) related to existence of self-similar solutions with finite moments of any order for a countable set of values of e from the interval set [, 1] (a possible logarithmic singularity was missing in the sketch of proof of Theorem 7. in (3) ). In fact the solution constructed in Ref. 3 has a finite number of moments for any e < 1, without exceptions. On the other hand, the above Proposition 8.1 can be applied to Maxwell mixtures (8.4). It shows that any physical (i.e. with a positive F( v ) ineq.(8.7)) solution of Eq. (8.6) corresponds to the distribution function F( v ) with a finite number of even integer moments. The exact solution constructed in Section 5 can be considered as a typical example. ACKNOWLEDGEMENTS The first author was supported by grant 61-3-5357 from the Swedish Research Council (NFR). The second author has been partially supported by NSF under grant DMS-4568. Support from the Institute for Computational Engineering and Sciences at the University of Texas at Austin is also gratefully acknowledged. REFERENCES 1. A. V. Bobylev, J. A. Carrillo and I. M. Gamba, On some properties of kinetic and hydrodynamic equations for inelastic interactions. J. Statist. Phys. 98(3 4):743 773 ().. A. V. Bobylev and C. Cercignani, Self-similar solutions of the Boltzmann equation and their applications. J. Statist. Phys. 16:139 171 (). 3. A. V. Bobylev and C. Cercignani, Self-similar asymptotics for the Boltzmann equation with inelastic and elastic interactions. J. Statist. Phys. 11:333 375 (3). 4. A. V. Bobylev, C. Cercignani and G. Toscani, Proof of an asymptotic property of self-similar solutions of the Boltzmann equation for granular materials. J. Statist. Phys. 111:43 417 (3). 5. A. V. Bobylev, The Fourier transform method for the Boltzmann equation for Maxwell molecules. Sov. Phys. Dokl. :8 8 (1976).

Bobylev and Gamba 6. A. V. Bobylev, The theory of the nonlinear spatially uniform Boltzmann equation for Maxwell molecules. Mathematical physics reviews, Vol. 7, 111 33, Soviet Sci. Rev. Sect. C Math. Phys. Rev., 7, Harwood Academic Publ., Chur, 1988. 7. A. Pulvirenti and B. Wennberg, A Maxwellian lower bound for solutions to the Boltzmann equation. Comm. Math. Phys. 183(1):145 16 (1997). 8. E. Gabetta, G. Toscani and B. Wennberg, Metrics for probability distributions and the trend to equilibrium for solutions of the Boltzmann equations. J. Statist. Phys. 81(5 6):91 934 (1995). 9. P. L. Krapivsky and E. Ben-Naim, Multiscaling in Infinite Dimensional Collision Processes. Phys. Rev. E 61: (). 1. T. Poschel and N. Brilliantov (Eds.), Granular Gas Dynamics, (Springer Berlin, 3). 11. M. H. Ernst and R. Brito, Scaling solutions of inelastic Boltzmann equations with overpopulated high energy tails. J. Stat. Phys. 19:47 43 (). 1. A. V. Bobylev, I. M. Gamba and V. A. Panferov, Moment inequalities and high-energy tails for Boltzmann equations with inelastic interactions. J. Statist. Phys. 116(5 6):1651 168 (4). 13. M. Krook and T. T. Wu, Exact solutions of the Boltzmann equation. Phys. Fl. :1589 1595 (1977). 14. A. V. Bobylev, C. Cercignani and I. M. Gamba, Generalized Maxwell models and self-similar asymptotics, to be published. 15. M. N. Kogan, Rarefied Gas Dynamics, (Plenum Press, 1969) 16. W. Feller, An Introduction to Probability Theory and Applications, Vol., (Wiley, N.Y., 1971). 17. A. V. Bobylev and C. Cercignani, The inverse Laplace transform of some analytic functions with an application to the eternal solutions of the Boltzmann equation. Appl. Math. Lett. 15(7):87 813 ().