Isoentropic equations of state of β -stable hadronic matter with a quark phase transition

Weconstructisoentropicequationsofstate(EOSs) of β -stable dense hadronic matter considering the possibility that a quark deconﬁnement phase transition can take place. These conditions can be actually realized in different astrophysical contexts like core-collapse supernovae (CCSNe), during the early stages of the evolution of a newly formed neutron star (protoneutron star, PNS) or in the postmerger compact object formed in binary neutron star (BNS) mergers. We consider four different EOSs to describe the hadronic phase: three EOSs from relativistic mean ﬁeld theory and one EOS recently derived from microscopic calculations in the framework of the Brueckner–Hartree–Fock approach. We combine these hadronic EOSs with a quark matter EOS obtained from a modiﬁed MIT-Bag model which takes into account some perturbative corrections in the grand canonical potential due to the quark–quark interaction. The two phases are then joined up through a Gibbs construction. For each model we study thermal and neutrino trapping effects on the matter composition and consequently on the EOS. We ﬁnally determine the PNS static structure integrating the Tolman– Oppenheimer–Volkoff equations. We ﬁnd that the thermal contribution and particularly the effect of neutrino trapping play an important role on the full EOS. The latter can get softer or stiffer according to the strangeness content in the hadronic phase. These effects are thus crucial to provide a proper description of the dynamical evolution of both the postmerger compact object formed in a BNS merger or the PNS formed in a CCSN.


Introduction
The physics of hot and dense matter is of interest for a large variety of different physical systems spanning from heavyion collisions to binary neutron star (BNS) mergers. The latter have received particular attention in the last years due to their strong connection with the physics of gravitational waves (GWs) and more in general of multimessenger astronomy [1][2][3][4]. The description of neutron stars (NSs) as well as of astrophysical events like core collapse supernovae (CCSNe) [5,6] and BNS mergers [7][8][9] requires the knowledge of the equation of state (EOS) of matter [10,11], namely the relation between pressure (P), energy density (ε) and temperature (T ) for given matter composition. One of the greatest challenge here is the large variety of thermodynamic conditions that matter experiences in these events, spanning several orders of magnitude in all the relevant variables [11,12]. However, according to the specific system under investigation, matter can be characterized by definite thermodynamical conditions like low temperature or approximately constant entropy per baryon (S/A), at least over some timescales of interest. In particular in the present work we consider the specific situation in which S/A can be considered constant and we accordingly study and build up the corresponding EOS. The latter condition is expected to be realized in good approximation in the core of protoneutron stars (PNSs) [10,[13][14][15][16] and in the compact object that forms after the merging of two neutron stars if a prompt collapse to a black hole (BH) does not take place [7,[17][18][19][20]. Such a postmerger remnant can be possibly gravitational unstable and survive only a few milliseconds before collapsing into a BH, but it can also evolve into a long-lived massive NS, with a collapse timescale changing by orders of magnitudes or eventually evolve to a stable NS. The precise fate of the postmerger remnant and the collapse timescale depend on several factors, including the masses and the spins of the colliding NSs, the properties of the EOS at supra-nuclear densities, the strength of the magnetic field and of all processes responsible for the redistribution of angular momentum inside the remnant. Indeed, several mechanisms can temporally prevent the collapse of the central object to a BH and the most relevant ones are the differential rotation of the remnant and the thermal support provided by dense matter in semi-degenerate conditions [7].
Another crucial point for the modeling of the above mentioned astrophysical phenomena is related to the composition of hot and dense matter. In fact due to the very large densities and/or temperatures reached in a PNS or in the postmerger remnant of a BNS merger, various hadronic species (in addition to nucleons) and phases of QCD matter are expected to appear. Hyperons are the first of such "exotic" constituents [21][22][23][24][25]. Moreover, a transition to a phase with deconfined quarks (quark matter) is also expected according to various QCD-inspired models [26][27][28][29][30][31][32][33][34]. The study of the quark deconfinement phase transition in hot and dense matter and of its implications for the structure and evolution of compact objects is one of the main objective of the present work.
Hot and dense matter in the center of CCSNe and BNS mergers produces copious quantities of neutrinos through weak interactions. Over the cooling timescale all neutrinos have escaped, the PNS or the long-lived BNS merger remnant has become transparent and cold and neutrino-less weak equilibrium (also called neutrino-less β-equilibrium) has been achieved. This is the case of cold, deleptonized NSs. However, during the early stages of the evolution of PNSs and BNS merger remnants, the neutrino diffusion time (∼ 30 s) is larger than any relevant dynamical timescale and neutrinos can be considered as trapped (or partially trapped) and possibly in weak equilibrium with matter [7,10,13,15,[35][36][37]). If this physical condition is actually realized for several seconds before neutrinos escape, the trapped neutrino component provides a contribution to the matter's pressure and alters the composition of the stellar material thus modifying the EOS with respect to the case of neutrino-free matter [10].
According to these considerations, it is relevant to ask if thermal and neutrino trapping effects act in direction to support the newly formed hot compact object against the collapse to BH or contrarily favor the collapse.
In the present work we thus study how these physical effects impact on the EOS of hot and dense matter, particularly how they influence the onset of the quark deconfinement phase transition, and as a whole how they affect the evolution of PNSs and BNS postmerger remnants. These effects may turn out to be important when encoded in numerical dynamical simulations of CCSNe and BNS mergers.
We want to emphasize that in the present work we assume that compact stars containing deconfined quarks are hybrid stars; in doing so we do not consider the possibility of the formation of the so called strange stars. The existence of the latter is based on the rightness of the hypothesis about the absolute stability of strange quark matter by Bodmer [38] and Witten [39], and it has been revised in light of the GW signal GW190814, whose light component was interpreted as a strange star of mass around 2.5 − 2.6 M [40] (but see [41] or [42,43] for alternative interpretations).
The paper is organized as follows: in Sect. 2, we briefly review the EOSs adopted for both the hadronic and the quark matter phases. In the Sect. 3, we show results about the composition, temperature profiles and structure of hot compact remnants using isoentropic EOSs and considering various stellar conditions expected during the PNS evolution and/or dynamical astrophysical processes like BNS mergers and CCSNe. In the last section we outline the main conclusions of our work.

Equations of State
In our study we have considered four different EOS models to describe the hadronic phase (i.e. the phase of matter with quarks confined within baryons and mesons). In the first case the EOS is calculated making use of the non-relativistic Brueckner-Hartree-Fock (BHF) many-body approach, whereas the other three EOSs are constructed according to a relativistic mean field (RMF) model. The EOS for the quark phase, i.e. the phase of matter composed of the three lightest quark flavors, namely up (u), down (d) and strange (s) quarks, is modeled using an extended version of the MIT bag model.
For both the hadronic and the quark phases we consider matter in equilibrium, with respect to the weak interaction, with electrons (e − ), muons (μ − ), and eventually with neutrinos (in the case of matter with trapped neutrinos) and we include the contributions of all the corresponding antiparticles.
We next assume a first order phase transition between the two phases and, following Glendenning [44], we require global electric charge neutrality of bulk stellar matter. An important consequence of imposing global charge neutrality is that the hadronic and the quark phases can coexist over a finite range of pressures. This treatment of the phase transition is known in the literature as the Gibbs construction for the hadronic-quark mixed phase.

Hadronic matter EOS
In the first case the hadronic phase is modeled as a uniform fluid of neutrons and protons (nuclear matter) in βequilibrium with electrons and muons, and possibly with neutrinos (in the case of matter with trapped neutrinos). The zero temperature version of this model (hereafter the BL EOS) was derived in Ref. [45] within the BHF manybody approach making use of realistic two-nucleon (NN) [46] and three-nucleon (NNN) [47] interactions derived in in the framework of chiral effective field theory (χ EFT) (e.g. [48,49]). These chiral nuclear interactions reproduce with high accuracy the nucleon-nucleon (NN) scattering data and the experimental binding energies of light (A = 3, 4) atomic nuclei. The BL EOS reproduces the empirical properties (energy per nucleon, symmetry energy and its slope parameter L, incompressibility) of nuclear matter at saturation density (n 0 = 0.16 fm −3 ; see [47]), it does not violate causality (i.e. v s < c, with v s being the speed of sound in the nuclear medium), and it is consistent (see Fig. 2 in [45]) with the measured elliptic flow of matter in heavy-ion collisions experiments [50]. When computing static neutron star configurations, the BL EOS (for the β-stable case) gives a maximum mass M max = 2.08 M (compatible with present measured NS masses [51][52][53][54]), and a quadrupolar tidal polarizability coefficient 1.4 = 385 (for the 1.4 M neutron star [55]) compatible with the constraints derived from GW170817 [2,56]. Based on the universal relation derived in Ref. [20], the threshold mass for the prompt collapse to a BH for equal mass BNS mergers is 2.79 M for this EOS, indicating that GW170817 is compatible with being a NS-NS system that did not undergo a prompt BH formation [18,20,57].
The BL EOS has been recently extended in Ref. [58] to finite temperature and out-of-β-equilibrium matter (hereafter the BLh EOS). This EOS model was successfully used in several numerical simulations of BNS mergers ( [59][60][61][62]).
Moreover, the zero temperature version of this EOS was extended in Ref. [23] to include, in addition to neutrons and protons, the hyperons. In particular, the authors of Ref. [23] have investigated the role of the three-baryon interaction (N N ) between two nucleon and a hyperon on the EOS of hyperonic matter. It was found that it is possible to have hyperonic stars with a maximum mass of ∼ 2 M when the NN interaction is included in the Hamiltonian indicating a possible solution of the so called hyperon puzzle in neutron stars [63,64].
The second, the third and the forth hadronic EOSs we consider are based on a RMF scheme which includes nucleons: DD2 ( [65]) and TM1-2 ( [66]), or nucleons and hyperons: TM1-2Y ( [66]). All the RMF models considered in this work are based on a quantum field theory Lagrangian density in which nucleons are assumed to interact via the exchange of effective σ , ω and ρ mesons. The equation of motion of the various fields are derived in the Hartree approximation. The TM1-2 model is a reparametrization of the original TM1 model ( [67]) where the parameters entering in the Lagrangian density were adjusted (see Ref. [66]) in order to give, for the calculated pressure in symmetric nuclear matter, a result which is consistent with the measured elliptic flow of matter in collision experiments between heavy atomic nuclei [50]. Specifically, in the present work we use the parameters set denoted as TM1-2 in Table 1 of Ref. [68]. One of the main difference between the TM1-2 and the DD2 models is that the latter contains some density dependent terms in the Lagrangian density while the TM1-2 one has none. Both these two EOS models when used in NS structure calculations produce a quadrupolar tidal polarizability coefficient 1.4 compatible with the constraints derived from GW170817 [1,56].
We remark that all the nucleonic EOS models used in our calculations, besides reproducing fairy well the empirical saturation density of nuclear matter n 0 = 0.16 ± 0.01 fm −3 and its corresponding energy per baryon E/A 0 = −16.0 ± 1.0 MeV (e.g., [66]), predict values of the nuclear symmetry energy at saturation density in the interval E sym (n 0 ) = 25 − 37 MeV, and slope parameter L = 30 − 90 MeV, in agreement with several microscopic calculations as well as with experimental results [69,70]. The DD2 EOS model has been also widely employed in literature in numerical simulations of CCSNe [5,6] and BNS mergers [7,8]. Finally, the TM1-2Y is an extension of the TM1-2 model which includes the hyperons of the baryonic octect. The coupling constants that determine the interactions between nucleons and hyperons and between hyperons for this model were set in Ref. [68] using SU(6) symmetry and to reproduce reasonable values of the hyperons' single particle potentials at saturation density in symmetric nuclear matter. The maximum mass of NSs predicted by the TM1-2Y is compatible with measured NS masses having M ∼ 2 M .
At baryon density below n ∼ 0.05 fm −3 one has a transition to non-uniform nuclear matter. In this low density regime the DD2 and the TM1-2 models have a consistent extension up to very low density according to the very same RMF approach. For the BLh model this extension is not straightforward and it would require the solution of the Bethe-Goldstone equation in finite nuclei. In order to avoid this complication, we have instead used the procedure discussed in Ref. [58], namely for the low density part of the EOS we have smoothly matched the microscopic BHF calculations with the SFHo EOS [71].

Quark matter EOS
To model the quark phase we use an extended version of the phenomenological bag model EOS which includes the effects of gluon mediated QCD interactions between quarks up to the second order in the QCD coupling α s [72][73][74]. The grand canonical potential per unit volume can be written as (we use units whereh = c = 1): The first term on the right hand side of Eq. (1) gives the contributions to originating from the fermionic constituents of matter, i.e. f = u, d, s, e − , μ − (and their antiparticles) with the possible addition of neutrinos and antineutrinos (in the case of matter with trapped neutrinos), all described as ideal relativistic Fermi gases. The second term on the right hand side of Eq. (1) accounts for the perturbative QCD corrections to O(α 2 s ) [72][73][74] and its value represents the degree of deviations from an ideal gas EOS, with a 4 = 1 corresponding to the ideal gas case. The chemical potential μ in Eq. (1) can be written in terms of the u, d and s quark chemical potentials as μ = (μ u +μ d +μ s )/3. The term B eff is an effective bag constant which takes into accounts in a phenomenological way nonperturbative aspects of QCD. We assume that both the a 4 and B eff coefficient do not depend on the temperature. The last two terms in Eq. (1) represents the contributions from gluons and photons, both calculated as ideal Bose gases.
In the following we refer to this quark matter EOS model as effective MIT bag model (eMIT). We adopt two different parametrizations of this EOS that we refer to as: eMIT1 and eMIT2. In the first case we set: B eff = 240 MeV and a 4 = 0.7 while in the second: B eff = 190 MeV and a 4 = 0.7. These values have been chosen in such a way that, when the quark EOS is combined with a given hadronic EOS the maximum mass configuration of the resulting hybrid star sequence satisfies the two solar mass requirement [51-54].

EOS for protoneutron stars and postmerger remnants
We construct finite temperature and composition dependent NS EOSs by combining a hadronic and a quark EOS. We assume a smooth Gibbs construction [44] to join the hadronic and the quark phases. We thus require global electric charge neutrality of bulk stellar matter. As a consequence of this, the hadronic and the quark phases can coexist for a finite range of pressures. We are aware that, depending on the nature of the hadron-quark matter interface [75,76], other possibilities exist like to consider a sharp Maxwell construction [77] or to analyze the formation of specific matter geometries in the mixed phase (spheres, rods, ...) that minimize the total energy of the system, as discussed for instance in Ref. [78][79][80].
Moreover, in the case of matter with trapped neutrinos, we include electron and muon neutrinos (and their antiparticles) and describe them as as ideal relativistic Fermi gases in thermal and weak equilibrium with matter. In addition, we have fixed the value of the electronic lepton fraction to Y L e = (n e + n ν e )/n to Y L e = 0.3 and the muonic lepton fraction Y L μ = 0, where n e and n ν e are the net electron and electron neutrino densities, respectively. These values for Y L i are typical ones obtained in hydrodynamical simulations of the early stages of PNS. Throughout the paper we have kept these values of Y L e and Y L μ in the case of trapped neutrinos, unless explicitly specified. In the subsequent evolution, the cooling process is associated with a significant deleptonization towards cold neutrino-less weak equilibrium. It is important to stress that in the case of BNS mergers, matter is characterized by an initially low Y L e ∼ 0.1, and experiences a leptonization process because of matter decompression and temperature increase. We emphasise that, in case of neutrino trapped matter, neutrinos have been included in the full EOS density range and not only in the high density regime.
We remark that we have included the contribution of the photons in all the density range of the EOSs considered in the present work. We additionally remark that even in the case of EOSs that take into account of a deconfinement phase transition, the assumption of constant entropy per baryon is kept. This last contest has been considered in a few works in the past [81][82][83]. A comment is in order: it is well known that during a first order phase transition the entropy per baryon is not constant but suffers a jump together with the total baryonic density of the system [84]. However, as pointed out in the previous works, stellar matter after phase transition can be considered again in a condition of approximately constant value of S/A. Such value will be different and in particular, due to neutrino diffusion, higher than the one before the phase transition. Clearly the precise value of the entropy variation should be calculated by a dynamical simulation while in the present work we can just address some general though qualitative behaviour.
Throughout the paper we adopt the following combinations of hadronic and quark EOS models: BLh+eMIT1, DD2+eMIT2, TM1-2+eMIT2 and TM1-2Y+eMIT2. Notice that when we combine the eMIT quark phase model to the DD2 and TM1-2 hadronic phase models various values of the (B eff , a 4 ) parameters can be adopted (in addition to the eMIT2 set) while keeping the two solar mass limit requirement for the resulting maximum mass hybrid star configuration. This is a consequence of the quite stiff behaviour of the EOS for these two models which allows for a large window of the quark matter parameters. The quark phase parameter space is instead strongly reduced in combination with the BLh model. This in turn reflects the softer character of the BLh model in comparison to the DD2 and TM1-2 ones. However for the analysis that we have carried out, the parameters choice is adequate to highlight the salient features of our discussion.

Results
In Fig. 1 we show the composition (i.e. the particle fractions Y i for the various matter constituents) of β-stable stellar matter as a function of the baryon density n. In the left panel of this figure we present the particle fractions for the BLh+eMIT1 EOS for matter with constant entropy per baryon S/A = 2 in the case of neutrino-free matter (con-  tinuous lines) and neutrino-trapped matter with Y L e = 0.3 (dashed lines). As we can see neutrino trapping shifts the onset of the hadron-quark mixed phase to a larger density with respect to the neutrino-free matter case [10]. This behaviour can be explained as a consequence of the softening of the EOS of the pure nucleonic phase with trapped neutrinos compared to the neutrino-free case. In fact, in neutrinotrapped matter the chemical equilibrium conditions between nucleons, electrons and electron neutrinos reads: to be compared with μ n − μ p = μ e in the case of neutrino free matter. Neutrino trapping, i.e. electronic and muonic lepton number conservation (Y L e = 0.3 and Y L μ = 0 for the corresponding lepton fractions) during the early stages of PNS evolution produces an increase of the net electron fraction and consequently (due to charge neutrality) of the proton fraction (compare the corresponding dashed and continuous lines in Fig. 1 left panel) and makes nuclear matter more proton rich (i.e. more symmetric) and thus softer with respect to the neutrino-free case. The upward pointing arrow in each panel of Fig. 1 denotes the central density of the stellar maximum mass configuration in the case of νfree matter. The corresponding arrow representing the same quantity for ν-trapped matter is very close (see 4th column in Table 1) to the one reported for the case of ν-free matter and it is not shown for clarity. The results reported in the left panel of Fig. 1 thus indicate that in the case of the BLh+eMIT1 EOS quark deconfinement occurs only in the core of stars with a mass close to the maximum mass configuration producing hybrid stars with a inner core with mixed hadron-quark phase. This outcome applies both to the case of neutrino free and neutrino trapped matter. Similar results are obtained for the composition of stellar matter in the case of the DD2+eMIT2 EOS model ( Fig. 1 central panel). Again neutrino trapping shifts the onset of the mixed hadron-quark phase to a larger density compared to the neutrino free matter. Now quark deconfinement occurs at a lower density as compared to the previous case (BLh+eMIT1 EOS), but again the pure quark phase is not present in hybrid stars described by this EOS model. The different values for the onset of the mixed hadron-quark phase for the two models (BLh+eMIT1 and DD2+eMIT2) is essentially a consequence of the larger stiffness of the DD2 EOS with respect to the BLh EOS; this feature favours the onset of the phase transition at lower densities even assuming the same quark matter EOS for the two models.
In the right panel of Fig. 1, we compare the composition of matter with S/A = 2 and trapped neutrinos (dashed lines) with that of cold (T = 0) neutrino free matter (continuous lines) using the DD2+eMIT2 EOS model. These two cases represent the composition of matter in PNSs and in cold deleptonized NSs respectively. From the results in the right panel of Fig. 1, we see that the combined thermal (S/A = 2) and neutrino trapping effects shift the onset of the mixed phase to a slightly higher density with respect to the cold neutrino free case. To disentangle the outcome of the thermal effects on matter composition we can compare the results for S/A = 2 neutrino free matter (continuous lines in the central panel of Fig. 1) with those relative to T = 0 neutrino free matter (continuous lines in the right panel of Fig. 1). From this comparison we see that thermal effects move the threshold density for quark appearance to a lower value with respect to the T = 0 case. Thus thermal and neutrino trap-ping effects have an opposite influence on the value of the baryon density where deconfined quarks appears in neutron star cores.
In Fig. 2 we show the temperature (T S ) of β-stable matter as a function of the baryon density along three different isoentropic paths (S/A = 1, 2, 3). Results in the left (middle) panel have been obtained using the BLh+eMIT1 (DD2+eMIT2) EOS model. The continuous (dashed) lines in both panels refer to the neutrino-free (neutrino trapped) matter case. For each curve in Fig. 2, T S increases up to a maximum value which occurs for a value of the baryon density that is close to the one for the onset of the hadron-quark mixed phase.
The increase with baryon density of the temperature along isoentropic paths in the pure nucleonic phase can be explained using arguments grounded on the properties of mixtures of ideal degenerate Fermi gases [85][86][87]. For example for a degenerate mixture of non-interacting neutrons and protons, with densities n n and n p respectively, neutron-proton asymmetry β = (n n − n p )/n = 1 − 2Y p and neglecting the neutron-proton mass difference (m in Eq. (4) denotes the nucleon mass), the temperature of the system along isoentropic paths can be writen [85,86] as: where Thus for fixed entropy per nucleon and fixed density, T S goes down by a factor 2 −2/3 ∼ 0.63 going from pure neutron matter (β = 1) to symmetric nuclear matter (β = 0).
This argument can be extended to the case of degenerate ideal relativistic Fermi gases [87], as electrons, muons and neutrinos, coming into play in dense stellar matter. The same Fermi gas argument can be used to explain the decrease, at fixed baryon density, of T S when neutrinos are trapped (dashed lines in Fig. 2) with respect to the case of neutrinofree matter (continuous lines in Fig. 2).
As the system enters in the mixed phase, T S decreases ( Fig. 2) with density. Again this behaviour can be understood using the argument based on the properties of mixtures of ideal degenerate Fermi gases and since the entropy is an extensive (i.e. additive) thermodynamic quantity. In fact, in the mixed phase one has the coexistence of nucleons (eventually hyperons), leptons and deconfined quarks and gluons. Thus the number of degrees of freedom in the mixed phase is larger than those in the pure hadronic (nucleonic) phase, accordingly for fixed S/A the temperature T S of the system is lower than the corresponding one for the pure hadronic (nucleonic) phase at the same density. As the density of the mixed phase rises, the quark Fermi seas are increasingly populated whereas the nucleons, electrons and muons concentrations (see Fig. 1) go down. As a result T S decreases with density in the mixed phase. Eventually, if stellar matter could enter in the pure quark phase, T S will rise again with density (see central panel in Fig. 2).
As already pointed out in connection with the matter composition results (Fig. 1), hybrid star configurations computed using the BLh+eMIT1 and the DD2+eMIT2 EOS models have inner cores made of mixed hadron-quark matter. To illustrate this outcome we plot in the left and central panels of Fig. 2  Thus the portion of the T S (n) curve on the right of the full circle can not be realized in hybrid stars described by the EOS models used in Fig. 2. In the right panel of Fig. 2 we show the effect of the variation of the electronic lepton fraction Y L e on T S (with S/A = 2) in the case of βstable matter with trapped neutrinos using the BLh+eMIT1 EOS. Increasing the value of Y L e makes stellar matter more neutrino rich and again, for fixed values of S/A, lower temperatures are obtained. Results similar to those reported in Fig. 2 have been obtained in the case of the TM1-2+eMIT2 and TM1-2Y+eMIT2 EOS models. In Fig. 3 we report the pressure as a function of the energy density for β-stable matter in the case of the BLh+eMIT1 EOS (left panel) and DD2+eMIT2 EOS (right panel) for three different values of the entropy per baryon (S/A = 1, 2, 3) in the case of neutrino-free matter (continuous lines) and matter with trapped neutinos (dashed lines). For both models the mixed hadron-quark phase is delimited by the two kinks in the pressure curves. The heavy dot in each of the continuous lines denotes the values of the central energy density and central pressure of the stellar maximum mass configuration in the case of ν-free matter. The largest outcome of thermal effects on the EOS is observed in the mixed hadron-quark phase. In this region the EOS becomes softer with increasing S/A both in the case of neutrino-free matter (continuous lines) and neutrino-trapped matter (dashed lines) [81]. This behaviour is opposite to what is obtained in the pure hadronic (nucleonic) phase, where thermal effects increase the pressure as the temperature or the entropy per baryon rises. We note however that in the latter case the effect of S/A on the EOS is small compared to that in the mixed phase and it can be barely be appreciated in Fig. 3. This reversal of the thermal effects in the mixed hadron-quark phase is mainly determined by the following two factors (see Fig. 3): (i) the onset of the mixed phase moves to smaller densities when S/A increases (as already noted analysing the results for matter's composition in Fig. 1); (ii) a larger values of S/A favour a larger quark fraction and, as a consequence, a smaller nucleonic content thus making the EOS of the mixed phase softer.
We now discuss the neutrino trapping effects on the EOS. Again the largest consequences of neutrino trapping on the EOS are observed in the mixed hadron-quark phase for all the considered values of S/A. As discussed in connection with the results reported in Fig. 1, neutrino trapping shifts the onset of the mixed phase to larger density with respect to the neutrino free matter case. As shown in Fig. 3 this density shift produces a significant net increase of the pressure in the mixed phase at fixed energy density [81]. In the pure hadronic phase the change of pressure due to neutrino trapping is small compared to the one in the mixed phase with neutrino trapping making the EOS stiffer (softer) at low (high) densities.
In the following lines we discuss the role of hyperons on the quark deconfinement phase transition and on the EOS of hot and dense stellar matter. We will pay special attention to the interplay between hyperon population and neutrino trapping to control the EOS of hybrid stars. For all the results reported in the following Figs. 4 and 5 the EOS for hadronic (hyperonic) phase has been derived using the TM1-2Y model whereas the EOS of the quark phase using the eMIT2 parametrization of the extended version of the MIT bag model.
In the left panel of Fig. 4 we show the composition of βstable hyperonic matter with S/A = 2 , in the case of neutrino free (continuous lines) and neutrino trapped (dashed lines) matter. As we can see neutrino trapping shifts the threshold densities for the various hyperon species to larger values 1 and makes the EOS of hyperonic matter stiffer (compare the dashed and continuous lines in the right panel of Fig. 4) in agreement with previous works (e.g. [10,88]. Again this outcome can be understood comparing the chemical equilibrium conditions in neutrino trapped and neutrino free matter.
In general the chemical potential μ i of a baryon species i in β-stable matter with trapped neutrinos can be written as [10]: where b i and q i are the baryon number and the electric charge (in unit of the elementary electric charge) of the baryon species i, μ n is the neutron chemical potential, and μ e and μ ν e are the chemical potentials for electrons and electron neutrinos respectively. In the case of hyperonic matter using Eq. (5), the chemical potential for the − hyperon in neutrino trapped matter can be written as Since (μ e − μ ν e ) in neutrino trapped matter is smaller than μ e in neutrino free matter [10] and since μ n decreases too, because matter becomes more proton rich (see Fig. 4 left panel), the right hand side of Eq. (6) decreases (at a given baryon density) when neutrinos are trapped. Thus the threshold density for the − hyperons is shifted to a larger density and their abundance decreases with respect to the neutrino free case (Fig. 4 left panel).
In the case of the hyperon (or other charge neutral hyperons) the chemical equilibrium condition (5) gives μ = μ n both in neutrino free and neutrino trapped matter. At first glance it would seem that neutrino trapping has no effect on the condition of chemical equilibrium and therefore on the 1 Looking at Fig. 4 (left panel), it seems that neutrino trapping moves the threshold densities for the − and 0 hyperons to smaller densities with respect to the neutrino free case. However this is a misleading visual effect due to the value for minimum particle fraction Y min = 10 −3 reported in the Y-axis of the figure. Using a smaller Y min (e.g. 10 −5 ) makes clear that neutrino trapping shifts the threshold densities of the − and 0 hyperons to larger values as for the case of the other hyperons of the bayon octet.  hyperons change if we consider the case of trapped or untrapped neutrinos. In fact, as can be seen in Fig. 4 (left panel), neutrino trapping produces an increase of the proton fraction Y p in β-stable nuclear matter and thus reduces the neutron fraction Y n and the neutron chemical potential μ n . Consequently neutrino trapping shifts the threshold density for the hyperons to larger values in agreemment with the results reported the left panel of Fig. 4. As a result (Fig. 4 left panel) neutrino trapping produces a stiffening of the EOS of hyperonic matter [10,88].
We next discuss the interplay between neutrino trapping and the presence of hyperons on the threshold density for the mixed hadron-quark phase and on the EOS of hybrid star matter. First we consider the case with no hyperons and report in the left panel of Fig. 5 the composition of β-stable matter with S/A = 2 using the TM1-2+eMIT2 EOS. In accordance with the results obtained for the BLh+eMIT1 (Fig. 1 left panel) and the DD2+eMIT2 (Fig. 1 central panel) EOS models, again we find that neutrino trapping shifts the onset of the mixed hadron-quark phase to a larger density with respect to the neutrino free case. Afterwards we include hyperons among the constituents of the hadronic phase (NY matter) and use the TM1-2Y model to calculate its EOS. The composition of β-stable matter with S/A = 2 using the TM1-2Y+eMIT2 EOS is depicted in the right panel of Fig.  5. One more neutrino trapping shifts the onset of the hadronquark mixed phase to a larger value compared to the neutrino free case.
To explore the role of hyperons on the quark deconfinement phase transition we compare the composition of β- stable stellar matter without and with hyperons depicted respectively in the left and right panels of Fig. 5. This comparison reveals that the presence of hyperons moves the threshold density for the hadron-quark mixed phase to a higher value with respect to the case where hyperons are not present [10]. This is true both in the case of neutrino free matter (continuous lines in the left and right panel in Fig. 5) and in the case of neutrino trapped matter (dashed lines in the left and central panel in Fig. 5).
The neutron star structure is determined by integrating the Tolman-Oppenheimer-Volkoff (TOV) equations which describe the hydrostatic equilibrium in general relativity for a spherical non-rotating body [26,27]. The calculated gravitational maximum mass (M G,max ), baryonic maximum mass (M B,max ), and the corresponding stellar radius (R) and central baryon density (n c for the various EOS models used in this work are reported in Tab. 1. The gravitational mass as a function of the radius (hereafter the mass-radius relation) is shown in the left and central panel of Fig. 6 for the BLh+eMIT1 and DD2+eMIT2 models respectively for different thermodynamic situations. First of all, it is interesting to note a different behaviour of the two EOS models concerning the dependence on the entropy per baryon of the gravitational maximum mass. In the case of the BLh+eMIT1 EOS model M G,max is almost constant with increasing S/A whereas in the case of the DD2+eMIT2 M G,max decreases with S/A. For example in the case of the DD2+eMIT2 model and neutrino free-matter M G,max decreases of about 5% passing from T = 0 to S/A = 2. As we already pointed out discussing the results in the left panel of Fig. 1, in the case of the BLh+eMIT1 EOS quark deconfinement occurs only in the core of stars with a mass close to the maximum mass configuration producing hybrid stars with a small inner core with mixed hadron-quark phase.
Thus for this EOS model deconfined quarks does not have a sizeable influence on the maximum mass. Consequently the present results for M G,max (S/A) are in line with the finding of Ref. [10,89] for nucleonic stars and accordingly the maximum mass does not change significantly as function of of S/A. We want to remark that adopting the BLh model for the confined phase and employing a different parametrization of the quark matter EOS that allows to obtain a larger quark content in the final hybrid star configurations, a behaviour similar of the DD2+eMIT2 model is recovered. However such alternative parametrizations (like the eMIT2 one for instance) do not allow to satisfy the two solar mass limit and have been therefore ruled out.
In the right panel of Fig. 6 we report the mass-radius relations for the TM1-2 and TM1-2Y models, at fixed entropy per baryon S/A = 2, considering different possible matter compositions as described in details in the figure caption. Concerning the cases in which only nucleonic degrees of freedom are considered the behaviours are very similar to the corresponding ones for the DD2 model. It is interesting to note instead that when hyperons are included in the system, together with a trapped neutrino component, the maximum mass supported is larger than the corresponding one for hyperonic neutrino free matter. This is clearly in agreement with the scenario discussed in Fig. 3. This shows that according to the different composition of hot neutron star matter and in particular to the strangeness content, neutrinos can act in such a way to stabilize the system against collapse to BH like in the case of hyperonic matter, or go in the opposite direction like in the case of nucleonic matter.
The early evolution of a PNS is driven by thermal and neutrino trapping effects on the EOS. The main features of this process can be schematically investigated considering the following two snapshots of the evolution process:  Since most of the matter accretion on the forming NS happens in the very early stages after birth (t < 3 s) [90], the neutron star baryonic mass M B [91] stays almost constant during the evolution between these two configurations. Thus the evolution of a PNS can be unambiguously analyzed in the M G -M B plane [92] (see also [82]). To this purpose we show in Fig. 7 Fig. 7, represents the stellar configuration with a central density equal to the onset of the mixed hadron-quark phase for the considered EOS model and stellar conditions. Thus all neutron star configurations with a baryonic mass larger than the one corresponding to the circle contains a mixed phase core (hybrid stars). We denote the baryonic masses of the stellar configurations associated to these circles as M (i) B,q (red circle) and M ( f ) B,q (black circle) with similar notation for the corresponding stellar gravitational masses.
To begin with we consider the BLh+eMIT1 EOS model (Fig. 7, left panel). A PNS born with a baryonic mass M As shown by our results in Fig. 7 (left panel), these stars will evolve to cold deleptonized nucleonic stars. In other words, in this stellar baryonic mass range neutrino escape and stellar cooling trigger a reconfinement phase transition (i.e. a transition from deconfined quark matter to nucleonic matter) in the stellar core. This peculiar behaviour of the evolution of PNSs, in the above quoted baryonic mass range, can be interpreted comparing the values of the central pressure in the PNS and in the final cold deleptonized NS. For example, considering M (i) B = 2.241 M the central pressure of the star goes down from 452 MeV/fm 3 (PNS) to 232 MeV/fm 3 (NS). At this "low" value of the pressure in cold neutrino free matter, the quark phase has a larger Gibbs energy per nucleon with respect to that of the nucleonic phase, thus a transition to the nucleonic phase take place [93]. We next consider a PNS with In this case the PNS can not be supported by the matter pressure against gravitational collapse since its baryonic mass M B max for the initial configuration. Thus the collapsing stellar core will collapse to a BH after reaching supranuclear densities. In the classical analysis à la Oppenheimer-Volkoff, where the dynamical evolution of the PNS is not taken into account, stars in this baryonic mass range will be considered to have a stable equilibrium configuration in the cold deleptonized stellar sequence. Thus the gravitational mass M * B max ) of the star corresponding to the evolution of the maximum mass PNS configuration, plays the role of an effective maximum mass for the cold deleptonized NS sequence [92]. Finally when M B max the stellar core will collapse to a BH. Consider now a NS in a binary stellar system in which the companion star is a normal star. If during the evolution of the binary system a common envelope is formed the NS can accrete matter from its companion star with a certain accretion rateṀ B . After a sufficiently long time the NS could increase its baryonic mass above the value M B max ] one can have both NSs and BHs. In the case of the BLh+eMIT1 EOS model (Fig. 7, left panel) this matter accretion mechanism is the only way to form hybrid stars, since PNS evolution will always produces nucleonic stars.
A similar analisys can be done in the case of the DD2+eMIT2 EOS model (Fig. 7, central panel). In particular, in the baryonic mass range [M B q ] the quark deconfinement phase transition will take place during the evolution of the newly born star, and a pure nucleonic PNS will evolve to a cold deleptonized hybrid star.
The evolution of PNSs described by TM1-2Y+eMIT2 EOS model (Fig. 7, right panel) shares many common features with those described by the DD2+eMIT2 model, with one major exception: in the case TM1-2Y+eMIT2 EOS model, PNSs with a baryonic mass in the range [M B max ] after neutrino escape and cooling will collapse to black holes. This outcome is related to the presence of hyperons among the stellar constituents in the TM1-2Y model and to role played by neutrino trapping on the EOS. In fact, neutrino trapping (as discussed in the previous pages) shifts the onset of hyperons (and of the quark deconfinement phase transition) to larger density values with respect to those for neutrino free matter, making as a consequence the EOS of PNS matter stiffer and thus resulting in M

Discussion and conclusions
In this work, we have carried out a study of hot β-stable neutrino trapped and neutrino free hadronic matter EOS considering the possibility that a quark deconfinement phase transition can take place. In this study we have considered the case in which the entropy per baryon, S/A, can be analyzed as constant since this condition is expected to be realized in the inner core of neutron stars during some stages of the evolution of high energy astrophysical phenomena like CCSNe or BNSs mergers. We have employed four different EOSs to describe the hadronic phase, namely the BLh, DD2, TM1-2 and TM1-2Y. We have then combined these hadronic EOSs with a quark matter EOS derived from a modified MIT-Bag model which takes into account some perturbative corrections in the grand-canonical potential due to the quark-quark interaction. We have finally performed a Gibbs construction to join up the two phases.
We have found that according to the strangeness content in the hadronic phase, the thermal contribution and the effect of neutrino trapping can get the EOS softer or stiffer. In particular the presence of hyperonic degrees of freedom in neutrino trapped matter produces a stiffer EOS compared to the case of neutrino free matter. In addition, the neutrino trapped component moves the onset of the phase transition to higher densities. The central density corresponding to the maximum mass configuration can accordingly lay in the mixed phase or in the quark phase. For the parametrizations adopted in the present work none of our models shows the formation of a pure quark phase in the final star configurations. Note however that a pure quark phase could be efficiently realized in dynamical BNS merger simulations [94][95][96][97]. Another very important point is that the metastable object that can be formed after the merging of two neutron stars is at least initially in a state of differential rotation. A more realistic description would indeed require the integration of the equations for differentially rotating neutron stars using stellar conditions similar to the ones that we have discussed in this paper. Differentially rotating neutron stars [98][99][100][101] represent indeed the previous stage to the one described in this paper and in some sense trigger the existence of the final star configurations obtained in the present work. These issues will be the topic of a forthcoming work.
Based on our results, we conclude that the effects described in the present paper are crucial for dynamical simulations of BNSMs and CCSNe and may affect the stability of the compact object formed after the neutron star-neutron star merging process.
Funding Open access funding provided by Università di Pisa within the CRUI-CARE Agreement.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data are available upon reasonable request to the corresponding author.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.