No Arabic abstract
Searching for spin liquids on the honeycomb J1-J2 Heisenberg model has been attracting great attention in the past decade. In this Paper we investigate the topological properties of the J1-J2 Heisenberg model by introducing nearest-neighbour and next-nearest-neighbour bond parameters. We find that there exist two topologically different phases in the spin disordered regime 0.2<J2/J1<0.5: for J2/J1<0.32, the system is a zero-flux spin liquid which is topological trivial and gapless; for J2/J1>0.32, it is a pi-flux chiral spin liquid, which is topological nontrivial and gapped. These results suggest that there exist two topologically different spin disorder phases in honeycomb J1-J2 Heisenberg model.
We perform an extensive density matrix renormalization group (DMRG) study of the ground-state phase diagram of the spin-1/2 J_1-J_2 Heisenberg model on the kagome lattice. We focus on the region of the phase diagram around the kagome Heisenberg antiferromagnet, i.e., at J_2=0. We investigate the static spin structure factor, the magnetic correlation lengths, and the spin gaps. Our results are consistent with the absence of magnetic order in a narrow region around J_2approx 0, although strong finite-size effects do not allow us to accurately determine the phase boundaries. This result is in agreement with the presence of an extended spin-liquid region, as it has been proposed recently. Outside the disordered region, we find that for ferromagnetic and antiferromagnetic J_2 the ground state displays signatures of the magnetic order of the sqrt{3}timessqrt{3} and the q=0 type, respectively. Finally, we focus on the structure of the entanglement spectrum (ES) in the q=0 ordered phase. We discuss the importance of the choice of the bipartition on the finite-size structure of the ES.
The two dimensional Heisenberg antiferromagnet on the square lattice with nearest (J1) and next-nearest (J2) neighbor couplings is investigated in the strong frustration regime (J2/J1>1/2). A new effective field theory describing the long wavelength physics of the model is derived from the quantum hamiltonian. The structure of the resulting non linear sigma model allows to recover the known spin wave results in the collinear regime, supports the presence of an Ising phase transition at finite temperature and suggests the possible occurrence of a non-magnetic ground state breaking rotational symmetry. By means of Lanczos diagonalizations we investigate the spin system at T=0, focusing our attention on the region where the collinear order parameter is strongly suppressed by quantum fluctuations and a transition to a non-magnetic state occurs. Correlation functions display a remarkable size independence and allow to identify the transition between the magnetic and non-magnetic region of the phase diagram. The numerical results support the presence of a non-magnetic phase with orientational ordering.
Based on the mapping between $s=1/2$ spin operators and hard-core bosons, we extend the cluster perturbation theory to spin systems and study the whole excitation spectrum of the antiferromagnetic $J_{1}$-$J_{2}$ Heisenberg model on the square lattice. In the Neel phase for $J_{2}lesssim0.4J_{1}$, in addition to the dominant magnon excitation, there is an obvious continuum close to $(pi,0)$ in the Brillouin zone indicating the deconfined spin-1/2 spinon excitations. In the stripe phase for $J_{2}gtrsim0.6J_{1}$, we find similar high-energy two-spinon continuums at $(pi/2,pi/2)$ and $(pi/2,pi)$, respectively. The intermediate phase is characterized by a spectrum with completely deconfined broad continuum, which is attributed to a $Z_{2}$ quantum spin liquid with the aid of a variational-Monte-Carlo analysis.
We study the plaquette valence-bond solid phase of the spin-1/2 J_1-J_2 antiferromagnet Heisenberg model on the square lattice within the bond-operator theory. We start by considering four S = 1/2 spins on a single plaquette and determine the bond operator representation for the spin operators in terms of singlet, triplet, and quintet boson operators. The formalism is then applied to the J_1-J_2 model and an effective interacting boson model in terms of singlets and triplets is derived. The effective model is analyzed within the harmonic approximation and the previous results of Zhitomirsky and Ueda [Phys. Rev. B 54, 9007 (1996)] are recovered. By perturbatively including cubic (triplet-triplet-triplet and singlet-triplet-triplet) and quartic interactions, we find that the plaquette valence-bond solid phase is stable within the parameter region 0.34 < J_2/J_1 < 0.59, which is narrower than the harmonic one. Differently from the harmonic approximation, the excitation gap vanishes at both critical couplings J_2 = 0.34 J_1 and J_2 = 0.59 J_1. Interestingly, for J_2 < 0.48 J_1, the excitation gap corresponds to a singlet-triplet excitation at the $Gamma$ point while, for J_2 > 0.48 J_1, it is related to a singlet-singlet excitation at the X = (pi/2,0) point of the tetramerized Brillouin zone.
We study the spin-1/2 Heisenberg model on the square lattice with first- and second-neighbor antiferromagnetic interactions J1 and J2, which possesses a nonmagnetic region that has been debated for many years and might realize the interesting Z2 spin liquid. We use the density matrix renormalization group approach with explicit implementation of SU(2) spin rotation symmetry and study the model accurately on open cylinders with different boundary conditions. With increasing J2, we find a Neel phase, a plaquette valence-bond (PVB) phase with a finite spin gap, and a possible spin liquid in a small region of J2 between these two phases. From the finite-size scaling of the magnetic order parameter, we estimate that the Neel order vanishes at J2/J1~0.44. For 0.5<J2/J1<0.61, we find dimer correlations and PVB textures whose decay lengths grow strongly with increasing system width, consistent with a long-range PVB order in the two-dimensional limit. The dimer-dimer correlations reveal the s-wave character of the PVB order. For 0.44<J2/J1<0.5, spin order, dimer order, and spin gap are small on finite-size systems and appear to scale to zero with increasing system width, which is consistent with a possible gapless SL or a near-critical behavior. We compare and contrast our results with earlier numerical studies.