Angular distribution and forward–backward asymmetry of the Higgs-boson decay to photon and lepton pair

The Higgs-boson decay h→γℓ+ℓ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h \rightarrow \gamma \ell ^+ \ell ^-$$\end{document} for various lepton states ℓ=(e,μ,τ)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ell = (e, \, \mu , \, \tau )$$\end{document} is analyzed. The differential decay width and forward–backward asymmetry are calculated as functions of the dilepton invariant mass in a model where the Higgs boson interacts with leptons and quarks via a mixture of scalar and pseudoscalar couplings. These couplings are partly constrained from data on the decays to leptons, h→ℓ+ℓ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h \rightarrow \ell ^+ \ell ^-$$\end{document}, and quarks h→qq¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h \rightarrow q \bar{q} $$\end{document} (where q=(c,b)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q = (c, \, b)$$\end{document}), while the Higgs couplings to the top quark are chosen from the two-photon and two-gluon decay rates. Nonzero values of the forward–backward asymmetry will manifest effects of new physics in the Higgs sector. The decay width and asymmetry integrated over the dilepton invariant mass are also presented.


Introduction
Since the discovery of the Higgs boson [1,2] its decay channels have been extensively studied. In general, the decay pattern and properties of the h boson are consistent [3] with the quantum numbers J PC = 0 ++ of the boson in the standard model (SM). Yet the nature of h needs to be clarified and will be investigated in detail in the next run of the LHC after its upgrade.
In many extensions of the SM a more complicated Higgs sector can exist, and some of the Higgs bosons may not have definite CP parity [4][5][6]. This aspect of the Higgs-boson physics is important for clarification of the origin of the CP violation, and possible additional mechanisms beyond the CP violation via the CKM matrix which can contribute to the observed matter-antimatter asymmetry in the Universe [7].
The CP properties of the Higgs boson were addressed for the two-photon decay h → γ γ in a model with veca e-mail: korchin@kipt.kharkov.ua b e-mail: koval@kipt.kharkov.ua torlike fermions [8]. It was shown that the mutual orientation of linear polarizations of the photons carries information on the CP violation. This idea was elaborated in Ref. [9], where the Bethe-Heitler conversion on nuclei of the two photons to electron-positron pairs was suggested as a means to probe the CP violation in the Higgs coupling to photons. In Ref. [10] the author analyzed possibilities of observation of the CP violation effects in the Higgs decays h → V 1 V 2 → ( f 1f2 ) ( f 3f4 ) to various final lepton and quark pairs, where V = W ± , Z .
In Refs. [11,12] the authors suggested to study the CP violation effects in the Higgs sector in the decay h → γ Z via the polarization parameters of the photon, or Z boson. A direct way for this is the forward-backward (FB) asymmetry in the decay h → γ Z → γ ff , where the Z boson on the mass shell decays to fermions. This observable vanishes in the SM and therefore it carries information on physics beyond the SM. Estimates of CP violation effects in some models of new physics were made in [11,12].
In Refs. [22,23] the authors discussed the angular distribution of the leptons = (e, μ, τ) in the decay h → γ + − in the framework of the SM. The importance of the FB asymmetry was emphasized, and its nonzero values were found. At the same time, beyond the SM, in Ref. [24] the FB asymmetry was proposed as a probe for CP-violating Higgs coupling to Z γ and γ γ states.
The FB asymmetry sure enough is an informative observable which can be of interest for future experiments at the LHC. In the present paper we address the decay h → γ + − in some detail. In addition to the loop mechanism h → γ Z * → γ + − considered in [11,12], we include here the photon bremsstrahlung off leptons, i.e. tree-level amplitudes for h → γ + − , and the loop amplitude h → γ γ * → γ + − . We remark that in the framework of the SM the FB asymmetry is equal to zero as a consequence of the scalar nature of the Higgs boson. This asymmetry can take nonzero values only in models beyond the SM and therefore this observable is sensitive to possible CP violation in the Higgs sector.
To estimate values of this asymmetry we apply a model in which the Higgs boson couples to fermions with a mixture of the scalar (S) and pseudoscalar (PS) interactions. The strength of the S and PS couplings, 1 + s f and p f , respectively, are partly constrained from the LHC measurements of the decay rates h → + − and h → qq (where q = (c, b)) [25,26]. As for the Higgs interaction with the top quark, the corresponding couplings are chosen from experimental information on the two-photon, h → γ γ , and twogluon, h → gg, decay widths.
In this model, for the decays h → γ + − we derive the distribution over the angle θ between the momentum of the lepton (in the rest frame of the pair + − ) and momentum of the photon (in the rest frame of h). The presence of the PS h ff coupling gives rise to the linear in cos θ terms in this distribution, and thereby to a FB asymmetry. We calculate the differential decay width and FB asymmetry as functions of the dilepton invariant mass squared q 2 = (q + + q − ) 2 (q + and q − are the four-momenta of leptons). The widths and FB asymmetries integrated over the invariant mass are also discussed.
The paper is organized as follows. In Sect. 2 amplitudes and angular distribution in h → γ + − are presented. The loop contributions are defined for the S and PS Higgs couplings to the fermions. The FB asymmetry is discussed. In Sect. 3 the differential decay width and FB asymmetry for various leptons are calculated. The results of the calculation are discussed. Section 4 contains the conclusions. In Appendix A the loop integrals are defined, and in Appendix B vanishing of the contribution from axial-vector Z ff coupling to the fermion-loop diagrams is shown.

Amplitudes and angular distribution
There are models with more than one Higgs doublet which induce CP violation due to the specific coupling of neutral Higgs bosons to fermions. We assume that the couplings of h boson to the fermion fields, ψ f , are given by the Lagrangian including both scalar and pseudoscalar parts, where v = ( √ 2G f ) −1/2 ≈ 246 GeV is the vacuum expectation value of the Higgs field, G f = 1.166378×10 −5 GeV −2 is the Fermi constant [27], m f is the fermion mass, and s f , p f are real parameters (s f = p f = 0 corresponds to the SM). Equation (1) can be considered as a phenomenological parametrization of effects of new physics. As for the Higgs interaction with the W ± and Z bosons, it is assumed to be the same as in the SM.
We consider the decay of the zero-spin Higgs h boson where the four-momenta of the h boson, photon, and leptons are p, k, q + , q − , respectively, and (k) is the polarization four-vector of the photon. The differential decay width can be written as where m h is the mass of the h boson, q ≡ q + + q − , q 2 is the invariant mass squared of the lepton pair, β = 1 − 4m 2 /q 2 is the lepton velocity in the rest frame of the lepton pair. The polar angle θ is defined in this frame and it is the angle between the momentum of lepton l + and the axis opposite to the direction of the Higgs-boson momentum.
The amplitude of the decay is where the tree-level amplitude ( Fig. 1) is where e = 4πα G f is the positron charge, Q = −1 (lepton charge in units of e) and m is the lepton mass. The electromagnetic coupling in the is the mass of the W (Z ) boson. For the rest we follow the standard definition of the γ matrices and lepton spinors (see, e.g. [29]).
The loop contributions h → γ γ * /Z * → γ + − (see Fig. 1) can be written in the form with coefficients c 1 , . . . , c 4 , which are specified below in terms of the loop functions, and 0123 = +1. Here we follow the notation in Refs. [22,23]. Note that we do not take into account the loop contributions of the type h → γ + − (the so-called box diagrams). The contribution of these diagrams to the considered decay is negligibly small in the SM [13,14]. Also the processes In particular, resonant production of the quarkonium states J/ψ (cc) and ϒ(1S) (bb) is of interest for studying the hqq coupling (see, for example [30][31][32][33]). The account of such mechanisms lies beyond the scope of the present work.
We evaluate the amplitude (4) squared, sum over the lepton and photon polarizations, and obtain the following result in the model (1): The fact that c 0 is real while c 1 , . . . , c 4 are generally complex-valued is used in derivation of (8).
The coefficients in Eq. (8) are defined as follows (we use the notation z ≡ cos θ ): The FB asymmetry is defined as (see, e.g. [11,12,22,23]) where As only the coefficients B, C, and F are linear in cos θ , then it is seen from Eqs. (9)-(17) that the numerator of the asymmetry (18) is determined by the imaginary part of the terms c 2 , c 4 , and the combination c 1 c * 4 + c 2 c * 3 : It may be instructive to analyze the asymmetry (18) in the limit of zero lepton masses. Putting m = 0 in (8), (11)- (17), and (20) one obtains the distribution over the + − invariant mass and the FB asymmetry Im c 1 c * 4 + c 2 c * For further reference we also introduce the integratedover-q 2 asymmetry [22]: for the appropriate integration limits q 2 min ≥ 4m 2 and q 2 max ≤ m 2 h .

Loop contributions
Let us specify the loop contributions in Fig. 1 to the coefficients c 1 , . . . , c 4 . We introduce below the Weinberg angle θ W and the notation s W ≡ sin θ W and c W ≡ cos θ W . We evaluate the loop diagrams using the Lagrangian (1) for the h ff vertex. The scalar coupling of the Higgs to fermions contributes to the coefficients c 1 , c 2 , which read where m Z ( Z ) is the mass (total decay width) of the Z boson, and Here g = 2m W ( √ 2G f ) 1/2 is the SU (2) L coupling, Q f is the charge of the fermion f in units of e, N f = 1(3) for leptons (quarks), g V, f = t 3L , f − 2Q f s 2 W and g A, f = t 3L , f are the vector and axial-vector couplings of Z boson to the fermion, where t 3L , f is the projection of the weak isospin, and The loop integrals for fermions, A f (λ f , λ f ), and W bosons, A W (λ W , λ W ), are expressed in terms of the loop functions I 1 (λ , λ) and I 2 (λ , λ) [34] (see Appendix A). The coefficients c 3 , c 4 in the amplitude (7) come only from the PS coupling of the Higgs boson to fermions in the loops. We obtain Of course, the sum over all fermions f = ( , q) in (26), (27), (31), and (32) is implied.

Forward-backward asymmetry in the SM
In the SM the angular distribution in Eq. (8) where (24) and (25).

In follows from Eq. (20) that in the SM
Therefore a nonzero value of the FB asymmetry can arise only in certain models beyond the SM. A similar conclusion for the decay h → γ Z → γ + − with on-mass-shell Z boson has been inferred in [11,12] and is used there to estimate the magnitude of possible CP violation effect.
The result (34) is at variance with the conclusion of Refs. [22,23], where the authors have found a nonzero FB asymmetry in the framework of the SM. The origin of a nonzero asymmetry in Ref. [22] is related to the axial-vector coupling of the Z boson to the fermions in the loop diagrams.
In fact, the axial-vector Z ff coupling does not contribute to the process h → γ Z * (for real or virtual Z ). This was noticed long ago in the framework of the SM in Refs. [35][36][37] on the basis of charge-conjugation parity arguments. As an alternative argument, in Appendix B we show in the model (1) and in the SM explicit cancelation of contributions from axial-vector Z ff coupling to the fermion-loop diagrams for h → γ * Z * .

Results of calculations and discussion
Let us discuss the choice of parameters s f and p f for the Higgs coupling to the fermions in (1). In terms of these parameters the decay width of the Higgs to fermions, except the top quark, is equal to where β f = 1 − 4m 2 f /m 2 h is the fermion velocity in the rest frame of h. Apparently, one can put β f ≈ 1. Then in order to keep the Higgs decay widths to fermions equal to their SM values we impose the following constraint on the parameters s f , p f : In this case, in order to ascertain the exact values of the parameters s f and p f one would need to measure polarization characteristics of the leptons, which is not accessible at present. Although Eq. (36) does not uniquely determine the parameters we choose the tentative values as in Ref. [11]: for all fermions. These values imply an equal weight of 1/2 of the S and PS couplings.
Regarding the Higgs couplings to the top quark, we will choose them by requiring that the ratios are consistent with the recent CMS results [38], This allows us to choose the following values of parameters s t and p t : With these parameters, values of μ gg and μ γ γ appear to be, respectively, 1.2 and 1.23.
As the interaction of the Higgs boson with the W ± and Z is not modified compared to the interaction in the SM, the observables in the decays h → Z Z → 4 and h → W W → ν ν , where = (e, μ), are consistent with the ATLAS and CMS data and spin-parity analyses [39,40].
Numerical values of the SM parameters are taken from [27], namely, the gauge boson masses, widths, and Z ff couplings. The quark masses are chosen according to [28,41], and sin 2 θ W = 1 − m 2 W /m 2 Z . In Fig. 2  The photon minimal energy in the Higgs-boson rest frame is taken E γ = 1 GeV in order to cut-off infrared divergence, so that q max = (m 2 h − 2m h E γ ) 1/2 ≈ m h − E γ . As is seen from Fig. 2, there is a deviation from the prediction of the SM with the chosen parameters s f , p f of new physics. Integration over the invariant mass within the interval [q min , q max ] leads to the widths shown in Table 1. The effect of new physics appears on the level 10-20%, if the invariant-mass interval lies below 30 GeV. Although the decay width in this interval is very small compared, for example, to the two-photon decay width of the Higgs boson in the SM (h → γ γ ) = 9.28 keV (Ref. [28]; see Table  A.10 therein). As one can also see from Fig. 2, for the decay h → γ e + e − the dominant contribution to the width in Table 1 comes from the loop amplitude. For the h → γ μ + μ − decay, the treelevel and loop contributions are comparable, while for the h → γ τ + τ − decay, the tree-level amplitude gives the dominant contribution.
In Fig. 3 the FB asymmetry (18) is presented as a function of q. As mentioned above, the FB asymmetry can take nonzero values only in models beyond the SM, though not all models of new physics lead to nonzero FB asymmetry. In the model (1), A FB (q 2 ) is proportional to Eq. (20). All terms in (20) are proportional to the parameters p f which characterize PS couplings of the Higgs to the fermions. So in this model the FB asymmetry is a direct measure of a possible CP violation in the h ff coupling.
For the light final leptons, e + e − and μ + μ − , the dominant contribution to A FB (q 2 ) comes from the term in (20) proportional to the imaginary part of the combination c 1 c * 4 + c 2 c * 3 . This imaginary part in turns originates from the Z -boson propagators in Eqs. (26), (27), (31), and (32), and the loop contributions Z , γ , Z and γ . The latter have small imaginary parts arising due to the intermediate on-mass-shell fermion-antifermion pairs with the masses m f ≤ m h /2. These imaginary parts come mainly from the bottom, charm quarks, and the τ lepton (this fact was also noticed in [11]).
As is seen from Fig. 3, for real values of the parameters s f , p f (see solid lines) the FB asymmetry takes values less than 1 % for the electrons and muons, with maximum value 0.8 % at the dilepton invariant mass around the Z -boson. For the τ leptons, the FB asymmetry is bigger, with a maximum value of about 2.5 %. In principle, observation of a nonzero FB asymmetry will point to CP violation in the Higgs coupling to fermions, though its small values make the corresponding experimental task difficult.
Let us emphasize that real parameters s f , p f follow from the requirement of Hermiticity of the Lagrangian L h f f in Eq. (1). Note that Hermiticity of Hamiltonian is a necessary condition, in addition to Lorentz invariance, locality, and the connection between spin and statistics, in the proof of the  (1). Solid lines correspond to parameters NP1 ( p t = +0.55), dashed lines to parameters NP2 ( p t = +0.55 i). The sign of the asymmetry is indicated at the curve C PT theorem in quantum field theory [42]. It is of interest to explore how a possible non-Hermiticity of the Lagrangian (1) will influence the FB asymmetry. Noticeable sensitivity of the FB asymmetry to non-Hermiticity of the h ff Lagrangian, in principle, can be used for testing the C PT symmetry.
For the purpose of this we change for the top quark the parameter p t from real value 0.55 to the imaginary value 0.55 i, while keeping the rest of parameters equal to their values in model NP1. This model is hereafter called NP2. Note that this choice of parameters does not affect the values of μ gg and μ γ γ calculated above.
As a result, the FB asymmetry increases substantially, up to 22 % for electron and 14 % for muon, while for τ lepton the maximal value of the asymmetry remains on the level of 2.5 % (see dashed lines in Fig. 3).
In general, A FB (q 2 ) changes sign as a function of the invariant mass, therefore the integrated FB asymmetry (23) over the whole interval of q is rather small and is not a suitable observable. In particular, in the model NP1 (NP2) the FB asymmetry integrated over the interval

Conclusions
The differential decay width and forward-backward asymmetry have been calculated for the decay of Higgs boson to the photon and lepton-antilepton pair, h → γ + − , where = (e, μ, τ). The calculations were performed in the framework of the SM and in a model of new physics, in which the Higgs boson interacts with fermions via a mixture of scalar and pseudoscalar couplings. Both the tree-level amplitudes and the one-loop h → γ Z * → γ + − and h → γ γ * → γ + − diagrams have been included.
We noted that the FB asymmetry vanishes identically in the SM. In models of new physics, which include effects of CP violation in the h ff interaction, this asymmetry takes nonzero values. The experimental study of the FB asymmetry is of interest in the search for effects of new physics in the Higgs-fermion interaction.
In numerical estimates of the decay width and FB asymmetry, the model parameters s f , p f have been chosen by requiring that the h → ff decay widths coincide with the widths in the SM for all leptons and quarks, except the top quark. For the latter the parameters s t , p t were constrained from the conditions that the rates of the h → γ γ and h → gg decays are consistent with the CMS data [38].
In the differential decay widths effects of new physics appear on the level of 10-20 %, especially at relatively small values of dilepton invariant mass 30 GeV.
As for the FB asymmetry, it takes nonzero values; however, these values are small. In particular, A FB (q 2 ) reaches 1 % for electrons and muons and 2.5 % for τ leptons in the region of invariant mass q ∼ m Z .
We have also shown that the FB asymmetry increases considerably if the parameter p t for the pseudoscalar htt coupling becomes complex. Specifically, for the imaginary value p t = 0.55 i the asymmetry rises up to 22 % for electrons and 14 % for muons in the region of invariant mass q ∼ 50-60 GeV. Hence the FB asymmetry for the e + e − and μ + μ − pairs turns out to be sensitive to the non-Hermiticity of the htt interaction Lagrangian. Since the requirement of Hermiticity underlies the proof of the C PT theorem, the FB asymmetry may also be used for testing the C PT symmetry.
For the τ leptons, the FB asymmetry is sensitive to non-Hermiticity of the htt interaction Lagrangian in the region of relatively small q ∼ 20 GeV, staying less than 1 %. Its maximal value 2.5 % remains the same with the real and imaginary parameter p t .
In our opinion experimental study of the differential decay width and FB asymmetry in the h → γ + − decays may give additional information on the couplings of the Higgs boson to fermions.
where the functions f (λ) and g(λ) can be expressed as