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.

Individual supermassive black hole binaries in non-circular orbits are possible nanohertz gravitational wave sources for the rapidly maturing Pulsar Timing Array experiments. We develop an accurate and efficient approach to compute Pulsar Timing Array signals due to gravitational waves from inspiraling supermassive black hole binaries in relativistic eccentric orbits. Our approach employs a Keplerian-type parametric solution to model third post-Newtonian accurate precessing eccentric orbits while a novel semi-analytic prescription is provided to model the effects of quadrupolar order gravitational wave emission. These inputs lead to a semi-analytic prescription to model such signals, induced by non-spinning black hole binaries inspiralling along arbitrary eccentricity orbits. Additionally, we provide a fully analytic prescription to model Pulsar Timing Array signals from black hole binaries inspiraling along moderately eccentric orbits, influenced by Boetzelet al.[Phys. Rev. D 96,044011(2017)]. These approaches are being incorporated into Enterprise and TEMPO2 for searching the presence of such binaries in Pulsar Timing Array datasets.

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.

Short-orbital period small-eccentricity binary pulsars can, in principle, experience substantial advance of periastron. We explore the possibility of measuring this effect by implementing a timing model, ELL1k, in the popular TEMPO2 pulsar timing package. True secular variations in the Laplace-Lagrange parameters, present in our ELL1k model, can lead to measurable timing residuals while pursuing decade-long timing campaigns using the existing ELL1 timing model of Lange et al. (2001), especially for binaries exhibiting significant periastron advance. We also list the main differences between our approach and various implementations of the ELL1 model present in both TEMPO and TEMPO2 packages. Detailed TEMPO2 simulations suggest the possibility of constraining the apsidal motion constant of pulsar companions in certain observed binary pulsars with minuscule eccentricities such as PSR J1719-1438. Fortunately, the ELL1k timing model does not pose any challenges to the on-going Pulsar Timing Array campaigns that routinely employ the ELL1 timing model. Detailed TEMPO2 simulations suggest the possibility of constraining the apsidal motion constant of pulsar companions in certain observed binary pulsars with minuscule eccentricities such as PSR J1719-1438. Fortunately, the ELL1k timing model does not pose any challenges to the on-going Pulsar Timing Array campaigns that routinely employ the ELL1 timing model.