Minimal flavor-changing models and muon after the measurement
There has been a steady interest in flavor anomalies and their global fits as ideal probes of new physics. If the anomalies are real, one promising explanation is a new gauge boson with a flavor-changing coupling to bottom and strange quarks and a flavor-conserving coupling to muons and, possibly, electrons. We point out that direct production of such a , emerging from the collision of and quarks, may offer a complementary window into these phenomena because collider searches already provide competitive constraints. On top of that, we analyze the same scenario in relation to another long-standing discrepancy between theory and experiment that concerns the anomalous magnetic moment of the muon. By scanning the allowed coupling strengths in the low-mass region, we assess the compatibility of the signals from LHCb with the searches in the high energy LHC data and the measurements of the anomalous magnetic moments of the involved leptons. We also argue that observations of the latter can break the degeneracy pattern in the Wilson coefficients and presented by LHCb data. The model we consider is compatible with the new measurement of , therefore it can potentially account for the long-standing deviations observed in -physics.
Processes involving flavor-changing neutral currents are sensitive probes of new physics. In the Standard Model (SM), transitions such as are loop-suppressed, but new particles can contribute at tree-level. The possible impact of new particles on these processes is usually analyzed by integrating out the heavy degrees of freedom and working with the effective Hamiltonian. For transitions we have,
which is expressed in terms of the effective operators , and the Wilson coefficients , . Several anomalies with respect to SM predictions have been measured, typically at the level, but with increasing statistical significance. One notable anomaly of lepton flavor universality is the measured ratio of branching fractions Aaij:2014ora (), as well as the recent measurement of by LHCb LHCB (), which has already prompted several studies Capdevila:2017bsm (); Altmannshofer:2017yso (); DAmico:2017mtc (); Hiller:2017bzc (); Geng:2017svp (); Ciuchini:2017mik (); 1592397 (); 1592392 (). Angular observables have become a popular testing ground, and recent updates ATLAS:2017dlm (); CMS:2017ivg () have confirmed previous measurements in the decay . Including other transitions and a large number of observables, global fits to the anomalies seem to be converging on preferred sets of Wilson coefficients Altmannshofer:2017fio (); Descotes-Genon:2013wba (); Descotes-Genon:2015uva (). However, it is not yet entirely clear which of these sets could be responsible for the anomalies. For example, refinement of the hadronic uncertainty in these analyses is an ongoing theoretical issue.
The general features of a gauge boson Langacker:2008yv () make it a good choice to generate transitions, as underlined by many previous studies Descotes-Genon:2014uoa (); Crivellin:2015era (); Buras:2012jb (); Gauld:2013qba (); Buras:2013qja (); Gauld:2013qja (); Buras:2013dea (); Allanach:2015gkd (); Crivellin:2015mga (); Crivellin:2015lwa (); Belanger:2015nma (); Megias:2017ove (); Megias:2016bde (); Megias:2017dzd (); GarciaGarcia:2016nvr (); Carmona:2015ena (). In this work, we follow a complementary approach and consider what else can be gained from colliders and other low-energy experiments. Specifically, we examine the sensitivity of the Wilson coefficients to constraints from the direct production of a suitable gauge boson that can also explain the discrepancy between the measured value of the muon magnetic moment and the SM prediction.
We will focus here on scenarios in which the supposed new physics contributions affect exclusively and , for . These coefficients are among those currently favored by global fits to the anomalies Altmannshofer:2017fio (), as well as by a recent analysis Altmannshofer:2017yso (), which examines compatibility with and the new measurement of . The relevant operators are
Another motivation for this choice is that although the measurements of and suggest there is new physics which discriminates between muons and electrons, the LHCb data currently exhibit an interesting degeneracy in and based on their statistical pulls Altmannshofer:2017yso (): we argue that such degeneracy can be broken by analyses of the magnetic moments of the involved leptons in scenarios where the considered boson is well below the TeV scale. In addition, our analysis will show that the LEP bounds strongly disfavor scenarios where the speculated boson couples purely to electrons, besides quarks.
2 A simple model
We assume a general interaction of a boson with quarks and leptons described by the Lagrangian
where . The companion interactions with neutrinos and up-type quarks are allowed by the SM gauge symmetry and may offer interesting features which, however, we do not pursue in this work. The quark flavor-violating couplings , and the lepton flavor-conserving couplings , are normalized relative to from the Standard Model for convenience, being the cosine of the Weinberg angle.
After integrating out the and performing tree-level matching, the Wilson coefficients bounded by the LHCb results are related to the couplings by
The Wilson coefficient encapsulates the vectorial couplings between the and leptons, whereas contains the axial couplings. As we can see, despite its simplicity, the model at hand can potentially account for deviations from lepton flavor universality.
The contributions to the magnetic moment of a charged lepton are given by
for scenarios, and by
for scenarios, where is the mass of . In our analysis we will refer to the following values for the measured discrepancies of the involved lepton magnetic moments: , Giudice:2012ms (); Jegerlehner:2009ry ().
The parameters to be scanned over are and the set of couplings
|, , and ,||(7)|
which are for instance subject to the LHC dilepton searches that provide upper bounds on the production cross section times branching ratio as a function of , . Using the interaction Lagrangian above with specified input values of the mass and couplings, we numerically simulate the dilepton searches ATLAS:2016cyf () yielding the constraint
where is the result of our simulation, after the usual kinematical cuts, and is the experimental upper bound.
The most stringent limits on the couplings of the to electrons originate from LEP electroweak precision measurements LEP-2 (). Integrating out the in our model generates the effective four-fermion operators
which are constrained by measurements of cross sections at LEP. The upper limits on the magnitudes of the electron- couplings are
which respectively hold in the case of vectorial and axial couplings.
We scan the allowed range of couplings for two different masses of , with the values GeV, 500 GeV for scenarios involving muons, and GeV, 500 GeV for those involving electrons. Our simulations were performed using MadGraph_aMC5@NLO Alwall:2014hca () and we validated our code in the special case of a sequential against the simulated cross section reported in Ref. ATLAS:2016cyf () as a function of .
In our analyses we consider also the constraint from – oscillations in terms of the measured mass difference given by Altmannshofer:2013foa ()
where the contribution from the second term on the right-hand side is bounded by experiment to be in magnitude below the level. The SM loop function is .
As for the low-energy observations, including recent measurements by LHCb, we remark that a comprehensive statistical analysis of the Wilson coefficients would involve about 100 observables and a careful treatment of about 100 nuisance parameters. This is so technically challenging that instead approximate “fast-fit” techniques are favoured (see e.g., DAmico:2017mtc (); Altmannshofer:2017yso ()), although they require considerable computing resources. Since a comprehensive analysis is beyond the purpose of this paper, we adopt a pragmatic approach to understand the impact of the latest data. We perform an independent fit of the Wilson coefficients to LHCb measurements of LHCB () and Aaij:2014ora (), Belle measurements of and Wehle:2016yoi (), and LHCb measurements of and Aaij:2017vad (). Recent studies DAmico:2017mtc (); Geng:2017svp () performed global frequentist statistical analysis, whereas Ciuchini:2017mik () performed a global Bayesian analysis with a subset of nuisance parameters and data with HEPfit HEPfit (). Our methodology and its results are detailed in the next section.
3 A first analysis
We calculated the profile likelihood with MultiNest Feroz:2008xx () interfaced with a modified version of flavio Altmannshofer:2017fio (). Confidence intervals were found by Wilks’ theorem (see Fowlie:2016hew () for conventions and definitions). All nuisance parameters (, , CKM matrix elements and form factors) were fixed to their central values. Wilson coefficients were varied between and . We find reasonable agreement for best-fits, confidence intervals and significances with similar flavio results with Markov Chain Monte Carlo algorithms and fast-fit techniques for nuisance parameters Altmannshofer:2017yso (). We estimate the significance of the preference for new physics versus the SM to be about , in agreement with recent literature DAmico:2017mtc (); Geng:2017svp (); Ciuchini:2017mik (); Altmannshofer:2017yso (); Capdevila:2017bsm (), but, as in the literature, there are important caveats about systematic uncertainties.
Since it is a product of our MultiNest fit, we note that data favour the model by a Bayes factor of about versus the SM and about versus , as well as versus or . Since, however, we omitted important nuisance parameters that may alleviate tension in the SM, should be regarded as an upper bound to the Bayes factor versus the SM (and is, of course, sensitive to priors). A Bayes factor preference for versus electron Wilson coefficients was also noted in an earlier study of Ghosh:2014awa ().
The results we obtained are presented in Figure 1, whereas Table 1 proposes the results from Ref. Altmannshofer:2017yso () for comparison. We find that the two analyses are in reasonable agreement given that the fits take into account different sets of observables.
We focus now on the constraints that collider and low energy experiments impose on the scenarios we delineated.
|[, ]||[, ]|
|[, ]||[, ]|
|[, ]||[, ]|
|[, ]||[, ]|
3.1 Scenarios with or only
We present in Figure 2 the results obtained for the scenario where new physics effects arise from the coupling of to electrons, and the relative effects are fully encapsulated in the Wilson coefficient (top panels), or (bottom panels). The different shaded areas indicate allowed regions of parameter space after taking into account the various constraints.
The green area represents the region allowed by the ATLAS searches in the dielectron channel. The blue band represents the values of the coupling selected by the LHCb measurements of the indicated Wilson coefficient. We refer here to the values for the credible regions from our analysis presented in Figure 1. The red band illustrates instead the region of the parameter space allowed by the constraints on the mass difference .
As mentioned before, we also checked that the contribution to the electron does not spoil the current agreement between theory and experiment. A plot of the new contribution to this quantity is presented in Figure 3 as a function of for different values of the relevant coupling. The results in the top panel holds for the scenario where all the new physics effects are contained in , the bottom one for the case of . We can see that values of GeV would negatively impact on the prediction for the anomalous magnetic moment of the electron.
We remark that these scenarios where the only couples to electrons, besides quarks, are strongly constrained by the bounds from LEP-II LEP-2 (). The gray shade in Figure 2 denotes the area of the parameter space that evades the latter. We can see that the scenario with non-zero generated from a vector coupling is more tightly constrained than by LEP electroweak precision measurements and – oscillations. Indeed in the case of non-zero a larger region of the parameter space is still allowed by all the considered experimental limits. A boson with axial vector couplings would also induce a negative contribution to the of the electron, and is therefore amenable to reducing the tension between the measurements and the SM prediction. Although in the simple models we are considering this contribution is not substantial, we remark that the preference for based on is consistent with the picture emerging from the other constraints analysed.
3.2 Scenarios with and only
We present in Figure 4 the analogous results obtained for the case in which, besides quarks, the speculated couples exclusively to muons.
Again we show in the top panel the case where new physics affects only, and the bottom panels cover the complementary case of . As before, the green area represents the region of the parameter space allowed by the ATLAS searches, here for the dimuon channel. The blue band represents the values of the coupling selected by the LHCb measurements of the indicated Wilson coefficient, according to the credible intervals shown in Figure 1. The red band is the region of the parameter space allowed by the constraints on the mass difference . We remark that LEP measurements do not constrain the couplings in these scenarios.
The contribution in the scenario can help to bridge the discrepancy between the measured value of the muon magnetic moment and the SM prediction. As we can see in Figure 5, the contribution is potentially able to reduce such a discrepancy below the level for moderate values of the -muon couplings for masses up to about GeV. The solid black line in this region with and GeV is compatibile with the parameter space selected by the LHCb data. In this case, the physical coupling is large, but marginally within the perturbative regime. The corresponding region in Figure 4 is allowed, when the top left panel is extrapolated to .
As for the case (not shown), the presence of a in the low-mass range can only worsen the theory prediction as the relative contribution, of negative sign, further lowers the value of the latter.
In this paper we considered four different new physics scenarios which could explain the anomalous -physics results from LHCb and argued that measurements of the anomalous magnetic moments of muon and electron could break the degeneracy in the – Wilson coefficients.
We find that a -boson which couples only to electrons can produce the correct values for and , but such a scenario is strongly constrained by the precision measurements from LEP. However, if the coupling to is purely axial, the allowed region of parameter space is larger. Such an axially coupled also generates a negative contribution to the anomalous magnetic moment of the electron and could therefore potentially accommodate both anomalies.
We also considered scenarios where the -boson couples only to muons. In this case, a vectorial coupling to the is favoured by the Bayesian approach employed for the performed fit. Such a scenario produces also a positive contribution to the of the muon and therefore could help to alleviate the tension between the measurement and the SM prediction. This scenario can also explain the LHCb anomalies while avoiding all the remaining constraints from high energy searches, provided the vectorial coupling to muons is large (). An axial coupling of the to muons can also explain the LHCb results, but increases the disagreement of the with measurements.
We thank Christian Veelken for useful discussions. This work was supported by the EU through the ERDF CoE program grant TK133. The work of CM is supported by the “Angelo Della Riccia” foundation.
- (1) LHCb collaboration, R. Aaij et al., Test of lepton universality using decays, Phys. Rev. Lett. 113 (2014) 151601, [1406.6482].
- (2) CMS collaboration, S. Bifani, Search for new physics with decays at LHCb, CERN Seminar, 18 April 2017 .
- (3) B. Capdevila, A. Crivellin, S. Descotes-Genon, J. Matias and J. Virto, Patterns of New Physics in transitions in the light of recent data, 1704.05340.
- (4) W. Altmannshofer, P. Stangl and D. M. Straub, Interpreting Hints for Lepton Flavor Universality Violation, 1704.05435.
- (5) G. D’Amico, M. Nardecchia, P. Panci, F. Sannino, A. Strumia, R. Torre et al., Flavour anomalies after the measurement, 1704.05438.
- (6) G. Hiller and I. Nisandzic, and beyond the Standard Model, 1704.05444.
- (7) L.-S. Geng, B. Grinstein, S. Jäger, J. Martin Camalich, X.-L. Ren and R.-X. Shi, Towards the discovery of new physics with lepton-universality ratios of decays, 1704.05446.
- (8) M. Ciuchini, A. M. Coutinho, M. Fedele, E. Franco, A. Paul, L. Silvestrini et al., On Flavourful Easter eggs for New Physics hunger and Lepton Flavour Universality violation, 1704.05447.
- (9) D. Bečirević and O. Sumensari, A leptoquark model to accommodate and , 1704.05835.
- (10) A. Celis, J. Fuentes-Martin, A. Vicente and J. Virto, Gauge-invariant implications of the LHCb measurements on Lepton-Flavour Non-Universality, 1704.05672.
- (11) ATLAS collaboration, Angular analysis of decays in collisions at TeV with the ATLAS detector, Tech. Rep. ATLAS-CONF-2017-023, CERN, Geneva, Apr, 2017.
- (12) CMS collaboration, Measurement of the and angular parameters of the decay in proton-proton collisions at , Tech. Rep. CMS-PAS-BPH-15-008, CERN, Geneva, 2017.
- (13) W. Altmannshofer, C. Niehoff, P. Stangl and D. M. Straub, Status of the anomaly after Moriond 2017, 1703.09189.
- (14) S. Descotes-Genon, J. Matias and J. Virto, Understanding the Anomaly, Phys. Rev. D88 (2013) 074002, [1307.5683].
- (15) S. Descotes-Genon, L. Hofer, J. Matias and J. Virto, Global analysis of anomalies, JHEP 06 (2016) 092, [1510.04239].
- (16) P. Langacker, The Physics of Heavy Gauge Bosons, Rev. Mod. Phys. 81 (2009) 1199–1228, [0801.1345].
- (17) S. Descotes-Genon, L. Hofer, J. Matias and J. Virto, On the impact of power corrections in the prediction of observables, JHEP 12 (2014) 125, [1407.8526].
- (18) A. Crivellin, L. Hofer, J. Matias, U. Nierste, S. Pokorski and J. Rosiek, Lepton-flavour violating decays in generic models, Phys. Rev. D92 (2015) 054013, [1504.07928].
- (19) A. J. Buras, F. De Fazio and J. Girrbach, The Anatomy of Z’ and Z with Flavour Changing Neutral Currents in the Flavour Precision Era, JHEP 02 (2013) 116, [1211.1896].
- (20) R. Gauld, F. Goertz and U. Haisch, On minimal explanations of the anomaly, Phys. Rev. D89 (2014) 015005, [1308.1959].
- (21) A. J. Buras and J. Girrbach, Left-handed and FCNC quark couplings facing new data, JHEP 12 (2013) 009, [1309.2466].
- (22) R. Gauld, F. Goertz and U. Haisch, An explicit -boson explanation of the anomaly, JHEP 01 (2014) 069, [1310.1082].
- (23) A. J. Buras, F. De Fazio and J. Girrbach, 331 models facing new data, JHEP 02 (2014) 112, [1311.6729].
- (24) B. Allanach, F. S. Queiroz, A. Strumia and S. Sun, models for the LHCb and muon anomalies, Phys. Rev. D93 (2016) 055045, [1511.07447].
- (25) A. Crivellin, G. D’Ambrosio and J. Heeck, Explaining , and in a two-Higgs-doublet model with gauged , Phys. Rev. Lett. 114 (2015) 151801, [1501.00993].
- (26) A. Crivellin, G. D’Ambrosio and J. Heeck, Addressing the LHC flavor anomalies with horizontal gauge symmetries, Phys. Rev. D91 (2015) 075006, [1503.03477].
- (27) G. Bélanger, C. Delaunay and S. Westhoff, A Dark Matter Relic From Muon Anomalies, Phys. Rev. D92 (2015) 055021, [1507.06660].
- (28) E. Megias, M. Quiros and L. Salas, Lepton-flavor universality violation in and from warped space, 1703.06019.
- (29) E. Megias, G. Panico, O. Pujolas and M. Quiros, A Natural origin for the LHCb anomalies, JHEP 09 (2016) 118, [1608.02362].
- (30) E. Megias, M. Quiros and L. Salas, from Vector-Like Leptons in Warped Space, 1701.05072.
- (31) I. Garcia Garcia, LHCb anomalies from a natural perspective, JHEP 03 (2017) 040, [1611.03507].
- (32) A. Carmona and F. Goertz, Lepton Flavor and Nonuniversality from Minimal Composite Higgs Setups, Phys. Rev. Lett. 116 (2016) 251801, [1510.07658].
- (33) G. F. Giudice, P. Paradisi and M. Passera, Testing new physics with the electron g-2, JHEP 11 (2012) 113, [1208.6583].
- (34) F. Jegerlehner and A. Nyffeler, The Muon g-2, Phys. Rept. 477 (2009) 1–110, [0902.3360].
- (35) ATLAS collaboration, Search for new high-mass resonances in the dilepton final state using proton-proton collisions at = 13 TeV with the ATLAS detector, Tech. Rep. ATLAS-CONF-2016-045, CERN, Geneva, Aug, 2016.
- (36) DELPHI, OPAL, LEP Electroweak, ALEPH, L3 collaboration, S. Schael et al., Electroweak Measurements in Electron-Positron Collisions at W-Boson-Pair Energies at LEP, Phys. Rept. 532 (2013) 119–244, [1302.3415].
- (37) J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer et al., The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations, JHEP 07 (2014) 079, [1405.0301].
- (38) W. Altmannshofer and D. M. Straub, New Physics in ?, Eur. Phys. J. C73 (2013) 2646, [1308.1501].
- (39) Belle collaboration, S. Wehle et al., Lepton-Flavor-Dependent Angular Analysis of , Phys. Rev. Lett. 118 (2017) 111801, [1612.05014].
- (40) LHCb collaboration, R. Aaij et al., Measurement of the branching fraction and effective lifetime and search for decays, 1703.05747.
- (41) J. de Blas et al., HEPfit: a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models, .
- (42) F. Feroz, M. P. Hobson and M. Bridges, MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics, Mon. Not. Roy. Astron. Soc. 398 (2009) 1601–1614, [0809.3437].
- (43) A. Fowlie and M. H. Bardsley, Superplot: a graphical interface for plotting and analysing MultiNest output, Eur. Phys. J. Plus 131 (2016) 391, [1603.00555].
- (44) D. Ghosh, M. Nardecchia and S. A. Renner, Hint of Lepton Flavour Non-Universality in Meson Decays, JHEP 12 (2014) 131, [1408.4097].