1 Introduction

The high energy limit of vector boson scattering in the electroweak sector of the standard model has been investigated in detail in [1]. The results are obtained in the multi-Regge limit of the theory and parallel those in QCD within the BFKL approach [2,3,4,5,6,7]. In [1] the main result is the derivation of an integral equation suited for the description of vector particle scattering in the vacuum exchange channel. The emerging picture is that of interacting reggeized charged particles and non-reggeizing neutral vector bosons. The consistency of these effective high energy degrees of freedom is sustained by the isospin-1 exchange bootstrap equations derived in [1].

Within this framework, our contribution in the present work is to offer an effective action useful to handle these high energy fields and evaluate efficiently scattering amplitudes for the interaction among electroweak vector bosons. We show how to reproduce the reggeization scheme above mentioned together with the calculation of effective production vertices needed to construct the integral evolution equations presented in [1]. Our calculational framework is based on that proposed by Lipatov [8,9,10] and our own work in the QCD context [11,12,13,14,15,16,17,18,19,20].

One of the virtues of this effective action approach is to set up a powerful tool to evaluate higher order corrections in the electroweak high energy effective theory. Special attention is put to the role of the scalar Higgs boson in the unitarity properties of these amplitudes and its contribution to the Pomeranchuk singularity governing the growth with energy of the total cross section. Our results should have an interesting interplay with the electroweak parton distributions discussed in [21]. They can be of importance also in the evolution of parton distribution and fragmentation functions up to scales beyond the electroweak scale (where the \(SU(2) \times U(1)\) symmetry is unbroken) [22,23,24,25,26]. Lipatov’s high energy effective action has had an abundant number of applications in QCD, see, e.g. [27,28,29,30,31,32,33,34,35,36], in supersymmetric theories [37, 38] and gravity [39,40,41,42,43,44].

In the next section we introduce the effective action for a general non-abelian SU(N) theory where the associated gauge bosons are known to reggeize in the leading logarithmic approximation (LLA) [2,3,4,5,6,7, 45]. In Sect. 3 we extend the formalism to the unbroken \(SU(2) \times U(1)\) case taking into account that the abelian gauge boson of QED does not reggeize in the LLA. We then introduce a vacuum expectation value for the Higgs field to break the symmetry and develop the high energy effective action for the different kinematical components of the gauge fields, treating the charged and neutral sectors on the same footing. In Sect. 4 we extract the associated Feynman rules which will allow us to evaluate the electroweak gauge bosons Regge trajectories in Sect. 5 and the effective production vertices, needed to construct the integral equations proposed in [1], in Sect. 6. We draw some Conclusions after this.

2 Effective action for the SU(N) sector

The analysis of the SU(N) sector is equivalent to that of QCD without fermions, which can be neglected in the LLA. This is also relevant in the electroweak case since in the Higgs sector of the standard model, with a spontaneously broken SU(2) symmetry, the non-abelian gauge bosons also reggeize. The corresponding Lagrangian reads

$$\begin{aligned} \mathcal{L}_\mathrm{SU(N)}[v_\mu ]= & {} \frac{1}{2} \mathrm{tr}\left[ G_{\mu \nu } G^{\mu \nu } \right] \, , \end{aligned}$$
(1)

where \(D_\mu = \partial _\mu + g v_\mu \), \(G_{\mu \nu } = [D_\mu , D_\nu ]/g\), and the vector fields \(v_\mu = -i t^a v_\mu (x)\) are considered as anti-hermitian matrices acting in the SU(N) fundamental representation with generators \(t^a\) fulfilling \(\mathrm{tr}(t^a t^b) = \delta ^{ab} / 2\).

To implement the effective description of multi-Regge kinematics (MRK), which is the dominant region governing high energy scattering in the leading logarithmic approximation (LLA), we first present Fig. 1 as a reference for the momenta conventions (in this case for a \(2 \rightarrow 4\) process).

Fig. 1
figure 1

Kinematics of the \(2 \rightarrow 4\) scattering process

Clusters of particles are represented by emissions with momenta \(k_i\) in the graph. In the Regge limit we have \(s=(p_A+p_B)^2 = (2E)^2 \gg -t\). In the LLA only one particle is allowed per cluster, with two in the next-to-leading order approximation (NLLA) and so on. Each cluster is emitted at rapidity

$$\begin{aligned} \eta _i&= \frac{1}{2} \ln \frac{k_i^+}{k_i^-} \end{aligned}$$
(2)

where \( p_{A,B}^\pm n_\pm /2 = p_{A,B}\) and \(k_\pm = k \cdot n_\pm = k^\mp = k \cdot n^\mp \); furthermore \((n^\pm )^2 = 0\) and \(n^\pm \cdot n^\mp = 2\). In the LLA we resum \(g^2 \ln {(s/s_0)}\) terms (with \(s_0\) a suitable, process dependent reference scale) and the exchange in the t-channel of vector bosons is the dominant contribution with fermions and scalars being relevant only in the NLLA and beyond.

The target is to describe in a gauge invariant way the interaction of particles locally in rapidity within each cluster and non-locally with reggeized vector bosons exchanged among different clusters. The separation into different modes is realised using the following decomposition into local and non-local components with each rapidity gap sector (with no emissions) denoted by s, i.e.

$$\begin{aligned} v_\mu (x)= & {} \sum _{s} v^{(s)}_\mu (x) ~=~ \sum _{s} \left[ V_\mu ^{\mathrm{local} (s)}(x) \right. \nonumber \\&\left. + \frac{n^+}{2} V_+^{\eta > \eta (s)} + \frac{n^-}{2} V_-^{\eta < \eta (s)} \right] \, . \end{aligned}$$
(3)

This separation translates into \(S^{(s)} = \int d^4 x \mathcal{L}_{SU(N)} [v^{(s)}_\mu (x)]\), the corresponding effective action for each local sector, where bilinear terms generating the transition from a single local into a single non-local field do not contribute. This cancelation can be made explicit in the effective action introducing the Sudakov decompositions

$$\begin{aligned}&v^{(s)}_\mu (x) = v^{(s)}_{-} (x) \frac{n^-_\mu }{2} + v^{(s)}_{+} (x) \frac{n^+_\mu }{2} + v^{(s)}_{\perp } (x) \, , \end{aligned}$$
(4)
$$\begin{aligned}&V_\mu ^{\mathrm{local} (s)}(x) = V_-^{\mathrm{local} (s)}(x) \frac{n^-_\mu }{2}\nonumber \\&+ V_+^{\mathrm{local} (s)}(x) \frac{n^+_\mu }{2} + V_\perp ^{\mathrm{local} (s)}(x) \, , \end{aligned}$$
(5)

and directly subtracting the bilinear terms using

$$\begin{aligned} S^{(s)}= & {} \int d^4 x \mathcal{L}_{SU(N)} [v^{(s)}_\mu (x)]\nonumber \\&+ \int d^4 x \, \mathrm{tr}\left[ \left( v_- - V_-^{\eta< \eta (s)} \right) \partial ^2 V_+^{\eta> \eta (s)} \right] \nonumber \\&+ \int d^4 x \, \mathrm{tr}\left[ \left( v_+ - V_+^{\eta > \eta (s)} \right) \partial ^2 V_-^{\eta < \eta (s)} \right] \, , \end{aligned}$$
(6)

where

$$\begin{aligned} v^{(s)}_{-} - V_-^{\eta < \eta (s)} ~=~ V_-^{\mathrm{local} (s)} \, \, v^{(s)}_{+} - V_+^{\eta > \eta (s)} ~=~ V_+^{\mathrm{local} (s)}.\nonumber \\ \end{aligned}$$
(7)

The next step is to replace the fields \( V_+^{\eta > \eta (s)} , V_-^{\eta < \eta (s)} \) by two new fields \(A_\pm \) which are invariant under local gauge transformations and subject to the constraints \(\partial _+ A_- = 0 = \partial _- A_+\) (for convenience we remove the s index from now on) to impose the \(\eta < \eta (s)\) \((\eta > \eta (s))\) rapidity modes separation. We can therefore write

$$\begin{aligned} S= & {} \int d^4 x \mathcal{L}_{SU(N)} [v_\mu (x)]\nonumber \\&+ \int d^4 x \, \mathrm{tr}\left[ \left( v_- - A_- \right) \partial ^2 A_+ \right] \nonumber \\&+ \int d^4 x \, \mathrm{tr}\left[ \left( v_+ - A_+ \right) \partial ^2 A_-\right] \, . \end{aligned}$$
(8)

The fields \(A_\pm \) correspond to the reggeized vector bosons in the theory. It is finally needed to introduce the replacement

$$\begin{aligned} v_\pm (x) \rightarrow v_\pm (x) U[v_\pm (x)] = -\frac{1}{g} \partial _\pm U[v_\pm (x)] \end{aligned}$$
(9)

where the path-ordered Wilson line

$$\begin{aligned} U[v(x)]&= \mathcal {P} \exp \left( -\frac{g}{2} \int _{-\infty }^{x^\pm } d {z}^\pm v_\pm (z) \right) \end{aligned}$$
(10)

ensures the full gauge invariance in the effective action of Eq. (8).

3 Effective action for the electroweak interaction

When constructing the high energy effective action for the electroweak sector of the standard model it is important to take into account that not all vector fields do actually reggeize. Besides this point, the procedure will be similar to the previous section and the vector fields \(V_\mu = W_\mu , B_\mu \) of the \(SU(2) \times U(1)\) theory will be expanded into local and non-local rapidity modes.

We will subtract those contributions coupling local and non-local fields by means of an induced term in the effective action. The non-local fields will be promoted to gauge invariant high energy fields where not all of them will correspond to reggeized fields, this is distinctly different to the QCD case where both quarks and gluons do reggeize in the high energy limit. Finally, suitable path-ordered Wilson lines will be introduced in order to complete the construction of a gauge invariant high energy effective action for the electroweak sector in the standard model.

To implement this program, let us first introduce the \(SU(2) \times U(1)\) gauge fields as

$$\begin{aligned} w_\mu (x)&= -i t^a W_\mu ^a(x) \, ,&b_\mu (x)&= -i B_\mu (x) \frac{1}{2}\begin{pmatrix} 1 &{} 0 \\ 0 &{} 1 \end{pmatrix} \, , \end{aligned}$$
(11)

where \(t^a\) are the generators of the fundamental representation of SU(2) with \(\mathrm{tr}(t^at^b)=1/2\). The covariant derivatives read

$$\begin{aligned} D[w]_\mu = \partial _\mu + g w_\mu \quad D[b]_\mu =\partial _\mu + g' b_\mu \, , \end{aligned}$$
(12)

and the field strength tensors are

$$\begin{aligned} F_{\mu \nu }^a[W]&= \partial _\mu W_\nu ^a - \partial _\nu W_\mu ^a + g f^{abc} W^b_\mu W^c_\nu \,\nonumber \\ F_{\mu \nu }[B]&= \partial _\mu B_\nu - \partial _\nu B_\mu \nonumber \\ F_{\mu \nu }^a[W]t^a&= \frac{1}{-ig}[D[w]_\mu , D[w]_\nu ] \nonumber \\ F_{\mu \nu }[B]\frac{1}{2}\begin{pmatrix} 1 &{} 0 \\ 0 &{} 1 \end{pmatrix}&= \frac{1}{-ig'}[D[b]_\mu , D[b]_\nu ] \, , \end{aligned}$$
(13)

where

$$\begin{aligned}&F_{\mu \nu }[w] \equiv \frac{1}{g}[D[w]_\mu , D[w]_\nu ],\nonumber \\&F_{\mu \nu }[b] \equiv \frac{1}{g'}[D[b]_\mu , D[b]_\nu ], \end{aligned}$$
(14)

are understood as \(2\times 2\) matrices. The kinetic terms for the gauge fields are

$$\begin{aligned}&-\frac{1}{4} F[W]^a_{\mu \nu }F[W]^{a,{\mu \nu }} -\frac{1}{4} F[B]_{\mu \nu }F[B]^{{\mu \nu }}\nonumber \\&\quad = \frac{1}{2}\mathrm{tr}\left( \frac{1}{2} F[w]_{\mu \nu }F[w]^{{\mu \nu }} + F[b]_{\mu \nu } F[b]^{\mu \nu } \right) \, . \end{aligned}$$
(15)

The covariant derivative associated to the Higgs field considered as a SU(2) scalar complex doublet, \(\phi \), reads

$$\begin{aligned} D_\mu \phi (x)&= \left( \partial _\mu + g w_\mu + g' b_\mu \right) \phi (x) \, . \end{aligned}$$
(16)

Together with the Higgs potential \(V(\phi ) = \mu ^2 \phi ^\dagger \phi - \lambda (\phi ^\dagger \phi )^2\) this yields the well-known electroweak Lagrangian (neglecting the coupling to fermions):

$$\begin{aligned} \mathcal {L}_\mathrm{EW}[w_\mu , b_\mu , \phi ]&= \frac{1}{2}\mathrm{tr}\left( F[b]_{\mu \nu } F[b]^{\mu \nu } + \frac{1}{2} F[w]_{\mu \nu }F[w]^{{\mu \nu }} \right) \nonumber \\&\quad + \left| D_\mu \phi (x)\right| ^2 + V(\phi ) \, . \end{aligned}$$
(17)

The separation into different rapidity modes can be achieved using

$$\begin{aligned} w_\mu (x)&= \sum _{s} \left[ W_\mu ^{local, s}(x) + \frac{n^+}{2} W_+^{\eta> \eta (s)} + \frac{n^-}{2} W_-^{\eta< \eta (s)} \right] \nonumber \\ b_\mu (x)&= \sum _{s} \left[ B_\mu ^{local, s}(x) + \frac{n^+}{2} B_+^{\eta > \eta (s)} + \frac{n^-}{2} B_-^{\eta < \eta (s)} \right] \, , \end{aligned}$$
(18)

where

$$\begin{aligned}&\partial _- W^{\eta> \eta (s)}_+(x) = 0 = \partial _+ W^{\eta< \eta (s)}_-(x),\nonumber \\&\partial _- B^{\eta > \eta (s)}_+(x) = 0 = \partial _+ B^{\eta < \eta (s)}_-(x) \, . \end{aligned}$$
(19)

Taking into account that the Higgs field is always local in the rapidity mode expansion, before spontaneous symmetry breaking the coupling between a local and a non-local field lies within the kinetic terms of the gauge fields. The situation is more complicated when the Higgs field acquires a vacuum expectation value of the form

$$\begin{aligned} \phi _0&= \langle \phi (x)\rangle = \frac{1}{\sqrt{2}} \begin{pmatrix} 0 \\ v \end{pmatrix} \end{aligned}$$
(20)

since in the kinetic term of the Higgs field we find

$$\begin{aligned} \mathrm{tr}\left[ w_\mu g^2M^2 w_\mu +b_\mu {g'}^2 M^2 b_\mu + gg' \left( w_\mu M^2 b_\mu + b_\mu M^2 w_\mu \right) \right] \end{aligned}$$
(21)

where

$$\begin{aligned} \phi _0 \otimes \phi _0^\dagger&= \frac{v^2}{2} \begin{pmatrix} 0 &{} 0 \\ 0 &{} 1 \end{pmatrix} = M^2 \, . \end{aligned}$$
(22)

This generates local to non-local field transitions after kinematical separation in rapidity of the electroweak boson fields of the form

$$\begin{aligned}&\frac{1}{2}\mathrm{tr}\bigg [ W^{local}_+ g^2M^2 W_-^{\eta< \eta (s)} + W_-^{\eta< \eta (s)} g^2M^2 W^{local}_+\nonumber \\&\quad + B^{local}_+ {g'}^2M^2 B_-^{\eta< \eta (s)} + B_-^{\eta< \eta (s)} {g'}^2 M^2 B^{local}_+ \nonumber \\&\quad + gg'\left( W^{local}_+ M^2 B_-^{\eta< \eta (s)} + B_-^{\eta< \eta (s)} g^2M^2 W^{local}_+ \right. \nonumber \\&\quad \left. + B^{local}_+ M^2 W_-^{\eta< \eta (s)} + W_-^{\eta< \eta (s)} g^2M^2 B^{local}_+ \right) \nonumber \\&\quad + \{ ``+\text {''} \leftrightarrow ``-\text {''} , ``>\text {''} \rightarrow ``<\text {''}\} \bigg ] \, . \end{aligned}$$
(23)

As explained in the previous section in the pure SU(N) theory, these terms need to be subtracted from the original electroweak effective action in order to correctly capture the multi-Regge limit of scattering amplitudes. Naming Eq. (23) as “Higgs subtraction” term, for the effective action valid in the multi-Regge limit \( {S}^\mathrm{EW}_\mathrm{eff}\) we would have

$$\begin{aligned} {S}^\mathrm{EW}_\mathrm{eff}= & {} \int d^4 x \, \mathcal {L}_{\text {EW}}[w_\mu , b_\mu , \phi ]\nonumber \\&+ \int d^4 x \, \mathrm{tr}\left[ \left( w_- - W_-^{\eta< \eta (s)} \right) \partial ^2 W_+^{\eta> \eta (s)} \right] \nonumber \\&+ \int d^4 x \, \mathrm{tr}\left[ \left( b_- - B_-^{\eta< \eta (s)} \right) \partial ^2 B_+^{\eta> \eta (s)} \right] \nonumber \\&- \text {Higgs subtraction}\nonumber \\&+ \{ ``+{\text {''}} \leftrightarrow ``-{\text {''}} , ``>{\text {''}} \rightarrow ``<{\text {''}}\} \, . \end{aligned}$$
(24)

As in Eq. (7) we now write \(W^{local}_- = w_- - W_-^{\eta < \eta (s)} \) and \(W^{local}_+ = w_+ - W_+^{\eta > \eta (s)}\) with equivalent expressions for the B field. The next step is to promote the non-local fields to gauge invariant high energy fields \(A_\pm ^{(w)}\) and \(A_\pm ^{(b)}\) subject to the kinematical constraints \(\partial _\mp A_\pm ^{(w)} = 0 = \partial _\mp A_\pm ^{(b)}\). We therefore have an effective action written as the sum of the usual standard model electroweak action plus additional “induced” terms (using Lipatov’s notation) where the Lagrangian for high energy fields reads

$$\begin{aligned} \mathcal {L}_\mathrm{induced}= & {} \frac{1}{2} \mathrm{tr}\bigg [ (w_+ - A_+^{(w)}) \left( \partial ^2 + {g^2}M^2\right) A_-^{(w)} \nonumber \\&+ A_-^{(w)} \left( \partial ^2 + {g^2}M^2\right) (w_+ - A_+^{(w)})\nonumber \\&+ (b_+ - A_+^{(b)}) \left( \partial ^2 +{g'}^2 M^2\right) A_-^{(b)}\nonumber \\&+ A_-^{(b)} \left( \partial ^2 +{{g'}^2} M^2\right) (b_+ - A_+^{(b)})\nonumber \\&+ {gg'} \bigg \{(w_+ - A_+^{(w)}) M^2 A_-^{(b)} + A_-^{(b)} M^2 (w_+ - A_+^{(w)})\nonumber \\&+ (b_+- A_+^{(b)}) M^2 A_-^{(w)} \nonumber \\&+ A_-^{(w)} M^2 (b_+- A_+^{(b)}) \bigg \} \bigg ] + \{``+{\text {''}} \leftrightarrow ``-{\text {''}} \} \, . \end{aligned}$$
(25)

This action is not yet fully gauge invariant since, even though the abelian b field is invariant under local gauge transformations, the non-abelian w field is not. Similarly to the SU(N) case in the previous section it is needed to replace \(w_\pm \) by a path-ordered Wilson line

$$\begin{aligned} w_\pm \rightarrow -\frac{1}{g}\partial _\pm U[w_\pm (x)] \equiv T[w_\pm ] \, . \end{aligned}$$
(26)

One should therefore write

$$\begin{aligned} \mathcal {L}_\mathrm{induced}= & {} \frac{1}{2}\mathrm{tr}\bigg [ (T[w_+] - A_+^{(w)}) \left( \partial ^2 + {g^2}M^2\right) A_-^{(w)} \nonumber \\&+ A_-^{(w)} \left( \partial ^2 + {g^2}M^2\right) (T[w_+] - A_+^{(w)})\nonumber \\&+ (b_+ - A_+^{(b)}) \left( \partial ^2 +{{g'}^2} M^2\right) A_-^{(b)}\nonumber \\&+ A_-^{(b)} \left( \partial ^2 +{{g'}^2} M^2\right) (b_+ - A_+^{(b)}) \nonumber \\&+ {gg'} \bigg \{(T[w_+] - A_+^{(w)}) M^2 A_-^{(b)} \nonumber \\&+ A_-^{(b)} M^2 (T[w_+] - A_+^{(w)}) \nonumber \\&+ (b_+- A_+^{(b)}) M^2 A_-^{(w)} + A_-^{(w)} M^2 (b_+- A_+^{(b)}) \bigg \} \bigg ] \nonumber \\&+ \{``+\text {''} \leftrightarrow ``-\text {''} \} \, . \end{aligned}$$
(27)

4 Effective Feynman rules

To study the high energy behaviour of electroweak scattering amplitudes in the multi-Regge limit two ingredients are needed in the LLA: effective emission vertices and electroweak boson Regge trajectories. We will now focus on identifying the induced effective vertices generated by the high energy action which will actually contribute in the calculation of those two fundamental components of the effective field theory.

In order to calculate the electroweak gauge bosons Regge trajectories it is needed to first find the propagator of the high energy fields, i.e. to evaluate the two-point vertex of a \(A_+\) and a \(A_-\) field, \(V_{A_+,A_-}(q^2)\). To this end we express the quadratic term of the induced Lagrangian in the basis of mass eigenstates (with \(M_Z^2 = \frac{M_W^2}{\cos ^2{(\theta _W)}} = \frac{ v^2}{4} (g^2+{g'}^2 )\) where \(\theta _W\) is Weinberg’s angle), i.e.

$$\begin{aligned} \mathcal {L}_\mathrm{induced}^\mathrm{quadratic}= & {} - \frac{1}{2} \bigg \{ \left( W_+^+ - A_+^{(W^+)} \right) \left( \partial ^2 + M_W^2 \right) A_-^{(W^-)}\nonumber \\&+ \left( W_+^- - A_+^{(W^-)} \right) \left( \partial ^2 + M_W^2 \right) A_-^{(W^+)} \nonumber \\&+ \left( Z_+ - A_+^{(Z)} \right) \left( \partial ^2 + M_Z^2 \right) A_-^{(Z)} \nonumber \\&+ \left( C_+ - A_+^{(C)} \right) \partial ^2 A_-^{(C)} \bigg \} + \{``+{\text {''}} \leftrightarrow ``-{\text {''}} \} \, , \end{aligned}$$
(28)

where \(C_\mu \) is the photon field and \(R_\pm \) denotes \(n_\pm ^\mu R_\mu \) (with \(R=W^\pm ,Z,C\)). The propagator in momentum space can then be obtained making use of

$$\begin{aligned} G_{A_+A_-}(q^2)&= \frac{-1}{V_{A_+A_-}(q^2)}\bigg |_{q^2 \rightarrow q^2 + i \epsilon } \end{aligned}$$
(29)

where the vertices receive both contributions from the conventional electroweak Lagrangian and the induced Lagrangian. The full set of relevant vertices and propagators for the high energy fields (not reggeized degrees of freedom yet) reads

$$\begin{aligned} V_{A_+^{(W^+)}A_-^{(W^+)}}(q^2)= & {} \frac{i}{2} (-q^2 + M_W^2) \, ,\nonumber \\ G_{A_+^{(W^+)}A_-^{(W^-)}}= & {} \frac{-1}{V_{A_+^{(W^+)}A_-^{(W^-)}} (q^2 + i\epsilon )} = \frac{-2i}{q^2 - M_W^2 + i \epsilon } \, , \nonumber \\ V_{A_+^{(W^-)}A_-^{(W^+)}}(q^2)= & {} \frac{i}{2} (-q^2 + M_W^2) \nonumber \\ G_{A_+^{(W^-)}A_-^{(W^+)}}= & {} \frac{-1}{V_{A_+^{(W^-)}A_-^{(W^+)}} (q^2 + i\epsilon )} = \frac{-2i}{q^2 - M_W^2 + i \epsilon } \, ,\nonumber \\ V_{A_+^{(Z)}A_-^{(Z)}}(q^2)= & {} \frac{i}{2} (-q^2 + M_Z^2) \nonumber \\ G_{A_+^{(Z)}A_-^{(Z)}}= & {} \frac{-1}{V_{A_+^{(Z)}A_-^{(Z)}} (q^2 + i\epsilon )} = \frac{-2i}{q^2 - M_Z^2 + i \epsilon } \, ,\nonumber \\ V_{A_+^{(\gamma )}A_-^{(\gamma )}}(q^2)= & {} \frac{i}{2} (-q^2) \nonumber \\ G_{A_+^{(\gamma )}A_-^{(\gamma )}}= & {} \frac{-1}{V_{A_+^{(\gamma )}A_-^{(\gamma )}} (q^2 + i\epsilon )} = \frac{-2i}{q^2 + i \epsilon } \, . \end{aligned}$$
(30)

As the next step we need to determine the first induced vertex, which corresponds to the interaction of three fields in the induced Lagrangian. With this target in mind, it is convenient to switch to a different SU(2) algebra basis. We define

$$\begin{aligned} t^\pm&= \frac{1}{\sqrt{2}} (t^1 \pm i t^2) \, ,&[t^+, t^-]&= -t^3 \, ,&[t^\pm , t^3]&= \pm t^\pm \, , \end{aligned}$$
(31)

which implies \(i W_\mu = W^+_\mu t^+ + W_\mu ^- t^- + W_\mu ^3 t^3\). Next, it is needed to study the first non-trivial term of the electroweak Wilson line \(T[w_+]\) in Eq. (26):

$$\begin{aligned} -gw_+ \frac{1}{\partial _+} w_+= & {} - \frac{g}{2} \left[ w_+ \frac{1}{\partial _+} w_+ - \left( \frac{1}{\partial _+} w_+ \right) w_+ \right] \nonumber \\= & {} \frac{g}{2} W_+^a \left( \frac{1}{\partial _+} W^b_+ \right) \cdot [t^a, t^b] \nonumber \\= & {} \frac{g}{2} \left[ W_+^+ \left( \frac{1}{\partial _+} W^-_+ \right) - W_+^- \left( \frac{1}{\partial _+} W^+_+ \right) \right] (-t^3) \nonumber \\&+\frac{g}{2} \left[ W_+^+ \left( \frac{1}{\partial _+} W^3_+ \right) - W_+^3 \left( \frac{1}{\partial _+} W^+_+ \right) \right] t^+ \nonumber \\&+ \frac{g}{2} \left[ W_+^3 \left( \frac{1}{\partial _+} W^-_+ \right) - W_+^- \left( \frac{1}{\partial _+} W^3_+ \right) \right] t^- \, , \end{aligned}$$
(32)

where \(W^3_+ = \left( g Z_+ + g' C_+ \right) / \sqrt{g^2 + {g'}^2} = \cos {(\theta _W)} Z_+ + \sin {(\theta _W)} C_+\). The relevant induced vertices stem from

$$\begin{aligned} \mathcal {L}_{ind}^{\mathcal {O}(g)}&= \frac{1}{2} \mathrm{tr}\left\{ -gw_+ \frac{1}{\partial _+} w_+ \left[ \partial ^2 A_-^{(w)} + M^2 \left( {g^2} A_-^{(w)} + {gg'} A_-^{(b)} \right) \right] \right\} \nonumber \\&\quad + \frac{1}{2} \mathrm{tr}\left\{ \left[ A_-^{(w)} \partial ^2 + \left( {g^2} A_-^{(w)} + {gg'} A_-^{(b)} \right) M^2\right] \left( -gw_+ \frac{1}{\partial _+} w_+ \right) \right\} . \end{aligned}$$
(33)

To evaluate the traces it is enough to study the projection of the term in the squared bracket on the three different generators. One finds

$$\begin{aligned}&i \, \mathrm{tr}\left[ \left\{ \left[ A_-^{(w)} \partial ^2 + \left( {g^2} A_-^{(w)} + {gg'} A_-^{(b)} \right) M^2\right] , t^+ \right\} \right] \nonumber \\&\quad = \left( \partial ^2 + M_W^2 \right) A_-^{{(W^-)}} i \, \mathrm{tr}\left[ \left\{ \left[ A_-^{(w)} \partial ^2 + \left( {g^2} A_-^{(w)} + {gg'} A_-^{(b)} \right) M^2\right] , t^- \right\} \right] \nonumber \\&\quad = \left( \partial ^2 + M_W^2 \right) A_-^{{(W^+)}} \, i \, \mathrm{tr}\left[ \left\{ \left[ A_-^{(w)} \partial ^2 + \left( {g^2} A_-^{(w)} + {gg'} A_-^{(b)} \right) M^2\right] , t^3 \right\} \right] \nonumber \\&\quad = \cos {(\theta _W)} \left( \partial ^2 + M_Z^2 \right) A_-^{{(Z)}} + \sin {( \theta _W)} \partial ^2 A_-^{{(\gamma )}} \, . \end{aligned}$$
(34)

The remaining task is to write down the corresponding momentum space expressions.

The two induced vertices of interest in our calculations can be depicted as follows,

figure a

where \(i,j,k = W^\pm , Z, \gamma ,\) and, from the induced Lagrangian, we obtain

$$\begin{aligned} K\left[ A_\pm ^{W^-}(q), W_{\nu _1}^+(k_1), Z_{\nu _2}(k_2) \right]&=K\left[ A_\pm ^{W^+}(q), Z_{\nu _1}(k_1), W^-_{\nu _2}(k_2) \right] \nonumber \\&=\cos \theta _W ({\varvec{q}}^2 + M_W^2) \, , \nonumber \\ K\left[ A_\pm ^{W^-}(q), W_{\nu _1}^+(k_1), \gamma _{\nu _2}(k_2) \right]&= K\left[ A_\pm ^{W^+}(q), \gamma _{\nu _1}(k_1), W^-_{\nu _2}(k_2) \right] \nonumber \\&= \sin \theta _W ({\varvec{q}}^2 + M_W^2),\nonumber \\ K\left[ A_\pm ^{Z}(q), W_{\nu _1}^+(k_1), W^-_{\nu _2}(k_2) \right]&=({\varvec{q}}^2 + M_Z^2) \cos \theta _W \, , \nonumber \\ K\left[ A_\pm ^{\gamma }(q), W_{\nu _1}^+(k_1), W^-_{\nu _2}(k_2) \right]&= {\varvec{q}}^2 \sin \theta _W \, . \end{aligned}$$
(35)

For clarity, field polarizations and momenta have been written explicitly.

5 Reggeization of the electroweak bosons

In the charged sector, the \(W^\pm \) bosons are known to lie on their Regge trajectories in the broken symmetry case, i.e. the Born amplitude for the scattering of W bosons of a given helicity and mass \(M_W\) has the all-orders form in the LLA (with \(1+ \omega _\mathrm{charged}(t)\) being the W boson Regge trajectory)

$$\begin{aligned} A \simeq \frac{A_\mathrm{Born}}{2} \left( \left( \frac{s}{M_W^2}\right) ^{\omega _\mathrm{charged}(t)}+\left( -\frac{s}{M_W^2}\right) ^{\omega _\mathrm{charged}(t)}\right) .\nonumber \\ \end{aligned}$$
(36)

To calculate the Regge trajectory it is needed to evaluate one-loop corrections to the amplitude for a charged isospin-1 exchange and look for rapidity divergent contributions (which are due to the kinematics, not to short distance dynamics). The starting point is the bare propagator in unitary gauge

$$\begin{aligned}&{G}_\mathrm{charged}^B\left( \rho \right) = \frac{2i}{{\varvec{q}}^2 + M_W^2}\nonumber \\&\qquad \times \left\{ 1 + \frac{2i}{{\varvec{q}}^2 + M_W^2} \Sigma _\mathrm{charged}\left( \rho , \epsilon , {\varvec{q}}^2 \right) + \cdots \right\} \, , \end{aligned}$$
(37)

where we focus on the high energy divergent piece, i.e.

$$\begin{aligned} \Sigma _\mathrm{charged}\left( \rho , \epsilon , {\varvec{q}}^2 \right)= & {} -i K_1^2 \frac{g^2}{(4 \pi )^{d/2}} \int \frac{d^d l}{i \pi ^{d/2}} \frac{n^\mu _+ n^\nu _+n^{\mu '}_+n^{\nu '}_+}{2l^+ 2l^-}\nonumber \\&\times G^{\mu \nu }_{W^\pm }(l) G^{\mu '\nu '}_{Z}(k+l) - i K_2^2 \frac{g^2}{(4 \pi )^{d/2}} \nonumber \\&\times \int \frac{d^d l}{i \pi ^{d/2}} \frac{n^\mu _+ n^\nu _+n^{\mu '}_+n^{\nu '}_+}{2l^+ 2l^-} G^{\mu \nu }_{W^\pm }(l) G^{\mu '\nu '}_{\gamma }(k+l) \nonumber \\\simeq & {} - \rho \frac{ ({\varvec{q}}^2 + M_W^2)^2}{2i} \left[ \beta _{WZ}({\varvec{q}}^2) \cos ^2{(\theta _W)}\right. \nonumber \\&+\left. \beta _{W\gamma }({\varvec{q}}^2) \sin ^2{(\theta _W)} \right] , \end{aligned}$$
(38)

where

$$\begin{aligned} G^{\mu \nu }_A(l)&= \frac{-i}{l^2 - m_A^2} \left( g^{\mu \nu } - \frac{l^\mu l^\nu }{m_A^2} \right) \, , \end{aligned}$$
(39)

\(p = n_+ + e^{-\rho } n_-\), \(n = n_- + e^{-\rho } n_+\), \(p \cdot n \simeq 2\), \(p^2 = n^2 = 4 e^{-\rho }\) and

$$\begin{aligned} \rho \simeq \frac{1}{2} \log \left( \frac{ 4 n \cdot p^2}{n^2 p^2}\right) \, . \end{aligned}$$
(40)

We have used the notation of [1] (together with \(g^2 = 4 \pi \alpha _{ew}, \bar{\alpha }_{ew} = 2 \alpha _{ew} / \pi \)) for

$$\begin{aligned} \beta _{ij}({\varvec{q}}^2)&= 2\frac{g^2}{(4\pi )^2} \int \frac{d^2{\varvec{l}}}{\pi } \frac{1}{[{\varvec{l}}^2 + M_i^2] [{\varvec{l}} + {\varvec{q}}^2 + M_j^2]} \, . \end{aligned}$$
(41)

We can then write

$$\begin{aligned} {G}^B_\mathrm{charged} \left( \rho \right)&= \frac{2i}{{\varvec{q}}^2 + M_W^2} - 2 i \rho \left[ \beta _{WZ}({\varvec{q}}^2) \cos ^2{(\theta _W)}\right. \nonumber \\&\quad +\left. \beta _{W\gamma }({\varvec{q}}^2) \sin ^2{(\theta _W)} \right] \, . \end{aligned}$$
(42)

In the LLA this expression is enough to show agreement with the results of [1] for the reggeization of the charged electroweak bosons with trajectory

$$\begin{aligned}&\omega _\mathrm{charged} (t) = (t-M_W^2)\nonumber \\&\times \left( \beta _{WZ}(-t) \cos ^2{(\theta _W)} + \beta _{W\gamma }(-t) \sin ^2{(\theta _W)} \right) \, . \end{aligned}$$
(43)

At higher orders it is needed to introduce the renormalization procedure proven to be valid at two loops in the case of the high energy QCD effective action and put forward in Refs. [11,12,13,14,15,16,17,18,19,20].

Fig. 2
figure 2

Production of a photon or Z boson (left) and a Higgs boson (right) from a charged W vector boson

The treatment of the neutral sector is more complicated since the Z and \(\gamma \) bosons do not reggeize within the Stantard Model (they do in grand unified theories [46, 47]) and a representation such as Eq. (36) does not apply. In this case both the Z and \(\gamma \) high energy fields can split into a \(W^+ W^-\) pair and mix in the t-channel. This fact generates four possible bare propagators which at one loop, and keeping only the high energy divergent contributions (proportional to \(\rho \sim \ln {s}\)), can be written (using the same notations as above) in the form

$$\begin{aligned} \widehat{G}_\mathrm{neutral}^B (\rho ) = \frac{1}{2i} \begin{pmatrix} {G}_{ZZ}^B &{} {G}_{Z\gamma }^B \\ {G}_{\gamma Z}^B &{} {G}_{\gamma \gamma }^B \end{pmatrix} = \widehat{G}_0 - \rho \, \beta _{WW}({\varvec{q}}^2) \, \widehat{G}_1 \, , \end{aligned}$$
(44)

where

$$\begin{aligned}&\widehat{G}_0 = \widehat{G}_\mathrm{neutral}^B (0) = \begin{pmatrix} \frac{1}{{\varvec{q}}^2 + M_Z^2} &{} 0 \\ 0 &{} \frac{1}{{\varvec{q}}^2} \end{pmatrix} \, , \end{aligned}$$
(45)

and

$$\begin{aligned} \widehat{G}_1= & {} \begin{pmatrix} \cos ^2{(\theta _W)} &{} \frac{\sin {(2 \theta _W)}}{2} \\ \frac{\sin {(2 \theta _W)}}{2} &{} \sin ^2{(\theta _W)} \end{pmatrix} \nonumber \\= & {} \frac{{\varvec{q}}^2 + M_W^2}{2} \left[ \widehat{S} \begin{pmatrix} 1&{} 0 \\ 0 &{} 0 \end{pmatrix} \widehat{S}^{-1} \widehat{G}_0 + \widehat{G}_0 \left( \widehat{S}^{-1}\right) ^T \begin{pmatrix} 1&{} 0 \\ 0 &{} 0 \end{pmatrix} \widehat{S}^{T} \right] ,\nonumber \\ \end{aligned}$$
(46)

with

$$\begin{aligned}&\widehat{S} = \begin{pmatrix} \cos {(\theta _W)}&{} - \frac{{\varvec{q}}^2}{{\varvec{q}}^2 + M_Z^2}\sin {(\theta _W)} \\ \sin {(\theta _W)} &{} \cos {(\theta _W)} \end{pmatrix} \nonumber \\&\widehat{S}^{-1} = \frac{{\varvec{q}}^2 + M_Z^2 }{{\varvec{q}}^2 + M_W^2 } \begin{pmatrix} \cos {(\theta _W)}&{} \frac{{\varvec{q}}^2}{{\varvec{q}}^2 + M_Z^2}\sin {(\theta _W)} \\ -\sin {(\theta _W)} &{} \cos {(\theta _W)} \end{pmatrix} \, . \end{aligned}$$
(47)

This result is in agreement with those presented in [1] where there exists a Regge trajectory,

$$\begin{aligned} \omega _\mathrm{neutral} (t) = (t-M_W^2)\beta _{WW} (-t) \, , \end{aligned}$$
(48)

in the neutral sector. However, this trajectory \(\alpha _\mathrm{neutral} (t) = 1+ \omega _\mathrm{neutral} (t)\) crosses one at \(t=M_W^2\) which implies that the neither the Z boson nor the photon lie on it.

6 Emission vertices

The second ingredient needed to construct the integral equations governing the high energy behaviour of electroweak scattering amplitudes in the multi-Regge limit is the set of real emission vertices with electroweak bosons or a Higgs in the final state.

The effective action proposed in this work allows to obtain the same central rapidity production vertices as those calculated by Bartels, Lipatov and Peters in the LLA in Ref. [1]. In the case of the exchange of a high energy W boson in the t-channel with associated production of a photon the three Feynman diagrams in the effective field theory are shown in Fig. 2 (left).

The sum of these diagrams renders the expression for the vertex

$$\begin{aligned} C_\mu ^{W^+W^-\gamma }&= \frac{ie}{2} \left[ -q_1^+ n_\mu ^- + q_2^- n_\mu ^+ + (q_1 - q_2)^\mu \right] \nonumber \\&\quad - \frac{ie}{2} \frac{q_1^2 - M_W^2}{q_2^-} n_\mu ^- + \frac{ie}{2} \frac{q_2^2 - M_W^2}{q_1^+} n_\mu ^+ \, . \end{aligned}$$
(49)

There is a one-to-one correspondence between the three terms in this vertex and each of the graphs in the figure.

The corresponding vertex for the exchange of a Z boson can be calculated by making the replacement \(e \rightarrow g \cos {(\theta _W)} = e \cot {(\theta _W)}\).

The graph for Higgs production shown in Fig. 2 (right) corresponds to the vertex

$$\begin{aligned} C^{W^+W^-H}&= 2i \frac{m_W^2}{v} \, . \end{aligned}$$
(50)

The remaining production vertices are calculated with the same method. Let us fix our notation for the above mentioned vertex \(C_\mu ^{W^+W^-\gamma }\) (see Fig. 2(left)) as \(r_+^{W^+}(q_1) + r_-^{W^-} (q_2)\rightarrow \gamma (q_1 + q_2)\), etc. (where “\(r_\pm ^A (q)\)” stands for a Reggeon A exchanged in the t-channel with + (-) lightcone projection and transverse momentum q). We then find for \(r_+^{\gamma }(q_1) + r_-^{W^+} (q_2)\rightarrow W^+ (q_1 + q_2) \):

$$\begin{aligned} C_\mu ^{\gamma W^+W^+}&= \frac{ie}{2} \left[ q_1^+ n_\mu ^- - q_2^- n_\mu ^+ + (q_2 - q_1)^\mu \right] \nonumber \\&+ \frac{ie}{2} \frac{q_1^2}{q_2^-}n_\mu ^- -\frac{ie}{2} \frac{q_2^2 - M_W^2}{q_1^+} n_\mu ^+ \, , \end{aligned}$$
(51)

and for \(r_+^{Z}(q_1) + r_-^{W^+} (q_2)\rightarrow W^+ (q_1 + q_2) \):

$$\begin{aligned} C_\mu ^{Z W^+W^+}&=\frac{ig \cos \theta _W}{2} \left[ q_1^+ n_\mu ^- - q_2^- n_\mu ^+ + (q_2 - q_1)^\mu \right] \nonumber \\&\quad + \frac{ig \cos \theta _W}{2} \frac{q_1^2 - M_Z^2}{q_2^-} n_\mu ^-\nonumber \\&-\frac{ig \cos \theta _W}{2} \frac{q_2^2 - M_W^2}{q_1^+} n_\mu ^+ . \end{aligned}$$
(52)

These vertices are equivalent to those calculated in [1]. This concludes our analysis of the effective action in the LLA since with the results of the last two sections it is possible to reproduce, using the product of pairs of these effective vertices and applying s-channel unitarity, both the vacuum and isospin-1 exchange integral equations derived in [1]. The detailed study of the solutions to those equations (see [48] for some studies in this direction) and the application of this formalism at NLL will be defered to a future work.

7 Conclusions

We have introduced an effective action valid to evaluate scattering amplitudes in the high energy limit of the electroweak sector of the standard model. We have made use of Wilson lines to deal with reggeized charged electroweak vector bosons and investigated their effective interactions among themselves and with the non-reggeizing neutral vector bosons. We have found agreement with the results of [1] for the charged electroweak vector boson Regge trajectories and production vertices needed to write the evolution equations governing the singlet and isospin-1 exchanges in the t-channel in multi-Regge kinematics. The results here presented belong to the leading logarithmic approximation but the proposed high energy effective action is valid beyond this. We will evaluate the next-to-leading order corrections to these results in future works making use of renormalized high energy propagators. This method has already been proven to be efficient at that order in the QCD case [17]. Probably the most interesting feature of the next-to-leading order calculation is to test if the reggeization in the neutral sector could already appear at that level of accuracy. In the effective theory context this would take place in the renormalization procedure mandatory for the associated propagators. A further interesting aspect of the calculation is to evaluate the likely presence of double logarithms related to the choice of energy scale in the resummed terms. The identification of these terms will, most likely, lead to their resummation to all orders using very similar arguments to those in [49].