A previously derived three-dimensional effective lattice theory describing the thermodynamics of QCD with heavy quarks in the cold and dense region is extended through order ∼ u 5 κ 8 in the combined character and hopping expansion of the original four-dimensional Wilson action. The systematics of the effective theory is investigated to determine its range of validity in parameter space. We demonstrate the severe cut-off effects due to lattice saturation, which afflict any lattice results at finite baryon density independent of the sign problem or the quality of effective theories, and which have to be removed by continuum extrapolation. We then show how the effective theory can be solved analytically by means of a linked cluster expansion, which is completely unaffected by the sign problem, in quantitative agreement with numerical simulations. As an application, we compute the cold nuclear equation of state of heavy QCD. Our continuum extrapolated result is consistent with a polytropic equation of state for non-relativistic fermions.

In this paper we test the complex Langevin algorithm for numerical simulations of a random matrix model of QCD with a first order phase transition to a phase of finite baryon density. We observe that a naive implementation of the algorithm leads to phase quenched results, which were also derived analytically in this article. We test several fixes for the convergence issues of the algorithm, in particular the method of gauge cooling, the shifted representation, the deformation technique and reweighted complex Langevin, but only the latter method reproduces the correct analytical results in the region where the quark mass is inside the domain of the eigenvalues. In order to shed more light on the issues of the methods we also apply them to a similar random matrix model with a milder sign problem and no phase transition, and in that case gauge cooling solves the convergence problems as was shown before in the literature.

No abstract

New results on bottomonium at nonzero temperature are presented, using the FASTSUM Generation 2L ensembles. Preliminary results for spectral function reconstruction using Kernel Ridge Regression, a machine learning technique, are shown as well and compared to results from the Maximum Entropy Method.

Abstract. We study a random matrix model for QCD at finite density via complex Langevin dynamics. This model has a phase transition to a phase with nonzero baryon density. We study the convergence of the algorithm as a function of the quark mass and the chemical potential and focus on two main observables: the baryon density and the chiral condensate. For simulations close to the chiral limit, the algorithm has wrong convergence properties when the quark mass is in the spectral domain of the Dirac operator. A possible solution of this problem is discussed.