Abstract
The initialâfinal mass relation (IFMR) links the birth mass of a star to the mass of the compact remnant left at its death. While the relevance of the IFMR across astrophysics is universally acknowledged, not all of its fine details have yet been resolved. A new analysis of a few carbonâoxygen white dwarfs in old open clusters of the Milky Way led us to identify a kink in the IFMR, located over a range of initial masses, 1.65ââ²âMi/Mâââ²â2.10. The kinkâs peak in white dwarf mass of about 0.70â0.75âMâ is produced by stars with Miâââ1.8â1.9âMâ, corresponding to ages of about 1.8â1.7âGyr. Interestingly, this peak coincides with the initial mass limit between low-mass stars that develop a degenerate helium core after central hydrogen exhaustion, and intermediate-mass stars that avoid electron degeneracy. We interpret the IFMR kink as the signature of carbon star formation in the Milky Way. This finding is critical to constraining the evolution and chemical enrichment of low-mass stars, and their impact on the spectrophotometric properties of galaxies.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 /Â 30Â days
cancel any time
Subscribe to this journal
Receive 12 digital issues and online access to articles
$119.00 per year
only $9.92 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request. Montreal WD cooling models are publicly available from http://www.astro.umontreal.ca/~bergeron/CoolingModels. The pulsation periods are computed with fitting relations based on publicly available models that can be found at http://starkey.astro.unipd.it/pulsation_models.html.
Code availability
The stellar evolution codes PARSEC and COLIBRI are not publicly available. The mass-loss routine for carbon stars can be found at https://www.astro.uu.se/coolstars/TOOLS/MLR-routines/C/. The code to compute the dust-grain growth in the outflows of AGB stars can be retrieved from http://www.ita.uni-heidelberg.de/~gail/agbdust/agbdust.html. The code used here to calculate photometry-based WD parameters is available from https://github.com/SihaoCheng/WD_models.
Change history
14 July 2020
In the version of this Article originally published, the Montreal WD cooling model link in the Data availability statement and the AGB stars link in the Code availability statement were incorrect. They have now been updated.
References
Cescutti, G., Matteucci, F., McWilliam, A. & Chiappini, C. The evolution of carbon and oxygen in the bulge and disk of the Milky Way. Astron. Astrophys. 505, 605â612 (2009).
Gustafsson, B., Karlsson, T., Olsson, E., Edvardsson, B. & Ryde, N. The origin of carbon, investigated by spectral analysis of solar-type stars in the Galactic Disk. Astron. Astrophys. 342, 426â439 (1999).
Mattsson, L. The origin of carbon: low-mass stars and an evolving, initially top-heavy IMF? Astron. Astrophys. 515, A68 (2010).
Bensby, T. & Feltzing, S. The origin and chemical evolution of carbon in the galactic thin and thick discs. Mon. Not. R. Astron. Soc. 367, 1181â1193 (2006).
Karakas, A. I. & Lattanzio, J. C. The Dawes Review 2: nucleosynthesis and stellar yields of low- and intermediate-mass single stars. Publ. Astron. Soc. Aus. 31, e030 (2014).
Herwig, F. Evolution of asymptotic giant branch stars. Annu. Rev. Astron. Astrophys. 43, 435â479 (2005).
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. 692, 1013â1032 (2009).
Kalirai, J. S., Marigo, P. & Tremblay, P.-E. The core mass growth and stellar lifetime of thermally pulsing asymptotic giant branch stars. Astrophys. J. 782, 17 (2014).
Bird, J. C. & Pinsonneault, M. H. A bound on the light emitted during the thermally pulsing asymptotic giant branch phase. Astrophys. J. 733, 81 (2011).
Marigo, P. & Girardi, L. Coupling emitted light and chemical yields from stars: a basic constraint to population synthesis models of galaxies. Astron. Astrophys. 377, 132â147 (2001).
Weidemann, V. Revision of the initial-to-final mass relation. Astron. Astrophys. 363, 647â656 (2000).
Kalirai, J. S. et al. The initialâfinal mass relation: direct constraints at the low-mass end. Astrophys. J. 676, 594â609 (2008).
Cummings, J. D., Kalirai, J. S., Tremblay, P.-E. & Ramirez-Ruiz, E. Two massive white dwarfs from NGC 2323 and the initialâfinal mass relation for progenitors of 4 to 6.5 Mâ. Astrophys. J. 818, 84 (2016).
Kalirai, J. S. et al. The masses of population II white dwarfs. Astrophys. J. 705, 408â425 (2009).
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. 866, 21 (2018).
Williams, K. A. et al. Ensemble properties of the white dwarf population of the old, solar metallicity open star cluster Messier 67. Astrophys. J. 867, 62 (2018).
Canton, P. The InitialâFinal Mass Relation Revisited: A Monte Carlo Approach with the Addition of the M67 White Dwarf Population. PhD thesis, Univ. Oklahoma (2018).
Tremblay, P.-E. et al. The field white dwarf mass distribution. Mon. Not. R. Astron. Soc. 461, 2100â2114 (2016).
El-Badry, K., Rix, H.-W. & Weisz, D. R. An empirical measurement of the initialâfinal mass relation with Gaia white dwarfs. Astrophys. J. Lett. 860, L17 (2018).
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. 623, A119 (2019).
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. 566, A95 (2014).
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. 509, A14 (2010).
Marigo, P. & Aringer, B. Low-temperature gas opacity. ÃSOPUS: a versatile and quick computational tool. Astron. Astrophys. 508, 1539â1569 (2009).
Ferrarotti, A. S. & Gail, H. P. Mineral formation in stellar winds. III. Dust formation in S stars. Astron. Astrophys. 382, 256â281 (2002).
Ferrarotti, A. S. & Gail, H.-P. Composition and quantities of dust produced by AGB-stars and returned to the interstellar medium. Astron. Astrophys. 447, 553â576 (2006).
DellâAgli, F. et al. Asymptotic giant branch and super-asymptotic giant branch stars: modelling dust production at solar metallicity. Mon. Not. R. Astron. Soc. 467, 4431â4440 (2017).
Nanni, A., Bressan, A., Marigo, P. & Girardi, L. Evolution of thermally pulsing asymptotic giant branch starsâII. Dust production at varying metallicity. Mon. Not. R. Astron. Soc. 434, 2390â2417 (2013).
Höfner, S. & Olofsson, H. Mass loss of stars on the asymptotic giant branch. Mechanisms, models and measurements. Astron. Astrophys. Rev. 26, 1 (2018).
Schöier, F. L. & Olofsson, H. Models of circumstellar molecular radio line emission. Mass loss rates for a sample of bright carbon stars. Astron. Astrophys. 368, 969â993 (2001).
McDonald, I., De Beck, E., Zijlstra, A. A. & Lagadec, E. Pulsation-triggered dust production by asymptotic giant branch stars. Mon. Not. R. Astron. Soc. 481, 4984â4999 (2018).
McDonald, I. & Trabucchi, M. The onset of the AGB wind tied to a transition between sequences in the period-luminosity diagram. Mon. Not. R. Astron. Soc. 484, 4678â4682 (2019).
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. 361, 641â659 (2000).
Cummings, J. D. et al. A novel approach to constrain rotational mixing and convective-core overshoot in stars using the initialâfinal mass relation. Astrophys. J. 871, L18 (2019).
Marigo, P., Bressan, A., Nanni, A., Girardi, L. & Pumo, M. L. Evolution of thermally pulsing asymptotic giant branch starsâI. The COLIBRI code. Mon. Not. R. Astron. Soc. 434, 488â526 (2013).
Wagenhuber, J. & Groenewegen, M. A. T. New input data for synthetic AGB evolution. Astron. Astrophys. 340, 183â195 (1998).
Ventura, P., Karakas, A., DellâAgli, F., GarcÃa-Hernández, D. A. & Guzman-Ramirez, L. Gas and dust from solar metallicity AGB stars. Mon. Not. R. Astron. Soc. 475, 2282â2305 (2018).
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. 219, 40 (2015).
Bloecker, T. Stellar evolution of low and intermediate-mass stars. I. Mass loss on the AGB and its consequences for stellar evolution. Astron. Astrophys. 297, 727â738 (1995).
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. 626, A100 (2019).
Girardi, L., Marigo, P., Bressan, A. & Rosenfield, P. The insidious boosting of thermally pulsing asymptotic giant branch stars in intermediate-age magellanic cloud clusters. Astrophys. J. 777, 142 (2013).
Maraston, C. et al. Evidence for TP-AGB stars in high-redshift galaxies, and their effect on deriving stellar population parameters. Astrophys. J. 652, 85â96 (2006).
Bruzual, A. & Charlot, G. Spectral evolution of stellar populations using isochrone synthesis. Astrophys. J. 405, 538â553 (1993).
Gaia Collaboration et al. The Gaia mission. Astron. Astrophys. 595, A1 (2016).
Gaia Collaboration et al. Gaia Data Release 2. Summary of the contents and survey properties. Astron. Astrophys. 616, A1 (2018).
Curtis, J. L., Wolfgang, A., Wright, J. T., Brewer, J. M. & Johnson, J. A. Ruprecht 147: the oldest nearby open cluster as a new benchmark for stellar astrophysics. Astron. J. 145, 134 (2013).
Tremblay, P.-E., Bergeron, P. & Gianninas, A. An improved spectroscopic analysis of DA white dwarfs from the Sloan Digital Sky Survey Data Release 4. Astrophys. J. 730, 128 (2011).
Bergeron, P. et al. A comprehensive spectroscopic analysis of DB white dwarfs. Astrophys. J. 737, 28 (2011).
Genest-Beaulieu, C. & Bergeron, P. A comprehensive spectroscopic and photometric analysis of DA and DB white dwarfs from SDSS and Gaia. Astrophys. J. 871, 169 (2019).
Fontaine, G., Brassard, P. & Bergeron, P. The potential of white dwarf cosmochronology. Publ. Astron. Soc. Pac. 113, 409â435 (2001).
Cukanovaite, E., Tremblay, P.-E., Freytag, B., Ludwig, H.-G. & Bergeron, P. Pure-helium 3D model atmospheres of white dwarfs. Mon. Not. R. Astron. Soc. 481, 1522â1537 (2018).
Tremblay, P. E., Cukanovaite, E., Gentile Fusillo, N. P., Cunningham, T. & Hollands, M. A. Fundamental parameter accuracy of DA and DB white dwarfs in Gaia data release 2. Mon. Not. R. Astron. Soc. 482, 5222â5232 (2019).
Cummings, J. D. & Kalirai, J. S. Improved main-sequence turnoff ages of young open clusters: multicolor UBV techniques and the challenges of rotation. Astron. J. 156, 165 (2018).
Bressan, A. et al. PARSEC: stellar tracks and isochrones with the PAdova and TRieste stellar evolution code. Mon. Not. R. Astron. Soc. 427, 127â145 (2012).
Marigo, P. et al. A new generation of PARSEC-COLIBRI stellar isochrones including the TP-AGB phase. Astrophys. J. 835, 77 (2017).
Asplund, M., Grevesse, N., Sauval, A. J. & Scott, P. The chemical composition of the Sun. Annu. Rev. Astron. Astrophys. 47, 481â522 (2009).
Reimers, D. Circumstellar absorption lines and mass loss from red giants. Mem. Soc. R. Sci. Liege 8, 369â382 (1975).
Pastorelli, G. et al. Constraining the thermally pulsing asymptotic giant branch phase with resolved stellar populations in the small magellanic cloud. Mon. Not. R. Astron. Soc. 485, 5666â5692 (2019).
Cranmer, S. R. & Saar, S. H. Testing a predictive theoretical model for the mass loss rates of cool stars. Astrophys. J. 741, 54 (2011).
Bowen, G. H. Dynamical modeling of long-period variable star atmospheres. Astrophys. J. 329, 299â317 (1988).
Bedijn, P. J. Pulsation, mass loss, and evolution of upper asymptotic giant branch stars. Astron. Astrophys. 205, 105â124 (1988).
Vassiliadis, E. & Wood, P. R. Evolution of low- and intermediate-mass stars to the end of the asymptotic giant branch with mass loss. Astrophys. J. 413, 641â657 (1993).
Lambert, D. L., Gustafsson, B., Eriksson, K. & Hinkle, K. H. The chemical composition of carbon stars. IâCarbon, nitrogen, and oxygen in 30 cool carbon stars in the galactic disk. Astrophys. J. Suppl. 62, 373â425 (1986).
Acknowledgements
P.M., S.B., Y.C., L.G., G.P., M.T. and B.A. acknowledge the support from the ERC Consolidator Grant funding scheme (project STARKEY, grant agreement number 615604). P.-E.T. has received ERC funding under the European Unionâs Horizon 2020 research and innovation programme (grant agreement number 677706 â WD3D).
Author information
Authors and Affiliations
Contributions
P.M. designed and performed the theoretical research, ran the TP-AGB models and the population synthesis simulations, and provided the interpretation of the new IFMR data in terms of stellar evolution. J.D.C. performed the Keck observations, processed the data, analysed the cluster parameters, spectroscopically analysed the DA WDs and determined memberships. J.L.C. identified the likely WD candidates for observations and assisted with the cluster-parameter analysis. J.K. coordinated the observational and theoretical work and provided expertise. P.-E.T. provided the DA WD atmospheric models and fitting program and his expertise. E.R.-R. assisted with Keck observations. P.B. provided the DB WD atmospheric models and fit the DB parameters. S.B. provided expertise and help in implementing the mass-loss grid of dynamical atmospheres for carbon stars in the COLIBRI code. Y.C., A.B., L.G., G.P. and M.T. contributed to the development of the stellar models and the discussion of the results. S.C. contributed his WD photometric analysis expertise and his publicly available Python 3 module was used for the photometric-based derivation of WD parameters. B.A. provided expertise and the molecular opacity data to model the atmospheres of carbon stars. P.D.T. implemented the WD models in the populations synthesis simulations.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Astronomy thanks Krzysztof Gesicki, Iain McDonald and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Comparison between the semi-empirical IFMR and model results.
The semi-empirical data are shown with diamonds and error bars covering the range of ± 1 Ï. Newly discovered and newly analysed WD data (see Table 1) are shown in green. a-b, Predictions for the whole (Mi, λ) grid of models. c-d, Selected models that are found to match the semi-empirical IFMR. The theoretical IFMR is colour-coded according to the values of the efficiency of the 3DU (a-c) and the photospheric C/O at the end of the TP-AGB phase (b-d).
Extended Data Fig. 2 Examples of theoretical IFMRs that fail to account for the kink in the semi-empirical IFMR.
a, Too high efficiency of the 3DU in low-mass stars: λ = 0.5 is assumed for all models that experience the 3DU. b, Mass loss insensitive to the photospheric chemical composition: the B95 mass-loss formula is applied to all models, irrespective of the photospheric C/O. The semi-empirical IFMR is the same as in Fig. 1, with error bars covering the range of ± 1 Ï.
Supplementary information
Supplementary Information
Supplementary Figs. 1â6 and sections âThe WD mass distributionâ and âOther supporting evidence: Galactic semi-regular variablesâ.
Rights and permissions
About this article
Cite this article
Marigo, P., Cummings, J.D., Curtis, J.L. et al. Carbon star formation as seen through the non-monotonic initialâfinal mass relation. Nat Astron 4, 1102â1110 (2020). https://doi.org/10.1038/s41550-020-1132-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41550-020-1132-1
This article is cited by
-
Alkali metals in white dwarf atmospheres as tracers of ancient planetary crusts
Nature Astronomy (2021)