Abstract We use the 4 − ϵ expansion to search for fixed points corresponding to 2 + 1 dimensional N $$ \mathcal{N} $$ =1 Wess-Zumino models of N Φ scalar superfields interacting through a cubic superpotential. In the N Φ = 3 case we classify all SUSY fixed points that are perturbatively unitary. In the N Φ = 4 and N Φ = 5 cases, we focus on fixed points where the scalar superfields form a single irreducible representation of the symmetry group (irreducible fixed points). For N Φ = 4 we show that the S5 invariant super Potts model is the only irreducible fixed point where the four scalar superfields are fully interacting. For N Φ = 5, we go through all Lie subgroups of O(5) and use the GAP system for computational discrete algebra to study finite subgroups of O(5) up to order 800. This analysis gives us three fully interacting irreducible fixed points. Of particular interest is a subgroup of O(5) that exhibits O(3)/Z2 symmetry. It turns out this fixed point can be generalized to a new family of models, with N Φ = N N − 1 2 $$ \frac{\mathrm{N}\left(\mathrm{N}-1\right)}{2} $$ − 1 and O(N)/Z2 symmetry, that exists for arbitrary integer N≥3.