Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrm{\mu}$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions...

Red supergiants are cool massive stars and are the largest and the most luminous stars in the universe. They are characterized by irregular or semi-regular photometric variations, the physics of which is not clearly understood. The paper aims at deriving the velocity field in the red supergiant star $\mu$ Cep and relating it to the photometric variability with the help of the tomographic method. The tomographic method allows to recover the line-of-sight velocity distribution over the stellar disk and within different optical-depth slices. The method is applied to a series of high-resolution spectra of $\mu$ Cep, and these results are compared to those obtained from 3D radiative-hydrodynamics CO5BOLD simulations of red supergiants. Fluctuations in the velocity field are compared with photometric and spectroscopic variations, the latter being derived from the TiO band strength and serving (at least partly) a proxy of the variations in effective temperature. The tomographic method reveals a phase shift between the velocity and spectroscopic/photometric variations. This phase shift results in a hysteresis loop in the temperature - velocity plane, with a timescale of a few hundred days, similar to the photometric one. The similarity between the hysteresis loop timescale measured in $\mu$ Cep and the timescale of acoustic waves disturbing the convective pattern suggests that such waves play an important role in triggering the hysteresis loops.

We perform a calibration of the mixing length parameter at the bottom boundary of the convection zone for helium-dominated atmospheres of white dwarfs. This calibration is based on a grid of 3D DB (pure-helium) and DBA (helium-dominated with traces of hydrogen) model atmospheres computed with the CO5BOLD code, and a grid of 1D DB and DBA envelope structures. The 3D models span a parameter space of hydrogen-to-helium abundances between -10.0 and -2.0, surface gravities between 7.5 and 9.0 and effective temperatures between 12000 K and 34000 K. The 1D envelopes cover a similar atmospheric parameter range, but are also calculated with different values of the mixing length parameter, namely ML2/alpha between 0.4 and 1.4. The calibration is performed based on two definitions of the bottom boundary of the convection zone, the Schwarzschild and the zero convective flux boundaries. Thus, our calibration is relevant for applications involving the bulk properties of the convection zone including its total mass, which excludes the spectroscopic technique. Overall, the calibrated ML2/alpha is smaller than what is commonly used in evolutionary models and theoretical determinations of the blue edge of the instability strip for pulsating DB and DBA stars. With calibrated ML2/alpha we are able to deduce more accurate convection zone sizes needed for studies of planetary debris mixing and dredge-up of carbon from the core. We highlight this by calculating examples of metal-rich 3D DBAZ models and finding their convection zone masses. Mixing length calibration represents the first step of in-depth investigations of convective overshoot in white dwarfs with helium-dominated atmospheres.

We present a theoretical description of macroscopic diffusion caused by convective overshoot in pure-hydrogen DA white dwarfs using three-dimensional (3D), closed-bottom, radiation hydrodynamics CO$^5$BOLD simulations. We rely on a new grid of deep 3D white dwarf models in the temperature range 11400 K $\leq T_{\mathrm{eff}} \leq$ 18000 K where tracer particles and a tracer density are used to derive macroscopic diffusion coefficients driven by convective overshoot. These diffusion coefficients are compared to microscopic diffusion coefficients from one-dimensional structures. We find that the mass of the fully mixed region is likely to increase by up to 2.5 orders of magnitude while inferred accretion rates increase by a more moderate order of magnitude. We present evidence that an increase in settling time of up to 2 orders of magnitude is to be expected which is of significance for time-variability studies of polluted white dwarfs. Our grid also provides the most robust constraint on the onset of convective instabilities in DA white dwarfs to be in the effective temperature range from 18000 to 18250 K.

Dust grains forming in the extended atmospheres of AGB stars are critical for the heavy mass loss of these cool luminous giants, as they provide radiative acceleration for the stellar winds. Characteristic mid-IR spectral features indicate that the grains consist mainly of silicates and corundum. The latter species seems to form in a narrow zone within about 2 stellar radii, preceding the condensation of silicate dust, which triggers the outflow. Recent high-angular-resolution observations show clumpy, variable dust clouds at these distances. We explore possible causes for the formation of inhomogeneous dust layers, using 3D dynamical simulations. We modeled the outer convective envelope and the dust-forming atmosphere of an M-type AGB star with the CO5BOLD radiation-hydrodynamics code. The simulations account for frequency-dependent gas opacities, and include a time-dependent description of grain growth and evaporation for corundum and olivine-type silicates. In the inner, gravitationally bound, and corundum-dominated layers of the circumstellar envelope, a patchy distribution of the dust emerges naturally, due to atmospheric shock waves that are generated by large-scale convective flows and pulsations. The formation of silicate dust at somewhat larger distances probably indicates the outer limit of the gravitationally bound layers. The current models do not describe wind acceleration, but the cloud formation mechanism should also work for stars with outflows. Timescales of atmospheric dynamics and grain growth are similar to observed values. In spherical averages of dust densities the variable 3D morphology manifests itself as cycle-to-cycle variations. Grain growth in the wake of large-scale non-spherical shock waves, generated by convection and pulsations, is a likely mechanism for producing the observed clumpy dust clouds, and for explaining their physical and dynamical properties.

A tomographic method, aiming at probing velocity fields at depth in stellar atmospheres, is applied to the red supergiant star {\mu} Cep and to snapshots of 3D radiative-hydrodynamics simulation in order to constrain atmospheric motions and relate them to photometric variability.

We present the first grid of 3D simulations for the pure-helium atmospheres of DB white dwarfs. The simulations were computed with the CO5BOLD radiation-hydrodynamics code and cover effective temperatures and surface gravities between 12000 K < Teff < 34000 K and 7.5 < log(g) (cgs units) < 9.0, respectively. In this introductory work, synthetic spectra calculated from the 3D simulations are compared to appropriate 1D model spectra under a differential approach. This results in the derivation of 3D corrections for the spectroscopically derived atmospheric parameters of DB stars with respect to the 1D ML2/alpha = 1.25 mixing-length parameterisation. No significant Teff corrections are found for the V777 Her instability strip region, and therefore no 3D revision is expected for the empirical blue and red edges of the strip. However, large log(g) corrections are found in the range 12000 K < Teff < 23000 K for all log(g) values covered by the 3D grid. These corrections indicate that 1D model atmospheres overpredict log(g), reminiscent of the results found from 3D simulations of pure-hydrogen white dwarfs. The next step will be to compute 3D simulations with mixed helium and hydrogen atmospheres to comprehend the full implications for the stellar parameters of DB and DBA white dwarfs.

The wind-driving mechanism of asymptotic giant branch (AGB) stars is commonly attributed to a two-step process: first, gas in the stellar atmosphere is levitated by shockwaves caused by stellar pulsation, then accelerated outwards by radiative pressure on newly formed dust, inducing a wind. Dynamical modelling of such winds usually assumes a spherically symmetric star. We explore the potential consequences of complex stellar surface structures, as predicted by three-dimensional (3D) star-in-a-box modelling of M-type AGB stars, on the resulting wind properties with the aim to improve the current wind models. Two different modelling approaches are used; the CO$^5$BOLD 3D star-in-a-box code to simulate the convective, pulsating interior and lower atmosphere of the star, and the DARWIN one-dimensional (1D) code to describe the dynamical atmosphere where the wind is accelerated. The gas dynamics of the inner atmosphere region at distances of $R\sim1-2R_\star$, which both modelling approaches simulate, are compared. Dynamical properties and luminosity variations derived from CO$^5$BOLD interior models are used as input for the inner boundary in DARWIN wind models in order to emulate the effects of giant convection cells and pulsation, and explore their influence on the dynamical properties. The CO$^5$BOLD models are inherently anisotropic, with non-uniform shock fronts and varying luminosity amplitudes, in contrast to the spherically symmetrical DARWIN wind models. DARWIN wind models with CO$^5$BOLD-derived inner boundary conditions produced wind velocities and mass-loss rates comparable to the standard DARWIN models, however the winds show large density variations on time-scales of 10-20 years.

Context. Asymptotic Giant Branch stars are characterised by complex stellar surface dynamics that affect the measurements and amplify the uncertainties on stellar parameters. The resulting motion of the stellar photo-center could have adverse effects on the parallax determination with Gaia. Aims. We explore the impact of the convection-related surface structure in AGBs on the photocentric variability. We quantify these effects to characterise the observed parallax errors and estimate fundamental stellar parameters and dynamical properties. Methods. We use 3D radiative-hydrodynamics simulations of convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the Gaia G band. From those maps, we calculate the intensity-weighted mean of all emitting points tiling the visible stellar surface and evaluate its motion as a function of time. We extract the parallax error from Gaia DR2 for a sample of semiregular variables in the solar neighbourhood and compare it to the synthetic predictions of photo-center displacements. Results. AGB stars show a complex surface morphology characterised by the presence of few large scale long-lived convective cells accompanied by short-lived and small scale structures. As a consequence, the position of the photo-center displays temporal excursions between 0.077 to 0.198 AU (5 to 11% of the corresponding stellar radius), depending on the simulation considered. We show that the convection-related variability accounts for a substantial part to the Gaia DR2 parallax error of our sample of semiregular variables. Finally, we put in evidence for a correlation between the mean photo-center displacement and the stellar fundamental parameters: surface gravity, radius, and pulsation. We denote that parallax variations could be exploited quantitatively using appropriate RHD simulations corresponding to the observed star.

The atmospheres of cool stars are temporally and spatially inhomogeneous due to the effects of convection.The influence of this inhomogeneity, referred to as granulation, on colours has never been investigated over a large range of effective temperatures and gravities. We aim to study, in a quantitative way, the impact of granulation on colours. We use the CIFIST grid of CO5BOLD hydrodynamical models to compute emerging fluxes. These in turn are used to compute theoretical colours in the UBVRI,2MASS,Hipparcos,Gaia and SDSS systems. Every CO5BOLD model has a corresponding one dimensional (1D) plane-parallel LHD model computed for the same atmospheric parameters, which we used to define a "3D correction" that can be applied to colours computed from fluxes computed from any 1D model atmosphere code. The 3D corrections on colours are generally small, of the order of a few hundredths of a magnitude, yet they are far from negligible. We find that ignoring granulation effects can lead to underestimation of Teff by up to 200K and overestimation of gravity by up to 0.5dex, when using colours as diagnostics. We have identified a major shortcoming in how scattering is treated in the current version of the CIFIST grid, which could lead to offsets of the order 0.01mag, especially for colours involving blue and UV bands.We have investigated the Gaia and Hipparcos photometric systems and found that the (G-H_p),(BP-RP) diagram is immune to the effects of granulation. In addition, we point to the potential of the RVS photometry as a metallicity diagnostic. Our investigation shows that the effects of granulation should not be neglected if one wants to use colours as diagnostics of the stellar parameters of F,G,K stars. A limitation is that scattering is treated as true absorption in our current computations, thus our 3D corrections are likely an upper limit to the true effect. (Abridged)

Standard spectroscopic analyses of variable stars are based on hydrostatic one-dimensional model atmospheres. This quasi-static approach has theoretically not been validated. We aim at investigating the validity of the quasi-static approximation for Cepheid variables. We focus on the spectroscopic determination of the effective temperature $T_\mathrm{eff}$, surface gravity $\log \,g$, microturbulent velocity $\xi_\mathrm{t}$, and a generic metal abundance $\log\,A$ -- here taken as iron. We calculate a grid of 1D hydrostatic plane-parallel models covering the ranges in effective temperature and gravity encountered during the evolution of a two-dimensional time-dependent envelope model of a Cepheid computed with the radiation-hydrodynamics code CO5BOLD. We perform 1D spectral syntheses for artificial iron lines in local thermodynamic equilibrium varying the microturbulent velocity and abundance. We fit the resulting equivalent widths to corresponding values obtained from our dynamical model. For the four-parametric case, the stellar parameters are typically underestimated exhibiting a bias in the iron abundance of $\approx-0.2\,\mbox{dex}$. To avoid biases of this kind it is favourable to restrict the spectroscopic analysis to photometric phases $\phi_\mathrm{ph}\approx0.3\ldots 0.65$ using additional information to fix effective temperature and surface gravity. Hydrostatic 1D model atmospheres can provide unbiased estimates of stellar parameters and abundances of Cepheid variables for particular phases of their pulsations. We identified convective inhomogeneities as the main driver behind potential biases. For obtaining a complete view on the effects when determining stellar parameters with 1D models, multi-dimensional Cepheid atmosphere models are necessary for variables of longer period than investigated here.

Cool giant and supergiant star atmospheres are characterized by complex velocity fields originating from convection and pulsation processes which are not fully understood yet. The velocity fields impact the formation of spectral lines, which thus contain information on the dynamics of stellar atmospheres. The tomographic method allows to recover the distribution of the component of the velocity field projected on the line of sight at different optical depths in the stellar atmosphere. The computation of the contribution function to the line depression aims at correctly identifying the depth of formation of spectral lines in order to construct numerical masks probing spectral lines forming at different optical depths. The tomographic method is applied to 1D model atmospheres and to a realistic 3D radiative hydrodynamics simulation performed with CO5BOLD in order to compare their spectral line formation depths and velocity fields. In 1D model atmospheres, each spectral line forms in a restricted range of optical depths. On the other hand, in 3D simulations, the line formation depths are spread in the atmosphere mainly because of temperature and density inhomogeneities. Comparison of CCF profiles obtained from 3D synthetic spectra with velocities from the 3D simulation shows that the tomographic method correctly recovers the distribution of the velocity component projected on the line of sight in the atmosphere.

Standard spectroscopic analyses of Cepheid variables are based on hydrostatic one-dimensional model atmospheres, with convection treated using various formulations of mixing-length theory. This paper aims to carry out an investigation of the validity of the quasi-static approximation in the context of pulsating stars. We check the adequacy of a two-dimensional time-dependent model of a Cepheid-like variable with focus on its spectroscopic properties. With the radiation-hydrodynamics code CO5BOLD, we construct a two-dimensional time-dependent envelope model of a Cepheid with $T_\mathrm{eff}= 5600$ K, $\log g=2.0$, solar metallicity, and a 2.8-day pulsation period. Subsequently, we perform extensive spectral syntheses of a set of artificial iron lines in local thermodynamic equilibrium. The set of lines allows us to systematically study effects of line strength, ionization stage, and excitation potential. We evaluate the microturbulent velocity, line asymmetry, projection factor, and Doppler shifts. The mean Doppler shift is non-zero and negative, -1 km/s, after averaging over several full periods and lines. This residual line-of-sight velocity (related to the "K-term") is primarily caused by horizontal inhomogeneities, and consequently we interpret it as the familiar convective blueshift ubiquitously present in non-pulsating late-type stars. Limited statistics prevent firm conclusions on the line asymmetries. Our two-dimensional model provides a reasonably accurate representation of the spectroscopic properties of a short-period Cepheid-like variable star. Some properties are primarily controlled by convective inhomogeneities rather than by the Cepheid-defining pulsations.

Context: Observations of asymptotic giant branch (AGB) stars with increasing spatial resolution reveal new layers of complexity of atmospheric processes on a variety of scales. Aim: To analyze the physical mechanisms that cause asymmetries and surface structures in observed images, we use detailed 3D dynamical simulations of AGB stars; these simulations self-consistently describe convection and pulsations. Methods: We used the CO5BOLD radiation-hydrodynamics code to produce an exploratory grid of global "star-in-a-box" models of the outer convective envelope and the inner atmosphere of AGB stars to study convection, pulsations, and shock waves and their dependence on stellar and numerical parameters. Results: The model dynamics are governed by the interaction of long-lasting giant convection cells, short-lived surface granules, and strong, radial, fundamental-mode pulsations. Radial pulsations and shorter wavelength, traveling, acoustic waves induce shocks on various scales in the atmosphere. Convection, waves, and shocks all contribute to the dynamical pressure and, thus, to an increase of the stellar radius and to a levitation of material into layers where dust can form. Consequently, the resulting relation of pulsation period and stellar radius is shifted toward larger radii compared to that of non-linear 1D models. The dependence of pulsation period on luminosity agrees well with observed relations. The interaction of the pulsation mode with the non-stationary convective flow causes occasional amplitude changes and phase shifts. The regularity of the pulsations decreases with decreasing gravity as the relative size of convection cells increases. The model stars do not have a well-defined surface. Instead, the light is emitted from a very extended inhomogeneous atmosphere with a complex dynamic pattern of high-contrast features.

Context. Small-scale bright features in the photosphere of the Sun, such as faculae or G-band bright points, appear in connection with small-scale magnetic flux concentrations. Aims. Here we report on a new class of photospheric bright points that are free of magnetic fields. So far, these are visible in numerical simulations only. We explore conditions required for their observational detection. Methods. Numerical radiation (magneto-)hydrodynamic simulations of the near-surface layers of the Sun were carried out. The magnetic field-free simulations show tiny bright points, reminiscent of magnetic bright points, only smaller. A simple toy model for these non-magnetic bright points (nMBPs) was established that serves as a base for the development of an algorithm for their automatic detection. Basic physical properties of 357 detected nMBPs were extracted and statistically evaluated. We produced synthetic intensity maps that mimic observations with various solar telescopes to obtain hints on their detectability. Results. The nMBPs of the simulations show a mean bolometric intensity contrast with respect to their intergranular surroundings of approximately 20%, a size of 60-80 km, and the isosurface of optical depth unity is at their location depressed by 80-100 km. They are caused by swirling downdrafts that provide, by means of the centripetal force, the necessary pressure gradient for the formation of a funnel of reduced mass density that reaches from the subsurface layers into the photosphere. Similar, frequently occurring funnels that do not reach into the photosphere, do not produce bright points. Conclusions. Non-magnetic bright points are the observable manifestation of vertically extending vortices (vortex tubes) in the photosphere. The resolving power of 4-m-class telescopes, such as the DKIST, is needed for an unambiguous detection of them.

Molecular features such as the G-band, CN-band and NH-band are important diagnostics for measuring a star's carbon and nitrogen abundances, especially in metal-poor stars where atomic lines are no longer visible in stellar spectra. Unlike atomic transitions, molecular features tend to form in bands, which cover large wavelength regions in a spectrum. While it is a trivial matter to compute carbon and nitrogen molecular bands under the assumption of 1D, it is extremely time consuming in 3D. In this contribution to the 2016 CO5BOLD workshop we review the improvements made to the 3D spectral synthesis code Linfor3D, and discuss the new challenges found when computing molecular features in 3D.

The $s$-process nucleosynthesis in Asymptotic Giant Branch (AGB) stars depends on the modeling of convective boundaries. We present models and s-process simulations that adopt a treatment of convective boundaries based on the results of hydrodynamic simulations and on the theory of mixing due to gravity waves in the vicinity of convective boundaries. Hydrodynamics simulations suggest the presence of convective boundary mixing (CBM) at the bottom of the thermal pulse-driven convective zone. Similarly, convection-induced mixing processes are proposed for the mixing below the convective envelope during third dredge-up where the 13C pocket for the s process in AGB stars forms. In this work we apply a CBM model motivated by simulations and theory to models with initial mass $M = 2$ and $M = 3M_\odot$, and with initial metal content Z = 0.01 and Z = 0.02. As reported previously, the He-intershell abundance of 12C and 16O are increased by CBM at the bottom of pulse-driven convection zone. This mixing is affecting the $^{22}Ne(\alpha,n)^{25}Mg$ activation and the s-process effciency in the 13C-pocket. In our model CBM at the bottom of the convective envelope during the third dredgeup represents gravity wave mixing. We take further into account that hydrodynamic simulations indicate a declining mixing efficiency already about a pressure scale height from the convective boundaries, compared to mixing-length theory. We obtain the formation of the 13C-pocket with a mass of $\approx 10^{-4}M_\odot$. The final $s$-process abundances are characterized by 0.36 < [s=Fe] < 0.78 and the heavy-to-light s-process ratio is 0.23 < [hs=ls] < 0.45. Finally, we compare our results with stellar observations, pre-solar grain measurements and previous work.

We obtained a total of 20 near-infrared K-band spectro-interferometric snapshot observations of the Mira variables o Cet, R Leo, R Aqr, X Hya, W Vel, and R Cnc with a spectral resolution of about 1500. We compared observed flux and visibility spectra with predictions by CODEX 1D dynamic model atmospheres and with azimuthally averaged intensities based on CO5BOLD 3D dynamic model atmospheres including convection. Our visibility data confirm the presence of spatially extended molecular atmospheres located above the continuum radii with large-scale inhomogeneities or clumps that contribute a few percent of the total flux. The detailed structure of the inhomogeneities or clumps show a variability on time scales of 3 months and above. Both modeling attempts provided satisfactory fits to our data. In particular, they are both consistent with the observed decrease in the visibility function at molecular bands of water vapor and CO, indicating a spatially extended molecular atmosphere. Observational variability phases are mostly consistent with those of the best-fit CODEX models, except for near-maximum phases, where data are better described by near-minimum models. Rosseland angular diameters derived from the model fits are broadly consistent between those based on the 1D and the 3D models and with earlier observations. We derived fundamental parameters including absolute radii, effective temperatures, and luminosities for our sources. Our results provide a first observational support for theoretical results that shocks induced by convection and pulsation in the 3D CO5BOLD models of AGB stars are roughly spherically expanding and of similar nature to those of self-excited pulsations in 1D CODEX models. Unlike for red supergiants, the pulsation- and shock-induced dynamics can levitate the molecular atmospheres of Mira variables to extensions that are consistent with observations.

Evolved stars are among the largest and brightest stars and they are ideal targets for the new generation of sensitive, high resolution instrumentation that provides spectrophotometric, interferometric, astrometric, and imaging observables. The interpretation of the complex stellar surface images requires numerical simulations of stellar convection that take into account multi-dimensional time-dependent radiation hydrodynamics with realistic input physics. We show how the evolved star simulations are obtained using the radiative hydrodynamics code CO5BOLD and how the accurate observables are computed with the post-processing radiative transfer code Optim3D. The synergy between observations and theoretical work is supported by a proper and quantitative analysis using these simulations, and by strong constraints from the observational side.

We present the first radiation magnetohydrodynamics simulations of the atmosphere of white dwarf stars. We demonstrate that convective energy transfer is seriously impeded by magnetic fields when the plasma-beta parameter, the thermal to magnetic pressure ratio, becomes smaller than unity. The critical field strength that inhibits convection in the photosphere of white dwarfs is in the range B = 1-50 kG, which is much smaller than the typical 1-1000 MG field strengths observed in magnetic white dwarfs, implying that these objects have radiative atmospheres. We have then employed evolutionary models to study the cooling process of high-field magnetic white dwarfs, where convection is entirely suppressed during the full evolution (B > 10 MG). We find that the inhibition of convection has no effect on cooling rates until the effective temperature (Teff) reaches a value of around 5500 K. In this regime, the standard convective sequences start to deviate from the ones without convection owing to the convective coupling between the outer layers and the degenerate reservoir of thermal energy. Since no magnetic white dwarfs are currently known at the low temperatures where this coupling significantly changes the evolution, effects of magnetism on cooling rates are not expected to be observed. This result contrasts with a recent suggestion that magnetic white dwarfs with Teff < 10,000 K cool significantly slower than non-magnetic degenerates.

We present an extended grid of mean three-dimensional (3D) spectra for low-mass, pure-hydrogen atmosphere DA white dwarfs (WDs). We use CO5BOLD radiation-hydrodynamics 3D simulations covering Teff = 6000-11,500 K and logg = 5-6.5 (cgs units) to derive analytical functions to convert spectroscopically determined 1D temperatures and surface gravities to 3D atmospheric parameters. Along with the previously published 3D models, the 1D to 3D corrections are now available for essentially all known convective DA WDs (i.e., logg = 5-9). For low-mass WDs, the correction in temperature is relatively small (a few per cent at the most), but the surface gravities measured from the 3D models are lower by as much as 0.35 dex. We revisit the spectroscopic analysis of the extremely low-mass (ELM) WDs, and demonstrate that the 3D models largely resolve the discrepancies seen in the radius and mass measurements for relatively cool ELM WDs in eclipsing double WD and WD + milli-second pulsar binary systems. We also use the 3D corrections to revise the boundaries of the ZZ Ceti instability strip, including the recently found ELM pulsators.

We present the atmospheric structure and the fundamental parameters of three red supergiants, increasing the sample of RSGs observed by near-infrared spectro-interferometry. Additionally, we test possible mechanisms that may explain the large observed atmospheric extensions of RSGs. We carried out spectro-interferometric observations of 3 RSGs in the near-infrared K-band with the VLTI/AMBER instrument at medium spectral resolution. To comprehend the extended atmospheres, we compared our observational results to predictions by available hydrostatic PHOENIX, available 3-D convection, and new 1-D self-excited pulsation models of RSGs. Our near-infrared flux spectra are well reproduced by the PHOENIX model atmospheres. The continuum visibility values are consistent with a limb-darkened disk as predicted by the PHOENIX models, allowing us to determine the angular diameter and the fundamental parameters of our sources. Nonetheless, in the case of V602 Car and HD 95686, the PHOENIX model visibilities do not predict the large observed extensions of molecular layers, most remarkably in the CO bands. Likewise, the 3-D convection models and the 1-D pulsation models with typical parameters of RSGs lead to compact atmospheric structures as well, which are similar to the structure of the hydrostatic PHOENIX models. They can also not explain the observed decreases in the visibilities and thus the large atmospheric molecular extensions. The full sample of our RSGs indicates increasing observed atmospheric extensions with increasing luminosity and decreasing surface gravity, and no correlation with effective temperature or variability amplitude, which supports a scenario of radiative acceleration on Doppler-shifted molecular lines.

A calibration of the mixing-length parameter in the local mixing-length theory (MLT) is presented for the lower part of the convection zone in pure-hydrogen atmosphere white dwarfs. The parameterization is performed from a comparison of 3D CO5BOLD simulations with a grid of 1D envelopes with a varying mixing-length parameter. In many instances, the 3D simulations are restricted to the upper part of the convection zone. The hydrodynamical calculations suggest, in those cases, that the entropy of the upflows does not change significantly from the bottom of the convection zone to regions immediately below the photosphere. We rely on this asymptotic entropy value, characteristic of the deep and adiabatically stratified layers, to calibrate 1D envelopes. The calibration encompasses the convective hydrogen-line (DA) white dwarfs in the effective temperature range 6000 < Teff (K) < 15,000 and the surface gravity range 7.0 < log g < 9.0. It is established that the local MLT is unable to reproduce simultaneously the thermodynamical, flux, and dynamical properties of the 3D simulations. We therefore propose three different parameterizations for these quantities. The resulting calibration can be applied to structure and envelope calculations, in particular for pulsation, chemical diffusion, and convective mixing studies. On the other hand, convection has no effect on the white dwarf cooling rates until there is a convective coupling with the degenerate core below Teff ~ 5000 K. In this regime, the 1D structures are insensitive to the MLT parameterization and converge to the mean 3D results, hence remain fully appropriate for age determinations.

Evolved stars are among the largest and brightest stars and they are ideal targets for the new generation of sensitive, high resolution instrumentation that pro- vides spectrophotometric, interferometric, astrometric, and imaging observables. The interpretation of the complex stellar surface images requires numerical simulations of stellar convection that take into account multi-dimensional time-dependent radiation hydrodynamics with realistic input physics. We show how the evolved star simulations are obtained using the radiative hydrodynamics code CO5BOLD and how the accurate observables are computed with the post-processing radiative transfer code Optim3D. The synergy between observations and theoretical work is supported by a proper and quantitative analysis using these simulations, and by strong constraints from the obser- vational side.

We have identified eight to ten new cool white dwarfs from the Large Area Survey (LAS) Data Release 9 of the United Kingdom InfraRed Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS). The data set was paired with the Sloan Digital Sky Survey (SDSS) to obtain proper motions and a broad ugrizYJHK wavelength coverage. Optical spectroscopic observations were secured at Gemini Observatory and confirm the degenerate status for eight of our targets. The final sample includes two additional white dwarf candidates with no spectroscopic observations. We rely on improved 1D model atmospheres and new multi-dimensional simulations with CO5BOLD to review the stellar parameters of the published LAS white dwarf sample along with our additional discoveries. Most of the new objects possess very cool atmospheres with effective temperatures below 5000 K, including two pure-hydrogen remnants with a cooling age between 8.5 and 9.0 Gyr, and tangential velocities in the range 40 km/s < vtan < 60 km/s. They are likely thick disk 10-11 Gyr-old objects. In addition we find a resolved double degenerate system with vtan ~ 155 km/s and a cooling age between 3.0 and 5.0 Gyr. These white dwarfs could be disk remnants with a very high velocity or former halo G stars. We also compare the LAS sample with earlier studies of very cool degenerates and observe a similar deficit of helium-dominated atmospheres in the range 5000 < Teff (K) < 6000. We review the possible explanations for the spectral evolution from helium-dominated towards hydrogen-rich atmospheres at low temperatures.

We present two epochs of MPG/ESO 2.2m GROND simultaneous 6-band ($r'i'z'JHK$) photometric monitoring of the closest known L/T transition brown dwarf binary WISE J104915.57-531906.1AB. We report here the first resolved variability monitoring of both the T0.5 and L7.5 components. We obtained 4 hours of focused observations on the night of UT 2013-04-22, as well as 4 hours of defocused (unresolved) observations on the night of UT 2013-04-16. We note a number of robust trends in our light curves. The $r'$ and $i'$ light curves appear to be anticorrelated with $z'$ and $H$ for the T0.5 component and in the unresolved lightcurve. In the defocused dataset, $J$ appears correlated with $z'$ and $H$ and anticorrelated with $r'$ and $i'$, while in the focused dataset we measure no variability for $J$ at the level of our photometric precision, likely due to evolving weather phenomena. In our focused T0.5 component lightcurve, the $K$ band lightcurve displays a significant phase offset relative to both $H$ and $z'$. We argue that the measured phase offsets are correlated with atmospheric pressure probed at each band, as estimated from 1D atmospheric models. We also report low-amplitude variability in $i'$ and $z'$ intrinsic to the L7.5 component.

We present the first grid of mean three-dimensional (3D) spectra for pure-hydrogen (DA) white dwarfs based on 3D model atmospheres. We use CO5BOLD radiation-hydrodynamics 3D simulations instead of the mixing-length theory for the treatment of convection. The simulations cover the effective temperature range of 6000 < Teff (K) < 15,000 and the surface gravity range of 7 < log g < 9 where the large majority of DAs with a convective atmosphere are located. We rely on horizontally averaged 3D structures (over constant Rosseland optical depth) to compute <3D> spectra. It is demonstrated that our <3D> spectra can be smoothly connected to their 1D counterparts at higher and lower Teff where the 3D effects are small. Analytical functions are provided in order to convert spectroscopically determined 1D effective temperatures and surface gravities to 3D atmospheric parameters. We apply our improved models to well studied spectroscopic data sets from the Sloan Digital Sky Survey and the White Dwarf Catalog. We confirm that the so-called high-log g problem is not present when employing <3D> spectra and that the issue was caused by inaccuracies in the 1D mixing-length approach. The white dwarfs with a radiative and a convective atmosphere have derived mean masses that are the same within ~0.01 Msun, in much better agreement with our understanding of stellar evolution. Furthermore, the 3D atmospheric parameters are in better agreement with independent Teff and log g values from photometric and parallax measurements.

3D model atmospheres for giants, dwarfs, and white dwarfs, computed with the CO5BOLD code and part of the CIFIST grid, have been used for spectroscopic and asteroseismic studies. Unlike existing plane-parallel 1D structures, these simulations predict the spatially and temporally resolved emergent intensity so that granulation can be analysed, which provides insights on how convective energy transfer operates in stars. The wide range of atmospheric parameters of the CIFIST 3D simulations (3600 < Teff (K) < 13,000 and 1 < log g < 9) allows the comparison of convective processes in significantly different environments. We show that the relative intensity contrast is correlated with both the Mach and Peclet numbers in the photosphere. The horizontal size of granules varies between 3 and 10 times the local pressure scale height, with a tight correlation between the factor and the Mach number of the flow. Given that convective giants, dwarfs, and white dwarfs cover the same range of Mach and Peclet numbers, we conclude that photospheric convection operates in a very similar way in those objects.

In this work we present first results of our current project aimed at combining the 3D hydrodynamical stellar atmosphere approach with non-LTE (NLTE) spectral line synthesis for a number of key chemical species. We carried out a full 3D-NLTE spectrum synthesis of the oxygen IR 777 nm triplet, using a modified and improved version of our NLTE3D package to calculate departure coefficients for the atomic levels of oxygen in a CO5BOLD 3D hydrodynamical solar model atmosphere. Spectral line synthesis was subsequently performed with the Linfor 3D code. In agreement with previous studies, we find that the lines of the oxygen triplet produce deeper cores under NLTE conditions, due to the diminished line source function in the line forming region. This means that the solar oxygen IR 777 nm lines should be stronger in NLTE, leading to negative 3D NLTE-LTE abundance corrections. Qualitatively this result would support previous claims for a relatively low solar oxygen abundance. Finally, we outline several further steps that need to be taken in order to improve the physical realism and numerical accuracy of our current 3D-NLTE calculations.

We review recent advancements in modeling the stellar to substellar transition. The revised molecular opacities, solar oxygen abundances and cloud models allow to reproduce the photometric and spectroscopic properties of this transition to a degree never achieved before, but problems remain in the important M-L transition characteristic of the effective temperature range of characterizable exoplanets. We discuss of the validity of these classical models. We also present new preliminary global Radiation HydroDynamical M dwarfs simulations.

A sequence of pure-hydrogen CO5BOLD 3D model atmospheres of DA white dwarfs is presented for a surface gravity of log g = 8 and effective temperatures from 6000 to 13,000 K. We show that convective properties, such as flow velocities, characteristic granulation size and intensity contrast of the granulation patterns, change significantly over this range. We demonstrate that these 3D simulations are not sensitive to numerical parameters unlike the 1D structures that considerably depend on the mixing-length parameters. We conclude that 3D spectra can be used directly in the spectroscopic analyses of DA white dwarfs. We confirm the result of an earlier preliminary study that 3D model spectra provide a much better characterization of the mass distribution of white dwarfs and that shortcomings of the 1D mixing-length theory are responsible for the spurious high-log g determinations of cool white dwarfs. In particular, the 1D theory is unable to account for the cooling effect of the convective overshoot in the upper atmospheres.

To derive space velocities of stars along the line of sight from wavelength shifts in stellar spectra requires accounting for a number of second-order effects. For most stars, gravitational redshifts, convective blueshifts, and transverse stellar motion are the dominant contributors. We provide theoretical corrections for the net velocity shifts due to convection expected for the measurements from the Gaia Radial Velocity Spectrometer (RVS). We used a set of three-dimensional time-dependent simulations of stellar surface convection computed with CO5BOLD to calculate spectra of late-type stars in the Gaia RVS range and to infer the net velocity offset that convective motions will induce in radial velocities derived by cross-correlation. The net velocity shifts derived by cross-correlation depend both on the wavelength range and spectral resolution of the observations. Convective shifts for Gaia RVS observations are less than 0.1 km/s for late-K-type stars, and they increase with stellar mass, reaching about 0.3 km/s or more for early F-type dwarfs. This tendency is the result of an increase with effective temperature in both temperature and velocity fluctuations in the line-forming region. Our simulations also indicate that the net RVS convective shifts can be positive (i.e. redshifts) in some cases. Overall, the blueshifts weaken slightly with increasing surface gravity, and are enhanced at low metallicity. Gravitational redshifts amount up to 0.7 km/s and dominate convective blueshifts for dwarfs, but become much weaker for giants.

The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: 1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; 2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image reconstruction process was "blind" in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with 6 to 10 telescopes, should be well suited to provide images of complex sceneries.

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.

Within the next few years, GAIA and several instruments aiming at imag- ing extrasolar planets will see first light. In parallel, low mass planets are being searched around red dwarfs which offer more favourable conditions, both for radial velocity de- tection and transit studies, than solar-type stars. Authors of the model atmosphere code which has allowed the detection of water vapour in the atmosphere of Hot Jupiters re- view recent advancement in modelling the stellar to substellar transition. The revised solar oxygen abundances and cloud model allow for the first time to reproduce the pho- tometric and spectroscopic properties of this transition. Also presented are highlight results of a model atmosphere grid for stars, brown dwarfs and extrasolar planets.

High-resolution images of the solar surface show a granulation pattern of hot rising and cooler downward-sinking material -- the top of the deep-reaching solar convection zone. Convection plays a role for the thermal structure of the solar interior and the dynamo acting there, for the stratification of the photosphere, where most of the visible light is emitted, as well as for the energy budget of the spectacular processes in the chromosphere and corona. Convective stellar atmospheres can be modeled by numerically solving the coupled equations of (magneto)hydrodynamics and non-local radiation transport in the presence of a gravity field. The CO5BOLD code described in this article is designed for so-called "realistic" simulations that take into account the detailed microphysics under the conditions in solar or stellar surface layers (equation-of-state and optical properties of the matter). These simulations indeed deserve the label "realistic" because they reproduce the various observables very well -- with only minor differences between different implementations. The agreement with observations has improved over time and the simulations are now well-established and have been performed for a number of stars. Still, severe challenges are encountered when it comes to extending these simulations to include ideally the entire star or substellar object: the strong stratification leads to completely different conditions in the interior, the photosphere, and the corona. Simulations have to account for a large range of spatial and time scales and also non-equilibrium processes. Last but not least, realistic simulations are based on detailed microphysics and depend on the quality of the input data. This article provides an overview of the physical problem and the numerical solution and the capabilities of CO5BOLD, illustrated with a number of applications.

Red supergiants are massive evolved stars that contribute extensively to the chemical enrichment of our Galaxy. It has been shown that convection in those stars gives rise to large granules that cause surface inhomogeneities and shock waves in the photosphere. The understanding of their dynamics is crucial to unveil the unknown mass-loss mechanism, their chemical composition and stellar parameters. We present a new generation of red supergiants simulations with a more sophisticated opacity treatment done with 3D radiative- hydrodynamics CO5BOLD. In the code, the coupled equations of compressible hydrodynamics and non-local radiation transport are solved in the presence of a spherical potential. The stellar core is replaced by a special spherical inner boundary condition, where the gravitational potential is smoothed and the energy production by fusion is mimicked by a simply producing heat corresponding to the stellar luminosity. The post-processing radiative transfer code OPTIM3D is used to extract spectroscopic and interferometric observables. We show that the relaxation of the assumption of frequency-independent opacities shows a steeper mean thermal gradient in the optical thin region that affect strongly the atomic strengths and the spectral energy distribution. Moreover, the weaker temperature fluctuations reduce the incertitude on the radius determination with interferometry. We show that 1D models of red supergiants must include a turbulent velocity calibrated on 3D simulations to obtain the effective surface gravity that mimic the effect of turbulent pressure on the stellar atmosphere. We provide an empirical calibration of the ad-hoc micro- and macroturbulence parameters for 1D models using the 3D simulations: we find that there is not a clear distinction between the different macroturbulent profiles needed in 1D models to fit 3D synthetic lines.

The surface gravities of cool (Teff < 13,000 K) hydrogen-atmosphere DA white dwarfs, determined from spectroscopic analyses, are found to be significantly higher than the canonical value of log g ~ 8 expected for these stars. It was recently concluded that a problem with the treatment of convective energy transport within the framework of the mixing-length theory was the most plausible explanation for this high-log g problem. We pursue the investigation of this discrepancy by computing model spectra of cool convective white dwarfs from a small sequence (11,300 K < Teff < 12,800 K) of 3D hydrodynamical model atmospheres, which feature a sophisticated treatment of convection and radiative transfer. Our approach is to proceed with a differential analysis between 3D and standard 1D models. We find that the 3D spectra predict significantly lower surface gravities, with corrections of the right amplitude as a function of effective temperature to obtain values of log g ~ 8 on average. We conclude that the surface gravity distribution of cool convective DA white dwarfs is much closer to that of hotter radiative objects when using, for the treatment of the convection, 3D models instead of the mixing-length framework.

Context. It has been shown that convection in red supergiant stars gives rise to large granules causing surface inhomogeneities together with shock waves in the photosphere. The resulting motion of the photocenter (on time scales ranging from months to years) could possibly have adverse effects on the parallax determination with Gaia. Aims. We explore the impact of the granulation on the photocentric and photometric variability. We quantify these effects in order to better characterize the error possibly altering the parallax. Methods. We use 3D radiative-hydrodynamics simulations of convection with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps and spectra in the Gaia G band [325-1030 nm]. Results. We provide astrometric and photometric predictions from 3D simulations of RSGs that are used to evaluate the degradation of the astrometric parameters of evolved stars derived by Gaia. We show from RHD simulations that a supergiant like Betelgeuse exhibits a photocentric noise characterised by a standard deviation of the order of 0.1 AU. The number of bright giant and supergiant stars whose Gaia parallaxes will be altered by the photocentric noise ranges from a few tens to several thousandths. The degradation of the astrometric fit due to the presence of this photocentric noise will be noticeable up to about 5 kpc for the brightest supergiants. Moreover, parallaxes of supergiants are affected by a error of the order of a few percents. We show that the photocentric noise, as predicted by the 3D simulation, does account for a substantial part of the supplementary 'cosmic noise' that affects Hipparcos measurements of Betelgeuse and Antares.

Since the discovery of brown dwarfs in 1994, and the discovery of dust cloud formation in the latest Very Low Mass Stars (VLMs) and Brown Dwarfs (BDs) in 1996, the most important challenge in modeling their atmospheres as become the understanding of cloud formation and advective mixing. For this purpose, we have developed radiation hydrodynamic 2D model atmosphere simulations to study the formation of forsterite dust in presence of advection, condensation, and sedimentation across the M-L-T VLMs to BDs sequence (Teff = 2800 K to 900 K, Freytag et al. 2010). We discovered the formation of gravity waves as a driving mechanism for the formation of clouds in these atmospheres, and derived a rule for the velocity field versus atmospheric depth and Teff , which is relatively insensitive to gravity. This rule has been used in the construction of the new model atmosphere grid, BT-Settl, to determine the microturbulence velocity, the diffusion coefficient, and the advective mixing of molecules as a function of depth. This new model grid of atmospheres and synthetic spectra has beencomputedfor100,000K>Teff >400K,5.5>logg>-0.5,and[M/H]=+0.5to -1.5, and the reference solar abundances of Asplund et al. (2009). We found that the new solar abundances allow an improved (close to perfect) reproduction of the photo- metric and spectroscopic VLMs properties, and, for the first time, a smooth transition between stellar and substellar regimes -- unlike the transition between the NextGen models from Hauschildt et al. 1999a,b, and the AMES-Dusty models from Allard et al. 2001). In the BDs regime, the BT-Settl models propose an improved explanation for the M-L-T spectral transition. In this paper, we therefore present the new BT-Settl model atmosphere grid, which explains the entire transition from the stellar to planetary mass regimes.

We present VLT-UVES Li abundances for 28 halo dwarf stars between [Fe/H]=-2.5 and -3.5, 10 of which have [Fe/H]<-3. Four different T_eff scales have been used. Direct Infrared Flux Method (IRFM) has been used on the basis of 2MASS infrared photometry. H_alpha wings have been fitted against synthetic grids computed by means of 1D LTE atmosphere models, assuming different self-broadening theories. Finally, a grid of H_alpha profiles has been computed by means of 3D hydrodynamical atmosphere models. The Li I doublet at 670.8 nm has been used to measure A(Li) by means of 3D hydrodynamical NLTE spectral syntheses. An analytical fit of A(Li)(3D, NLTE) as a function of equivalent width, T_eff, log g, and [Fe/H] has been derived and is made available. A(Li) does not exhibit a plateau below [Fe/H]=-3. A strong positive correlation with [Fe/H] appears, not influenced by the choice of the T_eff estimator. From a linear fit, we obtain a strong slope of about 0.30 dex in A(Li) per dex in [Fe/H], significant to 2-3 sigma, and consistent among all the four T_eff estimators. A significant slope is also detected in the A(Li)--T_eff plane, driven mainly by the coolest stars in the sample which appear Li-poor. Removing such stars does not alter the behavior in the A(Li)-[Fe/H] plane. The scatter in A(Li) increases by a factor of 2 towards lower metallicities, while the plateau appears very thin above [Fe/H]=-2.8. The meltdown of the Spite plateau below [Fe/H]\sim-3 is established, but its cause is unclear. If the primordial A(Li) is the one derived from WMAP, it appears difficult to envision a single depletion phenomenon producing a thin, metallicity independent plateau above [Fe/H]=-2.8, and a highly scattered, metallicity dependent distribution below. The fact that no star below [Fe/H]=-3 lies above the plateau suggests that they formed at plateau level and underwent subsequent depletion.

The measurement of oxygen lines in metal-poor unevolved stars, in particular near-UV OH lines, can provide invaluable information on the properties of the Early Galaxy. Near-UV OH lines constitute an important tool to derive oxygen abundances in metal-poor dwarf stars. Therefore, it is important to correctly model the line formation of OH lines, especially in metal-poor stars, where 3D hydrodynamical models commonly predict cooler temperatures than plane-parallel hydrostatic models in the upper photosphere. We have made use of a grid of 52 3D hydrodynamical model atmospheres for dwarf stars computed using the code CO5BOLD, extracted from the more extended CIFIST grid. The 52 models cover the effective temperature range 5000-6500K, the surface gravity range 3.5-4.5 and the metallicity range -3<[Fe/H]<0. We determine 3D-LTE abundance corrections in all the 52 3D models for several OH lines and FeI lines of different excitation potentials. These 3D-LTE corrections are generally negative reaching values of roughly -1 dex (for the OH 3167 with excitation potential of approximately 1 eV) for the higher temperatures and surface gravities. We apply these 3D-LTE corrections to the individual O abundances, derived from OH lines, of a sample the metal-poor dwarf stars reported in Israelian et al.(1998, 2001) and Boesgaard et al.(1999), by interpolating the stellar parameters of the dwarfs in the grid of 3D-LTE corrections. The new 3D-LTE [O/Fe] ratio still keeps a similar trend as the 1D-LTE, i.e, increasing towards lower [Fe/H] values. We applied 1D-NLTE corrections to 3D FeI abundances and we still see an increasing [O/Fe] ratio towards lower metallicites. However, the Galactic [O/Fe] ratio must be revisited once 3D-NLTE corrections become available for OH and Fe lines for a grid of 3D hydrodynamical model atmospheres.

We report on the development of an implicit multi-D hydrodynamic code for stellar evolution. We present two test-cases relevant for the first scientific goal of the code: the simulation of convection in pulsating stars. First results on a realistic stellar model are also presented.

Context. The red supergiant (RSG) Betelgeuse is an irregular variable star. Convection may play an important role in understanding this variability. Interferometric observations can be interpreted using sophisticated simulations of stellar convection. Aims. We compare the visibility curves and closure phases obtained from our 3D simulation of RSG convection with CO5BOLD to various interferometric observations of Betelgeuse from the optical to the H band in order to characterize and measure the convection pattern on this star. Methods. We use 3D radiative-hydrodynamics (RHD) simulation to compute intensity maps in different filters and we thus derive interferometric observables using the post-processing radiative transfer code OPTIM3D. The synthetic visibility curves and closure phases are compared to observations. Results. We provide a robust detection of the granulation pattern on the surface of Betelgeuse in the optical and in the H band based on excellent fits to the observed visibility points and closure phases. Moreover, we determine that the Betelgeuse surface in the H band is covered by small to medium scale (5-15 mas) convection-related surface structures and a large (30 mas) convective cell. In this spectral region, H2O molecules are the main absorbers and contribute to the small structures and to the position of the first null of the visibility curve (i.e. the apparent stellar radius).

In the last decade, the photospheric solar metallicity as determined from spectroscopy experienced a remarkable downward revision. Part of this effect can be attributed to an improvement of atomic data and the inclusion of NLTE computations, but also the use of hydrodynamical model atmospheres seemed to play a role. This "decrease" with time of the metallicity of the solar photosphere increased the disagreement with the results from helioseismology. With a CO5BOLD 3D model of the solar atmosphere, the CIFIST team at the Paris Observatory re-determined the photospheric solar abundances of several elements, among them C, N, and O. The spectroscopic abundances are obtained by fitting the equivalent width and/or the profile of observed spectral lines with synthetic spectra computed from the 3D model atmosphere. We conclude that the effects of granular fluctuations depend on the characteristics of the individual lines, but are found to be relevant only in a few particular cases. 3D effects are not reponsible for the systematic lowering of the solar abundances in recent years. The solar metallicity resulting from this analysis is Z=0.0153, Z/X=0.0209.

Observationally, spectra of brown dwarfs indicate the presence of dust in their atmospheres while theoretically it is not clear what prevents the dust from settling and disappearing from the regions of spectrum formation. Consequently, standard models have to rely on ad hoc assumptions about the mechanism that keeps dust grains aloft in the atmosphere. We apply hydrodynamical simulations to develop an improved physical understanding of the mixing properties of macroscopic flows in M dwarf and brown dwarf atmospheres, in particular of the influence of the underlying convection zone. We performed 2D radiation hydrodynamics simulations including a description of dust grain formation and transport with the CO5BOLD code. The simulations cover the very top of the convection zone and the photosphere including the dust layers for effective temperatures between 900K and 2800K, all with logg=5 assuming solar chemical composition. Convective overshoot occurs in the form of exponentially declining velocities with small scale heights, so that it affects only the region immediately above the almost adiabatic convective layers. From there on, mixing is provided by gravity waves that are strong enough to maintain thin dust clouds in the hotter models. With decreasing effective temperature, the amplitudes of the waves become smaller but the clouds become thicker and develop internal convective flows that are more efficient in mixing material than gravity waves. The presence of clouds leads to a highly structured appearance of the stellar surface on short temporal and small spatial scales. We identify convectively excited gravity waves as an essential mixing process in M dwarf and brown dwarf atmospheres. Under conditions of strong cloud formation, dust convection is the dominant self-sustaining mixing component.

The use of hydrodynamical simulations, the selection of atomic data, and the computation of deviations from local thermodynamical equilibrium for the analysis of the solar spectra have implied a downward revision of the solar metallicity. We are in the process of using the latest simulations computed with the CO5BOLD code to reassess the solar chemical composition. We determine the solar photospheric carbon abundance by using a radiation-hydrodynamical CO5BOLD model, and compute the departures from local thermodynamical equilibrium by using the Kiel code. We measure equivalent widths of atomic CI lines on high resolution, high signal-to-noise ratio solar atlases. Deviations from local thermodynamic equilibrium are computed in 1D with the Kiel code. Our recommended value for the solar carbon abundance, relies on 98 independent measurements of observed lines and is A(C)=8.50+-0.06, the quoted error is the sum of statistical and systematic error. Combined with our recent results for the solar oxygen and nitrogen abundances this implies a solar metallicity of Z=0.0154 and Z/X=0.0211. Our analysis implies a solar carbon abundance which is about 0.1 dex higher than what was found in previous analysis based on different 3D hydrodynamical computations. The difference is partly driven by our equivalent width measurements (we measure, on average, larger equivalent widths with respect to the other work based on a 3D model), in part it is likely due to the different properties of the hydrodynamical simulations and the spectrum synthesis code. The solar metallicity we obtain from the CO5BOLD analyses is in slightly better agreement with the constraints of helioseismology than the previous 3D abundance results. (Abridged)

Aims. We aim to explore the photosphere of the very cool late-type star VX Sgr and in particular the existence and characterization of molecular layers above the continuum forming photosphere. Methods. We obtained interferometric observations with the VLTI/AMBER interferometer using the fringe tracker FINITO in the spectral domain 1.45-2.50 micron with a spectral resolution of about 35 and baselines ranging from 15 to 88 meters.We perform independent image reconstruction for different wavelength bins and fit the interferometric data with a geometrical toy model.We also compare the data to 1D dynamical models of Miras atmosphere and to 3D hydrodynamical simulations of red supergiant (RSG) and asymptotic giant branch (AGB) stars. Results. Reconstructed images and visibilities show a strong wavelength dependence. The H-band images display two bright spots whose positions are confirmed by the geometrical toy model. The inhomogeneities are qualitatively predicted by 3D simulations. At about 2,00 micron and in the region 2,35 - 2,50 micron, the photosphere appears extended and the radius is larger than in the H band. In this spectral region, the geometrical toy model locates a third bright spot outside the photosphere that can be a feature of the molecular layers. The wavelength dependence of the visibility can be qualitatively explained by 1D dynamical models of Mira atmospheres. The best-fitting photospheric models show a good match with the observed visibilities and give a photospheric diameter of theta = 8,82+-0,50 mas. The H2O molecule seems to be the dominant absorber in the molecular layers. Conclusions. We show that the atmosphere of VX Sgr rather resembles Mira/AGB star model atmospheres than RSG model atmospheres. In particular, we see molecular (water) layers that are typical for Mira stars.

We present solar photospheric abundances for 12 elements from optical and near-infrared spectroscopy. The abundance analysis was conducted employing 3D hydrodynamical (CO5BOLD) as well as standard 1D hydrostatic model atmospheres. We compare our results to others with emphasis on discrepancies and still lingering problems, in particular exemplified by the pivotal abundance of oxygen. We argue that the thermal structure of the lower solar photosphere is very well represented by our 3D model. We obtain an excellent match of the observed center-to-limb variation of the line-blanketed continuum intensity, also at wavelengths shortward of the Balmer jump.

Grids of stellar atmosphere models and associated synthetic spectra are numerical products which have a large impact in astronomy due to their ubiquitous application in the interpretation of radiation from individual stars and stellar populations. 3D model atmospheres are now on the verge of becoming generally available for a wide range of stellar atmospheric parameters. We report on efforts to develop a grid of 3D model atmospheres for late-type stars within the CIFIST Team at Paris Observatory. The substantial demands in computational and human labor for the model production and post-processing render this apparently mundane task a challenging logistic exercise. At the moment the CIFIST grid comprises 77 3D model atmospheres with emphasis on dwarfs of solar and sub-solar metallicities. While the model production is still ongoing, first applications are already worked upon by the CIFIST Team and collaborators.

Context. It has been suggested that convection in Red Supergiant (RSG) stars gives rise to large-scale granules causing observable surface inhomogeneities. This convection is also extremely vigorous, and suspected to be one of the causes of mass-loss in RSGs. It must thus be understood in details. Evidence has been accumulated that there are asymmetries in the photospheres of RSGs, but detailedstudies of granulation are still lacking. Interferometric observations offer an exciting possibility to tackle this question, but they are still often interpreted using smooth symmetrical limb-darkened intensity distributions, or very simple spotted ad hoc models. Aims. We explore the impact of the granulation on visibility curves and closure phases using the radiative transfer code OPTIM3D. We simultaneously assess how 3D simulations of convection in RSG with CO5BOLD can be tested against these observations. Methods. We use 3D radiative-hydrodynamics (RHD) simulations of convection to compute intensity maps at various wavelengths and time, from which we derive interferometric visibility amplitudes and phases. We study their behaviour with time, position angle, and wavelength, and compare them to observations of the RSG alpha Ori Results. We provide average limb-darkening coefficients for RSGs. We detail the prospects for the detection and characterization of granulation (contrast, size) on RSGs. We demonstrate that our RHD simulations provide an excellent fit to existing interferometric observation of alpha Ori, contrary to limb darkened disks. This confirms the existence of large convective cells on the surface of Betelgeuse.

The solar oxygen abundance has undergone a major downward revision in the last decade, the most noticeable one being the update including 3D hydrodynamical simulations to model the solar photosphere. Up to now, such an analysis has been carried out only by one group using one radiation-hydrodynamics code. We investigate the photospheric oxygen abundance considering lines from atomic transitions. We also consider the relationship between the solar model used and the resulting solar oxygen abundance, to understand whether the downward abundance revision is specifically related to 3D hydrodynamical effects. We perform a new determination of the solar photospheric oxygen abundance by analysing different high-resolution high signal-to-noise ratio atlases of the solar flux and disc-centre intensity making use of the latest generation of CO5BOLD 3D solar model atmospheres. We find 8.73 < logNoxygen/Nhydrogen +12 < 8.79. The lower and upper value represent extreme assumptions on the role of collisional excitation and ionisation by neutral hydrogen for the NLTE level populations of neutral oxygen. The error of our analysis is +-0.04 +- 0.03 dex, the last being related to NLTE corrections, the first error to any other effect. 3D ``granulation effects'' do not play a decisive role in lowering the oxygen abundance. Our recommended value, considering our present ignorance of the role of collisions with hydrogen atoms on the NLTE level populations of oxygen, is log(Noxygen/Nhydrogen) = 8.76 +- 0.07. The reasons which have led to lower O abundances in the past are identified as (1) the lower equivalent widths adopted, and (2) the choice of neglecting collisions with hydrogen atoms in the statistical equilibrium calculations for oxygen.

Context. Europium is an almost pure r-process element, which may be useful as a reference in nucleocosmochronology. Aims. To determine the photospheric solar abundance using CO5BOLD 3-D hydrodynamical model atmospheres. Methods. Disc-centre and integrated-flux observed solar spectra are used. The europium abundance is derived from the equivalent width measurements. As a reference 1D model atmospheres have been used, in addition. Results. The europium photospheric solar abundance is 0.52 +- 0.02 in agreement with previous determinations. We also determine the photospheric isotopic fraction of Eu(151) to be 49 % +- 2.3 % from the intensity spectra and 50% +-2.3 from the flux spectra. This compares well to the the meteoritic isotopic fraction 47.8%. We explore the 3D corrections also for dwarfs and sub-giants in the temperature range ~5000 K to ~6500 K and solar and 1/10--solar metallicities and find them to be negligible for all the models investigated. Conclusions. Our photospheric Eu abundance is in good agreement with previous determinations based on 1D models. This is in line with our conclusion that 3D effects for this element are negligible in the case of the Sun.

We have written a 3D radiative transfer code that computes emerging spectra and intensity maps. We derive from radiative hydrodynamic (RHD) simulations of RSG stars carried out with CO5BOLD (Freytag et al. 2002) observables expected for red supergiant stars (RSG) especially for interferometric observations, with emphasis on small scale structures. We show that the convection-related surface structures are detectable in the H band with today's interferometers and that the diameter measurement should not be too dependent on the adopted model. The simulations are a great improvement over parametric models for the interpretation of interferometric observations.

We review the current state of modeling convective mixing in AGB stars. The focus is on results obtained through multi-dimensional hydrodynamic simulations of AGB convection, both in the envelope and the unstable He-shell. Using two different codes and a wide range of resolutions and modeling assumptions we find that mixing across convective boundaries is significant for He-shell flash convection. We present a preliminary quantitative analysis of this convectively induced extra mixing, based on a sub-set of our simulations. Other non-standard mixing will be discussed briefly.

The continuum intensity at wavelengths around 1 mm provides an excellent way to probe the solar chromosphere. Future high-resolution millimetre arrays, such as the Atacama Large Millimeter Array (ALMA), will thus produce valuable input for the ongoing controversy on the thermal structure and the dynamics of this layer. Synthetic brightness temperature maps are calculated on basis of three-dimensional radiation (magneto-)hydrodynamic (MHD) simulations. While the millimetre continuum at 0.3mm originates mainly from the upper photosphere, the longer wavelengths considered here map the low and middle chromosphere. The effective formation height increases generally with wavelength and also from disk-centre towards the solar limb. The average intensity contribution functions are usually rather broad and in some cases they are even double-peaked as there are contributions from hot shock waves and cool post-shock regions in the model chromosphere. Taking into account the deviations from ionisation equilibrium for hydrogen gives a less strong variation of the electron density and with it of the optical depth. The result is a narrower formation height range. The average brightness temperature increases with wavelength and towards the limb. The relative contrast depends on wavelength in the same way as the average intensity but decreases towards the limb. The dependence of the brightness temperature distribution on wavelength and disk-position can be explained with the differences in formation height and the variation of temperature fluctuations with height in the model atmospheres.

With numerical experiments we explore the feasibility of using high frequency waves for probing the magnetic fields in the photosphere and the chromosphere of the Sun. We track a plane-parallel, monochromatic wave that propagates through a non-stationary, realistic atmosphere, from the convection-zone through the photosphere into the magnetically dominated chromosphere, where it gets refracted and reflected. We compare the wave travel time between two fixed geometrical height levels in the atmosphere (representing the formation height of two spectral lines) with the topography of the surface of equal magnetic and thermal energy density (the magnetic canopy or beta=1 contour) and find good correspondence between the two. We conclude that high frequency waves indeed bear information on the topography of the `magnetic canopy'.

Broadening and asymmetry of spectral lines in slowly rotating late A-type stars provide evidence for high-amplitude convective motions. The properties of turbulence observed in the A-star atmospheres are not understood theoretically and contradict results of previous numerical simulations of convection. Here we describe an ongoing effort to understand the puzzling convection signatures of A stars with the help of 3-D hydrodynamic simulations. Our approach combines realistic spectrum synthesis and non-grey hydrodynamic models computed with the CO5BOLD code. We discuss these theoretical predictions and confront them with high-resolution spectra of A stars. Our models have, for the first time, succeeded in reproducing the observed profiles of weak spectral lines without introducing fudge broadening parameters.

The evolution of central stars of planetary nebulae can proceed in several distinct ways, either leading to H-deficiency or to H-normal surface composition. Several new simulations of the evolution channels that lead to H-deficiency are now available, mainly the born-again scenarios that are triggered by a He-shell flash during the hot pre-white dwarf evolution phase. A realistic AGB progenitor evolution is important for correct HRD tracks, that allow mass determinations. New hydrodynamic simulations of He-shell flash convection including cases with H-ingestion are now performed, and allow a determination of the convective extra-mixing efficiency. This has direct consequences for the intershell abundance distribution of AGB stars that can be observed in the H-deficient CSPN.

We present the first hydrodynamic, multi-dimensional simulations of He-shell flash convection. Specifically, we investigate the properties of shell convection at a time immediately before the He- luminosity peak during the 15th thermal pulse of a stellar evolution track with initially two solar masses and metallicity Z=0.01. This choice is a representative example of a low-mass asymptotic giant branch thermal pulse. We construct the initial vertical stratification with a set of polytropes to resemble the stellar evolution structure. Convection is driven by a constant volume heating in a thin layer at the bottom of the unstable layer. We calculate a grid of 2D simulations with different resolutions and heating rates. Our set of simulations includes one low-resolution 3D run. The computational domain includes 11.4 pressure scale heights. He-shell flash convection is dominated by large convective cells that are centered in the lower half of the convection zone. Convective rolls have an almost circular appearance because focusing mechanisms exist in the form of the density stratification for downdrafts and the heating of localized eddies that generate upflows. Nevertheless, downdrafts appear to be somewhat more focused. The He-shell flash convection generates a rich spectrum of gravity waves in both stable layers above and beneath the convective shell. The magnitude of the convective velocities from our 1D mixing-length theory model and the rms-averaged vertical velocities from the hydrodynamic model are consistent within a factor of a few. However, the velocity profile in the hydrodynamic simulation is more asymmetric, and decays exponentially inside the convection zone. [abbreviated]

Recent three-dimensional radiation hydrodynamic simulations by Wedemeyer et al. (2004) suggest that the solar chromosphere is highly structured in space and time on scales of only 1000 km and 20-25 sec, resp.. The resulting pattern consists of a network of hot gas and enclosed cool regions which are due to the propagation and interaction of shock fronts. In contrast to many other diagnostics, the radio continuum at millimeter wavelengths is formed in LTE, and provides a rather direct measure of the thermal structure. It thus facilitates the comparison between numerical model and observation. While the involved time and length scales are not accessible with todays equipment for that wavelength range, the next generation of instruments, such as the Atacama Large Millimeter Array (ALMA), will provide a big step towards the required resolution. Here we present results of radiative transfer calculations at mm and sub-mm wavelengths with emphasis on spatial and temporal resolution which are crucial for the ongoing discussion about the chromospheric temperature structure.

We present first realistic numerical simulations of 3D radiative convection in the surface layers of main sequence A-type stars with Teff = 8000 K and 8500 K, log g = 4.4 and 4.0, recently performed with the CO5BOLD radiation hydrodynamics code. The resulting models are used to investigate the structure of the H+HeI and the HeII convection zones in comparison with the predictions of local and non-local convection theories, and to determine the amount of "overshoot" into the stable layers below the HeII convection zone. The simulations also predict how the topology of the photospheric granulation pattern changes from solar to A-type star convection. The influence of the photospheric temperature fluctuations and velocity fields on the shape of spectral lines is demonstrated by computing synthetic line profiles and line bisectors for some representative examples, allowing us to confront the 3D model results with observations.

The radiation hydrodynamic code CO5BOLD has been supplemented with the time-dependent treatment of chemical reaction networks. Advection of particle densities due to the hydrodynamic flow field is also included. The radiative transfer is treated frequency-independently, i.e. grey, so far. The upgraded code has been applied to two-dimensional simulations of carbon monoxide (CO) in the non-magnetic solar photosphere and low chromosphere. For this purpose a reaction network has been constructed, taking into account the reactions which are most important for the formation and dissociation of CO under the physical conditions of the solar atmosphere. The network has been strongly reduced to 27 reactions, involving the chemical species H, H2, C, O, CO, CH, OH, and a representative metal. The resulting CO number density is highest in the cool regions of the reversed granulation pattern at mid-photospheric heights and decreases strongly above. There, the CO abundance stays close to a value of 8.3 on the usual logarithmic abundance scale with [H]=12 but is reduced in hot shock waves which are a ubiquitous phenomenon of the model atmosphere. For comparison, the corresponding equilibrium densities have been calculated, based on the reaction network but also under assumption of instantaneous chemical equilibrium by applying the Rybicki & Hummer (RH) code by Uitenbroek (2001). Owing to the short chemical timescales, the assumption holds for a large fraction of the atmosphere, in particular the photosphere. In contrast, the CO number density deviates strongly from the corresponding equilibrium value in the vicinity of chromospheric shock waves. Simulations with altered reaction network clearly show that the formation channel via hydroxide (OH) is the most important one under the conditions of the solar atmosphere.

Cool giant stars are highly dynamical objects, and complex micro-physical processes play an important role in their extended atmospheres and winds. The interpretation of observations, and in particular of high-resolution IR spectra, requires realistic self-consistent model atmospheres. Current dynamical models include rather detailed micro-physics, and the resulting synthetic spectra compare reasonably well with observations. A transition from qualitative to quantitative modelling is taking place at present. We give an overview of existing dynamical model atmospheres for AGB stars, discussing recent advances and current trends in modelling, including 3D 'star-in-a-box' models. When comparing synthetic spectra and other observable properties resulting from dynamical models with observations we focus on the near- and mid-IR wavelength range.

Three-dimensional numerical simulations with CO5BOLD, a new radiation hydrodynamics code, result in a dynamic, thermally bifurcated model of the non-magnetic chromosphere of the quiet Sun. The 3-D model includes the middle and low chromosphere, the photosphere, and the top of the convection zone, where acoustic waves are excited by convective motions. While the waves propagate upwards, they steepen into shocks, dissipate, and deposit their mechanical energy as heat in the chromosphere. Our numerical simulations show for the first time a complex 3-D structure of the chromospheric layers, formed by the interaction of shock waves. Horizontal temperature cross-sections of the model chromosphere exhibit a network of hot filaments and enclosed cool regions. The horizontal pattern evolves on short time-scales of the order of typically 20 - 25 seconds, and has spatial scales comparable to those of the underlying granulation. The resulting thermal bifurcation, i.e., the co-existence of cold and hot regions, provides temperatures high enough to produce the observed chromospheric UV emission and -- at the same time -- temperatures cold enough to allow the formation of molecules (e.g., carbon monoxide). Our 3-D model corroborates the finding by Carlsson & Stein (1994) that the chromospheric temperature rise of semi-empirical models does not necessarily imply an increase in the average gas temperature but can be explained by the presence of substantial spatial and temporal temperature inhomogeneities.

Recent numerical simulations by Freytag et al. of the outer convection envelope of the cool super-giant Betelgeuse, have shown that the fluctuations in the star's apparent luminosity may be caused by giant cell convection. These simulations bring forth the possibility of addressing another question; namely whether stars such as Betelgeuse may harbor magnetic activity. Taking the detailed numerical simulations of the star at face value, we have applied a kinematic dynamo analysis, to study whether the flow field of the super-giant may be able to amplify a weak seed magnetic field. We do indeed find a positive exponential growth rate of the magnetic energy. The possible Betelgeusian dynamo may be characterized as belonging to the class of so-called ``local small-scale dynamos'' (i.e. dynamos where rotation is not dominant), but this is a less meaningful designation in the case of Betelgeuse, since the field is both global and large-scale.

Based on detailed 2D numerical radiation hydrodynamics (RHD) calculations of time-dependent compressible convection, we have studied the dynamics and thermal structure of the convective surface layers of solar-type stars. The RHD models provide information about the convective efficiency in the superadiabatic region at the top of convective envelopes and predict the asymptotic value of the entropy of the deep, adiabatically stratified layers. This information is translated into an effective mixing-length parameter alpha suitable to construct standard stellar structure models. We validate the approach by a detailed comparison to helioseismic data. The grid of RHD models for solar metallicity comprises 58 simulation runs with a helium abundance of Y=0.28 in the range of effective temperatures 4300K < Teff < 7100K and gravities 2.54 < logg < 4.74. We find a moderate, nevertheless significant variation of alpha between about 1.3 for F-dwarfs and 1.75 for K-subgiants with a dominant dependence on Teff. In the close neighbourhood of the Sun we find a plateau where alpha remains almost constant. The internal accuracy of the calibration of alpha is estimated to be +/- 0.05 with a possible systematic bias towards lower values. An analogous calibration of the convection theory of Canuto & Mazzitelli gives a different temperature dependence but a similar variation of the free parameter. For the first time, values for the gravity-darkening exponent beta are derived independently of mixing-length theory: beta=0.07...0.10. We show that our findings are consistent with constraints from stellar stability considerations and provide compact fitting formulae for the calibrations.

One of the largest sources of uncertainty in the computation of globular cluster isochrones and hence in the age determination of globular clusters is the lack of a rigorous description of convection. Therefore, we calibrated the superadiabatic temperature gradient in the envelope of metal-poor low-mass stars according to the results from a new grid of 2D hydrodynamical models, which cover the Main Sequence and the lower Red Giant Branch of globular cluster stars. In practice, we still use for computing the evolutionary stellar models the traditional mixing length formalism, but we fix the mixing length parameter in order to reproduce the run of the entropy of the deeper adiabatic region of the stellar envelopes with effective temperature and gravity as obtained from the hydro-models. The detailed behaviour of the calibrated mixing length depends in a non-trivial way on the effective temperature, gravity and metallicity of the star. Nevertheless, the resulting isochrones for the relevant age range of galactic globular clusters have only small differences with respect to isochrones computed adopting a constant solar calibrated value of the mixing length. Accordingly, the age of globular clusters is reduced by 0.2 Gyr at most.

Based on radiation hydrodynamics modeling of stellar convection zones, a diffusion scheme has been devised describing the downward penetration of convective motions beyond the Schwarzschild boundary (overshoot) into the radiative interior. This scheme of exponential diffusive overshoot has already been successfully applied to AGB stars. Here we present an application to the Sun in order to determine the time scale and depth extent of this additional mixing, i.e. diffusive overshoot at the base of the convective envelope. We calculated the associated destruction of lithium during the evolution towards and on the main-sequence. We found that the slow-mixing processes induced by the diffusive overshoot may lead to a substantial depletion of lithium during the Sun's main-sequence evolution.

Results of nonradial, nonadiabatic pulsation calculations on hydrogen-rich white dwarf models are presented. In contrast to earlier attempts, the modeling builds on hydrodynamically simulated convective surface layers supplemented with standard interior models. Based on our stellar models and despite of various simple attempts to couple convection and pulsation we could not reproduce theoretically the presently adopted location of the observed blue edge of the ZZ Ceti variables. When the convective efficiency is high enough we found a sensitive dependence of the stability properties of the g-modes on the pulsational treatment of shear within the convection zone.