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

HD 51106 and HD 50747: an ellipsoidal binary and a triple system observed with CoRoT

We present an analysis of the observations of HD 51106 and HD 50747 by the satellite CoRoT, obtained during its initial run, and of the spectroscopic preparatory observations. AIMS: We complete an analysis of the light curve, extract the main frequencies observed, and discuss some preliminary interpretations about the stars. Methods: We used standard Fourier transform and pre-whitening methods to extract information about the periodicities of the stars. Results: HD 51106 is an ellipsoidal binary, the light curve of which can be completely explained by the tidal deformation of the star and smaller secondary effects. HD 50747 is a triple system containing a variable star, which exhibits many modes of oscillation with periods in the range of a few hours. On the basis of this period range and the analysis of the physical parameters of the star, we conclude that HD 50747 is a Gamma-Doradus star.Comment: 7 pages, 8 figures, use (Astronomy-Astrophysics format/macro LAtex

Asteroseismology of the planet-hosting star mu Arae. II. Seismic analysis

As most exoplanets host stars, HD 160691 (alias mu Ara) presents a metallicity excess in its spectrum compared to stars without detected planets. This excess may be primordial, in which case the star would be completely overmetallic, or it may be due to accretion in the early phases of planetary formation, in which case it would be overmetallic only in its outer layers. As discussed in a previous paper (Bazot and Vauclair 2004), seismology can help choosing between the two scenarios. This star was observed during eight nights with the spectrograph HARPS at La Silla Observatory. Forty three p-modes have been identified (Bouchy et al. 2005). In the present paper, we discuss the modelisation of this star. We computed stellar models iterated to present the same observable parameters (luminosity, effective temperature, outer chemical composition) while the internal structure was different according to the two extreme assumptions : original overmetallicity or accretion. We show that in any case the seismic constraints lead to models in complete agreement with the external parameters deduced from spectroscopy and from the Hipparcos parallax (L and Teff). We discuss the tests which may lead to a choice between the two typical scenarios. We show that the ``small separation'' seem to give a better fit for the accretion case than for the overmetallic case, but in spite of the very good data the uncertainties are still too large to conclude. We discuss the observations which would be needed to go further and solve this question.Comment: 16 pages, 8 figures, accepted in A&

The D/H Ratio in Interstellar Gas Towards G191-B2B

We reinvestigate the question of spatial variation of the local D/H abundance, using both archival GHRS spectra, and new echelle spectra of G191-B2B obtained with the Space Telescope Imaging Spectrograph (STIS) aboard HST. Our analysis uses stratified line-blanketed non-LTE model atmosphere calculations to determine the shape of the intrinsic WD Lyman-alpha profile and estimate the WD photospheric contamination of the interstellar lines. Although three velocity components were reported previously towards G191-B2B, we detect only two velocity components. The first component is at V(hel) ~ 8.6 km/s and the second at V(hel) ~ 19.3 km/s, which we identify with the Local Interstellar Cloud (LIC). From the STIS data we derive D/H = 1.60(+0.39,-0.27)X10^-5 for the LIC component, and D/H > 1.26X10^-5 for the 8.6 km/s component (uncertainties denote 2-sigma or 95% confidence limits). The STIS data provide no evidence for local or component-to-component variation in the D/H ratio. Despite using two velocity components for the profile fitting and using a more physically realistic WD Lyman-alpha profile for G191-B2B, our re-analysis of the GHRS data indicates a component-to-component variation as well as a variation of the D/H ratio in the LISM, neither of which are supported by the newer STIS data. We believe the most probable cause for this difference is the characterization of the background due to scattered light in the GHRS and STIS spectrographs. The two-dimensional MAMA detectors of STIS measure both the spatial and wavelength dependences of scattered light, allowing more accurate scattered light corrections than was possible with GHRS.Comment: Accepted for publication in Astrophysical Journal Letters. 10 pages + 3 figures. (Abstract is abridged.

Models of Metal Poor Stars with Gravitational Settling and Radiative Accelerations: I. Evolution and Abundance Anomalies

Evolutionary models have been calculated for Pop II stars of 0.5 to 1.0$M_\odot$ from the pre-main-sequence to the lower part of the giant branch. Rosseland opacities and radiative accelerations were calculated taking into account the concentration variations of 28 chemical species, including all species contributing to Rosseland opacities in the OPAL tables. The effects of radiative accelerations, thermal diffusion and gravitational settling are included. While models were calculated both for Z=0.00017 and 0.0017, we concentrate on models with Z=0.00017 in this paper. These are the first Pop II models calculated taking radiative acceleration into account. It is shown that, at least in a 0.8$M_\odot$ star, it is a better approximation not to let Fe diffuse than to calculate its gravitational settling without including the effects of $g_{rad}(Fe)$. In the absence of any turbulence outside of convection zones, the effects of atomic diffusion are large mainly for stars more massive than 0.7$M_\odot$. Overabundances are expected in some stars with \teff \ge 6000K. Most chemical species heavier than CNO are affected. At 12 Gyr, overabundance factors may reach 10 in some cases (e.g. for Al or Ni) while others are limited to 3 (e.g. for Fe). The calculated surface abundances are compared to recent observations of abundances in globular clusters as well as to observations of Li in halo stars. It is shown that, as in the case of Pop I stars, additional turbulence appears to be present.Comment: 40 pages, 17 color figures, to appear in The Astrophysical Journal, April 2002 (paper with original high resolution figures can be found at http://www.cerca.umontreal.ca/~richer/Fichiersps/popII_1.ps

Mass-Radius Relation for Magnetic White Dwarfs

Recently, several white dwarfs with very strong surface magnetic fields have been observed. In this paper we explore the possibility that such stars could have sufficiently strong internal fields to alter their structure. We obtain a revised white dwarf mass-radius relation in the presence of strong internal magnetic fields. We first derive the equation of state for a fully degenerate ideal electron gas in a magnetic field using an Euler-MacLaurin expansion. We use this to obtain the mass-radius relation for magnetic $^{4}$He, $^{12}$C, and $^{56}$Fe white dwarfs of uniform composition.Comment: 7 pages, 7 figures and 1 table, To appear in Ap

Low Mass Stars and the He3 Problem

The prediction of standard chemical evolution models of higher abundances of He3 at the solar and present-day epochs than are observed indicates a possible problem with the yield of He3 for stars in the range of 1-3 solar masses. Because He3 is one of the nuclei produced in Big Bang Nucleosynthesis (BBN), it is noted that galactic and stellar evolution uncertainties necessarily relax constraints based on He3. We incorporate into chemical evolution models which include outflow, the new yields for He3 of Boothroyd & Malaney (1995) which predict that low mass stars are net destroyers of He3. Since these yields do not account for the high \he3/H ratio observed in some planetary nebulae, we also consider the possibility that some fraction of stars in the 1 - 3 solar mass range do not destroy their He3 in theirpost main-sequence phase. We also consider the possibility that the gas expelled by stars in these mass ranges does not mix with the ISM instantaneously thus delaying the He3 produced in these stars, according to standard yields, from reaching the ISM. In general, we find that the Galactic D and He3 abundances can be fit regardless of whether the primordial D/H value is high (2 x 10^{-4}) or low (2.5 x 10^{-5}).Comment: 20 pages, latex, 9 ps figure

XMM−NewtonXMM-Newton Ω\Omega project: III. Gas mass fraction shape in high redshift clusters

We study the gas mass fraction, $f\_{\rm gas},$ behavior in $XMM-Newton$ $\Omega$ project. The typical $f\_{\rm gas}$ shape of high redshift galaxy clusters follows the global shape inferred at low redshift quite well. This result is consistent with the gravitational instability picture leading to self similar structures for both the dark and baryonic matter. However, the mean $f\_{\rm gas} in distant clusters shows some differences to local ones, indicating a departure from strict scaling. This result is consistent with the observed evolution in the luminosity-temperature relation. We quantitatively investigate this departure from scaling laws. Within the local sample we used, a moderate but clear variation of the amplitude of the gas mass fraction with temperature is found, a trend that weakens in the outer regions. These variations do not explain departure from scaling laws of our distant clusters. An important implication of our results is that the gas fraction evolution, a test of the cosmological parameters, can lead to biased values when applied at radii smaller than the virial radius. From our$XMM$ clusters, the apparent gas fraction at the virial radius is consistent with a non-evolving universal value in a high matter density model and not with a concordance.Comment: Accepted, A&A, in pres