No Arabic abstract
Using variational wave functions and Monte Carlo techniques, we study the antiferromagnetic Heisenberg model with first-neighbor $J_1$ and second-neighbor $J_2$ antiferromagnetic couplings on the honeycomb lattice. We perform a systematic comparison of magnetically ordered and nonmagnetic states (spin liquids and valence-bond solids) to obtain the ground-state phase diagram. Neel order is stabilized for small values of the frustrating second-neighbor coupling. Increasing the ratio $J_2/J_1$, we find strong evidence for a continuous transition to a nonmagnetic phase at $J_2/J_1 approx 0.23$. Close to the transition point, the Gutzwiller-projected uniform resonating valence bond state gives an excellent approximation to the exact ground-state energy. For $0.23 lesssim J_2/J_1 lesssim 0.4$, a gapless $Z_2$ spin liquid with Dirac nodes competes with a plaquette valence-bond solid. In contrast, the gapped spin liquid considered in previous works has significantly higher variational energy. Although the plaquette valence-bond order is expected to be present as soon as the Neel order melts, this ordered state becomes clearly favored only for $J_2/J_1 gtrsim 0.3$. Finally, for $0.36 lesssim J_2/J_1 le 0.5$, a valence-bond solid with columnar order takes over as the ground state, being also lower in energy than the magnetic state with collinear order. We perform a detailed finite-size scaling and standard data collapse analysis, and we discuss the possibility of a deconfined quantum critical point separating the Neel antiferromagnet from the plaquette valence-bond solid.
We consider the quasi-two-dimensional pseudo-spin-1/2 Kitaev - Heisenberg model proposed for A2IrO3 (A=Li, Na) compounds. The spin-wave excitation spectrum, the sublattice magnetization, and the transition temperatures are calculated in the random phase approximation (RPA) for four different ordered phases, observed in the parameter space of the model: antiferomagnetic, stripe, ferromagnetic, and zigzag phases. The N{e}el temperature and temperature dependence of the sublattice magnetization are compared with the experimental data on Na2IrO3.
We introduce a simple model of SO($N$) spins with two-site interactions which is amenable to quantum Monte-Carlo studies without a sign problem on non-bipartite lattices. We present numerical results for this model on the two-dimensional triangular lattice where we find evidence for a spin nematic at small $N$, a valence-bond solid (VBS) at large $N$ and a quantum spin liquid at intermediate $N$. By the introduction of a sign-free four-site interaction we uncover a rich phase diagram with evidence for both first-order and exotic continuous phase transitions.
We study the spin-$1/2$ Heisenberg model on the triangular lattice with the antiferromagnetic first ($J_1$) and second ($J_2$) nearest-neighbor interactions using density matrix renormalization group. By studying the spin correlation function, we find a $120^{circ}$ magnetic order phase for $J_2 lesssim 0.07 J_1$ and a stripe antiferromagnetic phase for $J_2 gtrsim 0.15 J_1$. Between these two phases, we identify a spin liquid region characterized by the exponential decaying spin and dimer correlations, as well as the 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 candidates in this region. While the sector with spinon is found to respect the time reversal symmetry, the even sector without a spinon 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 the gapped $Z_2$ spin liquid.
The honeycomb lattice material Li2RuO3 undergoes a dimerization of Ru4+ cations on cooling below 270C, where the magnetic susceptibility vanishes. We use density functional theory calculations to show that this reflects the formation of a valence bond crystal, with a strong bond disproportionation. On warming, x-ray diffraction shows that discrete three-fold symmetry is regained on average, and the dimerization apparently disappears. In contrast, local structural measurements using high-energy x-rays, show that disordered dimers survive at the nanoscale up to at least 650C. The high temperature phase of Li2RuO3 is thus an example of a valence bond liquid, where thermal fluctuations drive resonance between different dimer coverages, a classic analogue of the resonating valence bond state often discussed in connection with high T$_c$ cuprates.
Strongly correlated systems with geometric frustrations can host the emergent phases of matter with unconventional properties. Here, we study the spin $S = 1$ Heisenberg model on the honeycomb lattice with the antiferromagnetic first- ($J_1$) and second-neighbor ($J_2$) interactions ($0.0 leq J_2/J_1 leq 0.5$) by means of density matrix renormalization group (DMRG). In the parameter regime $J_2/J_1 lesssim 0.27$, the system sustains a N{e}el antiferromagnetic phase. At the large $J_2$ side $J_2/J_1 gtrsim 0.32$, a stripe antiferromagnetic phase is found. Between the two magnetic ordered phases $0.27 lesssim J_2/J_1 lesssim 0.32$, we find a textit{non-magnetic} intermediate region with a plaquette valence-bond order. Although our calculations are limited within $6$ unit-cell width on cylinder, we present evidence that this plaquette state could be a strong candidate for this non-magnetic region in the thermodynamic limit. We also briefly discuss the nature of the quantum phase transitions in the system. We gain further insight of the non-magnetic phases in the spin-$1$ system by comparing its phase diagram with the spin-$1/2$ system.