SciELO - Scientific Electronic Library Online

vol.38 issue4A multiparameter chaos control method applied to maps author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand




Related links


Brazilian Journal of Physics

Print version ISSN 0103-9733On-line version ISSN 1678-4448

Braz. J. Phys. vol.38 no.4 São Paulo Dec. 2008 

Leptonic decays of the B charged meson and B Xsγ in the Two Higgs Doublet Model type III



J. P. IdárragaI, II; R. MartinezII; J-Alexis RodríguezII; N. Poveda T.III

IDepartament de Physique, Université de Montréal, Montréal, Canada
IIDepartamento de Física, Universidad Nacional de Colombia, Bogota, Colombia
IIIEscuela de Física, Universidad Pedagógica y Tecnológica de Colombia, Grupo de Física Teórica y Computacional, Tunja, Colombia




We consider the Two Higgs Doublet Model (2HDM) of type III which leads to Flavour Changing Neutral Currents (FCNC) at tree level. In the framework of this model we calculate the NLO contribution for b sγ and the branchings for the meson decays B+ l+ν. We examine the limits on the new parameters λbb and M. We take into account the relationship between λtt and λbb coming from the validness of perturbation theory.

Keywords: B Meson decay; 2HDM type III; Charged Higgs mass



The Standard Model (SM) of particle physics based on the gauge group SU(3)c × SU(2)L × U(1)Y makes fit the symmetry breaking by including a fundamental weak doublet of scalar Higgs bosons with a scalar potential V() = λ(† - v2)2. However, the SM does not explain the dynamics responsible for the generation of masses. Between the models beyond the SM, many of them include more than one scalar Higgs doublet; for instance, the case of the minimal supersymmetric standard model (MSSM)). We consider one of these extensions which has a richer scalar sector, called generically as the Two Higgs Doublet Model (2HDM). There are several kinds of such 2HDM models. In the model called type I, one Higgs doublet provides masses to the up and down quarks, simultaneously. In the model type II, one Higgs doublet gives masses to the up type quarks and the other one to the down type quarks. These two models include a discrete symmetry to avoid flavour changing neutral currents (FCNC) at tree level [1]. However, the addition of these discrete symmetries is not compulsory and in this case both doublets are contributing to generate the masses for up-type and down-type quarks. In the literature, such a model is known as the 2HDM type III [2]. It has been used to search for physics beyond the SM and specifically for FCNC at tree level [3-5]. In general, both doublets can acquire a vacuum expectation value (VEV), but one of them can be absorbed redefining the Higgs boson fields properly. Nevertheless, other studies on 2HDM-III using different basis have been done and there is a case where both doublets get VEVs that allows to study the models type I and II in a specific limit [5, 6].

In the 2HDM models, the two complex Higgs doublets correspond to eigth scalar states. Spontanoues Symmetry breaking procedure leads to five Higgs fields: two neutral CP-even scalars h0 and H0, a neutral CP-odd scalar A0, and two charged scalars H±. While the neutral Higgs bosons may be difficult to distinguish from the one of the SM, the charged Higgs bosons would have a distinctive signal for physics beyond the SM. Therefore the direct or indirect evidence of a charged Higgs boson would play an important role in the discovery of an extended Higgs sector. Direct searches have carried out by LEP collaborations and they reported a combined lower limit on of 78.6 GeV [7] assuming H+ → τ+ντ(c). At the Tevatron, the direct searches for charged Higgs boson are based on where at least one top quark decay can be used in order to look for the channel t H+b. The CDF collaboration has reported a direct search for charged Higgs boson, setting an upper limit on B(t H+b) around 0.36 at 95 % C.L. for masses in the range of 60-160 GeV [8]. On the other hand, indirect and direct searches have been done by D0 looking for a decrease in the W+W&– signal expected from the SM and also the direct search for the decay mode H± → τ±ν. We should note that all these bounds have been gotten in the framework of the 2HDM type II. And, in the framework of the 2HDM type II and MSSM a full one loop calculation of Γ(t bH+) including all sources for large Yukawa couplings were presented in references [10, 11]. Other indirect bounds on the charged Higgs boson mass come from processes where the charged Higgs boson is a virtual particle which is the case of the process b sγ. However, the indirect limits which have been obtained from the branching ratio B Xsγ are strongly model dependent [12]. Finally, the search for the charged Higgs boson mass above the top quark mass will continue at LHC. The main production mechanisms would be the processes gg tbH+ and gb tH+ which have been studied using simulations of the LHC detectors [14].

The charged Higgs boson can also be revealed through contributions to low energy processes such as B0 &– , D0 &– and K0 &– where bounds on the charged Higgs sector have been found [2]. Moreover, there are other options through leptonic decays of the charged B mesons. They occur via the annihilation process B± W*(H*) l±νl. Then, it is possible to use the upper limits on these branching ratios obtained at CLEO [13], BELLE [15] and BABAR [16] in order to get bounds on the charged Higgs boson mass. Recent experimental results on B(Bu → τν) were reported by BELLE [17] and it is the first evidence of this kind of decays. The decays B± l±νl are sensitive at tree level to charged Higgs bosons and can be enhanced up to the current experimental limits by multi-Higgs models[18, 19]. On the other hand, the rare decay B Xsγ is sensitive to charged higgs bosons at one loop level through electromagnetic and chromomagnetic penguin diagrams, and therefore the decay B Xsγ can put strong constraints on the parameters of any charged Higgs sector because its high precision measurement done by CLEO [20].

In the present work, we study the processes B Xsγ and B+ l+ν in the framework of the 2HDM type III. And we concentrate on the charged Higgs boson sector of this model, with the relevant parameters being its mass and the Yukawa coupling.

The 2HDM type III is an extension of the SM which adds a new Higgs doublet and three new Yukawa couplings in the quark and leptonic sectors. The mass terms for the up-type or down-type sector depend on two Yukawa coupling matrices. The rotation of quarks and leptons to mass eigenstates allow us to diagonalize one of the matrices but not both simultaneously, then one of the Yukawa coupling matrix remains non-diagonal, generating FCNC at tree level. The most general structure for the Higgs-fermion Yukawa couplings in the so called 2HDM type-III is as follow[2]:

where Φ1,2 are the Higgs doublets, is the weak isospin quark doublet, and are weak isospin quark singlets, whereas and are non-diagonal 3 × 3 and i, j are family indices. The superscript 0 indicates that the fields are not mass eigenstates yet. In the so-called model type I, the discrete symmetry forbids the terms proportional to , meanwhile in the model type II the same symmetry forbids terms proportional to . On the other hand, we are considering the 2HDM-III in a basis where only one Higgs doublet acquire VEV and then it does not have the usual parameter tanβ = ν2/ν1 of the 2HDM type II. In this way we have the usual 2HDM type III [4], where the Lagrangian of the charged sector is given by

where K is the Cabbibo-Kobayashi-Maskawa (CKM) matrix and ξU,D the flavour changing matrices. In the framework of the 2HDM type III is useful to use the parameterization proposed by Cheng and Sher for the Yukawa couplings given by ξii = λii gmi/(2mW)[4].

The leptonic decays of the B± mesons are possible via annihilation processes into W± gauge bosons or H± scalar bosons. The first one is the usual SM contribution and the second one in our case is own to the 2HDM type III. We should mention that the branching fractions for and in the framework of the SM are helicity suppressed by factors of ~ 10-8 and ~ 10-3, respectively. But physics beyond the SM can enhance these branching fractions through the introduction of a charged Higgs boson, as we will notice. The decay width can be written as

where in the framework of the 2HDM-III, we have the factors

And R can be written as

where we have used the Cheng and Sher parametrization [4]. When R = 1, i.e., the new physics predicts the same value for the width Γ(B± l±νl)III that the SM, and there is an allow region in the plane mH &– λbb with the same SM values. In this form the decay width depends only on the free parameters λll, and λbb.

About the experimental data for the B meson decays B&– l&–, they are experimentally challenging because there are at least two undectetable neutrinos in the final state. These kind of decays has been searched at BELLE, BABAR and CLEO-b. Bounds on the braching fraction B(B → µν) have been reported. The stringent bounds come from BABAR measurements and they are B(B → µνµ) < 6.8 ×10-6 [16] and the SM prediction is B(B → µνµ) = 3.9 ×10-7. For the decay Bu → τν, the first evidence has been reported by BELLE [17], they report an experimental result of B(B → τντ) = ×10-4. In addition, the values predicted by the SM is B(B → τντ) = 1.59 ×10-4, which is consistent with the experiment value. This new measurement could guide to a deeper undertanding of flavour and electroweak dynamics, and it could provide evidence of a non-standard Higgs sector. As we already mentioned in the B meson decays is possible to reduce the number of parameters to λbb and the charged Higgs boson mass mH where we have used the flavour changing couplings for the leptonic sector from the literature[5]. These couplings are bounded by -0.12 < ξ22 < 0.12 and -1.8 ×10-2 < ξ33 < 2.2 ×10-2. These bounds can also be written as |λµµ| < 139.82 and |λττ| < 1.52. For our numerical analysis we always take the upper bounds λµµ= 139.82 and λττ= 1.52.

The Fig. 1 shows the allowed values (fullfilled) for mH vs λbb according to the experimental result from BELLE [17] for the B → µνµ (left) and B → τντ (right) decays where the allow regions correspond to the fullfilled regions. For the leptonic parameter we take the values λµµ= 139.82 and λττ= 1.52, respectively. For B → µνµ (left), the lower solid line corresponds to the experimental limit and the upper solid line corresponds to the solutions given by R = 1 which predicts for the new physics the same SM width. For B → τντ (right), the upper curve corresponds to experimental limit and the lower solid line corresponds to R = 1 solution which is excluded.



On the other hand, for the radiative decay B Xsγ we follow references [23, 24]. The B Xsγ process as any FCNC process does not arise at the tree level in the SM. In the framework of the SM it is generated by the one-lopp W-exchange diagrams but these contributions are small enough to be comparable to nonstandard contributions, in our case the exchange of a charged scalar Higgs boson. The branching ratio of the inclusive radiative decay B Xsγ is

at the leading order level, where (µb) is the effective coefficient at the scale µb,

f(z) = 1 - 8z2 + 8z6 - z8 - 24z4 log z is the phase space factor in the semileptonic b-decay parameterized in terms of z = and αem is the fine-structure constant. The coefficients (µb) have an important property and it is that they are quite similar in many interesting extensions of the SM, such as 2HDM or the MSSM [21-23] and therefore it is possible to parametrize the new contributions using new fuctions (µW) with i = 7,8 and j = YY, XY. These functions depend on the unknown parameter and also on the size and sign fo the couplings X and Y that in the case of the model III under study they are X = -λbb, and Y = λtt. To get these couplings we assume that the flavour changing parameters for the light quarks are negligible and λbb > 1, λtt < 1 which is the case disccussed by Atwood, Reina and Soni as their third case [3]. Then the LO Wilson coefficients at the matching energy scale mW are [23, 24],


where , and these leading order functions have no explicit µW dependence.

Now, at the next leading order level that is neccesary in order to use the experimental data, the branching ratio is

where BSL is the measured semileptonic branching ratio of B mesons, and κ(z) is the QCD correction for the semileptonic B decay. The term corresponds to the subprocesses b sγ which involves the NLO Wilson coefficient (µb), the virtual correction functions ri and the elements of the anomalous dimension matrix which govern the evolution of the Wilson coefficents from the matching scale µW to lower scale µb. The term A in equation (12) is the correction coming from the bremsstrahlung process b sγg and in the Δ have been included the nonperturbative corrections obtained with the method of the heavy-quark effective theory relating the actual hadronic process to the quark decay rate. The whole set of functions already mentioned have been given in references [23, 24]. With the set of above equations we can estimate the ratio B(B Xsγ) and use the experimental world averageB(B Xsγ)exp = (3.52 ± 0.30) ×10-4 [27].

From the perturbation theory considerations in the eq. (2) and using the Cheng Sher parametrization , we have already gotten the inequality [25]

which will be taken into account in the numerical evaluations for the B(B Xsγ) in order to reduce the number of variables. It defines an ellipse with |λbb| < 100 and |λtt| < bounds. In this case we consider the inequality from perturbation theory validness in order to reduce the space of parameters. This link between the parameter λtt and λbb allows to get the plane λbb versus mH using the experimental measurement for the branching ratio B(Bs Xsγ).

In Fig. 2, we present the allowed regions (fullfilled) for the B(B Xsγ) in the plane MH versus λbb for different values of λtt. And in Fig. 3, we present the allowed regions (fullfilled) for the B(B Xsγ) in the plane λbb-λtt for different values of the charged Higgs boson mass.





Finally, in Fig. 4 we show the case of the induced decay B Xsγ decay. The fullfilled regions are the allowed regions, it means these are the regions satisfying the experimental value and the perturbation theory constraint eq. (12). We notice that these regions in Fig. 4 correspond to a different choice of λtt as it was presented in Fig. 2.



In conclusion, We have studied in the framework of the 2HDM type III, the allowed region for the parameters λbb and mH using the processes B → τντ, B → µνµ and B Xsγ. We have used the condition on the parameter space coming from the fact that the Yukawa couplings should be perturbative, eq. (13), in to order to reduce the number of free parameters. Finally, we have compared the plots looking for the stringest regions in the plane λbb-mH and we have noticed that the B Xsγ decay is the most restrictive process constraining the parameters of the charged Higgs sector in the 2HDM-III. But however there are small regions for small values of λbb and light mH that leptonic decays can exclude. We also have found that in case of the leptonic decays B lν, there are values of the parameters λbb and mH given a 2HDM predition which cannot be distinguishable from the SM prediction. It is because the factor (1 - dbMB/)2 in eq. (3) could get the value equal one for some values of λbb and mH and then reach out the SM prediction. For the process B → µνµ, these values are in the allowed experimental region of the plane λbb-mH as it was showed in Fig. 1. But for the process B → τντ these ones are excluded.

This work has been supported by the Fundación Banco de la República.


[1] S. Glashow and S. Weinberg, Phys. Rev. D 15, 1958 (1977).         [ Links ]

[2] W. S. Hou, Phys. Lett. B 296, 179 (1992);         [ Links ] D. Cahng, W. S. Hou, and W. Y. Keung, Phys. Rev. D 48, 217 (1993);         [ Links ] S. Nie and M. Sher, Phys. Rev. D 58, 097701 (1998);         [ Links ] M. Sher and Y. Yuan, Phys. Rev. D 44, 1461 (1991);         [ Links ] D. Atwood, L. Reina, and A. Soni, Phys. Rev. D 55, 3156 (1997);         [ Links ] Phys. Rev. Lett. 75, 3800 (1995).         [ Links ]

[3] D. Atwood, L. Reina, and A. Soni, Phys. Rev. D 53, 1199 (1996);         [ Links ] Phys. Rev. D 54, 3296 (1996);         [ Links ] Phys. Rev. Lett. 75, 3800 (1993);         [ Links ] D. Atwood, L. Reina, and A. Soni, Phys. Rev. D 55, 3156 (1997);         [ Links ] G. Cvetic, S. S. Hwang, and C. S. Kim, Phys. Rev. D 58, 116003 (1998).         [ Links ]

[4] Marc Sher and Yao Yuan, Phys. Rev. D 44, 1461 (1991);         [ Links ] T. P. Cheng and M. Sher, Phys. Rev. D 35, 3490 (1987).         [ Links ]

[5] Rodolfo A. Diaz, R. Martinez, and J.-Alexis Rodriguez, Phys. Rev. D 64, 033004 (2001);         [ Links ] Phys. Rev. D 63, 095500 (2001).         [ Links ]

[6] S. Davidson and H. E. Haber, Phys. Rev. D 72, 035004 (2005) [arXiv:hep-ph/0504050]         [ Links ].

[7] LEP Collaborations, arXiv:hep-ex/0107031.         [ Links ]

[8] CDF Collaboration, F. Abe, et. al, Phys. Rev. Lett. 79, 357 (1997);         [ Links ] CDF Collaboration, T. Affolder, et. al, Phys. Rev. D 62, 012004 (2000).         [ Links ]

[9] D0 Collaboration, B. Abbot, et. al, Phys. Rev. Lett. 82, 4975 (1999);         [ Links ] V. Abazov, et. al, Phys. Rev. Lett. 88, 151803 (2002).         [ Links ]

[10] A. Mendez and A. Pomarol, Phys. Lett. B 360, 47 (1995);         [ Links ] C. Li and R. J. Oakes, Phys. Rev. D 43, 855 (1991);         [ Links ] A. Djouadi and P. Gambino, Phys. Rev. D 51, 218 (1995).         [ Links ]

[11] M. Carena, D. Garcia, U. Nierste, and C. Wagner, arXiv:hep-ph/9912516;         [ Links ] M. Carena, J. Conway, H. Haber, and J. Hobbs, arXiv:hep-ph/0010338.         [ Links ]

[12] W. S. Hou, Phys. Rev. D 48, 2342 (1993).         [ Links ]

[13] S. Chen et al. [CLEO Collaboration], Phys. Rev. Lett. 87, 251807 (2001) [arXiv:hep-ex/0108032]         [ Links ].

[14] K. A. Assamagan, Y. Coadou, and A. Deandrea, Eur. Phys. J. directC 4, 9 (2002) [arXiv:hep-ph/0203121]         [ Links ].K. A. Assamagan and N. Gollub, Eur. Phys. J. C 39S2, 25 (2005) [arXiv:hep-ph/0406013]         [ Links ].

[15] K. Abe et al. [Belle Collaboration], KEK-PREPRINT-2001-88 Prepared for 20th International Symposium on Lepton and Photon Interactions at High Energies (LP 01), Rome, Italy, 23-28 Jul 2001         [ Links ]

[16] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 95, 041804 (2005) [arXiv:hep-ex/0407038]         [ Links ].

[17] K. Ikado,[BELLE Collaboration],hep-ex/0604018.         [ Links ]

[18] A. G. Akeroyd and S. Recksiegel, J. Phys. G 29, 2311 (2003) [arXiv:hep-ph/0306037]         [ Links ].

[19] G. Isidori and P. Paradisi, hep-ph/0605012.         [ Links ]

[20] M. Artuso et al. [Cleo Collaboration], Phys. Rev. Lett. 75, 785 (1995).         [ Links ]

[21] B. Grinstein and M. B. Wise, Phys. Lett. B 201, 274 (1988);         [ Links ] P. L. Cho and B. Grinstein, Nucl. Phys. B 365, 279 (1991) [Erratum-ibid. B 427, 697 (1994)]         [ Links ]. B. Grinstein, R. P. Springer, and M. B. Wise, Phys. Lett. B 202, 138 (1988).         [ Links ]

[22] S. Bertolini, F. Borzumati, A. Masiero, and G. Ridolfi, Nucl. Phys. B 353, 591 (1991);         [ Links ] S. Bertolini and F. Vissani, Z. Phys. C 67, 513 (1995) [arXiv:hep-ph/9403397]         [ Links ].

[23] F. M. Borzumati and C. Greub, Phys. Rev. D 58, 074004 (1998) [arXiv:hep-ph/9802391]         [ Links ].

[24] Z. J. Xiao and L. Guo, Phys. Rev. D 69, 014002 (2004) [arXiv:hep-ph/0309103]         [ Links ].

[25] R. Martinez, J. A. Rodriguez, and M. Rozo, Phys. Rev. D 68, 035001 (2003) [arXiv:hep-ph/0212236]         [ Links ].

[26] R. A. Diaz, R. Martinez, and C. E. Sandoval, [arXiv:hep-ph/0311201] (2003).         [ Links ] R. A. Diaz, R. Martinez, and C. E. Sandoval, [arXiv:hep-ph/0406265] (2004).         [ Links ] R. A. Diaz, R. Martinez, and J-A. Rodriguez, [arXiv:hep-ph/0103050] (2001);         [ Links ] R. A. Diaz, R. Martinez, and J-A. Rodriguez, Phys. Rev. D 63, 095007 (2001), [arXiv:hep-ph/0010149]         [ Links ].

[27] E. Barberio, et. al [The Heavy Flavor Averaging Group], hep-ex/0603003.         [ Links ]



(Received on 15 February, 2008)

Creative Commons License All the contents of this journal, except where otherwise noted, is licensed under a Creative Commons Attribution License