Equation of state for hot QCD and compact stars from a mean field approach

The thermodynamic properties of high temperature and high density QCD-matter are explored within the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF. The quark sector of the CMF model is tuned to describe the $\mu_B=0$ thermodynamics data of lattice QCD. The resulting lines of constant physical variables as well as the baryon number susceptibilities are studied in some detail in the temperature/chemical potential plane. The CMF model predicts three consecutive transitions, the nuclear first-order liquid-vapor phase transition, chiral symmetry restoration, and the cross-over transition to a quark-dominated phase. All three phenomena are cross-over, for most of the $T-\mu_B$-plane. The deviations from the free ideal hadron gas baseline at $\mu_B=0$ and $T\approx 100-200$ MeV can be attributed to remnants of the liquid-vapor first order phase transition in nuclear matter. The chiral crossing transition determines the baryon fluctuations at much higher $\mu_B\approx1.5$ GeV, and at even higher baryon densities $\mu_B\approx2.4$ GeV, the behavior of fluctuations is controlled by the deconfinement cross-over. The CMF model also describe well the static properties of high $\mu_B$ neutron stars as well as the new neutron star merger observations. The effective EoS presented here describes simultaneously lattice QCD results at $\mu_B=0$, as well as observed physical phenomena (nuclear matter and neutron star matter) at $T\cong0$ and high densities, $\mu_B>1$ GeV.


I. INTRODUCTION
At low temperatures and high baryonic densities a transition from hadronic matter to a cold deconfined state with quark degrees of freedom is expected [1], that is supported by perturbative QCD calculations [2]. This transition is a particular case of appearance of a deconfined state of quark-gluon plasma (QGP) that expected at large densities and large temperatures [3]. The study of the unknown QCD phase diagram is one of the main motivations for the state-of-the-art research in the nuclear and particle physics community. A rich phase structure is conjectured for finite temperatures and chemical potentials [4][5][6][7][8]. The QCD phase structure is a most important ingredient for the understanding of the Early Universe, ultra relativistic heavy ion collisions, and of the evolution, structure and inspirals of neutron stars [9]. Even though QCD is a well-established theory with only a few parameters, perturbative calculations are inappropriate in the crossover regions of T and µ B discussed here due to the large values of the QCD coupling constant at scales relevant for most of these applications [10]. The infamous sign problem [11] prohibits lattice QCD (LQCD) calculations at finite densities.
The running experiments at the Large Hadron Collider (LHC), at the Relativistic Heavy Ion Collider (RHIC), at the Super Proton Synchrotron (SPS), and at the Heavy Ion Synchrotron (SIS) provide state-of-the-art data of the measurements of properties of matter produced in heavyion collisions (HIC). The detailed information of the particle production as measured in these experiments allow † Deceased.
to extract both thermodynamic and kinetic characteristics of the system that is created.
Astrophysical observations of compact stars, along with data from the recent gravitational-wave detection by LIGO provide an additional tool to probe the equation of state of dense nuclear and possible quark matter [12][13][14][15][16][17][18] in the region of moderate temperatures and high baryon densities, close to those as created in HIC.
First principle LQCD calculations suggest a smooth crossover transition at vanishing baryochemical potential µ B = 0 from hadronic to partonic degrees of freedom [19]. Although there is no indication of a first or a second order phase transition in the energy density, pressure, entropy density and speed of sound, there are other observable in LQCD -chiral susceptibilities which seem to indicate a chiral cross-over at a pseudocritical temperature T pc ≈ 155 MeV [20,21]. However, the extension of LQCD calculations to finite µ B is a difficult problem. There are approximate lattice methods, such as extrapolations by Taylor expansion and the analytic continuation from purely imaginary to real µ B are reasonable only for small baryon densities. The exploration of higher baryon densities require effective QCD models, which respect the known symmetries of QCD and describe appropriately the known phenomenology of strong interactions.
The current knowledge of the properties of strongly interacting matter suggests a number of features that ought to be incorporated in any reasonable effective QCD model: • First, nuclear matter and the nuclear liquid-vapor phase transition at moderate temperatures, close to the nuclear saturation density [22]; • Second, the chiral symmetry restoration should lead to the Stefan-Boltzmann limit for the thermodynamic properties at high temperature and at high chemical potential [23]; • The transition from hadronic-to quark-gluon degrees of freedom, at high temperatures and/or chemical potentials, is a crucial ingredient for the consistent description of QCD matter.
Often, these different aspects of QCD are modeled within separate frameworks, which are then merged through various constructions. The present work formulates a single combined framework to describing the QCD thermodynamics, which simultaneously satisfies all the constraints, from lattice QCD, and known nuclear matter properties, as well as neutron star observations.
The resulting equation of state is then used to estimate various properties of systems created in both heavyion collisions and neutron star physics. Sec. II presents the description of the Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean field model, CMF. Section III describes the fine-tuning of model parameters to the µ B = 0 LQCD data on the trace anomaly, and presents a comparison with the CMF model predictions for various conserved charge number fluctuations with the corresponding lattice data. The QCD phase diagram deduced from the CMF model is studied in Sec. IV. The creation of hot and dense QCD matter as created in heavy-ion collisions at various collision energies is studied using the 1-dimensional hydrodynamics in Sec. V, the respective trajectories along the QCD phase diagram are explored as well. Section VI presents the CMF model predictions for the observed neutron star properties. The Chiral SU(3)-flavor parity-doublet Polyakov-loop quark-hadron mean-field model, CMF, is an extension of the previously proposed σ-ω model with parity doubling for nuclear and hadron matter [24][25][26][27][28][29]. The CMF model was extended to include quark degrees of freedom [30][31][32][33][34]. This model is a phenomenological effective unified approach to describe interacting hadron-quark matter. The Lagrangian includes essential symmetries and features of QCD. These include: • Chiral symmetry restoration in the hadronic sector, in particular the baryon parity doubling so an explicit mass term for baryons is possible, even when chiral symmetry is restored. This leads to a restoration of mass degeneracy among baryons and their respective parity partners [35,36].
• Eigenvolume corrections for hadrons, which allow for an effective modeling of their repulsive interac-tions. This suppresses hadronic densities and ensures a transition to a parton-dominated matter at large densities, when quark and gluon d.o.f. appear.
• Chiral symmetry restoration for quarks and a dynamical generation of their masses.
• The Polyakov loop via a QCD-motivated potential incorporates the deconfinement transition.
A detailed description of the hadronic part of the CMF model can be found in the literature [30,31,33,37]. It is based on a realization of a σ-ω model in mean-field description. Here the relevant fermionic degrees of freedom are baryons that interact through mesonic meanfields. The version of this model used here includes all states in the SU(3) f baryon octet, together with their parity partners, i.e. states with the same quantum numbers but opposite parity. The LQCD data suggests that the same mechanism should be implemented for higher baryonic states -the baryon octet [36], this is a plan for future studies. In the limit of chiral symmetry restoration these parity partner states should be degenerate and their masses equal, which then serves as a signal for chiral symmetry restoration. To allow for such a behavior, the baryon masses are dynamically generated by their couplings to the scalar σ-field and the scalar strange ζfield, which serve as the order parameters for the chiral transition: σi σ ± g (2) ζi ζ .
Here + stand for positive and − for negative parity states, g (j) i are the coupling constants of baryons to the two scalar fields, m 0 = 759 MeV is the baryon mass at the restored phase. In addition, there is an SU(3) f symmetrybreaking mass term proportional to the strangeness content of the baryons, there n s is the number of strange quarks in the baryon, and m s = 130 MeV is the mass of the strange quark. The couplings g (j) i are tuned to reproduce the vacuum masses of baryons.
The mean-field values of the chiral fields are driven by the thermal contribution from baryons and quarks, and controlled by the scalar meson interaction, driving the spontaneous breaking of the chiral symmetry: with I 2 = (σ 2 + ζ 2 ) , I 4 = −(σ 4 /2 + ζ 4 ), where V 0 is fixed by demanding that the potential vanishes in the vacuum.
The parameters of the scalar and the vector interactions are fitted to describe nuclear matter properties [33,38]. Contributions of all established hadronic resonances are included here with their vacuum masses [39]. These states can be coupled to meson fields as parity doublets as well. However, this is not done in the current implementation, they only interact with the other particles via their excluded volume only.
The quark degrees of freedom are incorporated similarly with PNJL approach [40]. The appearance of quarks is controlled by the value of Polyakov loop Φ, which plays the role of the order parameter for the deconfinement transition. The coupling of the quarks to the Polyakov loop is introduced through the thermal energy of the quarks. Their thermal contribution to the grand canonical potential Ω is given by and (5) The sums run over all light quark flavors (u, d, and s), d i is the corresponding degeneracy factor, E * i = m * 2 i + p 2 is the energy, and µ * is the chemical potential of the quark. Note, two-and three-quark contributions to Ω are omitted in the CMF model since hadronic excitations are explicitly included.
The effective masses of the light quarks are generated by the σ field (non-strange chiral condensate) as well, the mass of the strange quark is generated by the ζ field (strange quark-antiquark state). The small explicit mass terms δm q = 5 MeV, and δm s = 150 MeV for the strange quark, and m 0 q = 253 MeV which corresponds to an explicit mass term which does not originate from chiral symmetry breaking: Dynamics of the Polyakov-loop is controlled by the effective Polyakov-loop potential U (Φ, Φ * , T ) [41]: The parameters of this potential can be fixed to the lattice QCD data in the pure gauge sector [41]. However, this yields an unsatisfactory description of the (2+1)-flavor QCD thermodynamics when the hadrons are explicitly included in the model. Therefore the parameters of the Polyakov Loop potential are adjusted in the present work to describe properly the (2+1)-flavor lattice data.
The CMF model incorporates excluded-volume effects, in order to suppress the hadronic degrees of freedom in the regions of the phase diagram where physically quarks and gluons dominate [31]. Consequently, all the thermodynamics densities, including the quark contribution, are reduced as parts of the system are occupied by EVhadrons: the v j are the eigenvolume parameters for the different species. p is a system pressure without contribution of mean fields, µ * is the chemical potential of the hadron.
The v is assumed to be v B = 1 fm 3 for (anti-) baryons, v M = 1/8 fm 3 for mesons, and is set to zero v q = 0 for quarks. To introduce the constraints on the CMF model parameters from lattice QCD at high temperature and zero net-baryon density we use the QCD trace anomaly I, "the interaction measure", as a reference:

III. CONSTRAINING THE CMF MODEL TO THE LATTICE DATA
The free parameters of the present model are tuned to reproduce this quantity. The LQCD trace anomaly permits to calculate all observable thermodynamic quantities. Analysis of other Lattice data (chiral susceptibility) seem to show that for chirally related observables there is a crossover transition with a "pseudo"-critical temperature at T ≈ 156 MeV. The analysis of this data by a phenomenological model suggests a half-hadron, half-quark composition in that region [42].
The parameters of the CMF model's quark sector needed to reproduce the trace anomaly data from LQCD are found by a least mean squares fitting procedure for the parameters of the Polyakov loop potential U (Φ, Φ * , T ) and for the coupling constants g qσ and g sζ of the quarks to the chiral condensates σ and ζ, respectively. All in all this fixes 5 model parameters, T 0 , a 1 , a 2 , b 3 , g qσ = g sζ (we set g qσ and g sζ to the same value). The quark parameter fitting is performed through a scan over the parameter space on a 8 × 6 × 7 × 6 × 6 sized grid, minimizing the root mean square deviation of the CMF model data on I/T 4 and computed on the lattice results. The resulting parameter values are presented in Table I. The comparison of the CMF model to the lattice data is shown in Fig. 1.  The values in Table I suggest couplings g qσ , g sζ of quarks to the chiral fields of about 1/3 of the baryons, this one may expect from the additive quark model. Larger values of quark couplings would significantly influence the size of the peak in interaction measure I/T 4 , as studied in Ref. [33]. In case of large values of g qσ , g sζ the strong interplay between the chiral symmetry restoration and the deconfinement transition would result in too large values of the interaction measure. The large peaks in the baryon number susceptibilities are in contrast to the lattice data.
Higher order baryon number susceptibilities χ B n that are a LQCD measure of the particle number fluctuations as well as the curvatures of various lines of constant physical quantities are also interesting in scope of LQCD data. The behavior of the χ B 4 /χ B 2 at µ B = 0 is presented in this section, the study for finite values of µ B is presented in Sec. IV. The comparison of the CMF model with the available LQCD data for the χ B 4 /χ B 2 is shown in Fig. 2, indicating a fair agreement of the CMF model with the lattice data.
Lattice QCD studies often explore regions of finite µ B by using the Taylor series expansion. The Taylor expansion in series of T and µ B up to O(µ 4 B ) was used in [45] to calculate "lines of constant physics": i.e. lines in the T − µ B plane where certain thermodynamic quantities like pressure, energy density, and entropy density P, ε, s, are constant. The coefficients κ f 2 and κ f 4 (f ≡ P, ε, s) represent these contour lines in the T − µ B plane using the following parametrization [45]: Here the coefficients κ f 2 and κ f 4 are calculated from Eqs. (12) and (13), see Ref. [45] for details: Here The coefficients κ f 2 and κ f 4 are calculated in the CMF model for the pressure P , the energy density ε, and the entropy density s as functions of the temperature T . The CMF model predictions are in a reasonable agreement with recent LQCD calculations [45]. The rather low values of κ f 2 and κ f 4 suggest also small curvatures of lines of constant physical observables in the temperature region studied here. Effects of the finite chemical potential are small, therefore these lines are almost horizontal in the T -µ B plane. The coefficients for the entropy and the energy density indicate that κ s 2 < κ ε 2 , meaning a decrease of the entropy density along the lines of constant energy density.

IV. THE CMF MODEL PHASE DIAGRAM
Two order parameters, the chiral condensate σ and the Polyakov loop Φ, plus the interacting baryon octet within the SU(3)-flavor σ-model permit four different phases within the CMF model. These phases are characterized as: • A dilute gas of interacting hadrons; • A hadronic liquid -a dense hadronic phase, the transition from the hadron gas to the hadronic liquid is the nuclear liquid-vapor phase transition. Quarks start to appear in the hadronic liquid, but they are negligible; • A chirally restored phase, where the mass symmetry between the parity partners is restored. Here the quark masses are decreased, hence quarks give a sizable contribution to the thermodynamics; • A quark dominated phase, where most of the energy density is carried by quark and gluon degrees of freedom. The gluon contribution is modeled by the Polyakov loop potential [41] and q / tot ≈ 1.
The baryon number susceptibilities χ B n , that can be calculated using Eq. 10, are proportional to the respective cumulants of the baryon number distribution. Higherorder baryon number susceptibilities do increase in proportion to the increasing powers of the correlation length [47]. Such an increase in correlation length would be reflected in large values of 2nd and higher-order susceptibilities in the vicinity of a critical point and in the region of phase transition. Hence, these quantities are useful indicators of a critical behavior in the CMF model. Deviations of χ B n from the corresponding baselines indicate a transformation between different phases, which is reflected usually in a non-monotonic behavior of these observables, e.g. skewness χ B 3 /χ B 2 and kurtosis χ B 4 /χ B 2 . The calculated skewness (χ B 3 /χ B 2 ) and kurtosis (χ B 4 /χ B 2 ) (Fig. 4) in the CMF model exhibit non-trivial structures in the T -µ B phase diagram. The regions of deviations from baseline separate regions with quantitatively different properties which are often dubbed -"phases". Note the sharp phase boundaries indicate first order phase transitions, FOPT, these are only observable at quite moderate temperatures T < 50 MeV. The hadron phase located at both low temperatures T and baryon chemical potential µ B represents a dilute gas of interacting hadrons. There, the fluctuation measures are quite close to unity, consistent with the Skellam distribution baseline. The system exhibit a FOPT to a dense hadronic liquid phase with rising chemical potential µ B ≈ 1 GeV. Here, fluctuations are reduced due to the repulsive interactions. Quarks start to appear in moderation. The liquid phase exhibits additional FOPT at µ B ≈ 1. lam baseline. Hence, the fluctuation measures are rather small 1 before and after the "transition" .
The chiral critical point of the CMF model is located at a rather low temperature T CP chiral ≈ 17 MeV. This value is close to the critical temperature of the nuclear liquidgas transition in the same model, the critical µ CP B is remarkably different, though: the appearance of the parity partners controls the dynamics of the chiral fields: as the parity partners -in the CMF model -obey the same repulsive interaction strength as the nucleons, the critical point appears at that low temperature. This phenomenon has been observed in various mean field models before.
The different phases shown in Fig. 5 in the T −µ B plane are related to the chiral field σ and the quark fraction. The chiral field is close to its vacuum value, σ = σ 0 , at the hadron gas region, here the quark fraction is close to zero, as expected. Both observables deviate from their vacuum values at higher densities and temperatures only.
The chiral field drops off more slowly or µ B = 0 than seen in lattice QCD calculations, where chiral field (the sigma field) (upper) and the quark fraction 1 3 nq/nB (lower) of the CMF model in baryon chemical potential µB and temperature T plane. Note that the rather fast change of the chiral condensate appears at moderate energy densities, while deconfinement appears only at much higher energy densities/chemical potentials. rapidly drops around T = 160 MeV. The reason for this discrepancy is due to the fact that in the present CMF model the thermodynamics at these temperatures are strongly influenced by baryonic resonances which are not coupled to the chiral fields. Baryon resonances like the ∆ ought to be coupled to chiral fields including their chiral partners. This brings down the chiral condensate of lower temperature as can be seen in [48].
The speed of sound is another important signature for the simulations of the dynamics of heavy ion collisions and neutron star mergers. c 2 s/n B presents a derivative c 2 s/n B = ∂P ∂ε at constant S/A = s/n B entropy per baryon, that allows to estimate a speed of propagation of soundlike excitations in non-dissipative hydrodynamic evolution. The isentropic speed of sound can be calculated as [49]:  Table II for details.
The partial derivatives with respect to the chemical potential and to the temperature are performed at constant temperature and at constant chemical potential, respectively. The calculated speed of sound shows three local minima which correspond to three locally softest points of the EoS. These three minima correspond to phase boundaries, where the baryon number susceptibilities present a non-monotonic behavior. Note, that the speed of sound reaches quite large values, c 2 s ≈ 0.7, in the higher density region of nuclear matter. This high speed of sound results due to the strong repulsion between the baryons, before the onset of deconfinement. Thereafter, the vector repulsion and baryon excluded volume cease, as such  [50,51]. The entropy per baryon S/A, temperature T , the initial baryon density nB/n0 and the quark fraction 1 3 nq/nB are presented for various colliding energies.
terms have not been predicted for the quarks.

V. APPLICATION TO HEAVY-ION COLLISIONS
The presented EoS is used as input for hydrodynamical simulations of both heavy-ion collisions and neutron star mergers. To illustrate which regions of the phase diagram can be reached in collisions at low and moderate collision energies the stationary 1-dimensional Taub adiabat model is used [50][51][52]. The expansion is described at lines of constant entropy per baryon S/A = const (isentropes). These lines depict the isentropic matter evolution of ideal fluid dynamics at different collision energies.
The entropy is produced in the earliest stage of a heavy ion collision by the shock violent compression [53]. During the system's expansion there is only a moderate increase of entropy due to the rather small viscosity [54,55], hence, an isentropic expansion scenario is a reasonable approximation [56].
The expansion of the equilibrated matter then continues until the system becomes so dilute that the chemical as well as the kinetic freeze-out occur and the chemical composition is fixed.
The entropy per baryon (S/A) is calculated in the 1dimensional stationary scenario of central heavy ion collisions -the two colliding slabs of cold nuclear matter [52,53,[57][58][59][60][61] conserve the baryon number, energy and momentum across the shock front in accord with the relativistic Rankine-Hugoniot equation (Taub adiabat), RRHT, [50,51]. Thus, the produced entropy is directly associated to the collision energy. The thermodynamic properties across the shock front are described by the RRHT-equation where P 0 , ε 0 and n 0 correspond to the initial pressure, energy density, and baryon density in the local rest frame of each of the two slabs. The two symmetric slabs consist of the nuclear matter in the ground state, P 0 = 0, ε 0 /n 0 −  Table II for  details. m N = −16 MeV and n 0 = 0.16 fm −3 . With any known relation P = P (ε, n), Eq. 15 can be solved. Furthermore, the collision energy is related to the created density as follows: Here γ CM is the Lorentz gamma factor in the center of mass frame of the heavy ion collisions and E lab is the beam energy per nucleon in the laboratory frame of a fixed target collision. This relation can be obtained from the full stopping condition [52,53,[58][59][60][61][62]. The initial state thermodynamics (density, temperature and entropy) of the hot, dense participant matter is obtained from Eqs. (15) and (16) as a function of the collision energy. The known initial entropy yields the lines of constant entropy which give the trajectories of the heavy ion collisions in the phase diagram. The predicted isentropic expansion trajectories are shown in the T − µ B plane phase diagram in Fig. 7.
Note that 1-dimensional stationary RRHT-adiabat scenario predicts a very strong compression and heating already at intermediate lab (fixed target) bombarding energies. The heavy ion participant system crosses the weak chiral transition predicted by the present CMF model already at E lab ≈ 2 A GeV, i.e. at GSI's SIS18 accelerator facility. Here the specific total entropy is predicted to reach S/A ≈ 3, in accord with previous RMFcalculations [61] which also used the 1-D RRHT-scenario. The T − µ B values, T ≈ 70 MeV, µ B ≈ 1.2 GeV, with net baryon densities n B /n 0 ≈ 3, reached here in heavy ion collisions, coincide with the T − µ B values reached in binary neutron star collisions, as recent general relativistic fully 3+1-dimensional hydrodynamical calculations have confirmed [15,17] for the gravitational wave event GW170817. At these temperatures and densities, T ≈ 70 MeV and n B /n 0 ≈ 3, the RRHT model predicts that there are abut 20% of the dense matter are already transformed to quarks.
At E lab = 5.6 A GeV, √ s NN = 3.5 A GeV roughly 40% of the CMF-matter is in the quark state in the RRHT model -a prerequisite for hot quarkyonic matter. Hence, this energy, which is presumably only reachable by the BMN detector at the Nuclotron at JINR Dubna, and the FXT fixed target setup at the STAR-detector at RHIC, is of great interest: here the matter starts to be dominated by quarks, rather then by in-medium baryons, at T > 100 MeV and µ B 1.5 GeV. This is predicted by the present CMF model, when using the 1-D RRHT ideal hydrodynamics. This model predicts that the quarkyonic transition is crossed also at higher energies, using the isentropic expansion of the matter at specific total entropy S/A > 6.
In fact, non-equilibrium viscous effects may increase the specific entropy of the system. However, pre-freezeout radiation, e.g. of Kaons and other hadrons with small scattering cross-sections which can escape early from the semi-equilibrated baryon-rich, dense system can considerably lower the specific entropy during the expansion. So an answer to the question whether the local entropy per baryon increases or decreases during the time evolution awaits more detailed microscopic/macroscopic modeling. Hence, heavy ions fixed target experiments of SIS at FAiR and SPS at CERN as well as STAR BES program at RHIC probe temperatures from 50 < T < 280 MeV and chemical potentials from 500 < µ B < 1700 MeV for the collision energy range √ s NN < 10 GeV considered here.
In this region the CMF model shows not an additional phase transition, but the remnants of the nuclear liquidvapor transition at T ≈ 20 MeV. The chiral transition at larger chemical potentials can influence the dynamical evolution, too. The present results suggest that heavyion collisions mostly probe regions where the nuclear matter liquid-vapor critical point dominates -hence, the observed baryon fluctuations are largely due to remnants of the nuclear liquid-vapor phase transition. This had been suggested also in previous works [33,[63][64][65]. The CP associated with the chiral symmetry restoration in the CMF model lies at µ B ≈ 1.5 GeV and T ≈ 20 MeV. This high density region is to the best of our knowledge is reachable in the interior of neutron stars, NS, and in binary general relativistic NS mergers [15,16,18,66,67]. Fig. 8 presents the square of the isentropic speed of sound, at fixed specific entropy, c 2 s and the quark frac-tion as function of the temperature for the studied collision energies, i.e. it shows how respective observable quantities evolve during the cooling of the system while it expands. Those isentropic lines which belong to √ s N N > 4.5 GeV probe the softest point of the EoS, which is attributed to the chiral symmetry restoration. At this energy region there are strong local maxima and minima of the speed of sound squared after which c 2 s rapidly increases during the expansion due to the decrease of the quark fraction, as a result of the rapid appearance of baryons, the EoS stiffens quickly doe to the hard-core baryon-baryon repulsion. For collision energies √ s N N < 4.5 GeV the initial state is not dominated by quarks, hence, the systems starts to expand at rather high values of c 2 s , which then monotonously decrease during the expansion, as a result of the diminishing repulsion between the baryons.

VI. APPLICATION TO NEUTRON STARS
The densities in the neutron star interiors also can exceed by several times the nuclear matter saturation density. At these high densities, the lack of a detailed knowledge of the equation of state and the appropriate microscopic degrees of freedom is similarly disturbing as in the relativistic heavy ion collisions discussed above. The discussion of the role of hyperonic, quarkyonic and strange quark degrees of freedom at these NS densities is ongoing.
The CMF model can be employed directly to describe the neutron star matter. Here we work out without any changes to the coupling constants and parameters used to describe the µ B = 0 LQCD results. The temperatures isolated in neutron star interiors are negligibly small, in comparison to what we observe in heavy ions collisions and hot QCD scales. The calculations here are done in the limit T = 0. In contrast to ordinary isospin symmetric nuclear matter, the neutron star matter is in βequilibrium which preserves the total electric neutrality of the NS matter and locally ensures stability with respect to β-decay. As a consequence, strangeness and hypercharges assume finite non-zero values. These constraints require the presence of leptons. In addition, one must allow for quark-hadron degrees of freedom. Fig. 9 depicts the CMF model predictions of the relative particle species number content of all different particle species present inside a CMF neutron star, at T = 0, as function of baryo-chemical potential. One feature of the present CMF model is the absence of baryon resonances (deltas etc. and of hyperons) even though they are included in the CMF model. Their total absence in the present CMF model calculation at T = 0 is due to the very strong hard-core repulsion by the excluded volume corrections.
The calculated EoS at T = 0 can be used as the input for the Tolman-Oppenheimer-Volkoff (TOV) equation, which allows to relate the mass and the radius of any static, spherical, gravitationally bound object [68,69], i.e. here a static neutron star, NS. The outer layer of neutron stars presumably consist of mostly neutron rich nuclei and clusters in chemical and β-equilibrium. Those nuclei are not yet included in the CMF model. Hence, another input for the EoS of the NS crust is needed. Here, we use the classical crust-EoS [70] matched to the CMF-EoS at n B ≈ 0.05 fm −3 . Figure 10 presents the results on the NS mass-radius relations obtained by solving the TOV equation with the present CMF EoS, matched just to that crust-EoS. The total fraction of the star's mass which consists of light and strange quarks is presented in color-code. The most massive stable solution of the TOV equation contains only < 30% deconfined quarks, i.e. for lighter NS only a small fraction of the star's mass originates from deconfined quark matter. If the quark fraction is increased above 30% the stars become unstable. The central density of the stable stars can never exceed n B = 6 n 0 , as shown in the lower part of Fig. 10. Here again the maximum mass indicates the "last stable star". The continuous slow transition from NS matter to a sizable deconfined quark phase implies a smooth appearance of quarks in the star structure and prevents a "second family" of stable solutions to appear. This does prohibit a strict separation between a quark core and the hadronic interior of the star. This is a CMF result due to the Polyakov loop implementation of the deconfinement mechanism. Similar results are obtained in the Quarkyonic Matter-model, where the deconfinement is realized by the appearance of the quarks from inside of the Fermi sea while the hadrons there reside exclusively on the surface shell in momentum space [1]. A similar approach to deconfinement was suggested in [71], however, there the produced mass-radius diagram differs from the CMF-model [34] due to the different realization of the chiral symmetry restoration. The response of a neutron star to non-spherical gravitational fields is reflected in the tidal deformability coefficient λ [72], which depends strongly on the EoS. During the inspiral phase of binary neutron star merger, both neutron stars experience tidal deformations induced by the other respective accompanying neutron star partner. The tidal deformability λ is a measure of the induced quadruple moment Q ij in a response to the external tidal field E ij : λ is directly proportional to the second Love number k 2 : For convenience, usually the dimensionless tidal deforma- bility Λ is presented as: Here, M and R are the mass and radius of the neutron star. A proper value of Λ is important for the description of the inspiral stage during the merger of two neutron stars. Various estimates of Λ emerged after the detection of GW170817 by the LIGO collaboration [73]. Ref [12] argued that for a 1.4M neutron star the tidal deformability and star radius are constrained to Λ 1.4M > 120 and R 1.4M < 13.6 km. It was concluded by means of a Bayesian analysis that for a 1.4M star the deformability should be 375.5 < Λ 1.4M < 800 and the radius is at 12.00 < R 1.4M < 13.45 km, with respective 2 σ confidence levels, see Ref. [13]. A recent analysis by the LIGO and Virgo collaborations [14,74] provides detailed constraints, by using a Bayesian analysis based on the reproducing of the details of the gravitational wave signal.

VII. SUMMARY
A unified and consistent approach to QCD thermodynamics is presented which is appropriate for essentially all temperatures and densities relevant for both heavy ion collisions and for neutron star matter. The Chiral SU(3)flavor parity-doublet Polyakov-loop quark-hadron meanfield model, CMF model, includes the main features of QCD hadron phenomenology, as well as a very good description of known QCD thermodynamics. The CMF model allows for a simultaneous description of many nuclear (astro-) physical data, consistent with astrophysical observations as well as heavy ion collisions and of compact stars. The CMF model is improved by the fixing of the relevant CMF parameters for the quark sector to the state-of-the-art LQCD data on the interaction measure. Here, the parameters of the Polyakov-loop potential and of the quark couplings to the chiral fields have been fixed. A good agreement is found between the CMF model predictions for LQCD data on both the baryon number susceptibilities and the "lines of constant physics". The CMF model is used to explore the phase diagram of strongly-interacting matter at a wide range of T and µ B . Three critical regions are found, which are connected to the nuclear liquid-vapor phase transition, to the chiral symmetry restoration, and to the deconfinement. The region of phase diagram accessible to experiments of high energy heavy-ions collisions is dominated by remnants of the nuclear liquid-vapor phase transition. Other critical regions may be probed by neutron star structure and in binary neutron star mergers. The calculated properties of neutron stars, like the mass-radius relation, the chemical composition of the stars and the tidal deformabilities are in a good agreement with recent experimental observations. The applicability of the improved CMF model to such a wide range of strongly interacting systems is impressive. For the first time, a QCD-motivated EoS is presented which precisely describes the thermodynamic observables for the whole QCD phase diagram.