Analytical results for large-$N_c$ scalar QCD$_2$

We study large-$N_c$ scalar QCD$_2$, a $1+1$-dimensional confining gauge theory with fundamental scalar quarks, whose meson spectrum is governed by a Bethe-Salpeter equation structurally parallel to the 't Hooft equation. Exploiting this structural analogy, we develop a nonperturbative analytic framework, based on integrability and inspired by the Fateev-Lukyanov-Zamolodchikov (FLZ) method, originally devised for the 't Hooft model and later extended in our previous works. Notably, the same Bethe-Salpeter equation also arises in the description of interchain mesons in the doubled Ising model coupled via a spin-spin interaction term. Within the FLZ approach, we find spectral sums and derive a systematic large-$n$ WKB expansion for the meson spectrum. The analytic results reproduce the expected behavior in key asymptotic regimes, such as the near-critical limit $m\to g/\sqrt{\pi}$ and the heavy-quark regime $m\gg g$, and are in good agreement with numerical data. Finally, by analytically continuing the mass parameter into the complex plane, we uncover two infinite families of singularities where individual mesons become massless, suggesting a hidden connection to nontrivial Conformal Field Theories.