• 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 four-quark 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 non-perturbatively 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 twist-2 operator in pion states from non-perturbative lattice calculations. We find that the non-perturbatively 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 non-perturbative renormalization, finite size effects and effects of a non-vanishing 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 non-singlet, twist-2 operator corresponding to the average momentum of non-singlet 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 non-perturbatively 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 one-loop perturbative calculation of the discretisation effects has been carried out for both the Wilson and the Clover-improved actions and for a large number of lattice resolutions. The non-perturbative 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 multi-teraflop 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 b-quark and m_c^{RGI} = 1.681(36) GeV for the c-quark, 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 non-perturbatively the evolution of the twist-2 operators corresponding to the average momentum of non-singlet quark densities. The calculation is based on a finite-size 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 by-product we apply the non-perturbatively 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 twist-2 operator corresponding to the average momentum of non-singlet 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 non-perturbatively improved clover actions in order to control better the extrapolation to the continuum limit. Moreover, we compute, fully non-perturbatively, 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 heavy--light 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 twist-2 operator corresponding to the average momentum of non-singlet 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 non-perturbative 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 two--scale 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 non--perturbative 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 two--scale 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 four-fermion operators which appear in the effective weak Hamiltonian of the Standard Model. Our calculations are done using Wilson fermions and the parity-odd 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 T-Flops processing power. This paper describes the main features of this new development.
  • We present the results of a non-perturbative determination of the pion matrix element of the twist-2 operator corresponding to the average momentum of non-singlet 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 non-perturbatively 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 non-singlet parton distribution function. The matching of the results of our non-perturbative calculation with the ones from hadronic matrix elements allows us to obtain eventually a renormalization group invariant average momentum of non-singlet 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 non-singlet parton densities. The evidence is provided by a non-perturbative 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 non-perturbatively improved clover action.
  • We determine non-perturbatively the anomalous dimensions of the second moment of non-singlet 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 twist-2 operator. In the region of renormalized couplings explored, we obtain a good description of our data in terms of a three-loop 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 non-perturbative 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.3-0.6% in r_0/a was achieved through good statistics, the application of a multi-hit 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 Yang--Mills 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 gluino-mass limit, we present results on the quenched low-lying spectrum of SU(2) N=1 Super-Yang--Mills, 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 non-perturbatively. The calculation is based on an extension of a finite-size technique previously employed to compute the running coupling in quenched QCD. As a by-product we obtain the $\Lambda$--parameter in this theory with completely controlled errors.
  • We discuss non-perturbative 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.