Search

In this study, we determine the fine-tuning of parameters in $\mathcal{N} = 1$ Supersymmetric QCD, discretized on a Euclidean lattice. Our focus is on the renormalization of Yukawa (gluino-quark-squark interactions) and quartic (four-squark interactions) couplings. Given that SUSY is broken on the lattice, non-supersymmetric counterterms must be added to the discretized Lagrangian, with coefficients which must be appropriately fine-tuned in order to recover SUSY in the continuum limit. To deduce the renormalization factors and the coefficients of the counterterms, we compute the relevant three-point (for the Yukawa couplings) and four-point Green's functions (for the quartic couplings) perturbatively to one-loop and to the lowest order in lattice spacing. Both dimensional and lattice regularizations are used to implement the Modified Minimal Subtraction scheme. Our lattice formulation employs the Wilson discretization for gluino and quark fields, the Wilson gauge action for gluons, and naive discretization for squark fields. One main difficulty in this study lies in the fact that different components of squark fields mix among themselves at the quantum level. Consequently, for an appropriate fine-tuning of the couplings, these mixings must be taken into account in the renormalization conditions. All Green's functions and renormalization factors exhibit an explicit dependence on the number of colors, $N_c$, the number of flavors, $N_f$, and the gauge parameter, $\alpha$, which are left unspecified. This work follows previous investigations on SQCD and finalizes the one-loop fine-tuning of the SQCD action on the lattice, paving the way for numerical simulations of SQCD., Comment: arXiv admin note: text overlap with arXiv:2409.09089

We study the renormalization of four-quark operators in one-loop perturbation theory. We employ a coordinate-space Gauge-Invariant Renormalization Scheme (GIRS), which can be advantageous compared to other schemes, especially in nonperturbative lattice investigations. From our perturbative calculations, we extract the conversion factors between GIRS and the modified Minimal Subtraction scheme ($\overline{\rm MS}$) at the next-to-leading order. A formidable issue in the study of the four-quark operators is that operators with different Dirac matrices mix among themselves upon renormalization. We focus on both parity-conserving and parity-violating four-quark operators, which change flavor numbers by two units ($\Delta F = 2$). The extraction of the conversion factors entails the calculation of two-point Green's functions involving products of two four-quark operators, as well as three-point Green's functions with one four-quark and two bilinear operators. The significance of our results lies in their potential to refine our understanding of QCD phenomena, offering insights into the precision of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and shedding light on the nonperturbative treatment of complex mixing patterns associated with four-quark operators., Comment: 27 pages, 3 figures, 13 tables, 1 supplemental Mathematica file. Version published in Physical Review D

We present preliminary results for the renormalization functions (RFs) of a number of quark and gluon operators studied in lattice QCD using a gauge-invariant renormalization scheme (GIRS). GIRS is a variant of the coordinate-space renormalization prescription, in which Green's functions of gauge-invariant operators are calculated in position space. A novel aspect is that summations over different time slices of the positions of the operators are employed in order to reduce the statistical noise in lattice simulations. We test the reliability of this scheme by calculating RFs for the vector one-derivative quark bilinear operator, which enters the average momentum fraction of the nucleon. We use $N_f=4$ degenerate twisted mass clover-improved fermion ensembles of different volumes and lattice spacings. We also present first results of applying GIRS when operator mixing occurs: the mixing coefficients of the gluon and quark singlet energy-momentum tensor operators are evaluated by imposing appropriate renormalization conditions on the lattice., Comment: 10 pages, 5 figures, 3 tables, Talk at the 39th International Symposium on Lattice Field Theory, LATTICE2022, 8th-13th August, 2022, Rheinische Friedrich-Wilhelms-Universit\"at Bonn, Bonn, Germany

In this work, we investigate the fine tuning of parameters in $\mathcal{N} = 1$ Supersymmetric QCD, discretized on a Euclidean lattice. Specifically, we study the renormalization of the Yukawa (gluino-quark-squark interactions) and the quartic (four-squark interactions) couplings. At the quantum level, these interactions suffer from mixing with other operators which have the same transformation properties. We exploit the symmetries of the action, such as charge conjugation and parity, in order to reduce the allowed mixing patterns. To deduce the renormalizations and the mixing coefficients we compute, perturbatively to one-loop and to the lowest order in the lattice spacing, the relevant three-point and four-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. We obtain analytic expressions for the renormalization and mixing coefficients of the Yukawa couplings; they are functions of the number of colors $N_c$, the gauge parameter $\alpha$, and the gauge coupling $g$. Furthermore, preliminary results on the quartic couplings are also presented., Comment: 9 pages, 2 figures, 2 tables, Proceedings of the 39th International Symposium on Lattice Field Theory, 8th-13th August, 2022, Rheinische Friedrich-Wilhelms-Universit\"at Bonn, Bonn, Germany

We study a gauge-invariant renormalization scheme (GIRS) for composite operators, regularized on the lattice, by extending the coordinate space (X-space) scheme proposed some years ago. In this scheme, Green's functions of products of gauge-invariant operators located at different spacetime points are considered. Due to the gauge-invariant nature of GIRS, gauge fixing is not needed in the lattice simulations. Also, when operator mixing occurs, the gauge-variant operators (BRST variations and operators which vanish by the equations of motion) can be safely excluded from the renormalization process. We propose a number of variants of GIRS, including integration over time slices of the operator insertion point in a Green's function, which may lead to reduced statistical noise in lattice simulations. We employ these variants in the renormalization of fermion bilinear operators and the study of mixing between the gluon and quark energy-momentum tensor operators. We extract the one-loop conversion factors relating the nonperturbative renormalization factors in different versions of GIRS to the reference scheme of $\overline{\rm MS}$., Comment: 10 pages, 6 figures, Talk at the 38th International Symposium on Lattice Field Theory, LATTICE2021, 26th-30th July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We consider a gauge-invariant, mass-independent prescription for renormalizing composite operators, regularized on the lattice, in the spirit of the coordinate space (X-space) renormalization scheme. The prescription involves only Green's functions of products of gauge-invariant operators, situated at distinct space-time points, in a way as to avoid potential contact singularities. Such Green's functions can be computed nonperturbatively in numerical simulations, with no need to fix a gauge: thus, renormalization to this "intermediate" scheme can be carried out in a completely nonperturbative manner. Expressing renormalized operators in the $\overline{\rm MS}$ scheme requires the calculation of corresponding conversion factors. The latter can only be computed in perturbation theory, by the very nature of the $\overline{\rm MS}$; however, the computations are greatly simplified by virtue of the following attributes: i) In the absense of operator mixing, they involve only massless, two-point functions; such quantities are calculable to very high perturbative order. ii) They are gauge invariant; thus, they may be computed in a convenient gauge. iii) Where operator mixing may occur, only gauge-invariant operators will appear in the mixing pattern: Unlike other schemes, involving mixing with gauge-variant operators (which may contain ghost fields), the mixing matrices in the present scheme are greatly reduced. Still, computation of some three-point functions may not be altogether avoidable. We exemplify the procedure by computing, to lowest order, the conversion factors for fermion bilinear operators of the form $\bar\psi\Gamma\psi$ in QCD. We also employ the gauge-invariant scheme in the study of mixing between gluon and quark energy-momentum tensor operators: We compute to one loop the conversion factors relating the nonperturbative mixing matrix to the $\overline{\rm MS}$ scheme., Comment: 22 pages, 9 figures, 3 tables

We evaluate the gluon and quark contributions to the spin of the proton using an ensemble of gauge configuration generated at physical pion mass. We compute all valence and sea quark contributions to high accuracy. We perform a non-perturbative renormalization for both quark and gluon matrix elements. We find that the contribution of the up, down, strange and charm quarks to the proton intrinsic spin is $\frac{1}{2}\sum_{q=u,d,s,c}\Delta\Sigma^{q^+}=0.191(15)$ and to the total spin $\sum_{q=u,d,s,c}J^{q^+}=0.285(45)$. The gluon contribution to the spin is $J^g=0.187(46)$ yielding $J=J^q+J^g=0.473(71)$ confirming the spin sum. The momentum fraction carried by quarks in the proton is found to be $0.618(60)$ and by gluons $0.427(92)$, the sum of which gives $1.045(118)$ confirming the momentum sum rule. All scale and scheme dependent quantities are given in the $\mathrm{ \overline{MS}}$ scheme at 2 GeV.

We study $4$-dimensional SQCD with gauge group $SU(N_c)$ and $N_f$ flavors of chiral super-multiplets on the lattice. We perform extensive calculations of matrix elements and renormalization factors of composite operators in Perturbation Theory. In particular, we compute the renormalization factors of quark and squark bilinears, as well as their mixing at the quantum level with gluino and gluon bilinear operators. From these results we construct correctly renormalized composite operators, which are free of mixing effects and may be employed in non-perturbative studies of Supersymmetry. All our calculations have been performed with massive matter fields, in order to regulate the infrared singularities which are inherent in renormalizing squark bilinears. Furthermore, the quark and squark propagators are computed in momentum space with nonzero masses. This work is a feasibility study for lattice computations relevant to a number of observables, such as spectra and distribution functions of hadrons, but in the context of supersymmetric QCD, as a forerunner to lattice investigations of SUSY extensions of the Standard Model., Comment: 21 pages, 5 figures, 2 tables

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, the calculation on the lattice of which allows the direct nonperturbative study of the corresponding physical parton distribution functions. 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 additional operator mixing which occurs in the presence of nonzero quark masses, both on the lattice and in dimensional regularization., Comment: 30 pages, 4 figures. Additional comments and clarifications. Version accepted for publication in Physical Review D

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., Comment: 20 pages, 4 figures, 2 table. Refined SU(3) ChPT analysis with no changes in the final result. Version to appear in PRD

We present results for the renormalization of gauge invariant nonlocal fermion operators which contain a Wilson line, to one loop level in lattice perturbation theory. Our calculations have been performed for Wilson/clover fermions and a wide class of Symanzik improved gluon actions. The extended nature of such `long-link' operators results in a nontrivial renormalization, including contributions which diverge linearly as well as logarithmically with the lattice spacing, along with additional finite factors. We present nonperturbative prescriptions to extract the linearly divergent contributions., Comment: 8 pages, 2 figures. Talk presented at the 35th International Symposium on Lattice Field Theory, 18-24 June 2017, Granada, Spain

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., Comment: 19 pages, 10 figures

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., Comment: 6 pages, 5 figure, XII Quark Confinement proceedings

In this paper we present the perturbative evaluation of the difference between the renormalization functions of flavor singlet and nonsinglet bilinear quark operators on the lattice. The computation is performed to two loops and to lowest order in the lattice spacing, for a class of improved lattice actions, including Wilson, tree-level (TL) Symanzik and Iwasaki gluons, twisted mass and SLiNC Wilson fermions, as well as staggered fermions with twice stout-smeared links. In the staggered formalism, the stout smearing procedure is also applied to the definition of bilinear operators., Comment: 20 pages, 11 figures. The only difference between v2 and v1 is the inclusion of an Addendum in which conversion factors for different renormalization schemes are individually listed

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., Comment: 28 pages, 13 figures, 7 tables

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., Comment: 7 pages, 1 figure, 3 tables, LATTICE2014 proceedings

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., Comment: 7 pages, 3 figures, LATTICE2014 proceedings

We compute the difference in the renormalization of flavor singlet and nonsinglet fermion bilinear operators, to two loops in perturbation theory. Our results are applicable to a rather wide class of lattice actions with Symanzik improved gluons, stout links and clover fermions, including the Twisted Mass and SLiNC actions. A more detailed presentation of our results, along with relevant references, will appear in a forthcoming publication., Comment: Talk presented at the 32nd International Symposium on Lattice Field Theory, 23-28 June, 2014, Columbia University. 7 pages, 3 figures, 2 tables

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., Comment: 33 pages, 27 figures, 6 tables

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]., Comment: 7 pages, 1 figure. Talk presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

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., Comment: 7 pages, 2 figures, Lattice2013

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., Comment: 7 pages, 2 figures, LATTICE 2013

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., Comment: 18 pages, 3 figures, 3 Appendices. Version accepted for publication in Phys. Rev. D

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., Comment: 20 pages, 11 figures, 4 tables

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., Comment: 13 pages, 9 figures, talks presented by FB and GE at Xth Quark Confinement and the Hadron Spectrum, 8-12 October 2012, TUM Campus Garching, Munich, Germany

The determination of renormalization factors is of crucial importance. They relate the observables obtained on finite, discrete lattices 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. They are always present because simulations are done at lattice spacings $a$ and momenta $p$ with $ap$ not necessarily small. In this paper we try to suppress these artifacts by subtraction of one-loop contributions in lattice perturbation theory. We compare results obtained from a complete one-loop subtraction with those calculated for a subtraction of $O(a^2)$., Comment: 7 pages, 8 figures

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., Comment: Version accepted for publication in PRD. 13 pages, 5 figures

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., Comment: 36 pages, 11 figures and 6 tables. The results are included in electronic form (Mathematica files)

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., Comment: 7 pages, 4 figures, 1 table, talk given at the XXVIII International Symposium on Lattice Field Theory, Lattice2010, June 14-19, 2010, Villasimius, Italy

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., Comment: Accepted for publication in Phys. Rev. D; systematic error discussion slightly modified; total uncertainty estimate on the BK result unchanged; 42 pages, 25 figures

Perturbative and non-perturbative results on the renormalization constants of the fermion field and the twist-2 fermion bilinears are presented with emphasis on the non-perturbative evaluation of the one-derivative twist-2 vector and axial vector operators. 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 have been performed for pion masses in the range of about 450-260 MeV and at three values of the lattice spacing $a$ corresponding to $\beta=3.9, 4.05, 4.20$. Subtraction of ${\cal O}(a^2)$ terms is carried out by performing the perturbative evaluation of these operators at 1-loop and up to ${\cal O}(a^2)$. The renormalization conditions are defined in the RI$'$-MOM scheme, for both perturbative and non-perturbative results. The renormalization 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 ${\bar{\rm MS}}$ at 2 GeV using 3-loop perturbative results for the conversion factors., Comment: ETM Collaboration, 25 pages, 15 figures. Improvement on the continuum extrapolation and clarification on the estimate of the systematic errors. Typos corrected

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., Comment: 31 pages, 11 figures

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., Comment: 7 pages, 5 figures. Presented at the "XXVII International Symposium on Lattice Field Theory", July 26-31 2009, Peking University, Beijing, China

We calculate corrections to the fermion propagator and to the Green's functions of all fermion bilinear operators of the form $\bar\Psi \Gamma \Psi$, to one-loop in perturbation theory. We employ the Wilson/clover action for fermions and a family of Symanzik improved actions for gluons. The novel aspect of our calculations is that they are carried out to second order in the lattice spacing, $O(a^2)$. Consequently, they have addressed a number of new issues, most notably the appearance of loop integrands with strong IR divergences (convergent only beyond 6 dimensions). Such integrands are not present in $O(a^1)$ improvement calculations; there, IR divergent terms are seen to have the same structure as in the $O(a^0)$ case, by virtue of parity under integration, and they can thus be handled by well-known techniques. We explain how to correctly extract the full $O(a^2)$ dependence; in fact, our method is generalizable to any order in $a$. The $O(a^2)$ corrections to the quark propagator and Green's functions computed in this paper are useful to improve the nonperturbative RI-MOM determination of renormalization constants for quark bilinear operators. Our results depend on a large number of parameters: coupling constant, number of colors, lattice spacing, external momentum, clover parameter, Symanzik coefficients, gauge parameter. To make these results most easily accessible to the reader, we have included them in the distribution package of this paper, as an ASCII file named: Oa2results.m ; the file is best perused as Mathematica input., Comment: 37 pages, 3 figures, 27 tables. Results are also presented online, in a Mathematica input file. Version 2 (to appear in JHEP): Added a comment regarding Fig.3, corrected typo on Eq.A3

We compute the two-loop renormalization functions, in the RI' scheme, of local bilinear quark operators $\bar{\psi}\Gamma\psi$, where $\Gamma$ corresponds to the Vector, Axial-Vector and Tensor Dirac operators, in the lattice formulation of QCD. We consider both the flavor nonsinglet and singlet operators. 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. Finally, we present our results in the MSbar scheme, for easier comparison with calculations in the continuum. The corresponding results, for fermions in an arbitrary representation, together with some special features of superficially divergent integrals, are included in the Appendices., Comment: 42 pages, 10 figures. Version accepted in PRD. Added comments and references, provided per diagram numerical values; final results and conclusions left unchanged. This paper is a sequel to arXiv:0707.2906 (Phys. Rev. D76 (2007) 094514), which regards the scalar and pseudoscalar cases.

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., Comment: 20 pages, 8 figures, 14 tables. Phys. Rev. D77, Number 1 (2008)

This work presents the calculation of the relation between the bare coupling constant g_0 and the MSbar-renormalized coupling g_MS, g_0 = Z_g(g_0,a\mu) g_MS, to 2 loops in perturbation theory, with fermions in an arbitrary representation of the gauge group SU(N). Our calculation is performed using overlap fermions and Wilson gluons, and the background field technique has been chosen for convenience. The corresponding results in the fundamental representation appear in our longer publication [arXiv:0709.4368]. The 3-loop coefficient of the bare beta-function, b_2^L, is extracted using the 2-loop expression for Z_g, and it is presented as a function of the overlap parameter rho, the number of fermion flavors (N_f) and the number of colors (N). We also provide the expression for the ratio Lambda_L/Lambda_MS, in an arbitrary representation. A plot of Lambda_L/Lambda_MS is given in the adjoint representation., Comment: 15 pages, 5 figures

We calculate the perturbative value of the free energy in Lattice QCD, up to three loops. Our calculation is performed using Wilson gluons and the Sheikholeslami - Wolhert (clover) improved action for fermions. The free energy is directly related to the average plaquette. To carry out the calculation, we compute all relevant Feynman diagrams up to 3 loops, using a set of automated procedures in Mathematica; numerical evaluation of the resulting loop integrals is performed on finite lattice, with subsequent extrapolation to infinite size. The results are presented as a function of the fermion mass m, for any SU(N_c) gauge group, and for an arbitrary number of fermion flavors. In order to enable independent comparisons, we also provide the results on a per diagram basis, for a specific mass value., Comment: 13 pages, 5 figures, 8 tables

We calculate the relation between the bare coupling constant g_0 and the MSbar-renormalized coupling g_MS: g_0 = Z_g(g_0,a \mu) g_MS, to 2 loops in perturbation theory. We employ the standard Wilson action for gluons and the overlap action for fermions. For convenience, we have worked with the background field technique, which only requires evaluation of 2-point Green's function for the problem at hand. Our results depend explicitly on the number of fermion flavors (N_f) and colors (N). Since the dependence of Z_g on the overlap parameter rho cannot be extracted analytically, we tabulate our results for different values in the allowed range of rho (0 < rho < 2), focusing on values which are being used most frequently in simulations. Knowledge of Z_g allows us to derive the 3-loop coefficient of the bare beta-function (beta_L(g_0)) which, unlike the 1- and 2-loop coefficients, is regularization-dependent. The nontrivial dependence of Z_g and of beta_L(g_0) on rho is plotted for various choices of N, N_f., Comment: 20 pages, 6 figures, 2 Tables

We compute the ratio Lambda_L/Lambda_MS, where the scale parameter Lambda_L is associated with a lattice formulation of QCD. We consider a 3-parameter family of gluon actions, which are most frequently used for O(a) improvement a` la Symanzik. The gluon action is put togeter with standard discretizations for fermions (Wilson/clover, overlap), to provide Lambda_L for several possible combinations of fermion and gluon actions. We employ the background field technique in order to calculate the 1PI 2-point function of the background field; this leads to the coupling constant renormalization function, Z_g, at 1-loop level. Our results are obtained for an extensive range of values for the Symanzik coefficients., Comment: 11 pages, 3 figures, 3 tables

We compute the two-loop renormalization functions, in the RI' 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. A longer write-up of the present work, including the conversion of our results to the MSbar scheme and a generalization to arbitrary fermion representations, can be found in arXiv:0707.2906 ., Comment: 7 pages, 4 figures, presented at the XXV International Symposium on Lattice Field Theory, July 30 - August 4, 2007, Regensburg, Germany

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., Comment: 30 pages, 15 figures. Correction of a minor mistake in the 2-loop result for $Z_S$ and $Z_P$. The mistake affects (very slightly) Eqs.(51-54); the change in numerical values is too small to alter the ensuing plots. All conclusions remain unchanged

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., Comment: 7 pages, 6 figures, Poster presented at XXIV International Symposium on Lattice Field Theory (Lattice 2006, Tucson, Arizona)

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., Comment: 17 pages, 3 tables, 6 figures

We calculate the 1-loop renormalization of the fermion self-energy, all local fermion bilinears, as well as a set of extended bilinears which form a basis corresponding to moments of the parton distribution functions. We use the overlap action for fermions and Symanzik improved action for gluons. Our results are presented as a function of the overlap parameter rho and the parameters entering the Symanzik action., Comment: 14 pages, 10 figures, 7 tables. v2: Expanded introduction, further description of the calculation, adjusted notation, no changes in results. Version accepted in Physical Review D

We calculate the perturbative value of the Free Energy in Lattice QCD in three dimensions, up to three loops. Our calculation is performed using the Wilson formulation for gluons in SU(N) gauge theories. The Free Energy is directly related to the average plaquette. To carry out the calculation, we compute the coefficients involved in the perturbative expansion of the Free Energy up to three loops, using an automated set of procedures developed by us in Mathematica. The dependence on N is shown explicitly in our results. For purposes of comparison, we also present the individual contributions from every diagram. These have been obtained by means of two independent calculations, in order to cross check our results., Comment: 12 pages, 3 figures. Expanded introduction and discussion, more details in presentation, no changes in results. Accepted in Phys. Rev. D

We calculate the 1-loop renormalization of a set of extended fermionic bilinears which form a basis corresponding to moments of the parton distribution functions. We use the overlap action for fermions and Luescher-Weisz (LW) action for gluons. Our results are presented as a function of the overlap parameter rho and the parameters entering the LW action., Comment: 6 pages, 6 figures, presented at the conference Lattice2005(renormalization and improvement)