# Low-Density Asymptotic Behavior of Observables of Hard Sphere Fluids.

1. Introduction

The main consistent approaches to the derivation of kinetic equations from underlying large particle dynamics were formulated by Bogolyubov [1] and Grad [2, 3]. For a hard sphere system Grad's method was developed by Cercignani [4] and Lanford [5, 6]. The rigorous results on the derivation of the Boltzmann equation with hard sphere collisions by methods of perturbation theory of the BBGKY hierarchy was proved in [7-10]. The most recent advances on the low-density (Boltzmann-Grad) scaling asymptotic behavior [11] of many-particle systems, in particular, systems with short-range interaction potentials, came in [12-24].

As is well known, many-particle systems are described by means of two objects: observables and states. A functional of the mean value of observables defines a duality between observables and states and as a consequence there exist two approaches to the description of the evolution within the framework of the evolution of observables and states, respectively [25]. Traditionally the evolution of many-particle systems is described within the framework of the evolution of states governed by the BBGKY hierarchy for marginal distribution functions. An equivalent approach to the description of the evolution of many-particle systems is given in terms of marginal observables governed by the dual BBGKY hierarchy [26].

The objective of the paper is to develop an approach to the description of the kinetic evolution of a hard sphere system within the framework of the evolution of observables. For this purpose in Section 2 we consider the microscopic description of the evolution of a hard sphere system within the framework of marginal observables governed by the dual BBGKY hierarchy. Then in Section 3 the origin of the dual kinetic evolution is stated; namely, a low-density (Boltzmann-Grad) limit of a nonperturbative solution of the Cauchy problem of the dual BBGKY hierarchy is established. In Sections 4 and 5 for initial states specified by means of a one-particle distribution function the link between the dual Boltzmann hierarchy for the limit marginal observables and the Boltzmann kinetic equation and the process of the propagation of initial chaos is established. In Sections 6 and 7 obtained results extended on hard spheres fluids, namely, for initial states specified by means of a one-particle distribution function and initial correlation functions, characterized condensed states. Finally, in Section 8 we conclude with some perspectives for future research.

2. The Dual BBGKY Hierarchy with Hard Sphere Collisions

As is well known, the evolution of many-particle systems can be described within the framework of a sequence of marginal (s-particle) distribution functions as well as in terms of a sequence of marginal observables. In this section we construct a nonperturbative solution of the Cauchy problem of a hierarchy of evolution equations for marginal observables of a hard sphere system.

We consider identical particles of a unit mass with a diameter [sigma] > 0, interacting as hard spheres with elastic collisions. Every particle is characterized by its phase coordinates ([q.sub.i], [p.sub.i]) = [x.sub.i] [member of] [R.sup.3] x [R.sup.3], i [greater than or equal to] 1. For configurations of such a system the following inequalities are satisfied: [absolute value of ([q.sub.i] - [q.sub.j])] [greater than or equal to] [sigma], i [not equal to] j [greater than or equal to] 1; that is, the set [W.sub.n] [equivalent to] {{[q.sub.1], ..., [q.sub.n]) [member of] [R.sup.3n] | [absolute value of ([q.sub.i] - [q.sub.j])] < [sigma] for at least one pair (i, j) : i [not equal to] j [member of] (1, ..., n)} is the set of forbidden configurations in the configuration space of n > 1 hard spheres. Let [C.sub.[gamma]] be the space of sequences b = ([b.sub.0], [b.sub.1], ..., [b.sub.n], ...) of bounded continuous functions on [R.sup.3n] x ([R.sup.3n] \ [W.sub.n]) which are symmetric with respect to permutations of the arguments [x.sub.1], ..., [x.sub.n], equal to zero on the set of forbidden configurations [W.sub.n] and equipped with the norm: [mathematical expression not reproducible], where 0 < [gamma] < 1.

If t [greater than or equal to] 0, the evolution of marginal observables B(t) = ([B.sub.0], [B.sub.1](t, [x.sub.1]), ..., [B.sub.s](t, [x.sub.1], ..., [x.sub.s]), ...) [member of] [C.sub.[gamma]] of a system of a nonfixed number of hard spheres is described by the Cauchy problem of the weak formulation of the following hierarchy of evolution equations [26]:

[mathematical expression not reproducible] (1)

[B.sub.s] (t, [x.sub.1], ..., [x.sub.s])[|.sub.t=0] = [B.sup.[epsilon],0.sub.2] ([x.sub.1], ..., [x.sub.s]), s [greater than or equal to] 1, (2)

where the coefficient [epsilon] > 0 is a scaling parameter (the ratio of the diameter [sigma] > 0 to the mean free path of hard spheres) and on the set [C.sub.s,0] [subset] [C.sub.s] of the continuously differentiable functions with compact supports the operators L(j) and [L.sub.int]([j.sub.1], [j.sub.2]) in a dimensionless form are defined by the formulas

[mathematical expression not reproducible] (3)

respectively. In (3) the symbol <x, x> denotes a scalar product, [delta] is the Dirac measure, [mathematical expression not reproducible], and the momenta [mathematical expression not reproducible], [mathematical expression not reproducible] are defined by the equalities

[mathematical expression not reproducible] (4)

We refer to recurrence evolution equations (1) as the dual BBGKY hierarchy for hard spheres in a dimensionless form. If t [less than or equal to] 0, a generator of the dual BBGKY hierarchy for hard spheres is defined by the expression of corresponding form [10].

To construct a solution of recurrence evolution equations (1) on the space [C.sub.n] = C([R.sup.3n] x ([R.sup.3n] \ [W.sub.n])) we introduce the group of operators [S.sub.n](t) that describes dynamics of n hard spheres. It is defined by means of the phase trajectories of a hard sphere system almost everywhere on the phase space [R.sup.3n] x ([R.sup.3n] \ [W.sub.n]), namely, beyond of the set [M.sup.0.sub.n] of the zero Lebesgue measure, as follows:

[mathematical expression not reproducible] (5)

where [X.sub.i](t) = [X.sub.i](t, [x.sub.1], ..., [X.sub.n]) is a phase trajectory of ith particle constructed in [7, 10], and the set [M.sup.0.sub.n] consists of the phase space points which are specified by initial data, generating multiple collisions of hard spheres in the evolutionary process, that is, collisions of more than two particles, more than one two-particle collision at the same instant, and infinite number of collisions on a finite time interval.

On the space [C.sub.n] one-parameter mapping (5) is an isometric *-weak continuous group of operators; that is, it is a [C.sup.*.sub.0]-group [27].

The infinitesimal generator [L.sub.n] of a group of operators (5) is defined in the sense of a *-weak convergence of the space [C.sub.n] and it has the structure [mathematical expression not reproducible], and the operators L(j) and [L.sub.int]([j.sub.1], [j.sub.2]) are defined by formulas (3).

A nonperturbative solution of the Cauchy problems (1) and (2) is determined by the following expansions:

[mathematical expression not reproducible] (6)

The generating operators of expansions (6) is the (1 + n)th-order cumulant of groups of operators (5) defined by the following expansion:

[mathematical expression not reproducible] (7)

where Y = (1, ..., s), Z = ([j.sub.1] ,..., [j.sub.n]) [subset] Y, the set, consisting of one element of the set of indices Y \ Z = (1, ..., [j.sub.1] - 1, [j.sub.1] + 1,..., jn - 1,jn + 1,..., s), we denoted by [Y \ Z}, the declusterization mapping [theta] is defined by the formula [theta]([Y \ Z}, Z) = Y, and the symbol [[summation].sub.p] means the sum over all possible partitions P of the set (1, ...,n) into |P| nonempty mutually disjoint subsets [X.sub.1] [subset] (1, ..., n).

The simplest examples of expansions (6) for marginal observables have the following form:

[mathematical expression not reproducible] (8)

On the space [C.sub.[gamma]] for the Cauchy problems (1) and (2) the following statement is true [28].

Theorem 1. For finite sequences of infinitely differentiable functions with compact supports B(0) = ([B.sub.0], [B.sup[epsilon],0.sub.1], ..., [B.sup.[epsilon],0.sub.s], ...) [member of] [C.sup.0.sub.[gamma]] [subset] [C.sub.[gamma]] a sequence of functions determined by expansions (6) is a classical solution and for arbitrary initial data B(0) [member of] [C.sub.[gamma]] it is a generalized solution.

Under the condition that [gamma] < [e.sup.-1], for a sequence of marginal observables (6), the estimate holds

[mathematical expression not reproducible] (9)

We remark that a one-component sequence of marginal observables corresponds to observables of certain structure; namely, the marginal observable [B.sup.(1)](0) = (0, [b.sup.[epsilon].sub.1] ([x.sub.1]), 0, ...) corresponds to the additive-type observable, and a one-component sequence of marginal observables [B.sup.(k)](0) = (0, ..., 0, [b.sup.[epsilon].sub.k] ([x.sub.1], ..., [x.sub.k]), 0, ...) corresponds to the k-ary-type observable [26]. If in capacity of initial data (2) we consider the additive-type marginal observables, then the structure of solution expansion (6) is simplified and attains the form

[mathematical expression not reproducible] (10)

where the generating operator of this expansion is the sth-order cumulant of groups of operators (7).

In the case of k-ary-type marginal observables solution expansion (6) has the form

[mathematical expression not reproducible] (11)

where the generating operator of this expansion is the (1 + s - k)th-order cumulant of groups of operators (7), and, if 1 [less than or equal to] s < k, we have [B.sup.(k).sub.s](t) = 0.

We remark also that expansion (6) can be also represented in the form of the perturbation (iteration) series [26] as a result of applying of analogs of the Duhamel equation to cumulants (7) of groups of operators (5).

Let Lin = Li([R.sup.3n] x ([R.sup.3n] \ [W.sub.n])) be the space of integrable functions that are symmetric with respect to permutations of the arguments [x.sub.1], ..., [x.sub.n], equal to zero on the set of forbidden configurations [W.sub.n] and equipped with the norm: [mathematical expression not reproducible]. A subspace of continuously differentiable functions with compact supports we denote by [L.sup.1.sub.n,0] [subset] [L.sup.1.sub.n].

The mean value of the marginal observable B(t) [member of] [C.sub.[gamma]] at t [member of] R is determined by the functional

[mathematical expression not reproducible] (12)

where initial state of finitely many hard spheres is described by means of a sequence of the marginal distribution functions F(0) = (1, [F.sup.[epsilon],0.sub.1], ..., [F.sup.[epsilon],0.sub.n], ...) [member of] [L.sup.1] = [[direct sum].sup.[infinity].sub.n=0] [L.sup.1.sub.n]. Owing to estimate (9), functional (12) exists under the condition that [gamma] < [e.sup.-1].

We remark that for mean value functional (12) the following equality holds:

(B(t), F(0)) = (B(0), F(t)), (13)

where the sequence F(t) = (1, [F.sub.1](t), ..., [F.sub.n](t), ...) [member of] [L.sup.1] = [[direct sum].sup.[infinity].sub.n=0] [L.sup.1.sub.n], is a solution of the BBGKY hierarchy for hard spheres. Generally such a solution is constructed by methods of perturbation theory [5-12, 29-32] (a nonperturbative solution was constructed in [33]). In case of infinitely many hard spheres [29, 30] a local in time solution of the Cauchy problem of the BBGKY hierarchy [7-12] is determined by perturbation series for arbitrary initial data from the space [L.sup.[infinity].sub.[xi]] of sequences of bounded functions equipped with the norm: [mathematical expression not reproducible]. In this case a local in time existence of the mean value functionals (B(0), F(t)) and (B(t), F(0)) was proved in papers [7], [10] and [34], [35], respectively.

3. The Kinetic Evolution of Hard Sphere Observables

We consider the problem of the rigorous description of the kinetic evolution of hard spheres within the framework of marginal observables by giving of a low-density (Boltzmann-Grad) asymptotic behavior of the Cauchy problem of the dual BBGKY hierarchy (1), (2).

Theorem 2. If for initial data [B.sup.[epsilon],0.sub.n] [member of] [C.sub.n], n [greater than or equal to] 1, there exists the limit [b.sup.0.sub.n] [member of] [C.sub.n]

[mathematical expression not reproducible], (14)

and then for arbitrary finite time interval there exists the Boltzmann-Grad limit of marginal observables (6) in the sense of a *-weak convergence of the space [C.sub.s]

[mathematical expression not reproducible], (15)

which is determined by the expansions

[mathematical expression not reproducible] (16)

where the operator [L.sup.0.sub.int]([j.sub.1], [j.sub.2]) is the collision operator of point particles, namely,

[mathematical expression not reproducible] (17)

Before proving this statement we give some comments.

We consider the Boltzmann-Grad limit of a special case of marginal observables, namely, the additive-type marginal observables. If for the initial additive-type marginal observable [b.sup.[epsilon].sub.1] the following condition is satisfied:

[mathematical expression not reproducible] (18)

then, according to statement (15), for additive-type marginal observables (10) we derive

[mathematical expression not reproducible] (19)

where the limit marginal observable [b.sup.(1).sub.s] (t) is determined as a special case of expansion (16):

[mathematical expression not reproducible] (20)

We make several examples of expansions (20) of the limit additive-type marginal observable:

[mathematical expression not reproducible] (21)

If for the initial k-ary-type marginal observable [b.sup.[epsilon].sub.k] the following condition is satisfied:

[mathematical expression not reproducible] (22)

then, according to statement (15), for k-ary-type marginal observables (11) we derive

[mathematical expression not reproducible] (23)

where the limit marginal observable [b.sup.(k).sub.s] (t) is determined as a special case of expansion (16):

[mathematical expression not reproducible] (24)

If [b.sup.0] [member of] [C.sub.[gamma]], then the sequence b(t) = ([b.sub.0], [b.sub.1](t), ..., [b.sub.s](t), ...) of limit marginal observables (16) is a generalized global solution of the Cauchy problem of the dual Boltzmann hierarchy with hard sphere collisions

[mathematical expression not reproducible] (25)

[b.sub.s] (t, [x.sub.1], ..., [x.sub.s])[|.sub.t=0] = [b.sup.0.sub.s] ([x.sub.1], ..., [x.sub.s]), s [greater than or equal to] 1, (26)

where the operator [L.sup.0.sub.nt] is defined by (17). This fact is proved similar to the case of an iteration series of the dual BBGKY hierarchy [26].

It should be noted that equations set (25) has the structure of recurrence evolution equations. Indeed, we make a few examples of the dual Boltzmann hierarchy with hard sphere collisions (25):

[mathematical expression not reproducible] (27)

Thus, in the Boltzmann-Grad scaling limit the kinetic evolution of hard spheres is described in terms of limit marginal observables (16) governed by the dual Boltzmann hierarchy with hard sphere collisions (25). Similar approach to the description of the mean field asymptotic behavior of quantum many-particle systems was developed in paper [36].

We outline the sketch of the proof of the limit theorem. For the group of operators (5) the analog of the Duhamel equation is valid [27]

[mathematical expression not reproducible] (28)

where the operator [L.sub.int](i, j) is defined by formula (3). Then for the (1 + n)th-order cumulant of groups of operators (5) the analog of the Duhamel equation holds

[mathematical expression not reproducible] (29)

where notations accepted above are used, [S.sub.s-n]([t.sub.n]) = [S.sub.s-n]([t.sub.n], (1, ..., s)\([j.sub.1], ..., [j.sub.n])), and we take into consideration the identity

[mathematical expression not reproducible] (30)

For arbitrary finite time interval *-weak continuous group of operators (5) has the following Boltzmann-Grad scaling limit in the sense of a *-weak convergence of the space [C.sub.s]

[mathematical expression not reproducible] (31)

Taking into account assumption (14) and an analog of the Duhamel equation (29), then in view of formula (31), for cumulants of asymptotically perturbed groups of operators we have

[mathematical expression not reproducible] (32)

As a result of equality (32) we establish the validity of statement (15) for nonperturbative solution (6) of the Cauchy problem of the dual BBGKY hierarchies (1) and (2).

4. The Derivation of the Boltzmann Kinetic Equation

We shall establish the link between the constructed asymptotic behavior of marginal observables of a hard sphere system (Theorem 2) and the description of kinetic evolution of states by means of a one-particle marginal distribution function governed by the Boltzmann kinetic equation.

In case of the absence of correlations between particles at initial time, that is, for initial states satisfying a chaos condition [10], the sequence of initial marginal distribution functions for a system of hard spheres has the form

[mathematical expression not reproducible] (33)

where [mathematical expression not reproducible] is the Heaviside step function of the allowed configurations. This assumption about initial state is intrinsic for the kinetic theory, because in this case all possible states of gases are described by means of a one-particle distribution function.

Let [F.sup.0,[epsilon].sub.1] [member of] [L.sup.[infinity].sub.[epsilon]] ([R.sup.3] x [R.sup.3]); that is, the inequality holds: [absolute value of ([F.sup.0,[epsilon].sub.1] ([x.sub.1]))] [less than or equal to] [xi] exp(- [beta]([p.sup.2.sub.i]/2)), where [xi] > 0, [beta] [greater than or equal to] 0 are parameters. We assume that the Boltzmann-Grad limit of the initial one-particle (marginal) distribution function [F.sup.0,[epsilon].sub.1] [member of] [L.sup.[infinity].sub.[xi]] ([R.sup.3] x [R.sup.3]) exists in the sense of a weak convergence of the space [L.sup.[infinity].sub.[xi]] ([R.sup.3] x [R.sup.3]), namely,

[mathematical expression not reproducible] (34)

and then the Boltzmann-Grad limit of the initial state (33) satisfies a chaos property too, that is, [f.sup.(c) [equivalent to] (1, [f.sup.0.sub.1] ([x.sub.1]), ..., [[PI].sup.s.sub.i=1] [f.sup.0.sub.1] ([x.sub.i), ...).

We note that assumption (34) with respect to the Boltzmann-Grad limit of initial states holds true for the equilibrium states [37].

If b(t) [member of] [C.sub.[gamma]] and [absolute value of ([f.sup.0.sub.1]([x.sub.i]))] [less than or equal to] [xi] exp(-[beta]([p.sup.2.sub.i]/2)), then the Boltzmann-Grad limit of mean value functional (12) exists under the condition that [7]: [mathematical expression not reproducible], and it is determined by the following series expansion:

[mathematical expression not reproducible] (35)

In consequence of the following equality for the limit additive-type marginal observables (20) (below it is proved in more general case)

[mathematical expression not reproducible] (36)

where function [b.sup.(1).sub.s](t) is given by expansion (20) and the distribution function [f.sub.1] (t, [x.sub.1]) is given by the series

[mathematical expression not reproducible] (37)

where the following operator was introduced:

[mathematical expression not reproducible] (38)

and the group of operators [S.sup.*.sub.1] (t) is a group of adjoint operators to operators (5) in the sense of mean value functional (12).

The distribution function [f.sub.1] (t) is a solution of the Cauchy problem of the Boltzmann kinetic equation

[mathematical expression not reproducible] (39)

[f.sub.1] (t, [x.sub.1])[|.sub.t=0] = [f.sup.0.sub.1] ([x.sub.1]). (40)

Thus, we establish that hierarchy (25) for additive-type marginal observables and initial state (34) describes the evolution of hard sphere systems just as the Boltzmann kinetic equation (39).

We differentiate over the time variable expression (37) in the sense of the pointwise convergence of the space [L.sup.[infinity].sub.[xi]] ([R.sup.3] x [R.sup.3])

[mathematical expression not reproducible] (41)

where the operator [L.sup.0,*.sub.int] ([k.sub.i], i + 2) is defined by formula (38).

Using the product formula for the one-particle marginal distribution function [f.sub.1](t, [x.sub.i]) defined by series expansion (37) in case of initial data (34)

[mathematical expression not reproducible] (42)

where the group property of one-parameter mapping (5) is applied, we express the second summand in the right-hand side of equality (41) in terms of [[PI].sup.2.sub.i=1] [f.sub.1](t, i), and, consequently, we get (39).

We remark that in a one-dimensional space the collision integral of the Boltzmann equation with elastic hard sphere collisions identically equals zero. In a one-dimensional space the Boltzmann-Grad limit is not trivial in case of hard sphere dynamics with inelastic collisions [38]. In paper [38] for one-dimensional granular gas the process of the creation of correlations in the Boltzmann-Grad limit was also described.

5. On Propagation of Initial Chaos in a Low-Density Limit

If the initial states of hard spheres are specified by a sequence of marginal distribution functions (33), then the property of the propagation of initial chaos holds in the Boltzmann-Grad limit. It is a result of the validity of the following equality for the limit k-ary marginal observables (24); that is, [b.sup.(k)](0) = (0, ..., [b.sup.0.sub.k]([x.sub.1], ..., [x.sub.k]), 0, ...),

[mathematical expression not reproducible] (43) where for finite time interval the limit one-particle marginal distribution function [f.sub.1](t) is defined by series expansion (37) and therefore it is governed by the Cauchy problem of the Boltzmann kinetic equations (39) and (40).

In fact, taking into account the validity of the following equality for expansion (16) of the function [b.sup.(k).sub.s](t)

[mathematical expression not reproducible] (44)

and product formula (42), for the limit one-particle marginal distribution function defined by series expansion (37), we finally verify the validity of equality (43).

Thus, in the Boltzmann-Grad scaling limit an equivalent approach to the description of the kinetic evolution of hard spheres within the framework of the Cauchy problem of the Boltzmann kinetic equations (39) and (40) is given by the Cauchy problem of the dual Boltzmann hierarchy with hard sphere collisions (25) and (26) for the additive-type marginal observables. In case of the nonadditive-type marginal observables a solution of the dual Boltzmann hierarchy with hard sphere collisions (25) is equivalent to the property of a propagation of initial chaos in the sense of equality (43).

6. The Boltzmann Equation for Hard Spheres Fluids

We consider initial states of a hard sphere system specified by the one-particle marginal distribution function [F.sup.0,[epsilon].sub.1] [member of] [L.sup.[infinity].sub.[xi]] ([R.sup.3] x [R.sup.3]) in the presence of correlations, that is, initial states defined by the following sequence of marginal distribution functions:

[mathematical expression not reproducible] (45)

where the functions [g.sup.[epsilon].sub.n] = [g.sup.[epsilon].sub.n] ([x.sub.1], ..., [x.sub.n]) [member of] [C.sub.n]([R.sup.3n] x ([R.sup.3n] \ [W.sub.n])), n [greater than or equal to] 2, are specified initial correlations. Since many-particle systems in condensed states are characterized by correlations sequence (45) describes the initial state of the kinetic evolution of hard sphere fluids.

We assume that the Boltzmann-Grad limit of initial one-particle marginal distribution function [F.sup.0,[epsilon].sub.1] [member of] [L.sup.[infinity].sub.[xi]] ([R.sup.3] x [R.sup.3]) exists in the sense as above; that is, in the sense of a weak convergence the equality holds: w - [lim.sub.[epsilon][right arrow]0]([[epsilon].sup.2] [F.sup.0,[epsilon].sub.1] - [f.sup.0.sub.1]) = 0, and in case of correlation functions let w - [lim.sub.[epsilon][right arrow]0]([g.sup.[epsilon].sub.n] - [g.sub.n]) = 0, n [greater than or equal to] 2; then in the Boltzmann-Grad limit initial state (45) is defined by the following sequence of the limit marginal distribution functions:

[mathematical expression not reproducible] (46)

We consider relationships of the constructed Boltzmann-Grad asymptotic behavior of marginal observables with the nonlinear Boltzmann-type kinetic equation in case of initial states (46).

For the limit additive-type marginal observables (20) and initial states (46) the following equality is true:

[mathematical expression not reproducible] (47)

where the functions [b.sup.(1).sub.s](t) are represented by expansions (20) and the limit marginal distribution function [f.sub.1](t) is represented by the following series expansion:

[mathematical expression not reproducible] (48)

Series (48) is uniformly convergent for finite time interval under the condition as above (37).

The function [f.sub.1](t) represented by series (48) is a weak solution of the following Cauchy problem of the Boltzmann kinetic equation with initial correlations [39, 40]

[mathematical expression not reproducible] (50)

This fact is proved similarly to the case of a perturbative solution of the BBGKY hierarchy for hard spheres represented by the iteration series [10, 29].

Thus, in case of initial states specified by one-particle marginal distribution function (46) we establish that the dual Boltzmann hierarchy with hard sphere collisions (25) for additive-type marginal observables describes the evolution of a hard sphere system just as the Boltzmann kinetic equation with initial correlations (49).

7. Propagation of Initial Correlations in a Low-Density Limit

The property of the propagation of initial correlations in a low-density limit is a consequence of the validity of the following equality for a mean value functional of the limit k-ary marginal observables:

[mathematical expression not reproducible] (51)

where the one-particle marginal distribution function [f.sub.1](t, [x.sub.j]) is solution (48) of the Cauchy problem of the Boltzmann kinetic equation with initial correlations (49) and (50), and the inverse group to the group of operators [S.sup.*.sub.1](t) we denote by [([S.sup.*.sub.1]).sup.-1](t) = [S.sup.*.sub.1](-t) = [S.sub.1](t).

This fact is proved similarly to the proof of a property of a propagation of initial chaos (43).

We note that, according to equality (51), in the Boltzmann-Grad limit the marginal correlation functions defined as cluster expansions of marginal distribution functions, namely,

[mathematical expression not reproducible] (52)

have the following explicit form:

[mathematical expression not reproducible] (53)

where for initial correlation functions (46) it is used the following notations:

[mathematical expression not reproducible] (54)

where the symbol [[summation].sub.p] means the sum over possible partitions P of the set of arguments ([x.sub.1], ..., [x.sub.s]) on [absolute value of (P)] nonempty subsets [X.sub.i], and the one-particle marginal distribution function [f.sub.1](t) is a solution of the Cauchy problem of the Boltzmann kinetic equation with initial correlations (49) and (50).

Thus, in case of the limit k-ary marginal observables a solution of the dual Boltzmann hierarchy with hard sphere collisions (25) is equivalent to a property of the propagation of initial correlations for the fc-particle marginal distribution function in the sense of equality (51) or in other words the Boltzmann-Grad scaling dynamics does not create new correlations except initial correlations.

8. Conclusion

In the paper a new approach to the problem of the rigorous description of the kinetic evolution of a system of hard spheres with elastic collisions was developed. For this purpose we established the low-density (Boltzmann-Grad) asymptotic behavior of a solution of the Cauchy problem of the dual BBGKY hierarchy for marginal observables of hard spheres (1) and (2). The constructed scaling limit is governed by the set of recurrence evolution equations, namely, by the dual Boltzmann hierarchy with hard sphere collisions (25).

Furthermore, it was established that for initial states specified by a one-particle distribution function the evolution of additive-type marginal observables is equivalent to a solution of the Boltzmann kinetic equation (39) and the evolution of nonadditive-type marginal observables is equivalent to the property of the propagation of initial chaos for states (43). In other words the Boltzmann-Grad dynamics does not create correlations.

One of the advantages of such an approach to the derivation of the Boltzmann equation is an opportunity to construct the kinetic equation, involving correlations at initial time, in particular, that can characterize the condensed states. Moreover, it gives opportunity to describe the propagation of initial correlations in the Boltzmann-Grad scaling limit (53).

Some applications ofthe developed method to the derivation of kinetic equations in scaling limits of large particle systems of different kinds, in particular, hard spheres with inelastic collisions [38], are considered in papers [38, 41, 42].

We note that one more approach to the description of the kinetic evolution of hard spheres is based on the non-Markovian generalization of the Enskog kinetic equation [43].

https://doi.org/10.1155/2018/6252919

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

[1] N. N. Bogolyubov, Problems of Dynamic Theory in Statistical Physics, Gostekhizdat, 1946.

[2] H. Grad, "On the kinetic theory of rarefied gases," Communications on Pure and Applied Mathematics, vol. 2, pp. 331-407, 1949.

[3] H. Grad, "Principles of the kinetic theory of gases," in Handbuch derPhysik, vol. 12, pp. 205-294, Springer, Berlin, Germany, 1958.

[4] C. Cercignani, "On the Boltzmann equation for rigid spheres," Transport Theory and Statistical Physics, vol. 2, no. 3, pp. 211-225, 1972.

[5] O. E. Lanford, "Time evolution of large classical systems," in Dynamical Systems, Theory and Applications, vol. 38 of Lecture Notes in Physics, pp. 1-111, Springer, Berlin, Germany, 1975.

[6] O. E. Lanford, "The hard sphere gas in the Boltzmann-Grad limit," Physica A: Statistical Mechanics and its Applications, vol. 106, no. 1-2, pp. 70-76, 1981.

[7] D. Y. Petrina and V I. Gerasimenko, "Mathematical problems of the statistical mechanics of a hard-sphere system," Russian Mathematical Surveys (Uspekhi Matematicheskikh Nauk), vol. 45, no. 3, pp. 135-182, 1990.

[8] H. Spohn, Large Scale Dynamics of Interacting Particles, Springer-Verlag, Berlin, Germany, 1991.

[9] C. Cercignani, R. Illner, and M. Pulvirenti, The Mathematical Theory of Dilute Gases, Springer-Verlag, Berlin, Germany, 1994.

[10] C. Cercignani, V. I. Gerasimenko, and D. Y. Petrina, Many-particle Dynamics and Kinetic Equations, Springer, Berlin, Germany, 2012.

[11] H. Spohn, "Kinetic equations from Hamiltonian dynamics: Markovian limits," Reviews of Modern Physics, vol. 52, no. 3, pp. 569-615, 1980.

[12] I. Gallagher, L. Saint-Raymond, and B. Texier, From Newton to Boltzmann: Hard Spheres and Short-range Potentials, Zurich Lectures in Advanced Mathematics, European Mathematical Society, Helsinki, Finland, 2014.

[13] T. Bodineau, I. Gallagher, and L. Saint-Raymond, "From hard sphere dynamics to the Stokes-Fourier equations: an [L.sup.2] analysis of the Boltzmann-Grad limit," Annals of PDE. Journal Dedicated to the Analysis of Problems from Physical Sciences, vol. 3, no. 2, 2017.

[14] M. Pulvirenti, C. Saffirio, and S. Simonella, "On the validity of the Boltzmann equation for short range potentials," Reviews in Mathematical Physics, vol. 26, no. 2, Article ID 1450001, 2014.

[15] S. Simonella, "Evolution of correlation functions in the hard sphere dynamics," Journal of Statistical Physics, vol. 155, no. 6, pp. 1191-1221, 2014.

[16] M. Pulvirenti and S. Simonella, "The Boltzmann-Grad limit of a hard sphere system: analysis of the correlation error," Inventiones Mathematicae, vol. 207, no. 3, pp. 1135-1237, 2017.

[17] M. Pulvirenti and S. Simonella, "On the evolution of the empirical measure for the hard-sphere dynamics," Bulletin of the Institute of Mathematics, Academia Sinica, vol. 10, no. 2, pp. 171-204, 2015.

[18] K. Aoki, M. Pulvirenti, S. Simonella, and T. Tsuji, "Backward clusters, hierarchy and wild sums for a hard sphere system in a low-density regime," Mathematical Models and Methods in Applied Sciences, vol. 25, no. 5, pp. 995-1010, 2015.

[19] R. I. Patterson, S. Simonella, and W. Wagner, "Kinetic theory of cluster dynamics," Physica D: Nonlinear Phenomena, vol. 335, pp. 26-32, 2016.

[20] Th. Paul, M. Pulvirenti, and S. Simonella, "On the size of chaos in the mean field dynamics," https://arxiv.org/abs/1708.07701.

[21] M. Tessarotto, "On the BoltzmannGrad limit for the classical hard-spheres systems," in Proceedings of the Proc. 26 Int. Symposium on Rarefied Gas Dynamics, pp. 194-199, AIP, NY, USA, 2008.

[22] R. Denlinger, "The propagation of chaos for a rarefied gas of hard spheres in the whole space," https://arxiv.org/abs/ 1605.00589.

[23] R. Denlinger, Observables and strong one-sided chaos in the Boltzmann-Grad limit, 2017, https://arxiv.org/abs/1706.09961.

[24] C. Saffirio, "Derivation of the Boltzmann equation: hard spheres, short-range potentials and beyond," in From Particle Systems to Partial Differential Equations III, vol. 162 of Springer Proceedings in Mathematics & Statistics, pp. 301-321, 2016.

[25] V. I. Gerasimenko, "Hierarchies of quantum evolution equations and dynamics of many-particle correlations," in in Statistical Mechanics and Random Walks: Principles, Processes and Applications, pp. 233-288, Nova Science Publishers, Hauppauge, NY, USA, 2012.

[26] G. Borgioli and V. I. Gerasimenko, "The dual BBGKY hierarchy for the evolution of observables," Rivista di Matematica della Universita di Parma, vol. 4, pp. 251-267, 2001.

[27] J. Banasiak and L. Arlotti, Perturbations of Positive Semi-groups with Applications, Springer-Verlag, Berlin, Germany, 2006.

[28] V. I. Gerasimenko, "On the approaches to the derivation of the Boltzmann equation with hard sphere collisions," Proc. Inst. Math. NASU, vol. 10, no. 2, pp. 71-95, 2013.

[29] V. I. Gerasimenko and D. Ya. Petrina, "Thermodynamical limit of nonequilibrium states of three-dimensional hard-sphere system," Theoretical and Mathematical Physics, vol. 64, no. 1, pp. 130-149, 1985.

[30] V. I. Gerasimenko, "On a solution of the BBGKY hierarchy for a one-dimensional hard sphere system," Theoretical and Mathematical Physics, vol. 91, no. 1, pp. 120-132, 1992.

[31] D. Y. Petrina, Stochastic Dynamics and Boltzmann Hierarchy, NASU Institute of Mathematics, Kiev, Ukraine, 2008.

[32] H. Spohn, "On the integrated form of the BBGKY hierarchy for hard spheres., arXiv," Article ID 0605068, https://arxiv.org/abs/math-ph/0605068.

[33] V. I. Gerasimenko, T. V. Ryabukha, and M. O. Stashenko, "On the structure of expansions for the BBGKY hierarchy solutions," Journal of Physics A: Mathematical and General, vol. 37, no. 42, pp. 9861-9872, 2004.

[34] V. I. Gerasimenko, Y. G. Kondratiev, T. Kuna, and D. Y. Petrina, "Evolution of observables and quasiobservables in classical statistical mechanics," Journal of Physics: Condensed Matter, vol. 3, no. 2, pp. 239-264, 2000.

[35] T. V. Ryabukha, "On mean-value ofobservalesfunctionals of one-dimensional infinite particle systems," Theoretical and Mathematical Physics, vol. 162, no. 3, pp. 422-438, 2010.

[36] V. Gerasimenko, "Heisenberg picture of quantum kinetic evolution in mean-field limit," Kinetic and Related Models, vol. 4, no. 1, pp. 385-399, 2011.

[37] V. I. Gerasimenko and D. Y. Petrina, "The Boltzmann-Grad limit of equilibrium states," Reports of the National Academy of Sciences of Ukraine SSR, vol. 12, pp. 17-19, 1988.

[38] M. S. Borovchenkova and V I. Gerasimenko, "On the non-Markovian Enskog equation for granular gases," Journal of Physics A: Mathematical and Theoretical, vol. 47, no. 3, Article ID 035001, 2014.

[39] I. V. Gapyak, "The kinetic equations of a hard sphere system with initial correlations," Proc. Inst. Math. NASU, vol. 11, no. 1, pp. 166-177, 2014.

[40] V. Gerasimenko and A. Kornienko, "The Boltzmann kinetic equation with correlations for hard sphere fluids," Reports of the National Academy of Sciences of Ukraine, no. 3, pp. 17-23, 2015.

[41] V. I. Gerasimenko, "New approach to derivation of quantum kinetic equations with initial correlations," Carpathian Journal of Mathematics, vol. 7, no. 1, pp. 38-48, 2015.

[42] V. I. Gerasimenko and Yu. Yu. Fedchun, "On semigroups of large particle systems and their scaling asymptotic behavior," in Semigroups of Operators - Theory and Applications, vol. 113 of Springer Proceedings in Mathematics & Statistics, pp. 165-182, Springer, Berlin, Germany, 2015.

[43] V. I. Gerasimenko and I. V. Gapyak, "Hard sphere dynamics and the Enskog equation," Kinetic and Related Models, vol. 5, no. 3, pp. 459-484, 2012.

Viktor Gerasimenko (iD) (1) and Igor Gapyak (2)

(1) Institute of Mathematics of the NAS of Ukraine, No. 3, Tereshchenkivs'ka Str., Kyiv 01004, Ukraine

(2) Department of Mechanics and Mathematics, Taras Shevchenko National University of Kyiv, No. 2, Academician Glushkov Av., Kyiv 03187, Ukraine

Correspondence should be addressed to Viktor Gerasimenko; gerasym@imath.kiev.ua

Received 15 November 2017; Revised 12 January 2018; Accepted 4 February 2018; Published 8 March 2018