Chiral magnetic properties of QCD phase-diagram

The QCD phase diagram is studied, at finite magnetic field. Our calculations are based on the QCD effective model, the SU($3$) Polyakov linear sigma model (PLSM), in which the chiral symmetry is integrated in the hadron phase and in the parton phase, the up-, down- and strange-quark degrees of freedom are incorporated besides the inclusion of Polyakov loop potentials in the pure gauge limit, which are motivated by various underlying QCD symmetries. The Landau quantization and the magnetic catalysis are implemented. The response of the QCD matter to an external magnetic field such as magnetization, magnetic susceptibility and permeability has been estimated. We conclude that the parton phase has higher values of magnetization, magnetic susceptibility, and permeability relative to the hadron phase. Depending on the contributions to the Landau levels, we conclude that the chiral magnetic field enhances the chiral quark condensates and hence the chiral QCD phase diagram, i.e. the hadron-parton phase transition likely takes place, at lower critical temperatures and chemical potentials.


Introduction
Exploring the quantum chromodynamic (QCD) phasediagram and studying the phase structures and the deconfinement phase-transitions of strongly interacting matter are among the fundamental issues in nuclear physics. Studying QCD matter in laboratory is one of the greatest challenges of the experimental nuclear physics as the parton matter is not directly accessible. There are different experimental methods implemented in accelerating and colliding ions (hadrons). Due to the center-of-mass energy of the collision achieved, different domains of the QCD phase-diagram could be drawn. The first prediction of end of the hadron domain, at high temperatures, was fomulated long time before the invention of QCD, where the partons are assumed as the effective degrees-of-freedom (dof), at temperatures larger than the Hagedorn temperature T H [1,2]. The hadron matter forms fireballs of new particles, which can again produce new fireballs. In 1975, Cabibbo proposed a QCD phase diagram in T − n B plane [3], where T H in the statistical Bootstrap model (SBM) [1] was interpreted as the critical temperature T c and is conjectured to be associated to second-order phase-transition into the deconfinement state. At large baryon density, n B , a weak interaction between quarks and gluons -due to the asymptotic dof -has been recognized [4]. A Historical suma Present address: Goethe University, Institute for Theoretical Physics (ITP), Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany mary on the QCD phase diagram and its investigation in the heavy-ion collisions (HIC) experiments are available, for instance, in refs. [5,6].
In statistical physics, the phase transitions are defined as singularities or non-analyticity in the free energy as a function of thermodynamic quantities. In lattice QCD simulations, the corresponding partition functions are taken as functional integrals over compact groups described and evaluated in dependence on the temperature T , the chemical potential µ, the volume V , the magnetic field eB, · · · , etc.
In HIC and due to oppositely ultra-relativistic motion of colliding heavy-ions, a huge magnetic field can be generated. Their motion generates an electric current which in turn induces magnetic field to the system. In a noncentral HIC, the two counter-propagating nuclei collide, at finite impact parameter b. In Fig. 1 the magnetic field in the center of the over-lapping surface in Au + Au collisions, for instance, at b = 10 fm and √ s N N = 200 GeV, is visualized. This is perpendicular to the reaction plane owing to the symmetry of the collision. Let us assume that all colliding nuclei are located, at the center of the nucleus, by applying Biot-Savart law, where the negative sign appears due to the assumption that the magnetic field is pointing in −y direction, v z = (1 − (2m N / √ s) 2 ) 1/2 ≈ 0.99 is the velocity of accelerated nuclei, m N is the nucleon mass, γ = 1/ 1 − v 2 z ≈ 100 is the Lorentz gamma, and Z Au = 79 is the change number of the gold nucleus. It is assumed that such a huge magnetic field would have great impacts on the dynamics of the parton (quarks and gluons) matter produced in HIC.
Over the last few decades, great efforts have been done to map out the QCD phase-diagram. Right panel of Fig.  1 illustrates a schematic version. At low temperatures and baryon chemical potentials, the quarks and gluons are still confinement forming colorless hadrons. At T ∼ 150 MeV, there is a crossover to the partonic colored phase; quarkgluon plasma. With increasing baryon-chemical potential, at low temperatures, the quarks shall be grouped in pairs known as correlated Cooper-pairs, which likely condense. Various experiments aim at studying the phase structures and the QCD phase-diagram such as RHIC at BNL √ s N N = 200 GeV and LHC at CERN with energies up √ s N N = 13 Tev. At √ s N N = 4 − 11 GeV, the large density regime shall be explored by CBM at FAIR and MPD at JINR.
In the heavy-ion experiments, as a result of the noncentral heavy-ion collisions, a huge magnetic field could be created. For instance, at RHIC and LHC energies, the magnetic field ranges between m 2 π and 10 − 15m 2 π , respectively [8,9], where m 2 π ∼ 10 8 Gauss. It is worth mentioning that this value of course is just a snapshot, since the magnetic field is strongly time-dependent. Detailed discussion on such a dynamical system does not lay within the scope of the present study. On the other hand, the proposed approach, the PLSM, take into account the evolution of such dynamical system. The magnetic field lifetime in HIC including the electric and chiral magnetic effects [10,11] and the electromagnetic impacts on the heavy-ion phenomenology are reviewed in ref. [12]. The features of the electromagnetic fields in HIC shall be addressed, quantitatively.
In 1960s, the linear sigma model (LSM), a low-energy model, was introduced by Gell-Mann and Levey [32], long before the invertion of the Quantum Chromdynamics (QCD), the theory of the strong interaction. Many studies have been performed on LSM O(4) at (non)-zero temperature [33,34] and for N f = 2, 3, and 4 quark flavors [35,36,37]. Moreover, the LSM is coupled with the Polyakov loop fields, known as the PLSM, to include the interaction and dynamics of the colored gluons. We have a solid term in developing the PLSM to obtain reliable results. For instance, in estimating the features of the moments in thermal QCD medium [38], obtaining the thermal spectrum for masses of meson states and QCD equation-s-tate (EoS) in thermal and dense QCD medium at (non)-zero magnetic back-grounds [39,40,41], Furthermore, the magnetic properties of the QCD matter such as the magnetization, magnetic susceptibility and inverse magnetic catalysis could be estimated [42,43,44]. We have also improved the PLSM to study the thermal structure of the transport properties, the bulk and the shear viscosity, the thermal and the electric conductivity of the QCD matter [45,46]. The extension to N f = 4)-PLSM [47,48] is an essential improvement in order to match with the recent lattice QCD simulations.
The present paper is divided into two main sections. Section 2 summarizes the main features of the effective QCD approach, the Polyakov linear-Sigma model. Section 3 outlines the results obtained.

Polyakov Linear-Sigma Model
Assuming that the hadronic degrees-of-freedom (dof) are the colors of their quark constituents, various effective models relying on the chiral symmetry of QCD have been proposed [49]. The linear sigma model (LSM) implements the chiral symmetry in the hadronic sector and incorporates additional partonic (quarks) dof. We briefly summary the basic concepts of the SU(N f ) Polyakov LSM (PLSM) in section 2.2. Various approaches for the Polyakov loop potentials shall be discussed in section 2.3. The inclusion of the magnetic effects in the mean field approximation by means of the Landau quantization shall be elaborated in section 2.4.

SU(N f ) Lagrangian
For LSM with SU(3) L × SU(3) R and N f = 2, 3, 4 quark flavors and N c color dof, the symmetric LSM Lagrangian L chiral = L f + L m , where the fermionic part is given as with µ is an additional Lorentz index [50], g is the Yukawa coupling of the quarks to the mesonic contributions L m = L SP + L V A + L Int + L U(1)A represented to L SP scalars (J P C = 0 ++ ) and pseudo-scalars (J P C = 0 −+ ), L V A to vectors (J P C = 1 − ) and axial-vectors (J P C = 1 ++ ) mesons and L Int being the interaction between them. The b is the impact parameter and R A is the radius of the nucleus. The magnetic field B is expected to be perpendicular to the reaction plane due to the left-right symmetry of the collision geometry. The figure is taken from [7]. Right panel illustrates a schematic QCD phase-diagram.
Equation (3) represents kinetic and potential terms for the scalar meson nonets. The third term -in this expression -gives the explicit symmetry breaking, which is defined in Eq. (11). This part of the Lagrangian creates scalar and pseudo-scalar mesonic states defined in Φ nonets, Eq. (10). Equation (4) represents the vector meson nonets involving explicit symmetry breaking as given in the second term in Eq. (11). It is obvious that the 3 × 3 matrix of the vector meson nonets involves vector and axial-vector fields, Eq. (10). Thus, this part creates vector and axialvector mesonic states and expresses the interactions between (pseudo)-scalar and (axial)-vector as outlined in Eq. (5). Because of the explicit and spontaneous symmetry breaking, an anomaly term L U(1)A in SU(3) r × SU(3) ℓ is included in the effective Lagrangian. The parameters c, c 0 , c 1 have to be determined, experimentally [56]. The first two terms in this Lagrangian approximate the original axial anomaly term [57,58], while the third term, which is proportional to the first term, is a mixed one. The first anomaly term, in which other terms are used to compare with other effects of different anomaly terms on the hadronic structure [55], is the one taken into account in the present calculations.
In order to repreoduce the related experimental results, the higher-order terms with local chiral symmetry have been included in [56]. It is worthy highlighting that L U(1)A symmetry is anomalous [59] and known as the QCD vacuum anomaly [59,60], i.e. broken by quantum effects. Without this term, a ninth pseudo-scalar Goldstone boson corresponding to the spontaneous breaking of the chiral U(3) ℓ × U(3) r symmetry likely unfold [59,60]. Therefore, the anomaly term is essential and the local chiral symmetry would not cause further numerical problems, at a mass scale of 1 − 2 GeV [56]. The constraint terms are conjetured to have great influences [56]. Thus, it is assumed that the L U(1)A problem is effectively controlled by the inclusion of the c-term. Also, it should be noticed that m 2 (squared tree-level masses of mesons) and m 2 1 have contributions from the spontaneous symmetry breaking [56].
The inclusion of scalar and vector meson nonets in the Lagrangian of PLSM is possible with a redefinition for the contra-covariant derivative of the quark-meson contributions, Eq. (7), where dof of scalar Φ and vector L µ and R µ meson nonets are coupled to A µ , the electromagnetic field. Eqs. (8) and (9), the left-and right-handed field strength tensors, respectively, represent self interaction between vector and axial-vector mesons A µ . Emerging from the globally invariant PLSM Lagrangian, the local chiral invariance requires that g 1 = g 2 = g 3 = g 4 = g 5 = g 6 = g [56] It is obvious that T a =λ a /2 with a = 0 . . . 8 are nine U(3) generators, whereλ a are Gell-Mann matrices with fields Φ of 3 × 3 complex matrix comprising of scalars σ a (J P C = 0 ++ ), pseudo-scalars π a (J P C = 0 −+ ), V µ a , vectors (J P C = 1 −− ), and A µ a axial-vectors (J P C = 1 ++ ) meson states, which are given as The chiral symmetry is explicitly broken by A non-vanishing vacuum expectation value for Φ, Φ = T a σ a breaks the chiral symmetry, spontaneously. Because of the parity is not broken in the vacuum, there are no non-vanishing vacuum expectation values for fields π a . In U (3) V × U (3) A symmetry, these patterns of the explicit symmetry breaking have been obtained as given in ref. [61].
Due to finite quark masses in the (pseudo)-scalar and (axial)-vector sectors, the breaking of U(3) A if H 0 , ∆ 0 = 0 and the symmetry breaking of [61]. The symmetry breaking terms are orig- They are proportional to the matrices H and ∆, Eq. (11). The spontaneous chiral symmetry breaking is conjectured to take place in vacuum state. Therefore, a finite vacuum expectation value for Φ andΦ are assumed to carry the quantum numbers of the vacuum, itself [51]. Thus, the explicit symmetry breaking components (diagonal) h 0 , h 3 and h 8 and δ 0 , δ 3 and δ 8 vanish [51], leading to extracting three finite condensatesσ 0 ,σ 3 andσ 8 . On the other hand,σ 3 breaks the isospin symmetry SU(2) [51]. To avoid this situation, we restrict ourselves to SU(3). This can be N f = 2 + 1 [60] flavors. Correspondingly, two degenerate light (upand down-quarks) and one heavier strange-quark are assumed, i.e. m u = m d = m s , where the violation of the isospin symmetry is neglected. This facilitates the choice of h a (h 0 = 0, h 3 = 0 and h 8 = 0) and for δ a (δ 0 = 0, δ 3 = 0 and δ 8 = 0).

Polyakov loops
The LSM Lagrangian can be coupled to the Polyakov loops [60,62], where the second term, U(φ, φ * , T ), represents the effective Polyakov loop potential [63]. There are various proposals motivated by different underlying QCD symmetries in the pure gauge limit. These are different parameterizations reproducing first-order transition, at T ∼ 187 MeV, N c = 3, and N f = 2 + 1 [20,64,65].
-A simple choice based on Ginzburg-Landau ansatz [20,66]. The underlying Z(3) center symmetry which is spontaneously broken should be conserved. Hence, an expansion in terms of the order-parameter can be expressed as U poly with the temperature-dependent coefficients b 2 (T ) = a 0 + a 1 (T 0 /T ) + a 2 (T 0 //T ) 2 + a 3 (T 0 /T ) 3  U log with the temperature-dependent coefficients a(T ) = a 0 +a 1 (T 0 /T )+a 2 (T 0 /T ) 2 and b(T ) = b 3 (T 0 /T ) 3 . This potential is qualitatively consistent with the leadingorder results from strong coupling expansion [67]. Equation (14) diverges as φ * → 1. This sets on limits to the Polyakov loop variables, i.e. they remain small, especially at high temperatures. a 0 , a 1 , a 2 , and b 3 can be determined from lattice QCD simulations, Tab. 1.
Higher-order terms have been included [68], It should be noticed that if c(T ) and d(T ) vanish, Eq. (15) is reduced to Eq. (14). The various coefficients in Eq. (15) are determined [68] x The different parameters are also summarized in Tab. 1. -A potential inspired by a strong-coupling analysis [65] U Fuku The nearest-neighbor interaction is given in the first term, while the logarithm term is the Haar measure, Eq. (14). There are only two parameters, a determines the deconfinement transition, i.e. the transition temperature in pure gauge theory and b controls the mixing of the chiral and the deconfinement transition. At deconfinement temperature T 0 ≃ 270MeV, a = 664MeV and b = 196.2MeV 3 .
In the present work, we have utilized the higher-order parameterization of the Polyakov loop fields based on the alternatively-improved extension of φ and φ * , for instance, the polynomial-logarithmic Polyakov loop potentials Eq. (15).

Landau quantization
The quantity 2n + 1 − σ can be replaced by a summation over the Landau Levels 0 ≤ ν ≤ ν max f . The earlier is the Lowest Landau Level (LLL), while the latter stands for Maximum Landau Level (MLL), i.e. ν max . For the sake of completeness, we recall that 2 − δ 0ν represents degenerate Landau Levels and ν max f contributes to the maximum quantization number, i.e. ν max f → ∞).
where the brackets represent floor of enclosed quantities and the parameter τ is related to µ, at varying T , i. e. for results given in varying T , we should take into consideration τ ≡ µ f . Concretely, when analyzing the results in thermal medium, µ f is given by τ f , while when analyzing the results in dense medium, T is given as function of τ . Equation (17) refers to the contribution of single and double degenerates for the upper Landau levels. For the present study, we merely need to highlight that various works should have been analyzed [69,70,71,72,73]. Their results in MLL, ν max f , at different temperatures and densities, for instance, can be proposed in terms of the medium parameters such as µ, T and eB. For the present calculations, we have to determine the maximum LL in order to infinity the contributions of MLL. We assume MLL as ν max f , where f → ∞. For details, interested readers are kindly advised to consult ref. [44].
At finite magnetic field, the dispersion relation contributes the Landau Level, so that Thus, the dispersion relation itself is modified by a quantization number, n, known as the Landau quantum number. σ is related to the spin quantum number, Σ = ±S z /2 and m f (q f ) is the quark flavor mass (charge). The present study counts for different contributions from Landau levels. The chiral condensates and deconfinement order-parameters shall be analyzed in a wide range of temperatures, baryon chemical potentials, magnetic fields, so that the chiral QCD phase-diagram could be mapped out in various directions.
A more challenging question is how the Landau quantization is fixed by the magnetic field, let us consider, for example, a Fermi sphere of quarks, where the system is considered as discrete energy levels with respect to the momentum space. At magnetic field background the all spin directions are aligned though the transverse magnetic field, B = Bê z "polarization". With increasing the magnetic field, the order of Fermi energy increases and the energy levels are discretizied. This is called the "Landau levels". Moreover, the phenomena of magnetic catalysis is mainly defined as an enhancement of the dynamical symmetry breaking by an external magnetic field. For instance, upon increasing the magnetic field tends to enhance or "catalysis" the quark-antiquark condensate. This chiral condensate is strongly associated with the breaking of chiral symmetries and also creates masses. Apparently, this is magnetic catalysis [44].
The phenomena of (inverse)-magnetic catalysis are reviewed in great detail in ref. [28,30], where the magnetic field enhances the spontaneous symmetry breaking and the chiral critical temperature decreases as increases the magnetic field strength, i.e. inverse magnetic catalysis.

Mean-field approximation
For a spatially uniform system in a thermal equilibrium, at finite temperature T and finite quark chemical potential µ f , where f stands for u, d and s quarks, the partition function can be constructed. The grand canonical partition function governs the change in numbers of particles and antiparticles. Then, a path integral over quark, antiquark and meson fields leads to [60] in which, we have abbreviated the space-time integration as The integration runs over imaginary time τ = i t from 0 to β = 1/T . For a symmetric quark matter, uniform independent chemical potentials are imposed µ f ≡ µ u = µ d = µ s [66,60,49]. where µ B , µ I and µ Y are the baryon, isospin and hypercharge chemical potentials, respectively. Converting the condensates σ 0 and σ 8 into a pure nonstrange σ l and a pure strange σ s quark flavor leads to [74] Then, the purely mesonic potential is given as The quarks and anti-quark contributions can be divided into two regimes: -At vanishing magnetic field (eB = 0) but finite temperature (T ) and baryon chemical potential (µ f ), where the E f = [p 2 + m 2 f ] 1/2 is the flavor-dependent single-particle energies of quark mass (f = l, s, c) as, where the σ-field coupled through Yukawa coupling g and the subscript l being the degenerate light up and down quarks. The SU(3) Polyakov linear-sigma model (PLSM) in mean-field approximation is utilized in analyzing the thermodynamic properties of quark matter in thermal and dense QCD medium at finite isospin asymmetry in [75]. -At finite magnetic background (eB = 0) where the magnetic field B = Bê z , all the spin directions should be aligned in parallel with the magnetic field, and the momentum directions are only determined according to chirality. The concepts of Landau quantization and magnetic catalysis, where the magnetic field is assumed to be oriented along the z-direction [76,77]. One can notice that the spin (Polarization) plays an essential role in the dimensional reduction of Dirac particles such as quarks. This latter process is known as dimension reduction or magnetic catalysis effect [28,30].
Moreover, the thermodynamic potential should be implemented in a finite temperature, T , chemical potential, µ B , and B as, Accordingly, one can solve the Dirac equation to obtain the modified dispersion relation, Eq. 18. At finite volume, V , and finite magnetic field, eB, the free energy is defined as in which,the LSM mesonic potential U (σ l , σ s ) counts for the contributions of the valence quarks. At very low temperatures, this part of the potential could be excluded, especially in the regime of temperatures typical for the QCD phase transition [38] At vanishing magnetic field B = 0, the second term vanishes, Eq. (26), while δ = 1 in the third term, i.e. Ωq q (T, µ f ), Eq. (24), should be evaluated with the standard dispersion-relation, E f . Therefore, Eq. (28), can be given as When the magnetic field is switched on, δ = 0 and therefore the fourth term vanishes, as well. In this case, Eq. (28) can be reduced to The potential of quark and antiquark contribution, at finite magnetic filed, Ωq q (T, µ f , B), can be divided into two regimes as shown in Eq. (26): 1. The first integral in Eq. (26) refers to the contribution of the magnetic field in zero Landau level, ν = 0. In this case, the modified dispersion-relation, Eq (18), tends to the standard one, E f . 2. The second integral in Eq. (26) gives the contribution of the magnetic field for the upper Landau levels, ν = 1 → ∞, where the dispersion relation will be modified, E B,f , as given in Eq. (18).
Furthermore, the excluding of U (σ l , σ s ) potential from the free energy, especially, at high temperature, tends to maintain dominant contributions of the sea quarks than the valence quarks. In order to evaluate the expectation values of the PLSM order-parameters In nonzero chemical-potential (µ = 0) and finite Ployakov loop variables, the PLSM free energy, at finite V , Eq. (28), is complex. A minimization of a such function would be seen as void of meaning. An analysis of the order parameters is given by minimizing the real part of the PLSM free energy (Re F ). In principle, the (thermal) expectation values of the Ployakov loopφ and its conjugateφ * must be real quantities as discussed in ref. [78]. The solutions of these equations can be determined by minimizing the real part pf the PLSM free energy (Re F ) at a saddle point. The remaining parameters are the chiral order-parametersσ l ,σ s and the Polyakov-loop expectation valuesφ,φ * as functions of T , µ and eB.

Chiral and deconfinement order-parameters
The estimation of the chiral condensates (σ l and σ s ) and the deconfinement order-parameters (φ and φ * ) in dense and thermal medium should be first computed by minimizing the free energy, Eq. (31). The parameters of PLSM, as discussed in sections 2.2 and 2.3, are estimated, at mass of (vacuum) sigma-meson σ = 800 MeV, are measured (vacuum) light and strange chiral condensates are σ lo = 92.5 MeV and σ so = 94.2 MeV, respectively. Figure 2 presents both PLSM chiral-condensates and order-parameters, at finite baryon chemical potentials and magnetic fields eB = m 2 π (upper panel) and eB = 10 m 2 π (bottom panel) as functions of T . In the left-hand panels, (a) and (c) for normalized chiral-condensates, σ l /σ lo and σ s /σ so , we notice that the chiral condensates are shifted to lower values with increasing baryon chemical potential. This means that the chiral critical temperature (T χ ) decreases with increasing eB and with increasing µ, as well. Therefore, we can draw a conclusion that the effect of strong magnetic field is almost the same as that of the baryon chemical potentials, namely on decreasing the phase transition (crossover).
In right-hand panels (b) and (d) for the Polyakov loop parameters, φ and φ * as functions of T , at diffident baryon chemical potentials, we find that both order parameters become differentiable, at finite µ; µ = 100 MeV (dashed curves) and µ = 200 MeV (dotted curves), while at µ = 0 MeV (solid curves), φ = φ * . The deconfinement critical temperature, T φ , is shifted to lower values as the magnetic field and baryon chemical potential increase. We thus draw the conclusion that the φ and φ * have opposite dependence on temperatures. The latter increases as the magnetic field and the baryon chemical potential µ increase, while the earlier decreases.  Fig. 2 but for different values of the magnetic field, at fixed values of baryon chemical potential. We notice that the chiral critical temperature deceases and the crossover phase-transition becomes sharper with increasing magnetic field. This can be interpreted due to the maximum occupation of the Landau levels (ν max → ∞). We conclude that the phase transition seems to be of first order, whenever the chiral condensate passes through a metastable phase, in which light quarks become massless and move freely.
It would be noticed that if one compares the curve for eB = 0.1 GeV 2 with the curve of eB = 10m 2 π ≃ 0.2 GeV 2 , which would be nearly similar values, one would find a slight difference. This is an artifact. It comes from the effect of the occupation Landau levels. The Landau levels and the quantization number are determined for the medium parameters such as T , µ, eB.
As discussed, PLSM is well-suited to study the chiral limit. The inclusion of magnetic field can be achieved by modifying the dispersion relation of quarks and antiquarks, Eq. (24). In doing this, the dimension of the momentum-space should be reduced from three to one and scaled via quark charge and magnetic field. This process is known as magnetic catalysis effect of dimension reduction [28]. Furthermore, the introduction of the mag-   Figure 4 shows the effects of the occupation number of the Landau levels on the temperature dependence of the chiral condensate, σ l and σ s , at a finite magnetic field eB = 0.2 GeV and a vanishing baryon chemical potential. We observe that the change in the Landau levels is relatively significant only in crossover phase-transition region and seems to disappear otherwise. At M LL = 14 (solid curves), 51 (dotted curves) and ∞ (dashed curves), the normalized chiral condensates for light and strange quarks are analyzed as functions of temperatures, at finite magnetic field and vanishing chemical potential. We conclude that the increase in the Landau levels very slightly sharp-ens the phase transition or the crossover and decreases the critical temperature T χ .
Furthermore, the temperature dependence of the chiral condensates, σ l and σ s , and the deconfinement orderparameters, φ and φ * , are estimated, at different values of magnetic fields (eB) and baryon chemical potential (µ). in the chiral condensates around the chiral phase-transition similar to what was observed in a previous study from PLSM without magnetic field [39]. Secondly, a small sudden drop in the chiral condensates referring to first-order phase-transition takes place. There is a gap difference between light and strange chiral condensates, at very high density. This could be understood because of the inclusion of the anomaly term in Eq. (23), where the fit parameters are accordingly modified [60,39]. This was conjectured as an evident on numerical estimation of the chiral condensates. The magnetic field seems to have a sudden drop in the chiral phase-structure. This causes an increase in the critical temperatures. We can draw a conclusion that increasing magnetic field tends to sharpen the phase transition and to accelerate the formation of metastable phase. The right-hand panels (b) and (d) give the deconfinement phase-transition in dependence on magnetic fields; eB = 1 (solid curves), 10 (dashed curves), 15 (dotted curves), 20 (dot-dashed curves), and 25 m 2 π (double dotted curves). The increase in temperature increases the deconfinement phase-transition to larger baryon chemical potential. The thermal and magnetic effects of the hadronic medium on the evolution of Polyakov loop parameters seem to be very smooth. The slope of φ and φ * to the baryon chemical potential can approximately be estimated. It is obvious that the magnetic field decreases these slopes (increases the critical temperature), while the temperature increases them (decreases the critical temperature).

Thermodynamics
We start this analysis with a brief introduction to the basic quantities of thermodynamics, at non-vanishing magnetic field strength. For a statistical system in equilibrium with volume V , temperature T and chemical potential µ, the grand-canonical density operator,ρ, the grandcanonical partition function, Z(T, V, µ), and the grand potential, Ω(T, V, µ) can be introduced in the natural units κ B = = c = 1. The pressure as a function of finite magnetic field eB reads P (T, µ, eB) = −Ω(T, µ, eB), which enables us to characterize the phenomenology to the strong interacting QCD matter, at finite magnetic fields, as the case in HIC, including magnetization, magnetic susceptibility and permeability. Then, the free energy density can be written as [79] where F = −T log ·Z, the total energy density ǫ tot = ǫ + ǫ field , which in turn is divided into two terms; one for the energy density of the medium ǫ and another one of the magnetic field ǫ field = eBM. We notice that the partition function in vanishing magnetic field is given by an integral over six-dimensional phase-space. The dispersion relations follow the Lorentz invariance principle. But, in finite magnetic field, the integral dimensional is reduced and simultaneously accompanied by a considerable modification in the dispersion relation.
The velocity of a test particle with momentum ∂ǫ tot /∂P in finite magnetic field B can be expressed as where σ = ±S/2. Then, the causality is guaranteed for v p not-exceeding the speed of light c, i.e. as long as the B-term is finite positive, which should be quantitatively estimated as a function of temperature and magnetic field. From Eq. (34), the entropy density and the magnetization, for instance, could be derived as Other thermodynamic quantities, such as, pressure can be derived, as well. Because the magnetic field is conjectured to mark a preferred direction, P i might be different along the geometry effect of the magnetic field. As V = L x L y L z and the magnetic field is distributed along z direction, we can distinguish between two different systems.
-B-scheme, in which the magnetic field remains constant in direction, results in an isotropic pressure ϕ-scheme, which sets up magnetic flux ϕ = eB · L x L y , remains constant and results in an anisotropic pressure Accordingly, the thermodynamic quantities should be modified. For instance, in Φ-scheme, the trace anomaly (the interaction measure) reads The speed of sound squared, at constant entropy, is given as [80], where the specific heat, c v , gives the thermal rate change of the energy density, at constant volume. In finite magnetic field, Stefan-Boltzmann (SB) limits can be deduced from lowest-order perturbation theory [81] T log Z(V, T, µ, eB) = 19 V π 2 36 T 4 The normalized entropy density, s/T 3 , which is derived from p z with respect to T , vanishes at T = 0. This may be understood from the fact that the vacuum contribution is conjectured to be a pure quantum effect. This emerges from the interaction of virtual quarks with the external field. Thus, it likely doesn't contribute to entropy [81]. We also notice that, at T > 0, the magnetic field changes the thermal distributions and is expected to modify entropy, We find that near the T c -regime, s/T 3 excellently agrees with the lattice QCD calculations. This might not be also the case, at higher temperatures. shows the same as in the left-hand panel but here P ǫ tot . 6. We find a reasonable agreement with the lattice QCD simulations The phase transition seems to smoothly take place. The temperature dependence continues even above T c . s/T 3 keeps its increase with increasing T /T c so that it becomes slightly lower than the lattice results.
The right-hand panel of Fig. 7 (b) presents the PLSM calculations for p/ǫ, at eB = 0.0 (solid curve), eB = 0.1 (dotted curve) and eB = 0.3 GeV 2 (double-dotted curve). We also compare with recent lattice QCD [81] (open square), (close square) and (open triangle), respectively. A fair agreement is also obtained, especially at eB = 0.1, at low temperature. It is obvious that such an agreement could be improved with increasing magnetic fields. Now, we are able to estimate some fundamental properties of the strongly interacting QCD matter in finite magnetic field such as the magnetization, the magnetic susceptibility, and the permeability. The response of QCD matter to an external magnetic field can be estimated from the free energy density. The magnetic susceptibility with proper renormalization has been introduced in litrature [82]. The magnetic permeability measures the ability of the QCD matter to generate magnetic field or the ability to store magnetic potential energy, which is proportionally constant for the magnetic flux. The magnetic flux -in turn -is produced from influences of the magnetic field. The magnetic permeability can be calculated along the magnetic field on the transverse direction to momentum space p z .
In thermal, dense and magnetic medium, the partition function ln Z is to be properly modified, from which the magnetization can be derived, Eq. (36). The sign of magnetization defines an essential magnetic property, namely whether QCD matter is para-or dia-magnetic, i.e. M > 0 (bara) or M < 0 (dia).
-For dia-magnetized QCD matter, the color charges align oppositely to the direction of the magnetic field and produce an induced current, which spreads as small loops attempting to cancel the effects of the applied magnetic field, and -For para-magnetized QCD matter, most color charges align towards the direction of the magnetic field. As discussed, the magnetization greatly affects the thermodynamic properties of the QCD matter, as the magne- tization measures the response of the system of interest to finite magnetic field. The latter is likely extremly generated in HIC, at least in very short time intervals. Aslo, because of the relativistic, off-center motion of spectators, i.e. peripheral collisions, the rapid motion of electric charges generates magnetic field perpendicular to plane of both motion direction and the electric field. Also, because of local imbalance in the momenta carried by the colliding nucleons in peripheral and central collisions. This local imbalance leads to an angular momentum and thus a magnetic field [83]. As discussed, such a magnetic field is typically very huge, O(m 2 π ). It largely exceeds the detector's magnet field.
The magnetization, M, can be derived from Eq. (28). The values obtained can be given in GeV 2 in the natural units. The sign of M refers to para-or dia-magnetic QCD matter. If M > 0 or M < 0, the QCD matter is either para-or dia-magnetized, respectively.
In Fig. 8, M is given as a function of T , at eB = 0.1 (dotted), 0.2 (dashed), 0.3 (double-dotted) and 0.4 GeV 2 (dash-dotted curve) at vanishing µ. The results are compared with recent lattice QCD [81] at eB = 0.1 (closed square), 0.2 (circle), 0.3 (triangle) and 0.4 GeV 2 (astride). It is obvious that M > 0 and increases as the magnetic field increases. This result indicates that paramagnetic properties of the QCD matter. At temperatures below the critical value, the PLSM results resemble the lattice data in an excellent way. At temperatures characterizing QGP (higher than the critical temperature), the PLSM curve becomes larger than the lattice data, especially at very high temperatures. In this range of temperatures, the colorless hadrons are conjectured to deconfine into colored quarks and gluons. The discrepancies suggest that the corresponding dof are not sufficient enough to achieve a good agreement, especially at very high temperature. Furthermore, we notice that the PLSM calculations give an evident on paramagnetic features of the hot QCD matter.
Also, the magnetic susceptibility and permeability reflect the magnetic response of the hot QCD matter. In other words, the response of QCD matter to the magnetic field can be determined by the slope of M with respect to the magnetic field. The second derivative of the free energy density with respect to the finite magnetic field is the magnetic susceptibility, In response to the magnetic field, the magnetic susceptibility is a dimensionless proportionality parameter indicating the degree of magnetization of the QCD matter. Furthermore, the relative magnetic permeability, µ r , normalized to the vacuum magnetic permeability µ 0 can be translated as the magnetic effect in thermal QCD medium. This can be determined by different methods. With a direct relation to the magnetic susceptibility, we have As shall be introduced, this relation agrees well with the lattice QCD simulations, in which the magnetic permeability is expressed in terms of the magnetic susceptibility, where α m = e 2 /4π is the fine structure constant. This expression distinguishes between dimensionless proportionality constants, namely external B ext and induced magnetic field B ind . We highlight that the higher-order permeability seems to be limited by the magnetic susceptibility, which is given by the reciprocal of the square of elementary charge e, i.e. χ B → 1/e 2 when µ → ∞.
Left-hand panel of Fig. 9 (a) shows the magnetic susceptibility as functions of temperature, at eB = 0.0 GeV 2 and vanishing µ. The results deduced from PLSM are compared with various lattice simulations (symbols), in which different calculation methods and algorithms are applied. We also compare with the resonance gas model (HRG) model. The recent lattice QCD simulations [81] (open circle) were estimated by using half-half method in 24 3 × 32 lattice (closed triangle) and integral method in 28 3 ×10 lattice (open triangle). The diamonds represent lattice simulations for N f = 2 + 1 and when using HISQ/tree action with light quark masses m l /m s = 0.05 and temporal dimension N τ = 8, the lattice results [14]. The closed circles stand for results obtained from isotropy lattice [84].
Right-panel of Fig. 9 (b) depicts the relative permeability with respect to that of the vacuum compared to recent lattice QCD calculations (open triangles) [81] in a wide range of temperatures, at eB = 0.0 GeV 2 and at vanishing µ. It is apparent that the agreement between PLSM and lattice QCD calculations is excellent.
Features of PLSM and lattice QCD results can be summarized as follows.
-The magnetic susceptibility obtained from the HRG model [81] (dashed curve) confirms the nature of the QCD matter as dia-magnetized, especially at low temperatures. Here, the free energy density is the sum over all contributions from the colorless hadrons and their resonances contributes to the hadron interactions in order to assure negative magnetic susceptibility [81]. -In PLSM, the free energy density is divided into three terms. The first one is the pure mesonic potential which is obtained from the Lagrangian for the pure gauge. The second one gives the contributions of the quark and antiquarks, which -as the name says -have mesonic fluctuations from the quarks and the antiquarks. The third term represents the interactions of the color charges and the gluons. This means that there two types of contributions to the hadronic fluctuations, while only one type contributes to the gluon interactions. -At very low temperatures, the slope of χ(T ) is apparently negative (inside-box in left-hand panel of Fig.  9) [85]. This gives a signature that the QCD matter is dia-magnetized and also well reproduces the different lattice simulations. When switching to the high temperature regime, i.e. restoring the broken chiral symmetry, we observe a transition between dia-and paramagnetic properties, where QCD para-magnetism is likely, at high temperature. -We notice that the PLSM results confirm that the strongly interacting QCD matter has para-magnetic properties. The magnetic susceptibility steeply increases when increasing temperatures towards the deconfinement phase-transition. These conclusions have been found in a wide range of temperatures 100 ≤ T ≤ 250 MeV [86,87].

QCD phase-diagram at finite magnetic field
There are two different mechanisms manifesting the influences of the magnetic field on the QCD phase-diagram. The first one is that the magnetic field improves the QCD phase-transition due to its contributions to the Landau quantizations or the Landau levels. The second one is that the magnetic field contributes to the suppression in the chiral condensate due to the restoration of the broken chiral symmetry. This suppression is known as inverse magnetic catalysis. It is a decrease in the chiral critical temperature with the increase in the magnetic field. The PLSM has two main types of the order parameters; the chiral condensates, which are connected with two light quarks, σ l , and one strange quark, σ s , and the Polyakov loop variables, φ and φ * . The intersects of φ and φ * with the chiral condensates are used in determining the quasi-critical temperatures (dotted curve) in Fig. 10. The solid curve is estimated from higher-order moments of quark multiplicity (not shown in the present paper), where the critical temperatures are estimated at the peak of normalized quark susceptibility χ q /T 2 . From both methods, we find that the critical temperatures decrease with increasing magnetic fields referring to an inverse magnetic catalysis.
In Fig. 10, the influences of the magnetic field on the chiral QCD phase-diagram are depicted. Left-hand panel shows the variation of the critical temperature with increasing magnetic field, at almost vanishing µ. The PLSM calculations are compared with recent lattice simulations given (circles with errorbars) [81]. The vertical bands refer to the magnetic fields generated, at RHIC (orders of per cent of GeV 2 or ∼ m 2 π ) and LHC energies (orders of per ten GeV 2 or ∼ 10 − 15 m 2 π ). The PLSM calculations take into consideration two methods in order to estimate the critical temperatures. One is based in normalized quark susceptibility χ q /T 2 and the other one is based on the freeze-out condition s/T 3 , where s is the entropy density.
There is an excellent agreement, especially at 0 ≤ eB [GeV 2 ] ≤ 0.2. The solid curve matches well with the lattice simulations, at a wider range of magnetic fields 0.13 ≤ eB [GeV 2 ] ≤ 0.55. The first method of determining the critical temperature apparently overestimates the lattice results, at low temperature, while the second method slightly underestimates these, especially at high temperatures. We conclude that the chiral magnetic field improves the agreement as it enhances the chiral condensates. This depends also on the type of contributions to the Landau levels introduced to the QCD effective approach.
The middle-panel draws the dependence of the critical baryon chemical potential on the magnetic fields, at T = 50 (solid curve) and 100 MeV (dashed curve). At constant magnetic field as that at RHIC or LHC energy, large critical baryon chemical potential can be reached, even at low temperatures, i.e. µ decreases with increasing eB. The intersection between the deconfinement phasetransition and the light quark chiral condensate is utilized to determine µ, at which the broken chiral symmetry is restored. There is no lattice QCD simulations to compare with. The right-hand panel presents T -µ phase-diagram, at eB = m 2 π (solid curve), 10 m 2 π (dashed curve) and 20 m 2 π (dotted curve). In determining the critical temperature, various methods have been utilized, for example, the intersection between the Polyakov loop variables φ and φ * , which are related to the deconfinement phase-transition, and the chiral condensates of light and strange quarks, σ l and σ s , respectively. The latter is related to the restoration of the broken chiral symmetry. For example, the critical temperature corresponding to the chiral restoration of light quark, T χ l c , can be determined by the intersection between φ and σ l , while the critical temperature corresponding to the chiral restoration of strange quark, T χs c , can be defined from the intersection between φ * and σ s . We observe that increasing magnetic field enhances the chiral QCD phase-diagram, i.e. the chiral phase-transition takes place at lower temperature. The estimation of freeze-out parameters, T and µ, in dependence on the heavy-ion centralities or the impact parameters allows to analyze the influence of the magnetic field, experimentally [9]. There are various experimental results on chemical and thermal freeze-out as reviewed in ref. [88].
In Fig. 11, the chemical freeze-out condition s(T, eB, µ)/T 3 = 7 is implemented in order to estimate the feeze-out temperature T . The entropy density is calculated, at different temperatures, baryon chemical potentials, and magnetic fields. At a given value of the entropy density normalized to T 3 , the related baryon chemical potential, µ, and the corresponding magnetic field, eB, are determined. These three parameters are then depicted in Fig. 11. It is a multi-dimensional chemical freeze-out boundary showing the dependence the freeze-out diagram (T − µ), which can directly be related to the one analyzed from ethe experimental measurements of various particle ratios [88], on the magnetic field. It is obvious that, at small µ, the effect of the magnetic field is almost negligible. At higher temperatures, the drop in T c around the chiral phase-transition moves to lower values with increasing eB. Again, this refers to an inverse catalysis. At very high temperatures,  Fig. 11. As in Fig. 10 but here the multi-dimensional chemical freeze-out boundary combines the chemical freeze-out parameters T with the magnetic field eB and the baryon chemical potentila µ.
there a slight increase in T with increasing eB. Increasing eB makes the chiral phase-transition smother. It is important to notice that the shape of T -µ phase-diagram looks different from the one at vanishing eB [88].

Conclusions
In relativistic HIC, the off-center motion of spectators in the peripheral collisions and the rapid motion of electric charges generate huge magnetic fields perpendicular to the plane of the motion direction and the electric field. Also, the local imbalance in momenta carried by the colliding nucleons in both peripheral and central collisions generates a huge magnetic field, as well. The response of QCD matter modeled in PLSM to these external magnetic fields such as magnetization, magnetic susceptibility and permeability is determined. The chiral magnetic properties of the QCD phase-diagram are analyzed. The inclusion of finite magnetic field in PLSM can be accomplished by firstly modifying the dispersion relation of the quarks and antiquarks, so that the dimension of the momentum-space is reduced from three to one and then scaled via quark charge and magnetic field (magnetic catalysis) and secondly Landau quantization should be integrated in. The latter modifies the dispersion relation by a quantization number; the Landau quantum number. σ, which is related to the spin quantum number and to the masses of quark flavors. Both chiral condensates and deconfinement order-parameters have been analyzed in a wide range of temperatures, baryon chemical potentials, magnetic fields, so that the chiral QCD phase-diagram could be mapped out in various directions.
The magnetization remarkably affects the thermodynamic properties of the QCD matter. The PLSM results give an evident on paramagnetic features of the hot QCD matter. The magnetic susceptibility, at low temperatures, is negative indicating that the QCD matter is dia-magnetized. At higher temperatures regime, i.e. restoring the broken chiral symmetry, there is a transition from dia-to paramagnetic properties. The permeability is a characterizing property of a magnetic material measuring the ability to create magnetic field and to store magnetic potential energy. The latter is proportionally constant for magnetic flux, which is produced from the influences of the magnetic field. At low temperatures (hadron phase) the QCD permeability is small but rapidly increases around the deconfinement phase-transition. At high temperatures (parton phase) the QCD permeability is large. We conclude that in a wide range of temperatures, the magnetic permeability normalized to the vacuum value agrees well with the lattice QCD simulations, so that this quantity can be utilized as a magnetic order-parameter.
To estimate the variation of the chiral critical temperature with the magnetic field, the PLSM results are confronted to the recent lattice QCD data. We found that the chiral critical temperature decreases as the magnetic field increases. We conclude that the magnetic catalysis of the thermal QCD medium is inverse and an inverse interrelation between the chiral chemical potential and the magnetic field is obtained. The chiral phase-diagram is shifted to lower temperatures due to the increase in the magnetic field. This result is confirmed by two different methods, one from the thermal and dense phase-transition of the PLSM parameters and the another one by applying a condition from freeze-out parameters s/T 3 .