1 INTRODUCTION
Neutron star interiors cover an incredible range of densities going from about 104 g cm−3 in the crust to about 1015 g cm−3 in the stellar core, corresponding to several times the nuclear saturation density. Although matter is reasonably understood until a little bit beyond saturation density, unfortunately, not much is known for larger densities. An exception is matter at extremely high densities, beyond the ones reached in the core of neutron stars, when perturbative QCD (PQCD) provides reliable results (Freedman & McLerran Reference Freedman and McLerran1978, Reference Freedman and McLerran1977a, Reference Freedman and McLerran1977b, Reference Freedman and McLerran1977c; Andersen & Strickland Reference Andersen and Strickland2002; Fraga, Kurkela, & Vuorinen Reference Fraga, Kurkela and Vuorinen2014). Since lattice QCD has not yet been extrapolated to high densities, a natural choice of description for such environments relies on effective models.
Effective models, after being calibrated to work on a certain regime of energies, can produce reliable results concerning the matter equation of state (EoS) and associated particle population, which in the case of neutron stars can be further used in dynamical simulations. These simulations include core collapse supernova explosions, star cooling, and compact star mergers and require tabulated data with information about the mycrophysics as input. In this article, I present tables calculated within the chiral mean field (CMF) model that can be used for this purpose.
2 THE CMF MODEL
Since the same physical laws govern particles under all conditions, it is only logical that a model used to describe cold and dense stars should also be compatible with the description of hot environments, such as the ones created in supernova explosions and compact star mergers, and ultimately in heavy-ion collisions or the early universe. After all, these environments are only different facets of a larger picture of matter with high energy, represented in the QCD phase diagram. Although fully evolved neutron stars have temperatures ≲1 MeV, proto-neutron stars can reach temperatures of about 30 MeV or more in their centres (Burrows & Lattimer Reference Burrows and Lattimer1986; Pons et al. Reference Pons, Reddy, Prakash, Lattimer and Miralles1999) and, in compact star mergers, temperatures of approximately 80 MeV can be reached (Galeazzi et al. Reference Galeazzi, Kastaun, Rezzolla and Font2013). Such temperatures are not significantly different from temperatures reached in high-energy heavy-ion collisions (such as the ones performed in the RHIC and LHC particle colliders (Schenke, Jeon, & Gale Reference Schenke, Jeon and Gale2011; Alqahtani et al. Reference Alqahtani, Nopoush, Ryblewski and Strickland2017). Having those numbers in mind, it becomes natural to use the same mathematical model, or at least the same kind of approach to describe all such systems.
The CMF model is based on a non-linear realization of the SU(3) sigma model. It is an effective quantum relativistic model that describes hadrons interacting via meson exchange and it is constructed in a chirally invariant manner, as the particle masses originate from interactions with the medium and, therefore, decrease at high densities/temperatures. The non-linear realization of the sigma model is an improvement over the widely used sigma model and it includes the pseudoscalar mesons as the angular parameters for the chiral transformation. In this case, the pseudoscalar mesons exhibit a pseudovector coupling to the baryons in agreement with the experimental finding of a vanishing π − N scattering length, chiral invariance for heavy particles is ensured if their coupling is invariant under local SU(3) vector transformations (allowing couplings between baryons and meson octets), a connection to the phenomenological Walecka model exists, the masses of the pseudoscalar mesons do not become imaginary at high densities, etc. As a consequence, in addition to reproducing chiral symmetry restoration, the model is in very good agreement with nuclear physics data (Papazoglou et al. Reference Papazoglou, Zschiesche, Schramm, Schaffner-Bielich, Stoecker and Greiner1999).
The Lagrangian density of the CMF model within the mean field approximation reads
$$\begin{eqnarray}
\mathcal {L} = \mathcal {L}_{\rm {Kin}}+\mathcal {L}_{\rm {Int}}+\mathcal {L}_{\rm {Self}}+\mathcal {L}_{\rm {SB}},
\end{eqnarray}$$
where, besides the kinetic energy term for baryons (entire octet) and free leptons (included to ensure charge neutrality), the terms
$$\begin{eqnarray}
\mathcal {L}_{\rm {Int}}&=&-\sum _i \bar{\psi _i}[\gamma _0(g_{i\omega }\omega +g_{i\phi }\phi +g_{i\rho }\tau _3\rho )+M_i^*]\psi _i,\nonumber \\
\mathcal {L}_{\rm {Self}}&=&+\frac{1}{2}(m_\omega ^2\omega ^2+m_\rho ^2\rho ^2+m_\phi ^2\phi ^2)\nonumber \\
&&+\,g_4 \left(\omega ^4 + 3 \omega ^2 \phi ^2 + \frac{\phi ^4}{4} + \frac{4 \omega ^3 \phi }{\sqrt{2}} + \frac{2 \omega \phi ^3}{\sqrt{2}}\right)\nonumber \\
&&-\,k_0(\sigma ^2+\zeta ^2+\delta ^2)-k_1(\sigma ^2+\zeta ^2+\delta ^2)^2\nonumber \\
&&-\,k_2\left(\frac{\sigma ^4}{2}+\frac{\delta ^4}{2}+3\sigma ^2\delta ^2+\zeta ^4\right)-k_3(\sigma ^2-\delta ^2)\zeta \nonumber \\
&&-\,k_4\ \ \ln {\frac{(\sigma ^2-\delta ^2)\zeta }{\sigma _0^2\zeta _0}}\,,\nonumber \\
\mathcal {L}_{\rm {SB}}&=&-m_\pi ^2 f_\pi \sigma -\left(\sqrt{2}m_k^ 2f_k-\frac{1}{\sqrt{2}}m_\pi ^ 2 f_\pi \right)\zeta \,,
\end{eqnarray}$$
represent the interactions between baryons and vector and scalar mesons, the self-interactions of scalar and vector mesons, and an explicit chiral symmetry breaking term, which is responsible for producing the masses of the pseudo-scalar mesons. The mesons included are the vector-isoscalars ω and ϕ (strange quark–antiquark state), the vector-isovector ρ, the scalar-isoscalars σ and ζ (strange quark–antiquark state) and the scalar-isovector δ, with τ3 being twice the isospin projection operator of each particle (±1). The isovector mesons affect isospin-asymmetric matter and, thus, are important for neutron star physics. Also, the δ meson has a contrary but complementary role to the ρ meson, much like the σ and ω mesons. The finite-temperature calculations include a heat bath of hadronic and quark quasiparticles and their antiparticles within the grand canonical potential of the system.
The effective masses for the baryons are generated by the scalar mesons:
$$\begin{eqnarray}
M_{i}^*&=&g_{i\sigma }\sigma +g_{i\delta }\tau _3\delta +g_{i\zeta }\zeta +M_{0_i}\,,
\end{eqnarray}$$
with the exception of small explicit mass terms M 0 N = 151.68 and M 0Λ, Σ, Ξ = 354.91 MeV.
The coupling constants of the model were presented in Dexheimer & Schramm (Reference Dexheimer and Schramm2008) and are shown in Table 1. The scalar sector was fitted to reproduce the vacuum masses of baryons and mesons and the pion and kaon decay constants. The vector sector was fitted to reproduce nuclear constraints for symmetric matter at saturation, such as baryon density (ρ0 = 0.15 fm−3), binding energy per nucleon (B/A = −16 MeV), compressibility (K = 300 MeV), as well as symmetry energy (E
sym = 30 MeV), and symmetry energy slope (L = 88 MeV). The reproduced pressure and compressibility for neutron matter at saturation are P = 4.56 MeV fm−3 and K = 870 MeV. The reproduced hyperon potentials at saturation are
$U_\Lambda =-28$
MeV,
$U_\Sigma =5$
MeV,
$U_\Xi =-18$
MeV. The reproduced critical point for the nuclear liquid–gas phase transition lies at T
c
= 16.4 MeV, n
B
c
= 0.05 fm−3, μ
B
c
= 910 MeV.
Coupling constants for the model, using χ = 401.93 MeV.
Regarding the constraints for the slope of the symmetry energy L, although a compilation of several studies indicates values lower than 60 MeV (Lattimer & Lim Reference Lattimer and Lim2013), there are other works that indicate that the values of such quantity should be larger than 90 MeV (Cozma et al. Reference Cozma, Leifels, Trautmann, Li and Russotto2013; Chen, Ko, & Li Reference Chen, Ko and Li2005; Sotani et al. Reference Sotani, Iida and Oyamatsu2015; Tsang et al. Reference Tsang2012; Wang et al. Reference Wang, Guo, Li, Zhang, Leifels and Trautmann2014). It is also important to note that the values suggested in Lattimer & Lim (Reference Lattimer and Lim2013) are a result of a compilation of experimental analyses that have different systematic and statistical errors, and hence should be interpreted carefully when it comes to excluding equations of state.
The numerical code for the CMF model solves a set of equations for each baryon chemical potential and temperature. Those include an equation of motion for each meson and one extra equation in the case that baryon number density is fixed (instead of chemical potential). Additional constraints such as, for example, charge neutrality, fixed charge fraction, fixed lepton fraction, fixed entropy per baryon and zero net isospin require additional equations.
In order to study neutron stars, charge neutrality and chemical equilibrium are required. As a result of the energy balance, highly isospin-asymmetric objects are formed. Nevertheless, in proto-neutron stars, the proton-to-neutron ratio is larger than in cold neutron stars, making these systems more similar to heavy-ion collision environments. Studies of the CMF model including fixed entropy per baryon together with trapped neutrinos were able to reproduce massive neutron stars (like the ones that have been recently observed (Antoniadis et al. Reference Antoniadis2013; Demorest et al. Reference Demorest, Pennucci, Ransom, Roberts and Hessels2010)) even taking hyperons into account (Dexheimer & Schramm Reference Dexheimer and Schramm2008). The CMF formalism was also used to study the effect of kaon condensation in neutron and proto-neutron stars (Mishra et al. Reference Mishra, Kumar, Sanyal, Dexheimer and Schramm2010), the inclusion of chiral partners in stars (Dexheimer, Schramm, & Zschiesche Reference Dexheimer, Schramm and Zschiesche2008b; Dexheimer et al. Reference Dexheimer, Pagliara, Tolos, Schaffner-Bielich and Schramm2008a; Dexheimer, Negreiros, & Schramm Reference Dexheimer, Negreiros and Schramm2015a; Mukherjee et al. in preparation), and the cooling profile of stars (Negreiros, Dexheimer, & Schramm Reference Negreiros, Dexheimer and Schramm2010; Dexheimer, Negreiros, & Schramm Reference Dexheimer, Negreiros and Schramm2015b), which was (and still is) in good agreement with observed data. For cold chemically equilibrated stars, a maximum mass star with 2.1 M⊙ and corresponding radius of 12 km is reproduced (1.93 M⊙ and 13 km when quarks are included). For the canonical star with mass 1.4 M⊙, a corresponding radius of 14 km is found.
Up, down, and strange quarks were introduced in the formalism within the same model (Dexheimer & Schramm Reference Dexheimer and Schramm2010; Negreiros et al. Reference Negreiros, Dexheimer and Schramm2010; Hempel et al. Reference Hempel, Dexheimer, Schramm and Iosilevskiy2013); however, the degrees of freedom, which are actually populated at a certain temperature and density, change from hadrons to quarks and vice-versa through the introduction of an extra field Φ in the effective masses of baryons and quarks (Dexheimer & Schramm Reference Dexheimer and Schramm2008, Equations (6) and (7)). The scalar field Φ was named in analogy with the Polyakov loop (Fukushima Reference Fukushima2004), since it also functions as the order parameter for deconfinement. The potential for Φ (see Dexheimer & Schramm Reference Dexheimer and Schramm2008, Equation (9) was modified from its original form in the Polyakov-loop-extended Nambu–Jona–Lasinio (PNJL) model (Ratti, Thaler, & Weise Reference Ratti, Thaler and Weise2006; Roessner, Ratti, & Weise Reference Roessner, Ratti and Weise2007) in order to be used to study low-temperature and high-density environments (in addition to high-temperature and low-density environments). It was shown in Fukushima (Reference Fukushima2011), Lourenco et al. (Reference Lourenco, Dutra, Delfino and Malheiro2011, Reference Lourenco, Dutra, Frederico, Delfino and Malheiro2012), and Blaschke et al. (Reference Blaschke, Alvarez Castillo, Benic, Contrera and Lastowiecki2013) that this choice for the potential U(Φ, T, μ B ) can also be used in the Polyakov–Nambu–Jona–Lasinio model, successfully reproducing QCD features. Nevertheless, the CMF model is significantly different from the widely used PNJL model (Fukushima Reference Fukushima2004). In our case, both hadronic and quark (in addition to leptonic) degrees of freedom are included. Note that although a relation between the effective mass for the fermions and the chiral condensate can be written using mean field theory in both models, such a derivation is quite different in each case.
The QCD phase diagram shown in Figure 1 was constructed by analysing the behaviour of the deconfinement to quark matter (and associated chiral symmetry restoration) at different temperatures and chemical potentials (or densities) using the formalism described above for isospin-symmetric matter (zero isospin and strangeness constraints) and neutron-star matter (charge neutral and in chemical equilibrium). The lines represent first-order phase transitions, meaning that for each baryon chemical potential and temperature, there are multiple metastable solutions in order-parameter space, although only one truly stable phase exists. The circles mark critical points, beyond which the deconfinement and chiral transitions become smooth crossovers. The shaded regions exemplify some of the different regimes that can be described within the CMF model (that can be applied to the entire μ B –T plane). Note that at low density/temperature, the nuclear physics liquid–gas phase transition is also reproduced.
QCD phase diagram resulting from the CMF model. The lines represent first-order transitions. The circles mark the critical end-points. Isospin-symmetric matter refers to zero isospin and strangeness constraints, while neutron-star matter stands for charged neutral matter in chemical equilibrium. The shaded regions exemplify some of the different regimes that can be described within the model.
Since the coupling constants related to quark matter cannot be related to nuclear properties, they were constrained using lattice QCD data as well as information about the QCD phase diagram from Ratti et al. (Reference Ratti, Thaler and Weise2006), Roessner et al. (Reference Roessner, Ratti and Weise2007), Aoki et al. (Reference Aoki, Endrodi, Fodor, Katz and Szabo2006), and Fodor & Katz (Reference Fodor and Katz2004) for symmetric matter (see Dexheimer & Schramm Reference Dexheimer and Schramm2010, Table 2). In this approach, the chiral symmetry restoration and deconfinement phase transitions range from sharp first-order phase transitions to smooth crossovers as the temperature increases. With the advent of RHIC and LHC, relativistic heavy-ion collision experiments have focused on the low-μ B part of the QCD phase diagram. The RHIC Beam Energy Scan is currently investigating the central region of the diagram, and future facilities, like FAIR in GSI and NICA in Dubna, will be capable of a more in-depth exploration of the high-μ B region. Lattice QCD calculations are also slowly advancing towards the middle of the diagram (from the left). Nevertheless, only neutron stars will be able to probe the very right side of the QCD phase diagram due to their incredibly high chemical potential to temperature ratio.
The most important and unique aspect of our description is that hadrons are included as quasi-particle degrees of freedom in a chemically equilibrated mixture with quarks. Therefore, the model gives a quasi-chemical representation of the deconfinement phase transition (so-called chemical picture in terms of electromagnetic non-ideal plasmas (Iosilevskiy Reference Iosilevskiy and Fortov2000)). The assumed full miscibility of hadrons and quarks is, for example, in contrast to the underlying picture of simple quark-bag models. At sufficiently high temperature, this will lead to the appearance of quarks in the hadronic sea. On the other hand, it is also possible that some hadrons survive in the quark sea. Nevertheless, quarks will always give the dominant contribution in the quark phase, and hadrons in the hadronic phase. The hadronic and the quark phases are characterized and distinguished from each other by the values of the order parameters, σ and Φ. The inter-penetration of quarks and hadrons in the two phases is physical, and is required to obtain a true crossover transition at low baryon chemical potential, which has been shown to exist by lattice QCD calculations.
It is important to note at this point, that astrophysical EoS’s including quark deconfinement are usually only carried out up to a few tens of MeV’s, if they include finite-temperature effects at all. Different EoS’s (representing different degrees of freedom) are put together ‘by hand’ in different regions of the neutron star, generating necessarily first-order phase transitions between them. This approach is not compatible with QCD calculations such as the one in Baym et al. (Reference Baym, Hatsuda, Tachibana and Yamamoto2008), Lourenco et al. (Reference Lourenco, Dutra, Frederico, Delfino and Malheiro2012), Kojo et al. (Reference Kojo, Powell, Song and Baym2015), and Masuda, Hatsuda, & Takatsuka (Reference Masuda, Hatsuda and Takatsuka2016), which suggests that the phase transition to quark matter might be a crossover, even at low temperatures. Unlike this approach, our description predicts different degrees of freedom appearing self-consistently at different densities and temperatures. This allowed for the first detailed comparison between the nature of the deconfinement phase transition and the one of the nuclear matter liquid–gas phase transition (Hempel et al. Reference Hempel, Dexheimer, Schramm and Iosilevskiy2013). It was found in this work (in agreement with Bombaci et al. Reference Bombaci, Logoteta, Panda, Providencia and Vidana2009) that the principle difference between both phase transitions is that, in contrast to the ordinary van der Waals like phase transition, the phase coexistence line of the deconfinement phase transition has a negative slope in the pressure–temperature plane. This feature is related to the quark phase having higher entropy per baryon than the hadronic phase. As another qualitative difference, we found that the non-congruent features of the deconfinement phase transition become vanishingly small around the critical point. Non-congruent phase transitions occur for first-order phase transitions with more than one globally conserved charge, allowing local concentrations of the charges to vary during a phase transition, i.e., the appearance of mixtures of phases.
As a final test of the validity of the CMF model, I present here for the first time a figure showing the speed of sound reproduced by the model when a mixture of hadronic and quark phases is allowed. Naturally, as in any relativistic formalism, the speed of sound does not go above 1 (the value of the speed of light in natural units), but note that neither does our speed of sound go above the speed of sound limit provided by kinetic theory (Moustakidis et al. Reference Moustakidis, Gaitanos, Margaritis and Lalazissis2017) given by
$$\begin{eqnarray}
v_s&=&\sqrt{\frac{\varepsilon -P/3}{P+\varepsilon }}\,.
\end{eqnarray}$$
Finally, it is important to note that the EoS presented in Figure 2 (black line) has already been successfully compared with PQCD calculations in Dexheimer, Franzon, & Schramm (Reference Dexheimer, Franzon and Schramm2017). Other models commonly used to describe quark matter in stars, such as the bag model, do not agree with PQCD calculations (Fraga et al. Reference Fraga, Kurkela and Vuorinen2014).
Equation of state for neutron-star matter at zero temperature within the CMF model, its derivative and square root of derivative (speed of sound). The kinetic limit for the speed of sound from Moustakidis et al. (Reference Moustakidis, Gaitanos, Margaritis and Lalazissis2017) is also shown. The arrows mark the threshold for the appearance of the Lambdas and strange quarks.
3 TABULATED EoS
The tables described in this section contain the EoS (a large set of thermodynamical quantities) and particle population for different sets of conditions calculated using the CMF model. These tables are available in two locations, in the CompOSE website http://compose.obspm.fr/ and in the following website http://personal.kent.edu/~vdexheim/. At this point in time, there are three different kinds of tables available produced using the CMF model, all containing only hadronic degrees of freedom (and in some cases leptons). In the future, equivalent tables will be uploaded that also contain quark degrees of freedom. The tables, produced under the assumption of different conditions, are the following:
-
• One-dimensional (1D) tables for chemically equilibrated neutron stars at zero temperature. In this case, the baryon number density ranges from 0.03 to 3.03 fm−3 with steps of 0.01 fm−3. They contain contributions from nucleons, hyperons, electrons, and muons.
-
• Three-dimensional (3D) tables for supernova/neutron-star merger simulations without leptons. In this case, the charge fraction ranges from 0 to 0.53 with steps of 0.01, the temperature ranges from 0 to 160 MeV with steps of 2 MeV, and the baryon number density ranges from 0.01 to 3.01 fm−3 with steps of 0.01 fm−3. They contain contributions from nucleons and hyperons.
-
• 3D tables for supernova/neutron-star merger simulations with electrons. In this case, the charge fraction ranges from 0 to 0.53 with steps of 0.01, the temperature ranges from 0 to 160 MeV with steps of 2 MeV, and the baryon number density ranges from 0.01 to 3.01 fm−3 with steps of 0.01 fm−3. They contain contributions from nucleons, hyperons, and electrons.
For each of those sets of conditions, five tables are provided. They are CompOSE standard data files (more information about them can be found in the following website http://compose.obspm.fr/manual/) and include the following:
-
• eos.thermo: a table with 2 dimensions, 1 316 574 grid points, and 1 316 575 lines (the first line gives the masses of the neutron and proton). The standard columns contain the pressure divided by the baryon number density (in MeV), entropy density divided by the baryon number density, scaled and shifted baryon chemical potential, scaled charge chemical potential, scaled effective lepton chemical potential, scaled free energy density divided by the baryon number density, and scaled internal energy density divided by the baryon number density. Scaled means divided by the nucleon mass, and shifted means it had 1 subtracted. The two extra columns contain the scaled enthalpy per baryon number density and strangeness number density (∑ i s i n i in fm−3). For the tables containing a quark phase, a third extra column containing the quantity ∑ i Q B i n i used to calculate the total hadronic electric charge will be added (see discussion of Equations (7) and (11)).
-
• eos.compo: a table with 2 dimensions, 1316 574 grid points, and 1 316 574 lines. The standard columns contain the particle fractions (n i /n B ) and the index encoding the type of phase, which in our case is 1 for hadronic matter and 2 for quark matter.
-
• eos.nb: a table with 1 dimension, 301 baryon number density grid values (in fm−3), and 303 lines (the first two lines contain the initial and final grid points).
-
• eos.t: a table with 1 dimension, 81 temperature grid values (in MeV), and 83 lines (the first two lines contain the initial and final grid points).
-
• eos.yq: a table with 1 dimension, 54 charge fraction grid values, and 56 lines (the first two lines contain the initial and final grid points).
In the case of the 1D tables for chemically equilibrated neutron stars at zero temperature, charge neutrality is fulfilled:
$$\begin{eqnarray}
\sum _i {Q_e}_i n_i=0,
\end{eqnarray}$$
where Q e i is the electric charge and n i the number density of each baryon, lepton, or quark and chemical equilibrium is ensured through
$$\begin{eqnarray}
\mu _i={Q_B}_i \mu _B+ {Q_e}_i \mu _q,
\end{eqnarray}$$
where μ i is the chemical potential of each baryon or quark, Q B i is the baryon number, μ B the baryon chemical potential, and μ q the charged chemical potential, which is equal to minus the electron/muon chemical potentials.
In the case of the 3D tables for supernova/neutron-star merger simulations without leptons, charge neutrality is not required and the leptons do not participate in chemical equilibrium, although Equation (6) still holds and μ q is only equal to zero in isospin-symmetric matter (as Coulomb interactions are not taken into account). The charge fraction is defined as the total charge over baryon number:
$$\begin{eqnarray}
Y_q=\frac{Q}{B}=\frac{\sum _i Q_{e_i} n_i}{\sum _i {Q_B}_i n_i} ,
\end{eqnarray}$$
where ∑Q B i n i is not the same as baryon number density n B , as the latter comes from the derivate of the pressure with respect to the baryon chemical potential and, therefore, also contains a contribution from the potential U for Φ (when quarks are included). The second law of thermodynamics can be written for the baryons and quarks as
$$\begin{eqnarray}
\sum _i \varepsilon _{i}=-\sum _i P_i+T\sum _i s_i+\sum _i \mu _i n_i,
\end{eqnarray}$$
where ε i , P i , and s i are the energy density, pressure, and entropy density of each baryon or quark and T is the temperature. This expression can be rewritten using Equation (6) and the definition of Y q :
$$\begin{equation}
\sum _i \varepsilon _{i} = -\sum _i P_i+T\sum _i s_i+\sum _i ({Q_B}_i \mu _B+ {Q_e}_i \mu _q) n_i,
\end{equation}$$\\
$$\begin{equation}
\sum \varepsilon _i = -\sum _i P_i+T\sum _i s_i+\widetilde{\mu }\sum _i {Q_B}_i n_i,
\end{equation}$$
$\widetilde{\mu }$
is the free energy of the system defined as
$\widetilde{\mu }=\mu _B+ \mu _q Y_q$
. For more details on this derivation, see Hempel et al. (Reference Hempel, Dexheimer, Schramm and Iosilevskiy2013). Note that the total energy density, pressure, and entropy density of the system also contain mesonic and Φ contributions.
In the case of the 3D tables for supernova/neutron-star merger simulations with electrons, the electrons still do not participate in chemical equilibrium, although, as before, Equation (6) still holds for the baryons and quarks and μ q is only equal to zero for isospin symmetric matter. The electron density is determined in order to establish charge neutrality, so from Equation (7) one can derive
$$\begin{eqnarray}
\sum _i {Q_e}_i n_i={\sum _i {Q_B}_i n_i}Y_q=n_e,
\end{eqnarray}$$
where n e is the number density of electrons.
Note that there are no nuclei included in our calculations within the CMF model, as its current parametrization describes only bulk matter. A version of the CMF model, which can describe hot and dense matter but also includes nuclei (as in Papazoglou et al. Reference Papazoglou, Zschiesche, Schramm, Schaffner-Bielich, Stoecker and Greiner1999; Beckmann et al. Reference Beckmann, Papazoglou, Zschiesche, Schramm, Stoecker and Greiner2002; Schramm Reference Schramm2002, Reference Schramm2015), is under way and will be reported in the near future. For this reason, a star crust should be added to the tables presented in this article before using them to perform realistic star simulations.
As an example of our results containing quarks (without the assumption of mixtures of phases), Figures 3 and 4 show the particle population for the extreme cases of charge-neutral mater (Y
q
= 0) and half-charged matter (Y
q
= 0.5) at zero temperature without the inclusion of leptons. In both figures, the y-axis contains the baryon number density, meaning that the quark densities were multiplied by 1/3. In Figure 3, one can see that all hyperons but the Lambdas are suppressed by the appearance of the quarks. In Figure 4, one can see that all hyperons are suppressed and that the quarks appear at larger free energy
$\widetilde{\mu }$
.
Particle population for charge-neutral matter without leptons (Y q = 0) at zero temperature within the CMF model.
Particle population for half-charged matter without leptons (Y q = 0.5) at zero temperature within the CMF model.
In Figure 3, the free energy equals the baryon chemical potential in each phase
$\widetilde{\mu }=\mu _{B_H}=\mu _{B_Q}$
, but, in Figure 4, the free energy is different from the baryon chemical potential in each phase
$\widetilde{\mu }=\mu _{B_H}+ 0.5 \mu _{q_H}=\mu _{B_Q}+ 0.5 \mu _{q_Q}$
. For more details on the correspondence of different chemical potentials in different phases, see Hempel et al. (Reference Hempel, Dexheimer, Schramm and Iosilevskiy2013). Finally, note that fixed charge fraction Y
q
= 0.5 does not correspond to isospin symmetric matter if net strangeness is not set to zero, which is the case here. This can easily be verified by the presence of strange quarks in Figure 4.
4 CONCLUSION
I presented in this article a review of some of the accomplishments of the CMF model in describing dense and/or hot matter and, in particular, neutron stars under different stages of evolution. These accomplishments also include a possible description of matter produced in heavy-ion collisions (Steinheimer et al. Reference Steinheimer, Dexheimer, Petersen, Bleicher, Schramm and Stoecker2010).
For the first time, a tabulated version of the EoS and particle population produced within the CMF model under different conditions was presented. These tables can be used in simulations of, for example, core-collapse stellar explosion, stellar cooling, and mergers.
Although the CMF tables that are already available online only contain hadronic matter (a version with quarks will be available in the near future), their results contain essential features of the description of hot and/or dense matter, such as chiral symmetry restoration, the inclusion of hyperons, and the inclusion of antiparticles for all fermions. In addition, since the CMF model results derive from a relativistic description, a subluminal speed of sound is guaranteed in any regime.
ACKNOWLEDGEMENTS
The author would like to thank Matthias Hempel for numerous helpful discussions. The author acknowledges the support from NewCompStar COST Action MP1304.
