M. H. Montgomery

Modules for Experiments in Stellar Astrophysics (MESA): Planets, Oscillations, Rotation, and Massive Stars

We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA star. Improvements in MESA stars ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 M ? stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA star solves the fully coupled stellar structure and composition equations, and we show how this has improved the scaling of MESAs calculational speed on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit that packages all the required components needed to form a unified, maintained, and well-validated build environment for MESA. We also highlight a few tools developed by the community for rapid visualization of MESA star results.

Whole earth telescope observations of the DBV white dwarf GD 358

We report on the analysis of 154 hours of early continuous high-speed photometry on the pulsating DB white dwarf (DBV) GD 358, obtained during the Whole Earth Telescope (WET) run of 1990 May. The power spectrum of the light curve is dominated by power in the range from 1000 to 2400 microHz with more than 180 significant peaks in the total spectrum. We identify all of the triplet frequencies as degree l = 1, and from the details of their spacings we derive the total stellar mass as 0.61 + or - 0.03 solar mass, the mass of the outer helium envelope as 2.0 + or - 1.0 x 10(exp -6) M(sub *), the absolute luminosity as 0.050 + or - 0.012 solar luminosity and the distance as 42 + or - 3 pc. We find strong evidence for differential rotation in the radial direction -- the outer envelope is rotating at least 1.8 times faster than the core -- and we detect the presence of a weak magnetic field with a strength of 1300 + or - 300 G. We also find a significant power at the sums and differences of the dominant frequencies, indicating nonlinear processes are significant, but they have a richness and complexity that rules out resonant mode coupling as a major cause.

Ensemble Characteristics of the ZZ Ceti Stars

We present the observed pulsation spectra of all known noninteracting ZZ Ceti stars (hydrogen atmosphere white dwarf variables [DAVs]) and examine changes in their pulsation properties across the instability strip. We confirm the well-established trend of increasing pulsation period with decreasing effective temperature across the ZZ Ceti instability strip. We do not find a dramatic order-of-magnitude increase in the number of observed independent modes in ZZ Ceti stars, traversing from the hot to the cool edge of the instability strip; we find that the cool DAVs have one more mode on average than the hot DAVs. We confirm the initial increase in pulsation amplitude at the blue edge and find strong evidence of a decline in amplitude prior to the red edge. We present the first observational evidence that ZZ Ceti stars lose pulsation energy just before pulsations shut down at the empirical red edge of the instability strip.

A New Technique for Probing Convection in Pulsating White Dwarf Stars

In this paper we demonstrate how pulsating white dwarfs can be used as an astrophysical laboratory for empirically constraining convection in these stars. We do this using a technique for fitting observed nonsinusoidal light curves, which allows us to recover the thermal response timescale of the convection zone (its depth), as well as demonstrating how this timescale changes as a function of effective temperature. We also obtain constraints on mode identifications for the pulsation modes, allowing us to use asteroseismology to study the interior structure of these stars. Aspects of this approach may have relevance for other classes of pulsators, including the Cepheids and RR Lyrae stars.

Binaries discovered by the SPY project. IV. Five single-lined DA double white dwarfs

We present results from our ongoing follow-up observations of double white dwarf binaries detected in the ESO SN Ia Progenitor SurveY (SPY). We discuss our observing strategy and data analysis and present the orbital solutions of five close double white dwarf binaries: HE0320−1917, HE1511−0448, WD0326−273, WD1013−010 and WD1210+140. Their periods range from 0.44 to 3.22 days. In none of these systems we find any spectral lines originating from the companion. This rules out main sequence companions and indicates that the companion white dwarfs are significantly older and cooler than the bright component. Infrared photometry suggests the presence of a cool, helium-rich white dwarf companion in the binary WD 0326−273. We briefly discuss the consequences of our findings for our understanding of the formation and evolution of double white dwarfs.

The effect of crystallization on the pulsations of white dwarf stars

We consider the pulsational properties of white dwarf star models with temperatures appropriate for the ZZ Ceti instability strip and with masses large enough that they should be substantially crystallized. Our work is motivated by the existence of a potentially crystallized DA variable (DAV), BPM 37093, and the expectation that digital surveys in progress will yield many more such massive pulsators. A crystallized core makes possible a new class of oscillations, the torsional modes, although we expect these modes to couple at most weakly to any motions in the fluid and therefore to remain unobservable. The p-modes should be affected at the level of a few percent in period, but are unlikely to be present with observable amplitudes in crystallizing white dwarfs any more than they are in the other ZZ Cetis. Most relevant to the observed light variations in white dwarfs are the g-modes. We find that the kinetic energy of these modes is effectively excluded from the crystallized cores of our models. As increasing crystallization pushes these modes farther out from the center, the mean period spacing ΔP between radial overtones increases substantially with the crystallized mass fraction, Mcr/M*. In addition, the degree and structure of mode trapping is affected. The fact that some periods are strongly affected by changes in the crystallized mass fraction while others are not suggests that we may be able to disentangle the effects of crystallization from those due to different surface layer masses.

New evolutionary models for massive ZZ Ceti stars. I. First results for their pulsational properties

We present new and improved evolutionary calculations for carbon-oxygen white dwarf (WD) stars appropriate for the study of massive ZZ Ceti stars. We follow the complete evolution of massive WD progenitors from the zero-age main sequence through the thermally pulsing and mass loss phases to the WD regime. Abundance changes are accounted for by means of a full coupling between nuclear evolution and time-dependent mixing due to diffusive overshoot, semiconvection and salt fingers. In addition, time-dependent element diffusion for multicomponent gases has been considered during the WD stage. We find that before the ZZ Ceti stage is reached, element diffusion has strongly smoothed out the chemical profile to such a degree that the resulting internal abundance distribution does not depend on the occurrence of overshoot episodes during the thermally pulsing phase. The mass of the hydrogen envelope left at the ZZ Ceti domain amounts to

Discovery of Pulsations, Including Possible Pressure Modes, in Two New Extremely Low Mass, He-core White Dwarfs

M_H \approx 2.3 \times 10^{-6}

A Strong Test of Electroweak Theory Using Pulsating DB White Dwarf Stars as Plasmon Neutrino Detectors

\msun. This is about half as large as for the case when element diffusion is neglected. The implications of our new models for the pulsational properties of massive ZZ Ceti stars are discussed. In this regard, we find that the occurrence of core overshooting during central helium burning leaves strong imprints on the theoretical period spectrum of massive ZZ Ceti stars. Finally, we present a simple new prescription for calculating the He/H profile which goes beyond the trace element approximation.

EVOLUTIONARY CALCULATIONS OF PHASE SEPARATION IN CRYSTALLIZING WHITE DWARF STARS

We report the discovery of the second and third pulsating extremely low mass white dwarfs (WDs), SDSS J111215.82+111745.0 (hereafter J1112) and SDSS J151826.68+065813.2 (hereafter J1518). Both have masses < 0.25 Msun and effective temperatures below 10,000 K, establishing these putatively He-core WDs as a cooler class of pulsating hydrogen-atmosphere WDs (DAVs, or ZZ Ceti stars). The short-period pulsations evidenced in the light curve of J1112 may also represent the first observation of acoustic (p-mode) pulsations in any WD, which provide an exciting opportunity to probe this WD in a complimentary way compared to the long-period g-modes also present. J1112 is a Teff = 9590 +/- 140 K and log(g) = 6.36 +/- 0.06 WD. The star displays sinusoidal variability at five distinct periodicities between 1792-2855 s. In this star we also see short-period variability, strongest at 134.3 s, well short of expected g-modes for such a low-mass WD. The other new pulsating WD, J1518, is a Teff = 9900 +/- 140 K and log(g) = 6.80 +/- 0.05 WD. The light curve of J1518 is highly non-sinusoidal, with at least seven significant periods between 1335-3848 s. Consistent with the expectation that ELM WDs must be formed in binaries, these two new pulsating He-core WDs, in addition to the prototype SDSS J184037.78+642312.3, have close companions. However, the observed variability is inconsistent with tidally induced pulsations and is so far best explained by the same hydrogen partial-ionization driving mechanism at work in classic C/O-core ZZ Ceti stars.