Testing clockwork axion with gravitational waves

We investigate the gravitational waves (GWs) produced from the Peccei-Quinn (PQ) phase transition associated with the clockwork axion. The PQ phase transition can be first-order when the dimension-6 operator is included into the scalar potential. The GWs from the PQ phase transition at scale in the range of $10^3-10^6$ GeV are detectable for the BBO and ALIA interferometers. The LISA and Taiji interferometers can probe the GWs from the PQ scale $f\lesssim 10^4$ GeV, while the GW signals from the scale $f\gtrsim 10^5$ GeV can be detected by the ground-based GW observatories ET and CE. We find that the parameter space $\kappa_m\sim 0.06-0.001$, $\kappa_l\sim 0.04-0.001$, and $\varepsilon\sim 0.1-0.01$ at the scale $f=10^5$ GeV and most of the parameter regions at the scale $f=10^6$ GeV have been excluded by the LIGO O2 run. The LIGO O3 and design phases can further probe the remaining parameter space. We show that the GWs from the annihilation of domain walls with a PQ scale $f\simeq 2\times 10^5$ GeV can induce the stochastic signals indicated by the 12.5-year observation of NANOGrav. The LIGO O3 run has the opportunity of detecting the GW signals from the first-order PQ phase transition around this scale.