We describe multiquark clusters in quark matter within a Beth-Uhlenbeck approach in a background gluon field coupled to the underlying chiral quark dynamics using the Polyakov gauge. An effective potential for the traced Polyakov loop is used to establish the center symmetry of the SU(3) color which suppresses colored states and its dynamical breaking as an aspect of the confinement/deconfinement transition. Quark confinement is modeled by a large quark mass in vacuum which is motivated by a confining density functional approach. A multiquark cluster containing n quarks and antiquarks is described as a binary composite of smaller subclusters $$n_1$$
n
1
and $$n_2$$
n
2
($$n_1+n_2=n$$
n
1
+
n
2
=
n
). It has a spectrum consisting of a bound state and a scattering state continuum. For the corresponding cluster-cluster phase shifts we use simple ansätze that capture the Mott dissociation of clusters as a function of temperature and chemical potential. We go beyond the simple “step-up-step-down” model that ignores continuum correlations and introduce an improved model that includes them in a generic form. In order to explain the model, we restrict ourselves here to the cases where the cluster size is $$1 \le n \le 6$$
1
≤
n
≤
6
. A striking result is the suppression of the abundance of colored multiquark clusters at low temperatures by the coupling to the Polyakov loop and their importance for a quantitative description of lattice QCD thermodynamics at non-vanishing baryochemical potentials. An important ingredient are Polyakov-loop generalized distribution functions of n-quark clusters which are derived here for the first time. Within our approach we calculate thermodynamic properties such as baryon density and entropy. We demonstrate that the limits of a hadron resonance gas at low temperatures and $$\mathcal {O}(g^2)$$
O
(
g
2
)
perturbative QCD at high temperatures are correctly reproduced. A comparison with lattice calculations shows that our model is able to give a unified, systematic approach to describe properties of the quark-gluon-hadron system.