Ectp201517
Wlcapp201517
QCD Phasetransition and chemical freezeout in nonzero magnetic field at NICA
Abstract
Because of relativistic offcenter motion of the charged spectators and the local momentumimbalance experienced by the participants, a huge magnetic field is likely generated in highenergy collisions. The influence of such shortlived magnetic field on the QCD phasetransition(s) shall be analysed. From Polyakov linearsigma model, we study the chiral phasetransition and the magnetic response and susceptibility in dependence on temperature, density and magnetic field strength. The systematic measurements of the phasetransition characterizing signals, such as the fluctuations, the dynamical correlations and the inmedium modifications of rhomeson, for instance, in different interacting systems and collision centralities are conjectured to reveal an almost complete description for the QCD phasestructure and the chemical freezeout. We limit the discussion to NICA energies.
Keywords: Chiral transition, Magnetic catalysis, Critical temperature,
Viscous properties of QGP
PACS Nos: 11.10.Wx,25.75.Nq,98.62.En,12.38.Cy
1 Introduction
In most of the experimental signals revealing the deconfinement phasetransition such as the correlations and the elliptic flow of produced particles, the increase in the collision size and centrality is conjectured to improve their detection. The heavyion collisions  on the other hand  introduce a remarkable effect; a nonzero magnetic field, which can be generated due to the relativistic offcenter motion of the spectators electric charges and due to the local imbalance in the momentum carried by colliding nucleons, which leads to local angularmomentum [1, 2]. Accordingly, the life span of the resulting magnetic field can approximatively be estimated from the duration of both mechanisms. The relativistic offcenter motion of the electric charges and the local momentum imbalance of colliding nuclei keep up with the strong interaction causing them, i.e., sec. Thus, the resulting magnetic field is shortlived but its strength can be very huge, [3]. We focus the discussion to NICA energies, where Gauss, Fig. 3.
Quantum electrodynamical (QED) phenomena such as magnetic catalysis which describes how the magnetic field generates masses, dynamically and the Meissner effect in which the magnetic field changes the order of the phase transition in typeI superconductor are borrowed from solidstate physics in order to study the possible influence of the generated magnetic field on the quantum choromodynamical (QCD) phasespace structure [4] and/or the response of the hadronic and partonic matter to nonzero magnetic field in thermal and dense medium.
We have utilized the Polyakov linearsigma model (PLSM) in characterizing various QCD phenomena. Thermodynamic quantities including the higherorder moments of particle multiplicities are studied and compared with lattice QCD results [5, 6]. The chiral phasestructure of (pseudo)scalar and (axial)vector meson masses in thermal and dense medium was analysed, as well [7]. Furthermore, we have reported on some properties of the quarkgluon plasma (QGP) in nonzero magnetic field [8]. From these studies, various parameters of PLSM can be fixed. Based on this PLSM can also be utilized at high density, where the nonperturbative lattice technique becomes no longer applicable.
In the present work, we shall describe how the chiral phasetransition and the magnetic response and susceptibility are affected at high density and nonzero magnetic field. In doing this, we utilize PLSM in the mean field approximation. Upon availability, our calculations shall be compared with lattice QCD simulations.
2 A short reminder to Polyakov linearsigma model in finite magnetic filed
By assuming that the magnetic field () is directed along axis of the collision plane and from magnetic catalysis and Landau quantization [9], the magnetic field strength causes a considerable modification to the quark (and equivalently in antiquark) dispersion relations,
(1) 
where subscript runs over all quarks (and antiquarks) flavors, is the component of the spin in the magnetic field direction, being index label for Landau levels and is the electric charge of th quark flavor. From the magnetic catalysis property, a dimensional reduction can be implemented [9],
(2) 
where stands for degenerate Landau levels and is replaced by the quantum number .
Further details about the PLSM formalism can be taken from Ref. [5, 6, 7]. In the mean field approximation, the PLSM grand potential can be summarized as [6],
(3) 
The mesonic potential is assumed to vanish at high temperatures . The purely mesonic potential is given as [10].
(4) 
Due to the Landau quantization and the magnetic catalysis, the quarks and antiquark contributions to the PLSM potential are given as [8]
(5)  
where and . is given in Eq. (1) and refers to the chemical potential of th quark flavor. In vanishing magnetic field, the quarks potential reads
(6)  
Landau quantization, which is applied in order to add restrictions to the quarks due to the existence of free charges in the plasma phase, affects the strongly interacting matter by reducing the electromagnetic interactions between the quarks. In nonzero magnetic field, one still has to sum over color and flavor degrees of freedom with a dispersion relation given in Eq. (1), because the electric charges of the quarks are not equal, Eq. (5). On the other hand, the zero modes, Eq. (6), are taken into consideration in the present calculations. Furthermore, we have estimated the occupation of Landau levels at different temperatures, densities (chemical potentials), and magnetic field strengths (not shown here). We have concluded that the population of Landau levels is most sensitive to the magnetic field and the quark charges, as well. Having PLSM potential, Eq. (3), we can then straightforwardly deduce the thermodynamics quantities of interest.
3 Results and discussion
3.1 Chiral phasetransition at high density and nonzero magnetic field
In Fig. 1, the normalized chiral condensates ( and ) are depicted in dependence on the temperature at different baryon chemical potentials; (solid), (dashed), MeV (dotted curves), and various magnetic field strengths; (lefthand panel) and (righthand panel). In these calculations, the PLSM parameters are determined at (vacuum) sigmameson mass MeV, where the vacuum light and strange chiral condensates are measured as MeV and MeV with and , respectively. We notice that the values of the chiral condensates decline (having smaller critical temperatures) with increasing the baryon chemical potential. This means that, the chiral critical temperature decreases with the increase in the magnetic field strength (). It is obvious that the increasing in the baryon chemical potential () has the same effect. In other words, we conclude that the influences of the magnetic field are almost the same as that of the baryon chemical potentials, especially on decreasing the phasetransition temperature. The influence of the corresponding magnetic field (Gauss) becomes stronger with increasing density, e.g. NICA, which operates at the highest baryon density.
Fig. 2 presents and as functions of the baryon chemical potential at fixed temperatures; MeV (lefthand panel) and MeV (righthand panel). Both chiral condensates are calculated in different magnetic field strengths; (solid), (dashed), (dotted), (dotdashed), and (double dotted curves). Although these strengths are much larger than the one at NICA top energy (), they reveal the tendency of with the collision energies and the significant influence of the magnetic field strength can determined.
We conclude that increasing temperature causes a rapid decrease in the chiral condensates, especially near the chiral phasetransition, Fig. 1. This is similar to the observation reported in a previous study without magnetic field [7]. At high density corresponding to NICA, there is a gap difference between light and strange chiral condensates, Fig. 2. This can be understood because of the anomaly term in the mesonic part of the Lagrangian [7]. The anomaly term was proposed [10] to enhance the numerical estimation for the chiral condensates. Originally, this proposed term models the axial U() anomaly of the QCD vacuum, which is likely broken by quantum effects. As the anomaly term explicitly appears in the strange quark condensate, we conclude that it considerably counts for its difference relative to the light quark condensate. At high temperature and density, it remains finite for the strange quark condensate while entirely vanishes for the light quark condensate. For the sake of completeness, we recall here that this anomaly term also appears in the meson masses. Its effects have been systematically analysed for sixteen meson states [7]. Also, here its effects at high density are remarkably large. From Fig. 3, it is obvious that the magnetic field causes a considerable effect on the chiral phasestructure with varying temperature. Reducing the magnetic fields leads to an increase in the chiral criticaltemperature, righthand panel of Fig 3. Furthermore, raising the magnetic field strength sharpens the QCD phasestructure and accelerates the formation of the metastable phase. The latter is to be understood due to early phasetransition, i.g. smaller critical temperature.
Now, we analyse the influence of the magnetic field strength on the chiral phasediagram. In Fig. 3, the chiral temperature (lefthand panel) and baryon chemical potential (righthand panel) are calculated in dependence on . is calculated according two criteria. The first one implements the secondorder moment of the quark multiplicity, in which is determined at the peak of the normalized quark susceptibility (). This criterion is frequently utilized in lattice QCD simulations. Depending on the order of the phase transition, the corresponding critical temperature can be determined. In crossover, which is likely for light and strange quarks at low density, the phasetransition takes place in a range of temperatures. Thus, an average value can be assigned to the quasicritical temperature. Our results are given by the solid curves. The second criterion defines as the intersect of the lightquark chiral condensate and the deconfinement orderparameter [6]. As illustrated in Figs. 1 and 2, the chiral quasicritical temperature can be determined from and , separately, which should become coincident with the deconfinement phasetransition. In PLSM, the latter is related to the temperature dependence of the order parameters and . Both condensates seem to have crossover phase transitions, where their values decrease from large values at low temperatures to low values at high temperatures. Accordingly, can be approximated as an average within the crossover region. In order to increase the certainty of the estimated , we utilize coincidence of both chiral condensates with and . The corresponding results are depicted by dashed curves. The vertical bands mark at GeV comparable to NICA (), RHIC beam energy scan ranging from to GeV () [11, 12] and GeV at LHC () calculated from PYTHIA with centrality and impact parameter fm.
We observe that decreases with increasing , i.e., inverse magnetic catalysis. A good agreement between dotted curve (second criterion) and the lattice results, especially at GeV, is obtained. At a wider range of ; GeV, the solid curve (first criterion) matches well with the lattice calculations [13]. estimated from the first criterion apparently overestimates the lattice simulations, especially at small . The PLSM calculations refer to larger than the one deduced from the lattice QCD. On the contrary, the second criterion seems to underestimates (relative to the lattice) at large .
In righthand panel (b), in a constant magnetic field such as the one marked by the vertical bands, we observe that, a larger can be reached at a lower temperature, i.e., decreases with increasing . In determining the critical baryon chemical potential at which the broken chiral symmetry should be restored, we implement the intersection between deconfinement phasetransition and the lightquark chiral condensate. At fixed temperatures, (solid curve) and MeV (dashed curve), the dependence of critical on is illustrated in the righthand panel of Fig. 3. In determining the critical chemical potential, a second criterion, the intersection between deconfinement orderparameter and lightquark chiral condensate, can be implemented.
3.2 Magnetic response and susceptibility
As mentioned in previous sections, the nonzero magnetic field comes up with modifications to the dispersion relation so that the response of the QCD matter can be determined from the derivative of the free energy density with respect to the magnetic field itself
(7) 
where is the elementary electric charge. The sign of the magnetization defines whether the QCD matter is para or diamagnetic, i.e. (bara), or (dia). Similar to solidstate physics,

if the QCD matter is diamagnetic, then the colorcharges are oriented oppositely to the magnetic field direction and an induced current shall be produced. The induced current spreads in form of small loops which  in tern  try to cancel out the magnetic effects, and

if the QCD matter is paramagnetic, then most colorcharges shall be aligned towards the magnetic field direction.
Recently, the magnetic susceptibility of the QCD matter has been reported [14]
(8) 
The magnetization is a dimensionless proportionality parameter indicating the degree of magnetization and depending on the magnetic permeability [15]. Thus, it is apparent that the temperature weakens the magnetization and at relativistic energies (or vanishing chemical potential), the magnetization likely vanishes.
In lefthand panel of Fig. 4 (a), the magnetization of the QCD matter due to the effects of nonvanishing magnetic field GeV is studied as a function of temperature at a vanishing baryon chemical potential and compared with recent lattice calculations (open triangles with errorbars) [13]. Increasing the magnetization with the temperature obviously refers to positive magnetization, , i.e., the paramagnetic property of the QCD matter becomes dominant. At temperatures smaller than the critical one, the calculations agree well with the lattice results. At temperatures larger than the critical one, we find that our PLSM calculations slightly overestimate the lattice results. In this limit, the degrees of freedom might not be sufficient to achieve a good agreement at very high temperature. Furthermore, such a discrepancy can be interpreted due to the PLSM applicability, which is determined up to the temperature range, where , , and remain finite. For the sake of completeness, we emphasize that the lattice results at GeV are obtained by using halfhalf method for three lattice temporal extensions; and and the continuum limit.
The righthand panel of Fig. 4 (b) presents the magnetic susceptibility as a function of temperature at GeV and a vanishing chemical potential. The results from PLSM are compared with various lattice simulations (symbols) in which different simulation methods and lattice configurations are applied. Features of PLSM and lattice QCD results can be summarized as follows.

The PLSM free energy, Eq. (3), consists of three parts. The first one is the gauge potential (pure meson). The second one represents the quark and antiquark potential, which includes fluctuations of quarks and antiquark. The third part takes into consideration the color and gluon interactions. This means that two types of contributions are not free of fluctuations, while the third one takes into consideration the gluon interactions.

At low temperatures, the slope of the magnetic susceptibility becomes negative (insidebox in the righthand panel). This signals a diamagnetic property and apparently confirms the lattice results, as well. Such result has been also obtained from partonhadronstring dynamics [16]. On the other hand, at high temperature, i.e., when the broken chiral symmetry is restored, a transition between dia and paramagnetic properties is observed. Similarly, from noninteracting MIT bag model, the phasetransition from dia to paramagnetism was observed [17]. QCDlike models with Polyakov potential agree well with the lattice calculations, especially at high temperatures [18].

The lattice simulations [13] implement halfhalf method in lattice (closed triangle) and integral method in lattice (open triangle). Thus, the continuum limit is likely achieved [19]. For degrees of freedom and by using HISQ/tree action with quark masses and , the lattice results are represented by diamonds [20]. The closed circles stand for isotropy lattice [21].
4 Conclusions
Although the magnetic field strength in heavyion collisions at NICA energies is small relative to top RHIC and LHC energies, Fig. 3, we observe that even such relative small considerably decreases the critical temperatures with increasing magnetic field strength, i.e., inverse magnetic catalysis. At finite magnetic field, we find that the phasetransition seems to take place at higher critical temperatures than that of the restoration of the chiral symmetry breaking. The estimation of both types of critical temperatures is effectively achieved from the location (in temperature axis) of the secondorder moment of the quark multiplicity, know as susceptibility and from the intersection of light and strange quark condensates with the deconfinement order parameters and .
Both chiral condensates ( and ), which can be characterized through thermal, dense and magnetic inmedium modifications of rhomeson, for instance [24]. We have presented results on their evolution with increasing temperature, chemical potential and magnetic field. We find that the chiral condensates rapidly decrease with increasing the baryon chemical potential, i.e., the chiral critical temperature decreases with the increase in the magnetic field strength. We conclude that the influences of the magnetic field are almost the same as that of the baryon chemical potentials, especially on decreasing the phasetransition temperature. Thus, the influence of the magnetic field at NICA energies () becomes stronger because of the high density. The susceptibilities and fluctuations of different particle yields are examples about possible experimental signatures for the deconfiment phasetransition and critical endpoint. The latter is conjectured to connect the crossover of the QCD boundary at low density (chemical potential) and the firstorder phase transition at high density [25]. In light of this, we recommend systematic measurements for the phasetransition characterizing signals such as the fluctuations and the dynamical correlations in various interacting systems and centralities. These shall reveal a complete description for the QCD phasestructure and for the chemical freezeout at high density.
The magnetization and magnetic susceptibility almost monotonically increase with increasing temperature. At low temperatures, the magnetic susceptibility has negative slope, i.e., the QCD matter possesses diamagnetic property. At high temperature, i.e., when the broken chiral symmetry is restored, the QCD matter becomes paramagnetic. In both temperature regimes, the PLSM calculations and the recent lattice calculations are in good agreement. In order to justify this agreement, we have shortly reviewed the configurations and the characteristics of the lattice QCD calculations.
Furthermore, we conclude that increasing critical temperature (as a result of decreasing magnetic field) is accompanied by a considerable decrease in the baryon chemical potential. We have analysed how at a fixed temperature the increase in decreases the chiral critical temperature, which  in turn  decreases further with increasing the temperature. Nearly the same effects has the baryon chemical potential.
References
 [1] ZT Liang and XN Wang Phys. Rev. Lett. 94 102301 (2005)
 [2] F Becattini, F Piccinini and J Rizzo Phys. Rev. C77 024906 (2008)
 [3] D E Kharzeev, L D McLerran and H J Warringa Nucl. Phys. A803 227 (2008)
 [4] E S Fraga and A J Mizher Phys. Rev. D78, 025016 (2008)
 [5] A Tawfik and N Magdy J. Phys. G42 015004 (2015)
 [6] A Tawfik, N Magdy and A Diab Phys. Rev. C89 055210 (2014)
 [7] Abdel Nasser Tawfik and Abdel Magied Diab Phys. Rev. C91 015204 (2015)
 [8] Abdel Nasser Tawfik and Niseem Magdy Phys. Rev. C90 015204 (2014)
 [9] I A Shovkovy Lect. Notes Phys. 871 13 (2013)
 [10] BJ Schaefer and M Wagner Phys. Rev. D79 014018 (2009)
 [11] D E Kharzeev, J Liao, S A Voloshin, and G Wang Prog. Part. Nucl. Phys. 88 1 (2016)
 [12] K Fukushima, D E Kharzeev and H J Warringa Phys. Rev. Lett. 104 212001 (2010)
 [13] G Bali, F Bruckmann, G Endrodi, S Katz and A Schaefer J. High Energy Phys. 1408 177 (2014)
 [14] K Kamikado and T Kanazawa J. High Energy Phys. 1501 129 (2015).
 [15] J R Hook and H E Hall Solid state physics (Chichester : Wiley) Vol 1 (1995)
 [16] T Steinert and W Cassing Phys. Rev. C89 035203 (2014)
 [17] N Agasian and S Fedorov Phys. Lett. B663 445 (2008)
 [18] V Orlovsky and Y A Simonov Int. J. Mod. Phys. A30 1550060 (2015)
 [19] C Bonati, M D’Elia, M Mariti, F Negro and F Sanfilippo Phys. Rev. D89 054506 (2014)
 [20] L Levkova and C DeTar Phys. Rev. Lett. 112 012002 (2014)
 [21] G Bali, F Bruckmann, G Endrodi and A Schafer Phys. Rev. Lett. 112 042301 (2014)
 [22] Y Aoki, S Borsanyi, S Durr, Z Fodor, S D Katz, S Krieg and K K Szabo J. High Energy Phys. 0906 088 (2009)
 [23] S Borsanyi, G Endrodi, Z Fodor, A Jakovac, S D Katz, S Krieg, C Ratti, and K K Szabo J. High Energy Phys. 1011 077 (2010)
 [24] A Tawfik Indian J. Phys. 85 755 (2011)
 [25] U S Gupta, R K Mohapatra, A M Srivastava and V K Tiwari Indian J. Phys. 85 115 (2011)