Combining the small-x evolution and DGLAP for description of inclusive photon induced processes

TL;DR

This paper develops a unified theoretical framework combining DGLAP and small-x resummation techniques to accurately describe inclusive photon-induced processes in DIS across all x and Q^2 ranges, improving the understanding of high-energy scattering.

Contribution

It introduces a novel approach that merges DGLAP evolution with double logarithm resummation, providing expressions valid at arbitrary x and Q^2, and offers a new perspective on Pomeron contributions in high-energy QCD.

Findings

Derived expressions for F_1 valid at all x and Q^2.

Unified description matching DGLAP and small-x limits.

Parton-parton amplitudes can replace BFKL Pomeron in models.

Abstract

Interest in studying the inclusive photon induced processes of Deep-Inelastic Scattering (DIS) and Diffractive DIS (DDIS) at high energies includes their experimental investigation and thorough theoretical description. The conventional instrument for theoretical investigation of DIS at large x is DGLAP. It describes the Q^2-evolution. Description of DIS in the small-x region requires accounting for the both Q^2 and x -evolution at the same time. Combining DGLAP with total resummation of double logarithms of x and Q^2, we present a description of structure function F_1 at large Q^2 and arbitrary x. Making use of the necessity of the shift of Q^2 in order to regulate infrared singularities, we obtain an expression for F_1 valid at arbitrary x and Q^2. The obtained expressions coincide with the DLA expressions at small x and at the same time coincides with the DGLAP result at large x and large Q^2. Accounting for virtualities k^2 of the external partons allows us to obtain expressions for F_1 in K_T-Factorization, which are valid at arbitrary Q^2 and arbitrary relations between Q^2 and k^2. Expressions for F_1 in all types of QCD factorization exhibit the small-x asymptotics of the Pomeron type. This Pomeron is generated by the high-energy asymptotics of amplitudes of the 2 -> 2 parton-parton forward scattering. These amplitudes are calculated in DLA, so their asymptotics have nothing to do with the BFKL Pomeron. We demonstrate that the parton-parton amplitudes can be used in the DDIS models instead of BFKL Pomeron or combinations of hard and soft Pomerons at non-asymptotic energies. On the other hand, this modification automatically ensures the Pomerons asymptotics, which simplifies and makes more consistent theory of DDIS.

Equations

$$T=T_q\otimes {\Phi }_q+T_g\otimes {\Phi }_g,$$

$$A_{\gamma }\equiv T{|}_{q^2=0}.$$

$$A_{\gamma }=A_{\gamma q}\otimes {\Phi }_q+A_{\gamma g}\otimes {\Phi }_g.$$

$$T^{(col)}={\int }_x^1\frac{d\beta }{\beta }\left[T_q^{(col)}(x/\beta ,Q^2/{\mu }^2){\Phi }_q^{(col)}(\beta ,{\mu }^2)+T_g^{(col)}(x/\beta ,Q^2/{\mu }^2){\Phi }_g^{(col)}(\beta ,{\mu }^2)\right],$$

$$T_q^{(col)}(x/\beta ,Q^2/{\mu }^2)$$

$$T_g^{(col)}(x/\beta ,Q^2/{\mu }^2)$$

$$y=\ln (Q^2/{\mu }^2)$$

$${\xi }^{(+)}(\omega )=-\left(e^{-ı\pi \omega }+1\right)/2.$$

$$F_q^{(col)}=C_q^{(+)}{e}^{{\Omega }_{(+)}y}+C_q^{(-)}{e}^{{\Omega }_{(-)}y},$$

$$F_g^{(col)}=C_g^{(+)}{e}^{{\Omega }_{(+)}y}+C_g^{(-)}{e}^{{\Omega }_{(-)}y}.$$

$$A_{r{r}^{\prime }}(s/{\mu }^2)={\int }_{-ı\infty }^{ı\infty }\frac{d\omega }{2ı\pi }{\left(s/{\mu }^2\right)}^{\omega }{\xi }^{(+)}(\omega )f_{r{r}^{\prime }}(\omega ).$$

$$h_{r{r}^{\prime }}(\omega )=f_{r{r}^{\prime }}(\omega )/(8{\pi }^2).$$

$$F_q^{DGLAP}$$

$$F_g^{DGLAP}$$

$${\tilde{\Omega }}_{(\pm )}={\Omega }_{(\pm )}+{\tilde{\gamma }}_{(\pm )}^{DGLAP}$$

$${\tilde{C}}_{q,g}^{(\pm )}=C_{q,g}^{(\pm )}+\Delta {\hat{C}}_{q,g}^{(\pm )}.$$

$${\tilde{T}}_q^{(col)}(x/\beta ,Q^2/{\mu }^2)$$

$$T_g^{(col)}(x/\beta ,Q^2/{\mu }^2)$$

$$Q^2\to Q^2+m^2,$$

$$m=\mu .$$

$$\mu =2.3{\Lambda }_{QCD}\approx 1GeV.$$

$${\tilde{T}}_q^{(col)}(x/\beta ,Q^2/{\mu }^2)$$

$${\tilde{T}}_g^{(col)}(x/\beta ,Q^2/{\mu }^2)$$

$$\tilde{x}=x+{\mu }^2/w\equiv x+x_0,\tilde{y}=\ln \left((Q^2+{\mu }^2)/{\mu }^2\right).$$

$$T_q^{(col)}$$

$${\omega }_H^{DL}=1.35,{\omega }_S^{DL}=\mathrm{1.07.}$$

$$x_{max}=10^{-6}.$$

$$T_q^{(col)}\sim \frac{{\Pi }_q({\omega }_0)}{{\ln }^{3/2}(1/x)}{\left(\zeta /{\mu }^2\right)}^{{\omega }_0/2},T_g^{(col)}\sim \frac{{\Pi }_g({\omega }_0)}{{\ln }^{3/2}(1/x)}{\left(\zeta /{\mu }^2\right)}^{{\omega }_0/2},$$

$$\zeta =Q^2/x^2.$$

$$T_{q,g}^{(DGLAP)}\sim x^{-a}(Q^2/{\mu }^2{)}^{{\gamma }_{DGLAP}/b},$$

$$\tilde{y}=y+\sum _{n=1}^{\infty }{c}_n{\left(\frac{{\mu }^2}{Q^2}\right)}^n.$$

$${\tilde{T}}_q^{(col)}(x/\beta ,{\tilde{Q}}^2)={\tilde{T}}_q^{(col)}(x/\beta ,Q^2)+\left(\frac{{\mu }^2}{Q^2}\right)C_1+{\left(\frac{{\mu }^2}{Q^2}\right)}^2{C}_2+...$$

$$C_1=\frac{\partial {\tilde{T}}_q^{(col)}(x/\beta ,y)}{\partial y},C_2=\frac{1}{2}\frac{{\partial }^2{\tilde{T}}_q^{(col)}(x/\beta ,y)}{\partial y^2}.$$

Full text

I Introduction

Photon induced inclusive processes such as Deep-Inelastic Scattering (DIS), Photoproduction and Diffractive Deep-Inelastic Scattering (DDIS) are the objects of intensive experimental and theoretical investigation. Theoretical studying these processes is based on the concept of QCD factorization which makes possible to separate perturbative and non-perturbative contributions. There are several forms/types of QCD factorization in the literature. Firstly, there is Collinear Factorizationcolfact . Secondly, there is the more general $K_{T}$ Factorization suggested in Refs. ktfact ; hefact . The third, the most general, type of QCD factorization is Basic Factorization suggested in Ref. bfact . In any of those factorizations the approximation of the single-parton photon-hadron interactions is used and as a result, the subjects of consideration are represented through convolutions of perturbative and non-perturbative objects connected by two-parton intermediate states. For instance, representation of the amplitude $T$ of the elastic forward Compton scattering off a hadron in any type of QCD factorization is depicted in Fig. 1 .

In the analytical form, QCD factorization of $T$ is generically written as follows:

[TABLE]

where $T_{q}$ and $T_{g}$ are the perturbative amplitudes of the elastic Compton scattering off a quark and gluon respectively. Non-perturbative blobs $\Phi_{q}$ and $\Phi_{g}$ denote distributions of those partons in the hadrons. The notations $\otimes$ refer to integrations over momentum $k$ of any parton in the two-quark and two-gluon intermediate states. The specific form of the integrations as well as a parametrization of momentum $k$ of the intermediate partons depend on the type of chosen factorization111See for detail Ref. bfact and refs therein.. Amplitudes $T_{q,g}$ as well as amplitudes $\Phi_{q,g}$ are also given by expressions depending on the type of factorization. In the present paper we focus on the perturbative amplitudes $T_{q,g}$. In any specific type of factorization, there are kinematic regions where amplitudes $T_{q,g}$ are given by essentially different expressions.

Such regions are :

Region A. Large $x$ and large $Q^{2}$: $(x\sim 1)\bigotimes(Q^{2}\gg\mu^{2})$.

Region B. Small $x$ and large $Q^{2}$: $(x\ll 1)\bigotimes(Q^{2}\gg\mu^{2})$.

Region C. Small $x$ and small $Q^{2}$: $(x\ll 1)\bigotimes(Q^{2}<\mu^{2})$.

Region D. Large $x$ and small $Q^{2}$: $(x\sim 1)\bigotimes(Q^{2}<\mu^{2})$.

We have used above the conventional notations: $Q^{2}=-q^{2}$, with $q$ being the external photon momentum, and $x=Q^{2}/w$, where $w=2pq$, with $p$ being the momentum of the incoming/outgoing parton. The parameter $\mu$ is associated with the factorization scale. Besides, it plays the role of the infrared (IR) cut-off in Collinear, $K_{T}$ and Basic Factorizations, when the double-logarithmic (DL) contributions to amplitudes $T_{q,g}$ are accounted for. The physical meaning of $\mu$ for the parton distributions $\Phi_{q,g}$ in $K_{T}$ Factorization and its role in transition from $K_{T}$ Factorization to Collinear Factorization were considered in Ref. bfact . Now let us consider the status of knowledge of amplitudes $T_{q,g}$ in the Regions A-D.

The region A is the DGLAP applicability region, so expressions for $T_{q,g}$ in this region are provided by the DGLAP evolution equationsdglap . Such expressions can easily be found in the literature. They account for logarithms of $Q^{2}$ but neglect total resummations of logarithms of $x$ because such logarithms are unessential in Region ** A**.

Description of $T_{q,g}$ in the small-$x$ region B in the double-logarithmic approximation (DLA) in the framework of Collinear Factorization was obtained in Ref. etf1 . In particular, it was shown in Ref. etf1 that the small-$x$ asymptotics of $T_{q,g}$ exhibits a new, DL contribution to Pomeron. In contrast, double-logarithmic expressions for $T_{q,g}$ in region **B ** in $K_{T}$ Factorization have not been obtained. We do it in the present paper, using the same method as in Ref. etf1 : constructing and solving Infrared Evolution Equations (IREE)222This method was suggested by L.N. Lipatov. History of this method and its development are discussed in detail in Ref. egtg1sum .. Then we obtain expressions for $T_{q,g}$ which combine DL and non-DL contributions available in DGLAP, so these expressions can universally be used at large $Q^{2}$ and arbitrary $x$ (i.e. in the region $\textbf{A}\oplus\textbf{B}$) in Collinear, $K_{T}$ and Basic Factorizations.

Extending expressions for $T_{qg}$ to low $Q^{2}$ was suggested in Ref. egtsmallq for Collinear Factorization. In the present paper we generalize this extension to the cases of $K_{T}$ and Basic Factorizations, obtaining thereby expressions which can universally be used at arbitrary $x$ and $Q^{2}$ in any form of QCD factorization.

Putting $q^{2}=0$ in Eq. (1), one arrives at the photoproduction amplitudes $A_{\gamma}$:

[TABLE]

Combining it with Eq. (2) leads to the factorized form of $A_{\gamma q}$:

[TABLE]

The Optical theorem relates the amplitude $A_{\gamma}$ to the total cross section of photoproduction. Using the expressions for $T_{qg}$ at low $Q^{2}$ allows us to obtain expressions for the perturbative amplitudes $A_{\gamma q},A_{\gamma g}$ first in Collinear and then in the other forms of QCD factorization.

Our paper is organized as follows: in Sect. II we remind results of Ref. etf1 for amplitudes $T_{q,g}$ in Region B in Collinear Factorization and extend $T_{q,g}$ to Regions A,C,D, obtaining explicit expressions for $T_{q,g}$ which can be used at any $x$ and $Q^{2}$. Then we briefly discuss the small-$x$ asymptotics of $T_{q,g}$, their applicability region the power $Q^{2}$-corrections to $T_{q,g}$ in the low-$Q^{2}$ region C. In Sect. III we study amplitudes $T_{q,g}$ in $K_{T}$-Factorization. This type of factorization involves dealing with essentially off-shell external partons, which brings a new technical problem which is absent in Collinear Factorization: there is no universal description of $T_{q,g}$ in Region B at at arbitrary virtualities of the external partons. We deal with this problem also by constructing and solving appropriate IREEs. In Sect. IV we consider the photoproduction amplitudes in Collinear, $K_{T}$ and Basic Factorizations. Finally, Sect. V is for concluding remarks.

II Compton amplitudes in Collinear Factorization

Eq. (1) for Collinear Factorization takes the following form:

[TABLE]

where the superscript $"col"$ refers to Collinear Factorization and $\beta$ is the longitudinal fraction of momentum $k$ in Fig. 1. We start with considering amplitudes $T^{(col)}_{q,g}$ in region B and then extend these expressions to regions A,C,D.

II.1 Amplitudes $T^{(col)}_{q.g}$ in region B

Perturbative amplitudes $T^{(col)}_{q,g}$ in region B were calculated in Ref. etf1 in the double-logarithmic approximation (DLA) and the cases of fixed and running $\alpha_{s}$ were investigated separately. Amplitudes $T^{(col)}_{q,g}$ were in Ref. etf1 represented as follows:

[TABLE]

with $\mu$ being the factorization scale. The logarithmic variable $y$ in Eq. (5) is related to $Q^{2}$:

[TABLE]

and $\xi^{(+)}(\omega)$ is the positive signature factor:

[TABLE]

The signature factor $\xi^{(+)}(\omega)$ guarantees that $T_{q,g}$ are invariant with respect to permutation of the incoming and outgoing external photons. The integral representation (5) is the asymptotic form of the Sommerfeld-Watson transform but it is frequently addressed as the Mellin transform. Throughout the paper we will name $T_{q,g}$ the Mellin amplitudes. Explicit expressions for the Mellin amplitudes $F^{(col)}_{q,g}$ in the region B are:

[TABLE]

Obtained in Ref. etf1 expressions for $\Omega_{(\pm)},C_{q}^{(\pm)}$ and $C_{g}^{(\pm)}$ can be found in Appendix B where they are expressed through the on-shell parton-parton amplitudes $A_{rr^{\prime}}$ ($r,r^{\prime}=q,g$). Graphs for amplitudes $A_{rr^{\prime}}$ are depicted in Fig. 2.

The Mellin transform for the parton-parton amplitudes is similar to the one in Eq. (5):

[TABLE]

Technically, it is more convenient to use the Mellin amplitudes $h_{rr^{\prime}}$ which are proportional to $f_{rr^{\prime}}$:

[TABLE]

Explicit expression for amplitudes $h_{rr^{\prime}}$ can be found in Appendix A. Let us discuss Eq. (8). It is easy to see that its structure exhibits a distinct similarity to the structure of the DGLAP description of $F^{(col)}_{q,g}$, which is especially obvious when the approximation of fixed $\alpha_{s}$ is used:

[TABLE]

Indeed, the factors $\Omega_{(\pm)}$ in the exponents of Eq. (8) are new anomalous dimensions instead of $\gamma^{DGLAP}_{(\pm)}$ in Eq. (11) while the factors $C_{g}^{(\pm)},C_{q}^{(\pm)}$ are new coefficient functions instead of the DGLAP coefficient functions $\hat{C}^{(\pm)}_{q,g}$. However in contrast to DGLAP, the coefficient functions and the anomalous dimensions in (8) are calculated with the same means and they include the total resummations of the DL contributions. The both coefficient functions and anomalous dimensions in (8) are made of amplitudes $h_{rr^{\prime}}$. When the partonic amplitudes $h_{rr^{\prime}}$ are replaced by their Born values, the integrands in Eq. (5) coincide with the integrands of the LO DGLAP expressions in which the most singular in $\omega$ terms are retained. Further expansions of $h_{rr^{\prime}}$ and substituting them in Eq. (8) lead to the NLO (as well as to NNLO DGLAP, etc), where the most singular terms only are accounted for in each order in $\alpha_{s}$.

II.2 Extending the small-x Eqs. (5,8) to region A

Expressions in Eqs. (5,8) are defined in the region B. Now we are going to obtain an interpolation formulae for amplitudes $T^{(col)}_{q,g}$, which would reproduce the DGLAP expressions at large $x$ and at the same time would coincide with Eqs. (5,8) at small $x$. Such extension can be obtained with the four steps:

Step 1: Subtract the most singular in $\omega$ terms (i.e. the terms $\sim\alpha_{s}/\omega$ in the case of LO DGLAP; $\alpha^{2}_{s}/\omega^{3}$ in the case of NLO DGLAP, etc.) from the DGLAP anomalous dimensions $\gamma_{(\pm)}^{(DGLAP)}(\omega,\alpha_{s})$. We denote $\widetilde{\gamma}^{DGLAP}_{(\pm)}$ the result of such amputation.

Step 2: Add $\widetilde{\gamma}^{DGLAP}_{(\pm)}$ to $\Omega_{(\pm)}$. The new anomalous dimensions

[TABLE]

contain the total resummations of the double-logarithmic contributions. They are essential at small $x$ and unimportant at large $x$. At the same time, (12) contains less singular in $1/\omega^{n}$ DGLAP contributions which are dominant at large $x$.

Step 3: Subtract the LO DGLAP term $=1$ and the most singular in $\omega$ terms (i.e. the term 1 $\sim\alpha_{s}/\omega^{2}$ in the case of NLO DGLAP, $\alpha^{2}_{s}/\omega^{4}$ in the case of NNLO DGLAP, etc.) from the DGLAP expressions for the coefficient functions. We denote $\Delta\hat{C}^{(\pm)}_{q,g}$ the result of such amputation.

Step 4: Add the results obtained in Step 3 to the DL expressions $C^{(\pm)}_{q,g}$, arriving thereby at new coefficient functions $\widetilde{C}^{(\pm)}_{q,g}$:

[TABLE]

The coefficient functions $\widetilde{C}^{(\pm)}_{q,g}$ include both DL contributions and the less singular DGLAP terms. The subtractions in Steps 1,3 are necessary to avoid the double counting.

Replacing $\Omega_{(\pm)}$ and $C^{(\pm)}_{q,g}$ in Eq. (5) by $\widetilde{\Omega}_{(\pm)}$ and $\widetilde{C}^{(\pm)}_{q,g}$, we obtain the interpolation expressions $\widetilde{T}^{(col)}_{q,g}$ for the Compton amplitudes :

[TABLE]

Eq. (14) combines the small-$x$ evolution in DLA with the DGLAP results for the coefficient functions and anomalous dimensions, which are important at large $x$. These expressions can universally be used as the interpolation formulae for $T^{(col)}_{q,g}$ in the region $\textbf{A}\oplus\textbf{B}$.

II.3 Extending Eq. (14) to the small-$Q^{2}$ region $\textbf{C}\oplus\textbf{ D}$

Extension of the expressions (14) to describe the amplitudes $\widetilde{T}^{(col)}_{q,g}$ to the region $\textbf{C}\oplus\textbf{D}$ is also not rigorous. The standard approach suggested in Ref. nacht and used since that in many papers (see e.g. Ref. bad ) to describe DIS at low $Q^{2}$ is to make a shift of $Q^{2}$:

[TABLE]

where the mass scale $m$ was fixed on basis of certain physical considerations, depending on specific situation. In contrast to preceding papers, we proved in Ref. egtsmallq (see also the overviewegtg1sum ) that the scale $m$ with DL accuracy can be unambiguously specified: the IR cut-off $\mu$ plays the role of the scale $m$. Our proof was based on the well-known fact that DL contributions arise from integrations $\sim dk^{2}_{i\perp}/k^{2}_{i\perp}$ over the transverse momenta $k_{i\perp}$ of virtual partons, which requires introducing an IR cut-off $\mu$. It can be introduced through the shift $k^{2}_{i\perp}\to k^{2}_{i\perp}+\mu^{2}$, which eventually leads to the shift (15), with

[TABLE]

The proof of (16) in Refs. egtsmallq ; egtg1sum was done in the context of the spin-dependent structure function $g_{1}$ but it holds for amplitudes $\widetilde{T}^{(col)}_{q,g}$ either. Applying the Principle of Minimal Sensitivitypms , we estimated in Ref. etf1 the value of $\mu$ of Eq. (16):

[TABLE]

So, the universal expression for the Compton amplitudes $\widetilde{T}^{(col)}_{q,g}$ valid in the region $\textbf{A}\oplus\textbf{B}\oplus\textbf{C}\oplus\textbf{D}$ (i.e. at arbitrary values of $x$ and $Q^{2}$) in the framework of Collinear Factorization is

[TABLE]

where we have denoted

[TABLE]

II.4 Small-$x$ asymptotics of the Compton amplitudes and comparison with the DGLAP asymptotics

Asymptotics of $T^{(col)}_{q,g}$ at $x\to 0$ and $Q^{2}>\mu^{2}$ were considered in detail in Ref. etf1 . In the present paper we briefly remind results of Ref. etf1 and compare these asymptotics with the asymptotics of the same amplitudes obtained in the DGLAP approach. The standard mathematical tool to calculate the small-$x$ asymptotics of $T^{(col)}_{q,g}$ is Saddle Point method. Applying this method to $T^{(col)}_{q,g}$, we obtain that the small-$x$ asymptotics of amplitudes $T^{(col)}_{q,g}$ are of the Regge type. They both have the same intercepts $\omega_{0}$:

[TABLE]

The factors $\Pi_{q,g}(\omega_{0})$ in Eq. (20) are called the impact factors333Explicit expressions for the impact factors can be found un Ref. etf1 . They are the only difference between the asymptotics of $T^{(col)}_{q}$ and $T^{(col)}_{g}$. The value of the intercept $\omega_{0}$ in Eq. (20) proved to be dependent on accuracy of the calculations. The maximal intercept $\omega^{DL}_{H}$ corresponded to the roughest approximation where the quark contributions were neglected and $\alpha_{s}$ was fixed while the minimal intercept $\omega^{DL}_{S}$ corresponded to the most accurate calculation where both gluon and quark contributions were accounted for and $\alpha_{s}$ was running:

[TABLE]

As the both intercepts $>1$, the Reggeons of Eq. (20) with the intercepts (21) are the DL contributions to Pomeron, or DL Pomerons. In accordance with the conventional Pomeron terminology we call them hard (with the subscript H) and soft (with the subscript S) DL Pomerons respectively. It is interesting fact that $\omega^{DL}_{H}$ is close to the LO BFLK Pomeron intercept and $\omega^{DL}_{S}$ almost coincides with the NLO BFLK Pomeron intercept despite that BFKL have nothing in common with resummation of DL contributions. The asymptotics in Eq. (20) are given by much simpler expressions than the parent amplitudes in Eq. (18). However, the asymptotics should not be used outside their applicability region. It was shown in Ref. etf1 that the asymptotics should be used at $x<x_{max}$ and estimates of $x_{max}$ at various values of $Q^{2}$ were obtained. In particular, at $Q^{2}=10$ GeV2

[TABLE]

It was shown in Ref. etf1 that the greater is $Q^{2}$, the less is $x_{max}$. When $x\geq x_{max}$ the parent amplitudes $T_{q,g}$ should be used instead of the asymptotics.

II.5 Remark on Asymptotic Scaling

Eq. (20) can be written in such a way:

[TABLE]

with

[TABLE]

Although both $T^{(col)}_{q}$ and $T^{(col)}_{g}$ by definition depend on two independent variables $x$ and $Q^{2}$, Eq. (23) manifests that $T^{(col)}_{q,g}\ln^{3/2}(1/x)$ asymptotically depend on the single variable $\zeta$ only. We name this remarkable property the Asymptotic Scaling. This property is absent in the DGLAP approach. Indeed, DGLAP predicts that the $x$- and $Q^{2}$- dependence of $T^{(col)}_{q,g}$ at $x\to 0$ are unrelated:

[TABLE]

where $\gamma_{DGLAP}$ is the anomalous dimension and $b$ is the first coefficient of the $\beta$-function. The intercept $a$ has the phenomenological origin: it is generated by the terms $x^{-a}$ in the fits for the initial parton densities.

II.6 Perturbative power $Q^{2}$-contributions

We begin studying the power $Q^{2}$-corrections with considering $T^{(col)}_{q}$ at the large-$Q^{2}$ region, where $Q^{2}\gg\mu^{2}$ and where the logarithmic variable $\widetilde{y}$ defined in (19) can be expanded in the series as follows:

[TABLE]

Being substituted in Eq. (18), it allows us to write $T^{(col)}_{q}$ as

[TABLE]

The first term in the r.h.s. of Eq. (27) describes the Compton amplitude in the region $\textbf{A}\otimes\textbf{B}$. At large $x$ it is given by the DGLAP formulae whereas the second and third terms can wrongly be interpreted as contributions of the higher twists, probably of the non-perturbative origin though actually they are altogether perturbative:

[TABLE]

This example demonstrates that ignoring the shift (15) and interpretating all terms inversely proportional to $Q^{2}$ as an impact of higher twist could be misleading. On the other, at $Q^{2}<\mu^{2}$ the expansion (26) is replaced by another one:

[TABLE]

It leads to the small-$Q^{2}$ expression for $\widetilde{T}^{(col)}_{q}$:

[TABLE]

which decreases down to zero when $Q^{2}\to 0$. It explains very well why the terms $\sim 1/(Q^{2})^{n}$ are never seen at small $Q^{2}$ where, naively, their impact could be great. We remind that our estimate of $\mu$ is $\mu\approx 1$ GeV (see Eq. (17)). Disappearing of the terms $\sim 1/(Q^{2})^{n}$ at $Q^{2}\lesssim 1$ GeV obtained from analysis of experimental data could confirm correctness of our reasoning.

III Calculating perturbative Compton amplitudes in $K_{T}$ Factorization

Originally, $K_{T}$ Factorization was introduced in Refs. ktfact ; hefact to make possible applying BFKL for description of hadronic reactions. It implied the kinematic region of very high, or asymptotic, energies, i.e. at $x\ll 1$. However, in the cases when BFKL is not involved one can use $K_{T}$ Factorization at $x\sim 1$ as well. In the present Sect. we generalize DL expressions for $T^{(col)}_{q,g}$ to $K_{T}$ Factorization. In the framework of $K_{T}$ Factorization Eq. (1) takes the following form:

[TABLE]

where the superscript $"KT"$ refers to $K_{T}$ -Factorization, $\mu$ is the factorization scale, $\beta$ is the longitudinal fraction of the external partons (see Fig. 1) and $k^{2}_{\perp}$ stands for their transverse momentum. $\Phi^{(KT)}_{q,g}$ denote initial parton distributions and amplitudes $T^{(KT)}_{q,g}$ are perturbative. The subscripts $"q,g"$ refer to quarks and gluons in the same way as in Sect. II. We are going to calculate $T^{(KT)}_{q,g}$ in DLA. DL contributions to them arrive from the kinematic region

[TABLE]

However, amplitudes $T^{(KT)}_{q,g}$ in DLA cannot be described by a single universal expression at different virtualities $k^{2}_{\perp}$ in (32). There are two regions, where they are given by different expressions:

Region E: Moderate-virtual $k^{2}_{\perp}$, where

[TABLE]

Region F: Deeply-virtual $k^{2}_{\perp}$, where

[TABLE]

In order to avoid confusions, we will denote $T^{(\textbf{E})}_{q,g}$ and $T^{(\textbf{F})}_{q,g}$ the amplitudes $T^{(KT)}_{q,g}$ in the regions **E **and F respectively.

We will use the Mellin transform for $T^{(\textbf{E,F)}}_{q,g}$ in the following form:

[TABLE]

where $y$ is given by Eq. (6) and the new variable

[TABLE]

describes dependence of $T^{(KT)}_{q,g}$ on $k_{\perp}$ in the kinematics (32). This dependence differs $T^{(KT)}_{q,g}$ from amplitudes $T^{(col)}_{q,g}$ considered in Sect. II in Collinear Factorization. In addition to $y$ and $z$, it is convenient to introduce one more logarithmic variable:

[TABLE]

In terms of the logarithmic variables regions E and F defined in Eqs. (33,34) look as follows:

Region E: $\rho>y+z$,

Region F: $\rho<y+z$.

III.1 Calculating the off-shell Compton amplitudes in region E

IREEs for $\varphi^{(\textbf{E})}_{q,g}$ have the following form:

[TABLE]

The l.h.s. of IREEs in (38) are obtained with applying the differential operator

[TABLE]

to Eq. (35). As a result, the l.h.s.of Eq. (38) contain explicitly the derivatives over $y$ and $z$ while the factor $\omega$ corresponds to differentiation of the Mellin factor $(w/\mu^{2})^{\omega}$. The convolutions in r.h.s. of (38) are similar in structure to the convolutions in the DGLAP evolution equations. They involve amplitudes $F^{(col)}_{q}$ which were considered in Sect. II and amplitudes $H_{rr^{\prime}}$ which are still unknown. We will calculate amplitudes $H_{rr^{\prime}}$ in the next Sect. In order to specify general solutions to Eq. (38), we will use the matching

[TABLE]

with $T^{(col)}_{q,g}$ defined in Eqs. (5). The matching (40) implies the following hierarchy between $y$ and $z$:

[TABLE]

The requirement (41) is temporary. The final expressions for $T^{(\textbf{E})}_{q,g}$ will be written in the form free of this hierarchy. Solving Eq. (38) goes easier when $y$ and $z$ are replaced by new variables $\xi,\eta$:

[TABLE]

Obviously, the inverse transform is

[TABLE]

In terms of $\xi,\eta$ Eq. (38) looks simpler:

[TABLE]

Solution to Eq. (44) is:

[TABLE]

with the variables $y^{\prime},z^{\prime}$ defined as follows:

[TABLE]

In order to lift the relation between $y$ and $z$ of Eq. (41) we replace $\eta$ with $|\eta|$ everywhere save Eq. (46). Doing so and substituting Eq. (45) in (35), we arrive at expressions for amplitudes $T^{(\textbf{E})}_{q,g}$ in region E:

[TABLE]

III.2 Calculating the off-shell Compton amplitudes in region F

IREEs for $T^{(\textbf{F})}_{q,g}$ are very simple:

[TABLE]

We do not use the Mellin transform for $T^{(\textbf{F})}_{q,g}$ and because of it the l.h.s. of (48) contain derivatives over $\rho$, $y$ and $z$. Integration over momenta of all virtual partons in the deeply-virtual region F do not bring any dependence444A detailed derivation of IREEs in the case of Deeply-Virtual kinematics can be found in Ref. eit . of $T^{(\textbf{F})}_{q,g}$on $\mu$. As a result, the r.h.s. of Eq. (48) are zeros. A general solution to Eq. (48) is

[TABLE]

with $\Psi_{q}$ and $\Psi_{g}$ being arbitrary functions. In order to specify them we use the following matching: Amplitudes $T^{(\textbf{F})}_{q,g}$ coincide with amplitudes $T^{(\textbf{E})}_{q,g}$ on the border surface between regions E and F, where $\rho=y+z$:

[TABLE]

with the new amplitudes $\bar{T}^{(\textbf{E})}_{q,g}$ defined as follows:

[TABLE]

According to Eq. (49) the expression for $T^{(\textbf{F})}_{q},T^{(\textbf{F})}_{g}$ in the whole region F can be obtained from amplitudes $\bar{T}^{(\textbf{E})}_{q,g}$ by the simple change of their arguments:

[TABLE]

Both $\rho-y$ and $\rho-z$ do not depend on the IR cut-off $\mu$, so amplitudes $T^{(\textbf{F})}_{q,g}$ are IR stable.

III.3 Extension of amplitudes $T^{(\textbf{E,F})}_{q,g}$ to the region of small $Q^{2}$

Extension of $T^{(\textbf{E})}_{q,g}$ to the small-$Q^{2}$ region can be done with the shift of $Q^{2}$ like it was done in Collinear Factorization. As a result, amplitudes $T^{(\textbf{E})}_{q,g}$ can be extended to the small-$Q^{2}$ region with replacement of $x,y$ in Eq. (47) by $\widetilde{x},\widetilde{y}$ defined in Eq. (19). A similar extension of $T^{(\textbf{F})}_{q,g}$ is impossible because it would involve partons with virtualities $k^{2}_{\perp}>w$ which contradicts to Eq. (32).

III.4 Extension of amplitudes $T^{(\textbf{E,F})}_{q,g}$ to Basic Factorization

It follows from Ref. bfact that extension of Compton amplitudes $T^{(\textbf{E})}_{q,g}$ and $T^{(\textbf{F})}_{q,p}$ defined in Eqs. (47,52) to Basic factorization can be done with the very simple replacement: it is necessary and sufficient to replace $k^{2}_{\perp}$ by $|k^{2}|$ in Eqs. (47,52).

III.5 Off-shell parton-parton amplitudes $H_{rr^{\prime}}$

Expressions for Compton amplitudes $T^{(KT)}_{q,g}$ in Eqs. (47,52) include off-shell Mellin amplitudes $H_{rr^{\prime}}$. Below we calculate them in DLA. They stem from the Mellin transform for off-shell parton-parton amplitudes $\widetilde{A}_{rr^{\prime}}$:

[TABLE]

and the following definition (cf. (10)):

[TABLE]

Amplitudes $H_{rr^{\prime}}(\omega,z)$ can also be found with constructing and solving appropriate IREEs. The IREEs for the off-shell $H_{rr^{\prime}}$ are quite similar to Eq. (38):

[TABLE]

where l.h.s. of each equation corresponds to applying operator $-\mu^{2}d/d\mu^{2}$ to Eq. (53) while each r.h.s involves convolutions of $H_{rr^{\prime}}$ and $h_{rr^{\prime}}$. Specifying the general solution to Eq. (55) should be done with using the matching:

[TABLE]

Solving (55) and using the matching (56), we obtain the following expressions for $H_{rr^{\prime}}$:

[TABLE]

Explicit expressions for the terms $\Omega_{(\pm)}$ are presented in Eq. (78) while $C_{1,2}$ and $C^{\prime}_{1,2}$ are defined in Eq. (82). The overall factor $e^{-\omega z}$ in Eq. (57) converts the IR-dependent Mellin factor $(w/\mu^{2})^{\omega}$ of Eq. (53) in the IR-stable factor $\left(w/k^{2}_{\perp}\right)^{\omega}$.

IV Photoproduction amplitudes

It follows from Eqs. (1,3) that the perturbative components $A_{\gamma q},A_{\gamma g}$ of the photoproduction are related to the perturbative Compton amplitudes $T_{q,g}$ in a simple manner:

[TABLE]

Eq. (58) holds in any form of QCD factorization but expressions for the photoproduction amplitudes are different in different forms of factorizations. We start with obtaining them in Collinear Factorization.

IV.1 Photoproduction amplitudes in Collinear kinematics

According to Eq. (58), putting $q^{2}=0$ in Eqs. (5, 8) should yield $A^{(col)}_{\gamma q}$ and $A^{(col)}_{\gamma q}$ in DLA. However, such procedure cannot be done in the straightforward way because the Mellin amplitudes $F^{(col)}_{q,g}$ in (8) were obtained in Ref. etf1 for $Q^{2}\gg\mu^{2}$ only and they cannot be used at $Q^{2}<\mu^{2}$. In order to describe $A^{(col)}_{q,g}$ at low $Q^{2}$, including $Q^{2}=0$, we use the prescription in the Sect. II and replace $x,y$ in Eqs. (5, 8) by $\widetilde{x},\widetilde{y}$ defined in Eq. (19). Then, putting $Q^{2}=0$ in $\widetilde{x},\widetilde{y}$, we obtain

[TABLE]

Replacement of $\widetilde{x},\widetilde{y}$ by $x_{0},y_{0}$ in Eqs. (5, 8) allows us to obtain photoproduction amplitudes in DLA in Collinear Factorization: Combining Eqs.  (5,8) and (59), we obtain:

[TABLE]

with

[TABLE]

where $a_{\gamma q}$ is the averaged photon-quark coupling $a_{\gamma q}=\bar{e_{q}^{2}}$, $x_{0}=\mu^{2}/w$ and

[TABLE]

Accounting for non-DL radiative correction in $A^{(col)}_{\gamma q},A^{(col)}_{\gamma g}$ can be done by the same way as it was done in Sec. IIB for the Compton amplitudes.

IV.2 Photoproduction amplitudes in $K_{T}$- Factorization

DL contributions to the perturbative photoproduction amplitudes $A^{(KT)}_{\gamma q}(s,k^{2}_{\perp})$ and $A^{(KT)}_{\gamma g}(s,k^{2}_{\perp})$ come from the region

[TABLE]

We will use the Mellin transform for $A^{(KT)}_{\gamma q}$ and $A^{(KT)}_{\gamma g}$ in the following form:

[TABLE]

IREE for the Mellin amplitudes $F_{\gamma q}(\omega,z),F_{\gamma g}(\omega,z)$ are similar to Eq. (55):

[TABLE]

Expressions for on-shell amplitudes $f^{(col)}_{\gamma q}$ and $f^{(col)}_{\gamma g}$ in the r.h.s. of Eq. (65) are given by Eq. (61) while new off-shell parton-parton amplitudes $H_{rr^{\prime}}(\omega,z)$ are unknown and should be specified. The general solution to Eq. (65) should be specified by the matching:

[TABLE]

Solving Eq. (65) and using the matching of Eq. (64) yields expressions for $A^{(KT)}_{\gamma q}$ and $A^{(KT)}_{\gamma g}$:

[TABLE]

We remind that expressions for amplitudes $f^{(col)}_{\gamma q}$ and $f^{(col)}_{\gamma g}$ can are obtained in Eq. (61) while $H_{rr^{\prime}}$ are defined in Eq. (57). To conclude this Sect. we notice that non-DL corrections can be incorporated in Eqs. (60, 67) absolutely similarly to the prescription of Sect. IIB. Transition from $K_{T}$ Factorization to Basic one is done with replacement of $k^{2}_{\perp}$ by $|k^{2}|$.

V Conclusion

In this paper we have studied first the perturbative amplitudes $T_{q,g}$ of elastic Compton scattering off quarks and gluons and then perturbative components $A_{\gamma q},A_{\gamma g}$ of the photoproduction amplitudes. $T_{q,g}$ were calculated in Ref. etf1 in DLA at small $x$ and large photon virtualities $Q^{2}$ (i.e. in kinematic region B) in the framework of Collinear Factorization . We converted these results in expressions which can be used at arbitrary $x$ and $Q^{2}$ in Collinear, $K_{T}$ and Basic Factorizations. Extension to Region A was done with combining the total resummation of DL contributions and DGLAP description of $T_{q,g}$. By doing so we obtained in Eq. (14) the interpolation expressions reproducing $T_{q,g}$ in DLA at small $x$ and coinciding with the DGLAP description of $T_{q,g}$ when $x$ are not small. Using the shift of Eq. (15), we extended description $T_{q,g}$ to small $Q^{2}$. As a result, we arrived at Eq. (18) to the expressions for $T_{q,g}$ in Collinear Factorization, which can be used at arbitrary $x$ and $Q^{2}$.

In contrast to Collinear Factorization, $T_{q,g}$ are the essentially off-shell in $K_{T}$ Factorization and they cannot be described in DLA by a single expression valid at arbitrary values of virtualities $k^{2}_{\perp}$ of the external partons. It made us consider separately regions of moderate (Eq. (33)) and deep (Eq. 34) virtualities and obtain expressions Eqs. (47) and (52) for $T_{q,g}$ in those regions. We obtained them by constructing appropriate IREEs and solving them.

After transition from Region B to the to low $Q^{2}$ Region C has been studied, we became able to obtain explicit expressions (60) and (67) for the photoproduction amplitudes $A_{\gamma q},A_{\gamma g}$ in Collinear and $K_{T}$ Factorizations respectively.

Small-$x$ asymptotics of amplitudes $T_{q,g}$ are of the Regge type with the same intercept in any form of QCD factorization. Because of it we considered such asymptotics in Collinear Factorization (see Eq. (20)) and discussed dependence of its intercept on accuracy of the calculations: the higher is the accuracy, the lesser is the intercept. In other words, the hard Pomeron becomes the soft one, when the accuracy grows. We think that further increasing the accuracy can lead to vanishing supercritical Pomeron(s) and restoration of the Unitarity.

The next interesting point is that when $T_{q,g}$ are calculated in DLA, their asymptotics depends on the single variable $\zeta=Q^{2}/x^{2}$ as shown in Eq. (23). Neither DGLAP nor BFKL cause such dependence.

It is important to use the Regge asymptotics within their applicability region, i.e. at $x<x_{max}$, with $x_{max}$ given by Eq. (22). When $x>x_{max}$ the asymptotic is considerably less than the parent amplitudes, so such amplitudes should be used instead of the asymptotics. Ignoring this point leads to various misconceptions: for instance, appearance of artificial/model (hard) Pomerons and spin-dependent Pomerons.

To conclude, let us notice that DL Pomeron can play an important role for description of various hadronic reactions where Pomerons are used and the diffractive DIS in the first place, see e.g. newman .

VI Acknowledgement

We are grateful to S.I. Alekhin, D.Yu. Ivanov, G.I. Lykasov, F. Olness and O.V. Teryaev for useful communications.

Appendix A Expressions for the parton-parton amplitudes

[TABLE]

where

[TABLE]

with

[TABLE]

and

[TABLE]

Eqs. (68,70,71) express $h_{rr^{\prime}}$ through terms $b_{rr^{\prime}}$. The terms $b_{rr^{\prime}}$ include the Born factors $a_{rr^{\prime}}$ and contributions of non-ladder graphs $V_{rr^{\prime}}$:

[TABLE]

The Born factors are (see Ref. egtg1sum for detail):

[TABLE]

where $A$ and $A^{\prime}$ stand for the running QCD couplings:

[TABLE]

with $\eta=\ln\left(\mu^{2}/\Lambda^{2}_{QCD}\right)$ and $b$ being the first coefficient of the Gell-Mann- Low function. When the running effects for the QCD coupling are neglected, $A(\omega)$ and $A^{\prime}(\omega)$ are replaced by $\alpha_{s}$. The terms $V_{rr^{\prime}}$ are represented in a similar albeit more involved way (see Ref. egtg1sum for detail):

[TABLE]

with

[TABLE]

and

[TABLE]

Appendix B Explicit expressions for ingredients of Eq. (8)

The factors $C^{(\pm)}_{q,g}$ and $\Omega_{(\pm)}$ of Eq. (5 were obtained in Ref. etf1 . We list them below. All of them are expressed though the parton-parton amplitudes $h_{rr^{\prime}}$ of Appendix A. The factors $\Omega_{(\pm)}$ are:

[TABLE]

and

[TABLE]

The factors $C_{(\pm)}$ are also expressed through the parton-parton amplitudes:

[TABLE]

[TABLE]

Appendix C Expressions for the factors $C_{1,2}$ and $C^{\prime}_{1,2}$ of Eq. (57)

The terms

[TABLE]

where $R$ is given by Eq. (79).