Axionic domain walls at Pulsar Timing Arrays: QCD bias and particle friction

The recent results from the Pulsar Timing Array (PTA) collaborations show the first evidence for the detection of a stochastic background of gravitational waves at the nHz frequencies. This discovery has profound implications for the physics of both the late and the early Universe. In fact, together with the interpretation in terms of supermassive black hole binaries, many sources in the early Universe can provide viable explanations as well. In this paper, we study the gravitational wave background sourced by a network of axion-like-particle (ALP) domain walls at temperatures around the QCD crossover, where the QCD-induced potential provides the necessary bias to annihilate the network. Remarkably, this implies a peak amplitude at frequencies around the sensitivity range of PTAs. We extend previous analysis by taking into account the unavoidable friction on the network stemming from the topological coupling of the ALP to QCD in terms of gluon and pion reflection off the domain walls at high and low temperatures, respectively. We identify the regions of parameter space where the network annihilates in the scaling regime ensuring compatibility with the PTA results, as well as those where friction can be important and a more detailed study around the QCD crossover is required.