We calculate the quantum phase diagram of an XXZ chain with nearest-neighbor (NN)
J1 and next-NN exchange
J2 with anisotropies ${\ensuremath{\Delta}}_{1}$ and ${\ensuremath{\Delta}}_{2},$ respectively. In particular we consider the case ${\ensuremath{\Delta}}_{1}=\ensuremath{-}{\ensuremath{\Delta}}_{2},$ to interpolate between the XX chain $({\ensuremath{\Delta}}_{i}=0)$ and the isotropic model with ferromagnetic phase
J2. For ${\ensuremath{\Delta}}_{1}<\ensuremath{-}1,$ a ferromagnetic phase and two antiferromagnetic phases exist. For $|{\ensuremath{\Delta}}_{i}|<1,$ the boundary between the dimer and spin fluid phases is determined by the method of crossing of excitation spectra. For large
J2/J1, this method seems to indicate the existence of a second spin fluid critical phase. However, an analysis of the spin stiffness and magnetic susceptibility for ${\ensuremath{\Delta}}_{1}={\ensuremath{\Delta}}_{2}=1$ suggests that a small gap is present.
The Hamiltonian of one-dimensional integrable multi-impurity XXZ spin chain with open boundary cases is constructed by using the method of adding parameters to the vertical matrices. This is the first approach to add multi-impurity and next-nearest-neighbor interaction to the integrable anisotropic Heisenberg spin chain. It is found that all of the impurity parameters in the bulk must be pure imaginary numbers if the cross parameter is real for getting the hermitian Hamiltonian, or vice versa, when the next-nearest-neighbor interactions are introduced. The Hamiltonian can be also applied to describe the many body systems with two kinds of particles if all of the impurity parameters are equal.
The one dimensional S=1 XXZ model with next-nearest-neighbor interaction
α and Ising-type anisotropy
Δ is studied by using a numerical diagonalization technique. We discuss the ground state phase diagram of this model numerically by the twisted-boundary-condition level spectroscopy method and the phenomenological renormalization group method, and analytically by the spin wave theory. We determine the phase boundaries among the XY phase, the Haldane phase, the ferromagnetic phase and the Néel phase, and then we confirm the universality class. Moreover, we map this model onto the non-linear
σ model and analyze the phase diagram in the
α ≪ -1 and
Δ ∼ 1 region by using the renormalization group method.