Abstract
We describe a stochastic series expansion quantum Monte Carlo method for a two-dimensional S = 1/2 XY model (or, equivalently, hard-core bosons at half filling) which in addition to the standard pair interaction J includes a four-site term K that flips spins on a square plaquette. The model has three ordered ground state phases; for K/J approximately < or = 8 it has long-range xy spin order (superfluid bosons), for K/J approximately > or = 15 it has staggered spin order in the z direction (charge-density wave), and between these phases it is in a state with columnar order in the bond and plaquette energy densities. We discuss an implementation of directed-loop updates for the SSE simulations of this model and also introduce a "multibranch" cluster update which significantly reduces the autocorrelation times for large K/J. In addition to the pure J-K model, which in the z basis has only off-diagonal terms, we also discuss modifications of the algorithm needed when various diagonal interactions are included.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
