We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.

OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.

Variations in the pulsation arrival time of five independent pulsation frequencies of the DB white dwarf EC 20058-5234 individually imitate the effects of reflex motion induced by a planet or companion but are inconsistent when considered in unison. The pulsation frequencies vary periodically in a 12.9 year cycle and undergo secular changes that are inconsistent with simple neutrino plus photon-cooling models. The magnitude of the periodic and secular variations increases with the period of the pulsations, possibly hinting that the corresponding physical mechanism is located near the surface of the star. The phase of the periodic variations appears coupled to the sign of the secular variations. The standards for pulsation-timing-based detection of planetary companions around pulsating white dwarfs, and possibly other variables such as subdwarf B stars, should be re-evaluated. The physical mechanism responsible for this surprising result may involve a redistribution of angular momentum or a magnetic cycle. Additionally, variations in a supposed combination frequency are shown to match the sum of the variations of the parent frequencies to remarkable precision, an expected but unprecedented confirmation of theoretical predictions.

We report on analysis of 308.3 hrs of high speed photometry targeting the pulsating DA white dwarf EC14012-1446. The data were acquired with the Whole Earth Telescope (WET) during the 2008 international observing run XCOV26. The Fourier transform of the light curve contains 19 independent frequencies and numerous combination frequencies. The dominant peaks are 1633.907, 1887.404, and 2504.897 microHz. Our analysis of the combination amplitudes reveals that the parent frequencies are consistent with modes of spherical degree l=1. The combination amplitudes also provide m identifications for the largest amplitude parent frequencies. Our seismology analysis, which includes 2004--2007 archival data, confirms these identifications, provides constraints on additional frequencies, and finds an average period spacing of 41 s. Building on this foundation, we present nonlinear fits to high signal-to-noise light curves from the SOAR 4.1m, McDonald 2.1m, and KPNO 2m telescopes. The fits indicate a time-averaged convective response timescale of 99.4 +/- 17 s, a temperature exponent 85 +/- 6.2 and an inclination angle of 32.9 +/- 3.2 degrees. We present our current empirical map of the convective response timescale across the DA instability strip.

We present 13 years of pulsation timing measurements of the DBV white dwarf EC 2005-5234. Each of the four O-C diagrams mimic the sinusoidal behavior typically attributed to a planet + WD system. However, the amplitude and phase of the O-C variations are inconsistent with each other. We discuss the impact of this result on timing based WD planet searches.

Convective driving, the mechanism originally proposed by Brickhill (1991, 1983) 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 (WET); 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.

We constrain the distribution of calcium across the surface of the white dwarf star G29-38 by combining time series spectroscopy from Gemini-North with global time series photometry from the Whole Earth Telescope. G29-38 is actively accreting metals from a known debris disk. Since the metals sink significantly faster than they mix across the surface, any inhomogeneity in the accretion process will appear as an inhomogeneity of the metals on the surface of the star. We measure the flux amplitudes and the calcium equivalent width amplitudes for two large pulsations excited on G29-38 in 2008. The ratio of these amplitudes best fits a model for polar accretion of calcium and rules out equatorial accretion.

We report on the analysis of 436.1 hrs of nearly continuous high-speed photometry on the pulsating DB white dwarf GD358 acquired with the Whole Earth Telescope (WET) during the 2006 international observing run, designated XCOV25. The Fourier transform (FT) of the light curve contains power between 1000 to 4000 microHz, with the dominant peak at 1234 microHz. We find 27 independent frequencies distributed in 10 modes, as well as numerous combination frequencies. Our discussion focuses on a new asteroseismological analysis of GD358, incorporating the 2006 data set and drawing on 24 years of archival observations. Our results reveal that, while the general frequency locations of the identified modes are consistent throughout the years, the multiplet structure is complex and cannot be interpreted simply as l=1 modes in the limit of slow rotation. The high k multiplets exhibit significant variability in structure, amplitude and frequency. Any identification of the m components for the high k multiplets is highly suspect. The k=9 and 8 modes typically do show triplet structure more consistent with theoretical expectations. The frequencies and amplitudes exhibit some variability, but much less than the high k modes. Analysis of the k=9 and 8 multiplet splittings from 1990 to 2008 reveal a long-term change in multiplet splittings coinciding with the 1996 "sforzando" event, where GD358 dramatically altered its pulsation characteristics on a timescale of hours. We explore potential implications, including the possible connections between convection and/or magnetic fields and pulsations. We suggest future investigations, including theoretical investigations of the relationship between magnetic fields, pulsation, growth rates, and convection.