  We consider a $\mathbb{Z}_3$ gauge theory coupled to three degenerate massive
flavors of fermions, which we term "QZD". The spectrum can be computed in $1+1$
dimensions using tensor networks. In weak coupling the spectrum is that of the
expected mesons and baryons, although the corrections in weak coupling are
nontrivial, analogous to those of non-relativistic QED in 1+1 dimensions. In
strong coupling, besides the usual baryon, the singlet meson is a baryon
anti-baryon state. For two special values of the coupling constant, the
lightest baryon is degenerate with the lightest octet meson, and the lightest
singlet meson, respectively.
