Studying the properties of plasma has received a great interest in both astrophysical and laboratory plasma application. The relativistic binary distribution function is one of the most important functions of statistical mechanics. The importance of the distribution function in statistical mechanics is due to the fact that all the thermodynamic quantities, such as the pressure, the internal energy and the free energies, can be calculated from it.

Relativistic statistical mechanics has a long story, but we may notice that, whereas the theory of relativistic ideal gases has received deep and detailed developments, little has been achieved in order to account for mutual interactions between particles [1] . In statistical physics, the BBGKY hierarchy (Bogoliubov-Born-Green- Kirkwood-Yvon hierarchy, sometimes called Bogoliubov hierarchy) is a set of equations describing the dynamics of a system of a large number of interacting particles. The equation for an s-particle distribution function (probability density function) in the BBGKY hierarchy includes the -particle distribution function thus forming a coupled chain of equations [2] . Their history of 70 years has brought enormous progress in the investigation of the transition from the microscopic to the macroscopic world, and they are still an attractive starting point for new developments. In particular, great advance has come by clever methods of truncation, approxi- mation and scaling limits of the hierarchy, providing in various cases a justification of the kinetic equations describing particle systems on mesoscopic (intermediate) scales [3] . Many authors studied the BBGKY hierarchy [4] and [5] ; Hussein and Eisa [6] [7] calculated the binary and triplet distribution function for one- and two- component plasmas for quantum and classical non-relativistic plasma. In this paper, we will calculate these distributions in the relativistic case in terms of plasma thermal parameter. Arendt and Eilek [8] showed that the

pair-plasma distribution functions could be described by a thermal parameter that was moderately

relativistic. The thermal parameter value plays a significant role for the stability of our system. Barcons and Lapiedra (1984) [9] gave explicit expressions for the thermodynamic functions of a high-temperature electron- positron plasma and gave expression for distribution functions for a classical dilute arbitrarily hot plasma in equilibrium which we compared our results with it.

Special relativity, however, does not permit velocities greater than the speed of light and is thus incompatible with a Maxwellian distribution that predicts a non-zero probability for every velocity [10] . It is possible to obtain the relativistic Maxwellian distribution (Jüttner distribution) of a moving gas when the so-called Planck-Einstein case is analyzed. It should be reiterated that the Planck-Einstein theory suffered a modification, principally in the transformation law of energy (the Planck-Einstein case) [11] . Treumann et al. [12] have also shown that the isotropic thermal equilibrium distribution function in relativistic plasma analytically has the form of a modified Jüttner distribution.

A little is known about relativistic distribution functions involving more than two particles, and in particular about the three-particle (or triplet) distribution function. This is of course due to the greater mathematical complexity of higher order correlation functions, and to a lack of a direct link with experiment. Although one can consider experimental determination of triplet distribution function from triple elastic scattering similar to that of the binary distribution function, but to our knowledge, the measurement of three-body correlation function requires the knowledge of the positions of three particles at the same time which is technically very demanding to obtain in 3D samples [13] ; this means that the experiment to measure the triplet distribution function directly requires high precision so that the multiple scattering can be differentiated from single scattering. Essentially, one follows the phase-space trajectory of the system as it evolves in time, and has thus the same amount of information as one obtains in a simulation. Similarly the particle positions have measured in dusty plasmas [14] . Lapiedra et al. [15] have undertaken an application of predictive methods to relativistic statistical mechanics. Coulomb forces between point charges are purely repulsive and charges approach very close to each other. Coulomb systems, such as plasma or electrolytes, are made of charged particles interacting through Coulomb’s law. The simplest model of a Coulomb system is the one-component plasma (OCP), also called jellium: an assembly of identical point charges, embedded in a neutralizing uniform background of the opposite sign. Here we consider the classical (i.e. non-quantum) equilibrium statistical mechanics of the OCP. It is rather straightforward to calculate higher-order correlation functions from the measured configurations. Moreover we study the model of two-component plasma (TCP) i.e. neutral system of point like particles of positive and negative charges such as electrons and ions. For the numerical calculation we restrict ourselves to the case of (TCP) which anti-symmetric with respect to the charges and therefore symmetrical with respect to the densities. To simplify the numeric investigations, we simulated so far only mass symmetrical electron-ion plasma with.

Homogeneous plasma is characterized by two parameters: the density of particles and the temperature; if there are several kinds of particles, it is also necessary to state their concentrations. There are different energies associated with the plasma, namely:

1) Energy of the rest mass per particle;

2) Kinetic energy per particle of order;

3) Coulomb energy of order per particle.

The ratios of these energies give us the two main dimensionless parameters of the plasma: the thermal para-

meter of plasma and the dilution parameter. The value of is very important. There

are four different regimes characterized by it: plasma with is relativistic, weakly relativistic plasma, low plasma temperature and ultra-relativistic high temperatures plasma. In our study

we note that the system is safely classical if where, , is the Planck constant, and

is the typical linear momentum of particles. A relativistic plasma with a thermal distribution function has temperatures greater than around 260 keV, or 3.0 GK (5.5 billion degrees Fahrenheit), where approximately 10% of

the electrons have where is the Lorentz factor. It is occurs in many environments in

astrophysics, including gamma-ray bursts, AGN jets, and pulsar winds. In physics, a particle is called ultrarelativistic when its speed is very close to the speed of light. Or, similarly, in the limit where the Lorentz factor is very large and. Our study is only valid for dilute plasmas. The plasma temperatures and velocities also plays a significant role for the stability of our system Figure 1 show qualitative sketch of the four different regimes characterized by the thermal parameter and Lorentz factor. Our model lies in the region of high speeds, and.

In the astrophysical environment, the fraction of ionized particles varies widely from nearly no ionization in cold regions to fully ionized in regions of high temperature. This leads to a wide range of parameters where astrophysical plasmas can exist. While the astrophysical environment is frequently dominated by the presence of the plasma, this plasma is often strongly influenced by and coupled to the presence of embedded particulates (i.e., dust). These dust grains which range in size from a few nanometers to micron-sized objects can become either positively or negatively charged due to interactions with the background plasma environment and ionizing radiation sources in the astrophysical environment [16] . Understanding the processes that govern these plasma particle interactions is critical to the study of astrophysics. The agglomeration and growth of larger particles from single atoms and dust grains leads to the eventual formation of objects so large that gravity becomes the dominant force controlling their subsequent evolution [17] .

Another important point is that of the walls. We shall study plasma which is homogeneous and isotropic but the plasma must be confined or it will expand. We may assume that the plasma is confined by some kind of walls which prevent the escape of particles, but that the container is so large that the effects of the walls are negligible. Our calculations are based on the phase-space distribution function; this is defined as the number of particles per unit volume of space per unit volume of velocity space: At time t, number of particles in elementary volume of space, with velocities in range.

The statistical state of a macroscopic system of N particles is in a complete―though in an intractably complex―

way described by the distribution function in dimensional phase space, which is spanned by the coordinates and velocities of all individual particles, , being the three posi-

tion of particle, and the spatial components of its four velocity in a given frame.

The s-particle reduced distribution is giving by:

The relativistic BBGKY hierarchy [15] is given by

where is the acceleration of the charge in the presence of the charge.

The Einstein summation convention is only valid for Greek labels. Therefore, for each value of s we have s equations, since. As in the non-relativistic case, the determination of the reduced generalized distribution function can only be made when the hierarchy is cut off somewhere, that is, when for some value s we give as a function of the other functions with.

According to what we done in non-relativistic case [7] we set, whatever particles 1, 2, 3 are,

Let us consider the case of homogeneous plasma in equilibrium. Then is the one-particle distribution function of an ideal gas see Figure 2 which shows one particle relativistic distribution function in the particle velocity interval (0, 0.9c) for different values of the thermal parameter μ. That is, in a frame relative to which the system is macroscopically at rest, we must set for the relativistic Maxwellian distribution [18]

where denotes the modified Bessel function.

At first sight, the calculation of the relativistic interaction between two charged particles seems rather involved, because the force on particle 1 at time would depend on the position and velocity of particle 2 at a retarded time. Also, the position and velocity of 1 at that time depends on the position and velocity of 1 at an earlier time, and so on. If one chooses a frame of reference, the acceleration on particle 1 due to the presence of 2 can be calculated from the position and velocity of 2 at that time. This does not mean that actions propagate instantaneously, but rather that there is a precise scheme to take retardation into account automatically, through the equations of the theory. All we shall need in that chapter is the time component of the acceleration of particle 2 due to particle 1 [19] .

To calculate the binary relativistic distribution function substituting from Equations (3) and (5) into (2) for we obtain

Then, by solving the integro-differential Equation (7) and by using the Fourier transform of

i.e.,. We get the relativistic binary

generalized distribution function at in the following form

Let us now study the two-particle relativistic distribution for the model of two-component plasma (TCP) i.e. neutral system of point like particles of positive and negative charges such as electrons and ions. For the numerical calculation we restrict ourselves to the case of two-component plasma which anti-symmetric with respect to the charges and therefore symmetrical with respect to the densities.

For two-component plasma we can use the two-particle correlation function which is given by

where is the Debye-Hiickel solution [20] for two-component plasma, by substituting Equation (9) into Equation (3) we get the binary distribution function for two-component plasma in the following form:

with given by Equation (5) (in fact only terms of up to must be retained). It can be verified that given by Equation (9), with satisfies identically the condition (6) for, to first order in.

The standard two-body distribution function for dilute slightly relativistic plasma has been calculated previously by Kosachev and Trubnikov [21] starting from the Darwin Lagrangian. Lapiedra and Santos [19] result agrees with theirs to order, but not to higher-order terms. We think that these terms are meaningless unless one goes beyond the Darwin Lagrangian, which is only correct to order. The Lagrangian approach introduces some modification of the expression for the statistical sum and what is more important the use of the restricted Breit-Darwin Hamiltonian leads to the wrong behavior of the pair correlation function and therefore to the incorrect expressions for various thermodynamical quantities [22] .

According to our knowledge few scientists studied the binary distribution function for relativistic dilute plasma [15] ; the new in our article is using effects of thermal parameter on the values of binary relativistic distribution function for one and two-component plasma.

The triplet distribution function is defined in such a way that, the calculation of the relativistic interaction between three charged particles seems rather involved, because the force on particle 1 at time t would depend on the position and velocity of particles 2 and 3 at a retarded time. Also, the position and velocity of 1 at that time depends on the position and velocity of 1 at an earlier time.

Substituting Equation (4) and (5) into (2) for we obtain

If we used the Kirkwood superposition approximation (KSA) [23] ; which is consisting of the assumption that the potential in a set of three particles is the sum of the three pair potentials, this is equivalent to assuming that the triplet distribution function is the product of the three radial distribution functions

It has been mentioned that before 2003 there is no direct measurement of three-body correlation function and this is the first study for the triplet distribution function in the case of dilute relativistic plasma. Such measurement requires the knowledge of the positions of three particles at the same time which is technically very demanding to obtain in 3D samples [13] ; this means that the experiment to measure the triplet distribution function directly requires high precision so that the multiple scattering can be differentiated from single scattering. Now because of the existence of video-microscopy, a modern experimental technique applied to colloidal systems to directly measure all particles’ positions at all times. Essentially, one follows the phase-space trajectory of the system as it evolves in time, and has thus the same amount of information as one obtains in a simulation. Similarly the particle positions have measured in dusty plasmas [24] .

From Equation (12) we can obtain the classical TDF for two-component plasma in the following form:

And we also can used (KSA) which is given in Equation (13) to get it in the form

The triplet and quadruple distribution functions as well as binary distribution function must be incorporated for a more accurate and complete discussion of macroscopic equilibrium properties. A little is known about distribution functions involving more than two particles, and in particular about the three-particle (or triplet) distribution functions. This is of course due to the greater mathematical complexity of higher order correlation functions, and to a lack of a direct link with experiment [25] . And this is considered the first study for the triplet distribution functions for dilute relativistic plasma.

In many physical systems, the description of a plasma as a Coulomb system is sufficient to reproduce most of the properties of interest. If the system is cold enough, the mean velocities of the particles are much smaller than the speed of light, and the charges may be assumed to interact via the instantaneous Coulomb potential. However, at sufficiently high temperatures, this approximation is no longer valid, and the contributions of the relativistic effects (which include, apart from the trivial kinetic corrections and of course all the retardation effects) must be incorporated when studying the equilibrium properties of the system [26] .

In the classical (non-quantum) case, the systematic approach adopted here follows the traditional route of the relativistic Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy for the reduced distribution functions by formal density expansion. Interestingly, this study is the first to display the effect of thermal parameter of plasma in the classical binary and triplet distribution function. Also the triplet relativistic distribution function for dilute plasma was calculated from the relativistic BBGKY hierarchy. We used the results to obtain the analytical forms of the classical triplet distribution functions for one- and two-component plasmas.

In Figure 3 and Figure 4 we noticed that the value of both binary and triplet relativistic distribution function increased when the value decreased at very high temperature. Physically these results seem acceptable because from the definition of the phase-space distribution function which means the number of particles per unit volume of space per unit volume of velocity space, and at high temperature, the velocity of the particles increases and as a result the number of particles increases per unit volume.

Our calculations are grounded in the classical relativistic statistical mechanics. Plasma is non-degenerate. The system is not dense, so one may neglect the contributions of higher order particle interactions. Figure 5 and Figure 6 show that the two- and three-particle distribution functions have become concentrated to ever-smaller region of speed, dramatically increasing the thermal parameter. When we made a comparison between the two-particle relativistic distribution function from Equation (11) and the result of Barcons and Lapiedra as shown in Figure 7, the results were nearly similar at very high velocities.

Figure 8 and Figure 9 show that the two- and three-particle distribution functions in the thermal parameter of plasma interval for one-component plasma for different values of speed, , and. We note that the curves are very close to each other, dramatically increasing the thermal parameter

and for three-particle distribution function more than two-particle distribution function.

Figure 10 shows the comparison between the three-particle relativistic distribution function from Equations (14) and (15). One of them is based on the Kirkwood superposition approximation (KSA) which is consisting of the assumption that the potential in a set of three particles is the sum of the three pair potentials. This is equivalent

to assume that the triplet distribution function is the product of the three radial distribution functions, and the other form is calculated by using the relativistic BBGKY hierarchy.

Finally, we can note that from Figures 11-14, the distribution function has become more concentrated to ever- smaller region of speed, dramatically increasing the thermal parameter and for three-particle distribution function more than two-particle distribution function.