Towards generating a new supernova equation of state: A systematic analysis of cold hybrid stars

The hadron-quark phase transition in core-collapse supernovae (CCSNe) has the potential to trigger explosions in otherwise nonexploding models. However, those hybrid supernova equations of state (EOS) shown to trigger an explosion do not support the observational 2 M$_\odot$ neutron star maximum mass constraint. In this work, we analyze cold hybrid stars by the means of a systematic parameter scan for the phase transition properties, with the aim to develop a new hybrid supernova EOS. The hadronic phase is described with the state-of-the-art supernova EOS HS(DD2), and quark matter by an EOS with a constant speed of sound (CSS) of $c_{QM}^2=1/3$. We find promising cases which meet the 2 M$_\odot$ criterion and are interesting for CCSN explosions. We show that the very simple CSS EOS is transferable into the well-known thermodynamic bag model, important for future application in CCSN simulations. In the second part, the occurrence of reconfinement and multiple phase transitions is discussed. In the last part, the influence of hyperons in our parameter scan is studied. Including hyperons no change in the general behavior is found, except for overall lower maximum masses. In both cases (with and without hyperons) we find that quark matter with $c_{QM}^2=1/3$ can increase the maximum mass only if reconfinement is suppressed or if quark matter is absolutely stable.Comment: 14 pages, 11 figures, v2: matches published versio

The r-process nucleosynthesis in the various jet-like explosions of magnetorotational core-collapse supernovae

The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., prompt-magnetic-jet and delayed-magnetic-jet, for which the magnetic fields of proto-neutron stars (PNSs) during collapse and the core-bounce are strong and comparatively moderate, respectively. Following the hydrodynamical trajectories of each explosion model, we confirmed that r-processes successfully occur in the prompt-magnetic-jets, which produce heavy nuclei including actinides. On the other hand, the r-process in the delayed-magnetic-jet is suppressed, which synthesizes only nuclei up to the second peak ($A \sim 130$). Thus, the r-process in the delayed-magnetic-jets could explain only "weak r-process" patterns observed in metal-poor stars rather than the "main r-process", represented by the solar abundances. Our results imply that core-collapse supernovae are possible astronomical sources of heavy r-process elements if their magnetic fields are strong enough, while weaker magnetic explosions may produce "weak r-process" patterns ($A \lesssim 130$). We show the potential importance and necessity of magneto-rotational supernovae for explaining the galactic chemical evolution, as well as abundances of r-process enhanced metal-poor stars. We also examine the effects of the remaining uncertainties in the nature of PNSs due to weak interactions that determine the final neutron-richness of ejecta. Additionally, we briefly discuss radioactive isotope yields in primary jets (e.g., $^{56}$Ni), with relation to several optical observation of SNe and relevant high-energy astronomical phenomena.Comment: 25 pages, 17 figures, 6 tables, ApJ in press, numerical data are available at http://www.astro.keele.ac.uk/~nobuya/mrsn and https://github.com/nnobuya/mrs

Explosive hydrogen burning during type I X-ray bursts

Explosive hydrogen burning in type I X-ray bursts (XRBs) comprise charged particle reactions creating isotopes with masses up to A~100. Since charged particle reactions in a stellar environment are very temperature sensitive, we use a realistic time-dependent general relativistic and self-consistent model of type I x-ray bursts to provide accurate values of the burst temperatures and densities. This allows a detailed and accurate time-dependent identification of the reaction flow from the surface layers through the convective region and the ignition region to the neutron star ocean. Using this, we determine the relative importance of specific nuclear reactions in the X-ray burst.Comment: 53 pages, 24 figures, submitted to Astrophys.

Origin of the elements

What is the origin of the oxygen we breathe, the hydrogen and oxygen (in form of water H2O) in rivers and oceans, the carbon in all organic compounds, the silicon in electronic hardware, the calcium in our bones, the iron in steel, silver and gold in jewels, the rare earths utilized, e.g. in magnets or lasers, lead or lithium in batteries, and also of naturally occurring uranium and plutonium? The answer lies in the skies. Astrophysical environments from the Big Bang to stars and stellar explosions are the cauldrons where all these elements are made. The papers by Burbidge (Rev Mod Phys 29:547–650, 1957) and Cameron (Publ Astron Soc Pac 69:201, 1957), as well as precursors by Bethe, von Weizsäcker, Hoyle, Gamow, and Suess and Urey provided a very basic understanding of the nucleosynthesis processes responsible for their production, combined with nuclear physics input and required environment conditions such as temperature, density and the overall neutron/proton ratio in seed material. Since then a steady stream of nuclear experiments and nuclear structure theory, astrophysical models of the early universe as well as stars and stellar explosions in single and binary stellar systems has led to a deeper understanding. This involved improvements in stellar models, the composition of stellar wind ejecta, the mechanism of core-collapse supernovae as final fate of massive stars, and the transition (as a function of initial stellar mass) from core-collapse supernovae to hypernovae and long duration gamma-ray bursts (accompanied by the formation of a black hole) in case of single star progenitors. Binary stellar systems give rise to nova explosions, X-ray bursts, type Ia supernovae, neutron star, and neutron star–black hole mergers. All of these events (possibly with the exception of X-ray bursts) eject material with an abundance composition unique to the specific event and lead over time to the evolution of elemental (and isotopic) abundances in the galactic gas and their imprint on the next generation of stars. In the present review, we want to give a modern overview of the nucleosynthesis processes involved, their astrophysical sites, and their impact on the evolution of galaxies

Massive Stars and their Supernovae

Massive stars and their supernovae are prominent sources of radioactive isotopes, the observations of which thus can help to improve our astrophysical models of those. Our understanding of stellar evolution and the final explosive endpoints such as supernovae or hypernovae or gamma-ray bursts relies on the combination of magneto-hydrodynamics, energy generation due to nuclear reactions accompanying composition changes, radiation transport, and thermodynamic properties (such as the equation of state of stellar matter). Nuclear energy production includes all nuclear reactions triggered during stellar evolution and explosive end stages, also among unstable isotopes produced on the way. Radiation transport covers atomic physics (e.g. opacities) for photon transport, but also nuclear physics and neutrino nucleon/nucleus interactions in late phases and core collapse. Here we want to focus on the astrophysical aspects, i.e. a description of the evolution of massive stars and their endpoints, with a special emphasis on the composition of their ejecta (in form of stellar winds during the evolution or of explosive ejecta). Low and intermediate mass stars end their evolution as a white dwarf with an unburned C and O composition. Massive stars evolve beyond this point and experience all stellar burning stages from H over He, C, Ne, O and Si-burning up to core collapse and explosive endstages. In this chapter we discuss the nucleosynthesis processes involved and the production of radioactive nuclei in more detail.Comment: 79 pages; Chapter of "Astronomy with Radioactivities", a book in Springer's 'lecture notes in physics series, Vol. 812, Eds. Roland Diehl, Dieter H. Hartmann, and Nikos Prantzos, to appear in summer 201

A full general relativistic neutrino radiation-hydrodynamics simulation of a collapsing very massive star and the formation of a black hole

We study the final fate of a very massive star by performing full general relativistic (GR), three-dimensional (3D) simulation with three-flavor multi-energy neutrino transport. Utilizing a 70 solar mass zero metallicity progenitor, we self-consistently follow the radiation-hydrodynamics from the onset of gravitational core-collapse until the second collapse of the proto-neutron star (PNS), leading to black hole (BH) formation. Our results show that the BH formation occurs at a post-bounce time of ~300 ms for the 70 Msun star. This is significantly earlier than those in the literature where lower mass progenitors were employed. At a few ~10 ms before BH formation, we find that the stalled bounce shock is revived by intense neutrino heating from the very hot PNS, which is aided by violent convection behind the shock. In the context of 3D-GR core-collapse modeling with multi-energy neutrino transport, our numerical results present the first evidence to validate a fallback BH formation scenario of the 70 Msun star.Comment: 5 pages, 4 figures, accepted for publication in MNRAS lette

Two-Dimensional Core-Collapse Supernova Simulations with the Isotropic Diffusion Source Approximation for Neutrino Transport

The neutrino mechanism of core-collapse supernova is investigated via non-relativistic, two-dimensional (2D), neutrino radiation-hydrodynamic simulations. For the transport of electron flavor neutrinos, we use the interaction rates defined by Bruenn (1985) and the isotropic diffusion source approximation (IDSA) scheme, which decomposes the transported particles into trapped particle and streaming particle components. Heavy neutrinos are described by a leakage scheme. Unlike the "ray-by-ray" approach in some other multi-dimensional supernova models, we use cylindrical coordinates and solve the trapped particle component in multiple dimensions, improving the proto-neutron star resolution and the neutrino transport in angular and temporal directions. We provide an IDSA verification by performing 1D and 2D simulations with 15 and 20 $M_\odot$ progenitors from Woosley et al.~(2007) and discuss the difference of our IDSA results with those existing in the literature. Additionally, we perform Newtonian 1D and 2D simulations from prebounce core collapse to several hundred milliseconds postbounce with 11, 15, 21, and 27 $M_\odot$ progenitors from Woosley et al.~(2002) with the HS(DD2) equation of state. General relativistic effects are neglected. We obtain robust explosions with diagnostic energies $E_{\rm dig} \gtrsim 0.1- 0.5$~B for all considered 2D models within approximately $100-300$ milliseconds after bounce and find that explosions are mostly dominated by the neutrino-driven convection, although standing accretion shock instabilities are observed as well. We also find that the level of electron deleptonization during collapse dramatically affect the postbounce evolution, e.g.~the ignorance of neutrino-electron scattering during collapse will lead to a stronger explosion.Comment: 23 pages. Accepted for publication in Ap

Multi-dimensional Core-Collapse Supernova Simulations with Neutrino Transport

We present multi-dimensional core-collapse supernova simulations using the Isotropic Diffusion Source Approximation (IDSA) for the neutrino transport and a modified potential for general relativity in two different supernova codes: FLASH and ELEPHANT. Due to the complexity of the core-collapse supernova explosion mechanism, simulations require not only high-performance computers and the exploitation of GPUs, but also sophisticated approximations to capture the essential microphysics. We demonstrate that the IDSA is an elegant and efficient neutrino radiation transfer scheme, which is portable to multiple hydrodynamics codes and fast enough to investigate long-term evolutions in two and three dimensions. Simulations with a 40 solar mass progenitor are presented in both FLASH (1D and 2D) and ELEPHANT (3D) as an extreme test condition. It is found that the black hole formation time is delayed in multiple dimensions and we argue that the strong standing accretion shock instability before black hole formation will lead to strong gravitational waves.Comment: 3 pages, proceedings for Nuclei in the Cosmos XIV, Niigata, Japan (2016

Silicon Burning II: Quasi-Equilibrium and Explosive Burning

Having examined the application of quasi-equilibrium to hydrostatic silicon burning in Paper I of this series, Hix & Thielemann (1996), we now turn our attention to explosive silicon burning. Previous authors have shown that for material which is heated to high temperature by a passing shock and then cooled by adiabatic expansion, the results can be divided into three broad categories; \emph{incomplete burning}, \emph{normal freezeout} and \emph{$\alpha$-rich freezeout}, with the outcome depending on the temperature, density and cooling timescale. In all three cases, we find that the important abundances obey quasi-equilibrium for temperatures greater than approximately 3 GK, with relatively little nucleosynthesis occurring following the breakdown of quasi-equilibrium. We will show that quasi-equilibrium provides better abundance estimates than global nuclear statistical equilibrium, even for normal freezeout and particularly for $\alpha$-rich freezeout. We will also examine the accuracy with which the final nuclear abundances can be estimated from quasi-equilibrium.Comment: 27 pages, including 15 inline figures. LaTeX 2e with aaspp4 and graphicx packages. Accepted to Ap

James W. Truran (1940–2022)

Truran made far-reaching contributions to the theory of nuclear reactions, stellar nucleosynthesis, big bang nucleosynthesis, stellar abundances, solar system formation, galaxy formation and galactic chemical evolution