An Introduction To Quantum Field Theory - M. Peskin, D. Schroeder (Perseus, 1995) · Problem FP
Ward Identities, Renormalization of QED — M. Peskin, D. Schroeder (Perseus, 1995)
Physics — QFT · decays-of-the-higgs-boson
Problem
Decays of the Higgs Boson.
Consider the minimal form of the Glashow-Weinberg-Salam electroweak theory with one Higgs scalar field ϕ. The physical Higgs boson h0 of this theory was discussed in Section 20.2. You can use that information to compute the amplitudes for the various possible decays of the h0 as a function of its mass mh. For the other parameters of weak-interaction theory, use the following values: mb=5 GeV, mt=175 GeV, mW=80 GeV, mZ=91 GeV, sin2θw=0.23, αs(mZ)=0.12.
(a) Compute first the rate for h0→ffˉ, where f is a quark or lepton of the standard model. After a completely trivial computation, you should find
Γ(h0→ffˉ)=(8sin2θwαmh)⋅mW2mf2(1−mh24mf2)3/2⋅Nc(f),
where Nc(f)=1 for leptons, 3 for quarks. If you have studied Chapter 18, you might improve this result for the case in which the fermion f is a quark, by computing the leading-log QCD corrections for the case mh≫mq. Remember that the quark mass mq is determined at the quark threshold M2∼mq2.
(b) If mh>2mW, the Higgs boson can decay to W+W−; if it is just a bit heavier, it can also decay to Z0Z0. Compute the decay widths to these final states. You can check your result in the following way: If mh≫mW, the dominant contribution to the decay comes from production of longitudinally polarized W or Z bosons, and this contribution can be estimated as follows:
Γ(h0→W+W−)≈Γ(h0→ϕ+ϕ−),Γ(h0→Z0Z0)≈Γ(h0→ϕ3ϕ3),
where ϕ±, ϕ3 are the Goldstone bosons of the Higgs sector and the quantities on the right-hand sides of these relations are computed in the ungauged Higgs theory. Explain why this statement should be true, and verify it explicitly.
(c) The third important decay mode of the Higgs is the decay to 2 gluons. The amplitude for this decay is generated by diagrams involving quark loops. Compute these diagrams, using dimensional regularization. The diagrams will be finite, but nevertheless there is a subtle contribution which apparently depends on the regulator. Check that you have computed the amplitude correctly by verifying that it is gauge invariant. Then square the amplitude and construct the decay rate. You should find
Γ(h0→2g)=(8sin2θwαmh)⋅mW2mh2⋅9π2αs2⋅q∑I(mq2mh2)2,
where the sum runs over all quark species and I(mh2/mq2) is a form factor with the property that I(x)→1 as x→0 and I(x)→0 as x→∞. This property implies that the dominant contribution to the decay rate comes from very heavy quarks. You need not evaluate I(x) explicitly at this stage; just leave it in the form of a Feynman parameter integral.
(d) The existence of the process h0→2g implies the existence of the inverse process g+g→h0, which is a mechanism for production of Higgs bosons in proton-proton collisions. Using the parton model, derive a relation between the partial width Γ(h0→2g) and the total cross section for pp→h0+X. Compute this cross section numerically (in nanobarns) for a 30 GeV Higgs for pp collisions of center of mass energy 1--40 TeV. (1 TeV =103 GeV.) For the purposes of this problem set (though this is not actually a good approximation) you may ignore the Q2 dependence of the gluon distribution function and take simply
fg(x)=8⋅x1(1−x)7.
You may also set I(mh2/mt2)=1; this is correct to about 10\%.
(e) The final decay mode that you should consider is h0→2γ. Consider first the contribution from the loop diagrams involving quarks and leptons. Show that the result is simply expressed in terms of the form factor I(mh2/m2) that you derived in part (c).
(f) Next, compute the contribution to this decay amplitude from the loop diagram involving W bosons, and the various diagrams one must add to this to obtain a gauge-invariant result. It is easiest to work in Feynman-'t Hooft gauge. Add this contribution to that of very heavy quarks and leptons, each with electric charge Qf. Your result should reduce to the following expression in the limit mh≪mW:
Γ(h0→2γ)=(8sin2θwαmh)⋅mW2mh2⋅18π2α2⋅f∑Qf2Nc(f)−4212.
(g) Now work out the detailed behavior of the form factor I(x) defined in part (c). Reduce your expression from part (c) to a one-parameter integral, then evaluate this integral numerically. Plot I(mh2/mt2) over the range 50 GeV <mh< 500 GeV, and compute the decay width Γ(h0→2g) numerically (in keV) over this range. The computation of part (f) introduces an additional form factor; compute this function in the same way.
(h) Finally, put together all the pieces. Find the branching fraction of the h0 into each of its five major decay modes bbˉ, ttˉ, gg, W+W−, Z0Z0, for Higgs bosons of mass 50 GeV -- 500 GeV.