Two-loop top and bottom Yukawa corrections to the Higgs-boson masses in the complex MSSM

Results for the two-loop corrections to the Higgs-boson masses of the MSSM with complex parameters of Oαt2+αtαb+αb2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {O}}{\left( \alpha _t^2+\alpha _t\alpha _b+\alpha _b^2\right) }$$\end{document} from the Yukawa sector in the gauge-less limit are presented. The corresponding self-energies and their renormalization have been obtained in the Feynman-diagrammatic approach. The impact of the new contributions on the Higgs spectrum is investigated. Furthermore, a comparison with an existing result in the limit of the MSSM with real parameters is carried out. The new results will be included in the public code FeynHiggs.


Introduction
After the discovery of a Higgs boson [1,2] with a mass around 125 GeV, intense studies were performed to reveal its nature. Although within the present experimental uncertainties the measured properties of this new boson are consistent with the expectations for the Higgs boson of the Standard Model (SM) [3,4], it could be part of an extended model like the theoretically well motivated minimal supersymmetric Standard Model (MSSM). In the MSSM the observed particle could in principle be interpreted as one of the three neutral physical Higgs bosons. At the tree level, the physical states are given by the neutral C P-even h, H and C P-odd A bosons, together with the charged H ± bosons, and can be parametrized in terms of the A-boson mass m A and the ratio of the two vacuum expectation values, tan β = v 2 v 1 . An admixture of these C P eigenstates is introduced to the Higgs sector via loop contributions involving complex parameters from other supersymmetric (SUSY) sectors [5][6][7][8].
Loop corrections to the masses of the Higgs bosons are sizable and therefore phenomenologically very important. Accordingly, numerous calculations for higher-order corrections to the mass spectrum within the MSSM for the a e-mail: sebastian.passehr@desy.de b e-mail: georg.weiglein@desy.de case where C P conservation has been assumed  as well as for the general case of the MSSM with complex parameters (cMSSM) [5][6][7][8][61][62][63][64][65][66][67][68][69][70][71][72] have already been performed. The largest one-loop contributions originate from the Yukawa sector due to the size of the top-quark Yukawa coupling h t , where α t = h 2 t /(4π) . For large values of tan β contributions of the order α b = h 2 b /(4π) , with the bottom Yukawa coupling h b , can become sizable. At the two-loop level both types of contributions receive further potentially large corrections. The dominant contribution is given by the leading O(α t α s ) terms [26][27][28]69] which are known in the MSSM with complex parameters. Additional corrections involving the strong coupling α s are known in the special case of the C P-conserving MSSM [49][50][51][52][53]. Another important class of two-loop corrections are Yukawa-coupling enhanced contributions of the order O α 2 t + α t α b + α 2 b which are known in the C P-conserving MSSM as well [36,39] (the corrections of O α t α b + α 2 b with on-shell parameters are only available in the approximation tan β → ∞ and m b → 0). A computation of the leading corrections of O α 2 t has been published for the general MSSM [70,71]. In this article also the other pieces of the two-loop Yukawa terms are obtained for the general case of the MSSM with complex parameters.
Following the usual convention, we choose to fix the phase of the mass of the electroweakinos, φ M 2 , to zero; then the phase of μ from the superpotential, φ μ , needs to be close to zero or π in order to be compatible with the experimental constraints. The other relevant parameters are the phase of the gluino mass parameter, φ M 3 , and the trilinear softbreaking parameters of the stops, φ A t , and sbottoms, φ A b . These phases, φ M 3 , φ A t and φ A b , are less constrained; especially the bounds on the phases of the trilinear soft-breaking parameters are weaker for the third generation than for the second and first generation.
The calculation presented here extends the Yukawa-type contributions of O α 2 t from Refs. [70,71], and profits from previously developed tools [80]. For this reason the theoretical framework is just briefly outlined and only new aspects are explained in detail in Sect. 2. The numerical analysis in Sect. 3 is focused on the impact of the new contributions on the Higgs masses, showing substantial shifts of 2 GeV and more in certain regions of parameter space. In the limit of vanishing phases of the parameters, our results agree well with previous results in the MSSM for the case where C P conservation has been assumed [39]; differences are shown for the comparison with an interpolation for non-zero phases which so far has been used in FeynHiggs [28,29,68,81,82]. The new results will become part of the public code FeynHiggs.

Higgs masses at higher orders in the complex MSSM
In this section, we briefly outline the theoretical framework for Higgs-mass predictions at higher orders in the MSSM. We introduce our conventions, explain some details of our chosen renormalization scheme, and comment on the gaugeless limit and the bottom mass resummation.

Notation and conventions at the tree level
The two scalar SU (2) Higgs doublets are expressed in terms of their components in the following way: (2.1) After rotation to mass eigenstates, the Higgs potential reads 3) The matrices M h H AG and M H ± G ± are diagonal at the tree level after minimizing the potential. Explicit expressions for the entries are given in Ref. [68].

Gauge-less limit
The gauge-less limit in our calculation is defined by neglecting all couplings proportional to g 1 or g 2 . As a consequence of this approximation the gauge-boson masses M W and M Z are equal to zero in the new two-loop contributions. Accordingly, the Higgs-boson masses entering the twoloop calculation take on the values (2.4) In this limit, the tree-level mixing angles α ∈ [−π/2, 0) and β ∈ [0, π/2) fulfill the relation (2.5)

Higgs masses at the two-loop order
The Higgs mass matrix elements at the two-loop order receive contributions from self-energies, leading in general to mixing of all neutral states. In this article the full oneloop corrections are used, while the O(α t α s ) and the new O α 2 t + α t α b + α 2 b terms are evaluated in the gauge-less limit and at zero external momentum. Therefore, the loopcorrected propagator h H AG is given by h H AG denotes the matrix of the renormalized diagonal and non-diagonal self-energies for the h, H, A, G fields at loop order k, and M (0) h H AG denotes the diagonal tree-level mass matrix.
Mixing of the Goldstone boson (and of the longitudinal Z boson) with the other Higgs bosons yields negligible effects to the propagators of the physical Higgs bosons [83][84][85]. Therefore, in the following we will only consider the (3 × 3) submatrix of h H AG involving the physical Higgs bosons. Though, Goldstone-Higgs mixing is taken into account in subloop renormalization terms of the type (one-loop) 2 [71].
The neutral Higgs masses are derived from the real parts of the complex poles of the h H A propagator matrix, obtained as the zeros of the determinant of the renormalized two-point function, (2.7)

Counterterms
The renormalized two-loop self-energies can be written aŝ with (2) h H A denoting the unrenormalized self-energies at the two-loop order, and δ (2) M Z h H A comprising all two-loop counterterms resulting from parameter and field renormalization. The notation follows [71], where the required expressions for δ (2) M Z h H A can be found. The Feynman-diagrammatic calculation of the self-energies has been performed with the help of FeynArts [86,87] for the generation of the Feynman diagrams, and TwoCalc [88] for the two-loop tensor reduction and trace evaluation. The one-loop renormalization constants have been obtained with the help of FormCalc [89].

Genuine two-loop renormalization
The two-loop counterterms for the Higgs self-energies given in Ref. [71] also apply to the corrections described in the present article. However, there is an interesting difference for the cancelation of the divergence in the self-energy (2) h H (0). The corresponding counterterm reads where terms with products of two one-loop counterterms have been omitted. In the gauge-less limit δ (2) m 2 h H is the only counterterm which contains δ (2) Here we define t β ≡ tan β, s β ≡ sin β and c β ≡ cos β. The two-loop field renormalization constant for the same matrix element is given by (2.11) The two-loop counterterm for t β in the DR scheme and in the gauge-less limit can be expressed as Combining Eqs. (2.9)-(2.12) yields The DR field-renormalization constant δ (1) Z H 1 is a pure UVdivergent term, calculated as the derivative with respect to the external momentum p 2 of the φ 1 Higgs self-energy. The only contribution in the gauge-less limit is a bottom loop, i.e. in the case of the previously calculated O α 2 t corrections [71], δ (1) Z H 1 was equal to zero due to the approximation m b = 0. The terms originating from two-loop fieldrenormalization and two-loop renormalization of t β canceled each other exactly. Now, for the more general case of a non-zero bottom mass, also δ (1) Z H 1 is non-zero and the cancelation is not complete anymore. The genuine two-loop parts of the fieldrenormalization constants, δ (2) Z H 1 and δ (2) Z H 2 , drop out in the gauge-less limit at zero external momentum in Eq. (2.13) because of a cancelation of the contributions in δ (2) t β and δ (2) Z h H . In principle, δ (2) Z H 1 and δ (2) Z H 2 could still appear as field-renormalization constants for the other Higgsmass counterterms. However, also there they drop out exactly (see Eq. (2.23) in [71]): • for δ (2) m Z h since m 2 h = 0 in the gauge-less limit, • for δ (2) H ± in the gauge-less limit, • for δ (2) m Z h A and δ (2) m Z H A since the Higgs sector is C P conserving at the tree level.

Resummation
Radiative corrections to the relation between the bottomquark mass and the Yukawa coupling of the bottom quark h b are proportional to t β . In order to resum the leading t βenhanced contributions, an effective bottom Yukawa coupling is used as described in Refs. [84,[90][91][92][93][94][95], leading to a UV finite and complex correction factor m b . Using a DR renormalization for m b in the MSSM, the largest contributions of this type are captured through an effective bottomquark mass which is given by 14) The symbol m DR,SM b (m os t ) denotes the bottom mass in the DR renormalization scheme, taking into account SM-type QCD corrections, evaluated at the on-shell top mass.
We use the correction factor m b at the one-loop order which is implemented in FeynHiggs. For illustrating the effects seen in our numerical analysis below, we give here the explicit form of the leading contributions: Further subleading contributions involve terms ∝ α b and ∝ α. With this definition a part of the considered two-loop cor- to the Higgs-boson masses are absorbed into an effective bottom-quark mass. In order to avoid a double counting of contributions from the bottomsbottom sector to the Higgs-boson self-energies, the bottommass is renormalized in the DR scheme as specified in Eq. (2.14).

Subloop renormalization
One-loop counterterms for subloop renormalization enter the self-energies (2) h H A in Eq. (2.8). In contrast to the previously calculated O α 2 t corrections, the approximation of massless bottom quarks is dropped in the present calculation. Accordingly, new counterterms for the bottom-sbottom sector are induced, which are specified in the following.
The squark mass matrices in the q L ,q R bases, q = t, b, in the gauge-less limit are given by The squark mass eigenvalues can be obtained by performing unitary transformations, The independent parameters entering the two-loop calculation via the quark-squark sector are the quark masses m q , the soft SUSY-breaking parameters mQ 3 and mq R ,q =t,b, the complex trilinear couplings A q = |A q |e i φ Aq , q = t, b, the complex μ parameter from the superpotential, and the ratio of the vacuum expectation values t β . All of them have to be renormalized at the one-loop level, Here δ (1) Mq denotes the matrix of counterterms after applying the renormalization transformation to the parameters in Eq. (2.16). The renormalization of the top-stop sector, as well as of μ and t β , is carried out as specified in Ref. [71].
For the renormalization of the bottom-sbottom sector, we refer to Refs. [39,50,96] where renormalization of m b and A b in the DR scheme has been proposed to avoid numerical instabilities. Also for the applied resummation of m b the DR scheme for m b is convenient, as explained above. The renormalization scale is chosen to be the on-shell top mass.
• The bottom-quark self-energy in a Lorentz decomposition is given by with the left-vector part L b , right-vector part R b , scalar part S b , and pseudo-scalar part PS b . The bottom-quark mass renormalization is fixed at the on-shell top-mass scale via (2.20) • With Eqs. (2.17) and (2.18) we define (2.21) The renormalization of the soft-breaking parameter A b follows from Eqs. (2.16) and (2.21) with q = b, yielding , where the self-energy is defined with an incomingb 2 and an outgoingb 1 . The renormalization conditions to fix the auxiliary counterterms δ (1) may be chosen analogously to the stop-sector counterterms in Ref. [71]. Again, the renormalization scale in Eq. (2.22) is the on-shell top mass.
• Invariance under SU (2) yields the following relation between the stop and sbottom sector: , and we apply on-shell conditions to both stop particles and the stop mixing angle (see Ref. [71]). Then we choose to make δ (1) (2.24) The quantity δ (1) is the off-diagonal entry of Eq. (2.21) for q = b. It is already fixed by the renormalization condition of Eq. (2.22) for the independent counterterm δ (1) A b .

Numerical results for the Higgs spectrum
In the following numerical analysis the new contributions of O α 2 t + α t α b + α 2 b are added to the known Higgs-mass corrections in the general case of the MSSM with complex parameters which are implemented in FeynHiggs (version 2.12.0). 1 While the improvement by resummation of leading logarithms as described in Refs. [57,58,60] can be applied also to the case of complex parameters (via an interpolation routine), we have not included contributions of this kind in our numerical results presented below. 2 The large impact of the O α 2 t terms has been investigated in Refs. [70,71] and is not presented here again. Instead the focus is set on the new corrections induced by the finite bottom mass. If not stated otherwise, we choose the following default setting for the parameters entering through the new contributions: The quantities in Eq. (3.1a) are real parameters. The charged Higgs mass m H ± is chosen as an input parameter, and its value is set to ensure the compatibility of scenarios with high t β with the current experimental constraints from searches for heavy MSSM Higgs bosons [97,98]. The parameters in Eq. (3.1b) are in general complex. Their respective phases φ A t , φ A b and φ M 3 are scanned in Sect. 3.2. Thereby the gluino mass parameter M 3 does not occur directly in the new Higgs self-energy contributions, but it appears in the leading term of the bottom-mass resummation. The parameter μ is also complex in general, but its phase is constrained to be very close to zero or π by EDM limits (see above). We remark that the phases φ M 3 , φ A t and φ A b are also constrained by EDM limits, but scenarios with large phases are possible (see e.g. Ref. [99]). We show results for the Higgs mass when varying two phases at the same time.
The absolute value of A t has been fixed to yield a lightest Higgs-boson mass close to 125 GeV which can then be identified with the Higgs signal discovered at ATLAS and CMS. Together with mQ 3 and mt R it determines the mass shift which is induced by the stop contributions. We choose different values for mQ 3 and mt R , mb R to avoid numerical instabilities due to degeneracies. Different setups for mt R and X t = A t − μ * /t β are possible to yield a lightest Higgs mass of 125 GeV as can be seen in Fig. 1. Therein the gray bar indicates the mass range 125.1 ± 0.21(stat) ± 0.11(syst) GeV as measured by ATLAS and CMS [100]. 2 For the scenarios in Fig. 1, for which in contrast to the other results shown below FeynHiggs version 2.13.0 has been used, the incorporation of higher-order leading and next-to-leading logarithmic contributions (FeynHiggs flag loglevel=1) would shift the displayed results for M h by 0.6 GeV (blue), 1.7 GeV (red) and 2.7 GeV (green) for X t = 0; −0.1 GeV (blue), 1.0 GeV (red) and 2.0 GeV (green) for |X t |/mt R = 2. Fig. 1 Dependence of the lightest Higgs-mass M h on the stop parameters mt R and X t = A t − μ * /t β as predicted by FeynHiggs-2.13.0 in the version for real parameters without contributions of resummed logarithms. We set mQ 3 = mt R + 100 GeV. The other parameters except mQ 3 , mt R and A t have been fixed to the values given in Eq. (3.1) with vanishing phases.
are comprised as implemented in the version of FeynHiggs for real parameters m tR 1 TeV 2 TeV 3 TeV The absolute value of A b is close to the upper limit from the approximate bound from the requirement of vacuum stability to avoid charge-and color-breaking minima [101,102] (see Refs. [103][104][105][106][107][108][109] for more detailed discussions of this issue).
In the following analysis we call M h the shift of the lightest Higgs-boson mass by the new Yukawa terms of O α t α b + α 2 b , i.e. excluding the previously analyzed contributions of O α 2 t . In Sect. 3.1, the impact of different parameters on the lightest Higgs boson mass in the C Pconserving case is investigated.
We have also investigated the mass shifts of the heavier neutral Higgs bosons. In general, the shifts are of the same absolute size as for the lightest Higgs but with opposite sign. However, since the tree-level input value m H ± needs to be large for high values of t β (where the O α t α b + α 2 b contributions are relevant) to be in agreement with experimental constraints, the relative mass shift for the heavy Higgs bosons is only ≈ 1 0 / 00 . Moreover, the two heavy Higgs bosons receive nearly identical corrections; in the investigated scenarios the largest difference was ≈ 0.1 GeV. For this reason we do not present numerical results for the mass shifts of the heavier Higgs bosons here. It should, however, be noted that even small mass shifts can have an important impact on the resonance-type behavior that typically occurs between the two heavy neutral Higgs states in C P-violating scenarios; see Refs. [110,111].

Scenarios with real parameters
We start to analyze our results by performing a comparison with the previously implemented two-loop corrections in FeynHiggs. The two-loop corrections of O α t α b + α 2 b were up to now only known for the MSSM with real parameters and m A being an input. We compare our new result with the predictions obtained so far with FeynHiggs from both the versions for real parameters 3 and for complex parameters (for the latter employing the same renormalization scheme with m A as input and in the limit of real parameters, but without terms of O α t α b + α 2 b ). Both versions contain shifts due to m b effects (see Eq. (2.15)), including contributions of O(α t α b ).
In Fig. 2  In our following analysis we choose m H ± as an input parameter. In this case the O α t α b + α 2 b terms are new contributions. We investigate the dependence of the prediction for M h on t β , μ and M 3 , whereby all parameters are still kept real. The results are depicted in Figs. 3, 4, 5 and 6.
As can be seen in Fig. 3 large contributions above 1 GeV are only visible at high values of t β . In this scenario M 3 is positive, leading to a much bigger M h if μ is negative, which can be understood from Eqs. (2.14) and (2.15). For later analysis we fix t β = 50.  [112][113][114]. Close to |μ| = mt 1,2 − m t ≈ 1.8 TeV one can see kinks which are induced by threshold effects from the higgsinotop-stop system.
In Fig. 5 the impact of the gluino mass parameter is depicted. This effect enters the Higgs self-energies at the investigated order purely via the employed effective bottom mass. We see a rising shift at growing |M 3 | for opposite signs of μ and M 3 (yielding the same enhancement in m b , see Eqs. (2.15)). At |M 3 | = mb 1,2 − m b ≈ 2 TeV (nearly invisible) threshold effects from the gluino-bottom-sbottom system appear. For our following analysis we fix |M 3 | above that region at 2.5 TeV.
Finally, in Fig. 6 the absolute value of A b is varied, and the resulting mass shift is plotted for positive sign (blue) and negative sign (red) of A b . The difference between both curves, i.e. the impact of the phase φ A b , is enhanced for larger absolute values. However, as too large values of |A b | lead to instable vacua according to the upper limit of Eq. (3.2), we set it to |A b | = 2.5 mb R in the scenarios of the following section.

Scenarios with complex parameters
Various phases enter the self-energies of the Higgs bosons Their impact on the Higgs sector is shown in Figs. 7, 8 and 9. Here we keep μ negative, i.e. sgn[μ] = −1, and M 3 positive, but we could also have chosen the opposite signs of both parameters to see enhanced effects for the phase dependent terms as has been shown in Fig. 4.
We start with the phases φ A t and φ A b . The results are depicted in Fig. 7, where mass shifts between 0.3 and 1.4 GeV can be seen. For φ A b = 0 the variation with respect to φ A t is maximal; the larger the phase of A b , the flatter the dependence on the phase of A t . Similarly, variation of φ A b yields  the largest effects for φ A t = 0. Also the signs of the phases matter, e.g. the mass shifts are different for φ A b = ± π 2 . In addition to the exact calculation (solid lines), FeynHiggs offers an implemented interpolation of the selfenergy corrections that have been known up to now for the case of real parameters but not for the complex case. Since the O α t α b + α 2 b terms were only available in FeynHiggs for the MSSM with real parameters in the limits t β → ∞ and m b → 0, deviations from the new mass shifts can be expected even for real parameters. Besides these relatively small differences, the linear interpolation can differ by ≈ 0.5 GeV from the full result in the investigated scenario. Also the asymmetric behavior for the change of two phases at the same time was not described correctly by the interpolation. Figures 8 and 9 show the influence of varying the gluino phase φ M 3 and in addition either φ A b or φ A t . These terms are induced by the correction factor m b as the investigated class of two-loop corrections does not contain the parameter M 3 . Also here the largest phase dependence is found when one phase is equal to zero. In Fig. 8 the mass shift is nearly symmetric in ±φ M 3 and ±φ A b , i.e. the red and yellow curves are lying on top of each other. Nevertheless, there are small asymmetries in the renormalized two-loop self-energiesˆ h A andˆ H A . On the contrary the mass shift M h in Fig. 9 shows a clear asymmetry similar to Fig. 7.
In summary, phase dependent contributions of O α t α b + α 2 b lead to mass shifts of the lightest Higgs boson of ≈ 1 GeV in the investigated scenarios. The sign of μ has been chosen to be negative in the considered scenarios, but similar effects can be found at positive large μ (and opposite sign of M 3 ).

Conclusions
The two-loop corrections of O α 2 t + α t α b + α 2 b to the Higgs-boson masses in the MSSM with complex parameters have been computed in the gauge-less limit at vanishing external momentum. The terms of O α t α b + α 2 b have only been known in the special case of the MSSM with real parameters before, and were incorporated in FeynHiggs in the limits t β → ∞ and m b → 0. The specific aspects related to the renormalization of these new contributions have been discussed, and their numerical impact on the Higgs spectrum has been investigated.
For the lightest Higgs boson mass at ≈ 125 GeV we have found shifts above 1 GeV at t β > 40 for different scenarios: moderate |μ| = 1 TeV with negative sign and positive M 3 , A t , A b , or with positive sign and negative M 3 , A t , A b . The reason for that enhancement can be found in the large correction factor m b yielding an enhancement of the bottom Yukawa coupling. The effect of varying the phases φ M 3 , φ A t and φ A b can be as large as 1 GeV. If one phase is set close to π , the dependence on the other phases is typically weakened; the largest effects are found when only one phase is varied with all others being zero. In FeynHiggs so far an interpolation of the corrections of O α t α b + α 2 b obtained for the case of real parameters is used for the case of complex parameters. We have found deviations with a size of ≈ 0.5 GeV from this approximation, especially when several phases are different from zero at the same time.
Mass shifts for the heavier neutral Higgs bosons have not been depicted. They are similar to the ones of the lightest Higgs boson; however, with opposite sign. Since we used a large value for t β in our scenarios, we need to choose a rather large input mass m H ± = 1.5 TeV to be consistent with existing experimental bounds. Therefore, the relative size of the mass shifts is small. Moreover, the two heavy Higgs bosons receive similar corrections with a maximal difference of ≈ 0.1 GeV in the investigated scenarios. Nevertheless, small mass shifts can be important to correctly describe the resonance-type behavior of nearly mass-degenerate mixed states like the two heavy Higgs bosons in the MSSM with complex parameters.
The new results will be implemented in the public code FeynHiggs.