Researcher profile

Leonardo Giusti

Leonardo Giusti contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 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.

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.

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_μ$.

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.

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.

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.

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.

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.

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).

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.

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.

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.

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.

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.

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.