No Arabic abstract
Puzzled by the indication of a new critical theory for the spin-1/2 Heisenberg model with a spatially staggered anisotropy on the square lattice as suggested in cite{Wenzel08}, we study a similar anisotropic spin-1/2 Heisenberg model on the honeycomb lattice. The critical point where the phase transition occurs due to the dimerization as well as the critical exponent $ u$ are analyzed in great detail. Remarkly, using most of the available data points in conjunction with the expected finite-size scaling ansatz with a sub-leading correction indeed leads to a consistent $ u = 0.691(2)$ with that calculated in cite{Wenzel08}. However by using the data with large number of spins $N$, we obtain $ u = 0.707(6)$ which agrees with the most accurate Monte Carlo O(3) value $ u = 0.7112(5)$ as well.
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.
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 investigate the spin-1/2 Heisenberg model on a rectangular lattice, using the Gutzwiller projected variational wave function known as the staggered flux state. Using Monte Carlo techniques, the variational parameters and static spin-structure factor for different coupling anisotropies $gamma=J_y/J_x$ are calculated. We observe a gradual evolution of the ground state energy towards a value which is very close to the 1D estimate provided by the Bethe ansatz and a good agreement between the finite size scaling of the energies. The spin-spin correlation functions exhibit a power-law decay with varying exponents for different anisotropies. Though the lack of Neel order makes the staggered flux state energetically unfavorable in the symmetric case $gamma=1$, it appears to capture the essence of the system close to 1D. Hence we believe that the staggered flux state provides an interesting starting point to explore the crossover from quantum disordered chains to the Neel ordered 2D square lattices.
We present a spin-rotation-invariant Green-function theory for the dynamic spin susceptibility in the spin-1/2 antiferromagnetic Heisenberg model on a stacked honeycomb lattice. Employing a generalized mean-field approximation for arbitrary temperatures, the thermodynamic quantities (two-spin correlation functions, internal energy, magnetic susceptibility, staggered magnetization, Neel temperature, correlation length) and the spin-excitation spectrum are calculated by solving a coupled system of self-consistency equations for the correlation functions. The temperature dependence of the magnetic (uniform static) susceptibility is ascribed to antiferromagnetic short-range order. The N{e}el temperature is calculated for arbitrary interlayer couplings. Our results are in a good agreement with numerical computations for finite clusters and with available experimental data on the beta-Cu2V2O2 compound.
We study the field dependence of the antiferromagnetic spin-1/2 Heisenberg model on the square lattice by means of exact diagonalizations. In a first part, we calculate the spin-wave velocity, the spin-stiffness, and the magnetic susceptibility and thus determine the microscopic parameters of the low-energy long-wavelength description. In a second part, we present a comprehensive study of dynamical spin correlation functions for magnetic fields ranging from zero up to saturation. We find that at low fields, magnons are well defined in the whole Brillouin zone, but the dispersion is substantially modified by quantum fluctuations compared to the classical spectrum. At higher fields, decay channels open and magnons become unstable with respect to multi-magnon scattering. Our results directly apply to inelastic neutron scattering experiments.