Extending the Higgs sector: an extra singlet

An extension of the Standard Model with an additional Higgs singlet is analyzed. Bounds on singlet admixture for the 125 GeV h boson from electroweak radiative corrections and data on h production and decays are obtained. The possibility of double h production enhancement at 14 TeV LHC due to a heavy Higgs contribution is considered.


Introduction
After the discovery of the Higgs boson [1,2], all fundamental particles of the Standard Model (SM) have finally been found, and now even passionate adepts of the SM should look for physics beyond it. The pattern of particles we have is rather asymmetric: there are 12 vector bosons, many leptons and quarks with spin 1/2, and only one scalar particle h with mass 125 GeV. Of course, there is only one particle with spin 2 as well, a graviton. However, unlike the spin 2 case, there are no fundamental principles according to which there should exist only one fundamental scalar particle. That is why it is quite probable that there are other still undiscovered fundamental scalar particles in Nature. The purpose of the present paper is to consider the simplest extension of the SM by adding one real scalar field to it. Such an extension of the SM has attracted considerable attention: relevant references can be found in recent papers [3][4][5][6]. An extra singlet can provide the first order electroweak phase transition needed for electroweak baryogenesis. It can act as a particle which connects SM particles to Dark Matter. Not going into these (very interesting) applications (see [7][8][9][10][11][12][13][14][15][16][17][18][19][20]), we will study the degree of enhancement of double Higgs production at LHC due to an extra singlet. To do this we should analyze the a e-mail: sgodunov@itep.ru b e-mail: vysotsky@itep.ru c e-mail: zhemchugov@itep.ru bounds on the mass of the additional scalar particle and its mixing with the isodoublet state.
An enhancement of hh production occurs due to the mixing of the SM isodoublet with an additional scalar field which is proportional to the vacuum expectation value (vev) of this field. Thus the isosinglet is singled out: its vev does not violate the custodial symmetry and can be large. For higher representations special care is needed; see Ref. [21] where an introduction of isotriplet(s) in the SM is discussed.
The paper is organized as follows: in Sect. 2 we describe the model and find the physical states. In Sect. 3 we get bounds on the model parameters of the scalar sector from the experimental data on h production and decays and from precision measurements of Z -and W -boson parameters and the t-quark and h masses. In Sect. 4 we discuss double h production at LHC Run 2.

The model
Adding to the SM a real field X , we take the scalar fields potential in the following form: where Φ is an isodoublet. 1 Terms proportional to X 3 , X 4 , and Φ † Φ X 2 are omitted despite that they are allowed by the demand of renormalizability: we always may assume that they are multiplied by small coupling constants. 2 Two combinations of the parameters entering (1) are known exper- Since we consider the possibility of double h production enhancement due to the H contribution, H → hh, we take m H > m h . Substituting in (1) at the minimum of the potential we get Footnote 2 continued always can choose the signs of the coefficients in front of the omitted terms so that the stability of the potential will not be violated. This is possible due to the freedom in the choice of the values of the bare coefficients.
so μ is negative. For the mass matrix using (4) we get where V φχ ≡ ∂ 2 V ∂φ∂χ , . . . The eigenvalues of (5) determine the masses of the scalar particles: where "−" corresponds to m h and "+" to m H . The eigenfunctions are determined by the mixing angle α: Equation (7) determines μ and λ for the given mixing angle α, while Eq. (6) determines m X for given α as well. Finally, Eq. (4) determines the values of m Φ and v X . Figure 1 demonstrates the dependencies just described for m H = 300 GeV.

ATLAS and CMS Collaborations had detected h production and decays in the reactions
where f i , i = 1, 2, . . . , 5 designate the so-called "Big five" final state channels: W W * , Z Z * , γ γ , ττ , bb. The cross sections of reactions (8) are equal to the Higgs production cross section times the branching ratio of the corresponding decay channel. The quantities μ i are introduced according to the following definition: According to the ATLAS and CMS results, all μ i are compatible with one within experimental and theoretical accuracy. It means that no New Physics are up to now observed in h production and decays.
In the model with an extra isosinglet, production and decay probabilities of h equal that in the SM multiplied by a factor cos 2 α, which is why we have and existing bounds on μ i are translated into bounds on the mixing angle α. Taking into account all measured production and decay channels, for the average values experimentalists obtain [23,24] ATLAS: μ = 1.30 +0.18 −0.17 , CMS: μ = 1.00 +0.14 −0.13 = ±0.09(stat.) +0.08 −0.07 (theor.) ± 0.07(syst.). (12) Let us stress that the theoretical uncertainty in the calculation of pp → h production cross section at LHC does not allow one to reduce substantially the uncertainty in the value of μ. The bounds from electroweak precision observables (EWPO) are not affected by this particular uncertainty. We fit experimental data with the help of LEPTOP program [25][26][27][28][29] using m h = 125.14 GeV. The result of the SM fit which accounts for the h mass measurement is shown in Table 1. The quality of the fit is characterized by the χ 2 value Higgs boson contributions to electroweak observables at one loop are described in LEPTOP by the functions In the case of an extra singlet the following substitution should be performed: The same substitution should be made for the functions δ 4 V i (t, h), t = m 2 t /m 2 Z , which describe two loops radiative corrections enhanced as m 4 t . In two loops a quadratic dependence on the Higgs mass appears which is described by the functions δ 5 V i . The calculation of these corrections in the case of an extra singlet Higgs is not easy. However, even for a 1000 GeV Higgs these corrections are small. We checked that corrections to the values of sin α due to δ 5 V i terms in Fig. 2 are less than 10 −3 . The bounds from EWPO on the singlet model parameters are presented in Fig. 2a. A χ 2 minimum is reached at sin α = 0, m H = 150 GeV, which is the minimum value allowed for m H in the fit. Experimental data are avoiding a heavy Higgs. The value of χ 2 at the minimum coincides with the SM result (13). Lines of constant χ 2 correspond to χ 2 = 1, 4, 9, . . . The probabilities that (sin α, m H ) values are below these lines are 39, 86, 98.9 %, … 3 Bounds accounting for both EWPO and direct h production data (11), (12), are shown in Fig. 2b. We see that for a heavy H the bounds from EWPO dominate, while for a light H the measurement of μ is more important.

h, H, and hh production at LHC
The main purpose of this section is to find what enhancement of double Higgs production cross section is possible with enlarged Higgs sector. Let us recall that in the SM the double h production cross section is very small. According to the recent result [30], at √ s = 14 TeV σ NNLO ( pp → hh) = 40 fb with a 10 ÷ 15 % accuracy. We will demonstrate that an enlarged Higgs sector allows one to strongly enhance double h production.
The cross section of H production at LHC equals that for the SM Higgs production (for (m h ) SM = m H ) multiplied by sin 2 α. We take the cross section of the SM Higgs production at NNLO from Table 3 of [31]. In order to obtain the cross section of resonant hh production in H decays we should multiply the cross section of H production by Br(H → hh).
Let us consider H decays. Decays to hh, W + W − , Z Z, and tt dominate. For the H hh coupling we obtain thus 3 Let us note that if a subset of the experimental data from Table 1 is fitted, then the allowed domains of the (sin α, m H ) values will be larger than those presented in Fig. 2a. Here we disagree with the statement made in [4] that the fit of only one observable (m W ) allows one to set the strongest constraint on (sin α, m H ). thus The dependences of the widths and branching ratios of H decays on the mixing angle α for m H = 300 GeV are shown in Fig. 3.
For the cross section of the reaction pp → H → hh we have Contour plot of R is presented in Fig. 5. Let us note that R does not depend on √ s. See also σ ( pp → H )·Br(H → Z Z) for different values of sin α in Fig. 6.

Conclusions
In the models with extended Higgs sector strong resonant enhancement of double Higgs production is possible, which makes the search of the pp → hh reaction at Run 2 LHC especially interesting. According to Fig. 4 the cross section of the pp → H → hh reaction can be as large as 0.5 pb, ten times larger than the SM value.
The search for a H boson can go in the same way as for the heavy SM boson h. The probability of H observation is diminished compared to that of h because of (a) suppression of H production cross section by the factor sin 2 α ≤ 0.2; and (b) suppression of Br(H → Z Z) because of an additional H → hh decay mode. Taking these two factors into account, we get a suppression of about a factor 10 of the pp → H → Z Z process probability compared to that for the SM Higgs boson (see Fig. 5).
Results for the search of the Higgs-like boson in Z Z decay mode by CMS can be found in [32], Figure 18, [33], Figure  7, and by ATLAS in [34], Figure 12. Comparing it with our Fig. 5, we observe that the experimental data start to be sensitive to the singlet model expectation for maximally allowed values of the mixing angle α.