where y ij is a Morse variable, given by y ij =1-exp(-r ij /a). The summation is over all powers of the y ij , the total polynomial order, subject to the constraint that the maximum order equals M. (For simplicity we denote the sum over powers by the single index m.)