Open Access
Issue
A&A
Volume 686, June 2024
Article Number L8
Number of page(s) 10
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202449994
Published online 31 May 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Massive stars (with initial masses > 8−10 M) end their lives with the formation and the collapse of a Fe core. The collapse halts and bounces back as soon as the central infalling material reaches the typical nuclear density. The bounce of the innermost regions of the star generates a shock wave that, in some cases, is able to emerge from the Fe core and propagate into the envelope of the star (Burrows & Vartanyan 2021). The shock wave induces local compression and heating, and this causes nucleosynthesis to occur within a few seconds rather than the days or even years it takes during the pre-supernova evolution. This explosive nucleosynthesis alters the chemical composition produced during the hydrostatic evolution (see, e.g., Rauscher et al. 2002; Pignatari et al. 2016a; Limongi & Chieffi 2018). The most abundant nuclear species in compositions of this so-called core-collapse supernova (CCSN) ejecta are the Fe-peak and α elements, mostly synthesized through (quasi)statistical equilibrium processes in the explosion, together with the stellar burning products from the hydrostatic evolution of the star (Woosley & Weaver 1995).

Core-collapse supernovae are also one of the production sites of the so-called p nuclei via the γ process (Rayet et al. 1995; Arnould & Goriely 2003; Pignatari et al. 2016b; Roberti et al. 2023). The ONe-rich layers in massive stars during a CCSN provide the necessary conditions (i.e., high temperatures and abundances of trans-iron seeds) to synthesize the p nuclei via a chain of photo-disintegrations (Woosley & Howard 1978). Several discrepancies still arise when comparing theoretical model predictions to each other and with solar abundances (Roberti et al. 2023). In particular, (i) the production is too low to match the solar abundances; (ii) the isotopic ratios of theoretical yields produced in the same temperature range cannot reproduce the solar ratios; and (iii) 92, 94Mo and 96, 98Ru show an even more severe underproduction compared to the other p nuclei (Rauscher et al. 2013; Pignatari et al. 2016b). While additional nucleosynthetic processes and astrophysical phenomena have been invoked to explain the solar abundances of p nuclei (see Roberti et al. 2023 for a recent overview of the production sites of p nuclei), this remains an open problem. For example, C−O shell mergers in the late stage of massive star evolution can play a crucial role in boosting the final yields of p nuclei, but this does not explain the underproduction of the Mo and Ru isotopes (Rauscher et al. 2002; Ritter et al. 2018a).

Choplin et al. (2022) show that fast rotation in a 25 M star at subsolar metallicity can enhance the efficiency of the γ-process nucleosynthesis, due to a higher production of trans-iron seeds via the slow neutron capture nucleosynthesis (the s process). The possibility of including such results in models of galactic chemical evolution needs to be explored as it may modify our current understanding, that the contribution from previous generations of massive stars to the Solar System abundances of p nuclei is negligible (Travaglio et al. 2018). Choplin et al. (2022) also briefly discuss how the different explosion energies affect the production of p nuclei for a wide range of energies (0.3−100 foe). They conclude that, at least in the specific case of their 25 M progenitor, the explosion energy has only a marginal impact on the total yields of p nuclei because the only effect of increasing the energy is the shift of the γ-process zone outward in mass. Nevertheless, a detailed study of the dependence of γ-process nucleosynthesis, and of the production of p nuclei in general, on different explosion energies and explosive prescriptions is still missing. With this work we aim to start an exploration in this direction.

In the first paper of this series (Roberti et al. 2023, hereafter Paper I), we presented the first step of our study of γ-process nucleosynthesis in CCSNe. We analyzed and compared the yields from five different sets of massive star models from the literature. In this second work, we discuss the effect of different explosive prescriptions on 15, 20, and 25 M star progenitors. In Sect. 2 we describe the different explosive prescriptions; Sect. 3 summarizes the conditions that affect the γ-process nucleosynthesis. In Sect. 4 we discuss how the different explosion prescriptions influence the average overproduction factors for p nuclei, and in Sect. 5 we discuss and summarize our results.

2. Datasets

We explored the effect of different supernova prescriptions on the γ process and on the p-nucleus yields in the two set of models described below, with initial masses of 15, 20, and 25 M and solar metallicity. We considered a broad range of explosion possibilities using different parameters that describe both the explosion energy and the amount of material that is ejected (versus the amount that is incorporated onto the compact object). Multiple parameters are currently required to capture the uncertainties in the explosive engine. The explosion energy depends on the growth time of the convection (Fryer et al. 2021, 2022; Boccioli et al. 2023) as well as on the existence of energy sources, such as continued accretion after the shock is launched (especially if there is sufficient angular momentum to form an accretion disk) or the development of a magnetar (Shankar et al. 2021; Soker 2022). The mass cut also depends on the convective growth time (Fryer et al. 2022). Fallback (which also depends on the explosion energy) can also impact the final remnant mass and can be particularly important in asymmetric explosions (Fryer 2009; Fryer et al. 2012). Our two supernova prescriptions help cover the broad uncertainties in the ejecta from supernovae.

2.1. Sedov blast-wave (SBW) explosion

We simulated a new set of CCSN explosions using the same approach and the same three massive star models calculated with the MESA code (Paxton et al. 2015, and references therein) by Ritter et al. (2018b) at solar metallicity (Z = 0.02, based on Grevesse & Noels 1993). We varied the initial shock velocity within the range 2 × 108 cm s−1 to 5 × 109 cm s−1, which correspond to a factor of 10 below and 2.5 above the Ritter et al. (2018b) standard supernova explosion setup. We briefly recall here the main features of the explosive prescription. The supernova explosion was modeled using a semi-analytical approach, with the prescription for the remnant mass from Fryer et al. (2012) and the Sedov (1946) blast wave (SBW) solution for the determination of the peak velocity of the shock throughout the stellar structure. The temperature and the density in each layer of the star after the shock were then obtained by imposing a strong shock limit (Chevalier 1989) and the radiation-dominated shock (Eqs. (3)–(5) of Pignatari et al. 2016b). After the passage of the shock, the temperature and the density decrease with time as a consequence of the cooling of the shocked material. This temporal evolution was described using a variant of the adiabatic exponential decay (Hoyle et al. 1964; Fowler & Hoyle 1964). As in Ritter et al. (2018b), we set an artificial cap to the maximum velocity of 5 × 109 cm s−1 to simulate the deceleration due to viscous forces. We note that in this set we kept the remnant mass fixed using the prescription from Fryer et al. (2012) and decided to study the effect on the γ-process nucleosynthesis of the explosion energy only.

Table 1.

Selected features of the CCSN models calculated using the SBW approach.

Both the hydrostatic and explosive nucleosynthesis were calculated using the same version of the JINA Reaclib database used by Ritter et al. (2018b) (V1.1; Cyburt 2011) and with an updated version of the Multi-zone Post-Processing Network – Parallel (MPPNP) code (Pignatari & Herwig 2012; Pignatari et al. 2016b; Ritter et al. 2018b, Paper I). Table 1 presents a summary of the main features of these models.

2.2. Lawson set (LAW)

We further used the CCSN yields from the models of Jones et al. (2019), Andrews et al. (2020), and Lawson et al. (2022, hereafter LAW) to study how different supernova parameterizations affect the production of p nuclei. The progenitors were computed with the 1D Kepler code (Heger & Woosley 2010) using the solar metallicity from Grevesse & Noels (1993), Z = 0.02. The explosions were calculated using a 1D hydrodynamic code that mimics the 3D convection-enhanced supernova engine (Fryer et al. 2018). Fryer et al. (2018) performed a parametric study of the explosion energy, compact remnant production, and basic nucleosynthetic yields of CCSNe by varying the power (Einj), duration (tinj), and extent of the energy injection region (Minj) in the stellar progenitors. On the basis of this structure and evolution, Jones et al. (2019), Andrews et al. (2020), and Lawson et al. (2022) calculated the hydrostatic nucleosynthesis using the MPPNP code and the explosive nucleosynthesis using the Tracer particle Post-Processing Network – Parallel (TPPNP) code (Jones et al. 2019). Both the hydrostatic and explosive nucleosynthesis were calculated using version V2.2 of the JINA Reaclib database. These calculations consist in a set of yields from 80 different explosive parametrizations, including 23 explosions of the 15 M model, 31 of the 20 M model, and 26 of the 25 M model1. The Einj ranged between 3 and 200 foe, Minj between 0.02 and 0.2 M, and tinj between 0.1 and 1 s. They produce explosion energies and remnant masses in the interval between 0.34 and 18.4 foe and 1.5 and 5.6 M, respectively.

3. The γ process and the CCSN explosion

The temperature required to produce p nuclei via the γ process ranges between 2 and 3.5 GK (see, e.g., Pignatari et al. 2016b), which in CCSNe roughly corresponds to the typical temperature of explosive Ne burning (Limongi & Chieffi 2008). Based on this temperature range, we defined the “hot” γ-process component as the light p nuclei that are produced at a temperature near ∼3.5 GK (for example, 74Se and 102Pd) and the “cold” γ-process component as the intermediate and heavy p nuclei that are produced at a temperature near ∼2 GK (for example, 130Ba and 196Hg).

thumbnail Fig. 1.

Peak shock temperature (solid red lines) and velocity (solid green lines) for different initial explosive conditions in the 20 M progenitor from Ritter et al. (2018b). The purple line represents the logarithm of the quantity ρr3 in the pre-supernova model in arbitrary units. The gray region represents the temperature range for the activation of the γ process, i.e., between 2 GK and 3.5 GK. Note that the saturation of the shock velocity at ∼4.2 M at the value of 5 × 109 cm s−1 is due to the artificial cap imposed for viscous forces (see Sect. 2.1).

thumbnail Fig. 2.

Effect of different explosion energies and prescriptions on three representative p-nuclei. Left panels: 92Mo (upper panel), 130Ba (central panel), and 196Hg (lower panel) abundances in mass fractions as a function of the internal mass coordinate (solid lines, left axis) and explosion peak temperature in T9 (dashed lines, right y-axis) in the 20 M SBW models. Different colors denote different shock velocities. Right panels: same as the left panels but for a sample of the 15 M LAW models, without the explosion peak temperature plotted. The horizontal dashed-dotted blue line in each panel represents the solar abundance of the representative isotope in mass fractions taken from Asplund et al. (2009).

The main effect of adopting different explosion prescriptions is to shift the mass coordinate in the structure of the supernova progenitor where the shock wave reaches a certain velocity and, therefore, a certain temperature. Consequently, a basic assumption that has often been made is that the production of p nuclei is almost independent of the properties of the explosion, since a different propagation of the shock would, in principle, only shift the location mass of the explosive Ne burning (see Choplin et al. 2022, and our Table 1). However, according to the Sedov solution2, a local decrease in the quantity ρr3 would facilitate the local propagation of the shock wave and its acceleration (Sedov 1946; Woosley et al. 2002). The acceleration results in a temperature increase, which in some cases (especially for high-energy explosions) can lead to the temperature threshold being exceed again for the γ process and can increase the extent in mass of the zone where the p nuclei are produced. Figure 1 shows this occurrence in the case of a 20 M model from the SBW set. In the case of the most energetic explosions, the peak temperature spends more time in the γ-process area (represented by the gray band in the figure) due to the acceleration of the shock wave where ρr3 decreases. The left panels of Fig. 2 show the effect of the late acceleration of the shock wave (in the same model as Fig. 1) on the final abundances of three representative p nuclei (92Mo, 130Ba, and 196Hg). The abundance peaks corresponding to the p nuclei are progressively more extended in mass with increasing initial shock velocity, at least in the models with vshock ≥ 2 × 109 cm s−1. However, in the most energetic explosions, the abundance of all the p nuclei drops abruptly at ∼4.5 M. This mass coordinate corresponds to the interface between the CO and He cores, where the temperature and the density of the shock sharply decrease. If the temperature is still high enough for an efficient γ-process nucleosynthesis when the shock wave reaches the CO–He core interface, the production of the cold γ-process components are generally suppressed because the peak temperature decrease is too sharp to extend farther into the γ-process region. We observe a qualitatively similar, albeit less intense, behavior in the LAW models (right panels of Fig. 2). In this case, the interplay among the multiple explosion parameters often leads to a non-monotonic increase in the γ-process yields with the final explosion energy.

The choice of mass cut (i.e., the mass coordinates that separate the ejecta from the remnant mass) plays a relevant role in the ejection of p nuclei. A more external mass cut will (partly) lock the ashes of the explosive Ne burning within the remnant; conversely, a deeper mass cut will allow the supernova to fully eject the γ-process products and, in some cases, to eject some material rich in the lightest p nuclei (74Se, 78Kr, 84Sr, 92 − 94Mo, and, possibly, 96 − 98Ru) produced by the α-rich freeze-out (see, e.g., Woosley et al. 2002; Rauscher et al. 2013; Pignatari et al. 2016b; Lugaro et al. 2016). The upper right panel of Fig. 2 shows the large abundances of the isotope 92Mo synthesized in the α-rich freeze-out in some of the 15 M models from the LAW set. To the contrary, no significant α-rich freeze-out component is ejected in any of the SBW models. We recall that in the SBW models the mass cut is independent of the explosion energy because it is predetermined by the adopted semi-analytical prescription, while for LAW models it is the result of the solution of the hydrodynamic equations.

In conclusion, different choices regarding the explosion prescription can significantly alter the amount of p-nucleus-rich material ejected by the CCSN. We analyze this effect in more detail in the next section.

4. Discussion: Effect of different explosion prescriptions on the production of p nuclei

We defined the overproduction factor, Fi, of a certain isotope, i, as the ratio between its mass fraction, Xi (which is equal to the total integrated yield divided by the total mass ejected) and the solar mass fraction, Xi, ⊙, from Asplund et al. (2009). In Paper I we discussed the use of the average overproduction factor, F0, defined as F 0 = ( i = 1 35 F i ) / 35 $ F_0 = (\sum_{i=1}^{35} F_i)/35 $, as a proxy for the overall p-nucleus production (see, e.g., Arnould & Goriely 2003; Pignatari et al. 2013). We also introduced alternative definitions of this parameter to disentangle the γ process from the other processes that can contribute to the production of p nuclei. In particular, we defined Fγ as the average overproduction factor of the three most produced γ-only nuclei. Since Fγ reflects the dominant γ-process production, we can consider the isotopes contributing to Fγ as “thermometers” for the γ process and use them to characterize the differences between the explosive prescriptions. To approximately estimate the contribution from CCSNe to the solar abundances of the p nuclei and to discuss how the p nucleus synthesis changes for different models, we compared F0 and Fγ with the overproduction of 16O (F[O16]), the most abundant metal ejected from massive stars. The F0, Fγ, and F[O16] of each model are given in Tables A.1A.3.

thumbnail Fig. 3.

F0/F[O16] (blue diamonds) and Fγ/F[O16] (green diamonds) in the SBW models as a function of the initial shock velocity (left panels) and in the LAW models as a function of the explosion energy (right panels). The dotted lines represent Fi = F[O16], and the shaded gray area represents a factor of 2 above and below the dotted line.

Figure 3 shows the F0/F[O16] and Fγ/F[O16] ratios in the SBW and LAW sets as a function of the initial mass and of their trend with the explosion energy or, equivalently, the initial shock velocity. In the 20 and 25 M SBW models, we find that that the increase in the explosion energy produces, on average, an increase within one order of magnitude of the p-nucleus yields relative to 16O, in contrast with Choplin et al. (2022). Almost no p nuclei are ejected for low explosion velocities (where the mass cut is larger than the outer mass coordinate for explosive Ne burning). At higher velocities, there is a plateau corresponding to F0 ∼ F[O16], caused by the fact that the external zone of γ-process production coincides with the interface between the CO and He cores, which suppresses the cold γ-process component (see the bottom- and central-left panels of Fig. 2, which show the representative cases of 130Ba and 196Hg). The production of the hot component is instead boosted by the acceleration of the shock wave in the outer region of the C shell, because the temperature is kept high enough to have a larger γ-process zone (the top-left panel of Fig. 2 shows the representative case of 92Mo). The Fγ increases more than F0 because it encompasses fewer isotopes that contribute the most to F0. The 15 M SBW models instead experience a C−O shell merger before the explosion. As discussed in Paper I, the yields of all the p nuclei heavier than 108Cd are almost exclusively dominated by the production obtained in the C−O shell mergers, with a negligible contribution from the explosive nucleosynthesis. Here we confirm this results: we find that, across a wide range of explosion energies, the average overproduction is constant and is equal to or higher than that of 16O (with F0/F[O16] ∼ 1, Fγ/F[O16] ∼ 2 − 3). Moreover, the most produced p nuclei do not change with increasing explosion energy, because the dominant contribution is always produced at the typical temperature of the O-burning shell in the pre-supernova stage. In the whole explosion velocity range, 130, 132Ba and 144Sm are the nuclei that contribute to Fγ. We note, however, that neither the more energetic explosions nor the C−O shell mergers can fully fill the production gap between the Mo−Ru isotopes and the other p nuclei in the SBW models.

In LAW models the explosion energy, Eexp, is the hydrodynamic result of a complex interplay between three initial parameters: the amount of injected energy (Einj) and the duration (tinj) and mass location (Minj) of the energy deposition (Sect. 2.2). Only models with Eexp >  2 foe have an efficient production of p nuclei (i.e., F0/F[O16] and/or Fγ/F[O16] ≥ 1). We note, however, that in the 15 M models the p nucleus production is often dominated by the ejection of some α-rich freeze-out material rather than the γ-process nucleosynthesis. In nine of the 15 M models, in fact, we find values of F0/F[O16] as high as 200, but with Fγ/F[O16] ≤ 1, due to the extremely large yields of some light p nuclei (from 74Se to 94Mo) that are produced very close to the mass cut, just below the abundance peak of 56Ni, as shown in the top-right panel of Fig. 2.

Finally, we find that in LAW models, Fγ receives the most contribution from intermediate and heavy p nuclei, while in SBW models it receives the most contribution from the light and intermediate p nuclei (see Tables A.1A.3). It follows that the SBW models experience hotter γ-process nucleosynthesis compared to LAW models. This is in qualitative agreement with what was already found by Woosley & Weaver (1995) and Limongi & Chieffi (2003), who observed that a radiation-dominated technique tends to result in higher peak temperatures than those estimated via a hydrodynamic calculation.

5. Summary and conclusions

We studied the effect of a wide range of explosion energies and parameters on the production of p nuclei in CCSNe. We compared the results obtained from a semi-analytical model (SBW; Pignatari et al. 2016a; Ritter et al. 2018b) with those from a set of hydrodynamic simulations (LAW; Fryer et al. 2018; Lawson et al. 2022), using two different definitions of the average p-nucleus overproduction factors (F0 and Fγ). As discussed in Paper I (see also Farmer et al. 2016; Boccioli & Roberti 2024), small differences in the input parameters used by different authors to calculate massive star models with the same initial masses and metallicities (and potentially rotation velocities) can lead to very different structures at the time of the collapse. The propagation of the shock wave and the explosive nucleosynthesis (and therefore the γ process) are then influenced by these differences. For this reason, finding that the overall production of p nuclei increases with the increase in the explosion power, regardless of the set and the explosion method used, is a remarkable result.

In the case of SBW models that do not experience a C−O shell merger (the 20 and 25 M models), the increase in the explosion energy produces an increase within one order of magnitude of the p-nucleus yields. Conversely, the p-nucleus yields from the model that experiences a C−O shell merger during the pre-supernova stage (15 M) are not significantly affected by the explosion because the pre-supernova production survives and dominates over the explosive nucleosynthesis for the wide range of explosion velocities we have explored. In any case, neither the increase in the explosion energy nor the occurrence of C−O shell mergers resolve the problem of the γ-process underproduction of Mo and Ru isotopes.

In the LAW models there is no monotonic trend with explosion energy, but there is an overall increase in the p-nucleus production in the case of more energetic supernovae. The ejection of some α-rich freeze-out material, in the case of the explosions above 2 foe in the 15 M models, increases the Mo−Ru overproduction relative to 16O by up to two orders of magnitude. Recently, Sasaki et al. (2022) suggested that high-energy hypernovae (Eexp ≥ 10 foe) could also provide the right conditions for an efficient νp process (Fröhlich et al. 2006), which may resolve the problem with Mo−Ru isotopes. In any case, high-energy supernovae could represent the key to resolving the long-pending problem about the underproduction of 92, 94Mo and 96Ru, although not via the γ process alone. We further note that the models included in our analysis did not consider neutrino wind nucleosynthesis, and therefore we are not able to estimate the impact at different explosion energies of the νp process nor the neutrino spallation contribution, which is crucial for the production of the p-nucleus 138La and 180Ta (see Cumming & Alburger 1985; Goriely & Siess 2001, and the detailed description in Paper I).

Finally, we stress that with this work we explored only a portion of the wide parameter space that determines the production of p nuclei. For instance, more recent and sophisticated multidimensional simulations of the final stages and explosion of massive stars show some features that are impossible to reproduce in 1D codes. Rizzuti et al. (2022) show that a significant entrainment is found in convective Ne-shell burning, confirming that the burning stages that lead to the collapse and to the explosion are still poorly constrained by 1D models. Moreover, most 3D CCSN simulations show that the dynamics of the innermost layers of the ejecta is strongly aspherical, with the formation of structures such as bubbles, plumes, and bullets that may even exceed the speed of the outer material and escape the fallback (Wongwathanarat et al. 2015; Burrows & Vartanyan 2021). Asymmetries in the explosion may eventually affect the spatial distribution of the free electron fraction, Ye, near the Fe core and, in turn, the production of the lightest p nuclei (up to Pd) via α-rich freeze-out or the νp process (Hoffman et al. 1996; Arcones & Thielemann 2013). A further exploration of the production of nuclei in multidimensional simulations will be needed in the future to confirm these predictions.


1

The full datasets of these models can be found at https://ccsweb.lanl.gov/astro/nucleosynthesis/nucleosynthesis_astro.html

2

From the Sedov solution: Eexp ∝ ρ r5 t−2 → vshock ∝ (ρ r3)−1.

Acknowledgments

We thank the support from the NKFI via K-project 138031 and the Lendület Program LP2023-10 of the Hungarian Academy of Sciences. L.R. and M.P. acknowledge the support to NuGrid from JINA-CEE (NSF Grant PHY-1430152) and STFC (through the University of Hull’s Consolidated Grant ST/R000840/1), and ongoing access to viper, the University of Hull High Performance Computing Facility. L.R. acknowledges the support from the ChETEC-INFRA – Transnational Access Projects 22102724-ST and 23103142-ST. This work was supported by the European Union’s Horizon 2020 research and innovation programme (ChETEC-INFRA – Project no. 101008324), and the IReNA network supported by US NSF AccelNet (Grant No. OISE-1927130). M.L. was also supported by the NKFIH excellence grant TKP2021-NKTA-64.

References

  1. Andrews, S., Fryer, C., Even, W., Jones, S., & Pignatari, M. 2020, ApJ, 890, 35 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arcones, A., & Thielemann, F. K. 2013, J. Phys. G Nucl. Phys., 40, 013201 [NASA ADS] [CrossRef] [Google Scholar]
  3. Arnould, M., & Goriely, S. 2003, Phys. Rep., 384, 1 [Google Scholar]
  4. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  5. Boccioli, L., & Roberti, L. 2024, Universe, 10, 148 [NASA ADS] [CrossRef] [Google Scholar]
  6. Boccioli, L., Roberti, L., Limongi, M., Mathews, G. J., & Chieffi, A. 2023, ApJ, 949, 17 [CrossRef] [Google Scholar]
  7. Burrows, A., & Vartanyan, D. 2021, Nature, 589, 29 [CrossRef] [PubMed] [Google Scholar]
  8. Chevalier, R. A. 1989, ApJ, 346, 847 [Google Scholar]
  9. Choplin, A., Goriely, S., Hirschi, R., Tominaga, N., & Meynet, G. 2022, A&A, 661, A86 [CrossRef] [EDP Sciences] [Google Scholar]
  10. Cumming, J. B., & Alburger, D. E. 1985, Phys. Rev. C, 31, 1494 [NASA ADS] [CrossRef] [Google Scholar]
  11. Cyburt, R. 2011, APS Div. Nucl. Phys. Meet. Abstr., 56, 7 [Google Scholar]
  12. Farmer, R., Fields, C. E., Petermann, I., et al. 2016, ApJS, 227, 22 [NASA ADS] [CrossRef] [Google Scholar]
  13. Fowler, W. A., & Hoyle, F. 1964, ApJS, 9, 201 [Google Scholar]
  14. Fröhlich, C., Hix, W. R., Martínez-Pinedo, G., et al. 2006, New Astron. Rev., 50, 496 [CrossRef] [Google Scholar]
  15. Fryer, C. L. 2009, ApJ, 699, 409 [NASA ADS] [CrossRef] [Google Scholar]
  16. Fryer, C. L., Belczynski, K., Wiktorowicz, G., et al. 2012, ApJ, 749, 91 [Google Scholar]
  17. Fryer, C. L., Andrews, S., Even, W., Heger, A., & Safi-Harb, S. 2018, ApJ, 856, 63 [CrossRef] [Google Scholar]
  18. Fryer, C. L., Karpov, P., & Livescu, D. 2021, Astron. Rep., 65, 937 [NASA ADS] [CrossRef] [Google Scholar]
  19. Fryer, C. L., Olejak, A., & Belczynski, K. 2022, ApJ, 931, 94 [NASA ADS] [CrossRef] [Google Scholar]
  20. Goriely, S., & Siess, L. 2001, A&A, 378, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Grevesse, N., & Noels, A. 1993, in Origin and Evolution of the Elements, eds. N. Prantzos, E. Vangioni-Flam, & M. Casse, 15 [Google Scholar]
  22. Heger, A., & Woosley, S. E. 2010, ApJ, 724, 341 [Google Scholar]
  23. Hoffman, R. D., Woosley, S. E., Fuller, G. M., & Meyer, B. S. 1996, ApJ, 460, 478 [NASA ADS] [CrossRef] [Google Scholar]
  24. Hoyle, F., Fowler, W. A., Burbidge, G. R., & Burbidge, E. M. 1964, ApJ, 139, 909 [NASA ADS] [CrossRef] [Google Scholar]
  25. Jones, S., Röpke, F. K., Fryer, C., et al. 2019, A&A, 622, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Lawson, T. V., Pignatari, M., Stancliffe, R. J., et al. 2022, MNRAS, 511, 886 [NASA ADS] [CrossRef] [Google Scholar]
  27. Limongi, M., & Chieffi, A. 2003, ApJ, 592, 404 [CrossRef] [Google Scholar]
  28. Limongi, M., & Chieffi, A. 2008, in EAS Publications Series, eds. C. Charbonnel, & J. P. Zahn, 32, 233 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Limongi, M., & Chieffi, A. 2018, ApJS, 237, 13 [NASA ADS] [CrossRef] [Google Scholar]
  30. Lugaro, M., Pignatari, M., Ott, U., et al. 2016, Proc. Natl. Acad. Sci., 113, 907 [NASA ADS] [CrossRef] [Google Scholar]
  31. Paxton, B., Marchant, P., Schwab, J., et al. 2015, ApJs, 220, 15 [NASA ADS] [CrossRef] [Google Scholar]
  32. Pignatari, M., & Herwig, F. 2012, Nucl. Phys. News, 22, 18 [NASA ADS] [CrossRef] [Google Scholar]
  33. Pignatari, M., Hirschi, R., Wiescher, M., et al. 2013, ApJ, 762, 31 [NASA ADS] [CrossRef] [Google Scholar]
  34. Pignatari, M., Göbel, K., Reifarth, R., & Travaglio, C. 2016a, Int. J. Mod. Phys. E, 25, 1630003 [NASA ADS] [CrossRef] [Google Scholar]
  35. Pignatari, M., Herwig, F., Hirschi, R., et al. 2016b, ApJS, 225, 24 [NASA ADS] [CrossRef] [Google Scholar]
  36. Rauscher, T., Heger, A., Hoffman, R. D., & Woosley, S. E. 2002, ApJ, 576, 323 [Google Scholar]
  37. Rauscher, T., Dauphas, N., Dillmann, I., et al. 2013, Rep. Progr. Phys., 76, 066201 [CrossRef] [Google Scholar]
  38. Rayet, M., Arnould, M., Hashimoto, M., Prantzos, N., & Nomoto, K. 1995, A&A, 298, 517 [Google Scholar]
  39. Ritter, C., Andrassy, R., Côté, B., et al. 2018a, MNRAS, 474, L1 [NASA ADS] [CrossRef] [Google Scholar]
  40. Ritter, C., Herwig, F., Jones, S., et al. 2018b, MNRAS, 480, 538 [NASA ADS] [CrossRef] [Google Scholar]
  41. Rizzuti, F., Hirschi, R., Georgy, C., et al. 2022, MNRAS, 515, 4013 [NASA ADS] [CrossRef] [Google Scholar]
  42. Roberti, L., Pignatari, M., Psaltis, A., et al. 2023, A&A, 677, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Sasaki, H., Yamazaki, Y., Kajino, T., et al. 2022, ApJ, 924, 29 [NASA ADS] [CrossRef] [Google Scholar]
  44. Sedov, L. I. 1946, Appl. Math. Mech., 10, 241 [NASA ADS] [Google Scholar]
  45. Shankar, S., Mösta, P., Barnes, J., Duffell, P. C., & Kasen, D. 2021, MNRAS, 508, 5390 [NASA ADS] [CrossRef] [Google Scholar]
  46. Soker, N. 2022, ApJ, 935, 108 [NASA ADS] [CrossRef] [Google Scholar]
  47. Travaglio, C., Rauscher, T., Heger, A., Pignatari, M., & West, C. 2018, ApJ, 854, 18 [NASA ADS] [CrossRef] [Google Scholar]
  48. Wongwathanarat, A., Müller, E., & Janka, H.-T. 2015, A&A, 577, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Woosley, S. E., & Howard, W. M. 1978, ApJS, 36, 285 [NASA ADS] [CrossRef] [Google Scholar]
  50. Woosley, S. E., & Weaver, T. A. 1995, ApJS, 101, 181 [Google Scholar]
  51. Woosley, S. E., Heger, A., & Weaver, T. A. 2002, Rev. Mod. Phys., 74, 1015 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Tables of average overproduction factors

Here we present the different average overproduction factors for all the models analyzed in this work. In particular: 15 M models in Table A.1; 20 M in Table A.2; and 25 M in Table A.3. We repeat below the different definitions for F0, F 0 I , II , III $ _{0}^{\mathrm{I,II,III}} $, and Fγ (see also Table 2 of Roberti et al. 2023). In the case of LAW, we use the same model names as in the online tables from Lawson et al. (2022) (Sect. 2.2). In the case of SBW, the name of each model is composed of the mass in solar masses (M), followed by the absolute metallicity (Z) and the velocity in cm/s (V). The column “Eexp or vshock” contains the final explosion energy in foe (1051 erg) in the case of the LAW models, or the initial shock velocity in 109 cm s−1 in the case of the SBW models.

  • F0: no p nuclei are excluded.

  • F 0 I $ _{0}^{\mathrm{I}} $: the excluded p nuclei are 74Se, 78Kr, 84Sr, 92, 94Mo, 96, 98Ru, 108Cd, 113In, 115Sn, 138La, 152Gd, 164Er, and 180Ta.

  • F 0 II $ _{0}^{\mathrm{II}} $: the excluded p nuclei are 74Se, 78Kr, 84Sr, 92, 94Mo, and 96, 98Ru.

  • F 0 III $ _{0}^{\mathrm{III}} $: the excluded p nuclei are 108Cd, 113In, 115Sn, 138La, 152Gd, 164Er, and 180Ta.

  • Fγ: all p nuclei are excluded except for the three most produced γ-only nuclei, which are listed in the last column.

Table A.1.

Average overproduction factors in the 15 M models.

Table A.2.

Average overproduction factors in the 20 M models.

Table A.3.

Average overproduction factors in the 25 M models.

All Tables

Table 1.

Selected features of the CCSN models calculated using the SBW approach.

Table A.1.

Average overproduction factors in the 15 M models.

Table A.2.

Average overproduction factors in the 20 M models.

Table A.3.

Average overproduction factors in the 25 M models.

All Figures

thumbnail Fig. 1.

Peak shock temperature (solid red lines) and velocity (solid green lines) for different initial explosive conditions in the 20 M progenitor from Ritter et al. (2018b). The purple line represents the logarithm of the quantity ρr3 in the pre-supernova model in arbitrary units. The gray region represents the temperature range for the activation of the γ process, i.e., between 2 GK and 3.5 GK. Note that the saturation of the shock velocity at ∼4.2 M at the value of 5 × 109 cm s−1 is due to the artificial cap imposed for viscous forces (see Sect. 2.1).

In the text
thumbnail Fig. 2.

Effect of different explosion energies and prescriptions on three representative p-nuclei. Left panels: 92Mo (upper panel), 130Ba (central panel), and 196Hg (lower panel) abundances in mass fractions as a function of the internal mass coordinate (solid lines, left axis) and explosion peak temperature in T9 (dashed lines, right y-axis) in the 20 M SBW models. Different colors denote different shock velocities. Right panels: same as the left panels but for a sample of the 15 M LAW models, without the explosion peak temperature plotted. The horizontal dashed-dotted blue line in each panel represents the solar abundance of the representative isotope in mass fractions taken from Asplund et al. (2009).

In the text
thumbnail Fig. 3.

F0/F[O16] (blue diamonds) and Fγ/F[O16] (green diamonds) in the SBW models as a function of the initial shock velocity (left panels) and in the LAW models as a function of the explosion energy (right panels). The dotted lines represent Fi = F[O16], and the shaded gray area represents a factor of 2 above and below the dotted line.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.