Unificaxion

Dark matter, gauge coupling unification, and the strong CP problem find a common and simple solution (in the absence of naturalness) within axion models. We show that such solution, even without specifying the details of the model implementation, makes testable predictions for the experimentally measurable axion parameters: the axion mass and its coupling to photons.


Introduction
Naturalness is considered the leading reason to believe that new physics beyond the Standard Model (SM) must exist at the weak scale. As no experimental evidence in favor of new physics at the weak scale has been observed yet, we explore the possibility of discarding the criterion of naturalness and following the lead of other arguments. Two interesting motivations to introduce new physics are Dark Matter (DM) and gauge coupling unification. Indeed these two arguments, in the absence of naturalness, have led to the intriguing hypothesis of Split Supersymmetry [1]. In this paper we will show that the same two arguments, together with the additional request of a solution to the strong CP problem, can lead to another interesting (and more minimal) hypothesis, which offers an experimentally testable prediction.
To account for DM we assume the existence of the axion, a particle which finds its motivation in the solution of the strong CP problem. Unlike the cases of the gauge hierarchy and the cosmological constant, there seems to be no anthropic explanation of the vanishingly small value of the QCD θ angle. Seeking a natural solution to the strong CP problem [2] while giving up naturalness on the other two problems appears therefore a logical option. A large class of invisible axion models (the so-called KSVZ [3]) make use of new matter, charged under color and PQ symmetry, with mass at about the same scale as the axion decay constant f a . It is possible that the new matter modifies the running of the SM coupling constants in such a way to achieve gauge coupling unification. We will show that this hypothesis, for which we coined the term unificaxion, leads to a prediction on the ratio between the axion coupling to photons and the the axion mass. Present and future axion experiments can test the hypothesis of unificaxion, although no new dynamics beyond the SM is predicted at the weak scale. This paper is organized as follows. In section 2 we discuss the requirements for new heavy fermionic particles to achieve unification of SM gauge couplings. These results are used in section 3 to predict the axion-photon coupling in unificaxion, as a function of the axion mass. We also extend our discussion to the case of supersymmetry. Finally our conclusions are drawn in section 4.

Unification with matter at an intermediate scale
Let us consider the existence of new fermions Ψ with common masses M Ψ . We consider fermions in real representations of the SM group (because we are interested in very massive particles) that can be embedded in SU(5) (because we have in mind a grand unified theory). We want to explore the restrictions on the possible quantum numbers of Ψ by requiring a  (5) multiplets. Thick red dots represent the contributions of single real representations, and the arrows with names represent the simplest of such cases, which are found inside the following SU (5) representations (see table 1): fundamental (+ conjugate) 5 ⊕ 5 = L ⊕ D, antisymmetric (+ conjugate) 10 ⊕ 10 = Q ⊕ U ⊕ D, symmetric (+ conjugate) 15 ⊕ 15 = Q ⊕ T ⊕ S, and adjoint 24 = 1 ⊕ V ⊕ G ⊕ X. Right: The green area shows the range of the differences between the beta function coefficients for the gauge couplings, ∆b 3 − ∆b 2 and ∆b 2 − ∆b 1 , that provide unification at a GUT scale between 10 16 and 10 18 GeV (red dotted lines) with intermediate scale indicated in blue (dashed lines). The point marked as represents the case of low-energy supersymmetry, and its nearest dot represents the case of split supersymmetry with new scalars at the unification scale.
successful unification of the SM gauge coupling constants at some scale M GUT in one loop approximation: Here b SM i = {41/10, −19/6, −7} and ∆b i are the contributions to the β-function coefficients due to SM particles and to new fermions associated with the axion sector, respectively.
Only the values of the differences ∆b i − ∆b j are relevant for unification and for fixing its scale, while adding a universal ∆b i (equal for any i) increases the value of α GUT . The arrows in fig. 1a show the values of (∆b 3 − ∆b 2 , ∆b 2 − ∆b 1 ) corresponding to real SM and SU(2) (n 2 ), and the index z for chiral irreducible representations of the SM group contained in the SU(5) representations 5 ⊕ 5, 10 ⊕ 10, 15 ⊕ 15, 24. The entries in the right-hand side give the contributions to ∆b i from fermions in real representations of the SM group, equal to the previous representations when they are real (Y = 0), or adding their conjugates when they are chiral (Y = 0). multiplets contained in SU (5) representations with low dimensions: 5 ⊕ 5, 10 ⊕ 10, 15 ⊕ 15, 24. The quantum numbers of these states are summarized in table 1. Within each complete SU(5) multiplet, arrows sum to zero. Next, by summing these arrows with integer non-negative coefficients we generate the most generic set of points produced by arbitrary combinations of these multiplets. The dots represent all possible combinations: we see that they form a sparse lattice, where each point can be produced in many different ways. In particular, arbitrary combinations of the SM representations contained in 5 ⊕5 and 10 ⊕ 10 are already enough to span the whole lattice, and nothing more is obtained by introducing the 24 or the 15 ⊕ 15.

Extension to arbitrary representations
Here we show that the lattice points shown in fig. 1a describe the most general case for fermionic matter and that no new points are added by including any other arbitrary irreducible representation R. In general, the representation R can be described by the SU(3) Dynkin label (n 3 ,n 3 ), the SU(2) Dynkin label (n 2 ), and hypercharge Y . The Dynkin labels count the differences between the number of boxes in successive rows of the corresponding Young tableau. So the indices n 3 ,n 3 , and n 2 are non-negative integers. The contributions of Ψ to the β-functions are where d is the dimensionality of R, given by The condition that R is embedded in an irreducible representation of SU(5) implies a constraint on the possible values of Y . By projecting the weights of a generic multiplet of SU(5) into the SM subgroup we find that the hypercharge Y must satisfy Here z is the component of a weight of the SU(5) irreducible representation, in the Dynkin basis, corresponding to the U(1) generator of the Cartan subalgebra. It can be proven that z can take all possible values in Z. As a result, ∆b i can only scan a discrete set of values, which is determined by eq.s (2)-(4), with n 2 , n 3 ,n 3 ∈ N and z ∈ Z. The points obtained by this procedure are shown in fig. 1a as red thick circles. It is easy to see that all these points can be generated by taking appropriate combinations of the representations contained in 5 ⊕5 and 10 ⊕ 10 of SU(5).

Unification
In fig. 1b we show the range of β-function coefficients compatible with unification of gauge couplings obtained by varying the intermediate mass M Ψ and the GUT scale M GUT as indicated. The range of acceptable M GUT is limited from below by the requirement of proton stability [4] An upper bound on M GUT of the order of the Planck mass comes from the request that unification be described at least in leading-log approximation within quantum field theory. Indeed the request M GUT ∼ < M Pl is also motivated by gauge/gravity unification. The green dot shown in fig. 1b represents the one-loop prediction of low-energy supersymmetry. Note that the green dot does not exactly correspond to any of the black dots because, in addition to the fermionic degrees of freedom of gauginos and higgsinos (G ⊕ V ⊕ H), low-energy supersymmetry introduces also some new scalars in incomplete GUT multiplets (the second Higgs doublet). Notice also that the green dot sits just outside the region where 1-loop unification is exactly achieved. This is a well know fact, but the mismatch is well within the size of plausible threshold corrections.  Table 2: Models of unificaxion with up to 3 fermion multiplets, intermediate mass between 10 3 and 10 14 GeV, and unification mass satisfying eq. (5). Their predictions for α GUT , M GUT , M Ψ , and E/N are shown.
particle content of the simplest models (containing at most 3 of the fermion representations listed in table 1) are summarized in table 2. These models are selected by requesting that 10 3 GeV < M Ψ < 10 14 GeV and that M GUT satisfies eq. (5). It is interesting that some of these models predict a value of M GUT close to its lower value and thus predict a rate for proton decay just beyond the present experimental sensitivity.
Furthermore, naïve gauge/gravity unification in 4 dimensions suggests the extra relation valid up to the model-dependent order-one factor k. This relation is shown as a green band in fig. 2b, where we consider the range 1 < k < 40 with the upper bound motivated by the heterotic string computation of [5]. It should be stressed that a value k 1 can be obtained in type I string theory or in M-theory as a result of a parametrically large volume of compactification. Therefore all values of M GUT ∼ < M Pl are in principle compatible with the unification of all forces, though M GUT on the high end does seem perhaps more plausible, in that it does not require additional very large (or very small) parameters.

Axions
In this paper we are focusing on KSVZ axion models [3] which, in addition to the SM particles, introduce a complex scalar A coupled to new Dirac fermions Ψ, in a representation of the gauge group Ψ = ⊕ r ψ r . Such models assume a Peccei-Quinn U(1) global symmetry which forbids the Dirac mass term MΨΨ but allows for Yukawa couplings r λ r Aψ r P L ψ r + A †ψ r P R ψ r .
The U(1) symmetry is spontaneously broken by the vacuum expectation value A = T 2 f a , where Tr T a T b = 1 2 T 2 δ ab for the QCD generators of the fermions Ψ, and f a is the decay constant of the light axion a = √ 2 Im A. We assume that all λ r couplings have comparable size, λ r ∼ λ Ψ . Consequently all fermions acquire roughly the same mass M Ψ = λ Ψ A , representing the only threshold between the weak and unification scales. Remarkably, as we shall discuss below, the overall size of λ Ψ does not affect our main prediction.
Non-observation of axion emission from stars and supernovae implies that the axion decay constant must be larger than f a > 10 9 GeV. Furthermore, requiring that the axion dark matter density generated by the initial misalignment mechanism [6], does not exceed the observed dark matter density Ω DM ≈ 0.23 implies the upper bound f a < 10 12 GeV, under the assumption that the axion vev a * in the early universe was of the order of f a [6]. Values of f a ∼ 10 12 GeV are therefore favored by the assumption that DM is made of axions. However, larger values of f a can be compatible with axion dark matter, if a * is sufficiently smaller than its natural value of order f a . Provided that λ Ψ = O(1), the heavy fermions Ψ that are associated with the axion dynamics have masses in the same range as f a but, for small λ Ψ , the intermediate scale M Ψ could be much less than 10 9 GeV, without conflict with axion bounds from stellar emission. Therefore, a wide range of intermediate scales M Ψ is compatible with axion DM.

Axion coupling to photons
The anomalous coupling of axion to photons is defined as whereF µν ≡ 1 2 µναβ F αβ . DM axion experiments are starting to probe the theoretically favored region of the axion mass m a and axion coupling to photons g aγγ : the ADMX experiment [7] obtained the limit |g aγγ | < 7 10 16 GeV as claimed in some analyses [8].
for m a in the range m a = 1.9 − 3.55 µeV. In eq. (11) we have shown the dependence on the local DM density, ρ DM ≈ 0.3 GeV/cm 3 , as unprecisely determined from rotation curves and halo dynamics. We used the most conservative limit corresponding to completely virialized axions in the galactic halo. Note that, in case of discovery, the uncertainty due to the galactic axion velocity distribution can be eliminated by studying the frequency dependence of the signal. A positive signal would allow to precisely measure m a and determine g aγγ . The experimental determination of m a and g aγγ can provide a crucial test of the idea of unificaxion, as we now discuss.
The measurable ratio g aγγ /m a is sensitive to the particle content of the theory: Here f π = 93 MeV and we defined E/N = r Q PQ q 2 / r Q PQ T 2 where the sums extend over all fermions ψ r with PQ charges Q PQ , electric charges q, and Tr T a T b = 1 2 T 2 δ ab for the QCD generators. As we already explained, we are considering models in which all fermions have the same PQ charges Q PQ . Adopting the quark masses given by chiral perturbation theory at lowest order, m u /m d = 0.56 [8], we get and fig. 3 shows the band induced by the uncertainty on quark masses. Concerning m u /m d , we recall that second-order effects in m q in chiral perturbation theory lead to an uncertainty on the value of m u /m d famously known as the Kaplan-Manohar (KM) ambiguity [8].
Conceivably the KM ambiguity could have made m u /m d = 0 compatible with experimental meson masses, thus disposing of the need for an axion. However that possibility is now disfavored by lattice simulations. Ref. [9] analyzed the issue in detail also considering the impact of the KM second-order effect in eq. (12). The uncertainty band shown in fig. 3 reflects the assessment of that study. It is also rather evident from eq. (12) that the leading source of uncertainty is given by the m d /m u term in the square root factor up front, which is singular as m u → 0. In that respect, the error in m s /m u = Q m 2 d /m 2 u − 1 has a minor impact and gives no additional uncertainty because Q = 22.7 ± 0.8 is precisely known and large.

Unificaxion prediction for g aγγ /m a
Coming to unificaxion, we find that the model dependent axion coupling coefficient in eq. (12) is related to the β-function coefficients as This result is pivotal for our analysis and so it is important to clarify the hypotheses upon which it rests. One simplifying characteristic of the class of models we are considering is, as previously mentioned, the existence of a single intermediate threshold at the energy scale M Ψ . Another condition is the existence of just one additional scalar describing the dynamics of PQ symmetry breaking and the physics at the intermediate threshold. Besides simplicity, the lack of multiple scalar fields may be justified in the multiverse by the condition of not exacerbating the naturalness problem. Note that scalars charged under the SM gauge group would affect gauge coupling unification in a way that is completely independent of axion couplings. Finally one last simplifying assumption, upon which our conclusions mostly rely, is that all fermions have the same PQ charge, as postulated in eq. (7). Algebraically this is the statement that the PQ charge matrix Q PQ is proportional to the identity and, in our normalization, Q PQ = 1. Because of that, the anomaly coefficients d QAB that control the effective coupling of the axion to the dual field strengths coincides with the contribution of fermions to the β function coefficients From this equation, the result in eq. (14) follows immediately. Note that with just one scalar A in principle one could accommodate fermions with PQ charge = −1 with Yukawa given by eq. (8), after the replacement A → A † . That would spoil the proportionality between contributions to the β function and the axion coupling, weakening our conclusions. This could not occur in a supersymmetric theory where, in the presence of just one scalar, holomorphy implies Q PQ = 1. In a non-supersymmetric theory the assumption Q PQ = 1 can be justified by assuming that all fermions Ψ sit in the same multiplet within a more fundamental description, but we will not try to construct explicit examples. It is also worth remarking that, under the assumption Q PQ = 1, the relation between the axion coupling and the β-function coefficients is preserved, regardless of the dynamics of the PQ breaking sector and, in particular, regardless of the number and PQ charges of the SM singlets in that sector.
Equations (13) and (14) provide the link between unification and axion phenomenology, which is the key feature of unificaxion. Gauge coupling unification selects a special range for ∆b i which, in turn, determines the measurable quantity g aγγ /m a . The prediction for g aγγ /m a is obtained only by the request of unification, with no need to specify the particular particle content of the model or their interactions. Fig. 4 shows the correlations between E/N , M GUT , α GUT , and M Ψ . In particular, fig. 4b illustrates how the prediction of unificaxion for E/N , which is directly related to g aγγ /m a through eq. (13), depends on the intermediate mass M Ψ . Under the simplifying assumption λ Ψ ≈ 1, M Ψ can be identified with f a and translated into a value of the axion mass m a (as shown in fig. 4b). However, the prediction E/N is independent of λ Ψ and thus more robust. Future experimental determinations of m a and g aγγ /m a will select a region in this plane, allowing for a test of unificaxion.
It is interesting that unificaxion allows for solutions with an intermediate mass M Ψ compatible with the favored range of f a , giving 1.0 < E/N < 2.5. Moreover, assuming that the unified coupling is large, as maybe suggested in simpler models of gauge/gravity unification, we obtain E/N > 1.6. Assuming that the unification scale is very close to the Planck scale leads to E/N < 2.1. A peculiar coincidence is that the predicted range of E/N is centered around the value for which g aγγ suffers a perfect cancellation. This is, of course, a worrisome result, because the axion would turn out to be literally invisible. Note also that E/N in unificaxion is always smaller than the value predicted by DFSZ axion models with fermions in complete SU(5) multiplets [10], which is E/N = 8/3, see eq. (14). Nevertheless, E/N = 8/3 gives a value of the axion-photon coupling which is inside the range predicted by unificaxion (as evident from eq. (13) and fig. 3, |g aγγ | turns out to be the same for E/N = 8/3 and E/N = 1.2). Thus, the DFSZ axion or any axion model with fermions in complete GUT multiplets are experimentally indistinguishable from unificaxion, unless one devises an experiment sensitive not only to the axion coupling to photons, but also to its coupling to hadrons.
The ADMX bound in eq. (11) is not yet strong enough to constrain unificaxion. However, an improved sensitivity is considered feasible, and furthermore new experimental techniques based on cold molecules could allow to probe also lighter axions with f a ∼ M GUT [11].

The supersymmetric case
It may be of interest to extend to the case of supersymmetry our considerations about the consequences of unification for the axion coupling to photons. Let us consider low-energy supersymmetry, with new matter in chiral supermultiplets at an intermediate scale M Ψ , associated with the dynamics of the axion sector. In this case, eq. (14) remains valid because heavy fermions and scalars contribute to ∆b i in the same proportion:  fig. 1, but now in the supersymmetric case. The region that allows unification collapses to a narrow strip because unification is already realized in low-energy supersymmetry without any extra field.
qualitative difference because, as well known, in low-energy supersymmetry unification is already achieved, at M GUT ≈ 2 × 10 16 GeV, without the need for extra heavy multiplets. As a consequence the most plausible scenario for the axion sector is to contain only complete SU(5) multiplets. In this case (∆b 1 = ∆b 2 = ∆b 3 ), eq. (14) gives the welldefined prediction E/N = 8/3. This value is larger than what expected in unificaxion but, as mentioned previously, it is indistinguishable through measurements of the effective axion-photon coupling. However, it is also possible that new matter at the scale M Ψ modifies the gauge coupling evolution, achieving unification at a scale M GUT different than the usual value of low-energy supersymmetry, 2 × 10 16 GeV. A solution of the one-loop renormalization-group equation shows that this happens when and that the new unification scale is given by where we worked in the limit of exact unification for low-energy supersymmetry. This is illustrated by the results shown in fig. 5, which is the analogous of fig. 1b in the case of supersymmetry. The region of ∆b i compatible with gauge coupling unification essentially collapses to a line, approximately described by eq. (17). A simple example of modified unification in low-energy supersymmetry is the addition of heavy chiral supermultiplets in the adjoint in the SM gauge group (G ⊕ V ), a case motivated by partial N = 2 supersymmetry. In this case, we find that the intermediate scale M Ψ is related to the unification scale by We also obtain E/N = 2/3, a value smaller than what expected in unificaxion (and thus a larger value of |g aγγ |). Another example, which gives a more accurate unification of gauge couplings, is the addition of chiral supermultiplets in the representation G ⊕ Q. This gives and, again, E/N = 2/3.

Conclusions
The exploration of the 'multiverse' prompted theoretical physicists to revisit their belief in the naturalness criterion and, in some cases, to replace it with a biased statistical approach. Although interesting from a theoretical point of view, this new paradigm suffers from a chronic lack of experimental predictions. For instance, when applied to the hierarchy problem, it leads to the conclusion that it is perfectly acceptable for the Higgs boson not to be accompanied by other new particles or new dynamics at LHC energies. This conclusion appears disheartening and it can be hardly used as evidence for the multiverse. The multiverse hypothesis is desperately in search for observational tests. One rare exception of this lack of experimental consequences is offered by Split Supersymmetry [1]. In this context, new physics should be present at the weak scale, not because of naturalness, but because of DM and gauge coupling unification. Discovery of a long-lived gluino at the LHC would provide crucial confirmation of this hypothesis. The Higgs searches have already narrowed down the possible range of the energy scale of Split Supersymmetry [12], stating that the gluino lifetime must be τg < ( TeV/mg) 5 10 −4 seconds. Future searches at the LHC with 14 TeV will probe the existence of a metastable gluino up to masses of about 2.5-3 TeV.
In this paper, we suggested an alternative approach -unificaxion -for predicting new physics without invoking naturalness. The idea of unificaxion is to introduce an invisible axion model with a single complex scalar field, which spontaneously breaks a PQ abelian symmetry at a scale f a and gives masses to a set of fermions. These fermions contribute to the QCD θ term and have the appropriate quantum numbers to achieve unification of the SM gauge couplings. In this way, the strong CP problem, DM, and unification find a common solution with new physics occurring at an intermediate scale M Ψ . While the axion decay constant f a must be larger than about 10 11 -10 12 GeV to account for DM, the intermediate scale M Ψ of the new fermions could conceivably take even smaller values. Nonetheless, the scale M Ψ is expected to be far from the Fermi mass. Although unificaxion gives no new physics at the LHC, it makes one experimentally testable prediction in the context of axion searches. The axion-photon coupling is determined by the particle content that is responsible for gauge coupling unification. In particular, we have found that 1.0 < E/N < 2.5, implying that |g aγγ | m a < 2 × 10 −16 GeV −1 µeV −1 .
This prediction is further narrowed down by the hypothesis that gauge unification is accompanied by an approximate gauge-gravity unification, which implies 1.6 < E/N < 2.1 and thus |g aγγ | m a < 1 × 10 −16 GeV −1 µeV −1 .
Present axion DM experiments are already probing an interesting range of parameters and future improvements can reach the sensitivity to test the prediction of unificaxion. However, it should be stressed that the discovery of an axion that satisfies eq. (21) would only provide an indication in favor of unificaxion, but not a definitive confirmation. For example, any axion model in which the new fermions form complete GUT multiplets predicts E/N = 8/3, leading to a value of |g aγγ | inside the range of eq. (21). Only a test of the hadronic coupling of the axion could disentangle the two cases. Moreover, a worrisome feature of the unificaxion prediction is that the axion-photon coupling g aγγ could be vanishingly small due to a fortuitous cancellation between short-distance and long-distance contributions. If this is the case, unificaxion would not be experimentally testable and it would remain forever buried in the obscurity of the multiverse.