$R_K$ and $R_{K^*}$ in an Aligned 2HDM with Right-Handed Neutrinos

We consider the possibility to explain the recent $R_K$ and $R_{K ^*}$ anomalies in a 2-Higgs Doublet Model, known as Aligned, combined with a low scale seesaw mechanism generating light neutrino masses and mixings. In this class of models, a large Yukawa coupling allows for significant non-universal leptonic contributions, through box diagrams mediated by charged Higgs bosons and right-handed neutrinos, to the $b \to s \ell^+ \ell^-$ transition that can then account for both $R_K$ and $R_{K^*}$ anomalies.

1 [3], by a ≈ 2.5σ deviation. Hence, they are considered as important hints of new physics that violates lepton universality.
Any signal of possible lepton non-universality would be striking evidence for physics Beyond the SM (BSM). Therefore, the R K and R K * anomalies attracted the attention of theoretical particle physicists and several new physics scenarios have been proposed to accommodate these results, for a review, see [4] and related references.
A non-trivial flavour structure in the lepton sector is already required, beyond the SM, in order to explain the observation of neutrino oscillations. Therefore, it is plausible and very attractive if the mechanism behind the B anomalies can be related to the same physics responsible for the non-zero neutrino masses and oscillations. In view of this, one then ought to explore extensions of the SM able to address both phenomena.
One of the possible scenarios for generating the observed lepton non-universality is to allow for large Yukawa couplings of the right-handed neutrinos with Higgs fields and charged leptons, as in a low scale seesaw mechanism, with the inverse seesaw being one of the most notorious examples, for generating light neutrino masses. In this case, a non-trivial neutrino Yukawa matrix may lead to different results for BR(B → Ke + e − ) and BR(B → Kµ + µ − ). This framework has been recently considered in the Supersymmetric (SUSY) B−L extension of the SM, where it was shown that that the box diagram mediated by a right-handed sneutrino, higgsino-like chargino and light stop can account simultaneously for both R K and R K * [5]. In non-SUSY models, a similar box diagram can be obtained through a charged Higgs, instead of a chargino, and a right-handed neutrino, instead of a right-handed sneutrino. Therefore, a rather minimal model that can account for these discrepancies is an extension of the SM with two Higgs doublets (so that we can have a physical charged Higgs boson) and right-handed neutrinos along with a low scale seesaw mechanism (to guarantee large neutrino Yukawa couplings). Note that there have been several attempts at explaining the above results through the penguin diagram as well, with a non-universal Z , and also through tree level mediation of flavour violating Z or leptoquark that induce a non-universal b → s + − transition. However, these types of flavour violating interactions are subject to severe experimental limits.
The 2-Higgs Doublet Model (2HDM), which is motivated by SUSY and Grand Unified Theories (GUTs), is the simplest model that includes charged Higgs bosons. According to the types of couplings of the two Higgs doublets to the SM fermions doublets and singlets, we may have different type of 2HDMs. For example, if only one Higgs doublet couples to the SM fermions one obtains the type I 2HDM. While, in the case of one Higgs doublet coupling to the up quarks and the second Higgs doublet coupling to the down quarks and charged leptons, one obtains the type II 2HDM. Also, we may have a type III or IV 2HDM if both Higgs doublets couple to both up and down quarks as well as charged leptons. However, severe constraints are imposed on 2HDMs due to their large contributions to Flavour Changing Neutral Currents (FCNCs) that contradict the current experimental limits. Therefore, assumptions on the Yukawa couplings generated by different Higgs doublets are usually imposed. One of these assumptions is the alignment between the Yukawa couplings generated by Φ 1 and Φ 2 , the two Higgs doublet fields. This class of models is called the Aligned 2HDM (A2HDM). It is worth mentioning that, as discussed below, other 2HDMs cannot account for the LHCb results of lepton non-universality.
In this paper, we emphasise that in the A2HDM, extended by (heavy) right-handed neutrinos in order to generate the (light) neutrino masses through a seesaw mechanism, interesting results can be obtained for several flavour observables, such as µ → eγ, B s → µµ and, indeed, R K ( * ) . In particular, one can account simultaneously for both aforementioned results on R K and R K * , through a box diagram mediated by a right-handed neutrino, top quark and charged Higgs boson.
The paper is organised as follows. In section 2 we introduce the main features of our A2HDM focusing on the neutrino sector and its interplay with the Higgs structures. (Here, we also briefly review some particular realisations of low scale seesaw models for generating light neutrino masses.) In section 3 we discuss the most relevant constraints from flavour physics that affect our A2HDM parameter space. The calculation of the A2HDM contributions to b → s + − transitions mediated by charged Higgs bosons and right-handed neutrinos is given in section 4. Our numerical results are presented in section 5. Finally, our conclusions and remarks are given in section 6.

The A2HDM
The 2HDM is characterised by two Higgs doublets with hypercharge Y = 1/2 which, in the Higgs basis, can be parameterised as where v is the Electro-Weak (EW) Vacuum Expectation Value (VEV) and G + , G 0 denote the Goldstone bosons. The two doublets describe five physical scalar degrees of freedom which are given by the two components of the charged Higgs H ± and three neutral states ϕ 0 i = {h, H, A}, the latter obtained from the rotation of the φ 0 i fields into the mass eigenstate basis. The scalar squared mass matrix M 2 S is determined by the structure of the 2HDM scalar potential, see for instance [6][7][8], and diagonalised by the orthogonal matrix R, where In general, the three mass eigenstates ϕ 0 i do not have definite CP transformation properties but in the CP-conserving scenario φ 0 3 does not mix with the other two neutral states and the scalar spectrum consists of a CP-odd field A = φ 0 3 and two CP-even fields h and H that are defined from the interaction eigenstates through the two-dimensional orthogonal matrix The most general Yukawa Lagrangian of the 2HDM can be written in the Higgs basis as where the quark Q L , u R , d R and lepton L R , R , ν R fields are defined in the weak interaction basis and we also included the couplings of the left-handed lepton doublets with the righthanded neutrinos. The Φ 1,2 fields are the two Higgs doublets in the Higgs basis and, as customary,Φ i = iσ 2 Φ * i . The Yukawa couplings Y 1j and Y 2j , with j = u, d, , are 3 × 3 complex matrices while Y 1ν and Y 2ν are 3×n R matrices, with n R being the number of righthanded neutrinos. In general, the Yukawas Y 1 and Y 2 cannot be simultaneously diagonalised in flavour space, so that, whilst the quark and the charged-lepton Y 1 can be recast into a diagonal form in the fermion mass eigenstate basis, namely, Y 1 = √ 2/vM , with M being the fermion mass matrix, Y 2 would remain non-diagonal and thus give rise to potentially dangerous tree-level FCNCs. This problem is usually solved by enforcing that only one of the two Higgs doublets couple to a given right-handed field. This requirement is satisfied by implementing a discrete Z 2 symmetry acting on the Higgs and fermion fields. There are four non-equivalent choices: type I, II, III and IV (as previously intimated). Another general way to avoid tree-level FCNCs in the Higgs sector is to require the alignment, in flavour space, of the two Yukawa matrices that couple to the same right-handed fermion [9], namely, where the proportionality constants ζ f are arbitrary family universal complex parameters. This scenario is dubbed A2HDM. The allowed sources of FCNCs at quantum level are highly constrained and the resulting structures are functions of the mass matrices and Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, so that this model provides an explicit example of the popular Minimal Flavour Violation (MFV) scenario [10]. Even though the alignment of the Yukawa matrices is strictly required, from observations, only in the quark and charged lepton sectors, we assume that the same mechanism which guarantees the aligned structure in the SM flavour space also holds in the neutrino sector and leads to In all sectors, the alignment is fixed to be exact at some specified scale µ 0 and subsequently will misalign due to radiative corrections, as discussed in [11,12]. However, the flavour structure of the model constrain the nature of the new sourced of FCNCs induced by Renormalisation Group effects. Quantitatively, in the quark sector, these FCNC contributions are suppressed by mass hierarchies m q m 2 q /v 3 and provide negligible effects [11,12]. We will not consider the impact of the misalignment in this work.
Interestingly, ζ f can provide new sources of CP violation but in this work we will consider only real values. Notice also that the usual 2HDMs in which tree-level FCNCs are removed by exploiting the discussed Z 2 discrete symmetry, namely the type I, II, III and IV, can be recovered for particular values of the proportionality constants ζ f as shown in Tab. 1.
Aligned Type I Type II Type III Type IV Table 1. Relation between the ζ f couplings of the A2HDM and the ones of the Z 2 symmetric scenarios.
The Yukawa Lagrangian in Eq. (2.4) generates a Dirac mass matrix for the standard neutrinos and can also be supplemented by a Majorana mass term M R for the right-handed where C is the charge-conjugation operator. In particular, by exploiting a bi-unitary transformation in the charged-lepton sector and a unitary transformation on the right-handed neutrinos L L = U L L L , R = U R R and ν R = U ν R ν R is always possible to diagonalise (with real eigenvalues) the charged-lepton and Majorana mass matrices at the same time, In this basis the neutrino mass matrix can be written as such that M ν = U T MU provides the masses of the three light active neutrinos ν l and of the remaining n R heavy sterile neutrinos ν h . The Yukawa interactions of the physical scalars with the mass eigenstate fermions are then described by where the couplings of the neutral Higgs states to the fermions are given by Due to the alignment of the Yukawa matrices all the couplings of the scalar fields to fermions are proportional to the corresponding mass matrices. Finally, the weak neutral and charged interactions of the neutrinos are In this paper, rather than presenting a complete model in the neutrino sector by specifying the structure and the hierarchies of the neutrino mass matrices, we work in a simplified framework that captures the interesting phenomenology whilst preserving a significant degree of model independence. In particular, we consider a single extra heavy neutrino despite the usual requirement of additional sterile states to fully accommodate the observed pattern of the light neutrino masses and mixing angles. Indeed, low scale right-handed neutrinos with sizeable mixings with the SM left-handed neutrino states, such that they may provide visible effects in physical observables at the EW scale, usually affect the light neutrino masses with inadmissible large contributions. This issue is nicely solved in extended seesaw models [13][14][15][16][17][18][19][20], as for example the linear or inverse seesaw mechanism, in which extra sterile neutrino states are introduced to allow for large mixings while correctly reproducing the observed smallness of the light neutrino masses. For the sake of definiteness, in the following section we briefly present some specific setup that could be employed to realise the phenomenological scenario described above.

Some explicit examples of low scale seesaw mechanism
As stated above, in order to achieve a sizeable mixing with the heavy sterile neutrino while avoiding, at the same time, large contributions to the light neutrino masses, it is necessary to require the neutrino Majorana mass matrix and the Dirac neutrino Yukawa coupling to realise a particular structure. In order to match the nomenclature usually employed in the literature, we split the set of right-handed neutrinos defined above into two classes, one of n S SM-singlet fermionic fields S i and another one of fermionic states that we continue to call right-handed neutrinos. The differences between the two will be clear in a moment.
In the basis N L = (ν L , ν c R , S) T , the mass matrix can be parameterised as which has the same structure of the one given in Eq. (2.9) provided that The m D and m S are, respectively, 3×n R and 3×n S mass matrices mediating the interactions between the charged leptons and the right-handed and sterile neutrinos while m N , m R and µ S are n R × n R , n R × n S and n S × n S mass matrices, respectively. As both right-handed ν R and sterile S neutrinos can be assigned lepton number L = 1, the mass terms m N , m S and µ S violate lepton number by two units. Two commonly studied mass patterns are the inverse and linear seesaws, which are characterised by m S = m N = 0 and µ S = m N = 0, respectively. In these cases, the vanishing of the Majorana mass µ S or m S would restore lepton number conservation and, as such, would increase the symmetry of the model. This feature makes the two masses naturally small accordingly to 't Hooft naturalness principle.
Following the standard seesaw calculation and by assuming the hierarchy µ S (m S ) m D , m R for the inverse (linear) seesaw scenario, the 3 × 3 light neutrino mass matrix is Differently from the standard type I seesaw case, in which m ν ∼ m 2 D /m R with m D m R , the lightness of the active neutrino masses is ensured in these low scale seesaw scenarios by the smallness of the µ S (m S ) parameters. This feature prevents m D to be extremely suppressed with respect to the Majorana mass and, as such, may allow for non-negligible couplings between the heavy neutrinos and the SM gauge bosons which are set by the mixing U Lh .
In order to understand the dependence of the mixing U Lh of the left-handed SM neutrinos with the extra sterile states, it is instructive to study the mixing matrix U in the limit of negligible µ S (m S ) and small m D /m R . While the first requirement is necessary to reproduce the lightness of the active neutrino states, the latter is used here only to simplify the structure of U which reads as Notice that the PMNS matrix has been set to the unit one, consistently with the approximation m ν i 0. From Eq. (2.18), one can immediately realise that U Lh is set, as naively expected from dimensional arguments, by the ratio m D /m R .
Once some specific inputs are provided for m D and m R , the corresponding µ S (m S ) matrix that ensures the agreement with the light neutrino mass splittings and mixing angles can always be reconstructed from Eq. (2.16). For instance, in the inverse seesaw case, we find where m ν and U PMNS are chosen in agreement with the bounds from the low-energy neutrino data which we report below for the sake of completeness. In particular, one should enforce the following constraints from the latest results of the Vfit group [21] extracted from the Vfit 3.2 (2018) data.

1) Neutrino mass squared differences
The 3σ Confidence Level (CL) ranges on the mass squared differences Notice also that m D and m R cannot be chosen freely since the U Lh block of the mixing matrix that they define is constrained by the unitarity requirement that directly affects the analysis presented in this work. The corresponding bound is discussed in the next section together with all the other relevant constraints.

Unitarity bounds on the neutrino mixing matrix
The 3 × 3 block of the mixing matrix U corresponds to a non-unitaryŨ PMNS matrix. The bounds on the deviation from unitarity ofŨ PMNS have been obtained in [22,23] using an effective field theory approach in which the masses of the heavy neutrinos lie above the EW scale. This constraint can be recast as follows

Lepton flavour violating processes
We consider the lepton flavour violating decays α → β γ induced at one-loop order by the sterile neutrinos and check their compatibility with the experimental upper bounds at 90% CL [24], The Branching Ratios (BRs) of the aforementioned decay rates are given by where Γ α is the total decay width of the lepton α and the loop functions are where we have neglected the mass of the lepton in the final state.
The G H ± can offer large contributions, larger than G W ± , for sizeable values of the couplings ζ ν , ζ l , which are, as such, strongly constrained by lepton flavour violating processes. These can be tamed by controlling the size of the mixing matrix elements (U Lh ) αi which should be highly suppressed in order to avoid any large contribution to these sensitive processes.

Flavour constraints from meson processes
Here we briefly mention the relevant constraints from measurement of flavour observables in meson mixing and decays. These have been studied in the context of general 2HDMs and the majority of the bounds extracted from these can be straightforwardly applied in our case since the presence of the sterile neutrinos does not add any significant contribution at leading order. In particular, the 2HDM with the alignment in the flavour sector has been scrutinised in [11,12,25,26] to which we refer for delineating the allowed parameter space spanned by the ζ u,d, couplings and the charged Higgs mass M H ± .
• Neutral meson mixing The meson mixing observables ∆M s , ∆M d and | K | constrain large value of the ζ u parameter. In turn, they do not significantly affect ζ d whose dependence is suppressed by m 2 b /M 2 W in the first two observables while does not appear at all in the third one. Obviously, there is no dependence on ζ and ζ ν .
• Radiative B s → X s γ decay The B s → X s γ decay rate represents one of the best measured observables and it is employed to constrains several new physics scenarios. The contribution of the charged Higgs boson is encoded in the Wilson coefficients C 7 and C 8 . At leading order, the corresponding new physics corrections are sensitive to ζ u and ζ d and are given by (3.7) As for the new physics contributions to the meson mixing observables, the values of the ζ and ζ ν parameters are completely irrelevant in the determination of the b → sγ transition.
• Leptonic decay of the neutral mesons B 0 q → µ + µ − The BR of the B 0 q → µ + µ − meson decay is given by The leading new physics contribution appears in the Wilson coefficient C 10 and it is mediated by the charged Higgs in 2HDMs. Interestingly, the same coefficient is also corrected by the presence of the heavy sterile neutrinos and it is sensitive to, besides ζ u,d, , the coupling ζ ν of the charged Higgs to the sterile neutrino states. The impact of the C P and C S coefficients is, in contrast, suppressed by the m 2 factor unless the former can be enhanced with respect to C SM 10 as, for instance, in a Z 2 symmetric model with large tan β. In the flavour aligned 2HDM and in the parameter space in which we are interested in, namely ζ u 1 and ζ d ζ 0, the bound from the measurement of the B 0 s → µ + µ − transition is usually weaker than the one from the B s → X s γ decay. Nevertheless, due to the contribution from the heavy sterile neutrinos which is proportional to the new parameter ζ ν , we recompute the corresponding constraint by employing the flavio package [27].
• Leptonic decay of the charged mesons M ± → τ ± ν The M ± → τ ± ν decay occurs at tree level through charged current processes and the corresponding BR is where f M and τ M are the decay constant and the lifetime, respectively, and the contribution of the charged Higgs boson is encoded in (3.11) The most constraining decay mode is found to be B → τ ν which is, however, only relevant for light charge Higgs masses.

Contributions to b → s + − processes
The effective Hamiltonian for the b → s + − transitions is given by where the relevant operators for the analysis of the R K and R K * anomalies are The new physics effects in the corresponding Wilson coefficients can be recasted as where ∆ Z , ∆ γ and denote the contributions from the Z penguins, the photon penguins and the box diagrams, respectively. Moreover, C

(•)
i (H ± ) represents the charged Higgs contribution typical of the 2HDM, C Penguins: contributions from the charged Higgs boson 10 (H ± ) = 0. Box: contributions from the heavy neutrinos Box: contributions from the charged Higgs boson and the heavy neutrinos In the previous equations we used the following mass ratios We now discuss some interesting features concerning the structure of these coefficients. Firstly, one can see for the box diagram contribution that the heavy neutrinos give C 9 = −C 10 . One can also see that the SM box contribution with the light neutrinos, not shown above, is rescaled by with the heavy neutrinos and virtual W ± s. This is proportional to n R i=1 |(U Lh ) i | 2 and we do not expect, as confirmed by the numerical analysis, that this contribution can provide large effects onto flavour observables. Indeed, the coupling of the heavy neutrinos to the leptons mediated by the W ± boson is fixed by the gauge invariance and proportional to the SU (2) weak gauge coupling. In order to allow for more freedom one has to rely on an extra charged degree of freedom with the simplest possibility being the charged scalar of a 2HDM extension. These contributions are encoded into C ( ) i (N R , H ± ). The Z 2 symmetric scenarios of the 2HDM are among the simplest ones but barely produce significant effects in the C 9,10 Wilson coefficients. This can be understood from Tab. 1, because the corrections to C 9,10 would be proportional to ζ 2 u = ζ 2 ν = cot 2 β, independently from the specific realisation, and thus relevant only for tan β < 1, which is severely constrained by b → sγ. The A2HDM allows to disentangle ζ u from ζ ν , such that, while the former is still bound from b → sγ, the latter can be varied freely thus providing significant contributions to the Wilson coefficients in some region of the parameter space. We recall again that the alignment in the neutrino sector is not strictly required by the flavour physics but we, nevertheless, impose it by assuming that the same mechanism ensuring the proportionality between the Yukawa couplings is in place in both the quark and lepton sectors.  Among the different lepton flavour violating processes, µ → eγ is the most constraining one. In Fig. 2 we show the corresponding BR generated by a single heavy neutrino as a function of the squared mixing angle of the same heavy neutrino with the electron one. The other parameters have been fixed as M H ± = 700 GeV, m ν h i = 500 GeV, the ζ = 0 and |(U Lh ) µi | 2 = 0.4 × 10 −3 . The latter corresponds to the maximum allowed value from unitarity constraints, see Eq. (3.1). Since in a model with a single Higgs doublet we would obtain BR/U 4 ∼ (6 − 7) × 10 −4 in the heavy neutrino mass range of 500 − 1000 GeV, it is clear that the diagrams with the charged Higgs give the dominant and a very large contribution to the BR. Using the largest values for the mixing angles allowed by unitarity we obtain a suppression of the BR of ∼ 10 −7 which cannot accommodate the strong bound on µ → eγ from the MEG experiment. This suggest that, if a low scale seesaw is embedded into a 2HDM framework, a given heavy neutrino (or, in a scenario with a large hierarchy between heavy neutrino states, the lightest one) may have a non-negligible mixing only with SM neutrinos of a given flavour eigenstate otherwise the charged Higgs boson would induce unacceptably large effects on lepton flavour violating processes. The two realisations would be (U Lh ) µi = 0, (U Lh ) ei 0 or (U Lh ) µi 0, (U Lh ) ei = 0. As we will see below, the first possibility can be also used to explain the deviation of R K * and R K from the SM prediction. These two conditions can be achieved by suitably choosing the mass matrices m D and m R . Another possibility, allowing one to control the large effects in BR(µ → eγ), which, anyway, we will not explore in this work, is to assume a ζ = 0 and tune it against the ζ 2 ν term in Eq. (3.5) to reduce the G H ± form factor with respect to G W ± .

Results
With the analytic expressions of the Wilson coefficients, we may turn to finding realistic benchmark scenarios. We explore the parameter space spanned by the three parameters: ζ ν , m ν i and (U Lh ) 2 µ,i considering the impact of a single heavy neutrino, thus assuming, for the sake of simplicity, a hierarchy in the neutrino mass spectrum. The general case can be obtained straightforwardly and does not add much to the present discussion. In particular, one finds that if all the neutrino masses are almost degenerate, the combination 6 i=1 (U Lh ) 2 µ,i can be factored out from the Wilson coefficients and can be treated as an independent parameter leading to the same conclusions of the single neutrino case. This combination of squared mixing angles is bound from the non-unitarity test of the PMNS matrix to be less than ∼ 0.4 × 10 −3 . Notice that we considered the scenario (U Lh ) µi = 0, (U Lh ) ei 0 in line with the previous discussion on lepton flavour violating processes. The performed scan takes the ranges: −80 < ζ ν < 80, 10 −5 < (U Lh ) 2 µ,i < 10 −3 and 200 GeV < m ν i < 2000 GeV. The other parameters are chosen as follows: M H ± = 550 GeV, ζ u = 1 and ζ d = ζ 0 which ensure that the flavour constraints discussed above, namely the ones which are not significantly affected by the presence of the heavy neutrinos, are all satisfied [26]. Other choices are obviously acceptable but not considered here.
The main result of the analysis is presented in Fig. 3. The numerical values of the flavour observables R K * and R K have been obtained from the evaluation of the Wilson coefficients computed above and from the flavio package [27]. In the interesting region of the parameter space, the predictions for the two ratios in the A2HDM are the same and, for the sake of simplicity we will only discuss R K * . The latter is presented as a function of the parameter ζ ν Y ν ≡ √ 2 ζ ν (U Lh ) µi (m ν h i /v) that mostly controls the C 9,10 (a) (b) Figure 3. R K * in the central bin and BR(B s → µ + µ − ) as a function of the combination of parameters that mostly controls the observable. The green bands are the 1σ and 2σ bands for the R K * measurement. Blue (red) points correspond to the scenario in which the heavy neutrino has a non-negligible coupling only to the muon (electron) flavour eigenstates.  Wilson coefficients. In particular, the blue points correspond to the configuration in which the heavy neutrino couples to the SM muon sector and has negligible mixing with the first family. As anticipated above, this setup allows to reproduce the measured reduction in the R K * ratio which is represented in the plot with a dashed horizontal line, together with 1σ and 2σ (green) bands. The red points, instead, are representative of the scenario in which the heavy neutrino has a mixing with the electrons. In this case the predicted R K * is above one and contradicts the LHCb observations. The extent of the reduction is mainly controlled by the parameter ζ ν Y ν . Interestingly, the same parameter defines the strength of the coupling among the charged Higgs, heavy neutrino and lepton and, as such, is subject to the perturbativity bound of [28]. The upper limit for the coupling is usually extracted from naive scaling arguments but also depends on the loop functions of the involved processes.
In the plot, to facilitate the reading, we have shown two commonly adopted upper bounds that can be used to estimate the region of tree-level perturbativity. With a ζ ν Y ν coupling below the most stringent perturbativity bound it is possible to explain the present R K * measurement and the BR(B s → µ + µ − ) one within the 2σ interval. Fig. 4(a) displays the 1σ bounds for both the R * K low and central bins. One can see that parameter configurations can satisfy either the low or central bins separately, but not simultaneously both. In Fig. 4(b) we consider the effects on the BR(B s → µ + µ − ). One can see that for several points both predictions are simultaneously compatible with the experimental measurements.

Conclusions
The LHCb experiment at CERN has recently reported the existence of some anomalies in their data, with respect to the predictions of the SM. Specifically, the measured values of the observables R K = BR(B + → K + µ + µ − )/BR(B + → K + e + e − ) and R K * = BR(B 0 → K * 0 µ + µ − )/BR(B 0 → K * 0 e + e − ) revealed a ≈ 2.5σ deviation when compared to the SM rates, which are essentially 1. In addition, the discrepancies occur in two di-lepton invariance mass bins. Therefore, these results must be taken as a serious hint of possible BSM physics.
Herein, we have considered the possibility of explaining such R K and R K * anomalies in an A2HDM, wherein an alignment is present between the Yukawa couplings generated by the two Higgs doublet fields, combined with a low scale seesaw mechanism generating light neutrino masses and mixings in compliance with current experimental measurements. Such a scenario allows for significant non-universal leptonic contributions, through box diagrams mediated by H ± and ν R states, which in turn alter the yield of the partonic decay b → s + − entering the definition of both the R K and R K * observables. In order to render our explanation phenomenologically viable, we have made sure to comply with both theoretical (chiefly, the unitarity bounds stemming from the neutrino mixing matrix) and experimental (the strongest being those due to lepton flavour violating processes and mesonic decay channels) constraints. In fact, the masses required for the charged Higgs and right-handed neutrino states entering the above transition are also well beyond their current direct limits. Furthermore, in order to make clear that our explanation of the R K and R K * anomalies is not particularly ad hoc, we have left the actual low scale dynamics onsetting the seesaw mechanism undetermined, by illustrating that this could be realised through different scenarios, e.g., the so-called inverse and linear seesaw cases. Therefore, our setup captures a variety of light neutrino masses and mixings that can be tuned to further experimental observation in the neutrino sector while leaving predictions in the B one unchanged. Finally, we have correlated our predictions for R K and R K * in both dilepton invariant mass bins to those for the highly constraining observable BR(B s → µ + µ − ) showing that a simultaneous solutions to both sets of measurements can be found in the envisioned A2HDM plus low scale seesaw scenario.