Open Access
Issue
A&A
Volume 667, November 2022
Article Number A155
Number of page(s) 18
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202244360
Published online 22 November 2022

© A. Choplin et al. 2022

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

Exploring the origins of elements heavier than Fe is among the most topical questions in modern astrophysics (e.g. Arnould & Goriely 2020). We do know that most of them come from neutron capture processes. Furthermore, besides the slow (s) and rapid (r) neutron capture processes, an intermediate (i) neutron capture process (first suggested by Cowan & Rose 1977) is expected to operate at neutron densities intermediate between the s- and r-processes.

Different observations support the existence of the i-process. The so-called carbon-enhanced metal-poor (CEMP)-r/s stars, whose chemical composition is difficult to reconcile with s- or r-process models, point towards the existence of an i-process (e.g. Jonsell et al. 2006; Lugaro et al. 2012; Dardelet et al. 2014; Roederer et al. 2016; Karinkuzhi et al. 2021; Goswami & Goswami 2022). Sometimes, the chemical composition is best reproduced if considering a combination of an s- and an i-process (Koch et al. 2019). Another clue suggesting the existence of the i-process is the puzzling Ba overabundance in open cluster stars (Mishenina et al. 2015). Several works have also shown that pre-solar grains may bear the isotopic signature of i-process nucleosynthesis (Fujiya et al. 2013; Jadhav et al. 2013; Liu et al. 2014).

The i-process is triggered when some hydrogen is mixed in a convective helium-burning zone. This mixing process is either referred to as proton ingestion event (PIE, e.g. Cristallo et al. 2009), flash-driven mixing (e.g. Lau et al. 2009), or dual shell flash (e.g. Campbell & Lattanzio 2008). In this work, we use the term PIE. When a PIE arises, protons are transported down in the convective helium-burning zone and burn on the way via the 12C(p, γ)13N reaction. The decay of 13N to 13C (in about 10 min) is followed by 13C(α, n)16O, operating mostly at the bottom of the convective helium-burning zone, where the temperature is high (> 2 × 108 K). Under these conditions, neutron densities as high as 1012 − 1016 cm−3 are reached for typically ∼1 yr, triggering an i-process nucleosynthesis. Because 12C(p, γ)13N releases a lot of energy in the middle of the convective helium-burning zone, a temperature inversion develops in the 1D models and the convective zone splits.

Proton ingestion events have been studied parametrically using one-zone models, for instance, Hampel et al. (2016, 2019), where the ability of the i-process to reproduce the chemical composition of CEMP-r/s stars is explained. Nevertheless, the astrophysical site(s) hosting PIEs, hence the i-process, remain(s) actively debated. Shortly after Schwarzschild & Härm (1965) first noticed the existence of thermal pulses in asymptotic giant branch (AGB) stars, it has been suggested that PIEs could appear during thermal pulses of these stars (e.g. Schwarzschild & Härm 1967; Despain & Scalo 1976; Fujimoto 1977; Scalo 1979; Fujimoto et al. 1984), despite the fact that the entropy barrier between the H- and He-rich layers tends to prevent such events from occurring (e.g. Iben 1976). Since then, various astrophysical sites hosting PIEs (i.e. the i-process) were suggested: the early AGB phase of metal-poor low-mass stars (Cassisi et al. 1996; Fujimoto et al. 2000; Chieffi et al. 2001; Siess et al. 2002; Iwamoto et al. 2004; Cristallo et al. 2009, 2016; Lau et al. 2009; Suda & Fujimoto 2010; Stancliffe et al. 2011; Gil-Pons et al. 2018; Choplin et al. 2021; Goriely et al. 2021), the core helium flash of very low-metallicity low-mass stars (Fujimoto et al. 1990, 2000; Schlattl et al. 2001; Suda & Fujimoto 2010; Campbell et al. 2010; Cruz et al. 2013), the very late thermal pulses of post-AGB stars (Herwig et al. 2011), rapidly accreting carbon-oxygen (C-O) or oxygen-neon (O-Ne) white dwarfs (RAWDs) in close binary systems (Denissenkov et al. 2017, 2019, 2021; Côté et al. 2018; Stephens et al. 2021), super-AGB stars (7 M ≲ Mini ≲ 10 M, Siess 2007; Jones et al. 2016), or the helium shell of very low- or zero-metallicity massive stars (Mini >  10 M, Banerjee et al. 2018; Clarkson et al. 2018; Clarkson & Herwig 2021).

In many of these studies, the i-process nucleosynthesis accompanying a PIE is not investigated and when it is, often just one model is considered. Indeed, such models can be difficult to compute and require large nuclear reaction networks. So far, only Denissenkov et al. (2019) have produced a grid of seven RAWD models with i-process nucleosynthesis calculated in post-processing. There is a clear need of computing grids of i-process models from the various possible sites.

In Choplin et al. (2021, hereafter Paper I) we studied in details the development of the i-process during the AGB phase of a a 1 M model at a metallicity of [Fe/H] = −2.5. In Goriely et al. (2021, hereafter Paper II) we focused on the nuclear physics uncertainties affecting the i-process during the AGB phase. In this third paper, we investigate the evolution and nucleosynthesis of a grid of 12 AGB models computed with various initial masses and metallicities. We highlight and explain the different evolutions, discuss the different nucleosynthesis processes at work, with a special emphasis on the i-process.

Section 2 presents the input physics. Section 3 discusses structure and evolution aspects, while Sect. 4 focuses on nucleosynthesis, surface enrichment, and stellar yields. Our summary and conclusions are given in Sect. 5.

2. Physical ingredients, numerical aspects, and important quantities

The models presented in this paper are computed with the stellar evolution code STAREVOL (Siess et al. 2000; Siess 2006; Goriely & Siess 2018, and references therein). We computed models with initial masses of 1, 2, 3 M. For each mass, we investigated the metallicities of [Fe/H] = −3, −2.5 and −2, corresponding to Z = 1.4 × 10−5, 4.3 × 10−5, and 1.4 × 10−4 in mass fractions, respectively. We also computed a model of 1 M at [Fe/H] = −2.3 (Z = 6.8 × 10−5). These ten models were computed with the solar mixture of Asplund et al. (2009). To investigate the effect of an α-enhancement, we computed a 1 M and a 2 M model at [Fe/H] = −3 with an α-enhanced mixture. We raised the initial abundance of the isotopes 12C, 16O, 20Ne, 24Mg, 28Si, and 32S following Ritter et al. (2018). Their Table 1 reports the mass fractions of α-enhanced isotopes according to halo and disc stars (Reddy et al. 2006, and references therein). For the elements Ne and S, they refer to Kobayashi et al. (2006). In the end, our α-enhanced models have a metallicity Z = 6.4 × 10−5 in mass fraction (i.e. similar to the [Fe/H] = −2.3 model).

Other input physics are the same as in Papers I and II. Especially, we use the mass-loss rate from Reimers (1975) from the main sequence up to the beginning of the AGB and then switch to the Vassiliadis & Wood (1993) prescription. The opacity change due to the formation of molecules when the star becomes carbon rich is taken into account (Marigo 2002). The mixing length parameter α is set to 1.75 and we do not consider extra mixing (e.g. overshoot, thermohaline). Table 1 summarizes the main characteristics of the models computed.

Table 1.

Main characteristics of the models computed in this work, with the initial mass Mini, [Fe/H] ratio, initial metallicity in mass fraction Z, initial chemical composition (solar-scaled or α-enhanced), total lifetime τstar, duration of the AGB phase τAGB, final mass Mfin, and final mass of the convective envelope , number of pulses Npulse, along with information on whether the model experiences a PIE or not and the pulse number during which the PIE arises.

2.1. The coupling of diffusion and nucleosynthesis

Particular attention is paid to the equation describing the abundances change of chemicals. The abundance Xi of a nucleus i is followed by the equation:

(1)

where D is the diffusion coefficient associated with convection only, since other mixing mechanisms are not included here.

The first and second terms in the right hand side of Eq. (1) account for the changes due to transport and nuclear burning, respectively. In general these two terms are treated separately, one after the other. However, during a PIE, the mixing timescale becomes comparable to the nuclear timescale, so that dealing with the transport and the nuclear burning independently may no longer be valid. Indeed, a chemical (especially protons) may be transported in a zone of the star where it should not exist since it would be burnt before reaching this zone. For this reason, it is necessary to fully couple the transport and nuclear burning. In this series of papers, during the PIE, the nucleosynthesis and transport equations are solved simultaneously once the structure has converged. This method requires many inversions of the Jacobian matrix of size (K × N) × (K × N), where K is the number of nuclei and N the number of shells. In our cases, it corresponds to matrices of about (2 × 106) × (2 × 106), which leads to a total of more than 1012 matrix elements. Fortunately, these matrices are very sparse. The fill factor can be defined as:

(2)

where Nnon_zero is the number of non zero value in the matrix and (K × N)2 the total number of elements. In our models, η ∼ 2 × 10−6. The sparsity of the matrix makes it interesting to use the compressed sparse row format, which reduces the matrix information into two one-dimensional (1D) arrays of size Nnon_zero and one 1D array of size K × N + 1. To solve the system of abundance equations (Eq. (1)), we use the PARDISO solver1 (Schenk & Gärtner 2004).

2.2. Convection

In a convective zone and in the present work, the coefficient D in Eq. (1) is equal to the convective diffusion coefficient Dconv. This coefficient is computed according to the mixing-length theory and can be expressed as:

(3)

where αmlt is the mixing length parameter, HP the pressure scale height, and vconv the convective velocity. The convective turnover timescale τconv is defined as:

(4)

with r1 and r2 the boundaries of the convective region. The turnover timescale expresses the approximate time it takes for a cell of matter to travel across the entire convection zone.

2.3. Nuclear reaction network

The STAREVOL stellar evolution code includes two reaction networks made of either 411 or 1160 nuclei. The largest network includes all species with a half-life greater than about 1 s and the 2123 nuclear reactions (n-, p-, α-captures and α-decays), weak (electron captures, β-decays), and electromagnetic interactions of relevance to properly follow neutron capture processes up to neutron densities of ∼1017 cm−3. In our models, as soon as the neutron density exceeds 1013 cm−3, the large network is adopted.

Nuclear reaction rates are taken from the Nuclear Astrophysics Library of the Université Libre de Bruxelles2 (Arnould & Goriely 2006). The latest experimental rates are extracted through the NETGEN interface tool (Xu et al. 2013). When not available experimentally, the Maxwellian-averaged cross sections are calculated with the TALYS reaction code (Koning & Rochman 2012; Goriely et al. 2008). Additional details can be found in Papers I and II.

2.4. Neutron exposure

Besides the neutron density, a useful indicator to quantify the efficiency of a neutron capture process is the neutron exposure. We define below some quantities that are relevant for the present work. The neutron exposure τ at a given mass coordinate Mr, between times t1 and t2 is defined as

(5)

where Nn is the neutron density and the neutron thermal velocity with kB the Boltzmann constant, T(t) the temperature at time, t, and mn the neutron mass. Since the neutron density Nn varies with time and mass coordinate in the model, it is also interesting to define the neutron exposure averaged over the convective pulse and over the duration of the PIE:

(6)

where the neutron density times thermal velocity averaged over the mass of the convective pulse, at time, t, is defined as:

(7)

with Mr, 1(t) and Mr, 2(t) the mass coordinate boundaries of the convective pulse at time, t, Nn(Mr), and vT(Mr) the neutron density and thermal velocity at mass coordinate, Mr, respectively. We also define τbot, which is the neutron exposure at the bottom of the convective pulse, where the neutron density is the highest. It is expressed as Eq. (5) but Mr can vary so as to follow the bottom of the convective pulse. In this case, the first temporal boundary t1 in Eq. (5) corresponds to the beginning of the PIE. The second temporal boundary t2 is chosen to be either just before the split, if any, or at the end of the PIE, when the maximal neutron density in the star has dropped below ∼1010 cm−3. These quantities are given in Table 2.

Table 2.

Important characteristics of our models experiencing a PIE.

2.5. Yields and overproduction factors

The yield 𝒴i of a nucleus i is computed according to the relation:

(8)

where τstar is total lifetime of the model star (given in Table 1), and Xi, s(t) and (t) are the surface mass fraction of nucleus, i, and the mass-loss rate at time, t, respectively.

We also use the overproduction factor, fi, associated with an element, i, defined as:

(9)

with Mini and Mfin the initial and final mass of the model, respectively, and Xi, ini the initial mass fraction of element, i.

2.6. Mass of hydrogen ingested during a PIE

The mass of hydrogen ingested during a PIE is critical. It can have dramatic consequences on the evolution, structure, and nucleosynthesis. In our models, protons are naturally ingested, so this mass is not a free parameter.

We define which is the mass of hydrogen engulfed from the start of the PIE to the point where the neutron density is maximum. It can be written as:

(10)

where MH is the total mass of hydrogen contained in the star, evaluated just before the start of the PIE at t = t0 and when the neutron density reaches its maximum at t = tNnmax. We typically have tNnmax − t0 ∼ 1 yr. Similarly, we define , which is the hydrogen mass difference between the start of the PIE and the time of splitting of the convective pulse at t = tsplit. We checked that the mass of hydrogen lost through winds during these intervals of time is negligible compared to the mass of the ingested hydrogen. We also considered that the H-burning shell is switched off just before the PIE, which is a reasonable assumption given that most of the energy comes from He-burning during a thermal pulse. We note that neither nor correspond to the total mass of hydrogen ingested, since more hydrogen is engulfed after the splitting of the convective pulse, when the upper part grows in mass before merging with the convective envelope.

3. Structure and evolution aspects

We routinely checked during the AGB phase of our models whether a PIE occurs by recomputing a given pulse with higher spatial and temporal resolutions. As mentioned in Paper I, we noticed in our models that a PIE can be missed if adopting a too large time step.

We can see from Table 1 that 6 out of our 12 models experience a PIE: the 1 M models with [Fe/H] = −3.0, −2.5 and −2.3, the 1 M model with α-enhancement and the 2 M models with [Fe/H] = −3.0 and −2.5. Such events are more likely to take place at lower initial masses and metallicities, as already noticed in previous works (e.g. Iwamoto et al. 2004; Campbell & Lattanzio 2008; Suda & Fujimoto 2010). Also, PIEs always take place during the first or second thermal pulses3 (Table 1). In the models with a PIE, two different categories emerge: (i) the models in which the TP-AGB phase ends after the PIE because of strong mass loss or (ii) the models in which the TP-AGB phase resumes after the PIE. The M1.0z2.5, M1.0z3.0α and M1.0z2.3 models belong to the first category while the M1.0z3.0, M2.0z3.0 and M2.0z2.5 to the second category. The physical reasoning behind to these different evolutionary pathways are discussed in the next sections.

3.1. The AGB phase of the 1 M, [Fe/H] = −3.0 model

The M1.0z3.0 model experiences a PIE during the very first thermal pulse and a normal AGB phase afterwards with 41 thermal pulses without any further PIE. We discuss below the peculiar structure and evolutionary aspects of this model during and after the PIE.

3.1.1. Evolution during the PIE

Once the PIE starts (Fig. 1 top panel, at t − t0 ∼ 0.4 yr), protons are transported down in the convective thermal pulse (the turnover timescale is about 1 h) and burnt on the fly by 12C(p, γ)13N. With a half-life of about 10 min, 13N decays into 13C which is transported at the bottom of the pulse where T >  200 MK (Fig. 1, second panel). At this temperature, 13C(α, n)16O is efficiently activated (Fig. 1, third panel) and releases neutrons. It gives rise to a maximum neutron density of about 1015 cm−3 (Fig. 1, fourth panel). Contrary to most of our models, this model does not experience a proper split of the pulse, as explained in Sect. 3.5. The energy released by nuclear burning produces an outward expansion of the convective pulse which engulfs more protons, as seen in the top panel of Fig. 1, for t − t0 ≳ 0.5 yr. At the same time, the pulse dilates (see the iso-mass contours in top panel of Fig. 1), leading to a decrease in the temperature at the bottom of the convective region. This reduces the efficiency of 13C(α,n)16O (blue line in Fig. 1, third panel) and the neutron production. Because a copious amount of protons has been engulfed, reactions from the CNO cycle become dominant (black, red, and green line in Fig. 1, third panel).

thumbnail Fig. 1.

Evolution of the M1.0z3.0 model during the PIE: Kippenhahn diagram (top panel), temperature at the bottom of the convective thermal pulse (or envelope for t − t0 ≳ 1 yr, second panel), energy released by the main reactions at the bottom of the pulse (or envelope, third panel), maximal neutron density (fourth panel). Bottom panel: 12C/13C, C/N, and O/N ratios at the bottom of the convective zone (solid lines) and at the stellar surface (dashed lines). The dashed green and blue lines in the Kippenhahn diagram (top panel) delineate the hydrogen and helium-burning zones (where the nuclear energy production by H- and He-burning exceeds 10 erg g−1 s−1). The thin light-grey lines indicate iso-masses expressed in M.

Reactions from the CNO cycle starts to be active in the thermal pulse as soon as protons are ingested. At first, they operate together with He-burning reactions when the temperature is high enough. In a second step they operate without He-burning when the temperature drops below about 150 MK. The operation of the CNO cycle in the pulse dramatically reduces the C/N, O/N, and 12C/13C ratios (Fig. 1, bottom panel). At the end of the PIE, at t − t0 ∼ 0.7 yr, the 12C/13C ratio at the bottom of the pulse is 3.4 which corresponds to the CNO equilibrium value. At t − t0 ∼ 1 yr, the pulse eventually merges with the convective envelope. The surface ratios (Fig. 1, bottom panel, dashed lines) are progressively modified by the large amount of metals coming from the pulse.

3.1.2. Evolution after the PIE

It’s reasonable to wonder why this model does not enter a standard thermally pulsating AGB phase (TP-AGB), while the M1.0z2.5 model does (cf. Paper I for more details on the M1.0z2.5 model). As explained below, the difference between the M1.0z2.5 and M1.0z3.0 models is due to the fact that the M1.0z3.0 PIE occurs in the very first thermal pulse, but not until the second pulse for the M1.0z2.5 model (Table 1).

During the early AGB phase, the maximal temperature at the bottom of the pulse quickly increases with the pulse number. The maximal temperature at the bottom of the pulse during the PIE reaches 213 MK and 257 MK for the M1.0z3.0 and M1.0z2.5 models, respectively. The lower temperature in the M1.0z3.0 models leads to slower nuclear burning, especially the 3α reaction. The amount of metals synthesized in the M1.0z3.0 model is consequently lower4 and the convective envelope is therefore less enriched in metals after the PIE. Just after the PIE, the surface metallicity in mass fraction rises to 7.5 × 10−3 (1.9 × 10−2) in the M1.0z3.0 (M1.0z2.5) model. In particular, the surface 12C and 16O mass fractions after the PIE are 0.7 × 10−3 and 1.6 × 10−3 for the M1.0z3.0 model while they are 7.7 × 10−3 and 4.9 × 10−3 for the M1.0z2.5 model. This leads to surface C/O ratios in number after the PIE of 0.75 and 2.55 for the for M1.0z3.0 and M1.0z2.5 models, respectively. This results in a less dramatic increase in CO molecular opacities in the stellar envelope of the M1.0z3.0 model, which, in turn, implies a weaker mass loss (more details in Sect. 3.3). In the end, because of the lower surface metallicity and C/O ratio, the M1.0z3.0 model does not lose its envelope as fast as the M1.0z2.5 model and the AGB phase consequently resumes after the PIE. We also note that the PIE starts 220 yr (305 yr) after the beginning of the thermal pulse in the M1.0z3.0 (M1.0z2.5) model. It means that helium-burning has less time to process in the M1.0z3.0 model, which also contributes to the smaller C and O production and surface enrichment in this star. The temperature at the bottom of the pulse and the interval of time between the start of the pulse and the PIE are two critical aspects that dictate the subsequent AGB evolution.

We confirmed our arguments by doing a numerical experiment on our M1.0z3.0 model. We intentionally missed the PIE happening in the first pulse by setting up time steps that were too large. Then, at the second pulse, a PIE happens. This time, the maximal temperature in the pulse is higher (around 250 MK), so that the model follows the same behaviour as the M1.0z2.5 model: high surface enrichment, strong mass loss, and the end to the AGB phase.

Ultimately, we see that PIEs tend to happen earlier in the AGB phase (during the very first thermal pulse) of lower metallicity models. If this happens, an AGB phase with (many) thermal pulses can follow the PIE, in contrast to higher metallicity models where the AGB phase ends right after the PIE.

We may also wonder why no other PIE happen in the M1.0z3.0 model after the first one. As mentioned in Iwamoto et al. (2004), once a PIE has taken place, it becomes more difficult to have a second one. Indeed, following a PIE, the envelope is enriched in metals which results in a larger entropy barrier between the He- and H-rich layers and, thus, this tends to prevent any further PIE. After a PIE, the situation resembles that of a star of higher metallicity, where PIEs are less prone to arise.

3.2. The AGB phase of the 2 M, [Fe/H] = −2.5, and −3.0 models

The M2.0z3.0 and M2.0z2.5 models experience a PIE followed by about 10 standard thermal pulses (Table 1). They follow a similar evolutionary pathway as the M1.0z3.0 model but the physical reasons leading to it are different.

The maximal temperature at the bottom of the PIE pulse is similar or even higher in the M2.0z3.0 and M2.0z2.5 models than in the M1.0z2.5 model. In the pulse, carbon and oxygen are produced in similar amounts than in the M1.0z2.5 model. An important difference however is higher dilution factors in the 2 M stars compared to the 1 M stars. The ratios between the pulse and envelope masses are about 0.1 and 0.01 for the 1 and 2 M, models, respectively5. The material from the pulse in the 2 M, models is therefore about ten times more diluted in the envelope than in the 1 M and the surface metallicity does not increase that much in this model. In the M1.0z2.5 model, the surface metallicity rises to 1.9 × 10−2 while it increases only to 1.2 × 10−3 and 2.6 × 10−3 in the M2.0z3.0 and M2.0z2.5 models, respectively. The high surface enrichment in the M1.0z2.5 model leads to high surface opacities and a large expansion of the stellar radius from ∼80 to ∼400 R. The effective temperature drops to 2700 K. In the M2.0z2.5 and M2.0z3.0 models, the radius ranges from ∼120 to ∼200 R and the effective temperature experiences a milder drop, down to about 4000 K. In the end, the mass loss rate after the PIE does not increase dramatically in the M2.0z2.5 and M2.0z3.0 models (unlike the M1.0z2.5 model), so the AGB phase resumes with additional thermal pulses.

3.3. The key role of molecular opacities

As mentioned in Sect. 2, molecular CO opacities are taken into account in our models. Molecular CO opacities become important when the effective temperature drops below 5000 K and when the C/O ratio becomes larger than about 0.8 in number. These opacities play a key role in the models that lose all their envelope quickly after the PIE (e.g. in the M1.0z2.5 model). The convective envelope of these models becomes highly enriched in C, N, and O after the PIE with a C/O ratio greater than 1. This leads to high CO molecular opacities and enhanced mass loss which prevents the further occurrence of thermal pulses.

We confirmed the critical role of CO molecular opacities by recomputing our M1.0z2.5 model starting from just after the PIE with the CO molecular opacities switched off. In this case, the surface opacity (hence, the stellar radius) does not rise dramatically, the mass loss rate modestly increases, and the star keeps its envelope for a much longer time. The PIE is consequently followed by a normal AGB phase with regular thermal pulses (as in the M1.0z3.0 model, see Sect. 3.1).

3.4. The entropy barrier

To our knowledge, there is no simple criterion at hand to predict whether a model will experience a PIE or not. However, as noticed in previous works (e.g. Fujimoto et al. 1990; Iwamoto et al. 2004), the height of the entropy barrier at the bottom of the H-rich layers plays an important role in the occurrence of a PIE. The higher the barrier, the less likely the convective pulse is to reach the H-rich layers, hence, the less likely it is for a PIE to occur.

In the present work, we estimate this barrier at the maximal development of each thermal pulse for each model. We define the entropy barrier as ds = s2 − s1 where s1 is the entropy at the top of the convective pulse and s2 is the entropy at the base of the H-rich zone (defined where the mass fraction of hydrogen drops below 10−1). For models experiencing a PIE, the entropy barrier was derived just before the PIE. The results are reported in Fig. 2.

thumbnail Fig. 2.

Entropy barrier between the top of the convective thermal pulse and the bottom of the H-rich layers as a function of the pulse number (see text for details on the computation of ds). The empty circles mark the occurrence of the PIE.

The entropy barrier overall increases with the pulse number. This is consistent with the fact that PIE preferentially happen during the first thermal pulses (Table 1). The entropy barrier tends to be higher in higher mass models which makes a PIE less likely to occur as mass increases (as recognized early by e.g. Iben 1977, especially his Fig. 3). This is clearly visible, at least, in our 3 M, models, as compared to the less massive models, especially during the first pulses. Higher entropy barriers in more massive AGB stars is consistent with the fact that PIEs take place in some of our 1 and 2 M, models, but not in our 3 M, models.

It was also recognized that PIEs preferentially occur at low metallicities because of the smaller entropy barrier at the base of the hydrogen-burning shell (e.g. Fujimoto et al. 1990). The metallicity dependence can be understood by the fact that less CNO catalysts are present in the H-burning shell of low metallicity AGB stars. The entropy barrier between the H- and He-burning zones consequently decreases with metallicity, thus facilitating PIEs. In our models, as soon as ds/NAkB ≲ 2.5, a PIE is triggered (Fig. 2). Nevertheless, the metallicity dependence is not clear as models with the same initial mass have overall similar entropy barriers (Fig. 2). To explain this, we first note that although the height of the entropy barrier gives a good hint on the possible occurrence of a PIE, it should be taken with caution because the definition of this barrier is not absolute and can be estimated using different criteria. Secondly, our models lie in the metallicity range separating models with and without PIE. Figure 3 shows a compilation of AGB models from different works investigating PIEs in AGB stars. All the models were computed without extra mixing processes except the model from Cristallo et al. (2009, 2016) that includes overshooting below the convective envelope following the time dependent formalism described in Chieffi et al. (2001). Overall, PIEs occur in all the models shown in the dark grey region, while the light grey region is subject to some uncertainties. Our models (especially the 1 and 2 M, models) lie close to the outer grey border, where very slight stellar structural differences can be responsible for the presence or absence of PIEs.

thumbnail Fig. 3.

Mass-metallicity diagram showing the occurrence of PIEs during the early AGB phase of models from various authors. Filled symbols show models experiencing a PIE while empty symbols are for models that do not experience a PIE. The dark grey zone shows the approximate region where PIEs happen in all of the models and the light grey one where PIEs happen in most of the models. The corresponding [Fe/H] ratios are indicated on the right axis assuming solar-scaled mixtures. Models are from Iwamoto et al. (2004, red triangles), Campbell & Lattanzio (2008, magenta circles), Cristallo et al. (2009, blue triangle), Lau et al. (2009, blue diamonds), Suda & Fujimoto (2010, orange triangles), Cristallo et al. (2016, black pentagons). Model results from this work are shown as green squares. All models were computed without extra mixing processes, except the models from Cristallo et al. (2009, 2016), which consider overshooting below the convective envelope.

3.5. The split of the convective pulse during a PIE

In 1D simulations, during a PIE, proton are transported down in the convective helium-burning zone until the timescale associated with the reaction 12C(p, γ)13N becomes similar to the local transport timescale of protons. At this location, protons accumulate, the nuclear energy production from 12C(p, γ)13N is maximal, a temperature inversion develops, and the helium-burning convection zone splits. From this point, protons cannot be mixed in the hottest bottom part of the pulse anymore and the chemicals synthesized in the lower part of the now-split convective region (especially through the i-process) remain locked in the star unless further thermal pulses and dredge-up events take place later on. The split is thus an important aspect of a PIE since it can prevent the enrichment of the AGB surface in heavy elements if it happens before neutrons are released.

In our M2.0z3.0 model the energy released from the main reactions just before (top panel) and after (bottom panel) the splitting is shown in Fig. 4. We clearly see the two peaks of nuclear energy at the bottom of the pulse and at Mr ≃ 0.642 M. The first peak is mainly due to the 3α reaction while the second one is attributed to 12C(p, γ)13N. A fraction of a day later, the pulse has split and the bottom part becomes much less active compared to the upper part, driven by hot hydrogen-burning. In this model, the split takes place after the neutron density has reached its maximum, so that the upper part of the pulse is fully enriched in i-process products before the splitting (Sect. 4.2 for a more detailed discussion).

thumbnail Fig. 4.

Energy generation from the main reactions during the PIE in the M2.0z3.0 model just before the splitting (top panel), and just after the splitting (bottom panel). The total nuclear energy production rate is shown by the thick grey line. The protons and neutrons mass fractions are shown.

All our models with a successful PIE (except for M1.0z3.0) follow a similar behaviour, although the split takes place at different depths in the pulse. Our M1.0z3.0 model experiences an early splitting of the convective pulse but both parts of the pulse merge again after a few days, as shown in Fig. 5. Also, contrary to the other models, this early splitting occurs before the neutron density reaches its maximum. Below we develop the reasoning behind this peculiar behaviour.

thumbnail Fig. 5.

Energy generation from the main reactions during the PIE in the M1.0z3.0 model just before the splitting (top panel), just after the splitting (middle panel) and just after both parts of the pulse have merged again (bottom panel). The total nuclear energy production rate is shown by the thick grey line. The protons and neutrons mass fractions are shown.

Firstly, as discussed in Sect. 3.1.2, the M1.0z3.0 model experiences a PIE during the very first thermal pulse, contrary to all other models (cf. Table 1). The temperature in the first pulse is lower than in subsequent pulses. In such conditions, protons can be transported deeper in the convective pulse before being completely burnt by 12C(p, γ)13N. This makes the two energy peaks closer in the M1.0z3.0 model compared to the M2.0z3.0 model (top panels of Figs. 4 and 5) and does not favor the split. The reaction 13N(β+)13C (pink) is mostly responsible for the smaller contrast between both peaks (this point is discussed in detail in Sect. 4.3).

Secondly, the M1.0z3.0 model is the model that ingests the smallest mass of protons (2 × 10−7 M, see Table 2). This makes a split less prone to happen. Hydrogen-burning in the pulse is consequently weaker in the M1.0z3.0 model while helium-burning operates similarly in all models (Figs. 4 and 5, during the split). This leads to a much lower entropy barrier in the M1.0z3.0 model: the entropy barrier at Mr = 0.515 M is s/NAkB = 10−3 while the barrier at Mr = 0.641 M in the M2.0z3.0 model is s/NAkB = 2.3. In the M1.0z3.0 model, the energy from the bottom part of the pulse is enough to make this zone fully convective again and overcomes the small entropy barrier separating both parts of the pulse. It results in the merging of the bottom and upper parts of the pulse (Fig. 5, bottom panel).

In the end, all our models with PIEs experience a split of the convective pulse but in one case (M1.0z3.0), both parts of the pulse quickly merge back again. This model behaves as if no split had developed. This can happen when a small amount of protons are engulfed or if the temperature in the pulse is sufficiently low (or both), so that the protons can be transported close enough to the bottom of the pulse.

As a final remark, we note that other PIE models in the literature are not always accompanied by a splitting. For instance, it is absent in the 1D RAWD models of Denissenkov et al. (2019) or in the 3D AGB simulation of Stancliffe et al. (2011) but reported in the 3D simulation of Herwig et al. (2011).

4. Nucleosynthesis, surface enrichment, and yields

4.1. Pulse splitting and the impact of the convective velocity

As discussed in Sect. 3.5, a split can take place at the location where τC the timescale associated with the reaction 12C(p, γ)13N is similar to the local transport timescale of protons by convection. To better understand this aspect, we recomputed the M1.0z3.0 model starting just before the PIE while arbitrarily changing the convective velocity vconv. This numerical experiment allows for the role of convective velocity on the splitting of the convective pulse to be understood. We divided vconv by a factor of 10, 100, and 1000.

The split occurs at Mr = 0.516 M for the vc/10 model and at Mr = 0.527 M for the vc/100 model (Fig. 6). The standard M1.0z3.0 model does not experience a proper split (as explained in Sect. 3.5) and the vc/1000 model follows a standard evolution. It is clear that the lower the convective velocity the higher in mass the split occurs. This is because a lower vconv leads to a higher τconv and since the split occurs where , if τconv is larger, this location corresponds to a higher mass coordinate, where the temperature is lower and thus the reaction rate of 12C(p, γ) weaker (i.e. τC is smaller).

thumbnail Fig. 6.

Kippenhahn diagram and maximal neutron density of the 1 M model with [Fe/H] = −3.0 during the PIE while considering different values for the convective velocity: the standard case (top left), vconv divided by 10 (top right), 100 (bottom left), and 1000 (bottom right). The dotted green and blue lines in the Kippenhahn diagram correspond to the maximum nuclear energy production by hydrogen- and helium-burning, respectively. The dashed green and blue lines delineate the H- and He-burning zones (when the production of energy by H- and He-burning exceeds 10 erg g−1 s−1).

In the vconv/1000 model, protons do not go deep enough in the pulse. They do not release enough energy through 12C(p, γ)13N to make a clear second energy peak in the pulse (as in Fig. 4, top panel). This prevents the splitting of the pulse. Also, the turnover timescale τconv is about 1 yr, which becomes non-negligible compared to the pulse evolution timescale (10 − 100 yr). The pulse recedes in mass before significant protons could be engulfed and burnt.

The maximal neutron densities obtained in the four cases mentioned above are shown in Fig. 7 and the resulting surface abundances after the PIE are shown in Fig. 8. The lower the convective velocity, the lower the neutron density and the less dramatic the surface enrichment in heavy elements. If vconv is too small, protons do not have time before the split occurs to reach the hot layers of the pulse where 13C(α,n)16O is efficiently activated. Nevertheless, the convective velocity has to be decreased by a factor of more than 10 before seeing any effect on the nucleosynthesis. The vconv/10 model (red pattern in Fig. 8) leads to a similar surface enrichment compared to the standard model (black pattern).

thumbnail Fig. 7.

Maximum neutron density for the M1.0z3.0 model during the PIE for different values for the convective velocity.

thumbnail Fig. 8.

Surface mass fractions after the PIE for the M1.0z3.0 model for different values for the convective velocity.

We also consider a case where vconv multiplied by a factor of 10, according to the results of a 3D simulation by Stancliffe et al. (2011). Increasing the convective velocity beyond the standard value has a small impact (Figs. 7 and 8). The maximal neutron density profile is similar to the standard case, although it reaches slightly higher values and the abundances are lower by a factor of typically 5 − 7 for 60 ≤ Z <  82 (except for Tantalum at Z = 73, where the abundance is lower by a factor of 14).

These results show that the i-process nucleosynthesis and surface enrichment following a PIE is a robust feature that does not depend strongly on the convective velocity. Nevertheless, we are aware that the physics of convection (hence, of PIE) cannot be fully captured by 1D models on the basis of the MLT formalism.

4.2. Surface enrichment after a PIE

Once a PIE has occurred, the important question of whether or not the nucleosynthesis products are brought up to the surface remains. The merging of the convective pulse with the convective envelope will bring some pulse material up to the surface. However, if the pulse splits, the lower part of the pulse after the split will remain locked deep into the star (this material can nevertheless reach the surface during subsequent thermal pulses, if any, followed by third dredge ups). In this case, one major question is whether or not i-process nucleosynthesis happens before the split, both before and after the split, or only after the split. If it happens only before the split for instance, the surface will be fully enriched by i-process products. To address this point, we define the following enrichment ratio, e:

(11)

where tsplit is the time of the split, tNnmax the time corresponding to the maximum neutron density, and ⟨τconv⟩ is the time-averaged convective turnover timescale between tNnmax and tsplit, which can be written as:

(12)

We assume that most of the i-process nucleosynthesis takes place at t = tNnmax. If e ≫ 1, the i-process material has enough time to homogenize in the thermal pulse before the split. In this case, the i-process products are fully transported to the surface. If 0 <  e <  1, the split would occur before the pulse homogenization. A negative e means that the neutron density peak takes place after the splitting. In this case, the i-process material is locked into the lower part of the pulse and will not enrich the surface.

In all our models experiencing a PIE and a split, e ≫ 1. For the M2.0z3.0 model for instance, tsplit − tNnmax = 0.23 yr = 1985 h and the mean turnover timescale in the convective pulse during this interval is 3.95 h. This leads to e = 503. For the M1.0z2.5 model, tsplit − tNnmax = 0.014 yr = 125 h, τconv ≈ 1.70 h leading to6e = 74. The enrichment factors e for our models are reported in Table 2.

4.3. The 13N(n,  p)13C reaction

Because protons mix into an He-burning region, PIEs involve different reactions from both H-burning at very high temperatures and overall He-burning. As can be seen in the top panels of Figs. 4 and 5, the proton abundance decreases with decreasing Mr until a point where it increases again, close to the bottom of the pulse. The protons here are not those that are ingested and transported downwards in the pulse. They come from the 13N(n, p)13C reaction and, to a smaller extent, from 14N(n, p)14C. At 200 − 250 MK, the rate of 13N(n, p)13C is about 500 times higher than the rate of 14N(n, p)14C. The neutrons required for these reactions mainly come from the 13C(α, n)16O reaction. The 13N isotope is synthesized upwards in the pulse by 12C(p, γ)13N and diffuses downwards until it reaches neutron-rich layers. Although 13N decays to 13C via β+ in only 10 min, this timescale is comparable to the time it takes for 13N to be transported from the middle to the bottom of the pulse. This can be seen in Fig. 5 (top and bottom panel) where it is clear that the 13N(β+)13C reaction is active, even at the bottom of the convective pulse. At the end, the production of protons at the bottom of the pulse reactivates the 12C(p, γ)13N reaction. Although the production of protons remains rather small, it is much stronger in the M1.0z3.0 model (Fig. 5, top panel) than in the M2.0z3.0 model (Fig. 4, top panel).

The rate of the 13N(n, p)13C reaction is derived as the reverse rate of 13C(p, n)13N, available in both the compilations of Caughlan & Fowler (1988) and Angulo et al. (1999, NACRE). In this work we use the rate from Angulo et al. (1999).

As a numerical test, we re-computed a PIE model with the rate of 13N(n, p)13C divided or multiplied by a factor of 10. We found that the surface abundances of trans-iron elements after the PIE are impacted by a factor of about 3 at most. Therefore, although this reaction appears to be important for i-process nucleosynthesis, both as a poisoning reaction and as an energy source, its variation by a factor of 10 does not impact the outcome of our models significantly. A more detailed analysis of the existing uncertainties and their subsequent impact on nucleosynthesis is required to draw firmer conclusions.

4.4. i-process nucleosynthesis

We highlight here the main features of the i-process nucleosynthesis in our models (see also Sect. 4.9 for a discussion on the overall i-process yields). For additional details, we refer to Papers I and II, where the i-process nucleosynthesis has already been extensively discussed.

Our models with PIE experience a maximum neutron density ranging between 6.8 × 1013 and 2.2 × 1015 cm−3 (Table 2). The i-process paths for these two extreme models (M2.0z3.0 and M1.0z2.3) are shown in Fig. 9. As expected, the M1.0z2.3 model with a higher neutron density follows a path which is farther away from the valley of β-stability. The i-process in our other models follows similar i-process path, in between the blue and green paths shown in Fig. 9.

thumbnail Fig. 9.

Main i-process path (starting from 56Fe) in the M1.0z2.3 (blue arrows) and M2.0z3.0 model (green arrows), at the bottom of the pulse, at the time of maximum neutron density (Nn = 2.19 × 1015 and 6.76 × 1013 cm−3, respectively). The dashed lines show the secondary paths, where at least 30% of the total flux goes. The four panels corresponds to four different zones in the (N, Z) plane. The black squares highlight the stable nuclei. The abundances of the M1.0z2.3 model are shown by the red colour scale in mass fraction.

Figure 10 shows that in the models experiencing a PIE (bold frames) the surface is enriched right after the PIE (when the pulse merges with the envelope) and is barely modified during the subsequent TP-AGB evolution. The rather weak convective s-process that may take place in the thermal pulses after the PIE (as discussed in Sect. 4.5) does not impact the surface abundances in these models.

thumbnail Fig. 10.

Surface elemental mass fractions normalized by their initial mass fractions for the nine models with Mini = 1, 2 and 3 M, and [Fe/H] = −2, −2.5 and −3. The initial mass is increasing from the top to the bottom panels and the initial metallicity is increasing from the left to the right panels. The patterns correspond to the surface abundances after the indicated pulse number. The bold frames highlights models experiencing a PIE. Unstable isotopes are not beta-decayed and the unstable Z = 43 and Z = 61 elements are not considered since Xini = 0 for those elements.

Our i-process models show, on average, a rising distribution from Fe to Pb. Tantalum (Z = 73) is largely overproduced except in the M2.0z3.0 model which experiences the weakest i-process (cf. Table 2). In Paper I, where only the M1.0z2.5 model was considered, we noticed that elements with Z ≲ 50 were underproduced by ∼0 − 0.5 dex compared to a standard s-process nucleosynthesis, while heavier elements with Z ≳ 50 were overproduced by ∼0 − 0.5 dex. The important production of heavier elements is visible in the M1.0z3.0 and M1.0z2.5 models (Fig. 10), where we notice a jump between I (Z = 53) and Xe (Z = 54). In the 2 M models, this feature is less pronounced. The difference between the elements with Z ≲ 50 and Z ≳ 50 is also smaller.

The progressive increase in Nb (Z = 41) in the M2.0z3.0 model for instance is due to the decay of 93Zr via 93Zr(β)93Nb. The lifetime of 93Zr is 1.61 Myr which is comparable to the AGB lifetime (1.09 Myr for the M2.0z3.0 model, cf. Table 2).

4.5. The convective s-process

After about ten pulses, the temperature at the bottom of the pulse of our 3 M models is high enough (at least 350 MK, Fig. 11, top panel) to efficiently activate the 22Ne(α,n)25Mg reaction. This leads to the convective s-process (e.g. Goriely & Siess 2005; Karakas & Lattanzio 2014). Our 2 M models barely reach the required temperature and therefore only experience a very weak convective s-process, as detailed below.

thumbnail Fig. 11.

Evolution of various quantities as a function of the pulse number. Top panel: maximum temperature reached at the bottom of the pulse. The horizontal dashed line shows the temperature above which the 22Ne(α, n)25Mg reaction becomes efficient. Middle panel: third dredge up efficiency parameter λ (see text for the definition). The λ parameter for the pulses where a PIE occurs is not shown since it cannot be properly defined. Bottom panel: surface 25Mg mass fraction.

The middle panel of Fig. 11 shows the λ parameter which indicates the third dredge-up efficiency. It is defined as the ratio of the mass extent of the envelope in the pulse region to the mass increase in the H-free region during the interpulse (e.g. Karakas et al. 2002). A high value corresponds to an efficient third dredge-up7. The λ parameter is greater than zero in our 2 and 3 M, models, which means that some of the pulse material reaches the stellar surface. In particular, this leads to an increase in 25Mg (Fig. 11, bottom panel) which is a signature of the operation of the 22Ne(α,n)25Mg reaction in the convective pulse. The surface enrichment in 25Mg is smaller for the 2 M, models (red patterns) compared to the 3 M, models (blue patterns) due to the lower pulse temperatures of 325 MK and smaller dredge-up efficiency (λ) in the 2 M, models.

In the 3 M, model, the convective s-process leaves a clear chemical signature at the surface. After about ten pulses, the envelope already shows enrichment in light s-elements (Fig. 10, bottom panels). The dilution of the pulse material in the large convective envelope does not lead to overproduction factor greater than about 50.

The situation is different in the 2 M, models. The M2.0z2.0 model (which does not experience any PIE) shows a very weak signature of the convective s-process at its surface. The overproduction factors does not exceed a factor of 2 − 3 (Fig. 10, middle right panel). The temperature of 325 MK at the bottom pulse in this model is too low to leave a noticeable s-process signature at the surface. The M2.0z3.0 and M2.0z2.5 models experience a PIE during the second pulse which strongly impacts the envelope composition. The subsequent dredge-up events have little effect on the surface abundances which remain almost unaffected (Fig. 10, middle panels).

4.6. Hot hydrogen-burning

Hydrogen-burning through the CNO cycle operates efficiently during the PIE in the thermal pulse as well as at the bottom of the convective envelope for our 3 M, models (hot bottom-burning). We briefly review its impact on the surface abundances below.

4.6.1. In the pulse during a PIE

The degree of CNO processing during a PIE varies from one model to another. In Table 3, we report some surface abundance ratios after the PIE. This pollution results from the interplay between He-burning, H-burning at high temperatures, and dilution with the material in the stellar envelope. All models but M2.0z3.0 show a 12C/13C ratio very close to the CNO equilibrium value of 3 − 4. The model that experiences the highest degree of CNO processing during the PIE is the M1.0z3.0 model (as seen from the low 12C/13C, C/N, and O/N ratios) because of the absence of a proper split in this model (cf. Sect. 3.5). As a result, the temperature at the bottom of the pulse remains high for a longer period allowing for further CNO processing (hence, the strength of H-burning reactions). In the other models, when the split occurs, the temperature at the bottom of the upper part of the pulse strongly decreases, so that H-burning is slowed down. The drop in temperature depends on where the splits occurs. The pulse of the M2.0z2.5 model splits close to the bottom of the pulse. Consequently, the CNO cycle remains efficient in the upper part of the pulse of this model. This results in a stronger CNO cycle signature at the surface (Table 3). In all the models (except M1.0z3.0 and M2.0z2.5), the C/N and O/N ratios are greater than 1, which means that the CNO cycle did not reach equilibrium otherwise 14N would have been the most abundant isotope and both C/N < 1 and O/N < 1. In the end, we found that models with no or deep split are more prone to show a very high degree of CNO processing at their surface.

Table 3.

Surface 12C/13C, C/N, O/N, and C/O ratios in number just after the PIE.

After a PIE, however, the AGB phase may resume and the surface composition can still change. In the M2.0z3.0 and M2.0z2.5 models for instance, after the decrease in the 12C/13C ratio caused by the PIE, the ratio progressively increases again as a result of the enrichment of the envelope with the pulse material (especially 12C), thanks to the operation of the third dredge up (Fig. 12).

thumbnail Fig. 12.

Surface 12C/13C ratio in number as a function of the pulse number. The ratio is shown at the middle of each pulse (i.e. before the third dredge up, if any).

4.6.2. Hot bottom-burning

The CNO cycle is also activated if the temperature at the base of the convective envelope rises above ∼50 MK (e.g. Karakas & Lattanzio 2014). In PIE models, right after the merging of the pulse with the envelope, the temperature at the bottom of the envelope can reach ∼90 MK (Fig. 13), but this lasts for a short period of time (less than 1 yr in our M1.0z3.0 model, Fig. 1), so that it barely affects the envelope composition.

thumbnail Fig. 13.

Temperature at the bottom of the convective envelope during the AGB phase. Models of 1, 2, and 3 M, are shown in the top, middle, and bottom panels, respectively.

Our 3 M, models experience hot bottom-burning, although the temperature at the base of its envelope barely reaches 50 MK (Fig. 13, bottom panel). It nevertheless leaves a clear signature at the stellar surface: a progressive decrease in the 12C/13C ratio (Fig. 12 top panel) and to a lesser extent, of the C/N ratio. The moderate envelope temperature prevents a full activation of the Ne-Na and Mg-Al cycles.

4.7. Effect of α-enhancement

The M1.0z3.0α model was computed with an α-enhanced mixture (cf. Sect. 2). This model does not behave similarly to the M1.0z3.0 model but, instead, it resembles the M1.0z2.5 and M1.0z2.3 models of higher [Fe/H]. It experiences a PIE during the second pulse and loses all its mass quickly, aborting a standard TP-AGB phase. This is because its metallicity in mass fraction is similar to those of the M1.0z2.5 and M1.0z2.3 models. The chemical yields of this model follow the same trend as the other 1 M, models experiencing a PIE (Fig. 14).

thumbnail Fig. 14.

Global overproduction factors (Eq. (9)) for our 1 M models (top panel), 2 M models (middle panel), and 3 M models (bottom panel).

Unlike the M2.0z3.0 model, the M2.0z3.0α model does not experience a standard PIE (as the models in Table 2), as the PIE is prevented due to its higher metallicity (Table 1). This model lies just above the light grey zone in Fig. 3. However, during the second and third thermal pulses, this model experiences weak proton ingestion events with maximal neutron densities of 1011 − 1012 cm−3. Unlike a standard PIE, the AGB structure and evolution is not affected. During these weak PIEs, first-peak heavy elements are synthesized in the pulse and later brought up to the surface by the third dredge-up. The resulting yields are shown by the dashed line in the middle panel of Fig. 14. We note that, like the other 2 M models, the M2.0z3.0α model does not experience an efficient convective s-process because the temperature at the bottom of the pulse is too low (about 325 MK at maximum; see also Sect. 4.5).

Our findings are consistent with those of Cristallo et al. (2016) who also concluded that the α-enhancement tends to suppress the PIE. The higher C and O abundances in the H-shell lead to lower H-burning temperatures and, thus, to an increase in the entropy difference between the H- and He-burning shells. In these conditions, the PIE is hampered.

4.8. Helium-rich ejecta

After a standard thermal pulse, some material of the pulse is mixed in the convective envelope provided the third dredge up occurs. In contrast, after a PIE, the whole helium-rich convective pulse merges with the hydrogen-rich convective envelope. The envelope He enrichment is clearly visible in some of our models (e.g. the M1.0z3.0 and M1.0z2.5 models, Fig. 15). Right after the PIE, the surface helium mass fraction rises up to about 0.34 (0.32) in the M1.0z3.0 and M1.0z2.5 models, respectively. Although present, this effect is much less pronounced in the M2.0z3.0 and M2.0z2.5 models because the helium contained in the pulse is diluted in a much bigger envelope and the pulse is less massive.

thumbnail Fig. 15.

Surface 4He mass fraction as a function of the pulse number. The mass fraction is shown after each thermal pulse, just after the third dredge up, if any.

Although the detection of helium is very challenging, there have been several recent direct detections of helium enrichment in stars. Pasquini et al. (2011) managed to estimate the abundance for two giants stars ([Fe/H] = −1.22 and −1.08) in the globular cluster NGC 2808. They found an enrichment of 0.17 in helium mass fraction, meaning that if one of the two star has a solar He abundance, the other would have X(4He) ∼ 0.42. Dupree & Avrett (2013) reported helium abundances of X(4He) < 0.22 and 0.39 < X(4He) < 0.44 in two giant stars ([Fe/H] = −1.86 and −1.79) belonging to the globular cluster ω Centauri. In NGC 2808, Marino et al. (2014) found X(4He) = 0.34 in a sample of horizontal branch stars with metallicities −1.5 <  [Fe/H] <  −1. The origin of this helium-rich population in Globular Clusters is debated. Possible polluters are intermediate-mass AGB stars (between ∼3 and ∼7 MVentura & D’Antona 2009), super-AGB stars (Pumo et al. 2008; Siess 2010) massive rotating stars (Decressin et al. 2007), or massive binary stars (de Mink et al. 2009). The present work suggests that low-mass AGB stars experiencing a PIE may also be interesting candidates. Nevertheless, PIEs do not seem to occur in AGB stars of metallicity higher than [Fe/H] ≃ −2.3 (Fig. 3), while some helium-rich stars in globular clusters have [Fe/H] ≃ −1 (cf. previous discussion). We recall, however, that Fig. 3 mostly considers models without extra mixing processes. It remains to be seen if AGB models including overshoot could experience PIEs at [Fe/H] ≃ −1. This possibility will be explored in a future work.

4.9. Yields and i-process signatures

We computed the yields and associated over-production factors of our models according to Eqs. (8) and (9). Results are displayed in Figs. 14 and 16.

thumbnail Fig. 16.

Global overproduction factors (Eq. (9)) for the nine models with M = 1, 2 and 3 M, and [Fe/H] = −3, −2.5 and −2. In each panel (corresponding to a given element), the x-axis (y-axis) shows the initial metallicity (initial mass) of the AGB model. The overproduction factors are color-coded according to the colorbar on the right.

At [Fe/H] = −2, neither the 1 M, nor the 2 M, show any enhancement in heavy elements (blue patterns). Unlike all other 1 and 2 M, models, these two models do not experience any PIE and their pulse temperature is two low to experience significant convective s-process (cf. Sect. 4.5). We also recall that in the absence of additional mixing below the convective envelope, our models do not experience radiative s-process nucleosynthesis during the interpulse phase.

As discussed in Sect. 4.5, the 3 M, models experience an efficient convective s-process which leads to the production of the light s-process elements Sr, Y, and Zr with similar overproduction factors (Fig. 14, bottom panel). The other models have indeed experienced a PIE and the yields result from a complex interplay between the characteristics of the PIE (mass of ingested protons, pulse temperature, maximum neutron density, split, etc.), the efficiency of the third dredge up (the λ parameter), and the dilution of the pulse material into the envelope of variable size and composition.

Despite these differences, all these i-process models show a similar rising distribution from Fe to Bi (Fig. 14, top, and middle panels). The overproduction factors of the 2 M models are lower on average than those of the 1 M models because of the larger dilution of the pulse material. Models at metallicity [Fe/H] = −3.0 (solid black) and −2.5 (solid red) give similar overproduction factors but the yields in M are ∼3 times higher in the [Fe/H] = −2.5 models because of the difference in metallicity. The M1.0z2.3 model (green) shows higher overproduction factors except for Pb and Bi. The highest overproduction factor is found for Tantalum (Z = 73) in this model.

Elements above Nb are only synthesized by low-mass low-metallicity models (Fig. 16). In contrast, lighter elements are also synthesized by 3 M models, thanks to the convective s-process (cf. Sect. 4.5). Nuclei between Ba and Bi are more produced in the 1 M models than in 2 M models. Xenon (Xe), ytterbium (Yb), tantalum (Ta), lead (Pb), and bismuth (Bi) are among the most produced species. A key feature of these models is the high xenon over iodine ratio (also visible in Fig. 14, top and middle panels, at Z = 53 − 54). Below iodine (Z = 53), the overproduction factors do not exceed 102 (blue and orange squares). Above xenon (Z = 54), the overproduction factors are often in the range of 102 − 103 (green) and sometimes higher (purple and grey). The high overproduction of elements with Z >  53 compared to elements with Z ≤ 53 is a significantly different feature that is observed between the i- and s-processes (as shown in Choplin et al. 2021, their Fig. 7).

5. Summary and conclusions

In this third paper of the series, we study the i-process originating from a PIE in AGB stars of 1, 2, and 3 M at metallicities of [Fe/H] = −3, −2.5, −2.3, and −2. We used the code STAREVOL with a network of 1160 species and provide the first i-process yields8 from a grid of AGB models experiencing PIEs. We note that extra mixing was not considered in these calculations.

We find that PIEs happen in six out of our 12 AGB stars models, preferentially at low-mass and low-metallicity. PIEs arise in 1 and 2 M AGB models with [Fe/H] = −3, −2.5, and −2.3 during the first or second thermal pulse. In these models, ≃10−6 − 10−4 M of hydrogen is ingested in the convective pulse. This activates hydrogen-burning at very high temperature, triggers i-process nucleosynthesis characterized by neutron densities of ≃1014 − 1015 cm−3. This event strongly impacts the stellar structure and subsequent AGB evolution.

In our 1 M models with [Fe/H] = −2.5 and −2.3, the AGB phase ends right after the PIE because of the strong effect of molecular CO opacities that leads to the quick loss of the whole convective envelope before any further thermal pulse can develop. Our 1 M model with [Fe/H] = −3 is special in the sense that the PIE happens during the very first thermal pulse, at relatively low temperature. This leads to a different evolutionary pathway: the surface is less enriched, the molecular CO opacities smaller and hence the mass loss weaker. This results in the resumption of the AGB phase with ∼40 further thermal pulses after the PIE. In our 2 M models, the PIE pulse material is more diluted in the more massive envelope, reducing the surface enrichment and allowing for the development of about ten subsequent thermal pulses.

All our PIE models, except the 1 M model with [Fe/H] = −3, experience a proper split of the convective pulse. We highlighted the dependence of the pulse splitting on the convective velocity (the lower the convective velocity, the higher in mass the split occurs) and we show that (at least in our models), the resulting surface enrichment does not depend on the convective velocity unless it is unrealistically low. This result suggests that PIEs can happen similarly under various convective recipes or assumptions. In addition, we showed that in all our PIE models, there are more than 50 turnover timescales (the e factor, cf. Eq. (11) and Table 2) between i-process nucleosynthesis and the split. This means that the surface of our AGB models is fully enriched in i-process products right after the PIE.

The yields of our PIE models show similar rising trends with overproduction factors of up to 105 for Pb and Bi. The elements Xe, Yb, Ta, Pb, and Bi show the highest overproduction factors. Also, overproduction factors of elements with 56 <  Z <  83 are higher in 1 M than in 2 M, AGB models. Our 1 M PIE models are also strong producers of helium. Finally, we found that our 1 M model with [Fe/H] = −3 and α-enhanced mixture behaves like our 1 M model with [Fe/H] = −2.5 and =−2.3 of similar metallicity (although a lower [Fe/H]).

Our 3 M models do not experience i-process but convective s-process and efficient third dredge up events. Their yields consequently show overproduction from the first neutron peak with factors of up to about 40.

This study has shown that AGB stars of low mass and low metallicity can undergo PIE and eject heavy elements made from i-process nucleosynthesis. The evolutionary pathway after a PIE depends on the details of the PIE itself (e.g. pulse where it occurs, mass of protons ingested). It is interesting to note that our AGB models make heavy elements without any parametrised extra mixing process such as overshoot and without the artificial inclusion of a 13C-pocket. Nevertheless, similar models with extra mixing deserve to be computed and analyzed. Overshooting could take place at the bottom of the envelope and at the bottom and top of the convective pulse. As it enhances mixing, it may facilitate ingestion of protons into the convective pulse. But the modelling of this process also include badly constrained parameters; the impact of overshoot on proton ingestion events is therefore complex and deserves a detailed study. This will be explored in a future work.


3

Pulses arising at the very early stage of the TP-AGB phase, that are not fully developed, are not counted as proper pulses. A pulse is referred to as such if the maximal temperature at the bottom of the He-driven convective zone is greater than 2 × 108 K.

4

The 12C mass fractions at the bottom of the PIE pulse and at peak of neutron density are 0.06 and 0.17 in the M1.0z3.0 and M1.0z2.5 models, respectively. For 16O, these are 4.1 × 10−3 and 0.016.

5

The envelope and pulse masses are about 0.35 and 0.04 M for the 1 M models and 1.31 and 0.017 M for the 2 M models, respectively.

6

In Paper I, we reported a value of tsplit − tNnmax = 0.002 yr for this model (the caption of Fig. 2). This leads to e ∼ 15. The difference is due to the fact that in the present paper, we recomputed this model with the latest version of STAREVOL for a sake of homogeneity. Slight changes were noticed between the old and new M1.0z2.5 model.

7

We note that the λ parameter is only defined for standard thermal pulses, not when a PIE occurs.

8

The yields are available online at http://www.astro.ulb.ac.be/~siess/Site/StellarModels

Acknowledgments

This work was supported by the Fonds de la Recherche Scientifique-FNRS under Grant No IISN 4.4502.19. L.S. and S.G. are senior FRS-F.N.R.S. research associates.

References

  1. Angulo, C., Arnould, M., Rayet, M., et al. 1999, Nucl. Phys. A, 656, 3 [Google Scholar]
  2. Arnould, M., & Goriely, S. 2006, Nucl. Phys. A, 777, 157 [Google Scholar]
  3. Arnould, M., & Goriely, S. 2020, Progr. Part. Nucl. Phys., 112, 103766 [NASA ADS] [CrossRef] [Google Scholar]
  4. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  5. Banerjee, P., Qian, Y.-Z., & Heger, A. 2018, ApJ, 865, 120 [Google Scholar]
  6. Campbell, S. W., & Lattanzio, J. C. 2008, A&A, 490, 769 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Campbell, S. W., Lugaro, M., & Karakas, A. I. 2010, A&A, 522, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Cassisi, S., Castellani, V., & Tornambe, A. 1996, ApJ, 459, 298 [NASA ADS] [CrossRef] [Google Scholar]
  9. Caughlan, G. R., & Fowler, W. A. 1988, At. Data Nucl. Data Tables, 40, 283 [NASA ADS] [CrossRef] [Google Scholar]
  10. Chieffi, A., Domínguez, I., Limongi, M., & Straniero, O. 2001, ApJ, 554, 1159 [Google Scholar]
  11. Choplin, A., Siess, L., & Goriely, S. 2021, A&A, 648, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Clarkson, O., & Herwig, F. 2021, MNRAS, 500, 2685 [Google Scholar]
  13. Clarkson, O., Herwig, F., & Pignatari, M. 2018, MNRAS, 474, L37 [Google Scholar]
  14. Côté, B., Denissenkov, P., Herwig, F., et al. 2018, ApJ, 854, 105 [Google Scholar]
  15. Cowan, J. J., & Rose, W. K. 1977, ApJ, 212, 149 [Google Scholar]
  16. Cristallo, S., Piersanti, L., Straniero, O., et al. 2009, PASA, 26, 139 [Google Scholar]
  17. Cristallo, S., Karinkuzhi, D., Goswami, A., Piersanti, L., & Gobrecht, D. 2016, ApJ, 833, 181 [Google Scholar]
  18. Cruz, M. A., Serenelli, A., & Weiss, A. 2013, A&A, 559, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Dardelet, L., Ritter, C., Prado, P., et al. 2014, XIII Nuclei in the Cosmos (NIC XIII), 145 [Google Scholar]
  20. de Mink, S. E., Pols, O. R., Langer, N., & Izzard, R. G. 2009, A&A, 507, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Decressin, T., Meynet, G., Charbonnel, C., Prantzos, N., & Ekström, S. 2007, A&A, 464, 1029 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Denissenkov, P. A., Herwig, F., Battino, U., et al. 2017, ApJ, 834, L10 [Google Scholar]
  23. Denissenkov, P. A., Herwig, F., Woodward, P., et al. 2019, MNRAS, 488, 4258 [Google Scholar]
  24. Denissenkov, P. A., Herwig, F., Perdikakis, G., & Schatz, H. 2021, MNRAS, 503, 3913 [NASA ADS] [CrossRef] [Google Scholar]
  25. Despain, K. H., & Scalo, J. M. 1976, ApJ, 208, 789 [NASA ADS] [CrossRef] [Google Scholar]
  26. Dupree, A. K., & Avrett, E. H. 2013, ApJ, 773, L28 [Google Scholar]
  27. Fujimoto, M. Y. 1977, PASJ, 29, 331 [NASA ADS] [Google Scholar]
  28. Fujimoto, M. Y., Iben, I., Jr., Chieffi, A., & Tornambe, A. 1984, ApJ, 287, 749 [NASA ADS] [CrossRef] [Google Scholar]
  29. Fujimoto, M. Y., Iben, I., Jr., & Hollowell, D. 1990, ApJ, 349, 580 [CrossRef] [Google Scholar]
  30. Fujimoto, M. Y., Ikeda, Y., & Iben, I., Jr. 2000, ApJ, 529, L25 [NASA ADS] [CrossRef] [Google Scholar]
  31. Fujiya, W., Hoppe, P., Zinner, E., Pignatari, M., & Herwig, F. 2013, ApJ, 776, L29 [Google Scholar]
  32. Gil-Pons, P., Doherty, C. L., Gutiérrez, J. L., et al. 2018, PASA, 35, e038 [NASA ADS] [CrossRef] [Google Scholar]
  33. Goriely, S., & Siess, L. 2005, in From Lithium to Uranium: Elemental Tracers of Early Cosmic Evolution, eds. V. Hill, P. Francois, & F. Primas, IAU Symp., 228, 451 [Google Scholar]
  34. Goriely, S., & Siess, L. 2018, A&A, 609, A29 [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., Siess, L., & Choplin, A. 2021, A&A, 654, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Goswami, P. P., & Goswami, A. 2022, A&A, 657, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Hampel, M., Stancliffe, R. J., Lugaro, M., & Meyer, B. S. 2016, ApJ, 831, 171 [Google Scholar]
  39. Hampel, M., Karakas, A. I., Stancliffe, R. J., Meyer, B. S., & Lugaro, M. 2019, ApJ, 887, 11 [Google Scholar]
  40. Herwig, F., Pignatari, M., Woodward, P. R., et al. 2011, ApJ, 727, 89 [Google Scholar]
  41. Iben, I., Jr. 1976, ApJ, 208, 165 [NASA ADS] [CrossRef] [Google Scholar]
  42. Iben, I., Jr. 1977, ApJ, 217, 788 [NASA ADS] [CrossRef] [Google Scholar]
  43. Iwamoto, N., Kajino, T., Mathews, G. J., Fujimoto, M. Y., & Aoki, W. 2004, ApJ, 602, 377 [Google Scholar]
  44. Jadhav, M., Pignatari, M., Herwig, F., et al. 2013, ApJ, 777, L27 [NASA ADS] [CrossRef] [Google Scholar]
  45. Jones, S., Ritter, C., Herwig, F., et al. 2016, MNRAS, 455, 3848 [Google Scholar]
  46. Jonsell, K., Barklem, P. S., Gustafsson, B., et al. 2006, A&A, 451, 651 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Karakas, A. I., & Lattanzio, J. C. 2014, PASA, 31, e030 [NASA ADS] [CrossRef] [Google Scholar]
  48. Karakas, A. I., Lattanzio, J. C., & Pols, O. R. 2002, PASA, 19, 515 [NASA ADS] [CrossRef] [Google Scholar]
  49. Karinkuzhi, D., Van Eck, S., Goriely, S., et al. 2021, A&A, 645, A61 [EDP Sciences] [Google Scholar]
  50. Kobayashi, C., Umeda, H., Nomoto, K., Tominaga, N., & Ohkubo, T. 2006, ApJ, 653, 1145 [NASA ADS] [CrossRef] [Google Scholar]
  51. Koch, A., Reichert, M., Hansen, C. J., et al. 2019, A&A, 622, A159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Koning, A. J., & Rochman, D. 2012, Nucl. Data Sheets, 113, 2841 [NASA ADS] [CrossRef] [Google Scholar]
  53. Lau, H. H. B., Stancliffe, R. J., & Tout, C. A. 2009, MNRAS, 396, 1046 [CrossRef] [Google Scholar]
  54. Liu, N., Savina, M. R., Davis, A. M., et al. 2014, ApJ, 786, 66 [NASA ADS] [CrossRef] [Google Scholar]
  55. Lugaro, M., Karakas, A. I., Stancliffe, R. J., & Rijs, C. 2012, ApJ, 747, 2 [Google Scholar]
  56. Marigo, P. 2002, A&A, 387, 507 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Marino, A. F., Milone, A. P., Przybilla, N., et al. 2014, MNRAS, 437, 1609 [NASA ADS] [CrossRef] [Google Scholar]
  58. Mishenina, T., Pignatari, M., Carraro, G., et al. 2015, MNRAS, 446, 3651 [Google Scholar]
  59. Pasquini, L., Mauas, P., Käufl, H. U., & Cacciari, C. 2011, A&A, 531, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Pumo, M. L., D’Antona, F., & Ventura, P. 2008, ApJ, 672, L25 [NASA ADS] [CrossRef] [Google Scholar]
  61. Reddy, B. E., Lambert, D. L., & Allende Prieto, C. 2006, MNRAS, 367, 1329 [Google Scholar]
  62. Reimers, D. 1975, Mem. Soc. R. Sci. Liege, 8, 369 [NASA ADS] [Google Scholar]
  63. Ritter, C., Herwig, F., Jones, S., et al. 2018, MNRAS, 480, 538 [NASA ADS] [CrossRef] [Google Scholar]
  64. Roederer, I. U., Karakas, A. I., Pignatari, M., & Herwig, F. 2016, ApJ, 821, 37 [Google Scholar]
  65. Scalo, J. M. 1979, A&A, 74, 6 [NASA ADS] [Google Scholar]
  66. Schenk, O., & Gärtner, K. 2004, Future Gen. Comput. Syst., 20, 475 [CrossRef] [Google Scholar]
  67. Schlattl, H., Cassisi, S., Salaris, M., & Weiss, A. 2001, ApJ, 559, 1082 [Google Scholar]
  68. Schwarzschild, M., & Härm, R. 1965, ApJ, 142, 855 [NASA ADS] [CrossRef] [Google Scholar]
  69. Schwarzschild, M., & Härm, R. 1967, ApJ, 150, 961 [NASA ADS] [CrossRef] [Google Scholar]
  70. Siess, L. 2006, A&A, 448, 717 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Siess, L. 2007, A&A, 476, 893 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Siess, L. 2010, A&A, 512, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 [Google Scholar]
  74. Siess, L., Livio, M., & Lattanzio, J. 2002, ApJ, 570, 329 [Google Scholar]
  75. Stancliffe, R. J., Dearborn, D. S. P., Lattanzio, J. C., Heap, S. A., & Campbell, S. W. 2011, ApJ, 742, 121 [Google Scholar]
  76. Stephens, D., Herwig, F., Woodward, P., et al. 2021, MNRAS, 504, 744 [NASA ADS] [CrossRef] [Google Scholar]
  77. Suda, T., & Fujimoto, M. Y. 2010, MNRAS, 405, 177 [NASA ADS] [Google Scholar]
  78. Vassiliadis, E., & Wood, P. R. 1993, ApJ, 413, 641 [Google Scholar]
  79. Ventura, P., & D’Antona, F. 2009, A&A, 499, 835 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Xu, Y., Goriely, S., Jorissen, A., Chen, G., & Arnould, M. 2013, A&A, 549, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1.

Main characteristics of the models computed in this work, with the initial mass Mini, [Fe/H] ratio, initial metallicity in mass fraction Z, initial chemical composition (solar-scaled or α-enhanced), total lifetime τstar, duration of the AGB phase τAGB, final mass Mfin, and final mass of the convective envelope , number of pulses Npulse, along with information on whether the model experiences a PIE or not and the pulse number during which the PIE arises.

Table 2.

Important characteristics of our models experiencing a PIE.

Table 3.

Surface 12C/13C, C/N, O/N, and C/O ratios in number just after the PIE.

All Figures

thumbnail Fig. 1.

Evolution of the M1.0z3.0 model during the PIE: Kippenhahn diagram (top panel), temperature at the bottom of the convective thermal pulse (or envelope for t − t0 ≳ 1 yr, second panel), energy released by the main reactions at the bottom of the pulse (or envelope, third panel), maximal neutron density (fourth panel). Bottom panel: 12C/13C, C/N, and O/N ratios at the bottom of the convective zone (solid lines) and at the stellar surface (dashed lines). The dashed green and blue lines in the Kippenhahn diagram (top panel) delineate the hydrogen and helium-burning zones (where the nuclear energy production by H- and He-burning exceeds 10 erg g−1 s−1). The thin light-grey lines indicate iso-masses expressed in M.

In the text
thumbnail Fig. 2.

Entropy barrier between the top of the convective thermal pulse and the bottom of the H-rich layers as a function of the pulse number (see text for details on the computation of ds). The empty circles mark the occurrence of the PIE.

In the text
thumbnail Fig. 3.

Mass-metallicity diagram showing the occurrence of PIEs during the early AGB phase of models from various authors. Filled symbols show models experiencing a PIE while empty symbols are for models that do not experience a PIE. The dark grey zone shows the approximate region where PIEs happen in all of the models and the light grey one where PIEs happen in most of the models. The corresponding [Fe/H] ratios are indicated on the right axis assuming solar-scaled mixtures. Models are from Iwamoto et al. (2004, red triangles), Campbell & Lattanzio (2008, magenta circles), Cristallo et al. (2009, blue triangle), Lau et al. (2009, blue diamonds), Suda & Fujimoto (2010, orange triangles), Cristallo et al. (2016, black pentagons). Model results from this work are shown as green squares. All models were computed without extra mixing processes, except the models from Cristallo et al. (2009, 2016), which consider overshooting below the convective envelope.

In the text
thumbnail Fig. 4.

Energy generation from the main reactions during the PIE in the M2.0z3.0 model just before the splitting (top panel), and just after the splitting (bottom panel). The total nuclear energy production rate is shown by the thick grey line. The protons and neutrons mass fractions are shown.

In the text
thumbnail Fig. 5.

Energy generation from the main reactions during the PIE in the M1.0z3.0 model just before the splitting (top panel), just after the splitting (middle panel) and just after both parts of the pulse have merged again (bottom panel). The total nuclear energy production rate is shown by the thick grey line. The protons and neutrons mass fractions are shown.

In the text
thumbnail Fig. 6.

Kippenhahn diagram and maximal neutron density of the 1 M model with [Fe/H] = −3.0 during the PIE while considering different values for the convective velocity: the standard case (top left), vconv divided by 10 (top right), 100 (bottom left), and 1000 (bottom right). The dotted green and blue lines in the Kippenhahn diagram correspond to the maximum nuclear energy production by hydrogen- and helium-burning, respectively. The dashed green and blue lines delineate the H- and He-burning zones (when the production of energy by H- and He-burning exceeds 10 erg g−1 s−1).

In the text
thumbnail Fig. 7.

Maximum neutron density for the M1.0z3.0 model during the PIE for different values for the convective velocity.

In the text
thumbnail Fig. 8.

Surface mass fractions after the PIE for the M1.0z3.0 model for different values for the convective velocity.

In the text
thumbnail Fig. 9.

Main i-process path (starting from 56Fe) in the M1.0z2.3 (blue arrows) and M2.0z3.0 model (green arrows), at the bottom of the pulse, at the time of maximum neutron density (Nn = 2.19 × 1015 and 6.76 × 1013 cm−3, respectively). The dashed lines show the secondary paths, where at least 30% of the total flux goes. The four panels corresponds to four different zones in the (N, Z) plane. The black squares highlight the stable nuclei. The abundances of the M1.0z2.3 model are shown by the red colour scale in mass fraction.

In the text
thumbnail Fig. 10.

Surface elemental mass fractions normalized by their initial mass fractions for the nine models with Mini = 1, 2 and 3 M, and [Fe/H] = −2, −2.5 and −3. The initial mass is increasing from the top to the bottom panels and the initial metallicity is increasing from the left to the right panels. The patterns correspond to the surface abundances after the indicated pulse number. The bold frames highlights models experiencing a PIE. Unstable isotopes are not beta-decayed and the unstable Z = 43 and Z = 61 elements are not considered since Xini = 0 for those elements.

In the text
thumbnail Fig. 11.

Evolution of various quantities as a function of the pulse number. Top panel: maximum temperature reached at the bottom of the pulse. The horizontal dashed line shows the temperature above which the 22Ne(α, n)25Mg reaction becomes efficient. Middle panel: third dredge up efficiency parameter λ (see text for the definition). The λ parameter for the pulses where a PIE occurs is not shown since it cannot be properly defined. Bottom panel: surface 25Mg mass fraction.

In the text
thumbnail Fig. 12.

Surface 12C/13C ratio in number as a function of the pulse number. The ratio is shown at the middle of each pulse (i.e. before the third dredge up, if any).

In the text
thumbnail Fig. 13.

Temperature at the bottom of the convective envelope during the AGB phase. Models of 1, 2, and 3 M, are shown in the top, middle, and bottom panels, respectively.

In the text
thumbnail Fig. 14.

Global overproduction factors (Eq. (9)) for our 1 M models (top panel), 2 M models (middle panel), and 3 M models (bottom panel).

In the text
thumbnail Fig. 15.

Surface 4He mass fraction as a function of the pulse number. The mass fraction is shown after each thermal pulse, just after the third dredge up, if any.

In the text
thumbnail Fig. 16.

Global overproduction factors (Eq. (9)) for the nine models with M = 1, 2 and 3 M, and [Fe/H] = −3, −2.5 and −2. In each panel (corresponding to a given element), the x-axis (y-axis) shows the initial metallicity (initial mass) of the AGB model. The overproduction factors are color-coded according to the colorbar on the right.

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.