Fighting off field dependence in MSSM Higgs-mass corrections of order $\alpha_t\,\alpha_s$ and $\alpha_t^2$

The connection between gauge and Higgs sectors makes supersymmetric extensions of the Standard Model predictive frameworks for the derivation of Higgs masses. In this paper, we study the contamination of such predictions by field-renormalization constants, in the MSSM with two-loop gaugeless corrections of $\mathcal{O}{\left(\alpha_{t,b}\,\alpha_s,\,\alpha_{t,b}^2\right)}$ and full momentum dependence, and demonstrate how strict perturbative expansions allow to systematically neutralize the dependence on such unphysical objects. On the other hand, the popular procedure consisting in an iterative pole search remains explicitly dependent on field counterterms. We then analyze the magnitude of the intrinsic uncertainty that this feature implies for the iterative method, both in non-degenerate and near-degenerate regimes, and conclude that this strategy does not improve on the predictions of the more straightforward expansion. We also discuss several features related to the inclusion of the orders $\alpha_{t,b}\,\alpha_s$ and $\alpha_{t,b}^2$ in the so-called 'fixed-order' approach, such as the resummation of UV-logarithms for heavy supersymmetric spectra.


Introduction
The discovery of a Standard-Model (SM)-like Higgs boson at the LHC [1][2][3] and the growingly precise measurements of its characteristics [4][5][6] make the need for controlled theoretical predictions in the Higgs sector of models of new physics quite clear, if one wishes to exploit associated observables and constrain the parameter space of beyond-the-SM (BSM) theories. In many cases, such precision calculations are difficult to conduct in the full model, e. g. due to the nonperturbative nature of the theory, and can only be performed in an effective context, in the limit of decouplingly heavy new physics. While the absence of conclusive signs of BSM physics at the LHC seems to justify the relevance of this type of approach anyway, the case of supersymmetric (SUSY) extensions of the SM [7,8] is still remarkable in that it provides a perturbative framework up to relatively high energies.
In Ref. [66], we argued that terms violating the EW symmetry in an uncontrolled fashion were introduced in the predictions of masses and decays at the 1L order unless one carefully expands and truncates the various contributions entering the amplitudes. In particular, momentumdependent corrections can lead to an explicit dependence on field-renormalization constants, thus producing UV-logarithms that have no quantitative meaning. While such artificial effects are formally of a higher order compared to that of the explicitly included self-energies, it is important to assess their numerical impact in order to estimate the actual gain in predictivity when incorporating radiative corrections in a fashion that explicitly depends on the regularization procedure. 2 Alternatively, a strict expansion and truncation order-after-order can limit the dependence on the field and gauge-fixing regulators and provide results where spurious symmetry-violating effects are systematically neutralized. The purpose of the present paper is to explain how such a formalism can be extended to address 2L mass corrections. We focus on the orders α q α s and α 2 q , since these are fully exploitable at 2L, while the generic EW corrections miss some contributions needed for a direct connection to observable quantities. We build upon the existing results of Refs. [13][14][15][16], but include full momentum dependence now, also making use of the relations derived in Ref. [31]. Below, we restrict ourselves to the simplest model where the considered orders matter, i. e. the Minimal Supersymmetric Standard Model (MSSM), although the method straightforwardly applies to extensions, e. g. the Next-to-MSSM (NMSSM)-yet more care is then needed in order to consistently account for the gaugeless limit and process the more motivated scenarios with large radiative Higgs mixing.
In parallel with the computations at FO, where self-energies are derived order after order in the full model considered in a given renormalization scheme, the mass-gap suggested by the absence of discovery of new-physics particles has encouraged the discussion of SUSY Higgs sectors in the context of EFTs, allowing for a resummation of the large UV-logarithms that develop between the SUSY and EW scales; see e. g. the recent work of Refs. [32][33][34][35][36][37][38][39][40][41][42][43][44][45] and a more complete list of references in the review [9]. While the effects that we discussed in Ref. [66] also formally apply in such EFTs, the weight of regulators in the definition of mass observables is less critical than in FO calculations. Indeed, the direct variation of the gauge-fixing parameters or the field counterterms in the EFT is necessarily reduced in such a formalism where large logarithms are forced to follow an SU (2) × U (1)-conserving pattern. On the other hand, such a variation represents only a partial picture in testing the connection between observables in the EFT context, as the MSSM-EFT matching procedure then plays a determining part in establishing this relation. This connection between the predicted Higgs masses and observable input is however at the center of what we probe through variations of regulators in the FO approach. We will not attempt to address the question of this relation in EFTs here. 1 By convention, we denote the considered classes of Feynman diagrams by the type of the involved couplings α f = Y 2 f (4 π), α s = g 2 s (4 π), and α = g 2 EW (4 π), where Y f is the Yukawa coupling of fermion f , g s is the strong gauge coupling, and g EW is any of the electroweak gauge couplings. In addition, we introduce α q ≡ α t,b when considering corrections controlled by the Yukawa couplings of both heavy quarks. 2 Here, we employ the term 'regularization' in reference to the procedure giving a meaning to ill-defined selfenergies away from their mass shell, not in reference to the counting of UV-divergences in Feynman amplitudes. Field-renormalization constants are the usual regulators-and the only ones that we consider in this paperallowing such an extrapolation, and the ambiguity rests with their arbitrary finite part, from which actual observables are supposedly independent.
In Sect. 2, we analyze how to avoid or minimize the dependence of 2L corrections to the Higgs masses on the choice of field renormalization, taking this criterion as our guideline for the inclusion of the 2L contributions of O α t α s and O α 2 t . We then numerically compare the corresponding mass predictions with those obtained through an iterative pole search, both in non-degenerate and near-degenerate scenarios, checking how variations of the field counterterms affect each in view of the magnitude of the 2L effects. We also probe how the effective-potential approximation performs for SM and non-SM Higgs states. This forms the content of Sect. 3 and Sect. 4. In Sect. 5, we perform a resummation of UV-logarithms of O α q , α q α s , α 2 q in the FO context, as the increasing weight of these effects for large SUSY scales otherwise limits the applicability of the FO approach. A brief summary is provided in Sect. 6.

Inclusion of 2L corrections to the Higgs mass observables
In this section, we present a brief description of the computation of Higgs masses at higher order, with reference to more detailed derivations in the appendix. We explicitly extract the dependence of the radiative contributions on Higgs-field counterterms and analyze the conditions for its cancellation at the level of Higgs-mass predictions at the 2L order.

Invariance under field-renormalization
As usual in Quantum Field Theory, the mass observables are calculated from the two-point truncated and connected correlators of the model. The complex zeroes of the characteristic polynomial in the Higgs sector should indeed correspond to the poles describing the Higgs resonances in particle scattering: a derivation is proposed in Appx. A. The defining equation for these poles M 2 k thus reads: where M 2 tree = diag m 2 i represents the tree-level mass matrix andΣ p 2 is the renormalized self-energy matrix for external momentum p.
A problematic feature in Eq. (1) is that the renormalized self-energiesΣ ij p 2 are not welldefined away from their mass-shell at p 2 = m 2 i + m 2 j 2. Indeed, absorption of the UV-divergences in such objects requires the introduction of field-renormalization constants, where the upper indices in parentheses refer to the loop order of the renormalization constant. However, since quantum fields are not measurable as such, these field counterterms should be strict bookkeeping quantities, which drop out at the level of observables. An iterative resolution of Eq. (1), as often advocated in the literature [17-19, 56, 57, 67-73], leads to an explicit dependence of the predicted masses on the field counterterms. This dependence should ideally become negligibly small when sufficiently high orders are included in the calculation and, in the meanwhile, the variation of the field counterterms could be seen as setting a lower bound on the higher-order uncertainty, by weighing such terms of higher order that are introduced in the pole search.
In Ref. [66], we argued against this iterative approach, in particular because it generates symmetry-violating terms that are not controlled by the vacuum-expectation values (v.e.v.-s) of the Higgs fields, but appear as pure artifacts of the formalism. Consequently, these partial higher-order, gauge-dependent terms introduced by the pole search do not represent a genuine estimate of higher-order corrections but amount to an intrinsic error of the procedure. Instead, we preferred to employ an explicit expansion and truncation where the dependence on field counterterms and gauge-fixing parameters vanishes or is minimized order-by-order. At 1L, this procedure leads us to systematically use diagonal self-energies evaluated on their mass-shell, while off-diagonal self-energies can be dismissed as contributions of higher order. Only in neardegenerate scenarios are the replacementΣ ij 2 and the inclusion of off-diagonal self-energies justified by the counting m 2 i − m 2 j = O(1L). Below, we discuss how to extend this procedure to include 2L mass corrections. As the considered orders do not involve EW gauge corrections, we must take solely independence from field-renormalization constants as our guiding criterion. This condition can be formulated in two different fashions: 1. the observables should be invariant under any choice of field counterterm; 2. the observables are UV-finite without need of considering field counterterms.
The second prescription proves to be a weaker requirement than the first one, in particular because the UV-divergent part of dΣ ij dp 2 is symmetric under the exchange i ↔ j, so that amplitudes may in general depend on antisymmetric finite contributions to the field counterterms while these are not needed to achieve UV-finiteness. In practice, the first criterion is satisfied in the non-degenerate case, while we must content ourselves with the second one in the neardegenerate scenario-at least in the strategy presented in Ref. [66]. We will look upon the second condition as being roughly equivalent to an invariance under symmetric choices of the field counterterms-in fact a more constraining requirement.

Non-degenerate case
For a non-degenerate state with tree-level mass m i , it is possible to extract the nearby pole M 2 i in the propagator matrix by solving the recursive Eq. (1); at the 2L order this amounts to solving M 2 i is the quantity defined in Eq. (40)-refer e. g. to the derivation in Appx. A.2. The right-hand side of this equation can be further expanded using the condition in the argument of the renormalized self-energies (Σ (k) represents the renormalized self-energy of k-loop order): The leading correctionΣ (1) ii M 2 i has been expanded so that it gives a contribution that is independent from field renormalization (and gauge-fixing) at strict O(1L).
Similarly, invariance under variations of the 2L field counterterm dictates the choice p 2 ! = m 2 i as the argument ofΣ (2) ii in Eq. (3). In addition, restriction to the orders α q α s and α 2 q entails a specific distinction between corrections of Yukawa-and gauge-type since both are otherwise mixed by the EW symmetry-breaking. A clean separation is only possible in the gaugeless limit (i. e. g EW ! = 0), where the tree-level spectrum slightly differs from the original one: the mass 'm i ' entering the self-energy of 2L order computed in the gaugeless limit and ensuring independence from the 2L field counterterm is not the original tree-level mass, but its gaugeless counterpartm i . The difference with m i corresponds to the order neglected in the gaugeless approximation. In the non-degenerate scenario, the identification of the 'original' Higgs state with its gaugeless counterpart is fortunately straightforward.
The choice of momenta in the 1L 2 terms of Eq. (3) should allow the cancellation of the dependence on the field counterterms of 1L order in the 2L contributions. While we already expanded the momenta in the off-diagonal terms of Eq. (3) and set them to m 2 i , any value p 2 ij = m 2 i + O(1L) is a priori as legitimate as long as dependence on the field counterterms is not considered. Let us therefore analyze the dependence ofΣ (2) ii on the 1L field counterterms. It originates in two pieces: 1L diagrams with counterterm insertion (1L×CT) and tree-level graphs with counterterm squared insertion (CT 2 ): 3 the 1L×CT and CT 2 contributions to the tadpoles cancel one another. Here, Σ ij (orΣ ij ) stands for the unrenormalized self-energy; δM 2 ij (or δM 2 ij ) denote the mass and δZ ji the field counterterms of 1L order. The˘ notation highlights that the corresponding quantities are obtained in the approximations applied in the 2L calculation (e. g. the gaugeless limit).
From the 1L 2 terms emerge the following contributions depending on field renormalization: The requirement for cancellation of the field-renormalization dependence between Eq. (4) and Eq. (5) dictates the following conditions: In other words, the self-energies employed in the 1L 2 terms should be calculated in the same approximation as in the 2L calculation in order to avoid the inclusion of arbitrary UV-logarithms in the calculation. 3 In this analysis of the dependence on 1L field counterterms δZ ij (for commodity we omit the superscript (1) ), we do not display the 2L term δZ (2) ii − k which separately cancels out for p 2 ! = m 2 i , as the vanishing condition for the dependence on the 2L field counterterm δZ (2) ii . This restores agreement between our expansion and the standard conventions for field counterterms-see Eq. (2).

Degenerate case
We assume the existence of a degenerate sector (denoted by 'D'). In this case, one should consider the effective mass matrix M 2 D (s) of Eq. (44), evaluated at a pole of the propagator-refer to Appx. A.2 for a derivation. For two field directions i, j ∈ D, we define the notation m 2 · U I (1) † now, while they could still amount to O(2L) in Ref. [66]; the reason is that we now also want to put 2L effects under control.
Let us now consider the 2L order: Once again, the dependence on field-renormalization constants of 2L order is neutralized by setting p 2 ! =m 2 ij = M 2 I + O(1L) in the argument ofΣ (2) ij . We still need to determine the momenta To this end, we consider the dependence on 1L 2 field counterterms, directly restricting ourselves to the symmetric case δZ ji = δZ ij . Then, instead of the diagonal matrix element of Eq. (4) one has to take the full self-energy matrix for the degenerate sector into account, depending on 1L field counterterms as follows: The dependence of W ij on δZ ij with i, j ∈ D almost cancels (provided all 1L 2 quantities are calculated in the same approximations as the self-energies of 2L order, e. g. in the gaugeless limit), up to a remainder of 3L order: We find no obvious method to absorb this piece by adding a finite term of 3L order to W ij ; thus, after applying the condition D I ! = O(3L), a subleading field-renormalization dependence remains in the determination of the pole of 2L order by This persisting dependence on field-renormalization constants in the degenerate case is intimately related to the form of the dependence of the off-diagonal self-energy of 2L order on 1L field counterterms-see Eq. (12). In the non-degenerate case, this object only intervenes in the mass corrections of 3L order, so that we can expect the field-renormalization dependence to be tackled by the inclusion of the diagonalΣ (3) ii (which, however, goes far beyond our purpose and our current means). The 2L terms of Eq. (15) (contributing at 3L order) have been chosen so as to avoid a quadratic dependence on the field counterterms. Other choices are of course possible. While the full cancellation of the dependence on field counterterms thus fails in the degenerate scenario, we still expect an improvement in the expansion method, as compared to the iterative pole search, due to the careful pairing of 2L and 1L 2 effects.
In case the degeneracy is not lifted at 1L order, we can still define M 2 (2) I as an eigenvalue of the effective mass matrix Indeed, similarly to the 1L piece in Eq. (10), the 2L corrections W ij of Eq. (16) are independent from the chosen pole, the selection being ensured by the projection via the eigenvectors U I . It is then convenient to directly consider these eigenvectors, albeit slightly dependent on the choice of field renormalization, as determining the mixing matrix S in the degenerate subspace at 2L order. In addition, the normalization k∈D S 2 Ik ! = 1 simplifies the application of this object, as can be read from Eq. (18) or the derivation of the Higgs-decay amplitudes in Appx. A.3.
Finally, we stress that it is necessary to compute all the 1L 2 pieces in the approximation of the 2L calculation, i. e. all self-energies, mass counterterms and tree-level masses in W are replaced byΣ (1) , δM andm. An apparent difficulty accompanies the observation that the identification of 'original' and 'gaugeless' states is no longer trivial in the degenerate scenario. Yet, as the full 2L + 1L 2 mass contribution is collected within the block W, it is possible to compute the latter in the gaugeless base and then rotate it to the 'original' base using the gauge eigenbase as reference. This ad-hoc procedure is however a sign that the combination of gaugeless 2L effects with a full 1L calculation is not defined in a completely consistent fashion in the near-degenerate case.
3. Field-dependence in the mass predictions for a non-degenerate scenario We first investigate the numerical significance of radiative corrections to the masses of MSSM Higgs bosons in the clearly defined configuration where all states are non-degenerate, and assess the weight of the dependence on field counterterms in various methods of evaluation. The required Feynman diagrams are computed with the help of FeynArts [74], FormCalc [75,76], TwoCalc [77] and TLDR [31]. The 1L integrals are implemented analytically, while the 2L integrals are numerically evaluated with the assistance of TSIL [78].

Preliminary considerations
We focus on a 'typical' MSSM scenario with squarks of third generation (Q 3 ) and gluinos (g) at the edge of the mass region probed by the LHC (mQ 3 ∼ 1.5 TeV and M 3 ∼ 2 TeV; A t = 2.3 TeV, A b = 1 TeV). EW-only-interacting SUSY particles are given a sub-TeV mass. The ratio of the doublet Higgs v.e.v.-s, t β , is set to 10 and we vary the charged-Higgs mass between 0.5 TeV and 4 TeV. In these conditions, the neutral SM-like Higgs (h), with mass at the EW scale, and the heavy-doublet states (H, A), with masses comparable to that of the charged Higgs (H ± ), are clearly non-degenerate. The CP-even (H) and CP-odd (A) heavy states do not mix in the absence of CP-violation, which we assume in this section. One can thus safely employ the formalism corresponding to the non-degenerate scenario. The lagrangian parameters in the Higgs sector are renormalized in the same scheme as employed in Refs. [66,79], i. e. with cancellation of the tadpoles, on-shell conditions for the EW gauge-boson, SM-fermion and charged-Higgs masses, while t β is defined in the DR scheme.
It is instructive to first consider the corrections of 1L order to the tree-level masses m h i of the neutral Higgs bosons h i ∈ {h, H, A}. In Fig. 1, we show the corresponding shifts in mass squared, The orange curves correspond to the gaugeless limit, while the full 1L shift is shown in blue. The situation of the SM-like Higgs is straightforward: the radiative corrections of Yukawa type O(α t,(b) ), captured in the gaugeless approach, indeed dominate the 1L shift. Contributions of EW type, beyond the gaugeless description, amount to only ∼ 15%. The gaugeless limit can thus be seen as predictive for this state. 5 Nevertheless, the impact of the various contributions follows a different pattern in the case of the heavy-doublet states. A first identifiable feature corresponds to the 'spikes' in the vicinity of M H ± = 3 TeV: these are associated with threshold effects in the loop integrals of the selfenergies, when M H,A ∼ m X + m Y for X and Y two particle species entering the loop diagram (m X and m Y , their masses). At M H ± ∼ 3 TeV, these internal 'on-shell' lines are squarks of the third generation, contributing at O(α q ). The sharp variation in the mass shift is consequently a physical effect, although it is not quantitatively described in the 'free-particle' expansion-we will make no attempt at addressing the threshold behavior by accounting for squark interactions in this paper: see e. g. Ref. [80]. Comparable features also appear close to M H ± ≈ 1 TeV for the blue curve (full 1L corrections) and correspond to electroweakino loops: these do not show in the gaugeless limit as EW interactions are turned off. Beyond this threshold behavior, the pattern of 1L corrections is essentially flat in the gaugeless limit while a slope is definitely identifiable in the presence of gauge effects. As explained in Sect. 3.2 of Ref. [66], the radiative contributions of gauge type to the squared mass-splitting between heavy-doublet states generate a term scaling linearly with the Higgs mass. This leading effect is shown as a dashed green line in Fig. 1 and indeed captures the slope of the full 1L result. The shift between the green and blue curves is due to corrections scaling like M 2 EW ln k M 2 H ± M 2 EW , where k ∈ {0, 1} and M 2 EW ∼ M 2 Z denotes the EW scale. Therefore, radiative corrections of Yukawa type do not capture the bulk of the contribution to the mass-splitting between heavy states and corrections of orders α q α s or α 2 q are not expected to dominate the 2L corrections either. While these orders are fully under control, their inclusion does not improve the numerical precision of the mass predictions to the heavy states as long as dominant corrections of EW gaugetype are not considered. The latter are known as far as scalar self-energies are concerned [31], but are currently not exploitable as 2L contributions to the vector self-energies that are needed for renormalization and the connection to observable input are still missing. Consequently, when we discuss the known O α q α s , α 2 q at the level of the heavy-doublet states below, it is purely for the sake of testing the formalism described in Sect. 2 on controlled orders: the theoretical prediction is actually not improved as compared to the strict 1L result.

Corrections of O(α q α s )
The (S)QCD corrections have a specific status in the contributions of 2L order to the Higgs masses, because no corresponding 1L 2 effects are associated. Thus, this order does not entail the inclusion of any off-diagonal self-energies-unless these are already needed at 1L, i. e. in the near-degenerate case. Given the large impact of QCD corrections to effects of Yukawa-type, they are particularly relevant for an accurate determination of the mass of the SM-like state.
In Fig. 2, we plot the shifts of the squared masses obtained at the considered order, i. e.
with corrections obtained in the gaugeless (gl) limit. The solid orange curves have momentum evaluated at the tree-level (gaugeless) mass of the Higgs states. The cyan dashed lines are derived in the effective-potential approximation, i. e. with momentum set to 0. This approximation is 'exact' in the case of the SM-like Higgs, because the corresponding tree-level mass in the gaugeless limit is indeed equal to 0 (this is not necessarily true in extensions of the MSSM), so that the orange and cyan curves overlap. For the heavy-doublet states, threshold effects originating in squark loops are again prominent (when considering momentum dependence). Even far from the threshold region, the effective-potential approximation does not appear as a particularly useful approach for the heavy states, as its predictions for the order α q α s are O(100%) away from the actual momentum-dependent contributions. This is not surprising for Higgs states with mass differing significantly from 0, even in the gaugeless limit. In this context, the self-energies of order α q α s with vanishing momentum can at best provide an estimate of the corresponding order in an assessment of the uncertainties, but are not predictive. In conformity with the discussion of Sect. 2, the mass shifts of O(α q α s ) presented in orange in Fig. 2 (or those of 1L order in Fig. 1) are independent from the field counterterms, because the self-energies have been expanded for momenta in the vicinity of the tree-level masses, and truncated at the strict desired order-meaning in practice that the self-energies are evaluated at the tree-level Higgs mass. For definiteness, we remind here the formal expression defining the pole mass in this expansion-and-truncation approach for a non-degenerate state h i : In competition with such an expansion, a popular approach consists in iteratively replacing the momenta in the self-energies by the value derived in the pole-mass determination. As off-diagonal effects do not matter at O(α q α s ) in a non-degenerate scenario, 6 we define such a pole search of order α q α s by simply considering the shift of the diagonal element of the effective mass matrix. In other words, the equation that we solve iteratively in this subsection reads denoting the (complex) pole mass-by an abuse of language, we will employ the same notation for its real part in the discussion below. In order to achieve UV-finite results for the renormalized self-energies of Eq. (22) away from p 2 = m 2 h i , the definition of these objects requires the introduction of field counterterms, thereby generating a dependence of higher order on field regulators in M To understand the differences between the expansion and iteration procedures, we investigate the mass shift between both through an expansion, where we only keep leading terms: with all self-energies evaluated at the (gaugeless) tree-level mass, indicating differentiation with respect to the external momentum squared, 'gl' and 'EW' referring to the gaugeless and electroweak contributions respectively. We assumedΣ h i h i (which may be seen as a definition of the EW piece). The systematic appearance of derivative self-energies makes evident the dependence on field counterterms.
Below, we study the stability of Higgs-mass predictions of O(α q α s ) under variations of the field-renormalization constants. For simplicity, we restrict ourselves to a 'minimal' form of the field counterterms, considering on-shell (OS) or DR field-renormalization constants. In the DR approach, the only loose parameter is the renormalization scale µ UV ; its variation between the physical scales of the model, i. e. the EW and the SUSY scales, offers a measurement of the arbitrariness introduced in the definition and UV-regularization of physical observables, hence a lower bound on the associated 'error'. 7 Yet, this very specific pattern where distinct counterterms are correlated by a common regulator may underestimate the actual theoretical uncertainty. In this picture, the MSSM Higgs-field counterterms can be written as where X R,I if encodes the decomposition of the tree-level neutral Higgs state h i on the gaugeeigenbase: Here, we neglect the Yukawa couplings of the first and second generation. The symbol ∆

−1 UV
represents the UV-divergence (including universal finite pieces). The finite piece can be viewed as resetting the renormalization scale of the fields from µ ren to µ UV . Similarly, the 2L fieldrenormalization constants of order α q α s read Only Yukawa and the QCD-gauge couplings appear in these expressions of the DR field counterterms. Since the EW corrections (gauge, gauginos) generate a vanishing UV-divergence at 1L, the corresponding pieces are insensitive to scale variations in the field counterterms. Accordingly, the variations of µ UV that we consider below probe all terms except for the last of Eq. (23). Thus, this type of scale variation is meaningful as long as Yukawa effects dominate and the gaugeless approximation holds. On the other hand, if corrections of EW type are large, the partial higher order appearing in the last term of Eq. (23) is not tested and the uncertainty from scale variation is only partial. The approach with OS-renormalized fields allows to take into account the impact of EW corrections to a certain extent, but it does not allow for variations. In this case, we simply express the field counterterms as cancelling the differentiated self-energies at the corresponding tree-level (OS) mass: which is a symmetric, but (for off-diagonal field counterterms) not fully conventional choice. We could define this object in the gaugeless limit, in which case EW corrections would still be overlooked. In practice, one then recovers results within the scope of the DR scale variation. Therefore we dismiss this choice. More usefully, we can define the 1L field counterterms in the full model, i. e. including EW effects. We can then compare the differences between this procedure and the DR evaluations. However, we do not attempt to express the 2L field counterterms in this scheme, keeping them DR with µ UV = m t , first, because the 2L contributions are explicitly calculated in the gaugeless limit, secondly because the evaluation of the differentiated 2L functions is technically involved.
In Fig. 3, we compare the Higgs-mass predictions of 1L order and α q α s obtained via an expansion-using Eq. (21)-or an iterative pole search-using Eq. (22)-for M H ± = 1 TeV. At the level of the SM-like state (upper row), the mass shift between the 1L and O(α q α s ) predictions is sizable (in our renormalization scheme), and much larger than the dispersion between the expansion and iteration methods. The dependence of the iteration method on the fieldrenormalization constants induces a variation with the regulator µ UV , when evaluating the latter between the EW and SUSY scales, which (in the considered scenario) amounts to about 5 GeV at 1L (long-dashed magenta curve) and about 3 GeV at O(α q α s ) (solid purple curve). This reduction of the scale dependence originates in the destructive interplay between the two orders, not in the completion of the partial order introduced by the pole search-which requires the terms of O α 2 q (see next subsection). Inclusion of the O(α q α s ) corrections in the effectivepotential approximation (short-dashed pink curve) leads to very comparable results, confirming the adequacy of this approach for the SM-like state: the 2L effects are in fact introduced at their correct tree-level gaugeless mass value. The evaluation with OS field counterterms, represented by a cross in the column on the left of each plot, returns results very near that of the expansion.
In the second and third rows of Fig. 3, we consider the shifts in squared mass for the heavy neutral states as compared to the charged one. The plots on the left show the dispersion of the mass predictions depending on the chosen order. Once again, the mass shift associated with the inclusion of O(α q α s ) corrections is larger than that induced by the choice of method (expansion vs. iteration). However, contrarily to the case of the SM-like state, the effectivepotential approximation ostensibly falls far away from the actual momentum-dependent result, implying that this approach is not predictive for the heavy-doublet states. Here, we should comment on the definition of this approach, since it is not completely straightforward. Indeed, beyond the zero-momentum assumption in the neutral Higgs self-energy of O(α q α s ), we also set the momentum equal to 0 in the corresponding charged-Higgs counterterm. 9 However, the charged Higgs self-energy at vanishing momentum requires the introduction of charged-Higgs field counterterms as well in order to obtain UV-finite results. Given that observables are supposed to be independent from field counterterms, the latter could be chosen a priori without connection to the neutral field counterterms: this is the choice governing the short-dashed pink curves, where DR conditions at the scale m t are employed. Alternatively, one can correlate the charged and neutral counterterms in the DR scheme, i. e. employ a common scale µ UV : this is shown in dotted pink, with much reduced variations along µ UV due to SU (2)-cancellations between neutral and charged counterterms-though in no way suggestive of an improved reliability of the effectivepotential method (the results for full momentum-dependent O(α q α s ) are significantly away). 10 We do not consider the OS field counterterms for the charged Higgs in the effective-potential approach, since these are infrared-divergent, which explains that only one pink cross appears in each plot. The relative failure of the effective-potential approximation to capture the O(α q α s ) corrections is not really surprising, as the corresponding Higgs masses are far from p 2 = 0, even in the gaugeless limit. Similarly, the large scale dependence in the dashed pink line proceeds from the necessary inclusion of a large correction proportional to the field counterterm to absorb the UV-divergences of the self-energies at p 2 = 0. [GeV] 200 400 600 800 1000 1200 1400 400 600 800 1000 1200 1400 400 600 800 1000 1200 1400  The plots on the right compare the momentum-dependent mass predictions of O(α q α s ) in the expansion (green) and iteration (purple) methods. We observe that the scale variations are not necessarily sufficient to allow both predictions to overlap-though the size of this dispersion remains compatible with the amplitude of the variations with µ UV . The origin of this separation between the expansion and iteration method is associated with the relevance of EW corrections at 1L for the heavy-doublet states. The last term of Eq. (23) indeed becomes sizable but, as we explained above, the scale variation does not probe the incompleteness of this order. As it is, these partial O α 2 effects are most certainly misleading and should not be interpreted as meaningful contributions of the iterative pole search. In fact, the evaluation with OS field counterterms typically pulls the iterative pole search in the direction of the expansion.
We complete this discussion about the mass prediction at O(α q α s ) for the heavy-doublet states with another example in Fig. 4: M H ± = 2.5 TeV, so that effects beyond the gaugeless approximation, both threshold and EW, are more relevant. The main difference with respect to the previous case is that the choice of procedure (expansion vs. iteration) leads to clearly separated predictions, not connected by the scale variation. As explained above, this shift originates in the last term of Eq. (23), which is not probed by the scale variation, though it is of partial higher order, hence not predictive. Considering the case of OS field counterterms-a setup that is sensitive to this term-extends the range of variation of the iterative pole search to engulf the prediction of the expansion. Once again, the effective-potential approximation provides no actual gain in precision with respect to the strict 1L calculation.

Corrections of
The inclusion of order α 2 q brings about the interplay between 2L and 1L 2 effects that we discussed in Sect. 2. Off-diagonal self-energies indeed become meaningful in the mass calculation even in the non-degenerate case-while they led to non-decoupling SU (2)-violating effects of higher order in a 1L calculation [66]. As explained in this reference, it is crucial, in the case of the heavydoublet states, to properly include 1L 2 off-diagonal contributions to the charged-Higgs mass in the 2L counterterms.  Once again, we first investigate the mass shifts generated at this order, in Fig. 5. Here, with all objects defined in the gaugeless limit, one has As before, the effective-potential approximation is only meaningful for the SM-like state, though it seems to work somewhat better for heavy states at low values of M H ± than at O(α q α s ). Threshold effects again appear in association with squark loops for M H ± ≈ 3 TeV, when the momentum dependence is accounted for. The background far from threshold contributions is essentially flat, as radiative The Higgs masses M 2(2,α 2 q ) h i defined via the expansion in Eq. (28) are by construction independent from field renormalization. Alternatively, one can numerically solve Eq. (1) with the momentum-dependent self-energy matrix including contributions up to O α 2 q , defining the poles M 2(2,α 2 q ) h i . As self-energies are evaluated away from their mass-shell, it is once again necessary to call upon field renormalization to neutralize UV-divergences and give a meaning to corresponding objects. In addition to the 1L field counterterms of Eq. (25), one should consider   the 2L field counterterms of order α 2 q , reading (in the DR scheme) In addition to these 2L field counterterms, the self-energies of order α 2 q depend on the 1L field counterterms of Eq. (25) according to Eq. (12)-even when the momentum is set to the treelevel Higgs mass. As explained in Sect. 2, this dependence on field-renormalization constants is mirrored by that of the 1L 2 terms resulting in a cancellation for the full order. It is also worth noticing that the dependence on the charged-Higgs field already vanishes separately within the 2L + 1L 2 terms forming the charged-Higgs mass counterterm of O α 2 q . In Fig. 6, we show the field-counterterm dependence of the mass prediction at O α 2 q for M H ± = 1 TeV. Once again, we compare the strict expansion-Eq. (28)-in solid green, with the iterative pole search, in solid purple. In this last case, the full CP-even mass matrix of O α 2 q is included in the pole search-with gaugeless approximation for the 2L pieces. For the CP-odd state, instead of considering the full matrix, we directly add the off-diagonal mixing contribution with the Goldstone boson Σ (1,gl) A 0 G 0 p 2 2 m 2 A to the diagonal element-omission of this term would otherwise induce an artificial SU (2)-breaking among the heavy-doublet states, see Ref. [66].
Finally, the dashed and the dotted pink curves include the 2L corrections of O α q α s , α 2 q in the effective-potential approximation before performing the pole search. There, we use two slightly different procedures for the CP-even and CP-odd sectors, as we add a Σ (1,gl) A 0 G 0 (0) 2 m 2 A term in the latter case, instead of keeping a full momentum-dependent 1L self-energy matrix. As before, the charged-Higgs field counterterms are either set to a fixed value (dashed curves) or correlated with the neutral sector (dotted).
At the level of the SM-like state, the corrections of O α 2 q are sufficiently significant to make them necessary in any attempt at precision predictions for the mass. Surprisingly, the variations with the field-renormalization scale are barely reduced as compared to the order α q α ssee Fig. 3-for the iterative approach. This situation originates in the inexact cancellation of the 1L 2 field counterterms that are generated in the variation of the 1L self-energy via the pole search, with those produced in the 2L corrections and assuming the gaugeless limit. A simple estimate actually suffices to recover the order of magnitude of these scale variations. The leading contribution from the pole search is a term − M h because all orders are involved in the mass shift. On the side of 2L self-energies, the dependence on field counterterms is dominated by the corresponding term −δZ hhΣ is large-of the order of 100% of the magnitude of radiative corrections to the mass of the SM-like Higgs-while, dominated by α t , δZ hh ≈ 0.1 for a variation between the EW and SUSY scales. In the aftermath, an effect of ∼ 2% is assessed on the field dependence at the level of the mass. Curiously enough, we can expect a comparable uncertainty from missing EW orders due to the gaugeless approximation. Indeed, as we observed in Fig. 1, the gaugeless approximation works at ∼ 15% at 1L for the SM-like state while the squared-mass shift of O α 2 q is of order 2 · 10 3 GeV 2 ; combining the two numbers, we arrive at a FO uncertainty of percent level from uncontrolled 2L orders at the level of a ∼ 125 GeV Higgs mass-corresponding to the squared logarithms of O(α t α). However, this coincidence does not establish the variation of field counterterms as a realistic uncertainty estimate for higher orders, as it strictly measures partial higher-order effects associated with the regularization of the pole-search procedure. We stress that the large dependence on field counterterms is intimately related to the size of the radiative corrections in the considered approach, i. e. to the FO procedure: in an EFT, contributions from the hierarchical spectrum would be absorbed within the tree-level couplings. Similar behaviours had been observed in Ref. [29]: see e. g. Fig. 3 of this reference. In addition, the mismatch between the total mass shift (Σ (1,αqαs,α 2 q ) hh m 2 h ) and the gaugeless 1L self-energy could certainly be minimized via an adequate choice of renormalization conditions for the Yukawa couplings, absorbing the impact of higher orders.
Concerning the heavy-doublet states, the field dependence induced by the pole search is in fact increased in the CP-even case with respect to O(α q α s )-see Fig. 3. In contrast, the dependence on DR field-renormalization constants remains mild in the CP-odd case, which is a consequence of our adding the off-diagonal self-energies in the gaugeless limit, hence satisfying the cancellation of the field counterterms from off-diagonal terms with the diagonal 2L contributions. Moreover, we observe that even the enlarged field-scale variations of the CP-even mass predictions do not capture the full magnitude of the dispersion between the expansion and iteration methods. Off-diagonal elements in the pole search indeed add further EW terms insensitive to the scale variation beyond that of the last term of Eq. (23); they correspond to partial higher-order α α q and α 2 terms and cannot be interpreted as a genuine physical effect, because the completion of the order α 2 is likely to sizably affect them. Comparing the DR regularization with that of OS field counterterms further extends the magnitude of the variations associated with the iterative approach, making it sensitive to EW effects and filling the gap with the expansion procedure. Finally, for 2L self-energies derived in the effective-potential approximation (dashed pink curves), the scale dependence caused by the approximation p 2 = 0 in the 2L self-energies reaches a magnitude much larger than the mass-squared shift induced by the inclusion of these contributions with respect to the strict 1L order. The predicted masses in the case of correlated field counterterms for the neutral-and charged-Higgs sectors (dotted lines) also fall significantly far away from the actual momentum-dependent 2L prediction-as compared to the shift from predictions of 1L order. This approach is thus meaningless for the heavy states.
To summarize, the inclusion of 2L corrections of orders α q α s and α 2 q is essential in a precise prediction of the mass of the SM-like Higgs state and the effective-potential approximation is reasonably predictive at this level. However we have observed that these orders are of little consequence for the heavy-doublet states in a non-degenerate scenario, since one expects them to be superseded by EW corrections. In addition, the effective-potential approximation works poorly for such states, as the obtained mass shift does not quantitatively improve the predictions as compared to O(1L). We also noted that the dependence on the field counterterm at the level of the SM-like state remains large in an iterative pole search including terms of O α 2 q due to the mismatch between the 1L 2 terms induced by the pole search and the corresponding 2L terms derived in the gaugeless limit: it thus appears that the iterative pole search at FO α 2 q contains an intrinsic uncertainty of ∼ GeV-size (depending on the hierarchy between the EW and SUSY scales). Concerning the heavy-doublet states, the dependence on DR field counterterms appears comparatively reduced for a full calculation of O α 2 q . However, this simply indicates that EW effects, not Yukawa, dominate the corrections to the corresponding masses: the DR renormalization is insensitive to these contributions, but the comparison with an OS regularization (as defined in Eq. (27)) allows to probe them and proves that the apparent gap between predictions from the expansion and the iterative pole search is strictly artificial. We thus conclude that there exists no advantage in employing the more costly iterative pole search in a non-degenerate scenario, as compared to the more straightforward expansion-and-truncation method. On the contrary, the partial higher orders introduced in the pole search may be unphysical and, in any case, they generate an 'uncertainty' intrinsic to the procedure and hardly representative of 'genuine' higher-order effects. Moreover, as the orders α q α s and α 2 q are the leading 2L corrections only for the SM-like Higgs, there is little significance-as long as 2L EW orders are not under control-in maintaining them for the mass determination of heavy-doublet states; in particular the inclusion in the effective-potential approximation induces uncertainties that are larger than the genuine mass shift.

Field-dependence in the mass predictions in near-degenerate scenarios
Higgs mixing in the near-degenerate scenario makes the situation more subtle for the mass determination, and we study its practical implementation in this section. As it is, the complications originate less in the general formalism described in Sect. 2.3 than in the features of the gaugeless approximation and its matching to the full model.

CP-violating mixing between heavy states
In the MSSM, the phenomenologically most relevant scenario with mass degeneracy involves CP-violating mixing between the neutral components of the heavy doublet. With the gaugeless description of O α 2 q , it is actually not possible, even for small mixing, to define loop-corrected masses according to the strict expansion of Eq. (3) since the CP-even and CP-odd components are exactly degenerate at tree level in the gaugeless limit. 11 However, as the mass-splitting between the two neutral components is of EW order, it can be regarded-consistently with the gaugeless counting-as numerically comparable to 1L effects and the near-degenerate formalism of Sect. 2.3 applies.
We turn to a scenario where CP-violation is induced at the loop-level-the tree-level MSSM Higgs sector is always CP-conserving-by the phase of the trilinear soft SUSY-breaking coupling in the stop sector, φ At . In practice, generating a sizable mixing (always scaling with an SU (2)breaking v.e.v.) in this fashion requires a rather large trilinear coupling |A t | as compared to the diagonal soft SUSY-breaking stop masses mQ 3 ,T , which could potentially produce chargeand color-breaking minima-see e. g. Ref. [83] for a recent reference. As our discussion is meant to be strictly illustrative of the mass calculation, we disregard this problem below, and take (dashed purple), for varying φ At . We see that the diagonal mass-splitting at the radiative level vanishes for φ At ≈ π 4 while the off-diagonal mass entry is comparatively large: we thus expect a sizable mass-mixing in this region. The plot in the upper right-hand corner shows the squared mass-splitting between neutral and charged states at order α 2 q in the expansion formalism for degenerate states. The blue and cyan dashed lines correspond to the diagonal entries of the effective mass matrix of Eq. (19), crossing at φ At ≈ π 4 . The other curves correspond to various evaluations of the masses corrected with the CP-violating mixing: eigenvalues of the effective mass matrix in solid lines, diagonal elements after rotation by the mixing matrix defined at order α q α s in dashed lines and diagonal elements after rotation by the mixing matrix defined at 1L order in dotted lines. We see that these various definitions give very close results, even though the maximal mixing is clearly shifted in phase at 1L. In the lower left-hand quadrant, we plot the squared mass-splitting between neutral states at 1L (green), O(α q α s ) (blue) and O α 2 q (purple, orange and red, depending on the definition through eigenvalues or rotation of the effective mass matrix) in the expansion formalism. The diagonal splitting is shown in dashed lines. Once again, we observe the good agreement among definitions at order α 2 q . Finally, in the lower right-hand panel of Fig. 7, we consider the dependence on field counterterms for φ At = π 4 in masses defined with the expansion formalism (red), in an iterative pole search with full momentum-dependence (blue) and in the effective-potential approximation for 2L self-energies, with correlated charged-Higgs DR field counterterms (dotted green) or independent charged-Higgs field counterterms (dashed brown). In fact, the masses obtained with the expansion method are fully independent from the field counterterms, because the remainder of Eq. (17) exactly vanishes due to the exact degeneracy of the two tree-level masses in the gaugeless approximation. For the iterative pole searches, the situation is largely comparable to what we discussed in the non-degenerate case. The predicted masses for the mixed (heavy) states show only a mild dependence on field-scale variations: we should stress here that the CP-violating self-energy is UV-finite without need of regularization by a field counterterm, so that the mixing α q α s , no mix. α q α s , eigenv. α q 2 , α q α s rot. α q 2 , no mix.   The size of the elements of the effective matrix elements is plotted against φ At ; more precisely, the mass-splitting between diagonal elements (dashed curves) is compared to the off-diagonal entry (solid curves) at various orders. Up Right: Predicted mass-splitting between neutral and charged states at O α 2 q , without accounting for the mixing (blue/cyan dashed), or in various evaluations of the 2L mixing. Down Left: Predicted mass-splitting between neutral states at various orders, without accounting for the mixing or including it. Down Right: Field dependence of the predicted masses at O α 2 q for φ At = π 4 , in the expansion formalism (red), with an iterative pole search retaining full momentum dependence (blue), or in the effective-potential description with correlated (green) or uncorrelated (brown) charged field counterterms. The crosses on the left correspond to the predictions with OS counterterms, while the scale variation employs a DR regularization of the self-energies. entry in the mass matrix is itself completely blind to the DR regularization of the fields-adding to the already noted insensitivity to EW effects. On the other hand, the OS field counterterms from Eq. (27) probe both effects: the corresponding mass predictions, indicated by crosses in the column at the left end of the plot, show a more significant deviation respective to those using a scale variation in the DR renormalization, which closes the apparent gap with the masses of the expansion formalism. Once again, this difference appears as a strictly artificial effect originating in the iterative method. Lastly, the effective-potential approximation performs somewhat better than in the non-degenerate scenario, which should be put in perspective with the small value of M H ± . The dependence on field counterterms of this description is made evident by the steeper variation with the scale (when keeping the charged-field counterterms uncorrelated).

CP-conserving mixing with the SM-like Higgs
As explained in the previous subsection, the CP-violating mixing between heavy-doublet states always falls in the degenerate limit, even for small mixing, as long as 2L effects are considered in the gaugeless approximation. Therefore, in order to study the transition between degenerate and nondegenerate regimes, we consider the scenario of degenerate CP-even states, i. e. M H ± = O(M Z ), even though it is now only marginally relevant from a phenomenological perspective due to tight experimental limits-refer e. g. to the related work in Ref. [84]. As the latter rather motivate weak-than strong-mixing scenarios, we do not bother and try to accommodate a Higgs boson with a mass of 125 GeV. This type of mixing also differs from the CP-violating mixing between heavy-doublet states in that the remainder of Eq. (17) no longer uniformally vanishes, meaning that the expansion formalism retains some amount of dependence on field counterterms, which we aim to quantify below.
In fact, the gaugeless approximation for 2L effects forbids an actual 'mass-crossing' in the offdiagonal 1L 2 term processed in the non-degenerate expansion formalism of Eq. (3), because the SM-like state takes a zero-mass in this limit. Therefore, this equation then returns a well-behaved mass-prediction, independent from field counterterms. On the other hand, this approach meets a first issue with the identification of gaugeless and 'actual' Higgs states, since the gaugeless state with mass equal to 0 is always SM-like, while the SU (2)-partner of the Goldstone bosons in the model with non-vanishing gauge couplings is distributed between both tree-level states, according to an angle α − β (α denoting the CP-even tree-level mixing angle with respect to gauge eigenstates): as M H ± narrows M Z , α departs from β − π/2 (for β > π 4 ; while α is fixed to this value in the gaugeless limit), underlining the irrelevance of a naive identification. However, performing a rotation of α − β of the total 2L effects obtained in the gaugeless limit (including 1L 2 contributions for controlled dependence on field counterterms) is ill-defined outside of the degenerate regime. In addition, the recourse to the near-degenerate formalism is originally motivated by the need to consistently process off-diagonal contributions to the Higgs masses intervening at 1L order, due to the near-degeneracy of diagonal terms. On the other hand, it is misleading to maintain this mixing-matrix description when the non-degenerate regime applies, since it generates SU (2) L -violating pieces that are not controlled by the EW-symmetry breaking, see Ref. [66]. It is thus legitimate to worry about defining the transition between both regimes.
In Fig. 8, we consider the same region in parameter space as in Sect. 3, but in the range M H ± ∈ [100, 500] GeV. The first row of plots shows the general perspective of mixing in the CPeven sector. In the plot on the left, we compare the off-diagonal mass-squared entry-we use DR counterterms with µ UV = m t -in the effective mass matrix of Eq. (19) (solid blue curve) to the diagonal splitting (dashed orange curve): in the range M H ± ≈ 120-150 GeV, the mass-squared splitting is smaller than the off-diagonal self-energy, highlighting the need for a near-degenerate formalism. We stress that we include the block of 2L corrections obtained in the gaugeless |M 22 2 (2,α q 2 , eff) -M 11 2 (2,α q 2 , eff) | |M 12 2 (2,α q 2 , eff) | Down: Dependence of the mass predictions at M H ± = 140 GeV on the field regularization for the expansion formalism with degeneracy (green), in an effective-potential description (magenta and purple, depending on whether charged-Higgs field counterterms are correlated with the neutral ones) and in an iterative pole search retaining full momentum dependence (blue). The plot on the left (resp. right) corresponds to the lightest (resp. heaviest) CP-even state. The crosses on the left-hand side of the plots correspond to the OS renormalization of the fields, and the curves to the DR renormalization with varying scale. limit after identification of the gauge eigenbasis, i. e. after rotation by an angle α − β. On the right, we plot the masses of the CP-even states obtained in various versions of the expansion formalism. The dashed purple and magenta lines correspond to the diagonal entries of the effective mass matrix of Eq. (19), crossing at M H ± ≈ 130 GeV. The dotted orange and dark-red curves represent the masses obtained in the non-degenerate formalism, i. e. according to Eq. (3). These noticeably depart from the diagonal entries of the mixing formalism at low M H ± , mainly due to the absence of corrections accounting for the tree-level mixing (i. e. no α − β rotation), which are not straightforward to include in a meaningful way in the non-degenerate scenario. Finally, the eigenvalues of the effective mass matrix are plotted with solid green lines. Obviously, the predictions from the mixing formalism barely differ from those of the non-degenerate one above M H ± ∼ 160 GeV.
In the second row of plots in Fig. 8, we study the difference between the degenerate and nondegenerate formalisms in the higher range of M H ± . The latter is obviously small-a few 100 MeV for M H ± 160 GeV. In addition, the dependence of the effective mass matrix on the field counterterms is made obvious by the dispersion among the various choices (DR with µ UV = M W , DR with µ UV = M SUSY , OS). However, this dispersion remains within O(10 MeV), although it tends to increase as M H ± reaches 500 GeV, then outstretching somewhat the regime of validity of the near-degenerate method. Given the good agreement between the predictions obtained with Eq. (3) and Eq. (19) in the range M H ± ≈ 200-400 GeV (corresponding to a 10%-1% mixing), there is no difficulty to extrapolate between the two. However, the combined results will still receive an uncertainty of order O(100 MeV) from this extrapolation (in addition to other sources of theoretical uncertainties). It is actually unclear whether the non-degenerate regime should not be altogether preferred in this intermediate regime. The difficulty here consists in estimating the uncertainty from neglected off-diagonal EW effects: the corresponding mixing entry in the effective mass matrix indeed generates partial effects of EW 2L order, which are not quantitatively reliable-in particular because they do not necessarily satisfy the symmetries of the system, see Ref. [66].
In the last row of Fig. 8, we compare the dependence on the field-renormalization scale in the expansion, effective-potential and iterative strategies for a point with near-maximal mixing (M H ± = 140 GeV). Expectedly, the field dependence in the expansion approach, of the order of ∼ 10 MeV, is much smaller than that obtained with the other approaches, of GeV order. As announced in Sect. 2, this is related to the careful pairing of 2L and 1L 2 terms in the expansion, which limits the contamination of the mass prediction by partial higher-order effects.
Finally, we explained in Sect. 3 that the 2L corrections of O α q α s , α 2 q are not really quantitatively meaningful for heavy-doublet states, since one then expects larger EW 2L effects. On the other hand, these orders are known to be dominant for the SM-like state. Therefore, with M H ± ∼ M Z , it appears necessary to keep these O α q α s , α 2 q corrections for the full effective mixing matrix. In Fig. 9, we compare the masses M (α 2 q ) obtained with full O α q α s , α 2 q corrections (as we always considered them till now) and those obtained in the approximation where these corrections of order α q α s and α 2 q are only applied in the SM-like direction, denoted as M (α 2 q ,SM) . The solid green lines correspond to the masses derived in the degenerate formalism, the dashed red and orange ones to those derived in the non-degenerate formalism. Obviously, the introduction of O α q α s , α 2 q for the non-SM states and mixings has little impact for all states in the parameter space above M H ± ∼ 300 GeV, and may as well be neglected since these orders are not quantitatively predictive. On the other hand, they contribute significantly in the mixing regime, justifying a complete inclusion at low M H ± .

Three-state mixing
For completeness, we study a scenario involving three-state mixing, although of little phenomenological relevance in the MSSM due to strong experimental constraints on the properties of the observed Higgs state. Such a setup has been considered in particular in Refs. [85,86], in view of studying interference effects close to degenerate scalar resonances, e. g. in the 'toy'scattering bb → h i → τ + τ − . The corresponding cross-section is then dominated by the S-channel exchange of Higgs bosons indeed. However, the formalism employed in these references fully disregards the dependence on field counterterms, both in the definition of the full propagator matrix and subsequent approximations, implying a sizable 'uncertainty' associated to regulators, as we show below. Instead, we prefer to describe such phenomena through the formalism derived in Appx. A, i. e. in a fashion minimizing the dependence on field counterterms (and linear gauge regulators). For simplicity, we focus on the scattering bb → h i → τ + τ − via a trio of near-degenerate Higgs states h i . A 'naive' Feynman-diagrammatic calculation would converge very slowly, due to the difference between the MSSM tree-level masses and the actual poles. It is thus useful to directly resum resonant effects. When performing this operation, one needs to extrapolate the form of the effective propagator away from tree-level Higgs masses. A direct resummation of self-energies as achieved in Ref. [85] then explicitly contaminates the pole values with gauge-dependent partial higher-order effects and field regulators. Instead, the identification of poles, residues and effective couplings via the expansion method, as proposed in Appx. A, minimizes the dependence on field counterterms and accordingly distributes the residues between pole and effective couplings, leading to an a priori more predictive result.
Below, we consider a scenario similar to that of Sect. 5.2 of Ref. [85] with M H ± = 175 GeV, tan β = 50, φ At = π/4. A first formal difficulty in such a setup with M H ± ∼ M EW originates in the mixing of (charged and neutral) Higgs states with Goldstone and gauge bosons: in order to avoid problems of consistency with the gaugeless limit, we assume that this mixing can be processed in the perturbative (non-degenerate) description. The latter is justified in the gaugeless evaluation where e. g. Σ (1,gl) H + G − M 2 H ± /2 M 2 H ± ≈ 5 · 10 −4 , but also in the full model, with Σ (1) We may then focus on the three-state mixing in the neutral sector and consider the effective mass matrix of Eq. (19) or, alternatively, a mass derivation through The effective couplings employed in the expansion approach are obtained at 1L from Eq. (53), which actually corresponds to a decay amplitude for the (loop-corrected) Higgs states. Such objects are independent on the choice of scheme for the fields. There exists no particular inconsistency-only added uncertainty-in working with different orders for the mass determination and the effective couplings, except possibly in the identification of the correspondence between poles and couplings. The latter offers no difficulty in the scenario with weak mixing considered here, and can always be made straightforward by using a mixing matrix defined at 2L order (at the cost of introducing a small dependence on field counterterms). QCD logarithms and t β -enhanced corrections to the Higgs-bottom couplings are factorized out and resummed. Turning to the poles obtained in the iterative procedure, we adjoin to them effective couplings that are derived according to the recipe of Ref. [69] for decay amplitudes, objects that are then also explicitly dependent on the choice of field renormalization. We then combine these objects to define the cross-section σ[bb → H 1,2,3 → τ + τ − ] in the Breit-Wigner description. These quantities are also gauge-dependent, see Ref. [66], but we only consider the 't Hooft-Feynman gauge here.
The results for the bb → τ + τ − scattering mediated by neutral Higgs states are shown in the left panel of Fig. 10. The cross-sections obtained with the expansion method, in green and cyan tones, are hardly distinguishable from one another, illustrating the weak dependence of this description on field renormalization. Conversely, the orange, red and magenta curves, obtained with different prescriptions for the field counterterms, demonstrate the strong dependence of the iterative method with off-shell momenta on these regulators: the predicted cross-sections then come with an in-built uncertainty of order 100%, strictly induced by the field dependence of the would-be (pseudo-)observables. Here, we note that this strong disparity among the cross-sections is mostly driven by the imaginary parts of the poles (Γ h i ): the amplitude of the resonances indeed scales like Γ −2 h i , so that moderate fluctuations of these quantities result in enhanced effects at the level of the scattering. The effective Higgs couplings to SM fermions, g H i f f eff with f ∈ {b, τ }, are shown in the right-hand panel of Fig. 10. With the expansion method (left-most column), the couplings of each Higgs resonance cluster at a definite point in the complex plane, irrespectively of the choice of scheme for the field renormalization. With the iteration method (right-most [fb] Hibb , exp. iter. column), this clustering is still perceptible but looser, again highlighting the dependence on the choice of scheme. We note that the property L cannot be maintained with a complex mixing matrix, as derived at 2L order, which reveals the presence of partial 2L pieces in the couplings thus defined. However, at 1L order, it would also be possible to restore this property by employing a real mixing matrix, as in Ref. [66].
To summarize on these mass predictions in near degenerate scenarios, we have seen how, in many cases, this formalism was rather forced upon us by the gaugeless approximation for 2L effects, than by the actual size of the Higgs mixing. The exact degeneracy among the SU (2) L partners of the charged Higgs at the tree level in the gaugeless limit thus forbids the description of 1L 2 CP-violating mixing terms in the formalism of Eq. (3). Similarly, the tree-level mixing angle α − β of the full MSSM (including gauge terms), relevant for M H ± ≈ M EW , can only be accounted for in an ad-hoc fashion in the gaugeless description, requiring a mixing-matrix formalism. This situation is not surprising as it naturally emerges from the neglected EW effects in the gaugeless description and can only be properly addressed by the inclusion of full EW 2L corrections. As it is, we observed that the expansion formalism for near-degenerate states, represented by the effective mass matrix of Eq. (19), provides predictive mass (pole) observables, showing little or no dependence on the choice of field renormalization, when the converse defect limits the usefulness of the iteration method.

5.
Resummation of UV-logarithms of O α q , α q α s , α 2 q and field dependence In the previous sections, we have derived the Higgs masses and discussed the field dependence in a strict expansion at FO. Nevertheless, the large size of the radiative corrections of O(2L) to the mass of the SM-like state in such an approach points at the slow convergence of the perturbative series in the presence of heavy (multi-TeV) SUSY particles. As is well-known, UV-logarithms of the type ln k M 2 SUSY M 2 EW -with M SUSY typically corresponding to the mass of the squarks of third generation-are responsible for these large effects and should be resummed for numerically meaningful predictions. One usually turns to the EFT framework-see e. g. Refs. [32][33][34][35][36][37][38][39][40][41][42][43][44][45]-in order to implement this resummation. However, for the corrections of O α q , α q α s , α 2 q that we discuss here, this operation can be performed directly in the FO context: all the relevant parameter input-strong gauge-and Yukawa couplings-is indeed accessible at low energy from SM observables, so that a matching at high scale is superfluous. Then, the (simple and squared) UV-logarithms contained in the FO expansion can be explicitly extracted and re-molded according to the flow of the Callan-Symanzik equations applying to the SM-like Higgs mass (or the mass-splitting between heavy states). While the method is thus formally distinct from that of an EFT, it is very similar at the technical level to the hybrid procedure described in Ref. [87], since both subtract UV-logarithms from the FO calculation to re-inject them in resummed form.
The relevant Renormalization-Group Equations (RGEs) are those obtained after screening-off the heavy fields-due to the large mass that cut these off in the emergence of UV-logarithms at the level of the loop functions: in other words, they match the field content of the SM, or of the Two-Higgs Doublet Model (THDM) for intermediate scales (if M H ± M SUSY ). Further thresholds could be considered, depending on the relative scales of the gluino and the squarks of third generation, but we will not discuss them here. For the considered orders, the RGEs are entirely determined by the running of the quartic Higgs couplings, the EW v.e.v., the strong gauge-and Yukawa couplings. For definiteness, we collect the resummation formulae for the SM-like state in the gaugeless limit below: The resummation of gaugeless orders has little significance for the heavy-doublet masses, where EW corrections are dominant, and we thus omit the corresponding formulae. Contrarily to the EFT description where the EFT parameters are run down from a UV-matching scale to the Higgs (EW) scale, the RGEs of the FO description are run up from the SM (or THDM) input scale M EW -set equal to m t in practice-towards the scale corresponding to the heavy screenedoff particles. The conversion of observable input (e. g. fermion and gauge pole masses) to running (MS) parameters generates further next-to-leading logarithms. The RGEs summarized in Eq. (30) incorporate all the logarithmically-enhanced corrections that one can derive in the gaugeless limit at 2L order. We thus explicitly restrict ourselves to the leading order (LO) and next-to-LO (NLO) that are explicitly contained in the diagrammatic calculation of Higgs self-energies at 2L FO. However, RGEs of higher order could also be employed-in that case, without subtracting the higher-order logarithms in expanded form, since these have no equivalent in the diagrammatic calculation of O(2L). Furthermore, a few subtleties associated with the Higgs self-coupling parameter (δλ) (αq) that are generated by loop effects of O(α q ) may be worth discussing. Indeed, the traditional (EFT) counting would associate to this object the boundary condition (δλ) (αq) (M SUSY ) . This is also the choice that allows to reproduce the explicit logarithmic expansion of the FO calculation. Nevertheless, the choice (δλ) (αq) (M EW ) ! = 0 is also a perfectly legitimate condition due to the fundamental ambiguities in order counting-α q ∼ α and α q ln M 2 SUSY M 2 EW ∼ 1-which allow the transmutation of the 2L quantity α 2 q ln 2 M 2 SUSY M 2 EW into the 1L object α ln M 2 SUSY M 2 EW : instead of being resummed as logarithms of Yukawa type, the squared logarithms thus overlooked would be left as simple EW logarithms-and possibly resummed as such when absorbed within the quartic Higgs coupling λ(M EW ). As we do not perform this EW resummation here, however, the numerical difference between the two procedures can be sizable at large squark masses. To keep the comparisons simple, we therefore adopt the traditional counting by default, i. e. set With LO and NLO logarithms of O α q , α q α s , α 2 q properly resummed (up to another caveat that we discuss below), the higher-order uncertainty, controlled by terms ∼ α q α 2 s ln 3 M 2 SUSY M 2 EW ( 10-30%, at the level of the squared mass, for M SUSY = 1.5-10 TeV) in the strict expansion, is now pushed back-as far as gaugeless orders are concerned-to ∼ α q α 2 s ln M 2 SUSY M 2 EW (< 1%). In fact, the higher-order uncertainty is now controlled by EW 2L effects ∼ α q α ln 2 M 2 SUSY M 2 EW (∼ 0.2-2% for M SUSY = 1.5-10 TeV). Furthermore, another term of gaugeless order appears in the RGEs, originating in Higgs self-interactions: As this is a contribution of 3L order, it does not appear in the 2L calculation at FO. In order to resum all NLO logarithms of Yukawa-type, we add this term to our calculation nonetheless, instead of leaving it as a contribution of O(α). Finally, it is fair to mention in this short uncertainty estimate that not only logarithmically-enhanced contributions may require the recourse to resummation techniques: see e. g. Ref. [88] for a resummation of squark-mixing terms. EW logarithms of LO could also be resummed, as is routinely performed in the EFT description. Our method with input at the low-energy boundary continues to apply, without need of a matching at high scale. EW gauge couplings are indeed well-defined at the low-energy end. As to the quartic Higgs couplings, they can be obtained from the observables (Higgs masses and decays) that are predicted in the FO calculation, without or with partial UV-resummation. This is a straightforward recipe for a low-energy effective SM-only the (preliminary) mass prediction for the SM-like state is then needed-but accounting for a THDM threshold becomes more intricate as a larger basis of observables is then needed in order to identify the more numerous parameters of the Higgs potential; yet, the situation is no different at this level in an EFT approach, in principle. In fact, the somewhat less robust boundary condition for the Higgs self-interactions is the only drawback as compared to the EFT procedure: we anticipate it to amount to little at the numerical level, albeit a detailed comparison (beyond our current scope) would be needed in order to ascertain this claim. In this paper, we choose not to carry out the resummation of EW logarithms and go through the intricacies of a THDM threshold, first because this would exceed the bounds of gaugeless orders that we meant to discuss here, then because these effects remain comparatively small (and not necessarily related to the scale of the sole gluino and squarks of third generation), finally because the actual focus of the paper concerns the dependence on field counterterms and that the gaugeless orders are quite sufficient to measure its impact in a resummed mass prediction for the SM-like Higgs.
In the upper row of Fig. 11, we compare the explicit calculation of Higgs self-energies at FO in the gaugeless limit and the corresponding expansion in ln k M 2 SUSY M 2 EW . This latter expansion is obtained from the integration of the RGEs of Eq. (30) in a linear way, but-as mentioned above-additional logarithms of NLO intervene from the conversion of parameters in our scheme to MS ones at the scale M EW -e. g. m OS t = m MS t (m t ) 1 + 16 3 αs 4 π + O(α q ) . The plots in the middle and lower rows show the impact of the resummation of the UV-logarithms for the prediction of the SM-like mass. Here, the UV-resummation is achieved by subtracting the identified UVlogarithms and substituting a resummed version where Eq. (30) is integrated numerically. The corresponding evaluations are shown in red and orange lines, while the blue and cyan lines correspond to the predictions without resummation. These calculations are conducted for both the expansion formalism (blue and red) and for the iterative pole search (cyan and orange) with DR regularization of the Higgs fields (µ UV = M W in solid and µ UV = M SUSY in dashed lines). In the lower row, the orange and cyan curves correspond to the approximation p 2 = 0 in the contributions of 2L order processed in an iterative fashion: their colors match those of the corresponding curves in the middle row, accounting for full momentum dependence.
On the left-hand side of Fig. 11, we consider an MSSM scenario with decoupling squarks and gluinos (mQ , keeping the mixing in the stop sector, X t ≡ A t − µ/t β , minimal, and setting M H ± = 1 TeV, t β = 10 for the THDM sector. This is the ideal setup for 1L, gaugeless full calc.
2L, α q α s , gaugeless log. 2L, α q 2 , gaugeless log.  checking the agreement between the logarithmic expansion and the FO calculations, as can be observed in the very narrow matching of the various self-energies. Turning to the mass predictions, the resummation of UV-logarithms accounts for a shift of already ∼ 1.5 GeV at M SUSY ∼ 1.5 TeV and over 10 GeV at 10 TeV. As already noted in earlier works-see e. g. Ref. [9]-the resummation counteracts the tendency of mass predictions of O α 2 q to fall at large M SUSY -due to large ln 2 M 2 SUSY M 2 EW added linearly in the strict FO expansion. On the right-hand side of Fig. 11, we maintain the stop mixing at |X t | ≈ M SUSY , which generates sizable shifts between the logarithmic expansion and the FO calculation (without endangering the relevance of the resummation of UV-logarithms). At least the leading orders from squark mixing are properly included within the non-logarithmic terms of the FO calculation.
As to the dependence of masses on field counterterms, in case an iterative pole search is employed, it is obvious that the resummation affects it only marginally, as the momenta evaluated in self-energies are still shifted with respect to the tree-level mass (the latter choice ensuring invariance in the expansion method). Correspondingly, the uncertainty associated with field variations steeply grows with increasing M SUSY in the FO approach with iterative pole search, reaching O(10%) at M SUSY ∼ 10 TeV-hence becoming larger than the actual higher-order uncertainty after resummation of the UV-logarithms. The approximation p 2 = 0 in pieces of 2L order of the iterative pole search (lower row) tends to systematically over-estimate the mass as compared to the prediction with full momentum dependence, but does not reduce the magnitude of the dependence on field-renormalization constants. Given the similarity of the procedure employed in the hybrid calculation of Ref. [87], then refined in Refs. [29,37,41,[46][47][48][49], the corresponding predictions of the public code FeynHiggs [46,47,57,60,67,87,89,90] are thus a priori subject to the large error associated with the field dependence. Yet, an ad-hoc choice of field counterterms, derived by comparison of the logarithms with those of an EFT (see Ref. [29]), then restores a more predictive behavior-largely compatible with that of the more straightforward expansion approach. We aim at a more detailed comparison in Fig. 12 for the same scenario as in the left-hand column of plots of Fig. 11 (X t ≈ 0).
FeynHiggs-2.18.0 delivers Higgs-mass predictions in the MSSM. The 2L corrections of order α q α s and α 2 q are included in the effective-potential approximation, which is a relevant choice for the SM-like state, as already discussed. UV-logarithms can optionally be resummed and included in the hybrid approach described in Refs. [46][47][48][49]87]. The default setup (labelled 'vanilla FeynHiggs' below) employs an iterative pole-search algorithm to find the loop-corrected Higgs masses. In order to neutralize the large unphysical scale dependence that is thus introduced by the field-renormalization constants, a special scheme has been devised in Ref. [29]; it includes all finite SUSY contributions in the field counterterms, thereby yielding the same logarithmic dependence as that of an EFT calculation. The predictions of 'vanilla FeynHiggs' are depicted as reference in all plots of Fig. 12 (blue and red short-dashed lines).
For consistency checks, FeynHiggs also provides some hidden flags that can be set through environment variables. Below, we make use of the following settings: • FHFOPOLEEQ=1 switches the determination of the loop-corrected masses to a partial perturbative expansion-only corrections for the diagonal self-energies are included; in the considered scenario with hierarchical Higgs sector, this is expected to agree relatively well with our calculation using the expansion formalism; • with FHFINFIELDREN=0 the finite terms of the field-renormalization constants are set to zero, hence reverting the field counterterms to a simple DR form; • FHTBSCALE=# interprets the input value for tan β at the scale # GeV, which can largely be seen as resetting the renormalization scale µ ren for this parameter.  In the upper-left plot of Fig. 12, we compare the mass predictions of the partial perturbative expansion available in FeynHiggs (solid curves) with the results using our own expansion formalism (disks). The plain FO determination is shown in blue, while a resummation of logarithms is included in the red curves. In the case of FeynHiggs, the setting loglevel=2 is employed, which corresponds to the resummation of next-to-leading logarithms (NLL) for the case where all SUSY masses are large. 12 We observe an almost exact agreement of FeynHiggs and our formalism for this expansion approach; small differences at the level of predictions with resummed UV-logarithms are of the expected percent-order magnitude for unresummed EW logarithms (in our setup) and can thus be attributed to additional classes of logarithms that are included in the resummation of FeynHiggs. In addition, the resummation of EW logarithms may not be completely suited to our scenario where electroweakino and slepton masses are not varied with M SUSY .
In the lower row of Fig. 12, we compare our predictions (disks and circles) to the corresponding ones obtained with FeynHiggs (solid and dashed curves) when using an iterative pole-search algorithm with DR-renormalized fields. In order to resemble the setup of FeynHiggs, our predictions include the 2L corrections in the effective-potential approximation. The plot on the left-hand side shows the FO results, while the curves on the right-hand side include resummed logarithms. In both plots we display predictions at the renormalization scales µ ren = m t (solid curves or disks) and µ ren = M SUSY (dashed curves or circles). The input value for the DR parameter tan β is interpreted at the scale m t before being run to the chosen renormalization scale and fed as input: we thus make sure that we are comparing the same points in parameter space. The running of the additional DR parameters m b and A b has a negligible impact in the chosen scenario. This whole procedure is meant to emulate the direct variation of finite field counterterms, as we considered it before, since an independent variation is not straightforwardly accessible in FeynHiggs. Again, we are able to recover the predictions of FeynHiggs at FO (left plot) to a good approximation, which confirms our observations as to the large inherent uncertainty associated with an iterative approach to the pole determination. A somewhat larger deviation is visible in the resummed prediction at the SUSY scale: we could not completely understand the origin of this discrepancy, but this may not matter much since the dependence on the choice of field renormalization in such a description spoils the logarithmic behavior anyway. In any case, it is obvious that the introduction of the UV-resummation does not neutralize the intrinsic uncertainty originating in the pole search. We note that the results obtained with µ UV = m t are very close to the-in our opinion more reliable-predictions of the expansion procedure. This good performance is not completely mysterious in the considered scenario. Indeed, with the external momentum set to 0 in the 2L self-energies and a sizable hierarchy between M EW and M H ± , the difference between the iteration and expansion procedures for the determination of the mass of the SM-like Higgs can be estimated as (0) (similarly to Eq. (23)). Then, in the gaugeless approximation,Σ SUSY , which vanishes for µ UV = O(m t ) at small X t . The inclusion of EW orders does not spoil this picture, as long as large logarithms do not develop from the electroweakino loops. Therefore, the apparent predictivity of the choice µ UV = m t is very specific to this scenario and not reliable on fundamental grounds.
Finally, we turn to the predictions by 'vanilla FeynHiggs'. These appear to be relatively close to our results obtained with the perturbative expansion (upper plot). Again, this comparative proximity is not really mysterious, as the injected counterterms in FeynHiggs have been designed such that one recovers the correct logarithmic behavior. Such a choice is comparable to our OS scheme for the field counterterms, for which we also observed comparative agreement with the expansion in similar setups: see the pink cross in the upper plot of Fig. 6. Therefore, although FeynHiggs employs an iterative pole search by default, it escapes the large uncertainty associated with this procedure through a judicious choice of field counterterms. Yet, a discrepancy reaching O(1 GeV) at M SUSY = 10 TeV is visible in the strict FO approach-blue curves; both calculations are comparable in that they consider exactly the same orders. Unquestionably, this difference originates in the choice of procedure (expansion vs. iteration), since the simplified expansion available in FeynHiggs agrees with our method. This is therefore the magnitude of the error-again, we stress that this shift has no predictive value-contained in this choice for the mass determination, which should be included as an irreducible uncertainty to the predictions (beyond estimated higher-order effects). This contribution is somewhat reduced after inclusion of the UV-resummation (in this scenario), because M 2 h then comes closer to m 2 h −Σ (1,gl) hh (0) (≈ 162 2 GeV 2 at M SUSY = 10 TeV)-the difference between these two quantities controls the leading dependence on the field counterterm, proportional to δZ hh .
In addition, the use of the technically more involved iteration method remains a choice of questionable efficiency since the predictivity of the FO calculation, directly accessible with the simple expansion procedure, is first wasted, then restored through the cross-reference of the logarithms with the EFT method. Finally, it is unlikely that this method can simultaneously produce predictive results at the level of the mass-splitting among heavy states. We therefore recommend the use of the more robust expansion and truncation procedure that we described in Sect. 2 after the principle of independence of observables from the choice of scheme for field renormalization.

Conclusions
In this paper, we investigated the dependence of MSSM Higgs masses on field counterterms in a FO approach at 2L. This dependence on regulators originates in the arbitrary regularization of Higgs self-energies away from their (tree-level) mass-shell, and is exacerbated when processing 2L and 1L 2 corrections in independent fashions. This situation is further complicated by the fact that-due to missing EW corrections in 2L vector self-energies-only effects of O α q α s , α 2 q are fully exploitable at 2L, rendering an evaluation in the gaugeless approximation necessary at the technical level. Masses derived in the strict expansion formalism evade these difficulties through a careful pairing of field-dependent pieces and the neutralization of field counterterms. On the other hand, the popular mass determination via an iterative pole search, computationally more costly, retains an explicit dependence on field regulators, which generates an irreducible uncertainty inherent to the method. In the case of the SM-like Higgs state, we have seen that this 'error' already amounts to a few GeV for a SUSY sector at the TeV scale. Concerning the heavydoublet states, unlike the SM-like one, the known orders α q α s and α 2 q are of limited relevance since EW corrections are expected to dominate. We still considered these contributions in order to test the impact of the regularization of Higgs self-energies away from their mass-shell. Then, a scale variation with DR counterterms is insufficient to capture the full extent of the field dependence, which is driven by leading EW effects, and we also introduced an OS regularization for comparison. In all cases, variation of the field regulators demonstrated that the mass-shifts generated by the iterative pole search with respect to the expansion approach are purely artificial in nature, and we thus believe it justified to prefer the simpler-and field independent-method. This argument adds to the one of the symmetry considerations that we raised at 1L order in Ref. [66].
In the presence of mass degeneracies, the expansion formalism can be extended to account for mixing effects in a fashion keeping the dependence on field regulators to a minimum. We studied several scenarios involving large mixing effects as well as the transition with the non-degenerate regime. The inclusion of 2L effects in the gaugeless approximation has various consequences at this level, such as imposing the degeneracy of heavy Higgs states in the presence of CP-violating mixing or complicating the connection between the gaugeless and 'full' tree-level Higgs states. At a technical level, it would thus seem desirable to put full 2L EW corrections under control, which would allow to escape the constraints of the gaugeless approximation and limit the use of a near-degenerate formalism strictly to scenarios with large mixing. As far as the dependence on field counterterms is concerned, the situation is very similar, however, in the non-degenerate and near-degenerate scenarios.
In addition, we observed that the effective-potential approximation for 2L self-energies leads to quantitatively reliable results only when applied to the mass of the SM-like state. Considering the numerical cost of evaluating 2L integrals at non-vanishing momentum, as well as the irrelevance (in the absence of EW 2L corrections) of the orders α q α s and α 2 q for heavy-doublet states, it makes limited sense, in the non-degenerate case, to consider 2L gaugeless self-energies for any other external state than the SM-like one. Furthermore, the limited pertinence of this approximation for non-vanishing tree-level masses questions as to the applicability of corresponding calculations to extensions of the MSSM, away from an MSSM-like regime, since tree-level masses do not necessarily vanish then, even in the gaugeless limit.
Given that large UV-logarithms develop with increasingly heavy SUSY spectrum, the strict FO formalism suffers from a slow convergence of the perturbative series. This issue can be evaded through an explicit resummation of logarithmic effects. This resummation can be directly included in the context of the FO calculation, without resorting to a matching scale, simply by exploiting low-energy observables as input, and we explicitly performed this operation for the orders α q , α q α s and α 2 q . As this resummation does not modify the problem of the regularization of Higgs self-energies away from their mass-shell, it does not affect our conclusions concerning the dependence of mass predictions at FO on field counterterms. This situation contrasts with the concurrent computation method through EFTs, where artificial field-dependent terms cannot receive large logarithmic enhancement-due to the very structure of the EFT that embeds UV-logarithms within effective tree-level couplings-but not with e. g. the hybrid approach of Ref. [87]. Nevertheless, we have also shown how the judicious choice of field renormalization devised in Ref. [29] for the public code FeynHiggs largely shields the latter from excessive uncertainties associated with the iterative mass determination.
i. e. theΩ i (s) 'rotation' isolates the direction i up to terms of 3L order, with the diagonal element ('eigenvalue') reading Coming back to the propagator matrix, we then havě where the ellipses represent non-singular pieces in the vicinity of s ∼ m 2 i . We thus obtain a pole