Influence of finite volume effect on the Polyakov Quark-Meson model

In the current work, we study the influence of a finite volume on $2+1$ $SU(3)$ Polyakov Quark-Meson model (PQM) order parameters, (fluctuations) correlations of conserved charges and the quark-hadron phase boundary. Our study of the PQM model order parameters and the (fluctuations) correlations of conserved charges indicates a sizable shift of the quark-hadron phase boundary to higher values of baryon chemical potential ($\mu_{B}$) and temperature ($T$) for decreasing the system volume. The detailed study of such effect could have important implications for the extraction of the (fluctuations) correlations of conserved charges of the QCD phase diagram from heavy ion data.


I. INTRODUCTION
One of the major aims of the current heavy ion collisions research is to study the properties of the strongly interacting matter created in such collisions theoretically and experimentally. On the experimental level, many facilities have been designed to investigate the strongly interacting matter phase diagram such as Relativistic Heavy Ion Collider (RHIC) program [1], and the Large Hadron Collider (LHC) [2]. Different studies suggest that the strongly interacting matter phase transition from hadronic phase to quark-gluon plasma (QGP) phase be a smooth crossover at low density and high-temperature [3], and first-order phase transition at high density and low temperature [4,5]. Both smooth crossover and first-order phase transitions expected to be connected by the critical endpoint (CEP), at which the phase transition is expected to be second order. One avenue to map out and study the QCD phase diagram is through the effective models such as the quark-meson (QM) model [6][7][8], the Nambu-Jona-Lasinio (NJL) model [9], and their Polyakov-loop extended versions [10].
Many studies have been devoted to investigating the QCD phase diagram, higher order moments and the thermodynamics of two [11,12] and three quark flavors [13] QM model and even PQM model with different Polyakovloop potentials. The thermodynamic properties (pressure, the equation of state, the speed of sound, specific heat, trace anomaly, and the bulk viscosity) have been evaluated at finite and vanishing chemical potential [13,14].
The effect of a finite-volume on the strongly interacting matter has been widely studied [15][16][17][18][19][20][21][22]. Those studies include that finite volume has a strong effect on the transition temperature (T c ), the location of the critical endpoint and other thermodynamic properties. In PQM model the transition temperature (T c ) shifted to large values as the volume decrease [21] and the location of the critical endpoint is shifted toward large µ and small T [17,18]. On another hand NJL [16] and PNJL [19] indicates that the transition temperature (T c ) shifted to small values as the volume decrease and the location of the critical endpoint is shifted toward large µ and small T for (2+1 flavor) and toward small µ and small T for (2 flavors).
In this work, we investigate the effect of the finite volumes on the PQM model order-parameters, phase-transition, and the conserved charges fluctuations and correlations. The present work is organized as follows. In section II we give a brief overview of the PQM model. The PQM model calculations of the order-parameters, thermodynamic properties, and the conserved charges fluctuations and correlations are compared with the LQCD [23,24], also, the influence of finite-volume effect on the PQM model conserved-quantities, baryon, charge, strangeness and them correlations will be presented in section III. We conclude with a summary and an outlook in section IV.

II. THE POLYAKOV QUARK MESON (PQM) MODEL
The SU (3) Quark Meson model with N f = 2 + 1 flavor quarks, coupled to Polyakov loop dynamics to formulate the Polyakov Quark Meson (PQM) model [13]. The related Lagrangian is given as; where the chiral part of the Lagrangian, L chiral = L quark + L meson , has SU (3) L × SU (3) R symmetry [25,26]. The first part provides the fermionic sector, and the second part represents the mesonic contribution, both contributions had been extensively discussed in Ref. [14].
The second term in Eq. (1), U(φ, φ * , T ), represents the Polyakov-loop effective potential [27], which is expressed by using the dynamics of the thermal expectation value of a color traced Wilson loop in the temporal direction Then, the Polyakov-loop potential and its conjugate read: where P is the Polyakov loop. This can be represented by a matrix in the color space [27] P where β = 1/T is the inverse temperature and A 4 = iA 0 is called Polyakov gauge [27,28].
In case of no quarks, zero quark chemical potential, then φ = φ * and the Polyakov loop is recognized as an order parameter for the deconfinement phase-transition [29]. In the present work, we used Polyakov loop effective potential U (φ, φ * , T ) as discussed in Refs. [29,30] but with a new dimensionless parameter K that help us get a better agreement with the LQCD. Other Polyakov loop potentials [31,32] were also examined in various work. However, the particular selection made for this work does not affect the main conclusions of our work.
where B and K are dimensionless parameters given as: The mean field approximation is used following Refs. [14,33] to obtain the grand potential as: where σ x and σ y are the non-strange and strange chiral condensates, the first term in Eq. (9) is a purely mesonic potential expressed as: Here, m 2 , h x , h y , λ 1 , λ 2 and c are model parameters as reported in Ref. [26]. The parameters values used in the current study, are listed in Table. I below. Different studies [34,35] indicate that extending the PQM model with the vector meson sector will help to accomplish better agreement with LQCD at T < T c . Such correction was not included in this work and shall be discussed in future work.
The third term in Eq. (9) Ωψ ψ (T, µ f ; φ, φ * ) which gives the quark and anti-quark contributions can be shown as [13], where N gives the number of the quark flavors, E f = p 2 + m 2 f (the index f runs over different quark flavors (u, d and s)) is the dispersion relation, energy, of the valence quark and antiquark. Assuming degenerate light quarks, q ≡ u, d, then we can give the masses as follows: The quark chemical potentials µ f are related to the baryon (µ B ), strange (µ S ) and charge (µ Q ) chemical potentials via the following transformations [36]; The influences of a finite volume are introduced in the PQM model by following the approximate method illustrated in [19,37] via a lower momentum cut-off p min [GeV ] = π/R[GeV ] = λ, where R is the length of a cubic volume. In this analysis, we are studying a simple situation (lower momentum cut-off ). A full implementation of the finite volume would require decent consideration of the effects of the surface and curvature, as well as boundary conditions which are periodic for bosons and anti-periodic for fermions. This full implementation of the boundary conditions leads to an infinite sum over discrete momentum values.
The PQM model has a set of parameters discussed in Refs [26,30] and listed in Tabs.(I,II).  In order to estimate the model different parameters, σ x , σ y , φ and φ * , we minimize the thermodynamic potential, Eq. (9), with respect to σ x , σ y , φ and φ * which gives us a set of four equations of motion: where σ x =σ x , σ y =σ y , φ =φ and φ * =φ * are the global minimum.

III. RESULTS
In this section, we will discuss our PQM model calculations using the parameters summarized in Tables. (I,II) to illustrate the effect of finite volume on the model order parameters, phase transition and fluctuations/correlations of the conserved charges.

A. Order parameters and phase transation
In the following, we present several calculations to illustrate the impacts of the finite volumes on the PQM model order parameters and the chiral phase transition. Figure. 1 shows the thermal dependence of the non-strange and strange chiral condensates (σ x , σ y ) panels (a, b and c) and the Polyakov loops (φ and φ * ) panels (d, e and f) for different volume selections and different µ B values. The upper panels show that both σ x and σ y increase as the system volume is decreased, with larger sensitivity for the non-strange chiral condensates (σ x ). The lower panels show very little if any, volume dependence for φ and φ * at different µ B .
As pointed out PQM model contains strange and non-strange chiral condensates which reflect the chiral phase transitions. Using both chiral condensates we can investigate the finite volume effects on the SU (3) 2 + 1 PQM model chiral phase transition via the normalized net-difference condensate ∆ ls (T ) as defined in Ref. [38], where h x (h y ) are non-strange (strange) explicit symmetry breaking parameters.    The study of the phase diagram of the PQM model for at fixed volume could be done through mapping out the µ B dependence of ∆ q,s (T ). For a fixed R and µ B values, d∆ ls /dT will peak up at a particular point expressing the phase transition. Therefore, the phase diagram can be studied by outlining such points for a wide range of baryon chemical potentials. Fig. 3 illustrates the effects of finite volume on the phase diagram. The parameters T c and µ Bc represent the transition temperature at µ B = 0.0 GeV and the transition chemical potential at low temperature respectively at R = ∞. Our calculations reveal that the PQM model phase diagram in the (µ B , T )-plane, increases with decreasing the system volume. For the R = 2.0 (fm) the µ B value at low temperature increased by about 30% and the T value at µ B = 0.0 GeV increased by about 19% from them values at R = ∞ (fm).

B. Fluctuations and correlations of conserved charges
The thermodynamics quantities and (diagonal) off-diagonal susceptibilities can be determined by using the thermodynamic pressure as [24], where superscripts i, j and k run over integers that represent the derivatives orders. The indexes B, Q and S represent the conserved-quantities, baryon, charge, and strangeness, respectively. Eq. (20) illustrates the dependence of the (fluctuations) correlations of conserved charges on the temperature, chemical potential, and system volume. The susceptibilities evaluated first by computing the thermodynamic potential at vanishing µ f and then expand the scaled thermodynamic potential in a Taylor series around µ f /T = 0. Before addressing the system volume effect, it is informative to contrast the PQM model thermodynamics quantities and (diagonal) off-diagonal susceptibilities calculations for (µ B = 0 and R = ∞), to similar results from LQCD calculations [23,24]. Such comparisons are presented in Figs. (4,5); which indicate a good agreement between the PQM model and LQCD [23,24]. These comparisons could be improved spatially at low temperature by including the vector mesons sector to the PQM model. The influence of the finite volume on the model thermodynamics quantities has been discussed in our previous study [21]. Figure 6 displays the temperature dependence of the normalized conserved-fluctuations, baryon (χ BB ), charge (χ QQ ) and strangeness (χ SS ), respectively. The results presented for several volume selections at three µ B values, µ B = 0.0, 0.4 and 0.6 GeV. Our results indicate that the normalized fluctuations decrease with the volume which quickly trends towards the infinite volume value at high temperature. The non-strange susceptibilities (χ 2 BB and χ 2 QQ ) shows a higher sensitivity to the volume change more than the strange susceptibility (χ 2 SS ). This weak sensitivity to the volume change of the strange quantities could be driving from the large mass of the strange quark. Figure 7 shows the temperature dependence of the off-diagonal susceptibilities, χ 2 BQ , χ 2 BS and χ 2 SQ for several volume selections and different µ B values. The net baryons show a high correlation to the net charge and less correlation to the net strange. Our results indicate that the normalized correlations decrease with the volume which quickly trends towards the infinite volume value at high temperature. Also, the non-strange correlation (χ BQ ) show a higher sensitivity to the volume change.
Also, the temperature dependence of the higher order baryon susceptibilities χ n B (n = 4, 6 and 8) for different volume selections at µ B values, µ B = 0.0, 0.4 and 0.6 GeV are shown in Figure 8. The n th -order susceptibilities decrease with the volume selections, and for n = 6, 8 they start to peak around the transition temperature T c . Also, we observed a stronger oscillation in all higher harmonics as we increase the µ B values.

IV. CONCLUSIONS
In this work, we have used the 2 + 1 SU (3) Polyakov Quark-Meson model (PQM) framework to study the properties of the QCD medium produced at finite volume in heavy ion collisions. This model framework provides several conserved-quantities, baryon, charge, and strangeness which compare well with those obtained in LQCD calculations for vanishing µ B . Our calculations indicate that the conserved-quantities (χ ijk BQS ) are significantly influenced by finite volume effects. The calculated conserved-quantities decrease with the volume which quickly trends towards the infinite volume value at high temperature. Also, the non-strange quantities show a higher sensitivity to the volume change more than the strange once. Finally, PQM model conserved-quantities, suggests that the quark-hadron phase boundary is shifted to higher values of µ B and T with decreasing system volume.