# Hadron-quark phase transition in a hadronic and Polyakov–Nambu–Jona-Lasinio models perspective

###### Abstract

In this work we study the hadron-quark phase transition matching relativistic hadrodynamical mean-field models (in the hadronic phase) with the more updated versions of the Polyakov-Nambu-Jona-Lasinio models (on the quark side). Systematic comparisons are performed showing that the predicted hadronic phases of the matching named as RMF-PNJL, are larger than the confined phase obtained exclusively by the Polyakov quark models. This important result is due to the effect of the nuclear force that causes more resistance of hadronic matter to isothermal compressions. For sake of comparison, we also obtain the matchings of the hadronic models with the MIT bag model, named as RMF-MIT, showing that it presents always larger hadron regions, while shows smaller mixed phases than that obtained from the RMF-PNJL ones. Thus, studies of the confinement transition in nuclear matter, done only with quark models, still need nuclear degrees of freedom to be more reliable in the whole phase diagram.

###### pacs:

12.38.Mh,25.75.Nq## I Introduction

The study of the properties of the strongly interacting matter is experimentally supported by the heavy-ion collisions at ultrarelativistic energies, accomplished in the most sophisticated accelerators such as the relativistic heavy-ion collider (RHIC) and the large hadron collider (LHC). The measurements coming from these experiments are indirectly used to furnish the basic information on the phase diagram concerning the state of matter, being the different regions related to the distinct hadronic (confined quarks), and the quark-gluon plasma (free quarks) phases. In order to cover the temperature-density/chemical potential plane of the hadron-quark phase diagram, the new experiments should be able to compress the baryonic matter even to higher densities compared to the well known nuclear saturation density, fm. The new facility for antiproton and ion research (FAIR) fair () at Darmstadt, and the nuclotron-based ion collider facility (NICA) nica () at the Joint Institute for Nuclear Research (JINR) in Dubna, will make possible such extreme conditions through collisions foreseen to reach an energy range of to GeV fair (). Therefore, it will be possible to test the predictions from the theoretical calculations about the order of hadron-quark phase transition (crossover, first order or both), and also the regime of high density and moderate temperatures.

On the theoretical side, the strongly interacting matter is described by Quantum Chromodynamics (QCD). In its nonperturbative regime of large distances, or equivalently low energies, the most important method to study the QCD structure is the numerical lattice calculations lattice (), with the Monte Carlo simulations mcarlo () being its mainly representative. Such techniques provide results as for the pure gluon sector, i.e., in the limit of infinitely heavy quarks, as for systems including dynamical quarks. The difficulties in the latter are in the fact that at finite quark chemical potential () regime, the numerical calculations face to the fermion sign problem sign (). However, alternatives to solve this question are addressed in many methods expansao (); reweighting (); review (); imaginary (); dos (). On the other hand, one can deal with the high energy QCD regime (asymptotically free region) freeqcd () in a complete theoretical way through the perturbative treatment.

To investigate the moderate and high density regime of the transition, it is needed effective models such as the MIT bag model mit () and the Nambu-Jona-Lasinio (NJL) njl (); buballa (); outros () one, that curiously were firstly proposed to describe the hadronic mass spectrum. The former treats gluons and massless quarks as free particles in which the confinement phenomenon is incorporated by including a bag constant in an ad hoc fashion. The latter presents further similarities with the full QCD theory but do not take into account the confinement, since the quarks interact each other via pointlike interactions with no mediator gluons. One of the mainly aims of these kind of models is to provide the QCD phase diagram in the plane in order to distinguish hadronic phase where the chiral symmetry is broken, and the quark one. The quark stars description with these effective models is in the same important way buballa (); stocker1 ().

The phenomenology of the hadron-quark phase transition is often studied with the mentioned effective QCD models astroj (); wang () and other. Recently, a study about the QCD phase diagram was accomplished in Ref. kapusta () where the linear model was used to construct the phase transition curve, in the two-flavor quark system, as a function of the vacuum pion mass. It was shown that the point in the plane where is changed the phase transition order, the critical end point, is sensitive to this parameter. In such cases, where there is only one model (the quark one), the identified phases in the transition curve are associated to those with broken and restored chiral symmetry.

To distinguish hadron and quark phases we have the following well accepted phenomenology. The strongly interacting matter at very high temperatures and baryon densities should have quarks and gluons as the degrees of freedom. Motivated QCD models, cited above mit (); njl (); buballa (); outros () take these degrees of freedom into account. On the other side, at low temperatures ( MeV) and moderate densities of the order of , the hadronic phase can not be described by quarks and gluons. Instead, nucleons and mesons are the relevant degrees of freedom. Here, different approaches as Brueckner-Hartree-Fock by using realistic two-nucleon interactions bhf (), or relativistic mean-field (RMF) models may describe quite well the nuclear matter and finite nuclei properties. Between these two phases a mixed region takes place.

In this work, we focus in two structurally distinct models to treat the hadron-quark phase transition, namely, the Polyakov Nambu-Jona-Lasinio (PNJL) models in the characterization of the quark phase, and the relativistic mean-field nonlinear Walecka ones in the description of the phase where the chiral symmetry is broken, i.e., the hadronic phase. The former is a generalization of the NJL model in the sense that the confinement phenomenon is incorporated in its structure. For sake of comparison, we also perform a hadron-quark phase transition using the MIT Bag model in the quark side, in order to see explicitly the effect of the dynamical confinement contained in the PNJL model compared with the MIT Bag one.

The motivation for such an approach matching two models with different degrees of freedom is the following. MIT bag and PNJL models do not describe any of the very well known properties of the ordinary nuclear matter while attempting only to achieve, through phase diagrams, the boundaries where this matter starts to be confined. In this aspect only with such models a more complete description is poor. On the other side, the RMF models constructed to well describe the bulk nuclear matter and finite nuclei properties, are often applied to investigate high density regimes as it is common for neutron stars. Whether in such regions the RMF models are still valid is questionable. The phase transition between the quark and the hadron phase helps understanding when one needs to start describing nuclear matter with other degrees of freedom than baryons and mesons employed by RMF models. Therefore, our study brings more information for the mixed phase between these two still disconnected quark and hadronic phases. The use of different PNJL and RMF models is needed if one wants to have a more precise conclusion about how the different families of hadronic and quark models predict mixed phases. If they are similar or not, or still whether the description of the confined/deconfined phase transition only via the PNJL model itself is totally trustable when compared with the approach matching the two models. As we will show, this phase transition obtained by only the PNJL model, is expected to occur in a much smaller chemical potential for the same temperature differing strongly of the transition obtained by RMF-PNJL models. Future experiments, discussed already before, will show what description is correct.

After the first version of this manuscript was submitted, we became aware of a very similar study in which the hadron-quark phase transition was investigated only with one parametrization for the PNJL and RMF model gyshao (). Our systematic study here goes beyond that work, in the following aspects. Our calculations will take into account here different parametrizations for both, the PNJL and RMF models. In the latter, we will analyze a class of hard and soft equations of state, in a very large range of the nuclear matter incompressibility, MeV, and show that, actually, the RMF models affect quite small the entire hadron-quark phase transition. In the former, the most up to date Polyakov potential parametrizations will be used and compared each other. Still, a comparative study involving the hadron-quark phase transition with the MIT bag model will be performed.

This paper is organized as follows. In Sec. II, we present the quark models used in this work, starting by the MIT Bag model and later with the PNJL one, where a comparison between the current parametrizations, regarding their thermodynamics, is also done. In Sec. III, the mainly features of the RMF models are presented, and in Sec. IV, the phase diagrams concerning the hadron-quark phase transitions are shown. Our summary and conclusions are in the last section.

## Ii Quark models

### ii.1 MIT Bag model

Possibly, the simplest model to describe an approximate physics describing the matter where quark and gluons are the proper degrees of freedom of the system is the MIT bag model mit (). Such a scenario is supposed to exist at very high density or temperature regimes. Back in the big bang theory, at the very beginning, the universe was very hot and dense before hadronization. Nowadays, very high energy heavy-ion collision experiments in LHC try to recreate the signatures to confirm this hypothesis. In the case it is confirmed, hot hadronic matter undergoes a phase transition to the quark-gluon plasma (QGP) at uncertain but very high and very low temperature and density, respectively. The phase diagram connecting the QGP phase to a pure hadronic phase from high to zero temperature is still a big challenge and has a model dependence exhibited in previous studies when some relativistic hadronic models and the MIT bag ones were investigated jphysg (); debora (). The MIT bag model itself does not describe the confined hadronic phase but rather the QGP. Therefore, the hadronic phase has to be represented by any hadronic model in which the degree of freedom are baryons and mesons as we will present later on. Once each phase is modeled, the phase diagram is obtained thermodynamically by the Gibbs criteria, matching the chemical potential and the pressure at a given temperature (), across the phase boundaries when a phase transition takes place. Explicitly, these criteria are

(1) | |||||

(2) |

As previously presented jphysg (), the set of critical ’s establishes curves in the plane. Below the curve the system may be interpreted as nuclear matter described by the hadronic models (hadronic sector). The region in between and corresponds to a mixed (H-QGP) coexistence phase. Above the curve, the system is in a pure QGP phase.

For the quark-gluon plasma phase in the MIT bag model, the pressure and baryon number density are given by heinz (); boqiang ()

(3) | |||||

(4) |

where is the bag constant, is the number of flavors and is the QCD running coupling constant, depending on the quark-gluon plasma temperature and the quark chemical potential through the first order perturbative expression

(5) |

In this paper, the MIT bag model will be used only as a comparison with the PNJL model, regarding the predictions for the boundaries of the hadronic and QGP phases for each specific RMF model describing the hadronic phase.

### ii.2 PNJL model

#### ii.2.1 Confinement

The confinement can be measured by the Polyakov loop, defined as a Wilson line in the Euclidean space-time by exp where is the gauge field. By this definition, it can be viewed as a closed path connecting the same point in the space, at the different times and . The traced Polyakov loop, with being the quark color number, plays an important role concerning the deconfinement since it is the order parameter of this transition. Actually, the deconfinement is a consequence of the global center symmetry breaking. In pure gauge systems at finite temperature, where the periodic boundary conditions are satisfied, this symmetry is realized by the gauge transformations due to the way how they affect the boson fields. The consequence in the traced Polyakov loop is that , with what means that the condition has to be fulfilled in order to keep the center symmetry preserved. Since is related with free energy of a color source, through , if the center symmetry is realized, then , meaning that the quark is confined. In dynamical quark systems where the fermionic field explicitly breaks the center symmetry, generating automatically , the traced Polyakov loop is considered as an approximated deconfinement order parameter. The situation is analog to the cases in which the system Lagrangian density presents a mass term that explicitly breaks the chiral symmetry. The condensate is also an approximated order parameter, in this case for the restored chiral symmetry phase. The confinement effect, associated with , was implemented originally in the NJL model by Fukushima fuku1 () and deeply studied in Refs. weise1 (); weise2 (); weise3 (); weise4 (); weise5 (); weise6 (); weise7 (). The PNJL model is able to describe as the broken-unbroken chiral symmetry as the confinement-deconfinement phase transitions even being this situations, respectively, in the opposite regimes of vanishing and infinite quark masses fuku1 ().

#### ii.2.2 Thermodynamical calculations at

The connection between the fermion and the gauge field in the PNJL model is done by making the substitution in the Lagrangian density. The after bosonization of the system, and the mean-field approximation lead to the following grand canonical potential per volume weise2 ():

(6) | |||||

being , in the isospin symmetric system, the effective Polyakov loop potential in terms of and its conjugate , that we will discuss later, and the degeneracy factor due to the spin, flavor, and color numbers ( and . The constituent quark mass should obey the autoconsistent equation

(7) |

where the quark condensate , determined by the condition , is given by

(8) |

with

(9) |

(10) |

being the generalized Fermi-Dirac distributions. As pointed out in Refs. ratti (); costa (), an important consequence of the structure of the coupling between and the quark sector, is the possibility to deal with the PNJL model in the same theoretical way as in the NJL one, regarding the statistical treatment. The modification in the equations, as, for example, in Eq. (8), is in the use of these new distribution functions for quarks and antiquarks.

Through , all the thermodynamical quantities can be obtained, namely, the pressure, given by that leads to

(11) | |||||

the quark density, ,

(12) |

and the energy density, ,

(13) | |||||

Notice here that we have subtracted the vacuum value of , resulting in the addition of the term in the expressions (11) and (13).

The entropy density is obtained through , or by the thermodynamical relationship . Therefore we have

(14) | |||||

It should be mention at this point that the prescription of generalization of the NJL model to the PNJL one concerning the Fermi-Dirac distributions functions of particles, , and antiparticles, , to those shown in Eqs. (9) and (10), is valid for the entropy density functional given in Eq. (14), and not for its usual expression .

To deal with the entire PNJL model, characterized by the above equations of state (EOS), it is needed specify the Polyakov loop potential, . Different versions were proposed in the literature, and following the language of Ref. fuku2 (), we refer two of them by RTW05 weise1 (), and RRW06 weise2 (); weise3 (); weise4 (); weise5 (); weise6 (); weise7 (); ratti (); costa (). We call the other two used in our work by FUKU08 fuku2 (), and DS10 schramm (). Their functional forms are given, respectively, by,

(15) | |||||

(16) | |||||

(17) | |||||

(18) |

where

(19) |

In a general way, the Polyakov potentials are constructed in order to reproduce the well established data from lattice calculations of the pure gauge sector (where ), concerning the temperature dependence of the traced Polyakov loop and its first order phase transition, characterized by the jump of from the vanishing to a finite value at MeV (see the dotted curve of Fig. 2 in Ref. weise2 (), for instance). Actually, in this work we reduced this value to MeV, following the rescaling adopted in Ref. weise1 (), in order to better reproduce the lattice QCD results that we will show in Fig. 2, and the transition temperature at vanishing chemical potential, , also predicted in the lattice, and that we will discuss in Sec. IV. By taking this rescaling into account, we match more realistic PNJL models with the RMF hadronic ones, that already describe quite realistically nuclear matter and finite nuclei.

The FUKU08 potential presents two free parameters, MeV and , and was derived from strong coupling lattice expansion. Also to correctly reproduce the lattice results for we changed here the parameter value to .

The RTW05 parametrization was based on a Ginzburg-Landau ansatz, presenting the polynomial form in terms of the order parameters, , and . The improved version RRW06 uses, instead the polynomial terms of third and fourth order, the logarithm of the Jacobi determinant. Differently from the FUKU08 potential, the dimensionless parameters , and of the RTW05, and RRW06 ones (see their values in Table 1) are found with the additional fit to the lattice data for the energy density, entropy density, and pressure of the gauge sector as a function of the temperature, including the proper Stefan-Boltzmann (SB) limit at high temperature regime. Fig. 1 of Refs. weise1 () and weise2 (), for example, shows the good agreement among these quantities and the lattice results.

Potentials | ||||||
---|---|---|---|---|---|---|

RTW05 | ||||||

RRW06 | - | - | ||||

DS10 | - | - |

The DS10 potential, proposed in Ref. schramm (), presents baryon chemical potential dependence and is used in a hybrid chiral model that has both hadrons and quarks as degrees of freedom. This parametrization is able to reproduce the gauge lattice results, just like in the potentials shown before, and also allows that the hadron-quark phase transition occurs, at zero temperature, at density of 4 times the saturation density, and that the critical end point (CEP) is situated at MeV, and MeV schramm (). By using this particular parametrization, we remark that the quark density, Eq. (12), and the energy density, Eq. (13), have to be modified, respectively, to and .

To define completely the PNJL model it is needed to determine the coupling constant , the cutoff , and the current quark mass . This is done by imposing the reproduction of the certain vacuum values, namely, the pion decay constant fixed in , the quark condensate , and the pion mass . These values, together with the zero temperature expressions of the PNJL model, are shown in Appendix A.

#### ii.2.3 Comparing the Different Polyakov Potentials

To evaluate the PNJL model one needs to solve simultaneously Eq. (7) and the minimization conditions for the thermodynamical potential relatively to and . Along all our study we will follow the lowest order approximation described in Refs. weise4 (); weise6 () that automatically leads to . This approach reduces the set of coupled equations to Eq. (7), and

(20) |

with

(21) |

coming from the condition .

To see how sensitive is the PNJL model regarding the different potentials, we present the behavior of some important quantities. First we show in Fig. 1 how the order parameters and vary as functions of the temperature for some fixed values of the quark chemical potential.

If we use only the PNJL model to construct the diagram, the maxima of , and at a given are used to define the transition temperature weise1 (); weise4 (); weise6 (), in this case a crossover transition. Indeed, the strength of this mixing is utilized to find the value of in the FUKU08 potential fuku2 (). Other similar criterium to find the transition point is to localize the maxima of the chiral, and Polyakov susceptibilities fuku1 (); weise6 (). Here we will also use the Gibbs criterium to find such transition points since we are dealing with two different models to treat the different phases.

In Fig. 2 we show, at vanishing chemical potential, the behavior of some important thermodynamical quantities. The pressure, energy density and entropy density as a function of the temperature, and scaled by the respective Stefan-Boltzmann values , , and , are displayed in Figs. 2a-2c. In units of , we furnish the called interaction measure given by . This quantity show us how the model deviates from the noninteracting massless quark system, since in such regime its value is vanishing. Finally, as a function of the energy density, we see the evolution of the ratio that is closely related with the sound velocity through . In the SB limit, one has .

As a reference the recent QCD lattice results with temporal extent lat1 (), and lat2 () for the flavor system, in calculations with the improved fermion actions asqtad, and p4, are also shown in Fig. 2. Rigorous comparisons among lattice data and quantities such as the difference and the moments of the pressure, can be found, for example, in Refs weise1 (); weise2 (); ghosh (), and weise2 (); weise4 (); weise5 (); weise6 (), respectively.

Notice here the good agreement of the PNJL parametrizations and the lattice data, even the latter being originated from the 3 flavor system. We highlight RTW05 and RRW06 as being the models that better agree with the data. These results can be explained by the rescaling procedure adopted in the value (from MeV to MeV).

As mentioned in Ref. fuku2 (), the FUKU08 model was not constructed to reproduce the Stefan-Boltzmann limit in high temperatures. This explains the deviation of the lattice data from temperatures higher than MeV even after our rescaling in the parameter of the model, from to .

## Iii Relativistic mean-field models

Different from the microscopic approach in which the nucleon-nucleon potentials are developed to reproduce the well established data of the few-nucleon physics, and where the basic informations about the many-body system are extracted, e.g., via Brueckner-Hartree-Fock methods, the Quantum Hadrodynamics (QHD), based on local Lagrangian densities, uses the nuclear matter bulk properties observables at zero temperature to adjust its free parameters and thereafter construct all the thermodynamics of the particular hadronic framework. One of the most used models coming from QHD is the nonlinear Walecka (or Boguta-Bodmer) boguta () model, that is given by the following renormalizable Lagrangian density

(22) | |||||

with , and the nucleon degree of freedom represented by the spinor . The responsible mesons for the attractive, and repulsive parts of the nuclear interaction are described in this formulation by the scalar (), and vector () neutral fields, respectively. The saturation in nuclear matter is understood in the QHD by the almost vanishing value of at the saturation density, with , and being the Lorentz meson potentials.. , , and are, respectively, the masses of the nucleon, scalar, and vector mesons.

In their model, Boguta and Bodmer boguta () considered the cubic and quartic self-interactions in the scalar field , in order to improve the original Walecka model walecka (). In this version, the models are able to control through the fitting of their coupling constants, the values of the saturation density, binding energy () as well as the incompressibility (), and the effective nucleon mass ().

Through the Dirac equation of the model, one identifies its nucleon effective mass as

(23) |

with , , , and . Let us remark here that this definition of effective mass, also called Dirac mass, is valid for the relativistic models. A deep discussion about other definitions and concepts of this physical quantity can be found in Ref. jaminon (). The scalar and vector densities are written, at finite temperature, as

(24) | |||||

(25) |

with for symmetric matter. The usual Fermi-Dirac distributions to particles and antiparticles are defined, respectively, by and , and the relation between the effective chemical potential, , and the baryon chemical potential is given by .

From the energy-momentum tensor obtained through Eq. (22), one can obtain the pressure of the system, that reads

(26) |

already written in terms of the mean-field approximation, as well as Eqs (23)-(25).

For the construction of the hadron-quark phase diagrams, we choose the following set of RMF parametrizations: Walecka (WAL) npb (), NLB serot (), NL2 reinhard (), NLSH lalazissis (), NLB1 reinhard (), NL3 lalazissis (), NLB2 reinhard (), NLC serot (), NL1 reinhard () and NLZ2 bender (). The WAL and NLZ2 models are chosen because their incompressibilities, MeV and MeV, respectively, indicate very hard and very soft EOS for infinite nuclear matter. The other models intermediate both. A full list of the nuclear matter saturation properties including , and is found in Table 2 of Appendix B.

Important features can be highlighted concerning these relativistic models. At moderate temperatures, MeV, a liquid-gas phase transition is predicted furn-serot (). In Figs. 3a-3c we show the typical van der Waals behavior of the used parametrizations in the pressure versus density curves.

The Maxwell construction can be used to find the coexistence points, and consequently the coexistence curve. It was shown in Ref. jbatista () that this curve, scaled by the critical parameters, is not universal in the liquid region due importance of the interactions in such phase. Another interesting feature of the RMF models occurs at high temperature regime ( MeV), where other kind of phase transition takes place. This scenario was studied by Theis et al theis () in the context of the linear Walecka model at . The authors have shown that the order of the phase transition depends on the values of the coupling constants used to fit , and . A signature of this transition is exhibited in the behavior of the nucleon effective mass as a function of the temperature, since an abrupt decreasing of characterizes a first order phase transition. In Figs. 3d-3f we present the high temperature regime of RMF models used in our work concerning the behavior of at zero baryon chemical potential. A recent study regarding the high temperature regime of RMF models at different number of nucleons and antinucleons () was done in Ref. npb (); ratios ().

## Iv Results and Discussions

Now we present our mainly results regarding the hadron-quark phase diagrams, constructed from the Gibbs criteria shown in Eqs (1)-(2), defining first our notation for the transitions curves we will exhibit. We call the matching between the RMF models and the PNJL ones as the “RMF-PNJL” transition curves. When the transition is done with the MIT Bag model on the quark side, the modification is straightforward in the sense that the curve is denoted by “RMF-MIT”. Thus, the Gibbs criteria applied, for example, in the NL3 model together with the RRW06 one, generate the curve called NL3-RRW06. If the PNJL model is replaced by the MIT Bag one in this case, we denote this transition by NL3-MIT.

### iv.1 RMF-PNJL transitions

Regarding the condition of equal pressures, we stress that the hadronic pressure, Eq. (26), was used to match the quark ones, Eq. (11), and Eq. (3) for the RMF-PNJL and RMF-MIT curves, respectively.

Let us remark here that for the RMF-PNJL transitions, we adopt the Gibbs criteria with the additional constraint that the RMF confined phase is actually matching the deconfined phase of the PNJL model. We only consider the solutions satisfying this condition. An example of this restriction is shown in Figs. 4a-4b for the NL3-RRW06, and NL3-FUKU08 matchings at MeV.

Notice that the FUKU08 and NL3 pressures cross each other at the temperature given by MeV, but since at this temperature the FUKU08 model is still in the confined phase (see the down panel of Fig. 4b), we did not consider this crossing for the NL3-FUKU08 phase diagram. Here we are following the same procedure used in Ref. weise1 () concerning the definition of the transition temperature. In the case of the adopted rescaling for (in RTW05, RRW06 and DS10) and for (in FUKU08), the almost perfect coincidence between the peaks of the temperatures derivatives of and is lost. To circumvent this problem, for each chemical potential, the transition temperature , as the average of the different temperatures related to the peaks of and . For the FUKU08 potential at MeV, this average is MeV.

For the RRW06 and NL3 models, its pressures cross at MeV, where the RRW06 model already predicts the deconfinement, since its transition temperature is given by MeV (see the up panel of Fig.. 4b). Therefore, the point in which MeV and MeV contributes to the NL3-RRW06 phase diagram.

In the next figure we present the hadron-quark phase diagrams for all the RMF-PNJL parametrizations used in this work.

In Fig. 5, each individual panel depicts the diagrams where the matching between the models is done by keeping fixed the indicated quark model, and changing the hadronic one. In the same figure, for comparison, we also give the transition line of the phase diagram for each PNJL model itself constructed, as aforementioned, from the average of the temperatures related to the peaks of and . The full (dotted) lines stand for the crossover (first order) transitions. It is important to stress here that there are different criteria to construct these PNJL transition curves, such as the choice between different peaks of plsm () and . This fact can favor the emergence of the quarkyonic phase, predicted in Ref. mclerran (). A detailed study about these different choices and its consequences in the hadron-quark phase diagram is under progress.

An important feature regarding the PNJL curves are the values of the transition temperatures at vanishing chemical potential, MeV, MeV, MeV, and MeV, respectively, for the RRW06, RTW05, FUKU08 and DS10 parametrizations. Notice that these values agree very well with the lattice QCD result for this quantity given by MeV star (), or even with the more stringent value of MeV stringent (). This nice agreement is completely destroyed if the original values of and are used in the construction of the PNJL curves. In this case, the values of are higher than MeV for all the Polyakov potentials.

We also remark that the PNJL phase diagrams could also furnish other phases beside the confined/broken chiral symmetry, and deconfined/restored chiral symmetry ones, delimited by the full and dotted lines in Fig. 5. These specific phases are closely related to the instabilities of the ground state, present in any fermionic system at sufficiently low temperatures, that are overcame by the formation of the fermion pairs (in BCS theory, these are the so-called Cooper pairs). In particular, the quark system described by the PNJL models should also be affected in this regime, giving rise to the emergence of the two-flavor superconducting color (2SC), and the color flavor locked (CFL) phases. This phenomenology can be incorporated in the PNJL model, via inclusion of the diquark condensate term in its Lagrangian density. However, in this work we focus exclusively in the nonsuperconducting phases of the hadron-quark diagrams by treating only the simplest version of the PNJL model, in which the color condensates are not being taken into account. Very interesting studies concerning the consequences of the 2SC and CFL phases in the description of quark stars were performed for the NJL sc1 () model, in which it is shown the role played by the diquark coupling strengths. Recently, this effect of the color superconductivity was also analyzed for the PNJL sc2 () model, where a modified version of the DS10 Polyakov potential is used in the description of the quark core of massive hybrid stars.

Regarding the RMF-PNJL transitions themselves, we stress the following interesting results. First of all, notice that for all the RMF-PNJL matchings, one can delimit very narrow bands that encompass all the transition curves, being the WAL-PNJL and NLB1-PNJL matchings their extremes. This is not an obvious result since we are dealing with a large set of RMF models, in the sense that they present a range of very soft to very hard equations of state, i. e., from the NLZ2 model with MeV, up to the WAL one in which MeV. This result indicates that at least for RMF models with scalar self-interactions of third and fourth order, all of them will predict very similar behavior concerning the hadron-quark phase transition, in connection with PNJL models.

Our systematic study still shows that the additional constraint exemplified in Figs. 4a and 4b, are useful to define the CEP’s of the RMF-PNJL transitions, given approximately by ( MeV, MeV), ( MeV, MeV), ( MeV, MeV), and ( MeV, MeV), respectively, for the RMF-RRW06, RMF-RTW05, RMF-FUKU08, and RMF-DS10 transitions. In a different way, the CEP is also defined in Ref. gyshao (). Notice how this restriction avoids the hadron-quark phase transition curves go inside the confined phases predicted by the PNJL transitions.

This situation is completely modified if we use the original values of and of the PNJL models. In this case there is no CEP for the RMF-PNJL transitions, in the sense that there is no transition inside the hadron phase of the PNJL models themselves. The only exception is for the RMF-FUKU08 transitions, that present the same qualitative behavior shown in Fig. 5 for the calculations done with the original value .

We still highlight here a very important physical consequence of the construction of the hadron-quark phase transitions via connection of the RMF and PNJL models. The panels in Fig. 5 clearly show us that the hadronic degrees of freedom present in the RMF models, make the system much more resistant to the quark liberation, in function of the density/chemical potential increasing, when compared with a system only described by the PNJL model. The physical origin of this important difference is due to the very known repulsive nature of the nuclear force, represented in the RMF models by the vector interaction, which strength is controlled by the coupling constant , see Eq. (22) (in nuclear matter equations of state, this strength is actually controlled by the ratio ). The repulsion makes the system support more strongly isothermal compressions. In other words, the RMF-PNJL transitions predict confined (hadron) phases larger than those obtained exclusively with the PNJL models. This region of highly compressed matter of the hadron-quark phase diagram, is expected to be reached in the new experiments, such as the planned to occur in the new facilities FAIR/GSI fair (), and NICA/JINR nica (). Therefore, it will be explicit what is the magnitude of the role played by the repulsive interaction part of the nuclear force, even guiding possible selections of better parametrizations.

The difference between the RMF-PNJL and the PNJL transitions can be decreased if we use PNJL models that also contain repulsive interactions, i.e., that present vector fields in its structure. Notice that the PNJL models used here are based on a structure that present only attractive interactions. There is no explicit terms proportional to the quark density in Eq. (6), coming from vector-type fields. Moreover, as pointed out in Ref. fuku2 (), there is no constraint at all for the choice of the strength of this kind of interaction in the PNJL model. A study about the determination of this magnitude, by making minimum the difference shown in Fig. 5, is underway.

### iv.2 Comparison with RMF-MIT transitions

We perform now, systematic comparisons between the RMF-PNJL transitions with RMF-MIT ones, in order to verify explicitly the role played by the dynamical confinement of the PNJL models in the hadron-quark phase diagrams. To construct such curves, we use the particular value of MeV, that furnish a critical temperature at vanishing chemical potential consistent with the lattice simulation results.

Firstly, we show in Fig. 6 the behavior of RMF-MIT curves in the plane for all the hadronic parametrizations.

Some interesting points about these particular transitions have to be mentioned here. Notice that differently from the RMF-PNJL diagrams, all the RMF-MIT curves start at the same temperature, MeV. There is no critical end points as shown in Fig. 5. Similar behavior between the curves can be seen at , since all used models are lied in a narrow band in this region. They lie in the range around MeV MeV ( MeV MeV) for the RMF-PNJL (RMF-MIT) curves. For the RMF-MIT transitions, this behavior is strongly changed if the RMF models with higher order terms in the vector field, or even mixing terms between and are used on the hadronic side. This is the case, for example, for the models used in Ref. jphysg (). We stress that the behavior shown in Figs. 5 and 6 is characteristic of the transitions in which the used RMF parametrizations only contain cubic and quartic self-coupling terms in the scalar field , see Eq. (22).

This almost model independent result for the hadron-quark transition at for the RMF-PNJL/MIT matchings, see Figs. 5 and 6, is not surprising from the point of view of the PNJL models, since the Polyakov potentials related to the RRW06, RTW05, and FUKU08 models vanish in the zero temperature regime, see Eqs. (15)-(17), and the modified Fermi-Dirac distributions, Eqs. (9) and (10), behave as the conventional ones in the limit. In particular, the DS10 Polyakov potential vanish at also due to , see Eq.. (18). That is, all the PNJL models used here are converted in the same model at , in this case, the NJL one on the quark side. This is also the case in the RMF-MIT transitions, i.e., one has only one parametrization of the quark model since we fixed and in Eq. (3) to construct the diagrams shown in Fig. 6. We remark that by construction, at all the RMF models lead to different bulk nuclear matter properties. We stress here that this almost model independence in the value at can be relevant to the study of hybrid stars, composed by a quark core and a hadron crust malheiro1 (); malheiro2 ().

In order to see the effect of the dynamical confinement present in the PNJL models in the RMF-PNJL transition curves, we explicitly compare the phase diagrams of both, RMF-PNJL/RMF-MIT phase diagrams in the same Fig. 7.

In this figure one can clearly see that the hadronic region predicted for the RMF-PNJL matchings is always smaller than that obtained from the RMF-MIT ones, i.e., the quark degrees of freedom emerge before in the RMF-PNJL curves. Although we have presented only the transition curves shown in Fig. 7, we streamline that all the other diagrams follow the same pattern, thus the respective RMF-PNJL/RMF-MIT curves can be considered as representatives of all the models treated in this work. An enlarged view of these features can be viewed in Fig. 8.

In this figure, despite present the larger hadron regions, that are delimited by the left branches of the curves, the RMF-MIT diagrams show that the mixed phase containing both hadrons and quarks is actually smaller in such diagrams. This difference, specifically at is extremely important to the study of hybrid stars, since its quark core is directly affected by the maximum density of the mixed phase malheiro1 (); malheiro2 ().

As our final remark, we point out that in the construction of the RMF-PNJL diagrams, we proceed in a different way that used in Ref. ciminale (). In this reference, the authors used a bag constant also in the PNJL pressure in order to ensure that the quark pressure is less than the hadronic one in the confined phase. Here we force , by subtracting from the grand-potential its vacuum value.

## V Summary and Conclusions

In this work we have studied the hadron-quark phase diagrams using two different models in the description of the distinct phases. The quark degrees of freedom were described via the recently suggested PNJL model fuku1 () which incorporates the confinement phenomenon in the previous NJL one. As one of our results, we compared the current PNJL parametrizations that differ each other by the Polyakov potential . We used the polynomial form RTW05 weise1 (), the logarithmic RRW06 weise2 (); weise3 (); weise4 (); weise5 (); weise6 (); weise7 (); ratti (); costa (), the FUKU08 fuku2 (), and also included the DS10 schramm () one that presents a chemical potential dependence. The Fig. 2 show the good agreement with the lattice data, specially for the RRW06 and RTW05 parametrizations even the PNJL models being treated in the two-flavor system. It was also shown that the parametrizations furnish similar results for the analyzed thermodynamical quantities. The good agreement with the lattice data remains valid even for the transition temperature at vanishing chemical potential, predicted to be given by MeV stringent (). Our calculations give the values of MeV, MeV, MeV, and MeV, respectively, for the RRW06, RTW05, FUKU08 and DS10 PNJL models used in this work. We still remark that these nice values are obtained with the rescaling of the parameters (RRW06, RTW05 and DS10) and (FUKU08).

On the hadronic side, we have used the well known RMF nonlinear Walecka models in its version that contains cubic and quartic self-coupling in the scalar field. We chose to deal with a set of these models that encompasses several incompressibilities values, representing very hard and soft equations of state.

Regarding the phase diagrams, we have constructed the RMF-PNJL curves using the Gibbs criteria and the additional constraint that the RMF pressure matches the PNJL one only in its deconfined phase. This condition is exemplified in Figs. 4a-4b, and played an important role in the final diagrams, since it determines the critical end points given approximately by ( MeV, MeV), ( MeV, MeV), ( MeV, MeV), and ( MeV, MeV), respectively, for the RMF-RRW06, RMF-RTW05, RMF-FUKU08, and RMF-DS10 transitions. Moreover, it is also important to stress that all these transitions furnish very similar results in a such way that one can define a very narrow band in the plane encompassing all the RMF-PNJL curves, being the WAL-PNJL and NLB1-PNJL the limiting curves of these bands, see Fig. 5. This is a surprising result since we are dealing with a very large class of RMF models. In principle, there is no reason to the hadron-quark phase diagram behaves in a very similar way with such variety of RMF parametrizations.

Other important result concerning the RMF-PNJL phase diagrams is their different predictions compared to those obtained exclusively with the PNJL quark models. We found that the hadron phase described by the RMF-PNJL transitions is meaningfully larger than that predicted by the PNJL ones, i.e., the RMF-PNJL hadron phase is more resistant to the isothermal compressions. Physically, such difference is due to the repulsive part of the nuclear force described in the RMF models by the vector field interaction. Therefore, one become clear that there are very different results in treating the hadron-quark phase transition via two distinct models, taking into account the different degrees of freedoms (hadrons and quarks), or only via quark models, even being the latter very realistic ones as the PNJL models that nicely agree with lattice QCD data, and where the dynamical confinement is considered in the NJL model through the Polyakov loop. We also remark that the region where the different descriptions do not agree each other will can be accessed in the future experiments planned to occur in the new facilities FAIR/GSI fair (), and NICA/JINR nica ().

As a final result, we compared the RMF-PNJL curves with that constructed by the matching between the hadronic models and the MIT Bag one, that incorporates the confinement via inclusion of the bag constant . The Figs. 7 and 8 show that the dynamical confinement predicted in the PNJL model force the hadronic phase of the RMF-PNJL diagrams be smaller than the RMF-MIT ones. Curiously, the opposite occurs for the mixed phase, where hadrons and quarks coexist, see Fig. 8. These comparisons were done by assuming the value of MeV for the MIT bag model, that nicely predicts a value of MeV for the transition temperature at vanishing chemical potential, see Fig. 6.

###### Acknowledgements.

This work was supported by the Brazilian agencies FAPERJ, CAPES, CNPq, and FAPESP. The authors gratefully acknowledge professor Tobias Frederico for stimulating discussions and valuable comments.## Appendix A Zero temperature expressions for the PNJL model: NJL sector

In this appendix we show the zero temperature expressions of the PNJL model, that actually are exactly the NJL ones, and give the remaining parameters needed to define the PNJL models presented in our work. We refer here only to that parametrizations that do not present contributions containing the traced Polyakov loop in the regime, i.e., we only consider the cases in which .

In the regime, the energy density given in Eq. (13), and the scalar density in Eq. (8), are replaced by

(27) | |||||

(28) |

with

(29) | |||||

(30) |

where is the Fermi momentum of the quark. The quark density is , and the pressure reads

(31) |

with the quark chemical potential given by .

Thus, the respective vacuum expressions obtained at are

(32) | |||||

(33) |

So, fixing the values MeV, MeV, and MeV, and using Eq. (33) together with the Gell-Mann-Oakes-Renner relation , and

(34) |

one obtains MeV, MeV, and GeV. The constituent vacuum quark mass obtained from these values is MeV.

The assumption leads to the following final expression for the energy density,

(35) |

where . Notice that the same condition also ensures that . For the parameters aforementioned one has that MeV. In this work we have considered this value of in the Polyakov potentials, even for the DS10 parametrization.

## Appendix B Parametrizations of the RMF models

Some important saturation quantities of the RMF models used in our work are listed in the next table.

Model | (fm) | (MeV) | (MeV) | |
---|---|---|---|---|

Walecka (WAL) | ||||

NLB | ||||

NL2 | ||||

NLSH | ||||

NLB1 | ||||

NL3 | ||||

NLB2 | ||||

NLC | ||||

NL1 | ||||

NLZ2 |

The binding energy is calculated from the energy density at ,