BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Leonardo Giusti

Leonardo Giusti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-thhep-phnucl-thhep-exmath-phmath.MP

Catalog footprint

What is connected

32works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

Four-dimensional factorization of the fermion determinant in lattice QCD

In the last few years it has been proposed a one-dimensional factorization of the fermion determinant in lattice QCD with Wilson-type fermions that leads to a block-local action of the auxiliary bosonic fields. Here we propose a four-dimensional generalization of this factorization. Possible applications are more efficient parallelizations of Monte Carlo algorithms and codes, master field simulations, and multi-level integration.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-perturbative thermal QCD at all temperatures: the case of mesonic screening masses

We present a strategy based on the step-scaling technique to study non-perturbatively thermal QCD up to very high temperatures. As a first concrete application, we compute the flavour non-singlet meson screening masses at 12 temperatures covering the range from $T \sim 1$ GeV up to $\sim 160$ GeV in the theory with three massless quarks. The calculation is carried out by Monte Carlo simulations on the lattice by considering large spatial extensions in order to have negligible finite volume effects. For each temperature we have simulated 3 or 4 values of the lattice spacing, so as to perform the continuum limit extrapolation with confidence at a few permille accuracy. Chiral symmetry restoration manifests itself in our results through the degeneracy of the vector and the axial vector channels and of the scalar and the pseudoscalar ones. In the entire range of temperatures explored, the meson screening masses deviate from the free theory result, $2 πT$, by at most a few percent. These deviations, however, cannot be explained by the known leading term in the QCD coupling constant $g$ up to the highest temperature, where other contributions are still very relevant. In particular the vector-pseudoscalar mass splitting turns out to be of $O(g^4)$ in the entire range explored, and it remains clearly visible up to the highest temperature, where the two screening masses are still significantly different within our numerical precision. The pattern of different contributions that we have found explains why it has been difficult in the past to match non-perturbative lattice results at $T \sim 1$ GeV with the analytic behaviour at asymptotically high temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A novel computational paradigm for a precise determination of the hadronic contribution to $(g_μ-2)$ from lattice QCD

The hadronic contribution to the muon anomalous magnetic moment $a_μ=(g_μ-2)/2$ has to be determined at the per-mille level for the Standard Model prediction to match the expected final uncertainty of the ongoing E989 experiment. That is 3 times better than the current precision from the dispersive approach, and 5-15 times smaller than the uncertainty based on the purely theoretical determinations from lattice QCD. So far the stumbling-block is the large statistical error in the Monte Carlo evaluation of the required correlation functions which can hardly be tamed by brute force. In this talk we present our proposal to solve this problem by multi-level Monte Carlo integration, a technique which reduces the variance of correlators exponentially in the distance of the fields. We report the results of our feasibility tests for the computation of the Hadronic Vacuum Polarization on a lattice with a linear extension of 3~fm, a spacing of 0.065 fm, and a pion mass of 270 MeV. Indeed the two-level integration makes the contribution to the statistical error from long-distances de-facto negligible by accelerating its inverse scaling with the cost of the simulation. These findings establish multi-level Monte Carlo as a solid and efficient method for a precise lattice determination of the hadronic contribution to $a_μ$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Frequency-splitting estimators for single-propagator traces

In these proceedings we address the computation of quark-line disconnected diagrams in lattice QCD. The evaluation of these diagrams is required for many phenomenologically interesting observables, but suffers from large statistical errors due to the vacuum and random-noise contributions to their variances. Motivated by a theoretical analysis of the variances, we introduce a new family of stochastic estimators of single-propagator traces built upon a frequency splitting combined with a hopping expansion of the quark propagator, and test their efficiency in two-flavour QCD with pions as light as 190 MeV. The use of these estimators reduces the cost of the computation by one to two orders of magnitude over standard estimators depending on the fermion bilinear. As a concrete application, we show the impact of these findings on the computation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-perturbative definition of the QCD energy-momentum tensor on the lattice

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 Poincare' 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Domain decomposition, multi-level integration and exponential noise reduction in lattice QCD

We explore the possibility of computing fermionic correlators on the lattice by combining a domain decomposition with a multi-level integration scheme. The quark propagator is expanded in series of terms with a well defined hierarchical structure. The higher the order of a term, the (exponentially) smaller its magnitude, the less local is its dependence on the gauge field. Once inserted in a Wick contraction, the gauge-field dependence of the terms in the resulting series can be factorized so that it is suitable for multi-level Monte Carlo integration. We test the strategy in quenched QCD by computing the disconnected correlator of two flavor-diagonal pseudoscalar densities, and a nucleon two-point function. In either cases we observe a significant exponential increase of the signal-to-noise ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The large $N$ limit of the topological susceptibility of Yang-Mills gauge theory

We present a precise computation of the topological susceptibility $χ_{_\mathrm{YM}}$ of SU$(N)$ Yang-Mills theory in the large $N$ limit. The computation is done on the lattice, using high-statistics Monte Carlo simulations with $N=3, 4, 5, 6$ and three different lattice spacings. Two major improvements make it possible to go to finer lattice spacing and larger $N$ compared to previous works. First, the topological charge is implemented through the gradient flow definition; and second, open boundary conditions in the time direction are employed in order to avoid the freezing of the topological charge. The results allow us to extrapolate the dimensionless quantity $t_0^2χ_{_\mathrm{YM}}$ to the continuum and large $N$ limits with confidence. The accuracy of the final result represents a new quality in the verification of large $N$ scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The topological susceptibility in the large-N limit of SU(N) Yang-Mills theory

We compute the topological susceptibility of the SU(N) Yang-Mills theory in the large-N limit with a percent level accuracy. This is achieved by measuring the gradient-flow definition of the susceptibility at three values of the lattice spacing for N=3,4,5,6. Thanks to this coverage of parameter space, we can extrapolate the results to the large-N and continuum limits with confidence. Open boundary conditions are instrumental to make simulations feasible on the finer lattices at the larger N.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A novel computation of the thermodynamics of the SU(3) Yang-Mills theory

We present an accurate computation of the Equation of State of the SU(3) Yang-Mills theory using shifted boundary conditions in the temporal direction. In this framework, the entropy density s can be obtained in a simple way from the expectation value of the space-time components T0k of the energy-momentum tensor. At each given value of the temperature, s is measured in an independent way at several values of the lattice spacing. The extrapolation to the continuum limit shows small discretization effects with respect to the statistical errors of approximatively 0.5%.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Energy-momentum tensor on the lattice: non-perturbative renormalization in Yang--Mills theory

We construct an energy-momentum tensor on the lattice which satisfies the appropriate Ward Identities (WIs) and has the right trace anomaly in the continuum limit. It is defined by imposing suitable WIs associated to the Poincare` invariance of the continuum theory. These relations come forth when the length of the box in the temporal direction is finite, and they take a particularly simple form if the coordinate and the periodicity axes are not aligned. We implement the method for the SU(3) Yang--Mills theory discretized with the standard Wilson action in presence of shifted boundary conditions in the (short) temporal direction. By carrying out extensive numerical simulations, the renormalization constants of the traceless components of the tensor are determined with a precision of roughly half a percent for values of the bare coupling constant in the range 0<= g^2_0<=1.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-Gaussianities in the topological charge distribution of the SU(3) Yang--Mills theory

We study the topological charge distribution of the SU(3) Yang--Mills theory with high precision in order to be able to detect deviations from Gaussianity. The computation is carried out on the lattice with high statistics Monte Carlo simulations by implementing a naive discretization of the topological charge evolved with the Yang--Mills gradient flow. This definition is far less demanding than the one suggested from Neuberger's fermions and, as shown in this paper, in the continuum limit its cumulants coincide with those of the universal definition appearing in the chiral Ward identities. Thanks to the range of lattice volumes and spacings considered, we can extrapolate the results for the second and fourth cumulant of the topological charge distribution to the continuum limit with confidence by keeping finite volume effects negligible with respect to the statistical errors. Our best results for the topological susceptibility is t_0^2*chi=6.67(7)*10^-4, where t_0 is a standard reference scale, while for the ratio of the forth cumulant over the second we obtain R=0.233(45). The latter is compatible with the expectations from the large Nc expansion, while it rules out the theta-behavior of the vacuum energy predicted by the dilute instanton model. Its large distance from 1 implies that, in the ensemble of gauge configurations that dominate the path integral, the fluctuations of the topological charge are of quantum non-perturbative nature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Renormalization of the energy-momentum tensor on the lattice

We present the calculation of the non-perturbative renormalization constants of the energy-momentum tensor in the SU(3) Yang-Mills theory. That computation is carried out in the framework of shifted boundary conditions, where a thermal quantum field theory is formulated in a moving reference frame. The non-perturbative renormalization factors are then used to measure the Equation of State of the SU(3) Yang-Mills theory. Preliminary numerical results are presented and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Chiral symmetry breaking in QCD Lite

A distinctive feature of the presence of spontaneous chiral symmetry breaking in QCD is the condensation of low modes of the Dirac operator near the origin. The rate of condensation must be equal to the slope of (Mpi^2 Fpi^2)/2 with respect to the quark mass m in the chiral limit, where Mpi and Fpi are the mass and the decay constant of the Nambu-Goldstone bosons. We compute the spectral density of the (Hermitian) Dirac operator, the quark mass, the pseudoscalar meson mass and decay constant by numerical simulations of lattice QCD with two light degenerate Wilson quarks. We use CLS lattices at three values of the lattice spacing in the range 0.05-0.08 fm, and for several quark masses corresponding to pseudoscalar mesons masses down to 190 MeV. Thanks to this coverage of parameters space, we can extrapolate all quantities to the chiral and continuum limits with confidence. The results show that the low quark modes do condense in the continuum as expected by the Banks-Casher mechanism, and the rate of condensation agrees with the Gell-Mann-Oakes-Renner (GMOR) relation. For the renormalisation-group-invariant ratios we obtain [Σ^RGI]^(1/3)/F =2.77(2)(4) and Lambda^MSbar/F = 3.6(2), which correspond to [Σ^\MSbar(2 GeV)]^(1/3) =263(3)(4) MeV and F=85.8(7)(20) MeV if FK is used to set the scale by supplementing the theory with a quenched strange quark.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Equation of state of a relativistic theory from a moving frame

We propose a new strategy for determining the equation of state of a relativistic thermal quantum field theory by considering it in a moving reference system. In this frame an observer can measure the entropy density of the system directly from its average total momentum. In the Euclidean path integral formalism, this amounts to compute the expectation value of the off-diagonal components T_{0k} of the energy-momentum tensor in presence of shifted boundary conditions. The entropy is thus easily measured from the expectation value of a local observable computed at the target temperature T only. At large T, the temperature itself is the only scale which drives the systematic errors, and the lattice spacing can be tuned to perform a reliable continuum limit extrapolation while keeping finite-size effects under control. We test this strategy for the four-dimensional SU(3) Yang-Mills theory. We present precise results for the entropy density and its step-scaling function in the temperature range 0.9 T_c - 20 T_c. At each temperature, we consider four lattice spacings in order to extrapolate the results to the continuum limit. As a byproduct we also determine the ultraviolet finite renormalization constant of T_{0k} by imposing suitable Ward identities. These findings establish this strategy as a solid, simple and efficient method for an accurate determination of the equation of state of a relativistic thermal field theory over several orders of magnitude in T.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-perturbative renormalization of the energy-momentum tensor in SU(3) Yang-Mills theory

We present a strategy for a non-perturbative determination of the finite renormalization constants of the energy-momentum tensor in the SU(3) Yang-Mills theory. The computation is performed by imposing on the lattice suitable Ward Identites at finite temperature in presence of shifted boundary conditions. We show accurate preliminary numerical data for values of the bare coupling g_0^2 ranging for 0 to 1.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectral density of the Dirac operator in two-flavour QCD

We compute the spectral density of the (Hermitean) Dirac operator in Quantum Chromodynamics with two light degenerate quarks near the origin. We use CLS/ALPHA lattices generated with two flavours of O(a)-improved Wilson fermions corresponding to pseudoscalar meson masses down to 190 MeV, and with spacings in the range 0.05-0.08 fm. Thanks to the coverage of parameter space, we can extrapolate our data to the chiral and continuum limits with confidence. The results show that the spectral density at the origin is non-zero because the low modes of the Dirac operator do condense as expected in the Banks-Casher mechanism. Within errors, the spectral density turns out to be a constant function up to eigenvalues of approximately 80 MeV. Its value agrees with the one extracted from the Gell-Mann-Oakes-Renner relation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Testing the Witten-Veneziano mechanism with the Yang-Mills gradient flow on the lattice

We present a precise computation of the topological charge distribution in the $SU(3)$ Yang-Mills theory. It is carried out on the lattice with high statistics Monte Carlo simulations by employing the clover discretization of the field strength tensor combined with the Yang-Mills gradient flow. The flow equations are integrated numerically by a fourth-order structure-preserving Runge-Kutta method. We have performed simulations at four lattice spacings and several lattice sizes to remove with confidence the systematic errors in the second (topological susceptibility $χ_t^\text{YM}$) and the fourth cumulant of the distribution. In the continuum we obtain the preliminary results $t_0^2χ_t^\text{YM}=6.53(8)\times 10^{-4}$ and the ratio between the fourth and the second cumulant $R=0.233(45)$. Our results disfavour the $θ$-behaviour of the vacuum energy predicted by dilute instanton models, while they are compatible with the expectation from the large-$N_c$ expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Chiral condensate from the Banks-Casher relation

We report on our ongoing project of determining the chiral condensate of two-flavor QCD from the Banks-Casher relation. We compute the mode number of the O(a)-improved Wilson-Dirac operator for several values of Λ, and we discuss different fitting strategies to extract the chiral condensate from its mass and Λdependence. Our preliminary results haven been obtained at two different lattice spacings by using CLS-configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Measuring the entropy from shifted boundary conditions

We explore a new computational strategy for determining the equation of state of the SU(3) Yang-Mills theory. By imposing shifted boundary conditions, the entropy density is computed from the vacuum expectation value of the off-diagonal components T_{0k} of the energy-momentum tensor. A step-scaling function is introduced to span a wide range in temperature values. We present preliminary numerical results for the entropy density and its step-scaling function obtained at eight temperature values in the range T_c - 15 T_c. At each temperature, discretization effects are removed by simulating the theory at several lattice spacings and by extrapolating the results to the continuum limit. Finite-size effects are always kept below the statistical errors. The absence of ultraviolet power divergences and the remarkably small discretization effects allow for a precise determination of the step-scaling function in the explored temperature range. These findings establish this strategy as a viable solution for an accurate determination of the equation of state in a wide range of temperature values.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermal field theories and shifted boundary conditions

The analytic continuation to an imaginary velocity 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. The Poincare' invariance underlying a relativistic theory implies a dependence of the free-energy on the compact length L_0 and the shift xi only through the combination beta=L_0(1+xi^2)^(1/2). This in turn implies that the energy and the momentum distributions of the thermal theory are related, a fact which is encoded in a set of Ward identities among the correlators of the energy-momentum tensor. The latter have interesting applications in lattice field theory: they offer novel ways to compute thermodynamic potentials, and a set of identities to renormalize non-perturbatively the energy-momentum tensor. At fixed bare parameters the shifted boundary conditions also provide a simple method to vary the temperature in much smaller steps than with the standard procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Implications of Poincare symmetry for thermal field theories in finite-volume

The analytic continuation to an imaginary velocity $iξ$ 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ξ.P)]$, the Poincare invariance underlying a relativistic theory implies a dependence of the free-energy on $L_0$ and the shift $ξ$ only through the combination $β= L_0 \sqrt{1+ξ^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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Glueball masses from ratios of path integrals

By generalizing our previous work on the parity symmetry, the partition function of a Yang-Mills theory is decomposed into a sum of path integrals each giving the contribution from multiplets of states with fixed quantum numbers associated to parity, charge conjugation, translations, rotations and central conjugations. Ratios of path integrals and correlation functions can then be computed with a multi-level Monte Carlo integration scheme whose numerical cost, at a fixed statistical precision and at asymptotically large times, increases power-like with the time extent of the lattice. The strategy is implemented for the SU(3) Yang-Mills theory, and a full-fledged computation of the mass and multiplicity of the lightest glueball with vacuum quantum numbers is carried out at two values of the lattice spacing (0.17 and 0.12 fm).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermal momentum distribution from path integrals with shifted boundary conditions

For a thermal field theory formulated in the grand canonical ensemble, the distribution of the total momentum is an observable characterizing the thermal state. We show that its cumulants are related to thermodynamic potentials. In a relativistic system for instance, the thermal variance of the total momentum is a direct measure of the enthalpy. We relate the generating function of the cumulants to the ratio of (a) a partition function expressed as a Matsubara path integral with shifted boundary conditions in the compact direction, and (b) the ordinary partition function. In this form the generating function 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermal momentum distribution from shifted boundary conditions

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermodynamic potentials from shifted boundary conditions: the scalar-field theory case

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 the Ward identities associated with the invariance of the theory under the Poincare' 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A novel approach for computing glueball masses and matrix elements in Yang-Mills theories on the lattice

We make use of the global symmetries of the Yang-Mills theory on the lattice to design a new computational strategy for extracting glueball masses and matrix elements which achieves an exponential reduction of the statistical error with respect to standard techniques. By generalizing our previous work on the parity symmetry, the partition function of the theory is decomposed into a sum of path integrals each giving the contribution from multiplets of states with fixed quantum numbers associated to parity, charge conjugation, translations, rotations and central conjugations Z_N^3. Ratios of path integrals and correlation functions can then be computed with a multi-level Monte Carlo integration scheme whose numerical cost, at a fixed statistical precision and at asymptotically large times, increases power-like with the time extent of the lattice. The strategy is implemented for the SU(3) Yang--Mills theory, and a full-fledged computation of the mass and multiplicity of the lightest glueball with vacuum quantum numbers is carried out at a lattice spacing of 0.17 fm.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Glueball masses with exponentially improved statistical precision

We briefly review the computational strategy we have recently introduced for computing glueball masses and matrix elements, which achieves an exponential reduction of statistical errors compared to standard techniques. The global symmetries of the theory play a crucial role in the approach. We show how our previous work on parity can be generalized to other symmetries. In particular we discuss how to extract the mass of the 0++, 2++ and 0-+ lightest glueballs avoiding the exponential degradation of the signal to noise ratio. We present new numerical results and update the published ones.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Towards a precise determination of the topological susceptibility in the SU(3) Yang-Mills theory

An ongoing effort to compute the topological susceptibility for the SU(3) Yang-Mills theory in the continuum limit with a precison of about 2% is reported. The susceptibility is computed by using the definition of the charge suggested by Neuberger fermions for two values of the negative mass parameter s. Finite volume and discretization effects are estimated to meet this level of precision. The large statistics required has been obtained by using PCs of the INFN-GRID. Simulations with larger lattice volumes are necessary in order to better understanding the continuum limit at small lattice spacing values.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Symmetries and exponential error reduction in YM theories on the lattice: theoretical aspects and simulation results

The path integral of a quantum system with an exact symmetry can be written as a sum of functional integrals each giving the contribution from quantum states with definite symmetry properties. We propose a strategy to compute each of them, normalized to the one with vacuum quantum numbers, by a Monte Carlo procedure whose cost increases power-like with the time extent of the lattice. This is achieved thanks to a multi-level integration scheme, inspired by the transfer matrix formalism, which exploits the symmetry and the locality in time of the underlying statistical system. As a result the cost of computing the lowest energy level in a given channel, its multiplicity and its matrix elements is exponentially reduced with respect to the standard path-integral Monte Carlo. We briefly illustrate the approach in the simple case of the one-dimensional harmonic oscillator and discuss in some detail its extension to the four-dimensional Yang Mills theories. We report on our recent new results in the SU(3) Yang--Mills theory on the relative contribution to the partition function of the parity-odd states.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Spontaneous chiral symmetry breaking in QCD:a finite-size scaling study on the lattice

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.

0 citations0 reviewsNo code linkedOpen access