A many-body generalization of quasi-solvable models with type C N-fold supersymmetry (I) regular cases

We make a generalization of the type C monomial space of a single variable, which was introduced in the construction of type C N-fold supersymmetry, to several variables. Then, we construct the most general quasi-solvable second-order operators preserving this multivariate type C space. These operators of several variables are characterized by the fact that two different polynomial type solutions are available. In particular, we investigate and classify all the possible Schrödinger operators realized as a subclass of this family. It turns out that the rational, hyperbolic, and trigonometric Calogero-Sutherland models as well as some particular type of the elliptic Inozemtsev system, all associated with the BCM root system, fall within the class.