skip to main content

Title: Non-LTE analysis of K I in late-type stars
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
Award ID(s):
1834740 1803844
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Astronomy & Astrophysics
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Aims. An analysis of the methylidyne (CH) radical in non-local thermodynamic equilibrium (NLTE) is performed for the physical conditions of cool stellar atmospheres typical of red giants (log ɡ = 2.0, T eff = 4500 K) and the Sun. The aim of the present work is to explore whether the G band of the CH molecule, which is commonly used in abundance diagnostics of carbon-enhanced metal-poor stars, is sensitive to NLTE effects. Methods. LTE and NLTE theoretical spectra were computed with the MULTI code. We used one-dimensional (1D) LTE hydrostatic MARCS model atmospheres with parameters representing eleven red giant stars with metallicities ranging from [Fe/H] = −4.0 to [Fe/H] = 0.0 and carbon-to-iron ratios of [C/Fe] = 0.0, +0.7, +1.5, and +3.0. The CH molecule model was represented by 1981 energy levels, 18 377 radiative bound-bound transitions, and 932 photo-dissociation reactions. The rates due to transitions caused by collisions with free electrons and hydrogen atoms were computed using classical recipes. Results. Our calculations suggest that NLTE effects in the statistical equilibrium of the CH molecule are significant and cannot be neglected for precision spectroscopic analysis of C abundances. The NLTE effects are mostly driven by radiative over-dissociation, owing to the very low dissociation threshold of the molecule and significant resonances in the photo-dissociation cross-sections. The NLTE effects in the G band increase with decreasing metallicity. When comparing the C abundances determined from the CH G band in LTE and in NLTE, we show that the C abundances are always under-estimated if LTE is assumed. The NLTE corrections to C abundance inferred from the CH feature range from +0.04 dex for the Sun to +0.21 dex for a red giant with metallicity [Fe/H] = −4.0. Conclusions. Departures from the LTE assumption in the CH molecule are non-negligible, and NLTE effects have to be taken into account in the diagnostic spectroscopy based on the CH lines. We show here that the NLTE effects in the optical CH lines are non-negligible for the Sun and red giant stars, but further calculations are warranted to investigate the effects in other regimes of stellar parameters. 
    more » « less
  3. Abstract Individual chemical abundances for 14 elements (C, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni) are derived for a sample of M dwarfs using high-resolution, near-infrared H -band spectra from the Sloan Digital Sky Survey-IV/Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. The quantitative analysis included synthetic spectra computed with 1D LTE plane-parallel MARCS models using the APOGEE Data Release 17 line list to determine chemical abundances. The sample consists of 11 M dwarfs in binary systems with warmer FGK dwarf primaries and 10 measured interferometric angular diameters. To minimize atomic diffusion effects, [X/Fe] ratios are used to compare M dwarfs in binary systems and literature results for their warmer primary stars, indicating good agreement (<0.08 dex) for all studied elements. The mean abundance difference in primaries minus this work’s M dwarfs is −0.05 ± 0.03 dex. It indicates that M dwarfs in binary systems are a reliable way to calibrate empirical relationships. A comparison with abundance, effective temperature, and surface gravity results from the APOGEE Stellar Parameter and Chemical Abundances Pipeline (ASPCAP) Data Release 16 finds a systematic offset of [M/H], T eff , log g = +0.21 dex, −50 K, and 0.30 dex, respectively, although ASPCAP [X/Fe] ratios are generally consistent with this study. The metallicities of the M dwarfs cover the range of [Fe/H] = −0.9 to +0.4 and are used to investigate Galactic chemical evolution via trends of [X/Fe] as a function of [Fe/H]. The behavior of the various elemental abundances [X/Fe] versus [Fe/H] agrees well with the corresponding trends derived from warmer FGK dwarfs, demonstrating that the APOGEE spectra can be used to examine Galactic chemical evolution using large samples of selected M dwarfs. 
    more » « less
  4. Abstract Precise fundamental atmospheric stellar parameters and abundance determination of individual elements in stars are important for all stellar population studies. Non–local thermodynamic equilibrium (non-LTE; hereafter NLTE) models are often important for such high precision, however, can be computationally complex and expensive, which renders the models less utilized in spectroscopic analyses. To alleviate the computational burden of such models, we developed a robust 1D, NLTE fundamental atmospheric stellar parameter derivation tool, LOTUS , to determine the effective temperature T eff , surface gravity log g , metallicity [Fe/H], and microturbulent velocity v mic for FGK-type stars, from equivalent width (EW) measurements of Fe i and Fe ii lines. We utilize a generalized curve of growth method to take into account the EW dependencies of each Fe i and Fe ii line on the corresponding atmospheric stellar parameters. A global differential evolution optimization algorithm is then used to derive the fundamental parameters. Additionally, LOTUS can determine precise uncertainties for each stellar parameter using a Markov Chain Monte Carlo algorithm. We test and apply LOTUS on a sample of benchmark stars, as well as stars with available asteroseismic surface gravities from the K2 survey, and metal-poor stars from the Gaia-ESO and R -Process Alliance surveys. We find very good agreement between our NLTE-derived parameters in LOTUS to nonspectroscopic values on average within T eff = ±30 K, and log g = ±0.10 dex for benchmark stars. We provide open access of our code, as well as of the interpolated precomputed NLTE EW grids available on Github (the software is available on GitHub 3 3 under an MIT License, and version 0.1.1 (as the persistent version) is archived in Zenodo) and documentation with working examples on the Readthedocs book. 
    more » « less
  5. null (Ed.)
    Context. Rubidium is one of the few elements produced by the neutron capture s - and r -processes in almost equal proportions. Recently, a Rb deficiency ([Rb/Fe] < 0.0), amounting to a factor of about two with respect to the Sun, has been found in M dwarfs of near-solar metallicity. This stands in contrast to the close-to-solar [Sr, Zr/Fe] ratios derived in the same stars. This deficiency is difficult to understand from the point of view of observations and of nucleosynthesis. Aims. To test the reliability of this Rb deficiency, we study the Rb and Zr abundances in a sample of KM-type giant stars across a similar metallicity range, extracted from the AMBRE Project. Methods. We used high-resolution and high signal-to-noise spectra to derive Rb and Zr abundances in a sample of 54 bright giant stars with metallicities in the range of −0.6 ≲ [Fe/H] ≲ +0.4 dex, via spectral synthesis in both local and non-local thermodynamic equilibrium (LTE and NLTE, respectively). We also studied the impact of the Zeeman broadening in the profile of the Rb  I at λ 7800 Å line. Results. The LTE analysis also results in a Rb deficiency in giant stars, however, it is considerably lower than that obtained in M dwarfs. However, once NLTE corrections are performed, the [Rb/Fe] ratios are very close to solar (average −0.01 ± 0.09 dex) in the full metallicity range studied here. This stands in contrast to the value found for M dwarfs. The [Zr/Fe] ratios derived are in excellent agreement with those obtained in previous studies in FGK dwarf stars with a similar metallicity. We investigate the effect of gravitational settling and magnetic activity as possible causes of the Rb deficiency found in M dwarfs. Although the former phenomenon has a negligible impact on the surface Rb abundance, the presence of an average magnetic field with an intensity that is typical of that observed in M dwarfs may result in systematic Rb abundance underestimations if the Zeeman broadening is not considered in the spectral synthesis. This may explain the Rb deficiency in M dwarfs, but not fully. On the other hand, the new [Rb/Fe] and [Rb/Zr] versus [Fe/H] relationships can be explained when the Rb production by rotating massive stars and low-to-intermediate mass stars (these latter also producing Zr) are considered, without the need to deviate from the standard s -process nucleosynthesis in asymptotic giant branch stars, as suggested previously. 
    more » « less