In this equation, each possible occupation of each active orbital n i is given its own matrix or vector A n i . The maximum inner dimension of these matrices is called the bond order M and is the number of states retained during the renormalization step. As M → ∞, results obtained with this method approach those obtained with full diagonalization (although useful results