An Alternative Approach to Particle-Particle Collisions

A new alternative approach to the statistical behavior of particle-particle collisions is introduced. The alternative approach is derived rigorously from well known mechanical laws; and the results given by it, quantitatively and qualitatively different from what the standard kinetic theory yields, can be directly checked with computer-simulated or realistic experiments. More importantly, from the introduction of it, a number of new concepts and new methodologies emerge, which might turn out to be very significant to the future development of nonequilibrium statistical mechanics.


Introduction
It is commonly believed that the statistical behavior of classical particle-particle collisions has been adequately codified into the equations of nonequilibrium gas dynamics (kinetic theory), and there is nothing truly new to today's physicists.However, as has been pointed out by us [1] [2], this belief is questionable.
What made us question the validity of the standard theory might be summarized as the following.In the context of gas dynamics, particle-particle collisions are supposed to be examined in the six-dimensional positionvelocity phase space, but the space, though superficially simple, is inherently counterintuitive.When a sixdimensional statistical dynamics is of concern, it is often the case that our attention is invited to certain incomplete and misleading pictures.To get a taste of such trickiness, let's briefly review a long-neglected fact associated with the Boltzmann equation [3]: while the differential operator on its left-hand side enjoys a relatively obvious position-velocity symmetry, the integral operator on its right-hand side does not exhibit any form of position-velocity symmetry.In this regard, one may reasonably wonder: How can the equation's two sides be always equal when they behave themselves so differently?
In this paper, after advancing several typical examples, an alternative approach is proposed, which takes care of the position space and the velocity space in a relatively balanced way.The proposed formalism is strictly based on well-known mechanical laws, and the results given by it can be directly compared with computer simulations.More importantly, from the introduction of this alternative approach, a number of new concepts and new methodologies emerge, which might turn out to be very significant to the future development of nonequilibrium statistical mechanics.
The structure of this paper is the following.In Section 2, several examples are advanced in which particleparticle collisions cannot be well treated by the standard kinetic theory.In Section 3, an alternative approach is proposed.In Section 4, the newly proposed alternative approach is extended to more general situations.In Section 5, a brief summary is provided.

Typical Examples of Particle-Particle Collisions
To investigate the collective behavior of particle-particle collisions, we shall in this paper concern ourselves with three different, but interconnected, situations shown in Figure 1.They are a) a parallel beam of particles scattered by a group of resting particles, b) two parallel beams of particles colliding in the head-on manner, and c) two ordinary gases colliding with each other.
To begin with, let's first look at the simplest situation illustrated in Figure 1(a).Suppose that the parallel beam on the left side is moving in the z -direction (rightwards), and the distribution function of it is ( ) ( ) ( ) ( ) where ( ) is considered to be time-independent and z -independent and The resting particles on the right side are described by the distribution function ) , , , 0 0 0 where ( ) There are several reasons why we start our investigation with this particular example.Firstly, this is the case to which the standard collision theory is presumably applicable.Secondly, this is the case that we, as physicists, can easily realize and easily monitor in computer simulation.Thirdly, this is the case in which the complexity of the subject is largely reduced but the statistical characteristics of the subject remain intact (so the insight gained will be generally instructive).
It is further assumed that a virtual detector, marked as 1 D in the figure, is placed at the position where the local distribution function is of interest.This arrangement possesses great importance in this paper; and the basic idea behind it is that the distribution function determines, almost uniquely, the local flux, and hence the distribution function should be directly computable if the flux information gathered by the detector situated there is somehow known.It will be seen that with certain improvement and refinement this simple idea works nicely.In contrast with that, the basic idea of the standard theory seems rather cumbersome.According to the standard theory, to determine the distribution function at a place we are supposed to investigate how the distribution function varies at the place, and to determine how the distribution function varies at the place we are supposed to investigate virtually all the fluxes around the place (eventually 12 coming-or-going fluxes are examined).
Before entering the next section, where a new alternative approach will be introduced, let's briefly go through what the Boltzmann equation has to say about the situation.Referring to Figure 1(a), suppose that the symmetry axis of 1 D is the positive ξ -coordinate axis and the origin of the ξ -coordinate starts from the bottom of the collision region.Thus, the distribution function of the collision-produced particles moving along the ξ -axis obeys (with no external forces) ( ) in which ′ v and 1 ′ v are the initial velocities of two colliding particles while v ξ and 1 v stand respectively for the final velocities of the two particles, σ is the cross section of the scattering and d c Ω is the solid angle of the scattering.(Notably, c σ and c Ω are defined in the so-called center-of-mass frame.)From Equation ( 3), ( ) for a definite v ξ can be determined by ( ) ( ) In writing these two formulas, it is understood that the two initial gases are dilute enough so that the collision probability between f ′ and 1 f ′ is a first-order quantity and the collision probability between f and f ′ (or 1 f ′ ) can be ignored.It is also understood that the relationship between v ξ , 1 v , ′ v and 1 ′ v in Equation (3), or (4), obeys the energy-momentum conservation law, which in general consists of four independent equations.
Surprisingly, although formally pertinent, expression (4) yields no meaningful result for the situation.Firstly, if we wish to compute the integral involved, we find no attainable way to get rid of the five δ -functions in f ′ and 1 f ′ .Secondly, if we assume that the five δ -functions disappear simply due to the five-fold velocity inte- gration of Equation ( 4), the role played by the four equations of the energy-momentum conservation law becomes unknown.Finally, if we somehow managed to make Equation (4) mathematically evaluable, the issue would become even more puzzling in the sense that the result of Equation ( 4) can be schematically represented by the dotted curve in Figure 2 (reflecting the fact that the integrand of the integral is positive-definite), while our physical and geometrical intuition, as well as a well-performed computer simulation, tells us that the real distribution function varies according to the solid curve given in Figure 2 (the asymptotic behavior of f must obey the inverse-square law).There has been a detailed investigation about the sources of the aforementioned problems [1].In this paper, we shall focus ourselves on proposing a new alternative approach.

An Alternative Approach
In connection with the situation shown in Figure 1   the aftermath of a collision is illustrated in the position space and in the velocity space respectively.The two figures are drawn quite generally and their full implications will manifest themselves in this and the next sections.
In Figure 3(a), suppose that ( ) is the position of the detector's inlet, and that the velocities of the particles recorded by the detector fall into a definite speed range v δ and a definite solid angle δΩ .Then, the distribution function at the inlet can be roughly expressed as where s δ is the area of the detector's inlet, t δ is the time interval of recording, N δ is the number of the particles recorded by the detector.Expression (5) can be interpreted as saying that there is a direct connection between the flux detected by the detector and the distribution function of our primary concern.It is this direct connection, though still primitive at this stage, that will guide us to the very end of our investigation.
Concerning expression (5), there are essential things worth discussing.To ensure that the right side of the expression stands for the "exact" distribution function at the position point r , the quantities t δ , s δ , v δ and δΩ have to be infinitely small. (By the term exact, we mean that the right side of the expression can be quali- fied as a strict mathematical limit.)However, allowing all these quantities to be simultaneously infinitesimal gives rise to practical problems.For instance, if both s δ and δΩ are indeed infinitely small, dN in expres- sion (5) receives contribution only from the collision region that is limited to a one-dimensional straight line, which is nothing but the very symmetry axis of the detector.Obviously, this type of approach can be accomplished only in mind.
To make expression (5) physically meaningful in rather general situations, we shall adopt the following two assumptions: 1) s δ is very small, infinitesimal in the theoretical sense, and 2) v δ and δΩ are finite and de- finite (still small though).It turns out, to describe force-free Boltzmann gases, accepting these two assumptions is convenient and largely necessary.
Under these two assumptions, we redefine the distribution function as a mathematical hybrid: ( ) where ( ) being respectively the lower and upper bounds of v δ , and dN includes all the particles recorded during t δ whose velocities fall in v δ and δΩ .Namely, this distri- bution function is a mathematical limit in the position space, but a mathematical average in the velocity space, averaged over the finite volume element ( ) At this point, one remark seems in order.In the standard theory the distribution function, as the primary concept of the theory, is often defined or interpreted in intentionally vague language.For instance, in one of the books on nonequilibrium statistical mechanics [4] it is stated that the distribution function, denoted as F therein, is the ratio δ and v δ are supposedly "large enough to contain a lot of par- ticles, but small compared to the range of variation of F ′′ .Then, an interesting question about expression (6) arises.Why does the approach herein define the distribution function in such clear and strict mathematical terms?To this question, the answer given by us is simply the following: to compute a distribution function certain limiting processes need to be invoked; but if there is not an adequate set of mathematically well-defined quantities those limiting processes may lose their validity explicitly or implicitly.
We take three steps to compute expression (6): 1) finding out the position region in which particle-particle collisions may possibly give contribution to dN ; 2) formulating the collision number in this position region; and 3) investigating how collision-produced particles spread in the position space and in the velocity space, and then determining in what portion those collision-produced particles will indeed be recorded by the detector.
The first two steps can be accomplished rather easily.Since the solid angle δΩ has been assumed to be de- finite and finite and the inlet area s δ has been assumed to be very small (infinitesimal), we find the trajectory cone denoted by { } r δ − Ω in Figure 3(a) to be the position region in which collisions may directly give contri- butions to dN .In this particular sense, the cone will be named as the upstream path zone, or up-zone for short.The number of the collisions taking place inside the up-zone can be expressed by ( ) ( ) where 0 dr stands for an infinitesimal region in { } r δ − Ω .Inserting Equations ( 1) and ( 2) into Equation (7), we obtain To accomplish the third step aforementioned, let's first note that the information concerning how the scattered particles spread in the position space and the velocity space is stored in c c σ δΩ .Since c c σ δΩ is conveniently defined in the center-of-mass frame, and the position and velocity spaces we care about are defined in the laboratory frame, an extra examination is needed.
Referring to Figure 3(a) and Figure 4(a), we find that there is an infinitely sharp cone formed by the point 0 r , as the apex, and the area element s δ , as the base.So, the solid angle bounded by the cone, δω , is equal to Since s δ is infinitely small, δω is infinitely small to δΩ .Thus, if a particle can enter the detector from the up-zone, the particle's velocity will definitely falls in δΩ .Namely, from now on, the only thing we need to make sure is that the collision-scattered particles emerge within v δ and δω in the laboratory frame.To connect the center-of-mass frame and the laboratory frame, we have to deal with the energy-momentum C. Y. Chen conservation law explicitly.Let all the particles in our consideration have the same mass (for simplicity); and let ′ v and 1 ′ v be the initial velocities of two colliding particles and v and 1 v be the final velocities of the two particles respectively.The energy-momentum conservation law then takes the form and u where The above expressions show that for a pair of ′ v and 1 ′ v , the final velocities v and 1 v have only two degrees of freedom, distributing on an infinitely thin energy-momentum (EM) membrane in the velocity space, as has already been illustrated in Figure 3 ) in which v = v and 0 2 in this section).This relationship defines the patch of the EM membrane on that the scattered particles will indeed be recorded by the detector.Finally, by virtue of Figure 4(b), we have ( ) With help of Equations ( 8), ( 9), ( 12) and ( 13), the limit-average-hybrid distribution function defined by Equation ( 6) is .In this formula only the particles coming from f ′ are taken into account.Unlike expression (4), expression ( 14) is obviously computable.It yields a finite and definite result no matter whether r is inside or outside the collision region.When r is located at a place remote to the collision region (which defines the asymptotic behavior of the particle-particle collisions), f in Equation ( 14) obeys the in- verse-square relationship, thus ensuring the conservation of particle number globally.

Extensions of the Alternative Approach
Now, we look at how to extend our proposed approach to the cases shown in Figure 1(b) and Figure 1(c).
Let the gas on the left-hand side of Figure 1(b) be represented by where n′ is time-and z -independent and ) ( ) where 1 n′ is also time-and z -independent and ( ) By defining the up-zone, { } r δ − Ω , around the symmetry axis of the detector, we know that the collision number taking place in the zone can still be expressed by Equation (7).So, the collision-produced particles that can enter the detector are

[
] [ ] in which d c Ω is limited to the patch on that the scattered particles will move almost exactly toward the detec- tor in the laboratory frame.Now, the energy-momentum conservation law shows us that and and Equation (17) becomes is the Jacobian of the variable transformation.Notice that Figure 4(a) makes sense generally as long as we disregard = c u in it.(Obviously, = c u holds only for the case discussed in the last It is then ob- vious that Equation (9) still holds.Furthermore, for the fixed ′ = c c in that figure, we have where 2 d d v v ω is the infinitesimal volume element into which the particles recorded by the detector are sup- posed to fall.Eventually, with help of Equations ( 9), ( 19) and ( 21) the distribution function defined by Equation (6) becomes ( ) ( ) in which the integration region of 0 dr covers the up-zone { } r δ − Ω , the integration range of d z c′ is from −∞ to +∞ and c σ is the function of ( ) , θ θ ′ = u u .Other quantities should be determined with care, and the fol- lowing diagram may provide an additional help ( ) g v′ can be determined.As for the general case shown in Figure 1(c), the formulation is about the same.Firstly, the collision number taking place in the up-zone can be represented by ( ) ( ) where f τ ′ and 1 f τ ′ are the distribution function at the collision position 0 r and at the collision time τ .By invoking the variable transformation defined by Equation (11), we arrive at Notice that in Equation (25) d ′ u can be replaced by 2 d d u u ′ Ω where ′ Ω stands for the direction of ′ u .Also notice that Equations ( 9) and (21) hold significance in this general situation.Eventually, we find that the distribution function defined by Equation (6), at a particular time t , is ( ) ( ) which is interpreted similarly to the diagram given in expression (23).
Obviously, both expressions ( 22) and ( 27) are directly computable, and ready to be checked with computer-simulated realistic experiments.
It is easy to see that if Ω δ in the formulas of this paper is allowed to approach zero, the whole formalism makes no sense, which again says that the "exact" distribution function is inherently indeterminable.More than that, this formalism has actually shown that, if we are a weightless and infinitely small observer moving together with a group of collision-produced particles, the distribution function observed by us is, in general, not pathinvariant, which in turn challenges the long-held tacit assumption that any distribution function, produced by whatever mechanisms, is a continuous one or can be treated as a continuous one.

Summary
In this paper, to formulate the statistical behavior of particle-particle collisions, a new alternative integral formalism has been introduced.The results given by the new formalism are quantitatively and qualitatively different from what the standard theory yields.If interested, readers may confirm or deny them by applying their own theoretical and/or numerical approaches.
More importantly, along with the introduction of the new approach, a set of new concepts and methodologies are proposed, which might turn out to be very significant to the future development of nonequilibrium statistical mechanics.Some of them are:  Due to the discontinuity concern, the distribution function in a nonequilibrium approach should be defined as an average in the position-velocity phase space, at least partially. Instead of using differential-integral equations, approaches of completely integral type should be employed. Instead of examining the events in a control volume element, what takes place in an upstream path zone should be investigated. Collisional effects should be studied in both the position space and the velocity space.The energymomentum conservation law should be fully incorporated.

Figure 1 .
Figure 1.Various setups: (a) A parallel beam of particles scattered by resting particles; (b) Two parallel beams colliding in the head-on manner; (c) Two ordinary gases colliding with each other.

Figure 1 (
b) and Figure1(c), similar difficulties will be met with.In particular, the difficulty shown in Figure2is always there.
(a), we now look at Figure 3(a) and Figure 3(b), in which

Figure 2 .
Figure 2. The qualitative illustration of the distribution function along the axis of the detector in Figure 1(a).

Figure 3 .
Figure 3. Aftermath of a collision: (a) In the position space one of the two colliding particles starts to move toward the detector; (b) in the velocity space the two particles get new velocities, v and 1 v respectively.

Figure 4 .
Figure 4. Schematic of (a) how the scattered particles with + = c u v enter the detector, and (b) how z v δ ′ in this section can be determined by v δ in view of 1 0 z v′ = .
(b).Referring to Figure 3(b) and Figure 4(a) simultaneously, we find out that in the velocity space, for the fixed ′ = c c , (

in which the integration region of 0 dr
is over the up-zone { } r δ − Ω , the integration range of d

,
+∞ is equal to 1; let the gas on the right- hand side of Figure 1(b) be represented by

r
defined by Equation (6), where r stands for the position of the detector's inlet in Figure 1(b).
v we obtain v , by geometrically subtracting c from v we get u , with help of u = u the function ( )

∫
r is over the upstream path zone, dv ∫ is over v .Other quantities in this expres- sion can be determined by