1. Introduction
As quoted long ago by Pisarski [
1], “It’s really surprising how difficult it is to calculate damping rates in hot gauge field theories”.
Even though a so-called Resummation Program of Hard Thermal Loops (that are leading order thermal fluctuations at high temperature, hereafter denoted by
) has been devised to obtain gauge-invariant complete results to damping rates (and has succeeded to some extent [
2,
3,
4]), serious difficulties have constantly bounced back and forth within the infrared sector of hot gauge field theories.
Literature testifies of two major obstructions to the Resummation Program that have become textbook material [
5]. The first one was discovered about the same time as the Resummation Program itself [
6], while using the latter in the situation to be described in the next section. The second obstruction was discovered a few years later, while using the Resummation Program to evaluate the soft photon emission rate out of a Quark–Gluon Plasma in thermal equilibrium [
7]. In both situations, the hot gauge theories infrared sector was recognized to be at the origin of two singular results.
In the latter case, the singular result reads (
),
where
is the space-like part of the fermionic spectral densities, and
with
, [
5]. The explicit Lorentz invariance breaking appears in the statistical factor
due to its dependence on
T, which associates with a preferred rest frame. In effect, in the context of quantum field theories, temperature at thermodynamical equilibrium can be implemented in a covariant way by means of a Lorentz 4-vector,
, of squared norm
. (We work in units where
, Boltzmann’s constant, is set equal to unity.) Only in the reference frame, where the plasma is at rest and in thermal equlibrium such that
(
u the 4-velocity of the plasma) does the invariant
allow for the interpretation of a temperature. The singularity in Equation (
1) is due to the following double entwined angular integrals. With
, the two light-like 4-vectors
and
,
which, to our knowledge, computers cannot evaluate. An estimated singular behaviour of
was accordingly retained, giving rise to the singular result above, whereas a full, exact and cross-checked calculation of
displays a series of mass singularities of strengths
and
, which exactly cancel out among themselves, leaving a regular result [
8,
9,
10].
In the former case, the current revisitation of the problem will conclude that its singular result has been induly derived as well.
As will be stated in the conclusion of this analysis, it turns out that the issue of infrared singularities has long masked a more fundamental difficulty, which appears to be inherent to the perturbative approach and has been loosely identified within the context of
-algebraic analyses [
11]. Elaborating on this fundamental difficulty is the objective of the current paper, organized as follows.
Section 2 is a quick reminder of the historical fast fermion damping rate calculation in a thermalized quantum electrodynamics (QED) or quantum chromodynamics (QCD) plasmas. The notations to be used throughout the article, as well as the process basic diagram are given.
Section 3 accounts for the photon polarization tensor at two-loop order, or equivalently, fast muon damping rate at three loops. The leading order contributions are retained in the limit of
where
is a component of the internal photonic line,
K. This is essentially technical, with some overall detailed balance examples of infrared (mass) singularity compensations displayed. The important point is comprised within the simple Equation (
55), which expresses the leading form of the photon polarization tensor at two-loop order.
Section 4 announces the point of the current paper, the peculiarity of the thermal case, by contrasting it to the customary
situation, as an independent resummation of leading fluctuations occurring at momentum scale
is shown to be necessary.
While it is believed that the ultrasoft scale is terminal, receiving contributions from higher number of loop diagrams, a three-loop evaluation of the photon polarization tensor appears to display an interplay of thermal and vacuum fluctuations leading to contributions of dominant order and calling for another independent resummation along the internal photonic line. Revealed by a high enough number of loop calculation, these leading mixed fluctuations of the photonic internal line take place at momentum scale
, are gauge invariant and contribute to the fast muon damping rate being on the same footing as the two first ones, at soft (
) and ultrasoft (
) momentum scales. This is what is dealt with in
Section 5.
Since only the contributions attached to the longitudinal degree of freedom of the internal photonic line were considered in the previous sections, the contributions of transverse degrees are examined in
Section 6 so as to make sure that there is not a full annihilation of the contributions due to the longitudinal leading order mixed fluctuations.
Section 7 is a discussion of our results in the light of previous
-algebraic analysis. The relevance of a perturbative approach bare or effective is questioned while possible non-perturbative alternative approaches are mentioned. Bare perturbation theory admits the same diagrammatic expansions as customary
quantum field theories, but now subject to the
-free field propagators and vertices [
5].
2. A Reminder of the Fast Muon Damping Rate Issue
The energy loss of a heavy fermion, a muon for example, propagating through a QED plasma is considered. It is assumed that the heavy muon has mass M much larger than the temperature T, and that, travelling through the plasma, it retains momenta and thus is not thermalised. In the plasma, electrons (positrons) are assumed to be massless, and, for the muon, energy loss is induced by Compton scattering off photons and Coulomb scattering off electrons and positrons of the plasma.
In the framework of real-time thermal perturbation theory of QED, at temperature
T, the damping rate is related to the fermionic self energy
by the relation [
5,
12],
The notation
will be used throughout. In order to get rid of inessential complications, it has been customary and convenient to assume the following ordering relations,
In view of
, in effect, Equation (
3) allows one to take the high velocity limit
in coincidence with the fermion mass shell limit, as in the massless case. We assume zero chemical potential and a high temperature, which is much bigger than the electron, positron or quark masses. One gets
where
, as depicted on
Figure 1. In Equation (
4), the
distribution is the muon
cut propagator taken at
, since, as stated above, this real particle is not at thermal equilibrium with the QED plasma. The function
instead is the thermalized cut photon propagator in Coulomb gauge. Whenever the
is used, quantities are denoted with a symbol
, and one has
The operators
are the
- longitudinal and transverse projectors [
5], and
the transverse and longitudinal spectral densities of the effective propagator. Eventually,
is the ordinary Bose–Einstein statistical factor, but in the
- real-time formalism being used [
5], without absolute value prescription, for example,
.
Summing Equation (
4) over the angles by taking advantage of the distribution
one gets, with
and
,
Next, since momentum
K is soft, on the order of
, one expands the statistical factor and writes,
The most infrared dangerous piece of the Bose–Einstein factor expansion is the term
which is hereafter retained. Appealing to the energy sum rule,
where
m is defined below, Equation (
14), one obtains, out of the transverse part, the well known mass shell singular result,
where the logarithmic diverging behaviour by the fermion mass shell is manifest. The upper bound,
, is some conveniently choosen intermediate scale between the soft,
, and hard
scales, usually the scale of order
. The
contribution, instead, is regular, the relevant sum rule being,
All of the proposed solutions rely on one or another use of the Resummation Program, that is a one-loop resummation of leading order thermal fluctuations. In Refs. [
13,
14], for example, a resummation of multi-loop exchanges is performed
à la Bloch–Nordsieck, within the one-loop effective theory provided by the Resummation Program. This approach, adequate in view of the ordering relation (
3), yields a finite, regular damping rate, both in QED and QCD. Calculations, though, clearly differentiate the abelian QED, from the non-abelian case of QCD, whereas both situations are formally identical at this order. At Hard Thermal Loop order in effect, the QCD effective Lagrangian density ‘abelianises’ in that it only differs from the QED effective lagrangian density, by the relevant Lie-algebra valuation of the gauge fields,
in QED versus
in QCD, generators
being taken in the appropriate representation [
15].
The current article aims at taking advantage of this long standing issue in order to point out a structural peculiarity of perturbative treatments when applied to the high temperature quantum field theory context. This is done by means of an up to four-loop analysis of the fast muon damping rate, where the number of loops refers not to the effective perturbation theory, i.e., to the
, but to bare perturbation theory and thus to the number of loops appearing in the basic diagram of
Figure 1.
3. Damping Rate at 3 Loop Order
Loop orders are hereafter referred to bare real-time perturbation theory endowed with Advanced and Retarded free field propagators. For example, four-loop order refers to that of diagram of
Figure 1, as the photon polarization tensor
is itself evaluated at three-loop order, as discussed in
Section 5.
Within the Resummation Program, one has
with the one-loop expressions of
and the thermal masses are [
16],
for QCD with
colors and
flavors of quark fields.
As displayed in Equation (
6), transverse and longitudinal degrees of freedom add up separately to the total damping rate (
6), and the latter will be considered now. One has [
5],
where
and
correspond to the two-loop pieces of
Figure 2 and
Figure 3, respectively (Retarded/Advanced indices,
, omitted for short).
One finds (with
),
and
where
is a light-like four-vector, and where as in Equation (
5), the usual Fermi–Dirac distribution
is deprived of any absolute value prescription. In (15) and (16), the mass squared term
stands for either of the two QED and QCD possibilities of (13), and the shorthand notation of
is used, as well as
for
. The
-fermionic self-energy components are
with the results,
and
being the Legendre functions,
In this way, the
prescriptions of Equation (
19) are encoded in the logarithmic determinations of the Legendre functions. For completeness, the
vertex reads as
It is convenient to choose the
z-axis along the vector
, so that
and
, with
. Out of Equations (
16) and (
17), several integrals come about,
and the leading order contributions to any of the seven above integrals will be retained in the limit of
. Note that, at this stage, the definite order of magnitude of
is not yet (and needs not being) determined. As it will turn out shortly, the specific order of magnitude of
is selected by the two-loop order fluctuations themselves.
The first three integrals correspond to the
self energy insertion depicted in
Figure 2 and yield
With , and accordingly , the parity properties of the Retarded determination of are met as they should.
An apparet leading behaviour shows up with the factor , but the whole contribution is ruled out by the odd -parity of the Legendre function .
For
, one finds
with
for the real part, and
for the imaginary part. Intermediate mass singularities develop at
, and cancel each other out in the end as further examples will display below.
For , contributions are found to cancel out, due to x- and -odd parities.
For
, the calculations of [
8,
9,
10] can be used (in particular Equation (3.15) in [
9]) to yield
On the basis of the leading order expansion,
it is immediate to check that the contributions attached to Equation (
34) are cancelled by the odd parity character of the integration on
.
Then, following [
10], to take the discontinuity of
, one has first to dispose of the result,
so as to proceed with the evaluation of
,
Considering the second term of the second line in curly brackets of (36), the sum on
yields
where
stands for
when letting
. At
, one obtains
Integrating on
x, the mass singularities, at
, are, therefore,
plus
The sum of the two previous expressions cancels out the
singular terms which are plaguing the real part of (34), while the imaginary part is finite and reads
With the overall multiplicative factor of , this produces an irrelevant imaginary part (even in ).
The regular (real) contributions that are left read
and add up to a null result,
For the first term in the second line of (36), one has at leading order
, after having integrated over
,
Integration on
x gives no real part,
and a non-zero imaginary one,
On the whole, at leading order, the evaluation of (36) yields no contribution to
other than Equation (
47) whose imaginary character comes from the only condition
and therefore is kinematically compatible with the overall constraint of Equation (
6), which is
.
Comments are in order. In Equation (
27), the fourth contribution is by a factor
- subleading and can be ignored. In the third line, a contribution of order
seems to dominate all the other contributions by a factor of
. One has in effect
where the QED value of
has been used for the sake of illustration, and where
denotes the function
with
such as defined after (28). The function
is infrared regular, but ultraviolet singular. At large
, it exhibits the ultraviolet logarithmic divergence. It arises from the
-part of the self-energy diagram of
Figure 2 and due to the factor 1 in
. One has, indeed,
One-loop renormalized at
, this ultraviolet logarithmic divergence is thus replaced by a (renormalization scheme dependent) constant
, and the corresponding contribution to
therefore is
. This also is the case of the non-zero
T-part of
. One finds [
17],
where the approximate equality refers to the substitution of
for
and to the neglect of
in the logarithm, and where
is the exponential integral function. Note that this is consistent in view of Equation (
6), which indicates that
is on the same order of magnitude as
k, and thus
. By preserving that term, however, the same result is reached.
Now, since
, with
the Euler’s constant, the two first terms of Equation (
51) yield
and are of order
, as are also the two remaining terms in the second line of Equation (
51). As displayed in
Figure 4, a numerical evaluation of the full, non-approximate
function of (48) confirms this conclusion, which can be summarised as follows:
To summarise, at two-loop order of bare perturbation theory and in the small
-limit, inspection of expressions (
22)–(
28) reveals that there are no leading order contribution to
other than those of form
, where
denotes some dimensionless function of
and
. In replacing
of QED by the one of QCD, the same conclusion applies. Two remarks are in order:
(i) The notation used above,
, instead of a would-be more natural notation like
for example, is to mean that as well known [
18,
19], not all of the leading order terms in the limit of
are accounted for by
. This is explained below, after Equation (
58).
(ii) For reasons that are made explicit in
Section 6, Equation (
104), in the present as well as in the following section, the contributions of transverse modes are not considered.
4. Perturbative Peculiarity of the Thermal Case
The diagrams of
Figure 2 and
Figure 3 involve the gauge invariant
self energy (
16) and vertex (
17). By changing the gauge, only the (external)
K-line propagator is changed, which doesn’t impact the damping rate
[
20]. One is therefore tempted to conclude that by construction these two contributions to the photonic polarisation tensor,
and
, are separately gauge invariant. This, however, needs to be checked. Calculations are long but straightforward and one finds effectively,
Moreover. these relations hold at leading -order, exemplifying that the principle of local gauge invariance is so deeply rooted in the theories that it is verified at any consistently determined momentum scale.
Taking advantage of the last results, Equation (12) could now be improved by writing
and used to calculate the longitudinal contribution to the moving fermion damping rate at this order of approximation. According to the previous section, one has (QED case chosen for the sake of illustration)
where
collects the various leading order contributions of (
22)–(
28) in the small
-limit, into some complex-valued function of
and
k. Like the Legendre function
of the one-loop order, (
13) and (
19),
is assumed to preserve the
- leading order of magnitude of two-loop contributions to
. Still, the rather
academic content of this assumption will be discussed in
Section 7. Re-writing Equation (
54),
it is made clear that, like in the
usual perturbative context, the two-loop contribution
represents but an order
-correction to
, the Hard Thermal Loop counterpart, and doesn’t bring the fermion damping rate anything else than a correction on the same
-order of magnitude, that is, using Equation (
56) instead of
one would get basically,
On the other hand, considering
in Equation (
55), it is clear that, over momenta of order of
, one has
, so that any number of
-insertions along the
K-line will contribute a similar amount to the longitudinal degree of the internal
K-line and must be resummed. For momenta of order
exclusively, this separated resummation leads to a spectral density of form,
and, as displayed in next
Section 5, the matter with Equation (
58), and further with Equation (
84), is that it contributes the fermionic damping rate
on the same order of magnitude as the Hard Thermal Loop spectral density itself,
, that is on the leading order of
.
This situation has no equivalent in the perturbative context of
theories and appears to be peculiar to the non-zero
T-case. In QCD, the necessity of resumming ‘ultra-soft’ order
-fluctuations into another effective perturbation theory was recognized years ago on the basis of an equivalent semi-classical approach [
19] to high-
T quantum field theories, and the new effective theory was shown to enjoy the same symmetry properties as those of the
-resummation program:
(1) For external momenta , ultra-soft amplitudes are as large as the corresponding and tree-level ones, (2) they are gauge fixing independent (covariant and Coulomb-like gauges …), and (3) ultrasoft amplitudes satisfy simple Ward identities. (4) Unlike the soft amplitudes controlled by the resummation program, though, ultra-soft amplitudes receive contributions from an infinite series of multi-loop diagrams (‘ladder diagrams’, to be specific).
It is therefore natural to speculate that a higher number of loop contributions could generate similar conclusions, with gauge invariant fluctuations at smaller momentum scales (orders
,
, etc.), generating for the cut photonic propagator
of Equation (
4) contributions all on the same order of magnitude as the tree-level ones.
The entire perturbative structure would therefore have to be considered along with layers of momentum scales, each of them controlled by a specific effective theory. Point (4), however, seems to cast some doubts on this expectation. In effect, the ultra-soft effective perturbation theory is taken to be terminal, as the ultra-soft denomination itself is intended to mean. In the next section, though, it is ventured that this could not be true, due to a specific interplay of T = 0-renormalized with Hard Thermal Loop vertices.
Though contributing equally to the net fermion damping rate (see Equation (
6)), transverse contributions have not been considered for reasons to be given explicitly in
Section 6.
5. Damping Rate at Four-Loop Order
In this section, three-loop contributions to
are considered such as pictured in
Figure 5, contributing
at four-loop order. One can start with two self energy insertions
and
, which at one-loop order and in real-time formalisms come about additive in their
and
parts,
and
Equations (
18) and (
19). Various terms contribute to
, whose leading terms in the small
-limit admit the generic form,
with some dimensionless associated functions
, possibly one- and/or two-loop renormalized at
. Clearly, in line with the expectations of the semi-classical analyses [
19], these contributions are to be considered as part of the ultra-soft amplitudes for
and do not define any new momentum scale of invariant fluctuations. Examples are as follows.
As bare power counting suggests, one may replace one of the two self energy insertions by its
T = 0-renormalized part only,
. At
, though, the fermionic self energy is not gauge invariant. In the Coulomb gauge being chosen, one has [
21],
so that one could expect contributions to
, on the order of
, to be a good candidate for a new relevant momentum scale at order
. Now, besides a lack of gauge invariance, explicit calculations of the various terms stick to the behaviour (
59) and hence, at best, would only be a part of the ultra-soft momentum scale fluctuations.
The same applies as well to two
self energy insertions
and
. For example, among other contributions, one has with the first term of Equation (
60) (in the small
limit always),
Explicit calculations give
for the leading term
in the second line of Equation (
61), while, for the second (sub-leading term)
, one finds
with
that is, besides an irrelevant
-even imaginary part, an order
subleading real contribution. The last term in the second line of Equation (
61),
, yields another contribution corresponding to a new function,
with correct parities but an order
sub-leading as well.
All in all, inspection of the various terms shows that the contributions to
due to the fluctuations depicted in
Figure 5 can be written in the form (
59) with related
functions
for the leading order ones. More precisely, within the topology of
Figure 5, contributions of order
are effectively met, but, compared to
-ones, they are sub-leading by a factor of
. Accordingly, the contributions attached to the topology of
Figure 4 (some of them retaining further dependences on the Coulomb gauge, what is more) could at best complete the contributions attached to fluctuations at ultra-soft momentum scale, but can in no way define any new gauge invariant fluctuation of the internal photonic line at a softer momentum scale, such as
.
It is with the diagram of
Figure 6 that things start to appear to be different. This is because the zero temperature one-loop renormalized vertex reads (QED)
with
,
representing the electric and magnetic form factors. Using this form of the vertex will be discussed below. The trace over Dirac indices, which are involved in this diagram, are
as the form factor
will not contribute at this number of loops because of an odd number of Dirac matrices. Evaluating Equation (
66) at
,
in order to get
, the result is
Integrated over
Q, Equation (
67) exhibits mass singularities whose cancellations can be seen as follows.
The first line of Equation (
67) is easily evaluated to give
Within the integration range of (5), a non-trivial imaginary part develops, which is essential to a non vanishing spectral density
,
This evaluation of the imaginary part is exactly what is obtained by putting the two fermionic lines on mass shell in agreement with the use of Equation (
65) for the
T = 0 part of the one-loop renormalized vertex correction.
Turning to the second line of Equation (
67), the first term in curly brackets (not proportional to
) contributes a term,
with an associated imaginary part of
For the second term in curly brackets in the second line of Equation (
67) (proportional to
), two non-vanishing
-discontinuities can be taken.
Taking the discontinuity of
, the first one gets
where as before,
, regularises mass singularities at
.
Then, taking the discontinuity of
(i.e.,
), one obtains
with
the function,
This function is a regular function of
at
, which allows for an isolation of singular behaviours in Equation (
73). At
, one obtains
On the other hand, at
, one has
Now, the first line of Equation (
75) and the first line of Equation (
76) compensate exactly the first line of Equation (
72), while the second lines of Equations (
75) and (
76) cancel each other out. In the end, the second line of Equation (
72) and the third line of Equation (
75) cancel each other out exactly.
The regular contribution of the second group of terms in the second line line of Equation (
67) thus is
and, due to
-improper parities, inducing cancellations in the integration over
, gets eventually reduced to
For the regular part of
, however, only those terms for which the
Q and
internal lines can be considered on mass shell (i.e.,
), or approximately (
), should be retained so as to be consistent with the use of Equation (
65) in lieu of the most general form of the
e.m. vertex (as will be argued later, 3rd and 6th items below, this limitation does not hamper the point of the current paper).
All in all, the contributions (
68) and (
70) only are retained for
and stand for a part of the leading order contributions generated by the quantum fluctuation depicted in
Figure 5 (remark (ii) below):
with the function
At spacelike
K-momenta, Equation (
80) develops a non-zero imaginary part that is essential to a non vanishing spectral density
.
Due to the form factor
a non-trivial effect comes into play between the invariant
T = 0-renormalized vertex (
65) and
- self energy
. When renormalized, in effect,
is gauge-invariant to all orders in perturbation theory [
22]. Since
is also gauge-invariant, the fluctuation of
Figure 5 brings to the cut photon propagator,
, a correction which is gauge-invariant by construction; or, more precisely, which doesn’t retain more gauge dependence than the original one-loop
expression
(as noticed before Equation (
53), this remaining gauge-dependence is not problematic). This gauge invariance, by construction, can be checked also by explicit calculation and is here a mere consequence of Equation (
53).
So far, as displayed for example in the ordering relation (
3), fermionic bare/renormalized masses have been ignored in front of the energy scales involved in the case considered. In the literature though, the one-loop form factor
is always given at non zero fermionic masses and reads [
22],
where
, and where
is a small photon mass introduced so as to regularise intermediate infrared singular behaviours (mass singularities indeed). For reasons to be explained shortly, in line with the use of Equation (
65), this logarithmic term can be consistently ignored. However, the
-dependence of
remains fairly involved. A small
limit of
exists, and when such a regime becomes relevant to the high temperature context of either QED or QCD, Equation (
81) simplifies to (QED case and ignoring the logarithmic term) [
22],
This form will be used to illustrate the main point of the current paper. Note that, in the case of QED, at
MeV and
MeV, one obtains, at
, a value of
,
, which renders admissible the use of Equation (
82) for
K-momenta on the order of
. The following remarks are in order:
(i) Using Equation (
65) for the
e.m. vertex amounts to reduce
in
Figure 6 to the contribution where the internal lines
Q and
are on mass shell, as illustrated in
Figure 7, whereas the full imaginary part of
entails also the cut diagram of
Figure 8.
Both contributions involve the mass singularity of Equation (
81), that is, the logarithmic term
, in such a way that these intermediate mass singularities cancel out in the sum, in agreement with standard
situation (theorem
[
23]). As intermediate terms, which are absent of the final result for
, the logarithmic term
can therefore be dropped out of Equation (
81) without prejudice. As has been known for a long time, it is only when one is interested in a physical process such as the one depicted in
Figure 8 that the intermediate logarithmic singularity of
is retained in a first step, and then consistently replaced by a regular term of
, where
stands for a parameter relevant to the case considered, such as an energy resolution,
, in a measure of the real (cut) photon of
Figure 8 [
24].
(ii) In the literature unfortunately, the general zero-temperature expression of
does not seem to exist in closed form [
25], and Equation (
65) is used instead. Accordingly, the polarisation tensor
of Equation (
79) is not the full one, but a part of it.
The situation can therefore be summarised as follows:
At soft momenta,
, the cut photon propagator
of Equation (
4) receives one-loop
contributions,
on the same order of magnitude as the cut bare one, evaluated at the same soft scale.
At ultra-soft momenta,
, the same applies also to
, or equivalently to the spectral density (
58) in the case of our illustration. To be complete, however, the ultra-soft scale would require the consideration of an infinite number of loop diagrams. In this diagrammatic analysis, this is displayed by Equation (
59) which corresponds to the three-loop topology of
Figure 4, while many more diagrams are known to contribute at this order as well, with contributions of form,
Such a case is not easily dealt with diagrammatically, but, hopefully, the semi-classical approaches [
19] may offer a sufficiently controlled alternative.
This is no longer possible at the hyper-soft momentum scale
, where the topology of
Figure 6 is able to generate a
, equivalently a
spectral density,
which, over this hyper-soft momentum scale, again is on the same order of magnitude as is the bare one evaluated on the same hyper-soft momentum scale. In deriving Equation (
84), Equation (
80) is used to identify the constant
as
and
stands for the expression,
where it is reminded that
parametrizes the ratio
, according to the relation
. A plot of
is found in
Figure 9.
In this last case, incompleteness comes from the use of Equation (
65) instead of the full
, as stated in remark (i) above. Unless explicitly demonstrated, it is not excluded either that higher loop diagrams contribute at the same hyper-soft scale. Completeness in this diagrammatic analysis is again extremely difficult to control [
26].
The point is therefore that the so far controlled situations of soft and ulta-soft scales doesn’t seem to exhaust all of the invariant photonic fluctuations contributing the leading order damping rate calculation, and that the contribution completeness to of the hyper-soft scale is not under full control; not to speak of other even softer invariant fluctuations, such as those arising at momentum scale , not dealt with in the present paper.
Now, there is more to it. Due to the final
k-momentum integration process, normally taken to be the range
, [
5,
13,
14], there remains the issue of partitioning the integration range into layers of every decreasing momenta, which contribute to the damping rate
by comparable amounts.
Each invariant fluctuation, occurring at momentum scale
, is assigned an effective range (a layer thickness, say) of,
in agreement with a common, still somewhat arbitrary usage. We know of no theoretical basis which could help with defining the layers in Equation (
87) with accuracy [
5]. As sketched in
Appendix A, a dramatic consequence of this is that even if, at any momentum scale
, each of the contributing invariant fluctuations yielded a perfectly controlled, complete
pieces (equivalently, a related
spectral densities), their leading order respective contributions to the damping rate,
say, would read typically as,
where the flexibility in the way Equation (
87) can be decided plagues the various values the constants
take.
This situation looks suggestive of two major points. First, several gauge-invariant quantum fluctuations of the internal photonic line, taking place at different momentum scales, may contribute at zeroth order to the fermionic damping rate, keeping aside the issue of sub-leading order corrections. Second, completeness of the leading order fermionic damping rate is out of reach in a (perturbative) diagrammatic approach.
A by-product of this result is that the famous infrared divergence, plaguing the rapid fermion damping rate by its mass shell, (
9), is an artefact of extending the
k-integration naively down to zero,
Namely, because of the existence of softer invariant fluctuations contributing at the same order of magnitude (see
Appendix A), the
effective perturbation theory cannot be extended down to
, but only to some lower limit, such as
, to comply with Equation (
87). Again, no real infrared singularity problem, but definitely something else to be discussed in
Section 7.
6. Transverse Degrees
So far, contributions from transverse degrees have been ignored. As displayed in Equation (
7) though, the contributions attached to
, [
5], should be considered as well, in the same limit of
and same topology (
Figure 6). Several cases are possible. The transverse contributions are ultra-soft ones, scaling like in Equation (
55) and nothing is changed of the previous sections, but for the fact that new contributions to the ultra-soft scale are identified. The transverse contributions are sub-leading with respect to the longitudinal ones (
79), so that they can be ignored. Transverse contributions annihilate the longitudinal ones and the matter of previous
Section 5 fades away. This is accordingly the point to be analysed.
As calculations display, however, supplants by a large factor of , doesn’t call for resummation and does not annihilate longitudinal contributions to the rapid fermion damping rate . This is seen as follows.
Getting back to the diagram of
Figure 6, two Traces are to be considered, which are (see Equation (
66)),
Calculations are a bit lengthy but straightforward. Up to a sub-leading contribution of,
which can obviously be discarded, the first trace of Equation (
90) yields
The second trace of Equation (
90) yields
where again,
regulates the mass singularities which appear at the boundaries of the integration range
. As usual, these singularities are exactly compensated for by those coming from taking the
-discontinuity of the
-self energy (
19),
In effect, integrating out the angular dependences, one obtains
The first term in the second line of Equation (
95) is made out of three contributions, attached to
,
and
, and only the two first ones contribute thanks to the according
-parities. The term
gives a real valued contribution,
The term
is responsible for a mass singularity,
and a regular null expression,
Note that, for the second term of Equation (
93), the piece coming from the constant term of
(i.e., the term
) is exactly cancelled by Equation (
97). Since
, the mass singularities attached to the
already cancelled out as just stated, it remains to prove that the mass singularities of the second term of Equation (
93) are compensated by those of Equation (
95). Now, out of the squared bracket in the second line of Equation (
95), one has the crossed term
and the following integration must be performed,
By the boundaries at
, it is then immediate to check that the mass singularities of Equation (
100) read as
and exactly cancel out the remaining mass singularities of Equation (
93).
What remains are contributions of types,
which, for the second type, are similar to those of the longitudinal case (
80) up to different overall multiplicative constants. The ensuing transverse spectral density would accordingly read as
with, like
,
developing an imaginary part at space-like momenta
K. Now, whatever
, the situation here differs from that of
in Equation (
84) because the real-valued contribution (
96) outweighs any other contribution to
by a large factor of
. This leads to a serious suppression of the transverse spectral density in comparison to the longitudinal one.
Over momenta of order
in effect, orders of magnitude are found to be in the relation,
Contrary to longitudinal contributions, the three-loop transverse contributions of
Figure 6 require no further resummation and preserve the peculiarity of the longitudinal ones, as it was necessary to check. Since Equation (
90) only differs from the two-loop case in the overall multiplicative form factor
, the same relation as (
104) applies to the case of
Figure 2 as well, with a milder depletion factor since in this case, one has
. This is why the transverse contributions were not considered in
Section 3 and
Section 4.
7. Discussion
In real-time formalisms, it has long been suspected, though poorly examined, that higher number of loop calculations could reveal some peculiarities. This was mainly expected of the so-called singularity mixing problem, that is, of some products of
ultraviolet and
-infrared singularities (which are leading over the
ones in the large-
T limit [
27]).
However intricate such a singularity mixing may be, there seems to be nothing really fundamental about it. Decades of efforts devoted to the matter in a large enough variety of cases have always proven that, when properly dealt with, all singularities cancel out in the end. As compared to the zero-temperature case, where theorems can be invoked, at non-zero temperature, patterns of singularity cancellations are only rendered more difficult to encompass into some -like general statements.
At the origin of this fate, two main reasons can be invoked. First, this is because of the explicit breaking of Lorentz invariance brought about by the thermodynamical equilibrium. Singularity cancellation patterns are not so constrained as they are in the
Lorentz symmetric case. In a second place, as is done throughout the current paper, one keeps analyzing
situations in terms of Feynman diagrams while, ...in contrast to vacuum observables, thermal observables are not themselves primitive blocks of Feynman diagrams [
11]. According to the observables under consideration, case-by-case treatments are therefore required, which render difficult attempts at generalizations (note that a somewhat related facet of this state of affairs is the one put forth in [
28]).
Now, as advertised above and supported by a large number of calculations, singularity cancellations should rather be regarded a technical matter rather than a real issue. This is why the current paper, already technical enough, doesn’t pay more attention to the (ultraviolet diverging) constants of Equation (
85) or Equation (
A5), which a relevant and canonical renormalization procedure will turn into finite numbers.
Rather, it is here posited that the singularity problem (infrared sector in particular) of non-zero temperature quantum field theories may not be that main issue we’ve long thought it is [
5]; and that it has long masked a more fundamental difficulty, as soon as one tries to enforce a perturbative treatment of these theories. The alluded more fundamental difficulty surfaces already at the level of the historical infrared singularity problem of Equation (
9) which, as it should be clear now, is artificial and arises from a too liberal continuation of the
k-integration in Equation (
9) down to zero within the
-resummation program.
The fundamental difficulty could be, instead, that any bare and/or effective perturbative treatment be simply unable to reach a complete answer. The example considered in this paper may be regarded as rather suggestive. The damping rate of a fast moving muon is calculated at leading order, as it travels through a QED (or QCD) plasma in thermal equilibrium. Not only the soft, order
, but ultra-soft, order
, and hyper-soft, order
gauge-invariant quantum fluctuations of the emitted and re-absorbed photon contribute the zeroth order muon damping rate, each. In this latter
instance, one observes that the gauge invariant character of the associated fluctuations results from an interplay of a
-self energy insertion with, at
, the gauge invariance of the one-loop renormalized vertex; that is, results from a peculiarity of a calculation performed at a high enough number of loops (of course, such a mixing of
with
parts wouldn’t come about in a Matsubara imaginary time formalism, but while this formalism is more fitted to the calculations of static quantities, it is not deprived of analogous difficulties either [
29]).
Now, aside from another quite obvious case at momentum scale , not addressed in the present paper, possible softer and softer gauge invariant fluctuations could contribute also on the same order. More precisely, on the basis of diagrammatics, it may be extremely difficult, if even doable, to control the possible number of similar gauge invariant fluctuations. By gauge invariant, what is understood is that no gauge dependence other than that of the original bare photonic line is retained within the higher number of loop fluctuations of this line.
As a result, even the zeroth order estimate of a physical quantity could be out of reach of Perturbation Theory bare and/or effective, not even speaking of further sub-leading corrections. More than the issue of singularity cancellations, which is of a technical nature, this difficulty could constitute a serious obstruction of principle to the application of Perturbation Theory in finite-temperature gauge-field theory.
It is instructive to compare this point of view to the ones suggested in the
-algebraic analysis of [
11]. In [
11], in effect, the so-called
-construction is used to build the thermal representation space of states
, and no state vector can be found in it out of which to devise a
-like Perturbation Theory. A perturbativization of the situation is therefore proposed, based on eigen-states of
which may be thought of in analogy to
unstable particle states. To some extent, these states may be viewed as corresponding to the quasi-particle excitations of the plasma, as long identified within the
-resummation program at scale
, and later on within the ultra-soft effective theory at scale
. In itself, using unstable/quasi-particle states instead of stable asymptotic ones doesn’t preclude a perturbative approach, in agreement with the perturbativization procedure evoked above, or the example explicitly worked out in [
30].
Now, it must also be realized that the point of view advocated in the current paper goes one step farther. Though possibly defined in a gauge invariant way, at least when explicitly constructed, there is no systematic control of the quantum fluctuations which may contribute to a given thermal observable at a given order of approximation; and thus there is no possibility to claim at a complete answer at any order of approximation. In the context of non-zero temperature quantum field theories, perturbative attempts could very well be doomed to failure because of this difficulty in which case the whole idea of diagrammatic perturbation theory at finite temperature may have to be abandoned, as quoted in [
11]. This is not too surprising in the end if one considers that Feynman diagrams are essentially tied to particle representations, while thermal representations of interacting relativistic field theories have long been known to possess no particle structure [
11,
31].
There is more on this: for the sake of illustration, the example of QED was used in the current paper and the case of QCD claimed to be the same up to constant redefinitions, like Equation (
14). Now, the point is that the
expansion parameter of QED has no equivalent in QCD. That is, if the QED coupling may be considered small enough to devise clear-cut separated momentum scales, such as
etc.., it is not so in QCD where a clear-cut separation of momentum scales
T,
,
,
, etc. is deprived of any physical realisation even at very high
T, as checked up to
in the pure Yang–Mills case [
32]. A consequence is that, as noticed by J.P. Blaizot in 1999, a condition necessary for a consistent implementation of the renormalisation group
à la Wilson fails to be met. The assumption on which one can rely in the case of QED doesn’t extend to the case of QCD where, at best, it defines an academic game.
At the very least, it appears that high-temperature quantum field theories call for non-perturbative methods [
33]. A new
formalism was born very recently [
34], which, by construction, implementing temperature by compactifying a spatial rather than the time direction, could perhaps disentangle the infrared sector from the high temperature limit of covariant (not invariant!) quantum field theories and avoid, in this way, the infrared enhancement mechanism, (i.e., the enhancement of infrared contributions remaining after infrared singularities cancellations have been taken into account: In our examples, the integration measures
appearing in Equations (
89) and (
A3) are responsible for this enhancement. They are induced by the high-temperature expansion of Equation (
7)). To our knowledge, however, the new formalism has not yet been explored in this respect.
Here may be the place where to recall another approach, functional and non-perturbative, concerning a similar calculation [
35,
36], that of the full fermionic 2- point function
under conditions identical to those considered in this paper. The initial energy of the fermion is
E, and
measures the travelling duration of the energetic fermion of mass
m entering at time
the QED plasma thermalized at temperature
T. A Bloch–Nordsieck estimate of
reads as
where
With
, and
the fermionic Compton wavelength, the momentum damping rate is
It is worth observing the large variety of different regimes and behaviours comprised within this non-perturbative Bloch–Nordsieck estimate of , according to the relative values of the involved parameters, and as time goes by. A richness that, at least, the perturbative attempts will have lead one to suspect.
It is worth mentioning also that another alternative to the perturbative loop expansion can be found within the approach developed in [
37], starting out with an explicit construction of the thermal ground state in terms of certain (anti) self-dual Yang–Mills configurations known as Harrington–Shepard (anti) calorons.