Renormalization of QED Near Decoupling Temperature

We study the effective parameters of QED near decoupling temperatures and show that the QED perturbative series is convergent, at temperatures below the decoupling temperature.The renormalization constant of QED acquires different values if a system cools down from a hotter system to the electron mass temperature or heats up from a cooler system to the same temperature. At T = m, the first order contribution to the electron self-mass, δm/m is 0.0076 for a heating system and 0.0115 for a cooling system and the difference between two values is equal to 1/3 of the low temperature value and 1/2 of the high temperature value around T∼m. This difference is a measure of hot fermion background at high temperatures. With the increase in release of more fermions at hotter temperatures, the fermion background contribution dominates and weak interactions have to be incorporated to understand the background effects.


Introduction
Renormalization techniques of perturbation theory are used to calculate temperature dependence of renormalization constants of QED (quantum electrodynamics) at finite temperatures [1][2][3][4][5][6][7][8][9][10][11]. The values of electron mass, charge, and wavefunction, at a given temperature, represent the effective parameters of QED at those temperatures. The magnetic moment of electrons, dynamically generated mass of photons, and QED coupling constants are estimated as functions of temperature. However, thermal contributions to electric permittivity, magnetic permeability, and dielectric constant of the medium are derived from the vacuum polarization. Some of the important parameters of QED plasma such as Debye shielding length, plasma frequency, and the phase transitions can be obtained from the properties of the medium itself.
In this paper, we reexamine the analytical results of temperature dependent renormalization constants and prove that QED can safely be renormalized at finite temperatures, using the perturbation theory in a vacuum, below the neutrino decoupling temperature. We use the renormalization scheme of QED in real-time formalism [1][2][3][4][5][6][7][8] to calculate the electron mass, wavefunction, and charge of electron as renormalization constants of QED [9][10][11]. It is now well known that the existing first order calculations of the renormalization constants, in the real-time formalism, give the quadratic dependence of QED parameters on temperature , expressed in units of electron mass . Renormalization constants of QED, using the perturbation theory, give effective parameters of QED in a hot and dense medium and are very useful to understand the physics of the universe. However, the renormalization scheme is fully reliable below the decoupling temperature only. It is explicitly checked that, in the existing scheme of calculations, the theory remains renormalizable at ≤ 4 ∼ 2 MeV. However, the perturbative corrections will exceed the original values of QED parameters at higher temperatures and hard thermal loops have to be dealt with, using already developed methods [12][13][14]. However, below the neutrino decoupling temperature, the real part of the propagators is enough to describe the perturbative behavior of the system and doubling of the field is not required.
We discuss here the distinct behavior of QED below the neutrino decoupling temperatures. At these temperatures, QED coupling starts to play its role in modifying QED parameters for nucleosynthesis. With the help of these effective parameters of QED, the abundance of helium in the early universe can be estimated [12] precisely at a given temperature. The temperature dependent QED corrections to the nucleosynthesis parameters improve the results of the 2 Physics Research International standard big bang model of cosmology [13][14][15][16] and they help to relate the observational data, for example, WMAP (Wilkinson Microwave Anisotropy Probe [15]) with the big bang theory. The same techniques can even be used to calculate perturbative effects in QCD [17] and electroweak processes [18] at low temperatures and estimate their contributions at high temperatures also. We assume that the change in the properties of neutrinos [19][20][21][22][23] does not affect the decoupling temperature significantly.
Without giving the calculational details, we have to briefly overview the existing form of QED renormalization constants in real-time formalism. The Feynman rules of vacuum theory are used with the statistically corrected propagators given as for bosons, and for fermions in a hot medium. Since the temperature corrections appear as additive contributions to the fermion and boson propagators in (1) and (2), temperature dependent terms can be handled independently of vacuum terms at the one-loop level. We restrict ourselves up to the two-loop level to show explicitly that second order thermal corrections are finite and smaller than the first order corrections [9], below the neutrino decoupling.
In the next section, we give the calculations of thermal corrections to the mass renormalization constant / of QED and the physical mass of electrons at finite temperatures. Section 3 is devoted to the calculations of electron wavefunction at finite temperature. Section 4 gives the calculations of the electron charge and the QED coupling constant up to the two-loop level, whereas Section 5, the last section, is devoted to the discussions of the results of all the calculations to explicitly prove the renormalizability of the theory up to the two-loop level, below the decoupling temperature.

Self-Mass of Electron
The renormalized mass of electrons can be represented as a physical mass phys of electrons and is defined in a hot and dense medium as where is the electron mass at zero temperature and ( = 0) represents the radiative corrections from a vacuum and ( ) are the contributions from the thermal background at nonzero temperature . The physical mass of electrons can then be represented as a perturbative series in and can be written as where 1 and 2 are the shifts in electron mass in the first and second order in , respectively. The physical mass is deduced by locating the pole of the fermion propagator ( ( ¡ + )/( 2 − 2 + )) in thermal background. For this purpose, we sum over all of the same order diagrams. Renormalization is established by demonstrating the orderby-order cancellation of singularities. All the finite terms from the same order in are combined together to evaluate the same order contribution to the physical mass given in (4). The physical mass in thermal backgrounds, up to order 2 [7,8], is calculated using the renormalization techniques of QED. Self-mass of electron is expressed as where ( ), ( ), and ( ) are the relevant coefficients that are functions of electron momentum only. We take the inverse of the propagator with momentum and mass terms separated as The temperature dependent radiative corrections to the electron mass up to the first order in are obtained from the temperature dependent propagator as giving where / is the relative shift in electron mass due to finite temperatures which was determined in [1,9] with The convergence of (4) can be ensured for < 2 MeV as \ is always smaller than unity within this limit. This scheme of calculations will not work for higher temperatures and the first order corrections may exceed the original values of QED parameters, after 5 MeV. At the low temperatures, < , the functions ( ), ( ), and ( ) fall off in powers of − in comparison with ( / ) 2 and can be neglected in the low-temperature limit, giving Physics Research International 3 In the high-temperature limit, ( ) and ( ) are vanishingly small and the total fermion contribution comes from ( ) → − 2 /12, yielding The above equations give \ = 7.647 × 10 −3 ( 2 \ 2 ) for the low temperature and \ = 1.147 × 10 −2 2 \ 2 for the high temperature, showing that the rate of change of mass \ is larger at > as compared to < . Subtracting (10) from (11), the change in \ between low-and hightemperature ranges can be written as showing that the Δ( \ ) = 3.8 × 10 −3 at = , \ = 0.0076 for a heating system, and \ = 0.0115 for a cooling system and the difference between two values Δ( \ ) = 0.0038 such that Δ( \ ) is equal to 1/3 of the low-temperature value and 1/2 of the high-temperature value at = . This difference is due to the photon background contributions at low temperatures and additional hot fermionic background at high temperatures. Therefore, the absence of hot fermion background contributes to a 50% decrease in self-mass as compared to cooling universe corresponding value. The high behavior will give 33% more self-mass as compared to the low behavior. Since Δ( \ ) quadratically grows with temperature, the fermion background contribution dominates over the hot boson background after the nucleosynthesis.
The temperature dependence of QED parameters is a little more complicated and significant during nucleosynthesis because of the change in matter composition during nucleosynthesis. Therefore (7) and (8) are required for the ∼ region and help to estimate the change in QED statistical behavior due to the change in composition. This difference becomes more visible when we plot (12) and (13) of self-mass at low temperature and high temperature showing that both values start to give a disconnected region near ∼ , that is, the nucleosynthesis temperature. Figure 1 shows that the slope of both graphs (corresponding to (10) and (11)) never meet at the common point = . The disconnected region around ∼ indicates a change in thermal properties for a heating and a cooling system. In a heating system, fermions start to produce around ∼ , whereas in a cooling system fermions are eliminating around these temperatures. Gap between two curves around ∼ is a measure of background fermion contribution.
That modification in the electron mass behavior in the range ∼ is estimated by (10). It is also clear from (11) that, after around 5 MeV, the temperature dependence correction term ( \ ) approaches unity or even bigger for higher temperatures, even at the one-loop level. Higher order corrections [7,8] will also grow rapidly at high temperatures, giving A change occurs around ∼ and is clearly related to nucleosynthesis where, during the cooling of the universe, right after decoupling, the beta decay processes, involving the electron mass, change the composition of matter and electrons pick up thermal mass from the hot fermion loop. At > 4 MeV, the renormalization scheme of perturbative QED may not be a very good theory as beta decay contributes through weak interactions, which may cause hard thermal loops and the associated singularities in QED. Therefore, all of our discussions in the following sections are referred to below decoupling temperature.

Second Order Contribution
Second order thermal corrections to the electron mass come from the two-loop diagrams. The overlapping diagrams give overlapping hot terms with divergent cold terms and the calculations become really cumbersome. However, in the limiting cases, simpler expressions can be obtained, both for low and for high limits. In the limit < , the second order thermal contribution to electron mass [9, 10] is This equation shows that the low-temperature expression for thermal contribution is very complicated, due to the overlapping hot and cold terms at the two-loop level, as compared to the thermal one-loop contribution which is just 2 / 2 . However, the leading order low-temperature second order contribution is simply whereas the first term indicates the contribution from the disconnected graph which is usually expected from the iteration method. The second term in this expression is clearly dominant for < . In the limit > , the electron mass contribution is a long expression and can be found in [10,11] in detail. Numerical evaluation is not simple and we postpone it for now. However, the leading order contribution at > can be written as The coefficients M's in (16) are complicated functions of electron mass, energy, and velocity of electron giving the leading order contribution as The second order contribution in the above equations is just a leading order contribution to prove that the second order contribution cannot be higher than the first order contribution below the decoupling temperature only. Higher order terms can blow up even at the lower temperatures.

Wavefunction Renormalization
The electron wavefunction in QED is related to the self-mass of electron through the Ward identity. The factor (1 − ) is required for renormalization, because then the propagator can also be renormalized by replacing Thus, for Lorentz invariant self-energy, the wavefunction renormalization constant can equivalently be expressed as The fermion wavefunction renormalization in the finite temperature field theory can be obtained in a similar way as discussed in vacuum theories. However, the Lorentz invariance in the finite temperature theory is imposed by setting = in (5). Thus using (14) and (5), one obtains [9] −1 giving the low-temperature values as and high-temperature value as For small values of the electron velocity V, thermal contributions to the wavefunction renormalization constant can be determined from (21) and (22) as for low temperature, and for high temperature. The finite part of (21) and (22) is equal to / at = 10 in the relevant temperature range. These terms are suppressed at large values of electron energy and they are suppressed by a factor 2 / 2 . However, the calculated value at that temperature is significantly different. The difference in the thermal contribution can easily be found to be about 50% of the low-temperature value and around 33% of the hightemperature value, just as in / . This difference can be mentioned as The finite part of the wavefunction renormalization constant can be obtained by finding a ratio of temperature with the Lorentz energy . The minimum value of this energy is equal to mass. Following (13), the higher order contributions to the wavefunction can then be written as Equation (26) indicates that thermal contributions to the wavefunction are always smaller than / because is always greater than . At = and even at the temperatures higher than nucleosynthesis, convergence of the series can be established as < at those temperatures. So the two interesting physical limits give smaller thermal contribution in electron wavefunction as the relevant temperature limits can be defined as < < and < < , which ensures the renormalizability of QED at comparatively higher temperature as compared to self-mass.

Second Order Thermal Corrections to the Electron Wavefunction
The renormalization of the wavefunction is directly related to the self-mass of electrons and the guaranteed finiteness of electron mass at finite temperatures, below the decoupling temperature, ensures the finiteness of the wavefunction automatically. However, the detailed expression for the wavefunction renormalization constant can be found in [10] and can be given as Similarly, the high-temperature limit for the wavefunction renormalization constant gives in terms of the one-loop integrals 's and 's which can be found in the original papers [9, 10].

Photon Self-Energy and QED Coupling Constant
The self-energy of photons and the electron charge also behave differently for a cooling and a heating system around = . It is well known that the electron charge and the coupling constant do not show significant temperature dependence for < . However, they have significant thermal contributions at high temperatures ( > ). Differences in the behavior of a cooling and a heating QED system start again near ∼ from the lower side and near the decoupling temperature from a high side. This difference in the coupling constant looks more natural due to the beta decay processes during nucleosynthesis.
Calculations of the vacuum polarization tensor show that the real part of the [9,10] longitudinal and transverse components of the polarization tensors can be evaluated, in the limit → 0, as giving the interaction potential, in the rest frame of the charged particle as [8] ( ) where is the renormalized charge. ( ) can be written at low temperature as: The constant in the longitudinal propagator is the plasma screening mass; therefore, the outside factor corresponds to the charge renormalization and in turn to the coupling constant. We may then write the coupling constant at low temperatures as The factor 1.55 × 10 3 ln( / ) is a slowly varying function of temperature and does not give any significant contribution near the decoupling temperature and remains insignificant for a large range of temperature, due to the absence of significant numbers of hot electrons in the background. Therefore, the coupling constant is not modified at < at all. The temperature dependent factor in the longitudinal propagator (4 2 / 2 ) is the plasma screening frequencies Also, the electric permittivity is and the magnetic permeability is In the limit > , the wavefunction renormalization constant can be written as giving the renormalized coupling constant as Equation (36) gives 3 = 1 + 1.2 × 10 3 2 / 2 and leaves the perturbation series valid for at least ≤ 4 .
It is clear from (32)-(36) that the coupling constant is basically changed due to the hot fermion loop contributions. Hot bosons do not change the coupling constant as the vacuum fluctuations occur due to fermion loops, at the first order in . However, the situation is different for higher order contributions.
At > 4 (the decoupling temperature), thermal contributions are significant enough to grow the coupling constant to the level where it can create a problem for the convergence of the perturbative series of QED. In that case, we need to use nonperturbative methods to establish the renormalization of QED at high temperatures. However, thermal corrections to the coupling constant are significant until the temperature is of the order of electron mass. When it is lower than the electron mass and the primordial nucleosynthesis almost stops attains the constant value (1/137). It happens because the constant thermal contribution from fermion background ( ( ) = − 2 /12) at > becomes negligible as soon as the universe cools down to < .

Second Order Correction to the QED Coupling Constant
The first order thermal corrections do not contribute to the coupling constant at < but they do not vanish at the twoloop level because of the overlap of the vacuum term and the thermal term. The low-temperature second order correction to coupling constant can be given as and the high-temperature contributions are given as Since the contribution to the coupling constant is always proportional to 2 / 2 and is sufficiently smaller than , that ensures the renormalizability at the temperatures below the decoupling temperature.
In the next section, we will summarize these results and discuss the behavior of QED in thermal background, below the decoupling temperature.
Physics Research International 7

Results and Discussion
A quantitative study of the QED renormalization constants, at finite temperatures, shows that all the renormalization constants are finite at ≤ 4 . Also the higher order radiative corrections are smaller than the lower order perturbative corrections in real-time formalism. In this range of temperature, the largest thermal contributions come from the mass renormalization constant. The wavefunction renormalization constant is suppressed, especially at low temperatures and relativistic energies. There is no low-temperature contribution to the electron charge and QED coupling constant, as the hot fermion loop contributions are ignorable at low temperatures ( < ). However, a rapid growth in the renormalization constants at large temperatures indicates increasing thermal corrections with the increase of hot electrons in the background. A comparison between the one-loop and twoloop contributions is shown to prove renormalizability of the theory, including thermal corrections. What happens around ≥ is expressed in terms of ( ), ( ), and ( ) functions, at the first loop level. These functions give vanishing contributions at low temperatures as they arise from the integration of the hot fermion propagator. Contributions of ( ) vanish at low and it sums up to (− 2 /12) for large values. However, at the two-loop level, thermal contributions overlap with the vacuum terms and are expressed in a much more complicated overlapping series as well as the same ( ), ( ), and ( ) functions. The problem arises when diverging vacuum terms overlap with the temperature dependent terms. However, they are renormalizable below decoupling temperatures. The low-temperature and hightemperature contributions are derivable from these general expressions.
We plot thermal contributions to electron mass ( / ), electron wavefunction renormalization constant ( −1 2 + (2 / ) ∫( / ) ( )), and the electron charge ( 3 − 1) that can be derived for < and > ranges from the same, (8), (19), and (33), respectively, at the one-loop level. Since we are dealing with the exponential functions in this study, even less than an order of magnitude difference is a safe limit to use low-temperature and high-temperature limits. Just to prove the renormalizability, we present the plots of low-temperature and high-temperature terms of renormalization constants, as they are derived from general expressions. We give a comparison between a heated and a cooled QED system, around the common temperature of the order of electron mass (see Figure 1). The difference between thermal background contributions of a heating and a cooling system changes due to the difference of background during the heating and cooling process. A cooling system such as the early universe starts off with more fermions and keeps losing them during the cooling process, whereas a heating system will start with a minimum number of fermions and they will be created during the heating process. We just quantitatively discuss the self-mass contribution as it changes the physically measurable mass, a very important parameter of the theory. Figure 1 indicates the difference between the self-masses of electron due to the low-and high-temperature first order corrections. It is found that the difference between two values is equal to 1/3 of the low-temperature value and 1/2 of the high-temperature value at = . Similarly the electron wavefunction and charge for a system that approaches from lower temperature to = will be different from the one reaching to the same temperature by cooling of a hotter system. It is understandable as the photons do not couple with each other directly. They only couple to the charged fermions and the presence of charge fermion is required to affect the QED coupling. The fermion distribution function (2) actually brings in the fermion loop contributions at high when more fermions are generated during the nucleosynthesis and are not ignorable any more in the medium.
We compare thermal contribution to the first order [9] and the second order [10] corrections to the self-mass of electron in Figure 2. The solid line corresponds to the first order corrections and the broken lines correspond to the second order corrections. Figure 3 gives a similar graph for the renormalized coupling constant of QED in a hot medium. A plot of these renormalization constants shows that the temperature corrections do not affect the renormalizability of the theory for sufficiently smaller than up to the two-loop 8 Physics Research International  3. This difference is significant at high temperatures. The plot of the QED coupling constant as a function of temperature shows the very small effect at low temperatures, at both loops. It only becomes significant for > . In the approximations, used in this paper, one-loop thermal contribution is zero at the first loop level at low temperature. However, at high temperatures, the coupling constant grows quadratically with temperature, expressed in units of electron mass. With a large coupling constant, the renormalizability of the theory cannot be guaranteed and nonperturbative methods have to be used to treat the hard thermal loops. However it is not needed under the decoupling temperatures at all. Table 1 shows the leading order thermal contributions to the numerical values of the renormalization constants for electron mass and the QED coupling constant near the decoupling temperatures. It also shows that the thermal contribution of the renormalization constants at temperatures around the decoupling temperature becomes significant. It is a table of the leading order contributions to prove the renormalizability. So we consider ≤ 4 , just to get an approximate thermal contribution to electron wavefunction. However, it is explicitly shown in the above figures and Table 1 that the renormalization scheme of QED works below the neutrino decoupling temperature. Above the decoupling temperature, a large number of hot electrons in the background lead to the failure of the QED renormalization scheme at larger temperatures and the electroweak theory has to be incorporated.
The cooling universe of the standard big bang model behaves differently after the neutrino decoupling. Nucleosynthesis starts right [17,18] after the neutrino decoupling and the helium synthesis takes place when the temperature of the universe is cooled down to the temperature of electron mass. This is actually a temperature, where the finite temperature corrections to QED parameters are significant but complicated enough to evaluate it numerically. However, the temperature dependent QED parameters are needed to describe the observations of WMAP data. After the nucleosynthesis is complete, the temperature dependence is back to the quadratic dependence on temperature, though it is not exactly the same as at the low temperature. Fermion background contribution can easily be seen.