skip to main content


Title: The chemical evolution of the solar neighbourhood for planet-hosting stars
ABSTRACT

Theoretical physical-chemical models for the formation of planetary systems depend on data quality for the Sun’s composition, that of stars in the solar neighbourhood, and of the estimated ’pristine’ compositions for stellar systems. The effective scatter and the observational uncertainties of elements within a few hundred parsecs from the Sun, even for the most abundant metals like carbon, oxygen and silicon, are still controversial. Here we analyse the stellar production and the chemical evolution of key elements that underpin the formation of rocky (C, O, Mg, Si) and gas/ice giant planets (C, N, O, S). We calculate 198 galactic chemical evolution (GCE) models of the solar neighbourhood to analyse the impact of different sets of stellar yields, of the upper mass limit for massive stars contributing to GCE (Mup) and of supernovae from massive-star progenitors which do not eject the bulk of the iron-peak elements (faint supernovae). Even considering the GCE variation produced via different sets of stellar yields, the observed dispersion of elements reported for stars in the Milky Way (MW) disc is not reproduced. Among others, the observed range of super-solar [Mg/Si] ratios, sub-solar [S/N], and the dispersion of up to 0.5 dex for [S/Si] challenge our models. The impact of varying Mup depends on the adopted supernova yields. Thus, observations do not provide a constraint on the Mup parametrization. When including the impact of faint supernova models in GCE calculations, elemental ratios vary by up to 0.1–0.2 dex in the MW disc; this modification better reproduces observations.

 
more » « less
Award ID(s):
1927130
NSF-PAR ID:
10439948
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
524
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
p. 6295-6330
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. Older models of Galactic chemical evolution (GCE) predict [K/Fe] ratios as much as 1 dex lower than those inferred from stellar observations. Abundances of potassium are mainly based on analyses of the 7698 Å resonance line, and the discrepancy between GCE models and observations is in part caused by the assumption of local thermodynamic equilibrium (LTE) in spectroscopic analyses. Aims. We study the statistical equilibrium of K  I , focusing on the non-LTE effects on the 7698 Å line. We aim to determine how non-LTE abundances of potassium can improve the analysis of its chemical evolution, and help to constrain the yields of GCE models. Methods. We construct a new model K  I atom that employs the most up-to-date atomic data. In particular, we calculate and present inelastic e+K collisional excitation cross-sections from the convergent close-coupling (CCC) and the B -Spline R -matrix (BSR) methods, and H+K collisions from the two-electron model (LCAO). We constructed a fine, extended grid of non-LTE abundance corrections based on 1D MARCS models that span 4000 < T eff ∕K < 8000, 0.50 < log g < 5.00, − 5.00 < [Fe/H] < + 0.50, and applied the corrections to potassium abundances extracted from the literature. Results. In concordance with previous studies, we find severe non-LTE effects in the 7698 Å line. The line is stronger in non-LTE and the abundance corrections can reach approximately − 0.7 dex for solar-metallicity stars such as Procyon. We determine potassium abundances in six benchmark stars, and obtain consistent results from different optical lines. We explore the effects of atmospheric inhomogeneity by computing for the first time a full 3D non-LTE stellar spectrum of K  I lines for a test star. We find that 3D modeling is necessary to predict a correct shape of the resonance 7698 Å line, but the line strength is similar to that found in 1D non-LTE. Conclusions. Our non-LTE abundance corrections reduce the scatter and change the cosmic trends of literature potassium abundances. In the regime [Fe/H] ≲−1.0 the non-LTE abundances show a good agreement with the GCE model with yields from rotating massive stars. The reduced scatter of the non-LTE corrected abundances of a sample of solar twins shows that line-by-line differential analysis techniques cannot fully compensate for systematic LTE modelling errors; the scatter introduced by such errors introduces a spurious dispersion to K evolution. 
    more » « less
  2. Context. Phosphorus (P) is considered to be one of the key elements for life, making it an important element to look for in the abundance analysis of spectra of stellar systems. Yet, only a select number of spectroscopic studies exist to estimate the phosphorus abundances and investigate its trend across a range of metallicities. This is due to the lack of good phosphorus lines in the optical wavelength region and the requirement of careful manual analysis of the blended phosphorus lines in near-infrared H-band spectra obtained with individual observations and surveys such as the Apache Point Observatory Galactic Evolution Experiment (APOGEE). Aims. Based on a consistent and systematic analysis of high-resolution, near-infrared Immersion GRating INfrared Spectrograph (IGRINS) spectra of 38 K giant stars in the Solar neighborhood, we present and investigate the phosphorus abundance trend in the metallicity range of −1.2 dex < [Fe/H] < 0.4 dex. Furthermore, we compare this trend with the available chemical evolution models to shed some light on the origin and evolution of phosphorus. Methods. We have observed full H - and K -band spectra at a spectral resolving power of R = 45 000 with IGRINS mounted on the Gemini South telescope, the Discovery Channel Telescope, and the Harlan J Smith Telescope at McDonald Observatory. Abundances were determined from spectral lines by modeling the synthetic spectrum that best matches the observed spectrum by χ 2 minimization. For this task, we used the Spectroscopy Made Easy (SME) tool in combination with one-dimensional (1D) Model Atmospheres in a Radiative and Convective Scheme (MARCS) stellar atmosphere models. The investigated sample of stars have reliable stellar parameters estimated using optical FIber-fed Echelle Spectrograph (FIES) spectra obtained in a previous study of a set of stars called Giants in the Local Disk (GILD). In order to determine the phosphorus abundances from the 16482.92 Å phosphorus line, we needed to take special care blending the CO( v = 7−4) line. With the stellar parameters known, we thus determined the C, N, and O abundances from atomic carbon and a range of nonblended molecular lines (CO, CN, and OH) which are plentiful in the H-band region of K giant stars, assuring an appropriate modeling of the blending CO( v = 7−4) line. Results. We present the [P/Fe] versus [Fe/H] trend for K giant stars in the metallicity range of −1.2 dex < [Fe/H] < 0.4 dex and enhanced phosphorus abundances for two metal-poor s-rich stars. We find that our trend matches well with the compiled literature sample of prominently dwarf stars and the limited number of giant stars. Our trend is found to be higher by ~0.05−0.1 dex compared to the theoretical chemical evolution trend resulting from the core collapse supernova (type II) of massive stars with the phosphorus yields arbitrarily increased by a factor of 2.75. Thus the enhancement factor might need to be ~0.05−0.1 dex higher to match our trend. We also find an empirically determined primary behavior for phosphorus. Furthermore, the phosphorus abundance is found to be elevated by ~0.6−0.9 dex in the two s-enriched stars compared to the theoretical chemical evolution trend. 
    more » « less
  3. Abstract

    We demonstrate that using up to seven stellar abundance ratios can place observational constraints on the star formation histories (SFHs) of Local Group dSphs, using Sculptor dSph as a test case. We use a one-zone chemical evolution model to fit the overall abundance patterns ofαelements (which probe the core-collapse supernovae that occur shortly after star formation),s-process elements (which probe AGB nucleosynthesis at intermediate delay times), and iron-peak elements (which probe delayed Type Ia supernovae). Our best-fit model indicates that Sculptor dSph has an ancient SFH, consistent with previous estimates from deep photometry. However, we derive a total star formation duration of ∼0.9 Gyr, which is shorter than photometrically derived SFHs. We explore the effect of various model assumptions on our measurement and find that modifications to these assumptions still produce relatively short SFHs of duration ≲1.4 Gyr. Our model is also able to compare sets of predicted nucleosynthetic yields for supernovae and AGB stars, and can provide insight into the nucleosynthesis of individual elements in Sculptor dSph. We find that observed [Mn/Fe] and [Ni/Fe] trends are most consistent with sub-MChType Ia supernova models, and that a combination of “prompt” (delay times similar to core-collapse supernovae) and “delayed” (minimum delay times ≳50 Myr)r-process events may be required to reproduce observed [Ba/Mg] and [Eu/Mg] trends.

     
    more » « less
  4. null (Ed.)
    ABSTRACT Stellar and supernova nucleosynthesis in the first few billion years of the cosmic history have set the scene for early structure formation in the Universe, while little is known about their nature. Making use of stellar physical parameters measured by GALAH Data Release 3 with accurate astrometry from the Gaia EDR3, we have selected ∼100 old main-sequence turn-off stars (ages ≳12 Gyr) with kinematics compatible with the Milky Way stellar halo population in the Solar neighbourhood. Detailed homogeneous elemental abundance estimates by GALAH DR3 are compared with supernova yield models of Pop III (zero-metal) core-collapse supernovae (CCSNe), normal (non-zero-metal) CCSNe, and Type Ia supernovae (SN Ia) to examine which of the individual yields or their combinations best reproduce the observed elemental abundance patterns for each of the old halo stars (‘OHS’). We find that the observed abundances in the OHS with [Fe/H] > −1.5 are best explained by contributions from both CCSNe and SN Ia, where the fraction of SN Ia among all the metal-enriching SNe is up to 10–20 per cent for stars with high [Mg/Fe] ratios and up to 20–27 per cent for stars with low [Mg/Fe] ratios, depending on the assumption about the relative fraction of near-Chandrasekhar-mass SNe Ia progenitors. The results suggest that, in the progenitor systems of the OHS with [Fe/H] > −1.5, ∼ 50–60 per cent of Fe mass originated from normal CCSNe at the earliest phases of the Milky Way formation. These results provide an insight into the birth environments of the oldest stars in the Galactic halo. 
    more » « less
  5. ABSTRACT

    The absence of Type IIP core-collapse supernovae arising from progenitors above 17 solar masses suggests the existence of another evolutionary path by which massive stars end their lives. The direct collapse of a stellar core to a black hole without the production of a bright, explosive transient is expected to produce a long-lived, dim, red transient known as a failed supernova. Despite the detection of a number of candidates for disappearing massive stars in recent years, conclusive observational evidence for failed supernovae remains elusive. A custom-built pipeline designed for the detection of faint transients is used to re-analyse 10 yr of observations of 231 nearby galaxies from the PTF/ZTF surveys. This analysis recovers known supernovae, and yields a number of interesting transients. However, none of these are consistent with a failed supernova. Through Monte Carlo tests the recovery efficiency of our pipeline is quantified. By assuming failed supernovae occur as a Poissonian process with zero detections in the data set, 95 per cent upper limits to the rate of failed supernovae are calculated as a function of failed supernova absolute magnitude. We estimate failed supernovae to be less than 0.61, 0.33, 0.26, or 0.23 of the core-collapse SN rate for absolute magnitudes of −11, −12, −13, and −14, respectively. Finally, we show that if they exist, the Vera C. Rubin Observatory will find 1.7–3.7 failed SNe per year for an absolute bolometric luminosity of ∼6 × 1039 erg s−1 out to distances of 33–43 Mpc, depending on their assumed spectral energy distribution.

     
    more » « less