We introduce pinta, a pipeline for reducing the uGMRT raw pulsar timing data developed for the Indian Pulsar Timing Array experiment. We provide a detailed description of the workflow and usage of pinta, as well as its computational performance and RFI mitigation characteristics. Furthermore, the results of a calibration experiment carried out to determine the relative time offsets between different back-end modes and the correct interpretation of the observatory frequency settings at the uGMRT, which are crucial for performing precision pulsar timing, are also discussed.

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.

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.

The bright blazar OJ 287 is the best-known candidate for hosting a nanohertz gravitational wave (GW) emitting supermassive binary black hole (SMBBH) in the present observable universe. The binary black hole (BBH) central engine model, proposed by Lehto and Valtonen in 1996, was influenced by the two distinct periodicities inferred from the optical light curve of OJ 287. The current improved model employs an accurate general relativistic description to track the trajectory of the secondary black hole (BH) which is crucial to predict the inherent impact flares of OJ 287. The successful observations of three predicted impact flares open up the possibility of using this BBH system to test general relativity in a hitherto unexplored strong field regime. Additionally, we briefly describe an on-going effort to interpret observations of OJ 287 in a Bayesian framework.

Blazar OJ287 exhibits large thermal flares at least twice every 12 years. The times of these flares have been predicted successfully using the model of a quasi-Keplerian eccentric black hole binary where the secondary impacts the accretion disk of the primary, creating the thermal flares. New measurements of the historical light curve have been combined with the observations of the 2015 November/December flare to identify the impact record since year 1886, and to constrain the orbit of the binary. The orbit solution shows that the binary period, now 12.062 year, is decreasing at the rate of 36 days per century. This corresponds to an energy loss to gravitational waves that is 6.5 +- 4 % less than the rate predicted by the standard quadrupolar gravitational wave (GW) emission. We show that the difference is due to higher order gravitational radiation reaction terms that include the dominant order tail contributions.

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 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 binary black hole 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 binary black hole 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.