Phase diagram of an extended classical dimer model

We present an extensive numerical study of the critical behavior of dimer models in three dimensions, focusing on the phase transition between Coulomb and crystalline columnar phases. The case of attractive interactions between parallel dimers on a plaquette was shown to undergo a continuous phase transition with critical exponents close to those of the O(N) tricritical universality class, a situation which is not easily captured by conventional field theories. That the dimer model is exactly fine-tuned to a highly symmetric point is a non trivial statement which needs careful numerical investigation. In this paper, we perform an extensive Monte Carlo study of a generalized dimer model with plaquette and cubic interactions and determine its extended phase diagram. We find that when both interactions favor alignment of the dimers, the phase transition is first order, in almost all cases. On the opposite, when interactions compete, the transition becomes continuous, with a critical exponent \eta ~ 0.2. The existence of a tricritical point between the two regimes is confirmed by simulations on very large size systems and a flowgram method. In addition, we find a highly-degenerate crystalline phase at very low temperature in the frustrated regime which is separated from the columnar phase by a first order transition.Comment: 12 pages, 13 figure

Valence Bond Entanglement Entropy

We introduce for SU(2) quantum spin systems the Valence Bond Entanglement Entropy as a counting of valence bond spin singlets shared by two subsystems. For a large class of antiferromagnetic systems, it can be calculated in all dimensions with Quantum Monte Carlo simulations in the valence bond basis. We show numerically that this quantity displays all features of the von Neumann entanglement entropy for several one-dimensional systems. For two-dimensional Heisenberg models, we find a strict area law for a Valence Bond Solid state and multiplicative logarithmic corrections for the Neel phase.Comment: 4 pages, 3 figures, v2: small corrections, published versio

Phase Diagram of Interacting Bosons on the Honeycomb Lattice

We study the ground state properties of repulsively interacting bosons on the honeycomb lattice using large-scale quantum Monte Carlo simulations. In the hard-core limit the half-filled system develops long ranged diagonal order for sufficiently strong nearest-neighbor repulsion. This staggered solid melts at a first order quantum phase transition into the superfluid phase, without the presence of any intermediate supersolid phase. Within the superfluid phase, both the superfluid density and the compressibility exhibit local minima near particle- (hole-) density one quarter, while the density and the condensate fraction show inflection points in this region. Relaxing the hard-core constraint, supersolid phases emerge for soft-core bosons. The suppression of the superfluid density is found to persist for sufficiently large, finite on-site repulsion.Comment: 4 pages with 5 figure

Universal Reduction of Effective Coordination Number in the Quasi-One-Dimensional Ising Model

Critical temperature of quasi-one-dimensional general-spin Ising ferromagnets is investigated by means of the cluster Monte Carlo method performed on infinite-length strips, L times infty or L times L times infty. We find that in the weak interchain coupling regime the critical temperature as a function of the interchain coupling is well-described by a chain mean-field formula with a reduced effective coordination number, as the quantum Heisenberg antiferromagnets recently reported by Yasuda et al. [Phys. Rev. Lett. 94, 217201 (2005)]. It is also confirmed that the effective coordination number is independent of the spin size. We show that in the weak interchain coupling limit the effective coordination number is, irrespective of the spin size, rigorously given by the quantum critical point of a spin-1/2 transverse-field Ising model.Comment: 12 pages, 6 figures, minor modifications, final version published in Phys. Rev.

Hybridization expansion impurity solver: General formulation and application to Kondo lattice and two-orbital models

A recently developed continuous time solver based on an expansion in hybridization about an exactly solved local limit is reformulated in a manner appropriate for general classes of quantum impurity models including spin exchange and pair hopping terms. The utility of the approach is demonstrated via applications to the dynamical mean field theory of the Kondo lattice and two-orbital models. The algorithm can handle low temperatures and strong couplings without encountering a sign problem.Comment: Published versio

Spin gap and string order parameter in the ferromagnetic Spiral Staircase Heisenberg Ladder: a quantum Monte Carlo study

We consider a spin-1/2 ladder with a ferromagnetic rung coupling J_\perp and inequivalent chains. This model is obtained by a twist (\theta) deformation of the ladder and interpolates between the isotropic ladder (\theta=0) and the SU(2) ferromagnetic Kondo necklace model (\theta=\pi). We show that the ground state in the (\theta,J_\perp) plane has a finite string order parameter characterising the Haldane phase. Twisting the chain introduces a new energy scale, which we interpret in terms of a Suhl-Nakamura interaction. As a consequence we observe a crossover in the scaling of the spin gap at weak coupling from \Delta/J_\| \propto J_\perp/J_\| for \theta < \theta_c \simeq 8\pi/9 to \Delta/J_\| \propto (J_\perp/J_\|)^2 for \theta > \theta_c. Those results are obtained on the basis of large scale Quantum Monte Carlo calculations.Comment: 4 page

Magnetization plateaus of an easy-axis Kagom\'e antiferromagnet with extended interactions

We investigate the properties in finite magnetic field of an extended anisotropic XXZ spin-1/2 model on the Kagome lattice, originally introduced by Balents, Fisher, and Girvin [Phys. Rev. B, 65, 224412 (2002)]. The magnetization curve displays plateaus at magnetization m=1/6 and 1/3 when the anisotropy is large. Using low-energy effective constrained models (quantum loop and quantum dimer models), we discuss the nature of the plateau phases, found to be crystals that break discrete rotation and/or translation symmetries. Large-scale quantum Monte-Carlo simulations were carried out in particular for the m=1/6 plateau. We first map out the phase diagram of the effective quantum loop model with an additional loop-loop interaction to find stripe order around the point relevant for the original model as well as a topological Z2 spin liquid. The existence of a stripe crystalline phase is further evidenced by measuring both standard structure factor and entanglement entropy of the original microscopic model.Comment: 14 pages, 14 figure

Neel order in square and triangular lattice Heisenberg models

Using examples of the square- and triangular-lattice Heisenberg models we demonstrate that the density matrix renormalization group method (DMRG) can be effectively used to study magnetic ordering in two-dimensional lattice spin models. We show that local quantities in DMRG calculations, such as the on-site magnetization M, should be extrapolated with the truncation error, not with its square root, as previously assumed. We also introduce convenient sequences of clusters, using cylindrical boundary conditions and pinning magnetic fields, which provide for rapidly converging finite-size scaling. This scaling behavior on our clusters is clarified using finite-size analysis of the effective sigma-model and finite-size spin-wave theory. The resulting greatly improved extrapolations allow us to determine the thermodynamic limit of M for the square lattice with an error comparable to quantum Monte Carlo. For the triangular lattice, we verify the existence of three-sublattice magnetic order, and estimate the order parameter to be M = 0.205(15).Comment: 4 pages, 5 figures, typo fixed, reference adde

Non-local updates for quantum Monte Carlo simulations

We review the development of update schemes for quantum lattice models simulated using world line quantum Monte Carlo algorithms. Starting from the Suzuki-Trotter mapping we discuss limitations of local update algorithms and highlight the main developments beyond Metropolis-style local updates: the development of cluster algorithms, their generalization to continuous time, the worm and directed-loop algorithms and finally a generalization of the flat histogram method of Wang and Landau to quantum systems.Comment: 14 pages, article for the proceedings of the "The Monte Carlo Method in the Physical Sciences: Celebrating the 50th Anniversary of the Metropolis Algorithm", Los Alamos, June 9-11, 200

Unconventional continuous phase transition in a three dimensional dimer model

Phase transitions occupy a central role in physics, due both to their experimental ubiquity and their fundamental conceptual importance. The explanation of universality at phase transitions was the great success of the theory formulated by Ginzburg and Landau, and extended through the renormalization group by Wilson. However, recent theoretical suggestions have challenged this point of view in certain situations. In this Letter we report the first large-scale simulations of a three-dimensional model proposed to be a candidate for requiring a description beyond the Landau-Ginzburg-Wilson framework: we study the phase transition from the dimer crystal to the Coulomb phase in the cubic dimer model. Our numerical results strongly indicate that the transition is continuous and are compatible with a tricritical universality class, at variance with previous proposals.Comment: 4 pages, 3 figures; v2: minor changes, published versio