Radiative corrections to the Triple Higgs Coupling in the Inert Higgs Doublet Model

We investigate the implication of the recent discovery of a Higgs-like particle in the first phase of the LHC Run 1 on the Inert Higgs Doublet Model (IHDM). The determination of the Higgs couplings to SM particles and its intrinsic properties will get improved during the new LHC Run 2 starting this year. The new LHC Run 2 would also shade some light on the triple Higgs coupling. Such measurement is very important in order to establish the details of the electroweak symmetry breaking mechanism. Given the importance of the Higgs couplings both at the LHC and $e^+e^-$ Linear Collider machines, accurate theoretical predictions are required. We study the radiative corrections to the triple Higgs coupling $hhh$ and to $hZZ$, $hWW$ couplings in the context of the IHDM. By combining several theoretical and experimental constraints on parameter space, we show that extra particles might modify the triple Higgs coupling near threshold regions. Finally, we discuss the effect of these corrections on the double Higgs production signal at the $e^+e^-$ LC and show that they can be rather important.


I. INTRODUCTION
The discovery of a new particle with a mass around 125-126 GeV in the search for the Standard Model (SM) Higgs boson [1] was announced simultaneously by the ATLAS and CMS collaborations in July 2012 [2][3][4][5][6][7]. Since then, more data has been taken and analyzed at the LHC. One of the primary goals of the Higgs groups at the LHC is now to study the properties of this new resonance and determine if it is indeed the state predicted by the SM.
With this new discovery, a program of precision measurement involving the Higgs boson has just started and will get improved with the new run of LHC and future e + e − Linear Collider (LC). In fact, the 7 ⊕ 8 TeV data allow the measurement of the Higgs couplings to gauge bosons and τ + τ − with about 20-30% of precision while the Higgs couplings to bb and tt still suffer large uncertainties of about 40-50%. All these measurements will be improved with the new run of LHC at [13][14] TeV and the future e + e − LC.
In order to confirm that the discovered Higgs-like particle is the SM Higgs responsible for the electroweak symmetry breaking, we need to know all its couplings to SM particles with accurate precision and also measure the trilinear and quartic self-couplings of the Higgs in order to be able to reconstruct the scalar potential. In this regards, the LHC with its high luminosity option have also the capability of measuring the SM triple Higgs couplings through one of the following channels gg → hh → bbγγ, bbτ + τ − , bbW ± W ∓ * [8,9]. The e + e − LC, which will provide some high precision measurement of the Higgs mass and its properties such as couplings to SM particles and quantum numbers, would also be able to perform SM triple Higgs coupling through e + e − → Zhh (double Higgs-Strahlung) and e + e − → ν eνe hh (WW fusion) with more than 700 GeV center of mass energy with better precision [10]. In the double Higgs-Strahlung process Zhh the Z boson will be reconstructed from l + l − or qq pairs, while for the WW fusion process the two Higgs can be reconstructed from bbbb or bb and W + W − .
The discovery of this Higgs-like particle resonance opens a new era in elementary particle physics and leads to several theoretical and phenomenological studies on Higgs physics both in the SM and beyond. One of the very simplest extension of SM is the IHDM proposed, more than three decade ago, by Deshpande and Ma [11] for electroweak symmetry breaking purpose. Recently, the IHDM model has been very attractive because it provides a dark matter candidate [12,13], generates tiny neutrino masses [14] and also solve the naturalness problem [15]. Phenomenology of IHDM have been extensively studied during last decade [13,16,17].
The aim of this paper is to study the effect of the one loop radiative corrections to the triple Higgs coupling hhh as well as hZZ coupling in the framework of the IHDM. We will also compute the well known SM radiative corrections to the triple Higgs coupling as check of our procedure. Once these effect are well studied, we then proceed to the evaluation of radiative corrections to the double Higgs Strahlung process e + e − → Zhh. For this purpose, we will apply an on-shell renormalization scheme to evaluate these one-loop corrections.
For the numerical evaluation, we will take into account all the theoretical and experimental constraints on the scalar sector of the Model.
The paper is organized as follow: in the second section we introduce the IHDM model and describe the theoretical and experimental constraints that the model is subject to.
In the third section we introduce the on-shell renormalization scheme for the triple Higgs coupling hhh and hZZ in the IHDM and present our numerical results in the fourth section.
Numerical analysis of the double Higgs-strahlung is presented in the fifth section. Our conclusion is given in the last section.

A. The Model
The IHDM is one of the most simplest models for the scalar dark matter, a version of a two Higgs double model with an exact Z 2 symmetry. The SM scalar sector is extended by an inert scalar doublet H 2 which can provide a stable dark matter candidate. Under Z 2 symmetry all the SM particles are even while H 2 is odd and it could mix with the SM-like Higgs doublet. We shall use the following parameterization of the two doublets : In the above potential there is no mixing terms like µ 2 12 (H † 1 H 2 +h.c) because of the unbroken Z 2 symmetry. By hermicity of the potential, all λ i , i = 1, · · · , 4 parameters are real. The phase of λ 5 can be absorbed by a suitable redefinition of the fields H 1 and H 2 , therefore the scalar sector is CP conserving. After spontaneous symmetry breaking of SU (2) L ⊗ U (1) Y down to U (1) Q , the spectrum of this potential will have five scalar particles: two CP even H 0 and h which will be identified as the SM Higgs boson, a CP odd A 0 and a pair of charged scalars H ± . Their masses are given by: where λ L,S are defined as: This model involves 8 independent parameters: five λ, two µ i and v. One parameter is eliminated by the minimization condition and the VEV is fixed by the W boson mass.
Finally, we are left with six independent parameters which we choose as follow :
In order to have an inert vacuum the following constraint should be satisfied [19]: The extra scalar particles affects quantum corrections to the W and Z bosons self energies. formulas for ∆S and ∆T in the IHDM can be found in Ref [15,22].
Searches of scalar particles 1 of the IHDM at colliders [23] is not directly performed yet.
However, several studies [24][25][26] applied SUSY searches involving two, three or multiple leptons with missing transverse energy E miss T to the case of IHDM and set some limits on the dark Higges. We choose in our analysis 2 : Finally, the magnitude of a possible Higgs boson signal at the LHC is characterized by the signal strength modifier, defined as R γγ by : h SM denotes a 125 GeV SM Higgs boson. In our analysis below, while we will show points which satisfy theoretical and experimental constraints from our scans, we will highlight the points for which R γγ is consistent with the measured µ γγ at the LHC. The latest publicly available measurements read [27,28] µ CMS γγ = 1.13 ± 0.24 (12) µ ATLAS γγ = 1.17 ± 0.27 (13) III section V. Those couplings are given at the tree-level by: As one can see, both couplings hhh and hZZ involve only SM parameters. Those couplings receive corrections from one-loop diagrams. The one-loop effects from the SM particles have been studied in [30][31][32] for hhh and in [29] for hZZ. These effects are dominated by the top quark loops which does not exceed 10% for hhh and 1.5% for hZZ.
New physics effects to hhh coupling have been analyzed in the context of the Two Higgs Doublet Model [30] and the MSSM [32]. It was found that these effects can enhance significantly this coupling in a wide range of parameter space. Furthermore, these corrections depend on the model and hence any deviation from the SM tree level relation (14) by more than 10% would be an evidence for the presence of new physics.
The coupling hZZ have been analyzed in the framework of the two Higgs doublet model [31] and it has been found that the effect is rather small 1% to 2%.
We have calculated the radiative corrections to the tree level triple Higgs coupling hhh and hZZ both in the SM and IHDM in the Feynman gauge including all the particles of the model in the loops. The Feynman diagrams from IHDM contributing to λ hhh coupling are shown in Fig.( 1).
The one-loop amplitude are calculated using dimensional regularization. The calculation was done with the help of FeynArts and FormCalc [33] packages. Numerical evaluation of the one-loop scalar integrals have been done with LoopTools [34]. We have checked both numerically and analytically the UV finiteness of the amplitudes.
In order to do that, we have considered hhh and hZZ at one-loop level: i) for hhh, we considered the decay of an off-shell Higgs boson into a Higgs boson pairs Where q, k 1 and k 2 are the 4-momenta of the three particles satisfying on shell conditions k 2 1 = k 2 2 = m 2 h for final state Higgs pairs and an off shell condition q 2 = m 2 h for the decaying Higgs.
ii) For hZZ, we follow ref [31] and write: where k 1 and k 2 are the momenta of outgoing Z bosons. We assume that the decaying Higgs and one of the Z boson are on-shell q 2 = m 2 h , k 2 1 = m 2 Z while the other Z boson is off-shell k 2 2 = (m h − m Z ) 2 . Using power counting arguments, it is expected that M hZZ 1 receives the highest power contribution of the heavy fermions masses. Therefore, in what follow we will take into account only the M hZZ 1 form-factor to hZZ coupling.
Since we are dealing with a processes at the one-loop level, a systematic treatment of the UV divergences have to be considered. We will use the on-shell renormalization scheme in which the input parameters coincide with the physical masses and couplings [35]. In the on shell scheme, a redefinition of the fields and parameters is performed. This redefinition cast the Lagrangian into a bare Lagrangian and counter-term. The counter-terms are calculated by specific renormalization conditions which allow us to cancel the UV divergences of the diagrams with loops. Furthermore, since we have three Higgs as external particles and there is no mixing between the SM doublet H 1 and the inert doublet H 2 , we do not need to renormalize the particle content of the scalar potential of the IHDM. The tree level coupling hhh eq. (14) depends only on Higgs mass and the vev as in the SM, then the renormalization procedure will be the same as in the SM [35]. We redefine the SM fields and parameters as follow: where s W = sin θ W is the Weinberg angle and t = v(µ 2 1 − λ 1 v 2 ) is the tadpole which is zero at tree level once the minimization condition is used but will receives again finite radiative corrections at the one-loop level. To ensure that the VEV is the same in all orders of perturbation theory, it is well known that one need to renormalize the Higgs tadpole: i.e, all Higgs tadpole amplitudes T are absorbed into the counter-term δt. Thus, we put the first condition:T = δt + T = 0 The mass counter-terms are fixed by the on shell conditions [35]: The field renormalization constants are fixed by imposing that the residue of the two point Green functions to be equal to unity and the mixing γ-Z vanish for k 2 = m 2 Z . While the electric charge renormalization constant δZ e is treated like in quantum electrodynamics and is fixed from the e + e − γ vertex. The renormalized three point functionΓ µ e + e − γ satisfy at the Thomson limit:Γ µ e + e − γ ( p 1 = p 2 = m, q 2 = 0) = e Furthermore, the counter-term δs W can be obtained from the on-shell definition s 2 as a function of δm W and δm Z .
Inserting these redefinitions into the Lagrangian, we find the following counter term for hhh and hZZ [35]: By adding the un-renormalized amplitude for hhh and hZZ to the above corresponding counter-terms, one find the renormalized amplitudeŝ which becomes UV finite. For our numerical illustrations, we define the following ratios: WhereΓ hhh is the renormalized vertex.

A. SM case
In our numerical analysis, the parameters are chosen as follow : In the SM, the dominant contribution to ∆Γ hhh (SM ) comes from top quark loops [30,32].
We have computed the top contribution and shown that it is in perfect agreement with Ref. [30]. We have also isolated and evaluated the other SM contributions without fermions.
It turn out that this bosonic contribution is of the order of 5% for large q.
In Fig. (2)(left), it is illustrated that the top contribution to ∆Γ hhh (SM ) is negative before the opening of h * → tt threshold and also for q ≥ 700 GeV. It is clear that for large q, ∆Γ hhh (SM ) is dominated by top-quark contribution.
In Fig. (2)(right), we show the radiative corrections to hZZ in the SM. We present separately the fermionic corrections which are dominated by the top contributions and the bosonic contributions. The total corrections to hZZ is of the order of 2%. In this plot, we also shift the triple Higgs SM coupling λ SM hhh by λ SM hhh (1 + ∆), where ∆ represent any deviation from SM coupling. As one can see from the green line, the sensitivity to ∆ is rather mild. Due to custodial symmetry, it is expected that hW W coupling will enjoy similar effect as hZZ and that is why we illustrate only the case of hZZ.

B. IHDM case
Here, we will show our numerical analysis for the triple coupling of the Higgs in the IHDM taking into account: unitarity, perturbativity, false vacuum as well as vacuum stability constraints described above. We take the mass of the SM Higgs m h = 125 GeV and the masses of the inert particles to be degenerate, i.e : parameters, we perform the following scan: We plot in Fig. (3)(left) the relative corrections to the triple coupling as a function of λ 2 . The theoretical constraints put a limit on λ 2 parameter which is λ 2 ≤ 4π 3 . One can see from Fig. 3 that the corrections are maximized for λ 2 ≤ 2 and decrease for λ 2 > 2. In our following analysis, We will take λ 2 = 2 in order to maximize the effect from λ 2 .
In Fig. (3)(right), we plot the relative corrections to the triple coupling hhh in the plane (m Φ , µ 2 2 ) for a fixed q = 300 GeV and λ 2 = 2. One can see that the corrections are very important in a large part of the parameter space with an enhancement up to 280% for large values of m Φ and negative µ 2 2 . Furthermore, these corrections are increasing, for a fixed value of m Φ , while µ 2 2 is decreasing. The maximum of the corrections is reached for µ 2 2 ≈ −30000 (GeV) 2 . Moderate or very small corrections which can be in the range [−50, 50]% are also possible for large area of parameter space with low m Φ ≤ 300 GeV and any positive µ 2 2 . It is also important to note that LHC constraint from diphoton at the 2σ level exclude light charged Higgs 100 < m H ± < 175 GeV and negative µ 2 2 : left-down corner of the scatter plot. In Fig. 4 we show the relative corrections ∆Γ hhh (IHDM ) as a function of the momentum of This is visible on the left panel of Fig. 4 where we can see a kink for q = 1000 GeV which correspond to threshold effect h * → ΦΦ with m Φ ≈ 500 GeV. As it is shown, negative values for µ 2 2 gives large corrections to the triple Higgs coupling. This is because in our assumption of taking degenerate Higges m H 0 = m A 0 = m H ± = m Φ one can show that λ 4 = λ 5 = 0, 2 ) and therefore the triple coupling are given by It is clear that those couplings gets stronger for negative µ 2 2 . We now examine the effect of the radiative corrections on the triple coupling in the case where the invisible decay h → HH is open. This is illustrated in Fig. 5(left) and right. In Fig. 5(left) we impose both |λ L | < 0.02 required by dark matter constraints as well as best fit limit on the invisible branching ratio Br(h → invisible) ≤ 10%(left) and Br(h → invisible) ≤ 20%(right) [36].
It is clear from left panel that with dark matter constraint, the size of the corrections and the range of µ 2 2 are smaller than in the previous case where the invisible decay was closed. The large corrections observed for high m Φ are mainly due to the charged Higgs loops.
In order to exibit the decoupling behavior on the triple Higgs coupling, We increase both the range of µ 2 2 to be [−10 6 , 10 7 ] GeV 2 as well as the range of m Φ ∈ [0.1, 3] TeV. We see that the decoupling effect occurs when appropriate combinations of the involved parameters are taken large compared to the electroweak scale. We find that ∆Γ hhh reaches its maximum for m Φ ≈ 500 GeV and decrease to SM value for large m Φ .
In Fig. (7) we illustrate the IHDM effect on hZZ coupling. Similar to the triple Higgs coupling, we fix λ 2 = 2 and scan over µ 2 2 and m Φ = m H = m A = m H± as in eq. (25). In Fig. (7)(left) we show scatter plot for ∆Γ hZZ in the plan (µ 2 2 , m Φ ), contrarily to the triple Note that e + e − → hhZ process arise in s−channel only and hence its cross section can be probed more efficiently at low energies above the threshold (typically between 350 GeV and 500 GeV). While, at high energies, for √ s 700 GeV, the trilinear Higgs-coupling is better probed through the process e + e − → ννhh assuming the SM. In Fig. (8)  We study in this section the effects of the one-loop radiative corrections to the Higgs trilinear self-coupling calculated in the previous section on the double Higgs-Strahlung process via Z boson exchange. In the context of the SM, the O(α) electroweak corrections have been studied in [37] and it was found that these corrections are of the order 10%. However, these loop effects can be very large in beyond SM enhancing the total cross section by about 2 orders of magnitude in particular in models with extended Higgs sector. As outlined above, at the tree level, e + e − → hhZ have four diagrams as depicted in Fig. (8).
Only the first diagram contribute to the signal while the others are considered as a background. In our analysis, we include one-loop correction only to the triple Higgs coupling hhh which is expected to give sizeable contribution. Only this correction contribute to the signal. Therefore, we did not include corrections to the initial state vertex e + e − Z, to the self energies of Z-Z and γ − Z mixing, to the hZZ coupling and also we did not include the initial state radiation. In fact these corrections are well known in the SM and are not expected to deviate that much in the IHDM as we already show in the previous section for hZZ coupling. Moreover, we will not include corrections to e + e − → hhZ coming from boxes and pentagon diagrams. The one-loop amplitude can be written as follow: The squared amplitude at the one-loop level is then : Thus, the cross section is written as: where q 1 and q 2 are the four momentum of the incoming electron and positron, p 1 , p 2 and p 3 are the four momentum of the outgoing particles, and the factor 1 (2π) 2 arises from the flux of the initial particles.
For our studies, we define the ratio ∆σ by : Where σ total = σ tree + σ loop . This ratio measure the relative correction of the IHDM to the cross section with respect to the tree level result. As stated before, in our analysis we will take into account the theoretical and experimental constraints discussed in the second section assuming the parameters to rely in the range given by eq. (25). The phase space and evaluation of the one-loop squared amplitude has been performed with FormCalc [33] with the help of LoopTools to evaluate numerically the oneloop scalar integrals. We have used the same on-shell renormalization scheme explained in the previous section.
In Fig. 10, we have plotted the ratio ∆σ versus m Φ for center of mass energy √ s = 500 GeV.
We assume again that: λ 2 = 2, the dark Higges to be degenerate m H = m A = m H± = m Φ and perform a scan over µ 2 2 and m Φ . From this plot one can see that the corrections can reach 160% for high dark Higgs masses m Φ ≈ 500 GeV and are negative for low dark Higgs masses 100 GeV ≤ m Φ ≤ 270 GeV depending on the value of µ 2 2 . One can see that the suppression of the total cross section can reach −15% for m Φ = 140 GeV while the enhancement is predominant in most part of the parameter space.
To understand this, let us remind first that in large area of parameter space the correction to the triple Higgs coupling ∆Γ hhh is positive (see section IV Fig. (3)). Moreover, according to the plot Fig. (9), if this correction is positive this lead to an enhancement of the total cross section and vice-versa.
This explain that in most of the case, the corrections to the cross section are positive and confirm that the behavior of the ∆σ is consistent with our analysis concerning the trilinear Higgs self coupling in the IHDM. We stress that the enhancement of ∆σ is observed in a large part of the parameter space and can exceed 100% only in the high mass region m Φ ≥ 400 GeV for µ 2 2 < 0. We plot in Fig. 10(right) the ratio R γγ as a function of λ 3 and showing the relative corrections in the left column. For small values of −1 < λ 3 < 2 (low m Φ ) the corrections are quite small. For λ 3 > 4 (high m Φ ), ∆σ becomes more important and exceed 100%, this region corresponds to R γγ ≈ 0.9 ± 0.02. Note that our results concerning R γγ are in agreement with the results of [22].

VI. CONCLUSION
We have computed the radiative corrections to triple Higgs coupling hhh, hZZ coupling as well as e + e − → Zhh in the framework of inert Higgs doublet model taking into account theoretical and experimental constraint on the parameter space of the model. The calculation was done in the Feynman gauge using dimensional regularization in the on-shell scheme.
In the SM it is known that the top contribution to hhh coupling is of the order 10%, we found that the bosonic contribution is somehow significant and goes up to 5%. In the IHDM, we found that the total radiative corrections to the triple Higgs coupling could be substantial exceeding 100% for heavy dark Higgs masses m H , m A and m H± . We also show that the corrections to hhh are decoupling for large m Φ and large µ 2 2 . In the case of hZZ coupling the effect is rather mild and do not exceed 2.5%. We also evaluate radiative corrections to the double Higgs strahlung process e + e − → Zhh by looking only to the correction to the diagram that contribute to the signal i.e the triple coupling hhh. We have shown that the correction are also very important. In general, the size of the loop effects, typically large, makes their proper inclusion in phenomenological analyses for future e + e − LC indispensable.  Γ loop hhh (q 2 , m 2 h , m 2 φ ) ≈ 1 8αeπ 3 q 2 3απq 2 m W s W 2λ 2 3 + 2λ 3 λ 4 + λ 2 4 + λ 2 5 x 1 log( −1 + x 1 x 1 ) + x 2 log( −1 + x 2 x 2 ) + 2λ 3 3 + 3λ 2 3 λ 4 + 3λ 3 λ 2 4 + λ 3 4 + 3(λ 3 + λ 4 )λ 2 5 m 2 W s 2 W log 2 (− where x 1,2 are given by: x 1,2 = 1 ∓ 1 − 4m 2 φ /q 2 2 (40) [1] F. Englert