BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
P. J. Silva contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
In this proceeding, we explain a few steps for an alternative extraction of the spectral density of a two-point function (propagator) based on a discrete set of data points. We present a so-called Tikhonov regularization of this particular inverse problem. We test it on 2 cases: lattice 0++} glueball data and mock gluon data.
The propagator of a physical degree of freedom ought to obey a Källén-Lehmann spectral representation, with positive spectral density. The latter quantity is directly related to a cross section based on the optical theorem. The spectral density is a crucial ingredient of a quantum field theory with elementary and bound states, with a direct experimental connection as the masses of the excitations reflect themselves into (continuum) $δ$-singularities. In usual lattice simulational approaches to the QCD spectrum the spectral density itself is not accessed. The (bound state) masses are extracted from the asymptotic exponential decay of the two-point function. Given the importance of the spectral density, each nonperturbative continuum approach to QCD should be able to adequately describe it or to take into proper account. In this work, we wish to present a first trial in extracting an estimate for the scalar glueball spectral density in SU(3) gluodynamics using lattice gauge theory.
In the past years a good comprehension of the infrared gluon propagator has been achieved, with a good qualitative agreement between lattice results and Dyson-Schwinger equations. However, lattice simulations have been performed at physical volumes which are close to 20 fm but using a large lattice spacing. The interplay between volume effects and lattice spacing effects has not been investigated. Here we aim to fill this gap and address how the two effects change the gluon propagator in the infrared region. Furthermore, we provide infinite volume extrapolations which take into account the finite volume and finite lattice spacing. We also report on preliminary results for the gluon propagator at finite temperature.
We show how to improve the molecular dynamics step of Hybrid Monte Carlo, both by tuning the integrator using Poisson brackets measurements and by the use of force gradient integrators. We present results for moderate lattice sizes.
The infrared behaviour of the lattice Landau gauge gluon propagator is discussed, combining results from simulations with different volumes and lattice spacings. In particular, the Cucchieri-Mendes bounds are computed and their implications for D(0) discussed.
The gluon and ghost propagators are computed using both quenched and dynamical configurations for the same lattice spacings. The Wilson fermions are simulated at several quark masses. Furthermore, the effect of the Gribov copies is evaluated for all sets of configurations.
The connection between QCD, a nonlocal Nambu--Jona-Lasinio type model and the Landau gauge gluon propagator is explored. This two point function is parameterized by a functional form which is compatible with Dyson-Schwinger and lattice QCD results. Demanding the nonlocal model to reproduce the experimental values for the pion mass, the pion decay constant, $Γ_{π\rightarrow γγ}$ and the light quark condensate we conclude that low energy physics does not distinguish between the so-called decoupling and scaling solutions of the Dyson-Schwinger equations. This result means that, provided that the model parameters are chosen appropriately, one is free to chose any of the above scenarios. Furthermore, the nonlocal Nambu--Jona-Lasinio quark model considered here is chiral invariant and satisfies the GMOR relation at the 1% level of precision.
We present initial results of the use of Force Gradient integrators for lattice field theories. These promise to give significant performance improvements, especially for light fermions and large lattices. Our results show that this is indeed the case, indicating a speed-up of more than a factor of two, which is expected to increase as the integration step size becomes smaller for larger lattices and smaller fermion masses.
The Landau gauge lattice gluon propagator is discussed for different sets of lattices. Particular attention is given to its infrared properties. Our results show that the lattice propagator can be made compatible with either the decoupling-like or the scaling-like solution of the Dyson-Schwinger equations. Furthermore, the analysis of the Cucchieri-Mendes bounds is performed onsidering large volume simulations and the Oliveira-Silva ratios are computed. If the first do not give a clear answer about the value of D(0), the second method favors a $D(0) = 0$. Finally, the SU(3) and SU(2) propagators are compared in the infrared. It comes out that the propagators are different although the infrared exponents seem to be similar. The analysis suggests a scaling behaviour $D(0) \sim N$ with the gauge group SU(N).
We report on the infrared behaviour of the lattice gluon propagator. The bounds on D(0) and their behaviour with the lattice volume are investigated, together with the full propagator. Moreover, a study of the gluon propagator using different lattice spacings is carried out.
We discuss how the integrators used for the Hybrid Monte Carlo (HMC) algorithm not only approximately conserve some Hamiltonian $H$ but exactly conserve a nearby shadow Hamiltonian (\tilde H), and how the difference $ΔH \equiv \tilde H - H $ may be expressed as an expansion in Poisson brackets. By measuring average values of these Poisson brackets over the equilibrium distribution $\propto e^{-H}$ generated by HMC we can find the optimal integrator parameters from a single simulation. We show that a good way of doing this in practice is to minimize the variance of $ΔH$ rather than its magnitude, as has been previously suggested. Some details of how to compute Poisson brackets for gauge and fermion fields, and for nested and force gradient integrators are also presented.