The excitation spectrum of the quantum antiferromagnetic Heisenberg model for a square lattice, as predicted from the in-phase flux state of Wigner-Jordan fermions, is calculated and compared with that calculated from conventional spin-wave theory. While the two spectra are overall similar, there are significant differences for short-wavelength excitations. The ratio of the excitation energy at wave vector (\ensuremath{\pi}/2,\ensuremath{\pi}/2) to that at (\ensuremath{\pi}/2,0) is predicted to be \ensuremath{\surd}2 for the Wigner-Jordan fermion excitations, and 2/ \ensuremath{\surd}3 for the spin-wave excitations. At long wavelengths, in contrast to the gapped spin-wave spectrum of Auerbach and Arovas, the excitation spectrum of Wigner-Jordan fermions is gapless at any finite temperature. Possible experimental examination of the excitation spectrum is discussed.
A Wigner-Jordan (WJ) transformation is used to study the one-dimensional (1D) and two-dimensional (2D) quantum antiferromagnetic Heisenberg model. The advantage of using the Wigner-Jordan transformation is that it preserves all spin-commutation relations as well as the spin on-site exclusion principle. In the 1D case a nearest-neighbor covalent-bonding state of the WJ fermions is found to have a ground-state energy (-0.4351J per site) comparable with that from the Bethe-ansatz solution (-0.4431J per site), and a linear energy spectrum at low energies with velocity 1.6366J, in close agreement with the velocity obtained by Haldane for the quantum antiferromagnetic Heisenberg model with 1/
d2 interaction (1.5708J). The method used for studying the 1D model is then applied to the 2D Heisenberg model in a square lattice. The resulting state at finite temperature is the in-phase flux state, i.e., a flux state of the Wigner-Jordan spinless fermions with an in-phase fermion orbital current circulating around each elementary plaquette.The single-particle excitation spectrum, i.e., the energy dispersion of reversing orientation of a spin in the system, of the in-phase flux state is overall similar to that of spin waves, with significant difference near the edge point k=(\ensuremath{\pi},0), at which the excitation energy of the WJ fermion is zero, whereas that of the spin-wave excitation is 2J. The specific heat of the in-phase flux state predicts a correct temperature dependence over the entire temperature range, namely a
T2 dependence at low temperature, a peak near T/J=0.6, and a 1/
T2 decreasing at high temperature. It also gives excellent agreement with the specific heat calculated from numerical methods. In contrast, the spin-wave theory only correctly predicts a
T2 dependence at low temperature. The Raman spectrum of the in-phase flux state is calculated and shows significant improvement over that from spin-wave theories compared with the experimental spectrum of
La2CuO4. The exchange parameter, J, obtained from the comparison is 1060 ${\mathrm{cm}}^{\mathrm{\ensuremath{-}}1}$, in agreement with that obtained from analyzing the neutron-scattering data.
The ground-state properties and spin excitation spectrum of the Heisenberg model are studied within the framework of the two-dimensional Jordan-Wigner theory. The mean-field state with the lowest energy is the magnetized \ensuremath{\pi}-flux state with a gap in the spin excitation spectrum. It is shown that this mean-field gap is decreased significantly by considering amplitude fluctuations in the random-phase approximation (RPA). The RPA ground-state energy and staggered magnetization are
E0=-0.332J per bond and M=0.327, respectively, which are very close to the best numerical estimates. \textcopyright{} 1996 The American Physical Society.
A disordered spin state, the nearest-neighbor virtual bonding state of Wigner-Jordan fermions, of the quantum antiferromagnetic Heisenberg model on a triangular lattice is found to have a zero-temperature energy -1.106J per spin, lower than that calculated from spin-wave theories, and in excellent agreement with that from numerical calculations, (-1.094\ifmmode\pm\else\textpm\fi{}0.018)J per spin. The excitations of the state are Wigner-Jordan fermions with a gapless spectrum. The low-temperature specific heat is predicted to be proportional to ${\mathit{T}}^{\ensuremath{\nu}}$, with \ensuremath{\nu} between 2/3 and 1, significantly deviating from the spin-wave prediction,
T2. At higher temperature, a broad peak in the specific heat near T/2J\ensuremath{\approxeq}0.75 is also predicted.
We investigate the quantum dynamics of the spins on different Heisenberg antiferromagnetic spin lattice systems. Firstly, we applied the coupled-cluster method to the spin-1/2 antiferromagnetic \emph{XXZ} model on a square lattice by employing an approximation which contains two-body long-range correlations and high-order four-body local correlations. Improvement is found for the ground-state energy, sublattice magnetization, and the critical anisotropy when comparing with the approximation including the two-body correlations alone. We also obtain the full excitation spectrum which is in good agreement with the quantum Monte Carlo results and the high-order spin-wave theory.Secondly, we study the longitudinal excitations of quantum antiferromagnets on a triangular lattice by a recently proposed microscopic many-body approach based on magnon-density waves. We calculate the full longitudinal excitation spectra of the antiferromagnetic Heisenberg model for a general spin quantum number in the isotropic limit. Similar to the square lattice model, we find that, at the center of the first hexagonal Brillouin zone
Γ(q=0) and at the magnetic ordering wavevectors
±[Q=(4π/3,0)], the excitation spectra become gapless in the thermodynamic limit, due to the slow, logarithmic divergence of the structure factor. However, these longitudinal modes on two-dimensional models may be considered as quasi-gapped, as any finite-size effect or small anisotropy will induce a large energy gap, when compared with the counterpart of the transverse spin-wave excitations. We have also investigated the excited states of the quasi-one-dimensional quantum antiferromagnets on hexagonal lattices, including the longitudinal modes based on the magnon-density waves. A model Hamiltonian with a uniaxial single-ion anisotropy is first studied by a spin-wave theory based on the one-boson method; the ground state thus obtained is employed for the study of the longitudinal modes. The full energy spectra of both the transverse modes (i.e., magnons) and the longitudinal modes are obtained as functions of the nearest-neighbor coupling and the anisotropy constants. We have found two longitudinal modes due to the non-collinear nature of the triangular antiferromagnetic order, similar to that of the phenomenological field theory approach by Affleck. The excitation energy gaps due to the anisotropy and the energy gaps of the longitudinal modes without anisotropy are then investigated. We then compares our results for the longitudinal energy gaps at the magnetic wavevectors with the experimental results for several antiferromagnetic compounds with both integer and non-integer spin quantum numbers, and we find good agreements after the higher-order contributions are included in our calculations.
We study the excitation spectra of two-dimensional spin-
1/2 Heisenberg model with checkerboard structures using stochastic analytic continuation of imaginary-time correlation function obtained from quantum Monte Carlo simulation. The checkerboard models have two different antiferromagnetic nearest-neighbor interactions
J1 and
J2, and the tuning parameter
g is defined as
J2/J1. The dynamic structure factors are systematically calculated in all phases of the models as well as at the critical points. To give a full understanding of the dynamic spectra, wave theory is employed to explain some features of numerical results, especially for the low-energy part. When
g is close to
1, the features of the excitation spectra of each checkerboard model are roughly the same as those of the original square lattice antiferromagnetic Heisenberg model, and the high-energy continuum among them are discussed. In contrast to the other checkerboard structures investigated in this paper, the
3×3 checkerboard model has distinctive excitation features, like a gap between a low-energy gapless branch and a gapped high-energy part exists when
g is small. The gapless branch in this case can be regarded as wave in N
eˊel order formed by block spin in each
3×3 plaquette with effective exchange interaction originated from renormalization. One finding is that the continuum also appears in this low-energy branch, which may originated from nearly deconfined process in block spin.
We study the spin excitation spectra of the two-dimensional spin-
21 Heisenberg model with a checkerboard structure using stochastic analytic continuation of the imaginary-time correlation function obtained from a quantum Monte Carlo simulation. The checkerboard models have two different antiferromagnetic nearest-neighbor interactions
J1 and
J2, and the tuning parameter
g is defined as
J2/J1. The dynamic spin structure factors are systematically calculated in all phases of the models as well as at the critical points. To give a full understanding of the dynamic spectra, spin wave theory is employed to explain some features of numerical results, especially for the low-energy part. When
g is close to 1, the features of the spin excitation spectra of each checkerboard model are roughly the same as those of the original square lattice antiferromagnetic Heisenberg model, and the high-energy continuum among them is discussed. In contrast to the other checkerboard structures investigated in this paper, the $3\ifmmode\times\else\texttimes\fi{}3$ checkerboard model has distinctive excitation features, such as a gap between a low-energy gapless branch and a gapped high-energy part that exists when
g is small. The gapless branch in this case can be regarded as a spin wave in N\'eel order formed by a ``block spin'' in each $3\ifmmode\times\else\texttimes\fi{}3$ plaquette with an effective exchange interaction originating from renormalization. One unexpected finding is that the continuum also appears in this low-energy branch.