format: text

text:
# Nonlinear Landau damping of a plasmino in the quark-gluon plasma

## I Introduction

This work deals with the study of the interactions between soft-fermion and gluon excitations in an off-equilibrium quark-gluon plasma (QGP) within the hard thermal loop effective theory. The specific problem of the damping of soft, collective, fermion excitations is ultimately considered. This problem has already been addressed in the literature (see the review below) in the formalism of thermal perturbation theory, but nevertheless it would be interesting to see how it emerges in a real-time formalism based on kinetic equations for soft excitations. In addition, the calculation of the damping rate for moving soft excitations is still incomplete.

The theoretical investigation of properties of plasma excitations bearing fermion quantum numbers in a quark-gluon plasma originated in the pioneering works by Klimov and Weldon klim . The gauge-independent dispersion law of fermion excitations was first derived in a high-temperature region within the framework of the imaginary time formalism (Matsubara technique) by Klimov and within the framework of the real-time formalism (Keldysh-Swinger technique) by Weldon at one-loop order. The quark propagator in the case of nonzero temperature and massless fermions was shown to have two poles corresponding to two different dispersion relations, both with a positive energy. The first branch describes normal-particle excitations with the relation between chirality and helicity at zero temperature. The second branch is collective excitations, where the usual relation between chirality and helicity is flipped pis1 . It has been called a plasmino bra  to emphasize that, like the plasmon mode of the gluons, it is a purely collective branch of plasma excitations that have no analog at zero temperature.

The damping rate of plasma modes is one of the most important characteristics of high-temperature plasma. The calculation of the damping rate of fermion excitations in the rigid one-loop approximation was made by Kapusta and Toimela kap . They showed that the standard one-loop calculations were incomplete and resulted in an explicit dependence of the damping rate on the gauge-fixing condition and were even negative in some gauges. Overcoming this problem both in the case of fermion and boson modes is known to lead to the development of an effective perturbation theory put forth by Braaten and Pisarski bra1  and by Frenkel and Taylor fre  that enables one to resum in a systematic way high-loop diagrams contributing to damping rates at the leading order in the coupling constant $`g`$ and gives formal proof that resummation produces gauge-invariant results for the damping rates of both quarks and gluons. The calculation of the damping rate of a heavy fermion pis2  by Pisarski is one of the first examples of the application of effective perturbation theory.

The damping rates depend in a nontrivial way on how fast the quasiparticle is moving through the thermal medium. Within the framework of the imaginary-time formalism using a dispersion relation method, Braaten and Pisarski bra2  and independently Kobes et al. kob  present by somewhat distinct approaches a complete calculation of the quark damping rate at zero momentum at leading order in $`g`$ by evaluating one-loop diagrams constructed out of the effective propagators and vertices. It is shown that the damping rate of quark modes at rest is some pure number times $`g^2T`$, where $`T`$ is a temperature.

Much interest has been taken in the so-called damping rate of moving (hard or soft) quarks in hot QCD plasma111Here, we have not discussed the problems associated with calculation of the damping rate of moving fermions in a hot QED plasma, which possesses proper specific features (see, e.g. Ref. leb ).. The hard thermal loop (HTL) resummation renders the infrared finite damping rate of plasma excitations at rest. However, the damping rate of fast moving collective excitations is infrared divergent, arising from the exchange of quasistatic, magnetic gluons, reflecting the breakdown of the Braaten-Pisarski resummation scheme in deriving the physical quantities that are sensitive to $`O(g^2T)`$ corrections. In hot QCD plasma this problem is avoided by the introduction of the magnetic screening mass of order $`g^2T`$, which is expected to come from soft-gluon mutual interactions lin . An accurate calculation of the damping rate for soft moving quarks and gluons has been done by Pisarski in Ref. pis3 . He found that the infrared-singular contributions were proportional to the group velocity of the respective quark and gluon modes times $`g^2T\mathrm{ln}(1/g)`$. In Ref. fle1  Flechsig et al. presented a unified treatment of all next-to-leading-order corrections to fermion and gluon dispersion laws in QGP’s, which are infrared singular due to mass-shell singularities. They reproduced the result of Pisarski pis3  for moving quasiparticle excitations.

The need for further resummation beyond the scheme of Braaten and Pisarski also arises in an attempt to calculate the production rate of soft real photons from equilibrium QGP’s bai . This rate was shown to be logarithmically divergent owing to the so-called collinear (mass) singularities of a HTL diagram. The reason for the collinear singularities is known to be the use of bare (massless) hard-quark propagators that are on the mass shell. In a somewhat different context a similar difficulty was encountered in the research of the behavior of the next-to-leading order longitudinal component of the polarization function of the hot gluon system near the light cone fle2 . Here, the singularity is even more strong and generated by a premature restriction to soft loop momenta (the detailed study of the light cone singularity in the simpler gauge theory of scalar electrodynamics was made by Kraemmer et al. kra ).

The way to solve the problems connected with collinear singularities in hot QCD plasma was proposed by Flechsig and Rebhan fle3 . They showed that the inclusion of asymptotic thermal masses for hard transverse gluons222Notice that the fact of acquiring the transverse gluons in hot QCD plasmas of the thermal mass was first mentioned by Shuryak shu . and hard quarks removed the collinear singularities of the HTL amplitudes without spoiling gauge invariance. The examples of employment of the effective hard-quark propagator and effective photon–(hard)-quark vertex for deriving the finite soft-photon production rate from the equilibrium QGP can be found in Refs. fle3 ; aur .

Here we would like to represent a somewhat different view about the study of the above-mentioned problems, based solely on the kinetic theory. More comprehensively, our purpose is to construct kinetic equations that describe the evolution of off-equilibrium deviations in the number densities of the soft excitations (fermions and gluons), and to obtain damping rates from the collision terms. As we have shown in the case of soft-boson modes mar1 ; mar2 , there is no direct connection between the results obtained within the framework of kinetic theory of soft excitations and the results derived within the framework of effective perturbation theory. Here, we would like to extend the analysis, carried out in Refs. mar1 ; mar2  on the fermion degree of freedom of plasma excitations in hot QCD plasma. Our approach is based on the fundamental system of equations derived by Blaizot and Iancu bla2 . These equations are obtained on the basis of a truncation of the Schwinger-Dyson hierarchy. They isolate consistently the dominant terms in the coupling constant in a set of equations that describe the response of the plasma to weak and slowly varying disturbances, and encompass all HTL’s. This system of equations is determined on the space-time scale $`(gT)^1`$. We show how to obtain the information connected with the damping of soft excitation QGP’s by using this system, i.e., with the processes determined on the other space-time scale.

In order to allow for the collective excitations with arbitrary quantum numbers, Blaizot and Iancu considered both fermion and boson soft fields. The pure gauge sector (when the influence of the soft-fermion field is neglected) of this system of equations has been studied in detail. The sector connected with a soft-fermion degree of freedom of plasma excitations is not actually studied. In the present paper we have shown that this sector of plasma excitations contains some extended information on the dynamics not only of soft-fermion modes, which are especially interesting for us, but also soft-boson modes interacting with them. Deriving a self-consistent system of the kinetic equations for plasminos and plasmons enables us to have a new look at the known results obtained from the diagrammatic technique. It also allows us to represent them in a more visual form and to obtain some new results, in particular, connected with the possibility to observe in more detail the pumping dynamics of plasma excitations energy over the oscillation spectrum and energy pumping from the fermion branch of excitation into the boson branch and vice versa. Our strategy for deriving the relevant kinetic equations is similar to that already used in connection with the purely (soft) gauge sector mar2 .

From most papers close to the subject of our research the recent work by Niégawa nieg  is to be mentioned. Within the framework of the closed-time-path formalism the generalized Boltzmann equation that describes the evolution of the number-density functions of fermion quasiparticles was derived by him. The transport equations for “normal and abnormal modes” emerge here from the requirement of the absence of large contributions (which is the result of pinch singularity) of perturbative scheme proposed in Ref. nieg . In spite of the generality of a suggested approach, particular expressions for collision terms unfortunately were not given. It makes it difficult to establish the connection with other approaches, also dealing with the problem of relaxation of plasma excitations, in particular, with hard thermal loop effective theory.

The paper is organized as follows. In Sec. II, after summarizing the conventions and the notations used in this paper, we write a starting coupling set of the soft-field equations and dynamical equations describing the motion of hard particles in the presence of the soft quark and gluon fields. The approach to the solution of a given nonlinear system based on the approximation scheme method–the weak-field expansion–is proposed. In Sec. III, the first-order approximation of the induced source is considered, and the correlation functions of fermion and boson excitations are introduced. In Sec. IV, the second- and third-order approximations of the induced source and the color current induced by the soft-fermion field are studied. In Sec. V, we discuss the consistency with gauge symmetry of the approximation scheme used. In Sec. VI we derive the generalized kinetic equation for soft-fermion excitations and supplement it with a generalized kinetic equation for soft-boson excitations. The right-hand side of the last equation contains the terms responsible for the nonlinear interaction between soft bosonic and fermionic modes. Sec. VII presents a detailed consideration of the terms on the right-hand side of the generalized kinetic equation obtained in the previous section, their identification with specific physical processes. Then, kinetic equations describing the nonlinear interactions between plasminos and plasmons are extracted. In Sec. VIII, by means of the effective Ward identities, we show that the function defining the plasmino–hard-particle scattering matrix element squared is the same in the covariant gauge as in the temporal one. In Sec. IX, the transformation of the nonlinear Landau damping rate for a plasmino similar to the transformation in a pure gauge sector of soft QGP excitations is performed. In Sec. X, on the basis of decomposition of the nonlinear Landau damping rate into positive and negative parts, the model problem of two interacting infinitely narrow packets with fermion and boson quantum numbers is considered, and the kinematic relations between wave vectors of excitations are defined, so that one or another process of pumping over of energy occurs. In Sec. XI, an explicit expression for the nonlinear Landau damping rate of the plasmino at rest is derived, and a comparison with a similar expression for the damping rate obtained within the framework of the HTL approach is carried out. In Sec. XII, the behavior of the nonlinear Landau damping rate for normal-particle excitations near the light cone is considered, and this damping rate is shown to diverge owing to the mass singularities of quarks. In the Conclusion we briefly discuss a more general statement of the research problem of soft excitation dynamics in hot QCD plasma within the framework of kinetic approach outlined in this paper.

## II Blaizot-Iancu equations

We adopt conventions of Blaizot and Iancu bla2 . We use the metric $`g^{\mu \nu }=\mathrm{diag}(1,1,1,1)`$, choose units such that $`c=k_B=1`$ and note $`X=(X_0,𝐗),p=(p_0,𝐩),\mathrm{}.`$ We consider an $`\mathrm{SU}(N_c)`$ gauge theory with $`N_f`$ flavors of massless quarks. The color indices for the adjoint representation, $`a,b,\mathrm{}`$ run from 1 to $`N_c^21`$, while those for the fundamental representation, $`i,j,\mathrm{},`$ run from 1 to $`N_c`$. The Greek indices $`\alpha ,\beta ,\mathrm{}`$ for the spinor representation run from 1 to 4.

The quark field $`\psi `$ and the gauge field $`A_\mu =A_\mu ^at^a`$ with $`N_c^21`$ Hermitian generators in the fundamental representation obey the field equations

$$i\overline{)}D\psi (X)=\eta (X),$$
(2.1)
$$[D^\nu ,F_{\mu \nu }(X)]\xi ^1_\mu ^\nu A_\nu (X)=g\overline{\psi }(X)\gamma _\mu t^a\psi (X)t^a+j_\mu (X),$$
(2.2)
where $`\overline{)}D=\gamma ^\mu D_\mu `$ is the covariant derivative in the fundamental representation:

$$D_\mu =_\mu +igA_\mu (X),$$

$`F_{\mu \nu }=F_{\mu \nu }^at^a`$ is field strength tensor with $`F_{\mu \nu }^a=_\mu A_\nu ^a_\nu A_\mu ^agf^{abc}A_\mu ^bA_\nu ^c`$; $`[,]`$ denotes the commutator and $`\xi `$ is a gauge parameter fixing a covariant gauge.

The induced source $`\eta (X)`$ on the right-hand side of Eq. (2.1) can be written as an integral over hard momentum $`𝐤`$ of a density $`\overline{)}\mathrm{\Lambda }(𝐤,X)t^a\gamma ^\mu \mathrm{\Lambda }_\mu ^a(𝐤,X)`$ bla2

$$\eta (X)=g\frac{\mathrm{d}𝐤}{(2\pi )^3}\frac{1}{ϵ_k}\overline{)}\mathrm{\Lambda }(𝐤,X),ϵ_k|𝐤|.$$
(2.3)
A density $`\mathrm{\Lambda }_\mu ^a(𝐤,X)`$ is a generalized one-body density matrix mixing fermion and boson degrees of freedom.

The total induced color current $`j_\mu `$ on the right-hand side of field equation (2.2) is a sum of two parts: the current $`j_\mu ^A(X)`$ induced by the soft-gauge field

$$j_\mu ^A(X)=g\frac{\mathrm{d}𝐤}{(2\pi )^3}v_\mu [N_f(\delta n_+^A(𝐤,X)\delta n_{}^A(𝐤,X))+2N_c\delta N^A(𝐤,X)],$$
(2.4)
and the current $`j_\mu ^\mathrm{\Psi }(X)`$ induced by the soft-quark field

$$j_\mu ^\mathrm{\Psi }(X)=g\frac{\mathrm{d}𝐤}{(2\pi )^3}v_\mu [\delta n_+^\mathrm{\Psi }(𝐤,X)\delta n_{}^\mathrm{\Psi }(𝐤,X)],v_\mu =(1,𝐯),𝐯=𝐤/|𝐤|,$$
(2.5)
where $`\delta n_\pm ^{A,\mathrm{\Psi }}(𝐤,X)=\delta n_\pm ^{A,\mathrm{\Psi }a}t^a`$ are soft fluctuations in the quark and antiquark color densities and $`\delta N^A(𝐤,X)=\delta N^{Aa}t^a`$ is a soft fluctuation in the gluon color density. The superscripts indicate the nature of the background field that induces fluctuations. On the space-time scale $`(gT)^1`$ these functions and density $`\overline{)}\mathrm{\Lambda }(𝐤,X)`$ satisfy the following system of equations:

$$[vD_X,\delta n_\pm ^A(𝐤,X)]=g𝐯𝐄(X)\frac{dn(ϵ_k)}{dϵ_k},$$
(2.6)
$$[vD_X,\delta N^A(𝐤,X)]=g𝐯𝐄(X)\frac{dN(ϵ_k)}{dϵ_k},$$
(2.7)
$$[vD_X,\delta n_\pm ^\mathrm{\Psi }(𝐤,X)]=\pm \frac{ig}{2ϵ_k}t^a(\overline{\psi }(X)t^a\overline{)}\mathrm{\Lambda }(𝐤,X)\overline{\overline{)}\mathrm{\Lambda }}(𝐤,X)t^a\psi (X)),$$
(2.8)
$$(vD_X)\overline{)}\mathrm{\Lambda }(𝐤,X)=igC_F[N(ϵ_k)+n(ϵ_k)]\overline{)}v\psi (X).$$
(2.9)
Here, $`𝐄(X)=𝐄^a(X)t^a`$ is the chromoelectric field, $`E^i=F^{i0};n(ϵ_k)=1/(\mathrm{exp}(ϵ_k/T)+1)`$ and $`N(ϵ_k)=1/(\mathrm{exp}(ϵ_k/T)1)`$ are the fermion and boson occupation factors, $`T`$ is the temperature of the plasma, and $`C_F`$ is the quadratic Casimir invariant of the fundamental representation. The function $`\overline{\overline{)}\mathrm{\Lambda }}(𝐤,X)`$ is related to $`\overline{)}\mathrm{\Lambda }(𝐤,X)`$ by $`\overline{\overline{)}\mathrm{\Lambda }}(𝐤,X)=\overline{)}\mathrm{\Lambda }^{}(𝐤,X)\gamma ^0`$, where the Hermitian conjugation refers to color and spinor indices.

The self-consistent system of Eqs. (2.1) – (2.9) for the soft fluctuating fields $`\psi `$ and $`A_\mu `$ and their induced sources was first derived in Ref. bla2 . In the subsequent discussion it will be called the Blaizot-Iancu equations. However, as distinct from the original paper bla2 , we take Eqs. (2.6) – (2.9) not as kinetic equations, i.e., time-irreversible ones, but as exact ”microscopic” dynamical equations coming from the hard thermal loop effective action and describing the evolution of the collisionless plasma with zero expectation values of the soft fields (or the associated induced sources).

The Blaizot-Iancu equations are solved by the approximation scheme method – the weak-field expansion. For this purpose first we expand soft fluctuations of the (anti)quark and gluon color densities as power series in the oscillation amplitudes of the functions $`A_\mu `$ and $`\psi `$,

$$\delta n_\pm ^A=\underset{n=1}{\overset{\mathrm{}}{}}\delta n_\pm ^{A(n)},\delta N^A=\underset{n=1}{\overset{\mathrm{}}{}}\delta N^{A(n)},$$
(2.10)
$$\delta n_\pm ^\mathrm{\Psi }=\underset{n=0}{\overset{\mathrm{}}{}}\delta n_\pm ^{\mathrm{\Psi }(n,2)},$$
(2.11)
$$\overline{)}\mathrm{\Lambda }=\underset{n=0}{\overset{\mathrm{}}{}}\overline{)}\mathrm{\Lambda }^{(n,1)},$$
(2.12)
where the index $`n`$ shows that $`\delta n_\pm ^{A(n)},\delta N^{A(n)},\delta n_\pm ^{\mathrm{\Psi }(n,2)}`$ and $`\overline{)}\mathrm{\Lambda }^{(n,1)}`$ are proportional to the $`n`$th power of $`A_\mu `$. Since the fermion fields appear in explicit form on the right-hand side of Eqs. (2.8) and (2.9), the functions $`\delta n_\pm ^{\mathrm{\Psi }(n,2)}`$ and $`\overline{)}\mathrm{\Lambda }^{(n,1)}`$ also depend on amplitudes $`\overline{\psi }`$ and $`\psi `$. By virtue of the structure of the right-hand side of Eqs. (2.8) and (2.9), we have

$$\delta n_\pm ^{\mathrm{\Psi }(n,2)}\overline{\psi }\psi ,\overline{)}\mathrm{\Lambda }^{(n,1)}\psi $$

for arbitrary values of $`n`$.

The induced color currents (2.4) and (2.5) and the induced source (2.3) are expressed as

$$j_\mu ^A=\underset{n=1}{\overset{\mathrm{}}{}}j_\mu ^{A(n)},$$
(2.13)
$$j_\mu ^\mathrm{\Psi }=\underset{n=0}{\overset{\mathrm{}}{}}j_\mu ^{\mathrm{\Psi }(n,2)},$$
(2.14)
where

$$j_\mu ^{A(n)}=g\frac{\mathrm{d}𝐤}{(2\pi )^3}v_\mu [N_f(\delta n_+^{A(n)}\delta n_{}^{A(n)})+2N_c\delta N^{A(n)}],$$
(2.15)
$$j_\mu ^{\mathrm{\Psi }(n,2)}=g\frac{\mathrm{d}𝐤}{(2\pi )^3}v_\mu [\delta n_+^{\mathrm{\Psi }(n,2)}\delta n_{}^{\mathrm{\Psi }(n,2)}],$$
(2.16)
and

$$\eta =\underset{n=0}{\overset{\mathrm{}}{}}\eta ^{(n,1)},$$
(2.17)
where

$$\eta ^{(n,1)}=g\frac{\mathrm{d}𝐤}{(2\pi )^3}\frac{1}{ϵ_k}\overline{)}\mathrm{\Lambda }^{(n,1)}.$$
(2.18)
Now we turn to the field equations (2.1) and (2.2), connecting the quark and gluon soft fields with induced source $`\eta `$ and total color current $`j_\mu `$. Let us rewrite these equations, explicitly separating the free parts in Eqs. (2.1) and (2.2) from interaction terms. Taking into account the expansion of induced current (2.13), (2.14) and induced source (2.17), we have

$$i\overline{)}\psi \eta ^{(0,1)}=g\overline{)}A\psi +\eta _{NL},$$
(2.19)
$$_\mu (F^{\mu \nu })_L\xi ^1^\nu ^\mu A_\mu j^{A(1)\nu }=j_{NL}^\nu +g\overline{\psi }\gamma ^\nu t^a\psi t^aig_\mu [A^\mu ,A^\nu ]ig[A_\mu ,(F^{\mu \nu })_L]+g^2[A_\mu ,[A^\mu ,A^\nu ]].$$
(2.20)
Here, the indices $`L`$ and $`NL`$ denote the linear and nonlinear parts of strength tensor, the induced source and the color-induced current with respect to $`A_\mu `$ and $`\psi `$. The approximation that will be made is the truncation of the series expansion (2.13), (2.14) and (2.17). The accuracy of this approximation is controlled by the characteristic amplitudes of the soft fields, which will be discussed in Sec. V. As will be shown to account for the nonlinear interaction between waves and hard particles in QGP at leading order in $`g`$, it is sufficient to restrict the consideration to the third order in total powers of $`\psi `$ and $`A_\mu `$ in expansions (2.13), (2.14), and (2.17).

## III Linear approximation of the induced source $`\eta `$. Correlation functions of the fermion and boson excitations

We now come to the derivation of the kinetic equation for the soft-fermion modes. The starting point is Eq. (2.19). Its left-hand side contains the linear approximation of the induced source $`\eta `$, whose explicit form can be easily found from Eq. (2.9). We prefer to work in momentum space. The corresponding equations are obtained by using

$$\psi (X)=𝑑q\psi (q)\mathrm{e}^{iqX},A_\mu (X)=𝑑pA_\mu (p)\mathrm{e}^{ipX}$$

and similar transformations for $`\delta n_\pm ^A,\delta N^A`$, etc. Here and in what follows we denote the momenta of the soft-quark fields by $`q,q^{},q_1,\mathrm{}`$ and the momenta of the soft-gauge fields by $`p,p^{},p_1,\mathrm{}`$.

Let us linearize Eq. (2.9) with respect to $`\psi `$ and $`A_\mu `$. The result of the Fourier transformation of the linearized equation is

$$\overline{)}\mathrm{\Lambda }^{(0,1)}(𝐤,q)=gC_F\frac{\overline{)}v\psi (q)}{vq+iϵ}[N(ϵ_k)+n(ϵ_k)],ϵ+0.$$
(3.1)
By substituting Eq. (3.1) into the Fourier transform of relation (2.18) (for $`n=0`$) and performing the radial integration over $`dϵ_k`$, we obtain the linear approximation with respect to the fermion field of the induced source,

$$\eta ^{(0,1)}(q)=\delta \mathrm{\Sigma }(q)\psi (q),$$
(3.2)
where

$$\delta \mathrm{\Sigma }(q)=\omega _0^2\frac{d\mathrm{\Omega }}{4\pi }\frac{\overline{)}v}{vq+iϵ}$$

is a well-known HTL expression for the soft-quark (retarded) self-energy, $`\omega _0^2=g^2C_FT^2/8`$ is the plasma frequency of the quark sector of plasma excitations, and $`d\mathrm{\Omega }`$ is an angular measure.

Furthermore, we rewrite Eq. (2.19) in momentum space. Taking into account Eq. (3.2) we obtain

$${}_{}{}^{}S_{}^{1}(q^{})\psi ^j(q^{})=g\overline{)}A(p_1)\psi ^j(q_1)\delta (q^{}q_1p_1)𝑑q_1𝑑p_1\eta ^{(1,1)j}(q^{})\eta ^{(2,1)j}(q^{}).$$
(3.3)
Here $`{}_{}{}^{}S(q)=[\overline{)}q+\delta \mathrm{\Sigma }(q)]^1`$ is an equilibrium propagator for a soft quark, corrected to leading order in $`g`$.

Let us multiply Eq. (3.3) by the Dirac conjugate amplitude $`\overline{\psi }^i(q)=(\psi ^i(q))^{}\gamma ^0`$ and take an expectation value over the off-equilibrium ensemble:

$$\overline{\psi }^i(q)^{}S^1(q^{})\psi ^j(q^{})=g\overline{\psi }^i(q)(\overline{)}A(p_1)\psi (q_1))^j\delta (q^{}q_1p_1)𝑑q_1𝑑p_1$$
(3.4)
$$\overline{\psi }^i(q)\eta ^{(1,1)j}(q^{})\overline{\psi }^i(q)\eta ^{(2,1)j}(q^{}).$$

We now introduce the correlation functions of soft-fermion and -boson excitations

$$\mathrm{{\rm Y}}_{\beta \alpha }^{ji}(q^{},q)=\overline{\psi }_\alpha ^i(q)\psi _\beta ^j(q^{}),I_{\mu \nu }^{ab}(p^{},p)=A_\mu ^a(p^{})A_\nu ^b(p),$$
(3.5)
respectively. The asterisk denotes a complex conjugate. The considered soft excitations are necessarily colorless and have zero fermion number by virtue of the fact that mean fields $`A_\mu ^a`$ and $`\psi _\alpha ^i`$ or the associated mean induced color current and source are assumed to be vanishing. Therefore, for the physical situation of interest, the off-equilibrium two-point functions (3.5) are diagonal in color space, which will be implied in what follows.

For the conditions of stationary and homogeneous QGP’s (i.e., when correlation functions (3.5) in the coordinate representation depend only on the difference of coordinates and time $`\mathrm{}X=X^{}X`$), we have

$$\mathrm{{\rm Y}}_{\beta \alpha }^{ji}(q^{},q)=\mathrm{{\rm Y}}_{\beta \alpha }^{ji}(q^{})\delta (q^{}q),I_{\mu \nu }^{ab}(p^{},p)=I_{\mu \nu }^{ab}(p^{})\delta (p^{}p).$$

The QGP state becomes slightly heterogeneous and nonstationary because of the effects of the nonlinear interaction between waves and particles. This leads to a $`\delta `$-function broadening, and $`\mathrm{{\rm Y}}_{\beta \alpha }^{ji}`$ and $`I_{\mu \nu }^{ab}`$ depend on both arguments.

Let us introduce $`\mathrm{{\rm Y}}_{\beta \alpha }^{ji}(q^{},q)=\mathrm{{\rm Y}}_{\beta \alpha }^{ji}(q,\mathrm{}q)`$, $`\mathrm{}q=q^{}q`$ with $`|\mathrm{}q/q|1`$ and $`I_{\mu \nu }^{ab}(p^{},p)=I_{\mu \nu }^{ab}(p,\mathrm{}p)`$, $`\mathrm{}p=p^{}p`$ with $`|\mathrm{}p/p|1`$, and insert the correlation functions in the Wigner form

$$\mathrm{{\rm Y}}_{\beta \alpha }^{ji}(q,x)=\mathrm{{\rm Y}}_{\beta \alpha }^{ji}(q,\mathrm{}q)\mathrm{e}^{i\mathrm{}qx}𝑑\mathrm{}q,I_{\mu \nu }^{ab}(p,x)=I_{\mu \nu }^{ab}(p,\mathrm{}p)\mathrm{e}^{i\mathrm{}px}𝑑\mathrm{}p,$$
(3.6)
slowly depending on $`x`$. In Eq. (3.4), we replace $`qq^{},ij`$, take a complex conjugation, and then subtract the resulting equation from Eq. (3.4), expanding beforehand the quark self-energy into “Hermitian” and “anti-Hermitian” parts:

$$\delta \mathrm{\Sigma }(q)=\delta \mathrm{\Sigma }^\mathrm{H}(q)+\delta \mathrm{\Sigma }^\mathrm{A}(q),\gamma ^0(\delta \mathrm{\Sigma }^\mathrm{H}(q))^{}\gamma ^0=\delta \mathrm{\Sigma }^\mathrm{H}(q),\gamma ^0(\delta \mathrm{\Sigma }^\mathrm{A}(q))^{}\gamma ^0=\delta \mathrm{\Sigma }^\mathrm{A}(q).$$

We assume that the anti-Hermitian part $`\delta \mathrm{\Sigma }^\mathrm{A}`$ is small relative to the Hermitian part $`\delta \mathrm{\Sigma }^\mathrm{H}`$ and is of the same order as the nonlinear terms on the right-hand side. We can therefore set $`\delta \mathrm{\Sigma }^\mathrm{A}(q)\delta \mathrm{\Sigma }^\mathrm{A}(q^{})`$ and move the term with $`\delta \mathrm{\Sigma }^\mathrm{A}`$ into the right-hand side of Eq. (3.4). We expand the remaining terms in the left-hand side with respect to $`\mathrm{}q`$ up to the first order. This corresponds to a gradient expansion procedure usually used in the derivation of kinetic equations. Multiplying the resulting equation by $`\mathrm{e}^{i\mathrm{}qx}`$ and integrating over $`\mathrm{}q`$, we obtain

$$\mathrm{tr}\left(\frac{}{q^\mu }[\overline{)}q+\delta \mathrm{\Sigma }^\mathrm{H}(q)]\frac{\mathrm{{\rm Y}}^{ji}(q,x)}{x_\mu }\right)=2i\mathrm{tr}(\delta \mathrm{\Sigma }^\mathrm{A}(q)\mathrm{{\rm Y}}^{ji}(q,x))$$
(3.7)
$$ig𝑑q^{}𝑑q_1𝑑p_1\{(\overline{\psi }(q_1)\overline{)}A(p_1))^i\psi ^j(q^{})\delta (qq_1p_1)\overline{\psi }^i(q)(\overline{)}A(p_1)\psi (q_1))^j\delta (q^{}q_1p_1)\}$$

$$i𝑑q^{}\{(\overline{\eta }^{(1,1)i}(q)\psi ^j(q^{})\overline{\psi }^i(q)\eta ^{(1,1)j}(q^{}))+(\overline{\eta }^{(2,1)i}(q)\psi ^j(q^{})\overline{\psi }^i(q)\eta ^{(2,1)j}(q^{}))\}.$$

Here, the Dirac trace is represented by $`\mathrm{tr}`$, $`\overline{\eta }(q)=\eta ^{}(q)\gamma ^0`$, and we take into account the reality of the gauge field: $`A_\mu ^a(p)=A_\mu ^a(p)`$. The linear term with $`\delta \mathrm{\Sigma }^\mathrm{A}`$ on the right-hand side corresponds to the linear Landau damping of soft-fermion excitations.

## IV Second and third approximations of the induced source $`\eta `$

Now we are concerned with computation of the nonlinear corrections to the induced source on the right-hand side of Eq. (3.7). At first, we define $`\overline{)}\mathrm{\Lambda }^{(1,1)}`$. Substituting the series expansion (2.12) into (2.9) and keeping only the terms of the second order in fields, we find

$$(v_X)\overline{)}\mathrm{\Lambda }^{(1,1)}(𝐤,X)=ig(vA(X))\overline{)}\mathrm{\Lambda }^{(0,1)}(𝐤,X).$$

Performing the Fourier transformation of the last equation and taking into account an explicit form of $`\overline{)}\mathrm{\Lambda }^{(0,1)}`$ (3.1), we derive

$$\overline{)}\mathrm{\Lambda }^{(1,1)}(𝐤,q)=g^2C_F[N(ϵ_k)+n(ϵ_k)]\frac{v^\mu \overline{)}v}{vq+iϵ}t^a\frac{1}{vq_1+iϵ}A_\mu ^a(p_1)\psi (q_1)\delta (qq_1p_1)𝑑q_1𝑑p_1.$$
(4.1)
Substituting the obtained expression into Eq. (2.18) (for $`n=1`$) and performing the radial integration over $`dϵ_k`$, we find the required induced source correction

$$\eta ^{(1,1)}(q)=g\omega _0^2t^a\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{v^\mu \overline{)}v}{(vq+iϵ)(vq_1+iϵ)}A_\mu ^a(p_1)\psi (q_1)\delta (qq_1p_1)𝑑q_1𝑑p_1.$$
(4.2)
The expression for the induced source of the third order in the fields is defined by means of reasoning similar to the previous one. Here, we have

$$\eta ^{(2,1)}(q)=g^2\omega _0^2t^at^b\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{v^\mu v^\nu \overline{)}v}{(vq+iϵ)(v(qp_1)+iϵ)(vq_1+iϵ)}A_\mu ^a(p_1)A_\nu ^b(p_2)\psi (q_1)\delta (qq_1p_1p_2)𝑑q_1𝑑p_1𝑑p_2.$$
(4.3)
Now we return to the initial equation for the soft-fermion field (3.7). We substitute nonlinear corrections to the induced source $`\eta `$ (Eqs. (4.2) and (4.3)) into the right-hand side of this equation. After simple algebraic transformations, instead of Eq. (3.7), we find

$$\mathrm{tr}\left(\frac{}{q^\mu }[\overline{)}q+\delta \mathrm{\Sigma }^\mathrm{H}(q)]\frac{\mathrm{{\rm Y}}^{ji}(q,x)}{x_\mu }\right)=2i\mathrm{tr}(\delta \mathrm{\Sigma }^\mathrm{A}(q)\mathrm{{\rm Y}}^{ji}(q,x))$$

$$igdq^{}dq_1dp_1\{\overline{\psi }^i(q)A^{a\mu }(p_1)(^{}\mathrm{\Gamma }_\mu ^{(Q)a}(p_1;q_1,q^{})\psi (q_1))^j\delta (q^{}q_1p_1)$$

$$(\overline{\psi }(q_1)^{}\mathrm{\Gamma }_\mu ^{(Q)a}(p_1;q_1,q))^iA^{a\mu }(p_1)\psi ^j(q^{})\delta (qq_1p_1)$$
(4.4)
$$+ig^2dq^{}dq_1dp_1dp_2\{A^{a\mu }(p_1)A^{b\nu }(p_2)\overline{\psi }^i(q)(\delta \mathrm{\Gamma }_{\mu \nu }^{(Q)ab}(p_1,p_2;q^{},q_1)\psi (q_1))^j\delta (q^{}q_1p_1p_2)$$

$$A^{a\mu }(p_1)A^{b\nu }(p_2)(\overline{\psi }(q_1)\delta \mathrm{\Gamma }_{\mu \nu }^{(Q)ba}(p_2,p_1;q,q_1))^i\psi ^j(q^{})\delta (qq_1p_1p_2)\}.$$

Here,

$${}_{}{}^{}\mathrm{\Gamma }_{\mu }^{(Q)a}(p;q_1,q_2)=t^a{}_{}{}^{}\mathrm{\Gamma }_{\mu }^{(Q)}(p;q_1,q_2)t^a(\gamma _\mu +\delta \mathrm{\Gamma }_\mu ^{(Q)}(p;q_1,q_2))$$
(4.5)
is an effective (i.e., HTL-resummed) vertex between a quark pair and a gluon, that is, a sum of the bare vertex $`\gamma _\mu `$ and a corresponding HTL correction bra1 ; tay ; bla2

$$\delta \mathrm{\Gamma }_\mu ^{(Q)}(p;q_1,q_2)=\omega _0^2\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{v_\mu \overline{)}v}{(vq_1+iϵ)(vq_2iϵ)},$$
(4.6)
and

$$\delta \mathrm{\Gamma }_{\mu \nu }^{(Q)ab}(p_1,p_2;q_1,q_2)=\omega _0^2\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{v_\mu v_\nu \overline{)}v}{(vq_1+iϵ)(vq_2iϵ)}\left(\frac{t^at^b}{v(q_1+p_1)+iϵ}+\frac{t^bt^a}{v(q_1+p_2)+iϵ}\right)$$
(4.7)
is an effective vertex between a quark pair and two gluons (this vertex does not exist at tree level, and in the leading order it arises entirely from the HTL bra1 ; tay ). The superscript $`(Q)`$ denotes that the vertex $`\delta \mathrm{\Gamma }_\mu (X;Y_1,Y_2)`$ corresponds to the function (4.6) in the coordinate representation, where the time arguments satisfy $`Y_1^0X^0Y_2^0`$ (boundary conditions), the vertex function (4.7) corresponds to $`\delta \mathrm{\Gamma }_\mu (X_1,X_2;Y_1,Y_2),`$ where the time arguments satisfy $`Y_1^0X_1^0X_2^0Y_2^0`$ for the first term in parentheses on the right-hand side of Eq. (4.7) and for the second term we have $`Y_1^0X_2^0X_1^0Y_2^0`$, i.e., the time argument of the external quark leg incoming in the vertex functions (4.6) and (4.7), is largest. In deriving Eq. (4.4) we have dropped the terms proportional to $`A_\mu \overline{\psi }A_\nu \psi `$ not contributing to the right-hand side of kinetic equation for the soft-fermion excitations at the leading order in $`g`$, in which we are interested.

Because of nonlinear wave interactions, phase correlation effects occur. By virtue of their smallness, the expectation value of a product of four quantities $`A_\mu A_\nu \overline{\psi }\psi `$ can be expressed approximately as a product of the expectation values of two fields $`A_\mu A_\nu `$ and $`\overline{\psi }\psi `$. For a product of three fields this approach yields zero; in this case, the weak correlation between the fields is to be taken into account. The third-order correlation functions on the right-hand side of Eq. (4.4) contain amplitudes of waves with different statistic, nonlinear corrections to that are defined by corresponding field equations. Let us consider first of all the nonlinear correction to a free quark field. For this purpose, we use field equation (3.3), keeping in the right-hand side only the second-order terms (with respect to $`A_\mu `$ plus $`\psi `$)

$${}_{}{}^{}S_{}^{1}(q)\psi (q)=g^{}\mathrm{\Gamma }^{(Q)\mu }(p_1;q_1,q)A_\mu (p_1)\psi (q_1)\delta (qq_1p_1)𝑑q_1𝑑p_1.$$

The approximate solution of this equation has the form

$$\psi (q)=\psi ^{(0)}(q)g^{}S(q)^{}\mathrm{\Gamma }^{(Q)\mu }(p_1;q_1,q)A_\mu ^{(0)}(p_1)\psi ^{(0)}(q_1)\delta (qq_1p_1)𝑑q_1𝑑p_1,$$
(4.8)
where $`A_\mu ^{(0)}`$ and $`\psi ^{(0)}`$ are solutions of the appropriate homogeneous field equations corresponding to free fields.

Now we consider the nonlinear correction to a free gauge field. For this purpose, we use field equation (2.20) is rewritten in momentum space, keeping on the right-hand side the terms of the second order in fields:

$${}_{}{}^{}𝒟_{}^{1\mu \nu }(p)A_\nu ^a(p)=j^{A(2)a\mu }(p)+j^{\mathrm{\Psi }(0,2)a\mu }(p)+g\overline{\psi }(q_1)\gamma ^\mu t^a\psi (q_2)\delta (p+q_1q_2)𝑑q_1𝑑q_2$$

$$\frac{i}{2}gf^{abc}\mathrm{\Gamma }^{\mu \nu \lambda }(p,p_1,p_2)A_\nu ^b(p_1)A_\lambda ^c(p_2)\delta (pp_1p_2)𝑑p_1𝑑p_2.$$
(4.9)
Here, $`{}_{}{}^{}𝒟_{}^{\mu \nu }(p)=[p^2g^{\mu \nu }(1+\xi ^1)p^\mu p^\nu +\delta \mathrm{\Pi }^{\mu \nu }(p)]^1`$ represents the medium-modified (retarded) gluon propagator with the soft-gluon polarization tensor,

$$\delta \mathrm{\Pi }^{\mu \nu }(p)=3\omega _{\mathrm{pl}}^2\left(g^{\mu 0}g^{\nu 0}p^0\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{v^\mu v^\nu }{vp+iϵ}\right),$$

corrected at the leading order in $`g`$, $`\omega _{\mathrm{pl}}^2=g^2T^2(N_f+2N_c)/18`$ is a plasma frequency of the gauge sector of plasma excitations, and $`\mathrm{\Gamma }^{\mu \nu \lambda }(p,p_1,p_2)`$ is a bare three-gluon vertex. The first and the last terms on the right-hand side of Eq. (4.9) give the contribution to a correlation function $`A\overline{\psi }\psi `$ proportional to $`AA\overline{\psi }\psi AA\overline{\psi }\psi `$. By replacing interacting fields by free ones and dividing the fourth-order correlation into a product of the second-order correlators, this contribution vanishes. Therefore these terms in this approximation can be dropped.

Now we consider the remaining terms on the right-hand side of Eq. (4.9). For determining the explicit form of $`j_\mu ^{\mathrm{\Psi }(0,2)a}`$ it is necessary to derive the expression $`\delta n_\pm ^{\mathrm{\Psi }(0,2)a}`$ according to Eq. (2.16). For this purpose, we substitute the expansion (2.11) into Eq. (2.8), keeping only the terms with quark fields $`\overline{\psi },\psi `$ on the right-hand side. Performing the Fourier transformation and taking into account Eq. (3.1), we find

$$\delta n_\pm ^{\mathrm{\Psi }(0,2)a}(𝐤,p)=\pm g^2C_F\frac{1}{2ϵ_k}[N(ϵ_k)+n(ϵ_k)]\frac{1}{(vq_1iϵ)(vq_2+iϵ)}\overline{\psi }(q_1)\overline{)}vt^a\psi (q_2)\delta (p+q_1q_2)𝑑q_1𝑑q_2.$$

Substituting the obtained expression into Eq. (2.16) and performing the radial integration over $`dϵ_k`$, we find the required current correction $`j_\mu ^{\mathrm{\Psi }(0,2)a}`$:

$$j_\mu ^{\mathrm{\Psi }(0,2)a}(p)=g\omega _0^2\frac{d\mathrm{\Omega }}{4\pi }\frac{v_\mu }{(vp+iϵ)}\left(\frac{1}{vq_2+iϵ}\frac{1}{vq_1iϵ}\right)\overline{\psi }(q_1)\overline{)}vt^a\psi (q_2)\delta (p+q_1q_2)𝑑q_1𝑑q_2.$$
(4.10)
Taking into account the discussion above, we obtain the following equation, instead of Eq. (4.9):

$${}_{}{}^{}𝒟_{}^{1\mu \nu }(p)A_\nu ^a(p)=g\overline{\psi }(q_1)^{}\mathrm{\Gamma }^{(G)\mu }(p;q_1,q_2)t^a\psi (q_2)\delta (p+q_1q_2)𝑑q_1𝑑q_2.$$
(4.11)
Here,

$${}_{}{}^{}\mathrm{\Gamma }_{\mu }^{(G)}(p;q_1,q_2)=\gamma _\mu +\delta \mathrm{\Gamma }_\mu ^{(G)}(p;q_1,q_2),$$
(4.12)
$$\delta \mathrm{\Gamma }_\mu ^{(G)}(p;q_1,q_2)=\omega _0^2\frac{d\mathrm{\Omega }}{4\pi }\frac{v_\mu \overline{)}v}{(vq_1iϵ)(vq_2iϵ)}$$

is the effective two-quarks–one-gluon vertex function, where now (in the coordinate representation) the time arguments satisfy $`X^0\mathrm{max}(Y_1^0,Y_2^0)`$ and the chronological order of $`Y_1`$ and $`Y_2`$ is arbitrary bla2 . The time argument of an external gluon leg incoming in the vertex is largest, as indicated by superscript $`(G)`$. The approximate solution of Eq. (4.11) (with accuracy required for our further calculations) is of the form

$$A_\mu ^a(p)=A_\mu ^{(0)a}(p)+g^{}𝒟_{\mu \nu }(p)\overline{\psi }^{(0)}(q_1)^{}\mathrm{\Gamma }^{(G)\nu }(p;q_1,q_2)t^a\psi ^{(0)}(q_2)\delta (p+q_1q_2)𝑑q_1𝑑q_2.$$
(4.13)
At the end of this section we present the expression for the next term in the expansion of current $`j_\mu ^\mathrm{\Psi }`$, which is needed for deriving general kinetic equation for soft-gluon modes in Sec. VI. Performing similar calculations, we obtain

$$j_\mu ^{\mathrm{\Psi }(1,2)a}(p)=g^2\omega _0^2\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{v^\mu v^\nu (\overline{)}v)_{\alpha \beta }}{vp+iϵ}[(T^a)^{bc}(t^c)^{ij}\frac{1}{(vq_1iϵ)(vq_2+iϵ)}$$

$$+(t^at^b)^{ij}\frac{1}{(v(q_2+p_1)+iϵ)(vq_2+iϵ)}(t^bt^a)^{ij}\frac{1}{(v(q_1p_1)iϵ)(vq_1iϵ)}]$$
(4.14)
$$\times A_\nu ^b(p_1)\overline{\psi }_\alpha ^i(q_1)\psi _\beta ^j(q_2)\delta (pp_1+q_1q_2)dp_1dq_1dq_2,$$

where $`(T^a)^{bc}=if^{abc}`$. The appearance of the term with the production of generators of the $`\mathrm{SU}(N_c)`$ group in other representations (the adjoint and the fundamental ones) is a special feature of the last expression. This term is not vanishing only in the processes of the higherorder in field powers, the processes of Boltzmann type, i.e., the scattering processes between the soft-fermion and -gluon excitations (see Conclusion).

## V Characteristic amplitudes of the soft fields

In this section we will estimate the typical amplitudes of the soft fields, both in coordinate and in momentum space, such that the truncation of the series expansions (2.13), (2.14), and (2.17) can be made and in the long run one can derive a closed system of gauge-invariant kinetic equations for soft plasma excitations.

For this purpose, let us discuss in more detail the approximation scheme used in this paper. In fact, the value $`\delta gRA_\mu `$ is a (dimensionless) parameter of the expansion in powers of nonlinearities in Eqs. (2.1) – (2.9), where $`R_X^1`$ is a typical length for off-equilibrium deviations. It follows, for example, from rewriting the equations of motion for the soft fields (2.1) and (2.2) in the form (2.19) and (2.20), where the gauge fields $`A_\mu (X)`$ are expanded from the covariant derivatives $`D_\mu =_\mu +igA_\mu (X)`$ and the field strength tensor $`F_{\mu \nu }(X)`$. The disadvantage of the proposed iterated scheme for deriving the terms in expansions (2.13), (2.14), and (2.17) is explicit breaking of the non-Abelian gauge symmetry of the dynamical equations for the soft fluctuations (2.6) – (2.9) at each step of the approximate calculations. Here, the answer to the nontrivial question on the recovery of a gauge symmetry depends directly on the order of the magnitudes of the soft fields, the dynamics of which we would like to describe. Below, we discuss two relevant cases.

Let us assume that the amplitudes of soft-boson and -fermion excitations are of orders $`g`$ coinciding with the orders of $`A_\mu `$ and $`\psi `$ presented in Ref. bla2 , i.e.,

$$|A_\mu (X)|T(|A_\mu (p)|1/g(gT)^3),$$

$$|\overline{\psi }(X)||\psi (X)|gT^3(|\overline{\psi }(q)||\psi (q)|1/g^2(gT)^5).$$

In this case, by using the obtained expressions for the terms in the expansion of the induced source $`\eta `$ (3.2), (4.2), and (4.3) and the color currents $`j_\mu ^A`$ mar2  and $`j_\mu ^\mathrm{\Psi }`$ (4.10) and (4.14), we have the following estimates:

$$\eta ^{(0,1)}(q)\eta ^{(1,1)}(q)\eta ^{(2,1)}(q)\mathrm{}gT|\psi (q)|,$$

$$j_\mu ^{A(1)}(p)j_\mu ^{A(2)}(p)j_\mu ^{A(3)}(p)\mathrm{}\frac{1}{g^2T},$$

$$j_\mu ^{\mathrm{\Psi }(0,2)}(p)j_\mu ^{\mathrm{\Psi }(1,2)}(p)\mathrm{}g(gT)^4|\overline{\psi }(q)||\psi (q)|\left(\frac{1}{g^2T}\right).$$

Thus, all terms in each expansion (2.13), (2.14), and (2.17) are of the same order of magnitude and the problem of resummation of all the relevant contributions arises. Thus, gauge symmetry is recovered. In a pure-gauge case such an approach is represented, for instance, in the work of Jackiw and Nair jac  dealing with the derivation of the non-Abelian version of the Kubo formula. Here, the explicit expression for the induced current $`j_\mu ^A(p)`$, including the contribution of all higher point functions, that actually coincide with the iterations $`j_\mu ^{A(1)},j_\mu ^{A(2)},\mathrm{}`$ mar2 , is defined. The requirement of gauge invariance leads to the completely nonlinear theory and principal impossibility of deriving closed kinetic equations for number densities (i.e. two-point functions) of the soft-fermion and -gluon excitations.

In this work, we restrict consideration only to a finite number of terms in expansions of the induced source and color currents. This imposes more rigid restrictions on the order of the magnitudes of soft fields. Blaizot and Iancu ian  (Sec. 1.1) showed in the special case when the soft fields were thermal fluctuations at the soft scale $`gT`$ (namely, this situation takes place in our case), their typical amplitudes would be of the orders

$$|A_\mu (X)|g^{1/2}T\left(|A_\mu (p)|\frac{1}{g^{1/2}(gT)^3}\right),$$

$$|\overline{\psi }(X)||\psi (X)|g^2T^3\left(|\overline{\psi }(q)||\psi (q)|\frac{1}{g(gT)^5}\right).$$

In this case the parameter of expansions $`\delta g^{1/2}1`$ and we have the estimates

$$\eta ^{(0,1)}(q)gT|\psi (q)|,\eta ^{(1,1)}(q)g^{3/2}T|\psi (q)|,\eta ^{(2,1)}(q)g^{5/2}T|\psi (q)|,\mathrm{},$$

$$j_\mu ^{A(1)}(p)\frac{1}{g^{3/2}T},j_\mu ^{A(2)}(p)\frac{1}{gT},j_\mu ^{A(3)}(p)\frac{1}{g^{1/2}T},\mathrm{},$$

$$j_\mu ^{\mathrm{\Psi }(0,2)}(p)g(gT)^4|\overline{\psi }(q)||\psi (q)|\left(\frac{1}{gT}\right),j_\mu ^{\mathrm{\Psi }(1,2)}(p)g^{3/2}(gT)^4|\overline{\psi }(q)||\psi (q)|\left(\frac{1}{g^{1/2}T}\right),$$

$$\mathrm{},$$

i.e., every successive term in the source and current expansions is suppressed by one power of $`g^{1/2}`$ in comparison with the preceding ones, and the use of the perturbation theory is therefore justified. By taking into account the nonlinear interaction between soft fields and hard particles at leading order in $`g`$, it is sufficient to keep only the first three terms in the expansions (2.13) and (2.17), and first two terms in expansion (2.14), as mentioned at the end of Sec. II. The recovery of a gauge symmetry here, occurs another way, by the account of weak correlation in the calculation of the expectation value of three soft fields in a slightly heterogeneous, slightly nonstationary quark-gluon plasma (see the following section). In this case we derive a closed system of the gauge-invariant kinetic equations that describe the evolution of the off-equilibrium deviations in the number densities, and obtain the damping rates from the collision terms, which are closely allied in form to the corresponding damping rates in the HTL approximation bra2 ; kob .

## VI Generalized kinetic equation for soft fermion excitations

Now we turn to the equation for the soft-fermion excitations (4.4). We substitute Eqs. (4.8) and (4.13) into third-order correlation functions entering into Eq. (4.4). Because of the fact that $`\overline{\psi }^{(0)},\psi ^{(0)}`$, and $`A_\mu ^{(0)}`$ represent the amplitudes of entirely uncorrelated waves, the correlation function $`A_\mu ^{(0)}\overline{\psi }^{(0)}\psi ^{(0)}`$ drops out. In this case, each term in $`\overline{\psi }^i(q)A^{a\mu }(p_1)(^{}\mathrm{\Gamma }_\mu ^{(Q)a}\psi (q_1))^j`$ and $`(\overline{\psi }(q_1)^{}\mathrm{\Gamma }_\mu ^{(Q)a})^iA^{a\mu }(p_1)\psi ^j(q^{})`$ should be defined more exactly. In the fourth-order correlation functions, within the accepted accuracy, we replace the fields $`\overline{\psi },\psi `$, and $`A_\mu `$ by $`\overline{\psi }^{(0)},\psi ^{(0)}`$, and $`A_\mu ^{(0)}`$.

Furthermore, we make the correlation decoupling of the fourth-order correlation functions in terms of the pair correlation functions by the rules

$$\overline{\psi }^{(0)}(q_1)\psi ^{(0)}(q_2)\overline{\psi }^{(0)}(q_3)\psi ^{(0)}(q_4)=\overline{\psi }^{(0)}(q_1)\psi ^{(0)}(q_2)\overline{\psi }^{(0)}(q_3)\psi ^{(0)}(q_4)$$

$$\overline{\psi }^{(0)}(q_1)\psi ^{(0)}(q_4)\overline{\psi }^{(0)}(q_3)\psi ^{(0)}(q_2),$$

$$A^{(0)}(p_1)A^{(0)}(p_2)\overline{\psi }^{(0)}(q_1)\psi ^{(0)}(q_2)=A^{(0)}(p_1)A^{(0)}(p_2)\overline{\psi }^{(0)}(q_1)\psi ^{(0)}(q_2).$$

Here, we suppress the color, spinor, and Lorentz indices, and take into consideration the Grassmanian nature of the quark field. Setting $`\mathrm{{\rm Y}}_{\beta \alpha }^{ij}=\delta ^{ji}\mathrm{{\rm Y}}_{\beta \alpha }`$ and $`I_{\mu \nu }^{ab}=\delta ^{ab}I_{\mu \nu }`$, taking into account $`\delta ^{ab}\delta \mathrm{\Gamma }_{\mu \nu }^{ab}(p_1,p_2;q_1,q_2)C_F\delta \mathrm{\Gamma }_{\mu \nu }(p_1,p_2;q_1,q_2)`$, after cumbersome calculations, we come to the following generalized kinetic equation333The term generalized here means that we have not yet restricted ourselves to any mass-shell conditions. for the soft-fermionic excitation QGP, instead of Eq.(4.4)

$$\mathrm{tr}\left(\frac{}{q^\mu }[\overline{)}q+\delta \mathrm{\Sigma }^\mathrm{H}(q)]\frac{\mathrm{{\rm Y}}(q,x)}{x_\mu }\right)=2i\mathrm{tr}(\delta \mathrm{\Sigma }^\mathrm{A}(q)\mathrm{{\rm Y}}(q,x))$$
(6.1)
$$+\mathrm{\hspace{0.17em}2}g^2C_F(dq_1\mathrm{Im}[^{}𝒟_{\mu \nu }(qq_1)\mathrm{tr}\left\{\mathrm{{\rm Y}}(q)^{}\mathrm{\Gamma }^{(Q)\mu }(qq_1;q_1,q)\mathrm{{\rm Y}}(q_1)^{}\mathrm{\Gamma }^{(G)\nu }(qq_1;q_1,q)\right\}]$$

$$dq_1dp_1I_{\mu \nu }(p_1)\mathrm{Im}\left[\mathrm{tr}\{^{}\mathrm{\Gamma }^{(Q)\mu }(p_1;q_1,q)^{}S(q)^{}\mathrm{\Gamma }^{(Q)\nu }(p_1;q_1,q)\mathrm{{\rm Y}}(q_1)\}\right]\delta (qq_1p_1)$$

$$dpI_{\mu \nu }(p)\mathrm{Im}[\mathrm{tr}\{(\delta \mathrm{\Gamma }^{(Q)\mu \nu }(p,p;q,q)^{}\mathrm{\Gamma }^{(Q)\mu }(p;qp,q)^{}S(qp)^{}\mathrm{\Gamma }^{(Q)\nu }(p;q,q+p)$$

$$^{}\mathrm{\Gamma }^{(Q)\nu }(p;q+p,q)^{}S(q+p)^{}\mathrm{\Gamma }^{(Q)\mu }(p;q,qp))\mathrm{{\rm Y}}(q)\}]).$$

On the right-hand side of Eq. (6.1) the $`x`$ dependence of $`\mathrm{{\rm Y}}_{\beta \alpha }(q)`$ and $`I_{\mu \nu }(p)`$ is understood, although not explicitly written.

In deriving Eq. (6.1) we use the following properties of effective two-quark–one-gluon and two-quark–two-gluon vertex functions, which immediately follows from initial definitions (4.5), (4.6), (4.7) and (4.12)

$$\gamma ^0{}_{}{}^{}\mathrm{\Gamma }_{\mu }^{(Q)}(p;q_1,q_2)\gamma ^0=^{}\mathrm{\Gamma }_\mu ^{(Q)}(p;q_2,q_1)=^{}\mathrm{\Gamma }_\mu ^{(Q)}(p;q_1,q_2),$$

$$\gamma ^0{}_{}{}^{}\mathrm{\Gamma }_{\mu }^{(G)}(p;q_1,q_2)\gamma ^0=^{}\mathrm{\Gamma }_\mu ^{(G)}(p;q_1,q_2)=^{}\mathrm{\Gamma }_\mu ^{(G)}(p;q_2,q_1),$$

$$\gamma ^0\delta \mathrm{\Gamma }_{\mu \nu }^{(Q)}(p_1,p_2;q_1,q_2)\gamma ^0=\delta \mathrm{\Gamma }_{\mu \nu }^{(Q)}(p_1,p_2;q_2,q_1)=\delta \mathrm{\Gamma }_{\mu \nu }^{(Q)}(p_1,p_2;q_1,q_2).$$

Moreover, we assume that under Hermitian conjugation the Wigner function $`\mathrm{{\rm Y}}(q,x)`$ behaves like an ordinary $`\gamma `$ matrix:

$$\gamma ^0\mathrm{{\rm Y}}^{}(q,x)\gamma ^0=\mathrm{{\rm Y}}(q,x).$$
(6.2)
Below this property is shown to ensure that the physical variable such as the fermion number density is real.

Equation (6.1) is incomplete, since the unknown function $`I_{\mu \nu }(p,x)`$ enters this equation. This function obeys the kinetic equation, which is similar to Eq. (6.1). The soft-field equation (2.20) is initial for derivation of this kinetic equation. We use the obtained expressions for the induced color current $`j_\mu ^\mathrm{\Psi }`$ (4.10) and (4.14) and corresponding expressions for the induced color current $`j_\mu ^A`$ defined in Ref. mar2 . Performing calculations similar to previous ones and keeping in the right-hand side only the terms responsible for the nonlinear interaction between soft-fermion excitations and soft-gluon ones (purely gauge sector was considered in detail in Ref. mar2 ), we come to the following generalized kinetic equation for soft-gluon excitation QGP’s (again, the $`x`$ dependence of most quantities is suppressed)

$$\frac{}{p_\lambda }[p^2g^{\mu \nu }(1+\xi ^1)p^\mu p^\nu +\delta \mathrm{\Pi }^{\mathrm{H}\mu \nu }(p)]\frac{I_{\mu \nu }(p,x)}{x^\lambda }=\mathrm{\hspace{0.17em}2}i\delta \mathrm{\Pi }^{\mathrm{A}\mu \nu }(p)I_{\mu \nu }(p,x)$$
(6.3)
$$\mathrm{\hspace{0.17em}2}g^2T_F(\mathrm{Im}(^{}𝒟_{\mu \nu }(p))dq_1dq_2\mathrm{tr}\left\{\mathrm{{\rm Y}}(q_1)^{}\mathrm{\Gamma }^{(G)\mu }(p;q_1,q_2)\mathrm{{\rm Y}}(q_2)^{}\mathrm{\Gamma }^{(G)\nu }(p;q_1,q_2)\right\}$$

$$+I_{\mu \nu }(p)dq\mathrm{Im}[\mathrm{tr}\{(\delta \mathrm{\Gamma }^{(G)\mu \nu }(p,p;q,q)+^{}\mathrm{\Gamma }^{(Q)\nu }(p;pq,q)^{}S(pq)^{}\mathrm{\Gamma }^{(G)\mu }(p;p+q,q)$$

$$^{}\mathrm{\Gamma }^{(G)\mu }(p;q,pq)^{}S(p+q)^{}\mathrm{\Gamma }^{(Q)\nu }(p;p+q,q))\mathrm{{\rm Y}}(q)\}]).$$

Here, $`T_F`$ is the index of the fundamental representation, $`\delta \mathrm{\Pi }^{\mathrm{H}\mu \nu }(p)`$ and $`\delta \mathrm{\Pi }^{\mathrm{A}\mu \nu }(p)`$ are the Hermitian and anti-Hermitian parts of the soft-gluon self-energy, respectively, and $`\delta \mathrm{\Gamma }^{(G)\mu \nu }`$ is the effective two-quark–two-gluon vertex function defined by an expression

$$\delta \mathrm{\Gamma }^{(G)\mu \nu }(p_1,p_2;q_1,q_2)=\omega _0^2\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{v^\mu v^\nu \overline{)}v}{(vq_1iϵ)(vq_2iϵ)}\left(\frac{1}{v(q_1+p_1)+iϵ}+\frac{1}{v(q_1+p_2)iϵ}\right),$$
(6.4)
where the first term in large parentheses satisfies the chronological order: $`X_1^0Y_2^0,X_1^0X_2^0Y_1^0`$ in the coordinate representation, and the second term satisfies the chronological order $`X_1^0Y_1^0,X_1^0X_2^0Y_2^0`$; i.e., the time argument of one of external gluon legs coming into the vertex is largest. The mean of the superscript $`(G)`$ is evident here. Notice that in deriving Eq. (6.3) we use the properties of Hermitian conjugate of the vertex function (6.4):

$$\gamma ^0\delta \mathrm{\Gamma }_{\mu \nu }^{(G)}(p_1,p_2;q_1,q_2)\gamma ^0=\delta \mathrm{\Gamma }_{\mu \nu }^{(G)}(p_2,p_1;q_1,q_2)=\delta \mathrm{\Gamma }_{\mu \nu }^{(G)}(p_1,p_2;q_2,q_1).$$

In closing this section we again call attention to the chronological order of arguments in the vertex functions entering into the right-hand side of generalized equations (6.1) and (6.3), which appear here in nonevident manner. Also, we call attention to the fact that the right-hand side of Eq. (6.3) is proportional to the third group invariant $`T_F`$ as distinct from the right-hand side of Eq. (6.1), proportional to the quark Casimir invariant $`C_F`$ and the kinetic equation for $`I_{\mu \nu }(p,x)`$ with purely gauge sector on the right-hand side mar2 , proportional to the gluon Casimir invariant $`C_A`$.

## VII Classification of types of nonlinear wave interactions. Plasmino and plasmon kinetic equations

Now we perform a preliminary analysis of the right-hand sides of the generalized kinetic equations (6.1) and (6.3), the purpose of which is the identification of the terms with specific physical processes described by them. To establish this connection it is convenient to use the Feynman diagrams, defining the contribution of the leading order in the coupling constant to the damping of soft-fermion excitations in framework of resummed perturbation theory. There are known bra1  to be two diagrams (Figs. 1(a) and 1(b))

with soft loop momentum that contribute at order $`g`$ to the effective self-energy for quarks: the usual self-energy graph at one-loop order (Fig. 1(a)), except that all of the vertices and propagators are effective, and the effective gluon tadpole on the soft-fermion line (Fig. 1(b)), which is special to the effective expansion.

Let us consider the second term444As was mentioned at the end of Sec. III, the first term on the right-hand side of Eq. (6.1), involving the imaginary part of the HTL self-energy, describes linear Landau damping. It vanishes, when one studies the decay of on-shell excitations and therefore can be omitted. on the right-hand side of Eq. (6.1). Physically, this contribution just corresponds to the stimulated scattering processes of soft-fermion excitations by hard-particle QGP through a soft-virtual gluon oscillation, without the change of statistic type of both soft and hard excitations, i.e.,

$$\mathrm{q}+\mathrm{Q}(\overline{\mathrm{Q}})\mathrm{q}+\mathrm{Q}(\overline{\mathrm{Q}}),\mathrm{q}+\mathrm{G}\mathrm{q}+\mathrm{G},$$

and similarly for $`\overline{\mathrm{q}}`$. Here, $`\mathrm{q}`$ represents fermionic collective excitations (we do not distinguish between normal-particle excitations and plasminos) and $`\mathrm{Q},\overline{\mathrm{Q}}`$, and $`\mathrm{G}`$ are quark, antiquark, and gluon excitations with characteristic momenta of order $`T`$. It can be easily inferred if cutting the usual self-energy graph before spelling out the effective gluon propagator, i.e., inserting a hard bubble along the gluon line, as depicted in Fig. 2.

The third term on the right-hand side of Eq. (6.1) is associated with the so-called decay processes depicted in Fig. 3 following from the cutting of the usual self-energy graph

$$\mathrm{q}\mathrm{q}+\mathrm{g},\mathrm{q}+\overline{\mathrm{q}}\mathrm{g},$$

where g is the gluonic collective excitations (here, we do not distinguish between longitudinal and transverse excitations). Similar decay processes have been studied in detail, since they are immediately connected with the processes of

soft-dilepton productions in hot QCD plasma. In the high-temperature QCD plasma in the approach based on using dispersion relation methods bra4 , these processes are connected with the so-called pole-pole terms.

Finally, the remaining terms grouped on the right-hand side of Eq. (6.1) are connected with the more interesting process of the nonlinear wave interaction with stimulated scattering of soft-fermion excitations by hard-particle QGP’s varying with the change of statistics of excitations:

$$\mathrm{q}(\overline{\mathrm{q}})+\mathrm{G}\mathrm{g}+\mathrm{Q}(\overline{\mathrm{Q}}),\mathrm{q}(\overline{\mathrm{q}})+\overline{\mathrm{Q}}(\mathrm{Q})\mathrm{g}+\mathrm{G}.$$

It is defined by two different physical processes. The first represents the Compton scattering type process and is connected with the effective two-quark–two-gluon vertex function (more exactly, with its imaginary part) on the right-hand side of Eq. (6.1). This corresponds diagrammatically to the appropriate cutting of the gluon “tadpole” graph, before illustrating the effective two-quark–two-gluon vertex function, as depicted in Fig. 4.

The second process connected with the remaining terms defines the scattering of a fermionic quantum oscillation through a virtual soft-fermion oscillation. Diagrammatically this process of scattering is defined by the corresponding cutting of the usual self-energy graph before illustrating the effective quark propagator as in Fig. 5.

As it will be shown in the next section, the terms corresponding to the processes drawn in Figs. 4 and 5 form a gauge-invariant function. The process, described by these terms, will be called further the nonlinear Landau damping of the soft-fermion excitations by virtue of the great similarity of the expression defining this process and the corresponding expression in a purely gauge case mar2 .

Here, it is necessary to note that the terms on the right-hand side of kinetic equation (6.1) corresponding to the process diagrammatically shown in Fig. 5 are very close in structure to the terms connected with stimulated scattering processes, without changing the statistic type and which are diagrammatically given in Fig. 2. However, in the last case the Compton scattering type contribution is absent, which, in turn, is associated with the absence of the effective quark tadpole on the soft-quark line. For this reason, we not assign this process to the one of the nonlinear Landau damping.

Now we turn to the right-hand side of generalized kinetic equation (6.3). The second term on the right-hand side of Eq. (6.3) is defined by time reversal of the decay processes that are drawn in Figs. 3(a) and 3(c). The remaining terms on the right-hand side of Eq. (6.3) are associated with the processes defined by the time reversal of the processes depicted in Figs. 4 and 5, i.e., with the process of the nonlinear Landau damping of the soft-gluon excitations.

In the remainder of this work, we restrict our attention to the detailed study of the process of nonlinear Landau damping, and moreover we restrict ourselves only to the process of the nonlinear interaction between plasminos and plasmons, i.e., purely collective excitations in hot QCD plasmas. The exception is only Sec. XII, where we consider the difficulties connected with the calculation of the nonlinear Landau damping rate for normal quark excitations in the vicinity of the light cone.

As it was mentioned in the Introduction in global equilibrium QGP’s for the case of massless quarks and for zero chemical potential there are two branches of fermion excitations with a positive energy. In this connection we define the Wigner function $`\mathrm{{\rm Y}}(q,x)`$ in the form of an expansion

$$\mathrm{{\rm Y}}(q,x)=h_+(\widehat{𝐪})\stackrel{~}{\mathrm{{\rm Y}}}_+(q,x)+h_{}(\widehat{𝐪})\stackrel{~}{\mathrm{{\rm Y}}}_{}(q,x),$$
(7.1)
where $`\stackrel{~}{\mathrm{{\rm Y}}}_\pm `$ are certain scalar functions and $`h_\pm (\widehat{𝐪})=(\gamma ^0\widehat{𝐪}\gamma )/2`$ with $`\widehat{𝐪}𝐪/|𝐪|.`$ By the condition (6.2) the functions $`\stackrel{~}{\mathrm{{\rm Y}}}_\pm `$ are real. Let us define dependence on variable $`q^0`$.

For this purpose we omit nonlinear terms and the anti-Hermitian part of the quark self-energy in Eq. (3.4). Substituting the function $`\delta ^{ji}(h_\pm (\widehat{𝐪}))_{\beta \alpha }\stackrel{~}{\mathrm{{\rm Y}}}_\pm (q,x)\delta (qq^{})`$ instead of $`\overline{\psi }_\alpha ^i(q)\psi _\beta ^j(q^{})`$, we find

$$\left[q^0\left(|𝐪|+\mathrm{Re}\delta \mathrm{\Sigma }_\pm (q)\right)\right]\stackrel{~}{\mathrm{{\rm Y}}}_\pm (q,x)=0,$$

where

$$\delta \mathrm{\Sigma }_\pm (q)=\pm \frac{1}{2}\mathrm{tr}\left\{h_\pm (\widehat{𝐪})\delta \mathrm{\Sigma }(q)\right\}=\frac{\omega _0^2}{|𝐪|}\left[\mathrm{\hspace{0.17em}1}\left(1\frac{|𝐪|}{q^0}\right)F\left(\frac{q^0}{|𝐪|}\right)\right]$$
(7.2)
with

$$F(z)=\frac{z}{2}\left[\mathrm{ln}\left|\frac{1+z}{1z}\right|i\pi \theta (1|z|)\right].$$

The solutions of these equations have the structure

$$\stackrel{~}{\mathrm{{\rm Y}}}_\pm (q,x)=\mathrm{{\rm Y}}_\pm (𝐪,x)\delta (q^0\omega _\pm (𝐪))+\mathrm{{\rm Y}}_{}(𝐪,x)\delta (q^0+\omega _{}(𝐪)).$$
(7.3)
Here, $`\mathrm{{\rm Y}}_\pm (𝐪,x)\mathrm{{\rm Y}}_𝐪^\pm `$ are certain functions of a wave vector $`𝐪`$ and $`\omega _\pm (𝐪)\omega _𝐪^\pm `$ are frequencies of the normal-particle excitations and plasminos, respectively. The solutions (7.3) describe what is called a quasiparticle approximation. In this approximation we assume that the off-equilibrium two-point functions $`\stackrel{~}{\mathrm{{\rm Y}}}_\pm (q,x)`$ have the same mass-shell conditions as in equilibrium, with an off-equilibrium deviation in the spectral densities $`\mathrm{{\rm Y}}_𝐪^\pm `$.

Furthermore, let us represent the function $`I_{\mu \nu }(p,x)=I_{\mu \nu }`$ in the form of an expansion

$$I_{\mu \nu }=Q_{\mu \nu }(p)I_p^l+P_{\mu \nu }(p)I_p^t,I_p^{l,t}I^{l,t}(p,x),$$
(7.4)
where

$$Q_{\mu \nu }(p)=\frac{\overline{u}_\mu (p)\overline{u}_\nu (p)}{\overline{u}^2(p)},P_{\mu \nu }(p)=g_{\mu \nu }D_{\mu \nu }(p)Q_{\mu \nu }(p),$$

$$D_{\mu \nu }=p_\mu p_\nu /p^2,\overline{u}_\mu (p)=p^2u_\mu p_\mu (pu),u_\mu =(1,0,0,0)$$

are longitudinal and transverse projections in the Lorentz covariant form, respectively, and we take the functions $`I_p^{l,t}`$ in the form of the quasiparticle approximation

$$I_p^l=I_𝐩^l\delta (p^0\omega _𝐩^l)+I_𝐩^l\delta (p^0+\omega _𝐩^l),I_p^t=I_𝐩^t\delta (p^0\omega _𝐩^t)+I_𝐩^t\delta (p^0+\omega _𝐩^t).$$
(7.5)
In Eq. (7.5) $`I_𝐩^{l,t}`$ are certain functions of a wave vector $`𝐩`$, and $`\omega _𝐩^{l,t}\omega ^{l,t}(𝐩)`$ are frequencies of the gluonic longitudinal and transverse eigenwaves in a QGP.

The equations describing the variation of the spectral densities $`\mathrm{{\rm Y}}_𝐪^{}`$ and $`I_𝐩^l`$ are obtained from Eqs. (6.1) and (6.3) by the replacements

$$\mathrm{{\rm Y}}(q,x)h_{}(\widehat{𝐪})\mathrm{{\rm Y}}_𝐪^{}\delta (q^0\omega _𝐪^{}),I_{\mu \nu }(p,x)Q_{\mu \nu }(p)[I_𝐩^l\delta (p^0\omega _𝐩^l)+I_𝐩^l\delta (p^0+\omega _𝐩^l)]$$
(7.6)
for Eq. (6.1), and

$$I_{\mu \nu }(p,x)Q_{\mu \nu }(p)I_𝐩^l\delta (p^0\omega _𝐩^l),\mathrm{{\rm Y}}(q,x)h_{}(\widehat{𝐪})\mathrm{{\rm Y}}_𝐪^{}\delta (q^0\omega _𝐪^{})+h_+(\widehat{𝐪})\mathrm{{\rm Y}}_𝐪^{}\delta (q^0+\omega _𝐪^{})$$

for Eq. (6.3). $`\delta `$ functions enable us to remove integration over $`q^0`$ and $`p^0`$. Retaining on the right-hand side the terms responsible for the process of nonlinear Landau damping only, we have instead of Eqs. (6.1) and (6.3)

$$2\left(\frac{}{q^\mu }[q^0+|𝐪|+\mathrm{Re}\delta \mathrm{\Sigma }_{}(q)]\right)_{q^0=\omega _𝐪^{}}\frac{\mathrm{{\rm Y}}_𝐪^{}}{x_\mu }=2g^2C_F\mathrm{{\rm Y}}_𝐪^{}d𝐩I_𝐩^l\mathrm{Im}T^{(Q)}(𝐪,𝐩),$$
(7.7)
where

$$T^{(Q)}(𝐪,𝐩)\mathrm{tr}\{(\delta \mathrm{\Gamma }^{(Q)\mu \nu }(p,p;q,q)^{}\mathrm{\Gamma }^{(Q)\mu }(p;qp,q)^{}S(qp)^{}\mathrm{\Gamma }^{(Q)\nu }(p;q,q+p)$$

$$^{}\mathrm{\Gamma }^{(Q)\nu }(p;q+p,q)^{}S(q+p)^{}\mathrm{\Gamma }^{(Q)\mu }(p;q,qp))h_{}(\widehat{𝐪})\}Q_{\mu \nu }(p)|_{q^0=\omega _𝐪^{},p^0=\omega _𝐩^l},$$
(7.8)
and

$$\left(\frac{}{p^\mu }[p^2+\mathrm{Re}\delta \mathrm{\Pi }^l(p)]\right)_{p^0=\omega _𝐩^l}\frac{I_𝐩^l}{x_\mu }=2g^2T_FI_𝐩^ld𝐪\mathrm{{\rm Y}}_𝐪^{}\mathrm{Im}T^{(G)}(𝐪,𝐩),$$
(7.9)
where

$$T^{(G)}(𝐪,𝐩)\mathrm{tr}\{(\delta \mathrm{\Gamma }^{(G)\mu \nu }(p,p;q,q)+^{}\mathrm{\Gamma }^{(Q)\nu }(p;pq,q)^{}S(pq)^{}\mathrm{\Gamma }^{(G)\mu }(p;p+q,q)$$

$$^{}\mathrm{\Gamma }^{(G)\mu }(p;q,pq)^{}S(p+q)^{}\mathrm{\Gamma }^{(Q)\nu }(p;p+q,q))h_{}(\widehat{𝐪})\}Q_{\mu \nu }(p)|_{q^0=\omega _𝐪^{},p^0=\omega _𝐩^l}.$$
(7.10)
The expressions $`T^{(Q)}(𝐪,𝐩)`$ and $`T^{(G)}(𝐪,𝐩)`$ contain the factors

$$\frac{1}{vq+iϵ},\frac{1}{vqiϵ},\frac{1}{v(qp)+iϵ},\mathrm{}.$$

The imaginary parts of the first two factors should be set equal to zero because they are connected with the linear Landau damping of a plasmino (which is absent in a QGP) and therefore, the imaginary part of the expressions (7.8) and (7.10) will be defined as

$$\mathrm{Im}\frac{1}{v(qp)+iϵ}|_{q^0=\omega _𝐪^{},p^0=\omega _𝐩^l}=\pi \delta (\omega _𝐪^{}\omega _𝐩^l𝐯(𝐪𝐩)),\mathrm{}.$$

Because of the absence of the linear Landau damping for the plasmino from the definitions of effective two-quark–two-gluon vertex functions (4.7) and (6.4) and two-quark–one-gluon vertex functions (4.6) and (4.12), the next equalities are as follows

$$\delta \mathrm{\Gamma }^{(Q)\mu \nu }(p,p;q,q)=\delta \mathrm{\Gamma }^{(G)\mu \nu }(p,p;q,q),^{}\mathrm{\Gamma }^{(Q)\mu }(p;qp,q)=^{}\mathrm{\Gamma }^{(G)\mu }(p;qp,q),\mathrm{},$$
(7.11)
and as consequence, we have

$$T^{(Q)}(𝐪,𝐩)=T^{(G)}(𝐪,𝐩).$$
(7.12)
Therefore, in the subsequent discussion the superscripts $`(Q)`$ and $`(G)`$ will be dropped. Let us emphasize that the equalities (7.11) and (7.12) hold only for on-shell plasma excitations.

## VIII Gauge invariance of Im $`T(𝐪,𝐩)`$

The issue of the gauge dependence of the plasmino nonlinear Landau damping rate is closely associated with the issue of the gauge dependence of $`\mathrm{Im}T(𝐪,𝐩)`$. In the next section the expression $`\mathrm{Im}T(𝐪,𝐩)`$ is shown to define the plasmino–hard-particle scattering matrix element squared, i.e., it has direct physical relevance. To establish the gauge invariance of the matrix element for plasmino–hard-particle scattering we need to show that the expression $`\mathrm{Im}T(𝐪,𝐩)`$ in the covariant gauge equals a similar one in the temporal gauge.

We prove the gauge invariance for a more general expression, which in a covariant gauge has the form

$${}_{}{}^{}\stackrel{~}{\mathrm{\Gamma }}(p_1,p_2;q_1,q_2)\{\delta \mathrm{\Gamma }^{\mu \nu }(p_1,p_2;q_1,q_2)^{}\mathrm{\Gamma }^\mu (p_1;q_1+p_2,q_2)^{}S(q_1+p_2)^{}\mathrm{\Gamma }^\nu (p_2;q_1,q_2+p_1)$$

$$^{}\mathrm{\Gamma }^\nu (p_2;q_1+p_1,q_2)^{}S(q_1+p_1)^{}\mathrm{\Gamma }^\mu (p_1;q_1,q_2+p_2)\}\overline{u}_\mu (p_1)\overline{u}_\nu (p_2)|_{\mathrm{on}\mathrm{shell}}.$$
(8.1)
Here, $`p_1+p_2+q_1+q_2=0`$. Notice that the order of the Lorentz indices of the last two terms in curly brackets on the right-hand side of Eq. (8.1) is important. The association of the expression (7.8) with (8.1) is given by

$$\mathrm{Im}T(𝐪,𝐩)=\frac{1}{\sqrt{\overline{u}^2(p)\overline{u}^2(p_1)}}\mathrm{Im}\mathrm{tr}\{^{}\stackrel{~}{\mathrm{\Gamma }}(p_1,p_2;q_1,q_2)h_{}(\widehat{𝐪}_1)\}|_{p_1=p_2=p,q_1=q_2=q}.$$

A similar expression (8.1) in the temporal gauge is obtained with the replacements

$$\overline{u}_\mu (p)\stackrel{~}{u}_\mu (p)p^2(p_\mu u_\mu (pu))/(pu),\text{etc}.$$

The gauge-invariance proof is based on using identities analogous to the effective Ward identities in hot gauge theory bra1 ; tay . It can be shown that the following equalities hold:

$$\delta \mathrm{\Gamma }^{\mu \nu }(p_1,p_2;q_1,q_2)p_{1\mu }=^{}\mathrm{\Gamma }^\nu (p_2;q_1,q_2+p_1)^{}\mathrm{\Gamma }^\nu (p_2;q_1+p_1,q_2),$$

$$\delta \mathrm{\Gamma }^{\mu \nu }(p_1,p_2;q_1,q_2)p_{2\nu }=^{}\mathrm{\Gamma }^\mu (p_1;q_1,q_2+p_2)^{}\mathrm{\Gamma }^\mu (p_1;q_1+p_2,q_2),$$
(8.2)
$${}_{}{}^{}\mathrm{\Gamma }_{}^{\mu }(p;q_1,q_2)p_\mu =^{}S^1(q_1)^{}S^1(q_1+p).$$

Initially we calculate the convolution with the effective two-quark–two-gluon vertex function $`\delta \mathrm{\Gamma }^{\mu \nu }`$. Slightly cumbersome but not complicated computations by using the effective Ward identities (8.2) and the mass-shell conditions lead to the following expression

$$\delta \mathrm{\Gamma }^{\mu \nu }(p_1,p_2;q_1,q_2)\overline{u}_\mu (p_1)\overline{u}_\nu (p_2)|_{\mathrm{onshell}}=p_1^2p_2^2\delta \mathrm{\Gamma }^{00}(p_1,p_2;q_1,q_2)+\mathrm{\Xi }(p_1,p_2;q_1,q_2),$$
(8.3)
where

$$\mathrm{\Xi }(p_1,p_2;q_1,q_2)=p_1^0p_2^2\{^{}\mathrm{\Gamma }^0(p_2;q_1+p_1,q_2)^{}\mathrm{\Gamma }^0(p_2;q_1,q_2+p_1)\}$$

$$+p_2^0p_1^2\{^{}\mathrm{\Gamma }^0(p_1;q_1+p_2,q_2)^{}\mathrm{\Gamma }^0(p_1;q_1,q_2+p_2)\}p_1^0p_2^0\{^{}S^1(q_1+p_2)+^{}S^1(q_1+p_1)\}.$$

Furthermore, we calculate the convolution with the terms containing the two-quark–one-gluon effective vertex in Eq. (8.1). Here, we derive the expression

$$\{^{}\mathrm{\Gamma }^\mu (p_1;q_1+p_2,q_2)^{}S(q_1+p_2)^{}\mathrm{\Gamma }^\nu (p_2;q_1,q_2+p_1)+(p_1p_2,\mu \nu )\}\overline{u}_\mu (p_1)\overline{u}_\nu (p_2)|_{\mathrm{on}\mathrm{shell}}$$
(8.4)
$$=p_1^2p_2^2\{^{}\mathrm{\Gamma }^0(p_1;q_1+p_2,q_2)^{}S(q_1+p_2)^{}\mathrm{\Gamma }^0(p_2;q_1,q_2+p_1)+(p_1p_2)\}+\mathrm{\Xi }(p_1,p_2;q_1,q_2).$$

Subtracting Eq. (8.4) from (8.3), we arrive at the desired expression

$${}_{}{}^{}\stackrel{~}{\mathrm{\Gamma }}(p_1,p_2;q_1,q_2)=p_1^2p_2^2\{\delta \mathrm{\Gamma }^{00}(p_1,p_2;q_1,q_2)^{}\mathrm{\Gamma }^0(p_1;q_1+p_2,q_2)^{}S(q_1+p_2)^{}\mathrm{\Gamma }^0(p_2;q_1,q_2+p_1)$$

$$^{}\mathrm{\Gamma }^0(p_2;q_1+p_1,q_2)^{}S(q_1+p_1)^{}\mathrm{\Gamma }^0(p_1;q_1,q_2+p_2)\}.$$
(8.5)
Now we consider the structure of $`{}_{}{}^{}\stackrel{~}{\mathrm{\Gamma }}`$ in the temporal gauge. For this purpose we replace $`\overline{u}_\mu `$ by $`\stackrel{~}{u}_\mu `$ in Eq. (8.1). The convolution with an effective two-quark–two-gluon vertex function leads to

$$\delta \mathrm{\Gamma }^{\mu \nu }(p_1,p_2;q_1,q_2)\stackrel{~}{u}_\mu (p_1)\stackrel{~}{u}_\nu (p_2)|_{\mathrm{onshell}}=p_1^2p_2^2\delta \mathrm{\Gamma }^{00}(p_1,p_2;q_1,q_2)+\stackrel{~}{\mathrm{\Xi }}(p_1,p_2;q_1,q_2),$$
(8.6)
where

$$\stackrel{~}{\mathrm{\Xi }}(p_1,p_2;q_1,q_2)=\frac{p_1^2p_2^2}{p_1^0p_2^0}\{\left(p_1^0{}_{}{}^{}\mathrm{\Gamma }_{}^{0}(p_1;q_1+p_2,q_2)^{}S^1(q_1+p_2)\right)+(p_1p_2)\}.$$

Convolution with the terms containing $`{}_{}{}^{}\mathrm{\Gamma }_{}^{\mu }`$ yields

$$\{^{}\mathrm{\Gamma }^\mu (p_1;q_1+p_2,q_2)^{}S(q_1+p_2)^{}\mathrm{\Gamma }^\nu (p_2;q_1,q_2+p_1)+(p_1p_2,\mu \nu )\}\stackrel{~}{u}_\mu (p_1)\stackrel{~}{u}_\nu (p_2)|_{\mathrm{on}\mathrm{shell}}$$
(8.7)
$$=p_1^2p_2^2\{^{}\mathrm{\Gamma }^0(p_1;q_1+p_2,q_2)^{}S(q_1+p_2)^{}\mathrm{\Gamma }^0(p_2;q_1,q_2+p_1)+(p_1p_2)\}+\stackrel{~}{\mathrm{\Xi }}(p_1,p_2;q_1,q_2).$$

Subtracting Eq. (8.7) from (8.6), we obtain a similar expression (8.5). Thus, we have shown that at least in the class of covariant and temporal gauges, the expression $`\mathrm{Im}T(𝐪,𝐩)`$ is gauge invariant.

## IX Nonlinear Landau damping rate of a plasmino

In our work mar2  the transformation of the nonlinear Landau damping rate for a plasmon in a purely gauge sector of soft QGP excitations was performed. An expression similar to $`\mathrm{Im}T(𝐪,𝐩)`$ was transformed to module squared of a sum of the terms, which are interpreted as the scattering amplitude of specific physical processes. In this section we represent a similar transformation for the function $`\mathrm{Im}T(𝐩,𝐪)`$. Since this function is not dependent on the choice of a gauge, we choose a temporal gauge for simplicity.

The right-hand side of Eq. (7.7) contains the contribution of two different processes. The first is associated with absorption of a plasmino by QGP particles with frequency $`\omega _𝐪^{}`$ and wave vector $`𝐪`$, and consequent radiation of a plasmon with frequency $`\omega _𝐩^l`$ and wave vector $`𝐩`$. It is defined by the second term on the right-hand side of Eq. (7.8). The frequencies and wave vectors of an incident plasmino and a recoil plasmon satisfy the conservation law

$$\omega _𝐪^{}\omega _𝐩^l𝐯(𝐪𝐩)=0.$$
(9.1)
The second process represents simultaneous radiation (or absorption) of a plasmino and plasmon with frequencies $`\omega _𝐪^{}`$, $`\omega _𝐩^l`$ and wave vectors $`𝐪`$, $`𝐩`$ satisfying the conservation law

$$\omega _𝐪^{}+\omega _𝐩^l𝐯(𝐪+𝐩)=0,$$
(9.2)
and is defined by the third term in Eq. (7.8). The contribution of the second process to the order of interest is not important. The first term on the right-hand side of Eq. (7.8), associated with an effective two-quark–two-gluon vertex function, contains both processes. Furthermore, we take into account the terms in Eq. (9.1) and we drop the terms that contain the $`\delta `$ function of Eq. (9.2). Notice that the conservation law (9.1) defines, in a very nontrivial manner, accessible kinematic regions of wave vectors $`𝐪`$ and $`𝐩`$ in the process of nonlinear Landau damping for a plasmino. In Sec. XI we consider only the simplest limiting case of accessible kinematic region of plasmon wave vector $`𝐩`$ for the case of the $`𝐪=0`$ mode.

Now we consider the expression with $`\delta \mathrm{\Gamma }^{\mu \nu }`$. With regard to the discussion above and to the definition (6.4), the contribution of a given term to $`\mathrm{Im}T(𝐪,𝐩)`$ can be represented as

$$2\pi \omega _0^2\left[\frac{p^2}{𝐩^2(p^0)^2}\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{(1+𝐯\widehat{𝐪})(𝐯𝐩)^2}{(vq)^2}\delta (\omega _𝐪^{}\omega _𝐩^l𝐯(𝐪𝐩))\right]_{\mathrm{on}\mathrm{shell}}.$$
(9.3)
Let us consider now a more complicated second term in Eq. (7.8) associated with the two-quark–one-gluon vertices. By the definition (4.5) and (4.6), the following equality is obeyed:

$${}_{}{}^{}\mathrm{\Gamma }_{}^{\mu }(p;l,q)=^{}\mathrm{\Gamma }^\mu (p;q,l)^{}\mathrm{\Gamma }^\mu $$

(hereafter, $`lqp`$). By using this relation and the decomposition of the effective quark propagator $`{}_{}{}^{}S(l)`$ onto $`h_\pm (\widehat{𝐥})`$,

$${}_{}{}^{}S(l)=h_+(\widehat{𝐥})^{}\mathrm{}_+(l)+h_{}(\widehat{𝐥})^{}\mathrm{}_{}(l),$$

where

$${}_{}{}^{}\mathrm{}_{\pm }^{}(l)=\frac{1}{l^0[|𝐥|+\delta \mathrm{\Sigma }_\pm (l)]},$$
(9.4)
the contribution to $`\mathrm{Im}T(𝐪,𝐩)`$ of the term associated with the effective two-quark–one-gluon vertex functions can be represented as

$$\left[\frac{p^2}{𝐩^2(p^0)^2}\mathrm{Im}(^{}\mathrm{}_+(l)\mathrm{tr}\left\{h_{}(\widehat{𝐪})(^{}\mathrm{\Gamma }^i\mathrm{p}^i)h_+(\widehat{𝐥})(^{}\mathrm{\Gamma }^j\mathrm{p}^j)\right\}+^{}\mathrm{}_{}(l)\mathrm{tr}\left\{h_{}(\widehat{𝐪})(^{}\mathrm{\Gamma }^i\mathrm{p}^i)h_{}(\widehat{𝐥})(^{}\mathrm{\Gamma }^j\mathrm{p}^j)\right\})\right]_{\mathrm{on}\mathrm{shell}}.$$
(9.5)
Let us compute the traces in the last expression. For this purpose we research the matrix structure of the function $`{}_{}{}^{}\mathrm{\Gamma }_{}^{i}\mathrm{p}^i`$ in more detail. From the analysis of the corresponding expressions for HTL functions, derived by Frenkel and Taylor (in particular, the expression (3.38) in Ref. fre ), it is easy to see that $`{}_{}{}^{}\mathrm{\Gamma }_{}^{i}\mathrm{p}^i`$ can be represented in the form of the expansion

$${}_{}{}^{}\mathrm{\Gamma }_{}^{i}\mathrm{p}^i=\gamma ^0\delta \Gamma _0+(𝐥\gamma )^{}\Gamma _{}+((𝐧\times 𝐥)\gamma )^{}\Gamma _{},$$
(9.6)
where $`𝐧𝐪\times 𝐩`$ and the coefficient functions are defined as

$$\delta \Gamma _0=\omega _0^2\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{𝐯𝐩}{(vl+iϵ)(vq)},$$

$${}_{}{}^{}\Gamma _{}^{}=\frac{𝐩𝐥}{𝐥^2}+\delta \Gamma _{}\frac{𝐩𝐥}{𝐥^2}\frac{\omega _0^2}{𝐥^2}\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{(𝐯𝐩)(𝐯𝐥)}{(vl+iϵ)(vq)},$$
(9.7)
$${}_{}{}^{}\Gamma _{}^{}=\frac{1}{𝐥^2}+\delta \Gamma _{}\frac{1}{𝐥^2}\frac{\omega _0^2}{𝐥^2𝐧^2}\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{(𝐯𝐩)(𝐯(𝐧\times 𝐥))}{(vl+iϵ)(vq)}.$$

The matrix basis in the expansion (9.6) is convenient in that it is “ortogonal” in trace computing. Substituting expression (9.6) into (9.5) we can compute desired traces in Eq. (9.5) in terms of the functions (9.7). However, this direct approach is not quite convenient in view of its nontransparency and the necessity of cumbersome calculations. We overdetermine the expansion (9.6), which enables us to represent the function $`{}_{}{}^{}\mathrm{\Gamma }_{}^{i}\mathrm{p}^i`$ in a more appropriate form for trace computing. Instead of expansion (9.6) we can write (the details of overdetermination are given in Appendix A )

$${}_{}{}^{}\mathrm{\Gamma }_{}^{i}\mathrm{p}^i=h_{}(\widehat{𝐥})^{}\Gamma _+h_+(\widehat{𝐥})^{}\Gamma _{}+2h_{}(\widehat{𝐪})𝐥^2|𝐪|^{}\Gamma _{},$$
(9.8)
where

$${}_{}{}^{}\Gamma _{\pm }^{}\delta \Gamma _0|𝐥|^{}\Gamma _{}+\frac{𝐧^2}{|𝐪|}\frac{1}{1\widehat{𝐪}\widehat{𝐥}}^{}\Gamma _{}.$$
(9.9)
The explicit selectivity of the matrix $`h_{}(\widehat{𝐪})`$ in the expansion (9.8) is connected with the fact that we have restricted our consideration to the study of a plasmino branch of fermion excitations only. In the case of the branch describing normal-particle excitations, it is necessary to use the following expansion instead of Eq. (9.8):

$${}_{}{}^{}\mathrm{\Gamma }_{}^{i}\mathrm{p}^i=h_{}(\widehat{𝐥})^{}\stackrel{´}{\Gamma }_+h_+(\widehat{𝐥})^{}\stackrel{´}{\Gamma }_{}2h_+(\widehat{𝐪})𝐥^2|𝐪|^{}\Gamma _{},$$
(9.10)
where now

$${}_{}{}^{}\stackrel{´}{\Gamma }_{\pm }^{}\delta \Gamma _0|𝐥|^{}\Gamma _{}\frac{𝐧^2}{|𝐪|}\frac{1}{1\pm \widehat{𝐪}\widehat{𝐥}}^{}\Gamma _{}.$$
(9.11)
Substituting expression (9.8) into (9.5), and using the identities for $`h_\pm (\widehat{𝐥})`$ matrices,

$$h_\pm (\widehat{𝐥})h_\pm (\widehat{𝐥})=0,h_\pm (\widehat{𝐥})h_{}(\widehat{𝐥})h_\pm (\widehat{𝐥})=h_\pm (\widehat{𝐥}),$$

after computing the trivial traces, we obtain

$$\left[\frac{p^2}{𝐩^2(p^0)^2}\{(1\widehat{𝐪}\widehat{𝐥})\mathrm{Im}(^{}\mathrm{}_+(l)(^{}\Gamma _+)^2)+(1+\widehat{𝐪}\widehat{𝐥})\mathrm{Im}(^{}\mathrm{}_{}(l)(^{}\Gamma _{})^2)\}\right]_{\mathrm{on}\mathrm{shell}}.$$
(9.12)
Furthermore, we use the relation

$$\mathrm{Im}(^{}\mathrm{}_\pm (l)(^{}\Gamma _\pm )^2)=\mathrm{Im}(^{}\mathrm{}_\pm ^1(l))|^{}\mathrm{}_\pm (l)^{}\Gamma _\pm |^2+2\mathrm{Im}(^{}\Gamma _\pm )\mathrm{Re}(^{}\mathrm{}_\pm (l)^{}\Gamma _\pm ).$$

Taking into consideration the last relation and the equalities

$$\mathrm{Im}^{}\mathrm{}_\pm ^1(l)=\pi \omega _0^2\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }(1𝐯\widehat{𝐥})\delta (vl),$$

$$\mathrm{Im}^{}\Gamma _\pm =\pi \omega _0^2\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\frac{𝐯𝐩}{vq}(1𝐯\widehat{𝐥}+\frac{𝐯(𝐧\times 𝐥)}{|𝐪|𝐥^2(1\widehat{𝐪}\widehat{𝐥})})\delta (vl),$$

subtracting Eq. (9.12) from (9.3), we define the desired expression for $`\mathrm{Im}T`$:

$$\mathrm{Im}T(𝐪,𝐩)=\pi \omega _0^2[\frac{p^2}{𝐩^2(p^0)^2}\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\{2\frac{(𝐯𝐩)^2}{(vq)^2}(1+𝐯\widehat{𝐪})$$

$$+(1\widehat{𝐪}\widehat{𝐥})(1𝐯\widehat{𝐥})|^{}\mathrm{}_+(l)^{}\Gamma _+|^2+2\frac{𝐯𝐩}{vq}\varrho _+(𝐯;\widehat{𝐪},\widehat{𝐩})\mathrm{Re}(^{}\mathrm{}_+(l)^{}\Gamma _+)$$
(9.13)
$$+(1+\widehat{𝐪}\widehat{𝐥})(1+𝐯\widehat{𝐥})|^{}\mathrm{}_{}(l)^{}\Gamma _{}|^2+2\frac{𝐯𝐩}{vq}\varrho _{}(𝐯;\widehat{𝐪},\widehat{𝐩})\mathrm{Re}(^{}\mathrm{}_{}(l)^{}\Gamma _{})\}\delta (vl)]_{\mathrm{on}\mathrm{shell}}.$$

Here,

$$\varrho _\pm (𝐯;\widehat{𝐪},\widehat{𝐩})(1\widehat{𝐪}\widehat{𝐥})(1𝐯\widehat{𝐥})+\frac{𝐯(𝐧\times 𝐥)}{|𝐪|𝐥^2}.$$
(9.14)
The expression (9.13) can be led to a more descriptive form. For this purpose we add to the expression in the curly brackets on the right-hand side of Eq. (9.13) the term

$$\frac{𝐯(𝐧\times 𝐥)}{|𝐪|𝐥^2}(|^{}\mathrm{}_+(l)^{}\Gamma _+|^2+|^{}\mathrm{}_{}(l)^{}\Gamma _{}|^2),$$

equals zero over solid integration. In addition, we rewrite the first term in curly brackets by using the relation

$$2(1+𝐯\widehat{𝐪})=\varrho _+(𝐯;\widehat{𝐪},\widehat{𝐩})+\varrho _{}(𝐯;\widehat{𝐪},\widehat{𝐩}),$$

which is a consequence of the definitions (9.14). Taking into account the discussion above, we rewrite expression (9.13) in the following form:

$$\mathrm{Im}T(𝐪,𝐩)=\pi \omega _0^2\frac{(\omega _𝐩^l)^2𝐩^2}{(\omega _𝐩^l)^2𝐩^2}\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\delta (\omega _𝐪^{}\omega _𝐩^l𝐯(𝐪𝐩))\{\varrho _+(𝐯;\widehat{𝐪},\widehat{𝐩})w_𝐯^+(𝐪,𝐩)+\varrho _{}(𝐯;\widehat{𝐪},\widehat{𝐩})w_𝐯^{}(𝐪,𝐩)\},$$
(9.15)
where

$$w_𝐯^\pm (𝐪,𝐩)=|_\pm (𝐪,𝐩)|^2\left|\frac{𝐯𝐩}{vq}+^{}\mathrm{}_\pm (l)^{}\Gamma _\pm \right|_{\mathrm{on}\mathrm{shell}}^2.$$
(9.16)
The denominator $`vq`$ in Eq. (9.16) is eikonal, which was expected in an approximation to the small-angle scattering of a high-energy particle. It defines the effective (here, effectiveness is not in terms of HTL resummation) propagator in the Compton type scattering process depicted in Fig. 4. The factor $`𝐯𝐩`$ is connected with the effective vertex with one soft external leg and two hard ones. Accordingly, the second term in the amplitudes $`_\pm (𝐪,𝐩)`$ represents the second type of scattering process depicted in Fig. 5. Here, the quark propagator transfers the soft momentum and is defined by the scalar functions $`{}_{}{}^{}\mathrm{}_{+}^{}`$ or $`{}_{}{}^{}\mathrm{}_{}^{}`$, and a vertex with soft external momenta is defined by the scalar functions $`{}_{}{}^{}\Gamma _{+}^{}`$ or $`{}_{}{}^{}\Gamma _{}^{}`$ for corresponding propagators.

Furthermore, we are relating the Wigner functions $`\mathrm{{\rm Y}}_𝐪^{}`$ and $`I_𝐩^l`$ to the plasmino and plasmon number densities, setting accordingly

$$n_𝐪^{}=2\mathrm{Z}_{}^1(𝐪)\mathrm{{\rm Y}}_𝐪^{},N_𝐩^l=\mathrm{\hspace{0.17em}2}\omega _𝐩^l\mathrm{Z}_l^1(𝐩)I_𝐩^l.$$
(9.17)
Here,

$$\mathrm{Z}_{}^1(𝐪)=1+\left(\frac{\mathrm{Re}\delta \mathrm{\Sigma }_{}(q)}{q^0}\right)_{q^0=\omega _𝐪^{}},\mathrm{Z}_l^1(𝐩)=1+\left(\frac{\mathrm{Re}\delta \mathrm{\Pi }^l(p)}{(p^0)^2}\right)_{p^0=\omega _𝐩^l}$$

are the residues of the effective quark and gluon propagators at the appropriate poles, respectively. The factor 2 in front of $`\mathrm{Z}_{}^1(𝐪)`$ takes into account the presence in QGP’s of antiplasminos. By using the definition (9.17), Eqs. (7.7) and (7.9) can be rewritten in the usual form of kinetic equations (containing on the left-hand side drift terms, and on the right-hand side the terms responsible for collisions), where the role distribution functions of quasiparticles–plasminos and plasmons–fulfill the functions (9.17)

$$\frac{n_𝐪^{}}{t}+𝐕_𝐪^{}\frac{n_𝐪^{}}{𝐱}=\gamma ^{}(𝐪)n_𝐪^{}g^2C_Fn_𝐪^{}d𝐩𝒬(𝐪,𝐩)N_𝐩^l,$$
(9.18)
$$\frac{N_𝐩^l}{t}+𝐕_𝐩^l\frac{N_𝐩^l}{𝐱}=+\gamma ^l(𝐩)N_𝐩^l+g^2T_FN_𝐩^ld𝐪𝒬(𝐪,𝐩)n_𝐪^{},$$
(9.19)
where $`𝐕_𝐪^{}\omega _𝐪^{}/𝐪`$ and $`𝐕_𝐩^l\omega _𝐩^l/𝐩`$ are group velocities of plasminos and plasmons respectively, and the kernel $`𝒬(𝐪,𝐩)`$ is defined by

$$𝒬(𝐪,𝐩)=\frac{1}{2\omega _𝐩^l}\mathrm{Z}_{}(𝐪)\mathrm{Z}_l(𝐩)\mathrm{Im}T(𝐪,𝐩).$$
(9.20)
The functions $`\gamma ^{}(𝐪)`$ and $`\gamma ^l(𝐩)`$ on the right-hand side of Eqs. (9.18) and (9.19) represent the nonlinear Landau damping rates for plasminos with momentum $`𝐪`$ and for plasmons with momentum $`𝐩`$, respectively.

The structure of kernel (9.20) is rather unexpected. As we see from the expression (9.15), this kernel is not reduced to the squared modulus of one scalar function, as this occurs in a purely gauge case mar2 . Here, $`𝒬(𝐪,𝐩)`$ is defined by a sum of the squared moduli of two independent scalar functions555As will be shown in the next section, the coefficient functions $`\varrho _\pm `$ in Eq. (9.15) are not, in general, case definite; therefore in principle it is impossible to express the expression $`\mathrm{Im}T`$ as the squared modulus of one scalar function–the total scattering amplitude.: amplitudes $`_+(𝐪,𝐩)`$ and $`_{}(𝐪,𝐩)`$. This is a point that we find difficult to interpret and therefore additional analysis of this problem is required. The only remark that may be made is that the interference between the scattering processes depicted in Fig. 5 proceeding through intermediate quark virtual oscillations with propagators $`{}_{}{}^{}\mathrm{}_{+}^{}(l)`$ and $`{}_{}{}^{}\mathrm{}_{}^{}(l)`$ accordingly, vanishes by the relation

$$h_{}(\widehat{𝐥})\left[\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\overline{)}v\delta (\omega _𝐪^{}\omega _𝐩^l𝐯(𝐪𝐩))\right]h_+(\widehat{𝐥})=0.$$

This relation is the analog of the relation in a purely gauge case mar2

$$(𝐩_2)^i\left[\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }v^iv^j\delta (\omega _𝐩^l\omega _{𝐩_1}^l𝐯(𝐩𝐩_1))\right](\stackrel{ˇ}{𝐧}\times 𝐩_2)^j=0,\stackrel{ˇ}{𝐧}𝐩\times 𝐩_1,𝐩_2𝐩𝐩_1,$$

that is responsible for the absence of interference between the scattering of plasmon by a QGP thermal particle through the longitudinal and transverse virtual gluon oscillations.

## X Decomposition of kernel $`𝒬(𝐪,𝐩)`$ into positive and negative parts

In this section we consider the problem of the direction of the effective pumping over of plasma excitation energy in the process of the nonlinear interaction of plasminos and plasmons. For this purpose, first, we study in more detail the structure of the $`\varrho _\pm (𝐯;\widehat{𝐪},\widehat{𝐩})`$ function.

By using the expansion (A4), we represent the functions $`\varrho _\pm `$ (9.14) in the form

$$\varrho _\pm (𝐯;\widehat{𝐪},\widehat{𝐩})=1+𝐯\widehat{𝐪}(\widehat{𝐪}\widehat{𝐥}+𝐯\widehat{𝐥}).$$
(10.1)
Let us introduce the coordinate system in which axis $`0Z`$ is aligned with vector $`\widehat{𝐥}`$; then the coordinates of vectors $`\widehat{𝐪}`$ and $`𝐯`$ are equal to $`\widehat{𝐪}=(1,\alpha ,\beta )`$ and $`𝐯=(1,\theta ,\phi ),`$ respectively. By $`\mathrm{\Phi }`$ we denote the angle between $`𝐯`$ and $`\widehat{𝐪}`$: $`𝐯\widehat{𝐪}=\mathrm{cos}\mathrm{\Phi }`$. The angle $`\mathrm{\Phi }`$ can be expressed as

$$\mathrm{cos}\mathrm{\Phi }=\mathrm{sin}\theta \mathrm{sin}\alpha \mathrm{cos}(\phi \beta )+\mathrm{cos}\theta \mathrm{cos}\alpha .$$
(10.2)
In the fixed coordinate system expression (10.1) reads

$$\varrho _\pm (\alpha ,\beta ;\theta ,\phi )=1+\mathrm{sin}\theta \mathrm{sin}\alpha \mathrm{cos}(\phi \beta )+\mathrm{cos}\theta \mathrm{cos}\alpha (\mathrm{cos}\alpha +\mathrm{cos}\theta ).$$
(10.3)
By using trigonometric formulas, it is easily to show that the last expression can be represented in the homogeneous quadratic form with respect to the variables $`\mathrm{cos}\frac{1}{2}(\theta \alpha )`$ and $`\mathrm{cos}\frac{1}{2}(\theta +\alpha )`$

$$\varrho _\pm (\alpha ,\beta ;\theta ,\phi )=𝐲^{}𝒜_\pm 𝐲,$$
(10.4)
where

$$𝐲𝐲(\alpha ,\theta )\left(\begin{array}{c}\mathrm{cos}\frac{1}{2}(\theta \alpha )\\ \mathrm{cos}\frac{1}{2}(\theta +\alpha )\end{array}\right),𝒜_\pm 𝒜_\pm (\beta ,\phi )\left(\begin{array}{cc}2\mathrm{cos}^2\frac{1}{2}(\phi \beta )& 1\\ 1& 2\mathrm{sin}^2\frac{1}{2}(\phi \beta )\end{array}\right).$$

The eigenvalues of matrices $`𝒜_\pm `$ are equal to

$$\lambda ^{(+)}=1+\sqrt{\mathrm{\hspace{0.17em}1}+\mathrm{cos}^2(\phi \beta )}>0,\lambda ^{()}=1\sqrt{\mathrm{\hspace{0.17em}1}+\mathrm{cos}^2(\phi \beta )}0.$$
(10.5)
From the theory of matrices gan  it is known that the real symmetric quadratic form is indefinite if and only if the eigenvalues $`\lambda `$ of matrix $`𝒜`$ have the different signs. Thus, by virtue of Eq. (10.5) we have proved that the coefficient functions $`\varrho _\pm `$ are indefinite. By the linear transformation of the variables $`\mathrm{cos}\frac{1}{2}(\theta \alpha )`$ and $`\mathrm{cos}\frac{1}{2}(\theta +\alpha )`$, the quadratic form (10.4) can be expressed as the canonical one, where the positive and negative parts of $`\varrho _\pm `$ are explicitly displayed

$$\varrho _\pm (\alpha ,\beta ;\theta ,\phi )=\lambda ^{(+)}(\chi _\pm ^{(+)})^2+\lambda ^{()}(\chi _\pm ^{()})^2.$$
(10.6)
Here, $`\chi _\pm ^{(+)}`$ and $`\chi _\pm ^{()}`$ are certain linear functions of $`\mathrm{cos}\frac{1}{2}(\theta \alpha )`$ and $`\mathrm{cos}\frac{1}{2}(\theta +\alpha )`$.

The explicit expressions for $`\chi _\pm ^{(+)}`$ and $`\chi _\pm ^{()}`$ are defined with the help of the transformations

$$\left(\begin{array}{c}\chi _\pm ^{(+)}\\ \chi _\pm ^{()}\end{array}\right)=𝒰_\pm \left(\begin{array}{c}\mathrm{cos}\frac{1}{2}(\theta \alpha )\\ \mathrm{cos}\frac{1}{2}(\theta +\alpha )\end{array}\right),$$
(10.7)
where, in turn, the transformation matrices $`𝒰_\pm `$ are defined by the solutions of the corresponding matrix equations

$$𝒰_\pm 𝒜_\pm =\mathrm{diag}(\lambda ^{(+)},\lambda ^{()})𝒰_\pm .$$

The solutions of the last equations have the following structure

$$𝒰_\pm =\left(\begin{array}{cc}\mathrm{cos}\vartheta & \pm \mathrm{sin}\vartheta \\ \mathrm{sin}\vartheta & \mathrm{cos}\vartheta \end{array}\right),\mathrm{tan}\vartheta =\frac{1}{\mathrm{cos}(\phi \beta )+\sqrt{\mathrm{\hspace{0.17em}1}+\mathrm{cos}^2(\phi \beta )}}.$$

By using a given expression and Eq. (10.7), we derive the desired expressions for $`\chi _\pm ^{(+)}`$ and $`\chi _\pm ^{()}`$. The functions $`\chi _\pm ^{(+)}`$ and $`\chi _\pm ^{()}`$ can be easily rewritten in the terms of initial vectors $`𝐥`$, $`𝐪`$ and $`𝐯`$. For example, from Eq. (10.2) it follows that

$$\mathrm{cos}(\phi \beta )=\frac{𝐧(𝐯\times 𝐥)}{|𝐧||𝐯\times 𝐥|}\text{etc}.$$

By substituting the expansion (10.6) into (9.15), we can represent the kernel $`𝒬(𝐪,𝐩)`$ in the form of a sum of positive and negative parts

$$𝒬(𝐪,𝐩)=𝒬^{(+)}(𝐪,𝐩)𝒬^{()}(𝐪,𝐩),𝒬^{(\pm )}(𝐪,𝐩)0,$$
(10.8)
where

$$𝒬^{(\pm )}(𝐪,𝐩)=\pm \pi \omega _0^2\left[\frac{(\omega _𝐩^l)^2𝐩^2}{2(\omega _𝐩^l)^3𝐩^2}\mathrm{Z}_{}(𝐪)\mathrm{Z}_l(𝐩)\right]$$
(10.9)
$$\times \frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\delta (\omega _𝐪^{}\omega _𝐩^l𝐯(𝐪𝐩))\{\lambda ^{(\pm )}[(\chi _+^{(\pm )})^2w_𝐯^+(𝐪,𝐩)+(\chi _{}^{(\pm )})^2w_𝐯^{}(𝐪,𝐩)]\}.$$

We use the decomposition (10.8) in general analysis of the problem of a direction of the pumping over of excitation energy in the study of the process of nonlinear interaction between plasminos and plasmons.

To consider the model problem of interaction of two infinitely narrow packets with typical wave vectors $`𝐪_0`$ and $`𝐩_0`$, let us introduce the number densities $`n_𝐪^{}`$ and $`N_𝐩^l`$ as follows:

$$n_𝐪^{}(t)=n^{}(t)\delta (𝐪𝐪_0),N_𝐩^l(t)=N^l(t)\delta (𝐩𝐩_0).$$

We have restricted ourselves to the spatially homogeneous case. Substituting the last expressions into Eqs. (9.18) and (9.19), we obtain the coupled nonlinear equations:

$$\frac{n^{}}{t}=g^2C_F𝒬n^{}N^l,n^{}(t_0)=n_0^{},$$
(10.10)
$$\frac{N^l}{t}=+g^2T_F𝒬n^{}N^l,N^l(t_0)=N_0^l.$$

Here, $`𝒬𝒬(𝐪_0,𝐩_0).`$ The system of the equations (10.10) possesses the integral of motion

$$𝒞T_Fn^{}(t)+C_FN^l(t)=T_Fn_0^{}+C_FN_0^l.$$

The general solution of this system, with regard to decomposition of the kernel (10.8), has the form

$$n^{}(t)=n_0^{}𝒞\frac{\mathrm{e}^{g^2𝒞𝒬^{(+)}(tt_0)}}{C_FN_0^l\mathrm{e}^{g^2𝒞𝒬^{()}(tt_0)}+T_Fn_0^{}\mathrm{e}^{g^2𝒞𝒬^{(+)}(tt_0)}},$$

$$N^l(t)=N_0^l𝒞\frac{\mathrm{e}^{g^2𝒞𝒬^{()}(tt_0)}}{C_FN_0^l\mathrm{e}^{g^2𝒞𝒬^{()}(tt_0)}+T_Fn_0^{}\mathrm{e}^{g^2𝒞𝒬^{(+)}(tt_0)}}.$$

Let us analyze the behavior of these solutions in the limit for $`t\mathrm{}`$.

1. Let us assume that the values of the wave vectors $`𝐪_0`$ and $`𝐩_0`$ are such that the following inequality

$$𝒬^{(+)}(𝐪_0,𝐩_0)>𝒬^{()}(𝐪_0,𝐩_0)$$
(10.11)
is true. Then in the limit for $`t\mathrm{}`$ we have

$$n^{}(t)0,N^l(t)N_0^l+(T_F/C_F)n_0^{}.$$
(10.12)
Thus we see that as a result of the nonlinear interaction of two infinitely narrow packets with fermion and boson quantum numbers, the effective pumping over of energy from the first packet to the second one takes place.
2. If the inverse inequality holds,

$$𝒬^{(+)}(𝐪_0,𝐩_0)<𝒬^{()}(𝐪_0,𝐩_0),$$
(10.13)
then in the limit for $`t\mathrm{}`$ we have

$$n^{}(t)n_0^{}+(C_F/T_F)N_0^l,N^l(t)0.$$

Here, we have an inverse case–the damping of bosonic excitations and the growth of fermionic excitations.
3. Finally, in the limiting case $`𝒬^{(+)}(𝐪_0,𝐩_0)=𝒬^{()}(𝐪_0,𝐩_0)`$, the interaction between wave packets is absent.

The inequalities (10.11) and (10.13) define the kinematic relations between wave vectors of excitations with different statistics, such that one or another process of pumping over of energy occurs. However, the general analysis of these inequalities is a complicated problem by virtue of the complexity of the expressions for positive and negative parts of the kernel $`𝒬`$. In the next section we shall study in detail only the simplest limiting case of the interaction of “standing” plasmino $`(𝐪=0)`$ with a plasmon. Here, the direction of the effective pumping over energy is easily defined.

## XI The damping rate of a plasmino at rest

The solution of the kinetic equation (9.18), defining a change of the plasmino number density, can be formally represented in the form (for a spatially homogeneous case)

$$n_𝐪^{()}(t)=n_𝐪^{()}(t_0)\mathrm{exp}\left\{g^2C_F\underset{t_0}{\overset{t}{}}dt^{}d𝐩𝒬(𝐪,𝐩)N_𝐩^l(t^{})\right\}.$$

Let us approximate the plasmon number density $`N_𝐩^l`$ by its equilibrium value or Planck distribution

$$N_𝐩^lN_{eq}^l(|𝐩|)=\frac{1}{(2\pi )^3}\frac{1}{\mathrm{e}^{\omega _𝐩^l/T}1}$$

and define the quasiparticle damping rate of the standing plasmino by means of the relation

$$\gamma ^{}(0)=\frac{1}{2}g^2C_F\underset{|𝐪|0}{lim}d𝐩𝒬(𝐪,𝐩)N_{eq}^l(|𝐩|).$$
(11.1)
Now we represent a complete calculation of $`\gamma ^{}(0)`$. We start with a representation of a kernel $`𝒬(𝐪,𝐩)`$ in the form (9.20), where the function $`\mathrm{Im}T(𝐪,𝐩)`$ is determined through Eqs. (9.15) and (9.16). Here, it is convenient for us to introduce a new coordinate system in which axis $`0Z`$ is aligned with the vector $`\widehat{𝐩}`$; then the coordinates of vectors $`\widehat{𝐪}`$ and $`𝐯`$ are equal to $`\widehat{𝐪}=(1,\alpha ^{},\beta ^{})`$ and $`𝐯=(1,\theta ^{},\phi ^{})`$, respectively. By $`\mathrm{\Phi }^{}`$ we denote the angle between $`𝐯`$, and $`\widehat{𝐪}`$. The angle $`\mathrm{\Phi }^{}`$ is connected with $`\alpha ^{}`$, $`\beta ^{}`$, $`\theta ^{}`$ and $`\phi ^{}`$ via the relation (10.2) with a corresponding change $`\alpha \alpha ^{}`$, etc.

In the limit for $`|𝐪|0`$, the kernel $`𝒬(𝐪,𝐩)`$ is reduced to

$$𝒬(0,𝐩)=\pi \omega _0^2\frac{(\omega _𝐩^l)^2𝐩^2}{8(\omega _𝐩^l)^3|𝐩|^3}\mathrm{Z}_l(𝐩)\frac{d\mathrm{\Omega }}{4\pi }\delta (\mathrm{cos}\theta ^{}\rho _𝐩^l)\{\varrho _+(𝐯;\widehat{𝐪},\widehat{𝐩})|_+(0,𝐩)|^2+\varrho _{}(𝐯;\widehat{𝐪},\widehat{𝐩})|_{}(0,𝐩)|^2\},$$
(11.2)
where $`\rho _𝐩^l(\omega _𝐩^l\omega _0)/|𝐩|`$; $`d\mathrm{\Omega }=\mathrm{sin}\theta ^{}\mathrm{d}\theta ^{}\mathrm{d}\phi ^{}.`$ The expressions for $`\varrho _\pm (𝐯;\widehat{𝐪},\widehat{𝐩})`$ are defined by Eq. (10.1). In a new coordinate system, instead of Eq. (10.3), we have the function representation $`\varrho _\pm `$ in terms of angles

$$\varrho _\pm (\alpha ^{},\beta ^{};\theta ^{},\phi ^{})=1+\mathrm{sin}\theta ^{}\mathrm{sin}\alpha ^{}\mathrm{cos}(\phi ^{}\beta ^{})+\mathrm{cos}\theta ^{}\mathrm{cos}\alpha ^{}\pm (\mathrm{cos}\alpha ^{}+\mathrm{cos}\theta ^{}).$$
(11.3)
The limits of the scattering amplitudes $`_\pm `$, by its definitions (9.16) equals

$$_\pm (0,𝐩)=\frac{|𝐩|}{\omega _0}\mathrm{cos}\theta ^{}+\underset{|𝐪|0}{lim}(^{}\mathrm{}_\pm (l)^{}\Gamma _\pm ).$$
(11.4)
Thus, the problem of computing of $`\gamma ^{}(0)`$ is reduced to the calculation of limits of scalar functions: the effective propagators $`{}_{}{}^{}\mathrm{}_{\pm }^{}`$ and the effective vertices $`{}_{}{}^{}\Gamma _{\pm }^{}`$.

By using the definitions of the $`{}_{}{}^{}\Gamma _{\pm }^{}`$ effective scalar vertex functions (9.9) and (9.7), after slightly cumbersome computations, we define

$$\underset{|𝐪|0}{\overset{}{lim}}\Gamma _\pm =\frac{|𝐩|}{\omega _0}\rho _𝐩^l\underset{|𝐪|0}{\overset{}{lim}}\mathrm{}_\pm ^1(l)[\omega _0|𝐩|\frac{𝐩^2}{\omega _0}\rho _𝐩^l(1\pm \rho _𝐩^l)].$$
(11.5)
In the last equality we use the $`{}_{}{}^{}\mathrm{}_{\pm }^{1}(q)`$ definition (9.4) as the $`F(q_0/|𝐪|)`$ function (7.2) (see Braaten and Pisarski bra2 ). Inserting Eq. (11.5) into (11.4), we reduce the scattering amplitudes to

$$_\pm (0,𝐩)=\frac{|𝐩|}{\omega _0}\mathrm{cos}\theta ^{}\frac{|𝐩|}{\omega _0}\rho _𝐩^l[\omega _0|𝐩|\frac{𝐩^2}{\omega _0}\rho _𝐩^l(1\pm \rho _𝐩^l)]\underset{|𝐪|0}{\overset{}{lim}}\mathrm{}_\pm (l).$$

By the $`\delta `$ function in Eq. (11.2), all terms in $`_\pm (0,𝐩)`$ not containing the propagators $`{}_{}{}^{}\mathrm{}_{\pm }^{}`$, are relatively reduced in the limit of the $`𝐪=0`$ mode.

The remaining terms, after substitution into Eq. (11.2) and integration over the solid angle with regard to Eq. (11.3), yield

$$𝒬(0,𝐩)=\pi \omega _0^2\frac{(\omega _𝐩^l)^2𝐩^2}{16(\omega _𝐩^l)^3|𝐩|^3}\mathrm{Z}_l(𝐩)\theta (1|\rho _𝐩^l|)\{(1+\mathrm{cos}\alpha ^{})(1+\rho _𝐩^l)[\omega _0|𝐩|\frac{𝐩^2}{\omega _0}\rho _𝐩^l(1+\rho _𝐩^l)]^2\left|\underset{|𝐪|0}{\overset{}{lim}}\mathrm{}_+(l)\right|^2$$

$$+(1\mathrm{cos}\alpha ^{})(1\rho _𝐩^l)[\omega _0+|𝐩|+\frac{𝐩^2}{\omega _0}\rho _𝐩^l(1\rho _𝐩^l)]^2|\underset{|𝐪|0}{\overset{}{lim}}\mathrm{}_{}(l)\left|^2\right\}.$$
(11.6)
We note that this expression is not dependent on the angle $`\beta ^{}`$. This enables us to represent the integration measure on the right-hand side of Eq. (11.1) in the form

$$d𝐩=2\pi \underset{0}{\overset{\mathrm{}}{}}𝐩^2\mathrm{d}|𝐩|\underset{1}{\overset{1}{}}\mathrm{d}(\mathrm{cos}\alpha ^{}).$$

Substituting Eq. (11.6) into (11.1) and performing an angular integration over $`\alpha ^{}`$, we finally obtain $`\gamma ^{}(0)`$

$$\gamma ^{}(0)=g^2C_F\underset{0}{\overset{\mathrm{}}{}}\mathrm{d}|𝐩|\theta (1|\rho _𝐩^l|)𝒬^{}(|𝐩|)N_{eq}^l(|𝐩|),$$
(11.7)
where the kernel $`𝒬^{}(|𝐩|)`$ has the form

$$𝒬^{}(|𝐩|)=\pi ^2\omega _0^2\frac{(\omega _𝐩^l)^2𝐩^2}{4(\omega _𝐩^l)^3|𝐩|}\mathrm{Z}_l(𝐩)\{(1+\rho _𝐩^l)[\omega _0|𝐩|\frac{𝐩^2}{\omega _0}\rho _𝐩^l(1+\rho _𝐩^l)]^2\left|\underset{|𝐪|0}{\overset{}{lim}}\mathrm{}_+(l)\right|^2$$
(11.8)
$$+(1\rho _𝐩^l)[\omega _0+|𝐩|+\frac{𝐩^2}{\omega _0}\rho _𝐩^l(1\rho _𝐩^l)]^2|\underset{|𝐪|0}{\overset{}{lim}}\mathrm{}_{}(l)\left|^2\right\}.$$

By virtue of the $`\theta `$ function in integrand (11.7), the kernel $`𝒬^{}(|𝐩|)`$ is positive; therefore

$$\gamma ^{}(0)>0,$$

i.e., the standing plasmino is damped. Thus, in the process of nonlinear Landau damping, the pumping over of the excitation energy of standing plasminos into plasmon branch of plasma excitations occurs, and therefore the first case (10.12) is true.

The function $`\rho _𝐩^l`$, entering the argument of $`\theta `$ function and in a kernel (11.8), decreases with momentum $`|𝐩|`$ from $`+\mathrm{}`$ reaching a minimum, and then monotonically increases at large $`|𝐩|`$, asymptotically tending to +1 from below. The equation $`\rho _𝐩^l=1`$ defines a lower limit of integration over the momentum of the recoil plasmon. The numerical solution of a given equation yields

$$|𝐩^{}|0.447\omega _{\mathrm{pl}}\text{for}N_f=2,|𝐩^{}|0.495\omega _{\mathrm{pl}}\text{for}N_f=3$$

(hereafter, the numerical estimates are presented for the $`\mathrm{SU}(3)`$ color group). The function $`\rho _𝐩^l`$ reaches an absolute minimum at

$$|𝐩^{}|1.236\omega _{\mathrm{pl}},\rho _𝐩^{}^l0.664\text{for}N_f=2,|𝐩^{}|1.305\omega _{\mathrm{pl}},\rho _𝐩^{}^l0.692\text{for}N_f=3.$$

Notice that the value of the function $`\rho _𝐩^l`$ at a minimum point coincides with the value of the plasmon group velocity $`\mathrm{d}\omega _𝐩^l/\mathrm{d}|𝐩|`$. Thus, if the incident plasmino is initially at rest, both the recoil plasmon and the effective quark that propagates below the light cone carry nonzero energies and momenta of order $`gT`$.

At the end of this section we compare derived expressions for the damping rate (11.7) and (11.8) of standing plasmino with the similar expression obtained in the framework of the HTL approximation bra2 ; kob . For this purpose in Eqs. (11.7) and (11.8) we rescale $`\omega _𝐩^l\omega _0\stackrel{~}{\omega }_𝐩^l`$ and $`|𝐩|\omega _0|\stackrel{~}{𝐩}|`$. The Planck distribution $`N_{\mathrm{eq}}^l(|𝐩|)`$ in the integrand in Eq. (11.7) should be set equal $`T/((2\pi )^3\omega _0\stackrel{~}{\omega }_𝐩^l)`$, since the energy is soft. Further rewriting the kernel (11.8) in terms of the functions accepted in the paper kob  and defining the damping constant as

$$\gamma ^{}(0)=\stackrel{ˇ}{a}(N_c,N_f)\frac{g^2TC_F}{4\pi },$$

we derive the expression for the coefficient $`\stackrel{ˇ}{a}(N_c,N_f)`$ required for comparison

$$\stackrel{ˇ}{a}(N_c,N_f)=\frac{1}{2}\underset{0}{\overset{\mathrm{}}{}}\frac{|\stackrel{~}{𝐩}|^2\mathrm{d}|\stackrel{~}{𝐩}|}{\stackrel{~}{\omega }_𝐩^l}\stackrel{~}{R}_l(\stackrel{~}{\omega }_𝐩^l,|\stackrel{~}{𝐩}|)\theta (|\stackrel{~}{𝐩}||\stackrel{~}{\omega }_𝐩^l1|)$$

$$\times \{[|\stackrel{~}{𝐩}|+\stackrel{~}{\omega }_𝐩^l2]^2\stackrel{~}{\beta }_+(1\stackrel{~}{\omega }_𝐩^l,|\stackrel{~}{𝐩}|)+[|\stackrel{~}{𝐩}|\stackrel{~}{\omega }_𝐩^l+2]^2\stackrel{~}{\beta }_{}(1\stackrel{~}{\omega }_𝐩^l,|\stackrel{~}{𝐩}|)\},$$

where

$$\stackrel{~}{R}_l(\stackrel{~}{\mathrm{p}}_0,|\stackrel{~}{𝐩}|)=\frac{\stackrel{~}{\mathrm{p}}_0}{\stackrel{~}{𝐩}^2}\frac{\stackrel{~}{\mathrm{p}}_0^2\stackrel{~}{𝐩}^2}{3r\stackrel{~}{\mathrm{p}}_0^2+\stackrel{~}{𝐩}^2},r=\frac{\omega _{pl}^2}{\omega _0^2}.$$

For the definitions of more cumbersome functions $`\stackrel{~}{\beta }_\pm (\stackrel{~}{\mathrm{p}}_0,|\stackrel{~}{𝐩}|)`$ see Ref. kob . The coefficient $`\stackrel{ˇ}{a}(N_c,N_f)`$ exactly coincides with the corresponding part in a similar coefficient $`a(N_c,N_f)`$ in the expression of the damping rate of the standing plasmino, responsible for the scattering process of plasminos by hard-particle QGP’s varying with the change of excitation statistic, derived by Kobes, Kunstatter, and Mak (kob , Eq. (5.11)). We have shown in doing so, by really remaining in the context of the Blaizot-Iancu equations (2.1)–(2.9), we are able to compute not only a spectrum of fermionic excitations in QGP’s, but also their (gauge-invariant) damping rate at the leading order in the coupling constant, which corresponds to the damping rate from the resummed perturbation theory bra2 ; kob .

## XII Light-cone singularity: Improved Blaizot-Iancu equations

In this section we discuss the behavior of the nonlinear Landau damping rate for normal-particle excitations, when the quark spectrum $`\omega _𝐪^+`$ approaches the light cone. Here, the initial expression is

$$\gamma ^+(𝐪)=\frac{1}{2}g^2C_Fd𝐩\stackrel{´}{𝒬}(𝐪,𝐩)N_{eq}^l(𝐩),\stackrel{´}{𝒬}(𝐪,𝐩)=\frac{1}{2\omega _𝐩^l}\mathrm{Z}_+(𝐪)\mathrm{Z}_l(𝐩)\mathrm{Im}\stackrel{´}{T}(𝐪,𝐩),$$
(12.1)
where

$$\mathrm{Im}\stackrel{´}{T}(𝐪,𝐩)=\pi \omega _0^2\frac{(\omega _𝐩^l)^2𝐩^2}{(\omega _𝐩^l)^2𝐩^2}\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\delta (\omega _𝐪^+\omega _𝐩^l𝐯(𝐪𝐩))\{\stackrel{´}{\varrho }_+(𝐯;\widehat{𝐪},\widehat{𝐩})\stackrel{´}{w}_𝐯^+(𝐪,𝐩)+\stackrel{´}{\varrho }_{}(𝐯;\widehat{𝐪},\widehat{𝐩})\stackrel{´}{w}_𝐯^{}(𝐪,𝐩)\},$$
(12.2)
and

$$\stackrel{´}{\varrho }_\pm (𝐯;\widehat{𝐪},\widehat{𝐩})(1\pm \widehat{𝐪}\widehat{𝐥})(1𝐯\widehat{𝐥})\frac{𝐯(𝐧\times 𝐥)}{|𝐪|𝐥^2},$$
(12.3)
$$\stackrel{´}{w}_𝐯^\pm (𝐪,𝐩)=|\stackrel{´}{}_\pm (𝐪,𝐩)|^2\left|\frac{𝐯𝐩}{vq}+^{}\mathrm{}_\pm (l)^{}\stackrel{´}{\Gamma }_\pm \right|_{\mathrm{onshell}}^2.$$
(12.4)
The scalar vertex functions $`{}_{}{}^{}\stackrel{´}{\Gamma }_{\pm }^{}`$ are defined by Eq. (9.11). On the right-hand side of Eq. (12.1) we have taken into account for simplicity only the contribution from the nonlinear interaction with longitudinal bosonic excitations, but all subsequent derivations easily extend to the case of the nonlinear interaction with transverse excitations.

Near the light cone, the HTL-vertex pieces $`\delta \Gamma _0`$, $`\delta \Gamma _{}`$, and $`\delta \Gamma _{}`$, the explicit forms of which are given by Eq. (9.7), are the origin of the strongest singularity of order $`1/\epsilon `$, where $`\epsilon ^2((\omega _𝐪^+)^2𝐪^2)/𝐪^2`$. The terms linear in $`\delta \Gamma `$ in the $`\stackrel{´}{w}_𝐯^\pm `$ probabilities (12.4) can lead to the logarithm of $`\epsilon `$ only fle2  (the first term in the amplitudes $`\stackrel{´}{}_\pm `$, defining Compton-like scattering processes, also results only in logarithmic divergence). Hence we restrict our consideration to the terms quadratic in $`\delta \Gamma `$.

Now let us single out in the HTL-vertex pieces the terms generating near-light-cone $`1/\epsilon `$-contributions to $`\gamma ^+(𝐪)`$. As an example we consider the scalar function $`\delta \Gamma _0`$. By using the explicit expression for HTL amplitudes derived by Frenkel and Taylor fre , we find

$$\delta \Gamma _0=\omega _0^2\frac{d\mathrm{\Omega }}{4\pi }\frac{𝐯𝐩}{(vl+iϵ)(vq)}\omega _0^2\frac{d\mathrm{\Omega }}{4\pi }\frac{𝒫}{vl}\frac{𝐯𝐩}{vq}$$

$$\frac{1}{𝐧^2}[l^0(𝐧\times 𝐪)𝐩q^0(𝐧\times 𝐥)𝐩]\mathrm{M}(l,q)=\omega _0^2p^0\mathrm{M}(l,q).$$

In the above expression, $`𝒫`$ indicates the principal-value prescription; here and in the following $``$ indicates that we have dropped less-divergent terms. The Lorentz-invariant function $`\mathrm{M}(l,q)`$ is defined by

$$\mathrm{M}(l,q)=\{\begin{array}{cc}\frac{1}{2\sqrt{\mathrm{\Delta }(l,q)}}\mathrm{ln}\left(\frac{lq+\sqrt{\mathrm{\Delta }(l,q)}}{lq\sqrt{\mathrm{\Delta }(l,q)}}\right),\hfill & \mathrm{\Delta }(l,q)l^2q^2(lq)^2<0\hfill \\ \frac{1}{\sqrt{\mathrm{\Delta }(l,q)}}\mathrm{arctan}\left(\frac{\sqrt{\mathrm{\Delta }(l,q)}}{lq}\right),\hfill & \mathrm{\Delta }(l,q)>0.\hfill \end{array}$$
(12.5)
Similarly, the terms containing the function $`\mathrm{M}(l,q)`$ in the HTL-vertex pieces $`\delta \Gamma _{}`$ and $`\delta \Gamma _{}`$ are singled out. Substituting next thus derived expressions into Eq. (9.11), we obtain

$${}_{}{}^{}\stackrel{´}{\Gamma }_{\pm }^{}\omega _0^2[(q^0l^0)\pm (q^0+l^0)\frac{q^0|𝐪|l^0|𝐥|}{|𝐪||𝐥|\pm 𝐪𝐥}]\mathrm{M}(l,q)\omega _0^2f_\pm (|𝐥|,|𝐪|,\widehat{𝐥}\widehat{𝐪})\mathrm{M}(l,q).$$

Furthermore, substituting the last expression into the $`\stackrel{´}{w}_𝐯^\pm `$ probabilities (12.4) and keeping only the terms quadratic in $`\mathrm{M}(l,q)`$, we have

$$\stackrel{´}{w}_𝐯^\pm (𝐪,𝐩)\omega _0^4|^{}\mathrm{}_\pm (l)|^2f_\pm ^2|\mathrm{M}(l,q)|^2,f_\pm f_\pm (|𝐥|,|𝐪|,\widehat{𝐥}\widehat{𝐪}).$$

The dependence on the unit vector $`𝐯`$ in the integrand of the expression (12.2) is defined only by the $`\delta `$ function and $`\stackrel{´}{\varrho }_\pm `$ function (12.3), that enables easily to perform the angular integration in Eq. (12.2) and thus to derive instead of Eq. (12.2)

$$\mathrm{Im}\stackrel{´}{T}(𝐪,𝐩)\pi \omega _0^6\frac{(\omega _𝐩^l)^2𝐩^2}{(\omega _𝐩^l)^2𝐩^2}\frac{1}{2|𝐥|}\theta (l^2)$$
(12.6)
$$\times \{(1\frac{l^0}{|𝐥|})(1+\widehat{𝐪}\widehat{𝐥})f_+^2|^{}\mathrm{}_+(l)|^2+(1+\frac{l^0}{|𝐥|})(1\widehat{𝐪}\widehat{𝐥})f_{}^2|^{}\mathrm{}_{}(l)|^2\}|\mathrm{M}(l,q)|^2.$$

The function $`\mathrm{Im}\stackrel{´}{T}(𝐪,𝐩)`$ is different from zero for $`l^2=(qp)^2<0`$. It is not difficult to show that the condition $`l^2<0`$ on the mass shell of plasma excitations leads to $`\mathrm{\Delta }(l,q)<0`$, and therefore in Eq. (12.5) it is necessary to take the first expression.

With $`q^20`$, the square of $`\mathrm{M}(l,q)`$ becomes $`1/\epsilon `$ times a representation of the $`\delta `$ function

$$|\mathrm{M}(l,q)|^2\stackrel{q^20}{}\frac{1}{\epsilon }\frac{\pi ^3}{|𝐪|\sqrt{l^2}}\delta (ql).$$
(12.7)
Substituting further the expression (12.6) into Eq. (12.1) and taking into account Eq. (12.7) and $`\mathrm{Z}_+(𝐪)\stackrel{q^20}{}1`$, we find the most singular contribution to the nonlinear Landau damping rate of the normal quark excitations near the light cone:

$$\gamma ^+(𝐪)\frac{1}{\epsilon }\frac{\pi ^3\omega _0^6}{2|𝐪|}g^2C_F\frac{\mathrm{d}\mathrm{\Omega }_{\widehat{𝐩}}}{4\pi }\underset{0}{\overset{\mathrm{}}{}}\mathrm{d}|𝐩|\frac{(\omega _𝐩^l)^2𝐩^2}{(\omega _𝐩^l)^3|𝐥|\sqrt{l^2}}\mathrm{Z}_l(𝐩)N_{eq}^l(𝐩)$$
(12.8)
$$\times \{(1\frac{l^0}{|𝐥|})(1+\widehat{𝐪}\widehat{𝐥})f_+^2|^{}\mathrm{}_+(l)|^2+(1+\frac{l^0}{|𝐥|})(1\widehat{𝐪}\widehat{𝐥})f_{}^2|^{}\mathrm{}_{}(l)|^2\}\theta (l^2)\delta (ql).$$

Here, the solid integral is over the directions of the unit vector $`\widehat{𝐩}`$. The $`\delta `$ function in integrand (12.8) enables us in principle to perform the angular integration. Because of the fact that the functions $`\mathrm{Z}_l(𝐩)`$ and $`N_{\mathrm{eq}}^l(𝐩)`$ vanish exponentially at large $`|𝐩|`$, the above expression is ultraviolet convergent. Notice also that by virtue of the $`\theta `$ function in the integrand, the damping rate (12.8) is positive. Thus we have shown that really the nonlinear Landau damping rate of the normal quark mode near the light cone diverges as $`1/\sqrt{q^2}`$, and thus signals the need to have further improvement of the Blaizot-Iancu equations, beyond the HTL approximation.

As mentioned in the Introduction, the light-cone singularities are associated with the massless basic constituents of a hot plasma, in our case, with massless hard quarks and hard transverse gluons. The inclusion of the asymptotic thermal masses for the basic constituents removes this type of singularity. The effective way of entering thermal masses into the resummed perturbation theory without spoiling gauge invariance was suggested by Flechsig and Rebhan fle3 . Below we use this approach, reformulating it into a “kinetic” language.

The light-cone singularity is generated by the HTL piece of the effective two-quark–one-gluon vertex function (4.6) (or (4.12)). In turn, the expression for $`\delta \mathrm{\Gamma }(p;q_1,q_2)`$ is defined with the help of the Blaizot-Iancu equation (2.8) for the function $`\delta n_\pm ^\mathrm{\Psi }`$ or by application of Eq. (2.9) for function $`\overline{)}\mathrm{\Lambda }`$. To derive the expression for $`\delta \mathrm{\Gamma }(p;q_1,q_2)`$, which is free from the light-cone singularity, we use improved Blaizot-Iancu equations instead of Eqs. (2.8) and (2.9) (the details of deriving these equations are given in Appendix B)

$$[vD_X,\delta n_\pm ^\mathrm{\Psi }(𝐤,X)]=\pm \frac{ig}{2ϵ_k}t^a(\overline{\psi }(X)t^a\overline{)}\mathrm{\Lambda }^\pm (𝐤,X)\overline{\overline{)}\mathrm{\Lambda }^\pm }(𝐤,X)t^a\psi (X)),$$
$`(2.8^{})`$
$$(vD_X)\overline{)}\mathrm{\Lambda }^\pm (𝐤,X)\pm \frac{\mathrm{}m_{\mathrm{}}^2}{2iϵ_k}\overline{)}\mathrm{\Lambda }^\pm (𝐤,X)=igC_F[N(ϵ_k)+n(ϵ_k)]\overline{)}v\psi (X),$$
$`(2.9^{})`$
where $`\overline{\overline{)}\mathrm{\Lambda }^\pm }(𝐤,X)=(\overline{)}\mathrm{\Lambda }^\pm (𝐤,X))^{}\gamma ^0`$, $`\mathrm{}m_{\mathrm{}}^2=m_{q,\mathrm{}}^2m_{g,\mathrm{}}^2`$, and $`m_{q,\mathrm{}}`$ and $`m_{g,\mathrm{}}`$ are asymptotic thermal masses for hard quarks and gluons, respectively. Instead of the expression (2.3) for the induced source now we have

$$\eta (X)=g\frac{\mathrm{d}𝐤}{(2\pi )^3}\frac{1}{2ϵ_k}\overline{)}\mathrm{\Lambda }^+(𝐤,X)+g\frac{\mathrm{d}𝐤}{(2\pi )^3}\frac{1}{2ϵ_k}\overline{)}\mathrm{\Lambda }^{}(𝐤,X).$$
$`(2.3^{})`$
The expression for the color current induced by the soft-quark field (2.5), is not changed.

Employing the improved Blaizot-Iancu equations $`(2.8^{})`$ and $`(2.9^{})`$, and expressions for induced source $`(2.3^{})`$ and the color current (2.5), it is not difficult to obtain the improved soft-quark self-energy and the effective two-quark–one-gluon vertex function used above. These expressions coincide with similar ones derived in the context of thermal field theory fle3 . Thus, for example, instead of Eq. (4.6) now we have

$$\delta \mathrm{\Gamma }_\mu ^{(Q)}(p;q_1,q_2)=\mathrm{\hspace{0.17em}2}\omega _0^2\frac{\mathrm{d}\mathrm{\Omega }}{4\pi }\{\frac{v_\mu \overline{)}v}{\left(vq_1+{\displaystyle \frac{\mathrm{}m_{\mathrm{}}^2}{2T\alpha }}+iϵ\right)\left(vq_2{\displaystyle \frac{\mathrm{}m_{\mathrm{}}^2}{2T\alpha }}iϵ\right)}$$

$$+\frac{v_\mu \overline{)}v}{\left(vq_1{\displaystyle \frac{\mathrm{}m_{\mathrm{}}^2}{2T\alpha }}+iϵ\right)\left(vq_2+{\displaystyle \frac{\mathrm{}m_{\mathrm{}}^2}{2T\alpha }}iϵ\right)}\}_\alpha .$$
(12.9)
Here,

$$𝒪(\alpha )_\alpha \frac{2}{\pi ^2}𝑑\alpha \frac{\alpha \mathrm{e}^\alpha }{\mathrm{e}^{2\alpha }1}𝒪(\alpha ),\alpha ϵ_k/T.$$

The use of expression (12.9) gives a finite nonlinear Landau damping rate $`\gamma ^+(𝐪)`$ near the light cone. The additional terms in denominators of Eq. (12.9) lead to the replacement of the divergence factor $`1/\epsilon `$ on the right-hand side of Eq. (12.8):

$$\frac{1}{\epsilon }\alpha ^{1/2}_\alpha \frac{1}{2}\left(\frac{T|𝐪|}{\mathrm{}m_{\mathrm{}}^2}\right)^{1/2},$$

where $`\alpha ^{1/2}_\alpha =(3/4\pi ^{3/2})(12^{2/3})\zeta (5/2)`$ and $`\zeta `$ is a Riemann zeta function.

At the end of this section notice that we run into similar divergence in research of a behavior near the light cone of the nonlinear Landau damping rate for the transverse bosonic excitations in purely gauge sector mar2 . The light cone singularity here is generated by the HTL piece of the effective three-gluon vertex function. It is clear that in this case for elimination of this singularity it is necessary to use the improved Blaizot-Iancu equations for functions $`\delta n_\pm ^A`$ and $`\delta N^A`$ instead of Eqs. (2.6) and (2.7). Here, we have not presented their explicit forms, noting only that their derivation is rather different from deriving improved equations $`(2.8^{})`$ and $`(2.9^{})`$.

## XIII Conclusion

In Sec. X the nonlinear interaction of plasminos with plasmons is shown to lead to the effective pumping over of plasma excitation energy from the fermionic branch of plasma oscillations to the bosonic branch and vice versa, and the kinematic relations between wave vectors of excitations are defined, such that one or another process of pumping over of energy occurs. However, it is clear that conclusions made in Sec. X are somewhat restricted since they were made with allowance only for one process–the process of nonlinear Landau damping. For precise study of nonlinear dynamics of soft-fermion and soft-boson excitations in hot QCD plasma it is necessary also to take into account in the plasmino and plasmon kinetic equations the remaining terms on the right-hand sides of generalized kinetic equations (6.1) and (6.3), responsible for the decay processes etc, and consider the existence of soft normal quark excitations and soft transverse gluon ones (the kinetic equations for them are defined from Eqs. (6.1) and (6.3) by the replacements similar to Eq. (7.6)). In addition to the right-hand side of the generalized kinetic equation (6.3) should be supplemented with the terms from purely gauge sector mar2 .

Thus obtained self-consistent nonlinear system of four kinetic equations will contain the maximum comprehensive information on soft excitations dynamics in hot QCD plasma, which may be obtained in the context of the initial Blaizot-Iancu equations (2.1)–(2.9) (or its improved variant) in the first nonvanishing approximation666The following term in the expansion over this small parameter leads to the processes of Boltzmann type:
$$\mathrm{q}+\mathrm{q}(\overline{\mathrm{q}})\mathrm{q}+\mathrm{q}(\overline{\mathrm{q}}),\mathrm{q}+\overline{\mathrm{q}}\mathrm{g}+\mathrm{g},\mathrm{g}+\mathrm{g}\mathrm{g}+\mathrm{g},\text{etc.},$$
going without exchange of energy between hard thermal particles and soft plasma waves. The example of the construction of the Boltzmann equation, describing the elastic scattering process of colorless plasmons between themselves can be found in Ref. mar3 . in the expansion over the small parameter $`gRA_\mu `$, where $`R(gT)^1`$ and $`|A_\mu |T`$. Certain nonlinear processes will be kinematically forbidden and therefore relevant contributions on the right-hand side of kinetic equations will drop out. However, even with allowance for the last case, this system of equations is especially complicated for any analytical research (similar to studies carried out in Sec. X of the present work and in Sec. 10 in Ref. mar2 ), and therefore here invoking numerical methods is required.

## Acknowledgment

This work was supported by an INTAS grant (No. 2000-15) and the Grant of the 6th Competition for Young Scientist of RAS (No. 1999-80).

## Appendix A

By using the relations $`\gamma ^0=h_{}(\widehat{𝐥})+h_+(\widehat{𝐥})`$ and $`𝐥\gamma =h_{}(\widehat{𝐥})h_+(\widehat{𝐥})`$, we rewrite the expansion (9.6) in the form

$${}_{}{}^{}\mathrm{\Gamma }_{}^{i}\mathrm{p}^i=h_{}(\widehat{𝐥})(\delta \Gamma _0+|𝐥|^{}\Gamma _{})+h_+(\widehat{𝐥})(\delta \Gamma _0|𝐥|^{}\Gamma _{})+((𝐧\times 𝐥)\gamma )^{}\Gamma _{}.$$
$`(\mathrm{A1})`$
Let us add to the expression in the first parentheses on the right-hand side of Eq. (A1) and then subtract from it the function $`𝐧^2/(|𝐪|(1\widehat{𝐪}\widehat{𝐥}))`$ and in the second parentheses add and subtract the function $`𝐧^2/(|𝐪|(1+\widehat{𝐪}\widehat{𝐥}))`$ correspondingly. Let us introduce the following functions:

$${}_{}{}^{}\Gamma _{\pm }^{}\delta \Gamma _0|𝐥|^{}\Gamma _{}+\frac{𝐧^2}{|𝐪|}\frac{1}{1\widehat{𝐪}\widehat{𝐥}}^{}\Gamma _{}.$$

Instead of Eq. (A1), now we have

$${}_{}{}^{}\mathrm{\Gamma }_{}^{i}\mathrm{p}^i=h_{}(\widehat{𝐥})^{}\Gamma _+h_+(\widehat{𝐥})^{}\Gamma _{}+\left[\frac{𝐧^2}{|𝐪|}\left(\frac{1}{1\widehat{𝐪}\widehat{𝐥}}h_{}(\widehat{𝐥})+\frac{1}{1+\widehat{𝐪}\widehat{𝐥}}h_+(\widehat{𝐥})\right)+((𝐧\times 𝐥)\gamma )\right]^{}\Gamma _{}.$$
$`(\mathrm{A2})`$
Furthermore, we transform the expression in the parentheses on the right-hand side of Eq. (A2). By using the definitions of matrices $`h_\pm (\widehat{𝐥})`$, this expression can be rewritten in the form

$$\frac{1}{2}\left[\gamma ^0\left(\frac{1}{1\widehat{𝐪}\widehat{𝐥}}+\frac{1}{1+\widehat{𝐪}\widehat{𝐥}}\right)+(\widehat{𝐥}\gamma )\left(\frac{1}{1\widehat{𝐪}\widehat{𝐥}}\frac{1}{1+\widehat{𝐪}\widehat{𝐥}}\right)\right]$$

$$=\frac{𝐪^2𝐥^2}{𝐧^2}(\gamma ^0+(\widehat{𝐥}\gamma )(\widehat{𝐪}\widehat{𝐥}))=2\frac{𝐪^2𝐥^2}{𝐧^2}h_{}(\widehat{𝐪})\frac{|𝐪|}{𝐧^2}((𝐧\times 𝐥)\gamma ).$$
$`(\mathrm{A3})`$
In deriving the last equality we have used the expansion of vector $`𝐪`$ into two mutually ortogonal vectors, $`𝐥`$ and $`𝐧\times 𝐥`$:

$$𝐪=\frac{𝐪𝐥}{𝐥^2}𝐥+\frac{1}{𝐥^2}(𝐧\times 𝐥).$$
$`(\mathrm{A4})`$
Substituting Eq. (A3) into (A2) and reducing the terms with vector product, finally we derive the expression (9.8) for convolution $`{}_{}{}^{}\mathrm{\Gamma }_{}^{i}\mathrm{p}^i`$, instead of Eq. (9.6).

## Appendix B

Here we show how one can obtain the improved Blaizot-Iancu equations, used in Sec. XII for deriving the finite nonlinear Landau damping rate for normal-particle excitations near the light cone. As the initial equations for obtaining an improved equation for a one-body density matrix mixing fermion and boson degrees of freedom, we take Eqs. (3.33) and (3.34) in the original paper of Blaizot and Iancu bla2  (see the accepted notations and definitions therein)

$$[_s^2+_s_X+2igA(X)_s]K_a^{\mu <}(s,X)=ig(\overline{)}_s\mathrm{\Delta })\gamma ^\mu t^a\psi (X),$$
$`(\mathrm{B1})`$
$$[\delta ^{ab}g_{\mu \nu }(_s^2_s_X)gf^{abc}\mathrm{\Gamma }_{\mu \nu \rho \lambda }A_c^\rho (X)_s^\lambda ]K_b^{\nu <}(s,X)=gS_0^<(s)\gamma _\mu t^a\psi (X).$$
$`(\mathrm{B2})`$
Furthermore, we expand $`K^<=K^{(0)}+K^{(1)}+\mathrm{}`$, with $`K^{(1)}gK^{(0)}`$, etc. The dominant terms in Eqs. (B1) and (B2) lead to the consistency condition

$$_s^2K^{(0)}(s,X)=0.$$
$`(\mathrm{B3})`$
Equation (B3) can be considered as the Klein-Gordon equation for hard, massless, free particles. However, due to interaction with the hot medium, the basic constituents of the plasma are known to acquire the dynamical (asymptotic) temperature-induced masses. To account for this fact we modify the condition (B3), replacing it in the case of Eq. (B1) by

$$_s^2K^{(0)}(s,X)m_{q,\mathrm{}}^2K^{(0)}(s,X),$$
$`(\mathrm{B4})`$
where $`m_{q,\mathrm{}}^2=g^2T^2C_F/4`$ is the asymptotic thermal mass for hard quark, and in the case of Eq. (B2) by

$$_s^2K^{(0)}(s,X)m_{g,\mathrm{}}^2K^{(0)}(s,X),$$
$`(\mathrm{B5})`$
where $`m_{g,\mathrm{}}^2=g^2T^2(N_f+2N_c)/12`$ is the asymptotic thermal mass for the hard transverse gluon. It is necessary to keep the terms $`_s^2`$ inside the brackets on the left-hand sides of Eqs. (B1) and (B2) setting them equals the right-hand side of Eqs. (B4) and (B5), accordingly. Further, subtracting Eq. (B2) from (B1) and performing transformations similar to those in Ref. bla2 , we come to improved equation for a density $`\overline{)}𝒦(k,X)=t^a\gamma ^\mu K_\mu ^a(k,X)`$:

$$(vD_X)\overline{)}𝒦(k,X)\frac{i}{2}\mathrm{}m_{\mathrm{}}^2\overline{)}𝒦(k,X)=igC_F\rho _0(k)[N(k_0)+n(k_0)]\overline{)}k\psi (X).$$

Here, $`\mathrm{}m_{\mathrm{}}^2m_{q,\mathrm{}}^2m_{g,\mathrm{}}^2.`$ Substituting into the last equation the expansion

$$\overline{)}𝒦(k,X)=2\pi \delta (k^2)\{\theta (k^0)\overline{)}\mathrm{\Lambda }^+(𝐤,X)+\theta (k^0)\overline{)}\mathrm{\Lambda }^{}(𝐤,X)\},$$
$`(\mathrm{B6})`$
replacing $`𝐤𝐤`$, where it is necessary, and dividing by $`ϵ_k`$, we obtain the improved Eq. $`(2.9^{})`$. By virtue of the fact that $`\mathrm{}m_{\mathrm{}}^20`$, the discrepancy $`\overline{)}\mathrm{\Lambda }^+(𝐤,X)\overline{)}\mathrm{\Lambda }^{}(𝐤,X)`$ will hold. The last case is an actual reflection at the kinetic description level of the discrepancy of the two terms in the integrands of the improved hard thermal loops with external fermions when exchanging $`𝐤𝐤`$, found by Flechsig and Rebhan fle3 .

Now we turn to deriving the improved kinetic equations for $`\delta n_\pm ^\mathrm{\Psi }`$. The equations for functions $`\delta n_\pm ^\mathrm{\Psi }`$ are defined from the equation for the current density $`J_\mu ^\mathrm{\Psi }(k,X)`$ bla2 :

$$[vD_X,J_\mu ^\mathrm{\Psi }(k,X)]=igk_\mu t^a\{\overline{\psi }(X)t^a\overline{)}𝒦(k,X)\overline{)}(k,X)t^a\psi (X)\}$$
$`(\mathrm{B7})`$
by substituting to it the expansion $`J_\mu ^\mathrm{\Psi }(k,X)=k_\mu 4\pi \delta (k^2)\{\theta (k_0)\delta n_+^\mathrm{\Psi }(𝐤,X)+\theta (k_0)\delta n_{}^\mathrm{\Psi }(𝐤,X)\}`$. The relevant modification of the equation for $`\delta n_\pm ^\mathrm{\Psi }`$ here is achieved by the simple requirement of using the expansion (B6) on the right-hand side of Eq. (B7), where the functions $`\overline{)}\mathrm{\Lambda }^\pm `$ obey the improved equations $`(2.8^{})`$.