Axion Quality from Superconformal Dynamics

We discuss a possibility that a superconformal dynamics induces the emergence of a global $U(1)_{\rm PQ}$ symmetry to solve the strong CP problem through the axion. Fields spontaneously breaking the $U(1)_{\rm PQ}$ symmetry couple to new quarks charged under the ordinary color $SU(3)_C$ and a new $SU(N)$ gauge group. The theory flows into an IR fixed point where the $U(1)_{\rm PQ}$ breaking fields hold a large anomalous dimension leading to the suppression of $U(1)_{\rm PQ}$-violating higher dimensional operators. The spontaneous breaking of the $U(1)_{\rm PQ}$ makes the new quarks massive. The $U(1)_{\rm PQ}$ symmetry is anomalous under the $SU(3)_C$ but not under the $SU(N)$ so that the axion couples to only the color $SU(3)_C$ and the usual axion potential is generated. We also comment on a model that the $U(1)_{\rm PQ}$ breaking fields are realized as meson superfields in a new supersymmetric QCD.

Introduction.-The strong CP problem is an intriguing puzzle to motivate physics beyond the Standard Model (SM). The current upper bound on the neutron electric dipole moment constrains the absolute value of the QCD vacuum angleθ to be smaller than 10 −11 [1,2]. Unlike other naturalness problems in the SM, some shifts ofθ would not provide a visible change in our world. The most common explanation for the strong CP problem is the introduction of a pseudo-Nambu-Goldstone boson, called axion a [3,4], associated with spontaneous breaking of a global U (1) Peccei-Quinn (U (1) PQ ) symmetry [5] (for reviews, see e.g. refs. [6,7]). Non-perturbative QCD effects break the U (1) PQ explicitly and generate a potential of the axion whose minimum setsθ to zero. Astrophysical observations provide a lower limit on the U (1) PQ breaking scale, f a 10 8 GeV [8].
In this letter, we explore an alternative approach to the axion quality problem that a superconformal dynamics induces the emergence of the U (1) PQ symmetry. Our model begins with the existence of a discrete Z N with N ∼ 5 which ensures that the model respects the U (1) PQ symmetry at the renormalizable level. We introduce a SU (N ) supersymmetric gauge theory with (anti-)fundamental quarks, some of which are also charged un-der the ordinary color SU (3) C . The Z N symmetry is anomaly-free under the SU (3) C as well as the SU (N ). All the new quarks couple to fields responsible for the spontaneous U (1) PQ breaking. The theory flows into an IR fixed point where the U (1) PQ breaking fields hold a large anomalous dimension. Then, even if there exist higher dimensional operators dangerously violating the U (1) PQ at the Planck scale, those operators are significantly suppressed at low-energies. The similar mechanism has been discussed in the context of the Nelson-Strassler model to realize quark and lepton mass hierarchies [44] (for a more recent development using the amaximization technique [45], see refs. [46,47]). According to the AdS/CFT correspondence [48], the approach is similar to that of the warped extra dimension model discussed in ref. [34]. However, to the best of our knowledge, our model is the first 4D calculable realization to utilize a conformal dynamics to suppress U (1) PQ -violating higher dimensional operators. The spontaneous breaking of the U (1) PQ makes all the new quarks massive. The new quarks leading to a large anomalous dimension of the U (1) PQ breaking fields also play the role of the so-called KSVZ quarks [49,50]. Since the U (1) PQ symmetry is anomalous under the SU (3) C but not under the SU (N ), the axion couples to only the color SU (3) C and the usual axion potential is generated. The SU (N ) finally confines and predicts the existence of SU (N ) glueballs.
While the U (1) PQ breaking fields are introduced as elementary fields in the main part of the present work, we will also comment on a possibility that they are realized as meson superfields in a new supersymmetric QCD (SQCD). Interestingly, in the magnetic picture of the theory [51,52], the coupling of the U (1) PQ breaking fields to dual quarks is automatic.
The model.-Let us consider a supersymmetric SU (N ) gauge theory with N f vector-like pairs of chiral superfields in the (anti-)fundamental representation, Q I ,Q I (I = 1, · · · , N f ). Here, N f is assumed to be even.
The charge assignments under the SU (N ) gauge group, the U (1)PQ (and ZN ) and the anomaly-free U (1)R which determines anomalous dimensions of the fields. Here, m = 1, · · · , N f /2 and k = N f /2+1, · · · , N f where N f is even.
We focus on 3 2 N < N f < 3N where the theory is in conformal window [52]. To implement the QCD axion, we introduce two SU (N ) singlet chiral superfields Φ,Φ charged under the U (1) PQ symmetry. They are coupled to the new SU (N ) quarks in the superpotential, where λ,λ denote dimensionless couplings, m runs from 1 to N f /2 and k runs from N f /2 + 1 to N f . These terms explicitly break the original Then, an anomaly-free discrete symmetry Z N ⊂ U (1) PQ is realized, which leads to the U (1) PQ symmetry at the renormalizable level. Below, we will discuss Planck-scale suppressed U (1) PQ -violating operators, but those operators must respect the Z N symmetry. The fields Φ,Φ obtain a non-zero vacuum expectation value (VEV) via the superpotential, which breaks the U (1) PQ symmetry spontaneously. Here, X is a singlet chiral superfield, κ is a dimensionless parameter and f is a constant with a mass dimension. 2 The gauge theory is in conformal window and believed to have a non-trivial IR fixed point. Here, let us assume the SU (N ) gauge coupling g, λ andλ approach values at the fixed point and the theory is in the conformal regime between the energy scales Λ and M c (Λ > M c ). We will demonstrate the existence of the IR fixed point later. In this regime, the conformal dynamics generates a large anomalous dimension of Φ,Φ through the superpotential terms of Eq. (1). The wave function renormalization factor of Φ (andΦ) at IR is given by where γ Φ = 6 N N f − 2 is the anomalous dimension of Φ which is exactly determined in terms of the anomaly-free U (1) R charges summarized in Tab. I. We now canonically normalize Φ as whose hatˆdenotes a field in the canonical normalization. Then, the superpotential (2) is rewritten in terms of the normalized fields, is determined by f which also gives the conformal breaking, M c ∼ f . The wave function renormalization factor of Eq. (3) will play a key role in suppressing U (1) PQ -violating higher dimensional operators as we will see below.
Once the U (1) PQ breaking fields Φ,Φ obtain the VEV, all the new quarks Q I ,Q I become massive, and then the axion-gluon coupling is generated in the effective Lagrangian after the integration of the new quarks, where a denotes the axion, G is the field strength of the gluon,G is its dual, g c is the QCD gauge coupling constant and F a /N = √ 2f /N is the axion decay constant. The same axion-gluon coupling is obtained in the KSVZ axion model [49,50] with N flavors of SU (3) C vector-like quarks. Since the U (1) PQ symmetry is not anomalous under the SU (N ), the terms in Eq. (1) do not lead to the axion-SU (N ) gluon coupling even after the integration of the quarks. The axion potential is obtained via the non-perturbative QCD effect, where m π and f π are the pion mass and the decay constant respectively and m 2 π f 2 π = (0.1 GeV) 4 . Then, the strong CP problem is solved in the ordinary way. After the decoupling of Q I ,Q I , the model becomes a SU (N ) pure Yang Mills theory. Because of a large gauge coupling of the SU (N ) at the fixed point, the theory confines just below the conformal breaking scale M c and predicts heavy SU (N ) glueballs and their superpartners.
Axion quality.-To address the axion quality problem, explicit U (1) PQ breaking terms must be highly suppressed compared to the axion potential generated by the non-perturbative QCD effect (7). The most dangerous Planck-scale suppressed operator respecting the Z N symmetry is the superpotential term such as which leads to the scalar potential in supergravity with e.g. the constant term Here, m 3/2 is the gravitino mass, κP Q is a model dependent coefficient, andΦ denotes the scalar component which is the same notation as the chiral superfield for notational simplicity. The U (1) PQ -violating axion potential is then obtained as where ϕ denotes a CP phase and γ Φ = 6 N N f − 2 has been used. We now define the axion quality factor Q by Assuming ϕ = O(1), the experimental upper bound on theθ parameter [1,2] requires Q 10 −10 to secure the axion quality. Fig. 1 shows the contours of Q calculated from the potential (10) in the m 3/2 − F a /N plane. Here, we take N f = 2 N , M c = F a , κP Q = 1, and Λ = 0.1 M Pl . The solid and dashed lines denote the quality factor Q = 10 −10 , 10 −8 , respectively. The axion decay constant F a /N is constrained from the supernova 1981A observation, F a /N 10 8 GeV [8]. We can see from the figure that there is a parameter space to solve the axion quality problem for N ≥ 5. While the case of N = 4 is not shown in the figure, Q = 10 −5 is obtained for F a /N ∼ 10 8 GeV and m 3/2 ∼ 1 eV.
Other potentially dangerous U (1) PQ -violating operators are with k = 0, · · · , N − 1. While these operators will not lead to the axion potential by themselves, we must be careful because they are enhanced at low-energies due to the negative anomalous dimension of QQ. However, for e.g. N f = 2N , Φ (Φ) and QQ have the same scaling dimension 3/2, and then Eq. (12) can be rewritten as which is suppressed compared to Eq. (8) for Φ ≈ M c and Λ < M Pl . We also note that U (1) PQ -violating operators in the Kähler potential are negligible compared to those in the superpotential.
The IR fixed point.-Let us now discuss the existence of the IR fixed point for the SU (N ) gauge coupling g and λ,λ in the superpotential (1). We first ignore the effect of the SU (3) C gauge coupling and solve the renormalization group equations (RGEs) for g, λ andλ, where t = ln(µ/Λ 0 ) with µ being the RG scale, C A = N and b 0 = 3N − N f . Here, we use the exact NSVZ β function [53][54][55] for the RGE of the gauge coupling, while the RGEs of λ andλ are shown at the one-loop level. The anomalous dimensions are given by with C 2 = N 2 −1 2N . We also calculate the RGEs for λ andλ at the two-loop level whose expressions are summarized in appendix. Fig. 2 shows the RG flows of g and λ from a scale Λ 0 to µ = 10 −9 Λ 0 for different initial values as a demonstration. We take N = 5, N f = 10 andλ = 2 at Λ 0 . Blue and red dots correspond to the cases using the one and two-loop RGEs for λ, respectively. The figure illustrates both couplings flow into a non-trivial IR fixed point. The blue circle around the center denotes the values of g and λ obtained by comparing the anomalous dimensions at one-loop (15) to those determined by the U (1) R charges in Tab. I, where γ Φ = 6 N N f − 2. The anomalous dimensions up to the two-loop order (24) are used to find the values of the couplings at the red circle. We also plot the anomalous dimension of Φ at two-loop γ 2 Φ in the left panel of Fig. 3  (black solid). We take N = 5, N f = 10 and g = λ =λ = 2 at the initial scale Λ 0 . The figure indicates that γ 2 Φ converges to γ Φ = 6 N N f − 2 = 1. Therefore, the theory is expected to enter the conformal regime in the IR region as we have assumed in the above discussion.
The IR fixed point can be disturbed by the SU (3) C gauge coupling. To discuss this effect, we first decompose the superpotential term in Eq. (1) as where Q a ,Q a (a = 1, 2, 3) denote the fundamental and anti-fundamental representations of the SU (3) C gauge group, Q α ,Q α (α = 4, · · · , N f /2) are the quarks that are not charged under the SU (3) C and λ 1,2 are dimensionless couplings. The anomalous dimensions including the SU (3) C effect at the two-loop level are summarized in appendix. We use the one-loop RGE for the SU (3) C gauge coupling, which is solved as 4π where we take b 3 = 3 − N for µ > M c by assuming all the new quarks have masses around M c . Here, the factor 3 is from the MSSM particles and the factor −N is from the Q a ,Q a quarks. For N = 5, the SU (3) C gauge coupling becomes asymptotic non-free. In this case, we obtain g c ≈ 1 around µ = 10 17 GeV for the spectrum of the MSSM particles at about 10 TeV and 4π/g 2 c | µ=Mc ≈ 20 at M c = 10 8 GeV. We numerically solve the twoloop RGEs from a scale Λ 0 to µ = 10 −9 Λ 0 . The left panel of Fig. 3 shows the flow of γ 2 Φ for g c = 1, 2 at Λ 0 denoted by the red dashed and green dotted lines, respectively. The initial values of the couplings at Λ 0 are g = λ 1 = λ 2 =λ = 2. We also plot the flow of the wave function renormalization factor Z Φ for g c = 0, 1, 2, 3 and g = λ 1 = λ 2 =λ = 2 at Λ 0 in the right panel of Fig. 3. From the figures, we can confirm that γ 2 Φ converges into the one without the SU (3) C effect and the smallness of Z Φ enables to solve the axion quality problem.
A model with the dual picture.-So far, we have discussed the model where the U (1) PQ breaking fields are introduced as elementary fields, but here let us comment on a possibility that they are realized as meson superfields in a new SQCD. Consider a SU (N f − N ) SQCD with N f vector-like pairs of quarks whose dual magnetic picture is given by a SU (N ) SQCD with the same number of flavors D i ,D i (i = 1, · · · , N f ) [51]. In the magnetic theory, there also exist meson chiral superfields M i j which are coupled to the dual quarks through the superpotential, where y is a dimensionless coupling. For 3 2 N ≤ N f ≤ 3N , where a, b (= 1, 2, 3) denote the color SU (3) C indices, α, β = 4, 5, 6 andī,j = 7, · · · , N f . The U (1) PQ charges are, for example, assigned as shown in Tab. II. With these assignments, the U (1) PQ symmetry is not anomalous under the SU (N ) but is anomalous under the SU (3) C .
With the decomposition of Eq. (21), we can see that the superpotential (20) contains the terms similar to those introduced in Eq. (1), Here, we have defined M 1 ≡ 1 3 M a 1a and M 2 ≡ 1 3 M α 2α . Note that M 1,2 are color singlet but U (1) PQ charged. Once they obtain non-zero VEVs, we get the axiongluon coupling (6). As before, the U (1) PQ symmetry at the renormalizable level is ensured by an anomaly-free Z N ⊂ U (1) PQ . Explicit U (1) PQ -violating higher dimensional operators are suppressed due to large anomalous dimensions of M 1,2 .
Several comments are in order. The IR fixed point can be disturbed by the SU (3) C gauge interaction. In order to keep the electric/magnetic duality reliable, the values of the couplings in both electric and magnetic pictures at the fixed point must be much larger than the QCD gauge coupling, which requires the theory to be near the middle of conformal window, N f ≈ 2 N . Extra meson and quark chiral superfields must get masses appropriately. In particular, SU (3) C -charged mesons must be stabilized at the origin to avoid the color breaking. If M 3 ≡ 1 N f −6 Mī 3ī obtains a non-zero VEV, all the quarks become massive. Below the scales of M 1,2,3 VEVs, the model becomes a confining SU (N ) pure Yang Mills theory. Further explorations of this model are left to a future study.
Conclusions and discussions.-We have considered a possibility that a superconformal dynamics helps to solve the strong CP problem through the axion with a sufficient quality. The U (1) PQ breaking fields are coupled to the new quarks charged under the SU (3) C and the new SU (N ). The theory flows into a non-trivial IR fixed point where the U (1) PQ breaking fields hold a large anomalous dimension leading to a strong suppression of explicit U (1) PQ breaking operators. The U (1) PQ is anomalous under the SU (3) C but not under the SU (N ) so that the usual axion potential is generated by non-perturbative QCD effects.
The model respects the anomaly-free Z N ⊂ U (1) PQ , which realizes the U (1) PQ symmetry at the renormalizable level. If the U (1) PQ is spontaneously broken after the end of inflation, cosmic strings are formed at a temperature close to the U (1) PQ breaking scale (see e.g. ref. [56] for a review on axion cosmology). Below around the QCD temperature, domain walls attached to the cosmic strings are formed. They are stable due to the Z N symmetry and cause a cosmological problem. In order to avoid this, the U (1) PQ symmetry must be broken before the end of inflation. In this case, the axion isocurvature perturbation is produced, which leads to a constraint on the Hubble scale of inflation, H inf 10 7 GeV. Cosmological aspects might be an interesting future direction.
We may be able to use the same superconformal dynamics to realize the quark and lepton mass hierarchies in the same way as the Nelson-Strassler model [44]. Such a possibility has been recently discussed in the 5D context [35]. One extra benefit of this scenario is that flavor-dependent soft scalar masses are automatically suppressed [57,58] (see also ref. [59]).

ACKNOWLEDGEMENTS
We would like to thank Ryosuke Sato for discussions and helpful comments on the manuscript. We are also grateful to Kavli IPMU for their hospitality during the COVID-19 pandemic.