An even lighter QCD axion

We explore whether the axion which solves the strong CP problem can naturally be much lighter than the canonical QCD axion. The ZN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {Z}_{\mathcal{N}} $$\end{document} symmetry proposed by Hook, with N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} mirror and degenerate worlds coexisting in Nature and linked by the axion field, is considered in terms of generic effective axion couplings. We show that the total potential is safely approximated by a single cosine in the large N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} limit, and we determine the analytical formula for the exponentially suppressed axion mass. The resulting universal enhancement of all axion interactions relative to those of the canonical QCD axion has a strong impact on the prospects of axion-like particle experiments such as ALPS II, IAXO and many others: experiments searching for generic axion-like particles have in fact discovery potential to solve the strong CP problem. The finite density axion potential is also analyzed and we show that the ZN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {Z}_{\mathcal{N}} $$\end{document} asymmetric background of high-density stellar environments sets already significant model-independent constraints: 3 ≤ N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} ≲ 47 for an axion scale fa ≲ 2.4 × 1015 GeV, with tantalizing discovery prospects for any value of fa and down to N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} ∼ 9 with future neutron star and gravitational wave data, down to the ultra-light mass region. In addition, two specific ultraviolet ZN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {Z}_{\mathcal{N}} $$\end{document} completions are developed: a composite axion one and a KSVZ-like model with improved Peccei-Quinn quality.