QCD Phase-transition and chemical freezeout in nonzero magnetic field at NICA
Because of relativistic off-center motion of the charged spectators and the local momentum-imbalance experienced by the participants, a huge magnetic field is likely generated in high-energy collisions. The influence of such short-lived magnetic field on the QCD phase-transition(s) shall be analysed. From Polyakov linear-sigma model, we study the chiral phase-transition and the magnetic response and susceptibility in dependence on temperature, density and magnetic field strength. The systematic measurements of the phase-transition characterizing signals, such as the fluctuations, the dynamical correlations and the in-medium modifications of rho-meson, for instance, in different interacting systems and collision centralities are conjectured to reveal an almost complete description for the QCD phase-structure 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
In most of the experimental signals revealing the deconfinement phase-transition 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 heavy-ion collisions - on the other hand - introduce a remarkable effect; a nonzero magnetic field, which can be generated due to the relativistic off-center motion of the spectators electric charges and due to the local imbalance in the momentum carried by colliding nucleons, which leads to local angular-momentum [1, 2]. Accordingly, the life span of the resulting magnetic field can approximatively be estimated from the duration of both mechanisms. The relativistic off-center 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 short-lived but its strength can be very huge, . 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 type-I superconductor are borrowed from solid-state physics in order to study the possible influence of the generated magnetic field on the quantum choromodynamical (QCD) phase-space structure  and/or the response of the hadronic and partonic matter to nonzero magnetic field in thermal and dense medium.
We have utilized the Polyakov linear-sigma model (PLSM) in characterizing various QCD phenomena. Thermodynamic quantities including the higher-order moments of particle multiplicities are studied and compared with lattice QCD results [5, 6]. The chiral phase-structure of (pseudo)-scalar and (axial)-vector meson masses in thermal and dense medium was analysed, as well . Furthermore, we have reported on some properties of the quark-gluon plasma (QGP) in nonzero magnetic field . From these studies, various parameters of PLSM can be fixed. Based on this PLSM can also be utilized at high density, where the non-perturbative lattice technique becomes no longer applicable.
In the present work, we shall describe how the chiral phase-transition 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 linear-sigma 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 , the magnetic field strength causes a considerable modification to the quark (and equivalently in antiquark) dispersion relations,
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 ,
where stands for degenerate Landau levels and is replaced by the quantum number .
The mesonic potential is assumed to vanish at high temperatures . The purely mesonic potential is given as .
Due to the Landau quantization and the magnetic catalysis, the quarks and antiquark contributions to the PLSM potential are given as 
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
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 phase-transition 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; (left-hand panel) and (right-hand panel). In these calculations, the PLSM parameters are determined at (vacuum) sigma-meson 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 phase-transition 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 (left-hand panel) and MeV (right-hand panel). Both chiral condensates are calculated in different magnetic field strengths; (solid), (dashed), (dotted), (dot-dashed), 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 phase-transition, Fig. 1. This is similar to the observation reported in a previous study without magnetic field . 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 . The anomaly term was proposed  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 . 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 phase-structure with varying temperature. Reducing the magnetic fields leads to an increase in the chiral critical-temperature, right-hand panel of Fig 3. Furthermore, raising the magnetic field strength sharpens the QCD phase-structure and accelerates the formation of the metastable phase. The latter is to be understood due to early phase-transition, i.g. smaller critical temperature.
Now, we analyse the influence of the magnetic field strength on the chiral phase-diagram. In Fig. 3, the chiral temperature (left-hand panel) and baryon chemical potential (right-hand panel) are calculated in dependence on . is calculated according two criteria. The first one implements the second-order 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 cross-over, which is likely for light and strange quarks at low density, the phase-transition takes place in a range of temperatures. Thus, an average value can be assigned to the quasi-critical temperature. Our results are given by the solid curves. The second criterion defines as the intersect of the light-quark chiral condensate and the deconfinement order-parameter . As illustrated in Figs. 1 and 2, the chiral quasi-critical temperature can be determined from and , separately, which should become coincident with the deconfinement phase-transition. In PLSM, the latter is related to the temperature dependence of the order parameters and . Both condensates seem to have cross-over 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 cross-over 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 . 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 right-hand 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 phase-transition and the light-quark chiral condensate. At fixed temperatures, (solid curve) and MeV (dashed curve), the dependence of critical on is illustrated in the right-hand panel of Fig. 3. In determining the critical chemical potential, a second criterion, the intersection between deconfinement order-parameter and light-quark 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
where is the elementary electric charge. The sign of the magnetization defines whether the QCD matter is para- or dia-magnetic, i.e. (bara-), or (dia-). Similar to solid-state physics,
if the QCD matter is dia-magnetic, then the color-charges 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 para-magnetic, then most color-charges shall be aligned towards the magnetic field direction.
Recently, the magnetic susceptibility of the QCD matter has been reported 
The magnetization is a dimensionless proportionality parameter indicating the degree of magnetization and depending on the magnetic permeability . Thus, it is apparent that the temperature weakens the magnetization and at relativistic energies (or vanishing chemical potential), the magnetization likely vanishes.
In left-hand 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) . 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 half-half method for three lattice temporal extensions; and and the continuum limit.
The right-hand 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 (inside-box in the right-hand panel). This signals a dia-magnetic property and apparently confirms the lattice results, as well. Such result has been also obtained from parton-hadron-string dynamics . On the other hand, at high temperature, i.e., when the broken chiral symmetry is restored, a transition between dia- and para-magnetic properties is observed. Similarly, from non-interacting MIT bag model, the phase-transition from dia- to para-magnetism was observed . QCD-like models with Polyakov potential agree well with the lattice calculations, especially at high temperatures .
The lattice simulations  implement half-half method in lattice (closed triangle) and integral method in lattice (open triangle). Thus, the continuum limit is likely achieved . For degrees of freedom and by using HISQ/tree action with quark masses and , the lattice results are represented by diamonds . The closed circles stand for isotropy lattice .
Although the magnetic field strength in heavy-ion 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 phase-transition 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 second-order 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 in-medium modifications of rho-meson, for instance . 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 phase-transition 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 phase-transition and critical endpoint. The latter is conjectured to connect the cross-over of the QCD boundary at low density (chemical potential) and the first-order phase transition at high density . In light of this, we recommend systematic measurements for the phase-transition 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 phase-structure 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 dia-magnetic 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.
-  Z-T Liang and X-N Wang Phys. Rev. Lett. 94 102301 (2005)
-  F Becattini, F Piccinini and J Rizzo Phys. Rev. C77 024906 (2008)
-  D E Kharzeev, L D McLerran and H J Warringa Nucl. Phys. A803 227 (2008)
-  E S Fraga and A J Mizher Phys. Rev. D78, 025016 (2008)
-  A Tawfik and N Magdy J. Phys. G42 015004 (2015)
-  A Tawfik, N Magdy and A Diab Phys. Rev. C89 055210 (2014)
-  Abdel Nasser Tawfik and Abdel Magied Diab Phys. Rev. C91 015204 (2015)
-  Abdel Nasser Tawfik and Niseem Magdy Phys. Rev. C90 015204 (2014)
-  I A Shovkovy Lect. Notes Phys. 871 13 (2013)
-  B-J Schaefer and M Wagner Phys. Rev. D79 014018 (2009)
-  D E Kharzeev, J Liao, S A Voloshin, and G Wang Prog. Part. Nucl. Phys. 88 1 (2016)
-  K Fukushima, D E Kharzeev and H J Warringa Phys. Rev. Lett. 104 212001 (2010)
-  G Bali, F Bruckmann, G Endrodi, S Katz and A Schaefer J. High Energy Phys. 1408 177 (2014)
-  K Kamikado and T Kanazawa J. High Energy Phys. 1501 129 (2015).
-  J R Hook and H E Hall Solid state physics (Chichester : Wiley) Vol 1 (1995)
-  T Steinert and W Cassing Phys. Rev. C89 035203 (2014)
-  N Agasian and S Fedorov Phys. Lett. B663 445 (2008)
-  V Orlovsky and Y A Simonov Int. J. Mod. Phys. A30 1550060 (2015)
-  C Bonati, M D’Elia, M Mariti, F Negro and F Sanfilippo Phys. Rev. D89 054506 (2014)
-  L Levkova and C DeTar Phys. Rev. Lett. 112 012002 (2014)
-  G Bali, F Bruckmann, G Endrodi and A Schafer Phys. Rev. Lett. 112 042301 (2014)
-  Y Aoki, S Borsanyi, S Durr, Z Fodor, S D Katz, S Krieg and K K Szabo J. High Energy Phys. 0906 088 (2009)
-  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)
-  A Tawfik Indian J. Phys. 85 755 (2011)
-  U S Gupta, R K Mohapatra, A M Srivastava and V K Tiwari Indian J. Phys. 85 115 (2011)