Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAdvancing Nucleosynthesis in Self-consistent, Multidimensional Models of Core-Collapse Supernovae
129
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We investigate core-collapse supernova (CCSN) nucleosynthesis in polar axisymmetric simulations using the multidimensional radiation hydrodynamics code CHIMERA. Computational costs have traditionally constrained the evolution of the nuclear composition in CCSN models to, at best, a 14-species $alpha$-network. Such a simplified network limits the ability to accurately evolve detailed composition, neutronization and the nuclear energy generation rate. Lagrangian tracer particles are commonly used to extend the nuclear network evolution by incorporating more realistic networks in post-processing nucleosynthesis calculations. Limitations such as poor spatial resolution of the tracer particles, estimation of the expansion timescales, and determination of the mass-cut at the end of the simulation impose uncertainties inherent to this approach. We present a detailed analysis of the impact of these uncertainties on post-processing nucleosynthesis calculations and implications for future models.
rate research
Read More
Multidimensional hydrodynamic simulations of shell convection in massive stars suggest the development of aspherical perturbations that may be amplified during iron core-collapse. These perturbations have a crucial and qualitative impact on the delayed neutrino-driven core-collapse supernova explosion mechanism by increasing the total stress behind the stalled shock. In this paper, we investigate the properties of a 15 msun model evolved in 1-,2-, and 3-dimensions (3D) for the final $sim$424 seconds before gravitational instability and iron core-collapse using MESA and the FLASH simulation framework. We find that just before collapse, our initially perturbed fully 3D model reaches angle-averaged convective velocity magnitudes of $approx$ 240-260 km s$^{-1}$ in the Si- and O-shell regions with a Mach number $approx$ 0.06. We find the bulk of the power in the O-shell resides at large scales, characterized by spherical harmonic orders ($ell$) of 2-4, while the Si-shell shows broad spectra on smaller scales of $ellapprox30-40$. Both convective regions show an increase in power at $ell=5$ near collapse. We show that the 1D texttt{MESA} model agrees with the convective velocity profile and speeds of the Si-shell when compared to our highest resolution 3D model. However, in the O-shell region, we find that texttt{MESA} predicts speeds approximately emph{four} times slower than all of our 3D models suggest. All eight of the multi-dimensional stellar models considered in this work are publicly available.
We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11-28 solar masses, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si-O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si-O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection time scales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.
Core-collapse supernovae (CCSNe) are the extremely energetic deaths of massive stars. They play a vital role in the synthesis and dissemination of many heavy elements in the universe. In the past, CCSN nucleosynthesis calculations have relied on artificial explosion methods that do not adequately capture the physics of the innermost layers of the star. The PUSH method, calibrated against SN1987A, utilizes the energy of heavy-flavor neutrinos emitted by the proto-neutron star (PNS) to trigger parametrized explosions. This makes it possible to follow the consistent evolution of the PNS and to ensure a more accurate treatment of the electron fraction of the ejecta. Here, we present the Iron group nucleosynthesis results for core-collapse supernovae, exploded with PUSH, for two different progenitor series. Comparisons of the calculated yields to observational metal-poor star data are also presented. Nucleosynthesis yields will be calculated for all elements and over a wide range of progenitor masses. These yields can be immensely useful for models of galactic chemical evolution.
We review some of the uncertainties in calculating nucleosynthetic yields, focusing on the explosion mechanism. Current yield calculations tend to either use a piston, energy injection, or enhancement of neutrino opacities to drive an explosion. We show that the energy injection, or more accurately, an entropy injection mechanism is best-suited to mimic our current understanding of the convection-enhanced supernova engine. The enhanced neutrino-opacity technique is in qualitative disagreement with simulations of core-collapse supernovae and will likely produce errors in the yields. But piston-driven explosions are the most discrepant. Piston-driven explosion severely underestimate the amount of fallback, leading to order-of-magnitude errors in the yields of heavy elements. To obtain yields accurate to the factor of a few level, we must use entropy or energy injection and this has become the NuGrid collaboration approach.
In a previously presented proof-of-principle study, we established a parametrized spherically symmetric explosion method (PUSH) that can reproduce many features of core-collapse supernovae for a wide range of pre-explosion models. The method is based on the neutrino-driven mechanism and follows collapse, bounce and explosion. There are two crucial aspects of our model for nucleosynthesis predictions. First, the mass cut and explosion energy emerge simultaneously from the simulation (determining, for each stellar model, the amount of Fe-group ejecta). Second, the interactions between neutrinos and matter are included consistently (setting the electron fraction of the innermost ejecta). In the present paper, we use the successful explosion models from Ebinger et al. (2018) which include two sets of pre-explosion models at solar metallicity, with combined masses between 10.8 and 120 M$_{odot}$. We perform systematic nucleosynthesis studies and predict detailed isotopic yields. The resulting $^{56}$Ni ejecta are in overall agreement with observationally derived values from normal core-collapse supernovae. The Fe-group yields are also in agreement with derived abundances for metal-poor star HD84937. We also present a comparison of our results with observational trends in alpha element to iron ratios.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
