  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.
