Non-perturbative Suppression of Chiral Vortical Effect in Hot (s)QGP for Hyperons Spin Polarization in Heavy Ion Collisions

With the Field Correlator Method (FCM) for QCD, we show that the Chiral Vortical Effect (CVE) in hot (strongly-interacting) Quark-Gluon Plasma ((s)QGP) is modified by non-perturbative interactions: by Color-Magnetic confinement, and by remnant Color-Electric interaction, which is encoded in the Polyakov line. The obtained result demonstrates numerical suppression of CVE comparable to the phenomenological suppression used for numerical simulations of RHIC-STAR data on hyperons spin polarization in non-central heavy ion collision (HIC). The parameters range in the temperature - quark chemical potential plane is expected to cover ALICE and RHIC data. The chiral current is calculated for the rigidly rotating model of (s)QGP in the linear order in angular velocity at the rotation axis with account of non-perturbative interactions.


I. INTRODUCTION
Hyperon spin polarization is a relatively new but a fine probe of strong-interaction theory [1][2][3][4][5][6][7][8][9].The chiral current in HIC leads to hyperons spin polarization asymmetry [2,3] (the Chiral Vortical Effect dominates, and the Chiral Separation Effect (CSE) [10][11][12][13][14][15][16][17], which is due to magnetic field, is at least an order of magnitude weaker [2][3][4][5][6][7][8]) Π µ is the spin polarization four-vector (values of its components depend on the choice of a reference frame).N Λ is the number of hyperons with transverse to the reaction plane momentum component p y (m Λ = 1116 MeV is the hyperons mass), v is the hydrodynamic velocity of the (s)QGP, and γ = (1 − v 2 ) −1/2 .µ s is the strange quarks chemical potential, which is conjugated to the net strange charge Q s (conserved up to absolutely negligible effects of the weak interaction).Since the initial state of the system of colliding ions contains only nucleons, the net strangeness of the system is zero at any moment.However, the strange quarks pairs production cross-section q q → ss, gg → ss is considerable during the collision and the (s)QGP evolution, so spatial fluctuations of the strange chemical potential in the fireball lead to the experimentally observable effect.
In the studies [3][4][5] the free value of the chiral vorticity coefficient c V = T 2 /6 + µ 2 /(2π 2 ) was used, which led to overestimation of the result by an order of magnitude in comparison to the signal found in RHIC-STAR data.The suppression of c V was attributed to the correlation effects in (s)QGP, and accounted for by introduction of a constant suppression factor.
In this paper we account for the leading nonperturbative interactions in (s)QGP with the Field Correlator Method, and obtain the formula c V = I 0 (T, µ) (25), which happens to provide a comparable suppression.
As a model of a relatively small (s)QGP grain that emerges in HIC in which the density fluctuation is seen as uniform, we consider rigidly rotating QCD in the deconfined phase -at temperatures above the crossover transition T c ∼ 160 MeV and relatively small quark densities (at least less than the nuclear density).Such a model is consistent with the HIC simulations, where graining (discretization in space and time) is usually larger than 1 fm(/c).
Since the mean free path l < T −1 [18][19][20][21] of a quark in (s)QGP is few times smaller than the typical size of the vortical structure [5,9,22], the rigid rotation may be a satisfactory model for the (s)QGP motion -the angular velocity corresponds to the vorticity at any given point inside the fireball.
The parameter range is expected to cover ALICE and RHIC data, see Fig. 1 [14,[23][24][25][26][27][28][29][30][31][32][33].SPS, NICA and FAIR might require a separate assessment.The axial current of quarks at the rotation axis is calculated in the leading order in the strong coupling and in the linear order in the angular velocity.Thus, only the non-perturbative interactions are accounted for, and the dynamical effect of rotation on the gluonic content of (s)QGP is neglected, which is an intriguing problem on its ownrigidly rotating Gluon Plasma was recently reported [34] to be unstable at temperatures T < 1.5T c where T c is the critical temperature of the Gluonic Plasma.However, in ongoing experiments, the vorticity of (s)QGP is at most ∼ 10 22 s −1 ∼ 10 MeV [1], which is a small scale in the FCM (the QCD vacuum correlation length is λ ∼ 1 GeV −1 [35]), so we expect the dependence of the field correlators on the angular velocity to be negligible.
As for the other regions of the phase diagram, in the hadronic phase we expect the CVE to vanish.Below the crossover temperature T c at relatively small densities µ ∼ T the quarks are confined.Bound states, below and above deconfinement, are heavy, and have spins of the constituent quarks correlated, so they barely contribute to the axial current (a small parameter that might describe the suppression is the ratio of temperature to the bound states mass splitting T /|m ρ − m π |).At large densities (µ ∼ µ qc ∼ 0.6 GeV the Field Correlators are less studied [33,36], so there is space for phenomenology.Even with the most optimistic assumption, that µ does not affect the gauge-invariant gluonic field strength correlators (except for suppressing the confining correlator D E ), for µ < M gb ∼ 1.5 GeV, analysis of bound states contributions might be mandatory.
This paper follows the outline of the recent paper on the Chiral Separation Effect [37] (where instead of rotation, a uniform magnetic field was considered).In the next Section II we calculate the axial current in the rigidly rotating hot QCD with account of nonperturbative interactions (a discussion of the gluonic background rotation may be found in the Appendix A), then we visualize the results in the Section III, and briefly sum up the results in the Section IV.

II. AXIAL CURRENT IN THE (THERMAL) WORLD-LINE FORMALISM
In this section we generally follow the second part of [37], and apply the powerful Method of Field Correlators (FCM) to deal with the (dominating even in the deconfined phase) non-perturbative interactions of QCD -for a review, see [35] and references therein.
We start with the definition of the axial current pro-duced by a given quark flavor, where all the perturbative corrections are neglected.In particular, the neglect of the fermionic determinant constrains the calculation to the Single Line Approximation (SLA).Intuitively, CVE is saturated by fermions with uncorrelated spin degrees of freedom, which are accounted for in SLA.tr c,D is the trace over color and Dirac spinor indices, respectively.⟨. ..⟩ stands for averaging over the thermodynamic ensemble with temperature T > T c (β = T −1 ) above the deconfinement transition crossover temperature T c ∼ 160 MeV, at non-zero baryon density that is defined with the given quark flavor chemical potential µ well below the critical chemical potential µ c = 0.6 GeV [33], in the gluonic background field B (B µ = B a µ t a , t a are the generators of the SU(3) algebra in fundamental representation).
The Matsubara formalism for free rigidly rotating fermions was derived in the original paper on CVE [18].The rotating propagator has a form of the free propagator, but 'twisted' with the rotation operator by the angle of rotation gained during the propagation time However, the rotational symmetry for a given quark trajectory may be broken by the background gluonic field.The symmetry is restored after the averaging, since the gluonic ensemble is isotropic.
Let us consider the quark propagator in the compactified Euclidean space, augmented with the parallel transporter (in the background field) up to the closed loop, and averaged over the background field where αγ are the spinor indices.Then, the proof of [18] applies verbatim for the (non-physical states) propagator S.
For a more concise prove of the form of S Ω (x, y), let us consider rotation around the x 3 -axis in the lab frame, and apply the rotation operator R3 (ϕ) = exp(ϕ Ĵ3 ).S is shift-invariant in temporal and angular directions, or, equivalently, the respective momentum conservation laws are fulfilled In this paper the current density on the rotation axis is calculated in the linear order in Ω, thus only the last operator, generated by the quark spin operator, matters.The quark trajectory on a closed path winds up on the temporal direction, so the parallel transporter reduces to the color trace, and the difference in the Euclidean time is proportional to the number of windings for a given trajectory.Dependence of the field correlators on the rotation is neglected, since the angular velocity (∼ 10 MeV [1]) is negligible on the QCD vacuum correlation length scale (λ ∼ 1 GeV −1 [35]).The 'kinematic' effect of the gluonic content rotation is estimated in Appendix A.
To calculate the average, we apply the world-line formalism -Fock-Feynman-Schwinger representation (FFSR), for quark motion, which is the standard way within FCM to separate the quark kinematics and the gluonic field dynamics.In the "quenched approximation" (with the fermion determinant dropped, thus perturbative quark loops disregarded) where accounts for the quark spin degrees of freedom, the world-line kinetic term is corresponds to the standard dimensional regularization d = 4 − ϵ with arbitrary energy scale ν [38].In the leading-order (LO) we disregard the perturbative gluons a as well A → B. As expected, the angular velocity appears in the expression "as the chemical potential for the angular momentum." The path integral (D 4 z) s xy describes the fermion motion from point x to y in world-line (proper) time s.The bar indicates the anti-periodic boundary conditions for the fermion as the trajectory wraps n times around the temporal direction in the Euclidean space R 3 × S 1 with the compactified temporal direction of length β = T −1 ).We used the plain discretization in order to work with the path integrals following [39] (D 4 z) s xy = lim The derivative acts on the path integral and affects quark kinematics as well as the Wilson loop endpoints.
As follows from the very definition of the parallel transporter (its closure forms the loop) for an arbitrary gauge field which we apply to the loop According to the Leibniz rule, the usual derivative ∂ λ also acts on the squared quark propagator, but not on the gauge factors.
The γ 5 properties simplify the Dirac trace, since the leading order in external field is dominant and the Dirac trace is provided by tr D γ 5 γ µ γ ν γ λ γ ρ = −4ϵ µνλρ (mind that in Euclidean space we use the convention {γ µ , γ ν } = 2δ µν , ϵ 1234 = 1).The term proportional to the quark mass m can not contribute -it involves an odd number of Dirac matrices in the trace and, therefore, vanishes in any case.The leading order in the angular velocity is provided by the first nonzero term in the 'rotation operator' expansion in power series To calculate the Dirac trace, we note that the gluonic FCs ⟨. ..⟩B have no preferred spatial direction (see Appendix A for a more thorough analysis), nor net color charge that could produce non-zero average color-magnetic field.
So the Dirac indices are fixed (µλ[ρσ] = 34 [12]) by the angular velocity direction At this point we apply (6).The usual derivative ∂ 4 also acts on the squared quark propagator as As a given quark trajectory z (n) wraps n times around the temporal direction, z 4 (s) − z 4 (0) = nβ.The Wilson loop in the deconfined phase approximately factorizes [40] (the correlator D EH 1 is neglected; D EH 1 was found so far to yield no pronounced physical effects, so the approximation seems to be numerically accurate) into a spatial loop and an n-times wrapped Polyakov line (in a certain, 'physical,' normalization, in which the additive constant in the Coulomb potential is equal to zero) The path integral factorizes exactly.The first factor -the spatial Wilson loop, provides the Color-Magnetic Confinement -the area law for spatial projection of the quark trajectory ⃗ z(τ ) with CM string tension σ H,f and the minimal area S 3 [⃗ z(τ )].CMC effectively suppresses large deviations of the trajectory z (n) (τ ) from the static trajectory ⃗ z(τ ) ∼ const, thus makes the second factor close to the Polyakov line (in power of winding number absolute value, since the winding direction is irrelevant for vacuum averaging).The approximation , where λ is the QCD vacuum correlation length [35].
The Polyakov line is normalized in a way that its po- is the energy required to overcome the remnant interaction that bounds the quark as a part of a color-singlet state.The potential was not calculated within FCM so far, so we rely on lattice data and use the fit [36] V V 1 (T c ) = 0.5 GeV, T c = 160 MeV.
The potential depends on quark chemical potential via quark loops, but, as a baseline approximation, we disregard the dependence.This allows us to evaluate the temporal path integral The spatial projection of the Wilson loop ) is defined by CMC, where S 3 is the area of the minimal area surface bounded by ⃗ z(τ ) and σ H is the Color-Magnetic string tension (for quarks -in fundamental representation).The Euclidean time dependence in the first factor of ( 10) is unimportant, because the spatial, or Color-Magnetic, string tension σ H is produced by short-distance correlations of the vacuum field B. CMC results in the Debye-like screening for quarks (and gluons) [40][41][42] where M is the screened quark mass.However, the analytical form ( 15) is a rough approximation that, on the other hand, allows us to proceed and obtain a manageable final result.With this approximation, the CVE suppression is overestimated.c D ≈ 2 and c σ ≈ 0.56 are numerical constants.c D ≈ 2 was calculated in the original paper [41].c σ was extracted from lattice data in [41], and has recently been calculated within FCM [43].The numerical constant entering the thermal running of the strong coupling, L σ ≈ 0.1, which we utilize in the one loop approximation, was extracted from lattice data in [41] CMC grows with temperature, which makes the Backgound Perturbation Theory for the thermal QCD selfconsistent in higher orders of perturbation theory: the thermal FCM (CMC) resolves the Linde problem [44].Meanwhile, the standard perturbation theory (the Hard Thermal Loop) provides no adequate suppression for high order diagrams.In particular, the perturbative Debye screening is weak in comparison to lattice data.On the other hand, account of CMC allowed the FCM to move from qualitative to quantitative agreement [42] on QCD thermodynamics with lattice data.Of course, the field correlators and the string tension depend on baryon density via string tension renormalization.Qualitatively, we expect the string tension to decrease with chemical potential.
However, CM interaction yields a considerable negative contribution to the effective quark mass via (nonperturbative) self-energy correction.Instead of direct calculation of the correction we adapt the result of [45] for the squared quark propagator G (in the non-perturbative gluonic field; of course, the self-energy is always a part of a gauge-invariant expression) to the deconfined phase (σ : F ≡ σ µν F µν ).Since the QCD vacuum correlation length fulfills λ ∼ 1 GeV −1 ≪ β (in the old paper [45] the length is denoted as T g ) in the temperature range T c < T < 3T c , T c ≈ 160 MeV, of our interest the integral for the non-perturbative self energy converges within one winding.Also, the current quark masses for the light flavors are small in comparison with the inverse correlation length.Thus we approximate the exact squared propagator G in (18) with the free scalar propagator.
In the deconfined phase the Color-Electric confining correlator is absent while the Color-magnetic is present (as explained near ( 16)), thus we disregard the non-perturbative part of Color-Electric correlator (to be specific, we neglect all the correlators D E,H,EH The consideration of [45] is applicable up to the overall spin-averaging factor: the factor σ µν σ µν = D(D − 1)/4 in the confined phase is to be replaced with σ ij σ ij = (D − 1)(D − 2)/4.Thus, the quark mass shift ∆m 2 q in the (s)QGP phase is twice smaller then in the hadronic phase.
Finally we conclude that in (s)QGP the shift and the resulting Debye-screened quark mass are The correction reduces the screened quark mass ( 16) by a factor 1 − 2 π ∼ 1 3 .The FC part of the problem is sorted out The non-winding (n = 0) trajectories drop out, so the regularization ζ(s) is no longer needed.Then we sequentially use the integral representations for the modified Bessel functions and their properties and substitute p = M sinh t To sum the series we reduce it to geometrical progression, which yields the Fermi-Dirac distribution where V 1 (T ) and M (T ) were defined in ( 13) and ( 19).The result is diagonal in flavor indices.

III. NUMERICAL RESULTS
Eq. ( 25) in the non-interacting limit (V 1 = 0, M = m = const) reproduces the standard results The massless case was studied in the original paper [18], while especially the latter limit might be formal due to finite size effects, see [18,46,47] for discussion (the temperature T is required to be much larger than the inverse size R −1 of the rigidly rotating system).The obtained chiral current suppression for the strange flavor is presented in Fig. 2. The ratio I 0 /I free ≡ ⟨j s,FCM 5 ⟩ / ⟨j free 5 ⟩ of (25) for the strange flavor (m s = 120 MeV) in N f = 2 + 1 QCD in an 'isospin-symmetric' medium µ ≡ µ s , µ I = 0, to (26) for a free massless fermion CVE is plotted as a function of temperature and the strange chemical potential.
Even if the rotating (s)QGP is stable in the deconfinement transition region, which is questionable [34], the result is qualitative due to heavy dependence on an 'unstable' numerical input -the Polyakov line potential (13).At large temperatures, on the other hand, the result might be quantitative.
Though in the 'realistic' parameter range the effect is suppressed, according to the formula (25) the effect may be enhanced in comparison to the free massless case (26) (for instance, if the CE interaction is neglected, which is The present calculation is only for the leading order, and the rotation of the gluonic content was neglected.Even though the strong coupling in the temperature range of interest might be numerically large α s ∼ 0.3, the non-perturbative background additionally suppresses the perturbative corrections with the non-perturbative small parameter σ H,f λ 2 ∼ 1/5 [35].
The suppression may be systematically overestimated due to the approximations used.Since the CMC string forms at distances above the vacuum correlation length λ ∼ 0.2 fm ≈ 1 GeV −1 , at smaller distances the effective screening is weaker than in our calculations.The same for the approximations -the baryon density might deteriorate the Color-Electric and CMC interactions.
Finally we did not consider finite size effects on our calculations which are irrelevant only for T ≫ R −1 .They imply a shift of the lower integration boundary in (25) such that R −1 ≲ p which results in a correction proportional to the small parameter (T R) −2 .
The accuracy of our calculations relies on assumptions on values of several parameters.On the one hand, we have α s , σλ 2 ≪ 1, while, on the other hand, we confine ourselves in the (µ, T )-plane to the region µ q , T ≪ M gb ∼ 1.5 GeV (more precisely, the chemical potential should be smaller than the nuclear chemical potential) from above and for fixed µ q to R −1 ≪ T c (µ q ) < T < λ −1 ∼ 1 GeV.

IV. CONCLUSIONS
Our calculations show that for the range of temperatures and densities experimentally accessible at ALICE and RHIC, the Chiral Vortical Effect is suppressed by non-perturbative interactions, -by the Color-Magnetic Confinement, which grows with temperature, and by the remnants of the Color-Electric interaction, which are encoded in the Polyakov line potential.
To deal with the strong interactions of QCD, which are dominant even in the deconfined phase, we used the Field Correlator Method [35].The method also relies on an expansion, but instead of the perturbative expansion in powers of α s , the non-perturbative cluster expansion in powers of a (eventually, defined by the gauge group algebraic structure) small parameter σλ 2 ∼ 1/5 (σ is the confining string tension at zero temperature and λ is the QCD vacuum correlation length).
Dependence of the field correlators on the rotation is neglected, since the characteristic angular velocity (∼ 10 MeV [1]) is negligible Ωλ ≪ 1.The 'kinematic' effect of the gluonic content rotation is non-perturbatively suppressed, as was shown in Appendix A.
In the non-perturbative background, the perturbative correction, which is naively of order α s , is additionally suppressed with the non-perturbative parameter σλ 2 .Without the non-perturbative suppression, the usual thermal perturbation theory yields the Linde problem, which is naturally resolved by the Color-Magnetic Confinement in the Background Perturbation Theory within FCM [44].Thus we expect the obtained (leading order) result to be numerically sufficient for analysis of the hyperons polarization.
We expect that the suppression provided by (25) (as compared to the free fermions, see Fig. 2) describes the suppression due to collective effects that was discussed in [3,5] (and accounted for with a phenomenological coefficient).Incorporation of our results to numerical simulations of HIC would provide an ultimate check.

V. ACKNOWLEDGMENTS
We are grateful to M.A.Zubkov for useful discussions and critical remarks.

Appendix A THE EFFECT OF GLUON BACKGROUND ROTATION ON THE AXIAL CURRENT
In this section we investigate if there are corrections to the CVE conductivity due to the rotation of the background gluonic fields B µ .We explicitly calculate their contribution at linear order in angular velocity Ω.
The rotation of the gluon background is incorporated by a change of the weight factor entering the ensemble average over is the gauge field (Yang-Mills) action of gluons, and S rot B is a modification that includes the effects of rotation.
The Euclidean metric in the rotating frame with an-gular velocity Ω around the x 3 -axis reads [34] The rotation amounts to imaginary contributions to the metric, as the change to imaginary time t → iτ implies Ω → −iΩ.To determine the contribution due to rotation, we change the sign of Ω and plug g rot µν into the action S B instead of g µν .Up to linear order in Ω we obtain for the gauge field action of gluons for rotation around the x 3 -axis where T µν denotes the gluon (Yang-Mills) energy momentum tensor.The temporal component of the energy momentum tensor also contributes an imaginary unit T 4i = 2i tr c F 4ν F ν i .Below we calculate the linear response of the quark axial current in Ω at the rotation axis δ ⟨j 5 µ (r = 0)⟩ in direction of the rotation axis to the energy-momentum flux of the background gluonic vacuum content.We use the notation of [48] for the field correlators form factors.
It is not only known how the vacuum background averaging of a single Wilson loop works but also for a pair of Wilson loops [49].We aim at bringing the linear response term into the latter form by replacing all non-Abelian field strength tensors by surface element variations (which are blind to the color trace structure and makes it necessary to consider a pair of Wilson loops).Due to the appearance of two traces in color space we introduce a Wilson loop over an infinitesimal contour C ′ within the trace comprising the Yang-Mills energy momentum tensor with the point limit of the contour implied throughout.
Let us denote the contribution to the axial current, introduce the auxiliary infinitesimal contour, apply the non-Abelian Stokes theorem, and replace all the field strength tensors by surface variations acting on the Wilson loop δ⟨j 5  3 (0)⟩ ∼⟨iΩ tr c,D (γ 5 γ 3 γ 4 exp(g =⟨2iΩ tr c,D (γ 5 γ 3 γ 4 exp(g The symbols P and P F as well as P ′ and P ′ F denote path and surface ordering, respectively.S min and S ′ min are the minimal areas, which perimeters are spanned by the closed quark path in the former case and by the auxiliary contour C ′ for the infinitesimal surface.x 0 is a reference point such that x 0 ∈ S min .Φ(x, y) is the parallel transporter associated with the background gluonic field from a point x to a point y.The field strength tensors supplemented with reference points are implicitly endowed with parallel transporters which we do not need to account for in the variational replacement since the transporters are equal to unity in the contour gauge.Nevertheless, we keep them for clarity.
The replacement of the field strength tensors comprising the energy momentum tensor by surface variations is only possible at points which lie on the mimimal surface S ′ min .We therefore require that x, x ′ ∈ C ′ .Since correlations between field strength tensors only persist up to distances of about one correlation length λ, correlation of the rotating white gluon contribution with the remaining terms implies that in the integral d 4 x only a volume comprising S min and its tubular neighborhood of thickness λ remains.The integration over d 4 x minus this tubular neighborhood of S min comprises factors ⟨T 4i (x)⟩ B = 0 with i = 1, 2, which disconnect from the remaining terms with respect to averaging over the gluonic background and vanish consequently (the gluonic background is isotropic).We will retain the full x-integration as long as possible and only later on keep the directions normal to S min by inclusion of the factor λ 2 and disregard corrections in the direction of continuation of the minimal surface (perimeter terms which are additionally suppressed).In order for the Dirac trace to be nonzero, the fermion spin interaction with the gluon background needs to be considered necessarily (contrary to the calculation of the CVE in the main text).
Then, as is customary within the FCM, we employ the cluster expansion, and retain only the Gaussian terms (quadratic in field strength tensors).We follow the notation of [49] for vacuum averaging of a pair of Wilson loops and obtain ))× (31) The omission of the reference point x 0 is exact up to corrections involving the parameter (λ/R) 2 , where λ ∼ 0.2f m ≈ 1 GeV −1 is the QCD vacuum correlation length, and R is some average size of an "overall white" system.We will work in this approximation which is commonly employed within the Field Correlator Method.The reference point is usually chosen such that the length over which parallel transporters act (and thereby the effect of choice of the reference point) is minimized.At this point we split the calculation which comprises a connected as well as a disconnected term After the differentiation we calculate the limit C ′ → 0, thus Λ 1 , Λ e → 0, and Λ 0 only comprises the term with the quark loop.
We start with the connected contribution, which is linear in the form factor D EH 1 .It embodies the colorelectric and color-magnetic cross correlation.To express the quadratic cumulant in terms of the form factor D EH 1 , which is defined as we keep the artificial point splitting z ∂z 2 e Λ 0 → 0.
We considered the limit x ′ → x by contracting C ′ to a point at the very end and made use of the freedom of contraction of C ′ to a point (first contract temporally and then spatially).Note that the even the singular terms in D EH 1 vanish, regardless of the limits order.
We now turn to the disconnected term which is quadratic in form factors.We neglect those terms containing the form factor D EH 1 , which contribution is small as compared to the main non-perturbative input provided by the color-electric vector-like form factor D E 1 (in form of the Polyakov loop) as well as the color-magnetic scalar form factor D H (in the form of color-magnetic confinement).This policy has been employed in the main text in (10), where it allows to split the Wilson loop into a product of the Polyakov line and a purely spatial loop subject to color-magnetic confinement.We still expect the following estimates to apply in general (that means taking also D EH We contracted the auxiliary Wilson loop to a point and thereby set x ′ = x.Further notice that S min comprises both a time-like as well as a space-like component.The first correlator in each term is determined by D E 1 , while the second one gives rise to D H . Thus, the u-integration extends in temporal direction while the v-integration is completely spatial.To sum up, a contribution from gluon rotation to the main term, which is ∼ T −1 ), of relative order O(λT (σ H,f λ 2 ) ).With this estimate, we conclude that the rotating gluon vacuum contribution is non-perturbatively suppressed and doesn't exceed 10% in comparison to the leading order contribution for T < 400 MeV.

2 Smin d 2 S 1 T 2 .
Both u-and v-integrations only contribute in the region |u − x|, |v − x| ≲ λ.As was already mentioned, this amounts to the replacement d 4 x → λ up to perimeter terms.The u-and x-integrations over D E 1 yield the contribution due to Polyakov lines, which may be shown to be of order O(α s λσ H,f ).The short-range correlation |u − x|, |v − x| ≲ λ implies that the winding around the Polyakov line is additionally suppressed by λT < 1.The remaining v-integration for a fixed x ∈ S min amounts at most to a contribution of order of σ H,f .The minimal surface in spatial direction is limited in size by |S min | ≲ 1 σ H,f due to the color magnetic confinement.Since we calculate the current at the rotation axis, the explicit x 1 and x 2 fulfill |x 1 |, |x 2 | ≲ 1 (σ H,f ) 1

3 2 )
. All explicit numerical factors combined are of order of unity.The omitted perimeter terms are additionally suppressed by O((σ H,f λ 2 ) 1 2