Scalar Quintuplet Minimal Dark Matter with Yukawa Interactions: Perturbative up to the Planck Scale

We confront the perturbativity problem in the real scalar quintuplet minimal dark matter model. In the original model, the quintuplet quartic self-coupling inevitably hits a Landau pole at a scale $\sim 10^{14}$ GeV, far below the Planck scale. In order to push up this Landau pole scale, we extend the model with a fermionic quintuplet and three fermionic singlets which couple to the scalar quintuplet via Yukawa interactions. Involving such Yukawa interactions at a scale $\sim 10^{10}$ GeV can not only keep all couplings perturbative up to the Planck scale, but can also explain the smallness of neutrino masses via the type-I seesaw mechanism. Furthermore, we identify the parameter regions favored by the condition that perturbativity and vacuum stability are both maintained up to the Planck scale.


I. INTRODUCTION
As one of the biggest mysteries of Nature, dark matter (DM) has drawn much attention from astrophysicists, cosmologists, and particle physicists. Among various guesses of the identity of the DM particle, the most extensively studied class of DM candidates is weakly interacting massive particles (WIMPs), because they can naturally explain the observed DM relic abundance via the thermal production mechanism in the early Universe [1][2][3][4]. WIMP models can be easily constructed by introducing a dark sector with electroweak SU(2) L multiplets. Introducing one nontrivial SU(2) L multiplet leads to the so-called minimal dark matter (MDM) models [5][6][7][8][9][10][11][12][13][14][15][16][17][18], which only involve the minimal content of new fields. Introducing more than one SU(2) L multiplet results in a richer phenomenology, but the models would be much more complicate .
The philosophy of the MDM models is to extend the standard model (SM) in a minimal way for involving dark matter [5]. For this purpose, a fermionic or scalar SU(2) L × U(1) Y multiplet in a representation (n, Y ) is introduced. The potential DM candidate would be the electrically neutral component that should be the lightest new state. If the dimension of the SU(2) L representation n is large enough to forbid dangerous decay operators, this neutral state would be able to play the role of the DM particle without imposing an artificial Z 2 symmetry. In other words, the stability of DM can be protected by an accidental Z 2 symmetry due to a large dimension of the SU(2) L representation and the Lorentz invariance. As proposed in the original paper [5], the minimal choices are a quintuplet (n = 5) for the fermionic case [6,8,[14][15][16] and a septuplet (n = 7) for the scalar case [9,11,12,16,17]. Adding new electroweak multiplets will push up the SU(2) L gauge coupling g 2 when it runs to high scales. Nevertheless, the MDM model with a quintuplet fermion or a septuplet scalar can keep g 2 perturbative up to ∼ 10 19 GeV [5] (∼ 10 16 − 10 17 GeV [43]) based on one-loop (two-loop) β functions. Scalar MDM models are quite different from fermionic ones, since scalars will bring in more coupling terms. Such a complexity has caused the neglect of a dangerous decay operator for the septuplet scalar model in the original consideration: the dimension-5 operator Φ 3 H † H will induce DM decay at loop level [43]. Therefore, the accidental Z 2 symmetry in the septuplet scalar model is not totally strict. Nonetheless, one can always introduce an artificial Z 2 symmetry to make the model work again, but considering n = 7 would not be special any more. In this case, discussing a triplet (n = 3) or quintuplet (n = 5) real 1 scalar multiplet with Y = 0 would be more economic. The real scalar triplet model have been well studied in Refs. [44][45][46], while the real scalar quintuplet model is less explored and will be the main topic of this paper.
In such scalar MDM models, scalar coupling terms may lead to another problem. Solutions to the renormalization group equations (RGEs) show that the scalar self-interaction couplings will go to infinity, i.e., a Landau pole (LP) will show up, at an energy scale far below the Planck scale [17,46,47]. With two independent septuplet self-interaction terms, the real scalar septuplet model hits a Landau pole at a scale around 10 8 GeV if the DM particle mass is fixed to satisfy the observed relic abundance. In our previous work [17], we attempted to push up the LP scale via introducing Yukawa couplings between the scalar septuplet and extra fermionic multiplets. A bonus of these extra fermions is to explain the smallness of neutrino masses through the type-III seesaw mechanism [48]. We found that such a setup can push up the LP scale to ∼ 10 14 GeV at best.
On the other hand, a real scalar quintuplet lives in a smaller representation and has only one independent self-interaction term. Consequently, the quintuplet couplings should evolve slower and reach a Landau pole at a higher scale. If extra fermionic multiplets are introduced, we may even push the LP scale above the Planck scale. Besides, such fermions could be used to explain the tiny neutrino masses via the type-I seesaw mechanism [49][50][51][52]. In this work, we will explore these possibilities in the real scalar quintuplet MDM model. For completeness, we will also discuss the constraint from the observed DM relic abundance, the bounds from direct and indirect detection experiments, and the stability of the electroweak vacuum.
The paper is organized as follows. In Sec. II, we introduce the quintuplet MDM model, and discuss its phenomenological constraints and the LP scale. In Sec. III, we study an extension with extra fermions for pushing up the LP scale, and discuss the constraints from perturbativity and vacuum stability. Conclusions and discussions are given in Sec. IV.

A. Model details
In the real scalar quintuplet MDM model, the dark sector only involves a real scalar quintuplet Φ with Y = 0, which can be expressed as The self-conjugate condition implies (∆ (Q) ) * = ∆ (−Q) . The real scalar ∆ (0) is a viable DM candidate. The gauge covariant derivative of Φ is where τ a are generators for the SU(2) L representation 5: Thus, the covariant kinetic term for Φ can be expanded as In order to protect the stability of ∆ (0) , we require that Φ is odd under a Z 2 symmetry, while all SM fields are even. The scalar potential is constructed by Φ and the SM Higgs doublet H. Since the operator Φ † τ a Φ vanishes due to the self-conjugation condition, the general form of the potential respecting the Z 2 symmetry is given by only five independent terms: Therefore, this model just brings in two couplings, λ 2 and λ 3 , and one mass parameter m as new free parameters. We assume that the vacuum expectation value (VEV) of the Higgs field is nonzero, while the VEV of Φ remains zero. Then the minimization of the potential implies two conditions, µ 2 < 0 and m 2 − λ 3 µ 2 /(2λ) ≥ 0. As in the SM, the VEV of the Higgs doublet is H After the Higgs field acquires a VEV, the λ 3 (H † H)(Φ † Φ) term equally contributes to the masses of all the Φ components. Therefore, at the tree level all components are degenerate with a shifted mass m 0 , given by Electroweak one-loop corrections break this degeneracy, making ∆ (1) and ∆ (2) slightly heavier than ∆ (0) . When m 0 m Z , the mass difference between ∆ (Q) and ∆ (0) is [5] m Q − m 0 = Q 2 ∆m, (2.8) where ∆m = α 2 m W sin 2 (θ W /2) 167 MeV, with α 2 ≡ g 2 2 /(4π). Vacuum stability (VS) sets a stringent constraint on the model. The philosophy is that the potential should keep bounded from below as the couplings evolving to high energies. The VS conditions can be obtained by means of the copositive criteria [53]: (2.9)

B. Experimental constraints
The observation of the DM relic abundance sets a constraint on the ∆ (0) mass m 0 . Assuming DM is thermally produced in the early Universe, its relic abundance can be expressed as [54] Ω DM h 2 1.07 where M Pl is the Planck mass, x F is the freeze-out parameter, g * is the total number of effectively relativistic degrees of freedom, and σ eff v the effective thermally averaged annihilation cross section accounting the coannihilation effect. For m 0 m h , annihilation and coannihilation into gauge and Higgs bosons in the s-wave are dominant, leading to the following result [9]: We take x F 25 and √ g * 10.33 for T ∼ O(TeV), and calculate the prediction to the relic abundance in the quintuplet MDM model. A more proper treatment is to include the Sommerfeld enhancement (SE) effect for DM annihilation, following the strategy in Ref. [6]. The annihilation channels in the dark sector are categorized by the electric charges of the two-body states. Enhancement factors can be computed in various categories via numerically solving the Schrödinger equations for the two-body states. The inclusion of the SE effect would increase the effective annihilation cross section, and hence reduce the relic abundance for fixed model parameters.

Vacuum Stability and Perturbativity
Relic Abundance(with SE Effect) Assuming thermally produced ∆ (0) particles in the early Universe fulfill the observed relic abundance, Ω DM h 2 = 0.1193 ± 0.0014 [55], the ∆ (0) mass m 0 can be constrained within a narrow range, depending on the coupling λ 3 . As shown in Fig. 1, a slightly larger m 0 is favored for a larger |λ 3 |. If the SE effect is not taken into account, the favored m 0 is ∼ 5 TeV. After including the SE effect, the favored value 2 is m 0 ∼ 7 TeV.

MAGIC and Fermi-LAT
Other constraints come from DM direct and indirect detection experiments. Direct detection is response to DM-nucleon scattering. The only tree-level diagram of the spinindependent ∆ (0) -nucleon scattering process is mediated by the Higgs boson, but its cross section is suppressed by m 2 0 and thus insignificant unless λ 3 is very large. On the other hand, according to the discussions in Ref. [58], the gauge loops lead to a DM-nucleon cross section of ∼ 2 × 10 −46 cm 2 for a fermionic quintuplet. As the gauge interactions of the scalar quintuplet is similar to a fermionic one, we may expect that the cross section for the scalar case would also be around this value 3 . Current direct detection experiments, such as PandaX-II [59] and XENON1T [60], have set an upper bound of ∼ 5 × 10 −45 cm 2 on the DM-nucleon cross section for a DM particle mass of 7 TeV. Thus, the scalar quintuplet model can evade current direct searches, but would be well tested in near future experiments.
For indirect detection of DM annihilation in space, the dominant process is ∆ (0) ∆ (0) → W + W − , whose cross section for m 0 m h without the SE effect is given by [9] But we must also involve the SE effect, which is even more significant as the velocities of Galactic DM particles are much lower than those at the freeze-out epoch. Using the 95% CL exclusion limit on DM annihilation into W + W − obtained by the combined MAGIC and Fermi-LAT observations of dwarf satellite galaxies [56], we find that indirect detection experiments have given very stringent constraints, as shown by the red regions in Fig. 1. Nonetheless, three windows around m 0 7 TeV, 8 − 10 TeV, and 13 − 14 TeV survive. As a result, the parameter regions suggested by the observation of relic abundance is still available.
C. β functions and the Landau pole scale RGE evolution of couplings are determined by β functions. In Appendix A, we list the β functions in the SM, as well as initial values of SM couplings. When the renormalization scale µ goes above a scale Λ s ∼ m 0 , the effect of the real scalar quintuplet should be involved. Hereafter we will adopt Λ s = 10 TeV, as suggested by the observed relic abundance. At oneloop level, the real scalar quintuplet MDM model introduces new terms in the β functions of the SU(2) L gauge coupling g 2 and the Higgs quartic coupling λ, while the β functions of the other SM couplings do no change. Here we list the β functions that are relevant in the following analysis: [47] 14) Note that g 1 is related to the U(1) Y gauge coupling g by g 1 ≡ 5/3g . The solution to in the beta function of λ 2 drive λ 2 going to divergence at high energy scales. At one-loop level, g 1 and y t do not enter the beta function of λ 2 , while λ 2 and λ 3 do not contribute to the beta functions of g 1 and y t . Therefore, g 1 and y t remain small at high energy scales. For an approximate analysis, we can safely neglect g 1 and y t in Eqs. (2.14) and (2.15). Thus, the RGEs for λ, λ 2 , and λ 3 become where t ≡ (4π) −2 ln(µ/Λ s ). For convenient, we define three functions f i by λ = f 1 g 2 2 , λ 2 = f 2 g 2 2 , and λ 3 = f 3 g 2 2 , and obtain the equations for them: where G(t ) = (b s 2 ) −1 ln[g 2 (t )/g 2 (0)]. In order to eliminate the linear terms, we further define functionsf 1 andf 2 by the shifts f 1 =f 1 + 1/8 and f 2 =f 2 + 69/52. Then we have Noting thatf 2 runs much faster thanf 1 and f 3 , we can simply neglect f 3 in the second equation and find its solution aŝ Settingf 2 (0) = −69/52, which is corresponding to λ 2 = 0 at µ = Λ s , we find that the maximal LP scale for λ 2 is Λ In this section, we will attempt to push up the LP scale obtained above. A lesson we can learn from the standard model is that the top Yukawa coupling gives a negative contribution to the self-coupling of the Higgs boson. As the Landau pole is induced by the self-coupling of the quintuplet, it is straightforward to introduce extra fermions with a Yukawa coupling to the quintuplet for postponing the Landau pole. Such a motivation leads to the 5-5-1 model studied below.

A. Yukawa interactions
There are three minimal ways to construct Yukawa interactions with the quintuplet scalar: introducing fermions living in (1, 0) ⊕ (5, 0), (3, 0) ⊕ (3, 0), and (4, ±1/2) ⊕ (2, ∓1/2). The first and second options have potential for explaining the tiny neutrino masses via the type-I and type-III seesaw mechanisms, respectively. In order to keep ∆ (0) stable, one of the two fermions participating the Yukawa interaction should be odd under the Z 2 symmetry. Particularly, in the second option, the two (3, 0) fermions should be different: one is Z 2 -odd, and the other one is Z 2 -even. Consequently, in order to give correct neutrino oscillation properties, we have to introduce at least one more triplet which is Z 2 -even. On the other hand, in the third option, the (4, ±1/2) representation should correspond to a Z 2 -odd Dirac fermion for avoiding anomalies [61][62][63][64]. The (2, ∓1/2) fermions can just be the SM lepton doublets which are Z 2 -even. However, this case is less interesting for us, as it cannot explain the neutrino masses.
For these reasons, in this work we only concentrate on the (1, 0) ⊕ (5, 0) case. The resulting model is dubbed as the "5-5-1" model. Minimally, we introduce a left-handed self-conjugated fermionic quintuplet Ψ L and several right-handed fermionic singlets N a,R . It is convenient to use the tensor notation for writing down the interaction terms. The tensor notation can be translated to the familiar vector notation using the following dictionaries of Φ and Ψ: (3.1) At the renormalizable level, the Yukawa interactions can be expressed as where a,L denotes the SM lepton doublets. i, j, k, l = 1, 2 are totally symmetric SU(2) L indices. a and b are family indices and at least two singlets are required for generating the realistic neutrino mixing. Eq. (3.2) respects the Z 2 symmetry, with Ψ L being Z 2 -odd and N a,R being Z 2 -even. Consequently, the new terms does not endanger stability of the scalar quintuplet, which has mass far below the fermionic one.
The above Yukawa interactions involve many new parameters. To illustrate the point, we adopt some working simplifications. In the first, we assume that there are three generations of N R , and their Yukawa couplings (y ν ) ab are close to those of the charged leptons. It means that the Yukawa coupling of the heavy neutrino is close to the τ Yukawa coupling, (y ν ) 33 y τ . Next, we assume all the new fermions, Ψ L and N a,R have the same mass M 51 , which is a characteristic scale of the 5-5-1 model. Thus, the neutrino masses given by the seesaw mechanism are ∼ m 2 τ /M 51 . According to the current cosmological constraint [65], m ν ∼ m 2 τ /M 51 0.2 eV, which implies that M 51 ∼ 10 10 GeV. Hereafter, M 51 = 10 10 GeV will be set as a benchmark scale of the 5-5-1 model. In the concrete numerical analysis, we will comment on the situation of deviation from this scale setup. Finally, in the 5-5-1 Yukawa interaction term, i.e., the second term of L yuk , we only consider the coupling to N 3,R , neglecting the other two couplings. Thus, we just need to deal with one 5-1-1 Yukawa coupling, y. This 5-5-1 Yukawa term can be expanded as where Ψ −Q = Ψ −Q,L + (Ψ +Q,L ) c and N 3 = N 3,R + (N 3,R ) c . The gauge couplings of the fermionic quintuplet Ψ are given by (3.4)
Note that the β function of g 2 becomes positive for µ > M 51 . Thus, one may worry about up to which scale g 2 can keep perturbative. By requiring α 2 = 4π, we find the non-perturbative scale of g 2 in the 5-5-1 model as where b tot 2 = b s 2 + 20/3 = 31/6. Setting Λ s = 10 TeV, we find that almost any M 51 > Λ s will give Λ NP g 2 > M Pl . Thus, g 2 would still be perturbative at the Planck scale. By solving the RGEs, we obtain the exact values of g 2 (µ) and y(µ) at µ > M 51 as .
Note that the perturbative condition y ≤ 4π will give a smaller upper bound on y(M 51 ).
As y grows at high scales, the effect of (y ν ) ab may become important. The β function of y ν with all Yukawa coupling included is where T = 3tr(y † u y u ) + 3tr(y † d y d ) + tr(y † e y e ) + tr(y † ν y ν ), and y u , y d , y e , y ν are the Yukawa coupling matrices for quarks and leptons. For simplicity we consider only one real element of y ν on the diagonal dominates, and denote itŷ ν . At high scales, y and g 2 are large, and the above equation can be approximated as (3.15) Asŷ ν contributes to the self-energy of N 3,R , it will modify β y by (3.16) Therefore, the growingŷ ν will boost the running of y, and hence the LP scale of y becomes lower. Nonetheless,ŷ ν is important only when y is very large (not far from its LP), so the LP scale would not change too much. Thus, it is still reasonable to neglect the effect of (y ν ) ab .
On the other hand, if y(M 51 ) < y c (M 51 ), i.e., F 0 > 0, as µ goes up y(µ) will increase at the beginning, and then turn its direction at some scale, and then exponentially drop down to zero. If the decreasing behavior happens at a scale lower than the Planck scale, the effect of the 5-5-1 Yukawa coupling would not be significant. In this case, the model is quite similar to the original quintuplet MDM model, where λ 2 blows up before the Planck scale. If the decrease happens at some scale higher than the Planck scale, then y might be large enough for keeping the λ 2 finite.

C. Numerical calculation
The above analysis is based on analytic calculations. Below we present the results obtained by numerically solving the RGEs.
Firstly, we investigate the impact of y(M 51 ) on the running of couplings. Fig. 2 shows the evolution of the couplings with different values of y(M 51 ) for fixing Λ s = 10 TeV, M 51 = 10 10 GeV, λ 2 (Λ s ) = 0.015, and λ 3 (Λ s ) = 0.2. We can see that the evolution behavior dramatically depends on the delicate input value of y(M 51 ).
• If y(M 51 ) = 0.662, it will be unable to slow down the growing of λ 2 , which reaches a Landau pole at a scale lower than the Planck scale.
• If y(M 51 ) = 0.666, the Landau pole scale of λ 2 will be push up to near the Planck scale.
• If y(M 51 ) = 0.667, all couplings will keep perturbative up to the Planck scale, and the VS conditions will be satisfied at the same time.    We have fixed Λ s = 10 TeV, M 51 = 10 10 GeV, λ 2 (Λ s ) = 0.015, and λ 3 (Λ s ) = 0.2. The red, blue, and green lines correspond to the evolution of λ 2 , λ 3 , and y, respectively.
• If y(M 51 ) = 0.672, although all couplings will keep perturbative up to the Planck scale, λ 2 will become negative and leads to an unstable vacuum before the Planck scale.
The above fine-tuning of y(M 51 ) is expected. In the β function of λ 2 one has to arrange a delicate cancellation between 108g 4 2 and 72y 4 at M 51 so that β λ 2 could be under well control in a wide energy region. Otherwise, the effect of Yukawa damping on λ 2 is insufficient for a slightly smaller y(M 51 ) while too much for a slightly larger y(M 51 ) that renders a negative β λ 2 too early.
Secondly, we study the parameter regions where the perturbativity and VS conditions are satisfied. We choose y(M 51 ) = 0.662, 0.666, 0.670 as three typical inputs, and preform scans in the λ 2 -λ 3 plane with Λ s = 10 TeV and M 51 = 10 10 GeV. The results are shown in Fig. 3. In the blue regions, all parameters can keep perturbative up to the Planck scale, while in the orange region, the vacuum keeps stable up to the Planck scale. The overlap regions simultaneously satisfied the perturbativity and VS conditions. As we vary the initial value of the Yukawa coupling y(M 51 ), the overlap region varies. We find that the perturbativity and VS conditions constrain λ 3 (Λ s ) within a range of 0.14 < λ 3 (Λ s ) < 0.4, and λ 2 (Λ s ) within a range of 0 < λ 2 (Λ s ) < 0.053. This favored range of λ 3 (Λ s ) is also indicated as a green strip in Fig. 1 for comparing with other phenomenological constraints.
The VS condition λ 3 + 2 √ λλ 2 ≥ 0 seems to allow a small negative λ 3 (Λ s ). However, the results in Fig. 3 exclude the whole regions with negative λ 3 (Λ s ). The reason is that the vacuum stability also requires the Higgs quartic coupling λ remaining positive when running to higher scales. In the SM, λ will turn to negative at a scale ∼ 10 9 GeV. The existence of the quintuplet scalar could change this behavior, because the λ 3 coupling has a positive contribution to the β function of λ, as shown in Eq. (2.14). If |λ 3 | is not large enough to turn around the trend of λ, however, the vacuum will be unstable at some scale. We set λ 3 (Λ s ) = −0.02, λ 2 (Λ s ) = 0.015, and y(M 51 ) = 0.666, and demonstrate the evolution of λ(µ) in Fig. 4. We can see that λ goes to negative values at scales ∼ 10 9 − 10 14 GeV. In this case, the vacuum stability cannot be ensured along the whole way to the Planck scale.
To end up this section, we would like to make a comment on the possible impacts from a different scale set for M 51 other than the benchmark scale 10 10 GeV. If M 51 lies below this scale, the Yukawa term can slow down the running of λ 2 since a lower scale, thus requiring a smaller y(M 51 ). But M 51 cannot be many orders of magnitude below, because β λ 2 is about to change the sign and render λ 2 negative even far below the Planck scale. M 51 10 10 GeV is disfavored neither, because λ 2 (M 51 ) would be too large to be stopped by a perturbative y. Anyway, one can employ a parallel analysis for other choice of M 51 , and the configuration of the viable parameter space will not change significantly as long as M 51 is not very far from the benchmark value.

IV. CONCLUSIONS AND DISCUSSIONS
Perturbativity puts a strong constraint on scalar MDM models, especially when the multiplet lives in a large SU(2) L representation. The scalar self-couplings usually reach a Landau pole at a energy scale far below the Planck scale, in spite of the initial values. There are two reasons leading to such a disaster. One is that the quadratic self-interaction Lagrangians result in terms with large coefficients in the β functions of the self-couplings. Once the self-couplings obtains a modest value, such terms will drive the self-couplings growing exponentially and soon violating the perturbativity. Another reason is that there are also significant terms contributed by other couplings in the β functions of the selfcouplings, e.g., a g 4 2 term with a large coefficient. These terms ensure that even when their initial values are very tiny, the self-couplings will quickly obtain modest values after a short journey of running.
In a previous work [17], we found that the perturbativity problem in the real septuplet scalar MDM model is quite stubborn: even after some fermions with Yukawa interactions are introduced to slow down the running of scalar self-couplings, the model is still unable to keep perturbative up to the Planck scale. Nonetheless, this may be achievable if the MDM scalar lives in a smaller representation. Therefore, we have studied the real scalar quintuplet MDM model in this work. The observed relic abundance implies that the scale for introducing such a quintuplet should be near 10 TeV. Our calculation suggests that the quintuplet self-coupling λ 2 will hit a Landau pole at a scale of 5.6 × 10 14 GeV, which is consistent with Ref. [47].
In order to push up this LP scale, we have extended the model with Yukawa couplings of the scalar quintuplet to a fermionic quintuplet and three fermionic singlets, resulting in the so-called 5-5-1 model. The new singlets can also play the role of right-handed neutrinos, explaining the smallness of neutrino masses by the type-I seesaw mechanism. We have found that if such Yukawa couplings are involved after a scale of M 51 ∼ 10 10 GeV, all couplings can remain perturbative up to the Planck scale, The reason is that the Yukawa couplings contribute a large negative term to the β function of λ 2 and hence slow down the growing of λ 2 at high scales.
We have also investigated the parameter regions favored by the perturbativity and vacuum stability conditions up to the Planck scale. It has been found that these conditions constrain the Higgs-quintuplet coupling λ 3 (Λ s ) at the quintuplet scale Λ s within a range of 0.14 < λ 3 (Λ s ) < 0.4, and the quintuplet self-coupling λ 2 (Λ s ) within a range of 0 < λ 2 (Λ s ) < 0.53. As most of the Yukawa couplings are negligible, only the top Yukawa coupling y t is considered in the above expressions.
In the RGE calculation, we use the following MS values at m Z as initial values for gauge couplings: [66]  The measured values of y t and λ are obtained from the pole masses of the top quark and the Higgs boson, m t and m h , respectively. Therefore, we need to derive their MS values y t (m Z ) and λ(m Z ) at m Z by the matching conditions [67]