Using combined strong coupling and hopping parameter expansions, we derive an effective three-dimensional theory from thermal lattice QCD with heavy Wilson quarks. The theory depends on traced Polyakov loops only and correctly reflects the centre symmetry of the pure gauge sector as well as its breaking by finite mass quarks. It is valid up to certain orders in the lattice gauge coupling and hopping parameter, which can be systematically improved. To its current order it is controlled for lattices up to N τ ∼ 6 at finite temperature. For nonzero quark chemical potentials, the effective theory has a fermionic sign problem which is mild enough to carry out simulations up to large chemical potentials. Moreover, by going to a flux representation of the partition function, the sign problem can be solved. As an application, we determine the deconfinement transition and its critical end point as a function of quark mass and all chemical potentials.

A three-dimensional effective lattice theory of Polyakov loops is derived from QCD by expansions in the fundamental character of the gauge action, u, and the hopping parameter, κ, whose action is correct to κ n u m with n + m = 4. At finite baryon density, the effective theory has a sign problem which meets all criteria to be simulated by complex Langevin as well as by Monte Carlo on small volumes. The theory is valid for the thermodynamics of heavy quarks, where its predictions agree with simulations of full QCD at zero and imaginary chemical potential. In its region of convergence, it is moreover amenable to perturbative calculations in the small effective couplings. In this work we study the challenging cold and dense regime. We find unambiguous evidence for the nuclear liquid gas transition once the baryon chemical potential approaches the baryon mass, and calculate the nuclear equation of state in the limit of heavy baryons. In particular, we find a negative binding energy per nucleon causing the condensation, whose absolute value decreases exponentially as mesons get heavier. For decreasing meson mass, we observe a first order liquid gas transition with an endpoint at some finite temperature, as well as a gap between the onset of isospin and baryon condensation.

We study lattice QCD with four flavors of staggered quarks. In the limit of infinite gauge coupling, "dual" variables can be introduced, which render the finite-density sign problem mild and allow a full determination of the μ-T phase diagram by Monte Carlo simulations, also in the chiral limit. However, the continuum limit coincides with the weak coupling limit. We propose a strong-coupling expansion approach towards the continuum limit. We show first results, including the phase diagram and its chiral critical point, from this expansion truncated at next-to-leading order.

The spectral function related to the correlator of two colour-electric fields along a Polyakov loop determines the momentum diffusion coefficient of a heavy quark near rest with respect to a heat bath. We compute this spectral function at next-to-leading order, O(alpha_s^2), in the weak-coupling expansion. The high-frequency part of our result (omega >> T), which is shown to be temperature-independent, is accurately determined thanks to asymptotic freedom; the low-frequency part of our result (omega << T), in which Hard Thermal Loop resummation is needed in order to cure infrared divergences, agrees with a previously determined expression. Our result may help to calibrate the overall normalization of a lattice-extracted spectral function in a perturbative frequency domain T << omega << 1/a, paving the way for a non-perturbative estimate of the momentum diffusion coefficient at omega -> 0. We also evaluate the colour-electric Euclidean correlator, which could be directly compared with lattice simulations. As an aside we determine the Euclidean correlator in the lattice strong-coupling expansion, showing that through a limiting procedure it can in principle be defined also in the confined phase of pure Yang-Mills theory, even if a practical measurement could be very noisy there.Comment: 38 page

Lattice QCD at finite density suffers from a severe sign problem, which has so far prohibited simulations of the cold and dense regime. Here we study the onset of nuclear matter employing a three-dimensional effective theory derived by combined strong coupling and hopping expansions, which is valid for heavy but dynamical quarks and has a mild sign problem only. Its numerical evaluations agree between a standard Metropolis and complex Langevin algorithm, where the latter is free of the sign problem. Our continuum extrapolated data approach a first order phase transition at μ(B) ≈ m(B) as the temperature approaches zero. An excellent description of the data is achieved by an analytic solution in the strong coupling limit.