Footprints of the QCD crossover on cosmological gravitational waves at pulsar timing arrays

Pulsar timing arrays (PTAs) have reported evidence for a stochastic gravitational wave (GW) background at nanohertz frequencies, possibly originating in the early Universe. We show that the spectral shape of the low-frequency (causality) tail of GW signals sourced at temperatures around T≳1 GeV is distinctively affected by confinement of strong interactions (QCD), due to the corresponding sharp decrease in the number of relativistic species, and significantly deviates from ∼f3 commonly adopted in the literature. Bayesian analyses in the NANOGrav 15 years and the previous international PTA datasets reveal a significant improvement in the fit with respect to cubic power-law spectra, previously employed for the causality tail. While no conclusion on the nature of the signal can be drawn at the moment, our results show that the inclusion of standard model effects on cosmological GWs can have a decisive impact on model selection.