We present a novel method to detect variable astrophysical objects and transient phenomena using anomalous excess scatter in repeated measurements from public catalogs of Gaia DR2 and Zwicky Transient Facility (ZTF) DR3 photometry. We first provide a generalized, all-sky proxy for variability using only Gaia DR2 photometry, calibrated to white dwarf stars. To ensure more robust candidate detection, we further employ a method combining Gaia with ZTF photometry and alerts. To demonstrate its efficacy, we apply this latter technique to a sample of roughly 12,100 white dwarfs within 200 pc centered on the ZZ Ceti instability strip, where hydrogen-atmosphere white dwarfs are known to pulsate. By inspecting the top 1% of the samples ranked by these methods, we demonstrate that both the Gaia-only and ZTF-informed techniques are highly effective at identifying known and new variable white dwarfs, which we verify using follow-up, high-speed photometry. We confirm variability in all 33 out of 33 (100%) observed white dwarfs within our top 1% highest-ranked candidates, both inside and outside the ZZ Ceti instability strip. In addition to dozens of new pulsating white dwarfs, we also identify five white dwarfs highly likely to show transiting planetary debris; if confirmed, these systems would more than triple the number of white dwarfs known to host transiting debris.

We report the discovery of pulsations in the extremely low-mass (ELM), likely helium-core white dwarf GD 278 via ground- and space-based photometry. GD 278 was observed by the Transiting Exoplanet Survey Satellite (TESS) in Sector 18 at a 2 minute cadence for roughly 24 days. The TESS data reveal at least 19 significant periodicities between 2447 and 6729 s, one of which is the longest pulsation period ever detected in a white dwarf. Previous spectroscopy found that this white dwarf is in a 4.61 hr orbit with an unseen >0.4 M ⊙ companion and has T eff = 9230 ± 100 K and log g = 6.627 ± 0.056 , which corresponds to a mass of 0.191 ± 0.013 M ⊙. Patterns in the TESS pulsation frequencies from rotational splittings appear to reveal a stellar rotation period of roughly 10 hr, making GD 278 the first ELM white dwarf with a measured rotation rate. The patterns inform our mode identification for asteroseismic fits, which, unfortunately, do not reveal a global best-fit solution. Asteroseismology reveals two main solutions roughly consistent with the spectroscopic parameters of this ELM white dwarf, but with vastly different hydrogen-layer masses; future seismic fits could be further improved by using the stellar parallax. GD 278 is now the tenth known pulsating ELM white dwarf; it is only the fifth known to be in a short-period binary, but is the first with extended, space-based photometry.

Hot subdwarf stars are mostly stripped red giants that can exhibit photometric variations due to stellar pulsations, eclipses, the reflection effect, ellipsoidal modulation, and Doppler beaming. Detailed studies of their light curves help constrain stellar parameters through asteroseismological analyses or binary light-curve modeling and generally improve our capacity to draw a statistically meaningful picture of this enigmatic stage of stellar evolution. From an analysis of Gaia DR2 flux errors, we have identified around 1200 candidate hot subdwarfs with inflated flux errors for their magnitudes—a strong indicator of photometric variability. As a pilot study, we obtained 2 minute cadence TESS Cycle 2 observations of 187 candidate hot subdwarfs with anomalous Gaia flux errors. More than 90% of our targets show significant photometric variations in their TESS light curves. Many of the new systems found are cataclysmic variables, but we report the discovery of several new variable hot subdwarfs, including HW Vir binaries, reflection-effect systems, pulsating sdBV s stars, and ellipsoidally modulated systems. We determine atmospheric parameters for select systems using follow-up spectroscopy from the 3 m Shane telescope. Finally, we present a Fourier diagnostic plot for classifying binary light curves using the relative amplitudes and phases of their fundamental and harmonic signals in their periodograms. This plot makes it possible to identify certain types of variables efficiently, without directly investigating their light curves, and may assist in the rapid classification of systems observed in large photometric surveys.