Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Sign-Problem-Free Monte Carlo Simulation of Certain Frustrated Quantum Magnets

Fabien Alet, Kedar Damle, and Sumiran Pujari
Phys. Rev. Lett. 117, 197203 – Published 4 November 2016
PDFHTMLExport Citation

Abstract

We introduce a quantum Monte Carlo (QMC) method for efficient sign-problem-free simulations of a broad class of frustrated S=1/2 antiferromagnets using the basis of spin eigenstates of clusters to avoid the severe sign problem faced by other QMC methods. We demonstrate the utility of the method in several cases with competing exchange interactions and flag important limitations as well as possible extensions of the method.

  • Figure
  • Figure
  • Figure
  • Figure
  • Received 21 January 2016

DOI:https://doi.org/10.1103/PhysRevLett.117.197203

© 2016 American Physical Society

Physics Subject Headings (PhySH)

Condensed Matter, Materials & Applied PhysicsStatistical Physics & Thermodynamics

Authors & Affiliations

Fabien Alet1, Kedar Damle2, and Sumiran Pujari3

  • 1Laboratoire de Physique Théorique, IRSAMC, Université de Toulouse, CNRS, 31062 Toulouse, France
  • 2Department of Theoretical Physics, Tata Institute of Fundamental Research, Mumbai 400 005, India
  • 3Department of Physics and Astronomy, University of Kentucky, Lexington, Kentucky 40506-0055, USA

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 117, Iss. 19 — 4 November 2016

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×

Images

  • Figure 1
    Figure 1

    Vertices that appear in the SSE operator string for Hbilayer, with corresponding weights in the canonical cluster basis. All other valid vertices are obtained by symmetry operations that exchange left and right, or upper and lower, legs (keeping the weight fixed). The constant C in the function f(lA,lB,mA,mB)=CJzmAmBζ(ΔzΔ)(mA2+mB2)ζΔ[lA(lA+1)+lB(lB+1)] is chosen to ensure that f0. Bottom right: Lattice structure and exchange couplings of Hbilayer. Vertices and lattice structure for Hmixed are detailed in Supplemental Material [56].

    Reuse & Permissions
  • Figure 2
    Figure 2

    Temperature (T) dependence of the susceptibility and specific heat of Hbilayer with Dz=D=1, J=K=1, Jz=1+Kz, and Kz=1Kz. Symbols display data for a sample with L=64 unit cells, plotted for a variety of values of Kz. The inset shows the perfect agreement between QMC data (symbols) and exact diagonalization results (lines) for L=6.

    Reuse & Permissions
  • Figure 3
    Figure 3

    QMC results (symbols) for Hbilayer in a field on the square lattice, with Jz=Kz=1, Dz=D=5, J=1+K, K=1K, and h=7. Left panel: Specific heat Cv for K=0 (inset zooms into the critical range). Right panels: Binder cumulant U=ms4/ms22 of the staggered magnetization ms=r()r(SIrz+SIIrz). The critical temperature Tc, estimated by the crossing point of U, decreases with K. U at the estimated Tc tends to the 2d Ising critical value U*=1.16793 [62] at large L for all K displayed. The solid lines in the K>0 panels are guides to the eye. At K=0, they denote results for the 2d classical Ising model at TIsing=4T.

    Reuse & Permissions
  • Figure 4
    Figure 4

    This operator string for a single plaquette of the Bravais lattice of the bilayer system illustrates the origin of the sign problem faced when simulating the general bilayer Hamiltonian Hbilayer: Its weight is negative, independent of the signs of the nonzero couplings J, Kz, and K.

    Reuse & Permissions
×

Sign up to receive regular email alerts from Physical Review Letters

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×