Issue 
A&A
Volume 535, November 2011



Article Number  A37  
Number of page(s)  13  
Section  Astrophysical processes  
DOI  https://doi.org/10.1051/00046361/201117225  
Published online  28 October 2011 
The equation of state and composition of hot, dense matter in corecollapse supernovae
^{1}
MaxPlanckInstitut für Astrophysik, KarlSchwarzschildStrasse 1, Postfach 1523, 85740 Garching, Germany
email seb@mpagarching.mpg.de
^{2}
Institute for Theoretical and Experimental Physics, B. Cheremushkinskaya St. 25, 117218 Moscow, Russia
^{3}
IPMU, University of Tokyo 515 Kashiwanoha, Kashiwa, Chiba 2778568, Japan
^{4}
Computer Science Department, TechnionIsrael Institute of Technology, 32000 Haifa, Israel
^{5}
Numazu College of Technology, Ooka 3600, Numazu, Shizuoka 4108501, Japan
^{6}
Theory Center, High Energy Accelerator Research Organization (KEK), Oho 11, Tsukuba, Ibaraki 3050801, Japan
Received: 24 April 2009
Accepted: 16 August 2011
The equation of state (EOS) and composition of matter are calculated for conditions typical for precollapse and early collapse stages in corecollapse supernovae. The composition is evaluated under the assumption of nuclear statistical equilibrium, when the matter is considered as an “almost” ideal gas with corrections owing to thermal excitations of nuclei, to free nucleon degeneracy, and to Coulomb and surfaceenergy corrections. The account of these corrections allows us to obtain the composition for densities that are slightly below the nuclear matter density. Through comparisons with the EOS developed by Shen et al., which is used in most of recent supernova simulations, we examine the differences of our EOS in the multicomposition with the EOS in the approximation of one representative nucleus. We find that widely distributed compositions in the nuclear chart are different because of the different mass formulae we used, while the thermodynamical quantities are quite close to those in Shen’s EOS.
Key words: nuclear reactions, nucleosynthesis, abundances / stars: neutron / supernovae: general
© ESO, 2011
1. Introduction
A proper description of the equation of state (EOS) for subnuclear and supernuclear densities is of vital importance for current studies on the explosion mechanism of corecollapse supernovae (Bethe 1990; Suzuki 1994; Janka et al. 2007). The success of the prompt shock propagation in the supernova mechanism depends on the size of the homologously contracting inner part of a collapsing core. The larger Y_{e} (the number of electrons per nucleon) at the time of bounce, the smaller the part of matter that will be dissociated and the stronger the prompt shock wave. Therefore, all factors that influence Y_{e} must be studied thoroughly.
One factor is the mass fraction of the free protons, X_{p}, in the dense matter of the supernova core. Because the rate of electron capture on a free proton is much higher than that on a nucleus, free protons play an important role in establishing the value of Y_{e} at the neutrino trapping. The abundance of free protons depends sensitively on the nuclear models of dense matter. The variations of X_{p} amount to more than an order of magnitude in the studies of supernova EOS so far (Cooperstein 1985; Cooperstein & Baron 1990; Hillebrandt & Wolf 1985; Hillebrandt 1991; Lattimer et al. 1985; Lattimer & Swesty 1991; Shen et al. 1998a,b) and can affect the initial strength of shock wave (Bruenn 1989a; Swesty et al. 1994; Sumiyoshi et al. 2004, 2005). As shown by Bruenn (1989a,b), the free proton fraction may differ considerably depending on the nuclear interaction and its model. Another factor is the composition of various nuclei, which exist under the nuclear statistical equilibrium (NSE). The total rate of electron captures on nuclei can dominate that on free protons, if the number of the free protons is small (Hix et al. 2003), and nuclei may affect the dynamics in this way as well. Therefore, it is crucial to evaluate the composition of dense matter in a precise manner.
Another astrophysical problem where an accurate chemical composition is important is the nucleosynthesis of heavy elements. The studies of various models of the explosive nucleosynthesis as well as the rapid neutron capture process (see, e.g. Ptitsyn & Chechetkin 1982; Woosley & Hoffman 1992; Nadyozhin et al. 1998; Arnould et al. 2007, and references therein) require knowledge of the chemical composition close to the NSE as an initial condition for nuclear reaction network calculations. Neutron star mergers and collapsar models for gammaray bursts may involve similar conditions as well.
The thermodynamic properties of hot, dense matter have been investigated in various approaches: Sahalike equations (Mazurek et al. 1979; El Eid & Hillebrandt 1980; Murphy 1980; Ishizuka et al. 2003; Nadyozhin & Yudin 2004), HartreeFock approaches (Bonche & Vautherin 1981; Wolff 1983; Hillebrandt & Wolf 1985; Lassaut et al. 1987), compressible liquiddrop model (Baym et al. 1971; Lattimer et al. 1985; Lattimer & Swesty 1991), and the relativistic mean field theory (Sutaria et al. 1999) with local density approximation (Shen et al. 1998a,b). The review of those approaches, beginning from Bethe et al. (1970), can be found in the context of cold neutron star matter in Rüster et al. (2006) and Haensel et al. (2007). In most of the sets of EOS used for recent supernova simulations, the approximation of one species of nuclei (in addition to neutrons, protons and alpha particles) has been adopted (Hillebrandt & Wolf 1985; Lattimer & Swesty 1991; Shen et al. 1998b). However, the advance in recent studies of electron capture rates on nuclei for a wide mass range (Langanke & MartínezPinedo 2003) necessitates the evaluation of multicomposition of nuclei to determine the total electron capture rate. Although the treatment with one nuclear species should be a good approximation to derive overall thermodynamical properties of dense matter (Burrows & Lattimer 1984), it is necessary to take into account the multicomposition at high densities ~ 10^{12} g/cm^{3} near the neutrino trapping regime, where nuclear interaction and multicompositional treatment may become more influential in determining the EOS. The description of the EOS with the multicomposition is, of course, crucial to predict the abundance of nuclei and free protons, see e.g. Murphy (1980).
In spite of an extensive use of the Saha equation in the literature, not all relevant physics was included in the published work within this approach from the very beginning. For example, Murphy (1980) has not taken into account the effects of nonideal nucleon gas. El Eid & Hillebrandt (1980) have consistently taken into account finite temperature effects in nucleon interactions following El Eid & Hilf (1977), but they have omitted the Coulomb corrections. Later Hillebrandt & Müller (1981), Bravo & GarcíaSenz (1999) and Nadyozhin & Yudin (2005) included them in different approximations and solved the implicit set of Saha equations. Moreover, Hillebrandt & Müller (1981) have taken into account an “excludedvolume” effect owing to the finite size of nuclei. It was noted by Hillebrandt (1991) that the number of nuclear species included in the Saha treatment strongly affects the mass fraction of free protons, which may change the effective adiabatic index through electron captures. New papers have appeared after the first submission of our manuscript which discuss many points touched upon below (Hempel & SchaffnerBielich 2010; Raduta & Gulminelli 2010; Odrzywolek 2010).
In this paper we describe our code, which is suitable for studying the properties of dense matter relevant for corecollapsing supernova conditions with an emphasis on the composition of nuclei. Our code is not new: we have already compared our NSE results with the kinetic approach in the paper (Panov et al. 2001), but the code has not been described in detail. Another motivation for our work is the availability of new mass formulae as in Koura et al. (2005) and Möller et al. (1995) and detailed EOS tables such as those given in Shen et al. (1998b). Our aim is to develop and to test a practical and reliable tool for predictions of NSE composition at subnuclear densities.
Following most closely the conventional Saha approach (Clifford & Tayler 1965) and the method by Mazurek et al. (1979), we extend this approach in the following points:

At a certain level of approximation we take into account theinfluence of free nucleon gas on the surface and Coulombenergies of nuclei. We retain some terms that were omitted byMazurek et al. (1979). However,we do not take into account the “excludedvolume” effect,because we do not pretend to reach very high densities.

We have an option to include various results for nuclear partition functions like those of Fowler et al. (1978) and more recent partition functions by Engelbrecht & Engelbrecht (1991) and Engelbrecht et al. (1990).

Our network is considerably more extensive than the ones used previously. The atomic mass table is updated using recent theoretical compilations of atomic masses. It covers ~20 000 nuclides (Koura 2007, priv. comm.) for the KTUY mass formula (Koura et al. 2005) and ~9000 nuclides for the FRDM mass formula (Möller et al. 1995) as an extra option.
To examine the basic properties of the EOS for supernovae in the current formulation, we calculate the properties of dense matter covering a wide range of density, ρ, electron fraction, Y_{e}, and temperature, T. We report here the results at the equilibrium for the nuclear and electromagnetic processes, i.e. the NSE, for fixed values of Y_{e} without imposing beta equilibrium. The beta equilibrium is easily incorporated in the Saha approach as in Mazurek et al. (1979) and El Eid & Hillebrandt (1980), and our code has this option. A brief discussion of the beta equilibrium is also given below. We show the general features of the dense matter (thermodynamical quantities and compositions) for the supernova environment. To assess the dependence on the nuclear mass formulae, we compare several choices. We also examine the similarity to and difference from the single nuclear species treatment by comparisons with the Shen EOS table (Shen et al. 1998b).
We explain the formulation of NSE in Sect. 2.1 with a detailed description of the spin of nuclei in Sect. 2.2. We briefly describe the atomic mass data used for the calculations in Sect. 2.3, the Shen equation of state in Sect. 2.4 and the beta equilibrium in Sect. 2.5. We show the properties of dense matter in the current framework in Sect. 3. The summary and discussion is given in Sect. 4.
2. Thermodynamics of interacting nuclides
2.1. Inclusion of nuclear and Coulomb contributions
The basis is the equation of NSE with respect to strong and electromagnetic processes for the chemical potentials of nuclei: (1)We denote i = 1 = p for proton and i = 2 = n for neutron. We use the index i mostly for nuclei (A_{i},Z_{i}) with A > 1, but we include nucleons (i = 1,2) in the summation over all species. To find the correct expressions for μ_{i} we suppose that the free energy of the set of { N_{i} } nuclides is (2)Here Φ_{i} is the free energy per nucleus i and it may depend not only on the concentration (3)of the ith nucleus, but also on concentrations of some other nuclei n_{k}. Then we have by definition (4)Thus we find for the pressure (5)which is consistent with the definition of the grand thermodynamic potential, PV. We recall these elementary details here because thereare different prescriptions in separating terms in F_{nuc} in the literature.
For example, let us start with the Coulomb contribution to F_{nuc}. As usual we introduce (6)with a_{Zi} being the radius of the ion sphere: (7)If Q_{i} denotes the Coulomb correction to the free energy and has the functional form (8)then we find the Coulomb contribution to the chemical potential as (9)See Sect. 3.4 below for a detailed discussion of nuances in the definition of .
We closely follow Mazurek et al. (1979) in the assumed form of the free energy, but with some modifications. The main results presented below use the Coulomb correction as in Eq. (8), as did Mazurek et al. (1979). In the end we check the effect of the additional term that appears in Eq. (31) in Sect. 3.4 below. Thus, we put (10)(11)Here and are the FermiDirac expressions for the free energy of nonrelativistic nucleons and they are expressed through the Fermi integrals of halfinteger index. (The superscript zero indicates that an expression is for noninteracting particles). To calculate the Fermi integrals, we use the code by Nadyozhin (1974a,b) – for an English version see Blinnikov et al. (1996). For the bulk interaction energy of free nucleons we take the zero temperature expression from Mazurek et al. (1979), where they have used the expression of Baym et al. (1971) with minor corrections from Mackie (1976): (12)with k^{3} = 1.5π^{2}(n_{p} + n_{n}) and β = n_{p}/(n_{p} + n_{n}). We used the same values of parameters entering in Eq. (12) as in formulae (A1–A5) from Mazurek et al. (1979): k_{0} = 1.34fm^{1}; W_{0} = 15.5 MeV; s = 27.1 MeV; K = 268 MeV.
For nuclei (i > 2) we assume (13)The binding energy in Eqs. (10), (11) and (13) is (14)with M_{i} being the mass of a nucleus number i, which depends on the choice of the reference nucleus m. We use ^{56}Fe as the reference nucleus in our calculations. In calculating the total energy we exclude electron mass m_{e} from the electron energy, then we have (15)where ΔM_{i} is the atomic mass excess, which accounts for the mass of electrons. To compare with Lattimer & Swesty (1991), one has to take into account that they assume B_{n} = 0, where m = 2 = n in our notation, and B_{i} = ΔM_{i} − A_{i}ΔM_{n}.
Considering the Coulomb corrections we follow Lattimer et al. (1985) and take for the total Coulomb energy of a nucleus (16)Here r_{A} is the nuclear radius, r_{A} = 1.2 × 10^{13}A^{1/3} cm. Accordingly the Coulomb correction is (17)In the leading term, which is what we need at high density, this is just − 0.9Γ_{i} with Γ_{i} defined in Eq. (6) when the form (8) is used.
The Boltzmann expression for the chemical potential satisfies the relation (18)One has to introduce the corrections Δ_{i} in addition to : (19)where (20)The corrections for degeneracy are (21)Here μ^{0} = k_{B}TΦ^{0} with Φ^{0} from Eqs. (10), (11), computed via Fermi integrals. For the Coulomb part we have a correction: (22)For nuclear interactions we take Eq. (26) from Mazurek et al. (1979) which is (23)Here B^{0} is the binding energy in vacuum and . This is a direct consequence of the definitions for NSE (1) and for the chemical potential (4). Mazurek et al. (1979) have argued that this long expression can be simplified: “The terms in the last brackets represent the effects of the changes in binding energies of nuclei due to the interactions with the free nucleons... For simplicity we neglect this term in our numerical calculations below”. Contrary to them we do not simplify the relation Eq. (23) and take it in full, which requires an additional loop of NSE iterations.
We assume the following correction to the surface energy: (24)The expression actually used in the mass formula for the surface energy may be much more complicated than just W_{sur}A^{2/3}, but because the factor in braces tends to zero for vanishing ρ this correction does not influence the vacuum values of nuclear masses.
2.2. Spin of ground state
To find the spin of ground states of exotic and experimentally unknown nuclei, we have the simple shell model with the WoodsSaxon potential for defining the scheme of levels and nuclear systematics. It is well known from experiments and calculations using uptodate versions of the shell model with residual interaction that the spin of ground states of eveneven nuclei is equal to zero. The spin of nuclear ground states of oddA nuclei is determined generally by the spin of the uncoupled nucleon (exceptions: ). That is why in the present paper we defined the groundstate spin of oddA nuclei by the spin of an uncoupled proton or neutron according to experimental data and shell model calculations with phenomenological coupling.
The spin of ground state of oddodd nuclei is defined by the spins of the uncoupled proton, j_{p}, and neutron, j_{n}, and could not be calculated precisely even on the base of the microscopic models of atomic nuclei. That is why we used the simple phenomenological Nordheim rules in the oddodd case (Nordheim 1950).
For opposite spins of proton and neutron (j_{n} = l_{n} ± 1/2,j_{p} = l_{p} ∓ 1/2) the spin of an oddodd nucleus is equal to (25)In the other case (j_{n} = l_{n} ± 1/2,j_{p} = l_{p} ± 1/2) it is quite difficult to find the correct value of the ground state spin because it can fall into the wide range (26)But it is well known from systematics that I ≈ j_{n} + j_{p} mainly. For simplicity we chose I = j_{p} + j_{n}.
The value of the ground state spin is unlikely to influence the results of the NSE calculations strongly because in our model the nuclei are mostly in the excited states. We calculated the ground state spin with another approach (Engelbrecht & Engelbrecht 1991) to check the influence. It is also based on the simple shell model level scheme, and significantly differs from ours only in the case of oddodd nuclei, for which the spin is defined as I = (j_{p} + j_{n} + 1)/2.
We find that the calculated results depend weakly on the different approaches to the groundstate spin calculation only in the region with the significant amount of “heavy” nuclei, where the role of oddodd nuclei is more important. Because the most abundant nuclei are eveneven and oddA ones under NSE conditions, one can use any simple estimation of the ground state spin of nuclei under these conditions. For T_{9} < 10 and low ρ one should prefer our model of the groundstate spin evaluation or an “exact” calculations of nuclear characteristics from first principles.
We agree with the conclusions of Liu et al. (2007) about the role of the excited states, the effect is not very large. Liu and collaborators used a reliable method of the nuclear partition function calculation based on the Fermi gas formula. However, the question deserves further investigation, especially considering the new data from Rauscher et al. (1997); Rauscher & Thielemann (2000); Rauscher (2003). The expression used in Liu et al. (2007) is related to an energydependent level density parameter with microscopic correction to a nuclear mass model (Rauscher et al. 1997; Rauscher & Thielemann 2000; Rauscher 2003). Our code can incorporate various options for the partition functions.
2.3. Nuclear masses
We adopt the table of nuclear masses (KTUY) by Koura et al. (2005) for our EOS calculations. The prediction of nuclear masses by Koura et al. (2005) is based on the mass formula composed of macroscopic and microscopic terms (Koura et al. 2000), which treat deformation, shell and evenodd energies. It covers ~9000 species of nuclei (Z ≥ 2 and N ≥ 2) in the nuclear chart. We use the extended mass table containing ~20 000 species (Koura 2007, priv. comm.) to cover the neutron and protonrich regions in the NSE. The rootmeansquare deviation from experimental masses is 666.7 keV. We examine the dependence on the nuclear mass by adopting other tables of the nuclear masses by Hilf et al. (1976). The mass data by Möller et al. (1995) will be implemented in the code. Those mass tables predict the experimental masses very well around the stability line in the nuclear chart, however, they provide different predictions in neutronrich region away from the stability. It would be interesting to examine whether these differences appear in the composition of dense matter in supernova cores.
2.4. Shen equation of state
We adopt the table of EOS by Shen et al. (1998a,b) to investigate the influence of mixture of nuclei on the equation of state. The Shen EOS is a set of thermodynamical quantities of dense matter for the wide range of density, proton fraction and temperature for supernovae. It has been widely used for numerical simulations of corecollapse supernovae (Sumiyoshi et al. 2005, 2006; Janka et al. 2005; Burrows et al. 2006) and other astrophysical phenomena (Rosswog & Liebendörfer 2003). The equation of state is calculated by the relativistic meanfield theory (Serot & Walecka 1986), which is based on the relativistic BrücknerHartreeFock theory (Brockmann & Machleidt 1990) and is constrained by the experimental data on neutronrich nuclei (Sugahara & Toki 1994). The nuclei in dense matter are described within the local density approximation assuming one species of nucleus surrounded by neutrons, protons, and alpha particles in the WignerSeitz cell. The basic behaviour of the Shen EOS in supernova core as well as the comparison with another EOS by Lattimer & Swesty (1991) can be found in Sumiyoshi et al. (2004, 2005).
Fig. 1 Isocontours of the mass fraction of heavy nuclei, X_{AZ}, (left) and free protons, X_{p}, (right) for Y_{e} = 0.316 are shown in the plane of density and temperature. 
2.5. Beta equilibrium
Although we fix the value of Y_{e} in the current calculations, we can impose the condition of beta equilibrium in the Saha approach. We add the formulation here for the applications to protoneutron stars.
The condition of beta equilibrium is governed by the relation (27)Therefore, one has to calculate chemical potentials of neutrinos and electrons in an additional loop of iterations.
For calculating lepton chemical potentials one can use simple relations derived by Nadyozhin; Nadyozhin (1974a,b; see also Blinnikov 1987; Blinnikov & Rudzsky 1988): (28)Here n is a fermion number in the unit volume, i.e., the fermion concentration, is an antifermion concentration. is a charge of the unit volume, for example, for the electrons and it is a lepton charge for neutrinos (g = 2 for electrons and g = 1 for neutrinos). It is possible to express the lepton chemical potentials μ through the radicals if and T are given. In actual calculations it is faster to do this by Newton iterations.
Fig. 2 Isocontours of the average mass number, A, of heavy nuclei, for Y_{e} = 0.316. 
3. Results
We report the numerical results in a set of selected conditions for density, ρ, electron fraction, Y_{e}, and temperature, T. We show the properties of dense matter at typical conditions (ρ = 10^{10}–10^{13} g/cm^{3}) during the core collapse, where the composition of nuclei and free protons are crucial to determine the electron capture. We remark that the present framework is generally applicable in the wide range of conditions that are necessary for corecollapse supernova simulations. However, it is limited from below by Y_{e} ~ 0.3 because of the coverage of nuclear mass models at neutronrich region. The limited range of low T ~ 2 MeV and high ρ ~ 10^{13} g/cm^{3} is because we consider the nuclear interaction and the convergence of iterations.
We show at first the general behaviour of equation of state in the current framework with the mass formula by Koura et al. (2005) in Sect. 3.1. Next, we compare the equation of state obtained by the current Sahaapproach with that of Shen EOS (Shen et al. 1998b) in Sect. 3.2. We examine the difference between the two nuclear mass models. We also compare this with the results calculated in the formulation of Mazurek et al. (1979). In Sect. 3.3 we show the composition of nuclei in the nuclear chart for supernova cores.
3.1. General behaviour
We show in Figs. 1–4 the general features of dense matter at Y_{e} = 0.316 for a wide range of density and temperature relevant to collapsing supernova cores. We also calculated the cases with Y_{e} = 0.398 and Y_{e} = 0.473, which are not shown in Figs. 1–4, but they are partly discussed below. The values of Y_{e} are selected to match the values in the original table of Shen EOS. We see very smooth variations of the calculated quantities in the range of density and temperature shown in the plot. There are some wiggles in the plots of the average of mass number, A, and proton number, Z in Figs. 2 and 3 because of the shell effects. One can recognize the appearance of nuclei at densities ~10^{10} g/cm^{3} at the lowest temperature. The appearance of nuclei occurs with the composition, which is different from that inthe Shen EOS, as we will see in the next subsection.
Fig. 3 Isocontours of the average proton number, Z, of heavy nuclei, for Y_{e} = 0.316. 
Fig. 4 Isocontours of the nuclear contribution to entropy, S_{nuc}, for Y_{e} = 0.316 in the plane of density and temperature. 
Fig. 5 Mass fraction of nuclei, X_{A}, (left) and free protons, X_{p}, (right) as a function of density for Y_{e} = 0.316 and T = 2 MeV. The solid and dashed lines denote the NSE calculation with massformulae by (Koura 2007, priv. comm.) and Hilf et al. (1976), respectively. The dashdotted lines are obtained using (Koura 2007, priv. comm.), but with a simplified treatment of relation Eq. (23) as in Mazurek et al. (1979). The dotted lines are the results from Shen et al. (1998b). 
3.2. Comparisons among models
We compare the current results with those in the Shen EOS to see the effect of multicomposition with respect to the onespecies treatment, which has been used in most of the supernova EOS’s. We also examine the dependence on the nuclear mass models by comparing the cases by Koura et al. (2005) and Hilf et al. (1976). We set here Y_{e} = 0.316 and T = 2 MeV, which roughly corresponds to the neutronrich supernova core. This T value also corresponds to the lower T boundary in the figures shown above.
Fig. 6 Mass fraction of alphaparticles, X_{α}, as a function of density for Y_{e} = 0.316 and T = 2 MeV. The notation is the same as in Fig. 5. 
We show the mass fractions of nuclei, free protons, and ^{4}He (alphaparticles) as a function of density in Figs. 5 and 6. As the density increases, we found that the light nuclei other than ^{4}He appear in the current approach. This behaviour takes place because we treat the multicomposition including light nuclei other than ^{4}He in addition to alphaparticles. In the Shen EOS, on the other hand, alphaparticles appear in a larger amount in the density region ~10^{11} g/cm^{3}, as can be seen in Fig. 6. Figure 7 elucidates the differences of composition. From ~10^{9} g/cm^{3}, the Sahatreatment provides nuclei of average mass number A ~ 10. At densities higher than ~3 × 10^{11} g/cm^{3}, the average mass number is smaller than that of representative nuclei in the Shen EOS. We note that the average mass number in the Sahatreatment does not include the contribution of alpha particles, which are treated separately in the Shen EOS.
The appearance of light clusters has been studied in microscopic approaches (for example, Röpke et al. 1982, 1983; Typel et al. 2010) and was recently applied to astrophysics. The mixture of light elements has been demonstrated in the profiles of supernova cores by O’Connor et al. (2007); Sumiyoshi & Röpke (2008); Arcones et al. (2008). Our findings in the current study agree with those recent results, but one has to remember that we have obtained them long ago, because our numerical code has been presented already in Panov et al. (2001). The appearance of light elements may influence the average reaction rate of neutrinos on nuclei, resulting in possible modifications (by factors) of the rates of scattering, emission, and absorption through nuclear structures different from nucleons or ^{4}He. It is necessary to evaluate the mean free path using the sum of individual crosssections weighted by the number density of corresponding target nuclei (see, e.g., Arcones et al. 2008; Nakamura et al. 2009; Odrzywolek 2009; Furusawa et al. 2011). Because the crosssections for the coherent scattering on nuclei is ∝ the result is quite different when only alphaparticles represent the whole set of light nuclei. For the given density the number of nuclei is ∝ A^{1}, so the extinction coefficient goes as A^{2}/A ∝ A, and the neutrino mean free path ∝ A^{1}. This has a substantial effect on the mean free path, e.g. for the conditions in Figs. 6 and 7 when in the realistic calculation the mass fraction of light nuclei with A_{i} > 4 at ρ ~ 10^{11} g/cm^{3} is larger than that of ^{4}He.
The thermodynamical quantities, i.e. entropy and pressure of nuclear contribution, are shown in Fig. 8. The entropy of the Shen EOS agrees very well with the current result. The pressure of the Shen EOS has a trough around 10^{13} g/cm^{3} because of the decrease of nuclear part of free energy. As already reported in Shen et al. (1998a), this is caused by the increasing binding energy of nuclei. The occurrence of the drop of nuclear pressure depends on the temperature and electron fraction. Because the pressure consists of contributions from the thermal part and the nuclear interactions, the nuclear pressure of the mixture of nucleons and nuclei is determined by a delicate sum. (We point out that the negative Coulomb corrections are not included into our plots for P_{nuc}, see the discussion in Sect. 3.4.) The pressure also depends on the treatment of the transition from the mixture of nuclei to the uniform matter (which we do not touch upon in the current work). This question is analyzed in a recent paper by Hempel & SchaffnerBielich (2010), who obtain a good agreement with negative P_{nuc} in Shen et al. (1998a) when they use the Maxwellconstruction in phase transition to uniform nuclear matter. They find that the use of a Gibbsconstruction would result in a strictly increasing pressure. In Fig. 22 of Raduta & Gulminelli (2010), the baryonic pressure stays positive and does not show any drop. Raduta & Gulminelli (2010) analysed in detail the difference of this property with single nucleus treatments like LattimerSwesty EOS (Lattimer & Swesty 1991). Our NSE approach is less sophisticated than that of Raduta & Gulminelli (2010), but of course it is much closer to them than to Lattimer & Swesty (1991); Shen et al. (1998b). We note that one should add lepton and radiation contributions, therefore, the total pressure is always positive, even considering negative Coulomb corrections (see also Fig. 23 in Raduta & Gulminelli 2010).
In Fig. 9 we show the average mass number and proton number of nuclei for Y_{e} = 0.473 and T=0.63 MeV, which is closer to the condition at the centre of the initial progenitors. We can see that the mass number (proton number also) in the current treatment islarger at low density and smaller at high density. This is quite similar to the more neutronrich and higher temperature case of Y_{e} = 0.316 and T = 2 MeV in Fig. 7.
Now let us comment on the difference in the calculated EOS if we use the treatment of by Eq. (23) suggested by Mazurek et al. (1979). The results of calculations made without the second line in Eq. (23), which expresses the dependence of binding energy of nuclei on the nucleon fractions, are shown by the dashdotted lines in Figs. 5 to 8. In general, the treatment without iterations gives similar results to the one with the iterations to calculate consistently the modified binding energy owing to the presence of nucleons at this condition. In Fig. 10 we show the mass fraction of free protons as a function of density for Y_{e} = 0.398 and T = 2 MeV. We found that the approximation used by Mazurek et al. (1979) for Eq. (23) gives similar results to those of the Shen EOS and our treatment yields a cosiderably lower proton fraction. This may influence the electron capture rate in collapsing supernova cores. It is worth exploring the reasons for the difference and its dependence on Y_{e}, but this is better conducted with a more modern approach to NSE as in Hempel & SchaffnerBielich (2010); Raduta & Gulminelli (2010), which is beyond the scope of the current paper.
We also explored the dependence of the bulk thermodynamics on the nuclear mass models. The calculations by the two mass models used by us are similar to each other and there is hardly any difference. Only in the prediction of mass number and proton number,the models differ slightly. In this respect the derivation of bulk quantities such as entropy and pressure, can safely be calculated irrespective of the uncertainties of mass formulae. We note, however, that the detailed composition depends on the adopted mass formula, as we will show below.
3.3. Nuclear composition
We explore the calculated composition using the mass formula by Koura et al. (2005) in Figs. 11–13 for three typical conditions of the supernova core. We took the conditions from the numerical simulation of corecollapse from a 15 M_{⊙} star by Sumiyoshi et al. (2005).
Fig. 7 Average mass number, A, (left) and average proton number, Z, (right) of nuclei as a function of density for Y_{e} = 0.316 and T = 2 MeV. The notation is the same as in Fig. 5. 
Fig. 8 Entropy, S_{nuc} (left), and pressure P_{nuc} (right), of nuclei as functions of density for Y_{e} = 0.316 and T = 2 MeV. The notation is the same as in Fig. 5. 
Fig. 9 Average mass number, A, (left) and average proton number, Z, (right) of nuclei as a function of density for Y_{e} = 0.473 and T = 0.63 MeV. The notation is the same as in Fig. 5. 
At ρ = 10^{11} g/cm^{3}, where electron captures continue during the gravitational collapse from the initial Fe core, nuclei up to A ~ 100 appear mostly with the peak abundance at Z ~ 34 and the magic number N = 50. By the time the density reaches ρ = 10^{12} g/cm^{3}, electron capture is about to cease and Y_{e} becomes smaller. Because of the higher temperature owing to the collapse, the distribution becomes wider, reaching A ~ 132 at the double magic nuclei, ^{132}Sn. Neutronrichness (small Y_{e}) makes the peak position at more neutronrich, but at N = 50. At high density of 10^{13} g/cm^{3} and high temperature of 3 MeV, which corresponds to the moment just before the core bounce, the distribution extends from nucleons to A ~ 160 continuously. We see the effect of magic numbers for neutron number 50 and 82 and evenodd numbers, which exist in nuclear mass data.
To demonstrate the difference of composition due to the mass models, we show the case of ρ = 10^{12} g/cm^{3} with the mass data by Hilf et al. (1976) in Fig. 14. Through the comparison with Fig. 12, the distribution is narrower both in N and A = N + Z directions. The peak positions are quite close to each other in the two models. We see similar differences also in the case of ρ = 10^{11} g/cm^{3} and ρ = 10^{13} g/cm^{3}, especially on the strength of magic numbers. Neutronrichness at high density causes more pronounced differences in composition among them.
3.4. Difference in Coulomb corrections
In this section we briefly discuss the different prescriptions on the form of Coulomb corrections and their influence on the composition of hot, dense matter. There are several expressions of the form of Coulomb corrections found in the literature. These expressions have been used for astrophysics as well as material science and are still under debate. Allmost all of them use a form for the Coulomb free energy like in Eq. (8). However, the expression for Γ_{i} may be understood differently. The most widely used expression is Eq. (9) based on Γ_{i} from Eq. (6). The question is how to treat the electron number density n_{e} that defines the radius of the ion sphere in (7).
The standard approach is to forget about the dependence of n_{e} in (7) on ion numbers and to treat n_{e} as an independent variable. Then we obtain the chemical potential for ions (9). The other approach substitutes from the electroneutrality into the expression for the free energy (8), which leads to another expression for the chemical potential that we will discuss below. The assumptions are different, however, their influence on supernova matter has not been studied so far.
We briefly illustrate the latter approach here in the simplest case of a fully ionized plasma with only one ion (A_{i},Z_{i}).Then we have n_{e} = Z_{i}n_{i}, and instead of (6) we may define a_{Zi} as (29)and then we write (30)and obtain (31)contrary to our Eq. (9), to Eq. (17) of Mazurek et al. (1979), and to Eq. (97) of Yakovlev & Shalybkov (1989). The expression μ_{i} = k_{B}Tf_{0}(Γ_{i}) is standard, but it leads to zero correction for P_{nuc} in Eq. (5). On the other hand, the expression Eq. (31) gives the expected correction to P_{nuc} after substitution into Eq. (5).
Basically, the formulae in Yakovlev & Shalybkov (1989) are correct, because one is free to separate the Coulomb contributions and ascribe them either to electrons or ions, but one has to avoid confusion with this separation. Haensel et al. (2007; Chabrier & Potekhin see also 1998) rely mostly on the free energy F instead of μ_{i} and avoid any confusion. Indeed, the expression that is equivalent to (31) appeared already in Dewitt et al. (1973). As we noticed after Eq. (17), we have in the leading term at high density just − 0.9Γ_{i} with Γ_{i} defined in Eq. (6) when the form (9) is used. With the modification as in Eq. (31) the correction will be − 1.2Γ_{i} – very close to the numbers in Eqs. (56) and (57) found by Dewitt et al. (1973) by two different methods. See the more detailed discussion of this difference in Glazyrin & Blinnikov (2010).
The difference between Eqs. (9) and (31) is explained by the neglect of electroneutrality in the standard Eq. (9): the latter assumes independence of n_{e} and n_{i}, which is not true in general (Glazyrin & Blinnikov 2010).
The argument on the role of the chemical potentials and on how to approach them still persists. For example, Brown et al. (2006) discuss the difference of their use of chemical potentials with that of Dewitt et al. (1973). Essentially, Brown et al. (2006) found that the Coulomb correction justified by Dewitt et al. (1973) for a classical plasma is applicable to a quantum plasma as well (see also Brydges & Martin 1999, for a detailed review on the modern perspective on Coulomb systems).
Fig. 10 Mass fraction of free protons, X_{p}, as a function of density for Y_{e} = 0.398 and T = 2 MeV. The notation is the same as in Fig. 5. 
We cannot use Eqs. (29, 31) directly in a multicomponent plasma because n_{e} depends on the sum of many ions, nevertheless, those expressions should be meaningful in many cases when the distribution of nuclei is dominated by a set of species with similar values of (A_{i},Z_{i}). We compare the effect of changing corrections as in formula (9) and in (31) on the proton fraction in Fig. 15. The effect is visible only at highest density, but it is still very small. An effect slightly more pronounced is obvious in Fig. 16 for the average mass of heavy nuclei again at highest density when the NSE model is near the border of applicability. A more detailed study is needed here along the lines undertaken recently by Nadyozhin & Yudin (2005).
Fig. 11 Mass fraction of nuclei, X_{A}, with the proton number, Z, and neutron number, N = A − Z in the nuclear chart for ρ = 10^{11}g/cm^{3},Y_{e} = 0.398 and T = 1 MeV. 
Fig. 12 Mass fraction of nuclei, X_{A}, with the proton number, Z, and neutron number, N = A − Z in the nuclear chart for ρ = 10^{12}g/cm^{3},Y_{e} = 0.316 and T = 2 MeV. 
Fig. 13 Mass fraction of nuclei, X_{A}, with the proton number, Z, and neutron number, N = A − Z in the nuclear chart for ρ = 10^{13}g/cm^{3},Y_{e} = 0.316 and T = 3 MeV. 
Fig. 14 Same as Fig. 12 but for Hilf et al. (1976). 
Fig. 15 Isocontours of the proton fraction, X_{p}, for Y_{e} = 0.473. Left panel: the Coulomb correction is taken as in Eq. (8). Right panel: as in Eq. (31). 
Fig. 16 Isocontours of the average mass number, A, of heavy nuclei, for Y_{e} = 0.473. Left panel: the Coulomb correction is taken as in Eq. (8). Right panel: as in Eq. (31). 
4. Summary and discussions
We have studied the properties of hot and dense matter at the conditions relevant to supernova cores in the extended Saha approach. The multicomposition of nuclei is taken into account by solving the chemical equilibrium among nuclei, neutrons, and protons. The contributions to chemical potentials owing to the degeneracy of nucleons, to Coulomb, nuclear, and surface effects are included in the Saha formulation. The modification of binding energies of nuclei through the existence of nucleons outside nucleus is solved selfconsistently. This is different from the treatment of Mazurek et al. (1979). We refer the reader to Nadyozhin & Yudin (2005), where the difficulty of selfconsistent treatment is pointed out. We have adopted the recent mass formula by Koura et al. (2005), covering ~20 000 nuclides, and compared the results with the ones by Hilf et al. (1976).
We have calculated extensive grids of density and temperature for several choices of Y_{e}. We found that the obtained quantities behave in general very smoothly. When we compared our results with the values in the Shen EOS, which is a popular set of EOS for supernova simulations, they coincide closely for thermodynamical quantities. This may suggest that the onespecies treatment of nuclei adopted in most of modern supernova EOS’s is a good approximation in this regard. However, the composition of hot and dense matter in the extended Saha approach is broader than just a single heavy nucleus and alpha particles in the onespecies approximation as in the Shen EOS. The appearance of nuclei starts with a mixture of various light nuclei at low density in the current approach. In the Shen EOS, ^{4}He is dominant at the same region and onespecies nuclei appear at higher density. The average mass and proton numbers of nuclei (apart from ^{4}He) are larger at low density and smaller at high density than those obtained from the Shen EOS. The former arises from the appearance of light nuclei around A ~ 10 in addition to alpha particles, which are only included in Shen EOS. We also explored the composition in the supernova core in the nuclear chart. We found that the distribution extends widely for high density and temperature and it depends on the choice of the mass formula, which provides different shell effects.
Our approach is fairly simple and may be compared with the more sophisticated modern treatment of nuclear reactions based on the statistical multifragmentation model (SMM): see Botvina & Mishustin (2010). Recently the difference between the treatments under single and multi nuclear species approximations with phenomenological mass formulae has been explored to demonstrate the overestimation of the mass number in the single nuclear species by Souza et al. (2009).
Our results assist the understanding of the gravitational collapse of massive stars for supernova explosions. The difference of composition from the onespecies treatment in the Shen EOS suggests that the proper treatment for a multicomposition of nuclei is important and may significantly affect the dynamics of core collapse. Because the electron capture on nuclei and free protons is crucial to determine the size of bounce core, i.e. the location of initial launch of shock wave, one should take the composition properly into account. These efforts have been made to evaluate the total electron capture rates on nuclei by averaging the distributions (Hix et al. 2003). One should, however, work carefully to predict the composition, because the mass and the nuclear effects may affect the detailed abundance. It is also necessary to calculate the equation of state as well as the thermodynamical quantities and compositions in a consistent manner. We are aiming to apply the extended Saha treatment to provide EOS sets for supernova simulations in the future.
Acknowledgments
A part of this work was made while S. Blinnikov, I. Panov, and K. Sumiyoshi were visiting the MaxPlanckInstitut für Astrophysik, Garching. These authors are especially thankful to W. Hillebrandt, E. Müller, and H.Th. Janka for the hospitality and to W.H. and E.M. for reading the manuscript at the early and final stages of the work. S.B. cordially thanks K. Nomoto, K. Sato, and H. Murayama for hospitality at RESCEU and IPMU in Tokyo, where he was able to complete the paper. S.B.’s work is supported by World Premier International Research Center Initiative (WPI), MEXT, Japan. We are grateful to D. Yakovlev for drawing our attention to the problem of Coulomb corrections in nuclear interactions and for his valuable consultations. The late C. Engelbrecht has given us his results on partition functions prior to publication. The work of S.B. and I.P. in Russia is partly supported by the grants RFBR 100200249, 100201398, by “Scientific School” grants 3458.2010.2, 3899.2010.2, by contract with Agency for Science and Innovation No. 02.740.11.0250, SNSF grant No. IZ73Z0128180/1 under the programme SCOPES. K.S. would like to thank C. Ishizuka, S. Furusawa, and H. Shen for helpful discussions. K.S. is partly supported by the GrantsinAid for the Scientific Research (18540291, 18540295, 19104006, 19540252, 22540296) and by the GrantsinAid for Scientific Research on Innovative Areas of the MEXT (20105004) in Japan. S.B. is grateful to Elena Blinnikova for typing many long formulae in TeX, and to Mona Frommert for her valuable help at one of the critical stages of the project.
References
 Arcones, A., MartínezPinedo, G., O’Connor, E., et al. 2008, Phys. Rev. C, 78, 015806 [NASA ADS] [CrossRef] [Google Scholar]
 Arnould, M., Goriely, S., & Takahashi, K. 2007, Phys. Rep., 450, 97 [NASA ADS] [CrossRef] [Google Scholar]
 Baym, G., Bethe, H. A., & Pethick, C. J. 1971, Nucl. Phys. A, 175, 225 [NASA ADS] [CrossRef] [Google Scholar]
 Bethe, H. A. 1990, Rev. Mod. Phys., 62, 801 [NASA ADS] [CrossRef] [Google Scholar]
 Bethe, H. A., Borner, G., & Sato, K. 1970, A&A, 7, 279 [NASA ADS] [Google Scholar]
 Blinnikov, S. I. 1987, Sov. Astron. Lett., C13, 346 [NASA ADS] [Google Scholar]
 Blinnikov, S. I., & Rudzsky, M. A. 1988, Astrofizika, 29, 385 [NASA ADS] [Google Scholar]
 Blinnikov, S. I., DuninaBarkovskaya, N. V., & Nadyozhin, D. K. 1996, ApJS, 106, 171 [NASA ADS] [CrossRef] [Google Scholar]
 Bonche, P., & Vautherin, D. 1981, Nucl. Phys. A, 372, 496 [NASA ADS] [CrossRef] [Google Scholar]
 Botvina, A., & Mishustin, I. 2010, Nucl. Phys. A, 843, 98 [NASA ADS] [CrossRef] [Google Scholar]
 Bravo, E., & GarcíaSenz, D. 1999, MNRAS, 307, 984 [NASA ADS] [CrossRef] [Google Scholar]
 Brockmann, R., & Machleidt, R. 1990, Phys. Rev. C, 42, 1965 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, L. S., Dooling, D. C., & Preston, D. L. 2006, Phys. Rev. E, 73, 056406 [NASA ADS] [CrossRef] [Google Scholar]
 Bruenn, S. W. 1989a, ApJ, 340, 955 [NASA ADS] [CrossRef] [Google Scholar]
 Bruenn, S. W. 1989b, ApJ, 341, 385 [NASA ADS] [CrossRef] [Google Scholar]
 Brydges, D. C., & Martin, P. A. 1999, J. Stat. Phys., 96, 1163 [NASA ADS] [CrossRef] [Google Scholar]
 Burrows, A., & Lattimer, J. M. 1984, ApJ, 285, 294 [Google Scholar]
 Burrows, A., Livne, E., Dessart, L., Ott, C. D., & Murphy, J. 2006, ApJ, 640, 878 [NASA ADS] [CrossRef] [Google Scholar]
 Chabrier, G., & Potekhin, A. Y. 1998, Phys. Rev. E, 58, 4941 [NASA ADS] [CrossRef] [Google Scholar]
 Clifford, F. E., & Tayler, R. J. 1965, MNRAS, 69, 21 [NASA ADS] [Google Scholar]
 Cooperstein, J. 1985, Nucl. Phys., A429, 527 [Google Scholar]
 Cooperstein, J., & Baron, E. A. 1990, in Supernovae, ed. A. Petschek (New York: SpringerVerlag), 213 [Google Scholar]
 Dewitt, H. E., Graboske, H. C., & Cooper, M. S. 1973, ApJ, 181, 439 [NASA ADS] [CrossRef] [Google Scholar]
 El Eid, M. F., & Hilf, E. R. 1977, A&A, 57, 243 [NASA ADS] [Google Scholar]
 El Eid, M. F., & Hillebrandt, W. 1980, A&AS, 42, 215 [NASA ADS] [Google Scholar]
 Engelbrecht, C. A., & Engelbrecht, J. R. 1991, Ann. Phys., 207, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Engelbrecht, C. A., Engelbrecht, J. R., & Joubert, D. P. 1990, unpublished [Google Scholar]
 Fowler, W. A., Engelbrecht, C. A., & Woosley, S. E. 1978, ApJ, 226, 984 [NASA ADS] [CrossRef] [Google Scholar]
 Furusawa, S., Yamada, S., Sumiyoshi, K., & Suzuki, H. 2011, ApJ, 738, 178 [Google Scholar]
 Glazyrin, S. I., & Blinnikov, S. I. 2010, J. Phys. A Math. Gen., 43, 075501 [NASA ADS] [CrossRef] [Google Scholar]
 Haensel, P., Potekhin, A. Y., & Yakovlev, D. G. 2007, Neutron stars, 1 Equation of state and structure, Astrophysics and Space Science Library (New York: SpringerVerlag), XXIII, 619 [NASA ADS] [Google Scholar]
 Hempel, M., & SchaffnerBielich, J. 2010, Nucl. Phys. A, 837, 210 [NASA ADS] [CrossRef] [Google Scholar]
 Hilf, E. R., von Groote, H., & Takahashi, K. 1976, in 3rd International Conference on Nuclei Far from Stability, Cargese, France, 19–26 May, CERN No. 76–13, 142 [Google Scholar]
 Hillebrandt, W. 1991, in HighPressure Equations of State: Theory and Applications, ed. S. Eliezer, & R. A. Ricci (Amsterdam: NorthHolland), 399 [Google Scholar]
 Hillebrandt, W., & Müller, E. 1981, A&A, 103, 147 [NASA ADS] [Google Scholar]
 Hillebrandt, W., & Wolf, R. G. 1985, in Nucleosynthesis: Challenges and New Developments, ed. W. D. Arnett, & J. W. Truran (Chicago: University of Chicago Press), 131 [Google Scholar]
 Hix, W. R., Messer, O. E. B., Mezzacappa, A., et al. 2003, Phys. Rev. Lett., 91, 201102 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Ishizuka, C., Ohnishi, A., & Sumiyoshi, K. 2003, Nucl. Phys. A, 723, 517 [NASA ADS] [CrossRef] [Google Scholar]
 Janka, H.T., Buras, R., Kitaura Joyanes, F. S., et al. 2005, Nucl. Phys. A, 758, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Janka, H.T., Langanke, K., Marek, A., MartinezPinedo, G., & Müller, B. 2007, Phys. Rep., 442, 38 [NASA ADS] [CrossRef] [Google Scholar]
 Koura, H., Uno, M., Tachibana, T., & Yamada, M. 2000, Nucl. Phys. A, 674, 47 [Google Scholar]
 Koura, H., Tachibana, T., Uno, M., & Yamada, M. 2005, Prog. Theor. Phys., 113, 305 [Google Scholar]
 Langanke, K., & MartínezPinedo, G. 2003, Rev. Mod. Phys., 75, 819 [NASA ADS] [CrossRef] [Google Scholar]
 Lassaut, M., Flocard, H., Bonche, P., Heenen, P. H., & Suraud, E. 1987, A&A, 183, L3 [NASA ADS] [Google Scholar]
 Lattimer, J. M., & Swesty, F. D. 1991, Nucl. Phys., A535, 331 [CrossRef] [Google Scholar]
 Lattimer, J. M., Pethick, C. J., Ravenhall, D. G., & Lamb, D. Q. 1985, Nucl. Phys. A, 432, 646 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, M.Q., Zhang, J., & Luo, Z.Q. 2007, Chin. Phys., 16, 3146 [NASA ADS] [CrossRef] [Google Scholar]
 Mackie, F. D. 1976, Ph.D. Thesis, University of Illinois at UrbanaChampaign [Google Scholar]
 Mazurek, T. J., Lattimer, J. M., & Brown, G. E. 1979, ApJ, 229, 713 [Google Scholar]
 Möller, P., Nix, J. R., & Swiatecki, W. J. 1995, Atomic Data Nucl. Data Tables, 59, 185 [Google Scholar]
 Murphy, M. J. 1980, APJS, 42, 385 [NASA ADS] [CrossRef] [Google Scholar]
 Nadyozhin, D. K. 1974a, Nauchnye informatsii, 32, 3, in Russian [NASA ADS] [Google Scholar]
 Nadyozhin, D. K. 1974b, Nauchnye informatsii, 33, 117, in Russian [NASA ADS] [Google Scholar]
 Nadyozhin, D. K., & Yudin, A. V. 2004, Astron. Lett., 30, 634 [Google Scholar]
 Nadyozhin, D. K., & Yudin, A. V. 2005, Astron. Lett., 31, 271 [NASA ADS] [CrossRef] [Google Scholar]
 Nadyozhin, D. K., Panov, I. V., & Blinnikov, S. I. 1998, A&A, 335, 207 [NASA ADS] [Google Scholar]
 Nakamura, S. X., Sumiyoshi, K., & Sato, T. 2009, Phys. Rev. C, 80, 035802 [NASA ADS] [CrossRef] [Google Scholar]
 Nordheim, L. W. 1950, Phys. Rev., 78, 294 [NASA ADS] [CrossRef] [Google Scholar]
 O’Connor, E., Gazit, D., Horowitz, C. J., Schwenk, A., & Barnea, N. 2007, Phys. Rev. C, 75, 055803 [NASA ADS] [CrossRef] [Google Scholar]
 Odrzywolek, A. 2009, Phys. Rev. C, 80, 045801 [NASA ADS] [CrossRef] [Google Scholar]
 Odrzywolek, A. 2010 [arXiv:1008.0798] [Google Scholar]
 Panov, I., Blinnikov, S., & Thielemann, F.K. 2001, Astron. Lett., 27, 239 [NASA ADS] [CrossRef] [Google Scholar]
 Ptitsyn, D. A., & Chechetkin, V. M. 1982, Soviet Astron. Lett., 8, 322 [NASA ADS] [Google Scholar]
 Raduta, A. R., & Gulminelli, F. 2010, Phys. Rev. C, 82, 065801 [NASA ADS] [CrossRef] [Google Scholar]
 Rauscher, T. 2003, ApJS, 147, 403 [Google Scholar]
 Rauscher, T., & Thielemann, F.K. 2000, Atomic Data Nucl. Data Tables, 75, 1 [Google Scholar]
 Rauscher, T., Thielemann, F.K., & Kratz, K.L. 1997, Phys. Rev. C, 56, 1613 [NASA ADS] [CrossRef] [Google Scholar]
 Röpke, G., Münchow, L., & Schulz, H. 1982, Nucl. Phys. A, 379, 536 [NASA ADS] [CrossRef] [Google Scholar]
 Röpke, G., Schmidt, M., Münchow, L., & Schulz, H. 1983, Nucl. Phys. A, 399, 587 [NASA ADS] [CrossRef] [Google Scholar]
 Rosswog, S., & Liebendörfer, M. 2003, MNRAS, 342, 673 [NASA ADS] [CrossRef] [Google Scholar]
 Rüster, S. B., Hempel, M., & SchaffnerBielich, J. 2006, Phys. Rev. C, 73, 035804 [NASA ADS] [CrossRef] [Google Scholar]
 Serot, B. D., & Walecka, J. D. 1986, in Advances in Nuclear Physics, ed. J. W. Negele, & E. Vogt (New York: Plenum Press), 16, 1 [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. Thoer. Phys., 100, 1013 [Google Scholar]
 Souza, S. R., Steiner, A. W., Lynch, W. G., Donangelo, R., & Famiano, M. A. 2009, ApJ, 707, 1495 [Google Scholar]
 Sugahara, Y., & Toki, H. 1994, Nucl. Phys. A, 579, 557 [NASA ADS] [CrossRef] [Google Scholar]
 Sumiyoshi, K., & Röpke, G. 2008, Phys. Rev. C, 77, 055804 [NASA ADS] [CrossRef] [Google Scholar]
 Sumiyoshi, K., Suzuki, H., Yamada, S., & Toki, H. 2004, Nucl. Phys. A, 730, 227 [NASA ADS] [CrossRef] [Google Scholar]
 Sumiyoshi, K., Yamada, S., Suzuki, H., et al. 2005, ApJ, 629, 922 [NASA ADS] [CrossRef] [Google Scholar]
 Sumiyoshi, K., Yamada, S., Suzuki, H., & Chiba, S. 2006, Phys. Rev. Lett., 97, 091101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Sutaria, F. K., Ray, A., Sheikh, J. A., & Ring, P. 1999, A&A, 349, 135 [NASA ADS] [Google Scholar]
 Suzuki, H. 1994, in Physics and Astrophysics of Neutrinos, ed. M. Fukugita, & A. Suzuki (Tokyo: SpringerVerlag), 763 [Google Scholar]
 Swesty, F. D., Lattimer, J. M., & Myra, E. S. 1994, ApJ, 425, 195 [NASA ADS] [CrossRef] [Google Scholar]
 Typel, S., Röpke, G., Klähn, T., Blaschke, D., & Wolter, H. H. 2010, Phys. Rev. C, 81, 015803 [NASA ADS] [CrossRef] [Google Scholar]
 Wolff, R. G. 1983, Ph.D. Thesis, Technische Universität, München [Google Scholar]
 Woosley, S. E., & Hoffman, R. 1992, ApJ, 395, 202 [NASA ADS] [CrossRef] [Google Scholar]
 Yakovlev, D. G., & Shalybkov, D. A. 1989, Astrophys. Space Phys. Rev., Soviet Sci. Rev. E., 7, 311 [Google Scholar]
All Figures
Fig. 1 Isocontours of the mass fraction of heavy nuclei, X_{AZ}, (left) and free protons, X_{p}, (right) for Y_{e} = 0.316 are shown in the plane of density and temperature. 

In the text 
Fig. 2 Isocontours of the average mass number, A, of heavy nuclei, for Y_{e} = 0.316. 

In the text 
Fig. 3 Isocontours of the average proton number, Z, of heavy nuclei, for Y_{e} = 0.316. 

In the text 
Fig. 4 Isocontours of the nuclear contribution to entropy, S_{nuc}, for Y_{e} = 0.316 in the plane of density and temperature. 

In the text 
Fig. 5 Mass fraction of nuclei, X_{A}, (left) and free protons, X_{p}, (right) as a function of density for Y_{e} = 0.316 and T = 2 MeV. The solid and dashed lines denote the NSE calculation with massformulae by (Koura 2007, priv. comm.) and Hilf et al. (1976), respectively. The dashdotted lines are obtained using (Koura 2007, priv. comm.), but with a simplified treatment of relation Eq. (23) as in Mazurek et al. (1979). The dotted lines are the results from Shen et al. (1998b). 

In the text 
Fig. 6 Mass fraction of alphaparticles, X_{α}, as a function of density for Y_{e} = 0.316 and T = 2 MeV. The notation is the same as in Fig. 5. 

In the text 
Fig. 7 Average mass number, A, (left) and average proton number, Z, (right) of nuclei as a function of density for Y_{e} = 0.316 and T = 2 MeV. The notation is the same as in Fig. 5. 

In the text 
Fig. 8 Entropy, S_{nuc} (left), and pressure P_{nuc} (right), of nuclei as functions of density for Y_{e} = 0.316 and T = 2 MeV. The notation is the same as in Fig. 5. 

In the text 
Fig. 9 Average mass number, A, (left) and average proton number, Z, (right) of nuclei as a function of density for Y_{e} = 0.473 and T = 0.63 MeV. The notation is the same as in Fig. 5. 

In the text 
Fig. 10 Mass fraction of free protons, X_{p}, as a function of density for Y_{e} = 0.398 and T = 2 MeV. The notation is the same as in Fig. 5. 

In the text 
Fig. 11 Mass fraction of nuclei, X_{A}, with the proton number, Z, and neutron number, N = A − Z in the nuclear chart for ρ = 10^{11}g/cm^{3},Y_{e} = 0.398 and T = 1 MeV. 

In the text 
Fig. 12 Mass fraction of nuclei, X_{A}, with the proton number, Z, and neutron number, N = A − Z in the nuclear chart for ρ = 10^{12}g/cm^{3},Y_{e} = 0.316 and T = 2 MeV. 

In the text 
Fig. 13 Mass fraction of nuclei, X_{A}, with the proton number, Z, and neutron number, N = A − Z in the nuclear chart for ρ = 10^{13}g/cm^{3},Y_{e} = 0.316 and T = 3 MeV. 

In the text 
Fig. 14 Same as Fig. 12 but for Hilf et al. (1976). 

In the text 
Fig. 15 Isocontours of the proton fraction, X_{p}, for Y_{e} = 0.473. Left panel: the Coulomb correction is taken as in Eq. (8). Right panel: as in Eq. (31). 

In the text 
Fig. 16 Isocontours of the average mass number, A, of heavy nuclei, for Y_{e} = 0.473. Left panel: the Coulomb correction is taken as in Eq. (8). Right panel: as in Eq. (31). 

In the text 
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.