Fast Parameter Inference on Pulsar Timing Arrays with Normalizing Flows

Pulsar timing arrays perform Bayesian posterior inference with expensive Markov chain Monte Carlo (MCMC) methods. Given a dataset of ∼10-100 pulsars and O(10^{3}) timing residuals each, producing a posterior distribution for the stochastic gravitational wave background (SGWB) can take days to a week. The computational bottleneck arises because the likelihood evaluation required for MCMC is extremely costly when considering the dimensionality of the search space. Fortunately, generating simulated data is fast, so modern simulation-based inference techniques can be brought to bear on the problem. In this Letter, we demonstrate how conditional normalizing flows trained on simulated data can be used for extremely fast and accurate estimation of the SGWB posteriors, reducing the sampling time from weeks to a matter of seconds.

A gamma-ray pulsar timing array constrains the nanohertz gravitational wave background

After large galaxies merge, their central supermassive black holes are expected to form binary systems. Their orbital motion should generate a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95% credible limit on the GWB characteristic strain of 1.0 × 10^-14 at a frequency of 1 yr^-1. The sensitivity is expected to scale with t _obs, the observing time span, as [Formula: see text]. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.

Neutron Stars with Baryon Number Violation, Probing Dark Sectors

The neutron lifetime anomaly has been used to motivate the introduction of new physics with hidden-sector particles coupled to baryon number, and on which neutron stars provide powerful constraints. Although the neutron lifetime anomaly may eventually prove to be of mundane origin, we use it as motivation for a broader review of the ways that baryon number violation, be it real or apparent, and dark sectors can intertwine and how neutron star observables, both present and future, can constrain them.

Constraining Cosmological Phase Transitions with the Parkes Pulsar Timing Array

A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave background with the latest data set of the Parkes Pulsar Timing Array. We find no evidence for a Hellings-Downs spatial correlation as expected for a stochastic gravitational wave background. Therefore, we present constraints on first-order phase transition model parameters. Our analysis shows that pulsar timing is particularly sensitive to the low-temperature (T∼1-100  MeV) phase transition with a duration (β/H_{*})^{-1}∼10^{-2}-10^{-1} and therefore can be used to constrain the dark and QCD phase transitions.

Searching for Gravitational Waves from Cosmological Phase Transitions with the NANOGrav 12.5-Year Dataset

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5-year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black hole mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.