#### Abstract

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 , the first order contribution to the electron self-mass, 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 . 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.

#### 1. Introduction

Renormalization techniques of perturbation theory are used to calculate temperature dependence of renormalization constants of QED (quantum electrodynamics) at finite temperatures [1–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–8] to calculate the electron mass, wavefunction, and charge of electron as renormalization constants of QED [9–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 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–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 standard big bang model of cosmology [13–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–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.

#### 2. Self-Mass of Electron

The renormalized mass of electrons can be represented as a physical mass of electrons and is defined in a hot and dense medium as where is the electron mass at zero temperature and 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 and 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 in thermal background. For this purpose, we sum over all of the same order diagrams. Renormalization is established by demonstrating the order-by-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 [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 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 and can be neglected in the low-temperature limit, giving

In the high-temperature limit, and are vanishingly small and the total fermion contribution comes from , yielding

The above equations give for the low temperature and 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 high-temperature ranges can be written as showing that the at , for a heating system, and for a cooling system and the difference between two values 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 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.

#### 3. 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 . 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 ’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.

#### 4. Wavefunction Renormalization

The electron wavefunction in QED is related to the self-mass of electron through the Ward identity. The factor 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] giving the low-temperature values as and high-temperature value as

For small values of the electron velocity , 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 in the relevant temperature range. These terms are suppressed at large values of electron energy and they are suppressed by a factor . 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 high-temperature 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.

#### 5. 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].

#### 6. 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 , 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 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 () is the plasma screening frequencies or self-mass of photons that contribute to the QED coupling constant at finite temperatures. For generalized temperatures, the charge renormalization constant can be written as [11]

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 and leaves the perturbation series valid for at least .

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 (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 at becomes negligible as soon as the universe cools down to .

#### 7. 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 two-loop 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 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.

#### 8. Results and Discussion

A quantitative study of the QED renormalization constants, at finite temperatures, shows that all the renormalization constants are finite at . 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 two-loop 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 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 high-temperature contributions are derivable from these general expressions.

We plot thermal contributions to electron mass , electron wavefunction renormalization constant (), and the electron charge 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 level. The distinction in behavior starts near without affecting the renormalizability. Due to the small contribution, the behavior of renormalization constants, for both orders, is comparable at low temperatures, indicated by Figures 2 and 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 , 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.

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.