Accretion from debris disks onto white dwarfs : Fingering (thermohaline) instability and derived accretion rates

Recent observations of a large number of DA and DB white dwarfs show evidence of debris disks, which are the remnants of old planetary systems. The infrared excess detected with \emph{Spitzer} and the lines of heavy elements observed in their atmospheres with high-resolution spectroscopy converge on the idea that planetary material accretes onto these stars. Accretion rates have been derived by several authors with the assumption of a steady state between accretion and gravitational settling. The results are unrealistically different for DA and DB white dwarfs. When heavy matter is accreted onto stars, it induces an inverse $\mu$-gradient that leads to fingering (thermohaline) convection. The aim of this letter is to study the impact of this specific process on the derived accretion rates in white dwarfs and on the difference between DA and DB. We solve the diffusion equation for the accreted heavy elements with a time-dependent method. The models we use have been obtained both with the IRAP code, which computes static models, and the La Plata code, which computes evolutionary sequences. Computations with pure gravitational settling are compared with computations that include fingering convection. The most important result is that fingering convection has very important effects on DAs but is inefficient in DBs. When only gravitational settling is taken into account, the time-dependent computations lead to a steady state, as postulated by previous authors. When fingering convection is added, this steady state occurs much later. The surprising difference found in the past for the accretion rates derived for DA and DB white dwarfs disappears. The derived accretion rates for DAs are increased when fingering convection is taken into account, whereas those for DBs are not modified. More precise and developed results will be given in a forthcoming paper

On the possible existence of a self-regulating hydrodynamical process in slowly rotating stars II. Lithium Plateau in Halo Stars and Primordial Abundance

The lithium plateau observed in halo stars has long appeared as a paradox in the general context of the lithium abundance behavior in stellar outer layers. First, the plateau is flat, second, the lithium abundance dispersion is extremely small. This seems in contradiction with the large lithium variations observed in younger stars. It is also difficult to understand theoretically as lithium nuclei are easily destroyed by nuclear reactions, and subject to microscopic diffusion which, in the case of halo stars, should also lead to depletion. Several ideas have been proposed to account for the lithium behavior in halo stars. The most promising possibilities were rotational-induced mixing, which could reduce lithium in the same way for all the stars, and mass-loss, which could oppose the lithium settling. In both cases however, the parameters should be tightly adjusted to prevent any dispersion in the final results. Vauclair 1999 (paper I) pointed out that the mu-gradient terms which appear in the computations of the meridional circulation velocity were not introduced in previous computations of rotationally-induced mixing. This can lead to a self-regulating process which reduces the efficiency of the meridional circulation as well as the microscopic diffusion. Here we present numerical computations of this process and its influence on the lithium abundance variations in halo stars. We show that in slowly rotating stars, under some conditions, lithium can be depleted by a factor of up to two with a dispersion smaller than 0.1 dex in the middle part of the lithium plateau. We derive a primordial lithium abundance of 2.5 +/- 0.1, consistent with the recent determinations of D/H and 4He/H.Comment: 15 pages, 10 figures. to be published in A&

Hydrodynamical instabilities induced by atomic diffusion in A stars and their consequences

Aims. Atomic diffusion, including the effect of radiative accelerations on individual elements, leads to important variations of the chemical composition inside the stars. The accumulation in specific layers of the elements, which are the main contributors of the local opacity, leads to hydrodynamical instabilities that modify the internal stellar structure and surface abundances. Our aim is to study these effects and compare the resulting surface abundances with spectroscopic observations Methods. We computed the detailed structure of A-type stars including these effects. We used the Toulouse-Geneva Evolution Code (TGEC), where radiative accelerations are computed using the Single Valued Parameter (SVP) method, and we added double-diffusive convection with mixing coefficients deduced from three-dimensional (3D) simulations. Results. We show that the modification of the initial chemical composition has important effects on the internal stellar mixing and leads to different surface abundances of the elements. The results fit the observed surface chemical composition well if the layers, which are individually mixed by double-diffusive convection, are connected.Comment: 7 pages, 7 figures, accepted in A&

The new Toulouse-Geneva Stellar Evolution Code including radiative accelerations of heavy elements

Atomic diffusion has been recognized as an important process that has to be considered in any computations of stellar models. In solar-type and cooler stars, this process is dominated by gravitational settling, which is now included in most stellar evolution codes. In hotter stars, radiative accelerations compete with gravity and become the dominant ingredient in the diffusion flux for most heavy elements. Introducing radiative accelerations into the computations of stellar models modifies the internal element distribution and may have major consequences on the stellar structure. Coupling these processes with hydrodynamical stellar motions has important consequences that need to be investigated in detail. We aim to include the computations of radiative accelerations in a stellar evolution code (here the TGEC code) using a simplified method (SVP) so that it may be coupled with sophisticated macroscopic motions. We also compare the results with those of the Montreal code in specific cases for validation and study the consequences of these coupled processes on accurate models of A- and early-type stars. We implemented radiative accelerations computations into the Toulouse-Geneva stellar evolution code following the semi-analytical prescription proposed by Alecian and LeBlanc. This allows more rapid computations than the full description used in the Montreal code. We present results for A-type stellar models computed with this updated version of TGEC and compare them with similar published models obtained with the Montreal evolution code. We discuss the consequences for the coupling with macroscopic motions, including thermohaline convection.Comment: 12 pages, 13 figures, published in A&