THE NANOGRAV NINE-YEAR DATA SET: OBSERVATIONS, ARRIVAL TIME MEASUREMENTS, AND ANALYSIS OF 37 MILLISECOND PULSARS

First-Order Phase Transition Interpretation of Pulsar Timing Array Signal Is Consistent with Solar-Mass Black Holes

We perform a Bayesian analysis of NANOGrav 15-yr and IPTA DR2 pulsar timing residuals and show that the recently detected stochastic gravitational-wave background is compatible with a stochastic gravitational-wave background produced by bubble dynamics during a cosmological first-order phase transition. The timing data suggest that the phase transition would occur around QCD confinement temperature and would have a slow rate of completion. This scenario can naturally lead to the abundant production of primordial black holes with solar masses. These primordial black holes can potentially be detected by current and advanced gravitational-wave detectors LIGO-Virgo-Kagra, Einstein Telescope, Cosmic Explorer, by astrometry with GAIA, and by 21-cm survey.

Stochastic Gravitational-Wave Backgrounds: Current Detection Efforts and Future Prospects

The collection of individually resolvable gravitational wave (GW) events makes up a tiny fraction of all GW signals that reach our detectors, while most lie below the confusion limit and are undetected. Similarly to voices in a crowded room, the collection of unresolved signals gives rise to a background that is well-described via stochastic variables and, hence, referred to as the stochastic GW background (SGWB). In this review, we provide an overview of stochastic GW signals and characterise them based on features of interest such as generation processes and observational properties. We then review the current detection strategies for stochastic backgrounds, offering a ready-to-use manual for stochastic GW searches in real data. In the process, we distinguish between interferometric measurements of GWs, either by ground-based or space-based laser interferometers, and timing-residuals analyses with pulsar timing arrays (PTAs). These detection methods have been applied to real data both by large GW collaborations and smaller research groups, and the most recent and instructive results are reported here. We close this review with an outlook on future observations with third generation detectors, space-based interferometers, and potential noninterferometric detection methods proposed in the literature.

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.

Point-to-point stabilized optical frequency transfer with active optics

Timescale comparison between optical atomic clocks over ground-to-space and terrestrial free-space laser links will have enormous benefits for fundamental and applied sciences. However, atmospheric turbulence creates phase noise and beam wander that degrade the measurement precision. Here we report on phase-stabilized optical frequency transfer over a 265 m horizontal point-to-point free-space link between optical terminals with active tip-tilt mirrors to suppress beam wander, in a compact, human-portable set-up. A phase-stabilized 715 m underground optical fiber link between the two terminals is used to measure the performance of the free-space link. The active optical terminals enable continuous, cycle-slip free, coherent transmission over periods longer than an hour. In this work, we achieve residual instabilities of 2.7 × 10-6 rad2 Hz-1 at 1 Hz in phase, and 1.6 × 10-19 at 40 s of integration in fractional frequency; this performance surpasses the best optical atomic clocks, ensuring clock-limited frequency comparison over turbulent free-space links.

Universality of free fall from the orbital motion of a pulsar in a stellar triple system

Einstein's theory of gravity-the general theory of relativity¹-is based on the universality of free fall, which specifies that all objects accelerate identically in an external gravitational field. In contrast to almost all alternative theories of gravity², the strong equivalence principle of general relativity requires universality of free fall to apply even to bodies with strong self-gravity. Direct tests of this principle using Solar System bodies^3,4 are limited by the weak self-gravity of the bodies, and tests using pulsar-white-dwarf binaries^5,6 have been limited by the weak gravitational pull of the Milky Way. PSR J0337+1715 is a hierarchical system of three stars (a stellar triple system) in which a binary consisting of a millisecond radio pulsar and a white dwarf in a 1.6-day orbit is itself in a 327-day orbit with another white dwarf. This system permits a test that compares how the gravitational pull of the outer white dwarf affects the pulsar, which has strong self-gravity, and the inner white dwarf. Here we report that the accelerations of the pulsar and its nearby white-dwarf companion differ fractionally by no more than 2.6 × 10^-6. For a rough comparison, our limit on the strong-field Nordtvedt parameter, which measures violation of the universality of free fall, is a factor of ten smaller than that obtained from (weak-field) Solar System tests^3,4 and a factor of almost a thousand smaller than that obtained from other strong-field tests^5,6.

Observations of the missing baryons in the warm-hot intergalactic medium

It has been known for decades that the observed number of baryons in the local Universe falls about 30-40 per cent short^1,2 of the total number of baryons predicted ³ by Big Bang nucleosynthesis, as inferred^4,5 from density fluctuations of the cosmic microwave background and seen during the first 2-3 billion years of the Universe in the so-called 'Lyman α forest'^6,7 (a dense series of intervening H I Lyman α absorption lines in the optical spectra of background quasars). A theoretical solution to this paradox locates the missing baryons in the hot and tenuous filamentary gas between galaxies, known as the warm-hot intergalactic medium. However, it is difficult to detect them there because the largest by far constituent of this gas-hydrogen-is mostly ionized and therefore almost invisible in far-ultraviolet spectra with typical signal-to-noise ratios^8,9. Indeed, despite large observational efforts, only a few marginal claims of detection have been made so far^2,10. Here we report observations of two absorbers of highly ionized oxygen (O VII) in the high-signal-to-noise-ratio X-ray spectrum of a quasar at a redshift higher than 0.4. These absorbers show no variability over a two-year timescale and have no associated cold absorption, making the assumption that they originate from the quasar's intrinsic outflow or the host galaxy's interstellar medium implausible. The O VII systems lie in regions characterized by large (four times larger than average ¹¹ ) galaxy overdensities and their number (down to the sensitivity threshold of our data) agrees well with numerical simulation predictions for the long-sought warm-hot intergalactic medium. We conclude that the missing baryons have been found.

Testing the Universality of Free Fall with the Triple System J0337+1715

The millisecond pulsar PSR J0337+1715 is in a mildly relativistic hierarchical triple system with two white dwarfs. This offers the possibility of testing the universality of free fall: does the neutron star fall with the same acceleration as the inner white dwarf in the gravity of the outer white dwarf? We have carried out an intensive pulsar timing campaign, yielding some 27000 pulse time-of-arrival (TOA) measurements with a median uncertainty of 1.2 μs. Here we describe our analysis procedure and timing model.

Gravitational wave research using pulsar timing arrays

A pulsar timing array (PTA) refers to a program of regular, high-precision timing observations of a widely distributed array of millisecond pulsars. Here we review the status of the three primary PTA projects and the joint International Pulsar Timing Array project. We discuss current results related to ultra-low-frequency gravitational wave searches and highlight opportunities for the near future.

Constraints on the Dynamical Environments of Supermassive Black-Hole Binaries Using Pulsar-Timing Arrays

We introduce a technique for gravitational-wave analysis, where Gaussian process regression is used to emulate the strain spectrum of a stochastic background by training on population-synthesis simulations. This leads to direct Bayesian inference on astrophysical parameters. For pulsar timing arrays specifically, we interpolate over the parameter space of supermassive black-hole binary environments, including three-body stellar scattering, and evolving orbital eccentricity. We illustrate our approach on mock data, and assess the prospects for inference with data similar to the NANOGrav 9-yr data release.

Detection methods for stochastic gravitational-wave backgrounds: a unified treatment

We review detection methods that are currently in use or have been proposed to search for a stochastic background of gravitational radiation. We consider both Bayesian and frequentist searches using ground-based and space-based laser interferometers, spacecraft Doppler tracking, and pulsar timing arrays; and we allow for anisotropy, non-Gaussianity, and non-standard polarization states. Our focus is on relevant data analysis issues, and not on the particular astrophysical or early Universe sources that might give rise to such backgrounds. We provide a unified treatment of these searches at the level of detector response functions, detection sensitivity curves, and, more generally, at the level of the likelihood function, since the choice of signal and noise models and prior probability distributions are actually what define the search. Pedagogical examples are given whenever possible to compare and contrast different approaches. We have tried to make the article as self-contained and comprehensive as possible, targeting graduate students and new researchers looking to enter this field.

A repeating fast radio burst

Fast radio bursts are millisecond-duration astronomical radio pulses of unknown physical origin that appear to come from extragalactic distances. Previous follow-up observations have failed to find additional bursts at the same dispersion measure (that is, the integrated column density of free electrons between source and telescope) and sky position as the original detections. The apparent non-repeating nature of these bursts has led to the suggestion that they originate in cataclysmic events. Here we report observations of ten additional bursts from the direction of the fast radio burst FRB 121102. These bursts have dispersion measures and sky positions consistent with the original burst. This unambiguously identifies FRB 121102 as repeating and demonstrates that its source survives the energetic events that cause the bursts. Additionally, the bursts from FRB 121102 show a wide range of spectral shapes that appear to be predominantly intrinsic to the source and which vary on timescales of minutes or less. Although there may be multiple physical origins for the population of fast radio bursts, these repeat bursts with high dispersion measure and variable spectra specifically seen from the direction of FRB 121102 support an origin in a young, highly magnetized, extragalactic neutron star.