Axion Quality Straight from the GUT

Composite axion scenarios offer a robust field theoretic justification for the existence of a Peccei-Quinn symmetry of high quality. We present a new class of realizations that are naturally embedded in Grand-Unified Theories, retain asymptotic freedom for all gauge groups, and protect the axion symmetry up to operators of dimension 12. Our setup leads to a number of distinctive signatures at low energies. First, additional composite scalars are predicted; some of these are viable dark matter candidates for values of the axion decay constant that are too low for the QCD axion abundance to be relevant. Second, an approximate unification of the Standard Model gauge couplings takes place at the axion scale, while leaving the actual quark-lepton unification at much higher energies as usual. This suggests the existence of GUT relics with Standard Model gauge quantum numbers at potentially accessible scales.


Introduction
The non-observation of the electric dipole moment of the neutron sets an upper bound on the QCD topological angle: The smallness of |θ| can be explained by postulating non-generic properties of the UV completion of the Standard Model (SM). The identification of such properties is a challenging theoretical and experimental task, called Strong CP Problem. Among the various solutions of the Strong CP Problem are UV completions that invoke a global U (1) PQ anomalous under QCD.
Within the new physics sector the symmetry U (1) PQ may be linearly realized, explicitly broken, or spontaneously broken. When the U (1) PQ is linearly realized,θ can be rotated away becoming unphysical. By anomaly matching, once the massive new physics states are integrated out the low energy spectrum necessarily includes (composite or fundamental) massless colored fermions. Unfortunately, experiments firmly preclude the possibility of exotic massless quarks, and also the massless up-quark hypothesis [1] appears to be rather strongly rejected (see, e.g., [2] and more recently [3]). If instead the U (1) PQ is explicitly broken by the new physics, the full theory has a CP-odd angleθ full =θ −θ constructed as a linear combination of the QCD topological angle and some other anomalous parameterθ of the new physics. Postulating a softly-broken Z 2 mirror symmetry such thatθ =θ makes the overall phaseθ full vanish. [4] Because this is the quantity we actually measure below the mirror sector scale, the strong CP problem is solved. To be compatible with collider searches, explicit realizations of this idea must ensure the scale of soft Z 2 breaking lies within a certain range.
In the most popular U (1) PQ -based approach to the Strong CP Problem, the anomalous Peccei-Quinn symmetry is spontaneously broken and an approximate Nambu-Goldston boson, the QCD axion a, is predicted. [5] An anomalous rotation of the colored fermions does not remove the topological angle; rather, it shifts the axion in such a way as to keep the dynamical fieldā/f a =θ + a/f a invariant. Remarkably, the vacuum energy V QCD (ā) of QCD gets minimized at ā = 0. This means that the value of the observable topological angle gets automatically relaxed to zero,θ obs ≡ ā = 0, provided the axion potential is entirely dominated by non-perturbative QCD effects.
The QCD axion has been sought after for many years. No evidence has ever been found and the experimental searches are interpreted as constraints on the parameter space. Because all axion couplings are controlled by the scale f a of spontaneous U (1) PQ breaking, data translate into lower bounds on this quantity. For the invisible axion models considered here, the most stringent constraints read f a 10 8 GeV. Such large scales make the axion potential extremely sensitivity to tiny effects from unknown physics at the shortest imaginable distances. [6] Being an anomalous symmetry, U (1) PQ cannot be exact in the far UV. At best, we may expect it to be broken solely by quantum gravity effects parametrized at low energies via Planck-scale suppressed interactions. These induce corrections ∆V to the axion potential that are generically not aligned with V QCD (ā), the origin of field space is destabilized and a non-vanishing observableθ obs ≡ ā /f a is re-introduced. The full potential is now the sum of a term due entirely to QCD and one from tiny effects at the Planck scale: At leading order in ∆V we findθ The question of why |θ obs | is so small becomes, in axion models, the question of why the U (1) PQ is so accurate despite it being only an approximate symmetry of the UV dynamics. This "U (1) PQ quality problem" is so significant because the QCD potential is controlled by physics at the GeV scale, V QCD ∼ m 2 π f 2 π , whereas ∆V is a function of the large f a and the Planck scale f Pl = 2.4 × 10 18 GeV. Using naive dimensional analysis we expect ∆V ∼ c 16π 2 f 4 a (f a /f Pl ) d−4 , with c some coefficient and d the dimension of the associated U (1) PQviolating operator. For a generic, maximally strong UV dynamics c ∼ 1. As long as f a f Pl one can hope to suppress the destabilizing effect by making sure the operator is sufficiently irrelevant. Plugging these estimates in (3), the experimental bound (1) translates into d 9 for f a ∼ 10 8 GeV. A scale of order, say, f a ∼ 10 11 GeV would instead require d 12. As f a approaches the Planck scale the only way to prevent |θ obs | from exceeding the experimental limit is to have |c| 10 −10 . Unfortunately, to a low energy observer such a condition is no different than the experimental input (1) we started with. The necessary non-generic UV conditions to ensure |c| 1 may be somehow realized in the correct theory of quantum gravity, for example in String Theory. However, in this paper we take the point of view of the low energy observer and look for explanations within the realm of quantum field theory.
We see that robust realizations of the QCD axion within a field theory approach require two fundamental ingredients: an axion decay constant parametrically smaller than the Planck scale and an accurate U (1) PQ . Composite axion models can feature both ingredients. In such a framework the axion arises as a Nambu-Goldstone boson (NGB) of some exotic strong dynamics that undergoes chiral symmetry breaking at a scale ∼ f a . [7] As long as the axion constituents are chiral under the new gauge interactions, a number of accidental chiral symmetries, among which U (1) PQ , may be very accurate. Furthermore, a value of f a that is Table 1: Axion constituents, and their accidental U (1) symmetries, in our toy-model. simultaneously parametrically larger than the weak scale and lower than the Planck scale is completely natural in composite axion models, where f a gets generated via dimensional transmutation. In non-supersymmetric models with fundamental scalars, on the other hand, the value of f a results from a fine-tuning that is much more severe than just setting |θ| < 10 −10 by hand in a given UV completion of the SM.
The first composite axion model with chiral constituents was proposed in [8]. In those types of models the fermion constituents are charged under three non-abelian gauge symmetries and to have a sufficiently accurate U (1) PQ some of the gauge groups must have rank bigger than a certain lower bound. As a consequence there is a structural tendency to developing low energy Landau poles. [9] One can evade this problem by replacing one or more of these nonabelian factors with abelian groups [10] (one should also worry about cubic anomalies in this case). The price to pay is however a loss of asymptotic freedom. Alternatively, one can construct asymptotically-free models using a "moose" gauge structure [11] that resembles axion scenarios in a local extra dimension [12].
The purpose of the present paper is proposing an alternative class of asymptotically-free composite axion models with a natural embedding in some Grand-Unified Theory (GUT) and a U (1) PQ protected up to operators of dimension d = 12.
2 Warming up with SU (N c ) × SU (5) 2 We begin our discussion with a very simple model. Unfortunately its protection of the U (1) PQ is at most barely sufficient. Yet, this toy model allows us to illustrate some of the key features that characterize the more efficient scenarios discussed in Section 3.
Our toy-model is based on an SU (N c ) × SU (5) 1 × SU (5) 2 gauge symmetry. The axion constituents are shown in Table 1. They are all Weyl fermions, chiral under SU (5) 2 . The SU (N c ) is the dynamics that becomes strong and generates a composite axion at scales of order f . The SU (5) 1 is to be identified with the GUT, with the SM fermions belonging to three families in the 10 ⊕ 5 ∈ SU (5) 1 and being singlets of SU (N c ) × SU (5) 2 . The SU (5) 1 GUT is not instrumental to the realization of our composite axion scenarios. Yet, being one of our motivations we find it useful to entertain the possibility of grand-unification. In this case scalar fields at the GUT scale should be added to obtain a realistic model. The axion physics in our models is not going to be affected by them provided they do not carry SU (N c )×SU (5) 2 charges, which is what we assume. No other fundamental scalar is introduced. 1 At scales of order f GUT the GUT is broken into the SM as usual: At f f GUT , when the strong dynamics of the axion becomes important, it is appropriate to decompose ψ, ψ in irreducible G SM,1 representations: ψ = ψ q ⊕ ψ u ⊕ ψ e , ψ = ψ d ⊕ ψ , where ψ q,u,e,d, are in the fundamental of SU (N c ), singlets under SU (5) 2 and transform under G SM,1 as the SM quarks and leptons q, u, e, d, .
Let us next discuss under which conditions our gauge interactions remain asymptotically free. The requirement of asymptotic freedom for SU (5) 2 is satisfied for any N c ≤ 13. This is a very mild constraint. A much more significant upper bound on N c is found demanding that SU (5) 1 remains UV-free. This bound is model-dependent because it depends on the additional fields necessary to explain (4) and the SM Yukawa couplings. For definiteness we consider a minimal (non-supersymmetric) scenario in which the scalar responsible for breaking SU (5) 1 into the SM is in the representation 24. The observed pattern of Yukawa couplings for the SM fermions may be reproduced with a Higgs in the 5. Rather than introducing additional fundamental fields, which would quickly take our model away from asymptotic freedom, the necessary quark-lepton splitting is achieved via effective dimension-5 operators with two SM fermions and two powers of the scalars φ H ∼ 24 ⊕ 5. These interactions are large enough provided With this minimal GUT-scale field content the upper bound reads N c ≤ 9. The key assumption that SU (N c ) confines and breaks its chiral symmetries provides instead a lower bound on N c . Ignoring the weak gauging of SU (5) 2 , the SU (N c ) dynamics behaves as if it was a QCD-like theory with N c colors and N f = 15 massless flavors. According to a number of independent arguments (for an incomplete list see, e.g., [13][14] [15][16] [17]), confinement and dynamical chiral symmetry breaking occurs when N f /N c 3 − 4. We will work in the regime 5 ≤ N c ≤ 9. Under our assumptions, at f ∼ 10 8 − 10 9 GeV the SU (N c ) becomes strong and non-trivial vacuum condensates form: where we borrowed the factors of 4π from QCD and the powers of N c from large N c scaling. 2 The particle content of Table 1 features a number of accidental symmetries. We spare the reader of the details, which will instead be presented for the models of Section 3. The bottom line is that (6) break the gauge symmetry G SM,1 × SU (5) 2 into the diagonal G SM , which is to be identified with the SM gauge symmetry we experimentally observe. In the process, many Nambu-Goldstone bosons (NGBs) are generated. Some are eaten via the Higgs mechanism, some acquire large masses squared ∼ g 2 SM f 2 /N c , and finally four SM-neutral NGBs remain massless at the renormalizable level. Among these, three are associated to chiral symmetries that appear below the GUT scale but are in fact explicitly broken by the weak gauging of 2 For QCD in the chiral limit we have ψψ c (240 MeV) 3 at µ ∼ 1 GeV (no summation over flavor). This agrees with the naive dimensional analysis estimate ψψ c ∼ 4πf 3 π up to a number close to 1.
SU (5) 2 , which explains why these are not shown in Table 1. As a consequence, these three NGBs acquire masses from dimension-6 interactions below the GUT scale (we will study them in Section 4.2). The actual QCD axion a is associated to the breaking of the anomalous U (1) PQ of Table 1.
The lowest dimensional operators that violate U (1) PQ are formed by the product of ψ †σµ ψ and appropriate SM fermion currents, and have dimension d = 6. Yet, such interactions are innocuous. Indeed, to generate an axion potential ∆V in our QCD-like models one must break In the presence of explicit breaking of SU (N f ) L only, as it is the case for the dimension-6 interactions just mentioned, one could remove completely from the Lagrangian the non-derivative couplings of the axion performing a global SU (N f ) R rotation. This demonstrates that the above dimension-6 interactions are not harmful unless they are accompanied by operators that also violate the U (1) PQ ⊂ SU (N f ) R carried by ψ c , ψ c . Fortunately, the couplings of the latter fields are more severely constrained by invariance under SU (5) 2 and SU (N c ), as we now argue.
The leading SU (5) The former requires three powers of ψ, ψ to make an SU (N c )-singlet, leading us to operators of at least dimension d = 9. The other option, i.e. ψ c ψ c ( ψ c ) † , may be combined with ψψ( ψ) † , which again gives d = 9 operators. Hence, we conclude that ∆V is first induced by operators of the type Strictly speaking, one may also build an SU (N c ) singlet from ψ c ψ c ( ψ c ) † by adding either ψ or ψ (or even their conjugates if we are willing to pay the additional price of a gauge field strength), but the resulting operator would not be invariant under G SM,1 unless at least an additional SU (N c )-singlet, say a Higgs and/or an even number of SM fermions, are added. One can verify that the effect of these operators on the axion potential is negligible compared to the one coming directly from (7). Having established that the dominant contribution to ∆V arises from the set of operators in eq. (7), we can finally estimate the expected size of the observed topological angle (3). For later convenience we analyze the impact of a general operator of dimension d of the type (ΨΨ c ) d/3 , with Ψ, Ψ c the axion constituents. We expect the operator to be suppressed by the ultimate cutoff, say the Planck scale f Pl 2.4 × 10 18 . Naive dimensional analysis is used to estimate the Wilson coefficient. Denoting by g Pl the typical coupling at the cutoff, e.g. the string coupling, the operator appears in the effective Lagrangian as where c Pl is of order unity if we conservatively assume the operator is generated at tree-level by exotic states of mass ∼ g Pl f Pl . For any d > 6 the Wilson coefficient gets smaller as the coupling increases (we have fixed f Pl ). Large-N counting suggests that the maximal allowed value satisfies g 2 Pl N dof ∼ 16π 2 , with N dof 1 the number of degrees of freedom. In our estimates we will saturate the bound with N dof ∼ N f N c . Fortunately, our numerical results are modestly affected by the unknown parameter g Pl .
One should then run the Wilson coefficient down to energies of order the confinement scale, i.e. c Pl → c IR . We may estimate this effect assuming (as in the large N c approximation) the operator anomalous dimension factorizes as d/3 times the anomalous dimension of the quark bilinear ΨΨ c . Limiting our analysis to a domain g 2 IR N c /(4π) < 0.5 − 1 in which a perturbative calculation is expected to be reliable, we observe an enhancement of order c IR /c Pl ∼ 10 − 10 2 . The contribution to the axion potential due to (8) can thus be written as where q a is the U (1) PQ charge of (8) and is a numerical factor that depends on the number of condensates (6) that appear in the operator (8). For example, the second operator in (7) contributes with n O = 10 2 × 5 1 and q a = 4/ √ 15. For generic complex coefficients c IR , uncorrelated with the bareθ, ∆V implies a contribution to the QCD topological angle (see eq. (3)) .
Here we employed the leading order expression for the axion mass m 2 where A G a is the QCD anomaly coefficient determining the axion coupling to gluons (see Section 3.2 for details). Note that the powers of , as well as the RG effects, go in the direction of enhancing |θ obs |. Other than that, the estimate (10) agrees very well with the one given in Section 1.
Employing typical values of n O and q a (such as those that have been mentioned above (10)), the constraint (1) reads f a 10 8 GeV, which falls too short to convince us this is a viable model. One can certainly hope that the actual value of the unknown coefficient c Pl may be suppressed in a realistic theory of gravity. But this is not what we do here. Our point of view is the one of an agnostic low energy observer who would like to find an answer to a puzzling result, i.e. (1), without having to resort to unjustified assumptions about the UV. In other words, we are interested in finding scenarios that work with coefficients c Pl of order unity. Such realizations are presented in Section 3.
Before concluding this discussion we would like to stress that our toy model differs significantly from the one of Ref. [18], where the axion constituents are also in a chiral 10 ⊕ 5 representation of some SU (5). In that paper the gauge group SU (5) is identified with the strong axion dynamics; the constituents are therefore chiral under the confining force but vector-like under the SM. Here the situation is completely reversed. Our axion constituents carry vector-like representations under the confining SU (N c ) dynamics and are instead chiral under the weakly-gauged, grand-unified SU (5) 2 . The dimensionality of the leading operators that violate U (1) PQ is inherited from the chiral 10 ⊕ 5 ∈ SU (5), and is thus the same here as well as [18], though.
The previous model does not quite reach the desired level of protection of the axion U (1) PQ . The ultimate reason is that it is sufficient to have three powers of 10 ⊕ 5 ∈ SU (5) to build a complete gauge singlet under the chiral SU (5) 2 . Our next models prevent this from happening by embedding the axion constituents in the spinorial representation of SO (10). Apart from a more significant suppression of ∆V , in many respects the new models are qualitatively similar to our earlier toy scenario.

Model Setup and Symmetry Breaking
Consider then a model based on the gauge symmetry SU (N c ) × SO(10) 4 . The representations of the axion constituents are shown in Table 2, and are all chiral under the weakly-gauged SO(10) 4 . SO(10) 1 is the GUT and the SM fermions (plus the right-handed neutrinos) are as usual embedded in three generations of 16 i . They are singlets of SU (N c ) and SO(10) 2,3,4 . As in Section 2, other fields charged under SO(10) 1 , but not under SU (N c ) and SO(10) 2,3,4 , can and must be included to obtain a realistic GUT. Specifically, the minimal (non-supersymmetric) set of Higgses belong to the representation φ H ∼ 10 ⊕ 16 ⊕ 45 ∈ SO(10) 1 . The observed pattern of Yukawa couplings for the SM fermions may be reproduced, as long as (5) As in the toy model, at f f GUT the SU (N c ) becomes strong and a non-trivial vacuum for the fermion bilinears form: Switching off the gauge SO(10) 4 couplings the particle content of Table 2 Table 2. Below f GUT we can decompose ψ in irreducible representations of G SM,1 : with ψ N a singlet of G SM,1 and SO(10) 2,3,4 . The global SO(10) 1 is thus further broken into G SM,1 × U (1) 4 . Without loss of generality the abelian factor U (1) 4 can be chosen to act only on the components ψ q,u,e,d, ,N . It is very simple to identify it. Of the five independent nonanomalous rotations that these fields can be subject to, one is hyper-charge and is already included in G SM,1 ; the remaining four are shown explicitly in Table 3. Table 3: Accidental U (1) global symmetries of the fields of Table 2 arising below the GUT scale.
The symmetry breaking pattern induced by the condensates (12) turns out to be where for brevity we introduced the notation for the broken axial charges. The condensate ψψ c breaks G SM,1 × SO (10) Table 3. The vectorial combination is left unbroken by the condensates (12). This is what we denoted by U (1) S,V in (13). Without the gauged SO (10)  There are then four complete singlets that remain massless at the renormalizable level. As we anticipated for the model in Section 2, these arise from the breaking of (14). One is obviously the axion in the last line of (15), and is associated to the breaking of U (1) PQ . The other three are contained in the SM-singlet components of (210, 1). They arise from the breaking of the accidental U (1) A1 × U (1) A2 × U (1) A3 appearing below the GUT scale. A completely analogous thing happens in the toy model of the previous section, as the reader can verify.
The four singlet NGBs will be denoted by a, a 1 , a 2 , a 3 , respectively. They can be described by a 32-dimensional special unitary matrix Σ = e 2iaαTα/f transforming as Σ → LΣR * under SU (32) L × SU (32) R , where the generators, normalized such that Tr[T α T β ] = δ αβ /2, are diagonal and can be readily seen from Table 2 and 3. In particular T PQ = 1 8 diag(1 16×16 , −1 16×16 ) and similarly for T A1,A2,A3 . The factor of 2 in Σ is chosen to conform to the standard definition 0|J µ PQ (0)|a(p) = if p µ . Some of these pseudo-scalars have anomalous couplings to the SM gauge fields proportional to the anomaly coefficients 3 where T A gauge are the SM generators, with A, B indices in the adjoint, and the trace runs over all the flavors of the axion constituents. For example, those of a read . Similar considerations apply to a 1,2,3 . Despite a unique linear combination of a, a 1,2,3 can be chosen to couple to the SM gluons, the actual QCD axion is to be identified with a. Indeed, as anticipated before and shown explicitly in Section 4.2, the scalars a 1,2,3 acquire masses above the TeV from dimension-6 interactions. Hence they cannot be relevant to the Strong CP Problem. The only potential impact they can have is to shiftθ, an effect which is completely irrelevant to our purposes. On the other hand the QCD axion a remains massless up to a much higher degree, see Section 3.2. Our focus for the rest of the section will therefore be on a. Conventionally the scale f a is defined as Before moving on let us note we can verify a posteriori the vectorial gauge group G SM × SO(10) 3+4 remains unbroken. All the physical NGBs charged under it have large positive masses squared and, at the renormalizable level, vanishing vacuum expectation values. [19][20] The (misaligning) effect induced by unavoidable higher-dimensional operators is suppressed by at least powers of f 2 /f 2 GUT g 2 unbroken /16π 2 and hence is much smaller than the (aligning) contributions due to the gauge couplings.

UV Contributions to the Axion Potential
To relaxθ to zero, and thus solve the Strong CP Problem, we have to make sure that the NGB associated to the breaking of U (1) PQ , namely the QCD axion a, has a potential dominated by non-perturbative QCD effects. The UV correction ∆V of eq. (2) originates from local operators that explicitly break the chiral symmetry U (1) PQ . In these models such effects are severely constrained by the local SO(10)'s.
The most minimal gauge singlets built out of a 16 ∈ SO(10) are either 16⊗16⊗16⊗16 or 16 ⊗ 16. With the particle content of Table 2 each fermion is charged under a single SO(10), so the combination 16 ⊗ 16 is automatically invariant under U (1) PQ . We are thus left with the former structure, which by construction cannot be SU (N c ) invariant, see (11). To form a complete gauge and Lorentz singlet, four powers of ψ, ψ should be accompanied by at least four powers of ψ c , ψ c , and vice-versa. It follows that the most relevant operators that violate U (1) PQ are of the type and thus belong to the class in eq. (8) with dimension d = 12. This is a qualitative improvement compared to the model of Section 2. We can next employ the general expression (10) to estimate the observable vacuum angle in this model. One can check that the operator of type (18) with the largest n O q a is ( ψ ψ c ) 4 , which contributes with n O = 16 4 = 65 536 (!) and |q a | = 1, see Table 3. Plugging these values in (10), together with A G a = N c from (17) and |c Pl | = 1 (the RG effect again leads to an enhancement of a factor of 10 2 ), the experimental constraint (1) becomes f a (3 − 4) × 10 9 GeV, with N c varying in the range (11). Had we not included the large factor of n O = 16 4 , nor the enhancement (A G a ) d−1 (3/N c ) d/6 ∼ N 9 c , nor the RG enhancement, we would have approached the optimistic f a 10 11 GeV mentioned in the introduction. Unfortunately such enhancements are there.
Finally we compare the above "stability" bound on f a to the lower bound from experimental searches. The low energy phenomenology of our models is in the KSVZ class. [21] The effective coupling of the axion to photons, after the mixing with π 0 is taken into account, is given by The coupling to neutrons and protons is g ann = −0.02(3)m n /f a and g app = −0.47(3)m p /f a (see, e.g., [23]). In these scenarios the astrophysics bounds are currently dominated by the supernova constraint on the proton coupling (see e.g. [24]) and read f a 1.3 × 10 8 GeV. We conclude that these scenarios have an allowed region of parameter space of one to two orders of magnitude in which the Strong CP Problem is robustly solved.

Alternative Realizations
One can obtain alternative models adding some or removing one of the weakly gauged SO(10) groups. Consider for instance removing one of the gauge factors in Table 2, say SO(10) 4 . The resulting scenario has a smaller SU (N c ) × SO(10) 3 gauge and the axion constituents are: 1, 1, 1).
The physics of the QCD axion is virtually unchanged compared to that of Section 3.2. In particular, also here the leading operators violating U (1) PQ have dimension 12.
The main difference is that now no additional gauge symmetry survives after chiral symmetry breaking besides the SM, i.e. SO(10) 3 is completely broken. Furthermore, as opposed to the model of Section 3 where only four neutral NGBs arise, here many more are present. These reside in the SM-singlets in (16,16) ⊕ (16,16) and of course in all of the (1, 210) of (15). Interestingly, among these are potential dark matter candidates, as discussed in the Section 4.3.2.

Signatures
Our models predict a number of generic signatures: an approximate gauge-coupling unification, heavy unstable ALPs, and dark matter candidates. The first two are essentially the same in the models of Section 2 and 3 as well as the one explicitly mentioned Section 3.3 and will be studied in Sections 4.1 and 4.2. Which stable relic can play the role of dark matter depends on the model, though. In the scenario of Section 3 natural candidates are the lightest glueball of the unbroken SO(10) 3+4 and/or heavy SM-neutral NGBs; in scenarios of the type of Section 3.3 dark matter could be in the form of a neutral GeV-scale NGB and/or a heavy hadron again. We will discuss these in Sections 4.3.1 and 4.3.2 respectively.

Gauge-Coupling Unification at f a ?
The axion constituents are chiral under the weakly-coupled gauge groups and the SM gauge symmetry observed at low energy, G SM , is always in the unbroken vectorial subgroup of the chiral symmetry. An approximate unification of the SM gauge couplings g G , g W , g Y ≡ 5/3g Y at the axion scale follows as a generic prediction. Let us discuss this feature in the context of the SU (N c ) × SO(10) 4 model.
Below the GUT scale the three G SM,1 couplings g G,1 , g W,1 , and g Y ,1 ≡ 5/3 g Y,1 start to depart from each other. When we reach µ ∼ µ gauge ≡ g GUT f / √ N c = g GUT √ N c f a we find that the observed SM gauge symmetry is in the vectorial subgroup of G SM,1 × SO(10) 2 . Eq. (13) induces a new threshold, which at tree-level (13) reads with g 2 the SO(10) 2 and all couplings are renormalized at µ ∼ µ gauge . Now, for generic values g 2 2 g 2 i,1 the threshold correction (20) is very mild, if relevant at all. In the complementary regime g 2 2 g 2 i,1 (for earlier speculations on this possibility see [25]) 4 the linear combination of vector fields that describes the low energy SM gauge bosons is mostly oriented along the SO(10) 2 vectors and the couplings observed at low energy thus have a comparable value This results in an approximate gauge coupling unification at scales µ ∼ µ gauge g GUT f GUT parametrically smaller than those at which true unification of quarks and leptons actually occurs. Gauge coupling unification is only approximate because an accurate one would require g 2 2 g 2 i,1 , which can only be realized if the couplings g 2 i,1 are non-perturbative at threshold.
Physics in the approximately unified regime g 2 2 g 2 i,1 is quite interesting. The gauge couplings almost unify at the axion threshold, and yet proton decay is not induced because the SO(10) 2 gauge bosons do not carry baryon nor lepton charges. One has to probe physics at the GUT scale in order to discover that the SM fermions fill complete SO(10) 1 GUT representations. We have thus separated gauge-coupling unification from proton decay. As an extreme limit, one can even push f GUT to the Planck scale, thus suppressing proton decay beyond any foreseeable experimental reach, while preserving an approximate gauge-coupling unification at relatively low scales. [25] Looking from bottom-up we see that approximate gauge coupling unification at the axion scale leads to a clear low-energy prediction: there must exist new light particles charged under the SM that make the gauge couplings g i unify at much lower scales than in usual GUTs. It is not hard to find relics of grand-unified multiplets that can serve our purpose. For example, four Dirac electroweak doublets with hypercharge ±1/2 plus two Dirac singlets with hypercharge ±1, all with masses of order the TeV, would be enough to reach SM gauge coupling unification at 10 11 GeV.

Unstable ALPs at the TeV
These models feature a unique axial symmetry at the renormalizable level, namely U (1) PQ . Yet, below the GUT scale an additional approximate one arises. For definiteness we focus on the model of Section 3, where the extra broken generators are those of U (1) A1 × U (1) A2 × U (1) A3 . Completely analogous considerations apply to the scenario of Sections 2 and 3.3. We are thus interested in the physics of the axion-like particles (ALPs) a ALP = {a 1 , a 2 , a 3 }, which we recall are the SM-singlets in (210, 1) of (15).
The key observation is that U (1) A1 × U (1) A2 × U (1) A3 is violated by the gauge couplings. Indeed, as opposed to T PQ , the corresponding generators do not commute with the generators of the heavy vectors, T L ∈ SO(10) 1 /G SM,1 . This latter fact explains why that symmetry did not appear above the GUT scale in Table 2. The tree exchange of the heavy vectors therefore leads to dimension-6 interactions of the type (ψ † T L ψ) 2 /f 2 GUT that break U (1) A1 × U (1) A2 × U (1) A3 completely. As we emphasized in Section 2, this is not sufficient to generate a potential for a ALP , though, because such interactions only break SU (N f ) L . The necessary breaking of SU (N f ) R comes from the gauging of T R ∈ SO(10) 3,4 . Combining the dimension-6 operators (ψ † T L ψ) 2 /f 2 GUT with a loop of SO(10) 3,4 vectors we get a contribution to the potential of a ALP , which technically can be parametrized by chiral invariants of the form . This translates into a mass matrix (including mass mixing terms) of the order Nothing like this can be written down for a, since U (1) PQ is not broken by the gauge interactions.
One can readily see from Table 3 that these extra NGBs have anomalous couplings ∝ A SM ALP to the SM gauge fields (B µ , W µ , G µ ), see above (16) for an explicit expression of A SM ALP , and are therefore unstable. Up to unknown mixing angles, their lifetimes are with n SM = 1, 3, 8 for SM = Y, W, G a multiplicity factor. Given typical values of the parameters the decay safely takes place at temperatures T 10 1−2 GeV. Hence, the heavy partners of the QCD axion disappear well before nucleosynthesis and lead to no obvious signatures.

Long-lived Relics
A generic problem of axion models is the production of topological defects via the Kibble-Zurek mechanism. To avoid conflicts with cosmology we assume the highest temperature reached by the Universe, T max , never exceeded the critical temperature of the SU (N c ) dynamics, expected to be of order Λ ≡ 4πf / √ N c = 4πf a √ N c by analogy with large N QCD: T max 4 × 10 9 GeV N c 12 Under this condition the SU (N c ) dynamics underwent its phase transition during or before inflation (we will see below that also H I 4πf a holds) and topological defects were subsequently inflated away.
Composite axion models also predict a whole tower of heavy hadrons. The lightest baryons have mass ∼ N c Λ and are cosmologically stable. Because of (23) their abundance is completely negligible throughout the thermal history of the Universe. Heavy mesons with masses of order ∼ Λ quickly decay into NGBs as in ordinary QCD. The physics of the NGBs is more modeldependent.

Heavy NGBs
In the scenario of Section 3 all NGBs besides the QCD axion and the three unstable singlets of Section 4.2 are heavy. They decay into SM fermions, gauge bosons, or into fermions and other NGBs. 5 The lightest NGB carrying SO(10) 3+4 charge, we call it Π s , is stable. We expect it to be the SM-neutral component of (16,16) ⊕ (16,16) in the third line of (15), say Π s ∼ ψ N ψ c , since the masses of all its charged partners receive additional positive contributions from the SM gauge loops. A robust way to ensure its abundance does not exceed the density of the observed dark matter, in conflict with the well established hot big-bang theory, is to postulate that the reheating temperature T RH ≤ T max is sufficiently smaller than its mass. This constraint roughly reads T RH < m Πs /30, where m 2 Πs ∼ (g 2 unbroken /16π 2 )C 2 Λ 2 . Hence we demand with g unbroken the coupling of SO(10) 3+4 and C 2 = 45/8 the quadratic Casimir of the spinorial of SO(10) 3+4 , the representation carried by Π s . With (24) satisfied, the model of Section 3 has a negligible abundance of SU (N c ) hadrons. Potential cosmological signatures may come from the SO(10) 3+4 gauge bosons, as discussed in Section 4.3.1. The situation is a bit different in the scenario of Section 3.3. In that case the NGBs in the (210, 1) of (15) again decay into SM fermions, and those in the third line of (15) again decay into the lightest of them, i.e. Π s . However, Π s and all the states in (1,210) are now massless at tree-level, and couple very weakly to the bath via higher-dimensional operators. Still, Π s may be abundantly produced via the decay of its heavy partners. To avoid this we impose the same condition as (24), with g unbroken replaced by the SM couplings. As a consequence of this assumption, also in the scenario of Section 3.3 the present-day population of heavy hadrons is negligible. The physics of the light, decoupled singlets Π s and (1, 210) will be analyzed in Section 4.3.2.

Axion Density
Eq. (23) implies the axion abundance is set by vacuum misalignment. Assuming for definiteness that the temperature-dependence of the axion mass is approximately given by m a ∝ T −δ for temperatures above the critical one, the fraction of dark matter in the form of cold QCD axions turns out to scale as where θ init,a is the initial misalignment angle. The result (25) depends uniquely on the initial condition and the (temperature-dependent) QCD potential, and is therefore governed by f a = f /N c rather than f , see (17). The main unknown is δ. Taking as a reference the result of a dilute instanton gas approximation, where δ ∼ 4, we see that for the relatively low values f a ∼ 10 8−10 GeV considered here the axion constitutes at most a small fraction of the observed dark matter. With f a fixed, the abundance decreases as δ decreases. During inflation the axion field undergoes quantum fluctuations of typical size of the Hubble scale H I in that epoch. The resulting isocurvature perturbations are constrained by cosmological observations. This latter constraint can be expressed as an upper bound on H I [29]: H I 10 7 GeV (f a /10 11 GeV) (1+δ)/(4+2δ) ρ DM /ρ a , which gets weaker as the QCD axion abundance decreases. Combining with (25) we obtain This bound applies to all pre-inflationary QCD axion scenarios. In our case it is a bit relaxed compared to what is often quoted because the axion is a small component of the dark matter.

Glueball Dark Matter
The model of Section 3 predicts an unbroken hidden non-abelian gauge group, i.e. SO(10) 3+4 . Below the scale Λ = 4πf a √ N c the associated vectors are coupled to ordinary gluons very weakly via higher dimensional operators of the form 6 and similarly for the electroweak bosons. The interaction is so irrelevant that they never equilibrate with the plasma. The hidden gluons form a gas of self-interacting particles living with its own entropy s (T ) at its own temperature T , which differs from the temperature T of radiation. At the end of inflation the hidden gluons are populated dominantly by the annihilation of SM gluons according to dn /dt + 3Hn = Γ, with H the Hubble scale and Γ the Boltzmann collision term. (By the assumption (24) the amount of vectors produced by the annihilation of the hadrons at the scale ∼ Λ is instead completely negligible.) After being produced, the hidden gas thermalizes right away with entropy s (T ) = (2π 4 /45ζ 3 )n (T ). The non-renormalizable nature of (27) implies Γ ∝ T 12 , so that the production is mostly active at the highest accessible temperatures T ∼ T RH . The entropy of the hidden gluons over the radiation entropy may be estimated to be of order where we find The quantity ξ then remains constant at all T T RH . At later times the SO(10) 3+4 force confines and the entropy carried by the relativistic hidden gluons gets re-distributed among non-relativistic glueballs. The heavy SO(10) 3+4 glueballs quickly decay into the lightest, Φ, of mass m Φ . A simple estimate shows that Φ has a lifetime longer than the age of the Universe for m Φ 5 × 10 5 GeV, c ∼ 1, N c in the range (11), and f a = 10 8−10 GeV.
The present-day abundance of Φ can be calculated as follows (see ref. [30]). After confinement the new gas of non-relativistic Φ's remains self-interacting for a while until the number-changing processes go out of equilibrium at T d m Φ . We can relate the temperatures of the hidden gas to that of ordinary radiation via s Φ (T d ) = s (T d ) = ξs(T d ). For the range of parameters we are interested in, T d is always many orders of magnitude above the keV and the energy density ρ Φ of the glueballs at decoupling is much smaller than that of radiation. For any T < T d the non-interacting glueball gas behaves as ordinary cold relics: ρ Φ /s = ξρ Φ /s stays constant and progressively becomes more and more relevant than radiation. The present-day glueball energy density is of order Using (29) and saturating the upper bound (24) the abundance of the stable glueballs reaches a value comparable to the observed dark matter. This conclusion is however extremely sensitive to T RH because of the high temperature-dependence in Γ. A dedicated study of the reheating process would be needed to clarify if Φ can be a realistic dark matter candidate in a regime in which the heavy NGBs have a negligible relic density.

NGB Dark Matter
In the scenario of Section 3.3 there is no hidden gauge symmetry but there is a whole host of light SM-neutral NGBs, as we have seen at the end of Section 4.3. In particular, Π s and those in the representation (1, 210) of (15) are SM singlets, massless at the renormalizable level, and have no anomalous couplings to the SM gauge fields. Some of these are exactly stable in the renormalizable Lagrangian, namely those that are protected by an "isospin" symmetry under which the SO(10) 3 gauge bosons transform in the adjoint.
contributions to the axion potential from operators of dimension 9, which may be too little to obtain a realistic axion model. An axion symmetry with higher quality is obtained exploiting the spinorial 16 of SO(10). In models of this type all SO(10) singlets with U (1) PQ charge contain at least four powers of the axion constituents. As a result the minimal Lagrangian operator that can contribute to the axion potential has dimension 12. We analyzed two models of this class, but many others can be built along these lines.
A rough estimate naively implies that Planck-scale suppressed operators of dimension 12 could be compatible with f a ∼ 10 11 GeV. However, composite axion models tend to have large color and flavor multiplicities -a must if we want to embed our scenarios in GUTs in which all gauge groups are asymptotically free -as well as renormalization group effects which go in the direction of enhancing the sensitivity to UV contributions to the potential. Employing naive dimensional analysis, the best tool the effective field theorist has at her or his disposal to estimate the Wilson coefficient of the Planck-scale suppressed operator, we find that our SO(10) models are left with a window of about one to two orders of magnitude 10 8 GeV f a 10 9−10 GeV in which the Strong CP Problem is robustly solved consistently with all experimental data. Our conservative estimate seems to us more justified than the naive one.
Interestingly, the very structure of our scenarios implies a number of very distinctive predictions at low energies. Because of the particular embedding of the SM characterizing these models an approximate gauge coupling unification of (all or just two of) the SM gauge couplings generically occurs at around the axion scale. From a bottom-up perspective this requires new light particles, potentially accessible to colliders, that make the SM gauge couplings unify at 10 10−11 GeV. Above that scale one discovers that the SM gauge fields are mostly the vectors of a grand-unified group, and yet the SM fermions do not fill complete GUT multiplets, so baryon number violating operators do not appear at that threshold. Rather, the actual lepton-quark unification occurs at much larger scales, as usual.
Another generic prediction of our scenarios is the existence of several light SM-neutral bosonic states beyond the QCD axion. Some of these can be potential dark matter candidates precisely in the regime f a ∼ 10 8−10 GeV expected in our models, i.e. when the QCD axion contributes negligibly. In the models of Section 3 these dark matter candidates are in the form of long-lived glueballs of an unbroken hidden gauge sector populated at re-heating, whereas in the models of Section 3.3 they correspond to GeV-scale NGBs with a cold abundance set by vacuum misalignment. In addition we always find unstable ALPs with anomalous couplings to the SM gauge bosons and masses of order a few TeV generated at the GUT scale.
In one way or another, all field theoretic explanations of the axion quality problem look a bit ad hoc. Our SO(10) models are no exception. Yet, taking into account the good quality of their U (1) PQ , the natural predisposition to accommodate grand-unified scenarios, and their distinctive phenomenological signatures, we think that our models do not compare too bad with the competition either. Perhaps the ideas developed in this paper will find a more compelling incarnation in the future.
for conversations on axions in string theory, M. Redi for sharing some of his expertise in composite axion models, I. Shoemaker for suggestions, and especially R. Contino, A. Podo, and F. Revello for comments on the draft and for discussions about their independent forthcoming paper on axions.