Review Article  Open Access
Anomalies in Transitions and Dark Matter
Abstract
Since 2013, the LHCb collaboration has reported on the measurement of several observables associated with transitions, finding various deviations from their predicted values in the Standard Model. These include a set of deviations in branching ratios and angular observables, as well as in the observables and , specially built to test the possible violation of Lepton Flavor Universality. Even though these tantalizing hints are not conclusive yet, the anomalies have gained considerable attention in the flavor community. Here we review new physics models that address these anomalies and explore their possible connection to the dark matter of the Universe. After discussing some of the ideas introduced in these works and classifying the proposed models, two selected examples are presented in detail in order to illustrate the potential interplay between these two areas of current particle physics.
1. Introduction
The Standard Model (SM) of particle physics provides an excellent description for a vast amount of phenomena and can be regarded as one of the most successful scientific theories ever built. In fact, with the recent discovery of the last missing piece, the Higgs boson, the particle spectrum is finally complete and the SM looks stronger than ever. However, and despite its enormous success, there are several indications that clearly point towards the existence of a more complete theory, with neutrino masses and the baryon asymmetry of the Universe as the most prominent examples.
Another open question is the nature of the dark matter (DM) that accounts for of the energy density of the Universe [1]. Several ideas have been proposed to address this fundamental problem in current physics. Under the hypothesis that the DM is composed of particles, these cannot be identified with any of the states in the SM, hence demanding an extension of the model with new states and, possibly, new dynamics. Again, many directions exist. Interestingly, in scenarios involving New Physics (NP) at the TeV scale, the first signals from the new DM sector might be found in experiments not specially designed to look for them.
Rare decays stand among the most powerful tests of the SM. Since 2013, results obtained by the LHCb collaboration have led to an increasing interest in B physics, particularly in processes involving transitions. Deviations from the SM expectations have been reported in several observables, some of them hinting at the violation of Lepton Flavor Universality (LFU), a central feature in the SM. Even though these anomalies could be caused by a combination of unfortunate fluctuations and, perhaps, a poor theoretical understanding of some processes, it is tempting to speculate about their possible origin in terms of NP models, in particular models linking them to other open problems.
This minireview will pursue this goal, focusing on NP scenarios that relate the anomalies to DM. Several works [2–20] have already explored this connection, mostly by means of specific models that accommodate the observations in transitions with a new dark sector. We will review some of the ideas introduced in these works and highlight those that deserve further exploration. We will also classify the proposed models into two general categories: (i) models in which the NP contributions to transitions and DM production in the early Universe share a common mediator, and (ii) models with the DM particle running in loop diagrams that contribute to the solution of the anomalies. After a general discussion, a selected example of each class will be presented in detail.
The rest of the manuscript is organized as follows. First, we review the anomalies in transitions in Section 2 and interpret the experimental results in a model independent way in Section 3. In Section 4 we discuss and classify the proposed New Physics explanations to these anomalies that involve a link to the dark matter problem. Sections 5 and 6 present two simple example models that illustrate this connection. Finally, we summarize and draw our conclusions in Section 7.
2. Experimental Situation
We begin by discussing the present experimental situation. The observed anomalies in transitions can be classified into two classes: (1) branching ratios and angular observables and (2) lepton flavor universality violating (LFUV) anomalies. Although they might be related (and caused by the same NP), they are conceptually different.
Branching Ratios and Angular Observables. using stateoftheart computations of the hadronic form factors involved, one can compute branching ratios and angular observables for processes such as and look for deviations from the SM predictions. For the comparison to be meaningful, one must have a good knowledge of all possible Quantum ChromoDynamics (QCD) effects that might pollute the theoretical calculation and we currently have at our disposal several methods to minimize or at least estimate the uncertainties. (The size of the hadronic uncertainties in different calculations is a matter of hot debate nowadays. We will not discuss this issue here but just refer to the recent studies regarding form factors [21, 22] and nonlocal contributions [23–28] for extended discussions.) In particular, a basis of optimized observables for the decay , specially designed to reduce the hadronic uncertainties, was introduced in [29]. In 2013, the LHCb collaboration published results on these observables using their fb^{−1} dataset, finding a deviation between the measurement and the SM prediction for the angular observable in one dimuon invariant mass bin [30]. A systematic deficit with respect to the SM predictions in the branching ratios of several processes, mainly , was also reported by LHCb [31]. These discrepancies have been found later in other datasets. In 2015, LHCb confirmed these anomalies using their full Run 1 dataset with fb^{−1} [32, 33], whereas in 2016 the Belle collaboration presented an independent measurement of , compatible with the LHCb result [34, 35]. More recently, both ATLAS [36] and CMS [37] have also presented preliminary results on the angular observables, with relatively good agreement with LHCb.
LFUV Anomalies. One of the central features of the SM is that gauge bosons coupled with the same strength to all three families of leptons. This prediction can be tested by measuring observables such as the ratios, defined as [38]measured in specific dilepton invariant mass squared ranges . In the SM, these ratios should be very approximately equal to one. Furthermore, hadronic uncertainties are expected to cancel very good approximation in these ratios, which implies that, in contrast to the previous class of anomalies, deviations in these observables cannot be explained by uncontrolled QCD effects and would be a clear indication of NP at work. For this reason, they are sometimes referred to as clean observables. Interestingly, in 2014, the LHCb collaboration measured in the region GeV^{2} [39], finding a value significantly lower than one, while in 2017 similar measurements of the ratio in two bins [40] were also found to depart from their SM expected values:The comparison between these experimental results and the SM predictions [41, 42],shows deviations from the SM at the level in the case of , for in the low region, and for in the central region. Finally, Belle has recently measured the LFUV observable , with the observable defined for analogously to for [43]. The result, although statistically not very significant, also points towards the violation of LFU [35].
Summarizing, there are at present two sets of experimental anomalies in processes involving transitions at the quark level. While the relevance of the first set is currently a matter of discussion due to the possibility of unknown QCD effects faking the deviations from the SM, the second can only be explained by NP violating LFU. In principle, these two classes of anomalies can be completely unrelated but, as we will see in the next section, global analyses of all experimental data in transitions indicate that a common explanation (in terms of a single effective operator) can address both sets in a satisfactory and economical way. This intriguing result has made the anomalies a topic of great interest currently.
Finally, it is very interesting to note the existence of an independent set of anomalies in transitions. Several experimental measurements of the ratios and have been found to depart from their SM predictions, with a global discrepancy at the level [44]. Recently, the ratio has also been measured by the LHCb collaboration, finding again a deviation from the SM expected value [45]. Compared to the anomalies, the anomalies are of a different nature and, if real, they could have a completely different origin. For instance, they would involve a new charged current, instead of a neutral one, hence requiring the new mediators to be much lighter to be able to compete with the SM boson treelevel exchange. However, many authors have proposed models that can simultaneously address both sets of anomalies. We refer to [46] for a general discussion on combined explanations and ignore the anomalies for the rest of this paper.
3. Model Independent Interpretation
The experimental tensions discussed in the previous section must be properly quantified and interpreted. Quantification is crucial to determine whether the anomalies can be explained by fluctuations in the data or they truly indicate a statistically significant deviation from the SM. Assuming that these tensions are caused by genuine NP, the ultimate goal is to construct a specific model in which they are solved. However, the first step in this direction must be a model independent interpretation of the experimental data in order to identify the ingredients that this new scenario must include. This is achieved by adopting an approach based on effective operators, valid under the assumption that all NP degrees of freedom lie at energies well above the relevant energy scales for the observables of interest.
The effective Hamiltonian for transitions is usually written asHere is the Fermi constant, the electric charge, and the CabibboKobayashiMaskawa (CKM) matrix. and are the effective operators that contribute to transitions, and and their Wilson coefficients. The most relevant operators for the interpretation of the anomalies areHere . In fact, the operators and Wilson coefficients carry flavor indices and we are omitting them to simplify the notation. When necessary, we will denote a particular lepton flavor with a superscript, e.g., and , for muons. It is also convenient to split the Wilson coefficients in two pieces: the SM contributions and the NP contributions, defining the following (similar splittings could be defined for the Wilson coefficients of the primed operators, and , but in this case the SM contributions are suppressed and one has and ):The SM contributions have been computed at NNLO at GeV, obtaining and (see [29] and references therein), leaving the NP contributions as parameters to be determined (or at least constrained) by using experimental data.
It is in principle possible to derive limits for the NP contributions considering each observable independently, but this approach would completely miss the global picture. The effective operators in (4) contribute to several observables and one expects the presence of NP to be revealed by a pattern of deviations from the SM expectations, rather than by a single anomaly. For this reason, global fits constitute the best approach to analyze the available experimental data. Interestingly, several independent fits [47–54] have found a remarkable tension between the SM and experimental data on transitions which is clearly reduced with the addition of NP contributions. Although the numerical details (such as statistical significance) differ among different analyses, there is a general consensus on the following qualitative results:(i)Global fits improve substantially with a negative contribution in , with , leading to a total Wilson coefficient significantly smaller than the one in the SM.(ii)NP contributions in other Wilson coefficients can also improve the fit, but only in a subdominant way. For instance, the anomalies can also be accommodated in scenarios with or , without a clear statistical preference with respect to the scenario with NP only in . (Such patterns for the Wilson coefficients are automatically obtained if the NP states coupled to SM fermions with specific chiralities. For instance, the relation is obtained in models where the NP states only couple to the lefthanded muons. Two examples of this class of models are shown in Sections 5 and 6.)(iii)Other operators involving muons are perfectly compatible with their SM values. Similarly, no NP is required for operators involving electrons or tau leptons.
Armed with these results, model builders can construct specific models where all requirements are met and the anomalies explained. Similarly, one can extract interesting implications for model building by explaining the anomalies in terms of gauge invariant effective operators; see [55] for a recent analysis. Either way, the resulting profile of NP contributions reveals a pattern that was not predicted by any theoretical framework, such as supersymmetry, and many new models have been put forward. In the next Section we will discuss some of these models, in particular those linking the anomalies to dark matter.
4. Linking the Anomalies to Dark Matter
After discussing the current experimental situation in transitions, let us focus on possible connections to the dark matter problem. These have been explored in [2–20]. In general, the proposed models that solve the anomalies and explain the origin of the dark matter of the Universe can be classified into two principal categories:(i)Portal models: models in which the mediator responsible for the NP contributions to transitions also mediates the DM production in the early Universe.(ii)Loop models: models that induce the required NP contributions to transitions with loops containing the DM particle.
In the case of portal models, the usual scenario considers a gauge extension of the SM that leads to the existence of a new massive gauge boson after spontaneous symmetry breaking. The resulting boson induces a new neutral current contribution in transitions and mediates the production of DM particles in the early Universe via a portal interaction. This setup was first considered in [2]. In this particular realization of the general idea, the SM fermions were assumed to be neutral under the new gauge symmetry and the couplings to quarks () and leptons (), necessary to explain the anomalies, are generated at treelevel via mixing with new vectorlike (VL) fermions. Additionally, the boson also couples to the scalar field , the DM candidate in this model, automatically stabilized by an remnant symmetry after breaking. This model will be reviewed in more detail in Section 5. Variants of this setup with fermionic DM also exist. In [6], a horizontal gauge symmetry is introduced, with being the baryon number of the ith fermion family. The resulting boson couples directly to the SM quarks, while the coupling to muons is obtained by introducing a VL lepton. This allows accommodating the anomalies in transitions. Furthermore, the model also contains a Dirac fermion that is stable due to a remnant symmetry, in a similar fashion as in [2], which becomes the DM candidate. Similarly, [7] builds on the wellknown model of [56] and extends it to include a stable Dirac fermion with a relic density also determined by portal interactions, while [20] considers a similar model but makes use of kinetic mixing between the and the SM neutral gauge bosons. Reference [19] considers vectorlike neutrino DM in a setup analogous to [2] extended with additional VL fermions. Reference [10] explored a pair of scenarios based on a gauge symmetry supplemented with VL fermions and a fermionic DM candidate, of Dirac or Majorana nature. This paper focuses on effects in indirect detection experiments, aiming at an explanation of the excess of events in antiproton spectra reported by the AMS02 experiment in 2016 [57]. Other works that adopt the standard portal setup are [4, 12]. Finally, [11] considers a light mediator (not the usual heavy ) that contributes to transitions and decays predominantly into invisible final states, possibly made of light DM particles.
It is important to note that the phenomenology of these portal models differs substantially from the standard portal phenomenology. This is due to the fact that the bosons in these models couple with different strengths to different fermion families, as required to accommodate the LFUV hints observed by the LHCb collaboration ( and ). For instance, DM annihilation typically yields muon and tau lepton pairs, but not electrons and positrons. Direct detection experiments are also more challenging in the standard portal scenario, since the DM candidate typically does not couple to first generation quarks, more abundant in the nucleons.
In what concerns loop models, many variations are possible. To the best of our knowledge, the first model of this type that appeared in the literature is [13], based on previous work on loop models for the anomalies, without connecting to the DM problem, in [58, 59]. In this model the SM particle content is extended with two VL pairs of doublets, with the same quantum numbers as the SM quark and lepton doublets, but charged under a global symmetry. The model also contains the complex scalar , singlet under the SM gauge symmetry and also charged under . With these states, one can draw a 1loop diagram contributing to the observables relevant to explaining the anomalies. Furthermore, if the global is conserved, the lightest charged state becomes stable. In this work, this state is assumed to be , hence the DM candidate in the model. A more detailed discussion about this model can be found in Section 6. Two similar setups can be found in [18], where a different set of global symmetries are considered ( and ) in order to stabilize a scalar DM candidate. This paper also includes righthanded neutrinos in order to accommodate nonzero neutrino masses with the typeI seesaw mechanism and explores the lepton flavor violating phenomenology of the model in detail. A Majorana fermionic DM candidate was considered in [16]. Similarly to the previously mentioned models, this scenario also addresses the anomalies at 1loop level introducing a minimal number of fields: just a VL quark () and an inert scalar doublet (), in addition to the fermion singlet that constitutes the DM candidate. The model is supplemented with a discrete symmetry to ensure the stability of the DM particle. Interestingly, the model can be tested in direct DM detection experiments as well as at the Large Hadron Collider (LHC), where the states and can be pairproduced and lead to final states with hard leptons and missing energy. Finally, an extended loop model for the anomalies, which also has an additional gauge symmetry, contains a scalar DM candidate and explains neutrino masses that can be found in [17].
Finally, let us comment on other models and works that do not easily fit within any of the two categories mentioned above. The model in [3] is very similar to the model in [2]. It also extends the SM with a complex scalar, VL quarks, and leptons and a new gauge symmetry that breaks down to a parity. However, the VL leptons carry different charges, leading to a loopinduced coupling. This changes the DM phenomenology dramatically. The dominant mechanism for the DM production in the early Universe is not a portal interaction, but tchannel exchange of VL leptons. The model in [9] can be regarded as a hybrid model, with features from both portal and loop models. The SM symmetry group is extended with a new piece. The first factor leads to the existence of a massive boson while the second one stabilizes a scalar DM candidate. The coupling is generated with a loop containing the odd fields and the dominant DM production mechanism is a portal interaction, mainly with leptons. The coupling is also loopgenerated in [14], but in this case production of DM particles takes place via a Higgs portal. Reference [8] proposes an extended Scotogenic model for neutrino masses [60] supplemented with a nonuniversal gauge group. The DM candidate in this case is the lightest fermion singlet and is produced by Yukawa interactions. Finally, two models that address the anomalies with leptoquarks and include DM candidates were introduced in [5, 15]. In the former the DM candidate is a component of an multiplet introduced to enhance the diphoton rate of a scalar in the model, whereas in the latter the DM candidate is a baryonlike composite state in a model with strong dynamics at the TeV scale.
Having reviewed and classified the proposed models, we now proceed to discuss in some detail two specific examples. These illustrate the main features of portal and loop models.
5. An Example Portal Model
We will now review the model introduced in [2], arguably one of the simplest scenarios to account for the anomalies with a dark sector. Some of the ingredients of this model were already present in the model of [56], which is extended in the quark sector (following the same lines as in the lepton sector). It also includes a dark matter candidate that couples to the SM fields via the same mediator that leads to an explanation of the anomalies, a heavy boson. A variation of this scenario with a loopinduced coupling to muons appeared afterwards in [3], whereas the phenomenology of an extension to account for neutrino masses will be discussed in [61].
The model extends the SM gauge group with a new dark factor, under which all the SM particles are assumed to be singlets. The dark sector contains two pairs of vectorlike fermions, and , as well as the complex scalar fields, and . Tables 1 and 2 show all the details about the gauge sector and the new scalars and fermions in the model. and have the same representation under the SM gauge group as the SM doublets and , and they can be decomposed under asWith the electric charges of , , , and being , , , and , respectively. In contrast to their SM counterparts, and are vectorlike fermions charged under the dark . In addition to the usual canonical kinetic terms, the new vectorlike fermions and have Dirac mass terms,as well as Yukawa couplings with the SM doublets and and the scalar ,where and are 3 component vectors. The scalar potential of the model can be split into different pieces, is the usual SM scalar potential containing quadratic and quartic terms for the Higgs doublet . The new terms involving the scalars and areandAll couplings are dimensionless, whereas has dimensions of mass and and have dimensions of mass^{2}. We will assume that the scalar potential parameters allow for the vacuum configurationwhere is the neutral component of the Higgs doublet . The scalar does not get a vacuum expectation value (VEV). Therefore, the scalar will be responsible for the spontaneous breaking of . This automatically leads to the existence of a new massive gauge boson, the boson, with mass . In the absence of mixing between the gauge bosons, the boson can be identified with the original boson in Table 1. We note that a Lagrangian term of the form , where are the usual field strength tensors for the groups, would induce this mixing. In order to avoid phenomenological difficulties associated with this mixing we will assume that . Moreover, it can be shown that loop contributions to this mixing are kept under control if [2].


Let us now discuss how this model solves the anomalies. After the spontaneous breaking of , the Yukawa interactions in (11) lead to mixings between the vectorlike fermions and their SM counterparts. This mixing results in effective couplings to the SM fermions. If these are parametrized as [62, 63]and one assumes for the sake of simplicity, the couplings to and , necessary to solve the anomalies, are found to beThese couplings induce a treelevel contribution to the semileptonic fourfermion operators in (5) and (6), as shown schematically in Figure 1. More specifically, given that the SM fermions participating in the effective vertices are lefthanded, see (11), the operators and are generated simultaneously, with [63]where we have introducedwith being the electromagnetic fine structure constant. and the CKM elements appear in (18) in order to normalize the Wilson coefficients as defined in (5) and (6). By taking proper ranges for the model parameters, the required values for these Wilson coefficients, previously identified by the global fits to flavor data, can be easily obtained.
Finally, we move to the discussion on the Dark Matter phenomenology of the model. We note that the model does not include any ad hoc stabilizing symmetry for the DM candidate . However, this state is perfectly stable. This is due to the fact that the continuous symmetry leaves a remnant parity, under which is odd, after spontaneous symmetry breaking [64–66]. Therefore, the same symmetry that leads to the dynamics behind the anomalies is also at the origin of the DM stabilization mechanism. Furthermore, the DM production in the early Universe can take place via processes mediated by the boson, thus establishing another link with the anomalies. Indeed, purely gauge interactions open a portal that induces annihilation processes, with , as shown in Figure 2. (Another possibility is the socalled Higgs portal, activated in this model with the scalar potential term , which induces processes. This DM production mechanism will be subdominant for sufficiently small ). We notice that these processes match those in Figure 1 if one trades one of the fermion pairs for . Therefore, one can establish an interplay between flavor and DM physics in this scenario. Figure 3 illustrates this connection displaying contours for constant (the DM relic density) and the ratio in the plane. This figure has been obtained with fixed , , , and . The calculation of the flavor observables has been performed with FlavorKit [67], whereas the DM relic density has been evaluated with MicrOmegas [68]. We see that there is a region in parameter space, with moderately large , where the observed DM relic density can be reproduced and a ratio in agreement with the global fits is obtained. We also note that the DM relic density tends to be large. In fact, in order to obtain a numerical value in the ballpark of one has to be rather close to the resonant region with , which in this plot is located around TeV.
Besides flavor and DM physics, the model has rich phenomenological prospects in other fronts. The new states can be discovered at the LHC in large portions of the parameter space. Although one typically assumes that the boson couples predominantly to the second and thirdgeneration quarks (), the resulting suppressed production crosssections at the LHC can still be sufficient for a discovery; see, for instance, [69]. Furthermore, the new VL fermions can also be produced and detected. In particular, the heavy VL quarks masses are already pushed beyond the TeV scale due to their efficient production in collisions. In what concerns direct and indirect DM detection, scenarios with a dark portal have been discussed in [70, 71]. For more details about this model, its predictions, and the most relevant experimental constraints we refer to [2].
6. An Example Loop Model
We now turn our attention to the second class of models, those that explain the anomalies via loop diagrams including DM particles. A simple but illustrative example of this category is that presented in [13]. Previous work on loop models for the anomalies, without connecting to the DM problem, can be found in [58, 59].
The model introduces two VL fermions, and , with the same gauge quantum numbers as the SM quark and lepton doublets, respectively. It also adds the complex scalar , a complete singlet under the SM gauge symmetry. The new fields are charged under a global Abelian symmetry, , under which all SM fields are assumed to be singlets. As we will see below, this particle content is sufficient to address the anomalies. Table 3 details the new fields and their charges under the gauge and global symmetries of the model.

The VL fermions and can be decomposed as in (9), with their components having exactly the same electric charges as in that case. Therefore, the same Dirac mass terms as in (10) can be written. In addition, the symmetries of the model allow for the Yukawa couplings with the SM doublets and and the scalar where and are 3 component vectors. The scalar potential of the model contains the following terms:All couplings are dimensionless, whereas has dimensions of mass^{2}. In the following, possible effects due to the coupling will be ignored, assuming . We will also assume that the scalar potential parameters allow for a vacuum configuration with . In this case, the global symmetry is conserved and the lightest state with a nonvanishing charge under this symmetry is completely stable. Moreover, we note that the conservation of prevents the VL fermions from mixing with the SM ones.
We move now to the solution of the anomalies in the context of this model. It is straightforward to check that no NP contributions to transitions are generated at treelevel in this model. (For instance, in contrast to the model discussed in Section 5, there is no SMVL mixing, nor a boson that can mediate these transitions at treelevel.) However, the semileptonic operators and are generated at the 1loop level as shown in Figure 4. This diagram leads towhere was introduced in (19) and is the loop function
This looplevel solution to the anomalies corresponds to scenario AI, model class b), in [59]. Figure 5 shows that the model can accommodate the and measurements by the LHCb collaboration. This figure has been obtained with fixed , TeV, and TeV. One finds that the and regions for and almost overlap, and thus they can be accommodated in the same region of parameter space. Furthermore, in order to be compatible with the bounds coming from mixing one needs , which implies a relatively large value of , . This feature, a hierarchy between the NP couplings to quarks and leptons, is shared by most models addressing the anomalies. For a general discussion about the mixing constraint in the context of the anomalies we refer to [72].
Finally, let us discuss the Dark Matter phenomenology of the model. As explained above, the global symmetry is assumed to be conserved, and this implies that a stable state must exist. Assuming that the lightest state charged under is the neutral scalar , it constitutes the DM candidate in the model. One then needs to determine whether the observed DM relic density can be achieved in the region of parameter space where the anomalies are solved, without conflict with other experimental constraints. This is shown in Figure 6, where contours of are shown in the plane. This figure has been obtained with fixed and TeV. For each parameter point, the value of is chosen to reproduce the observed DM relic density, which is calculated using MicrOmegas [68]. Large values of are obtained in this way. For this reason, the most relevant DM annihilation channels for the determination of the relic density with these parameter values are . Even though the experimental constraints, in particular those from direct LHC searches for extra quarks, reduce the allowed parameter space substantially, one finds valid regions with . This value would explain the anomalies at ; see, for instance, [48]. Interestingly, the model is testable in future direct DM detection experiments, such as XENON1T, as shown in Figure 6. In the region of parameter space selected for this figure, the dominant process leading to DMnucleon scattering is 1loop photon exchange, with leptons running in the loop. The loop suppression is compensated by the large coupling.
The new states in this model can be discovered at the LHC. For instance, the heavy VL charged lepton can be produced in DrellYan processes. Due to the required large values for the coupling, this exotic state is expected to decay mainly to a DM particle (invisible at the LHC) and a muon. Since conservation requires the particles to be produced in pairs, the expected signature is the observation of two energetic muons and large missing energy. Similar events replacing the muons by jets (mainly jets) are expected for the VL quarks. We conclude the discussion of this model by referring for more details to the original work in [13].
7. Summary and Discussion
In this minireview we have discussed New Physics models that address the anomalies and link them to the dark matter of the Universe. The interplay between these two areas of particle physics may offer novel model building directions as well as additional phenomenological tests for the proposed scenarios. We have shown that most of the proposed models can be classified into two categories: (i) models in which the anomalies and the DM production mechanism share a common mediator (such as a heavy boson), and (ii) models that induce the NP contributions to explain the anomalies via loops including the DM particle. These generic ideas have been illustrated with two particular realizations (the models introduced in [2, 13]), which clearly show that the combination of flavor physics and dark matter leads to new scenarios with a rich phenomenology.
The introduction of a dark sector in a model for the anomalies can also have phenomenological consequences besides the existence of a DM candidate. For instance, both problems, the dark matter of the Universe and the anomalies, might be connected to another longstanding question in particle physics: the muon anomalous magnetic moment [4, 73–75]. Furthermore, it is interesting to note that the introduction of dark matter in some models can help alleviate some of the most stringent constraints. Indeed, the LHC bounds on some mediators become weaker if they have invisible decay channels [76]. We believe that this is a promising line of research to be pursued in order to fully assess the validity of some scenarios that are currently under experimental tension.
We are living an exciting moment in flavor physics, with several interesting anomalies in Bmeson decays. Whether real or not, only time can tell. New LHCb analyses based on larger datasets are expected to appear in the near future, possibly shedding new light on these anomalies. In the longer term, fundamental contributions from the Belle II experiment will also be crucial to settle the issue [77]. In the meantime, an intense model building effort is opening new avenues with rich phenomenological scenarios. The possible connection to one of the central problems in current physics, the nature of the dark matter of the Universe, would definitely be a fascinating outcome of this endeavour.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The author is grateful to Junichiro Kawamura, Shohei Okawa, and Yuji Omura for clarifications regarding [13] and to Javier Virto for comments about the manuscript. They are also very grateful to their collaborators in the subjects discussed in this review and acknowledge financial support from the Grants FPA201785216P and SEV20140398 (MINECO) and PROMETEOII/2014/084 (Generalitat Valenciana).
References
 Planck Collaboration, P. A. R. Ade et al., “Planck 2015 results. XIII. Cosmological parameters,” Astronomy & Astrophysics, vol. 594, 2016. View at: Publisher Site  Google Scholar
 D. A. Sierra, F. Staub, and A. Vicente, “Shedding light on the b→s anomalies with a dark sector,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 92, no. 1, Article ID 015001, 2015. View at: Publisher Site  Google Scholar
 G. Bélanger, C. Delaunay, and S. Westhoff, “A Dark matter relic from muon anomalies,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 92, no. 5, Article ID 055021, 2015. View at: Publisher Site  Google Scholar
 B. Allanach, F. S. Queiroz, A. Strumia, and S. Sun, “Z' models for the LHCb and g2 muon anomalies,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 93, no. 5, Article ID 055045, 2016, Erratum: Phys. Rev.D95,no.11,119902(2017). View at: Publisher Site  Google Scholar
 M. Bauer and M. Neubert, “Flavor anomalies, the 750 GeV diphoton excess, and a dark matter candidate,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 93, no. 11, 2016. View at: Publisher Site  Google Scholar
 A. Celis, W. Feng, and M. Vollmann, “Dirac dark matter and b→sl+l with U(1) gauge symmetry,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 95, no. 3, Article ID 035018, 2017. View at: Publisher Site  Google Scholar
 W. Altmannshofer, S. Gori, S. Profumo, and F. S. Queiroz, “Explaining dark matter and B decay anomalies with an LμLτ model,” Journal of High Energy Physics, vol. 12, no. 106, 2016. View at: Publisher Site  Google Scholar
 P. Ko, T. Nomura, and H. Okada, “A flavor dependent gauge symmetry, predictive radiative seesaw and LHCb anomalies,” Physics Letters. B. Particle Physics, Nuclear Physics and Cosmology, vol. 772, pp. 547–552, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 P. Ko, T. Nomura, and H. Okada, “Explaining B→K(*)l+l anomaly by radiatively induced coupling in U(1)μτ gauge symmetry,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 95, no. 11, Article ID 111701, 2017. View at: Publisher Site  Google Scholar
 J. M. Cline, J. M. Cornell, D. London, and R. Watanabe, “Hidden sector explanation of,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 95, no. 9, Article ID 095015, 2017. View at: Publisher Site  Google Scholar
 F. Sala and D. M. Straub, “A new light particle in B decays?” Physics Letters B, vol. 774, pp. 205–209, 2017. View at: Publisher Site  Google Scholar
 J. Ellis, M. Fairbairn, and P. Tunney, “Anomalyfree models for flavour anomalies,” The European Physical Journal C, vol. 78, no. 3, 2018. View at: Publisher Site  Google Scholar
 J. Kawamura, S. Okawa, and Y. Omura, “Interplay between the b→sll anomalies and dark matter physics,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 96, no. 7, Article ID 075041, 2017. View at: Publisher Site  Google Scholar
 S. Baek, “Dark matter contribution to b→sμ+μ anomaly in local U(1)LμLτ model,” Physics Letters B, vol. 781, pp. 376–382, 2018. View at: Publisher Site  Google Scholar
 J. M. Cline, “B decay anomalies and dark matter from vectorlike confinement,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 97, no. 1, Article ID 015013, 2018. View at: Publisher Site  Google Scholar
 J. M. Cline and J. M. Cornell, “R(K(⁎)) from dark matter exchange,” Physics Letters B, vol. 782, pp. 232–237, 2018. View at: Publisher Site  Google Scholar
 L. Dhargyal, “A simple model to explain the observed muon sector anomalies, small neutrino masses, baryongenesis and darkmatter,” https://arxiv.org/abs/1711.09772. View at: Google Scholar
 C.W. Chiang and H. Okada, “A simple model for explaining muonrelated anomalies and dark matter,” https://arxiv.org/abs/1711.07365. View at: Publisher Site  Google Scholar
 A. Falkowski, S. F. King, E. Perdomo, and M. Pierre, “Flavourful Z′ portal for vectorlike neutrino Dark Matter and RK(∗),” https://arxiv.org/abs/1803.04430. View at: Google Scholar
 G. Arcadi, T. Hugle, and F. S. Queiroz, “The Dark Lμ−Lτ Rises via Kinetic Mixing,” https://arxiv.org/abs/1803.05723. View at: Google Scholar
 S. Jäger and J. M. Camalich, “On B → Vℓℓ at small dilepton invariant mass, power corrections, and new physics,” Journal of High Energy Physics, vol. 2013, no. 5, 2013. View at: Publisher Site  Google Scholar
 S. DescotesGenon, L. Hofer, J. Matias, and J. Virto, “On the impact of power corrections in the prediction of B → K *μ+μ− observables,” Journal of High Energy Physics, vol. 2014, no. 12, 2014. View at: Publisher Site  Google Scholar
 A. Khodjamirian, T. Mannel, A. Pivovarov, and Y. Wang, “Charmloop effect in B → K (*) ℓ + ℓ − and B → K * γ,” Journal of High Energy Physics, vol. 2010, no. 9, 2010. View at: Publisher Site  Google Scholar
 M. Ciuchini, M. Fedele, E. Franco et al., “B → K∗ ℓ + ℓ − decays at large recoil in the Standard Model: a theoretical reappraisal,” Journal of High Energy Physics, vol. 2016, no. 6, 2016. View at: Publisher Site  Google Scholar
 B. Capdevila, S. DescotesGenon, L. Hofer, and J. Matias, “Hadronic uncertainties in B→K*μ+μ: a stateoftheart analysis,” Journal of High Energy Physics, vol. 2017, no. 4, 2017. View at: Publisher Site  Google Scholar
 V. G. Chobanova, T. Hurth, F. Mahmoudi, D. M. Santos, and S. Neshatpour, “Large hadronic power corrections or new physics in the rare decay B → K ∗ μ + μ −?” Journal of High Energy Physics, vol. 2017, no. 7, 2017. View at: Publisher Site  Google Scholar
 C. Bobeth, M. Chrzaszcz, D. van Dyk, and J. Virto, “Longdistance effects in B→K∗ℓℓ from Analyticity,” https://arxiv.org/abs/1707.07305. View at: Google Scholar
 T. Blake, U. Egede, P. Owen, G. Pomery, and K. A. Petridis, “An empirical model of the longdistance contributions to B0→K*0μ+μ transitions,” The European Physical Journal C, vol. 78, no. 453, 2018. View at: Publisher Site  Google Scholar
 S. DescotesGenon, T. Hurth, J. Matias, and J. Virto, “Optimizing the basis of B → K * ℓ + ℓ − observables in the full kinematic range,” Journal of High Energy Physics, vol. 2013, no. 5, 2013. View at: Publisher Site  Google Scholar
 LHCb Collaboration, R. Aaij et al., “Measurement of FormFactorIndependent Observables in the Decay B0→K*0μ+μ,” Physical Review Letters, vol. 118, no. 12, 2017. View at: Publisher Site  Google Scholar
 LHCb Collaboration, R. Aaij et al., “Differential branching fraction and angular analysis of the decay Bs0→ϕμ+μ,” JHEP, vol. 7, no. 84, 2013. View at: Google Scholar
 LHCb Collaboration, R. Aaij et al., “Angular analysis of the B0→K*0μ+μ decay using 3 fb 1 of integrated luminosity,” JHEP, vol. 2, no. 104, 2016. View at: Google Scholar
 LHCb Collaboration, R. Aaij et al., “Angular analysis and differential branching fraction of the decay Bs0→ϕμ+μ,” JHEP, vol. 09, no. 179, 2015. View at: Google Scholar
 Belle Collaboration, A. Abdesselam et al., “Angular analysis of B0→K*(892)0l+l,” in in Proceedings, LHCSki  A First Discussion of 13 TeV Results: Obergurgl, Austria, 2016, arXiv:1604.04042. View at: Google Scholar
 Belle Collaboration, S. Wehle et al., “LeptonFlavorDependent Angular Analysis of B→K*l+l,” Phys. Rev. Lett, vol. 118, no. 11, Article ID 111801, 2017. View at: Google Scholar
 ATLAS Collaboration, “Angular analysis of Bd0→K*μ+μ decays in pp collisions at s=8 TeV with the ATLAS detector,” in Proceedings of the collision8 TeV with the ATLAS detector, pp. 2017–023, 2017. View at: Google Scholar
 CMS Collaboration, “Measurement of the P1 and P5' angular parameters of the decay B0?K*0µ+µ in protonproton collisions at s=8TeV,” CMSPASBPH15008, 2017. View at: Google Scholar
 G. Hiller and F. Krüger, “More modelindependent analysis of b→s processes,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 69, no. 7, 2004. View at: Publisher Site  Google Scholar
 LHCb Collaboration, R. Aaij et al., “Test of lepton universality using B+→K+l+l decays,” Physical Review Letters, vol. 113, Article ID 151601, 2014. View at: Publisher Site  Google Scholar
 LHCb Collaboration, R. Aaij et al., “Test of lepton universality with B 0 → K *0 ℓ + ℓ − decays,” Journal of High Energy Physics, vol. 2017, no. 55, 2017. View at: Publisher Site  Google Scholar
 S. DescotesGenon, L. Hofer, J. Matias, and J. Virto, “Global analysis of b→sll anomalies,” Journal of High Energy Physics, vol. 1606, article 092, no. 6, 2016. View at: Publisher Site  Google Scholar
 M. Bordone, G. Isidori, and A. Pattori, “On the Standard Model predictions for RK and RK*,” The European Physical Journal C, vol. 76, no. 8, 2016. View at: Publisher Site  Google Scholar
 B. Capdevila, S. DescotesGenon, J. Matias, and J. Virto, “Assessing leptonflavour nonuniversality from B →K ∗ ℓℓ angular analyses,” Journal of High Energy Physics, vol. 2016, no. 10, 2016. View at: Publisher Site  Google Scholar
 Heavy flavor averaging group, http://www.slac.stanford.edu/xorg/hfag/semi/fpcp17/RDRDs.html.
 LHCb Collaboration, R. Aaij et al., “Measurement of the Ratio of Branching Fractions B(B+c→J/ψτ+ντ)/B(B+c→J/ψμ+νμ),” Physical Review Letters, vol. 120, Article ID 121801, 2018. View at: Publisher Site  Google Scholar
 D. Buttazzo, A. Greljo, G. Isidori, and D. Marzocca, “Bphysics anomalies: a guide to combined explanations,” Journal of High Energy Physics, vol. 2017, no. 11, 2017. View at: Publisher Site  Google Scholar
 B. Capdevila, A. Crivellin, S. DescotesGenon, J. Matias, and J. Virto, “Patterns of New Physics in b → sℓ+ℓ− transitions in the light of recent data,” Journal of High Energy Physics, vol. 2018, no. 1, 2018. View at: Publisher Site  Google Scholar
 W. Altmannshofer, P. Stangl, and D. M. Straub, “Interpreting hints for lepton flavor universality violation,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 96, no. 5, 2017. View at: Publisher Site  Google Scholar
 G. D’Amico, M. Nardecchia, P. Panci et al., “Flavour anomalies after the RK* measurement,” Journal of High Energy Physics, vol. 2017, no. 9, 2017. View at: Publisher Site  Google Scholar
 G. Hiller and I. Nišandžić, “RK and RK* beyond the standard model,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 96, no. 3, Article ID 035003, 2017. View at: Publisher Site  Google Scholar
 L. Geng, B. Grinstein, S. Jäger, J. Martin Camalich, X. Ren, and R. Shi, “Towards the discovery of new physics with leptonuniversality ratios of b→sll decays,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 96, no. 9, 2017. View at: Publisher Site  Google Scholar
 M. Ciuchini, A. M. Coutinho, M. Fedele et al., “On flavourful Easter eggs for New Physics hunger and lepton flavour universality violation,” The European Physical Journal C, vol. 77, no. 10, 2017. View at: Publisher Site  Google Scholar
 A. K. Alok, B. Bhattacharya, A. Datta, D. Kumar, J. Kumar, and D. London, “New Physics in b→sμ+μ after the Measurement of RK*,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 96, no. 9, 2017. View at: Publisher Site  Google Scholar
 T. Hurth, F. Mahmoudi, D. Martínez Santos, and S. Neshatpour, “Lepton nonuniversality in exclusive b→sll decays,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 96, no. 9, 2017. View at: Publisher Site  Google Scholar
 A. Celis, J. FuentesMartín, A. Vicente, and J. Virto, “Gaugeinvariant implications of the LHCb measurements on leptonflavor nonuniversality,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 96, no. 3, 2017. View at: Publisher Site  Google Scholar
 W. Altmannshofer, S. Gori, M. Pospelov, and I. Yavin, “Quark flavor transitions in LμLτ models,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 89, no. 9, 2014. View at: Publisher Site  Google Scholar
 AMS Collaboration, M. Aguilar et al., “Antiproton Flux, AntiprotontoProton Flux Ratio, and Properties of Elementary Particle Fluxes in Primary Cosmic Rays Measured with the Alpha Magnetic Spectrometer on the International Space Station,” Physical Review Letters, vol. 117, no. 9, Article ID 091103, 2016. View at: Google Scholar
 B. Gripaios, M. Nardecchia, and S. Renner, “Linear flavour violation and anomalies in B physics,” Journal of High Energy Physics, vol. 2016, no. 6, 2016. View at: Publisher Site  Google Scholar
 P. Arnan, L. Hofer, F. Mescia, and A. Crivellin, “Loop effects of heavy new scalars and fermions in b→sμ+μ,” Journal of High Energy Physics, vol. 2017, no. 4, 2017. View at: Publisher Site  Google Scholar
 E. Ma, “Verifiable radiative seesaw mechanism of neutrino mass and dark matter,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 73, Article ID 077301, 2006. View at: Publisher Site  Google Scholar
 P. RochaMoran and A. Vicente, “Lepton flavor violation in a model for the anomalies,” In press. View at: Google Scholar
 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,” Journal of High Energy Physics, vol. 2013, no. 2, 2013. View at: Publisher Site  Google Scholar
 W. Altmannshofer and D. M. Straub, “New physics in b→s transitions after LHC run 1,” The European Physical Journal C, vol. 75, no. 8, 2015. View at: Publisher Site  Google Scholar
 L. M. Krauss and F. Wilczek, “Discrete gauge symmetry in continuum theories,” Physical Review Letters, vol. 62, no. 11, pp. 1221–1223, 1989. View at: Publisher Site  Google Scholar
 B. Petersen, M. Ratz, and R. Schieren, “Patterns of remnant discrete symmetries,” Journal of High Energy Physics, vol. 2009, no. 08, article 111, 2009. View at: Publisher Site  Google Scholar
 D. A. Sierra, M. Dhen, C. S. Fong, and A. Vicente, “Dynamical flavor origin of ZN symmetries,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 91, no. 9, 2015. View at: Publisher Site  Google Scholar
 W. Porod, F. Staub, and A. Vicente, “A flavor kit for BSM models,” The European Physical Journal C, vol. 74, no. 8, article 2992, 2014. View at: Publisher Site  Google Scholar
 G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, “micrOMEGAs_3: a program for calculating dark matter observables,” Computer Physics Communications, vol. 185, no. 3, pp. 960–985, 2014. View at: Publisher Site  Google Scholar
 A. Greljo, G. Isidori, and D. Marzocca, “On the breaking of lepton flavor universality in B decays,” Journal of High Energy Physics, vol. 2015, no. 7, 2015. View at: Publisher Site  Google Scholar
 A. Alves, S. Profumo, and F. S. Queiroz, “The dark Z' portal: direct, indirect and collider searches,” Journal of High Energy Physics, vol. 2014, article 63, 2014. View at: Publisher Site  Google Scholar
 A. Alves, A. Berlin, S. Profumo, and F. S. Queiroz, “Dark Matter Complementarity and the Z' Portal,” Physical Review D, vol. 92, no. 8, Article ID 083004, 2015. View at: Google Scholar
 L. Di Luzio, M. Kirk, and A. Lenz, “Updated Bsmixing constraints on new physics models for b→sℓ+ℓ− anomalies,” Physical Review D, vol. 97, Article ID 095035, 2018. View at: Publisher Site  Google Scholar
 G. Bélanger and C. Delaunay, “A Dark Sector for gμ2, RK and a Diphoton Resonance,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 94, no. 7, Article ID 075019, 2016. View at: Publisher Site  Google Scholar
 S. Di Chiara, A. Fowlie, S. Fraser et al., “Minimal flavorchanging Z′ models and muon g2 after the RK* measurement,” Nuclear Physics B, vol. 923, pp. 245–257, 2017. View at: Publisher Site  Google Scholar
 K. Kowalska and E. M. Sessolo, “Expectations for the muon g2 in simplified models with dark matter,” Journal of High Energy Physics, vol. 2017, no. 9, article 112, 2017. View at: Publisher Site  Google Scholar
 D. A. Faroughy, A. Greljo, and J. F. Kamenik, “Confronting lepton flavor universality violation in B decays with high pT tau lepton searches at LHC,” Physics Letters B, vol. 764, pp. 126–134, 2017. View at: Publisher Site  Google Scholar
 J. Albrecht, F. Bernlochner, M. Kenzie, S. Reichert, D. Straub, and A. Tully, “Future prospects for exploring present day anomalies in flavour physics measurements with Belle II and LHCb,” https://arxiv.org/abs/1709.10308. View at: Google Scholar
Copyright
Copyright © 2018 Avelino Vicente. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.