Showing 1–50 of 86 results for author: Barklem, P S

3D NLTE modelling of Y and Eu. Centre-to-limb variation and solar abundances

Authors: N. Storm, P. S. Barklem, S. A. Yakovleva, A. K. Belyaev, P. Palmeri, P. Quinet, K. Lodders, M. Bergemann, R. Hoppe

Abstract: Context. Abundances of s- and r-process elements in Sun-like stars constrain nucleosynthesis in extreme astrophysical events, such as compact binary mergers and explosions of highly magnetised rapidly rotating massive stars. Aims. We measure solar abundances of yttrium (Y) and europium (Eu) using 3D non-local thermal equilibrium (NLTE) models. We use the model to determine the abundance of Y, an… ▽ More Context. Abundances of s- and r-process elements in Sun-like stars constrain nucleosynthesis in extreme astrophysical events, such as compact binary mergers and explosions of highly magnetised rapidly rotating massive stars. Aims. We measure solar abundances of yttrium (Y) and europium (Eu) using 3D non-local thermal equilibrium (NLTE) models. We use the model to determine the abundance of Y, and also explore the model's ability to reproduce the solar centre-to-limb variation of its lines. In addition, we determine the Eu abundance using solar disc-centre and integrated flux spectra. Methods. We developed an NLTE model of Eu and updated our model of Y with collisional data from detailed quantum-mechanical calculations. We used the IAG spatially resolved high-resolution solar spectra to derive the solar abundances of Y across the solar disc and of Eu for integrated flux and at disc centre using a set of carefully selected lines and a 3D radiation-hydrodynamics model of the solar atmosphere. Results. We find 3D NLTE solar abundances of A(Y)$_{\textrm{3D NLTE}}$=$2.30 \pm 0.03_{\textrm{stat}} \pm 0.07_{\textrm{syst}}$ dex based on observations at all angles and A(Eu)$_{\textrm{3D NLTE}}$=$0.57 \pm 0.01_{\textrm{stat}} \pm 0.06_{\textrm{syst}}$ dex based on the integrated flux and disc-centre intensity. 3D NLTE modelling offers the most consistent abundances across the solar disc, and resolves the problem of severe systematic bias in Y and Eu abundances inherent to 1D LTE, 1D NLTE, and 3D LTE modelling. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 8 pages, 4 figures, accepted by A&A

Detailed $α$ abundance trends in the inner Galactic bulge

Authors: N. Nieuwmunster, G. Nandakumar, E. Spitoni, N. Ryde, M. Schultheis, R. M. Rich, P. S. Barklem, O. Agertz, F. Renaud, F. Matteucci

Abstract: In this paper, we aim to derive high-precision alpha-element abundances using CRIRES high-resolution IR spectra of 72 cool M giants of the inner Galactic bulge. Silicon, magnesium, and calcium abundances were determined by fitting a synthetic spectrum for each star. We also incorporated recent theoretical data into our spectroscopic analysis (i.e. updated K-band line list, better broadening parame… ▽ More In this paper, we aim to derive high-precision alpha-element abundances using CRIRES high-resolution IR spectra of 72 cool M giants of the inner Galactic bulge. Silicon, magnesium, and calcium abundances were determined by fitting a synthetic spectrum for each star. We also incorporated recent theoretical data into our spectroscopic analysis (i.e. updated K-band line list, better broadening parameter estimation, non-local thermodynamic equilibrium (NLTE) corrections). We compare these inner bulge alpha abundance trends with those of solar neighbourhood stars observed with IGRINS using the same line list and analysis technique; we also compare our sample to APOGEE DR17 abundances for inner bulge stars. We investigate bulge membership using spectro-photometric distances and orbital simulations. We construct a chemical-evolution model that fits our metallicity distribution function (MDF) and our alpha-element trends. Among our 72 stars, we find four that are not bulge members. [Si/Fe] and [Mg/Fe] versus [Fe/H] trends show a typical thick disc alpha-element behaviour, except that we do not see any plateau at supersolar metallicities as seen in other works. The NLTE analysis lowers [Mg/Fe] typically by $\sim$0.1 dex, resulting in a noticeably lower trend of [Mg/Fe] versus [Fe/H]. The derived [Ca/Fe] versus [Fe/H] trend has a larger scatter than those for Si and Mg, but is in excellent agreement with local thin and thick disc trends. With our updated analysis, we constructed one of the most detailed studies of the alpha abundance trends of cool M giants in the inner Galactic bulge. We modelled these abundances by adopting a two-infall chemical-evolution model with two distinct gas-infall episodes with timescales of 0.4 Gyr and 2 Gyr, respectively. Based on a very meticulous spectral analysis, we have constructed detailed and precise chemical abundances of Mg, Si, and Ca for cool M giants. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 19 pages, 11 figures

Journal ref: A&A 671, A94 (2023)

The solar photospheric silicon abundance according to CO5BOLD: Investigating line broadening, magnetic fields, and model effects

Authors: S. A. Deshmukh, H. -G. Ludwig, A. Kučinskas, M. Steffen, P. S. Barklem, E. Caffau, V. Dobrovolskas, P. Bonifacio

Abstract: In this work, we present a photospheric solar silicon abundance derived using CO5BOLD model atmospheres and the LINFOR3D spectral synthesis code. Previous works have differed in their choice of a spectral line sample and model atmosphere as well as their treatment of observational material, and the solar silicon abundance has undergone a downward revision in recent years. We additionally show the… ▽ More In this work, we present a photospheric solar silicon abundance derived using CO5BOLD model atmospheres and the LINFOR3D spectral synthesis code. Previous works have differed in their choice of a spectral line sample and model atmosphere as well as their treatment of observational material, and the solar silicon abundance has undergone a downward revision in recent years. We additionally show the effects of the chosen line sample, broadening due to velocity fields, collisional broadening, model spatial resolution, and magnetic fields. CO5BOLD model atmospheres for the Sun were used in conjunction with the LINFOR3D spectral synthesis code to generate model spectra, which were then fit to observations in the Hamburg solar atlas. We present a sample of 11 carefully selected lines (from an initial choice of 39 lines) in the optical and infrared, made possible with newly determined oscillator strengths for the majority of these lines. Our final sample includes seven optical Si I lines, three infrared Si I lines, and one optical Si II line. We derived a photospheric solar silicon abundance of $\log ε_\mathrm{Si} = 7.57 \pm 0.04$, including a $-0.01$ dex correction from Non-Local Thermodynamic Equilibrium (NLTE) effects. Combining this with meteoritic abundances and previously determined photospheric abundances results in a metal mass fraction Z/X = $0.0220 \pm 0.0020$. We found a tendency of obtaining overly broad synthetic lines. We mitigated the impact of this by devising a de-broadening procedure. The over-broadening of synthetic lines does not substantially affect the abundance determined in the end. It is primarily the line selection that affects the final fitted abundance. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: 24 pages

Journal ref: A&A 668, A48 (2022)

Titanium abundances in late-type stars I. 1D non-local thermodynamic equilibrium modelling in benchmark dwarfs and giants

Authors: J. W. E. Mallinson, K. Lind, A. M. Amarsi, P. S. Barklem, J. Grumer, A. K. Belyaev, K. Youakim

Abstract: The titanium abundances of late-type stars are important tracers of Galactic formation history. However, abundances inferred from Ti I and Ti II lines can be in stark disagreement in very metal-poor giants. Departures from local thermodynamic equilibrium (LTE) have a large impact on the minority neutral species and thus influences the ionisation imbalance, but satisfactory non-LTE modelling for bo… ▽ More The titanium abundances of late-type stars are important tracers of Galactic formation history. However, abundances inferred from Ti I and Ti II lines can be in stark disagreement in very metal-poor giants. Departures from local thermodynamic equilibrium (LTE) have a large impact on the minority neutral species and thus influences the ionisation imbalance, but satisfactory non-LTE modelling for both dwarfs and giants has not been achieved in previous literature. The reliability of titanium abundances is reassessed in benchmark dwarfs and giants using a new non-LTE model and one-dimensional (1D) model atmospheres. A comprehensive model atom was compiled with a more extended level structure and newly published data for inelastic collisions between Ti I and neutral hydrogen. In 1D LTE, the Ti I and Ti II lines agree to within $0.06$ dex for the Sun, Arcturus, and the very metal-poor stars HD84937 and HD140283. For the very metal-poor giant HD122563, the Ti I lines give an abundance that is $0.47$ dex lower than that from Ti II. The 1D non-LTE corrections can reach $+0.4$ dex for individual Ti I lines and $+0.1$ dex for individual Ti II lines, and reduce the overall ionisation imbalance to $-0.17$ dex for HD122563. However, it also increases the imbalance for the very metal-poor dwarf and sub-giant to around $0.2$ dex. Using 1D non-LTE reduces the ionisation imbalance in very metal-poor giants but breaks the balance of other very metal-poor stars, consistent with the conclusions in earlier literature. To make further progress, consistent 3D non-LTE models are needed. △ Less

Submitted 7 February, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 9 pages plus appendix, 6 figures; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 668, A103 (2022)

The Gaia-ESO Public Spectroscopic Survey: Motivation, implementation, GIRAFFE data processing, analysis, and final data products

Authors: G. Gilmore, S. Randich, C. C. Worley, A. Hourihane, A. Gonneau, G. G. Sacco, J. R. Lewis, L. Magrini, P. Francois, R. D. Jeffries, S. E. Koposov, A. Bragaglia, E. J. Alfaro, C. Allende Prieto, R. Blomme, A. J. Korn, A. C. Lanzafame, E. Pancino, A. Recio-Blanco, R. Smiljanic, S. Van Eck, T. Zwitter, T. Bensby, E. Flaccomio, M. J. Irwin , et al. (143 additional authors not shown)

Abstract: The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending a… ▽ More The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: 38 pages. A&A in press

Journal ref: A&A 666, A120 (2022)

Non-LTE abundance corrections for late-type stars from 2000Å to 3μm: I. Na, Mg, and Al

Authors: K. Lind, T. Nordlander, A. Wehrhahn, M. Montelius, Y. Osorio, P. S. Barklem, M. Afsar, C. Sneden, C. Kobayashi

Abstract: It is well known that cool star atmospheres depart from local thermodynamic equilibrium (LTE). Accurate abundance determination requires taking those effects into account, but the necessary non-LTE calculations are often lacking. Our goal is to provide detailed estimates of NLTE effects for FGK type stars for all spectral lines from the ultraviolet to the infrared that are potentially useful as ab… ▽ More It is well known that cool star atmospheres depart from local thermodynamic equilibrium (LTE). Accurate abundance determination requires taking those effects into account, but the necessary non-LTE calculations are often lacking. Our goal is to provide detailed estimates of NLTE effects for FGK type stars for all spectral lines from the ultraviolet to the infrared that are potentially useful as abundance diagnostics. The first paper in this series focusses on the light elements Na, Mg and Al. The code PySME is used to compute curves-of-growth for 2158 MARCS model atmospheres in a wide parameter range. Nine abundance points are used to construct individual line curves-of-growth by calculating the equivalent widths of 35 Na lines, 134 Mg lines, and 34 Al lines. The lines are selected from the ultra-violet to the near infrared wavelength range. We demonstrate the power of the new grids with LTE and NLTE abundance analysis by means of equivalent width measurements of five benchmark stars; the Sun, Arcturus, HD84937, HD140283 and HD122563. For Na, the NLTE abundances are lower than in LTE and show markedly reduced line-to-line scatter in the metal-poor stars. For Mg, we confirm previous reports of a significant 0.25 dex LTE ionization imbalance in metal-poor stars that is only slightly improved in NLTE (0.18 dex). LTE abundances based on Mg II lines agree better with models of Galactic chemical evolution. For Al, NLTE calculations strongly reduce a 0.6 dex ionization imbalance seen in LTE for the metal-poor stars. The abundance corrections presented in this work are in good agreement with previous studies for the subset of lines that overlap, except for strongly saturated lines. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 12 pages main text, 7 figures, 14 pages appendix. Accepted for publication in A&A

Journal ref: A&A 665, A33 (2022)

The Gaia-ESO Public Spectroscopic Survey: Implementation, data products, open cluster survey, science, and legacy

Authors: S. Randich, G. Gilmore, L. Magrini, G. G. Sacco, R. J. Jackson, R. D. Jeffries, C. C. Worley, A. Hourihane, A. Gonneau, C. Viscasillas Vàzquez, E. Franciosini, J. R. Lewis, E. J. Alfaro, C. Allende Prieto, T. Bensby R. Blomme, A. Bragaglia, E. Flaccomio, P. François, M. J. Irwin, S. E. Koposov, A. J. Korn, A. C. Lanzafame, E. Pancino, A. Recio-Blanco, R. Smiljanic , et al. (139 additional authors not shown)

Abstract: In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars… ▽ More In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Astronomy and Astrophysics. 30 pages, 30 figures, 4 tables

The molecular chemistry of Type Ibc supernovae, and diagnostic potential with the James Webb Space Telescope

Authors: S. Liljegren, A. Jerkstrand, P. S. Barklem, G. Nyman, R. Brady, S. N. Yurchenko

Abstract: We aim to investigate how the molecular chemistry in stripped-envelope supernovae (SESNe) affect physical conditions and optical spectra, and produce ro-vibrational emission in the mid-infrared (MIR). We also aim to assess the diagnostic potential of observations of such MIR emission with JWST. We coupled a chemical kinetic network including carbon, oxygen, silicon, and sulfur-bearing molecules in… ▽ More We aim to investigate how the molecular chemistry in stripped-envelope supernovae (SESNe) affect physical conditions and optical spectra, and produce ro-vibrational emission in the mid-infrared (MIR). We also aim to assess the diagnostic potential of observations of such MIR emission with JWST. We coupled a chemical kinetic network including carbon, oxygen, silicon, and sulfur-bearing molecules into the nonlocal thermal equilibrium (NLTE) spectral synthesis code SUMO. We let four species - CO, SiO, SiS, and SO - participate in NLTE cooling of the gas to achieve self-consistency between the molecule formation and the temperature. We applied the new framework to model the spectrum of a Type Ic SN in the 100-600d time range. Molecules are predicted to form in SESN ejecta in significant quantities (typical mass $10^{-3}$ $M_\odot$) throughout the 100-600d interval. The impact on the temperature and optical emission depends on the density of the oxygen zones and varies with epoch. For example, the [O I] 6300, 6364 feature can be quenched by molecules from 200 to 450d depending on density. The MIR predictions show strong emission in the fundamental bands of CO, SiO, and SiS, and in the CO and SiO overtones. Type Ibc SN ejecta have a rich chemistry and considering the effect of molecules is important for modeling the temperature and atomic emission in the nebular phase. Observations of SESNe with JWST hold promise to provide the first detections of SiS and SO, and to give information on zone masses and densities of the ejecta. Combined optical, near-infrared, and MIR observations can break degeneracies and achieve a more complete picture of the nucleosynthesis, chemistry, and origin of Type Ibc SNe. △ Less

Submitted 25 June, 2023; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: 34 pages, 26 figures

Journal ref: A&A 674, A184 (2023)

Correspondence between the surface integral and linear combination of atomic orbitals methods for ionic-covalent interactions in mutual neutralisation processes involving H$^-$/D$^-$

Authors: Paul S. Barklem

Abstract: The surface integral method for estimating ionic-covalent interactions in diatomic systems been successful in producing cross sections for mutual neutralisation (MN) in reasonable agreement with experimental results for branching fractions between final states in systems such as O$^+$/O$^-$ and N$^+$/O$^-$. However, for simpler cases of MN involving H$^-$ or D$^-$, such as Li$^+$/D$^-$ and Na$^+$/… ▽ More The surface integral method for estimating ionic-covalent interactions in diatomic systems been successful in producing cross sections for mutual neutralisation (MN) in reasonable agreement with experimental results for branching fractions between final states in systems such as O$^+$/O$^-$ and N$^+$/O$^-$. However, for simpler cases of MN involving H$^-$ or D$^-$, such as Li$^+$/D$^-$ and Na$^+$/D$^-$, it has not produced results that are in agreement with experiments and other theoretical calculations; in particular, for Li$^+$/D$^-$ calculations predict the wrong ordering of importance of final channels, including the incorrect most populated channel. The reason for this anomaly is investigated and a leading constant to the asymptotic H$^-$ wavefunction is found that is different by roughly a factor $1/\sqrt{2}$ to that which has been used in previous calculations with the surface integral method involving H$^-$ or D$^-$. With this correction, far better agreement with both experimental results and with calculations with full quantum and LCAO methods is obtained. Further, it is shown that the surface integral method and LCAO methods have the same asymptotic behaviour, in contrast to previous claims. This result suggests the surface integral method, which is comparatively easy to calculate, has greater potential for estimating MN processes than earlier comparisons had suggested. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: Accepted by Physical Review A

Diagnostic capabilities of spectropolarimetric observations for understanding solar phenomena I. Zeeman-sensitive photospheric lines

Authors: C. Quintero Noda, P. S. Barklem, R. Gafeira, B. Ruiz Cobo, M. Collados, M. Carlsson, V. Martínez Pillet, D. Orozco Suárez, H. Uitenbroek, Y. Katsukawa

Abstract: Future ground-based telescopes will expand our capabilities for simultaneous multi-line polarimetric observations in a wide range of wavelengths, from the near-ultraviolet to the near-infrared. This creates a strong demand to compare candidate spectral lines to establish a guideline of the lines that are most appropriate for each observation target. We focused in this first work on Zeeman-sensitiv… ▽ More Future ground-based telescopes will expand our capabilities for simultaneous multi-line polarimetric observations in a wide range of wavelengths, from the near-ultraviolet to the near-infrared. This creates a strong demand to compare candidate spectral lines to establish a guideline of the lines that are most appropriate for each observation target. We focused in this first work on Zeeman-sensitive photospheric lines in the visible and infrared. We first examined their polarisation signals and response functions using a 1D semi-empirical atmosphere. Then we studied the spatial distribution of the line core intensity and linear and circular polarisation signals using a realistic 3D numerical simulation. We ran inversions of synthetic profiles, and we compared the heights at which we obtain a high correlation between the input and the inferred atmosphere. We also used this opportunity to revisit the atomic information we have on these lines and computed the broadening cross-sections due to collisions with neutral hydrogen atoms for all the studied spectral lines. The results reveal that four spectral lines stand out from the rest for quiet-Sun and network conditions: Fe I 5250.2, 6302, 8468, and 15648 A. The first three form higher in the atmosphere, and the last line is mainly sensitive to the atmospheric parameters at the bottom of the photosphere. However, as they reach different heights, we strongly recommend using at least one of the first three candidates together with the Fe I 15648 A line to optimise our capabilities for inferring the thermal and magnetic properties of the lower atmosphere. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 15 pages, 12 figures, regular publication

Journal ref: A&A 652, A161 (2021)

Mutual neutralisation in Li$^+$+H$^-$/D$^-$ and Na$^+$+H$^-$/D$^-$ collisions: Implications of experimental results for non-LTE modelling of stellar spectra

Authors: Paul S. Barklem, Anish M. Amarsi, Jon Grumer, Gustav Eklund, Stefan Rosén, MingChao Ji, Henrik Cederquist, Henning Zettergren, Henning T. Schmidt

Abstract: Advances in merged-beams instruments have allowed experimental studies of the mutual neutralisation (MN) processes in collisions of both Li$^+$ and Na$^+$ ions with D$^-$ at energies below 1 eV. These experimental results place constraints on theoretical predictions of MN processes of Li$^+$ and Na$^+$ with H$^-$, important for non-LTE modelling of Li and Na spectra in late-type stars. We compare… ▽ More Advances in merged-beams instruments have allowed experimental studies of the mutual neutralisation (MN) processes in collisions of both Li$^+$ and Na$^+$ ions with D$^-$ at energies below 1 eV. These experimental results place constraints on theoretical predictions of MN processes of Li$^+$ and Na$^+$ with H$^-$, important for non-LTE modelling of Li and Na spectra in late-type stars. We compare experimental results with calculations for methods typically used to calculate MN processes, namely the full quantum (FQ) approach, and asymptotic model approaches based on the linear combination of atomic orbitals (LCAO) and semi-empirical (SE) methods for deriving couplings. It is found that FQ calculations compare best overall with the experiments, followed by the LCAO, and the SE approaches. The experimental results together with the theoretical calculations, allow us to investigate the effects on modelled spectra and derived abundances and their uncertainties arising from uncertainties in the MN rates. Numerical experiments in a large grid of 1D model atmospheres, and a smaller set of 3D models, indicate that neglect of MN can lead to abundance errors of up to 0.1 dex (26\%) for Li at low metallicity, and 0.2 dex (58\%) for Na at high metallicity, while the uncertainties in the relevant MN rates as constrained by experiments correspond to uncertainties in abundances of much less than 0.01~dex (2\%). This agreement for simple atoms gives confidence in the FQ, LCAO and SE model approaches to be able to predict MN with the accuracy required for non-LTE modelling in stellar atmospheres. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: Accepted by ApJ

Atomic data for the Gaia-ESO Survey

Authors: Ulrike Heiter, Karin Lind, Maria Bergemann, Martin Asplund, Šarunas Mikolaitis, Paul S. Barklem, Thomas Masseron, Patrick de Laverny, Laura Magrini, Bengt Edvardsson, Henrik Jönsson, Juliet C. Pickering, Nils Ryde, Amelia Bayo Arán, Thomas Bensby, Andrew R. Casey, Sofia Feltzing, Paula Jofré, Andreas J. Korn, Elena Pancino, Francesco Damiani, Alessandro Lanzafame, Carmela Lardo, Lorenzo Monaco, Lorenzo Morbidelli , et al. (5 additional authors not shown)

Abstract: We describe the atomic and molecular data that were used for the abundance analyses of FGK-type stars carried out within the Gaia-ESO Survey. We present an unprecedented effort to create a homogeneous line list, which was used by several abundance analysis groups to calculate synthetic spectra and equivalent widths. The atomic data are accompanied by quality indicators and detailed references to t… ▽ More We describe the atomic and molecular data that were used for the abundance analyses of FGK-type stars carried out within the Gaia-ESO Survey. We present an unprecedented effort to create a homogeneous line list, which was used by several abundance analysis groups to calculate synthetic spectra and equivalent widths. The atomic data are accompanied by quality indicators and detailed references to the sources. The atomic and molecular data are made publicly available in electronic form. In general experimental transition probabilities were preferred but theoretical values were also used. Astrophysical gf-values were avoided due to the model-dependence of such a procedure. For elements whose lines are significantly affected by hyperfine structure or isotopic splitting a concerted effort has been made to collate the necessary data for the individual line components. We also performed a detailed investigation of available data for line broadening due to collisions with neutral hydrogen atoms. Synthetic spectra calculated for the Sun and Arcturus were used to assess the blending properties of the lines. Among a subset of over 1300 lines of 35 elements in the wavelength ranges from 475 nm to 685 nm and from 850 nm to 895 nm we identified about 200 lines of 24 species which have accurate gf-values and are free of blends in the spectra of the Sun and Arcturus. For the broadening due to collisions with neutral hydrogen we recommend data based on Anstee-Barklem-O'Mara theory, where available, and to avoid lines of neutral species otherwise. Theoretical broadening data by R.L. Kurucz should be used for Sc II, Ti II, and Y II lines. For ionised rare-earth species the Unsöld approximation with an enhancement factor of 1.5 for the line width can be used. Desirable improvements in atomic data were identified for a number of species, including Al I, S I, Cr II, Na I, Si I, Ca II, and Ni I. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: Accepted for publication in section 12. Atomic, molecular, and nuclear data of Astronomy and Astrophysics; main part 25 pages, 8 tables, 5 figures; appendices 53 pages, 23 tables, 29 figures

Journal ref: A&A 645, A106 (2021)

The GALAH Survey: Non-LTE departure coefficients for large spectroscopic surveys

Authors: A. M. Amarsi, K. Lind, Y. Osorio, T. Nordlander, M. Bergemann, H. Reggiani, E. X. Wang, S. Buder, M. Asplund, P. S. Barklem, A. Wehrhahn, Á. Skúladóttir, C. Kobayashi, A. I. Karakas, X. D. Gao, J. Bland-Hawthorn, G. M. De Silva, J. Kos, G. F. Lewis, S. L. Martell, S. Sharma, J. D. Simpson, D. B. Zucker, K. Čotar, J. Horner , et al. (1 additional authors not shown)

Abstract: Massive sets of stellar spectroscopic observations are rapidly becoming available and these can be used to determine the chemical composition and evolution of the Galaxy with unprecedented precision. One of the major challenges in this endeavour involves constructing realistic models of stellar spectra with which to reliably determine stellar abundances. At present, large stellar surveys commonly… ▽ More Massive sets of stellar spectroscopic observations are rapidly becoming available and these can be used to determine the chemical composition and evolution of the Galaxy with unprecedented precision. One of the major challenges in this endeavour involves constructing realistic models of stellar spectra with which to reliably determine stellar abundances. At present, large stellar surveys commonly use simplified models that assume that the stellar atmospheres are approximately in local thermodynamic equilibrium (LTE). To test and ultimately relax this assumption, we have performed non-LTE calculations for $13$ different elements (H, Li, C, N, O, Na, Mg, Al, Si, K, Ca, Mn, and Ba), using recent model atoms that have physically-motivated descriptions for the inelastic collisions with neutral hydrogen, across a grid of $3756$ 1D MARCS model atmospheres that spans $3000\leq T_{\mathrm{eff}}/\mathrm{K}\leq8000$, $-0.5\leq\log{g/\mathrm{cm\,s^{-2}}}\leq5.5$, and $-5\leq\mathrm{[Fe/H]}\leq1$. We present the grids of departure coefficients that have been implemented into the GALAH DR3 analysis pipeline in order to complement the extant non-LTE grid for iron. We also present a detailed line-by-line re-analysis of $50126$ stars from GALAH DR3. We found that relaxing LTE can change the abundances by between $-0.7\,\mathrm{dex}$ and $+0.2\,\mathrm{dex}$ for different lines and stars. Taking departures from LTE into account can reduce the dispersion in the $\mathrm{[A/Fe]}$ versus $\mathrm{[Fe/H]}$ plane by up to $0.1\,\mathrm{dex}$, and it can remove spurious differences between the dwarfs and giants by up to $0.2\,\mathrm{dex}$. The resulting abundance slopes can thus be qualitatively different in non-LTE, possibly with important implications for the chemical evolution of our Galaxy. △ Less

Submitted 21 August, 2020; originally announced August 2020.

Comments: 19 pages, 25 figures, 2 tables, arXiv abstract abridged; accepted for publication in A&A

Journal ref: A&A 642, A62 (2020)

The 3D non-LTE solar nitrogen abundance from atomic lines

Authors: A. M. Amarsi, N. Grevesse, J. Grumer, M. Asplund, P. S. Barklem, R. Collet

Abstract: Nitrogen is an important element in various fields of stellar and Galactic astronomy, and the solar nitrogen abundance is crucial as a yardstick for comparing different objects in the cosmos. In order to obtain a precise and accurate value for this abundance, we carried out N i line formation calculations in a 3D radiative-hydrodynamic STAGGER model solar atmosphere, in full 3D non-local thermodyn… ▽ More Nitrogen is an important element in various fields of stellar and Galactic astronomy, and the solar nitrogen abundance is crucial as a yardstick for comparing different objects in the cosmos. In order to obtain a precise and accurate value for this abundance, we carried out N i line formation calculations in a 3D radiative-hydrodynamic STAGGER model solar atmosphere, in full 3D non-local thermodynamic equilibrium (non-LTE), using a model atom that includes physically-motivated descriptions for the inelastic collisions of N i with free electrons and with neutral hydrogen. We selected five N i lines of high excitation energy to study in detail, based on their strengths and on their being relatively free of blends. We found that these lines are slightly strengthened from non-LTE photon losses and from 3D granulation effects, resulting in negative abundance corrections of around $-0.01$ dex and $-0.04$ dex respectively. Our advocated solar nitrogen abundance is $\logε_{\mathrm{N}} = 7.77$, with the systematic $1σ$ uncertainty estimated to be $0.05$ dex. This result is consistent with earlier studies after correcting for differences in line selections and equivalent widths. △ Less

Submitted 31 March, 2020; v1 submitted 27 March, 2020; originally announced March 2020.

Comments: 11 pages, 7 figures, 3 tables; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 636, A120 (2020)

Excitation and charge transfer in low-energy hydrogen atom collisions with neutral manganese and titanium

Authors: Jon Grumer, Paul S. Barklem

Abstract: Data for inelastic processes due to hydrogen atom collisions with manganese and titanium are needed for accurate modeling of the corresponding spectra in late-type stars. In this work excitation and charge transfer in low-energy Mn+H and Ti+H collisions have been studied theoretically using a method based on an asymptotic two-electron linear combination of an atomic orbitals model of ionic-covalen… ▽ More Data for inelastic processes due to hydrogen atom collisions with manganese and titanium are needed for accurate modeling of the corresponding spectra in late-type stars. In this work excitation and charge transfer in low-energy Mn+H and Ti+H collisions have been studied theoretically using a method based on an asymptotic two-electron linear combination of an atomic orbitals model of ionic-covalent interactions in the neutral atom-hydrogen-atom system, together with the multichannel Landau-Zener model to treat the dynamics. Extensive calculations of charge transfer (mutual neutralization, ion-pair production), excitation and de-excitation processes in the two collisional systems are carried out for all transitions between covalent states dissociating to energies below the first ionic limit and the dominating ionic states. Rate coefficients are determined for temperatures in the range 1000 - 20 000 K in steps of 1000 K. Like for earlier studies of other atomic species, charge transfer processes are found to lead to much larger rate coefficients than excitation processes. △ Less

Submitted 8 April, 2020; v1 submitted 22 March, 2020; originally announced March 2020.

Comments: 9 pages, 3 figures, 5 tables; accepted for publication in A&A

Journal ref: A&A 637, A28 (2020)

Non-LTE analysis of K I in late-type stars

Authors: Henrique Reggiani, Anish M. Amarsi, Karin Lind, Paul S. Barklem, Oleg Zatsarinny, Klaus Bartschat, Dmitry V. Fursa, Igor Bray, Lorenzo Spina, Jorge Meléndez

Abstract: Older GCE models 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 models and observations is in part caused by the LTE assumption. We study the statistical equilibrium of KI, focusing on the non-LTE effects on the $7698 \ Å$ line. We aim to det… ▽ More Older GCE models 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 models and observations is in part caused by the LTE assumption. We study the statistical equilibrium of KI, focusing on the non-LTE effects on the $7698 \ Å$ line. We aim to determine how non-LTE abundances of K can improve the analysis of its chemical evolution, and help to constrain the yields of models. We construct a model atom that employs the most up-to-date data. In particular, we calculate and present inelastic e+K collisional excitation cross-sections from the convergent close-coupling and the $B$-Spline $R$-matrix methods, and H+K collisions from the two-electron model. We constructed a fine grid of non-LTE abundance corrections that span $4000<\teff / \rm{K}<8000$, $0.50<\lgg<5.00$, $-5.00<\feh<+0.50$, and applied the corrections to abundances from the literature. In concordance with previous studies, we find severe non-LTE effects in the $7698 \ Å$ line, which is stronger in non-LTE with abundance corrections that can reach $\sim-0.7\,\dex$. We explore the effects of atmospheric inhomogeneity by computing a full 3D non-LTE stellar spectrum of KI for a test star. We find that 3D 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. Our non-LTE abundance corrections reduce the scatter and change the cosmic trends of literature K abundances. In the regime [Fe/H]$\lesssim-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 modelling errors. △ Less

Submitted 19 June, 2019; originally announced June 2019.

Comments: 18 pages, 15 Figures

Journal ref: A&A 627, A177 (2019)

Excitation and charge transfer in low-energy hydrogen atom collisions with neutral carbon and nitrogen

Authors: A. M. Amarsi, P. S. Barklem

Abstract: Low-energy inelastic collisions with neutral hydrogen atoms are important processes in stellar atmospheres, and a persistent source of uncertainty in non-LTE modelling of stellar spectra. We have calculated and studied excitation and charge transfer of C i and of N i due to such collisions. We used a previously presented method that is based on an asymptotic two-electron linear combination of atom… ▽ More Low-energy inelastic collisions with neutral hydrogen atoms are important processes in stellar atmospheres, and a persistent source of uncertainty in non-LTE modelling of stellar spectra. We have calculated and studied excitation and charge transfer of C i and of N i due to such collisions. We used a previously presented method that is based on an asymptotic two-electron linear combination of atomic orbitals (LCAO) model of ionic-covalent interactions for the adiabatic potential energies, combined with the multichannel Landau-Zener model for the collision dynamics. We find that charge transfer processes typically lead to much larger rate coefficients than excitation processes do, consistent with studies of other atomic species. Two-electron processes were considered and lead to non-zero rate coefficients that can potentially impact statistical equilibrium calculations. However, they were included in the model in an approximate way, via an estimate for the two-electron coupling that was presented earlier in the literature: the validity of these data should be checked in a future work. △ Less

Submitted 23 May, 2019; v1 submitted 30 March, 2019; originally announced April 2019.

Comments: 11 pages, 5 figures, 4 tables; accepted for publication in A&A

Journal ref: A&A 625, A78 (2019)

3D non-LTE line formation of neutral carbon in the Sun

Authors: A. M. Amarsi, P. S. Barklem, R. Collet, N. Grevesse, M. Asplund

Abstract: Carbon abundances in late-type stars are important in a variety of astrophysical contexts. However C i lines, one of the main abundance diagnostics, are sensitive to departures from local thermodynamic equilibrium (LTE). We present a model atom for non-LTE analyses of C i lines, that uses a new, physically-motivated recipe for the rates of neutral hydrogen impact excitation. We analyse C i lines i… ▽ More Carbon abundances in late-type stars are important in a variety of astrophysical contexts. However C i lines, one of the main abundance diagnostics, are sensitive to departures from local thermodynamic equilibrium (LTE). We present a model atom for non-LTE analyses of C i lines, that uses a new, physically-motivated recipe for the rates of neutral hydrogen impact excitation. We analyse C i lines in the solar spectrum, employing a three-dimensional (3D) hydrodynamic model solar atmosphere and 3D non-LTE radiative transfer. We find negative non-LTE abundance corrections for C i lines in the solar photosphere, in accordance with previous studies, reaching up to around 0.1 dex in the disk-integrated flux. We also present the first fully consistent 3D non-LTE solar carbon abundance determination: we infer log $ε_{\text{C}}$ = $8.44\pm0.02$, in good agreement with the current standard value. Our models reproduce the observed solar centre-to-limb variations of various C i lines, without any adjustments to the rates of neutral hydrogen impact excitation, suggesting that the proposed recipe may be a solution to the long-standing problem of how to reliably model inelastic collisions with neutral hydrogen in late-type stellar atmospheres. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 14 pages, 8 figures, 2 tables; accepted for publication in A&A

Journal ref: A&A 624, A111 (2019)

Ca line formation in late-type stellar atmospheres: I. The model atom

Authors: Y. Osorio, K. Lind, P. S. Barklem, C. Allende Prieto, Oleg Zatsarinny

Abstract: Context. Departures from local thermodynamic equilibrium (LTE) distort the calcium abundance derived from stellar spectra in various ways, depending on the lines used and the stellar atmospheric parameters. The collection of atomic data adopted in non-LTE (NLTE) calculations must be sufficiently complete and accurate. Aims. We derive NLTE abundances from high-quality observations and reliable stel… ▽ More Context. Departures from local thermodynamic equilibrium (LTE) distort the calcium abundance derived from stellar spectra in various ways, depending on the lines used and the stellar atmospheric parameters. The collection of atomic data adopted in non-LTE (NLTE) calculations must be sufficiently complete and accurate. Aims. We derive NLTE abundances from high-quality observations and reliable stellar parameters using a model atom built afresh for this work, and check the consistency of our results over a wide wavelength range with transitions of atomic and singly ionised calcium. Methods. We built and tested Ca i and Ca ii model atoms with state-of-the-art radiative and collisional data, and tested their performance deriving the Ca abundance in three benchmark stars: Procyon, the Sun, and Arcturus. We have excellent-quality observations and accurate stellar parameters for these stars. Two methods to derive the LTE / NLTE abundances were used and compared. The LTE / NLTE centre-to-limb variation (CLV) of Ca lines in the Sun was also investigated. Results. The two methods used give similar results in all three stars. Several discrepancies found in LTE do not appear in our NLTE results; in particular the agreement between abundances in the visual and infra-red (IR) and the Ca i and Ca ii ionisation balance is improved overall, although substantial line-to-line scatter remains. The CLV of the calcium lines around 6165 Angs can be partially reproduced. We suspect differences between our modelling and CLV results are due to inhomogeneities in the atmosphere that require 3D modelling. △ Less

Submitted 31 January, 2019; originally announced January 2019.

Comments: 17 pages, 11 Figures

Journal ref: A&A 623, A103 (2019)

Carbon and oxygen in metal-poor halo stars

Authors: A. M. Amarsi, P. E. Nissen, M. Asplund, K. Lind, P. S. Barklem

Abstract: Carbon and oxygen are key tracers of the Galactic chemical evolution; in particular, a reported upturn in [C/O] towards decreasing [O/H] in metal-poor halo stars could be a signature of nucleosynthesis by massive Population III stars. We reanalyse carbon, oxygen, and iron abundances in thirty-nine metal-poor turn-off stars. For the first time, we take into account three-dimensional (3D) hydrodynam… ▽ More Carbon and oxygen are key tracers of the Galactic chemical evolution; in particular, a reported upturn in [C/O] towards decreasing [O/H] in metal-poor halo stars could be a signature of nucleosynthesis by massive Population III stars. We reanalyse carbon, oxygen, and iron abundances in thirty-nine metal-poor turn-off stars. For the first time, we take into account three-dimensional (3D) hydrodynamic effects together with departures from local thermodynamic equilibrium (LTE) when determining both the stellar parameters and the elemental abundances, by deriving effective temperatures from 3D non-LTE H$β$ profiles, surface gravities from Gaia parallaxes, iron abundances from 3D LTE Feii equivalent widths, and carbon and oxygen abundances from 3D non-LTE Ci and Oi equivalent widths. We find that [C/Fe] stays flat with [Fe/H], whereas [O/Fe] increases linearly up to $0.75$ dex with decreasing [Fe/H] down to $-3.0$ dex. As such [C/O] monotonically decreases towards decreasing [O/H], in contrast to previous findings, mainly by virtue of less severe non-LTE effects for Oi at low [Fe/H] with our improved calculations. △ Less

Submitted 1 February, 2019; v1 submitted 11 January, 2019; originally announced January 2019.

Comments: 5 pages, 2 figures; published in A&A Letters

Journal ref: A&A 622, L4 (2019)

A collection of model stellar spectra for spectral types B to early-M

Authors: Carlos Allende Prieto, Lars Koesterke, Ivan Hubeny, Manuel A. Bautista, Paul S. Barklem, Sultana N. Nahar

Abstract: Models of stellar spectra are necessary for interpreting light from individual stars, planets, integrated stellar populations, nebulae, and the interstellar medium. We provide a comprehensive and homogeneous collection of synthetic spectra for a wide range of atmospheric parameters and chemical compositions. We compile atomic and molecular data from the literature. We adopt the largest and most re… ▽ More Models of stellar spectra are necessary for interpreting light from individual stars, planets, integrated stellar populations, nebulae, and the interstellar medium. We provide a comprehensive and homogeneous collection of synthetic spectra for a wide range of atmospheric parameters and chemical compositions. We compile atomic and molecular data from the literature. We adopt the largest and most recent set of ATLAS9 model atmospheres, and use the radiative code ASS$ε$T.The resulting collection of spectra is made publicly available at medium and high-resolution ($R\equivλ/δλ$ = 10,000, 100,000 and 300,000 spectral grids, which include variations in effective temperature between 3500 K and 30,000 K, surface gravity ($0\le \log g \le 5$), and metallicity ($-5 \le$[Fe/H]$\le +0. 5$), spanning the wavelength interval 120-6500 nm. A second set of denser grids with additional dimensions, [$α$/Fe] and micro-turbulence, are also provided (covering 200-2500 nm). We compare models with observations for a few representative cases. △ Less

Submitted 16 July, 2018; originally announced July 2018.

Comments: 7 pages, 4 figures. To appear in A&A

Journal ref: A&A 618, A25 (2018)

Chemical Characterization of the Inner Galactic bulge:North-South Symmetry

Authors: G. Nandakumar, N. Ryde, M. Schultheis, B. Thorsbro, H. Jönsson, P. S. Barklem, R. M. Rich, F. Fragkoudi

Abstract: While the number of stars in the Galactic bulge with detailed chemical abundance measurements is increasingly rapidly, the inner Galactic bulge ( |b| < 2$^\circ$) remains poorly studied, due to heavy interstellar absorption and photometric crowding. We have carried out a high-resolution IR spectroscopic study of 72 M giants in the inner bulge using the CRIRES (ESO/VLT) facility. Our spectra cover… ▽ More While the number of stars in the Galactic bulge with detailed chemical abundance measurements is increasingly rapidly, the inner Galactic bulge ( |b| < 2$^\circ$) remains poorly studied, due to heavy interstellar absorption and photometric crowding. We have carried out a high-resolution IR spectroscopic study of 72 M giants in the inner bulge using the CRIRES (ESO/VLT) facility. Our spectra cover the wavelength range of 2.0818 - 2.1444 $μ$m with the resolution of R$\sim$50,000 and have signal-to-noise ratio of 50 - 100. Our stars are located along the bulge minor axis at l = 0$^\circ$, b = $\pm$0$^\circ$, $\pm$1$^\circ$, $\pm$2$^\circ$ and +3$^\circ$, selected with the aim of investigating any North-South asymmetries in the metallicity distribution function and composition, and comparing them to the outer bulge fields. Our sample was analysed in a homogeneous way using the most current K-band line list. We clearly detect a bimodal MDF with a metal-rich peak at $\rm \sim +0.3\,dex$ and a metal-poor peak at $\rm \sim -0.5\,dex$. Only a single star is found to exceed $\rm [Fe/H]=+0.5\,dex$. They show a symmetric behaviour along the $\pm$1$^\circ$, $\pm$2$^\circ$ fields. The Galactic Center field reveals in contrast a mainly metal-rich population with a mean metallicity of $\rm +0.3\,dex$. We derived $\rm [Mg/Fe]$ and $\rm [Si/Fe]$ abundances which are consistent with trends from the outer bulge, with a gradually decreasing trend with increasing metallicity. We confirm for the supersolar metallicity stars the decreasing trend in \mgfe\ and \sife\ as expected from chemical evolution models. With the caveat of a relatively small sample, we do not find significant differences in the chemical abundances between the Northern and the Southern fields, hence the evidence is consistent with symmetry in chemistry between North and South. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Comments: 17 pages, 15 figures, Accepted for publication in MNRAS

Gaia Data Release 2: Observational Hertzsprung-Russell diagrams

Authors: Gaia Collaboration, C. Babusiaux, F. van Leeuwen, M. A. Barstow, C. Jordi, A. Vallenari, D. Bossini, A. Bressan, T. Cantat-Gaudin, M. van Leeuwen, A. G. A. Brown, T. Prusti, J. H. J. de Bruijne, C. A. L. Bailer-Jones, M. Biermann, D. W. Evans, L. Eyer, F. Jansen, S. A. Klioner, U. Lammers, L. Lindegren, X. Luri, F. Mignard, C. Panem, D. Pourbaix , et al. (428 additional authors not shown)

Abstract: We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the select… ▽ More We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies. △ Less

Submitted 13 August, 2018; v1 submitted 25 April, 2018; originally announced April 2018.

Comments: Published in the A&A Gaia Data Release 2 special issue. Tables 2 and A.4 corrected. Tables available at http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/616/A10

Journal ref: A&A 616, A10 (2018)

Effective temperature determinations of late-type stars based on 3D non-LTE Balmer line formation

Authors: A. M. Amarsi, T. Nordlander, P. S. Barklem, M. Asplund, R. Collet, K. Lind

Abstract: Hydrogen Balmer lines are commonly used as spectroscopic effective temperature diagnostics of late-type stars. However, the absolute accuracy of classical methods that are based on one-dimensional (1D) hydrostatic model atmospheres and local thermodynamic equilibrium (LTE) is still unclear. To investigate this, we carry out 3D non-LTE calculations for the Balmer lines, performed, for the first tim… ▽ More Hydrogen Balmer lines are commonly used as spectroscopic effective temperature diagnostics of late-type stars. However, the absolute accuracy of classical methods that are based on one-dimensional (1D) hydrostatic model atmospheres and local thermodynamic equilibrium (LTE) is still unclear. To investigate this, we carry out 3D non-LTE calculations for the Balmer lines, performed, for the first time, over an extensive grid of 3D hydrodynamic STAGGER model atmospheres. For H$α$, H$β$, and H$γ$, we find significant 1D non-LTE versus 3D non-LTE differences (3D effects): the outer wings tend to be stronger in 3D models, particularly for H$γ$, while the inner wings can be weaker in 3D models, particularly for H$α$. For H$α$, we also find significant 3D LTE versus 3D non-LTE differences (non-LTE effects): in warmer stars ($T_{\text{eff}}\approx6500$K) the inner wings tend to be weaker in non-LTE models, while at lower effective temperatures ($T_{\text{eff}}\approx4500$K) the inner wings can be stronger in non-LTE models; the non-LTE effects are more severe at lower metallicities. We test our 3D non-LTE models against observations of well-studied benchmark stars. For the Sun, we infer concordant effective temperatures from H$α$, H$β$, and H$γ$; however the value is too low by around 50K which could signal residual modelling shortcomings. For other benchmark stars, our 3D non-LTE models generally reproduce the effective temperatures to within $1σ$ uncertainties. For H$α$, the absolute 3D effects and non-LTE effects can separately reach around 100K, in terms of inferred effective temperatures. For metal-poor turn-off stars, 1D LTE models of H$α$ can underestimate effective temperatures by around 150K. Our 3D non-LTE model spectra are publicly available, and can be used for more reliable spectroscopic effective temperature determinations. △ Less

Submitted 2 August, 2018; v1 submitted 6 April, 2018; originally announced April 2018.

Comments: 19 pages, 10 figures, abstract abridged; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 615, A139 (2018)

Inelastic O+H collisions and the OI 777nm solar centre-to-limb variation

Authors: A. M. Amarsi, P. S. Barklem, M. Asplund, R. Collet, O. Zatsarinny

Abstract: The OI 777 nm triplet is a key diagnostic of oxygen abundances in the atmospheres of FGK-type stars; however it is sensitive to departures from local thermodynamic equilibrium (LTE). The accuracy of non-LTE line formation calculations has hitherto been limited by errors in the inelastic O+H collisional rate coefficients: several recent studies have used the so-called Drawin recipe, albeit with a c… ▽ More The OI 777 nm triplet is a key diagnostic of oxygen abundances in the atmospheres of FGK-type stars; however it is sensitive to departures from local thermodynamic equilibrium (LTE). The accuracy of non-LTE line formation calculations has hitherto been limited by errors in the inelastic O+H collisional rate coefficients: several recent studies have used the so-called Drawin recipe, albeit with a correction factor $\mathrm{S_{H}}$ that is calibrated to the solar centre-to-limb variation of the triplet. We present a new model oxygen atom that incorporates inelastic O+H collisional rate coefficients using an asymptotic two-electron model based on linear combinations of atomic orbitals, combined with a free electron model, based on the impulse approximation. Using a 3D hydrodynamic stagger model solar atmosphere and 3D non-LTE line formation calculations, we demonstrate that this physically-motivated approach is able to reproduce the solar centre-to-limb variation of the triplet to 0.02 dex, without any calibration of the inelastic collisional rate coefficients or other free parameters. We infer $\logε_{\mathrm{O}}=8.69\pm0.03$ from the triplet alone, strengthening the case for a low solar oxygen abundance. △ Less

Submitted 28 October, 2018; v1 submitted 28 March, 2018; originally announced March 2018.

Comments: 13 pages, 8 figures; published in Astronomy & Astrophysics

Journal ref: A&A 616, A89 (2018)

A New Test of Copper and Zinc Abundances in Late-Type Stars Using Ultraviolet Cu II and Zn II Lines

Authors: Ian U. Roederer, Paul S. Barklem

Abstract: We present new abundances derived from Cu I, Cu II, Zn I, and Zn II lines in six warm (5766 < Teff < 6427 K), metal-poor (-2.50 < [Fe/H] < -0.95) dwarf and subgiant (3.64 < log g < 4.44) stars. These abundances are derived from archival high-resolution ultraviolet spectra from the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope and ground-based optical spectra from several… ▽ More We present new abundances derived from Cu I, Cu II, Zn I, and Zn II lines in six warm (5766 < Teff < 6427 K), metal-poor (-2.50 < [Fe/H] < -0.95) dwarf and subgiant (3.64 < log g < 4.44) stars. These abundances are derived from archival high-resolution ultraviolet spectra from the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope and ground-based optical spectra from several observatories. Ionized Cu and Zn are the majority species, and abundances derived from Cu II and Zn II lines should be largely insensitive to departures from local thermodynamic equilibrium (LTE). We find good agreement between the [Zn/H] ratios derived separately from Zn I and Zn II lines, suggesting that departures from LTE are, at most, minimal (< 0.1 dex). We find that the [Cu/H] ratios derived from Cu II lines are 0.36 +/- 0.06 dex larger than those derived from Cu I lines in the most metal-poor stars ([Fe/H] < -1.8), suggesting that LTE underestimates the Cu abundance derived from Cu I lines. The deviations decrease in more metal-rich stars. Our results validate previous theoretical non-LTE calculations for both Cu and Zn, supporting earlier conclusions that the enhancement of [Zn/Fe] in metal-poor stars is legitimate, and the deficiency of [Cu/Fe] in metal-poor stars may not be as large as previously thought. △ Less

Submitted 26 March, 2018; originally announced March 2018.

Comments: Accepted for publication in the Astrophysical Journal (15 pages, 5 figures, 8 tables)

Excitation and charge transfer in low-energy hydrogen atom collisions with neutral iron

Authors: P. S. Barklem

Abstract: Data for inelastic processes due to hydrogen atom collisions with iron are needed for accurate modelling of the iron spectrum in late-type stars. Excitation and charge transfer in low-energy Fe+H collisions is studied theoretically using a previously presented method based on an asymptotic two-electron linear combination of atomic orbitals (LCAO) model of ionic-covalent interactions in the neutral… ▽ More Data for inelastic processes due to hydrogen atom collisions with iron are needed for accurate modelling of the iron spectrum in late-type stars. Excitation and charge transfer in low-energy Fe+H collisions is studied theoretically using a previously presented method based on an asymptotic two-electron linear combination of atomic orbitals (LCAO) model of ionic-covalent interactions in the neutral atom-hydrogen-atom system, together with the multi-channel Landau-Zener model. An extensive calculation including 166 covalent states and 25 ionic states is presented and rate coefficients are calculated for temperatures in the range 1000 - 20000 K. The largest rates are found for charge transfer processes to and from two clusters of states around 6.3 and 6.6 eV excitation, corresponding in both cases to active 4d and 5p electrons undergoing transfer. Excitation and de-excitation processes among these two sets of states are also significant. △ Less

Submitted 22 January, 2018; originally announced January 2018.

Comments: Accepted by A&A

Journal ref: A&A 612, A90 (2018)

Excitation and charge transfer in low-energy hydrogen atom collisions with neutral oxygen

Authors: P. S. Barklem

Abstract: Excitation and charge transfer in low-energy O+H collisions is studied; it is a problem of importance for modelling stellar spectra and obtaining accurate oxygen abundances in late-type stars including the Sun. The collisions have been studied theoretically using a previously presented method based on an asymptotic two-electron linear combination of atomic orbitals (LCAO) model of ionic-covalent i… ▽ More Excitation and charge transfer in low-energy O+H collisions is studied; it is a problem of importance for modelling stellar spectra and obtaining accurate oxygen abundances in late-type stars including the Sun. The collisions have been studied theoretically using a previously presented method based on an asymptotic two-electron linear combination of atomic orbitals (LCAO) model of ionic-covalent interactions in the neutral atom-hydrogen-atom system, together with the multichannel Landau-Zener model. The method has been extended to include configurations involving excited states of hydrogen using an estimate for the two-electron transition coupling, but this extension was found to not lead to any remarkably high rates. Rate coefficients are calculated for temperatures in the range 1000 - 20000 K, and charge transfer and (de)excitation processes involving the first excited S-states, 4s.5So and 4s.3So, are found to have the highest rates. △ Less

Submitted 5 December, 2017; v1 submitted 4 December, 2017; originally announced December 2017.

Comments: Accepted for A&A. Data will be made available at CDS. Is available here: https://github.com/barklem/public-data. Replaced version corrects url and adds it to the paper

Journal ref: A&A 610, A57 (2018)

Inelastic e+Mg collision data and its impact on modelling stellar and supernova spectra

Authors: P. S. Barklem, Y. Osorio, D. V. Fursa, I. Bray, O. Zatsarinny, K. Bartschat, A. Jerkstrand

Abstract: Results of calculations for inelastic e+Mg effective collision strengths for the lowest 25 physical states of Mg I (up to 3s6p 1P), and thus 300 transitions, from the convergent close-coupling (CCC) and the B-spline R-matrix (BSR) methods are presented. At temperatures of interest, ~5000 K, the results of the two calculations differ on average by only 4%, with a scatter of 27%. As the methods are… ▽ More Results of calculations for inelastic e+Mg effective collision strengths for the lowest 25 physical states of Mg I (up to 3s6p 1P), and thus 300 transitions, from the convergent close-coupling (CCC) and the B-spline R-matrix (BSR) methods are presented. At temperatures of interest, ~5000 K, the results of the two calculations differ on average by only 4%, with a scatter of 27%. As the methods are independent, this suggests that the calculations provide datasets for e+Mg collisions accurate to this level. Comparison with the commonly used dataset compiled by Mauas et al. (1988), covering 25 transitions among 12 states, suggests the Mauas et al. data are on average ~57% too low, and with a very large scatter of a factor of ~6.5. In particular the collision strength for the transition corresponding to the Mg I intercombination line at 457 nm is significantly underestimated by Mauas et al., which has consequences for models that employ this dataset. In giant stars the new data leads to a stronger line compared to previous non-LTE calculations, and thus a reduction in the non-LTE abundance correction by ~0.1 dex (~25%). A non-LTE calculation in a supernova ejecta model shows this line becomes significantly stronger, by a factor of around two, alleviating the discrepancy where the 457 nm line in typical models with Mg/O ratios close to solar tended to be too weak compared to observations. △ Less

Submitted 11 June, 2017; originally announced June 2017.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 606, A11 (2017)

Non-LTE line formation of Fe in late-type stars IV: Modelling of the solar centre-to-limb variation in 3D

Authors: K. Lind, A. M. Amarsi, M. Asplund, P. S. Barklem, M. Bautista, M. Bergemann, R. Collet, D. Kiselman, J. Leenaarts, T. M. D. Pereira

Abstract: Our ability to model the shapes and strengths of iron lines in the solar spectrum is a critical test of the accuracy of the solar iron abundance, which sets the absolute zero-point of all stellar metallicities. We use an extensive 463-level Fe atom with new photoionisation cross-sections for FeI as well as quantum mechanical calculations of collisional excitation and charge transfer with neutral h… ▽ More Our ability to model the shapes and strengths of iron lines in the solar spectrum is a critical test of the accuracy of the solar iron abundance, which sets the absolute zero-point of all stellar metallicities. We use an extensive 463-level Fe atom with new photoionisation cross-sections for FeI as well as quantum mechanical calculations of collisional excitation and charge transfer with neutral hydrogen; the latter effectively remove a free parameter that has hampered all previous line formation studies of Fe in non-local thermodynamic equilibrium (NLTE). For the first time, we use realistic 3D NLTE calculations of Fe for a quantitative comparison to solar observations. We confront our theoretical line profiles with observations taken at different viewing angles across the solar disk with the Swedish 1-m Solar Telescope. We find that 3D modelling well reproduces the observed centre-to-limb behaviour of spectral lines overall, but highlight aspects that may require further work, especially cross-sections for inelastic collisions with electrons. Our inferred solar iron abundance is log(eps(Fe))=7.48+-0.04. △ Less

Submitted 11 March, 2017; originally announced March 2017.

Comments: 13 pages, 8 figures, accepted for publication in MNRAS

3D NLTE analysis of the most iron-deficient star, SMSS0313-6708

Authors: Thomas Nordlander, Anish M. Amarsi, Karin Lind, Martin Asplund, Paul S. Barklem, Andy R. Casey, Remo Collet, Jorrit Leenaarts

Abstract: Models of star formation in the early universe depend on the details of accretion, fragmentation and radiative feedback. Different simulations predict different initial mass functions of the first stars, ranging from predominantly low mass (0.1-10 Msol), to massive (10-100 Msol), or even supermassive (100-1000 Msol). The mass distribution of the first stars should lead to unique chemical imprints… ▽ More Models of star formation in the early universe depend on the details of accretion, fragmentation and radiative feedback. Different simulations predict different initial mass functions of the first stars, ranging from predominantly low mass (0.1-10 Msol), to massive (10-100 Msol), or even supermassive (100-1000 Msol). The mass distribution of the first stars should lead to unique chemical imprints on the low-mass second and later generation metal-poor stars still in existence. The chemical composition of SMSS0313-6708, which has the lowest abundances of Ca and Fe of any star known, indicates it was enriched by a single massive supernova. However, even weak spectral lines may be affected by strong 3D and NLTE effects in metal-poor stars. If these effects are ignored or treated incorrectly, errors in the inferred abundances may significantly bias the inferred properties of the polluting supernovae. We redetermine the chemical composition of SMSS0313-6708 using 3D NLTE radiative transfer to obtain accurate abundances for Li, Na, Mg, Al, Ca and Fe. The model atoms employ realistic collisional rates, with no calibrated free parameters. We find significantly higher abundances in 3D NLTE than 1D LTE by 0.8 dex for Fe, and 0.5 dex for Mg, Al and Ca, while Li and Na are unaffected to within 0.03 dex. In particular, our upper limit for [Fe/H] is now a factor ten larger, at [Fe/H] < -6.53 (3 sigma), than previous estimates based on <3D> NLTE (i.e. using averaged 3D models). This higher estimate is due to a conservative upper limit estimation, updated NLTE data, and 3D-<3D> NLTE differences, all of which lead to a higher abundance determination. We find that the revised chemical composition of SMSS0313-6708 matches supernova yields for massive progenitors of 20-60 Msol exploding with low energies (1-2 x 10^51 erg), as well as progenitors of 10 Msol with very low explosion energies (< 10^51 erg). △ Less

Submitted 23 September, 2016; originally announced September 2016.

Comments: 17 pages, 10 pages; Accepted for publication in A&A

Journal ref: A&A 597, A6 (2017)

The Hamburg/ESO R-process Enhanced Star survey (HERES) XI. The highly $r$-process-enhanced star CS 29497-004

Authors: V. Hill, N. Christlieb, T. C. Beers, P. S. Barklem, K. -L. Kratz, B. Nordström, B. Pfeiffer, K. Farouqi

Abstract: We report an abundance analysis for the highly r-process-enhanced (r-II) star CS 29497-004, a very metal-poor giant with Teff = 5013K and [Fe/H]=-2.85, whose nature was initially discovered in the course of the HERES project. Our analysis is based on high signal-to-noise, high-resolution (R~75000) VLT/UVES spectra and MARCS model atmospheres under the assumption of local thermodynamic equilibrium,… ▽ More We report an abundance analysis for the highly r-process-enhanced (r-II) star CS 29497-004, a very metal-poor giant with Teff = 5013K and [Fe/H]=-2.85, whose nature was initially discovered in the course of the HERES project. Our analysis is based on high signal-to-noise, high-resolution (R~75000) VLT/UVES spectra and MARCS model atmospheres under the assumption of local thermodynamic equilibrium, and obtains abundance measurements for a total of 46 elements, 31 of which are neutron-capture elements. As is the case for the other 25 r-II stars currently known, the heavy-element abundance pattern of CS 29497-004 well-matches a scaled Solar System second peak r-process-element abundance pattern. We confirm our previous detection of Th, and demonstrate that this star does not exhibit an "actinide boost". Uranium is also detected (log e(U) =-2.20+/-0.30), albeit with a large measurement error that hampers its use as a precision cosmo-chronometer. Combining the various elemental chronometer pairs that are available for this star, we derive a mean age of 12.2+/-3.7 Gyr using the theoretical production ratios from waiting-point approximation models (Kratz et al. 2007). We further explore the high-entropy wind model (Farouqi et al. 2010) production ratios arising from different neutron richness of the ejecta (Y_e), and derive an age of 13.7+/-4.4 Gyr for a best-fitting Y_e=0.447. The U/Th nuclei-chronometer is confirmed to be the most resilient to theoretical production ratios and yields an age of 16.5+:-6.6 Gyr. Lead (Pb) is also tentatively detected in CS 29497-004, at a level compatible with a scaled Solar r-process, or with the theoretical expectations for a pure r-process in this star. △ Less

Submitted 26 August, 2016; originally announced August 2016.

Comments: 24 pages, 15 figures, accepted in Astronomy and Astrophysics

Journal ref: A&A 607, A91 (2017)

Non-LTE line formation of Fe in late-type stars - III. 3D non-LTE analysis of metal-poor stars

Authors: A. M. Amarsi, K. Lind, M. Asplund, P. S. Barklem, R. Collet

Abstract: As one of the most important elements in astronomy, iron abundance determinations need to be as accurate as possible. We investigate the accuracy of spectroscopic iron abundance analyses using archetypal metal-poor stars. We perform detailed 3D non-LTE radiative transfer calculations based on 3D hydrodynamic Stagger model atmospheres, and employ a new model atom that includes new quantum-mechanica… ▽ More As one of the most important elements in astronomy, iron abundance determinations need to be as accurate as possible. We investigate the accuracy of spectroscopic iron abundance analyses using archetypal metal-poor stars. We perform detailed 3D non-LTE radiative transfer calculations based on 3D hydrodynamic Stagger model atmospheres, and employ a new model atom that includes new quantum-mechanical neutral hydrogen collisional rate coefficients. With the exception of the red giant HD122563, we find that the 3D non-LTE models achieve Fe i/Fe ii excitation and ionization balance as well as not having any trends with equivalent width to within modelling uncertainties of 0.05 dex, all without having to invoke any microturbulent broadening; for HD122563 we predict that the current best parallax-based surface gravity is over-estimated by 0.5 dex. Using a 3D non-LTE analysis, we infer iron abundances from the 3D model atmospheres that are roughly 0.1 dex higher than corresponding abundances from 1D MARCS model atmospheres; these differences go in the same direction as the non-LTE effects themselves. We make available grids of departure coefficients, equivalent widths and abundance corrections, calculated on 1D MARCS model atmospheres and horizontally- and temporally-averaged 3D Stagger model atmospheres. △ Less

Submitted 13 September, 2016; v1 submitted 23 August, 2016; originally announced August 2016.

Comments: 18 pages, 5 figures; Accepted for publication in MNRAS

Accurate abundance analysis of late-type stars: advances in atomic physics

Authors: Paul S. Barklem

Abstract: The measurement of stellar properties such as chemical compositions, masses and ages, through stellar spectra, is a fundamental problem in astrophysics. Progress in the understanding, calculation and measurement of atomic properties and processes relevant to the high-accuracy analysis of F-, G-, and K-type stellar spectra is reviewed, with particular emphasis on abundance analysis. This includes f… ▽ More The measurement of stellar properties such as chemical compositions, masses and ages, through stellar spectra, is a fundamental problem in astrophysics. Progress in the understanding, calculation and measurement of atomic properties and processes relevant to the high-accuracy analysis of F-, G-, and K-type stellar spectra is reviewed, with particular emphasis on abundance analysis. This includes fundamental atomic data such as energy levels, wavelengths, and transition probabilities, as well as processes of photoionisation, collisional broadening and inelastic collisions. A recurring theme throughout the review is the interplay between theoretical atomic physics, laboratory measurements, and astrophysical modelling, all of which contribute to our understanding of atoms and atomic processes, as well as to modelling stellar spectra. △ Less

Submitted 9 May, 2016; v1 submitted 26 April, 2016; originally announced April 2016.

Comments: Review accepted for publication by The A&A Review (Springer), 9 Figures, 3 Tables, 56 pages

Excitation and charge transfer in low-energy hydrogen atom collisions with neutral atoms: Theory, comparisons, and application to Ca

Authors: Paul S. Barklem

Abstract: A theoretical method for the estimation of cross sections and rates for excitation and charge transfer processes in low-energy hydrogen atom collisions with neutral atoms, based on an asymptotic two-electron model of ionic-covalent interactions in the neutral atom-hydrogen atom system, is presented. The calculation of potentials and non-adiabatic radial couplings using the method is demonstrated.… ▽ More A theoretical method for the estimation of cross sections and rates for excitation and charge transfer processes in low-energy hydrogen atom collisions with neutral atoms, based on an asymptotic two-electron model of ionic-covalent interactions in the neutral atom-hydrogen atom system, is presented. The calculation of potentials and non-adiabatic radial couplings using the method is demonstrated. The potentials are used together with the multi-channel Landau-Zener model to calculate cross sections and rate coefficients. The main feature of the method is that it employs asymptotically exact atomic wavefunctions, which can be determined from known atomic parameters. The method is applied to Li+H, Na+H, and Mg+H collisions, and the results compare well with existing detailed full-quantum calculations. The method is applied to the astrophysically important problem of Ca+H collisions, and rate coefficients are calculated for temperatures in the range 1000-20000 K. △ Less

Submitted 20 April, 2016; v1 submitted 23 March, 2016; originally announced March 2016.

Comments: Accepted by Phys. Rev. A, small corrections made

Partition functions and equilibrium constants for diatomic molecules and atoms of astrophysical interest

Authors: Paul S. Barklem, Remo Collet

Abstract: Partition functions and dissociation equilibrium constants are presented for 291 diatomic molecules for temperatures in the range from near absolute zero to 10000 K, thus providing data for many diatomic molecules of astrophysical interest at low temperature. The calculations are based on molecular spectroscopic data from the book of Huber and Herzberg with significant improvements from the litera… ▽ More Partition functions and dissociation equilibrium constants are presented for 291 diatomic molecules for temperatures in the range from near absolute zero to 10000 K, thus providing data for many diatomic molecules of astrophysical interest at low temperature. The calculations are based on molecular spectroscopic data from the book of Huber and Herzberg with significant improvements from the literature, especially updated data for ground states of many of the most important molecules by Irikura. Dissociation energies are collated from compilations of experimental and theoretical values. Partition functions for 284 species of atoms for all elements from H to U are also presented based on data collected at NIST. The calculated data are expected to be useful for modelling a range of low density astrophysical environments, especially star-forming regions, protoplanetary disks, the interstellar medium, and planetary and cool stellar atmospheres. The input data, which will be made available electronically, also provides a possible foundation for future improvement by the community. △ Less

Submitted 10 February, 2016; originally announced February 2016.

Comments: 13 pages, 8 figures, 8 tables. Full tables 1, 2, 4, 5, 6, 7 and 8 to be made available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/

Journal ref: A&A 588, A96 (2016)

Mg line formation in late-type stellar atmospheres: II. Calculations in a grid of 1D models

Authors: Yeisson Osorio, Paul S. Barklem

Abstract: Mg is the alpha element of choice for Galactic population and chemical evolution studies as it is easily detectable in all late-type stars. Such studies require precise elemental abundances, and thus departures from LTE need to be accounted for. Our goal is to provide reliable departure coefficients and equivalent widths in non-LTE, and for reference in LTE, for diagnostic lines of Mg studied in l… ▽ More Mg is the alpha element of choice for Galactic population and chemical evolution studies as it is easily detectable in all late-type stars. Such studies require precise elemental abundances, and thus departures from LTE need to be accounted for. Our goal is to provide reliable departure coefficients and equivalent widths in non-LTE, and for reference in LTE, for diagnostic lines of Mg studied in late-type stars. These can be used e.g., to correct LTE spectra and abundances. Using the model atom built and tested in the preceding paper in this series, we performed non-LTE radiative transfer calculations in a grid of 3945 stellar 1D atmospheric models. We used a sub-grid of 86 models to explore the propagation of errors in the recent atomic collision calculations to the radiative transfer results. We obtained departure coefficients for all the levels and equivalent widths (in LTE and non-LTE) for all the radiative transitions included in the "final" model atom of Osorio et al.. We present and describe our results and show some examples of applications of the data. The errors due to uncertainties in the collisional data are investigated and tabulated. The results for equivalent widths and departure coefficients are made freely available. Giants tend to have negative abundance corrections while dwarfs have positive, though small, corrections. Error analysis results show that uncertainties related to the atomic collision data are typically of order 0.01 dex or less, although for few stellar models in specific lines uncertainties can be as large as 0.03 dex. As these errors are less than or of the same order as typical corrections, we expect that we can use these results to extract Mg abundances from high quality spectra more reliably than from classical LTE analysis. △ Less

Submitted 29 January, 2016; v1 submitted 17 October, 2015; originally announced October 2015.

Comments: 11 pages

Journal ref: A&A 586, A120 (2016)

A unified numerical model of collisional depolarization and broadening rates due to hydrogen atom collisions

Authors: M. Derouich, A. Radi, P. S. Barklem

Abstract: Interpretation of solar polarization spectra accounting for partial or complete frequency redistribution requires data on various collisional processes. Data for depolarization and polarization transfer are needed but often missing, while data for collisional broadening are usually more readily available. Recent work by Sahal-Bréchot and Bommier concluded that despite underlying similarities in th… ▽ More Interpretation of solar polarization spectra accounting for partial or complete frequency redistribution requires data on various collisional processes. Data for depolarization and polarization transfer are needed but often missing, while data for collisional broadening are usually more readily available. Recent work by Sahal-Bréchot and Bommier concluded that despite underlying similarities in the physics of collisional broadening and depolarization processes, relationships between them are not possible to derive purely analytically. We aim to derive accurate numerical relationships between the collisional broadening rates and the collisional depolarization and polarization transfer rates due to hydrogen atom collisions. Such relationships would enable accurate and efficient estimation of collisional data for solar applications. Using earlier results for broadening and depolarization processes based on general (i.e. not specific to a given atom), semi-classical calculations employing interaction potentials from perturbation theory, genetic programming (GP) has been used to fit the available data and generate analytical functions describing the relationships between them. The predicted relationships from the GP-based model are compared with the original data to estimate the accuracy of the method. △ Less

Submitted 26 August, 2015; originally announced August 2015.

Comments: 10 pages, 7 figures, Accepted for publication in Astronomy & Astrophysics

Atomic and Molecular Data for Optical Stellar Spectroscopy

Authors: U. Heiter, K. Lind, M. Asplund, P. S. Barklem, M. Bergemann, L. Magrini, T. Masseron, Š. Mikolaitis, J. C. Pickering, M. P. Ruffoni

Abstract: High-precision spectroscopy of large stellar samples plays a crucial role for several topical issues in astrophysics. Examples include studying the chemical structure and evolution of the Milky Way galaxy, tracing the origin of chemical elements, and characterizing planetary host stars. Data are accumulating from instruments that obtain high-quality spectra of stars in the ultraviolet, optical and… ▽ More High-precision spectroscopy of large stellar samples plays a crucial role for several topical issues in astrophysics. Examples include studying the chemical structure and evolution of the Milky Way galaxy, tracing the origin of chemical elements, and characterizing planetary host stars. Data are accumulating from instruments that obtain high-quality spectra of stars in the ultraviolet, optical and infrared wavelength regions on a routine basis. These instruments are located at ground-based 2- to 10-m class telescopes around the world, in addition to the spectrographs with unique capabilities available at the Hubble Space Telescope. The interpretation of these spectra requires high-quality transition data for numerous species, in particular neutral and singly ionized atoms, and di- or triatomic molecules. We rely heavily on the continuous efforts of laboratory astrophysics groups that produce and improve the relevant experimental and theoretical atomic and molecular data. The compilation of the best available data is facilitated by databases and electronic infrastructures such as the NIST Atomic Spectra Database, the VALD database, or the Virtual Atomic and Molecular Data Centre (VAMDC). We illustrate the current status of atomic data for optical stellar spectra with the example of the Gaia-ESO Public Spectroscopic Survey. Data sources for 35 chemical elements were reviewed in an effort to construct a line list for a homogeneous abundance analysis of up to 100000 stars. △ Less

Submitted 22 June, 2015; originally announced June 2015.

Comments: Published 30 April 2015 in Physica Scripta

Journal ref: Phys. Scr. 90 (2015) 054010

Mg line formation in late-type stellar atmospheres: I. The model atom

Authors: Y. Osorio, P. S. Barklem, K. Lind, A. K. Belyaev, A. Spielfiedel, M. Guitou, N. Feautrier

Abstract: Mg is often traced in late-type stars using lines of neutral magnesium, which is expected to be subject to departures from LTE. The astrophysical importance of Mg as well as its relative simplicity from an atomic physics point of view, makes it a prime target and test bed for detailed ab initio non-LTE modelling in stellar atmospheres. For the low-lying states of Mg i, electron collision data were… ▽ More Mg is often traced in late-type stars using lines of neutral magnesium, which is expected to be subject to departures from LTE. The astrophysical importance of Mg as well as its relative simplicity from an atomic physics point of view, makes it a prime target and test bed for detailed ab initio non-LTE modelling in stellar atmospheres. For the low-lying states of Mg i, electron collision data were calculated using the R-matrix method. Calculations for collisional broadening by neutral hydrogen were also performed where data were missing. These calculations, together with data from the literature, were used to build a model atom. First, the modelling was tested by comparisons with observed spectra of benchmark stars with well-known parameters. Second, the spectral line behaviour and uncertainties were explored by extensive experiments in which sets of collisional data were changed or removed. The modelled spectra agree well with observed spectra. The line-to-line scatter in the derived abundances shows improvements compared to LTE. The observed Mg emission features at 7 and 12 microns in the spectra of the Sun and Arcturus are reasonably well reproduced. The modelling predicts non-LTE abundance corrections in dwarfs, both solar metallicity and metal-poor, to be very small (of order 0.01 dex), even smaller than found in previous studies. In giants, corrections vary greatly between lines, but can be as large as 0.4 dex. Our results emphasise the need for accurate data of Mg collisions with both electrons and H atoms for precise non-LTE predictions of stellar spectra, but demonstrate that such data can be calculated and that ab initio non-LTE modelling without resort to free parameters is possible. Grids of departure coefficients and abundance corrections for a range of stellar parameters are planned for a forthcoming paper. △ Less

Submitted 28 April, 2015; originally announced April 2015.

Comments: 23 pages, 11 figures

The Gaia-ESO Survey: A globular cluster escapee in the Galactic halo

Authors: K. Lind, S. E. Koposov, C. Battistini, A. F. Marino, G. Ruchti, A. Serenelli, C. C. Worley, A. Alves-Brito, M. Asplund, P. S. Barklem, T. Bensby, M. Bergemann, S. Blanco-Cuaresma, A. Bragaglia, B. Edvardsson, S. Feltzing, P. Gruyters, U. Heiter, P. Jofre, A. J. Korn, T. Nordlander, N. Ryde, C. Soubiran, G. Gilmore, S. Randich , et al. (15 additional authors not shown)

Abstract: A small fraction of the halo field is made up of stars that share the light element (Z<=13) anomalies characteristic of second generation globular cluster (GC) stars. The ejected stars shed light on the formation of the Galactic halo by tracing the dynamical history of the clusters, which are believed to have once been more massive. Some of these ejected stars are expected to show strong Al enhanc… ▽ More A small fraction of the halo field is made up of stars that share the light element (Z<=13) anomalies characteristic of second generation globular cluster (GC) stars. The ejected stars shed light on the formation of the Galactic halo by tracing the dynamical history of the clusters, which are believed to have once been more massive. Some of these ejected stars are expected to show strong Al enhancement at the expense of shortage of Mg, but until now no such star has been found. We search for outliers in the Mg and Al abundances of the few hundreds of halo field stars observed in the first eighteen months of the Gaia-ESO public spectroscopic survey. One halo star at the base of the red giant branch, here referred to as 22593757-4648029 is found to have [Mg/Fe]=-0.36+-0.04 and [Al/Fe]=0.99+-0.08, which is compatible with the most extreme ratios detected in GCs so far. We compare the orbit of 22593757-4648029 to GCs of similar metallicity and find it unlikely that this star has been tidally stripped with low ejection velocity from any of the clusters. However, both chemical and kinematic arguments render it plausible that the star has been ejected at high velocity from the anomalous GC omega Centauri within the last few billion years. We cannot rule out other progenitor GCs, because some may have disrupted fully, and the abundance and orbital data are inadequate for many of those that are still intact. △ Less

Submitted 13 February, 2015; originally announced February 2015.

Comments: 4 pages, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 575, L12 (2015)

Atomic diffusion and mixing in old stars IV: Weak abundance trends in the globular cluster NGC 6752

Authors: Pieter Gruyters, Andreas J. Korn, Olivier Richard, Frank Grundahl, Remo Collet, Lyudmila I. Mashonkina, Yeisson Osorio, Paul S. Barklem

Abstract: Atomic diffusion in stars can create systematic trends of surface abundances with evolutionary stage. Globular clusters offer useful laboratories to put observational constraints on this theory as one needs to compare abundances in unevolved and evolved stars, all drawn from the same stellar population. In this paper, we show the results of an abundance study of stars in the globular cluster NGC67… ▽ More Atomic diffusion in stars can create systematic trends of surface abundances with evolutionary stage. Globular clusters offer useful laboratories to put observational constraints on this theory as one needs to compare abundances in unevolved and evolved stars, all drawn from the same stellar population. In this paper, we show the results of an abundance study of stars in the globular cluster NGC6752 which shows weak but systematic abundances trends with evolutionary phase for Fe, Sc, Ti and Ca. The trends are best explained by a stellar structure model including atomic diffusion with efficient additional mixing. The model allows to correct for sub-primordial stellar lithium abundances of the stars on the Spite plateau, and to match it to the WMAP-calibrated Big-Bang nucleosynthesis predictions to within the mutual 1-sigma errors. △ Less

Submitted 10 May, 2013; v1 submitted 8 May, 2013; originally announced May 2013.

Comments: 15 pages, 4 figures and 8 tables

The Most Metal-Poor Stars. IV. The Two Populations With [Fe/H] < -3.0

Authors: John E. Norris, David Yong, M. S. Bessell, N. Christlieb, M. Asplund, Gerard Gilmore, Rosemary F. G. Wyse, Timothy C. Beers, P. S. Barklem, Anna Frebel, S. G. Ryan

Abstract: We discuss the carbon-normal and carbon-rich populations of Galactic halo stars having [Fe/H] < -3.0, utilizing chemical abundances from high-resolution, high-S/N model-atmosphere analyses. The C-rich population represents ~28% of stars below [Fe/H] = -3.1, with the present C-rich sample comprising 16 CEMP-no stars, and two others with [Fe/H] ~ -5.5 and uncertain classification. The population is… ▽ More We discuss the carbon-normal and carbon-rich populations of Galactic halo stars having [Fe/H] < -3.0, utilizing chemical abundances from high-resolution, high-S/N model-atmosphere analyses. The C-rich population represents ~28% of stars below [Fe/H] = -3.1, with the present C-rich sample comprising 16 CEMP-no stars, and two others with [Fe/H] ~ -5.5 and uncertain classification. The population is O-rich ([O/Fe] > +1.5); the light elements Na, Mg, and Al are enhanced relative to Fe in half the sample; and for Z > 20 (Ca) there is little evidence for enhancements relative to solar values. These results are best explained in terms of the admixing and processing of material from H-burning and He-burning regions as achieved by nucleosynthesis in zero-heavy-element models in the literature of "mixing and fallback" supernovae (SNe); of rotating, massive and intermediate mass stars; and of Type II SNe with relativistic jets. The available (limited) radial velocities offer little support for the C-rich stars with [Fe/H] < -3.1 being binary. More data are required before one could conclude that binarity is key to an understanding of this population. We suggest that the C-rich and C-normal populations result from two different gas cooling channels in the very early Universe, of material that formed the progenitors of the two populations. The first was cooling by fine-structure line transitions of CII and OI (to form the C-rich population); the second, while not well-defined (perhaps dust-induced cooling?), led to the C-normal group. In this scenario, the C-rich population contains the oldest stars currently observed. △ Less

Submitted 13 November, 2012; originally announced November 2012.

Comments: Accepted for publication in ApJ

The Most Metal-Poor Stars. III. The Metallicity Distribution Function and CEMP Fraction

Authors: David Yong, John E. Norris, M. S. Bessell, N. Christlieb, M. Asplund, Timothy C. Beers, P. S. Barklem, Anna Frebel, S. G. Ryan

Abstract: We examine the metallicity distribution function (MDF) and fraction of carbon-enhanced metal-poor (CEMP) stars in a sample that includes 86 stars with [Fe/H] < -3.0, based on high-resolution, high-S/N spectroscopy, of which some 32 objects lie below [Fe/H] = -3.5. After accounting for the completeness function, the "corrected" MDF does not exhibit the sudden drop at [Fe/H] = -3.6 that was found in… ▽ More We examine the metallicity distribution function (MDF) and fraction of carbon-enhanced metal-poor (CEMP) stars in a sample that includes 86 stars with [Fe/H] < -3.0, based on high-resolution, high-S/N spectroscopy, of which some 32 objects lie below [Fe/H] = -3.5. After accounting for the completeness function, the "corrected" MDF does not exhibit the sudden drop at [Fe/H] = -3.6 that was found in recent samples of dwarfs and giants from the Hamburg/ESO survey. Rather, the MDF decreases smoothly down to [Fe/H] = -4.1. Similar results are obtained from the "raw" MDF. We find the fraction of CEMP objects below [Fe/H] = -3.0 is 23 +/- 6% and 32 +/- 8% when adopting the Beers et al. and Aoki et al. CEMP definitions, respectively. The former value is in fair agreement with some previous measurements, which adopt the Beers et al. criterion. △ Less

Submitted 14 August, 2012; originally announced August 2012.

Comments: Accepted for publication in ApJ

The Most Metal-Poor Stars. II. Chemical Abundances of 190 Metal-Poor Stars Including 10 New Stars With [Fe/H] < -3.5

Authors: David Yong, John E. Norris, M. S. Bessell, N. Christlieb, M. Asplund, Timothy C. Beers, P. S. Barklem, Anna Frebel, S. G. Ryan

Abstract: We present a homogeneous chemical abundance analysis of 16 elements in 190 metal-poor Galactic halo stars (38 program and 152 literature objects). The sample includes 171 stars with [Fe/H] < -2.5, of which 86 are extremely metal poor, [Fe/H] < -3.0. Our program stars include ten new objects with [Fe/H] < -3.5. We identify a sample of "normal" metal-poor stars and measure the trends between [X/Fe]… ▽ More We present a homogeneous chemical abundance analysis of 16 elements in 190 metal-poor Galactic halo stars (38 program and 152 literature objects). The sample includes 171 stars with [Fe/H] < -2.5, of which 86 are extremely metal poor, [Fe/H] < -3.0. Our program stars include ten new objects with [Fe/H] < -3.5. We identify a sample of "normal" metal-poor stars and measure the trends between [X/Fe] and [Fe/H], as well as the dispersion about the mean trend for this sample. Using this mean trend, we identify objects that are chemically peculiar relative to "normal" stars at the same metallicity. These chemically unusual stars include CEMP-no objects, one star with high [Si/Fe], another with high [Ba/Sr], and one with unusually low [X/Fe] for all elements heavier than Na. The Sr and Ba abundances indicate that there may be two nucleosynthetic processes at lowest metallicity that are distinct from the main r-process. Finally, for many elements, we find a significant trend between [X/Fe] versus Teff which likely reflects non-LTE and/or 3D effects. Such trends demonstrate that care must be exercised when using abundance measurements in metal-poor stars to constrain chemical evolution and/or nucleosynthesis predictions. △ Less

Submitted 14 August, 2012; originally announced August 2012.

Comments: Accepted for publication in ApJ

The Most Metal-Poor Stars. I. Discovery, Data, and Atmospheric Parameters

Authors: John E. Norris, M. S. Bessell, David Yong, N. Christlieb, P. S. Barklem, M. Asplund, Simon J. Murphy, Timothy C. Beers, Anna Frebel, S. G. Ryan

Abstract: We report the discovery of 34 stars in the Hamburg/ESO Survey for metal-poor stars and the Sloan Digital Sky Survey that have [Fe/H] < -3.0. Their median and minimum abundances are [Fe/H] = -3.1 and -4.1, respectively, while 10 stars have [Fe/H] < -3.5. High-resolution, high-S/N spectroscopic data - equivalent widths and radial velocities - are presented for these stars, together with an additiona… ▽ More We report the discovery of 34 stars in the Hamburg/ESO Survey for metal-poor stars and the Sloan Digital Sky Survey that have [Fe/H] < -3.0. Their median and minimum abundances are [Fe/H] = -3.1 and -4.1, respectively, while 10 stars have [Fe/H] < -3.5. High-resolution, high-S/N spectroscopic data - equivalent widths and radial velocities - are presented for these stars, together with an additional four objects previously reported or currently being investigated elsewhere. We have determined the atmospheric parameters, effective temperature (Teff) and surface gravity (logg), which are critical in the determination of the chemical abundances and the evolutionary status of these stars. Three techniques were used to derive these parameters. Spectrophotometric fits to model atmosphere fluxes were used to derive Teff, logg, and an estimate of E(B-V); Halpha, Hbeta, and Hgamma profile fitting to model atmosphere results provided the second determination of Teff and logg; and finally, we used an empirical Teff-calibrated Hdelta index, for the third, independent Teff determination. The three values of Teff are in good agreement, although the profile fitting may yield systematically cooler Teff values, by ~100K. This collective data set will be analyzed in future papers in the present series to utilize the most metal-poor stars as probes of conditions in the early Universe. △ Less

Submitted 14 August, 2012; originally announced August 2012.

Comments: Accepted for publication in ApJ

Precise Effective Temperatures of Solar Analog Stars

Authors: D. Cornejo, I. Ramirez, P. S. Barklem

Abstract: We perform a study of 62 solar analog stars to compute their effective temperatures (Teff) using the Balmer line wing fitting procedure and compare them with Teff values obtained using other commonly employed methods. We use observed H-alpha spectral lines and a fine grid of theoretical LTE model spectra calculated with the best available atomic data and most recent quantum theory. Our spectroscop… ▽ More We perform a study of 62 solar analog stars to compute their effective temperatures (Teff) using the Balmer line wing fitting procedure and compare them with Teff values obtained using other commonly employed methods. We use observed H-alpha spectral lines and a fine grid of theoretical LTE model spectra calculated with the best available atomic data and most recent quantum theory. Our spectroscopic data are of very high quality and have been carefully normalized to recover the proper shape of the H-alpha line profile. We obtain Teff values with internal errors of about 25 K. Comparison of our results with those from other methods shows reasonably good agreement. Then, combining Teff values obtained from four independent techniques, we are able to determine final Teff values with errors of about 10 K. △ Less

Submitted 4 June, 2012; originally announced June 2012.

Comments: 4 pages, 4 figures, Comparative Magnetic Minima: Charactering proceedings of the IAU Symposium 286 on Comparative Magnetic Minima: Characterizing quiet times in the Sun and Stars, held in Mendoza, Argentina, October 3-7, 2011

Inelastic Mg+H collision data for non-LTE applications in stellar atmospheres

Authors: P. S. Barklem, A. K. Belyaev, A. Spielfiedel, M. Guitou, N. Feautrier

Abstract: Rate coefficients for inelastic Mg+H collisions are calculated for all transitions between the lowest seven levels and the ionic state (charge transfer). The rate coefficients are based on cross-sections from full quantum scattering calculations, which are themselves based on detailed quantum chemical calculations for the MgH molecule. The data are needed for non-LTE applications in cool astrophys… ▽ More Rate coefficients for inelastic Mg+H collisions are calculated for all transitions between the lowest seven levels and the ionic state (charge transfer). The rate coefficients are based on cross-sections from full quantum scattering calculations, which are themselves based on detailed quantum chemical calculations for the MgH molecule. The data are needed for non-LTE applications in cool astrophysical environments, especially cool stellar atmospheres, and are presented for a temperature range of 500-8000 K. From consideration of the sensitivity of the cross-sections to various uncertainties in the calculations, most importantly input quantum chemical data and the numerical accuracy of the scattering calculations, a measure of the possible uncertainties in the rate coefficients is estimated. △ Less

Submitted 21 March, 2012; originally announced March 2012.

Comments: Accepted by A&A. 5 pages, 1 figure

On inelastic hydrogen atom collisions in stellar atmospheres

Authors: P. S. Barklem, A. K. Belyaev, M. Guitou, N. Feautrier, F. X. Gadea, A. Spielfiedel

Abstract: The influence of inelastic hydrogen atom collisions on non-LTE spectral line formation has been, and remains to be, a significant source of uncertainty for stellar abundance analyses, due to the difficulty in obtaining accurate data for low-energy atomic collisions either experimentally or theoretically. For lack of a better alternative, the classical "Drawin formula" is often used. Over recent de… ▽ More The influence of inelastic hydrogen atom collisions on non-LTE spectral line formation has been, and remains to be, a significant source of uncertainty for stellar abundance analyses, due to the difficulty in obtaining accurate data for low-energy atomic collisions either experimentally or theoretically. For lack of a better alternative, the classical "Drawin formula" is often used. Over recent decades, our understanding of these collisions has improved markedly, predominantly through a number of detailed quantum mechanical calculations. In this paper, the Drawin formula is compared with the quantum mechanical calculations both in terms of the underlying physics and the resulting rate coefficients. It is shown that the Drawin formula does not contain the essential physics behind direct excitation by H atom collisions, the important physical mechanism being quantum mechanical in character. Quantitatively, the Drawin formula compares poorly with the results of the available quantum mechanical calculations, usually significantly overestimating the collision rates by amounts that vary markedly between transitions. △ Less

Submitted 7 April, 2011; originally announced April 2011.

Comments: 9 pages, 6 figures, accepted for A&A

Non-LTE calculations for neutral Na in late-type stars using improved atomic data

Authors: K. Lind, M. Asplund, P. S. Barklem, A. K. Belyaev

Abstract: Neutral sodium is a minority species in the atmospheres of late-type stars, and line formation in local thermodynamic equilibrium (LTE) is often a poor assumption, in particular for strong lines. We present an extensive grid of non-LTE calculations for several NaI lines in cool stellar atmospheres, including metal-rich and metal-poor dwarfs and giants. For the first time, we constructed a Na model… ▽ More Neutral sodium is a minority species in the atmospheres of late-type stars, and line formation in local thermodynamic equilibrium (LTE) is often a poor assumption, in particular for strong lines. We present an extensive grid of non-LTE calculations for several NaI lines in cool stellar atmospheres, including metal-rich and metal-poor dwarfs and giants. For the first time, we constructed a Na model atom that incorporates accurate quantum mechanical calculations for collisional excitation and ionisation by electrons as well as collisional excitation and charge exchange reactions with neutral hydrogen. Similar to LiI, the new rates for hydrogen impact excitation do not affect the statistical equilibrium calculations, while charge exchange reactions have a small but non-negligible influence. The presented LTE and non-LTE curves-of-growth can be interpolated to obtain non-LTE abundances and abundance corrections for arbitrary stellar parameter combinations and line strengths. The typical corrections for weak lines are -0.1...-0.2dex, whereas saturated lines may overestimate the abundance in LTE by more than 0.5dex. The non-LTE Na abundances appear very robust with respect to uncertainties in the input collisional data. △ Less

Submitted 10 February, 2011; originally announced February 2011.

Comments: 9 pages, 8 figures, accepted for publication in A&A