Reconciling EFT and hybrid calculations of the light MSSM Higgs-boson mass

Various methods are used in the literature for predicting the lightest CP-even Higgs boson mass in the Minimal Supersymmetric Standard Model (MSSM). Fixed-order diagrammatic calculations capture all effects at a given order and yield accurate results for scales of supersymmetric (SUSY) particles that are not separated too much from the weak scale. Effective field theory calculations allow a resummation of large logarithmic contributions up to all orders and therefore yield accurate results for a high SUSY scale. A hybrid approach, where both methods have been combined, is implemented in the computer code FeynHiggs. So far, however, at large scales sizeable differences have been observed between FeynHiggs and other pure EFT codes. In this work, the various approaches are analytically compared with each other in a simple scenario in which all SUSY mass scales are chosen to be equal to each other. Three main sources are identified that account for the major part of the observed differences. Firstly, it is shown that the scheme conversion of the input parameters that is commonly used for the comparison of fixed-order results is not adequate for the comparison of results containing a series of higher-order logarithms. Secondly, the treatment of higher-order terms arising from the determination of the Higgs propagator pole is addressed. Thirdly, the effect of different parametrizations in particular of the top Yukawa coupling in the non-logarithmic terms is investigated. Taking into account all of these effects, in the considered simple scenario very good agreement is found for scales above 1 TeV between the results obtained using the EFT approach and the hybrid approach of FeynHiggs.


Introduction
The properties of the Higgs boson that has been discovered by the ATLAS and CMS collaborations at the CERN Large Hadron Collider [1,2] are compatible with those predicted for the Higgs boson of the Standard Model (SM) at the present level of accuracy. Despite this apparent success of the SM, there are several open questions that cannot be answered by the SM and ask for extended or alternative theoretical concepts. Supersymmetry is one of best motivated frameworks for physics beyond the Standard Model (BSM), and in particular the Minimal Supersymmetric Standard Model (MSSM) is the most intensively studied scenario providing precise predictions for experimental phenomena in the LHC era.
Apart from associating a superpartner to each SM degree of freedom, the MSSM extends the Higgs sector of the SM by a second complex doublet. Consequently, the MSSM employs two Higgs-boson doublets, denoted by H 1 and H 2 , with hypercharges −1 and +1, respectively. After minimizing the scalar potential, the neutral components of H 1 and H 2 acquire vacuum expectation values (vevs), v 1 and v 2 . Without loss of generality, one can assume that the vevs are real and non-negative, yielding The two Higgs doublets in the MSSM accommodate five physical Higgs bosons. In lowest order these are the light and heavy CP-even Higgs bosons, h and H, the CP-odd Higgs boson, A, and two charged Higgs bosons, H ± . Two parameters are required to describe the Higgs sector at the tree level (conventionally chosen as tan β and the mass M A of the CP-odd Higgs particle); masses and couplings, however, are substantially affected by higher-order contributions. Until now, experiments have not found direct evidence for supersymmetric (SUSY) particles. On the other hand, precision observables provide an indirect access to the MSSM parameter space from which significant constraints on the allowed parameter regions can be obtained. On top of the classical set of electroweak precision observables, the mass of the detected Higgs boson constitutes an additional important precision observable, M exp h = 125.09 ± 0.24 GeV [3]. If the measured value is associated with the mass M h of the lightest CP-even Higgs boson within the MSSM (for a recent discussion of the viability of the interpretation in terms of the heavy CP-even Higgs boson H, see [4]), the comparison of the predicted value with the measurement constitutes an important test of the model with high sensitivity to the SUSY mass scales (see e.g. [5][6][7] for reviews). In order to fully exploit the high precision of the experimental measurement for constraining the SUSY parameter space the accuracy of the theoretical prediction for M h has to be improved very significantly.
So far, the full one-loop corrections [8][9][10][11], dominant two-loop corrections  and partial threeloop results [36][37][38] for the light MSSM Higgs-boson mass have been calculated diagrammatically. Besides fixed-order calculations, effective field theory (EFT) methods have been used to resum large logarithmic contributions in case of a large mass hierarchy between the electroweak and the SUSY scale [39][40][41][42][43]. These EFT calculations, however, are less accurate for relatively low SUSY mass scales owing to terms suppressed by the SUSY scale(s) which correspond to higher-dimensional operators in the EFT framework (see [44] for recent work in this direction).
In order to profit from the advantages of both methods -high accuracy for relatively low SUSY scales in the case of the diagrammatic approach versus high accuracy for a high SUSY scale in the case of the EFT approach -a hybrid method combining both approaches has been developed [45,46], see also [47,48] for different implementations. The method introduced in [45,46] has been implemented into the publicly available code FeynHiggs [11,17,[49][50][51] such that the fixed-order result is supplemented with higher-order logarithmic contributions.
Comparisons between FeynHiggs and pure EFT codes in the literature [43,47,48] have revealed nonnegligible differences between the predicted values for M h . In particular, deviations have been observed for large SUSY scales, where terms not captured in the EFT framework are supposed to be negligible. At first glance, such differences appear to be unexpected since the resummation of logarithms included in FeynHiggs is at the same level of accuracy as in pure EFT calculations.
In order to clarify the situation, it is the purpose of this work to perform an in-depth comparison of the various approaches to explain the origin of the observed differences. For simplicity, we choose a single-scale scenario, where M soft are the soft SUSY-breaking masses and µ is the Higgsino mass parameter. Furthermore, all parameters are assumed to be real, i.e. we work in the CP-conserving MSSM with real parameters. 1 While the chosen single-scale scenario is particularly suitable for the EFT approach, it should be noted that in realistic cases the actual task is to provide the most accurate prediction (together with a reliable estimate of the remaining theoretical uncertainties) for the Higgs-boson masses of the model for a given SUSY mass spectrum which may contain a variety of SUSY scales. We leave an investigation of such multi-scale scenarios for future work. We shall explain that there are essentially three sources of the observed differences. In a first step, we show that the usual scheme conversion of input parameters is not suitable for the comparison of results containing a series of higher-order logarithms. Such a scheme conversion can lead to large shifts corresponding to formally uncontrolled higher-order terms. Secondly, we analytically identify specific terms arising through the determination of the Higgs propagator pole which cancel with subloop renormalization contributions in the irreducible self-energies of the diagrammatic approach for a large SUSY scale. We develop an improved treatment where unwanted effects from incomplete cancellations are avoided. Thirdly, we show how different parametrizations of non-logarithmic terms can explain remaining differences between the results of FeynHiggs and pure EFT codes for high scales.
The paper is organized as follows. In Section 2, we review the different approaches with a particular focus on how the Higgs pole mass is extracted. In Section 3, we compare the results of the various approaches for the Higgs pole mass to each other. In Section 4, we discuss the issue of using DR input parameters as input of an OS calculation. In Section 5, we give a brief overview about the levels of accuracy of the M h evaluation implemented in various codes. In Section 6, we present a numerical analysis showing the impact of the effects discussed in the previous Sections and numerically compare FeynHiggs to other codes. The conclusions can be found in Section 7. Two appendices provide additional details.

Calculating the Higgs mass
In this Section, we shortly review how the pole mass of the lightest CP-even Higgs boson of the MSSM is calculated in a pure diagrammatic calculation, in a pure EFT calculation, and in the hybrid approach of FeynHiggs.

Diagrammatic fixed-order calculation
A well-established way to calculate corrections to the mass of the SM-like Higgs of the MSSM, as well as to the mass of the heavier CP-even neutral Higgs boson and the charged Higgs boson, is a fixed-order Feynman diagrammatic (FD) calculation. The prediction is based on the calculation of Higgs self-energies involving contributions from SM particles, extra Higgs bosons, as well as their corresponding superpartners. In this approach the contributions from all sectors of the model and of all particles in the loop can be incorporated at a given order. The mass effects of all particles in the loop can be taken into account for any pattern of the mass spectrum. If there is however a large splitting between the relevant scales, in particular a large mass hierarchy between the electroweak and the scale of some or all of the SUSY particles, the fixed-order result will contain numerically large logarithms that can spoil the convergence of the perturbative expansion.
In the MSSM with real parameters, after calculating the renormalized Higgs-boson self-energies, the physical masses of the CP-even Higgs bosons h, H can be obtained by finding the poles of their propagator matrix, whose inverse is given by where m h (m H ) denotes the tree-level mass of the h (H) boson andΣ hh,hH,HH are the corresponding selfenergies. We introduced the label "MSSM" to indicate that the corresponding self-energy contains SM-type contributions as well as non-SM contributions.
Concerning the renormalization, we follow here the approach used in the program FeynHiggs. Accordingly, the circumflexˆindicates that the self-energies have been renormalized using the mixed on-shell (OS) and DR-scheme of [11]. In particular, the A-boson mass is renormalized on-shell, whereas the Higgs field renormalization and the renormalization of tan β is performed using the DR scheme.
The masses of the weak gauge bosons (M Z , M W ) and the electromagnetic charge e are renormalized on-shell, and the tadpole renormalization is carried out such that the tadpole contributions are cancelled by their respective counterterms. The OS vev is a dependent quantity, which is given in terms of the OS values of the observables M W , s w and e by where s w denotes the sine of the weak mixing angle. The renormalization of this quantity at the one-loop level is therefore given in terms of the OS counterterms of M W , s w and e, where δM 2 W,Z are the mass counterterms of the W and Z bosons, respectively, and δe 2 is the counterterm of the electromagnetic charge (c 2 w = 1 − s 2 w ). Motivated by the fact that the renormalization of the vev receives a contribution from the field renormalization of the Higgs doublet, we identify the counterterm given in Eq. (5) with δv 2 OS /v 2 OS + δZ hh , where δZ hh is the field renormalization counterterm of the SM-like Higgs field fixed in the DR scheme. 2 Accordingly, the OS counterterm of the vev defined in this way reads The results for the self-energies in FeynHiggs have been reparametrized in terms of the Fermi constant G F instead of the electric charge e. The corresponding vev MSSM predictions for the quantity ∆r can be found in [52][53][54][55]. The effect of this reparametrization in the one-loop self-energies is formally of two-loop order. Furthermore (in the default choice), the stop sector is renormalized using the OS scheme, which is defined by applying on-shell conditions for the respective masses: the top-quark mass M t , and the top-squark masses Mt 1 and Mt 2 . A fourth renormalization condition fixes the mixing of the stops and can be identified with a condition for the top-squark mixing angle.
Employing this scheme, in FeynHiggs the full one-loop corrections to the Higgs self-energies as well as two- [11,17,22,23,25,27,30,33,34,[49][50][51]. While those two-loop corrections in the gauge-less limit have been obtained for vanishing external momentum, there is futhermore an option to incorporate the momentum dependence of the corrections at O(α t α s ) [56,57] (see also [58]). Finding the (complex) poles for the case where CP conservation is assumed corresponds to solving the equation In the decoupling limit, M A M Z , the physical mass of the lightest Higgs boson can approximately be obtained as solution of the simpler equation up to corrections from the hH and HH self-energies, which are suppressed by powers of M A . In the following discussion we will for simplicity use Eq. (9) for determining the pole of the propagator and we will furthermore neglect the imaginary parts of the self-energies. In FeynHiggs the complex poles of the propagator are obtained from the full propagator matrix, taking into account the real and imaginary parts of the Higgs-boson self-energies. Solving Eq. (9) iteratively for the case where imaginary parts are neglected yields an expression for the Higgs pole mass, where the prime denotes the derivative of the self-energy with respect to the momentum squared. The ellipsis stands for terms involving higher-order derivatives and products of differentiated self-energies. In App. B we provide a formula from which these terms can be derived recursively. The Higgs pole mass at a given order is obtained from Eq. (10) via a loop expansion to the appropriate order.

Effective Field Theory calculation
Another approach to calculate the mass of the SM-like Higgs boson in the MSSM is using effective field theory (EFT) methods. These allow the resummation of large logarithmic contributions, so that higherorder contributions beyond the order of fixed-order diagrammatic calculations can be incorporated. Without including higher dimensional operators in the effective Lagrangian, contributions suppressed by a heavy scale are however not captured.
with the MS vev (at the scale M t ). The MS vev can be related to the on-shell vev via the finite part of δv 2 OS defined in Eq. (6), It should be noted that since the quantity in Eq. (11) is the SM MS vev, in Eq. (12) only SM-type contributions have to be considered in δv 2 OS . Getting from the running mass (11) to the physical Higgs mass one has to solve the pole equation for the Higgs-boson propagator, involving the renormalized SM Higgs boson self-energy (denoted by a tilde) which is renormalized accordingly in the MS scheme at the scale M t but with the Higgs tadpoles renormalized to zero, i.e. the tadpole counterterm is chosen to cancel the sum of the tadpole diagrams, T SM h , for the Higgs field, With all these ingredients, the Higgs pole mass is now obtained as the solution of the equation Expanding the Higgs self-energy perturbatively around the tree-level mass m 2 h of the MSSM yields where the ellipsis indicates higher-order terms in the expansion. We discuss the current status of EFT calculations in Section 5.

Hybrid calculation
In FeynHiggs, the fixed-order approach is combined with the EFT approach in order to supplement the full diagrammatic result with leading higher-order contributions [45,46]. The logarithmic contributions resummed using the EFT approach are incorporated into Eq. (9), The quantity ∆Σ hh contains all logarithmic contributions obtained via the EFT approach as well as subtraction terms compensating the logarithmic terms already present in the diagrammatic fixed-order result forΣ MSSM hh , The subscript 'log' indicates that we take only logarithmic contributions into account. Note that in Σ MSSM hh (m 2 h ) the logarithms appear only explicitly when expanding in v/M SUSY . For more details on the combination of the fixed-order and the EFT result, we refer to [45,46].
Plugging the expression for ∆Σ hh into Eq. (18), we obtain for the physical Higgs mass We use the label 'nolog' to indicate that we take only terms not involving large logarithms into account for the labelled quantity. We again would like to stress that the large logarithms (and thereby the meant non-logarithmic terms) appear only explicitly inΣ MSSM hh (m 2 h ) when expanding in v/M SUSY . Before comparing the various approaches in depth, we also shortly comment on the renormalization scheme conversion needed for the combination of the fixed-order and the EFT calculation. As mentioned before, in FeynHiggs (in the default choice) the stop sector is renormalized using the OS scheme. In contrast, in the EFT calculation, i.e. the calculation of λ(M t ), all SUSY parameters enter in DR-renormalized form. As argued in [46], it is sufficient to convert only the stop mixing parameter X t using only the one-loop large logarithmic terms, where M 2 S = Mt 1 Mt 2 , α s = g 2 3 /(4π) (with g 3 being the strong gauge coupling) and α t = y 2 t /(4π) (with y t being the top Yukawa coupling).

Comparison of the different approaches
In the following we will discuss the differences between the various approaches. It is obvious from the discussion of the previous section that the diagrammatic fixed-order result and the pure EFT result differ by higher-order logarithmic terms that are contained in the EFT result but not in the diagrammatic fixed-order result as well as by non-logarihmic terms that are contained in the diagrammatic fixed-order result but not in the pure EFT result. In the hybrid approach the diagrammatic fixed-order result is supplemented by the higher-order logarithmic terms obtained by the EFT approach. We focus in the following on the comparison between the hybrid approach and the pure EFT result. In the present section we leave aside issues related to the used renormalization schemes, which will be addressed in Section 4.
While the hybrid approach and the pure EFT approach both incorporate the higher-order logarithmic terms obtained by the EFT approach, this does not necessarily imply that all logarithmic terms in the two results are the same. This is due to the fact that the determination of the Higgs-boson mass from the pole of the progagator within the hybrid approach is performed in the full model (in the example considered here the MSSM, incorporating loop contributions from all SUSY particles), while in the EFT approach it is determined in the effective low-scale model (in the considered example the SM). We will demonstrate below that the determination of the propagator pole in the hybrid approach generates logarithmic terms beyond the ones contained in the EFT approach at the two-loop level and beyond which actually cancel in the limit of a heavy SUSY scale with contributions from the subloop renormalization. This cancellation is explicitly demonstrated at the two-loop level. We will furthermore discuss the difference in non-logarithmic terms between the results of the hybrid and the EFT approach.

Higher-order logarithmic terms from the determination of the pole of the propagator
In the EFT approach where the Higgs boson mass is determined as the pole of the propagator in the SM as the effective low-scale model, while the SUSY particles have been integrated out, the logarithmic terms are given by (see Eq. (17)) The logarithmic terms contained in the result of the hybrid approach implemented in FeynHiggs are given by (see Eq. (20)) In the decoupling limit (M SUSY = M A M t , where in particular the light CP-even Higgs boson has SM-like couplings), we can split up the MSSM Higgs self-energy into a SM part and a non-SM part, In the mixed OS/DR scheme of the full diagrammatic calculation, the Higgs field renormalization constants are fixed in the DR scheme. For scalar propagators, there is no difference between the DR and the MS scheme at the one-loop level. Consequently, holds.
Using this relation, we obtain for the difference between the higher-order logarithmic terms from the determination of the pole of the propagator obtained in the EFT and the hybrid approach = : ∆ log p 2 . Since this difference, which is of two-loop order and beyond, results only from the momentum dependence of the non-SM contributions to the Higgs self-energy, we call it ∆ log p 2 in the following. We give analytic expressions for ∆ log p 2 in App. B. In Section 3.3 we will demonstrate at the two-loop level that in the limit of a heavy SUSY scale the quantity ∆ log p 2 consisting of "momentum-dependent non-SM contributions" as given in Eq. (26) cancels out with contributions of the Higgs self-energy's subloop renormalization. Before we address this issue we first compare the non-logarithmic terms in the two approaches.

Non-logarithmic terms
In the EFT approach, the non-logarithmic terms are given by (see Eq. (17)) By construction, all non-logarithmic terms contained in the result of the hybrid approach originate from the fixed-order diagrammatic calculation (see Eq. (20)), In this way one-and two-loop terms that are suppressed by the SUSY scale, ∆ nolog v/MSUSY , are included in the result of the hybrid approach. Terms of this kind would result from higher-dimensional operators in the EFT approach. Those terms that are included in the hybrid result as implemented in FeynHiggs but not in the publicly available pure EFT results constitute an important source of difference between the corresponding results, which is expected to be sizeable if some or all SUSY particles are relatively light (see also [44] for a recent discussion of contributions of this kind in the EFT approach). It should be noted that in general terms of O(v/M SUSY ) also originate from solving the full pole mass equation, Eq. (8), rather than the approximated one, Eq. (9).
At zeroth order in v/M SUSY , the non-logarithmic terms of the EFT approach contained in λ(M t ) in Eq. (27) At higher orders, all terms involving a derivative ofΣ nonSM hh are affected. As we will demonstrate in the following section, also the non-logarithmic non-SM contributions arising from the determination of the pole of the propagator cancel out with contributions of the subloop renormalization in the limit of a high SUSY scale.
Apart from these terms and from the non-logarithmic terms of O(v/M SUSY ) discussed above, ∆ nolog v/MSUSY , a further difference between the hybrid approach and the EFT aproach is due to the parametrization of the non-logarithmic terms. In the EFT approach all low-scale parameters are MS quantities. The results of FeynHiggs, on the other hand, are expressed in terms of physical, i.e. on-shell, parameters. For the top-quark mass both the results expressed in terms of the pole mass, M t , and the running mass at the scale M t , m t (M t ) (see [59] for details on the involved reparametrization) have been implemented (the applied renormalization schemes for SUSY parameters will be discussed below). The Higgs vev is a dependent quantity in FeynHiggs which is expressed in terms of the physical observables M W , s w and e according to Eq. (4) (where e is furthermore reparametrized in terms of the Fermi constant, see Eq. (7)). Accordingly, if choosing low-energy SM parameters to express the EFT result, the non-logarithmic terms in this result are parametrized in terms of the MS quantities m t (M t ) and v MS (M t ), while depending on the option chosen for the top-quark mass the non-logarithmic terms in FeynHiggs are expressed in terms of either m t (M t ) and v G F or M t and v G F . Those parametrizations differ from each other by higher-order terms. The observed differences are therefore related to the remaining uncertainties of unknown higher-order corrections.
It should be noted that also within the EFT approach there is a certain freedom for choosing different parametrizations. For instance, the threshold corrections at the matching scale can be expressed in terms of the SM MS top Yukawa coupling or in terms of the MSSM DR top Yukawa coupling.
As a result, the deviations ∆ nolog between the non-logarithmic terms in the hybrid approach and the EFT approach arise from the following sources, 4 Two-loop corrections controlled by the bottom and tau Yukawa couplings have recently been derived in [44].
Here ∆ nolog v/MSUSY are terms present in the hybrid approach that would correspond to higher-dimensional operators in the EFT approach. The term ∆ nolog para indicates the differences in the parametrization of the non-logarithmic terms, and ∆ nolog are terms arising from the different determination of the propagator poles, as discussed above.

Terms arising from the determination of the propagator pole at the twoloop level
We saw in Section 3.1 and Section 3.2 that the different determination of the propagator pole in the hybrid approach and the EFT approach gives rise to both logarithmic and non-logarithmic contributions in which the expressions given for the two approaches in the previous sections differ from each other. We will now explicitly demonstrate at the two-loop level that those differences in fact cancel out in the limit of a heavy SUSY scale if all the relevant terms at this order are taken into account. As a first step, we write down the correction to M 2 h , derived by an explicit diagrammatic calculation. At strict two-loop order, we obtain The superscripts indicate the loop-order of the corresponding self-energy. 5 We obtain the renormalized two-loop self-energy from the unrenormalized one viâ The subloop-renormalization can be derived from the one-loop self-energy via a counterterm-expansion. Expressing all couplings appearing in the one-loop self-energy through masses divided by v G F (for the remainder of this section we drop the subscript "G F ", i.e. we use the shorthand v ≡ v G F ), we can write where we used in the last line thatΣ We are interested in terms involving the finite parts of the derivative of the Higgs self-energy, i.e. terms which could potentially cancel the term proportional toΣ . At first sight it would seem that terms of this kind could arise from an on-shell field renormalization of the Higgs field. It is wellknown, however, that those field renormalization constants drop out of the prediction of the mass parameter order by order in perturbation theory (in FeynHiggs, a DR renormalization is employed for the Higgs fields). Also the mass counterterms as well as the genuine two-loop counterterms do not contribute terms that are 5 In our discussion here we treat the two-loop self-energy as the full result containg all contributions that appear at this order. The specific approximations that have been made at the two-loop level in FeynHiggs will be discussed below.
proportional toΣ nonSM,(1) hh (m 2 h ). The only remaining term is the vev counterterm. According to Eq. (6) and Eq. (7) it is given at the one-loop level by, having the same form in the SM and the MSSM, The renormalization constant δZ hh represents within the MSSM the DR field renormalization constant of the SM-like Higgs field, while in the SM it is understood to be the MS field renormalization constant of the Higgs field. We verified by explicit calculation that in the limit of a large SUSY scale the following relation holds Using this relation, we can rewrite the two-loop self-energy (omitting terms of O(v/M SUSY )), where the subscript '(δv 2 ) MSSM → (δv 2 ) SM ' is used to indicate that the MSSM vev counterterm, appearing in the subloop renormalization, is replaced by its SM counterpart. Plugging this expression back into Eq. (33) and Eq. (31), we obtain We observe that the corresponding subloop renormalization term cancels in Eq. (31) the termΣ nonSM,(1) hh (m 2 h ) involving the non-SM contributions to the Higgs self-energy by which the determination of the propagator pole in the hybrid approach differs from the EFT approach.
The origin of Eq. (35) is the different normalization of the SM-like MSSM Higgs doublet Φ MSSM and the SM Higgs doublet Φ SM . Comparing the derivative of the two-point function, appearing in the LSZ factor of amplitudes with external Higgs fields, 6 we obtain in the limit of a heavy SUSY scale or equivalently Expressed in terms of a relation between the counterterms of the vevs, this implies Eq. (35). While as mentioned above the Higgs field renormalization constant drops out in the Higgs mass prediction order by order, it is nevertheless noteworthy that the introduction of an OS field renormalization constant would lead toΣ and implying that no terms involvingΣ nonSM hh appear in the subloop renormalization at the two-loop level. While we have demonstrated this cancellation at the two-loop level, it is to be expected that it would also occur at higher orders. Explicit formulas for higher-order terms of this kind are given in App. B. While the described cancellation occurs at the full two-loop level, only partial cancellations occur between the full one-loop self-energy times its derivative and the two-loop self-energy if for the latter certain approximations are made.
In FeynHiggs, the two-loop self-energies are derived in the gaugeless limit (i.e., two-loop corrections of [22,23,25,27,33,34]), 7 and by default the external momentum of the two-loop graphs is neglected. There is, however, an option to include momentum dependence at O(α t α s ) (see [56,57]). Accordingly, all O(α 2 t , α t α b , α 2 b ) non-SM terms arising through the determination of the propagator pole at the two-loop level are cancelled in the limit of a large SUSY scale by corresponding subloop renormalization contributions within the diagrammatic calculation (the determination of the propagator pole obviously does not give rise to terms of O(α t α s , α b α s )). In previous versions of FeynHiggs, we have already taken care when constructing the subtraction terms according to Eq. (19) that we do not subtract logarithmic contributions that are needed for the cancellation with the corresponding terms arising from the determination of the propagator poles. For terms arising through the determination of the propagator pole beyond O(α 2 t , α t α b , α 2 b ), however, so far the cancellation in the limit of a large SUSY scale did not occur because the corresponding contributions in the irreducible self-energies at the two-loop level and beyond are not incorporated. In order to avoid unwanted effects from an incomplete cancellation, we have removed the uncompensated terms arising from the determination of the propagator pole in FeynHiggs.

DR parameters as input for an OS calculation
In this section we discuss issues related to the conversion between parameters of OS and DR renormalization schemes. While the discussion will focus on the case where DR input parameters are converted into OS ones that are then inserted into a result in the OS scheme, it should be stressed that the related problems are not intrinsic to the OS approach. The same problems would occur if a DR result were used with OS input parameters. The discussed problems are also not specific to Higgs mass predictions in SUSY models, but would appear whenever there are numerically large higher-order logarithms arising from a large splitting between the relevant scales of the considered quantity. In predictions for the mass of the SM-like Higgs boson within the MSSM, the result is however particularly sensitive to higher-order effects of this kind through the pronounced dependence on the stop mixing parameter X t , which receives large corrections when converting from the DR to the OS scheme or vice versa.
In the case where fixed-order results at the n-loop level obtained in two different renormalization schemes are compared with each other, and higher-order logarithms are unknown and not expected to be particularly enhanced, it is well known that the results based on the same type of corrections in two schemes differ by terms that are of O(n + 1). The same is true for different options regarding how to perform the parameter conversion that differ from each other by higher-order contributions. The numerical differences observed in such a comparison can therefore be used as an indication of the possible size of unknown higher-order corrections.
The situation is different, however, in the case that we are considering here, since the comparison is not performed between fixed-order results but between results incorporating a series of (resummed) higher-order logarithms. It is crucial in such a case that the correct form of the higher-order logarithms that can be derived via EFT methods, which in our case arise from the large splitting between the assumed SUSY scale and the weak scale, is maintained in the parameter conversion. We will demonstrate below that the parameter conversion that is usually applied for a comparison of renormalization schemes in fixed-order results does not maintain the correct form of the higher-order logarithms. Since those higher-order logarithms are numerically important, a conversion carried out in the described way leads to very large numerical discrepancies for large values of the SUSY scale.

Conversion between DR and OS parameters applicable to fixed-order results
The most straightforward method used for the conversion of DR input parameters to OS parameters in fixedorder results is to derive the shift between a parameter p in the two schemes according to p OS = p DR + ∆p at the considered loop order, see e.g. [60]. Accordingly, at the full one-loop level, including logarithmic as well as non-logarithmic terms, the conversion from DR to OS parameters for the stop mixing parameter and the stop masses, which are particularly relevant in the context of MSSM Higgs mass predictions, reads (for explicit formulas see [22,23,25,59]) Here ∆mt 1,2 is given by the corresponding difference of the DR and the OS counterterm. In FeynHiggs, the shift of X t is obtained by first calculating the OS stop masses and the OS stop mixing angle θ OS t . These are then used to obtain X OS Relating this prescription for X OS t to the DR input parameters X DR t , m DR t1 , m DR t2 , one can see that Eq. (45) contains products of one-loop contributions and therefore involves higher-order terms. Alternatively one could have used an expression for the conversion that is truncated at the one-loop level. The difference between the two prescriptions would be of the order of unknown higher-order corrections in a fixed-order comparison. The on-shell parameters obtained as described above are then used as input of the fixed-order OS renormalized calculation. This means in particular that the knowledge of the initial DR parameters is not used any further once the conversion to OS parameters has been carried out. While this procedure is suitable for fixed-order results, it leads to problems if results containing a series of higher-order logarithms are meant to be converted. Indeed, applying the described parameter conversion to the case of a DR result that incorporates higherorder logarithms generates additional higher-order terms causing a deviation in the logarithmic contributions. This can be seen by investigating the Higgs self-energy up to the two-loop level where the parameter X OS t obtained from the conversion has been inserted, Using instead Eq. (42) to write X OS t in terms of X DR t , and performing an expansion in ∆X t yieldŝ Thus, the obtained expression obviously differs from the original DR result by terms of 3-loop order and beyond. One would furthermore need to convert also all other parameters entering the self-energy to the DR scheme in order to exactly recover the DR renormalized self-energy.

The case of large higher-order logarithms
The higher-order terms in Eq. (49) that are not present in the original DR result contain in general logarithmic contributions which for a result containing a series of higher-order logarithms cause a deviation from the logarithmic corrections determined via the RGE. In our numerical discussion in Section 6 below we will demonstrate that those higher-order contributions that are induced by the parameter conversion are indeed numerically sizeable. Another issue that is relevant in a hybrid approach, as pursued in FeynHiggs, where a fixed-order result in the OS scheme is combined with higher-order logarithmic expressions that are expressed in the DR scheme concerns the DR value of X t that is used in the EFT part of the calculation. Only logarithmic terms are kept in the relation between X DR,EFT t and X OS t , see Eq. (21). If instead an input value for X DR t were converted to X OS t using the full one-loop contributions according to Eq. (42), the stop mixing parameter used in the EFT calculation of FeynHiggs, X DR,EFT t , would differ from the input parameter X DR t . In order to properly address the case where DR parameters associated with a result containing a series of higher-order logarithms are used as input for FeynHiggs, we follow the strategy to perform the parameter conversion in the fixed-order result rather than in the infinite series of higher-order logarithms. For this purpose we have extended FeynHiggs such that the incorporated fixed-order result is given in terms of the DR parameters X DR t , m DR t1 , m DR t2 (the actual input parameters are the soft-breaking parameters of the stop sector). This new result complements the existing result that is given in terms of the on-shell parameters The reparametrisation on which the new result is based can be viewed as the parameter conversion described in the example of the previous section, but truncated at the two-loop level,Σ We have used the same procedure as the one described here for the stop mixing parameter also for the stop masses. The two-loop terms that are induced by the conversion at the one-loop level have been added to the two-loop result derived in the on-shell scheme in order to arrive at the corresponding expression in the DR scheme. Compact expressions for these additional terms valid in the case M SUSY M t and degenerate M L = Mt R = M SUSY can be found in App. A. It should be noted that we would have obtained the same result if we had performed the diagrammatic calculation with a DR renormalization of the respective parameters instead of reparametrizing the final result. Using the above result given in terms of DR parameters, the value of X t that is used in the EFT part of the calculation equals the DR input parameter, X DR,EFT t = X DR t . For this setting in FeynHiggs with DR input parameters the subtraction terms have been adjusted such that the logarithms already contained in the fixed-order result for the DR renormalized self-energy are subtracted (rather than the ones contained in the OS renormalized self-energy, as it is the case for OS input parameters).
Accordingly, depending on the provided input parameters the evaluation of the prediction for the mass of the SM-like Higgs boson in FeynHiggs proceeds in the following ways: • For on-shell input parameters the on-shell fixed-order result is combined with the higher-order logs obtained in the EFT approach, where X DR,EFT t is related to X OS t as specified in Eq. (21).
• For DR input parameters in the stop sector associated with a result containing a series of higher-order logarithms the DR fixed-order result is combined with the higher-order logs obtained in the EFT approach, where X DR,EFT t = X DR t .
• For DR input parameters in a low-scale SUSY scenario where the impact of higher-order logarithms is expected to be small, both the fixed-order DR result and the fixed-order on-shell result can be employed, where for the latter the parameter conversion described in the previous section is used.

Comparison of FeynHiggs to other codes
In the previous sections, we investigated methodical differences between the different approaches for predicting the lightest CP-even Higgs boson mass in the MSSM, focusing in particular on the comparison of the hybrid approach implemented in FeynHiggs with a pure EFT calculation. In the following, we compare FeynHiggs numerically to other codes. Publicly available codes based on diagrammatic fixed-order results or effective potential methods include CPSuperH [61][62][63], SoftSUSY [64], SPheno [65,66] and SUSPECT [67]. Publicly available pure EFT calculations are SUSYHD [43] and MhEFT [40,42,68]. FlexibleSUSY [69], based on SARAH [70][71][72][73], includes both a diagrammatic and an EFT result. Furthermore, it also has the option to use a hybrid method different from the one pursued in FeynHiggs, called FlexibleEFTHiggs [47]. Its basic idea is to include terms suppressed by the SUSY scale into the matching conditions in order to obtain accurate results for both low and high scales. Recently, the same approach has been included into SPheno [48].
The different levels of higher-order corrections implemented in the various diagrammatic codes are listed in [74]. A detailed numerical comparison between various diagrammatic and EFT codes can be found in [47]. In there, it is also discussed in detail how FlexibleEFTHiggs compares to other codes. We therefore focus in this work on a comparison of FeynHiggs to SUSYHD as an exemplary EFT calculation.
Before we can investigate the impact of the effects discussed in the previous Sections on the comparison of FeynHiggs and SUSYHD, we have to ensure that the RGE results, i.e. the results for λ(M t ), of both codes are compatible with each other. Both codes implement full leading and next-to-leading resummation and O(α s α t , α 2 t ) next-to-next-to-leading resummation of large logarithms. So the levels of accuracy are basically identical. There are however several differences which are listed below.
• SUSYHD by default uses the top-Yukawa coupling extracted at the NNNLO level. FeynHiggs instead uses the NNLO value by default, which is formally the appropriate setting for the resummation of NNLL contributions. For all numerical results shown in this work, we deactivate the NNNLO corrections to the top-Yukawa coupling in SUSYHD.
• SUSYHD includes the bottom-and tau-Yukawa couplings in the renormalization group running and also includes corresponding one-loop threshold corrections. In FeynHiggs, the bottom and tau Yukawa couplings are set to zero in the EFT calculation. In the fixed-order diagrammatic calculation, however, terms proportional to the bottom Yukawa coupling are included at the one-and two-loop level (at the one-loop level for the case of the tau Yukawa coupling).
• SUSYHD includes the electroweak gauge couplings in the running up to the three-loop level. FeynHiggs takes them into account up to the two-loop level. At the three-loop level, they are set to zero.
• FeynHiggs includes a one-loop running of tan β to relate tan β(M t ), which is used as input of FeynHiggs, to tan β(M SUSY ), which enters through the matching at the SUSY scale. In contrast, SUSYHD uses tan β(M SUSY ) as input.
• Similarly, FeynHiggs uses a DR renormalized Higgsino mass parameter µ at the scale M t . The running to the scale M SUSY , at which it enters the EFT calculation via the matching conditions at the SUSY scale, is neglected. SUSYHD uses µ(M SUSY ) as input.
More details on the implemented EFT calculations are given in [43,46]. Despite the listed differences including the different treatment of the renormalization scales of tan β and µ, we find excellent agreement between the results of the RGE running of both codes. The numerical difference of the quantity v 2 λ(M t ) calculated using the two codes is always 50 GeV 2 for the single scale scenario defined in Eq. (2) and tan β ∼ O(10). This translates into a difference in M h of 0.1 GeV.

Numerical results
In this Section, we present a numerical investigation of the effects discussed in the previous Sections and compare the result obtained by FeynHiggs to SUSYHD as an exemplary pure EFT code. We restrict ourselves to the single scale scenario defined in Eq. (2). Apart from the parameters of the stop sector, we neglect all renormalization scheme conversions necessary to relate the parameters of Eq. (2) as defined in FeynHiggs to the parameters as defined in SUSYHD. We furthermore set tan β = 10, i.e. we use tan β(M t ) = 10 as input for FeynHiggs and tan β(M SUSY ) = 10 as input for SUSYHD. As mentioned in Section 5, the difference in the renormalization scales is negligible for the considered scenario. All soft-breaking trilinear couplings except the one of the scalar top quarks are choosen to be A e,µ,τ,u,d,c,s,b = 0.
For all soft-breaking parameters (i.e. those of the stop sector), we use the DR scheme with the renormalization scale being M SUSY . If not stated otherwise, we use a parametrization of the non-logarithmic contributions in terms of the SM MS NNLO top mass and v G F (see Section 3.2), corresponding to choosing runningMT = 1 as FeynHiggs flag. We first look at the numerical difference between employing the type of conversion from DR to OS input parameters which is suitable for the comparison of fixed-order results ("FH 2.13.0 param conv") and using a DR renormalized fixed-order result ("FH 2.13.0 DR"), see the discussion in Section 4. The left plot of Fig. 1 shows the corresponding results for X DR t /M SUSY = 0 (2) as solid (dashed) lines as a function of M SUSY . One can see that for M SUSY 5 TeV the difference between the two methods leads to an approximately constant shift in the prediction for M h . For vanishing mixing the prediction obtained by using a DR renormalized fixed-order result is ∼ 0.5 GeV higher than the one obtained by a naive scheme conversion of the input parameters. For X t /M SUSY = 2, the shift is larger. The result obtained using a DR fixed-order result is ∼ 1 − 1.5 GeV smaller than the one obtained by the naive conversion of the input parameters. The shifts occur not only for scales of a few TeV, but also for very low scales (M SUSY 0.3 TeV). Therefore, we conclude that at low scales the observed shifts are mainly caused by non-logarithmic higher-order terms by which the DR result and the result involving a parameter conversion differ from each other. As usual, non-logarithmic terms tend to be larger for |X DR t /M SUSY | ∼ 2 than for vanishing stop mixing. For M SUSY 5 TeV, we observe that the difference be-tween the two results is increasing rapidly to up to 10 GeV for vanishing mixing and up to 5 GeV for |X DR t /M SUSY | ∼ 2 in the region up to M SUSY ∼ 20 TeV. This behavior is mainly due to the fact that the parameter conversion that is used for the comparison of fixedorder results induces higher-order logarithmic contributions that are not compatible with the implemented resummation of logarithms to all orders (see the discussion in Section 4.1). For high SUSY scales, where the higher-order logarithmic contributions become numerically large, this mismatch leads to the observed large deviations. To a lesser extent, also the deviation between the input X DR t and the X DR,EFT t used in the EFT calculation plays a role in this context, see Section 4.2.
In the right plot of Fig. 1 the two results are compared as a function of X DR t /M SUSY for M SUSY = 1, 5, 20 TeV, shown as solid, dashed and dot-dashed lines, respectively. For M SUSY = 1 TeV and M SUSY = 5 TeV the deviations stay relatively small except for the highest values of |X DR t /M SUSY |. In contrast, for M SUSY = 20 TeV the uncontrolled higher-order contributions induced by the naive conversion of the input parameters are seen to have a huge effect which even reverts the usual pattern of the dependence on |X DR t /M SUSY |, giving rise to local minima at |X DR t /M SUSY | ±2.3. We emphasize again that the same kind of uncontrolled higher-order effects would occur if a naive conversion of OS to DR parameters would be used as input for a DR result containing a series of numerically large higher-order logarithms. Fig. 1 shows that numerical instabilities noticed in comparisons of EFT results with FeynHiggs carried out in the literature are a consequence of an inappropriate application of the conversion of input parameters between the OS and the DR schemes. The higher-order contributions implemented in FeynHiggs are seen to be numerically stable up to very high SUSY scales in the considered scenario.
For the further FeynHiggs results shown below we use the DR renormalization of the parameters in the stop sector. As a next step we investigate the impact of the terms arising from the determination of the propagator pole. As explained in Section 3, there occurs a cancellation in the limit of a large SUSY scale between non-SM terms arising through the determination of the propagator pole and contributions from the subloop renormalization of the irreducible self-energy diagrams. While up to the version FeynHiggs2.13.0 this cancellation was incomplete for terms beyond O(α 2 t , α t α b , α 2 b ) (see Eq. (20)), we have modified the determination of the propagator poles in the new version of FeynHiggs such that terms are omitted that would not cancel because their counterpart in the irreducible self-energies is not incorporated at present. In Fig. 2 FeynHiggs2.13.0 DR is compared with the new version, which is labelled as FeynHiggsnew DR. The difference between the two results corresponds essentially to the terms ∆ logs p 2 and ∆ nolog p 2 given in Eqs. (26) and (29). In the left plot of Fig. 2, we show the results as a function of M SUSY for X DR t = 0 and X DR t /M SUSY = 2. One observes that the difference grows nearly logarithmically with M SUSY . This is expected since the Having investigated the numerical impact of the scheme conversion of the input parameters as well as of the terms arising from the determination of the propagator pole, we now turn to a direct comparison of FeynHiggs with SUSYHD. 8 The FeynHiggs results in this comparison are the ones of the new version, FeynHiggsnew DR, where the stop sector is renormalized in the DR scheme and terms arising from the determination of the propagator pole are omitted that go beyond the level of the corrections implemented in the irreducible self-energies, as described above.
The left plot of Fig. 3 shows M h as a function of M SUSY for X DR t /M SUSY = 0 (2) as solid (dashed) lines. For vanishing stop mixing and M SUSY 1 TeV, we observe an excellent agreement of the FeynHiggs curve with the SUSYHD result. Even for very large scales M SUSY 20 TeV, we find agreement within ∼ 0.5 GeV in the considered simple numerical scenario, in which all SUSY scales are chosen to be equal to each other. For low scales (M SUSY 1 GeV), it can be seen that the FeynHiggs result is higher by up to ∼ 1.7 GeV compared to the SUSYHD result. The origin of this difference are terms suppressed by the SUSY scale, which are included in FeynHiggs but not in SUSYHD, as will be discussed below. For X DR t /M SUSY = 2, we basically observe the same behavior as in case of vanishing stop mixing. The overall agreement in the simple numerical scenario is very good (within ∼ 0.7 GeV for M SUSY 0.5 TeV). For low scales (M SUSY 0.5 GeV), the FeynHiggs result is lower compared to the SUSYHD result by up to ∼ 1 GeV. As in the case of vanishing stop mixing, this can be traced back to terms suppressed by the SUSY scale. We will discuss this and investigate the remaining differences in more detail below.
In the right plot of Fig. 3 3 TeV), the difference is only slowly increasing when raising M SUSY . For vanishing stop mixing, the difference is growing by ∼ 50 GeV 2 when raising M SUSY from 3 TeV to 20 TeV. For X DR t /M SUSY = 2, similarly a growth of ∼ 50 GeV 2 is recognizable. This behavior is mostly due to the differences in the EFT calculations implemented in FeynHiggs and SUSYHD discussed in Section 5. In addition, however, we observe an offset relative to the zero axis for M SUSY 3 TeV. For vanishing stop mixing, it is small (∼ 50 GeV 2 ), whereas for X DR t /M SUSY = 2, the shift is more significant (∼ 150 GeV 2 ). The nearly constant offset between the two codes can be traced back to the different parametrization of the non-logarithmic terms discussed in Section 3.2.
We further analyse the influence of the different ways to parametrize the non-logarithmic terms in the right plot of Fig. 4. It shows the difference in M 2 h obtained from a diagrammatic calculation of O(α t α s , α 2 t ) using different parametrizations of the vev for the non-logarithmic one-and two-loop terms (see Section 3.2 for more details). Note that these non-logarithmic terms, apart of O(v/M SUSY ) contributions, stay constant when varying M SUSY . For X DR t /M SUSY ∼ 2 the difference between parametrizations in terms of v G F and v MS (both using the SM NNLO MS top-quark mass) amounts to ∼ 170 GeV 2 . Such a shift accounts for the main part of the nearly constant offset observed in the left plot of Fig. 4. For X DR t /M SUSY ∼ 0 the difference between the parametrizations in terms of v G F and v MS is seen to become very small. The nearly constant offset for vanishing stop mixing observed in the left plot of Fig. 4 can be explained in a similar way by different parameterization of terms that are not of O(α t α s , α 2 t ). Finally, we briefly comment on the differences between FeynHiggs and other codes that have been reported in the literature. In [43] it was claimed that differences between FeynHiggs and SUSYHD of up to ∼ 9 GeV would occur for M SUSY = 2 TeV and X DR t /M SUSY ∼ √ 6. As already noted in [43], this difference was somewhat reduced if the NNLO MS top mass was employed in the calculation of FeynHiggs. 9 While at the time of the comparison carried out in [43] the EFT calculation of FeynHiggs was not yet at the same level of accuracy as the one of SUSYHD, the differences claimed in [43] were in fact primarily caused by an inappropriate application of the conversion of input parameters between the DR and the OS scheme. The inappropriate parameter conversion, for which the authors of [43] used their own routine, caused a deviation of 3-4 GeV for M SUSY = 2 TeV and X DR t /M SUSY ∼ √ 6 and was also responsible for the apparent numerical instability at large SUSY scales of the FeynHiggs curve with X DR t /M SUSY = 0 shown in [43]. The numerical effect of this deviation was larger than the shift caused by employing the NNLO or NNNLO MS top-quark mass in FeynHiggs, in contrast to the claim made in [43].
Also the comparison figures shown in [47,48] are plagued by deficiencies arising from an inappropriate application of the parameter conversion between the DR and the OS scheme. We stress again that such a parameter conversion would give rise to the same kind of problems when starting from OS parameters and converting to DR ones.

Conclusions
We have presented a detailed comparison between various approaches used to predict the mass of the SM-like Higgs boson in the MSSM in a scenario in which all SUSY mass scales are chosen equal to each other. In particular we have compared pure EFT calculations with the hybrid approach, in which an explicit Feynmandiagrammatic fixed-order result is combined with the leading higher-order contributions obtained from EFT methods. In the literature significant deviations between the results obtained via the two approaches have been reported especially at large SUSY scales. In this work, we have identified three sources of the observed differences.
We could show that a large part of the reported discrepancies can be traced back to parameter conversions between different renormalization schemes. In EFT calculations typically the DR scheme is used for the renormalization of SUSY breaking parameters, e.g. the stop mixing parameter. In the diagrammatic calculation of FeynHiggs (in the default case) however, the OS scheme is employed in the scalar top sector. We have demonstrated that the usual scheme conversion of input parameters used for the comparison of fixed-order results is not suitable for the comparison of results containing a series of higher-order logarithms. This kind of parameter conversion would induce higher-order logarithmic contributions that are not compatible with the implemented resummation of logarithms to all orders. We have shown that the form of the higher-order logarithms obtained in one scheme can manifestly be maintained if the fixed-order part of the calculation is consistently reparametrized to this scheme. In order to enable this approach for DR input parameters, we have extended FeynHiggs such that the results are provided both in terms of the on-shell parameters X OS t , Mt 1 ≡ m OS t1 , Mt 2 ≡ m OS t2 (as before) and the DR parameters X DR t , m DR t1 , m DR t2 . In practice, this was achieved by reparametrizing the existing OS fixed-order result. We have demonstrated that many of the apparent discrepancies reported in the literature have mainly been caused by an inappropriate application of the conversion of input parameters between the OS and the DR schemes. It should be emphasized that this issue is not a problem of the OS renormalization, but analogously appears if OS parameters are used as input for codes employing the DR scheme.
Another difference between pure EFT calculations and the hybrid approach arises from the determination of the poles of the Higgs propagator matrix. We have shown explicitly at the two-loop level that there occurs a cancellation in the limit of a large SUSY scale between non-SM terms arising through the determination of the propagator pole and contributions from the subloop renormalization of the irreducible self-energy diagrams. Since we expect that similar cancellations will happen at higher loops, we have modified the determination of the propagator poles in the new version of FeynHiggs such that terms are omitted that would not cancel because their counterpart in the irreducible self-energies is not incorporated at present. Unless otherwise stated, the numerical results presented in this paper have been obtained with this new version of FeynHiggs. Numerically, we found that the terms beyond O(α 2 t , α t α b ) for which in previous versions of FeynHiggs the cancellation was incomplete are negligible for low scales (M SUSY 0.5 TeV). They can be more significant for high scales (∼ 1.5 GeV for M SUSY ∼ 20 TeV). Furthermore, we investigated the impact of different parametrizations of the non-logarithmic one-and two-loop terms. In this context, we found the top-quark quark mass and the vev to be especially relevant. Despite the results being formally identical at the strict two-loop level, using e.g. a SM NNLO MS top quark mass instead of the OS top quark mass induces changes in the higher-order non-logarithmic contributions.
In our numerical comparison, we focused on a single scale scenario with a moderate tan β, which is particularly suited for an EFT calculation. We specifically compared the results of FeynHiggs and the EFT code SUSYHD. Using the NNLO value of the MS top Yukawa coupling in SUSYHD (by default the NNNLO value is used in SUSYHD, which leads to a downward shift by ∼ 0.5 GeV in M h ), we find very good agreement between the new version of FeynHiggs and SUSYHD for scales M SUSY 1 TeV. Such a good agreement is in fact expected for high SUSY scales since the hybrid approach of FeynHiggs incorporates essentially the same logarithmic contributions as pure EFT calculations. For M SUSY 1 TeV we observe significant differences between FeynHiggs and SUSYHD due to terms suppressed by the SUSY scale that are not incorporated in the EFT calculation of SUSYHD. The observed differences stay relatively small for the considered simple scenario with a single SUSY scale, reaching ∼ 1 GeV for M SUSY ∼ 300 GeV. Larger deviations can be expected in SUSY scenarios with non-negligible mass splittings between the various SUSY particles. Such kind of mass patterns are accounted for in the diagrammatic fixed-order part of the hybrid approach.
The new version of FeynHiggs described in this paper, comprising an improvement in the determination of the propagator poles and an option for using the DR scheme for the renormalization of the stop sector, will be made public soon.
The results obtained in this paper provide important input for an improved estimate of the remaining theoretical uncertainties from unknown higher-order corrections. In this context, we would like to stress once more that for the numerical evaluations in this paper we have used a rather simple scenario where all SUSY masses have been set to be equal to each other. Having reconciled the hybrid approach of FeynHiggs with pure EFT calculations for this simple single scale scenario, we are now in a position to assess the accuracy of the theoretical predictions also for more general scenarios with different hierarchies of scales. This will be analysed in a fothcoming publication. drop the subscript "G F ", i.e. we use the shorthand v ≡ v G F ) where an n-tuple of non negative integers (a 1 , ..., a n ) is an element of T n if 1 · a 1 + 2 · a 2 + ... + n · a n = n.
The zeroth order correction serves as starting point of the recursion. We split ∆ logs p 2 into a leading, a next-to-leading and a next-to-next-to-leading logarithm piece, In FeynHiggs, the full momentum dependence by default is taken into account only at the one-loop level. At the two-loop level, the external momentum is set to zero (see [56,57] for a discussion of the momentum dependence at the two-loop level). We can therefore split up the non-SM contributions to the Higgs selfenergy into a one-and a two-loop piece, To shorten the expressions for the individual contributions, we first introduce abbreviations. We write the non-SM contributions to the Higgs self-energy aŝ Σ nonSM,(1) hh (m 2 h ) = k c χ 1,1 L χ + c A 1,1 L A + cf 1,1 L S + c 1,0 , Σ nonSM,(2) hh (0) = k 2 c 2,2 L 2 S + c 2,1 L S + c 2,0 , where k ≡ (4π) −2 is used to keep track of the loop order and Here it should be noted that in this work we set In this Appendix, however, we keep them independent to be able to use the expressions also for more general cases.
The subscript of a coefficient c a,b indicates that it is the prefactor of the term k a L b (L = L χ , L A , L S ). The corresponding superscript marks the origin of the respective term (from EWinos χ, from heavy Higgses A or from sfermionsf ). These superscripts are used only at the one-loop level to be able to differentiate between the different types of appearing logarithms (L χ , L A and L S ). In the DR scheme, the appearing coefficients up to O(v 2 /M 2 heavy ) (M heavy = M χ , M A , M SUSY ) are given by (for the remainder of this section we drop the subscript "G F ", i.e. we use the shorthand v ≡ v G F ) where all terms in the c coefficients proportional to g or g are set to zero. Correspondingly, the derivatives of the light self-energy only include terms proportional to y t . These are listed in Eqs. (79)-(83). This loop expansion quickly converges such that we can safely drop higher-order contributions (8-loop and beyond).
We find the electroweak contributions at the NNLL level and even higher-order logarithms (N n L with n > 2) to be completely negligible. Similar expressions can easily be obtained for the non-logarithmic terms of the same origin (see Eq. (30)).