Binary black hole (BH) central engine description for the unique blazar OJ 287 predicted that the next secondary BH impact-induced bremsstrahlung flare should peak on 2019 July 31. This prediction was based on detailed general relativistic modeling of the secondary BH trajectory around the primary BH and its accretion disk. The expected flare was termed the Eddington flare to commemorate the centennial celebrations of now-famous solar eclipse observations to test general relativity by Sir Arthur Eddington. We analyze the multi-epoch Spitzer observations of the expected flare between 2019 July 31 and 2019 September 6, as well as baseline observations during 2019 February-March. Observed Spitzer flux density variations during the predicted outburst time display a strong similarity with the observed optical pericenter flare from OJ 287 during 2007 September. The predicted flare appears comparable to the 2007 flare after subtracting the expected higher base-level Spitzer flux densities at 3.55 and 4.49 $\mu$m compared to the optical R-band. Comparing the 2019 and 2007 outburst lightcurves and the previously calculated predictions, we find that the Eddington flare arrived within 4 hours of the predicted time. Our Spitzer observations are well consistent with the presence of a nano-Hertz gravitational wave emitting spinning massive binary BH that inspirals along a general relativistic eccentric orbit in OJ 287. These multi-epoch Spitzer observations provide a parametric constraint on the celebrated BH no-hair theorem.

We report on the results of a multiwavelength monitoring campaign of the bright, nearby Seyfert galaxy, Ark 120 using a ~50-day observing programme with Swift and a ~4-month co-ordinated ground-based observing campaign, predominantly using the Skynet Robotic Telescope Network. We find Ark 120 to be variable at all optical, UV, and X-ray wavelengths, with the variability observed to be well-correlated between wavelength bands on short timescales. We perform cross-correlation analysis across all available wavelength bands, detecting time delays between emission in the X-ray band and the Swift V, B and UVW1 bands. In each case, we find that the longer-wavelength emission is delayed with respect to the shorter-wavelength emission. Within our measurement uncertainties, the time delays are consistent with the \tau ~ \lambda^{4/3} relation, as predicted by a disc reprocessing scenario. The measured lag centroids are \tau_{cent} = 11.90 +/- 7.33, 10.80 +/- 4.08, and 10.60 +/- 2.87 days between the X-ray and V, B, and UVW1 bands, respectively. These time delays are longer than those expected from standard accretion theory and, as such, Ark 120 may be another example of an active galaxy whose accretion disc appears to exist on a larger scale than predicted by the standard thin-disc model. Additionally, we detect further inter-band time delays: most notably between the ground-based I and B bands (\tau_{cent} = 3.46 +/- 0.86 days), and between both the Swift XRT and UVW1 bands and the I band (\tau_{cent} = 12.34 +/- 4.83 and 2.69 +/- 2.05 days, respectively), highlighting the importance of co-ordinated ground-based optical observations.

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes.

Here we explore the disk-jet connection in the broad-line radio quasar 4C+74.26, utilizing the results of the multiwavelength monitoring of the source. The target is unique in that its radiative output at radio wavelengths is dominated by a moderately-beamed nuclear jet, at optical frequencies by the accretion disk, and in the hard X-ray range by the disk corona. Our analysis reveals a correlation (local and global significance of 96\% and 98\%, respectively) between the optical and radio bands, with the disk lagging behind the jet by $250 \pm 42$ days. We discuss the possible explanation for this, speculating that the observed disk and the jet flux changes are generated by magnetic fluctuations originating within the innermost parts of a truncated disk, and that the lag is related to a delayed radiative response of the disk when compared with the propagation timescale of magnetic perturbations along relativistic outflow. This scenario is supported by the re-analysis of the NuSTAR data, modelled in terms of a relativistic reflection from the disk illuminated by the coronal emission, which returns the inner disk radius $R_{\rm in}/R_{\rm ISCO} =35^{+40}_{-16}$. We discuss the global energetics in the system, arguing that while the accretion proceeds at the Eddington rate, with the accretion-related bolometric luminosity $L_{\rm bol} \sim 9 \times 10^{46}$ erg s$^{-1}$ $\sim 0.2 L_{\rm Edd}$, the jet total kinetic energy $L_\textrm{j} \sim 4 \times 10^{44}$ erg s$^{-1}$, inferred from the dynamical modelling of the giant radio lobes in the source, constitutes only a small fraction of the available accretion power.

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.

We report the results of spectroscopic and photometric observations of the emission-line object AS 386. For the first time, we found that it exhibits the B[e] phenomenon and fits the definition of an FS CMa type object. The optical spectrum shows the presence of a B-type star with the following properties: T_ eff = 11000+/-500 K, log L/L_sun = 3.7+/-0.3, a mass of 7+/-1 M_sun, and a distance D = 2.4+/-0.3 kpc from the Sun. We detected regular radial velocity variations of both absorption and emission lines with the following orbital parameters: P_orb = 131.27+/-0.09 days, semi-amplitude K_1 = 51.7+/-3.0 km/s, systemic radial velocity gamma = -31.8+/-2.6 km/s, and a mass function of f(m) = 1.9+/-0.3 M_sun. AS 386 exhibits irregular variations of the optical brightness (V=10.92+/-0.05 mag), while the near-IR brightness varies up to ~0.3 mag following the spectroscopic period. We explain this behavior by a variable illumination of the dusty disk inner rim by the B-type component. Doppler tomography based on the orbital variations of emission-line profiles shows that the material is distributed near the B-type component and in a circumbinary disk. We conclude that the system has undergone a strong mass transfer that created the circumstellar material and increased the B-type component mass. The absence of any traces of a secondary component, whose mass should be >= 7 M_sun, suggests that it is most likely a black hole.

We present results from a search for short time-scale white dwarf variability using \texttt{gPhoton}, a time-tagged database of \textit{GALEX} photon events and associated software package. We conducted a survey of $320$ white dwarf stars in the McCook-Sion catalogue, inspecting each for photometric variability with particular emphasis on variability over time-scales less than $\sim 30$ minutes. From that survey, we present the discovery of a new pulsating white dwarf: WD 2246-069. A Ca II K line is found in archival ESO spectra and an IR excess is seen in WISE $W1$ and $W2$ bands. Its independent modes are identified in follow-up optical photometry and used to model its interior structure. Additionally, we detect UV pulsations in four previously known pulsating ZZ Ceti-type (DAVs). Included in this group is the simultaneous fitting of the pulsations of WD 1401-147 in optical, near-ultraviolet and far-ultraviolet bands using nearly concurrent Whole Earth Telescope and \textit{GALEX} data, providing observational insight into the wavelength dependence of white dwarf pulsation amplitudes.

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.

The eclipsing and strongly interacting binary star system R Arae (HD149730) is in a very active and very short-lived stage of its evolution. R Ara consists of a B9V primary and an unknown secondary. We have collected the International Ultraviolet Explorer (IUE) archival data on R Ara, with most of the data being studied for the first time. There are 117 high resolution IUE spectra taken in 1980, 1982, 1985, 1989, and 1991. We provide photometric and spectroscopic evidence for mass transfer and propose a geometry for the accretion structure. We use colour scale radial velocity plots to view the complicated behavior of the blended absorption features and to distinguish the motions of hotter and cooler regions within the system. We observed a primary eclipse of R Ara in 2008 and have verified that its period is increasing. A model of the system and its evolutionary status is presented.

Spectroscopic measurements of Mira variable stars, as a function of phase, probe the stellar atmospheres and underlying pulsation mechanisms. For example, measuring variations in TiO, VO, and ZrO with phase can be used to help determine whether these molecular species are produced in an extended region above the layers where Balmer line emission occurs or below this shocked region. Using the same methods, the Balmer-line increment, where the strongest Balmer line at phase zero is H-delta and not H-alpha can be measured and explanations tested, along with another peculiarity, the absence of the H-epsilon line in the spectra of Miras when other Balmer lines are strong. We present new spectra covering the spectral range from 6200 Angstroms to 9000 Angstroms of 20 Mira variables. A relationship between variations in the CaII IR triplet and H-alpha as a function of phase support the hypothesis that H-epsilon's observational characteristics result from an interaction of H-epsilon photons with the CaII H line. New periods and epochs of variability are also presented for each star.