# Double Logarithms Resummation in Exclusive Processes: The Surprising Behavior of DVCS^{§}

**T. Altinoluk**

^{1, *}**,**

**B. Pire**

^{2}**,**

**L. Szymanowski**

^{3}**,**

**S. Wallon4**

^{4}^{1}Departamento de Fisica de Particulas and IGFAE, Universidade de Santaigo de Compostela, E-15782 Santiago de Compostela, Galicia, Spain

^{4}LPT, Université Paris-Sud, CNRS, 91405, Orsay, France & UPMC Univ. Paris 06, faculté de physique, 4 place Jussieu, 75252 Paris Cedex 05, France

### Article Information

#### Identifiers and Pagination:

**Year:**2014

**Volume:**1

**Issue:**Suppl 1: M2

**First Page:**28

**Last Page:**32

**Publisher Id:**PHY-1-28

**DOI:**10.2174/1874843001401010028

#### Article History:

**Received Date:**15/6/2013

**Revision Received Date:**27/9/2013

**Acceptance Date:**25/8/2013

**Electronic publication date:**31/12/2014

**Collection year:**2014

**© Altinoluk**

*etal.*; Licensee*Bentham Open*.open-access license: This is an open access article licensed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted, non-commercial use, distribution and reproduction in any medium, provided the work is properly cited.

## Abstract

Double logarithms resummation has been much studied in inclusive as well as exclusive processes. The
Sudakov mechanism has often be the crucial tool to exponentiate potentially large contributions to amplitudes or crosssections
near phase-space boundaries. We report on a recent work where a very different pattern emerges : the DVCS
quark coefficient function C^{q} (x,ξ) develops near the particular point x =ξ a non-alternate series in ${a}_{s}^{n}\text{}{\mathrm{log}}^{2n}\text{}\left(x\u2013\xi \right)$
which may be resummed in a cosh $\left(k\sqrt{{a}_{s}}\right)\text{}\mathrm{log}\text{}\left(\left(x\u2013\xi \right)\right)$ factor. This result is at odds with the known result for the
corresponding coefficient function for the pion transition form factor near the end point C^{q} (z) although they are much
related through a *C ^{q}(Z)* correspondence.

**Keywords:**Collinear gluon radiation, DVCS, quark coefficient function, resummation, semi-eikonal approximation, soft gluon radiation..

## 1. INTRODUCTION

While perturbative calculations are widely used in quantum field theory, their summation is always a formidable task, unreachable but in the simplest, not to say most simplistic, occurences. >From elementary particle to atomic and to solid state physics, resummation techniques have been developed to go beyond a fixed order perturbation estimate through the sampling and evaluation of an infinite class of diagrams shown to dominate in a limited kinematical region for some given observables. The result of such procedures is often an exponentiated factor as in the famous Sudakov case [1], exp(*-Kg _{2}log_{2}z*) where

*g*is the coupling constant and

*z*is the ratio of two different characteristic scales, which governs both QED and QCD calculations of exclusive form factors. As explained in detail in [2] we obtain for a specific case of exclusive scattering amplitude, namely the deeply virtual Compton scattering in the generalized cosh(

*kg log z*) orken regime, a very different resummed result of the form cosh(

*kg log z*) where

*z*is a momentum fraction. To our knowledge, this form never previously emerged in field theoretic calculations. The process that we focus on is the most studied case of a class of reactions - exclusive hard hadronic processes - which are under intense experimental investigation. The result presented here provides an important stepping-stone for further developments enabling a consistent extraction of the quantities describing the 3-dimensional structure of the proton.

In the collinear factorization framework the scattering amplitude for exclusive processes such as deeply virtual Compton scattering (DVCS) has been shown [3] to factorize in specific kinematical regions, provided a large scale controls the separation of short distance dominated partonic subprocesses and long distance hadronic matrix elements, the generalized parton distributions (GPDs) [4]. The amplitude for the DVCS process

1 ${y}^{*}{\left(q\right)}^{N}\left(p\right)\text{}\to y{\left({q}^{\text{'}}\right)}^{{N}^{\text{'}}}\left({p}^{\text{'}}\right)$

with a large virtuality *q ^{2}=-Q^{2}*, factorizes in terms of perturbatively calculable coefficient functions $C\left(x,\xi ,{\alpha}_{s}\right)$ and GPDs H(x,ξ,t), where the scaling variable in the generalized Bjorken limit is the skewness ξ defined as
$\xi =\frac{{Q}^{2}}{\left(p+{p}^{\text{'}}\right).\left(q+{q}^{\text{'}}\right)}.$The calculation of first order perturbative corrections to the partonic amplitude has shown that terms of order $\frac{{\mathrm{log}}^{2}\left(x\pm \xi \right)}{x\pm \xi}$ play an important role in the region of small (x±ξ) i.e. in the vicinity of the boundary between the domains where the QCD evolution equations of GPDs take distinct forms (the so-called ERBL and DGLAP domains). We scrutinize these regions and demonstrate that they are dominated by soft fermion and gluon propagation. This explains why they can be exponentiated using quasi-eikonal techniques.

## 2. MAIN STEPS OF OUR ANALYSIS

To set up our notations, let us remind the reader of the known results for the NLO corrections to the DVCS amplitude (1), specializing to the quark contribution to its symmetric part. After proper renormalization, it reads

2 ${A}^{\mu v}\text{}=\text{}{g}_{T}^{\mu v}\text{}{\int}_{\u20131}^{1}\mathrm{dx}\text{}\left(\left(\sum _{q}^{{n}_{F}}{T}^{q}\left(x\right){H}^{q}\left(x,\xi ,t\right)\right)\right),$

where the quark coefficient function T^{q} reads [5]

(3) ${T}^{q}={C}_{0}^{q}+{C}_{1}^{q}+{C}_{\mathrm{coll}}^{q}\mathrm{log}\frac{\left({Q}^{2}\right)}{{}_{F}^{2}},$

4 ${C}_{0}^{q}={e}_{q}^{2}\left(\left(\frac{1}{x\u2013\xi +i\epsilon}\u2013x\to \u2013x\right)\right),$

(5) ${C}_{1}^{q}=\frac{{e}_{q}^{2}{}_{s}{C}_{F}}{4\left(x+i\right)}\left({\mathrm{log}}^{2}\left(\frac{x}{2}i\right)93\frac{x}{+x}\mathrm{log}\left(\frac{x}{2}i\right)\right)\left(xx\right).$

The first (resp. second) terms in Eqs. (4) and (5) correspond to the *s*-channel (resp. *u*-channel) class of diagrams. One goes from the *s*-channel to the *u-*channel by the interchange of the photon attachments. Since these two contributions are obtained from one another by a simple (x↔-x) interchange, we now restrict mostly to the discussion of the former class of diagrams.

Let us first point out that in the same spirit as for evolution equations, the extraction of the soft-collinear singularities which dominate the amplitude in the limit x → ± ξ is made easier if one uses the light-like gauge p_{1}. *A* = 0 with p_{1} = q^{'} We argue (and verified) that in this gauge the amplitude is dominated by ladder-like diagrams. We expand any momentum in the Sudakov basis *p _{1}p_{2}*, as $k=\alpha {p}_{1}+\beta {p}_{2}+{k}_{}$ where

*p*is the light-cone direction of the two incoming and outgoing partons $\left({p}_{1}^{2}={p}_{2}^{2}=0,2{p}_{1}.\text{}{p}_{2}=s={Q}^{2}2\xi \right)$. In this basis, ${q}_{\stackrel{\text{\u2217}}{y}}={p}_{1}\u20132\text{}\xi {p}_{2}.$.

_{2} We now restrict our study to the limit x → + ξ . The dominant kinematics is given by a strong ordering both in longitudinal and transverse momenta, according to (see Fig. **1**) :

$x\xi \left(\left(\beta \right)\left(\right)x\u2013\left(\xi \right)\left(\right)x\u2013\xi +\left({\beta}_{1}\right)\left(\right)\left({\beta}_{2}\right)...\right)$

6 $...\left(x\u2013\xi +{\beta}_{1}+{\beta}_{2}+...+{\beta}_{n\u20131}\left(\right)\left({\beta}_{n}\right),\right)$

7 $\left({k}_{1}^{2}\left(\left({k}_{n}^{2}\left(s\text{}{Q}^{2}\right)\right)\right)\right)$

8 $\left({a}_{1}\left(...\left({a}_{n}\left(1.\right)\right)\right)\right)$

This ordering is related to the fact that the dominant double logarithmic contribution for each loop arises from the region of phase space where both soft and collinear singularities manifest themselves. When x→ξ the left fermionic line is a hard line, from which the gluons are emitted in an eikonal way (which means that these gluons have their all four-components neglected in the vertex w.r.t. the momentum of the emitter), with an ordering in *p _{2}* direction and a collinear ordering. For the right fermionic line, eikonal approximation is not valid, since the dominant momentum flow along

*p*is from gluon to fermion, nevertheless the collinear approximation can still be applied.

_{2} When computing the coefficient functions, one faces both UV and IR divergencies. On the one hand, the UV divergencies are taken care of through renormalization, which manifest themselves by a renormalization scale *µ _{R}* dependency. On the other hand, the IR divergencies remain, but factorization proofs at any order for DVCS justify the fact that they can be absorbed inside the generalized parton distributions and result in finite coefficient functions. In our study, we are only interested into finite parts. Thus, using dimensional regularization, in a factorization scheme like $\stackrel{}{\mathrm{MS}}$, any scaleless integral can be safely put to zero although it contains both UV and IR divergencies. Following this line of thought, we can thus safely deal with DVCS on a quark for our resummation purpose.

Finally, the issue related to the iε prescription in Eq. (5) is solved by computing the coefficient function in the unphysical region ξ > 1. After analytical continuation to the physical region 0≤ ξ ≤ 1, the final result is then obtained through the shift ξ→ξ-iε.

We define *K _{n}* as the contribution of a

*n*-loop ladder to the coefficient function. Let us sketch the main steps of the derivation of

*K*and then generalize it for

_{1}*K*.

_{n}### 2.1. The Ladder Diagram at Order α_{s}

A careful analysis [2] shows that among the one loop diagrams and in the light-like axial gauge *p _{1}*.A=0, the box diagram is dominant for x→ξ. Starting from the dominant part of the numerator of the Born term which is θ=-2

*p*, the numerator of the box diagram is

_{1}9 $\mathrm{tr}\left({p}_{2}{y}^{\mu}\left(k+\left(x\u2013\xi \right)\left({p}_{2}\right)\left(k+\right)\right)\left(x+\xi \right)\left({p}_{2}\right)\left({y}^{v}\right){d}^{\mu v}.\right)$

In the limit x → ξ, while the left fermionic line is hard with a large *p _{2}* momentum, the gluonic line is soft with respect to the left fermionic line. So we perform soft gluon approximation in the numerator by taking $k+\left(x+\xi \right){p}_{2}\to \left(x+\xi \right){p}_{2}.$. The dominant contribution comes from the residue of the gluonic propagator. Thus, the numerator of the on-shell gluon propagator, d

^{µv}, is expressed in terms of transverse polarizations, i.e.

${d}^{\mu v\text{}}\u2013{\sum}_{}{}_{\left(\right)}^{\mu}\text{}{}_{\left(\right)}^{v}$

Writing the gluon polarization vectors in the light-like *P _{1}*.A=0 gauge through their Sudakov decomposition

10 ${}_{\left(\right)}^{\mu}={}_{\left(\right)}^{\mu}\u20132\text{}\frac{{}_{\left(\right).\text{}{k}_{}}}{\beta s}{p}_{1}^{\mu}$

allows us to define an effective vertex for the gluon and outgoing quark through the polarization sum

11 $\sum _{}{}_{\left(\right).\text{}}{k}_{}{}_{\left(\right)}^{\mu}\text{}=\text{}\u2013\text{}{k}_{}^{\mu}+2\text{}\frac{{k}_{}^{2}}{\beta s}{p}_{1}^{\mu}$

The numerator, (Num)_{1}, is α - independent and reads

$\frac{\u20134\left(x+\xi \right)}{\beta}\mathrm{tr}\left(\left({p}_{2}\left({k}_{}\u20132\right)\text{}\frac{{k}_{}^{2}}{\beta s}\left({p}_{1}\right)\left(k+\left(x\u2013\xi \right)\left({p}_{2}\right){p}_{1}\right)\right)\right)$

12 $=\u20134\text{}\left(x+\xi \right)s\text{}\frac{2{k}_{}^{2}}{\beta}\left(1+\frac{\left(x\u2013\xi \right)}{\beta}\right)$

We now calculate the integral over the gluon momentum *k*, using dimensional regularization

$\int {d}^{d}\text{}k\to \frac{s}{2}\int d\alpha d\beta {d}^{d\u20132}\underset{}{k},$

13 $\u20132i\frac{s}{2}{\int}_{0}^{\xi \u2013x}\frac{d\beta}{s\beta \text{}}{\int}_{0}^{}{d}^{d\u20132}\text{}\underset{\u2013}{k}\text{}{\left(\frac{{\left(\mathrm{Num}\right)}_{1}}{{L}_{1}^{2}{R}_{1}^{2}{S}^{2}}\right)}_{\alpha =\frac{{k}^{2}}{{\beta}^{s}}}$

with the denominators

${L}_{1}^{2}=-{\underset{\_}{k}}^{2}+\alpha (\beta +x+\xi )s$

${R}_{1}^{2}=-{\underset{\_}{k}}^{2}+\alpha (\beta +x-\xi )s$

${S}^{2}=-{\underset{\_}{k}}^{2}+(\beta +x-\xi )s$

${k}^{2}=-{\underset{\_}{k}}^{2}+\alpha \beta s.$

$|\beta +x-\xi |$

14 ${k}_{1\text{}}=\frac{i}{4}{e}_{q}^{2}\left(\u2013{c}_{F}{\alpha}_{s}\frac{1}{{\left(2\right)}^{2}}\right)\frac{4}{x\u2013\xi}\frac{2i}{2!}{\mathrm{log}}^{2}\left(a\left(x\u2013\xi \right)\right),$

where we kept only the most singular terms in the x→ξ region and have no control of the value of *a* within our approximation. To fix *a*, we match our approximated one-loop result with the full one-loop result (5). This amounts to cut the
${\underset{\_}{k}}^{2}$ integral at Q^{2}. The ίε term is included according to the same matching. This leads to

15 ${k}_{1\text{}}=\frac{i}{4}{e}_{q}^{2}\left(\u2013{c}_{F}{\alpha}_{s}\frac{1}{{\left(2\right)}^{2}}\right)\frac{4}{x\u2013\xi +i}\frac{2i}{2!}{\mathrm{log}}^{2}\left(\frac{\xi \u2013x}{2\xi}\text{}\u2013i\right),\text{}$

which is the known result. This is a positive test of the validity of our approximation procedure that we now generalize to the *n*-rung ladder.

### 2.2. The Ladder Diagram at Order ${\alpha}_{s}^{n}$

Let us now turn to the estimation of all ${log}^{2n}(x-\xi )$ terms in the diagram shown on Fig. (**1**). Assuming the strong ordering (7, 8) in K_{⊥} and α, the distribution of the poles generates nested integrals in β_{ί} as :

16 ${\int}_{0}^{\xi \u2013x}d{\beta}_{1}{\int}_{0}^{\xi \u2013x\u2013{\beta}_{1}}{d\beta}_{2}...{\int}_{0}^{\xi \u2013x\u2013{\beta}_{1}\u2013..\u2013{\beta}_{n\u20131}}d\text{}{\beta}_{n}$

The numerator for the n^{th} order box diagram is obtained as:

17 ${\left(\mathrm{Num}\right)}_{2}=4s{\left(x+\xi \right)}^{n}\frac{2{k}_{1}^{2}}{{\beta}_{1}}\left(1+\left(\frac{2\left(x\u2013\xi \right)}{{\beta}_{1}}\right)\frac{2{k}_{2}^{2}}{{\beta}_{2}}\text{}\left(1+\left(\frac{2\left({\beta}_{1}+x\u2013\xi \right)}{{\beta}_{2}}\right)...\frac{{2k}_{n}^{2}}{{\beta}_{n}}\right)\text{}\left(1+\left(\frac{2\left({\beta}_{n\u20131}+...+{\beta}_{1}+x\u2013\xi \right)}{{\beta}_{n}}\right)\right)\right)$

and the denominators of propagators are, for ί=1...n,

18 ${R}_{i}^{2}\mathrm{}=\u2013{\underset{}{k}}_{i}^{2}{\alpha}_{i}\mathrm{}({\beta}_{1}\dots {\beta}_{i}x\u2013\xi \mathrm{})s\text{}\mathrm{,}{S}^{2}\mathrm{}=\u2013{\underset{}{k}}_{n}^{2}\mathrm{}({\beta}_{1}\dots {\beta}_{n}x\u2013\xi \mathrm{})s\text{}\mathrm{.}$

Using dimensional regularization and omitting scaleless integrals, the integral reads:

19 $\begin{array}{c}\text{}{\int}_{0}^{\xi \u2013x}\text{}d{\beta}_{1}\text{}\dots \text{}{\int}_{0}^{\xi \u2013x\u2013\dots \u2013{\beta}_{n\u20131}}\text{}d{\beta}_{n}\text{}{\int}_{0}^{}\text{}{d}^{d\u20132}{\underset{}{k}}_{n}\text{}\dots \text{}{\int}_{0}^{{\underset{}{k}}_{2}^{2}}\text{}{d}^{d\u20132}{\underset{}{k}}_{1}{\mathrm{}(\u20131)}^{n}\\ \frac{4\text{}s{\mathrm{}\left(2i\mathrm{}\right)}^{n}}{x\u2013\xi}\frac{1}{{\beta}_{1}x\u2013\xi}\dots \frac{1}{{\beta}_{1}\text{}\text{}\dots \text{}\text{}{\beta}_{n\u20131}\text{}\text{}x\text{}\u2013\text{}\xi}\\ \frac{1}{{\underset{}{k}}_{1}^{2}}\dots \frac{1}{{\underset{}{k}}_{n}^{2}}\frac{1}{{\underset{}{k}}_{n}^{2}\u2013\mathrm{}({\beta}_{1}\dots {\beta}_{n}x\u2013\xi \mathrm{})s}\text{}\end{array}$

The integrals over ${\underset{\_}{k}}_{1}\cdots {\underset{\_}{k}}_{n}$ are performed similarly as in the one-loop case, resulting in:

20 ${k}_{n}=\frac{i}{4}{e}_{q}^{2}{\left(\u2013i{C}_{F}{a}_{s}\frac{1}{{\left(2\right)}^{2}}\right)}^{n}\frac{4}{x\u2013\xi +i}\frac{{\left(2i\right)}^{n}}{{\left(2n\right)}^{!}}{\mathrm{log}}^{2n}\left(\frac{\xi \u2013x}{2\xi}\u2013i\xi \right),$

where the matching condition introduced in one-loop case is extended to *n*-loops.

## 3. THE RESUMMED FORMULA

Based on the results Eqs. (15, 20), one can build the resummed formula for the complete amplitude; we get with

$D\mathrm{}=\sqrt{\frac{{\alpha}_{s}{C}_{F}}{2}}$

21 $\sum _{n\mathrm{}=0}^{}{K}_{n}\mathrm{}=\frac{{e}_{q}^{2}}{x\u2013\xi i}cosh\left[Dlog\left(\frac{\xi \u2013x}{2\xi}\u2013i\right)\right]\mathrm{}=\frac{1}{2}\frac{{e}_{q}^{2}}{x\u2013\xi i}\left[{\left(\frac{\xi \u2013x}{2\xi}\u2013i\right)}^{D}{\left(\frac{\xi \u2013x}{2\xi}\u2013i\right)}^{\u2013D}\right]\text{}\mathrm{.}$

In the absence of a next to leading logarithmic calculation, the minimal and most natural resummed formula which has the same *O(α _{s})* expression as the full NLO result, reads, :

22 ${\mathrm{}({C}_{0}{C}_{1}\mathrm{})}^{res}\mathrm{}=\frac{{e}_{q}^{2}}{x\u2013\xi i}\mathrm{}\&\mathrm{lcub};cosh\mathrm{}\left[Dlog\mathrm{}\right(\frac{\xi \u2013x}{2\xi}\u2013i\mathrm{}\left)\right]\u2013\frac{{D}^{2}}{2}\mathrm{}[93\frac{\xi \u2013x}{x\xi}log\mathrm{}(\frac{\xi \u2013x}{2\xi}\u2013i\mathrm{}\left)\right]\&\mathrm{rcub};\u2013\mathrm{}(x\u2013x\mathrm{})\text{}\mathrm{.}$

## 4. THE GLUON COEFFICIENT FUNCTION

For several decades the effects of gluons on many high energy processes has been widely studied. Specifically the theory of "Color Glass Condensate" shows that at very high energies the behavior of the scattering amplitudes are dominated by gluons [6]. Recently it was also shown that even at moderate energies, there are significant O(α_{s}) corrections to scattering amplitudes due to gluonic contributions for spacelike and timelike virtual Compton scatterings [7]. With the above mentioned motivations performing a similar resummation procedure for gluon coefficient function of DVCS and TCS would result in a more trustful extraction gluon GPDs.

### SUMMARY AND OUTLOOK

We have demonstrated that resummation of soft-collinear gluon radiation effects can be performed in hard exclusive reactions amplitudes. The resulting formula for coefficient function stabilizes the perturbative expansion, which is crucial for a trustful extraction of GPDs from experimental data. A related expression should emerge in various reactions, such as the crossed case of timelike Compton scattering [8] and exclusive meson electroproduction.

Giving these results, a question should be raised : what is the physics beyond this result, or in other words, why is the Sudakov resummation [9] familiar to experts of hard exclusive processes not applicable here? An even more precise question may be: how is our analysis compatible with the discussion of soft effects in the pion transition form factor, a quantity which has been much discussed [10] recently thanks to the experimental results of BABAR and BELLE? Let us stress that the coefficient function of this quantity is identical to the ERBL part of the coefficient function of the DVCS amplitude after a rescaling →x/ξ. Our result thus may be applied to the transition form factor. In [11], it has been argued that the ${\alpha}_{s}{log}^{2}\mathrm{}(1\u2013z\mathrm{})$ factor in the one loop expression of the coefficient function had to be understood as the sum of two very distinct terms, one of them exponentiating in a Sudakov form factor. To advocate this fact, the authors allow themselves an excursion outside the colinear factorization framework and use the familiar detour into the coordinate space framework. Our procedure is different and we resum the complete one loop result. In other words, one may ask to the authors of [11]: what happens to the remnant term proportional to ${\alpha}_{s}{log}^{2}\mathrm{}(1\u2013z\mathrm{})$ If indeed the usual resummation procedure of the transition pion form factor must be revised following our new results, one may ask whether the understanding of the meson form factor [12] should also be reconsidered.

## CONFLICT OF INTEREST

The authors confirm that this article content has no conflict of interest.

## ACKNOWLEDGEMENTS

We thank the organisors and the French CEA (IPhT and DSM) for support. This work is supported by the P2IO consortium, the Polish Grant NCN No. DEC-2011/01/B/ ST2/03915, the French grant grant ANR PARTONS (ANR-12-MONU-0008-01), the Joint Research Activity "Study of Strongly Interacting Matter" (HadronPhysics3, Grant Agree-ment no. 283286) under the 7^{th} Framework Programme of the European Community, the European Research Council grant HotLHC ERC-2001- StG-279579, Ministerio de Ciencia e Innovac on of Spain grants FPA2009-06867-E, Consolider-Ingenio 2010 CPAN CSD2007-00042 and FEDER.