Matrix product states and exactly solvable spin 1/2 Heisenberg chains with nearest neighbor interactions

Using the matrix product formalism, we introduce a two parameter family of exactly solvable $xyz$ spin 1/2 Heisenberg chains in magnetic field (with nearest neighbor interactions) and calculate the ground state and correlation functions in compact form. The ground state has a very interesting property: all the pairs of spins are equally entangled with each other. Therefore it is possible to engineer long-range entanglement in experimentally realizable spin systems on the one hand and study more closely quantum phase transition in such systems on the other.