The hMSSM approach for Higgs self-couplings revisited

We compare the decay of the heavy Higgs boson into two SM-like Higgs bosons, $H\to hh$, calculated in a Feynman-diagrammatic approach at the one-loop level based on the one hand on the full effective potential involving the top quark and stops in the Minimal Supersymmetric Standard Model (MSSM) accompanied by the matched Two-Higgs-Doublet Model (2HDM) as its low-energy limit and on the other hand on the hMSSM approximation. We identify missing contributions due to the top quark in the Higgs self-couplings of the hMSSM, that - when taken into account - lead to a good agreement between the hMSSM and a full MSSM calculation, at least in the limit of the Higgsino mass parameter $\mu$ being small compared to the stop spectrum. We also thoroughly analyze momentum-dependent and kinetic corrections intrinsic to the Feynman-diagrammatic approach and the matching to the effective Lagrangian, respectively, for both our calculation in the MSSM and the hMSSM and for the latter suggest to include additional corrections from the top quark, which are independent of the unknown supersymmetric spectrum.


Introduction
The discovery of the Higgs boson with a mass of (125.09 ± 0.24) GeV [1] in 2012 by the LHC experiments ATLAS [2] and CMS [3] has marked a milestone for particle physics. While this structurally completes the Standard Model (SM) the SM itself leaves open many questions that require extensions of the model. The SM is therefore considered as an effective low-energy description of a more complete model valid at high-energy scales. Since the discovered Higgs boson behaves very SM-like any such beyond-the-SM theory has to contain a SM-like Higgs boson with a mass of about 125 GeV.
The tree-level Higgs sector can be described by two parameters, usually chosen to be the mass of the CP-odd Higgs boson, M A , and the ratio of the two vacuum expectation values of the two Higgs doublets, tan β = v 2 /v 1 , in the case of real supersymmetric parameters. Supersymmetry restricts the tree-level mass of the lightest CP-even scalar h to values below the Z boson mass M Z . This constraint is relaxed by the inclusion of radiative corrections in the Higgs sector that can shift its value to the measured 125 GeV. The dominant corrections originate from third generation quark/squark loops. Depending on the parameter choices, the squark masses must be quite large in order to match the observed Higgs mass value for small values of tan β. Moreover, in a significant part of the MSSM parameter space the limits on the squark masses are pushed into the TeV range by the unsuccessful LHC searches for supersymmetric (SUSY) particles so far. The loop-corrected Higgs sector depends on many SUSY parameters so that the investigation of the MSSM parameter space becomes a complicated task. This triggered the introduction of benchmark scenarios that are used by the experimental collaborations for the interpretation of their results. Among these, the hMSSM presented in Refs. [17][18][19][20] exploits the fact that the dominant corrections to the lightest CP-even Higgs mass and the mixing parameters that enter the Higgs couplings have a common origin and that the dominant corrections stem from the top-quark and its supersymmetric partners, the stops.
In the hMSSM the measured Higgs mass value M h is taken as an input parameter in addition to M A and tan β. This removes the explicit dependence of the Higgs sector on other SUSY parameters through the radiative corrections. In its region of applicability, the hMSSM approach has been shown to describe the MSSM Higgs mass spectrum and mixing angle α of the CP-even sector very well [20][21][22]. In particular, it allows to probe the low tan β regime where a very high SUSY scale is required for the radiative corrections to be large enough to achieve 125 GeV for the light CP-even Higgs mass. The Higgs self-couplings that are related to the Higgs masses through the Higgs potential are also affected by large radiative corrections. In order to make reliable predictions, the large logarithms that appear in the corrections in case of very large SUSY masses have to be resummed using effective field theory (EFT) methods. In physical processes containing the trilinear Higgs self-couplings, like Higgs decays into a pair of lighter Higgs bosons, momentum-dependent corrections to the vertex and to the kinetic factors can become important. These are not taken into account in the EFT approach, however, and have to be computed through a diagrammatic fixed-order calculation.
In this paper, we revisit the hMSSM approach with focus on the Higgs-to-Higgs decay of the heavier H into two SM-like Higgs bosons, H → hh. We compute the decay at next-toleading order (NLO) taking into account the dominant radiative corrections from the top-quark and stop sector. The calculation is performed in an effective low-energy 2HDM with MSSMlike quartic couplings that are properly matched to the MSSM and in the MSSM itself. In both cases the calculation is performed in the Feynman-diagrammatic approach thus including momentum-dependent corrections. Moreover, radiative corrections to the Higgs self-couplings from the top-quark contributions in the 2HDM, and the top-quark and stop contributions in the MSSM, are taken into account through effective couplings. By choosing appropriate counterterms according to the low-energy limit, double counting is avoided when including the diagrammatic NLO corrections. By plugging in the effective trilinear Higgs self-coupling of the hMSSM and comparing with the full MSSM result, we are able to disentangle the deviations due to the hMSSM approximation of the coupling from those originating from momentum-dependent contributions. In this way, we are able to properly dissect the Higgs self-coupling of the effective hMSSM approximation and to propose improvements that allow to better approximate the full result. It turns out that the bulk of the improvement does not introduce additional parameters so that it is appropriate to dub it "improved hMSSM".
The outline of the paper is as follows. In Section 2 we introduce the MSSM, the 2HDM and the hMSSM approach and define our notation. In particular in Subsection 2.4 we discuss the effective potential as part of the low-energy effective Lagrangian, while Subsection 2.8 presents our proposed improvement of the hMSSM approach. Section 3 contains the explicit computation of the NLO decay width for the decay H → hh with effective couplings. Section 4 is dedicated to the presentation of our numerical results. Our conclusions are given in Section 5.

The MSSM Higgs sector as an effective 2HDM
In this section we first introduce our notation of the MSSM Higgs sector. If the supersymmetric spectrum is heavy, the MSSM Higgs sector can be understood as an effective (properly matched) low-energy 2HDM, where all heavy (s)particles are integrated out. This approach allows to resum logarithms as they appear e.g. in the derivation of the Higgs masses at loop-level. In this section we use the Effective Potential Approach (EPA) [23,24] for the matching, define the relevant effective potential, kinetic corrections in the effective Lagrangian and provide results for loop corrections to the Higgs masses and the Higgs self-couplings at order O(α t ), where α t = y 2 t /(4π) with the top-Yukawa coupling y t = √ 2m t /v (v ≈ 246.22 GeV is the SM vacuum expectation value and m t the top-quark mass.). We explain the hMSSM approach and identify terms that are missing in the Higgs self-couplings of the hMSSM approach.

The MSSM Higgs and squark sectors at tree-level
Supersymmetry requires the introduction of at least two complex Higgs doublets, which is realized in the minimal SUSY version, the MSSM. In addition, the Adler-Bell-Jackiw anomaly contributions [25,26] due to the Higgsino doublet states cancel each other thanks to their opposite hypercharges. The two doublets H u and H d with hypercharge Y = 1 and −1, respectively, can be expressed in terms of the charged and neutral components φ ± i and φ 0 i (i = u, d), as The tree-level potential for the MSSM Higgs fields is derived from the F -and D-term contributions and the soft-SUSY breaking Lagrangian. The most general supersymmetric potential for two Higgs doublets at tree level reads where 12 = − 21 = 1 and m 2 H d , m 2 Hu and Bµ denote the corresponding soft-SUSY breaking mass parameters. The SU(2) L and U(1) Y gauge couplings are given by g and g , respectively. After electroweak symmetry breaking we expand the neutral fields around the vacuum expectation values (VEVs) according to The ratio of the VEVs v u and v d is defined as tan β = vu v d while obeying the sum rule v 2 u +v 2 d = v 2 . Throughout our work we fix m 2 H d and m 2 Hu through the radiatively corrected tadpole equations.
The mass matrix of the stop sector using left-and right-handed stopst L andt R is given by where Mt L and Mt R are the left-and right-handed soft-SUSY breaking mass terms, respectively. Mt L equals the soft-SUSY breaking mass term MQ L of the doublet of the third generation squarks. The stop-mixing parameter X t is defined through X t = A t − µ/t β involving the trilinear soft-SUSY breaking stop parameter A t and the µ-term, which has already been part of the tree-level potential in Eq. (2). In the subsequent derivation of the effective potential we work in the gaugeless limit, which is why we also omit D-terms proportional to M 2 Z in the stop mass matrix. Its diagonalization yields the stop masses mt 1 and mt 2 , which for M S := Mt L = Mt R are given by m 2

The 2HDM Higgs sector at tree-level
Being the low-energy limit of the MSSM, we work only in the type-II 2HDM, where H 1 and H 2 couple to down-type and up-type quarks, respectively. The tree-level Higgs potential including only terms that arise in the MSSM at tree-level takes the form We label the fields H 1 and H 2 instead of H d and H u to follow the standard notation. Both field conventions follow the MSSM field H u , i.e. 2). For the relation between the Higgs fields in the MSSM and the 2HDM we also refer to Ref. [27]. The parameters m 2 1 and m 2 2 can again be fixed through the tadpole equations, such that in a generic 2HDM the potential has six free parameters, which in the above λ-basis are m 2 3 , λ 1 − λ 4 and the ratio of the VEVs, tan β. The trilinear Higgs self-couplings within the λ-basis are given by In the MSSM at tree level the couplings λ i are given by , so that λ 34 ≡ λ 3 + λ 4 = −λ 1 = −λ 2 (12) and m 2 3 = Bµ which is thus related to M 2 A , see Eq. (4).

The effective 2HDM Higgs sector beyond tree-level
Beyond tree level the masses and couplings of the 2HDM discussed in the previous section receive quantum corrections from heavier (s)particles 1 and thus form an effective low-energy 2HDM. This effective model is obtained through a proper matching at the scale, where the heavier (s)particles are integrated out. Practically, this matching can be performed by a Feynmandiagrammatic calculation, in which all squared external momenta p 2 are strictly set to zero. This limit of vanishing momenta defines the EPA [23,24]. This approach includes all one-particle irreducible diagrams and allows to define an effective 2HDM as the low-energy limit of the MSSM.
Within the rest of this section we will define the effective 2HDM relevant for our purposes and later supplement it with a Feynman-diagrammatic calculation of the decay H → hh.
The EPA allows to calculate the corrections to the potential V LO 2HDM (= V LO MSSM ) in Eq. (10) and thus defines the effective potential V eff 2HDM . 2 A detailed discussion follows in the next section. However, also the kinetic term of the effective Lagrangian receives corrections, such that we can decompose the effective low-energy 2HDM Lagrangian as which involves an additional kinetic matrix Z eff . The matrix Z eff can be obtained, as in Ref. [27], by an expansion of the (off-)diagonal H i H j (i, j = 1, 2) two-point functions in their external momenta. In our Feynman diagrammatic calculation Z eff appears as an effective wave-function renormalization and is named kinetic correction. We will present its detailed form in Section 3.1.
We note that beyond tree level also couplings λ 5 to λ 7 are generated, see e.g. in Ref. [27]. Their generation is connected to a non-vanishing value of µ [28]. However, since we work with physical quantities as the Higgs self-couplings, masses and mixings derived from the full effective Higgs potential we do not need their explicit form in this work since they are taken into account intrinsically. The MSSM Higgs masses in the EPA were first calculated at order O(α t ) in Refs. [29][30][31][32], followed by calculations including O(α b ) and electroweak corrections in Refs. [33][34][35][36]. In this approach also the large radiative corrections to the Higgs self-interactions are known at one-loop order O(α t,b ) [37][38][39] and two-loop order O(α t α s ) [40]. In the hMSSM approach, which is anyhow agnostic to the exact form of the higher-order corrections, the corrections to the Higgs masses and the Higgs self-couplings are reexpressed in terms of the light Higgs mass M h only, see Section 2.7. In both cases, the proper calculation of the 2HDM effective potential and the hMSSM approach, we work with an effective 2HDM.

The 2HDM effective potential
As we argued beforehand the application of the EPA defines the Higgs potential of an effective low-energy 2HDM, which we provide explicitly in this section. We consistently distinguish quark from squark contributions, in order to later omit the effect of stops.
We are only interested in the O(α t ) corrections to the effective potential of the 2HDM in view of the corrections that are implicitly taken into account in the hMSSM approach. To stay strictly at O(α t ) in perturbation theory, we are working in the gaugeless limit and thus drop supersymmetric D-term contributions beyond tree-level. This limit captures the dominant corrections and has the advantage that our discussion of H → hh remains independent of the renormalization of electroweak parameters, in particular the vacuum expectation value(s), see Section 3.1. We also omit corrections from the bottom and sbottom sector. We split the O(α t ) corrections to the effective potential in top-and stop contributions as follows: The effective potential of the 2HDM is then given by The field-dependent mass parameters are defined as Therein h t denotes the top-Yukawa coupling h t = √ 2m t /v u . If not mentioned otherwise all parameters are running parameters evaluated at the scale Q. The coefficient C ε = Γ(1+ε)(4π) ε expands to 1 + (−γ E + log(4π))ε for small ε and results in the ultraviolet divergent term The renormalization scale Q is a priori not fixed but should be of O(M S ). It should be noted that the scale Q represents the matching scale between the full MSSM and the low-energy 2HDM. Below this scale the top quark and the stop states are integrated out and thus do not contribute to the low-energy running of the Higgs self-interactions, i.e. they are decoupled. On the other hand the low-energy self-couplings λ ijk develop a residual scale dependence due to the light particles still present in the low-energy 2HDM spectrum.
We are interested both in the corrections to the Higgs masses and to the triple Higgs selfcouplings in order to compare differences directly. Some of our subsequent discussion is a historical review, but allows to understand the underlying basis of the additional terms that we add to the hMSSM approach.

Corrections to the Higgs masses
In this section we discuss the corrections to the Higgs masses obtained from the 2HDM effective potential. Within the discussion we keep top-quark and stop contributions separated. In order to obtain the CP-even Higgs masses, the 2HDM effective potential is expanded in the CP-even and CP-odd components of the neutral fields φ 0 d and φ 0 u , see Eq.
(3). 3 The mass corrections are then obtained as second derivatives with respect to these components. A nice explanation of this procedure is given in Ref. [41], which reproduces the corrections presented in the original publications [29][30][31][32]. The symmetric mass corrections ∆M 2 ij in gauge eigenstates are added to the tree-level mass matrix M 2 tree in Eq. (5). We split them according to in top-quark and stop-induced corrections originating from V NLO (t) and V NLO (t), respectively. The top-induced correction contributes only to ∆M 2 uu (t). The corrections are given by where we have used the abbreviations Two remarks are in order: Only through a non-vanishing µ-term also the corrections ∆M 2 dd and ∆M 2 du are non-zero. The ultraviolet divergent terms and the renormalization scale dependence cancel in the sum of top-quark and stop-induced corrections: It is therefore obvious that an effective 2HDM, where both the top quark and stops are taken into account in the effective potential, yields finite corrections at NLO. In contrast in an effective 2HDM, which only involves tops, also ultraviolet divergences appear, which can only be cured through a proper renormalization. We will discuss this when we will later calculate our partial decay width.

The approximation for the light Higgs mass
Our previous discussion of the loop corrections was carried out in gauge eigenstates. We can rotate to mass eigenstates to obtain the correction to the light Higgs mass, ∆M h . We now consider the case where left-and right-handed soft-SUSY breaking mass parameters M S = Mt L = Mt R are identical. We expand the mass corrections in inverse powers of M S and subsequently present the result for ∆M 2 uu and the total mass correction ∆M 2 h : Therein, we expressed the SM VEV through the Fermi constant G F . Moreover for ∆M 2 h we employed the decoupling limit α → β − π/2, which is associated with M A M Z . The last well-known relation in Eq. (21), see Refs. [29][30][31][32] and later updates in Refs. [42,43], shows that the light Higgs mass grows with the fourth power of the top-quark mass and logarithmically with the stop masses associated with the SUSY scale M S . We define the parameter = ∆M 2 uu . Setting µ = 0 and thus X t = A t we see that exactly corresponds to the correction ∆M 2 h /s 2 β in accordance with the fact that all other elements of ∆M 2 ij in gauge eigenstates vanish. The correction increases the upper mass bound to The mass of the lightest Higgs boson is explicitly given by Within this approximation the masses of the heavy neutral and charged Higgs bosons are obtained by sum rules, and the effective mixing parameter α between the CP-even scalars is given by We will later come back to the approximation also for the Higgs self-couplings, but before introduce the hMSSM approach.

The hMSSM approach
As argued before, in the limit of a µ-term, which is small compared to the stop spectrum, the corrections and ∆M 2 h /s 2 β in Eq. (21) are identical. Then, the complete correction to the light Higgs mass originating from the top-quark and stop sector only enters the element ∆M 2 uu . The hMSSM approach [17][18][19][20] assumes all supersymmetric particles to be heavy and is agnostic for what concerns the origin of the mass corrections. Instead in the hMSSM one obtains ∆M 2 uu by inverting Eq. (23) using M h as an input parameter, which yields The heavy Higgs mass M H and the mixing angle α are then in turn also fixed to the values in Eq. (24) and Eq. (25), respectively. The limitations of this procedure are rather obvious: First, the corrections of the top-quark and stop sector are assumed to be dominant, which is true at low values of tan β, while at larger values of tan β also corrections from the bottom and sbottom sector provide a potentially relevant (subleading) contribution. Second, neglecting the corrections ∆M 2 dd and ∆M 2 du is only compatible with a µ-term, which is small compared to the stop spectrum and thus implies relatively light Higgsinos. This means not all supersymmetric particles are necessarily heavy. They can influence Higgs physics mainly through decays, either by allowing for additional decay channels for heavy Higgs bosons or through contributions to loop-induced decays of the Higgs bosons.

The improved hMSSM for Higgs self-couplings
We are interested in the corrections to the triple Higgs self-couplings presented in Eq. (8) that emerge from the effective potential in Eq. (14). Corrections are obtained by performing the third derivatives with respect to the corresponding Higgs fields. Again, we only focus on corrections from the top-quark and stop sector and thus split the individual contributions as follows both in gauge eigenstates {i, j, k} ∈ {d, u} and in mass eigenstates {i, j, k} ∈ {h, H}. We present them subsequently in gauge eigenstates. The top-quark contribution enters ∆λ uuu , while for a non-vanishing µ-term the stop contributions yield a correction to all couplings in accordance with the original results in Ref. [37]: The ultraviolet divergences, which are only part of ∆λ uuu , cancel between the top-quark and stop contributions ∆λ uuu (t) and ∆λ uuu (t), respectively. So does the renormalization scale dependence related to the top-quark and stop contributions. It is apparent that in contrast to the calculation of the Higgs mass terms the top-quark induced correction develops an additional constant factor (− 2 3 ) at the one-loop level. Since we are interested in the decay H → hh we combine the results above to obtain ∆λ Hhh by rotating with the CP-even mixing angle α. This yields Therein we introduced the additional abbreviations For degenerate soft-SUSY breaking masses M S = Mt L = Mt R and in the limit µ M S the expansion in inverse powers of the SUSY scale M S implies Here our expansion in inverse powers of M S adds terms proportional to m 6 t and all relevant terms that are of order 1/M 2 S for A t ∼ M S . If we perform a similar expansion in M 2 uu in Eq. (20), we obtain In ∆λ Hhh we can identify the corrections that are part of in Eq. (32), and therefore write The hMSSM approach advocates to just use the first term in Eq. (33) as a correction not only for the Higgs masses, but also for the Higgs self-couplings and thus misses the second bracket. However, the second bracket includes a purely top-induced contribution, which originates from the top-quark correction to ∆λ uuu (t) and should not be missed in the Higgs self-couplings, neither in the approximation nor in the hMSSM approach. The other terms of the second bracket are instead well suppressed for heavy squark masses and in the spirit of the hMSSM approach can be neglected. This constant correction beyond the terms comprised in the approximation also appears in all other Higgs self-couplings. We therefore define and obtain effective couplings of the form where the tree-level couplings are taken from Eq. (8). The usage of the effective couplings λ ijk can be considered an improvement of the original hMSSM approach, which is why we dub it "improved hMSSM". Additionally, we will later use the couplings λ ijk , which are defined as in Eq. (35), but with instead of . Thus, they correspond to the original hMSSM approach. Note that and therefore also can be calculated either from the actual correction = ∆M 2 uu , which equals the approximation, or according to the hMSSM approach. We follow the hMSSM approach and thus obtain from the right-hand side of Eq. (26) and from Eq. (34).
We summarize that in our subsequent calculation of the partial decay width H → hh we work with effective low-energy 2HDM couplings and mixing angles and thus do not only employ the tree-level coupling λ Hhh in the leading-order amplitude. We instead use the effective couplings λ Hhh and λ Hhh of the original and the improved hMSSM approach, respectively, and, moreover, make use of the complete correction in the effective potential given by that corresponds to the proper matching of the low-energy 2HDM to the MSSM. We add the arguments (t,t) to λ eff Hhh , which have to be understood as flags, i.e. we allow to add the topand stop-induced contribution to ∆λ Hhh separately. In practice at LO in the decay width we will always use λ eff Hhh (1, 1), which is ultraviolet finite in contrast to λ eff Hhh (1, 0).

The partial decay width H → hh
For now we have discussed both the MSSM and its approximation, the hMSSM, in the effective potential approach, which allows us to match them to an effective 2HDM. In our subsequent discussion the couplings and masses of the 2HDM include the previously mentioned one-loop corrections. We explicitly provided formulas for the corrected Higgs self-couplings and note that the Higgs masses are obtained by diagonalizing the one-loop corrected mass matrices in the EPA, see also Ref. [40]. For the scope of this work we neglect the contributions of other particles than the top quark and squarks that, however, could be taken into account in a full diagrammatic calculation in a straightforward manner. In the same context we will also neglect the residual RGE-evolution of the Higgs self-couplings λ ijk within the low-energy 2HDM in the following, i.e. work with the values obtained at the matching scale. In a full calculation the running due to the light degrees of freedom would have to be taken into account for consistency. Since the dominant radiative corrections to the decay H → hh are known to emerge from topquark and stop loops, these neglected effects are only subleading and do not contribute to the mismatch between the full MSSM and the 2HDM as the low-energy limit.
In this section we calculate the partial decay width H → hh. We perform a Feynmandiagrammatic calculation at the one-loop level including the full momentum-dependent corrections. We denote the momentum of the incoming Higgs boson by q H and the momenta of the outgoing Higgs bosons by q 1 and q 2 . Ultimately, we perform an on-shell calculation and thus set q 2 H = M 2 H and q 2 1 = q 2 2 = M 2 h . In the Feynman-diagrammatic approach we obtain the partial decay width Γ(H → hh) according to where A = A LO + A NLO denotes the amplitude. At pure tree-level the contribution A LO equals the Higgs self-coupling λ Hhh as given in Eq. (8). The corresponding Feynman diagram is shown in Fig. 1 (a). However, we want to work in the effective 2HDM, i.e. we employ the effective couplings λ eff Hhh and λ Hhh and thus incorporate higher-order effects already in A LO . We also apply the previously discussed one-loop corrections in the EPA, see also Ref. [40], to obtain the external Higgs masses.
The loop-corrected amplitude A NLO can be split into the following pieces where A virt (t,t) denotes the momentum-dependent one-particle irreducible Feynman diagrams. They are depicted in Fig. 1 (b)-(f). A ext (t,t) are external self-energy corrections adjusted to the amputated Green's functions including a mixing between H and h, see Fig. 1 (g)-(i). The contribution A ext,eff (t,t) originates from the kinetic mixing Z eff already discussed in Section 2.3 and thus provides the proper normalization of the (effective) Higgs fields in the effective lowenergy 2HDM. Finally A δλ (t,t) comprises additional counterterms induced by the effective couplings of the EPA and is generically depicted in Fig. 1 (j). A virt can be easily expressed in terms of Passarino-Veltman integrals [44,45]. We present the corresponding analytic expression in Appendix A. On the other hand we have to define a renormalization scheme to fix all the remaining counterterms. By adding the arguments (t,t), which have to be understood as a flag to include or to disregard top-quark and stop contributions, we emphasize that we can add both contributions separately to all individual ingredients of the one-loop amplitude, see Appendix A for an explanation.

Self-energy corrections and renormalization
We employ effective couplings and effective masses in an effective 2HDM at tree-level and also use them in our one-loop Feynman-diagrammatic calculation. This implies that beyond the corrections already implemented through the effective potential, we add only momentumdependent terms. First we discuss the external self-energies, which boil down to the following contributions where Σ ij (p 2 ) is the self-energy involving the two mass eigenstates i, j ∈ {h, H} and Σ ij (p 2 ) is its derivative with respect to the squared external momentum p 2 . The self-energy corrections δZ all enter A ext . Again these corrections can be expressed in terms of Passarino-Veltman integrals and can be split into top-quark and stop-induced corrections, see Appendix A. The self-energy corrections enter the amplitude as follows: A ext (t,t) =  Additionally we have to take into account the kinetic corrections explained in Section 2.3. They imply an additional contribution at the one-loop level that corresponds to the matching of the kinetic term of the full theory to the kinetic term of the effective low-energy 2HDM after integrating out the top-quark (stops). Since our calculation is based on the field definition of the full MSSM we have to divide through Z eff to obtain the field normalization in the effective low-energy 2HDM, which explains the subtraction of the corresponding terms δZ eff . The effective Z eff matrix is given by Also for the self-energies presented here, we can take into account top-quark and stop contributions separately.
The mixing angle α is renormalized by promoting the tree-level relation for α 0 = α + δα, see Eq. (6), to the one-loop level, which results in Using the explicit expressions for the mass corrections in the EPA, see Eq. (17), implies that the mixing angle is renormalized for vanishing external momenta in accordance with the consistent definition of effective low-energy parameters.
The tree-level coupling λ Hhh in Eq. (8) suggests that we need to renormalize not only the mixing angle α, but also the angle β as well as the electroweak sector, i.e. M Z and v. However, since we are working with effective parameters defined by the radiatively corrected EPA in our one-loop corrections, the renormalization of the parameters β, M Z and v is already part of the renormalization of the effective potential intrinsically, since the full counterterm δλ Hhh is defined by the EPA that determines the full mismatch between the MSSM and the lowenergy 2HDM. This is different for the mixing angle α that enters as the external rotation of the current eigenstates to the mass eigenstates applied to the fully corrected and renormalized effective potential in the current-eigenstate basis. We can therefore write The renormalization of the mixing angle α is formally part of A δλ , but adds to the nondiagonal renormalization factor δZ Hh , which is obvious from the relation ∂λ Hhh ∂α = λ hhh −2λ HHh . Since we employ an effective coupling in the tree-level amplitude A LO , we need to adjust the counterterm A eff accordingly to avoid double-counting. We obtain We explicitly checked the last relation, which is in accordance with Ref. [46]. Again all these relations hold for top-quark and stop contributions separately. It is obvious that the combination of the self-energy corrections A ext and A ext,eff and the counterterm A δλ only leaves momentum-dependent corrections in the amplitude A NLO .

Combining the results
We have now presented all relevant ingredients for the calculation of the partial decay width H → hh at the one-loop level within the effective 2HDM. In the following we will present our numerical results for different combinations of the effective or tree-level couplings and including top-quark and squark contributions separately. We emphasize that also just taking the top-quark contribution into account leads to an ultraviolet finite result in accordance with the Appelquist-Carazzone decoupling theorem [47]. According to Eq. (45) the contribution A eff (t,t) cancels all divergences of the vertex corrections of Fig. 1 (b)-(f). The contributions of the counterterm δα and the counterterms δZ eff ij together cancel the divergences of the external self-energies of Fig. 1 (g)-(i). All these cancellations emerge for the top-quark and stop contributions separately. On the other hand the sum of the top-quark and stop contributions leads to finite vertex corrections already before renormalization and in the same way to a finite counterterm A eff .
We summarize the options for our numerical analysis in Tab. 1. Option 1 is the MSSM calculation involving top-quark (and stop) loops without the use of the effective Higgs selfcoupling λ eff Hhh , but employing the tree-level values at LO and NLO. We name the corresponding decay widths Γ LO and Γ NLO (t,t). Through the arguments (t,t) we indicate whether only top quarks (1, 0) or additionally stops (1, 1) are included in the Feynman-diagrammatic calculation at the one-loop level. Option 2 is the same calculation using the effective Higgs self-couplings λ Hhh and λ Hhh within the original and our improved hMSSM approach, respectively. We make use of these effective couplings both at LO and NLO, which results in Γ LO / and Γ NLO / (t,t), respectively. For the latter to avoid double-counting of contributions, we subtract the vertex correction using the exact value of λ eff Hhh (t,t) in A eff in Eq. (45). This can be understood from the fact that the Feynman-diagrammatic calculation of the partial decay width H → hh adds the exact corrections due to the top quark and the stops. Lastly, in option 3 the calculation is performed in a consistently matched 2HDM to the MSSM with effective Higgs self-couplings and subleading terms in the matching beyond the hMSSM approach. This results in Γ LO eff and Γ NLO eff (t,t). Both for option 2 and option 3 we use λ ijk , λ ijk and λ eff ijk (1, 1), respectively, in all occurrences of Higgs self-couplings in Eq. (40) and Eq. (41). This choice is taking finite higher-order effects into account and does not harm the cancellation of ultraviolet divergences.

Comparison with earlier work
Our method of calculating the decay width for H → hh deviates from Ref. [46], since the latter performed a DR renormalization of all parameters involved, i.e. did not include an explicit decoupling of the heavy top-quark and stop states involving a proper matching to the full MSSM. In this way the work of Ref. [46] could not isolate the pure momentum-dependent contributions beyond the effective parameters rigorously.
The work of Ref. [48] presented the full NLO results within the complex MSSM with conventional SUSY-electroweak renormalization, i.e. without the introduction of effective Higgs self-couplings and mixing angles. This work finds large radiative corrections that, however, should be explainable as the missing contributions to the effective parameters at LO to a large extent. We will demonstrate this effect in our subsequent numerical analysis. Our approach can be extended to the full calculation within the MSSM. This, however, is beyond the scope of our work.
The recent work of Ref. [49] renormalized the Higgs self-coupling in the MS-scheme so that also this work did not perform an explicit matching of the low-energy 2HDM to the full MSSM. On the other hand the authors used the hMSSM approach to approximate the dominant radiative corrections within the MSSM Higgs sector. The residual effects beyond the use of their "effective" parameters, however, range at the same level as the consistent momentum-dependent contributions beyond the effective couplings and mixing angles as obtained in this work.

Numerical Results
We perform our comparisons in one benchmark scenario, in which we allow for different values of the Higgsino mass parameter µ. We will use the results of the work performed in Ref. [40] for the effective Higgs self-couplings and masses at the one-loop level. The most important parameter settings are as follows, where M S is the soft-SUSY breaking mass used for all left-and right-handed squarks. We set α s (M Z ) = 0.118, M Z = 91.15449 GeV, m t = 173.2 GeV and G F = 1.166378 · 10 −5 GeV −2 . The latter fixes the vacuum expectation value v. The top-quark mass is understood to be on-shell and internally transformed into a DR top-quark mass at NLO evaluated at Q = M S . This transformation as well as the running and matching of α s , see Ref. [40] for details, induces a dependence on the gluino mass and the other squark masses, which are fixed by the choices M 3 = 2500 GeV, m b = 4.84 GeV and A b = A t . The quark masses of the first two generations are set to zero. In the calculation of the Higgs masses and the effective Higgs self-couplings λ eff ijk (1, 1) at the one-loop level we also incorporate subleading D-terms in the squark masses to match the results of Ref. [40]. The determination of the D-terms also need the weak mixing angle, fixed through M Z above and M W = 80.36951 GeV. We checked that the inclusion of D-terms in the squark masses has almost no impact on our findings due to the relatively heavy SUSY scale, i.e. M Z M S .
The choice of X t allows to keep the light CP-even Higgs mass close to 125 GeV at least for values of tan β > 4. A Higgs mass of 125 GeV can also be reached for lower values of tan β, but only in combination with larger values of M S . Larger values motivate an RGErunning of the effective couplings. Since such a discussion is beyond our scope, we stick to M S = 1500 GeV and emphasize that our findings can be considered quite general, despite the fact that we "undershoot" the experimental value of 125 GeV significantly at low values of tan β. For a suitable and conclusive comparison we use the light CP-even Higgs mass M h obtained in the EPA at the one-loop level as input to the hMSSM approach, see Eq. (26). This implies that our values of and are based on a value of M h that differs from 125 GeV for low values of M A or tan β. In the following we choose three different values of µ. Our first choice is µ = 0 GeV, since only a Higgsino mass, which is small compared to the stop masses, allows the hMSSM approximation to be valid. We are aware that µ = 0 GeV results in very light Higgsinos, which are not compatible with chargino mass bounds [50]. Thus, the scenario with µ = 0 GeV is considered pedagogical to set up the consistency of our analysis with the hMSSM approach. The second and third choice are µ = 400 GeV and µ = 1000 GeV. For the latter choice, due to µ and M S being rather close, we will see remaining differences between the hMSSM approximation and the exact one-loop results in the Higgs masses, self-couplings and the partial decay width H → hh. They arise, since we intrinsically violate an assumption of the hMSSM, being µ M S .

Light Higgs mass and improved Higgs self-couplings
We display the value of the light Higgs mass M h in Fig. 2 for three values of µ = 0 GeV, 400 GeV and 1000 GeV. The Higgs mass is calculated at NLO taking into account the corrections described in Section 2.5. Since the value of X t is equal in all scenarios and thus are the stop masses, the light Higgs mass is almost identical in the three scenarios. We emphasize that for µ = 0 GeV the exact values of the heavy Higgs mass m H and the Higgs mixing angle α can be obtained through Eq. (24) and Eq. (25). For µ = 400 GeV and µ = 1000 GeV we depict the differences between the predictions of the hMSSM approach through Eq. (24) and Eq. (25) and their exact determination at NLO in Fig. 3. For the heavy Higgs mass such differences are mostly below 1% throughout the parameter plane. For the Higgs mixing angle α instead differences rise to a couple of percent, in particular for a very large value of µ = 1000 GeV, which is close to the soft-SUSY breaking mass M S = 1500 GeV in our example. Ref. [21] showed partially smaller discrepancies for the Higgs mixing angle α due to the smaller ratio µ/M S 1. We emphasize again that for µ = 0 GeV we find perfect agreement between the two approaches. Thus, we note that such discrepancies in m H and α are not the dominant source of differences observed in the partial decay width of H → hh in previous studies. On the other hand, as for the subsequently discussed Higgs self-couplings, the differences in α for larger µ/M S point towards the limitations of the hMSSM approach, see below.
We show the Higgs self-coupling λ Hhh in Fig. 4 again for µ = 0 GeV, 400 GeV and 1000 GeV. While the Higgs self-coupling obtained in the standard hMSSM approach λ Hhh shows large deviations from the exact value λ eff Hhh (1, 1), the improved version λ Hhh , performs significantly better, compare Fig. 4 (a) and (d). The kinks at tan β = 4 in Fig. 4 (d) are induced by our choice of X t , which is constant for tan β ≤ 4 and rescaled for tan β > 4. For non-vanishing µ the improved coupling shows remaining differences of a few percent, see

Momentum-dependent and kinetic corrections to H → hh
Before we discuss the quality of the hMSSM approximation for what concerns H → hh we focus on the relevance of momentum-dependent and kinetic corrections for the partial decay width H → hh in the benchmark scenario with µ = 0 GeV. They arise from our Feynmandiagrammatic calculation with non-zero external momenta and the kinetic mixing discussed in Section 2.3, respectively. The total partial decay width for H → hh, Γ NLO eff (1, 1), including all top-quark and stop corrections is depicted in Fig. 5 (a). First it is apparent that for low values of M A the heavy Higgs boson H is too light to decay into a pair of two CP-even light Higgs bosons hh, such that the on-shell decay is kinematically closed. Only for low values of tan β, where M h is decreasing rapidly, this decay mode is again of relevance. In all subsequent figures the region, where the on-shell decay H → hh is kinematically closed, is shown in white. We show the NLO partial decay width Γ NLO eff relative to the pure tree-level partial width Γ LO including the treelevel Higgs self-coupling λ Hhh in Fig. 5 (b). We observe large corrections in accordance with the literature, see e.g. Ref. [48]. However, we point out that using the effective coupling λ eff Hhh (1, 1) in the tree-level decay width Γ LO eff diminishes the one-loop corrections substantially, see Fig. 5 (c). Accordingly, Fig. 5 (c) demonstrates the relevance of the remaining momentum-dependent and kinetic corrections, which we incorporated through a Feynman-diagrammatic calculation and through the inclusion of the external kinetic mixing described in Section 2.3, respectively. Exactly such momentum-dependent and kinetic contributions were missing in the comparison of H → hh decay widths performed in Ref. [21]: Therein the partial decay width Γ LO was compared against a full Feynman-diagrammatic one-loop calculation including an additional resummation of large logarithms. The discrepancies found at the level of the partial decay width in Ref. [21] were between 15 − 25% in a large part of the M A -tan β-plane and thus in the same ballpark as the missing momentum-dependent and kinetic contributions. Keep in mind that the partial decay width Γ LO and thus the comparison in Ref. [21] missed also the constant factor between and . We leave a detailed comparison employing the newest predictions for a future study. We emphasize that the momentum-dependent and kinetic corrections in our benchmark scenario are mostly induced by the top-quark contribution, whereas the stops only yield a tiny fraction due to the relatively large stop masses. Given that in Fig. 5 we are not yet applying the hMSSM approximation similar results are obtained for larger values of µ. We therefore refrain from showing the corresponding results for µ = 400 GeV or µ = 1000 GeV.

H → hh in the hMSSM
We finally focus on the performance of the hMSSM approximation in the description of the partial decay width H → hh. In Fig. 6 we compare the partial decay widths Γ NLO (1, 1) and Γ NLO (1, 1), both named option 2, against the full result Γ NLO eff (1, 1), named option 3 in Section 3.2. In a first comparison we keep both top-quark and stop corrections in our calculation. For µ = 0 GeV it is apparent, see Fig. 6 (a), that Γ NLO (1, 1) deviates from the exact result by more than 20% for small values of tan β. Using instead the improved Higgs self-coupling λ Hhh Choosing a non-vanishing value of µ, see Fig. 6 (b) and (c) for µ = 400 GeV and µ = 1000 GeV, respectively, leads to deviations between the hMSSM prediction Γ NLO (1, 1) and the exact result Γ NLO eff . Such deviations increase with the value of µ and follow a similar pattern as observed for λ Hhh when compared against λ eff Hhh , see Fig. 4. Keep in mind, that the squared value of λ enters the LO partial decay width. Again the usage of λ Hhh , depicted in Fig. 6 (e) and (f), makes the overall spread of the discrepancies between the hMSSM prediction and the exact result smaller, but in particular for larger values of tan β they can still exceed 10 %.
Lastly, we suggest to calculate Γ NLO (1, 0) in the hMSSM approach, since this calculation of the decay width incorporates the most dominant one-loop corrections from top quarks, but on the other hand is not sensitive to the actual supersymmetric spectrum and is thus in the spirit of the hMSSM approach. This decay width is compared against the exact result including top-quark and stop corrections in Fig. 7. From the fact that for small µ values there is only a small difference to the results including one-loop stop corrections Γ NLO (1, 1), compare e.g. Fig. 6 (d) against Fig. 7 (a), we again conclude that momentum-dependent stop corrections are subdominant, at least for our benchmark scenario with stop masses above 1 TeV and a heavy Higgs mass m H well below the TeV scale. The comparison of Fig. 6 (e) and (f) against Fig. 7 (b) and (c), respectively, shows that also for larger µ the momentum-dependent corrections from the stop sector can well be neglected.

Conclusions
We revisited the calculation of the decay of the heavy Higgs boson H into two SM-like Higgs bosons h in the hMSSM approach. For this purpose we considered the full effective potential, in which the top quark and stops are integrated out, allowing to match the MSSM on a 2HDM as its low-energy limit. By carefully integrating out the top quark and stops separately, we identified missing contributions in the Higgs self-couplings of the hMSSM approach and suggested the definition of improved Higgs self-couplings for the hMSSM approximation. In particular in the limit of a small Higgsino mass parameter µ M S , which is an intrinsic assumption of the hMSSM approach, these improved Higgs self-couplings yield an excellent approximation for the actual Higgs self-couplings calculated in our effective potential approach. In this context the hMSSM approach can and has to be understood as an approximation to a full effective low-energy 2HDM matched to the MSSM. Since the effective potential is evaluated for zero external momenta, we included kinetic corrections at the level of the effective Lagrangian of the 2HDM involving dimension-4 operators. These corrections allow us to perform a oneloop Feynman-diagrammatic calculation of H → hh in the 2HDM-setup, where again we can consider the top quark and stops separately. By this procedure we discussed the relevance of momentum-dependent and kinetic corrections to the decay H → hh at the one-loop level, that were missing in previous comparisons. As for the Higgs self-couplings themselves using our improved Higgs self-couplings for the decay H → hh leads to a much better agreement between the calculation performed in the 2HDM-setup and the hMSSM approximation, in particular for small µ M S . We advocate to supplement the calculation of H → hh with momentum-dependent and kinetic corrections from the top quark to obtain a reliable prediction at one-loop level without introducing a dependence on the actual supersymmetric particle spectrum. Our study misses an inclusion of RGE effects, such that also heavier supersymmetric spectra can be investigated. We leave such comparisons to future work.