Pair production of the Elementary Goldstone Higgs boson at the LHC

The Elementary Goldstone Higgs (EGH) model is a perturbative extension of the standard model (SM), which identifies the EGH boson as the observed Higgs boson. In this paper, we study pair production of the EGH boson via gluon fusion at the LHC and find that the resonant contribution of the heavy scalar is very small and the SM-like triangle diagram contribution is strongly suppressed. The total production cross section mainly comes from the box diagram contribution and its value can be significantly enhanced with respect to the SM prediction.


Introduction
The discovery of a Higgs boson at the LHC [1,2], which in principle seems to be the Higgs boson of the standard model (SM), opens a new window to physics related to the Higgs boson within and beyond the SM. The possibility of detailed and careful study of Higgs physics in the TeV range is higher than past years with the upgrade of the CMS and ATLAS detectors [3][4][5][6]. For the current and future high-energy experiments, such as the LHC Run II and e + e − colliders, one of their main goals is to measure the processes which give information on the Higgs boson couplings to fermions, gauge bosons and its self-coupling and to compare with the most accurate SM predictions. Any deviation from the SM predictions could be valuable information about physics beyond the SM.
To address the SM gauge hierarchy problem, several alternative paradigms have been put forward. The Elementary Goldstone Higgs (EGH) scenario [7,8] is one possible paradigm, which is based on an elementary scalar sector with a global symmetry larger than that in the SM. This symmetry is explicitly broken by the couplings with the electroweak (EW) gauge currents and SM Yukawa interactions. Under radiative corrections, this symmetry breaking will align the vacuum with respect to the EW symmetry. The observed Higgs boson is a fundamental pseudo Nambu-Goldstone (PNG) boson of the global symmetry breaking. It can obtain a light mass through radiative corrections which could cause the symmetry breaking and explain the origin of the known fermion masses.
In the EGH model [7,8], the elementary Higgs sector of the SM is enhanced to an SU (4) symmetry that breaks spontaneously to Sp(4). The embedding of the EW gauge sector is parameterized by an angle θ. Its value can be fixed by minimizing the quantum corrected effective potential of this model in the presence of the EW and top corrections, which has been dynamically determined to be centered around θ≈0.018 [8].
The EGH model is a perturbative extension of the SM, in which the Higgs boson is a fundamental particle, like in the SM, but the mechanism of symmetry breaking is completely different. Thus, it might produce rich new phenomenology at the LHC Run II or future high energy collider experiments, which has not been fully studied in the literature. The relations of the EGH idea with unification scenario, the relaxation leptogenesis mechanism, and supersymmetry have been studied in Refs. [9,10] and [11], respectively. SM Higgs inflation has also been discussed in the context of the EGH model [12].
Pair production of the Higgs boson is well known for its sensitivity to the tri-Higgs coupling, providing a way to test the structure of Higgs potential and further EW symmetry breaking mechanism. In the SM, di-Higgs production at the LHC, mainly from gluon fusion, arises from both triangle and box loop contributions, which interfere destructively, causing a suppression of the total production rate from the naive estimate [13,14]. Thus, sizable production of di-Higgs directly implies a new physics signature [15]. The main goal of this paper is to examine pair production of the EGH boson at the LHC. We find that, although its trilinear coupling with respect to the SM one is strongly suppressed, the box loop contributions, the resonant contributions of the heavy scalar and the constructive interference effects can make up this suppression and enhance the production rate of EGH boson to a sizeable level.
This paper is structured as follows. In Section 2, we summarize the main features of the EGH model and show the relevant couplings. The possible decay channels of the EGH boson H 1 , which is identified as the observed Higgs boson, and the heavy scalar H 2 , are discussed in Section 3. Section 4 studies pair production of the EGH boson at the LHC, and our conclusions are given in Section 5.

Main features of the EGH model
For the EGH model [7,8], the Higgs sector is embedded into a SU (4) → Sp(4) pattern of chiral symmetry breaking. The SM Higgs doublet is part of the SU (4)/Sp(4) coset, while the EW symmetry, SU (2) L × U (1) Y , is embedded in SU (4). The SM Higgs boson is identified with one of the Goldstone bosons which acquires mass via a slight vacuum misalignment mechanism induced by quantum corrections. The most general vacuum E θ of this model can be expressed as where 0 θ π/2 and the two independent vacua E B and E H are defined as Here σ 2 is the second Pauli matrix. The Higgs sector of this model strictly depends on the choice of the vacuum E θ . The alignment angle θ is completely determined by the radiative corrections and the requirement that the EGH model reproduces the phenomenological success of the SM, which prefers small values of θ.
The above vacuum structure can be realised by introducing the scalar matrix M , which is the two-index antisymmetric irrep ∼6 of SU (4) where X i θ (i = 1,...,5) are the broken generators associated with the breaking of SU (4) to Sp(4), reported in Appendix A of Ref. [8].
In order to embed the EW gauge sector of the SM into the SU (4) group, the EGH model gauges the SU (2) L ×U (1) Y part of the chiral symmetry group SU (2) L ×SU (2) R ⊂ SU (4). In this case, the scalars are minimally coupled to the EW gauge bosons via the covariant derivative of M R . The kinetic and EW gauge interaction Lagrangian of the scalar sector can be written as which explicitly breaks the global SU (4) symmetry. In the EGH model [7,8], the EGH boson is one of the two linear combinations of the PNG bosons σ and Π 4 at low energy, that is where H 1 and H 2 are mass eigenstates, α is the mixing angle and its value taken as 0 < α < π/2. The lightest of H 1 and H 2 is identified as the observed Higgs boson with mass m h =125.09±0.24 GeV [16].
The renormalizability of the EGH model together with the perturbative corrections determine dynamically the direction of the vacuum E θ . By investigating the available parameter space of the scalar sector, Ref. [8] has shown that the preferred value of the vacuum alignment angle θ is θ = 0.018 +0.004 −0.003 , corresponding to the SU (4) spontaneous symmetry breaking scale of f =13.9 +2.9 −2.1 TeV via the phenomenological constant f sinθ = ν = 246 GeV and the mixing angle α = 1.57. This means that the EGH boson H 1 is taken as the observed Higgs boson and is expressed by h, like the SM Higgs boson. It is mainly comprised of the PNG boson Π 4 with a tiny admixture of σ, while the heavier scalar H 2 is mainly made up of the PNG boson σ, which is taken as H.
The first three of the five PNG bosons Π i (i=1,...,5) become the longitudinal components of the EW gauge bosons W and Z, while the fourth is used to constitute the EGH boson h. Reference [7] has shown that Π 5 is a stable massive particle and provides a viable dark matter candidate. For M Π 5 = M DM m h , the EGH model is compatible with the experimental constraints. In our following numerical calculation, we will take M DM =m h .
The SM normalised coupling strength of the scalars predicted by the EGH model can be written as [7,8]: Here f f denote all of the fermion pairs, V V =WW and ZZ, and λ SM hhh = 3m 2 h /ν is the SM trilinear self-coupling constant of the Higgs boson. In the following section, we will use these relations to consider the possible decay channels of the scalar particles h and H, and focus our attention on the branching ratios of the heavy scalar H.

Decays of scalars h and H
From the discussions given in Section 2, we can see that, except the decay mode Π 5 Π 5 , the decay modes of the EGH boson h are same as those of the SM Higgs boson and its partial decay widths are universally shifted from the SM predictions by a factor sin 2 (α+θ). Since the decay channel h→Π 5 Π 5 is kinematically prohibited, the values of the branching ratios BR(h→XX) are also same as those of the SM Higgs boson.
For a specific production process and decay mode i→h→f , the Higgs signal strength is defined as Here σ i (i = ggF, V BF, Wh, Zh, tth) and BR f (f = ZZ, WW, γγ, ττ, bb, µµ) are respectively the production cross section for i → h and the branching ratio for the decay process h → f . The subscript "SM" refers to their respective SM predictions. The values of µ f i are the same for all production processes i and decay channels f in the EGH model, and µ f i =µ=sin 2 (α+θ). A fit to the combined ATLAS and CMS data at the center-of-mass (c.m.) energies √ s=7 and 8 TeV give the best fit value of the Higgs signal strength µ as µ = 1.09 +0.11 −0.1 [17,18]. The preferred values θ = 0.018 +0.004 −0.003 and α = 1.57 given by Ref. [8] satisfy this experimental constraint.
The heavy scalar H can decay to the SM gauge bosons and fermions with partial widths of where Γ SM (H→XX) is the total decay width of the SM Higgs boson into the XX final states evaluated at the H mass M H . At tree level, the partial decay widths for the processes H→hh and Π 5 Π 5 can be written as In the above equation we have taken M Π 5 ≈ m h . Then the total decay width of the heavy scalar H is given by Considering the heavy scalar H mainly from the PNG boson σ, we take M H ≈ M σ and assume its values in the range of 1 TeV ∼ 3 TeV, using the Mad-graph5/aMC@NLO program [19] to give our numerical results. In Fig. 1 we plot the branching ratio BR(H→hh) as a function of the mass parameter M H for α=1.57 and θ = 0.018 +0.004 −0.003 , where the breaking scale f is taken as f = ν/sinθ. In most of the parameter space, the decay channel H→hh is one of the dominant ones of the heavy Higgs boson and its branching ratio value is about 20%.  For α=1.57 and θ =0.018 +0.004 −0.003 , the allowed range of the ratio Γ H /M H is shown as a function of M H in Fig. 2. In almost all the parameter space, the value of Γ H /M H is smaller than 0.5%. Thus, we can safely use the narrow width approximation method to calculate the production cross section of the process gg→H→hh. Pair production of the EGH boson h at the LHC is mainly induced by two sources: one comes from the diagrams depicted in Fig. 3(a) and 3(b), which are similarly to the SM di-Higgs production, while the other comes from resonant process through the heavy scalar H decay H → hh, as shown in Fig. 3(c). The particles in the loops are the heavy SM fermions, such as top and bottom quarks.
In the EGH model, the trilinear coupling of the EGH boson and its couplings to the SM fermions are modified from the relevant SM couplings. Thus, compared with the SM prediction, the production cross section induced by Figs. 3(a) and 3(b) is changed as where σ T , σ B and α I respectively represent the cross sections only from the triangle and box graphs, and the interference angle α I . K F h = µ h = 1 for the SM. For the SM, Ref. [20] has calculated the values of σ T and σ B for various c.m. energy √ s and shown that cosα I is almost independent of the c.m. energy √ s, the scale and the parton distribution functions (PDFs), and cosα I = −0.898 for √ s=13 TeV. The second source contributing to EGH pair production at the LHC is resonant contribution from Fig. 3(c). If a new scalar is sufficiently heavy and can decay on-shell into two SM-like Higgs bosons, it has been shown that the new scalar might significantly enhance the di-Higgs production rate over the SM prediction via its resonant production by gluon-gluon fusion [21][22][23][24][25][26][29][30][31][32][33]. Using the narrow width approximation, the resonant production cross section can be approximately written as where σ(gg → h) is the cross section of single production for the SM Higgs boson via the gluon-gluon fusion process evaluated at the H mass M H .
The interference effects between the two kinds of contributions, the SM-like and heavy scalar contributions, might be significantly large [32,33]. So, in our numerical calculation, we will include these effects. Then the total cross section of pair production of the EGH boson at the LHC is written as where σ 12 represents the cross section from the interference contributions. Using the PDFs of CT14 [34] and the NNLO κ-factor as κ=2.30 for √ s=13 TeV [35], we plot in Fig. 4 the cross section σ 2 contributed by the heavy scalar H resonant contributions as a function of the mass parameter M H for α=1.57 and θ=0.018 +0.004 −0.003 . From Fig. 4, the production cross section is very small and σ 2 0.01 fb for M H 1 TeV, which cannot be de-tected at the LHC in the near future. This is unlike the resonant enhancement arising from a new scalar, which is because of the suppression factor cos 2 (α+θ).  In the EGH model, the trilinear self-coupling of the SM-like Higgs boson with respect to the SM one is strongly suppressed, while its couplings with the SM fermions are almost the same as those for the SM Higgs boson. Thus, the cross section of pair production of the EGH boson at the LHC is dominated by the box contributions [ Fig. 3(b)] and its destructive interference with the triangle diagram [ Fig. 3(a)] is suppressed. Furthermore, the interference contribution between Fig. 3(b) and Fig. 3(c) is positive. So, in the context of the EGH model, the production cross section is larger than that in the SM. Our numerical results show that it is indeed this case. Fig. 5 shows the ratio R = σ EGH /σ SM as a function of the mass parameter M H , where the value of the SM di-Higgs production cross section σ SM is taken as 37.91 fb calculated at NNLL level for m h =125.09 GeV and the c.m. energy √ s=13 TeV [35]. One can see from this figure that the production rate for pair production of the SM-like Higgs boson in the EGH model is significantly larger than the SM prediction.  The triangle and box diagrams contributing to di-Higgs production have different phase space dependencies and distinct kinematic features. The triangle contribution would peak around √ŝ =m h and become subleading at larger invariant masses. A modified triple Higgs self-coupling of the SM-like Higgs boson will mostly reveal itself at low invariant masses, while shifted Yukawa couplings will typically become more apparent in the larger invariant mass region. The small change in the triple Higgs self-coupling will mainly affect the low invariant mass region which will be removed by the cuts. However, when the triple Higgs self-coupling is nearly turned off, the high invariant mass region signal event number will not increase as much as the total cross section [36]. In the EGH model, the trilinear self-coupling of the SM-like Higgs boson is strongly suppressed, while its couplings with the SM fermions are almost the same as those for the SM Higgs boson. Furthermore, the production cross section contributed by the heavy scalar H is very small for M H 1 TeV. Thus, in the context of the EGH model, the production rates of the final states, such as bbγγ and bbττ, can be enhanced with respect to the SM ones. As a reasonable estimation, the value of the ratio R for the signal event number is in the range of 1.6 ∼ 1.7 after the cuts effect is considered.
Currently, the ATLAS and CMS collaborations have been searching for the signal of pair production of the Higgs boson at the LHC and also focusing on resonanceenhanced production mechanisms [37][38][39][40][41][42][43][44]. Upper limits on the production cross sections for some categories of signal final states have been obtained. In the near future, the LHC will give more meaningful data, which can be used to test the SM and further probe new physics beyond the SM. Our work will help to examine the EGH model via pair production of the SM-like Higgs boson at the LHC.

Conclusions
The EGH model is a perturbative extension of the SM featuring an EGH boson and dark matter particles, which is theoretically well-motivated and phenomenologically viable. In this model, the EGH boson is taken as the observed Higgs boson, which has almost the same couplings as the SM fermions and the EW gauge bosons as those for the SM Higgs boson. Thus, this model can easily satisfy the experimental constraints from the Higgs signal data at the LHC.
The trilinear self-coupling of the EGH boson is suppressed with respect to the SM one, which modifies the cross section for pair production of the EGH boson from the SM prediction at the LHC. In addition, the existence of the heavy scalar H in the EGH model gives an additional contribution to this cross section. So, in this paper, we have studied pair production of the EGH boson via the gluon-gluon fusion process at the LHC.
Our numerical results show that, since the couplings of the heavy scalar H to the SM fermion pairs are strongly suppressed by the factor cos(α+θ) for α = 1.57 and θ = 0.018 +0.004 −0.003 , its resonant contribution to the di-Higgs production cross section is negligible and its value is smaller than 0.1 fb in most of the parameter space, while its interference contribution with the box diagram is positive and cannot be neglected. In the EGH model, the total production cross section of the SM-like Higgs boson pair is larger than the SM prediction. The value of the ratio R=σ EGH /σ SM is in the range of 2.04∼1.91 for α = 1.57, θ = 0.018 +0.004 −0.003 and 1 TeV M H 3 TeV. The EGH model might be probed or ruled out by the LHC via the di-Higgs production process in the near future.