No Arabic abstract
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.
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.
We use the state-of-the-art tensor network state method, specifically, the finite projected entangled pair state (PEPS) algorithm, to simulate the global phase diagram of spin-$1/2$ $J_1$-$J_2$ Heisenberg model on square lattices up to $24times 24$. We provide very solid evidences to show that the nature of the intermediate nonmagnetic phase is a gapless quantum spin liquid (QSL), whose spin-spin and dimer-dimer correlations both decay with a power law behavior. There also exists a valence-bond solid (VBS) phase in a very narrow region $0.56lesssim J_2/J_1leq0.61$ before the system enters the well known collinear antiferromagnetic phase. We stress that our work gives rise to the first solid PEPS results beyond the well established density matrix renormalization group (DMRG) through one-to-one direct benchmark for small system sizes. Thus our numerical evidences explicitly demonstrate the huge power of PEPS for solving long-standing 2D quantum many-body problems. The physical nature of the discovered gapless QSL and potential experimental implications are also addressed.
We have explored the magnetic excitation spectrum of the S=1/2 square lattice Heisenberg antiferromagnet, K2V3O8 using both triple-axis and time-of-flight inelastic neutron scattering. The long-wavelength spin waves are consistent with the previously determined Hamiltonian for this material. A small energy gap of 72+/-9 micro-eV is observed at the antiferromagnetic zone center and the near-neighbor exchange constant is determined to be 1.08+/-0.03 meV. A finite ferromagnetic interplanar coupling is observed along the crystallographic c-axis with a magnitude of Jc=-0.0036+/-0.006 meV. However, upon approaching the zone boundary, the observed excitation spectrum deviates significantly from the expectation of linear spin wave theory resulting in split modes at the (pi/2,pi/2) zone boundary point. The effects of magnon-phonon interaction, orbital degrees of freedom, multimagnon scattering, and dilution/site randomness are considered in the context of the mode splitting. Unfortunately, no fully satisfactory explanation of this phenomenon is found and further theoretical and experimental work is needed.
We study the quantum phase diagram and excitation spectrum of the frustrated $J_1$-$J_2$ spin-1/2 Heisenberg Hamiltonian. A hierarchical mean-field approach, at the heart of which lies the idea of identifying {it relevant} degrees of freedom, is developed. Thus, by performing educated, manifestly symmetry preserving mean-field approximations, we unveil fundamental properties of the system. We then compare various coverings of the square lattice with plaquettes, dimers and other degrees of freedom, and show that only the {it symmetric plaquette} covering, which reproduces the original Bravais lattice, leads to the known phase diagram. The intermediate quantum paramagnetic phase is shown to be a (singlet) {it plaquette crystal}, connected with the neighboring Neel phase by a continuous phase transition. We also introduce fluctuations around the hierarchical mean-field solutions, and demonstrate that in the paramagnetic phase the ground and first excited states are separated by a finite gap, which closes in the Neel and columnar phases. Our results suggest that the quantum phase transition between Neel and paramagnetic phases can be properly described within the Ginzburg-Landau-Wilson paradigm.
The properties of ground state of spin-$frac{1}{2}$ kagome antiferromagnetic Heisenberg (KAFH) model have attracted considerable interest in the past few decades, and recent numerical simulations reported a spin liquid phase. The nature of the spin liquid phase remains unclear. For instance, the interplay between symmetries and $Z_2$ topological order leads to different types of $Z_2$ spin liquid phases. In this paper, we develop a numerical simulation method based on symmetric projected entangled-pair states (PEPS), which is generally applicable to strongly correlated model systems in two spatial dimensions. We then apply this method to study the nature of the ground state of the KAFH model. Our results are consistent with that the ground state is a $U(1)$ Dirac spin liquid rather than a $Z_2$ spin liquid.