Non-perturbative gluon-hadron inputs for all available forms of QCD factorization

Description of hadronic reactions at high energies is conventionally done on basis of QCD factoriza- tion so that factorization convolutions involve non-perturbative inputs mimicking non-perturbative contributions and perturbative evolution of those inputs. We construct the inputs for the gluon- hadron scattering amplitudes in the forward kinematics and, using the Optical theorem, convert them into inputs for gluon distributions in the both polarized and unpolarized hadrons. Firstly, we derive general mathematical criteria which any model for the inputs should obey and then suggest a Resonance Model satisfying those criteria. This model is inspired by a simple observation: after emitting an active parton off the hadron, the remaining ensemble of spectators becomes unstable and therefore it can be described through factors of the resonance type. Exploiting Resonance Model, we obtain non-perturbative inputs for gluon distributions in unpolarized and polarized hadrons for all available forms of QCD factorization: Basic, KT - and Collinear Factorizations. We compare the obtained inputs to the inputs available in the literature.


I. INTRODUCTION
Description of hadronic reactions at high energies requires the use of QCD in the perturbative and non-perturbative domains. It is well-known that such calculations cannot be carried out in the straightforward way because the nonperturbative QCD is poorly-known. The most efficient approximation approach to theory of hadronic reactions is QCD factorization, where perturbative QCD calculations are complemented by phenomenological fits or model expressions which mimic non-perturbative QCD contributions. Throughout the present paper we will address those expressions as non-perturbative inputs. Technically speaking, such inputs act as initial conditions for the evolution equations which account for perturbative contributions. The inputs are supposed to approximately describe short-time dissociations of each of the interacting hadrons into active partons and spectators. Depending on the number of the active partons emitted by each hadron, there can be Single-Parton Scattering (SPS) and Multi-Parton Scattering (MPS). An example of the factorized amplitude for the scattering of hadrons under the SPS approximation is depicted in Fig. 1. When such amplitudes are conjugated with the mirror amplitudes, the intermediate states consist of two partons. For instance, QCD factorization of the hadronic tensor W µν for DIS off a hadron with momentum p is represented in Fig. 2, where the s-cut is implied. The lowest blob stands for a non-perturbative input describing emition/absorption of the active partons with momentum k from the initial hadron while the upper blob corresponds to DIS off the active partons. The Optical theorem relates the elastic scattering amplitudes in the forward kinematics to cross-sections. For example, it relates A µν to W µν : So, in order to calculate cross-sections or parton distributions under the SPS approximation, it is sufficient to calculate the related elastic amplitudes with two-parton intermediate partons in the t-channel only. In order to avoid misunderstanding, we note that t-channel intermediate states inside the upper (perturbative) blob can involve unlimited number of partons. Despite that there is a quite extensive literature on MPS scenario, the SPS scenario is still the most popular. By this reason we focus on it in the present paper.
There are different kinds of QCD factorization in the literature and each of them is tailored to a specific perturbative approach. For example, when DGLAP [1] or its generalization [2] to the small-x region are used to calculate W µν , the first step is to introduce an arbitrary mass scale µ acting as a starting point of the perturbative k 2 ⊥ -evolution from µ 2 to Q 2 , with Q 2 ≫ µ 2 . This scale is called the factorization scale. At the same moment, µ can act as a cut-off for infrared-divergent perturbative contributions. In Collinear Factorization [3], transverse components of momentum k in Fig. 2 are neglected, so these partons are regarded as collinear to the incoming hadrons and therefore one-dimensional non-perturbative inputs can be used.
In contrast, BFKL [4] is free of infrared divergences and because of that integrations over transverse momenta in the BFKL ladders run down to zero. This approach operates with essentially non-collinear initial partons, which excludes neglecting their transverse components and as a result, it excludes a simple matching of BFKL with Collinear Factorization. By this reason, K T -Factorization [5] named also High-Energy Factorization [6] was suggested. These factorizations use different ways of parametrization of momenta k of the connecting partons. In Collinear Factorization, the parametrization of k is one-dimensional: with β being the longitudinal momentum fraction. K T -Factorization involves the same longitudinal parameter β and, in addition, accounts for the transverse momentum k ⊥ : There are different ways to construct the non-perturbative inputs in Collinear and K T -Factorizations. Quite often, see e.g. Ref. [7], the inputs applied in the context in Collinear Factorization are introduced entirely on basis of practical reasons without any theoretical grounds. In contrast, there also are the models with solid theoretical background. In particular, the models in Ref. [8] are actually based on various theoretical approaches to approximate the problem of confinement: the chiral quark-soliton model, diquark model, etc. An overview of the most popular models of hadrons can be found in Ref. [9]. The features of the saturation model [10] are used in Refs. [11,12] for modeling the inputs in the context of K T -Factorization while the model in Ref. [13] combines features of a variety of models for the Fock space wave function on the light cone. Some features of this model are similar to the results obtained in Ref. [14]. Another interesting approach is the lattice calculations. They are a complementary method to study the non-perturbative inputs for parton distributions. Some of the recent results of the lattice calculations can be found is Ref. [15].
Although the parametrization of momentum k in K T -Factorization is more general than in Collinear Factorization, it misses one more parameter for the longitudinal component of k, which can be seen from comparison of Eq. (3) to the standard Sudakov parametrization [16]: where the light-cone (i.e. p ′2 = q ′2 = 0) momenta p ′ , q ′ are made of the external momenta p and q satisfying the inequality |pq| ≫ |p 2 |, |q 2 | (for instance, they can be the momenta p, q introduced in Fig. 2): In Ref. [17] we presented a new kind of QCD factorization: Basic Factorization. It accounts for dependence of the factorized blobs on all components of momentum k and because of that it is the most general form of factorization. We proved that Basic Factorization can step-by-step be reduced to K T -Factorization and to Collinear Factorization. Doing so, we considered the non-perturbative inputs in the most general form, without specifying them. In Ref. [17] we inferred general theoretical constraints on the non-perturbative inputs. They follow from the obvious requirement: despite the integrands in the factorization convolutions for forward scattering amplitudes have both infra-red (IR) and ultra-violet (UV) singularities, integrations over k in the factorization convolutions must yield finite results. This is possible only if the non-perturbative inputs kill all divergences, which leads to restrictions on the non-perturbative inputs. In Ref. [18] we suggested a model for the non-perturbative inputs to quark distributions in hadrons. This model is based on a simple observation: after emitting an active parton off the hadron, the set of remaining quarks and gluons (usually named spectators) becomes unstable, so it can be described by expressions of the resonance type. Because of that we named this model in Ref, [18] the Resonance Model. We used the Resonance Model model so as to obtain the quark distributions in the hadrons whereas the very important case of gluon distribution was left uninspected in Ref. [18].
There is a certain similarity between handling quarks and gluons in Perturbative QCD. It allows us to believe that main features of non-perturbative inputs for gluons and quarks can be alike. However, the quark and gluon inputs cannot coincide. Firstly, they have different polarization structure: the quark input should be a spinor whereas the gluon one should be a tensor. Secondly, the well-known difference between the high-energy behavior of the gluon and quark perturbative amplitudes can lead to an essential difference between the quark and gluon inputs.
We think that this issue needs a thorough investigation and do it in the present paper. In the first place we obtain restrictions on the gluon inputs which guarantee both IR and UV stability of the factorization convolutions in Basic Factorization and then use those restrictions so as to extend the Resonance Model to description of the non-perturbative inputs for the gluon distributions in the hadrons. We consider here both polarized and unpolarized hadrons. Our strategy is to calculate the gluon-hadron elastic scattering amplitudes in the forward kinematics and, using the Optical theorem, to arrive at the gluon distributions in the hadrons in Basic Factorization. Then we reduce them down to the expressions which can be used in K T -and Collinear Factorizations. As is well-known, K T -Factorization (High-Energy Factorization) by definition can be used in the Regge (small-x) kinematics, so throughout the present paper we will consider the parton distributions in this region only.
Our paper is organized as follows: In Section II we study the elastic gluon-hadron amplitudes for the forward kinematic region in the Born approximation, and then analyze the impact of radiative corrections. We investigate the convergence of the factorization convolutions, using a general form for the non-perturbative inputs, to determine constraints on the non-perturbative inputs. In Section III we formulate general criteria for the non-perturbative inputs in Basic Factorization and use them to construct a Resonance Model. Transition from the Basic Factorization to K T -Factorization is considered in Sect. IV. In this Sect. we also compare our results with the ones available in the literature. In Sect. V we reduce the inputs in K T -Factorization down to the ones in Collinear Factorization. We focus here on comparing our results with the standard DGLAP-fits. At last, Sect. VI is for our concluding remarks.

II. ELASTIC GLUON-HADRON SCATTERING AMPLITUDES IN THE FORWARD KINEMATICS
In this Sect. we consider the elastic gluon-hadron amplitude A in the forward kinematics and inspect integrability of convolutions for such amplitude in Basic Factorization. We study the convolutions involving two-gluon intermediate t-channel states only because we presume the approximation of Single-Parton Scattering for the gluon distributions. The convolutions describing the gluon-hadron amplitude are depicted by the sum of two graphs 1 : the graph in Fig. 3 and a similar graph, where q is replaced by −q. Our aim here is to obtain mathematical restrictions on the non-perturbative inputs. In order to do it, we use the obvious reasoning: Integration over momentum k in the factorization convolutions (see Fig. 3) runs over the whole phase space and it should yield a finite result. Features of perturbative components of the convolutions are known, at least qualitatively, so integrability requirement can be used to establish necessary restrictions on the non-perturbative inputs in a general form. Then any model for the inputs should meet these restrictions. Therefore, the restrictions can be used as criteria for accepting or rejecting the models. To begin with, we consider amplitude A in the Born approximation and then examine an impact of the radiative corrections. We show that, in contrast to the concepts of K T -and Collinear Factorizations, Basic Factorization allows the simple scenario, where the upper blob in in Fig. 3 can be regarded as altogether perturbative one while the lowest blob includes non-perturbative contributions only.

A. Gluon-hadron scattering amplitudes in the Born approximation
Let us consider the simplest case of the factorization convolutions for the gluon-hadron amplitude, where the nonperturbative input is convoluted with the gluon-gluon perturbative amplitude in the Born approximation. In this case the upper blob in Fig. 3 is represented by two simple graphs, each with the single-gluon exchange. The first graph, with a non-zero imaginary part in s, is depicted in Fig. 4 and the second graph, with a non-zero imaginary part in u, can be obtained from the first one by replacement q → −q. A remarkable feature of Basic Factorization is that the analytic expression corresponding to the convolution graphs can be obtained with the use of the standard Fetnman rules. Applying these rules to the graph in Fig. 4 and complementing the result by the exchange q → −q, we arrive at the following expression for the Born gluon-hadron amplitude A B in Basic Factorization: where N = 3 is the color factor and S denotes the hadron spin. The term in squared brackets corresponds to the Born amplitude for gluon-gluon scattering, l µ (q) and l ′ ν (q) denote the polarization vectors of the external gluons, the terms 1/k 2 correspond to propagators of the connecting gluons and the input T λρ contains non-pertrurbative contributions only. Throughout the paper we use the Feynman gauge. Let us notice that both T λρ and the term in the squared brackets are dimensionless. The perturbative term H µνλρ in Eq. (6) is We remind that in the present paper we consider the gluon-hadron scattering amplitudes and parton distributions in the small-x region. The point is that we are going to reduce Basic Factorization to K T -Factorization which is defined in the small-x region only. It is easy to check that A B with H µνλρ given by Eq. (8) is gauge-invariant in the small-x, though approximately 2 In order to make further progress we should specify T λρ . First of all, we should fix its polarization (tensor) structure. For instance, when the hadron and the gluons are unpolarized, T λρ should be symmetric in λ, ρ. The general form of such a tensor, T but it involves four arbitrary invariant amplitudes A, .., D. Trying to diminish their number, we recall an instance from perturbative QCD, namely that, when the hadron is replaced by a bare quark, T λρ is replaced by the quark-gluon Born amplitude Q λρ . It consists of the unpolarized part Q U λρ and the spin-dependent part Q S λρ : with where m q is the quark mass, S q is the quark spin and C F = 4/3. We suggest the simplest generalization of Eqs. (10,11), for the hadrons: we presume that T λρ keeps the polarization structure of Q λρ . It means that with where M h , S are the hadron mass and spin respectively, and The invariant amplitudes M U and M S may either coincide or differ from each other. Substituting T U,S λρ in Eq. (6), we obtain: with We have done summation over the gluon polarizations in the expression for N (1) s in Eq. (17). In what follows we will use the Sudakov variables defined in Eq. (4). In their terms Let us consider integrability of the Born amplitudes in Eqs. (15,16). Integration over momentum k in both convolutions runs over the whole phase space and result of the integration must be finite despite the integrands may have singularities which must be regulated. First of all, let us note that the expressions in parentheses in Eqs. (15,16) stand for the gluon perturbative amplitudes in the Born approximation, which are obviously free of divergences. In contrast, the factors 1/k 2 in Eqs. (15,16) become divergent at k 2 = 0. We cannot implement any IR cut-off to regulate those singularities because there is no physical reason to restrict the integrations regions. It leaves us with the only option: amplitudes M U,S in Eqs. (15,16) at small k 2 must decrease rapidly enough in order to kill these IR singularities. Therefore we infer that M U,S should satisfy the following restrictions at k 2 → 0: with η > 0. Eq. (19) sets IR stability for the integrals in Eqs. (15,16). Then, there may be UV singularities in Eqs. (15,16) at large |k|, i.e. at large |α| in terms of the Sudakov variables, when the integration over α goes first. It is easy to see that the integrands in Eqs. (15,16) at large |α| behave as which can be UV-divergent. In order to kill this divergence, M U,S should decrease at large α: with χ > 0. In other words, the non-perturbative inputs should decrease, when their invariant energy s 1 grows. Also, we note that the restrictions in Eqs. (19,21) are valid at any choice of the polarization structure of T λρ , including the most general parametrization in Eq. (9).

C. Gluon-hadron scattering amplitudes beyond the Born approximation
Diagrammatically, accounting for the radiative corrections to the Born amplitudes A B U,S of Eqs. (15,16) means adding more quark and gluon propagators to the single-gluon exchange in Fig. 4 and summing up contributions of such graphs. This procedure converts the Born perturbative amplitudes of Eqs. (15,16) into amplitudes A (pert) U,S (the upper blob in Fig. 3) but keeps unchanged the non-perturbative inputs M U,S (the lowest blob in Fig. 3). Such natural separation of the perturbative and non-perturbative contributions is a remarkable feature of our approach and perfectly agrees with the factorization concept.
Replacing the Born amplitudes in Eqs. (15,16) by the amplitudes A U,S , we arrive at the gluon-hadron amplitudes A U , A S beyond the Born approximation: Applying the Optical theorem to Eq. (22), we obtain the spin-independent, D U , and polarized, D S , distributions of gluons in the hadrons: where We remind that the invariant sub-energies s 1,2 are defined as follows: depends on s 2 , k 2 , q 2 , so, being dimensionless, it can be parameterized, for instance, as follows: Perturbative higher-loop contributions can induce both IR and UV divergences in A (pert) U,S . All those divergences should be regulated before using A (pert) U,S in the factorization convolutions (22). We start with considering regularization of IR divergences in A (pert) U,S . They appear from integrations of gluon loops over the regions, where some of the virtual gluons are soft. In addition, the problem of IR divergences is essential for soft quarks when their masses are neglected. The ways of regulating IR divergences are different for the amplitudes with on-shell and off-shell external partons. We will address such amplitudes as the on-shell and off-shell amplitudes respectively. For the on-shell amplitudes, the IR divergences are regulated by suitable IR cut-offs, while virtualities of the external partons play the role of natural IR cut-offs for the off-shell amplitudes. The latter was shown first in Ref. [16] and then was used in numerous calculations. Therefore, q 2 and k 2 can be used as IR cut-offs for A are ∼ ln n (wβ/k 2 ) (or ∼ ln n (q 2 /k 2 )), so these contributions are IR stable at k 2 = 0 but become divergent at the point k 2 = 0. Integration in Eqs. (22) covers the whole phase space, so the integration region includes the point k 2 = 0. Such logarithmic divergences were addressed as mass singularities in the pioneer studies [3] of QCD hard processes in the framework of Collinear Factorization, where the mass singularities were moved from the hard part to universal structure functions. However, such a treatment of the mass singularities becomes unproductive in the Regge (small-x) kinematics which we pursuit throughout this paper. Furthermore, the logarithmic singularities are the only kind of IR divergences in Collinear Factorization but in Basic Factorization we should handle both them and the power singularity ∼ 1/(k 2 k 2 ) (see Eqs. (22)). In Sect. IIB we demonstrated that the power singularity 1/(k 2 k 2 ) can be regulated by the condition (19). Obviously, this condition persists when the singularity 1/(k 2 k 2 ) is accompanied by the logarithmically divergent terms. It is worth noticing here that the well-known power factor s/k 2 in the spin-independent amplitude appears in the lowest order in α s , borrowing k 2 from the denominator in Eqs. (22), so its impact on IR stability was taken into account in Eq. (19).
Finally, it is easy to prove that accounting for the perturbative radiative corrections does not change the condition (21) of UV stability of the factorization convolutions at large |α|. Indeed, the fact that the perturbative QCD is renormalizable means that the fastest-growing with energy terms in amplitudes A (pert) U,S are of logarithmic kind, i.e. ∼ ln n wα/q 2 , ln n wα/k 2 , so such logarithmic divergences can be regulated by the power restriction in Eq. (21).

III. MODELING THE NON-PERTURBATIVE INVARIANT AMPLITUDES MU,S
In this Sect. we consider first the Resonance Model for non-perturbative amplitudes M U,S and then proceed to the non-perturbative inputs for the gluon distributions. Our models is based on simple and clear-cut physical arguments which we list below and call them criteria which any model for M U,S should satisfy. Although constructing such a model cannot be done unambiguously, any proposed model should satisfy the following criteria: A. General criteria for modeling non-perturbative inputs Criterion (i) Expressions for M U,S should satisfy the requirements of IR and UV stability of the factorization convolutions given in Eqs. (19,21).

Criterion (ii)
Expressions for M U,S should have non-zero imaginary parts in s 1 to make possible the use of the Optical theorem. It is necessary in order to proceed from elastic gluon-hadron amplitudes to gluon distributions in the hadrons. Criterion (iii) Expressions for M U,S should enable the step-by-step reduction of Basic Factorization down to K Tand Collinear Factorizations. Such a reduction was suggested in general in Ref. [17].

B. Specifying non-perturbative gluon inputs in Basic Factorization
We suggest the following model expressions for the amplitudes M U,S (s 1 , k 2 ) : Eq. (26) reads that dependences on s 1 and k 2 are separated. This is done for the sake of simplicity. We could assume that R = R(k 2 , s 1 ) but we find it unnecessary. As handling M U and M S is much the same, we skip the subscripts U, S and will work with the generic notations This is the most UV divergent case which should not be ignored, so in order to satisfy Eqs. (21,29), Z should behave as follows: at large |α| If R does not grow at large |α|, we can all the same use Z obeying Eq. (30). In order to specify Z, we suggest the guiding idea: After emitting an active parton by the initial hadron, the remaining set of quark and gluons pick up a color and therefore it cannot be stable. Unstable states, in general, are known to be often modeled by expressions of the resonance type. In particular, we can approximate Z in Eq. (27) by the following expression: with arbitrary positive integer n. In order to get Z n (s 1 ) satisfying Eq. (30), we should choose n > 1. For the sake of simplicity, we consider the minimal value n = 2. Obviously, Z 2 can also be written as an interference of two resonances: In terms of the Sudakov variables, where C Z = ∆M 2 12 + ı∆Γ 12 and µ 2 1,2 = M 2 1,2 − p 2 . Applying the Optical theorem to Eq. (33) allows us to obtain the non-perturbative contribution G to the gluon distributions in the hadrons: Obviously, this expression is of the Breit-Wigner type.

IV. NON-PERTURBATIVE GLUON INPUTS IN KT -FACTORIZATION
The expression for the gluon-hadron scattering amplitude M (α, β, k 2 ) in Eq. (33) corresponds to Basic Factorization. In order to reduce M (α, β, k 2 ) to the amplitude M KT (β, k 2 ⊥ ) of the same process in K T -Factorization, we should integrate out the α-dependence of M (α, β, k 2 ). In principle, as soon as However, in factorization convolutions (see e.g. Figs. 2,3,4 ) both the upper and the lowest blobs depend on α, so one cannot integrate M (α, β, k 2 ) only, without integration of the perturbative blob. Therefore, M (β, k 2 ⊥ ) cannot be derived from M (α, β, k 2 ) in the straightforward way, which makes impossible a straightforward reduction of Basic Factorization to K T -Factorization. Nevertheless, it can be done approximately. The perturbative blobs (the upper blob in Figs. 2,3,4) approximately do not depend on α in the region and the non-perturbative blobs in this region are independent of β. Indeed, the upper blob depends on α through k 2 only and in the region (35) k 2 ≈ −k 2 ⊥ . At the same time, the region (35) is known to be the region for gaining the perturbative contributions most essential at small x. For instance, all BFKL logarithms come from this region. So, Eq. (35) makes possible to integrate M α, β, k 2 independently of perturbative contributions. It is convenient to define the non-perturbative input M KT in K T Factorization as follows: with So, Eq. (37) is satisfied when with a positive ξ obeying the inequality ξ ≪ 1. Combining Eqs. (33,36,38) we arrive at the following expression for the non-perturbative input M KT in K T -Factorization: where R KT = R/(ξk 2 ⊥ C Z ) and We also used the well-known observation that the most essential region in the factorization convolutions is the small-β region, where ξ/β is not far from x. Eq. (39) is valid when i.e. when k 2 ⊥ /β is pretty far from the resonance region k 2 ⊥ /β ∼ µ ′2 1,2 . The closer k 2 ⊥ /β is to µ ′2 1,2 the greater are corrections to Eq. (39). It means that transition from Basic Factorization to K T -Factorization considered in Ref. [17] is straightforwardly done outside the resonance region. After that, Eq. (39) can be analytically continued in the resonance region. As a result, we arrive at the gluon-hadron amplitude A KT and the gluon distribution D KT in K T -Factorization: where A  24)) are the non-perturbative inputs. It is convenient to represent G KT in the following form: with and G R and G B obviously consist of expressions of the Breit-Wigner type. Signs of µ ′ 2 1 and µ ′ 2 2 cannot be fixed a priory, but G R ⇆ G B when µ ′ 2 1,2 → −µ ′ 2 1,2 , so we can consider the case of positive µ ′ 2 1,2 without loss of generality. We remind that k 2 ⊥ /β > 0, so G R is within the resonant region k 2 ⊥ /β ∼ µ ′2 1,2 while G B is out of that region. Eqs. (45,46) allows us to represent the gluon distribution D KT in K T -Factorization as the sum of its resonance part D R KT and the background contribution D B KT : with and Representation of D KT through the resonance D R KT and background D B KT contributions is similar to the structure of expressions in the Duality concept.
A. Discussion of the distribution GKT First of all, let us compare the mass parameters µ 1,2 of the gluon distribution D of Eq. (34) in Basic Factorization and the parameters µ ′ 1,2 of the gluon distribution D KT of Eq. (43) in K T -Factorization. They are related by Eq. (40). As ξ ≪ 1, Eq. (40) means that µ ′ 1 ≫ µ 1 , µ ′ 2 ≫ µ 2 . It means that despite the non-perturbative input G is defined in Basic Factorization within the non-perturbative domain, the resonance maximums of G KT are located at perturbative values of k 2 ⊥ /β > 0. Then, we would like to comment on the factor R KT in Eqs. (39,45,46). This factor describes a dependence of G KT on k 2 ⊥ . Refs. [10][11][12][13] suggest the exponential form of R KT : with N and λ in Eq. (50) being parameters. However, Eq. (50) contradicts to the general requirement of IR stability (cf. Eq. (28)) stating that at k 2 ⊥ → 0, so R ′ in Eq. (50) must be modified. The simplest modification compatible with Eq. (51) is which agrees with the expressions used in Refs. [11,12]. So,Eqs,(45,46,52) present the complete set of formulas to describe the non-perturbative input G KT for the gluon distributions in K T Factorization. It is interesting to note a further similarity of our expressions in Eqs. (45,46) and the model suggested in Ref. [13], where non-perturbative distributions for the five-quark state are studied. The quark distributions in Ref. [13] involve both a Gaussian (though without the power factor k 2 ⊥ η ) and propagators ∼ (k 2 ⊥ /β − m 2 + ıε) −1 of the five quarks. The latter resembles our factor Z in Eq. (32) save the difference between the factors ıε and ıΓ 1,2 . The quarks in Ref. [13] are assumed to be free and stable, the gluon content is dropped. In contrast to Ref. [13], we do not specify the content of the spectators but describe the whole set of spectators through the resonances and background contributions.

V. GLUON DISTRIBUTIONS IN COLLINEAR FACTORIZATION
Our next aim is to reduce D KT of Eq. (43) to the gluon distribution D col in Collinear Factorization and bring the gluon distribution to the following form: with D (pert) col being the perturbative contribution and ϕ being the non-perturbative input; µ F is the factorization scale. To this end, the k ⊥ -dependence in Eq. (43) should be integrated out. At the first sight, such a reduction is impossible: D (pert) KT and G KT in Eq. (43) explicitly depend on k ⊥ , so they both must be integrated and because of that the integration can yield some entangled mixture of the perturbative and non-perturbative terms instead of their product as represented in Eq. (53). However, the specific form of the inputs in Eq. (45,46) allows us to reduce D KT to D col though approximately. According to Eq. (47), D KT consists of D R KT and D B KT . Let us consider first the resonance component, D R KT given by Eq. (48). Integration over ζ therein runs in the resonance region, so only G R must be integrated: where ζ 0 = w(1 − x/β) ≈ w, with w = 2pq as defined in Eq. (5). We presume that the resonances are narrow: In Eq. (54) the notations D (pert) col (x/β, µ ′2 r ), with r = 1, 2, stand for the perturbative components while non-perturbative inputs ϕ R (β, µ ′2 r ) are In Eq. (56) we have denoted N r = N µ 2 r η−1 and C r = λµ 2 r . Besides, we presumed that λµ 2 r < 1 to truncate the series of the power expansion of the exponential. In contrast to D R KT , integration over k ⊥ (over ζ in fact) in Eq. (49) runs outside the resonance region and therefore it necessitates integration of both G B and D (perp) KT . Because of that the outcome scarcely can be represented in the factorized form of Eq. (53), which violates the reduction to Collinear Factorization. Fortunately, this violation is small. Indeed, so when the resonance input G R KT is narrow, i.e. when Γ ′ 1,2 ≪ µ ′2 1,2 , we can neglect all in Eqs. (54,57), thereby arriving at the following expression for the gluon distribution D col in Collinear Factorization: The integrand in Eq. (58) looks similar to the conventional one in Eq. (53) but does not coincide with it. It is easy to bring Eq. (58) to the conventional form (53) with the perturbative evolution in the k ⊥ -space. It can be done in two steps: First, let us introduce a scale µ F so that |q 2 | > µ F > µ ′ 1,2 . Then, using DGLAP (or another evolution approach) we evolve ϕ R (β, µ ′2 1 ) and ϕ R (β, µ ′2 2 ) in Eq. (58) from their scales up to the scale µ F , keeping β fixed. This procedure automatically sets D (pert) col (x/β, µ ′2 1,2 ) in Eq. (58) on the scale µ F . As a result, we convert Eq. (58) into Eq. (53), where ϕ col is related to the non-perturbative input ϕ R by the operator E of the DGLAP-evolution: The evolution operator E in the Mellin (momentum) space is expressed in terms of the DGLAP anomalous dimensions. In contrast to the non-perturbative inputs ϕ R (β, µ ′2 r ), the input ϕ col comprises both non-perturbative (through operators E) and perturbative (through ϕ R ) contributions.

A. Comparison of Eq. (56) to the standard DGLAP fits
Let us compare the non-perturbative input ϕ R with the standard DGLAP fit ϕ. Quite often (see e.g. Ref. [7]) ϕ is chosen in the following form: where N, a, b, c, d are phenomenological parameters defined from analysis of experimental data. All these parameters implicitly depend on the factorization scale µ F . Its value can be chosen arbitrary but factorization convolutions, where ϕ participates, do not depend on µ F . First of all, let us note that the singular factor β −a in Eq. (60) provides the fast growth of the parton distributions at small x and leads to their Regge asymptotics. We proved (see e.g. the overview [2]) that the singular factors in the DGLAP fits mimic accounting for the total resummation of ln n (1/x). When such resummation is taken into account, the factors β −a become irrelevant and should be dropped. On the other hand, our analysis does not exclude the factor β η in Eq. (56) because η is positive.
The factor (1 − β) b of Eq. (60) looks as the large-x asymptotics of of the parton distributions at x → 1. By this reason, we suggest that the origin of this factor also is perturbative, so the factor (1 − β) b can be excluded from the fit when the total resummation of ln n (1 − x) is taken into account.
After the factors β −a and (1 − β) b have been excluded from Eq. (60), it becomes quite similar to the input ϕ R .

VI. SUMMARY
In the present paper we have performed the detailed study of the general structure of the gluon non-perturbative inputs for the gluon distributions in hadrons and applied the Resonance Model to construct the gluon inputs. These inputs can also apply to description of various hadronic high-energy processes, including the DIS structure functions. We constructed the gluon inputs in Basic Factorization and then calculated the gluon inputs in the reduced forms so as to use them in K T -and Collinear Factorizations.