Issue 
A&A
Volume 518, JulyAugust 2010
Herschel: the first science highlights



Article Number  A17  
Number of page(s)  8  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201014308  
Published online  24 August 2010 
The maximum and minimum mass of protoneutron stars in the Brueckner theory
G. F. Burgio and H.J. Schulze
INFN Sezione di Catania and Dipartimento di Fisica, Universitá di Catania, via Santa Sofia 64, 95123 Catania, Italy
Received 24 February 2010 / Accepted 17 April 2010
Abstract
We study the structure of protoneutron stars within the
finitetemperature BruecknerBetheGoldstone theoretical approach,
paying particular attention to how it is joined to a lowdensity
nuclear equation of state. We find a slight sensitivity of the minimum
value of the protoneutron star mass to the lowdensity equation of
state, whereas the maximum mass is hardly affected.
Key words: dense matter  equation of state  stars: interiors  stars: neutron
1 Introduction
A lot of effort is currently dedicated to dynamical simulations of supernovae explosions (Burrows & Lattimer 1986; Pons et al. 1999; Villain et al. 2004; Liebendörfer et al. 2005; Fischer et al. 2009) and the subsequent formation and evolution of a protoneutron star (PNS). The fundamental input to these calculations is the nuclear equation of state (EOS) over a wide range of densities, apart from microscopic information regarding diffusion and cooling processes. The output are timedependent radial profiles of the thermodynamic quantities of interest, such as temperature, entropy, particle fractions, etc.
According to current understanding, one can very roughly identify different stages of this process (Burrows & Lattimer 1986; Pons et al. 1999; Strobel et al. 1999). (i) For about 1 s after supernova core bounce, the system consists of a relatively cool central region surrounded by a hot mantle, collapsing and radiating off neutrinos quickly, while still also accreting material. (ii) For the next 20 or so seconds one can identify a slowly developing state of the proper PNS, the system first deleptonizing and heating up the interior parts of the star in the process, and then beginning to cool down by further neutrino diffusion. (iii) After several tens of seconds, the final state of a cold NS has basically been formed, which continues to cool down slowly first by neutrino and later by photon emission.
In this article we do not intend to perform dynamical simulations, but focus on the consistent construction and evaluation of basic consequences of the temperaturedependent nuclear EOS during the prominent PNS stage (ii). We therefore assume strongly idealized, static profiles representing this environment. For simplicity we use a constant entropy per baryon throughout the star and investigate the sensitivity of the results to the chosen value of entropy S/A, as is usually done (Gondek et al. 1997; Goussard et al. 1997).
Another important aspect is the treatment of neutrino trapping. Gravitational collapse calculations of the whitedwarf core of massive stars indicate that, at the onset of trapping, the electron lepton fraction is , the precise value depending on the efficiency of electron capture reactions during the initial collapse stage. Also, because no muons are present when neutrinos become trapped, the constraint is imposed. We fix the at these values in our calculations for neutrinotrapped matter. Furthermore, since lowdensity nuclear matter loses neutrino very fast, the concept of a neutrino sphere, inside which neutrinos are trapped, is often introduced during the PNS stage (Gondek et al. 1997; Fischer et al. 2009). Our treatment of the neutrino sphere is explained later in Sect. 4.
One of the most advanced microscopic approaches to the EOS of nuclear matter is the Brueckner theory. In recent years, it has made rapid progress in several aspects. (i) The convergence of the BruecknerBetheGoldstone (BBG) expansion has been firmly established (Song et al. 1998; Baldo et al. 2000b, 2001). (ii) The addition of phenomenological threebody forces (TBF) based on the Urbana model (Carlson et al. 1983; Schiavilla et al. 1986), permitted improving the agreement to a large extent with the empirical saturation properties (Baldo et al. 1997; Zhou et al. 2004; Li et al. 2006, 2008a). (iii) The BruecknerHartreeFock (BHF) approach has been extended in a fully microscopic and selfconsistent way to describe nuclear matter also containing hyperons (Schulze et al. 1995, 1998). (iv) It has also been extended to the finitetemperature regime within the BlochDe Dominicis formalism with a realistic nucleonnucleon interaction (Baldo & Ferreira 1999; Baldo 1999).
We used the finitetemperature BHF EOS to model PNSs in our previous papers (Nicotra et al. 2006; Burgio & Schulze 2009). The present article represents further evolution in this approach, in particular improving two aspects of the previous work: (a) an exact calculation of the finitetemperature EOS, in contrast to the socalled ``frozen correlations approximation'' used before; (b) particular attention to how to join the BHF highdensity EOS to different EOSs describing the lowdensity domain of clustered nuclear matter.
The structure of this article is as follows. The BHF approach at finite temperature is reviewed in Sect. 2, and the resulting composition and EOS of betastable matter are presented in Sect. 3. The joining to different lowdensity EOSs is discussed in Sect. 4, and the final properties of PNSs are shown in Sect. 5. Conclusions are drawn in Sect. 6.
2 Brueckner theory at finite temperature
The formalism that is closest to the BBG expansion,
and that actually reduces to it in the zerotemperature limit,
is the one formulated by Bloch & De Dominicis (1958, 1959a, b).
In this approach, the essential ingredient is the twobody scattering
matrix K, which, along with the singleparticle potential U,
satisfies the selfconsistent equations
and
where k_{i} generally denote momentum, spin, and isospin. Here V is the twobody interaction, W = E_{k1} + E_{k2} represents the starting energy, and E_{k} = k^{2}/2m + U(k) the singleparticle energy. Equation (1) coincides with the Brueckner equation for the K matrix at zero temperature if the singleparticle occupation numbers n(k) are taken at T = 0. At finite temperature, n(k) is a Fermi distribution. For a given density and temperature, Eqs. (1) and (2) have to be solved selfconsistently along with the following equation for the auxiliary chemical potential ,
The grandcanonical potential density in the BlochDe Dominicis framework can be written as the sum of a meanfield term and a correlation contribution (Baldo & Ferreira 1999; Baldo 1999),
where the trace, , is taken in the space of antisymmetrized twobody states, and the twobody scattering matrix is defined by
(5) 
If one only considers the first term in the expansion of the arctan function (Baldo & Ferreira 1999, see Appendix A), then the correlation term reduces to
which defines the grandcanonical potential in total analogy with the BBG binding potential, just using Fermi functions instead of the usual step functions at zero temperature. In this way one neglects a series of terms proportional to n(k)[1n(k)] (or powers of it), which turn out to be negligible in the temperature and density ranges relevant for neutron and protoneutron stars.
A further simplification can be achieved by disregarding the effects of finite temperature on the singleparticle potential U, and using the T=0 results in order to speed up the calculations (frozen correlations approximation). This was the procedure followed in our previous publications (Nicotra et al. 2006; Burgio & Schulze 2009). On the contrary, all calculations shown in this paper are obtained without introducing any approximation on the singleparticle potential U, and using the exact expression for the grandcanonical potential, Eq. (4).
Furthermore, in the previous work the calculations were performed using the Gmatrix instead of the Kmatrix, [i.e., taking the real part in Eq. (2) instead of Eq. (1)], which produces a more repulsive EOS at high density. Strictly speaking, the BlochDe Dominicis formalism has been derived for the Kmatrix; on the other hand, the use of the Gmatrix ensures the compatibility with the usually performed continuouschoice BHF calculations at zero temperature.
In this framework, the free energy density is then
and all other thermodynamic quantities of interest can be computed from it; namely, the ``true" chemical potential , pressure p, entropy density s, and internal energy density read as
=  (8)  
p  =  (9)  
s  =  (10)  
=  f + Ts.  (11) 
We stress that this procedure permits to fulfill the HugenholtzVan Hove theorem in the calculation of thermodynamical quantities in the Brueckner theory. For an extensive discussion of this topic, the reader is referred to Baldo (1999), and references therein.
It is well known that at zero temperature the nonrelativistic microscopic approaches do not correctly reproduce the nuclear matter saturation point , MeV. One common way of correcting this deficiency is to introduce threebody forces among nucleons. Since a complete microscopic theory of TBF is not available yet, we have adopted the phenomenological Urbana model (Carlson et al. 1983; Schiavilla et al. 1986), which consists of an attractive term due to twopion exchange with excitation of an intermediate resonance and a repulsive phenomenological central term. In the BBG approach, the TBF is reduced to a densitydependent twobody force by averaging over the position of the third particle, assuming that the probability of having two particles at a given distance is reduced according to the twobody correlation function. The corresponding EOS reproduces the nuclear matter saturation point correctly (Baldo et al. 1997; Zhou et al. 2004; Li et al. 2008a) and fulfills several requirements from the nuclear phenomenology. In all calculations presented in this paper we used the Argonne V_{18} nucleonnucleon potential (Wiringa et al. 1995) along with the phenomenological Urbana TBF.
Figure 1: Free energy per nucleon for symmetric ( upper panel) and pure neutron ( lower panel) matter (solid lines). The temperatures vary from 0 to 50 MeV in steps of 10 MeV. See text for details. 

Open with DEXTER 
Figure 2: The entropy per particle for symmetric ( upper panel) and pure neutron ( lower panel) matter at different temperatures. 

Open with DEXTER 
In Fig. 1 we display the free energy vs. the nucleon density, obtained following the procedure discussed above, both for symmetric and pure neutron matter, and for several values of the temperature between 0 and 50 MeV. The dots represent the numerical results, and the curves are empirical fits, which can be easily implemented in numerical simulations. The thin dashed curves represent the calculations obtained with the Gmatrix and the frozen correlations approximation for the singleparticle potential U, as adopted in our previous papers (Nicotra et al. 2006; Burgio & Schulze 2009). The discrepancy at high density is mainly due to the use of the Gmatrix instead of the Kmatrix, whereas the differences observed at low density and high temperature are consequences of the frozen correlations approximation.
We find that the following functional forms provide excellent parametrizations
of the new numerical results in the required ranges of density
(
)
and temperature (
)
for symmetric (SNM) and pure neutron matter (PNM):
where and F and are given in MeV and fm^{3}, respectively.
We notice that the free energy of symmetric matter shows a typical Van der Waals behavior (with MeV, ) and is a monotonically decreasing function of the temperature. At T=0 the free energy coincides with the total energy and the corresponding curve is just the usual nuclear matter saturation curve.
From the free energy we can calculate the entropy (in units of the Boltzmann constant)
from the thermodynamical relationship
(14) 
and this is displayed in Fig. 2 for both SNM and PNM.
It turns out that the dependence of the free energy on the proton
fraction can be approximated very well by a quadratic dependence, as at
zero temperature (Bombaci & Lombardo 1991; Baldo et al. 1998, 2000a):
where the symmetry energy can be expressed in terms of the difference of the energy per particle between pure neutron (x=0) and symmetric (x=0.5) matter:
Therefore, in order to treat the betaequilibrated case, it is only necessary to use the parametrizations Eqs. (12, 13) of the free energy for SNM and PNM.
Figure 3: The free symmetry energy as a function of the nucleon density at different temperatures. 

Open with DEXTER 
For completeness, we display in Fig. 3 the behavior of the free symmetry energy for several values of the temperature. Full dots represent the numerical values, and the solid curves are the results of the polynomial fits. Even at finite temperature, the symmetry energy increases nearly linearly with density, as is generally the case in the BHF approach (Li et al. 2006, 2008a). However, at low density there are significant deviations due to finite temperature, which might have experimental implications. The symmetry energy plays a significant role not only in the cooling of protoneutron stars, but also in the dynamics of heavy ion collisions induced by radioactive beams and the structure of exotic nuclei. For a review, see Li et al. (2008b).
3 Composition and EOS of hot stellar matter
Figure 4: Relative populations as functions of baryon density in neutrinofree ( upper panels), and neutrinotrapped ( lower panels) betaequilibrated matter at entropies S/A = 0, 1, 2. 

Open with DEXTER 
For stellar matter containing nucleons and leptons as relevant degrees of freedom, the chemical equilibrium conditions read explicitly as
At given baryon density , these equations have to be solved with the charge neutrality condition
and those expressing conservation of lepton numbers
As discussed in Sect. 1, we fix the lepton fractions to and for neutrinotrapped matter. When the neutrinos have left the system, their partial densities and chemical potentials vanish and the above equations simplify accordingly.
The nucleon chemical potentials are obtained from the
free energy density f, Eq. (7),
and the chemical potentials of the noninteracting leptons are obtained by solving numerically the free Fermi gas model at finite temperature. Once the hadronic and leptonic chemical potentials are known, one can proceed to calculate the composition of the stable stellar matter, and then the total pressure p through the usual thermodynamical relation
(23) 
Following this procedure, we first discuss the populations of betaequilibrated stellar matter. In Fig. 4 we display the relative particle fractions as a function of the baryon density for cold (left panels) and hot betastable matter characterized by entropy values S/A=1 (central panels) and S/A=2 (right panels). The upper panels show the particle fractions when stellar matter contains only neutrons, protons, electrons, and muons, whereas the lower panels show the particle fractions in neutrinotrapped matter. We notice that the electron fraction is greater in neutrinotrapped matter than in the neutrinofree case; therefore, the nuclear matter is more symmetric, and the resulting EOS will be softer. The appearance of muons in trapped matter is shifted to higher density values, because their onset is determined by the difference between the neutron and proton chemical potentials, Eq. (19).
We observe that thermal effects influence the populations mainly in the lowdensity region. In fact, the presence of tails in the Fermi distribution makes it possible to create (anti)particles at any density and thus typical production thresholds, like for muon creation, disappear at finite entropy. At high density, thermal effects are less important and do not change the composition appreciably when increasing entropy. As already found in (Nicotra et al. 2006) in the isothermal case and in (Burgio & Schulze 2009), the proton fraction increases slightly at high density with increasing entropy, thus causing a softening of the EOS. This is at variance with the results shown in (Prakash et al. 1997), and is presumably caused by the different temperature dependence of the potential part of the EOS.
Figure 5: Pressure as a function of energy density for betaequilibrated cold and hot matter with and without neutrino trapping. 

Open with DEXTER 
In Fig. 5, we display the pressure for betastable asymmetric matter with and without neutrinos, as a function of the energy density at entropy S/A=2 (red curves), in comparison with cold betastable matter (black curves). As discussed above, we notice that both thermal effects and neutrino trapping produce a softer EOS than in the cold case.
4 Lowdensity EOS
The Brueckner approach provides a realistic modeling of nuclear matter only at densities above about half normal nuclear matter density. Below this threshold, clusterization sets in and the system becomes inhomogeneous. Another theoretical approach therefore has to be used in this ``lowdensity'' regime, and we employ two widely used liquiddroptype models at finite temperature, namely the one of Lattimer & Swesty (1997) and the one of Shen et al. (1998a, b). In the first case we in fact compare a model with a (too) large compressibility, K=370 MeV, here denoted by ``LS'' and a second one with a lower compressibility, K=263 MeV, denoted by ``SKa''. The Shen EOS is characterized by a compressibility and a symmetry energy at saturation of , while the BHF EOS has and .
Of course, since no phase transition is involved, but only two different theoretical descriptions of the same state of matter, the joining of the two EOSs requires the thermodynamical observables, i.e., free energy, internal energy, pressure, and chemical potentials, to be continuous functions of the baryon density. Of particular interest are therefore (i) the consistent joining of the low and highdensity EOS, and (ii) eventual differences caused by the use of different lowdensity models. These aspects are studied in detail in the following.
Figure 6: The lepton fractions in betastable matter as function of baryon density for the different EOSs. 

Open with DEXTER 
Another important feature of the lowdensity domain is the treatment of neutrino trapping. Physically, neutrinos escape rapidly from the lowdensity matter during the PNS evolution, and so the trapping condition does not apply anymore. This effect can be roughly modeled by the concept of a neutrino sphere inside which the neutrinos are trapped. Typical modeldependent values for the location of the neutrino sphere found in the literature are (Gondek et al. 1997), (Strobel et al. 1999), and (Fischer et al. 2009). Given these variations, we choose the following ``natural cutoff'' procedure: by imposing the condition at any density, at a certain threshold density , the electron fraction becomes 0.4, and the neutrinos disappear naturally. For lower densities we consider the matter untrapped. This is illustrated in Fig. 6, showing the electron and neutrino fractions following this recipe with the different EOSs. This simple procedure avoids making assumptions about the neutrino sphere, but clearly a satisfactory treatment of neutrino trapping would require coupled dynamical calculations of PNS evolution and the temperaturedependent EOS, which might be possible in the future.
Figure 7: The EOS for cold untrapped matter (black curves) and for neutrinotrapped matter with entropy S/A=1 (red curves) and S/A=2 (blue curves). 

Open with DEXTER 
In Fig. 7, we display the joint lowdensity LS, SKa, Shen + BHF EOSs over a wide range of densities. The black lines represent the calculations performed at T=0, whereas the red (blue) lines are the ones for fixed entropy in neutrinotrapped matter. The transition densities between lowdensity inhomogeneous regime and homogeneous BHF regime in the different configurations are displayed in Fig. 7 as vertical bars, and summarized in Table 1. One observes in general a reduction of the transition density with decreasing compressibility of the lowdensity EOS and with increasing temperature. In all configurations, a consistent joining of the two segments of the EOS is possible, fulfilling the abovementioned continuity requirements. In practice we performed a Maxwell construction by equating pressure and chemical potentials of the low and highdensity sectors, and verified that the other thermodynamic variables do not exhibit significant discontinuities at the transition point.
As a further general prominent feature, we notice a strong stiffening of the EOS in neutrinotrapped matter as compared to the cold untrapped case in the domain (black vs. colored curves in Fig. 7), which comes from the greater lepton Fermi pressure in trapped matter (larger electron fraction and additional neutrino contribution, cf., Fig. 6).
Table 1: Transition densities (in fm^{3}) between lowdensity (LS, SKa, Shen) and BHF EOS.
5 (Proto)Neutron star configurations
The stable configurations of a (P)NS can be obtained from the
wellknown hydrostatic equilibrium equations
of Tolman, Oppenheimer, and Volkov (Shapiro & Teukolsky 1983)
for pressure p(r) and enclosed mass m(r),
once the EOS is specified, with the total internal energy density (G is the gravitational constant). For a given central value of the energy density, the numerical integration of Eqs. (24) and (25) provides the massradius relation.
Figure 8: The gravitational mass as a function of the radius (left panels), and the central energy density (right panels), obtained with the LS (top panels), SKa (central panels), and Shen (bottom panels) EOSs. 

Open with DEXTER 
Our results, using the different EOSs introduced in the previous section, are displayed in Fig. 8, which shows the gravitational mass (in units of the solar mass g) as a function of the radius (left panels), and the central density (right panels) (normalized with respect to the saturation value ). The black curves represent the calculations performed at T=0, whereas the thick (thin) colored lines denote stable configurations in neutrinofree (neutrinotrapped) matter at constant entropy S/A=1 (dashed lines) and S/A=2 (dotdashed lines).
We notice that the value of the maximum mass decreases slightly in neutrinotrapped matter, because of the softening of the EOS discussed before. The value of the maximum mass turns out to be practically independent of the entropy, whereas the value of the radius for a fixed gravitational mass depends strongly on the entropy and on the presence of neutrinos. The dots mark the values of the minimum mass and the corresponding radius a PNS may have, with typical values in the range 0.60.77 and about 4050 km, according to the chosen lowdensity EOS. The proper values are summarized in Table 2. One notes that the minimum mass configurations only probe the lowdensity part of the EOS. As expected, the minimum mass increases with decreasing compressibility of the EOS, and with increasing entropy of the matter.
Our values of the minimum mass are very similar to the ones found in Gondek et al. (1997), where the LS EOS was used over the full density range of the PNS. On the other hand, our values turn out to be lower than those calculated in Strobel et al. (1999), the reason being that our minimum mass is calculated during stage (ii) in the PNS evolution, when neutrinos are still trapped, and the core and the mantle are characterized by the same entropy values as described in the introduction.
Table 2: Properties of (P)NS minimum and maximum mass configurations.
6 Conclusions
We have presented a microscopic calculation of the EOS of hot nuclear matter in the BHF approach and provided convenient parametrizations of the free energy as a function of baryon density, proton fraction, and temperature. All other thermodynamic quantities of interest can be derived from these. This EOS was joined to different standard EOSs describing the lowdensity, inhomogeneous domain of hot nuclear matter, and we roughly modeled the structure of PNSs, using strongly idealized temperature and neutrinotrapping profiles.
Altogether we find only small variations in our results for the minimum mass configurations when different lowdensity EOSs are employed, while maximum masses are practically unaffected. Both neutrino trapping and (to a minor degree) finite temperature decrease the asymmetry of nuclear matter, and thus soften the EOS and decrease the maximum mass.
However, to provide more quantitative results, it will be necessary to perform consistent dynamical simulations of the PNS evolution, using selfconsistent temperature and neutrinotrapping profiles, along with the temperature dependent EOS. This is a difficult task for the future.
AcknowledgementsWe thank M. Baldo, V. Ferrari, and L. Gualtieri for useful discussions. We acknowledge the support of COMPSTAR, a research and training program of the European Science Foundation.
References
 Baldo, M. 1999, Nuclear Methods and the Nuclear Equation of State (Singapore: World Scientific) [Google Scholar]
 Baldo, M., & Ferreira, L. S. 1999, Phys. Rev. C, 59, 682 [NASA ADS] [CrossRef] [Google Scholar]
 Baldo, M., Bombaci, I., & Burgio, G. F. 1997, A&A, 328, 274 [NASA ADS] [Google Scholar]
 Baldo, M., Burgio, G. F., & Schulze, H.J. 1998, Phys. Rev. C, 58, 3688 [NASA ADS] [CrossRef] [Google Scholar]
 Baldo, M., Burgio, G. F., & Schulze, H.J. 2000a, Phys. Rev. C, 61, 055801 [NASA ADS] [CrossRef] [Google Scholar]
 Baldo, M., Song, H. Q., Giansiracusa, G., & Lombardo, U. 2000b, Phys. Lett. B, 473, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Baldo, M., Fiasconaro, A., Song, H. Q., Giansiracusa, G., & Lombardo, U. 2001, Phys. Rev. C, 65, 017303 [NASA ADS] [CrossRef] [Google Scholar]
 Bloch, C., & De Dominicis, C. 1958, Nucl. Phys., 7, 459 [CrossRef] [Google Scholar]
 Bloch, C., & De Dominicis, C. 1959a, Nucl. Phys., 10, 181 [CrossRef] [Google Scholar]
 Bloch, C., & De Dominicis, C. 1959b, Nucl. Phys., 10, 509 [CrossRef] [Google Scholar]
 Bombaci, I., & Lombardo, U. 1991, Phys. Rev. C, 44, 1892 [NASA ADS] [CrossRef] [Google Scholar]
 Burgio, G. F., & Schulze, H.J. 2009, Phys. Atomic Nuclei, 72, 1197 [Google Scholar]
 Burrows, A., & Lattimer, J. M. 1986, ApJ, 307, 178 [NASA ADS] [CrossRef] [Google Scholar]
 Carlson, J., Pandharipande, V. R., & Wiringa, R. B. 1983, Nucl. Phys. A, 401, 59 [NASA ADS] [CrossRef] [Google Scholar]
 Fischer, T., Whitehouse, S. C., Mezzacappa, A., Thielemann, F.K., & Liebendörfer, M. 2009, A&A, 499, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gondek, D., Haensel, P., & Zdunik, J. L. 1997, A&A, 325, 217 [NASA ADS] [Google Scholar]
 Goussard, J. O., Haensel, P., & Zdunik, J. L. 1997, A&A, 321, 822 [NASA ADS] [Google Scholar]
 Lattimer, J. M., & Swesty, F. D. 1997, Nucl. Phys. A, 535, 331, http://www.astro.sunysb.edu/lattimer/EOS/main.html [Google Scholar]
 Li, Z. H., Lombardo, U., Schulze, H.J., et al. 2006, Phys. Rev. C, 74, 047304 [NASA ADS] [CrossRef] [Google Scholar]
 Li, Z. H., Lombardo, U., Schulze, H.J., & Zuo, W. 2008a, Phys. Rev. C, 77, 034316 [NASA ADS] [CrossRef] [Google Scholar]
 Li, B. A., Chen, L. W., & Ko, C. M. 2008b, Phys. Rep., 464, 113 [NASA ADS] [CrossRef] [Google Scholar]
 Liebendörfer, M., Rampp, M., Janka, H.Th., & Mezzacappa, A. 2005, ApJ, 620, 840 [NASA ADS] [CrossRef] [Google Scholar]
 Nicotra, O. E., Baldo, M., Burgio, G. F., & Schulze, H.J. 2006, A&A, 451, 213 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pons, J. A., Reddy, S., Prakash, M., Lattimer, J. M., & Miralles, J. A. 1999, ApJ, 513, 780 [NASA ADS] [CrossRef] [Google Scholar]
 Prakash, M., Bombaci, I., Prakash, M., et al. 1997, Phys. Rep., 280, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Schiavilla, R., Pandharipande, V. R., & Wiringa, R. B. 1986, Nucl. Phys. A, 449, 219 [NASA ADS] [CrossRef] [Google Scholar]
 Schulze, H.J., Lejeune, A., Cugnon, J., Baldo, M., & Lombardo, U. 1995, Phys. Lett. B, 355, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Schulze, H.J., Lejeune, A., Cugnon, J., Baldo, M., & Lombardo, U. 1998, Phys. Rev. C, 57, 704 [NASA ADS] [CrossRef] [Google Scholar]
 Shapiro, S. L., & Teukolsky, S. A. 1983, Black Holes, White Dwarfs, and Neutron Stars (New York: John Wiley & Sons) [Google Scholar]
 Shen, H., Toki, H., Oyamatsu, K., & Sumiyoshi, K. 1998a, Nucl. Phys. A, 637, 435 [NASA ADS] [CrossRef] [Google Scholar]
 Shen, H., Toki, H., Oyamatsu, K., & Sumiyoshi, K. 1998b, Prog. Theor. Phys., 100, 1013, http://user.numazuct.ac.jp/ sumi/eos/index.html [Google Scholar]
 Song, H. Q., Baldo, M., Giansiracusa, G., & Lombardo, U. 1998, Phys. Rev. Lett., 81, 1584 [NASA ADS] [CrossRef] [Google Scholar]
 Strobel, K., Schaab, C., & Weigel, M. K. 1999, A&A, 350, 497 [NASA ADS] [Google Scholar]
 Villain, L., Pons, J. A., CerdáDurán, P., & Gourgoulhon, E. 2004, A&A, 418, 283 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wiringa, R. B., Stoks, V. G. J., & Schiavilla, R. 1995, Phys. Rev. C, 51, 38 [NASA ADS] [CrossRef] [Google Scholar]
 Zhou, X. R., Burgio, G. F., Lombardo, U., Schulze, H.J., & Zuo, W. 2004, Phys. Rev. C, 69, 018801 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: Transition densities (in fm^{3}) between lowdensity (LS, SKa, Shen) and BHF EOS.
Table 2: Properties of (P)NS minimum and maximum mass configurations.
All Figures
Figure 1: Free energy per nucleon for symmetric ( upper panel) and pure neutron ( lower panel) matter (solid lines). The temperatures vary from 0 to 50 MeV in steps of 10 MeV. See text for details. 

Open with DEXTER  
In the text 
Figure 2: The entropy per particle for symmetric ( upper panel) and pure neutron ( lower panel) matter at different temperatures. 

Open with DEXTER  
In the text 
Figure 3: The free symmetry energy as a function of the nucleon density at different temperatures. 

Open with DEXTER  
In the text 
Figure 4: Relative populations as functions of baryon density in neutrinofree ( upper panels), and neutrinotrapped ( lower panels) betaequilibrated matter at entropies S/A = 0, 1, 2. 

Open with DEXTER  
In the text 
Figure 5: Pressure as a function of energy density for betaequilibrated cold and hot matter with and without neutrino trapping. 

Open with DEXTER  
In the text 
Figure 6: The lepton fractions in betastable matter as function of baryon density for the different EOSs. 

Open with DEXTER  
In the text 
Figure 7: The EOS for cold untrapped matter (black curves) and for neutrinotrapped matter with entropy S/A=1 (red curves) and S/A=2 (blue curves). 

Open with DEXTER  
In the text 
Figure 8: The gravitational mass as a function of the radius (left panels), and the central energy density (right panels), obtained with the LS (top panels), SKa (central panels), and Shen (bottom panels) EOSs. 

Open with DEXTER  
In the text 
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.