Recently there has been extensive theoretical interest in the remarkable Haldane prediction that spin-1 antiferromagnetic Heisenberg-like chains should display an energy gap between singlet ground state and excited continua. This behavior has been confirmed numerically, and it has been shown that the energy gap is of a magnitude to be observable experimentally. The phenomenon has captured the interest of experimental groups in neutron scattering and magnetism. A recent experimental study obtains very good experimental agreement with the theoretical T=0 magnetization curve obtained numerically by Parkinson and Bonner.1,2 An explanation will be provided for the unusual phase diagram and other anomalous features observed by the authors of Ref. 2. Numerical susceptibility studies of Heisenberg antiferromagnetic and general XXZ spin-1 chains will be presented, and discussed in terms of a novel quantum-classical cross-over phenomenon.2 Finally, a comparison will be made of T=0 magnetization isotherms for several spin-1 chains, including the Takhtajan–Babujian chain and the chain with pure biquadratic interactions.
Recently there has been extensive theoretical interest in the remarkable Haldane prediction that spin-1 antiferromagnetic Heisenberg-like chains should display an energy gap between singlet ground state and excited continua. This behavior has been confirmed numerically, and it has been shown that the energy gap is of a magnitude to be observable experimentally. The phenomenon has captured the interest of experimental groups in neutron scattering and magnetism. A recent experimental study obtains very good experimental agreement with the theoretical T=0 magnetization curve obtained numerically by Parkinson and Bonner.1,2 An explanation will be provided for the unusual phase diagram and other anomalous features observed by the authors of Ref. 2. Numerical susceptibility studies of Heisenberg antiferromagnetic and general XXZ spin-1 chains will be presented, and discussed in terms of a novel quantum-classical cross-over phenomenon.2 Finally, a comparison will be made of T=0 magnetization isotherms for several spin-1 chains, including the Takhtajan–Babujian chain and the chain with pure biquadratic interactions.
Since the Heisenberg spin chain can be considered the simplest realistic model of magnetism, surprise and some degree of controversy have resulted from recent work of Haldane. The prediction is that quantum spin chains with half-integer spin should all display T=0 phase behavior equivalent to that of the Bethe Ansatz integrable (solvable) spin-1/2 quantum chain. More remarkably, the class of integer spin chains is predicted to show very different phase behavior. In particular, a gap should be present in the spectrum of a Heisenberg antiferromagnetic chain. This remarkable feature is counterintuitive in terms of accepted wisdom in magnetism (spin-wave theory, spin-Peierls theory) and critical phenomena. Consequently the vertification of the prediction is of great interest. A considerable amount of numerical work has been done, involving finite-chain, finite-size scaling, variational, Monte Carlo and other calculations, which will be reviewed here. The present consensus is that the weight of numerical evidence supports the prediction, although puzzling features still remain. Adding additional interactions to the basic Heisenberg Hamiltonian such as spin (XXZ) anisotropy, single-ion anisotropy, biquadratic exchange, and an applied magnetic field, generates a rich and complicated phase diagram for chains with spin >1/2, particularly for the case of integer spin. The s=1 phase diagram seems to display critical behavior of a type not previously encountered. A theoretical appraisal of the Haldane phenomenon will include a discussion of the possible role of nonintegrability. Mention will also be made of current progress in experimental investigation of the phenomenon, including problems that might be encountered. More recent work of Affleck has greatly generalized the field-theoretic mappings which underlay the original work of Haldane. A number of interesting problems have been mapped into quantum spin chains of various types, including field theoretic phenomena and the localization problem of the quantum Hall effect.
The spectral properties, and associated critical singularities, of spin-1 Heisenberg antiferromagnetic chains are of high current interest, both theoretical and experimental. Theoretical studies using the single mode approximation and also a Monte Carlo calculation indicate that the spectrum has a gap at k=π, as predicted by Haldane. They also predict the unusual feature that the excited state continuum is asymmetric about the mid-point of the Brillouin zone, k=π/2. However, neither approach is sufficiently accurate to predict reliably the magnitude of the gap at k=0 which was not discussed by Haldane. Our approach is to use static (spectral) calculations of excited states of chains of N=4 through 14 spins for all values of k, a novel feature not associated with previous finite-size calculations, to discuss the dispersion spectra of the S=1 Heisenberg chain. In addition, we have calculated dynamical properties, in particular S(q,ω), exactly for finite chains up to N=14 spins. This allows us to estimate the gap at k=0 to good accuracy, and we find that to within this accuracy the gap at k=0 is twice the Haldane gap at k=π. This result is important in connection with very recent neutron scattering and magnetization measurements. In addition, we have calculated the correlation length and critical singularities associated with dynamical quantities.
The antiferromagnetic (AFM) s=1 chain with biquadratic exchange coupling has recently been the subject of considerable attention, and some controversy. Theoretical arguments by Affleck that the model has dimerized character associated with a small spectral excitation gap have not received unanimous support from numerical calculations. For finite systems the field-dependent dispersion spectra are interesting and informative concerning the nature of the model. The excitations divide into two classes having very different character, corresponding to SzT even and SzT odd. Reflection symmetry about k=π/2 is apparent if states with a given total spin ST, rather than SzT, are plotted as a function of wavevector k, consistent with the Affleck prediction of dimerized character. The high degeneracy of the eigenstates reflects the existence of new symmetry operators which (a) shift the SzT value, (b) change wave vector k to π/2−k, and (c) shift the total spin quantum ST, of the energy eigenstates. Parkinson has recently observed that finite N states with SzT even and SzT≥2 map into states of the finite N s= 1/2 XXZ model at a special value of uniaxial anisotropy. Barber and Batchelor (unpublished) have shown the existence of an analytic mapping between states of the s=1 biquadratic chain and the 2D 9-state Potts model. They further show that all states of the Parkinson s= 1/2 XXZ model map into states of the s=1 biquadratic model. Hence, they find the ground state and first excited state energies exactly, and verify the Affleck picture. Their approach falls short of establishing complete integrability, however. Numerical studies on finite chains indicate that the T=0 magnetization curve is determined by the lowest-energy states for given SzT even, which are states of the s= 1/2 XXZ chain. Hence the magnetization curve of the s=1 AFM biquadratic chain is the same as the analytically known magnetization curve of the s= 1/2 XXZ chain!
Extensive numerical studies have been performed on Heisenberg antiferromagnetic chains of spin (1/2) (up to N=20), spin 1 (N=14), spin (3/2) (N=10), and spin 2 (N=8). With use of the Lancz\"os technique, primarily, the two lowest-lying eigenvalues have been calculated for all values of wave vector q and all values of magnetization (STz) up to saturation for each chain. From a knowledge of the eigenvector corresponding to the lowest eigenstate for each STz, the T=0 spin-pair correlation functions have also been calculated as a function of field. We find a most unusual quantum-classical crossover phenomenon. It shows up in greatest detail in the field-dependent dispersion spectra, but the consequences are consistently manifested in the behavior of the T=0 magnetization isotherms and in the correlations both in real space and the Fourier transforms in q space. The additional data relevant to behavior in a field have allowed us to extend previous numerical studies of Heisenberg antiferromagnetic chains with higher spin whose purpose was to examine the validity of the Haldane conjecture. The Haldane conjecture implies that Heisenberg antiferromagnetic chains with integer spin have a gap in their excitation spectrum whereas chains with half-integer spin do not. While no feature of our extended investigations is in conflict with the conjecture, unusual features associated apparently with very slow convergence make the outcome less than conclusive. It appears that calculations on significantly longer chains are required to observe with confidence the large-N asymptotic limiting behavior.
In a previous paper (Liu et al 2012 Phys. Rev. B 85 195144 ), a variational Monte Carlo method (based on Gutzwiller projected states) was generalized to S = 1 systems. This method provided very good trial ground states for the gapped phases of an S = 1 bilinear-biquadratic (BLBQ) Heisenberg chain. In this paper, we extend the approach to study the low-lying elementary excitations in S = 1 chains. We calculate the one-magnon and two-magnon excitation spectra of the BLBQ Heisenberg chain and the results agree very well with recent data in the literature. In our approach, the difference of the excitation spectrum between the Haldane phase and the dimer phase (such as the even/odd size effect) can be understood from their different topologies of the corresponding mean field theory. We especially study the Takhtajan–Babujian critical point. Despite the fact that the 'elementary excitations' are spin-1 magnons, which are different from the spin-1/2 spinons in Bethe solution, we show that the excitation spectrum, critical exponent ( ) and central charge ( c = 1.45) calculated from our theory agree well with the Bethe ansatz solution and conformal field theory predictions.