The nature of quantum spin liquids is studied for the spin-
1/2 antiferromagnetic Heisenberg model on a square lattice containing exchange interactions between nearest-neighbor sites,
J1, and those between next-nearest-neighbor sites,
J2. We perform variational Monte Carlo simulations together with the quantum-number-projection technique and clarify the phase diagram in the ground state together with its excitation spectra. We obtain the nonmagnetic phase in the region
0.4<J2/J1≤0.6 sandwiched by the staggered and stripe antiferromagnetic (AF) phases. Our direct calculations of the spin gap support the notion that the triplet excitation from the singlet ground state is gapless in the range of
0.4<J2/J1≤0.5, while the gapped valence-bond-crystal (VBC) phase is stabilized for
0.5<J2/J1≤0.6. The VBC order is likely to have the columnar symmetry with a spontaneous symmetry breaking of the
C4v symmetry. The power-law behaviors of the spin-spin and dimer-dimer correlation functions in the gapless region are consistent with the emergence of the algebraic quantum-spin-liquid phase (critical phase). The exponent of the spin correlation
⟨S(0)S(r)⟩∝1/rz+η at a long distance
r appears to increase from
z+η∼1 at
J2/J1∼0.4 toward the continuous transition to the VBC phase at
J1/J1∼0.5. Our results, however, do not fully exclude the possibility of a direct quantum transition between the staggered AF and VBC phases with a wide critical region and deconfined criticality.
We study the spin-
21 Heisenberg model on the triangular lattice with antiferromagnetic first-
(J1) and second-
(J2) nearest-neighbor interactions using density matrix renormalization group. By studying the spin correlation function, we find a ${120}^{\ensuremath{\circ}}$ magnetic order phase for ${J}_{2}\ensuremath{\lesssim}0.07{J}_{1}$ and a stripe antiferromagnetic phase for ${J}_{2}\ensuremath{\gtrsim}0.15{J}_{1}$. Between these two phases, we identify a spin-liquid region characterized by exponential decaying spin and dimer correlations, as well as large spin singlet and triplet excitation gaps on finite-size systems. We find two near degenerating ground states with distinct properties in two sectors, which indicates more than one spin-liquid candidate in this region. While the sector with spinons is found to respect time reversal symmetry, the even sector without spinons breaks such a symmetry for finite-size systems. Furthermore, we detect the signature of the fractionalization by following the evolution of different ground states with inserting spin flux into the cylinder system. Moreover, by tuning the anisotropic bond coupling, we explore the nature of the spin-liquid phase and find the optimal parameter region for gapped
Z2 spin liquids.
The linear spin-
21 Heisenberg antiferromagnet with exchanges
J1 and
J2 between first and second neighbors has a bond-order wave (BOW) phase that starts at the fluid-dimer transition at
J2/J1=0.2411 and is particularly simple at
J2/J1=1/2. The BOW phase has a doubly degenerate singlet ground state, broken inversion symmetry, and a finite-energy gap
Em to the lowest-triplet state. The interval
0.4<J2/J1<1.0 has large
Em and small finite-size corrections. Exact solutions are presented up to
N=28 spins with either periodic or open boundary conditions and for thermodynamics up to
N=18. The elementary excitations of the BOW phase with large
Em are topological spin-
21 solitons that separate BOWs with opposite phase in a regular array of spins. The molar spin susceptibility ${\ensuremath{\chi}}_{M}(T)$ is exponentially small for
T⪡Em and increases nearly linearly with
T to a broad maximum.
J1 and
J2 spin chains approximate the magnetic properties of the BOW phase of Hubbard-type models and provide a starting point for modeling alkali-tetracyanoquinodimethane salts.
We study Heisenberg spin-1/2 and spin-1 chains with alternating ferromagnetic (
J1F) and antiferromagnetic (
J1A) nearest-neighbor interactions and ferromagnetic next-nearest-neighbor interaction (
J2F). In this model frustration is induced by the non-zero
J2F and the effect of dimerization is realized by changing the strength of
J1A for a fixed value of
J1F. The model with site spin
s behaves like a Haldane spin chain with site spin 2
s in the limit of vanishing
J2F and large
J1F/J1A. We show that the exact ground state of the model can be found along a line in the parameter space. The phase diagram in the space of
J1A−J2F is determined using numerical techniques complemented by some analytical calculations. A number of quantities, including the ground state spin, energy gap, structure factor, entanglement entropy, and zero temperature magnetization, are studied to understand the complete phase diagram. An interesting and application-wise important feature of this model is that it exhibits a macroscopic magnetization jump in the presence of a magnetic field; we study this using an effective Hamiltonian.
The density matrix renormalization group and quantum Monte Carlo methods are used to describe coupled trimer chains in a magnetic field
h. The Hamiltonian contains exchange terms involving the intratrimer coupling
J1 (taken as the unit of energy) and the intertrimer coupling
J2, plus the Zeeman interaction for a magnetic field
h along the
z direction. Results for the magnetization per trimer
m are calculated in regimes of positive and negative values of the ratio
J=J2/J1, from which the rich field-induced ground state phase diagram
h versus
J is derived, with the presence of a Luttinger liquid, the 1/3 plateau (
m=1/2), and the one of fully polarized magnetization (
m=3/2). Also, zero-field Lanczos calculation of the spin-wave dispersion from the 1/3 plateau for
Sz=1 is shown at the previous regimes of
J values. In addition, we also report on the decay of correlation functions of trimers along open chains, as well as the average two-magnon distribution. The ground state is ferrimagnetic for $0<J\ensuremath{\le}1$, and is a singlet for $\ensuremath{-}1\ensuremath{\le}J<0$. In the singlet phase, the spin correlation functions along the legs present an antiferromagnetic power-law decay, similar to the spin-1/2 linear chain, thus suggesting that the ground state is made of three coupled antiferromagnetically oriented chains. In the singlet phase, the dimensionless thermal magnetic susceptibility per site normalized by
1/∣J∣ gets closer to $1/{\ensuremath{\pi}}^{2}$ as the temperature $T\ensuremath{\rightarrow}0$. For the ferrimagnetic phase, we fit the susceptibility to the experimental data for the compound
Pb3Cu3(PO4)4 and estimate the model exchange couplings:
J1=74.8K and
J=0.4. These values imply a range of energies for the magnon excitations that are in accord with the data from neutron scattering experiments on
Pb3Cu3(PO4)4 for two excitation modes. The 1/3 plateau closes only at
1/∣J∣=0 with
J<0.
We analyze the stability of the spontaneously dimerized phase of the frustrated Heisenberg antiferromagnet---the ${J}_{1}{\ensuremath{-}J}_{2}$ model---in two dimensions. The lowest triplet excitation, corresponding to breaking of a singlet bond, is found to be stable in the region $0.38\ensuremath{\lesssim}{J}_{2}{/J}_{1}\ensuremath{\lesssim}0.62.$ In addition we find a stable low-energy collective singlet mode, which reflects the spontaneous violation of the discrete symmetry. The spontaneous dimerization vanishes at the point of second-order quantum transition into the N\'eel ordered phase $[{(J}_{2}{/J}_{1}{)}_{c1}\ensuremath{\approx}0.38].$ We argue that the disappearance of dimer order is related to the vanishing of the singlet energy gap at the transition point.