Search

We present Vela, an efficient, modular, easy-to-use Bayesian pulsar timing and noise analysis package written in Julia. Vela provides an independent, efficient, and parallelized implementation of the full non-linear pulsar timing and noise model along with a Python binding named pyvela. One-time operations such as data file input, clock corrections, and solar system ephemeris computations are performed by pyvela with the help of the PINT pulsar timing package. Its reliability is ensured via careful design utilizing Julia's type system, strict version control, and an exhaustive test suite. This paper describes the design and usage of Vela focusing on the narrowband paradigm., Comment: Accepted for publication in ApJ

Pulse profile stability is a central assumption of standard pulsar timing methods. Thus, it is important for pulsar timing array experiments such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) to account for any pulse profile variability present in their data sets. We show that in the NANOGrav 15-yr data set, the integrated pulse profile of PSR J1022+1001 as seen by the Arecibo radio telescope at 430, 1380, and 2030 MHz varies considerably in its shape from observation to observation. We investigate the possibility that this is due to the "ideal feed assumption" (IFA), on which NANOGrav's routine polarization calibration procedure relies. PSR J1022+1001 is $\sim 90\%$ polarized in one pulse profile component, and also has significant levels of circular polarization. Time-dependent deviations in the feed's polarimetric response (PR) could cause mixing between the intensity I and the other Stokes parameters, leading to the observed variability. We calibrate the PR using a mixture of Measurement Equation Modeling and Measurement Equation Template Matching techniques. The resulting profiles are no less variable than those calibrated using the IFA method, nor do they provide an improvement in the timing quality of this pulsar. We observe the pulse shape in 25-MHz bandwidths to vary consistently across the band, which cannot be explained by interstellar scintillation in combination with profile evolution with frequency. Instead, we favor phenomena intrinsic to the pulsar as the cause., Comment: 18 pages, 16 figures, 4 tables. Submitted to ApJ

Evidence has emerged for a stochastic signal correlated among 67 pulsars within the 15-year pulsar-timing data set compiled by the NANOGrav collaboration. Similar signals have been found in data from the European, Indian, Parkes, and Chinese PTAs. This signal has been interpreted as indicative of the presence of a nanohertz stochastic gravitational wave background. To explore the internal consistency of this result we investigate how the recovered signal strength changes as we remove the pulsars one by one from the data set. We calculate the signal strength using the (noise-marginalized) optimal statistic, a frequentist metric designed to measure correlated excess power in the residuals of the arrival times of the radio pulses. We identify several features emerging from this analysis that were initially unexpected. The significance of these features, however, can only be assessed by comparing the real data to synthetic data sets. After conducting identical analyses on simulated data sets, we do not find anything inconsistent with the presence of a stochastic gravitational wave background in the NANOGrav 15-year data. The methodologies developed here can offer additional tools for application to future, more sensitive data sets. While this analysis provides an internal consistency check of the NANOGrav results, it does not eliminate the necessity for additional investigations that could identify potential systematics or uncover unmodeled physical phenomena in the data., Comment: 21 pages, 11 figures, 2 tables

Pulsar timing array observations have found evidence for an isotropic gravitational wave background with the Hellings-Downs angular correlations, expected from general relativity. This interpretation hinges on the measured shape of the angular correlations, which is predominately quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlations into a sum of Legendre polynomials and use a Bayesian analysis to constrain their coefficients with the 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). When including Legendre polynomials with multipoles $\ell \geq 2$, we only find a significant signal in the quadrupole with an amplitude consistent with general relativity and non-zero at the $\sim 95\%$ confidence level and a Bayes factor of 200. When we include multipoles $\ell \leq 1$, the Bayes factor evidence for quadrupole correlations decreases by more than an order of magnitude due to evidence for a monopolar signal at approximately 4 nHz which has also been noted in previous analyses of the NANOGrav 15-year data. Further work needs to be done in order to better characterize the properties of this monopolar signal and its effect on the evidence for quadrupolar angular correlations., Comment: 15 pages, 6 figures

The detection of a stochastic gravitational wave background by pulsar timing arrays suggests the presence of a supermassive black hole binary population. Although the observed spectrum generally aligns with predictions from orbital evolution driven by gravitational wave emission in circular orbits, there is a discernible preference for a turnover at the lowest observed frequencies. This turnover could indicate a significant hardening phase, transitioning from early environmental influences to later stages predominantly influenced by gravitational wave emission. In the vicinity of these binaries, the ejection of stars or dark matter particles through gravitational three-body slingshots efficiently extracts orbital energy, leading to a low-frequency turnover in the spectrum. By analyzing the NANOGrav 15-year data, we assess how the gravitational wave spectrum depends on the initial inner galactic profile prior to disruption by binary ejections, accounting for a range of initial binary eccentricities. Our findings suggest a parsec-scale galactic center density around $10^6\,M_\odot/\textrm{pc}^3$ across most of the parameter space, offering insights into the environmental effects on black hole evolution and combined matter density near galaxy centers., Comment: 15 pages, 5 figures

The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index \gamma. In this paper, we consider the next logical step beyond this minimal spectral model, allowing for a running (i.e., logarithmic frequency dependence) of the spectral index, \gamma_run(f) = \gamma + \beta \ln(f/f_ref). We fit this running-power-law (RPL) model to the NANOGrav 15-year data and perform a Bayesian model comparison with the minimal constant-power-law (CPL) model, which results in a 95% credible interval for the parameter \beta consistent with no running, \beta \in [-0.80,2.96], and an inconclusive Bayes factor, B(RPL vs. CPL) = 0.69 +- 0.01. We thus conclude that, at present, the minimal CPL model still suffices to adequately describe the NANOGrav signal; however, future data sets may well lead to a measurement of nonzero \beta. Finally, we interpret the RPL model as a description of primordial GWs generated during cosmic inflation, which allows us to combine our results with upper limits from big-bang nucleosynthesis, the cosmic microwave background, and LIGO-Virgo-KAGRA., Comment: 17 pages, 4 figures, 2 tables

Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Those analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residuals across pulsars. These assumptions are idealizations that may not be realized in actuality. We test them in the NANOGrav 15 yr data set using Bayesian posterior predictive checks: after fitting our fiducial model to real data, we generate a population of simulated data-set replications, and use them to assess whether the optimal-statistic significance, inter-pulsar correlations, and spectral coefficients assume extreme values for the real data when compared to the replications. We confirm that the NANOGrav 15 yr data set is consistent with power-law and Hellings--Downs assumptions. We also evaluate the evidence for the stochastic background using posterior-predictive versions of the frequentist optimal statistic and of Bayesian model comparison, and find comparable significance (3.2\ $\sigma$ and 3\ $\sigma$ respectively) to what was previously reported for the standard statistics. We conclude with novel visualizations of the reconstructed gravitational waveforms that enter the residuals for each pulsar. Our analysis strengthens confidence in the identification and characterization of the gravitational-wave background as reported by NANOGrav., Comment: 20 pages, 14 Figures

Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4\sigma$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in \textit{Gaia} DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6\%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 80 kyr, a companion mass of 0.3 $M_\odot$ (indicative of a white dwarf or low-mass star), and a semi-major axis of 2000 AU. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios.

Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect achromatic noise properties in each pulsar. The models we compare include existing DM models used by NANOGrav and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar dataset. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from $\log_{10}A_{\text{RN}} = -14.1^{+0.1}_{-0.1}$ to $-14.7^{+0.3}_{-0.5}$, and the spectral index broadens from $\gamma_{\text{RN}} = 2.6^{+0.5}_{-0.4}$ to $\gamma_{\text{RN}} = 3.5^{+1.2}_{-0.9}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs.

PINT is a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework within PINT to characterize the single-pulsar noise processes present in pulsar timing datasets. This framework enables the parameter estimation for both uncorrelated and correlated noise processes as well as the model comparison between different timing and noise models in a computationally inexpensive way. We demonstrate the efficacy of the new framework by applying it to simulated datasets as well as a real dataset of PSR B1855+09. We also describe the new features implemented in PINT since it was first described in the literature., Comment: Accepted for publication in ApJ

The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 \sigma - 1.9 \sigma$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 \sigma - 2.1 \sigma$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26~\mathrm{nHz}$ break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe., Comment: 13 pages, 8 figures, 2 appendices, in press at ApJ

Burst with memory events are potential transient gravitational wave sources for the maturing pulsar timing array (PTA) efforts. We provide a computationally efficient prescription to model pulsar timing residuals induced by supermassive black hole pairs in general relativistic hyperbolic trajectories employing a Keplerian-type parametric solution. Injection studies have been pursued on the resulting bursts with linear GW memory (LGWM) events with simulated datasets to test the performance of our pipeline, followed by its application to the publicly available NANOGrav 12.5-year (NG12.5) dataset. Given the absence of any evidence of LGWM events within the real NG12.5 dataset, we impose $95\%$ upper limits on the PTA signal amplitude as a function of the sky location of the source and certain characteristic frequency ($n$) of the signal. The upper limits are computed using a signal model that takes into account the presence of intrinsic timing noise specific to each pulsar, as well as a common, spatially uncorrelated red noise, alongside the LGWM signal. Our investigations reveal that the $95\%$ upper limits on LGWM amplitude, marginalized over all other parameters, is 3.48 $\pm 0.51 \ \mu$s for $n>3.16$ nHz. This effort should be relevant for constraining both burst and memory events in the upcoming International Pulsar Timing Array data releases., Comment: 20 pages, 11 figures

Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis., Comment: 11 pages, 5 figures

The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity $e_0<0.5$ and the symmetric mass ratio $\eta>0.1$. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of $88.1\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $\eta>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations., Comment: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

A pulsar's pulse profile gets broadened at low frequencies due to dispersion along the line of sight or due to multi-path propagation. The dynamic nature of the interstellar medium makes both of these effects time-dependent and introduces slowly varying time delays in the measured times-of-arrival similar to those introduced by passing gravitational waves. In this article, we present an improved method to correct for such delays by obtaining unbiased dispersion measure (DM) measurements by using low-frequency estimates of the scattering parameters. We evaluate this method by comparing the obtained DM estimates with those, where scatter-broadening is ignored using simulated data. A bias is seen in the estimated DMs for simulated data with pulse-broadening with a larger variability for a data set with a variable frequency scaling index, $\alpha$, as compared to that assuming a Kolmogorov turbulence. Application of the proposed method removes this bias robustly for data with band averaged signal-to-noise ratio larger than 100. We report the measurements of the scatter-broadening time and $\alpha$ from analysis of PSR J1643$-$1224, observed with upgraded Giant Metrewave Radio Telescope as part of the Indian Pulsar Timing Array experiment. These scattering parameters were found to vary with epoch and $\alpha$ was different from that expected for Kolmogorov turbulence. Finally, we present the DM time-series after application of this technique to PSR J1643$-$1224., Comment: 10 pages, 8 figures, Accepted for publication in MNRAS

Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them., Comment: 14 pages, 8 figures, version matching published paper

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\sigma$. 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 and 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 as single PTA., Comment: 21 pages, 9 figures, submitted to ApJ

The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy., Comment: 30 pages, 23 figures, replaced to match the version published in Astronomy & Astrophysics, note the change in the numbering order in the series (now paper IV)

We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favor of the presence of the CGW over a common uncorrelated noise process, while the frequentist analysis estimates the p-value of the candidate to be 1%, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only 0.7. Therefore, we cannot conclusively determine the origin of the observed feature, but we cannot rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations of the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature., Comment: 13 figures, 15 pages, accepted

The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the subtle imprints that the GWB induces on pulsar timing data are obscured by many sources of noise. These must be carefully characterized to increase the sensitivity to the GWB. In this paper, we present a novel technique to estimate the optimal number of frequency coefficients for modelling achromatic and chromatic noise and perform model selection. We also incorporate a new model to fit for scattering variations in the pulsar timing package temponest and created realistic simulations of the European Pulsar Timing Array (EPTA) datasets that allowed us to test the efficacy of our noise modelling algorithms. We present an in-depth analysis of the noise properties of 25 millisecond pulsars (MSPs) that form the second data release (DR2) of the EPTA and investigate the effect of incorporating low-frequency data from the Indian PTA collaboration. We use enterprise and temponest packages to compare noise models with those reported with the EPTA DR1. We find that, while in some pulsars we can successfully disentangle chromatic from achromatic noise owing to the wider frequency coverage in DR2, in others the noise models evolve in a more complicated way. We also find evidence of long-term scattering variations in PSR J1600$-$3053. Through our simulations, we identify intrinsic biases in our current noise analysis techniques and discuss their effect on GWB searches. The results presented here directly help improve sensitivity to the GWB and are already being used as part of global PTA efforts., Comment: 20 pages, 6 figures, 9 tables

This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data., Comment: 30 pages, 10 figures, 1 table; Companion paper to "The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background"; For questions or comments, please email comments@nanograv.org

Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of $\sim$4 nHz, and another at $\sim$170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings-Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems., Comment: 23 pages, 13 figures, 2 tables. Accepted for publication in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data., Comment: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations., Comment: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons., Comment: 74 pages, 31 figures, 4 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white noise parameters and two red noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of $7\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals., Comment: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by three years, now spanning nearly 16 years for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic gravitational-wave background., Comment: 90 pages, 74 figures, 6 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on modern observing systems, (iii) the combination of the full data set with the first data release of the InPTA for ten commonly timed millisecond pulsars, and (iv) the combination of the 10.3-year subset with the InPTA data. These combinations allowed us to probe the contributions of instrumental noise and interstellar propagation effects. With the full data set, we find marginal evidence for a GWB, with a Bayes factor of four and a false alarm probability of $4\%$. With the 10.3-year subset, we report evidence for a GWB, with a Bayes factor of $60$ and a false alarm probability of about $0.1\%$ ($\gtrsim 3\sigma$ significance). The addition of the InPTA data yields results that are broadly consistent with the EPTA-only data sets, with the benefit of better noise modelling. Analyses were performed with different data processing pipelines to test the consistency of the results from independent software packages. The inferred spectrum from the latest EPTA data from new generation observing systems is rather uncertain and in mild tension with the common signal measured in the full data set. However, if the spectral index is fixed at 13/3, the two data sets give a similar amplitude of ($2.5\pm0.7)\times10^{-15}$ at a reference frequency of $1\,{\rm yr}^{-1}$. By continuing our detection efforts as part of the International Pulsar Timing Array (IPTA), we expect to be able to improve the measurement of spatial correlations and better characterise this signal in the coming years., Comment: 21 pages, 14 figures, 4 appendix figures, accepted for publication in A&A

We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3\sigma$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4\sigma$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal., Comment: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

Close hyperbolic encounters of black holes (BHs) generate certain Burst With Memory (BWM) events in the frequency windows of the operational, planned, and proposed gravitational wave (GW) observatories. We present detailed explorations of the detectable parameter space of such events that are relevant for the LIGO-Virgo-KAGRA and the International Pulsar Timing Array (IPTA) consortia. The underlying temporally evolving GW polarization states are adapted from Cho et al. [Phys. Rev. D 98, 024039 (2018)] and therefore incorporate general relativistic effects up to the third post-Newtonian order. Further, we provide a prescription to ensure the validity of our waveform family while describing close encounters. Preliminary investigations reveal that optimally placed BWM events should be visible to megaparsec distances for the existing ground-based observatories. In contrast, maturing IPTA datasets should be able to provide constraints on the occurrences of such hyperbolic encounters of supermassive BHs to gigaparsec distances., Comment: 19 pages, 11 figures, accepted for publication in Phys. Rev. D

The wideband timing technique enables the high-precision simultaneous estimation of pulsar Times of Arrival (ToAs) and Dispersion Measures (DMs) while effectively modeling frequency-dependent profile evolution. We present two novel independent methods that extend the standard wideband technique to handle simultaneous multi-band pulsar data incorporating profile evolution over a larger frequency span to estimate DMs and ToAs with enhanced precision. We implement the wideband likelihood using the libstempo python interface to perform wideband timing in the tempo2 framework. We present the application of these techniques to the dataset of fourteen millisecond pulsars observed simultaneously in Band 3 (300 - 500 MHz) and Band 5 (1260 - 1460 MHz) of the upgraded Giant Metrewave Radio Telescope (uGMRT) with a large band gap of 760 MHz as a part of the Indian Pulsar Timing Array (InPTA) campaign. We achieve increased ToA and DM precision and sub-microsecond root mean square post-fit timing residuals by combining simultaneous multi-band pulsar observations done in non-contiguous bands for the first time using our novel techniques., Comment: Published in MNRAS

The Indian Pulsar Timing Array (InPTA) collaboration has recently made its first official data release (DR1) for a sample of 14 pulsars using 3.5 years of uGMRT observations. We present the results of single-pulsar noise analysis for each of these 14 pulsars using the InPTA DR1. For this purpose, we consider white noise, achromatic red noise, dispersion measure (DM) variations, and scattering variations in our analysis. We apply Bayesian model selection to obtain the preferred noise models among these for each pulsar. For PSR J1600$-$3053, we find no evidence of DM and scattering variations, while for PSR J1909$-$3744, we find no significant scattering variations. Properties vary dramatically among pulsars. For example, we find a strong chromatic noise with chromatic index $\sim$ 2.9 for PSR J1939+2134, indicating the possibility of a scattering index that doesn't agree with that expected for a Kolmogorov scattering medium consistent with similar results for millisecond pulsars in past studies. Despite the relatively short time baseline, the noise models broadly agree with the other PTAs and provide, at the same time, well-constrained DM and scattering variations., Comment: Accepted for publication in PRD, 30 pages, 17 figures, 4 tables

Pulsar Timing Array (PTA) experiments are expected to be sensitive to gravitational waves (GWs) emitted by individual supermassive black hole binaries (SMBHBs) inspiralling along eccentric orbits. We compare the computational cost of different methods of computing the PTA signals induced by relativistic eccentric SMBHBs, namely approximate analytic expressions, Fourier series expansion, post-circular expansion, and numerical integration. We show that the fastest method for evaluating PTA signals is by using the approximate analytic expressions, which provides up to ~50 times improvement in computational speed over the alternative methods. We investigate the accuracy of the approximate analytic expressions by employing a mismatch metric valid for PTA signals. We show that this method is accurate within the region of the binary parameter space that is of interest to PTA experiments. We introduce a spline-based method to further accelerate the PTA signal evaluations for narrowband PTA datasets. The efficient methods for computing the eccentric SMBHB-induced PTA signals were implemented in the GWecc.jl package and can be readily accessed from the popular ENTERPRISE package to search for such signals in PTA datasets. Further, we simplify the eccentric SMBHB PTA signal expression for the case of a single-pulsar search and demonstrate our computationally efficient methods by performing a single-pulsar search in the 12.5-year NANOGrav narrowband dataset of PSR J1909-3744 using the simplified expression. These results will be crucial for searching for eccentric SMBHBs in large PTA datasets., Comment: Accepted for publication in Classical and Quantum Gravity

Decades long monitoring of millisecond pulsars, which exhibit highly stable rotational periods, in pulsar timing array experiments is on the threshold of discovering nanohertz stochastic gravitational wave background. This paper describes the Indian Pulsar timing array (InPTA) experiment, which employs the upgraded Giant Metrewave Radio Telescope (uGMRT) for timing an ensemble of millisecond pulsars for this purpose. We highlight InPTA's observation strategies and analysis methods, which are relevant for a future PTA experiment with the more sensitive Square Kilometer Array (SKA) telescope. We show that the unique multi-sub-array multi-band wide-bandwidth frequency coverage of the InPTA provides Dispersion Measure estimates with unprecedented precision for PTA pulsars, e.g., ~ 2 x 10{-5} pc-cm{-3} for PSR J1909-3744. Configuring the SKA-low and SKA-mid as two and four sub-arrays respectively, it is shown that comparable precision is achievable, using observation strategies similar to those pursued by the InPTA, for a larger sample of 62 pulsars requiring about 26 and 7 hours per epoch for the SKA-mid and the SKA-low telescopes respectively. We also review the ongoing efforts to develop PTA-relevant general relativistic constructs that will be required to search for nanohertz gravitational waves from isolated super-massive black hole binary systems like blazar OJ 287. These efforts should be relevant to pursue persistent multi-messenger gravitational wave astronomy during the forthcoming era of the SKA telescope, the Thirty Meter Telescope, and the next-generation Event Horizon Telescope., Comment: Accepted for publication in Journal of Astronomy and Astrophysics for Special Issue on Indian Participation in the SKA (Editors : Abhirup Datta, Nirupam Roy, Preeti Kharb and Tirthankar Roy Choudhury)

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 (uGMRT). 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 uGMRT observations in Band 3 and Band 5 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., Comment: 23 pages, 21 figures, 3 tables. Published in PASA

High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the Indian Pulsar Timing Array (InPTA) campaign. We have used the observations with the 300-500 MHz band of the uGMRT for this purpose. We obtain high precision in DM measurements with precisions of the order 10^{-6}cm^{-3}pc. The ToAs obtained have sub-{\mu}s precision and the root-mean-square of the post-fit ToA residuals are in the sub-{\mu}s range. We find that the uncertainties in the DMs and ToAs obtained with this wideband technique, applied to low-frequency data, are consistent with the results obtained with traditional pulsar timing techniques and comparable to high-frequency results from other PTAs. This work opens up an interesting possibility of using low-frequency wideband observations for precision pulsar timing and gravitational wave detection with similar precision as high-frequency observations used conventionally., Comment: Accepted for publication in MNRAS

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 April 16, 2021 (MJD 59320) and April 17, 2021 (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 timescale of $\sim 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 $\sim\nu^{+1.34}$., Comment: Accepted for publication in MNRAS-Letters

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., Comment: 12 pages, 7 figures, 1 table. Accepted for publication in MNRAS

Pulsar radio emission undergoes dispersion due to the presence of free electrons in the interstellar medium (ISM). The dispersive delay in the arrival time of pulsar signal changes over time due to the varying ISM electron column density along the line of sight. Correcting for this delay accurately is crucial for the detection of nanohertz gravitational waves using Pulsar Timing Arrays. In this work, we present in-band and inter-band DM estimates of four pulsars observed with uGMRT over the timescale of a year using two different template alignment methods. The DMs obtained using both these methods show only subtle differences for PSR 1713+0747 and J1909$-$3744. A considerable offset is seen in the DM of PSR J1939+2134 and J2145$-$0750 between the two methods. This could be due to the presence of scattering in the former and profile evolution in the latter. We find that both methods are useful but could have a systematic offset between the DMs obtained. Irrespective of the template alignment methods followed, the precision on the DMs obtained is about $10^{-3}$ pc cm$^{-3}$ using only BAND3 and $10^{-4}$ pc cm$^{-3}$ after combining data from BAND3 and BAND5 of the uGMRT. In a particular result, we have detected a DM excess of about $5\times10^{-3}$ pc cm$^{-3}$ on 24 February 2019 for PSR J2145$-$0750. This excess appears to be due to the interaction region created by fast solar wind from a coronal hole and a coronal mass ejection (CME) observed from the Sun on that epoch. A detailed analysis of this interesting event is presented., Comment: 11 pages, 6 figures, 2 tables. Accepted by A&A

We introduce pinta, a pipeline for reducing the upgraded Giant Metre-wave Radio Telescope (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. We also discuss a novel and independent determination of the relative time offsets between the different back-end modes of uGMRT and the interpretation of the uGMRT observation frequency settings, and their agreement with results obtained from engineering tests. Further, we demonstrate the capability of pinta to generate data products which can produce high-precision TOAs using PSR J1909-3744 as an example. These results are crucial for performing precision pulsar timing with the uGMRT.

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., Comment: Invited review article; 16 pages, 4 figures, 2 tables, to appear in Universe

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.

We provide an elegant way of solving analytically the third post-Newtonian (3PN) accurate Kepler equation, associated with the 3PN-accurate generalized quasi-Keplerian parametrization for compact binaries in eccentric orbits. An additional analytic solution is presented to check the correctness of our compact solution and we perform comparisons between our PN-accurate analytic solution and a very accurate numerical solution of the PN-accurate Kepler equation. We adapt our approach to compute crucial 3PN-accurate inputs that will be required to compute analytically both the time and frequency domain ready-to-use amplitude-corrected PN-accurate search templates for compact binaries in inspiralling eccentric orbits., Comment: 20 pages, 4 figures