Full text

Translate

Turn on search term navigation

http://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-015-3625-0&domain=pdf

Web End = http://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-015-3625-0&domain=pdf

Web End = Eur. Phys. J. C (2015) 75:392DOI 10.1140/epjc/s10052-015-3625-0

Regular Article - Theoretical Physics

http://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-015-3625-0&domain=pdf

Web End = http://crossmark.crossref.org/dialog/?doi=10.1140/epjc/s10052-015-3625-0&domain=pdf

Web End = Double vector meson production in the International Linear Collider

F. Carvalho1, V. P. Gonalves2, B. D. Moreira3, F. S. Navarra3,a

1 Departamento de Cincias Exatas e da Terra, Universidade Federal de So Paulo, Campus Diadema, Rua Prof. Artur Riedel, 275, Jd. Eldorado, Diadema, SP 09972-270, Brazil

2 High and Medium Energy Group, Instituto de Fsica e Matemtica, Universidade Federal de Pelotas, Caixa Postal 354, Pelotas, RS 96010-900, Brazil

3 Instituto de Fsica, Universidade de So Paulo, C.P. 66318, So Paulo, SP 05315-970, Brazil

Received: 26 April 2015 / Accepted: 17 August 2015 The Author(s) 2015. This article is published with open access at Springerlink.com

Abstract In this paper we study double vector meson production in interactions at high energies and estimate, using the color dipole picture, the main observables which can be probed at the International Linear Collider (ILC). The total (Q21) + (Q22) V1 + V2 cross sections for Vi = ,

, J/, and are computed and the energy and virtuality dependencies are studied in detail. Our results demonstrate that the experimental analysis of this process is feasible at the ILC and it can be useful to constrain the QCD dynamics at high energies.

1 Introduction

There is an increasing interest in the construction of a high energy electronpositron collider [13]. The primary goal of this new facility will be to carry out precision measurements of electroweak physics, including the Higgs boson properties. An important byproduct of this program will be the study of high energy photonphoton collisions [4] and the continuation, at energies one order of magnitude higher, of the measurements performed at CERN-LEP, almost 15 years ago. Photonphoton collisions are a very clean laboratory for the theory of strong interactionsquantum chromodynamics (QCD)where we can test details of the QCD dynamics at high energies, such as the evolution both in virtuality (Q2) and in energy (1/x) (for a review see, e.g. Ref. [5]). It has motivated the development of a large number of phenomenological studies in the last two decades [638]. In particular, several authors have discussed the possibility of use the scattering of two off-shell photons at high energy in e+ e colliders as a probe of the parton saturation effects in the QCD dynamics, which are predicted to be present in

a e-mail: mailto:[email protected]

Web End [email protected]

the high energy regime [3943]. Although the experimental results on several inclusive and diffractive observables measured in ep scattering at HERA and hadronhadron collisions at RHIC and LHC suggest that these effects are already observed in the energy regime probed by current colliders, these observations still need further conrmation.

The state-of-art framework to treat QCD at high energies is the color glass condensate (CGC) formalism [4452], which predicts the saturation of the growth of parton distributions, with the evolution with the energy being described by an innite hierarchy of coupled equations for the correlators of Wilson linesthe Balitsky-JIMWLK hierarchy (for recent reviews see [3943]). In the mean eld approximation, this set of equations can be approximated by the Balitsky Kovchegov (BK) equation [5355]. As emphasized in Ref. [23], in general, the applications of the CGC formalism to scattering problems require an asymmetric frame, in which the projectile has a simple structure and the evolution occurs in the target wave function, as it is the case in deep inelastic scattering. Therefore the extension of the BK equation to the calculation of the scattering cross section is not a trivial task. In Ref. [23] we have discussed this generalization in order to use the solution of the BK equation as input of our calculations of the total cross sections and photon structure functions, which were compared with the LEP data.

In particular, in Ref. [23] we have improved the treatment of the dipoledipole cross section, which is the main ingredient of the description of the interactions in the dipole picture. Differently from previous phenomenological studies, which disregarded the impact parameter dependence, we have proposed an educated guess for this dependence and demonstrated that the LEP data can be described in this approach. The high energy behavior of the observables predicted in Ref. [23] is largely different from those obtained in previous

123

392 Page 2 of 11 Eur. Phys. J. C (2015) 75:392

studies. This conclusion motivates us to review the analysis of other observables which could be measured at the ILC. One promising observable is double vector meson production in collisions, which has attracted the attention of several theoretical groups in the last years, with the cross section being estimated in different theoretical frameworks [2438], as, for instance, the solution of the BFKL equation and impact factors at leading and next-to-leading orders. In this paper we will estimate the total (Q21)+ (Q22) V1 + V2 cross sec

tions for Vi = , , J/, and considering the improved

treatment of the dipoledipole cross section and the energy and virtuality dependencies of the total cross sections will be analyzed in detail. Our analysis is strongly motivated by the fact that our knowledge about vector meson wave functions has improved considerably over the last years with the progress of phenomenological studies of vector meson production at HERA. As a consequence the main ingredients of our calculations are constrained by LEP and HERA data and hence our predictions for the ILC energies have only one free parameterthe slope parameter BV1V2which determines the t-dependence of the cross sections. The magnitude of this parameter for different combinations of vector mesons is still an open issue that deserves more detailed studies.

This paper is organized as follows. In the next section we present a brief review of the formalism, discussing in more detail the vector meson wave functions and the dipoledipole cross section, which are the main inputs of our calculations. In Sect. 3 we present our predictions for the production of different combinations of vector mesons. In particular, the dependencies of the cross sections on the energy and photon virtualities are analyzed in detail. Finally, in Sect. 4 we present our summary.

2 Double vector meson production

2.1 The cross section

Let us review the main formulas of vector meson production in the color dipole picture (for more details see, e.g. Ref. [34]). The relevant scattering process is V1 V2,

where Vi stands for both light and heavy vector mesons. At high energies, this scattering can be seen as a succession in time of three factorizable subprocesses (see Fig. 1): (i) the photons uctuate into quarkantiquark pairs (the dipoles), (ii) these color dipoles interact and, (iii) the pairs convert into the vector meson nal states. Using as kinematic variables the c.m.s. energy squared s = W2 = (p + q)2, where p

and q are the photon momenta, the photon virtualities squared given by Q21 = q2 and Q22 = p2, and t, the squared

momentum transfer, the total cross section for double vector meson production is given by

Fig. 1 Double vector meson production in interactions at high energies in the color dipole picture

( V1 V2) = dt

d( V1 V2) dt

1 BV1 V2

d( V1 V2) dt

tmin=0

m A(s, t = 0)]2

16 BV1 V2 , (1)

where we have approximated the t-dependence of the differential cross section by an exponential with BV1 V2 being the slope parameter. The imaginary part of the amplitude at zero momentum transfer A(s, t = 0) reads

Im A ( V1 V2)

h, h

n,n

= [I

dz1 d2r1 h,h(z1, r1, Q21) V1h,h(z1, r1)

dz2 d2r2 n,n(z2, r2, Q22) V2n,n(z2, r2)

dd(r1, r2, Y ), (2)

where and Vi are the light-cone wave functions of the photon and vector meson, respectively. The quark and anti-quark helicities are labeled by h, h, n, and n and reference

to the meson and photon helicities are implicitly understood. The variable r1 denes the relative transverse separation of the pair (dipole) and z1 (1 z1) is the longitudinal momen

tum fraction of the quark (antiquark). Similar denitions are valid for r2 and z2. The variable Y is the rapidity and it will be dened later. The basic blocks are the photon wave function, , the meson wave function, VT, L , and the dipoledipole cross section, dd.

123

Eur. Phys. J. C (2015) 75:392 Page 3 of 11 392

2.2 Wave functions

In the dipole formalism, the light-cone wave functions h,h(z, r) in the mixed representation (r, z) are obtained

through two dimensional Fourier transform of the momentum space light-cone wave functions h,h(z, k) [56]. This sub

ject has been intensely discussed in several references (see e.g. Refs. [5761]). In what follows we present, for completeness, some of the main formulas. The normalized light-cone wave functions for longitudinally (L) and transversely (T ) polarized photons are given by

Lh,h(z, r) =

Nc4 h,h e e

where the effective charge f = 1/3, 2/3, 1/3, or 1/2,

for , J/, , or mesons, respectively. The assumption that the quantum numbers of the meson are saturated by the quarkantiquark pair and that the possible contributions of gluon or sea-quark states to the wave function may be neglected, allows the normalization of the vector meson wave functions to unity. The normalization conditions for the scalar parts of the wave functions are then

1 =

Nc 2

dzz2(1 z)2

d2r

m2f

0 2T + z2 +

(1 z)2 (r

T )2

, (9)

f 2z(1 z) Q

K0(r)

2 , (3)

T(=)h,h(z, r)

1 =

Nc 2

0 dz d2r MV L +

m2f 2r

MV z(1 z)

Nc2 e e f [ieir (zh,h (1 z)h,h)r

+m f h,h]

K0(r)

2 , (4)

where 2 = z(1 z)Q2 + m2f . The quark mass m f plays

the role of a regulator when the photoproduction regime is reached. The electric charge of the quark of avor f is given by e e f .

One simple way to model the vector meson wave function is to assume, following Refs. [57,5961], that the vector meson is a quarkantiquark state and that the spin and polarization structure is the same as in the photon case. The transversely polarized vector meson wave function is then given by

Vhh,=1

(r, z)

= s 2N

(10)

Another important constraint on the vector meson wave functions is obtained from the decay width. It is commonly assumed that the decay width can be described in a factorized way: the perturbative matrix element q q l+l fac

torizes out from the details of the wave function, which contributes only through its properties at the origin. The decay widths are then given by

fV,T = f

Nc 2 MV

dz z2(1 z)2

m2f

z2 + (1 z)2 2r

T (r, z)|r=0, (11)

fV,L = f

0 dz MV +

m2f 2r

MV z(1 z)

c 1z(1 z) {

ieir [zh,h,

(1 z)h,h,]r + mf h,h,} T (r, z). (5)

and the longitudinally polarized wave function is given by

Vhh,=0

(r, z)

= N

L(r, z)|r=0.

(12)

The coupling of the meson to the electromagnetic current, fV , is obtained from the measured electronic decay width by

L(r, z), (6)

where 2r (1/r)r + 2r and MV is the meson mass. The

overlaps between the photon and the vector meson wave functions read then

( V )T = f e

Ncz(1 z)

Ve+e =

42em f 2V

3MV . (13)

We need now to specify the scalar parts of the wave functions, T,L(r, z). Dosch, Gousset, Kulzinger and Pirner (DGKP)

[57] made the assumption that the longitudinal momentum fraction z uctuates independently of the transverse quark momentum k, where k is the Fourier conjugate variable to the dipole vector r. In the DGKP model one chooses = 0

in Eqs. (6), (8), (10), and (12). The DGKP model was further simplied by Kowalski and Teaney [60,61], who assumed that the z dependence of the wave function for the longitudinally polarized meson is given by the short-distance limit of z(1 z). For the transversely polarized meson they set

T (r, z) [z(1 z)]2 in order to suppress the contribution

from the end-points (z 0, 1). This leads to the Gauss-LC

c h,h MV +

m2f 2r

MV z(1 z)

m2f K0( r)T (r, z)

z2 + (1 z)2

K1( r)rT (r, z)

, (7)

( V )L = f e

2Qz(1 z) K0( r)

MV L(r, z) +

m2f 2r

MV z(1 z)

L(r, z)

, (8)

123

392 Page 4 of 11 Eur. Phys. J. C (2015) 75:392

Table 1 Parameters of the Gauss-LC vector meson wave functions

Meson MV (GeV) fV m f (GeV) NT R2T (GeV2) NL R2L (GeV2)

9.460 0.236 4.5 0.67 2.16 0.47 1.01

J/ 3.097 0.274 1.4 1.23 6.5 0.83 3.0

1.019 0.076 0.14 4.75 16.0 1.41 9.7

0.776 0.156 0.14 4.47 21.9 1.79 10.4

[60,61] wave functions given by

T (r, z) = NT [z(1 z)]2 exp(r2/2R2T ), (14) L(r, z) = NLz(1 z) exp(r2/2R2L). (15)

The values of the constants NT,L and RT,L in Eqs. (14) and (15), determined by requiring the correct normalization and by the condition fV = fV,T = fV,L, are given in Table 1. It

is important to emphasize that this model allows to describe the HERA data and the recent LHC data for the exclusive vector meson photoproduction in hadronhadron collisions (see, e.g. Refs. [6264]).

2.3 The dipoledipole scattering cross section

At lowest order, the dipoledipole interaction can be described by the two-gluon exchange between the dipoles, with the resulting cross section being energy independent (see, e.g. Ref. [65]). Taking into account the leading corrections associated to terms log(1/x), as described by the

BFKL equation, implies a power-law energy behavior for the cross section, which violates the unitarity at high energies. These unitarity corrections were addressed in Ref. [66], considering the color dipole picture and independent multiple scatterings between the dipoles, and in Refs. [67,68] considering the Color Glass Condensate formalism. In the eikonal approximation the dipoledipole cross section can be expressed as follows:

dd(r1, r2, Y ) = 2 d2b N (r1, r2, b, Y ) (16)

where N (r1, r2, b, Y ) is the scattering amplitude for the

two dipoles with transverse sizes r1 and r2, relative impact parameter b and rapidity separation Y . How to write N for

the interaction of two dipoles of similar sizes is still an open question (see, e.g. Ref. [69]). In a rst approximation, it is useful to express N in terms of the solution of the BK equa

tion (obtained disregarding the b dependence), which has been derived considering an asymmetric frame where the projectile has a simple structure and the evolution occurs in the target wave function. A shortcoming of this approach, used in a previous analysis of double vector meson production, is that although the unitarity of the S-matrix (N 1)

is respected by the solution of the BK equation, the associated dipoledipole cross section can still rise indenitely

with the energy, even after the black disk limit (N = 1)

has been reached at central impact parameters, due to the non-locality of the evolution. In Ref. [23] we have proposed a more elaborated model for the impact parameter dependence in order to obtain more realistic predictions for the dipoledipole cross section. Basically, we assumed that only the range b < R, where R = Max(r1, r2), contributes to

the dipoledipole cross section, i.e., we assumed that N is

negligibly small when the dipoles have no overlap with each other (b > R). Therefore the dipoledipole cross section can be expressed as follows [23]:

dd(r1, r2, Y ) = 2 N(r, Y )

0 d2b = 2 R2N(r, Y ), (17) where N(r, Y ) is the dipole scattering amplitude. The explicit form of dd reads

dd(r1, r2, Y ) = 2r21N(r2, Y2) (r1 r2)

+2r22N(r1, Y1) (r2 r1), (18) where Yi = ln(1/xi) and

xi =

Q2i + 4m2f

W2 + Q2i

. (19)

Following Ref. [23] this model for dd will be called Model 2. As input for this model we will use two forward dipole scattering amplitudes. The rst one is the solution of the BK equation obtained in Ref. [70], which we call rcBK. The second one is the phenomenological model of the forward dipole scattering N(r, Y ) proposed in Ref. [71] and updated in [72], which was constructed so as to reproduce two limits of the LO BK equation analytically under control: the solution of the BFKL equation for small dipole sizes, r 1/Qs(x), and the LevinTuchin law for larger ones,

r 1/Qs(x). In the updated version of this parametrization

[72], the free parameters were obtained by tting the new H1 and ZEUS data. In this parametrization the forward dipole scattering amplitude is given by

N(r, Y ) =

2 s+ ln(2/r Qs) Y , for r Qs(x) 2, 1 expa ln2 (b r Qs), for r Qs(x) > 2,

r Qs 2

(20)

where a and b are determined by continuity conditions at r Qs(x) = 2, s = 0.6194, = 9.9, = 0.2545, Q20 =

123

Eur. Phys. J. C (2015) 75:392 Page 5 of 11 392

1.0 GeV2, x0 = 0.2131 104, and N0 = 0.7. Hereafter,

we shall call the model above IIM-S. The rst line from Eq. (20) describes the linear regime whereas the second one includes saturation effects. One of the main motivations to use this model in our analysis is that it allows one to estimate the magnitude of the saturation effects, by the comparison between the predictions of the full model with those obtained considering only the linear term. As demonstrated in Ref. [23], using this model we can describe the LEP data for the total cross sections and photon structure functions.

Finally, for the sake of comparison with a previous analysis, in what follows we also will present the predictions obtained using the phenomenological model for the dipole dipole cross section proposed in [22], which disregards the impact parameter dependence and, consequently, presents the shortcoming discussed above. The inclusion of these predictions in our analysis allows us to estimate the theoretical uncertainty present in ILC predictions, as well as to make comparisons with existing results in the literature. The dipoledipole cross section proposed in Ref. [22] is the following:

dda,b(r1, r2, Y ) = a,b0 N(r1, r2, Y ) (21)

with a,b0 = (2/3)0, where 0 is a free parameter in the

saturation model considered, xed by tting the DIS HERA data. In the above equation N(r1, r2, Y ) = N(reff, Y =

ln(1/ xab)), where

r2eff =

r21r22

r21 + r22

and xab =

Q21 + Q22 + 4m2a + 4m2b W2 + Q21 + Q22

(22)

Keeping the notation introduced in [23], this model will be called Model 1. As input for the calculations, we will use the same two dipole amplitudes: rcBK and IIM-S.

3 Results

In what follows, as in Ref. [23], we will denote the predictions obtained using the dipoledipole cross section given by Eq. (21) by model 1 and those using Eq. (18) as input by model2. Moreover, we will consider the rcBK and IIM-S models for the scattering amplitude N(r, Y ). The parameters of our calculations are the same used in Ref. [23] and this implies that our model gives a good description of the LEP data. As the value of the slope BV1V2 for the different combinations of vector mesons in the nal state is still poorly known, we will, in almost all cases, present our predictions for the product BV1V2 ( V1V2), which can be estimated without

free parameters, since all parameters are constrained by the LEP and HERA data.

In Fig. 2 we present our predictions for the energy dependence of the product BV1V2 ( V1V2) assuming

V1 = V2 (Vi = , , J/, ) and considering Q21 = Q22 =

0. In this case only the transverse photon polarizations contribute to the total cross sections. It is important to emphasize

104

103

IIM-S Model 1 IIM-S Model 2 rcBK Model 1 rcBK Model 2

2 ]

B V 1 V 2[nb/GeV

102

101

10 0 500 1000 1500 2000

100

10-4

10-5

2 ]

10-1

B V 1 V 2[nb/GeV

10-6

10-2

J/J/

10-7

Upsilon Upsilon

10-3

10-8

0 500 1000 1500 2000W(GeV)

Fig. 2 Energy dependence of the product BV1V2 [ (Q21) (Q22) V1V2] assuming V1 = V2 (Vi = , , J/, ) and considering Q21 =

Q22 = 0

123

392 Page 6 of 11 Eur. Phys. J. C (2015) 75:392

102

IIM-S Model 1 IIM-S Model 2 rcBK Model 1 rcBK Model 2

2 ]

101

B V 1 V 2[nb/GeV

100

-1

10 0 500 1000 1500 2000

10-1

2 ]

10-2

B V 1 V 2[nb/GeV

10-3

10-4

10-5

0 500 1000 1500 2000W(GeV)

Fig. 3 Energy dependence of the product BV1V2 [ (Q21) (Q22) V1V2] assuming V1 = V2 (V1V2 = J/, J/, J/ , ) and

considering Q21 = Q22 = 0

that the color dipole picture allows us to treat simultaneously double production by real photons, which is a typical soft process, and double production, which is the ideal laboratory to study the basic example of a hard process at high energies: the oniumonium scattering. Moreover, it allows us to study the transition between these two regimes, where we expect to see nonlinear (saturation) effects in the QCD dynamics. In our calculations we consider the two different models for the dipoledipole cross section as well as the two models for the forward dipole scattering amplitude. We can observe that the main distinction is associated to the choice of the dipoledipole cross sections. The predictions obtained using model 2 are always smaller than those from model 1.Consequently, assuming that model 2 is more adequate to describe the dipoledipole interaction, we see that the previous estimates of the double vector production have overestimated the magnitude of the total cross sections. This behavior was expected from our previous results for the total

cross section [23]. We also see that the difference between the predictions increases with the quark masses, going from a factor 4, in the case, to almost two orders of magnitude in the case of production. This suggests that the experimental analysis of double vector production at ILC can, in principle, constrain the model for the dipoledipole interaction. Moreover, in the case of model 1, the IIM-S and rcBK predictions are almost identical for all combinations of vector mesons in the nal state. In model 2 the IIM-S predictions are smaller than the rcBK ones for light vector meson production

and larger for heavy vector meson production. Such behavior is directly associated to the distinct transition between small and large dipoles predicted by these two models for the forward dipole scattering amplitude (see Fig. 2 in Ref. [23]). As expected, we nd that our predictions are strongly dependent on the quark mass, with the cross sections being smaller for the production of heavier vector mesons. Similar conclusions are obtained in the analysis shown in Fig. 3, where we present our predictions for the energy dependence of the product BV1V2 [ (Q21) (Q22) V1V2] assuming

V1 = V2 (V1V2 = J/, J/, J/ , ) and consider

ing Q21 = Q22 = 0.

In Fig. 4 we present our predictions for the dependence on the photon virtualities Q21 = Q22 = Q2 of the product

BV1V2 [ (Q21) (Q22) V1V2] for different combina

tions of vector mesons in the nal state and xed center-of-mass energy (W = 500 GeV). In this case we take into

account the contributions of the transverse and longitudinal photon polarizations, i.e.:

[ (Q21) (Q22) V1V2]= [ L(Q21) L(Q22) V1V2]

+[ L(Q21) T(Q22) V1V2] +[ T(Q21) L(Q22) V1V2] +[ T(Q21) T(Q22) V1V2]. (23) For Q2 = 0 we have two hard scales present in the pro

cess: the mass of the quarks (vector mesons) and the photon

123

Eur. Phys. J. C (2015) 75:392 Page 7 of 11 392

102

101

2 ]

100

B V 1 V 2[nb/GeV

10-1

10-2

W = 500 GeV

10-3

W = 500 GeV

10-3

10-4

5 10 15 20

100

10-1

10-2

Upsilon-Upsilon

W = 500 GeV

IIM-S Model 1 IIM-S Model 2 rcBK Model 1 rcBK Model 2

2 ]

B V 1 V 2[nb/GeV

10-3

10-4

10-5

10-6

J/J/

W = 500 GeV

10-6

10-7

10-8

5 10 15 20Q2(GeV2)

Fig. 4 Dependence on the photon virtualities Q21 = Q22 = Q2 of the product BV1V2 [ (Q21) (Q22) V1V2] assuming V1 = V2 (Vi =

, , J/, ) for a xed center-of-mass energy (W = 500 GeV)

102

101

W = 500 GeV

101

W = 500 GeV

2 ]

B V 1 V 2[nb/GeV

100

10-1

10-2

10-3

10-4

5 10 15 20

100

10-1

W = 500 GeV

10-1

W = 500 GeV

IIM-S Model 1 IIM-S Model 2 rcBK Model 1 rcBK Model 2

2 ]

10-2

B V 1 V 2[nb/GeV

10-3

10-4

10-5

10-6

10-7

10-8

5 10 15 20Q2(GeV2)

Fig. 5 Dependence on the photon virtualities Q21 = Q22 = Q2 of the product BV1V2 [ (Q21) (Q22) V1V2] assuming V1 = V2 (V1V2 =

J/, J/, J/ , ) for a xed center-of-mass energy (W = 500 GeV).

virtualities. For the double light vector meson (, ) production, the dominant scale is the photon virtuality. In this case our predictions strongly decrease with Q2. On the other

hand, for the double production, our predictions are almost Q2-independent in the range considered, since the dominant scale that denes the size of the two interacting dipoles

123

392 Page 8 of 11 Eur. Phys. J. C (2015) 75:392

Model 1

Q2 = 0

Model 1Q2 = 20 GeV2

0 100 200 300 400 500 600 700 800 900 1000

Model 2

Q2 = 0

J/ J/J/

Model 2Q2 = 20 GeV2

J/ J/J/

0 100 200 300 400 500 600 700 800 900 1000

W(GeV)

0 100 200 300 400 500 600 700 800 900 1000

W(GeV)

(a) (b)

Fig. 6 Energy dependence of the normalized cross sections (see text) for different nal states and different values of Q21 = Q22 = Q2. a Q2 = 0

and b Q2 = 20 GeV2

is the bottom quark mass. In contrast, for the double J/ production, the characteristic dipole sizes are determined at small Q2 by the charm quark mass and at medium Q2 by the photon virtualities. Consequently, we observe a mild Q2 dependence in the corresponding predictions. Moreover, we observe that the difference between model 1 and model 2 predictions increases at larger Q2 and for heavier vector mesons.In Fig. 5 we present our predictions for the production of two different vector mesons, which are similar to those observed in the production of identical vector mesons. Basically, the Q2 dependence is reduced for larger values of the sum of the masses of the vector mesons in the nal state.

In order to illustrate how the energy behavior depends on the masses of the nal state mesons, on the photon virtualities Q21 = Q22 = Q2 and on the choice of the model for the

dipoledipole cross section, in Fig. 6 we present our predictions for the normalized cross sections. The different cross sections were all normalized to the unity at W = 100 GeV

to better exhibit the different trends. For Q2 = 0 we observe

a clear transition between the soft and hard regimes, with the growth with the energy being faster for heavier mesons in the nal state. Moreover, we nd that model 2 predicts a smaller slope than model 1. For Q2 = 20 GeV2, a similar behavior

is observed, but in this case already for the production we see a steep rise of the cross section with the energy, which is directly associated to the presence of the hard scale Q2.

In certain cases, where the slope parameters are phenomenologically known, it is possible to make denite predictions. In Fig. 7 we show the cross sections calculated with models 1 and 2 as a function of the energy W with the proper slope coefcients, taken from [34]: B = 10

GeV2, B = 5 GeV2, and B = 0.44 GeV2. Our

predictions with model 1 are similar to those obtained in

[34], with small differences mainly associated to the different forward dipole scattering amplitude and to the treatment of the vector meson wave functions. In contrast, with model 2, we predict that at W = 1 TeV the cross sections

are ( ) 15 nb, ( J/ ) 1.2 nb

and ( J/ J/ ) 0.25 nb, which are a factor 4 smaller than previous estimates in the literature [27,28]

obtained using the color dipole picture. In our calculation we have been using a formalism in which the real part of the amplitude is neglected. This might not be a good approximation for heavy mesons. Nevertheless, since one of our primary goals in this work is to investigate the predictions of our dipoledipole scattering model, we keep using this approximation to compare our numbers with those obtained in [27,28], where the same approximation was made.

A similar analysis can be performed for the interaction of virtual photons. In this case we can compare our predictions with those obtained in Refs. [3133] for the double production. In [3133] the cross section

was obtained using the kT -factorization formalism, with the scattering amplitude being given in terms of the convolution between the impact factors, BFKL kernel and the corresponding leading twist distribution amplitude, which describes the hadronization into the nal state mesons. The authors have considered the BFKL evolution in the leading logarithm (LL) approximation as well as the renormalization group improved BFKL kernel, which gives good agreement with the next-to-leading logarithm (NLL) evolution derived in Refs. [35,36]. In comparison with [3133], our results are a factor 103

smaller than the LL BFKL predictions and similar to the NLL BFKL one.

Although saturation effects of QCD dynamics are expected to take over at higher energies, this change of

123

Eur. Phys. J. C (2015) 75:392 Page 9 of 11 392

Fig. 7 Energy dependence of the V1V2 cross section

for different nal states considering Q21 = Q22 = 0

200

IIM-S Model 1 IIM-S Model 2 rcBK Model 1 rcBK Model 2

J/J/

150

1.5

(nb)

100

0.5

0 0 500 1000 1500 2000

W(GeV)

0 0 500 1000 1500 2000

W(GeV)

0 0 500 1000 1500 2000

W(GeV)

300

IIM-S Model 1 - full IIM-S Model 1 - linear IIM-S Model 2 - full IIM-S Model 2 - linear

250

200

(nb)

150

100

0 0 500 1000 1500 2000

W(GeV)

Fig. 8 Comparison between the linear and full IIM-S predictions for the energy behavior of the cross section

dynamics can manifest itself with different strength for different observables. For this reason, when working with the dipole approach, which is naturally prepared to incorporate nonlinear corrections, it is always interesting to quantify the importance of the saturation effects. In Fig. 8 we show our results for the energy dependence of the cross sec

tion, where we compare the full IIM-S model predictions with those obtained considering the linear regime of this model [rst line in Eq. (20)]. We see that the high energy behavior of the cross section is strongly modied by saturation effects. This conclusion was already obtained in [34] and we see that it remains valid even after updating the dipole cross sections and model for the dipoledipole interaction.

A nal comment is in order. Recently, in Refs. [73,74], the dipole representation of vector meson electroproduction

was studied beyond the standard twist-2 level. The authors of these references have considered the description of hard exclusive processes beyond the leading twist approximation derived in Refs. [75,76] and demonstrated that the helicity amplitudes, expressed in the impact parameter space and computed in the collinear factorization scheme, factorize into the dipole cross section and the wave functions of the virtual photon combined with the rst moments of the meson wave functions parametrized by the distribution amplitudes. In Ref. [74] the predictions of this approach were compared with the HERA ep data and a good description was obtained for Q2 5 GeV2. The extension of this approach to

interactions and the comparison of its predictions with those obtained in this paper is an important issue to be addressed in the future.

4 Summary

The scattering of two off-shell photons at high energy in e+ e colliders is an interesting process to look for parton saturation efffects. In these two-photon reactions, the photon virtualities can be made large enough to ensure the applicability of perturbative methods or can be varied in order to test the transition between the soft and hard regimes of the QCD dynamics. In recent years, a series of studies have discussed in detail the treatment of the total cross section and the exclusive production of different nal states in interactions considering very distinct theoretical approaches. One great motivation for these works is the possibility that in a near future interactions may be investigated at the

123

392 Page 10 of 11 Eur. Phys. J. C (2015) 75:392

International Linear Collider (ILC). In particular, in Ref. [23] we presented a detailed analysis of the cross section at high energies using the color dipole picture and taking into account saturation effects, which are expected to be visible at high energies. In this paper we extended our approach to double vector meson production, improving the previous analysis in three important aspects: (i) the theoretical treatment of the dipoledipole cross section; (ii) the forward scattering amplitude, considering the solution of the running coupling BK equation (which is the state-of-art of the CGC formalism); and (iii) the treatment of the vector meson wave functions. Considering that all parameters of our approach have been xed by tting HERA and LEP data, our predictions for double vector meson production at ILC are parameter free, except for the only unknown parameter: the slope parameter BV1V2, which deserves a more detailed analysis. Our main conclusion is that the improvement of the theoretical framework for double vector meson production in interactions resulted in a reduction of the previously estimated cross sections at ILC energies. However, our results indicate that the experimental analysis at ILC is feasible and may be useful to constrain the QCD dynamics at high energies. As a nal remark we would like to say that understanding vector meson production in collisions is very important not only for the phenomenology of future electronpositron colliders but also for exclusive double vector meson production in hadronhadron collisions, which have been studied at the LHC and that could be further studied in future hadronic colliders.

Acknowledgments This work was partially nanced by the Brazilian funding agencies CNPq, CAPES, FAPERGS and FAPESP.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecomm http://creativecommons.org/licenses/by/4.0/

Web End =ons.org/licenses/by/4.0/ ), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Funded by SCOAP3.

References

1. G. Moortgat-Pick, H. Baer, M. Battaglia, G. Belanger, K. Fujii,J. Kalinowski, S. Heinemeyer, Y. Kiyo et al. http://arxiv.org/abs/1504.01726

Web End =arXiv:1504.01726 [hep-ph]2. H. Baer, T. Barklow, K. Fujii, Y. Gao, A. Hoang et al. http://arxiv.org/abs/1306.6352

Web End =arXiv:1306.6352 [hep-ph]

3. M. Bicer et al., JHEP 1401, 164 (2014) and references therein4. V.M. Budnev, I.F. Ginzburg, G.V. Meledin, V.G. Serbo, Phys. Rept. 15, 181 (1975)

5. R. Nisius, Phys. Rept. 332, 165 (2000)6. I.F. Ginzburg, S.L. Panl, V.G. Serbo, Nucl. Phys. B 284, 685 (1987)

7. J. Bartels, A. De Roeck, H. Lotter, Phys. Lett. B 389, 742 (1996)

8. A. Donnachie, H.G. Dosch, M. Rueter, Phys. Rev. D 59, 074011 (1999)

9. J. Bartels, C. Ewerz, R. Staritzbichler, Phys. Lett. B 492, 56 (2000)10. A. Bialas, W. Czyz, W. Florkowski, Eur. Phys. J. C 2, 683 (1998)11. J. Kwiecinski, L. Motyka, Eur. Phys. J. C 18, 343 (2000)12. N.N. Nikolaev, B.G. Zakharov, V.R. Zoller, JETP 93, 957 (2001)13. S.J. Brodsky, F. Hautmann, D.E. Soper, Phys. Rev. D 56, 6957 (1997)

14. S.J. Brodsky, F. Hautmann, D.E. Soper, Phys. Rev. Lett. 78, 803 (1997)

15. M. Boonekamp, A. De Roeck, C. Royon, S. Wallon, Nucl. Phys. B 555, 540 (1999)

16. S.J. Brodsky, V.S. Fadin, V.T. Kim, L.N. Lipatov, G.B. Pivovarov, Pisma ZHETF 76, 306 (2002) [JETP Lett. 76, 249 (2002)]

17. M. Kozlov, E. Levin, Eur. Phys. J. C 28, 483 (2003)18. V.P. Goncalves, M.V.T. Machado, W.K. Sauter, J. Phys. G 34, 1673 (2007)

19. F. Caporale, D.Y. Ivanov, A. Papa, Eur. Phys. J. C 58, 1 (2008)20. D.Y. Ivanov, B. Murdaca, A. Papa, JHEP 1410, 58 (2014)21. G.A. Chirilli, Y.V. Kovchegov, JHEP 1405, 099 (2014)22. N. Timneanu, J. Kwiecinski, L. Motyka, Eur. Phys. J. C 23, 513 (2002)

23. V.P. Goncalves, M.S. Kugeratski, E.R. Cazaroto, F. Carvalho, F.S. Navarra, Eur. Phys. J. C 71, 1779 (2011)

24. I.F. Ginzburg, S.L. Panl, V.G. Serbo, Nucl. Phys. B 296, 569 (1988)

25. J. Kwiecinski, L. Motyka, Phys. Lett. B 438, 203 (1998)26. C.F. Qiao, Phys. Rev. D 64, 077503 (2001)27. V.P. Goncalves, M.V.T. Machado, Eur. Phys. J. C 28, 71 (2003)28. V.P. Goncalves, M.V.T. Machado, Eur. Phys. J. C 29, 271 (2003)29. V.P. Goncalves, W.K. Sauter, Eur. Phys. J. C 44, 515 (2005)30. V.P. Goncalves, W.K. Sauter, Phys. Rev. D 73, 077502 (2006)31. B. Pire, L. Szymanowski, S. Wallon, Eur. Phys. J. C 44, 545 (2005)32. R. Enberg, B. Pire, L. Szymanowski, S. Wallon, Eur. Phys. J. C 45, 759 (2006) [Erratum-ibid. C 51, 1015 (2007)]

33. M. Segond, L. Szymanowski, S. Wallon, Eur. Phys. J. C 52, 93 (2007)

34. V.P. Goncalves, M.V.T. Machado, Eur. Phys. J. C 49, 675 (2007)35. D.Y. Ivanov, A. Papa, Eur. Phys. J. C 49, 947 (2007)36. D.Y. Ivanov, A. Papa, Nucl. Phys. B 732, 183 (2006)37. M. Klusek, W. Schafer, A. Szczurek, Phys. Lett. B 674, 92 (2009)38. S. Baranov, A. Cisek, M. Klusek-Gawenda, W. Schafer, A. Szczurek, Eur. Phys. J. C 73, 2335 (2013)

39. F. Gelis, Int. J. Mod. Phys. A 28, 1330001 (2013)40. F. Gelis, E. Iancu, J. Jalilian-Marian and R. Venugopalan. http://arxiv.org/abs/1002.0333

Web End =arXiv:1002.0333

41. E. Iancu, R. Venugopalan. http://arxiv.org/abs/hep-ph/0303204

Web End =arXiv:hep-ph/0303204 42. H. Weigert, Prog. Part. Nucl. Phys. 55, 461 (2005)43. J. Jalilian-Marian, Y.V. Kovchegov, Prog. Part. Nucl. Phys. 56, 104 (2006)

44. L.D. McLerran, R. Venugopalan, Phys. Rev. D 49, 2233 (1994)45. E. Iancu, A. Leonidov, L. McLerran, Nucl. Phys. A 692, 583 (2001)46. E. Ferreiro, E. Iancu, A. Leonidov, L. McLerran, Nucl. Phys. A 703, 489 (2002)

47. J. Jalilian-Marian, A. Kovner, L. McLerran, H. Weigert, Phys. Rev. D 55, 5414 (1997)

48. J. Jalilian-Marian, A. Kovner, H. Weigert, Phys. Rev. D 59, 014014 (1999)

49. J. Jalilian-Marian, A. Kovner, H. Weigert, Phys. Rev. D 59, 014015 (1999)

50. J. Jalilian-Marian, A. Kovner, H. Weigert, Phys. Rev. D 59, 034007 (1999)

51. A. Kovner, J. Guilherme Milhano, H. Weigert, Phys. Rev. D 62, 114005 (2000)

52. H. Weigert, Nucl. Phys. A 703, 823 (2002)53. I. Balitsky, Nucl. Phys. B 463, 99 (1996)

123

Eur. Phys. J. C (2015) 75:392 Page 11 of 11 392

54. Y.V. Kovchegov, Phys. Rev. D 60, 034008 (1999)55. Y.V. Kovchegov, Phys. Rev. D 61, 074018 (2000)56. V. Barone, E. Predazzi, High-Energy Particle Diffraction (Springer, Berlin, 2002)

57. H.G. Dosch, T. Gousset, G. Kulzinger, H.J. Pirner, Phys. Rev. D 55, 2602 (1997)

58. S. Munier, A.M. Stasto, A.H. Mueller, Nucl. Phys. B 603, 427 (2001)

59. J.R. Forshaw, R. Sandapen, G. Shaw, Phys. Rev. D 69, 094013 (2004)

60. H. Kowalski, D. Teaney, Phys. Rev. D 68, 114005 (2003)61. H. Kowalski, L. Motyka, G. Watt, Phys. Rev. D 74, 074016 (2006)62. V.P. Goncalves, B.D. Moreira, F.S. Navarra, Phys. Rev. C 90, 015203 (2014)

63. V.P. Goncalves, B.D. Moreira, F.S. Navarra, Phys. Lett. B 742, 172 (2015)

64. N. Armesto, A.H. Rezaeian, Phys. Rev. D 90, 054003 (2014)

65. H. Navelet, S. Wallon, Nucl. Phys. B 522, 237 (1998)66. A.H. Mueller, G.P. Salam, Nucl. Phys. B 475, 293 (1996)67. E. Iancu, A.H. Mueller, Nucl. Phys. A 730, 460 (2004)68. G.P. Salam, Nucl. Phys. B 461, 512 (1996)69. Y.V. Kovchegov, Phys. Rev. D 72, 094009 (2005)70. J.L. Albacete, N. Armesto, J.G. Milhano, C.A. Salgado, Phys. Rev. D 80, 034031 (2009)

71. E. Iancu, K. Itakura, S. Munier, Phys. Lett. B 590, 199 (2004)72. G. Soyez, Phys. Lett. B 655, 32 (2007)73. A. Besse, L. Szymanowski, S. Wallon, Nucl. Phys. B 867, 19 (2013)74. A. Besse, L. Szymanowski, S. Wallon, JHEP 1311, 062 (2013)75. I.V. Anikin, D.Y. Ivanov, B. Pire, L. Szymanowski, S. Wallon, Phys. Lett. B 682, 413 (2010)

76. I.V. Anikin, D.Y. Ivanov, B. Pire, L. Szymanowski, S. Wallon, Nucl. Phys. B 828, 1 (2010)

123

Word count: 7022

Show less

The Author(s) 2015

Abstract

Translate

(ProQuest: ... denotes formulae and/or non-USASCII text omitted; see image)

In this paper we study double vector meson production in ...... interactions at high energies and estimate, using the color dipole picture, the main observables which can be probed at the International Linear Collider (ILC). The total ...... cross sections for ......, ......, ......, and ...... are computed and the energy and virtuality dependencies are studied in detail. Our results demonstrate that the experimental analysis of this process is feasible at the ILC and it can be useful to constrain the QCD dynamics at high energies.

Details

Title

Double vector meson production in the International Linear Collider

Author

Carvalho, F; Gonçalves, V P; Moreira, B D; Navarra, F S

Pages

1-11

Publication year

2015

Publication date

Aug 2015

Publisher

Springer Nature B.V.

ISSN

14346044

e-ISSN

14346052

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1140/epjc/s10052-015-3625-0

ProQuest document ID

1732878420

The Author(s) 2015

Double vector meson production in the International Linear Collider

Jump to:

Full text

Abstract

Details

Suggested sources