Update of the flavour-physics constraints in the NMSSM

We consider the impact of several flavour-changing observables in the B- and the Kaon sectors on the parameter space of the NMSSM, in a minimal flavour violating version of this model. Our purpose consists in updating our previous results in [4] and designing an up-to-date flavour test for the public package NMSSMTools. We provide details concerning our implementation of the constraints in a series of brief reviews of the current status of the considered channels. Finally, we present a few consequences of these flavour constraints for the NMSSM, turning to two specific scenarios: one is characteristic of the MSSM-limit and illustrates the workings of charged-Higgs and genuinely supersymmetric contributions to flavour-changing processes; the second focus is a region where a light CP-odd Higgs is present. Strong limits are found whenever an enhancement factor – large tanβ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tan \beta $$\end{document}, light H±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H^{\pm }$$\end{document}, resonant pseudoscalar – comes into play.


Introduction
Flavour-changing rare decays and oscillation parameters are known as uncircumventable tests of the Standard Model (SM) and its new-physics extensions. In the quark sector of the SM, flavour-violation is induced by the non-alignment of the Yukawa matrices, resulting in a Cabbibbo-Kobayashi-Maskawa (CKM) mixing matrix, and conveyed only by charged currents at tree-level. While tensions are occasionally reported -see e.g. [1] and references therein for a recent example -this minimal picture seems globally consistent with the current experimental status of flavour observables [2], such that new sources or new mediators of flavour violation are relevantly constrained by these measurements. Yet, a proper confrontation of a new-physics model to such experimental results does not depend exclusively on the accuracy of the measurements or of the theoretical predictions in the a e-mail: florian.domingo@csic.es SM, but also on the magnitude of the effects induced beyond the SM (BSM).
In this paper, we will consider the well-motivated supersymmetric (SUSY) extension of the SM known as the Nextto-Minimal Supersymmetric Standard Model (NMSSM)see [3] for a review -in a minimal flavour-violating version: we will assume that the squark sector is aligned with the mass-states in the quark sector, so that, at tree-level, only charged particles convey flavour-violating effects, which are always proportional to the CKM matrix. Our aims consist in updating our previous work in [4] to the current status of flavour observables and accordingly designing a tool for a test in the flavour-sector which will be attached to the public package NMSSMTools [5][6][7][8]. Beyond ours, several projects for the study of flavour observables in the NMSSM, or more generally in SUSY extensions of the SM, have been presented in the literature: see e.g. [9][10][11][12][13][14][15][16].
Our original work dealing with B-physics in the NMSSM [4] discussed the processes B R(B → X s γ ), B R(B 0 s → μ + μ − ), B R(B + → τ + ν τ ) as well as the oscillation parameters M d,s . These processes had been implemented in the Fortran code bsg.f at (grossly) leading order (LO) in terms of the BSM contributions, using the NLO formalism for the SM and locally correcting it to account for NNLO effects in an ad-hoc fashion: in other words, this analysis essentially compiled results of the late 90'/early 2000's [17][18][19][20][21][22][23][24][25]. In doing so, it ignored existing NLO results in the MSSM [26][27][28], focussed instead, at the loop level, on tan β-enhanced Higgs-penguin contributions [22,25] and only caught the early developments of the NNLO calculation in the SM [29,30].
The SM analysis of B R(B → X s γ ) at NNLO has been recently updated in [31,32]: the corresponding results account for significant progress since [29] and shift the SM expectation ∼1σ upwards, very close to the experimental measurement. Similarly, B R(B 0 s → μ + μ − ) has been considered up to three-loop order in the SM [33][34][35], shifting the result upwards with respect to the LO. Moreover, LHCb and CMS now provide an actual measurement of this process [36], which tightens the associated constraint significantly with respect to the previous upper limits. The SM status of B → X s l + l − has also received some attention lately [37]. Finally, several other channels -e.g. B + → D ( * ) τ + ν τ , the b → sνν or the s → dνν transitions -have been suggested as complementary probes of new physics.
In addition to these recent developments concerning the SM and experimental status of flavour processes, we note that, as the NLO contributions in supersymmetric extensions of the SM can be extracted from e.g. [26][27][28], it is scientifically sound to include them into our implementation of the observables, so as to reduce the associated uncertainty in the test. This is particularly true in the case of B R(B 0 s → μ + μ − ), since this process is now measured and no longer simply bounded two orders of magnitude from above. The substantial shift in the SM estimate for B R(B → X s γ ) also tightens the constraint on BSM effects, so that enhanced precision is relevant.
Our purpose in this paper consists in describing the new status of flavour observables within the package NMSSM-Tools. In this implementation, we update the SM estimates and experimental measurements to current values. We also aim at an inclusion at NLO 1 of the SUSY contributions (whenever these results are accessible in the literature). We stress, however, that our study is restricted to a specific pattern of flavour violation -determined by the simplifying assumption in NMSSMTools of a squark sector strictly aligned with the quarks -where flavour-changing effects are only mediated by charged particles (i.e. not by gluinos and neutralinos) and in proportion to the CKM matrix. This 'minimal' version of flavour violation in a SUSY model appears as a phenomenologically justified assumption, since no compelling flavour-violating effects beyond the SM have been reported. We note however that this formulation is not stable under radiative corrections and that higher orders would restore neutral flavour-changing mediators. Moreover, as gluino vertices involve the strong coupling constant (though at high-energy), gluino corrections to LO diagrams -intervening e.g. at two-loop in the b → sγ transition -can be viewed as NLO. We still neglect such effects under the assumption that they are suppressed by the large mass of the gluinos (again a phenomenologically justified approximation). Moreover, we will generally assume that the two first generations of sfermions are degenerate. Another situation which escapes our implementation of flavour-transitions (and more generally the partonic description of quarks employed in NMSSMTools) appears in the presence of very light Higgs states with mass below ∼1 GeV, where the interplay between the Higgs state and the strongly-interacting sector is non-trivial. Finally, sources of CP-violation beyond the SM are not considered here, i.e. the NMSSM parameters are assumed real (except for the CKM matrix). Concerning the observables, we retrieve and extend the list of [4], including few rare Kaon decays -K → πνν -and oscillation parameters -M d,s , M K . On the other hand, we omit the exclusive channels B → K ( * ) l + l − , which have grown popular since the LHC Run-I due to possible hints of tensions with the SM interpretation [1]. A first reason is that a clearer understanding of the reported anomalies may be in order, a second, that the SUSY framework may not be suited to resolving them (see e.g. [38] for a discussion in the MSSM). Other transitions -e.g. K → μν -or observables could also be considered and we have no claim to exhaustivity, though we believe to have included most of the classical flavour tests for BSM physics.
In the following sections, we will first remind succinctly of the formalism employed to account for modified Higgs couplings at large tan β. We will then review briefly each observable and refer explicitly to the literature that we use in their implementation: in a first step, we shall focus on processes in the B sector before turning to Kaon physics. Finally, we will illustrate the workings of the new flavour constraints on the parameter space of the NMSSM, comparing the results of our new implementation with the former ones in a few scenarios, and discussing the relevance of the new observables that have been included. We also include a brief comparison with SuperIso [11][12][13].

tan β-Enhanced corrections to the Higgs-quark couplings
The Higgs sector of the NMSSM consists of two doublets as well as a singlet S. As in the MSSM, the tree-level couplings to quarks involve H u and H d in a Type-II 2-Higgs-Doublet-Model (2HDM) fashion: where the diagonal Yukawa parameters can be written in terms of the tree level quark masses and the Higgs vac- Yet, radiative corrections, particularly those driven by the SUSY sector, spoil this Type-II picture and generate effective terms such as -in the SU (2) × U (1)-conserving approximation: While in principle a higher-order concern, such terms may be enhanced for large values of tan β, so that a resummation becomes necessary for a consistent evaluation.
Here, as in our original work, Buras et al. [25] remains our main guide. This paper shows that the corrections to the Higgs-quark couplings driven by supersymmetric loops are well approximated in an effective SU (2) × U (1)-conserving theory. Corrections to the down-Yukawa couplings and the associated Higgs-quark vertices are dominated by the loopinduced and tan β-enhanced contributions to the H † u q L d c R operator, which, in turn, can be encoded as the corrections to the down-type quark mass-matrix m Corrections to the up-Yukawa are somewhat more subtle, as tan β-enhanced terms do not appear at the level of the quark masses, but only at that of some -e.g. charged -Higgs couplings. There, Buras et al. [25] shows that the parametrisation of [22], complemented by additional corrections, gives competitive numerical results. We will be working within these approximations.
In this framework, all the tan β-enhanced Higgs-quark vertices can be encoded in terms of the 'apparent' quark massesm f q and CKM matrix elements V e f f f f , as well as a bunch of 'ε-parameters' which parametrise: • corrections to the down-type masses: Such diagonal contributions are mediated by gluino-sdown 2 /neutralino-sdown or charginosup loops and can be extracted from Eqs. (2.5) and (A.2) in [25]. Note that, here as below, it is quite straightforward to infer the relevant NMSSM expressions from their couterparts in the MSSM, since the NMSSM only differs in the presence of singlet/singlino states, which are sterile in the couplings to SM fermions.
• off-diagonal corrections to the down-type mass matrix: ε These, in our minimal flavour violation approximation, are exclusively mediated by chargino-sup loops. Eqs. (2.5) and (A.2) of [25] again provide explicit expressions in terms of the supersymmetric spectrum.
• for the up-type couplings, ε u (d) is defined as the effective correction to the H + d u c R d L vertex (see [22]): It is computed according to Eqs. (5.6)-(5.8) of [25], i.e. including relevant electroweak-gauge effects.
• corrections to the CKM matrix elements can be encoded in The relevant flavour-changing Higgs-quark couplings are then given in Eqs. (3.55)-(3.61) and (5.8) of [25]. We note that, in this approach, the couplings of the Goldstone bosons expressed in terms of the effective quark masses and CKM elements are formally identical to the tree-level vertices expressed in terms of the tree-level masses and CKM matrix, so that the Goldstone bosons do not convey explicit tan β-enhanced terms.
Another remark addresses the explicit calculation of the ε-parameters: we neglect the Yukawa couplings of the two first generations and assume degeneracy of the corresponding sfermions. Consequently, the unitarity of the CKM matrix can be invoked in order to include the contributions of both generations at once as, e.g.,

Observables in the B-sector
In this section, we intend to present which B-physics observables are included within the new NMSSMTools flavour subroutines, and how. Let us first observe that a substantial part of the literature that we refer to below considers the MSSM, not the NMSSM: yet, as we mentioned in the previous section, generalising these results to the NMSSM offers no particular difficulty. Indeed, at the formal level, the contributions in the NMSSM differ from their analogues in the MSSM only by the inclusion of the singlet and singlino components -not directly coupled to SM fermions and gauge bosons -in the Higgs and neutralino sectors. The NMSSM thus generates no new topology in flavour-transitions (at least at NLO order) and it consequently proves sufficient to extend the sums over one CP-odd Higgs, two CP-even Higgs or four neutralino components of the MSSM to the two CP-odd Higgs, three CP-even Higgs or five neutralino mass-states of the NMSSM, whenever these particles intervene. In this sense, the impact of the NMSSM on flavour-observables with respect to the MSSM is essentially 'kinematical' in nature, meaning that in this extended model, the MSSM contributions are split and projected on various mass-states, with possibly uncon-ventional characteristics from the point of view of the MSSM (e.g. light singlet-dominated Higgs states). Additionally, the couplings involving the Higgs and neutralino states should be replaced by the correct Feynman rules in the NMSSM (see e.g. appendix A.2 of [3] for the Higgs couplings), which typically consists in promoting trigonometric functions of the MSSM mixing angles α and β to elements of the NMSSM rotation matrices in the CP-even and CP-odd Higgs sectors, S and P. We will not detail this straightforward procedure any further in the following and we will thus direct the reader to the MSSM literature without further comments. Then, we also wish to remind the reader that we work under the strict assumption that the squark sector is aligned with the quarks at the SUSY scale (defined by the squark masses), so that flavour-changing gluino or neutralino vertices vanish. This hypothesis is dictated by consistency with the general processing of squarks within NMSSMTools, which neglects any mixing of the generations at all scales. We note that the minimal flavour violation condition can be formulated in a much milder way, so that the Renormalisation Group Equations (RGE) would result in a non-alignment of quarks and squarks at the SUSY scale: such effects are necessarily beyond the scope of this work, or indeed the program it aims at completing.

B → X s(d) γ
As mentioned earlier, the status of B → X s γ in the SM has substantially evolved since the analyses of [29,30]. The new NNLO SM estimate for E 0 = 1.6 GeV [31,32] (where E 0 is the cut on the photon energy), shifted ∼1σ upwards with respect to the older estimate, is indeed very close to the experimental measurement [39] (combining results from CLEO, Belle and BABAR): Trying to account for this result by tuning the c-quark mass in the NLO formalism, as we did in [4], potentially opens new sources of uncertainty. On the other hand, employing the full NNLO formalism is an effort-consuming task of limited interest (in our position), considering that BSM effects will be included at NLO only (see below). We therefore settled for a 'middle-way', using the NNLO formalism but encoding the pure SM NNLO effects within free parameters which are numerically evaluated by comparison with the numbers provided in [31,32].
In the NNLO formalism, one can write [40]: where: • C accounts for the normalisation using semi-leptonic decays: a current estimate can be extracted from Eqs. in [40]: Beyond the Wilson coefficients at LO i j play a central role in this expression (for simplicity of notations, we factor out 2 in the case of K (1) ii ). They can be extracted from Eqs. (3.2)-(3.13) of [40] as well as Eqs. (3.1) and (6.3) of [41] and convey the NLO corrections to the partonic process b → sγ as well as the associated Bremsstrahlung contributions. All the i j (δ, z) functions entering K (1) i j are fitted numerically. The only real difficulty lies in incorporating the third line of Eq. 5, which contains the NNLO Wilson coefficient C the NNLO coefficients K (2) i j . However, if we confine to the NLO for BSM contributions, we see that these missing quantities originate purely from the SM and may be parametrised as: P only include the BSM effects. Using the numerical input from [31] -see Eqs. (6) and (10) -it will be possible to identify the coefficients P (2) 2 (SM) and Q (2) 7,8 (SM). Several other (numerically small) effects are actually absorbed within this parametrisation, though not formally NNLO. However, we include the contributions from 4-body final states (b → sqqγ ) at NLO explicitly, 3 following [42,43]. N (E 0 ) stands for non-perturbative corrections. From Eqs. (3.9) and (3.14) of [44]: where one can use the input from Appendix D -Eq. (D.1) -of [32], with the dictionary: . This procedure is aimed at parametrising phenomenologically the non-perturbative effects, the parameters being determined in a fit of the semi-leptonic B decays. [31,32] then invoke [46] to estimate the irreducible uncertainties. We come to the Wilson coefficients at the low scale. These are connected to the Wilson coefficients at the high scale via the RGE. For the LO coefficients, the solution to the RGE's -provided C (0)eff i (μ 0 ) = δ i2 , for i = 1, · · · , 6 -can be found in Appendix E of [23]. Alternatively, one may directly use [47], which allows to derive the NLO coefficients as well: with the 'magic numbers' a i , m (00) kl,i , m (10) kl,i and m (11) kl,i of Tables 1, 3, 4 in the cited reference. Finally, the QED coefficient can be obtained from Eqs. (27), (85) and (86) of [24] -multiplied with a factor α em 4π −1 -and proceeds originally from the study in Ref. [48]. 3 We thank J. Virto for his suggestion and assistance.
Having sketchily described the general formalism in the previous lines, we are left with the sole remaining task of defining the Wilson coefficients at the matching scale. From the discussion above, it should be clear that, at the considered order, we need only the C Still, this improves on the treatment in [4] where only 2HDM effects were included at NLO.
(52)-(64) of [20]. • LO contributions from chargino/stop loops were given in [21] -Eqs. (4)-(7) -but the NLO effects in Eqs. (9)-(27) (of the same reference) are not straightforward. Instead, we prefer to use [26,28]. In order to avoid the explicit ln μ 0 mT in the NLO coefficient, we take care of defining the LO coefficients at the stop (or scharm/sup) scale directly, then running it down to μ 0 via the RGE's -and taking into account the flavour-dependence in the running, i.e. the anomalous dimensions 14 23 and 16 23 for five flavours become 14 21 and 16 21 for six flavours. • Finally, for tan β-enhanced two-loop effects at the level of the Higgs-quark couplings, we no longer follow [22], Eqs. (18)-(19) -which are phrased in terms of the treelevel, and not of the apparent, parameters -but Eqs. (6.51) and (6.53) of [25], the former amounting to 0 for the G ± contribution in the SU (2) × U (1)-conserving limit. As before, effective neutral Higgs/bottom quark flavourchanging loops are included -see Eq. (6.61) in [25] and the correction for C 8 in Eq. (79) of [49].
We note that, at two-loop order, diagrams involving one gluino propagator formally contribute at NLO: such effects were considered in e.g. [50]. We shall not include such contributions under the assumption that they are suppressed in the heavy-gluino limit. At this point, the implementation at SM NNLO + BSM NLO is almost complete. The only remaining task consists in identifying the NNLO coefficients P numerically. For this, we take good care of employing the input parameters described in Appendix D of [32] and turning off the BSM contributions. To recover the branching ratio in [31] -see Eq. (6) of this reference -we determine a correction P (2) 2 (SM) of the order of 5 % of the total P(E 0 ). Then, linearising Eq. 5 in terms of LO BSM coefficients at the matching scale, we find that our implementation should be supplemented with coefficients Q (2) 7,8 (SM) at the permil level in order to recover the numbers appearing in Eq. (10) of [31]. These numbers are of the expected order of magnitude.
Let us finally comment on the error estimate. The SM + CKM + non-perturbative uncertainties have been combined in quadrature in Eq. (6) of [31] and we simply double the resulting number 0.23 in order to obtain 2σ bounds. On top of this SM + CKM + non-perturbative error, we add linearly a higher-order uncertainty of 10 % on the LO and 30 % on the NLO new-physics contributions, each type -namely 2HDM, SUSY, neutral Higgs -being added separately in absolute value. To incorporate this uncertainty, we simply use the linearisation that has been employed to determine the NNLO parameters just before.
was originally considered in [24] at NLO and then, in view of the BABAR measurement [51], by [52]. Finally, [31] extended the analysis to NNLO. Beyond the trivial substitution s → d in CKM matrix elements, the chief difference with B R[B → X s γ ] originates in sizable contributions from the partonic is not negligible. The latter can be sampled in several ways -see e.g. [53] -which provides some handle on the associated error estimate. We will be content with the evaluation using constituent quark masses given in Eq. (3.1) of [53], setting the ratio m m b -with m standing for the mass of the light quarks -in such a way as to recover, in the SM limit, the central value of [31], Eq. (8): We can then check the consistency with Eq. (10) of [31] for the new physics contributions. As before, the SM + CKM + non-perturbative uncertainties are taken over from [31], Eq. (8) -again we double the error bands to test the observable at the 2σ level -and we add linearly the new physics uncertainties. On the experimental side, the BABAR measurement [51] has to be extrapolated to the test region, leading to the estimate [52]: is the observable where the evolution since [4] has been the most critical. The experimental status has seen the upper bound replaced by an actual measurement at LHCb and CMS [36,39]: The corresponding value agrees well with the recent SM calculation [33]: It is thus no longer sufficient to consider effects driven by tan β-enhanced Higgs penguins only, and we therefore design a more complete test at NLO in the new version of bsg.f. The general formalism remains unchanged and the master formula can be recovered e.g. using Eqs. (5.15)-(5.16) of [27]: As before, we shall neglect effects from the 'mirror operators' -which are suppressed as m s /m b -and focus on the leading coefficients c A (pseudovector), c S (scalar) and c P (pseudoscalar) of the (bs)(μμ) system. The analysis is simplified by the fact that -provided the corresponding operators have been suitably normalised -these semi-leptonic coefficients have a trivial running.
is known up to three-loop QCD [34] and leading QED order [35]. It projects on the pseudovector operator exclusively. We shall use the numerical parametrisation of [33], Eq. (4), to account for it.  [27]. Instead of using Eqs.
(3.50)-(3.58) of that same reference for the neutral-Higgs penguins, we resort to [25], Eqs. (6.35) and (6.36). As in [4] -see also [54] -we replace the squared Higgs mass in the denominator by a Breit-Wigner function, so as to account for potentially light Higgs states.
The prefactor induces an uncertainty related to CKM, lattice (hadronic form factor) and B-width measurement. These are combined in quadrature at the 2σ level. In practice, we use: m B s = 5.36677 [55], τ B s = (1.607 ± 0.010) ps [39], f B s = (226 ± 6) MeV -which is an ad-hoc combination of the various results presented in [56] -and |V tb V * ts | = (41.3 ± 1.4) · 10 −3 [57]. Then, a higher-order uncertainty of 2.2 % for the SM [33] and 10 % for new physics contributions of each type are added linearly. We note that, together with the parametric (CKM, form factor) uncertainty, our error estimate in the SM limit is somewhat more conservative than that quoted in Eq. 11 and would amount to ∼8.5 % (at 1 σ ).
A similar analysis can be conducted for The experimental measurement [39]: combines the LHCb and CMS limits. The formalism is the same as for the B 0 s decay up to the trivial replacement s → d. In practice, we use the quantities m B d = 5.27958 [55], τ B d = (1.520 ± 0.004) ps [39], f B d = (188.5 ± 5.25) MeV -again an ad-hoc combination of the various results presented in [56] -and |V tb V * td | = (8.6 ± 2.8) · 10 −3 [57]. With this choice, we obtain a total (CKM-dominated) uncertainty of ∼35 % (at 1 σ ) in the SM limit. Due to larger uncertainties, one expects milder limits than in the B 0 s case.
The processB → X s l + l − was not considered in [4] but had been added in bsg.f later, including only the scalar contributions from tan β-enhanced Higgs penguins. Here we aim at a more complete analysis. The study in [37] provides a recent overview of the observables that can be extracted fromB → X s l + l − . We will confine to the branching fractions in the low -  [37], while the prefactor in Eq. (4.6) of [37] can be evaluated separately to allow for a different choice of the central values of CKM/non-perturbative contributions: we choose to take the latter from [32] since the normalisation coincides with that of B R[B → X s γ ].
The computation of the Wilson coefficients for the (chromo)-magnetic operators -C eff 7 , C eff 8 -has already been described in connection withB → X s γ : we simply run these coefficients down to the matching scale of [37], μ 0 = 120 GeV. Moreover, C 10 coincides with c A -discussed in the context of B 0 s → μ + μ − -up to a normalisation factor. Only C 9 is thus missing: it can be obtained in [28] -see Eq.
(3.6) and Appendix A of this reference. Although the lepton flavour has very little impact on C 9,10 -it intervenes only via the lepton Yukawa couplings in subleading terms -we still distinguish among C e 9,10 and C μ 9,10 . While this is ignored by [37],B → X s l + l − could also be mediated by scalar operators as shown in Eq. (2.5) of [58] note that the coefficients C Q 1,2 there coincide, up to a normalisation factor, with c S,P introduced before. Therefore, we add these contributions accordingly, estimating the integrals over m 2 l + l − numerically. To account for possibly light Higgs states, the Higgs-penguin contributions from SUSY loops are isolated in c S,P and receive denominators of the form m 2 , which are then integrated. Note that the scalar coefficients c S,P depend linearly on the lepton mass, so that they matter only in the case of the muonic final state. Finally, we come to the error estimate: the SM uncertainties (including e.g. CKM effects) are extracted from Eqs. (5.13)-(5.15) of [37]. Linearising Eqs. (B.33) and (B.36)-(B.38) (of that same reference) in terms of C BSM 7, ...,10 , we associate a 10 % uncertainty to these new-physics contributions and add it linearly. As the contributions from scalar operators is added 'by hand', we use a larger uncertainty of 30 %. The experimental values relevant for theB → X s l + l − transition are extracted from Eqs. (1.1) and (1.2) in [37].
The normalised FB asymmetryĀ F B [B → X s l + l − ] could also be implemented using the results in [37]: with the quantities H A , H T and H L explicited in Appendix B of [37]. Note that the contributions from scalar operators are suppressed as m l m b 2 [58] and may thus be neglected. However, the only experimental source (Belle) [59] chose a different binning, so that the results cannot be compared.
Beyond the inclusive decay rates, much effort has been mobilised in the study of the B → K ( * ) l + l − exclusive modes in the last few years. The full angular analysis of these modes provide two dozen independent observables [60]. Tensions with the SM estimates have been reported in some of these channels, however, leading to a substantial literature (see e.g. [1,61,62]). In this context, we choose to disregard these exclusive modes for the time being, waiting for a clearer understanding of the reported anomalies (we refer the reader to [38] for a discussion within the MSSM).

The b → sνν transition
The b → sνν transition is known to provide theoretically clean channels. While ignored in our original work, we decide to include the three following observables in the new version of the code: We follow the analysis of [63] (Section 5.9), updated in [64]. The Wilson coefficients are provided at NLO in Section 3.2 of [27]. Under our assumption of minimal flavour violation, with no flavour-changing gluinos or neutralinos, and neglecting the masses of the light quarks, only the coefficient C L (or X L in the notations of [63]) receives contributions in the model. The relation between the branching ratios in the NMSSM and that in the SM thus becomes particularly simple: see Eqs. (229-232) of [63]. We employ the updated SM evaluations in Eqs. (10), (11) and (23) of [64]: We also note that the ratio of the B + /B 0 lifetimes controls that of the B + → K ( * )+ /B 0 → K ( * )0 transitions.  [66] and BELLE [67] respectively (see also the compilation in [39]). At 90 % CL: Generalising to the b → dνν transition is trivial, though not competitive at the moment.

Flavour transitions via a charged current
The central observable in the b → u transition is B R[B + → τ + ν τ ]. Here, we perform little modification of the original implementation in [4]. In other words, we follow [68], where the effects of the W and charged-Higgs exchanges at treelevel, corrected by tan β-enhanced supersymmetric loops, appear in Eq. (5-7) of the quoted paper. The uncertainty is assumed dominated by V ub [55] and the hadronic form factor [56].
These contributions are mediated by a charged Higgs and we can easily translate, for the charged-Higgs/quark couplings, the notations of [77] to ours (see Sect. 2). We find the following Wilson coefficients: We assume a 30 % uncertainty on these new-physics coefficients, which we add linearly to the SM uncertainty quoted above. Since these observables are only marginally compatible with the SM prediction, we do not devise an actual test for them, but simply propose an evaluation.

B 0 d,s oscillation parameters M d,s
The old version of the code used the formalism of [25] to encode the SM -see Eq. As before, SM and double-penguin contributions are also included within this formalism. • We follow Section 3.1 of [79] to connect the matching scale to the low-energy matrix elements -5-flavour running. • The low-energy physics is described by the so-called 'Bag' parameters -matrix elements of the operators. We rely essentially on the lattice calculations of [80], except for the operator Q V L L , which receives the SM contributions and has thus attracted more recent attention. In this later case, we use the current FLAG average [56] for B B d,s -which coincides with the Bag parameter up to a rescaling.
These ingredients allow to derive a prediction for M d,s using the master formula of [25], Eqs. (6.6)-(6.8). The hadronic form factors f B s,d are taken from [56], where we combine the various results. For the CKM elements, we continue to rely on the evaluation from tree-level processes proposed in [57]. Note that our central value for M s in the SM limit is somewhat higher than the latest estimates [81,82]. This is essentially due to the choice for the lattice input: [81,82] have their own averaging, leading to a smaller form factor, while we follow [56]. We come to the error estimate. The uncertainty associated to the SM contributions to the operator Q V L L is often neglected in the literature. Equation (11) of [81] shows however that there could be an error of at least a few permil. We therefore associate a 1 % uncertainty to this contribution, which we add linearly to a 30 % uncertainty on each type -charged Higgs Box/SUSY Box/double penguin -of new-physics contribution. Then, the uncertainties on the Bag parameters are taken from [56,80] and combined in quadrature at the 2σ level. Finally, we factor out the uncertainties on the CKM [57] and lattice form factor [56], adding them in quadrature at the 2σ level. Note that the CKM uncertainty dominates the total error on M d (at the level of 60 %) and is actually of the order of magnitude of the central value, so that it is important not to linearise the associated error. In the SM limit, we obtain an uncertainty (at the 1 σ level) of ∼11 % for M s and ∼40 % for M d .
These results are then confronted to the experimental measurements [39]:

The s → dνν transition
The physics of Kaons also provides limits on new physics, one example being the s → dνν transition. We again follow [63] (Section 5.8), together with the updated SM results of [83]. The Wilson coefficients mediating the transition [84] are very similar to those intervening in the b → sνν transition, except that the interplay of CKM matrix elements is formally different. The normalisation to V td V * ts , instead of V tb V * ts , gives more weight to the two first generations: for completion, we thus incorporate the effects proportional to the charm Yukawa coupling. Note that we continue to neglect the quark masses of the first generation and that tree-level neutralino and gluino couplings do not mediate flavour transitions (per assumption), so that only the X L (in the notations of [63] and equivalent to the C L of Sect. 3.4) coefficient is relevant.
The decay of a charged kaon to π + νν, as well as that of the neutral K L to π 0 νν, can then be encoded in terms of this Wilson coefficient: see Eqs. (213-214) of [63], where, however, we substitute the more recent SM input of Eqs. (2.2), (2.9), (2.11) of [83].

K −K mixing
As for the B mesons, one can consider the mixing of K andK mesons. Associated quantities are very precisely measured experimentally [55]: where M K stands for the mass-difference between K L and K S and ε K measures indirect CP-violation in the K −K system. However, the theoretical evaluations of these quantities of the K −K system suffer from a substantial uncertainty associated to long-distance effects. Several estimates based on representations of large N QCD had been proposed in the 80 ies -see e.g. [87]. Lately, lattice collaborations have been emphasising the possibility to perform an evaluation in a realistic kinematical configuration in the near future: see e.g. [88]. We will follow [89] (see also discussion and literature therein) in estimating the long-distance contribution to M K at (20±10) % of the experimental value, while [90] provides some lattice input for ε K : we take over the quantity ξ 0 from Eq. (74) and the error estimate on ξ L D , Eq. (75) (of that reference) -these values were originally computed in [91] and [92], and Eq. (67) of [90] explicits their impact on ε K .
We now turn to the short-distance contributions in the K −K system. The discussion is very similar to the case of the B −B mixing in Sect. 3.6. We follow [93] for the SM part: this paper performed a NNLO evaluation of the charm contribution -see Eq. (15) in that work -completing earlier results for the mixed charm-top [94] and top [95] contributions. Note that [89,90] also propose recent evaluations of the quantity η cc , with slightly lower central value, but we choose to stick to the conservative estimate of [93]. The master formulae for M K and ε K are provided in Eqs. (18) and (16) of this reference, respectively -see also Eqs. (XVIII. [6][7][8][9] of [17], as well as Eqs. (XII. [3][4][5] for the expression of the functions S 0 . The kaon mass of 0.4976 GeV and the form factor f K = 0.1563 ± 0.0009 GeV are taken from [55] and [56] respectively. The inclusion of BSM contributions follows the general NLO formalism of [79]: see Eqs. (7.24-7.32). Eqs. (3.20-3.38) (of this reference) explicit the running between the matching-and the low-energy scales. However, we will be using more recent 'Bag parameters' for a low-energy scale of 3 GeV: Table XIII of [96] compiles several recent lattice calculations, which we put to use. As in the B −B case, the Wilson coefficients account for Higgs/quark and chargino/squark box diagrams as well as Higgs doublepenguin contributions and we follow appendix A of [25]. Yet, we also include effects associated to the charm Yukawa, as the interplay of CKM elements gives more weight to such terms than for the B −B mixing. Finally, we use appendix C of [79] to run each new-physics contribution from the relevant BSM scale (charged-Higgs or squark mass) down to the matching scale of 166 GeV.
The SM uncertainty, where the uncertainty on the η cc parameter and the long distance effects dominate, is added linearly to the uncertainty driven by the bag parameters and a 30 % error on higher-order contributions to the BSM Wilson coefficients. Moving to the SM limit and considering the theoretical error at 1 σ , we obtain a total uncertainty of about 35 % for M K and 15 % for ε K .
Note that one can lead a similar analysis for the D −D mixing, which we do not consider here, however.

Sampling the impact of the flavour constraints
Based on the discussion of the previous sections, we design two Fortran subroutines bsg.f and Kphys.f for the evaluation in the NMSSM of the considered observables in the B-and the Kaon-sectors (respectively), as well as a confrontation to experimental results. These subroutines are then attached to the public tool NMSSMTools [5][6][7][8].

Comparing the new and the old codes
In order to test the differences between the new and the former implementations of B-physics observables in bsg.f, we perform a scan over the plane defined by m H ± -the charged Higgs mass -and tan β and display the exclusion contours associated with flavour constraints in Figs. 1 and 2. The chosen region in the NMSSM parameter space corresponds to the MSSM-limit, with degenerate sfermions and hierarchical neutralinos. Note that we disregard the phenomenological limits from other sectors (e.g. Higgs physics, Dark matter, etc.), which are largely orthogonal to the flavour-constraints and would unnecessarily complicate the picture. We consider a large value of the trilinear stop coupling |A t | = 2.5 TeV, which is known to enhance effects driven by supersymmetric loops, and study separately the two opposite signs -a negative value of A t , when μ > 0, typically triggers destructive interferences among the SUSY and 2HDM contributions tō B → X s γ .
The general appearance of the exclusion contours in Figs. 1 and 2 remains qualitatively similar, when comparing the results obtained with the new (plots on the top) and old (plots on the bottom) versions of the code. 4 Yet, quantitatively, one witnesses a few deviations: • The limits fromB → X s γ are more severe in the new version, which is mostly apparent in Fig. 1 (A t > 0): this is not unexpected since the larger SM central value -closer to the experimental measurement -correspondingly disfavours new physics effects that interfere constructively with the SM contribution (2HDM effects or supersymmetric loops for μ, A t > 0). Consequently, the areas with a light charged Higgs or large tan β receive excessive BSM contributions in view of the experimental measurement and are thus disfavoured. Moreover, note that the NLO implementation reduces somewhat the error bar associated to higher-order new-physics contributions, which also results in tighter bounds for the more recent code. For A t < 0, one observes two separate exclusion regions: for low values of tan β and m H ± , the 2HDM contribution is large (excessive) while the negative SUSY effect is too small to balance it; on the contrary, with large tan β and heavy H ± , the SUSY contribution dominates and is responsible for the mismatch with the experimental measurement. In between, the destructive interplay between the SM and 2HDM effects on one side and the SUSY loops on the other succeeds in keeping B R[B → X s γ ] within phenomenologically acceptable values. • Limits from B 0 s → μ + μ − used to be little sensitive to the sign of A t in the older implementation. This is no longer true, the reason being that the scalar coefficients c S,P receive new contributions, namely 2HDMand SUSY-mediated Z -penguins and boxes, which (had been neglected in the previous version of the code and) may interfere constructively or destructively with the Higgs-penguin effects. This channel appears as the most sensitive one, together withB → X s γ , in the considered scenario. Given the shape of the exclusion regions driven byB → X s γ , however, B 0 s → μ + μ − seems Points excluded by the various constraints are coloured in red (horizontal lines) -B → X s γblue (vertical lines) -B 0 s → μ + μ −dark green (circles) -B → X s l + l −light green (oblique lines) -B + → τ + ν τyellow (diamonds) -M s -and orange (triangles) -M d -while remaining gray points satisfy all these limits most relevant for A t < 0 (Fig. 2). Expectedly, the limits are tighter for large tan β, where SUSY contributions are enhanced.
• Limits fromB → X s l + l − differ more significantly between the two implementations -although they remain subleading. In particular, an excluded region appears at low tan β: it is largely driven by the 2HDM contributions to the semi-leptonic vector coefficients C 9,10 -which indeed involve terms in tan −1 β. On the other hands, the exclusion region at low m H ± is largely unchanged: it is associated with the enhancement of the Higgs-penguin contributions for a light Higgs sector. • Despite the corrections to the tan β-enhanced Higgs/quark vertices, the constraints from B + → τ + ν τ , M s and M d are little affected by the modernisation of the code and remain subleading.
We observe thatB → X s γ and B 0 s → μ + μ − intervene as the determining limits from the flavour sector in the considered scenario: they exclude all the region beyond tan β > ∼ 20 (for our specific choice of parameters in Figs. 1, 2). The low m H ± -region is in tension with most of the observables in the B-sector (unsurprisingly), thoughB → X s γ and B 0 s → μ + μ − again appear as the limiting factors at low-tomoderate tan β. Interestingly,B → X s l + l − seems to offer a competitive test for tan β < ∼ 2.
We perform a second test in a region involving a light CPodd Higgs state with mass below 15 GeV -still presuming nothing of the limits from other sectors: note that this is a phenomenologically viable scenario in the NMSSM, although the limits on unconventional decays of the SM-like Higgs state at ∼125 GeV place severe constraints on the properties of the light pseudoscalar. The results are displayed in Fig. 3 Fig. 1, but with A t = −2.5 TeV -in terms of the mass of the pseudoscalar m A 1 and tan βand confirm the trends that we signaled before:

Fig. 2 Similar plots and scan as in
• Limits fromB → X s γ intervene here at low tan βwhere the supersymmetric contributions cannot balance the effect triggered by the charged-Higgs (note that A t < 0). A few points are also excluded for low m A 1 and large tan β: these result from the two-loop effect mediated by a neutral Higgs. They prove subleading in the considered region. • Limits from B 0 s → μ + μ − appear somewhat tighter in the new implementation. In particular, a narrow corridor where the new physics effects reverse the SM contribution is visible in the plot on the bottom of Fig. 3 (which corresponds to the older implementation of the limits)from (m A 1 ∼ 6, tan β ∼ 2) to (m A 1 ∼ 15, tan β ∼ 5); this region is no longer accessible with the more recent code (it is, in fact, shifted to lower values of tan β). This channel is the main flavour limit in the considered region, due to the large contribution mediated by an almost onshell Higgs penguin.
• Limits fromB → X s l + l − intervene in two fashions. One is the exclusion driven by an almost-resonant pseudoscalar and the associated bounds are essentially unchanged with respect to the older implementation. Additionally, a new excluded area appears at low tan β. • Limits from M d,s are qualitatively unchanged among the two versions, though the bounds associated with M d seem somewhat more conservative in the new implementation. These constraints remain subleading however, in view of the more efficient B R(B 0 s → μ + μ − ), and confine to the resonant regime -note e.g. the allowed 'corridor' where new-physics contributions reverse the SM effect -or the very-low range tan β < ∼ 1.
B 0 s → μ + μ − thus appears as the constraint that is most sensitive to the enhancement-effect related to a near-resonant below, the results of the old version. The colour code remains the same as before. The symbols are also unchanged, except forB → X s l + l − , with horizontal lines, and,B → X s γ , with circles (for reasons of visibility) pseudoscalar. The exclusion effects are most severe for larger tan β as the Higgs-penguin is correspondingly enhanced. For tan β < ∼ 2,B → X s l + l − proves a sensitive probe in its new implementation.
Note that, in the two scenarios that we discussed here, the precise limits on the {m H ± , tan β} or {tan β, m A 1 } planes of course depend on the details of the parameters. In particular, the large value of |A t | triggers enhanced SUSY effects, resulting in severe bounds on the considered planes. We thus warn the reader against over-interpreting the impression that only corners of the parameter space of the NMSSM are in a position to satisfy B-constraints at 95 % CL, as Figs. 1, 2 and 3 might lead one to believe. To counteract this effect, we present in Fig. 4 the limits from flavour processes obtained with the new implementation, for A t = 500 GeV and a somewhat heavier chargino/neutralino sector. The plot on the top again considers the plane {m H ± , tan β} in the MSSM limit. SUSY contributions are suppressed by the choice of low A t .
Correspondingly, limits fromB → X s γ only intervene in the region with low m H ± < 300 GeV. The constraints driven by B 0 s → μ + μ − eventually exclude the large tan β > ∼ 45 range but are obviously weaker than before. On the other hand, the exclusion contour associated with B + → τ + ν τ andB → X s l + l − remain largely unaffected. The plot on the bottom part of Fig. 4 addresses the scenario with a light pseudoscalar: contrarily to the case of Fig. 3, CP-odd masses above m A 1 ∼6 GeV are left unconstrained by the flavour test, with exclusions intervening only at very-low tan β or for m A 1 in the immediate vicinity of a resonant energy (for B 0 s → μ + μ − , M s orB → X s l + l − ).

Impact of the new flavour tests
Beyond the observables that had been considered in [4], we have extended our analysis to several new channels. We now wish to discuss their impact on the NMSSM parameter space. In Fig. 5, we consider the scenario of Fig. 1 once more and present the exclusion limits driven by the newly implemented channels. Note that the constraints considered in the previous section form the black exclusion zone on the background. The limits from the various channels shown in this plane seem to be essentially subleading in view of these previous constraints of Fig. 1.
• Limits fromB → X d γ intervene essentially for A t > 0 (i.e. for constructive SUSY contributions), large tan β (driving large SUSY contributions) and light H ± (driving large 2HDM contributions). Yet the corresponding bounds are superseded byB → X s γ . • B 0 d → μ + μ − intervenes in the large tan β/low m H ± corner as well, e.g. for A t < 0, but seems less sensitive than B 0 s → μ + μ − , except in the low tan β < ∼ 2 region. • The processes of the b → sνν and s → dνν transitions are found to be well under the current experimental upper bounds.
• The K −K mixing excludes a few points (driven by ε K where the SM is slightly off, with respect to the experimental results) but is not competitive in view of the, admittedly conservative, uncertainties.
Note that the limits induced by the b → cτ ν τ channels have been omitted in Fig. 5. Given the current data, this transition would exclude the whole {m H ± , tan β} plane, with the exception of the large tan β/low m H ± corner -which is excluded by most of the other flavour constraints: the significant discrepancy of the SM estimate with the experimental measurement, especially for B → D * τ ν τ , explains this broad exclusion range. SUSY 2HDM effects cannot reduce the gap much, except in already excluded regions of the parameter space.
Then, we return to the light pseudoscalar scenario of Fig. 3 and display the constraints associated with the new channels in Fig. 6. Again, these limits are found to be weaker than those shown in the previous section. Limits from B 0 d → μ + μ −  prove the most constraining of the new channels in this regime: this again results from the enhancement of the Higgspenguin mediated by a resonant A 1 . Subleading constraints from the K −K mixing also intervene at low tan β < ∼ 1 and for very light CP-odd Higgs with m A 1 < ∼ 2 GeV. Again, the discrepancy among SM predictions and experimental measurements for the b → cτ ν τ transition cannot be interpreted in this scenario, so that applying a 95 % CL test for the ratios R D ( * ) would lead to the exclusion of the whole portion of parameter space displayed in Fig. 6.
Finally, we complete this discussion by considering the parameter sets of Fig. 4, where the flavour limits discussed in the previous section were found weaker. The impact of the new channels can be read in Fig. 7. The corresponding exclusion regions in the considered regime with A t = 500 GeV again prove narrower than those considered in Fig. 4. (Note again that we have omitted the b → cτ ν τ channels, however.) Therefore, we find that the new channels tested in bsg.f and Kphys.f are typically less constraining than the older ones, which we discussed before. Limits fromB → X d γ and B 0 d → μ + μ − are found to be significant, however, and an evolution of the experimental limits or an improvement in understanding the SM uncertainties may provide them with more relevance in the future. The b → cτ ν τ transition stands apart as the tension between SM and experiment resists an NMSSM interpretation, at least in the scenarios that we have been considering here.
As a final remark, let us mention that the flavour routines that we discussed are not designed to address very light CPodd states, below the threshold for decays to pions or even muons, as e.g. [97] considered -note that the decay widths computed within NMSSMTools are not reliable for masses below ∼1 GeV and that the interplay of the light Higgs with the hadronic sector is non-trivial. In such cases, the light pseudoscalar could be long-lived and decay outside of the detectors, hence behave like an invisible final state. Moreover large branching fractions for the decay A 1 → e + e − appear at low mass and could induce sizable contributions to the b → se + e − and s → de + e − transitions: this will not be properly tested.

A short comparison with SuperIso
SuperIso [11][12][13] is a public tool devoted to the calculation of flavour observables in a number of models, among which the NMSSM. The list of channels that are considered there differs somewhat from ours, including e.g. observables in the and B R(B + → Dτ + ν τ )/B R(B + → Dl + ν l ) -coincide with those considered in the present work, so that we may attempt a brief (and necessarily limited) comparison. For this, we will employ the current version of SuperIso -i.e. 3.5 -and feed it the slha output of NMSSMTools as input for several points (to ensure that we perform the comparison for exactly the same spectra).
In Fig. 8, we consider the MSSM-like scenario of Fig. 1, where we fix M A = 800 GeV and restrict to eight values of tan β. The six common observables of bsg.f and SuperIso are displayed, with our estimate appearing as a black dot accompanied by its 2σ theoretical uncertainty; the central value delivered by SuperIso (to our knowledge, SuperIso does not provide an error estimate) is shown as red diamonds and the experimental measurements is plotted in green in the background (for reference). Let us first emphasize that the theoretical estimates roughly agree, in the sense that the central value computed by SuperIso generally falls within the uncertainty range of our prediction. Furthermore, the variations of the central values share a common trend (i.e. the estimates of both codes rise or go down simultaneously) for all observables and the results coincide almost exactly in the case of B R(B → X s l + l − ) 1 GeV 2 <m 2 l + l − <6 GeV 2 , B R(B → X s l + l − ) | m 2 l + l − >14.4 GeV 2 and B R(B + → Dτ + ν τ )/B R(B + → Dl + ν l ).
Looking more closely, however, we observe some discrepancies. Their origins are potentially multiple -use of distinct experimental (e.g. CKM) or lattice input, inclusion of the contributions at different orders, diverging treatment of SUSY corrections to the Higgs/down quark system, etc. Actually, given the large parametric uncertainties of the input, the results may differ without any inconsistency. Characterizing these discrepancies with accuracy goes beyond the scope of this paper or indeed the skills of its author, since a detailed comparison would require access to intermediate quantities in SuperIso. Here we will merely state what differences appear and propose an interpretation based on what superficial means stay at our disposal.
Let us first consider B R(B → X s γ ), where the two predictions coincide at low tan β but tend to diverge apart at large tan β, resulting in only a marginal agreement at tan β = 40. This effect develops far from the phenomenologically relevant regime -the predicted B R(B → X s γ ) cannot be reconciled with its experimental value then -and the question of concordance between the two codes may thus appear of limited importance. Still, we notice in the documentation of SuperIso that large-tan β corrections to the SM contributions for the Wilson coefficients C 7,8 -corresponding to the re-definition of the Goldstone boson couplings to down-type quarks, subsequent to the inclusion of radiative corrections -are included in this tool, according to Eq. (18) of [22] or Eq. (B.72) of [9]. We checked that a naive implementation of this contribution would account numerically for the difference between our results and those of SuperIso in Fig. 8 -although this could also be a mere coincidence. Including this contribution would be inappropriate, at least in our case. Indeed, as was shown in [25] and as we reminded in Sect. 2, such contributions are only consistent with an uncorrected description, at LO, of the Higgs/down-quark system -the b-mass and CKM elements would then be pure lagrangian parameters which do not coincide with the measured b mass or CKM elements.
For B R(B 0 s → μ + μ − ), the small deviations between the predictions of SuperIso and bsg.f are more difficult to trace. At low tan β, they simply result from distinct choices for the estimate of the SM contribution -see Eq. 11. Then another likely contributor is the divergent treatment of Higgs-penguin contributions: SuperIso employs the NLO expansion of [27] while we rely exclusively on the effective lagrangian (at LO) describing the Higgs / down-type couplings in terms of the physical quark masses (see Sect. 2). Finally, for B R(B + → τ + ν τ ), the central values are almost parallel and it is easy to identify the source of the discrepancy as the choice of distinct values of V ub : as such, this difference is accounted for in the parametric uncertainty that we implemented.
In Fig. 9, we turn to the low-m A 1 scenario of Fig. 3, setting tan β = 6 and choosing five values of m A 1 . The results for B R(B → X s γ ) and B R(B + → Dτ + ν τ )/B R(B + → Dl + ν l ) look very compatible between the two codes. Again, the deviation for B R(B + → τ + ν τ ) -within theoretical uncertainties -can be interpreted in terms of the choice of V ub . On the other hand, for B R(B 0 s → μ + μ − ), the two codes only agree as far as concluding that the points give a branching fraction far beyond the experimentally acceptable range. The two estimates actually differ by orders of magnitude. First, we note that, so far from the physical region (where BSM effects should remain small), the question of explaining this difference could appear as largely academical. Then, a fair guess would assign this discrepancy to the divergent descriptions of the Higgs penguin contributions. More precisely, SuperIso employs the formalism of [98] in the presence of a light pseudoscalar. Instead of charac- The output of SuperIso is plotted as red diamonds. The experimental measurement with 2σ error bars is also shown in the background, as green lines terising the pseudoscalar -and in particular its couplingsstrictly in terms of the associated rotation matrix, the latter reference directly injects lagrangian parameters (using treelevel relations; e.g. via the quantity δ − ) into the expressions of the Wilson coefficients. While receivable in a description of the NMSSM Higgs sector at tree-level, these formulae will not properly account for a radiatively-corrected Higgs sector (as e.g. delivered by NMSSMTools): the properties of a light state may vary sizably in this fashion. Addition-ally, the conventions of [98] differ somewhat from those usually employed in the NMSSM -particularly a factor √ 2 for the singlet v.e.v. We checked that a naive use of the expressions in [98] could easily account for a factor O(10) at the level of the Wilson coefficients, hence might be the source of the disparity of the estimates performed by SuperIso and bsg.f for B R(B 0 s → μ + μ − ) in Fig. 9. If so, we believe that our description of the Higgs penguins is better suited, as it properly implements the properties of the loop-corrected Higgs sector. The predictions for B R(B → X s l + l − ) in the low and high m 2 l + l − ranges should be considered under the same light: our estimates only depart from the experimentally allowed region when the pseudoscalar falls under the tested range in m 2 l + l − -i.e. for m A 1 = 2 GeV in the low range and m A 1 = 4 GeV in the high range -which seems physically intuitive and is consistent with our relatively mild excess in B 0 s → μ + μ − (the effective A 1 − b − s coupling is not extremely large in our approach). On the other hand, considering the predictions of SuperIso, the presence of a pole in the integral over m 2 l + l − does not appear to make much of a difference and the B R(B → X s l + l − ) are invariably excessive (the effective A 1 −b−s coupling computed by SuperIso is much larger than the one in bsg.f as could already be inferred from B 0 s → μ + μ − ). Another difference is that our result is averaged over l = e and l = μ final states, while SuperIso only considers l = μ: this is relevant in view of the Higgs couplings being related to m l .
We shall close the comparison with SuperIso here. We found a qualitative agreement between the estimates of this code and ours in most cases. As for the discrepancies, we can understand them most of the time, either in terms of different choices of the parameters -for B R(B + → τ + ν τ ) -or divergent processing of the effective Higgs/down quark couplings -for B R(B → X s γ ) or B R(B 0 s → μ + μ − ). In the latter case, we believe that we have employed a consistent formalism. We also stress that these effects appear for large new-physics contributions, i.e. far from the phenomenologically relevant regime. Yet a full understanding of these differences would require a much more detailed analysis and goes beyond our ambitions for this paper.

Conclusions
We have considered a set of flavour observables in the NMSSM, updating and extending our former analysis in [4]. These channels have been implemented in a pair of Fortran subroutines, which allow for both the evaluation of the observables in the NMSSM and confrontation with the current experimental results. We have taken into account the recent upgrades of the SM status of e.g. B R[B → X s γ ] or B R[B 0 s → μ + μ − ] and included BSM effects at NLO. The tools thus designed will be/have been partially made public within the package NMSSMTools [5][6][7][8].
We observe that the bounds on the NMSSM parameter space driven by B R[B → X s γ ] or B R[B 0 s → μ + μ − ] have become more efficient, which should be considered in the light of the recent evolution of the SM status and/or the experimental measurement for both these channels. In particular, the large tan β region is rapidly subject to constraints originating from the flavour sector. Similarly, the light pseudoscalar scenario is tightly corseted due to the efficiency of Higgs-penguins in the presence of such a light mediator.
Among the new channels that we have included, we note the specific status of the b → cτ ν τ transition, where the discrepancy between SM and experiment seems difficult to address in a SUSY context.
Other channels of the flavour-changing sector may prove interesting to include in the future. Note e.g. the current evolution in the B → K ( * ) l + l − observables.