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

TL;DR

This paper investigates whether the lightest Higgs boson in the supersymmetric minimal 3-3-1 model can replicate the observed 125 GeV Higgs, analyzing its mass, couplings, and decay signatures in comparison with experimental data.

Contribution

It develops the Higgs potential in the supersymmetric minimal 3-3-1 model and evaluates the Higgs-like boson's properties, including mass, couplings, and flavor-changing processes, aligning them with experimental results.

Findings

Higgs mass calculated up to one loop matches 125 GeV.

Higgs couplings with standard particles are consistent with CMS and ATLAS data.

Decay $t ightarrow h+c$ has a calculable branching ratio within bounds.

Abstract

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.

Equations

$$M_h^2\le 4M_Z^2{\sin }^2({\theta }_w)(1+{\cos }^2(\beta ){)}^2,$$

$$\displaystyle\hat{L}_{a}=\left(\begin{array}[]{c}\hat{\nu}_{a}\\ \hat{e}_{a}\\ \hat{e}_{a}^{c}\end{array}\right)_{L}\sim(1\,,\,3\,,\,0)\,,$$

$$\displaystyle\hat{Q}_{\alpha_{L}}=\left(\begin{array}[]{c}\hat{d}_{\alpha}\\ \hat{u}_{\alpha}\\ \hat{j}^{\prime}_{\alpha}\end{array}\right)_{L}\sim(3,3^{*},-\frac{1}{3}),$$

$${\hat{u}}_{{\alpha }_L}^c\sim (3^{\ast },1,-\frac{2}{3}),{\hat{d}}_{{\alpha }_L}^c\sim (3^{\ast },1,1/3),{\hat{j}}_{{\alpha }_L}^{\prime c}\sim (3^{\ast },1,\frac{4}{3}),$$

$$\displaystyle\hat{Q}_{3L}=\left(\begin{array}[]{c}\hat{u}_{3}\\ \hat{d}_{3}\\ \hat{j}^{\prime}_{3}\end{array}\right)_{L}\sim(3,3,\frac{2}{3}),$$

$${\hat{u}}_{3_L}^c\sim (3^{\ast },1,-\frac{2}{3}),{\hat{d}}_{3_L}^c\sim (3^{\ast },1,1/3),{\hat{j}}_{3L}^c\sim (3^{\ast },1,-\frac{5}{3}),$$

$$\displaystyle\hat{\eta}=\left(\begin{array}[]{c}\hat{\eta}\\ \hat{\eta}_{1}^{-}\\ \hat{\eta}_{2}^{+}\\ \end{array}\right),\hat{\chi}=\left(\begin{array}[]{c}\hat{\chi}^{-}\\ \hat{\chi}^{--}\\ \hat{\chi}\\ \end{array}\right),\hat{\rho}=\left(\begin{array}[]{c}\hat{\rho}^{+}\\ \hat{\rho}\\ \hat{\rho}^{++}\\ \end{array}\right),$$

$$\displaystyle\hat{\eta}^{\prime}=\left(\begin{array}[]{c}\hat{\eta}^{\prime}\\ \hat{\eta}_{1}^{\prime+}\\ \hat{\eta}_{2}^{\prime-}\\ \end{array}\right),\hat{\chi^{\prime}}=\left(\begin{array}[]{c}\hat{\chi}^{\prime+}\\ \hat{\chi}^{\prime++}\\ \hat{\chi}^{\prime}\\ \end{array}\right),\hat{\rho}^{\prime}=\left(\begin{array}[]{c}\hat{\rho}^{\prime-}\!\!\!\\ \hat{\rho}^{\prime}\\ \hat{\rho}^{\prime--}\!\!\!\\ \end{array}\right),$$

$$⟨\eta ⟩=\frac{v_{\eta }}{\sqrt{2}},⟨{\eta }^{\prime }⟩=\frac{v_{{\eta }^{\prime }}}{\sqrt{2}},⟨\rho ⟩=\frac{v_{\rho }}{\sqrt{2}},⟨{\rho }^{\prime }⟩=\frac{v_{{\rho }^{\prime }}}{\sqrt{2}},⟨\chi ⟩=\frac{v_{\chi }}{\sqrt{2}},⟨{\chi }^{\prime }⟩=\frac{v_{{\chi }^{\prime }}}{\sqrt{2}}.$$

$$SU(3{)}_C\otimes SU(3{)}_L\otimes U(1{)}_X\stackrel{v_{\chi },v_{{\chi }^{\prime }}}{⟹}SU(3{)}_C\otimes SU(2{)}_L\otimes U(1{)}_Y\stackrel{v_{\eta },v_{{\eta }^{\prime }},v_{\rho },v_{{\rho }^{\prime }}}{⟹}SU(3{)}_C\otimes U(1{)}_{QED}.$$

$$M_{W^{\pm }}$$

$$M_{V^{\pm }}$$

$$M_{U^{\pm \pm }}$$

$$M_Z$$

$$M_{Z^{\prime }}$$

$$W$$

$${\mathcal{L}}_{soft}$$

$$b_{\rho ,\eta ,\chi }$$

$$v_{331}$$

$$\tan {\beta }_{\eta }$$

$$X_t$$

$$V_F$$

$$V_D$$

$$V_{soft}$$

$${\rho }^0,{\rho }^{\prime 0},{\eta }^0,{\eta }^{\prime 0}{\chi }^0,{\chi }^{\prime 0}\to \frac{1}{\sqrt{2}}(v_{\rho ,{\rho }^{\prime },\eta ,{\eta }^{\prime }\chi ,{\chi }^{\prime }}+R_{\rho ,{\rho }^{\prime },\eta ,{\eta }^{\prime }\chi ,{\chi }^{\prime }}+i{I}_{\rho ,{\rho }^{\prime },\eta ,{\eta }^{\prime }\chi ,{\chi }^{\prime }}),$$

$$\text{centering}-0.03\le f_1\le 0.2\text{@add@centering}$$

$$-2000\mathrm{GeV}\le k_1\le 2000\mathrm{GeV}$$

$$-800\mathrm{GeV}\le {\mu }_{\eta }\le -500\mathrm{GeV},-1800\mathrm{GeV}\le$$

$$2\le \tan \beta \le 4,8\le (\tan {\beta }_{\eta }$$

$$v_{\rho }^2+v_{\eta }^2+v_{{\eta }^{\prime }}^2+v_{{\rho }^{\prime }}^2$$

$$\tan (\beta )=2.7,\tan ({\beta }_{\eta })=10,\tan ({\beta }_{\rho })=10,\tan ({\beta }_{\chi })=1.2,v_{\eta }=22.3\mathrm{GeV},v_{{\eta }^{\prime }}=223.2\mathrm{GeV},$$

$$v_{\rho }=8.1\mathrm{GeV},v_{{\rho }^{\prime }}=81\mathrm{GeV},v_{\chi }=1763.6\mathrm{GeV},v_{{\chi }^{\prime }}=2116.4\mathrm{GeV},f_1=0.01,f_1^{\prime }=0.6,$$

$$k_1=-500\mathrm{GeV},k_2=-600\mathrm{GeV},b_{\eta }=-55000{\mathrm{GeV}}^2,b_{\rho }=-1.7\times 10^6{\mathrm{GeV}}^2,$$

$$b_{\chi }=-1.4\times 10^6{\mathrm{GeV}}^2,{\mu }_{\rho }=-1700\mathrm{GeV},{\mu }_{\eta }=-550\mathrm{GeV},{\mu }_{\chi }=-1400\mathrm{GeV}.$$

$$h=-0.93R_{{\eta }^{\prime }}-0.34R_{{\rho }^{\prime }}-0.09R_{\eta }-0.03R_{\rho }-0.05R_{{\chi }^{\prime }}-0.06R_{\chi }.$$

$$\text{centering}{f}_1=0.001\text{@add@centering}$$

$$k_1=-500\mathrm{GeV}$$

$${\mu }_{\eta }=-550\mathrm{GeV},{\mu }_{\rho }=$$

$$\tan \beta =2,\tan {\beta }_{\eta }=$$

$$A_0,B_0$$

$$f_1^{\prime }=0.6,-0.03\le f_1\le 0.2,$$

Full text

I Introduction

The discovery of a scalar state by the ATLAS and CMS at CERN Aad et al. (2012); Chatrchyan et al. (2012) 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\pm 0.21(stat.)\pm 0.11(syst.)$ GeV, but also know that its signal strengths are found consistent with the standard model (SM) predictionsAad et al. (2016). 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 principleCarena and Haber (2003). 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}\times SU(3)_{L}\times U(1)_{N}$ (331) model. The motivations to study 331 models are various as, for example, explanation of family replicationPisano and Pleitez (1992)Frampton (1992), electric charge quantizationde Sousa Pires and Ravinez (1998)de Sousa Pires (1999), strong CP-problemPal (1995)Dias and Pleitez (2004)Dias et al. (2003), incorporation of inflationFerreira et al. (2017) 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 Duong and Ma (1993). For subsequent works, see Montero et al. (2002)Huong et al. (2013)Ferreira et al. (2013). As the Higgs sector of the minimal 331 modelPisano and Pleitez (1992)Frampton (1992) requires at least three Higgs triplets to generate correctly the masses of all massive particles and be phenomenologically viableMontero et al. (2002)Queiroz et al. (2010)Dong and Si (2014), 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 $\chi$ and $\chi^{\prime}$ decouple from the other four $\eta$, $\rho$, $\eta^{\prime}$ and $\rho^{\prime}$, and that $\eta$ and $\rho$ 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

[TABLE]

where $M_{Z}$ is the mass of the standard neutral gauge boson, $\theta_{w}$ is the electroweak mixing angle and the angle $\beta$ is such that $\tan\beta=\frac{v_{\eta^{\prime}}}{v_{\rho^{\prime}}}$. According to this expression the tree level contribution to the lightest Higgs may attain 174 GeV for $\cos(\beta)=1$.This is so because the superpotential naturally contains cubic invariant terms that act like in the NMSSMDuong and Ma (1993). 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\rightarrow h+c$ and $t\rightarrow 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 symmetry

[TABLE]

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 by

[TABLE]

where $\alpha=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 superfields

[TABLE]

where $\hat{\eta}\sim(1,3,0)\,,\hat{\chi}\sim(1,3,-1)\,,\,\,\,$ $\hat{\rho}\sim(1,3,1)$, and

[TABLE]

where $\hat{\eta}^{\prime}\sim(1,3^{*},0)\,,\hat{\chi}^{\prime}\sim(1,3^{*},1)\,,\,\,\,$ $\hat{\rho}^{\prime}\sim(1,3^{*},-1)$.

We assume that the neutral scalars $\eta$, $\eta^{\prime}$, $\rho$, $\rho^{\prime}$, $\chi$ and $\chi^{\prime}$ develop nonzero VEV according to

[TABLE]

These VEVs lead to the following gauge symmetry breaking pattern

[TABLE]

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, $\gamma$, $Z^{0}$ and $W^{\pm}$, one new neutral massive gauge bosons $Z^{\prime}$, two doubly charged gauge bosons $U^{\pm\pm}$, and two simply charged gauge bosons $V^{\pm}$. The charged gauge bosons gain the following mass expressions

[TABLE]

while the neutral gauge bosons have masses

[TABLE]

where $t=g_{N}/g$, $v_{\rho}^{2}+v^{2}_{\rho^{\prime}}+v_{\eta}^{2}+v_{\eta^{\prime}}^{2}=v^{2}_{ew}$ and $v_{\chi}^{2}+v_{\chi^{\prime}}^{2}\equiv v^{2}_{331}$ 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

[TABLE]

where $\alpha,\beta=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 $\Lambda=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 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

[TABLE]

where $\alpha,\beta=1,2$ and $i=1,2,3$. $\lambda_{c}^{a}$ are the gluinos, $\lambda^{a}$ are gauginos associated to $SU(3)_{L}$ (in both cases $a=1,...,8$ is the gauge group index) and $\lambda$ is the gaugino associated to $U(1)_{N}$. The scalar supersymmetric partners of fermion fields, $f$, are denoted by $\tilde{f}$, while the remaining fields are self-evident. For simplicity, we will take $m_{L}^{2}=m_{Q_{\alpha}}^{2}=m_{Q_{3}}^{2}=m_{SUSY_{L}}^{2}$, $m_{u_{i}}^{2}=m_{d_{i}}^{2}=m_{j\beta}^{2}=m_{SUSY_{R}}^{2}$ and $m_{J}^{2}=10\times m_{SUSY_{R}}^{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

[TABLE]

where $\kappa$’s are the corresponding Yukawa couplings, $M_{s}$ is the SUSY scale, $m_{\tilde{t}}$ 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_{SUSY_{L}}$ and $m_{SUSY_{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_{soft}$. These contributions are, respectively

[TABLE]

where $t_{A}$ are the Gell-Mann matrices.

The fields are assumed to be shifted in the usual way

[TABLE]

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 Higgs in what concerns its mass, branching ratios and signal strengths. For this we first obtain the $6\times 6$ mass matrix, $M^{2}_{R_{H}}$, associated to the CP-even scalars $R_{\eta}$, $R_{\eta^{\prime}}$, $R_{\rho}$, $R_{\rho^{\prime}}$, $R_{\chi}$ and $R_{\chi^{\prime}}$. 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}_{H_{D}}$ relates to the $M^{2}_{R_{H}}$ through the relation $M^{2}_{H_{D}}=U^{\dagger}_{R}M^{2}_{R_{H}}U_{R}$. The physical eigenstates relates to the symmetrical ones by $H=U_{R}R$ where $H=\left(h_{1},h_{2},h_{3},h_{4},h_{5},h_{6}\right)^{T}$ and $R=\left(R_{\eta},R_{\eta^{\prime}},R_{\rho},R_{\rho^{\prime}},R_{\chi},R_{\chi^{\prime}}\right)^{T}$. 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 Staub (2008, 2015); Vicente (2015), Spheno Porod (2003); Porod and Staub (2012) and SSP Staub et al. (2012). 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.

III.1 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\times 6$ mass matrix, agrees with the prediction in Eq. (1).

The results displayed in FIG. 1 is for the range of parameters

[TABLE]

Observe that the tree level contribution may attain $115$ GeV for $\tan(\beta)=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(\beta)$ such that the smaller $\tan(\beta)$, 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(\beta)\geq 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

[TABLE]

In this case, the eigenvector corresponding to the 125 GeV Higgs is

[TABLE]

See that the eigenvector $h$ is composed mainly by $R_{\eta^{\prime}}$ and $R_{\rho^{\prime}}$. For any other choice of the set of values of the parameters, the eigenvector will keep being dominantly a composition of $R_{\eta^{\prime}}$ and $R_{\rho^{\prime}}$.

III.2 One Loop Level Contribution

In the MSSM a 125 GeV Higgs demands robust loop corrections which requires stop heavy enough such that threaten the naturalness principle. In the minimal SUSY331 model things are a little different once tree level contribution to the Higgs-like boson mass may surpass $100$ GeV.

In what follows we calculate numerically, by employing the package Spheno, the Higgs-like boson mass but taking into account one loop corrections. Our results are displayed in graphics showing the behavior of our Higgs-like boson mass with the parameters $M_{s}$, $X_{t}$ and $v_{331}$. For the other parameters, we take the values

[TABLE]

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 may obtain a 125 GeV Higgs in the minimal SUSY331 model even for $X_{t}=0$. We highlight this result because it is not allowed in the MSSM.

We also obtained the $2\times 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 we show the behavior of $M_{h}$ with $M_{s}$ for the specific case $X_{t}=0$, and in FIG. 5 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 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 Andersen et al. (2013). 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

[TABLE]

As we know, the Higgs prefers to decay into pairs of $b\bar{b}$, $WW^{*}$, $\tau\bar{\tau}$, $ZZ^{*}$ and $\gamma\gamma$. 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 ATLASAad et al. (2016) and CMSCollaboration (2014). 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.

V 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 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_{\eta^{\prime}}$ and $R_{\rho^{\prime}}$, 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\rightarrow bs\,,\,bd$ are very surppressed even in relation to the loop predictions of the SMBenitez-Guzm n et al. (2015). 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 Khachatryan et al. (2017). The behavior of the branching ratio of these decays as function of $v_{331}$ is presented in FIG. 7. Our calculations were done for the following region of the parameter space

[TABLE]

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 Camargo-Molina et al. (2014, 2013). 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(\beta)$ versus $v_{331}$ displayed in FIG. 8 .

[TABLE]

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\rightarrow h+u$ and $t\rightarrow 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.

Acknowledgements.

We thank Felipe Ferreira for useful discussions. This work was supported by Conselho Nacional de Pesquisa e Desenvolvimento Científico- CNPq (C.A.S.P, P.S.R.S.). CS is supported by CAPES/PDSE Process 88881.134759/2016-01.

Appendix A

[TABLE]