On the Higgs-like boson in the minimal supersymmetric 3-3-1 model

It is imperative that any proposal of new physics beyond the standard model possesses a Higgs-like boson with 125 GeV of mass and couplings with the standard particles that recover the branching ratios and signal strengths as measured by CMS and ATLAS. We address this issue within the supersymmetric version of the minimal 3-3-1 model. For this we develop the Higgs potential with focus on the lightest Higgs provided by the model. Our proposal is to verify if it recovers the properties of the Standard Model Higgs. With respect to its mass, we calculate it up to one loop level by taking into account all contributions provided by the model. In regard to its couplings, we restrict our investigation to couplings of the Higgs-like boson with the standard particles, only. We then calculate the dominant branching ratios and the respective signal strengths and confront our results with the recent measurements of CMS and ATLAS. As distinctive aspects, we remark that our Higgs-like boson intermediates flavor changing neutral processes and has as signature the decay t→h+c\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t \rightarrow h+c$$\end{document}. We calculate its branching ratio and compare it with current bounds. We also show that the Higgs potential of the model is stable for the region of parameter space employed in our calculations.


Introduction
The discovery of a scalar state by the ATLAS and CMS at CERN [1,2] compatible with the standard model Higgs boson was a milestone in particle physics since it was the missing part of the theory. Nowadays we not only know the Higgs mass with precision, M h = 125.09 ± 0.21(stat.) ± 0.11(syst.) GeV, but also know that its signal strengths are found consistent with the standard model (SM) predictions [3]. However, the present experimental uncertainties allow a e-mail: cpires@fisica.ufpb.br us to interpret such observed state as belonging to a more complex theory that encompasses the standard one.
Supersymmetry is the most popular proposal of new physics and all supersymmetric models predict at least one light Higgs-like boson. For example, the minimal supersymmetric standard model (MSSM) provides a tree level contribution to the mass of the Higgs-like boson that does not surpass the neutral gauge boson Z 0 mass. Consequently the amount of radiative corrections required to complete 125 GeV demands stops with mass in a region of values that threaten the naturalness principle [4]. In view of this it turns out imperative to examine the Higgs sector of any supersymmetric phenomenological model with the aim of investigating the lightest Higgs concerning its mass and couplings with the standard particles.
In this work we develop the scalar sector of the supersymmetric version of the minimal SU (3) C × SU (3) L × U (1) N (331) model. The motivations to study 331 models are various as, for example, explanation of family replication [5,6], electric charge quantization [7][8][9], strong CP-problem [10][11][12], incorporation of inflation [13] and so forth. As we show in this paper, its SUSY versions, besides solving the hierarchy problem, provides a tree level contribution to the lightest Higgs boson that may surpass 100 GeV. This is nice in what concerns the naturalness principle since now 125 GeV Higgs is compatible with stop with mass below the TeV scale.
The minimal supersymmetric 331 model was firstly developed in [14]. For subsequent works, see [15][16][17]. As the Higgs sector of the minimal 331 model [5,6] requires at least three Higgs triplets to generate correctly the masses of all massive particles and be phenomenologically viable [18][19][20], then its supersymmetric version, what we call the minimal SUSY331 model, must involve six Higgs triplets as required to cancel anomalies. On assuming that the two triplets χ and χ decouple from the other four η, ρ, η and ρ , and that η and ρ are inert, then after spontaneous breaking of the 331 symmetry to the standard one, we obtain the potential for two Higgs doublets modified by cubic invariant terms. On developing such potential, we obtain the following approximate expression for the lightest scalar, the Higgs-like boson candidate, of the model where M Z is the mass of the standard neutral gauge boson, θ w is the electroweak mixing angle and the angle β is such that tan β = v η v ρ . According to this expression the tree level contribution to the lightest Higgs may attain 174 GeV for cos(β) = 1. This is so because the symmetries of the model allow the superpotential to have cubic invariant terms formed with the scalar triplets of the model. This acts like in the NMSSM [14]. In view of this, it is compulsory to review the Higgs sector of this model in the light of the recent discoveries of the LHC, and verifying how realistic the expression above is and define precisely the profile of the lightest Higgs in the minimal SUSY331 model.
Our contribution to this discussion restricts to obtain the mass matrix that contains the neutral Higgs, diagonalize it in the most general way and obtain the eigenvalue and the eigenvector that corresponds to the lightest Higgs. We investigate the behavior of its mass by taking into account one loopcorrections. We also obtain the branching ratios for the dominant processes and calculate the respective signal strengths. It is notorious that our Higgs-like boson candidate intermediates flavor changing processes. We discuss such processes but focusing on the dominant ones that are t → h + c and t → h +u. Finally, we guarantee that the region of parameter space employed in our calculations leads to a stable potential.
The paper is divided in the following way: In Sect. 2 we present the main ingredients of the model. Next, in Sect. 3, we develop the Higgs sector with focus on the lightest Higgs where we calculate its mass up to one-loop level. In Sect. 4 we calculate its branching ratios and the respective signal strengths. In Sect. 5 we calculate the flavor changing processes. In Sect. 6 we address the stability of the potential and, finally, we conclude in Sect. 7.

The main ingredients of the model
In this section we present the core of the minimal SUSY331 model relevant for what follows. In regard to the leptonic sector, the superfields of each generation are arranged in triplet of superfields according to the following transformation by the 3-3-1 symmetrŷ where a = 1, 2, 3 represents the family index for the usual three generations of leptons.
In the Hadronic sector, the superfields of the third generation comes in the triplet representation and the superfields of the other two come in anti-triplet representations of SU (3) L , as a requirement for anomaly cancellation. They are given bŷ where α = 1, 2.
The scalar sector of the 3-3-1 model, responsible for the spontaneously broken gauge symmetry, is composed by three scalar triplets. In the supersymmetric version, anomaly cancellation requires we double these fields. Thus, the scalar sector of the minimal SUSY331 is composed by the following superfieldŝ whereη ∼ (1, 3, 0) ,χ ∼ (1, 3, −1) ,ρ ∼ (1, 3, 1), and We assume that the neutral scalars η, η , ρ, ρ , χ and χ develop nonzero VEV according to These VEVs lead to the following gauge symmetry breaking pattern With the breaking of the gauge symmetry by this set of VEVs all the massive particles, including the supersymmetric ones, receive mass. What matters for us here are the scalar and gauge boson masses. Concerning the gauge bosons, they are composed by the standard gauge boson, γ, Z 0 and W ± , one new neutral massive gauge bosons Z , two doubly charged gauge bosons U ±± , and two simply charged gauge bosons V ± . The charged gauge bosons gain the following mass expressions while the neutral gauge bosons have masses with v 331 lying in the TeV scale. After all this, we are ready to build up the superpotential of the model. The superpotential that respects the gauge symmetry and R-parity is composed by the following terms where α, β = 1, 2 and i = 1, 2, 3. We could be economical and resort to a Z 3 symmetry so as to avoid the bilinear terms in the superpotential above. However, in order to be general enough, we do not follow this path here. We call the attention to the fact that the last term in the first line of the superpotential is an effective 5-D operator. It will generate the masses of the charged leptons. This point has been well developed in many previous papers, but it is appropriate to recall it here. The highest energy scale where the model is found to be perturbatively reliable is about = 4 − 5 TeV. Hence, effective operators may be required to generate corrections to the mass of some charged fermions. We make use of this here with the aim of avoiding the addition of scalars sextets to the model.
Up to this point the masses of the standard particles are equal to the masses of their superpartners. As usual, in phenomenological supersymmetric models SUSY must be broken so as to provide a reasonable shift between ordinary par-ticles and their supersymmetric partners. In this work we assume that SUSY is broken explicitly through the following set of soft breaking terms that are invariant under the symmetries assumed here where α, β = 1, 2 and i = 1, 2, 3. λ a c are the gluinos, λ a are gauginos associated to SU (3) L (in both cases a = 1, . . . , 8 is the gauge group index) and λ is the gaugino associated to U (1) N . The scalar supersymmetric partners of fermion fields, f , are denoted byf , while the remaining fields are self-evident. For simplicity, we will take m 2 As it is usual in SUSY models, which involves a large number of free parameters, simplifications are necessary in order to easily get a better view of the big picture. As simplification we use the following parametrization where κ's are the corresponding Yukawa couplings, M s is the SUSY scale, mt is the mass of the stop and X t is the equivalent of the soft trilinear coupling of the stops. We also assume that all soft left and right masses, m SU SY L and m SU SY R , are equals. With all this we are ready to further explore the Higgs physics.

Higgs physics in the minimal SUSY331 model
In supersymmetric models the Higgs potential receives contributions from three different sources, F-term, D-term and the soft SUSY breaking terms, that adds up to form the scalar potential V = V F + V D + V so f t . These contributions are, respectively where t A are the Gell-Mann matrices. The fields are assumed to be shifted in the usual way and the set of minimum conditions are given in Appendix A. That set of constraint equations enable us to obtain the mass matrix of the scalars of the model. In this work we focus exclusively in the lightest neutral scalar which is expected to recover the properties of the Standard Model Higgs in what concerns its mass, branching ratios and signal strengths. For this we first obtain the 6 × 6 mass matrix, M 2 R H , associated to the CP-even scalars R η , R η , R ρ , R ρ , R χ and R χ . It is a very complex, and not illuminating, mass matrix, then we do not show it here. Next, we diagonalize it by a rotating mixing matrix U R such that the diagonal mass matrix M 2 The lightest eigenstate, let us call it h i , must be identified as the Higgs-like boson. This means that it must have a mass of 125 GeV and its eigenvector, h i = (U R ) i j R j , must recover the standard Higgs couplings. In practical terms, to know the eigenvector h i means to obtain the set of entries (U R ) i j . Such entries, when combined with the adequate couplings, must recover the existing branching ratios and signal strengths of the Higgs. From now on we refer to the eigenvector associated to the Higgs-like boson as h.
Due to the large number of parameters involved in the diagonalization of M 2 R H , an analytical approach is unviable. In view of this our approach will be completely numerical.
In what follows we make use of powerful tools like the packages Sarah 4.10.2 [21][22][23], Spheno 4.0.2 [24,25] and SSP 1.2.3 [26] ( We include in Appendix B and C the implemented model in Sarah and the Spheno input, respectively). Throughout this work we employ the following routine: we implement the model in the Sarah 10.2 and export it to the Spheno 4.0.2, where we make all the numerical calculations. The SSP helps in making the scan of the parameters. Hereafter all numerical computations done in this paper follow this routine, including the diagonalization of this mass matrix.

Tree level contribution
Differently from the MSSM, where the tree level contribution to the Higgs mass does not surpass 91 GeV, in the minimal SUSY331 model things are a little different. The expression in Eq. (1) predicts that the lightest scalar, the Higgs-like boson, may attain a mass of 125 GeV already at tree level. It is important to note that the approximations assumed in Eq. (1) do not take into account the stability of the potential. Thus, it becomes mandatory, and challenging, to clarify if the general case, which involves the diagonalization of a 6 × 6 mass matrix, agrees with the prediction in Eq. (1).
The results displayed in Fig. 1 is for the range of parameters Observe that the tree level contribution may attain 115 GeV for tan(β) = 2 which is in a good agreement with the predictions of the Eq. (1). In other words, our results confirm that Eq. (1) is a valid approximation for the tree level contribution to the Higgs-like boson mass in the minimal SUSY331 model. Thus, the robustness of the Higgs-like boson mass at tree level is dictated by tan(β) such that the smaller tan(β), the larger tree level Higgs-like boson mass.
Nevertheless, we know that loop correction to the Higgs mass are considerable, enforcing us to conclude that tan(β) ≥ 2 and that the model may support a stop with mass below the TeV scale as required by naturalness principle.
For sake of completeness, we obtain the eigenvector associated to the eigenvalue 125 GeV for the following set of values of the parameters tan(β) = 2.7, tan(β η ) = 10, tan(β ρ ) = 10, In this case, the eigenvector corresponding to the 125 GeV Higgs is See that the eigenvector h is composed mainly by R η and R ρ . For any other choice of the set of values of the parameters, the eigenvector will keep being dominantly a composition of R η and R ρ .

One loop level contribution
In In Fig. 2 we plot M h versus X t . As we can see in that plot, small X t is allowed for large value of v 331 , while large X t is allowed by small values of v 331 . Moreover, we highlight that We also obtained the 2 × 2 mass matrix for the stops and diagonalized it with the Spheno package. The result is displayed in Fig. 3 where we show the behavior of M s with X t . Note that the larger v 331 , the smaller M s . In Fig. 4 (left) we show the behavior of M h with M s for the specific case X t = 0, and in Fig. 4 (right) we show the behavior of M h with v 331 for X t = 0. Perceive that the model provides easily a 125 GeV Higgs boson for M s and X t , both, below the TeV scale.
As we can see, the minimal SUSY331 model naturally recovers the Standard Model Higgs boson mass, but we know that it is different from the MSSM one because of the peculiarities of the minimal SUSY331 model. Thus, in order to conclude definitely that our Higgs-like boson recovers all the observed properties of the Standard Model Higgs, we need to calculate its branching ratios and the respective signal strengths and confront the results with the experimental data measured by LHC.

Branching ratios and signal strengths
In what concerns the Standard Model Higgs, any extension of the SM must possess a scalar with 125 GeV of mass and couplings with the Standard Model particles that fit the measured branching ratios and signal strengths [27]. This is the reason why we wonder if the Higgs-like boson discussed here recovers the Higgs branching ratios and signal strengths as observed by LHC.
As we know, the Higgs prefers to decay into pairs of bb, W W * , ττ , Z Z * and γ γ . We restrict our investigation for these channels. Our results are shown in Fig. 5. Perceive that the predictions are in perfect agreement with the values measured by ATLAS [28] and CMS [29]. In this point of the work, we are secure in affirming that the minimal SUSY331 model is privileged in what concerns Higgs physics since it possesses a Higgs-like boson with mass and couplings with the standard particles that recovers its observed properties and respecting the naturalness principle.

Flavor changing neutral current as signature of our Higgs-like boson
As we saw above, the lightest Higgs of the minimal SUSY331 model fulfills the conditions to be a SM-like Higgs. However, it is important to recall that our Higgs-like bosom is different from the SM-like one in many aspects with the major one being that it intermediates flavor changing processes involv- ing quarks already at tree level. This is a consequence of quark masses having more than one source, as depicted in the superpotential of the model. According to the last line of the superpotential, and the fact that our Higgs-like boson is mostly a composition of R η and R ρ , the dominant flavor changing processes are those involving the third family of quarks. We studied all the possible processes here and obtained that the decays h → bs , bd are very suppressed even in relation to the loop predictions of the SM [30]. Thus, the significant processes are those involving the Yukawa interactions t − h − u and t − h − c. For this case, the signature of our Higgs-like boson is through top quark decays via Higgs-mediated flavor-changing processes [31]. The behavior of the branching ratio of these decays as function of v 331 is presented in Fig. 6. Although our results are far below the excluded region, we hope this to be probed in the next generation of collider.

Stability of the potential
In this section we check if that region of parameter space considered throughout this paper is so that the potential of the model is stable for such value of the parameters. In practical terms, we have to be sure that the minimum of the potential we are considering is in fact the global one. We check this by means of the Vevacious 1.2.02 Package [32,33]. We export the model implemented in the Sarah package to the Vevacious one and scanned that region of the parameter space used in the results above. Our result is presented in the plane tan(β) versus v 331 displayed in Fig. 7 .
Perceive that, except for a small region of the parameter space presenting long lived behavior, but with decay time larger than the age of our Universe (blue region), major part of the region of the parameter space considered in this work is associated to a stable minimum of the potential.

Conclusions
In this work we developed the Higgs sector of the minimal SUSY331 model with the focus on the lightest Higgs. We obtained its mass and couplings with the standard particles. With regard to its mass, tree level contribution may surpass the usual MSSM prediction and attain more than 100 GeV. Consequently, a 125 GeV Higgs will demand feeble loop corrections. In our calculations we showed that a 125 GeV Higgs is compatible with stop with mass below the TeV scale running in the loops. This is a remarkable result concerning the naturalness principle. Although Higgs mass is the most important aspect of the Higgs physics, a complete work demands we extend our investigation to the Higgs couplings with the standard particles. We performed such an analysis and confirmed that they recover the observed pattern of branching ratios and signal strengths for its dominant processes. Finally we discussed the signature of our Higgs-like boson which manifests in the form of flavor changing processes intermediated by the Higgs. The most relevant ones are the top decay channels t → h + u and t → h + c. We studied the behavior of these processes and concluded that they are out of reach of the LHC but, perhaps, may be probed in future colliders. We also checked if the region of parameter space considered in this work is compatible with the stability of the potential. We obtained that, except for a small region of the parameter space presenting long lived behavior, but with decay time larger than the age of our Universe, major part of it belongs to those points where the potential has a stable global minimum. All of this leads us to conclude that the minimal SUSY331 model is an interesting supersymmetric model for particle physics.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm 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 .