$Chempy$: A flexible chemical evolution model for abundance fitting - Do the Sun's abundances alone constrain chemical evolution models?

Abstract

Elemental abundances of stars are the result of the complex enrichment history of their galaxy. Interpretation of observed abundances requires flexible modeling tools to explore and quantify the information about Galactic chemical evolution (GCE) stored in such data. Here we present Chempy, a newly developed code for GCE modeling, representing a parametrized open one-zone model within a Bayesian framework. A Chempy model is specified by a set of 5-10 parameters that describe the effective galaxy evolution along with the stellar and star-formation physics: e.g. the star-formation history, the feedback efficiency, the stellar initial mass function (IMF) and the incidence of supernova type Ia (SN Ia). Unlike established approaches, Chempy can sample the posterior probability distribution in the full model parameter space and test data-model matches for different nucleosynthetic yield sets. We extend Chempy to a multi-zone scheme. As an illustrative application, we show that interesting parameter constraints result from only the ages and elemental abundances of Sun, Arcturus and the present-day interstellar medium (ISM). For the first time, we use such information to infer IMF parameter via GCE modeling, where we properly marginalize over nuisance parameters and account for different yield sets. We find that of the IMF $11.6_{-1.6}^{+2.1}$ % explodes as core-collapse SN, compatible with Salpeter 1955. We also constrain the incidence of SN Ia per 10^3 Msun to 0.5-1.4. At the same time, this Chempy application shows persistent discrepancies between predicted and observed abundances for some elements, irrespective of the chosen yield set. These cannot be remedied by any variations of Chempy's parameters and could be an indication for missing nucleosynthetic channels. Chempy should be a powerful tool to confront predictions from stellar nucleosynthesis with far more complex abundance data sets.

19 pages, 17 figures, accepted for publication in A&A, python code: https://github.com/jan-rybizki/Chempy

References in corpus (6)

• The NumPy array: a structure for efficient numerical computation
• Galactic chemical evolution: Carbon through Zinc
• Impact of radial migration on stellar and gas radial metallicity distribution
• Gaia FGK Benchmark stars: Opening the black box of stellar element abundance determination
• Towards a fully consistent Milky Way disc model - III. Constraining the initial mass function
• Fluorine in the solar neighborhood - is it all produced in AGB-stars?

Cited by in corpus (12)

• The Origin of r-Process Elements in the Milky Way
• An Empirical Measurement of the Initial-Final Mass Relation with Gaia White Dwarfs
• APOGEE DR16: a multi-zone chemical evolution model for the Galactic disc based on MCMC methods
• Spatial variations in the Milky Way disc metallicity-age relation
• Age-resolved chemistry of red giants in the solar neighbourhood
• Bimodality of [α/Fe]-[Fe/H] distributions is a natural outcome of dissipative collapse and disc growth in Milky Way-type galaxies
• Abundance ratios in GALAH DR2 and their implications for nucleosynthesis
• The chemical evolution of r-process elements from neutron star mergers: the role of a 2-phase interstellar medium
• The Milky Way's bulge star formation history as constrained from its bimodal chemical abundance distribution
• On the Optimal Choice of Nucleosynthetic Yields, IMF and Number of SNe Ia for Chemical Evolution Modelling
• BEDE: Bayesian Estimates of Dust Evolution For Nearby Galaxies
• Towards a fully consistent Milky Way disk model -- IV. The impact of Gaia DR2 and APOGEE