Dark Stars: Improved Models and First Pulsation Results

We use the stellar evolution code MESA to study dark stars. Dark stars (DSs), which are powered by dark matter (DM) self-annihilation rather than by nuclear fusion, may be the first stars to form in the Universe. We compute stellar models for accreting DSs with masses up to 10^6 M_{sun}. The heating due to DM annihilation is self-consistently included, assuming extended adiabatic contraction of DM within the minihalos in which DSs form. We find remarkably good overall agreement with previous models, which assumed polytropic interiors. There are some differences in the details, with positive implications for observability. We found that, in the mass range of 10^4 -10^5 M_{sun}, our DSs are hotter by a factor of 1.5 than those in Freese et al.(2010), are smaller in radius by a factor of 0.6, denser by a factor of 3 - 4, and more luminous by a factor of 2. Our models also confirm previous results, according to which supermassive DSs are very well approximated by (n=3)-polytropes. We also perform a first study of dark star pulsations. Our DS models have pulsation modes with timescales ranging from less than a day to more than two years in their rest frames, at z ~ 15, depending on DM particle mass and overtone number. Such pulsations may someday be used to identify bright, cool objects uniquely as DSs; if properly calibrated, they might, in principle, also supply novel standard candles for cosmological studies.Comment: 17 pages; 11 figures; revised version; accepted by Ap

The White Dwarf Luminosity Function from Sloan Digital Sky Survey Imaging Data

A sample of white dwarfs is selected from the Sloan Digital Sky Survey (SDSS) Data Release 3 using their reduced proper motions, based on improved proper motions from combined SDSS and USNO-B data. Numerous SDSS and follow-up spectra (Kilic and coworkers) are used to quantify completeness and contamination of the sample; kinematicsmodels are used to understand and correct for velocity-dependent selection biases.A luminosity function is constructed covering the range 7 \u3c Mbol \u3c 16, and its sensitivity to various assumptions and selection limits is discussed. The white dwarf luminosity function based on 6000 stars is remarkably smooth and rises nearly monotonically to Mbol=15.3. It then drops abruptly, although the small number of low-luminosity stars in the sample and their unknown atmospheric composition prevent quantitative conclusions about this decline. Stars are identified that may have high tangential velocities, and a preliminary luminosity function is constructed for them

Whole Earth Telescope Observations of the Helium Interacting Binary PG 1346+082 (CR Bootis)

We present our analysis of 240 hr of white-light, high-speed photometry of the dwarf nova-like helium variable PG 1346+082 (CR Boo). We identify two frequencies in the low-state power spectrum, at 679.670 ± 0.004 μHz and 669.887 ± 0.008 μHz. The 679.670 μHz variation is coherent over at least a 2 week time span, the first demonstration of a phase-coherent photometric variation in any dwarf nova-like interacting binary white dwarf system. The high-state power spectrum contains a complex fundamental with a frequency similar, but not identical, to the low-state spectrum, and a series of harmonics not detected in low state. We also uncover an unexpected dependence of the high-frequency power\u27s amplitude and frequency structure on overall system magnitude. We discuss these findings in light of the general AM CVn system model, particularly the implications of the high-order harmonics on future studies of disk structure, mass transfer, and disk viscosity

The White Dwarf Luminosity Function from SDSS Imaging Data

A sample of white dwarfs is selected from SDSS DR3 imaging data using their reduced proper motions, based on improved proper motions from SDSS plus USNO-B combined data. Numerous SDSS and followup spectra (Kilic et al. 2005) are used to quantify completeness and contamination of the sample; kinematic models are used to understand and correct for velocity-dependent selection biases. A luminosity function is constructed covering the range 7 < M_bol < 16, and its sensitivity to various assumptions and selection limits is discussed. The white dwarf luminosity function based on 6000 stars is remarkably smooth, and rises nearly monotonically to M_bol = 15.3. It then drops abruptly, although the small number of low-luminosity stars in the sample and their unknown atmospheric composition prevent quantitative conclusions about this decline. Stars are identified that may have high tangential velocities, and a preliminary luminosity function is constructed for them.Comment: Accepted for AJ (Jan 2006). 35 pages (includes 10 figures

A New Very Cool White Dwarf Discovered by the Sloan Digital Sky Survey

Early data taken during commissioning of the SDSS have resulted in the discovery of a very cool white dwarf. It appears to have stronger collision induced absorption from molecular hydrogen than any other known white dwarf, suggesting it has a cooler temperature than any other. While its distance is presently unknown, it has a surprisingly small proper motion, making it unlikely to be a halo star. An analysis of white dwarf cooling times suggests that this object may be a low-mass star with a helium core. The SDSS imaging and spectroscopy also recovered LHS 3250, the coolest previously known white dwarf, indicating that the SDSS will be an effective tool for identifying these extreme objects.Comment: 15 pages, including 5 figures. Accepted for Astrophysical Journal Letter

Constraining the Evolution of Zz Ceti

We report our analysis of the stability of pulsation periods in the DAV star (pulsating hydrogen atmosphere white dwarf) ZZ Ceti, also called R548. On the basis of observations that span 31 years, we conclude that the period 213.13 s observed in ZZ Ceti drifts at a rate dP/dt ≤ (5:5 ± 1:9) x 10-15 s s-1, after correcting for proper motion. Our results are consistent with previous Ṗ values for this mode and an improvement over them because of the larger time base. The characteristic stability timescale implied for the pulsation period is ⎸P / Ṗ ⎸=⎹≥ 1:2 Gyr, comparable to the theoretical cooling timescale for the star. Our current stability limit for the period 213.13 s is only slightly less than the present measurement for another DAV, G117-B15A, for the period 215.2 s, establishing this mode in ZZ Ceti as the second most stable optical clock known, comparable to atomic clocks and more stable than most pulsars. Constraining the cooling rate of ZZ Ceti aids theoretical evolutionary models and white dwarf cosmochronology. The drift rate of this clock is small enough that we can set interesting limits on reflex motion due to planetary companions

Limits on planets around pulsating white dwarf stars

We present limits on planetary companions to pulsating white dwarf stars. A subset of these stars exhibit extreme stability in the period and phase of some of their pulsation modes; a planet can be detected around such a star by searching for periodic variations in the arrival time of these pulsations. We present limits on companions greater than a few Jupiter masses around a sample of 15 white dwarf stars as part of an ongoing survey. One star shows a variation in arrival time consistent with a 2MJ planet in a 4.5 yr orbit.We discuss other possible explanations for the observed signal and conclude that a planet is the most plausible explanation based on the data available

Evidence for temperature change and oblique pulsation from light curve fits of the pulsating white dwarf GD 358

Convective driving, the mechanism originally proposed by Brickhill for pulsating white dwarf stars, has gained general acceptance as the generic linear instability mechanism in DAV and DBV white dwarfs. This physical mechanism naturally leads to a nonlinear formulation, reproducing the observed light curves of many pulsating white dwarfs. This numerical model can also provide information on the average depth of a star’s convection zone and the inclination angle of its pulsation axis. In this paper, we give two sets of results of nonlinear light curve fits to data on the DBV GD 358. Our first fit is based on data gathered in 2006 by the Whole Earth Telescope; this data set was multiperiodic containing at least 12 individual modes. Our second fit utilizes data obtained in 1996, when GD 358 underwent a dramatic change in excited frequencies accompanied by a rapid increase in fractional amplitude; during this event it was essentially monoperiodic. We argue that GD 358’s convection zone was much thinner in 1996 than in 2006, and we interpret this as a result of a short-lived increase in its surface temperature. In addition, we find strong evidence of oblique pulsation using two sets of evenly split triplets in the 2006 data. This marks the first time that oblique pulsation has been identified in a variable white dwarf star

Eleven new DA white dwarf variable stars from the Sloan Digital Sky Survey

We report the discovery of 11 new variable DAwhite dwarf (ZZ Ceti) stars. Candidates were selected by deriving temperatures from model fits to spectra obtained from the Sloan Digital Sky Survey (SDSS). We also find objects whose temperatures and gravities indicate they lie within the instability strip for pulsation but were not observed to vary. Although the temperatures are based on relatively low signal-to-noise ratio spectra, an impure strip is unexpected, which if confirmed, has implications for DA asteroseismology. This work brings the total number of published variable DA white dwarf stars to 82