
It is shown that the algorithm introduced in [1] and conceived to deal with
continuous degrees of freedom models is well suited to compute the density of
states in models with a discrete energy spectrum too. The q=10 D=2 Potts model
is considered as a test case, and it is shown that using the Maxwell
construction the interface free energy can be obtained, in the thermodynamic
limit, with a good degree of accuracy.

We define a family of Schroedinger Functional renormalization schemes for the
fourquark multiplicatively renormalizable operators of the $\Delta F = 1$ and
$\Delta F = 2$ effective weak Hamiltonians. Using the lattice regularization
with quenched Wilson quarks, we compute nonperturbatively the renormalization
group running of these operators in the continuum limit in a large range of
renormalization scales. Continuum limit extrapolations are well controlled
thanks to the implementation of two fermionic actions (Wilson and Clover). The
ratio of the renormalization group invariant operator to its renormalized
counterpart at a low energy scale, as well as the renormalization constant at
this scale, is obtained for all schemes.

We give a continuum limit value of the lowest moment of a twist2 operator in
pion states from nonperturbative lattice calculations. We find that the
nonperturbatively obtained renormalization group invariant matrix element is
<x>_{RGI} = 0.179(11), which corresponds to <x>^{MSbar}(2 GeV) = 0.246(15). In
obtaining the renormalization group invariant matrix element, we have
controlled important systematic errors that appear in typical lattice
simulations, such as nonperturbative renormalization, finite size effects and
effects of a nonvanishing lattice spacing. The crucial limitation of our
calculation is the use of the quenched approximation. Another question that
remains not fully clarified is the chiral extrapolation of the numerical data.

We investigate finite size effects of the pion matrix element of the
nonsinglet, twist2 operator corresponding to the average momentum of
nonsinglet quark densities. Using the quenched approximation, they come out to
be surprisingly large when compared to the finite size effects of the pion
mass. As a consequence, simulations of corresponding nucleon matrix elements
could be affected by finite size effects even stronger which could lead to
serious systematic uncertainties in their evaluation.

The renormalisation group running of the quark mass is determined
nonperturbatively for a large range of scales, by computing the step scaling
function in the Schroedinger Functional formalism of quenched lattice QCD both
with and without O(a) improvement. A oneloop perturbative calculation of the
discretisation effects has been carried out for both the Wilson and the
Cloverimproved actions and for a large number of lattice resolutions. The
nonperturbative computation yields continuum results which are regularisation
independent, thus providing convincing evidence for the uniqueness of the
continuum limit. As a byproduct, the ratio of the renormalisation group
invariant quark mass to the quark mass, renormalised at a hadronic scale, is
obtained with very high accuracy.

We present the current status of the apeNEXT project. Aim of this project is
the development of the next generation of APE machines which will provide
multiteraflop computing power. Like previous machines, apeNEXT is based on a
custom designed processor, which is specifically optimized for simulating QCD.
We discuss the machine design, report on benchmarks, and give an overview on
the status of the software development.

We compute charm and bottom quark masses in the quenched approximation and in
the continuum limit of lattice QCD. We make use of a step scaling method,
previously introduced to deal with two scale problems, that allows to take the
continuum limit of the lattice data. We determine the RGI quark masses and make
the connection to the MSbar scheme. The continuum extrapolation gives us a
value m_b^{RGI} = 6.73(16) GeV for the bquark and m_c^{RGI} = 1.681(36) GeV
for the cquark, corresponding respectively to m_b^{MSbar}(m_b^{MSbar}) =
4.33(10) GeV and m_c^{MSbar}(m_c^{MSbar}) = 1.319(28) GeV. The latter result,
in agreement with current estimates, is for us a check of the method. Using our
results on the heavy quark masses we compute the mass of the Bc meson, M_{Bc} =
6.46(15) GeV.

We compute nonperturbatively the evolution of the twist2 operators
corresponding to the average momentum of nonsinglet quark densities. The
calculation is based on a finitesize technique, using the Schr\"odinger
Functional, in quenched QCD. We find that a careful choice of the boundary
conditions, is essential, for such operators, to render possible the
computation. As a byproduct we apply the nonperturbatively computed
renormalization constants to available data of bare matrix elements between
nucleon states.

Some new results on nonperturbative renormalisation of quark bilinears in
quenched QCD with Schroedinger Functional techniques are presented. Special
emphasis is put on a study of the universality of the continuum limit for step
scaling functions computed with different levels of O(a) improvement.

We report on recent results for the pion matrix element of the twist2
operator corresponding to the average momentum of nonsinglet quark densities.
For the first time finite volume effects of this matrix element are
investigated and come out to be surprisingly large. We use standard Wilson and
nonperturbatively improved clover actions in order to control better the
extrapolation to the continuum limit. Moreover, we compute, fully
nonperturbatively, the renormalization group invariant matrix element, which
allows a comparison with experimental results in a broad range of energy
scales. Finally, we discuss the remaining uncertainties, the extrapolation to
the chiral limit and the quenched approximation.

We compute the decay constants for the heavylight pseudoscalar mesons in
the quenched approximation and continuum limit of lattice QCD. Within the
Schrodinger Functional framework, we make use of the step scaling method, which
has been previously introduced in order to deal with the two scale problem
represented by the coexistence of a light and a heavy quark. The continuum
extrapolation gives us a value $f_{B_s} = 192(6)(4)$ MeV for the $B_s$ meson
decay constant and $f_{D_s} = 240(5)(5)$ MeV for the $D_s$ meson.

We discuss the usage of continuous external momenta for computing
renormalization factors as needed to renormalize operator matrix elements.
These kind of external momenta are encoded in special boundary conditions for
the fermion fields. The method allows to compute certain renormalization
factors on the lattice that would have been very difficult, if not impossible,
to compute with standard methods. As a result we give the renormalization group
invariant step scaling function for a twist2 operator corresponding to the
average momentum of nonsinglet quark densities.

A novel method to calculate f_B on the lattice is introduced, based on the
study of the dependence of finite size effects upon the heavy quark mass of
flavoured mesons and on a nonperturbative recursive finite size technique. We
avoid the systematic errors related to extrapolations from the static limit or
to the tuning of the coefficients of effective Lagrangian and the results admit
an extrapolation to the continuum limit. We perform a first estimate at finite
lattice spacing, but close to the continuum limit, giving f_B = 170(11)(5)(22)
MeV. We also obtain f_{B_s} = 192(9)(5)(24) MeV. The first error is
statistical, the second is our estimate of the systematic error from the method
and the third the systematic error from the specific approximations adopted in
this first exploratory calculation. The method can be generalized to twoscale
problems in lattice QCD.

We discuss a novel method to calculate fB on the lattice based on the study
of the dependence of finite size effects upon the heavy quark mass of flavoured
mesons and on a nonperturbative recursive finite size technique. This method
avoids the systematic errors related to extrapolations from the static limit or
to the tuning of the coefficients of effective Lagrangian and the results admit
an extrapolation to the continuum limit. We show the results of a first
estimate at finite lattice spacing, but close to the continuum limit, giving fB
= 170(11)(5)(22)MeV. We also obtain fBs = 192(9)(5)(24)MeV. The first error is
statistical, the second is our estimate of the systematic error from the method
and the third the systematic error from the specific approximations adopted in
this first exploratory calculation. The method can be generalized to twoscale
problems in lattice QCD.

In this paper we extend the estimate of the value of the lattice spacing a in
units of the r0 scale at values of the bare coupling larger than those
available. By using results from the computation of the renormalised coupling
in the Schroedinger functional formalism we find that from beta \simeq 7 onward
the behaviour predicted by asymptotic freedom at tree loop describes very well
the data if a value of (r0 Lambda) slightly lower then the latest available is
used. We also show that, by sticking to the current value, an effective four
loop term can describe as well the data. The systematic relative error on the
lattice spacing induced by the choice of the procedure is between 1% and 3% for
6.92 < beta < 8.5.

We apply the Schroedinger Functional (SF) formalism to determine the
renormalisation group running of fourfermion operators which appear in the
effective weak Hamiltonian of the Standard Model. Our calculations are done
using Wilson fermions and the parityodd components of the operators.
Preliminary results are presented for the operator $O_{VA}=(\bar s \gamma_\mu
d)(\bar s \gamma_\mu \gamma_5 d)$.

APENEXT is a new generation APE processor, optimized for LGT simulations. The
project follows the basic ideas of previous APE machines and develops simple
and cheap parallel systems with multi TFlops processing power. This paper
describes the main features of this new development.

We present the results of a nonperturbative determination of the pion matrix
element of the twist2 operator corresponding to the average momentum of
nonsinglet quark densities. The calculation is made within the Schroedinger
functional scheme. We report the results of simulations done with the standard
Wilson action and with the nonperturbatively improved clover action and we
show that their ratio correctly extrapolates, in the continuum limit, to a
value compatible with the residual correction factor expected from perturbation
theory.

We compute, within the Schr\"odinger functional scheme, a renormalization
group invariant renormalization constant for the first moment of the
nonsinglet parton distribution function. The matching of the results of our
nonperturbative calculation with the ones from hadronic matrix elements allows
us to obtain eventually a renormalization group invariant average momentum of
nonsinglet parton densities, which can be translated into a preferred scheme
at a specific scale.

We present evidence for the universality of the continuum limit of the scale
dependence of the renormalization constant associated with the operator
corresponding to the average momentum of nonsinglet parton densities. The
evidence is provided by a nonperturbative computation in quenched lattice QCD
using the Schr\"odinger Functional scheme. In particular, we show that the
continuum limit is independent of the form of the fermion action used, i.e. the
Wilson action and the nonperturbatively improved clover action.

We determine nonperturbatively the anomalous dimensions of the second moment
of nonsinglet parton densities from a continuum extrapolation of results
computed in quenched lattice simulations at different lattice spacings. We use
a Schr\"odinger functional scheme for the definition of the renormalization
constant of the relevant twist2 operator. In the region of renormalized
couplings explored, we obtain a good description of our data in terms of a
threeloop expression for the anomalous dimensions. The calculation can be used
for exploring values of the coupling where a perturbative expansion of the
anomalous dimensions is not valid a priori. Moreover, our results provide the
nonperturbative renormalization constant that connects hadron matrix elements
on the lattice, renormalized at a low scale, with the experimental results,
renormalized at much higher energy scales.

We present results for the reference scale r_0 in SU(3) Lattice Gauge Theory
for beta = 6/g_0^2 in the range 5.7 <= beta <= 6.57. The high relative accuracy
of 0.30.6% in r_0/a was achieved through good statistics, the application of a
multihit procedure and a variational approach in the computation of Wilson
loops. A precise definition of the force used to extract r_0 has been used
throughout the calculation which guarantees that r_0/a is a smooth function of
the bare coupling and that subsequent continuum extrapolations are possible.
The results are applied to the continuum extrapolations of the energy gap Delta
in the static quark potential and the scale L_max/r_0 used in the calculation
of the running coupling constant.

We consider the lattice regularization of N=1 supersymmetric YangMills
theory with Wilson fermions. This formulation breaks supersymmetry at any
finite lattice spacing; we discuss how Ward identities can be used to define a
supersymmetric continuum limit, which coincides with the point where the gluino
becomes massless. As a first step towards the understanding of the zero
gluinomass limit, we present results on the quenched lowlying spectrum of
SU(2) N=1 SuperYangMills, at $\beta=2.6$ on a $V=16^3 \times 32$ lattice, in
the OZI approximation. Our results, in spite of the quenched and OZI
approximations, are in remarkable agreement with theoretical predictions in the
supersymmetric theory, for the states with masses which are not expected to get
a large contribution from fermion loops.

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.

We discuss nonperturbative improvement of the vector current, using the
Schroedinger Functional formalism. By considering a suitable Ward identity, we
compute the improvement coefficient which gives the O(a) mixing of the tensor
current with the vector one.