skip to main content


Title: LOTUS: A (Non-) LTE Optimization Tool for Uniform Derivation of Stellar Atmospheric Parameters
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 https://github.com/Li-Yangyang/LOTUS 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
Award ID(s):
2206263
NSF-PAR ID:
10435845
Author(s) / Creator(s):
;
Date Published:
Journal Name:
The Astronomical Journal
Volume:
165
Issue:
4
ISSN:
0004-6256
Page Range / eLocation ID:
145
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. Large spectroscopic surveys of the Milky Way must be calibrated against a sample of benchmark stars to ensure the reliable determination of atmospheric parameters. Aims. Here, we present new fundamental stellar parameters of seven giant and subgiant stars that will serve as benchmark stars for large surveys. The aim is to reach a precision of 1% in the effective temperature. This precision is essential for accurate determinations of the full set of fundamental parameters and abundances for stars observed by the stellar surveys. Methods. We observed HD 121370 ( η Boo), HD 161797 ( μ Her), HD 175955, HD 182736, HD 185351, HD 188512 ( β Aql), and HD 189349, using the high angular resolution optical interferometric instrument PAVO at the CHARA Array. The limb-darkening corrections were determined from 3D model atmospheres based on the STAGGER grid. The T eff were determined directly from the Stefan-Boltzmann relation, with an iterative procedure to interpolate over tables of bolometric corrections. We estimated surface gravities from comparisons to Dartmouth stellar evolution model tracks. The spectroscopic observations were collected from the ELODIE and FIES spectrographs. We estimated metallicities ([Fe/H]) from a 1D non-local thermodynamic equilibrium (NLTE) abundance analysis of unblended lines of neutral and singly ionised iron. Results. For six of the seven stars, we measured the value of T eff to better than 1% accuracy. For one star, HD 189349, the uncertainty on T eff is 2%, due to an uncertain bolometric flux. We do not recommend this star as a benchmark until this measurement can be improved. Median uncertainties for all stars in log  g and [Fe/H] are 0.034 dex and 0.07 dex, respectively. Conclusions. This study presents updated fundamental stellar parameters of seven giant and subgiant stars that can be used as a new set of benchmarks. All the fundamental stellar parameters were established on the basis of consistent combinations of interferometric observations, 3D limb-darkening modelling, and spectroscopic analysis. This paper in this series follows our previous papers featuring dwarf stars and stars in the metal-poor range. 
    more » « less
  2. Context. Stellar models applied to large stellar surveys of the Milky Way need to be properly tested against a sample of stars with highly reliable fundamental stellar parameters. We have established a programme aiming to deliver such a sample of stars. Aims. Here we present new fundamental stellar parameters of nine dwarf stars that will be used as benchmark stars for large stellar surveys. One of these stars is the solar-twin 18 Sco, which is also one of the Gaia -ESO benchmarks. The goal is to reach a precision of 1% in effective temperature ( T eff ). This precision is important for accurate determinations of the full set of fundamental parameters and abundances of stars observed by the surveys. Methods. We observed HD 131156 ( ξ Boo), HD 146233 (18 Sco), HD 152391, HD 173701, HD 185395 ( θ Cyg), HD 186408 (16 Cyg A), HD 186427 (16 Cyg B), HD 190360, and HD 207978 (15 Peg) using the high angular resolution optical interferometric instrument PAVO at the CHARA Array. We derived limb-darkening corrections from 3D model atmospheres and determined T eff directly from the Stefan–Boltzmann relation, with an iterative procedure to interpolate over tables of bolometric corrections. Surface gravities were estimated from comparisons to Dartmouth stellar evolution model tracks. We collected spectroscopic observations from the ELODIE spectrograph and estimated metallicities ([Fe/H]) from a 1D non-local thermodynamic equilibrium (NLTE) abundance analysis of unblended lines of neutral and singly ionised iron. Results. For eight of the nine stars we measure the T eff ⪅ 1%, and for one star better than 2%. We determined the median uncertainties in log  g and [Fe/H] as 0.015 dex and 0.05 dex, respectively. Conclusions. This study presents updated fundamental stellar parameters of nine dwarf stars that can be used as a new set of benchmarks. All the fundamental stellar parameters were based on consistently combining interferometric observations, 3D limb-darkening modelling, and spectroscopic analysis. The next paper in this series will extend our sample to giants in the metal-rich range. 
    more » « less
  3. Context. Benchmark stars are crucial as validating standards for current as well as future large stellar surveys of the Milky Way. However, the number of suitable metal-poor benchmark stars is currently limited, owing to the difficulty in determining reliable effective temperatures ( T eff ) in this regime. Aims. We aim to construct a new set of metal-poor benchmark stars based on reliable interferometric effective temperature determinations and a homogeneous analysis. The aim is to reach a precision of 1% in T eff , as is crucial for sufficiently accurate determinations of the full set of fundamental parameters and abundances for the survey sources. Methods. We observed ten late-type metal-poor dwarfs and giants: HD 2665, HD 6755, HD 6833, HD 103095, HD 122563, HD 127243, HD 140283, HD 175305, HD 221170, and HD 224930. Only three of them (HD 103095, HD 122563, and HD 140283) have previously been used as benchmark stars. For the observations, we used the high-angular-resolution optical interferometric instrument PAVO at the CHARA array. We modelled angular diameters using 3D limb-darkening models and determined effective temperatures directly from the Stefan-Boltzmann relation, with an iterative procedure to interpolate over tables of bolometric corrections. Surface gravities (log( g )) were estimated from comparisons to Dartmouth stellar evolution model tracks. We collected spectroscopic observations from the ELODIE and FIES spectrographs and estimated metallicities ([Fe/H]) from a 1D non-local thermodynamic equilibrium (NLTE) abundance analysis of unblended lines of neutral and singly ionised iron. Results. We inferred T eff to better than 1% for five of the stars (HD 103095, HD 122563, HD 127243, HD 140283, and HD 224930). The effective temperatures of the other five stars are reliable to between 2 and 3%; the higher uncertainty on the T eff for those stars is mainly due to their having a larger uncertainty in the bolometric fluxes. We also determined log( g ) and [Fe/H] with median uncertainties of 0.03 dex and 0.09 dex, respectively. Conclusions. This study presents reliable and homogeneous fundamental stellar parameters for ten metal-poor stars that can be adopted as a new set of benchmarks. The parameters are based on our consistent approach of combining interferometric observations, 3D limb-darkening-modelling and spectroscopic observations. The next paper in this series will extend this approach to dwarfs and giants in the metal-rich regime. 
    more » « less
  4. Physically realistic models of stellar spectra are needed in a variety of astronomical studies, from the analysis of fundamental stellar parameters, to studies of exoplanets and stellar populations in galaxies. Here we present a new version of the widely used radiative transfer code Turbospectrum, which we update so that it is able to perform spectrum synthesis for lines of multiple chemical elements in non-local thermodynamic equilibrium (NLTE). We use the code in the analysis of metallicites and abundances of the Gaia FGK benchmark stars, using 1D MARCS atmospheric models and the averages of 3D radiation-hydrodynamics simulations of stellar surface convection. We show that the new more physically realistic models offer a better description of the observed data, and we make the program and the associated microphysics data publicly available, including grids of NLTE departure coefficients for H, O, Na, Mg, Si, Ca, Ti, Mn, Fe, Co, Ni, Sr, and Ba. 
    more » « less
  5. 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