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.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We analyze a large sample of 885 GCs giants from the APOGEE survey. We used the Cannon results to separate the red giant branch and the asymptotic giant branch stars, not only allowing for a refinement of surface gravity from isochrones, but also providing an independent H-band spectroscopic method to distinguish stellar evolutionary status in clusters. We then use the BACCHUS code to derive metallicity, microturbulence, acroturbulence and many light-element abundances as well as the neutron-capture elements Nd and Ce for the first time from the APOGEE GCs data. Our independent analysis helped us to diagnose issues regarding the standard analysis of the APOGEE DR14 for low-metallicity GC stars. Furthermore, while we confirm most of the known correlations and anti-correlation trends (Na-O, Mg-Al,C-N), we discover that some stars within our most metal-poor clusters show an extreme Mg depletion and some Si enhancement but at the same time show some relative Al depletion, displaying a turnover in the Mg-Al diagram. These stars suggest that Al has been partially depleted in their progenitors by very hot proton-capture nucleosynthetic processes. Furthermore, we attempted to quantitatively correlate the spread of Al abundances with the global properties of GCs. We find an anti-correlation of the Al spread against clusters metallicity and luminosity, but the data do not allow to find clear evidence of a dependence of N against metallicity in the more metal-poor clusters. Large and homogeneously analyzed samples from on-going spectroscopic surveys unveil unseen chemical details for many clusters, including a turnover in the Mg-Al anti-correlation, thus yielding new constrains for GCs formation evolution models.

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.

Data and analysis methodology used for the SDSS/APOGEE Data Releases 13 and 14 are described, highlighting differences from the DR12 analysis presented in Holtzman (2015). Some improvement in the handling of telluric absorption and persistence is demonstrated. The derivation and calibration of stellar parameters, chemical abundances, and respective uncertainties are described, along with the ranges over which calibration was performed. Some known issues with the public data related to the calibration of the effective temperatures (DR13), surface gravity (DR13 and DR14), and C and N abundances for dwarfs (DR13 and DR14) are highlighted. We discuss how results from a data-driven technique, The Cannon (Casey 2016), are included in DR14, and compare those with results from the APOGEE Stellar Parameters and Chemical Abundances Pipeline (ASPCAP). We describe how using The Cannon in a mode that restricts the abundance analysis of each element to regions of the spectrum with known features from that element leads to Cannon abundances can lead to significantly different results for some elements than when all regions of the spectrum are used to derive abundances.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

The Gaia-ESO Survey (GES) is a large public spectroscopic survey at the European Southern Observatory Very Large Telescope. A key aim is to provide precise radial velocities (RVs) and projected equatorial velocities (v sin i) for representative samples of Galactic stars, that will complement information obtained by the Gaia astrometry satellite. We present an analysis to empirically quantify the size and distribution of uncertainties in RV and v sin i using spectra from repeated exposures of the same stars. We show that the uncertainties vary as simple scaling functions of signal-to-noise ratio (S/N) and v sin i, that the uncertainties become larger with increasing photospheric temperature, but that the dependence on stellar gravity, metallicity and age is weak. The underlying uncertainty distributions have extended tails that are better represented by Student's t-distributions than by normal distributions. Parametrised results are provided, that enable estimates of the RV precision for almost all GES measurements, and estimates of the v sin i precision for stars in young clusters, as a function of S/N, v sin i and stellar temperature. The precision of individual high S/N GES RV measurements is 0.22-0.26 km/s, dependent on instrumental configuration.

The Sloan Digital Sky Survey--III (SDSS--III) Apache Point Observatory Galactic Evolution Experiment (APOGEE) has obtained high resolution (R $\sim$ 22,500), high signal-to-noise ratio ($>$ 100) spectra in the H$-$band ($\sim$1.5$-$1.7 $\mu$m) for about 146,000 stars in the Milky Way galaxy. We have computed spectral libraries with effective temperature ($T\rm{_{eff}}$) ranging from 3500 to 8000 K for the automated chemical analy\-sis of the survey data. The libraries, used to derive stellar parameters and abundances from the APOGEE spectra in the SDSS--III data release 12 (DR12), are based on ATLAS9 model atmospheres and the ASS$\epsilon$T spectral synthesis code. We present a second set of libraries based on MARCS model atmospheres and the spectral synthesis code Turbospectrum. The ATLAS9/ASS$\epsilon$T ($T\rm{_{eff}}$ = 3500$-$8000 K) and MARCS/Turbospectrum ($T\rm{_{eff}}$ = 3500$-$5500 K) grids cover a wide range of metallicity ($-$2.5 $\leq$ [M/H] $\leq$ $+$0.5 dex), surface gravity (0 $\leq$ log $g$ $\leq$ 5 dex), microturbulence (0.5 $\leq$ $\xi$ $\leq$ 8 km~s$^{-1}$), carbon ($-$1 $\leq$ [C/M] $\leq$ $+$1 dex), nitrogen ($-$1 $\leq$ [N/M] $\leq$ $+$1 dex), and $\alpha$-element ($-$1 $\leq$ [$\alpha$/M] $\leq$ $+$1 dex) variations, having thus seven dimensions. We compare the ATLAS9/ASS$\epsilon$T and MARCS/Turbospectrum libraries and apply both of them to the analysis of the observed H$-$band spectra of the Sun and the K2 giant Arcturus, as well as to a selected sample of well-known giant stars observed at very high-resolution. The new APOGEE libraries are publicly available and can be employed for chemical studies in the H$-$band using other high-resolution spectrographs.

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.

This paper describes the analysis of UVES and GIRAFFE spectra acquired by the Gaia-ESO Public Spectroscopic Survey in the fields of young clusters whose population includes pre-main sequence (PMS) stars. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia-ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing the results obtained by such different methods, while the accuracy is estimated by comparison with independent external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. Specific strategies are implemented to deal with fast rotation, accretion signatures, chromospheric activity, and veiling. The analysis carried out on spectra acquired in young clusters' fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be $\approx$ 120 K r.m.s. in Teff, $\approx$0.3 dex r.m.s. in logg, and $\approx$0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups.

Canonical grids of stellar evolutionary sequences have been computed for the helium mass-fraction abundances $Y = 0.25$, $0.29$, and $0.33$, and for iron abundances that vary from $-2.4$ to $+0.4$ (in 0.2 dex increments) when [$\alpha$/Fe] $= +0.4$, or for the ranges $-2.0 \le$ [Fe/H] $\le +0.6$, $-1.8 \le$ [Fe/H] $\le +0.6$ when [$\alpha$/Fe] $= 0.0$ and $-0.4$, respectively. The grids, which consist of tracks for masses from $0.12 {{\cal M}_\odot}$ to $1.1$-$1.5 {{\cal M}_\odot}$ (depending on the metallicity) are based on up-to-date physics, including the gravitational settling of helium (but not metals diffusion). Interpolation software is provided to generate isochrones for arbitrary ages between $\approx 5$ and $15$ Gyr and any values of $Y$, [$\alpha$/Fe], and [Fe/H] within the aforementioned ranges. Comparisons of isochrones with published color-magnitude diagrams (CMDs) for the open clusters M67 ([Fe/H] $\approx 0.0$) and NGC 6791 ([Fe/H] $\approx 0.3$) and for four of the metal-poor globular clusters (47 Tuc, M3, M5, and M92) indicate that the models for the observed metallicities do a reasonably good job of reproducing the locations and slopes of the cluster main sequences and giant branches. The same conclusion is reached from a consideration of plots of nearby subdwarfs that have accurate Hipparcos parallaxes and metallicities in the range $-2.0 \le$ [Fe/H] $\le -1.0$ on various CMDs and on the ($\log T_{\rm eff},\,M_V$)-diagram. A relatively hot temperature scale similar to that derived in recent calibrations of the infrared flux method is favored by both the isochrones and the adopted color transformations, which are based on the latest MARCS model atmospheres.

The Gaia-ESO Survey is obtaining high-quality spectroscopic data for about 10^5 stars using FLAMES at the VLT. UVES high-resolution spectra are being collected for about 5000 FGK-type stars. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO 2nd internal release and will be part of its 1st public release of advanced data products. The final parameter scale is tied to the one defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. A set of open and globular clusters is used to evaluate the physical soundness of the results. Each methodology is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted-medians of those from the individual methods. The recommended results successfully reproduce the benchmark stars atmospheric parameters and the expected Teff-log g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log g, and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50-100 K for Teff, 0.10-0.25 dex for log g, and 0.05-0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex.

We present a new grid of model photospheres for the SDSS-III/APOGEE survey of stellar populations of the Galaxy, calculated using the ATLAS9 and MARCS codes. New opacity distribution functions were generated to calculate ATLAS9 model photospheres. MARCS models were calculated based on opacity sampling techniques. The metallicity ([M/H]) spans from -5 to 1.5 for ATLAS and -2.5 to 0.5 for MARCS models. There are three main differences with respect to previous ATLAS9 model grids: a new corrected H2O linelist, a wide range of carbon ([C/M]) and alpha element [alpha/M] variations, and solar reference abundances from Asplund et al. 2005. The added range of varying carbon and alpha element abundances also extends the previously calculated MARCS model grids. Altogether 1980 chemical compositions were used for the ATLAS9 grid, and 175 for the MARCS grid. Over 808 thousand ATLAS9 models were computed spanning temperatures from 3500K to 30000K and log g from 0 to 5, where larger temperatures only have high gravities. The MARCS models span from 3500K to 5500K, and log g from 0 to 5. All model atmospheres are publically available online.

S-type stars are late-type giants whose atmosphere is enriched in carbon and s-process elements because of either extrinsic pollution by a binary companion or intrinsic nucleosynthesis and dredge-up on the thermally-pulsing AGB. A large grid of S-star model atmospheres has been computed covering the range 2700 < Teff < 4000 K with 0.5 < C/O < 0.99. ZrO and TiO band strength indices as well as VJHKL photometry are needed to disentangle Teff, C/O and [s/Fe]. A "best-model finding tool" was developed using a set of well-chosen indices and checked against photometry as well as low- and high-resolution spectroscopy. It is found that applying M-star model atmospheres (i.e., with a solar C/O ratio) to S stars can lead to errors on Teff up to 400K. We constrain the parameter space occupied by S stars of the vast sample of Henize stars in terms of Teff, [C/O] and [s/Fe].

The Hipparcos orbiting observatory has revealed a large number of helium-core-burning "clump" stars in the Galactic field. These low-mass stars exhibit signatures of extra-mixing processes that require modeling beyond the first dredge-up of standard models. The 12C/13C ratio is the most robust diagnostic of deep mixing, because it is insensitive to the adopted stellar parameters. In this work we present 12C/13C determinations in a sample of 34 Galactic clump stars as well as abundances of nitrogen, carbon and oxygen. Abundances of carbon were studied using the C2 Swan (0,1) band head at 5635.5 A. The wavelength interval 7980-8130 A with strong CN features was analysed in order to determine nitrogen abundances and 12C/13C isotope ratios. The oxygen abundances were determined from the [O I] line at 6300 A. Compared with the Sun and dwarf stars of the Galactic disk, mean abundances in the investigated clump stars suggest that carbon is depleted by about 0.2 dex, nitrogen is enhanced by 0.2 dex and oxygen is close to abundances in dwarfs. Comparisons to evolutionary models show that the stars fall into two groups: the one is of first ascent giants with carbon isotope ratios altered according to the first dredge-up prediction, and the other one is of helium-core-burning stars with carbon isotope ratios altered by extra mixing. The stars investigated fall to these groups in approximately equal numbers.

The first stars in the history of the Universe are likely to form in the dense central regions of 10^5-10^6 Msolar cold dark matter halos at z=10-50. The annihilation of dark matter particles in these environments may lead to the formation of so-called dark stars, which are predicted to be cooler, larger, more massive and potentially more long-lived than conventional population III stars. Here, we investigate the prospects of detecting high-redshift dark stars with the upcoming James Webb Space Telescope (JWST). We find that dark stars at z>6 are intrinsically too faint to be detected by JWST. However, by exploiting foreground galaxy clusters as gravitational telescopes, certain varieties of cool (Teff < 30000 K) dark stars should be within reach at redshifts up to z=10. If the lifetimes of dark stars are sufficiently long, many such objects may also congregate inside the first galaxies. We demonstrate that this could give rise to peculiar features in the integrated spectra of galaxies at high redshifts, provided that dark stars make up at least 1 percent of the total stellar mass in such objects.

It is debated whether the Milky Way bulge has the characteristics of a classical bulge sooner than those of a pseudobulge. Detailed abundance studies of bulge stars is a key to investigate the origin, history, and classification of the bulge. The aim is to add to the discussion on the origin of the bulge and to study detailed abundances determined from near-IR spectra for bulge giants already investigated with optical spectra, the latter also providing the stellar parameters which are very significant for the results of the present study. Especially, the important CNO elements are better determined in the near-IR. High-resolution, near-infrared spectra in the H band are recorded using the CRIRES spectrometer on the Very Large Telescope. The CNO abundances can all be determined from the numerous molecular lines in the wavelength range observed. Abundances of the alpha elements are also determined from the near-IR spectra. [O/Fe], [Si/Fe] and [S/Fe] are enhanced up to metallicities of at least [Fe/H]=-0.3, after which they decline. This suggests that the Milky Way bulge experienced a rapid and early star-formation history like that of a classical bulge. However, a similarity between the bulge trend and the trend of the local thick disk seems present. Such a similarity could suggest that the bulge has a pseudobulge origin. Our [C/Fe] trend does not show any increase with [Fe/H] which could have been expected if W-R stars have contributed substantially to the C abundances. No "cosmic scatter" can be traced around our observed abundance trends; the scatter found is expected, given the observational uncertainties.

We determine the metallicity distribution function (MDF) of the Galactic halo by means of a sample of 1638 metal-poor stars selected from the Hamburg/ESO objective-prism survey (HES). The sample was corrected for minor biases introduced by the strategy for spectroscopic follow-up observations of the metal-poor candidates, namely "best and brightest stars first". [...] We determined the selection function of the HES, which must be taken into account for a proper comparison between the HES MDF with MDFs of other stellar populations or those predicted by models of Galactic chemical evolution. The latter show a reasonable agreement with the overall shape of the HES MDF for [Fe/H] > -3.6, but only a model of Salvadori et al. (2007) with a critical metallicity for low-mass star formation of Z_cr = 10^{-3.4} * Z_Sun reproduces the sharp drop at [Fe/H] ~-3.6 present in the HES MDF. [...] A comparison of the MDF of Galactic globular clusters and of dSph satellites to the Galaxy shows qualitative agreement with the halo MDF, derived from the HES, once the selection function of the latter is included. However, statistical tests show that the differences between these are still highly significant. [ABSTRACT ABRIDGED]

The formation and evolution of the Milky Way bulge is not yet well understood and its classification is ambiguous. Constraints can, however, be obtained by studying the abundances of key elements in bulge stars. The aim of this study is to determine the chemical evolution of CNO, and a few other elements in stars in the Galactic bulge, and to discuss the sensitivities of the derived abundances from molecular lines. High-resolution, near-IR spectra in the H band were recorded using VLT/CRIRES. Due to the high and variable visual extinction in the line-of-sight towards the bulge, an analysis in the near-IR is preferred. The CNO abundances can all be determined simultaneously from the numerous molecular lines in the wavelength range observed. The three giant stars in Baade's window presented here are the first bulge stars observed with CRIRES. We have especially determined the CNO abundances, with uncertainties of less than 0.20 dex, from CO, CN, and OH lines. Since the systematic uncertainties in the derived CNO abundances due to uncertainties in the stellar fundamental parameters, notably Teff, are significant, a detailed discussion of the sensitivities of the derived abundances is included. We find good agreement between near-IR and optically determined O, Ti, Fe, and Si abundances. Two of our stars show a solar [C+N/Fe], suggesting that these giants have experienced the first dredge-up and that the oxygen abundance should reflect the original abundance of the giants. The two giants fit into the picture, in which there is no significant difference between the O abundance in bulge and thick-disk stars. Our determination of the S abundances is the first for bulge stars. The high [S/Fe] values for all the stars indicate a high star-formation rate in an early phase of the bulge evolution.

We aim to test the power of theoretical calibrations based on a new generation of MARCS models by comparisons with observational photomteric data. We calculate synthetic uvby-Hbeta colour indices from synthetic spectra. A sample of 388 field stars as well as stars in globular clusters is used for a direct comparison of the synthetic indices versus empirical data and for scrutinizing the possibilities of theoretical calibrations for temperature, metallicity and gravity. We show that the temperature sensitivity of the synthetic (b-y) colour is very close to its empirical counterpart, whereas the temperature scale based upon Hbeta shows a slight offset. The theoretical metallicity sensitivity of the m1 index (and for G-type stars its combination with c1) is somewhat larger than the empirical one, based upon spectroscopic determinations. The gravity sensitivity of the synthetic c1 index shows a satisfactory behaviour when compared to obervations of F stars. For stars cooler than the sun a deviation is significant in the c1-(b-y) diagram. The theoretical calibrations of (b-y), (v-y) and c1 seem to work well for Pop II stars and lead to effective temperatures for globular cluster stars supporting recent claims by Korn et al. (2007) that atomic diffusion occurs in stars near the turnoff point of NGC 6397. Synthetic colours of stellar atmospheres can indeed be used, in many cases, to derive reliable fundamental stellar parameters. The deviations seen when compared to observational data could be due to incomplete linelists but are possibly also due to effects of assuming plane-parallell or spherical geometry and LTE.

High resolution synthetic stellar libraries are of fundamental importance for the preparation of the Gaia Mission. We present new sets of spectral stellar libraries covering two spectral ranges: 300 --1100 nm at 0.1 nm resolution, and 840 -- 890 nm at 0.001 nm resolution. These libraries span a large range in atmospheric parameters, from super-metal-rich to very metal-poor (-5.0 $<$[Fe/H]$<$+1.0), from cool to hot (\teff=3000--50000 K) stars, including peculiar abundance variations. The spectral resolution, spectral type coverage and number of models represent a substantial improvement over previous libraries used in population synthesis models and in atmospheric analysis.

We have constructed a grid of about 10,000 spherically symmetric and plane-parallel models with the MARCS program, and make it available for public use. Parameter ranges are: Teff=2500 to 8000 K, log g =log(GM/R2)= -1 to 5 (cgs) with various masses and radii, [Me/H]=-5 to +1, with [Alpha/Fe] = 0.0 and 0.4 and different choices of C and N abundances to also represent stars of types R, S and N, and with microturbulence parameters from 1 to 5 km/s. We also list fluxes in approximately 108,000 wavelength points. Underlying assumptions in addition to 1D stratification include hydrostatic equilibrium, MLT convection and LTE. A number of general properties of the models are discussed, in relation to the effects of changing blanketing and sphericity. Models are compared with other available grids and excellent agreement is found with plane-parallel models of Castelli and Kurucz within the overlapping parameter range. Although there are departures from the spherically symmetric NextGen models, the agreement with more recent PHOENIX models is gratifying. The models of the grid show regularities, but some interesting departures from general patterns occur for the coolest models due to the molecular opacities. We have tested rules of thumb concerning effects of blanketing and sphericity and found them to often be astonishingly accurate. Some interesting new phenomena have been discovered, such as the intricate coupling between blanketing and sphericity, and the strong effects of carbon enhancement on metal-poor models. We give further details of models and comparisons with observations in subsequent papers.

The mass at which a transition is made between stars that have radiative or convective cores throughout the core H-burning phase is a fairly sensitive function of Z (particularly the CNO abundances). As a consequence, the ~4 Gyr, open cluster M67 provides a constraint on Z_\odot (and the solar heavy-element mixture) because (i) high-resolution spectroscopy indicates that this system has virtually the same metal abundances as the Sun, and (ii) its turnoff stars have masses just above the lower limit for sustained core convection on the main sequence. In this study, evolutionary tracks and isochrones using the latest MARCS model atmospheres as boundary conditions have been computed for 0.6-1.4 solar masses on the assumption of a metals mix (implying Z_\odot = 0.0125) based on the solar abundances derived by M. Asplund and collaborators using 3-D model atmospheres. These calculations do not predict a turnoff gap where one is observed in M67. No such difficulty is found if the analysis uses isochrones for Z_\odot = 0.0165, assuming the Grevesse & Sauval (1998) mix of heavy elements. Our findings, like the inferences from helioseismology, indicate a problem with the Asplund et al. abundances. However, it is possible that low-Z models with diffusive processes taken into account will be less problematic.

Stellar models have been computed for stars having [Fe/H] = 0.0 and -2.0 to determine the effects of using boundary conditions derived from the latest MARCS model atmospheres. The latter were fitted to the interior models at both the photosphere and at tau = 100, and at least for the 0.8-1.0 solar mass stars considered here, the resultant evolutionary tracks were found to be nearly independent of the chosen fitting point. Particular care was taken to treat the entire star as consistently as possible; i.e., both the interior and atmosphere codes assumed the same abundances and the same treatment of convection. Tracks were also computed using either the classical gray T(tau,T_eff) relation or that derived by Krishna Swamy (1966) to derive the boundary pressure. The latter predict warmer giant branches (by ~150 K) at solar abundances than those based on gray or MARCS atmospheres, which happens to be in good agreement with the inferred temperatures of giants in the open cluster M67 from the latest (V-K)-T_eff relations. Most of the calculations assumed Z=0.0125 (Asplund et al.), though a few models were computed for Z=0.0165 (Grevesse & Sauval) to determine the dependence of the tracks on Z_\odot. Grids of "scaled solar, differentially corrected" (SDC) atmospheres were also computed to try to improve upon theoretical MARCS models. When they were used as boundary conditions, the resultant tracks agreed very well with those based on a standard scaled-solar (e.g., Krishna Swamy) T(tau,T_eff) relation, independently of the assumed metal abundance. Fits of isochrones to the C-M diagram of the [Fe/H] = -2 globular cluster M68 were examined, as was the possibility that the mixing-length parameter varies with stellar parameters.

We present our on-going work on the determination of elemental abundances of giants in the Galactic Bulge by means of infrared spectroscopy. We show a preliminarily reduced spectrum and a synthetic spectrum fit of the Bulge giant Arp 4203 recorded with the near-infrared, high-resolution Crires spectrograph mounted on the VLT during its science verification run in August 2006. Abundances derived from this spectrum are discussed.

We have derived abundances of O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Fe, Ni and Ba for 43 metal-poor field stars in the solar neighbourhood with iron abundances [Fe/H] ranging from -0.4 to -3.0. Total abundance errors are estimated to about 0.1-0.2 dex for most elements, while relative abundances are most probably more accurate. We have made an extensive comparison with other similar studies and traced the reasons for the, in most cases, relatively small differences. Among the results we note the following: We find [O/Fe] from the O I IR- triplet lines to possibly increase beyond [Fe/H]=-1.0, though considerably less than obtained by others from abundances based on OH lines. We find correlations among the deviations for individual stars from the trends of individual alpha elements relative to Fe. This may support the idea that the formation of the Halo stars occurred in smaller systems with different star formation rates. We verify the finding by Gratton et al. (2003b) that stars that do not participate in the rotation of the galactic disk show a lower mean and larger spread in [alpha/Fe] than stars participating in the general rotation. The latter stars also seem to show some correlation between [alpha/Fe] and rotation speed. We finally advocate that a spectroscopic study of a larger sample of halo stars with well-defined selection criteria is very important, in addition to the very considerable efforts that various groups have already made. (Abridged)

Chemical abundances for 23 candidate subgiant stars have been derived with the aim at exploring their usefulness for studies of galactic chemical evolution. High-resolution spectra from ESO CAT-CES and NOT-SOFIN covered 16 different spectral regions in the visible part of the spectrum. Some 200 different atomic and molecular spectral lines have been used for abundance analysis of about 30 elemental species. The wings of strong, pressure-broadened metal lines were used for determination of stellar surface gravities, which have been compared with gravities derived from Hipparcos parallaxes and isochronic masses. Stellar space velocities have been derived from Hipparcos and Simbad data, and ages and masses were derived with recent isochrones. Only 12 of the stars turned out to be subgiants, i.e. on the "horizontal" part of the evolutionary track between the dwarf- and the giant stages. The abundances derived for the subgiants correspond closely to those of dwarf stars. With the possible exceptions of lithium and carbon we find that subgiant stars show no "chemical" traces of post-main-sequence evolution and that they are therefore very useful targets for studies of galactic chemical evolution.

The sensitivity of stellar spectra to alpha/Fe abundance changes is investigated with the aim to be detected photometricly and employed for scientific goals of the GAIA mission. A grid of plane parallel, line blanketed, flux constant, LTE model atmospheres with different [alpha/Fe] ratios was calculated. As a first step, the modelled stellar energy fluxes for solar-type stars and giants were computed and intercompared. The spectral sensitivity talpha/Fe abundance changes is noticeable and has to be taken into account when selecting photometric filters for GAIA. The Ca II H and K lines and Mg I b triplet are the most sensitive direct indicators of alpha/Fe abundance changes.

Europium abundances for 74 F and G dwarf stars of the galactic disk have been determined from the 4129.7 A Eu II line. The stars were selected from the sample of Edvardsson et al. (1993) and [Eu/Fe] shows a smaller scatter and a slightly weaker trend with [Fe/H] than found by Woolf et al. (1995). The data of the two analyses are homogenized and merged. We also discuss the adopted effective temperature scale.

High-resolution spectra of 13 core helium-burning stars in the thick disk of the Galaxy have been obtained with the SOFIN spectrograph on the Nordic Optical Telescope to investigate abundances of up to 22 chemical elements. Abundances of carbon were studied using the C_2 Swan (0,1) band head at 5635.5 A. The wavelength interval 7980-8130 A with strong CN features was analysed in order to determine nitrogen abundances and 12C/13C isotope ratios. The oxygen abundances were determined from the O I line at 6300 A. Abundances in the investigated stars suggest that carbon is depleted by about 0.3 dex, nitrogen is enhanced by more than 0.4 dex and oxygen is unaltered. The 12C/13C ratios are lowered and lie between values 3 and 7 which is in agreement with "cool bottom processing" predictions (Boothroyd & Sackmann 1999). The C/N ratios in the investigated stars are lowered to values between 0.7 and 1.2 which is less than present day theoretical predictions and call for further studies of stellar mixing processes. Abundance ratios of O, Mg, Eu and other heavy chemical elements to iron in the investigated stars show a pattern characteristic of thick disk stars. The results provide evidence that the thick disk population has a distinct chemical history from the thin disk. The onset of the bulk of SN Ia is suggested to appear around [Fe/H] = -0.6 dex.

We have observed the B I line at 2089.6 AA in the metal poor star HD 140283 with HST/STIS. The observation was an attempt to confirm the B abundances derived from the 2497 AA B I line in earlier works in general, and for for this star in particular (Edvardsson et al. 1994 A&A 290, 176; Kiselman & Carlsson 1996, A&A 311, 680). The resulting spectrum gained from 8 orbits of observations is hardly consistent with the boron abundance derived in the earlier works. A pure 10^B line could produce the feature observed but no reasonable process can produce such a ratio between 10^B and 11^B. More likely the analysed feature is affected by a statistical fluctuation. A conservative upper limit estimate of the NLTE corrected boron abundance gives a value marginally consistent with that derived from the 2497 AA line. We are proposing further observations to derive a more definite abundance from the line.

High-resolution spectra of six core helium-burning `clump' stars and three giants in the open cluster M 67 have been obtained with the SOFIN spectrograph on the Nordic Optical Telescope to investigate abundances of up to 25 chemical elements. Abundances of carbon were studied using the C2 Swan (0,1) band head at 5635.5 A. The wavelength interval 7980-8130 A with strong CN features was analysed in order to determine nitrogen abundances and 12C/13C isotope ratios. The oxygen abundances were determined from the [O I] line at 6300 A. The overall metallicity of the cluster stars was found to be close to solar ([Fe/H]=-0.03+-0.03). Compared with the Sun and other dwarf stars of the Galactic disk, as well as with dwarf stars of M 67 itself, abundances in the investigated stars suggest that carbon is depleted by about 0.2 dex, nitrogen is enhanced by about 0.2 dex and oxygen is unaltered. Among other mixing- sensitive chemical elements an overabundance of sodium may be suspected. The mean C/N and 12C/13C ratios are lowered to the values of 1.7+-0.2 and 24+-4 in the giants and to the values of 1.4+-0.2 and 16+-4 in the clump stars. These results suggest that extra mixing of CN-cycled material to the stellar surface takes place after the He-core flash. Abundances of heavy chemical elements in all nine stars were found to be almost identical and close to solar.

Abundance analysis of carbon has been performed in a sample of 80 late F and early G type dwarf stars in the metallicity range -1.06 < [Fe/H] < 0.26 using the forbidden [C I] line at 8727 A. This line is presumably less sensitive to temperature, atmospheric structure and departures from LTE than alternative carbon criteria. We find that [C/Fe] decreases slowly with increasing [Fe/H] with an overall slope of -0.17 +- 0.03. Our results are consistent with carbon enrichment by superwinds of metal-rich massive stars but inconsistent with a main origin of carbon in low-mass stars. This follows in particular from a comparison between the relation of [C/O] with metallicity for the Galactic stars and the corresponding relation observed for dwarf irregular galaxies. The significance of intermediate-mass stars for the production of carbon in the Galaxy is still somewhat unclear.

We have observed the B I 2497 A line to derive the boron abundances of two very metal-poor stars selected to help in tracing the origin and evolution of this element in the early Galaxy: BD +23 3130 and HD 84937. The observations were conducted using the Goddard High Resolution Spectrograph on board the Hubble Space Telescope. A very detailed abundance analysis via spectral synthesis has been carried out for these two stars, as well as for two other metal-poor objects with published spectra, using both Kurucz and OSMARCS model photospheres, and taking into account consistently the NLTE effects on the line formation. We have also re-assessed all published boron abundances of old disk and halo unevolved stars. Our analysis shows that the combination of high effective temperature (Teff > 6000 K, for which boron is mainly ionized) and low metallicity ([Fe/H]<-1) makes it difficult to obtain accurate estimates of boron abundances from the B I 2497 A line. This is the case of HD 84937 and three other published objects (including two stars with [Fe/H] ~ -3), for which only upper limits can be established. BD +23 3130, with [Fe/H] ~ -2.9 and logN(B)_NLTE=0.05+/-0.30, appears then as the most metal-poor star for which a firm measurement of the boron abundance presently exists. The evolution of the boron abundance with metallicity that emerges from the seven remaining stars with Teff < 6000 K and [Fe/H]<-1, for which beryllium abundances were derived using the same stellar parameters, shows a linear increase with a slope ~ 1. Furthermore, the B/Be ratio found is constant at a value ~ 20 for stars in the range -3<[Fe/H]<-1. These results point to spallation reactions of ambient protons and alpha particles with energetic particles enriched in CNO as the origin of boron and beryllium in halo stars.

We have made new abundance determinations for a sample of NaMgAl stars. These stars, which are a subgroup of the nearby metal-rich field F and G disk dwarfs, were first identified by Edvardsson et al. (1993) on the basis of their apparent enrichment in Na, Mg and Al relative to other elements. The discovery of a planetary companion to the nearby solar type star 51 Peg (Mayor & Queloz 1995) combined with Edvardsson et al.'s earlier identification of 51 Peg as a NaMgAl star highlighted the group's potential importance. Our new analysis, which uses new spectra of higher resolution and better wavelength coverage than the analysis of Edvardsson et al., shows that the Na, Mg and Al abundances of the NaMgAl stars are indistinguishable from those of non-NaMgAl stars with otherwise similar properties. The group thus appears to be spurious. Our study, which includes 51 Peg, also provides the most complete set of abundances for this star available to date. The new Fe abundance, $[Fe/H] = +0.20 \pm 0.07$, of 51 Peg confirms earlier measurements of its metal richness. Abundances for 19 other elements, including C, N and O, reveal a fairly uniform enrichment similar to that of Fe and show no evidence of abnormality compared to other metal rich stars of similar spectral type.

We have determined carbon abundances or upper limits from the forbidden [C I] line at 8727.13 A for 85 F and early G type main sequence stars, with metallicities ranging between [Fe/H] = -1.0 and [Fe/H] = +0.25. The [C I] line has not been previously used for a study of this size. We find that the C/Fe abundance ratio is slowly decreasing with time and increasing metallicity in the disk, but with considerable observational scatter. A discussion of results in the literature supports this result. Our data do not reveal any gradients with galactocentric birth distance for our stars, but more accurate observations are desirable since possible gradients should be important to the question of the synthesis of C and Fe in the disk. (Accepted for publication in A&A. Postscript text and figures available via anonymous ftp to "chiron.astro.uu.se", in directory "pub/articles/atmos/P89")

We have determined chemical abundances and radial velocities for stars in the field of the zeta Sculptoris cluster. We find that the cluster metal deficiency previously found from UBV photometry is too high; the cluster overall metallicity, [Fe/H]=+0.24, is about 70% higher than the solar value. The chemical abundance pattern is unusual: the Ni/Fe ratio is significantly lower than typical for field stars with similar Fe and Ca abundances, and the cluster may also be deficient in the alpha element Si. For its age, the cluster is unusually far from the galactic disk, approximately 240 pc. The adopted heliocentric cluster radial velocity, +3.9 km/s, shows that it is close to its maximum distance from the galactic plane. We suggest that the zeta Sculptoris cluster was formed in the galactic disk 45 Myr ago by the interaction of a high velocity cloud with the interstellar medium, and that its formation may be connected with that of Gould's Belt. (Table3, postscript figures and text available by anonymous ftp from "chiron.astro.uu.se" in directory pub/articles/atmos/P90)

Using observations of the 2496.7 A B I line with the HST GHRS at a nominal resolution of 90,000, we have found the abundance of boron of HD 140283 to be log epsilon_B (= 12 + log (N_B/N_H)) = 0.34 \pm 0.20. This result is found when a significant non-LTE effect in the formation of the B I line is taken into account. The resulting N_B/N_Be ratio is about 17 (in the range 9 -- 34), which is in very good agreement with what is expected from spallation by cosmic rays. We conclude that this origin of Be and B in the Early Galaxy is the most probable of recently suggested formation mechanisms. (Submitted to A&A. Postscript figures available via annymous ftp to "chiron.astro.uu.se" in directory pub/Preprints/P87)