Source author record

H. -G. Ludwig

H. -G. Ludwig appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

88works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

88 published item(s)

preprint2025arXiv

Unveiling the nature of HE 0107-5240

The vast majority of the most iron-poor stars in the Galaxy exhibit a strong carbon enhancement, with C/H ratios only about two orders of magnitude below solar. This unusual chemical composition likely reflects the properties of the gas cloud from which these stars formed, having been enriched by one, or at most a few, supernovae. A remarkable member of this stellar class, HE 0107-5240 with [Fe/H]=-5.56, has been identified as part of a binary system. To constrain its orbital parameters, radial velocity monitoring has been carried out using the ESPRESSO spectrograph. Radial velocities were derived using cross-correlation with a template, taking advantage of the strong G-band feature. Combining all observations yielded a high signal-to-noise spectrum, which has been used to refine our understanding of the stellar chemical composition. Additionally, a co-added UVES spectrum in the blue was used to complement the wavelength coverage of ESPRESSO. Observations of HE 0107-5240 over a span of more than four years have yielded a revised orbital period of about 29 years. Updated elemental abundances have been determined for Sc, Cr, Co, and, tentatively, Al, along with a new upper limit for Be. The iron abundance has been derived from ionised Fe lines. Significant upper limits have been established for Li, Si, and Sr. The star is confirmed to be a long-period binary. Iron abundances derived from neutral and ionised lines are consistent with local thermodynamical equilibrium (LTE) assumption, casting doubt on published deviation from LTE corrections for Fe for this star. The heavy elements Sr and Ba remain undetected, confirming the classification of HE 0107-5240 as a carbon enhanced metal-poor and non enhanced in heavy elements (CEMP-no) star and supporting the absence of an n-capture element plateau at the lowest metallicities.

preprint2022arXiv

Interpolation of spectra from 3D model atmospheres

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

preprint2020arXiv

Modelling the asymmetries of the Sun's radial $p$-mode line profiles

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 $ν< ν_\text{max}$ and positive asymmetry for $ν> ν_\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.

preprint2020arXiv

Velocity-intensity asymmetry reversal of solar radial p-modes

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.

preprint2016arXiv

An in-depth spectroscopic examination of molecular bands from 3D hydrodynamical model atmospheres I. Formation of the G-band in metal-poor dwarf stars

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.

preprint2016arXiv

Enhanced methods for computing spectra from CO5BOLD models using Linfor3D. Molecular bands in metal-poor stars

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.

preprint2016arXiv

Solar science with the Atacama Large Millimeter/submillimeter Array - A new view of our Sun

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

preprint2016arXiv

The chemical composition of red giants in 47 Tucanae. II. Magnesium isotopes and pollution scenarios

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

preprint2016arXiv

TOPoS: III. An ultra iron-poor multiple CEMP system

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

preprint2015arXiv

3D Model Atmospheres for Extremely Low-Mass White Dwarfs

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.

preprint2015arXiv

A three-dimensional hydrodynamical line profile analysis of iron lines and barium isotopes in HD140283

Heavy-elements, i.e. those beyond the iron peak, mostly form via two neutron capture processes: the s- and r-process. Metal-poor stars should contain fewer isotopes that form via the s-process, according to currently accepted theory. It has been shown in several investigations that theory and observation do not agree well, raising questions on the validity of either the methodology or the theory. We analyse the metal-poor star HD140283, for which we have a high quality spectrum. We test whether a 3D LTE stellar atmosphere and spectrum synthesis code permits a more reliable analysis of the iron abundance and barium isotope ratio than a 1D LTE analysis. Using 3D model atmospheres, we examine 91 iron lines of varying strength and formation depth. This provides us with the star's rotational speed. With this, we model the barium isotope ratio by exploiting the hyperfine structure of the singly ionised 4554 resonance line, and study the impact of the uncertainties in the stellar parameters. HD140283's vsini = 1.65 +/- 0.05 km/s. Barium isotopes under the 3D paradigm show a dominant r-process signature as 77 +/- 6 +/- 17% of barium isotopes form via the r-process, where errors represent the assigned random and systematic errors, respectively. We find that 3D LTE fits reproduce iron line profiles better than those in 1D, but do not provide a unique abundance (within the uncertainties). However, we demonstrate that the isotopic ratio is robust against this shortcoming. Our barium isotope result agrees well with currently accepted theory regarding the formation of the heavy-elements during the early Galaxy. The improved fit to the asymmetric iron line profiles suggests that the current state of 3D LTE modelling provides excellent simulations of fluid flows. However, the abundances they provide are not yet self-consistent. This may improve with NLTE considerations and higher resolution models.

preprint2015arXiv

GIANO Y-band spectroscopy of dwarf stars: Phosphorus, Sulphur, and Strontium abundances

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.

preprint2015arXiv

Lithium spectral line formation in stellar atmospheres. The impact of convection and NLTE effects

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.

preprint2015arXiv

On The Evolution of Magnetic White Dwarfs

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.

preprint2015arXiv

Stellar science from a blue wavelength range - A possible design for the blue arm of 4MOST

From stellar spectra, a variety of physical properties of stars can be derived. In particular, the chemical composition of stellar atmospheres can be inferred from absorption line analyses. These provide key information on large scales, such as the formation of our Galaxy, down to the small-scale nucleosynthesis processes that take place in stars and supernovae. By extending the observed wavelength range toward bluer wavelengths, we optimize such studies to also include critical absorption lines in metal-poor stars, and allow for studies of heavy elements (Z>38) whose formation processes remain poorly constrained. In this context, spectrographs optimized for observing blue wavelength ranges are essential, since many absorption lines at redder wavelengths are too weak to be detected in metal-poor stars. This means that some elements cannot be studied in the visual-redder regions, and important scientific tracers and science cases are lost. The present era of large public surveys will target millions of stars. Here we describe the requirements driving the design of the forthcoming survey instrument 4MOST, a multi-object spectrograph commissioned for the ESO VISTA 4m-telescope. We focus here on high-density, wide-area survey of stars and the science that can be achieved with high-resolution stellar spectroscopy. Scientific and technical requirements that governed the design are described along with a thorough line blending analysis. For the high-resolution spectrograph, we find that a sampling of >2.5 (pixels per resolving element), spectral resolution of 18000 or higher, and a wavelength range covering 393-436 nm, is the most well-balanced solution for the instrument. A spectrograph with these characteristics will enable accurate abundance analysis (+/-0.1 dex) in the blue and allow us to confront the outlined scientific questions. (abridged)

preprint2015arXiv

Surface-effect corrections for solar-like oscillations using 3D hydrodynamical simulations

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 $ν_{\rm max}$.

preprint2015arXiv

The photospheric solar oxygen project: III. Investigation of the centre-to-limb variation of the 630nm [OI]-NiI blend

The solar photospheric abundance of oxygen is still a matter of debate. For about ten years some determinations have favoured a low oxygen abundance which is at variance with the value inferred by helioseismology. Among the oxygen abundance indicators, the forbidden line at 630nm has often been considered the most reliable even though it is blended with a NiI line. In Papers I and Paper II of this series we reported a discrepancy in the oxygen abundance derived from the 630nm and the subordinate [OI] line at 636nm in dwarf stars, including the Sun. Here we analyse several, in part new, solar observations of the the centre-to-limb variation of the spectral region including the blend at 630nm in order to separate the individual contributions of oxygen and nickel. We analyse intensity spectra observed at different limb angles in comparison with line formation computations performed on a CO5BOLD 3D hydrodynamical simulation of the solar atmosphere. The oxygen abundances obtained from the forbidden line at different limb angles are inconsistent if the commonly adopted nickel abundance of 6.25 is assumed in our local thermodynamic equilibrium computations. With a slightly lower nickel abundance, A(Ni)~6.1, we obtain consistent fits indicating an oxygen abundance of A(O)=8.73+/-0.05. At this value the discrepancy with the subordinate oxygen line remains. The derived value of the oxygen abundance supports the notion of a rather low oxygen abundance in the solar hotosphere. However, it is disconcerting that the forbidden oxygen lines at 630 and 636nm give noticeably different results, and that the nickel abundance derived here from the 630nm blend is lower than expected from other nickel lines.

preprint2015arXiv

The photospheric solar oxygen project: IV. 3D-NLTE investigation of the 777 nm triplet lines

The solar photospheric oxygen abundance is still widely debated. Adopting the solar chemical composition based on the "low" oxygen abundance, as determined with the use of three-dimensional (3D) hydrodynamical model atmospheres, results in a well-known mismatch between theoretical solar models and helioseismic measurements that is so far unresolved. We carry out an independent redetermination of the solar oxygen abundance by investigating the center-to-limb variation of the OI IR triplet lines at 777 nm in different sets of spectra with the help of detailed synthetic line profiles based on 3D hydrodynamical CO5BOLD model atmospheres and 3D non-LTE line formation calculations with NLTETD. The idea is to simultaneously derive the oxygen abundance,A(O), and the scaling factor SH that describes the cross-sections for inelastic collisions with neutral hydrogen relative the classical Drawin formula. The best fit of the center-to-limb variation of the triplet lines achieved with the CO5BOLD 3D solar model is clearly of superior quality compared to the line profile fits obtained with standard 1D model atmospheres. Our best estimate of the 3D non-LTE solar oxygen abundance is A(O) = 8.76 +/- 0.02, with the scaling factor SH in the range between 1.2 and 1.8. All 1D non-LTE models give much lower oxygen abundances, by up to -0.15 dex. This is mainly a consequence of the assumption of a $μ$-independent microturbulence.

preprint2015arXiv

Three-dimensional hydrodynamical CO5BOLD model atmospheres of red giant stars. IV. Oxygen diagnostics in extremely metal-poor red giants with infrared OH lines

Context. Although oxygen is an important tracer of Galactic chemical evolution, measurements of its abundance in the atmospheres of the oldest Galactic stars are still scarce and rather imprecise. At the lowest end of the metallicity scale, oxygen can only be measured in giant stars and in most of cases such measurements rely on a single forbidden [O I] 630 nm line that is very weak and frequently blended with telluric lines. Although molecular OH lines located in the ultraviolet and infrared could also be used for the diagnostics, oxygen abundances obtained from the OH lines and the [O I] 630 nm line are usually discrepant to a level of ~0.3-0.4 dex. Aims. We study the influence of convection on the formation of the infrared (IR) OH lines and the forbidden [O I] 630 nm line in the atmospheres of extremely metal-poor (EMP) red giant stars. Methods. We used high-resolution and high signal-to-noise ratio spectra of four EMP red giant stars obtained with the VLT CRIRES spectrograph. For each EMP star, 4-14 IR OH vibrational-rotational lines located in the spectral range of 1514-1548 and 1595-1632 nm were used to determine oxygen abundances by employing standard 1D LTE abundance analysis methodology. We then corrected the 1D LTE abundances obtained from each individual OH line for the 3D hydrodynamical effects. Results. We find that the influence of convection on the formation of [O I] 630 nm line in the atmospheres of EMP giants studied here is minor, which leads to very small 3D-1D abundance corrections (< -0.01 dex). On the contrary, IR OH lines are strongly affected by convection and thus the abundance corrections for these lines are significant, 3D-1D ~ -0.2 ... -0.3 dex. These abundance corrections do indeed bring the 1D LTE oxygen abundances of EMP red giants obtained using IR OH lines into better agreement with those determined from the [O I] 630 nm line.

preprint2014arXiv

Abundances of lithium, oxygen, and sodium in the turn-off stars of Galactic globular cluster 47 Tuc

We aim to determine abundances of Li, O and Na in a sample of of 110 turn-off (TO) stars, in order to study the evolution of light elements in this cluster and to put our results in perspective with observations of other globular and open clusters, as well as with field stars. We use medium resolution spectra obtained with the GIRAFFE spectrograph at the ESO 8.2m Kueyen VLT telescope and use state of the art 1D model atmospheres and NLTE line transfer to determine the abundances. We also employ CO5BOLD hydrodynamical simulations to assess the impact of stellar granulation on the line formation and inferred abundances. Our results confirm the existence of Na-O abundance anti-correlation and hint towards a possible Li-O anti-correlation in the TO stars of 47 Tuc. We find no convincing evidence supporting the existence of Li-Na correlation. The obtained 3D NLTE mean lithium abundance in a sample of 94 TO stars where Li lines were detected reliably, $\langle A({\rm Li})_{\rm 3D~NLTE}\rangle = 1.78 \pm 0.18$ dex, appears to be significantly lower than what is observed in other globular clusters. At the same time, star-to-star spread in Li abundance is also larger than seen in other clusters. The highest Li abundance observed in 47 Tuc is about 0.1 dex lower than the lowest Li abundance observed among the un-depleted stars of the metal-poor open cluster NGC 2243. The lithium abundances in 47 Tuc, when put into context with observations in other clusters and field stars, suggest that stars that are more metal-rich than [FeH] \sim -1.0 experience significant lithium depletion during their lifetime on the main sequence, while the more metal-poor stars do not. Rather strikingly, our results suggest that initial lithium abundance with which the star was created may only depend on its age (the younger the star, the higher its Li content) and not on its metallicity.

preprint2014arXiv

Calibration of the Mixing-Length Theory for Convective White Dwarf Envelopes

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.

preprint2014arXiv

Clues on the Galactic evolution of sulphur from star clusters

(Abridged) The abundances of alpha-elements are a powerful diagnostic of the star formation history and chemical evolution of a galaxy. Sulphur, being moderately volatile, can be reliably measured in the interstellar medium (ISM) of damped Ly-alpha galaxies and extragalactic HII regions. Measurements in stars of different metallicity in our Galaxy can then be readily compared to the abundances in external galaxies. Such a comparison is not possible for Si or Ca that suffer depletion onto dust in the ISM. Furthermore, studying sulphur is interesting because it probes nucleosynthetic conditions that are very different from those of O or Mg. The measurements in star clusters are a reliable tracers of the Galactic evolution of sulphur. We find <A(S)>NLTE=6.11+/-0.04 for M 4, <A(S)>NLTE=7.17+/-0.02 for NGC 2477, and <A(S)>NLTE=7.13+/-0.06 for NGC 5822. For the only star studied in Trumpler 5 we find A(S)NLTE=6.43+/-0.03 and A(S)LTE=6.94+/-0.05. Our measurements show that, by and large, the S abundances in Galactic clusters trace reliably those in field stars. The only possible exception is Trumpler 5, for which the NLTE sulphur abundance implies an [S/Fe] ratio lower by roughly 0.4 dex than observed in field stars of comparable metallicity, even though its LTE sulphur abundance is in line with abundances of field stars. Moreover the LTE sulphur abundance is consistent only with the abundance of another alpha-element, Mg, in the same star, while the low NLTE value is consistent with Si and Ca. The S abundances in our sample of stars in clusters imply that the clusters are chemically homogeneous for S within 0.05 dex.

preprint2014arXiv

Oxygen in the early Galaxy: OH lines as tracers of oxygen abundance in extremely metal-poor giant stars

Oxygen is a powerful tracer element of Galactic chemical evolution. Unfortunately, only a few oxygen lines are available in the ultraviolet-infrared stellar spectra for the reliable determination of its abundance. Moreover, oxygen abundances obtained using different spectral lines often disagree significantly. In this contribution we therefore investigate whether the inadequate treatment of convection in 1D hydrostatic model atmospheres used in the abundance determinations may be responsible for this disagreement. For this purpose, we used VLT CRIRES spectra of three EMP giants, as well as 3D hydrodynamical CO$^5$BOLD and 1D hydrostatic LHD model atmospheres, to investigate the role of convection in the formation of infrared (IR) OH lines. Our results show that the presence of convection leads to significantly stronger IR OH lines. As a result, the difference in the oxygen abundance determined from IR OH lines with 3D hydrodynamical and classical 1D hydrostatic model atmospheres may reach -0.2 ... -0.3 dex. In case of the three EMP giants studied here, we obtain a good agrement between the 3D LTE oxygen abundances determined by us using vibrational-rotational IR OH lines in the spectral range of 1514-1626 nm, and oxygen abundances determined from forbidden [O I] 630 nm line in previous studies.

preprint2014arXiv

The chemical composition of red giants in 47 Tucanae I: Fundamental parameters and chemical abundance patterns

Context: The study of chemical abundance patterns in globular clusters is of key importance to constrain the different candidates for intra-cluster pollution of light elements. Aims: We aim at deriving accurate abundances for a large range of elements in the globular cluster 47 Tucanae (NGC 104) to add new constraints to the pollution scenarios for this particular cluster, expanding the range of previously derived element abundances. Methods: Using tailored 1D LTE atmospheric models together with a combination of equivalent width measurements, LTE, and NLTE synthesis we derive stellar parameters and element abundances from high-resolution, high signal-to-noise spectra of 13 red giant stars near the tip of the RGB. Results: We derive abundances of a total 27 elements (O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ru, Ba, La, Ce, Pr, Nd, Eu, Dy). Departures from LTE were taken into account for Na, Al and Ba. We find a mean [Fe/H] = $-0.78\pm0.07$ and $[α/{\rm Fe}]=0.34\pm0.03$ in good agreement with previous studies. The remaining elements show good agreement with the literature, but the inclusion of NLTE for Al has a significant impact on the behaviour of this key element. Conclusions: We confirm the presence of an Na-O anti-correlation in 47 Tucanae found by several other works. Our NLTE analysis of Al shifts the [Al/Fe] to lower values, indicating that this may be overestimated in earlier works. No evidence for an intrinsic variation is found in any of the remaining elements.

preprint2014arXiv

White Dwarfs In The UKIRT Infrared Deep Sky Survey Data Release 9

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.

preprint2013arXiv

Carbon-enhanced metal-poor stars: the most pristine objects?

Carbon-enhanced metal poor stars (CEMP) form a significant proportion of the metal-poor stars, their origin is not well understood. Three very metal-poor C-rich turnoff stars were selected from the SDSS survey, observed with the ESO VLT (UVES) to precisely determine the element abundances. In turnoff stars (unlike giants) the carbon abundance has not been affected by mixing with deep layers and is therefore easier to interpret. The analysis was performed with 1D LTE static model atmospheres. When available, non-LTE corrections were applied to the classical LTE abundances. The 3D effects on the CH and CN molecular bands were computed using hydrodynamical simulations of the stellar atmosphere (CO5BOLD) and are found to be very important. To facilitate a comparison with previous results, only 1D abundances are used in the discussion. The abundances (or upper limits) of the elements enable us to place these stars in different CEMP classes. The carbon abundances confirm the existence of a plateau at A(C)= 8.25 for [Fe/H] \geq -3.4. The most metal-poor stars ([Fe/H] < -3.4) have significantly lower carbon abundances, suggesting a lower plateau at A(C) \approx 6.5. Detailed analyses of a larger sample of very low metallicity carbon-rich stars are required to confirm (or refute) this possible second plateau and specify the behavior of the CEMP stars at very low metallicity.

preprint2013arXiv

Convective line shifts for the Gaia RVS from the CIFIST 3D model atmosphere grid

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.

preprint2013arXiv

Fluorine Abundances of Galactic Low-Metallicity Giants

With abundances and 2σ upper limits of fluorine (F) in seven metal-poor field giants, nucleosynthesis of stellar F at low metallicity is discussed. The measurements are derived from the HF(1-0) R9 line at 23358Å using nearinfrared K-band high-resolution spectra obtained with CRIRES at the Very Large Telescope. The sample reaches lower metallicities than previous studies on F of field giants, ranging from [Fe/H] = -1.56 down to -2.13. Effects of three-dimensional model atmospheres on the derived F and O abundances are quantitatively estimated and shown to be insignificant for the program stars. The observed F yield in the form of [F/O] is compared with two sets of Galactic chemical evolution models, which quantitatively demonstrate the contribution of Type II supernova (SN II) ν-process and asymptotic giant branch/Wolf-Rayet stars. It is found that at this low-metallicity region, models cannot well predict the observed distribution of [F/O], while the observations are better fit by models considering an SN II ν-process with a neutrino energy of E_ν = 3 x 10^53 erg. Our sample contains HD 110281, a retrograde orbiting low-α halo star, showing a similar F evolution as globular clusters. This supports the theory that such halo stars are possibly accreted from dwarf galaxy progenitors of globular clusters in the halo.

preprint2013arXiv

Granulation properties of giants, dwarfs, and white dwarfs from the CIFIST 3D model atmosphere grid

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.

preprint2013arXiv

Is the Sun Lighter than the Earth? Isotopic CO in the Photosphere, Viewed through the Lens of 3D Spectrum Synthesis

We consider the formation of solar infrared (2-6 micron) rovibrational bands of carbon monoxide (CO) in CO5BOLD 3D convection models, with the aim to refine abundances of the heavy isotopes of carbon (13C) and oxygen (18O,17O), to compare with direct capture measurements of solar wind light ions by the Genesis Discovery Mission. We find that previous, mainly 1D, analyses were systematically biased toward lower isotopic ratios (e.g., R23= 12C/13C), suggesting an isotopically "heavy" Sun contrary to accepted fractionation processes thought to have operated in the primitive solar nebula. The new 3D ratios for 13C and 18O are: R23= 91.4 +/- 1.3 (Rsun= 89.2); and R68= 511 +/- 10 (Rsun= 499), where the uncertainties are 1 sigma and "optimistic." We also obtained R67= 2738 +/- 118 (Rsun= 2632), but we caution that the observed 12C17O features are extremely weak. The new solar ratios for the oxygen isotopes fall between the terrestrial values and those reported by Genesis (R68= 530, R6= 2798), although including both within 2 sigma error flags, and go in the direction favoring recent theories for the oxygen isotope composition of Ca-Al inclusions (CAI) in primitive meteorites. While not a major focus of this work, we derive an oxygen abundance of 603 +/- 9 ppm (relative to hydrogen; 8.78 on the logarithmic H= 12 scale). That the Sun likely is lighter than the Earth, isotopically speaking, removes the necessity to invoke exotic fractionation processes during the early construction of the inner solar system.

preprint2013arXiv

Isotope Spectroscopy

The measurement of isotopic ratios provides a privileged insight both into nucleosynthesis and into the mechanisms operating in stellar envelopes, such as gravitational settling. In this article, we give a few examples of how isotopic ratios can be determined from high-resolution, high-quality stellar spectra. We consider examples of the lightest elements, H and He, for which the isotopic shifts are very large and easily measurable, and examples of heavier elements for which the determination of isotopic ratios is more difficult. The presence of 6Li in the stellar atmospheres causes a subtle extra depression in the red wing of the 7Li 670.7 nm doublet which can only be detected in spectra of the highest quality. But even with the best spectra, the derived $^6$Li abundance can only be as good as the synthetic spectra used for their interpretation. It is now known that 3D non-LTE modelling of the lithium spectral line profiles is necessary to account properly for the intrinsic line asymmetry, which is produced by convective flows in the atmospheres of cool stars, and can mimic the presence of 6Li. We also discuss briefly the case of the carbon isotopic ratio in metal-poor stars, and provide a new determination of the nickel isotopic ratios in the solar atmosphere.

preprint2013arXiv

Micro- and macroturbulence predictions from CO5BOLD 3D stellar atmospheres

We present an overview of the current status of our efforts to derive the microturbulence and macroturbulence parameters (ximic and ximac) from the CIFIST grid of CO5BOLD 3D model atmospheres as a function of the basic stellar parameters Teff, log g, and [M/H]. The latest results for the Sun and Procyon show that the derived microturbulence parameter depends significantly on the numerical resolution of the underlying 3D simulation, confirming that `low-resolution' models tend to underestimate the true value of ximic. Extending the investigation to twelve further simulations with different Teff, log g, and [M/H], we obtain a first impression of the predicted trend of ximic over the Hertzsprung-Russell diagram: in agreement with empirical evidence, microturbulence increases towards higher effective temperature and lower gravity. The metallicity dependence of ximic must be interpreted with care, since it also reflects the deviation between the 1D and 3D photospheric temperature stratifications that increases systematically towards lower metallicity.

preprint2013arXiv

On the seismic scaling relations $Δν- \barρ$ and $ν_{\rm max}-ν_{\rm c}$

Scaling relations between asteroseismic quantities and stellar parameters are essential tools for studying stellar structure and evolution. We will address two of them, namely, the relation between the large frequency separation ($Δν$) and the mean density ($\barρ$) as well as the relation between the frequency of the maximum in the power spectrum of solar-like oscillations ($ν_{\rm max}$) and the cut-off frequency ($ν_{\rm c}$). For the first relation, we will consider the possible sources of uncertainties and explore them with the help of a grid of stellar models. For the second one, we will show that the basic physical picture is understood and that departure from the observed relation arises from the complexity of non-adiabatic processes involving time-dependent treatment of convection. This will be further discussed on the basis of a set of 3D hydrodynamical simulation of surface convection.

preprint2013arXiv

Oxygen spectral line synthesis: 3D non-LTE with CO5BOLD hydrodynamical model atmospheres

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.

preprint2013arXiv

Pure-hydrogen 3D model atmospheres of cool white dwarfs

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.

preprint2013arXiv

Signs of atmospheric inhomogeneities in cool stars from 1D-NLTE analysis of iron lines

For the well studied halo star HD 122563 and the four stars in the globular cluster NGC 6397, we determine NLTE abundances of iron using classical plane-parallel model atmospheres. Each star reveals a discrepancy in abundances between the Fe I lines arising from the ground state and the other Fe I lines, in qualitative agreement with the 3D-LTE line formation predictions, however, the magnitude of the observed effect is a factor of 2 smaller compared with the predicted one. When ignoring the Fe I low-excitation lines, the NLTE abundances from the two ionization stages, Fe I and Fe II, are consistent in each investigated star. For the subgiants in NGC 6397, this is only true when using the cooler effective temperature scale of Alonso et al. (1999). We also present full 3D-LTE line formation calculations for some selected iron lines in the solar and metal-poor 4480/2/-3 models and NLTE calculations with the corresponding spatial and temporal average <3D> models. The use of the <3D> models is justified only for particular Fe I lines in particular physical conditions. Our NLTE calculations reproduce well the centre-to-limb variation of the solar Fe I 7780 A line, but they are unsuccessful for Fe I 6151 A. The metal-poor <3D> model was found to be adequite for the strong Fe I 5166 A (Eexc = 0) line, but inadequite in all other investigated cases.

preprint2013arXiv

Spectral line asymmetries in the metal-poor red giant HD 122563: CO5BOLD predictions versus observations

We study the influence of convection on the asymmetries and Doppler shifts of Fe I spectral lines in the metal-poor red giant HD 122563. To this end, we compute theoretical Fe I line shifts and line bisectors using 3D hydrodynamical model atmosphere of HD 122563 calculated with the CO5BOLD code. We then make a detailed comparison of the theoretical line shifts and bisectors with those derived from the high quality HARPS spectrum of HD 122563 taken from the ESO Science Archive Facility (R = 115 000, average signal-to-noise ratio, S/N = 310). In general, we find a good agreement between the theoretically predicted and observed Doppler shifts of Fe I line cores, with somewhat larger discrepancies seen in the case of weaker (equivalent width W < 5 pm) and stronger lines (W > 11 pm). Both observed and theoretical coreshifts cover a range between 0 and -1 km/s, with increasingly stronger blueshifts for weaker lines and slight hints of a coreshift dependence on wavelength. Theoretical bisectors reproduce the observed ones reasonably well too, however, theoretical bisectors of the weak red (lambda > 600 nm) Fe I lines have blueshifts that are by up to 200 m/s larger than observed. The obtained results therefore suggest that the current CO5BOLD models are capable of reproducing the large-scale velocity fields in the atmosphere of HD 122563 suficiently well. Nevertheless, further efforts are needed in order to understand the physical reasons behind the discrepancies in theoretical predictions and observed properties of the weakest and strongest Fe I lines.

preprint2013arXiv

Spectroscopic analysis of DA white dwarfs with 3D model atmospheres

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.

preprint2013arXiv

Stellar granulation as seen in disk-integrated intensity. I. Simplified theoretical modeling

The solar granulation is known for a long time to be a surface manifestation of convection. Thanks to the current space-borne missions CoRoT and Kepler, it is now possible to observe in disk-integrated intensity the signature of this phenomena in a growing number of stars. The space-based photometric measurements show that the global brightness fluctuations and the lifetime associated with granulation obeys characteristic scaling relations. We thus aim at providing a simple theoretical modeling to reproduce these scaling relations and subsequently at inferring the physical properties of granulation properties across the HR diagram. We develop a simple 1D theoretical model that enable us to test any prescription concerning the time-correlation between granules. The input parameters of the model are extracted from 3D hydrodynamical models of the surface layers of stars, and the free parameters involved in the model are calibrated with solar observations. Two different prescriptions for representing the eddy time-correlation in the Fourier space are compared: a Lorentzian and an exponential form. Finally, we compare our theoretical prediction with a 3D radiative hydrodynamical (RHD) numerical modeling of stellar granulation (ab-initio approach). Provided that the free parameters are appropriately adjusted, our theoretical model satisfactorily reproduces the shape and the amplitude of the observed solar granulation spectrum. The best agreement is obtained with an exponential form. Furthermore, our theoretical model results in granulation spectra that consistently agree with the these calculated on the basis of the ab-initio approach with two 3D RHD models. Comparison between theoretical granulation spectra calculated with the present model and high precision photometry measurements of stellar granulation is undertaken in a companion paper.

preprint2013arXiv

Stellar granulation as seen in disk-integrated intensity. II. Theoretical scaling relations compared with observations

A large set of stars observed by CoRoT and Kepler shows clear evidence for the presence of a stellar background, which is interpreted to arise from surface convection, i.e., granulation. These observations show that the characteristic time-scale (tau_eff) and the root-mean-square (rms) brightness fluctuations (sigma) associated with the granulation scale as a function of the peak frequency (nu_max) of the solar-like oscillations. We aim at providing a theoretical background to the observed scaling relations based on a model developed in the companion paper. We computed for each 3D model the theoretical power density spectrum (PDS) associated with the granulation as seen in disk-integrated intensity on the basis of the theoretical model. For each PDS we derived tau_eff and sigma and compared these theoretical values with the theoretical scaling relations derived from the theoretical model and the Kepler measurements. We derive theoretical scaling relations for tau_eff and sigma, which show the same dependence on nu_max as the observed scaling relations. In addition, we show that these quantities also scale as a function of the turbulent Mach number (Ma) estimated at the photosphere. The theoretical scaling relations for tau_eff and sigma match the observations well on a global scale. Our modelling provides additional theoretical support for the observed variations of sigma and tau_eff with nu_m max. It also highlights the important role of Ma in controlling the properties of the stellar granulation. However, the observations made with Kepler on a wide variety of stars cannot confirm the dependence of our scaling relations on Ma. Measurements of the granulation background and detections of solar-like oscillations in a statistically sufficient number of cool dwarf stars will be required for confirming the dependence of the theoretical scaling relations with Ma.

preprint2013arXiv

The influence of convection on the atmospheric structures and observable properties of red giant stars

During the recent years significant progress has been made in the modeling of red giant atmospheres with the aid of 3D hydrodynamical model atmosphere codes. In this contribution we provide an overview of selected results obtained in this context by utilizing 3D hydrodynamical CO5BOLD stellar model atmospheres. Hydrodynamical simulations show that convective motions lead to significant differences in the atmospheric structures of red giants with respect to those predicted by the classical 1D model atmospheres. Results of these simulations also show that in certain cases 1D models fail to reproduce even the average properties of the 3D hydrodynamical models, such as P-T profiles. Large horizontal temperature fluctuations in the 3D model atmospheres, as well as differences between the temperature profiles of the average <3D> and 1D models, lead to large discrepancies in the strengths of spectral lines predicted by the 3D and 1D model atmospheres. This is especially important in models at lowest metallicities ([M/H]<-2.0) where the 3D-1D abundance differences may reach (or even exceed) -0.6 dex for lines of neutral atoms and molecules. We also discuss several simplifications and numerical aspects involved in the present 3D hydrodynamical modeling of red giant atmospheres, and briefly address several issues where urgent progress may be needed.

preprint2013arXiv

Three-dimensional hydrodynamical CO5BOLD model atmospheres of red giant stars. III. Line formation in the atmospheres of giants located close to the base of RGB

We utilize state-of-the-art 3D hydrodynamical and classical 1D stellar model atmospheres to study the influence of convection on the formation properties of various atomic and molecular spectral lines in the atmospheres of four red giant stars, located close to the base of the red giant branch, RGB ($T_{\mathrm eff}\approx5000$ K, $\log g=2.5$), and characterized by four different metallicities, [M/H] = 0.0, -1.0, -2.0, -3.0. The role of convection in the spectral line formation is assessed with the aid of abundance corrections, i.e., the differences in abundances predicted for a given equivalent width of a particular spectral line with the 3D and 1D model atmospheres. We find that for lines of certain neutral atoms the abundance corrections strongly depend both on metallicity of a given model atmosphere and the line excitation potential. While abundance corrections for all lines of both neutral and ionized elements tend to be small at solar metallicity, for lines of neutral elements with low ionization potential and low-to-intermediate $χ$ they quickly increase with decreasing metallicity, reaching in their extremes to -0.6...-0.8 dex. In all such cases the large abundance corrections are due to horizontal temperature fluctuations in the 3D hydrodynamical models. Abundance corrections of molecular lines are very sensitive to metallicity of the underlying model atmosphere and may be larger (in absolute value) than -0.5 dex at [M/H] = -3.0 (-1.5 dex in the case of CO). We also find that an approximate treatment of scattering in the 3D model calculations leads to the abundance corrections that are altered by less than ~0.1 dex, both for atomic and molecular (CO) lines, with respect to the model where scattering is treated as true absorption throughout the entire atmosphere, with the largest differences for the resonance and low-excitation lines.

preprint2013arXiv

Three-dimensional magnetohydrodynamic simulations of M-dwarf chromospheres

We present first results from three-dimensional radiation magnetohydrodynamic simulations of M-type dwarf stars with CO5BOLD. The local models include the top of the convection zone, the photosphere, and the chromosphere. The results are illustrated for models with an effective temperature of 3240 K and a gravitational acceleration of log g = 4.5, which represent analogues of AD Leo. The models have different initial magnetic field strengths and field topologies. This first generation of models demonstrates that the atmospheres of M-dwarfs are highly dynamic and intermittent. Magnetic fields and propagating shock waves produce a complicated fine-structure, which is clearly visible in synthetic intensity maps in the core of the Ca II K spectral line and also at millimeter wavelengths. The dynamic small-scale pattern cannot be described by means of one-dimensional models, which has important implications for the construction of semi-empirical model atmospheres and thus for the interpretation of observations in general. Detailed three-dimensional numerical simulations are valuable in this respect. Furthermore, such models facilitate the analysis of small-scale processes, which cannot be observed on stars but nevertheless might be essential for understanding M-dwarf atmospheres and their activity. An example are so-called "magnetic tornadoes", which have recently been found on the Sun and are presented here in M-dwarf models for the first time.

preprint2013arXiv

TOPoS: I. Survey design and analysis of the first sample

The metal-weak tail of the metallicity distribution function (MDF) of the Galactic Halo stars contains crucial information on the formation mode of the first generation of stars. To determine this observationally, it is necessary to observe large numbers of extremely metal-poor stars. We present here the Turn-Off Primordial Stars survey (TOPoS) that is conducted as an ESO Large Programme at the VLT. This project has {four} main goals: (i) to understand the formation of low-mass stars in a low-metallicity gas: determine the metal-weak tail of the halo MDF below [M/H]=-3.5. In particular, we aim at determining the critical metallicity, that is the lowest metallicity sufficient for the formation of low-mass stars; (ii) to determine the relative abundance of the elements in extremely metal-poor stars, which are the signature of the massive first stars; (iii) to determine the trend of the lithium abundance at the time when the Galaxy formed; and (iv) to derive the fraction of C-enhanced extremely metal-poor stars with respect to normal extremely metal-poor stars. The large number of stars observed in the SDSS provides a good sample of candidates of stars at extremely low metallicity. Candidates with turn-off colours down to magnitude g=20 were selected from the low-resolution spectra of SDSS by means of an automated procedure. X-Shooter has the potential of performing the necessary follow-up spectroscopy, providing accurate metallicities and abundance ratios for several key elements for these stars. We here present the stellar parameters of the first set of stars. The nineteen stars range in iron abundance between -4.1 and -2.9 dex relative to the Sun. Two stars have a high radial velocity and, according to our estimate of their kinematics, appear to be marginally bound to the Galaxy and are possibly accreted from another galaxy.

preprint2013arXiv

X-Shooter GTO: evidence for a population of extremely metal-poor, alpha-poor stars

The extremely metal-poor stars are the direct descendants of the first generation stars. They carry the chemical signature of the pristine Universe at the time they formed, shortly after the Big Bang. We aim to derive information about extremely metal-poor stars from their observed spectra. Four extremely metal-poor stars were selected from the Sloan Digital Sky Survey (SDSS) and observed during the guaranteed observing time of X-Shooter. The X-Shooter spectra were analysed using an automatic code, MyGIsFOS, which is based on a traditional analysis method. It makes use of a synthetic grid computed from one-dimensional, plane-parallel, hydrostatic model atmospheres. The low metallicity derived from the SDSS spectra is confirmed here. Two kinds of stars are found. Two stars are confirmed to be extremely metal-poor, with no evidence of any enhancement in carbon. The two other stars are strongly enhanced in carbon. We could not derive iron abundance for one of them, while [Ca/H] is below -4.5. Two of the stars are members of the rare population of extremely metal-poor stars low in alpha elements.

preprint2012arXiv

3D hydrodynamical CO5BOLD model atmospheres of red giant stars: I. Atmospheric structure of a giant located near the RGB tip

We investigate the character and role of convection in the atmosphere of a prototypical red giant located close to the red giant branch (RGB) tip with atmospheric parameters, Teff=3660K, log(g)=1.0, [M/H]=0.0. Differential analysis of the atmospheric structures is performed using the 3D hydrodynamical and 1D classical atmosphere models calculated with the CO5BOLD and LHD codes, respectively. All models share identical atmospheric parameters, elemental composition, opacities and equation-of-state. We find that the atmosphere of this particular red giant consists of two rather distinct regions: the lower atmosphere dominated by convective motions and the upper atmosphere dominated by wave activity. Convective motions form a prominent granulation pattern with an intensity contrast (~18%) which is larger than in the solar models (~15%). The upper atmosphere is frequently traversed by fast shock waves, with vertical and horizontal velocities of up to Mach ~2.5 and ~6.0, respectively. The typical diameter of the granules amounts to ~5Gm which translates into ~400 granules covering the whole stellar surface. The turbulent pressure in the giant model contributes up to ~35% to the total (i.e., gas plus turbulent) pressure which shows that it cannot be neglected in stellar atmosphere and evolutionary modeling. However, there exists no combination of the mixing-length parameter and turbulent pressure that would allow to satisfactorily reproduce the 3D temperature-pressure profile with 1D atmosphere models based on a standard formulation of mixing-length theory.

preprint2012arXiv

A primordial star in the heart of the Lion

Context: The discovery and chemical analysis of extremely metal-poor stars permit a better understanding of the star formation of the first generation of stars and of the Universe emerging from the Big Bang. aims: We report the study of a primordial star situated in the centre of the constellation Leo (SDSS J102915+172027). method: The star, selected from the low resolution-spectrum of the Sloan Digital Sky Survey, was observed at intermediate (with X-Shooter at VLT) and at high spectral resolution (with UVES at VLT). The stellar parameters were derived from the photometry. The standard spectroscopic analysis based on 1D ATLAS models was completed by applying 3D and non-LTE corrections. results: An iron abundance of [Fe/H]=--4.89 makes SDSS J102915+172927 one of the lowest [Fe/H] stars known. However, the absence of measurable C and N enhancements indicates that it has the lowest metallicity, Z<= 7.40x10^{-7} (metal-mass fraction), ever detected. No oxygen measurement was possible. conclusions: The discovery of SDSS J102915+172927 highlights that low-mass star formation occurred at metallicities lower than previously assumed. Even lower metallicity stars may yet be discovered, with a chemical composition closer to the composition of the primordial gas and of the first supernovae.

preprint2012arXiv

Amplitudes of solar-like oscillations in red-giant stars: Evidences for non-adiabatic effects using CoRoT observations

A growing number of solar-like oscillations has been detected in red giant stars thanks to CoRoT and Kepler space-crafts. The seismic data gathered by CoRoT on red giant stars allow us to test mode driving theory in physical conditions different from main-sequence stars. Using a set of 3D hydrodynamical models representative of the upper layers of sub- and red giant stars, we computed the acoustic mode energy supply rate (Pmax). Assuming adiabatic pulsations and using global stellar models that assume that the surface stratification comes from the 3D hydrodynamical models, we computed the mode amplitude in terms of surface velocity. This was converted into intensity fluctuations using either a simplified adiabatic scaling relation or a non-adiabatic one. From L and M (the luminosity and mass), the energy supply rate Pmax is found to scale as (L/M)^2.6 for both main-sequence and red giant stars, extending previous results. The theoretical amplitudes in velocity under-estimate the Doppler velocity measurements obtained so far from the ground for red giant stars by about 30%. In terms of intensity, the theoretical scaling law based on the adiabatic intensity-velocity scaling relation results in an under-estimation by a factor of about 2.5 with respect to the CoRoT seismic measurements. On the other hand, using the non-adiabatic intensity-velocity relation significantly reduces the discrepancy with the CoRoT data. The theoretical amplitudes remain 40% below, however, the CoRoT measurements. Our results show that scaling relations of mode amplitudes cannot be simply extended from main-sequence to red giant stars in terms of intensity on the basis of adiabatic relations because non-adiabatic effects for red giant stars are important and cannot be neglected. We discuss possible reasons for the remaining differences.

preprint2012arXiv

An extremely primitive halo star

The early Universe had a chemical composition consisting of hydrogen, helium and traces of lithium1, almost all other elements were created in stars and supernovae. The mass fraction, Z, of elements more massive than helium, is called "metallicity". A number of very metal poor stars have been found some of which, while having a low iron abundance, are rich in carbon, nitrogen and oxygen. For theoretical reasons and because of an observed absence of stars with metallicities lower than Z=1.5E-5, it has been suggested that low mass stars (M<0.8M\odot, the ones that survive to the present day) cannot form until the interstellar medium has been enriched above a critical value, estimated to lie in the range 1.5E-8\leqZ\leq1.5E-6, although competing theories claiming the contrary do exist. Here we report the chemical composition of a star with a very low Z\leq6.9E-7 (4.5E-5 of that of the Sun) and a chemical pattern typical of classical extremely metal poor stars, meaning without the enrichment of carbon, nitrogen and oxygen. This shows that low mass stars can be formed at very low metallicity. Lithium is not detected, suggesting a low metallicity extension of the previously observed trend in lithium depletion. Lithium depletion implies that the stellar material must have experienced temperatures above two million K in its history, which points to rather particular formation condition or internal mixing process, for low Z stars.

preprint2012arXiv

Barium abundance in red giants of NGC 6752. Non-local thermodynamic equilibrium and three-dimensional effects

(Abridged) Aims: We study the effects related to departures from non-local thermodynamic equilibrium (NLTE) and homogeneity in the atmospheres of red giant stars in Galactic globular cluster NGC 6752, to assess their influence on the formation of Ba II lines. Methods: One-dimensional (1D) local thermodynamic equilibrium (LTE) and 1D NLTE barium abundances were derived using classical 1D ATLAS stellar model atmospheres. The three-dimensional (3D) LTE abundances were obtained for 8 red giants on the lower RGB, by adjusting their 1D LTE abundances using 3D-1D abundance corrections, i.e., the differences between the abundances obtained from the same spectral line using the 3D hydrodynamical (CO5BOLD) and classical 1D (LHD) stellar model atmospheres. Results: The mean 1D barium-to-iron abundance ratios derived for 20 giants are <[Ba/Fe]>_{1D NLTE} = 0.05 \pm0.06 (stat.) \pm0.08 (sys.). The 3D-1D abundance correction obtained for 8 giants is small (~+0.05 dex), thus leads to only minor adjustment when applied to the mean 1D NLTE barium-to-iron abundance ratio for the 20 giants, <[Ba/Fe]>_{3D+NLTE} = 0.10 \pm0.06(stat.) \pm0.10(sys.). The intrinsic abundance spread between the individual cluster stars is small and can be explained in terms of uncertainties in the abundance determinations. Conclusions: Deviations from LTE play an important role in the formation of barium lines in the atmospheres of red giants studied here. The role of 3D hydrodynamical effects should not be dismissed either, even if the obtained 3D-1D abundance corrections are small. This result is a consequence of subtle fine-tuning of individual contributions from horizontal temperature fluctuations and differences between the average temperature profiles in the 3D and 1D model atmospheres: owing to the comparable size and opposite sign, their contributions nearly cancel each other.

preprint2012arXiv

Chemical abundances of distant extremely metal-poor unevolved stars

Aims: The purpose of our study is to determine the chemical composition of a sample of 16 candidate Extremely Metal-Poor (EMP) dwarf stars, extracted from the Sloan Digital Sky Survey (SDSS). There are two main purposes: in the first place to verify the reliability of the metallicity estimates derived from the SDSS spectra; in the second place to see if the abundance trends found for the brighter nearer stars studied previously also hold for this sample of fainter, more distant stars. Methods: We used the UVES at the VLT to obtain high-resolution spectra of the programme stars. The abundances were determined by an automatic analysis with the MyGIsFOS code, with the exception of lithium, for which the abundances were determined from the measured equivalent widths of the Li I resonance doublet. Results: All candidates are confirmed to be EMP stars, with [Fe/H]<= -3.0. The chemical composition of the sample of stars is similar to that of brighter and nearer samples. We measured the lithium abundance for 12 stars and provide stringent upper limits for three other stars, for a fourth star the upper limit is not significant, owing to the low signal-to noise ratio of the spectrum. The "meltdown" of the Spite plateau is confirmed, but some of the lowest metallicity stars of the sample lie on the plateau. Conclusions: The concordance of the metallicities derived from high-resolution spectra and those estimated from the SDSS spectra suggests that the latter may be used to study the metallicity distribution of the halo. The abundance pattern suggests that the halo was well mixed for all probed metallicities and distances. The fact that at the lowest metallicities we find stars on the Spite plateau suggests that the meltdown depends on at least another parameter, besides metallicity. (abridged)

preprint2012arXiv

Lithium and sodium in the globular cluster M4. Detection of a Li-rich dwarf star: preservation or pollution?

Context. The abundance inhomogeneities of light elements observed in Globular Clusters (GCs), and notably the ubiquitous Na-O anti-correlation, are generally interpreted as evidence that GCs comprise several generations of stars. There is an on-going debate as to the nature of the stars which produce the inhomogeneous elements, and investigating the behavior of several elements is a way to shed new light on this problem. Aims. We aim at investigating the Li and Na content of the GC M 4, that is known to have a well defined Na-O anti-correlation. Methods. We obtained moderate resolution (R=17 000-18 700) spectra for 91 main sequence (MS)/sub-giant branch stars of M 4 with the Giraffe spectrograph at the FLAMES/VLT ESO facility. Using model atmospheres analysis we measured lithium and sodium abundances. Results. We detect a weak Li-Na anti-correlation among un-evolved MS stars. One star in the sample, # 37934, shows the remarkably high lithium abundance A(Li)=2.87, compatible with current estimates of the primordial lithium abundance. Conclusions. The shallow slope found for the Li-Na anti-correlation suggests that lithium is produced in parallel to sodium. This evidence, coupled with its sodium-rich nature, suggests that the high lithium abundance of star # 37934 may originate by pollution from a previous generations of stars. The recent detection of a Li-rich dwarf of pollution origin in the globular cluster NGC 6397 may also point in this direction. Still, no clear cut evidence is available against a possible preservation of the primordial lithium abundance for star # 37934.

preprint2012arXiv

Simulations of the solar near-surface layers with the CO5BOLD, MURaM, and Stagger codes

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.

preprint2012arXiv

Spectroscopic and photometric studies of white dwarfs in the Hyades

The Hyades cluster is known to harbour ten so-called classical white dwarf members. Numerous studies through the years have predicted that more than twice this amount of degenerate stars should be associated with the cluster. Using the PPMXL catalog of proper motions and positions, a recent study proposed 17 new white dwarf candidates. We review the membership of these candidates by using published spectroscopic and photometric observations, as well as by simulating the contamination from field white dwarfs. In addition to the ten classical Hyades white dwarfs, we find six white dwarfs that may be of Hyades origin and three more objects that have an uncertain membership status due to their unknown or imprecise atmospheric parameters. Among those, two to three are expected as field stars contamination. Accurate radial velocity measurements will confirm or reject the candidates. One consequence is that the longstanding problem that no white dwarf older than ~340 Myr appears to be associated with the cluster remains unsolved.

preprint2012arXiv

Three-dimensional hydrodynamical CO5BOLD model atmospheres of red giant stars II. Spectral line formation in the atmosphere of a giant located near the RGB tip

We investigate the role of convection in the formation of atomic and molecular lines in the atmosphere of a red giant star. For this purpose we study the formation properties of spectral lines that belong to a number of astrophysically important tracer elements, including neutral and singly ionized atoms, and molecules. We focus our investigation on a prototypical red giant located close to the red giant branch (RGB) tip (Teff=3660K, logg=1.0, [M/H]=0.0). We used two types of model atmospheres, 3D hydrodynamical and classical 1D, calculated with the CO5BOLD and LHD stellar atmosphere codes, respectively. Both codes share the same atmospheric parameters, chemical composition, equation of state, and opacities, which allowed us to make a strictly differential comparison between the line formation properties predicted in 3D and 1D. The influence of convection on the spectral line formation was assessed with the aid of 3D-1D abundance corrections, which measure the difference between the abundances of chemical species derived with the 3D hydrodynamical and 1D classical model atmospheres. We find that convection plays a significant role in the spectral line formation in this particular red giant. The derived 3D-1D abundance corrections rarely exceed \pm0.1 dex when lines of neutral atoms and molecules are considered, which is in line with the previous findings for solar-metallicity red giants located on the lower RGB. The situation is different with lines that belong to ionized atoms, or to neutral atoms with high ionization potential. In both cases, the corrections for high-excitation lines (χ>8 eV) may amount to Δ_3D-1D ~ -0.4 dex. The 3D--1D abundance corrections generally show a significant wavelength dependence; in most cases they are smaller in the near-infrared, at 1600-2500 nm.

preprint2011arXiv

3D hydrodynamical CO5BOLD model atmospheres of late-type giants: chemical abundances from molecular lines

We investigate the influence of convection on the formation of molecular spectral lines in the atmospheres of late-type giants. For this purpose we use the 3D hydrodynamical CO5BOLD and classical 1D LHD stellar atmosphere codes and synthesize a number of fictitious lines belonging to a number of astrophysically relevant molecules, C2, CH, CN, CO, NH, OH. We find that differences between the abundances obtained from molecular lines using the 3D and 1D model atmospheres are generally small at [M/H]=0.0, but they quickly increase at sub-solar metallicities where for certain molecules they may reach -2.0 dex. The 3D-1D abundance differences show a significant dependence on the spectral line parameters, such as wavelength and excitation potential. Our comparison, therefore, reveals a complex interplay between the spectral line formation and convection that can not be properly accounted for with the classical 1D model atmospheres.

preprint2011arXiv

Chemical abundances in metal-poor giants: limitations imposed by the use of classical 1D stellar atmosphere models

In this work we have used 3D hydrodynamical (CO5BOLD) and 1D hydrostatic (LHD) stellar atmosphere models to study the importance of convection and horizontal temperature inhomogeneities in stellar abundance work related to late-type giants. We have found that for a number of key elements, such as Na, Mg, Si, Ca, Ti, Fe, Ni, Zn, Ba, Eu, differences in abundances predicted by 3D and 1D models are typically minor (< 0.1 dex) at solar metallicity. However, at [M/H] = -3 they become larger and reach to -0.5...-0.8 dex. In case of neutral atoms and fixed metallicity, the largest abundance differences were obtained for the spectral lines with lowest excitation potential, while for ionized species the largest 3D-1D abundance differences were found for lines of highest excitation potential. The large abundance differences at low metallicity are caused by large horizontal temperature fluctuations and lower mean temperature in the outer layers of the 3D hydrodynamical model compared with its 1D counterpart.

preprint2011arXiv

First Stars XIV. Sulfur abundances in extremely metal-poor (EMP) stars

Sulfur is important: the site of its formation is uncertain, and at very low metallicity the trend of [S/Fe] against [Fe/H] is controversial. Below [Fe/H]=-2.0, [S/Fe] remains constant or it decreases with [Fe/H], depending on the author and the multiplet used in the analysis. Moreover, although sulfur is not significantly bound in dust grains in the ISM, it seems to behave differently in DLAs and in old metal-poor stars. We aim to determine precise S abundance in a sample of extremely metal-poor stars taking into account NLTE and 3D effects. NLTE profiles of the lines of the multiplet 1 of SI have been computed using a new model atom for S. We find sulfur in EMP stars to behave like the other alpha-elements, with [S/Fe] remaining approximately constant for [Fe/H]<-3. However, [S/Mg] seems to decrease slightly as a function of [Mg/H]. The overall abundance patterns of O, Na, Mg, Al, S, and K are best matched by the SN model yields by Heger & Woosley. The [S/Zn] ratio in EMP stars is solar, as found also in DLAs. We obtain an upper limit on the abundance of sulfur, [S/Fe] < +0.5, for the ultra metal-poor star CS 22949-037. This, along with a previous reported measurement of zinc, argues against the conjecture that the light-element abundances pattern in this star, and, by analogy, the hyper metal-poor stars HE 0107-5240 and HE 1327-2326, are due to dust depletion.

preprint2011arXiv

Solution to the problem of the surface gravity distribution of cool DA white dwarfs from improved 3D model atmospheres

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.

preprint2011arXiv

X-Shooter GTO: Chemical analysis of a sample of EMP candidates

Context: Extremely metal-poor stars (EMP) are very rare objects that hold in their atmospheres the fossil record of the chemical composition of the early phases of Galactic evolution. Finding these objects and determining their chemical composition provides important constraints on these early phases. Aims: Using a carefully designed selection method, we chose a sample of candidate EMP stars from the low resolution spectra of the Sloan Digital Sky Survey and observed them with X-Shooter at the VLT to confirm their metallicities and determine abundances for as many elements as possible. Method: The X-Shooter spectra are analysed by means of one-dimensional, plane-parallel, hydrostatic model atmospheres. Corrections for the granulation effects are computed using CO5BOLD hydrodynamical simulations. Results: All the candidates are confirmed to be EMP stars, proving the efficiency of our selection method within about 0.5 dex. The chemical composition of this sample is compatible with those of brighter samples, suggesting that the stars in the Galactic halo are well mixed. Conclusions These observations show that it is feasible to observe, in a limited amount of time, a large sample of about one hundred stars among EMP candidates selected from the SDSS. Such a size of sample will allow us, in particular, to confirm or refute the existence of a vertical drop in the Galactic Halo metallicity distribution function around [Fe/H] ~ -3.5.

preprint2010arXiv

Asteroseismology of Solar-type Stars with Kepler I: Data Analysis

We report on the first asteroseismic analysis of solar-type stars observed by Kepler. Observations of three G-type stars, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation can clearly be distinguished in each star. We discuss the appearance of the oscillation spectra, including the presence of a possible signature of faculae, and the presence of mixed modes in one of the three stars.

preprint2010arXiv

Convection and 6Li in the atmospheres of metal-poor halo stars

Based on 3D hydrodynamical model atmospheres computed with the CO5BOLD code and 3D non-LTE (NLTE) line formation calculations, we study the effect of the convection-induced line asymmetry on the derived 6Li abundance for a range in effective temperature, gravity, and metallicity covering the stars of the Asplund et al. (2006) sample. When the asymmetry effect is taken into account for this sample of stars, the resulting 6Li/7Li ratios are reduced by about 1.5% on average with respect to the isotopic ratios determined by Asplund et al. (2006). This purely theoretical correction diminishes the number of significant 6Li detections from 9 to 4 (2 sigma criterion), or from 5 to 2 (3 sigma criterion). In view of this result the existence of a 6Li plateau appears questionable. A careful reanalysis of individual objects by fitting the observed lithium 6707 A doublet both with 3D NLTE and 1D LTE synthetic line profiles confirms that the inferred 6Li abundance is systematically lower when using 3D NLTE instead of 1D LTE line fitting. Nevertheless, halo stars with unquestionable 6Li detection do exist even if analyzed in 3D-NLTE, the most prominent example being HD 84937.

preprint2010arXiv

Galactic Evolution of Oxygen: OH lines in 3D hydrodynamical model atmospheres

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.

preprint2010arXiv

Gravitational redshifts in main-sequence and giant stars

Gravitational redshifts in solar-type main-sequence stars are expected to be some 500 ms$^{-1}$ greater than those in giants. Such a signature is searched for between groups of open-cluster stars which share the same average space motion and thus have the same average Doppler shift. 144 main-sequence stars and cool giants were observed in the M67 open cluster using the ESO FEROS spectrograph, obtaining radial velocities by cross correlation with a spectral template. M67 dwarf and giant radial-velocity distributions are well represented by Gaussian functions, sharing the same apparent average radial velocity within $\simeq$ 100 ms$^{-1}$. In addition, dwarfs in M67 appear to be dynamically hotter ($σ$ = 0.90 kms$^{-1}$) than giants ($σ$ = 0.68 kms$^{-1}$). Explanations for the lack of an expected signal are sought: a likely cause is the differential wavelength shifts produced by different hydrodynamics in dwarf and giant atmospheres. Radial-velocity differences measured between unblended lines in low-noise averaged spectra vary with line-strength: stronger lines are more blushifted in dwarfs than in giants, apparently 'compensating' for the gravitational redshift. Synthetic high-resolution spectra are computed from 3-dimensional hydrodynamic model atmospheres for both giants and dwarfs, and synthetic wavelength shifts obtained. In agreement with observations, 3D models predict substantially smaller wavelength-shift differences than expected from gravitational redshift only. The procedures developed could be used to test 3D models for different classes of stars, but will ultimately require high-fidelity spectra for measurements of wavelength shifts in individual spectral lines.

preprint2010arXiv

Radiative hydrodynamics simulations of red supergiant stars. III. Spectro-photocentric variability, photometric variability, and consequences on Gaia measurements

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.

preprint2010arXiv

The asteroseismic potential of Kepler: first results for solar-type stars

We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies and frequency separations to provide first results on the radii, masses and ages of the stars, and comment in the light of these results on prospects for inference on other solar-type stars that Kepler will observe.

preprint2010arXiv

The metal-poor end of the Spite plateau. 1: Stellar parameters, metallicities and lithium abundances

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.

preprint2010arXiv

The solar photospheric abundance of carbon.Analysis of atomic carbon lines with the CO5BOLD solar model

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)

preprint2010arXiv

Three carbon-enhanced metal-poor dwarf stars from the SDSS - Chemical abundances from CO^5BOLD 3D hydrodynamical model atmospheres

The origin of carbon-enhanced metal-poor stars enriched with both s and r elements is highly debated. Detailed abundances of these types of stars are crucial to understand the nature of their progenitors. The aim of this investigation is to study in detail the abundances of SDSS J1349-0229, SDSS J0912+0216 and SDSS J1036+1212, three dwarf CEMP stars, selected from the Sloan Digital Sky Survey. Using high resolution VLT/UVES spectra (R ~ 30 000) we determine abundances for Li, C, N, O, Na, Mg, Al, Ca, Sc, Ti, Cr, Mn, Fe, Co, Ni and 21 neutron-capture elements. We made use of CO^5BOLD 3D hydrodynamical model atmospheres in the analysis of the carbon, nitrogen and oxygen abundances. NLTE corrections for C I and O I lines were computed using the Kiel code. We classify SDSS J1349-0229 and SDSS J0912+0216 as CEMP-r+s stars. SDSS J1036+1212 belongs to the class CEMP-no/s, with enhanced Ba, but deficient Sr, of which it is the third member discovered to date. Radial-velocity variations have been observed in SDSS J1349-0229, providing evidence that it is a member of a binary system. The chemical composition of the three stars is generally compatible with mass transfer from an AGB companion. However, many details remain difficult to explain. Most notably of those are the abundance of Li at the level of the Spite plateau in SDSS J1036+1212 and the large over-abundance of the pure r-process element Eu in all three stars.

preprint2009arXiv

3D simulations of M star atmosphere velocities and their influence on molecular FeH lines

We present an investigation of the velocity fields in early to late M-type star hydrodynamic models, and we simulate their influence on FeH molecular line shapes. The M star model parameters range between log g of 3.0 - 5.0 and Teff of 2500 K and 4000 K. Our aim is to characterize the Teff- and log g -dependence of the velocity fields and express them in terms of micro- and macro-turbulent velocities in the one dimensional sense. We present also a direct comparison between 3D hydrodynamical velocity fields and 1D turbulent velocities. The velocity fields strongly affect the line shapes of FeH, and it is our goal to give a rough estimate for the log g and Teff parameter range in which 3D spectral synthesis is necessary and where 1D synthesis suffices. In order to calculate M-star structure models we employ the 3D radiative-hydrodynamics (RHD) code CO5BOLD. The spectral synthesis on these models is performed with the line synthesis code LINFOR3D. We describe the 3D velocity fields in terms of a Gaussian standard deviation and project them onto the line of sight to include geometrical and limb-darkening effects. The micro- and macro-turbulent velocities are determined with the "Curve of Growth" method and convolution with a Gaussian velocity profile, respectively. To characterize the log g and Teff dependence of FeH lines, the equivalent width, line width, and line depth are regarded. The velocity fields in M-stars strongly depend on log g and Teff. They become stronger with decreasing log g and increasing Teff.

preprint2009arXiv

6Li in metal-poor halo stars: real or spurious?

The presence of convective motions in the atmospheres of metal-poor halo stars leads to systematic asymmetries of the emergent spectral line profiles. Since such line asymmetries are very small, they can be safely ignored for standard spectroscopic abundance analysis. However, when it comes to the determination of the 6Li/7Li isotopic ratio, q(Li)=n(6Li)/n(7Li), the intrinsic asymmetry of the 7Li line must be taken into account, because its signature is essentially indistinguishable from the presence of a weak 6Li blend in the red wing of the 7Li line. In this contribution we quantity the error of the inferred 6Li/7Li isotopic ratio that arises if the convective line asymmetry is ignored in the fitting of the lithium blend at 6707 A. Our conclusion is that 6Li/7Li ratios derived by Asplund et al. (2006), using symmetric line profiles, must be reduced by typically dq(Li) ~ 0.015. This diminishes the number of certain 6Li detections from 9 to 4 stars or less, casting some doubt on the existence of a 6Li plateau.

preprint2009arXiv

Accuracy of spectroscopy-based radioactive dating of stars

Combined spectroscopic abundance analyses of stable and radioactive elements can be applied for deriving stellar ages. The achievable precision depends on factors related to spectroscopy, nucleosynthesis, and chemical evolution. We quantify the uncertainties arising from the spectroscopic analysis, and compare these to the other error sources. We derive formulae for the age uncertainties arising from the spectroscopic abundance analysis, and apply them to spectroscopic and nucleosynthetic data compiled from the literature for the Sun and metal-poor stars. We obtained ready-to-use analytic formulae of the age uncertainty for the cases of stable+unstable and unstable+unstable chronometer pairs, and discuss the optimal relation between to-be-measured age and mean lifetime of a radioactive species. Application to the literature data indicates that, for a single star, the achievable spectroscopic accuracy is limited to about +/- 20% for the foreseeable future. At present, theoretical uncertainties in nucleosynthesis and chemical evolution models form the precision bottleneck. For stellar clusters, isochrone fitting provides a higher accuracy than radioactive dating, but radioactive dating becomes competitive when applied to many cluster members simultaneously, reducing the statistical errors by a factor sqrt(N). Spectroscopy-based radioactive stellar dating would benefit from improvements in the theoretical understanding of nucleosynthesis and chemical evolution. Its application to clusters can provide strong constraints for nucleosynthetic models.

preprint2009arXiv

Can we trust elemental abundances derived in late-type giants with the classical 1D stellar atmosphere models?

We compare the abundances of various chemical species as derived with 3D hydrodynamical and classical 1D stellar atmosphere codes in a late-type giant characterized by T_eff=3640K, log g = 1.0, [M/H] = 0.0. For this particular set of atmospheric parameters the 3D-1D abundance differences are generally small for neutral atoms and molecules but they may reach up to 0.3-0.4 dex in case of ions. The 3D-1D differences generally become increasingly more negative at higher excitation potentials and are typically largest in the optical wavelength range. Their sign can be both positive and negative, and depends on the excitation potential and wavelength of a given spectral line. While our results obtained with this particular late-type giant model suggest that 1D stellar atmosphere models may be safe to use with neutral atoms and molecules, care should be taken if they are exploited with ions.

preprint2009arXiv

Detailed analyses of three neutron-capture-rich carbon-enhanced metal-poor stars

Approximately 20% of very metal-poor stars ([Fe/H] < -2.0) are strongly enhanced in carbon ([C/Fe] > +1.0). Such stars are referred to as carbon-enhanced metal-poor (CEMP) stars. We present a chemical abundance analysis based on high resolution spectra acquired with UVES at the VLT of three dwarf CEMP stars: SDSS J1349-0229, SDSS J0912+0216 and SDSS J1036+1212. These very metal-poor stars, with [Fe/H] < -2.5, were selected from our ongoing survey of extremely metal-poor dwarf candidates from the SDSS. Among these CEMPs, SDSS J1349-0229 has been identified as a carbon star ([C/O] > +1.0). First and second peak s-process elements, as well as second peak r-process elements have been detected in all stars. In addition, elements from the third r-process peak were detected in one of the stars, SDSS J1036+1212. We present the abundance results of these stars in the context of neutron-capture nucleosynthesis theories.

preprint2009arXiv

Hydrodynamical simulations of convection-related stellar micro-variability. II. The enigmatic granulation background of the COROT target HD49933

Local-box hydrodynamical model atmospheres provide statistical information about a star's emergent radiation field which allows one to predict the level of its granulation-related micro-variability. Space-based photometry is now sufficiently accurate to test model predictions. We aim to model the photometric granulation background of HD49933 as well as the Sun, and compare the predictions to the measurements obtained by the COROT and SOHO satellite missions. We construct hydrodynamical model atmospheres representing HD49933 and the Sun, and use a previously developed scaling technique to obtain the observable disk-integrated brightness fluctuations. We further performed exploratory magneto-hydrodynamical simulations to gauge the impact of small scale magnetic fields on the synthetic light-curves. We find that the granulation-related brightness fluctuations depend on metallicity. We obtain a satisfactory correspondence between prediction and observation for the Sun, validating our approach. For HD49933, we arrive at a significant over-estimation by a factor of two to three in total power. Locally generated magnetic fields are unlikely to be responsible, otherwise existing fields would need to be rather strong to sufficiently suppress the granulation signal. Presently suggested updates on the fundamental stellar parameters do not improve the correspondence; however, an ad-hoc increase of the HD49933 surface gravity by about 0.2dex would eliminate most of the discrepancy. We diagnose a puzzling discrepancy between the predicted and observed granulation background in HD49933, with only rather ad-hoc ideas for remedies at hand.

preprint2009arXiv

Impact of granulation effects on the use of Balmer lines as temperature indicators

Balmer lines serve as important indicators of stellar effective temperatures in late-type stellar spectra. One of their modelling uncertainties is the influence of convective flows on their shape. We aim to characterize the influence of convection on the wings of Balmer lines. We perform a differential comparison of synthetic Balmer line profiles obtained from 3D hydrodynamical model atmospheres and 1D hydrostatic standard ones. The model parameters are appropriate for F,G,K dwarf and subgiant stars of metallicity ranging from solar to 1/1000 solar. The shape of the Balmer lines predicted by 3D models can never be exactly reproduced by a 1D model, irrespective of its effective temperature. We introduce the concept of a 3D temperature correction, as the effective temperature difference between a 3D model and a 1D model which provides the closest match to the 3D profile. The temperature correction is different for the different members of the Balmer series and depends on the adopted mixing-length parameter in the 1D model. Among the investigated models, the 3D correction ranges from -300K to +300K. Horizontal temperature fluctuations tend to reduce the 3D correction. Accurate effective temperatures cannot be derived from the wings of Balmer lines, unless the effects of convection are properly accounted for. The 3D models offer a physically well justified way of doing so. The use of 1D models treating convection with the mixing-length theory do not appear to be suitable for this purpose. In particular, there are indications that it is not possible to determine a single value of the mixing-length parameter which will optimally reproduce the Balmer lines for any choice of atmospheric parameters.

preprint2009arXiv

Lithium abundances of Main Sequence and Sub-Giant stars in the Globular Cluster NGC 6397

We present FLAMES/GIRAFFE spectroscopy obtained at the Very Large Telescope (VLT). Using these observations we have been able for the first time to observe the Li I doublet in the Main Sequence stars of a Globular Cluster. We also observed Li in a sample of Sub-Giant stars of the same B-V colour. Our final sample is composed of 84 SG stars and 79 MS stars. In spite of the fact that SG and MS span the same temperature range we find that the equivalent widths of the Li I doublet in SG stars are systematically larger than those in MS stars, suggesting a higher Li content among SG stars. This is confirmed by our quantitative analysis which makes use of both 1D and 3D model atmospheres. We find that SG stars show, on average, a Li abundance higher by 0.1 dex than MS stars. We also detect a positive slope of Li abundance with effective temperature, the higher the temperature the higher the Li abundance, both for SG and MS stars, although the slope is slightly steeper for MS stars. These results provide an unambigous evidence that the Li abundance changes with evolutionary status. The physical mechanisms that contribute to this are not yet clear, since none of the proposed models seem to describe accurately the observations. Whether such mechanism can explain the cosmological lithium problem, is still an open question.

preprint2009arXiv

Lithium in the Globular Cluster NGC 6397: Evidence for dependence on evolutionary status

Most Globular Clusters are believed to host a single stellar populations. They can thus be considered a good place to study the Spite plateau and probe for possible evolutionary modifications of the Li content. We want to determine the Li content of subgiant (SG) and Main Sequence (MS) stars of the old, metal-poor globular cluster NGC 6397. This work was aimed not only at studying possible Li abundance variations but to investigate the cosmological Li discrepancy. Here, we present FLAMES/GIRAFFE observations of a sample of 84 SG and 79 MS stars in NGC 6397 selected in a narrow range of B-V colour and, therefore, effective temperatures. We determine both Teff and A(Li) using 3D hydrodynamical model atmospheres for all the MS and SG stars of the sample. We find a significant difference in the Li abundance between SG stars and MS stars, the SG stars having an A(Li) higher by almost 0.1 dex on average. We also find a decrease in the A(Li) with decreasing Teff, both in MS and SG stars, albeit with a significantly different slope for the two classes of stars. This suggests that the lithium abundance in these stars is, indeed, altered by some process, which is Teff-dependent. The Li abundance pattern observed in NGC 6397 is different from what is found among field stars, casting some doubt on the use of Globular Cluster stars as representative of Population II with respect to the Li abundance. None of the available theories of Li depletion appears to satisfactorily describe our observations.

preprint2009arXiv

Main-Sequence and sub-giant stars in the Globular Cluster NGC6397: The complex evolution of the lithium abundance

Thanks to the high multiplex and efficiency of Giraffe at the VLT we have been able for the first time to observe the Li I doublet in the Main Sequence (MS) stars of a Globular Cluster. At the same time we observed Li in a sample of Sub-Giant (SG) stars of the same B-V colour. Our final sample is composed of 84 SG stars and 79 MS stars. In spite of the fact that SG and MS span the same temperature range we find that the equivalent widths of the Li I doublet in SG stars are systematically larger than those in MS stars, suggesting a higher Li content among SG stars. This is confirmed by our quantitative analysis. We derived the effective temperatures, from H$α$ fitting, and NLTE Li abundances of the stars in our the sample, using 3D and 1D models. We find that SG stars have a mean Li abundance higher by 0.1dex than MS stars, using both 1D and 3D models. We also detect a positive slope of Li abundance with effective temperature. These results provide an unambiguous evidence that the Li abundance changes with evolutionary status. The physical mechanisms responsible for this behaviour are not yet clear, and none of the existing models seems to describe accurately these observations. Based on these conclusions, we believe that the cosmological lithium problem still remains an open question.

preprint2009arXiv

Solar abundances and 3D model atmospheres

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.

preprint2009arXiv

Sulphur in the Globular Clusters 47 Tucanae and NGC 6752

The light elements Li, O, Na, Al, and Mg are known to show star to star variations in the globular clusters 47 Tuc and NGC 6752. Such variations are interpreted as due to processing in a previous generation of stars. In this paper we investigate the abundances of the alpha-element sulphur, for which no previous measurements exist. In fact this element has not been investigated in any Galactic globular cluster so far. The only globular cluster for which such measurements are available is Terzan 7, which belongs to the Sgr dSph. We use high resolution spectra of the S I Mult. 1, acquired with the UVES spectrograph at the 8.2m VLT-Kueyen telescope, for turn-off and giant stars in the two globular clusters. The spectra are analyzed making use of ATLAS static plane parallel model atmospheres and SYNTHE spectrum synthesis. We also compute 3D corrections from CO5BOLD hydrodynamic models and apply corrections due to NLTE effects taken from the literature. In the cluster NGC 6752 sulphur has been measured only in four subgiant stars. We find no significant star to star scatter and a mean <[S/Fe]> = +0.49 +- 0.15, consistent with what observed in field stars of the same metallicity. In the cluster 47 Tuc we measured S in 4 turn-off and 5 subgiant stars with a mean <[S/Fe]> = +0.18 +- 0.14. While this result is compatible with no star to star scatter we notice a statistically significant correlation of the sulphur abundance with the sodium abundance and a tentative correlation with the silicon abundance.

preprint2009arXiv

Teff and log g dependence of velocity fields in M-stars

We present an investigation of velocity fields in early to late M-type hydrodynamic stellar atmosphere models. These velocities will be expressed in classical terms of micro- and macro-turbulent velocities for usage in 1D spectral synthesis. The M-star model parameters range between log g of 3.0 - 5.0 and Teff of 2500 K - 4000 K. We characterize the Teff- and log g-dependence of the hydrodynamical velocity fields in these models with a binning method, and for the determination of micro-turbulent velocities, the Curve of Growth method is used. The macro-turbulent velocities are obtained by convolutions with Gaussian profiles. Velocity fields in M-stars strongly depend on log g and Teff. Their velocity amplitudes increase with decreasing log g and increasing Teff. The 3D hydrodynamical and 1D macro-turbulent velocities range from ~100 m/s for cool high gravity models to ~ 800 m/s - 1000 m/s for hot models or models with low log g. The micro-turbulent velocities range in the order of ~100 m/s for cool models, to ~600 m/s for hot or low log g models. Our M-star structure models are calculated with the 3D radiative-hydrodynamics (RHD) code CO5BOLD. The spectral synthesis on these models is performed with the line synthesis code LINFOR3D.

preprint2009arXiv

The CoRoT target HD 49933: 1- Role of the metal abundance

Solar-like oscillations are stochastically excited by turbulent convection at the surface layers of the stars. We study the role of the surface metal abundance on the efficiency of the stochastic driving in the case of the CoRoT target HD 49933. We compute two 3D hydrodynamical simulations representative -- in effective temperature and gravity -- of the surface layers of the CoRoT target HD 49933, a star that is rather metal poor and significantly hotter compared to the Sun. One 3D simulation has a solar metal abundance and the other has a surface iron-to-hydrogen, [Fe/H], abundance ten times smaller. For each 3D simulation we match an associated global 1D model and we compute the associated acoustic modes using a theoretical model of stochastic excitation validated in the case of the Sun and Alpha Cen A. The rate at which energy is supplied per unit time into the acoustic modes associated with the 3D simulation with [Fe/H]=-1 are found about three times smaller than those associated with the 3D simulation with [Fe/H]=0. As shown here, these differences are related to the fact that low metallicity implies surface layers with a higher mean density. In turn, a higher mean density favors smaller convective velocities and hence less efficient driving of the acoustic modes. Our result shows the importance of taking the surface metal abundance into account in the modeling of the mode driving by turbulent convection. A comparison with observational data is presented in a companion paper using seismic data obtained for the CoRoT target HD 49933.

preprint2009arXiv

The CoRoT target HD 49933: 2- Comparison of theoretical mode amplitudes with observations

From the seismic data obtained by CoRoT for the star HD 49933 it is possible, as for the Sun, to constrain models of the excitation of acoustic modes by turbulent convection. We compare a stochastic excitation model described in Paper I (arXiv:0910.4027) with the asteroseismology data for HD 49933, a star that is rather metal poor and significantly hotter than the Sun. Using the mode linewidths measured by CoRoT for HD 49933 and the theoretical mode excitation rates computed in Paper I, we derive the expected surface velocity amplitudes of the acoustic modes detected in HD 49933. Using a calibrated quasi-adiabatic approximation relating the mode amplitudes in intensity to those in velocity, we derive the expected values of the mode amplitude in intensity. Our amplitude calculations are within 1-sigma error bars of the mode surface velocity spectrum derived with the HARPS spectrograph. The same is found with the mode amplitudes in intensity derived for HD 49933 from the CoRoT data. On the other hand, at high frequency, our calculations significantly depart from the CoRoT and HARPS measurements. We show that assuming a solar metal abundance rather than the actual metal abundance of the star would result in a larger discrepancy with the seismic data. Furthermore, calculations that assume the ``new'' solar chemical mixture are in better agreement with the seismic data than those that assume the ``old'' solar chemical mixture. These results validate, in the case of a star significantly hotter than the Sun and Alpha Cen A, the main assumptions in the model of stochastic excitation. However, the discrepancies seen at high frequency highlight some deficiencies of the modelling, whose origin remains to be understood.

preprint2009arXiv

The Solar Photospheric Nitrogen Abundance: Determination with 3D and 1D Model Atmospheres

We present a new determination of the solar nitrogen abundance making use of 3D hydrodynamical modelling of the solar photosphere, which is more physically motivated than traditional static 1D models. We selected suitable atomic spectral lines, relying on equivalent width measurements already existing in the literature. For atmospheric modelling we used the co 5 bold 3D radiation hydrodynamics code. We investigated the influence of both deviations from local thermodynamic equilibrium (non-LTE effects) and photospheric inhomogeneities (granulation effects) on the resulting abundance. We also compared several atlases of solar flux and centre-disc intensity presently available. As a result of our analysis, the photospheric solar nitrogen abundance is A(N) = 7.86 +/- 0.12.

preprint2009arXiv

The solar photospheric nitrogen abundance. Analysis of atomic transitions with 3D and 1D model atmospheres

CONTEXT: In recent years, the solar chemical abundances have been studied in considerable detail because of discrepant values of solar metallicity inferred from different indicators, i.e., on the one hand, the "sub-solar" photospheric abundances resulting from spectroscopic chemical composition analyses with the aid of 3D hydrodynamical models of the solar atmosphere, and, on the other hand, the high metallicity inferred by helioseismology. AIMS: After investigating the solar oxygen abundance using a CO5BOLD 3D hydrodynamical solar model in previous work, we undertake a similar approach studying the solar abundance of nitrogen, since this element accounts for a significant fraction of the overall solar metallicity, Z. METHOD: We used a selection of atomic spectral lines to determine the solar nitrogen abundance, relying mainly on equivalent width measurements in the literature. We investigate the influence on the abundance analysis, of both deviations from local thermodynamic equilibrium ("NLTE effects") and photospheric inhomogeneities ("granulation effects"). RESULTS: We recommend use of a solar nitrogen abundance of A(N)=7.86+-0.12 whose error bar reflects the line-to-line scatter. CONCLUSION: The solar metallicity implied by the CO5BOLD-based nitrogen and oxygen abundances is in the range 0.0145<= Z <= 0.0167. This result is a step towards reconciling photospheric abundances with helioseismic constraints on Z. Our most suitable estimates are Z=0.0156 and Z/X=0.0213.

preprint2009arXiv

Theoretical amplitudes and lifetimes of non-radial solar-like oscillations in red giants

Solar-like oscillations have been observed in numerous red giants from ground and from space. An important question arises: could we expect to detect non-radial modes probing the internal structure of these stars? We investigate under what physical circumstances non-radial modes could be observable in red giants; what would be their amplitudes, lifetimes and heights in the power spectrum (PS)? Using a non-radial non-adiabatic pulsation code including a non-local time-dependent treatment of convection, we compute the theoretical lifetimes of radial and non-radial modes in several red giant models. Next, using a stochastic excitation model, we compute the amplitudes of these modes and their heights in the PS. Very distinct results are found depending on the evolutionary status of the star.

preprint2002arXiv

A simulation of solar convection at supergranulation scale

We present here numerical simulations of surface solar convection which cover a box of 30$\times30\times$3.2 Mm$^3$ with a resolution of 315$\times315\times$82, which is used to investigate the dynamics of scales larger than granulation. No structure resembling supergranulation is present; possibly higher Reynolds numbers (i.e. higher numerical resolution), or magnetic fields, or greater depth are necessary. The results also show interesting aspects of granular dynamics which are briefly presented, like extensive p-mode ridges in the k-$ω$ diagram and a ringlike distribution of horizontal vorticity around granules. At large scales, the horizontal velocity is much larger than the vertical velocity and the vertical motion is dominated by p-mode oscillations.