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

TL;DR

This paper investigates how finite volume effects influence the phase transition and conserved charge correlations in the Polyakov Quark-Meson model, with implications for interpreting heavy ion collision data.

Contribution

It provides a detailed analysis of finite volume impacts on the PQM model's order parameters and phase boundary shifts, a novel aspect in this context.

Findings

Finite volume causes the quark-hadron phase boundary to shift to higher chemical potential and temperature.

Finite volume effects are significant for interpreting conserved charge fluctuations in heavy ion collisions.

The study highlights the importance of volume considerations in QCD phase diagram analysis.

Abstract

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.

Tables

Table 1. Table 1: Summary of the QM model parameters employed in the presented calculations.

$c$ (MeV)	${\lambda }_1$	$m^2$ ($Me{V}^2$)	${\lambda }_2$	$h_x$ ($Me{V}^3$)	$h_y$ ($Me{V}^3$)
$4807.84$	$1.40$	${\left(342.52\right)}^2$	$46.48$	${\left(120.73\right)}^3$	${\left(336.41\right)}^3$

Table 2. Table 2: Summary of the Polyakov loop potential parameters employed in the presented calculations.

$a_1$	$a_2$	$b_0$	$b_1$	$b_2$	$b_3$
$0.75$	$7.5$	$6.75$	$-1.95$	$2.625$	$-7.44$
$k_1$	$k_2$	$k_3$	$k_4$
$0.30$	$0.25$	$0.24$	$0.20$

Equations

$$\mathcal{L}={\mathcal{L}}_{\mathrm{chiral}}-\mathcal{U}(\varphi ,{\varphi }^{\ast },T),$$

$$\Phi (\vec{x})=\frac{1}{N_c}⟨\mathcal{P}(\vec{x})⟩,$$

$$\varphi$$

$${\varphi }^{\ast }$$

$$\mathcal{P}(\vec{x})=\mathcal{P}\exp \left[i{\int }_0^{\beta }d\tau A_4(\vec{x},\tau )\right],$$

$$\frac{\mathcal{U}(\varphi ,{\varphi }^{\ast },T)}{T^4}$$

$$B$$

$$K$$

$$\Omega (T,{\mu }_f)$$

$$U({\sigma }_x,{\sigma }_y)$$

$${\Omega }_{\bar{\psi }\psi }(T,{\mu }_f;\varphi ,{\varphi }^{\ast })$$

$$m_q$$

$$m_s$$

$${\mu }_u$$

$${\mu }_d$$

$${\mu }_s$$

$${\frac{\partial \Omega }{\partial {\sigma }_x}=\frac{\partial \Omega }{\partial {\sigma }_y}=\frac{\partial \Omega }{\partial \varphi }=\frac{\partial \Omega }{\partial {\varphi }^{\ast }}|}_{min}=0,$$

$${\Delta }_{ls}(T)$$

$$p$$

$$s$$

$$ϵ$$

$${\chi }_{BQS}^{ijk}$$

Full text

Chiral Lagrangian, Quark confinement, Quark-Gluon Plasma

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  Fachini (2006), and the Large Hadron Collider (LHC)  Monteno (2007). 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 Aoki et al. (2006), and first-order phase transition at high density and low temperature Ejiri (2008); Pisarski and Wilczek (1984). 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 Lee (1972); Kovacs and Szep (2007, 2008), the Nambu-Jona- Lasinio (NJL) model Nambu and Jona-Lasinio (1961), and their Polyakov-loop extended versions Fukushima (2004).

Many studies have been devoted to investigating the QCD phase diagram, higher order moments and the thermodynamics of two Kahara and Tuominen (2008); Wambach et al. (2010) and three quark flavors Schaefer and Wagner (2009a) QM model and even PQM model with different Polyakov–loop 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 Schaefer and Wagner (2009a); Mao et al. (2010).

The effect of a finite-volume on the strongly interacting matter has been widely studied  Fisher and Barber (1972); Abreu et al. (2006); Palhares et al. (2011); Fraga et al. (2011); Bhattacharyya et al. (2013, 2015a); Magdy et al. (2017); Almasi et al. (2017). 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 Magdy et al. (2017) and the location of the critical endpoint is shifted toward large $\mu$ and small $T$ Palhares et al. (2011); Fraga et al. (2011). On another hand NJL Abreu et al. (2006) and PNJL Bhattacharyya et al. (2013) 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 $\mu$ and small $T$ for (2+1 flavor) and toward small $\mu$ 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 Borsanyi et al. (2014); Bazavov et al. (2012), 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 Schaefer and Wagner (2009a). The related Lagrangian is given as;

[TABLE]

where the chiral part of the Lagrangian, $\mathcal{L}_{\rm chiral}=\mathcal{L}_{quark}+\mathcal{L}_{meson}$, has $SU(3)_{L}\times SU(3)_{R}$ symmetry Lenaghan et al. (2000); Schaefer and Wagner (2009b). The first part provides the fermionic sector, and the second part represents the mesonic contribution, both contributions had been extensively discussed in Ref. (Mao et al., 2010).

The second term in Eq. (1), $\mathbf{\mathcal{U}}(\phi,\phi^{*},T)$, represents the Polyakov–loop effective potential Polyakov (1978), which is expressed by using the dynamics of the thermal expectation value of a color traced Wilson loop in the temporal direction

[TABLE]

Then, the Polyakov–loop potential and its conjugate read:

[TABLE]

where $\mathcal{P}$ is the Polyakov loop. This can be represented by a matrix in the color space Polyakov (1978)

[TABLE]

where $\beta=1/T$ is the inverse temperature and $A_{4}=iA^{0}$ is called Polyakov gauge Polyakov (1978); Susskind (1979).

In case of no quarks, zero quark chemical potential, then $\phi=\phi^{*}$ and the Polyakov loop is recognized as an order parameter for the deconfinement phase-transition  Ratti et al. (2006). In the present work, we used Polyakov loop effective potential $U(\phi,\phi^{*},T)$ as discussed in Refs. Ratti et al. (2006); Ghosh et al. (2008) but with a new dimensionless parameter $K$ that help us get a better agreement with the LQCD. Other Polyakov loop potentials Haas et al. (2013); Schaefer et al. (2009) were also examined in various work. However, the particular selection made for this work does not affect the main conclusions of our work.

[TABLE]

where $B$ and $K$ are dimensionless parameters given as:

[TABLE]

The mean field approximation is used following Refs. (Mao et al., 2010; Tawfik et al., 2014) to obtain the grand potential as:

[TABLE]

where $\sigma_{x}$ and $\sigma_{y}$ are the non-strange and strange chiral condensates, the first term in Eq. (9) is a purely mesonic potential expressed as:

[TABLE]

Here, $m^{2}$, $h_{x}$, $h_{y}$, $\lambda_{1}$, $\lambda_{2}$ and $c$ are model parameters as reported in Ref. Schaefer and Wagner (2009b). The parameters values used in the current study, are listed in Table. 1 below. Different studies Kovács and Wolf (2015); Kovács et al. (2016) 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) $\Omega_{\bar{\psi}\psi}(T,\mu_{f};\phi,\phi^{*})$ which gives the quark and anti-quark contributions can be shown as Schaefer and Wagner (2009a),

[TABLE]

where $N$ gives the number of the quark flavors, $E_{f}=\sqrt{\vec{p}^{2}+m_{f}^{2}}$ (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\equiv u,d$, then we can give the masses as follows:

[TABLE]

The quark chemical potentials $\mu_{f}$ are related to the baryon ($\mu_{B}$), strange ($\mu_{S}$) and charge ($\mu_{Q}$) chemical potentials via the following transformations Fu (2013);

[TABLE]

The influences of a finite volume are introduced in the PQM model by following the approximate method illustrated in Bhattacharyya et al. (2015b, 2013) via a lower momentum cut-off $p_{min}[GeV]=\pi/R[GeV]=\lambda$, 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 Schaefer and Wagner (2009b); Ghosh et al. (2008) and listed in Tabs.(1,2).

In order to estimate the model different parameters, $\sigma_{x}$, $\sigma_{y}$, $\phi$ and $\phi^{*}$, we minimize the thermodynamic potential, Eq. (9), with respect to $\sigma_{x}$, $\sigma_{y}$, $\phi$ and $\phi^{*}$ which gives us a set of four equations of motion:

[TABLE]

where $\sigma_{x}=\bar{\sigma_{x}}$, $\sigma_{y}=\bar{\sigma_{y}}$, $\phi=\bar{\phi}$ and $\phi^{*}=\bar{\phi^{*}}$ are the global minimum.

III Results

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

III.1 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 ($\sigma_{x}$, $\sigma_{y}$) panels (a, b and c) and the Polyakov loops ($\phi$ and $\phi^{*}$) panels (d, e and f) for different volume selections and different $\mu_{B}$ values. The upper panels show that both $\sigma_{x}$ and $\sigma_{y}$ increase as the system volume is decreased, with larger sensitivity for the non-strange chiral condensates ($\sigma_{x}$). The lower panels show very little if any, volume dependence for $\phi$ and $\phi^{*}$ at different $\mu_{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 $\Delta_{ls}(T)$ as defined in Ref. Schaefer et al. (2010),

[TABLE]

where $h_{x}$ ($h_{y}$) are non-strange (strange) explicit symmetry breaking parameters.

Figure. 2 shows the thermal dependence of the normalized net-difference condensate $\Delta_{ls}$ panels (a, b and c) and $d\Delta_{ls}/dT$ panels (d, e and f) for different volume selections and different $\mu_{B}$ values. The upper panels indicate an increase in $\Delta_{ls}(T)$ as the system volume is decreased. The lower panels show that for fixed values of $R$ and $\mu_{B}$, the $d\Delta_{ls}/dT$ is peaking up at a specific point indicating the phase transition. The peak position is shifted toward lower temperature as the $\mu_{B}$ value increase.

The study of the phase diagram of the PQM model for at fixed volume could be done through mapping out the $\mu_{B}$ dependence of $\Delta_{q,s}(T)$. For a fixed $R$ and $\mu_{B}$ values, $d\Delta_{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 $\mu_{Bc}$ represent the transition temperature at $\mu_{B}$ = 0.0 GeV and the transition chemical potential at low temperature respectively at $R$ = $\infty$. Our calculations reveal that the PQM model phase diagram in the $(\mu_{B},T)$-plane, increases with decreasing the system volume. For the $R$ = 2.0 (fm) the $\mu_{B}$ value at low temperature increased by about $30\%$ and the $T$ value at $\mu_{B}$ = 0.0 GeV increased by about $19\%$ from them values at $R$ = $\infty$ (fm).

III.2 Fluctuations and correlations of conserved charges

The thermodynamics quantities and (diagonal) off-diagonal susceptibilities can be determined by using the thermodynamic pressure as Bazavov et al. (2012),

[TABLE]

[TABLE]

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 $\mu_{f}$ and then expand the scaled thermodynamic potential in a Taylor series around $\mu_{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 ($\mu_{B}=0$ and $R=\infty$), to similar results from LQCD calculations Bazavov et al. (2012); Borsanyi et al. (2014). Such comparisons are presented in Figs. (4,5); which indicate a good agreement between the PQM model and LQCD Bazavov et al. (2012); Borsanyi et al. (2014). 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 Magdy et al. (2017).

Figure 6 displays the temperature dependence of the normalized conserved-fluctuations, baryon ($\chi_{BB}$), charge ($\chi_{QQ}$) and strangeness ($\chi_{SS}$), respectively. The results presented for several volume selections at three $\mu_{B}$ values, $\mu_{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 ($\chi^{2}_{BB}$ and $\chi^{2}_{QQ}$) shows a higher sensitivity to the volume change more than the strange susceptibility ($\chi^{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, $\chi^{2}_{BQ}$, $\chi^{2}_{BS}$ and $\chi^{2}_{SQ}$ for several volume selections and different $\mu_{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 ($\chi_{BQ}$) show a higher sensitivity to the volume change.

Also, the temperature dependence of the higher order baryon susceptibilities $\chi^{n}_{B}$ ($n$ = 4, 6 and 8) for different volume selections at $\mu_{B}$ values, $\mu_{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 $\mu_{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 $\mu_{B}$. Our calculations indicate that the conserved-quantities ($\chi^{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 $\mu_{B}$ and $T$ with decreasing system volume.

Acknowledgements

This work is supported by the US Department of Energy under contract DE-FG02-94ER40865.