Search

We present a strategy to define non-perturbatively the energy-momentum tensor in Quantum Chromodynamics (QCD) which satisfies the appropriate Ward identities and has the right trace anomaly. The tensor is defined by regularizing the theory on a lattice, and by fixing its renormalization constants non-perturbatively by suitable Ward identities associated to the Poincaré invariance of the continuum theory. The latter are derived in thermal QCD with a non-zero imaginary chemical potential formulated in a moving reference frame. A renormalization group analysis leads to simple renormalization- group-invariant definitions of the gluonic and fermionic contributions to either the singlet or the non-singlet components of the tensor, and therefore of their form factors among physical states. The lattice discussion focuses on the Wilson discretization of quark fields but the strategy is general. Specific to that case, we also carry out the analysis for the on-shell O(a)-improvement of the energy-momentum tensor. The renormalization and improvement programs profit from the fact that, as shown here, the thermal theory enjoys de-facto automatic O(a)-improvement at finite temperature. The validity of the proposal is scrutinized analytically by a study to 1-loop order in lattice perturbation theory with shifted and twisted (for quarks only) boundary conditions. The latter provides also additional useful insight for a precise non-perturbative calculation of the renormalization constants. The strategy proposed here is accessible to Monte Carlo computations, and in this sense it provides a practical way to define non-perturbatively the energy-momentum tensor in QCD.

The analytic continuation to an imaginary velocity $i\xi$ of the canonical partition function of a thermal system expressed in a moving frame has a natural implementation in the Euclidean path-integral formulation in terms of shifted boundary conditions. Writing the Boltzmann factor as $exp[-L_0(H-i\xi.P)]$, the Poincare invariance underlying a relativistic theory implies a dependence of the free-energy on $L_0$ and the shift $\xi$ only through the combination $eta= L_0 sqrt{1+\xi^2}$. This in turn implies a set of Ward identities, some of which were previously derived by us, among the correlators of the energy-momentum tensor. In the infinite-volume limit they lead to relations among the cumulants of the total energy distribution and those of the momentum, i.e. they connect the energy and the momentum distributions in the canonical ensemble. In finite volume the Poincare symmetry translates into exact relations among partition functions and correlation functions defined with different sets of (generalized) periodic boundary conditions. They have interesting applications in lattice field theory. In particular, they offer Ward identities to renormalize non-perturbatively the energy-momentum tensor and novel ways to compute thermodynamic potentials. At fixed bare parameters they also provide a simple method to vary the temperature in much smaller steps than with the standard procedure

In a thermal field theory, the cumulants of the momentum distribution can be extracted from the dependence of the Euclidean path integral on a shift in the fields built into the temporal boundary condition. When combined with theWard identities associated with the invariance of the theory under the Poincaré group, thermodynamic potentials such as the entropy or the pressure can be directly inferred from the response of the system to the shift. Crucially the argument holds, up to harmless finite-size and discretization effects, even if translational and rotational invariance are broken to a discrete subgroup of finite shifts and rotations such as in a lattice box. The formulas are thus applicable at finite lattice spacing and volume provided the derivatives are replaced by their discrete counterpart, and no additive or multiplicative ultraviolet-divergent renormalizations are needed to take the continuum limit. In this paper we present a complete derivation of the relevant formulas in the scalar field theory, where several technical complications are avoided with respect to gauge theories. As a by-product we obtain a recursion relation among the cumulants of the momentum distribution, and formulae for finite-volume corrections to several well-known thermodynamic identities.

At finite temperature the distribution of the total momentum is an observable characterizing the thermal state of a field theory, and its cumulants are related to thermodynamic potentials. In a relativistic system at zero chemical potential, for instance, the thermal variance of the total momentum is a direct measure of the entropy. We relate the generating function of the cumulants to the ratio of a path integral with properly shifted boundary conditions in the compact direction over the ordinary partition function. In this form it is well suited for Monte-Carlo evaluation, and the cumulants can be extracted straightforwardly. We test the method in the SU(3) Yang–Mills theory, and obtain the entropy density at three different temperatures.

The Banks-Casher relation links the spontaneous breaking of chiral symmetry in QCD to the presence of a non-zero density of quark modes at the low end of the spectrum of the Dirac operator. Spectral observables like the number of modes in a given energy interval are renormalizable and can therefore be computed using the Wilson formulation of lattice QCD even though the latter violates chiral symmetry at energies on the order of the inverse lattice spacing. Using numerical simulations, we find (in two-flavour QCD) that the low quark modes do condense in the expected way. In particular, the chiral condensate can be accurately calculated simply by counting the low modes on large lattices. Other spectral observables can be considered as well and have a potentially wide range of uses. © 2009 SISSA.

We investigate two-point correlation functions of left-handed currents computed in quenched lattice QCD with the Neuberger-Dirac operator. We consider two lattice spacings a~0.09,0.12 fm and two different lattice extents L~ 1.5, 2.0 fm; quark masses span both the p- and the epsilon-regimes. We compare the results with the predictions of quenched chiral perturbation theory, with the purpose of testing to what extent the effective theory reproduces quenched QCD at low energy. In the p-regime we test volume and quark mass dependence of the pseudoscalar decay constant and mass; in the epsilon-regime, we investigate volume and topology dependence of the correlators. While the leading order behaviour predicted by the effective theory is very well reproduced by the lattice data in the range of parameters that we explored, our numerical data are not precise enough to test next-to-leading order effects.

Spontaneous chiral symmetry breaking in QCD with massless quarks at infinite volume can be seen in a finite box by studying, for instance, the dependence of the chiral condensate from the volume and the quark mass. We perform a feasibility study of this program by computing the quark condensate on the lattice in the quenched approximation of QCD at small quark masses. We carry out simulations in various topological sectors of the theory at several volumes, quark masses and lattice spacings by employing fermions with an exact chiral symmetry, and we focus on observables which are infrared stable and free from mass-dependent ultraviolet divergences. The numerical calculation is carried out with an exact variance-reduction technique, which is designed to be particularly efficient when spontaneous symmetry breaking is at work in generating a few very small low-lying eigenvalues of the Dirac operator. The finite-size scaling behaviour of the condensate in the topological sectors considered agrees, within our statistical accuracy, with the expectations of the chiral effective theory. Close to the chiral limit we observe a detailed agreement with the first Leutwyler-Smilga sum rule. By comparing the mass, the volume and the topology dependence of our results with the predictions of the chiral effective theory, we extract the corresponding low-energy constant.

In this second report on our recent numerical simulations of two-flavour QCD, we provide further technical details on the simulations and describe the methods we used to extract the meson masses and decay constants from the generated ensembles of gauge fields. Among the topics covered are the choice of the DD-HMC parameters, the issue of stability, autocorrelations and the statistical error analysis. Extensive data tables are included as well as a short discussion of the quark-mass dependence in partially quenched QCD, supplementing the physics analysis that was presented in the first paper in this series.

Recent conceptual, algorithmic and technical advances allow numerical simulations of lattice QCD with Wilson quarks to be performed at significantly smaller quark masses than was possible before. Here we report on simulations of two-flavour QCD at sea-quark masses from slightly above to approximately 1/4 of the strange-quark mass, on lattices with up to 64x32^3 points and spacings from 0.05 to 0.08 fm. Physical sea-quark effects are clearly seen on these lattices, while the lattice effects appear to be quite small, even without O(a) improvement. A striking result is that the dependence of the pion mass on the sea-quark mass is accurately described by leading-order chiral perturbation theory up to meson masses of about 500 MeV.

We study the spectral gap of the Wilson--Dirac operator in two-flavour lattice QCD as a function of the lattice spacing $a$, the space-time volume $V$ and the current-quark mass $m$. It turns out that the median of the probability distribution of the gap scales proportionally to $m$ and that its width is practically equal to $a/sqrt{V}$. In particular, numerical simulations are safe from accidental zero modes in the large-volume regime of QCD.

We present a strategy designed to separate several possible origins of the well-known enhancement of the Delta{I}=1/2 amplitude in non-leptonic kaon decays. In particular, we seek to disentangle the contribution of physics at the typical QCD scale (soft-gluon 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 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 low-lying eigenmodes of the Dirac operator, are reduced by an exact low-mode averaging. To demonstrate the feasibility of the technique, we compute the two-point correlator of the left-handed 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 epsilon-regime. By matching the results with the corresponding (quenched) chiral perturbation theory expressions, an estimate of (quenched) low-energy constants can be obtained. We find agreement between the quenched values of F extrapolated from the p-regime and extracted in the epsilon-regime.

In the epsilon-regime of QCD the main features of the spectrum of the low-lying eigenvalues of the (euclidean) Dirac operator are expected to be described by a certain universality class of random matrix models. In particular, the latter predict the joint statistical distribution of the individual eigenvalues in any topological sector of the theory. We compare some of these predictions with high-precision numerical data obtained from lattice QCD for a range of lattice spacings and volumes. While no complete matching is observed, the results agree with theoretical expectations at volumes larger than about 5 fm^4

By combining the first unquenched lattice computation of the B-meson binding energy and the two-loop contribution to the lattice HQET residual mass, we determine the (ar{{MS}}) (b)-quark mass, (ar{m}_{b}(ar{m}_{b})). The inclusion of the two-loop corrections is essential to extract (ar{m}_{b}(ar{m}_{b})) with a precision of ({cal O}(Lambda^{2}_{QCD}/m_{b})), which is the uncertainty due to the renormalon singularities in the perturbative series of the residual mass. Our best estimate is (ar{m}_{b}(ar{m}_{b}) = (4.26 pm 0.09) { m GeV}), where we have combined the different errors in quadrature. A detailed discussion of the systematic errors contributing to the final number is presented. Our results have been obtained on a sample of (60) lattices of size (24^{3} imes 40) at (eta =5.6), using the Wilson action for light quarks and the lattice HQET for the (b) quark, at two values of the sea quark masses. The quark propagators have been computed using the unquenched links generated by the T(chi)L Collaboration.

In this poster we present a few preliminary results obtained using our method to fix generic covariant gauges on the lattice. We have computed the gluon propagator and we have found a sensitive dependence on the gauge parameter

We perform a Next-to-Leading order analysis of Delta S=2 processes beyond the Standard Model. Combining the recently computed NLO anomalous dimensions and the B parameters of the most general Delta S=2 effective Hamiltonian, we give an analytic formula for Delta M_K and epsilon_K in terms of the Wilson coefficients at the high energy scale. This expression can be used for any extension of the Standard Model with new heavy particles. Using this result, we consider gluino-mediated contributions to Delta S=2 transitions in general SUSY models and provide an improved analysis of the constraints on off-diagonal mass terms between the first two generations of down-type squarks. Finally, we improve the constraints on R-violating couplings from Delta M_K and epsilon_K