Next Article in Journal
Nonsingular Black Holes in 4D Einstein–Gauss–Bonnet Gravity
Next Article in Special Issue
Correction: Buompane et al. Recent Achievements of the ERNA Collaboration. Universe 2022, 8, 135
Previous Article in Journal
Analogue Quantum Gravity in Hyperbolic Metamaterials
Previous Article in Special Issue
Post-AGB Stars as Tracers of AGB Nucleosynthesis: An Update
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Initial-Final Mass Relation of White Dwarfs: A Tool to Calibrate the Third Dredge-Up

Department of Physics and Astronomy, University of Padova, 35122 Padova, Italy
Submission received: 21 February 2022 / Revised: 9 April 2022 / Accepted: 11 April 2022 / Published: 14 April 2022

Abstract

:
The initial mass-final mass relationship (IFMR) of white dwarfs (WD) represents a crucial benchmark for stellar evolution models, especially for the efficiency of mixing episodes and mass loss during the asymptotic giant branch (AGB) phase. In this study, we argue that this relation offers the opportunity to constrain the third dredge-up (3DU), with important consequences for chemical yields. The results are discussed in light of recent studies that have identified a kink in the IFMR for initial masses close to 2 M . Adopting a physically-sound approach in which the efficiency λ of the 3DU varies as a function of core and envelope masses, we calibrate λ in solar-metallicity TP-AGB models in order to reproduce the final masses of their WD progeny, over the range of initial masses 0.9 M i / M 6 . In particular, we find that in low-mass stars with 1.4 M i / M 2.0 the efficiency is small, λ 0.3 , it steeply rises to about λ 0.65 in intermediate-mass stars with 2.0 M i / M 4.0 , and then it drops in massive TP-AGB stars with 4.0 M i / M 6.0 . Our study also suggests that a second kink may show up in the IFMR at the transition between the most massive carbon stars and those that are dominated by hot-bottom burning.

1. Introduction

The initial-final mass relation (IFMR) of white dwarfs plays a key role in several astrophysical applications [1,2]. A recent analysis [3] of a few carbon-oxygen white dwarfs in old open clusters of the Milky Way (MW) identified a kink in the initial-final mass relation (IFMR), located over a range of initial masses, 1.65 M i / M 2.10 , which unexpectedly interrupts the commonly assumed monotonic trend. The proposed interpretation links this observational fact to the formation of carbon stars and the modest outflows (with mass loss rate < 10 7 M / yr ) that are expected as long as the carbon excess, C–O, remains too low to produce carbonaceous dust grains in a sufficient amount. Under these conditions the mass of the carbon-oxygen core can grow more than is generally predicted by stellar models. A new systematic follow-up investigation ([4], to which we refer for all details), based on GAIA DR2 and EDR3, of TP-AGB stars belonging to open clusters of known age support the findings of [3].
The key point in explaining the kink is the interplay between the third dredge-up (3DU), which determines the surface enrichment in carbon, and the dependence of stellar winds for C stars on the excess of carbon with respect to oxygen, C–O [5,6,7]. We leverage this mutual interconnection to constrain the efficiency of 3DU, λ , in solar metallicity TP-AGB stars, using semi-empirical IFMR as the primary calibrator. A similar approach was already introduced by Kalirai et al. [8] and Marigo et al. [3]. In this work we improve the methodology, as we take into account the dependence of λ on main stellar parameters (core mass and envelope mass), on the basis of the indications of complete TP-AGB models [9,10].
The paper is structured as follows. Section 2 details the novel approach to calibrating λ and recall the basic ingredients of our COLIBRI code for the TP-AGB phase. In Section 3 we analyze and discuss the results, with a particular focus on the impact of carbon enrichment on the mass-loss rate (Section 2). We present our derived λ -law at solar metallicity in Section 3.2, and compare it with the predictions of full TP-AGB models in the literature. Section 3.3 analyses the semi-empirical IFMR, compares the data with models, and points to a possible secondary kink drawn by the WDs produced by high-mass AGB stars (Section 3.4). In Section 4 we discuss other relevant quantities that derive from the λ -calibration, namely the TP-AGB lifetimes and the carbon ejecta. We compare our findings with other published models. Finally, Section 5 closes the paper.

2. Materials and Methods

The 3DU influences the IFMR in at least two ways: 1. The mass of the core, M c , is instantaneously reduced at each mixing episode; 2. It affects the mass-loss rate. The higher the efficiency λ 1 of the dredge-up episode, the greater the decrease in M c , equal to λ × Δ M c . The effect on mass loss is a consequence of the variation of the surface chemical composition which induces notable changes in atmospheric opacity [11,12] and hence in the effective temperature [13], a key parameter for the pulsation periods [14], and the winds of AGB stars [15,16]. Furthermore, the increase of C and O in the atmospheres of carbon stars controls the excess of carbon relative to oxygen, C–O, which in turn sets the budget of C atoms available for the formation of carbonaceous dust [5,7,17].
In Marigo et al. [3], we calibrated the average efficiency of the 3DU as a function of the initial mass to reproduce the IFMR and its kink at M i 2 M , under the simple hypothesis that λ remains constant during the TP-AGB phase, at given M i . In this study we aim to improve the description of the 3DU, by relaxing the assumption of constant λ . In principle, the depth of the 3DU depends primarily on the mass of the core, M c , the mass of the envelope, M env , and the chemical composition (e.g., [9,18]). These parameters, in turn, affect other key quantities for the 3DU, such as, for instance, the flash-luminosity peak L P . At this stage the envelope convection zone reaches its maximum inward extension in mass fraction and temperature while, at the same time, the nuclear-processed material is pushed out to its lowest temperature during the thermal-pulse cycle [19]. In principle, the larger L P , the deeper the inward envelope extension.
An important step forward to include these dependencies was already addressed by [20,21], in their population synthesis simulations aimed at reproducing the luminosity functions of AGB stars in the Magellanic Clouds. In this paper, we move to solar metallicities, and adopt the IFMR relation as the main calibrator. To take a more physical approach we rely on a study based on complete TP-AGB models of the FRUITY database2 ([9], see their Equation (3) and Figure 1). There, the authors proposed a fitting formula that accurately reproduces their predictions for the mass of dredged-up material, Δ M 3 DU , at each pulse, as a function of the minimum mass of the core, M c min , for the onset 3DU, the current envelope and core masses, M env and M c , and the metallicity, Z.
We adopt the same functional form introduced by [9], with some modifications, and fit Δ M 3 DU of the entire set of FRUITY TP-AGB models with metallicity Z = 0.014 . Using the original fitting relation in combination with our mass-loss prescriptions in COLIBRI, we fail to reproduce the IFMR. Therefore, we opt to maintain the dependence on M c min , M c , and M env , and introduce a free parameter that allows us to vary the relation to best reproduce the observational IFMR data.
In practice, we apply a corrective factor, f cor , to the FRUITY Δ M 3 DU relation, and make it vary iteratively until we find a good match between models and observations. This procedure is carried out for every COLIBRI stellar model of given initial mass M i , calculating, for each iteration with a test f cor , the entire TP-AGB phase up to the complete ejection of the envelope. When the final M c of the model and the measured WD mass coincide (within a tolerance of about 0.01–0.05 M ), the loop closes and the appropriate value of f cor is found. The results of the calibration are discussed in Section 3.

TP-AGB Models

At this point it is worth mentioning here the main input prescriptions of the COLIBRI code, focusing on those that are particularly relevant for the present study. For more details we refer to [3,13,20,21].
  • Equation of state and opacity.COLIBRI contains the ÆSOPUS3 and Opacity Project4 software packages [11,22] as internal routines, to compute on-the-fly both the equation of state (for ≳ 800 atomic and molecular species) and Rosseland mean opacities at temperatures 10 3 T / K 10 8 . This allows us to consistently describe the impact on the stellar structure of variations in abundances due to nucleosynthesis and mixing.
  • Nucleosynthesis. Hot-bottom burning is followed with a complete nuclear network coupled to a diffusive description of convection. Nuclear reactions of the p-p chains, CNO cycles, Ne-Na and Mg-Al chains are included.
  • Mass loss. Before the onset of dust-driven mass loss we assume that stellar winds are sustained by low-amplitude pulsation [23,24,25]. The resulting rates are typically low, from ≈ 10 9 10 7 M yr−1. In the dust-driven regime the mass-loss rate is computed with a different formalism depending on the photospheric C / O   ratio, namely: a widely-used relation ([26], with an efficiency parameter η B = 0.01 ) for the O-rich stages when C / O < 1 , and a routine based on state-of-the-art pulsating atmosphere models for carbon stars, when C / O > 1 [5,6]. These models, in particular, predict that radiation-driven winds are effectively triggered by carbonaceous dust only when (i) the amount of free carbon, C–O, overcomes a threshold (C–O)min, and (ii) suitable conditions for dust formation (in terms of stellar luminosity, effective temperature and mass) are met in the extended atmospheres.
  • Intershell composition. A nuclear network with the main α -capture reactions is included. The typical abundance pattern is characterized by He / C / O = 0.70–0.78/0.20–0.25/0.005–0.02 (abundances in mass fraction). No overshoot is applied to the bottom of the pulse-driven convection zone.

3. Results

3.1. The Impact of Carbon Enrichment on the Mass-Loss Rate

Two key aspects characterize our COLIBRI models, namely: the computation on-the-fly of gas opacities correctly coupled to the changes in the envelope chemical composition, and the adoption of a mass-loss prescription for the carbon-rich phase that is explicitly dependent on carbon excess, C–O (see Section 2). This allows us to explore the feedback due to carbon enrichment both on the stellar structure and on mass loss, and their mutual interaction as well.
Figure 1 shows the evolution of luminosity, effective temperature and mass-loss rate, in two models that become both carbon stars and end their TP-AGB phase with a bare core of about M f 0.7 M . First, we note that when C / O increases due to the 3DU and enters the narrow range 0.98 C / O 1.05 (shades from light pink to white), which characterizes the stars belonging to the advanced S and SC spectral classes, the tracks in the HR diagram tend to heat up slightly (due to a minimum in opacity, since most of C and O atoms are locked in the CO molecule; see Marigo and Aringer [11]). Then, when C / O exceeds 1.05 and attains larger values as a consequence of repeated 3DU episodes, the atmospheric opacity increases making the effective temperature decrease sizeably, a well-known property that distinguishes the class of carbon stars [12]. We see that the model with M i = 2.8 M achieves a higher surface enrichment in carbon, and therefore the track cools more than the M i = 1.8 M model. In turn, the increase in C/O and the consequent decrease in T eff favor the increase in mass loss (bottom panels). In both tracks, the minimum in T eff is reached when a large part of the envelope has been ejected, i.e., as soon as the envelope mass falls below 35–40% of the total mass ( M env / M < 0.35–0.40). From that stage on, the tracks reverse their evolutionary direction, and begin to heat up as the last thermal pulses occur. In the present models the evolutionary calculations stop when M env < 0.05 M .
Let us now move to discuss the evolution of the mass-loss rate (bottom panels of Figure 1). During the O-rich phases it gradually increases, as expected from the [26] formula. Then, as soon as the model passes through the transition from C / O < 1 to C / O > 1 , the mass-loss rate drops and remains moderate ( M ˙ < 5 × 10 7 M yr 1 ) until the carbon enrichment, resulting from the 3DU, reaches and exceeds the limit necessary to activate a dust-driven wind.
This threshold is not constant, but depends on stellar parameters: in general it increases with the current mass of the star and its effective temperature. In other words, more massive TP-AGB stars must enrich more in carbon to form the minimum amount of carbonaceous dust necessary to trigger the radiative wind. This is evident from the comparison of panels c and d: the model with M i = 2.8 M must attain C / O 1.8 to power a dust-driven wind, while in the model with M i = 1.8 M this happens when C / O 1.3 . The fact that the model with M i = 2.8 M must exceed a higher C–O threshold explains why the dust-driven wind regime (with M ˙ 10 6 M yr 1 ) is activated about 13 TPs after the transition to the C-star regime, while in the model with M i = 1.8 M the same mass-loss rate is reached after 9 TPs. In fact, although the efficiency of 3DU for M i = 2.8 M is greater, the carbon surface enrichment is slower from pulse to pulse, as the dredged-up material is diluted in a more massive envelope. Overall, the model with M i = 2.8 M achieves a higher carbon enrichment and ends the TP-AGB phase with C / O = 2.84 , while the model with M i = 1.8 M terminates its life with C / O = 1.40 .
In both cases the main consequences are: 1. A transient period of modest mass loss during the C-rich phase in which the star evolves along its Hayashi line at approximately constant mass and without significant dust production; 2. A prolongation of the TP-AGB lifetime compared to the predictions obtained with mass loss formulas that do not depend on the surface carbon excess. In particular, for the model with M i = 1.8 M , this circumstance allows the core to grow more than usually predicted by stellar models in the literature. As extensively discussed by Marigo et al. [3], this feature is of key importance to explain the IFMR kink around M i 2 M .

3.2. A Calibrated λ -Law at Solar Metallicity

Taking advantage of the fact that our mass-loss prescription for carbon stars is sensitive to the surface carbon and oxygen abundances, we have the opportunity to introduce a methodology to calibrate the efficiency parameter, λ , of the 3DU, using the IFMR as an observable reference. The details of the procedure have already been illustrated in Section 2. We only recall that, compared to the approach used in [3,8], in this study we make a substantial improvement and include the dependence of λ on the main stellar parameters ( M c min , M c , M env ), based on complete TP-AGB models [9].
The results of the 3DU calibration are shown in Figure 2 (left panel). As expected, for each model, the efficiency of the 3DU first increases (following the growth of the pulse strength), attains a maximum value, λ max , and then it decreases eventually reaching zero, due to the reduction of the envelope mass by stellar winds.
Some interesting aspects emerge from the analysis:
  • Low-mass stars with M i 2 M —those experimenting with the He-flash at the tip of the red giant branch—have a low 3DU efficiency, with λ 0.3 .
  • Intermediate-mass stars with 2 M i 6 —those that avoid electronic degeneracy in their helium cores—undergo a more efficient dredge-up. The λ parameter becomes more than twice as large ( λ max 0.6–0.7) in TP-AGB stars with 2 M i / M 4.5 .
  • At higher masses, 4.5 M i M 6.0 , the efficiency decreases reaching zero at the largest M c .
With respect to the latter point, we note that the decrease of λ in the most massive TP-AGB stars follows from the constraint of obtaining sufficiently high final masses, in agreement with those measured for WDs members of young open clusters ( M f > 0.9–1.0  M ). From a physical point of view this could be explained by the fact that relatively lower efficiencies of the TDU are predicted at increasing core mass and hence shorter interpulse periods, as a consequence of the weaker thermal pulses [10,27]. Furthermore, following the analysis of [28] the combination of hot dredge-up ([29]) and HBB limits the occurrence of the third dredge-up in stars with M i > 4–5 M . By extrapolating the λ -calibration to higher masses ( M i > 6 M ), we could speculate that 3DU might be virtually absent in Super-AGB stars. However, this prediction requires a further observational verification, such as, for example, the absence or scarcity of heavy elements produced by the s-process in the atmosphere of these very luminous red super-giants.
We can now compare our results with predictions from full TP-AGB models. Panel c of Figure 2 shows the trend of the 3DU maximum efficiency, λ max , as a function of the initial mass of the star, as predicted by different sets of models. We note that our calibration has a qualitative trend similar to the FRUITY and Ventura models: it increases with the initial mass, touches a maximum at M i 2.5–3.0 M , and then drops at the highest masses. On the other hand, when we quantitatively compare the results, substantial differences appear, with the largest deviations in λ max showing up for M i > 3 M . For instance, taking the M i = 6 M model, Ventura [30] predicts λ max 0.1 , FRUITY [10] yields λ max 0.4 , while the Monash model of Karakas [31] has λ max 0.93 . In general, our calibrated λ max is closer to FRUITY and Ventura models, while it deviates sizeably from the Monash models, being much lower for M i > 2.5 M .

3.3. The IFMR: Data vs. Stellar Models

Figure 3 (panel a) shows the results of our λ -calibration at solar metallicity: the semi-empirical IFMR is well reproduced by the COLIBRI models all over the range 0.9 M i / M 6 .
The minimum initial mass for carbon star formation is M i 1.5 M , in agreement with the observational evidence from Galactic open clusters [4] and inference from studies on galactic enrichment and solar system formation [33]. The WDs that populate the IFMR kink should be the progeny of low-mass carbon stars with 1.5 M i / M 2.0 , which are characterized by a shallow 3DU and a mild carbon enrichment. The final C / O lies in the range of 1.2–1.6 in these models.
Figure 3. Comparison between theoretical and observed IFMR [3,34]. The black solid line is a fit to the data. (a) Results of this work, obtained with the calibrated λ values, color-coded as a function of the final C / O (right color bar). The core mass at the first TP is also shown (magenta line). (b): Predictions from various authors, namely: ([30], Ventura, purple stars); ([31,32], Monash, red circles); ([10], FRUITY, blue squares); ([35], Miller–Bertolami, green triangles).
Figure 3. Comparison between theoretical and observed IFMR [3,34]. The black solid line is a fit to the data. (a) Results of this work, obtained with the calibrated λ values, color-coded as a function of the final C / O (right color bar). The core mass at the first TP is also shown (magenta line). (b): Predictions from various authors, namely: ([30], Ventura, purple stars); ([31,32], Monash, red circles); ([10], FRUITY, blue squares); ([35], Miller–Bertolami, green triangles).
Universe 08 00243 g003
In the domain of intermediate-mass stars ( 2 M i / M 4.2 ), at increasing M i the 3DU is more efficient, and the surface composition becomes more enriched in carbon. These models are expected to become carbon stars, with a final carbon-to-oxygen ratio of 2 C / O 3.4 . In our grid, the M i = 4.4 M model represents the most massive model that makes the transition to a carbon star and remains in this condition until the end of its evolution. This result is in line with the most massive carbon star found in open clusters: BM IV 34, member of Haffner 14, has an estimated initial mass of ≃3.4–4.0 M . At higher initial masses, ( 4.4 M i / M 6.0 ), the models end the TP-AGB phase as O-rich stars. They all undergo hot-bottom burning, which becomes stronger at increasing M i . However, we note that there is a narrow mass interval, 4.6 M i / M 4.8 , in which the models experience a temporary C-rich phase, before being converted back to O-rich stars by HBB. This circumstance may leave a fingerprint in the IFMR, as discussed in Section 3.4.
We can now examine the theoretical IFMRs published by different authors (panel b of Figure 3). We see that models generally tend to underestimate the final masses of white dwarfs. The Miller–Bertolami relation [35] is the one that most differs from the observations, while the Ventura relation [30] is the one that comes closest to the data. No model reproduces the kink observed around M i 2 M . The reason is that none of the evolutionary calculations adopts a mass-loss prescription for carbon stars that is explicitly dependent on the carbon excess.

3.4. A Second IFMR Kink at Higher Masses?

Looking at Figure 3 (panel a) we note an interesting peculiarity: around the transition mass between carbon stars and O-rich stars experiencing HBB, a second IFMR kink could appear. As can be seen, the M i = 4.4 M model ends its life as a carbon star and produces a white dwarf with M f 1.06 M , while the M i = 4.6 , 4.8 , 5.0 M models terminate the TP-AGB phase as O-rich stars and leave progressively less massive white dwarfs, with M f 1.00 , 0.98 , 0.95 M , respectively. The reason lies in the different efficiency of stellar winds as a function of the surface chemical composition, in particular of the C / O ratio. As to the M i = 4.4 M model, it is necessary that the 3DU enriches the envelope with enough carbon to start the dust-driven wind. In general, the greater the mass of the star, the stronger the inward gravitational force, and therefore the higher the minimum carbon threshold that must be reached. To meet this condition the star needs to experience a relatively large number of mixing episodes, and therefore, in the mean time, the masser of the core can grow.
Moving at a higher stellar mass, the two models with M i = 4.6 , 4.8 M , before leaving the AGB as O-rich stars, experience a transient and short C-rich phase, with C / O > 1 . This circumstance can be appreciated in Figure 4. The M i = 4.6 M model is massive and bright enough to experience HBB, as well as a modest 3DU. Both processes can affect the surface C / O in opposite ways: The 3DU which makes it grow instantaneously at each dredge-up episode, while HBB tends to lower it during the brighter phases of the interpulse periods. Initially the competition is won by the 3DU, so that the star enters the C-star domain after about 40 TPs. During these C-rich phases, lasting ≃0.33 Myr, the mass-loss rate undergoes a substantial slowdown for the reasons already discussed, and the evolution proceeds at an almost constant mass (middle panel).
In the mean time, the growth of the core mass together with the increase in luminosity make the temperature of the deeper layers of the convective envelope increase. HBB becomes stronger and the surface C / O ratio, after reaching a maximum of about 1.62, begins to decline due to the conversion of 12C and 13C into 14N by the CN cycle (middle panel). At the 102 nd TP the C / O falls below one and the mass-loss rate switches to the Bloecker [26] prescription, which depends significantly on luminosity. At this stage, L 34,000 L and the mass-loss rate increases suddenly, leading to a rapid termination of the TP-AGB phase with a final mass M f 1.00 M .
The M i = 4.8 M model experiences a similar evolution, but being more massive and brighter, the HBB process is more efficient and therefore the duration of the C-rich phase, when mass loss undergoes a stop, is shorter (≃0.18 Myr). As a consequence, the final mass is smaller, M f 0.98 M . Finally, starting from the M i = 5 M model and upwards in mass, the C-rich transient phase disappears as HBB dominates, and therefore the theoretical IFMR resumes an increasing monotonic trend.
In short, our study indicates that the application of a mass-loss prescription dependent on carbon excess may produce two kinks in the IFMR of white dwarfs:
  • A first kink starts near the minimum mass for star formation at carbon, peaking at the transition mass between low- and intermediate-mass stars, covering approximately the 1.6 M i / M 2.2 range;
  • A second kink shows up around the transition mass between the highest-mass carbon stars and those that remain O-rich due to the HBB process. Our models indicate a range 4.2 M i / M 4.8 , but the exact prediction may change depending on the assumed mass-loss prescriptions, and on the efficiencies of the 3DU and HBB (hence on the adopted convection theory and/or mixing-length parameter). A thorough study that aims to explore the fine details of the second IFMR kink is postponed for a future work.
While there is observational evidence for the existence of the first kink [3,4], the situation is less clear for the second kink. We only note that the data in Figure 3 show a large dispersion in WD masses of between 4.5 M and 6 M , with some white dwarfs reaching M f 1.1 M at M i 5 M . Were the second massive kink confirmed, it could serve as a powerful calibrator of the HBB efficiency, which is currently very uncertain in AGB star modelling. At the same time, we should also recall that the scatter in WD masses at large M i may also be attributed to a dispersion of initial rotational velocities and convective core overshoot in the evolutionary phases prior the AGB ([36], see the discussion in). Both processes in main-sequence stars, rotational mixing and convective overshoot, lead to creating more massive cores and hence more massive WDs.

4. TP-AGB Lifetimes and Chemical Ejecta

The λ -calibration affects the entire TP-AGB evolution. Here we focus on lifetimes and chemical yields. Figure 5 (panel a) compares the duration of the TP-AGB phase, as a function of the initial mass for solar metallicity. We find that all models predict τ TP AGB to be an initially increasing function with M i reaching a peak, and then decreasing to larger masses. The position of the peak lies in the range 1.8 M i / M 2.5 , while its height, that is, the maximum TP-AGB lifetime, can vary by several factors, from about 1.6 Myr [35] to ≃5.2 Myr (this work).
We note that the relation τ TP AGB ( M i ) predicted by our calibration shows a primary peak of ≃5.2 Myr at M i 1.8 M , and a smaller peak of ≃0.9 Myr for M i = 4.4 M . The latter is produced by the temporary phase C-rich already discussed in Section 3.4, which is responsible for the second kink in the theoretical IFMR report (Figure 3, panel a). These predictions should be tested through population synthesis simulations to be compared with observational data (star counts, luminosity functions), similarly to the studies of Pastorelli et al. [20,21] for the Magellanic Clouds.
The number of dredge-up events, their efficiency and the total duration of the TP-AGB phase impact on chemical yields. In panel b of Figure 5 we compare the 12C ejecta calculated by different authors. As can be seen, for M i 2 M differences exist but are relatively moderate, while large discrepancies among the authors occur for M i > 2 , M , due to different efficiencies of 3DU, HBB and different mass loss prescriptions adopted in the computations. A more in-depth discussion of chemical yields, in light of our new models, will be carried out in a future study.

5. Concluding Remarks

Benefiting from recent systematic studies of the IFMR [3,34] and of state-of the art mass loss prescriptions for carbon stars [5,6], in this work we have introduced a methodology to calibrate the efficiency of the 3DU at solar metallicity, as a function of the parameters M c min , M c , and M env . Our calibration indicates that low-mass stars are characterized by low λ 0.3 , while the efficiency increases significantly in the class of intermediate-mass stars, and eventually decreases in massive AGB stars. We emphasize that the calibration depends on the mass loss used in the models, which underlines the utmost importance of having reliable and up-to-date mass loss prescriptions to follow the TP-AGB phase.
Our models are able to reproduce well the IFMR kink around M i 2 , M , first identified by Marigo et al. [3]. We also suggest the possible existence of a secondary kink, located in the approximate range 4 M i / M 5 , at the transition between the most massive carbon stars and those that remain O-rich due to the onset of HBB. This feature needs further theoretical and observational investigation to be confirmed or rejected.
Finally, we remark again that these results apply to solar-like metallicities. In fact, one should consider that, in principle, the IFMR depends on metallicity (see, for example, [20]) since the efficiencies of both the 3DU and AGB winds are affected by the envelope chemical composition. We will investigate this important aspect in a follow-up study.

Funding

This research was funded by the ERC Consolidator Grant, Project STARKEY, grant agreement n. 615604, and from PRD 2021, University of Padova.

Data Availability Statement

In this study we made use of data available through the FRUITY public website at http://fruity.oa-teramo.inaf.it/ (accessed on 9 April 2022).

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
IFMRinitial-final mass relation
TP-AGBthermally pulsing asymptotic giant branch
TPthermal pulse
WDwhite dwarf
3DUthird dredge-up
Δ M 3 DU dredged-up mass
HBBhot-bottom burning
M i initial mass
M f final mass
M c core mass
M c min minimum core mass for the 3DU
M env envelope mass
T bce temperature at the base of the convective envelope
C / O surface carbon-to-oxygen ratio
C–O carbon excess with respect to oxygen

Notes

1
The efficiency of a 3DU event is commonly described by the dimensionless parameter λ = Δ M 3 DU / Δ M c , defined as the amount of dredged-up material, Δ M 3 DU , relative to the growth of the core mass, Δ M c , during the previous inter-pulse period.
2
Full-network Repository of Updated Isotopic Tables and Yields: http://fruity.oa-teramo.inaf.it/ (accessed on 9 April 2022).
3
http://stev.oapd.inaf.it/cgi-bin/aesopus (accessed on 9 April 2022).
4

References

  1. Marigo, P. Asymptotic Giant Branch Evolution and the Initial-Final Mass Relation of Single CO White Dwarfs. In Binary Paths to Type Ia Supernovae Explosions; Di Stefano, R., Orio, M., Moe, M., Eds.; Cambridge University Press: Cambridge, UK, 2013; Volume 281, pp. 36–43. [Google Scholar] [CrossRef] [Green Version]
  2. Salaris, M.; Serenelli, A.; Weiss, A.; Miller Bertolami, M. Semi-empirical White Dwarf Initial-Final Mass Relationships: A thorough Analysis of Systematic Uncertainties Due to Stellar Evolution Models. Astrophys. J. 2009, 692, 1013–1032. [Google Scholar] [CrossRef] [Green Version]
  3. Marigo, P.; Cummings, J.D.; Curtis, J.L.; Kalirai, J.; Chen, Y.; Tremblay, P.E.; Ramirez-Ruiz, E.; Bergeron, P.; Bladh, S.; Bressan, A.; et al. Carbon star formation as seen through the non-monotonic initial-final mass relation. Nat. Astron. 2020, 4, 1102–1110. [Google Scholar] [CrossRef]
  4. Marigo, P.; Bossini, D.; Trabucchi, M.; Addari, F.; Girardi, L.; Cummings, J.D.; Pastorelli, G.; Dal Tio, P.; Costa, G.; Bressan, A. A Fresh Look at AGB Stars in Galactic Open Clusters with Gaia: Impact on Stellar Models and the Initial-Final Mass Relation. Astrophys. J. Suppl. 2022, 258, 43. [Google Scholar] [CrossRef]
  5. Mattsson, L.; Wahlin, R.; Höfner, S. Dust driven mass loss from carbon stars as a function of stellar parameters. I. A grid of solar-metallicity wind models. Astron. Astrophys. 2010, 509, A14. [Google Scholar] [CrossRef]
  6. Eriksson, K.; Nowotny, W.; Höfner, S.; Aringer, B.; Wachter, A. Synthetic photometry for carbon-rich giants. IV. An extensive grid of dynamic atmosphere and wind models. Astron. Astrophys. 2014, 566, A95. [Google Scholar] [CrossRef] [Green Version]
  7. Bladh, S.; Eriksson, K.; Marigo, P.; Liljegren, S.; Aringer, B. Carbon star wind models at solar and sub-solar metallicities: A comparative study. I. Mass loss and the properties of dust-driven winds. Astron. Astrophys. 2019, 623, A119. [Google Scholar] [CrossRef] [Green Version]
  8. Kalirai, J.S.; Marigo, P.; Tremblay, P.E. The Core Mass Growth and Stellar Lifetime of Thermally Pulsing Asymptotic Giant Branch Stars. Astrophys. J. 2014, 782, 17. [Google Scholar] [CrossRef] [Green Version]
  9. Straniero, O.; Domínguez, I.; Cristallo, S.; Gallino, R. Low-Mass AGB Stellar Models for 0.003 <= Z <= 0.02: Basic Formulae for Nucleosynthesis Calculations. Publ. Astron. Soc. Aust. 2003, 20, 389–392. [Google Scholar] [CrossRef] [Green Version]
  10. Cristallo, S.; Straniero, O.; Piersanti, L.; Gobrecht, D. Evolution, Nucleosynthesis, and Yields of AGB Stars at Different Metallicities. III. Intermediate-mass Models, Revised Low-mass Models, and the ph-FRUITY Interface. Astrophys. J. Suppl. 2015, 219, 40. [Google Scholar] [CrossRef]
  11. Marigo, P.; Aringer, B. Low-temperature gas opacity. ÆSOPUS: A versatile and quick computational tool. Astron. Astrophys. 2009, 508, 1539–1569. [Google Scholar] [CrossRef]
  12. Marigo, P. Asymptotic Giant Branch evolution at varying surface C/O ratio: Effects of changes in molecular opacities. Astron. Astrophys. 2002, 387, 507–519. [Google Scholar] [CrossRef]
  13. Marigo, P.; Bressan, A.; Nanni, A.; Girardi, L.; Pumo, M.L. Evolution of thermally pulsing asymptotic giant branch stars–I. The COLIBRI code. Month. Not. RAS 2013, 434, 488–526. [Google Scholar] [CrossRef]
  14. Trabucchi, M.; Wood, P.R.; Montalbán, J.; Marigo, P.; Pastorelli, G.; Girardi, L. Modelling long-period variables–I. A new grid of O-rich and C-rich pulsation models. Month. Not. RAS 2019, 482, 929–949. [Google Scholar] [CrossRef]
  15. Höfner, S.; Olofsson, H. Mass loss of stars on the asymptotic giant branch. Mechanisms, models and measurements. Astron. Astrophys. Rev. 2018, 26, 1. [Google Scholar] [CrossRef] [Green Version]
  16. Bladh, S.; Liljegren, S.; Höfner, S.; Aringer, B.; Marigo, P. An extensive grid of DARWIN models for M-type AGB stars. I. Mass-loss rates and other properties of dust-driven winds. Astron. Astrophys. 2019, 626, A100. [Google Scholar] [CrossRef] [Green Version]
  17. Ferrarotti, A.S.; Gail, H.P. Composition and quantities of dust produced by AGB-stars and returned to the interstellar medium. Astron. Astrophys. 2006, 447, 553–576. [Google Scholar] [CrossRef]
  18. Herwig, F. The evolution of AGB stars with convective overshoot. Astron. Astrophys. 2000, 360, 952–968. [Google Scholar]
  19. Wood, P.R. The conditions for dredge-up of carbon during the helium shell flash and the production of carbon stars. In Physical Processes in Red Giants; Astrophysics and Space Science Library; Iben, I.J., Renzini, A., Eds.; Springer: Dordrecht, The Netherlands, 1981; Volume 88, pp. 135–139. [Google Scholar] [CrossRef]
  20. Pastorelli, G.; Marigo, P.; Girardi, L.; Chen, Y.; Rubele, S.; Trabucchi, M.; Aringer, B.; Bladh, S.; Bressan, A.; Montalbán, J.; et al. Constraining the thermally pulsing asymptotic giant branch phase with resolved stellar populations in the Small Magellanic Cloud. Month. Not. RAS 2019, 485, 5666–5692. [Google Scholar] [CrossRef]
  21. Pastorelli, G.; Marigo, P.; Girardi, L.; Aringer, B.; Chen, Y.; Rubele, S.; Trabucchi, M.; Bladh, S.; Boyer, M.L.; Bressan, A.; et al. Constraining the thermally pulsing asymptotic giant branch phase with resolved stellar populations in the Large Magellanic Cloud. Month. Not. RAS 2020, 498, 3283–3301. [Google Scholar] [CrossRef]
  22. Seaton, M.J. Opacity Project data on CD for mean opacities and radiative accelerations. Month. Not. RAS 2005, 362, L1–L3. [Google Scholar] [CrossRef]
  23. Winters, J.M.; Le Bertre, T.; Jeong, K.S.; Helling, C.; Sedlmayr, E. A systematic investigation of the mass loss mechanism in dust forming long-period variable stars. Astron. Astrophys. 2000, 361, 641–659. [Google Scholar]
  24. Bowen, G.H. Dynamical modeling of long-period variable star atmospheres. Astrophys. J. 1988, 329, 299–317. [Google Scholar] [CrossRef]
  25. Bedijn, P.J. Pulsation, mass loss, and evolution of upper asymptotic giant branch stars. Astron. Astrophys. 1988, 205, 105–124. [Google Scholar]
  26. Bloecker, T. Stellar evolution of low and intermediate-mass stars. I. Mass loss on the AGB and its consequences for stellar evolution. Astron. Astrophys. 1995, 297, 727–738. [Google Scholar]
  27. Ventura, P.; D’Antona, F. The self-enrichment scenario in intermediate metallicity globular clusters. Astron. Astrophys. 2008, 479, 805–816. [Google Scholar] [CrossRef] [Green Version]
  28. Straniero, O.; Cristallo, S.; Piersanti, L. Heavy Elements in Globular Clusters: The Role of Asymptotic Giant Branch Stars. Astrophys. J. 2014, 785, 77. [Google Scholar] [CrossRef] [Green Version]
  29. Goriely, S.; Siess, L. S-process in hot AGB stars: A complex interplay between diffusive mixing and nuclear burning. Astron. Astrophys. 2004, 421, L25–L28. [Google Scholar] [CrossRef]
  30. Ventura, P.; Karakas, A.; Dell’Agli, F.; García-Hernández, D.A.; Guzman-Ramirez, L. Gas and dust from solar metallicity AGB stars. Month. Not. RAS 2018, 475, 2282–2305. [Google Scholar] [CrossRef]
  31. Karakas, A.I. Helium enrichment and carbon-star production in metal-rich populations. Month. Not. RAS 2014, 445, 347–358. [Google Scholar] [CrossRef] [Green Version]
  32. Karakas, A.I.; Lugaro, M. Stellar Yields from Metal-rich Asymptotic Giant Branch Models. Astrophys. J. 2016, 825, 26. [Google Scholar] [CrossRef]
  33. Busso, M.; Gallino, R.; Wasserburg, G.J. Nucleosynthesis in Asymptotic Giant Branch Stars: Relevance for Galactic Enrichment and Solar System Formation. Ann. Rev. Astron. Astrophys. 1999, 37, 239–309. [Google Scholar] [CrossRef] [Green Version]
  34. Cummings, J.D.; Kalirai, J.S.; Tremblay, P.E.; Ramirez-Ruiz, E.; Choi, J. The White Dwarf Initial-Final Mass Relation for Progenitor Stars from 0.85 to 7.5 M. Astrophys. J. 2018, 866, 21. [Google Scholar] [CrossRef] [Green Version]
  35. Miller Bertolami, M.M. New models for the evolution of post-asymptotic giant branch stars and central stars of planetary nebulae. Astron. Astrophys. 2016, 588, A25. [Google Scholar] [CrossRef]
  36. Cummings, J.D.; Kalirai, J.S.; Choi, J.; Georgy, C.; Tremblay, P.E.; Ramirez-Ruiz, E. A Novel Approach to Constrain Rotational Mixing and Convective-core Overshoot in Stars Using the Initial-Final Mass Relation. Astrophys. J. Lett. 2019, 871, L18. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Impact of carbon enrichment on the evolution of TP-AGB stars. Two models are shown as an example, namely M i = 1.8 M (left panels) and M i = 2.8 M (right panels), which reach approximately the same final mass, M f 0.70 0.71 M at the end of the TP-AGB phase. (a,b) Evolution on the H-R diagram during the whole TP-AGB evolution, color-coded according to the current photospheric C / O ratio. (c,d) Evolution of the mass-loss rate. Time is set to zero at the first TP. The tracks are color-coded according to the current photospheric C / O ratio. The drop in mass loss during the carbon-star phase occurs when the carbon excess, C–O, is below the current threshold to trigger a dust-driven wind. See the text for more details.
Figure 1. Impact of carbon enrichment on the evolution of TP-AGB stars. Two models are shown as an example, namely M i = 1.8 M (left panels) and M i = 2.8 M (right panels), which reach approximately the same final mass, M f 0.70 0.71 M at the end of the TP-AGB phase. (a,b) Evolution on the H-R diagram during the whole TP-AGB evolution, color-coded according to the current photospheric C / O ratio. (c,d) Evolution of the mass-loss rate. Time is set to zero at the first TP. The tracks are color-coded according to the current photospheric C / O ratio. The drop in mass loss during the carbon-star phase occurs when the carbon excess, C–O, is below the current threshold to trigger a dust-driven wind. See the text for more details.
Universe 08 00243 g001
Figure 2. Properties of the third dredge-up in TP-AGB stars at solar metallicity. The IFMR-based calibration of the dredged-up mass Δ M 3 DU (a), and efficiency λ (b), as a function of the current core mass. COLIBRI TP-AGB models with solar metallicity cover the initial mass range 1.4 M i / M 6.0 . (c) The maximum efficiency of the 3DU as a function of the initial stellar mass, predicted from various authors, namely: This work [magenta diamonds]; ([30], Ventura, purple stars); ([31,32], Monash, red circles); ([10], FRUITY, blue squares).
Figure 2. Properties of the third dredge-up in TP-AGB stars at solar metallicity. The IFMR-based calibration of the dredged-up mass Δ M 3 DU (a), and efficiency λ (b), as a function of the current core mass. COLIBRI TP-AGB models with solar metallicity cover the initial mass range 1.4 M i / M 6.0 . (c) The maximum efficiency of the 3DU as a function of the initial stellar mass, predicted from various authors, namely: This work [magenta diamonds]; ([30], Ventura, purple stars); ([31,32], Monash, red circles); ([10], FRUITY, blue squares).
Universe 08 00243 g002
Figure 4. Structural properties of the Z = 0.014 models with M i = 4.6 M (a), and M i = 4.8 M (b). Moving downward, the panels show the evolution during the whole TP-AGB phase of the temperature, T bce , at the bottom the convective envelope (top panel), the current mass being reduced by stellar winds, and the photospheric C/O (middle panel), the surface abundance (in molar fraction X / A , where X is the mass fraction abundance, and A denotes the mass number of the species) of a few elements (bottom panel). The horizontal dotted line corresponds to C / O = 1 . Note that the during the temporary C-rich stages the evolution proceeds at almost constant mass. See the text for more details.
Figure 4. Structural properties of the Z = 0.014 models with M i = 4.6 M (a), and M i = 4.8 M (b). Moving downward, the panels show the evolution during the whole TP-AGB phase of the temperature, T bce , at the bottom the convective envelope (top panel), the current mass being reduced by stellar winds, and the photospheric C/O (middle panel), the surface abundance (in molar fraction X / A , where X is the mass fraction abundance, and A denotes the mass number of the species) of a few elements (bottom panel). The horizontal dotted line corresponds to C / O = 1 . Note that the during the temporary C-rich stages the evolution proceeds at almost constant mass. See the text for more details.
Universe 08 00243 g004
Figure 5. Comparison of TP-AGB models at solar metallicity, from various sources. References and symbols are the same as in Figure 3. (a) TP-AGB phase lifetime as a function of the stellar initial mass. (b) 12C ejecta contributed by TP-AGB stars as a function of the stellar initial mass.
Figure 5. Comparison of TP-AGB models at solar metallicity, from various sources. References and symbols are the same as in Figure 3. (a) TP-AGB phase lifetime as a function of the stellar initial mass. (b) 12C ejecta contributed by TP-AGB stars as a function of the stellar initial mass.
Universe 08 00243 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Marigo, P. The Initial-Final Mass Relation of White Dwarfs: A Tool to Calibrate the Third Dredge-Up. Universe 2022, 8, 243. https://0-doi-org.brum.beds.ac.uk/10.3390/universe8040243

AMA Style

Marigo P. The Initial-Final Mass Relation of White Dwarfs: A Tool to Calibrate the Third Dredge-Up. Universe. 2022; 8(4):243. https://0-doi-org.brum.beds.ac.uk/10.3390/universe8040243

Chicago/Turabian Style

Marigo, Paola. 2022. "The Initial-Final Mass Relation of White Dwarfs: A Tool to Calibrate the Third Dredge-Up" Universe 8, no. 4: 243. https://0-doi-org.brum.beds.ac.uk/10.3390/universe8040243

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop