スピンハミルトニアンの計算

＃見苦しくてすみません。

等幅のフォントで見て下さい。''''は桁揃えのために入れてます。

ブラを行ベクトル、演算子を行列、ケットを列ベクトルで表します。

<+1| = [ 1, 0, 0]

< 0| = [ 0, 1, 0]

<-1| = [ 0, 0, 1]

''''┌ 0, 1, 0 ┐
Sx =│ 1, 0, 1 │/√2
''''└ 0, 1, 0 ┘

''''┌ 0,-i, 0 ┐
Sy =│ i, 0,-i │/√2
''''└ 0, i, 0 ┘

''''┌ 1, 0, 0 ┐
Sz =│ 0, 0, 0 │
''''└ 0, 0,-1 ┘

''''''┌ 1 ┐
|+1> =│ 0 │
''''''└ 0 ┘

''''''┌ 0 ┐
|0 > =│ 1 │
''''''└ 0 ┘

''''''┌ 0 ┐
|-1> =│ 0 │
''''''└ 1 ┘

そうすれば、<+1|H|+1>などの行列要素は、行列のかけ算から求まります。
例えば <+1|DzSz^2|+1> なら
''''''''''''''''''''''''''''''┌ 1, 0, 0 ┐┌ 1, 0, 0 ┐┌ 1 ┐
<+1|DzSz^2|+1> = Dz [ 1, 0, 0]│ 0, 0, 0 ││ 0, 0, 0 ││ 0 │= Dz
''''''''''''''''''''''''''''''└ 0, 0,-1 ┘└ 0, 0,-1 ┘└ 0 ┘
です。
あとはハミルトニアンを3×3の行列で書いて、
S・H = Sx Hx + Sy Hy + Sz Hz と Dx + Dy + Dz = 0
の関係を使えば解けると思います。