# Towards a viable scalar interpretation of

###### Abstract

Recent measurements of semileptonic B-meson decays seemingly imply violations of lepton flavor universality beyond the Standard Model predictions. With three-level explanations based on extended Higgs sectors being strongly challenged by the measurements of the lifetime, new theories invoking leptoquark or vector fields appear as the only feasible answer. However, in this work we show that simple scalar extensions of the Standard Model still offer a possible solution to the physics puzzle, owing to sizeable loop-level corrections which mimic the effects of new vector contributions. Considering a simplified model characterised by a charged and a neutral scalar particle, we verify the compatibility of the observed signal with the relevant collider bounds. We also study an embedding of the simplified model into a three-Higgs-doublet framework, and investigate its main phenomenological constraints.

^{†}

^{†}affiliationtext: Laboratory of High Energy and Computational Physics, National Institute of Chemical Physics and Biophysics, Rävala pst. 10, 10143 Tallinn, Estonia.

## 1 Introduction

Recent results of the LHC experiment [1, 2] highlight a significant amount of tension between the measured properties of the -meson and their Standard Model (SM) predictions, confirming previous investigations by the BaBar [3, 4] and Belle [5, 6, 7] collaborations. In more detail, anomalies of lepton flavour universality have been measured for the underlying transition in the ratio of branching fractions

(1) |

for . Including the most recent LHC measurement, the Heavy Flavour Average Group (HFAG) reports the world averages [8]

(2) |

which exceed the SM predictions

(3) |

with a combined significance of about . Being defined as ratios, and are quantities particularly suitable for investigating lepton universality regardless of hadronic uncertainties in the involved SM parameters. All independent probes of transitions have exhibited an increased affinity for third generation leptons, resulting in a departure from universality that exceeds the SM predictions. The present situation is summarised in Fig. 1, where the elliptical areas shaded in grey represent different joint confidence regions for the mentioned observables. The horizontal band indicates instead the range probed by the latest LHC result. The SM prediction, marked by a red dot, falls just on the outside of the 4 confidence region.

New particles could clearly induce extra contributions to , however most of the proposed interpretations are immediately falsified by correlated measurements in different experiments. For instance, scalar extensions of the SM that seek to simultaneously explain through a tree-level contribution inevitably leads to an enhancement of the decay , which is incompatible with the parameter space indicated by the anomaly [9, 10, 11, 12, 13, 14]. Once the strictest experimental bounds are considered, this incompatibility apparently rules out such interpretations of the signal for extended Higgs sector models, opening the way to alternative solutions involving leptoquark fields [15], extra-dimensions [16], gauge extensions [17, 18, 19] or specific supersymmetric models [20].

However, in this paper we demonstrate that scalar extensions of the SM are still viable and can relax, if not completely solve, the tensions due to the -physics measurements. In essence, we propose the generation of effective operators typically associated with vector currents via scalar mediators at the loop level. Our results indicate that if such loop effects are considered on top of the SM contribution, the anomaly can be explained via a sizeable but perturbative correction.

The paper is organized as follows. In the upcoming section, we briefly review the operator basis used in the literature to discuss flavour anomalies. After that, in Sec. 3 we introduce the simplified model and investigate its possible role in the measured flavour anomalies. In addition, we discuss the effects of the bounds cast by measurements of the and lifetimes, as well as the collider signatures of the proposed framework. In Sec. 4 we identify a possible high-energy completion and review the main phenomenological constraints implied. Finally, we summarise our findings in Sec. 5.

## 2 Flavour anomalies and Wilson coefficients

The low-energy behaviour investigated at LHC and in -factories can be described by matching the framework at hand with an effective theory at the -quark mass-scale. The procedure results in a class of dimension six operators^{1}^{1}1Tensor operators play no role in our investigation and will therefore be neglected for sake of simplicity. that preserve color and electric charge, with coefficients to be determined from the dynamics of the considered model. Including right handed neutrinos, the effective Lagrangian for the process takes the form [21, 22, 23, 24]

(4) | |||||

where the eight operators are given by

(5) | |||||

(6) | |||||

(7) | |||||

(8) |

for and with being the corresponding chirality projector.

As mentioned before, the observables are defined as the ratios of the total branching fractions for

(9) |

where the meson branching ratios are functions of the coefficients of the adopted operator basis: .

In order to illustrate the reach of the approach, we schematise in Table 1 the dependence of selected flavour observables on the Wilson coefficient of the effective Lagrangian. The reader is referred to Refs [24, 25] for the exact expressions and the relative hadronic form factors.

Yes | No | Yes | No | Yes | No | Yes | No | |

No | Yes | No | Yes | No | Yes | No | Yes | |

No | Yes | No | Yes | No | Yes | No | Yes |

The case of is shown in more detail in Fig. 1. Clearly, non-vanishing values of the combination - (indicated by a purple solid line) allow to ameliorate the agreement with observations. Notice that fitting the anomaly involves only modest deviations of the relevant Wilson coefficient from the corresponding SM values =.

Traditional scalar extensions of the SM yield tree-level contributions to and , and, provided right handed neutrinos are considered, to their right-handed counterpart and . These coefficients are however severely constrained by measurements of the lifetime. In particular, the dependence of the branching ratio

(10) |

for on is enhanced by the mass ratio. The impact of the lifetime measurements is illustrated in more detail in Fig. 1(a), where the red, orange and yellow areas show the regions for the joint fit of and . The shaded areas indicate the values of the Wilson coefficients causing deviations larger than (light gray) or (dark gray) from the measured lifetime. As we can see, the latter is fully consistent with the SM prediction, = 0, and precludes explanations of the anomaly based on the operator.

Measurements of the decay spectrum provide instead a complementary way to probe the simplest scalar extensions [4]. As shown in Fig. 1(b), this observable constrains the scalar Wilson coefficient but is essentially independent of the remaining operators in the effective Lagrangian (4). Hence, models that for instance induce deviations only in the coefficient result in predictions for indistinguishable from the SM one.

Given the importance of the constraints cast by the measurements of the lifetime and by the distribution of , in the following we propose an effective scalar extension of the SM that explains the anomalies by generating sizeable contributions to and keeps the coefficients and at their SM values. In our analyses we refer to the 10% limit on the lifetime when assessing the viability of the proposed scheme. To conclude the section, we mention that angular and differential distributions have shown further concrete possibilities to rule out some of the solutions of the puzzle [26].

## 3 A simplified scalar model for flavour anomalies

As mentioned in the previous section, scalar extensions of the SM with sizeable tree-level contributions cannot account for the current flavour anomalies. In order to investigate the effect of one-loop processes in the B-meson decay, we then focus on the following simplified model

(11) | |||||

which introduces both an electrically neutral and a charged complex scalar field. The interactions contained in Eq. (11) have a non-zero projection on the vector and axial vector operators and , providing new contributions to the relative Wilson coefficients via the diagrams shown in Fig. 3:

(12) | |||

(13) |

The model includes also a sterile right-handed neutrino that plays an active role in the decay of the meson. We assume for our calculation that the masses of these particles are negligible but still large enough to evade the cosmological bounds on additional relativistic species. We remark that the results we present are independent of the details of the neutrino mass generation mechanism as long as the interactions of are not significantly diminished by mixing effects.

Given the structure of the amplitudes for the diagrams in Fig. 3 and the overall normalization of the SM contribution proportional to , we can produce percent-level variations in for new scalar fields with masses in the range GeV and perturbative values of the couplings. Larger values would require the loss of perturbativity when fitting the anomaly.

In order to assess the viability of the simplified model at hand, we investigate its collider phenomenology in conjunction with the signal produced from meson decays. After that, we discuss the details of a possible high-energy completion of the framework and the potential additional constraints it entails.

### 3.1 Collider phenomenology of the simplified model

The proposed model contains neutral and charged scalars, which can be singly produced at the LHC in association with top quarks [27]. The initial state for such topologies contains and quarks, which are suppressed by their small contributions to the proton parton distribution functions at large momentum fractions. Because of this, the production cross section for single scalars and top quarks falls off steeply with increasing scalar mass. At the same time, however, the events become more visible as the transverse momentum of the decay products, as well as the missing energy from neutrinos, increases. In order to derive experimental limits for our model, we therefore scan the relevant parameter space to estimate the excluded region.

In addition to single scalar/mono-top events, the charged scalars can be pair produced from quark/antiquark initial states through s-channel photon and t-channel top quark diagrams. The neutral scalars can also be pair produced via the same t-channel topology, resulting in particular in and final states, the latter being identifiable in mono-jet searches. Because the final states involve tau leptons, the suppressed electroweak processes are difficult to identify in the strongly interacting hadronic environment of the LHC. We therefore included the pair production events in the general scan over the model parameters.

The topologies considered for the production of the charged and neutral scalars are shown in Fig. 4. We decided to disregard the hadronic decay channels of the scalar bosons since events with decays to top quarks are kinematically suppressed, and decays to charm quarks cannot be identified over the QCD background. The analysed final states are reported in Table 2.

Charged scalar field | Neutral scalar field | |
---|---|---|

Single production | , | , |

Double production | , |

In order to derive the bounds presented below, we implemented the Lagrangian for the proposed simplified model with FeynRules [28] and exported the resulting model files to MadGraph5_aMC@NLO [29] via the UFO interface. For each of the topologies, we then performed a grid scan of the parameter space via Monte-Carlo simulations, applying the default kinematic cuts and collecting events per point. Finally, the partonic LHE files generated were compared against the 8 TeV analysis available in CheckMATE 2.0.14 [30]. In addition, we seperately included the constraints presented in the dedicated CMS search [31] for the single production of a charged Higgs boson.

Based on our scan of the parameter space, we now present three benchmark points.

#### Benchmark 1: GeV GeV

In Fig. 5 we show the bounds that LHC searches impose once we assume a mass of 200 GeV for both the charged and neutrals scalar fields of the simplified model. As we can see in the top panels, the most stringent bounds on the couplings of are due to the single production of the particle in association with a top quark. The allowed leptonic final states result in comparable 95% exclusion regions, which force the transverse coupling to and below unity. The double production of results in a weaker constraint that is only marginally relevant once the decay channel is predominant. The bottom panels refer instead to the charged scalars . Also in this case we find that the single production in concomitance with an (anti)top quark casts the strongest bound, limiting the value of the coupling to and quarks below about 0.7. In these plots we show also the confidence regions designated by the measurements of , where we assume the values of the neutral scalar couplings indicated by a red star in the top panels.

#### Benchmark 2: GeV GeV

Considering larger values for the scalar masses, the constraints in the top panels of Fig. 5 significantly relax and the simplified model in Eq. (11) is compatible with the measured values of . As shown in Fig. 6, while the properties of the charged scalar are still constrained by the dedicated searches, the current physics measurements can be explained for values of the couplings smaller than 2. For the computation of the confidence regions we here assumed .

#### Benchmark 3: GeV GeV

## 4 A possible high-energy completion

We investigate now a possible high-energy completion of our simplified model for flavour anomalies proposed in Eq. (11). As a first attempt, we could identify the new scalar fields with the components of an doublet and immediately recover the invariance of the Lagrangian under the full SM gauge group. In this case, however, the SM quark mixing induced by the CKM matrix inevitably results in the dilution of the quark couplings to either the charged or the neutral scalar fields, diminishing the net contribution of the diagrams in Fig. 3. This problem can be circumvented by assigning the new scalars to two different doublets, effectively embedding the framework into a three-Higgs-doublets model (3HDM). Proceeding along the lines of Ref. [32], we therefore consider

(14) |

where plays the role of the SM Higgs doublet and is solely responsible for spontaneous symmetry breaking with a VEV of GeV. The doublets instead provide the new scalar degrees of freedom which interact sizeably with the top quark, lepton and neutrinos.

In order to embed the Lagrangian of Eq. (11) in a 3HDM, we need to address the presence of two additional complex fields that can cause further conflicts with the LHC and low energy measurements, even though not directly involved in the process. To estimate the mass scale of these states, we consider the potential

(15) | |||||

resulting in

(16) |

The mass-splitting between the neutral and charged elements of a doublet is controlled by , therefore the mass degeneracy can be broken by a factor of about .

The interactions of the new states follow form the Lagrangian

(17) |

where we identified the lightest components of each doublet with the scalar fields used in our simplified model: , . The interactions that characterise the latter can then be recovered by appropriately choosing Yukawa couplings of the form

(18) |

written on the basis of the quark mass eigenstates. After having performed the CKM rotation, the scalars and couple then as

(19) |

matching the interactions in Eq. (11). The heavier components of the doublets interact instead according to

(20) |

presenting suppressed couplings with respect to those of and .

Although the collider signatures of the heavy components are thus less prominent than those of their light counterparts^{2}^{2}2The presence of additional gauge interaction does not worsen the bounds found in the previous section, which are mainly driven by single production processes that gauge interactions do not enhance. , the new couplings contained in Eq. (20) potentially induce flavour changing neutral currents and source the process. In the following we investigate the constraints that the corresponding searches yield, and use these results to bound the mass spectrum of the full model.

### 4.1

The charged current interactions in Eq. (20) induce new couplings of the top quark to the down-type quarks, resulting in a non-zero contribution to the photon and gluon magnetic dipole operators^{3}^{3}3We neglect the operators and with opposite chirality because of the overall suppression factor.

(21) |

By matching Eq. (20) with the effective Lagrangian [33]

(22) |

and requiring that the new physics contributions fall in the range [34, 35]

(23) |

we obtain the exclusion bound on the coupling and on the mass of shown in Fig. 8.

We remark that similar to , the interactions of with the up-type quarks involve only the top quark. Any sizeable contribution to due to the exchange of light quarks is therefore absent.

### 4.2 Meson mixing

Measurements of the mass splitting in the neutral system provide bounds on and the mass of which may compete with the constraint from . In regard of this, the interactions in Eq. (17) contribute to via the two one-loop diagrams presented in Fig. 9. The value of the corresponding amplitude is to be matched with the effective Lagrangian

(24) |

Bounds on are then derived from measurements of the mass splitting and of the involved CP-violating phase , which respectively provide and [36]. In our model we find

(25) |

which is well in agreement with the considered limits on the considered parameter space.

### 4.3 Lepton universality

The coupling involving the third lepton generation in Eq. (17) is subject to tight bounds from the decay of to lighter leptons mediated by gauge interactions in the full theory, shown in Fig. 10. We then repeat for the model at hand the analysis in Ref. [37], assuming the data therein. The observables of reference are the two decay rates

(26) |

where and a suitable normalization is used to cancel the dependence on . The Wilson coefficient in Eq. (26) are linked to the corresponding four leptons operator so that, in absence of new physics, the full SM effect is included in the Fermi constant. New physics contributions are then quantified by computing the induced finite shift of the weak coupling due to the extra scalar fields. At the leading order in the new Yukawa couplings, the divergences of the vertex are cancelled by the field strength renormalization of the and fields, leading to the finite contribution

(27) |

and the Wilson coefficient

(28) |

The conventions for the Passarino-Veltman integrals used in Eq. (27) are taken from Ref. [38], as provided by Ref.s [39, 40]. In the loop computation we neglect all external momenta and the mass of the lepton and neutrinos.

Observable | Experiment | SM prediction |
---|---|---|

By comparing the result obtained with Eq. (26) against the average of Table 3, we find that our benchmark points are subject to the exclusion regions presented in Fig. 11, which constrain the mass splitting amongst the members of the two extra doublets. The more than discrepancy of the SM is always worsened by the interactions of non-degenerate doublets, therefore we refer to bounds in order to assess the exclusion regions of the model in Eq. (17).

### 4.4 Electroweak precision tests

We finally compute the contributions of the extra doublets to the parameter and compare it to the results of electroweak precision tests:

(29) |

where [41]

(30) |

As shown in Fig. 12, these measurements bound the mass splitting of the new doublets to GeV, regardless of the mass of the lightest state.

## 5 Conclusions

If confirmed, the present B-physics anomalies would imply the existence of new physics beyond the Standard Model. Current provisional explanations seem to require exotic new physics, such as leptoquarks characterised by very large couplings close to the perturbative limit.

As an alternative, we propose a simple scalar extension of the SM that explains the observed deviations from lepton universality via a loop contribution that mimics the effect of vector degrees of freedom. The simplified scalar model we consider induces deviations only in vector Wilson coefficients, avoiding in this way the phenomenological bounds that affect traditional scalar extensions. We have shown that our proposal is phenomenologically consistent with all relevant collider as well as precision physics bounds. Our solution requires new scalars with couplings of order unity, thereby avoiding new non-perturbative physics close to the electroweak scale.

Considering the resulting loop contributions on top of the Standard Model prediction reduces the tension with the physics measurements below the level. We speculate that, once these loop effects are considered on top of a small tree level contribution, agreement with the anomalous signal could reach well into the range.

We also investigated a possible high-energy completion of the proposed simplified scalar model, identified in a three-Higgs doublet framework. In this case, measurements of lepton universality and of the parameters strongly constrain the mass splitting between the elements of the new scalar doublets. Given that explanation of the measurements require scalars slightly above the electroweak scale, we expect that future collider searches will exhaustively probe the proposed scalar extension.

## Acknowledgement

It is a pleasure to acknowledge A. Crivellin, A. Djouadi and J. Martin Camalich for useful discussions. We also thank S. Di Chiara for participating in the early stage of this project. This work is supported by the Estonian Research Council grants MOBTT5, IUT23-6, PUT799, PUT808, the ERDF Centre of Excellence project TK133 and the European Research Council grant NEO-NAT.

## References

- [1] Roel Aaij et al. Measurement of the ratio of branching fractions . Phys. Rev. Lett., 115(11):111803, 2015, 1506.08614. [Erratum: Phys. Rev. Lett.115,no.15,159901(2015)].
- [2] Roel Aaij et al. Measurement of the ratio of the and branching fractions using three-prong -lepton decays. 2017, 1708.08856.
- [3] J. P. Lees et al. Evidence for an excess of decays. Phys. Rev. Lett., 109:101802, 2012, 1205.5442.
- [4] J. P. Lees et al. Measurement of an Excess of Decays and Implications for Charged Higgs Bosons. Phys. Rev., D88(7):072012, 2013, 1303.0571.
- [5] M. Huschle et al. Measurement of the branching ratio of relative to decays with hadronic tagging at Belle. Phys. Rev., D92(7):072014, 2015, 1507.03233.
- [6] Y. Sato et al. Measurement of the branching ratio of relative to decays with a semileptonic tagging method. Phys. Rev., D94(7):072007, 2016, 1607.07923.
- [7] S. Hirose et al. Measurement of the lepton polarization and in the decay . Phys. Rev. Lett., 118(21):211801, 2017, 1612.00529.
- [8] Heavy Flavour Average Group. Average of R(D) and R(D*) for FPCP 2017.
- [9] Andreas Crivellin, Christoph Greub, and Ahmet Kokulu. Explaining , and in a 2HDM of type III. Phys. Rev., D86:054014, 2012, 1206.2634.
- [10] Syuhei Iguro and Kazuhiro Tobe. in a general two Higgs doublet model. 2017, 1708.06176.
- [11] Jong-Phil Lee. in the 2HDM with an anomalous coupling. Phys. Rev., D96(5):055005, 2017, 1705.02465.
- [12] Rodrigo Alonso, Benjamín Grinstein, and Jorge Martin Camalich. Lifetime of Constrains Explanations for Anomalies in . Phys. Rev. Lett., 118(8):081802, 2017, 1611.06676.
- [13] A. G. Akeroyd and Chuan-Hung Chen. Constraint on the branching ratio of from LEP1 and consequences for anomaly. Phys. Rev., D96(7):075011, 2017, 1708.04072.
- [14] R. Martinez, C. F. Sierra, and German Valencia. Beyond with the general 2HDM-III for . 2018, 1805.04098.
- [15] Xin-Qiang Li, Ya-Dong Yang, and Xin Zhang. Revisiting the one leptoquark solution to the R(D) anomalies and its phenomenological implications. JHEP, 08:054, 2016, 1605.09308.
- [16] Eugenio Megias, Giuliano Panico, Oriol Pujolas, and Mariano Quiros. A Natural origin for the LHCb anomalies. JHEP, 09:118, 2016, 1608.02362.
- [17] Pouya Asadi, Matthew R. Buckley, and David Shih. It’s all right(-handed neutrinos): a new model for the anomaly. 2018, 1804.04135.
- [18] Admir Greljo, Dean J. Robinson, Bibhushan Shakya, and Jure Zupan. from and right-handed neutrinos. 2018, 1804.04642.
- [19] Mohammad Abdullah, Julian Calle, Bhaskar Dutta, Andrés Flórez, and Diego Restrepo. Probing a simplified, model of anomalies using -tags, leptons and missing energy. 2018, 1805.01869.
- [20] Wolfgang Altmannshofer, P. S. Bhupal Dev, and Amarjit Soni. anomaly: A possible hint for natural supersymmetry with -parity violation. 2017, 1704.06659.
- [21] Walter D. Goldberger. Semileptonic B decays as a probe of new physics. 1999, hep-ph/9902311.
- [22] Vincenzo Cirigliano, Martin Gonzalez-Alonso, and Michael L. Graesser. Non-standard Charged Current Interactions: beta decays versus the LHC. JHEP, 02:046, 2013, 1210.4553.
- [23] Rupak Dutta. Exploring , and anomalies. 2017, 1710.00351.
- [24] Debjyoti Bardhan, Pritibhajan Byakti, and Diptimoy Ghosh. A closer look at the R and R anomalies. JHEP, 01:125, 2017, 1610.03038.
- [25] Ryoutaro Watanabe. New Physics effect on in relation to the anomaly. 2017, 1709.08644.
- [26] Marat Freytsis, Zoltan Ligeti, and Joshua T. Ruderman. Flavor models for . Phys. Rev., D92(5):054018, 2015, 1506.08896.
- [27] A. G. Akeroyd et al. Prospects for charged Higgs searches at the LHC. Eur. Phys. J., C77(5):276, 2017, 1607.01320.
- [28] Adam Alloul, Neil D. Christensen, Céline Degrande, Claude Duhr, and Benjamin Fuks. FeynRules 2.0 - A complete toolbox for tree-level phenomenology. Comput. Phys. Commun., 185:2250–2300, 2014, 1310.1921.
- [29] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer, H. S. Shao, T. Stelzer, P. Torrielli, and M. Zaro. The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations. JHEP, 07:079, 2014, 1405.0301.
- [30] Daniel Dercks, Nishita Desai, Jong Soo Kim, Krzysztof Rolbiecki, Jamie Tattersall, and Torsten Weber. CheckMATE 2: From the model to the limit. Comput. Phys. Commun., 221:383–418, 2017, 1611.09856.
- [31] V. Khachatryan et al. Search for a charged Higgs boson in pp collisions at TeV. Journal of High Energy Physics, 2015(11):18, Nov 2015.
- [32] A. Cordero, J. Hernandez-Sanchez, V. Keus, S. F. King, S. Moretti, D. Rojas, and D. Sokolowska. Dark Matter Signals at the LHC from a 3HDM. JHEP, 05:030, 2018, 1712.09598.
- [33] Christoph Bobeth, Mikolaj Misiak, and Joerg Urban. Matching conditions for and in extensions of the standard model. Nucl. Phys., B567:153–185, 2000, hep-ph/9904413.
- [34] Quan-Yi Hu, Xin-Qiang Li, and Ya-Dong Yang. decay in the Aligned Two-Higgs-Doublet Model. Eur. Phys. J., C77(3):190, 2017, 1612.08867.
- [35] M. Misiak et al. Updated NNLO QCD predictions for the weak radiative B-meson decays. Phys. Rev. Lett., 114(22):221801, 2015, 1503.01789.
- [36] Kfir Blum, Yuval Grossman, Yosef Nir, and Gilad Perez. Combining K0 - anti-K0 mixing and D0 - anti-D0 mixing to constrain the flavor structure of new physics. Phys. Rev. Lett., 102:211802, 2009, 0903.2118.
- [37] Sofiane M. Boucenna, Alejandro Celis, Javier Fuentes-Martin, Avelino Vicente, and Javier Virto. Phenomenology of an model with lepton-flavour non-universality. JHEP, 12:059, 2016, 1608.01349.
- [38] Ansgar Denner. Techniques for calculation of electroweak radiative corrections at the one loop level and results for W physics at LEP-200. Fortsch. Phys., 41:307–420, 1993, 0709.1075.
- [39] Thomas Hahn. Generating Feynman diagrams and amplitudes with FeynArts 3. Comput. Phys. Commun., 140:418–431, 2001, hep-ph/0012260.
- [40] T. Hahn and M. Perez-Victoria. Automatized one loop calculations in four-dimensions and D-dimensions. Comput. Phys. Commun., 118:153–165, 1999, hep-ph/9807565.
- [41] W. Grimus, L. Lavoura, O. M. Ogreid, and P. Osland. A Precision constraint on multi-Higgs-doublet models. J. Phys., G35:075001, 2008, 0711.4022.