The extended BLMSSM with a 125 GeV Higgs boson and dark matter

To extend the BLMSSM, we not only add exotic Higgs superfields $(\Phi_{NL},\varphi_{NL})$ to make the exotic lepton heavy, but also introduce the superfields($Y$,$Y^\prime$) having couplings with lepton and exotic lepton at tree level. The obtained model is called as EBLMSSM, which has difference from BLMSSM especially for the exotic slepton(lepton) and exotic sneutrino(neutrino). We deduce the mass matrices and the needed couplings in this model. To confine the parameter space, the Higgs boson mass $m_{h^0}$ and the processes $h^0\rightarrow \gamma\gamma$, $h^0\rightarrow VV, V=(Z,W)$ are studied in the EBLMSSM. With the assumed parameter space, we obtain reasonable numerical results according to data on Higgs from ATLAS and CMS. As a cold dark mater candidate, the relic density for the lightest mass eigenstate of $Y$ and $Y'$ mixing is also studied.


I. INTRODUCTION
The total lepton number (L) and baryon number (B) are good symmetries because neutrinoless double beta decay or proton decay has not been observed. In the standard model(SM), L and B are global symmetries [1]. However, the individual lepton numbers L i = L e , L µ , L τ are not exact symmetries at the electroweak scale because of the neutrino oscillation and the neutrinos with tiny masses [2]. In the Universe, there is matter-antimatter asymmetry, then the baryon number must be broken.
With the detection of the light Higgs h 0 (m 0 h = 125.1GeV) [3], the SM succeeds greatly and the Higgs mechanism is compellent. Beyond the SM, supersymmetry [4] provides a possibility to understand the light Higgs. The minimal supersymmetric extension of the SM(MSSM) [5] is one of the favorite models, where the light Higgs mass at tree level is m tree h = m Z | cos 2β| [6]. The one loop corrections to Higgs mass mainly come from fermions and sfermions, that depend on the virtual particle masses and the couplings with the Higgs.
There are many papers about the gauged B and L models, although most of them are non-supersymmetric [7]. Extending MSSM with the local gauged B and L, one obtains the so called BLMSSM, which was proposed by the authors in Ref. [6]. The proton remains stable, as B and L are broken at the TeV scale. Therefore, a large desert between the electroweak scale and grand unified scale is not necessary. In BLMSSM, the baryon number is changed by one unit, at the same time the lepton number is broken in an even number. R-parity in BLMSSM is not conserved, and it can explain the matter-antimatter asymmetry in the Universe. There are some works for Higgs and dark matters [8] in the BLMSSM [9,10]. In the framework of BLMSSM, the light Higgs mass and the decays h 0 → γγ and h 0 → V V, V = (Z, W ) are studied in our previous work [10]. Some lepton flavor violating processes and CP-violating processes are researched with the new parameters in BLMSSM [11].
In BLMSSM, the exotic leptons are not heavy, because their masses just have relation with the parameters Y e 4 υ d , Y e 5 υ u . Here υ u and υ d are the vacuum expectation values(VEVs) of two Higgs doublets H u and H d . In general, the Yukawa couplings Y e 4 and Y e 5 are not large parameters, so the exotic lepton masses are around 100 GeV. The light exotic leptons may lead to that the BLMSSM is excluded by high energy physics experiments in the future. To obtain heavy exotic leptons, we add two exotic Higgs superfields to the BLMSSM, and they are SU (2) singlets Φ N L and ϕ N L , whose VEVs are υ N L andῡ N L [12]. The exotic leptons and the superfields Φ N L , ϕ N L have Yukawa couplings, then υ N L andῡ N L give contributions to the diagonal elements of the exotic lepton mass matrix. So the exotic leptons turn heavy and should be unstable. In the end, the super fields Y and Y ′ are also introduced. At tree level, there are couplings for lepton-exotic lepton-Y (Y ′ ). It is appealing that this extension of BLMSSM produces some new cold dark matter candidates, such as the lightest mass egeinstate of Y and Y ′ mixing. The four-component spinorỸ is made up of the superpartners of Y and Y ′ . In this extended BLMSSM(EBLMSSM), we study the lightest CP even Higgs mass with the one loop corrections. The Higgs decays h 0 → γγ and h 0 → V V, V = (Z, W ) are also calculated here. Supposing the lightest mass eigenstate of Y and Y ′ mixing as a cold dark matter candidate, we study the relic density.
After this introduction, in Section 2, we introduce the EBLMSSM in detail, including the mass matrices and the couplings different from those in the BLMSSM. The mass of the lightest CP-even Higgs h 0 is deduced in the Section 3. The Section 4 is used to give the formulation of the Higgs decays h 0 → γγ, h 0 → V V, V = (Z, W ) and dark matter relic density. The corresponding numerical results are computed in Section 5. The last section is used for the discussion and conclusion.

II. EXTEND THE BLMSSM
The local gauge group of the BLMSSM [6] In the BLMSSM, the exotic lepton masses are obtained from the Yukawa couplings with the two Higgs doublets H u and H d . The VEVs of H u and H d are υ u and υ d with the relation GeV. Therefore, the exotic lepton masses are not very heavy, though they can satisfy the experiment bounds at present. In the future, with the development of high energy experiments, the experiment bounds for the exotic lepton masses can improve in a great possibility. Therefore, we introduce the exotic Higgs superfields Φ N L and ϕ N L with nonzero VEVs to make the exotic lepton heavy. The heavy exotic leptons should be unstable, then the superfields Y, Y ′ are introduced accordingly. These introduced superfields lead to tree level couplings for lepton-exotic lepton-Y (Y ′ ).
In EBLMSSM, we show the superfields in the Table I. The superpotential of EBLMSSM is shown here W M SSM is the superpotential of MSSM. W B and W X are same as the terms in BLMSSM [10].
W Y includes the terms beyond BLMSSM, and they include the couplings of lepton-exotic lepton-Y (l I − L ′ − Y ). Therefore, the heavy exotic leptons can decay to leptons and mass eigenstates of Y and Y ′ mixing whose lighter one can be a dark matter candidate. From c 5Ŷ is a matrix and has non-zero elements relating with lepton flavor, this term can enhance the lepton flavor violating effects. In the whole, W Y enriches the lepton physics to a certain degree, and these subjects will be researched in our latter works.
Because of the introduction of the superfields Φ N L , ϕ N L , Y and Y ′ , the soft breaking terms are written as Here L BLM SSM sof t is the soft breaking terms of BLMSSM, whose concrete form is in our previous work [10]. The SU(2) L doublets H u , H d acquire the nonzero VEVs υ u , υ d . The SU(2) L singlets Φ B , ϕ B , Φ L , ϕ L , Φ N L , ϕ N L obtain the nonzero VEVs υ B , υ B , υ L , υ L , υ N L , υ N L respectively.
Here, we define tan The VEVs of the Higgs satisfy the following equations Here, the Eqs. (8) and (9) are similar as the corresponding equations in BLMSSM, but Eqs. (8) and (9) have relation with the new parameters υ N L and υ N L . We obtain the new Eqs.(10) and (11) through ∂V ∂Φ NL and ∂V ∂ϕ NL , with V denoting the Higgs scalar potential.
Here we deduce the mass matrices in the EBLMSSM. Compared with BLMSSM, the superfields Φ N L and ϕ N L are introduced and they give corrections to the mass matrices of the slepton, sneutrino, exotic lepton, exotic neutrino, exotic slepton and exotic sneutrino.
That is to say, in EBLMSSM, the mass matrices of squark, exotic quark, exotic squark, baryon neutralino, MSSM neutralino, X andX are same as those in the BLMSSM, and their concrete forms can be found in our previous works [13,14]. Though the mass squared matrices of slepton and sneutrino in EBLMSSM are different from those in BLMSSM, we can obtain the slepton and sneutrino mass squared matrices in EBLMSSM easily just using the replacement υ 2 L − υ 2 L → V 2 L for the BLMSSM results. In the BLMSSM, the issue of Landau pole has been discussed in detail by the authors of Ref. [6]. Their conclusion is that there are no Landau poles at the low scale due to the new families. In EBLMSSM, the parts of quark (squark), exotic quark (exotic squark) are same as those in BLMSSM. Therefore, the Landau pole conditions for the Yukawa couplings of quark (squark), exotic quark (exotic squark) have same behaviors of BLMSSM.
The added superfields (Φ N L , ϕ N L , Y, Y ′ ) do not have couplings with the gauge fields of So the characters of gauge couplings g 1 , g 2 , g 3 and g B in BLMSSM and EBLMSSM are same.
The different parts between BLMSSM and EBLMSSM are the terms including In BLMSSM, the exotic lepton masses are not heavy, because they obtain masses only from H u and H d . The VEVs of Φ N L and ϕ N L are υ N L andῡ N L , that can be large parameters.
So, the EBLMSSM exotic leptons are heavier than those in BLMSSM.
The mass matrix for the exotic leptons reads as Obviously, υ N L and υ N L are the diagonal elements of the mass matrix in the Eq. (12). It is easy to obtain heavy exotic lepton masses with large υ N L and υ N L . If we take υ N L and υ N L as zero, the mass matrix is same as that in BLMSSM. In fact, our used values of υ N L and υ N L are at TeV order, which produce TeV scale exotic leptons. Heavy exotic leptons have strong adaptive capacity to the experiment bounds. The exotic neutrinos are four-component spinors, whose mass matrix is Similar as the exotic lepton condition, heavy exotic neutrinos are also gotten.
In BLMSSM, the exotic sleptons of 4 generation and 5 generation do not mix, and their mass matrices are both 2 × 2. In EBLMSSM, the exotic sleptons of 4 generation and 5 generation mix together, and their mass matrix is 4 × 4. With the base (ẽ 4 ,ẽ c * 4 ,ẽ 5 ,ẽ c * 5 ), we show the elements of exotic slepton mass matrix M 2Ẽ in the following form.
In Eq. (14), the non-zero terms and M 2Ẽ (ẽ c 4ẽ 5 ) are the reason for the exotic slepton mixing of generations 4 and 5. These mixing terms all include the parameters υ N L andῡ N L . It shows that this mixing is caused basically by the added Higgs superfields Φ N L and ϕ N L . Using the matrix ZẼ, we obtain mass eigenstates with the In the same way, the exotic sneutrino mass squared matrix is also obtained For the exotic sneutrino, the mixing of generations 4 and 5 is similar as that of exotic slepton. In the base (ν 4 ,ν c * 4 ,ν 5 ,ν c * 5 ), we get the mass squared matrix of the exotic sneutrino, and obtain the mass eigenstates by the matrix ZÑ through the formula

B. The lepton neutralino mass matrix in EBLMSSM
In EBLMSSM, the superfields (Φ L , ϕ L , Φ N L , ϕ N L ) have their SUSY superpartners . They mix with λ L , which is the superpartner of the new lepton type gauge boson Z µ L . Therefore, we deduce their mass matrix in the base The lepton neutralino mass egeinstates are four-component spinors X 0 In BLMSSM, there are no ψ Φ NL , ψ ϕ NL , and the base of lepton neutralino is (iλ L , ψ Φ L , ψ ϕ L ), whose mass matrix is 3 × 3. EBLMSSM extends this matrix to 5 × 5 including the BLMSSM results.

C. The Higgs superfields and Y in EBLMSSM
The superfields Φ L , ϕ L , Φ N L , ϕ N L mix together and form 4×4 mass squared matrix, which is larger than the corresponding 2 × 2 mass matrix in the BLMSSM. Diagonalizing the mass squared matrix, four CP even exotic Higgs are obtained.
We use Zφ L to diagonalize the mass squared matrix in Eq. (18), and the relation between mass eigenstates and the comments are In EBLMSSM, the conditions for the exotic CP odd Higgs P 0 L ,P 0 L are same as those in BLMSSM, and they do not mix with the added exotic CP odd Higgs P 0 N L ,P 0 N L . Here, we show the mass squared matrix for the added exotic CP odd Higgs P 0 N L ,P 0 N L .
The scalar superfields Y and Y ′ mix, and their mass squared matrix is deduced here.
This condition is similar as that of X and X ′ , then the lightest mass egeinstate of Y and Y ′ can be a candidate of the dark matter. With S Y = g 2 L (2 + L 4 )V 2 L , the concrete form for the mass squared matrix is shown here. To obtain mass eigenstates, the matrix Z Y is used through the following formula, with the supposition m 2 The superpartners of Y and Y ′ form four-component Dirac spinors, and the mass term for superfieldsỸ is shown as The spinorỸ and the mixing of superfields Y, Y ′ are all new terms beyond BLMSSM, that add abundant contents to lepton physics and dark matter physics.
In Eq.(23), different from BLMSSM, there are new terms ( In this coupling, the new terms beyond BLMSSM are −( The h 0 −L −L coupling has the same form as that in BLMSSM. While, the h 0 −ν −ν coupling gets corrected terms, but these terms are suppressed by the tiny neutrino Yukawa Here, s W (c W ) denotes sin θ W (cos θ W ), with θ W representing the weak-mixing angle. The concrete form of B 2 R is in Ref. [5].

E. The couplings with Y
For the dark matter candidate Y 1 , the necessary tree level couplings are deduced in EBLMSSM. We show the couplings (lepton-exotic lepton-Y ) and (neutrino-exotic neutrino- The new gauge boson Z L couples with leptons, neutrinos and Y , whose concrete forms are ϕ L gives masses to the light neutrinos trough the see-saw mechanism and Φ L , ϕ L , Φ N L , ϕ N L mix together producing lepton Higgs H 0 L . Then the couplings of H 0

III. THE MASS OF h 0
Similar as BLMSSM, in EBLMSSM the mass squared matrix for the neutral CP even Higgs are studied, and in the basis (H 0 d , H 0 u ) it is written as where M 2 11 , M 2 12 , M 2 22 are the tree level results, whose concrete forms can be found in Ref. [10] The   [15,16] with x a = m 2 h 0 /(4m 2 a ). Here, q and q ′ are quark and exotic quark. While,q andq ′ denote squark and exotic squark. The concrete expressions for g h 0 qq , g h 0 q ′ q ′ , g h 0qq , g h 0q′q′ (i = 1, 2) are in literature [10]. The functions A 1/2 (x) and A 0 (x) are [16] The decay h 0 → γγ obtains contributions from loop diagrams, and the leading order contributions are from the one loop diagrams. In the EBLMSSM, the exotic quark(squark) and exotic lepton (slepton) give new corrections to the decay width of h 0 → γγ. Different from BLMSSM, the exotic leptons in EBLMSSM are more heavy and the exotic sleptons of the 4 and 5 generations mix together. These parts should influence the numerical results of the EBLMSSM theoretical prediction to the process h 0 → γγ to some extent.
The decay width of h 0 → γγ can be expressed as [17] where The formulae for h 0 → ZZ, W W are with g h 0 ZZ = g h 0 W W and F (x) is given out in Ref [18,19]. The observed signals for the diphoton and ZZ, W W channels are quantified by the ratios R γγ and R V V , V = (Z, W ), whose current values are R γγ = 1.16 ± 0.18 and R V V = 1.19 +0.22 −0.20 [20].

B. Dark matter Y
In BLMSSM, there are some dark matter candidates such as: the lightest mass eigenstate of XX ′ mixing,X the four-component spinor composed by the super partners of X and X ′ .
They are studied in Ref. [9]. In EBLMSSM, the dark matter candidates are more than those in BLMSSM, because the lightest mass eigenstate of Y Y ′ mixing andỸ are dark matter candidates. After U(1) L is broken by Φ L and Φ N L , Z2 symmetry is left, which guarantees the stability of the dark matters. There are only two elements (1,-1) in Z2 group. This symmetry eliminates the coupling for the mass eigenstates of Y Y ′ mixing with two SM particles. The condition for X is similar as that of Y , and it is also guaranteed by the Z2 symmetry.
In this subsection, we suppose the lightest mass eigenstate of Y Y ′ mixing in Eq.(21) as a dark matter candidate, and calculate the relic density. So we summarize the relic density constraints that any WIMP candidate has to satisfy. The interactions of the WIMP with SM particles are deduced from the EBLMSSM, then we study its annihilation rate and its relic density Ω D by the thermal dynamics of the Universe. The annihilation cross section σ(Y 1 Y * 1 → anything) should be calculated and can be written as σv rel = a + bv 2 rel in the Y 1 Y * 1 center of mass frame. v rel is the twice velocity of Y 1 in the Y 1 Y * 1 c.m. system frame. To a good approximation, the freeze-out temperature(T F ) can be iteratively computed from [8] x with x F ≡ m D /T F and m D = m Y 1 representing the WIMP mass. M P l = 1.22 × 10 19 GeV is the Planck mass and g * is the number of the relativistic degrees of freedom with mass less than T F . The density of cold non-baryonic matter is Ω D h 2 = 0.1186 ± 0.0020 [20], whose formula is simplified as To obtain a and b in the σv rel , we study the Y 1 Y * 1 dominate decay channels whose final states are leptons and light neutrinos: Nα=νe,νµ,ντ Nα=νe,νµ,ντ 1 96π

V. NUMERICAL RESULTS
A. h 0 decays and m A 0 , m H 0 In this section, we research the numerical results. For the parameter space, the most strict constraint is that the mass of the lightest eigenvector for the mass squared matrix in Eq.(29) is around 125.1 GeV. To satisfy this constraint, we use m h 0 = 125.1 GeV as an input parameter. Therefore, the CP odd Higgs mass should meet the following relation.
To embody the exotic squark corrections, we calculate the results versus A

VI. DISCUSSION AND CONCLUSION
Considering the light exotic lepton in BLMSSM, we add exotic Higgs superfields Φ N L and ϕ N L to BLMSSM in order to make the exotic leptons heavy. Light exotic leptons may be excluded by the experiment in the future. On the other hand, heavy exotic leptons should not be stable. So we also introduce the superfields Y and Y ′ to make exotic leptons decay quickly. The lightest mass eigenstate of Y and Y ′ mixing mass matrix can be a dark matter candidate. Therefore, the exotic leptons are heavy enough to decay to SM leptons and Y at tree level. We call this extended BLMSSM as EBLMSSM, where the mass matrices for the particles are deduced and compared with those in BLMSSM. Different from BLMSSM, the exotic sleptons of 4 and 5 generations mix together forming 4 × 4 mass squared matrix.
EBLMSSM has more abundant content than BLMSSM for the lepton physics.
To confine the parameter space of EBLMSSM, we study the decays h 0 → γγ and Based on the supposition that the lightest mass eigenstate Y 1 of Y and Y ′ mixing possesses the character of cold dark matter, we research the relic density of Y 1 . In our used parameter space, Ω D h 2 of Y 1 can match the experiment bounds. EBLMSSM has a bit more particles and parameters than those in BLMSSM. Therefore, EBLMSSM possesses stronger adaptive capacity to explain the experiment results and some problems in the theory. In our later work, we shall study the EBLMSSM and confine its parameter space to move forward a single step.