Free Access
Issue
A&A
Volume 661, May 2022
Article Number A86
Number of page(s) 11
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202243331
Published online 05 May 2022

© ESO 2022

1. Introduction

Despite tremendous progress during the past decades, the origin of the trans-iron chemical elements is still debated and not yet fully understood (e.g. Arnould & Goriely 2020). The slow (s) and rapid (r) neutron capture processes are the two main processes that have each forged about the half of the trans-iron nuclides. The s-process (e.g. the review of Käppeler et al. 2011) operates during the late life of asymptotic giant branch (AGB) stars (main s-process; e.g. Gallino et al. 1998; Herwig 2005; Cristallo et al. 2011; Karakas & Lattanzio 2014) and during the core helium-burning and shell carbon-burning stages of massive stars (weak s-process; e.g. Langer et al. 1989; Prantzos et al. 1990; Raiteri et al. 1991; The et al. 2007). The r-process is associated with explosive events such as neutron star mergers (e.g. Arnould et al. 2007; Goriely et al. 2011; Wanajo et al. 2014; Just et al. 2015), magnetorotational supernovae (Winteler et al. 2012; Nishimura et al. 2015), or collapsars (Siegel et al. 2019). At neutron densities in between the s- and r-processes, the existence of an intermediate (i) neutron capture process (first named by Cowan & Rose 1977) is expected. Its astrophysical site(s) is (are) actively debated (see Sect. 1 of Choplin et al. 2021, for a list of possible sites). Other nuclear processes also include short but possibly intense neutron bursts taking place in the helium shell of exploding massive stars (Blake & Schramm 1976; Thielemann et al. 1979; Meyer et al. 2004; Choplin et al. 2020). This process is thought to cause the anomalous isotopic signatures found in meteorites (Meyer et al. 2000; Pignatari et al. 2015, 2018), and possibly the abundances of some metal-poor r/s-stars (Choplin et al. 2020).

Finally, the p-process (Arnould 1976; Woosley & Howard 1978; Meyer 1994; Arnould & Goriely 2003) is thought to take place in the hydrostatic oxygen- and neon-burning shells of massive stars (Arnould 1976; Rauscher et al. 2002), during core-collapse supernovae (CCSNe; Rayet et al. 1995), or during type Ia supernovae (Travaglio et al. 2015). The p-process occurs through combinations of (γ,n), (γ,p), and (γ,α) reactions affecting pre-existing s- or r-nuclides. The p-process provides a possible explanation for the neutron-deficient nuclides in Solar System abundances (e.g. Arnould & Goriely 2003; Lugaro et al. 2016). Some discrepancies remain, however. In particular, 92, 94Mo and 96, 98Ru p-isotopes are systematically underproduced (Rayet et al. 1995; Arnould & Goriely 2003, 2020, their Fig. 22). This fact motivates the search for alternative or additional ways to produce these nuclides. In particular, the so-called pn process, i.e. a proton-poor neutron-boosted rp-process, which is encountered during He detonation, has been suggested as a promising nucleosynthesis source (Goriely et al. 2005). Such an object is made of a carbon-oxygen white dwarf with a sub-Chandrasekhar mass (M <  1.4 M) that accumulates a He-rich layer at its surface. An alternative site proposed to explain the origin of the Mo and Ru p-nuclei is the p-rich neutrino-driven wind in CCSNe, where antineutrino absorptions in the proton-rich environment produce neutrons that are immediately captured by neutron-deficient nuclei (the νp-process; Fröhlich et al. 2006; Ghosh et al. 2022). Recent results of Bliss et al. (2018) suggest that these proton-rich winds can make dominant contributions to the solar abundance of 98Ru, but that additional astrophysical sources are likely required to account for 92, 94Mo and 96Ru. The proton richness of the wind also remains a highly debated question.

Costa et al. (2000) have shown that an increase in 22Ne(α,n)25Mg reaction during stellar evolution leads to a pre-supernova seed distribution that can ultimately provide enough 92, 94Mo and 96, 98Ru p-isotopes during the explosion. Interestingly, the 22Ne(α,n)25Mg reaction can be naturally boosted in rotating massive stars through additional production of 22Ne during core-helium burning (Meynet et al. 2006; Hirschi 2007; Hirschi et al. 2008). The efficiency of the weak s-process is thus significantly boosted if the massive star is rotating (Pignatari et al. 2008; Frischknecht et al. 2016; Choplin et al. 2018; Limongi & Chieffi 2018; Banerjee et al. 2019). In rotating massive stars, the pre-supernova distribution of trans-iron elements is therefore different from that in non-rotating massive stars. This will impact the p-process that can take place during the explosion of such stars, and which largely depends on the trans-iron seeds before the explosion.

In this paper we explore how rotation impacts the p-process nucleosynthesis in exploding rotating massive stars and whether these stars can account for the abundances of solar p-nuclides. Section 2 presents the input physics relevant to the present work, Sect. 3 discusses the results obtained in terms of s-process during the hydrostatic burning stages and of p-process during both the hydrostatic and explosive evolutions. Section 4 investigates the contribution of massive rotating stars to the Galactic p-nuclei. Conclusions are given in Sect. 5.

2. Input physics

2.1. Stellar evolution models

We started from two 25 M models at a metallicity of Z = 10−3 in mass fraction, which were computed with the Geneva stellar evolution code (Eggenberger et al. 2008) and published in Choplin et al. (2018, models labelled 25S0 and 25S4 in their Table 1). The first model was non-rotating, and the second model was rotating with1υini/υcrit = 0.4. All details of the input physics can be found in Choplin et al. (2018). We recall the main physical ingredients here.

During stellar evolution, a nuclear reaction network of 737 isotopes (from hydrogen to polonium), coupled to the structure equations, was used. Opacity tables were computed with the OPAL tool2 and complemented at low temperature by opacities from Ferguson et al. (2005). Radiative mass-loss rates were taken from Vink et al. (2001) if log Teff ≥ 3.9 and from de Jager et al. (1988) otherwise. For convection, the Schwarzschild criterion was used. During the H- and He-burning phases, the size of the convective core is extended by dover = αHP, with HP the pressure scale height estimated at the Schwarzschild boundary, and α = 0.1. Rotation was treated according to the shellular scheme (Zahn 1992; Chaboyer & Zahn 1992; Maeder & Zahn 1998). The diffusion coefficients for horizontal and vertical shear were taken from Zahn (1992) and Talon & Zahn (1997), respectively. The rates of 22Ne(α, n)25Mg and 22Ne(α, γ)26Mg were taken from Longland et al. (2012). The rates of 17O(α, n)20Ne and 17O(α, γ)21Ne were taken from Best et al. (2013) if not stated otherwise.

Choplin et al. (2018) evolved the stellar models up to core O-burning ignition with the 737 isotopes network. At core O-burning ignition, most of the stellar layers ejected at the time of the supernova have reached their final state and will not be impacted during O- and Si-burning phases. However, the innermost layers (below a mass coordinate of about Mr = 3 M for the non-rotating 25 M model; Fig. 1, left panel) are impacted during these short ultimate stages. The p-process occurs in relatively deep layers into the star, and for the present work, it is therefore preferable to go further in the evolution to have a more reliable pre-supernova structure of the inner layers. We thus continued the evolution until the end of the core Si-burning stage. During these ultimate stages, the effect of rotation was switched off. This was shown to be a good approximation because the evolutionary timescale is far shorter than the rotational mixing timescale (Choplin et al. 2017, their Fig. 2). We also used the standard minimum network of the Geneva code (e.g. Hirschi 2007; Ekström et al. 2012) instead of the full 737 isotope network. This minimum network ensures a proper description of the energetics and keeps track of the main isotopes. The pre-supernova abundance profiles of some specific nuclei are illustrated in Fig. 1.

thumbnail Fig. 1.

Pre-supernova abundance profiles for the M25Zm3V0 (non-rotating, left panel) and M25Zm3V4 (rotating, υini/υcrit = 0.4, right panel) models computed with the rate of 17O(α,γ)21Ne from Best et al. (2013). The abundances of trans-iron elements are not explicitly followed from core O-burning ignition (see text for details). The dashed grey line shows the pre-supernova temperature profile, the scale of which is indicated on the right axis. The shaded grey area represents the remnant. The cyan area represents the zone in which the maximum temperature during a spherical 1051 erg explosion ranges between 1.8 and 3.7 GK (i.e. the zone in which the explosive p-process takes place).

In addition to the two reference 25 M models (with υini/υcrit = 0 and 0.4), we investigated the case of a model with the same initial properties that was a faster rotator, however, with υini/υcrit = 0.7. Unlike the first two models, the evolution of this model was not followed after core O-burning ignition. We adopted the same explosion model as for the υini/υcrit = 0.4 model, but with the important difference that the initial abundances of trans-iron elements prior to the explosion were not the same. This means that the same mass cut3 as for the υini/υcrit = 0.4 model was adopted and the temperature and density histories of the stellar layers during the explosion were the same as for the υini/υcrit = 0.4 model (cf. Sect. 2.2 for more details of the explosive modelling). This is a good approximation because the υini/υcrit = 0.4 and 0.7 models have very similar characteristics at core oxygen-burning ignition (in particular, similar temperature and density profiles, and similar CO-core masses). The only important difference between the υini/υcrit = 0.4 and 0.7 models is the abundance of trans-iron elements prior to the explosion (cf. Sect. 3.1).

For each of our three model stars, we also considered two alternative models for which the rate of the 17O(α,γ)21Ne reaction was changed, as illustrated in Fig. 2. In one case (models B), we used the rate from Best et al. (2013) divided by 10 instead of the original rate from Best et al. (2013). In a second case (models C) we adopted the lower limit from Taggart et al. (2019) instead of Best et al. (2013). Although theoretical and experimental works were carried out to study this rate (Descouvemont 1993; Best et al. 2011, 2013; Taggart et al. 2019), significant uncertainties still remain in the temperature range of interest for the s-process. The uncertainty on this rate was shown to dramatically affect the s-process efficiency in rotating massive stars (e.g. Taggart et al. 2019). A low 17O(α,γ)21Ne rate enhances the s-process efficiency because in this case, the competing 17O(α,n)20Ne reaction becomes dominant and recycles neutrons (more details in Sect. 3.1, also in Sect. 3.5 of Choplin et al. 2018). Only a theoretical estimate of 17O(α,γ)21Ne is available at low temperature (Caughlan & Fowler 1988; Best et al. 2013, black and red lines in Fig. 2). Taggart et al. (2019) constrained its lower limit experimentally. In the temperature range of interest for the s-process, it is 10 to 1000 times lower than the recommended rate of Best et al. (2013, green line in Fig. 2).

thumbnail Fig. 2.

Ratio of 17O(α,γ)21Ne reaction rates to the rate of Caughlan & Fowler (1988, cf88 label). The red curve (label bg13) is the recommended rate of Best et al. (2013). The green curve (label ta19) is the experimental lower limit from Taggart et al. (2019). The shaded area indicates the approximate range of temperatures in the helium-burning core of massive stars.

As mentioned above, in comparison with the standard models, the sets of models computed with the Best et al. (2013) rate divided by 10 or the Taggart et al. (2019) rate experience a more efficient s-process during the evolution because the 17O(α,γ)21Ne rate is lower (Fig. 2). Like for the υini/υcrit = 0.7 model with the rate of Best et al. (2013), the evolution after core O-burning ignition is not followed for these six additional models. We rely on the explosion of the corresponding models computed with the Best et al. (2013) rate. Here again, this is a good approximation because for a given υini/υcrit ratio, the models with different 17O(α,γ)21Ne rates behaves similarly. The pre-supernova abundance distributions for the different models are compared in Figs. 3 and 4 and are further discussed in Sect. 3.1. Table 1 summarizes the main characteristics of the nine models considered in this work.

thumbnail Fig. 3.

Pre-supernova mass fractions of elements heavier than iron as a function of the mass number A at a mass coordinate of 4 M for the models computed with the rate of 17O(α,γ)21Ne from Best et al. (2013, left panel) and Taggart et al. (2019, right panel). The four vertical dashed lines show the location of the s-nuclides 56Fe, 88Sr, 138Ba, and 208Pb.

thumbnail Fig. 4.

Elemental mass fractions of all models normalized by the mass fractions of the reference M25Zm3V0 model. Like in Fig. 3, the abundances are shown at a mass coordinate of 4 M. The four vertical dashed lines show the location of Fe (Z = 26), Sr (Z = 38), Ba (Z = 56), and Pb (Z = 82).

Table 1.

Characteristics of our 25 M, Z = 10−3 stellar models at the end of their evolution.

2.2. Explosion models

Two common and simple approximations for simulating the explosion of massive stars are the piston model (e.g. Woosley & Weaver 1995; Limongi et al. 2003; Heger & Woosley 2010) and the thermal bomb (e.g. Thielemann et al. 1996; Tominaga et al. 2007). The piston model usually imitates the collapse and bounce of the inner boundary. The radius of inner boundary is first reduced and then increased. In the second method, the energy is deposited at the mass cut of the exploding star in the form of internal energy. More recent methods are the so-called PUSH (e.g. Perego et al. 2015; Ebinger et al. 2019; Curtis et al. 2019) or P-HOTB (e.g. Sukhbold et al. 2016; Ertl et al. 2020) approaches, in which the mass cut and explosion energy can emerge from the simulation instead of being assumed.

In this work, the explosion energy is deposited as kinetic energy at the chosen mass cut. We used the relativistic hydrodynamical code from Tominaga et al. (2007) and Tominaga (2009). Although this code allows the computation of two dimensional explosions, we considered one-dimensional spherical explosions. A code like this was already used in Choplin et al. (2020) to investigate the effect of a jet-like explosion hitting the helium-burning shell of rotating models computed with the Geneva code.

In the present work, a spherical explosion was triggered in the non-rotating and the υini/υcrit = 0.4 models with the Best et al. (2013) rate (M25Zm3V0 and M25Zm3V4 in Table 1). The total energy deposited was Etot = 1051 erg in the standard case. Different explosion energies are considered in Sect. 3.5. The energy was deposited at a mass coordinate of 1.2 and 1.4 M for the non-rotating and rotating model, respectively. These mass coordinates correspond to the top of the 56Ni core (Fig. 1). The hydrodynamics was followed for 500 s. The temperature and density histories are recorded by mass particles representing Lagrangian mass elements of the stellar mantle.

2.3. p-process nucleosynthesis

During the last hydrostatic burning phases of massive stars and in the most inner layers, the trans-iron elements progressively photodisintegrate through a combination of (γ, n), (γ, α), and (γ, p) reactions. It was shown that some p-nuclides could form before the supernova explosion in the oxygen-burning shell and that this process is very sensitive to initial mass or to the convection model (Arnould 1976; Rauscher et al. 2002). The production of p-nuclides during the explosion takes place at higher radii than during the hydrostatic evolution (in the layers that are sufficiently heated by the supernova wave). We considered both production channels (during the evolution and during the explosion). The p-process nucleosynthesis was treated in post-processing calculations during the late hydrostatic evolution (oxygen-burning onwards) and during the explosion.

The nucleosynthesis was calculated with a code that was especially designed to follow all reactions of relevance during the p-process nucleosynthesis (Arnould & Goriely 2003). Changes in composition were followed by a full network calculation including all 2200 species from protons up to Z = 84 that lie between the proton-drip line and the neutron-rich region that may be populated. All neutron and charged-particle fusion reactions as well as their reverse reactions on elements up to Po isotopes were included. The reaction rates on light species were taken from the NETGEN library, which includes all the latest compilations of experimentally determined reaction rates (Xu et al. 2013). Experimentally unknown reaction rates were estimated with the TALYS code (Koning & Rochman 2012; Goriely et al. 2008) on the basis of the Skyrme Hartree-Fock-Bogolyubov (HFB) nuclear mass model, HFB-24 (Goriely et al. 2013), when they were not available experimentally. In addition to these reactions, electron captures and β-decays were also included. The corresponding rates were taken from experimental data (Kondev et al. 2021) when available, and from the gross theory (Tachibana et al. 1990), otherwise.

p-process nuclei are produced exclusively in layers heated at temperatures ranging between typically 1.8 and 3.5 × 109 K (Rayet et al. 1990, 1995; Arnould & Goriely 2003). Below 1.8 GK, the heavy seeds do not efficiently photodisintegrate during the evolutionary timescale. Above 3.5 GK, all seeds are photodisintegrated into iron-peak elements. For this reason, only layers with 1.8 ≲ Tmax ≲ 3.7 GK (where Tmax is the maximum temperature reached during the hydrostatic evolution and/or explosion) were post-processed in the present study. These layers are referred to as p-process layers (PPLs).

In this work, the production of a p-nucleus i is described by ⟨Fi⟩, which is the overproduction factor of this p-nucleus averaged over the PPLs, i.e.

F i = 1 M p X i , ini M 1 M 2 X i ( M r ) d M r , $$ \begin{aligned} \langle F_i \rangle = \frac{1}{M_p \, X_{i,\mathrm{ini}}} \, \int _{M_{1}}^{M_{2}} X_{i}(M_{r}) \, \mathrm{d}M_r, \end{aligned} $$(1)

where Mp = M2 − M1 is the total mass of the PPLs (these layers are delimited by the lower Lagrangian mass coordinate M1 and the upper one, M2), Xi, ini the initial abundance of the p-nucleus i and Xi(Mr) the mass fraction of p-nucleus i at mass coordinate Mr (after the explosion and beta-decays, if not stated otherwise).

3. Results

Figure 1 shows that the non-rotating model (left panel) only lost ∼0.5 M through winds during the evolution, while the rotating model ejected about ∼8 M. The main reason for this is that rotation produces larger helium-burning cores, which boosts the stellar luminosity and hence the mass loss. A more luminous star is also more likely to enter the supra-Eddington regime (cf. Sect. 3.2 in Choplin et al. 2018, for more details of this model).

3.1. s-process nucleosynthesis during hydrostatic burning

s-process nuclei represent the seed from which p-process elements are made. The s-process has been shown to be affected by rotation during core helium-burning. Consequently, the p-process will also be impacted later in the evolution or during the explosion.

As shown in Fig. 3, the s-process becomes more efficient with increasing initial rotation. This is due to the stronger operation of the rotational mixing during the core-helium burning phase. It first transports 12C and 16O from the He-core to the H-shell, which creates primary 14N. The 14N diffuses backward and penetrates the growing convective He-core. This makes primary 22Ne through the 14N(α, γ)18F(β+)18O(α, γ)22Ne chain. The neutron source 22Ne(α,n)25Mg is then boosted, and so is the s-process (e.g. Pignatari et al. 2008; Frischknecht et al. 2012).

The s-process also becomes more efficient with the 17O(α,γ)21Ne reaction rate from Taggart et al. (2019, right panel in Fig. 3) compared to the rate of Best et al. (2013, left panel). The reason is that the 16O(n,γ)17O(α,γ)21Ne chain is weaker than the competing 16O(n,γ)17O(α,n)20Ne neutron recycling chain. The latter gives the neutrons back to the medium and hence favours neutron captures by heavy seeds.

In the rotating models computed with the rate of Best et al. (2013), rotation boosts the s-process by up to a factor of 10−103 for 30 <  Z <  60 with a peak at Z ∼ 45 (solid lines in Fig. 4). For a given rotation, the uncertainties associated with the 17O(α,γ)21Ne rate change the production by a factor of typically 10, except when fast rotation and heavy elements with Z >  55 are considered, where the differences reach a factor of about 100 (the solid and dashed red patterns in Fig. 4).

3.2. p-process nucleosynthesis during hydrostatic burning

At solar metallicity, some p-nuclides can already form during the ultimate hydrostatic burning stages in the oxygen-burning shell, as initially suggested by Arnould (1976). Due to the relatively long evolution of the ultimate hydrostatic stages (with respect to the explosion timescale), the p-process during hydrostatic burning already takes place at a minimum temperature of about 1.5 GK (Rayet et al. 1990). In our M25Zm3V4 model, this temperature was reached at some point during the evolution in the layers with Mr <  3.9 M. This is illustrated by the light blue area in Fig. 5. We evaluated the p-process contribution during hydrostatic burning by post-processing the M25Zm3V4 model. We selected the layers that reached a temperature of at least 1 GK before the end of the hydrostatic evolution. These layers were post-processed with our p-process nucleosynthetic code from the onset of core oxygen-burning to the pre-supernova stage.

thumbnail Fig. 5.

Temperature profiles at different stages of the evolution (blue lines) and maximum temperature reached during the explosion (magenta line) for the M25Zm3V4 model. The abundance profile of the s-nuclide 88Sr is shown in black with the scale on the right. The blue and magenta rectangle show where the p-process can take place (in terms of temperature range and mass coordinate range) during the hydrostatic and explosive burning, respectively (see text for details). The dark grey area shows the extent of the remnant (corresponding to the size of the 56Ni core; cf. Sect. 2.2). The three horizontal dashed lines correspond to temperatures of 1.5, 1.8, and 3.7 GK.

Figure 6 shows the abundance profiles of the s-nuclide 88Sr (black) and p-nuclide 74Se (red) at different burning stages. At the start of core oxygen burning, only 88Sr is abundant because 74Se was fully destroyed by the s-process during core helium-burning (dot-dashed black line in Fig. 6). During the ultimate hydrostatic phases, 88Sr is partially photodisintegrated at Mr <  2.4 M (dashed black line) and 74Se is synthesized around Mr ∼ 2.2 M (dashed red line) through the photodisintegration of s-nuclides.

thumbnail Fig. 6.

Abundance profiles of the s-nuclide 88Sr (black) and the p-nuclide 74Se (red) in the region of the star that is relevant for the p-process. The abundance profiles are shown at four different stages: at the birth of the star (or zero-age main-sequence, dotted lines), at the start of core oxygen burning (dash-dot lines, at this point, the 74Se abundance is about zero, hence not visible on the plot), at the end of the hydrostatic post-processing calculation (dashed lines) and after the supernova explosion (solid lines).

Figure 7 shows the results of p-process nucleosynthesis when only the hydrostatic evolution (green), only the explosion (blue, see Sect. 3.3 for details; in this case, the abundances of heavy seeds are those at the start of the core oxygen-burning phase), or both the hydrostatic plus explosive burning (red) are considered. The contribution stemming from the hydrostatic burning phase is seen to be significant and smaller than the explosive contribution by a factor of 2−5 for most nuclides. This matter is further processed by the explosive nucleosynthesis. Considering only the explosive nucleosynthesis is seen to be an excellent approximation to the total hydrostatic plus explosive processing. It shows that at least in this model, p-nuclides are mostly synthesized during the explosion. In other conditions, and more particularly, in very massive (M ≳ 100 M) stars leading to pair creation supernovae, hydrostatic Ne/O burning has been shown to provide a dominant component to the total production of p-nuclei (Rayet et al. 1993). As proposed by Arnould & Goriely (2003), the p-process may also be found to develop in multi-dimensionally simulated pre-supernova O-rich shells where additional convective mixing is induced by gravity waves (Arnett 2001). In this case, this extended mixing may also move the p-nuclides that are synthesized before the explosion to external stellar regions, in which they might survive the supernova explosion, in contrast to the pattern found in one-dimensional simulations.

thumbnail Fig. 7.

Overproduction factors ⟨Fi⟩ (Eq. (1)) of p-nuclei mass-averaged over the 4.25 M PPLs (corresponding to the union of the blue and magenta boxes in Fig. 5) in the M25Zm3V4 model at the end of the hydrostatic burning pre-supernova phase (green diamonds), only due to the explosive burning (blue) or including both the hydrostatic and explosive burning phases (red; see text for more details).

3.3. p-process nucleosynthesis during the explosion

We considered explosions of total energy Etot = 1051 erg (the impact of the explosion energy is discussed in Sect. 3.5). The energy was deposited at a mass cut located at the top of the iron (56Ni) core. During the explosion, the p-process took place in PPLs reaching a peak temperature between ∼1.8 and ∼3.5 GK and corresponding to mass coordinates of 2.59 <  Mr <  5.44 M in our M25Zm3V4 model (magenta area in Fig. 5). The temperature and density evolution of several PPLs verifying this criteria is shown in Fig. 8. The relatively high abundance of heavy s-process seeds in this region (represented by the s-nucleus 88Sr, black line in Fig. 5) gives rise to an efficient p-process during the explosion. This is visible in Fig. 6 for the specific cases of 88Sr and 74Se. As other s-nuclides, the pre-supernova 88Sr abundance (dashed black line) is high above Mr = 2.5 M. This enables an abundant production of p-nuclides such as 74Se during the explosion (solid red line). We note that the zone in the star in which the explosive p-process takes place is above the zone in which the hydrostatic p-process takes place (compare the blue and magenta area in Fig. 5) because the maximum temperature reached during the explosion is higher than in the hydrostatic evolution. The overlap between these regions is about 1.3 M. This region contains p-nuclides produced by both hydrostatic and explosive burning. In the inner blue zone (around Mr = 2 M for instance), the p-nuclides built during hydrostatic burning are all destroyed by photodisintegration due to the high temperatures (Tex, max >  3.7 GK) encountered during the explosion.

thumbnail Fig. 8.

Temperature (left panel) and density (right panel) evolution of several tracer particles coming from the M25Zm3V4 explosion model. The maximum temperature of these particles during the explosion verifies 1.8 <  Tex, max <  3.7 GK. These particles are initially located in the cyan area in Fig. 1 (right panel).

The s-process efficiency during the evolution increases with initial rotation and with decreasing 17O(α, γ) rate. Consequently, the p-process efficiency, which very strongly depends on the initial trans-iron seeds content, also increases with initial rotation and when the 17O(α, γ) rate is lowered (Fig. 9). For p-nuclides with 90 ≲ A ≲ 150, the overproduction differences between the non-rotating and the fast-rotating models reach about 3 dex at most (Fig. 9, black and red curves). This is similar for the differences shown in Figs. 3 and 4 (black and red curves) for the s-process. For p-nuclides with A ≲ 90 and A ≳ 150, the differences decrease, as expected from the pre-supernova s-process yields shown in Figs. 3 and 4: they reach about 1 dex at most. A lower 17O(α, γ) rate increases the s-process efficiency (Figs. 3 and 4), and thus also the p-process (Fig. 9). In particular, it leads to the production of a substantial amount of p-nuclides with A ≳ 100, as expected from the s-process yields shown in Figs. 3 and 4.

thumbnail Fig. 9.

Overproduction factors ⟨Fi⟩ (Eq. (1)) mass-averaged over all PPLs for all p-nuclides when considering the 17O(α,γ)21Ne rate of Best et al. (2013, top panel), Best et al. (2013) divided by 10 (middle), and Taggart et al. (2019, bottom panel). The filled symbols highlight the npp-nuclides (Table 2) that have the highest overabundances, i.e. those with ⟨Fi⟩ higher than the highest value divided by an arbitrary factor of 20. The green curves show the velocity-averaged yields over the three 25 M models assuming they follow the rotational distribution of observed young B stars from Huang et al. (2010, their Fig. 6). See text for more details.

Figure 10 illustrates the overproduction factors in the rotating 25 M models with υini/υcrit = 0.4 and 0.7, not only for p-nuclei, but also for all stable nuclei. The p-nuclei appear to be rather well co-produced with the α-elements such as 16O, 20Ne, 24Mg, or 28Si. Although the overabundances shown in Fig. 10 only represent those characterizing the PPLs, the integration over the entire stellar mass gives an overproduction of light p-elements (A = 74−84 for M25Zm3V4 and A = 74−132 for M25Zm3V7) similar to those of the light α-elements. By contrast, the non-rotating model gives an overproduction factor of about 100 in the PPLs, i.e. about ten times smaller than the factor obtained in such stars for 16O (black symbols in Fig. 10). This difference is due to the low initial metallicity of about Z/10 adopted here.

thumbnail Fig. 10.

Overproduction factors ⟨Fi⟩ (Eq. (1)) mass-averaged over all PPLs for all stable nuclei (open symbols) and p-nuclides (full symbols) for the 25 M models with the 17O(α,γ)21Ne rate from Best et al. (2013) and with υini/υcrit = 0 (black diamonds), 0.4 (blue squares) and 0.7 (red circles). The overproduction factors of some α-elements, including 16O, are highlighted.

Finally, we note that the emitted flux of (anti)neutrinos is not followed in our explosion simulation. For this reason, neutrino-induced nucleosynthesis that can boost the production of some rare p-nuclei, such as the odd-odd 138La and 180Ta (Goriely et al. 2001; Sieverding et al. 2018), is not described in the present work.

3.4. 92, 94Mo and 96, 98Ru p-nuclides

As shown in Fig. 9 and found systematically by previous calculations (see e.g. Rayet et al. 1995; Arnould & Goriely 2003; Travaglio et al. 2018), the p-nuclides 92, 94Mo and 96, 98Ru are always underproduced by non-rotating models. In particular, the overproduction factors of Mo-Ru p-isotopes are smaller by a factor of about 100 than the overproduction factors of the light p-nuclei 74Se, 78Kr, and 84Sr. For rotating models, this deviation is seen to be substantially reduced to a ratio of about 20 and 3 for the models with υini/υcrit = 0.4 and 0.7, respectively. In the M25Zm3V7B model with a lower 17O(α,γ)21Ne rate (Fig. 9, middle panel, red pattern), the Mo-Ru isotopes are overproduced at the same level as the light p-nuclei, but the 100 ≲ A ≲ 132 p-nuclides are now overproduced with respect to the light nuclides by a factor of about 10. The increased production of Mo-Ru p-nuclides is directly linked to the s-enrichment in heavy seeds such as Ba during the evolution of rotating models (cf. Sect. 3.1).

3.5. Impact of explosion energy on the p-process

To obtain an efficient explosive p-process, a region in the star must be enriched in trans-iron seeds and experience temperatures between 1.8 and 3.7 GK during the explosion. In the M25Zm3V4 model, the region enriched in trans-iron seeds extends to about Mr = 8 M (Figs. 5 and 11). Above this, the abundances of trans-iron elements are very similar to their initial abundances. The extension of the 1.8 <  Tex, max <  3.7 GK region is shown in Fig. 11 for various explosion energies. As expected, this zone shifts upward in mass for more energetic explosions. Nevertheless, for explosion energies 0.3 × 1051 <  Etot <  5 × 1052 erg, this zone stays in the stellar region that is enriched in trans-iron seeds. This implies that an efficient p-process takes place for a whole range of explosion energies in this model. Moreover, the yields will be weakly impacted by the explosion energy because the pre-supernova trans-iron seed abundances vary weakly (by less than a factor of 2) up to Mr = 8 M.

thumbnail Fig. 11.

Mass coordinate as a function of explosion energy for the M25Zm3V4 model. The color code indicates the mass fraction of the s-nuclide 88Sr prior to the explosion. The two black lines delimit the zone in which the maximum temperature during the explosion verifies 1.8 <  Tex, max <  3.7 GK. The dark grey rectangular zone shows the extent of the remnant. The dashed line shows the mass coordinate below which some photodisentegration of trans-iron seeds can have occurred during the last hydrostatic burning stages (zone in which T >  1.5 GK at some point during hydrostatic burning).

The resulting abundance distribution of p-nuclides assuming two different values of the explosion energy, namely Etot = 1051 and 1052 erg, are illustrated in Fig. 12. Both distributions are globally identical, except for the lightest 74Se and 78Kr and the A = 180 Ta and W p-isotopes. As shown in Fig. 11, for Etot = 1052 erg, the PPLs extend up to the region in which s-process overabundances are highest. For these outer layers, the peak temperature remains relatively low (Tex, max ≃ 1.8 GK), and this consequently essentially boosts the production of the heaviest p-nuclei.

thumbnail Fig. 12.

Overproduction factors ⟨Fi⟩ (Eq. (1)) mass-averaged over all PPLs for all p-nuclides for the M25Zm3V4 model with two values of the total explosion energy Etot = 1051 and 1052 erg.

We note that non-spherical explosions can lead to a wide variety of temperature and density conditions compared to the spherical explosions considered in this work. This could ultimately impact the yields of p-nuclides. We plan to investigate the effect of the explosion geometry on the p-process nucleosynthesis in a future work.

3.6. Velocity-averaged yields for a 25 M star

To make a first rough estimate of the integrated p-nuclide yields of a population of rotating massive stars, we can use the velocity distribution derived by Huang et al. (2010, their Fig. 6). Their velocity distribution was obtained from the observation of 220 young main-sequence B-type stars. We divided their υini/υcrit probability density function into three intervals corresponding to our three 25 M models with υini/υcrit = 0, 0.4, and 0.7. The first interval extends from υini/υcrit = 0 to 0.2, the second from 0.2 to 0.55 (midway between 0.4 and 0.7), and the third from 0.55 to 1. We then integrated their probability density function over these three intervals to obtain the relative weights of 0.14, 0.45, and 0.41 for the υini/υcrit = 0, 0.4, and 0.7 models, respectively. The green curves in Fig. 9 (also reported in Fig. 13) show the resulting overproduction factors when weighting the yields of our three models by the derived coefficients. In each of the three panels of Fig. 9, the final yields are mainly influenced by the fast-rotating υini/υcrit = 0.7 model, which indeed has (i) high overproduction factors and (ii) a relatively high weight of 0.41. Overall, at this mass and metallicity, the p-nuclide yields of a fast-rotating (typically υini/υcrit = 0.7) massive star are likely well representative of a population of massive stars with various initial rotation rates. A more detailed estimation including a finer initial rotational grid, different masses, and metallicities, is deferred to a future work.

thumbnail Fig. 13.

Same as Fig. 9, where only the overproduction factors for the three velocity-averaged models are compared.

4. Contribution of massive rotating stars to Galactic p-nuclei

Since oxygen is known to be essentially produced by CCSNe, the contribution of massive stars to the Galactic p-nuclei enrichment can be estimated through the oxygen to p-nuclide yield ratio normalized to the solar abundances, (O/p), as defined by Rayet et al. (1995). The calculation of this quantity requires the estimate of the average overproduction factor F0 over the 35 p-nuclides given by

F 0 = i F i / 35 , $$ \begin{aligned} F_0 = \sum _i \, \langle F_i \rangle /35, \end{aligned} $$(2)

where ⟨Fi⟩ is the overproduction factor of p-nuclide i averaged over the PPLs (Eq. (1)). We note that we considered only the explosive contribution of the p-process to estimate the (O/p) ratios (which is a good approximation, as shown in Fig. 7). Therefore, the PPLs mass Mp used in Eq. (1) and given in Table 2 corresponds to the zone in the star that experiences p-process only during the explosion (magenta zone in Fig. 5).

Table 2.

Quantities required to estimate the (O/p) ratio corresponding to the oxygen to p-nuclide yield ratio normalized to the solar abundances (Eq. (6)).

The net yield yi of a p-nuclide i is defined as the difference between the mass of a p-nuclide i returned by the star to the interstellar medium and the mass of that same p-nuclide engulfed at the birth of the star. In our case, it can be simplified as follows:

y i = X i , ini F 0 M p X i , ini M CO , $$ \begin{aligned} y_{i} = X_{i,\mathrm{ini}} \, F_0 \, M_p - X_{i,\mathrm{ini}} \, M_{\rm CO}, \end{aligned} $$(3)

where MCO is the carbon-oxygen core mass, which also corresponds to the helium-free core (Table 2). As in Rayet et al. (1995), Eq. (3) assumes that the p-nuclei initially present in the helium-burning core have been destroyed by the s-process and that their production in the PPLs can be averaged by the same overproduction factor F0.

To derive the (O/p) ratio, a similar estimate of the net oxygen yield can be derived by

y O = Y O X O , ini M ini , $$ \begin{aligned} y_{\rm O} = \mathcal{Y} _{\rm O} - X_{\rm O,ini} \, M_{\rm ini}, \end{aligned} $$(4)

where 𝒴O is the absolute oxygen yield (i.e. the total mass of oxygen ejected through stellar winds and supernova, as reported in Table 2), XO, ini is the initial oxygen mass fraction, and Mini is the initial mass. The supernova ejecta enrich the interstellar medium with oxygen and p-nucleus i in solar proportion if

y O y i = X O , X i , , $$ \begin{aligned} \frac{y_{\rm O}}{y_{i}} = \frac{X_{\rm O,\odot }}{X_{i,\odot }}, \end{aligned} $$(5)

where XO, ⊙ and Xi, ⊙ are the solar mass fraction of oxygen and p-nuclide i, respectively. Combining Eq. (5) with Eqs. (3) and (4) gives the (O/p)0 ratio for a given averaged p-enrichment, given by F0 as

( O / p ) 0 y O y i X i , X O , = Y O / X O , f Z M ini f Z F 0 M p f Z M CO , $$ \begin{aligned} (\mathrm{O}/p)_0 \equiv \frac{y_{\rm O}}{y_{i}}\frac{X_{i,\odot }}{X_{\rm O,\odot }} = \frac{\mathcal{Y} _{\rm O}/X_{\rm O,\odot } - f_Z M_{\rm ini}}{f_Z F_0 M_p - f_Z M_{\rm CO}}, \end{aligned} $$(6)

where fZ = Z/Z = 0.001/0.014 = 0.071 in our case. We note that considering a simple average overproduction factor F0 may not be optimal in our case because of the relatively large dispersion of the overproduction factors ⟨Fi⟩ in some of our models (Fig. 9). Therefore we also define an averaged overproduction factor F1, which in contrast to F0 only considers the most abundant p-nuclei. This new factor F1 is defined by an equation similar to Eq. (2), but the sum runs only over the npp-nuclides (given in Table 2), the mass fraction of which is greater than Xmax/20, where Xmax is the mass fraction of the most abundant p-nucleus. These npp-nuclides are highlighted with filled symbols in Fig. 9. The factor of 20 considered here is purely arbitrary and aims at considering the various uncertainties associated with both the astrophysical modelling (including those affecting the s-process in massive stars) and the nuclear ingredients (see in particular Fig. 35 of Arnould & Goriely 2003). Therefore, we define two ratios (O/p)0 and (O/p)1 depending on whether the (O/p) ratio in Eq. (6) is calculated with F0 or F1 (Table 2). We also computed the (O/p) ratios for the velocity-averaged models, as reported in the fourth row of each series of models in Table 2.

If (O/p) = 1, the p-nuclides are produced in solar proportion with oxygen on average. However, if our models predict (O/p) > 1, p-nuclides are predicted to be underproduced compared to oxygen. This would imply that our models alone cannot fully account for the p-nuclide enrichment in the Solar System, and, more generally, in the Galaxy. As shown in Table 2, 0.17 < (O/p)0 <  61.3 and 0.12 < (O/p)1 <  36.9. The (O/p) ratios rapidly drop with increasing rotation and with decreasing 17O(α, γ) rate, meaning that more p-nuclei are produced compared to oxygen if massive stars are rotating or if the 17O(α, γ) rate is lower (as already discussed in Sect. 3).

Rayet et al. (1995) considered solar metallicity non-rotating massive stars with initial masses ranging between 13 and 25 M. They found 1.8 < (O/p)0 <  8.4 and a ratio of 4.2 when they integrated over the initial mass function from Kroupa et al. (1993). In our case, the (O/p)0 ratios for non-rotating models are significantly higher, with values ranging between 6.7 and 61.3 (Table 2). Nevertheless, their results are compatible with ours in view of the different metallicities considered. Oxygen is a primary product, in contrast to s- and p-nuclei, which are secondary in non-rotating models and are hence metallicity dependent. Rayet et al. (1995) concluded that about one quarter of the solar system p-nuclei could be attributed to supernovae from massive stars (assuming that the entire solar oxygen is coming from such events). Similarly, on the basis of a Galactic chemical evolution model, Travaglio et al. (2018) estimated that CCSNe from populations of non-rotating massive stars with different initial masses and metallicities could contribute no more than 10% to the Galactic p-enrichment, with only a few exceptions (e.g. the light p-nuclides 74Se, 78Kr, and 84Sr). The contribution stemming from models with sub-solar metallicities was also found to be negligible due to the secondary nature of the s- and p-processes.

Our results with (O/p) ratios close to or even below one (Table 2) suggest that sub-solar-metallicity rotating massive stars can co-produce or even overproduce p-nuclei with respect to oxygen. The impact of rotation on the s-process was found to be so large in stars with a metallicity Z ≃ 0.001 (Frischknecht et al. 2016; Choplin et al. 2018, 2020) that it gives rise to large overproduction factors of p-nuclei (see Table 2 and Fig. 9) that can compensate for their sub-solar metallicity. Table 5 in Frischknecht et al. (2016) shows that the overproduction factors of light s-nuclides for their solar metallicity models are ∼50 at most and that rotation increases them by no more than a factor of about 2. The overproduction factors of p-nuclides for these models therefore is not significantly affected by rotation and is characterized by averaged values of F0 ≃ 100 (see e.g. Table 3 in Rayet et al. 1995). The averaged overproduction factors of our models with respect to the solar abundances F0, ⊙ and F1, ⊙ can also be estimated from our Table 2 by simply dividing the F0 and F1 factors by the Z/Z = 14 ratio. Considering the velocity-averaged yields of p-nuclides, overproduction factors of F0, ⊙ ≃ 50 and F1, ⊙ ≃ 70 are found when considering the rates of Best et al. (2013). When using the rate of Best et al. (2013) divided by 10, we find F0, ⊙ ≃ 280 and F1, ⊙ ≃ 350. Finally, F0, ⊙ ≃ 520 and F1, ⊙ ≃ 700 if using the rate of Taggart et al. (2019). These value are similar to or higher than those expected from solar metallicity non-rotating and rotating massive stars. On this basis, sub-solar-metallicity rotating stars may even be dominant contributors to the Galactic enrichment of p-nuclei (and s-nuclei) compared to solar metallicity (rotating) stars. Detailed Galactic evolution models are needed to quantify their global contribution. Such a quantitative study is postponed to a future work.

5. Summary and conclusions

We studied the p-process in rotating massive stars during their ultimate hydrostatic burning stages and during their explosions. These stars can experience an enhanced s-process during their evolution due to the effect of rotational mixing. Consequently, they are enriched in trans-iron seeds in the last evolutionary stages and during the explosion. Because of the higher abundance in heavy seeds, they are expected to experience an efficient p-process nucleosynthesis during their explosion and possibly during their last hydrostatic burning stages.

We computed 25 M stellar models at a metallicity of Z = 10−3 with different initial rotation velocities and 17O(α,γ)21Ne rates. We found that most of the p-nuclides are synthesized during the explosion. The impact of rotation on the p-process follows the effect of rotation on the s-process. From no to fast rotation, both the s- and p-process efficiencies are boosted by about 3−4 dex for nuclides with A <  140. The impact remains small for nuclides with A ≥ 140, unless both fast rotation and a lower 17O(α,γ)21Ne reaction rate are considered. In this case, a significant number of heavy p-nuclides (with A ≥ 140) is produced. The impact of the explosion energy on the p-process yields remains weak because a different explosion energy just shifts the zone in the star in which the p-process takes place. Because this zone is similarly enriched in s-nuclides (built during the previous evolutionary stages), the p-process yields are weakly impacted.

By considering a population of solar-metallicity non-rotating massive stars, Rayet et al. (1995) and Travaglio et al. (2018) concluded that CCSNe from populations of non-rotating massive stars could contribute no more than 10−25% to the Galactic p-enrichment, and as in previous works, they reported that sub-solar-metallicity stars play a minor role. However, their studies only considered non-rotating models and may consequently have underestimated the role of CCSNe in the Galactic enrichment of p-nuclei. Our present results suggest that rotating massive stars with sub-solar metallicity may substantially contribute to the Galactic p-enrichment, and that the global contribution stemming from CCSNe needs to be revisited. In particular, rotating massive stars with sub-solar metallicity may be the dominant contributors compared to solar metallicity massive stars because the effect of rotational mixing at sub-solar metallicity is stronger. However, before drawing quantitative conclusions, more detailed studies are required that use a larger grid of masses, metallicities, and rotation velocities, supplemented by Galactic chemical evolution simulations.

The p-process efficiency in rotating massive stars also remains very strongly impacted by the uncertainty associated with the 17O(α,γ)21Ne reaction rate and its effects on the s-process seed distribution. An accurate determination of this rate will greatly help constrain the contribution of massive rotating stars to the Galactic content not only in s-nuclei, but also in p-nuclei.

As mentioned at the end of Sect. 3.5, the p-process nucleosynthesis in self-consistent multi-dimensional explosions from rotating progenitors also remains to be assessed. These models are thought to experience r-process nucleosynthesis (Winteler et al. 2012; Nishimura et al. 2015; Mösta et al. 2018; Reichert et al. 2021). Our simulations suggest that exploding rotating massive stars can be rich nucleosynthesis astrophysical sites for the production of elements heavier than iron, in particular by enriching the interstellar medium in (1) s-process material from their outer He-rich layers, (2) p-process material from deep O-rich layers, and (3) some r-process material from the innermost ejecta.

Finally, our study once again highlights the important impact of mixing mechanisms in deep stellar interiors on the nucleosynthesis in general, and of the p-process in particular. The only way to substantially improve the current situation in this respect is to master multi-dimensional star simulations in future developments.


1

The critical velocity υcrit is reached when gravitational acceleration is counterbalanced by centrifugal force. In the Roche approximation, it is expressed as υ crit = 2 3 GM R p , c $ \upsilon_{\mathrm{crit}} = \sqrt{\frac{2}{3}\frac{GM}{R_{\mathrm{p,c}}}} $, with Rp, c the polar radius at the critical limit.

3

At the time of the explosion, the mass cut is the mass coordinate that delimits the part of the star that is expelled from the part that is locked into the remnant.

Acknowledgments

We wish to thank M. Arnould for his careful reading and valuable remarks. This work was supported by the Fonds de la Recherche Scientifique-FNRS under Grant No IISN 4.4502.19. S.G. is FRS-F.N.R.S. research associate. R.H. acknowledges support from the World Premier International Research Centre Initiative (WPI Initiative, MEXT, Japan), STFC UK, the European Union’s Horizon 2020 research and innovation programme under grant agreement No 101008324 (ChETEC-INFRA) and the IReNA AccelNet Network of Networks, supported by the National Science Foundation under Grant No. OISE-1927130. R.H. also acknowledges support from the ChETEC COST Action (CA16117), supported by COST (European Cooperation in Science and Technology). G.M. has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 833925, project STAREX).

References

  1. Arnett, W. 2001, Nucl. Phys. A, 688, 177 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arnould, M. 1976, A&A, 46, 117 [Google Scholar]
  3. Arnould, M., & Goriely, S. 2003, Phys. Rep., 384, 1 [Google Scholar]
  4. Arnould, M., & Goriely, S. 2020, Progr. Part. Nucl. Phys., 112, 103766 [NASA ADS] [CrossRef] [Google Scholar]
  5. Arnould, M., Goriely, S., & Takahashi, K. 2007, Phys. Rep., 450, 97 [Google Scholar]
  6. Banerjee, P., Heger, A., & Qian, Y.-Z. 2019, ApJ, 887, 187 [Google Scholar]
  7. Best, A., Görres, J., Couder, M., et al. 2011, Phys. Rev. C, 83, 052802 [NASA ADS] [CrossRef] [Google Scholar]
  8. Best, A., Beard, M., Görres, J., et al. 2013, Phys. Rev. C, 87, 045805 [NASA ADS] [CrossRef] [Google Scholar]
  9. Blake, J. B., & Schramm, D. N. 1976, ApJ, 209, 846 [Google Scholar]
  10. Bliss, J., Arcones, A., & Qian, Y. Z. 2018, ApJ, 866, 105 [NASA ADS] [CrossRef] [Google Scholar]
  11. Caughlan, G. R., & Fowler, W. A. 1988, At. Data Nucl. Data Tab., 40, 283 [NASA ADS] [CrossRef] [Google Scholar]
  12. Chaboyer, B., & Zahn, J.-P. 1992, A&A, 253, 173 [Google Scholar]
  13. Choplin, A., Ekström, S., Meynet, G., et al. 2017, A&A, 605, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Choplin, A., Hirschi, R., Meynet, G., et al. 2018, A&A, 618, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Choplin, A., Tominaga, N., & Meyer, B. S. 2020, A&A, 639, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Choplin, A., Siess, L., & Goriely, S. 2021, A&A, 648, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Costa, V., Rayet, M., Zappalà, R. A., & Arnould, M. 2000, A&A, 358, L67 [NASA ADS] [Google Scholar]
  18. Cowan, J. J., & Rose, W. K. 1977, ApJ, 212, 149 [Google Scholar]
  19. Cristallo, S., Piersanti, L., Straniero, O., et al. 2011, ApJS, 197, 17 [NASA ADS] [CrossRef] [Google Scholar]
  20. Curtis, S., Ebinger, K., Fröhlich, C., et al. 2019, ApJ, 870, 2 [Google Scholar]
  21. de Jager, C., Nieuwenhuijzen, H., & van der Hucht, K. A. 1988, A&AS, 72, 259 [NASA ADS] [Google Scholar]
  22. Descouvemont, P. 1993, Phys. Rev. C, 48, 2746 [NASA ADS] [CrossRef] [Google Scholar]
  23. Ebinger, K., Curtis, S., Fröhlich, C., et al. 2019, ApJ, 870, 1 [NASA ADS] [Google Scholar]
  24. Eggenberger, P., Meynet, G., Maeder, A., et al. 2008, Ap&SS, 316, 43 [Google Scholar]
  25. Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146 [Google Scholar]
  26. Ertl, T., Woosley, S. E., Sukhbold, T., & Janka, H. T. 2020, ApJ, 890, 51 [CrossRef] [Google Scholar]
  27. Ferguson, J. W., Alexander, D. R., Allard, F., et al. 2005, ApJ, 623, 585 [Google Scholar]
  28. Frischknecht, U., Hirschi, R., & Thielemann, F.-K. 2012, A&A, 538, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Frischknecht, U., Hirschi, R., Pignatari, M., et al. 2016, MNRAS, 456, 1803 [Google Scholar]
  30. Fröhlich, C., Martínez-Pinedo, G., Liebendörfer, M., et al. 2006, Phys. Rev. Lett., 96, 142502 [CrossRef] [Google Scholar]
  31. Gallino, R., Arlandini, C., Busso, M., et al. 1998, ApJ, 497, 388 [Google Scholar]
  32. Ghosh, S., Wolfe, N., & Fröhlich, C. 2022, ApJ, 929, 43 [NASA ADS] [CrossRef] [Google Scholar]
  33. Goriely, S., Arnould, M., Borzov, I., & Rayet, M. 2001, A&A, 375, L35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Goriely, S., Garcia-Senz, D., Bravo, E., & José, J. 2005, A&A, 444, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Goriely, S., Hilaire, S., & Koning, A. J. 2008, A&A, 487, 767 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Goriely, S., Bauswein, A., & Janka, H.-T. 2011, ApJ, 738, L32 [Google Scholar]
  37. Goriely, S., Chamel, N., & Pearson, J. 2013, Phys. Rev. C, 88, 024308 [NASA ADS] [CrossRef] [Google Scholar]
  38. Heger, A., & Woosley, S. E. 2010, ApJ, 724, 341 [Google Scholar]
  39. Herwig, F. 2005, ARA&A, 43, 435 [NASA ADS] [CrossRef] [Google Scholar]
  40. Hirschi, R. 2007, A&A, 461, 571 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Hirschi, R., Chiappini, C., Meynet, G., Maeder, A., & Ekström, S. 2008, in Massive Stars as Cosmic Engines, eds. F. Bresolin, P. A. Crowther, & J. Puls, 250, 217 [NASA ADS] [Google Scholar]
  42. Huang, W., Gies, D. R., & McSwain, M. V. 2010, ApJ, 722, 605 [Google Scholar]
  43. Just, O., Bauswein, A., Ardevol Pulpillo, R., Goriely, S., & Janka, H.-T. 2015, MNRAS, 448, 541 [Google Scholar]
  44. Käppeler, F., Gallino, R., Bisterzo, S., & Aoki, W. 2011, Rev. Mod. Phys., 83, 157 [Google Scholar]
  45. Karakas, A. I., & Lattanzio, J. C. 2014, PASA, 31, e030 [NASA ADS] [CrossRef] [Google Scholar]
  46. Kondev, F., Wang, M., Huang, W., Naimi, S., & Audi, G. 2021, Chin. Phys. C, 45, 030001 [NASA ADS] [CrossRef] [Google Scholar]
  47. Koning, A. J., & Rochman, D. 2012, Nucl. Data Sheets, 113, 2841 [NASA ADS] [CrossRef] [Google Scholar]
  48. Kroupa, P., Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545 [NASA ADS] [CrossRef] [Google Scholar]
  49. Langer, N., Arcoragi, J.-P., & Arnould, M. 1989, A&A, 210, 187 [Google Scholar]
  50. Limongi, M., & Chieffi, A. 2018, ApJS, 237, 13 [NASA ADS] [CrossRef] [Google Scholar]
  51. Limongi, M., Chieffi, A., & Bonifacio, P. 2003, ApJ, 594, L123 [CrossRef] [Google Scholar]
  52. Longland, R., Iliadis, C., & Karakas, A. I. 2012, Phys. Rev. C, 85, 065809 [CrossRef] [Google Scholar]
  53. Lugaro, M., Pignatari, M., Ott, U., et al. 2016, Proc. Natl. Acad. Sci., 113, 907 [NASA ADS] [CrossRef] [Google Scholar]
  54. Maeder, A., & Zahn, J.-P. 1998, A&A, 334, 1000 [NASA ADS] [Google Scholar]
  55. Meyer, B. S. 1994, ARA&A, 32, 153 [NASA ADS] [CrossRef] [Google Scholar]
  56. Meyer, B. S., Clayton, D. D., & The, L.-S. 2000, ApJ, 540, L49 [NASA ADS] [CrossRef] [Google Scholar]
  57. Meyer, B. S., The, L. S., Clayton, D. D., & El Eid, M. F. 2004, in Lunar and Planetary Science Conference, eds. S. Mackwell, & E. Stansbery, 1908 [Google Scholar]
  58. Meynet, G., Ekström, S., & Maeder, A. 2006, A&A, 447, 623 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Mösta, P., Roberts, L. F., Halevi, G., et al. 2018, ApJ, 864, 171 [Google Scholar]
  60. Nishimura, N., Takiwaki, T., & Thielemann, F.-K. 2015, ApJ, 810, 109 [Google Scholar]
  61. Perego, A., Hempel, M., Fröhlich, C., et al. 2015, ApJ, 806, 275 [NASA ADS] [CrossRef] [Google Scholar]
  62. Pignatari, M., Gallino, R., Meynet, G., et al. 2008, ApJ, 687, L95 [Google Scholar]
  63. Pignatari, M., Zinner, E., Hoppe, P., et al. 2015, ApJ, 808, L43 [NASA ADS] [CrossRef] [Google Scholar]
  64. Pignatari, M., Hoppe, P., Trappitsch, R., et al. 2018, Geochim. Cosmochim. Acta, 221, 37 [NASA ADS] [CrossRef] [Google Scholar]
  65. Prantzos, N., Hashimoto, M., & Nomoto, K. 1990, A&A, 234, 211 [NASA ADS] [Google Scholar]
  66. Raiteri, C. M., Busso, M., Picchio, G., & Gallino, R. 1991, ApJ, 371, 665 [Google Scholar]
  67. Rauscher, T., Heger, A., Hoffman, R. D., & Woosley, S. E. 2002, ApJ, 576, 323 [Google Scholar]
  68. Rayet, M., Arnould, M., & Prantzos, N. 1990, A&A, 227, 271 [NASA ADS] [Google Scholar]
  69. Rayet, M., Eid, M. E., & Arnould, M. 1993, in Nuclei in the Cosmos, eds. F. Käppeler, & K. Wisshak (Bristol: Institute of Physics Publishing), 613 [Google Scholar]
  70. Rayet, M., Arnould, M., Hashimoto, M., Prantzos, N., & Nomoto, K. 1995, A&A, 298, 517 [Google Scholar]
  71. Reichert, M., Obergaulinger, M., Eichler, M., Aloy, M. Á., & Arcones, A. 2021, MNRAS, 501, 5733 [NASA ADS] [Google Scholar]
  72. Siegel, D. M., Barnes, J., & Metzger, B. D. 2019, Nature, 569, 241 [Google Scholar]
  73. Sieverding, A., Martinez-Pinedo, G., Huther, L., Langanke, K., & Heger, A. 2018, ApJ, 865, 143 [NASA ADS] [CrossRef] [Google Scholar]
  74. Sukhbold, T., Ertl, T., Woosley, S. E., Brown, J. M., & Janka, H. T. 2016, ApJ, 821, 38 [NASA ADS] [CrossRef] [Google Scholar]
  75. Tachibana, T., Yamada, M., & Yoshida, Y. 1990, Prog. Theor. Phys., 84, 641 [CrossRef] [Google Scholar]
  76. Taggart, M. P., Akers, C., Laird, A. M., et al. 2019, Phys. Lett. B, 798, 134894 [NASA ADS] [CrossRef] [Google Scholar]
  77. Talon, S., & Zahn, J.-P. 1997, A&A, 317, 749 [Google Scholar]
  78. The, L.-S., El Eid, M. F., & Meyer, B. S. 2007, ApJ, 655, 1058 [Google Scholar]
  79. Thielemann, F.-K., Arnould, M., & Hillebrandt, W. 1979, A&A, 74, 175 [Google Scholar]
  80. Thielemann, F.-K., Nomoto, K., & Hashimoto, M.-A. 1996, ApJ, 460, 408 [NASA ADS] [CrossRef] [Google Scholar]
  81. Tominaga, N. 2009, ApJ, 690, 526 [NASA ADS] [CrossRef] [Google Scholar]
  82. Tominaga, N., Maeda, K., Umeda, H., et al. 2007, ApJ, 657, L77 [NASA ADS] [CrossRef] [Google Scholar]
  83. Travaglio, C., Gallino, R., Rauscher, T., Röpke, F. K., & Hillebrandt, W. 2015, ApJ, 799, 54 [Google Scholar]
  84. Travaglio, C., Rauscher, T., Heger, A., Pignatari, M., & West, C. 2018, ApJ, 854, 18 [NASA ADS] [CrossRef] [Google Scholar]
  85. Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Wanajo, S., Sekiguchi, Y., Nishimura, N., et al. 2014, ApJ, 789, L39 [Google Scholar]
  87. Winteler, C., Käppeli, R., Perego, A., et al. 2012, ApJ, 750, L22 [Google Scholar]
  88. Woosley, S. E., & Howard, W. M. 1978, ApJS, 36, 285 [NASA ADS] [CrossRef] [Google Scholar]
  89. Woosley, S. E., & Weaver, T. A. 1995, ApJS, 101, 181 [Google Scholar]
  90. Xu, Y., Goriely, S., Jorissen, A., Chen, G., & Arnould, M. 2013, A&A, 549, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  91. Zahn, J.-P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]

All Tables

Table 1.

Characteristics of our 25 M, Z = 10−3 stellar models at the end of their evolution.

Table 2.

Quantities required to estimate the (O/p) ratio corresponding to the oxygen to p-nuclide yield ratio normalized to the solar abundances (Eq. (6)).

All Figures

thumbnail Fig. 1.

Pre-supernova abundance profiles for the M25Zm3V0 (non-rotating, left panel) and M25Zm3V4 (rotating, υini/υcrit = 0.4, right panel) models computed with the rate of 17O(α,γ)21Ne from Best et al. (2013). The abundances of trans-iron elements are not explicitly followed from core O-burning ignition (see text for details). The dashed grey line shows the pre-supernova temperature profile, the scale of which is indicated on the right axis. The shaded grey area represents the remnant. The cyan area represents the zone in which the maximum temperature during a spherical 1051 erg explosion ranges between 1.8 and 3.7 GK (i.e. the zone in which the explosive p-process takes place).

In the text
thumbnail Fig. 2.

Ratio of 17O(α,γ)21Ne reaction rates to the rate of Caughlan & Fowler (1988, cf88 label). The red curve (label bg13) is the recommended rate of Best et al. (2013). The green curve (label ta19) is the experimental lower limit from Taggart et al. (2019). The shaded area indicates the approximate range of temperatures in the helium-burning core of massive stars.

In the text
thumbnail Fig. 3.

Pre-supernova mass fractions of elements heavier than iron as a function of the mass number A at a mass coordinate of 4 M for the models computed with the rate of 17O(α,γ)21Ne from Best et al. (2013, left panel) and Taggart et al. (2019, right panel). The four vertical dashed lines show the location of the s-nuclides 56Fe, 88Sr, 138Ba, and 208Pb.

In the text
thumbnail Fig. 4.

Elemental mass fractions of all models normalized by the mass fractions of the reference M25Zm3V0 model. Like in Fig. 3, the abundances are shown at a mass coordinate of 4 M. The four vertical dashed lines show the location of Fe (Z = 26), Sr (Z = 38), Ba (Z = 56), and Pb (Z = 82).

In the text
thumbnail Fig. 5.

Temperature profiles at different stages of the evolution (blue lines) and maximum temperature reached during the explosion (magenta line) for the M25Zm3V4 model. The abundance profile of the s-nuclide 88Sr is shown in black with the scale on the right. The blue and magenta rectangle show where the p-process can take place (in terms of temperature range and mass coordinate range) during the hydrostatic and explosive burning, respectively (see text for details). The dark grey area shows the extent of the remnant (corresponding to the size of the 56Ni core; cf. Sect. 2.2). The three horizontal dashed lines correspond to temperatures of 1.5, 1.8, and 3.7 GK.

In the text
thumbnail Fig. 6.

Abundance profiles of the s-nuclide 88Sr (black) and the p-nuclide 74Se (red) in the region of the star that is relevant for the p-process. The abundance profiles are shown at four different stages: at the birth of the star (or zero-age main-sequence, dotted lines), at the start of core oxygen burning (dash-dot lines, at this point, the 74Se abundance is about zero, hence not visible on the plot), at the end of the hydrostatic post-processing calculation (dashed lines) and after the supernova explosion (solid lines).

In the text
thumbnail Fig. 7.

Overproduction factors ⟨Fi⟩ (Eq. (1)) of p-nuclei mass-averaged over the 4.25 M PPLs (corresponding to the union of the blue and magenta boxes in Fig. 5) in the M25Zm3V4 model at the end of the hydrostatic burning pre-supernova phase (green diamonds), only due to the explosive burning (blue) or including both the hydrostatic and explosive burning phases (red; see text for more details).

In the text
thumbnail Fig. 8.

Temperature (left panel) and density (right panel) evolution of several tracer particles coming from the M25Zm3V4 explosion model. The maximum temperature of these particles during the explosion verifies 1.8 <  Tex, max <  3.7 GK. These particles are initially located in the cyan area in Fig. 1 (right panel).

In the text
thumbnail Fig. 9.

Overproduction factors ⟨Fi⟩ (Eq. (1)) mass-averaged over all PPLs for all p-nuclides when considering the 17O(α,γ)21Ne rate of Best et al. (2013, top panel), Best et al. (2013) divided by 10 (middle), and Taggart et al. (2019, bottom panel). The filled symbols highlight the npp-nuclides (Table 2) that have the highest overabundances, i.e. those with ⟨Fi⟩ higher than the highest value divided by an arbitrary factor of 20. The green curves show the velocity-averaged yields over the three 25 M models assuming they follow the rotational distribution of observed young B stars from Huang et al. (2010, their Fig. 6). See text for more details.

In the text
thumbnail Fig. 10.

Overproduction factors ⟨Fi⟩ (Eq. (1)) mass-averaged over all PPLs for all stable nuclei (open symbols) and p-nuclides (full symbols) for the 25 M models with the 17O(α,γ)21Ne rate from Best et al. (2013) and with υini/υcrit = 0 (black diamonds), 0.4 (blue squares) and 0.7 (red circles). The overproduction factors of some α-elements, including 16O, are highlighted.

In the text
thumbnail Fig. 11.

Mass coordinate as a function of explosion energy for the M25Zm3V4 model. The color code indicates the mass fraction of the s-nuclide 88Sr prior to the explosion. The two black lines delimit the zone in which the maximum temperature during the explosion verifies 1.8 <  Tex, max <  3.7 GK. The dark grey rectangular zone shows the extent of the remnant. The dashed line shows the mass coordinate below which some photodisentegration of trans-iron seeds can have occurred during the last hydrostatic burning stages (zone in which T >  1.5 GK at some point during hydrostatic burning).

In the text
thumbnail Fig. 12.

Overproduction factors ⟨Fi⟩ (Eq. (1)) mass-averaged over all PPLs for all p-nuclides for the M25Zm3V4 model with two values of the total explosion energy Etot = 1051 and 1052 erg.

In the text
thumbnail Fig. 13.

Same as Fig. 9, where only the overproduction factors for the three velocity-averaged models are compared.

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.