Muon Anomalous Magnetic Moment in SUSY B-L Model with Inverse Seesaw

Motivated by the tension between the Higgs mass and muon g-2 in minimal supersymmetric standard model (MSSM), we analyze the muon g-2 in supersymmertic B-L extension of the standard model (BLSSM) with inverse seesaw mechanism. In this model, the Higgs mass receives extra important radiative corrections proportional to large neutrino Yukawa coupling. We point out that muon g-2 also gets significant contribution, due to the constructive interferences of light neutralino effects. The light neutralinos are typically the MSSM Bino like and the supersymmetric partner of U(1)_B-L gauge boson (B'-ino). We show that with universal soft supersymmetry breaking terms, the muon g-2 resides within 2 sigma of the measured value, namely ~17 x 10^-10, with Higgs mass equal to 125 GeV .

Motivated by the tension between the Higgs mass and muon g − 2 in minimal supersymmetric standard model (MSSM), we analyze the muon g − 2 in supersymmertic B − L extension of the standard model (BLSSM) with inverse seesaw mechanism. In this model, the Higgs mass receives extra important radiative corrections proportional to large neutrino Yukawa coupling. We point out that muon g − 2 also gets significant contribution, due to the constructive interferences of light neutralino effects. The light neutralinos are typically the MSSM Bino like and the supersymmetric partner of U (1)B−L gauge boson (B ino). We show that with universal soft supersymmetry breaking terms, the muon g − 2 resides within 2σ of the measured value, namely ∼ 17 × 10 −10 , with Higgs mass equal to 125 GeV .
This discrepancy has survived after performing highly accurate theoretical calculations [1] within the SM framework and experimental analyses [2]; and hence, it can be resolved or ameliorated by contributions from new physics beyond the SM (BSM). If supersymmetry (SUSY), as one of the forefront candidates for the BSM physics, is a solution to the muon g − 2, the SUSY particles, namely, smuon and weak gaugino (Bino or Wino) masses should be around a few hundred GeV, in order to utilize the supersymmetric contributions [3]. However, the observation of the Higgs boson of mass about 125 GeV requires rather heavy sparticle spectrum within the MSSM framework, and it results in a strong tension in simultaneous resolution for both the 125 GeV Higgs boson and the muon g−2 problem since SUSY contributions to muon g − 2 is suppressed by the heavy spectrum. Non-universality in gaugino and/or scalar masses may remove this tension [4], nevertheless in this case SUSY models will have plenty of free parameters and will lose their productivity.
In this article we show that this tension can be alleviated in the U (1) B−L extended Supersymmetric Standard Model (BLSSM), which is one of the interesting nonminimal realizations of supersymmetry. The BLSSM is also well motivated by the established existence of nonzero neutrino masses [5]. It turns out that, in this class of model, the scale of B − L symmetry breaking can be related to the scale of SUSY breaking [6], therefore, a TeV scale type I or inverse seesaw mechanism can be naturally implemented [7]. The one-loop radiative corrections to the SM-like Higgs boson mass, due to the right-handed (s)neutrinos in BLSSM, with inverse seesaw mechanism, provide new contribution to the Higgs boson mass in addition to the stop sector of MSSM [8]. Thus the lower bound imposed by Higgs mass on the uni-versal gaugino soft masses m 1/2 is reduced, which makes possible to find solutions with the light weak gauginos. Moreover, in BLSSM with inverse seesaw the g−2 may receive new contributions, in addition to the usual MSSM ones, due to the extension of the neutralino sector by SM singlet (B − L) Higgsino and B -ino and also due to the possibility that one of the right-handed sneutrinos is light (due to large mixing between right-handed sneutrinos and right-handed anti-sneutrinos). We emphasize that within BLSSM with inverse seesaw the g − 2 resides within 2σ of the measured value. We show that in this model, even if we are restricted to the universal boundary conditions at the GUT scale, we can account for the Higgs mass and keep the lightest neutralino to be as light as about 200 GeV.
The BLSSM is based on the gauge group SU One can implement the seesaw mechanism by adding the right-handed neutrino one per family as required by the anomaly cancellation. In order to count for the right-handed neutrino contributions, we employ the inverse seesaw model in our study. The BLSSM with inverse seesaw includes, in addition to the MSSM particle contents, two SM singlet chiral Higgs superfields χ 1,2 and three sets of SM singlet chiral superfields ν i , s 1i , s 2i [9]. The Superpotential of this model is given by [10] The definition of the parameters appear in W , the corresponding soft SUSY breaking terms and the details of the associate spectrum can be found in Refs. [9,10]. Note that the U (1) Y and U (1) B−L gauge kinetic mixing can be absorbed in the covariant derivative redefinition. Therefore, this mixing provides another source of coupling between MSSM and B − L sectors. Here, we will focus only on the particles involved in the g − 2 loops, namely light neutralino, smuon and chargino, sneutrino.
Considering the additional SM singlet fields of BLSSM with inverse seesaw mentioned above the B − L exten-arXiv:1509.05391v1 [hep-ph] 17 Sep 2015 sion may modify only the neutral sectors of the MSSM only, and hence the chargino and slepton mass matrices remain intact. The 7 × 7 neutralino mass matrix, in the basis: B ,W 0 ,H 0 d ,H 0 u ,B ,χ 1 ,χ 2 , can be found in Ref. [10]. One can easily show that depending on the ratio of the couplings g 1 and g BL , the lightest neutralino could be B-ino (B) or B -ino (B ) like. It is worth noting that in order to account for the Higgs mass, m 0 should be of order TeV, thus µ and µ η , which are determined from electroweak and B − L symmetry breaking conditions, respectively, are also of order TeV scale. Hence, the lightest eigenstates of neutralinos are mostly formed by the neutral gauginos, since the neutral Higgsinos are quite heavy. Now we turn to the sneutrino mass matrix.
, then we can get the CP-odd/even sneutrinos matrices as given in [10].
The supersymmetric contributions to a µ in BLSSM can be split into neutralino and chargino parts as for the MSSM [3], where Z µ L,R , Uμ are the rotation matrices which diagonalize the muon and smuon mass matrices respectively, while Uνi and UνR diagonalize the CP-odd and CP-even sneutrino mass matrices. Note that one can neglect the mixing between slepton families, and consider the smuon mass matrix as 2 × 2 matrix separately from the first and third families. In addition, the mixing between two smuons is proportional to the Yukawa coupling associated with muon, which is of the order ∼ 10 −4 , and hence left-and right-handed smuons are approximately match with the mass eigenstates, and hence, the rotation matrix for the smuons, Uμ can be set to unity in a good approximation. A similar discussion holds for the muon mass matrix diagonalized by Z µ L,R . As seen from Eqs. (5,6), the Bino contribution is in a similar form as obtained in the MSSM, but in BLSSM it has an enhancement by the gauge mixing between U (1) Y and U (1) B−L characterized by the couplingg. It is worth emphasizing the contribution from B −ino (B ). It contributes to a µ through interactions with muon governed by B − L gauge group. In addition, its contribution exhibits an enhancement with the gauge kinetic mixing. Moreover, since it is allowed to be as light as Bino, and even lighter, the lightest neutralino can be formed to be mostlyB orB −B mixing. Thus, one can expect its contribution to be comparable with that from Bino in BLSSM; i.e. N a1 ≈ N a5 numerically. We also present the contribution from the Higgsino component of the Neutralino; however, due to the smallness of Y µ , its contribution is strongly suppressed.
Similarly Eqs. (7,8) reveal the contribution from the chargino expressed in terms of CP-odd and CP-even sneutrino sectors separately. Note that Eqs. (7,8) hold approximately, and one can combine these two sectors in the case of strong mixing between them by summing over j = 1, . . . , 9. The contribution denoted by C L bij are mostly suppressed because of Y µ . On the other hand, C R bij is expected to dominate in the chargino contribution to a µ . It arises from the interactions between muon and sneutrinos through SU(2) interactions as shown in the first and third terms in Eq. (8). The remaining terms are from the Yukawa interactions, and they cannot be neglected, since Y ν is allowed to be of the order O(1) when one employs the inverse seesaw mechanism in BLSSM.
In this regard, the parameters relevant to muon g−2 in BLSSM with inverse seesaw can be listed as MB , mνR 1 , g, and Y ν in addition to those present in MSSM, where MB is the soft SUSY breaking mass term for theB , mνi 1 and mνR 1 are the masses of the lightest CP-even and CP-odd sneutrinos respectively. Besides the mass eigenstates,g, which quantifies the kinetic gauge mixing between U (1) Y and U (1) B−L , and Dirac Yukawa coupling associated with the neutrinos denoted by Y ν are also effective in supersymmetric contributions to a µ . In scanning the parameter space, we have employed the Metropolis-Hastings algorithm as described in [11], and used SARAH [12] and SPheno [13] for the numerical results. The data points collected all satisfy the require- ment of radiative electroweak symmetry breaking. After collecting the data, we impose the mass bounds on all the particles [14]. We have employed the Higgs mass bound as 123 GeV ≤ m h ≤ 127 GeV [15,16], where we take into account about 2 GeV uncertainty in Higgs boson mass due to the theoretical uncertainties in calculation of the minimum of the scalar potential, and the experimental uncertainties in m t and α s . We also employ the gluino mass bound: mg ≥ 1 TeV [17] and the neutral gauge boson Z mass bound: M Z ≥ 2.5 TeV [18].
We present our results first in Fig. 1, where we plot the muon g − 2 as function of lightest neutralino (top) and lightest chargino (bottom). The green points satisfy the mass bounds. As can be seen from the panel ∆a µ − mχ0 1 , muon g − 2 puts a sharp upper bound on the neutralino as mχ0 1 400 GeV. Similarly, ∆a µ − mχ± 1 exhibits a sharp bound on the lightest chargino mass at about 600 GeV. Also the chargino cannot be lighter than about 500 GeV because of the mass bounds we employed in our analysis. Therefore, the sneutrino-chargino contribution is rather suppressed by the chargino mass. Both planes indicate that the best result for muon g−2 is found about 17 × 10 −10 consistently with the constraints. In Fig. 2 we display the correlation between the masses of the light neutralinos;B andB (top) and also the dependence of ∆a µ on tan β (bottom). The color coding is the same as Figure 1 except the plane MB − M 1 in which the yellow points are a subset of green and they satisfy the muon g − 2 condition within 2σ. The unit line in this plane indicates the region with MB = M 1 . As can be seen from the MB −M 1 plane, MB is lighter than M 1 . Hence, Bino and B −ino dominate in muon g − 2 results through the mixing of neutralinos. According to the results, MB is bounded in the interval ∼ 200 − 250 GeV, while M 1 in ∼ 250 − 360 GeV. One can obtain some significant contributions to muon g − 2 in such a narrow mass intervals forB andB with the help of tan β enhancement as shown in the plot of ∆a µ − tan β that muon g − 2 within 2σ requires tan β 45.
It is worth noting that the contributions from MB and M B is further enhanced by their gauge kinetic mixing in smuon-neutralino channel as given in Eqs. (5,6). In Fig.  3 we show the dependence of muon g − 2 on the gauge kinetic mixing characterized withg. The color coding is the same as Figure 1. The figure indicates thatg ∼ −0.1 is favorably for enhancing ∆a µ . In conclusion, we have found that the supersymmetric contribution to muon g − 2 in BLSSM with inverse seesaw mostly relies on the light Bino with M 1 250 GeV, and even lighterB of mass about 180 GeV. In this case the lightest neutralino mass eigenstate is either mostlyB or a linear superposition ofB andB. The tan β dependence of ∆a µ is also represented, and we have found that tan β needs to be 45 in order to raise the supersymmetric contribution such that muon g − 2 results satisfy the measurement within 2σ. In addition, the kinetic gauge mixing couplingg enhance the supersymmetric contribution to muon g − 2. We have revealed in our analysis that muon g − 2 favorsg ≈ −0.1. The gauge kinetic mixing affects the contribution from smuon-neutralino channel, and with its enhancement, the smuons can have the mass about a TeV, and hence the higher m 0 values can still yield results compatible with muon g − 2 as well as consistent with the 125 GeV Higgs boson mass.