Four-loop lattice-regularized vacuum energy density of the three-dimensional SU(3) + adjoint Higgs theory

The pressure of QCD admits at high temperatures a factorization into purely perturbative contributions from "hard" thermal momenta, and slowly convergent as well as non-perturbative contributions from "soft" thermal momenta. The latter can be related to various effective gluon condensates in a dimensionally reduced effective field theory, and measured there through lattice simulations. Practical measurements of one of the relevant condensates have suffered, however, from difficulties in extrapolating convincingly to the continuum limit. In order to gain insight on this problem, we employ Numerical Stochastic Perturbation Theory to estimate the problematic condensate up to 4-loop order in lattice perturbation theory. Our results seem to confirm the presence of "large" discretization effects, going like $a\ln(1/a)$, where $a$ is the lattice spacing. For definite conclusions, however, it would be helpful to repeat the corresponding part of our study with standard lattice perturbation theory techniques.Comment: 35 pages. v2: minor corrections, published versio

Renormalization of infrared contributions to the QCD pressure

Thanks to dimensional reduction, the infrared contributions to the QCD pressure can be obtained from two different three-dimensional effective field theories, called the Electrostatic QCD (Yang-Mills plus adjoint Higgs) and the Magnetostatic QCD (pure Yang-Mills theory). Lattice measurements have been carried out within these theories, but a proper interpretation of the results requires renormalization, and in some cases also improvement, i.e. the removal of terms of O(a) or O(a^2). We discuss how these computations can be implemented and carried out up to 4-loop level with the help of Numerical Stochastic Perturbation Theory.Comment: 7 pages, 4 figures, talk presented at Lattice 2006 (High temperature and density

Four-loop pressure of massless O(N) scalar field theory

Inspired by the corresponding problem in QCD, we determine the pressure of massless O(N) scalar field theory up to order g^6 in the weak-coupling expansion, where g^2 denotes the quartic coupling constant. This necessitates the computation of all 4-loop vacuum graphs at a finite temperature: by making use of methods developed by Arnold and Zhai at 3-loop level, we demonstrate that this task is manageable at least if one restricts to computing the logarithmic terms analytically, while handling the ``constant'' 4-loop contributions numerically. We also inspect the numerical convergence of the weak-coupling expansion after the inclusion of the new terms. Finally, we point out that while the present computation introduces strategies that should be helpful for the full 4-loop computation on the QCD-side, it also highlights the need to develop novel computational techniques, in order to be able to complete this formidable task in a systematic fashion.Comment: 34 page

The leading non-perturbative coefficient in the weak-coupling expansion of hot QCD pressure

Using Numerical Stochastic Perturbation Theory within three-dimensional pure SU(3) gauge theory, we estimate the last unknown renormalization constant that is needed for converting the vacuum energy density of this model from lattice regularization to the MSbar scheme. Making use of a previous non-perturbative lattice measurement of the plaquette expectation value in three dimensions, this allows us to approximate the first non-perturbative coefficient that appears in the weak-coupling expansion of hot QCD pressure.Comment: 16 pages. v2: published versio

The lattice ghost propagator in Landau gauge up to three loops using Numerical Stochastic Perturbation Theory

We complete our high-accuracy studies of the lattice ghost propagator in Landau gauge in Numerical Stochastic Perturbation Theory up to three loops. We present a systematic strategy which allows to extract with sufficient precision the non-logarithmic parts of logarithmically divergent quantities as a function of the propagator momentum squared in the infinite-volume and $a\to 0$ limits. We find accurate coincidence with the one-loop result for the ghost self-energy known from standard Lattice Perturbation Theory and improve our previous estimate for the two-loop constant contribution to the ghost self-energy in Landau gauge. Our results for the perturbative ghost propagator are compared with Monte Carlo measurements of the ghost propagator performed by the Berlin Humboldt university group which has used the exponential relation between potentials and gauge links.Comment: 8 pages, 6 figures, XXVII International Symposium on Lattice Field Theory - LAT2009, Beijin

Two-point functions of quenched lattice QCD in Numerical Stochastic Perturbation Theory

We summarize the higher-loop perturbative computation of the ghost and gluon propagators in SU(3) Lattice Gauge Theory. Our final aim is to compare with results from lattice simulations in order to expose the genuinely non-perturbative content of the latter. By means of Numerical Stochastic Perturbation Theory we compute the ghost and gluon propagators in Landau gauge up to three and four loops. We present results in the infinite volume and $a \to 0$ limits, based on a general fitting strategy.Comment: 3 pages, 5 figures, talk at conference QCHS-IX, Madrid 201

Two-point functions of quenched lattice QCD in Numerical Stochastic Perturbation Theory. (I) The ghost propagator in Landau gauge

This is the first of a series of two papers on the perturbative computation of the ghost and gluon propagators in SU(3) Lattice Gauge Theory. Our final aim is to eventually compare with results from lattice simulations in order to enlight the genuinely non-perturbative content of the latter. By means of Numerical Stochastic Perturbation Theory we compute the ghost propagator in Landau gauge up to three loops. We present results in the infinite volume and $a \to 0$ limits, based on a general strategy that we discuss in detail.Comment: 27 pages, 11 figure

An efficient method to compute the residual phase on a Lefschetz thimble

We propose an efficient method to compute the so-called residual phase that appears when performing Monte Carlo calculations on a Lefschetz thimble. The method is stochastic and its cost scales linearly with the physical volume, linearly with the number of stochastic estimators and quadratically with the length of the extra dimension along the gradient flow. This is a drastic improvement over previous estimates of the cost of computing the residual phase. We also report on basic tests of correctness and scaling of the code.Comment: New simulations, new plot, new appendix added. To appear in PRD. 9 pages, 3 figure