Searching for Domain Wall Network by Parkes Pulsar Timing Array

We search for stochastic gravitational wave background generated by domain wall networks in the Data Release-2 of Parkes Pulsar Timing Array and find that the observed strong common power-law process can be explained by domain wall networks for the wall tension $\sigma_{\rm{DW}}\sim (29-414~\rm{TeV})^3$ and the wall-decay temperature $T_d\sim 26-363~\rm{MeV}$. Interestingly, the same parameter region can largely alleviate the Hubble tension, if the free particles generated from domain wall networks further decay into dark radiation. In addition, the preferred parameter space corresponds to the axion mass range $m_a \sim 10^{-13}-10^{-8}\ {\rm eV}$ for QCD axion. On the other hand, assuming that the common power-law process is not due to domain wall networks, we can put stringent constraints on the wall tension and decay temperature around the energy scale of QCD phase transition.