Helicity amplitudes without gauge cancellation for electroweak processes

In the 5-component representation of weak bosons, the first four components make a Lorentz four vector, representing the transverse and longitudinal polarizations excluding the scalar component of the weak bosons, whereas its fifth component corresponds to the Goldstone boson. We obtain the 5×5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$5\times 5$$\end{document} component propagators of off-shell weak bosons, proposed previously and named after the Goldstone boson equivalence theorem, by starting from the unitary-gauge representation of the tree-level scattering amplitudes, and by applying the BRST (Becchi–Rouet–Stora–Tyutin) identities to the two sub-amplitudes connected by each off-shell weak-boson line. By replacing all weak boson vertices with those among the off-shell 5-component wavefunctions, we arrive at the expression of the electroweak scattering amplitudes, where the magnitude of each Feynman amplitude has the correct on-shell limits for all internal propagators, and hence with no artificial gauge cancellation among diagrams. Although our derivation is limited to the tree-level only, it allows us to study the properties of each Feynman amplitude separately, and then learn how they interfere in the full amplitudes. We implement the 5-component weak boson propagators and their vertices in the numerical helicity amplitude calculation code HELAS (Helicity Amplitude Subroutines), so that an automatic amplitude generation program such as MadGraph can generate the scattering amplitudes without gauge cancellation. We present results for several high-energy scattering processes where subtle gauge-theory cancellation among diagrams takes place in all the other known approaches.


Introduction
It has been customary that scattering amplitudes in QED, QCD, and those in the standard model (SM) are expressed in terms of the Feynman amplitudes, where the gaugeboson propagators are expressed in a given gauge, such as in the Feynman gauge, Landau gauge, light-cone gauge, axial gauge, unitary gauge or renormalizable covariant R ξ gauge, etc.It has been shown that scattering amplitudes, or S-matrix elements of external on-shell particles, are independent of the gauge which we need to fix to quantize the gauge bosons.The gauge invariance of the amplitudes has thus been used to test our analytic or numerical computations.
Recently, a new type of gauge-boson propagators is proposed for massless gauge bosons [1], such as photons and gluons: a chenjm@jnu.edu.cnb kaoru.hagiwara@kek.jpc junichi.kanzaki@ipmu.jpd mawatari@iwate-u.ac.jp which has the form of the light-cone gauge propagator, gµν = g µν − n µ q ν + q µ n ν n • q , where the light-like vector n µ is chosen for each propagator as depending on the three momentum of the off-shell gauge bosons.Here, sgn(q 0 ) = 1 when q 0 ≥ 0, while sgn(q 0 ) = −1 otherwise.It has been shown in Ref. [1] that tree-level scattering amplitudes calculated with the above photon and gluon propagators agree exactly with the known gauge invariant helicity amplitudes.The agreement has been shown to follow from the BRST (Becchi-Rouet-Stora-Tyutin) identities [2] ⟨phys|{Q BRST , c}|phys⟩ = ⟨phys|∂ µ A µ |phys⟩ = 0, (4) where the BRST transformation of the anti-ghost field c gives the covariant gauge fixing term both in QED and QCD.The physical states are defined to be BRST invariant, or When an arbitrary scattering amplitude is connected by a photon/gluon propagator as with P R ξ µν = −g µν + (1 − ξ) q µ q ν q 2 (7) in the covariant (R ξ ) gauge, the theorem can be applied for the two sub-amplitudes T µ 1 and T ν 2 , giving q µ T µ k = 0 (8) for both k = 1 and 2. Because of the BRST identities for the two off-shell currents in Eq. ( 8), we can drop all terms proportional to q µ and q ν in the polarization tensor as Pµν = λ=±1 ϵ µ (q, λ) * ϵ ν (q, λ) + sgn(q 2 )ε µ (q, 0)ε ν (q, 0) where εµ (q, 0) = ϵ µ (q, 0) − with Q = |q 2 |, when q 0 > 0. The essence of the propagator in Eqs. ( 1) and ( 2) is hence that the component of the longitudinally polarized mode of the virtual photon/gluon which grows with q µ at high energies is not allowed to propagate.Indeed, εµ (q, 0) vanishes at high energies.Our notation for the polarization vectors is summarized in Appendix C. The striking property of the QED/QCD scattering amplitudes calculated with the propagators of the form (1) is that there is no subtle cancellation among Feynman amplitudes and that the individual Feynman amplitude shows all the on-shell singularities at all energies in an arbitrary Lorentz frame.Accordingly, the absolute value squared of the Feynman amplitudes calculated with the propagator (1) reproduces the well-known DGLAP splitting functions [3] in the collinear limit.Because of this property, the propagator of Eqs. ( 1) and (2) has been named 'parton-shower (PS) gauge' in Ref. [1].
In this paper, we extend the study of Ref. [1] to massive vector bosons in the SM of particle physics, W and Z bosons.The BRST identities for the physical states now read ⟨phys|{Q BRST , c}|phys⟩ = ⟨phys|(∂ µ V µ − ξm V π V )|phys⟩ = 0 (11) in the covariant R ξ gauge [4], where π V corresponds to the Goldstone boson associated with the weak bosons V = W ± and Z.We can apply the BRST identities (11) to a pair of an arbitrary tree-level amplitudes which are connected by a weak boson propagator, e.g. as with in the unitary gauge.We can express the unitary-gauge polarization tensor (13) as where gµν denotes the polarization tensor (2) or (9) of QED and QCD [1].As we show below in Sect.2, the BRST identities (11) for the sub-amplitudes where T π V k (k = 1, 2) are the Goldstone-boson amplitudes.The weak-boson propagators are then expressed as a summation over the contribution of the QED/QCD-like term, the mixing terms between the reduced longitudinal polarization vector (10) and the Goldstone boson, and the Goldstone-boson exchange term.All the unitary-gauge propagators are hence expressed in the 5 × 5 matrix form in the 5-component representation of the weak bosons.
The 5-component representation of the weak boson wavefunctions was introduced by Chanowitz and Gaillard [5] in their derivation of the Goldstone boson equivalence theorem.It was adopted by the authors of Refs.[6][7][8][9] in obtaining the equivalent distributions of the weak bosons in massless quarks and leptons.The 5 × 5 matrix form of the weak-boson propagators, Eq. ( 14), was introduced in Ref. [9] in their derivation of the electroweak (EW) splitting functions, and derived rigorously in Ref. [10] as a consistent covariant quantization of the EW theory.Although our derivation is limited to the tree-level amplitudes only, as explained in the following section, we keep it in the manuscript because it was inspired by the success of obtaining 'PS-gauge' amplitudes in unbroken gauge theories [1] and also because its simplicity has allowed us to develop an automatic amplitude generation code in the SM.
In this paper, we present explicit expressions for all the weak boson propagators and their vertices in the SM, and prepare numerical codes to calculate arbitrary treelevel amplitudes by updating HELAS (HELicity Amplitude Subroutines) [11,12] for automatic amplitude generation programs such as MadGraph [13][14][15][16].
The rest of the paper is organized as follows.In Sect.2, we derive the 5-component weak boson wavefunctions and their 5 × 5 propagators, starting from the unitary-gauge amplitudes, and discuss the relationships with the amplitudes in the R ξ gauge.We also explain how to generate an arbitrary tree-level helicity amplitude with the 5component representation of weak bosons automatically by using MadGraph5 aMC@NLO [16].In Sect.3, we present numerical results for weak boson scattering processes V V → V V for V = W ± or Z, W-boson pair production processes l − l + → W − W + → f 1 f2 f 3 f4 , and l − l + → l ′l′ V V (l = e or µ, l ′ = l or ν l ).Section 4 summarizes our findings.
We give three appendices.Appendix A presents the new HELAS codes which are necessary to obtain helicity amplitudes in the 5-component weak-boson formalism.Appendix B gives all the Goldstone-boson couplings in the SM, while Appendix C summarizes our notation for polarization vectors.

Derivation
As stated above, the derivation of the 5 × 5 component weak-boson propagator given in this section is valid only for the tree-level scattering amplitudes.For a more general derivation, we refer the readers to Ref. [10].

5-component polarization vector of weak bosons
Let us denote a helicity amplitude with an external weak boson V (V = W ± or Z) as where q and λ are the four momentum and the helicity of the vector boson.When the vector boson is in the final state, the polarization vector should be replaced by its complex conjugate.The three polarization vectors take a simple form in the frame where the three momentum − → q is along the z-axis, Here, y is the rapidity of the vector boson along its three momentum direction.
Because the polarization vector for longitudinally polarized (λ = 0) vector boson grows with E/m V or as e y at high energies, the individual Feynman amplitude for longitudinally polarized vector bosons often grows with energy, violating perturbative unitarity.In theories where the vector boson masses are obtained by spontaneous breaking of the gauge symmetry, all these unitarity violating terms cancel after summing over all the interfering Feynman amplitudes.The subtle gauge-theory cancellation among interfering amplitudes can lead to the loss of terms of order (E/m V ) n in numerical calculation of amplitudes for processes with n external weak bosons.The factor can be orders of 10 8 for weak boson scattering processes (n = 4) at E = 10 TeV.
The BRST invariance of quantized EW theory allows us to obtain the helicity amplitudes without subtle cancellation.We first split the λ = 0 polarization vector as where in the frame where the three momentum of the vector boson is taken along the z-axis as in Eq. ( 18).The BRST identities (11) for the helicity amplitudes M λ (16) give where T π V is the amplitude of the Goldstone boson π V which gives the mass to the gauge boson V when the gauge invariance is broken spontaneously. 1 The sign of the lefthand-side of Eq. ( 21) should be reversed when the vector boson is in the final state.
An arbitrary amplitude of the longitudinally polarized massive gauge bosons (16) can hence be expressed as where the index M runs from 0 to 4, M = {µ, 4}, with and when the weak boson V is in the initial state.When the weak boson V is in the final state, all the polarization 'vectors', εM V (q, λ), should be replaced by their complex conjugates.Since all the components of εµ V (q, 0) in Eq. ( 20) vanish at high energies as m V /(E + | − → q |) = e −y , and since the Goldstone boson amplitudes T π V do not grow with energy, we can tell that the λ = 0 vector boson scattering amplitudes should reduce to the associate Goldstone boson amplitudes at high energies, leading to the Goldstone boson equivalence theorem [5,18].As we show in Subsection 3.1, we study all 2 → 2 weak-boson scattering amplitudes, and verify numerically that the BRST identity (21) holds even when one or more external weak-boson polarization vectors are replaced by the 5-component wavefunctions (23).
If we keep the finite contribution of εµ V (q, 0), which can be significant at low energies, we have an exact expression of the amplitudes without gauge-theory cancellation.This has been shown first in the axial gauge in obtaining the equivalent real weak-boson distributions [6,7], and the weak-boson splitting amplitudes [9], and also in the covariant gauge in Refs.[8,10]. 1 It should be noted that the BRST identities, Eqs.(15a) and (15b), are valid in an arbitrary R ξ gauge, including the unitary gauge (ξ → ∞) in the tree level.This is because the sub-amplitudes T µ V and T π V in Eq. ( 21) are truncated, and all the Goldstone-boson couplings are ξ independent.The BRST identities, Eq. ( 21), in the unitary gauge has been used in checking the correctness of numerically evaluated scattering amplitudes [12,17].
At this stage, we have obtained the expression of the weak-boson scattering amplitudes which have no external wavefunction components that grow with energy.However, because we adopt the unitary gauge for the internal propagators, individual Feynman amplitude can still have components which grow with energy when the energy of the exchanged off-shell weak boson is large.This problem can be solved by replacing all the unitary-gauge propagators by the 5 × 5 component propagators of the 5-component weak bosons [9,10].

Propagators of the 5-component weak bosons
The BRST identities, Eqs.(15a) and (15b), allow us to replace one unitary-gauge propagator in an arbitrary treelevel scattering amplitude by the 5×5 component propagator of the 5-component representation of the weak bosons: with for k = 1 and 2, and or in the matrix representation.
In the first step, the subamplitudes T µ 1 and T ν 2 are both off-shell weak-boson currents connected by all Feynman diagrams to on-shell states in the unitary gauge.What we verified by using our numerical code is that the identities (26) hold even if the off-shell weak-boson line is connected by the 5 × 5 propagator (28).Because of this, all the unitary-gauge propagators can be replaced by the 5×5 propagators.
Let us give a few key properties of the 4×4 components of the polarization tensor PV µν , which is exactly the same as the special light-cone gauge propagator introduced in Ref. [1] for the photons and gluons, Eqs. ( 2) and ( 9).The light-cone vector n µ (3) has its three-vector component along − − → q and the time component is fixed to give for an arbitrary four momentum q µ .The reduced polarization vector for the λ = 0 state (10) can be expressed as εµ (q, 0) = sgn(q 2 ) e −y (−1, 0, 0, 1), (30) when the three momentum − → q is along the positive z-axis with q z = Q sinh y for the time-like, and q z = Q cosh y for the space-like four momentum.The reduced polarization tensor becomes both for the time-like and space-like momentum.Therefore, the longitudinal polarization mode is suppressed by the e −2y factor as compared to the transverse polarization modes at high energies.
For massive vector bosons, the remaining part of the longitudinally polarized weak bosons and also the scalar (j = 0) component can propagate.For instance, the polarization transfer tensor in the unitary gauge can be expressed as The first term gives the propagation of the spin-1 (j = 1) components just as in QED and QCD [1], while the second term gives transitions between the reduced longitudinal polarization state and the spin-0 (j = 0) component, and the last part gives the propagation of the spin-0 component.This is manifest in the equivalent expression given in Eq. (14).Applying the BRST identities (15), which are valid at an arbitrary q µ , we arrive at the 5-component weak boson propagator in the matrix representation.Note that the expression in Eq. (33) agrees with Eq. ( 28) when q 2 > 0, because Eq. ( 28) is obtained from Eq. ( 15) when the weak boson is produced by T µ 1 with the time-like momentum and then decays into T ν 2 .The form of our 5 × 5 polarization transfer matrix (33) is valid for an arbitrary four momentum q µ , including the space-like momentum with q 2 < 0. This can be confirmed by replacing q µ by −q µ in the propagator; the diagonal terms in Eq. (33) are invariant under the replacement.In the off-diagonal terms, n • q, given in Eq. (29), is invariant, while n µ changes sign.At the same time the Goldstone-boson flow is reversed, resulting in the exchange between i and −i in the BRST identities (15).We hence recover the same form (33).

Implementation in MadGraph5 aMC@NLO
The 5-component representation of the weak bosons and their 5×5 polarization transfer matrix propagators can be a powerful tool in the study of the SM physics if they can be implemented into an automatic amplitude generation program, such as MadGraph5 aMC@NLO (MG5) [16].In this paper, we develop a working model which can be applied to calculate tree-level helicity amplitudes of an arbitrary processes in the SM.We first generate the amplitudes in the unitary gauge of the SM weak bosons by using MG5 [16].MG5 produces numerical codes to compute helicity amplitudes corresponding to each Feynman diagram.The numerical codes calculate the helicity amplitudes by successively calling HELAS subroutines [11], which calculate off-shell fermion and vector boson wavefunctions starting from the external onshell particle wavefunctions.The amplitude is obtained as a complex number when all the off-shell wavefunctions meet at one vertex. 2 We modify the program as follows: 1.When an external particle is a weak boson (W ± or Z), we replace the external wavefunction with our 5component wavefunction in Eq. ( 23). 2. All the vertices with external weak bosons should be extended to contain Goldstone boson contributions.3.All the unitary-gauge propagator of weak bosons should be replaced by the 5×5 matrix propagator in Eq. (33).4. All photon and gluon propagators should be given in the PS gauge in Eqs. ( 1) and (2).
The new HELAS subroutines for the 5-component representation of a massive vector boson (step 1), those for all the weak boson vertices in the SM (step 2), and those for the 5-component off-shell vector-boson wavefunctions (step 3) have been created, and they are listed in Appendix A. The off-shell photon and gluon propagators in the PS gauge (step 4) have been given in Ref. [1].
We note here that all the off-shell HELAS subroutines for a n-point vertex have n − 1 inputs and one off-shell output.When all the n − 1 inputs are on-shell states, the subroutine should hence satisfy the BRST identities (15).All our new HELAS subroutines have been tested explicitly by using the BRST, with the double complex accuracy of 12 to 13 digits.The BRST identities turn out to be valid when one or all n − 1 inputs are replaced by off-shell HELAS subroutines which are obtained from the external 2 While MadGraph/MadEvent v4 [14] employs the original HELAS library [11], MadGraph5( aMC@NLO) [15]( [16]) adopts ALOHA [19], which automatically generates the HELAS library by using the UFO format [20].
on-shell states.In this way, we have checked step-by-step the BRST invariance of all off-shell vertices which appear in the tree-level SM amplitudes.
It turns out that the above four steps are not sufficient to generate all the SM amplitudes.This is because there are four 4-point vertices in the SM which couple only to the Goldstone boson(s) (the fifth component of our weakboson wavefunctions) but not to the corresponding weak boson(s) in the unitary gauge.They are listed in Fig. 1: (a) The quartic coupling of the Goldstone boson π Z = π 0 .(b) The quartic gauge coupling of the Goldstone boson π 0 , which couples to Z-boson pair.(c) The quartic coupling among the photon (γ) or Z, H, W ± and the Goldstone boson of W ∓ (π Because the only way that we could identify these four missing vertices in the unitary gauge has been to check all Goldstone vertices in the SM, we present our parametrization of the SM Higgs sector in Appendix B. The quartic coupling of π 0 (Fig. 1(a)) appears in the Higgs potential, as a consequence of its custodial SU(2) invariance, under which the three Goldstone bosons π k with k = 1, 2, 3 transform as a triplet.All the other quartic vertices, Fig. 1(b,c), appear in the gauge coupling of the Higgs doublet, as so-called seagull terms.
Those four vertices do not appear among the weak boson vertices in the unitary gauge, and hence we should add them as additional vertices in the SM.This can easily be achieved in the SM UFO model by introducing the following three additional interactions, ZZZZ, W W HZ and W W Hγ. MG5 then generates helicity amplitudes with the HELAS subroutines with the above vertices.It is then straightforward to prepare the HELAS subroutines by using the Goldstone boson couplings of the SM.These additional codes are also given in Appendix A.
It should be noted that these quartic vertices are needed in our representation of the scattering amplitudes with the 5-component weak bosons, because the BRST identities, Eqs.(15a) and (15b), relate the total sum of all the contributing amplitudes rather than individual Feynman amplitudes.For example, the ZZZZ sub-amplitudes with four unitary gauge Z-boson propagators generate the Goldstone vertices like Fig. 1(a) and (b).Therefore, the scalar (j = 0) polarization component of the physical (unitary-gauge) weak boson interacts with the Goldstoneboson vertices which are gauge (ξ) independent as shown explicitly in the Higgs Lagrangian (given in Appendix B for the minimal SM).

Comments on amplitudes in the R ξ gauge
The tree-level amplitudes expressed in the unitary gauge have a special property that when we cut one gauge-boson propagator of a given four momentum in all the Feynman amplitudes, the full helicity amplitude is split into two pieces, each satisfying the BRST invariance.This is because when the cut vector-boson four momentum is set on-shell, each sub-amplitude can be regarded as an independent scattering amplitude.
In general renormalizable covariant gauge, the R ξ gauge with finite ξ [4], there can appear both a vector-boson and its associate Goldstone-boson propagator which share the same four momentum.In our 5-component representation of the weak bosons, its fifth component is interacts with the Goldstone boson couplings.It might therefore seem that our 5 × 5 matrix propagator of the weak bosons contains the associated Goldstone boson propagators in the Feynman gauge (ξ = 1).However, this is not the case, as explained below.
This can be seen clearly by denoting the massive weakboson propagator in the R ξ gauge in the following form: where D U V µν (q) is the unitary-gauge propagator.Upon making use of the BRST identities, the second term in Eq. ( 35) gives the negative of the Goldstone-boson contribution in the R ξ gauge with the propagator whenever the same four-momentum is shared by the gauge boson and its associated Goldstone boson.We can hence conclude that their sum should agree with the contribution of the unitary-gauge propagator, and hence that of our 5 × 5 propagator.This statement is valid for an arbitrary R ξ gauge, including the Feynman gauge (ξ = 1).Because of this cancellation, Goldstone-boson propagators cannot appear inside the tree-level scattering amplitudes.Therefore even if we start from the scattering amplitude expression in the R ξ gauge, the vector-boson and its associated Goldstone-boson contributions sum up to be replaced by the unitary-gauge propagator contribution.We find that in the R ξ gauge amplitudes, we always have both a gauge-boson and its associated Goldstone-boson lines which share the same four-momentum, and this cancellation always takes place. 3lthough it may seem that the Goldstone boson propagates inside our amplitudes, none of them are the Goldstone bosons of the R ξ gauge.We can verify this by keeping the ξ factor in the mass term of the Goldstone-boson and the weak-boson propagators of the R ξ gauge, Eqs. ( 36) and (35), respectively.The two contributions cancel out exactly in the scattering amplitudes, and only the unitarygauge part of the R ξ -gauge propagators, the first term in Eq. (35), survives.It is the scalar component of this remaining unitary-gauge polarization tensor, Eq. ( 14), which interacts with the Goldstone-boson vertices by the BRST identities, Eqs.(15a) and (15b).The propagating degrees of freedom are just that of the unitary-gauge weak bosons, Eq. ( 12).
We can therefore arrive at the same representation of the amplitudes by starting from the R ξ gauge with finite ξ, and replacing all the pair of the gauge-boson and its associated Goldstone-boson propagators by our 5×5 component propagators.The derivation of our amplitudes given in this section tells that the resulting amplitudes should agree with the one which we obtain from the unitary-gauge expression, with the additional four types of vertices, given in Fig. 1.In the SM, we identify the four vertices of Fig. 1 which do not have weak-boson counterparts in the unitary gauge.In models beyond the SM, we should start from a new set of Feynman rules which treat the Goldstone bosons as the 5th component of the weak bosons.In the SM effective field theory (SMEFT), the Goldstone bosons can have additional interactions in higher dimensional operators, whereas in general models with spontaneously broken gauge symmetry, the Goldstone bosons of the SM and additional weak bosons can have components in additional Higgs multiplets.The Feynman diagrams with the Goldstone-boson couplings are then generated automatically as the vertices among the fifth component of the weak bosons.

Sample results
In order to demonstrate the power of the 5-component description of weak bosons in numerical calculation of the amplitudes, we study several high-energy EW scattering processes.In the following, we denote "New HELAS" as the 5-component weak-boson description introduced in this paper, and "HELAS" as the original HELAS [11], in which the unitary-gauge and Feynman-gauge propagators are employed for massive and massless gauge bosons, respectively, and compare the two results.

Weak boson scattering processes
First, we study weak boson scattering processes: < l a t e x i t s h a 1 _ b a s e 6 4 = " p h s 4 Figure 5 shows the total cross section of W − W + → ZZ as a function of the total energy √ ŝ from 200 GeV up to 100 TeV.Solid lines show the cross sections computed by the new HELAS subroutines, while dashed ones are by the original HELAS.Lines with black squares denote helicity-summed cross sections, while those with red circles represent cross sections among longitudinally polarized weak bosons.They are both calculated with the W-boson propagator with its physical width as , since it was a default setting of MadGraph until recently. 4 It is well known that the width of the weak bosons spoils gauge cancellation.Although the uncanceled terms are of higher order in the EW perturbation theory, when the amplitudes are calculated in the unitary gauge, violation of subtle cancellation among individual Feynman amplitudes spoils the accuracy of the cross section at high energies.On the other hand, the Feynman amplitudes calculated with new HELAS do not have subtle gauge cancellation at high energies, and the weak boson width effect, whether it is needed or not, remains small at all energies.
We confirmed that, if we set Γ V = 0 in the t/u-channel propagators of the weak bosons in the original HELAS computation, the total cross sections recover the perturbative unitarity and agree with the new HELAS results.In the following studies, we always set the width of weak bosons to zero in the t/u-channel propagators, while we should introduce the width term to unitarize the amplitudes near the weak-boson mass shell, as we study in Sect.3.2. 54 The weak boson widths in the space-like propagators are automatically set to zero in MG5 v2.8.0 or later [16]. 5Issues on gauge invariance of the amplitudes around the weak-boson mass-square pole have been studied in Ref. [10].In Fig. 6 we show the scattering angle distribution for TeV by new HELAS (left) and by original HELAS (right), where all the helicities for the external weak bosons are taken to be longitudinal.A solid line denotes the total differential distribution, while dashed lines with circles, triangles and squares show the distributions of the amplitude squared of the s-channel, t/u-channel and the contact diagrams, in Fig. 2(a), (b),  (c) and (d), respectively.
The physical distributions, calculated from the absolute value square of the total sum of all the Feynman amplitudes, are identical between the new and original HELAS calculations as they should be.However, the contribution from each amplitude is quite different between the two calculations.While there is no energy growth of each amplitude at all in the new HELAS, one can observe about O(10 15 ) cancellation among the amplitudes in the original HELAS, which is computationally expensive to obtain physical cross sections accurately at very high energies.
In the new HELAS amplitudes with the 5-component weak-boson description, no term which grows with energy (|q 0 | ≫ Q) appears, and the absolute value of each Feynman amplitude is expressed as a product of the Lorentz invariant propagator factor D V (q 2 ) = (q 2 −m 2 V ) −1 .As expected, the contributions from the t/u-channel amplitudes are dominant especially for the forward region, while the ones from the contact amplitude are constant and significant only in the central region.We find that the contribution from the s-channel Higgs boson exchange amplitude is completely negligible at this energy.
In contrast, the original HELAS amplitudes given in the right panel of Fig. 6 show that the square of t/u-channel exchange amplitudes (blue triangles) and the contact term (red squares) are both almost flat in cos θ, with 10 to 15 orders of magnitude larger than the physical distribution (black solid).Likewise, in Fig. 7 we show the angular distributions for W − W + → W − W + .In the new HELAS results given in the left panel, the t-channel γ and Z exchange amplitudes dominate the cross section at cos θ ∼ 1, while we can study the interference pattern among the different amplitudes quantitatively in the central (cos θ ∼ 0) and the backward (cos θ < 0) regions.On the other hand, each amplitude in the unitary gauge, i.e. in the original HELAS calculation, apparently does not give any useful information about physical properties of the scattering process.As text books describe, only after summing all the amplitudes, the energy-growing terms are cancelled by the gauge symmetry, and then the physically observable distributions given by black solid curves are the same between new HELAS (left panel) and the original HELAS in the unitary gauge (right panel).
Figure 8 shows the cos θ distribution for ZZ → ZZ.The physical distribution is flat in both new HELAS (left) and original HELAS (right).The origin of the flat distribution is the dominance of the contact interaction term (red squares) from the Goldstone boson coupling in Fig. 1(a the s-channel (green circles) and t/u-channel (blue triangles) Higgs exchange amplitudes, each giving three orders of magnitude larger absolute values.

W-boson pair production at lepton colliders
Next, we study W-boson pair production in lepton collisions, whose Feynman diagrams are depicted in Fig. 9.
Figure 10 shows the total cross section for e − e + → W − W + including the subsequent leptonic decays as a function of the collision energy.In the left (right) panel, the collision with the left-handed (right-handed) electron and the right-handed (left-handed) positron is given.
The total cross section, denoted by a black solid line, falls as 1/s and agree between the two calculations as expected.However, the contributions from each amplitude are completely different between new HELAS (solid curves) and the original unitary-gauge HELAS (dashed curves).In the 5-component description with new HELAS, we find that the t-channel neutrino-exchange amplitude is dominant for the left-handed electron beam at all the energies.For the right-handed electron shown in the right panel, on the other hand, we can identify the interference between the s-channel photon and Z-boson amplitudes.The sum of the square of the photon and Z-boson exchange diagrams, as shown by the green solid curves, is about ten times larger than the observable cross section, given by the black solid curve, at all energies.This has a clear physics interpretation.At energies high above the Z-boson mass, only the U(1) gauge boson B µ couples to the right-handed electron.B µ deos not couple to the W 's because the W hypercharge is zero.It is only the Goldstone boson component of the W , π ± , which couples to the B µ boson since they have hypercharge.For the individual γ and Z exchange amplitudes, both the SU(2) L triplet W ± and the doublet π ± contributes, while in the sum, for the e − R e + L annihilation, the W ± components cancel out at high energies.The produced π ± decays via longitudinal W mixing into a lepton pair.
On the other hand, the unitary-gauge amplitudes with original HELAS give contribution of each Feynman diagram which grows with energy, as shown by dashed curves in both left (e − L ) and right (e − R ) plots in Fig. 10.At √ s = 100 TeV, the cancellation in the amplitudes give 10 −9 and 10 −12 times the individual contribution for e − L and e − R , respectively.This 'fake' cancellation between the s-channel γ and Z exchange amplitudes (with γW W and ZW W couplings) and the t-channel ν exchange amplitudes gave us a wrong hope that we could have a better sensitivity to the triple gauge boson coupling at high energies.In order to overcome such unjustified expectation, gauge invariant operator formalism which may be a prototype of SMEFT was introduced in Ref. [23].With new HELAS amplitudes, we can tell that the magnitudes of the s-channel γ and Z exchange amplitudes (given by green solid curves in Fig. 10) are essentially the same for e − L (left panel) and e − R (right panel).Because of their interference with the dominant ν exchange amplitude (whose square is given by a blue solid curve with triangles), e − L annihilation process is more sensitive to the triple weak boson couplings.
In the new HELAS calculation, on the other hand, there is no subtle cancellation among the amplitudes.Moreover, each squared amplitude describes can be interpreted as parton distribution of the corresponding parton-shower history.For instance in the left-panel of Fig. 11(b), we observe that the physical cross section for e − L e + R → W − W + is dominated by the t-channel neutrino-exchange contribution in the forward direction, where the collinear limit is reached at cos θ = 1, whereas at cos θ < 0, where lower partial wave contributions dominate, and we observe destructive interference between the t-and s-channel amplitudes.We believe that it is the merit of the new HELAS amplitudes which allows us to make such qualitative statements on the magnitude and the interference patterns among interfering Feynman amplitudes.

Weak boson scattering at lepton colliders
Finally, we demonstrate weak boson scattering processes, discussed in Sect.3.1, in a more realistic setup, i.e. in lepton collisions.We study l − l + → l ′l′ V V (l = e or µ, l ′ = l or ν l , V = W ± , Z), specifically, for Z-boson and W-boson pair productions associated with the leptons.Here, in order to make physics discussions simpler, we consider e − µ + collisions in order to remove annihilation amplitudes.While these processes include the following weak boson scatterings: respectively, the majority of the Feynman diagrams consists of Z-or W-boson emissions from the external lepton currents.
Figure 13 shows total cross sections for each process in (40) as a function of the collision energy, where all the helicities in the initial and final states are summed.Similar to the previous sub-sections, we separate the contributions from each Feynman amplitude and show the sum of each amplitude squared by lines with symbols.Solid lines show the results computed by the new HELAS, while dashed ones are by the original HELAS.Since the processes (40c) and (40d) contain the t-channel photon exchange diagrams, we impose p T (l) > 10 GeV to avoid divergence in the forward region.
The common features for the four processes (40) in Fig. 13 are as follows.In the original HELAS computation, each amplitude squared grows as the collision energy is increased, which results in subtle cancellation among the amplitudes to obtain the physical cross section.In the new HELAS computation, on the other hand, the energy dependence of each amplitude squared behaves as the physical cross section, and subtle cancellations do not exist. 6 We note, however, that even in new HELAS amplitudes the cancellation can be seen most clearly in Fig. 13(c) for e − µ + → e − µ + ZZ (40c).Hints of similar cancellation can also be seen in other processes (40a,b,d), respectively shown in Fig. 13(a), (b) and (d), at very high energies.We investigate the origin of the cancellations among new HELAS amplitudes below.
In Fig. 14 we present p T distributions of the Z boson with smaller p T for e − µ + → e − µ + ZZ at √ s = 20 TeV, where p T (l) > 10 GeV is imposed as in Fig. 13(c).The total distribution is shown by a thick solid line.A solid line with squares denotes the contribution from the ZZ → ZZ 6 In Ref. [24] the authors report suppression of cancellation among the amplitudes for the process γγ → W − W + in the axial gauge [25].
scattering, which is computed by new HELAS as corresponding to the sum of the four amplitude squares, i.e. the s, t, u, and contact amplitudes, shown in Fig. 4. A solid line with crosses gives the contributions from the sum of the Z emission amplitude squares: as in Fig. 13(c).We can observe cancellation among the Z emission amplitudes, whose square sum in Eq. ( 43) is a factor of 10 3 larger than the physical cross section at the smallest pT (Z) bin of p T (Z) < 0.5 TeV at √ s = 20 TeV.From our experiences in QED, we examine the possibility that the cancellation is a consequence of soft Zboson emission from the initial and final charged leptons.
As a test of this postulate, we separate the Z emission amplitudes into three sub-groups such as Z emissions from the electron current, from the muon current, and one Z from the electron current and the other from the muon current, because we know that soft-photon emission amplitudes in QED are factorized into the three groups.A blue dashed line with circles in Fig. 14 gives the sum of the contributions from the three sub-groups, where the individual Feynman amplitudes (M k ) are summed within each sub-group before taking the absolute value square as No more subtle cancellation remains after the above 'soft-Z' emission amplitudes are summed up before squaring.As a further confirmation of the validity of our interpretation, we confirm that the majority of the small-p T Z boson from individual diagrams is transversely polarized.
As a final remark of our sample studies, the new HELAS amplitudes with the 5-component wavefunction and the 5 × 5 component propagators of the weak bosons are free from unphysical 'gauge theory cancellation' among individual Feynman amplitudes, and hence we can study physical origin of destructive or constructive interference among contributing amplitudes.

Summary
In this paper we successfully implement the 5 × 5 representation of the weak boson propagators [9,10] in an arbitrary tree-level helicity amplitudes of the electroweak theory which is free from subtle gauge-theory cancellation among interfering amplitudes.
In this representation, the longitudinally polarized weak boson states are expressed as a quantum mechanical superposition of the vector boson with reduced longitudinal polarization vector and its associate Goldstone boson with the same four momenta.The reduced polarization vector for the longitudinal polarization state εµ (q, 0) is defined as the difference between the longitudinal and the scalar polarization vectors as in Eq. ( 10).The amplitudes of the scalar polarization component are expressed as those of the Goldstone boson amplitudes, by making use of the BRST identities, Eq. ( 21), which are satisfied for the physical states in the S-matrix elements.We show that the same BRST identities hold for the two subamplitudes which are connected by a gauge boson propagator of an arbitrary tree-level scattering amplitudes.By applying the identities successively to all the gauge boson propagators in a scattering amplitude, we arrive at the expression of the amplitude in which all the off-shell weak boson propagators have the 5×5 polarization transfer matrix of Refs.[9,10], between the 5-component off-shell wavefunctions.The propagator reduces to the 'partonshower gauge' form for QED and QCD [1], in which only the transverse and the reduced longitudinal polarization states are allowed to propagate.
In this representation of the amplitudes, all the components of individual Feynman amplitudes which grow with four momentum of on-shell and off-shell gauge bosons are subtracted systematically, while the total sum of all the contributing amplitudes, the S-matrix element, remains intact by the BRST invariance.As a consequence of the removal of all components of the sub-amplitudes which grow with the weak boson energies, all the sub-amplitudes, those amplitudes associated with an individual Feynman diagram, are expressed as a product of the invariant propagator factors, D V (q 2 ) = (q 2 − m 2 V ) −1 , with the spin polarization sum at each vertex giving the splitting amplitudes [9,10,26].This allows us to interpret the full scattering amplitude as a superposition of all contributing subamplitudes (Feynman diagrams), and we can study interference among them quantitatively.
We believe that the above properties of the amplitudes in the 5 × 5 representation of weak-boson propagators [9,10] are valuable in the study of scattering amplitudes at high energies, and prepared a set of numerical codes, new HELAS codes, which allow us to obtain the tree-level amplitudes of an arbitrary SM processes by using an automatic Feynman amplitude generation code like MadGraph aMC@NLO [16].Some of the representative codes are presented in Appendix A, and sample results are reported in Sect.3, for the weak boson scattering processes (W − W + → ZZ, W − W + → W − W + , ZZ → ZZ), and for the processes where the above 2 → 2 weak boson scattering appears in processes with massless leptons in the initial state.All our findings look encouraging, allowing us to interpret the property of each contributing Feynman amplitude and the interference among them.
The gauge boson propagators are available only after gauge fixing, and their explicit forms depend on the gauge-fixing condition.A particular form of the gauge boson propagators has hence been named after the chosen gauge-fixing condition, such as the Feynman gauge, Landau gauge, axial or light-cone gauge, covariant renormalizable R ξ gauge, etc.The particular form of the massless gauge-boson propagators with the reduced longitudinal polarization vector has been named as 'parton shower gauge' [1], because of the property that the absolute value square of each sub-amplitude give the lepton, quark, photon and gluon splitting functions [3,26].For the EW gauge bosons, the 5-component representation of the weak bosons and their 5×5 propagators are called 'equivalent gauge' [8,10] or as 'Goldstone equivalence gauge' [9], referring to the original 5-component description of the weak bosons adopted in deriving the Goldstone boson equivalence theorem [5].That this particular form of the propagator is useful in obtaining the splitting functions among weak bosons and the Higgs boson has been clearly observed in Refs.[9,10].Since weak bosons can be treated as partons at high energies, the expression like 'parton shower' may be a good common description of our propagators and amplitudes.However, the name 'parton shower gauge' was adopted by Nagy and Soper [27,28] for a specific lightcone gauge in their study of parton-shower properties with quantum interference.We may call the gauge-boson propagator representations of unbroken [1] and broken [9,10] gauge theories as 'Feynman diagram (FD) gauge', because of the common property that the individual Feynman amplitude can be expressed as a product of invariant propagator factors for the internal lines and the splitting amplitudes at each vertex.Zheng for sending us valuable comments on the use of preliminary versions of new HELAS codes.All the Feynman diagrams were drawn with TikZ-FeynHand [29].JMC is supported by Fundamental Research Funds for the Central Universities of China NO.11620330.The work was supported in part by World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan, and also by JSPS KAKENHI Grant No. 18K03648, 20H 05239, 21H01077 and 21K03583.

A HELAS subroutines with 5-component weak bosons
In this appendix we list and present new HELAS subroutines needed to compute amplitudes and off-shell currents with 5-component weak bosons.We refrain from presenting all of them, but we show the representative subroutines for each type of vertices, i.e. those to compute amplitudes.Only for the fermion-fermion-vector boson vertex, we also present the subroutines to compute the off-shell currents as examples.They are publicly available on the web site. 7e already modified the internal structure of HELAS wavefunctions in our previous paper [1] to make it easier to extend them.In the new HELAS, which includes the new 5-component weak boson calculation, the size of the complex array of the vector wavefunction becomes seven to have the fifth component of wavefunctions.
The coupling constants of the new HELAS subroutines become the array of double-precision complex numbers to include new 5-component weak-boson interactions.Tables in each subsection show the coupling constants of each interaction vertex.The content of the array of coupling constants should be permutated depending on the order of input/output particles because the vertices of the new HELAS library with 5-component weak bosons are not symmetric for some interaction vertices.Other Tables describe the order of permutations in subsections if necessary.
We note that the original HELAS library [12] was coded by Fortran77, while the new one was done by Fortran90.

A.1 Parameter file and utility subroutine
In the module, helas params, in helas modules.f90, a new value of helas mode=3 is introduced for the 5-component weak boson calculations.The subroutine, hlmode, can change this internal flag by calling it with a mode number, helas mode.The subroutine call, call hlmode(3), changes gauge definitions of gauge-boson propagators.
In the module, helas utils, some utility subroutines are also included: define gauge dir calculates n-vector (Eq.( 3)) from particle four-momentum.

A.2 Wavefunction subroutines
The function, vxxxxx (List 1), computes the vector particle wavefunction of the 5-component vector boson representation.
The coupling constants of the fermion-fermion-vector boson vertices are listed in Table 1.

A.4 VVV vertex
Subroutines in this section compute amplitudes (vvvxxx) and off-shell currents (jvvxxx) of the three-point vector boson vertex.The subroutine for the amplitudes is shown explicitly in List 5.
The coupling constants and their permutations according to the order of input particles for the three-point vector boson vertices are listed in Table 2 and Table 3.

A.9 Interface to MadGraph5
The new HELAS library contains the subroutines to interface the HELAS subroutines to the MadGraph5 [15] subroutines.Amplitude subroutines generated by MadGraph5 can calculate helicity amplitude with 5-components vector bosons by linking with the new HELAS library.

B Goldstone boson couplings in the SM
In this appendix, we present all the Goldstone-boson couplings of the SM explicitly, because not only the magnitude but also the relative signs of all the gauge-boson and the corresponding Goldstone-boson couplings should be kept exact in order to keep the BRST invariance of the amplitudes.The Goldstone bosons appear in the four sector of the SM Lagrangian, the Higgs potential V H , the Higgs gauge interactions K H , the gauge-fixing term L GF , and the Yukawa term L Y .We show all of them explicitly below.
The minimal SM has just one SU(2) L doublet Higgs field ϕ(x).The Lagrangian of the Higgs field is given by with We parametrize the Higgs doublet field with hyper-  Although the Goldstone-boson interactions in the Higgs potential (B.12) take a very simple form due to its custodial SU(2) symmetry, their gauge couplings look more complicated because the U(1) Y and hence U(1) EM gauge interactions violate the custodial symmetry explicitly.
In the renormalizable covariant (R ξ ) gauge [4], the gauge-fixing terms are which allow us to quantize the gauge fields, in a covariant manner, while removing the weak-boson-Goldstoneboson mixing terms in (B.15c) and (B.15d), and giving the Goldstone-boson masses ξ Z m 2 Z and ξ W m 2 W for π 0 and π ± , respectively.We may call the choice ξ A = ξ Z = ξ W = 1 as the Feynman-'t Hooft gauge.

Fig. 1 .
Fig. 1.Contact interactions which do not appear in the unitary-gauge weak-boson vertices.

whose
Feynman diagrams are shown in Figs. 2, 3 and 4, respectively.As discussed in Sect.2.3, the 4-point Z-boson diagram, Fig. 4(d), has been introduced to account for the Goldstone boson amplitudes, Fig. 1(a) and 1(b), via the fifth component of the Z boson wavefunction.

Fig. 5 .
Fig.5.Historical plot of the total cross sections for W − W + → ZZ as a function of the colliding W W energy, when W-boson propagators exchanged in the t/u-channels are given the physical width.Solid lines show the cross sections computed by the new HELAS subroutines, while dashed ones are by the original HELAS.Lines with squares (circles) denote helicity-summed (only longitudinally polarized) cross sections.

Fig. 6 .Fig. 7 .Fig. 8 .
Fig.6.Distributions of the scattering angle for W − W + → ZZ at √ ŝ = 10 TeV by using the new (original) HELAS in left (right) panel, where all the helicities for the external bosons are taken to be longitudinal.A solid line denotes the total distribution, while dashed lines with circles, triangles and squares show the distributions of the amplitude squared of the s-channel, t/uchannel and the contact diagrams, respectively.

Fig. 10 .Fig. 11 .Fig. 12 .
Fig. 10.Total cross section for e − e + → W − W + → µ − νµe + νe with the left-handed (right-handed) electron and the righthanded (left-handed) positron as a function of the collision energy in the left (right) panel.A black solid line shows the total cross section, while lines with symbols denote contributions from the amplitude squared of each diagram.Solid lines are the results computed by the new HELAS subroutines, while dashed ones are by the original HELAS.
11 and 12, we present differential distributions for e − L e + R → W − W + → µ − νµ e + ν e at √ s = 250 GeV in the left two panels (a) and at 10 TeV in the right two panels (b).Fig. 11 gives cos θ distribution of W − , whereas Fig. 12 gives p T (µ) distribution in the lab (colliding e + e − center-of-mass) frame.For the angular distributions in Fig. 11, in the unitary gauge by the original HELAS, the cancellation among the amplitudes is still mild at √ s = 250 GeV, especially at cos θ ≳ 0.5, but it becomes very large at √ s = 10 TeV,

Fig. 13 .
Fig.13.Total cross sections for l − l + → ν νZZ (a), l − l + → ν νW − W + (b), l − l + → l − l + ZZ (c), and l − l + → l − l + W − W + (d) as a function of the collision energy.A Black solid line shows the total cross section, while lines with symbols denote contributions from the amplitude squared of each Feynman diagram in different categories.Solid lines show the results computed by the new HELAS, while dashed ones are by the original HELAS.

2 Fig. 14 .
Fig. 14. pT distributions of the Z boson with smaller pT for e − µ + → e − µ + ZZ at √ s = 20 TeV.A thick solid line shows the total distribution, while other lines with symbols are explained in the text.

2 π
+ π − , (B.15j) where we drop repeated Lorentz indices.Let us briefly explain each term.The first line (B.15a)gives kinetic terms of the Higgs boson and the Goldstone bosons.The terms with (v + H) 2 in (B.15b) give W and Z boson masses and their Higgs-boson couplings, while the last term with (π 0 ) 2 contributes to W W ZZ and ZZZZ vertices via the Goldstone (fifth) component of the Z boson.The next terms in (B.15c) and (B.15d) contain the bilinear Goldstone-boson-gauge-boson mixing terms which are removed by the R ξ gauge fixing, as given below, and the associate ZHH, ZZH, W W H coupling terms through the Z and W Goldstone-boson components.All the remaining terms in (B.15e)-(B.15j)contribute to three-or four-point vertices among the Higgs boson and the weak bosons with one or two Goldstone-boson components.