Investigation of CP-even Higgs bosons decays $H\to \mu \tau$ within constraints of $l_a\to l_b\gamma$ in a 3-3-1 model with inverse seesaw neutrinos

In a 3-3-1 model with inverse seesaw neutrinos, we use a simple form of Higgs potential to give four CP-even Higgs bosons ($H\equiv h_1^0,h_2^0,h_3^0,h_4^0$). We investigate $H\to \mu \tau$ decays in the parameter space regions satisfying the experimental limits of $l_a\to l_b\gamma$ with running parameters being the mass of the charged Higgs boson ($m_{H_1^{\pm }}$) and the mixing matrix of the heavy neutrinos ($M_R$). We show that there exist regions of parameter space where all partial widths $\mathrm{\Gamma }(H\to \mu \tau )$ are less than the current experimental limit ($4.1\times {10}^{-6}GeV$). Analyzing the contributing components to $\mathrm{\Gamma }(H\to \mu \tau )$, we also compare the mass of the SM-like Higgs boson with the corresponding ones of the other CP-even Higgs bosons in this model.