• In this paper, we examine the effect of nonzero quark masses on the renormalization of gauge-invariant nonlocal quark bilinear operators, including a finite-length Wilson line (called Wilson-line operators). These operators are relevant to the definition of parton quasi-distribution functions (quasi-PDFs), whose calculation on the lattice allows the direct nonperturbative study of the corresponding physical parton distribution functions (PDFs). We present our perturbative calculations of the bare Green's functions, the renormalization factors in RI' and MSbar schemes, as well as the conversion factors of these operators between the two renormalization schemes. Our computations have been performed in Dimensional Regularization at one-loop level, using massive quarks. The conversion factors can be used to convert the corresponding lattice nonperturbative results to the MSbar scheme, which is the most widely used renormalization scheme for the analysis of experimental data in high-energy physics. Also, our study is relevant for disentangling the observed operator mixing on the lattice, in the presence of nonzero quark masses.
  • We present the results of the first lattice QCD calculation of the $K \to \pi$ matrix elements of the chromomagnetic operator $O_{CM} = g\, \bar s\, \sigma_{\mu\nu} G_{\mu\nu} d$, which appears in the effective Hamiltonian describing $\Delta S = 1$ transitions in and beyond the Standard Model. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with operators of equal and lower dimensionality. The multiplicative renormalization factor as well as the mixing coefficients with the operators of equal dimension have been computed at one loop in perturbation theory. The power divergent coefficients controlling the mixing with operators of lower dimension have been determined non-perturbatively, by imposing suitable subtraction conditions. The numerical simulations have been carried out using the gauge field configurations produced by the European Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing. Our result for the B-parameter of the chromomagnetic operator at the physical pion and kaon point is $B_{CMO}^{K \pi} = 0.273 ~ (70)$, while in the SU(3) chiral limit we obtain $B_{CMO} = 0.072 ~ (22)$. Our findings are significantly smaller than the model-dependent estimate $B_{CMO} \sim 1 - 4$, currently used in phenomenological analyses, and improve the uncertainty on this important phenomenological quantity.
  • We perform a pilot study of the perturbative renormalization of a Supersymmetric gauge theory with matter fields on the lattice. As a specific example, we consider Supersymmetric ${\cal N}{=}1$ QCD (SQCD). We study the self-energies of all particles which appear in this theory, as well as the renormalization of the coupling constant. To this end we compute, perturbatively to one-loop, the relevant two-point and three-point Green's functions using both dimensional and lattice regularizations. Our lattice formulation involves the Wilson discretization for the gluino and quark fields; for gluons we employ the Wilson gauge action; for scalar fields (squarks) we use na\"ive discretization. The gauge group that we consider is $SU(N_c)$, while the number of colors, $N_c$, the number of flavors, $N_f$, and the gauge parameter, $\alpha$, are left unspecified. We obtain analytic expressions for the renormalization factors of the coupling constant ($Z_g$) and of the quark ($Z_\psi$), gluon ($Z_u$), gluino ($Z_\lambda$), squark ($Z_{A_\pm}$), and ghost ($Z_c$) fields on the lattice. We also compute the critical values of the gluino, quark and squark masses. Finally, we address the mixing which occurs among squark degrees of freedom beyond tree level: we calculate the corresponding mixing matrix which is necessary in order to disentangle the components of the squark field via an additional finite renormalization.
  • We study the self energies of all particles which appear in a lattice regularization of supersymmetric QCD (${\cal N}=1$). We compute, perturbatively to one-loop, the relevant two-point Green's functions using both the dimensional and the lattice regularizations. Our lattice formulation employs the Wilson fermion acrion for the gluino and quark fields. The gauge group that we consider is $SU(N_c)$ while the number of colors, $N_c$ and the number of flavors, $N_f$, are kept as generic parameters. We have also searched for relations among the propagators which are computed from our one-loop results. We have obtained analytic expressions for the renormalization functions of the quark field ($Z_\psi$), gluon field ($Z_u$), gluino field ($Z_\lambda$) and squark field ($Z_{A_\pm}$). We present here results from dimensional regularization, relegating to a forthcoming publication our results along with a more complete list of references. Part of the lattice study regards also the renormalization of quark bilinear operators which, unlike the non-supersymmetric case, exhibit a rich pattern of operator mixing at the quantum level.
  • We present our study of the renormalization of the chromomagnetic operator,O(CM), which appears in the effective Hamiltonian describing Delta S = 1 transitions in and beyond the Standard Model. We have computed, perturbatively to one-loop, the relevant Green's functions with two (quark-quark) and three (quark-quark-gluon) external fields, at nonzero quark masses, using both the lattice and dimensional regularizations. The perturbative computation on the lattice is carried out using the maximally twisted-mass action for the fermions, while for the gluons we employed the Symanzik improved gauge action for different sets of values of the Symanzik coefficients. We have identified all the operators which can possibly mix with O(CM), including lower dimensional and non gauge invariant operators, and we have calculated those elements of the mixing matrix which are relevant for the renormalization of O(CM). We have also performed numerical lattice calculations to determine non-perturbatively the mixings of the chromomagnetic operator with lower dimensional operators, through proper renormalization conditions. For the first time the 1/a**2-divergent mixing of the chromomagnetic operator with the scalar density has been determined non-perturbatively with high precision. Moreover, the 1/a-divergent mixing with the pseudoscalar density, due to the breaking of parity within the twisted-mass regularization of QCD, has been calculated non-perturbatively and found to be smaller than its one-loop perturbative estimate. The QCD simulations have been carried out using the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 2 + 1 + 1 dynamical quarks, which include in the sea, besides two light mass degenerate quarks, also the strange and charm quarks with masses close to their physical values.
  • The Chromomagnetic operator (CMO) mixes with a large number of operators under renormalization. We identify which operators can mix with the CMO, at the quantum level. Even in dimensional regularization (DR), which has the simplest mixing pattern, the CMO mixes with a total of 9 other operators, forming a basis of dimension-five, Lorentz scalar operators with the same flavor content as the CMO. Among them, there are also gauge noninvariant operators; these are BRST invariant and vanish by the equations of motion, as required by renormalization theory. On the other hand using a lattice regularization further operators with $d \leq 5$ will mix; choosing the lattice action in a manner as to preserve certain discrete symmetries, a minimul set of 3 additional operators (all with $d<5$) will appear. In order to compute all relevant mixing coefficients, we calculate the quark-antiquark (2-pt) and the quark-antiquark-gluon (3-pt) Green's functions of the CMO at nonzero quark masses. These calculations were performed in the continuum (dimensional regularization) and on the lattice using the maximally twisted mass fermion action and the Symanzik improved gluon action. In parallel, non-perturbative measurements of the $K-\pi$ matrix element are being performed in simulations with 4 dynamical ($N_f = 2+1+1$) twisted mass fermions and the Iwasaki improved gluon action.
  • We present preliminary results of the first lattice QCD calculation of the K -> pi matrix elements of the chromomagnetic operator O_{CM}=g sbar sigma_{munu} G_{munu} d, which appears in the effective Hamiltonian describing Delta S=1 transitions in and beyond the Standard Model. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with operators of equal and lower dimensionality. The multiplicative renormalization factor as well as the mixing coefficients with the operators of equal dimension have been computed at one-loop in perturbation theory. The power divergent coefficients controlling the mixing with operators of lower dimension have been computed non-perturbatively, by imposing suitable subtraction conditions. The numerical simulations have been carried out using the gauge field configurations produced by the European Twisted Mass Collaboration with N_f=2+1+1 dynamical quarks at three values of the lattice spacing. Our preliminary result for the B-parameter of the chromomagnetic operator is B_{CMO}=0.29(11), which can be compared with the estimate B_{CMO}~1-4 currently used in phenomenological analyses.
  • The renormalization factors of local quark-bilinear operators are computed non-perturbatively for $N_f=3$ flavors of SLiNC fermions, with emphasis on the various procedures for the chiral and continuum extrapolations. The simulations are performed at a lattice spacing $a=0.074$ fm, and for five values of the pion mass in the range of 290-465 MeV, allowing a safe and stable chiral extrapolation. Emphasis is given in the subtraction of the well-known pion pole which affects the renormalization factor of the pseudoscalar current. We also compute the inverse propagator and the Green's functions of the local bilinears to one loop in perturbation theory. We investigate lattice artifacts by computing them perturbatively to second order as well as to all orders in the lattice spacing. The renormalization conditions are defined in the RI$'$-MOM scheme, for both the perturbative and non-perturbative results. The renormalization factors, obtained at different values of the renormalization scale, are translated to the ${\bar{\rm MS}}$ scheme and are evolved perturbatively to 2 GeV. Any residual dependence on the initial renormalization scale is eliminated by an extrapolation to the continuum limit. We also study the various sources of systematic errors. Particular care is taken in correcting the non-perturbative estimates by subtracting lattice artifacts computed to one loop perturbation theory using the same action. We test two different methods, by subtracting either the ${\cal O}(g^2\,a^2)$ contributions, or the complete (all orders in $a$) one-loop lattice artifacts.
  • We calculate the fermion propagator and the quark-antiquark Green's functions for a complete set of ultralocal fermion bilinears, ${{\cal O}_\Gamma}$ [$\Gamma$: scalar (S), pseudoscalar (P), vector (V), axial (A) and tensor (T)], using perturbation theory up to one-loop and to lowest order in the lattice spacing. We employ the staggered action for fermions and the Symanzik Improved action for gluons. From our calculations we determine the renormalization functions for the quark field and for all ultralocal taste-singlet bilinear operators. The novel aspect of our calculations is that the gluon links which appear both in the fermion action and in the definition of the bilinears have been improved by applying a stout smearing procedure up to two times, iteratively. Compared to most other improved formulations of staggered fermions, the above action, as well as the HISQ action, lead to smaller taste violating effects. The renormalization functions are presented in the RI$'$ scheme; the dependence on all stout parameters, as well as on the coupling constant, the number of colors, the lattice spacing, the gauge fixing parameter and the renormalization scale, is shown explicitly. We apply our results to a nonperturbative study of the magnetic susceptibility of QCD at zero and finite temperature. In particular, we evaluate the "tensor coefficient", $\tau$, which is relevant to the anomalous magnetic moment of the muon.
  • We study matrix elements of the "chromomagnetic" operator on the lattice. This operator is contained in the strangeness-changing effective Hamiltonian which describes electroweak effects in the Standard Model and beyond. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with other operators of equal and lower dimensionality, including also non gauge invariant quantities; it is thus quite a challenge to extract from lattice simulations a clear signal for the hadronic matrix elements of this operator. We compute all relevant mixing coefficients to one loop in lattice perturbation theory; this necessitates calculating both 2-point (quark-antiquark) and 3-point (gluon-quark-antiquark) Green's functions at nonzero quark masses. We use the twisted mass lattice formulation, with Symanzik improved gluon action. For a comprehensive presentation of our results, along with detailed explanations and a more complete list of references, we refer to our forthcoming publication [1].
  • Renormalization factors relate the observables obtained on the lattice to their measured counterparts in the continuum in a suitable renormalization scheme. They have to be computed very precisely which requires a careful treatment of lattice artifacts. In this work we present a method to suppress these artifacts by subtracting one-loop contributions proportional to the square of the lattice spacing calculated in lattice perturbation theory.
  • In this paper we present the perturbative computation of the renormalization functions for the quark field and for a complete set of ultra-local fermion bilinears. The computation of the relevant Green's functions was carried out at 1-loop level for the staggered action using massive fermions. The gluon links which appear both in the fermion action and in the definition of the bilinears have been improved by applying a stout smearing procedure up to 2 times, iteratively. In the gluon sector we employed the Symanzik improved gauge action for different sets of values of the Symanzik coefficients. The renormalization functions are presented in (two variants of) the RI' and in the MSbar renormalization scheme; the dependence on all stout parameters, as well as on the fermion mass, the gauge fixing parameter and the renormalization scale, is shown explicitly. This work is related to our recent paper [Phys. Rev. D86 (2012) 094512, arXiv:1209.6015]. To make our results easily accessible to the reader, we include them in the distribution package of this paper, as a Mathematica input file, Staggered.m.
  • The determination of renormalization factors is of crucial importance in lattice QCD. They relate the observables obtained on the lattice to their measured counterparts in the continuum in a suitable renormalization scheme. Therefore, they have to be computed as precisely as possible. A widely used approach is the nonperturbative Rome-Southampton method. It requires, however, a careful treatment of lattice artifacts. In this paper we investigate a method to suppress these artifacts by subtracting one-loop contributions to renormalization factors calculated in lattice perturbation theory. We compare results obtained from a complete one-loop subtraction with those calculated for a subtraction of contributions proportional to the square of the lattice spacing.
  • We consider the effect of strong external electromagnetic fields on thermodynamic observables in QCD, through lattice simulations with 1+1+1 flavors of staggered quarks at physical quark masses. Continuum extrapolated results are presented for the light quark condensates and for their tensor polarizations, as functions of the temperature and the magnetic field. We find the light condensates to undergo inverse magnetic catalysis in the transition region, in a manner that the transition temperature decreases with growing magnetic field. We also compare the results to other approaches and lattice simulations. Furthermore, we relate the tensor polarization to the spin part of the magnetic susceptibility of the QCD vacuum, and show that this contribution is diamagnetic.
  • The response of the QCD vacuum to a constant external (electro)magnetic field is studied through the tensor polarization of the chiral condensate and the magnetic susceptibility at zero and at finite temperature. We determine these quantities using lattice configurations generated with the tree-level Symanzik improved gauge action and N_f=1+1+1 flavors of stout smeared staggered quarks with physical masses. We carry out the renormalization of the observables under study and perform the continuum limit both at T>0 and at T=0, using different lattice spacings. Finite size effects are studied by using various spatial lattice volumes. The magnetic susceptibilities \chi_f reveal a diamagnetic behavior; we obtain at zero temperature \chi_u=-(2.08 +/- 0.08) 1/GeV^2, \chi_d=-(2.02 +/- 0.09) 1/GeV^2 and \chi_s=-(3.4 +/- 1.4) 1/GeV^2 for the up, down and strange quarks, respectively, in the MSBar scheme at a renormalization scale of 2 GeV. We also find the polarization to change smoothly with the temperature in the confinement phase and then to drastically reduce around the transition region.
  • Perturbative and non-perturbative results are presented on the renormalization constants of the quark field and the vector, axial-vector, pseudoscalar, scalar and tensor currents. The perturbative computation, carried out at one-loop level and up to second order in the lattice spacing, is performed for a fermion action, which includes the clover term and the twisted mass parameter yielding results that are applicable for unimproved Wilson fermions, as well as for improved clover and twisted mass fermions. We consider ten variants of the Symanzik improved gauge action corresponding to ten different values of the plaquette coefficients. Non-perturbative results are obtained using the twisted mass Wilson fermion formulation employing two degenerate dynamical quarks and the tree-level Symanzik improved gluon action. The simulations are performed for pion masses in the range of 480 MeV to 260 MeV and at three values of the lattice spacing, a, corresponding to beta=3.9, 4.05, 4.20. For each renormalization factor computed non-perturbatively we subtract its perturbative O(a^2) terms so that we eliminate part of the cut-off artifacts. The renormalization constants are converted to MS-bar at a scale of mu=2 GeV. The perturbative results depend on a large number of parameters and are made easily accessible to the reader by including them in the distribution package of this paper, as a Mathematica input file.
  • We present an unquenched $N_f=2$ lattice computation of the $B_{K}$ parameter which controls $K^0-\bar K^0$ oscillations. A partially quenched setup is employed with two maximally twisted dynamical (sea) light Wilson quarks, and valence quarks of both the maximally twisted and the Osterwalder--Seiler variety. Suitable combinations of these two kinds of valence quarks lead to a lattice definition of the $B_{K}$ parameter which is both multiplicatively renormalizable and O($a$) improved. Employing the non-perturbative RI-MOM scheme, in the continuum limit and at the physical value of the pion mass we get $B^{\rm RGI}_K=0.729\pm 0.030$, a number well in line with the existing quenched and unquenched determinations.
  • We present perturbative and non-perturbative results on the renormalization constants of the local and one-derivative vector and axial vector operators. Non-perturbative results are obtained using the twistedmassWilson fermion formulation employing two degenerate dynamical quarks and the tree-level Symanzik improved gluon action for pion masses in the range of about 450-260 MeV and at there values of the lattice spacing, namely 0.055 fm, 0.070 fm and 0.089 fm. Subtraction of O(a^2) terms is carried out by performing the perturbative evaluation of these operators at 1- loop and up to O(a^2). The renormalization conditions are defined in the RI'-MOM scheme, for both perturbative and non-perturbative results. The Z-factors, obtained for different values of the renormalization scale, are evolved perturbatively to a reference scale set by the inverse of the lattice spacing. In addition, they are translated to MS-bar at 2 GeV using 3-loop perturbative results for the conversion factors.
  • We present results for the renormalization constants of bilinear quark operators obtained by using the tree-level Symanzik improved gauge action and the Nf=2 twisted mass fermion action at maximal twist, which guarantees automatic O(a)-improvement. Our results are also relevant for the corresponding standard (un-twisted) Wilson fermionic action since the two actions only differ, in the massless limit, by a chiral rotation of the quark fields. The scale-independent renormalization constants ZV, ZA and the ratio ZP/ZS have been computed using the RI-MOM approach, as well as other alternative methods. For ZA and ZP/ZS, the latter are based on both standard twisted mass and Osterwalder-Seiler fermions, while for ZV a Ward Identity has been used. The quark field renormalization constant Zq and the scale dependent renormalization constants ZS, ZP and ZT are determined in the RI-MOM scheme. Leading discretization effects of O(g^2 a^2), evaluated in one-loop perturbation theory, are explicitly subtracted from the RI-MOM estimates.
  • We compute the renormalization functions on the lattice, in the RI' scheme, of local bilinear quark operators $\bar{\psi}\Gamma\psi$, where $\Gamma= 1, \gamma_5, \gamma_\mu, \gamma_5\gamma_\mu, \gamma_5\sigma_{\mu\nu}$. This calculation is carried out to two loops for the first time. We consider both the flavor non-singlet and singlet operators. As a prerequisite for the above, we compute the quark field renormalization, $Z_\psi$, up to two loops. We also compute the 1-loop renormalization functions for the gluon field, $Z_A$, ghost field, $Z_c$, gauge parameter, $Z_\alpha$, and coupling constant $Z_g$. We use the clover action for fermions and the Wilson action for gluons. Our results are given as an explicit function of the coupling constant, the clover coefficient $c_{SW}$, and the number of fermion colors ($N_c$) and flavors ($N_f$), in the renormalized Feynman gauge. All 1-loop quantities are evaluated in an arbitrary gauge. Finally, we present our results in the MS-bar scheme, for easier comparison with calculations in the continuum. We have generalized to fermionic fields in an arbitrary representation. Some special features of superficially divergent integrals, obtained from the evaluation of two-loop Feynman diagrams, are presented in detail in Ref. 1.
  • We compute the two-loop renormalization functions, in the RI $^\prime$ scheme, of local bilinear quark operators $\bar{\psi}\Gamma\psi$, where $\Gamma$ denotes the Scalar and Pseudoscalar Dirac matrices, in the lattice formulation of QCD. We consider both the flavor non-singlet and singlet operators; the latter, in the scalar case, leads directly to the two-loop fermion mass renormalization, $Z_m$. As a prerequisite for the above, we also compute the quark field renormalization, $Z_{\psi}$, up to two loops. We use the clover action for fermions and the Wilson action for gluons. Our results are given as a polynomial in $c_{SW}$, in terms of both the renormalized and bare coupling constant, in the renormalized Feynman gauge. We also confirm the 1-loop renormalization functions, for generic gauge. Finally, we present our results in the $\bar{MS}$ scheme, for easier comparison with calculations in the continuum. The corresponding results, for fermions in an arbitrary representation, are included in an Appendix.
  • We calculate the critical value of the hopping parameter, $\kappa_c$, in Lattice QCD, up to two loops in perturbation theory. We employ the Sheikholeslami-Wohlert (clover) improved action for fermions and the Symanzik improved gluon action with 4- and 6-link loops. The quantity which we study is a typical case of a vacuum expectation value resulting in an additive renormalization; as such, it is characterized by a power (linear) divergence in the lattice spacing, and its calculation lies at the limits of applicability of perturbation theory. Our results are polynomial in $c_{SW}$ (clover parameter) and cover a wide range of values for the Symanzik coefficients $c_i$. The dependence on the number of colors N and the number of fermion flavors $N_f$ is shown explicitly. In order to compare our results to non perturbative evaluations of $\kappa_c$ coming from Monte Carlo simulations, we employ an improved perturbation theory method for improved actions.
  • We calculate the critical value of the hopping parameter, kappa_c, in Lattice QCD, up to two loops in perturbation theory. We employ the Sheikholeslami-Wohlert (clover) improved action for Wilson fermions and the Symanzik improved gluon action for 4- and 6-link loops. The quantity which we study is a typical case of a vacuum expectation value resulting in an additive renormalization; as such, it is characterized by a power (linear) divergence in the lattice spacing, and its calculation lies at the limits of applicability of perturbation theory. Our results are polynomial in c_{SW} (clover parameter) and cover a wide range of values for the Symanzik coefficients c_i. Furthermore, the dependence on the number of colors N and the number of fermionic flavors N_f is shown explicitly. In order to compare our results to non perturbative evaluations of kappa_c coming from Monte Carlo simulations, we employ an improved perturbation theory method applied to improved actions.
  • We study a systematic improvement of perturbation theory for gauge fields on the lattice; the improvement entails resumming, to all orders in the coupling constant, a dominant subclass of tadpole diagrams. This method, originally proposed for the Wilson gluon action, is extended here to encompass all possible gluon actions made of closed Wilson loops; any fermion action can be employed as well. The effect of resummation is to replace various parameters in the action (coupling constant, Symanzik coefficients, clover coefficient) by ``dressed'' values; the latter are solutions to certain coupled integral equations, which are easy to solve numerically. Some positive features of this method are: a) It is gauge invariant, b) it can be systematically applied to improve (to all orders) results obtained at any given order in perturbation theory, c) it does indeed absorb in the dressed parameters the bulk of tadpole contributions. Two different applications are presented: The additive renormalization of fermion masses, and the multiplicative renormalization Z_V (Z_A) of the vector (axial) current. In many cases where non-perturbative estimates of renormalization functions are also available for comparison, the agreement with improved perturbative results is significantly better as compared to results from bare perturbation theory.
  • We present a scheme for the analytic computation of renormalization functions on the lattice, using a symbolic manipulation computer language. Our first nontrivial application is a new three-loop result for the topological susceptibility.