Next Article in Journal
Entropy Generation of Desalination Powered by Variable Temperature Waste Heat
Next Article in Special Issue
From Lattice Boltzmann Method to Lattice Boltzmann Flux Solver
Previous Article in Journal
Analytical Solutions of the Black–Scholes Pricing Model for European Option Valuation via a Projected Differential Transformation Method
Previous Article in Special Issue
Extension of the Improved Bounce-Back Scheme for Electrokinetic Flow in the Lattice Boltzmann Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Truncation Scheme for the BBGKY2 Equation

by
Gregor Chliamovitch
1,2,*,
Orestis Malaspinas
1 and
Bastien Chopard
1
1
Department of Computer Science, University of Geneva, Route de Drize 7, 1227 Geneva, Switzerland
2
Department of Theoretical Physics, University of Geneva, Quai Ernest-Ansermet 24, 1211 Geneva, Switzerland
*
Author to whom correspondence should be addressed.
Entropy 2015, 17(11), 7522-7529; https://0-doi-org.brum.beds.ac.uk/10.3390/e17117522
Submission received: 28 September 2015 / Revised: 23 October 2015 / Accepted: 26 October 2015 / Published: 30 October 2015
(This article belongs to the Special Issue Non-Linear Lattice)

Abstract

:
In recent years, the maximum entropy principle has been applied to a wide range of different fields, often successfully. While these works are usually focussed on cross-disciplinary applications, the point of this letter is instead to reconsider a fundamental point of kinetic theory. Namely, we shall re-examine the Stosszahlansatz leading to the irreversible Boltzmann equation at the light of the MaxEnt principle. We assert that this way of thinking allows to move one step further than the factorization hypothesis and provides a coherent—though implicit—closure scheme for the two-particle distribution function. Such higher-order dependences are believed to open the way to a deeper understanding of fluctuating phenomena.

1. Introduction

While the formulation of equilibrium statistical mechanics in terms of the maximum entropy (“MaxEnt”) principle goes back to Jaynes’ seminal work [1,2] in the 50s, the last decade has seen a spectacular revival of this approach. In particular, the MaxEnt-based characterization of complex systems presented in [3] paved the way to applications in a variety of fields ranging from linguistics to biology [4,5,6,7,8]. First focussed on equilibrium situations, these works soon turned their attention to non-equilibrium properties as well [9,10,11,12].
However, it seems that these authors paid comparatively little attention to more “fundamental” issues. In the present letter, we would like to reconsider some aspects of the kinetic theory of gases at the light of the MaxEnt philosophy. More precisely, kinetic theory relies heavily on the so-called Stosszahlansatz which asserts that, before colliding, particles are uncorrelated. While this bold assumption can be motivated physically, it is not completely clear how it should be generalized when considering higher-order distribution functions. Our point is that if one considers the Stosszahlansatz as a heuristically motivated assumption, it generalizes naturally to higher-order cases—even though this raises extra mathematical challenges!
We start with a brief reminder on MaxEnt distributions as well as on the BBGKY hierarchy leading to kinetic equations, and bridge both in the last two sections.

2. Maximum Entropy Distributions

The Shannon entropy H ( X ) = d x P ( x ) ln P ( x ) has all properties one would expect from an uncertainty measure [13,14]. In other terms, among a set of distributions, the least biased guess an observer can make is the one having the largest entropy while still satisfying available observational constraints. Assume for instance we try to maximize H ( X ) under the constraint O = d x O ( x ) P ( x ) = μ . Introducing a multiplier for the constraint and another for the probabilistic normalization, one has to compute
P d x P ( x ) ln 1 P ( x ) + λ 0 d x P ( x ) + λ d x O ( x ) P ( x ) = d x ln 1 P ( x ) 1 + λ 0 + λ O ( x ) .
Letting this expression vanish yields
P ( x ) = exp 1 + λ 0 + λ O ( x ) ,
where the multipliers have to be determined to match the constraints. This result extends straightforwardly to the case of several constraints, namely O k = d x O k ( x ) P ( x ) = μ k for k = 1 , 2 , , K . Then
P ( x ) = 1 Z exp k = 1 K λ k O k ( x ) .
In what follows we shall be concerned primarily with constraints on marginals instead of averages. An appropriate use of δ functions allows to generalize the previous result. As an example consider the case of a quadrivariate variable x = ( w , x , y , z ) the marginal P 123 ( a , b , c ) of which is assumed to be known. Putting O ( x ) = δ ( w , a ) δ ( x , b ) δ ( y , c ) we can write
O = x O ( x ) P ( x ) = w , x , y δ ( w , a ) δ ( x , b ) δ ( y , c ) z P ( x ) = w , x , y δ ( w , a ) δ ( x , b ) δ ( y , c ) P 123 ( w , x , y ) = P 123 ( a , b , c ) .
Using the result (3) derived for the averages one gets
P ( x ) = 1 Z exp a , b , c λ ( a , b , c ) O ( x ) = 1 Z exp λ ( w , x , y ) ,
λ now denoting a function. This result can be extended to any number of marginals of any order. If for instance besides P 123 the marginals P 124 and P 34 are given we get
P ( x ) = 1 Z exp λ 1 ( w , x , y ) + λ 2 ( w , x , z ) + λ 3 ( y , z )
for functions λ 1 , λ 2 and λ 3 to be determined. Unfortunately, this determination is difficult except in the case of univariate constraints. Then λ 1 ( w ) = ln P 1 ( w ) , etc., obviously solves the problem, so that the corresponding MaxEnt distribution is the factorized distribution. When turning our attention to the applicability of this result to the realm of kinetic theory, it will appear that in that context the problem can be slightly simplified due to the structure of reduced distributions.

3. The BBGKY Hierarchy

Let us consider N particles of mass m, the coordinates of which in phase space are their positions q i and momenta p i . It will be convenient to define a condensed notation x i = ( q i , p i ) . Let f N ( x 1 , , x N , t ) denote the joint distribution function characterizing the system, which obeys Liouville’s equation
d f N d t = f N t + i = 1 N p i m f N q i + i = 1 N F i f N p i = 0 ,
where F i denotes the force exerted on particle i. We shall restrict ourselves to the case without external force and where particles interact pairwise through some radial potential V ( | q i q j | ) = V i j . Then F i = j i V i j q i .
Reminding that f N itself is normalized to N ! , we now introduce the reduced s-particle distribution f s ( x 1 , , x s , t ) = N ! ( N s ) ! d x s + 1 d x N f N ( x 1 , , x N , t ) . The standard result [15] is that by integrating Liouville’s equation, one obtains a dynamical equation for f s given by the so-called BBGKY hierarchy (from the non-chronological list of its co-discoverers’ names : Bogoliubov, Born, Green, Kirkwood, Yvon) :
f s t + i = 1 s p i m f s q i i = 1 s j i s V i j q i f s p i d x s + 1 i = 1 s V i , s + 1 q i f s + 1 p i = 0 .
This expression forms a hierarchy since the dynamics for f s is expressed in terms of the higher-order distribution f s + 1 . Of course each equation can be deduced from its higher-order precursor by integration, at the cost of an information loss. In what follows we shall denote the s-th equation of the hierarchy as BBGKY-s. BBGKY1 and BBGKY2, in which we are primarily interested here, are
f 1 t + p 1 m f 1 q 1 = d x 2 V 12 q 1 f 2 p 1
and
f 2 t + p 1 m f 2 q 1 + p 2 m f 2 q 2 V 12 q 1 p 1 p 2 f 2 = d x 3 V 13 q 1 f 3 p 1 + V 23 q 2 f 3 p 2 ,
where of course f 1 = f 1 ( p 1 , q 1 , t ) , f 2 = f 2 ( p 1 , q 1 , p 2 , q 2 , t ) and f 3 = f 3 ( p 1 , q 1 , p 2 , q 2 , p 3 , q 3 , t ) . The purpose of this paper is to investigate the second of these equations. As stressed above, we shall not try to express BBGKY1 and BBGKY2 as a set of coupled equations relating f 1 and f 2 , since such an approach would not “fit” nicely in the spirit of the BBGKY approach. Instead, we shall manage to truncate BBGKY2 in order to obtain a single, self-standing equation for f 2 .

4. The Stosszahlansatz for BBGKY2

The procedure leading from the BBGKY1 equation to a consistent kinetic equation for f 1 is standard: the Stosszahlansatz asserts that before colliding two particles are uncorrelated, i.e., f 2 factorizes as f 2 ( x 1 , x 2 ) = f 1 ( x 1 ) f 1 ( x 2 ) . This allows us to express the collision integral in terms of f 1 , so that BBGKY1 becomes a closed equation for f 1 . Since this factorization hypothesis may be supported from a physical standpoint, it is tempting to use this ansatz in the collision term for BBGKY2 as well. But this raises an issue: if BBGKY2 can be cast into an equation relating a streaming term expressed in terms of f 2 to a collision term expressed in terms of f 1 , then this equation is obviously not consistent by itself and has to be supplemented, so as to obtain a system of coupled equations.
Our point is that this issue vanishes if the Stosszahlansatz is reconsidered as a heuristic ansatz instead of a physically-grounded assumption. We propose to reformulate it as follows: since the exact codependence of particles entering the collision range is unknown, one has to make a reasonable guess on it, and the MaxEnt distribution steps out at this point. The MaxEnt guess for f 2 , compatible with the univariate distribution appearing in the streaming term, is the factorized one, but on the contrary the guess for f 3 , compatible with the f 2 appearing in the left-hand side, will be quite different from a factorized distribution (as exemplified by (6)).
Let us now see to what extent the result (6) obtained for the MaxEnt distribution may be particularized to our current purpose. We showed that, given bivariate marginals, the MaxEnt estimate for f 3 ( x 1 , x 2 , x 3 ) was given by
f 3 M E ( x 1 , x 2 , x 3 ) = 1 Z exp λ 1 ( x 1 , x 2 ) + λ 2 ( x 1 , x 3 ) + λ 3 ( x 2 , x 3 )
for some λ 1 , λ 2 and λ 3 . The point is that these marginals are the same for each pair by definition of the reduced distribution f 2 , and accordingly all three λs are the same. Absorbing the normalization, one is therefore allowed to write that
f 3 M E ( x 1 , x 2 , x 3 ) = g ( x 1 , x 2 ) g ( x 1 , x 3 ) g ( x 2 , x 3 )
for a function g that is nevertheless unknown, except for the fact that it has to satisfy the marginal constraint
g ( x , y ) d z g ( x , z ) g ( y , z ) = f 2 ( x , y ) .

5. The Collision Term

Using this ansatz we now proceed to write down the kinetic equation for f 2 . All through, we shall retain the usual assumptions of kinetic theory, leading us to neglect triple collisions: the streaming term for the two-particle distribution characterizing particles “1” and “2” will thus be altered by (1) binary collisions between “1”and another particle, “2” being spectator, and (2) binary collisions between “2”and another particle, “1” being spectator.
The binary interaction is defined as occurring when two particles meet in a ball B of radius R. Defining ternary interactions is more subtle since, inasmuch as the interaction potential is the same whatever the order of the interaction, it seems artificial to introduce a specific cutoff. We shall therefore define the range of triple collisions as the lenticular overlap of balls B R ( 1 ) and B R ( 2 ) characterizing the domain of interaction with “1” and “2” respectively. Neglecting triple collisions thus amounts to assuming that | q 1 q 2 | > 2 R .
We first compute the contribution of collisions of “2” with “3”, “1” being left aside. Let us recall that the collision term is given by
f 2 t c o l l = d x 3 V 13 q 1 f 3 p 1 + V 23 q 2 f 3 p 2 .
In the usual derivation of the Boltzmann equation from the BBGKY hierarchy, the right-hand side of BBGKY1 is transformed using BBGKY2. Similarly, we can transform ( t f 2 ) c o l l using BBGKY3, namely
f 3 t + p 1 m f 3 q 1 + p 2 m f 3 q 2 + p 3 m f 3 q 3 V 12 q 1 p 1 p 2 f 3 V 13 q 1 p 1 p 3 f 3 V 23 q 2 p 2 p 3 f 3 = f 3 t c o l l
(we do not make explicit the collision term ( t f 3 ) c o l l since we shall cancel it soon anyway). Under usual dimensional assumptions, we can write t f 3 0 and ( t f 3 ) c o l l 0 , so that, substituting in the collision term, t f 2 c o l l is rewritten as
f 2 t c o l l = d x 3 p 1 m f 3 q 1 + p 2 m f 3 q 2 + p 3 m f 3 q 3 V 12 q 1 p 1 p 2 f 3 + V 13 q 1 + V 23 q 2 f 3 p 3 = d x 3 p 2 m f 3 q 2 + p 3 m f 3 q 3
(the last term vanishes due to the boundary condition f 3 ( | p | ) = 0 , the penultimate since “1” and “2” are supposed far apart from each other and the first because f 3 depends but weakly on q 1 ). More precisely,
f 2 t c o l l = q 3 B R ( 2 ) d q 3 d p 3 p 2 m f 3 q 2 + p 3 m f 3 q 3 .
The following is standard [15]. We introduce the relative coordinate r 23 = q 3 q 2 and use Gauss’ theorem in order to rewrite t f 2 c o l l as a surface integral, so that
f 2 t c o l l = r 23 B R d r 23 d p 3 p 3 p 2 m r 23 f 3 ( q 1 , p 1 , q 2 , p 2 , q 3 , p 3 , t ) = S R d p 3 d Σ · p 3 p 2 m f 3 ( q 1 , p 1 , q 2 , p 2 , q 3 , p 3 , t ) = S R S R + d p 3 d Σ · p 3 p 2 m f 3 ( q 1 , p 1 , q 2 , p 2 , q 3 , p 3 , t ) ,
where d Σ denotes the surface element of the sphere S R such that | r 23 | = R . The southern hemisphere is interpreted as the contribution of oncoming collisions since ( p 3 p 2 ) · d Σ < 0 , while the northern one is the contribution of ending collisions since ( p 3 p 2 ) · d Σ > 0 .
Orienting the polar axis along p 3 p 2 , we have d Σ · ( p 3 p 2 ) = | p 3 p 2 | R 2 sin θ cos θ d θ d ϕ . This can be re-expressed in terms of the surface element of the azimuthal plane such that θ = π / 2 . Letting r denote the radial component on the plane, we obviously have r = R sin θ , whence d r = ± R cos θ d θ (depending on θ being lesser or larger than π / 2 ) and d Σ · ( p 3 p 2 ) = ± | p 3 p 2 | d ω . The collision term can thus be rewritten as (approximating q 3 q 2 )
f 2 t c o l l = a f t e r d p 3 d ω | p 3 p 2 | m f 3 ( q 1 , p 1 , q 2 , p 2 , q 2 , p 3 , t ) b e f o r e d p 3 d ω | p 3 p 2 | m f 3 ( q 1 , p 1 , q 2 , p 2 , q 2 , p 3 , t ) .
Now comes the ansatz. Before the collision we obviously have
f 3 ( q 1 , p 1 , q 2 , p 2 , q 2 , p 3 , t ) = g ( q 1 , p 1 , q 2 , p 2 , t ) g ( q 1 , p 1 , q 2 , p 3 , t ) g ( q 2 , p 2 , q 2 , p 3 , t ) .
The ansatz may be extended after the collision using the fact that, by Liouville’s equation, f 3 ( q 1 , p 1 , q 2 , p 2 , q 2 , p 3 , t ) = f 3 ( q 1 τ , p 1 , q 2 τ , p 2 , q 2 τ , p 3 , t τ ) , where τ is the retardation such that at t τ the particles are entering the collision range with momenta p 2 , p 3 . Since q i τ q i and t t τ , and since p 2 , p 3 are pre-collisional momenta, the ansatz may be introduced in the first integral as well with arguments p 2 , p 3 . We are therefore eventually led to the following Boltzmann-like form for the BBGKY2 equation:
f 2 t + p 1 m f 2 q 1 + p 2 m f 2 q 2 V 12 q 1 p 1 p 2 f 2 = d p 3 d ω | p 3 p 2 | m ( g ( q 1 , p 1 , q 2 , p 2 ) g ( q 1 , p 1 , q 2 , p 3 ) g ( q 2 , p 2 , q 2 , p 3 ) g ( q 1 , p 1 , q 2 , p 2 ) g ( q 1 , p 1 , q 2 , p 3 ) g ( q 2 , p 2 , q 2 , p 3 ) ) + ( 1 2 ) .
The last term accounts for the contribution of collisions undergone by particle “1”.

6. Final Remarks

As promised, equation (21) is coherent for f 2 since g can—in principle—be solved in terms of f 2 . Unfortunately, we are not aware of any readily available solution to the integral equation (13). It is interesting to note that kinetic theory provides a motivation for studying the mathematical object (12), which proves surprisingly involved in spite of its deceptive apparent simplicity. In the authors’ opinion, it might turn sound to tackle the problem from a linearized vantage point, considering situations where the particles are almost uncorrelated, that is where g is almost factorizable. Such an approach would also be in line with usual methods of kinetic theory [16].
This current implicit form of the collision term bears a close ressemblance with the one appearing in the standard Boltzmann equation. This ressemblance might however turn deceptive since g is likely to be a complicated functional of f 2 . We hope that the present letter can prompt further work on this collision term, which might eventually lead to a form susceptible of a hydrodynamical treatment. Our hope is that such a treatment could lead to a deeper understanding of fluctuating phenomena, for which the first-order theory provides only lacunary insights. In particular, it seems reasonable to expect that (21) can be cast in the lattice-based formalism which proved so successful for the usual Boltzmann equation.

Acknowledgments

The authors would like to acknowledge funding from the European Union Seventh Framework Programme, under grant agreement 317534 (Sophocles).

Author Contributions

G.C., O.M. and B.C. performed the research; B.C. supervised the project; G.C. wrote the manuscript. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jaynes, E.T. Information Theory and Statistical Mechanics. Phys. Rev. 1957, 106. [Google Scholar] [CrossRef]
  2. Jaynes, E.T. Information Theory and Statistical Mechanics. II. Phys. Rev. 1957, 108. [Google Scholar] [CrossRef]
  3. Schneidman, E.; Still, S.; Berry, M.J.; Bialek, W. Network Information and Connected Correlations. Phys. Rev. Lett. 2003, 91, 238701. [Google Scholar] [CrossRef] [PubMed]
  4. Schneidman, E.; Berry, M.J.; Segev, R.; Bialek, W. Weak Pairwise Correlations Imply Strongly Correlated Network States in a Neural Population. Nature 2006, 440, 1007–1012. [Google Scholar] [CrossRef] [PubMed]
  5. Stephens, G.J.; Bialek, W. Statistical Mechanics of Letters in Words. Phys. Rev. E 2010, 81, 066119. [Google Scholar] [CrossRef]
  6. Mora, T.; Bialek, W. Are Biological Systems Poised at Criticality? J. Stat. Phys. 2011, 144, 268–302. [Google Scholar] [CrossRef]
  7. Bialek, W.; Cavagna, A.; Giardina, I.; Mora, T.; Silvestri, E.; Viale, M.; Walczak, A.M. Statistical Mechanics for Natural Flocks of Birds. Proc. Natl. Acad. Sci. USA 2012, 109, 4786–4791. [Google Scholar] [CrossRef] [PubMed]
  8. Stephens, G.J.; Mora, T.; Tkačik, G.; Bialek, W. Statistical Thermodynamics of Natural Images. Phys. Rev. Lett. 2013, 110, 018701. [Google Scholar] [CrossRef] [PubMed]
  9. Van der Straeten, E. Maximum Entropy Estimation of Transition Probabilities of Reversible Markov Chains. Entropy 2009, 11, 867–887. [Google Scholar] [CrossRef]
  10. Marre, O.; El Boustani, S.; Frégnac, Y.; Destexhe, A. Prediction of Spatiotemporal Patterns of Neural Activity from Pairwise Correlations. Phys. Rev. Lett. 2009, 102, 138101. [Google Scholar] [CrossRef] [PubMed]
  11. Cavagna, A.; Giardina, I.; Ginelli, F.; Mora, T.; Piovani, D.; Tavarone, R.; Walczak, A.M. Dynamical Maximum Entropy Approach to Flocking. Phys. Rev. E 2014, 89, 042707. [Google Scholar] [CrossRef]
  12. Chliamovitch, G.; Dupuis, A.; Golub, A.; Chopard, B. Improving Predictability of Time Series Using Maximum Entropy Methods. Europhys. Lett. 2015, 110. [Google Scholar] [CrossRef]
  13. Shannon, C.E. A Mathematical Theory of Communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  14. Khinchin, A.Y. Mathematical Foundations of Information Theory; Dover: Mineola, NY, USA, 1957. [Google Scholar]
  15. Kreuzer, H.J. Nonequilibrium Thermodynamics and its Statistical Foundations; Oxford University Press: Oxford, UK, 1984. [Google Scholar]
  16. Liboff, R.L. Kinetic Theory; Springer: New York, NY, USA, 2003. [Google Scholar]

Share and Cite

MDPI and ACS Style

Chliamovitch, G.; Malaspinas, O.; Chopard, B. A Truncation Scheme for the BBGKY2 Equation. Entropy 2015, 17, 7522-7529. https://0-doi-org.brum.beds.ac.uk/10.3390/e17117522

AMA Style

Chliamovitch G, Malaspinas O, Chopard B. A Truncation Scheme for the BBGKY2 Equation. Entropy. 2015; 17(11):7522-7529. https://0-doi-org.brum.beds.ac.uk/10.3390/e17117522

Chicago/Turabian Style

Chliamovitch, Gregor, Orestis Malaspinas, and Bastien Chopard. 2015. "A Truncation Scheme for the BBGKY2 Equation" Entropy 17, no. 11: 7522-7529. https://0-doi-org.brum.beds.ac.uk/10.3390/e17117522

Article Metrics

Back to TopTop