The quantum thermodynamic functions of plasma in terms of the Green ’ s function

The objective of this paper is to calculate the third virial coefficient in Hartree approximation, Hartree-Fock approximation and the MontrollWard contribution of plasma by using the Green’s function technique in terms of the interaction parameter ξ , and used the result to calculate the quantum thermodynamic functions for one and two component plasma in the case of ab nλ 3 1  , where ab ab h m KT λ 2 2 = is the thermal De Broglie wave-length. We compared our results with others.


INTRODUCTION
The thermodynamic functions are of great interest for understanding the properties of plasma such as the excess free energy and the pressure.An accurate description of the pressure, volume and temperature (P-V-T) behavior represents one of the most important goals of statistical thermodynamics.Besides simulations, the equilibrium behavior can be determined either from the equation of state (EOS) or from virial expansion.The properties and the behavior of many particle systems with Coulomb interactions are essentially determined by the long range character of the forces between the charged particles.Therefore, systems with Coulomb interactions are of special interest and importance in statistical physics [1].
Many authors have calculated the excess free energy, Ebeling et al. [2] have calculated the virial expansion of the excess free energy until the second virial coefficient from the pair distribution function.Hussein and Eisa [1] calculated the quantum excess free energy for two component plasma.All these researches are used the method of Slater Sum in quantum statistical mechanics.Also, the thermodynamic functions which are calculated from these researches near to the classical limit.In this paper we use the Green's function technique, which means that we have to start from the grand canonical ensemble.To our knowledge, no attempt to find the thermodynamic function until third virial coefficient by using Green's function technique for one and two component plasma.The method of Green's functions is one of the most powerful techniques in quantum statistical physics and quantum chemistry.The main advantage of Green's functions is that one can formulate exact equations of approximations and set up consequent schemes to gain any higher approximation [3].Another advantage of such Green's functions is the existence of highly effective methods for their determination, such as Feynman diagram techniques, functional techniques and the formulation of equations of motion [4].The single particle Green's function in fact contains more detailed information than the total energy alone, to the extent that the local Slater correlation potential can be obtained from it [5].Many authors used the method of Green's functions such as Dewitt et al. [6] and have calculated the low density expansion of the equation of state for a quantum electron gas.Riemann et al. [7] have calculated the equation of state of the weakly degenerate one-component plasma (OCP).
The model under consideration is the two-component plasma (TCP); i.e. neutral system of point like particles of positive and negative charges which antisymmetric with respect to the charges e i e e e = − = − and therefore symmetrical with respect to the densities e i n n n = = .Also, we used the one component plasma model (the model of identical point charges immersed in a uniform background, while the continuous charge density of the background is chosen to be equal and opposite to the average charge density of the point charges, so that the system as a whole is electrically neutral) for example the electron gas i e e e e = = − and therefore symmetrical with respect to the mass e i m m m = = .This paper is organized as follows: In Section 2, we present the Green's function.In Section 3, we calculate the excess free energy until the third virial coefficient for one and two component plasma in quantum form.Also, we calulate the general formula of the third virial coefficient in Hartree-Fock approximation.Finally, in Section 4, we calculate the pressure for one and two component plasma until the third virial coefficient.

THE GREEN'S FUNCTION
The n-particle Green's function is defined by [8] in the following form where, T is the time-ordering operator and So we can defined the greater and lesser Green's function by In the Hartree and the Hartree-Fock approximation the two particle Green's function is given by Or the variational representation [8] ( ) ( ) ( ) ( ) ( ) where the polarization function ( ) ( ) ( ) and ( ) is the effective potential.Also the three particle Green's function is defined by [8]

THE EXCESS FREE ENERGY
The excess free energy corresponds to the part of the free energy change in the real system that arises from interactions among ions [9].The excess free energy ex F can be calculated from the mean value of the potential energy according to a general quantum statistical formula with the coupling parameter λ in the following form [10] ( ) . , 0 1, lim , , ( ) abc v t are the mean interaction potential for two and three particle respectively.
After performing the integration with respect to λ in Equation ( 8) we can get 3 , 12π where the second virial coefficient is given by [4] ( ) , Now we will calculate the third virial coefficient; for Coulomb systems it is useful to apply, instead of v , the screened potential s v in Eq.14 for three particle which is given by the following form , v v v G G G B B B are the binary potential, the triplet potential, the triplet screened potential, the triplet polarization function, one particle Green's function, the two particle Green's function in Hartree approximation, the three particle Green's function, the triplet Hartree term, the triplet Hartree-Fock, the triplet screened respectively.

The Hartree Term of B abc
By substituting from Eq.5 and 7 into Eq.17 then the Hartree term of the quantum third virial coefficient By taking the inverse Fourier transformation of the above equation and making use of the Wigner distribution function [4]; and using Equation (20) we obtain where the number density Eq.23 is vanishes for example for one component plasma such as the electron gas.

The Hartree-Fock Term of B abc
The polarization function is defined random phase approximation by ( ) ( ) By substituting from Eq.25 into Eq.18 we can get the Hartree-Fock term of the quantum third virial HF abc B as follow The inverse Fourier transform of this equation with the help of the Eq.22 gives , We assume that The calculation of ( ) By substituting from Eq.33 into Eq.31 then we get where, ( ) , , , , and T are the Fermi integral, the volume, the spin projection, the degeneracy parameter, the chemical potential, the Boltzmann's constant and the absolute temperature respectively.
Then we can written the third virial coefficient in the following form The analytical calculation of the integral from Eq.38 is evaluated by solving this integral by parts and using Gamma functions in the regions of small κλ (low densities, high temperatures) then we have By substituting from Eqs.13, 23, 35 and 39 into Eq.12 we get the excess free energy until the third virial coefficient for one component plasma as follow where K  is the quantum virial function which given by [4] and is the normalized thermal wave length.
Similarly, we can write the excess free energy until the third virial coefficient for two component plasma as follow

THE EQUATION OF STATE
Following the method of effective potentials developed by [9] we get the pressure from the excess free energy as follow , .
where id p is the ideal pressure and .
By substituting from Eq.40 into Eq.42 we can get the equation of state until third virial coefficient for one component plasma; ( Also, by substituting from Eq.41 into Eq.42 we get the equation of state until third virial coefficient for two component plasma In the numerical calculation, we let

DISCUSSIONS
To our knowledge there is no paper to calculate the third virial coefficient by using Green's function technique until now; this paper is the first paper to calculate the third virial coefficient in Hartree, Hartree-Fock approximation and the Montroll-Ward contribution by using the Green's function technique, and used it to calculate the quantum thermodynamic functions.In past the potential was used as the mean potential for two particles only so their results were until the second virial coefficient, but in this paper we used the potential as the sum of the mean potential of two and three particles so our results were evaluated until the third virial coefficient.Also the quantum thermodynamic functions until the third virial coefficient which are calculated by using the binary Slater sum are near at the classical limit only; they used the potential as the pair potential only and neglected the triplet potential so there results Is not exactly correct results.We considered only the thermal equilibrium plasma in the case of one and two component plasma by using Green's function method.We obtained the general formula of the third virial coefficient in Hartree-Fock approximation analytically (Eq.34).
As shown in Figures 1-3, we plotted the comparison between the excess free energy until the second virial coefficient for one and two component plasma of Ebeling et al. [11], Hussein et al. [1] until the third virial coefficient and our results until the third virial coefficient at .We noticed that the curves are small nearly for small values of ξ , this is due to the difference between Green's function technique which was used here and Slater Sum technique in Hussein et al. [2] and Ebeling et al. [11] which are given in Figures 1-3.In these figures we observe that there exist difference for small values of ξ for two comopnent plasma especially order 4  e .Also, we plotted the comparison between the pressure until the second virial coefficient for one and two component plasma of Kremp et al. [4], Ebeling et al. [11], Hussein et al. [1] until the third virial coefficient and our results until the third virial coefficient up to order 2 4  , e e and 6  e as shown in Figures 4-6.In Figures 5 and 6, the curves are far from each other for one and two component plasma.But in Figure 4, we observed that the red curve for Ebeling et    al. [11] and green curve for our result are nearly for some values of ξ for two component plasma.
v r r r G = − ∑ ∫ where, a f and aa v are Fermi functions and the Fourier transform of Coulomb potential which is defined by ( ) the dielectric constant where , r ε ε  are the vacuum and the relative dielectric constant.
wers of e α and using the spherical polar coordinates (see Appdenix A) then

35) 3 . 3 .is the inverse Debye length and 1 1 F
The Screened Contribution of B abc By substituting from Eqs.6 and 7 into Eqs.19 we can rewrite the screened quantum third virial S abc Buntil first three terms in the following form transform of binary potential then we get in the weakly degenerate case or low degeneracy limit denotes the confluent hypergeometric function.

Figure 1 .
Figure1.The comparison between the quantum excess free energy until the second virial coefficient[11] (black solid line),[1] (blue short dashed dotted line) until the third virial coefficient and our results (green long dashed line) up to order e 2 .

Figure 2 .Figure 3 .
Figure 2.The comparison between the quantum excess free energy until the second virial coefficient[11] (black solid line),[1] (blue short dashed dotted line) until the third virial coefficient and our results (green long dashed line) up to order e 4 .

Figure 4 .Figure 5 .
Figure 4.The comparison between the quantum pressure until the second virial coefficient [11] (black solid line), [4] (red long dashed dotted line) and [1] (blue short dashed dotted line) until the third virial coefficient and our results (green long dashed line) up to order e 2 .

Figure 6 .
Figure 6.The comparison between the quantum pressure until the second virial coefficient (red dotted line), [1] (blue short dashed dotted line) until the third virial and our result (green long dashed line) up to order e 6 . )