# Perturbative renormalization of GPDs to , for various fermion/gluon actions

###### Abstract

We present a 1-loop perturbative calculation of the fermion propagator, up to (: lattice spacing). The fermions are described by Wilson, clover and twisted-mass actions; for gluons we use Symanzik improved actions (Plaquette, Tree-level Symanzik, Iwasaki, TILW, DBW2). Our results are given in a general covariant gauge, and their dependence on the coupling constant, the external momentum, the masses and the clover parameter is shown explicitly. We also study the corrections to matrix elements of unpolarized/polarized fermion bilinear operators, which include up to one derivative. These corrections are essential ingredients for improving, to , the renormalization constants of the operators under study. In addition, they can be used to minimize lattice artifacts in non-perturbative studies.

Perturbative renormalization of GPDs to , for various fermion/gluon actions

Martha Constantinou^{†}^{†}thanks: Speaker.
^{†}^{†}thanks: Work supported in part by the Research Promotion Foundation
of Cyprus (Proposal Nr: TEXN/0308/17) ,
Haralambos Panagopoulos, Fotos Stylianou

Department of Physics, University of Cyprus,

P.O.Box 20537, Nicosia CY-1678, Cyprus

E-mail: marthac@ucy.ac.cy, haris@ucy.ac.cy, fstyli01@ucy.ac.cy

\abstract@cs

PACS: 11.15.Ha, 12.38.Gc, 11.10.Gh, 12.38.Bx

## 1 Introduction

A major issue of Lattice Gauge Theory has been the reduction of effects which are due to the finite lattice spacing , in order to better approach the elusive continuum limit. Over the years, many efforts have been made for improvement in lattice observables, which in many cases is automatic by virtue of symmetries of the fermion action. According to Symanzik’s program [2], one can improve the action by a judicious addition of irrelevant operators. Also, in the twisted mass formulation of QCD [3] at maximal twist, certain observables are improved, by symmetry considerations. The first 1-loop perturbative computation of effects was recently performed by our group [4]. This regards the evaluation of the fermion propagator and Green’s functions for ultralocal bilinears of the form , using the Wilson/clover fermions and Symanzik improved gluons. Extending a calculation up to brings in new difficulties, compared to lower order in ; for instance, there appear new types of singularities. The procedure to address this issue is extensively described in Ref. [4]. The terms of such perturbative computations are of great utility since they can be subtracted from non-perturbative estimates to minimize their lattice artifacts.

The generalized parton distributions (GPDs) of the nucleon determine non-forward matrix elements of gauge invariant light cone operators; the moments of such operators can be evaluated on the lattice. GPDs also give information on interesting quantities such as the quark orbital angular momentum contribution to proton spin [5]. Moreover, the generalized parton distributions can be measured in high energy scattering experiments, for instance the deeply virtual Compton scattering of virtual photons off a nucleon. Thus, by computing moments of GPDs in lattice QCD one can explore many aspects of the nucleon structure. The operators that are related to the GPDs must be renormalized, before one compares results from simulations to physical, experimentally measurable quantities.

In this work, we investigate the perturbative renormalization of the fermion propagator, local and twist-2 fermion operators. We compute all matrix elements for the amputated Green’s functions of the inverse fermion propagator, ultralocal bilinears and twist-2 one-derivative operators, defined as (symmetrized and traceless). Although our expressions for all these matrix elements are as general as possible (the dependence on almost all parameters is shown explicitly), these are extremely lengthy and complicated to be presented here. Thus, we show our results only for the renormalization constants and for particular choices of the various parameters. In Section 2 we briefly describe the procedure for these computations and in Sections 3 - 4 we present our results for the renormalization functions. We also compare with non-perturbative estimates provided by the ETM Collaboration. Due to space limitations the renormalization of the ultralocal operators is not discussed here, but will appear in a longer write-up.

## 2 Description of the calculation

In the framework of this computation we employ a fermion action which includes the clover term, and also an additional mass parameter,

(2.0) |

A summation over the Dirac indices is implied. The advantage of the above action is that it combines different actions that are used widely: Wilson (), clover () and twisted mass (). It is important to mention that for the case of twisted mass fermions, Eq. (2) corresponds to the twisted mass action in the so called twisted basis (). The appropriate form of this action in physical basis has a rotation of the Wilson term, . This is crucial because the definition of the -factors in physical and twisted bases are different, for instance and , for the axial and vector currents with or without a derivative. The results presented here are all in the physical basis.

For the gluon part we use a family of Symanzik improved gluons and more precisely 10 different sets of Symanzik parameters that are widely used in simulations; these are shown in Table 1 of Ref. [4]. In the same reference, the reader can find all details for the calculation procedure and the evaluation of the superficially divergent integrals. The substantial difference with our previous computations is the presence of masses, and , in the fermion propagators. The dependence of superficially divergent integrals on is nontrivial and it is a rather complicated task to extract this dependence explicitly. In fact, there are 11 "strongly" IR divergent integrals (they are convergent only at dimensions) and a few hundred of divergent integrals at dimensions.

The 1-loop diagrams that enter our three calculations are shown in Fig. 1. Diagrams 1a and 1b regard the correction to the inverse fermion propagator, diagram 2 corrects the local bilinears and finally diagrams 3a-3d contribute to the calculation of the twist-2 bilinears. The appropriate operator employed in each computation is represented by a cross in diagrams 2, 3a-3d.

Figure 1: One-loop diagrams contributing to the correction of the amputated Green’s functions of the propagator (1a, 1b), local bilinears (2) and twist-2 operators (3a-3d). A wavy (solid) line represents gluons (fermions). A cross denotes an insertion of the operator under study.

## 3 Renormalization of the fermion field to

For all the renormalization factors used here we employ the RI’-MOM scheme. The fermion field renormalization constant, , can be obtained using alternative renormalization conditions that differ in their lattice artifacts. The most widely used conditions are the following two

(3.0) | |||||

(3.0) |

with the trace taken over the Dirac indices. The is the amputated Green’s function for the inverse propagator up to 1-loop; it was computed from diagrams 1a and 1b of Fig. 1. Indeed, the constants and differ only in their lattice artifacts, since the expansion of is simply .

In the most general expression for , we show explicitly the dependence on the action parameter , the coupling constant , the number of colors , the gauge fixing parameter , the masses , the lattice spacing and the external momentum . On the contrary, we cannot express in a closed form as a function of the Symanzik parameters, ; we have computed for each of the 10 sets of Symanzik coefficients shown in Table 1 of Ref. [4]. Next we provide our result for using the condition in Eq. (3), for the Landau gauge, , and tree-level Symanzik improved gluons

(3.0) | |||||

In the above expression , , and is the renormalization scale. No terms appear in Eq. (3.0) since we set equal to zero. Although the mass dependence shown in Eq. (3.0) is very complicated, the numerical result for the renormalization constant depends mildly on (for the values used in numerical simulations by ETMC, ).

A very important issue is that the terms depend not only on , but also on the direction of the momentum , as manifested by the presence of . As a consequence, different renormalization prescriptions, involving the same renormalization scale but different directions of , lead to a different renormalization constant. The strategy one may follow is to average over the results for which are obtained using these different directions.

In Table 1 we present our 1-loop results for , up to , for the two renormalization conditions: Eq. (3) () and Eq. (3) (). The results are given in the chiral limit for the Landau gauge, and for 3 values of the coupling constant . Instead of the bare coupling we have used the tadpole-improved coupling [7], defined as , to achieve further improvement. The renormalization scale was set to , that is . This is a well-defined choice, since all possible lattice momenta (permutations of ) have the same . Thus, is the same for all these permutations. and are the effects of and respectively. The last column gives the non-perturbative estimates of ETMC [6], using the renormalization condition of Eq. (3) and employing the same parameters as we did. In our corrections were already subtracted to reduce lattice artifacts. In fact, ETMC applied the subtraction procedure to all their data; therefore, the behavior of against the renormalization scale becomes flatter in the energy region where perturbation theory is valid. This indicates that the lattice artifacts are suppressed.

3.80 0 .655255014(4) -0 .039554503(1) 0 .771056601(4) 0 .076247083(1) 0 .755(4) 3.90 0 .663871997(4) -0 .045705298(1) 0 .782134880(4) 0 .072557584(1) 0 .757(3) 4.05 0 .675276827(4) -0 .053846057(1) 0 .796797308(3) 0 .067674424(1) 0 .777(5) Table 1: Perturbative and non-perturbative estimates of for different coupling constants, in the chiral limit and using two renormalization conditions, for tree-level Symanzik gluons, .The numbers is parenthesis denote the systematic errors coming from our procedures for numerical integration over loop momenta. Since we have used two alternative renormalization conditions, it is worth comparing their lattice artifacts. For and the same parameters of Table 1, we isolated the contribution of the terms and tabulated them in Table 2. In the same Table we present the averaged results, . For the averaging we used all 24 lattice momenta with : 16 permutations of and 8 permutations of . One observes that, although the condition in Eq. (3) gives correction at tree level for , the combined effects in are much smaller than the ones of , which is exactly 1 at tree level. It is also evident that averaging over momenta with different directions helps to reduce the overall contributions.

3.80 0 .35679267(9) 0 .29105315(4) 0 .45876742(9) 0 .34204053(4) 3.90 0 .26014128(8) 0 .23208860(4) 0 .41544327(8) 0 .30973960(4) 4.05 0 .13958988(7) 0 .15854326(3) 0 .36140590(7) 0 .26945127(3)Table 2: The contributions up to 1-loop in the perturbative results for the at .

## 4 corrections to the renormalization of twist-2 fermion bilinears

In this section we present the computation of the amputated Green’s functions for the following two twist-2 operators

(4.0) | |||||

(4.0) |

which are symmetrized and traceless, to avoid mixing with lower dimension operators. We have computed, to , the matrix elements of these operators for general external indices , mass, and . The final results are available for the 10 sets of Symanzik coefficients we have used in the calculation of .

The one-derivative operators fall into two different irreducible representations of the hypercubic group, depending on the choice of the external indices . Hence, we distinguish between

(4.0) |

The renormalization conditions from which and are obtained, is defined in the RI’-MOM scheme, as

(4.0) |

where . The quantity is our result up to 1-loop for the amputated Green’s function for each operator. In the renormalization condition there is a choice of using instead of . The final -factors of the two choices have different lattice artifacts. Applying the renormalization condition and using the that we computed perturbatively with its corrections, we have the results for and . In Table 3 we present these results in the chiral limit, for tree-level Symanzik gluons, , and Landau gauge. We chose six different values of the momentum, thus six renormalization scales. For we employed Eq. (3), and we used in Eq. (4.0). The extrapolation errors are smaller then the last digit shown in Table 3.

(4,2,2,2) 1.04129 1.05738 1.08270 1.06309 (5,2,2,2) 1.04411 1.04331 1.08772 1.04669 (6,2,2,2) 1.04881 1.02365 1.09312 1.02496 (3,3,3,3) 1.00343 1.05427 1.03778 1.03607 (4,3,3,3) 0.99660 1.05237 1.04456 1.03289 (5,3,3,3) 0.99731 1.04736 1.05543 1.02660 Table 3: Perturbative results for and in the chiral limit, for various renormalization scales. The action parameters are: tree-level Symanzik gluons, Landau gauge, .For the case of the vector operator with we give the result for which contributes to , for the Landau gauge, tree-level Symanzik gluons,

(4.0) | |||||

It is now interesting to use non-perturbative estimates for and , combined with our terms, in order to see if the subtraction procedure works well for these operators. Our collaborators in ETMC provided us with data for the physical [8] corresponding to the six momenta of Table 3. Note that the momenta can be grouped into two categories, according to their spatial values: (2,2,2) or (3,3,3). Non-perturbatively it is observed that these groups have different lattice artifacts and a discontinuity appears in the -factors with one derivative, as shown in Fig. 2. We also plot the subtracted non-perturbative data, which exhibit a smooth behavior of . The errors are too small to be visible.

Figure 2: Physical non-perturbative with subtractions of terms.

Another plot demonstrating the achieved improvement of non-perturbative estimates for the renormalization constants using our correction terms, is shown in Ref. [9]. In particular, for the local axial renormalization constant it is shown that, by converting the subtracted numbers into scheme and setting the renormalization scale to 2GeV, one obtains a good plateau when plotting for each of the momentum 4-vectors.

## References

- [1]
- [2] K. Symanzik, Continuum limit and improved action in lattice theories, Nucl. Phys. B226 (1983) 187; Nucl. Phys. B226 (1983) 205.
- [3] R. Frezzotti, P. Grassi, S. Sint, P. Weisz, Lattice QCD with a chirally twisted mass term, JHEP 08 (2001) 058, [hep-lat/0101001].
- [4] M. Constantinou, V. Lubicz, H. Panagopoulos, F. Stylianou, corrections to the one-loop propagator and bilinears of clover fermions with Symanzik improved gluons, JHEP, to appear, [arXiv:0907.0381].
- [5] X. Ji, Gauge-Invariant Decomposition of Nucleon Spin and Its Spin-Off, Phys. Rev. Lett. 78 (1997) 610, [hep-ph/9603249].
- [6] ETM Collaboration, in preparation.
- [7] G. Lepage and P. Mackenzie, On the Viability of Lattice Perturbation Theory, Phys. Rev. D48 (1993) 2250, [hep-lat/9209022].
- [8] C. Alexandrou, M. Constantinou, T. Korzec, in preparation.
- [9] C. Alexandrou, M. Constantinou, T. Korzec, Generalized Parton Distributions of the Nucleon from twisted mass QCD, PoS(LAT2009)136.