“…(5). For this purpose, 11 it is quite advisable to pass to the matrix representation of the cell HamiltonianĤ k in the basis spanned over four available 12 states of two Heisenberg spins σ z k,1 and σ z k,2 : 13 | ↑, ↑⟩ k = |↑⟩ k,1 |↑⟩ k,2 , | ↓, ↓⟩ k = |↓⟩ k,1 |↓⟩ k,2 , | ↑, ↓⟩ k = |↑⟩ k,1 |↓⟩ k,2 , | ↓, ↑⟩ k = |↓⟩ k,1 |↑⟩ k,2 , 14 whereas | ↑⟩ k,i and | ↓⟩ k,i denote two eigenvectors of the spin operatorσ z k,i with the respective eigenvalues σ z k,i = 1/2 and 15 −1/2. After a straightforward diagonalization of the cell HamiltonianĤ k one obtains the following four eigenvalues:…”