PHYSICAL REVIEW D 102, 015023 (2020)
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, s/n, E-38205, La Laguna, Tenerife, Spain 2Departamento de Astrofísica, Universidad de La Laguna, E-38206, La Laguna, Tenerife, Spain 3School of Physics and Astronomy, University of Minnesota, Minneapolis, Minnesota 55455, USA 4William I. Fine Theoretical Physics Institute, School of Physics and Astronomy, University of Minnesota, Minneapolis, Minnesota 55455, USA 5Laboratoire de Physique Subatomique et de Cosmologie, Universit´e Grenoble-Alpes, CNRS/IN2P3, Grenoble INP, 38000 Grenoble, France 6Theoretical Physics Department, CERN, 1 Esplanade des Particules, 1211 Geneva 23, Switzerland 7Institut für Theoretische Teilchenphysik (TTP), Karlsruher Institut für Technologie (KIT), 76131 Karlsruhe, Germany 8Department of Physics, University of Cincinnati, Cincinnati, Ohio 45221, USA
(Received 20 February 2020; accepted 24 June 2020; published 27 July 2020)
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.
DOI: 10.1103/PhysRevD.102.015023
I. INTRODUCTION gluons [30]. Much activity is being devoted to devise new experimental techniques, and ambitious projects have been The QCD axion is arguably one of the best-motivated proposed for the upcoming years (for a review, see Ref. [31]). particles beyond the Standard Model (SM). Originally In this paper, we reiterate the importance of flavor- predicted by the Peccei-Quinn (PQ) solution to the strong violating transitions for axion searches. Indeed, a QCD CP problem [1–4], i.e., the apparent absence of large CP axion with flavor-violating couplings could well be first violation in strong interactions [5–13], the axion is also an discovered in flavor-physics experiments. While this pos- excellent cold dark matter (DM) candidate in large parts of sibility was already contemplated in the literature for rare the parameter space [14–16]. meson decays and neutral meson mixing [32–35],wego The original electroweak axion models [3,4] were well beyond the state of the art. We show that, in light of strongly disfavored by the nonobservation of flavor- upcoming experiments, there are a number of additional violating Kþ → πþa decays (where a is the axion) [17]. dedicated searches that could and should be performed. Since then, experimental searches for axions have been Besides the two-body meson decays, also the three-body primarily focused on the flavor-conserving axion couplings. and baryonic decays can provide the best sensitivity to Searches using haloscopes [18–25] and helioscopes specific axion couplings. We provide a careful and rigorous [26–29] rely on axion couplings to photons, while searches analysis of the resulting constraints by exploiting the entire using precision magnetometry rely on axion couplings to set of current experimental information. We also improve the predictions for axion-induced neutral-meson oscilla- tions by using effective-field theory methods and derive a Published by the American Physical Society under the terms of new bound from supernova cooling. the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to Predictions for axion couplings to fermions are model the author(s) and the published article’s title, journal citation, dependent. The only requirement for a successful axion and DOI. Funded by SCOAP3. model is an almost exact global Uð1Þ PQ symmetry that is
2470-0010=2020=102(1)=015023(24) 015023-1 Published by the American Physical Society JORGE MARTIN CAMALICH et al. PHYS. REV. D 102, 015023 (2020) spontaneously broken and anomalous with respect to QCD. from mixing of neutral mesons, Sec. V reviews bounds Therefore, the only coupling shared by all axion models is on flavor-diagonal couplings, and Sec. VI discusses axion the axion coupling to the CP-odd gluon operator, GG˜ , couplings involving the top quark. Finally, in Sec. VII,we arising from the QCD anomaly and responsible for the present the results and experimental projections. Details solution to the strong CP problem. Axion models divide about renormalization of effective axion couplings, exper- into two classes, depending on how the color anomaly imental recasts of two-body meson decays, and hadronic arises. In the KSVZ-type models [36,37] the color anomaly inputs are deferred to the Appendix. is due to a set of heavy fermions that are vectorlike under the SM but chiral under the PQ symmetry. In this case, the II. AXION COUPLINGS TO FERMIONS axion does not couple to elementary SM fermions at tree level. In the DFSZ-type models [38,39], on the other hand, The Lagrangian describing the most general interactions 1 the SM fermions carry PQ charges, and the axion always of the axion with the SM fermions is given by (see also couples to the fermionic currents. Whether these couplings Appendix A) are flavor conserving or flavor violating is a model- ∂ dependent choice. μa ¯ μ V A L ¼ f γ ðc þ c γ5Þf ; ð1Þ aff i fifj fifj j In the original DFSZ model, the PQ charges are taken 2fa to be flavor universal [38,39], so that flavor-violating axion couplings arise only at loop level. In general, though, where f is the axion decay constant; cV;A are Hermitian a fifj flavor-violating axion couplings are present already at matrices in flavor space; and the sum over repeated tree level. This is the case, for instance, in generalized generational indices, i; j ¼ 1; 2; 3, is implied. For future DFSZ-type models with generation-dependent PQ charges convenience, we define the effective decay constants as [40–44], which can also allow one to suppress the axion 2 couplings to nucleons [45–47]. Particularly motivated fa FV;A ≡ : ð2Þ scenarios, which lead to flavor-violating axion couplings fifj V;A cf f at tree level, arise when the PQ symmetry is part of a i j flavor group that shapes the structure of the Yukawa In general, FV;A , i ≠ j, are complex, with ðFV;A Þ ¼ FV;A . sector [32,48]. The PQ symmetry could enforce texture fifj fifj fjfi zeros in the Yukawa matrices [49–51] or be responsible Throughout the paper, we take a to be the QCD axion, so for their hierarchical structure ala` Froggatt-Nielsen that its mass is inversely proportional to fa [83], (FN) [52]. While in the simplest scenario PQ and FN 1012 symmetries are identified [53–55], PQ could also be GeV ma ¼ 5.691ð51Þ μeV : ð3Þ a subgroup of a larger flavor symmetry; see, e.g., fa Refs. [56–65]. Finally, flavored PQ symmetries can arise also in the context of minimal-flavor violation (MFV) For the “invisible” axion, the decay constant is fa ≫ [66,67] or as accidental symmetries in models with gauged 106 GeV [84], in which case the axion is much lighter flavor symmetries [68–71]. than an eV and essentially decoupled from the SM. We will In our analysis, we remain, for the most part, agnostic always be working in this limit, so that in the flavor about the origin of the flavor and chiral structure of axion transitions the axion can be taken as massless for all couplings to SM fermions and simply treat axion couplings practical purposes. to fermions as independent parameters in an effective In this mass range, the axion has a lifetime that is larger Lagrangain. For related studies of axionlike particles with than the age of the Universe and therefore is a suitable DM flavor-violating couplings, see Refs. [72–74] (for loop- candidate. If the PQ symmetry is broken before inflation, induced transitions, see [75–82]). We restrict the analysis to axions are produced near the QCD phase transition and the case of the (practically) massless QCD axion, but our yield the observed DM abundance for axion decay con- 11 13 results can be repurposed for any other light scalar or stants of the order fa ∼ ð10 ÷ 10 Þ GeV [14–16], assum- pseudoscalar with flavor-violating couplings to the SM ing natural values of the misalignment angle. Other fermions, as long as the mass of the (pseudo)scalar is much production mechanisms, e.g., via parametric resonance, smaller than the typical energy release in the flavor allow for axion DM also for smaller decay constants, down 8 transition. to fa ∼ 10 GeV [85]. We will see below that precision The paper is structured as follows. In Sec. II,we flavor experiments are able to test this most interesting introduce our notation for the axion couplings to fermions region of the QCD axion parameter space. and comment on their flavor structure. In Sec. III, we derive the bounds on these couplings from two-body and three- 1Note that diagonal vector couplings are unphysical up to body meson decays, from baryon decays, and from baryon electroweak anomaly terms, which are irrelevant for the purpose transitions in supernovae. Section IV contains bounds of this paper.
015023-2 QUARK FLAVOR PHENOMENOLOGY OF THE QCD AXION PHYS. REV. D 102, 015023 (2020)
The axion couplings to the SM fermions in the mass V;A simply take cf f to be unknown parameters in an effective basis, cV and cA , are related to the PQ charge matrices i j fifj fifj Lagrangian, which will be constrained solely from data. in the flavor basis, Xf, through III. BOUNDS FROM HADRON DECAYS 1 V;A † † c ¼ ðV Xf Vf V Xq Vf Þ ; ð4Þ Bounds on the vector and axial-vector parts of the flavor- fifj 2 fR R R fL L L ij N violating axion couplings, Eqs. (1) and (2), can be derived from searches for hadron decays with missing energy. where N is the QCD anomaly coefficient of the PQ sym- In this section, we consider the bounds from two-body metry. The unitary rotations V diagonalize the appro- fL;fR decays of pseudoscalar mesons to pseudoscalar and vector † priate SM fermion Yukawa matrices, V y V ¼ ydiag, for mesons, respectively, from three-body meson decays and fL f fR f the “up” and “down” quark flavors, f ¼ u, d. We focus on from decays of baryons. In each case, we first review the axion couplings to quarks and refer the reader to Ref. [86] available and planned experimental measurements and then for present and future prospects for testing lepton flavor– interpret them in terms of the bounds on flavor-violating violating axion couplings. Off-diagonal couplings arise axion couplings. whenever PQ charges, XqL and XfR , are not diagonal in the A. Bounds from two-body meson decays same basis as the Yukawa matrices, yf. Their sizes depend on the misalignment between the two bases, parametrized Because of parity conservation of strong interactions, → by the unitary rotations VfL and VfR (taking XqL and XfR to the P1 P2a decays of a pseudoscalar meson P1 to a be diagonal). pseudoscalar meson P2 are only sensitive to the vector Very different flavor textures of cV;A are possible. couplings of the axion, while the P1 → V2a decays, where fifj V2 is the vector meson, are only sensitive to the axial-vector Provided a suitable set of PQ charges and appropriate flavor couplings (see Appendix C). Searches targeting specifically structures of the SM Yukawa matrices, it is possible for just a þ → πþ V;A the massless axion were performed in K a [89], single off-diagonal coupling to be large c ∼ Oð1Þ, i ≠ j, þ þ þ þ fifj B → K a [90], and B → π a [90] decays. In addition, with all the other off-diagonal couplings zero or very small. searches for SM decays where the invisible final state is V A ∼ 1 For example, one can realize a situation where cbs ¼ cbs a νν¯ pair can be recast to derive limits on the axion V;A 0 1 couplings. This requires that the two-body kinematics of an and all the other cf ≠f ¼ ,bytakingXqL ¼ XuR ¼ , i j (essentially) massless axion is included in the search ðX Þ2 ≠ ðX Þ3, with down Yukawa matrix y such that the dR dR d region, as was the case in the BABAR and CLEO searches only nonvanishing rotation is in the 2-3 right-handed sector, → ð Þνν¯ → πνν¯ → τ → πν¯ ν sRd ∼ 1. Moreover, while one would generically expect axial for B K [91], B [92], and D ð Þ 23 [93]. Note that the corresponding Belle datasets analyzed in couplings cA and vector couplings cV to be of the same fifj fifj Refs. [94,95] cannot be readily used to set bounds on axion X order, the latter can be suppressed in a situation where fR ¼ couplings, because the analyses either cut out two-body − XqL and VfR ¼ VfL , which can arise in models where PQ decays with a massless axion [95] or used multivariate charges are compatible with a grand unified structure [see methods [94] that are difficult to reinterpret for different Ref. [87] for a recent example in SO(10)], and Yukawas are purposes [96] (see also Ref. [97]). Hermitian, positive-definite matrices [see, e.g., Ref. [88] for The available experimental information is summarized in a realization of this scenario in SO(10)]. Table I, where we give the 90% C.L. upper limits on the In the absence of a theory of flavor, we will be agnostic branching ratios for decays involving neutrinos or invisible about the origin of the possible flavor misalignment and massless axions. We include the limits on the decays
TABLE I. Experimental inputs for meson decays; see text for details. We show the 90% C.L. upper bounds on the branching ratios of a pseudoscalar meson P1 to another pseudoscalar (P2) or vector (V2) meson (for sd transitions V2 ¼ ππ instead). The bounds shown are for decays to neutrinos or massless invisible axions. In the latter case, we also show our bounds obtained by recasting related searches for invisible decays (subscript “recast”).
Decay sd cu bd bs −11 −5 −5 BRðP1 → P2 þ aÞ 7.3 × 10 [89] No analysis 4.9 × 10 [90] 4.9 × 10 [90] → −6 −5 −6 BRðP1 P2 þ aÞrecast No need 8.0 × 10 [93] 2.3 × 10 [92] 7.1 × 10 [91] → νν¯ 1 47þ1.30 10−10 0 8 10−5 1 6 10−5 BRðP1 P2 þ Þ . −0.89 × [89] No analysis . × [94] . × [94] −5 BRðP1 → V2 þ aÞ 3.8 × 10 [98] No analysis No analysis No analysis → −5 BRðP1 V2 þ aÞrecast No need No data No data 5.3 × 10 [91] −5 −5 −5 BRðP1 → V2 þ νν¯Þ 4.3 × 10 [98] No analysis 2.8 × 10 [94] 2.7 × 10 [94]
015023-3 JORGE MARTIN CAMALICH et al. PHYS. REV. D 102, 015023 (2020) involving axions that we obtained by recasting the exper- D0 → ρ0a, which results in three charged pions þ MET imental searches for decays involving neutrinos (“recast”). and two displaced vertices. The lack of such analyses Table I shows bounds on the branching ratios for decays to means that there is at present no bound from meson decays þ þ pseudoscalar mesons, P1 → P2a, for K → π a (s → d on axial cu couplings to the axion. Similarly, there is at transition, experimental analysis in Ref. [89]), Dþ → πþa present no publicly available experimental analysis that (c → u transition, our recast of Ref. [93]), Bþ → πþa bounds the B → ρa decays (as discussed above, one cannot (b → d transition, experimental analysis in Ref. [90] and readily use for that purpose the B → ρνν¯ Belle data from our recast of Ref. [92]), and Bþ → Kþa (b → s transition, Ref. [94], while BABAR has not performed such an experimental analysis in Ref. [90] and our recast of analysis). Finally, our recast bounds on B → Kð Þa, B → Ref. [91]). For decays to vector mesons, P1 → V2a, there πa could be easily improved by dedicated experimental is experimental information on the B → K a decay (b → s searches using already collected data. At LHCb, one could transition, our recast of Ref. [91]). In the same section of measure the B → K a and B → ρa branching ratios using 0 Table I, we also include the bounds on Kþ → πþπ a decay ¯0 þ − ¯0 þ − the decay chains such as Bs → K B or B → π B used in Sec. III B below (s → d transition, experimental − − − ¯0 ¯0 followed by B → K ð→ KSπ Þa,orBs → KSB fol- analysis in Ref. [98]). For details on the recast, see lowed by B¯0 → K¯ 0a, ρ0a [105]. One could also attempt Appendix B. more challenging decay chain measurements such as The above analyses could be improved with dedicated Bs → Bsγ, followed by Bs → ϕa or Bs → K a. axion searches applied to existing data or to ongoing We now convert the bounds on the branching ratios in þ → πþ experiments. Sensitivity to the K a decay better Table I to bounds on flavor-violating couplings of axions to than the present world best limit can be achieved at the quarks, Eqs. (1) and (2). The corresponding partial decay NA62 experiment. For a massless axion, an improvement widths are given by by an order of magnitude compared to the BNL result, −11 þ → πþ 7 3 10 8 0 2 BRðK aÞ < . × [89], is expected by 2025 fþð Þ < 2 ;P1 → P2a jFV j [99]. The same flavor transition is probed by KOTO Γ κ ij → π0 ¼ 12 2 ð5Þ searching for the neutral decay mode KL a. The : A0ð0Þ 2 ;P1 → V2a jFA j current limit using data collected in 2015 is BRðKL → ij π0aÞ < 2 × 10−9 and is in the same ballpark as for the decay to the neutrino pair [100]. The KOTO Collaboration with the kinematic prefactor expects the sensitivity to be improved down to the 10−11 M3 M2 3 level [101]. Improved sensitivity in the neutral-kaon κ 1 1 − 2 12 ¼ 16π 2 ; ð6Þ mode can also be expected from the proposed KLEVER M1 experiment [102]. To the best of our knowledge, the future prospects for the where M1 (M2) is the mass of the parent (daughter) meson. → π0 heavy-meson axion decays P1 → P2a or P1 → V2a have Since KL a decay is CP violating, the partial decay not yet been estimated by the experimental collaborations. width in that case is given by Nevertheless, it is clear that improvements over the present 2 V 2 Γ →π0 ¼ κ12f ð0Þ ½Imð1=F Þ ð7Þ situation are justifiably expected. For instance, the exper- KL a þ sd imental error on the Dþ → τþν branching ratio is projected −1 V 0 to be reduced by a factor of 3 with 20 fb integrated and thus vanishes in the CP-conserving limit, ImFsd ¼ ; → π0 þ → πþ luminosity at BESIII [103] compared to the present value cf. Eq. (2). The KL a and K a decay rates obey → π0 ≤ 4 3 þ → [104], indicating a potential significant improvement in the Grossman-Nir bound BRðKL aÞ . BRðK sensitivity to the c → ua transition. The reach on the πþaÞ [106,107]. 2 2 branching ratios for the axionic decay modes of B mesons The form factors fþðq Þ and A0ðq Þ are defined in will be improved with the amount of data expected to be Appendix C, where we also collect the numerical values collected at Belle II, which is roughly a factor of 50 larger used as inputs in the numerical analysis. The resulting V;A compared to the BABAR and Belle samples. bounds on axion couplings Fij are shown in Table III. The Several potentially interesting channels are lacking any implications of these results and future projections will be experimental analyses so far. For example, there is no discussed in Sec. VII. experimental analysis of c → ua transitions that are sensi- tive to the axial-vector coupling; i.e., there are no D → ππ → ρ νν¯ B. Bounds from three-body meson decays Xinv or D Xinv, Xinv ¼ ;a, searches. One could also search for a c → ua signal in D → Ka, D → K a The E787 experiment at Brookhaven performed a search s s 0 decays, all of which could be performed at Belle II and for the three-body Kþ → π πþa decay mediated by the BESIII. Potentially, LHCb could also probe these couplings s → da transition and set the bound BRðKþ → π0πþaÞ ≤ using decay chains, such as B− → D0π− followed by 3.8 × 10−5 at 90% C.L. [98]. The related decay mode
015023-4 QUARK FLAVOR PHENOMENOLOGY OF THE QCD AXION PHYS. REV. D 102, 015023 (2020)
→ π0π0 1 A → 1 A KL a has also been searched for, resulting in the obtained by replacing Reð =FsdÞ Imð =FsdÞ in the rates 0 0 upper limit for light massive axions BRðKL → π π aÞ ≲ of KL. 0.7 × 10−6 [108]. However, this analysis excluded the Better sensitivity to the axial sd axion coupling can be → π0π0 ma ¼ 0 kinematic region and is thus not applicable to achieved by searching for KL a at KOTO [109], the case of the QCD axion [109]. Other decay modes such while a search for Kþ → πþπ0a could be attempted at þ − as KL → π π a or those involving the decays of KS have NA62. Note that Eqs. (8) and (9) can be combined to the not been investigated experimentally. inequality Parity conservation implies that the K → ππa decays are 0 0 þ þ 0 sensitive only to the axial-vector couplings of the axion to BRðKL → π π aÞ ≤ 31BRðK → π π aÞ; ð11Þ quarks (see Appendix C). The form factors entering the τ τ 4 1 predictions are related via isospin symmetry to the form which, besides the ratio of kaon lifetimes L= þ ¼ . [that factors measured in Kþ → πþπ−eþν [110–113], making gives the main effect in the Grossman-Nir bound for precise predictions for K → ππa decay rates possible. The BRðK → πaÞ], also includes the effects of different com- two final-state pions can be only in the total isospin I ¼ 0 binations of Clebsch-Gordan coefficients, form factors, and or the I ¼ 1 state, since the s → da Lagrangian is phase space factors. jΔIj¼1=2, while the initial kaon is part of an isodoublet. The numerical input values for the hadronic matrix Bose symmetry demands the decay amplitude to be elements are described in Appendix C, while the resulting symmetric with respect to the exchange of the two pions. bound on the axion coupling obtained from the current The I ¼ 0 (I ¼ 1) amplitude is even (odd) under this bound on Kþ → πþπ0a [98], is included in Table III. The permutation. The form factors must therefore enter in implications of these results and prospects are discussed combinations which are even (odd) with respect to the further in Sec. VII. exchange of pion momenta, pπ ↔ pπ . The two pions in Other searches using three-body decays could be of 1 2 0 0 0 0 the decay K0 → π0π0a (Kþ → πþπ0a) are in a pure I ¼ 0 interest. Decays such as Bþ → ρ πþa, B → ρ ρ a or → ρ0 → 0ρ0 → ϕρ0 (I ¼ 1) state, and one obtains Bs KS a, Bs K a, Bs a could potentially be even attempted at LHCb, since they result in three or 0 0 2 3 dΓðK → π π aÞ ðm 0 − sÞ four charged particles and a massless invisible particle. L ¼½Reð1=FA Þ 2 K βF2 ð8Þ ds sd 1024π3m5 s Other possible decays of theoretical interest, but probably K only measurable at Belle II, are Bþ → π0πþa, B0 → − 0 − 0 and πþπ a, B → ρþρ a, and Bþ → ρþρ a. Belle II could also access from the ϒð5SÞ run the Bs decays with neutral 0 0 0 þ þ 0 2 − 3 → π → π dΓðK → π π aÞ 1 ðm þ sÞ pions such as Bs KS a, Bs K a, etc. Providing ¼ K β ds jFA j2 1536π3m5 predictions for these decays lies beyond the scope of the sd K present paper, though controlled calculations using QCD 2 β2 2 2β × ðFp þ Gp þ FpGpÞ; ð9Þ factorization and soft-collinear effective theory may be possible [116,117] and could be attempted in the future (see 2 β2 1– 4 2 where s ¼ðpπ1 þ pπ2 Þ and ¼ mπ=s. The functions also related work on two-body form factors [118–121]). Fs, Fp, and Gp correspond to the moduli of the coefficients Until then, one can use as a rough guide the bounds that in the partial-wave expansion of the form factors (which were obtained for the related two-body decays, i.e., Bþ → are complex functions of the kinematic variables of the πþa (on FV Þ and Bþ → ρþa (on FA ), as a ballpark figure − bd bd two-pion system) in the Kþ → πþπ eþν decay channel for what an interesting experimental reach for Bþ → ρ0πþa [110,111,114,115] (see Appendix C). In the first equation, V A (bounding a combination of Fbd and Fbd) might be. we also used the fact that the final state is CP odd, and Finally, one can also use LHCb dimuon data collected therefore the decay occurs through the predominantly for the Bq → μμ analysis to constrain the axial couplings CP-odd component of KL. A A Fbd and Fbs. As long as no vetoes on extra particles in the One can also obtain expressions for other decay modes event are applied in the LHCb analysis, their datasets can be → πþπ− from Eqs. (8) and (9). The amplitude of KL a used to constrain decays with additional particles in the 0 1 contains both I ¼ and I ¼ isospin components, and the final state such as axions. In Ref. [74], the present data decay rate in the isospin limit is A were used to derive constraints on Fbq of the order of 5 −1 − 0 10 GeV, cf. Table III. With 300 fb , the bounds can be ΓðK → πþπ aÞ¼Γ˜ðKþ → πþπ aÞ L strengthened by about an order of magnitude, and could be 0 0 þ 2ΓðKL → π π aÞ; ð10Þ further improved by using also the ATLAS and CMS data on Bq → μμ. The same strategy might also be applied to Γ˜ þ → πþπ0 1 A → where ðK aÞ is obtained by replacing =jFsdj sd transitions using KS → μμ decays, as proposed in 1 A Γ þ → πþπ0 Reð =FsdÞ in ðK aÞ. Rates for the KS decays are Ref. [122].
015023-5 JORGE MARTIN CAMALICH et al. PHYS. REV. D 102, 015023 (2020)
C. Bounds from baryon decays D. Supernova bound
Baryon decays, B1 → B2a, are sensitive to both axial and In the core of neutron stars (NS), hyperons coexist vector couplings of the axion. The decay rates are given by in equilibrium with neutrons, protons, and electrons – Λ → [130 133]. The decay na would represent a new 2 2 f1ð0Þ g1ð0Þ cooling mechanism for NS and can thus be constrained ΓðB1 → B2aÞ¼κ12 þ ; ð12Þ jFV j2 jFA j2 by stellar structure calculations and observations. At ij ij exactly zero temperature, the degenerate Λ and neutron distributions must have the same Fermi energy, leaving with the kinematic prefactor κ12 given in Eq. (6). The Λ → 2 2 no phase space for the na decays to occur. The form factors f1ðq Þ and g1ðq Þ are discussed in detail in degeneracy is partially lifted at finite temperature allowing Appendix C. for the Λ → na transitions with a rate that increases At present, there are no published experimental searches with the temperature. The impact of this new cooling → for B1 B2a decays. We therefore set 90% C.L. upper mechanism is maximal during the few seconds after → limits on BRðB1 B2aÞ indirectly. For hyperons, we the supernova (SN) explosion, when a protoneutron star subtract from unity the branching ratios for all channels (PNS) reaches temperatures of several tens of MeV that have been measured so far, adding the experimental [134,135]. Λ errors in quadrature [10].For b decays, we use the SM To estimate the cooling facilitated by the sd-axion prediction for its lifetime at next-to-leading order (NLO) in interaction in this early phase of the supernova evolution, α O 1 s and at ð =mbÞ in heavy quark expansion [123], com- we assume that the PNS is a system of noninteracting (finite pare it to the experimental measurements [10], and ascribe temperature) Fermi gases of neutrons, protons, electrons, Λ → the difference to the allowed value for BRð b B2aÞ.For and Λ baryons that are in thermal and chemical equilib- Λ Λþ → c, we saturate the observed lifetime with the c pa rium. Furthermore, we assume that the neutrinos are decay width. The resulting upper bounds on the branching trapped inside the PNS, while the lepton fraction number, → ratios BRðB1 B2aÞ are collected in Table II and used to relative to baryon number density, is taken to be Y ¼ 0.3 V;A L derive the limits on axion couplings Fij shown in [136]. The occupancy of Λ states is distributed according Table III. Λ 1 1 EΛ−μΛ to the Fermi distribution fp ¼ =ð þ expð T ÞÞ, where The sensitivity could be improved substantially with p is the Λ 3-momentum in the star’s rest frame, EΛ is its dedicated searches for B1 → B2a decays. The BESIII 2 2 2 energy, E p m , and μΛ is its chemical potential. → νν¯ Λ ¼ þ Λ Collaboration plans to measure B1 B2 decays of Neutrons are distributed following an analogous distribu- hyperons by using a sample of hyperon-antihyperon pairs n tion, f 0 , characterized by μ and labeled by the corre- collected at the eþe− → J=ψ peak [127] (see also p n sponding neutron 3-momentum p0, also in the star’s frame. Ref. [128]). At BESIII, one could also study the decays Antiparticles follow identical distributions with the replace- of charmed baryons into axions. Namely, approximately μ → −μ 104 ΛþΛ¯ − 567 −1 ment , so that for the temperatures expected in a c c pairs have been collected in a run with pb ¯ − PNS the densities of Λ and n¯ are negligible. of integrated luminosity at eþe collisions just above the The volume emission rate Q inside the PNS due to the pair-production threshold [129]. Note that bottom baryons − process Λ → na is given by are not produced in the modern eþe machines, since they run at energies below the corresponding pair-production Z 3 Γ Λ → ∞ thresholds, so only LHCb, while challenging, could have mΛ ð naÞ Q ¼ 2 2 2 pdp access to these decays. π ðmΛ − m Þ 0 Z n 0 pmax EΛ − E 0 0 n Λ 1 − n × p dp fp ð fp0 Þ; ð13Þ 0 EΛE TABLE II. The 90% C.L. upper bounds on the branching pmin n fractions for the baryon B1 → B2a decays obtained by adding up the measured branching fractions of the exclusive modes (hyper- 0 0 ons) or by comparing theory predictions for lifetimes with the where pmax (pmin) is the maximal (minimal) neutron measurements (heavy baryons). momentum in the Λ → na decay, if Λ has momentum p ¼jpj (all in the PNS’s rest frame).2 Notice that in the Baryon B1 BRðB1 → B2aÞ90% 0 nonrelativistic limit where p; p ≪ mΛ ∼ mn, and in the Λ 8.5 × 10−3 limit of no Pauli blocking of the final state neutrons this Σþ 4.9 × 10−3 formula reduces to a more familiar form, Ξ0 2.3 × 10−4 Ξ− 6 4 10−4 . × 2In the star’s rest frame Λ is moving in the direction of pˆ. The ’ Λc 1 maximum (minimum) 3-momentum of the neutron in the PNS s −2 rest frame is reached when the neutron recoils in the Λ’s rest Λb 4.1 × 10 frame in the direction (in the direction opposite) of pˆ.
015023-6 QUARK FLAVOR PHENOMENOLOGY OF THE QCD AXION PHYS. REV. D 102, 015023 (2020)
−mΛ−mn O 2 Q ≃ nnðmΛ − mnÞΓðΛ → naÞe T ; ð14Þ Notice that ða Þ terms that also appear in the trans- formation from Eq. (1) to Eq. (17) do not affect ΔF ¼ 2 where nn is the number density of neutrons. processes and are omitted for that reason. The above L Evaluating the distributions (chemical potentials) for form of aff simplifies somewhat the calculations of the benchmark conditions of T ¼ 30 MeV and nuclear density nonlocal meson mixing matrix elements, Eq. (16). For this, ρ ¼ ρ , using Eq. (13), we obtain for the energy loss per we utilize the appropriate effective-field theories; we use nuc ¯ unit mass ϵ ¼ Q=ρ, chiral perturbation theory (ChPT) for contributions to K−K ¯ mixing, while for the heavy-quark systems, Bd;s − Bd;s and 0 2 0 2 ¯ ϵ 3 6 1038 erg 2 f1ð Þ g1ð Þ D − D, we use the operator product expansion (OPE), ¼ . × GeV V 2 þ A 2 : ð15Þ sg jFsdj jFsdj matching onto local four-quark operators describing the meson mixing. In the following, we use the relativistic Setting as the maximal limit on ϵ the energy lost through 0 0 0 2 δ ⃗− ⃗0 neutrino emission 1 sec after the collapse of the supernova normalization of states, hP ðkÞjP ðk Þi ¼ E ðk k Þ, and 0 − ¯0 SN 1987A, ϵ ≲ 1019 erg=sg[136,137], one obtains bounds the phase convention CPjP i¼ jP i. on jFA j and jFV j in the range 109—1010 GeV. sd sd ¯ Our estimates are afflicted by significant uncertainties. A. K−K mixing Nuclear interactions induce important corrections in the Since ma ≪ 1 GeV, we can use ChPT to describe calculation of the number densities [131,133], and there are contributions from axion exchanges in K−K¯ mixing. For considerable stellar uncertainties stemming from the com- axial couplings, cA , the leading contributions arise at tree plex physics at work in the supernova. Note that the energy fifj level, while for vector couplings, cV , the first nonzero loss per unit mass obtained using the approximate formula fifj in Eq. (14) is independent of the structural details of the contributions are at one loop. PNS, except for the temperature. At T ¼ 30 MeV, this To construct the ChPT Lagrangian in the presence leads to an emission rate that is approximately 40% larger of flavor-violating axions, we use a spurion analysis than in Eq. (15). More than anything, the emission rate [139,140]. In terms of the (pseudo)scalar interactions, suffers from the uncertainty in the temperature of the the Lagrangian for QCD with a flavor-violating axion central region. Variation of this quantity from 20 to can be conveniently written as 40 MeV changes Q by 2 orders of magnitude. Finally, L ¯ ∂ a a − ¯M our bound crucially relies on the validity of the standard QCDþa ¼ qði= þ gs=G T Þq q qq scenario for the SN explosion as applied to SN 1987A, − aq¯ðχ − iχ γ5Þq; ð19Þ which was disputed in a recent publication [138]. S P where we keep only light quarks, q ¼ðu; d; sÞ. The IV. BOUNDS FROM MESON MIXING diagonal mass matrix is Mq ¼ diagðmu;md;msÞ, while χ 3 3 The exchanges of axions with flavor-violating couplings S;P are × Hermitian matrices describing the quark- contribute to ΔF ¼ 2 transitions and can modify meson axion couplings, 0 1 mixing rates from the SM predictions. The contribution from 00 0 axion exchanges to the mixing amplitude of the P0 − P¯0 B − C 00md ms neutral meson system is given by the time-ordered correlator χ B V C S ¼ i@ Fds A; ð20Þ Z − 0 ms md 0 i 4 0 0 FV − L L 0 ¯ sd M12 ¼ 4 d xhP jTf affðxÞ; affð ÞgjP i; ð16Þ mP 0 1 2mu A 00 L Fuu where affðxÞ is the axion-fermion interaction Lagrangian B C B 2 C in Eq. (1). χ − 0 md mdþms P ¼ B FA FA C: ð21Þ @ dd ds A In this section, it will prove useful to use the following 2 0 msþmd ms form of the axion-fermion interaction Lagrangian, A A Fsd Fss a L ¼ −i f¯ ½ðm − m ÞcV aff i fi fj fifj The off-diagonal couplings in (20) and (21) induce kaon 2fa oscillations. The mixing matrix element M12 follows from A þðm þ m Þc γ5 f ; ð17Þ L fi fj fifj j a double insertion of the interaction Lagrangian aff ¼ −aq¯ðχS − iχPγ5Þq, where q ¼ u, d, s; cf. Eq. (16). which is obtained from Eq. (1) with the axion-dependent The Lagrangian for QCD with the flavor-violating field transformation L 3 axion, QCDþa, is formally invariant under SUð ÞR × 3 → χ a V − A a V A SUð ÞL transformations, qR;L gR;LðxÞqR;L,ifa S;P i2f ðcij cijÞ i2f ðcijþcijÞ fi → ½e a PL þ e a PR fj: ð18Þ and Mq are promoted to spurions that transform as
015023-7 JORGE MARTIN CAMALICH et al. PHYS. REV. D 102, 015023 (2020)
→ † s þ ip gRðs þ ipÞgL ð22Þ Expanding (24) in meson fields gives
2 a and that take the values Lð Þ ⊃− f m2 K¯ 0 ChPTþa FA K K ds M χ χ η s ¼ q þ Sa; p ¼ Pa: ð23Þ − iffiffiffi a 2 − 2 ¯ 0 π0 ffiffiffi p V ðmK mπÞK þ p 2 Fds 3 We also define the spurion χ ¼ χ þ iχ that does not 2 S P a 2m contain the axion field and that transforms similarly as þ K K0ðK¯ 0Þ2 þ H:c: ð26Þ † FA 3f χ → g χg . The identification of this symmetry structure ds K R L Here, we kept only the terms relevant for K-K¯ mixing and allows one to build the ChPT Lagrangian, including the pffiffiffi 2 155 6 chiral-symmetry breaking terms. The introduction of the replaced f with the kaon decay constant, fK ¼ . 0 4 spurion χ is needed because the axion cannot be treated . MeV [10], thus capturing part of the SU(3) breaking merely as an external field and enters in the chiral loops that corrects our results at higher orders in the ChPT expansion [139]. We also traded the products B0m for the with two insertions of χS;P. Thus, the ChPT Lagrangian is q also invariant under a Z2 symmetry that transforms meson masses squared [139,140,142]. The two terms in χ → −χ, with all the other fields taken to be Z2 even. (25), expanded in meson fields, are Up to overall normalization factors, the axion-induced ¯ 2 2 contributions to the K−K mixing amplitude have the Lð4Þ ⊃ 2 mK α 2α ν ν ChPTþa A ð 0 þ 1Þ scalings 2m M12 ∼ ðp=ΛχÞ ð1=FÞ F . The integer ν char- F K ds acterizes the usual chiral scaling [139,140] where the 2 − 2 2 mK mπ ¯ 0 2 derivatives of meson fields (and the axion) count as − ðα0 − 2α1Þ ðK Þ þ ; ð27Þ FV OðpÞ and the quark masses count as Oðp2Þ and where ds Λχ ≃ 4πfπ ∼ O 1 fπ the UV cutoff of ChPT is ð GeVÞ ( is where again we only keep the terms relevant for K-K¯ ν the pion decay constant). On the other hand, F counts the mixing. 1 V;A ¯ number of =Fds insertions in the amplitude. Thus, the The axion exchange contributions to the K-K mixing 2 chiral counting of the spurions is Mq ∼ Oðp Þ and νF ¼ 0 amplitude due to axial vector couplings are proportional χ ∼ O 2 ν 1 1 A 2 O 4 and S;P ðp Þ and F ¼ . In the following, we use to ð =FdsÞ and are, up to ðp Þ in the chiral counting, ChPT to calculate the leading axion-exchange contribu- given by ¯ tions to the K−K mixing amplitude including corrections ν ≤ 4 2 2 2 up to NLO corrections in the chiral counting, . This A fK mK mK M ¼ 1 − ðα0 þ 2α1Þ requires two insertions of 1=FV;A and thus ν ¼ 2. 12 FA 2 f2 ds F ds K The leading-order (LO) ChPT Lagrangian, including a 8 2 2 mK 1 − mK as the light degree of freedom, is given by þ 3 16π2 2 log μ2 : ð28Þ fK 2 ð2Þ f μ † L Tr ∂μU∂ U ChPTþa ¼ 4 ð Þ The first term in the parentheses is due to the tree-level 2 axion exchange, Fig. 1(a), and is induced by the first term f † þ B0 Tr½ðs − ipÞU þðs þ ipÞU in the expanded LO ChPT Lagrangian, Eq. (26). The NLO 2 correction is due to the loop diagram in Fig. 1(b), induced 2 1 μ ma 2 by the last term in Eq. (26). We use dimensional regulari- þ ∂μa∂ a − a ; ð24Þ 2 2 zation in the MS scheme and fix the renormalization scale to μ ¼ mK in order to minimize the size of the chiral while the relevant terms in the NLO ChPT Lagrangian are logarithm. The counterterms in Fig. 1(d) which cancel the μ ð4Þ 2 2 † † dependence of the loops are provided by the two terms in L ⊃ B f α0Tr χU Tr χ U 4 ChPTþa 0 ð ½ ½ the expanded Oðp Þ Lagrangian, Eq. (27). While the † 2 † 2 þ α1Tr½ðχU Þ þðχ UÞ Þ: ð25Þ numerical values of the low-energy constants α0;1 are not known, we can estimate their sizes by varying μ in Here, UðxÞ¼expðiλaπa=fÞ is the unitary matrix para- the loop contribution around its nominal value by a factor metrizing the meson fields [139,140]; B0 is a constant of 2, while artificially setting α0;1 to zero. This estimates the related to the quark condensate; B0ðμ ¼ 2 GeVÞ¼ Oðp4Þ contribution in (28) to be ð17% 23%Þ of the 2.666 57 GeV; f is related to the pion decay constant 2 ð pÞffiffiffi Oðp Þ one. f ≃ fπ= 2 ¼ 92.2ð1Þ MeV [141], with normalization The vector couplings of the axions first contribute at − 2 2 4 h0ju¯γμdð0Þjπ ðpÞi ¼ ipμfπ; and α0;1 ∼ f =Λχ are the Oðp Þ through the one loop diagram in Fig. 1(c) with the (unknown) low energy constants. counterterms in Fig. 1(d), resulting in
015023-8 QUARK FLAVOR PHENOMENOLOGY OF THE QCD AXION PHYS. REV. D 102, 015023 (2020)
(a) (b) (c) (d)
FIG. 1. Contributions to K0 − K¯ 0 mixing from exchanges of the flavor-violating axion, up to and including one loop in ChPT.