No Arabic abstract
Lattice QCD with heavy quarks reduces to a three-dimensional effective theory of Polyakov loops, which is amenable to series expansion methods. We analyse the effective theory in the cold and dense regime for a general number of colours, $N_c$. In particular, we investigate the transition from a hadron gas to baryon condensation. For any finite lattice spacing, we find the transition to become stronger, i.e. ultimately first-order, as $N_c$ is made large. Moreover, in the baryon condensed regime, we find the pressure to scale as $psim N_c$ through three orders in the hopping expansion. Such a phase differs from a hadron gas with $psim N_c^0$, or a quark gluon plasma, $psim N_c^2$, and was termed quarkyonic in the literature, since it shows both baryon-like and quark-like aspects. A lattice filling with baryon number shows a rapid and smooth transition from condensing baryons to a crystal of saturated quark matter, due to the Pauli principle, and is consistent with this picture. For continuum physics, the continuum limit needs to be taken before the large $N_c$ limit, which is not yet possible in practice. However, in the controlled range of lattice spacings and $N_c$-values, our results are stable when the limits are approached in this order. We discuss possible implications for physical QCD.
After combined character and hopping expansions and integration over the spatial gauge links, lattice QCD reduces to a three-dimensional $SU(3)$ Polyakov loop model with complicated interactions. A simple truncation of the effective theory is valid for heavy quarks on reasonably fine lattices and can be solved by linked cluster expansion in its effective couplings. This was used ealier to demonstrate the onset transition to baryon matter in the cold and dense regime. Repeating these studies for general $N_c$, one finds that for large $N_c$ the onset transition becomes first-order, and the pressure scales as $psim N_c$ through three consecutive orders in the hoppoing expansion. These features are consistent with the formal definition of quarkyonic matter given in the literature. We discuss the implications for $N_c=3$ and physical QCD.
Combining strong coupling and hopping expansion one can derive a dimensionally reduced effective theory of lattice QCD. This theory has a reduced sign problem, is amenable to analytic evaluation and was successfully used to study the cold and dense regime of QCD for sufficiently heavy quarks. We show results from the evaluation of the effective theory for arbitrary $N_c$ up to $kappa^4$. The inclusion of gauge corrections is also investigated. We find that the onset transition to finite baryon number density steepens with growing $N_c$ even for $T eq 0$. This suggests that in the large $N_c$ limit the onset transition is first order up to the deconfinement transition. Beyond the onset, the pressure is shown to scale as $p sim N_c$ through three orders in the hopping expansion, which is characteristic for a phase termed quarkyonic matter in the literature.
During the last years it has become possible to address the cold and dense regime of QCD directly for sufficiently heavy quarks, where combined strong coupling and hopping expansions are convergent and a 3d effective theory can be derived, which allows to control the sign problem either in simulations or by fully analytic calculations. In this contribution we review the effective theory and study the $N_c$-dependence of the nuclear liquid gas transition, as well as the equation of state of baryonic matter in the strong coupling limit. We find the transition to become more strongly first order with growing $N_c$, suggesting that in the large $N_c$ limit its critical endpoint moves to high temperatures to connect with the deconfinement transition. Furthermore, to leading and next-to-leading order in the strong coupling and hopping expansions, respectively, the pressure is found to scale as $psim N_c$. This suggests that baryonic and quarkyonic matter might be the same at nuclear densities. Further work is needed to see whether this result is stable under gauge corrections.
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, kappa, whose action is correct to kappa^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 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 gap between the onset of isospin and baryon condensation.
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 $sim u^5kappa^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.