- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Advances in High Energy Physics

Volume 2013 (2013), Article ID 631841, 13 pages

http://dx.doi.org/10.1155/2013/631841

## Anatomy of the Higgs Boson Decay into Two Photons in the Unitary Gauge

Division of Theoretical Physics, University of Ioannina, 45110 Ioannina, Greece

Received 10 December 2012; Revised 4 April 2013; Accepted 17 May 2013

Academic Editor: Alexey Petrov

Copyright © 2013 Athanasios Dedes and Kristaq Suxho. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We review and clarify computational issues about the -gauge boson one-loop contribution to the decay amplitude, in the unitary gauge and in the Standard Model. We find that highly divergent integrals depend upon the choice of shifting momenta with arbitrary vectors. One particular combination of these arbitrary vectors reduces the superficial divergency down to a logarithmic one. The remaining ambiguity is then fixed by exploiting gauge invariance and the Goldstone Boson Equivalence Theorem. Our method is strictly realised in four dimensions. The result for the amplitude agrees with the “famous” one obtained using dimensional regularisation (DR) in the limit , where is the number of spatial dimensions in Euclidean space. At the exact equality , a three-sphere surface term appears that renders the Ward Identities and the equivalence theorem inconsistent. We also examined a recently proposed four-dimensional regularisation scheme and found agreement with the DR outcome.

#### 1. Introduction

Today one of the main focal points at the Large Hadron Collider (LHC) is to search for the Higgs boson () [1–3] through its decay into two photons, (for reviews see [4, 5]). Indeed, the recent [6, 7] observation by ATLAS and CMS experiments of a resonance, which could be the Standard Model (SM) Higgs boson, is based on data mainly driven by . In the SM [8–10], this particular decay process goes through loop-induced diagrams involving either charged fermions or -gauge bosons. Their calculation was first performed in [11] in the limit of light Higgs mass , using dimensional regularisation in the ’t Hooft-Feynman gauge. Since then, there are numerous works spent on improving this calculation including finite Higgs mass effects in linear and nonlinear gauges [12–14], different regularisation schemes [15–18], and/or different gauge choices [19].

The amplitude is originated, in broken (unbroken) phase, by a dimension-5 (dimension-6) SM gauge invariant operator(s) and, therefore, its expression, within a renormalizable theory, must be finite, gauge invariant and independent of any gauge choice. The amplitude should also be consistent with the Goldstone Boson Equivalence Theorem (GBET) [20–22] since the SM is a renormalizable, spontaneously broken, gauge field theory.

A problem arises when the -gauge boson contribution (see Figure 1) to produces “infinite” results at intermediate steps. These problems are usually treated by using a gauge invariant regulator method, for example, dimensional regularization. In the unitary gauge [23], this indeterminacy is more pronounced and more difficult to handle with (however, using DR and unitary gauge with modern computer algorithms, this may not be a hard problem today [19]) due to the particular form of the -gauge boson propagator. On the other hand it is much simpler to work with only few diagrams that involve physical particle masses, rather than many.

More specifically, in the unitary gauge, one encounters divergencies up to the sixth power. It is well known that, in four dimensions, shifting momenta in integrals that are more than logarithmically divergent is a “tricky business”—recall the calculation of linearly divergent fermion triangles in chiral anomalies [24, 25]—that requires keeping track of several “surface” terms for these integrals. There is also the situation we face here where apparent logarithmically divergent integrals turn out to be finite but discontinuous at .

We would like to bypass those ambiguities and at the same time to present a “regularisation” method, by performing the calculation for the amplitude strictly in 4 dimensions and in the physical unitary gauge. Our method is similar to the one used elsewhere for calculating triple gauge boson amplitudes [26, 27], or Lorentz noninvariant amplitudes [28], and consists of three steps:(1)We write down the most general Lorentz invariant amplitude.(2)We introduce arbitrary vectors that account for the “shifting momentum” indeterminacy. We show that a particular choice of those “shifting vectors” cancels higher powers of infinities leaving still behind at most logarithmically divergent integrals that are treated as undetermined variables.(3)We exploit physics, that is, gauge invariance (Ward Identities) and the GBET, in order to fix the last undetermined variables.This method is quite general within a renormalizable theory and can be applied to other observables too. Following these steps we arrive at the same result for the amplitude obtained by Ellis et al. [11] and by Shifman et al. [13] almost 35 years ago. Our analysis, among other issues, highlights that the recent observation [6, 7] of the at the LHC signifies the validity of the Goldstone Boson Equivalence Theorem. As a further clarification we also make a remark on the direct calculation in the following three cases: we first perform the integrals in exactly (with no regularisation method beyond the one discussed in point 2 above), second, by exploiting Dimensional Regularisation (DR) as defined in [29, 30] and then taking the limit , and finally third by using a four-dimensional regularization scheme introduced in [31].

Our calculation is complementary to, but somewhat different from, the two existing ones [19, 32, 33] performed in the unitary gauge. It is not our intend to redo the calculation in unitary gauge with DR as in [19]. On the contrary, we want to clarify subtle issues related to the amplitude in unitary gauge and raised in part by [32, 33]. We find, using arbitrary vectors, that divergencies (up to the 6th power) are reduced down to logarithmic ones. This is a new result that is not obvious when working in unitary gauge and cannot be seen when using dimensional regularisation. This fact was stated incorrectly in [32, 33].

The outline of the paper is as follows: in Section 2 we present the calculation of the -loop contribution (note that the calculation of the fermion triangle contribution is well defined; that is, it is independent of arbitrary vectors and finite. We are not going to repeat this calculation here and refer the reader to the reviews in [4, 5]) to amplitude, its ambiguities, and the resolution within physics arising from GBET. Next, in Section 3 we examine details of the amplitude calculation within an alternative and recently proposed four-dimensional regularization scheme [31], the one that resembles most closely the symmetry approach taken here. In Section 4 we discuss other possible physical setups plus experiment that may help to resolve inconsistencies. We conclude in Section 5. There are two appendices: Appendix A that contains some intermediate formulae and Appendix B where we present the details about the origin of surface terms in four dimensions.

#### 2. The -Loop Contribution to in SM

The most general, Lorentz and CP-invariant, form of the of-shell amplitude is, where and are the outgoing photon momenta shown in Figure 1 and the coefficients are scalar functions of , and . By considering that all particles are on-mass-shell, that is, , we obtain an amplitude with only two undetermined (for the time being) coefficients In unitary gauge, the Feynman diagrams that contribute to and are displayed in Figure 1. In order to calculate them, we introduce three arbitrary four-vectors , and , one for each diagram. These vectors shift the integration momentum, that is, for the first diagram, for the second diagram, and for the third diagram. As we will see, these arbitrary vectors operate as regulators capable to handle highly divergent integrals related to unitary gauge choice. Furthermore, the vectors , , and linearly depend upon the external momenta and . Hence , , and are not linearly independent (cf. (4)). This is an important fact leading to the cancellation of infinities.

We first calculate the less divergent part of in (2) which is the coefficient. (The coefficient will be fixed later on by the requirement of gauge invariance.) By naive power counting, we see that diverges by at most four powers. Then we perform the Feynman integral calculations strictly in 4 dimensions. For reasons that will become clear later, we shall keep the number of dimensions general in all intermediate steps of the calculation, that is, . As we will see, contributes only to finite pieces of (cf. (9)). (On the contrary, we shall see that there are nontrivial -contributions into -coefficient.)

With all the above definitions, we can write down the total amplitude in the form where the coefficients with are given in Appendix A, and the sign is the proportionality factor: . Note that is a (superficially) 6th power divergent amplitude in the unitary gauge. in (3) solely contributes to while all other -elements contribute to both and/or in (2).

First we focus on the calculation of the “less divergent” coefficient of (2). Based on naive power counting, we observe that the terms in (3) lead to and the most quartic divergent integrals. However, when adding all these pieces together, we find that quartic divergent integrals vanish for every arbitrary vectors , and leaving behind an expression with integrals of third power (in momenta) plus integrals with smaller divergencies. Then the cubically divergent integrals are proportional to all possible Lorentz invariant combinations like , , and . Therefore, choosing
we ensure that third-order divergent integrals related to terms vanish identically. In the same way, by naive power counting, and terms—these terms in (3) together with contribute solely to —lead again to at most third-order divergent integrals. However, in the sum of and terms in (3), third-order divergent integrals vanish for arbitrary , and , leading to an expression that when added to term consists of at most quadratically divergent integrals, proportional to and . We choose
for the quadratically divergent integrals to vanish. Likewise, when we add and terms—these terms, together with , solely contribute to in (2)—the third-order divergent integrals vanish for every choice of leading to an expression that when added to consists of at most quadratically divergent integrals proportional to and . Therefore, we choose
for infinities to vanish identically. From (4), (5), and (6) we arrive at the final relation among the three introduced vectors:
Equation (7) suggests that the rest of the divergent integrals depend on, at most, one arbitrary vector, say the -vector. Note that contributes only to the finite part of . Now, if we impose conditions (7) onto the remaining expressions for terms of (3), we find that *all* quadratically and linearly divergent integrals vanish, independently of the direction of the -vector. We stress here on the fact that *the cancellation of divergencies down to logarithmic ones is a highly non-trivial, almost “miraculous” result. These cancellations only take place for a particular choice of the momentum-variable shift vectors (7).* (As a corollary, if for instance, we had split the -vertex into three pieces, each one associated with three different arbitrary vectors, then the generalised condition (7) would again downgrade the divergency of the amplitude to a logarithmic one.) Of course this is an expected outcome for an observable in a renormalizable theory.

Our final result contains *at most logarithmically divergent integrals*. (This result is different from the one obtained in [32, 33], where there are remaining quadratically and linearly divergent terms. Following (11) in Gastmans et al. paper [32] or (3.36) in their sequel paper [33], and unless there is a typo in both their formulae, we find that there is a missing quadratically divergent term proportional to . Even in the case that this is a typo, their formulae contain a linearly divergent term that is also referred to by the authors claiming that this last term reduces to a logarithmically divergent integral by changing the internal momentum and further manipulating the integral. Note here that in our calculation all divergent integrals are reduced to at most logarithmically divergent ones *without* any further manipulation nor any assumption other than (7) for every arbitrary vector: , , and . Finally, the authors in [32, 33] have made a specific choice for routing the internal loop momentum corresponding to , and therefore their result should be the same as ours.) Despite of the fact that the resulting expressions so far contain the shift instead of with an arbitrary vector , its presence is irrelevant since logarithmically divergent integrals are momentum-variable shift independent [34]. Summing up all the above contributions to , we find a particularly nice and symmetric form for :Introducing Feynman parameters, shifting momentum variable from to , and ignoring all terms (these terms will be used later in arriving at (24)) that contribute to , we find that the contribution to in (2) arises solely from the term
with . Obviously, the first integral in (9) is (superficially) logarithmically divergent while the second one is finite. The number of dimensions appears only at the finite integral, and therefore we can fearlessly set everywhere. This means that we do not use dimensional regularisation in what follows (see however the discussion below). We state here few additional remarks to be exploited later on: (a) we observe that the top line in the integrand of (9) *does not vanish* in the limit and (b) despite of the appearances in (8), there is no in the numerators of the subsequent expression equation (9). The whole contribution arises from the denominator’s -term.

Our next step is to parametrize the logarithmically divergent integral in (9) by an unknown, dimensionless parameter to be determined later by a physical argument. So we defineA parenthesis must be added at this point. We could of course promote and use dimensional regularisation [29] by exploiting symmetric integration in -dimensions. In this case, and after taking the limit , one finds which is finite and nonzero and agrees with the one we find below in (20) after imposing the GBET condition. This is also the result found in the original references [11–13]. However, according to [32, 33], the integral in (10) is discontinuous at ; in fact, when symmetric integration in is used, one finds instead . This is also understood in a slightly different context. It has long been known [34–36] that shifts of integration variables in linearly (and above) divergent integrals are accompanied by “surface” terms that appear only in four dimensions, a famous example being the integrals in chiral anomaly triangle graphs. For our purpose here let us start with the following shift of variables in a linearly divergent integral that has been generalised [36] to work in -dimensions following the expression
that is valid for and constant, possibly dependent on Feynman parameters, like the one given below (9), and is an arbitrary constant four-vector. By taking the derivative, of both sides in (11) shifting the integration variable for the logarithmically divergent integral encountered, and evaluating the finite one, we easily arrive at
For an alternative and detailed proof of (12), see Appendix B. (The same result is obtained by standard algebraic tricks. We would like to thank R. Jackiw for communicating his calculation to us.) Applying (12) to terms of (10) with , we find
This is consistent with the symmetric integration in 4 dimensions () *but*, is also consistent with the usual tabulated textbook result [37] from dimensional regularisation in -dimensions (). Equation (13) shows that is discontinuous at . Then the question arises: *which ** to believe in?* Answer: *the one that is indicated by well defined, calculable, boundary conditions and symmetries of the underlying theory*.

The above parenthesis to our calculation motivates us to avoid the direct calculation of integral (10), but set everywhere and treat as an unknown parameter to be defined later within a physical context or experiment. Substituting (10) into (9), we arrive at Evaluating the double finite integral in (14) and restoring the proportionality factor given below (3), we obtain where

Our final step is to determine the unknown parameter in (15). For this we need physics that reproduces in a different and unambiguous way. One choice, probably not the only one, is to adopt the Goldstone Boson Equivalence Theorem (GBET) [20–22] which states that the amplitude for emission or absorption of a longitudinally polarised at high energy becomes equivalent to the emission or absorption of the Goldstone boson that was eaten. Mathematically, this is written by an [38]
which says that the -matrix elements for the scattering of the *physical* longitudinal vector bosons with other physical particles are the same as the -matrix elements of the theory where the ’s have been replaced by *physical* Goldstone bosons . We are not going to get into details here; apart from the original literature the reader is also referred to the articles [19, 38–40]. Following [38], within perturbation theory and in the limit of high energies, , GBET can be expressed with physics in two different limits of the theory: (a) or (b) .

The limit (b) is irrelevant (the limit (b) simply says that matrix elements for the theory which contains the physical ’s and zero v.e.v are equal to those produced by scattering of massless physical Goldstone bosons (instead of 's) at high energies. We have checked that (17) is satisfied in this limit) for defining in (15) so we completely focus on the limit (a). It is very easy to see that, in the unitary gauge, the ’s do not decouple (this is another advantage of calculating in the unitary gauge. Note that here the word “decoupling” is stated in a different context than is usually quoted following the Appelquist-Carazzone [41] theorem, where particle masses circulating in loops are much heavier than the external ones) for vanishing gauge coupling . Consider for example the diagrams in Figure 1: there is always an from the -vertex that cancels another sitting in the denominator of the longitudinal part for the internal -boson propagator expression written in the unitary gauge. So, as it was already noted in the paragraph below (9), in the limit , there are remaining nondecoupled terms. Unfortunately, these effects may be obscured or misjudged by the regularisation method needed to handle divergent, intermediate loop integrals. This is exactly what happens here when trying to calculate directly from its ambiguous form (10). On the other hand, however, at the exact , with fixed v.e.v and Higgs quartic coupling , (17) suggests that the Goldstone bosons () should reappear at the physical spectrum of the theory while the longitudinal components of ’s become unphysical. At this limit, the SM is a spontaneously broken global -symmetry that couples, minimally, to electromagnetism. The interactions between the Higgs and photon with the Goldstone bosons are simply those of a spontaneously broken scalar QED with :

Armed with these Feynman rules we calculate the diagrams in Figure 2. By doing so, we introduce again three momentum variable shift vectors, one for each diagram, exactly in the same way we did for the calculation of the diagrams in Figure 1. The Lorentz structure of the amplitude is completely analogous to (2) with , but now due to the scalar propagators, the superficial degree of divergence, for diagrams contributing to , is . Hence, all integrals involved in are finite and in addition, they are independent of any momentum integration shift vector variable. As a consequence, is well defined, calculable, independent of any regularisation method, and at the limit of (or ) is By equating (15) (in the limit ) and (19) which represent the l.h.s and r.h.s of the GBET condition (17), respectively, we find This value agrees with dimensional regularization [29, 30] in the limit (see (13)). The final form of the in (2) is with and defined in (16).

To complete the picture there is still the coefficient in (2) to be calculated. Naive power counting says that this is by two powers more divergent than and, in general, undetermined. It can be fixed however by using quantum gauge invariance, that is, conservation of charge, for the : and thus from (2), Equation (23) is substituted to (2) with read by (21). This is exactly the same result for the -boson contribution to amplitude that has been obtained in [11–14, 19] using dimensional regularisation in -gauges.

It is interesting here to note the result from the explicit algebraic manipulation of in the unitary gauge and check the validity of gauge invariance (see (23)). Exactly as for , the condition for the arbitrary vectors given in (7) is crucial in reducing the divergence of down to a logarithmic one (see expression equation ). In -dimensions the expression for is finally independent of any arbitrary vector and, up to a proportionality factor, readsClearly the first integral in (24) is ill-defined in four dimensions. If, however, we insist on doing the calculation of (24) in with symmetric integration, like in [32, 33], we find that gauge invariance (see (23)) is *not* satisfied. This is of course unacceptable. By going a little bit deeper, gauge invariance is lost because of the term proportional to in (24) which vanishes when . Quite the contrary in DR, this term results in a non-zero contribution when since the (log divergent) integral in front of contains a simple pole at . This changes the final result and renders (21), (24), and (23) consistent only if . This outcome is in agreement with [19].

Few remarks are worth mentioning here. Had we started first calculating , there would be no possibility of defining unambiguously without using a gauge invariant regulator: the part of the amplitude at involving Goldstone bosons (see diagrams Figure 2) is not well defined—an integral as the one in (12) appears again. Another remark is that the same expressions for the coefficients displayed in Appendix A in the unitary gauge appear also when one exploits the -gauge. In the latter there are in addition -dependent terms [19] that vanish in the end from unphysical scalar contributions. Therefore, the logarithmic ambiguity in (12), found here in the unitary gauge, is similar in every other gauge.

The ambiguity of the integral in (12) has been discussed by many articles in the recent literature. References [15, 16] have used gauge-invariant regulators alá Pauli-Villars. Most notably, Piccinini et al. [18] showed that the unitary gauge with a cut-off regularisation scheme turns out to be nonpredictive: new physics input, along the lines of [28] also followed in our article, is needed. As already noted, Marciano et al. [19] were the first to make the calculation of in unitary gauge with DR. Furthermore, [15, 18, 19, 39, 40] showed that the “decoupling limit” must hold because of the GBET. All these articles together with the one at hand conclude that the result in [32, 33] is incorrect. Furthermore, the aim of our paper is not to redo the calculation in the unitary gauge with DR as in [19]; this is only a by-product of our analysis. On the contrary, we want to clarify subtle issues related to this calculation in the unitary gauge and in raised in part by [32, 33]. We find by working strictly in , using arbitrary vectors, that divergencies (up to the 6th power) are reduced down to logarithmic ones. This is a new result that is not obvious when working in the unitary gauge and cannot be seen when using dimensional regularisation.

#### 3. Four-Dimensional Regularization (FDR)

So far we have proposed a regularization scheme which is four-dimensional and uses the basic symmetries and underlying physics of the SM. However, in more complicated models or observables with more parameters to adjust, such a scheme can become cumbersome. For example, it is not always obvious which physics argument will fix undefined integrals.

Very recently, Pittau [31] proposed a scheme that is fairly easy to handle and, to the best of our knowledge, is the closest to four-dimensional calculations and thereby coined four-dimensional regularisation/renormalization scheme or just FDR. According to this scheme, infinite bubble graph contributions, that is, large loop momenta contributions that do not depend upon external momenta, are absorbed into the shift of the vacuum while the remaining finite corrections are calculable in four dimensions in addition to being Lorentz and gauge invariant.

We have applied FDR to the calculation of the amplitude and found agreement with our physics approach and with DR results. In FDR one introduces an arbitrary scale which is considered to be much smaller than internal momenta and particle masses in loops. Self-contracted loop momenta quantities like become , while, for gauge invariance to hold, vector momenta, , remain untouched. For example the integral of (12) becomes
where and stands for integration over , dropping all divergent terms from the integrand (see below) and taking the limit . In going from l.h.s to r.h.s of (25) the symmetry property has been used in four dimensions. Then, using the partial fractions identity,
the term in square bracket is recognised as divergent and therefore removed, and integrating the r.h.s of (25) over one obtains
that is, exactly the same result as in DR which eventually leads to consistent with gauge invariance and GBET. What in fact FDR scheme does is to restate the correct DR answer through the regulator keeping (12) correct in . We therefore understand that the constant (-independent) term of (21) in FDR arises from the fact that the arbitrary scale, , *must* disappear from physical observables.

#### 4. Discussion

It is evident that our calculation for the amplitude incorporates two physical inputs: one is the conservation of charge and the other is the equivalence theorem. They are both direct consequences of the gauge invariance of the underlying physical theory. The first one is experimentally indisputable while the second one is theoretical (this is not entirely correct. There is of course the high energy behaviour of found at LEP [42] consistent with the GBET.) and has been proven in [43], which is valid in any spontaneously broken renormalizable theory, like, for example, the SM. One may think however that there is a loophole in our use of this second argument: so far, and, to our knowledge, the replacement of the -bosons with Goldstone bosons at high energies has been proven to be valid only for external -bosons [38, 44, 45] and *not* for internal ones which is the case exploited here. Although it has been tested in several phenomenological examples [46], a formal, to all orders, proof is still missing. Although this may be true, it is difficult to argue against the validity of decoupling limit (with fixed v.e.v and Higgs quartic coupling) discussed in the paragraph above (18).

Is there another physics context from which one can define ? One possibility is to exploit the low-energy Higgs theorem [11, 47–49] instead. Although this may serve as a consistency check and indeed is compatible with , we cannot use it to define . The reason here is threefold: first, when treating the Higgs field as an external background field with zero momentum one needs to take partial derivative w.r.t of the 2-point photon vacuum polarization amplitude, . The latter is notoriously difficult, if meaningful, to be calculated in the unitary gauge. Second, according to [11], we know that to the lowest order in weak coupling, the amplitude for the process is proportional to , where is the improved energy momentum tensor [50]. However, the calculation of goes through the same steps as for the calculation for the amplitude and therefore involves the same ambiguity for calculating . Third, one could examine the -contribution to within the dispersion relation approach. It can be shown [51] that the nonvanishing limit at is due to a finite subtraction induced by the corresponding trace anomaly [52]. However, in order to calculate unambiguously this finite piece, one has to make full use of a (physical) boundary condition of the theory.

As a final remark, suppose that we did not know DR and wanted to calculate a certain observable in 4 dimensions. In this observable we encounter singularities, that is, undefined and undetermined integrals. Then we use physics arguments to fix these ambiguities. However, we can always question whether we are using the right physics setup or not. In that sense the final judgement should come from the experiment. Therefore it may be not academic to ask whether LHC could see the difference between and . Setting the SM Higgs mass GeV and including the top-loop contribution, we find
This is certainly within LHC’s sensitivity for 14 TeV c.m energy and luminosity of fb^{−1} (see, e.g., Figure 3 in [53]). In fact, the recent observation by LHC experiments [6, 7] indicates a value [54] which highly disfavours the case by almost four standard deviations. We can turn this around and state that this is an indirect hint towards the validity of the equivalence theorem.

#### 5. Conclusions

In this work we review the -gauge boson loop contribution to the amplitude in the unitary gauge. Our objective is to fix intermediate step indeterminacies arising from divergent diagrams by making full use of physics at much in the same way as in the calculation of the chiral anomaly triangle.

We anticipate a finite result for the loop-induced -amplitude in the renormalizable SM. Therefore the amplitude has to be independent of any shifting momentum variables we have originally introduced. But finite or even log divergent integrals are independent of these vectors, so the vectors have to be accompanied only by infinite contributions, if at all. Therefore, infinities and arbitrary vectors are eliminated altogether by a certain combination among them (see (7)).

The whole calculation in the unitary gauge boils down to a logarithmically divergent integral (10). We find that this integral results in two different values depending on whether or . This is due to a surface term remaining at the exact case after the part-by-part integration in -dimensions (see Appendix B). To proceed, we identify this integral with an undefined parameter (see (12)). This parameter is then fixed unambiguously by assuming the validity of the Goldstone Boson Equivalence Theorem (GBET). Its value is consistent with DR in the limit .

In our calculation we are very careful not to perform shifting of integration variables for highly divergent integrals by introducing three arbitrary momentum variable shift vectors straight from the beginning. Divergencies and arbitrariness from these unknown vectors are altogether removed, leaving behind a log-like divergent integral in of (2). This is defined by a physical input taken from the GBET and is connected to by electromagnetic charge conservation.

As noted many times in the text, the key point towards deriving an unambiguous amplitude for in the unitary gauge is the limit of vanishing gauge couplings; this is an aspect of GBET (see (17)). In this limit, the Goldstone boson loop contributions to the coefficient finite, that is, independent of any regularisation scheme.

We also saw that DR (FDR), a regularisation scheme introduced to maintain Ward Identities at intermediate steps of a calculation, supports the GBET in the limit (). On the contrary, we find that performing the integrals in with symmetric integration is not a good choice because it leads to the violation of gauge invariance (see (23) and the discussion below). The main reason is due to surface terms that are developed in exactly dimensions (see discussion below (10) and Appendix B). The latter are axiomatically discarded in DR [29, 30]. Another reason is the appearance of the -term in the numerator of (24).

In conclusion, the four-dimensional calculation of amplitude in the unitary gauge is ambiguous without introduction of a physics input beyond gauge invariance. As we have demonstrated, this physics, which uniquely defines the amplitude, may arise from the Goldstone Boson Equivalence Theorem (GBET). This effectively proves that GBET comprises an additional important pillar of the Standard Model dynamics.

#### Appendices

#### A. The Coefficients

We append here the integrand expressions for the coefficients in (3). The corresponding formula for is quite long and is not included here. It can be provided by the authors upon request. Note that the number of dimensions has been kept arbitrary throughout, and on-shell conditions for the external particles have been imposed:
It is straightforward, but long and tedious, to show that after implementing the condition (7) to coefficients in (A.1) we arrive at (8) which is *at most logarithmically divergent. *

For complementarity reasons, it is useful in deriving (24) to present the expression for the coefficient after the imposition of the arbitrary vector relation equation (7):
This integrand expression, under , is obviously *at most logarithmically divergent*.

#### B. Dimensional Regularization and the Surface Term

We would like to examine the surface terms arising in when calculating the integral on the l.h.s of (12). This integral after Wick rotation into Euclidean space reads where and drop for clarity the subscript from now on. We follow very closely’t Hooft and Veltman’s seminal paper in [29]. In our calculation for a physical process we should notice first that are strictly 4-vectors since they are contracted with physical external momenta or . On the other hand, the loop momentum in has components in all, , dimensions. We write as a sum of a vector which has non-zero components in dimensions 0, 1, 2, and 3 and a vector which has nonzero components in -dimensions: With this definition, the integral (B.1) reduces to where the components in the numerator of (B.1) vanish thanks to symmetric integration formula, . In order not to carry the in all formulae below we just concentrate on the integral Integrating over the extra-dimensional solid angle we arrive at where is the Euler -function and is the length of the vector. This integral is UV divergent for and IR divergent for . Therefore, the region of convergence, , is finite, but it does not yet include the point . In order to enlarge the region of convergence to include one has to analytically continue by inserting the identity in (B.5). After integrating by parts in the region of convergence, rewriting the r.h.s in terms of from (B.5), and keeping only, potentially, nonvanishing surface terms, we arrive at where the first integral is over the Euclidean spatial components of a 4-vector on a three-sphere. The surface integral converges in while the other in . By taking the surface integral on a three-sphere with radius and eventually taking the limit we find which now converges in the region that is, it includes the point . For this surface term vanishes while for there is a finite piece, , remaining. This is exactly the term that spoils gauge invariance and the equivalence theorem. In DR this term is axiomatically absent—allowing the shifting of integral momenta is among DR’s main properties.

Turning into the second integral of (B.7) we note first that the region of convergence includes now . It gives By placing (B.8) and (B.9) into (B.7) we finally arrive at (12).

#### Acknowledgments

The authors are grateful to Adrian Signer for illuminating discussions and critical comments on the manuscript. The authors would also like to thank R. Pittau, R. Jackiw, F. Piccinini, P. Kanti, S. Martin, and K. Tamvakis for discussions and comments. The authors also thank C. Coriano for bringing to our attention [51]. This research Project is cofinanced by the European Union—European Social Fund (ESF) and National Sources, in the framework of the program “ARISTEIA” of the “Operational Program Education and Lifelong Learning” of the National Strategic Reference Framework (NSRF) 2007–2013. Kristaq Suxho acknowledges full the financial support from the Greek State Scholarships Foundation (I.K.Y).

#### References

- P. W. Higgs, “Broken symmetries and the masses of gauge bosons,”
*Physical Review Letters*, vol. 13, no. 16, pp. 508–509, 1964. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - F. Englert and R. Brout, “Broken symmetry and the mass of gauge vector mesons,”
*Physical Review Letters*, vol. 13, no. 9, pp. 321–323, 1964. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - G. S. Guralnik, C. R. Hagen, and T. W. B. Kibble, “Global conservation laws and massless particles,”
*Physical Review Letters*, vol. 13, no. 20, pp. 585–587, 1964. View at Publisher · View at Google Scholar · View at Scopus - J. F. Gunion, H. E. Haber, G. L. Kane, and S. Dawson, “The Higgs hunter's guide,”
*Frontiers of Physics*, vol. 80, pp. 1–448, 2000. View at Google Scholar - A. Djouadi, “The anatomy of electroweak symmetry breaking. I: the Higgs boson in the standard model,”
*Physics Reports*, vol. 457, no. 1–4, pp. 1–216, 2008. View at Publisher · View at Google Scholar · View at Scopus - G. Aad, T. Abajyan, B. Abbott, et al., “Observation of a new particle in the search for the standard model Higgs boson with the ATLAS detector at the LHC,”
*Physics Letters B*, vol. 716, no. 1, pp. 1–29, 2012. View at Publisher · View at Google Scholar - S. Chatrchyan, V. Khachatryan, A. M. Sirunyan, et al., “Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC,”
*Physics Letters B*, vol. 716, no. 1, pp. 30–61, 2012. View at Publisher · View at Google Scholar - S. Weinberg, “A model of leptons,”
*Physical Review Letters*, vol. 19, no. 21, pp. 1264–1266, 1967. View at Publisher · View at Google Scholar · View at Scopus - S. L. Glashow, “Partial-symmetries of weak interactions,”
*Nuclear Physics*, vol. 22, no. 4, pp. 579–588, 1961. View at Publisher · View at Google Scholar · View at Scopus - A. Salam, “Elementary particle theory,” in
*Proceedings of the 8th Nobel Symposium*, N. Svartholm, Ed., p. 367, Wiley, New York, NY, USA, 1968. View at Google Scholar - J. Ellis, M. K. Gaillard, and D. V. Nanopoulos, “A phenomenological profile of the Higgs boson,”
*Nuclear Physics B*, vol. 106, pp. 292–340, 1976. View at Google Scholar · View at Scopus - B. Ioffe and V. A. Khoze, “What can be expected from experiments on colliding ${e}^{+}{e}^{-}$ beams with energy approximately equal to 100-GeV?”
*Soviet Journal of Particles and Nuclei*, vol. 9, p. 50, 1978. View at Google Scholar - M. A. Shifman, A. Vainshtein, M. Voloshin, and V. I. Zakharov, “Low-energy theorems for Higgs boson couplings to photons,”
*Soviet Journal of Nuclear Physics*, vol. 30, pp. 711–716, 1979. View at Google Scholar - M. B. Gavela, G. Girardi, C. Malleville, and P. Sorba, “A non-linear ${R}_{\zeta}$ gauge condition for the electroweak $SU(2)\times U(1)$ model,”
*Nuclear Physics B*, vol. 193, no. 1, pp. 257–268, 1981. View at Publisher · View at Google Scholar · View at Scopus - D. Huang, Y. Tang, and Y.-L. Wu, “Note on Higgs Decay into Two Photons $H\to \gamma \gamma $,”
*Communications in Theoretical Physics*, vol. 57, pp. 427–434, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H. Shao, Y. Zhang, and K. Chao, “Reduction schemes in cutoff regularization and Higgs decay into two photons,”
*Journal of High Energy Physics*, vol. 2012, no. 1, article 053, 2012. View at Publisher · View at Google Scholar · View at Scopus - F. Bursa, A. Cherman, T. C. Hammant, R. R. Horgan, and M. Wingate, “Calculation of the one W loop $H\to \gamma \gamma $ decay amplitude with a lattice regulator,”
*Physical Review D*, vol. 85, Article ID 093009, 2012. View at Google Scholar - F. Piccinini, A. Pilloni, and A. Polosa, “$H\to \gamma \gamma $: a comment on the Indeterminacy of non-gauge-invariant integrals,”
*Chinese Physics C*, vol. 37, no. 4, Article ID 043102. View at Publisher · View at Google Scholar - W. J. Marciano, C. Zhang, and S. Willenbrock, “Higgs decay to two photons,”
*Physical Review D*, vol. 85, no. 1, Article ID 013002, 10 pages, 2012. View at Publisher · View at Google Scholar - J. M. Cornwall, D. N. Levin, and G. Tiktopoulos, “Derivation of gauge invariance from high-energy unitarity bounds on the S matrix,”
*Physical Review D*, vol. 10, no. 4, pp. 1145–1167, 1974. View at Publisher · View at Google Scholar · View at Scopus - C. E. Vayonakis, “Born helicity amplitudes and cross-sections in non-Abelian gauge theories,”
*Lettere Al Nuovo Cimento Series 2*, vol. 17, no. 11, pp. 383–387, 1976. View at Publisher · View at Google Scholar · View at Scopus - B. W. Lee, C. Quigg, and H. B. Thacker, “Weak interactions at very high energies: the role of the Higgs-boson mass,”
*Physical Review D*, vol. 16, no. 5, pp. 1519–1531, 1977. View at Publisher · View at Google Scholar · View at Scopus - S. Weinberg, “General theory of broken local symmetries,”
*Physical Review D*, vol. 7, no. 4, pp. 1068–1082, 1973. View at Publisher · View at Google Scholar · View at Scopus - S. L. Adler, “Axial-vector vertex in spinor electrodynamics,”
*Physical Review*, vol. 177, no. 5, pp. 2426–2438, 1969. View at Publisher · View at Google Scholar · View at Scopus - J. Bell and R. Jackiw, “A PCAC puzzle: ${\pi}^{0}\to \gamma \gamma $ in the sigma model,”
*Il Nuovo Cimento A*, vol. 60, pp. 47–61, 1969. View at Publisher · View at Google Scholar - L. Rosenberg, “Electromagnetic interactions of neutrinos,”
*Physical Review*, vol. 129, no. 6, pp. 2786–2788, 1963. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - A. Dedes and K. Suxho, “Heavy fermion non-decoupling effects in triple gauge boson vertices,”
*Physical Review D*, vol. 85, no. 9, Article ID 095024, 31 pages, 2012. View at Publisher · View at Google Scholar - R. Jackiw, “When radiative corrections are finite, but undetermined,”
*International Journal of Modern Physics B*, vol. 14, no. 19-20, pp. 2011–2021, 2000. View at Publisher · View at Google Scholar · View at Scopus - G. 't Hooft and M. Veltman, “Regularization and renormalization of gauge fields,”
*Nuclear Physics B*, vol. 44, no. 1, pp. 189–213, 1972. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - J. C. Collins,
*Renormalization*, Cambridge Monographs on Mathematical Physics, Cambridge University Press, Cambridge, UK, 1984. View at Publisher · View at Google Scholar · View at MathSciNet - R. Pittau, “A four-dimensional approach to quantum field theories,”
*Journal of High Energy Physics*, vol. 2012, p. 151, 2012. View at Publisher · View at Google Scholar - R. Gastmans, S. L. Wu, and T. T. Wu, “Higgs decay $H\to \gamma \gamma $ through a W loop: difficulty with dimensional regularization,” http://arxiv.org/abs/1108.5322.
- R. Gastmans, S. L. Wu, and T. T. Wu, “Higgs decay into two photons, revisited,” http://arxiv.org/abs/1108.5872.
- J. M. Jauch and F. Rohrlich,
*The Theory of Photons and Electrons*, Springer, New York, NY, USA, 1976. View at MathSciNet - R. E. Pugh, “Origin shifts in divergent Feynman integrals,”
*Canadian Journal of Physics*, vol. 47, pp. 1263–1269, 1969. View at Publisher · View at Google Scholar · View at MathSciNet - V. Elias, G. McKeon, and R. B. Mann, “Shifts of integration variable within four- and n-dimensional Feynman integrals,”
*Physical Review D*, vol. 28, no. 8, pp. 1978–1992, 1983. View at Publisher · View at Google Scholar · View at Scopus - M. E. Peskin and D. V. Schroeder,
*An Introduction to Quantum Field Theory*, Addison-Wesley, Reading, Mass, USA, 1995. View at MathSciNet - J. Bagger and C. Schmidt, “Equivalence theorem redux,”
*Physical Review D*, vol. 41, no. 1, pp. 264–270, 1990. View at Publisher · View at Google Scholar · View at Scopus - M. Shifman, A. Vainshtein, M. B. Voloshin, and V. Zakharov, “Higgs boson decay into two photons through the W-boson loop: No decoupling in the ${m}_{W}\to 0$ limit,”
*Physical Review D*, vol. 85, no. 1, Article ID 013015, 2012. View at Publisher · View at Google Scholar · View at Scopus - F. Jegerlehner, “Comment on $H\to \gamma \gamma $ and the role of the decoupling theorem and the equivalence theorem,” http://arxiv.org/abs/1110.0869.
- T. Appelquist and J. Carazzone, “Infrared singularities and massive fields,”
*Physical Review D*, vol. 11, no. 10, pp. 2856–2861, 1975. View at Publisher · View at Google Scholar · View at Scopus - J. Abdallah, P. Abreu, W. Adam, et al., “Measurements of CP-conserving trilinear gauge boson couplings WWV ($\text{V}\equiv \gamma ,\text{Z}$) in ${e}^{+}{e}^{-}$ collisions at LEP2,”
*The European Physical Journal C*, vol. 66, pp. 35–56, 2010. View at Publisher · View at Google Scholar - J. M. Cornwall, D. N. Levin, and G. Tiktopoulos, “Uniqueness of spontaneously broken gauge theories,”
*Physical Review Letters*, vol. 30, no. 25, pp. 1268–1270, 1973. View at Publisher · View at Google Scholar · View at Scopus - M. S. Chanowitz and M. K. Gaillard, “The TeV physics of strongly interacting W's and Z's,”
*Nuclear Physics B*, vol. 261, pp. 379–431, 1985. View at Publisher · View at Google Scholar · View at Scopus - G. J. Gounaris, R. Kögerler, and H. Neufeld, “Relationship between longitudinally polarized vector bosons and their unphysical scalar partners,”
*Physical Review D*, vol. 34, no. 10, pp. 3257–3259, 1986. View at Publisher · View at Google Scholar · View at Scopus - S. Dawson and S. Willenbrock, “Radiative corrections to longitudinal-vector-boson scattering,”
*Physical Review D*, vol. 40, no. 9, pp. 2880–2887, 1989. View at Publisher · View at Google Scholar · View at Scopus - A. Vainshtein, V. I. Zakharov, and M. A. Shifman, “Higgs particles,”
*Soviet Physics Uspekhi*, vol. 23, pp. 429–449, 1980. View at Publisher · View at Google Scholar - B. A. Kniehl and M. Spira, “Low-energy theorems in Higgs physics,”
*Zeitschrift für Physik C*, vol. 69, no. 1, pp. 77–87, 1995. View at Publisher · View at Google Scholar · View at Scopus - A. Pilaftsis, “Higgs-boson low-energy theorem and compatible gauge-fixing conditions,”
*Physics Letters B*, vol. 422, no. 1–4, pp. 201–211, 1998. View at Publisher · View at Google Scholar · View at Scopus - C. G. Callan Jr., S. Coleman, and R. Jackiw, “A new improved energy-momentum tensor,”
*Annals of Physics*, vol. 59, no. 1, pp. 42–73, 1970. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - J. Horejsi and M. Stohr, “Higgs decay into two photons, dispersion relations and trace anomaly,”
*Physics Letters B*, vol. 379, pp. 159–162, 1996. View at Publisher · View at Google Scholar - S. L. Adler, J. C. Collins, and A. Duncan, “Energy-momentum-tensor trace anomaly in spin-1/2 quantum electrodynamics,”
*Physical Review D*, vol. 15, no. 6, pp. 1712–1721, 1977. View at Publisher · View at Google Scholar · View at Scopus - M. Klute, R. Lafaye, T. Plehn, M. Rauch, and D. Zerwas, “Measuring higgs couplings from LHC data,”
*Physical Review Letters*, vol. 109, no. 10, Article ID 101801, 10 pages. View at Publisher · View at Google Scholar - P. P. Giardino, K. Kannike, M. Raidal, and A. Strumia, “Is the resonance at 125 GeV the Higgs boson?”
*Physics Letters B*, vol. 718, no. 2, pp. 469–474, 2012. View at Publisher · View at Google Scholar