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.

Andrea Chiavassa1, L. Casagrande2, R. Collet3, Z. Magic, L. Bigot1, F. Thévenin1 and M. Asplund2 Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Lagrange, CS 34229, Nice, France email: andrea.chiavassa@oca.eu Research School of Astronomy & Astrophysics, Australian National University, Cotter Road, Weston ACT 2611, Australia Stellar Astrophysics Centre, Department of Physics and Astronomy, Ny Munkegade 120, Aarhus University, DK-8000 Aarhus C, Denmark

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.

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 interferometer. 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 (HR 7633) = 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.

From exploratory studies and theoretical expectations it is known that simplifying approximations in spectroscopic analysis (local thermodynamic equilibrium (LTE), 1D) lead to systematic biases of stellar parameters and abundances. These biases depend strongly on surface gravity, temperature and, in particular, for LTE versus non-LTE (NLTE), on metallicity of the stars. Here we analyze 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, the other three methods display increasing biases toward 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 disk stars show a clearer behavior 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 disk stars and this population extends down to at least [Fe/H] similar to -1.6 dex, further than previously proven. The thick disk stars display a constant [Mg/Fe] approximate to 0.3 dex, with a small intrinsic dispersion in [Mg/Fe] that suggests that a fast SN Ia channel is not relevant for the disk 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 six Gaia benchmark stars using theoretical one-dimensional (1D) hydrostatic model atmospheres, as well as temporally and spatially averaged three-dimensional ( ) model atmospheres. The stars cover a range of from 4700 to 6500 K, from 1.6 to 4.4 dex, and 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 in and in 1D, (c) the solar Mg abundance is (total error), in excellent agreement with the Mg abundance measured in CI chondritic meteorites, (d) the 1D NLTE and NLTE approaches can be used with confidence to analyze optical Mg i lines in spectra of dwarfs and sub-giants, but for red giants the Mg i 5711 Å line should be preferred, (e) low-excitation Mg i lines are sensitive to the atmospheric structure; for these lines, LTE calculations with 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.

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.

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

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.

Aims. 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.Methods. The adiabatic entropy value of the deep convection zone, s(bot), and the entropy jump,Delta s, determined from the 3D RHD models, were matched with the mixing length parameter, alpha(MLT), from 1D hydrostatic atmosphere models with identical microphysics (opacities and equation-of-state). We also derived the mass mixing length parameter, 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.Results. The calibrated mixing length parameter for the Sun is alpha(circle dot)(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(circle dot)(MLT) (Delta s). In addition, 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 parameter, alpha(m), and find that it qualitatively shows a similar variation with stellar parameter (between 1.6 and 2.4) with the solar value of alpha(circle dot)(m) = 1.83. The vertical correlation length scaled with the pressure scale height yields 1.71 for the Sun, but only displays a small systematic variation with stellar parameters, the correlation length slightly increases with T-eff.Conclusions. 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 strongly. Our results will help to replace a constant alpha(MLT).

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.

Aims. We present a theoretical study of the effects and signatures of realistic velocity field and atmosphe ric inhomogeneities associated with convective motions at the surface of cool late-type stars on the emergent profiles of iron spectral lines for a la rge range in stellar parameters. Methods. 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-hydrod ynamical atmosphere models from the Stagger-grid. A set of 35 real unblended, optical Fei and Feii lines of varying excitation potential are considered. Addi tionally, fictitious Fe i and Feii lines (5000 ˚ A 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 il lustrate systematical trends stemming from the intrinsic structural differences among 3D model atmospheres with different stellar parameters. Theoretical line shifts and bisectors are derived to analyz e 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. Results. We find systematic variations in line strength, shift, width , and bisectors, that can be related to the respective physic al conditions at the height of the line formation in the stellar atmospheric environment, in particular the amplitude of the vertical velocity field. Conclusions. 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 star s.

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.

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.

The element abundance ratios of four low-mass stars with extremely low metallicities (abundances of elements heavier than helium) 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 seems surprising, because it had 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 has been 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 about 60 times that of the Sun (and that the supernova left behind a black hole). Taken together with the four previously mentioned low-metallicity stars, we conclude that low-energy supernovae were common in the early Universe, and that such supernovae yielded 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.