Dear itensor

I am working on finite temperature spin-1 Heisenberg model including single-ion anisotropy term

H=J\sum*{i,j} S}*i . S*j -K \sum*i (S_i^z)^2

To catch the correct physics at temperature at T-->0, it seems one needs to use a kind of "soft edges" trick, to approximate the infinite system.

I mean, I want to have a chain of spin-1 excepts those are sitting at both edges should be spin1/2.

S1/2---S1---S1---S1. . . S1---S1--S1/2

any suggestions are highly appreciated.

Regards

Javad