Nonperturbative QCD thermodynamics in the external magnetic field

The thermodynamics of quarks and gluons strongly depends on the vacuum colormagnetic field, which grows with the temperature T, so that spatial string tension σs=constg4(T)T2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _s =\mathrm{const}~ g^4 (T) T^2$$\end{document}. We investigate below what happens when one imposes in addition constant magnetic field and discover remarkable structure of the resulting thermodynamic potential.


Introduction
The problem of quark gluon thermodynamics in magnetic field is of the high interest in the modern physics, since heavyion experiments produce important data on the properties of resulting hadron yields and hadron interactions, which might be influenced by the strong magnetic fields (MF) created during the collision process [1][2][3]. For a recent review on the effects of MF see [4].
On the theoretical side the problem of MF in the quark gluon plasma (qgp) was studied in different aspects, e.g. in the NJL-type models [5] and in the holographic approach [6,7]. Within the nonperturbative QCD the theory of qgp in MF was developed in [8,9,11], where the general form of the thermodynamic potentials was found in MF with zero or nonzero baryon density, summing over all Landau levels including LLL.
In this approach the only nonperturbative interaction, which was taken into account, reduced to the inclusion of Polyakov lines in the resulting expression for the pressure.
The resulting expressions for magnetic susceptibilitieŝ χ q (T ), obtained in [8], were used in [9,11] to compare with the lattice data from [12,13], and a reasonable agreement was found forχ q (T )with different q = u, d, s, in [9,11] as well as for the sum [11], however somewhat renormalised values of effective quark masses were used. a e-mail: andreichicov@mail.ru Recently in [14,15] a new step in the development of the np QCD thermodynamics was made, where the colormagnetic confinement (CC) was included in the dynamics of the qgp. This interaction with the spatial string tension σ s grows with temperature,σ s ∼ g 4 T 2 and is important in the whole region T c < T < 10 GeV. A concise form of the final expression was found in [14,15] in the case of an oscillatory type CC and an approximate one in the realistic case of the linear CC. The resulting behavior both in the SU(3) case, found in [16,17], and in the qgp case [14,15] agrees well with the corresponding lattice data.
It is the purpose of the present paper to extend our previous analysis of the qgp thermodynamics in MF, done in [8,9,11], including the dynamics of CC with the explicit form of the magnetic screening mass m D , generated by CC.
As will bee seen, we propose a simple generalization of the results [8,9,11], where the CC produces the mass M D , entering the final expressions in the combination m 2 q + m 2 D ≡ M instead of m q . We check the limiting cases and compare the result with lattice data.
In the next section the general analysis of the MF effects in thermodynamics is explained, in Sect. 3 the magnetic susceptibility is defined, in Sect. 4 the results are compared to the lattice data.

General structure of the pressure with and without magnetic field
Field correlator method (FCM) [18,19], based on the path integral formalism and Feynman-Fock-Schwinger proper time formalism (FSF) [20,22], provides a relativistic-and gauge-invariant approach to QCD. Both perturbative and nonperturbative (np) gauge interactions between quarks could be written in a form of expansion in Wilson loops. A detailed review of the FCM technique and background field theory applied to the quark-gluon thermodynamics could be found in [21,23]. Here we emphasize that two basic interactions -colorelectric (CE) and colormagnetic (CM) are defined by two pairs of CE (D E , D E 1 ) and CM (D H , D H 1 ) gauge-invariant correlators obtained through the averaging over the stochastic background color field in FCM where u = |x − y|, the parallel transporter is = e ig y x A μ dx μ , and its contour is parametrized with a straight line.
In the confinement region, T < T c , spatial CE string tension is defined by scalar Since the D E vanishes beyond the deconfinement transition region, the CE interaction is provided only by the D E 1 corre- , via the np saturating potential V 1 (r, T ), and perturbative part yields color Coulomb potential. Despite the vanishing of the D E , the scalar part of the CM interaction D H still holds in the deconfinement region, yielding the spatial string tension σ S (T ) = 1 2 D H (z)d 2 z. As it was shown in [16,24], the magnetic string tension grows with temperature as σ S (T ) = cg 4 (T )T 2 , c = 0.566 [16] and therefore is important in the whole region T > T c .
Using T -dependent path integral formalism, one can express quark and gluon pressures via the Wilson loop series where s is Schwinger proper time amd G (n) (s) are the winding (Matsubara) path integrals dτ is a kinetic kernel in the Schwingr proper time method, m q is a current quark mass, and W (C n ) is the Wilson loop defined by the quark path C n . It is convenient to make a redefinition of the kinetic kernel splitting it into spatial and temporal components K = K 4 + K 3d =  [21] CE and CM fields in T > T c region correlate very weakly, since the field correlator E i (x) H k 0 and therefore temporal and spatial Wilson loops can be factorized where, after the applying FCM background averaging porcedure, the temporal projection L (n) defines a Polyakov loop and the dynamics of the spatial projection W 3 is defined by the CM string tension σ s . As a result, one can integrate out the z 4 part of the path integral in (4) and write the result as Here S 3 (s) is a 3d proper time Green's function with a closed path, which can be expanded in a series of eigenstates of the quark-antiquark pair in FSF formalism with the CM confinement interaction between quarks [15].
where the ground state energy for the 3d spatial motion is defined by a color-magnetic Debye mass m 2 D . As it was argued in [15,16], in the case of linear potential for the CMC and in the absence of the external MF, one obtains for S 3 (s) an approximate form and the pressure can be written as where One can see that the Debye mass m D growing with T provides an effective damping of the quark pressure. Due to this effect the high-T limit of P q never tends to Stefan-Boltzmann value. The same behaviour was observed in lattice calculations [10] where the difference of P T 4 with Stefan-Boltzmann limit is ∼ 25% at T ≥ 500 MeV. In the presence of the external magnetic field (MF) B||z, the discrete spectrum of the spatial Green's function (7) modifies [8] as and the summation over eigenstates ν is replaced by the summation over Landau levels and spin projections ν → n ⊥ ,σ . An influence of the MF in |eB| ≥ m 2 D 0.3 GeV 2 (Fig. 1) region dominates over the CMC interaction for quarks. It os sufficient to assume that the quark spectrum comprises Landau levels and self-energy quark mass m 2 D . However, the detailed behaviour in 0 < |eB| < m 2 D ∼ g 2 (T )T 2 region in B − T plane requires to solve the spectral problem in (7) with both CM and MF interactions. One should also take into account an additional multiplier proportional to eB and s due to the |ψ ν | 2 normalization in the presence of the MF during the phase space modification V 3 As a result, one can rewrite (9) as The total quark pressure is defined as a sum over flavors P q = f P f q . We have neglected above the influence of MF on the gluonic component of pressure, since in the deconfined region it is of pertubative origin only. As was shown in [8], the form (12) can be summed up over N ⊥ and σ to obtain the following result Note, that the first term in (13) is provided by the lowest Landau levels (LLL). For these levels it is known from analysis in [26], that the asymptotic quark energy values do not depend on eB and are equal to √ σ for small quark mass.
This agrees with our valuesM = m 2 q + σ s ≈ √ 4σ s and supports our expression (13) at least in the high e q B limit, e q B M , when the second and the third term in (13) tend to zero. Hence one obtains in the limit |e q B| M , T Note, that the factor |e q B| appears due to the phase space factor in MF, discussed above.. Now we turn to the limit of small MF, |e q B| M , T .
One obtains from (13) the contribution of the second term only One can compare (14) with (9), obtained in the case of zero MF, and insertion of (10) in (9) yields the same answer as in (15).

Magnetic susceptibility of the quark matter
Using general expression for the quark pressure (13), one can define a more convenient quantity, the magnetic susceptibilitiesχ To this end one expands the Mc Donald functions K n ( √ n 2 + b 2 ) entering in (13) in powers of b, following [8], and one obtains As a consequence, one has forχ q It is possible to sum up the series over n in (17), when one exploits the representation As a result one obtains for the quadratic magnetic susceptibility (ms) The ms in (21) is defined for a given quark flavor q, and the total m.s. for the quark ensemble, e.g. for 2 + 1 species of quarks can be written as

Results and discussions
We have presented above the thermodynamic theory of quarks in the magnetic field, when quarks are affected also by Polyakov line interaction and the CMC interaction, generalizing in this way our old results of [8,9,11], where the CMC part was absent.
We have included the CMC interaction in the energy eigenvalues ε σ n ⊥ , Eq. (7) tentatively via the replacement (11). This substitute can be corroborated in the case of lowest Landau levels withσ = 1, n ⊥ = 0, where magnetic field does not eneter, and the effective quark mass in subject to the CMC interaction only. In the general case one can expect possible interference of eB and CC terms, which can spoil the suggested replacement.
To make this first analysis more realistic, we have checked the limits of small and large values of eB. In the first case we have shown the correct correspondence with the eB = 0 result of [14,15], and in the second case of large eB, the leading linear in eB term is just LLL term, which is not influenced by magnetic fields, except for the phase space redefinition. These results enable us to proceed with the analysis and comparisons of obtained equations with lattice data.
We present below an analysis of the MF influence on the quark thermodynamics. It is interesting, that the basic expression for the pressure P ( f ) q in (14) has the property, that at small eB < (eB) crit , the dependence of P(B, T ) on eB is quadratic with a good accuracy, according to Eq. (17). At larger eB, eB > (eB) crit , one has the linear dependence of P(B, T ) on eB, given in (15). Figure 1 illustrates this behaviour for thee fixed temperatures T = 0.15, 0.2, 0.4 GeV. One can see from (14) and Fig. 1 that (eB) crit >M, and actually is around 0.5 GeV 2 for T 0.2 GeV. We have computed analytically the difference P(T ) = P(T, eB) − P(T, 0) using Eq. (14) and compare with the lattice data from [27] for averaged u, d, s quark ensemble. One can see in Fig. 2 the normalized pressure P(B, T ) for T > 0.135 GeV and eB = 0.2, 0.4 GeV 2 for the averaged quark ensemble of u,d,s quarks. At this point one should notice, that the calculations in FCM formalism were performed for the quark-gluon plasma in the infinite volume, which yields to the isotropic (thermodynamic) pressure distribution in the presence of the external MF. The effect of the pressure anisotropy observed on the lattice [27] is provided by the finite volume effects caused by the surface currents on the volume boundary, which are induced by the bulk magnetization of the quark-gluon plasma, as was discussed in [29,30]. Neglecting these effects we are comparing  [27]. For lattice results, the line thickness corresponds to the estimated error of the calculation our results with the lattice scheme-independent longitudial pressure from [27]. In Fig. 2 are shown our calculations of P(B, T ) for eB = 0.2, 0.4 GeV 2 and m 2 D = 3σ s using (14) are compared with the lattice data from [27]. We have used in (14) the Polyakov line L(T ) obtained in [28]. One can see a reasonable agreement within the accuracy of the lattice data, which supports the main structure of the theory used in the paper. Detailed analysis for larger intervals of eB and T and for specific quark flavours is possible within the approach and is planned for next publications. and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .