In this fourth paper of the series, we use the parametrized, spherically symmetric explosion method PUSH to perform a systematic study of two sets of non-rotating stellar progenitor models. Our study includes pre-explosion models with metallicities Z=0 and Z=Z$_{\odot}\times 10^{-4}$ and covers a progenitor mass range from 11 up to 75 M$_\odot$. We present and discuss the explosion properties of all models and predict remnant (neutron star or black hole) mass distributions within this approach. We also perform systematic nucleosynthesis studies and predict detailed isotopic yields as function of the progenitor mass and metallicity. We present a comparison of our nucleosynthesis results with observationally derived $^{56}$Ni ejecta from normal core-collapse supernovae and with iron-group abundances for metal-poor star HD~84937. Overall, our results for explosion energies, remnant mass distribution, $^{56}$Ni mass, and iron group yields are consistent with observations of normal CCSNe. We find that stellar progenitors at low and zero metallicity are more prone to BH formation than those at solar metallicity, which allows for the formation of BHs in the mass range observed by LIGO/VIRGO.

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. The present paper goes beyond a specific application that is able to reproduce observational properties of SN1987A and performs a systematic study of the explosion properties for an extensive set of non-rotating, solar metallicity stellar progenitor models in the mass range from 10.8 to 120 M$_\odot$.This includes the transition from neutron stars to black holes as the final result of the collapse of massive stars, and the relation of the latter to supernovae and faint/failed supernovae. The present paper provides the basis for extended nucleosynthesis predictions in a forthcoming paper to be employed in galactic evolution models.

Dependable predictions of the X-ray burst ashes and light curves require a stringent constraint on the composition of the accreted material as an input parameter. Lower metallicity models are generally based on a metal deficient donor and all metals are summed up in CNO abundances or solar metal distribution is assumed. In this work, we study the alteration of accreted composition due to spallation in the atmosphere of accreting neutron stars considering a cascading destruction process. We find that the inclusion of the cascading process brings the replenishment of CNO elements and overall survival probability is higher compared to isolated destruction of CNO elements. Spallation model provides the distribution of metals as a function of mass accretion rate. Multi-zone X-ray burst models calculated with reduced metallicities have enhanced abundances for high-mass nuclei in X-ray burst ashes. The increased metallicity due to the replenishment of CNO elements changes the composition of burst ashes compared to lower metallicity conditions. This will modify the thermal and compositional structure of accreted neutron star crusts.

Modeling core-collapse supernovae (CCSNe) with neutrino transport in three dimensions (3D) requires tremendous computing resources and some level of approximation. We present a first comparison study of CCSNe in 3D with different physics approximations and hydrodynamics codes. We aim to assess the impact of the hydrodynamics code, approximations for the neutrino and gravity treatments, and rotation on the simulation of CCSNe in 3D. We use four different hydrodynamics codes in this work (ELEPHANT, FLASH, fGR1, and SPHYNX) in combination with two different neutrino treatments, the isotropic diffusion source approximation (IDSA) and two-moment M1, and three different gravity treatments: Newtonian, 1D General Relativity (GR) correction, and full GR). Additional parameters discussed in this study are the inclusion of neutrino-electron scattering via a parametrized deleptonization (PD) and the influence of rotation. The four codes compared in this work include Eulerian and fully Lagrangian (smoothed particle hydrodynamics) codes for the first time. They show agreement in the overall evolution of the collapse phase and early post-bounce within the range of 10% (20% in some cases). The comparison of the different neutrino treatments highlights the need to further investigate the antineutrino luminosities in IDSA, which tend to be relatively high. We also demonstrate the requirement for a more detailed heavy-lepton neutrino leakage. When comparing with a full GR code, including an M1 transport method, we confirm the influence of neutrino-electron scattering during the collapse phase, which is adequately captured by the PD scheme. Also, the effective GR potential reproduces the overall dynamic evolution correctly in all Newtonian codes. Additionally, we verify that rotation aids the shock expansion and estimate the overall angular momentum losses for each code in rotating scenarios.

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.

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 report on a method, PUSH, for triggering core-collapse supernova (CCSN) explosions of massive stars in spherical symmetry. This method provides a framework to study many important aspects of core collapse supernovae: the effects of the shock passage through the star, explosive supernova nucleosynthesis and the progenitor-remnant connection. Here we give an overview of the method, compare the results to multi-dimensional simulations and investigate the effects of the progenitor and the equation of state on black hole formation.

We report on a method, PUSH, for triggering core-collapse supernova explosions of massive stars in spherical symmetry. We explore basic explosion properties and calibrate PUSH such that the observables of SN1987A are reproduced. Our simulations are based on the general relativistic hydrodynamics code AGILE combined with the detailed neutrino transport scheme IDSA for electron neutrinos and ALS for the muon and tau neutrinos. To trigger explosions in the otherwise non-exploding simulations, we rely on the neutrino-driven mechanism. The PUSH method locally increases the energy deposition in the gain region through energy deposition by the heavy neutrino flavors. Our setup allows us to model the explosion for several seconds after core bounce. We explore the progenitor range 18-21M$_{\odot}$. Our studies reveal a distinction between high compactness (HC) and low compactness (LC) progenitor models, where LC models tend to explore earlier, with a lower explosion energy, and with a lower remnant mass. HC models are needed to obtain explosion energies around 1 Bethe, as observed for SN1987A. However, all the models with sufficiently high explosion energy overproduce $^{56}$Ni. We conclude that fallback is needed to reproduce the observed nucleosynthesis yields. The nucleosynthesis yields of $^{57-58}$Ni depend sensitively on the electron fraction and on the location of the mass cut with respect to the initial shell structure of the progenitor star. We identify a progenitor and a suitable set of PUSH parameters that fit the explosion properties of SN1987A when assuming 0.1M$_{\odot}$ of fallback. We predict a neutron star with a gravitational mass of 1.50M$_{\odot}$. We find correlations between explosion properties and the compactness of the progenitor model in the explored progenitors. However, a more complete analysis will require the exploration of a larger set of progenitors with PUSH.

The smoothed-particle hydrodynamics (SPH) technique is a numerical method for solving gas-dynamical problems. It has been applied to simulate the evolution of a wide variety of astrophysical systems. The method has a second-order accuracy, with a resolution that is usually much higher in the compressed regions than in the diluted zones of the fluid. In this work, we propose and check a scheme to balance and equalize the resolution of SPH between high- and low-density regions. This method relies on the versatility of a family of interpolators called Sinc kernels, which allows increasing the interpolation quality by varying only a single parameter (the exponent of the Sinc function). The scheme is checked and validated through a number of numerical tests, from standard one-dimensional Riemann problems in shock tubes, to multidimensional simulations of explosions, hydrodynamic instabilities and the collapse of a sun-like polytrope. The analysis of the hydrodynamical simulations suggests that the scheme devised to equalizing accuracy improves the treatment of the post-shock regions and, in general, of the rarefacted zones of fluids while causing no harm to the growth of hydrodynamic instabilities. The method is robust and easy to implement with a low computational overload. It conserves mass, energy, and momentum and reduces to the standard SPH scheme in regions of the fluid that have smooth density gradients.