
We review recent investigations on the short distance behaviors of potentials
among baryons, which are formulated through the NambuBetheSalpeter (NBS) wave
function. After explaining the method to define the potentials, we analyze the
short distance behavior of the NBS wave functions and the corresponding
potentials by combining the operator product expansion and a renormalization
group analysis in the perturbation theory of QCD. These analytic results are
compared with numerical results obtained in lattice QCD simulations.

We investigate the cutoff effects in 2d lattice O(N) models for a variety of
lattice actions, and we identify a class of very simple actions for which the
lattice artifacts are extremely small. One action agrees with the standard
action, except that it constrains neighboring spins to a maximal relative angle
delta. We fix delta by demanding that a particular value of the step scaling
function agrees with its continuum result already on a rather coarse lattice.
Remarkably, the cutoff effects of the entire step scaling function are then
reduced to the per mille level. This also applies to the thetavacuum effects
of the step scaling function in the 2d O(3) model. The cutoff effects of other
physical observables including the renormalized coupling and the mass in the
isotensor channel are also reduced drastically. Another choice, the mixed
action, which combines the standard quadratic with an appropriately tuned large
quartic term, also has extremely small cutoff effects. The size of cutoff
effects is also investigated analytically in 1d and at N = infinity in 2d.

Lectures given at the Summer School on "Modern perspectives in lattice QCD",
Les Houches, August 328, 2009

Noncompact SO(1,N) sigmamodels are studied in terms of their large N
expansion in a lattice formulation in dimensions d \geq 2. Explicit results for
the spin and current twopoint functions as well as for the Binder cumulant are
presented to next to leading order on a finite lattice. The dynamically
generated gap is negative and serves as a couplingdependent infrared regulator
which vanishes in the limit of infinite lattice size. The cancellation of
infrared divergences in invariant correlation functions in this limit is
nontrivial and is in d=2 demonstrated by explicit computation for the above
quantities. For the Binder cumulant the thermodynamic limit is finite and is
given by 2/(N+1) in the order considered. Monte Carlo simulations suggest that
the remainder is small or zero. The potential implications for ``criticality''
and ``triviality'' of the theories in the SO(1,N) invariant sector are
discussed.

Multiparticle form factors of local operators in integrable models in two
dimensions seem to have the property that they factorize when one subset of the
particles in the external states are boosted by a large rapidity with respect
to the others. This remarkable property, which goes under the name of form
factor clustering, was first observed by Smirnov in the O(3) nonlinear
sigmamodel and has subsequently found useful applications in integrable models
without internal symmetry structure. In this paper we conjecture the nature of
form factor clustering for the general O(n) sigmamodel and make some tests in
leading orders of the 1/n expansion and for the special cases n=3,4.

The quantum field theory describing the massive O(2) nonlinear sigmamodel is
investigated through two nonperturbative constructions: The form factor
bootstrap based on integrability and the lattice formulation as the XY model.
The Smatrix, the spin and current twopoint functions, as well as the 4point
coupling are computed and critically compared in both constructions. On the
bootstrap side a new parafermionic super selection sector is found; in the
lattice theory a recent prediction for the (logarithmic) decay of lattice
artifacts is probed.

In a recent paper Stevenson claimed that analysis of the data on the wave
function renormalization constant near the critical point of the 4d Ising model
is not consistent with analytical expectations. Here we present data with
improved statistics and show that the results are indeed consistent with
conventional wisdom once one takes into account the uncertainty of lattice
artifacts in the analytical computations.

We present a strategy designed to separate several possible origins of the
wellknown enhancement of the Delta{I}=1/2 amplitude in nonleptonic kaon
decays. In particular, we seek to disentangle the contribution of physics at
the typical QCD scale (softgluon exchange) from the effects at the scale of
the charm quark mass. This is achieved by considering QCD with an unphysically
light charm quark, so that the theory possesses an approximate SU(4)_L x
SU(4)_R chiral symmetry. By computing the relevant operator matrix elements and
monitoring their values as the charm quark mass departs from the
SU(4)symmetric situation, the role of the charm quark can be assessed. We
study the influence of the charm quark mass in Chiral Perturbation Theory.
First results from lattice simulations in the SU(4)symmetric limit are also
discussed.

We report on recent work on the determination of lowenergy constants
describing Delta{S}=1 weak transitions, in order to investigate the origins of
the Delta{I}=1/2 rule. We focus on numerical techniques designed to enhance the
statistical signal in threepoint correlation functions computed with overlap
fermions near the chiral limit.

We investigate a new numerical procedure to compute fermionic correlation
functions at very small quark masses. Large statistical fluctuations, due to
the presence of local ``bumps'' in the wave functions associated with the
lowlying eigenmodes of the Dirac operator, are reduced by an exact lowmode
averaging. To demonstrate the feasibility of the technique, we compute the
twopoint correlator of the lefthanded vector current with Neuberger fermions
in the quenched approximation, for lattices with a linear extent of L~1.5 fm, a
lattice spacing a~0.09 fm, and quark masses down to the epsilonregime. By
matching the results with the corresponding (quenched) chiral perturbation
theory expressions, an estimate of (quenched) lowenergy constants can be
obtained. We find agreement between the quenched values of F extrapolated from
the pregime and extracted in the epsilonregime.

By matching 1/m^2 divergences in finitevolume twopoint correlation
functions of the scalar or pseudoscalar densities with those obtained in chiral
perturbation theory, we derive a relation between the Dirac operator zeromode
eigenfunctions at fixed nontrivial topology and the lowenergy constants of
QCD. We investigate the feasibility of using this relation to extract the pion
decay constant, by computing the zeromode correlation functions on the lattice
in the quenched approximation and comparing them with the corresponding
expressions in quenched chiral perturbation theory.

The 7particle form factors of the fundamental spin field of the O(3)
nonlinear $\sigma$model are constructed. We calculate the corresponding
contribution to the spinspin correlation function, and compare with
predictions from the spectral density scaling hypothesis. The resulting
approximation to the spinspin correlation function agrees well with that
computed in renormalized (asymptotically free) perturbation theory in the
expected energy range. Further we observe simple lower and upper bounds for the
sum of the absolute square of the form factors which may be of use for analytic
estimates.

The contributions of the improved fermion action of Sheikholeslami and
Wohlert to the two loop coefficient of the expansion of the Schroedinger
functional coupling in terms of the lattice coupling are calculated for the
gauge group SU(3). These coefficients are required for the second order
relation of lattice data to the MSbarcoupling. By taking into account all
improvement coefficients we are able to improve the Schroedinger functional to
two loop order.

The intrinsic 4point coupling, defined in terms of a truncated 4point
function at zero momentum, provides a wellestablished measure for the
interaction strength of a QFT. We show that this coupling can be computed
nonperturbatively and to high accuracy from the form factors of an
(integrable) QFT. The technique is illustrated and tested with the Ising model,
the XYmodel and the O(3) nonlinear sigmamodel. The results are compared to
those from high precision lattice simulations.

The clover coefficient $\csw$ is computed at one loop order of perturbation
theory for improved gluon actions including sixlink loops. The O(a)
improvement coefficients for the dimension three isovector composite operators
bilinear in the quark fields are also calculated.

We compute the Schroedinger functional (SF) for the case of pure SU(3) gauge
theory at twoloop order in lattice perturbation theory. This allows us to
extract the threeloop betafunction in the SFscheme. These results are
required to compute the running coupling, the Lambdaparameter and quark masses
by finite size techniques with negligible systematic errors. In addition, we
can now implement twoloop O(a) improvement in simulations and extend and study
series in alternative ("tadpoleimproved") bare couplings.

We show that the renormalization factor relating the renormalization group
invariant quark masses to the bare quark masses computed in lattice QCD can be
determined nonperturbatively. The calculation is based on an extension of a
finitesize technique previously employed to compute the running coupling in
quenched QCD. As a byproduct we obtain the $\Lambda$parameter in this theory
with completely controlled errors.

It is shown how onshell O(a) improvement can be implemented
nonperturbatively in lattice QCD with Wilson quarks. Improvement conditions
are obtained by requiring the PCAC relation to hold exactly in certain matrix
elements. These are derived from the QCD Schr\"odinger functional which enables
us to simulate directly at vanishing quark masses. In the quenched
approximation and for bare couplings in the range $0\leq g_0\leq 1$, we
determine the improved action, the improved axial current, the additive
renormalization of the quark mass and the isospin current normalization
constants Z_A and Z_V.

The universality of the continuum limit and the applicability of renormalized
perturbation theory are tested in the SU(2) lattice gauge theory by computing
two different nonperturbatively defined running couplings over a large range
of energies. The lattice data (which were generated on the powerful APE
computers at Rome II and DESY) are extrapolated to the continuum limit by
simulating sequences of lattices with decreasing spacings. Our results confirm
the expected universality at all energies to a precision of a few percent. We
find, however, that perturbation theory must be used with care when matching
different renormalized couplings at high energies.