Spinring:  N=6, spin=1/2

Heisenberg-Modell,
vollständige Diagonalisierung,
allgemeine Wechselwirkungsmatrix

(* <<DiscreteMath`KroneckerDelta` *)

Initialisierung

N, s, Dimension von H und Basis

Setze -J = 1;

NS = 6;
spin = (1/2);
\(\(DimH = \((2*spin + 1)\)\^NS;\)
\)
Basis = Table[{0, 0, 0, 0, 0, 0}, {n, 1, DimH}];

Kopplungsmatrix für Ringe

J = 
        Table[If[Abs[m - n] == 1, 1, 0], {n, 1, NS}, {m, 1, 
            NS}];
J[[1, NS]] = 1;
J[[NS, 1]] = 1;

MatrixForm[J]

Hamiltonmatrix erzeugen

Eigenwerte der Hamiltonmatrix

EV = Sort[Chop[Eigenvalues[HMatrix]]]