We investigate the possibility of a supernova in supermassive ($5 \times 10^4 \;M_\odot$) population III stars induced by a general relativistic instability occurring in the helium burning phase. This explosion could occur via rapid helium burning during an early contraction of the isentropic core. Such an explosion would be visible to future telescopes and could disrupt the proposed direct collapse formation channel for early universe supermassive black holes. We simulate first the stellar evolution from hydrogen burning using a 1D stellar evolution code with a post Newtonian approximation; at the point of dynamical collapse, we switch to a 1D (general relativistic) hydrodynamics code with the Misner-Sharpe metric. In opposition to a previous study, we do not find an explosion in the non rotating case, although our model is close to exploding for a similar mass to the explosion in the previous study. When we include slow rotation, we find one exploding model, and we conclude that there likely exist additional exploding models, though they may be rare.

We report on the core-collapse supernova simulation we conducted for a 11.2M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at ~10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is ~0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and possibly improves the closure relation there.

Using the Boltzmann-radiation-hydrodynamics code which solves the Boltzmann equation for the neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer-Swesty (LS) and Furusawa-Shen (FS). We extend the simulation time of the LS model and conducted thorough investigations, though the previous paper reported some of the results briefly. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by the photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The proto-neutron star seen in the FS model is more compact than that in the LS model since the existence of multi-nuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the matter motion determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semi-transparent region.

A dense neutrino gas, such as the one anticipated in the supernova environment, can experience fast neutrino flavor conversions on scales much shorter than those expected in vacuum probably provided that the angular distributions of $\nu_e$ and $\bar\nu_e$ cross each other. We perform a detailed investigation of the neutrino angular distributions obtained by solving the Boltzmann equations for fixed matter profiles of some representative snapshots during the post-bounce phase of core-collapse supernovae in multidimensional calculations of an $11.2\mathrm{M}_{\odot}$ and a $27\mathrm{M}_{\odot}$ progenitor models. Although the $11.2\mathrm{M}_{\odot}$ model features $\nu_e - \bar\nu_e$ angular crossings and the associated fast modes at different time snapshots, the $27\mathrm{M}_{\odot}$ model does not show any crossings within the decoupling region. We show that this can be understood by studying the multipole components of the neutrino distributions. In fact, there is a higher chance for the occurrence of $\nu_e - \bar\nu_e$ angular crossings for the zones where the multipole components of the neutrino distributions are strong enough. We also show that there can exist more than one crossings between the angular distributions of $\nu_e$ and $\bar{\nu}_e$. In addition, apart from the crossings within the neutrino decoupling region, there is a class of $\nu_e-\bar\nu_e$ angular crossings which appears very deep inside the proto-neutron star.

We construct a new equation of state (EOS) for numerical simulations of core-collapse supernovae and neutron-star mergers based on an extended relativistic mean-field model with a small symmetry energy slope $L$, which is compatible with both experimental nuclear data and recent observations of neutron stars. The new EOS table (EOS4) based on the extended TM1 (TM1e) model with $L=40$ MeV is designed in the same tabular form and compared with the commonly used Shen EOS (EOS2) based on the original TM1 model with $L=110.8$ MeV. This is convenient and useful for performing numerical simulations and examining the influences of symmetry energy and its density dependence on astrophysical phenomena. In comparison with the TM1 model used in EOS2, the TM1e model provides a similar maximum neutron-star mass but smaller radius and tidal deformability for a $1.4 M_\odot$ neutron star, which is more consistent with current constraints. By comparing the phase diagram and thermodynamic quantities between EOS4 and EOS2, it is found that the TM1e model predicts relatively larger region of nonuniform matter and softer EOS for neutron-rich matter. Significant differences between EOS4 and EOS2 are observed in the case with low proton fraction, while the properties of symmetric matter remain unchanged.

Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough survey of the ELN crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $\nu_e$ and $\bar\nu_e$ angular distributions, or the electron lepton number (ELN) crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $r$ $\approx$ 16 - 21 km. We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $\nu_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy and nucleosynthesis of heavy nuclei.

The dense neutrino medium in a core-collapse supernova or a neutron-star merger event can experience fast flavor conversions on time/distance scales that are much smaller than those of vacuum oscillations. It is believed that fast neutrino flavor transformation occurs in the region where the angular distributions of $\nu_e$ and $\bar\nu_e$ cross each other. We present the first study of this crossing phenomenon and the fast neutrino flavor conversions in multidimensional (multi-D) supernova models. We examine the neutrino distributions obtained by solving the Boltzmann transport equation for several fixed profiles which are representative snapshots taken from separate 2D and 3D supernova simulations with an $11.2 M_\odot$ progenitor model. Our research shows that the spherically asymmetric patterns of the $\nu_e$ and $\bar\nu_e$ fluxes in multi-D models can assist the appearance of the crossing between the $\nu_e$ and $\bar\nu_e$ angular distributions. In the models that we have studied, there exist unstable neutrino oscillation modes in and beyond the neutrino decoupling region which have amplitude growth rates as large as an $e$-fold per nanosecond. This finding can have important consequences for the explosion mechanism, nucleosynthesis, and neutrino signals of core-collapse supernovae.

We study the influence of density-dependent symmetry energy at high densities in simulations of core-collapse supernovae, black hole formation and proto-neutron star cooling by extending the relativistic mean field (RMF) theory used for the Shen EOS table. We adopt the extended RMF theory to examine the density dependence of the symmetry energy with a small value of the slope parameter $L$, while the original properties of the symmetric nuclear matter are unchanged. In order to assess matter effects at high densities, we perform numerical simulations of gravitational collapse of massive stars adopting the EOS table at high densities beyond $10^{14}$ g/cm$^3$ with the small $L$ value, which is in accord with the experimental and observational constraints, and compare them with the results obtained by using the Shen EOS. Numerical results for 11.2M$_{\odot}$ and 15M$_{\odot}$ stars exhibit minor effects around the core bounce and in the following evolution for 200 ms. Numerical results for 40M$_{\odot}$ and 50M$_{\odot}$ stars reveal a shorter duration toward the black hole formation with a smaller maximum mass for the small $L$ case. Numerical simulations of proto-neutron star cooling over 10 s through neutrino emissions demonstrate increasing effects of the symmetry energy at high densities. Neutrino cooling drastically proceeds in a relatively long timescale with high luminosities and average energies with the small symmetry energy. Evolution toward the cold neutron star is affected because of the different behavior of neutron-rich matter while supernova dynamics around core bounce remains similar in less neutron-rich environments.

We examine the frequencies of the fundamental ($f$-) and the gravity ($g_i$-) mode in gravitational waves from accreting protoneutron stars (PNSs) toward black hole formation. For this purpose, we analyze numerical results of gravitational collapse of massive stars for two different progenitors with three equations of state. We adopt profiles of central objects obtained from the numerical simulations by solving general relativistic neutrino-radiation hydrodynamics under the spherical symmetry. Using a series of snapshots as a static configuration at each time step, we solve the eigenvalue problem to determine the specific frequencies of gravitational waves from the evolving PNSs with accretion from core bounce to black hole formation by the relativistic Cowling approximation. We find that the frequency of the $f$-mode gravitational waves can be expressed as a function of the average density of the PNS at each time step almost independently of the progenitor models in accord with the previous studies on ordinary PNSs. We also find that the ratio of the $g_1$-mode frequency to the $f$-mode frequency is insensitive to the progenitor, but strongly depends on the compactness of the PNSs. Having a simultaneous observation of the $f$- and $g_1$-mode gravitational waves, it would provide the information of the average density and compactness of PNS, which can reconstruct the evolution of mass and radius of the PNSs as a constraint on the equation of state at high densities. Furthermore, we show that the highest average density of the PNS with the maximum mass can be estimated with the detection of the $f$-mode gravitational waves together with the help of the neutrino observations.

Supernova neutrinos are crucially important to probe the final phases of massive star evolution. As is well known from observations of SN1987A, neutrinos provide information on the physical conditions responsible for neutron star formation and on the supernova explosion mechanism. However, there is still no complete understanding of the long-term evolution of neutrino emission in supernova explosions, although there are a number of modern simulations of neutrino radiation hydrodynamics, which study neutrino emission at times less than one second after the bounce. In the present work we systematically calculate the number of neutrinos that can be observed in Super-Kamiokande over periods longer than ten seconds using the database of Nakazato et al. (2013) anticipating that neutrinos from a Galactic supernova can be detected for several tens of seconds. We find that for a supernova at a distance of 10 kpc, neutrinos remain observable for longer than 30 s for a low-mass neutron star ($1.20M_\odot$ gravitational mass) and even longer than 100 s for a high-mass neutron star ($2.05M_\odot$). These scenarios are much longer than the observations of SN1987A and longer than the duration of existing numerical simulations. We propose a new analysis method based on the cumulative neutrino event distribution as a function of reverse time from the last observed event, as a useful probe of the neutron star mass. Our result demonstrates the importance of complete modeling of neutrino light curves in order to extract physical quantities essential for understanding supernova explosion mechanisms, such as the mass and radius of the resulting neutron star.

In this {\it Letter}, we present the result of an axisymmetric CCSN simulation conducted with appropriate treatments of neutrino transport and proper motions of proto-neutron star (PNS), in which a remarkable PNS acceleration is observed in association with asymmetric neutrino emissions by the time of $300$ms after bounce. We find that these asymmetric neutrino emissions play important roles in the acceleration of PNS in this phase. The correlation between the PNS proper motion and the asymmetric ejecta is similar to that in the neutron star (NS) kick of hydrodynamic origin. Both electron-type neutrinos ($\nu_{\rm e}$) and their anti-particles ($\bar{\nu}_{\rm e}$) have $\sim 10\%$ level of asymmetry between the northern and southern hemispheres, while other heavy-leptonic neutrinos ($\nu_x$) have much smaller asymmetry of $\sim 1 \%$. The emissions of $\bar{\nu}_{\rm e}$ and $\nu_x$ are higher in the hemisphere of stronger shock expansion whereas the $\nu_{\rm e}$ emission is enhanced in the opposite hemisphere: in total the neutrinos carry some linear momentum to the hemisphere of the stronger shock expansion. The asymmetry is attributed to the non-spherical distribution of electron-fraction ($Y_e$) in the envelope of PNS. Although it is similar to LESA (lepton-emission self-sustained asymmetry), the $Y_e$ asymmetry seems to be associated with the PNS motion: the latter triggers lateral circular motions in the envelope of PNS by breaking symmetry of the matter distribution there, which are then sustained by a combination of convection, lateral neutrino diffusion and matter-pressure gradient. Our findings may have an influence on the current theories on the NS kick mechanism although long-term simulations are required to assess their impact on the later evolution.

We present a new method for neutrino-matter coupling in multi-dimensional radiation-hydrodynamic simulations of core-collapse supernova (CCSN) with the full Boltzmann neutrino transport. The development is motivated by the fact that the accurate conservation of momentum is required for reliable numerical modelings of CCSN dynamics including a recoil of proto-neutron star (PNS). The new method is built on a hybrid approach, in which we use the energy-momentum tensor of neutrinos to compute the momentum feedback from neutrino to matter in the optically thick region while we employ the collision term in the optically thin region. In this method we utilize a general relativistic description of radiation-hydrodynamics with angular moments, which allows us to evaluate the momentum feedback from neutrino to matter without inconsistency with our Boltzmann solver. We demonstrate that the new method substantially improves the accuracy of linear momentum conservation in our CCSN simulations under reasonable angular resolutions in momentum space, alleviating the difficulty in giving the diffusion limit precisely with the discrete ordinate (Sn) method. It is the first-ever demonstration that the PNS kick can be handled directly and properly in multi-dimensional radiation-hydrodynamic simulations with the full Boltzmann neutrino transport.

Neutrinos are densely populated deep inside the core of massive stars after their gravitational collapse to produce supernova explosions and form compact stars such as neutron stars (NS) and black holes (BH). It has been considered that they may change their flavor identities through so-called fast-pairwise conversions induced by mutual forward scatterings. If that is really the case, the dynamics of supernova explosion will be influenced, since the conversion may occur near the neutrino sphere, from which neutrinos are effectively emitted. In this paper, we conduct a pilot study of such possibilities based on the results of fully self-consistent, realistic simulations of a core-collapse supernova explosion in two spatial dimensions under axisymmetry. As we solved the Boltzmann equations for neutrino transfer in the simulation not as a post-process but in real time, the angular distributions of neutrinos in momentum space for all points in the core at all times are available, a distinct feature of our simulations. We employ some of these distributions extracted at a few selected points and times from the numerical data and apply linear analysis to assess the possibility of the conversion. We focus on the vicinity of the neutrino sphere, where different species of neutrinos move in different directions and have different angular distributions as a result. This is a pilot study for a more thorough survey that will follow soon. We find no positive sign of conversion unfortunately at least for the spatial points and times we studied in this particular model. We hence investigate rather in detail the condition for the conversion by modifying the neutrino distributions rather arbitrarily by hand.

With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of $1\,{\rm rad\,s^{-1}}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches $\sim 20\%$ in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.

We perform an extensive study of the influence of nuclear weak interactions on core-collapse supernovae (CCSNe), paying particular attention to consistency between nuclear abundances in the equation of state (EOS) and nuclear weak interactions. We compute properties of uniform matter based on the variational method. For inhomogeneous nuclear matter, we take a full ensemble of nuclei into account with various finite-density and thermal effects and directly use the nuclear abundances to compute nuclear weak interaction rates. To quantify the impact of a consistent treatment of nuclear abundances on CCSN dynamics, we carry out spherically symmetric CCSN simulations with full Boltzmann neutrino transport, systematically changing the treatment of weak interactions, EOSs, and progenitor models. We find that the inconsistent treatment of nuclear abundances between the EOS and weak interaction rates weakens the EOS dependence of both the dynamics and neutrino signals. We also test the validity of two artificial prescriptions for weak interactions of light nuclei and find that both prescriptions affect the dynamics. Furthermore, there are differences in neutrino luminosities by ~10% and in average neutrino energies by 0.25-1 MeV from those of the fiducial model. We also find that the neutronization burst neutrino signal depends on the progenitor more strongly than on the EOS, preventing a detection of this signal from constraining the EOS.

A critical mass ONe core with a high ignition density is considered to end in gravitational collapse leading to neutron star formation. Being distinct from a Fe core collapse, the final evolution involves combustion flame propagation, in which complex phase transition from ONe elements into the nuclear-statistical-equilibrium (NSE) state takes place. We simulate the core evolution from the O+Ne ignition until the bounce shock penetrates the whole core, using a state-of-the-art 1D Lagrangian neutrino-radiation-hydrodynamic code, in which important nuclear burning, electron capture, and neutrino reactions are taken into account. Special care is also taken in making a stable initial condition by importing the stellar EOS, which is used for the progenitor evolution calculation, and by improving the remapping process. We find that the central ignition leads to intense $\nu_e$ radiation with $L_{\nu_e} \gtrsim 10^{51}$ erg s$^{-1}$ powered by fast electron captures onto NSE isotopes. This pre-bounce $\nu_e$ radiation heats the surroundings by the neutrino-electron scattering, which acts as a new driving mechanism of the flame propagation together with the adiabatic contraction. The resulting flame velocity of $\sim10^8$ cm s$^{-1}$ will be more than one-order-of-magnitude faster than that of laminar flame driven by heat conduction. We also find that the duration of the pre-bounce $\nu_e$ radiation phase depends on the degree of the core hydrostatic/dynamical stability. Therefore, the future detection of the pre-bounce neutrino is important not only to discriminate the ONe core collapse from the Fe core collapse but also to constrain the progenitor hydrodynamical stability.

We investigate the collective neutrino oscillations under the three flavor multi-angle approximation in a spherically symmetric simulation of failed supernovae. A failed supernova emits high neutrino fluxes in a short time, while intense accretion proceeds with high electron density enough to experience re-collapse into a black hole. Our results show that matter-induced effects completely dominate over neutrino self-interaction effects and multi-angle matter suppression occurs at all time snapshots we studied. These facts suggest us that only MSW resonances affect the neutrino flavor conversions in failed supernovae and simple spectra will be observed at neutrino detectors. We also estimate the neutrino event rate in current and future neutrino detectors from a source at 10 kpc as a Galactic event. The time evolution of neutrino detection could provide information about the dense and hot matter and constrain the neutrino mass ordering problem.

We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing different nuclear equations-of-state (EOS's): Lattimer and Swesty's EOS with the incompressibility of K = 220MeV (LS EOS) and Furusawa's EOS based on the relativistic mean field theory with the TM1 parameter set (FS EOS). In the LS EOS the shock wave reaches ~700km at 300ms after bounce and is still expanding whereas in the FS EOS it stalled at ~200km and has started to recede by the same time. This seems to be due to more vigorous turbulent motions in the former during the entire post-bounce phase, which leads to higher neutrino-heating efficiency in the neutrino-driven convection. We also look into the neutrino distributions in momentum space, which is the advantage of the Boltzmann transport over other approximate methods. We find non-axisymmetric angular distributions with respect to the local radial direction, which also generate off-diagonal components of the Eddington tensor. We find that the r {\theta}-component reaches ~10% of the dominant rr-component and, more importantly, it dictates the evolution of lateral neutrino fluxes, dominating over the {\theta}{\theta}-component, in the semi-transparent region. These data will be useful to further test and possibly improve the prescriptions used in the approximate methods.

We revisit the diffuse supernova neutrino background in light of recent systematic studies of stellar core collapse that reveal the quantitative impacts of the progenitor conditions on the collapse process. In general, the dependence of the progenitor on the core-collapse neutrino emission is not monotonic in progenitor initial mass, but we show that it can, at first order, be characterized by the core compactness. For the first time, we incorporate the detailed variations in the neutrino emission over the entire mass range $8$-$100 {\rm M}_\odot$, based on (i) a long-term simulation of the core collapse of a $8.8 {\rm M}_\odot$ O-Ne-Mg core progenitor, (ii) over 100 simulations of iron core collapse to neutron stars, and (iii) half a dozen simulations of core collapse to black holes (the "failed channel"). The fraction of massive stars that undergo the failed channel remains uncertain, but in view of recent simulations which reveal high compactness to be conducive to collapse to black holes, we characterize the failed fraction by considering a threshold compactness above which massive stars collapse to black holes and below which the final remnant is a neutron star. We predict that future detections of the diffuse supernova neutrino background may have the power to reveal this threshold compactness, if its value is relatively small as suggested by interpretations of several recent astronomical observations.

The mechanism driving core-collapse supernovae is sensitive to the interplay between matter and neutrino radiation. However, neutrino radiation transport is very difficult to simulate, and several radiation transport methods of varying levels of approximation are available. We carefully compare for the first time in multiple spatial dimensions the discrete ordinates (DO) code of Nagakura, Yamada, and Sumiyoshi and the Monte Carlo (MC) code Sedonu, under the assumptions of a static fluid background, flat spacetime, elastic scattering, and full special relativity. We find remarkably good agreement in all spectral, angular, and fluid interaction quantities, lending confidence to both methods. The DO method excels in determining the heating and cooling rates in the optically thick region. The MC method predicts sharper angular features due to the effectively infinite angular resolution, but struggles to drive down noise in quantities where subtractive cancellation is prevalent, such as the net gain in the protoneutron star and off-diagonal components of the Eddington tensor. We also find that errors in the angular moments of the distribution functions induced by neglecting velocity dependence are sub-dominant to those from limited momentum-space resolution. We briefly compare directly computed second angular moments to those predicted by popular algebraic two-moment closures, and find that the errors from the approximate closures are comparable to the difference between the DO and MC methods. Included in this work is an improved Sedonu code, which now implements a fully special relativistic, time-independent version of the grid-agnostic Monte Carlo random walk approximation.

We have constructed a nuclear equation of state (EOS) that includes a full nuclear ensemble for use in core-collapse supernova simulations. It is based on the EOS for uniform nuclear matter that two of the authors derived recently, applying a variational method to realistic two- and there-body nuclear forces. We utilize an extended liquid drop model of heavy nuclei and a quantum-theoretical mass evaluation for light nuclei. In addition to realistic nuclear forces, the inclusion of in-medium effects on the full ensemble of nuclei makes the new EOS one of the most realistic EOS's for supernova simulations. We make comparisons with the FYSS EOS, which is based on the same formulation for the nuclear ensemble but adopts the relativistic mean field (RMF) theory with the TM1 parameter set for uniform nuclear matter. The new EOS is softer than the FYSS EOS around and above nuclear saturation densities. We find that neutron-rich nuclei with small mass numbers are more abundant in the new EOS than in the FYSS EOS because of the larger saturation densities and smaller symmetry energy of nuclei in the former. We apply the two EOS's to 1D supernova simulations and find that the new EOS gives lower electron fractions and higher temperatures in the collapse phase owing to the smaller symmetry energy. As a result, the inner core has smaller masses for the new EOS. It is more compact, on the other hand, due to the softness of the new EOS and bounces at higher densities. The ensuing outward propagations of the shock wave in the outer core are very similar in the two simulations, which may be an artifact, though, caused by the use of the same tabulated electron capture rates for heavy nuclei ignoring differences in the nuclear composition between the two EOS's in these computations.

We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse supernova simulations. The abundance of various nuclei is obtained together with thermodynamic quantities. The formulation is an extension of the previous model, in which we adopted the relativistic mean field theory with the TM1 parameter set for nucleons, the quantum approach for $d$, $t$, $h$ and $\alpha$ as well as the liquid drop model for the other nuclei under the nuclear statistical equilibrium. We reformulate the model of the light nuclei other than $d$, $t$, $h$ and $\alpha$ based on the quasi-particle description. Furthermore, we modify the model so that the temperature dependences of surface and shell energies of heavy nuclei could be taken into account. The pasta phases for heavy nuclei and the Pauli- and self-energy shifts for $d$, $t$, $h$ and $\alpha$ are taken into account in the same way as in the previous model. We find that nuclear composition is considerably affected by the modifications in this work, whereas thermodynamical quantities are not changed much. In particular, the washout of shell effect has a great impact on the mass distribution above $T \sim 1$ MeV. This improvement may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores.

We investigate the influences of the nuclear composition on the weak interaction rates of heavy nuclei during the core collapse of massive stars. The nuclear abundances in nuclear statistical equilibrium (NSE) are calculated by some equation of state (EOS) models including in-medium effects on nuclear masses. We systematically examine the sensitivities of electron capture and neutrino-nucleus scattering on heavy nuclei to the nuclear shell effects and the single nucleus approximation. We find that the washout of shell effects at high temperatures brings significant change to weak rates by smoothing the nuclear abundance distribution: the electron capture rate decreases by $\sim$20$\%$ in the early phase and increases by $\sim$40$\%$ in the late phase at most, while the cross section for neutrino-nucleus scattering is reduced by $\sim$15$\%$. This is because the open-shell nuclei become abundant instead of those with closed neutron shells as the shell effects disappear. We also find that the single-nucleus description based on the average values leads to underestimations of weak rates. Electron captures and neutrino coherent scattering on heavy nuclei are reduced by $\sim$80$\%$ in the early phase and by $\sim$5$\%$ in the late phase, respectively. These results indicate that NSE like EOS accounting for shell washout is indispensable for the reliable estimation of weak interaction rates in simulations of core-collapse supernovae.

We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3+1 formalism of general relativity (GR), making use of the shift vector to specify the time evolution of the coordinate system. This means that the transport part of our code is essentially general relativistic although in this paper it is applied only to the moving curvilinear coordinates in the flat Minknowski spacetime, since the gravity part is still Newtonian. The numerical aspect of the implementation is also described in detail. Employing the axisymmetric two-dimensional version of the code, we conduct two test computations: oscillations and runaways of proto-neutron star (PNS). We show that our new method works fine, tracking the motions of PNS correctly. We believe that this is a major advancement toward the realistic simulation of CCSNe.

Energetics of nuclear reaction is fundamentally important to understand the mechanism of pair instability supernovae (PISNe). Based on the hydrodynamic equations and thermodynamic relations, we derive exact expressions for energy conservation suitable to be solved in simulation. We also show that some formulae commonly used in the literature are obtained as approximations of the exact expressions. We simulate the evolution of very massive stars of ~100-320 Msun with zero- and 1/10 Zsun, and calculate further explosions as PISNe, applying each of the exact and approximate formulae. The calculations demonstrate that the explosion properties of PISN, such as the mass range, the 56Ni yield, and the explosion energy, are significantly affected by applying the different energy generation rates. We discuss how these results affect the estimate of the PISN detection rate, which depends on the theoretical predictions of such explosion properties.

We study the multi-dimensional properties of neutrino transfer inside supernova cores by solving the Boltzmann equations for neutrino distribution functions in genuinely six dimensional (6D) phase space. Adopting representative snapshots of the post-bounce core from other supernova simulations in three dimensions, we solve the temporal evolutions to stationary states of neutrino distribution functions by our Boltzmann solver. Taking advantage of the multi-angle and multi-energy feature realized by the S$_n$ method in our code, we reveal the genuine characteristics of spatially three dimensional (3D) neutrino transfer such as non-radial fluxes and non-diagonal Eddington tensors. In addition, we assess the ray-by-ray approximation, turning off the lateral-transport terms in our code. We demonstrate that the ray-by-ray approximation tends to propagate fluctuations in thermodynamical states around the neutrino-sphere along each radial ray and overestimate the variations between the neutrino distributions on different radial rays. We find that the difference in the densities and fluxes of neutrinos between the ray-by-ray approximation and the full Boltzmann transport becomes ~20%, which is also the case for the local heating rate, whereas the volume-integrated heating rate in the Boltzmann transport is found to be only slightly larger (~2%) than the counterpart in the ray-by-ray approximation due to cancellation among different rays. These results suggest that we had better assess carefully the possible influences of various approximations in the neutrino transfer employed in the current simulations on supernova dynamics. Detailed information on the angle and energy moments of neutrino distribution functions will be profitable for the future development of numerical methods in neutrino-radiation hydrodynamics.

Nuclear equation of state (EoS) plays an important role in understanding the formation of compact objects such as neutron stars and black holes. The true nature of the EoS has been a matter of debate at any density range not only in the nuclear physics but also in the astronomy and astrophysics. We have constructed a database of EoSs by compiling data from the literature. Our database contains the basic properties of the nuclear EoS of symmetric nuclear matter and of pure neutron matter. It also includes the detailed information about the theoretical models, for example the adopted methods and assumptions in individual models. The novelty of the database is to consider new experimental probes such as the symmetry energy, its slope relative to the baryon density, and the incompressibility, which enables the users to check their model dependences. We demonstrate the performance of the EOSDB through the examinations of the model dependence among different nuclear EoSs. It is reveled that some theoretical EoSs, which is commonly used in astrophysics, do not satisfactorily agree with the experimental constraints.

We propose a novel numerical method for solving multi-dimensional, special relativistic Boltzmann equations for neutrinos coupled to hydrodynamics equations. It is meant to be applied to simulations of core-collapse supernovae. We handle special relativity in a non-conventional way, taking account of all orders of v/c. Consistent treatment of advection and collision terms in the Boltzmann equations is the source of difficulties, which we overcome by employing two different energy grids: Lagrangian remapped and laboratory fixed grids. We conduct a series of basic tests and perform a one-dimensional simulation of core-collapse, bounce and shock-stall for a 15M_{sun} progenitor model with a minimum but essential set of microphysics. We demonstrate in the latter simulation that our new code is capable of handling all phases in core-collapse supernova. For comparison, a non-relativistic simulation is also conducted with the same code, and we show that they produce qualitatively wrong results in neutrino transfer. Finally, we discuss a possible incorporation of general relativistic effects in our method.

The equation of state (EOS) for hot asymmetric nuclear matter which is constructed with the variational method starting from the Argonne v18 and Urbana IX nuclear forces is applied to spherically symmetric core-collapse supernovae (SNe). We first investigate the EOS of isentropic beta-stable SN matter, and find that the matter with the variational EOS is more neutron-rich than that with the Shen EOS. Using the variational EOS for uniform matter supplemented by the Shen EOS of non-uniform matter at low densities, we perform general-relativistic spherically symmetric simulations of core-collapse SNe with and without neutrino transfer, starting from a presupernova model of 15 solar mass. In the adiabatic simulation without neutrino transfer, the explosion is successful, and the explosion energy with the variational EOS is larger than that with the Shen EOS. In the case of the simulation with neutrino transfer, the shock wave stalls and then the explosion fails, as in other spherically symmetric simulations. The inner core with the variational EOS is more compact than that with the Shen EOS, due to the relative softness of the variational EOS. This implies that the variational EOS is more advantageous for SN explosions than the Shen EOS.

Neutrino emissions from electron/positron capture on the deuteron and the nucleon-nucleon fusion processes in the surface region of a supernova core are studied. These weak processes are evaluated in the standard nuclear physics approach, which consists of one-nucleon and two-nucleon-exchange currents and nuclear wave functions generated by a high precision nucleon-nucleon potential. In addition to the cross sections for these processes involving the deuteron, we present neutrino emissivities due to these processes calculated for typical profiles of core-collapsed supernovae. These novel neutrino emissivities are compared with the standard neutrino emission mechanisms. We find that the neutrino emissivity due to the electron capture on the deuteron is comparable to that on the proton in the deuteron abundant region. The implications of the new channels involving deuterons for the supernova mechanism are discussed.

Hadronic matter undergoes a deconfinement transition to quark matter at high temperature and/or high density. It would be realized in collapsing cores of massive stars. In the framework of MIT bag model, the ambiguities of the interaction are encapsulated in the bag constant. Some progenitor stars that invoke the core collapses explode as supernovae, and other ones become black holes. The fates of core collapses are investigated for various cases. Equations of state including the hadron-quark phase transition are constructed for the cases of the bag constant B=90, 150 and 250 MeV fm^{-3}. To describe the mixed phase, the Gibbs condition is used. Adopting the equations of state with different bag constants, the core collapse simulations are performed for the progenitor models with 15 and 40Msolar. If the bag constant is small as B=90 MeV fm^{-3}, an interval between the bounce and black hole formation is shortened drastically for the model with 40Msolar and the second bounce revives the shock wave leading to explosion for the model with 15Msolar.

We perform numerical experiments to investigate the influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability (SASI). The time evolution of shock waves is calculated with a simple light-bulb approximation for the neutrino transport and a multi-nuclei equation of state. The neutrino absorptions and inelastic interactions with deuterons, tritons, helions and alpha particles are taken into account in the hydrodynamical simulations in addition to the ordinary charged-current interactions with nucleons. Axial symmetry is assumed but no equatorial symmetry is imposed. We show that the heating rates of deuterons reach as high as 10% of those of nucleons around the bottom of the gain region. On the other hand, alpha particles are heated near the shock wave, which is important when the shock wave expands and the density and temperature of matter become low. It is also found that the models with heating by light nuclei have different evolutions from those without it in the non-linear phase of SASI. This results is because matter in the gain region has a varying density and temperature and there appear sub-regions that are locally rich in deuterons and alpha particles. Although the light nuclei are never dominant heating sources and they work favorably for shock revival in some cases and unfavorably in other cases, they are non-negligible and warrant further investigation.

We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with thermodynamic quantities. A model free energy is constructed, based on the relativistic mean field theory for nucleons and the mass formula for nuclei with the proton number up to ~ 1000. The formulation is an extension of the previous model, in which we adopted the liquid drop model to all nuclei under the nuclear statistical equilibrium. We reformulate the new liquid drop model so that the temperature dependences of bulk energies could be taken into account. Furthermore, we extend the region in the nuclear chart, in which shell affects are included, by using theoretical mass data in addition to experimental ones. We also adopt a quantum theoretical mass evaluation of light nuclei, which incorporates the Pauli- and self-energy shifts that are not included in the ordinary liquid drop model. The pasta phases for heavy nuclei are taken into account in the same way as in the previous model. We find that the abundances of heavy nuclei are modified by the shell effects of nuclei and temperature dependence of bulk energies. These changes may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. The abundances of light nuclei are also modified by the new mass evaluation, which may affect the heating and cooling rates of supernova cores and shocked envelopes.

We compare three different statistical models for the equation of state (EOS) of stellar matter at subnuclear densities and temperatures (0.5-10 MeV) expected to occur during the collapse of massive stars and supernova explosions. The models introduce the distributions of various nuclear species in nuclear statistical equilibrium, but use somewhat different nuclear physics inputs. It is demonstrated that the basic thermodynamical quantities of stellar matter under these conditions are similar, except in the region of high densities and low temperatures. We demonstrate that mass and isotopic distributions have considerable differences related to the different assumptions of the models on properties of nuclei at these stellar conditions. Overall, the three models give similar trends, but the details reflect the uncertainties related to the modeling of medium effects, such as the temperature and density dependence of surface and bulk energies of heavy nuclei, and the nuclear shell structure effects. We discuss importance of new physics inputs for astrophysical calculations from experimental data obtained in intermediate energy heavy-ion collisions, in particular, the similarities of the conditions reached during supernova explosions and multifragmentation reactions.

We present a new series of supernova neutrino light curves and spectra calculated by numerical simulations for a variety of progenitor stellar masses (13-50Msolar) and metallicities (Z = 0.02 and 0.004), which would be useful for a broad range of supernova neutrino studies, e.g., simulations of future neutrino burst detection by underground detectors, or theoretical predictions for the relic supernova neutrino background. To follow the evolution from the onset of collapse to 20 s after the core bounce, we combine the results of neutrino-radiation hydrodynamic simulations for the early phase and quasi-static evolutionary calculations of neutrino diffusion for the late phase, with different values of shock revival time as a parameter that should depend on the still unknown explosion mechanism. We here describe the calculation methods and basic results including the dependence on progenitor models and the shock revival time. The neutrino data are publicly available electronically.

Massive stars (M> 10Msun) end their lives with spectacular explosions due to gravitational collapse. The collapse turns the stars into compact objects such as neutron stars and black holes with the ejection of cosmic rays and heavy elements. Despite the importance of these astrophysical events, the mechanism of supernova explosions has been an unsolved issue in astrophysics. This is because clarification of the supernova dynamics requires the full knowledge of nuclear and neutrino physics at extreme conditions, and large-scale numerical simulations of neutrino radiation hydrodynamics in multi-dimensions. This article is a brief overview of the understanding (with difficulty) of the supernova mechanism through the recent advance of numerical modeling at supercomputing facilities. Numerical studies with the progress of nuclear physics are applied to follow the evolution of compact objects with neutrino emissions in order to reveal the birth of pulsars/black holes from the massive stars.

This is a status report on our endeavor to reveal the mechanism of core-collapse supernovae (CCSNe) by large-scale numerical simulations. Multi-dimensionality of the supernova engine, general relativistic magnetohydrodynamics, energy and lepton number transport by neutrinos emitted from the forming neutron star as well as nuclear interactions there, are all believed to play crucial roles in repelling infalling matter and producing energetic explosions. These ingredients are nonlinearly coupled with one another in the dynamics of core-collapse, bounce, and shock expansion. Serious quantitative studies of CCSNe hence make extensive numerical computations mandatory. Since neutrinos are neither in thermal nor in chemical equilibrium in general, their distributions in the phase space should be computed. This is a six dimensional (6D) neutrino transport problem and quite a challenge even for those with an access to the most advanced numerical resources such as the "K computer". To tackle this problem, we have embarked on multi-front efforts. In particular we report in this paper our recent progresses in the treatments of multi-dimensional (multi-D) radiation-hydrodynamics. We are currently proceeding on two different paths to the ultimate goal; in one approach we employ an approximate but highly efficient scheme for neutrino transport and treat 3D hydrodynamics and/or general relativity rigorously; some neutrino-driven explosions will be presented and comparisons will be made between 2D and 3D models quantitatively; in the second approach, on the other hand, exact but so far Newtonian Boltzmann equations are solved in two and three spatial dimensions; we will show some demonstrative test simulations. We will also address the perspectives of exa-scale computations on the next generation supercomputers.

We develop a numerical code to calculate the neutrino transfer with multi-energy and multi-angle in three dimensions (3D) for the study of core-collapse supernovae. The numerical code solves the Boltzmann equations for neutrino distributions by the discrete-ordinate (S_n) method with a fully implicit differencing for time advance. The Boltzmann equations are formulated in the inertial frame with collision terms being evaluated to the zeroth order of v/c. A basic set of neutrino reactions for three neutrino species is implemented together with a realistic equation of state of dense matter. The pair process is included approximately in order to keep the system linear. We present numerical results for a set of test problems to demonstrate the ability of the code. The numerical treatments of advection and collision terms are validated first in the diffusion and free streaming limits. Then we compute steady neutrino distributions for a background extracted from a spherically symmetric, general relativistic simulation of 15Msun star and compare them with the results in the latter computation. We also demonstrate multi-D capabilities of the 3D code solving neutrino transfers for artificially deformed supernova cores in 2D and 3D. Formal solutions along neutrino paths are utilized as exact solutions. We plan to apply this code to the 3D neutrino-radiation hydrodynamics simulations of supernovae. This is the first article in a series of reports on the development.

We construct the equation of state (EOS) of dense matter covering a wide range of temperature, proton fraction, and density for the use of core-collapse supernova simulations. The study is based on the relativistic mean-field (RMF) theory, which can provide an excellent description of nuclear matter and finite nuclei. The Thomas--Fermi approximation in combination with assumed nucleon distribution functions and a free energy minimization is adopted to describe the non-uniform matter, which is composed of a lattice of heavy nuclei. We treat the uniform matter and non-uniform matter consistently using the same RMF theory. We present two sets of EOS tables, namely EOS2 and EOS3. EOS2 is an update of our earlier work published in 1998 (EOS1), where only the nucleon degree of freedom is taken into account. EOS3 includes additional contributions from $\Lambda$ hyperons. The effect of $\Lambda$ hyperons on the EOS is negligible in the low-temperature and low-density region, whereas it tends to soften the EOS at high density. In comparison with EOS1, EOS2 and EOS3 have an improved design of ranges and grids, which covers the temperature range $T=0.1$--$10^{2.6}$ MeV with the logarithmic grid spacing $\Delta \log_{10}(T/\rm{[MeV]})=0.04$ (92 points including T=0), the proton fraction range $Y_p=0$--0.65 with the linear grid spacing $\Delta Y_p = 0.01$ (66 points), and the density range $\rho_B=10^{5.1}$--$10^{16}\,\rm{g\,cm^{-3}}$ with the logarithmic grid spacing $\Delta \log_{10}(\rho_B/\rm{[g\,cm^{-3}]}) = 0.1$ (110 points).

We investigate the emergence of hyperons in black-hole-forming failed supernovae, which are caused by the dynamical collapse of nonrotating massive stars. We perform neutrino-radiation hydrodynamical simulations in general relativity adopting realistic hyperonic equation-of-state (EOS). Attractive and repulsive cases are examined for the potential of \Sigma\ hyperons. Since hyperons soften the EOS, they shorten the time interval from the bounce to black hole formation, which corresponds to the duration of neutrino emission. This effect is larger for the attractive case than the repulsive case because \Sigma\ hyperons appear more easily. In addition, we investigate the impacts of pions to find that they also promotes the recollapse towards the black hole formation.

We discuss the possibility to probe the QCD critical point during the dynamical black hole formation from a gravitational collapse of a massive star, where the temperature and the baryon chemical potential become as high as T ~ 90 MeV and $\mu_B$ ~ 1300 MeV. Comparison with the phase boundary in chiral effective models suggests that quark matter is likely to be formed before the horizon is formed. Furthermore, the QCD critical point may be probed during the black hole formation. The critical point is found to move in the lower temperature direction in asymmetric nuclear matter, and in some of the chiral models it is found to be in the reachable region during the black hole formation processes.

We calculate a new equation of state for baryons at sub-nuclear densities meant for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with the thermodynamic quantities. The formulation is the NSE description and the liquid drop approximation of nuclei. The model free energy to minimize is calculated by relativistic mean field theory for nucleons and the mass formula for nuclei with the atomic number up to ~ 1000. We have also taken into account the pasta phase, thanks to which the transition to uniform nuclear matter in our EOS occurs in the conventional manner: nuclei are not dissociated to nucleons but survive right up to the transition to uniform nuclear matter. We find that the free energy and other thermodynamical quantities are not very different from those given in the Shen's EOS, one of the standard EOS's that adopt the single nucleus approximation. The average mass is systematically different, on the other hand, which may have an important ramification to the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. It is also interesting that the root mean square of the mass number is not very different from the average mass number, since the former is important for the evaluation of coherent scattering rates on nuclei but has been unavailable so far. The EOS table is currently under construction, which will include the weak interaction rates.

In the formation process of black holes, the density and temperature of matter become sufficiently high for quarks and pions to appear. In this study we numerically investigate stellar core collapse and black hole formation taking into account the equations of state involving quarks and/or pions. In our simulations, we utilize a code that solves the general relativistic hydrodynamics and neutrino transfer equations simultaneously, treating neutrino reactions in detail under spherical symmetry. Initial models with three different masses, namely, 40, 100 and 375Msolar, are adopted. Our results show that quarks and pions shorten the duration of neutrino emission if the collapse bounces before black hole formation. In addition, pions increase the luminosity and average energy of neutrinos before black hole formation. We also find that the hadron-quark phase transition leads to an interesting evolution of temperature. Moreover, the neutrino event number is evaluated for the currently operating neutrino detector, SuperKamiokande, to confirm that it is not only detectable but also affected by the emergence of quarks and pions for Galactic events. While there are some issues, such as hyperons, beyond the scope of this study, this is the first serious attempt to assess the impact of quarks and pions in dynamical simulations of black hole formation and will serve as an important foundation for future studies.

The detection of neutrinos from massive stellar collapses can teach us a lot not only about source objects but also about microphysics working deep inside them. In this study we discuss quantitatively the possibility to extract information on the properties of dense and hot hadronic matter from neutrino signals coming out of black-hole-forming collapses of non-rotational massive stars. Based on our detailed numerical simulations we evaluate the event numbers for SuperKamiokande with neutrino oscillations being fully taken into account. We demonstrate that the event numbers from a Galactic event are large enough not only to detect it but also to distinguish one hadronic equation of state from another by our statistical method assuming the same progenitor model and non-rotation. This means that the massive stellar collapse can be a unique probe into hadron physics and will be a promising target of the nascent neutrino astronomy.

We study how the influence of the shock wave appears in neutrino oscillations and the neutrino spectrum using density profile of adiabatic explosion model of a core-collapse supernova which is calculated in an implicit Lagrangian code for general relativistic spherical hydrodynamics. We calculate expected event rates of neutrino detection at SK and SNO for various theta_{13} values and both normal and inverted hierarchies. The predicted event rates of bar{nu}_e and nu_e depend on the mixing angle theta_{13} for the inverted and normal hierarchies, respectively, and the influence of the shock appears for about 2 - 8 s when sin^2 2 theta_{13} is larger than 10^{-3}. These neutrino signals for the shock propagation is decreased by < 30 % for bar{nu}_e in inverted (SK) or by < 15 % for nu_e in normal hierarchy (SNO) compared with the case without shock. The obtained ratio of the total event for high-energy neutrinos (20 MeV < E_{nu} < 60 MeV) to low-energy neutrinos (5 MeV < E_{nu} < 20 MeV) is consistent with the previous studies in schematic semi-analytic or other hydrodynamic models of the shock propagation. The time dependence of the calculated ratio of the event rates of high-energy to low-energy neutrinos is a very useful observable which is sensitive to theta_{13} and hierarchies. Namely, time-dependent ratio shows clearer signal of the shock propagation that exhibits remarkable decrease by at most factor \sim 2 for bar{nu}_e in inverted (SK), whereas it exhibits smaller change by \sim 10 % for nu_e in normal hierarchy (SNO). Observing time-dependent high-energy to low-energy ratio of the neutrino events thus would provide a piece of very useful information to constrain theta_{13} and mass hierarchy, and eventually help understanding the propagation how the shock wave propagates inside the star.

We examine a potential role of the neutrino deuteron reactions in the mechanism of supernova explosion by evaluating the energy transfer cross section for the neutrino heating. We calculate the energy loss rate due to the neutrino absorptions through the charged-current process as well as the neutrino scattering through the neutral-current process. In so doing, we adopt a detailed evaluation of cross sections for the neutrino deuteron reactions with the phenomenological Lagrangian approach. We find the energy transfer cross section for the deuteron is larger than those for $^{3}$H, $^{3}$He and $^{4}$He for neutrino temperatures (T$_\nu$ $\sim 4$ MeV) relevant to supernova core. Because of the low energy threshold for the deuteron breakup, the energy transfer rate rapidly increases from low temperature, T$_\nu$ $\sim 1$ MeV. This suggests that the neutrino deuteron reactions may contribute effectively to the heating mechanism during the dissociation of irons into light elements and nucleons in the shocked material of supernova core.

The equation of state and composition of matter are calculated for conditions typical for pre-collapse and early collapse stages in core collapse supernovae. The composition is evaluated under the assumption of nuclear statistical equilibrium, when the matter is considered as an `almost' ideal gas with corrections due to thermal excitations of nuclei, to free nucleon degeneracy, and to Coulomb and surface energy corrections. The account of these corrections allows us to obtain the composition for densities a bit below the nuclear matter density. Through comparisons with the equation of state (EOS) developed by Shen et al. we examine the approximation of one representative nucleus used in most of recent supernova EOS's. We find that widely distributed compositions in the nuclear chart are different, depending on the mass formula, while the thermodynamical quantities are quite close to those in the Shen's EOS.

We investigate the emergence of strange baryons in the dynamical collapse of a non-rotating massive star to a black hole by the neutrino-radiation hydrodynamical simulations in general relativity. By following the dynamical formation and collapse of nascent proto-neutron star from the gravitational collapse of a 40Msun star adopting a new hyperonic EOS table, we show that the hyperons do not appear at the core bounce but populate quickly at ~0.5-0.7 s after the bounce to trigger the re-collapse to a black hole. They start to show up off center owing to high temperatures and later prevail at center when the central density becomes high enough. The neutrino emission from the accreting proto-neutron star with the hyperonic EOS stops much earlier than the corresponding case with a nucleonic EOS while the average energies and luminosities are quite similar between them. These features of neutrino signal are a potential probe of the emergence of new degrees of freedom inside the black hole forming collapse.

Recently, stellar collapse involving black hole formation from massive stars is suggested to emit enormous fluxes of neutrinos on par with ordinary core-collapse supernovae. We investigate their detectability for the currently operating neutrino detector, SuperKamiokande. Neutrino oscillation is also taken into account for the evaluation. We find that the event number is larger than or comparable to that of supernova neutrinos and, hence, black hole formation is also a candidate for neutrino astronomy. Moreover, we find that the event number depends dominantly on the equation of state used in the computations of the black hole formation. This fact implies that the detection of neutrinos emitted from the black hole progenitors is very valuable to probe the properties of the equation of state for hot and/or dense matter.

We study the progenitor dependence of the black hole formation and its associated neutrino signals from the gravitational collapse of non-rotating massive stars, following the preceding study on the single progenitor model in Sumiyoshi et al. (2007). We aim to clarify whether the dynamical evolution toward the black hole formation occurs in the same manner for different progenitors and to examine whether the characteristic of neutrino bursts is general having the short duration and the rapidly increasing average energies. We perform the numerical simulations by general relativistic neutrino-radiation hydrodynamics to follow the dynamical evolution from the collapse of pre-supernova models of 40Msun and 50Msun toward the black hole formation via contracting proto-neutron stars. For the three progenitor models studied in this paper, we found that the black hole formation occurs in ~0.4-1.5 s after core bounce through the increase of proto-neutron star mass together with the short and energetic neutrino burst. We found that density profile of progenitor is important to determine the accretion rate onto the proto-neutron star and, therefore, the duration of neutrino burst. We compare the neutrino bursts of black hole forming events from different progenitors and discuss whether we can probe clearly the progenitor and/or the dense matter.

We explore the abundance of light clusters in core-collapse supernovae at post-bounce stage in a quantum statistical approach. Adopting the profile of a supernova core from detailed numerical simulations, we study the distribution of light bound clusters up to alpha particles (A=2-4) as well as heavy nuclei (A > 4) in dense matter at finite temperature. Within the frame of a cluster-mean field approach, the abundances of light clusters are evaluated accounting for self-energy, Pauli blocking and effects of continuum correlations. We find that deuterons and tritons, in addition to 3He and 4He, appear abundantly in a wide region from the surface of the proto-neutron star to the position of the shock wave. The appearance of light clusters may modify the neutrino emission in the cooling region and the neutrino absorption in the heating region, and, thereby, influence the supernova mechanism.

We construct an equation of state including the hadron-quark phase transition. The mixed phase is obtained by the Gibbs conditions for finite temperature. We adopt the equation of state based on the relativistic mean field theory for the hadronic phase taking into account pions. As for the quark phase, the MIT bag model of the deconfined 3-flavor strange quark matter is used. As a result, our equation of state is thermodynamically stable and exhibits qualitatively the desired properties of hadron-quark mixed matter, such as the temperature dependence of the transition density. The pions raise the transition density because they make the equation of state softer. Using the equation of state constructed here, we study its astrophysical implications. The maximum mass of compact stars is investigated, and our equation of state is consistent with recent observations. We also compute the collapse of a massive star with 100 solar masses using our equation of state and find that the interval time from the bounce to the black hole formation becomes shorter for the model with pions and quarks. The pions and quarks affect the total energy of the emitted neutrinos because the duration time of the neutrino emission becomes shorter. The neutrino luminosity rises under the effect of pions since the density of the proto-neutron star becomes high.

We investigate the role of Alfven waves in the core-collapse supernova (SN) explosion. We assume that Alfven waves are generated by convections inside a proto-neutron star (PNS) and emitted from its surface. Then these waves propagate outwards, dissipate via nonlinear processes, and heat up matter around a stalled prompt shock. To quantitatively assess the importance of this process for the revival of the stalled shock, we perform 1D time-dependent hydrodynamical simulations, taking into account the heating via the dissipation of Alfven waves that propagate radially outwards along open flux tubes. We show that the shock revival occurs if the surface field strength is larger than ~2e15 G and if the amplitude of velocity fluctuation at the PNS surface is larger than 20% of the local sound speed. Interestingly, the Alfven wave mechanism is self-regulating in the sense that the explosion energy is not very sensitive to the surface field strength and initial amplitude of Alfven waves as long as they are larger than the threshold values given above.

We study the black hole formation and the neutrino signal from the gravitational collapse of a non-rotating massive star of 40 Msun. Adopting two different sets of realistic equation of state (EOS) of dense matter, we perform the numerical simulations of general relativistic neutrino-radiation hydrodynamics under the spherical symmetry. We make comparisons of the core bounce, the shock propagation, the evolution of nascent proto-neutron star and the resulting re-collapse to black hole to reveal the influence of EOS. We also explore the influence of EOS on the neutrino emission during the evolution toward the black hole formation. We find that the speed of contraction of the nascent proto-neutron star, whose mass increases fast due to the intense accretion, is different depending on the EOS and the resulting profiles of density and temperature differ significantly. The black hole formation occurs at 0.6-1.3 sec after bounce when the proto-neutron star exceeds its maximum mass, which is crucially determined by the EOS. We find that the average energies of neutrinos increase after bounce because of rapid temperature increase, but at different speeds depending on the EOS. The duration of neutrino emission up to the black hole formation is found different according to the different timing of re-collapse. These characteristics of neutrino signatures are distinguishable from those for ordinary proto-neutron stars in successful core-collapse supernovae. We discuss that a future detection of neutrinos from black-hole-forming collapse will contribute to reveal the black hole formation and to constrain the EOS at high density and temperature.

An equation of state (EOS) for uniform nuclear matter is constructed at zero and finite temperatures with the variational method starting from the realistic nuclear Hamiltonian composed of the Argonne V18 and UIX potentials. The energy is evaluated in the two-body cluster approximation with the three-body-force contribution treated phenomenologically so as to reproduce the empirical saturation conditions. The obtained energies for symmetric nuclear matter and neutron matter at zero temperature are in fair agreement with those by Akmal, Pandharipande and Ravenhall, and the maximum mass of the neutron star is 2.2 Msolar. At finite temperatures, a variational method by Schmidt and Pandharipande is employed to evaluate the free energy, which is used to derive various thermodynamic quantities of nuclear matter necessary for supernova simulations. The result of this variational method at finite temperatures is found to be self-consistent.

The existence of various anomalous stars, such as the first stars in the universe or stars produced by stellar mergers, has been recently proposed. Some of these stars will result in black hole formation. In this study, we investigate iron core collapse and black hole formation systematically for the iron-core mass range of 3 - 30Msun, which has not been studied well so far. Models used here are mostly isentropic iron cores that may be produced in merged stars in the present universe but we also employ a model that is meant for a Population III star and is obtained by evolutionary calculation. We solve numerically the general relativistic hydrodynamics and neutrino transfer equations simultaneously, treating neutrino reactions in detail under spherical symmetry. As a result, we find that massive iron cores with ~10Msun unexpectedly produce a bounce owing to the thermal pressure of nucleons before black hole formation. The features of neutrino signals emitted from such massive iron cores differ in time evolution and spectrum from those of ordinary supernovae. Firstly, the neutronization burst is less remarkable or disappears completely for more massive models because the density is lower at the bounce. Secondly, the spectra of neutrinos, except the electron type, are softer owing to the electron-positron pair creation before the bounce. We also study the effects of the initial density profile, finding that the larger the initial density gradient is, the more steeply the neutronization burst declines. Further more, we suggest a way to probe into the black hole progenitors from the neutrino emission and estimate the event number for the currently operating neutrino detectors.

The gravitational collapse of a non-rotating, black-hole-forming massive star is studied by neutrino-radiation-hydrodynamical simulations for two different sets of realistic equation of state of dense matter. We show that the event will produce as many neutrinos as the ordinary supernova, but with distinctive characteristics in luminosities and spectra that will be an unmistakable indication of black hole formation. More importantly, the neutrino signals are quite sensitive to the difference of equation of state and can be used as a useful probe into the properties of dense matter. The event will be unique in that they will be shining only by neutrinos (and, possibly, gravitational waves) but not by photons, and hence they should be an important target of neutrino astronomy.

Pop III stars are the first stars in the universe. They do not contain metals and their formation and evolution may be different from that of stars of later generations. In fact, according to the theory of star formation, Pop III stars might have very massive components ($\sim 100 - 10000M_\odot$). In this paper, we compute the spherically symmetric gravitational collapse of these Pop III massive stars. We solve the general relativistic hydrodynamics and neutrino transfer equations simultaneously, treating neutrino reactions in detail. Unlike supermassive stars ($\gtrsim 10^5 M_\odot$), the stars of concern in this paper become opaque to neutrinos. The collapse is simulated until after an apparent horizon is formed. We confirm that the neutrino transfer plays a crucial role in the dynamics of gravitational collapse, and find also that the $\beta$-equilibration leads to a somewhat unfamiliar evolution of electron fraction. Contrary to the naive expectation, the neutrino spectrum does not become harder for more massive stars. This is mainly because the neutrino cooling is more efficient and the outer core is more massive as the stellar mass increases. Here the outer core is the outer part of the iron core falling supersonically. We also evaluate the flux of relic neutrino from Pop III massive stars. As expected, the detection of these neutrinos is difficult for the currently operating detectors. However, if ever observed, the spectrum will enable us to obtain the information on the formation history of Pop III stars. We investigate 18 models covering the mass range of $300 - 10^4 M_\odot$, making this study the most detailed numerical exploration of spherical gravitational collapse of Pop III massive stars. This will also serve as an important foundation for multi-dimensional investigations.

We discuss effects of magnetic fields on proto-neutron star winds by performing numerical simulation. We assume that the atmosphere of proto-neutron star has a homogenous magnetic field (ranging from ~10^{12} G to ~10^{15} G) perpendicular to the radial direction and examine the dependence of the three key quantities (dynamical time scale, electron fraction, and entropy per baryon) for the successful r-process on the magnetic field strength. Our results show that even with a magneter-class field strength, ~10^{15} G, the feature of the wind dynamics varies only little from that of non-magnetic winds, and that the condition for successful r-process is not realized.

We study the evolution of supernova core from the beginning of gravitational collapse of a 15Msolar star up to 1 second after core bounce. We present results of spherically symmetric simulations of core-collapse supernovae by solving general relativistic neutrino-radiation-hydrodynamics in the implicit time-differencing. We aim to explore the evolution of shock wave in a long term and investigate the formation of protoneutron star together with supernova neutrino signatures. These studies are done to examine the influence of equation of state (EOS) on the postbounce evolution of shock wave in the late phase and the resulting thermal evolution of protoneutron star. We make a comparison of two sets of EOS, that is, by Lattimer and Swesty (LS-EOS) and by Shen et al.(SH-EOS). We found that, for both EOSs, the core does not explode and the shock wave stalls similarly in the first 100 milliseconds after bounce. The revival of shock wave does not occur even after a long period in either cases. However, the recession of shock wave appears different beyond 200 milliseconds after bounce, having different thermal evolution of central core. A more compact protoneutron star is found for LS-EOS than SH-EOS with a difference in the central density by a factor of ~2 and a difference of ~10 MeV in the peak temperature. Resulting spectra of supernova neutrinos are different to the extent that may be detectable by terrestrial neutrino detectors.

We study the explosion mechanism of collapse-driven supernovae by numerical simulations with a new nuclear EOS based on unstable nuclei. We report new results of simulations of general relativistic hydrodynamics together with the Boltzmann neutrino-transport in spherical symmetry. We adopt the new data set of relativistic EOS and the conventional set of EOS (Lattimer-Swesty EOS) to examine the influence on dynamics of core-collapse, bounce and shock propagation. We follow the behavior of stalled shock more than 500 ms after the bounce and compare the evolutions of supernova core.

We perform a series of two-dimensional, axisymmetric, magnetohydrodynamic simulations of the rotational collapse of a supernova core. In order to calculate the waveforms of the gravitational wave, we derive the quadrupole formula including the contributions from the electromagnetic fields. Recent stellar evolution calculations imply that the magnetic fields of the toroidal components are much stronger than those of the poloidal ones at the presupernova stage. Thus, we systematically investigate the effects of the toroidal magnetic fields on the amplitudes and waveforms. Furthermore, we employ the two kinds of the realistic equation of states, which are often used in the supernova simulations. Then, we investigate the effects of the equation of states on the gravitational wave signals. With these computations, we find that the peak amplitudes are lowered by an order of 10% for the models with the strongest toroidal magnetic fields. However, the peak amplitudes are mostly within sensitivity range of laser interferometers such as TAMA and the first LIGO for a source at a distance of 10 kpc. Furthermore, we point out that the amplitudes of second peaks are still within the detection limit of the first LIGO for the source, although the characteristics of second peaks are reduced by the magnetic fields. We stress the importance of the detection, since it will give us information about the angular momentum distribution of massive evolved stars. When we compare the gravitational waves from the two realistic equation of states, significant differences are not found, except that the typical frequencies of the gravitational wave become slightly higher for the softer equation of state.

We study the nucleosyntheses of the light elements 7Li and 11B and the r-process elements in Type II supernovae from the point of view of supernova neutrinos and Galactic chemical evolution. We investigate the influence of the luminosity and average energy (temperature) of supernova neutrinos on these two nucleosynthesis processes. Common models of the total energy E_nu and the decay time of the neutrino luminosity tau_nu of supernova neutrinos are adopted to understand both processes. We adopt the model of the supernova explosion of a 16.2 M_odot star, which corresponds to SN 1987A, and calculate the nucleosynthesis of the light elements by the postprocessing. We find that the ejected masses of 7Li and 11B are roughly proportional to the total neutrino energy and weakly dependent on the decay time of the neutrino luminosity. As for the r-process nucleosynthesis, we adopt the same models of the neutrino luminosity in the neutrino-driven wind models of a 1.4 M_odot neutron star. We find that the r-process nucleosynthesis is affected through the peak neutrino luminosity, which depends on E_nu/tau_nu. We also discuss the unresolved problem of the overproduction of 11B in the Galactic chemical evolution of the light elements. We have to assume E_{nu} < 1.2 \times 10^{53} ergs to avoid the overproduction of 11B, which is too small to accept in comparison to 3.0 \times 10^{53} ergs deduced from the observation of SN1987A. We propose to reduce the temperature of nu_{mu, tau} and their anti-particles to T_{nu_{mu, tau}} = T_{bar{nu}_{mu,tau}} = 6.0 MeV/k. This modification of the neutrino temperature is shown to resolve the overproduction problem of 11B while still keeping a successful r-process abundance pattern.

We study characteristics of the relativistic equation of state (EOS) for collapse-driven supernovae, which is derived by relativistic nuclear many body theory. Recently the relativistic EOS table has become available as a new complete set of physical EOS for numerical simulations of supernova explosion. We examine this EOS table by using general relativistic hydrodynamics of the gravitational collapse and bounce of supernova cores. In order to study dense matter in dynamical situation, we perform simplified calculations of core collapse and bounce by following adiabatic collapse with the fixed electron fraction for a series of progenitor models. This is intended to give us ``approximate models'' of prompt explosion. We investigate the profiles of thermodynamical quantities and the compositions during collapse and bounce. We also perform the calculations with the Lattimer-Swesty EOS to compare the properties of dense matter. As a measure of the stiffness of the EOS, we examine the explosion energy of the prompt explosion with electron capture totally suppressed. We study the derivative of the thermodynamical quantities obtained by the relativistic EOS to discuss the convective condition in neutron-rich environment, which may be important in the delayed explosion.

We study the astrophysical reaction rate for the formation of $^{9}$Be through the three body reaction $\alpha(\alpha n,\gamma)$. This reaction is one of the key reactions which could bridge the mass gap at A = 8 nuclear systems to produce intermediate-to-heavy mass elements in alpha- and neutron-rich environments such as r-process nucleosynthesis in supernova explosions, s-process nucleosynthesis in asymptotic giant branch (AGB) stars, and primordial nucleosynthesis in baryon inhomogeneous cosmological models. To calculate the thermonuclear reaction rate in a wide range of temperatures, we numerically integrate the thermal average of cross sections assuming a two-steps formation through a metastable $^{8}$Be. Off-resonant and on-resonant contributions from the ground state in $^{8}$Be are taken into account. As input cross section, we adopt the latest experimental data by photodisintegration of $^{9}$Be with laser-electron photon beams, which covers all relevant resonances in $^{9}$Be. We provide the reaction rate for $\alpha(\alpha n,\gamma)^{9}$Be in the temperature range from T$_{9}$=10$^{-3}$ to T$_{9}$=10$^{1}$ both in the tabular form and in the analytical form. The calculated reaction rate is compared with the reaction rates of the CF88 and the NACRE compilations. The CF88 rate is valid at $T_{9} > 0.028$ due to lack of the off-resonant contribution. The CF88 rate differs from the present rate by a factor of two in a temperature range $T_{9} \geq 0.1$. The NACRE rate, which adopted different sources of experimental information on resonance states in $^{9}$Be, is 4--12 times larger than the present rate at $T_{9} \leq 0.028$, but is consistent with the present rate to within $\pm 20 %$ at $T_{9} \geq 0.1$.

We investigate the liquid-gas phase transition of dense matter in supernova explosion by the relativistic mean field approach and fragment based statistical model. The boiling temperature is found to be high (T_{boil} >= 0.7 MeV for rho_B >= 10^{-7} fm^{-3}), and adiabatic paths are shown to go across the boundary of coexisting region even with high entropy. This suggests that materials experienced phase transition can be ejected to outside. We calculated fragment mass and isotope distribution around the boiling point. We found that heavy elements at the iron, the first, second, and third peaks of r-process are abundantly formed at rho_B = 10^{-7}, 10^{-5}, 10^{-3} and 10^{-2} fm^{-3}, respectively.

We study the effects of the outer boundary conditions in neutrino-driven winds on the r-process nucleosynthesis. We perform numerical simulations of hydrodynamics of neutrino-driven winds and nuclear reaction network calculations of the r-process. As an outer boundary condition of hydrodynamic calculations, we set a pressure upon the outermost layer of the wind, which is approaching toward the shock wall. Varying the boundary pressure, we obtain various asymptotic thermal temperature of expanding material in the neutrino-driven winds for resulting nucleosynthesis. We find that the asymptotic temperature slightly lower than those used in the previous studies of the neutrino-driven winds can lead to a successful r-process abundance pattern, which is in a reasonable agreement with the solar system r-process abundance pattern even for the typical proto-neutron star mass Mns ~ 1.4 Msun. A slightly lower asymptotic temperature reduces the charged particle reaction rates and the resulting amount of seed elements and lead to a high neutron-to-seed ratio for successful r-process. This is a new idea which is different from the previous models of neutrino-driven winds from very massive (Mns ~ 2.0 Msun) and compact (Rns ~ 10 km) neutron star to get a short expansion time and a high entropy for a successful r-process abundance pattern. Although such a large mass is sometimes criticized from observational facts on a neutron star mass, we dissolve this criticism by reconsidering the boundary condition of the wind. We also explore the relation between the boundary condition and neutron star mass, which is related to the progenitor mass, for successful r-process.

We propose a new nucleosynthesis process in which various nuclei are formed through the liquid-gas phase transition of supernova matter as a preprocess of the standard r-process, and study its possibility of realization qualitatively. In the relativistic mean field and statistical model calculations, we find that this process may take place and help the later r-process to proceed.

We reanalyze $r$-process nucleosynthesis in the neutron-rich ejecta from a prompt supernova explosion of a low-mass (11 M$_\odot$) progenitor. A pompt explosion is not yet ruled out as a possibility for low-mass supernova progenitors. Moreover, there is mounting evidence that a new $r$-process site may be required. Hence, we assume that a prompt explosion can occur and make a study of r-process nucleosynthesis in the supernova ejecta. To achieve a prompt explosion we have performed a general relativistic hydrodynamic simulation of adiabatic collapse and bounce using a relativistic nuclear-matter equation of state. The electron fraction $Y_e$ during the collapse was fixed at the initial-model value. The size of the inner collapsing core was then large enough to enable a prompt explosion to occur. Adopting the calculated trajectories of promptly ejected material, we explicitly computed the burst of neutronization due to electron captures on free protons in the photodissociated ejecta after the passage of the shock. The thermal and compositional evolution of the resulting neutron-rich ejecta originating from near the surface of proto-neutron star was obtained. These were used in nuclear reaction network calculations to evaluate the products of $r$-process nucleosynthesis. We find that, unlike earlier studies, the amount of $r$-process material ejected per supernova is quite consistent with observed galactic $r$-process abundances. Furthermore, the computed $r$-process abundances are in good agreement with Solar abundances of $r$-process elements for A$> 100$. This suggests that prompt supernovae are still a viable $r$-process site. Such events may be responsible for the abundances of the heaviest $r$-process nuclei.

We study the role of light neutron-rich nuclei during r-process nucleosynthesis in supernovae. Most previous studies of the r-process have concentrated on the reaction flow of heavy unstable nuclei. Although the nuclear reaction network includes a few thousand heavy nuclei, only limited reaction flow through light-mass nuclei near the stability line has been used in those studies. However, in a viable scenario of the r-process in neutrino-driven winds, the initial condition is a high-entropy hot plasma consisting of neutrons, protons, and electron-positron pairs experiencing an intense flux of neutrinos. In such environments light-mass nuclei as well as heavy nuclei are expected to play important roles in the production of seed nuclei and r-process elements. Thus, we have extended our fully implicit nuclear reaction network so that it includes all nuclei up to the neutron drip line for Z $ \leq 10$, in addition to a larger network for Z $ \geq 10$. In the present nucleosynthesis study, we utilize a wind model of massive SNeII explosions to study the effects of this extended network. We find that a new nuclear-reaction flow path opens in the very light neutron-rich region. This new nuclear reaction flow can change the final heavy-element abundances by as much as an order of magnitude.

We study the neutrino-driven wind from the proto-neutron star by the general relativistic hydrodynamical simulations. We examine the properties of the neutrino-driven wind to explore the possibility of the r-process nucleosynthesis. The numerical simulations with the neutrino heating and cooling processes are performed with the assumption of the constant neutrino luminosity by using realistic profiles of the proto-neutron star (PNS) as well as simplified models. The dependence on the mass of PNS and the neutrino luminosity is studied systematically. Comparisons with the analytic treatment in the previous studies are also done. In the cases with the realistic PNS, we found that the entropy per baryon and the expansion time scale are neither high nor short enough for the r-process within the current assumptions. On the other hand, we found that the expansion time scale obtained by the hydrodynamical simulations is systematically shorter than that in the analytic solutions due to our proper treatment of the equation of state. This fact might lead to the increase of the neutron-to-seed ratio, which is suitable for the r-process in the neutrino-driven wind. Indeed, in the case of massive and compact proto-neutron stars with high neutrino luminosities, the expansion time scale is found short enough in the hydrodynamical simulations and the r-process elements up to A ~ 200 are produced in the r-process network calculation.

We construct the equation of state (EOS) of nuclear matter at finite temperature and density with various proton fractions within the relativistic mean field (RMF) theory for the use in the supernova simulations. The Thomas-Fermi approximation is adopted to describe the non-uniform matter where we consider nucleus, alpha-particle, proton and neutron in equilibrium. We treat the uniform matter and non-uniform matter consistently using the RMF theory. We tabulate the outcome as the pressure, free energy, entropy etc, with enough mesh points in wide ranges of the temperature, proton fraction, and baryon mass density.

We construct the equation of state (EOS) of nuclear matter using the relativistic mean field (RMF) theory in the wide density, temperature range with various proton fractions for the use of supernova simulation and the neutron star calculations. We first construct the EOS of homogeneous nuclear matter. We use then the Thomas-Fermi approximation to describe inhomogeneous matter, where heavy nuclei are formed together with free nucleon gas. We discuss the results on free energy, pressure and entropy in the wide range of astrophysical interest. As an example, we apply the resulting EOS on the neutron star properties by using the Oppenheimer-Volkoff equation.

First results of numerical simulations are presented which compute the dynamical evolution of a neutron star with a mass slightly below the minimum stable mass by means of a new implicit (general relativistic) hydrodynamic code. We show that such a star first undergoes a phase of quasi-static expansion, caused by slow nuclear $\beta$-decays, lasting for about 20 seconds, but then explodes violently. The kinetic energy of the explosion is around $10^{49} erg$, the peak luminosity in electron anti-neutrinos is of order $10^{52} erg/s$, and the thermodynamic conditions of the expanding matter are favorable for r-process nucleosynthesis. These results are obtained for the Harrison-Wheeler equation of state and a simple and, possibly, unrealistic treatment of $\beta$-decay rates and nuclear fission, which were adopted for comparison with previous works. However, we do not expect that the outcome will change qualitatively if more recent nuclear input physics used. Although our study does not rely on a specific scenario ofhow a neutron star starting from a bigger (and stable) masscan reach the dynamical phase, we assume that the final mass-loss event happens on a very short time scale, i.e., on a time scale shorter than a sound-crossing time, by removing a certain amount of mass as an initial perturbation. This assumption implies that the star has no time to adjust its nuclear composition to the new mass mass through a sequence of quasi-equilibria. In the latter case, however, there exists no stable configuration below the minimum mass, because the equation of state of fully catalyzed matter is too soft. Therefore, the dynamics of the explosion will not be too different from what have obtained if different initial perturbations are assumed.

We study the influence of the symmetry energy of the equation of state on the thermal evolution of protoneutron stars and the properties of supernova neutrinos by the numerical simulations after the protoneutron stars are formed. As for the equation of state (EOS) of nuclear matter, we take two EOS's with different symmetry energies obtained by the relativistic mean field theory. We find the symmetry energy plays the essential role on the evolution of lepton profiles and the neutrino fluxes.