We have been monitoring the metal-poor Population II giant, HD 122563, for radial velocity variations since 2016 using the SONG telescope on Tenerife. We have detected the global seismic quantity, numax, which provides information related to the stellar parameters. By combining these data with complementary data, we derive a new precise surface gravity, radius and distance to the star. Our results are corroborated by using the parallax from Gaia DR2. We present these results and some of their implications.

The last decade has seen a rapid development in asteroseismology thanks to the CoRoT and Kepler missions. With more detailed asteroseismic observations available, it is becoming possible to infer exactly how oscillations are driven and dissipated in solar-type stars. We have carried out three-dimensional (3D) stellar atmosphere simulations together with one-dimensional (1D) stellar structural models of key benchmark turn-off and subgiant stars to study this problem from a theoretical perspective. Mode excitation and damping rates are extracted from 3D and 1D stellar models based on analytical expressions. Mode velocity amplitudes are determined by the balance between stochastic excitation and linear damping, which then allows the estimation of the frequency of maximum oscillation power, $\nu_{\max}$, for the first time based on ab initio and parameter-free modelling. We have made detailed comparisons between our numerical results and observational data and achieved very encouraging agreement for all of our target stars. This opens the exciting prospect of using such realistic 3D hydrodynamical stellar models to predict solar-like oscillations across the HR-diagram, thereby enabling accurate estimates of stellar properties such as mass, radius and age.

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\epsilon_{\mathrm{N}} = 7.77$, with the systematic $1\sigma$ uncertainty estimated to be $0.05$ dex. This result is consistent with earlier studies after correcting for differences in line selections and equivalent widths.

Context. The pursuit of more realistic spectroscopic modelling and consistent abundances has led us to begin a new series of papers designed to improve current solar and stellar abundances of various atomic species. To achieve this, we have began updating the three-dimensional (3D) non-local thermodynamic equilibrium (non-LTE) radiative transfer code, Multi3D, and the equivalent one-dimensional (1D) non-LTE radiative transfer code, MULTI. Aims. We examine our improvements to these codes by redetermining the solar barium abundance. Barium was chosen for this test as it is an important diagnostic element of the s-process in the context of galactic chemical evolution. New Ba II + H collisional data for excitation and charge exchange reactions computed from first principles had recently become available and were included in the model atom. The atom also includes the effects of isotopic line shifts and hyperfine splitting. Method. A grid of 1D LTE barium lines were constructed with MULTI and fit to the four Ba II lines available to us in the optical region of the solar spectrum. Abundance corrections were then determined in 1D non-LTE, 3D LTE, and 3D non-LTE. A new 3D non-LTE solar barium abundance was computed from these corrections. Results. We present for the first time the full 3D non-LTE barium abundance of $A({\rm Ba})=2.27\pm0.02\pm0.01$, which was derived from four individual fully consistent barium lines. Errors here represent the systematic and random errors, respectively.

Manganese (Mn) is a key Fe-group elements, commonly employed in stellar population and nucleosynthesis studies to explore the role of SN Ia. We have developed a new non-local thermodynamic equilibrium (NLTE) model of Mn, including new photo-ionisation cross-sections and new transition rates caused by collisions with H and H- atoms. We applied the model in combination with 1-dimensional (1D) LTE model atmospheres and 3D hydrodynamical simulations of stellar convection to quantify the impact of NLTE and convection on the line formation. We show that the effects of NLTE are present in Mn I and, to a lesser degree, in Mn II lines, and these increase with metallicity and with effective temperature of a model. Employing 3D NLTE radiative transfer, we derive new abundance of Mn in the Sun, A(Mn)=5.52 +/- 0.03 dex, consistent with the element abundance in C I meteorites. We also apply our methods to the analysis of three metal-poor benchmark stars. We find that 3D NLTE abundances are significantly higher than 1D LTE. For dwarfs, the differences between 1D NLTE and 3D NLTE abundances are typically within 0.15 dex, however, the effects are much larger in the atmospheres of giants owing to their more vigorous convection. We show that 3D NLTE successfully solves the ionisation and excitation balance for the RGB star HD 122563 that cannot be achieved by 1D LTE or 1D NLTE modelling. For HD 84937 and HD 140283, the ionisation balance is satisfied, however, the resonance Mn I triplet lines still show somewhat lower abundances compared to the high-excitation lines. Our results for the benchmark stars confirm that 1D LTE modelling leads to significant systematic biases in Mn abundances across the full wavelength range from the blue to the IR. We also produce a list of Mn lines that are not significantly biased by 3D and can be reliably, within the 0.1 dex uncertainty, modelled in 1D NLTE.

The amplitude of solar p-mode oscillations is governed by stochastic excitation and mode damping, both of which take place in the surface convection zone. However, the time-dependent, turbulent nature of convection makes it difficult to self-consistently study excitation and damping processes through the use of traditional one-dimensional hydrostatic models. To this end, we carried out \textit{ab initio} three-dimensional, hydrodynamical numerical simulations of the solar atmosphere to investigate how p-modes are driven and dissipated in the Sun. The description of surface convection in the simulations is free from the tuneable parameters typically adopted in traditional one-dimensional models. Mode excitation and damping rates are computed based on analytical expressions whose ingredients are evaluated directly from the three-dimensional model. With excitation and damping rates both available, we estimate the theoretical oscillation amplitude and frequency of maximum power, $\nu_{\max}$, for the Sun. We compare our numerical results with helioseismic observations, finding encouraging agreement between the two. The numerical method presented here provides a novel way to investigate the physical processes responsible for mode driving and damping, and should be valid for all solar-type oscillating stars.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many 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 $\epsilon_{\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.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

The nearby metal-poor giant HD122563 is an important astrophysical laboratory for which to test stellar atmospheric and interior physics. It is also a benchmark star for which to calibrate methods to apply to large scale surveys. Recently it has been remeasured using various methodologies given the new high precision instruments at our disposal. However, inconsistencies in the observations and models have been found. In order to better characterise this star we have been measuring its radial velocities since 2016 using the Hertzsprung telescope (SONG network node). In this work we report the first detections of sun-like oscillations in this star, and to our knowledge, a detection in the most metal-poor giant to date. We apply the classical seismic scaling relation to derive a new surface gravity of $\log g_{\nu} = 1.39 \pm 0.01$ dex. Constraints on the mass imposed by its PopII giant classification then yield a radius of $30.8 \pm 1.0$ R$_{\odot}$. By coupling this with recent interferometric measurements we infer a distance to the star of 306 $\pm$ 9 pc. Data from the Gaia mission corroborates the distance hypothesis ($d_{\rm GDR2}$ = 290 $\pm$ 5 pc), and thus the updated fundamental parameters. We confirm the validity of the seismic scaling relation without corrections for surface gravity in metal-poor and evolved star regimes. The small discrepancy of 0.04 dex reduces to 0.02 dex by applying corrections to the scaling relations. The new constraints on the HR diagram ($L_{\odot} = 381 \pm 26$) reduce the disagreement between the stellar parameters and evolution models, however, a discrepancy still exists. Fine-tuned stellar evolution calculations show that this can be reconciled by changing the mixing-length parameter by an amount (--0.35) that is in agreement with predictions from recent 3D simulations and empirical results.

Detailed chemical abundances of Galactic stars are needed in order to improve our knowledge of the formation and evolution of our galaxy, the Milky Way. We took advantage of the GIANO archive spectra to select a sample of Galactic disc stars in order to derive their chemical inventory and to compare the abundances we derived from these infrared spectra to the chemical pattern derived from optical spectra. We analysed high-quality spectra of 40 stars observed with GIANO. We derived the stellar parameters from the photometry and the Gaia data-release 2 (DR2) parallax; the chemical abundances were derived with the code MyGIsFOS. For a subsample of stars we compared the chemical pattern derived from the GIANO spectra with the abundances derived from optical spectra. We derived P abundances for all 40 stars, increasing the number of Galactic stars for which phosphorus abundance is known. We could derive abundances of 14 elements, 8 of which are also derived from optical spectra. The comparison of the abundances derived from infrared and optical spectra is very good. The chemical pattern of these stars is the one expected for Galactic disc stars and is in agreement with the results from the literature. GIANO is providing the astronomical community with an extremely useful instrument, able to produce spectra with high resolution and a wide wavelength range in the infrared.

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\epsilon_{\mathrm{O}}=8.69\pm0.03$ from the triplet alone, strengthening the case for a low solar oxygen abundance.

GALAH is a large-scale magnitude-limited southern stellar spectroscopic survey. Its second data release (GALAH DR2) provides values of stellar parameters and abundances of 23 elements for 342,682 stars (Buder et al.). Here we add a description of the public release of radial velocities with a typical accuracy of 0.1 km/s for 336,215 of these stars, achievable due to the large wavelength coverage, high resolving power and good signal to noise ratio of the observed spectra, but also because convective motions in stellar atmosphere and gravitational redshift from the star to the observer are taken into account. In the process we derive medians of observed spectra which are nearly noiseless, as they are obtained from between 100 and 1116 observed spectra belonging to the same bin with a width of 50 K in temperature, 0.2 dex in gravity, and 0.1 dex in metallicity. Publicly released 1181 median spectra have a resolving power of 28,000 and trace the well-populated stellar types with metallicities between -0.6 and +0.3. Note that radial velocities from GALAH are an excellent match to the accuracy of velocity components along the sky plane derived by Gaia for the same stars. The level of accuracy achieved here is adequate for studies of dynamics within stellar clusters, associations and streams in the Galaxy. So it may be relevant for studies of the distribution of dark matter.

Female role models reduce the impact on women of stereotype threat, i.e., of being at risk of conforming to a negative stereotype about one's social, gender, or racial group. This can lead women scientists to underperform or to leave their scientific career because of negative stereotypes such as, not being as talented or as interested in science as men. Sadly, history rarely provides role models for women scientists; instead, it often renders these women invisible. In response to this situation, we present a selection of twelve outstanding women who helped to develop nuclear astrophysics.

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$\alpha$, H$\beta$, and H$\gamma$, 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$\gamma$, while the inner wings can be weaker in 3D models, particularly for H$\alpha$. For H$\alpha$, 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$\alpha$, H$\beta$, and H$\gamma$; 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\sigma$ uncertainties. For H$\alpha$, 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$\alpha$ 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.

Context. The surface structures and dynamics of cool stars are characterized by the presence of convective motions and turbulent flows which shape the emergent spectrum. Aims. We used realistic three-dimensional radiative hydrodynamical simulations from the Stagger-grid to calculate synthetic spectra with the radiative transfer code Optim3D for stars with different stellar parameters to predict photometric colors and convective velocity shifts. Methods. We calculated spectra from 1000 to 200 000 A with a constant resolving power of 20 000 and from 8470 and 8710 A (Gaia Radial Velocity Spectrometer - RVS - spectral range), with a constant resolving power of 300 000. Results. We used synthetic spectra to compute theoretical colors in the Johnson-Cousins UBV(RI), SDSS, 2MASS, Gaia, SkyMapper, Stromgren systems, and HST-WFC3. We showed that 1D versus 3D differences are limited to a small percent except for the narrow filters that span the optical and UV region of the spectrum. In addition, we derived the effect of the convective velocity fields on selected Fe I lines. We found the overall convective shift for 3D simulations with respect to the reference 1D hydrostatic models, revealing line shifts of between -0.235 and +0.361 km/s. We showed a net correlation of the convective shifts with the effective temperature: lower effective temperatures denote redshifts and higher effective temperatures denote blueshifts. We conclude that the extraction of accurate radial velocities from RVS spectra need an appropriate wavelength correction from convection shifts. Conclusions. The use of realistic 3D hydrodynamical stellar atmosphere simulations has a small but significant impact on the predicted photometry compared with classical 1D hydrostatic models for late-type stars. We make all the spectra publicly available for the community through the POLLUX database.

We present an abundance analysis of the low-metallicity benchmark red giant star HD 122563 based on realistic, state-of-the-art, high-resolution, three-dimensional (3D) model stellar atmospheres including non-grey radiative transfer through opacity binning with four, twelve, and 48 bins. The 48-bin 3D simulation reaches temperatures lower by ~ 300 - 500 K than the corresponding 1D model in the upper atmosphere. Small variations in the opacity binning, adopted line opacities, or chemical mixture can cool the photospheric layers by a further ~ 100 - 300 K and alter the effective temperature by ~ 100 K. A 3D local thermodynamic equilibrium (LTE) spectroscopic analysis of Fe I and Fe II lines gives discrepant results in terms of derived Fe abundance, which we ascribe to non-LTE effects and systematic errors on the stellar parameters. We also determine C, N, and O abundances by simultaneously fitting CH, OH, NH, and CN molecular bands and lines in the ultraviolet, visible, and infrared. We find a small positive 3D-1D abundance correction for carbon (+0.03 dex) and negative ones for nitrogen (-0.07 dex) and oxygen (-0.34 dex). From the analysis of the [O I] line at 6300.3 {\AA}, we derive a significantly higher oxygen abundance than from molecular lines (+0.46 dex in 3D and +0.15 dex in 1D). We rule out important OH photodissociation effects as possible explanation for the discrepancy and note that lowering the surface gravity would reduce the oxygen abundance difference between molecular and atomic indicators.

HARPS-N spectra with S/N > 250 and MARCS model atmospheres were used to derive abundances of C, O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Zn, and Y in ten stars from the Kepler LEGACY sample (including the binary pair 16 Cyg A and B) selected to have metallicities in the range -0.15 < [Fe/H] < +0.15 and ages between 1 and 7 Gyr. Stellar gravities were obtained from seismic data and effective temperatures were determined by comparing non-LTE iron abundances derived from FeI and FeII lines. Available non-LTE corrections were also applied when deriving abundances of the other elements. The results support the [X/Fe]-age relations previously found for solar twins. [Mg/Fe], [Al/Fe], and [Zn/Fe] decrease by ~0.1 dex over the lifetime of the Galactic thin disk due to delayed contribution of iron from Type Ia supernovae relative to prompt production of Mg, Al, and Zn in Type II supernovae. [Y/Mg] and [Y/Al], on the other hand, increase by ~0.3 dex, which can be explained by an increasing contribution of s-process elements from low-mass AGB stars as time goes on. The trends of [C/Fe] and [O/Fe] are more complicated due to variations of the ratio between refractory and volatile elements among stars of similar age. Two stars with about the same age as the Sun show very different trends of [X/H] as a function of elemental condensation temperature Tc and for 16 Cyg, the two components have an abundance difference, which increases with Tc. These anomalies may be connected to planet-star interactions.

From exploratory studies and theoretical expectations it is known that simplifying approximations in spectroscopic analysis (LTE, 1D) lead to systematic biases of stellar parameters and abundances. These biases depend strongly on surface gravity, temperature, and, in particular, for LTE vs. non-LTE (NLTE) on metallicity of the stars. Here we analyse the [Mg/Fe] and [Fe/H] plane of a sample of 326 stars, comparing LTE and NLTE results obtained using 1D hydrostatic models and averaged <3D> models. We show that compared to the <3D>NLTE benchmark, all other three methods display increasing biases towards lower metallicities, resulting in false trends of [Mg/Fe] against [Fe/H], which have profound implications for interpretations by chemical evolution models. In our best <3D> NLTE model, the halo and disc stars show a clearer behaviour in the [Mg/Fe] - [Fe/H] plane, from the knee in abundance space down to the lowest metallicities. Our sample has a large fraction of thick disc stars and this population extends down to at least [Fe/H] ~ -1.6 dex, further than previously proven. The thick disc stars display a constant [Mg/Fe] ~ 0.3 dex, with a small intrinsic dispersion in [Mg/Fe] that suggests that a fast SN Ia channel is not relevant for the disc formation. The halo stars reach higher [Mg/Fe] ratios and display a net trend of [Mg/Fe] at low metallicities, paired with a large dispersion in [Mg/Fe]. These indicate the diverse origin of halo stars from accreted low-mass systems to stochastic/inhomogeneous chemical evolution in the Galactic halo.

We determine Mg abundances in 6 Gaia benchmark stars using theoretical one-dimensional (1D) hydrostatic model atmospheres, as well as temporally- and spatially-averaged 3D model atmospheres (<3D>). The stars cover a range of Teff from 4700 to 6500 K, log g from 1.6 to 4.4 dex, and [Fe/H] from -3.0 dex to solar. Spectrum synthesis calculations are performed in local thermodynamic equilibrium (LTE) and in non-LTE (NLTE) using the oscillator strengths recently published by Pehlivan Rhodin et al. We find that: a) Mg abundances determined from the infrared spectra are as accurate as the optical diagnostics, b) the NLTE effects on Mg I line strengths and abundances in this sample of stars are minor (although for a few Mg I lines the NLTE effects on abundance exceed 0.6 dex in <3D> and 0.1 dex in 1D, c) the solar Mg abundance is 7.56 +/- 0.05 dex (total error), in the excellent agreement with the Mg abundance measured in CI chondritic meteorites, d) the 1D NLTE and <3D> NLTE approach can be used with confidence to analyse optical Mg I lines in spectra of dwarfs and sub-giants, but for red giants the Mg I 5711 A line should be preferred, e) low-excitation Mg I lines are sensitive to the atmospheric structure; for these lines, LTE calculations with <3D> models lead to significant systematic abundance errors. The methods developed in this work will be used to study Mg abundances of a large sample of stars in the next paper in the series.

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.

Context. Red giant branch (RGB) stars are very bright objects in galaxies and are often used as standard candles. Interferometry is the ideal tool to characterize the dynamics and morphology of their atmospheres. Aims. We aim at precisely characterising the surface dynamics of a sample of RGB stars. Methods. We obtained interferometric observations for three RGB stars with the MIRC instrument mounted at the CHARA interfer- ometer. We looked for asymmetries on the stellar surfaces using limb-darkening models. Results. We measured the apparent diameters of HD197989 (Epsilon Cyg) = 4.61+-0.02 mas, HD189276 (HR7633) = 2.95+-0.01 mas, and HD161096 (Beta Oph) = 4.43+-0.01 mas. We detected departures from the centrosymmetric case for all three stars with the tendency of a greater effect for lower logg of the sample. We explored the causes of this signal and conclude that a possible explanation to the interferometric signal is the convection-related and/or the magnetic-related surface activity. However, it is necessary to monitor these stars with new observations, possibly coupled with spectroscopy, in order to firmly establish the cause.

During primary transits, the spectral signatures of an exoplanet atmosphere can be measured using transmission spectroscopy. The goal of this work is to accurately measure the atomspheric sodium absorption light curve in HD189733b, correcting for the effects of stellar differential limb-darkening, stellar activity and a "bump" caused by the changing radial velocity of the exoplanet. In fact, due to the high cadence and quality of our data, it is the first time that the last feature can be detected even by visual inspection. We use 244 high-resolution optical spectra taken by the UVES instrument mounted at the VLT. Our observations cover a full transit of HD 189733b, with a cadence of 45 seconds. To probe the transmission spectrum of sodium we produce excess light- curves integrating the stellar flux in passbands of 1 \AA, 1.5 \AA, and 3 \AA inside the core of each sodium D-line. We model the effects of external sources on the excess light-curves, which correspond to an observed stellar flare beginning close to mid-transit time and the wavelength dependent limb-darkening effects. In addition, by characterizing the effect of the changing radial velocity and Doppler shifts of the planetary sodium lines inside the stellar sodium lines, we estimate the depth and width of the exoplanetary sodium feature. We estimate the shape of the planetary sodium line by a Gaussian profile, with an equivalent width of 0.0023 \pm 0.0010 \AA, thereby confirming the presence of sodium in the atmosphere of HD189733b with excess absorption levels of 0.72 \pm 0.25%, 0.34 \pm 0.11%, and 0.20 \pm 0.06% for the integration bands of 1 \AA, 1.5 \AA, and 3 \AA, respectively. From these, we produce a first order estimate of the number density of sodium in the exoplanet atmosphere.

We present the scientific case for a future space-based telescope aimed at very high spatial and temporal resolution imaging of the solar photosphere and chromosphere. Previous missions (e.g., HINODE, SUNRISE) have demonstrated the power of observing the solar photosphere and chromosphere at high spatial resolution without contamination from Earth's atmosphere. We argue here that increased spatial resolution (from currently 70 km to 25 km in the future) and high temporal cadence of the observations will vastly improve our understanding of the physical processes controlling solar magnetism and its characteristic scales. This is particularly important as the Sun's magnetic field drives solar activity and can significantly influence the Sun-Earth system. At the same time a better knowledge of solar magnetism can greatly improve our understanding of other astrophysical objects.

We analyse the effect on adiabatic stellar oscillation frequencies of replacing the near-surface layers in 1D stellar structure models with averaged 3D stellar surface convection simulations. The main difference is an expansion of the atmosphere by 3D convection, expected to explain a major part of the asteroseismic surface effect; a systematic overestimation of p-mode frequencies due to inadequate surface physics. We employ pairs of 1D stellar envelope models and 3D simulations from a previous calibration of the mixing-length parameter, alpha. That calibration constitutes the hitherto most consistent matching of 1D models to 3D simulations, ensuring that their differences are not spurious, but entirely due to the 3D nature of convection. The resulting frequency shift is identified as the structural part of the surface effect. The important, typically non-adiabatic, modal components of the surface effect are not included in the present analysis, but relegated to future papers. Evaluating the structural surface effect at the frequency of maximum mode amplitude, $\nu_{\rm max}$, we find shifts from $\delta\nu$=-0.8 microHz for giants at $\log g$=2.2 to -35 microHz for a ($T_{\rm eff}=6901$ K, $\log g$=4.29) dwarf. The fractional effect $\delta\nu(\nu_{\rm max})/\nu_{\rm max}$, ranges from -0.1% for a cool dwarf (4185 K, 4.74) to -6% for a warm giant (4962 K, 2.20).

Stellar activity and convection-related surface structures might cause bias in planet detection and characterization that use these transits. Surface convection simulations help to quantify the granulation signal. We used realistic three-dimensional radiative hydrodynamical simulations from the Stagger grid and synthetic images computed with the radiative transfer code Optim3D to model the transits of three prototype planets: a hot Jupiter, a hot Neptune, and a terrestrial planet. We computed intensity maps from RHD simulations of the Sun and a K-dwarf star at different wavelength bands from optical to far-infrared. We modeled the transit using synthetic stellar-disk images and emulated the temporal variation of the granulation intensity. We identified two types of granulation noise that act simultaneously during the planet transit: (i) the intrinsic change in the granulation pattern with timescales smaller than the usual planet transit, and (ii) the fact that the transiting planet occults isolated regions of the photosphere that differ in local surface brightness. Our modeling approach shows that the granulation pattern has a non-negligible effect on the light curve depth during the transit, and, consequentially on the determination of the radius of the planet transiting. The granulation noise appears to be correlated among the different wavelength ranges either in the visible or in the infrared regions. The granulation has to be considered as an intrinsic uncertainty (as a result of stellar variability) on the precise measurements of exoplanet transits of planets. The full characterization of the granulation is essential for determining the degree of uncertainty on the planet parameters. In this context, the use of 3D RHD simulations is important to measure the convection-related fluctuations.

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).

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.

The globular cluster 47 Tuc exhibits a complex sub-giant branch (SGB) with a faint-SGB comprising only about the 10% of the cluster mass and a bright-SGB hosting at least two distinct populations.We present a spectroscopic analysis of 62 SGB stars including 21 faint-SGB stars. We thus provide the first chemical analysis of the intriguing faint-SGB population and compare its abundances with those of the dominant populations. We have inferred abundances of Fe, representative light elements C, N, Na, and Al, {\alpha} elements Mg and Si for individual stars. Oxygen has been obtained by co-adding spectra of stars on different sequences. In addition, we have analysed 12 stars along the two main RGBs of 47 Tuc. Our principal results are: (i) star-to-star variations in C/N/Na among RGB and bright-SGB stars; (ii) substantial N and Na enhancements for the minor population corresponding to the faint-SGB; (iii) no high enrichment in C+N+O for faint-SGB stars. Specifically, the C+N+O of the faint-SGB is a factor of 1.1 higher than the bright-SGB, which, considering random (+-1.3) plus systematic errors (+-0.3), means that their C+N+O is consistent within observational uncertainties. However, a small C+N+O enrichment for the faint-SGB, similar to what predicted on theoretical ground, cannot be excluded. The N and Na enrichment of the faint-SGB qualitatively agrees with this population possibly being He-enhanced, as suggested by theory. The iron abundance of the bright and faint-SGB is the same to a level of ~0.10 dex, and no other significant difference for the analysed elements has been detected.

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.

The phenomenon of multiple populations in globular clusters is still far from understood, with several proposed mechanisms to explain the observed behaviour. The study of elemental and isotopic abundance patterns are crucial for investigating the differences among candidate pollution mechanisms. We derive magnesium isotopic ratios for 13 stars in the globular cluster 47 Tucanae (NGC 104) to provide new, detailed information about the nucleosynthesis that has occurred within the cluster. For the first time, the impact of 3D model stellar atmospheres on the derived Mg isotopic ratios is investigated. Using both tailored 1D atmospheric models and 3D hydrodynamical models, we derive magnesium isotopic ratios from four features of MgH near 5135{\AA} in 13 giants near the tip of the RGB, using high signal-to-noise, high-resolution spectra. We derive the magnesium isotopic ratios for all stars and find no significant offset of the isotopic distribution between the pristine and the polluted populations. Furthermore, we do not detect any statistically significant differences in the spread in the Mg isotopes in either population. No trends were found between the Mg isotopes and [Al/Fe]. The inclusion of 3D atmospheres has a significant impact on the derived 25Mg/24Mg ratio, increasing it by a factor of up to 2.5, compared to 1D. The 26Mg/24Mg ratio, on the other hand, essentially remains unchanged. We confirm the results seen from other globular clusters, where no strong variation in the isotopic ratios is observed between stellar populations, for observed ranges in [Al/Fe]. We see no evidence for any significant activation of the Mg-Al burning chain. The use of 3D atmospheres causes an increase of a factor of up to 2.5 in the fraction of 25Mg, resolving part of the discrepancy between the observed isotopic fraction and the predictions from pollution models.

The OI 777 nm lines are among the most commonly used diagnostics for the oxygen abundances in the atmospheres of FGK-type stars. However, they form in conditions that are far from local thermodynamic equilibrium (LTE). We explore the departures from LTE of atomic oxygen, and their impact on OI lines, across the Stagger-grid of three-dimensional hydrodynamic model atmospheres. For the OI 777 nm triplet we find significant departures from LTE. These departures are larger in stars with larger effective temperatures, smaller surface gravities, and larger oxygen abundances. We present grids of predicted 3D non-LTE based equivalent widths for the OI616nm, [OI] 630 nm, [OI] 636 nm, and OI 777 nm lines, as well as abundance corrections to 1D LTE based results.

We revisit the Galactic chemical evolution of oxygen, addressing the systematic errors inherent in classical determinations of the oxygen abundance that arise from the use of one dimensional hydrostatic (1D) model atmospheres and from the assumption of local thermodynamic equilibrium (LTE). We perform detailed 3D non-LTE radiative transfer calculations for atomic oxygen lines across a grid of 3D hydrodynamic stag- ger model atmospheres for dwarfs and subgiants. We apply our grid of predicted line strengths of the [OI] 630 nm and OI 777 nm lines using accurate stellar parameters from the literature. We infer a steep decay in [O/Fe] for [Fe/H] $\gtrsim$ -1.0, a plateau [O/Fe] $\approx$ 0.5 down to [Fe/H] $\approx$ -2.5 and an increasing trend for [Fe/H] $\lesssim$ -2.5. Our 3D non-LTE calculations yield overall concordant results from the two oxygen abundance diagnostics.

We present a high resolution spectroscopic analysis of 62 red giants in the Milky Way globular cluster NGC5286. We have determined abundances of representative light proton-capture, alpha, Fe-peak and neutron-capture element groups, and combined them with photometry of multiple sequences observed along the colour-magnitude diagram. Our principal results are: (i) a broad, bimodal distribution in s-process element abundance ratios, with two main groups, the s-poor and s-rich groups; (ii) substantial star-to-star Fe variations, with the s-rich stars having higher Fe, e.g. <[Fe/H]>_s-rich - <[Fe/H]>_s-poor ~ 0.2~dex; and (iii) the presence of O-Na-Al (anti-)correlations in both stellar groups. We have defined a new photometric index, c_{BVI}=(B-V)-(V-I), to maximise the separation in the colour-magnitude diagram between the two stellar groups with different Fe and s-element content, and this index is not significantly affected by variations in light elements (such as the O-Na anticorrelation). The variations in the overall metallicity present in NGC5286 add this object to the class of "anomalous" GCs. Furthermore, the chemical abundance pattern of NGC5286 resembles that observed in some of the anomalous GCs, e.g. M22, NGC1851, M2, and the more extreme Omega Centauri, that also show internal variations in s-elements, and in light elements within stars with different Fe and s-elements content. In view of the common variations in s-elements, we propose the term s-Fe-anomalous GCs to describe this sub-class of objects. The similarities in chemical abundance ratios between these objects strongly suggest similar formation and evolution histories, possibly associated with an origin in tidally disrupted dwarf satellites.

Stellar activity and, in particular, convection-related surface structures, potentially cause fluctuations that can affect the transit light curves. Surface convection simulations can help the interpretation of ToV. We used realistic three-dimensional radiative hydrodynamical simulation of the Sun from the Stagger-grid and synthetic images computed with the radiative transfer code Optim3D to provide predictions for the transit of Venus in 2004 observed by the satellite ACRIMSAT. We computed intensity maps from RHD simulation of the Sun and produced synthetic stellar disk image. We computed the light curve and compared it to the ACRIMSAT observations and also to the light curves obtained with solar surface representations carried out using radial profiles with different limb-darkening laws. We also applied the same spherical tile imaging method to the observations of center-to-limb Sun granulation with HINODE. We managed to explain ACRIMSAT observations of 2004 ToV and showed that the granulation pattern causes fluctuations in the transit light curve. We evaluated the contribution of the granulation to the ToV. We showed that the granulation pattern can partially explain the observed discrepancies between models and data. This confirms that the limb-darkening and the granulation pattern simulated in 3D RHD Sun represent well what is imaged by HINODE. In the end, we found that the Venus's aureole contribution during ToV is less intense than the solar photosphere, and thus negligible. Being able to explain consistently the data of 2004 ToV is a new step forward for 3D RHD simulations that are becoming essential for the detection and characterization of exoplanets. They show that the granulation have to be considered as an intrinsic incertitude, due to the stellar variability, on precise measurements of exoplanet transits of, most likely, planets with small diameters.

We compute the emergent stellar spectra from the UV to far infrared for different viewing angles using realistic 3D model atmospheres for a large range in stellar parameters to predict the stellar limb darkening. We have computed full 3D LTE synthetic spectra based on 3D radiative hydrodynamic atmosphere models from the Stagger-grid. From the resulting intensities at different wavelength, we derived coefficients for the standard limb darkening laws considering a number of often-used photometric filters. Furthermore, we calculated theoretical transit light curves, in order to quantify the differences between predictions by the widely used 1D model atmosphere and our 3D models. The 3D models are often found to predict steeper limb darkening compared to the 1D models, mainly due to the temperature stratifications and temperature gradients being different in the 3D models compared to those predicted with 1D models based on the mixing length theory description of convective energy transport. The resulting differences in the transit light curves are rather small; however, these can be significant for high-precision observations of extrasolar transits, and are able to lower the residuals from the fits with 1D limb darkening profiles. We advocate the use of the new limb darkening coefficients provided for the standard four-parameter non-linear power law, which can fit the limb darkening more accurately than other choices.

Context: The study of chemical abundance patterns in globular clusters is of key importance to constrain the different candidates for intra-cluster pollution of light elements. Aims: We aim at deriving accurate abundances for a large range of elements in the globular cluster 47 Tucanae (NGC 104) to add new constraints to the pollution scenarios for this particular cluster, expanding the range of previously derived element abundances. Methods: Using tailored 1D LTE atmospheric models together with a combination of equivalent width measurements, LTE, and NLTE synthesis we derive stellar parameters and element abundances from high-resolution, high signal-to-noise spectra of 13 red giant stars near the tip of the RGB. Results: We derive abundances of a total 27 elements (O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ru, Ba, La, Ce, Pr, Nd, Eu, Dy). Departures from LTE were taken into account for Na, Al and Ba. We find a mean [Fe/H] = $-0.78\pm0.07$ and $[\alpha/{\rm Fe}]=0.34\pm0.03$ in good agreement with previous studies. The remaining elements show good agreement with the literature, but the inclusion of NLTE for Al has a significant impact on the behaviour of this key element. Conclusions: We confirm the presence of an Na-O anti-correlation in 47 Tucanae found by several other works. Our NLTE analysis of Al shifts the [Al/Fe] to lower values, indicating that this may be overestimated in earlier works. No evidence for an intrinsic variation is found in any of the remaining elements.

The composition of the Sun is an essential piece of reference data for astronomy, cosmology, astroparticle, space and geo-physics. This article, dealing with the intermediate-mass elements Na to Ca, is the first in a series describing the comprehensive re-determination of the solar composition. In this series we severely scrutinise all ingredients of the analysis across all elements, to obtain the most accurate, homogeneous and reliable results possible. We employ a highly realistic 3D hydrodynamic solar photospheric model, which has successfully passed an arsenal of observational diagnostics. To quantify systematic errors, we repeat the analysis with three 1D hydrostatic model atmospheres (MARCS, MISS and Holweger & M\"{u}ller 1974) and a horizontally and temporally-averaged version of the 3D model ($\langle$3D$\rangle$). We account for departures from LTE wherever possible. We have scoured the literature for the best transition probabilities, partition functions, hyperfine and other data, and stringently checked all observed profiles for blends. Our final 3D+NLTE abundances are: $\log\epsilon_{\mathrm{Na}}=6.21\pm0.04$, $\log\epsilon_{\mathrm{Mg}}=7.59\pm0.04$, $\log\epsilon_{\mathrm{Al}}=6.43\pm0.04$, $\log\epsilon_{\mathrm{Si}}=7.51\pm0.03$, $\log\epsilon_{\mathrm{P}}=5.41\pm0.03$, $\log \epsilon_{\mathrm{S}}=7.13\pm0.03$, $\log\epsilon_{\mathrm{K}}=5.04\pm0.05$ and $\log\epsilon_{\mathrm{Ca}}=6.32\pm0.03$. The uncertainties include both statistical and systematic errors. Our results are systematically smaller than most previous ones with the 1D semi-empirical Holweger & M\"uller model. The $\langle$3D$\rangle$ model returns abundances very similar to the full 3D calculations. This analysis provides a complete description and a slight update of the Na to Ca results presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), with full details of all lines and input data.

Aims. We used realistic three-dimensional (3D) radiative hydrodynamical (RHD) simulations from the Stagger-grid and synthetic images computed with the radiative transfer code Optim3D to provide interferometric observables to extract the signature of stellar granulation and transiting planets. Methods. We computed intensity maps from RHD simulations for twelve interferometric instruments covering wavelengths ranging from optical to infrared. The stellar surface asymmetries in the brightness distribution mostly affect closure phases. We compared the closure phases of the system star with a transiting planet and the star alone and considered the impact of magnetic spots constructing a hypothetical starspots image. Results. All the simulations show departure from the axisymmetric case at all wavelengths. We presented two possible targets (Beta Com and Procyon) and found that departures up to 16 deg can be detected on the 3rd lobe and higher. In particular, MIRC is the most appropriate instrument because it combines good UV coverage and long baselines. Moreover, we explored the impact of convection on interferometric planet signature for three prototypes of planets. It is possible to disentangle the signature of the planet at particular wavelengths (either in the infrared or in the optical) by comparing the closure phases of the star at difference phases of the planetary transit. Conclusions. The detection and characterisation of planets must be based on a comprehensive knowledge of the host star; this includes the detailed study of the stellar surface convection with interferometric techniques. In this context, RHD simulations are crucial to reach this aim. We emphasize that interferometric observations should be pushed at high spatial frequencies by accumulating observations on closure phases at short and long baselines.

We present a theoretical study of the effects and signatures of realistic velocity field and atmospheric inhomogeneities associated with convective motions at the surface of cool late-type stars on the emergent profiles of iron spectral lines for a large range in stellar parameters. We compute 3D spectral line flux profiles under the assumption of local thermodynamic equilibrium (LTE) by employing state-of-the-art, time-dependent, 3D, radiative-hydrodynamical atmosphere models from the Stagger-grid. A set of 35 real unblended, optical FeI and FeII lines of varying excitation potential are considered. Additionally, fictitious Fe i and Fe ii lines (5000A and 0, 2, 4 eV) are used to construct general curves of growth and enable comparison of line profiles with the same line strength to illustrate systematical trends stemming from the intrinsic structural differences among 3D model atmospheres with different stellar parameters. Theoretical line shifts and bisectors are derived to analyze the shapes, shifts, and asymmetries imprinted in the full 3D line profiles emerging self-consistently from the convective simulations with velocity fields and atmospheric inhomogeneities. We find systematic variations in line strength, shift, width, and bisectors, that can be related to the respective physical conditions at the height of the line formation in the stellar atmospheric environment, in particular the amplitude of the vertical velocity field. Line shifts and asymmetries arise due to the presence of convective velocities and the granulation pattern that are ubiquitously found in observed stellar spectra of cool stars.

Theoretical atmosphere models provide the basis for a variety of applications in astronomy. In simplified one-dimensional (1D) atmosphere models, convection is usually treated with the mixing length theory despite its well-known insufficiency, therefore, the superadiabatic regime is poorly rendered. Due to the increasing computational power over the last decades, we are now capable to compute large grids of realistic three-dimensional (3D) hydrodynamical model atmospheres with the realistic treatment of the radiative transfer. We have computed the Stagger-grid, a comprehensive grid of 3D atmosphere models for late-type stars. In the presented contribution, we discuss initial results of the grid by exploring global properties and mean stratifications of the 3D models. Furthermore, we also depict the differences to classic 1D atmosphere models.

We study the implications of averaging methods with different reference depth scales for 3D hydrodynamical model atmospheres computed with the Stagger-code. The temporally and spatially averaged (hereafter denoted as <3D>) models are explored in the light of local thermodynamic equilibrium (LTE) spectral line formation by comparing spectrum calculations using full 3D atmosphere structures with those from <3D> averages. We explore methods for computing mean <3D> stratifications from the Stagger-grid time-dependent 3D radiative hydro- dynamical atmosphere models by considering four different reference depth scales (geometrical depth, column-mass density, and two optical depth scales). Furthermore, we investigate the influence of alternative averages (logarithmic or enforced hydrostatic equilibrium, flux-weighted temperatures). For the line formation we compute curves of growth for Fe i and Fe ii lines in LTE . The resulting <3D> stratifications for the four reference depth scales can be considerably different. We find typically that in the upper atmosphere and in the superadiabatic region just below the optical surface, where the temperature and density fluctuations are highest, the differences become considerable and increase for higher Teff, lower logg, and lower [Fe/H]. The differential comparison of spectral line formation shows distinctive differences depending on which <3D> model is applied. The averages over layers of constant column-mass density yield the best mean <3D> representation for LTE line formation, while the averages on layers at constant geometrical height are the least appropriate. Unexpectedly, the usually preferred averages over layers of constant optical depth are prone to the increasing interference of the reversed granulation towards higher effective temperature, in particular at low metallicity.

We study the effects of different descriptions of the solar surface convection on the eigenfrequencies of p-modes. 1-D evolution calculations of the whole Sun and 3-D hydrodynamic and magnetohydrodynamic simulations of the current surface are performed. These calculations rely on realistic physics. Averaged stratifications of the 3-D simulations are introduced in the 1-D solar evolution or in the structure models. The eigenfrequencies obtained are compared to those of 1-D models relying on the usual phenomenologies of convection and to observations of the MDI instrument aboard SoHO. We also investigate how the magnetic activity could change the eigenfrequencies and the solar radius, assuming that, 3 Mm below the surface, the upgoing plasma advects a 1.2 kG horizontal field. All models and observed eigenfrequencies are fairly close below 3 mHz. Above 3 mHz the eigenfrequencies of the phenomenological convection models are above the observed eigenfrequencies. The frequencies of the models based on the 3-D simulations are slightly below the observed frequencies. Their maximum deviation is ~ 3 \mu Hz at 3 mHz but drops below 1 \mu Hz at 4 mHz. Replacing the hydrodynamic by the magnetohydrodynamic simulation increases the eigenfrequencies. The shift is negligible below 2.2 mHz and then increases linearly with frequency to reach ~ 1.7 \mu Hz at 4 mHz. The impact of the simulated activity is a 14 milliarcsecond shrinking of the solar layers near the optical depth unity.

We report on a strictly differential line-by-line analysis of high quality UVES spectra of bright giants in the metal-poor globular cluster NGC 6752. We achieved high precision differential chemical abundance measurements for Fe, Na, Si, Ca, Ti, Cr, Ni, Zn, Y, Zr, Ba, La, Ce, Pr, Nd, Sm, Eu and Dy with uncertainties as low as ~0.01 dex (~2%). We obtained the following main results. (1) The observed abundance dispersions are a factor of ~2 larger than the average measurement uncertainty. (2) There are positive correlations, of high statistical significance, between all elements and Na. (3) For any pair of elements, there are positive correlations of high statistical significance, although the amplitudes of the abundance variations are small. Removing abundance trends with effective temperature and/or using a different pair of reference stars does not alter these results. These abundance variations and correlations may reflect a combination of (a) He abundance variations and (b) inhomogeneous chemical evolution in the pre- or proto-cluster environment. Regarding the former, the current constraints on Delta Y from photometry likely preclude He as being the sole explanation. Regarding the latter, the nucleosynthetic source(s) must have synthesised Na, alpha, Fe-peak and neutron-capture elements and in constant amounts for species heavier than Si; no individual object can achieve such nucleosynthesis. We speculate that other, if not all, globular clusters may exhibit comparable abundance variations and correlations to NGC 6752 if subjected to a similarly precise analysis.

We present the Stagger-grid, a comprehensive grid of time-dependent, 3D hydrodynamic model atmospheres for late-type stars with realistic treatment of radiative transfer, covering a wide range in stellar parameters. This grid of 3D models is intended for various applications like stellar spectroscopy, asteroseismology and the study of stellar convection. In this introductory paper, we describe the methods used for the computation of the grid and discuss the general properties of the 3D models as well as their temporal and spatial averages (<3D>). All our models were generated with the Stagger-code, using realistic input physics for the equation of state (EOS) and for continuous and line opacities. Our ~220 grid models range in Teff from 4000 to 7000K in steps of 500K, in log g from 1.5 to 5.0 in steps of 0.5 dex, and [Fe/H] from -4.0 to +0.5 in steps of 0.5 and 1.0 dex. We find a tight scaling relation between the vertical velocity and the surface entropy jump, which itself correlates with the constant entropy value of the adiabatic convection zone. The range in intensity contrast is enhanced at lower metallicity. The granule size correlates closely with the pressure scale height sampled at the depth of maximum velocity. We compare the <3D> models with widely applied 1D models, as well as with theoretical 1D hydrostatic models generated with the same EOS and opacity tables as the 3D models, in order to isolate the effects of using self-consistent and hydrodynamic modeling of convection, rather than the classical mixing length theory approach. For the first time, we are able to quantify systematically over a broad range of stellar parameters the uncertainties of 1D models arising from the simplified treatment of physics, in particular convective energy transport. In agreement with previous findings, we find that the differences can be significant, especially for metal-poor stars.

Un-evolved, very metal-poor stars are the most important tracers of the cosmic abundance of lithium in the early universe. Combining the standard Big Bang nucleosynthesis model with Galactic production through cosmic ray spallation, these stars at [Fe/H]<-2 are expected to show an undetectably small 6Li/7Li isotopic signature. Evidence to the contrary may necessitate an additional pre-galactic production source or a revision of the standard model of Big Bang nucleosynthesis. We revisit the isotopic analysis of four halo stars, two with claimed 6Li-detections in the literature, to investigate the influence of improved model atmospheres and line formation treatment. For the first time, a combined 3D, NLTE (non-local thermodynamic equilibrium) modelling technique for Li, Na, and Ca lines is utilised to constrain the intrinsic line-broadening and to determine the Li isotopic ratio. We discuss the influence of 3D NLTE effects on line profile shapes and assess the realism of our modelling using the Ca excitation and ionisation balance. By accounting for NLTE line formation in realistic 3D hydrodynamical model atmospheres, we can model the Li resonance line and other neutral lines with a consistency that is superior to LTE, with no need for additional line asymmetry caused by the presence of 6Li. Contrary to the results from 1D and 3D LTE modelling, no star in our sample has a significant 2-sigma detection of the lighter isotope in NLTE. Over a large parameter space, NLTE modelling systematically reduces the best-fit Li isotopic ratios by up to five percentage points. As a bi-product, we also present the first ever 3D NLTE Ca and Na abundances of halo stars, which reveal significant departures from LTE. The observational support for a significant and non-standard 6Li production source in the early universe is substantially weakened by our findings.

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.

Recently, new solar model atmospheres have been developed to replace classical 1D LTE hydrostatic models and used to for example derive the solar chemical composition. We aim to test various models against key observational constraints. In particular, a 3D model used to derive the solar abundances, a 3D MHD model (with an imposed 10 mT vertical magnetic field), 1D models from the PHOENIX project, the 1D MARCS model, and the 1D semi-empirical model of Holweger & M\"uller. We confront the models with observational diagnostics of the temperature profile: continuum centre-to-limb variations (CLV), absolute continuum fluxes, and the wings of hydrogen lines. We also test the 3D models for the intensity distribution of the granulation and spectral line shapes. The predictions from the 3D model are in excellent agreement with the continuum CLV observations, performing even better than the Holweger & M\"uller model (constructed largely to fulfil such observations). The predictions of the 1D theoretical models are worse, given their steeper temperature gradients. For the continuum fluxes, predictions for most models agree well with the observations. No model fits all hydrogen lines perfectly, but again the 3D model comes ahead. The 3D model also reproduces the observed continuum intensity fluctuations and spectral line shapes very well. The excellent agreement of the 3D model with the observables reinforces the view that its temperature structure is realistic. It outperforms the MHD simulation in all diagnostics, implying that recent claims for revised abundances based on MHD modelling are premature. Several weaknesses in the 1D models are exposed. The differences between the PHOENIX LTE and NLTE models are small. We conclude that the 3D hydrodynamical model is superior to any of the tested 1D models, which gives further confidence in the solar abundance analyses based on it.

Present grids of stellar atmosphere models are the workhorses in interpreting stellar observations, and determining their fundamental parameters. These models rely on greatly simplified models of convection, however, lending less predictive power to such models of late type stars. We present a grid of improved and more reliable stellar atmosphere models of late type stars, based on deep, 3D, convective, stellar atmosphere simulations. This grid is to be used in general for interpreting observations, and improve stellar and asteroseismic modeling. We solve the Navier Stokes equations in 3D and concurrent with the radiative transfer equation, for a range of atmospheric parameters, covering most of stellar evolution with convection at the surface. We emphasize use of the best available atomic physics for quantitative predictions and comparisons with observations. We present granulation size, convective expansion of the acoustic cavity, asymptotic adiabat, as function of atmospheric parameters. These and other results are also available in electronic form.

We have determined the angular diameters of two metal-poor stars, HD 122563 and Gmb 1830, using CHARA and Palomar Testbed Interferometer observations. For the giant star HD 122563, we derive an angular diameter theta_3D = 0.940 +- 0.011 milliarcseconds (mas) using limb-darkening from 3D convection simulations and for the dwarf star Gmb 1830 (HD 103095) we obtain a 1D limb-darkened angular diameter theta_1D = 0.679 +- 0.007 mas. Coupling the angular diameters with photometry yields effective temperatures with precisions better than 55 K (Teff = 4598 +- 41 K and 4818 +- 54 K --- for the giant and the dwarf star, respectively). Including their distances results in very well-determined luminosities and radii (L = 230 +- 6 L_sun, R = 23.9 +- 1.9 R_sun and L = 0.213 +- 0.002 L_sun, R = 0.664 +- 0.015 R_sun, respectively). We used the CESAM2k stellar structure and evolution code in order to produce models that fit the observational data. We found values of the mixing-length parameter alpha (which describes 1D convection) that depend on the mass of the star. The masses were determined from the models with precisions of <3% and with the well-measured radii excellent constraints on the surface gravity are obtained (log g = 1.60 +- 0.04, 4.59 +- 0.02, respectively). The very small errors on both log g and Teff provide stringent constraints for spectroscopic analyses given the sensitivity of abundances to both of these values. The precise determination of Teff for the two stars brings into question the photometric scales for metal-poor stars.

We investigate departures from LTE in the line formation of Fe for a number of well-studied late-type stars in different evolutionary stages. A new model of Fe atom was constructed from the most up-to-date theoretical and experimental atomic data available so far. Non-local thermodynamic equilibrium (NLTE) line formation calculations for Fe were performed using 1D hydrostatic MARCS and MAFAGS-OS model atmospheres, as well as the spatial and temporal average stratifications from full 3D hydrodynamical simulations of stellar convection computed using the Stagger code. It is shown that the Fe I/Fe II ionization balance can be well established with the 1D and mean 3D models under NLTE including calibrated inelastic collisions with H I calculated from the Drawin's (1969) formulae. Strong low-excitation Fe I lines are very sensitive to the atmospheric structure; classical 1D models fail to provide consistent excitation balance, particularly so for cool metal-poor stars. A better agreement between Fe I lines spanning a range of excitation potentials is obtained with the mean 3D models. Mean NLTE metallicities determined for the standard stars using the 1D and mean 3D models are fully consistent. Also, the NLTE spectroscopic effective temperatures and gravities from ionization balance agree with that determined by other methods, e.g., infrared flux method and parallaxes, if one of the stellar parameters is constrained independently.

We used a new realistic 3D radiative-hydrodynamical model atmosphere of Procyon generated with the Stagger Code and synthetic spectra computed with the radiative transfer code Optim3D to re-analyze interferometric and spectroscopic data from the optical to the infrared of Procyon. We compute intensity maps in two optical filters centered at 500 and 800 nm (MARK III) and one infrared filter centered at 2200 nm (VINCI). We constructed stellar disk images accounting for the center-to-limb variations and used them to derive visibility amplitudes and closure phases. We provide 3D limb-darkening coefficients in the optical as well as in the infrared. We show that visibility curves and closure phases show clear deviations from circular symmetry from the 3rd lobe on. These deviations are detectable with current interferometers using closure phases. We derive new angular diameters at different wavelengths with two independent methods based on 3D simulations. We find a diameter_Vinci = 5.390 \pm 0.03 mas that this is confirmed by an independent asteroseismic estimation. The resulting Teff is 6591 K, which is consistent with the infrared flux method determinations. We find also a value of the surface gravity log g = 4.01 \pm 0.03 that is larger by 0.05 dex from literature values. Spectrophotometric comparisons with observations provide very good agreement with the spectral energy distribution and photometric colors, allowing us to conclude that the thermal gradient of the simulation matches fairly well Procyon. Finally, we show that the granulation pattern of a planet hosting Procyon-like star has a non-negligible impact on the detection of hot Jupiters in the infrared using interferometry closure phases. It is then crucial to have a comprehensive knowledge of the host star to directly detect and characterize hot Jupiters. In this respect, RHD simulations are very important to reach this aim.

Radiative hydrodynamic simulations of solar and stellar surface convection have become an important tool for exploring the structure and gas dynamics in the envelopes and atmospheres of late-type stars and for improving our understanding of the formation of stellar spectra. We quantitatively compare results from three-dimensional, radiative hydrodynamic simulations of convection near the solar surface generated with three numerical codes CO5BOLD, MURaM, and STAGGER and different simulation setups in order to investigate the level of similarity and to cross-validate the simulations. For all three simulations, we considered the average stratifications of various quantities (temperature, pressure, flow velocity, etc.) on surfaces of constant geometrical or optical depth, as well as their temporal and spatial fluctuations. We also compared observables, such as the spatially resolved patterns of the emerging intensity and of the vertical velocity at the solar optical surface as well as the center-to-limb variation of the continuum intensity at various wavelengths. The depth profiles of the thermodynamical quantities and of the convective velocities as well as their spatial fluctuations agree quite well. Slight deviations can be understood in terms of differences in box size, spatial resolution and in the treatment of non-gray radiative transfer between the simulations. The results give confidence in the reliability of the results from comprehensive radiative hydrodynamic simulations.

Convection plays an essential role in the emerging intensity for many stars that will be observed by Gaia. Convective-related surface structures affect the shape, shift, and asymmetry of absorption lines, the phocentric and photometric variability causing bias in Gaia measurements. Regarding the importance of Gaia mission and its goals, it is mandatory to have the best models of the observed stars. 3D time-dependent hydrodynamical simulations of surface convection are crucial to model the photosphere of late type stars in a very realistic way. These simulations are an important tool to correct the radial velocities and better estimates the parallaxes and photometric variability.

In this contribution, we present the StaggerGrid, a collaborative project for the construction of a comprehensive grid of time-dependent, three-dimensional (3D), hydrodynamic model atmospheres of solar- and late-type stars with different effective temperatures, surface gravities, and chemical compositions. We illustrate the main characteristics of these 3D models and their effects on the predicted strengths, wavelength-shifts, and shapes of spectral lines, highlighting the differences with respect to calculations based on classical, one-dimensional, hydrostatic models, and discuss some of their possible applications to elemental abundance analysis of stellar spectra in the context of large observational surveys.

Bisectors of the HARPS cross-correlation function (CCF) can discern between planetary radial-velocity (RV) signals and spurious RV signals from stellar magnetic activity variations. However, little is known about the effects of the stellar atmosphere on CCF bisectors or how these effects vary with spectral type and luminosity class. Here we investigate the variations in the shapes of HARPS CCF bisectors across the HR diagram in order to relate these to the basic stellar parameters, surface gravity and temperature. We use archive spectra of 67 well studied stars observed with HARPS and extract mean CCF bisectors. We derive previously defined bisector measures (BIS, v_bot, c_b) and we define and derive a new measure called the CCF Bisector Span (CBS) from the minimum radius of curvature on direct fits to the CCF bisector. We show that the bisector measures correlate differently, and non-linearly with log g and T_eff. The resulting correlations allow for the estimation of log g and T_eff from the bisector measures. We compare our results with 3D stellar atmosphere models and show that we can reproduce the shape of the CCF bisector for the Sun.

We describe the first results from our project aimed at large-scale calculations of NLTE abundance corrections for important astrophysical atoms and ions. In this paper, the focus is on Fe which is a proxy of stellar metallicity and is commonly used to derive effective temperature and gravity. We present a small grid of NLTE abundance corrections for Fe I lines and discuss how NLTE effects influence determination of effective temperature, surface gravity, and metallicity for late-type stars.

We investigate the effects of coherent isotropic continuum scattering on the formation of spectral lines in local thermodynamic equilibrium (LTE) using 3D hydrodynamical and 1D hydrostatic model atmospheres of red giant stars. Continuum flux levels, spectral line profiles and curves of growth for different species are compared with calculations that treat scattering as absorption. Photons may escape from deeper, hotter layers through scattering, resulting in significantly higher continuum flux levels beneath a wavelength of 5000 A. The magnitude of the effect is determined by the importance of scattering opacity with respect to absorption opacity; we observe the largest changes in continuum flux at the shortest wavelengths and lowest metallicities; intergranular lanes of 3D models are more strongly affected than granules. Continuum scattering acts to increase the profile depth of LTE lines: continua gain more brightness than line cores due to their larger thermalization depth in hotter layers. We thus observe the strongest changes in line depth for high-excitation species and ionized species, which contribute significantly to photon thermalization through their absorption opacity near the continuum optical surface. Scattering desaturates the line profiles, leading to larger abundance corrections for stronger lines, which reach -0.5 dex at 3000 A for Fe II lines in 3D with excitation potential 2 eV at [Fe/H]=-3.0. The corrections are less severe for low-excitation lines, longer wavelengths, and higher metallicity. Velocity fields increase the effects of scattering by separating emission from granules and intergranular lanes in wavelength. 1D calculations exhibit similar scattering abundance corrections for weak lines, but those for strong lines are generally smaller compared to 3D models and depend on the choice of microturbulence.

It is still debated whether or not the Galactic chemical evolution of sulphur in the halo followed the constant or flat trend with [Fe/H], ascribed to the result of explosive nucleosynthesis in type II SNe. The aim of this study is to try to clarify this situation by measuring the sulphur abundance in a sample of halo giants using two diagnostics; the S I triplet around 1045 nm and the [S I] line at 1082 nm. The latter of the two is not believed to be sensitive to non-LTE effects. We can thereby minimize the uncertainties in the diagnostic used and estimate the usefulness of the triplet in sulphur determination in halo K giants. We will also be able to compare our sulphur abundance differences from the two diagnostics with the expected non-LTE effects in the 1045 nm triplet previously calculated by others. High-resolution near-infrared spectra of ten K giants were recorded using the spectrometer CRIRES mounted on VLT. Two standard settings were used; one covering the S I triplet and one covering the [S I] line. The sulphur abundances were determined individually with equivalent widths and synthetic spectra for the two diagnostics using tailored 1D model atmospheres and relying on non-LTE corrections from the litterature. Effects of convective inhomogeneities in the stellar atmospheres are investigated. We corroborate the flat trend in the [S/Fe] vs. [Fe/H] plot for halo stars found in other works and cannot find a scatter nor a rise in [S/Fe] obtained in some other previous studies. We find the sulphur abundances deduced from the non-LTE corrected triplet somewhat lower than the abundances from the [S I] line, possibly indicating too large non-LTE corrections. Considering 3D modeling, however, they might instead be too small. Further we show that the [S I] line is possible to use as a sulphur diagnostic down to [Fe/H] = -2.3 in giants.

Three-dimensional (3D) radiative hydrodynamic model atmospheres of metal-poor late-type stars are characterized by cooler upper photospheric layers than their 1D counterparts. This property of 3D models can dramatically affect elemental abundances derived from temperature-sensitive spectral lines. We investigate whether the cool surface temperatures predicted by metal-poor 3D models can be ascribed to the approximated treatment of scattering in the radiative transfer. We use the Bifrost code to test three different ways to handle scattering in 3D model atmospheres of metal-poor stars. First, we solve self-consistently the radiative transfer equation for a source function with a coherent scattering term. Second, we solve the radiative transfer equation for a Planckian source function, neglecting the contribution of continuum scattering to extinction in the optically thin layers; this has been the default mode in previous models of ours. Third, we treat scattering as pure absorption everywhere, which is the standard case in CO5BOLD models. We find that the second approach produces temperature structures with cool upper photospheric layers very similar to the correct coherent scattering solution. In contrast, treating scattering as pure absorption leads to significantly hotter and shallower temperature stratifications. The main differences in temperature structure between our published models and those generated with the CO5BOLD code can be traced to the different treatments of scattering. Neglecting the contribution of continuum scattering to extinction in optically thin layers provides a good approximation to the full radiative transfer solution for metal-poor stars. Our results demonstrate that the cool temperature stratifications predicted for metal-poor late-type stellar atmospheres by previous models of ours are not an artifact of the approximated treatment of scattering.

A very high resolution (R=200,000), high signal-to-noise ratio (S/N=340) blue-green spectrum of the very metal-poor ([Fe/H]=-2.6) red giant star HD122563 has been obtained by us at McDonald Observatory. We measure the asymmetries and core wavelengths of a set of unblended FeI lines covering a wide range of line strength. Line bisectors exhibit the characteristic C-shape signature of surface convection (granulation) and they span from about 100 m/s in the strongest FeI features to 800 m/s in the weakest ones. Core wavelength shifts range from about -100 to -900 m/s, depending on line strength. In general, larger blueshifts are observed in weaker lines, but there is increasing scatter with increasing residual flux. Assuming local thermodynamic equilibrium (LTE), we synthesize the same set of spectral lines using a state-of-the-art three-dimensional hydrodynamic simulation for a stellar atmosphere of fundamental parameters similar to those of HD122563. We find good agreement between model predictions and observations. This allows us to infer an absolute zero-point for the line shifts and radial velocity. Moreover, it indicates that the structure and dynamics of the simulation are realistic, thus providing support to previous claims of large 3D-LTE corrections, based on the hydrodynamic model used here, to elemental abundances and fundamental parameters of very metal-poor red giant stars obtained with standard 1D-LTE spectroscopic analyses.

Context. Theoretical predictions from models of red giant branch stars are a valuable tool for various applications in astrophysics ranging from galactic chemical evolution to studies of exoplanetary systems. Aims. We use the radiative transfer code OPTIM3D and realistic 3D radiative-hydrodynamical (RHD) surface convection simulations of red giants to explore the impact of granulation on interferometric observables. Methods. We compute intensity maps for the 3D simulation snapshots in two filters: in the optical at 5000 \pm 300 {\AA} and in the K band 2.14 $\pm$ 0.26 {\mu}m FLUOR filter, corresponding to the wavelength-range of instruments mounted on the CHARA interferometer. From the intensity maps, we construct images of the stellar disks, accounting for center-to-limb variations. We then derive interferometric visibility amplitudes and phases. We study their behavior with position angle and wavelength. Results. We provide average limb-darkening coefficients for different metallicities and wavelength-ranges. We detail the prospects for the detection and characterization of granulation and center-to-limb variations of red giant stars with today's interferometers. We find that the effect of convective-related surface structures depends on metallicity and surface gravity. We provided theoretical closure phases that should be incorporated into the analysis of red giant planet companion closure phase signals. We estimate 3D-1D corrections to stellar radii determination: 3D models are ~ 3.5% smaller to ~ 1% larger in the optical with respect to 1D, and roughly 0.5 to 1.5% smaller in the infrared. Even if these corrections are small, they are important to properly set the zero point of effective temperature scale derived by interferometry and to strengthen the confidence of existing red giant catalogues of calibrating stars for interferometry.

We present the implementation of a radiative transfer solver with coherent scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD) simulations of stellar surface convection. The code is fully parallelized using MPI domain decomposition, which allows for large grid sizes and improved resolution of hydrodynamical structures. We apply the code to simulate the surface granulation in a solar-type star, ignoring magnetic fields, and investigate the importance of coherent scattering for the atmospheric structure. A scattering term is added to the radiative transfer equation, requiring an iterative computation of the radiation field. We use a short-characteristics-based Gauss-Seidel acceleration scheme to compute radiative flux divergences for the energy equation. The effects of coherent scattering are tested by comparing the temperature stratification of three 3D time-dependent hydrodynamical atmosphere models of a solar-type star: without scattering, with continuum scattering only, and with both continuum and line scattering. We show that continuum scattering does not have a significant impact on the photospheric temperature structure for a star like the Sun. Including scattering in line-blanketing, however, leads to a decrease of temperatures by about 350\,K below log tau < -4. The effect is opposite to that of 1D hydrostatic models in radiative equilibrium, where scattering reduces the cooling effect of strong LTE lines in the higher layers of the photosphere. Coherent line scattering also changes the temperature distribution in the high atmosphere, where we observe stronger fluctuations compared to a treatment of lines as true absorbers.

One of the most important factors in determining the stellar lithium abundance is the effective temperature. In a previous study by the authors, new effective temperatures Teff for sixteen metal-poor halo dwarfs were derived using a local thermodynamic equilibrium (LTE) description of the formation of Fe lines. This new Teff scale reinforced the discrepancy. For six of the stars from our previous study we calculate revised temperatures using a non-local thermodynamic equilibrium (NLTE) approach. These are then used to derive a new mean primordial lithium abundance in an attempt to solve the lithium discrepancy. Using the code MULTI we calculate NLTE corrections to the LTE abundances for the Fe I lines measured in the six stars, and determine new Teff's. We keep other physical parameters, i.e. log g, [Fe/H] and xi, constant at the values calculated in Paper I. With the revised Teff scale we derive new Li abundances. We compare the NLTE values of Teff with the photometric temperatures of Ryan et al. (1999, ApJ, 523, 654), the infrared flux method (IRFM) temperatures of Melendez & Ramirez (2004, ApJ, 615, 33), and the Balmer line wing temperatures of Asplund et al. (2006, ApJ, 644, 229). We find that our temperatures are hotter than both the Ryan et al. and Asplund et al. temperatures by typically ~ 110 K - 160 K, but are still cooler than the temperatures of Melendez & Ramirez by typically ~ 190 K. The temperatures imply a primordial Li abundance of 2.19 dex or 2.21 dex, depending on the magnitude of collisions with hydrogen in the calculations, still well below the value of 2.72 dex inferred from WMAP + BBN. We discuss the effects of collisions on trends of 7Li abundances with [Fe/H] and Teff, as well as the NLTE effects on the determination of log g through ionization equilibrium, which imply a collisional scaling factor SH > 1 for collisions between Fe and H atoms.

We present a preliminary local thermodynamic equilibrium (LTE) abundance analysis of the template halo red giant HD122563 based on a realistic, three-dimensional (3D), time-dependent, hydrodynamical model atmosphere of the very metal-poor star. We compare the results of the 3D analysis with the abundances derived by means of a standard LTE analysis based on a classical, 1D, hydrostatic model atmosphere of the star. Due to the different upper photospheric temperature stratifications predicted by 1D and 3D models, we find large, negative, 3D-1D LTE abundance differences for low-excitation OH and Fe I lines. We also find trends with lower excitation potential in the derived Fe LTE abundances from Fe I lines, in both the 1D and 3D analyses. Such trends may be attributed to the neglected departures from LTE in the spectral line formation calculations.

We derive the mixture of odd to even barium isotopes in the atmosphere of the metal-poor subgiant HD140283 from the analysis of the Ba II transition at 455.4 nm in a high-resolution high signal-to-noise spectrum of the star. The detailed shape of this spectral line depends on the relative contributions of odd and even isotopes via isotopic and hyperfine splitting. We measure the fractional abundance of odd Ba isotopes by modelling the formation of the Ba II 455.4 nm line profile with the use of both a classical 1D hydrostatic and a 3D hydrodynamical model atmosphere of HD140283. We interpret the results in terms of contributions by the slow (s-) and rapid (r-) neutron-capture processes to the isotopic mix. While the result of the 1D analysis of the Ba II feature indicates a (64 +/- 36)% contribution of the r-process to the isotopic mix, the 3D analysis points toward a mere (15 +/- 34)% contribution from this process, that is consistent with a solar-like mixture of barium isotopes.

I present here the main results of recent realistic, 3D, hydrodynamical simulations of convection at the surface of metal-poor red giant stars. I discuss the application of these convection simulations as time-dependent, 3D, hydrodynamical model atmospheres to spectral line formation calculations and abundance analyses. The impact of 3D models on derived elemental abundances is investigated by means of a differential comparison of the line strengths predicted in 3D under the assumption of local thermodynamic equilibrium (LTE) with the results of analogous line formation calculations performed with classical, 1D, hydrostatic model atmospheres. The low surface temperatures encountered in the upper photospheric layers of 3D model atmospheres of very metal-poor stars cause spectral lines of neutral metals and molecules to appear stronger in 3D than in 1D calculations. Hence, 3D elemental abundances derived from such lines are significantly lower than estimated by analyses with 1D models. In particular, differential 3D$-$1D LTE abundances for C, N, and O derived from CH, NH, and OH lines are found to be in the range -0.5 to -1 dex. Large negative differential 3D-1D corrections to the Fe abundance are also computed for weak low-excitation neutral Fe lines. The application of metal-poor 3D models to the spectroscopic analysis of extremely iron-poor halo stars is discussed.

We present a new abundance analysis of HE 1327-2326, the currently most iron-poor star, based on observational data obtained with VLT/UVES. We correct the 1D LTE abundances for 3D effects to provide an abundance pattern that supersedes previous works, and should be used to observationally test current models of the chemical yields of the first-generation SNe. Apart from confirming the 1D LTE abundances found in previous studies before accounting for 3D effects, we make use of a novel technique to apply the 3D-1D corrections for CNO which are a function of excitation potential and line strength for the molecular lines that comprise the observable CH, NH, and OH features. We find that the fit to the NH band at 3360 A is greatly improved due to the application of the 3D-1D corrections. This may indicate that 3D effects are actually observable in this star. We also report the first detection of several weak Ni lines. The cosmologically important element Li is still not detected; the new Li upper limit is extremely low, A(Li)<0.62, and in stark contrast with results not only from WMAP but also from other metal-poor stars. We also discuss how the new corrected abundance pattern of HE 1327-2326 is being reproduced by individual and integrated yields of SNe.

We present the results of realistic, 3D, hydrodynamical, simulations of surface convection in red giant stars with varying effective temperatures and metallicities. We use the convection simulations as time-dependent, hydrodynamical, model atmospheres to compute spectral line profiles for a number of ions and molecules under the assumption of local thermodynamic equilibrium (LTE). We compare the results with the predictions of line formation calculations based on 1D, hydrostatic, model stellar atmospheres in order to estimate the impact of 3D models on the derivation of elemental abundances. We find large negative 3D-1D LTE abundance corrections (typically -0.5 to -1 dex) for weak low-excitation lines from molecules and neutral species in the very low metallicity cases. Finally, we discuss the extent of departures from LTE in the case of neutral iron spectral line formation.

We present a homogeneous photometric and spectroscopic analysis of 18 stars along the evolutionary sequence of the metal-poor globular cluster NGC 6397 ([Fe/H] = -2), from the main-sequence turnoff point to red giants below the bump. The spectroscopic stellar parameters, in particular stellar-parameter differences between groups of stars, are in good agreement with broad-band and Stroemgren photometry calibrated on the infrared-flux method. The spectroscopic abundance analysis reveals, for the first time, systematic trends of iron abundance with evolutionary stage. Iron is found to be 31% less abundant in the turnoff-point stars than in the red giants. An abundance difference in lithium is seen between the turnoff-point and warm subgiant stars. The impact of potential systematic errors on these abundance trends (stellar parameters, the hydrostatic and LTE approximations) is quantitatively evaluated and found not to alter our conclusions significantly. Trends for various elements (Li, Mg, Ca, Ti and Fe) are compared with stellar-structure models including the effects of atomic diffusion and radiative acceleration. Such models are found to describe the observed element-specific trends well, if extra (turbulent) mixing just below the convection zone is introduced. It is concluded that atomic diffusion and turbulent mixing are largely responsible for the sub-primordial stellar lithium abundances of warm halo stars. Other consequences of atomic diffusion in old metal-poor stars are also discussed.

We investigate the impact of 3D hydrodynamical model atmospheres of red giant stars at different metallicities on the formation of spectral lines of a number of ions and molecules. We carry out realistic 3D simulations of surface convection in red giant stars with varying stellar parameters. We use the simulations as time-dependent hydrodynamical model stellar atmospheres to compute atomic (Li, O, Na, Mg, Ca, Fe) and molecular (CH, NH, OH) spectral lines under the assumption of local thermodynamic equilibrium (LTE). We compare the line strengths computed in 3D with the results of analogous line formation calculations for 1D, hydrostatic, plane-parallel MARCS model atmospheres in order to estimate the impact of 3D models on the derivation of elemental abundances. The temperature and density inhomogeneities and correlated velocities in 3D models, as well as the differences between the 1D and mean 3D structures significantly affect the predicted line strengths. Under the assumption of LTE, the low atmospheric temperatures of very metal-poor 3D model atmospheres cause the lines from neutral species and molecules to appear stronger than in 1D. Therefore, elemental abundances derived from these lines using 3D models are significantly lower than according to 1D analyses. Differences between 3D and 1D abundances of C, N, and O derived from CH, NH, and OH weak low-excitation lines are found to be in the range -0.5 dex to -1.0 dex for the the red giant stars at [Fe/H]=-3 considered here. At this metallicity, large negative corrections (about -0.8 dex) are also found for weak low-excitation Fe I lines. We caution, however, that departures from LTE might be significant for these and other elements and comparable to the effects due to stellar granulation.

This paper presents abundances for 12 red giants of the old open cluster Collinder 261 based on spectra from VLT/UVES. Abundances were derived for Na, Mg, Si, Ca, Mn, Fe, Ni, Zr and Ba. We find the cluster has a solar-level metallicity of [Fe/H] = -0.03 dex. However some alpha elements were found to be enhanced. The star-to-star scatter was consistent with the expected measurement uncertainty for all elements. The observed rms scatter is as follows: Na = 0.07, Mg = 0.05, Si = 0.06, Ca = 0.05, Mn = 0.03, Fe = 0.02, Ni = 0.04, Zr = 0.12, and Ba = 0.03 dex. The intrinsic scatter was estimated to be less than 0.05 dex. Such high levels of homogeneity indicate that chemical information remains preserved in this old open cluster. We use the chemical homogeneity we have now established in Cr 261, Hyades and the HR1614 moving group to examine the uniqueness of the individual cluster abundance patterns, ie. chemical signatures. We demonstrate that the three studied clusters have unique chemical signatures, and discuss how other such signatures may be searched for in the future. Our findings support the prospect of chemically tagging disk stars to common formation sites in order to unravel the dissipative history of the Galactic disk.

A homogeneous spectroscopic analysis of unevolved and evolved stars in the metal-poor globular cluster NGC 6397 with FLAMES-UVES reveals systematic trends of stellar surface abundances that are likely caused by atomic diffusion. This finding helps to understand, among other issues, why the lithium abundances of old halo stars are significantly lower than the abundance found to be produced shortly after the Big Bang.

The measurement of the cosmic microwave background has strongly constrained the cosmological parameters of the Universe. When the measured density of baryons (ordinary matter) is combined with standard Big Bang nucleosynthesis calculations, the amounts of hydrogen, helium and lithium produced shortly after the Big Bang can be predicted with unprecedented precision. The predicted primordial lithium abundance is a factor of two to three higher than the value measured in the atmospheres of old stars. With estimated errors of 10 to 25%, this cosmological lithium discrepancy seriously challenges our understanding of stellar physics, Big Bang nucleosynthesis or both. Certain modifications to nucleosynthesis have been proposed, but found experimentally not to be viable. Diffusion theory, however, predicts atmospheric abundances of stars to vary with time, which offers a possible explanation of the discrepancy. Here we report spectroscopic observations of stars in the metalpoor globular cluster NGC 6397 that reveal trends of atmospheric abundance with evolutionary stage for various elements. These element-specific trends are reproduced by stellar-evolution models with diffusion and turbulent mixing. We thus conclude that diffusion is predominantly responsible for the low apparent stellar lithium abundance in the atmospheres of old stars by transporting the lithium deep into the star.

We investigate the impact of realistic 3D hydrodynamical model stellar atmospheres on the determination of elemental abundances in the carbon-rich, hyper iron-poor stars HE0107-5240 and HE1327-2326. We derive the chemical compositions of the two stars by means of a detailed 3D analysis of spectral lines under the assumption of local thermodynamic equilibrium (LTE). The lower temperatures of the line-forming regions of the hydrodynamical models cause changes in the predicted spectral line strengths. In particular we find the 3D abundances of C, N, and O to be lower by ~ -0.8 dex (or more) than estimated from a 1D analysis. The 3D abundances of iron-peak elements are also decreased but by smaller factors (~ -0.2 dex). We caution however that the neglected non-LTE effects might actually be substantial for these metals. We finally discuss possible implications for studies of early Galactic chemical evolution.

The effects of background line opacity (line-blocking) in statistical equilibrium calculations for Fe in late-type stellar atmospheres have been investigated using an extensive and up-to-date model atom with radiative data primarily from the IRON Project. The background metal line opacities have been computed using data from the MARCS stellar model atmospheres. While accounting for this line opacity is important at solar metallicity, the differences between calculations including and excluding line-blocking at low metallicity are insignificant for the non-local thermodynamic equilibrium (non-LTE) abundance corrections for Fe I lines. The line-blocking has no impact on the non-LTE effects of Fe II lines. The dominant uncertainty in Fe non-LTE calculations for metal-poor stars is still the treatment of the inelastic H I collisions, which have here been included using scaling factors to the classical Drawin formalism, and whether or not thermalisation of the high Fe I levels to Fe II ground state should be enforced. Without such thermalisation, the Fe I non-LTE abundance corrections are substantial in metal-poor stars: about 0.3 dex with efficient (i.e. Drawin-like) H I collisions and about 0.5 dex without. Without both thermalisation and H I collisions, even Fe II lines show significant non-LTE effects in such stars.

Recently 3D hydrodynamical simulations of stellar surface convection have become feasible thanks to advances in computer technology and efficient numerical algorithms. Available observational diagnostics indicate that these models are highly realistic in describing the topology of stellar granulation and for spectral line formation purposes. The traditional free parameters (mixing length parameters, micro- and macroturbulence) always inherent in standard 1D analyses have thus finally become obsolete. These 3D models can therefore both shed light on the elusive nature of stellar convection as well as be employed in element abundance analyses. In the present contribution we will describe some aspects of the models and possible applications of them in terms of radiative transfer.