
In this paper, we examine the effect of nonzero quark masses on the
renormalization of gaugeinvariant nonlocal quark bilinear operators, including
a finitelength Wilson line (called Wilsonline operators). These operators are
relevant to the definition of parton quasidistribution functions (quasiPDFs),
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 oneloop 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 highenergy
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 nonperturbatively, 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 Bparameter 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 modeldependent 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
selfenergies of all particles which appear in this theory, as well as the
renormalization of the coupling constant. To this end we compute,
perturbatively to oneloop, the relevant twopoint and threepoint 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 oneloop, the relevant twopoint 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 oneloop 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 nonsupersymmetric 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 oneloop, the relevant Green's functions with two
(quarkquark) and three (quarkquarkgluon) 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
twistedmass 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 nonperturbatively the mixings of the chromomagnetic operator with
lower dimensional operators, through proper renormalization conditions. For the
first time the 1/a**2divergent mixing of the chromomagnetic operator with the
scalar density has been determined nonperturbatively with high precision.
Moreover, the 1/adivergent mixing with the pseudoscalar density, due to the
breaking of parity within the twistedmass regularization of QCD, has been
calculated nonperturbatively and found to be smaller than its oneloop
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 dimensionfive, 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 quarkantiquark (2pt) and the
quarkantiquarkgluon (3pt) 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, nonperturbative
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 oneloop in perturbation theory. The power divergent
coefficients controlling the mixing with operators of lower dimension have been
computed nonperturbatively, 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 Bparameter of the chromomagnetic operator is
B_{CMO}=0.29(11), which can be compared with the estimate B_{CMO}~14 currently
used in phenomenological analyses.

The renormalization factors of local quarkbilinear operators are computed
nonperturbatively 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 290465 MeV, allowing a safe and stable chiral
extrapolation. Emphasis is given in the subtraction of the wellknown 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 nonperturbative 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 nonperturbative 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$)
oneloop lattice artifacts.

We calculate the fermion propagator and the quarkantiquark 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 oneloop 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 tastesinglet 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 strangenesschanging 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 2point
(quarkantiquark) and 3point (gluonquarkantiquark) 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 oneloop 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 ultralocal fermion
bilinears. The computation of the relevant Green's functions was carried out at
1loop 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 RomeSouthampton method. It requires, however,
a careful treatment of lattice artifacts. In this paper we investigate a method
to suppress these artifacts by subtracting oneloop contributions to
renormalization factors calculated in lattice perturbation theory. We compare
results obtained from a complete oneloop 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 treelevel 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 nonperturbative results are presented on the
renormalization constants of the quark field and the vector, axialvector,
pseudoscalar, scalar and tensor currents. The perturbative computation, carried
out at oneloop 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. Nonperturbative results are
obtained using the twisted mass Wilson fermion formulation employing two
degenerate dynamical quarks and the treelevel 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 nonperturbatively we
subtract its perturbative O(a^2) terms so that we eliminate part of the cutoff
artifacts. The renormalization constants are converted to MSbar 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 OsterwalderSeiler
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 nonperturbative RIMOM
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 nonperturbative results on the renormalization
constants of the local and onederivative vector and axial vector operators.
Nonperturbative results are obtained using the twistedmassWilson fermion
formulation employing two degenerate dynamical quarks and the treelevel
Symanzik improved gluon action for pion masses in the range of about 450260
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 nonperturbative results. The Zfactors, 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 MSbar at 2 GeV using 3loop perturbative results for the
conversion factors.

We present results for the renormalization constants of bilinear quark
operators obtained by using the treelevel 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 (untwisted) Wilson fermionic action since the two actions only
differ, in the massless limit, by a chiral rotation of the quark fields. The
scaleindependent renormalization constants ZV, ZA and the ratio ZP/ZS have
been computed using the RIMOM approach, as well as other alternative methods.
For ZA and ZP/ZS, the latter are based on both standard twisted mass and
OsterwalderSeiler 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 RIMOM scheme. Leading
discretization effects of O(g^2 a^2), evaluated in oneloop perturbation
theory, are explicitly subtracted from the RIMOM 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 nonsinglet 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 1loop 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 1loop quantities are evaluated
in an arbitrary gauge.
Finally, we present our results in the MSbar 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 twoloop Feynman diagrams, are
presented in detail in Ref. 1.

We compute the twoloop 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 nonsinglet and singlet
operators; the latter, in the scalar case, leads directly to the twoloop
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 1loop 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
SheikholeslamiWohlert (clover) improved action for fermions and the Symanzik
improved gluon action with 4 and 6link 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
SheikholeslamiWohlert (clover) improved action for Wilson fermions and the
Symanzik improved gluon action for 4 and 6link 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 nonperturbative 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 threeloop result for the topological
susceptibility.