New gauge boson W′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W^\prime $$\end{document} and radiative leptonic decays of charged B and D mesons

Considering the constraints on model parameters from the latest experimental data, we study the contributions of new gauge boson W′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W^\prime $$\end{document} predicted by the general 221 models to the radiative leptonic decays of charged B and D mesons. It is found that the largest contributions of W′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W^\prime $$\end{document} in the non-universal 221 model to the branching fractions occur in the decay modes B-→γμ-ν¯μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^{-} \rightarrow \gamma \mu ^{-} \bar{\nu }_\mu $$\end{document} and B-→γτ-ν¯τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^{-} \rightarrow \gamma \tau ^{-} \bar{\nu }_\tau $$\end{document}, which can reach to about 8%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8 \%$$\end{document} of the corresponding SM values in the up-quark mixing scenario. In contrast, the relative corrections of W′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W^\prime $$\end{document} predicted by the un-unified 221 model to the branching fractions are very tiny and identical for all the decay modes.


Introduction
In spite of the tremendous success of the standard model (SM), there are still open questions which imply the existence of new physics beyond it. One of the most common ways of building new physics models is to extend the gauge structure of the SM. An additional U (1) X gauge symmetry is introduced in the simplest extension. One of the next simplest extensions involves an additional SU (2) such as the left-right [1][2][3], leptophobic, hadrophobic, fermiophobic [4][5][6], nonuniversal [7][8][9][10][11] and un-unified models [12,13], which are called the general 221 models in the literature [14]. The most important feature of these models is the existence of new heavy gauge bosons W and Z which are relevant to the low energy charged and neutral current processes, respectively. These new gauge bosons can have rich phenomenology in the current and future collider experiments.
The study of heavy meson decays is one of the most interesting and challenging fields in high energy physics, which involves both strong and weak interactions. In recent years, both experimental and theoretical studies have been improved greatly. However, the limitation in understanding a e-mail: zuoyabing@lnnu.edu.cn b e-mail: cxyue@lnnu.edu.cn and controlling the strong interactions in non-perturbative region is so far still a problem. Compared with hadronic decays, the leptonic decays are simpler due to the fact that strong interaction only occurs in the initial state. The pure leptonic decays of heavy mesons H → lν l (H = B, B s , D, D s and l = e, μ, τ ) can be used to determine the decay constants f H if the CKM matrix elements |V Qq | (Q = b, c and q = u, d, s) are known. Conversely, if we know the values of decay constants from other methods, e.g. the lattice calculation, these decays can also be used to extract the CKM matrix elements. In order to cancel the infrared divergences appearing in the loop calculations, the corresponding radiative leptonic decays H → γ lν l should be considered. More importantly, the radiative leptonic decays are of interest in themselves, which may be separated properly and compared with measurements directly as long as the theoretical softness of the photon corresponds to the experimental resolutions. As well known, the measurement of pure leptonic decays is very difficult because of the helicity suppression and only one detected final state particle. However, having an extra real photon emitted in the final state, the radiative leptonic decays escape from the helicity suppression and therefore can be as large as, or even larger than the pure leptonic mode. Furthermore, with the extra photon to identify, the reconstruction of these decays is easier to do. For these reasons, the radiative leptonic decays have so far received a great deal of attention and been studied with various methods in the literature including the constituent quark model [15][16][17][18][19][20][21][22], light front quark model [23,24], QCD factorization [25][26][27][28][29][30], pQCD [31], light cone sum rules [32,33], dispersion approach [34][35][36][37] and so on, as a means of probing aspects of the strong and weak interactions in a heavy meson system. Concretely, the radiative leptonic decays have been used as an alternative and more viable way for the determination of f H and |V Qq |. In addition, the branching fractions of these decays are sensitive to the inverse moment λ H of the leading twist light cone distribution amplitude of heavy mesons which also enters the QCD factorization formulae for hadronic decays. A measurement of Br(H → γ lν l ) thus provides a determination of λ H free of hadronic final state uncertainties. Moreover, the radiative leptonic decays play a peculiar role in understanding the dynamics of annihilation processes occurring in heavy meson decays which are relevant for many important quantities. Apart from the above, these decays also provide a good opportunity to study the factorization approach, where many ingredients required for the hadronic decays already appear but in a simplified setting.
Besides those aspects in the SM, these decays provide a way of probing possible new physics effects related to the charged current interactions. In Ref. [38], the effects of charged Higgs bosons in the two Higgs doublet model type II (2HDMII) on the radiative leptonic decays of charged B and D mesons were studied based on the factorization procedure. Using the constraints on the free parameters of the 2HDMII given in previous works, it was found that the corrections of charged Higgs boson to the decay B → γ τ ν τ could be as large as 13%. To our knowledge, there is little information available in the literature about the contributions of new charged gauge boson W to the radiative leptonic decays of heavy mesons. In the present paper, we intend to study the corrections of W predicted by the general 221 models to the radiative leptonic decays of charged B and D mesons. The extensions to the contributions of W in other models to these decays are straightforward with simply changing the mass and relevant couplings. The remaining part of this paper is organized as follows. In Sect. 2, we describe the general 221 models briefly and give the corresponding couplings related to these decays. In Sect. 3, we give the formulae for the decay rates including the contributions of W and define the relative corrections to the branching fractions. The numerical results and discussions are given in Sect. 4. Section 5 is our summary.

General 221 models and couplings relevant for the radiative leptonic decays of charged B and D mesons
As mentioned above, the general 221 models are extensions of the SM by introducing an additional SU (2) gauge symmetry. Here we only briefly revisit the basics and some aspects of particular importance to the present work and refer the reader to the original papers for a more thorough description [1][2][3][4][5][6][7][8][9][10][11][12][13][14]. The most important feature of these models is the existence of new heavy gauge bosons W and Z which are relevant to the low energy charged and neutral current processes, respectively. According to the patterns of gauge symmetry breaking, the 221 models can be divided into two classes. The first stage gauge symmetry breaking in the models of class I and II are SU (2) 2 × U (1) X → U (1) Y and SU (2) 1 × SU (2) 2 → SU (2) L , respectively. The charge assignments of the SM fermions and the Higgs representation used to achieve the gauge symmetry breaking can be found in Ref. [14]. After gauge symmetry breaking, the Lagrangian involving the gauge boson masses and fermionic gauge interactions can be written as [14] with In Eqs. (1)-(5), the L 1 and L 2 correspond to the mass terms of the ordinary (Z , W ) and heavy (Z , W ) gauge bosons, respectively. The L 3 denotes the neutral current fermionic gauge interactions and L 4 represents the charged current gauge interactions.M 2 Z ,W are the mass squared of heavy gauge bosons obtained after the first stage symmetry breaking. The parameters ΔM 2 Z ,W denote the additional contributions to the mass squared of heavy gauge bosons after the second stage symmetry breaking, while δM 2 Z ,W are related to the mixing effects between the ordinary and heavy gauge bosons. Comparing withM 2 Z ,W , these four parameters are all at the next leading order of v 2 /u 2 expansion where v and u are the scales of the first and second stage symmetry breaking, respectively. Therefore, at the leading order of v 2 /u 2 expansion, the heavy gauge bosons W and Z in the 221 model of class II are degenerate in mass, The charged current K + μ for the LR, LP, HP, FP, NU and UU models. For the NU and UU models, the mixing angle of the extended gauge group θ E is defined as tan θ E = g 1 /g 2

Quark part
Lepton part (ν e γ μ P R e +ν μ γ μ P R μ +ν τ γ μ P R τ ) The J ± μ , J 0 μ and J μ mainly couple to the ordinary gauge bosons (Z , W , A) and are corresponding to the weak and electromagnetic currents in the SM. The K ± μ and K 0 μ mainly couple to the heavy gauge bosons, the explicit expressions of which depend on the specific models and can be found in Ref. [14].
Note that the subject of this work is to investigate the contributions of heavy gauge boson W to the radiative leptonic decays of charged B and D mesons. We list the charged current K + μ for the left-right (LR), leptophobic (LP), hadrophobic(HP), fermiophobic (FP), non-universal (NU) and ununified (UU) models in Table 1. The current K − μ is just the Hermitian conjugate of K + μ . It is easily seen from Eqs. (1)-(5) and Table 1 that the heavy gauge boson W in the LR, LP, HP, FP models does not couple to the ordinary left-handed neutrinos and therefore makes no contributions to these decays at the leading order of v 2 /u 2 expansion. Therefore, we only need to consider the contributions from W in the NU and UU models.
In the fermion mass eigenstate basis, the gauge couplings of W with fermions in the NU and UU models have the following forms at the leading order of v 2 /u 2 expansion [10,11,14], where For the quark sector, the unitary transformation matrices As analyzed in Ref. [11], the flavor changing neutral currents must occur in the interaction of quarks to gauge bosons, which can be realized in three possible ways, i.e. in the downquark sector only, in the up-quark sector only and in both sectors. Here we only consider the first two ways for simplicity, known as the down-quark mixing scenario and up-quark mixing scenario, respectively. In the down-quark mixing scenario, and the additional charged currents mixing matrix for quarks In the up-quark mixing scenario, and correspondingly Fig. 1 The Feynman diagrams of the radiative leptonic decays of charged B and D mesons at tree level in the SM For the lepton sector, the unitary transformation matrix T E can be written as It is easy to show that the currents mixing matrix for leptons Thus, the interaction of leptons to gauge bosons only depends on the third row of unitary matrix T E , i.e. T 31 , T 32 and T 33 .

Decay rates and relative corrections to the branching fractions
In the radiative leptonic decays of charged B and D mesons, the strong interactions are involved only in the initial states. The Feynman diagrams of these decays at tree level in the SM are given in Fig. 1. The tree level diagrams corresponding to the contributions of W are similar to those in Fig. 1 with simply replacing W with W . Extending the expression for the amplitude in the SM given in Ref. [38] to include the contributions of W by using Eqs. (7)-(10), we have with and P denotes the charged B or D meson and ε is the polarization vector of the photon in final state. p P and p γ are the 4-momentum of P and the photon, respectively. The form The contribution of (d) is suppressed by a factor of 1/m 2 W and thus is neglected here. Among these quantities, For the form factors F V , F A , the fitted numerical results based on the factorization procedure [38] are where E γ is the energy of the photon and Λ QC D denotes the characteristic energy scale of strong interaction.
In the present work, we only consider the contributions of W to the leading order of v 2 /u 2 expansion. Correspondingly, according to Eqs. (18) and (20), (21), the amplitude squared is In the rest frame of P, the initial and final state 4-momenta can be written as where the momentum of the photon has been chosen to be in the − z direction as shown in Fig. 2. E l is the energy of the charged lepton in final state. Then, we have Note that the sum of momenta is conserved in the decay and the particles in the initial and final state are on mass shell, the p lz can be expressed in terms of E γ and E l , i.e.
With the amplitude squared in Eq. (26) and relevant expressions given above, we can calculate the differential decay rates for radiative leptonic decays of charged B and D mesons by using the formula [39] The total decay rates can be obtained by integrating the differential decay rates over E γ and E l , with the integral limits Note that since there are IR divergences in the radiative leptonic decays when the photon is soft or collinear with the emitted lepton, the total decay rates depend on the experimental resolution to the energy of the photon ΔE γ . By using the lifetimes of charged B and D mesons, we can obtain the branching fractions from the total decay rates for all relevant decay modes. To see the corrections of W predicted by the NU and UU models more clearly and reduce the possible uncertainties stemming from the CKM matrix elements and form factors, we define the relative corrections of W to the branching fractions of these decays where Br total and Br SM are the branching fractions of these decays with and without the contributions of W , respectively.

Numerical results and discussions
Before performing the actual calculations, we need to consider the constraints on the model parameters from the latest experimental data. As stated above, the mixing angle of the extended gauge group θ E is defined as tan θ E = g 1 /g 2 .
Assuming the perturbativity of all the gauge couplings, g 2 1,2 /4π < 1, we have 0.03 < sin 2 θ E < 0.96 [40,41]. It is easily seen from Eq. (6) that the mass of heavy gauge boson m W is relevant to sin 2 θ E . Considering the latest constraints from the CKM unitarity, lepton flavor violating processes and direct search for heavy gauge bosons at the LHC [40][41][42], m W can be chosen as m W > 1800 GeV and m W > 3000 GeV in the small and large sin 2 θ E region, respectively. The T 31 , T 32 , T 33 in the currents mixing matrix for leptons T † E ΔT E are constrained by considering the unitarity of T E and lepton flavor violating processes [40]. It is concluded that one of the matrix elements T 31 , T 32 , T 33 is close to 1 and the others are sufficiently small. It is more natural to assume the mixing strength between leptons to be directly related to their masses. So we only consider the cases that T 33 ≈ 1 or T 32 ≈ 1. Actually, we choose T 33 = 1, T 31 = T 32 = 0 or T 32 = 1, T 31 = T 33 = 0 for simplicity in our calculations.
Other input parameters are as follows [39], First of all, we give the branching fractions of these decays in the SM and corresponding experimental upper limits [39] in Table 2. Since the branching fraction of D + → γ τ + ν τ is too small due to the phase space suppression, we only consider the decay modes D + → γ e + ν e and D + → γ μ + ν μ for the D decays. We can see that the branching fractions are all within their experimental upper limits and decrease slightly with the increase of ΔE γ . The theoretical uncertainties can come from the CKM matrix elements |V ub |, |V cd | and the form factors F V , F A , which are about 10-20% via roughly estimation. Considering that the main concern of this work is the relative corrections of W and the uncertainties mentioned above can be canceled to a large extent in this quantity, we only give the central values for the branching fractions here. In the NU model, it is easily seen from Eqs. (18) and (19) that the contributions of W to the radiative leptonic decays of charged B and D mesons depend on the mixing angle of the extended gauge group θ E and the charged currents mixing matrices T † U ΔT D , T † E ΔT E . Considering the relevant arguments given above, we obtain the relative corrections of W to the branching fractions of these decays, which are all independent of the experimental resolution to the energy of the photon ΔE γ due to the fact that the relative corrections is defined via the ratio of branching fractions, c.f. Eq. (33), and therefore the effects of ΔE γ are canceled.
For the decay B − → γ e −ν e , we can see from Eqs. (17) Table 2 The branching fractions of the decays B − → γ l −ν l (l = e, μ, τ ) and D + → γ l + ν l (l = e, μ) in the SM and corresponding experimental upper limits [ are related with T 31 which we have chosen to be zero as stated above. In the down-quark mixing scenario, the relative corrections R as functions of sin 2 θ E are shown in Fig. 3. Apparently, the contributions of W to this decay mode are positive and become more significant in the large sin 2 θ E region. Hence we choose 0.5 < sin 2 θ E < 0.96 and m W > 3000 GeV. The numerical results for R corresponding to different values of m W and sin 2 θ E in the above regions are given in Table 3. We can see that the relative corrections of W to this decay mode in the down-quark mixing scenario are insignificant, 2% for the largest possible case. In the upquark mixing scenario, the corrections of W to the branching fractions of B − → γ e −ν e are independent of sin 2 θ E which can be illuminated by scrutinizing Eqs. (15), (17)- (19), and we choose relatively low values for m W , requiring m W > 1800 GeV. The corresponding numerical results for R are given in Table 4. We can observe that the corrections of W to the branching fractions of this decay mode in the up-quark mixing scenario are negative and tiny, around 0.3% of the values in the SM for the largest possible case.
For the decay B − → γ μ −ν μ , the corrections of W in the NU model to the branching fractions depend on the parameter T 32 which can be chosen to be 0 or 1 as stated above. Choosing T 32 = 0, it is found that the relative corrections of W to the branching fractions of this decay mode are identical to those of B − → γ e −ν e which are shown in Tables 3 and 4 for the down-quark and up-quark mixing scenario, respectively. This is due to the fact that the relative corrections R is defined via the ratio of branching fractions and thus the effects of the mass of charged lepton in the final state are canceled. For the same reason, we shall have the same results for the decay B − → γ τ −ν τ under the condition that T 33 = 0. Under the condition that T 32 = 1, the relative corrections of W to the branching fractions are independent of sin 2 θ E in the downquark mixing scenario and this can be seen by inspecting  Eqs. (13), (17)- (19). Requiring m W > 1800 GeV, we obtain the same numerical results for R as those shown in Table 4. It is worth mentioning that not only the relative corrections to the branching fractions but also the genuine ones are identical for the case T 32 = 0 in the up-quark mixing scenario and T 32 = 1 in the down-quark mixing scenario. In contrast, the relative corrections of W to the branching fractions of B − → γ μ −ν μ decrease with the increase of sin 2 θ E in the up-quark mixing scenario and therefore are more significant in the small sin 2 θ E region as shown in Fig.4. Correspondingly, we choose 0.03 < sin 2 θ E < 0.5 and m W > 1800 GeV. The numerical results for R are given in Table 5. We can observe that the largest possible corrections of W to the decay B − → γ μ −ν μ are about 8% of the corresponding values in the SM. For the decay B − → γ τ −ν τ , the corrections of W predicted by the NU model to the branching fractions of this decay mode depend on the parameter T 33 . As inferred above, under the condition that T 33 = 0, the relative corrections of W to the branching fractions are the same as those of B − → γ e −ν e , c.f. Tables 3 and 4 for the down-quark and up-quark mixing scenario, respectively. By the same token, under the condition that T 33 = 1, the relative corrections of W to this decay mode are identical to those of B − → γ μ −ν μ with T 32 = 1 which are shown in Tables 4 and 5 for the downquark and up-quark mixing scenario, respectively. Similar to the decay B − → γ μ −ν μ , the genuine corrections to the branching fractions of B − → γ τ −ν τ are the same for the case T 33 = 0 in the up-quark mixing scenario and T 33 = 1 in the down-quark mixing scenario.
For the decays D + → γ l + ν l (l = e, μ), it is found that the relative corrections of W predicted by the NU model to the branching fractions are the same in the down-quark and up-quark mixing scenarios, owing to the fact that the relevant matrix elements of T † U ΔT D for the D decays are both zero in these two scenarios, c.f. Eqs. (13) and (15). With T 31 = 0, we have the same numerical results for the relative corrections of W to the branching fractions of D + → γ e + ν e as those shown in Table 3. The effects of the mass difference between D and B mesons in the initial state are canceled due to the fact that the relative corrections R is defined via the ratio of branching fractions. As mentioned above, the effects of the mass of charged lepton in the final state are also canceled in the relative corrections. And indeed the numerical results for the relative corrections to the decay D + → γ μ + ν μ with T 32 = 0 are identical to those shown in Table 3 as well. Under the condition that T 32 = 1, the contributions of W to the branching fractions of this decay mode are independent of sin 2 θ E . Choosing m W > 1800 GeV, we obtain the same numerical results for the relative corrections as those given in Table 4. Now we calculate the contributions of W predicted by the UU model to the branching fractions of B − → γ l −ν l (l = e, μ, τ ) and D + → γ l + ν l (l = e, μ). For this model, it is easily seen from Eqs. (18) and (19) that the corrections of W to the branching fractions of these decays are independent of sin 2 θ E and the charged current mixing matrices In other words, our results in this model only depend on m W . Moreover, the relative corrections of W predicted by this model to the branching fractions are all the same for these decay modes owing to the fact that the effects of the mass difference between B and D in the initial state as well as those of the mass of charged lepton in the final state are all canceled. Choosing m W > 1800 GeV, we have the numerical results for R as shown in the Table 6. We can observe that the corrections of W in the UU model to the branching fractions of B − → γ l −ν l (l = e, μ, τ ) and D + → γ l + ν l (l = e, μ) are negative and very insignificant, about 0.3% of the corresponding values in SM for the largest possible case.

Summary
In the present work, we study the contributions of new heavy gauge boson W predicted by the general 221 models to the radiative leptonic decays of charged B and D mesons. Due to the fact that the heavy gauge boson W in the LR, LP, HP and FP models does not couple to the ordinary left-handed neutrinos and therefore makes no contributions to these decays at the leading order of v 2 /u 2 expansion, we only calculate the contributions of W in the NU and UU models. To see the corrections of W more clearly and reduce the uncertainties stemming from the CKM matrix elements and form factors, we define the relative corrections R via the ratio of branching fractions. Considering the constraints on the model parameters from the latest experimental data, we give the relative corrections of W to the branching fractions of these decays. Due to the fact that the effects on R coming from the mass difference between B and D mesons in the initial state as well as those of the mass of charged lepton in the final state are all canceled, the numerical results for many decay modes are identical as noted above. Specially, the relative corrections of W predicted by the NU model to the branching fractions of B − → γ l −ν l (l = e, μ, τ ) are all the same under the condition that T 31 = 0, T 32 = 0 and T 33 = 0, respectively. For the decays D + → γ l + ν l (l = e, μ), the contributions of W in the NU model are equal in the down-quark and up-quark mixing scenarios. The largest corrections of W in the NU model to the branching fractions occur in the decay modes B − → γ μ −ν μ and B − → γ τ −ν τ , which can reach to about 8% of the values in the SM in the up-quark mixing scenario under the condition that T 32 = 1 and T 33 = 1, respectively. For the time being, these effects are still smaller compared with the corresponding theoretical uncertainties of the SM branching fractions and thus difficult to discern in experiments. It is noteworthy that the near future BELLE II experiment at SuperB will be able to improve the accuracy of measurement on the CKM matrix element |V ub | significantly [43,44]. In addition, it is possible to achieve a better control on the uncertainties of the form factors in the future. Under these circumstances, the theoretical uncertainties of branching fractions in the SM could be reduced to below 5% and our results may be tested by future more precise experiments such as BELLE II, LHCb upgraded. In contrast, the contributions of W predicted by the UU model only depend on m W and the relative corrections to the branching fractions are the same for all the decay modes which are negative and very tiny, about 0.3% for the largest possible case.
ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .