Ultrarelativistic Electron-Positron Plasma

# Ultrarelativistic Electron-Positron Plasma

Markus H. Thoma Max-Planck-Institut für extraterrestrische Physik, Gießenbachstr., 85748 Garching, Germany
Received: date / Revised version: date
###### Abstract

Ultrarelativistic electron-positron plasmas can be produced in high-intensity laser fields and play a role in various astrophysical situations. Their properties can be calculated using QED at finite temperature. Here we will use perturbative QED at finite temperature for calculating various important properties, such as the equation of state, dispersion relations of collective plasma modes of photons and electrons, Debye screening, damping rates, mean free paths, collision times, transport coefficients, and particle production rates, of ultrarelativistic electron-positron plasmas. In particular, we will focus on electron-positron plasmas produced with ultra-strong lasers.

###### pacs:
52.27.EpElectron-positron plasmas and 11.10.WxFinite-temperature field theory

## 1 Introduction

Plasmas, i.e. (partly) ionized gases, are considered to be the fourth state of matter after the solid, liquid, and gaseous states. The plasma state dominates in the observable Universe: 99% of the visible matter is in the plasma state, namely in the form of stars and hot interstellar and intergalactic gases. Plasmas emit light due to the excitation of atoms and ions and recombination. Plasmas can be produced by high temperatures, such as in the sun or fusion reactors, by electric fields (discharges), as used for illumination in neon tubes or in lightening, or by radiation, such as in the Crab nebula where the pulsar in the center emits synchrotron radiation.

Plasmas can be classified according to various aspects: Relativistic plasmas, e.g. the electron-positron plasma in a supernova explosion, are plasmas in which the thermal energy of the plasma particles is of the order of their rest mass energy or larger. Quantum plasmas, e.g. the degenerate electron component in a white dwarf, are plasmas in which the thermal de Broglie wave length is of the order of the interparticle distance or larger. Here is the thermal velocity of the particles. Strongly coupled plasmas, e.g. the ion component in white dwarfs, are plasmas in which the interaction energy between the particles is larger than their thermal energy. In non-relativistic plasmas this corresponds to a Coulomb coupling parameter .

Plasmas in Nature, like in comets, in aurorae, in the corona of the sun, in lightening, in flames, or in the sun, and in technology, like in discharges or fusion reactors, cover a wide range of pressures and temperatures. All of these plasmas are non-relativistic, classical, and weakly coupled systems.

Electron-positron plasmas (EPPs) are created in the presence of strong electric or magnetic fields or extremely high temperatures, where massive pair production sets in. For example in a supernova explosion temperatures up to K corresponding to will lead to an ultrarelativistic EPP. Also in the vicinity of magnetars, i.e. neutron stars with magnetic fields G, and in accretion disks around black holes EPPs show up.

Recently the possibility to create ultrarelativistic EPPs with high-intensity lasers ( W/cm) have been discussed. For example two opposite laser pulses hitting a thin gold foil will heat up the electrons in the foil up to several MeV leading to pair creation Shen ().

In the following we will discuss the properties of an ultrarelativistic EPP using quantum field theory (QED) at finite temperature. We will follow closely the review article Thoma1 (). We will not discuss the production mechanism and equilibration of the EPP here.

## 2 Field theoretic description of an electron-positron plasma

Throughout the paper we will use natural units, i.e. , as usual in quantum field theory, in which all units are given in powers of MeV. The conversion to conventional units can be achieved by MeV m from which follows. In these units the electron charge corresponding to a fine structure constant .

### 2.1 Equation of state

We will start with the equation of state of an EPP and compute it under the following assumptions:

1. ultrarelativistic EPP, i.e. ,

2. thermal and chemical equilibrium,

3. equal electron and positron density, i.e. vanishing chemical potential,

4. ideal gas, i.e. no interactions in the plasma,

5. infinitely extended, homogeneous, and isotropic EPP.

We will relax some of these assumptions in the following sections. According to these assumptions the distribution function of the electrons and positrons is given by the Fermi-Dirac distribution

 nF(p)=1ep/T+1 (1)

and of the photons by the Bose-Einstein distribution

 nF(p)=1ep/T−1, (2)

where the momentum is identical to the energy of the particles in the ultrarelativistic case. It should be noted that the photons are in equilibrium with electrons and positrons under the above assumptions, i.e. the system is actually an electron-positron-photon gas.

The particle and energy density can be calculated by integrating over the distribution functions. The particle number density of the electrons and positrons follows from integrating over the Fermi-Dirac distribution as

 ρeqe=gF∫d3p(2π)3nF(p)=3π2ζ(3)T3=0.37T3, (3)

where is the number of degrees of freedom corresponding to the electrons and positrons in the two spin states. Assuming a temperature of MeV, we find MeV m.

The photon density follows accordingly by integrating over the Bose-Einstein distribution with degrees of freedom corresponding to the two polarization states as .

The energy density of the electron-positron-photon gas is obtained from

 ϵeq = gF∫d3p(2π)3pnF(p)+gB∫d3p(2π)3pnB(p) (4) = 11π260T4=1.81T4,

where the photons contribute 36 % to the energy density. Here the Boltzmann law, , holds also for the fermions because we neglected their masses.

For MeV we find J m. In a volume of m (corresponding to the size of a neutron star) the total thermal energy of the EPP is J, which corresponds to about 10% of the entire energy (without neutrinos) released in a supernova type II explosion. In a volume of 1 m there is still an energy of J contained.

The Coulomb coupling parameter of the EPP, which is a measure for the non-ideal behavior of a plasma Ichimaru (), is given by , where m is the interparticle distance. For MeV we find which shows that the EPP is a weakly coupled plasma. Therefore the ideal gas results for the equation of state derived above are a good approximation. After all, interactions in the EPP play an important role, for example, for the collective behavior of the plasma as discussed in the next section and for equilibration of the plasma. Obviously, the interactions can be treated by perturbation theory.

### 2.2 Collective Phenomena

Collective effects in a plasma are associated with long-range interactions within the plasma. Important examples are Debye screening and plasma waves. In non-relativistic ion-electron plasmas Lifshitz () plasmas these phenomena can be described by classical transport theory (Vlasov equation). For example, the electron plasma frequency reads

 ωpl=√4πe2ρeme (5)

and the Debye screening length due to the electrons in the plasma

 λD=√kBTe4πe2ρe, (6)

where is the electron number density, the temperature of the electron component, and the electron mass. In an ultrarelativistic plasma with the masses can be neglected and the important scales are the temperature , called the hard scale, and the soft scale , which determines the collective phenomena as we will see below.

Interactions between relativistic electrons and positrons can be treated by using perturbative QED. This corresponds to an expansion in the fine structure constant . Most conveniently Feynman diagrams are considered from which via Feynman rules quantities such as scattering cross sections, decay and production rates, or life times can be calculated directly. In an EPP the interactions take place in the presence of a heat bath. Hence we have to consider QED at finite temperature. For this purpose the Feynman rules are extended to finite temperature which can be achieved by using the imaginary or real time formalism Kapusta (); Landsmann (). The calculations are similar to the ones done already within the last 30 years for the properties of the quark-gluon plasma using perturbative QCD at finite temperatures. As a matter of fact, many results from the quark-gluon plasma (see e.g. Ref.Thoma2 ()) can be directly carried over to the EPP.

An important quantity is the polarization tensor or photon self-energy. The lowest order diagram for the polarization tensor is shown in Fig.1. Assuming the external four momentum to be soft, i.e. the frequency and to be much smaller than , and the internal loop momenta to be hard, an analytic result can be found using the real or imaginary time formalism Klimov (); Weldon1 ()

 ΠL(ω,k) = −3m2ph(1−ω2klnω+kω−k), ΠT(ω,k) = 32m2phω2k2[1−(1−k2ω2)ω2klnω+kω−k] (7)

where is called the effective photon mass. For MeV we get MeV.

The crucial quantity from which the collective phenomena are derived is the dielectric tensor relating the macroscopic electric field in the medium to the external field (), i.e. in momentum space

 Di(ω,k)=∑jϵij(ω,k)Ej(ω,k). (8)

In the case of an isotropic medium it depends only on and on and has two independent components

 ϵij(ω,k)=ϵT(ω,k)(δij−kikjk2)+ϵL(ω,k)kikjk2. (9)

The dielectric tensor is closely related to the polarization tensor or photon self-energy by Elze ()

 ϵL(ω,k) = 1−ΠL(ω,k)k2, ϵT(ω,k) = 1−ΠT(ω,k)ω2, (10)

where and are the longitudinal and transverse components of the polarization tensor, respectively.

The dielectric functions following from (10) and (7) can also be derived from the classical Vlasov equation together with the Maxwell equations Silin (), since the high-temperature limit corresponds to the classical limit.

The dispersion relations of collective plasma modes, i.e. propagation of electromagnetic waves in the plasma, can be found by using the Maxwell equation, leading to

 ϵL(ω,k) = 0, ϵT(ω,k) = k2ω2. (11)

Combining (10), (7), and (11) gives the dispersion relations of the transverse as well as longitudinal plasma waves as shown in Fig.2. The longitudinal branch, which does not exist in vacuum, is called plasmon as in the case of non-relativistic plasmas. The transverse branch does not play a role in non-relativistic plasmas but is equally important as the longitudinal one in relativistic plasmas. Both branches start at the plasma frequency . Consequently the collective plasma waves have soft momenta of the order . At high momenta the transverse mode approaches the free dispersion , corresponding to a real photon in vacuum, whereas the longitudinal mode disappears, i.e. its spectral strength is exponentially suppressed. For MeV we find s.

Another important quantity which can be derived from the polarization or dielectric tensor is the Debye screening length, entering the Yukawa potential of a heavy, non-relativistic test charge in the EPP. The Debye screening length is given by the static limit of the longitudinal component of the polarization tensor , leading to , which is m for MeV.

Finally from (7) we see that the polarization tensor and the dielectric function become imaginary for , i.e. below the light cone , corresponding to Landau damping. We also observe that the plasma waves calculated at lowest order perturbation theory are undamped since they are located at .

A complete new phenomenon that does not appear in non-relativistic plasmas is the existence of fermionic plasma waves because all fermion masses are much too large in the non-relativistic case. Their dispersion relations follow from the pole of the electron propagator containing the electron self-energy. Using again the high temperature approximation for the one-loop electron self-energy of Fig.3 leads to (, ) Klimov (); Weldon2 ()

 Σ(P)=−a(p0,p)Pμγμ−b(p0,p)γ0 (12)

with

 a(p0,p) = 14p2[tr(PμγμΣ)−p0tr(γ0Σ)], b(p0,p) = 14p2[P2tr(γ0Σ)−p0tr(PμγμΣ)], (13)

where the traces over the matrices are given by

 tr(PμγμΣ) = 4m2F, tr(γ0Σ) = 2m2F1plnp0+pp0−p (14)

with the effective electron mass , which is 1.1 MeV at MeV.

The full electron propagator in the helicity representation is given by Braaten1 ()

 S⋆(P)=12D+(P)(γ0−^p⋅γ)+12D−(P)(γ0+^p⋅γ), (15)

where

 D±(P)=−p0±p+14p[±tr(PμγμΣ)−(±p0−p)tr(γ0Σ)]. (16)

The dispersion relations following from the pole of this propagator are shown in Fig.4. Again two branches show up, one with a positive ratio of the helicity to chirality () following from , the other one with a negative ratio () following from , called plasmino Braaten1 (). The plasmino branch , which does not exist in vacuum, shows an interesting behavior, namely a minimum at , which may lead to van Hove singularities Braaten1 (); Peshier (). Whether these van Hove singularities will lead to observable effects in the EPP, e.g. in the electron spectrum, is a very interesting question which should be investigated in detail. It could open the exciting possibility to observe a new collective plasma wave, the plasmino, experimentally in a laser induced EPP.

### 2.3 Transport properties

Now we want to consider the interaction and properties of particles in the plasma with hard momenta, i.e. of the order of or larger. In particular we are interested in damping and transport rates, mean free paths, collision times, energy losses of these particles and other transport properties such as the shear viscosity of the EPP.

It was shown by Braaten and Pisarski Braaten2 () that a consistent treatment of gauge theories such as QED at finite temperature, i.e. for obtaining results that are gauge independent, infrared finite, and complete to leading order, require the use of an effective perturbation theory using resummed Green functions based on the hard thermal loop (HTL) approximation (HTL resummation technique).

The damping rate of an electron or positron in the EPP is defined as the imaginary part of the dispersion relation . To lowest order it follows from the elastic scattering diagram of Fig.5. In the case of a hard electron or positron with momenta of the order of or higher it exhibits a quadratic infrared (IR) divergence which can be reduced to a logarithmic one using a HTL resummed photon propagator. This logarithmic singularity is expected to be cut-off by higher order contributions leading to Thoma2 ()

 γe=e2T4πln1e (17)

within logarithmic accuracy, i.e. the constant under the logarithm is not determined. For MeV we obtain keV, which is much smaller than MeV, showing that the EPP is not overdamped.

Physically more important are the transport rates which are related to the mean free path and collision time of electrons and positrons in the EPP. They differ from the damping rate in cutting off the long range interactions with small scattering angles by a factor under the integral defining the rate Lifshitz (). This leads to an improvement of the IR behavior (logarithmic instead of quadratic singularity in perturbation theory) and a finite result using the HTL method. Logarithmic divergent quantities can be treated consistently by splitting them into a soft part and a hard part, where the soft part is calculated using the HTL resummation technique Braaten3 (). For the transport rate we find to logarithmic accuracy

 Γe=e4T33πsln1e, (18)

where the Mandelstam variable is the square of the sum of the four momenta of the incoming particles in the scattering diagram of Fig.5. For thermal particles we replace by its thermal average , where . Assuming again MeV, we get keV.

The mean free path and collision time of the plasma particles (electrons and positrons) are given by the inverse of the transport rate , leading to nm and s at MeV.

In a non-relativistic plasma the shear viscosity can be estimated from elementary kinetic theory as Reif ()

 ηi=13∑iρi⟨pi⟩λmfpi (19)

where the sum is performed over the various components of the system. In an relativistic plasma the coefficient 1/3 should be replaced by 4/15 deGroot (). Using the mean free path following from (18), the density of (3), and the thermal momentum , the shear viscosity is given by (within logarithmic accuracy)

 ηe=55.8T3e4ln(1/e). (20)

At MeV the shear viscosity coefficient is Pa s.

Another quantity of interest in a plasma is its stopping power or the energy loss of an energetic particle in the plasma. There are two contributions, namely the energy loss by collisions and the radiative one by bremsstrahlung. In a relativistic plasma the latter one becomes important. The collisional energy loss is given by the mean energy transfer divided by the mean free path leading to Braaten4 ()

 dEdx=1v∫dγω (21)

where is the particle velocity, the damping or interaction rate proportional to the plasma density and the collision cross section, and the energy transfer from the energetic particle to the plasma particle in the collision. Using for the collision cross section the lowest order diagrams in Fig.6, the collisional energy of a muon with mass in an EPP has been calculated by Braaten and Thoma Braaten4 () applying the HTL resummation technique

 dEdx=e4T224π(1v−1−v22v2ln1+v1−v)(lnEM+ln1e+A(v)), (22)

where is a slowly varying function of the muon velocity between 1.3 and 1.5.

The collisional energy loss of an electron with energy is approximately given by Braaten4 ()

 dEdx=e4T248πln15.3Ee2T. (23)

This leads to an energy loss of 200 MeV/nm for an electron or positron with an energy of MeV at MeV, showing that such an electron is stopped (thermalized) within a fraction of a nanometer. So far no calculations of the radiative energy loss in an EPP have been performed to our knowledge.

The damping rate of a photon in an EPP follows from the diagram in Fig.7, where a HTL resummed electron propagator has to be used in case of soft momenta of the exchanged electron (positron). In contrast to the electron damping rate, the photon rate is infrared finite using the HTL method due to the presence of an electron propagator in Fig.7 instead of the photon propagator in Fig.5. Hence there is no need to cut off the long range interaction introducing a transport cross section. The result for a photon with energy reads Thoma3 ()

 Γph=e4T264πEln3.88Ee2T. (24)

The mean free path and the collision time of photons in an EPP are given by . For a thermal photon with the mean momentum at MeV the mean free path nm and the collision time s follow. Actually the damping rate given in (24) is a lower limit as higher order effects will enlarge it. As a matter of fact, the photon production rate in a QGP, which is the inverse process of the damping rate Thoma3 (), was shown to be about a factor of 2 larger taking bremsstrahlung into account Arnold1 ().

For the viscosity of the photon component using the above mean free path, the photon density (see above), and the mean photon energy we find

 ηph=48.7T3e4ln(3.27/e) (25)

corresponding to Pa s at MeV. Hence the viscosity of the EPP has similar contributions from the electrons and photons.

A more advanced calculation of the total viscosity of the EPP based on the Kubo formula yields within logarithmic accuracy Arnold2 ()

 η=188T3e4ln(1/e). (26)

This result is about a factor of 1.5 larger than the one presented here based on the elementary kinetic theory, which is typically valid within a factor of 2 Reif ().

A summary of the QED results of the EPP properties discussed above is presented in the table below.

Quantity Formula Value at MeV
Electron-Positron Density m
Photon Density m
Electron-Positron Energy Density J m
Photon Energy Density J m
Total Energy Density J m
Thermal electron momentum 31 MeV
Thermal photon momentum 28 MeV
Interparticle distance m
Coulomb Coupling Parameter
Effective Photon Mass 1 MeV
Plasma Frequency s
Debye Screening Length m
Effective Electron Mass 1.1 MeV
Electron Damping Rate 86 keV
Electron Transport Rate 0.54 keV for
Photon Damping Rate 0.70 keV for
Electron Mean Free Path 0.37 nm
Photon Mean Free Path nm
Electron Collision Time s
Photon Collision Time s
Electron Viscosity Pa s
Photon Viscosity Pa s
Total Viscosity Pa s
Electron Energy Loss 200 MeV/nm for MeV

### 2.4 Non-equilibrium and finite chemical potential

EPPs produced in strong laser fields are probably not in complete equilibrium. For example, it has been predicted by Shen and Meyer-ter-Vehn Shen () that a positron density of about m at a temperature of 10 MeV can be reached. This density deviates from the equilibrium density (3) by 12 orders of magnitude. In the following we will therefore assume that the EPP produced by lasers is in thermal but not in chemical equilibrium. Then we can replace the distribution functions for the electrons and positrons by Fermi-Dirac distributions multiplied by a fugacity factor describing the deviation from chemical equilibrium, . This assumption has been used for example for describing the chemical equilibration of the QGP in ultrarelativistic heavy-ion collisions Biro (). The fugacity is given by the ratio of the experimental to equilibrium particle density, since the experimental density follows from integrating over the non-equilibrium distribution, i.e.

 ρexp=gF∫d3p(2π)3λnF(p)=λρeq⇒λ=10−12. (27)

Using the real time formalism, QED perturbation theory and the HTL method can also be extended to non-equilibrium situations like the one discussed above Carrington (). For example, the effective photon mass is given now by

 m2ph=4e23π2∫∞0dppfF(p). (28)

For MeV we then find for the non-equilibrium photon mass eV and the plasma frequency Hz. The Debye screening length in such an EPP is m. In order to speak of a plasma with collective behavior its dimension should be much larger than , i.e. at least of the order of 1 m.

Furthermore, an anisotropic EPP can also be described by quantum field theoretic methods Mrowczynski (). In this case instabilities can occur Randrup ().

Finally a possible difference between the positron density and the electron density can be treated by introducing a finite chemical potential , i.e. using the distribution

 nF(p)=1e(p±μ)/T+1 (29)

for the electrons (negative sign) and positrons (positive sign). Such a difference comes from the fact that the laser produced EPP is embedded in a hot electron and cold ion background of the target. Therefore there will be an excess of electrons over positrons in the hot, relativistic EPP. The methods described above, such as the HTL resummation, can be generalized easily to this case Vija (). For example, the energy density is given by

 ϵeq=11π260T4+12T2μ2+14π2μ4 (30)

or the effective photon energy by

 m2ph=e2T29(1+3μ2π2T2). (31)

### 2.5 Particle Production

At high temperatures above 10 MeV also other particle species will be produced, e.g. muons with a mass of MeV. Their rate follows to lowest order from the diagram in Fig.8 (Born term). We assume that holds. The first inequality implies that the electron mass can be put to zero and the latter inequality implies that muons are not equilibrated. Then the muon production rate to lowest order () is given by (for more details see Ref.Thoma1 ())

 dNd4xd4p=α224π4(1+2m2μM2)(1−4m2μM2)1/2Tp 1exp(E/T)−1ln1+exp[−(E+p)/(2T)]1+exp[−(E−p)/(2T)], (32)

where is the invariant mass of the virtual photon , its energy and its momentum. Because of the rate is suppressed exponentially for temperatures below .

In order to obtain the spectrum from this formula one has to integrate over the space-time volume, taking into account the space-time evolution by using, for example, a hydrodynamical model. The total muon yield then follows from integrating the spectrum over the energy and momentum of the virtual photon.

At temperatures above 10 MeV also hadron production becomes important, in particular pion production Kuznetsova ().

## 3 Conclusions

The aim of this presentation is the prediction of properties of ultrarelativistic thermalized EPPs produced in laser fields or supernovae. The ultrarelativistic EPP is a weakly coupled system. Therefore its equation of state can be described in first approximation by an ideal ultrarelativistic gas. Interactions within the EPP can be described by perturbative QED at finite temperature. In this way collective phenomena (plasma waves, Debye screening) and transport properties (damping rates, mean fee paths, relaxation times, production rates, viscosity, energy loss) can be computed. A complete new phenomenon, fermionic plasma waves (plasmino), which is absent in non-relativistic plasmas, might be observable by van Hove singularities. The deviation from chemical equilibrium, as expected for laser produced EPPs, can also be treated by extending QED to the non-equilibrium case.

Acknowledgment: I would like to thank D. Habs, J. Rafelski, I. Kouznetsova, and G. Moore for helpful discussions and hints. The author acknowledges the support by the European Commission under contract ELI pp 212105 in the framework of the program FP7 Infrastructures-2007-1.

## References

• (1) B. Shen and J. Meyer-ter-Vehn, Phys. Rev. E 65, 016405 (2001)
• (2) M.H. Thoma, arXiv:0801.0956, Rev. Mod. Phys. (to be published)
• (3) S. Ichimaru, 1982, Rev. Mod. Phys. 54, 1017 (1982)
• (4) E.M. Lifshitz and L.P. Pitaevskii, Physical Kinetics (Pergamon Press, Oxford, 1981)
• (5) J.I. Kapusta, Finite Temperature Field Theory, (Cambridge University Press, New York, 1989).
• (6) N.P. Landsmann and C.G. van Weert, Phys. Rep. 145, 141 (1987)
• (7) M.H. Thoma, in Quark-Gluon Plasma 2, edited by R.C. Hwa, (World Scientific, Singapore, 1995), p.51, hep-ph/9503400
• (8) V.V. Klimov, Zh. Eksp. Teor. Fiz. 82, 336 (1982) [Sov. Phys. JETP 55, 199 (1982)]
• (9) H.A. Weldon, Phys. Rev. D 26, 1394 (1982)
• (10) H.T. Elze and U. Heinz, Phys. Rep. 183, 81 (1989)
• (11) V.P. Silin, J. Exptl. Theoret. Phys. (U.S.S.R.) 38, 1577 (1960) [Sov. Phys. JETP 11, 1136 (1960)]
• (12) H.A. Weldon, Phys. Rev. D 26, 2789 (1982)
• (13) E. Braaten, R.D. Pisarski, and T.C. Yuan, Phys. Rev. Lett. 64, 2242 (1990)
• (14) A. Peshier and M.H. Thoma, Phys. Rev. Lett. 84, 841 (2000)
• (15) E. Braaten and R.D. Pisarski, Nucl. Phys. B 337, 569 (1990)
• (16) E. Braaten and T.C. Yuan, Phys. Rev. Lett. 66, 2183 (1991)
• (17) F. Reif, Fundamentals of Statistical and Thermal Physics, (McGraw-Hill, New York, 1965)
• (18) S.R. de Groot, W.A. van Leeuwen, and C.G. van Weert, Relativistic Kinetic Theory, (North-Holland, Amsterdam, 1980)
• (19) E. Braaten and M.H. Thoma, Phys. Rev. D 44, 1298 (1991)
• (20) M.H. Thoma, Phys. Rev D 51, 862 (1995)
• (21) P. Arnold, G.D. Moore, and L.G. Yaffe, JHEP 0112, 009 (2001)
• (22) P. Arnold, L.G. Yaffe, and G.D. Moore, JHEP 0011, 001 (2000)
• (23) T.S. Biró, E. van Doorn, B. Müller, M.H. Thoma, and X.N. Wang, Phys. Rev. C 48, 1275 (1993)
• (24) M.E. Carrington, H. Defu, and M.H. Thoma, Eur. Phys. J. C 67, 347 (1999)
• (25) S. Mrówczyński and M. H. Thoma, Phys. Rev. D 62, 036011 (2000)
• (26) J. Randrup and S. Mrówczyński, Phys. Rev. C 68, 034909, (2003)
• (27) H. Vija and M.H. Thoma, Phys. Lett. B 342, 212 (1995)
• (28) I. Kuznetsova, D. Habs, and J Rafelski, Phys. Rev. D 78, 014027 (2008)
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters