Novel tests of gravity using nano-Hertz stochastic gravitational-wave background signals

Gravity theories that modify General Relativity in the slow-motion regime can introduce nonperturbative corrections to the stochastic gravitational-wave background~(SGWB) from supermassive black-hole binaries in the nano-Hertz band, while remaining perturbative in the highly-relativistic regime and satisfying current post-Newtonian~(PN) constraints. We present a model-agnostic formalism to map such theories into a modified tilt for the SGWB spectrum, showing that negative PN corrections (in particular -2PN) can alleviate the tension in the recent pulsar-timing-array data if the detected SGWB is interpreted as arising from supermassive binaries. Despite being preliminary, current data have already strong constraining power, for example they set a novel (conservative) upper bound on theories with time-varying Newton's constant at least at the level of $\dot{G}/G \lesssim 10^{-5} \text{yr}^{-1}$ for redshift $z=[0.1\div1]$. We also show that NANOGrav data are best fitted by a broken power-law interpolating between a dominant -2PN or -3PN modification at low frequency, and the standard general-relativity scaling at high frequency. Nonetheless, a modified gravity explanation should be confronted with binary eccentricity, environmental effects, nonastrophysical origins of the signal, and scrutinized against statistical uncertainties. These novel tests of gravity will soon become more stringent when combining all pulsar-timing-array facilities and when collecting more data.