B. Freytag

Solar Chemical Abundances Determined with a CO5BOLD 3D Model Atmosphere

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

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

Context. The primordial nature of the Spite plateau is at odds with the WMAP satellite measurements, implying a primordial Li production at least three times higher than observed. It has also been suggested that A(Li) might exhibit a positive correlation with metallicity below [Fe/H] ~ -2.5. Previous samples studied comprised few stars below [Fe/H] = -3. Aims. We present VLT-UVES Li abundances of 28 halo dwarf stars between [Fe/H] = -2.5 and -3.5, ten of which have [Fe/H] < -3. Methods. We determined stellar parameters and abundances using four different T eff scales. The direct infrared flux method was applied to infrared photometry. Hα wings were fitted with two synthetic grids computed by means of 1D LTE atmosphere models, assuming two different self-broadening theories. A grid of Hα profiles was finally 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. Results. We confirm previous claims that A(Li) does not exhibit a plateau below [Fe/H] = -3. We detect a strong positive correlation with [Fe/H] that is insensitive to the choice of T eff estimator. From a linear fit, we infer a steep slope of about 0.30 dex in A(Li) per dex in [Fe/H], which has a significance of 2-3σ. The slopes derived using the four T eff estimators are consistent to within 1σ. A significant slope is also detected in the A(Li)-T eff plane, driven mainly by the coolest stars in the sample (T eff < 6250), which appear to be Li-poor. However, when we remove these stars the slope detected in the A(Li)-[Fe/H] plane is not altered significantly. When the full sample is considered, 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. At this metallicity, the plateau lies at 〈A(Li) 3D,NLTE 〉 = 2.199 ± 0.086. Conclusions. The meltdown of the Spite plateau below [Fe/H] ~ -3 is established, but its cause is unclear. If the primordial A(Li) were that derived from standard BBN, 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. That no star below [Fe/H] = -3 lies above the plateau suggests that they formed at plateau level and experienced subsequent depletion.

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 CO 5 BOLD 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 ... ... ... ... ... ...

The role of convection, overshoot, and gravity waves for the transport of dust in M dwarf and brown dwarf atmospheres

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

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

Context. 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. Aims. We quantitatively compare results from three-dimensional, radiative hydrodynamic simulations of convection near the solar surface generated with three numerical codes (CO 5 BOLD, MURaM, and Stagger) and different simulation setups in order to investigate the level of similarity and to cross-validate the simulations. Methods. 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. Results. 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. Conclusions. The results give confidence in the reliability of the results from comprehensive radiative hydrodynamic simulations.

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

Context. The analysis of the solar spectra using hydrodynamical simulations, with a specific selection of lines, atomic data, and method for computing deviations from local thermodynamical equilibrium, has led to a downward revision of the solar metallicity, Z. We are using the latest simulations computed with the CO5BOLD code to reassess the solar chemical composition. Our previous analyses of the key elements, oxygen and nitrogen, have not confirmed any extreme downward revision of Z, as derived in other works based on hydrodynamical models. Aims. We determine the solar photospheric carbon abundance with a radiation-hydrodynamical CO5BOLD model and compute the departures from local thermodynamical equilibrium by using the Kiel code. Methods. We measured equivalent widths of atomic C i lines on high-resolution, high signal-to-noise ratio solar atlases of disccentre intensity and integrated disc flux. These equivalent widths were analysed with our latest solar 3D hydrodynamical simulation computed with CO5BOLD. Deviations from local thermodynamic equilibrium we computed in 1D with the Kiel code, using the average temperature structure of the hydrodynamical simulation as a background model. Results. 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 errors. 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. Conclusions. Our analysis implies a solar carbon abundance that is about 0.1 dex higher than what was found in previous analyses based on different 3D hydrodynamical computations. The difference is partly driven by our equivalent width measurements (we measure, on average, larger equivalent widths than the other work based on a 3D model), in part because of 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.

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

The surface gravities of cool (T eff < 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 < T eff < 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 ID 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.

Atmospheres from very low-mass stars to extrasolar planets

Within the next few years, several instruments aiming at imaging extrasolar planets will see first light. In parallel, low mass planets are being searched around red dwarfs which offer more favorable conditions, both for radial velocity detection and transit studies, than solar-type stars. We review recent advancements in modeling the stellar to substellar transition. The revised solar oxygen abundances and cloud models allow to reproduce the photometric and spectroscopic properties of this transition to a degree never achieved before, but problems remain in the important M-L transition characteristic of the effective temperature range of characterizable exoplanets.

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.

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.