On the Higgs-like boson in the Minimal Supersymmetric 3-3-1 Model

It is imperative that any proposal of new physics 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 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 then argue that its signature is the decay $t \rightarrow h+c$. We calculate its branching ratio and compare it with current bounds. We also show that the potential is stable for the region of parameter space employed in our calculations.


I. 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 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. Consequentely 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 [10] [11], strong CP-problem [12][13] [14], incorporation of inflation [15] 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 [16]. For subsequent works, see [17][18] [19]. 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 [7][8] [9], 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 superpotential naturally contains cubic invariant terms that act like in the NMSSM [16]. 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 loop-corrections. We also obtain the branching ratios for the dominant processes and calculate the respective signal strengths. It is notorious that our Higgs-like boson intermediates flavor changing processes. We discuss such processes but focus 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 Sec. II we present the main ingredients of the model. Next, in Sec. III, we develop the Higgs sector with focus on the lightest Higgs where we calculate its mass up to one-loop level. In Sec. IV we calculate its branching ratios and the respective signal strengths. In Sec. V we calculate the flavor changing processes. In Sec. VI we address the stability of the potential and, finally, we conclude in Sec. VII.

II. 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.
We assume that the neutral scalars η, η , ρ, ρ , χ and χ develop nonzero VEV according 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, 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 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 particles 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.

III. 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 sof t . These contributions are, respectively where t A are the Gell-Mann matrices.
In this work we focus exclusively in the lightest neutral scalar which is expected to recover the properties of the standard 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 The physical eigenstates relates to the symmetrical ones by 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 ) ij 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 ) ij . 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 diagonalisation of M 2 R H , an analytical approach is unviable, and we turn it completely numerical. In what follows we make use of powerful tools like the packages Sarah [20][21][22], Spheno [23,24] and SSP [25]. Throughout this work we employ the following routine: we implement the model in the Sarah and export it to the Spheno, where we make all the numerical calculations. The SSP helps in making the scan of the parameters. Hereafter all numerical computatiions done in this paper will follow this routine, including the diagonalisation of this mass matrix.

A. 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 diagonalisation 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 In this case, the eigenvector corresponding to the 125 GeV Higgs is 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, observe that we  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 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 Higgs, we need to calculate its branching ratios and the respective signal strengths and confront the results with the experimental data measured by LHC.

IV. BRANCHING RATIOS AND SIGNAL STRENGTHS
In what concerns the standard Higgs, any extension of the SM must possess a scalar with 125 GeV of mass and couplings with the standard particles that fit the measured branching ratios and signal strengths [26]. 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.
Our analysis is done for the following set of parameters As we know, the Higgs prefers to decay into pairs of bb, W W * , ττ , ZZ * and γγ. We restrict our investigation for these channels. Our results are shown in FIG. 6. Perceive that the predictions are in perfect agreement with the values measured by ATLAS [27] and CMS [28].
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 posseses a Higgs-like boson with mass and couplings with the standard particles that recovers its observed properties and respecting the naturalness principle. different from the SM-like one in many aspects with the major one being that itintermediates flavor changing processes involving 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 surppressed even in relation to the loop predictions of the SM [29]. Thus, the significant processes are those involving the Yukawa interactions Although our results are far below the excluded region, we hope this to be probed in the next generation of collider.

VI. 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 Package [31,32]. We 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.

VII. 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.