  We numerically simulate a non-Abelian lattice gauge theory in two spatial
dimensions, with tensor networks (TN), up to intermediate sizes (>30 matter
sites) well beyond exact diagonalization. We focus on the SU(2) Yang-Mills
model in Hamiltonian formulation, with dynamical matter and minimally truncated
gauge field (hardcore gluon). Thanks to the TN sign-problem-free approach, we
characterize the phase diagram of the model at zero and finite baryon number as
a function of the quark bare mass and color charge. At intermediate system
sizes, we detect a liquid phase of quark-pair bound-state quasiparticles
(baryons), whose mass is finite towards the continuum limit. Interesting
phenomena arise at the transition boundary where color-electric and
color-magnetic terms are maximally frustrated: For low quark masses, we see
traces of potential deconfinement, while for high masses, signatures of a
possible topological order.
