Arxiv:1407.4336V2 [Hep-Ph] 10 May 2016 I.Laigtrelo Orcin Otehgsmass Higgs the to Corrections Three-Loop Leading III
Total Page:16
File Type:pdf, Size:1020Kb
Higgs boson mass in the Standard Model at two-loop order and beyond Stephen P. Martin1,2 and David G. Robertson3 1Department of Physics, Northern Illinois University, DeKalb IL 60115 2Fermi National Accelerator Laboratory, P.O. Box 500, Batavia IL 60510 3Department of Physics, Otterbein University, Westerville OH 43081 We calculate the mass of the Higgs boson in the Standard Model in terms of the underlying Lagrangian parameters at complete 2-loop order with leading 3-loop corrections. A computer program implementing the results is provided. The program also computes and minimizes the Standard Model effective po- tential in Landau gauge at 2-loop order with leading 3-loop corrections. Contents I. Introduction 1 II. Higgs pole mass at 2-loop order 2 III. Leading three-loop corrections to the Higgs mass 11 IV. Computer code implementation and numerical results 13 V. Outlook 19 arXiv:1407.4336v2 [hep-ph] 10 May 2016 Appendix: Some loop integral identities 20 References 22 I. INTRODUCTION The Large Hadron Collider (LHC) has discovered [1] a Higgs scalar boson h with mass Mh near 125.5 GeV [2] and properties consistent with the predictions of the minimal Standard Model. At the present time, there are no signals or hints of other new elementary particles. In the case of supersymmetry, the limits on strongly interacting superpartners are model dependent, but typically extend to over an order of magnitude above Mh. It is therefore quite possible, if not likely, that the Standard Model with a minimal Higgs sector exists 2 as an effective theory below 1 TeV, with all other fundamental physics decoupled from it to a very good approximation. Within this model, precision calculations can help to relate observable quantities to the underlying Lagrangian parameters, as well as help to constrain new physics models, including those for which decoupling may not hold. One such observable quantity is the physical mass Mh itself. At tree level, Mh is directly proportional to the square root of the Higgs field self interaction coupling, λ. One important question has to do with the stability of the Standard Model vacuum [4–15]. The observed value of Mh is in the range that would apparently correspond to metastability of the vacuum [16–20], assuming that there is no new physics between the electroweak scale and the Planck scale. It is therefore important to pin down the relationship between λ and Mh as accurately as possible. Parametric uncertainties, notably the dependences on the top-quark mass and the QCD coupling, are not insignificant, and will likely remain so for some time. However, our attitude is that theoretical calculations should, to the extent possible, be pushed to the point that all limitations of our understanding can be reliably and unambiguously blamed on experimental error. The purpose of this paper is to present a full 2-loop calculation of the minimal Standard Model Higgs boson pole mass Mh, in terms of the MS Lagrangian parameters v,λ,yt,g,g′,g3, with the leading 3-loop corrections in the limit g ,y λ,g,g′. The relations between these 3 t ≫ parameters and other observables, such as the physical masses of the top quark and the Z and W bosons, are left to separate calculations. The result for Mh is probably too long to present as an analytical formula in print without forfeiting the goodwill of the reader, and in any case evaluation of it will necessarily rely on numerical work done by computer. We therefore present most of our results in the form of an electronic file, and as a public computer code. The computer code also performs the related task of minimizing the 2-loop effective potential [21] of the Standard Model with leading 3-loop corrections [22], implementing the form of the minimization condition given recently in [23, 24], which resummed Goldstone contributions to eliminate spurious imaginary parts and potentially infrared singular contributions. Previous work on the 2-loop contributions to the relation between λ and Mh includes the QCD corrections [18, 19], which can be obtained from the 2-loop QCD correction [25, 26] to the Higgs self-energy function. The non-QCD corrections have been obtained by [19] and [20] but were given there only in the form of simple interpolating formulas. II. HIGGS POLE MASS AT 2-LOOP ORDER To fix our conventions and notation, we write the Higgs kinetic and self-interaction La- grangian as µ 2 2 = ∂ Φ†∂ Φ Λ m Φ†Φ λ(Φ†Φ) , (2.1) L − µ − − − 3 where we use the metric with signature ( ,+,+,+), and m2 < 0, and the complex doublet − Higgs field is 1 [v + h(x)+ iG0(x)] Φ(x)= √2 . (2.2) G+(x) Here v is the Higgs vacuum expectation value (VEV), which we take to be evaluated at the minimum of the effective potential evaluated at 2-loop order with leading 3-loop corrections. This means that the sum of tadpole diagrams (including the tree-level one) vanishes at that same order, and so need not be included. Because the Landau gauge is used for the evaluation of the effective potential in [21]-[24], our calculation also is restricted to that gauge-fixing scheme. The other relevant couplings in the theory are the top-quark Yukawa coupling yt and the SU(3) SU(2) U(1) gauge couplings g , g, g′. In principle, the bottom quark c × L × Y 3 and other fermion Yukawa couplings can also be included, but they make only a very tiny difference even at 1-loop order, where their inclusion is straightforward (see below). All of 2 the couplings λ, m , yt, g3, g, g′, and the VEV v, are running parameters in the MS scheme. In order to obtain the Higgs boson physical mass Mh, we calculate the self-energy function 1 1 Π(s)= Π(1)(s) + Π(2)(s)+ ... (2.3) 16π2 (16π2)2 consisting of the sum of all 1-particle-irreducible 2-point Feynman diagrams, in the regulated theory in d =4 2ǫ dimensions. In this paper, factors of 1/(16π2)ℓ are extracted as a way − of signifying the loop order ℓ. Rather than including counterterm diagrams separately, we 2 found it more efficient to do the calculation in terms of the bare quantities λB, mB, ytB, g3B, gB, gB′ , and VEV vB, and then re-express the results in terms of the MS quantities. The complex pole squared mass is the solution of 1 1 M 2 iΓ M s = m2 +3λ v2 + Π(1)(s ) + Π(2)(s ), (2.4) h − h h ≡ pole B B B 16π2 pole (16π2)2 pole where 3-loop order effects are consistently neglected in this section. We then apply the MS relations between bare and renormalized parameters: φ φ φ 2 2ǫ 2 1 c1,1 1 c2,2 c2,1 v = µ− v 1+ + + + ... , (2.5) B 16π2 ǫ (16π2)2 ǫ2 ǫ h i 1 cλ 1 cλ cλ λ = µ2ǫ λ + 1,1 + 2,2 + 2,1 + ... , (2.6) B 16π2 ǫ (16π2)2 ǫ2 ǫ h i 4 2 2 2 1 cm 1 cm cm m2 = m2 + 1,1 + 2,2 + 2,1 + ..., (2.7) B 16π2 ǫ (16π2)2 ǫ2 ǫ yt 1 c y = µǫ y + 1,1 + ... , (2.8) tB t 16π2 ǫ g h 1 c i g = µǫ g + 1,1 + ... , (2.9) B 16π2 ǫ h g′ i ǫ 1 c1,1 g′ = µ g′ + + ... , (2.10) B 16π2 ǫ ǫ g3B = µ h[g3 + ...] , i (2.11) to obtain spole in terms of the renormalized parameters. Here µ is the regularization scale, related to the MS renormalization scale Q by 2 γE 2 Q =4πe− µ , (2.12) where γE = 0.5772 ... is the Euler-Mascheroni constant, and the counterterm coefficients are, to the orders required for this paper: φ 2 9 2 3 2 c = 3y + g + g′ , (2.13) 1,1 − t 4 4 φ 2 2 9 4 27 2 2 1 2 2 33 4 27 2 2 91 4 c = 12g y y y g y g′ g + g g′ + g′ , (2.14) 2,2 3 t − 4 t − 8 t − 8 t − 32 16 32 φ 2 2 27 4 45 2 2 85 2 2 271 4 9 2 2 431 4 2 c = 10g y + y y g y g′ + g g g′ g′ 3λ , (2.15) 2,1 − 3 t 8 t − 16 t − 48 t 64 − 32 − 192 − λ 4 2 2 9 2 3 2 9 4 3 2 2 3 4 c = 3y +6λy + 12λ λg λg′ + g + g g′ + g′ , (2.16) 1,1 − t t − 2 − 2 16 8 16 λ 2 4 2 2 45 6 27 4 2 13 4 2 9 4 2 2 135 2 2 c = 24g y 24g y λ y + y g + y g′ y λ + 108y λ y λg 2,2 3 t − 3 t − 2 t 2 t 2 t − 2 t t − 4 t 53 2 2 27 2 4 9 2 2 2 9 2 4 3 2 2 2 2 y λg′ + y g + y g g′ + y g′ + 144λ 81λ g 27λ g′ − 4 t 16 t 8 t 16 t − − 4 45 2 2 7 4 195 6 119 4 2 37 2 4 73 6 +24λg + λg g′ λg′ g g g′ + g g′ + g′ , (2.17) 4 − 4 − 64 − 64 64 64 λ 2 4 2 2 15 6 2 4 2 3 4 2 2 45 2 2 c = 8g y + 20g y λ + y y g′ y λ 36y λ + y λg 2,1 − 3 t 3 t 2 t − 3 t − 4 t − t 8 t 85 2 2 9 2 4 21 2 2 2 19 2 4 3 2 2 2 2 + y λg′ y g + y g g′ y g′ 78λ + 27λ g +9λ g′ 24 t − 16 t 8 t − 16 t − 73 4 39 2 2 629 4 305 6 289 4 2 559 2 4 379 6 λg + λg g′ + λg′ + g g g′ g g′ g′ , (2.18) −32 16 96 64 − 192 − 192 − 192 m2 2 2 9 2 3 2 c = m 3y +6λ g g′ , (2.19) 1,1 t − 4 − 4 m2 2h 2 2 9 4 2i 81 2 2 35 2 2 2 2 2 c = m 12g y + y + 36y λ y g y g′ + 54λ 27λg 9λg′ 2,2 − 3 t 4 t t − 8 t − 8 t − − 249h 4 45 2 2 55 4 + g + g g′ g′ , (2.20) 32 16 − 32 i 5 x x x x y y y z z A B I S x x y u y y z v x z z u T U M FIG.