Minimal axion model from flavor

Lorenzo Calibbi, Florian Goertz, Diego Redigolo, Robert Ziegler, and Jure Zupan CAS Key Laboratory of Theoretical Physics and Kavli Institute for Theoretical Physics China (KITPC), Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, People’s Republic of China Theory Division, CERN, CH-1211 Geneva 23, Switzerland Raymond and Beverly Sackler School of Physics and Astronomy, Tel-Aviv University, Tel-Aviv 69978, Israel Department of Particle Physics and Astrophysics, Weizmann Institute of Science, Rehovot 7610001, Israel Institute for Theoretical Particle Physics (TTP), Karlsruhe Institute of Technology, Engesserstrasse 7, D-76128 Karlsruhe, Germany Department of Physics, University of Cincinnati, Cincinnati, Ohio 45221,USA (Received 23 January 2017; published 15 May 2017)


I. INTRODUCTION
Three of the major open questions in particle physics are (i) the strong CP problem-why is the QCD θ angle so small, (ii) what is the origin of Dark Matter (DM), and (iii) the Standard Model (SM) flavor puzzle-why are the masses of fermions so hierarchical.The first problem can be elegantly addressed by the QCD axion: the pseudo-Goldstone boson of an approximate global Uð1Þ symmetry that has a color anomaly [1][2][3].The two main classes of axion models based on this mechanism are usually referred to as the KSVZ [4,5] and the DFSZ [6,7] axion solutions.It is well known [8][9][10] that in most regions of the parameter space the QCD axion serves as a viable DM candidate.Finally, the SM flavor problem can be elegantly resolved by introducing approximate flavor symmetries, which are spontaneously broken at large scales as in the original Froggatt-Nielsen (FN) mechanism [11].
In this paper we propose a unified framework where the approximate symmetry of the QCD axion is identified with the simplest flavor symmetry of the FN mechanism (the setup is the minimal realization of an old idea by Wilczek [12] that axion and flavor physics could be connected).The structure of quark and lepton masses and mixings follows from a spontaneously broken Uð1Þ H flavor symmetry which generically has a QCD anomaly.The resulting Nambu-Goldstone boson, the axiflavon, solves automatically the strong CP problem by dynamically driving the theory to a CP conserving minimum [13].Nonthermal production of the axiflavon from the misalignment mechanism can then reproduce the observed DM relic density, provided that the Uð1Þ H breaking scale is sufficiently large.
This simultaneous solution of flavor, strong CP and DM problem leads to sharp predictions for the properties of the axiflavon that can be tested experimentally. 1 Of particular importance is the axion coupling to photons that is determined by the ratio E=N, i.e., the ratio of the electromagnetic over the QCD anomaly coefficient.This ratio is essentially a free parameter in generic axion models (see Refs. [15,16] for a recent discussion).In the axiflavon setup E=N is directly related to the Uð1Þ H charges of SM fermions and thus to the hierarchies between SM fermion masses.Despite the considerable freedom of choosing these charges in the simplest Uð1Þ H model, we find a surprisingly sharp prediction for E=N centered around 8=3, the prediction of the simplest DFSZ model, This result is a direct consequence of the strong hierarchies in up-and down-type quark masses and only weak hierarchies in the ratio of down-quark to charged lepton masses.A similarly restrictive range for E=N can be found also in a broad class of models with nonminimal flavor symmetries like Uð2Þ (which are more predictive in the fermion sector).The above range for E=N can be translated into a prediction for the ratio of axion-photon coupling g aγγ and axion mass m a , g aγγ m a ∈ ½1.0; 2.2 10 16 GeV For axion masses in the natural range for axion DM, m a ≈ ð10 −3 ÷ 0.1Þ meV, this region will be tested in the near future by the ADMX experiment.
The axiflavon can also be tested by precision flavor experiments looking for the decay K þ → π þ a. Indeed the flavor violating couplings of the axiflavon to quarks are also related to quark masses, but in contrast to E=N are more sensitive to model-dependent Oð1Þ coefficients: where κ sd =N ∼ Oð1Þ.In the natural range of axion DM this decay can be within the reach of the NA62 and ORKA experiments, depending on the model-dependent coefficients.We summarize our results along with the present and expected experimental constraints in Fig. 1 at the end of this paper.

II. SETUP
We assume that the masses of the SM fermions come from the vacuum expectation value (vev) v ¼ 174 GeV of the SM Higgs H, while the hierarchies of the Yukawa couplings are due to a global horizontal symmetry Uð1Þ H .The SM Weyl fermion fields Q i ; U c i ; D c i ; L i ; E c i have positive flavor-dependent charges ½q i ; ½u i ; ½d i ; ½l i ; ½e i , respectively.Here Q i and L i are the quark and lepton electroweak doublets, the remaining fields are SUð2Þ L singlets, and i ¼ 1, 2, 3 is the generation index.For simplicity we assume that the Higgs does not carry a Uð1Þ H charge, so that the flavor hierarchies are explained entirely by the fermion sector.This assumption will be relaxed below.The Uð1Þ H symmetry is spontaneously broken at a very high scale by the vev V Φ of a complex scalar field Φ with Uð1Þ H charge of −1.All other fields in the model, the FN messengers, have masses of OðΛÞ ≳ V Φ ≫ v and can be integrated out.Note that Λ is a scale above Uð1Þ H breaking, implying that fermionic FN messengers are vectorlike under the Uð1Þ H .The Yukawa sector in the resulting effective theory is then given by where a u;d;e ij are complex numbers, assumed to be Oð1Þ.
, gives the SM Yukawa couplings with where ½L i ¼ ½q i ; ½R i ¼ ½u i ; ½d i in the quark sectors, ½L i ¼ ½l i ; ½R i ¼ ½e i in the charged lepton sector and we have defined the small parameter ϵ ≡ V Φ =ð ffiffi ffi 2 p ΛÞ.The hierarchy of masses follows from Uð1Þ H charge assignments, giving Here V is the Cabibbo-Kobayashi-Maskawa (CKM) matrix in the quark sector and the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix in the charged lepton sector.The observed CKM structure is typically obtained for ϵ of the order of the Cabibbo angle, ϵ ∼ 0.23.The exact values of Uð1Þ H charges can be obtained from a fit to fermion masses and mixings, and are subject to the uncertainties in the unknown Oð1Þ numbers a u;d;e ij .As we are going to demonstrate, these uncertainties will only weakly influence the main phenomenological predictions.Note that the pattern of masses and mixings in the neutrino sector can also be explained in this setup, however, this sector of the SM is irrelevant for the prediction of color and electromagnetic Uð1Þ H anomalies.
The field Φ contains two excitations, the CP-even flavon, ϕ, and the CP-odd axiflavon, a, The flavon field ϕ has a mass m ϕ ∼ OðV Φ Þ, and thus is not directly relevant for low energy phenomenology, and can be integrated out.The axiflavon, a, is a Nambu-Goldstone boson.It is massless at the classical level, but receives a nonzero mass from the breaking of Uð1Þ H by the QCD anomaly.Its couplings to SM fermions F i are given by The couplings of the axiflavon to the SM fermions are in general not diagonal in the fermion mass eigenstate basis due to the generation dependency of charges ½q i ; ½u i ; ½d i ; ½l i ; ½e i .This induces flavor changing neutral currents, which are experimentally well constrained and will be discussed in the next section. 2Note that several axion models with flavor-violating couplings to fermions have been proposed in the literature, see e.g.[18][19][20][21][22][23][24][25][26].In the axiflavon setup they are directly related to the SM fermion masses and thus predicted up to Oð1Þ uncertainties.

2
For flavor constraints on a heavy CP-odd flavon and possible collider signatures see Ref. [17].
The axiflavon couplings to gluons and photons are controlled by the color and electromagnetic anomalies, where Gμν ¼ 1 2 ϵ μνρσ G ρσ and we have switched to the standard axion notation introducing f a ¼ V Φ =2N.The two anomaly coefficients, N, E, are completely determined by the Uð1Þ H charges of SM fermions, in the minimal scenario where these are the only states with chiral Uð1Þ H charge assignments (see a more detailed discussion below).Interestingly, these coefficients can be directly related to the determinants of the fermion mass matrices as [27][28][29] det where the quantities α ud ¼ det a u det a d and α de ¼ det a d = det a e contain the Oð1Þ uncertainties, given by the anarchical coefficients in Eq. ( 4).Taking fermion masses at 10 9 GeV from Ref. [30], one finds det m u det m d =v 6 ≈ 5 × 10 −20 and det m d = m e ≈ 0.7, which makes it clear that up to small model-dependent corrections we have E ¼ 8=3N and so are close to the simplest DFSZ axion solution [31].Indeed the phenomenologically relevant ratio E=N is independent of ϵ and given by The most natural values for the coefficients are α ud ¼ α de ¼ 1, in the sense that Yukawa hierarchies are entirely explained by Uð1Þ H charges, giving E=N ≈ 2.7.To estimate the freedom from Oð1Þ uncertainties, we simply take flatly distributed Oð1Þ coefficients in the range ½1=3; 3 with random sign.This results in a 99.9% range given in Eq. (1) or j E N −1.92j∈½0.5;1.1 (allowing for loop induced textures, a ij ∈ ½1=ð4πÞ 2 ; 3, with a u 33 ∈ ½1=3; 3, still gives a relatively narrow range j E N − 1.92j ∈ ½0.3; 1.7).This should be compared with the usual KSVZ/DFSZ axion window j E N − 1.92j ∈ ½0.07; 7 [32].Note that the restricted range is due to the suppression of the second term in Eq. ( 14) since the denominator is dominated by log det m u det m d =v 6 ≈ −44, while the first term in the numerator is log det m d = det m e ≈ −0.36.Following Ref. [33], we therefore obtain a quite sharp prediction for the axion-photon coupling, 1  4 g aγγ aF F, as while the axion mass induced by the QCD anomaly is given by [33] m a ¼ 5.7 μeV It is remarkable that the prediction for E=N in Eq. ( 1) is largely insensitive on the details of the underlying flavor model.We therefore briefly review the underlying assumptions that lead to the above results and discuss their relevance and generality.For simplicity we assumed in the above discussion positive fermion charges.This assumption can be relaxed to the extent that just the sums of charges in each Yukawa entry are positive.This assumption can be even further relaxed, if we embed the setup into a supersymmetric model in order to address also the hierarchy problem.The negative sums of charges then give holomorphic zeros in the mass textures, while our result ( 14) remains unchanged.Our second assumption was that only the fermion fields and the flavon carry the Uð1Þ H charges.This assumption can be easily dropped since a possible Uð1Þ H charge for the Higgs, [h], would simply drop out of Eq. ( 1), as it would enter as det m u → det m u ϵ 3½h and det m d;e → det m d;e ϵ −3½h .Finally we have assumed that only light fermions contribute to the QCD and electromagnetic anomalies, i.e., that all the other fields in the model are either bosons or vectorlike fermions under Uð1Þ H .This is a natural feature of the FN messengers needed to UV complete the effective setup in Eq. ( 4), see also the explicit UV completions in Refs.[34,35].We also note that the same prediction for E=N holds in any flavor model where a global, anomalous Uð1Þ factor determines exclusively the determinant of the SM Yukawa matrices.For example in Uð2Þ flavor models [36][37][38][39], where the three fermion generations transform as 2 þ 1, one has a SUð2Þ breaking flavon and a Uð1Þ breaking flavon.In the supersymmetric realization, or upon imposing positive charge sums in the non-SUSY realizations, one finds texture zeros for the 11, 13 and 31 entries of the Yukawa matrices.The determinant is therefore given by the 12, 21 and 33 entries which are SUð2Þ singlets and therefore depend only on Uð1Þ charges, resulting in the same prediction for E=N when the Uð1Þ breaking flavon contains the axiflavon [and the SUð2Þ is gauged].
Finally we comment on the modification for the E=N range in the context of an additional light Higgs doublet, restricting for simplicity to the case of a 2HDM of type-II.Then Eq. ( 12) is modified by the rescaling v 6 → sin 3 β cos 3 β v 6 where tan β ¼ v u =v d is the ratio of Higgs vevs.Large values of tan β can reduce the suppression of the model-dependent in Eq. ( 14), and we find essentially the same 99.9% ranges for tan β ¼ 20, while for tan β ¼ 50 the range is slightly increased, E=N ∈ ½2.3; 3.0.

III. PHENOMENOLOGY
Being a QCD axion, the axiflavon is a very light particle with a large decay constant making it stable on cosmological scales.Assuming that the phase transition corresponding to the Uð1Þ H breaking happens before inflation, the energy density stored in the axion oscillations can be easily related to the present Dark Matter (DM) abundance [8][9][10]: For a given axion mass below roughly ≲10 −5 -10 −4 eV it is then always possible to choose a misalignment angle θ to get the correct dark matter abundance Ω DM h 2 ≈ 0.12.The axion domain wall problem is automatically solved in this setup, but interesting constraints can arise from isocurvature perturbations [40].
We show in Fig. 1 present and future bounds on the axiflavon both from axion searches and from flavor experiments in terms of its mass m a and its coupling to photons g aγγ .In this plane one can appreciate how the allowed range of E=N is considerably reduced compared to the standard axion window [32].Assuming that the axiflavon is also accounting for the total DM abundance we give the corresponding value of θ for a given mass.
In the high mass region with m a ∼ 0.1À10 meV stringent bounds on the axiflavon come from its coupling to fermions and are hence independent of g aγγ .A mild lower bound on the axiflavon decay constant f a can be derived from axiflavon coupling to electrons which affects white dwarf cooling [41].This bound cuts off our parameter space at around m a ≲ 10 meV.
A stronger bound comes from the flavor-violating coupling of the axiflavon to down and strange quarks, asd, leading to (bounds from kaon decays are more restrictive than the bounds from kaon mixing) where m K;π are the kaon and pion masses, and B s ¼ 4.6ð8Þ is the nonperturbative parameter related to the quark condensate [42].The 90% C.L. combined bound from E787 and E949, BRðK þ → π þ aÞ < 7.3 × 10 −11 [43], gives where κ sd =N ∼ Oð1Þ are model-dependent coefficients controlled by the particular flavor charge assignments, and quark masses are taken at μ ∼ 2 GeV.Similarly in the B sector we find with f K 0 ð0Þ ¼ 0.331 [44] and the shorthand notation where again κ bs =N ∼ Oð1Þ.A bound BRðB þ → K þ aÞ < 10 −6 ÷ 10 −8 , potentially in the reach of Belle II, would translate into m a < ð8 ÷ 80Þ meV × N=κ bs .A careful experimental analysis of this decay would be very interesting, as suggested also in Ref. [23].The solid blue line in Fig. 1 shows the lower bound on m a from flavor-violating kaon decays for κ sd =N ¼ 1.The reach on BRðK þ → π þ aÞ is expected to be improved by a factor ∼70 by NA62 [45,46] (and possibly also ORKA [47]   and KOTO [48]), giving sensitivity to scales as high as f a ≳ κ sd =N 6.3 × 10 11 GeV.The expected sensitivity on the axion mass for κ sd =N ¼ 1 is shown by the dashed blue line in Fig. 1.Therefore the future flavor experiment will probe the axiflavon parameter space in the interesting region where it can account for the dark matter relic abundance with θ ∼ Oð1Þ.
Going to lower axiflavon masses, below 0.1 keV, the phenomenology becomes essentially identical to the one of the original DFSZ model but with a sharper prediction for the value of E=N, given in Eq. ( 1).This corresponds to the brown band in Fig. 1.
The gray shaded regions in Fig. 1 summarize the present constraints on axionlike particles.An upper bound on the photon coupling for the full range of masses of our interest comes from its indirect effects on stellar evolution in Globular Clusters [49].A comparable bound is set by the CAST experiment [50].Stronger constraints for axions lighter than 0.1 μeV can be derived from the lack of gamma-ray signals emitted from the supernova SN1987A [51] and from the bounds on spectral irregularities of the Fermi-LAT and H.E.S.S. telescopes [52,53].The region of very low axion masses below 10 −5 μeV is disfavored by black hole superradiance independently on the photon coupling [54].In the axion mass region between and 100 μeV present bounds from the ADMX experiment [55] do not put yet a constraint on the axiflavon band.This is a well-known feature of the original DFSZ model with E=N ¼ 8=3 that is shared by the axiflavon and further motivates future developments in microcavity experiments.
In Fig. 1 we also display the projections for the different axion future experiments.The combination of the upgraded ADMX experiment and its high frequency version [56] can probe a wide range of the axiflavon parameter space in the mass window between 1 and 100 μeV.This region is strongly preferred because the correct axion abundance can be obtained without a tuning of the initial misalignment angle.Dielectric Haloscopes [57] have a similar reach of ADMX-HF and are not displayed in the plot.The IAXO experiment [58] gives instead a bound only at large axiflavon masses m a ≳ meV.Such large masses are already robustly ruled out by flavor-violating kaon decays.The low mass window of the axiflavon band for m a ≲ 0.1 μeV will be probed by the resonant ABRACADABRA experiment and its upgrade [59].Interestingly, the axiflavon band lives below the reach of the first phase of the broadband ABRACADABRA experiment.Axiflavon masses below 10 −3 μeV will eventually be probed in the final phase of the CASPEr experiment [60].
In conclusion, the axiflavon parameter space is considerably narrower than that of KSVZ/DFSZ models, as visible in Fig. 1, and will be covered in a wide range of masses by a combination of future axion searches and kaon experiments.In the high mass window with 10 −6 eV ≲ m a ≲ 10 −4 eV the comparable projected reaches of ADMX-HF and future kaon experiments leave the exciting possibility to tell apart the axiflavon scenario from other QCD axions.

FIG. 1 .
FIG.1.The axiflavon band (light brown) projected on the axion parameter space: mass vs photon coupling defined in Eq.(15).The standard KSVZ/DFSZ band is shown in light yellow.The grey exclusion region is obtained from the combination of various axion constraints that are summarized in the legend.The dashed colored lines show the projected reach of future axion experiments.The solid blue line is the exclusion reach from current flavor experiments for an axiflavon model with κ sd =N ¼ 1 [cf.Eq. (20)].The dashed blue line depicts the expected reach of future flavor experiments for the same choice of parameters.