Unintegrated gluon distribution from forward polarized $\rho$-electroproduction

We present here some arguments to support our suggestion that data on the helicity structure for the hard exclusive electroproduction of $\rho$ mesons at HERA (and in possible future high-energy electron-proton colliders) provide useful information to constrain the $\kappa$-shape of the unintegrated gluon distribution in the proton.


I. INTRODUCTION
Our ability to find new Physics and understand the dynamics of strong interactions at the LHC strongly relies on getting a more and more precise knowledge of the structure of the proton. In general, the latter is encoded in different types of partonic distribution functions that enter the factorization formalism for the description of the hard processes.
Collinear factorization is the most developed approach to calculate cross sections of inclusive reactions as a power expansion over the hard-scale parameter. A prominent example here is the deep inelastic scattering (DIS) of an electron off a proton. Its cross section, at the leading order in the power expansion over the virtuality Q 2 of the exchanged photon γ * , is factorized as a convolution of a hard cross sections (calculable in perturbation theory) with parton distribution functions (PDFs) of quarks and gluons, q i (ζ, µ F ) and g(ζ, µ F ), that depend on the longitudinal momentum fraction of the proton carried by the parton, ζ, and on the factorization scale µ F , and obey DGLAP evolution equations [1]. At the leading order (LO) of perturbation theory the variable ζ coincides with the Bjorken variable x = Q 2 /(W 2 + Q 2 ), where W 2 is the squared center-of-mass energy of the γ * p system.
The collinear factorization scheme can be also applied to the amplitudes of hard exclusive processes, where the nonperturbative part is factorized in generalized parton distributions [2,3].
At high energy, W Q Λ QCD , the application of collinear factorization is limited because the perturbative expansion includes in this kinematics large logarithms of the energy that have to be resumed. Such a resummation is incorporated in the κ-factorization.
The scattering amplitudes are basically written as a convolution of the unintegrated gluon distribution (UGD) in the proton with the impact factor (IF) that depends on the considered process. In the DIS case the γ * → γ * IF is calculated fully in perturbation theory.
The UGD is a nonperturbative quantity, function of x and κ, where the latter represents the gluon momentum transverse to the direction of the proton and is the Fourier-conjugate variable of the transverse separation r d of the color dipole into which the virtual photon splits. Therefore small values of r d mean large values of κ and vice versa. The UGD, in its original definition, obeys the BFKL [4] evolution equation in the x variable. Differently from collinear PDFs, the UGD is not well known and several types of models for it do exist, which lead to very different shapes in the (x, κ)-plane (see, for instance, Ref. [5]).
The aim of this paper is to present our arguments that HERA data on polarization observables in vector meson (VM) electroproduction can be used to constrain the κ-dependence of the UGD in the HERA energy range. In particular, we will focus our attention on the ratio of the two dominant amplitudes for the polarized electroproduction of ρ mesons, i.e. the longitudinal VM production from longitudinally polarized virtual photons and the transverse VM production from transversely polarized virtual photons.
The H1 and ZEUS collaborations performed a complete analysis [6,7]  with VM distribution amplitudes (DAs) of twist-2 and twist-3 [8]. Such approach to helicity amplitudes of VM electroproduction was used earlier in Ref. [9], where a rather simple model for UGD was used.
Note that the κ-dependence of the IFs is different in the cases of longitudinal and transverse polarizations and this poses a strong constraint on the κ-dependence of the UGD in the HERA energy range. The main point of our work will be to demonstrate, considering different models for UGD, that the uncertainties of the theoretical description do not prevent us from some, at least qualitative conclusions about the κ-shape of the UGD.
In this paper we concentrate on κ-factorization method. The dipole approach is based on similar physical ideas, but formulated not in κ-but in the transverse coordinate space; this scheme is especially suitable to account for nonlinear evolution and gluon saturation effects. Interesting developments are the results of the papers [10,11], where the helicity amplitudes of VM production were considered in the dipole approach.
The paper is organized as follows: in Section II we will present the expressions for the amplitudes of interest here, discuss the sources of theoretical uncertainties and sketch the main properties of a few models for the UGD; in Section III we compare theoretical predictions from the different models of UGD with HERA data; in Section IV we draw our conclusions and give some outlook.

II. THEORETICAL SETUP
The H1 and ZEUS collaborations have provided extended analyses of the helicity structure in the hard exclusive production of the ρ meson in ep collisions through the subprocess Here λ ρ and λ γ represent the meson and photon helicities, respectively, and can take the values 0 (longitudinal polarization) and ±1 (transverse polarizations). The helicity amplitudes T λρλγ extracted at HERA [6] exhibit the following hierarchy, that follows from the dominance of a small-size dipole scattering mechanism, as discussed first in Ref. [12]: The H1 and ZEUS collaborations have analyzed data in different ranges of Q 2 and W . In what follows we will refer only to the H1 ranges, 2.5 GeV 2 < Q 2 < 60 GeV 2 , and will concentrate only on the dominant helicity ratio, T 11 /T 00 .

B. Theoretical uncertainties and approximations
There are four sources of uncertainty and/or approximation in our analysis, as based on the above expressions for the helicity amplitudes.
(i) The γ * → V IF is a function of Q 2 and κ which is not fully perturbative and includes also physics of large distances. Here we use collinear factorization to express the IF as a convolution -integration over longitudinal momentum fraction -of the nonperturbative twist-2 and -3 DAs and a perturbative hard part.
In the region of large κ, κ ∼ Q, which corresponds to the range of small dipole sizes, the IFs for the production of both the longitudinally and transversely polarized meson are well described in our collinear factorization scheme. The neglected contributions are relatively suppressed as powers of Λ QCD /Q and are therefore neglected.
The region of small κ, κ Q, is also present in our κ-factorization formulas and corresponds to the range of larger dipole sizes r d . Can we calculate also in this case our IFs as convolution of the perturbative hard part with the meson DAs?
The situation here is different in the cases of longitudinal and transverse polarizations.
In the longitudinal case, we have, in fact, small r d dominance in the region of all κ. Note that, as calculated in our scheme, the longitudinal IF divided by κ 2 is finite for κ → 0, and, therefore, the neglected contributions for our longitudinal IF are suppressed by powers of Λ QCD /Q in the region of all κ. In relation with that, we note here that the longitudinal VM electroproduction can be described not only in κ-, but also fully in QCD collinear factorization (in terms of generalized parton distributions).
For the transverse polarization the situation is different, since the transverse IF divided by κ 2 behaves like log(κ 2 /Q 2 ), which means that the collinear limit, κ → 0, is not safe and one cannot describe the transverse VM electroproduction fully in QCD collinear factorization.
One can easily trace that this behavior ∼ log(κ 2 /Q 2 ) for κ → 0 of the transverse IF appears due to the integration over the longitudinal fraction z and its logarithmic divergence near the end points z → 0, 1. The light-cone wave function of the virtual photon, that controls the hard part of the IF, has a scale Q 2 r 2 d z(1 − z). This means that, at κ → 0, when the endpoint region of the z-integration is important, large values of Q do not mean automatically that the small-r d region is dominant, but instead both large and small r d contribute. Therefore we do not control the accuracy of our IF calculation for κ → 0 in the transverse polarization case.
However, in κ-factorization the small-κ region is only a corner of the κ integration domain and the importance of this corner is a matter of investigation. On the experimental side we do not see indications that large r d 's, and therefore small κ's, are dominant; indeed, HERA data show similar t-and W -dependence for the both longitudinal and transverse helicity amplitudes. In our phenomenological analysis we will check the importance of the region of small κ's by studying the dependence of our predictions on the κ lower cut value.
(ii) Another source of uncertainty comes from the adopted form of the light-cone DAs.
We considered, for the sake of simplicity, the so called asymptotic choice for the twist-2 DA given in Eq. (6), corresponding to fixing a 2 (µ 2 ) = 0. The impact of this approximation was estimated by letting a 2 (µ 2 0 ) take a non-zero value as large as 0.6 at µ 2 0 = 1 GeV 2 in the analysis with one specific model for UGD.
We used typically twist-3 DAs in the Wandzura-Wilczek approximation, but considered in one case the effect of the inclusion of the genuine twist-3 contribution to check the validity of this approximation.
(iii) We calculate the forward amplitudes for both longitudinal and transverse case. The experimental analysis showed that the t-dependence is similar for the two helicity amplitudes, the measured values of the slope parameter have, within errors, the same values for the both polarizations cases. Therefore in T 11 /T 00 ratio considered here the t-dependence drops.
(iv) The expression in Eq. (4) represents, as a matter of fact, the imaginary part of the amplitude and not the full amplitude. The real parts of the amplitudes at high energy are smaller: they are suppressed in comparison to the imaginary parts by the factor ∼ 1/ log s, and they are related to the latter by dispersion relations. Here again we appeal to the results of the experimental analysis, that showed similar W -dependence for both helicity amplitudes, that means the effective cancellation of the contribution from the real parts of the amplitudes in the ratio T 11 /T 00 .

C. Models of Unintegrated Gluon Distribution
In this work we have considered a selection of six models of UGD, without pretension to exhaustive coverage, but with the aim of comparing (sometimes radically) different approaches. We refer the reader to the original papers for details on the derivation of each model and limit ourselves to presenting here just the functional form F(x, κ 2 ) of the UGD as we implemented it in the numerical analysis.
1. An x-independent model (ABIPSW) The simplest UGD model is x-independent and merely coincides with the proton impact factor [9]: where M corresponds to the non-perturbative hadronic scale. The constant A is unessential since we are going to consider the ratio T 11 /T 00 .

Gluon momentum derivative
This UGD is given by and encompasses the collinear gluon density g(x, µ 2 F ), taken at µ 2 F = κ 2 . It is based on the obvious requirement that, when integrated over κ 2 up to some factorization scale, the UGD must give the collinear gluon density. We have employed the CT14 parametrization [14], using the appropriate cutoff κ min = 0.3 GeV (see Section III A for further details).

Ivanov-Nikolaev' (IN) UGD: a soft-hard model
The UGD proposed in Ref. [15] is developed with the purpose of probing different regions of the transverse momentum. In the large-κ region, DGLAP parametrizations for g(x, κ 2 ) are employed. Moreover, for the extrapolation of the hard gluon densities to small κ 2 , an Ansatz is made [16], which describes the color gauge invariance constraints on the radiation of soft gluons by color singlet targets. The gluon density at small κ 2 is supplemented by a non-perturbative soft component, according to the color-dipole phenomenology.
The soft term reads where C F = N 2 c − 1 2N c and µ soft = 0.1 GeV. The parameter a soft = 2 gives a measure of how important is the soft part compared to the hard one. On the other hand, the hard component reads where F pt (x, κ 2 ) is related to the standard gluon parton distribution as in Eq. (10) and Q 2 c = 3.26 is the soft-hard interface (see Section IIIA for further details). We refer to Ref. [15] for the expressions of the vertex function V N (κ) and of F We stress that this model was successfully tested on the unpolarized electroproduction of VMs at HERA.

Hentschinski-Salas-Sabio Vera' (HSS) model
This model, originally used in the study of DIS structure functions [17], takes the form of a convolution between the BFKL gluon Green's function and a LO proton impact factor. It has been employed in the description of single-bottom quark production at LHC in Ref. [18] and to investigate the photoproduction of J/Ψ and Υ in Ref. [19]. We implemented the formula given in Ref. [18] (up to a κ 2 overall factor needed to match our definition), which reads , with N f the number of active quarks (put equal to four in the following),ᾱ s = α s (µ 2 ) N c π , with µ 2 = Q 0 M h , and χ 0 ( 1 2 + iν) ≡ χ 0 (γ) = 2ψ(1) − ψ(γ) − ψ(1 − γ) is the LO eigenvalue of the BFKL kernel and ψ(γ) is the logarithmic derivative of Euler Gamma function. Here, M h plays the role of the hard scale which can be identified with the photon virtuality, Q 2 . In Eq. (15), χ(γ) (with γ = 1 2 + iν) is the NLO eigenvalue of the BFKL kernel, collinearly improved and BLM optimized; it reads with χ 1 (γ) and χ RG (ᾱ s , γ) given in Section 2 of Ref. [18].
This UGD model is characterized by a peculiar parametrization for the proton impact factor, whose expression is which depends on three parameters Q 0 , δ and C which were fitted to the combined HERA data for the F 2 (x) proton structure function. We adopted here the so called kinematically improved values (see Section III A for further details) and given by Q 0 = 0.28 GeV, δ = 6.5, C = 2.35 .

Golec-Biernat-Wüsthoff ' (GBW) UGD
This UGD parametrization derives from the effective dipole cross sectionσ(x, r) for the scattering of a qq pair off a nucleon [20], through a reverse Fourier transform of the expression with (22) and the following values σ 0 = 23.03 mb, λ p = 0.288, The normalization σ 0 and the parameters x 0 and λ p > 0 of R 2 0 (x) have been determined by a global fit to F 2 (x) in the region x < 0.01.

Watt-Martin-Ryskin' (WMR) model
The UGD introduced in Ref. [21] reads where the term gives the probability of evolving from the scale κ to the scale µ without parton emission.
Here z max ≡ 1 − z min = µ/(µ + κ t ); N f is the number of active quarks. This UGD model depends on an extra-scale µ, which we fixed at Q. The splitting functions P qg (z) and P gg (z) are given by with the plus-prescription defined as . (26)

III. NUMERICAL ANALYSIS
In this Section we present our results for the helicity-amplitude ratio T 11 /T 00 , as obtained with the six UGD models presented above, and compare them with HERA data.
We preliminarily present a plot, Fig. 1, with the κ 2 -dependence of all the considered UGD models, for two different values of x. The plot clearly exhibits the marked difference in the κ 2 -shape of the six UGDs.
In Fig. 2 we compare the Q 2 -dependence of T 11 /T 00 for all six models at W = 100 GeV, together with the experimental result. We used here the asymptotic twist-2 DA (a 2 (µ 2 ) = 0) and the WW approximation for twist-3 contributions. Theoretical results are spread over a large interval, thus supporting our claim that the observable T 11 /T 00 is potentially able to strongly constrain the κ-dependence of the UGD. None of the models is able to reproduce data over the entire Q 2 range; the x-independent ABIPSW model and the GBW model seem to better catch the intermediate-Q 2 behavior of data.
To gauge the impact of the approximation made in the DAs, we calculated the T 11 /T 00 ratio with the GBW model, at W = 35 and 180 GeV, by varying a 2 (µ 0 = 1 GeV) in the range 0. to 0.6 and properly taking into account its evolution. Moreover, for the same UGD model, we relaxed the WW approximation in T 11 and considered also the genuine twist-3 contribution. All that is summarized in Fig. 3, which indicates that the approximations made are quite reliable.
The stability of T 11 /T 00 under the lower cut-off for κ, in the range 0 < κ min < 1 GeV, has been investigated. This is a fundamental test since, if passed, it underpins the main underlying assumption of this work, namely that both the helicity amplitudes considered here are dominated by the large κ region. In Fig. 4  gives only a marginal contribution.

A. Tools and systematics
All numerical calculations we done in Fortran, making use of specific CERNLIB routines [22] to perform numerical integrations and the computation of (poly-)gamma functions.
In order to deal with all the considered UGD models, we found advantageous to create a modular library which allowed us to bind and call all sets via a unique and simple interface, serving at the same time as a working environment for the creation of new, user-customized UGD parametrizations.
The uncertainty coming from the numerical 2-dimensional over κ and y in Eqs. (4), (5) and (7) was directly estimated by the Dadmul integrator [22] and it was constantly kept below 0.5%. In the case of the HSS and WMR models, one should take into account of an extra-source of systematic uncertainties coming from the integration on ν (Eq. (15)) and on z and κ 2 t (Eqs. (24) and (25)), respectively. Even in this case, we managed to keep the numerical error very small. Furthermore, it is worth to note that for all the UGD models involving the use of standard PDF parametrizations, it was needed to put a lower cut-off in κ, in order to respect the kinematical regime where each set has been extracted. We gauged the effect of use different PDF parametrizations by making tests with the most popular sets extracted from global fits, namely MMHT14 [23], CT14 [14] and NNPDF3.0 [24], as provided by the LHAPDF Interface 6.2.1 [25], after imposing a provisional cut-off of κ (test) min = 1 GeV. We checked that the discrepancy among the various cases is small or negligible. Then, we did the final calculations by using the CT14 parametrization, which allowed us to integrate over κ up to κ min = 0.3 GeV. Results with the gluon momentum derivative and the WMR model were obtained by imposing such a cut-off, while we adopted a dynamic strategy as for the IN one: we cut the contribution coming from the IN hard component (Eq. (13)), at κ min , while no cut-off was imposed for the soft component (Eq. (12)). With respect to this model, we made further tests by considering the effect of using different DGLAP inputs (see Table I of Ref. [15]) for the parameters Q c , κ h and µ soft entering Eqs. (12) and (13). No significant discrepancy among them was found, so we gave our results for the IN model by using the so-called CTEQ4L DGLAP input.
Following the definition of IN and WMR models, the PDF set was taken at LO, while the NLO one was employed in the gluon momentum derivative parametrization. Moreover, as for the HSS UGD parametrization, we checked that the discrepancy between the so-called improved setup given in Eq. (18) and the standard one (see, e.g. Ref. [18]) is negligible when considering the helicity-amplitude ratio T 11 /T 00 .

IV. DISCUSSION
In this paper we have proposed an observable that is well measured in the experiments at HERA (and could be studied in possible future electron-proton colliders) -the dominant helicity amplitudes ratio for the electroproduction of vector mesons -as a nontrivial testfield to discriminate the models for the unintegrated gluon distribution in the proton.
The main motivation of our study are the features, observed at HERA, of polarization observables for exclusive vector meson electroproduction. In the cases of both longitudinal and transverse polarizations, the measured cross sections demonstrate similar dependencies on kinematic variables: specific Q 2 scaling, t-and W -dependencies that are distinct from the ones seen in soft diffractive exclusive processes. This indicates that the dominant physical mechanism in both cases is the scattering of a small transverse-size, ∼ 1/Q, dipole on a proton.
On the theoretical side we have a description in κ-factorization, where the nonperturbative physics is encoded in the unintegrated gluon distribution F(x, κ 2 ) and in the vector meson twist-2 and twist-3 DAs (which includes both WW and genuine twist-3 contributions), that parameterize the probability amplitudes for the transition of 2-and 3-parton small-transverse-size colorless states to the vector meson.
In our analysis we have considered six models for F(x, κ 2 ), which exhibit rather different shape of κ-dependence in the region, κ 2 ∼ few GeV 2 , relevant for the kinematic of the ρ-meson electroproduction at HERA, as shown in Fig. 1.
In our numerical study we have found rather weak sensitivity of our predictions for the helicity amplitudes ratio to the physics encoded in the meson DAs (though values of longitudinal and transverse amplitudes separately depend strongly on the model for DAs).
As an example, in Fig. 3 we have presented results for the GBW model of F(x, κ 2 ). Here the dominance of the WW contribution over the genuine twist-3 one is clearly seen. Besides, we have found rather moderate dependence of our observable on the shape of twist-2 DA.
Indeed, varying the value of a 2 in a wide range in comparison to the value a 2 (µ 0 ) = 0.18±0.10 obtained from the QCD sum rules [13], we have found small variation of the amplitude ratio, on the level of the experimental errors.
Another important issue is the small-size color dipole dominance that allows us to use results for the γ * → ρ IF calculated unambiguously in terms of the meson DAs. To clarify this question we have introduced a cut-off in the κ-integration and studied the stability of our predictions on the excluded region of small gluon transverse momenta. In Fig. 4, considering again GBW model as an example, we have shown that the sensitivity of our predictions to the region of small κ is indeed not strong, the variation of our results is less or comparable to the data errors.
In this way we have seen that the dominance of the small-size dipole production mechanism is supported both by the qualitative features of the data and by the theoretical calculations in κ-factorization. This gives evidence to our main statement that, having precise HERA data on the helicity amplitude ratio, one can obtain important information about the κ-shape of the UGD. To demonstrate this in Fig. 2 we have confronted HERA data with the predictions calculated with six different UGD models. We have seen that none of the models is able to reproduce data over the entire Q 2 range and that HERA data on the transverse to longitudinal amplitudes ratio are really precise enough to discriminate predictions of different UGD models.
We recommend that further tests of models for the unintegrated gluon distribution, as well as possible new model proposals, take into due account our suggestion to utilize important information encoded in the HERA data on the helicity structure in the light vector meson electroproduction.

V. ACKNOWLEDGMENTS
We thank M. Hentschinski for fruitful discussions.
FGC acknowledges support from the Italian Foundation "Angelo della Riccia". In the twist-2 DA we have put a 2 (µ 0 = 1 GeV) = 0 and the T 11 amplitude has been calculated in the WW approximation.