9512.net

甜梦文库

甜梦文库

当前位置：首页 >> >> # QCD Corrections to Electroweak Vector Boson Scattering at Small Scattering Angles

DESY 03-053

ISSN 0418–9833

arXiv:hep-ph/0305137v2 26 Jun 2003

QCD Corrections to Electroweak Vector Boson Scattering at Small Scattering Angles

Kriszti? an Peters(a)1and Gian Paolo Vacca(b)2

(a) II.

Institut f¨ ur Theoretische Physik, Universit¨ at Hamburg, Luruper Chaussee 149, D-22761 Hamburg

(b)

Dipartimento di Fisica - Universit` a di Bologna and INFN - Sezione di Bologna, via Irnerio 46, 40126 Bologna, Italy

Abstract We investigate the role of a certain class of QCD corrections to electroweak vector boson scattering at small scattering angles and large energies. These are present since, from the perturbative analysis, the vector bosons may dissociate into quark-antiquark pairs giving rise to colour dipoles interacting through gluon exchanges. After the computation of the vector boson impact factors, we present expressions for the lowest order QCD scattering amplitude and for the leading logarithmic BFKL amplitude. Particularly we discuss numerical results for the process γγ → ZZ . The QCD corrections to the cross section resulting from the interference with the electroweak ones are estimated and compared with the leading pure electroweak part. Corrections resulting from the leading log BFKL amplitude are of the order of few percent already at the 0.5 ? 1 TeV energy range.

1 2

email: krisztian.peters@desy.de email: vacca@bo.infn.it

1

Introduction

The search for physics beyond the Standard Model (SM) and the determination of the crucial and still not understood mechanism of the electroweak symmetry breaking in the SM itself are major subjects of investigation in particle physics. There is the hope to increase the understanding with the help of the next generation of supercolliders, and major e?orts in this direction have currently been made. First insights are expected from the future large hadron collider (LHC) runs. On the other hand high precision measurements at the Next Linear Collider (NLC) are considered to be of great importance, with features complementary to LHC. Next Linear Collider is a generic name for a e+ e? machine operating in the energy regime up to 1 TeV, providing a very clean environment, and making available observables for many processes at an accuracy better than the percent level. In addition to the e+ e? collider mode, one may have the capability of running the NLC in a γγ collision mode via Compton backscattering of laser photons o? the linear collider electrons. One of the important processes one may consider in a linear collider is the scattering of vector bosons. Indeed, at the NLC the electron and positron bunches can radiate vector bosons which interact thereafter. Diverse scenarios are considered also in the Compton collider mode. Each incoming photon may turn into a virtual WW pair, followed by the scattering of one W of each pair to end with WW or ZZ states [1]. There are also the processes where initial photons scatter directly into arbitrary vector bosons. Vector boson scattering is related to the open question of mass generation. Since the longitudinal modes of the massive vector bosons re?ect the nature of the Goldstone modes of the unbroken theory, these may give important informations about the nature of electroweak symmetry breaking. Without the discovery of a scalar Higgs particle the vector bosons have to interact strongly at TeV scales, in order to preserve unitarity. In this scenario the vector boson scattering will be one of the key processes to be studied in the nearest future [2]. Another important question one may pose while searching for new physics is the existence of anomalous triple and quartic vector boson couplings which may a?ect vector boson scattering. Contact interactions are regarded as a possible extension of the SM, which can be well investigated with vector boson scattering [3]. The natural order of magnitude of the anomalous couplings [4] is small and, in order to eventually extract them, one needs to both measure and compute theoretically the SM cross sections with a precision better than 1%. In this context the γγ → ZZ process has a special role. Since there is no perturbative contribution at the tree level, such a process is sensitive to the particles and the new physics phenomena which contribute through radiative corrections [5]. This leads to an independent and complementary method to the direct production of new particles. In the same process the detection of CP violating phases [6] or additional e?ects due to the exchange of Kaluza-Klein gravitons in large extra-dimension scenarios [7] have been discussed. In summary, vector boson scattering processes may be seen as a central tool to probe physics beyond the SM and to reveal the nature of symmetry breaking. In order to exploit these features, a high precision is needed, both on experimental and on theoretical side. Calculations of these processes at the lowest available order may not be accurate enough.

1

One loop corrections have been calculated for all those vector bosonic processes [8] where a tree level contribution is present. For the γγ → ZZ process, which starts at one loop in αEW , no higher order corrections are available yet. The question naturally arise, if this accuracy is high enough to disentangle new physics from the SM, providing a window to non-accessible energy regimes. The aim of this paper is to address the question if the QCD corrections could play a role in vector boson scattering processes. Naively this seems not to be the case, since they arise only at higher loop corrections. On the other hand, even if the lowest order nonsuppressed QCD correction to the amplitude for this process is smaller by O(α2 s ) w.r.t. the electroweak one, in the cross section the interference among the two can reach the percent level. Moreover at high energies there are kinematical regimes wherein the emission of soft gluons generates large logs in the energy which make the QCD amplitude rising with energy. The electroweak cross section of V V → V V (with V = W, Z, γ ) gains its biggest contribution from small scattering angles, and with rising scattering energy the cross sections are approaching a constant in computations up to the one loop level. This contribution is coming from diagrams where only spin-1 particles are exchanged in the t-channel. At leading order and O(α2 EW ) these exchanged particles are only electroweak vector bosons. The same con?guration arises for QCD corrections with gluons in the t-channel. The ?rst gluonic contributions, which arise at the two loop level, are suppressed since these gluons are always accompanied by fermions in the t-channel. The non-suppressed QCD contributions occur at higher orders. These emerge when the vector bosons ?uctuate into a quark-antiquark pair and these dipoles interact through gluons. At the lowest order, when two gluons are exchanged in the t-channel, one obtains a contribution constant in energy while higher order corrections provide large logarithms in the energy. These contributions cannot be neglected at large energies. The extreme attitude is to consider a resummation described by the BFKL equation [9] which gives an upper bound estimate of these e?ects. We study in detail these QCD corrections for the reaction γγ → ZZ , whose cross section up to now was only studied at leading order (one loop in the electroweak sector) by G. Jikia [10]. His results show that the cross section is approaching a constant for high energies and its main contribution comes from the helicity conserving amplitudes at small scattering angles. This is exactly the kinematical region where the above discussed QCD contributions have their biggest in?uence. We begin with the analytic calculations for an arbitrary bosonic process. In the next section the two gluon exchange amplitude is computed, including the calculation of the non-forward impact factor. Since we consider the high energy limit the amplitude factorise into two independent impact factors, built from the bosonic wave functions which describe the probability that a vector Boson ?uctuates into a quark anti-quark pair. In the third section the cross section of the resummed gluonic contribution is presented, i.e. the BFKL amplitude. In order to give some numerical estimates, we concentrate on the de?nite process γγ → ZZ . In the fourth section we calculate the EW and QCD parts of the amplitude in the Regge limit. Next we evaluate the corrections numerically. Because of the strong scale dependence we can give only estimates of the corrections, which we ?nd to be of the order of a few percent. Finally we present our conclusions. Our result suggest that, in particular in the small angle region, possible deviations from EW calculations cannot be interpreted as signals of new physics without taking into account the QCD corrections.

2

??

??? ?·?

¨ ·?

??

¨

???

??

??

???

?·?

?

Figure 1: The two gluon exchange of the V V → V V elastic cross section.

2

Impact factors and lowest order QCD contribution

At high energies the lowest order leading QCD contribution to the vector boson scattering is pictured in Fig.1 where the lowest order diagrams with two-gluon exchange are illustrated. The vector boson ?uctuation into a quark anti-quark pair is described by the product of the two boson wave functions. This quantity is called the (non-forward) vector boson impact factor. Due to the high-energy factorization one can compute ?rst the impact factors and then integrate these with the gluon propagators. In the high-energy limit, corresponding to s going to in?nity and t ?xed (Regge limit), we calculate ?rst the s-channel discontinuity and restore after the amplitude through dispersion relations. Thus the intermediate quark lines inside the impact factors are taken on-shell. This procedure and the high-energy factorization simpli?es the calculation signi?cantly. The s-channel discontinuity of the scattering amplitude takes the form: Discs A = is d2 k Φ(q, k, l) Φ(p, ?k, ?l) , (2π )3 (k + l)2 (k ? l)2 (1)

where Φ is the impact factor and the boldface marked letters are the transverse part of the euclidean momenta. As usual, the integration of the two longitudinal components of the loop momenta are absorbed into the de?nition of the impact factors. The remaining integrals may be computed numerically. The momenta of the incoming bosons are q + l and p ? l, while the outgoing bosons have momenta q ? l and p + l. We use Sudakov variables with the light cone vectors p′ = p + yq ′ , y ? ?p2 /2p′ · q ′ and q ′ = q + xp′ , x ? ?q 2 /2p′ · q ′ , while s = (p + q )2 ? 2p′ · q ′ denotes as usual the squared center of mass energy of the V V scattering process. Then one parametrizes any 4-momentum u as u = αu q ′ + βu p′ + u⊥ . (2)

The impact factor for incoming and outgoing photons has been calculated recently [11] and since we adopt the momentum assignment which was given there, we refer the reader to this paper for further details of the calculation. 3

The incoming and outgoing bosons are characterized by the invariants Q2 + = ?(q + 2 2 2 l) , Q? = ?(q ? l) and by similar expressions associated to the other two bosons involved in the process. For the full process the Mandelstam variable t = (2l)2 has to be considered. 2. In the case of on-shell bosons the virtuality reduces to the boson mass, i.e. Q2 = ?MV We remind that in the limit of Regge kinematics one can neglect terms which are suppressed by powers of t/s and Q2 ± /s, leading to important simpli?cations. Let us consider the polarization vectors for the upper incoming and outgoing bosons. Their longitudinal polarization in terms of the Sudakov variables are given by ?L (q ± l) = 1 2l2 (1 ± αl ) q ′ + x ? βl ? ⊥ |Q± | s p′ ± l⊥ , (3)

while the transverse polarisation vectors are ?T (q ± l) = ?⊥ ±

(h) (h)

2l⊥ · ?⊥ (q ′ ? p′ ± l⊥ ) , s

(h)

(4)

where h = ± denotes two helicity states and 1 (h) ?⊥ = √ (0, 1, ± i, 0) . 2 (5)

Again, for the bosons appearing in the lower impact factor the same expression holds with an exchange of q and p and l → ?l. A di?erence between the photon and the other bosons is given by the chiral coupling to the quarks of the latters. Consequently extra terms in the traces appear and the ones proportional to odd powers of γ 5 do not contribute because of charge conjugation invariance. Thus, regarding the impact factors, the W and Z bosons behave both in the same way as the photon, i.e. we can use the result from [11], and only the magnitude of the coupling to the quarks depends on the nature of the boson. The impact factor where both of the bosons have longitudinal polarisations takes the form:

2?1 ΦLL (q, k, l) = αs Nc

f

√ C f 4 4π |Q+ ||Q? | ×α2 (1 ? α)2

1

dα

0

d2 u × (2π )2 1 1 ? ? ? D1 D2 ,

(6)

1 1 + ? + D1 D2

while for transverse polarized bosons we have: ΦT T (q, k, l) = αs ×

(ij ) 2?1 Nc

Cf

f

√

1

4π

0

dα

d2 u (2π )2 · ?? j 1 1 + ? + D1 D2 1 1 ? ? ? D1 D2 .

(7)

?4α(1 ? α) ?i · N+ N+ 2 1 + ? + D1 D2

N+ N+ 1 2 ? + + D1 D2 ·

N? N? 1 2 ? ? ? D1 D2 + m2 f

+ ?i · ?? j

N? N? 1 2 ? ? ? D1 D2

4

Finally, if the incoming boson is transverse polarized while the outgoing one is longitudinal, the impact factor reads:

2?1 ΦLT (q, k, l) = αs Nc (j )

f

√ Cf 2 4π |Q+ |

1

dα

0

d2 u × (2π )2 N? N? 1 2 ? ? ? D1 D2 · ?? i ,

(8)

× α(1 ? α)(1 ? 2α)

1 1 + ? + D1 D2

and viceversa for the incoming boson longitudinal and the outgoing one transverse, one has

2?1 ΦT L (q, k, l) = αs Nc (i)

f

√ Cf 2 4π |Q? |

1

dα

0

d2 u × (2π )2 N+ N+ 1 2 ? + + D1 D2 1 1 ? ? ? D1 D2 ,

(9)

× α(1 ? α)(1 ? 2α) ?i ·

√ where ?j = 1/ 2(1, ±i) are two-dimensional polarization vectors corresponding to the two transverse polarizations j = ±. We have also used:

± 2 D1 = (u ± (1 ? α)l)2 + α(1 ? α) Q2 ± + mf ± 2 D2 = (u ? k ? αl)2 + α(1 ? α) Q2 ± + mf

N± = u ? k ? αl . 2

N± = u ± (1 ? α)l 1

(10)

2 Marked in boldface are the Euclidean form of the transverse momenta, i.e. u2 ⊥ = ?u < 0. The constant Cf is the product of the coupling of the incoming and outgoing boson to the quarks of ?avor f . The impact factors have been written in a factorised form where the two parts correspond to the boson wave functions, describing the probability that a boson dissociates into a quark-antiquark pair. In order to obtain the lowest order QCD amplitude, one has to substitute these impact factors into eq.(1) and integrate over the gluon momenta. The two helicity-?ip impact factors are in general non-zero. They vanish only for forward scattering, where l = 0. This is due to the di?erent symmetry behaviour of the transverse and longitudinal boson wave functions. This property is more trasparent in the coordinate space formulation, eq.(34) of [11]. Here the impact factor is written as a convolution of the two boson wave functions and the dipol interaction. The dipole interaction is invariant under the transformation of the dipole size vector r → ?r, and so is the longitudinal wave function, eq.(22) of [11]. In contrast to this the transverse wave function, eq.(28-31) of [11], is antisymmetric under this transformation, since it is proportional to r. The convolution of two transverse wave functions is therefore symmetric under this transformation, whereas the convolution of a transverse and a longitudinal one is antisymmetric, leading to the vanishing result for l = 0. If l is non-zero, this symmetry properties are broken, resulting in a non-vanishing ΦT L and ΦT L .

5

??

?

? ·?

¨? ? ?? ·?

?

? ? ??

?

? ??

??

??

?

·?

?

?

??

?

??

? ?· ?

¨?

?

? ?? ?

?

? ? ??

?

? ??

?

Figure 2: The BFKL Pomeron exchange in the V V → V V elastic cross section.

3

Leading log (BFKL) contribution

At high energies and close to the forward region (Regge limit), higher order contributions in αs start to be important. Diagrams of ladder topology, built with (non elementary) n reggeized gluons, can contribute up to the order αn s ln (s/s0 ). On retaining the leading logarithmic (LL) contributions and performing a resummation a rough estimate (overestimate) of the total QCD corrections is obtained. Nonetheless one expects to ?nd the QCD contribution to the process in the window constrained by the simple two gluon exchange on one side and the resummed LL BFKL Pomeron exchange on the other. This possibility is related to the property of the factorisation of the QCD scattering amplitude, in the Regge limit, in terms of the impact factors associated to the external particles and the BFKL Green’s function, as illustrated in Fig. 2. For the physical process we want to study the impact factors have been calculated in the preceeding section and the discontinuity of the LL resummed amplitude is simply obtained by inserting the non-forward BFKL Pomeron Green’s function in the place of the two gluon propagators. Indeed one has Discs A = is d2 k d2 k′ Φ(q, k, l)G(y |k + l, k ? l; k′ + l, k′ ? l)Φ(p, ?k′ , ?l), (2π )3 (2π )3 (11)

where y = ln s/s0 is the rapidity variable and s0 is a typical scale in the process. Clearly, for αs → 0, when all rungs of the BFKL resummation decouple, the Green’s function reduces to describing the two gluon exchange, which means that we use the following normalisation:

αs →0

lim G(y |k + l, k ? l; k′ + l, k′ ? l) =

(2π )3 δ(2) (k ? k′ ) 2 2 (k + l) (k ? l)

(12)

In the following we brie?y give the ingredients needed to construct the BFKL Green’s function. The eigenfunctions Eh,h ? of the BFKL kernel, which appears in the integral Bethe-Salpeter like equation resumming all the LL terms, are well known in coordinate space, wherein their form is dictated by conformal invariance. They are given by Eh,h ? (r 10 , r 20 ) = r12 r10 r20 6

h ? r12 ? r? r10 20 ? h

,

(13)

? = (1 ? n)/2 + iν (h? = 1 ? h ?, h ? ? = 1 ? h). with r ij = r i ? r j etc., h = (1 + n)/2 + iν , h Let us note that standard complex notation for two-dimensional vectors is used on the right-hand side. The conformal weights h, decomposed in n (integer conformal spin) and ν (real), label the principal series of unitary representations of the Moebius group. Fourier transforming to momentum space (we use the Lipatov’s convention which assigns a 1/(2π )2 to any coordinate integration) one ?nds [15]: d2 r 1 d2 r 2 E ? (r 1 , r 2 )ei(k 1 ·r1 +k2 ·r2 ) = C X (k1 , k2 )+(?1)n X (k2 , k1 ) . (2π )2 (2π )2 h,h (14) The coe?cient C is given by ? ? (k1 , k2 ) = E hh C= (?i)n ? ? )Γ(1 ? h)Γ(1 ? h ? ). hh(1 ? h)(1 ? h (4π )2 (15)

The functions X in complex notation can be expressed in terms of hypergeometric functions: X (k1 , k2 ) = k1 2

? ?2 h ? k2 2 h ?2

F 1 ? h, 2 ? h; 2; ?

(16) This analytic form does not contain any term of the type δ2 (k1 ) or δ2 (k2 ) which are present in the coordinate representation eq.(13). The impact factor of a colorless external particle vanishes for zero momentum of any gluon exchanged and therefore the delta-function type terms do not contribute. The Pomeron intercept has the form αP (0) = 1 + χ(ν, n), where χ(ν, n) = χh = α ? s 2ψ (1) ? ψ ( 1 + |n| 1 + |n| + iν ) ? ψ ( ? iν ) 2 2 , α ?s = Nc αs , π (17)

? k1 ? k2

? 2?h ? ; 2; ? k2 F 1 ? h, k1

.

are the kernel eigenvalues corresponding to the eigenfunctions in equations (13)-(14). We shall use the Pomeron Green’s function in the momentum representation. Starting from the form in coordinate representation [16] one can perform a Fourier transform [15] and making use of the Casimir operator properties of the M¨ obius group, one ?nally obtains [17] for the non-amputated Green’s function: ? 2 (y |k1 , k2 ; k1′ k2′ ) = (2π )3 G ? ? (k1 , k2 )E ? ? ? (k1′ , k2′ ) , d?(h) ey χh Nh × (2π )2 E h,h h,h d?(h) =

n

(18)

where we use the measure in the conformal weight space following normalization factor Nh : Nh =

dν with the

(ν 2 + n2 /4) . [ν 2 + (n ? 1)2 /4][ν 2 + (n + 1)2 /4]

(19)

The (2π )3 factor in front of the integral comes from the normalization ?xed by eq.(12). Let us ?nally rewrite the expression of the amplitude discontinuity in eq.(11) by inserting in it the relation (18). On de?ning the conformal weight representation for the impact factors according to Φh (q, l) = d2 k ? ? (k + l, k ? l) , Φ(q, k, l)E h,h (2π )3 7 (20)

together with the primed Φ′ h (q, l) computed with the conjugate eigenstate, one has Discs A = is(2π )5 d?(h) Nh ey χh Φh (q, l)Φ′ h (p, l) . (21)

We stress that due to high energy factorization the energy dependence is encoded in the LL BFKL Green’s function and the only other dependence in the momenta of the incoming particles is through the invariant external masses on which the impact factors depend. In the last expression the dominant contribution for large y can be easily extracted with a saddle point analysis for n = 0 around ν = 0 and therefore reduces to √ 2π 5 y χ0 Discs A = is(2π ) 16e Φ0 (q, l)Φ′ (22) 0 (p, l) , (2ay )3/2 where χ0 = 4 ln(2)? αs is the leading BFKL intercept and a = 14ζ (3)? αs . The forward limit presents a well known cusp behavior and it must converge to the result obtained on using a forward BFKL Green’s function GF (y |k, k′ ) = (2π )3 ? ? (k′ ) , ? ? (k)E d?(h) ey χh E h,h h,h

(F ) (F )?

(23)

which is written in the spectral represention with the forward BFKL amputated eigen√ ? ?1 ? (F?) (k) = 1/(π 2)kh?1 k? h . functions given by E h,h The saddle point approximations can be carried on in a way similar to the non forward case. We shall de?ne Φh (q ) =

′ (F )

(F )

1 ? (F ) d2 k Φ(q, k, 0) 2 E ? (k) , 3 (2π ) k h,h

(24)

as well as a Φh (q ) for the scalar product with the conjugate Pomeron state and write the discontinuity in the forward case as √ ′ (F ) 2π (F ) (F ) 3 y χ0 Φ0 (q )Φ0 (p) . (25) Discs A = is(2π ) e 1 / 2 (2ay ) Let us note that in the transition from the nonforward to forward physics the perturbative QCD analysis is increasingly a?ected by the long distance interactions (contributions for all values of the impact parameter), since gluons are massless. This implies that some care is required to treat the general non-forward expression in this limit. The saddle point approximation for n = 0 around ν = 0 leading to eq.(22) is accurate only for l down to values of the order of one. Strictly speaking, performing the saddle point approximation of eq.(21) in the region of very small l one needs to shift the saddle point which becomes 2, l dependent. In our numerical calculation we have used eq.(22) down to |2l|2 = 1/2 MZ and then we have extrapolated to the forward value, eq.(25). A more precise analysis, avoiding the use of the saddle point approximation is of course possible but it is not needed for the order of magnitude estimate which is pursued here. By the way few checks have been done with an exact Green’s function computation, con?rming such an expectation.

8

? ?

?

? ?

?

Figure 3: The surviving contributions to the leading electroweak amplitude in the Regge limit for real γ and transverse polarised Z.

Finally, in order to evaluate the elastic cross sections we need to obtain the amplitude from the computed discontinuity. Instead of using the full dispersion relation technique, a simple evaluation compatible with the approximation adopted is given by A? 1?i π ? 2 ?y 1 Discs A , 2 (26)

as provided by the signature factor in the partial wave amplitude analysis.

4

The process γγ → ZZ in the small angle limit

We shall henceforth concentrate on a de?nite process, γγ → ZZ , and compute for it a numerical estimate of the QCD contributions. We remind that we are interested in the kinematical region of small scattering angles, i.e. scattering close to the forward region, so that one can do the calculation in the high-energy approximation suited for the Regge limit: terms which are suppressed by powers of t/s can be neglected, with a signi?cant simpli?cation in the computation. We shall in fact use the results presented in section 2 and 3.

4.1

The leading electroweak contribution

Let us ?rst analyze the pure electroweak contributions to the process. This process starts at one loop and was calculated by Jikia [10]. In order to calculate the QCD corrections to γγ → ZZ one has to know the leading electroweak contribution at the amplitude level, add to it the leading QCD amplitude and take the modulus squared. The main contribution to the corrections will come from the interference term. In the high energy limit where t is kept ?xed most of the contributions are suppressed by powers of s. The surviving graphs are pictured in Fig. 3. Also for these only the helicity conserving amplitudes survive, with both of the bosons transverse polarised. One may con?rm this by taking the limit s to in?nity and keeping t ?xed within the expressions of [10]. In order to simplify the calculation we compute ?rst the s-channel discontinuity, i.e. the s-channel internal W lines are taken on-shell. Here the s-channel discontinuity coincides with twice the imaginary part of the amplitude. With the same assignment of the momenta already used in the last sections the imaginary part of the box contribution reads: Im Aγγ →ZZ =

(1)

1 2

d4 k (0) (0)? 2 2 δ[(q ? k)2 ? MW ] δ[(p + k)2 ? MW ] Aγγ →W W AW W →ZZ , (27) (2π )2 9

where a summation over the W boson helicities is implied. Performing the loop integration and retaining only the leading power of s we get: Im ABox = c2 w s α2 em 2 sw 8π t2

2 ? 4tMW

where |2l|2 ≈ ?t was used. Equation (28) yields immediately the real part by just noting that ln(?s) = ln s ? iπ . Thus the reconstructed (leading) full box amplitude reads: ABox = ?s s ln π t c2 w α2 em 2 sw 8π

2 t2 ? 4tMW

ln ?

??

? ? ?t +

2 t2 ? 4tMW 2 t2 ? 4tMW

?t ?

?2 ? ? ? ?,

(28)

The amplitude of the crossed diagram is obtained from the box contribution (29) by replacing s by u ≈ ?s. Adding the box and crossed contributions, at the end one is left with the imaginary part of the box diagram:

?+? A++++ = A+ = i Im ABox EW EW

?? ? ?t + ln ??

t2

?

2 4tMW

?t ?

2 t2 ? 4tMW

?2 ? ? ? ?.

(29)

(30)

Thus in the Regge limit the leading order electroweak contribution is purely imaginary (s ? ?t), and the cross section is constant in s. Considering the limiting behavior of 2 and for the case t = 0 one obtains the expressions eq.(30) in the region s ? ?t ? MW presented in [10].

4.2

QCD corrections

The amplitude of this process in the high energy limit is illustrated in Fig. 2. The ingredients needed have been discussed in sections 2 and 3. To proceed with an explicit calculation we have to specify the couplings Cf , which depend on the external bosons, and the virtualities appearing in the impact factors given in equations (7) and (9). The starting expressions are for massive quarks. For practical purposes only the top quark mass has to be taken into account because of the large Z mass. The integrals which correspond to the quark loop have to be performed. The last integration we compute is the scalar product of the impact factors with the BFKL states eq.(20) and (24). We calculate this process for all external bosons on-shell. In such a case we have Q2 + 2 identically zero and Q2 ? = ?MZ . The coupling constants Cf - the product of the γ -quark and Z-quark couplings - are given by: Cf = 4παem 3 2 qf Tf ? 2 qf sin2 θw 2 sin θw cos θw . (31)

First we consider the scattering of the photons into transverse Z bosons. To proceed we have to solve the remaining integrations for the impact factor ΦT T in eq.(7). In order to compute the integral over the two dimensional euclidean vector u we employ a standard Feynman parametrization and rewrite eq.(7) as: ΦT T (q, k, l) = αs

(ij ) 2?1 Nc

Cf

f

√

1

4π

0

dα

d2 u (2π )2

1 0

dx ×

(32)

10

× ?

2 2 2 ? ? 4α(1 ? α)[?i · u ?? j · u ? x(1 ? x)?i · pa ?j · pa ] ? (u ? x(1 ? x)pa + mf )(?i · ?j ) 2 2 2 ? ? 4α(1 ? α)[?i · u ?? j · u ? x(1 ? x)?i · pb ?j · pb ] ? (u ? x(1 ? x)pb + mf )(?i · ?j ) 2 2 [u2 + x(1 ? x)p2 b ? xm + mf ] 2 2 [u2 + x(1 ? x)p2 a ? xm + mf ]

? .

2 + i?). Strictly Here we have de?ned pa = ?k ? (1 ? 2α)l , pb = 2αl and m = α(1 ? α)(MZ speaking, at this stage of the calculation m is taken negative and subsequently, after having performed the integrations, analytically continued to positive values obtaining the above de?nition. The ?rst fraction corresponds to the two diagrams where one of the gluons is attached to the quark and the other to the antiquark. The second fraction corresponds to the two diagrams where the two gluons are attached to the same (anti)quark. Each of them would lead to a divergent integral but, or course, the sum in eq.(32) gives a ?nite result. On performing the u integration we end up with:

ΦT T (q, k, l)

(ij )

= αs × + +

2?1 Nc

Cf

√

1

4π

0

x(1 ?i · ?? j

f ?i · pa ?? j · pa 2 ? x)pa ? xm + 2

dα 2π

1

dx

0

? 4α(1 ? α) x(1 ? x) × + + . (33)

m2 f

2

?

[α + (1 ? α) ] ln ?

2 x(1 ? x)p2 a ? xm + mf 2 x(1 ? x)p2 b ? xm + mf 2 x(1 ? x)p2 b ? mf

?i · pb ?? j · pb 2 x(1 ? x)pb ? xm + m2 f

2 x(1 ? x)p2 a ? xm + mf

2 x(1 ? x)p2 a ? mf

2 x(1 ? x)p2 b ? xm + mf

The computation of the x-integration is straightforward. Due to the lenght of the result we present here only the limit mf → 0 which preserves a signi?cant simpli?cation compared to the massive case. The corresponding helicity conserving result with massive fermions is presented in the appendix. After the integration in the Feynman parameter variable, one obtains in the zero quark mass limit for the separate helicity states: ΦT T (q, k, l) ΦT T (q, k, l) ΦT T (q, k, l)

(+?) (??) (++)

= αs

2?1 Nc

Cf

f

√

1

4π

0

p2 ? m dα 2 [α + (1 ? α)2 ] ln a 2π p2 b ?m dα 4α(1 ? α) × 2π ?

y px b pb p2 b

(34) (35)

= ΦT T (q, k, l) = i αs ×

2?1 Nc

(++)

Cf

f

√ 4π

0

1

(36)

y px a pa p2 a (+?)

1+

ΦT T (q, k, l)

(?+)

p2 m a?m ln 2 pa ?m

1+

p2 m b ?m ln 2 ?m pb (37)

= ?ΦT T (q, k, l) ,

with (px , py ) = p. For the helicity conserving parts in the massless quark limit the remaining α integration was also performed analytically and is presented in the appendix. For the result with massive quarks a numerical integration has been performed.

11

In the same way one can calculate the scattering of real photons into longitudinal Z-bosons and obtain (in the mf → 0 limit): ΦT L (q, k, l)

(i)

= αs ×

2?1 Nc

√ Cf 2 4π MZ

f

1 0

dα α(1 ? α)(1 ? 2α) × 2π . (38)

?i · pb p2 ?i · pa p2 a?m b ?m ln ? ln 2 2 ?m pa ?m pb

In order to compute the two-gluon or the BFKL amplitude one has to substitute these expressions into eq.(1) or eqs. (20, 24).

5

Numerical results

In this section we present numerical results on the QCD corrections for ZZ pair production in γγ collisions. Here we consider full circular polarization of the incoming photons. The relative QCD contributions are presented for transverse Z bosons only, because in the high-energy approximation the helicity changing parts of the electroweak amplitude are supressed by powers of s and thus are negligible. Nevertheless we shall display the helicity changing part of the BFKL amplitude as the pure QCD contribution alone. The parameters αW = α/s2 W , α = 1/128, mZ = 91.2 GeV, mt = 174.3 GeV and αs (s0 ) have been used throughout the numerical computations. We have been neglecting all quark masses apart from the top quark mass in the numerical calculation, since they are neglible compared to the Z boson mass. The inclusion of the top quark mass reduces the QCD amplitude by 25%. In terms of the amplitudes (30,21,26) the cross section reads: 1 dσ |AEW + AQCD |2 = 2 dp2 16 πs T (39)

where pT = |2 l| is the exchanged transverse momenta. The angular distribution of the cross section is obtained by the relation: s dt = β d cos θ , 2 (40)

2 /s. For small scattering angles in the high-energy approximation the with β = 1 ? 4MZ two incoming and two outgoing particles do not ’communicate’ with each other because of the huge separation in rapidity. As a result, the amplitude with the helicity ++ → ++ is equal to the amplitude with helicity +? → +? etc. Since the helicity conserving amplitudes completely dominate the cross section, after the summation of outgoing helicities the result is the same for all the possibilities of incoming helicities, i.e. ++ → TT is equal to +? → TT etc. In the following γγ → ZT ZT stands for all these possibilities as well as for unpolarised scattering. In Fig. 4 the angular distribution of the leading electroweak contribution for transverse √ Z pair production for polarised γγ collision at sγγ = 1 TeV is given. The solid line is the result of the calculation in the high-energy approximation, presented in section 4.1. √ The dotted line is the result of the exact calculation, presented also in [10] for sγγ = 500

12

3000

γγ → ZT ZT

2500

EW approx. EW exact

dσEW / dcos(θZ) [fb]

2000

√sγγ = 1 TeV

1500

1000

500 1 0.99 0.98 0.97 0.96 0.95 0.94 0.93 0.92 0.91 0.9

cos(θZ)

Figure 4: Angular distribution of the leading electroweak contribution for transverse Z pair pro√ duction in polarised γγ collision at s = 1 TeV. The solid line is the result of the high-energy approximation, the dotted line is the result of the exact calculation.

GeV. For cos θ ≈ 0.9 the result of the high-energy approximation has a deviation of ? 15% from the exact calculation, for smaller scattering angles the deviation is decreasing. With rising center of mass energy the accuracy of the Regge calculation is continously increasing. Let us consider the QCD corrections. In the Born analysis, from the naive counting of the couplings (see Fig. 1), one expects in the interference term a suppression, with respect to the one loop EW case, due to the factor α2 s and an additional suppression of a factor ≈ 4 due to the di?erence in the EW couplings. This would lead to corrections of less than a percent. At high energies in the Regge limit we expect an enhancement due to the appearance of some large logarithms, which may be resummed in the BFKL scheme. We stress that the resummed leading log QCD corrections at the lower energies we consider are overshooting what is expected from the real contribution, nevertheless they provide a ?rst estimate. We remind, as discussed previously, that the BFKL resummation is evaluated in the saddle point approximation. Fig. 5 displays the di?erential cross section with respect √ to the exchanged transverse momenta at a center of mass energy of sγγ = 1 TeV. The solid line is the pure electroweak cross section, and the dashed line is the one with the QCD corrections included, eq.(39). The cross section has a slight rise approaching the forward region. This slight rise comes from the electroweak part since it completely dominates the amplitude. As in the following it will be visible, for the QCD part of the cross section the rise approaching the forward region is much stronger. Therefore the QCD corrections are enhanced close to the forward region. This behavior is shown in Fig. 6, where the QCD corrections to the di?erential cross

13

x 10 0.6 0.55 0.5

2 dσ / dpT [fb/GeV ] 2

-2

√sγγ = 1 TeV s0 = (1/2 MZ) EW EW + QCD

2

0.45 0.4

0.35 0.3

0.25 0.2 0.15 0.1 0 0.5 1 1.5 2

2 2 pT / MZ

2.5

3

3.5

4

Figure 5: The pure √ EW di?erential cross section (solid line) versus the exchanged transverse momenta at s = 1 TeV and the one which include the QCD corrections (dashed line). For the BFKL contribution the scale s0 = (1/2MZ )2 was used.

section relative to the pure electroweak contribution at center of mass energies of 500 GeV 2 (a), 800 GeV (b), 1.5 TeV (c) and 3 TeV (d) are shown as functions of p2 T /MZ . At these 2 energies p2 T /MZ = 4 corresponds to values of cos θ (where θ is the scattering angle) of 0.75 (a), 0.90 (b), 0.97 (c) and 0.99 (d) respectively. This relative correction is de?ned as ?= dσQCD+EW dσEW ? 2 dpT dp2 T dσEW . dp2 T (41)

Since for leading order BFKL the scale is not ?xed, there is an uncertainty which a?ects the magnitude of the QCD corrections. The solid lines in Fig. 6 are results for the scale 2 . At the scale s = M 2 s0 = (1/2 MZ )2 while the dashed lines correspond to s0 = MZ 0 Z √ the minimal center of mass energy s, for which the saddle point approximation for the leading logarithmic calculations gives a reliable result, is 800 GeV. Corrections resulting from leading order BFKL contribution are at the percent level. For higher pT they are about one percent and approaching the forward region (smaller pT values) they rise and lead to corrections of a few percent. Close to the forward region the sign of the corrections is changing as a result of the fact that for small l the imaginary part of Φh (eq.20) gets an enhancement. Fig. 6 (a-d) summarizes the dependencies of the corrections resulting √ from leading order BFKL calculations with respect to the center of mass energy s, the exchanged transverse momenta pT and the scale s0 . For rising center of mass energy the corrections increase. The solid lines (s0 = (1/2 MZ )2 ) in Fig. 6(a) and (d) show that the corrections may rise in this regime by a factor of 2 if the energy increases from 500 GeV to 3 TeV. The scale s0 implicates also a strong dependence. In fact close to the forward 14

6

6

5

√sγγ = 500 GeV

5

√sγγ = 800 GeV

4

4

s0 = (1/2 Mz) s0 = M2 z

2

3

3

? [%]

2

? [%]

(a)

2

1

1

0

0

-1

-1

(b)

1 2

2 2 pT / MZ

-2

0

3

4

-2

0

1

2

2 2 pT / MZ

3

4

6

6

5

√sγγ = 1.5 TeV

5

√sγγ = 3 TeV

4

4

3

3

? [%]

2

? [%]

(c)

2

1

1

0

0

-1

-1

(d)

1 2

2 2 pT / MZ

-2

0

3

4

-2

0

1

2

2 2 pT / MZ

3

4

Figure 6: QCD corrections to the di?erential cross section relative to the pure EW contribution √ for di?erent scales s0 and center of mass energies s. The relative correction is de?ned dσ +EW EW EW as ? = ( QCD ? dσ )/ dσ . dp2 dp2 dp2

T T T

15

4 3.5

γγ → ZT ZT

(σEW+QCD - σEW) / σEW [%]

3 2.5 2 1.5 1 0.5 0

s0 = (1/2 Mz) s0 = M2 z

2

1000

1500

2000

2500

3000

3500

4000

4500

5000

√sγγ [GeV]

Figure 7: The QCD corrections to the integrated cross section relative to the pure EW contribution 2 2 for di?erent scales. The integration region was choosen as 1/4 MZ < pT < 4 M Z .

2 region the corrections also grow by a factor of 2 if the scale is decreasing from s0 = MZ 2 to s0 = (1/2 MZ ) . The biggest dependence is coming from the exchanged momenta or equivalently from the scattering angle. The corrections have a strong rise with decreasing scattering angle or pT and change sign close to the forward region. In Fig. 7 relative corrections to the integrated cross section (σEW +QCD ? σEW )/σEW are plotted as a function of the center of mass energy. These are also at the order of percent level. The corrections rise with increasing center of mass energy. The magnitude of the corrections is again scale dependent and the rise is stronger for smaller scales, as expected. Since the calculations are done in the Regge limit we integrated eq.(39) in the 2 < p < 4 M 2 . An integration down to p = 0 would change the result region of 1/4 MZ T T Z by less then 0.1 %, since the sign of the corrections changes very close to the forward region, (see Fig. 6). If one integrates over all scattering angles the corrections become less important, but since the main contribution of the cross section is coming from small scattering angles the result would not be a?ected strongly. In the following ?gures some properties of the pure QCD part of the amplitude are discussed. In Fig. 8 the absolute value of the QCD amplitudes is displayed for the helicity conserving case. The solid line is the BFKL (resummed) contribution while the dashed line √ is the Born (two-gluon) contribution. At s = 1 TeV the gluon resummation is enhancing the pure QCD Born amplitude by one order of magnitude. The same enhancement is visible in the QCD corrections to the electroweak cross section. Both amplitudes are rising with decreasing pT , whereas for higher energies the rise of the resummed amplitude will be stronger as the Born one.

16

10 5

γγ → ZT ZT √sγγ = 1 TeV s0 = (1/2 MZ) BFKL Born

2

2 |AQCD(pT )| [√fb GeV]

10 4

10

3

10 2

0

0.5

1

1.5

2

2 2 pT / MZ

2.5

3

3.5

4

Figure 8: The pure QCD part of the amplitude versus the exchanged transverse momenta. The solid line is the BFKL contribution and the dashed line is the QCD Born approximation (two-gluon amplitude).

In Fig. 9 the absolute value of the two helicity changing parts of the amplitude at s = 1TeV are plotted. The most noticable property is that the helicity changing parts are suppressed by more than two orders of magnitude with respect to the helicity conserving part. The ++ → ?? cross section has, apart from its magnitude, properties similar to ++ → ++, for example its slope as a function of pT . The ++ → LL part is also rising when pT is getting smaller but in the region of pT = 1/2 MZ it has an extremum and is steeply going to zero for pT = 0. This behavior comes from eq.(9), where the integrand is antisymmetric in α ? (1 ? α) for l = 0. √

6

Conclusions

In this paper we have presented formulae for the leading QCD corrections to vector boson scattering at small scattering angles. These emerge since, like the photon, all vector bosons ?uctuate into quark-antiquark pairs, which are strongly interacting colour dipoles. The nonforward boson impact factors describe the probability that the boson ?uctuates into a quark-antiquark dipole. First the Born and after the leading logarithmic resummed (BFKL) contributions have been investigated for an arbitrary process of elastic scattering 17

10 2

γγ → Z Z

++ → – – ++ → LL

2 |AQCD(pT )| [√fb GeV]

√sγγ = 1 TeV s0 = (1/2 MZ)

2

10

1

0

0.5

1

1.5

2

2 2 pT / MZ

2.5

3

3.5

4

Figure 9: The same as Fig. 5 for the helicity changing parts of the BFKL contribution.

of vector bosons. In particular the process γγ → ZZ has been considered in detail. Both, the electroweak and the QCD parts of the amplitude have been calculated in the high-energy approximation. The results are displayed for the di?erential and integrated cross section, comparing QCD corrections relative to the one loop elektroweak calculation. We discussed some properties of the pure QCD part of the cross section as well, including also nonleading contributions. The general result of the numerical calculation is that at s = O(1 TeV) the QCD corrections a?ect the cross section of the process by a few percents. The corrections are smaller for larger pT values and rise with falling pT . Moreover they rise with the center of mass energy as expected. Since the leading log BFKL contribution ends up with a noticable scale dependence this results should not be regarded as a precision calculation. Therefore, we interpret this result as an indication that, in particular in the small angle region, possible deviations from EW calculations cannot be interpreted as signals of new physics without taking into account the QCD corrections. The same reasoning is valid for the integrated cross sections, since the biggest contribution to the cross section is coming from scattering at small angles. Acknowledgements We wish to thank J. Bartels for many helpful discussions. Furthermore we wish to thank G. Chachamis for giving us data regarding exact EW calculations, and S. Dittmaier for helpful discussions. K. P. is supported by the Graduiertenkolleg ”Zuk¨ unftige Entwicklungen in der Teilchenphysik”. G.P. Vacca thanks J. Bartels and the II. Institut f¨ ur Theoretische Physik for the warm hospitality. 18

7

Appendix

In this appendix we present further formulas for the impact factor calculated in chapter 4.2. In the case of massive fermions eq.(34) has the form: ΦT T (q, k, l) = αs × ? p2 a

(++) 2?1 Nc

f

1 Cf √ π log log

1 0

dα [α2 + (1 ? α)2 ] ×

?m+

2 m2 f

2 ?a

2 p2 b ? m + 2 mf 2 ?b

m2 f ? ?a + m2 f ?b

log log

?m ? p2 ?m + p2 a ? ?a a ? ?a ? log 2 +? ?m ? p2 + ? ? m + p a a a a

?m ? p2 ?m + p2 b ? ?b b ? ?b ? log 2 ?m ? pb + ?b ?m + p2 b + ?b

?m ? p2 ?m + p2 a ? ?a a ? ?a ? log 2+? ?m ? p2 ? m + p + ? a a a a

?m ? p2 ?m + p2 b ? ?b b ? ?b ? log 2 2 ?m ? pb + ?b ?m + pb + ?b

,

(42)

while eq.(38) becomes:

(i) 2?1 Nc

ΦT L (q, k, l) = αs ?i · pb 2 p2 b ? ?i · pa 2 p2 a p2 b ?m ?b

f

2 Cf √ MZ π

1 0

dα α(1 ? α)(1 ? 2α) ×

(43)

m2 p2 ?m ? p2 ?m + p2 f ?m a?m a ? ?a a ? ?a log ? log ? log 2 2 ?a ?m ? pa + ?a ?m + pa + ?a m2 f ?a = ?b =

2 2 2 (p2 a ? m) + 4mf pa , 2 2 2 (p2 b ? m) + 4mf pb ,

m2 ?m ? p2 ?m + p2 f ?m b ? ?b b ? ?b log ? log ? log ?m ? p2 ?m + p2 m2 b + ?b b + ?b f

with

(44)

2 and p2 a , pb and m was already de?ned in section 4.2. Taking the limit mf → 0 one obtains eq.(34) and eq.(38) respectively. This formulas were used for the top quark contribution, where the remaining α integration was done numerically. In the case of massless fermions for the helicity conserving part the α integration may be performed analytically, giving the result:

ΦT T (q, k, l) = αs

(++)

2?1 Nc

f

1 2 Cf √ 2 x1 ? x2 1 + x2 ? x2 + π

3 2 3 + 2 ? 3 x1 + 3 x2 1 ? 2 x1 log(1 ? x1 ) + 3 x1 ? 3 x1 + 2 x1 log(?x1 ) + 3 2 3 + 2 ? 3 x2 + 3 x2 2 ? 2 x2 log(1 ? x2 ) + 3 x2 ? 3 x2 + 2 x2 log(?x2 ) +

+

d 1 d ?3 a2 + 3 a d + 2 d2 log( )+ 3 3a a 19

d + (a + d) 2 a d ? 2 a2 + a d + 2 d2 log(1 + ) a with

2 a = 4 l2 + MZ + i? 2 b = ?4(l2 + kl) ? MZ ? i?

,

(45)

c = (k + l)2

2 d = ?MZ ? i? √ ?b ± b2 ? 4ac . x1,2 = 2a

(46)

References

[1] M. Baillargeon, G. Belanger and F. Boudjema, hep-ph/9405359. [2] B.W. Lee, C. Quigg and H.B. Thacker, Phys. Rev. D 16 (1977) 1519. [3] G. B? elanger and F. Boudjema, Phys. Lett. B288 (1992) 210; K. A. Peterson and S. ? Godfrey, hep-ph/9302297; O. J. P. Eobli, M. C. Gonzal? ez-Garc? ?a and S. F. Novaes, ? Nucl. Phys. B411 (1994) 381, hep-ph/9306306; O. J. P. Eobli, M. B. Magro, P. G. Mercandante and S. F. Novaes, Phys.Rev. D52 (1995) 15, hep-ph/9503432; S. Godfrey, hep-ph/9505252; G. J. Gounaris, F. M. Renard, Z. Phys. C69 (1996) 513, hep-ph/9505429; G. J. Gounaris, J. Layssac, F. M. Renard, Z. Phys. C69 (1996) 505, hep-ph/9505430; M. Baillargeon, G. B? elanger and F. Boudjema, Nucl.Phys. B500 (1997) 224, hep-ph/9701372; E. Boos. et. al., hep-ph/0103090. [4] F. Boudjema, in Physics and Experiments with Linear e+ e? Colliders, Waikoloa, Hawaii, 1993, Ed. F.A. Harris et al., Word Sci., vol.II, p. 712. [5] G.J. Gounaris, J. Layssac, P.I. Porfyriadis and F.M. Renard, Eur. Phys. J. C 13 (2000) 79. [6] J.L. Hewett and F.J. Petriello, Phys. Rev. D 64 (2001) 095017. [7] T. G. Rizzo, Phys. Rev. D60 (1999) 115010, hep-ph/9904380; K. Cheung, Phys. Rev. D61 (2000) 015005, hep-ph/9904266; H. Davoudiasl, Phys. Rev. D60 (1999) 084022, hep-ph/9904425; Phys. Rev. D61 (2000) 044018, hep-ph/9907347. [8] A. Denner, S. Dittmaier and T. Hahn, Phys. Rev. D 56 (1997) 117; A. Denner and T. Hahn, Nucl. Phys. B 525 (1998) 27; A. Denner, S. Dittmaier and R. Schuster, Nucl. Phys. B 452 (1995) 80. [9] E.A. Kuraev, L.N. Lipatov, V.S.Fadin, Sov. Phys. JETP 44 (1976) 443, ibid. 45 (1977) 199, Ya.Ya. Balitskii, L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. [10] G. Jikia, Nucl. Phys. B 405 (1993) 24. [11] J. Bartels, K. Golec-Biernat and K. Peters, Acta Phys. Polon. B 34 (2003) 3051. [12] S.J. Brodsky, F. Hauptmann and D.E. Soper, Phys. Rev. D 56 (1997) 6957. 20

[13] H. Cheng and T.T. Wu, Phys. Rev. 182 (1969) 1852, 1868, 1873, 1899. [14] V. S. Fadin and L. N. Lipatov, Phys. Lett B429 (1998) 127, and references therein; M. Ciafaloni and G. Camici, Phys. Lett. B430 (1998) 430. [15] J. Bartels, M. A. Braun, D. Colferai and G. P. Vacca, Eur.Phys.J.C 20 (2001) 323 [arXiv:hep-ph/0102221]. [16] L.N. Lipatov, Pomeron in quantum chromodynamics, in “Perturbative QCD”, pp. 411-489, ed. A. H. Mueller, World Scienti?c, Singapore, 1989; Phys. Rep. 286 (1997) 131. [17] J. Bartels, M. G. Ryskin and G. P. Vacca, Eur. Phys. J. C 27 (2003) 101 [arXiv:hep-ph/0207173].

21

赞助商链接

- NLO QCD corrections to tri-boson production
- QCD Corrections and Non-standard Three Vector Boson Couplings in $W^+W^-$ Production at Had
- One-loop Electroweak and QCD corrections to vector boson scattering into top pairs and appl
- QCD CORRECTIONS TO HIGGS BOSON DECAYS
- QCD corrections to tri-boson production

更多相关文章：
更多相关标签：