# Non-perturbative renormalization of tensor bilinears in Schrödinger Functional schemes^{1}^{1}1
Ift-Uam/csic-15-10, Ftuam-15-41

###### Abstract

We present preliminary result for the study of the renormalization group evolution of tensor bilinears in Schrödinger Functional (SF) schemes for and QCD with non-perturbatively -improved Wilson fermions. First results (proceeding in parallel with the ongoing computation of the running quark masses [1]) are also discussed. A one-loop perturbative calculation of the discretisation effects for the relevant step scaling functions has been carried out for both Wilson and -improved actions and for a large number of lattice resolutions. We also calculate the two-loop anomalous dimension in SF schemes for tensor currents through a scheme matching procedure with RI and . Thanks to the SF iterative procedure the non-perturbative running over two orders of magnitude in energy scales, as well as the corresponding Renormalization Group Invariant operators, have been determined.

Non-perturbative renormalization of tensor bilinears in Schrödinger Functional schemes^{†}^{†}thanks:
IFT-UAM/CSIC-15-10, FTUAM-15-41

Patrick Fritzsch, Carlos Pena, David Preti^{†}^{†}thanks: Speaker.

Instituto de Física Teórica UAM/CSIC, Universidad Autónoma de Madrid,

C/ Nicolás Cabrera 13-15, Cantoblanco, Madrid 28049

Departamento de Física Teórica, Universidad Autónoma de Madrid,

Cantoblanco, Madrid 28049

E-mail: p.fritzsch@csic.es, carlos.pena@uam.es, david.preti@csic.es

\abstract@cs

## 1 Introduction

Tensor currents play an important rôle in interesting processes through which the consistency of the Standard Model (SM) is being currently probed, as e.g. rare meson decays (see e.g. [2, 3] and references therein) or precision measurements of -decays and limits on the neutron electric dipole moment (see e.g. [4]). Considering the tensor operator for two generical (and formally distinct) flavours , as its improvement in the chiral limit is achieved by considering the combination

(1.0) |

where the coefficient was computed at 1-loop within the Schrödinger Functional (SF) in [5] and reads , for and N colours. At vanishing spatial momentum, the only non-vanishing two-point functions with boundary operators allowed by the SF boundary conditions involve the "electric" components (see Eq. (2) below)

(1.0) |

where the second term in the parenthesis of the rhs vanishes when inserted in Eq. (2). This operator renormalizes multiplicatively; i.e. the corresponding operator insertion in any on-shell renormalized correlation function is given by

(1.0) |

where , are the bare coupling and the lattice spacing respectively and is the renormalization scale. The renormalization group running is described by the Callan-Symanzik equations

(1.0) |

Since we will work in a mass independent scheme (i.e. renormalization conditions will be imposed in the chiral limit), the -function and all anomalous dimensions will only depend on the renormalized coupling and they can be expanded perturbatively in powers of as

(1.0) |

with universal coefficients , , given by

(1.0) |

All the other coefficients of the expansions are scheme dependent.

## 2 Renormalization in SF schemes

In order to impose renormalization conditions we introduce a SF two-point function of the tensor current with boundary sources of the form

(2.0) |

and the respective improved correlator is then In order to avoid extra divergences arising from the boundary the correlator can be normalized with boundary-to-boundary correlators

(2.0) |

(2.0) |

The (mass independent) renormalization conditions then read

(2.0) |

The freedom in the choice of in (2) and in the angle entering spatial boundary conditions [6] define a class of renormalization schemes. In the present work we have considered and . Following the standard SF iterative renormalization procedure [7], we define step scaling functions (SSF) as

(2.0) |

The continuum SSF for a bilinear correlator [7, 8] and for the coupling [9, 10] are defined respectively by

(2.0) |

(2.0) |

where the index has been suppressed.

## 3 Perturbative one-loop computation

We can expand all the correlators entering Eq. (2), as well as renormalization constants, in powers of

(3.0) |

where in Eq. (3) can be either , , , and . The one-loop improvement for the now reads

(3.0) |

and the renormalization constant for is given by

(3.0) |

where the notation stands for one-loop coefficients where the contribution from boundary counter terms, as well as the one related to the critical mass, have been subtracted [11]. The 1-loop critical mass is taken from [12]. Following [13] and [11], in order to study the approach to the continuum limit of the SSFs we define the relative deviation

(3.0) |

where the one-loop coefficient (see Fig. (1)) is given by

(3.0) |

The dependence on of the 1-loop renormalization constant can be described according to [14] as

(3.0) |

In order to assess the systematic uncertainties in fit coefficients, we tried several fit ansaetze by truncating the series at different orders and changing the minimal value in the fit range . We checked that, within systematic uncertainties, the correct value of the LO anomalous dimension is reproduced; after that, the term with can be subtracted to improve the fitting precision. In the errors on the fit parameters the systematics related to the choice of the ansatz has been taken into account in a conservative way. In order to extract the two-loops anomalous dimension in the SF scheme for the tensor bilinear, we used a scheme-matching procedure with a given reference scheme where is known. This bypasses a direct two-loop calculation in the SF. Since tensor bilinears renormalize multiplicatively, the same strategy employed for quark masses [13] can be used. SF NLO anomalous dimension can be written as

(3.0) |

where is related to a scheme matching for the coupling, and is related to the finite part of Eq. (3). In fact

(3.0) |

and . Since Eq. (3) is independent on the choice of the “ref” scheme, we have crosschecked our results using both and schemes ([15], [16] for finite parts, and [17] for the two-loops anomalous dimension in both schemes).

scheme s | SF () | SF () | |
---|---|---|---|

## 4 Non-perturbative renormalization and running

Our non-perturbative computation has been carried out for both and is ongoing for in parallel with the mass [1]. Once the renormalization constants given by imposing the condition (2) have been computed on a given lattice of size and the double lattice of size , eq. (2) is used to obtain the non-perturbative value of the SSF. We have computed here the SSFs for values for the coupling in the range for quenched data, and for couplings in the range for . Since we did not implement improvement of the tensor operator in the quenched case, the continuum extrapolation is linear in and has been performed using lattices with and . In the final analysis the smallest lattice has been discarded because it is affected by large cutoff effects. In order to reduce the uncertainty on the extrapolation we have performed a constrained fit between data from an -improved action and an unimproved one after testing universality of the continuum limit. Regarding , since both action and operator are -improved, the continuum limit is approached quadratically. In this case the extrapolation has been performed using and . Once in the continuum we adopted a polynomial fit ansatz for the SSFs of the form

(4.0) |

where the first two coefficients are kept fixed to their perturbative values

(4.0) |

and the cubic coefficient has been fitted for both quenched and two-flavour data.

The running between two scale and is given by eq. (1) and can be written as

(4.0) |

Once SSF has been fitted on the range of couplings, the non-perturbative running can be obtained. The evolution coefficient is computed non-perturbatively as a product of SSFs with . For both and we have computed non-perturbative steps (i.e. achieving a factor in the ratio between and ), connecting an hadronic scale ( for ) up to () an high energy scale, where perturbation theory is supposed to be safe. At those scales (computed with for and from [10, 18] respectively) the NP evolution is matched with perturbation theory at NLO, where is defined as

(4.0) |

Equivalently, the total RGI renormalization constant is defined as

(4.0) |

## 5 Conclusions

On the perturbative side we have analysed cutoff effects of the SSF for various SF schemes, and from the finite parte of the 1-loop renormalization constant we have been able to give the first preliminary determination of the NLO anomalous dimension in SF schemes for the tensor currents, which are the only quark bilinears with a non-trivial anomalous dimension independent from that of quark masses. Moreover, thanks to the non-perturbative lattice computation for both and we have computed the non-perturbative SSF in the continuum through which with 7 recursion steps the running over more than 2 orders of magnitude is computed, from an hadronic scale up to a perturbative one. Despite the dependence of the running on the scheme and on , since the correction given to the running by the NLO anomalous dimension respect to the LO is large, we still observe large systematics due to the matching with perturbation theory on the scale of -. The same strategy, briefly explained here, is being applied to simulations that will allow a more physical and precise determination of both the running and the RGI.

## References

- [1] I. Campos at al., PoS LATTICE2015 (2015) [arXiv:1508.06939 [hep-lat]].
- [2] M. Artuso et al., Eur. Phys. J. C 57, 309 (2008) [arXiv:0801.1833 [hep-ph]].
- [3] M. Antonelli et al., Phys. Rept. 494, 197 (2010) [arXiv:0907.5386 [hep-ph]].
- [4] T. Bhattacharya et al. Phys. Rev. D 85, 054512 (2012) [arXiv:1110.6448 [hep-ph]].
- [5] S. Sint and P. Weisz, Nucl. Phys. Proc. Suppl. 63, 856 (1998) [hep-lat/9709096].
- [6] S. Sint, Nucl. Phys. B 421, 135 (1994) [hep-lat/9312079].
- [7] S. Capitani et al. [ALPHA Collaboration], Nucl. Phys. B 544, 669 (1999) [hep-lat/9810063].
- [8] M. Della Morte et al. [ALPHA Collaboration], Nucl. Phys. B 729, 117 (2005) [hep-lat/0507035].
- [9] A. Bode et al. [ALPHA Collaboration], Phys. Lett. B 515, 49 (2001) [hep-lat/0105003].
- [10] M. Lüscher et al. Nucl. Phys. B 413, 481 (1994) [hep-lat/9309005].
- [11] F. Palombi, C. Pena and S. Sint, JHEP 0603, 089 (2006) [hep-lat/0505003].
- [12] H. Panagopoulos and Y. Proestos, Phys. Rev. D 65, 014511 (2002) [hep-lat/0108021].
- [13] S. Sint et al. [ALPHA Collaboration], Nucl. Phys. B 545, 529 (1999) [hep-lat/9808013].
- [14] A. Bode et al. [ALPHA Collaboration], Nucl. Phys. B 576, 517 (2000) [hep-lat/9911018].
- [15] A. Skouroupathis and H. Panagopoulos, Phys. Rev. D 79, 094508 (2009) [arXiv:0811.4264 [hep-lat]].
- [16] S. Capitani et al. Nucl. Phys. B 593, 183 (2001) [hep-lat/0007004].
- [17] J. A. Gracey, Nucl. Phys. B 667, 242 (2003) [hep-ph/0306163].
- [18] P. Fritzsch et al. Nucl. Phys. B 865, 397 (2012) [arXiv:1205.5380 [hep-lat]].