Prepared for submission to JHEP

## Hot and Dense QCD Shear Viscosity at (almost) NLO

## Isabella Danhoni, Guy D. Moore

Institut f¨r Kernphysik, Technische Universit¨t Darmstadt u a Schlossgartenstraße 2, D-64289 Darmstadt, Germany E-mail:

idanhoni@theorie.ikp.physik.tu-darmstadt.de,guy.moore@physik.tu-darmstadt.de

Abstract: The next-to-leading order weak-coupling shear viscosity of QCD was computed 6 years ago. However, these results have never been applied at finite baryon chemical potential µ , even though intermediate-energy heavy ion collisions and merging neutron stars may explore the Quark-Gluon Plasma in a regime where baryon chemical potentials are large. Here, we extend the next-to-leading order shear viscosity calculations to finite µ , and we show that, while the convergence of the weak-coupling expansion is questionable for achievable plasmas, it is somewhat better at µ &gt; T than at µ = 0.

Keywords: Quark-Gluon Plasma, QCD, Viscosity, Chemical Potential

## Contents

| 1   | Introduction                                             | Introduction                                             |   1 |
|-----|----------------------------------------------------------|----------------------------------------------------------|-----|
| 2   | Kinetic Theory                                           | Kinetic Theory                                           |   2 |
| 3   | Collision operators for an NLO evaluation                | Collision operators for an NLO evaluation                |   5 |
|     | 3.1                                                      | Collision Integrals at LO                                |   6 |
|     | 3.2                                                      | Next-to-Leading Order corrections                        |   9 |
|     |                                                          | 3.2.1 Momentum diffusion                                 |   9 |
|     |                                                          | 3.2.2 Collinear and semi-collinear contributions         |  12 |
|     |                                                          | 3.2.3 Summary                                            |  14 |
| 4   | From zero to nonzero chemical potential                  | From zero to nonzero chemical potential                  |  14 |
| 5   | Results                                                  | Results                                                  |  15 |
| 6   | Summary and Discussion                                   | Summary and Discussion                                   |  18 |
| A   | Self-energies beyond the hard thermal loop approximation | Self-energies beyond the hard thermal loop approximation |  18 |

## 1 Introduction

The deconfined phase of hot and dense QCD is known as the quark-gluon plasma (QGP). Evidence shows that this phase was present in the early universe and that it has been recreated at microscopic scales in heavy ion collision experiments[1]. These results have sparked large interest in this subject over the past few decades, and experiments indicate that the QGP in such collisions is a strongly coupled fluid[2-4] with behavior close to that of a perfect fluid[5-7]. On the other hand, sufficiently early in the history of the early Universe, the running coupling should have been small enough that a perturbative treatment may be valid. For thermodynamical quantities, suitably resummed perturbative calculations work starting at temperatures a few times that of the chiral symmetry restoration transition[8, 9].

̸

In lower-energy collisions, quark-gluon plasma is generated with a significant net baryon number, such that the baryon chemical potential µ B is several times higher than the temperature. This regime cannot be directly studied on the lattice due to the sign problem (or complex action problem). For a chemical potential µ B = 0 the fermion determinant (or Monte Carlo action) becomes complex, which ruins the probabilistic interpretation of this quantity used in lattice QCD calculations[10]. For this reason, the high-density region of the QCD phase diagram is much less well understood than the µ B = 0 axis. As experimental efforts focus on bringing clarity to this region, such as RHIC-BES [11], HADES[12] and CBM at FAIR[13], the need for parallel theoretical development becomes clear.

At high temperatures, we claim that the perturbative series should work better at a high chemical potential, such that µ B &gt; T &gt; T c , than in the case of µ B = 0[14]. This is a consequence of the fact that scattering from quarks is enhanced by a factor of µ 2 B /T 2 . The failure of perturbation theory in hot QCD occurs before the coupling becomes large, due to the strongly-coupled nature of soft gluons. These play a central role in the breakdown of the perturbative expansion for the shear viscosity [15]. But at high chemical potential, where quarks play a larger role, one expects the perturbative series to be better behaved. This means that perturbation theory has a better chance of working at high µ/T than at µ = 0. Although the ( µ, T ) values which can be explored in experiment are unlikely to extend into the region where perturbation theory is directly predictive, the perturbative results can still be used in modeling the ( µ, T ) behavior of the shear viscosity phenomenologically, as was recently done using leading-log results in Ref. [16]. Essentially, perturbation theory serves as a limiting behavior so that η/s can be interpolated based on high-density perturbative, low-density hadron-resonance gas, and intermediate density data-driven estimates.

Previously Arnold, Moore, and Yaffe have computed the shear viscosity at leading perturbative order in the absence of a chemical potential [17], and Ghiglieri et al have extended this calculation to 'almost' next-to-leading order (NLO) [15]. Separately, we have presented the shear viscosity at finite chemical potential at leading-log order [14]. The goal of this paper is to work to almost next-to-leading order at finite chemical potential, both to obtain the most reliable estimates for the viscosity in this regime that we can, and to test how perturbation theory converges in the highµ regime. We start in Section 2 by reviewing the kinetic theory and computational approach that Arnold, Moore, and Yaffe used in [17, 18]. In Section 3, we review how to extend the collision operator to next-toleading order and describe the computational tools required. In Section 4, we show what changes must be made to the calculation to include chemical potentials. The analysis of our calculations and final results are shown in Section 5. We end in Section 6 with a short discussion and summary. But we can present the main conclusions very briefly now. The shear viscosity, rescaled as ηT/ ε ( + P ), is a weak function of µ /T q (with µ q = µ B / 3 the chemical potential of a quark). The convergence of the perturbative series improves when N µ f 2 q &gt; N π T c 2 2 , but the improvement is somewhat modest. For more complete details, see Section 5.

## 2 Kinetic Theory

In this section, we briefly summarize how shear viscosity is defined and how it can be computed to leading order (LO) in the Effective Kinetic Theory (EKT). Shear viscosity is a property relevant for systems with space-nonuniform flow velocity, meaning that the energy and momentum densities vary through space, but in a manner which is sufficiently smooth that the system is locally near equilibrium. In this case we can exchange the conserved densities T µ 0 for the local 4-velocity and energy density, defined with the Landau-Lifshitz convention [19]: u µ ( x, t ) is the timelike unit eigen(four)vector of T µν and ε is its eigenvalue,

T µν u ν = εu µ with 1 u u µ µ = -1. For uniform u , ε µ the system will be in equilibrium with T µν = εu u µ ν + Pg µν , where P the pressure is determined from ε by the equilibrium equation of state.

However, for nonuniform u µ the system will depart from equilibrium and T µν will not follow its ideal form. Defining the projection operator into the local rest frame,

$$\Delta ^ { \mu \nu } \equiv g ^ { \mu \nu } + u ^ { \mu } u ^ { \mu } \,,$$

the shear tensor is defined as

$$\sigma ^ { \mu \nu } \equiv \nabla ^ { \mu } u ^ { \nu } + \nabla ^ { \nu } u ^ { \mu } - \frac { 2 } { 3 } \Delta ^ { \mu \nu } \nabla _ { \alpha } u ^ { \alpha } \,.$$

For small and slowly varying σ ij , linear response theory predicts that the stress tensor should change by an amount linear in the shear tensor: in the local rest frame,

$$T _ { i j } = P g _ { i j } - \eta \, \sigma _ { i j }$$

with P the pressure and η , the shear viscosity, defined as the coefficient of the linear response of T ij to nonvanishing σ ij . If the fluid flow also has a divergence θ = ∂ u i i , then an additional term T ij = . . . -ζθ ∆ ij would be present, with ζ the bulk viscosity. Bulk viscosity requires a different analysis [20] which we will not attempt here.

Many (not all) theories have a kinetic description, and this should include nonabelian gauge theories in the weak coupling regime. In this case one can give a coarse-grained description of the state of the system in terms of the occupancies of each field type at each point in time and phase space, f a ( ⃗, ⃗, t p x ). In equilibrium and at lowest order in the coupling, we have

$$f _ { 0 } ^ { a } ( \vec { p }, \vec { x }, t ) = \frac { 1 } { \exp ( \gamma \beta ( p ^ { 0 } - u _ { i } p ^ { i } - \mu _ { a } ) \mp 1 } = \frac { 1 } { \exp ( \beta ( - u _ { \mu } P ^ { \mu } - \mu _ { a } ) ) \mp 1 }.$$

Here u µ = ( γ, γ⃗ v ) is the local velocity 4-vector as before, γ is the associated relativistic gamma-factor, and β = β ⃗ x, t ( ) = 1 /T is the local temperature. The sign ∓ 1 is -for bosons and + for fermions (the upper sign will always apply to bosons), and µ a is the chemical potential for species a , which will be nonzero for quark species and opposite for a particle and its antiparticle. In this work, we will consider SU (3) gauge theory with N f vectorlike quark species each with the same chemical potential.

Shear viscosity, as well as other transport coefficients, is theoretically defined in a Kubo-type formulae relating them to the zero-frequency transport limit of the spectral functions of two-point correlators of the appropriate operators. Performing a perturbative evaluation at leading order using the diagrammatic approach based directly on these Kubo formulas would require complex resummations of numerous sub-diagrams to all orders. This has been done successfully for scalar field theory [21] and for QED for electrical conductivity [22] and shear viscosity [23], and in each case the result proved to be the same

1 We use a mostly-positive metric tensor g µν = Diag[ -1 +1 +1 +1], and natural units, , , , c = 1, ℏ = 1. Capital letters P, K are 4-vectors, lower case letters ⃗, ⃗ p k or p , k i i are the space components, p , k 0 0 are the time components, and p, k are the magnitudes of the space components.

as a kinetic description. In the case of full QCD such a diagrammatic calculation has never been successfully attempted, but physical arguments indicate that the result should again justify a kinetic description [24]. Therefore, we adopt the linearized version of the EKT developed in [24].

One starts as usual with a Boltzmann equation

$$\left [ \vec { v } _ { p } \cdot \frac { \partial } { \partial \vec { x } } \right ] f ^ { a } ( \vec { p }, \vec { x }, t ) = - \mathcal { C } ^ { a } [ f ],$$

where f a refers to the nonequilibrium form of the statistical function, which we linearize about the equilibrium form

$$f ^ { a } ( \vec { p }, \vec { x } ) = f _ { 0 } ^ { a } ( \vec { p }, \vec { x } ) + f _ { 0 } ( 1 \pm f _ { 0 } ) \, f _ { 1 } ^ { a } ( \vec { p }, \vec { x } ) \,,$$

according to the Landau-Lifshitz conventions under which f 0 contributes all of the energy, momentum density, and net particle numbers in the system. The choice to normalize f 1 with a factor of f 0 (1 ± f 0 ) will be convenient in what follows. For a more detailed description one can refer to [25] which developed the leading-log application of Boltzmann equations to the evaluation of the QCD shear viscosity, [17], which extended the treatment to leading order and established most of our notation and techniques, and [14], which extended the treatment to finite chemical potentials, albeit at leading-logarithmic order. We will save the discussion of the collision operator C a [ f ] for the next section.

The departure from equilibrium must first arise at linear order in σ ij . Using rotational invariance, one can then write the departure from equilibrium as

$$f _ { 1 } ^ { a } ( \vec { p }, \vec { x } ) = \beta ^ { 2 } X _ { i j } ( \vec { x } ) \chi _ { i j } ^ { a } ( \vec { p } ) = \beta ^ { 2 } X _ { i j } ( \vec { x } ) I _ { i j } ( \vec { p } ) \chi ^ { a } ( p ) \,.$$

Here χ a ( p ) is a species-dependent, pure scalar function of p which characterizes the departure from equilibrium in the presence of shear stress, and

$$X _ { i j } \equiv \frac { 1 } { \sqrt { 6 } } \sigma _ { i j }$$

$$I _ { i j } ( \vec { p } ) \equiv \sqrt { \frac { 3 } { 2 } } \left ( \hat { p } _ { i } \hat { p } _ { j } - \frac { 1 } { 3 } \delta _ { i j } \right ) \,.$$

The normalization is explained in Ref. [17] and is chosen so that I ij is a properly normalized projection operator. Evaluating the space derivatives on the lefthand side of Eq. (2.5) one finds to lowest order,

$$\vec { v } _ { p } \cdot \frac { \partial } { \partial \vec { x } } f _ { 0 } ^ { a } ( \vec { p }, \vec { x }, t ) = \beta ^ { 2 } X _ { i j } S _ { i j } ^ { a }$$

with

$$S _ { i j } ^ { a } ( \vec { p } ) \equiv p ^ { 0 } \, T \, f _ { 0 } ^ { a } ( 1 \pm f _ { 0 } ^ { a } ) I _ { i j } ( \vec { p } ) \,.$$

Using this and inserting Eq. (2.7) into the RHS, Eq. (2.5) reduces to

$$\beta ^ { 2 } X _ { i j } \ S _ { i j } ^ { a } ( \vec { p } ) = \mathcal { C } ^ { a } [ \beta ^ { 2 } X _ { i j } \chi _ { i j } ^ { a } ( \vec { p } ) ] \ \text{ and therefore} \ \ S _ { i j } ^ { a } ( \vec { p } ) = \mathcal { C } \chi _ { i j } ^ { a } ( \vec { p } ), \quad ( 2. 1 2 )$$

where C χ a ij ( ⃗ p ) means the result of factoring β X 2 ij out of the collision operator.

This equation is solved using a variational method. This strategy was used in [17, 18] in leading log and leading order calculations. First one defines an inner product on departures from equilibrium:

$$( f, g ) = \beta ^ { 3 } \sum _ { a } ^ { \tilde { \tt f h c } } \int _ { \vec { p } } f ^ { a } ( \vec { p } ) g ^ { a } ( \vec { p } )$$

It can be shown that the collision operator is Hermitian with respect to this inner product. We define

$$\mathcal { Q } [ \chi ] = \left ( \chi _ { i j }, S _ { i j } \right ) - \frac { 1 } { 2 } \left ( \chi _ { i j }, \mathcal { C } \chi _ { i j } \right ) \,,$$

which can be viewed as a functional over the space of χ a ( p ) values. That is, the functional Q takes its maximum value when the Boltzmann equation is satisfied. Furthermore, its value at this maximum,

$$Q _ { \max } = \frac { 1 } { 2 } ( \chi _ { i j }, S _ { i j } ) = \frac { 1 } { 2 } ( \chi _ { i j }, \mathcal { C } \chi _ { i j } ) = \frac { 1 } { 2 } ( S _ { i j }, \mathcal { C } ^ { - 1 } S _ { i j } ) \,,$$

determines the viscosity:

$$\eta = \frac { 2 } { 1 5 } Q _ { \max }.$$

As known from other variational problems, one can obtain highly accurate approximate results by performing a restricted extremization within a well-chosen finite dimensional subspace. For this purpose, we expand the χ p ( ) functions into a finite basis set:

$$\chi ( p ) = \sum _ { m = 1 } ^ { N } a _ { m } \phi ^ { ( m ) } ( p ).$$

We will use the variational Ansatz proposed 2 in Ref. [18, 26]:

$$\phi ^ { ( m ) } = \frac { p ( p / T ) ^ { m } } { ( 1 + p / T ) ^ { N - 1 } }, \ \ m = 1, \dots, N \,.$$

The subspace spanned by this basis strictly increases as N is increased.

## 3 Collision operators for an NLO evaluation

At vanishing chemical potential, the perturbative framework suggests that the QGP predominantly comprises quasiparticle excitations with momenta of the order of the temperature. Their dispersion and occupancy are corrected only at O ( g 2 ), which for an O ( g )-subleading calculation means that we can ignore corrections both to the lowest-order entropy density s and to the lowest-order particle dispersion relations (except for splitting, see below). Therefore NLO corrections will only be needed in the collision integral, which we turn to next. The LO calculation [17] requires two types of scattering processes, the 2 ↔ 2 scatterings with all hard ( O ( T )) external particles and 1 ↔ 2 effective splitting processes between hard participants.

2 To our knowledge this variational Ansatz was first used in Ref. [26], but it is not actually presented there; we thank Larry Yaffe for private communications on this point.

At NLO, all the scatterings and corrections to be added are the ones suppressed by a single power of g . As shown in detail in [15], there are only a few such O ( g ) subleading effects, which can be described as follows. The rate of soft 2 ↔ 2 scattering is modified to be considered as an additional momentum-diffusion coefficient δq ˆ. This along with an O ( g ) correction to the medium corrections to dispersion provides an O ( g ) shift in the 1 ↔ 2 splitting rate. Next, the 1 ↔ 2 splitting rate must be corrected wherever one participant of the process is 'soft', p ∼ gT , or the opening angle is less collinear. Finally, the numerical implementation of the LO scattering kernel [17] already resums a small amount of these NLO effects, requiring a subtraction to δq ˆ and the longitudinal momentum diffusion, δq ˆ . L This section will briefly review the collision operator at LO and NLO at vanishing chemical potential.

## 3.1 Collision Integrals at LO

The leading-order collision operator encodes the contribution of tree-level 2 ↔ 2 scattering processes, with Hard-Loop resummed propagators in the soft-sensitive channels, as well as collinear, effective 1 ↔ 2 processes resumming the effect of splitting caused by any number of soft scatterings. Both processes contribute to order g T 4 to the collision operator; a subset of C 2 ↔ 2 a [ f ] is logarithmically enhanced, g T 4 ln(1 /g ), due to the sensitivity to the soft scale gT . C 2 ↔ 2 a [ f ] and C 1 ↔ 2 a [ f ] are described in detail in [17, 24]. Therefore, the collision term has the form:

$$\mathcal { C } _ { a } [ f ] = \mathcal { C } _ { a } ^ { 2 \leftrightarrow 2 } [ f ] + \mathcal { C } _ { a } ^ { 1 \leftrightarrow 2 } [ f ]$$

The matrix elements that contribute to 2 ↔ 2 processes at leading order come from the diagrams shown in Fig. 1. In terms of the collision operator, these contributions can be written as:

$$\infty \text{ window } \infty. \\ \mathcal { C } ^ { a } [ f ] ( \vec { p } ) = \frac { 1 } { 2 } \sum _ { b c d } \int _ { \vec { k }, \vec { p }, \vec { k ^ { \prime } } } & \frac { | \mathcal { M } _ { a b c d } ( P, K, P ^ { \prime }, K ^ { \prime } ) | ^ { 2 } } { 2 p ^ { 0 } \, 2 k ^ { 0 } \, 2 p ^ { \prime 0 } \, 2 k ^ { \prime 0 } } ( 2 \pi ) ^ { 4 } \delta ^ { 4 } ( P + K - P ^ { \prime } - K ^ { \prime } ) \\ \times & \left \{ f ^ { a } ( \vec { p } ) f ^ { b } ( \vec { k } ) [ 1 \pm f ^ { c } ( \vec { p ^ { \prime } } ) ] [ 1 \pm f ^ { d } ( \vec { k ^ { \prime } } ) ] \\ & \quad - f ^ { c } ( \vec { p ^ { \prime } } ) f ^ { d } ( \vec { k ^ { \prime } } ) [ 1 \pm f ^ { a } ( \vec { p } ) ] [ 1 \pm f ^ { b } ( \vec { k } ) ] \right \}$$

where the incoming/outgoing momenta ⃗, ⃗ p k and ⃗ , ⃗ p ′ k ′ are all on shell, p 0 = p . The squared matrix element ∣ ∣ M ab cd ∣ ∣ is hard thermal loop (HTL) resummed [27, 28] and is summed over all spin polarizations and colors. We write ∫ ⃗ k = ∫ d ⃗ 3 k (2 π ) 3 to simplify the notation.

Diagrams (D) and (E) are processes involving t (or u) channel fermion exchange. These processes have a logarithmic infrared divergence when computed using free propagators. To cut off this divergence, one can include the retarded HTL fermionic thermal self-energy, as was originally done in Ref. [17]. It was later shown in Ref. [29] that there is a simpler treatment which is also correct at the leading order. The most important feature of the HTL self-energy insertion is that it correctly reproduces the total scattering rate for this process; since the process is nearly isotropic, the detailed angle dependence is not relevant. A mass-like insertion, with the correctly chosen coefficient, reproduces this rate and is

Figure 1 . The diagrams relevant to leading-order scattering in QCD, showing our labeling conventions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00524v1-figure-a61dc969-f129-4f2b-afb2-7e27da9e1c4d.jpg)

much simpler to use in calculation. Specifically for a soft fermion exchange in the t or u channel, one can simply use:

$$\frac { u - s } { t } \rightarrow \frac { u - s } { t } \frac { q ^ { 2 } } { q ^ { 2 } + \xi _ { q } ^ { 2 } m _ { \infty, q } ^ { 2 } }$$

where q is the exchanged momentum, ξ q = e/ 2 ≃ 1 359, and . m 2 ∞ ,q is the medium-induced effective quark mass-squared for quarks, given by

$$\cdots \frac { \mu \dots \dots \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

We have allowed ourselves to make this simplification - also because the cumulative effect of these diagrams is rather small and is exponentially suppressed ∝ e -µ/T in the large µ/T region.

On the other hand, diagrams (A),(B), and (C) represent processes involving t or uchannel gauge boson exchange that require the inclusion of the thermal gauge boson selfenergy on the internal propagator to cut off the infrared sensitivity of these processes. The matrix element of interest is ( s 2 + u 2 ) /t 2 which can be written as:

$$\frac { 1 } { 8 } D _ { \mu \alpha } ^ { \text{ret} } D _ { \nu \beta } ^ { \text{adv} } \text{Tr} \left ( \psi \gamma ^ { \mu } k \gamma ^ { \nu } \right ) \text{Tr} \left ( \psi ^ { \prime } \gamma ^ { \alpha } k ^ { \prime \prime } \gamma ^ { \beta } \right )$$

where D ret µα and D adv µα , are the retarded and advanced boson propagators, respectively. Following the same notation as [30], the plasma frame frequency and momentum carried by the gauge boson propagator are denoted as ω and q . The retarded gauge boson propagator can be conveniently expressed in Coulomb gauge as:

$$\ u y \exp r e s s e { \emph { m o u m o g a g u e s } } ; \\ D _ { 0 0 } ^ { \text{ret} } ( \omega, \mathbf q ) & = \frac { - 1 } { q ^ { 2 } - \Pi _ { 0 0 } ^ { \text{ret} } ( \omega, \mathbf q ) } \\ D _ { i j } ^ { \text{ret} } ( \omega, \mathbf q ) & = \frac { \delta _ { i j } - \hat { \mathbf q } _ { i } \hat { \mathbf q } _ { j } } { q ^ { 2 } - \omega ^ { 2 } + \Pi _ { T } ^ { \text{ret} } ( \omega, \mathbf p ) } \\ D _ { 0 i } ^ { \text{ret} } ( \omega, \mathbf q ) & = D _ { i 0 } ^ { \text{ret} } ( \omega, \mathbf q ) = 0.$$

The equilibrium transverse and longitudinal gauge boson self-energies in the HTL approximation are [31, 32]

$$\Pi _ { \text{T} } ^ { \text{ret} } ( \omega, q ) = m _ { D } ^ { 2 } \left \{ \frac { \omega ^ { 2 } } { 2 q ^ { 2 } } + \frac { \omega \left ( q ^ { 2 } - \omega ^ { 2 } \right ) } { 4 q ^ { 3 } } \left [ \ln \left ( \frac { q + \omega } { q - \omega } \right ) - i \, \pi \right ] \right \} \,,$$

$$\Pi ^ { \text{ret} } ( \omega, q ) = m _ { D } ^ { 2 } \left \{ 1 - \frac { \omega } { 2 q } \left [ \ln \left ( \frac { q + \omega } { q - \omega } \right ) - i \, \pi \right ] \right \},$$

where we have assumed | ω &lt; q | , which is the only case of relevance. Here m 2 D is the squared Debye mass, which we will give for the case of nonzero µ, T in Eq. (4.1). The advanced propagators are the complex conjugates of the retarded ones.

The remaining diagrams and the interference terms are all finite and can be computed directly. One should refer to [17] for more details on these calculations.

The second type of contribution to leading order comes from 1 ↔ 2 nearly collinear processes. These can be approximated as strictly collinear, which only produces relative changes suppressed by at least one power of g. Corrections at this order will be treated in the next subsection. The leading-order computation proceeds as follows [24]:

$$\text{the next subsection. } I \text{$ne$ leading-order computation proceeds as follows [24}$} \colon \\ \mathcal { C } _ { a } ^ { 1 \mapsto 2 } [ f ] ( \vec { p } ) = & \frac { ( 2 \pi ) ^ { 3 } } { 2 | \vec { p } | ^ { 2 } \nu _ { a } } \sum _ { b c } \int _ { 0 } ^ { \infty } d p ^ { \prime } d k \delta ( | \vec { p } | - p ^ { \prime } - k ^ { \prime } ) \gamma _ { b c } ^ { a } ( \vec { p } ; p ^ { \prime } \hat { p }, k ^ { \prime } \hat { p } ) \\ & \left \{ f ^ { a } ( \vec { p } ) [ 1 \pm f ^ { b } ( p ^ { \prime } \hat { p } ) ] [ 1 \pm f ^ { c } ( k \hat { p } ) ] - f ^ { b } ( p ^ { \prime } \hat { p } ) f ^ { c } ( k \hat { p } ) [ 1 \pm f ^ { a } ( \vec { p } ) ] \right \} \\ + & \frac { ( 2 \pi ) ^ { 3 } } { 2 | \vec { p } | ^ { 2 } \nu _ { a } } \sum _ { b c } \int _ { 0 } ^ { \infty } d k d p ^ { \prime } \delta ( | \vec { p } | + k - p ^ { \prime } ) \gamma _ { a b } ^ { c } ( p ^ { \prime } \hat { p } ; p, k \hat { p } ) \\ & \left \{ f ^ { a } ( \vec { p } ) f ^ { b } ( k \hat { p } ) [ 1 \pm f ^ { c } ( p ^ { \prime } \hat { p } ) ] - f ^ { c } ( p ^ { \prime } \hat { p } ) [ 1 \pm f ^ { a } ( \hat { p } ) ] [ 1 \pm f ^ { b } ( k \hat { p } ) ] \right \}, \quad ( 3. 9 ) \\ \cdot \quad a \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \end{cases}$$

where γ a bc represents the splitting rate, and is given by [24]:

$$\gamma _ { b c } ^ { a } ( p ; p - k, k ) = \frac { g ^ { 2 } d _ { R _ { b } } C _ { R _ { b } } } { 6 4 \pi ^ { 4 } } \left \{ \frac { \frac { p ^ { 4 } + k ^ { 4 } + ( p - k ) ^ { 4 } } { p ^ { 3 } _ { k ^ { 3 } } ( p - k ) ^ { 2 } } } { p ^ { 2 } _ { k ^ { 3 } } ( p - k ) ^ { 2 } } } \, \begin{array} { c } g \leftrightarrow g g \\ \frac { p ^ { 2 } + ( p - k ) ^ { 2 } } { p ^ { 2 } _ { k ^ { 3 } } ( p - k ) ^ { 2 } } } & q \leftrightarrow q g \\ \frac { p ^ { 2 } + ( p - k ) ^ { 2 } } { p ^ { 3 } _ { k ^ { 2 } } ( p - k ) ^ { 2 } } & g \leftrightarrow q \bar { q } \end{array} \right \} \times \int \frac { d \mathbf h ^ { 2 } } { ( 2 \pi ) ^ { 2 } } 2 \mathbf h \cdot \text{Re} \mathbf B,$$

where h = p × k is a transverse, two-dimensional vector related to the transverse momentum picked up during the splitting process. And F b can be obtained through the following integral equation:

$$2 = & i \delta E ( \mathbf h ) \mathbf F _ { b } ( \mathbf h ) + \int \frac { d q _ { \perp } ^ { 2 } } { ( 2 \pi ) ^ { 2 } } \bar { C } ( q _ { \perp } ) \{ ( C _ { R _ { b } } - C _ { A } / 2 ) \left [ \mathbf F _ { b } ( \mathbf h ) - \mathbf F _ { b } ( \mathbf F - k \mathbf q _ { \perp } ) \right ] } \\ & C _ { A } / 2 [ \mathbf F _ { b } ( \mathbf h ) - \mathbf F _ { b } ( \mathbf h + p \mathbf q _ { \perp } ) ] + C _ { A } / 2 [ \mathbf F _ { b } ( \mathbf h ) - \mathbf F _ { b } ( \mathbf h - ( p - k ) \mathbf q _ { \perp } ) ] \}.$$

For the case of g ↔ qq ¯, C R b -C / A 2 = C F -C / A 2 multiplies the term with F b ( h -p q ⊥ ) rather than F b ( h -k q ⊥ ). Here ¯ ( C q ⊥ ) represents the leading-order transverse scattering kernel:

$$\bar { C } ( q _ { \perp } ) = \frac { g ^ { 2 } T m _ { D } ^ { 2 } } { \mathbf q _ { \perp } ^ { 2 } ( \mathbf q _ { \perp } ^ { 2 } + m _ { D } ^ { 2 } ) },$$

and δE ( h ) is the energy difference between initial and final collinear particles, expanded to first nontrivial order in thermal masses and transverse momenta:

$$\delta E ( \mathbf h ) = \frac { h ^ { 2 } } { 2 p k ( p - k ) } + \frac { m _ { \infty, c } ^ { 2 } } { 2 k } + \frac { m _ { \infty, b } ^ { 2 } } { 2 ( p - k ) } - \frac { m _ { \infty, a } ^ { 2 } } { 2 p }.$$

We encountered m 2 ∞ ,q already in Eq. (3.4). For gluons, m 2 ∞ ,g = m / 2 D 2.

## 3.2 Next-to-Leading Order corrections

We now focus on the next-to-leading order corrections. Treating kinetic theory to NLO was pioneered in Ref. [33, 34] and applied to QCD shear viscosity in Ref. [15]. In this work, the authors observed that the analytically complex parts of the collision operator do not receive O ( g ) corrections, which instead only modify collinear exchange, incorporated by changing γ a bc , and soft exchange, which can be modeled as momentum-diffusion and identity-changing processes. We briefly review the extended treatment that captures these O g ( ) corrections, and then we summarize it and present an explicit form of the collision operator at NLO.

## 3.2.1 Momentum diffusion

The small momentum exchange part of diagrams (A), (B), (C) can be well approximated by momentum diffusion for the particles on each line. Beyond leading order, there are additional small-momentum-exchange effects which are only suppressed by O ( g ) relative to LO, which must be incorporated. One such effect is the possibility for the gluon in diagrams (A), (B) to carry a soft momentum. Other effects involve the interference between multiple soft scatterings, soft gluon emission during scattering, and nonabelian interactions between scattering processes, such as those depicted in Figure 2.

Figure 2 . Some scattering effects which affect momentum diffusion at NLO.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00524v1-figure-8623b670-16ff-4a0e-b684-3c0fc25b6ffa.jpg)

Because the momentum exchanges involved in these processes are small (that is the reason they can occur at O ( g ) and not O ( g 2 )), these terms can be adequately described in terms of Langevin dynamics, that is, in terms of momentum drag and longitudinal and transverse momentum diffusion. This is analyzed in the context of heavy quark dyamics in Ref. [35] and jet energy loss in Ref. [34], and applied to shear viscosity in Ref. [15]. As shown there, the effective collision term due to soft scatterings, cut off at a scale µ ⊥ , is

$$( C _ { \text{diff} } ^ { 2 \leftrightarrow 2 } f _ { 1 } ) ^ { a } ( \boldsymbol p ) = - \frac { 1 } { 2 } \frac { \partial } { \partial p ^ { j } } \, \hat { q } _ { a } ^ { i j } \, f _ { 0 } ^ { a } ( p ) ( 1 \pm f _ { 0 } ^ { a } ( p ) ) \frac { \partial f _ { 1 } ^ { a } ( \boldsymbol p ) } { \partial p ^ { j } } + \text{gain-terms} \,,$$

where

$$\hat { q } _ { a } ^ { i j } = \hat { q } _ { a, L } \hat { p } ^ { i } \hat { p } ^ { j } + \frac { 1 } { 2 } \hat { q } _ { a } \left ( \delta ^ { i j } - \hat { p } ^ { i } \hat { p } ^ { j } \right )$$

represents the momentum diffusion parallel and perpendicular to the particle's momentum through ˆ q L and ˆ respectively. q

Contracting with the departure from equilibrium and integrating over p , the contribution to the inner product of Eq. (2.14) takes the form

$$( f _ { 1 }, C _ { \text{diff} } ^ { 2 \rightarrow 2 } f _ { 1 } ) \equiv ( f _ { 1 }, C _ { \text{diff} } ^ { 2 \rightarrow 2 } f _ { 1 } ) \Big | _ { \text{loss} } + ( f _ { 1 }, C _ { \text{diff} } ^ { 2 \rightarrow 2 } f _ { 1 } ) \Big | _ { \text{gain} } \,,$$

where

$$( f _ { 1 }, C _ { \text{diff} } ^ { 2 \to 2 } f _ { 1 } ) | _ { \text{loss} } = \frac { 1 } { 2 } \beta ^ { 3 } \sum _ { a } \nu _ { a } \, \int _ { p } f _ { 0 } ^ { a } ( p ) ( 1 \pm f _ { 0 } ^ { a } ( p ) ) \, \dot { q } _ { a } ^ { i j } \, \frac { \partial f _ { 1 } ^ { a } ( p ) } { \partial p ^ { i } } \frac { \partial f _ { 1 } ^ { a } ( p ) } { \partial p ^ { j } } \,,$$

The gain collision term takes the form

$$\text{$megan common term cases use inm} \\ \text{gain-terms} & = \frac { 1 } { 2 } \frac { \partial } { \partial p ^ { i } } \left ( f _ { 0 } ^ { a } ( p ) ( 1 \pm f _ { 0 } ^ { a } ( p ) ) \, \sum _ { b } \nu _ { b } \int _ { k } \mathcal { C } _ { a b } ^ { i j } ( \hat { p } \cdot \hat { k } ) \, f _ { 0 } ^ { b } ( k ) ( 1 \pm f _ { 0 } ^ { b } ( k ) ) \frac { \partial f _ { 1 } ^ { b } ( k ) } { \partial k ^ { j } } \right ) \\ ( f _ { 1 }, C _ { \text{diff} } ^ { 2 \to 2 } f _ { 1 } ) | _ { \text{gain} } & = \, - \frac { \beta ^ { 3 } } { 2 } \sum _ { a b } \nu _ { a } \nu _ { b } \int _ { p k } f _ { 0 } ^ { a } ( p ) ( 1 \pm f _ { 0 } ^ { a } ( p ) ) f _ { 0 } ^ { b } ( k ) ( 1 \pm f _ { 0 } ^ { b } ( k ) ) \\ & \times \mathcal { C } _ { a b } ^ { i j } ( \hat { p } \cdot \hat { k } ) \frac { \partial f _ { 1 } ^ { a } ( p ) } { \partial p ^ { i } } \frac { \partial f _ { 1 } ^ { b } ( k ) } { \partial k ^ { j } }, \\ \text{where the anular function is}$$

where the angular function is

$$\mathcal { C } ^ { i j } _ { a b } ( \hat { p } \cdot \hat { k } ) = \frac { g ^ { 4 } C _ { R _ { a } } C _ { R _ { b } } } { d _ { A } } \int \frac { d ^ { 4 } Q } { ( 2 \pi ) ^ { 4 } } | G ^ { R } _ { \mu \nu } ( Q ) v ^ { \mu } _ { \mathbf p } v ^ { \nu } _ { \mathbf k } | ^ { 2 } \, 2 \pi \delta ( v _ { \mathbf p } \cdot Q ) 2 \pi \delta ( v _ { \mathbf k } \cdot Q ) q ^ { i } q ^ { j } \,,$$

which is related to ˆ q ij a through

$$\hat { q } _ { a } ^ { i j } = \sum _ { b } \nu _ { b } \int _ { k } f _ { 0 } ^ { b } ( k ) ( 1 \pm f _ { 0 } ^ { b } ( k ) ) \, \mathcal { C } _ { a b } ^ { i j } ( \hat { \mathbf p } \cdot \hat { k } ) \,.$$

Beyond leading order it turns out to be much easier to determine ˆ q ij a than C ab ij ( ˆ p · ˆ ), k since light-cone techniques work for the former but not the latter. Therefore, in practice we will not be able to compute this gain term. Since our best estimates based on leading-order analyses indicate that angular cancellations make this term small, we will instead drop it. This is why we generally refer to our results as 'almost next-to-leading order.'

Examining the expression for ( f , C 1 2 ↔ 2 diff f 1 ) ∣ ∣ loss , it becomes evident that it involves a single light-like vector denoted as v p . Using the causality and KMS properties of such lightlike correlators [36], these leading-order soft contributions to ˆ and ˆ q q L can be evaluated in closed form [33, 34, 37],

$$\hat { q } ^ { a } \Big | _ { \text{soft} } & = \frac { g ^ { 2 } C _ { R _ { a } } T m _ { D } ^ { 2 } } { 2 \pi } \ln \frac { \mu _ { \perp } } { m _ { D } }, \quad \hat { q } ^ { a } _ { L } \Big | _ { \text{soft} } = \frac { g ^ { 2 } C _ { R _ { a } } T m _ { D } ^ { 2 } } { 4 \pi } \ln \frac { \sqrt { 2 } \mu _ { \perp } } { m _ { D } }. \quad \quad ( 3. 2 1 )$$

In these expressions, gT ≪ µ ⊥ ≪ T is a cutoff on the q ⊥ ≡ √ q 2 -ω 2 integration, which separates the soft from the hard scale. In the region where ω, q ≳ T , and the bare matrix elements can be used to compute the hard contribution to the collision operator C 2 ↔ 2 , there is a cancellation on this cutoff dependence. Since light-like separated points are effectively causally disconnected as far as the soft gauge fields are concerned, ˆ and ˆ q q L have a simple form. Using the explicit form of the angular dependence of f a 1 ( p ) = β χ 2 ij ( p ) X ij , straightforward analysis shows that the loss term reduces to

$$\cdots \circ \cdots \circ \cdots \circ \cdots \circ \cdots \circ \cdots \\ \left ( \chi _ { i j }, C _ { \text{lossshort} } ^ { 2 \leftrightarrow 2 } \chi _ { i j } \right ) \right | _ { \text{loss} } = \frac { \beta ^ { 3 } } { 2 } \sum _ { a } \nu _ { a } \int _ { p } f _ { 0 } ^ { a } ( p ) ( 1 \pm f _ { 0 } ^ { a } ( p ) ) \times \\ \left [ ( \chi ^ { a } ( p ) ^ { \prime } ) ^ { 2 } q _ { L } ^ { a } \right | _ { \text{soft} } + \frac { \ell ( \ell + 1 ) \chi ^ { a } ( p ) ^ { 2 } } { 2 p ^ { 2 } } q ^ { a } \right | _ { \text{soft} } \right ], \ \ ( 3. 2 2 ) \\ \text{which is the most useful form for evaluation} \text{the transform coefficients numerically}$$

which is the most useful form for evaluating the transport coefficients numerically.

The next-to-leading order value of ˆ q a has been computed [36]. This result can be simply included by adding it to the above results [15]:

$$\delta \hat { q } ^ { a } = \frac { g ^ { 4 } C _ { R _ { a } } C _ { A } m _ { D } T ^ { 2 } } { 3 2 \pi ^ { 2 } } \left ( 3 \pi ^ { 2 } + 1 0 - 4 \ln 2 \right ) \,.$$

the shift included in the soft gluon exchange diffusion rate gives

$$\left ( \chi _ { i j }, C ^ { \delta \bar { q } } \chi _ { i j } \right ) = \frac { \beta ^ { 3 } } { 2 } \sum _ { a } \nu _ { a } \, \delta \dot { q } ^ { a } \int _ { p } f _ { 0 } ^ { a } ( p ) [ 1 \pm f _ { 0 } ^ { a } ( p ) ] \frac { \ell ( \ell + 1 ) \chi ^ { a } ( p ) ^ { 2 } } { 2 p ^ { 2 } }. \quad \ \ ( 3. 2 4 )$$

As discussed in Ref. [15], the precise procedure described in the last subsection for including HTL effects in the leading-order calculation also introduces some NLO effects in that calculation. In particular, the phase space and statistical functions for scattering from soft gluons incorporate some of the O ( g ) effects which we have just described. Therefore, the leading-order procedure we use will not in fact reproduce Eq. (3.21), but will reproduce that expression plus some O ( g ) additional terms. These contributions must then be subtracted from the O ( g ) contribution we include here, in order to avoid overcounting such effects. Fortunately, since the effects arise from scattering off of gluons, they are unmodified by the inclusion of quark chemical potentials, and we can take over the analysis of Ref. [15]. The resulting counterterm can be written as

$$\int \colon [ \infty ] {. } \, \text{$1\infty$ rounding counter} \, \text{and} \, \infty \, \text{when} \, \infty \\ \left ( \chi _ { i j }, C _ { \mathcal { O } ( g ) } ^ { 2 \leftrightarrow 2 } \, & \text{finite} \, g \, \chi _ { i j } \right ) = \frac { d _ { A } C _ { A } g ^ { 4 } m _ { D } } { 3 2 \pi ^ { 5 } T } \sum _ { a } T _ { R _ { a } } \int _ { 0 } ^ { \infty } d p \, p ^ { 2 } \, f _ { 0 } ^ { a } ( p ) [ 1 \pm f _ { 0 } ^ { a } ( p ) ] \\ \times & \left \{ 4. 2 6 9 5 [ ( \chi ^ { a } ( p ) ^ { \prime } ) ^ { 2 } + ( \chi ^ { g } ( 0 ) ^ { \prime } ) ^ { 2 } ] + 7. 1 7 6 9 \frac { \ell ( \ell + 1 ) } { 2 p ^ { 2 } } \left [ ( \chi ^ { a } ( p ) ) ^ { 2 } + ( p \chi ^ { g } ( 0 ) ^ { \prime } ) ^ { 2 } \right ]$$

$$+ 1 8. 0 6 6 9 \, \delta _ { \ell 2 } [ \chi ^ { g } ( 0 ) ^ { \prime } ] ^ { 2 } \Big \}, \quad ( 3. 2 5 )$$

and the soft fermion exchange counterterm is given by

$$\left ( \chi _ { i j }, C _ { \mathcal { O } ( g ) } ^ { 2 \leftrightarrow 2 } \chi _ { i j } \right ) = \stackrel { \rangle } { \frac { d _ { F } C _ { F } ^ { 2 } N _ { \mathfrak { t } } g ^ { 4 } m _ { \infty } } { 3 2 \pi ^ { 5 } T ^ { 2 } } } 9. 9 5 2 6 8 \int _ { 0 } ^ { \infty } d k \, k \, f _ { 0 } ^ { q } ( k ) [ 1 + f _ { 0 } ^ { g } ( k ) ] \\ & \times \left [ ( \chi ^ { q } ( k ) - \chi ^ { g } ( k ) ) ^ { 2 } + ( \chi ^ { q } ( 0 ) ) ^ { 2 } \right ].$$

One should refer to [15] for a detailed derivation of these counterterms.

## 3.2.2 Collinear and semi-collinear contributions

The collinear splitting process, Eq. (3.9) to Eq. (3.13), involves the same momentum diffusion effects we just described. And as we just described, the strength of momentum diffusion receives O ( g ) corrections. Furthermore, Eq. (3.13) shows that collinear splitting is sensitive to thermal effective masses. These also turn out to receive O ( g ) corrections. Therefore we need to analyze how these two corrections impact the collinear splitting rate. To conduct the analysis, we will follow Ref. [37] and Fourier transform Eq. (3.9) from transverse-momentum q ⊥ to transverse impact parameter b ⊥ space. This not only diagonalizes the collision kernel C k ( ⊥ ), but also turns an integral equation into a differential equation. Furthermore, the source on the left-hand side becomes a boundary condition at b = 0 and the desired final integral becomes a boundary value of the ODE solution. Hence, defining

$$\mathbf F ( \mathbf b ) = \int \frac { d ^ { 2 } h } { ( 2 \pi ) ^ { 2 } } e ^ { i \mathbf b \cdot \mathbf h } \mathbf F ( \mathbf h ) \,,$$

we have

$$\text{Re} \, \int \frac { d ^ { 2 } h } { ( 2 \pi ) ^ { 2 } } 2 \cdot \mathbf F ( \mathbf q ) = \text{Im} \left ( 2 \nabla _ { \mathbf b } \cdot \mathbf F ( 0 ) \right ),$$

and Eq. (3.11) becomes

$$- 2 i \nabla \delta ^ { 2 } ( \mathbf b ) & = \frac { i } { 2 p \omega ( p - \omega ) } \left ( p ( p - \omega ) m _ { \infty } ^ { 2 } \omega + p \omega \, m _ { \infty } ^ { 2 } \, _ { p - \omega } - \omega ( p - \omega ) m _ { \infty } ^ { 2 } \, _ { p } - \nabla _ { \mathbf b } ^ { 2 } \right ) \mathbf F ( \mathbf b ) \\ & \quad + \left ( \mathcal { C } _ { R } ^ { \prime } ( | \omega | \, b ) - \frac { \mathcal { C } _ { A } ^ { \prime } ( | \omega | \, b ) } { 2 } + \frac { \mathcal { C } _ { A } ^ { \prime } ( | p | \, b ) } { 2 } + \frac { \mathcal { C } _ { A } ^ { \prime } ( | p - \omega | b ) } { 2 } \right ) \mathbf F ( \mathbf b ),$$

with

$$\mathcal { C } _ { R } ^ { \prime } ( | \omega | \, b ) \equiv \int \frac { d ^ { 2 } { \mathbf k } _ { \perp } } { ( 2 \pi ) ^ { 2 } } \left ( 1 - e ^ { i \omega b \cdot k _ { \perp } } \right ) \mathcal { C } _ { R } ( { \mathbf k } _ { \perp } ) \,.$$

As noted above, we need to incorporate the change in the collision operator,

$$\mathcal { C } ^ { \prime } _ { R, \text{LO} + \text{NLO} } ( b ) = \mathcal { C } ^ { \prime } _ { R } ( b ) + \delta \mathcal { C } ^ { \prime } _ { R } ( b ),$$

where the q ⊥ -space correction of C R is the unintegrated version of δq ˆ a detailed in Eq. (3.23). The details appear in Appendix C of Ref. [38]. We must also incorporate a shift in the thermal masses:

$$m _ { \infty \, a, \mathrm L O + \text{NLO} } ^ { 2 } = m _ { \infty \, a } ^ { 2 } + \delta m _ { \infty \, a } ^ { 2 }, \quad \delta m _ { \infty q } ^ { 2 } = - g C _ { \mathrm F } \frac { T m _ { D } } { 2 \pi }, \quad \delta m _ { \infty g } ^ { 2 } = - g C _ { A } \frac { T m _ { D } } { 2 \pi }. \quad ( 3. 3 2 )$$

These shifts are incorporated into Eq. (3.29) by solving for their effects perturbatively, by treating F b ( ) formally as an expansion in powers of δm ∞ , δ C ; F b ( ) = F 0 ( b ) + F 1 ( b ) + . . . to first order. For details on the boundary conditions and the numerical evaluation of these equations, one should refer to [38-40]

The above analysis considers O ( g ) corrections to the regime in which the external states are collinear up to O ( g 2 ) corrections. This allowed for simplifications in the analysis, in particular in expressing Eq. (3.13) in terms of the leading mass and transverse-momentum corrections. But this is not always the case. Integrating over opening angles, Eq. (3.13)

breaks down if the external momenta P, K obey P · K ∼ gT 2 , and this regime represents an O ( g ) part of the splitting processes considered. Similarly, for hard scattering, if the transfer momentum Q becomes nearly lightlike, our treatment also breaks down, and this again represents an O ( g ) effect. Both regimes feature small angles, but not as small as in the collinear splitting processes; so we refer to them as semi-collinear . The analysis of this kinematic region was carried out in Ref. [33], and we will follow the development there and in Ref. [15].

The good news is that, in treating this regime, the kinematic variable h becomes large, namely h ≫ gT 2 . The collision operator then becomes subdominant, C ≪ δE , and one can perturb in it to first order. The lowest-order and first nontrivial-order solutions are

$$& \mathbf F _ { b 1 } ( \mathbf h ) = 2 \mathbf h / i \delta E ( \mathbf h ) \,, \\ & \mathbf F _ { b 2 } ( \mathbf h ) = \frac { i } { \delta E ( \mathbf h ) } \int \frac { d ^ { 2 } q _ { \perp } } { ( 2 \pi ) ^ { 2 } } \bar { C } ( q _ { \perp } ) \left \{ \left ( C _ { R _ { b } } - \frac { C _ { A } } { 2 } \right ) \left [ \mathbf F _ { b 1 } ( \mathbf h ) - \mathbf F _ { b 1 } ( \mathbf h - k \mathbf q _ { \perp } ) \right ] \\ & \quad + \frac { C _ { A } } { 2 } [ \mathbf F _ { b 1 } ( \mathbf h ) - \mathbf F _ { b 1 } ( \mathbf h + p \mathbf q _ { \perp } ) ] + \frac { C _ { A } } { 2 } [ \mathbf F _ { b 1 } ( \mathbf h ) - \mathbf F _ { b 1 } ( \mathbf q - ( p - k ) \mathbf q _ { \perp } ) ] \right \}. \\ - - \quad \cdot \quad \cdot \quad \cdot \quad - \quad - \quad \cdot \quad - \quad \cdot \quad - \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

This is the same as treating the emission in the Bethe-Heitler limit. However, at large h , q -= δE ∼ gT ∼ q ⊥ , which changes the kinematics of the scattering. Therefore, ¯ ( C q ⊥ ) must be recomputed. A more accurate form for ¯ ( C q ⊥ ) in this regime is [33, 38]

$$\bar { C } _ { N L O } ( q _ { \perp }, \delta E ) = \frac { g ^ { 2 } T m _ { D } ^ { 2 } } { ( q _ { \perp } ^ { 2 } + \delta E ^ { 2 } ) ( q _ { \perp } ^ { 2 } + \delta E ^ { 2 } + m _ { D } ^ { 2 } ) } + \frac { 2 g ^ { 2 } T \delta E ^ { 2 } } { q _ { \perp } ^ { 2 } ( q _ { \perp } ^ { 2 } + \delta E ^ { 2 } ) } \,. \quad \quad ( 3. 3 4 )$$

Physically, this is a consequence of the splitting being induced by an elastic scattering or by the absorption of a soft on-shell particle, as observed in [15]. Therefore, the alreadycomputed LO 1 ↔ 2 contribution, corresponding to the small δE limit, and the alreadyincluded LO 2 ↔ 2 contribution [38] need to be subtracted as:

$$\delta \bar { C } ( q _ { \perp }, \delta E ) = \frac { g ^ { 2 } T m _ { D } ^ { 2 } } { ( q _ { \perp } ^ { 2 } + \delta E ^ { 2 } ) ( q _ { \perp } ^ { 2 } + \delta E ^ { 2 } + m _ { D } ^ { 2 } ) } - \frac { g ^ { 2 } T m _ { D } ^ { 2 } } { q _ { \perp } ^ { 2 } ( q _ { \perp } ^ { 2 } + m _ { D } ^ { 2 } ) } \,.$$

The second term is the LO collinear form for ¯ (the small C δE limit of Eq. (3.34)).

The semi-collinear contribution is found by substituting Eq. (3.35) into Eq. (3.33) and using it to evaluate Eq. (3.11) and finally Eq. (3.9), which results in

$$\sin & \left | \psi _ { b c } ^ { a } \right | _ { \text{semi} } ( p ; p - k, k ) = \frac { g ^ { 2 } } { 3 2 \pi ^ { 4 } } \left [ \begin{array} { c c } d _ { A } C _ { A } \frac { p ^ { 4 } + k ^ { 4 } + ( p - k ) ^ { 4 } } { p ^ { 3 } k ^ { 3 } ( p - k ) ^ { 3 } } & g \leftrightarrow g g \\ d _ { F } C _ { F } \frac { p ^ { 2 } + ( p - k ) ^ { 2 } } { p ^ { 2 } ( p - k ) ^ { 2 } k ^ { 3 } } & q \leftrightarrow q g \end{array} \right ] \left ( \begin{array} { c } d ^ { 2 } h \\ d _ { F } C _ { F } \frac { ( p - k ) ^ { 2 } + k ^ { 2 } } { p ^ { 3 } } & g \leftrightarrow q ^ { 2 } \end{array} \right ] \left ( \begin{array} { c } d ^ { 2 } h \\ d _ { F } C _ { F } \frac { ( p - k ) ^ { 2 } k ^ { 2 } p ^ { 3 } } { p ^ { 2 } ( p - k ) ^ { 2 } k ^ { 2 } p ^ { 3 } } & g \leftrightarrow q ^ { 2 } \end{array} \right ] \\ & \times \left [ \left ( C _ { R } - \frac { C _ { A } } { 2 } \right ) \left ( \frac { h } { \delta E ( h ) } - \frac { h - k \mathbf q _ { \perp } } { \delta E ( h - k \mathbf q _ { \perp } ) } \right ) ^ { 2 } + \frac { C _ { A } } { 2 } \left ( \frac { h } { \delta E ( h ) } - \frac { h } { \delta E ( h + p \mathbf q _ { \perp } } \right ) ^ { 2 } \\ & + \frac { C _ { A } } { 2 } \left ( \frac { h } { \delta E ( h ) } - \frac { h - ( p - k ) \mathbf q _ { \perp } } { \delta E ( h - ( p - k ) \mathbf q _ { \perp } ) } \right ) ^ { 2 } \right ],$$

with ( C R -C / A 2) appearing on the h + p q p term for g ↔ qq ¯ processes. This is then inserted into Eq. (3.9), resulting in

$$( f _ { 1 }, \mathcal { C } ^ { s e m i } f _ { 1 } ) & \equiv \frac { 2 \pi } { T ^ { 3 } } \sum _ { a b c } \int _ { 0 } ^ { \infty } d p \int _ { 0 } ^ { p } d k \, \gamma _ { b c } ^ { a } \Big | _ { s e m i } ( p ; p - k, k ) f _ { 0 } ^ { a } ( p ) [ 1 \pm f _ { 0 } ^ { b } ( k ) ] [ 1 \pm f _ { 0 } ^ { c } ( p - k ) ] \\ & \quad \times \Big [ f _ { 1 } ^ { a } ( p ) - f _ { 1 } ^ { b } ( k \hat { p } ) - f _ { 1 } ^ { c } ( ( p - k ) \hat { p } ) \Big ] ^ { 2 }.$$

## 3.2.3 Summary

In summary, the corrections to be added to the strict leading-order contributions of the previous section produce a collision operator that is fully NLO. At NLO, the collision term will be given as

$$\Phi ^ { \Phi \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Psi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Philon _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi } \Phi _ { \Phi }$$

where

$$\left ( f _ { 1 }, \delta C \, f _ { 1 } \right ) \equiv & \left ( f _ { 1 }, C ^ { \delta \hat { q } } \, f _ { 1 } \right ) - \left ( f _ { 1 }, C _ { \mathcal { O } ( g ) \, \text{finite} } ^ { 2 \leftrightarrow 2 } \, f _ { 1 } \right ) + \left ( f _ { 1 }, C ^ { \text{semi} } \, f _ { 1 } \right ) + \left ( f _ { 1 }, \delta C ^ { 1 \leftrightarrow 2 } \, f _ { 1 } \right ). \ \ ( 3. 3 9 )$$

Here, the first two terms are the contribution from the soft gluon exchange and its counterterm. The last two terms are the contribution from semi-collinear processes and the correction to the splitting rate, respectively. All of these corrections were carefully studied in [15], and what was found is that the main difference between LO and NLO results arises from the soft gluon exchange.

## 4 From zero to nonzero chemical potential

In this section, we delve into the differences between the description of the QGP at vanishing chemical potential, µ = 0, and at high densities, µ &gt; 0. Temperature, density, and screening enter the calculation, and we must make sure that we include them correctly each place they both arise. The structure fits within the general effective kinetic framework of Ref. [24], for which finite µ and T case is a special case.

When introducing a chemical potential we must shift the equilibrium occupancies from their simple thermal values; we have already done this in Eq. (2.4). This change propagates into the Boltzmann equation and the definition of the extremization problem, Eq. (2.14), in a straightforward way. There is no need to change the basis of variational functions, but more terms may be needed because the most important momentum range is likely to be higher, p ∼ µ rather than p ∼ T .

The matrix elements depend on the chemical potential indirectly through the strength of screening effects. As shown in Ref. [41], the form of the HTLs are identical for high density as for high temperature. The two effects simply combine in the Debye screening mass squared,

$$m _ { D } ^ { 2 } = \frac { N _ { c } } { 3 } g ^ { 2 } T ^ { 2 } + \frac { N _ { f } } { 6 } g ^ { 2 } \left ( T ^ { 2 } + \frac { 3 } { \pi ^ { 2 } } \mu ^ { 2 } \right ).$$

The last term represents the increase in the screening due to the quark abundance induced by the chemical potential. The strength of quark screening is similarly enhanced; one uses

the occupancies with chemical potentials in evaluating Eq. (3.4). The resulting expression is already written in that equation.

In evaluating splitting rates, a pivotal role was played by the transverse scattering rate, described by the kernel ¯ ( C q ⊥ ). Remarkably, Eq. (3.12) still holds at finite chemical potential, provided one uses the adjusted value of the Debye mass squared. This is discussed in Ref. [24]. Essentially, the shifted m 2 D value accounts for the change in the number of scattering 'targets' in the medium and therefore correctly adjusts the strength of medium scattering which induces the splitting. Similarly, the calculation of δE for splitting, Eq. (3.13), must use the form of the asymptotic masses which accounts for chemical potentials.

Similarly, at NLO we must correctly account for the strength of scattering in the medium, but generally the calculations presented so far remain true after one does so. The light-cone methods which were used to derive Eq. (3.21) and Eq. (3.23) are based on integrating out the quarks and incorporating their effects in the HTLs or equivalently within the effective theory of EQCD. This proceeds the same with µ as without, except that the value of m 2 D must be adjusted using Eq. (4.1). By similar reasoning, Eq. (3.32) and Eq. (3.34) apply as written, but using the corrected form of m 2 D .

Next, let us briefly analyze the formal convergence of perturbation theory. At weak coupling α s ≪ 1 and µ = 0 hot QCD has three scales: the hard scale πT , the soft scale gT , and the ultrasoft scale g T 2 . At leading order the scale πT participates as the dominant propagating degrees of freedom, while the gT scale is relevant as the momentum-exchange scale for soft scattering relevant in both 2 ↔ 2 and 1 ↔ 2 processes. NLO corrections occur primarily because these gT scale degrees of freedom have occupancies f ( q ) ∼ 1 /g and loop corrections are only suppressed by a single power of g . Therefore the structure of NLO O ( g ) corrections relative to LO is often as a ratio g T/m 2 D . Including color and other numerical factors, these corrections often prove to be quite large.

The inclusion of a chemical potential with µ &gt; πT raises the hard scale from πT to µ , and it increases the Debye scale as one sees in Eq. (4.1): m D ∼ gµ rather than gT . Therefore the separation between the soft and ultrasoft scales is increased, and O ( g ) corrections are expected to scale as g T/gµ 2 instead of g T/gT 2 . This underlies the expectation that NLO effects should be milder in the largeµ/T case.

In practice, the largest NLO effect is the shift in the transverse momentum diffusion coefficient, Eq. (3.23). This effect is linear in m D , while the leading-order part is proportional to m 2 D . The ratio δq/q ˆ ˆ ∼ g m 4 D T /g m 2 2 2 D T is therefore of order g T/m 2 D , as anticipated above. For parametrically large µ/T , one indeed finds an improvement in the NLO convergence by a single factor of T/µ . In the next section we will see to what extent this really improves the convergence in practice.

## 5 Results

In the last sections, we have presented the ingredients necessary for the computation of shear viscosity at leading order, discussed the relevant next-to-leading-order corrections, and presented a good basis set for these calculations. Using this computational apparatus,

we compute the shear viscosity as a function of µ/T at both leading order and (almost) next-to-leading order.

Figure 3 . η/s as a function of the temperature for µ = 0 3 6. , , The dashed lines represent NLO order results, while continuous lines are LO results. Both were calculated using the renormalization point µ EQCD = 2 7 . T

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00524v1-figure-9879d087-e365-44a4-8a5a-dcf3a4aa0174.jpg)

In figure 3 we plot the shear-viscosity to entropy density ratio η/s for LO and NLO, and in figure 4 we show the kinematic viscosity, that is, the ratio of the shear viscosity to the enthalpy density, ηT/ e ( + P ). 3 Note that the chemical potential µ we refer to is that for a quark, and we treat it to be the same for all quark species. The baryon number chemical potential is then 3 times larger. We have used the EQCD-inspired scale setting by choosing the MS renormalization point to be µ EQCD = 2 7 . T , following the treatment from [15], and we run the coupling using the 2-loop beta function and the value α s = 0 118 at the . Z -pole. We assume that when µ EQCD reaches a quark mass, one crosses the quark-mass threshold and the quark becomes active both in the medium and in the running of the coupling constant. So for temperatures below 2 7 . T = m b we use N f = 4 and for temperatures above 2 7 . T = m b we use N f = 5 (and similarly for the charm and top quarks). The plots clearly show that the difference between LO and NLO is smaller at large µ than at µ = 0, particularly at high temperatures where the coupling is smaller and N f is larger.

Since one motivation to study high-density regimes comes from the expectation that perturbation theory should show better convergence than at µ = 0, we present a more detailed comparison between LO and NLO results in Fig. 5. We plot the ratio between LO and NLO calculations as a function of the temperature for each fixed chemical potential. We reproduce the previous result [15] that the NLO corrections at µ = 0 are large, more than a factor of 2, up to 10 3 GeV temperatures. Once one goes to higher densities, the NLO effects become less dominant, falling below a factor of 2 by T = 500 MeV for µ/T = 6.

3 Raw data for all plots are included in the ancillary materials.

Figure 4 . Kinematic shear viscosity as a function of the temperature for µ = 0 3 6. , , The dashed lines represent NLO order results, while continuous lines are LO results. Both were calculated using the coupling µ EQCD = 2 7 . T

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00524v1-figure-0923c6e5-424b-4e9b-8e82-ae8596aae29d.jpg)

Figure 5 . Ratio of the NLO shear viscosity to the LO shear viscosity, indicating the importance of subleading effects, for three chemical potentials over a range of temperatures. The renormalization scale was set at µ EQCD = 2 7 . T in all cases.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00524v1-figure-4bcf9662-14dd-49d7-974a-98532c2b5977.jpg)

This suggests that the expansion in g is better behaved at large µ/T . Nevertheless, NLO effects are rather important at all temperatures and chemical potentials we consider.

## 6 Summary and Discussion

We have extended the investigation of shear viscosity for hot QCD with chemical potentials to both leading and (almost) next-to-leading order. We find that the viscosity-to-entropy ratio η/s is larger at finite µ than at zero µ for all temperatures we consider. However, the kinematic viscosity ηT/ e ( + P ) is more weakly dependent on µ , falling below the µ = 0 result except at NLO below T = 1 GeV, where the µ/T dependence almost vanishes. We argued previously that, at a given temperature or coupling strength, the convergence of perturbation theory should be better in the large µ/T regime than for µ = 0. Our results show that this is indeed the case, but the improvement is not dramatic. NLO effects are a factor of 2 at µ/T = 0 for temperatures of hundreds of GeV, whereas for µ = 6 T they are a factor of 2 at about 500 MeV temperature. Perturbation theory does not converge at experimentally achievable combinations of ( µ, T ), but it comes closer at high µ than at small µ .

It would be straightforward to apply the techniques developed here to the more physical case in which different quark species experience different chemical potentials. This could allow the development of an interpolating estimate for the viscosity based on experimental data, hadron-resonance modeling, and NLO perturbative input, similar to what was recently developed by one of us [16].

## Acknowledgements

The authors acknowledge the support by the State of Hesse within the Research Cluster ELEMENTS (Project ID 500/10.006), and support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through the CRC-TR 211 'Strong-interaction matter under extreme conditions'- project number 315477589 - TRR 211. The authors also thank Nicolas Wink, Jacopo Ghiglieri, and Kevin Ingles for instructive conversations.

## A Self-energies beyond the hard thermal loop approximation

As discussed in [24], the t-channel exchange process is subtle, since the use of vacuum matrix elements leads to logarithmic divergences in the 'transport' cross-section and in the inverse viscosity. Baym et al already identified the correct procedure to fix this problem in the context of viscosity [25]; one should amend the matrix element through the inclusion of the thermal self-energy. At the leading-log level, relevant in Refs. [18, 25], it is sufficient to identify the energy scale where this occurs. At leading order it is necessary to really include the self-energy, at least within an approximation which is adequate in the formal kinematic regime ( p, k ) ∼ T but ( q, ω ) ∼ gT ≪ T . This is precisely the regime where the HTL approximation should be valid.

At large chemical potentials we can try to do better. In this regime, the self-energy is dominated by fermionic loops. For these, the self-energy is gauge invariant and the full expression, not just the hard-loop limit, is known [31, 32]. Specifically, the medium-induced

part of the fermion-loop contribution to the gluonic self-energy is [31, 32]

$$\Pi _ { 0 0 } = \frac { g ^ { 2 } N _ { f } T ^ { 2 } } { 3 } H ( \omega, q ), \quad \Pi _ { T } = \frac { g ^ { 2 } N _ { f } T ^ { 2 } } { 3 } \left ( - \frac { 1 } { 2 } H ( \omega, q ) + \frac { q ^ { 2 } } { 2 ( q ^ { 2 } - \omega ^ { 2 } ) } G ( \omega, q ) \right ) \quad ( \text{A.1} )$$

where H ω,q ( ) and G ω,q ( ) are, writing 2 ω + = ω + q and 2 ω -= ω -q :

$$\text{Re} \, G ( \omega, q ) = \frac { 3 } { \pi ^ { 2 } T ^ { 2 } } \int _ { 0 } ^ { \infty } d k f ( k ) \left ( 4 k + \frac { q ^ { 2 } - \omega ^ { 2 } } { 2 q } \ln \left [ \frac { ( k + \omega _ { - } ) ( k - \omega _ { + } ) } { ( k - \omega _ { - } ) ( k + \omega _ { + } ) } \right ] \right ) \quad \text{.}$$

$$\intertext { w e } \cup ( \omega, q ) & = \frac { 3 } { \pi ^ { 2 } T ^ { 2 } } \int _ { 0 } ^ { u \kappa _ { J } \left ( \kappa \right ) } \left \lfloor \kappa + \frac { 2 q } { 2 q } \right \| \left \lfloor \overline { ( k - \omega _ { - } ) ( k + \omega _ { + } ) } \right \rfloor \right ) \\ \text{Re} \, H ( \omega, q ) & = \frac { 3 } { \pi ^ { 2 } T ^ { 2 } } \int _ { 0 } ^ { \infty } d k f ( k ) \left ( 2 k - \frac { 4 k ^ { 2 } + \omega ^ { 2 } - q ^ { 2 } } { 4 q } \ln \left [ \frac { ( k + \omega _ { - } ) ( k - \omega _ { + } ) } { ( k - \omega _ { - } ) ( k + \omega _ { + } ) } \right ] \\ & - \frac { 2 k \omega } { q } \ln \frac { \omega _ { + } } { \omega _ { - } } - \frac { k \omega } { q } \ln \left [ \frac { ( k + \omega _ { - } ) ( k - \omega _ { + } ) } { ( k - \omega _ { - } ) ( k + \omega _ { + } ) } \right ] \right ) \\ \text{Im} \, G ( \omega, q ) & = \frac { 3 } { \pi ^ { 2 } T ^ { 2 } } \int _ { 0 } ^ { \infty } d k f ( k ) \frac { q ^ { 2 } - \omega ^ { 2 } } { 2 q } [ \Theta ( k + \omega _ { + } ) \Theta ( - k - \omega _ { - } ) - \Theta ( - k + \omega _ { + } ) \Theta ( k - \omega _ { - } ) ] \\ \text{Im} \, H ( \omega, q ) & = \frac { 3 } { \pi ^ { 2 } T ^ { 2 } } \int _ { \Lambda } ^ { \infty } d k f ( k ) \left [ - \frac { ( 2 k + \omega ) ^ { 2 } - q ^ { 2 } } { 4 a } \Theta ( k + \omega _ { + } ) \Theta ( - k - \omega _ { - } )$$

$$\frac { 3 } { \pi ^ { 2 } T ^ { 2 } } & \int _ { 0 } ^ { \infty } d k f ( k ) \left ( 4 k + \frac { q ^ { 2 } - \omega ^ { 2 } } { 2 q } \ln \left [ \frac { ( k + \omega _ { - } ) ( k - \omega _ { + } ) } { ( k - \omega _ { - } ) ( k + \omega _ { + } ) } \right ] \right ) \\ & 3 \quad \int _ { \ n, \ e t r \searrow \infty } ^ { \infty } 4 k ^ { 2 } + \omega ^ { 2 } - q ^ { 2 } \, _ { \iota _ { - } } \left [ ( k + \omega _ { - } ) ( k - \omega _ { + } ) \right ]$$

$$\text{Im} \, H ( \omega, q ) & = \frac { \overset { n } { 3 } } { \pi ^ { 2 } T ^ { 2 } } \int _ { 0 } ^ { \infty } d k f ( k ) \left [ - \frac { ( 2 k + \omega ) ^ { 2 } - q ^ { 2 } } { 4 q } \Theta ( k + \omega _ { + } ) \Theta ( - k - \omega _ { - } ) \\ & + \frac { ( 2 k - \omega ) ^ { 2 } - q ^ { 2 } } { 4 q } \Theta ( k - \omega _ { - } ) \Theta ( - k + \omega _ { + } ) + \frac { 2 k \omega } { q } \Theta ( \omega _ { + } ) \Theta ( - \omega _ { - } ) \right ] \\ \text{Im} \, \quad \, \quad \, \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \end{cases}$$

The references only present the real parts, but the imaginary parts are easily obtained by using the retarded prescription ω → ω + 0 and carefully interpreting the logarithms. i The hard thermal loop approximation corresponds to expanding the k integrals in the approximation k ≫ ( q, ω ), in which case the functions simplify:

$$H _ { H T L } ( \omega / q ) = 1 - \frac { \omega } { 2 q } \ln \frac { \omega _ { + } } { \omega _ { - } } + i \frac { \pi \omega } { 2 q } \Theta ( \omega _ { + } ) \Theta ( - \omega _ { - } )$$

$$G _ { H T L } ( \omega / q ) & = 1.$$

One concern regarding the AMY leading-order calculation [17] has always been the use of hard thermal loop self-energies, rather than full thermal self-energies. At high chemical potential, where the hard loop is dominated by fermions, the above self-energies accurately represent the correct loop screening effect at any energy and momentum, without the need for an HTL approximation. Therefore we can find out how well the HTL approximation works by comparing how the leading-order shear viscosity changes when we 'improve' the screening effects by replacing the HTL approximation, Eq. (A.3), with the full self-energy, Eq. (A.2).

We make this comparison in figure 6. The figure considers a single collision process, diagram (C), and computes the collision integral ( χ ij , C χ ij ) for the single-parameter Ansatz. We then take the ratio of this collision rate for the case of the full self-energy over the collision rate using the HTL self-energy. The Debye screening strength is chosen to vary with µ as it would if the coupling strength g 2 is held fixed. The figure shows, first and foremost, that the effect of replacing the HTL self-energy with the full self-energy is not that large. Also, the HTL approximation actually becomes more reliable as the chemical potential is increased. This indicates that the use of the HTL approximation for the self-energies is justified and represents a smaller effect than the NLO (and, presumably, unknown NNLO) corrections.

Figure 6 . Ratio of diagram (C), integrated over the 1-parameter Ansatz, for the full vs the HTL self-energy. The Debye mass is varied such that the coupling g 2 is kept fixed.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00524v1-figure-effa81bf-6862-4b66-91dc-d9540ef4fb19.jpg)

## References

- [1] J. Grefa, M. Hippert, J. Noronha, J. Noronha-Hostler, I. Portillo, C. Ratti et al., Transport coefficients of the quark-gluon plasma at the critical point and across the first-order line , Phys. Rev. D 106 (2022) 034024 [ 2203.00139 ].
- [2] PHENIX collaboration, Formation of dense partonic matter in relativistic nucleus-nucleus collisions at RHIC: Experimental evaluation by the PHENIX collaboration , Nucl. Phys. A 757 (2005) 184 [ nucl-ex/0410003 ].
- [3] STAR collaboration, Experimental and theoretical challenges in the search for the quark gluon plasma: The STAR Collaboration's critical assessment of the evidence from RHIC collisions , Nucl. Phys. A 757 (2005) 102 [ nucl-ex/0501009 ].
- [4] ALICE collaboration, Elliptic flow of charged particles in Pb-Pb collisions at 2.76 TeV , Phys. Rev. Lett. 105 (2010) 252302 [ 1011.3914 ].
- [5] U. Heinz and R. Snellings, Collective flow and viscosity in relativistic heavy-ion collisions , Ann. Rev. Nucl. Part. Sci. 63 (2013) 123 [ 1301.2826 ].
- [6] M. Luzum and H. Petersen, Initial State Fluctuations and Final State Correlations in Relativistic Heavy-Ion Collisions , J. Phys. G 41 (2014) 063102 [ 1312.5503 ].
- [7] JETSCAPE collaboration, Multisystem Bayesian constraints on the transport coefficients of QCD matter , Phys. Rev. C 103 (2021) 054904 [ 2011.01430 ].
- [8] M. Laine, What is the simplest effective approach to hot QCD thermodynamics? , in 5th Internationa Conference on Strong and Electroweak Matter , pp. 137-146, 2003, hep-ph/0301011 , DOI.
- [9] A. Hietanen, K. Kajantie, M. Laine, K. Rummukainen and Y. Schroder, Three-dimensional physics and the pressure of hot QCD , Phys. Rev. D 79 (2009) 045018 [ 0811.4664 ].

- [10] V. A. Goy, V. Bornyakov, D. Boyda, A. Molochkov, A. Nakamura, A. Nikolaev et al., Sign problem in finite density lattice QCD , PTEP 2017 (2017) 031D01 [ 1611.08093 ].
- [11] STAR collaboration, Beam energy dependence of moments of the net-charge multiplicity distributions in Au+Au collisions at RHIC , Phys. Rev. Lett. 113 (2014) 092301 [ 1402.1558 ].
- [12] HADES collaboration, Probing dense baryon-rich matter with virtual photons , Nature Phys. 15 (2019) 1040.
- [13] V. Friese, The CBM experiment at GSI/FAIR , Nucl. Phys. A 774 (2006) 377.
- [14] I. Danhoni and G. D. Moore, Hot and dense QCD shear viscosity at leading log , JHEP 02 (2023) 124 [ 2212.02325 ].
- [15] J. Ghiglieri, G. D. Moore and D. Teaney, QCD Shear Viscosity at (almost) NLO , JHEP 03 (2018) 179 [ 1802.09535 ].
- [16] I. Danhoni, J. S. S. Martin and J. Noronha-Hostler, Shear viscosity from perturbative Quantum Chromodynamics to the hadron resonance gas at finite baryon, strangeness, and electric charge densities , 2406.04968 .
- [17] P. B. Arnold, G. D. Moore and L. G. Yaffe, Transport coefficients in high temperature gauge theories. 2. Beyond leading log , JHEP 05 (2003) 051 [ hep-ph/0302165 ].
- [18] P. B. Arnold, G. D. Moore and L. G. Yaffe, Transport coefficients in high temperature gauge theories. 1. Leading log results , JHEP 11 (2000) 001 [ hep-ph/0010177 ].
- [19] L. D. Landau and E. M. Lifshitz, Fluid Mechanics, Second Edition: Volume 6 (Course of Theoretical Physics) , Course of theoretical physics / by L. D. Landau and E. M. Lifshitz, Vol. 6. Butterworth-Heinemann, 1987.
- [20] P. B. Arnold, C. Dogan and G. D. Moore, The Bulk Viscosity of High-Temperature QCD , Phys. Rev. D 74 (2006) 085021 [ hep-ph/0608012 ].
- [21] S. Jeon, Hydrodynamic transport coefficients in relativistic scalar field theory , Phys. Rev. D 52 (1995) 3591 [ hep-ph/9409250 ].
- [22] J.-S. Gagnon and S. Jeon, Leading order calculation of electric conductivity in hot quantum electrodynamics from diagrammatic methods , Phys. Rev. D 75 (2007) 025014 [ hep-ph/0610235 ].
- [23] J.-S. Gagnon and S. Jeon, Leading Order Calculation of Shear Viscosity in Hot Quantum Electrodynamics from Diagrammatic Methods , Phys. Rev. D 76 (2007) 105019 [ 0708.1631 ].
- [24] P. B. Arnold, G. D. Moore and L. G. Yaffe, Effective kinetic theory for high temperature gauge theories , JHEP 01 (2003) 030 [ hep-ph/0209353 ].
- [25] G. Baym, H. Monien, C. J. Pethick and D. G. Ravenhall, Transverse Interactions and Transport in Relativistic Quark - Gluon and Electromagnetic Plasmas , Phys. Rev. Lett. 64 (1990) 1867.
- [26] S. Jeon and L. G. Yaffe, From quantum field theory to hydrodynamics: Transport coefficients and effective kinetic theory , Phys. Rev. D 53 (1996) 5799 [ hep-ph/9512263 ].
- [27] E. Braaten and R. D. Pisarski, Soft Amplitudes in Hot Gauge Theories: A General Analysis , Nucl. Phys. B 337 (1990) 569.
- [28] J. Frenkel and J. C. Taylor, High Temperature Limit of Thermal QCD , Nucl. Phys. B 334 (1990) 199.

- [29] A. Kurkela and A. Mazeliauskas, Chemical equilibration in weakly coupled QCD , Phys. Rev. D 99 (2019) 054018 [ 1811.03068 ].
- [30] G. D. Moore, Transport coefficients in large N(f) gauge theory: Testing hard thermal loops , JHEP 05 (2001) 039 [ hep-ph/0104121 ].
- [31] V. V. Klimov, Collective Excitations in a Hot Quark Gluon Plasma , Sov. Phys. JETP 55 (1982) 199.
- [32] H. A. Weldon, Covariant Calculations at Finite Temperature: The Relativistic Plasma , Phys. Rev. D 26 (1982) 1394.
- [33] J. Ghiglieri, G. D. Moore and D. Teaney, Jet-Medium Interactions at NLO in a Weakly-Coupled Quark-Gluon Plasma , JHEP 03 (2016) 095 [ 1509.07773 ].
- [34] J. Ghiglieri and D. Teaney, Parton energy loss and momentum broadening at NLO in high temperature QCD plasmas , Int. J. Mod. Phys. E 24 (2015) 1530013 [ 1502.03730 ].
- [35] G. D. Moore and D. Teaney, How much do heavy quarks thermalize in a heavy ion collision? , Phys. Rev. C 71 (2005) 064904 [ hep-ph/0412346 ].
- [36] S. Caron-Huot, O(g) plasma effects in jet quenching , Phys. Rev. D 79 (2009) 065039 [ 0811.1603 ].
- [37] P. Aurenche, F. Gelis and H. Zaraket, A Simple sum rule for the thermal gluon spectral function and applications , JHEP 05 (2002) 043 [ hep-ph/0204146 ].
- [38] J. Ghiglieri, J. Hong, A. Kurkela, E. Lu, G. D. Moore and D. Teaney, Next-to-leading order thermal photon production in a weakly coupled quark-gluon plasma , JHEP 05 (2013) 010 [ 1302.5970 ].
- [39] I. Ghisoiu and M. Laine, Interpolation of hard and soft dilepton rates , JHEP 10 (2014) 083 [ 1407.7955 ].
- [40] J. Ghiglieri and G. D. Moore, Low Mass Thermal Dilepton Production at NLO in a Weakly Coupled Quark-Gluon Plasma , JHEP 12 (2014) 029 [ 1410.4203 ].
- [41] E. Braaten and R. D. Pisarski, Simple effective Lagrangian for hard thermal loops , Phys. Rev. D 45 (1992) R1827.
Isabella Danhoni Guy D. Moore
Hot and Dense QCD Shear Viscosity at (almost) NLO
The next-to-leading order weak-coupling shear viscosity of QCD was computed 6
years ago.
  However, these results have never been applied at finite baryon chemical
potential $\mu$, even though intermediate-energy heavy ion collisions and
merging neutron stars may explore the Quark-Gluon Plasma in a regime where
baryon chemical potentials are large.
  Here, we extend the next-to-leading order shear viscosity calculations to
finite $\mu$, and we show that, while the convergence of the weak-coupling
expansion is questionable for achievable plasmas, it is somewhat better at $\mu
> T$ than at $\mu=0$.