BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kohsuke Sumiyoshi

Kohsuke Sumiyoshi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEnucl-thastro-ph.SRhep-phastro-phastro-ph.COgr-qchep-thphysics.comp-ph

Catalog footprint

What is connected

35works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2022arXiv

Effects of nuclear matter and composition in core-collapse supernovae and long-term proto-neutron star cooling

We study the influence of hot and dense matter in core-collapse supernovae by adopting up-to-date nuclear equation of state (EOS) based on the microscopic nuclear many-body frameworks. We explore effects of EOS based on the Dirac Brueckner Hartree-Fock theory through comparisons with those based on the variational method. We also examine effects of the differences in the composition of nuclei and nucleons by using the same EOS by the variational method but employing two different treatments in computations of nuclear abundances. We perform numerical simulations of core-collapse supernovae adopting the three EOSs. We also perform numerical simulations of the long-term evolution over 70 s of the proto-neutron star cooling. We show that impacts by different modeling of composition are remarkable as in those by different treatments of uniform matter in the gravitational collapse, bounce, and shock propagation. The cooling of proto-neutron star and the resulting neutrino emission are also affected by the compositional difference even if the same treatment in computing uniform matter of EOS.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Neutrino emission from the collapse of $\sim 10^4$ $M_\odot$ population III supermassive stars

We calculate the neutrino signal from Population III supermassive star collapse using a neutrino transfer code originally developed for core collapse supernovae and massive star collapse. Using this code, we are able to investigate the supermassive star mass range thought to undergo neutrino trapping ($\sim 10^4$ M$_\odot$), a mass range which has been neglected by previous works because of the difficulty of neutrino transfer. For models in this mass range, we observe a neutrino-sphere with a large radius and low density compared to typical massive star neutrino-spheres. We calculate the neutrino light-curve emitted from this neutrino-sphere. The resulting neutrino luminosity is significantly lower than the results of a previous analytical model. We briefly discuss the possibility of detecting a neutrino burst from a supermassive star or the neutrino background from many supermassive stars and conclude that the former is unlikely with current technology, unless the SMS collapse is located as close as 1 Mpc, while the latter is also unlikely even under very generous assumptions. However, the supermassive star neutrino background is still of interest as it may serve as a source of noise in proposed dark matter direct detection experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Observing Supernova Neutrino Light Curves with Super-Kamiokande. II. Impact of the Nuclear Equation of State

The late-time evolution of the neutrino event rate from supernovae is evaluated for Super-Kamiokande using simulated results of proto-neutron star (PNS) cooling. In the present work we extend the result of Suwa et al. (2019) [arXiv:1904.09996], which studied the dependence on the PNS mass, but focus on the impact of the nuclear equation of state (EOS). We find that the neutrino event rate depends on both the high-density and low-density EOS, where the former determines the radius of the PNS and the latter affects its surface temperature. Based on the present evaluation of the neutrino event rate we propose a new analysis method to extract the time variability of the neutrino average energy taking into account the statistical error in the observation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Observing Supernova Neutrino Light Curves with Super-Kamiokande. III. Extraction of Mass and Radius of Neutron Stars from Synthetic Data

Neutrinos are guaranteed to be observable from the next galactic supernova (SN). Optical light and gravitational waves are also observable, but may be difficult to observe if the location of the SN in the galaxy or the details of the explosion are unsuitable. The key to observing the next supernova is to first use neutrinos to understand various physical quantities and then link them to other signals. In this paper, we present Monte Carlo sampling calculations of neutrino events from galactic supernova explosions observed with Super-Kamiokande. The analytical solution of neutrino emission, which represents the long-term evolution of neutrino-light curve from supernovae, is used as a theoretical template. It gives the event rate and event spectrum through inverse beta decay interactions with explicit model parameter dependence. Parameter estimation is performed on these simulated sample data by fitting least squares using the analytical solution. The results show that the mass, radius and total energy of a remnant neutron star produced by a SN can be determined with an accuracy of $\sim 0.1M_\odot$, $\sim 1$ km, and $\sim 10^{51}$ erg, respectively, for a galactic SN at 8 kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Principal-Axis Analysis of the Eddington Tensor for the Early Post-Bounce Phase of Rotational Core-Collapse Supernovae

Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. We employed the results obtained in the early post-bounce phase (t < 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment method. We first made a comparison in 1D under spherical symmetry, though, of the Eddington factor p defined in the fluid rest frame (FR). We confirmed that the maximum entropy closure for the Fermionic distribution (MEFD) performs better than others near the proto-neutron star surface, where p < 1/3 occurs, but does not work well even in 1D when the phase space occupancy satisfies e < 0.5 together with p < 1/3, the condition known to be not represented by MEFD. For the 2D models with the rapid rotation, we employed the principal-axis analysis of the Eddington tensor. We paid particular attention to the direction of the longest principal axis. We observed in FR that it is aligned neither with the radial direction nor with the neutrino flux in 2D, particularly so in convective and/or rapidly rotating regions, the fact not accommodated in the moment method. We repeated the same analysis in the laboratory frame (LB) and found again that the direction of the longest principal axis is not well reproduced by MEFD because the interpolation between the optically thick and thin limits is not very accurate in this frame.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stability analysis of supermassive primordial stars: a new mass range for general relativistic instability supernovae

Observed supermassive black holes in the early universe have several proposed formation channels, in part because most of these channels are difficult to probe. One of the more promising channels, the direct collapse of a supermassive star, has several possible probes including the explosion of a helium-core supermassive star triggered by a general relativistic instability. We develop a straightforward method for evaluating the general relativistic radial instability without simplifying assumptions and apply it to population III supermassive stars taken from a post Newtonian stellar evolution code. This method is more accurate than previous determinations and it finds that the instability occurs earlier in the evolutionary life of the star. Using the results of the stability analysis, we perform 1D general relativistic hydrodynamical simulations and we find two general relativistic instability supernovae fueled by alpha capture reactions as well as several lower mass pulsations, analogous to the puslational pair instability process. The mass range for the events (2.6-3.0 $\times 10^4$ ${\rm M}_\odot$) is lower than had been suggested by previous works (5.5 $\times 10^4$ ${\rm M}_\odot$) because the instability occurs earlier in the star's evolution. The explosion may be visible to, among others, JWST, while the discovery of the pulsations opens up additional possibilities for observation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Analytic solutions for neutrino-light curves of core-collapse supernovae

Neutrinos are a guaranteed signal from supernova explosions in the Milky Way, and a most valuable messenger that can provide us with information about the deepest parts of supernovae. In particular, neutrinos will provide us with physical quantities, such as the radius and mass of protoneutron stars (PNS), which are the central engine of supernovae. This requires a theoretical model that connects observables such as neutrino luminosity and average energy with physical quantities. Here, we show analytic solutions for the neutrino-light curve derived from the neutrino radiation transport equation by employing the diffusion approximation and the analytic density solution of the hydrostatic equation for a PNS. The neutrino luminosity and the average energy as functions of time are explicitly presented, with dependence on PNS mass, radius, the total energy of neutrinos, surface density, and opacity. The analytic solutions provide good representations of the numerical models from a few seconds after the explosion and allow a rough estimate of these physical quantities to be made from observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Developing an end-to-end simulation framework of supernova neutrino detection

Massive stars can explode as supernovae at the end of their life cycle, releasing neutrinos whose total energy reaches $10^{53}$ erg. Moreover, neutrinos play key roles in supernovae, heating and reviving the shock wave as well as cooling the resulting protoneutron star. Therefore, neutrino detectors are waiting to observe the next galactic supernova and several theoretical simulations of supernova neutrinos are underway. While these simulation concentrate mainly on only the first one second after the supernova bounce, the only observation of a supernova with neutrinos, SN 1987A, revealed that neutrino emission lasts for more than 10 seconds. For this reason, long-time simulation and analysis tools are needed to compare theories with the next observation. Our study is to develop an integrated supernova analysis framework to prepare an analysis pipeline for treating galactic supernovae observations in the near future. This framework deals with the core-collapse, bounce and proto-neutron star cooling processes, as well as with neutrino detection on earth in a consistent manner. We have developed a new long-time supernova simulation in one dimension that explodes successfully and computes the neutrino emission for up to 20 seconds. Using this model we estimate the resulting neutrino signal in the Super-Kamiokande detector to be about 1,800 events for an explosion at 10 kpc and discuss its implications in this paper. We compare this result with the SN 1987A observation to test its reliability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of symmetry energy on equation of state for simulations of core-collapse supernovae and neutron-star mergers

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of symmetry energy on the radius and tidal deformability of neutron stars in relativistic mean-field model

The radii and tidal deformabilities of neutron stars are investigated in the framework of relativistic mean-field (RMF) model with different density-dependent behaviors of symmetry energy. To study the effects of symmetry energy on the properties of neutron stars, an $ω$ meson and $ρ$ meson coupling term is included in a popular RMF Lagrangian, i.e. the TM1 parameter set, which is used for the widely used supernova equation of state (EoS) table. The coupling constants relevant to the vector-isovector meson, $ρ$, are refitted by a fixed symmetry energy at subsaturation density and its slope at saturation density, while other coupling constants remain the same as the original ones in TM1 so as to update the supernova EoS table. The radius and mass of maximum neutron stars are not so sensitive to the symmetry energy in these family TM1 parameterizations. However, the radii at intermediate mass region are strongly correlated with the slope of symmetry energy. Furthermore, the dimensionless tidal deformabilities of neutron stars are also calculated within the associated Love number. We find that its value at $1.4 M_\odot$ has a linear correlation to the slope of symmetry energy being different from the previous studied. With the latest constraints of tidal deformabilities from GW170817 event, the slope of symmetry energy at nuclear saturation density should be smaller than $60$ MeV in the family TM1 parameterizations. This fact supports the usage of lower symmetry energy slope for the update supernova EoS, which is applicable to simulations of neutron star merger. Furthermore, the analogous analysis are also done within the family IUFSU parameter sets. It is found that the correlations between the symmetry energy slope with the radius and tidal deformability at $1.4 M_\odot$ have very similar linear relations in these RMF models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fast collective neutrino oscillations inside the neutrino sphere in core-collapse supernovae

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 $ν_e$ and $\barν_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 $ν_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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fast Neutrino Flavor Conversion Modes in Multidimensional Core-collapse Supernova Models: the Role of the Asymmetric Neutrino Distributions

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 $ν_e$ and $\barν_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 $ν_e - \barν_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 $ν_e - \barν_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 $ν_e$ and $\barν_e$. In addition, apart from the crossings within the neutrino decoupling region, there is a class of $ν_e-\barν_e$ angular crossings which appears very deep inside the proto-neutron star.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the occurrence of fast neutrino flavor conversions in multidimensional supernova models

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 $ν_e$ and $\barν_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 $ν_e$ and $\barν_e$ fluxes in multi-D models can assist the appearance of the crossing between the $ν_e$ and $\barν_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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Properties of neutrino transfer in a deformed remnant of neutron star merger

We study properties of neutrino transfer in a remnant of neutron star merger, consisting of a massive neutron star and a surrounding torus. We perform numerical simulations of the neutrino transfer by solving the Boltzmann equation with momentum-space angles and energies of neutrinos for snapshots of the merger remnant having elongated shapes. The evaluation of the neutrino distributions in the multi-dimensions enable us to provide the detailed information of angle and energy spectra and neutrino reaction rates. We demonstrate features of asymmetric neutrino fluxes from the deformed remnant and investigate the neutrino emission region by determining the neutrinosphere for each energy. We examine the emission and absorption of neutrinos to identify important ingredients of heating rates through neutrino irradiation. We show that the contributions of $μ$- and $τ$-types neutrinos are important for the heating in the region above the massive neutron star. We also examine the angle moments and the Eddington tensor calculated directly by the neutrino distribution functions and compare them with those obtained by a moment closure approach, which is often used in the study of neutrino-radiation hydrodynamics. We show that the components of the Eddington tensor have non-monotonic behaviors and the approximation of the closure relation may become inaccurate for high energy neutrinos, whose fluxes are highly aspherical due to the extended merger remnant.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Final Fate of Supermassive $M \sim 5 \times 10^4 \; M_\odot$ Pop III Stars: Explosion or Collapse?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Influence of density dependence of symmetry energy in hot and dense matter for supernova simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exact and approximate expressions of energy generation rates and their impact on the explosion properties of Pair Instability Supernovae

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

EOSDB: The Database for Nuclear EoS

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Three-dimensional Boltzmann-Hydro code for core-collapse in massive stars I. special relativistic treatments

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New equations of state based on the liquid drop model of heavy nuclei and quantum approach to light nuclei for core-collapse supernova simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

QCD phase diagram at finite baryon and isospin chemical potentials in Polyakov loop extended quark meson model with vector interaction

We investigate the QCD phase diagram of isospin asymmetric matter using the Polyakov loop extended quark meson (PQM) model with vector interaction. The critical point temperature is found to decrease in isospin asymmetric matter and disappear at large isospin chemical potential. We also discuss the QCD phase transition in the neutron star core. From comparison of the QCD phase diagram in PQM and corresponding baryon and isospin chemical potentials of neutron star matter in relativistic mean field models, we show that the order of the chiral phase transition in the neutron star core could be crossover because of large isospin chemical potential.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stellar Core Collapse with Hadron-Quark Phase Transition

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Supernova Neutrino Light Curves and Spectra for Various Progenitor Stars: From Core Collapse to Proto-neutron Star Cooling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability in core-collapse supernovae

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Core-Collapse Supernovae as Supercomputing Science: a status report toward 6D simulations with exact Boltzmann neutrino transport in full general relativity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Neutrino Transfer in Three Dimensions for Core-Collapse Supernovae. I. Static Configurations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A new baryonic equation of state at sub-nuclear densities for core-collapse simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hyperon Matter and Black Hole Formation in Failed Supernovae

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exploring Hadron Physics in Black Hole Formations: a New Promising Target of Neutrino Astronomy

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impact of Quarks and Pions on Dynamics and Neutrino Signal of Black Hole Formation in Non-rotating Stellar Core Collapse

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Oscillation and Future Detection of Failed Supernova Neutrinos from 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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Numerical Study on Stellar Core Collapse and Neutrino Emission: Probe into the Spherically Symmetric Black Hole Progenitors with 3 - 30Msun Iron Cores

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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Effects of Magnetic Fields on Proto-Neutron Star Winds

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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Gravitational Collapse and Neutrino Emission of Population III Massive Stars

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

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Gravitational Radiation from Rotational Core Collapse: Effects of Magnetic Fields and Realistic Equation of States

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.

0 citations0 reviewsNo code linkedOpen access