Effective field theory versus UV-complete model: vector boson scattering as a case study

Effective field theories (EFT) are commonly used to parameterize effects of BSM physics in vector boson scattering (VBS). For Wilson coefficients which are large enough to produce presently observable effects, the validity range of the EFT represents only a fraction of the energy range covered by the LHC, however. In order to shed light on possible extrapolations into the high energy region, a class of UV-complete toy models, with extra SU(2) multiplets of scalars or of fermions with vector-like weak couplings, is considered. By calculating the Wilson coefficients up to energy-dimension eight, and full one-loop contributions to VBS due to the heavy multiplets, the EFT approach, with and without unitarization at high energy, is compared to the perturbative prediction. For high multiplicities, e.g. nonets of fermions, the toy models predict sizable effects in transversely polarized VBS, but only outside the validity range of the EFT. At lower energies, dimension-eight operators are needed for an adequate description of the models, providing another example that dimension-eight can be more important than dimension-six operators. A simplified VBFNLO implementation is used to estimate sensitivity of VBS to such BSM effects at the LHC. Unitarization captures qualitative features of the toy models at high energy but significantly underestimates signal cross sections in the threshold region of the new particles.


Introduction
With over 300 fb −1 of data collected at the LHC, vector boson scattering processes (VBS) have become important testing grounds for the non-abelian structure of electroweak interactions, as described by the Standard Model (SM). While early measurements of the CMS and ATLAS collaborations firmly established the presence of these electroweak qq → qqV V signals [1,2], at a level as predicted by the SM, the focus is now shifting to precise comparisons between data and theoretical predictions and the search for possible signals beyond the SM. Possible deviations in the three-and four-vector-boson-couplings as compared to SM predictions, so called anomalous triple gauge couplings (aTGC) or anomalous quartic gauge couplings (aQGC) have received particular attention [2,3].
Anomalous couplings are conveniently described within an effective field theory (EFT) framework [4][5][6][7] with operators of energy dimension d = 6 and 8 such as Here Φ represents the Higgs doublet field (Φ = (0, v + h)/ √ 2 in the unitary gauge),Ŵ µν is the SU(2) field strength tensor, and the scale of new physics is generically called Λ. For a process with typical momentum transfer Q 2 , deviations from SM amplitudes will generally scale like (Q/Λ) d−4 or (v/Λ) d−4 within the validity range of the EFT, which is set by Q 2 < Λ 2 , and, thus, dimension-6 operators are expected to dominate at the LHC. The dynamics of the new physics may, however, suppress the Wilson coefficients f i /Λ 2 of dimension-6 operators such that the leading effects arise from dimension-8 operators. A simple example is the Euler-Heisenberg Lagrangian [8], where a term analogous to f T 0 in Eq. (1) exists, due to the electron loop producing a 4-photon interaction, while a f W W W -like term is forbidden by Furry's theorem, i.e. the fact that the photon is C-odd. As we will see later, such a suppression of aTGCs vs. aQGCs, due to Furry's theorem, does indeed occur for simple extensions of the SM.
A second consideration is whether the operators in the EFT Lagrangian are loop induced or whether they may occur due to tree level effects of the underlying UV-complete model. While terms such as f S 0 in Eq. (1) can arise from tree level scalar exchange in e.g. a two-Higgs-doublet model, the appearance of field strength tensors in the effective Lagrangian generally points to a loop origin of the effective interaction. The concomitant suppression by loop factors g 2 /(16π 2 ) tends to substantially reduce the phenomenological importance of such operators [9], at least when they are competing with SM tree level effects. 1 One might thus be tempted to ignore any loop induced operators in the study of VBS, arguing that loop induced effects are too small to be observed above sizable SM contributions to VBS. We will show that this assumption, while warranted for a large class of models, can be false when high multiplicities of the additional fields occur in the UV physics.
New charged particles are expected to first appear in direct pair production at colliders, or, indirectly, in the simplest n-point functions which can be probed. Already three decades ago, de Rujula and collaborators asked "Does LEP-1 obviate LEP-2?" [10], i.e. does the nonobservation of new physics in 2-point functions, via the S, T , U parameters of LEP-1, exclude the observation of anomalies in W -pair production at LEP-2. Translated to the present, the question arises whether the measurement of quartic couplings in VBS can lead to a first BSM signal, even though theqq → V V vector boson pair production processes are observed at the LHC with much higher statistics and are also known better theoretically, with the availability of e.g. NNLO QCD corrections [11][12][13][14][15]. An additional scalar resonance would be an obvious example for the power of VBS. However, we will show that for large isospin multiplets also purely transverse aQGCs, like the f T term in Eq. (1), may be more sensitive probes of new physics than aTGCs, because aQGCs are enhanced by an additional factor of isospin-squared. Of course, very large isospin representations will push models out of the perturbative domain, a tension which we will probe by unitarity considerations.
Typically, Wilson coefficients of dimension-8 operators, which are large enough to be observable in present data, will produce scattering amplitudes which exceed unitarity bounds within the energy range of the LHC. Not taking into account the resulting unitarity limits on cross sections results in overly stringent constraints on Wilson coefficients. Accordingly, recent LHC analyses take unitarity constraints into account in their extraction of limits on aQGC (see e.g. Refs. [2,3]). Unfortunately, unitarization procedures for VBS amplitudes such as the T-matrix approach of Refs. [16][17][18] or the T u -model of Ref. [19] are arbitrary to a significant degree, introducing model dependence in the measurement of the Wilson coefficients of an EFT. A comparison of unitarized cross sections with predictions of UV-complete toy models can help to make the unitarization procedures more realistic.
In this paper we address the above questions by confronting the EFT description of the transverse, dimension-6 f W W W -and dimension-8 f T -type operators in Eq. (1) with a fairly generic class of UV-complete models with additional matter fields in arbitrary weak isospin representations. The model, given by the addition to the SM of (2J R + 1)-dimensional SU(2) multiplets of heavy scalars and/or heavy non-chiral Dirac fermions, is further specified in Section 2.1. At the one-loop level, which is summarized in Section 2.2, its low energy manifestation is described by an EFT Lagrangian involving (derivatives of) SU(2) field strengths, i.e. only transverse operators contribute. The relevant operators are introduced in Section 2.3. In Section 2.4 we present the Wilson coefficients of all EFT terms up to energy dimension eight. Many of them have been measured in various experiments and the existing bounds are translated into exclusions of (J R , M R ) combinations in Section 2.5, with M R denoting the mass of the heavy scalars or fermions. The resulting amplitudes for VBS processes diverge at high energy, proportional up to s 2 /M 4 R where √ s is the c.m. energy of the V V → V V VBS process. In order to assess the validity range of the EFT description, we also calculate, in Section 3, the full one-loop corrections to the V V → V V scattering amplitudes which arise from the heavy new matter fields. Performing a partial wave analysis, their trajectory around the Argand circle is analyzed in Section 3.1 and we are thus able to specify for which isospin J R a perturbative description is still justified. At the limit of this perturbative range, we compare the full one-loop corrected cross sections for V V → V V scattering with the EFT approximation in Section 3.2. We find good agreement of the two approaches only for √ s 1.3M R , as is to be expected. While the above results are obtained in a simplified setting, neglecting the hypercharge coupling g and considering on-shell scattering only, an approximate implementation for the full qq → qqV V processes within the VBFNLO Monte Carlo program [20,21] is presented in Section 4. On the theoretical side, this allows for a more realistic comparison of the full model with its EFT approximation, at dimension-8 level, and with the unitarized replacement of the EFT in the T u model of Ref. [19]. Also here, good agreement between the three approaches is found for √ s < 1.3M R , very large deviations are seen in the threshold region for pair production of the heavy particles (M R < √ s < 2M R ), and qualitatively acceptable agreement of the T umodel with the full one-loop corrected results occurs well above pair production threshold. The VBFNLO implementation also allows to compare the expected VBS signals of high isospin fermion or scalar multiplets with LHC data. Finally, conclusions are presented in Section 5.

The toy model and its EFT
As a set of minimal, ultraviolet-complete models, we consider the addition of heavy fermionic or scalar SU(2) L multiplets to the SM. These new degrees of freedom are assumed to couple through the SU(2) L gauge interaction only, i.e. they are color singlets and have vanishing hypercharge. In addition, they are assumed to have negligible Yukawa or Φ 4 couplings to the SM Higgs doublet field and, thus, they do not mix with the other particles of the SM.

Definition and description of the toy model
To be specific, and in order to fix the notation, we consider Dirac fermions Ψ with mass M F or complex scalars Φ with mass M S , which transform under a generic SU(2) L representation R described by the generators t a R . The representation R is specified by its isospin, J R , i.e. each multiplet has (2J R + 1) components and t a R t a R = J R (J R + 1). Since the covariant derivatives of the new, hypercharge-zero fields contain SU(2) L gauge fields only, no relevant information will be lost in our loop calculations and EFT considerations by simplifying the model to a pure SU(2) L gauge theory, i.e. we restrict most of the discussion to the SU(2) L limit of the electroweak sector of the SM, where the hypercharge coupling g of U(1) Y is set to zero.
Within this approximation, the covariant derivative of the Higgs doublet field isD µ = ∂ µ +ig τ a 2 W a µ , with the Pauli matrices τ a , and the SU(2) L field-strength tensor is defined througĥ Suppressing the SM fermions and the Higgs potential, the relevant toy-model Lagrangian can be written as The first line corresponds to the well-known SM Lagrangian of the gauge and Higgs sector in the unitary gauge, albeit with equal W and Z masses, 2 m W = m Z = gv/2. The second line of Eq. (4) comprises the new particles and their interactions with the gauge fields, which give rise to familiar expressions for the Feynman rules. Our loop calculations will involve internal lines of only the new fermions or scalars. Thus, results in the unitary gauge, outlined above, are identical to those of a more general R ξ gauge. The new fermions' gauge couplings are assumed to be non-chiral and therefore their mass, M F , can be chosen arbitrarily (and large) prior to spontaneous symmetry breaking of the SU(2) L . In principle our toy model can host various fermionic and scalar multiplets without mass mixing. However, this will be irrelevant for our phenomenological discussion. Our toy model is closely related to new-physics solutions of problems that remain open in the SM, like an explanation of neutrino masses or dark matter: Large fermionic and scalar SU(2) L multiplets up to quintets that couple among each other and to the lepton sector of the SM are known to provide an explanation of neutrino masses [22][23][24][25][26][27][28]. However, we are interested in the impact of such SU(2) L multiplets on vector boson scattering only and thus allow the new particles to couple to the SM only through the gauge coupling of SU(2) L . Such large multiplets can also contain a dark matter candidate, as the lightest particle in the multiplet spectrum can be stable. Our model coincides with a class of minimal dark matter models [29][30][31], which range from SU(2) L triplets up to septets (see Ref. [32] for a review).
Embedding the heavy matter fields into the full SM, i.e. including the U(1) Y gauge field, the mass degeneracy of the additional SU(2) L multiplets at lowest-order in perturbation theory is lifted at the one-loop level, due to the mass splitting between photon, W and Z. As an example, and in agreement with the literature [29], for a J F = 4 fermion nonet of zero hypercharge, the differences in mass of the charged states with respect to the neutral state are given by ∆m 4± = m Ψ 4± − m Ψ 0 = 2.7 GeV, ∆m 3± = 1.5 GeV, ∆m 2± = 0.66 GeV and ∆m ± = 0.166 GeV. Thus, cascade decays from the heavier, charged states to the neutral, stable state proceed through far off-shell gauge bosons that result in very soft leptons and pions. If on the one hand the decay length of the charged states is short enough that they decay within the interaction 2 The absence of a massless photon in our toy model and the common W and Z masses significantly simplify the analytical results for VBS amplitudes as well as the unitarity considerations in Section 3.1.
region and do not result in a charged track in the detector, but on the other hand the final state leptons and pions are still soft enough, then such dark matter models are very hard to be probed at the LHC.
A thorough discussion of the current exclusions of a minimal dark matter quintet, based on LHC searches, can be found in Ref. [33], including projections for the high-luminosity runs. Expected exclusion bounds range to quintet masses up to 750 GeV for integrated luminosities of 3 ab −1 . Going beyond quintets, i.e. beyond J R = 2, results in the existence of additional charged states, that are pair-produced with large cross sections that rise proportional to J R (J R +1) [29]. It is also known that the minimal dark matter model needs masses of the dark matter candidate in the range of 10 TeV [34] in order to achieve the correct relic abundance, such that for lighter masses, as discussed in this paper, only a fraction of the observed abundance can be explained by the neutral, stable particle of the minimal dark matter model. Alternatively, extensions with additional singlet states [35] are discussed in the literature, which, however, change the collider phenomenology substantially without impacting too much the discussion carried out in this paper. Because the collider constraints can always be avoided by somewhat increasing the mass parameters M F or M S , we omit a more detailed discussion of direct exclusion limits from existing collider searches, as this is not the focus of the present paper.

One-loop corrections to gauge-boson vertices, isospin factors and renormalization
The additional SU(2) L multiplets contribute to various observables at the one-loop level. Here we focus on contributions to VBS, which enter through propagator, three-and four-boson vertex corrections. In practice we use FeynCalc [36,37] for the handling of Dirac structures, we regularize ultraviolet divergences through dimensional regularization and express our findings in terms of Passarino-Veltman integrals [38,39] to obtain numerical results. We refrain from showing the explicit, rather lengthy expressions for the corrections to the three-and fourgauge-boson vertices. 3 Instead, we here concentrate on the general structure of the one-loop contributions, including their representation dependence. We start with the vacuum polarization. For the fermion case only a single diagram contributes. For the scalar case a second diagram involving the four-particle vertex needs to be added. As the propagator corrections are purely transverse, they generically take the form for the scalar case and similarly, replacing T S Π S (p 2 , M 2 S ) by T F Π F (p 2 , M 2 F ) for the fermion loop. For the diagrams involving two three-particle vertices the isospin factor is Tr t a R t b R = T R δ ab while the scalar diagram with the four-particle vertex has Tr {t a R , t b R } = 2T R δ ab . Therein T R = T F or T S is the index of the representation, given by T R = 1 3 [J R (J R + 1)(2J R + 1)] for an isospin J R representation. The explicit form of the propagator corrections is given by These expressions can be expanded in small momenta p 2 motivated by p 2 ∼ m 2 W M 2 F , M 2 S . As we employ dimensional regularization, working in d = 4 − 2 dimensions, we parameterize the divergences in terms of ∆ = 1 − γ E + log(4π). The propagator corrections then turn into where µ denotes the renormalization scale. It is therefore clear that choosing an MS scheme with µ = M F , M S leads practically to an on-shell scheme with almost vanishing propagator corrections, such that electroweak precision data are not significantly impacted by the extra degrees of freedom. We now turn our attention to the three-boson vertex. The result is generically of the form for heavy fermions, and analogously for scalars. Here "permutations" of the external gauge bosons refers to the Feynman graph with the opposite direction of the fermion arrow. For the three-boson vertex correction in the scalar case, all diagrams involving the W W ΦΦ seagull vertex vanish individually because, among other reasons, their isospin factor Tr t a R {t b R , t c R } is zero. The sum of the depicted diagrams, on the other hand, carries an isospin factor Tr t a For the scalar case, the one-loop corrections to the four-gauge boson vertex take the form while the corresponding heavy fermion contribution,Γ µναβ,abcd 4,F (p 1 , p 2 , p 3 , p 4 , M 2 F , J F ) is given by a box-diagram analogous to the one shown, and its permutations. Above, the three permutations for the scalar one-loop graphs involving two four-particle vertices are labeled s-, t-and u-type diagram. For all diagrams we get traces over four generators, which generically result in Here, ] is the index, and C 2,R = J R (J R +1) is the quadratic casimir of a representation. Eq. (10) separates the isospin symmetric part proportional to T R (3C 2,R −1), which is the only contribution to pure ZZ → ZZ scattering, from terms coming from products of tensors. Alternatively, all one-loop corrections can be sorted into terms proportional to T R only and terms proportional to C 2,R T R , which provides the full isospin dependence. As a consequence, numerical evaluations can be performed for a fixed mass choice and can then easily be reweighted to any representation R. We use this procedure in VeBoS and our implementation in VBFNLO.
We now turn our attention to the renormalization of ultraviolet divergences. Since we are only considering n-point functions of the SU(2) L gauge fields which are transverse, as in Eq. (5), only the Tr Ŵ µνŴ µν term in Eq. (4) can produce the necessary counter terms, via the renormalization of the gauge fields and the coupling constant. Denoting the renormalization constants between the bare Lagrangian parameters (with subscript 0) and the physical parameters as the renormalization constants of the three and four gauge-boson vertices are given by Z 3W = Z g Z 3/2 3 and Z 4W = Z 2 g Z 2 3 , respectively. Writing Z i = 1 + δ i , we obtain δ 3W = δ g + 3 2 δ 3 and δ 4W = 2δ g + 2δ 3 , which defines the counter-terms for the three and four gauge-boson vertices, respectively. The counter-term of the propagator reads δ 2W = δ 3 . It directly follows from the vacuum polarization in Eq. (7) that in the MS scheme the field renormalization constant δ 3 is given by assuming n F and n S copies of a fermionic and scalar multiplet with isospin J F and J S , respectively. From the ultraviolet divergence of the three gauge-boson vertex we obtain and, therefore, the renormalization constant of the SU(2) L coupling is given by δ g = − 1 2 δ 3 . It is easy to check that indeed the now fully determined counter term to the four gauge-boson vertex, namely δ 4W = 2(δ g + δ 3 ) = δ 3 , also cancels the divergence of the corresponding one-loop corrections to gauge-boson four-point functions.

Effective-field theory description
Considering the low-energy expansion of the now finite two-, three-and four-point functions of the gauge bosons, we can derive the contributions of additional heavy SU(2) L multiplets to the EFT Lagrangian i.e. we can now determine the Wilson coefficients f i /Λ 2 and f i /Λ 4 of the dimension-6 and dimension-8 operators which are relevant for VBS in our concrete model. By construction, the SU(2) L gauge bosons are the only SM particles which couple via the heavy fermions or scalars and which can, at the one-loop level, appear as fields in the effective Lagrangian. Gauge invariance then implies that we can restrict the EFT Lagrangian to operators which contain the SU(2) L field strength tensorŴ µν of Eq. (3) or its commutator with covariant derivatives. For a description of vector boson scattering we need vacuum polarization effects (dimension-4 operators and higher), anomalous triple-gauge couplings (dimension-6 and -8 operators) and anomalous quartic-gauge couplings (dimension-8 operators).
The dimension-4 operator is unique and given by It is already part of the SM kinetic term. Since it receives contributions from the additional SU(2) L multiplets it is explicitly listed here. However, as discussed in the last section, it disappears as a result of renormalization. As a basis for the dimension-6 operators, following Ref. [5], we choose Our focus is on the dimension-8 operators which first appear as aQGC in vector boson scattering, and which involve four field-strength tensors, namely Here we follow the notation of Refs. [7,19,41,42]. However, since our definition ofŴ µν includes the coupling factor ig, the operators in Eq. (16) contain an overall factor g 4 as compared to thé Eboli convention of Refs. [7,41]. Please note that O T 3 is a non-vanishing, linearly independent operator which was missing in previous discussions and, therefore, is also not experimentally constrained explicitly. We refer to the four operators in Eq. (16) as T -operators. 4 For three field-strength tensors in combination with covariant derivatives, seven operators can be written down which are, however, related through integration-by-parts for the covariant derivative, the Jacobi identity for covariant derivatives, andŴ µν = D µ ,D ν . Our set of independent operators is then given by For the case of two field-strength tensors in combination with covariant derivatives we pick as, again, other combinations are not linearly independent. The above choice of basis is by no means unique but simply driven by the physics of the underlying theory. Looking at the dimension-6 operators one could make an argument for using the Warsaw Basis [6] as this could make comparisons easier and gives more direct access to exclusion limits for the Wilson coefficients. Rewriting the O DW operator using integration-byparts, the commutator for the covariant derivatives as well as the Jacobi identity for covariant derivatives, one arrives at 5 4 Similarly, the operators named O T5 , O T6 and O T7 in Ref. [7] should be supplemented by O T4 = Tr Ŵ µνŴ αβ B ναBβµ when discussing mixed SU(2) L and U(1) Y gauge boson scattering. 5 The operator O W W W is equivalent to the corresponding operator in the Warsaw Basis up to a factor, At this point we can use the equations of motion for the W -field twice: Here the last term accounts for corrections to the equations of motion due to the presence of O DW and O W W W which are of order 1 Λ 2 . This leaves us with four types of terms: operators containing four Higgs fields, the ones containing fermionic currents, and also terms of order 1 Λ 2 and 1 Λ 4 of which only the latter may be neglected as they are of order 1 Λ 6 when including the Wilson coefficient of O DW . The terms with four Higgs fields can again be rewritten by using relations for the SU(2) L generators as well as the equations of motion for the Higgs field and one finds where the operators on the r.h.s., with φ representing the SM Higgs doublet field, are written in the Warsaw basis. Here ψ stands for either leptons or quarks. The last two operators are just symbolizing the class of additional operators that can appear following the conventions in Ref. [6]. Individual terms in this new set now give local corrections to Higgs and fermionic observables that, at the one-loop level, are not present for the BSM heavy matter fields which we consider. Therefore, these corrections have to cancel each other when calculating amplitudes for physical processes in our new physics model. This also implies that any exclusion limits for Wilson coefficients are only useful when provided with full error correlations. With this in mind it becomes apparent that there is no advantage in using the Warsaw basis over the basis defined by Eqs. (15)-(18).

Wilson coefficients due to additional SU(2) L multiplets
After presenting the full EFT Lagrangian we can map our concrete model onto the Wilson coefficients at low energies by integrating out the heavy new degrees of freedom. Again we assume to have n F and n S copies of fermionic and scalar multiplets with isospin J F and J S , respectively. We expand the one-loop corrections in small momenta with respect to the mass scale of the new heavy degrees of freedom, i.e.
1, which for the Wilson coefficients originating from propagator corrections results in We continue to expand the three gauge boson-vertex corrections in a similar manner in small momenta, which leads to Finally, the Wilson coefficients of the T -operators are given by As discussed previously, all UV divergences of propagator-, triangle-, and box-diagrams are contained in the coefficient of the field-strength-square operator O W W and, hence, disappear upon renormalization of the SU(2) L gauge coupling g and the gauge boson fields. Some additional remarks are in order: • The isospin factors of propagator and three-boson-vertex corrections are given by T R and thus rise with J 3 R . The Wilson coefficients of the T -operators, on the other hand, are enhanced by an additional casimir factor, i.e. they rise as J 5 R , enhancing their importance over the dimension-6 operators for high isospin.
• Individual operators in the scalar case are more suppressed than for heavy fermions of the same mass and isospin.
• The Wilson coefficients of the T -operators f T 0 and f T 1 receive opposite sign contributions from heavy scalars vs. heavy fermions, which can induce cancellations between T -operator contributions to e.g. VBS cross sections. f T 2 changes its sign at J R ∼ 1 for both the fermionic and scalar multiplets.

Confronting Wilson coefficients with experiment
Having obtained the low-energy EFT of the general model, we next investigate whether specific model realizations are compatible with experimental limits on four-fermion contact terms, anomalous triple-and anomalous quartic-gauge couplings. Since we are only interested in rough estimates for our toy model, we restrict the discussion to the single multiplet realization and determine bounds in the (J R , M R ) parameter space. The insertion of O DW into the gauge boson propagator of a qq → l + l − amplitude scales like at high energy and, thus, reduces to a flavor universal, dimension-6 four-fermion interaction which affects left-chiral amplitudes only. Writing the coefficient of theQ L γ µ Q LL γ µ L EFT operator as 2π η LL /Λ 2 LL one finds η LL = −1 and for a single fermion or scalar multiplet in the loop, respectively. In the EFT approximation, experimental bounds on four-fermion contact terms, such as recent ATLAS measurements [43], can thus be turned into the lower limits on M 2 R /T R which are shown in Fig. 1. aTGCs are typically measured in the LEP parameterization, where the effective W W Vcoupling parameters ∆g 1 , ∆κ and λ are used [44]. λ parameterizes transverse weak boson interactions and thus is the one of primary interest here. The loop-corrections of our model contribute both to propagators, via the dimension-6 operator O DW , and to triple gauge-boson vertices, via both O DW and O W W W . In theqq → W W amplitude, in the SU(2) L limit, the former is equivalent to an energy dependent ∆g 1 = ∆κ contribution, while the latter are best captured after the operator redefinition in Eq. (19). One finds Going beyond the SU(2) L limit, the λ-aTGC can easily be compared to existing data, via λ γ = λ Z = λ. The ∆g 1 and ∆κ contributions would appear as process dependent form-factors, however, which are constrained at a similar level as the induced λ-aTGC. Table 1 summarizes relevant bounds on four-fermion couplings, aTGCs, and aQGCs. For the bounds on T -operators, we use the recent CMS results of Refs. [46,47] on leptonic decay modes in VBS processes, which also take into account limitations from unitarity. Since these aQGC bounds are only specified for insertions of single EFT operators, Fig. 1 depicts the limits on the (J R , M R ) parameter space for each Wilson coefficient separately. Although the strong correlations between the Wilson coefficients in our model are not taken into account by this Coefficient limit without unitar. channel limit with unitarity cut channel  [43], aTGCs [45] and aQGCs [46,47]. The values for aQGCs from the CMS analyses are derived from leptonic decay channels in VBS processes only. aQGC bounds are presented with and without a unitarization in form of a cut at the perturbative unitarity limit of the EFT calculation. In addition, we name the channels from which the presented bounds are derived.
procedure, it suffices here, since we are only aiming at a rough estimate of allowed regions in the (J R , M R ) plane. 6 The figures indicate that a measurement of aQGCs has higher potential to find effects of the BSM physics considered here than aTGC measurements. To be more concrete, we consider the cases of a fermion with J F = 4 and of a scalar with J S = 6 in the following, values which represent maximal multiplet sizes compatible with unitarity and perturbativity limits to be discussed in Section 3.1.
For the fermion multiplet, the strongest bound is associated with searches for qqll contact terms, disfavoring masses below 700 GeV. Somewhat weaker, at M F > 540 GeV, are bounds derived from limits on the dimension-8 T 1 -operator, which do take into account unitarity constraints. Restrictions via the other T -operators are qualitatively similar. The aQGC bounds in turn are much more stringent than the bound from the λ-aTGC, which only excludes masses below 100 GeV in the EFT dimension-6 approximation.
The discussion is comparable in the scalar case. The bounds in the (J S , M S ) parameter plane are less stringent, however, as anticipated from the higher numerical suppression in the Wilson coefficients. Again all T -operator bounds appear to be stronger than the aTGC bounds, even when including unitarity considerations only for the former (see Fig. 1, right), making VBS more promising for the investigation of this model than vector boson pair production. For J S = 6, the bound on the coefficient of the T 1 operator disfavors masses below 470 GeV, while searches for qqll contact terms appear to imply a slightly weaker bound of M S > 430 GeV.
Guided by these constraints, we will use J F = 4, M F = 600 GeV and J S = 6, M S = 600 GeV as benchmark points in the remainder of the paper to illustrate consequences for VBS at the LHC. For fermion multiplets, M F = 600 GeV appears to be in slight tension with qqll contact term limits. However, additional BSM effects which we do not consider, like an extra Z , might ameliorate the tension. Since we are not explicitly proposing a viable BSM model, but are rather interested in generic features for VBS, we prefer a common mass for fermions and scalars in our subsequent discussions.
Bounds on (J S , M S ) with cut at unitarity limit for aQGCs Figure 1: Translation of the experimental bounds on Wilson coefficients as in Table 1 to the (J R , M R ) parameter space: masses above the curves are allowed by the individual measurements. Shown are the cases of a fermion multiplet (left panel) and a scalar multiplet (right panel). For aQGCs, experimental results are derived with a step-function form-factor, cutting off the EFT prediction at the unitarity limit (see Ref. [3] for details).
Another reason for taking the (J R , M R ) bounds derived from EFT considerations only as indicative is the concern that the experimental EFT analyses, especially for the qqll contact terms and the aQGC, are largely based on data with dilepton or V V invariant masses above 1 TeV. For our large multiplet model the EFT approximation cannot be expected to hold for invariant masses of the order of the threshold energy of 2 M R or larger.
As a case in point, let us have a brief look at qq → l + l − production at the LHC, from which qqll contact term bounds are derived. In this process and at 1-loop level, our model only influences the vector boson propagator. The Dyson-resummation of the on-shell renormalized propagator corrections, which follow from Eqs. (5) and (6), leads to for W ± exchange between purely left-chiral fermions. For comparison with l + l − data, γ-Z propagator mixing needs to be considered and the BSM contributions are somewhat diluted by U(1) Y hypercharge contributions which, in our toy model, are left unchanged with respect to the SM. The results are shown in Fig. 2 where the ratio of the dilepton invariant mass distribution in our toy model as compared to the SM is displayed for M R = 600 GeV and a selection of isospins. Also shown, for the J F = 4 and J S = 6 benchmark points, is the BSM to SM ratio in the EFT approximation, considering the equivalent dimension-6 qqll contact term only. While BSM effects are negligible in the m ll < m R /2 region (below 1.5% for J F = 4 and well below 1% for J S = 6), sizable deviations from the SM are expected at higher energies, for large isospin representations. The EFT dimension-6 approximation reproduces the full oneloop result considerably below production threshold, but it fails spectacularly above threshold. Since the bounds on qqll contact terms in Table 1 heavily depend on data at m ll > 2.5 TeV, one concludes that the contact term bounds on the (J R , M R ) parameter plane cannot be taken at face value. However, around the production threshold, the deviations from the SM are even bigger than suggested by the EFT approximation, and a full analysis of the model should be performed in this region to extract bounds on M R .
As a first rough approximation, we have translated the results of a recent CMS analysis of dilepton production [48] to 1-and 2-σ error bands in Fig. 2, assuming fully correlated systematic errors between electron and muon channels and adding systematic and statistical errors in quadrature, otherwise. Inspection of Fig. 2 then suggests that for M R = 600 GeV and at 95% CL the scalar benchmark point with J S = 6 is still allowed, while strong evidence for the fermionic case with J F = 4 should have been seen already. Since we consider our setup as a toy model, and since we are primarily interested in its qualitative predictions for VBS, we postpone a more detailed analysis to future work and use the benchmark point merely for illustration in the following.

On-shell vector boson scattering at one-loop
The large deviations between EFT results and the full model prediction for lepton pairproduction provide a strong motivation to also consider the full model for VBS predictions. For a study of its main phenomenological features, we start with a discussion of on-shell VBS. In the following we neglect the normal SM electroweak corrections, which are parametrically small compared to the contributions from extra SU(2) L multiplets, due to the J 3 R to J 5 R enhancement of the latter. Corrections from the new heavy degrees of freedom originate from the previously discussed one-loop contributions to the gauge-boson propagators and the three-and four-gauge-boson vertex functions. In terms of Feynman diagrams, for external gauge bosons in the adjoint basis (W µ a i with a i = 1, 2, 3) the scattering amplitudes can be depicted as The last equality follows from isospin symmetry, because all isospin combinations for the four indices a i associated with the external gauge bosons W µ a i can be expressed by Kronecker deltas in the adjoint space. Obviously, the reduced amplitudes M 1 , M 2 , and M 3 can be obtained from each other by crossing.
In this section we consider on-shell vector boson scattering in the SU(2) L limit, i.e. we set g = 0, which eliminates any photon contributions 7 and leads to Z = W 3 and m Z = m W = gv 2 , thus simplifying the scattering kinematics. Projecting the amplitude of Eq. (36) onto Z bosons (with a factor δ 3a i ) or W bosons (with a factor 1 7 As a result one does not need to contend with any infrared or collinear singularities due to the massless photon. Also the unphysical Rutherford singularity (due to t-channel photon exchange) in charged W -scattering is eliminated, which in real life is regularized by the space-like virtuality of the incoming weak bosons.
As indicated by the absence of double vertex insertions in Eq. (36), our full amplitudes correspond to a coupling expansion to order g 4 , with MS renormalization. Before proceeding, we need to address whether such a perturbative expansion is justified for large multiplicities, 2J R + 1, which we do by a unitarity analysis of partial wave amplitudes.

Partial wave analysis of VBS amplitudes
The notation for VBS partial wave amplitudes follows Ref. [19], the results of which are implemented in VBFNLO. However, the present case is simplified by taking the SU(2) L limits and by considering on-shell scattering only. Exploiting angular momentum conservation, the helicity amplitudes M λ 1 λ 2 λ 3 λ 4 of the various V V → V V processes are decomposed into coefficients of Wigner d-functions where λ ij := λ i − λ j , and the normalization factor is given by W /s and statistics factors S f or S i = 1/2 in case of two identical weak bosons. N f i is chosen such as to simplify the constraints on the partial wave amplitudes A j λ 1 λ 2 λ 3 λ 4 (s) which follow from S-matrix unitarity, namely Here Im A j λ 1 λ 2 λ 3 λ 4 corresponds to the anti-hermitian part of the VBS scattering amplitude, understood as a matrix in isospin and helicity space. The sum over n on the right hand side includes all relevant di-boson intermediate states. Equality in Eq. (39) would only be reached when summing over a complete set of intermediate states, including multi-particle states. For any normalized linear combination of diboson states, |n , the expectation value of A j will lie within the Argand circle. In particular, this will be true for the eigenvector |n max of the largest eigenvalue of Re (A j ) or Im (A j ). Setting which corresponds to the largest eigenvalue for a normal scattering matrix, we thus get the Argand circle condition Re a j 2 + Im a j − 1 2 ≤ 1 .
Since a finite perturbative expansion is only expected to fulfill this requirement approximately, e.g. an on-shell tree-level calculation remains always on the real axis and is thus never compatible with the strict requirement, the perturbative unitarity bounds are given by which should be understood as demanding that the eigenvalues of A j remain close to the Argand circle in a reasonable perturbative calculation. When investigating the behavior of the partial wave amplitude contributions from single (J R , M R ) multiplets, one finds sizable effects only for the transverse helicity amplitudes M 1111 , M 11-1-1 , M 1-11-1 , M 1-1-11 (and with corresponding parity flipped helicities). Thus it is sufficient to diagonalize the partial wave amplitudes in this restricted helicity space 8 . This is made easier by the fact that the VBS amplitudes are also block-diagonal in the basis of total isospin of the two particle states, |J 2 , J 2,z , as we work in the g = 0 limit of the SM and our BSM contributions do not break the SU(2) L symmetry. The dominant eigenvalue is found in the j = 0 partial wave and corresponds to vanishing total isospin of the weak boson pair, J 2 = 0. Its BSM contribution is called a 0 in the following.
Results are presented in Fig. 3 for both the fermionic multiplet (left panel) and the scalar case (right panel), for a range of isospin choices, J R , for a multiplet mass of M F = 600 GeV = M S , and for 200 GeV < √ s < 6000 GeV. For J F ≤ 4 and J S ≤ 5, the path of a 0 J R (s) in the complex plane is compatible with the unitarity bounds of Eq. (42), while the scalar case with J S = 6 can be considered marginally perturbative. The strong dependence on multiplet size was to be expected because the four-boson vertex function at 1-loop order grows like J 5 R . In the following we would like to estimate maximal deviations from the SM which can be expected in VBS. For this purpose, the multiplet realizations of (J F = 4, M F = 600 GeV) and (J S = 6, M S = 600 GeV), chosen earlier, appear well suited for a qualitative discussion. Fig. 4 shows the dependence of both the real and imaginary parts of the dominant partial wave eigenvalue, a 0 J R (s) on the center of mass energy of the V V → V V process. Here we fix J F = 4 for the fermion model (left panel) and J S = 6 for the extra scalar multiplet (right panel). Also the EFT prediction for a 0 J R (s) is shown, which is purely real and which, sufficiently below production threshold at √ s = 2 M R , agrees well with the full calculation. The LO SM contribution to the dominant partial wave amplitude is also depicted (gray solid line), for comparison, and because it should be added to the anomalous contribution, a 0 J R (s). Adding these two contributions, the resulting full partial wave amplitude stays fully within the perturbative unitarity bounds of Eq. (42) for the J F = 4 fermion multiplet and agreement for the J S = 6 scalar case is also improved at high energy. Violation of the |Re (a j ) | 1 bound for the scalar case is actually limited to a small region around threshold, while the imaginary part of the amplitude, which dominates somewhat above threshold, is consistent with the unitarity limit. The problem with the real part of the amplitude around √ s = 2 M S could be ameliorated by SU(2) L -breaking effects, i.e. by mass splitting the multiplet and thus distributing its threshold effects over a larger energy range, by binding effects due to additional (strong) interactions of the scalar multiplet(s) or by other modifications of our toy model. This provides another argument to not discard the J S = 6 case.
In order to demonstrate the effect of different mass values, Fig. 4 shows results for M R = 600 GeV and M R = 1000 GeV. One finds excellent scaling: the amplitude effectively only depends on s/M 2 R because the electroweak scale merely enters via the mass, m W , of the external gauge bosons, and this gives tiny (m W /M R ) 2 corrections. This scaling behavior implies that the effects which we will demonstrate with M R = 600 GeV in the following can simply be shifted to higher energy for larger multiplet masses, with an approximately invariant ratio of BSM to SM VBS cross sections.

Cross sections for on-shell vector boson scattering
With limited statistics and sizable backgrounds for VBS events at the LHC, integrated, unpolarized VBS cross sections are sufficient for a first survey of LHC capabilities. According to Eq. (37), the information is contained in three crossing related helicity amplitudes which can be disentangled by measuring cross sections for different combinations of weak bosons: samesign W W scattering, i.e. W ± W ± → W ± W ± gives access to M 2 + M 3 , W Z → W Z depends on M 2 only, and the ZZ production processes W W → ZZ and ZZ → ZZ provide separate information on M 1 .
Separating SM and new physics (NP) contributions, the squared amplitude, which determines the cross section, is given by Our new physics amplitude is the 1-loop contribution due to the additional heavy fermion or scalar multiplet. In a typical NLO calculation, only the interference term of the NP would be considered, together with the SM Born contribution. On the other hand, when considering anomalous couplings, frequently also the |M NP | 2 term is included in the cross section prediction, and this is indeed the default VBFNLO implementation of VBS. Since there are no infrared divergences in the present case, and because the purely anomalous part is necessary for a prediction of the ZZ → ZZ cross section, which due to the small Higgs boson mass has a tiny SM amplitude, we do include the |M NP | 2 term in our calculation, but we also show the NLO-type result in the following, dubbed "full model (int)" in the figures. An analysis of the parametrically leading terms in a perturbative expansion support the inclusion of the quadratic NP term. At large J R , the dominant two-loop contribution is given by diagrams with an additional virtual vector boson propagator within the the multiplet loop of the four-particle vertex functions, e.g. and in the fermion case, which have representation factors of T R C 2 2,R ∼ J 7 R . On the other hand, the square of the one-loop four-particle vertex function, i.e. the |M NP | 2 term above, is enhanced by a factor T 2 R C 2 2,R ∼ J 10 R (see Eq. (10)), which would not arise until the three-loop level in the interference term.
Subsequently we present the integrated, unpolarized cross sections for scattering angles 5 • < θ < 175 • in the center of mass frame for the still academic W ± W ± → W ± W ± , W Z → W Z, W W → ZZ, and ZZ → ZZ scattering processes. A more complete LHC simulation of the full qq → qqV V processes will be postponed until Section 4. The limited angular range reduces the forward and backward enhancement in the SM contribution, which is due to t-or uchannel Higgs or W, Z exchange. For the new-physics contribution, we use the single multiplet realization with parameter choices discussed in Section 2.4.
Cross section results, as a function of the center-of-mass energy √ s = m V V , are shown in Fig. 5 for the case of a fermion nonet and in Fig. 6 for a single scalar multiplet with J S = 6. A striking and well-known feature of the SM prediction is the strong suppression of the ZZ → ZZ cross section, due to the small value of m H = 125 GeV. According to Eq. (37) this implies a relation which does not carry over to the NP contributions. Since M 1 describes W + W − → ZZ while M 2 + M 3 is the amplitude for same-sign W W scattering, one finds opposite signs for the interference terms, M SM M NP * , when comparing these two processes. M 2 is the amplitude for W Z scattering, which is related to M 3 by θ → π − θ crossing. As a result one finds the same interference characteristics for W Z scattering as for same-sign W W scattering. Effects are larger for the latter, however, by roughly a factor of two, because the important NP amplitudes are fairly independent of scattering angle and thus add up, while the SM t-and u-channel amplitudes M 2;SM and M 3;SM peak in opposite hemispheres.
Comparing the J F = 4 fermion case in Fig. 5 with the J S = 6 scalar multiplet in Fig. 6, the effect of the scalar multiplet on cross-sections is more pronounced due to the higher representation. The constructive interference between the SM and the anomalous contribution in the W W → ZZ process (lower left) is particularly striking around the threshold for fermion or scalar pair production, which results in a peak structure centered around √ s = 2M R . For same-sign W W and for W Z scattering (upper row) we see destructive interference between the SM and the anomalous contribution in the low-energy regime. This leads to a dip in the cross section at pair production threshold which, however, is soon eclipsed by the strong cross section enhancement above threshold, in particular for the scalar case. This peak above threshold is driven to a substantial degree by a large contribution from the imaginary part of the scattering amplitude, which was already evident in Fig. 4.
Comparing the behavior of the total amplitude squared with the NLO approximation, we observe discrepancies in the high energy regime, starting at threshold, especially for the cases of W Z and same-sign W W scattering. This demonstrates that our model realizations, with large isospin representations, are only marginally perturbative and might exhibit sizable higher order corrections. This is also evident for ZZ scattering (lower right panels) which is clearly dominated by the anomalous contribution, i.e. the one-loop amplitudes effectively provide LO estimates, similar to gg → ZZ in the SM.
In all cases one finds good accordance between the full model and its complete EFT realization in the low-energy regime. For the EFT, four different results are shown. The complete "EFT" curves (orange dotted) contain the contributions from all dimension-6 and dimension-8 operators in Section 2.4. They are virtually indistinguishable from the orange dashed lines,     For a better evaluation of the validity of the EFT and its subsets in the low-energy region, we present in Fig. 7 and Fig. 8 the cross sections, normalized to the SM prediction, for same-sign W W scattering and the W W → ZZ process. As discussed above, results for W Z scattering are qualitatively the same as for same-sign W W scattering.
The EFT agrees well with the full model for energies below ∼ 800 GeV, which is consistent with its naive validity range, given by the mass scale of the loop particle. In the fermion case, the discrepancy stays below ∼ 1 % up to 1000 GeV, but would grow with J 5 F . We also note that the EFT dimension-6 contribution only accounts for a slowly varying offset of less than 1 % to the pure SM cross section (gray dash-dotted line).

Implication for V V jj events at the LHC
The analysis of the last section provided cross sections for on-shell 2 → 2 VBS processes, in the g = 0 limit, and thus was largely academic. A more realistic simulation for the LHC needs to consider the full qQ → qQV V processes, including decay of the produced electroweak gauge bosons, it needs to consider the off-shell nature of both the space-like initial and the time-like final gauge bosons, and it needs to go beyond the SU(2) L limit and incorporate that the W 3 is actually a linear combination of photon and Z as mass eigenstates. While a full implementation of the model into VBFNLO is foreseen for the future, here we only make a first, approximate assessment of consequences for the LHC.
A prominent feature of typical VBS kinematics consists of small virtualities for all four electroweak bosons in the V 1 V 2 → V 3 V 4 subgraphs, |p 2 i | s = m 2 V V . At such small p 2 i , modifications of the weak boson propagators are tiny, as evidenced by Fig. 2. Also corrections to W W Z or W W γ vertices, which appear for a high virtuality gauge boson which is radiated off a quark line in W Zjj and W + W − jj production, are sub-dominant as was seen in the discussion of aTGCs in Section 2.5. The remaining BSM effects in our model then contain V 1 V 2 → V 3 V 4 subgraphs with off-shell helicity amplitudes M V BS λ 1 λ 2 λ 3 λ 4 ({p 2 i }; s, t, u). Following Ref. [19], the BSM part, which needs to be added to the SM qq → qqf 1 f 2f3 f 4 amplitude is then given by Here the J (p i , λ i ) are off-shell polarization vectors (as defined in Ref. [19]) and the J q→qV i and J V i →f k f l are quark-and lepton currents. A final approximation now consists in replacing the off-shell helicity amplitudes M V BS λ 1 λ 2 λ 3 λ 4 ({p 2 i }; s, t, u) = M V BS λ 1 λ 2 λ 3 λ 4 ({p 2 i }; s, θ), determined in the VBS center-of-mass frame, by their on-shell values M V BS λ 1 λ 2 λ 3 λ 4 ({m 2 W }; s, θ) = M λ 1 λ 2 λ 3 λ 4 (s, θ) of Eq. (38). We have verified that the above approximations are good at the 10% level by comparing on-and off-shell VBFNLO results for the dimension-8 T -operators [40].
In the next step we use our new VBFNLO implementation of full on-shell model amplitudes to study fiducial cross sections for the various VBS processes. We want to find out to which degree ATLAS and CMS strategies for measuring anomalous quartic gauge couplings 9 would also reveal high multiplicity extra matter multiplets, via their impact on weak boson 4-point functions. All results are produced for a proton-proton collider with a center-of-mass energy of 13 TeV using the MMHT2014 [50] parton distribution functions at leading order, which are linked through LHAPDF [51]. The Higgs-boson mass is set to 125.09 GeV. For the W and Z boson mass VBFNLO uses M W = 80.398 GeV and M Z = 91.1876 GeV, respectively. The weak mixing angle is set to sin 2 θ W = 0.222646. However we remind the reader that the inserted on-shell amplitudes are derived in the SU(2) L limit, in which M W = M Z = 81.18 GeV. When we depict the number of events we assume an integrated luminosity of 137 fb −1 .

Destructive interference in same-sign W W and W Z scattering
We use VBFNLO to calculate the cross sections for the electroweak processes with l, l ∈ {e, µ}, which were exploited in Refs. [3,52] to set bounds on anomalous quartic gauge couplings. In order to make a qualitative comparison with the results of Ref. [3], in particular their Fig. 6, we have implemented a very similar cut-flow. For the W Z final state this cut-flow includes: Therein the abbreviation "3l" denotes the three-lepton system, i.e. the sum of the three lepton 4-momenta, and its invariant mass m 3l . In contrast to the CMS analysis we only simulate the flavor combination l = e, l = µ and as an estimate for all flavor combinations multiply our results by a factor of 4. As a consequence we have just one cut on the transverse momentum p l T , which is chosen equal for all leptons. The cuts depicted in the last two lines of Eq. (47) are typical vector boson scattering cuts, which enhance the contribution of electroweak VBS over QCD-induced V V jj events and other SM backgrounds.
The cut-flow for the same-sign W W final state is very similar: for charged leptons we merely replace the m 3l cut by m ll > 20 GeV. Furthermore we set max(z * l ) < 0.75. In contrast Events/bin pp → W ± Zjj → lνl l jj, l, l ∈ {e, µ} J S = 6, M S = 600 GeV Events/bin pp → W ± Zjj → lνl l jj, l, l ∈ {e, µ} J F = 4, M F = 600 GeV Figure 10: Same as Fig. 9, but for a fermionic multiplet with J F = 4 and M F = 600 GeV.
to Ref. [3] we do not have a cut on |m ee − m Z |, as we generate our result for the flavor combination l = e, l = µ and multiply by 2.
We present our results primarily for the scalar case with the combination J S = 6 and M S = 600 GeV which, according to the discussion in the previous sections, may be considered a maximal plausible signal within our framework. The results for the fermionic case with J F = 4 and M F = 600 GeV, which is disfavored by Drell-Yan data, do not provide much additional insight and differ visibly only in the vicinity of the threshold peak, which is more pronounced in the scalar case. We present the number of events for an integrated luminosity of 137 fb −1 as a function of the invariant mass of the final state W W and W Z system which, however, is not accessible experimentally. Therefore, also for a direct comparison with Fig. 6 of Ref. [3], we add figures showing the transverse mass as defined there.
For the W ± Z final state the results for the scalar case are depicted in Fig. 9 and for the fermionic case in Fig. 10. The m(W ± Z) invariant mass distribution in Fig. 9 (upper left) and its ratio to the SM contribution (upper right) show destructive interference between the multiplet contribution and the SM well below the threshold at m(W ± Z) ≈ 2M F = 1200 GeV. The EFT description follows this behavior until around m(W ± Z) ≈ 1000 GeV, after which it clearly deviates from the full-model prediction. In this region of reliable EFT description the BSM effects are tiny, however, deviating from the SM by at most 1%, which renders their observation hopeless. The figures also depict the contribution from the combined dimension-8 T i operators only, which well coincides with the full EFT description. This exacerbates the observation in Section 3.2 that the contribution of the dimension 6 operators is sub-dominant and irrelevant for VBS in practice. The lower panels of Fig. 9 show the transverse momentum distribution, in which the destructive interference is not visible anymore, due to the migration of excess events at higher invariant mass into the low transverse mass region. The lower right panel reproduces the binning of the events in Fig. 6 of Ref. [3]. Comparison of the last two bins with the data (which are compatible with the SM) reveals that the full fermion model in Fig. 10 cannot be excluded, whereas scalar model at J S = 6 in Fig. 9 should have been seen with considerable significance. Similarly, the sizable signals of the non-unitarized EFT's in the last (overflow) bin are excluded. Also shown are unitarized versions of the EFT description, following the unitarization procedure of Ref. [19], which stays closer to the full model description in both the invariant and the transverse mass distributions.
We continue with a description of the W ± W ± final state, for which we depict the scalar case in Fig. 11 and the fermionic case in Fig. 12. The destructive interference in the invariant mass distribution is more pronounced than in the W ± Z case. Unfortunately, the actually observable transverse mass distribution suffers from two final-state neutrinos which leave the detector unobserved. This again leads to the migration of a substantial fraction of high-energy events in the invariant-mass distribution to lower values of the transverse mass, which wipes out the destructive interference signal. This migration also leads to significant deviations between the full-model description and the EFT setup for all values of the transverse mass, even far below threshold, m T (W ± W ± ) 2M R . Lastly, in the transverse-mass distribution also the unitarized EFT case is far from both the full-model description and the EFT curve, as the different behavior at high invariant-mass has a clear impact at lower values of the transverse mass. The choice of the unitarization procedure, which is arbitrary to a considerable extent, unfortunately impacts the transverse-mass distribution substantially. We note that the dilepton invariant mass distribution would show a behavior very similar to the transverse mass, i.e. throughout the whole range of masses the full model, the EFT description and its unitarized version differ. Again the lower right panel shows the binning as performed in the left panel of Fig. 6 in Ref. [3]. In the next-to-last bin, from 850 to 1050 GeV in m T (W ± W ± ), all curves Events/bin pp → W ± W ± jj → lνl νjj, l, l ∈{e, µ} J S = 6, M S = 600 GeV Figure 11: W W invariant mass (upper row) and llp miss T transverse mass (lower row) distributions for VBS production of same-sign W W jj events at the LHC. Three panels show expected event numbers per bin for a run 2 integrated luminosity of 137 fb −1 while the upper right panel is normalized to the SM. Bin width in the two left panels is 100 GeV, while on the lower right bin size is chosen as in Ref. [3]. The considered scenario is a scalar multiplet with J S = 6 and M S = 600 GeV. Events/bin pp → W ± W ± jj → lνl νjj, l, l ∈{e, µ} J F = 4, M F = 600 GeV Figure 12: Same as Fig. 11, but for a fermionic multiplet with J F = 4 and M F = 600 GeV.
for the J F = 4 fermionic case are compatible with the data, whereas the measurement, with a single event in the m T (W ± W ± ) > 1050 GeV bin, is visibly under tension with all non-SM curves depicted in the lower right panels of Fig. 11 and Fig. 12. It is again apparent that it is much easier to exclude the EFT description compared to the full model or the unitarized EFT description.

Constructive interference in ZZjj production
Within the VBFNLO framework we next consider the electroweak process pp → ZZjj → l + l − l + l − jj with l, l ∈ {e, µ}, which was used in Refs. [46,47] to bound aQGCs. For a qualitative comparison we again use a cut-flow similar to the experimental analysis, namely We generate the flavor combination l = l ∈ {e, µ} and multiply our results with a factor of 2, thus ignoring the problem of correctly assigning the leptons to the two Z boson candidates in 4e or 4µ events. Note that the actual experimental analysis in Ref. [46] is using a boosted decision tree, which also takes into account the z * l variable of Eq. (48) to enhance the fraction of vector boson scattering events over QCD background events. The more recent analysis in Ref. [47] is moreover extracting the bounds not from a VBS-enhanced region, but from a larger phase space region with only m jj > 100 GeV. We stick to the VBS-enhanced region defined above for our qualitative discussion, because otherwise triple gauge boson production processes would have to be considered as well. As a consequence we cannot directly compare against Ref. [47], which would anyhow need a detailed Monte Carlo simulation of all underlying, dominant background processes.
In Fig. 13 we show the ZZ invariant-mass distribution for the scalar case with J S = 6 and M S = 600 GeV which, in contrast to W W or W Z production, is fully accessible experimentally for ZZ → 4l. While the left panel is for a constant 100 GeV bin width, the central panel uses the binning employed in Refs. [46,47], and the ratio to the SM is shown in the right panel. In contrast to the same-sign W W and W ± Z final states the interference with the SM is constructive (see right panel) and the peak itself is more pronounced (see left panel). We emphasize again that a direct comparison to the experiments is non-trivial due to different regions in the phase space of the two jets. Still both the full model as well as the unitarized EFT description yield about three extra events in the last bin (m(ZZ) > 1200 GeV), which seems well compatible with the experimental data. The qualitative behavior of the J F = 4, M F = 600 GeV model is similar to the scalar model shown. However, deviations from the SM are somewhat smaller, with a single excess event expected in the last bin.

Implications for experimental analysis
The various cases discussed in the previous sections differ in details like constructive vs. destructive interference below threshold, or how large a threshold peak may be expected. However, there are a number of generic observations which can be made and which are important for BSM searches in VBS at the LHC. 1. One finds a very limited energy range, m V V 1.3M R for the diboson invariant masses, where the EFT description (even at dimension-8 level) adequately approximates the underlying UV-complete model. Within this validity range of the EFT description, cross section deviations from the SM stay below 10% even for isospins as large as J F = 4 or J S = 6, which are at the edge of perturbative behavior of the electroweak interactions. For less extreme isospin choices, deviations are even smaller, scaling like J 5 R .
2. In spite of small BSM effects in the EFT validity range, the overall BSM signal can be sizable. Enhancements by a factor of 10 in invariant mass or transverse mass distributions are possible, starting in the threshold region of the new physics. However, they require large isospin, J R , close to the perturbativity limit. For more "reasonable" multiplet assignments, the cross section increase at high m V V is substantially smaller. For example, three replicas of J F = 2, M F = 600 GeV Dirac fermions lead to an increase by about a factor 1.3 around m W ± W ± ≈ 3 TeV in the W ± W ± jj cross section, as compared to the factor 7 visible in Fig. 12 6. In the threshold region, 1.5M R m V V 3M R , the EFT description strongly underestimates the BSM signal of our toy model. This intermediate energy range thus provides an attractive opportunity to search for extra isospin multiplets in VBS. Because of the growth of amplitudes with J 5 R in VBS vs. J 3 R in dilepton production (as discussed in Section 2.5), VBS can actually be competitive with and complementary to analyses of Drell-Yan dilepton events in finding BSM hints in the (J R , M R ) parameter plane.
The above observations were made for a toy model with a large SU(2) L multiplet of otherwise inert heavy scalars or fermions. Many modifications of the model can be contemplated. For example, the new heavy matter fields might come in multiplets of an additional (confining) gauge interaction, which would add non-perturbative effects and resonances, similar to the formation of heavy quarkonium states, like the J/ψ or Υ, in QCD. Borrowing from quarkhadron duality in QCD, however, we would still expect the above features to approximately hold when smearing V V invariant mass distributions over energy intervals which are larger than the spacing of the resonances which are induced by the new confining interaction. In contrast, going to higher SU(2) L multiplets than contemplated in e.g. Fig. 3, i.e. entering the nonperturbative realm for electroweak interactions, might lead to new features, not discussed in this paper. Even then, the bounds established by the T u -model unitarization of the dimension-8 EFT would still present upper limits for observable cross sections in VBS.

Discussion and conclusions
In this paper we have investigated possible sources of anomalous quartic gauge boson couplings and their impact on vector boson scattering processes at the LHC. The T -operators of Eq. (16) arise naturally, at the one-loop level, in any BSM extension with extra fields which carry weak isospin, and they modify the scattering of transversely polarized weak bosons. Conversely, the appearance of field strength tensors in the T -operators signifies that they must be loopinduced [9]. Within the setting of renormalizable field theories, and ignoring the possibility of an extended electroweak gauge group, 10 the extra scalar and spin 1 /2 matter fields, which we have considered, constitute the most general source of these operators. We have avoided operators involving the SM Higgs-doublet field, like the M 0 -or S 0 -operators of Eq. (1), by only considering isospin and hypercharge representations which cannot form gauge invariant Yukawa type interactions with the Higgs doublet field (and possibly SM fermions) and by assuming small H † HΦ † Φ couplings in the UV-complete model.
Relaxing these model constraints will not significantly alter our results for the T -operators, and transverse VBS in general, as long as extra interactions do not induce very large mass splitting within the new matter multiplets. Our toy model represents a variant of natural dark matter models in the spirit of Ref. [29]. However, we have only considered relatively modest masses, of order 1 TeV, which would provide only a fraction of the observed dark matter in the universe. In a less constrained UV-complete model, in particular when allowing mixing with SM matter fields, direct collider searches for the extra multiplets as well as their dark matter impact would be strongly affected by additional interactions, leading to a vast and rich phenomenology. We were not interested in such issues here but rather have concentrated on the generic loop-induced effects of extra matter multiplets for weak boson interactions.
The one-loop effects of generic, degenerate scalar or Dirac multiplets have been calculated up to weak boson four-point functions, in the zero-hypercharge limit. In addition, the EFT representation of these results was derived, including all necessary operators up to energy dimension 8. Both the full calculation as well as the EFT low-energy approximation were implemented for on-shell VBS in the SU(2) L limit of the SM. In addition, an approximated version was added into VBFNLO [20,21] for the estimation of its impact in current experiments at the LHC. Assuming a single multiplet for simplicity, we have tuned the parameters of the multiplet, i.e. its isospin J R and its mass M R , to produce sizable, but not yet excluded deviations from the SM in VBS. Surprisingly, present VBS constraints are of comparable strength as bounds from dilepton production at the LHC, and considerably stronger than existing constraints from aTGCs as measured via vector boson pair production.
This somewhat surprising result (given the much higher event rate for Drell-Yan or qq → V V production at the LHC as compared to VBS) is due to the fact that large isospin multiplets induce one-loop effects in weak boson two-and three-point functions which rise as J 3 R only, while a J 5 R rise is found for four-point functions and, thus, also aQGCs. In addition, an accidental cancellation between the two relevant dimension-6 operators leads to a particularly small value for the λ-aTGC (see Eq. (34)). Degenerate, non-mixing matter multiplets at one-loop level thus provide another example that dimension-8 operators can be more important than dimension-6 operators in LHC phenomenology. One should keep in mind, however, that for small isospin multiplets (like J R = 1 /2 or 1) the one-loop BSM effects on weak boson vertex-functions are tiny, and models with large isospin matter fields are required to produce noticeable effects in VBS.
An upper bound on reasonable values of J R is provided by unitarity considerations, combined with perturbativity of the model. As shown in Section 3.1, the J 5 R growth of the one-loop VBS amplitude leaves the allowed range of the Argand circle for fermion multiplets with J F 5 and scalars with J S > 6. We have therefore considered a fermion or scalar with the largest isospin J R which is compatible with the unitarity bounds (i.e. J F = 4 and J S = 6) and a mass value such that at least one of the model's Wilson coefficients has a value in the ball-park of current EFT operator bounds obtained from aQGC measurements at the LHC [2, 3, 46, 49]. Both particle species show comparable impact in VBS processes.
Turning to individual VBS processes, we observe constructive interference with the SM for the case of ZZ production and destructive interference for W Z and same-sign W W scattering. However, the characteristic destructive interference in W Z and same-sign W W scattering is only visible as a function of the invariant mass of the diboson pair, which, unfortunately, is not a kinematic variable that is directly accessible in purely leptonic decays involving neutrinos. These interference effects, which remain at the 10% level even for isospins at the perturbativity limit, are thus hard to detect. Large BSM signals are possible in VBS at and above pair production threshold of the new particles, however, and such effects should be searched for at the LHC in the 1 to 2 TeV region. As discussed in the last section, an EFT approximation to the loop effects is only justified in the interference region, well below threshold. Dissecting the effects of the individual dimension-8 operators further, one finds large destructive interference between the different T i operators in W Z and same-sign W W scattering for the fermionic case. This implies that the experimental bounds on individual Wilson coefficients for these operators are only of limited use and tend to overly restrain the allowed parameter space of our fermionic model.
Because the validity region of the EFT approximation is confined to m V V 1.3M R , where the BSM effects are still small, the EFT is not the appropriate tool to search for extra large multiplets of Dirac fermions or scalars. Rather a full one-loop simulation of the extra matter fields should be used. An approximate implementation in VBFNLO, as discussed in Section 4 is available upon request and a full implementation is envisaged for the future.