Gravitational waves from domain walls in pulsar timing array datasets

We present a model-independent search for a gravitational wave background from cosmic domain walls (DWs) in the NANOGrav 12.5 years dataset and International PTA Data Release 2. DWs that annihilate at temperatures " 20-50 MeV with tensions " (40-100 TeV)3 provide as good a fit to both datasets as the astrophysical background from supermassive black hole mergers. DWs may decay into the Standard Model (SM) or a dark sector. In the latter case we predict an abundance ANeff of dark radiation well within the reach of upcoming CMB surveys. Complementary signatures at colliders and laboratories can arise if couplings to the SM are present. As an example, we discuss heavy axion scenarios, where DW annihilation may interestingly be induced by QCD confinement.