Abstract:The NANOGrav, Parkes and European Pulsar Timing Array (PTA) experiments have collected strong evidence for a stochastic gravitational wave background in the nHz-frequency band. In this work we perform a detailed statistical analysis of the signal in order to elucidate its physical origin. Specifically, we test the standard explanation in terms of supermassive black hole mergers against the prominent alternative explanation in terms of a first-order phase transition. By means of a frequentist hypothesis test we find that the observed gravitational wave spectrum prefers a first-order phase transition at $2-3\sigma$ significance compared to black hole mergers (depending on the underlying black hole model). This mild preference is linked to the relatively large amplitude of the observed gravitational wave signal (above the typical expectation of black hole models) and to its spectral shape (which slightly favors the phase-transition spectrum over the predominantly single power-law spectrum predicted in black hole models). The best fit to the combined PTA data set is obtained for a phase transition which dominantly produces the gravitational wave signal by bubble collisions (rather than by sound waves). The best-fit (energy-density) spectrum features, within the frequency band of the PTA experiments, a crossover from a steeply rising power law (causality tail) to a softly rising power law; the peak frequency then falls slightly above the PTA-measured range. Such a spectrum can be obtained for a strong first-order phase transition in the thick-wall regime of vacuum tunneling which reheats the Universe to a temperature of $T_*\sim \text{GeV}$. A dark sector phase transition at the GeV-scale provides a comparably good fit.

Induced Gravitational Wave interpretation of PTA data: a complete study for general equation of state

Induced gravitational wave interpretation of PTA data: a complete study for general equation of state

Probing the equation of state of the early Universe with pulsar timing arrays

Unraveling the early universe's equation of state and primordial black hole production with PTA, BBN, and CMB observations

High Frequency Gravitational Waves from Pulsar Timing Arrays

Constraints On Scalar-Induced Gravitational Waves Up To Third Order From Joint Analysis of BBN, CMB, And PTA Data

Induced Gravitational Waves with Kination Era for Recent Pulsar Timing Array Signals

Simultaneously probing the sound speed and equation of state of the early Universe with pulsar timing arrays

Untangling PBH overproduction in $w$-SIGWs generated by Pulsar Timing Arrays for MST-EFT of single field inflation

Scrutinizing the primordial black hole interpretation of PTA gravitational waves and JWST early galaxies

From formation to evaporation: Induced gravitational wave probes of the primordial black hole reheating scenario

The Gravitational Wave Background from Massive Black Hole Binaries in Illustris: spectral features and time to detection with pulsar timing arrays

Scrutinizing the Primordial Black Holes Interpretation of PTA Gravitational Waves and JWST Early Galaxies

Explaining PTA Data with Inflationary GWs in a PBH-Dominated Universe

Scalar-induced gravitational wave interpretation of PTA data: the role of scalar fluctuation propagation speed

Origin of the Stochastic Gravitational Wave Background: First-Order Phase Transition vs. Black Hole Mergers

Search for scalar induced gravitational waves in the International Pulsar Timing Array Data Release 2 and NANOgrav 12.5 years datasets

Obviating PBH overproduction for SIGWs generated by Pulsar Timing Arrays in loop corrected EFT of bounce