Statistical models corresponding to a new class of braid matrices (ôN;N⩾3) presented in a previous paper are studied. Indices labeling states spanning the Nr dimensional base space of T(r)(θ), the rth order transfer matrix are so chosen that the operators W (the sum of the state labels) and (CP) (the circular permutation of state labels) commute with T(r)(θ). This drastically simplifies the construction of eigenstates, reducing it to solutions of relatively small number of simultaneous linear equations. Roots of unity play a crucial role. Thus for diagonalizing the 81 dimensional space for N=3, r=4, one has to solve a maximal set of five linear equations. A supplementary symmetry relates invariant subspaces pairwise [W=(r,Nr) and so on] so that only one of each pair needs study. The case N=3 is studied fully for r=(1,2,3,4). Basic aspects for all (N,r) are discussed. Full exploitation of such symmetries lead to a formalism quite different from, possibly generalized, algebraic Bethe ansatz. Chain Hamiltonians are studied. The specific types of spin flips they induce and propagate are pointed out. The inverse Cayley transform of the YB matrix giving the potential leading to factorizable S matrix is constructed explicitly for N=3 as also the full set of R̂tt relations. Perspectives are discussed in a final section.