Neutron stars are formed in core-collapse supernova explosions, where a large number of neutrinos are emitted. In this paper, supernova neutrino light curves are computed for the cooling phase of protoneutron stars, which lasts a few minutes. In the numerical simulations, 90 models of the phenomenological equation of state with different incompressibilities, symmetry energies, and nucleon effective masses are employed for a comprehensive study. It is found that the cooling timescale is longer for a model with a larger neutron star mass and a smaller neutron star radius. Furthermore, a theoretical expression of the cooling timescale is presented as a function of the mass and radius and it is found to describe the numerical results faithfully. These findings suggest that diagnosing the mass and radius of a newly formed neutron star using its neutrino signal is possible.

Whether fast cooling processes occur or not is crucial for the thermal evolution of neutron stars. In particular, the threshold of the direct Urca process, which is one of the fast cooling processes, is determined by the interior proton fraction $Y_p$, or the nuclear symmetry energy. Since recent observations indicate the small radius of neutron stars, a low value is preferred for the symmetry energy. In this study, simulations of neutron star cooling are performed adopting three models for equation of state (EoS): Togashi, Shen, and LS220 EoSs. The Togashi EoS has been recently constructed with realistic nuclear potentials under finite temperature, and found to account for the small radius of neutron stars. As a result, we find that, since the direct Urca process is forbidden, the neutron star cooling is slow with use of the Togashi EoS. This is because symmetry energy of Togashi EoS is lower than those of other EoSs. Hence, in order to account for observed age and surface temperature of isolated neutron stars (INS) with use of the Togashi EoS, other fast cooling processes are needed regardless of the surface composition.

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.

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.

The cooling process of a protoneutron star is investigated with focus on its sensitivity to properties of hot and dense matter. An equation of state, which includes the nucleon effective mass and nuclear symmetry energy at twice the saturation density as control parameters, is constructed for systematic studies. The numerical code utilized in this study follows a quasi-static evolution of a protoneutron star solving the general-relativistic stellar structure with neutrino diffusion. The cooling timescale evaluated from the neutrino light curve is found to be longer for the models with larger effective masses and smaller symmetry energies at high densities. The present results are compared with those for other equations of state and it is found that they are consistent in terms of their dependences on the effective mass and neutron star radius.

Event spectra of the neutrino-$^{16}$O charged-current reactions in Super-Kamiokande are evaluated for a future supernova neutrino burst. Since these channels are expected to be useful for diagnosing a neutrino spectrum with high average energy, the evaluations are performed not only for an ordinary supernova neutrino model but also for a model of neutrino emission from a black-hole-forming collapse. Using shell model results, whose excitation energies are consistent with the experimental data, the cross sections of the $^{16}$O($\nu_e, e^-$)X and $^{16}$O($\bar\nu_e, e^+$)X reactions for each nuclear state with a different excitation energy are employed in this study. It is found that, owing to the components of the reaction with higher excitation energy, the event spectrum becomes 4-7 MeV softer than that in the case without considering the excitation energies. In addition, a simplified approach to evaluate the event spectra is proposed for convenience and its validity is examined.

A proto-neutron star (PNS) is a newly formed compact object in a core collapse supernova. In this Letter, the neutrino emission from the cooling process of a PNS is investigated using two types of nuclear equation of state (EOS). It is found that the neutrino signal is mainly determined by the high-density EOS. The neutrino luminosity and mean energy are higher and the cooling time scale is longer for the softer EOS. Meanwhile, the neutrino mean energy and the cooling time scale are also affected by the low-density EOS because of the difference in the population of heavy nuclei. Heavy nuclei have a large scattering cross section with neutrinos owing to the coherent effects and act as thermal insulation near the surface of a PNS. The neutrino mean energy is higher and the cooling time scale is longer for an EOS with a large symmetry energy at low densities, namely a small density derivative coefficient of the symmetry energy, $L$.

A new table of the nuclear equation of state (EOS) based on realistic nuclear potentials is constructed for core-collapse supernova numerical simulations. Adopting the EOS of uniform nuclear matter constructed by two of the present authors with the cluster variational method starting from the Argonne v18 and Urbana IX nuclear potentials, the Thomas-Fermi calculation is performed to obtain the minimized free energy of a Wigner-Seitz cell in non-uniform nuclear matter. As a preparation for the Thomas-Fermi calculation, the EOS of uniform nuclear matter is modified so as to remove the effects of deuteron cluster formation in uniform matter at low densities. Mixing of alpha particles is also taken into account following the procedure used by Shen et al. (1998, 2011). The critical densities with respect to the phase transition from non-uniform to uniform phase with the present EOS are slightly higher than those with the Shen EOS at small proton fractions. The critical temperature with respect to the liquid-gas phase transition decreases with the proton fraction in a more gradual manner than in the Shen EOS. Furthermore, the mass and proton numbers of nuclides appearing in non-uniform nuclear matter with small proton fractions are larger than those of the Shen EOS. These results are consequences of the fact that the density derivative coefficient of the symmetry energy of our EOS is smaller than that of the Shen EOS.

The cosmological evolution of the binary black hole (BH) merger rate and the energy density of the gravitational-wave (GW) background are investigated. To evaluate the redshift dependence of the BH formation rate, BHs are assumed to originate from low-metallicity stars, and the relations between the star formation rate, metallicity and stellar mass of galaxies are combined with the stellar mass function at each redshift. As a result, it is found that when the energy density of the GW background is scaled with the merger rate at the local Universe, the scaling factor does not depend on the critical metallicity for the formation of BHs. Also taking into account the merger of binary neutron stars, a simple formula to express the energy spectrum of the GW background is constructed for the inspiral phase. The relation between the local merger rate and the energy density of the GW background will be examined by future GW observations.

The spectrum of the supernova relic neutrino (SRN) background from past stellar collapses including black hole formation (failed supernovae) is calculated. The redshift dependence of the black hole formation rate is considered on the basis of the metallicity evolution of galaxies. Assuming the mass and metallicity ranges of failed supernova progenitors, their contribution to SRNs is quantitatively estimated for the first time. Using this model, the dependences of SRNs on the cosmic star formation rate density, shock revival time and equation of state are investigated. The shock revival time is introduced as a parameter that should depend on the still unknown explosion mechanism of core collapse supernovae. The dependence on equation of state is considered for failed supernovae, whose collapse dynamics and neutrino emission are certainly affected. It is found that the low-energy spectrum of SRNs is mainly determined by the cosmic star formation rate density. These low-energy events will be observed in the Super-Kamiokande experiment with gadolinium-loaded water.

A new cellular automaton (CA) model is presented for the self-organized criticality (SOC) in recurrent bursts of soft gamma repeaters (SGRs), which are interpreted as avalanches of reconnection in the magnetosphere of neutron stars. The nodes of a regular dodecahedron and a truncated icosahedron are adopted as spherically closed grids enclosing a neutron star. It is found that the system enters the SOC state if there are sites where the expectation value of the added perturbation is nonzero. The energy distributions of SOC avalanches in CA simulations are described by a power law with a cutoff, which is consistent with the observations of SGR 1806-20 and SGR 1900+14. The power-law index is not universal and depends on the amplitude of the perturbation. This result shows that the SOC of SGRs can be illustrated not only by the crust quake model but also by the magnetic reconnection model.

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.

The spectrum and event rate of supernova relic neutrinos are calculated taking into account the dependence on the time it takes for the shock wave in supernova cores to revive. The shock revival time should depend on the still unknown explosion mechanism of collapse-driven supernovae. The contribution of black-hole-forming failed supernovae is also considered. The total event rate is higher for models with a longer shock revival time and/or a failed-supernova contribution. The hardness of the spectrum does not strongly depend on the shock revival time, but the spectrum becomes hard owing to the failed supernovae. Therefore, the shock-revival-time dependence of supernova relic neutrinos has different systematics from the fractions of failed supernovae.

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.

The equation of state for neutron stars in a wide-density range at zero temperature is constructed. The chiral quark-meson coupling model within relativistic Hartree-Fock approximation is adopted for uniform nuclear matter. The coupling constants are determined so as to reproduce the experimental data of atomic nuclei and hypernuclei. In the crust region, nuclei are taken into account within the Thomas-Fermi calculation. All octet baryons are considered in the core region, while only $\Xi^{-}$ appears in neutron stars. The resultant maximum mass of neutron stars is $1.95M_\odot$, which is consistent with the constraint from the recently observed massive pulsar, PSR J1614-2230.

We systematically examine the fundamental frequencies of shear torsional oscillations in neutron star crusts in a manner that is dependent on the parameter $L$ characterizing the poorly known density dependence of the symmetry energy. The identification of the lowest quasiperiodic oscillation (QPO) among the observed QPOs from giant flares in soft-gamma repeaters as the $/ell=2$ fundamental torsional oscillations enables us to constrain the parameter $L$ as $L/ge 47.4$ MeV, which is the most conservative restriction on $L$ obtained in the present work that assumes that the mass and radius of the flaring neutron stars range 1.4-1.8 $M_{/odot}$ and 10-14 km. Next, we identify one by one a set of the low-lying frequencies observed in giant flares as the fundamental torsional oscillations. The values of $L$ that can reproduce all the observed frequencies in terms of the torsional oscillations coupled with a part of dripped neutrons via entrainment effects are then constrained as 101.1 MeV $/le L /le$ 131.0 MeV. Alternatively, if only the second lowest frequency observed in SGR 1806$-$20 has a different origin, one obtains relatively low $L$ values ranging 58.0 MeV $/le L /le$ 85.3 MeV, which seem more consistent with other empirical constraints despite large uncertainties.

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.

We consider how superfluidity of dripped neutrons in the crust of a neutron star affects the frequencies of the crust's fundamental torsional oscillations. A nonnegligible superfluid part of dripped neutrons, which do not comove with nuclei, act to reduce the enthalpy density and thus enhance the oscillation frequencies. By assuming that the quasi-periodic oscillations observed in giant flares of soft gamma repeaters arise from the fundamental torsional oscillations and that the mass and radius of the neutron star is in the range of 1.4 < M/M_SUN < 1.8 and 10 km < R < 14 km, we constrain the density derivative of the symmetry energy as 100 MeV < L < 130 MeV, which is far severer than the previous one, L > 50 MeV, derived by ignoring the superfluidity.

We general relativistically calculate the frequency of fundamental torsional oscillations of neutron star crusts, where we focus on the crystalline properties obtained from macroscopic nuclear models in a way depending on the equation of state of nuclear matter. We find that the calculated frequency is sensitive to the density dependence of the symmetry energy, but almost independent of the incompressibility of symmetric nuclear matter. By identifying the lowest-frequency quasi-periodic oscillation in giant flares observed from soft gamma-ray repeaters as the fundamental torsional mode and allowing for the dependence of the calculated frequency on stellar models, we provide a lower limit of the density derivative of the symmetry energy as $L\simeq 50$ MeV.

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.

In supernova cores and neutron star crusts, nuclei are thought to deform to rodlike and slablike shapes, which are often called nuclear pasta. We study the equilibrium properties of the nuclear pasta by using a liquid drop model with curvature corrections. It is confirmed that the curvature effect acts to lower the transition densities between different shapes. We also examine the gyroid structure, which was recently suggested as a different type of nuclear pasta by analogy with the polymer systems. The gyroid structure investigated in this paper is approximately formulated as an extension of the periodic minimal surface whose mean curvature vanishes. In contrast to our expectations, we find from the present approximate formulation that the curvature corrections act to slightly disfavor the appearance of the gyroid structure. By comparing the energy corrections in the gyroid phase and the hypothetical phases composed of d-dimensional spheres, where d is a general dimensionality, we show that the gyroid is unlikely to belong to a family of the generalized dimensional spheres.

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.

Nuclear matter is considered to be inhomogeneous at subnuclear densities that are realized in supernova cores and neutron star crusts, and the structures of nuclear matter change from spheres to cylinders, slabs, cylindrical holes and spherical holes as the density increases. In this letter, we discuss other possible structures, that is, gyroid and double-diamond morphologies, which are periodic bicontinuous structures discovered in a block copolymer. Utilizing the compressible liquid drop model, we evaluate their surface and Coulomb energies and show that there is a chance of gyroid appearance near the transition point from a cylinder to a slab. This interesting analogy between nuclear and polymer systems is not merely qualitative. The volume fraction at the phase transition is also similar for the two systems. Although the five shapes listed initially have been long thought to be the only major constituents of so-called nuclear pasta at subnuclear densities, our findings imply that this may not be the case and suggest that more detailed studies on nuclear pasta including the gyroid phase are needed.

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 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.

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.

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.