The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational-wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within 1σ. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we “extended” each PTA’s data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings–Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA’s Data Release 3, which will involve not just adding in additional pulsars but also including data from all three PTAs where any given pulsar is timed by more than a single PTA.

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 pulse arrival times and high-precision dispersion measure estimates for 14 millisecond pulsars observed simultaneously in the 300 $-$ 500 MHz and 1260 $-$ 1460 MHz frequency bands using the upgraded Giant Metrewave Radio Telescope. The data spans over a baseline of 3.5 years (2018-2021), and is the first official data release made available by the Indian Pulsar Timing Array collaboration. This data release presents a unique opportunity for investigating the interstellar medium effects at low radio frequencies and their impact on the timing precision of pulsar timing array experiments. In addition to the dispersion measure time series and pulse arrival times obtained using both narrowband and wideband timing techniques, we also present the dispersion measure structure function analysis for selected pulsars. Our ongoing investigations regarding the frequency dependence of dispersion measures have been discussed. Based on the preliminary analysis for five millisecond pulsars, we do not find any conclusive evidence of chromaticity in dispersion measures. Data from regular simultaneous two-frequency observations are presented for the first time in this work. This distinctive feature leads us to the highest precision dispersion measure estimates obtained so far for a subset of our sample. Simultaneous multi-band upgraded Giant Metrewave Radio Telescope observations in 300 $-$ 500 MHz and 1260 $-$ 1460 MHz are crucial for high-precision dispersion measure estimation and for the prospect of expanding the overall frequency coverage upon the combination of data from the various Pulsar Timing Array consortia in the near future. Parts of the data presented in this work are expected to be incorporated into the upcoming third data release of the International Pulsar Timing Array.

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 μ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 hr 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 present the first polarimetric space very long baseline interferometry (VLBI) observations of OJ 287, observed with RadioAstron at 22 GHz during a perigee session on 2014 April 4 and five near-in-time snapshots, together with contemporaneous ground VLBI observations at 15, 43, and 86 GHz. Ground-space fringes were obtained up to a projected baseline of 3.9 Earth diameters during the perigee session, and at a record 15.1 Earth diameters during the snapshot sessions, allowing us to image the innermost jet at an angular resolution of ∼50μ as, the highest ever achieved at 22 GHz for OJ 287. Comparison with ground-based VLBI observations reveals a progressive jet bending with increasing angular resolution that agrees with predictions from a supermassive binary black hole model, although other models cannot be ruled out. Spectral analyses suggest that the VLBI core is dominated by the internal energy of the emitting particles during the onset of a multiwavelength flare, while the parsec-scale jet is consistent with being in equipartition between the particles and magnetic field. Estimated minimum brightness temperatures from the visibility amplitudes show a continued rising trend with projected baseline length up to 10¹³ K, reconciled with the inverse-Compton limit through Doppler boosting for a jet closely oriented to the line of sight. The observed electric vector position angle suggests that the innermost jet has a predominantly toroidal magnetic field, which, together with marginal evidence of a gradient in rotation measure across the jet width, indicates that the VLBI core is threaded by a helical magnetic field, in agreement with jet formation models.

PSR J1713+0747 is one of the most precisely timed pulsars in the international pulsar timing array experiment. This pulsar showed an abrupt profile shape change between 2021 April 16, (MJD 59320) and 2021 April 17 (MJD 59321). In this paper, we report the results from multi-frequency observations of this pulsar carried out with the upgraded Giant Metrewave Radio Telescope (uGMRT) before and after the event. We demonstrate the profile change seen in Band 5 (1260 MHz–1460 MHz) and Band 3 (300 MHz–500 MHz). The timing analysis of this pulsar shows a disturbance accompanying this profile change followed by a recovery with a time-scale of ∼159 days. Our data suggest that a model with chromatic index as a free parameter is preferred over models with combinations of achromaticity with DM bump or scattering bump. We determine the frequency dependence to be ∼ν^+1.34.