TTP-20-005 Quark Flavor Phenomenology of the QCD Axion Jorge Martin Camalich,1, 2 Maxim Pospelov,3, 4 Pham Ngoc Hoa Vuong,5 Robert Ziegler,6, 7 and Jure Zupan8 1Instituto de Astrof´ısica de Canarias, C/ V´ıaL´actea, s/n E38205 - La Laguna, Tenerife, Spain 2Universidad de La Laguna, Departamento de Astrof´ısica, La Laguna, Tenerife, Spain 3School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455, USA 4William I. Fine Theoretical Physics Institute, School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455, USA 5Laboratoire de Physique Subatomique et de Cosmologie, Universit´eGrenoble-Alpes, CNRS/IN2P3, Grenoble INP, 38000 Grenoble, France 6Theoretical Physics Department, CERN, Geneva, Switzerland 7Institut f¨urTheoretische Teilchenphysik (TTP), Karlsruher Institut f¨urTechnologie (KIT), 76131 Karlsruhe, Germany 8Department of Physics, University of Cincinnati, Cincinnati, Ohio 45221,USA Axion models with generation-dependent Peccei-Quinn charges can lead to flavor-changing neutral currents, thus motivating QCD axion searches at precision flavor experiments. We rigorously derive limits on the most general effective flavor-violating couplings from current measurements and assess their discovery potential. For two-body decays we use available experimental data to derive limits on q ! q0a decay rates for all flavor transitions. Axion contributions to neutral-meson mixing are calculated in a systematic way using chiral perturbation theory and operator product expansion. We also discuss in detail baryonic decays and three-body meson decays, which can lead to the best search strategies for some of the couplings. For instance, a strong limit on the Λ ! na transition can be derived from the supernova SN 1987A. In the near future, dedicated searches for q ! q0a decays at ongoing experiments could potentially test Peccei-Quinn breaking scales up to 1012 GeV at NA62 or KOTO, and up to 109 GeV at Belle II or BES III. I. INTRODUCTION be performed. Besides the two-body meson decays, also the three-body and baryonic decays can provide the best sensitivity to specific axion couplings. We The QCD axion is arguably one of the best- provide a careful and rigorous analysis of the result- motivated particles beyond the Standard Model ing constraints by exploiting the entire set of current (SM). Originally predicted by the Peccei-Quinn experimental information. We also improve the pre- (PQ) solution to the strong CP problem [1{4], dictions for axion induced neutral-meson oscillations i.e., the apparent absence of large CP violation in by using effective-field theory methods, and derive strong interactions [5{12], the axion is also an ex- a new bound from supernova cooling. cellent cold Dark Matter (DM) candidate in large parts of the parameter space [13{15]. Predictions for axion couplings to fermions are The original electroweak axion models [3, 4] were model-dependent. The only requirement for a suc- strongly disfavored by the non-observation of flavor- cessful axion model is an almost exact global U(1) violating K+ π+a decays (where a is the ax- PQ symmetry that is spontaneously broken and ion) [16]. Since! then, experimental searches for anomalous with respect to QCD. Therefore, the only axions have been primarily focused on the flavor- coupling shared by all axion models is the axion ~ conserving axion couplings. Searches using halo- coupling to the CP-odd gluon operator, GG, aris- scopes [17{24] and helioscopes [25{28] rely on ax- ing from the QCD anomaly and responsible for the ion couplings to photons, while searches using pre- solution to the strong CP problem. Axion models divide into two classes, depending on how the color arXiv:2002.04623v1 [hep-ph] 11 Feb 2020 cision magnetometry rely on axion couplings to glu- ons [29]. Much activity is being devoted to de- anomaly arises. In the KSVZ-type models [35, 36] vise new experimental techniques, and ambitious the color anomaly is due to a set of heavy fermions projects have been proposed for the upcoming years that are vectorlike under the SM but chiral under (for a review see Ref. [30]). the PQ symmetry. In this case the axion does not couple to elementary SM fermions at tree level. In In this paper we reiterate the importance of fla- the DFSZ-type models [37, 38], on the other hand, vor violating transitions for axion searches. In- the SM fermions carry PQ charges, and the axion deed, a QCD axion with flavor violating couplings always couples to the fermionic currents. Whether could well be first discovered in flavor-physics ex- these couplings are flavor conserving or flavor vio- periments. While this possibility was already con- lating is a model-dependent choice. templated in the literature for rare meson decays and neutral meson mixing [31{34], we go well be- In the original DFSZ model the PQ charges are yond the state of the art. We show that, in light taken to be flavor universal [37, 38], so that flavor- of upcoming experiments, there are a number of ad- violating axion couplings arise only at loop level. ditional dedicated searches that could and should In general, though, flavor violating axion couplings 2 are present already at tree level. This is the case, where f is the axion decay constant, cV;A are her- a fifj for instance, in generalized DFSZ-type models with mitian matrices in flavor space, and the sum over generation dependent PQ charges [39{43], which can repeated generational indices, i; j = 1; 2; 3, is im- also allow to suppress the axion couplings to nucle- plied. For future convenience we define effective de- ons [44{46]. Particularly motivated scenarios, which cay constants as lead to flavor violating axion couplings at tree level, arise when the PQ symmetry is part of a flavor 2fa F V;A : (2) group that shapes the structure of the yukawa sec- fifj ≡ cV;A tor [31, 47]. The PQ symmetry could enforce texture fifj zeros in the Yukawa matrices [48{50], or be respon- ∗ In general F V;A; i = j, are complex, with F V;A = sible for their hierarchical structure `ala Froggatt- fifj 6 fifj F V;A. Throughout the paper we take a to be the Nielsen (FN) [51]. While in the simplest scenario fj fi PQ and FN symmetries are identified [52{54], PQ QCD axion, so that its mass is inversely propor- could also be a subgroup of a larger flavor symmetry, tional to fa [79], see e.g. Refs. [55{61]. Finally, flavored PQ symme- tries can arise also in the context of Minimal Fla- 1012 GeV ma = 5:691(51) µeV : (3) vor Violation [62, 63] or as accidental symmetries in fa models with gauged flavor symmetries [64{67]. In our analysis we remain, for the most part, ag- For the \invisible" axion the decay constant is fa nostic about the origin of the flavor and chiral struc- 106 GeV [80], in which case the axion is much lighter ture of axion couplings to SM fermions, and sim- than an eV and essentially decoupled from the SM. ply treat axion couplings to fermions as independent We will always be working in this limit, so that in parameters in an effective Lagrangain. For related the flavor transitions the axion can be taken as mass- studies of axion-like particles with flavor violating less for all practical purposes. couplings, see [68{70] (for loop induced transitions In this mass range the axion has a lifetime that see [71{78]). We restrict the analysis to the case of is larger than the age of the universe, and there- the (practically) massless QCD axion, but our re- fore is a suitable DM candidate. If the PQ symme- sults can be repurposed for any other light scalar or try is broken before inflation, axions are produced pseudoscalar with flavor violating couplings to the near the QCD phase transition and yield the ob- SM fermions, as long as the mass of the (pseudo- served DM abundance for axion decay constants of 11 13 )scalar is much smaller than the typical energy re- the order fa (10 10 ) GeV [13{15], assum- lease in the flavor transition. ing natural values∼ of the÷ misalignment angle. Other The paper is structured as follows. In Section II production mechanisms, e.g., via parametric reso- we introduce our notation for the axion couplings nance, allow for axion DM also for smaller decay 8 to fermions and comment on their flavor structure. constants, down to fa 10 GeV [81]. We will see In Section III we derive the bounds on these cou- below that precision flavor∼ experiments are able to plings from two-body and three-body meson decays, test this most interesting region of the QCD axion from baryon decays and from baryon transitions in parameter space. supernovae. Section IV contains bounds from mix- The axion couplings to the SM fermions in the ing of neutral mesons, Section V reviews bounds on mass basis, cV and cA , are related to the PQ fifj fifj flavor-diagonal couplings, and Section VI discusses charge matrices in the flavor basis, Xf , through axion couplings involving the top quark. Finally, in Section VII we present the results and experi- V;A 1 y y cf f = Vf XfR VfR Vf XqL VfL ; (4) mental projections. Details about renormalization i j 2N R ± L ij of effective axion couplings, experimental recasts of two-body meson decays and hadronic inputs are de- where N is the QCD anomaly coefficient of the ferred to the Appendix. PQ symmetry. The unitary rotations VfL;fR diag- onalize the appropriate SM fermion yukawa matri- ces, V y y V = ydiag, for the \up" and \down" fL f fR f II.