Susceptibility of the 2D spin-1 / 2 Heisenberg antiferromagnet with an impurity

We use a quantum Monte Carlo method (stochastic series expansion) to study the effects of a magnetic or nonmagnetic impurity on the magnetic susceptibility of the two-dimensional Heisenberg antiferromagnet. At low temperatures, we find a log-divergent contribution to the transverse susceptibility. We also introduce an effective few-spin model that can quantitatively capture the differences between magnetic and nonmagnetic impurities at high and intermediate temperatures.     

Valence bond solid phases in a cubic antiferromagnet

We report on a valence bond projector Monte Carlo simulation of the cubic lattice quantum Heisenberg model with additional higher-order exchange interactions in each unit cell. The model supports two different valence bond solid (VBS) ground states. In one of these states, the dimer pattern is a three-dimensional analogue of the columnar pattern familiar from two dimensions. In the other, the dimers are regularly arranged along the four main diagonals in 1/8 of the unit cells. The phases are separated from one another and from a Néel phase by strongly first-order boundaries. Our results strengthen the case for exotic transitions in two dimensions, where no discontinuities have been detected at the Heisenberg Néel-VBS transition driven by four-spin plaquette interactions.     

Impurity induced spin texture in quantum critical 2D antiferromagnets

We describe the uniform and staggered magnetization distributions around a vacancy in a quantum critical two-dimensional S=1/2 antiferromagnet. The distributions are delocalized across the entire sample with a universal functional form arising from an uncompensated Berry phase. The numerical results, obtained using quantum Monte Carlo simulations of the Heisenberg model on bilayer lattices with up to approximately equal to 10(5) spins, are in good agreement with the proposed scaling structure. We determine the exponent eta' = 0.40+/-0.02, which governs both the staggered and uniform magnetic structure away from the impurity and also controls the impurity spin dynamics.     

Method to characterize spinons as emergent elementary particles

We develop a technique to directly study spinons (emergent spin S=1/2 particles) in quantum spin models in any number of dimensions. The size of a spinon wave packet and of a bound pair (a triplon) are defined in terms of wave-function overlaps that can be evaluated by quantum Monte?Carlo simulations. We show that the same information is contained in the spin-spin correlation function as well. We illustrate the method in one dimension. We confirm that spinons are well-defined particles (have exponentially localized wave packet) in a valence-bond-solid state, are marginally defined (with power-law shaped wave packet) in the standard Heisenberg critical state, and are not well defined in an ordered Néel state (achieved in one dimension using long-range interactions).