Promise of Persistent Multi-Messenger Astronomy with the Blazar OJ 287

Successful observations of the seven predicted bremsstrahlung flares from the unique bright blazar OJ 287 firmly point to the presence of a nanohertz gravitational wave (GW) emitting supermassive black hole (SMBH) binary central engine. We present arguments for the continued monitoring of the source in several electromagnetic windows to firmly establish various details of the SMBH binary central engine description for OJ 287. In this article, we explore what more can be known about this system, particularly with regard to accretion and outflows from its two accretion disks. We mainly concentrate on the expected impact of the secondary black hole on the disk of the primary on 3 December 2021 and the resulting electromagnetic signals in the following years. We also predict the times of exceptional fades, and outline their usefulness in the study of the host galaxy. A spectral survey has been carried out, and spectral lines from the secondary were searched for but were not found. The jet of the secondary has been studied and proposals to discover it in future VLBI observations are mentioned. In conclusion, the binary black hole model explains a large number of observations of different kinds in OJ 287. Carefully timed future observations will be able to provide further details of its central engine. Such multi-wavelength and multidisciplinary efforts will be required to pursue multi-messenger nanohertz GW astronomy with OJ 287 in the coming decades.

MOMO. IV. The Complete Swift X-Ray and UV/Optical Light Curve and Characteristic Variability of the Blazar OJ 287 during the Last Two Decades

We have been carrying out a dense monitoring of the blazar OJ 287 with Swift since late 2015 as part of our project MOMO (Multiwavelength Observations and Modeling of OJ 287). This is the densest existing monitoring of OJ 287 involving X-ray/UV data. In this latest publication of a sequence, we characterize the multiwavelength variability of OJ 287 based on >4000 Swift single-wave-band data sets including archival data since 2005. A structure function analysis reveals a characteristic timescale of ∼5 days in the optical–UV at epochs of low-level activity and larger during outbursts. The discrete correlation function shows zero lag between optical and UV, with τ = 0 ± 1 day at the epoch of densest cadence. During outbursts (in 2016/17 and 2020) the X-rays follow the UV with near-zero lags. However, during quiescence, the delay is 7–18 days with X-rays leading or lagging, interpreted as due to a different X-ray component dominated by inverse Compton emission. Scaling relations are used to derive the characteristic length scales of the broad-line region and torus in OJ 287. A remarkable, symmetric UV–optical deep fade is identified in late 2017, lasting 2 months. We rule out occultation from the passage of a dusty cloud and a model where the secondary black hole deflects the jet between the primary and observer. We speculate about a temporary dispersion or jet swing event in the core or in a bright quasi-stationary jet feature. The deep fade reveals an additional, spatially distinct X-ray component. The epoch 2020.9–2021.1 was searched for precursor flare activity predicted by the binary black hole model of OJ 287.

The Blazar OJ 287 Jet from Parsec to Kiloparsec Scales

The curved shape of the kiloparsec-scale jet of the blazar OJ 287 is analyzed in the framework of the precession of the central engine, on the existence on which a large number of studies over the past decades are based. The data necessary for the analysis on the kiloparsec-scale jet velocity and angle with the line of sight are obtained based on two competing assumptions about the X-ray emission mechanism of the OJ 287 jet. Namely, there were both the inverse Compton scattering of the cosmic microwave background under the assumption of relativistic kiloparsec-scale jet and the inverse Compton scattering of the central source radiation. For the latter one, we showed that the expected flux from the kiloparsec-scale jet in the gamma range does not exceed the limit set for it according to Fermi-LAT data. We found that only the period of the kiloparsec-scale jet helix, estimated in the framework of the inverse Compton scattering of the central source radiation, agrees with the precession period of the central engine, determined from the modulation of the peak values of 12‑year optical flares.

Explaining temporal variations in the jet position angle of the blazar OJ 287 using its binary black hole central engine model

The bright blazar OJ 287 is the best-known candidate for hosting a supermassive black hole binary system. It inspirals due to the emission of nanohertz gravitational waves (GWs). Observations of historical and predicted quasi-periodic high-brightness flares in its century-long optical lightcurve, allow us to determine the orbital parameters associated with the binary black hole (BBH) central engine. In contrast, the radio jet of OJ 287 has been covered with Very Long Baseline Interferometry (VLBI) observations for only about 30 years and these observations reveal that the position angle (PA) of the jet exhibits temporal variations at both millimetre and centimetre wavelengths. Here we associate the observed PA variations in OJ 287 with the precession of its radio jet. In our model, the evolution of the jet direction can be associated either with the primary black hole (BH) spin evolution or with the precession of the angular momentum direction of the inner region of the accretion disc. Our Bayesian analysis shows that the BBH central engine model, primarily developed from optical observations, can also broadly explain the observed temporal variations in the radio jet of OJ 287 at frequencies of 86, 43, and 15 GHz. Ongoing Global mm-VLBI Array (GMVA) observations of OJ 287 have the potential to verify our predictions for the evolution of its 86 GHz PA values. Additionally, thanks to the extremely high angular resolution that the Event Horizon Telescope (EHT) can provide, we explore the possibility to test our BBH model through the detection of the jet in the secondary black hole .

A possible γ-ray quasi-periodic oscillation of ∼314 days in the blazar OJ 287

ABSTRACT We report the detection of a probable γ-ray quasi-periodic oscillation (QPO) of around 314 d in the monthly binned 0.1–300 GeV γ-ray Fermi-Large Area Telescope light curve of the well-known BL Lacertae blazar OJ 287. To identify and quantify the QPO nature of the γ-ray light curve of OJ 287, we used the Lomb–Scargle periodogram (LSP), REDFIT, and weighted wavelet Z-transform (WWZ) analyses. We briefly discuss possible emission models for radio-loud active galactic nuclei that can explain a γ-ray QPO of such a period in a blazar.

A hadronic emission model for black hole-disc impacts in the blazar OJ 287

ABSTRACT A supermassive black hole (SMBH) binary in the core of the blazar OJ 287 has been invoked in previous works to explain its observed optical flare quasi-periodicity. Following this picture, we investigate a hadronic origin for the X-ray and γ-ray counterparts of the November 2015 major optical flare of this source. An impact outflow must result after the lighter SMBH (the secondary) crosses the accretion disc of the heavier one (the primary). We then consider acceleration of cosmic ray (CR) protons in the shock driven by the impact outflow as it expands and collides with the active galactic nucleus (AGN) wind of the primary SMBH. We show that the emission of these CRs can reproduce the X-ray and γ-ray flare data self-consistently with the optical component of the 2015 November major flare. The derived emission models are consistent with a magnetic field B ∼ 5 G in the emission region and a power-law index of q ∼ 2.2 for the energy distribution of the emitting CRs. The mechanical luminosity of the AGN wind represents $\lesssim 50{{\ \rm per\ cent}}$ of the mass accretion power of the primary SMBH in all the derived emission profiles.