Results from regular monitoring of relativistic compact binaries like PSR 1913+16 are consistent with the dominant (quadrupole) order emission of gravitational waves (GWs). We show that observations ...associated with the binary black hole (BBH) central engine of blazar OJ 287 demand the inclusion of gravitational radiation reaction effects beyond the quadrupolar order. It turns out that even the effects of certain hereditary contributions to GW emission are required to predict impact flare timings of OJ 287. We develop an approach that incorporates this effect into the BBH model for OJ 287. This allows us to demonstrate an excellent agreement between the observed impact flare timings and those predicted from ten orbital cycles of the BBH central engine model. The deduced rate of orbital period decay is nine orders of magnitude higher than the observed rate in PSR 1913+16, demonstrating again the relativistic nature of OJ 287's central engine. Finally, we argue that precise timing of the predicted 2019 impact flare should allow a test of the celebrated black hole "no-hair theorem" at the 10% level.
We present the results of our power spectral density analysis for the BL Lac object OJ 287, utilizing the Fermi-LAT survey at high-energy γ-rays, Swift-XRT in X-rays, several ground-based telescopes ...and the Kepler satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time autoregressive moving average (CARMA) processes. Owing to the inclusion of the Kepler data, we were able to construct for the first time the optical variability power spectrum of a blazar without any gaps across ∼6 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 γ-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 Fermi-LAT data, corresponding to ∼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.
ABSTRACT OJ 287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts that are predictable in a binary black hole model. The model predicted a major optical ...outburst in 2015 December. 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, . 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 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 an analysis of 308.3 hr of high-speed photometry targeting the pulsating DA white dwarf EC14012-1446. The data were acquired with the Whole Earth Telescope 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 mu Hz. 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.1 m, McDonald 2.1 m, and KPNO 2 m telescopes. The fits indicate a time-averaged convective response timescale of tau sub(0) = 99.4 + or - 17 s, a temperature exponent N = 85 + or - 6.2, and an inclination angle of theta sub(i) = 32degrees.9 + or - 3degrees.2. We present our current empirical map of the convective response timescale across the DA instability strip.
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.
The achievement of suitable toroidal-current-density profiles in tokamak plasmas plays an important role in enabling high fusion gain and noninductive sustainment of the plasma current for ...steady-state operation with improved magnetohydrodynamic stability. The evolution in time of the current profile is related to the evolution of the poloidal magnetic flux, which is modeled in normalized cylindrical coordinates using a partial differential equation (PDE) usually referred to as the magnetic flux diffusion equation. The dynamics of the plasma current density profile can be modified by the total plasma current and the power of the noninductive current drive. These two actuators, which are constrained not only in value and rate but also in their initial and final values, are used to drive the current profile as close as possible to a desired target profile at a specific final time. To solve this constrained finite-time open-loop PDE optimal control problem, model reduction based on proper orthogonal decomposition is combined with sequential quadratic programming in an iterative fashion. The use of a low-dimensional dynamical model dramatically reduces the computational effort and, therefore, the time required to solve the optimization problem, which is critical for a potential implementation of a real-time receding-horizon control strategy.
We analyse the light curve in the R band of the blazar OJ287, gathered during the 2015/2016 observing season. We did a search for quasi-periodic oscillations (QPOs) using several methods over a wide ...range of timescales. No statistically significant periods were found in the high-frequency domain both in the ground-based data and in Kepler observations. In the longer-period domain, the Lomb-Scargle periodogram revealed several peaks above the 99% significance level. The longest one-about 95 days-corresponds to the innermost stable circular orbit (ISCO) period of the more massive black hole. The 43-day period could be an alias, or it can be attributed to accretion in the form of a two-armed spiral wave.
We report on the analysis of 436.1 hr 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 and 4000 mHz, with the dominant peak at 1234 mHz. 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.
PDF
