Seismological studies of ZZ Ceti stars – II. Application to the ZZ Ceti class

We used the detected pulsation modes and adiabatic pulsation models to do seismology of the class of ZZ Ceti stars and measure the H layer mass for 83 stars. We found the surface hydrogen layer to be within the range  10^−9.5≥ M _H/ M _*≥ 10⁻⁴  , with an average of   M _H/ M _*= 10^−6.3  , which is thinner than the predicted value of   M _H/ M _*= 10⁻⁴  , indicating that the stars lose more mass during their evolution than previously expected. These results are preliminary and do not include the possible effects of realistic C/O profiles on the fits.     

Evolution of a 3-M⊙ star from the main sequence to the ZZ Ceti stage: the role played by element diffusion

The purpose of this paper is to present new full evolutionary calculations for DA white dwarf stars with the major aim of providing a physically sound reference frame for exploring the pulsation properties of the resulting models in future communications. Here, white dwarf evolution is followed in a self-consistent way with the predictions of time-dependent element diffusion and nuclear burning. In addition, full account is taken of the evolutionary stages prior to white dwarf formation. In particular, we follow the evolution of a 3-M_⊙ model from the zero-age main sequence (the adopted metallicity is   Z =0.02)  , all the way from the stages of hydrogen and helium burning in the core up to the thermally pulsing phase. After experiencing 11 thermal pulses, the model is forced to evolve towards its white dwarf configuration by invoking strong mass loss episodes. Further evolution is followed down to the domain of the ZZ Ceti stars on the white dwarf cooling branch.
Emphasis is placed on the evolution of the chemical abundance distribution caused by diffusion processes and the role played by hydrogen burning during the white dwarf evolution. We find that discontinuities in the abundance distribution at the start of the cooling branch are considerably smoothed out by diffusion processes by the time the ZZ Ceti domain is reached. Nuclear burning during the white dwarf stage does not represent a major source of energy, as expected for a progenitor star of initially high metallicity. We also find that thermal diffusion lessens even further the importance of nuclear burning.
Furthermore, the implications of our evolutionary models for the main quantities relevant for adiabatic pulsation analysis are discussed. Interestingly, the shape of the Ledoux term is markedly smoother compared with previous detailed studies of white dwarfs. This is translated into a different behaviour of the Brunt–Väisälä frequency.     

The ages and colours of cool helium-core white dwarf stars

The purpose of this work is to explore the evolution of helium-core white dwarf stars in a self-consistent way with the predictions of detailed non-grey model atmospheres and element diffusion. To this end, we consider helium-core white dwarf models with stellar masses of 0.406, 0.360, 0.327, 0.292, 0.242, 0.196 and 0.169 M_⊙ and follow their evolution from the end of mass-loss episodes, during their pre-white dwarf evolution, down to very low surface luminosities.
We find that when the effective temperature decreases below 4000 K, the emergent spectrum of these stars becomes bluer within time-scales of astrophysical interest. In particular, we analyse the evolution of our models in the colour–colour and in the colour–magnitude diagrams and find that helium-core white dwarfs with masses ranging from ∼0.18 to 0.3 M_⊙ can reach the turn-off in their colours and become blue again within cooling times much less than 15 Gyr and then remain brighter than M _V ≈16.5 . In view of these results, many low-mass helium white dwarfs could have had enough time to evolve to the domain of collision-induced absorption from molecular hydrogen, showing blue colours.     

BVR observations and period variation of the near-contact binary ZZ Aurigae

BVR light curves of ZZ Aurigae were obtained with the 60-cm Cassegrain reflector at the Sobaek Observatory, Korea, between 2000 February and 2001 February. All collected times of minimum light, including our observations, were used for the period study. The period variation could be of quasi-sinusoidal form superposed on an upward parabola. A continuous period increase of  d P /d t =+2.3 × 10⁻⁸ d yr⁻¹  was determined for ZZ Aur. The period of quasi-sinusoidal variation is about ∼26–31 yr. Photometric solutions were found using the Wilson–Devinney method. The Roche configuration of ZZ Aur is that of an Algol-type semidetached system where the primary star nearly fills its Roche lobe and the secondary star fills its lobe. The spot model was used to explain the asymmetry in the light curve known as the O'Connell effect.     

Evolutionary calculations of carbon dredge-up in helium envelope white dwarfs

We investigate the evolution of cooling helium atmosphere white dwarfs using a full evolutionary code, specifically developed to follow the effects of element diffusion and gravitational settling on white dwarf cooling. The major difference between this work and previous work is that we use more recent opacity data from the OPAL project. Since, in general, these opacities are higher than those available 10 years ago, at a given effective temperature, convection zones go deeper than in models with older opacity data. Thus convective dredge-up of observationally detectable carbon in helium atmosphere white dwarfs can occur for thicker helium layers than found by Pelletier et al. We find that the range of observed C to He ratios in different DQ white dwarfs of similar effective temperature is well explained by a range of initial helium layer mass between 10⁻³ and 10⁻² M⊙, in good agreement with stellar evolution theory, assuming a typical white dwarf mass of 0.6 M⊙. We also predict that oxygen will be present in DQ white dwarf atmospheres in detectable amounts if the helium layer mass is near the lower limit compatible with stellar evolution theory. Determination of the oxygen abundance has the potential of providing information on the profile of oxygen in the core and hence on the important ¹²C(α,γ)¹⁶O reaction rate.     

The evolution of low-metallicity asymptotic giant branch stars and the formation of carbon-enhanced metal-poor stars

We investigate the behaviour of asymptotic giant branch (AGB) stars between metallicities   Z = 10⁻⁴  and 10⁻⁸. We determine which stars undergo an episode of flash-driven mixing, where protons are ingested into the intershell convection zone, as they enter the thermally pulsing AGB phase and which undergo third dredge-up. We find that flash-driven mixing does not occur above a metallicity of   Z = 10⁻⁵  for any mass of star and that stars above  2 M_⊙  do not experience this phenomenon at any metallicity. We find carbon ingestion (CI), the mixing of carbon into the tail of hydrogen-burning region, occurs in the mass range  2 M_⊙  to around  4 M_⊙  . We suggest that CI may be a weak version of the flash-driven mechanism. We also investigate the effects of convective overshooting on the behaviour of these objects. Our models struggle to explain the frequency of Carbon-Enhanced Metal-Poor (CEMP) stars that have both significant carbon and nitrogen enhancement. Carbon can be enhanced through flash-driven mixing, CI or just third dredge-up. Nitrogen can be enhanced through hot bottom burning and the occurrence of hot dredge-up also converts carbon into nitrogen. The C/N ratio may be a good indicator of the mass of the primary AGB stars.     

The eclipsing binary PV Cas: a case of strong differential rotation in depth

In this study, we show that only models with differential rotation, as opposed to non-rotation or solid-body rotation for the components of PV Cas, are in satisfactory agreement with the observations (including constraints on the apsidal advance rate and the synchronous rotation of the components in addition to their luminosities and radii). Internal rotation profiles found for solar and metal-rich chemical compositions are similar to each other in that only a small part of the outermost regions is rotating very slowly  ( v _eq∼ 65 km s ⁻¹)  and the rest is rotating very rapidly. Confirmation of Ap-like variations in the light curve of the system leads us to search for a connection between this type (and similar types) of internal rotation and chemically peculiar stars. The temperature difference between the blue sides of the main sequence for the normal and for the magnetic Ap stars may arise from such a connection, since we find a similar difference between the models with this kind of rotation and the non-rotating models of somewhat lower mass.     

Characterizing the pulsations of the ZZ Ceti star KUV 02464+3239

We present the results on period search and modelling of the cool DAV star KUV 02464+3239. Our observations resolved the multiperiodic pulsational behaviour of the star. In agreement with its position near the red edge of the DAV instability strip, it shows large amplitude, long-period pulsation modes, and has a strongly non-sinusoidal light curve. We determined six frequencies as normal modes and revealed remarkable short-term amplitude variations. A rigorous test was performed for the possible source of amplitude variation: beating of modes, effect of noise, unresolved frequencies or rotational triplets. Among the best-fitting models resulting from a grid search, we selected three that gave   l = 1  solutions for the largest amplitude modes. These models had masses of 0.645, 0.650 and  0.680  M _⊙  . The three 'favoured' models have   M _H  between  2.5 × 10⁻⁵ and 6.3 × 10⁻⁶  M _*  and give 14.2–14.8 mas seismological parallax. The  0.645  M _⊙  (11 400 K) model also matches the spectroscopic  log  g   and   T _eff  within 1σ. We investigated the possibility of mode trapping and concluded that while it can explain high amplitude modes, it is not required.     

Non-adiabatic oscillations of red giants

Using our non-local time-dependent theory of convection, the linear non-adiabatic oscillations of 10 evolutionary model series with masses of  1–3 M_⊙  are calculated. The results show that there is a red giant instability strip in the lower temperature side of the Hertzsprung–Russell diagram which goes along the sequences of the red giant branch and the asymptotic giant branch. For red giants of lower luminosities, pulsation instability is found at high order overtones; the lower order modes from the fundamental to the second overtone are stable. Towards higher luminosity and lower effective temperature, instability moves to lower order modes, and the amplitude growth rate of oscillations also grows. At the high luminosity end of the strip, the fundamental and the first overtone become unstable, while all the modes above the fourth order become stable. The excitation mechanisms have been studied in detail. It is found that turbulent pressure plays a key role for excitation of red variables. The frozen convection approximation is unavailable for the low temperature stars with extended convective envelopes. In any case, this approximation can explain neither the red edge of the Cepheid instability strip, nor the blue edge of the pulsating red giant instability strip. An analytic expression of a pulsation constant as a function of stellar mass, luminosity and effective temperature is presented from this work.     

Burial of the polar magnetic field of an accreting neutron star – II. Hydromagnetic stability of axisymmetric equilibria

The theory of polar magnetic burial in accreting neutron stars predicts that a mountain of accreted material accumulates at the magnetic poles of the star, and that, as the mountain spreads equatorward, it is confined by, and compresses, the equatorial magnetic field. Here, we extend previous, axisymmetric, Grad–Shafranov calculations of the hydromagnetic structure of a magnetic mountain up to accreted masses as high as   M _a= 6 × 10⁻⁴ M_⊙  , by importing the output from previous calculations (which were limited by numerical problems and the formation of closed bubbles to   M _a < 10⁻⁴ M_⊙  ) into the time-dependent, ideal-magnetohydrodynamic code zeus-3d and loading additional mass on to the star dynamically. The rise of buoyant magnetic bubbles through the accreted layer is observed in these experiments. We also investigate the stability of the resulting hydromagnetic equilibria by perturbing them in zeus-3d . Surprisingly, it is observed that the equilibria are marginally stable for all   M _a≤ 6 × 10⁻⁴ M_⊙  ; the mountain oscillates persistently when perturbed, in a combination of Alfvén and acoustic modes, without appreciable damping or growth, and is therefore not disrupted (apart from a transient Parker instability initially, which expels <1 per cent of the mass and magnetic flux).     

The discovery of photospheric nickel in the hot DO white dwarf REJ 0503−289

We present the first evidence for the direct detection of nickel in the photosphere of the hot DO white dwarf REJ 0503−289. While this element has been seen previously in the atmospheres of hot H-rich white dwarfs, this is one of the first similar discoveries in a He-rich object. Intriguingly, iron, which is observed to be more abundant than Ni in the hot DA stars, is not detected, the upper limit to its abundance (Fe/He=10⁻⁶) implying an Fe/Ni ratio a factor of 10 lower than seen in the H-rich objects (Ni/He=10⁻⁵ for REJ 0503−289). The abundances of nickel and various other elements heavier than He were determined from Goddard High Resolution Spectrograph spectra. We used two completely independent sets of non-local thermodynamic equilibrium model atmospheres, which both provide the same results. This not only reduces the possibility of systematic errors in our analysis, but is also an important consistency check for both model atmosphere codes.
We have also developed a more objective method of determining T _eff and log  g , from the He lines in the optical spectrum, in the form of a formal fitting of the line profiles to a grid of model spectra, an analogue of the standard procedure utilizing the Balmer lines in DA white dwarfs. This gives the assigned uncertainties in T _eff and log  g a firm statistical basis and allows us to demonstrate that inclusion of elements heavier than H, He and C in the spectral calculations, exclusively considered in most published optical analyses, yields a systematic downward shift in the measured value of T _eff.     

The effect of the 19F(α, p)22Ne reaction rate uncertainty on the yield of fluorine from Wolf–Rayet stars

In the light of recent recalculations of the  ¹⁹F(α, p)²²Ne  reaction rate, we present results of the expected yield of ¹⁹F from Wolf–Rayet (WR) stars. In addition to using the recommended rate, we have computed models using the upper and lower limits for the rate, and hence we constrain the uncertainty in the yield with respect to this reaction. We find a yield of  3.1 × 10⁻⁴ M_⊙  of ¹⁹F with our recommended rate, and a difference of a factor of 2 between the yields computed with the upper and lower limits. In comparison with previous work we find a difference in the yield of a factor of approximately 4, connected with a different choice of mass loss. Model uncertainties must be carefully evaluated in order to obtain a reliable estimate of the yield, together with its uncertainties, of fluorine from WR stars.     

Modelling the components of binaries in the Hyades: the dependence of the mixing-length parameter on stellar mass

We present our findings based on a detailed analysis of the binaries of the Hyades, in which the masses of the components are well known. We fit the models of the components of a binary system to observations so as to give the observed total V and B − V of that system and the observed slope of the main sequence in the corresponding parts. According to our findings, there is a very definite relationship between the mixing-length parameter and the stellar mass. The fitting formula for this relationship can be given as  α= 9.19( M /M_⊙− 0.74)^0.053− 6.65  , which is valid for stellar masses greater than  0.77 M_⊙  . While no strict information is gathered for the chemical composition of the cluster, as a result of degeneracy in the colour–magnitude diagram, by adopting   Z = 0.033  and using models for the components of 70 Tau and θ² Tau we find the hydrogen abundance to be   X = 0.676  and the age to be 670 Myr. If we assume that   Z = 0.024  , then   X = 0.718  and the age is 720 Myr. Our findings concerning the mixing-length parameter are valid for both sets of the solution. For both components of the active binary system V818 Tau, the differences between radii of the models with   Z = 0.024  and the observed radii are only about 4 per cent. More generally, the effective temperatures of the models of low-mass stars in the binary systems studied are in good agreement with those determined by spectroscopic methods.     

Spitzer observations of acetylene bands in carbon-rich asymptotic giant branch stars in the Large Magellanic Cloud

We investigate the molecular bands in carbon-rich asymptotic giant branch (AGB) stars in the Large Magellanic Cloud (LMC), using the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope ( SST ) over the 5–38 μm range. All 26 low-resolution spectra show acetylene (C₂H₂) bands at 7 and 14 μm. The hydrogen cyanide (HCN) bands at these wavelengths are very weak or absent. This is consistent with low nitrogen abundances in the LMC. The observed 14 μm C₂H₂  band is reasonably reproduced by an excitation temperature of 500 K. There is no clear dilution of the 14 μm C₂H₂  band by circumstellar dust emission. This 14-μm band originates from molecular gas in the circumstellar envelope in these high mass-loss rate stars, in agreement with previous findings for Galactic stars. The C₂H₂ column density, derived from the 13.7 μm band, shows a gas mass-loss rate in the range 3 × 10⁻⁶ to 5 × 10⁻⁵ M_⊙ yr⁻¹. This is comparable with the total mass-loss rate of these stars estimated from the spectral energy distribution. Additionally, we compare the line strengths of the 13.7 μm C₂H₂  band of our LMC sample with those of a Galactic sample. Despite the low metallicity of the LMC, there is no clear difference in the C₂H₂  abundance among LMC and Galactic stars. This reflects the effect of the third dredge-up bringing self-produced carbon to the surface, leading to high carbon-to-oxygen ratio at low metallicity.     

Hydrogen-model atmospheres for white dwarf stars

We present a detailed calculation of model atmospheres for DA white dwarfs. Our atmosphere code solves the atmosphere structure in local thermodynamic equilibrium with a standard partial linearization technique, which takes into account the energy transfer by radiation and convection. This code incorporates recent improved and extended data base of collision-induced absorption by molecular hydrogen. We analyse the thermodynamic structure and emergent flux of atmospheres in the range 2500 T _eff60 000 K and 6.5log  g 9.0. Bolometric correction and colour indices are provided for a subsample of the model grid. Comparison of the colours is made with published observational material and results of other recent model calculations.
Motivated by the increasing interest in helium-core white dwarfs, we analyse the photometric characteristics of these stars during their cooling, using evolutionary models recently available. Effective temperatures, surface gravities, masses and ages have been determined for some helium-core white dwarf candidates, and their possible binary nature is briefly discussed.     

A-star envelopes: a test of local and non-local models of convection

We present results of a fully non-local, compressible model of convection for A-star envelopes. This model quite naturally reproduces a variety of results from observations and numerical simulations which local models based on a mixing length do not. Our principal results, which are for models with T _eff between 7200 and 8500 K, are the following. First, the photospheric velocities and filling factors are in qualitative agreement with those derived from observations of line profiles of A-type stars. Secondly, the He  ii and H  i convection zones are separated in terms of convective flux and thermal interaction, but joined in terms of the convective velocity field, in agreement with numerical simulations. In addition, we attempt to quantify the amount of overshooting in our models at the base of the He  ii convection zone.     

An alternative symbiotic channel to Type Ia supernovae

By assuming an aspherical stellar wind with an equatorial disc from a red giant, we investigate the production of Type Ia supernovae (SNe Ia) via a symbiotic channel. We estimate that the Galactic birthrate of SNe Ia via the symbiotic channel is between  1.03 × 10⁻³  and  2.27 × 10⁻⁵ yr⁻¹  , while the delay time of SNe Ia has a wide range from ∼0.07 to 5 Gyr. The results are greatly affected by the outflow velocity and mass-loss rate of the equatorial disc. Using our model, we discuss the progenitors of SN 2002ic and SN 2006X.     

On the late spectral types of cataclysmic variable secondaries

We investigate why the spectral type of most cataclysmic variable (CV) secondaries is significantly later than that of a zero-age main-sequence (ZAMS) star with the same mean density. Using improved stellar input physics, tested against observations of low-mass stars at the bottom of the main sequence, we calculate the secular evolution of CVs with low-mass donors. We consider sequences with different mass transfer rates and with a different degree of nuclear evolution of the donor prior to mass transfer.
Systems near the upper edge of the gap ( P ∼3–6 h) can be reproduced by models with a wide range of mass transfer rates from 1.5×10⁻⁹ M_⊙ yr⁻¹ to 10⁻⁸ M_⊙ yr⁻¹. Evolutionary sequences with a small transfer rate and donors that are substantially evolved off the ZAMS (central hydrogen content 0.05–0.5) reproduce CVs with late spectral types above P ≳6 h. Systems with the most discrepant (late) spectral type should have the smallest donor mass at any given P .
Consistency with the period gap suggests that the mass transfer rate increases with decreasing donor mass for evolved sequences above the period gap. In this case, a single-parameter family of sequences with varying X _c and increasing mass transfer rate reproduces the full range of observed spectral types. This would imply that CVs with such evolved secondaries dominate the CV population.