Two-loop improved predictions for $M_W$ and $\sin^2\theta_{eff}$ in Two-Higgs-Doublet Models

We present the currently most precise predictions for the $W$-boson mass and the leptonic effective mixing angle $\sin^2\theta_{\text{eff}}$ in the aligned Two-Higgs-Doublet Model. The evaluation includes the full one-loop result, all known higher-order corrections of the Standard Model, and the non-standard two-loop contributions that increase with mass splittings between charged and neutral THDM Higgs bosons. They depend on $\tan\beta$ and the soft $Z_2$-symmetry breakingparameter $m_{12}^2$ of the scalar potential, in addition to the non-standard boson masses. Whenever the one-loop corrections become large, the two-loop contributions yield substantial modifications of the predictions, which is of particular importance for the $W$ mass where large mass shifts are required to reach the recently published final resultof the CDF collaboration. Numerical results are shown for the dependence on the various non-standard parameters and in comparison with experimental data.


Introduction
Precise measurements of electroweak observables at lepton and hadron colliders provide important tests of the electroweak Standard Model (SM) as well as of possible physics beyond the Standard Model (BSM). Model parameters entering the theoretical predictions via loop corrections can be constrained by comparison with the measured data provided the theoretical precision can compete with the experimental accuracy. Prominent cases of precision observables, measured with highest accuracy and very sensitive to virtual heavy particles, are the effective leptonic mixing angle in terms of sin 2 θ eff , and the correlation of the gauge boson masses M W , M Z via the Fermi constant G F . The latter allows a precise prediction of M W from M Z , the electromagnetic fine-structure constant and the Fermi constant, in combination with additional parameters in the higher-order contributions. With the discovery of a SM-like Higgs boson by the LHC experiments ATLAS [1] and CMS [2] all the required input parameters of the SM are determined, and together with the adequate calculation of the loop corrections, acccurate predictions for the electroweak observables have become available.
Recently the final result of the W mass measurement by the CDF Collaboration [7] has been presented, M W,CDF = 80.4335 ± 0.0094 GeV (2) which is in obvious tension with the average (1) from previous measurements and with the SM prediction.
The effective leptonic mixing angle at the Z resonance has been measured in electronpositron reactions e + e − → ff at the colliders LEP and SLC via the forward-backward and τ -polarization asymmetries as well as the left-right asymmetry accessible by longitudinally polarized initial-state electrons that were available at the SLC. The average is given by [8] sin 2 θ eff,exp = 0.23153 ± 0.00016.
The mere extension of the SM by a second Higgs doublet, the Two-Higgs-Doublet Model (THDM), can provide a new type of large contributions to ∆ρ arising fom the Higgs sector, in contrast to the SM where the Higgs sector respects the custodial symmetry and hence deviations of ρ from unity can only originate from other sources, like gauge and Yukawa couplings. The THDM Higgs sector accommodates three neutral and a pair of charged Higgs particles, with masses as independent free parameters; in general, custodial symmetry is broken and thus ρ different from unity can occur for mass splittings between different isospin states of the scalar sector. The non-standard one-loop corrections to precision observables are very sensitive to mass differences between charged and neutral Higgs bosons [68][69][70][71][72][73][74][75][76] and in particular have a large impact on the prediction of M W . Because of the potentially large quantum effects, the THDM has become popular to provide an explanation for the larger value of the W mass given by the new CDF result; accordingly, a series of recent papers deal with analyses of electroweak precision data in various versions of the THDM [77][78][79][80][81][82][83][84][85][86][87][88][89].
The one-loop corrections in the THDM for scenarios with large mass splittings, as required for M W , are quite sizeable; hence calculations beyond one loop are needed in order to get reliable theoretical predictions for precision observables. Since the dominating quantum effects are associated with the ρ-parameter, the leading two-loop effects can be embedded via ∆ρ in case of broken custodial symmetry. In [90], we presented the two-loop corrections to the ρparameter in the CP -conserving THDM originating from the top-Yukawa and the scalar self interactions. This was done under the assumption that one of the neutral CP -even states of the THDM can be identified with the scalar boson observed at the LHC with properties like the Higgs particle of the SM.
In this paper we provide the currently most precise predictions for M W and the leptonic sin 2 θ eff in the CP -conserving THDM combining the two-loop corrections of [90] with the THDM one-loop corrections and the complete set of available SM loop contributions. The outline is as follows: Section 2 specifies the CP -conserving THDM and gives a list of theoretical constraints on the model parameters. A phenomenologically interesting special case is the Inert Higgs Doublet Model (IHDM), which is briefly desribed as well. Section 3 reviews the status of the W -boson mass prediction in the SM and explains the incorporation of the available nonstandard corrections for the W -mass prediction in the THDM. Analogously, Section 4 deals with the effective mixing angle. Numerical results are presented in Section 5, with conclusions given in Section 6.

The Two-Higgs-Doublet model
The Higgs sector of the THDM contains two SU (2) L doublets of complex scalar fields with hypercharge Y = 1, Under the assumption of a discrete Z 2 symmetry ( which is only softly broken by a non-diagonal mass term, the general gauge invariant potential is given by (see e.g. [91]) yields the mass terms quadratic in the fields, which are diagonalized by unitary transformations with in order to disentangle the Goldstone bosons G 0 and G ± from the physical mass eigenstates: two CP -even states h 0 and H 0 , one CP -odd state A 0 , and a pair H ± of charged Higgs bosons. The gauge-boson masses are determined by v and the electroweak gauge couplings g 1 , g 2 , and provide the doublet-specific relation to the electroweak mixing angle, Whereas the combination v 2 = v 2 1 + v 2 2 is fixed by the vector boson masses, there are seven other free parameters of the Higgs potential, which can be expressed in terms of the Higgs boson masses, the mixing angles α and β, and one remaining independent parameter (see [91] for the general relations). Conventionally the soft Z 2 -breaking mass parameter m 2 12 is chosen. In our previous calculation of the two-loop contributions to the ρ-parameter [90], we used the dimensionless quantity from the parameterisation given in [92] instead, which is employed in the FeynArts model-file of the THDM. For consistency we present the results of the W boson mass and the effective leptonic mixing angle also using the combination λ 5 in (15). Furthermore, two additional assumptions have been made in [90], which are adopted also for the present study. First the CP -even state h 0 is identified with the scalar particle observed by the LHC experiments. Second, the alignment limit [93] is applied where the couplings of h 0 to the gauge bosons and fermions are identical to the corresponding couplings of the SM Higgs boson. Formally this can be achieved by setting Flavour-changing neutral currents (FCNC) through neutral Higgs exchange at the tree level can be avoided by the arrangement that not more than one of the doublets couples to fermions of a given charge [94,95], This has lead to four different model classes in the literature which go by the names type-I, type-II, type-X and type-Y (see for example the review [96] for more details). In all of them, however, the couplings between up-type quarks and Higgs bosons are the same. In our calculation of the non-standard two-loop corrections in the THDM all the Yukawa couplings except those of the top quark are neglected, hence the results of [90] are valid independent of the model classification. In the alignment limit, the coupling between the SM-like scalar h 0 and the top quark is identical to the top-Yukawa coupling in the SM, while the couplings between the top quark and the non-standard Higgs bosons A 0 , H 0 and H ± are modified by an additional factor of t −1 β .

The Inert-Higgs-Doublet model
A particular interesting version of the THDM is the Inert-Higgs-Doublet-Model (IHDM) [97] which is distinguished by an unbroken Z 2 symmetry. Under this symmetry all the SM fields and the doublet Φ 1 are even, whereas Φ 2 is odd. The IHDM has received attention in the context of radiative neutrino masses [98] or as a solution to the naturalness problem [99]. Moreover, since the Z 2 symmetry is unbroken, the lightest particle that is oddly charged can provide a dark matter candidate [100]. An overview with more details on phenomenology and additional The couplings of the scalar h 0 to fermions and gauge bosons are identical to the corresponding couplings of the Higgs boson in the SM. Owing to the Z 2 symmetry, the additional particles H 0 , A 0 and H ± do not couple to fermions. Moreover, they can appear only pairwise in interaction vertices and the lightest of these non-standard particles will be stable. If this is one of the neutral bosons, the IHDM provides a suitable dark matter candidate. In order to specify the free parameters of the IHDM, it is convenient to choose, besides m h 0 = M HSM , the masses m H 0 , m A 0 , m H ± , the quartic coupling Λ 2 of the non-standard scalars, and the quantity which is of special interest since Λ 345 describes the coupling of the standard-like Higgs particle to a potential dark matter candidate.

Theoretical constraints
The parameters of the potential given in (6) are subject to various restrictions. A stable vacuum requires the potential to be bounded from below. In the THDM this requirement has to be fulfilled for all possible directions along which the component fields of Φ 1,2 go to large values. As explained in [97,102,103] the conditions ensure that the quartic terms in the potential are positive for large values of the field components in all directions. These tree-level bounds can be improved by considering higher-order corrections to the potential. For more details see the discussion in [96] and references therein. For our analysis we employ the tree-level bounds as an estimate for the allowed parameter range. The unitarity requirement for the scattering matrix puts additional constraints on the parameters of the Higgs potential. Due to the optical theorem, the s-wave scattering length a 0 is restricted to |a 0 | ≤ 1/2. For scattering processes with four scalars in the high-energy limit, a 0 is directly proportional to the scalar couplings. Moreover, due to the Goldstone boson equivalence theorem the scattering of longitudinal gauge bosons can be calculated as scalar-scalar scattering by replacing the gauge boson with the corresponding Goldstone bosons. In the SM, the constraints from tree-level unitarity gives an upper bound on the Higgs mass [105,106]. The application of the analysis in the THDM [107][108][109][110][111][112][113] is more complicated due to the larger number of possible scattering processes and the involved structure of the scalar quartic couplings. With the help of a unitary tranformation, the scattering matrix of the coupled scalar-scalar channels can be simplified by using the original fields φ + i , η i and χ i from (10) instead of the mass eigenstates. The restrictions on the s-wave scattering length constrain the eigenvalues of the scattering matrix at tree-level (see for example [96]) to fulfill |e i | ≤ 8π. We are employing these tree-level bounds as an estimate of the validity of perturbativity. For more accurate restrictions higher-order corrections have to be considered for the scattering processes. A one-loop analysis of the unitarity bounds can be found in [114,115]. The constraints from vacuum stability and S-matrix unitarity are identical in the aligned THDM and the IHDM. In the IHDM the parameters are further constrained by the condition to ensure that v in (17) corresponds to the global minimum of the potential [116].

The M W -M Z interdependence
The correlation between the gauge boson masses M W and M Z can be established via the Fermi constant G F , which is determined with high accuracy from precise measurements of the muon lifetime and the calculation of the muon decay width within the low-energy effective Fermi model including QED corrections up to O(α 2 ) for the point-like interactions [117][118][119][120][121][122]. Comparison of the muon-decay amplitude calculated in electroweak theories like the SM or THDM with the Fermi model result yields the relation where the non-QED loop corrections are summarized in the quantity ∆r. Since it depends on all the virtual particles in the loop contributions, is a model-dependent quantity, and the relation (24) provides the prediction of M W in specific models in terms of the model parameters and the highly accurate input quantitites M Z , G F and the electromagnetic fine-structure constant α.

One-loop calculations
At the one-loop level, the contributions to ∆r consist of the W self-energy, vertex and box diagrams, and the related counterterms. In the on-shell renormalization scheme, the required counterterms arise from charge and mass renormalization (see for example [123,124]), according to with the transverse part Σ V of the gauge-boson self-energies, the photon vacuum polarization and the non-diagonal photon-Z mixing self-energy. In the on-shell scheme the electroweak mixing angle is a derived quantity, yielding and the one-loop counterterm With the notation δ vertex+box for the vertex and box diagram contributions including external wave-function renormalization, one can write At the one-loop level, ∆r can be split into three parts, containing • the shift of the fine-structure constant ∆α from charge renormalization, originating from the light-fermion contribution to the photon vacuum polarization; • the loop correction ∆ρ (1) to the ρ-parameter, which and can be written as • the remainder part ∆r (1) rem summarizing all the other terms.
The loop correction to the ρ-parameter is sensitive to the mass splitting between the partners within an isospin doublet [125]. In the SM, this yields a sizeable contribution from the topbottom quark doublet [125][126][127]. In the THDM, it can moreover get large contributions from the non-standard Higgs bosons in case of mass splittings between neutral and charged scalars, yielding the dominant non-standard loop corrections in ∆r.
In the alignment limit of the THDM, one can identify the loop contributions of the scalars h 0 , G 0 , G ± with the standard scalar contributions to the SM part of ∆r. Consequently, the non-standard contributions to ∆r arise from the scalars H 0 , A 0 , and H ± . Thus, one can write at the one-loop level Since the scalar contributions to the vertex and box corrections are negligible due to the small Yukawa couplings to electron and muon, the non-standard scalars contribute to ∆r only through the gauge boson self-energies, yielding the non-standard part ∆r (1) The subindex NS indicates that only the non-standard parts of the self-energies and counterterms, respectively, have to be taken. The corresponding expressions are listed in the Appendix. The dominant effect on ∆r, as noted above, can be traced back to the additional non-standard correction to the ρ-parameter, entering (33) with The origin of ∆ρ (1) NS are the couplings of the non-standard Higgs sector when they violate the custodial symmetry. ∆ρ

Two-loop corrections in the THDM
The THDM prediction for ∆r beyond the one-loop level can be decomposed according to where ∆r SM contains all the known SM corrections mentioned above and ∆r NS comprises the additional non-standard contribution, which in the alignment case originates from the nonstandard bosons H 0 , A 0 and H ± . The expansion up to two-loop order, accommodates the complete non-standard one-loop part ∆r (1) NS from (36) and the two-loop part ∆r (2) NS which we approximate by including the potentially large terms associated with the ρparameter. These contain products of ∆ρ NS , ∆α as well as the non-standard two-loop irreducible corrections to the ρ-parameter calculated in [90]. Technically, they are obtained in the gauge-less limit (vanishing gauge couplings g 1 , g 2 while keeping the ratio M W /M Z constant) and in the top-Yukawa approximation where only the top-quark mass is kept different from zero. This yields a significant step of improvement at the two-loop level taking into account those contributions that can become sizeable whenever the one-loop contribution in ∆r is large.
In analogy to the SM, the non-standard two-loop contribution can be written in terms of a one-particle irreducible part ∆r (2) NS,irr and a reducible part ∆r (2) NS,red , ∆r (2) NS = ∆r (2) NS,red + ∆r The reducible non-standard contribution, ∆r NS,red , arises from the non-standard one-loop correction to the ρ-parameter (37) in the expansion of ∆r up to two-loop order, in the on-shell scheme given by with The two-loop terms of ∆r which contain only ∆α and ∆ρ (1) t are already included in the SM part; therefore, the non-standard part remains as follows, As a side note, we mention that a reparametrization of ∆ρ (1) → ∆ρ (1) replacing the basic on-shell parameters in the normalization of ∆ρ (1) by the Fermi constant with the help of induces additional factorized two-loop terms in ∆r that allow to rewrite (42) as follows, which corresponds to the expansion of the resummed form found in [14] for the SM.
The irreducible non-standard two-loop contribution ∆r (2) NS,irr originates from the W selfenergy and the counterterms for the parameters in (24). In our approximation, the two-loop charge renormalization constant is zero and the two-loop self-energy Σ (2) W (0) is canceled by the two-loop W mass counterterm. Hence, the only remaining quantity is the two-loop counterterm for s 2 W , obtained from expanding the bare relation (29), which in the gauge-less limit reduces to where ∆ρ (2) is the two-loop correction to the ρ parameter as derived in [90]. The W, Z self-energies Σ (2) W,Z correspond to the set of genuine two-loop diagrams and diagrams with subloop renormalization, i.e. one-loop diagrams with insertion of mass countertems and of δ (1) s 2 W . The two-loop terms with δ (1) s 2 W factorize and can be extracted; in combination with the factorized first term in (49) one obtains The residual part δρ (2) contains the genuine two-loop self-energies combined with the part of the subloop renormalization from the mass counterterms of the internal particles. The quadratic term in ∆ρ (1) from (50) is already included in ∆r 2 NS,red (in the G F -parametrization of ∆ρ it would vanish). The irreducible two-loop part of ∆r is thus given by ∆r (2) it still contains the SM contribution. In the notation of [90], δρ (2) is divided into four finite parts δρ (2) = δρ (2) t,SM + δρ (2) t,NS + δρ H,NS + δρ where δρ (2) t,SM contains the SM-like scalars h 0 , G 0 , G ± and the top-Yukawa coupling, whereas the non-standard scalars H 0 , A 0 , H ± appear in the residual three parts.
The various entries of ∆r have the following properties.
• ∆r (2) t,NS incorporates the coupling of the top quark to the non-standard scalars, which is proportional to t −1 β and enters all the diagrams for this contribution quadratically.
• ∆r (2) H,NS incorporates the non-standard scalar interactions. It is proportional to the squared coupling between three non-standard scalars, which has the form in the alignment limit. Consequently, ∆r H,NS vanishes for Moreover, it is zero for m A 0 = m H ± , since this mass configuration restores the custodial symmetry of the scalar potential. However, differently from ∆ρ (1) NS , this contribution does not vanish for • ∆r (2) H,Mix incorporates the interaction between the standard and the non-standard scalars. Similarly to ∆ρ (1) NS , it vanishes for m H 0 = m H ± or m A 0 = m H ± due to the restoration of the custodial symmetry. Differently from ∆ρ (1) NS , this contribution contains additional couplings between h 0 and the non-standard scalars of the form which can be enhanced by λ 5 . In contrast to ∆r (2) t,NS and ∆r (2) H,NS , the contribution ∆r

Loop corrections in the IHDM
The non-standard one-loop corrections to the gauge-boson self-energies are the same in the IHDM and the general THDM in the alignment limit, yielding ∆r (1) NS for both cases which depends only on the masses of the H 0 , A 0 , H ± bosons. Differences arise at the two-loop level. Since Yukawa interactions of the non-standard scalar doublet are absent, there is no ∆r (2) t,NS term in the IHDM. Moreover, as discussed in [90], the only non-standard two-loop contribution to the ρ parameter, δρ (2) IHDM , follows from the interaction of h 0 , G 0 , G ± with H 0 , A 0 , H ± and is thus equivalent to δρ (2) H,Mix ; it vanishes for equal charged and non-standard neutral Higgs-boson masses and contains the coupling between h 0 and the non-standard scalars, which is determined by the combination Λ 345 of coefficients in the scalar potential, as specified in (20) of Section 2.1. The corresponding contribution to ∆r is given by which introduces an additional dependence on the IHDM parameter Λ 345 , which is not present in ∆r (1) NS and ∆r (2) NS,red .

Incorporation of the higher-order corrections
The SM part, ∆r SM , contains, besides the one-loop result, the complete two-loop and the partial higher-order contributions listed in (38). Since the two-loop electroweak part is quite an inconvenient expression which involves furthermore numerical integrations, we make use of the parametrization given in [43] in terms of the SM input parameters.
Summarizing the various standard and non-standard contributions to ∆r we write (58) as follows, ∆r = ∆r SM + ∆r (1) NS + ∆r (2) NS,red + ∆r where the three non-standard terms are specified in the equations (36), (44), and (53). In case of the IHDM, ∆r (2) NS,irr has to replaced by ∆r (2) IHDM,irr in (57). The predicted value of M W fulfills the relation (24) for a given specific set of parameters. Since ∆r itself does depend on M W , the solution of (24) for M W has to be determined numerically.

The effective electroweak mixing angle
The electroweak mixing angle in the effective leptonic vertex of the Z boson [35] is another important precision observable, measured with high accuracy at the Z resonance by the LEP and SLC experiments [8]. Theoretically it is derived from the ratio of the dressed leptonic vector and axial vector couplings g V,A of the Z boson, keeping the notation s 2 l as a short term. The effective couplings contain the lowest-order leptonic Z couplings and the corresponding loop contributions ∆g V,A evaluated at the Z-mass scale (without the QED corrections associated with virtual photons). It is convenient to relate s 2 l to the on-shell quantity s 2 W in (30) by a factor κ wich incorporates the loop corrections in terms of ∆κ, Analogously to ∆r in (25), ∆κ = ∆κ(M W , M Z , m t , . . . ) is a model-dependent quantity, depending on all parameters of the standard and non-standard particles which enter the loop contributions. An additional model dependence in the prediction of s 2 l occurs via s 2 W through the W mass obtained by means of ∆r, as described in the previous Section 3.

One-loop calculations
Loop corrections lead to deviations ∆κ from the relation κ = 1. The expansion of (60) yields the one-loop contribution which can be written as follows, The first two terms originate from photon-Z mixing and from renormalization of the on-shell weak mixing angle (31) The corrections from the vertex diagrams and the wave-function renormalization of the external fermions are collected in the vector and axial vector form factors F V,A , with lepton masses neglected.
In the alignment limit of the THDM, one can identify the one-loop contribution of the scalars h 0 , G 0 , G ± with the standard scalar contribution to the SM part of ∆κ. Consequently, the non-standard contribution to ∆κ arise from the scalars H 0 , A 0 , and H ± . Thus, one can write at the one-loop level with the SM part where the subindex indicates that only the standard particles are kept in the one-loop self-energy and vertex corrections. The contributions from the non-standard scalars to the vertex corrections and to the lepton self-energies are suppressed owing to the small Yukawa couplings. Accordingly, we can neglect the non-standard contributions to the form factors F V,A . The one-loop non-standard contribution to ∆κ is thus given by and is obtained via the non-standard content of the vector-boson self-energies from Appendix A. A sizeable part of the one-loop contribution to ∆κ is associated with ∆ρ (1) , Eq. (34), via involving the non-standard term (37), which becomes large in case of a large mass splitting in the non-standard scalar spectrum.

Higher order corrections in the SM
The one-loop calculation of the effective mixing angle was performed in [11] for the neutralcurrent vertex in neutrino scattering. The first results for s 2 l at the Z resonance [35,36] are based on one-loop calculations improved by higher orders from ∆α and ∆ρ, and an expansion in the top-quark mass [37]. QCD corrections are known at the two-loop order [12,13,[15][16][17][18], together with the leading three-loop [26][27][28] and four-loop [32][33][34] terms from the top quark. The two-loop electroweak SM contributions have been obtained in [38][39][40][41][42][43]. The leading threeloop corrections via the ρ parameter at were calculated for a massless Higgs boson in [31] and with Higgs mass dependence in [30].
In [43] a simple parametrization is given, which incorporates the complete electroweak one-and two-loop corrections together with the QCD corrections of O (αα s ) [12,13,[15][16][17][18] and O (αα s ) [26][27][28] and the leading electroweak three- 30,31]. An update is found in [44] which includes also the leading O G F m 2 t α 3 s terms [32][33][34]; the differences, however, are very small, not more than 2 · 10 −5 for s 2 l . The required input parameters are the masses of the Higgs particle, the top quark and the Z boson, together with the strong coupling constant α s M 2 Z .

Two-loop corrections in the THDM
The one-loop corrections can be supplemented by the leading reducible and irreducible twoloop contributions, containing products of ∆α and ∆ρ (1) and the irreducible two-loop correction δρ (2) to ρ, see Eq. (52). The resulting two-loop contribution to ∆κ is given by As a side-remark, we note that the appearance of the reducible products in ∆κ is again a consequence of the on-shell parameters in the normalization of ∆ρ. After a reparametrization ∆ρ (1) → ∆ρ (1) by the Fermi constant with the help of Eqs. (45) and (33) in the expression the resulting two-loop shift cancels the reducible term in (70) and we obtain as the two-loop correction to ∆κ. Next we have to identify the non-standard two-loop corrections. Separating ∆ρ (1) t and ∆ρ (1) NS in the reducible term of (70) yields Since the terms without ∆ρ (1) NS are already incorporated in ∆κ SM , we retain the non-standard reducible part ∆κ (2) NS,red = ∆α The irreducible part δρ (2) can be separated into the different non-standard contributions according to (52), which leads to the following contributions to ∆κ, • contribution from the non-standard top-Yukawa couplings ∆κ (2) t,NS = • contribution from the exclusively non-standard scalar corrections ∆κ (2) H,NS = • contribution from the mixed standard-non-standard scalar corrections ∆κ (2) In the IHDM, the only irreducible non-standard two-loop contribution to the ρ parameter is δρ (2) IHDM (which is equivalent to δρ With these specifications, the two-loop irreducible contribution to ∆κ is given by ∆κ (2) NS,irr = ∆κ (2) t,NS + ∆κ (2) H,NS + ∆κ in the aligned THDM, and ∆κ (2) NS,irr = ∆κ (2) IHDM,irr (80) in the IHDM.

Incorporation of the higher-order corrections
For an accurate prediction of s 2 l in the THDM we have to combine all the by now available loop contributions from the SM and beyond in the quantity ∆κ, which depends on the SM input quantitites and on the non-standard parameters of the THDM, For a given set of parameters, the predicted value of s 2 l fulfills the relation (63) with the corresponding M W calculated from (24). Since the W mass varies with the model parameters, ∆κ is needed as an expression where M W appears as an independent variable In the following we write only M W as an explicit variable in ∆κ, dropping the ellipsis. The SM part, contains besides the one-loop result the complete two-loop and the partial higher-order contributions listed above. Since the full SM result is an inconvenient expression, we make use of the parametrization for the SM prediction of s 2 l given in [43] for the set of SM input parameters M Z , m t , m h 0 , α s , yielding the value s 2 l,SM . By inverting (63) one obtains while M W is simultaneously fixed by the standard value M W,SM via (24) with ∆r SM . Since the W mass varies with the non-standard model parameters, we need ∆κ as an expression where M W appears as an independent variable. For the one-loop correction of ∆κ in the SM we use the explicit result (67) with M W as one of the input quantities. From (83) we obtain ∆κ The prediction of the effective leptonic mixing angle can now be written in the following way, NS + ∆κ (2) NS,red + ∆κ evaluated for M W as obtained from (24) with the full ∆r = ∆r SM + ∆r NS .
For comparison, we will also need the effective mixing angle including just the one-loop non-standard corrections, calculated from where M W is used as the corresponding input for the W -boson mass obtained from (24) with ∆r = ∆r SM + ∆r

Numerical results
For the analysis of the precision observables the following input [6] for the SM parameters is used, together with the mass of the SM-like Higgs boson m h 0 = 125.1 GeV, The shift ∆α = 0.05907 in the electromagnetic fine structure constant, ∆α = ∆α lept + ∆α had , contains the leptonic contribution up to four-loop order [133,134] ∆α lept = 0.031498, and the hadronic part extracted from experimental data with the help of dispersion relations, ∆α had = 0.027572 taken from [135]; this yields ∆α as the default value in the SM parametrizations for M W , ∆r and s 2 l given in [132] and [43].

Parameter dependence
Since the non-standard loop contribution is very sensitive to the difference between the masses of the charged and neutral scalars, we illustrate in Fig. 1 with ∆r (2) x = ∆r (2) t,NS (red line), ∆r (2) x = ∆r (2) H,NS (orange line) and ∆r (2) H,Mix (green line).
The one-loop non-standard contribution ∆r (1) NS depends only on the scalar masses and is independent t β and λ 5 . It is dominated by the non-standard correction ∆ρ (1) NS to the ρ parameter, which increases quadratically with the mass difference between charged and neutral scalars and is small for m H ± m H 0 or m H ± m A 0 . Adding the non-standard two-loop corrections leads to notable deviations. Similiar to the one-loop contribution, the mass shifts from ∆r (2) H,Mix are small for m H 0 m H ± or m A 0 m H ± and grow with the mass difference between charged and neutral scalars. The shifts from ∆r (2) H,NS are zero for m H ± = m A 0 but not for m H ± = m H 0 . These terms are additionally enhanced for large values of t β and vanish for t β = 1. The mass shifts from ∆r (2) t,NS originating from the top-Yukawa contribution to ρ are substantially smaller than those from the scalar sector. They are of some influence only for t β 1 and are suppressed for larger values of t β . Since the SM prediction of M W is below the experimental 1 σ limits, also for the current PDG world average, mass splitting between charged and neutral scalars can improve the agreement between the measurement and the theoretical prediction, in which the growing significance of the two-loop corrections is clearly visible.
Whereas for the PDG value moderate non-standard corrections are needed, a substantial upward shift of M W is required in order to reach the 1σ interval for the new CDF value, where the two-loop corrections are important: they enhance the one-loop result sizeably and hence have a considerable impact on constraining the appropriate range of the corresponding THDM parameters.
In a similar way, Fig. 2 illustrates the effect of a mass splitting between charged and neutral scalars, showing the impact of a variation of m H ± on the leptonic sin 2 θ eff (= s 2 l ), for the same parameters as in Fig. 1 to allow a direct comparison of the two figures. The masses of the neutral Higgs bosons are kept fixed, assuming equal masses of H 0 and A 0 on the left side and a small mass difference on the right hand side.    Differently to the W mass, increasing the mass splitting between neutral and charged Higgs bosons distorts the good agreement between theory and the world-average sin 2 θ eff,exp in (3) for the effective mixing angle. The values of m H ± required for getting to the CDF result for M W pull sin 2 θ eff away from the world average towards the 2 σ boundary. It is interesting to observe that instead full agreement is achieved with the SLD measurement (4) from the left-right asymmetry. In order to illustrate the influence of the additional parameters t β and λ 5 , we display in Fig. 3 regions in the m H ± -t β plane where the predicted M W is in agreement with the experimental results, either with the PDG value in (1) or with the new CDF result in (2). We present results for λ 5 = 1 (left side) and λ 5 = 6 (right side) with neutral masses m H 0 = 500 GeV and m A 0 = 550 GeV. The parameter configurations that lead to a prediction in accordance with the 1σ range of the experimental values are illustrated by the coloured band (PDG) and by the dashed band (CDF). Moreover, the areas shaded in red display the parameters which are excluded by the theoretical bounds from vacuum stability and tree-level unitarity. Around t β = 1 the contribution via δρ (2) H,NS is close to zero and the main two-loop effects arise from δρ (2) H,Mix , which are additionally enhanced for smaller values of λ 5 . Larger t β values yield additional contributions fom δρ (2) H,NS , especially for m H ± close to m H 0 when both ∆ρ (1) NS and δρ (2) H,Mix are small. For λ 5 = 1 the theoretical bounds restrict the allowed values of t β quite strongly; nevertheless, the contributions from δρ (2) H,NS can be sizeable, giving rise to the curvature of the bands. For λ 5 = 6 the theoretical constraints give room to a wider range of t β , simultaneously the influence of δρ (2) H,NS is weaker, owing to the factor (54), such that the dependence on t β becomes rather flat for larger values. For very small values of t β , on the other hand, the top-Yukawa contribution δρ (2) t,NS can become important due to its enhancement by t −2 β . However, such small values are strongly restricted from flavour physics (see for example [136]).

Inert-Higgs-Doublet Model
The predictions of M W and s 2 l in the IHDM are dominantly influenced by mass differences between the neutral and charged Higgs bosons as well. When only the one-loop corrections ∆r (1) NS and ∆κ (1) NS are considered, the results are equivalent to those in the aligned THDM, depending only on the non-standard particle masses. Also the reducible corrections from products of ∆ρ (1) NS and ∆α will be identical. Special features occur at the two-loop level. Differences arise through the irreducible corrections: contributions corresponding to δρ (2) t,NS and δρ (2) H,NS are absent in the IHDM, and the corrections result from the remaining contribution δρ (2) IHDM , which corresponds to δρ (2) H,Mix . This term contains both standard and non-standard scalars with couplings involving the parameter combination Λ 345 , Eq. (20). Hence, the two-loop corrections induce a dependence of the precision observables on Λ 345 as the only model parameter in addition to the masses entering the predictions.   In Fig. 4 we exemplify the dependence of M W and s 2 l on the mass difference between charged and neutral Higgs bosons for different values of Λ 345 , in terms of their shifts with respect to the SM predictions. In order to emphasize the effects of mass splitting we allow also the case with H ± as the lightest scalar, although in such a scenario the IHDM does not provide a dark matter candidate. For comparison, the shifts in M (1) W and sin 2 θ (1) eff , calculated with one-loop non-standard corrections only, are displayed by the blue dashed line. As expected, the two-loop corrections show a significant dependence on the value of Λ 345 which can enhance or diminish the shifts in an asymmetric way with respect to the minimum.
In Fig. 5 we visualize values in the m H ± -Λ 345 plane that lead to an agreement between the theoretical prediction and the experimental 1σ limits for M W , again for the PDG result (coulored band) and the recent CDF result (dashed band). The areas shaded in red display the theoretical constraints on the parameter space. Note that, differently to the two-loop corrections to the ρ parameter, the theoretical constraints are affected also by the choice of the coefficient Λ 2 in the scalar potential. For the selected value Λ 2 = 1 the bounds are rather loose, and the allowed regions are most significantly constrained via the W mass. The growing influence of Λ 345 with the mass splitting is visible by the curvature of the dashed band.

Computational tools
For computation of the electroweak precision observables the code THDM EWPOS has been developed, which contains all the entries described in this paper and in [88]. For the development, the calculation of the one-loop corrections was done with the help of the Mathematica packages FeynArts [137] and FormCalc [138]. The non-standard two-loop contributions were calculated using the Mathematica package TwoCalc [139,140]. The one-and two-loop results have been implemented as Fortran routines with the help of FormCalc and the methods from [141], which were developed for the implementation of the two-loop Higgs-mass corrections at O α 2 t in FeynHiggs. These routines are then called in THDM EWPOS in the calculation of the various precision observables. For the numerical evaluation of the one-loop integrals the program uses the library LoopTools [138]. For the numerical evaluation of the required two-loop integrals the Fortran routine implemented in the program FeynHiggs [142][143][144] is used.

Conclusions
We have given an overview over the currently most precise calculation of the electroweak precision observables M W and sin 2 θ eff in the aligned Two-Higgs-Doublet Model where one of the CP -even scalars (h 0 ) is identified with the scalar boson at 125 GeV observed by the LHC experiments ATLAS and CMS. It includes, besides all the known SM loop corrections, the full non-standard one loop contributions and the leading non-standard two-loop contributions from the top-Yukawa coupling and the self-couplings of the Higgs bosons. As already at the one-loop level, the non-standard two-loop contributions from the scalar self-interactions are particularly sensitive to mass splittings between neutral and charged scalars and become sizeable whenever the one-loop corrections are large. Hence, they turn out to be of special importance when a large loop-induced shift in the predicted W -boson mass is required, as it is the case for the recently published final result of the CDF experiment. Simultaneously, the predictions for sin 2 θ eff move away from the experimental world average by more than one standard deviation, towards the 2 σ limit, but show full agreement with the individual value of the SLD experiment measured via the left-right asymmetry.
As a new feature, the two-loop contributions have a significant dependence on the parameters tan β and λ 5 (or m 2 12 , respectively), a coefficient of the THDM scalar potential that is not fixed by the masses of the neutral and charged Higgs bosons. 3 In the Inert-Higgs-Doublet Model, the two-loop contributions depend on a specific combination Λ 345 of quartic coefficients in the scalar potential and are independent of tan β. In both cases, they can modify the one-loop predictions substantially and are relevant for phenomenological analyses.
A compact analytic expression for B 0 can be found in [123].