• Using images from the Hubble Space Telescope Advanced Camera for Surveys, we measure the rate of cooling of white dwarfs in the globular cluster 47 Tucanae and compare it to modelled cooling curves. We examine the effects of the outer convective envelope reaching the nearly isothermal degenerate core and the release of latent heat during core crystallisation on the white dwarf cooling rates. For white dwarfs typical of 47 Tuc, the onset of these effects occur at similar times. The latent heat released during crystallisation is a small heat source. In contrast, the heat reservoir of the degenerate core is substantially larger. When the convective envelope reaches the nearly isothermal interior of the white dwarf, the star becomes brighter than it would be in the absence of this effect. Our modelled cooling curves that include this convective coupling closely match the observed luminosity function of the white dwarfs in 47 Tuc.
  • The broadening of the hydrogen lines during flares is thought to result from increased charge (electron, proton) density in the flare chromosphere. However, disagreements between theory and modeling prescriptions have precluded an accurate diagnostic of the degree of ionization and compression resulting from flare heating in the chromosphere. To resolve this issue, we have incorporated the unified theory of electric pressure broadening of the hydrogen lines into the non-LTE radiative transfer code RH. This broadening prescription produces a much more realistic spectrum of the quiescent, A0 star Vega compared to the analytic approximations used as a damping parameter in the Voigt profiles. We test recent radiative-hydrodynamic (RHD) simulations of the atmospheric response to high nonthermal electron beam fluxes with the new broadening prescription and find that the Balmer lines are over-broadened at the densest times in the simulations. Adding many simultaneously heated and cooling model loops as a "multithread" model improves the agreement with the observations. We revisit the three-component phenomenological flare model of the YZ CMi Megaflare using recent and new RHD models. The evolution of the broadening, line flux ratios, and continuum flux ratios are well-reproduced by a multithread model with high-flux nonthermal electron beam heating, an extended decay phase model, and a "hot spot" atmosphere heated by an ultrarelativistic electron beam with reasonable filling factors: 0.1%, 1%, and 0.1% of the visible stellar hemisphere, respectively. The new modeling motivates future work to understand the origin of the extended gradual phase emission.
  • We present a robust statistical analysis of the white dwarf cooling sequence in 47 Tucanae. We combine HST UV and optical data in the core of the cluster, Modules for Experiments in Stellar Evolution (MESA) white dwarf cooling models, white dwarf atmosphere models, artificial star tests, and a Markov Chain Monte Carlo (MCMC) sampling method to fit white dwarf cooling models to our data directly. We use a technique known as the unbinned maximum likelihood to fit these models to our data without binning. We use these data to constrain neutrino production and the thickness of the hydrogen layer in these white dwarfs. The data prefer thicker hydrogen layers $(q_\mathrm{H}=3.2\e{-5})$ and we can strongly rule out thin layers $(q_\mathrm{H}=10^{-6})$. The neutrino rates currently in the models are consistent with the data. This analysis does not provide a constraint on the number of neutrino species.
  • Gaia will identify several 1e5 white dwarfs, most of which will be in the solar neighborhood at distances of a few hundred parsecs. Ground-based optical follow-up spectroscopy of this sample of stellar remnants is essential to unlock the enormous scientific potential it holds for our understanding of stellar evolution, and the Galactic formation history of both stars and planets.
  • Characterizing the local space density of double degenerate binary systems is a complementary approach to broad sky surveys of double degenerates to determine the expected rates of white dwarf binary mergers, in particular those that may evolve into other observable phenomena such as extreme helium stars, Am CVn systems, and supernovae Ia. However, there have been few such systems detected in local space. We report here the discovery that WD 1242$-$105, a nearby bright WD, is a double-line spectroscopic binary consisting of two degenerate DA white dwarfs of similar mass and temperature, despite it previously having been spectroscopically characterized as a single degenerate. Follow-up photometry, spectroscopy, and trigonometric parallax have been obtained in an effort to determine the fundamental parameters of each component of this system. The binary has a mass ratio of 0.7 and a trigonometric parallax of 25.5 mas, placing it at a distance of 39 pc. The system's total mass is 0.95 M$_\odot$ and has an orbital period of 2.85 hours, making it the strongest known gravitational wave source ($\log h = -20.78$) in the mHz regime. Because of its orbital period and total mass, WD 1242$-$105 is predicted to merge via gravitational radiation on a timescale of 740 Myr, which will most likely not result in a catastrophic explosion.
  • Using images from the Hubble Space Telescope Wide-Field Camera 3, we measure the rate of diffusion of stars through the core of the globular cluster 47 Tucanae using a sample of young white dwarfs identified in these observations. This is the first direct measurement of diffusion due to gravitational relaxation. We find that the diffusion rate $\kappa\approx 10-13$ arcsecond$^2$ Myr$^{-1}$ is consistent with theoretical estimates of the relaxation time in the core of 47 Tucanae of about 70 Myr.
  • We establish new constraints on the intermediate-mass range of the initial-final mass relation by studying white dwarfs in four young star clusters, and apply the results to study the evolution of stars on the thermally pulsing asymptotic giant branch (TP-AGB). We show that the stellar core mass on the AGB grows rapidly from 10% to 30% for stars with $M_{\rm initial}$ = 1.6 to 2.0 $M_\odot$. At larger masses, the core-mass growth decreases steadily to $\sim$10% at $M_{\rm initial}$ = 3.4 $M_\odot$. These observations are in excellent agreement with predictions from the latest TP-AGB evolutionary models in Marigo et al. (2013). We also compare to models with varying efficiencies of the third dredge-up and mass loss, and demonstrate that the process governing the growth of the core is largely the stellar wind, while the third dredge-up plays a secondary, but non-negligible role. Based on the new white dwarf measurements, we perform an exploratory calibration of the most popular mass-loss prescriptions in the literature. Finally, we estimate the lifetime and the integrated luminosity of stars on the TP-AGB to peak at $t$ $\sim$ 3 Myr and $E$ = 1.2 $\times$ 10$^{10}$ $L_\odot$ yr for $M_{\rm initial}$ $\sim$ 2 $M_\odot$ ($t$ $\sim$ 2 Myr for luminosities brighter than the RGB tip at $\log(L/L_{\odot})$ $>$ 3.4), decreasing to $t$ = 0.4 Myr and $E$ = 6.1 $\times$ 10$^{9}$ $L_\odot$ yr for stars with $M_{\rm initial}$ $\sim$ 3.5 $M_\odot$. The implications of these results are discussed with respect to general population synthesis studies that require correct modeling of the TP-AGB phase of stellar evolution.
  • We present a grid of optical (3200--7200 \AA) synthetic spectra calculated with Tlusty/Synspec. The new NLTE model atmospheres include the most recent hydrogen Stark broadening profiles; were calculated in opacity sampling and limited to pure H/He composition. The grid covers the observed parameter space of (He-)sdB and (He-)sdO stars, therefore it is suitable for the homogeneous spectral analyses of such evolved stars.