Two-color QCD with heavy quarks at finite densities

We extend the density-of-states approach to gauge systems (LLR method [K. Langfeld, B. Lucini, and A. Rago, Phys. Rev. Lett. 109, 111601 (2012)]) to QCD at finite temperature and density with heavy quarks. The approach features an exponential error suppression and yields the Polyakov loop probability distribution function over a range of hundreds of orders of magnitude. $SU(2)$ gauge theory is considered in the confinement and the high-temperature phase, and at finite densities of heavy quarks. In the latter case, a smooth rise of the density with the chemical potential is observed, and no critical phenomenon associated with deconfinement due to finite densities is found.