Search

Context. The complex dynamics of asymptotic giant branch (AGB) stars and the resulting stellar winds have a significant impact on the measurements of stellar parameters and amplify their uncertainties. Three-dimensional (3D) radiative hydrodynamic (RHD) simulations of convection suggest that convection-related structures at the surface of AGB star affect the photocentre displacement and the parallax uncertainty measured by Gaia. Aims. We explore the impact of the convection on the photocentre variability and aim to establish analytical laws between the photocentre displacement and stellar parameters to retrieve such parameters from the parallax uncertainty. Methods. We used a selection of 31 RHD simulations with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the Gaia G band [320-1050 nm]. From these maps, we calculated the photocentre position and temporal fluctuations. We then compared the synthetic standard deviation to the parallax uncertainty of a sample of 53 Mira stars observed with Gaia. Results. The simulations show a displacement of the photocentre across the surface ranging from 4 to 13 % of the corresponding stellar radius, in agreement with previous studies. We provide an analytical law relating the pulsation period of the simulations and the photocentre displacement as well as the pulsation period and stellar parameters. By combining these laws, we retrieve the surface gravity, the effective temperature, and the radius for the stars in our sample. Conclusions. Our analysis highlights an original procedure to retrieve stellar parameters by using both state-of-the-art 3D numerical simulations of AGB stellar convection and parallax observations of AGB stars. This will help us refine our understanding of these giants., Comment: 14 pages, 17 figures

Stellar convection plays an important role in atmospheric dynamics, wind formation and the mass-loss processes in Asymptotic Giant Branch (AGB) stars. However, a direct characterization of convective surface structures in terms of size, contrast, and life-span is quite challenging. Spatially resolving these features requires the highest angular resolution. In this work, we aim at characterizing the size of convective structures on the surface of the O-rich AGB star R Car to test different theoretical predictions, based on mixing-length theory from solar models. We used infrared low-spectral resolution (R~35) interferometric data in the H-band (~1.76 $\mu$m) with the instrument PIONIER at the Very Large Telescope Interferometer (VLTI) to image the star's surface at two epochs separated by ~6 years. Using a power spectrum analysis, we estimate the horizontal size of the structures on the surface of R Car. The sizes of the stellar disk, at different phases of a pulsation cycle, were obtained using parametric model-fitting in the Fourier domain. Our analysis supports that the sizes of the structures in R Car are correlated with variations of the pressure scale height in the atmosphere of the target, as predicted by theoretical models based on solar convective processes. We observe that these structures grow in size when the star expands within a pulsation cycle. While the information is still scarce, this observational finding highlights the role of convection in the dynamics of those objects. New interferometric imaging campaigns with the renewed capabilities of the VLTI are envisioned to expand our analysis to a larger sample of objects., Comment: Accepted by A&A on 15/05/2024, 11 pages, 7 figures

Red supergiant (RSGs) are cool massive stars in a late phase of their evolution when the stellar envelope becomes fully convective. They are the brightest stars in the universe at infrared light and can be detected in galaxies far beyond the Local Group, allowing for accurate determination of chemical composition of galaxies. The study of their physical properties is extremely important for various phenomena including the final fate of massive stars as type II supernovae and gravitational wave progenitors. We explore the well-studied nearby young stellar cluster chi Per. Using Gaia EDR3 data, we find the distance of the cluster (d = 2.260+-0.020 kpc). We then investigate the variability of the convection-related surface structure as a source for parallax measurement uncertainty. We use state-of-the-art 3D radiative hydrodynamics simulations with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps in the Gaia G photometric system. We calculate the variabiltiy, as a function of time, of the intensity-weighted mean from the synthetic maps. We then select the RSG stars in the cluster and compare their uncertainty on parallaxes to the predictions of photocentre displacements. The synthetic maps of RSG show extremely irregular and temporal variable surfaces due to convection-related dynamics. Consequentially, the position of the photo-center varies during Gaia measurements between 0.033 and 0.130 AU (up to 5% of the corresponding simulation stellar radius). We argue that the variability of the convection-related surface structures accounts for a substantial part of the Gaia EDR3 parallax error of the RSG sample. We suggest that the variation of the uncertainty on Gaia parallax could be exploited quantitatively using appropriate 3D simulations to extract, in a unique way, important information about the stellar dynamics and parameters of RSG stars., Comment: Press-release: https://www.mpa-garching.mpg.de/1054448/news20220530

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3 to 4 {\mu}m), M (4.5 to 5 {\mu}m), and N (8 to 13 {\mu}m) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics (RHD) simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes Optim3D and Radmc3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of the stellar surface and circumstellar environment across a very large spectral domain for the first time., Comment: Accepted for publication on Astronomy and Astrophysics

Despite being the best studied red supergiant star in our Galaxy, the physics behind the photometric variability and mass loss of Betelgeuse is poorly understood. Moreover, recently the star has experienced an unusual fading with its visual magnitude reaching a historical minimum. We investigate the nature of this event with the help of a recently developed tomographic technique. Tomography allows us to probe different depths in the stellar atmosphere and to recover the corresponding disk-averaged velocity field. We apply the tomographic method to a few-year time-series high-resolution spectroscopic observations of Betelgeuse in order to relate its atmospheric dynamics to the photometric variability. Our results show that a sudden increase of the molecular opacity in the cooler upper atmosphere of Betelgeuse is likely the reason of the observed unusual decrease of the star's brightness., Comment: A&A Letters, accepted for publication

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., Comment: Accepted for publication on Astronomy and Astrophysics

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..., Comment: 16 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

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., Comment: Accepted by A&A; 17 pages, 17 figures

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., Comment: 17 pages (main text), 12 pages (appendix), 15 figures, accepted for publication in MNRAS

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., Comment: 2 pages, 2 figures, accepted as Proceedings of IAU Symposium No. 343, 2018

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., Comment: 17 pages, 18 figures, accepted for publication in MNRAS

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., Comment: 13 pages, 12 figures, Accepted for publication in A&A

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., Comment: Letter to the Editor. Accepted for publication on Astronomy & Astrophysics

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), Comment: Accepted for publication by A&A, typo in the title corrected, cross reference to Paper II

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., Comment: 12 pages, 17 figures

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., Comment: accepted for publication in Astronomy & Astrophysics

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.

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., Comment: To be published as a proceeding of the 2016 CO5BOLD workshop in Memorie della SAIt Supplementi. 8 pages. 6 figures

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., Comment: Submitted to ApJ on 11-24-2015. Accepted on 5-17-2016 (Manuscript #: ApJ101257)

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., Comment: Accepted for publication in A&A. 16 pages, 10 figures

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., Comment: 4 pages, proceedings of the conference: "The Physics of Evolved Stars: A Conference Dedicated to the Memory of Olivier Chesneau", EAS Publications Series, 71-72, pages 237-242. arXiv admin note: text overlap with arXiv:1410.3868

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., Comment: 11 pages, 12 figures, accepted for publication in the Astrophysical Journal

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., Comment: 11 pages, 8 figures, accepted for publication in the Astrophysical Journal

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., Comment: Accepted for publication in A&A

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., Comment: 15 pages, 14 figures, accepted for publication in the Astrophysical Journal

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., Comment: Proceedings of the conference "Why Galaxies Care About AGB Stars III", Vienna July 28 - August 1, 2014. 10 pages, to appear on ASP Conference Series

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., Comment: 40 pages, 14 figures, accepted for publication in the Astrophysical Journal

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., Comment: 15 pages, 18 figures, 10 pages online appendix, accepted for publication in Astronomy and Astrophysics

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., Comment: 16 pages, 17 figures, 37 pages online appendix, accepted for publication in Astronomy and Astrophysics

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., Comment: 9 pages, 5 figures. Accepted for publication in Mem. S.A.It., proceedings of the 2nd CO5BOLD Workshop (Heidelberg, 2012)

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., Comment: Submitted to Mem. S. A. It. Suppl

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., Comment: 14 pages, 17 figures, accepted for publication in Astronomy and Astrophysics

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., Comment: 13 pages, to appear in A&A; model fluxes available from ftp://leda.as.utexas.edu/pub/callende/Gaia3D and soon from CDS

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., Comment: Acccepted to Astronomy and Astrophysics Review

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., Comment: Astronomy & Astrophysics, in press

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., Comment: Accepted for publication of A&A, 15 pages, 12 figures

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., Comment: 5 pages, 5 figures, accepted for publication in the Astronomy & Astrophysics Letters

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., Comment: 17 pages, accepted Astronomy & Astrophysics

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., Comment: Cool Star 16 conference paper submitted to the proceedings on Nov. 23rd 2010

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., Comment: Accepted for publication in Astronomy and Astrophysics

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., Comment: 22 pages, 15 figures, accepted for publication in Astronomy and Astrophysics, replaced with version with final format for table 3, typos corrected in table 4 and in equation B.1

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)., Comment: 11 pages, Accepted for publication on A&A

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), Comment: Astronomy and Astrophysics, accepted

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., Comment: 9 Pages, Accepted for publication on Astronomy & Astrophysics, two references updated

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., Comment: Contributed paper, to be published in the proceedings of IAU Symposium 265, eds. K. Cunha, M. Spite, and B. Barbuy, Cambridge University Press (CUP). 2 figures, 4 pages

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., Comment: Paper presented at IAU General Assembly Joint Discussion 10 "3D views on cool stellar atmospheres: theory meets observation"; to appear in a dedicated volume of the "Memorie della Societa Astronomica Italiana"; 4 pages, 3 figures

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., Comment: 16 pages with 18 Figures. Accepted for publication on Astronomy and Astrophysics

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., Comment: Accepted for publication on A&A