The spectral analysis and data products in Data Release 16 (DR16; December 2019) from the high-resolution near-infrared APOGEE-2/SDSS-IV survey are described. Compared to the previous APOGEE data release (DR14; July 2017), APOGEE DR16 includes about 200000 new stellar spectra, of which 100000 are from a new southern APOGEE instrument mounted on the 2.5 m du Pont telescope at Las Campanas Observatory in Chile. DR16 includes all data taken up to August 2018, including data released in previous data releases. All of the data have been re-reduced and re-analyzed using the latest pipelines, resulting in a total of 473307 spectra of 437445 stars. Changes to the analysis methods for this release include, but are not limited to, the use of MARCS model atmospheres for calculation of the entire main grid of synthetic spectra used in the analysis, a new method for filling "holes" in the grids due to unconverged model atmospheres, and a new scheme for continuum normalization. Abundances of the neutron capture element Ce are included for the first time. A new scheme for estimating uncertainties of the derived quantities using stars with multiple observations has been applied, and calibrated values of surface gravities for dwarf stars are now supplied. Compared to DR14, the radial velocities derived for this release more closely match those in the Gaia DR2 data base, and a clear improvement in the spectral analysis of the coolest giants can be seen. The reduced spectra as well as the result of the analysis can be downloaded using links provided in the SDSS DR16 web page.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

The new version of TLUSTY allows for the calculation of restricted NLTE in cool stars using pre-calculated opacity tables. We demonstrate that TLUSTY gives consistent results with MULTI, a well-tested code for NLTE in cool stars. We use TLUSTY to perform LTE and a series of NLTE calculations of Na, Mg, K and Ca using all combinations of 1, 2, 3 and the 4 elements mentioned above simultaneously in NLTE. In this work we take into account how departures from LTE in one element can affect others through changes in the opacities due to NLTE. We find that atomic Mg, which provides strong UV opacity, and shows departures from LTE in the low-energy states, can impact the NLTE populations of Ca, leading to abundance corrections as large as 0.07 dex. The differences in the derived abundances between the single-element and the multi-element cases can exceed those between the single-element NLTE and an LTE analysis, warning that this is not always a second-order effect. By means of detailed tests for three stars with reliable atmospheric parameters (Arcturus, Procyon and the Sun) we conclude that our NLTE calculations provide abundance corrections in the optical up to 0.1, 0.2 and 0.7 dex for Ca, Na and K, but LTE is a good approximation for Mg. In the H-band, NLTE corrections are much smaller, under 0.1 dex. The derived NLTE abundances in the optical and in the IR are consistent. For all four elements, in all three stars, NLTE line profiles fit better the observations than the LTE counterparts. For elements where over-ionisation is an important NLTE mechanism are likely affected by departures from LTE in Mg . Special care must be taken with the collisions adopted for high-lying levels when calculating NLTE profiles of lines in the H-band. The derived NLTE corrections in the optical and in the H-band differ, but the derived NLTE abundances are consistent between the two spectral regions.

The orbits of the least chemically enriched stars open a window on the formation of our Galaxy when it was still in its infancy. The common picture is that these low-metallicity stars are distributed as an isotropic, pressure-supported component since these stars were either accreted from the early building blocks of the assembling Milky Way, or were later brought by the accretion of faint dwarf galaxies. Combining the metallicities and radial velocities from the Pristine and LAMOST surveys and Gaia DR2 parallaxes and proper motions for an unprecedented large and unbiased sample of very metal-poor stars at $[Fe/H]\leq-2.5$ we show that this picture is incomplete. This sample shows strong statistical evidence (at the $5.0\sigma$ level) of asymmetry in their kinematics, favouring prograde motion. Moreover, we find that $31\%$ of the stars that currently reside in the disk do not venture outside of the disk plane throughout their orbit. The discovery of this population implies that a significant fraction of stars with iron abundances $[Fe/H]\leq-2.5$ formed within or concurrently with the Milky Way disk and that the history of the disk was quiet enough to allow them to retain their disk-like orbital properties.

We present the results of a 3-year long, medium-resolution spectroscopic campaign aimed at identifying very metal-poor stars from candidates selected with the CaHK, metallicity-sensitive Pristine survey. The catalogue consists of a total of 1007 stars, and includes 146 rediscoveries of metal-poor stars already presented in previous surveys, 707 new very metal-poor stars with [Fe/H] < -2.0, and 95 new extremely metal-poor stars with [Fe/H] < -3.0. We provide a spectroscopic [Fe/H] for every star in the catalogue, and [C/Fe] measurements for a subset of the stars (10% with [Fe/H] < -3 and 24% with -3 < [Fe/H] < -2) for which a carbon determination is possible, contingent mainly on the carbon abundance, effective temperature and S/N of the stellar spectra. We find an average carbon enhancement fraction ([C/Fe] >= +0.7) of 41 +- 4% for stars with -3 < [Fe/H] < -2 and 58 +- 14% for stars with [Fe/H] < -3, and report updated success rates for the Pristine survey of 56 % and 23 % to recover stars with [Fe/H] < -2.5 and [Fe/H] < -3, respectively. Finally, we discuss the current status of the survey and its preparation for providing targets to upcoming multi-object spectroscopic surveys such as WEAVE.

We previously explored the circumstellar effects on the Rb and Zr abundances in massive Galactic O-rich AGB stars. Here we are interested in the role of the extended atmosphere in the case of Li and Ca. Li is an important indicator of HBB while the total Ca in these stars could be affected by neutron captures. The Li abundances were previously studied with hydrostatic models, while the Ca abundances have been determined for the first time. We use a modified version of Turbospectrum and consider the presence of a gaseous circumstellar envelope and radial wind. The new Li abundances derived with the pseudo-dynamical models are very similar to those obtained from hydrostatic models (the average difference is 0.18 dex), with no difference for Ca. The Li and Ca content in these stars is only slightly affected by the presence of a circumstellar envelope. We also found that the Li I and Ca I line profiles are not very sensitive to variations of the model wind parameters. The new Li abundances confirm the Li-rich nature of the sample stars, supporting the activation of HBB in massive Galactic AGB stars. This is in good agreement with the theoretical predictions for solar metallicity AGB models from ATON, Monash, and NuGrid/MESA but is at odds with the FRUITY database, which predicts no HBB leading to the production of Li. Most sample stars display nearly solar Ca abundances that are consistent with the available s-process nucleosynthesis models for solar metallicity massive AGB stars, which predict overproduction of 46Ca relatively to the other Ca isotope and the creation of the radiactive isotope 41Ca but no change in the total Ca abundance. A minority of the stars seem to show a significant Ca depletion (by up to 1.0 dex). Possible explanations are offered to explain their apparent and unexpected Ca depletion.

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.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

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.

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.

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.

We study the Kepler object Kepler-432, an evolved star ascending the red giant branch. By deriving precise radial velocities from multi-epoch high-resolution spectra of Kepler-432 taken with the CAFE spectrograph at the 2.2m telescope of Calar Alto Observatory and the FIES spectrograph at the Nordic Optical Telescope of Roque de Los Muchachos Observatory, we confirm the planetary nature of the object Kepler-432 b, which has a transit period of 52 days. We find a planetary mass of Mp=5.84 +\- 0.05 Mjup and a high eccentricity of e=0.478 +\- 0.004. With a semi-major axis of a=0.303 +\- 0.007 AU, Kepler-432 b is the first bona fide warm Jupiter detected to transit a giant star. We also find a radial velocity linear trend of 0.44 +\- 0.04 m s$^{-1}$ d$^{-1}$, which suggests the presence of a third object in the system. Current models of planetary evolution in the post-main-sequence phase predict that Kepler-432 b will be most likely engulfed by its host star before the latter reaches the tip of the red giant branch.

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.

To date, the CoRoT space mission has produced more than 124,471 light curves. Classifying these curves in terms of unambiguous variability behavior is mandatory for obtaining an unbiased statistical view on their controlling root-causes. The present study provides an overview of semi-sinusoidal light curves observed by the CoRoT exo-field CCDs. We selected a sample of 4,206 light curves presenting well-defined semi-sinusoidal signatures. The variability periods were computed based on Lomb-Scargle periodograms, harmonic fits, and visual inspection. Color-period diagrams for the present sample show the trend of an increase of the variability periods as long as the stars evolve. This evolutionary behavior is also noticed when comparing the period distribution in the Galactic center and anti-center directions. These aspects indicate a compatibility with stellar rotation, although more information is needed to confirm their root-causes. Considering this possibility, we identified a subset of three Sun-like candidates by their photometric period. Finally, the variability period versus color diagram behavior was found to be highly dependent on the reddening correction.

The influence of the uncertainties in the rate coefficient data for electron-impact excitation and ionization on non-LTE Li line formation in cool stellar atmospheres is investigated. We examine the electron collision data used in previous non-LTE calculations and compare them to recent calculations that use convergent close-coupling (CCC) techniques and to our own calculations using the R-matrix with pseudostates (RMPS) method. We find excellent agreement between rate coefficients from the CCC and RMPS calculations, and reasonable agreement between these data and the semi-empirical data used in non-LTE calculations up to now. The results of non-LTE calculations using the old and new data sets are compared and only small differences found: about 0.01 dex (~ 2%) or less in the abundance corrections. We therefore conclude that the influence on non-LTE calculations of uncertainties in the electron collision data is negligible. Indeed, together with the collision data for the charge exchange process Li(3s) + H <-> Li^+ + H^- now available, and barring the existence of an unknown important collisional process, the collisional data in general is not a source of significant uncertainty in non-LTE Li line formation calculations.