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.

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.

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.

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.

Solar p-mode oscillations exhibit a systematic offset towards higher frequencies due to shortcomings in the 1D stellar structure models, especially, the lack of turbulent pressure in the superadiabatic layers just below the optical surface, arising from the convective velocity field. We study the influence of the turbulent expansion, chemical composition, and magnetic fields on the stratification in the upper layers of the solar models in comparison with solar observations. Furthermore, we test alternative <3D> averages for improved results on the oscillation frequencies. We appended temporally and spatially averaged <3D> stratifications to 1D models to compute adiabatic oscillation frequencies that we then tested against observations. We also developed depth-dependent corrections for the solar 1D model, for which we expanded the geometrical depth to match the pressure stratification of the solar <3D> model, and we reduced the density that is caused by the turbulent pressure. We obtain the same results with our <3D> models as have been reported previously. Our depth-dependent corrected 1D models match the observations to almost a similar extent as the <3D> model. We find that correcting for the expansion of the geometrical depth and the reducing of the density are both equally necessary. Interestingly, the influence of the adiabatic exponent Gam1 is less pronounced than anticipated. The turbulent elevation directly from the <3D> model does not match the observations properly. Considering different reference depth scales for the <3D> averaging leads to very similar frequencies. Solar models with high metal abundances in their initial chemical composition match the low-frequency part much better. We find a linear relation between the p-mode frequency shift and the vertical magnetic field strength with dvnl = 26.21Bz microHz/kG, which is able to render the solar activity cycles correctly.

Stellar structure calculations are able to predict precisely the properties of stars during their evolution. However, convection is still modelled by the mixing length theory; therefore, the upper boundary conditions near the optical surface do not agree with asteroseismic observations. We want to improve how the outer boundary conditions are determined in stellar structure calculations. We study realistic 3D stellar atmosphere models to find alternative properties. We find that the asymptotic entropy run of the superadiabatic convective surface layers exhibit a distinct universal stratification when normalised by the entropy minimum and jump. The normalised entropy can be represented by a 5th order polynomial very accurately, and a 3rd order polynomial also yields accurate coefficients. This generic entropy stratification or the solar stratification, when scaled by the entropy jump and minimum, can be used to improve the modelling of superadiabatic surface layers in stellar structure calculations. Furthermore, this finding indicates that surface convection operates in the same way for all cool stars, but requires further scrutiny in order to improve our understanding of stellar atmospheres.

SMSS J031300.36-670839.3 (hereafter SM0313-6708) is a sub-giant halo star, with no detectable Fe lines and large overabundances of C and Mg relative to Ca. We obtained VLT-UVES spectra extending to 3060 Angstroms showing strong OH A-X band lines enabling an oxygen abundance to be derived. The OH A-X band lines in SM0313-6708 are much stronger than the CH C-X band lines. Spectrum synthesis fits indicate an [O/C] ratio of 0.02 +- 0.175. Our high S/N UVES data also enabled us to lower the Fe abundance limit to [Fe/H]{3D},NLTE < -7.52 (3 sigma). These data support our previous suggestion that the star formed from the iron-poor ejecta of a single massive star Population III supernova.

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.

We investigate the relation between 1D atmosphere models that rely on the mixing length theory and models based on full 3D radiative hydrodynamic (RHD) calculations to describe convection in the envelopes of late-type stars. The adiabatic entropy value of the deep convection zone, s_bot, and the entropy jump, {\Delta}s, determined from the 3D RHD models, are matched with the mixing length parameter, {\alpha}_MLT, from 1D hydrostatic atmosphere models with identical microphysics (opacities and equation-of-state). We also derive the mass mixing length, {\alpha}_m, and the vertical correlation length of the vertical velocity, C[v_z,v_z], directly from the 3D hydrodynamical simulations of stellar subsurface convection. The calibrated mixing length parameter for the Sun is {\alpha}_MLT (s_bot) = 1.98. For different stellar parameters, {\alpha}_MLT varies systematically in the range of 1.7 - 2.4. In particular, {\alpha}_MLT decreases towards higher effective temperature, lower surface gravity and higher metallicity. We find equivalent results for {\alpha}_MLT ({\Delta}s). Also, we find a tight correlation between the mixing length parameter and the inverse entropy jump. We derive an analytical expression from the hydrodynamic mean field equations that motivates the relation to the mass mixing length, {\alpha}_m, and find that it exhibits qualitatively a similar variation with stellar parameter (between 1.6 and 2.4) with a solar value of {\alpha}_m = 1.83. The vertical correlation length scaled with the pressure scale height yields for the Sun 1.71, but displays only a small systematic variation with stellar parameters, the correlation length slightly increasing with Teff. We derive mixing length parameters for various stellar parameters that can be used to replace a constant value. Within any convective envelope, {\alpha}_m and related quantities vary a lot.

In the surface layers of late-type stars, stellar convection is manifested with its typical granulation pattern due to the presence of convective motions. The resulting photospheric up- and downflows leave imprints in the observed spectral line profiles. We perform a careful statistical analysis of stellar granulation and its properties for different stellar parameters. We employ realistic 3D radiative hydrodynamic (RHD) simulations of surface convection from the Stagger-grid, a comprehensive grid of atmosphere models that covers a large parameter space in terms of Teff, logg, and [Fe/H]. Individual granules are detected from the (bolometric) intensity maps at disk center with an efficient granulation pattern recognition algorithm. From these we derive their respective properties: diameter, fractal dimension (area-perimeter relation), geometry, topology, variation of intensity, temperature, density and velocity with granule size. Also, the correlation of the physical properties at the optical surface are studied. We find in all of our 3D RHD simulations stellar granulation patterns imprinted, which are qualitatively similar to the solar case, despite the large differences in stellar parameters. The granules exhibit a large range in size, which can be divided into two groups - smaller and larger granules - by the mean granule size. These are distinct in their properties: smaller granules are regular shaped and dimmer, while the larger ones are increasingly irregular and more complex in their shapes and distribution in intensity contrast. This is reflected in their fractal dimensions, which is close to unity for the smaller granules, and close to two for larger granules, which is due to the fragmentation of granules. Stellar surface convection seems to operate scale-invariant over a large range in stellar parameters, which translates into a self-similar stellar granulation pattern.

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.

The element abundance ratios of four low-mass stars with extremely low metallicities indicate that the gas out of which the stars formed was enriched in each case by at most a few, and potentially only one low-energy, supernova. Such supernovae yield large quantities of light elements such as carbon but very little iron. The dominance of low-energy supernovae is surprising, because it has been expected that the first stars were extremely massive, and that they disintegrated in pair-instability explosions that would rapidly enrich galaxies in iron. What has remained unclear is the yield of iron from the first supernovae, because hitherto no star is unambiguously interpreted as encapsulating the yield of a single supernova. Here we report the optical spectrum of SMSS J031300.36- 670839.3, which shows no evidence of iron (with an upper limit of 10^-7.1 times solar abundance). Based on a comparison of its abundance pattern with those of models, we conclude that the star was seeded with material from a single supernova with an original mass of ~60 Mo (and that the supernova left behind a black hole). Taken together with the previously mentioned low-metallicity stars, we conclude that low-energy supernovae were common in the early Universe, and that such supernovae yield light element enrichment with insignificant iron. Reduced stellar feedback both chemically and mechanically from low-energy supernovae would have enabled first-generation stars to form over an extended period. We speculate that such stars may perhaps have had an important role in the epoch of cosmic reionization and the chemical evolution of early galaxies.

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.

Exoplanetary transmission spectroscopy in the near-infrared using Hubble/NICMOS is currently ambiguous because different observational groups claim different results from the same data, depending on their analysis methodologies. Spatial scanning with Hubble/WFC3 provides an opportunity to resolve this ambiguity. We here report WFC3 spectroscopy of the giant planets HD209458b and XO-1b in transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decorrelating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels collected in spatial scan mode. Our errors are within 6-percent (XO-1) and 26-percent (HD209458b) of the photon-limit at a spectral resolving power of 70, and are better than 0.01-percent per spectral channel. Both planets exhibit water absorption of approximately 200 ppm at the water peak near 1.38 microns. Our result for XO-1b contradicts the much larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for HD209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of an exoplanet atmosphere by Charbonneau et al. (2002). Model atmospheres having uniformly-distributed extra opacity of 0.012 cm^2 per gram account approximately for both our water measurement and the sodium absorption in this planet. Our results for HD209458b support the picture advocated by Pont et al. (2013) in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by continuous opacity due to haze and/or dust. However, the extra opacity needed for HD209458b is grayer than for HD189733b, with a weaker Rayleigh component.

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.

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.

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.

The open cluster M67 has solar metallicity and an age of about 4Gyr. The turn-off mass is close to the minimum mass for which solar metallicity stars develop a convective core during main sequence evolution as a result of the development of hydrogen burning through the CNO-cycle. The morphology of the color-magnitude-diagram (CMD) of M67 around the turn-off shows a clear hook-like feature, direct sign that stars close to the turn-off have convective cores. VandenBerg et al. investigated the possibility of using the morphology of the M67 turn-off to put constraints on the solar metallicity, particularly CNO elements, for which solar abundances have been revised downwards by more than 30% over the last few years. Here, we extend their work filling in the gaps in their analysis. To this aim, we compute isochrones appropriate for M67 using new (low metallicity) and old (high metallicity) solar abundances and study whether the characteristic turn-off in the CMD of M67 can be reproduced or not. We also study the importance of other constitutive physics on determining the presence of such a hook, particularly element diffusion, overshooting and nuclear reaction rates. We find that using the new solar abundance determinations, with low CNO abundances, makes it more difficult to reproduce the characteristic CMD of M67. This result is in agreement with results by VandenBerg et al. However, changes in the constitutive physics of the models, particularly overshooting, can influence and alter this result to the extent that isochrones constructed with models using low CNO solar abundances can also reproduce the turn-off morphology in M67. We conclude that only if all factors affecting the turn-off morphology are completely under control (and this is not the case), M67 could be used to put constraints on solar abundances.