In the binary black hole model of OJ 287 the secondary black hole orbits a much more massive primary, and impacts on the primary accretion disk at predictable times. We update the parameters of the disk, the viscosity $\alpha$ and the mass accretion rate $\dot m$. We find $\alpha=0.26 \pm 0.1$ and $\dot m = 0.08 \pm 0.04$ in Eddington units. The former value is consistent with Coroniti (1981) and the latter with Marscher and Jorstad (2011). Predictions are made for the 2019 July 30 superflare in OJ 287. We expect that it will take place simultaneously at the Spitzer infrared channels as well as in the optical and that therefore the timing of the flare in optical can be accurately determined from Spitzer observations. We also discuss in detail the light curve of the 2015 flare and find that the radiating volume has regions where bremsstrahlung dominates as well as regions that radiate primarily in synchrotron radiation. The former region produces the unpolarised first flare while the latter region gives rise to a highly polarized second flare.

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.

The cyclic behaviour of (O-C) residuals of eclipse timings in the sdB+M eclipsing binary NSVS 14256825 was previously attributed to one or two Jovian-type circumbinary planets. We report 83 new eclipse timings that not only fill in the gaps in those already published but also extend the time span of the (O-C) diagram by three years. Based on the archival and our new data spanning over more than 17 years we re-examined the up to date system (O-C). The data revealed systematic, quasi-sinusoidal variation deviating from an older linear ephemeris by about 100 s. It also exhibits a maximum in the (O-C) near JD 2,456,400 that was previously unknown. We consider two most credible explanations of the (O-C) variability: the light propagation time due to the presence of an invisible companion in a distant circumbinary orbit, and magnetic cycles reshaping one of the binary components, known as the Applegate or Lanza-Rodono effect. We found that the latter mechanism is unlikely due to the insufficient energy budget of the M-dwarf secondary. In the framework of the third-body hypothesis, we obtained meaningful constraints on the Keplerian parameters of a putative companion and its mass. Our best-fitting model indicates that the observed quasi-periodic (O-C) variability can be explained by the presence of a brown dwarf with the minimal mass of 15 Jupiter masses rather than a planet, orbiting the binary in a moderately elliptical orbit (~ 0.175) with the period of ~ 10 years. Our analysis rules out two planets model proposed earlier.

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.

We study the mid-egress eclipse timing data gathered for the cataclysmic binary HU Aquarii during the years 1993-2014. The (O-C) residuals were previously attributed to a single ~7 Jupiter mass companion in ~5 au orbit or to a stable 2-planet system with an unconstrained outermost orbit. We present 22 new observations gathered between June, 2011 and July, 2014 with four instruments around the world. They reveal a systematic deviation of ~60 - 120 seconds from the older ephemeris. We re-analyse the whole set of the timing data available. Our results provide an erratum to the previous HU Aqr planetary models, indicating that the hypothesis for a third and fourth body in this system is uncertain. The dynamical stability criterion and a particular geometry of orbits rule out coplanar 2-planet configurations. A putative HU Aqr planetary system may be more complex, e.g., highly non-coplanar. Indeed, we found examples of 3-planet configurations with the middle planet in a retrograde orbit, which are stable for at least 1Gyr, and consistent with the observations. The (O-C) may be also driven by oscillations of the gravitational quadrupole moment of the secondary, as predicted by the Lanza et al. modification of the Applegate mechanism. Further systematic, long-term monitoring of HU Aqr is required to interpret the (O-C) residuals.