  We compute the spectra of flux tubes formed between a static quark antiquark
pair up to a significant number of excitations and for eight symmetries of the
flux tubes, up to $\Delta_u$, using pure $SU(3)$ gauge lattice QCD in $3+1$
dimensions. To accomplish this goal, we use a large set of appropriate
operators, an anisotropic tadpole improved action, smearing techniques, and
solve a generalized eigenvalue problem. Moreover, we compare our results with
the Nambu-Goto string model to evaluate possible tensions which could be a
signal for novel phenomena.
