Search

The high-precision measurements of exoplanet transit light curves that are now available contain information about the planet properties, their orbital parameters, and stellar limb darkening (LD). Recent 3D magneto-hydrodynamical (MHD) simulations of stellar atmospheres have shown that LD depends on the photospheric magnetic field, and hence its precise determination can be used to estimate the field strength. Among existing LD laws, the uses of the simplest ones may lead to biased inferences, whereas the uses of complex laws typically lead to a large degeneracy among the LD parameters. We have developed a novel approach in which we use a complex LD model but with second derivative regularisation during the fitting process. Regularisation controls the complexity of the model appropriately and reduces the degeneracy among LD parameters, thus resulting in precise inferences. The tests on simulated data suggest that our inferences are not only precise but also accurate. This technique is used to re-analyse 43 transit light curves measured by the NASA Kepler and TESS missions. Comparisons of our LD inferences with the corresponding literature values show good agreement, while the precisions of our measurements are better by up to a factor of 2. We find that 1D non-magnetic model atmospheres fail to reproduce the observations while 3D MHD simulations are qualitatively consistent. The LD measurements, together with MHD simulations, confirm that Kepler-17, WASP-18, and KELT-24 have relatively high magnetic fields ($>200$ G). This study paves the way for estimating the stellar surface magnetic field using the LD measurements., Comment: 14 pages, 12 figures (including 2 in the appendix), 2 tables, accepted for publication in MNRAS

The Atacama Large Millimeter/submillimeter Array (ALMA) allows for solar observations in the wavelength range of 0.3$-$10 mm, giving us a new view of the chromosphere. The measured brightness temperature at various frequencies can be fitted with theoretical models of density and temperature versus height. We use the available ALMA and Mets\"ahovi measurements of selected solar structures (quiet sun (QS), active regions (AR) devoid of sunspots, and coronal holes (CH)). The measured QS brightness temperature in the ALMA wavelength range agrees well with the predictions of the semiempirical Avrett$-$Tian$-$Landi$-$Curdt$-$W\"ulser (ATLCW) model, better than previous models such as the Avrett$-$Loeser (AL) or Fontenla$-$Avrett$-$Loeser model (FAL). We scaled the ATLCW model in density and temperature to fit the observations of the other structures. For ARs, the fitted models require 9%$-$13% higher electron densities and 9%$-$10% higher electron temperatures, consistent with expectations. The CH fitted models require electron densities 2%$-$40% lower than the QS level, while the predicted electron temperatures, although somewhat lower, do not deviate significantly from the QS model. Despite the limitations of the one-dimensional ATLCW model, we confirm that this model and its appropriate adaptations are sufficient for describing the basic physical properties of the solar structures., Comment: 11 pages, 3 figures, 2 tables, accepted and published in Astronomische Nachrichten/Astronomical Notes

The Gaia satellite has already provided the astronomical community with three data releases, and the Radial Velocity Spectrometer (RVS) on board Gaia has provided the radial velocity for 33 million stars. When deriving the radial velocity from the RVS spectra, several stars are measured to have large values. To verify the credibility of these measurements, we selected some bright stars with the modulus of radial velocity in excess of 500\ to be observed with SOPHIE at OHP and UVES at VLT. This paper is devoted to investigating the chemical composition of the stars observed with UVES. We derived atmospheric parameters using Gaia photometry and parallaxes, and we performed a chemical analysis using the code. We find that the sample consists of metal-poor stars, although none have extremely low metallicities. The abundance patterns match what has been found in other samples of metal-poor stars selected irrespective of their radial velocities. We highlight the presence of three stars with low Cu and Zn abundances that are likely descendants of pair-instability supernovae. Two stars are apparently younger than 1\,Ga, and their masses exceed twice the turn-off mass of metal-poor populations. This makes it unlikely that they are blue stragglers because it would imply they formed from triple or multiple systems. We suggest instead that they are young metal-poor stars accreted from a dwarf galaxy. Finally, we find that the star RVS721 is associated with the Gjoll stream, which itself is associated with the Globular Cluster NGC\,3201., Comment: Astronomy and Astrophysics - A\&A, In press

The modeling of convection is a long standing problem in stellar physics. Up-to-now, all ad hoc models rely on a free parameter alpha (among others) which has no real physical justification and is therefore poorly constrained. However, a link exists between this free parameter and the entropy of the stellar adiabat. Prescriptions, derived from 3D stellar atmospheric models, are available that provide entropy as a function of stellar atmospheric parameters (effective temperature, surface gravity, chemical composition). This can provide constraints on alpha through the development of entropy-calibrated models. Several questions arise as these models are increasingly used. Which prescription should be used? How do uncertainties impact entropy-calibrated models? We aim to study the three existing prescriptions and determine which one should be used, and how. We implemented the entropy-calibration method into the stellar evolution code Cesam2k20 and performed comparisons with the Sun and the alpha Cen system. In addition, we used data from the CIFIST grid of 3D atmosphere models to evaluate the accuracy of the prescriptions. Of the three entropy prescriptions available, we determine which one best reproduces the entropies of the 3D models. We also demonstrate that the entropy obtained from this prescription should be corrected for the evolving chemical composition and for an entropy offset different between various EoS tables, following a precise procedure, otherwise classical parameters obtained from the models will be strongly biased. Finally, we also provide table with entropy of the adiabat of the CIFIST grid, as well as fits of these entropies. We performed a precise examination of entropy-calibrated modelling, and gave recommendations on which adiabatic entropy prescription to use, how to correct it and to implement the method into a stellar evolution code., Comment: 18 pages, 11 figures. Submitted to Astronomy and Astrophysics, recommended for publication with minor revisions

Binary molecules such as CO, OH, CH, CN, and C$_2$ are often used as abundance indicators in stars. These species are usually assumed to be formed in chemical equilibrium. The time-dependent effects of hydrodynamics can affect the formation and dissociation of these species and may lead to deviations from chemical equilibrium. We aim to model departures from chemical equilibrium in dwarf stellar atmospheres by considering time-dependent chemical kinetics alongside hydrodynamics and radiation transfer. We examine the effects of a decreasing metallicity and an altered C/O ratio on the chemistry when compared to the equilibrium state. We used the radiation-(magneto)hydrodynamics code CO5BOLD, and its own chemical solver to solve for the chemistry of 15 species and 83 reactions. The species were treated as passive tracers and were advected by the velocity field. The steady-state chemistry was also computed to isolate the effects of hydrodynamics. In most of the photospheres in the models we present, the mean deviations are smaller than $0.2$ dex, and they generally appear above $\log{\tau} = -2$. The deviations increase with height because the chemical timescales become longer with decreasing density and temperature. A reduced metallicity similarly results in longer chemical timescales and in a reduction in yield that is proportional to the drop in metallicity; a decrease by a factor $100$ in metallicity loosely corresponds to an increase by factor $100$ in chemical timescales. As both CH and OH are formed along reaction pathways to CO, the C/O ratio means that the more abundant element gives faster timescales to the constituent molecular species. Overall, the carbon enhancement phenomenon seen in very metal-poor stars is not a result of an improper treatment of molecular chemistry for stars up to a metallicity as low as [Fe/H] = $-3.0$., Comment: 29 pages

We use the ESPRESSO spectrograph at the Very Large Telescope to measure velocity shifts and gravitational redshifts of eight bona fide Hyades white dwarfs, with an accuracy better than 1.5 percent. By comparing the gravitational redshift measurements of the mass-to-radius ratio with the same ratios derived by fitting the \textit{Gaia} photometry with theoretical models, we find an agreement to better than one per cent. It is possible to reproduce the observed white dwarf cooling sequence and the trend of the mass-to-radius ratios as a function of colour using isochrones with ages between 725 and 800 Myr, tuned for the Hyades. One star, EGGR\,29, consistently stands out in all diagrams, indicating that it is possibly the remnant of a blue straggler. We also computed mass-to-radius ratios from published gravities and masses, determined from spectroscopy. The comparison between photometric and spectroscopic stellar parameters reveals that spectroscopic effective temperature and gravity are systematically larger than the photometric values. Spectroscopic mass-to-radius ratios disagree with those measured from gravitational redshift, indicating the presence of systematics affecting the white dwarf parameters derived from the spectroscopic analysis., Comment: MNRAS, accepted for pubblication

This study investigates and compares brightness and area of coronal bright points (CBPs) inside and outside of coronal holes (CHs) using the single-dish Band 6 observations by the Atacama Large Millimeter/submillimeter Array (ALMA), combined with extreme-ultraviolet (EUV) 193 $\overset{\circ}{\mathrm{A}}$ filtergrams obtained by the Atmospheric Imaging Assembly (AIA) and magnetograms obtained by the Helioseismic and Magnetic Imager (HMI), both on board Solar Dynamics Observatory (SDO). The CH boundaries were extracted from the SDO/AIA images using the Collection of Analysis Tools for Coronal Holes (CATCH) and CBPs were identified in the SDO/AIA, SDO/HMI, and ALMA data. Measurements of brightness and areas in both ALMA and SDO/AIA images were conducted for CBPs within CHs and quiet Sun regions outside CHs. A statistical analysis of the measured physical properties resulted in a lower average CBP brightness in both ALMA and SDO/AIA data for CBPs within the CHs. Depending on the CBP sample size, the difference in intensity for the SDO/AIA data, and brightness temperature for the ALMA data, between the CBPs inside and outside CHs ranged from between 2$\sigma$ and 4.5$\sigma$, showing a statistically significant difference between those two CBP groups. For CBP areas, CBPs within the CH boundaries showed smaller areas on average, with the observed difference between the two CBP groups between 1$\sigma$ and 2$\sigma$ for the SDO/AIA data, and up to 3.5$\sigma$ for the ALMA data, indicating that CBP areas are also significantly different. Given the measured properties, we conclude that the CBPs inside CHs tend to be less bright on average, but also smaller in comparison to those outside of CHs. This conclusion might point to the specific physical conditions and properties of the local CH region around a CBP limiting the maximum achievable intensity (temperature) and size of a CBP., Comment: 27 pages, 28 figures, accepted in Astronomy and Astrophysics

In this work, we present a photospheric solar silicon abundance derived using CO5BOLD model atmospheres and the LINFOR3D spectral synthesis code. Previous works have differed in their choice of a spectral line sample and model atmosphere as well as their treatment of observational material, and the solar silicon abundance has undergone a downward revision in recent years. We additionally show the effects of the chosen line sample, broadening due to velocity fields, collisional broadening, model spatial resolution, and magnetic fields. CO5BOLD model atmospheres for the Sun were used in conjunction with the LINFOR3D spectral synthesis code to generate model spectra, which were then fit to observations in the Hamburg solar atlas. We present a sample of 11 carefully selected lines (from an initial choice of 39 lines) in the optical and infrared, made possible with newly determined oscillator strengths for the majority of these lines. Our final sample includes seven optical Si I lines, three infrared Si I lines, and one optical Si II line. We derived a photospheric solar silicon abundance of $\log \epsilon_\mathrm{Si} = 7.57 \pm 0.04$, including a $-0.01$ dex correction from Non-Local Thermodynamic Equilibrium (NLTE) effects. Combining this with meteoritic abundances and previously determined photospheric abundances results in a metal mass fraction Z/X = $0.0220 \pm 0.0020$. We found a tendency of obtaining overly broad synthetic lines. We mitigated the impact of this by devising a de-broadening procedure. The over-broadening of synthetic lines does not substantially affect the abundance determined in the end. It is primarily the line selection that affects the final fitted abundance., Comment: 24 pages

The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022., Comment: 38 pages. A&A in press

The use of 3D hydrodynamical simulations of stellar surface convection for model atmospheres is computationally expensive. Although these models have been available for quite some time, their use is limited because of the lack of extensive grids of simulations and associated spectra. Our goal is to provide a method to interpolate spectra that can be applied to both 1D and 3D models, and implement it in a code available to the community. This tool will enable the routine use of 3D model atmospheres in the analysis of stellar spectra.} We have developed a code that makes use of radial basis functions to interpolate the spectra included in the CIFIST grid of 84 three-dimensional model atmospheres. Spectral synthesis on the hydrodynamical simulations was previously performed with the code ASS$\epsilon$T. We make a tool for the interpolation of 3D spectra available to the community. The code provides interpolated spectra and interpolation errors for a given wavelength interval, and a combination of effective temperature, surface gravity, and metallicity. In addition, it optionally provides graphical representations of the RMS and mean ratio between 1D and 3D spectra, and maps of the errors in the interpolated spectra across the parameter space., Comment: 14 pages, 9 figures, accepted for publication in A&A

Context. The TOPoS project has the goal to find and analyse Turn-Off (TO) stars of extremely low metallicity. To select the targets for spectroscopic follow-up at high spectral resolution, we have relied on low-resolution spectra from the Sloan Digital Sky Survey. Aims. In this paper we use the metallicity estimates we have obtained from our analysis of the SDSS spectra to construct the metallicity distribution function (MDF) of the Milky Way, with special emphasis on its metal-weak tail. The goal is to provide the underlying distribution out of which the TOPoS sample was extracted. Methods. We make use of SDSS photometry, Gaia photometry and distance estimates derived from the Gaia parallaxes to derive a metallicity estimate for a large sample of over 24 million TO stars. This sample is used to derive the metallicity bias of the sample for which SDSS spectra are available. Results. We determined that the spectroscopic sample is strongly biased in favour of metal-poor stars, as intended. A comparison with the unbiased photometric sample allows to correct for the selection bias. We select a sub-sample of stars with reliable parallaxes for which we combine the SDSS radial velocities with Gaia proper motions and parallaxes to compute actions and orbital parameters in the Galactic potential. This allows us to characterize the stars dynamically, and in particular to select a sub-sample that belongs to the Gaia-Sausage-Enceladus (GSE) accretion event. We are thus able to provide also the MDF of GSE. Conclusions. The metal-weak tail derived in our study is very similar to that derived in the H3 survey and in the Hamburg/ESO Survey. This allows us to average the three MDFs and provide an error bar for each metallicity bin. Inasmuch the GSE structure is representative of the progenitor galaxy that collided with the Milky Way, that galaxy appears to be strongly deficient in metal-poor stars compared to the Milky Way, suggesting that the metal-weak tail of the latter has been largely formed by accretion of low mass galaxies rather than massive galaxies, such as the GSE progenitor., Comment: Astronomy and Astrophysics - A\&A, EDP Sciences, in press

Aims. The main aim of the present analysis is to decipher (i) the small-scale bright features in solar images of the quiet Sun and active regions obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) and (ii) the ALMA correspondence of various known chromospheric structures visible in the H-alpha images of the Sun. Methods. Small-scale ALMA bright features in the quiet Sun region were analyzed using single-dish ALMA observations (1.21 mm, 248 GHz) and in an active region using interferometric ALMA measurements (3 mm, 100 GHz). With the single-dish observations, a full-disk solar image is produced, while interferometric measurements enable the high-resolution reconstruction of part of the solar disk, including the active region. The selected quiet Sun and active regions are compared with the H-alpha (core and wing sum), EUV, and soft X-ray images and with the magnetograms. Results. In the quiet Sun region, enhanced emission seen in the ALMA is almost always associated with a strong line-of-sight (LOS) magnetic field. Four coronal bright points were identified, while other small-scale ALMA bright features are most likely associated with magnetic network elements and plages. In the active region, in 14 small-scale ALMA bright features randomly selected and compared with other images, we found five good candidates for coronal bright points, two for plages, and five for fibrils. Two unclear cases remain: a fibril or a jet, and a coronal bright point or a plage. A comparison of the H-alpha core image and the 3 mm ALMA image of the analyzed active region showed that the sunspot appears dark in both images (with a local ALMA radiation enhancement in sunspot umbra), the four plage areas are bright in both images and dark small H-alpha filaments are clearly recognized as dark structures of the same shape also in ALMA., Comment: 12 pages, 4 figures, to be published in Astronomy and Astrophysics

The development of space-borne missions has significantly improved the quality of the measured spectra of solar-like oscillators. Their $p$-mode line profiles can now be resolved, and the asymmetries inferred for a variety of stars other than the Sun. However, it has been known for a long time that the asymmetries of solar $p$-modes are reversed between the velocity and the intensity spectra. Understanding the origin of this reversal is necessary in order to use asymmetries as a tool for seismic diagnosis. For stars other than the Sun, only the intensity power spectrum is sufficiently resolved to allow for an estimation of mode asymmetries. We recently developed an approach designed to model and predict these asymmetries in the velocity power spectrum of the Sun and to successfully compare them to their observationally derived counterparts. In this paper we expand our model and predict the asymmetries featured in the intensity power spectrum. We find that the shape of the mode line profiles in intensity is largely dependent on how the oscillation-induced variations of the radiative flux are treated, and that modelling it realistically is crucial to understanding asymmetry reversal. Perturbing a solar-calibrated grey atmosphere model, and adopting the quasi-adiabatic framework as a first step, we reproduce the asymmetries observed in the solar intensity spectrum for low-frequency modes. We conclude that, unlike previously thought, it is not necessary to invoke an additional mechanism (e.g. non-adiabatic effects, coherent non-resonant background signal) to explain asymmetry reversal. This additional mechanism is necessary, however, to explain asymmetry reversal for higher-order modes., Comment: 9 pages, 4 figures. Accepted for publication in A&A

The General Theory of Relativity predicts the redshift of spectral lines in the solar photosphere, as a consequence of the gravitational potential of the Sun. This effect can be measured from a solar disk-integrated flux spectrum of the Sun's reflected light on solar system bodies. The laser frequency comb (LFC) calibration system attached to the HARPS spectrograph offers the possibility to perform an accurate measurement of the solar gravitational redshift (GRS) by observing the Moon or other solar system bodies. We have analysed the line shift observed in Fe absorption lines from five high-quality HARPS-LFC spectra of the Moon. We select an initial sample of 326 photospheric Fe lines in the spectral range 476-585 nm and measure their line positions and equivalent widths (EWs). Accurate line shifts are derived from the wavelength position of the core of the lines compared with the laboratory wavelengths. We fit the observed spectral Fe lines using CO$^5$BOLD 3D synthetic profiles. Convective motions in the solar photosphere do not affect the line cores of Fe lines stronger than about $\sim 150$ mA. In our sample, only 15 FeI lines have EWs in the range $150 <$ EW(mA) $< 550$, providing a measurement of the solar GRS at $639\pm14$ ${\rm m\;s^{-1}}$, consistent with the expected theoretical value on Earth of $\sim 633.1$ ${\rm m\;s^{-1}}$. A final sample of about 97 weak Fe lines with EW $<180$ mA allows us to derive a mean global line shift of $638\pm6$ ${\rm m\;s^{-1}}$ in agreement with the theoretical solar GRS. These are the most accurate measurements of the solar GRS so far. Ultrastable spectrographs calibrated with the LFC over a larger spectral range, such as HARPS or ESPRESSO, together with a further improvement on the laboratory wavelengths, could provide a more robust measurement of the solar GRS and further tests for the 3D hydrodynamical models., Comment: Accepted for publication in Astronomy & Astrophysics

In this paper, we aim to develop a predictive model for solar radial $p$-mode line profiles in the velocity spectrum. Unlike the approach favoured by prior studies, this model is not described by free parameters and we do not use fitting procedures to match the observations. Instead, we use an analytical turbulence model coupled with constraints extracted from a 3D hydrodynamic simulation of the solar atmosphere. We then compare the resulting asymmetries with their observationally derived counterpart. We find that stochastic excitation localised beneath the mode upper turning point generates negative asymmetry for $\nu < \nu_\text{max}$ and positive asymmetry for $\nu > \nu_\text{max}$. On the other hand, stochastic excitation localised above this limit generates negative asymmetry throughout the $p$-mode spectrum. As a result of the spatial extent of the source of excitation, both cases play a role in the total observed asymmetries. By taking this spatial extent into account and using a realistic description of the spectrum of turbulent kinetic energy, both a qualitative and quantitative agreement can be found with solar observations perfoemed by the GONG network. We also find that the impact of the correlation between acoustic noise and oscillation is negligible for mode asymmetry in the velocity spectrum., Comment: 26 pages (19 without appendices), 14 figures. Accepted for publication in A&A

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

We present high-resolution and high-quality UVES spectroscopic data of the metal-poor double-lined spectroscopic binary CS 22876--032 ([Fe/H] $\sim -3.7$ dex), with the goal to derive the $^6$Li/$^7$Li isotopic ratio by analysing the \ion{Li}{i} $\lambda$~670.8~nm doublet. We coadd all 28 useful spectra normalised and corrected for radial velocity to the rest frame of the primary star. We fit the Li profile with a grid of the 3D-NLTE synthetic spectra, to take into account the line profile asymmetries induced by stellar convection, and perform Monte Carlo simulations to evaluate the uncertainty of the fit of the Li line profile. We check that the veiling factor does not affect the derived isotopic ratio, $^6$Li/$^7$Li, and only modifies the Li abundance, A(Li), by about 0.15~dex. The best fit of the Li profile of the primary star provides A(Li)~$ = 2.17 \pm 0.01$~dex and $^6$Li/$^7$Li~$=8^{+2}_{-5}$\% at 68\% confidence level. In addition, we improve the Li abundance of the secondary star at A(Li)~$= 1.55 \pm 0.04$~dex, which is about 0.6~dex lower than that of the primary star. The analysis of the Li profile of the primary star is consistent with no detection of $^6$Li and provides an upper-limit to the isotopic ratio of $^6$Li/$^7$Li~$< 10$\% at this very low metallicity, about 0.5~dex lower in metallicity than previous attempts for detection of $^6$Li in extremely metal poor stars. These results do not solve or worsen the cosmological $^7$Li problem, nor support the need for non standard $^6$Li production in the early Universe., Comment: Accepted for publication in A&A

Context. It is possible to accurately measure the masses of the white dwarfs (WDs) in the Hyades cluster using gravitational redshift, because the radial velocity of the stars can be obtained independently of spectroscopy from astrometry and the cluster has a low velocity dispersion. Aims. We aim to obtain an accurate measurement of the Hyades WD masses by determining the mass-to-radius ratio (M/R) from the observed gravitational redshift, and to compare them with masses derived from other methods. Methods. We analyse archive high-resolution UVES-VLT spectra of six WDs belonging to the Hyades to measure their Doppler shift, from which M/R is determined after subtracting the astrometric radial velocity. We estimate the radii using Gaia photometry as well as literature data. Results. The M/R error associated to the gravitational redshift measurement is about 5%. The radii estimates, evaluated with different methods, are in very good agreement, though they can differ by up to 4% depending on the quality of the data. The masses based on gravitational redshift are systematically smaller than those derived from other methods, by a minimum of $\sim 0.02$ up to 0.05 solar masses. While this difference is within our measurement uncertainty, the fact that it is systematic indicates a likely real discrepancy between the different methods. Conclusions. We show that the M/R derived from gravitational redshift measurements is a powerful tool to determine the masses of the Hyades WDs and could reveal interesting properties of their atmospheres. The technique can be improved by using dedicated spectrographs, and can be extended to other clusters, making it unique in its ability to accurately and empirically determine the masses of WDs in open clusters. At the same time we prove that gravitational redshift in WDs agrees with the predictions of stellar evolution models to within a few percent., Comment: Accepted for publication in Astronomy & Astrophysics

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

We will study the formation history of the Milky Way, and the earliest phases of its chemical enrichment, with a sample of more than 1.5 million stars at high galactic latitude. Elemental abundances of up to 20 elements with a precision of better than 0.2 dex will be derived for these stars. The sample will include members of kinematically coherent substructures, which we will associate with their possible birthplaces by means of their abundance signatures and kinematics, allowing us to test models of galaxy formation. Our target catalogue is also expected to contain 30,000 stars at a metallicity of less than one hundredth that of the Sun. This sample will therefore be almost a factor of 100 larger than currently existing samples of metal-poor stars for which precise elemental abundances are available (determined from high-resolution spectroscopy), enabling us to study the early chemical evolution of the Milky Way in unprecedented detail., Comment: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Oxygen abundance measurements are important for understanding stellar structure and evolution. Measured in Cepheids, they further provide clues on the metallicity gradient and chemo-dynamical evolution in the Galaxy. However, most of the abundance analyses of Cepheids to date have been based on one-dimensional (1D) hydrostatic model atmospheres. Here, we test the validity of this approach for the key oxygen abundance diagnostic, the \ion{O}{I} $777\,\mathrm{nm}$~triplet lines. We carry out 2D non-LTE radiative transfer clculations across two different 2D radiation hydrodynamics simulations of Cepheid atmospheres, having stellar parameters of $T_\mathrm{eff}= 5600$ K, solar chemical compositions, and $\log\,g= 1.5$ and $2.0$, corresponding to pulsation periods of 9 and 3 days, respectively. We find that the 2D non-LTE versus 1D LTE abundance differences range from $-1.0$~dex to $-0.25$~dex depending on pulsational phase. The 2D non-LTE versus 1D non-LTE abundance differences range from $-0.2$~dex to $0.8$~dex. The abundance differences are smallest when the Cepheid atmospheres are closest to hydrostatic equilibrium, corresponding to phases of around $0.3$ to $0.8$, and we recommend these phases for observers deriving the oxygen abundance from \ion{O}{I} $777\,\mathrm{nm}$ triplet with 1D hydrostatic models., Comment: 9 pages, 10 figures; Published in Astronomy and Astrophysics

In this work we use solar observations with the ALMA radio telescope at the wavelength of 1.21 mm. The aim of the analysis is to improve understanding of the solar chromosphere, a dynamic layer in the solar atmosphere between the photosphere and corona. The study has an observational and a modeling part. In the observational part full-disc solar images are analyzed. Based on a modified FAL atmospheric model, radiation models for various observed solar structures are developed. Finally, the observational and modeling results are compared and discussed., Comment: 4 pages, presented at The 20th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, July 29 - Aug 3 2018, Boston / Cambridge, USA

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

We study the accuracy of spectroscopic RVs by comparing spectroscopic and astrometric RVs for stars of the Hyades open cluster. Rather accurate astrometric RVs are available for the Hyades' stars, based on Hipparcos and on the first Gaia data release. We obtained HARPS spectra for a large sample of Hyades stars, and homogeneously analysed them. After cleaning the sample from binaries, RV variables, and outliers, 71 stars remained. The distribution of the observed RV difference (spectroscopic -- astrometric) is skewed and depends on the star right ascension. This is consistent with the Hyades cluster rotating at 42.3 +/- 10 m/s/pc. The two Hyades giants in the sample show, as predicted by gravitational redshift (GR), a spectroscopic RV that is blue-shifted with respect to the dwarfs, and the empirical GR slope is of 626 +/- 131 m/s, in agreement with the theoretical prediction. The difference between spectroscopic and astrometric RVs is very close to zero, within the uncertainties. In particular, the mean difference is of -33 m/s and the median is of -16 m/s when considering the Gaia-based RVs (corrected for cluster rotation), with a sigma of 347 m/s, very close to the expected cluster velocity dispersion. We also determine a new value of the cluster centroid spectroscopic RV: 39.36 +/- 0.26 km/s. The spectroscopic RV measurements are expected, from simulations, to depend on stellar rotation, but our data do not confirm these predictions. We finally discuss other phenomena that can influence the RV difference, such as cluster expansion, stellar activity, general relativity, and Galactic potential. Clusters within the reach of current telescopes are expected to show differences of several hundreds m/s, depending on their position in the Galaxy., Comment: 18 pages, 9 figures, accepted for publication in MNRAS

The CoRoT and Kepler missions provided a wealth of high-quality data for solar-like oscillations. To make the best of such data for seismic inferences, we need theoretical models with precise near-surface structure, which has significant influence on solar-like oscillation frequencies. The mixing-length parameter, $\alpha$, is a key factor for the near-surface structure. In the convection formulations used in evolution codes, the $\alpha$ is a free parameter that needs to be properly specified. We calibrated $\alpha$ values by matching entropy profiles of 1D envelope models with those of 3D CO$^5$BOLD models. For such calibration, previous works concentrated on the classical mixing-length theory (MLT). Here we also analyzed the full spectrum turbulence (FST) models. For the atmosphere part in the 1D models, we use the Eddington grey $T(\tau)$ relation and the one with the solar-calibrated Hopf-like function. For both the MLT and FST models with a mixing length $l=\alpha H_p$, calibrated $\alpha$ values increase with increasing $g$ or decreasing $T_{\rm eff}$. For the FST models, we also calibrated values of $\alpha^*$ defined as $l=r_{\rm top}-r+\alpha^*H_{p,{\rm top}}$. $\alpha^*$ is found to increase with $T_{\rm eff}$ and $g$. As for the correspondence to the 3D models, the solar Hopf-like function gives a photospheric-minimum entropy closer to a 3D model than the Eddington $T(\tau)$. The structure below the photosphere depends on the convection model. However, not a single convection model gives the best correspondence since the averaged 3D quantities are not necessarily related via an EOS. Although the FST models with $l=r_{\rm top}-r+\alpha^*H_{p,{\rm top}}$ are found to give the frequencies closest to the solar observed ones, a more appropriate treatment of the top part of the 1D convective envelope is necessary., Comment: 14 pages, 18 figures, accepted to Astronomy and Astrophysics

Extremely metal-poor stars are keys to understand the early evolution of our Galaxy. The ESO large programme TOPoS has been tailored to analyse a new set of metal-poor turn-off stars, whereas most of the previously known extremely metal-poor stars are giant stars. Sixty five turn-off stars (preselected from SDSS spectra) have been observed with the X-Shooter spectrograph at the ESO VLT Unit Telescope 2, to derive accurate and detailed abundances of magnesium, silicon, calcium, iron, strontium and barium. We analysed medium-resolution spectra (R ~ 10 000) obtained with the ESO X-Shooter spectrograph and computed the abundances of several alpha and neutron-capture elements using standard one-dimensional local thermodynamic equilibrium (1D LTE) model atmospheres. Our results confirms the super-solar [Mg/Fe] and [Ca/Fe] ratios in metal-poor turn-off stars as observed in metal-poor giant stars. We found a significant spread of the [alpha/Fe] ratios with several stars showing sub-solar [Ca/Fe] ratios. We could measure the abundance of strontium in 12 stars of the sample, leading to abundance ratios [Sr/Fe] around the Solar value. We detected barium in two stars of the sample. One of the stars (SDSS J114424-004658) shows both very high [Ba/Fe] and [Sr/Fe] abundance ratios (>1 dex)., Comment: Accepted for publication in A&A : 17 pages, 6 figures

The CoRoT and Kepler missions have provided high-quality measurements of the frequency spectra of solar-like pulsators, enabling us to probe stellar interiors with a very high degree of accuracy by comparing the observed and modeled frequencies. However, the frequencies computed with 1D models suffer from systematic errors related to the poor modeling of the uppermost layers of stars. These biases are what is commonly named the near surface effect. The dominant effect is related to the turbulent pressure that modifies the hydrostatic equilibrium and thus the frequencies. This has already been investigated using grids of 3D RMHD simulations, which also were used to constrain the parameters of the empirical correction models. However, the effect of metallicity has not been considered so far. We study the impact of metallicity on the surface effect across the HR diagram, and provide a method for accounting for it when using the empirical correction models. We computed a grid of patched 1D stellar models with the stellar evolution code CESTAM in which poorly modeled surface layers have been replaced by averaged stratification computed with the 3D RMHD code CO5BOLD. We found that metallicity has a strong impact on the surface effect: keeping T_eff and log g constant, the frequency residuals can vary by up to a factor two. Therefore, the influence of metallicity cannot be neglected. We found that a correct way of accounting for it is to consider the surface Rosseland mean opacity. It allowed us to give a physically-grounded justification as well as a scaling relation for the frequency differences at nu_max as a function of T_eff, log g and kappa. Finally, we provide prescriptions for the fitting parameters of the correction functions. We show that the impact of metallicity through the Rosseland mean opacity must be taken into account when studying and correcting the surface effect., Comment: 15 pages, 7 figures. Accepted for publication in A&A

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

We investigate possible relations between the abundances of zinc and the light elements sodium, magnesium, and potassium in the atmospheres ofred giant branch (RGB) stars of the Galactic globular cluster 47~Tuc and study connections between the chemical composition and dynamical properties of the cluster RGB stars. The abundance of zinc was determined in 27 RGB stars of 47~Tuc using 1D local thermal equilibrium (LTE) synthetic line profile fitting to the high-resolution 2dF spectra obtained with the Anglo-Australian Telescope (AAT). Synthetic spectra used in the fitting procedure were computed with the SYNTHE code and 1D ATLAS9 stellar model atmospheres. The average 1D~LTE zinc-to-iron abundance ratio and its RMS variations due to star-to-star abundance spread determined in the sample of 27 RGB stars is $\langle{\rm[Zn/Fe]}\rangle^{\rm 1D~LTE}=0.11\pm0.09$. We did not detect any statistically significant relations between the abundances of zinc and those of light elements. Neither did we find any significant correlation or anticorrelation between the zinc abundance in individual stars and their projected distance from the cluster center. Finally, no statistically significant relation between the absolute radial velocities of individual stars and the abundance of zinc in their atmospheres was detected. The determined average [Zn/Fe]$^{\rm~1DLTE}$ ratio agrees well with those determined in this cluster in earlier studies and nearly coincides with that of Galactic field stars at this metallicity. All these results suggest that nucleosynthesis of zinc and light elements proceeded in separate, unrelated pathways in 47~Tuc., Comment: 8 pages, 8 figures. Accepted to A&A

We determined abundances of Mg and K in the atmospheres of 53 (Mg) and 75 (K) turn-off (TO) stars of the Galactic globular cluster 47 Tuc. The obtained abundances, together with those of Li, O, and Na, were used to search for possible relations between the abundances of K and other light elements, Li, O, Na, and Mg, as well as the connections between the chemical composition of TO stars and their kinematical properties. Abundances of Mg and K were determined using archival VLT FLAMES/GIRAFFE spectra, in combination with the 1D NLTE spectral synthesis methodology. Spectral line profiles were computed with the MULTI code, using ATLAS9 stellar model atmospheres. We also utilized 3D hydrodynamical CO5BOLD and 1D hydrostatic LHD model atmospheres for computing 3D-1D LTE abundance corrections for the spectral lines of Mg and K, in order to assess the influence of convection on their formation in the atmospheres of TO stars. The determined average abundance-to-iron ratios and their RMS variations due to star-to-star abundance spreads were $\langle{\rm[Mg/Fe]}\rangle^{\rm 1D NLTE}=0.47\pm0.12$, and $\langle{\rm [K/Fe]}\rangle^{\rm 1D NLTE}=0.39\pm0.09$. Although the data suggest the possible existence of a weak correlation in the [K/Fe]-[Na/Fe] plane, its statistical significance is low. No statistically significant relations between the abundance of K and other light elements were detected. Also, we did not find any significant correlations or anti-correlations between the [Mg/Fe] and [K/Fe] ratios and projected distance from the cluster center. Similarly, no relations between the absolute radial velocities of individual stars and abundances of Mg and K in their atmospheres were detected. The 3D-1D abundance corrections were found to be small ($\leq 0.1$\,dex) for the lines of Mg and K used in this study, thus indicating that the influence of convection on their formation is small., Comment: 13 pages, 7 figures. Accepted to A&A

We studied the influence of convection on the spectral energy distributions, photometric magnitudes, and colour indices of different types of stars across the H-R diagram. The 3D hydrodynamical CO5BOLD, averaged <3D>, and 1D hydrostatic LHD model atmospheres were used to compute spectral energy distributions of stars on the main sequence (MS), main sequence turn-off (TO), subgiant branch (SGB), and red giant branch (RGB), in each case at two different effective temperatures and two metallicities, [M/H]=0.0 and -2.0. Using the obtained spectral energy distributions, we calculated photometric magnitudes and colour indices in the broad-band Johnson-Cousins $UBVRI$ and 2MASS $JHK_{\rm s}$, and the medium-band Str\"{o}mgren $uvby$ photometric systems. The 3D-1D differences in photometric magnitudes and colour indices are small in both photometric systems and typically do not exceed $\pm0.03$ mag. Only in the case of the coolest giants located on the upper RGB are the differences in the $U$ and $u$ bands able reach $\approx-0.2$ mag at [M/H]=0.0 and $\approx-0.1$ mag at [M/H]=-2.0. Generally, the 3D-1D differences are largest in the blue-UV part of the spectrum and decrease towards longer wavelengths. They are also sensitive to the effective temperature and are significantly smaller in hotter stars. Metallicity also plays a role and leads to slightly larger 3D-1D differences at [M/H]=0.0. All these patterns are caused by a complex interplay between the radiation field, opacities, and horizontal temperature fluctuations that occur due to convective motions in stellar atmospheres. Although small, the 3D-1D differences in the magnitudes and colour indices are nevertheless comparable to or larger than typical photometric uncertainties and may therefore cause non-negligible systematic differences in the estimated effective temperatures., Comment: 12 pages, 14 figures, accepted for publication by A&A

Extremely metal-poor stars provide us with indirect information on the first generations of massive stars. The TOPoS survey has been designed to increase the census of these stars and to provide a chemical inventory that is as detailed as possible. Seven of the most iron-poor stars have been observed with the UVES spectrograph at the ESO VLT Kueyen 8.2m telescope to refine their chemical composition. We analysed the spectra based on 1D LTE model atmospheres, but also used 3D hydrodynamical simulations of stellar atmospheres. We measured carbon in six of the seven stars: all are carbon-enhanced and belong to the low-carbon band, defined in the TOPoS II paper. We measured lithium (A(Li)=1.9) in the most iron-poor star (SDSS J1035+0641, [Fe/H] < -5.2 ). We were also able to measure Li in three stars at [Fe/H]~ -4.0, two of which lie on the Spite plateau. We confirm that SDSS J1349+1407 is extremely rich in Mg, but not in Ca. It is also very rich in Na. Several of our stars are characterised by low alpha-to-iron ratios. The lack of high-carbon band stars at low metallicity can be understood in terms of evolutionary timescales of binary systems. The detection of Li in SDSS J1035+0641 places a strong constraint on theories that aim at solving the cosmological lithium problem. The Li abundance of the two warmer stars at [Fe/H]~ -4.0 places them on the Spite plateau, while the third, cooler star, lies below. We argue that this suggests that the temperature at which Li depletion begins increases with decreasing [Fe/H]. SDSS J1349+1407 may belong to a class of Mg-rich EMP stars. We cannot assess if there is a scatter in alpha-to-iron ratios among the EMP stars or if there are several discrete populations. However, the existence of stars with low alpha-to-iron ratios is supported by our observations., Comment: A&A, in press

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

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.

We present our latest 3D model atmospheres for carbon-enhanced metal-poor (CEMP) stars computed with the CO5BOLD code. The stellar parameters are representative of hot turn-off objects (Teff ~ 6250 K, log g=4.0, [Fe/H]=-3.0). The main purpose of these models is to investigate the role of 3D effects on synthetic spectra of the CH G-band (4140-4400 A), the CN BX-band (3870-3890 A), and several UV OH transitions (3122-3128 A). By comparison with the synthetic spectra from standard 1D model atmospheres (assuming local thermodynamic equilibrium, LTE), we derive 3D abundance corrections for carbon and oxygen of up to -0.5 and -0.7 dex, respectively., Comment: poster paper presented at IAU Symposium 334

Due to their volatile nature, when sulfur and zinc are observed in external galaxies, their determined abundances represent the gas-phase abundances in the interstellar medium. This implies that they can be used as tracers of the chemical enrichment of matter in the Universe at high redshift. Comparable observations in stars are more difficult and, until recently, plagued by small number statistics. We wish to exploit the Gaia ESO Survey (GES) data to study the behaviour of sulfur and zinc abundances of a large number of Galactic stars, in a homogeneous way. By using the UVES spectra of the GES sample, we are able to assemble a sample of 1301 Galactic stars, including stars in open and globular clusters in which both sulfur and zinc were measured. We confirm the results from the literature that sulfur behaves as an alpha-element. We find a large scatter in [Zn/Fe] ratios among giant stars around solar metallicity. The lower ratios are observed in giant stars at Galactocentric distances less than 7.5 kpc. No such effect is observed among dwarf stars, since they do not extend to that radius. Given the sample selection, giants and dwarfs are observed at different Galactic locations, and it is plausible, and compatible with simple calculations, that Zn-poor giants trace a younger population more polluted by SN Ia yields. It is necessary to extend observations in order to observe both giants and dwarfs at the same Galactic location. Further theoretical work on the evolution of zinc is also necessary., Comment: 15 pages, 17 figures, accepted for publication in A&A

We study the abundances of Na, Mg, and K in the atmospheres of 32 RGB stars in the Galactic globular cluster (GGC) 47 Tuc, with the goal to investigate the possible existence of Na-K and Mg-K correlations/anti-correlations, similar to those that were recently discovered in two other GGCs, NGC 2419 and 2808. The abundances of K, Na, and Mg were determined using high-resolution 2df spectra obtained with the AAT. The 1D NLTE abundance estimates were obtained using 1D hydrostatic ATLAS9 model atmospheres and spectral line profiles synthesized with the MULTI package. We also used 3D hydrodynamical CO5BOLD and 1D hydrostatic LHD model atmospheres to compute 3D-1D LTE abundance corrections, $\Delta_{\rm 3D-1D~LTE}$, for the spectral lines of Na, Mg, and K used in our study. These abundance corrections were used to understand the role of convection in the formation of spectral lines, as well as to estimate the differences in the abundances obtained with the 3D hydrodynamical and 1D hydrostatic model atmospheres. The average element-to-iron abundance ratios and their RMS variations due to star-to-star abundance spreads determined in our sample of RGB stars were $\langle{\rm [Na/Fe]}\rangle^{\rm 1D~NLTE}=0.42\pm0.13$, $\langle{\rm [Mg/Fe]}\rangle^{\rm 1D~NLTE}=0.41\pm0.11$, and $\langle{\rm [K/Fe]}\rangle^{\rm 1D~NLTE}=0.05\pm0.14$. We found no statistically significant relations between the abundances of the three elements studied here. Also, there were no abundance trends with the distance from the cluster center, nor any statistically significant relations between the abundance/abundance ratios and absolute radial velocities of individual stars. All these facts suggest the similarity of K abundance in stars that belong to different generations in 47 Tuc which, in turn, may hint that evolution of K in this particular cluster was unrelated to the nucleosynthesis of Na and/or Mg., Comment: 12 pages, 6 figures. Accepted to A&A. Minor adjustments in the text and the Figs. 4-6

Context. Tighter constraints on metal-poor stars we observe are needed to better understand the chemical processes of the early Universe. Computing a stellar spectrum in 3D allows one to model complex stellar behaviours, which cannot be replicated in 1D. Aims. We examine the effect that the intrinsic CNO abundances have on a 3D model structure and the resulting 3D spectrum synthesis. Methods. Model atmospheres were computed in 3D for three distinct CNO chemical compositions using the CO 5 BOLD model atmosphere code, and their internal structures were examined. Synthetic spectra were computed from these models using Linfor3D and they were compared. New 3D abundance corrections for the G-band and a selection of UV OH lines were also computed. Results. The varying CNO abundances change the metal content of the 3D models. This had an effect on the model structure and the resulting synthesis. However, it was found that the C/O ratio had a larger effect than the overall metal content of a model. Conclusions. Our results suggest that varying the C/O ratio has a substantial impact on the internal structure of the 3D model, even in the hot turn-off star models explored here. This suggests that bespoke 3D models, for specific CNO abundances should be sought. Such effects are not seen in 1D at these temperature regimes., Comment: Accepted for publication as a letter in A&A. 4 figures, 2 tables, 6 pages

The space-borne missions have provided a wealth of highly accurate data. However, our inability to properly model the upper-most region of solar-like stars prevents us from making the best of these observations. This problem is called "surface effect" and a key ingredient to solve it is turbulent pressure for the computation of both the equilibrium models and the oscillations. While 3D hydrodynamic simulations help to include properly the turbulent pressure in the equilibrium models, the way this surface effect is included in the computation of stellar oscillations is still subject to uncertainties. We aim at determining how to properly include the effect of turbulent pressure and its Lagrangian perturbation in the adiabatic computation of the oscillations. We also discuss the validity of the gas-gamma model (GGM) and reduced gamma model (RGM) approximations, which have been used to compute adiabatic oscillations of equilibrium models including turbulent pressure. We use a patched model of the Sun with an inner part constructed by a 1D stellar evolution code (CESTAM) and an outer part by the 3D hydrodynamical code (CO$^5$BOLD). Then, the adiabatic oscillations are computed using the ADIPLS code for the GGM and RGM and with the MAD code imposing the adiabatic condition on an existing time-dependent convection (TDC) formalism. We show that the computation of the oscillations using the TDC formalism in the adiabatic limit improves significantly the agreement with the observed frequencies compared to the GGM and RGM. Of the components of the turbulent pressure perturbation, the perturbation of the density and advection term is found to contribute most to the frequency shift. We propose a formalism to evaluate the frequency shift due to the inclusion of the term with the turbulent pressure perturbation in the variational principle in order to extrapolate our result to other stars., Comment: 11 pages, 7 figures, accepted for publication in Astronomy and Astrophysics

Although oxygen is an important tracer of the early Galactic evolution, its abundance trends with metallicity are still relatively poorly known at [Fe/H] < -2.5. This is in part due to a lack of reliable oxygen abundance indicators in the metal-poor stars, in part due to shortcomings in 1D LTE abundance analyses. In this study we determined the oxygen abundance in the metal-poor halo giant HD 122563 using a 3D hydrodynamical CO5BOLD model atmosphere. Our main goal was to understand whether a 3D LTE analysis may help to improve the reliability of oxygen abundances determined from OH UV lines in comparison to those obtained using standard 1D LTE methodology. The oxygen abundance in HD 122563 was determined using 71 OH UV lines located in the wavelength range between 308-330 nm. The analysis was done using a high-resolution VLT UVES spectrum with a 1D LTE spectral line synthesis performed using the SYNTHE package and classical ATLAS9 model atmosphere. Subsequently, a 3D hydrodynamical CO5BOLD, and 1D hydrostatic LHD model atmospheres were used in order to compute 3D-1D abundance corrections. For this, the microturbulence velocity used with the 1D LHD model atmosphere was derived from the hydrodynamical CO5BOLD model atmosphere. As in previous studies, we found trends of the 1D LTE oxygen abundances determined from OH UV lines with line parameters, such as the line excitation potential and the line strength. These trends become significantly less pronounced in 3D LTE. Using OH UV lines we determined a 3D LTE oxygen abundance in HD 122563 of A(O) = 6.23 +/- 0.13. This is in fair agreement with the oxygen abundance obtained from OH IR lines, A(O) = 6.39 +/- 0.11, but it is noticeably lower than that determined using the forbidden [OI] line, A(O) = 6.53 +/- 0.15. While the exact cause for this discrepancy remains unclear, it is very likely that non-LTE effects may play a decisive role here., Comment: 18 pages, 13 figures, 6 tables. Accepted for publication in A&A

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

One of the primary objectives of the TOPoS survey is to search for the most metal-poor stars. Our search has led to the discovery of one of the most iron-poor objects known, SDSS\,J092912.32+023817.0. This object is a multiple system, in which two components are clearly detected in the spectrum. We have analysed 16 high-resolution spectra obtained using the UVES spectrograph at the ESO 8.2m VLT telescope to measure radial velocities and determine the chemical composition of the system. Cross correlation of the spectra with a synthetic template yields a double-peaked cross-correlation function (CCF) for eight spectra, and in one case there is evidence for the presence of a third peak. Chemical analysis of the spectrum obtained by averaging all the spectra for which the CCF showed a single peak found that the iron abundance is [Fe/H]=-4.97. The system is also carbon enhanced with [C/Fe]=+3.91 (A(C)=7.44). From the permitted oxygen triplet we determined an upper limit for oxygen of [O/Fe]<+3.52 such that C/O>1.3. We are also able to provide more stringent upper limits on the Sr and Ba abundances ([Sr/Fe]<+0.70, and [Ba/Fe]<+1.46, respectively)., Comment: Accepted for publication in A&A Letter

Recent developments in the three-dimensional (3D) spectral synthesis code Linfor3D have meant that, for the first time, large spectral wavelength regions, such as molecular bands, can be synthesised with it in a short amount of time. A detailed spectral analysis of the synthetic G-band for several dwarf turn-off-type 3D atmospheres (5850 <= T_eff [K] <= 6550, 4.0 <= log g <= 4.5, -3.0 <= [Fe/H] <= -1.0) was conducted, under the assumption of local thermodynamic equilibrium. We also examine carbon and oxygen molecule formation at various metallicity regimes and discuss the impact it has on the G-band. Using a qualitative approach, we describe the different behaviours between the 3D atmospheres and the traditional one-dimensional (1D) atmospheres and how the different physics involved inevitably leads to abundance corrections, which differ over varying metallicities. Spectra computed in 1D were fit to every 3D spectrum to determine the 3D abundance correction. Early analysis revealed that the CH molecules that make up the G-band exhibited an oxygen abundance dependency; a higher oxygen abundance leads to weaker CH features. Nitrogen abundances showed zero impact to CH formation. The 3D corrections are also stronger at lower metallicity. Analysis of the 3D corrections to the G-band allows us to assign estimations of the 3D abundance correction to most dwarf stars presented in the literature. The 3D corrections suggest that A(C) in CEMP stars with high A(C) would remain unchanged, but would decrease in CEMP stars with lower A(C). It was found that the C/O ratio is an important parameter to the G-band in 3D. Additional testing confirmed that the C/O ratio is an equally important parameter for OH transitions under 3D. This presents a clear interrelation between the carbon and oxygen abundances in 3D atmospheres through their molecular species, which is not seen in 1D., Comment: 19 pages, 13 figures, 4 tables. Accepted for publication in A&A

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

Different simplified approaches are used to account for the non-local thermodynamic equilibrium (NLTE) effects with 3D hydrodynamical model atmospheres. In certain cases, chemical abundances are derived in 1D NLTE and corrected for the 3D effects by adding 3D-1D LTE abundance corrections (3D+NLTE approach). Alternatively, average <3D> model atmospheres are sometimes used to substitute for the full 3D hydrodynamical models. We tested whether the results obtained using these simplified schemes (i.e., 3D+NLTE, <3D> NLTE) may reproduce those derived using the full 3D NLTE computations. The tests were made using 3D hydrodynamical CO5BOLD model atmospheres of the main sequence (MS), main sequence turn-off (TO), subgiant (SGB), and red giant branch (RGB) stars, all at [M/H]=0.0 and -2.0. Our goal was to investigate the role of 3D and NLTE effects on the formation of the 670.8 nm lithium line by assessing strengths of synthetic 670.8 nm line profiles, computed using 3D/1D NLTE/LTE approaches. Our results show that Li 670.8 nm line strengths obtained using different methodologies differ only slightly in most of the models at solar metallicity. However, the line strengths predicted with the 3D NLTE and 3D+NLTE approaches become significantly different at subsolar metallicities. At [M/H]=-2.0, this may lead to (3D NLTE)-(3D+NLTE) differences in the predicted lithium abundance of ~0.46 and ~0.31 dex in the TO and RGB stars, respectively. On the other hand, NLTE line strengths computed with the average <3D> and 1D model atmospheres are similar to those obtained with the full 3D NLTE approach for MS, TO, SGB, and RGB stars, at all metallicities; 3D-<3D> and 3D-1D differences in the predicted abundances are always less than ~0.04 dex and ~0.08 dex, respectively. However, neither of the simplified approaches can reliably substitute 3D NLTE spectral synthesis when precision is required., Comment: 11 pages, 6 figures, accepted for publication in A&A

In recent years a number of poorly studied chemical elements, such as phosphorus, sulphur, and strontium, have received special attention as important tracers of the Galactic chemical evolution. By exploiting the capabilities of the infrared echelle spectrograph GIANO mounted at the Telescopio Nazionale Galileo, we acquired high resolution spectra of four Galactic dwarf stars spanning the metallicity range between about one-third and twice the solar value. We performed a detailed feasibility study about the effectiveness of the P, S, and Sr line diagnostics in the Y band between 1.03 and 1.10 microm. Accurate chemical abundances have been derived using one-dimensional model atmospheres computed in local thermodynamic equilibrium (LTE). We computed the line formation assuming LTE for P, while we performed non-LTE analysis to derive S and Sr abundances. We were able to derive phosphorus abundance for three stars and an upper limit for one star, while we obtained the abundance of sulphur and strontium for all of the stars. We find [P/Fe] and [S/Fe] abundance ratios consistent with solar-scaled or slightly depleted values, while the [Sr/Fe] abundance ratios are more scattered (by +/-0.2 dex) around the solar-scaled value. This is fully consistent with previous studies using both optical and infrared spectroscopy. We verified that high-resolution, Y-band spectroscopy as provided by GIANO is a powerful tool to study the chemical evolution of P, S, and Sr in dwarf stars., Comment: To appear on A&A

The space-borne missions have provided us with a wealth of high-quality observational data that allows for seismic inferences of stellar interiors. This requires the computation of precise and accurate theoretical frequencies, but imperfect modeling of the uppermost stellar layers introduces systematic errors. To overcome this problem, an empirical correction has been introduced by Kjeldsen et al. (2008, ApJ, 683, L175) and is now commonly used for seismic inferences. Nevertheless, we still lack a physical justification allowing for the quantification of the surface-effect corrections. We used a grid of these simulations computed with the CO$^5$BOLD code to model the outer layers of solar-like stars. Upper layers of the corresponding 1D standard models were then replaced by the layers obtained from the horizontally averaged 3D models. The frequency differences between these patched models and the 1D standard models were then calculated using the adiabatic approximation and allowed us to constrain the Kjeldsen et al. power law, as well as a Lorentzian formulation. We find that the surface effects on modal frequencies depend significantly on both the effective temperature and the surface gravity. We further provide the variation in the parameters related to the surface-effect corrections using their power law as well as a Lorentzian formulation. Scaling relations between these parameters and the elevation (related to the Mach number) is also provided. The Lorentzian formulation is shown to be more robust for the whole frequency spectrum, while the power law is not suitable for the frequency shifts in the frequency range above $\nu_{\rm max}$., Comment: 11 pages, 14 figures, 4 tables; accepted for publication in Astronomy & Astrophysics

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