Principles of Gravitational-Wave Detection with Pulsar Timing Arrays

Pulsar timing uses the highly stable pulsar spin period to investigate many astrophysical topics. In particular, pulsar timing arrays make use of a set of extremely well-timed pulsars and their time correlations as a challenging detector of gravitational waves. It turns out that pulsar timing arrays are particularly sensitive to ultra-low-frequency gravitational waves, which makes them complementary to other gravitational-wave detectors. Here, we summarize the basics, focusing especially on supermassive black-hole binaries and cosmic strings, which have the potential to form a stochastic gravitational-wave background in the pulsar timing array detection band, and the scientific goals on this challenging topic. We also briefly outline the recent interesting results of the main pulsar timing array collaborations, which have found strong evidence of a common-spectrum process compatible with a stochastic gravitational-wave background and mention some new perspectives that are particularly interesting in view of the forthcoming radio observatories such as the Five hundred-meter Aperture Spherical Telescope, the MeerKAT telescope, and the Square Kilometer Array.

Mode changing in J1909−3744: the most precisely timed pulsar

We present baseband radio observations of the millisecond pulsar J1909−3744, the most precisely timed pulsar, using the MeerKAT telescope as part of the MeerTime pulsar timing array campaign. During a particularly bright scintillation event the pulsar showed strong evidence of pulse mode changing, among the first millisecond pulsars and the shortest duty cycle millisecond pulsar to do so. Two modes appear to be present, with the weak (lower signal-to-noise ratio) mode arriving 9.26 ±3.94 μs earlier than the strong counterpart. Further, we present a new value of the jitter noise for this pulsar of 8.20 ± 0.14 ns in one hour, finding it to be consistent with previous measurements taken with the MeerKAT (9 ± 3 ns) and Parkes (8.6 ± 0.8 ns) telescopes, but inconsistent with the previously most precise measurement taken with the Green Bank telescope (14 ± 0.5 ns). Timing analysis on the individual modes is carried out for this pulsar, and we find an approximate $10\%$ improvement in the timing precision is achievable through timing the strong mode only as opposed to the full sample of pulses. By forming a model of the average pulse from templates of the two modes, we time them simultaneously and demonstrate that this timing improvement can also be achieved in regular timing observations. We discuss the impact an improvement of this degree on this pulsar would have on searches for the stochastic gravitational wave background, as well as the impact of a similar improvement on all MeerTime PTA pulsars.

Search for Continuous Gravitational-wave Signals in Pulsar Timing Residuals: A New Scalable Approach with Diffusive Nested Sampling

Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov Chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution is to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semianalytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling has shown capability in coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including  ( 10 2 ) pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also highly beneficial to include DNest in current toolboxes of PTA analysis.

Probing dark contents in globular clusters with timing effects of pulsar acceleration

We investigate pulsar timing residuals due to the coupling effect of the pulsar transverse acceleration and the Römer delay. The effect is relatively small and usually negligible. Only for pulsars in globular clusters, it is possibly important. The maximum residual amplitude, which is from the pulsar near the surface of the core of the cluster, is about tens of nanoseconds, and may hardly be identified for most globular clusters currently. However, an intermediate-mass black hole in the center of a cluster can apparently increase the timing residual magnitudes. Particularly for pulsars in the innermost core region, their residual magnitudes may be significant. The high-magnitude residuals, which are above critical lines of each cluster, are strong evidence for the presence of a black hole or dark remnants of comparable total mass in the center of the cluster. We also explored the timing effects of line-of-sight accelerations for the pulsars. The distribution of measured line-of-sight accelerations are simulated with a Monte Carlo method. Two-dimensional Kolmogorov-Smirnov tests are performed to reexamine the consistency of distributions of the simulated and reported data for various values of parameters of the clusters. It is shown that the structure parameters of Terzan 5 can be constrained well by comparing the distribution of measured line-of-sight accelerations with the distributions from Monte Carlo simulations. We find that the cluster has an upper limit on the central black hole/dark remnant mass of ∼ 6000 M ⊙.

Exploring the early Universe with Gaia and Theia

It has recently been pointed out that Gaia is capable of detecting a stochastic gravitational wave background in the sensitivity band between the frequency of pulsar timing arrays and LISA. We argue that Gaia and Theia have great potential for early universe cosmology, since such a frequency range is ideal for probing phase transitions in asymmetric dark matter, SIMP and the cosmological QCD transition. Furthermore, there is the potential for detecting primordial black holes in the solar mass range produced during such an early universe transition and distinguish them from those expected from the QCD epoch. Finally, we discuss the potential for Gaia and Theia to probe topological defects and the ability of Gaia to potentially shed light on the recent NANOGrav results.

Probing Ultralight Bosons with Compact Eccentric Binaries

Ultralight bosons can be abundantly produced through superradiance process by a spinning black hole and form a bound state with hydrogen-like spectrum. We show that such a gravitational atom typically possesses anomalously large mass quadrupole and leads to significant orbital precession when it forms an eccentric binary with a second compact object. Dynamically formed black hole binaries or pulsar-black hole binaries are typically eccentric during their early inspirals. We show that the large orbital precession can generate distinct and observable signature in their gravitational wave or pulsar timing signals.

Gravitational waves and dark photon dark matter from axion rotations

An axion rotating in field space can produce dark photons in the early universe via tachyonic instability. This explosive particle production creates a background of stochastic gravitational waves that may be visible at pulsar timing arrays or other gravitational wave detectors. This scenario provides a novel history for dark photon dark matter. The dark photons may be warm at a level detectable in future 21-cm line surveys. For a consistent cosmology, the radial direction of the complex field containing the axion must be thermalized. We explore a concrete thermalization mechanism in detail and also demonstrate how this setup can be responsible for the generation of the observed baryon asymmetry.

Stochastic gravitational-wave background from metastable cosmic strings

A metastable cosmic-string network is a generic consequence of many grand unified theories (GUTs) when combined with cosmic inflation. Metastable cosmic strings are not topologically stable, but decay on cosmic time scales due to pair production of GUT monopoles. This leads to a network consisting of metastable long strings on superhorizon scales as well as of string loops and segments on subhorizon scales. We compute for the first time the complete stochastic gravitational-wave background (SGWB) arising from all these network constituents, including several technical improvements to both the derivation of the loop and segment contributions. We find that the gravitational waves emitted by string loops provide the main contribution to the gravitational-wave spectrum in the relevant parameter space. The resulting spectrum is consistent with the tentative signal observed by the NANOGrav and Parkes pulsar timing collaborations for a string tension of G μ ∼ 10-11…-7 and has ample discovery space for ground- and space-based detectors. For GUT-scale string tensions, G μ ∼ 10-8…-7, metastable strings predict a SGWB in the LIGO-Virgo-KAGRA band that could be discovered in the near future.