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 variational study of the ground state of the spin
−1/2 anisotropic Heisenberg antiferromagnet has been revisited on a square lattice by improving and correcting past numerical results found in Sol. State. Comm. 165, 33 (2013). The Hamiltonian has been implemented on a square lattice with antiferromagnetic interactions between nearest- and next-nearest neighbors. The nearest-neighbor couplings have different strengths, namely,
J1 and
J1′, for the x and y directions, respectively. These couplings compete with the next-nearest ones denoted by
J2. We obtained a new phase diagram in the
λ−α plane, where
λ=J1′/J1 and
α=J2/J1, whose topology is slightly different of that previously found. There is no direct frontier dividing the collinear (CAF) and the antiferromagnetic order (AF), rather, the quantum paramagnetic phase (QP) separates these two phases for all positive values of
λ and
α. The true nature of the frontiers has been obtained by scanning rigorously the relevant points of the
λ−α plane.
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 zero-temperature phase diagram and the low-lying excitations of a square-lattice spin-half Heisenberg antiferromagnet with two types of regularly distributed nearest-neighbor exchange bonds
[J>0 (antiferromagnetic) and $\ensuremath{-}\ensuremath{\infty}<{J}^{\ensuremath{'}}<\ensuremath{\infty}]$ using the coupled cluster method (CCM) for high orders of approximation (up to LSUB8). We use a N\'eel model state as well as a helical model state as a starting point for the CCM calculations. We find a second-order transition from a phase with N\'eel order to a finite-gap quantum disordered phase for sufficiently large antiferromagnetic exchange constants ${J}^{\ensuremath{'}}>0.$ For frustrating ferromagnetic couplings ${J}^{\ensuremath{'}}<0$ we find indications that quantum fluctuations favor a first-order phase transition from the N\'eel order to a quantum helical state, by contrast with the corresponding second-order transition in the corresponding classical model. The results are compared to those of exact diagonalizations of finite systems (up to 32 sites) and those of spin-wave and variational calculations. The CCM results agree well with the exact diagonalization data over the whole range of the parameters. The special case of ${J}^{\ensuremath{'}}=0,$ which is equivalent to the honeycomb lattice, is treated more closely.
We study the spin-1/2 and spin-1 [Formula: see text] Heisenberg antiferromagnets (HAFs) on an infinite, anisotropic, two-dimensional triangular lattice, using the coupled cluster method. With respect to an underlying square-lattice geometry the model contains antiferromagnetic (J1 > 0) bonds between nearest neighbours and competing [Formula: see text] bonds between next-nearest-neighbours across only one of the diagonals of each square plaquette, the same diagonal in each square. In a topologically equivalent triangular-lattice geometry the model has two sorts of nearest-neighbour bonds, with [Formula: see text] bonds along parallel chains and with J1 bonds providing an interchain coupling. The model thus interpolates between an isotropic HAF on the square lattice at one extreme (κ = 0) and a set of decoupled chains at the other (κ → ∞), with the isotropic HAF on the triangular lattice in between at κ = 1. For the spin-1/2 [Formula: see text] model, we find a weakly first-order (or possibly second-order) quantum phase transition from a Néel-ordered state to a helical state at a first critical point at κc1 = 0.80 ± 0.01, and a second critical point at κc2 = 1.8 ± 0.4 where a first-order transition occurs between the helical state and a collinear stripe-ordered state. For the corresponding spin-1 model we find an analogous transition of the second-order type at κc1 = 0.62 ± 0.01 between states with Néel and helical ordering, but we find no evidence of a further transition in this case to a stripe-ordered phase.
The ruby lattice has been the subject of much interest recently due its realization in Rydberg atom arrays, where its rich variety of frustrated interactions gives rise to topologically ordered quantum spin liquids. Similarly, numerical studies of ruby-lattice spin models, with both isotropic and anisotropic interactions, have provided evidence of gapped and gapless spin-liquid ground states with different low-energy gauge structures. Motivated by these findings, we perform a projective symmetry group (PSG) classification of U(1) and
Z2 fermionic spinon mean-field theories$\unicode{x2014}$respecting space-group and time-reversal symmetries$\unicode{x2014}$for
S=1/2 spins. We obtain a total of 50 U(1) and 64
Z2 PSGs, and upon restricting their realization via mean-field
Ansa¨tze with up to second-nearest-neighbor singlet amplitudes (relevant to the models studied here), only 8 U(1) and 18
Z2 distinct states are obtained. We present the singlet fields for all
Ansa¨tze up to third-nearest-neighboring bonds and discuss their spinon dispersions as well as their dynamical spin structure factors. Building on this information, we also obtain the phase diagram of the Heisenberg model in the presence of first (
J1), second (
J1′), and third (
J2) neighbor antiferromagnetic couplings within a self-consistent mean-field approximation.