Search

The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravitational waves (GWs). Here we describe the noise models developed for each of the MSPs in the Parkes Pulsar Timing Array (PPTA) third data release, which have been used as the basis of a search for the isotropic stochastic GW background. We model pulsar spin noise, dispersion measure variations, scattering variations, events in the pulsar magnetospheres, solar wind variability, and instrumental effects. We also search for new timing model parameters and detected Shapiro delays in PSR~J0614$-$3329 and PSR~J1902$-$5105. The noise and timing models are validated by testing the normalized and whitened timing residuals for Gaussianity and residual correlations with time. We demonstrate that the choice of noise models significantly affects the inferred properties of a common-spectrum process. Using our detailed models, the recovered common-spectrum noise in the PPTA is consistent with a power law with a spectral index of $\gamma=13/3$, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission., Comment: 18 pages, 10 figures. Accepted for publication in ApJL

Pulsar timing arrays (PTAs) and the Laser Interferometer Space Antenna (LISA) will open complementary observational windows on massive black-hole binaries (MBHBs), i.e., with masses in the range $\sim 10^6 - 10^{10}\,$ M$_{\odot}$. While PTAs may detect a stochastic gravitational-wave background from a population of MBHBs, during operation LISA will detect individual merging MBHBs. To demonstrate the profound interplay between LISA and PTAs, we estimate the number of MBHB mergers that one can expect to observe with LISA by extrapolating direct observational constraints on the MBHB merger rate inferred from PTA data. For this, we postulate that the common noise currently detected in PTAs is an astrophysical background sourced by a single MBHB population. We then constrain the LISA detection rate, $\mathcal{R}$, in the mass-redshift space by combining our Bayesian-inferred merger rate with LISA's sensitivity to spin-aligned, inspiral-merger-ringdown waveforms. Using an astrophysically-informed formation model, we predict a 95$\%$ upper limit on the detection rate of $\mathcal{R} < 134\,{\rm yr}^{-1}$ for binaries with total masses in the range $10^7 - 10^8\,$ M$_{\odot}$. For higher masses, i.e., $>10^8\,$ M$_{\odot}$, we find $\mathcal{R} < 2\,(1)\,\mathrm{yr}^{-1}$ using an astrophysically-informed (agnostic) formation model, rising to $11\,(6)\,\mathrm{yr}^{-1}$ if the LISA sensitivity bandwidth extends down to $10^{-5}$ Hz. Forecasts of LISA science potential with PTA background measurements should improve as PTAs continue their search.

The observation of massive black hole binary systems is one of the main science objectives of the Laser Interferometer Space Antenna (LISA). The instrument's design requirements have recently been revised: they set a requirement at $0.1\,\mathrm{mHz}$, with no additional explicit requirements at lower frequencies. This has implications for observations of the short-lived signals produced by the coalescence of massive and high-redshift binaries. Here we consider the most pessimistic scenario: the (unlikely) case in which LISA has no sensitivity below $0.1\,\mathrm{mHz}$. We show that the presence of higher multipoles (beyond the dominant $\ell = |m| = 2$ mode) in the gravitational radiation from these systems, which will be detectable with a total signal-to-noise ratio $\sim 10^3$, allows LISA to retain the capability to accurately measure the physical parameters, the redshift, and to constrain the sky location. To illustrate this point, we consider a few select binaries in a total (redshifted) mass range of $4 \times10^6 - 4 \times 10^7\,M_\odot$ whose ($\ell = |m| = 2$) gravitational-wave signals last between $\approx 12$ hours and $\approx 20$ days in band. We model the emitted gravitational radiation using the highly accurate (spin-aligned) waveform approximant IMRPhenomXHM and carry out a fully coherent Bayesian analysis on the LISA noise-orthogonal time-delay-interferometry channels., Comment: 13 pages, 6 figures, comments welcome!