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

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 spectra. It is demonstrated that our 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 spectra and that the issue was caused by inaccuracies in the 1D mixing-length approach. The white dwarfs with a radiative and a convective atmosphere have derived mean masses that are the same within ~0.01 Msun, in much better agreement with our understanding of stellar evolution. Furthermore, the 3D atmospheric parameters are in better agreement with independent Teff and log g values from photometric and parallax measurements.Comment: 15 pages, 18 figures, 10 pages online appendix, accepted for publication in Astronomy and Astrophysic

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

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

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

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

Stellar Envelope Convection calibrated by Radiation Hydrodynamics Simulations: Influence on Globular Clusters Isochrones

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