Color unified dynamical axion

We consider an enlarged color sector which solves the strong CP problem via new massless fermions. The spontaneous breaking of a unified color group into QCD and another confining group provides a source of naturally large axion mass ma\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_a$$\end{document} due to small size instantons. This extra source of axion mass respects automatically the alignment of the vacuum, ensuring a low-energy CP-conserving vacuum. The mechanism does not appeal to a Z2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_2$$\end{document} “mirror” copy of the SM, nor does it require any fine-tuning of the axion-related couplings at the unification scale. There is no very light axion and uncharacteristically the lighter spectrum contains instead sterile fermions. The axion scale fa\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_a$$\end{document} can be naturally brought down to a few TeV, with an exotic spectrum of colored pseudoscalars lighter than this scale, observable at colliders exclusively via strong interactions. The {ma,fa}\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{m_a, f_a\}$$\end{document} parameter space which allows a solution of the strong CP problem is thus enlarged well beyond that of invisible axion models.


Introduction
Phenomenological analyses based on chiral perturbation theory and supported by lattice computations indicate that all Standard Model (SM) quarks have non-zero masses. This disfavors the solution to the strong CP problem via one massless SM quark, which automatically guarantees a U (1) axial invariance at the classical level. The interesting possibility of having a massless up quark in the microscopic theory which appears as massive at QCD scales due to non-perturbative instanton contributions [1] does not seem to be realized in nature, even if this option is not completely excluded [2][3][4][5][6][7].
It is still possible to solve the strong CP problem using massless fermions if the SM up quark is massive. The idea is to enlarge the SM gauge group with a new confining sector [8][9][10], whose scale is much larger than that of the QCD group, SU (3) c . Extra massless quarks charged under both QCD and the new confining sector may realistically solve the problem [11]. A new spectrum of confined states results.
The term "axion" denotes any (pseudo)Goldstone boson (pGB) of a global chiral U (1) symmetry which is exact at the classical level but has anomalous couplings to the field strength of a confining group. 1 Axions are characteristic of solutions to the strong CP problem based on an anomalous U (1) axial symmetry, usually called Peccei-Quinn (PQ) symmetry [12]. When the number of axions in a given theory outnumbers the total number of distinct instanton-induced scales other than QCD to which they couple, one (or more) light axions remain. Axions which are not elementary but composed of fermions are typical of theories with an enlarged confining sector and are often referred to as "dynamical" or composite axions. Very heavy dynamical axions made out of massless fermions will typically acquire the bulk of their mass from the largest instanton-induced scale to which they exhibit anomalous couplings, For a very light axion coupled only to QCD, the mixing with the η pseudoscalar is relevant instead, and that axion obeys the usual relation [13,14] where m π , f π , m u , m d denote the pion mass and coupling constant, and the up and down quark masses, respectively. The first step in the direction of solving the strong CP problem with exotic massless quarks was the proposal by K. Choi and E. Kim [11] to enlarge the confining gauge sector of the SM to SU (3) c × SU (Ñ ), with the latter having a scale larger than that of QCD,˜ QCD . Two confined charges would then exist in nature, color and axicolor respectively, and correspondingly two distinct sources of instanton potentials. A massless color-triplet quark Q, charged also under axicolor and singlet under the SM electroweak symmetry, would solve the QCD strong CP problem. The fact that the axicolor scale is very large would explain the non-observation of exotic bound states at low energies. An issue arises because there are now two potentially harmful vacuum angles to absorb: θ c of QCD andθ . Only one combination of them would be redefined away by a chiral rotation of Q. This was easily remedied by adding a second exotic quark χ charged only under SU (Ñ ).
In the limit of vanishing QCD coupling, the SU (Ñ ) sector described four flavors which seed two singlet pseudoscalars with anomalous couplings: two dynamical axions. Finally, taking into account the SM quark sector and thus the SM η , three flavor-singlet pseudoscalars result in the low-energy spectrum for only two instanton sources of masses: the η , a very heavy axion with mass ∼˜ and a second axion almost massless and obeying Eq. (2). Because of this last axion, the axicolor construction can be seen as an ultraviolet dynamical completion of the invisible axion paradigm. As usual, a very large f a scale is required to be orders of magnitude larger than the electroweak (EW) one, albeit with the advantage of being free from scalar potential fine-tunings.
In a different and recent attempt [15] to solve the strong CP problem with extra massless quarks, the same SU (3) c × SU (Ñ ) confining sector is considered. No light axion remains in the low-energy spectrum, though, as only two pseudoscalar mesons are present which couple to the two anomalous currents: the customary η meson and one axion. This is achieved by assuming only one exotic massless quark Q instead of two. In the limit of vanishing QCD coupling, the SU (Ñ ) sector then describes only three flavors, resulting in only one gauge singlet pseudoscalar with anomalous SU (Ñ ) couplings: a dynamical axion. Both the η and the axion thus acquire a mass, and the axion mass is induced by SU (Ñ ) instantons, m a , f a ∼˜ . For such a heavy axion the f a scale can be as low as the TeV range without incurring unacceptable phenomenological consequences [15][16][17][18][19][20][21][22][23][24]27]. The issue of the two θ parameters is solved in this proposal by imposing a discrete Z 2 symmetry relating the two sectors. Unfortunately, the practical implementation of this idea requires a complete Z 2 -"mirror" copy of the SM. This Z 2 symmetry is explicitly broken by a scalar potential which gives the second Higgs field a very large vacuum expectation value, e.g. 10 14 GeV, in order to sufficiently modify the running of the two confining scales. This is overall a quite complex, tuned, and large structure.
Even more recently, the SM massless quark avenue has been revisited in a theory that involves a product of SU (3) groups which break spontaneously to QCD [22]. A very interesting aspect developed by the same authors in a previous work [23] is the impact of small size instantons of Table 1 The massless fermion sector of the SU (6) construction above the unification scale SU (6) SU (2) L U (1) L 20 1 0 the group which undergoes spontaneous symmetry breaking. These instantons are shown to provide a possible extra source of large masses for the putative axions of the model. We will develop in what follows a new solution to the strong CP problem via massless fermions, in which the issue of the different θ parameters that arise in the presence of two or more confining groups is solved via color unification. Color unification with massless quarks is attempted here for the first time. This path is an alternative to the axicolortype constructions and will lead to different phenomenology. QCD will be unified with another confining sector singlet under the electroweak gauge symmetry. The color unified theory (CUT) breaks spontaneously to QCD and another confining group. The small-size instantons of the unified color group provide an extra source of high masses for the axions of the theory, and it will be shown that typically no axion remains at low scales. The exotic low-energy spectrum is instead fermionic. Furthermore, it will be shown that interesting new phenomenological signals can be explored at colliders. The complete ultraviolet completion of this idea will be developed, implementing two different scenarios: in one of them the two resulting heavy axions are dynamical, while in the other one axion is elementary.
The structure of the paper can be easily inferred from the Table of Contents.

SU(6) color unification
We propose a scenario in which QCD is unified with another confining group into SU (6), and a single, strictly massless SU (6) fermion rotates away simultaneously all θ parameters. The unification path in the context of an extended strong sector to solve the strong CP problem was first proposed by Rubakov long ago [24], in a Grand Unification construction that relied on traditional models à la DFSZ [25,26] with massive exotic fields, and required a Z 2 mirror copy of the complete SM field content. Another recent attempt [27] using unification ideas also relied on massive exotic fermions à la DFSZ. Here we instead consider color unification in the presence of massless fermions. The massless SU (6) fermion belongs to the 20 representation of the SU (6) CUT, having a definite chirality (e.g. left-handed) while being a singlet of the SM SU (2) L × U (1) Y gauge symmetry, see Table 1.
At a color unification scale CUT much higher than the EW one, the SU (6) group breaks into (3) Table 2 The massless fermion sector of the SU (6) construction below the unification scale. The notation is such that The parameters θ c of SU (3) c andθ of SU (3) are necessarily equal and unphysical down to the unification scale, and will remain so even below the unification scale as long as remains massless, protected by chiral symmetry. Under spontaneous symmetry breaking of the CUT symmetry, decomposes as where the charges under the U (1) group in Eq. (3) are shown in parenthesis for completeness. If the components of the L field are to remain massless under the CUT scale, SU (3) must confine. The two confining scales QCD and˜ need to be separated with˜ QCD , as no bound states are observed other than those compatible with QCD. The 20dimensional representation is thus advantageous because all its components charged under QCD are also charged under SU (3), and so will form bound states at the higher scale˜ . This representation is also pseudo-real, and so the theory is anomaly free. The non-trivial issue of how to separate˜ and QCD is discussed further below. The colored-axicolored massless fermions in Eq. (4) will be denoted ψ L ,R , see Table 2, while ψ ν will refer to the singlet massless fermions to convey that they act like sterile 2 neutrinos. The ψ ν fields only connect to the other fields through the unified strong forces, and thus their couplings to the visible universe will be safely suppressed by CUT , provided the U (1) gauge group in (3) is also broken near that scale. 3 SU (6) color unification is thus a successful path to solve the strong CP problem, and this fact will remain at the heart of the developments in this paper. The remaining problem is to obtain a low-energy spectrum which is fully compatible with observations.

The SM fermions
Because of color unification, the SM quarks must belong to SU (6) multiplets. The simplest option is to include them in six-dimensional fundamental representations. For each fermion generation, where q L , u R and d R denote the SM quarks, while their SU (6) partners are signaled by tildes. Theq L fields are necessarily electroweak doublets, and this character turns out to be the major practical issue of this model: -Leaving the tilde-quark sector massless but confined is unacceptable, as the condensate -assuming chiral symmetry breaking of SU (3) -would typically break the SM EW symmetry at the large˜ scale. -Alternatively, giving much larger masses (≥˜ ) to the tilde quarks is not viable either in this SU (6) setup without spoiling SM quark masses, since they belong to the same multiplet. If a scalar field gave high masses to the tilde quarks 4 by obtaining a high vacuum expectation value (vev), that scalar field would have to be an SU (2) L doublet. Then its large vev would spontaneously break SM EW symmetry, giving gigantic masses to the W and Z boson.
The main problem of this model is then the unacceptably light tilde-fermion sector. We will develop next an extension whose only purpose is precisely to achieve high masses for the tilde-sector quarks, decoupling them from the low-energy spectrum. By the same token, the necessary separation of QCD and a larger confining scale will naturally follow. 5 We will develop in detail two realistic ultraviolet (UV) completions.

The realistic color unified theory: SU(6) × SU()
It is necessary to give large masses to the tilde-quark sector without giving masses to the SM quarks, a challenging enterprise as explained above due to the SU (6) unification. An external mechanism is ideal for this task. The color unified SU (6) group which contains QCD is enlarged via an external non-abelian SU (3 ) group with additional fermions charged only under the latter. In fact, all fermions in the theory will be charged under only one of the two groups, SU (6) or SU (3 ). L will be thus taken to be a singlet of SU (3 ) and the same applies to the multiplets in Eq. (5) which contain the SM quarks. The two sectors are connected exclusively via a new scalar . QCD remains a subgroup of SU (6), whose θ -parameter is rotated away by the massless L fermion in Eq. (4). This type of auxiliary extension was suggested in Ref. [27] to give high masses to exotic fermions in a different context. The field content of our model is summarized in Table 3, in which all fermions except L will become massive. It is easy to see that the theory with this matter content is anomaly free. The scalar appearing in the table belongs to the bifundamental of SU (6) × SU (3 ), and its vev breaks color unification at a scale CUT , taken to be much larger than all SM scales, The fermion quantum numbers under the two resulting groups are also shown in Table 3.
A simple CUT-invariant Yukawa Lagrangian which connects the SU (6) and the auxiliary SU (3 ) extension reads The CUT symmetry is spontaneously broken upon taking a vev of the order of the CUT breaking scale CUT Unless otherwise stated, we will assume in what follows that all κ i Yukawa couplings are O(1), meaning all tilde and prime fermion masses are of order CUT . Some tuning of the κ i values could be acceptable, though, as discussed further below. The SM fermions get their masses through the usual SM Higgs doublet , which in this model is a singlet of SU where Y SM i denote the SM Yukawa couplings. Analogously, the most general Lagrangian compatible with all symmetries discussed above allows us to write Yukawa couplings of the Higgs field to the prime-sector fermions, 6 Note that we take to be real. The phases of the nonvanishing entries in Eq. (8) can all be made equal by an SU (6) × SU (3 ) transformation; the remaining phase can be removed by a transformation under the U (1) defined in Eq. (12). 7 In this notation taken from unified models the contraction of the spinor indices is implicit, more precisely the first term would read Q T L C U c L , where C = iγ 2 γ 0 is the charge conjugation matrix.
Eqs. (10) and (11) induce contributions to the tilde fermion masses which are quantitatively irrelevant in comparison with those from Eq. (9). In addition, the couplings in Eq. (11) will be absent for symmetry reasons in one of the models to be developed in this paper (Model II in Sect. 3.2). Both SU (3) c and SU (3) diag can now remain unbroken and confine at two different scales, QCD and diag , with diag QCD . The task of achieving different values for the two confining scales and getting rid of the tilde sector or any other dangerous exotic sector is thus accomplished.
Note that the Yukawa-type Lagrangian in Eq. (7) has an inherent global U (1) symmetry under which only the prime fermions and the field would transform -a generalized Baryon number symmetry in the prime sector, with charges This symmetry is not chiral, thus not anomalous under SU (3 ), and irrelevant to the strong CP problem. An associated pGB 8 results after spontaneous breaking, albeit with its couplings safely suppressed by the CUT scale and interacting only with the very heavy prime and tilde sectors. 9 Finally, the theory below CUT contains phenomenologically interesting bound states formed from the massless ψ L ,R fermions, to be studied below. The spectrum of free eigenstates below the EW scale contains the usual SM spectrum, plus a harmless pGB and sterile neutrinos. 8 This symmetry is broken at loop level by SU (2) L sphalerons, in the same way that in the SM baryon number current is anomalous. For our purposes this effect is negligible. 9 As suggested in Ref. [23], this type of pGB could be entirely removed by gauging the U(1) group. There is no real need to implement this procedure in our case, though, given the strongly suppressed couplings of this pGB.

θ issue
The extension of the strong sector by the auxiliary external group SU (3 ) brings a new θ parameter into the game: where G i denote gauge field strengths with tensorial indices omitted. G c , G 6 , G and G diag correspond respectively to the SM QCD gauge group, SU (6), SU (3 ) and SU (3) diag . While the rotation of the massless field was designed to reabsorb θ 6 and ultimately θ c , θ may source back a SM strong CP problem through the contamination to the visible sector via the scalar. Indeed, at low energies the massless quark ψ transforms as a (3,3), therefore the phase θ 6 cannot be fully reabsorbed in the Lagrangian since the chiral rotation that removes the SU (3) c θ -term generates a new contribution to the SU (3) diag topological term. Reference [27] acknowledges this issue (in the context of a different model which does not rely on massless fermions) and leaves it unsolved hoping that some UV completion solves it. In what follows, we will determine and exhaustively analyze two UV solutions, via the simple addition of either -An extra massless fermion transforming only under SU (3 ). -A second bifundamental scalar field, which automatically endows PQ invariance to the above extension procedure.
The first solution is more in line with the spirit of the present paper, as all θ parameters inducing a strong CP problem are made unphysical via massless fermions, and it is developed next.

Model I: adding a massless fermion charged under SU (3 ).
The θ parameter of the auxiliary SU (3 ) gauge group can be made unphysical by the addition of a massless fermion field χ that transforms as a fundamental of SU (3 ) and is an EW and SU (6) singlet. In other words, the field content for this solution is that previously shown in Table 3 plus the massless fermion χ with quantum numbers shown in Table 4. Additional composite bound states will result from χ , among them composite pseudoscalars with anomalous couplings -dynamical axions -whose masses are discussed further below.

Running of the coupling constants
The CUT breaking pattern in Eq. (6) imposes the following relations among the gauge couplings with the constraint where α c , α diag , α and α 6 denote respectively the coupling strength of QCD, SU (3) diag , SU (3 ) and SU (6). As shown in Fig. 1, there is a discontinuity in the running of the coupling constants at the CUT-breaking scale that allows α to have large values while reproducing the known QCD running at low scales. 10 Those α values will seed a source of large axion masses, as discussed in Sect. 3.1.3 further below. Although the relation in Eq. (14) imposes α diag ( CUT ) < α c ( CUT ), the presence of the SM q L , u R , and d R quarks at energies well below CUT slows down the running of QCD with respect to that of SU (3) diag . In this regime ψ and χ are the only fields left charged under SU (3) diag [theq sector generically decouples as their mass scale is set by CUT , see Eq. (9)]. As a consequence, α diag runs faster and thus the SU (3) diag group confines at a higher scale than QCD , see Fig. 1. This mechanism easily achieves the separation of the two confining scales. We computed both the one-and twoloop running and the latter actually reinforces the pattern, as illustrated in the figure for the choice diag = 4 TeV. Lower values of diag are also phenomenologically acceptable, see Sect. 4.1 below.

Confinement of SU (3) diag and pseudoscalar anomalous couplings to the confining interactions
At the scale diag , SU (3) diag confines and the remaining massless fermions will form massive QCD-colored bound states.
Assuming that both ψ and χ in Table 5 remain massless after SU (3) diag confinement, the SU (3) diag Lagrangian exhibits, at the classical level and in the limit in which α c is switched off, a global flavor symmetry U (4) L × U (4) R −→ U (4) V . 11 The chiral symmetry is then spontaneously broken by the quark condensates ψ L ψ R and χ L χ R . This results in 16 (p)GBs, decomposed here in terms of their QCD charges. There is a QCD octet plus a singlet with flavor contentψψ (3 c × 3 c = 8 c + 1 c ). The two QCD triplets, 3 c and3 c , correspond to the combinationsψχ andχψ. Finally, a color-singlet composite state is made out ofχχ. The fourteen colored mesons in Eq. (16) acquire large masses induced by gluon loops that are quadratically divergent and therefore sensitive to the cutoff scale diag , 12 The remaining two QCD singlets will be denoted here by η ψ and η χ and are shown next to be dynamical axions. The associated currents are where f d denotes the SU (3) diag pGB scale, with diag ≤ 4π f d . These classically conserved currents are broken at the quantum level by the SU (6) and SU (3 ) instantons, and so the currents are anomalous. The anomalous terms are These anomalous terms modify the classical equations of motion of the η ψ and η χ , 11 U (4) V remains unbroken and contains as a subgroup the SU (3) c QCD gauge group. 12 They contribute to the running of the QCD coupling constant, but given their high masses their impact is unnoticeable. and give rise to an effective Lagrangian, η ψ and η χ are thus two dynamical axions. It is to be stressed that the PQ scale in this model is f d ∼ diag for both axions and not the much larger CUT scale. When the SM quarks are taken into account, the η QCD pseudoscalar meson is also present at energies below the QCD confinement scale, and the effective Lagrangian of anomalous couplings reads where QCD ≤ 4π f π .
The scale E QCD can be expressed in terms of QCD observables from the chiral effective Lagrangian according to the results from Refs. [28][29][30] obtained in the large N limit. In the two-quark approximation, E QCD reads and E diag is obtained by rescaling the previous value As a consequence, the two instanton-induced scales QCD and diag provide a contribution to the masses of the pseudoscalars which have anomalous couplings. It follows that there are only two sources of mass (disregarding corrections from SM quark masses) for three states coupling to anomalous currents: η QCD , η ψ and η χ . In the absence of supplementary mass sources, one axion would get a mass of order diag while another one would have remained almost massless, see Eqs. (1) and (2), as often happens in models with dynamical axions. The model would be simply an ultraviolet implementation of the invisible axion paradigm.
As we will see next, an additional and important instanton source of mass for the axions is present, though, which lifts the light axion mass. In fact, depending on the model parameters, the χ fermion mass itself may: (1) still be lighter than diag , in which case the above discussed U (4) pattern of global symmetry holds; (2) alternatively, become more massive than diag and thus be decoupled from the spectrum above diag . In the latter case, the approximate global flavor symmetry of the SU The mass of this QCD-colored octet pseudoscalar is as previously given in Eq. (17) above, while the η χ would then disappear from the spectrum at energies above diag : the presence of η χ is to be then disregarded in Eqs. (19)-(26).

The impact of small-size instantons on the dynamical axion mass
An additional and putatively large contribution to the axion mass(es) applies in the presence of a spontaneously broken theory: the small-size instantons (SSI) of the theory at the breaking scale, as pointed out long ago in Refs. [31][32][33] and very recently in Ref. [23]. SSI can induce a large mass even for perturbative theories if the breaking scale is large enough to overcome the exponential suppression of instanton effects. In our model, the instantons of the color-unified theory in Eq. (6) near the CUT scale provide automatically this third source of axion mass. The SU (6) SSI can be neglected because of the smallness of α 6 at CUT (e.g. see Fig. 1) and the analysis below will focus on the SU (3 ) SSI contribution.
At the scales where we will compute the SSI effects, the SU (3 ) gauge coupling is perturbative. Therefore, for these instantons, the dilute gas approximation [34] gives a reliable estimate of the effective potential for the pseudoscalars. This was not the case for the previous instanton effects which, corresponding to the confinement scales, required the use of the effective chiral Lagrangian to obtain the potential from QCD observables.
It is well known [34][35][36][37] that, in the absence of fermions, the effective Lagrangian that describes instanton configurations for a pure Yang-Mills theory SU (N c ) induces a scale inst in the instanton potential given by where ρ is the instanton size, D[α ] is the dimensionless instanton density, The constant C inst reads [38,39] and the function c(x) is defined in Ref. [35] such that c( 1 /2) = 0.145873 and c(1) = 0.443307. For the SU (3 ) instantons of our model C inst = 0.0015. 13 In order to compute the integral in Eq. (30), the running of the coupling constant α (μ) must be included. At one loop this reads where α CUT ≡ α ( CUT ) and b is the one-loop β-function coefficient. For the spontaneously broken theory, only the SSI instantons with size ≤ 1/ CUT are relevant [40], This has the form where f (α , b) is given by For instance for the benchmark value α CUT = 0.3, the value of SU (3 ) SSI-induced scale in the absence of fermions (for which b = 10) is

Small-size instantons with small Yukawa couplings
The massless χ fermions may acquire an effective mass due to the instantons of SU (3) diag and SU (3 ) SSI, similar to the effective mass in QCD for a hypothetically massless SM up quark. In this section, we assume the y i couplings in Eq. (11) to be small, e.g O(0.2) for illustration. This is a regime of moderate small-size instanton effects, in which the instantoninduced effective χ mass is smaller than the confinement scale. As a consequence, a composite axion made out of this fermion can be considered to be a pGB. 14 The impact of the y i Yukawa couplings is illustrated by the one-instanton "flower" contribution in the left side of Fig. 2. The pure gauge results in Eq. (30) are now suppressed by two factors: the χ chiral condensate and the y i dependence, χχ is the order parameter controlling SU (3) diag chiral symmetry breaking and thus expected to be χχ − 3 diag . For ρ ≤ 1/ CUT , the product ρ 3 χχ 1 reduces the SSI-induced scale by orders of magnitude, with where, in the presence of N f Dirac fermions, The integral in Eq. (39) can be computed exactly, although a good estimation follows from the approximation which leads to 14 In Sect. 3.1.5, we will consider y i couplings of O(1) which corresponds to the regime of large small-size instanton effects, translating to a very heavy χ fermion which decouples from the spectrum well above the confinement regime.
The complete computation including fermions can be compared with the one in in Eq. (37): a strong suppression by a factor of order ( diag / CUT ) 3 is now present.
There is an additional contribution to the SSI scale independent of the the y i couplings, but suppressed by the χ chiral condensate, times the product of κ i Yukawa coupling of the prime-fermion sector in Eq. (7) and the product of SM Yukawa couplings. This contribution is illustrated by the instanton "double flower" in the right side of Fig. 2, and given by where the power of the 4π factor results from the 6 SM Yukawa couplings and the 12 κ i couplings in the product. For ρ ≤ 1/ CUT , this contribution is well approximated by In summary, putting together the two instanton contributions discussed, the SSI scale is given by Even with O(1) κ i Yukawa couplings, the second term is strongly suppressed by the SM Yukawa couplings. Overall, the size of the new scale SSI is quite sensitive to the value of the SU (3 ) coupling constant at the CUT-breaking scale.  Yukawa couplings in the prime sector, the diagram to the left represents the leading order contribution, whereas the one to the right is subdominant since it is further suppressed by the SM Yukawa couplings A mass is thus generated for the η χ axion, given by where the replacement diag 4π f d has been used. The lighter axion may thus be as light as O(GeV). The other dynamical axion of the theory has been shown to acquire a mass of order diag , see Eq. (26). Both dynamical axions have thus acquired masses far above the typical invisible axion regime, as a direct and unavoidable consequence of the instanton potentials inherent to the theory.

How light can the axion that couples to SSI become?
The y i values are very relevant for the size of SSI scale. Nevertheless, the mass spectrum and thus the running of coupling constants is basically unaffected by them. Should the y i couplings be negligible, SSI would be determined by the second term in Eq. (46). For vanishing y i values and generic κ i couplings of O(1), the product of Y SM i /4π factors in the second term in Eq. (46) would suppress the η χ axion mass to values of order keV. Such low masses are excluded up to the GeV range [41] for axion scales not far from a TeV, as is the case here. The allowed range is illustrated in Fig. 3. 15

Solution to the strong CP problem
It is pertinent to briefly re-check the status of the strong CP problem after taking into account the impact of the SSI of the spontaneously broken symmetry discussed above. Any new mass term for the axions breaks the PQ symmetry and therefore perturbs the axion potential; it is then important to verify that the vevs of the axions remain in the CP-conserving minimum, solving the strong CP problem. Indeed, this is the case with our color-unified proposal as, according to Eq. (13), the potential including θ i dependencies explicitly reads at second order, for the case discussed above in which both dynamical axions are present at the SU (3) diag confinement scale, For this potential, the minimum is CP-conserving: since all θ i dependences cancel. A word of caution is pertinent as the exact dependence of the potential on the phases of the different couplings in the Lagrangian which participate in fermion mass generation (κ i , Y SM i , y i , . . . ) remains to be computed. Nevertheless, the two massless fermions and χ guarantee that at energies above CUT the two parameters θ 6 andθ are unphysical. Below CUT, the spectrum of the theory is exclusively the SM one plus massless fermions, and the EW SM contributions are known to be negligible [44], even if a mismatch remained in spite of the low-energy presence of the massless quarks. An explicit computation of the threshold effects is elaborate and it is left for future work. Note that without the presence of the second PQ mechanism, that is, without the presence of the η χ field and its vev, all θ i in Eq. (13) would not have been reabsorbed, while Eq. (50) demonstrates that its inclusion does ensure a CP-conserving minimum.
It is very positive that in this model there is no contribution to the EW hierarchy problem coming from axion physics. No potential connects the EW and axion scales: the PQ scale f a 16 is set by diag and not CUT , and all axions are dynamically generated. This is a feature that our Model I shares with the original axicolor model, and in general with models of composite dynamical axion(s). There remains instead the customary fine-tuning in spontaneously broken unified theories, as CUT and the EW scale are connected via the scalar potential, but the latter does not communicate to our PQ mechanism.
The demonstrated possibility of lowering the PQ scale towards the electroweak one raises the question of the compatibility of the setup presented here with attempts to solve the EW hierarchy problem, e.g. via compositeness or supersymmetry. This is an interesting question which deserves future work. One could probably build supersymmetric or techincolor version of the models presented here. Nevertheless, the need to strongly separate the CUT scale from the EW scale is a non-trivial source of instability in the scalar potential.
Computation of the pseudoscalar mass matrix: η χ , η ψ , η QCD and light spectrum After the replacement of the pGBs with anomalous couplings by their physical excitations, η χ −→ η χ + η χ , η ψ −→ η ψ + η ψ , the effective low-energy potential for the axions and the SM η QCD field is given at second order by (disregarding the effects of SM quark masses) Expanding to second order in the fields yields the following mass matrix: Assuming the hierarchy of scales As advertised, the usual QCD η QC D phys remains as a light eigenstate, while one of the two composite axions acquires a mass in the range GeV to tens of TeV. The other composite 17 The physical states correspond to the following combinations: η ψ, phys 1/ f d (2 η χ + √ 6 η ψ ) and η χ, axion will have a mass orders of magnitude larger and will be out of collider reach.

Low energy spectrum and observable effects
In addition to the dynamical axion with mass O (GeV) or larger, under a few TeV the spectrum of the theory contains: -The SM pseudoscalar meson η QC D phys , plus the rest of the SM hadronic spectrum. -The exotic QCD-colored "pions" -color octets and color triplets -whose masses are given in Eq. (17) as With masses naturally lighter than the TeV scale, these QCD-colored pions can be easily produced at the LHC. -The two sterile fermions stemming from the 20-representation . They are basically invisible as their interactions with the visible world are suppressed by CUT , which is much larger than diag without any tuning. -Possibly, a GB associated with generalized baryon number. This GB is harmless as its interactions are suppressed by CUT . It can also easily be made arbitrarily heavy by gauging that global symmetry.
The very interesting phenomenological bounds and detection prospects for the exotic QCD-colored mesons will be quite similar to those applying to the next model. The ensemble will then be briefly developed in Sect. 4 further below. The same applies to the cosmological consequences of the two color-unified UV completions developed in this paper.

Small-size instantons with O(1) Yukawa couplings
If the Yukawa couplings in the prime sector y u i , y d i are O(1), the instanton-induced mass of χ is larger than diag . The mass of the bound state η χ , phys will be then dominated by the χ mass instead of the confinement scale. It is thus relevant to consider the instanton effects on the constituent quarks (i.e. χ ), as opposed to the discussion in the previous section in which the instantons contributed directly to the mass of the bound state η χ . As a consequence, the relevant diagram (Fig. 4) generates an effective mass m χ for the χ quark induced by SSI instantons: which can be computed using the dilute instanton gas approximation:  The slight kinks in the η ψ curves mark the χ threshold. The η χ axion mass, dominated by the instanton effects on the constituent quarks, is shown together with the mass for the dynamical η ψ axion sourced by diag instantons. The α CUT values corresponding to solid lines are allowed, while those for dashed lines correspond to excluded values, for which SU (3) diag would confine below 2.9 TeV, ruled out by searches for scalar octets at LHC [45] (see Sect. 4.1) For the benchmark values used for CUT this indicates a mass of tens of TeV or above. As expected, the mass of the χ thus lies well above the confinement scale d , and therefore the mass of the η χ will be dominated by the mass of its constituent fermion: m η χ ∼ 2m χ The two axion masses are represented in the Fig. 5, where the mass of the η ψ is taken from Eq. (53) and that of η χ from Eq. (55). Both composite axions have masses typically above the TeV regime and are out of present collider reach.

Solution to the strong CP problem
Upon SU (3) diag confinement, only one massless fermion remains in the spectrum, ψ. Therefore, the potential contains only one pseudoscalar meson that inherits a shift symmetry: η ψ . This does not invalidate the solution to the strong CP problem, since the other phaseθ in Eq. (49) is reabsorbed when the χ is integrated out of the spectrum. This is a consequence of the mass of the χ being generated by instanton effects: the phase of this mass term will be exactly that of the topological term, and it will be completely removed after integrating out the fermion. After taking all instanton effects into account, the low energy potential for the dynamical axion now reads: and the minimum is CP-conserving:

Low energy spectrum and observable effects
For most of the parameter space, the two dynamical axions η ψ, phys and η χ, phys are typically heavier than the TeV scale and thus very difficult to observe at LHC. Otherwise, under a few TeV the spectrum is the same as that itemized in Sect. 3.1.4, except that there are no color triplet "axi-pions" because of the χ absence at the relevant energies, see Eq. (29). Therefore only the color octet "axi-pion" can be searched for at the LHC.

Model II: Addition of a second scalar.
This solution to the θ problem is an alternative to extending the spectrum by a massless fermion, discussed in the previous subsection. In this second model no extra fermion is added to the SU (6) × SU (3 ) Lagrangian, while a second field will be considered instead. The spectrum is that in Table 3 albeit with the scalar line duplicated, → { 1 , 2 }. This simple extension allows the implementation of a PQ symmetry which reabsorbs the θ contribution to the strong CP problem. The corresponding PQ symmetry is automatic if the terms in Eq. (11) are omitted and Eq. (7) is replaced by This Lagrangian is invariant under two independent abelian global symmetries; one of them is anomalous with respect to SU (3 ) and corresponds to the PQ charge assignment 18 P Q{ 1 , 2 , q L , u c L , d c L } = {+1, −1, +1, +1, +1} . (59) The vevs of 1 and 2 generalize the CUT spontaneous breaking in Eq. (8) and at the same time break spontaneously the PQ symmetry; therefore, this PQ scale coincides with the CUT scale. This distinguishes model 2 from Model I, as in the latter the PQ scale coincided with diag . A pGB -an elementary axion-is generated at this stage. The corresponding PQ conserved current is given by The i fields are parameterized as where v 1 and v 2 denote respectively the 1 and 2 vevs, both of which we take to be real for simplicity. Decoupling the heavy radial modes, the PQ current reads where the elementary axion field a(x) corresponds to the GB combination with This classically exact PQ symmetry is broken at the quantum level by the SU (3 ) anomaly, which at lower energies translates into an anomalous current for the SU (3) diag gauge theory.
where N and N diag are the group factors, and where T a P Q corresponds to the PQ generator and t b = λ b 2 to the Gell-Mann matrices for the SU (3) generators. The anomalous term modifies the classical equations of motion of the axion, and gives rise to an effective Lagrangian, The impact of the SSI of the spontaneously broken theory will again add further contributions, inducing a putatively high mass for the elementary axion as discussed further below.

Running of the coupling constants
The matter content allows SU (3) diag to confine at higher scales than the QCD group SU (3) c as in Model I. The separation of both scales is made even sharper in the Model II because α diag runs faster. In Model II, only one massless fermion charged under SU (3) diag is ever present under the CUT scale (compare Tables 5 and 6). We have estimated both the one and two-loop running, as illustrated in Fig. 6 for diag = 4 TeV.

Confinement of SU (3) diag and pseudoscalar anomalous couplings to the confining interactions.
At the scale diag the QCD coupling constant is small, and the SU (3) diag spectrum with only one massless fermion in Table 6 has an approximate classical global symmetry The gauge QCD group SU (3) c is again a subgroup of U (3) V which remains unbroken. The octet of pGBs colored under QCD will acquire large masses due to gluon loops, which at one-loop is given by The QCD singlet 1 c , denoted η ψ , is a dynamical axion. Note that it has the same quark composition as the η ψ meson in Model I. The η ψ couples to both the SU (3) diag and SU (3) c anomalies, resulting in a low-energy effective Lagrangian for this axion given by In summary, this solution to the strong CP problem is a hybrid one with two axions: a heavy dynamical axion η ψ with mass of order diag stemming from a PQ symmetry which reabsorbs the original θ SU (6) (and thus θ QCD ) parameter as in the previous section, and a second elementary axion a resulting from solving the external SU (3 ) sector à la PQWW [12,13,46]. Up to now, only two sources of masses have been identified for the ensemble of three pseudoscalars with anomalous couplings (η QCD , η ψ and a). We analyze next the SSI of this model which provide a large source of axion mass for the elementary axion a.

Small-size instantons of the spontaneously broken CUT
The analysis of SSI for Model II under discussion is simpler than that for Model I developed in the previous subsection.
No massless fermions charged under SU (3 ) are present in Model II (in contrast with Model I). Furthermore, PQ symmetry forbids here the y i Yukawa couplings which gave the dominant contribution in Model I. In consequence, the terms proportional to κ i and mediated by the 1 and 2 scalars and the Higgs field will dominate SSI . This is illustrated by the instanton "flower" in Fig. 7. It results in SSI given by which can be written as where For the benchmark values used for CUT , this implies an elementary axion mass in the range 10 −5 − 1 GeV, see Fig. 8 for illustration. In summary, the dynamical axion has a mass of order diag and thus of a few TeV or above, while the elementary axion is light although typically heavier than the usual invisible axion.

Solution to the strong CP problem
The minimum of the axion potential can be easily shown to remain CP-conserving after including all contributions to the axion masses. Indeed, the θ i dependence of the potential can be again read off of Eq. (13), In this model both the usual QCD η and the axion a phys remain as light eigenstates, while the other eigenstate η ψ, phys will have a mass generically above the TeV scale. Apart from the axion, the lowest set of exotic states is an octet of exotic "pions" whose masses are TeV, see Eq. (70). A similar colored octet is present in Model I discussed in Sect. 3.1, although Model I contains an additional set of color-triplet pseudoscalars in the case of small y i Yukawa couplings. In contrast, no color-triplet is ever expected here as the exotic classical flavor symmetry is This Model II with an additional scalar may be less appealing than than Model I with an extra massless fermions for two reasons: (a) its axion sector contributes directly to the EW hierarchy problem, as its elementary axion results from a scalar potential which a priori communicates with the Higgs potential; (b) it is a hybrid model with both one elementary and one dynamical axion, while Model I is more aligned with the spirit of solving fully the strong CP problem via massless fermions.

Phenomenological and cosmological limits on the lightest exotic states
A common feature of the ultraviolet complete models constructed above is that the generic spectrum under the EW scale includes, in addition to the SM spectrum, sterile fermions, in contrast with usual axion models. An axion may also be present in this range depending on the model parameters.

Collider observable signals
A set of observable exotic states are expected to be the exotic SU (3) c -colored "pions" whose masses may lie under the TeV scale. These resulted from the chiral symmetry breaking of the confining group SU (3) diag . All models exhibit as a common characteristic QCD color-octet meson bound states made out of their massless fermions, shown in Table 5 and Eq. (17) for Model I and Table 6 and Eq. (70) for Model II.
In addition, QCD color-triplet meson bound states may be observable for Model I with small Yukawa couplings in the primed sector.
In this color-unified axion solution, the exotic fundamental fermions have no SM SU (2)×U (1) charges. The heavy pions will be produced in colliders only via QCD interactions, e.g. gluon-gluon couplings, through which they also presumably decay before they can hadronize to make color neutral states. As they are colored, they do not mix with ordinary pions or other visible matter.

Color-octet pions from SU (3) diag confinement
The lightest scalar octets in Eqs. (17) and (70) are denoted by π d ≡ 8 c . Their effective coupling to QCD gluons can be written as where D μ denotes the SU (3) c covariant derivative and d abc is the corresponding symmetric group structure constant. The second term in Eq. (83) results at one-loop from the triangle diagram with the fermions χ running in the loop. The coloroctet exotic pions can thus be either produced in pairs through the gluon-gluon-π d -π d coupling in the kinetic term, or singly produced through their anomalous coupling to gluons. The kinetic term dominates the scalar octet production channels, while the second term allows the π d decay, yielding a dijet final state. Experimental limits on scalar octet pair production via the gluonic interactions in the kinetic term can be inferred by recasting searches of sgluons. A recent search of sgluon pair production by ATLAS using 36.7 fb −1 of √ s = 13 TeV data [45], whose prediction was obtained from the NLO computation in Ref. [47] for √ s = 8 TeV (rescaled to 13 TeV according to Ref. [48]), sets a bound on the octet scalar exotic pions given by From Eqs. (17) and (70), this translates into diag 2.9 TeV .
For Model I and in the particular case of small Yukawa couplings in the primed sector, an alternative bound may be inferred from the limits on color-triplet scalars, which can be produced via their color interactions. In the absence of couplings which mediate their decay (as the second type of coupling in Eq. (83) is not possible for scalar triplets), they will bind with SM quarks to form stable hadrons. This search is expected to result in a sensitivity similar to the one above [49]. It is very interesting to pursue the experimental search for colored pseudoscalars and stable exotic hadrons. Their detection would be a powerful indication of the dynamical solution to the strong CP problem proposed here.

Dynamical axion and exotic fermions
The dynamical axion denoted above by η ψ, phys with instanton-induced mass of order diag , Eqs. (53) and (82), can a priori be either pair-produced through the kinetic coupling or singly produced through the dimension five anomalous operator. It would decay dominantly to two back-to-back jets and can be searched for in dijet resonance searches. Its production, however, may be suppressed by its high mass. For instance, for diag 2.9 TeV, m η ψ, phys ∼ 90 TeV, which is beyond the reach of present collider searches. The other mesons and baryons resulting from the SU (3) diag confinement have masses in the TeV range and above; they would lead to collider signatures similar to those for the exotic pseudo-goldstone bosons.

The axion coupled to SSI
In Model I, the second axion is also dynamical and it was denoted by η χ, phys . Its mass can span a large range, depending on the strength of the Yukawa couplings in the prime sector: -For y i couplings of O(1), this second axion becomes typically heavier than the TeV scale and is thus out of LHC reach. -For small y i couplings (O(0.2) have been illustrated), an observable axion with mass in the GeV-tens of TeV regime is realistic. Given the low PQ scale ( f d , close to a TeV) lower masses are excluded.
The second axion in Model II, denoted a and of elementary nature, is expected to be light although heavier than the invisible axion, with a mass in the 10 eV -10 GeV range for O(1) Yukawa couplings in the prime-fermion sector. By lowering the values of the latter couplings, it can even become as light as the usual invisible axion and with similar phenomenology, as its associated axion scale is the color-unified scale, f a ∼ CUT .

Cosmological and gravitational aspects
We briefly discuss next the cosmological aspects of the models constructed above, as well as the putative instability threat from gravitational non-perturbative effects.

Stable particles and cosmological structures
Stable particles with masses higher than about 10 5 GeV may lead to cosmological problems, dominating the mass density and overclosing the universe [50]. This is often a problem in previous models of composite axions because exotic stable baryons bound by the extra confining force [11,51] are expected.
However, as pointed out in Refs. [51,52], if axicolor can be unified with a SM gauge group, then the unified forces could mediate the decay of axihadrons into lighter states, and the model would be cosmologically safe. The color unification of this work automatically employs this mechanism. The heavy exotic hadrons decay to the sterile fermions ψ ν , which are part of the CUT massless multiplet in Eq. (4) and Table 2. It remains to be determined whether the lifetime for these CUT-induced decay channels is too large to avoid problems from stable exotic hadrons. If the decays are made to be fast enough, the resulting sterile fermions may induce in turn other cosmological problems. This analysis is left to a future work.
A similar concern pertains to the domain walls which may arise due to the spontaneous breaking of the discrete shift symmetry in the instanton potential, Eqs. (51) and (80). The walls need to disappear before they dominate the matter density of the universe, or else other mechanisms must be applied to solve the domain wall problem [53][54][55][56][57].
In any case, if the universe went through an inflation phase any relic previously present will be wiped out. If the reheating temperature is lower than the PQ scales, meaning lower than diag here, neither the heavy stable particles nor any putative domain walls are produced again after inflation, and the problems mentioned above are avoided altogether. We defer to a future work the in-depth study of the cosmological aspects of Model I and Model II, with either high-or low-scale inflation.

The question of gravitational quantum effects
Gravitational quantum corrections have been suggested to be relevant and dangerous for axion models in which the PQ scale is not far from the Planck scale. In Model I, both PQ scales correspond to diag , which is much lower than the Planck scale, and no instability resulting from gravitational quantum effects is at stake. In Model II instead, while the dynamical PQ scale is analogously low, the second PQ symmetry is realized at the CUT scale and gravitational quantum effects could be relevant.
It has often been argued that all global symmetries may be violated by non-perturbative quantum gravitational effects, see for instance Refs. [54,[58][59][60][61]. For instance, a black hole can eat global charges and subsequently evaporate. Similar effects may exist with virtual black holes. Another indication that gravity might not respect global symmetries comes from wormhole physics [62][63][64][65]. The natural scale of violation in this case is the wormhole scale, usually thought to be very near (within an order of magnitude or so) the Planck mass M Planck .
For axion models with high PQ scales, such as the typical scale of invisible axion models f a ∼ 10 9 − 10 12 GeV, it has been argued that those non-perturbative quantum gravitational effects could lead to extreme fine-tunings. The authors of Ref. [54,[58][59][60] concentrated on the simplest (and the most dangerous) hypothetical dimension five effective operator where g 5 is a dimensionless coefficient and would be a field whose vev breaks the PQ invariance. This term threatens the standard invisible axion solutions to the strong CP problem, as it would change the shape of the effective potential. The minimum moves unacceptably away from a CPconserving solution unless the coefficient is strongly finetuned, for instance g 5 < 10 −54 for f a ∼ 10 12 GeV. These potentially dangerous terms can be avoided if the global PQ symmetry arises accidentally as a consequence of other gauge groups [66,67]. Nevertheless, the idea that gravity breaks all global symmetries is indeed an assumption -and sometimes an incorrect one -at least at the Lagrangian level. 19 Furthermore, very recently the impact of non-perturbative effects has been clarified and quantified in Ref. [69]. The effects may be extremely suppressed by an exponential dependence on the gravitational instanton action, and be harmless even with high axion scales. The demonstration relied on assuming that the spontaneous breaking of the PQ symmetry is implemented through the vev of a scalar field, and thus it directly applies to our Model II. 20 In summary, Model I is safe from instabilities induced by gravitational quantum corrections and plausibly this also applies to Model II.

Conclusions
Color unification with massless quarks has been proposed and developed here for the first time. As a simple implementation of the idea, the SM color group has been embedded in SU (6), which is spontaneously broken to QCD and a second confining and unbroken gauge group. An exactly massless SU (6) fermion multiplet solves the strong CP problem. We have fully developed two ultraviolet completions of the mechanism.
In order to implement this idea successfully, it is necessary to give satisfactorily high masses to the SU (6) partners of the SM quarks to achieve a separation between the QCD scale and that of the second confining group. For this purpose, an auxiliary SU (3 ) gauge group is introduced under which the aforementioned massless fermion is a singlet. SU (6) × SU (3 ) −→ SU (3) c × SU (3) diag is a simple and realistic option. Both final groups remain unbroken and confine at two different scales, QCD and diag , with diag ∼ O(# TeV) QCD . The scale diag then gives the order of magnitude of the mass of the dynamical composite axion inherent to the color-unified mechanism. Further-more, massless (or almost massless) sterile fermions are a low-energy trademark remnant of the massless multiplet that solves the SM strong CP problem.
In order to avoid the SU (3 ) sector sourcing back an extra contribution to the strong CP problem, a minimal extension of its matter sector suffices. Two examples of ultraviolet complete models have been explored in this work: in Model I an extra SU (3 ) massless fermion is added, while Model II includes instead a second scalar with the same quantum numbers as the color-unification breaking scalar. From the point of view of the strong CP problem, those two models are very different. Model I features a second dynamical axion with a second PQ scale which is also of order diag and thus low. In Model II, this second PQ scale coincides instead with the much larger color-unification scale, and the associated axion is elementary. We computed the two-loop running of all coupling constants involved, showing that the desired separation of all relevant scales is achieved naturally: a color-unification scale much larger than the two confining ones, diag and QCD , and the subsequent separation of the last two. This separation of scales is robust and stable over a wide range of parameter values.
We have found that regardless of the details of the ultraviolet implementation, generically there are the three sources of anomalous currents: the instantons of the confining SU (3) c , the instantons of the confining SU (3) diag , and finally the small-size instantons of the spontaneously broken colorunified theory. There are thus three diverse sources of mass for the three pseudoscalars in the theory which couple to anomalous currents: the QCD η , the dynamical axion inherent to color unification, and the second axion (either dynamical or elementary) associated to the solution of the θ problem. These three bosons then acquire masses of order QCD , diag and SSI , respectively, and no standard invisible axion (coupling anomalously only to QCD instantons) is left in the low-energy spectrum. This is generically a very interesting mechanism from the point of view of solving the strong CP problem with heavy axions and scales around a TeV. The mechanism allows a wide extension beyond the invisible axion range of axion parameter space which solves the strong CP-problem. Typical axion masses are in the TeV regime and above, although strictly speaking one of the axions can become as light as the usual invisible axion.
With axion scales around the TeV, observable signals at colliders are expected, as well as other rich phenomenology. Generically, the lightest exotic bound states are colored pseudoscalars (QCD octets and in some cases also triplets in Model I, and only octets in Model II). We have recast the results from present experimental searches for heavy colored mesons to infer a 2.9 TeV bound on the confinement scale of the second confining group, SU (3) diag , which is directly related in Model I to the axion scale.
Overall, Model I may be preferred as: (i) it is exclusively based on solving the strong CP problem dynamically via massless quarks; (ii) from the point of view of naturalness it does not require any fine-tuning to ensure the hierarchy between the PQ and the electroweak scales, as no PQ field is involved in the scalar potential. In Model II instead, one PQ field participates in the color-unification scalar potential, and furthermore this model is a hybrid dynamical-elementary axion solution to the strong CP problem.
Model I is also unquestionably safe from the point of view of stability of the axion solution with respect to nonperturbative effects of quantum gravity, as its PQ scales are diag ∼ TeV. Furthermore, recent advances suggest that the quantum gravity threat should not be considered a risk even for Model II. The other issue of the cosmological impact of (quasi) stable heavy exotic hadrons and of the (almost) massless sterile fermion remnants can be simply avoided by introducing an inflation scale and reheating temperature lower than diag . This last subject deserves future detailed attention in particular in view of the dark matter puzzle.