Exact solutions of an Ising spin chain with a spin-1 impurity
Huang Xuchu
       

(color online) The curve of specific heat C vs temperature T for different effective magnetic field λ = hx/J = 0, 5, 1.0, and 1.5. In the computation, we chose x = Jimp/J = 1, Dz = 0.4J, Dx = 1.0J, and N = 2000.