Schwinger Mechanism with Energy Dissipation in Glasma
Initial states of glasma in high energy heavy ion collisions are longitudinal classical color electric and magnetic fields. Assuming finite color electric conductivity, we show that the color electric field decays by quark pair production with the life time of the order of , i.e. the inverse of the saturation momentum. Quarks and anti-quarks created in the pair production are immediately thermalized as far as their temperature is lower than . Namely, a relaxation time of the quarks to be thermalized is much shorter than when . We also show that the quarks acquire longitudinal momentum of the order of by the acceleration of the electric field. To discuss the quark pair production, we use chiral anomaly which has been shown to be very powerful tool in the presence of strong magnetic field.
pacs:12.38.-t, 24.85.+p, 12.38.Mh, 25.75.-q, 12.20.-m, 12.20.Ds
Schwinger mechanism, Chiral Anomaly, Color Glass Condensate
Initial states of color gauge fields (glasma) produced immediately after high energy heavy-ion collisions have recently received much attention. The gauge fields are longitudinal classical color electric and magnetic fields; their field strengths are given by the square of the saturation momentum . The presence of such classical gauge fields has been discussed on the basis of a fairly reliable effective theory of QCD at high energies, that is, a model of color glass condensate (CGC)gl (); cgc (). It is expected that the decay of the glasma leads to thermalized quark gluon plasma (QGP).
The glasma is homogeneous in the longitudinal direction and inhomogeneous in the transverse directions. Hence, we may view that it forms electric and magnetic flux tubes extending in the longitudinal direction. In the previous papersiwa (); itakura (); hii () we have shown a possibility that the famous Nielsen-Olesen instabilitynielsen () makes the color magnetic field decay. The possibility has been partially confirmed by the comparison between our results and numerical simulationsvenugopalan (); berges (). Such decay is a first step toward the generation of QGP.
On the other hand we have also discussediwazaki1 (); iwazaki2 (); old () the decay of the color electric field; the decay is caused by Schwinger mechanismschwinger (), that is, the pair production of quarks and anti-quarks. The mechanism has been extensively exploredtanji () since the discovery of Klein paradox. Among them, the pair production in the expanding glasma has been discussedlap (). A new feature in the glasma is that it is composed of not only electric field but also magnetic field parallel to the electric field. Such a feature has also been explored. In particular, recently there are studies of the back reaction of the particles on the electric field under the presence of the magnetic fieldiwazaki2 (); tanji (); suga (). The back reaction is essentially important for the decay of the electric field.
Our originalityiwazaki2 () for the investigation of the decay is to use chiral anomaly. As is well known, the anomaly is effective when collinear magnetic and electric fields are present. This is the very situation in the glasma. When we use the chiral anomaly, we can discuss Schwinger mechanism without detail calculationstanji (); lap (); suga () of wave functions but simply by solving classical anomaly equation and Maxwell equations. In particular, when the strong magnetic field is present, the anomaly is much simplified because the quarks are effectively massless and only relevant states are ones in the lowest Landau level. ( Both and in the glasma are much larger than mass of quarks. ) Since the motions of the quarks in transverse directions are frozen, only possible motion is along the longitudinal direction. Thus, the anomaly equation is reduced to the one in two dimensional space-time. With the simplification, we can find important quantities even in complicated realistic situations for which the investigations have not yet performed. Actually, we have discussed the decay of axial symmetric electric flux tube by taking account of the azimuthal magnetic field around the tube. The field is generated by the current carried by the pair created quarks and anti-quarks. Although the electric field loses its energy by the pair creation and the generation of the azimuthal magnetic field, it never decays. The field oscillating with time propagates to the infinity in the axial symmetric wayiwazaki2 (). This is because the quarks are free particles and there is no energy dissipation. ( In the case of homogeneous electric field, the field simply oscillates with time. )
In this paper we wish to discuss the decay of the color electric field by taking account of energy dissipation in heat baths. The dissipation arises due to the presence of finite electric conductivity. Namely, the pair production generates color electric current , which dissipates its energy owing to the interaction between the quarks and the surrounding; the surrounding is composed of quarks and gluons. Actually, we assume Ohm law with electric conductivity . The conductivity is calculated by using Boltzmann equation in the relaxation time approximation. In the approximation a relaxation time is obtained by calculating electron’s scattering rates. Then, we can show that the quarks are thermalized immediately after their production as far as their temperature is much smaller than ; the relaxation time of a slightly deformed momentum distribution of quarks becoming the equilibrium Fermi distribution is much shorter than the life time of the field. As numerical calculations have showntanji (), the longitudinal momentum distribution of the free particles produced in the vacuum is almost equal to the equilibrium one, that is Fermi distribution at zero temperature. Thus, even in non zero but low temperature, the momentum distribution is nearly equal to the equilibrium one. Our relaxation time approximation in Boltzmann equation may holds in such a situation. Therefore, owing to the energy dissipation by the scattering between electrons and positrons, the electric field decays and never oscillates.
For simplicity, we examine homogeneous electric and magnetic fields of U(1) gauge theory instead of SU(3) gauge theory. Thus, we use terminology of electron or positron instead of quark or anti-quark. The generalization of our results to the case of SU(3) gauge theory is straightforward done simply by assuming maximal Abelian group of SU(3)tanji (). We assume that both the electric and magnetic fields are much larger than the square of the electron mass. Thus, they are taken to be massless.
In the next section we explain how the chiral anomaly is useful for the discussion of Schwinger mechanism. We apply the anomaly to the discussion of the pair production with energy dissipation in the section . In the discussion we use the conductivity of electron-positron gas under strong magnetic field. The detail calculation of it is shown in the appendix.
Ii Chiral anomaly and pair production
We first make a brief review of the utility of the chiral anomaly. As we will show below, only by simple calculations we can obtain most of important results of Schwinger mechanism, which was derived previously by detail calculations.
Under the homogenous magnetic field , massless fermions with charge occupy the states characterized by Landau level and momentum parallel to the magnetic field. Their energies are given by
where integer denotes Landau level. The term of “parallel” (“anti-parallel”) implies magnetic moment parallel (anti-parallel) to . The magnetic moment of electrons (positrons) is antiparallel (parallel) to their spin. Thus, electrons (positrons) with spin anti-parallel (parallel) to can have zero energy states in the lowest Landau level; their energy spectrum is given by . On the other hand, the other states cannot be zero energy states. Hence only relevant states for the pair production are states with energy in the limit of strong magnetic field. Hereafter, we consider only such states.
The equation of the chiral anomaly for the states leads to
where ( ) denotes number density of right ( left ) handed chiral fermions; in which the expectation value is taken by using the states in the lowest Landau level. We have used spatial uniformness of the chiral current , that is, .
We note that electrons move to the direction anti-parallel to while positrons move to the direction parallel to after their pair production. Therefore, both the electrons and positrons have right handed helicity when being parallel to . ( Right or left handed helicity means the state with momentum parallel or antiparallel to spin, respectively.) Since chirality is equivalent to helicity of massless fermions, and where is the number density of electrons. ( When is anti-parallel to , and . ) Therefore, the anomaly equation describes the pair production rate of electron and positrons under and ,
For example, we suppose that the electric field is switched on at . Then, solving this equation with initial condition , we can obtain the number density of electronstanji (); . When is constant, increases linearly with time. In this estimation, the back reaction of the produced particles to the electric field is not taken into account.
When we take account of the back reaction, we can find that the electric field oscillates with time. The back reaction is that the energy of electric field decreases owing to the acceleration of the charged particles. It can be included by solving a Maxwell equation, along with the anomaly equation. To solve the equations, we need to express the electric current in terms of . This can be done by considering the energy conservation. Namely, the total energy of the electric field and the charged particles is conserved,
where the energy density of electrons or positrons is given such that with Fermi momentum for . Namely, we assume the energy of free electrons and positrons. Their momentum distributions are supposed to be for positrons or for electrons. We have used the formula in eq(4). Thus, from the equation we obtain the electric current for . Obviously, there is no energy dissipation.
Using the Maxwell equation and the anomaly equation we derive the equation of the number density,
We can see that after the switch on of , gradually decreases with time, while increases. Eventually, these quantities oscillate with the frequency ( plasma oscillation ) as shown in Fig.1. The oscillation can be understood in the following.
When ( ) parallel to is switched on at with , the number of the particles increases while decreases. When vanishes, takes the maximum. After that changes its sign, namely, becomes negative ( anti-parallel to .) Then, the particles are forced to move to the direction inverse to the one in the case of . Thus, the pair annihilation occurs and begins to decrease. The energy of the particles is transmitted to the energy of the electric field. Actually, according to the anomaly equation, becomes large with time. When vanishes, takes the maximum ( ). Subsequently, the pair creation occurs and increases according to the anomaly equation(3) but with . This is a plasma oscillation shown in the figure. The energy of the electric field is transmitted to the particles and then the energy of the particles goes back to the electric field. The plasma oscillation is a natural consequence of the absence of the dissipation.
We should make a comment on the scales involved in the discussion. The above result ( plasma oscillation ) holds in the limit of where only the states in the lowest Landau level are relevant. On the other hand when , states in higher Landau levels become important since the energy of electrons can be much larger than . In the glasma we have both and ( square of electron mass ). But states in higher Landau levels are still not important according to the following reason. When some of electrons in the lowest Landau level are accelerated and their energies reach the energy , they may make a transition to a state in a higher Landau level with energy . Such electrons are the ones produced in the very early stage of the pair production. The number of the electrons is still a fraction. Most of electrons whose energies are less than , occupy the lowest Landau level. It has recently been shownsuga () that the decay of the electric field is mainly caused by particles in the lowest Landau level. Thus, our result is approximately correct.
In this way, we obtain the number density of electrons simply by using the chiral anomaly. It is important to note that the anomaly equation involves all of quantum effects. Thus, by simply solving the classical anomaly equation, we can obtain the quantities resulting from the purely quantum effect, i.e. Schwinger mechanism. In the standard manipulationtanji () of Schwinger mechanism, we need to find wave functions of electrons under the electric and magnetic fields in order to obtain such quantities. When the time dependence of the electric field is complicated, it is very difficult to obtain the wave functions. On the other hand, the anomaly equation is a very simple and powerful tool for investigating Schwinger mechanism.
Iii Pair production with energy dissipation
Up to now, we have not taken into account energy dissipation of electrons and positrons produced by the electric field. We proceed to analyse the Schwinger mechanism in heat bath, especially, the case with energy dissipation. We assume that the particles are immediately thermalized after their production. ( Later we see that relaxation time, within which slightly deformed particle distribution becomes equilibrium distribution, is much shorter than the typical scale as far as the system is in low temperature. As numerical calculations have showntanji (), the momentum distribution of the produced particles is nearly equal to the Fermi distribution at zero temperature. ) Thus, the number density is given by
with , where is the inverse of temperature and is chemical potential of electrons or positrons.
The number density depends on two parameters, and , which may be functions of time. Hereafter, we consider only the case that the chemical potential vanishes due to non-conservation of electron number. It is caused by frequent pair annihilation or creation of electrons and positrons after their production. Since initial states have no lepton number, the chemical potential vanishes. ( This corresponds to the situation in the glasma where the pair annihilation or creation of quarks and anti-quarks frequently take place after the pair creation by Schwinger mechanism. ) Thus, the number density is only a function of the temperature and is given by . Similarly, the energy density is given by
with . Since the temperature represents an average energy of the particles, it is reasonable that we have ; the energy density is equal to the number density times the average energy of a particle.
The electric field produces electrons and positrons and generates their electric currents. When the electric current flows, the current dissipates its energy due to the finite electric conductivity of these particles. Thus, the electric field dissipates its energy. Here we may assume that the electric current satisfies the Ohm law where denotes electric conductivity. Obviously, the Ohm law represents the energy dissipation of the electric field .
As we will show in the appendix, the electric conductivity of the charged particles in the lowest Landau level is given by
where denotes the momentum distribution of the charged particles. The conductivity depends on the relaxation time ; the time needed for a momentum distribution to become the equilibrium one . The relaxation is realized by scattering of the charged particles so that it can be expressed by their scattering amplitudes. The detail calculation shows that in the limit of low temperature . ( This fact makes valid our assumption of the immediate thermalization in the low temperature where the relaxation time is small. ) Thus, the conductivity depends on the temperature such as in the limit ( is a numerical constant ).
Thus, using the Maxwell equation and the anomaly equation ( ), we find that
Setting , we obtain the equation of the electric field
with initial conditions and . The number density can be obtained by integrating the solution such that .
The equation (10) describes a motion of a classical particle with the potential energy ; denotes the coordinate of the particle. Starting at , the particle goes down to along the slope of the potential. It means that . Namely, the electric field vanishes with the life time of the typical scale . ( The life time of the electric field is sufficiently short to be consistent with phenomenologyhirano () of QGP. ) Since , the number density or the temperature monotonically increases with time and reaches its maximum. This is contrasted with the case of no dissipation in which both and oscillate. These behaviors are shown in the Fig.2.
We can roughly estimate the attainable number density , temperature and average energy of the particles in the pair creation,
where we set which corresponds to gauge coupling constant in QCD with the energies at RICH or LHC.
We should make a comment. The result has been derived in the limit of the low temperature. The relaxation time is sufficiently short in the limit so that we can use the equilibrium momentum distribution. The temperature increases with time from zero after the switch-on of the electric field. Therefore, the result is reliable in the early stage of the pair creation. However, when it becomes comparable to , the assumption does not hold. Then, we cannot use the formula , and the equilibrium momentum distribution. Hence, the result of the saturation in the number density at the late stage of the evolution cannot be reliable. But, as the temperature increases and becomes comparable to , the conductivity is almost constant in time, Then, we find that from the equation and that from the anomaly equation. The number density is saturated and never oscillate. ( To obtain this result we do not use the equilibrium momentum distribution. Namely, the assumption of the particles being immediately thermalized is not necessary. ) Actually, when the temperature or the average energy of the particles is comparable with , the electric conductivity may depends weakly on the time because only dimensional quantity in is or the average energy. Consequently, we find that the number density ( temperature ) gradually increases in the early stage and is saturated at the late stage. It never oscillates unlike the case with no energy dissipation.
To put it another way, in the early stage when the pair creation starts, the particles are thermalized soon because the relaxation time is short. Thus, the temperature can be defined. The temperature or the average energy of the particles is much less than . But in the stage when relaxation time is not sufficiently short for the particles to be thermalized soon, the temperature cannot be defined. In the stage all of the typical energy scales are of the order of . Thus, the electric conductivity is also of the order of . This leads to the above result of the saturation . The stage is not early stage but may be late stage in the pair creation. In this way we understand that the results shown in Fig.2 are reliable.
Iv discussion and summary
We have shown that the color electric field decays into the quarks and anti-quarks under the color magnetic field. They occupy only the states in the lowest Landau level. We have also discussed that the average energy of the particles is of the order of the saturation momentum . Namely, the Fermi momentum, which is the longitudinal momentum, is given by . The particles obtain the longitudinal momentum with the acceleration by the electric field. We remind you that the glasma has only transverse momentum of the order of and has no longitudinal momentum owing to the homogeneity in the direction. Thus, initially there is anisotropy in momentum distribution. But owing to the decay of the electric field, the particles have longitudinal momenta of the order of . Thus, the longitudinal and transverse momenta become almost identical to each other. This is a step toward the isotropy of the momenta or thermalization of the quarks.
The life time of the electric field has been found to be of the order of , while that of the magnetic field by a decay mechanism of Nielsen Olesen instability is much longervenugopalan () than . We may speculate that the magnetic field also decay rapidly due to Schwinger mechanism. This is because azimuthal electric field is generated by the expansion of the longitudinal magnetic flux tube and this electric field decays by Schwinger mechanism. The expansion in the transverse direction of the magnetic flux tube is simply a result of Maxwell equationsiwa (); itakura (). Similarly, azimuthal magnetic field is also generated by the expansion of the longitudinal electric flux tube. Thus, both azimuthal magnetic and electric fields are generated around the longitudinal flux tubes. Then, the azimuthal electric field decays in the way as we have shown. The field strengths of these longitudinal and transverse ( azimuthal ) fields have been shownlappi () to become the same order of magnitudes at about the time after the high energy heavy ion collisions. In this way, the longitudinal magnetic flux may decay rapidly through the decay of the azimuthal electric field generated by the expansion of the magnetic flux. Such a rapid decay of the magnetic field as well as the electric field is favorable for establishing early thermalizationhirano () of QGP.
To summarize, we have shown how the electric field decays by Schwinger mechanism in the presence of the energy dissipation. Namely, the electric current carried by electrons and positrons dissipate its energy owing to the Ohm law, which results in the increase of the temperature with time. Our result holds only in the limit of low temperature, in which the relaxation time has been shown to be much short; the thermalization of produced particles can be realized immediately after their production. However, we have also discussed that our result is approximately correct even when the temperature or the average energy of the particles is of the order of .
Appendix A electric conductivity
We show that the electric conductivity behaves as in the limit of the low temperature . We calculate in the relaxation time approximation of Boltzmann equation,
with , where we have used classical equation of motion, . denotes relaxation time and does Fermi momentum distribution with the vanishing chemical potential; . represents a deformed momentum distribution caused by the effect of the electric field. Since the electric current is given by
we obtain the electric conductivity
The relaxation time is given by the inverse of the scattering rate between electrons and positrons,
where we denote longitudinal momentum of each particle ( parallel to ) as , , ,and . We use the notations like , etc. to distinguish the contributions of different Feymann diagrams, which was used in the paper W (). We also use the notations and to distinguish the contributions of the processes and .
Our concern is scattering of the massless fermions occupying only the states in the lowest Landau level. Their wave functions are given by
with and Dirac spinor ,
where we use the gauge potential of the magnetic field . and denote spatial lengths in the direction of and axis, respectively. The states are characterized by two momenta and . Using the formulae, the scattering amplitudes can be calculated in the following way. For example, the scattering amplitude of associated with is obtained by estimating the formula,
with and . Here we also used another notations in eq(17). Using the amplitude , we define the scattering rate ,
The other scattering rates can be obtained in the similar way. Consequently, we find the formulae of which are fairly simplified since the fermions are massless,
where the delta functions represent momentum and energy conservation. All of the formulae involve the integration variables, and .
Now we explain that the electric conductivity behaves as in the limit of . In the integrals of , we rescale all momenta such as e.g. . Then, we find in which dimensionless quantity depends on and through the dimensionless quantity . Using the delta functions we can trivially perform the integration of the variables , and . Among remaining integrations of , and , the integration of is finite at due to the factor for . On the other hand, the integrations of and are finite at or owing to the factors or and or