BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Anders W. Sandvik

Anders W. Sandvik appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.stat-mechhep-latquant-phcond-mat.dis-nncond-matcond-mat.mtrl-scicond-mat.quant-gasmath-phmath.MPphysics.comp-ph

Catalog footprint

What is connected

60works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

60 published item(s)

preprint2023arXiv

Progress on stochastic analytic continuation of quantum Monte Carlo data

We report multipronged progress on the stochastic averaging approach to numerical analytic continuation of quantum Monte Carlo data. With the sampled spectrum parametrized with delta-functions in continuous frequency space, a calculation of the configurational entropy lends support to a simple goodness-of-fit criterion for the optimal sampling temperature. To further investigate entropic effects, we compare spectra sampled in continuous frequency with results of amplitudes sampled on a fixed frequency grid. We demonstrate equivalences between sampling and optimizing spectral functions with the maximum-entropy approach with different forms of the entropy. These insights revise prevailing notions of the maximum-entropy method and its relationship to stochastic analytic continuation. We further explore various adjustable (optimized) constraints that allow sharp spectral features to be resolved, in particular at the lower frequency edge. The constraints, e.g., the location of the edge or the spectral weight of a quasi-particle peak, are optimized using a statistical criterion. We show that this method can correctly reproduce both narrow and broad quasi-particle peaks. We next introduce a parametrization for more intricate spectral functions with sharp edges, e.g., power-law singularities. Tests with synthetic data as well as with real simulation data for the spin-1/2 Heisenberg chain demonstrate that constrained sampling methods can reproduce spectral functions with sharp edge features at unprecedented fidelity. We present new results for S=1/2 Heisenberg 2-leg and 3-leg ladders to illustrate the ability of the methods to resolve spectral features arising from both elementary and composite excitations. Finally, we also propose how the methods developed here could be used as "pre processors" for analytic continuation by machine learning.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fractional and composite excitations of antiferromagnetic quantum spin trimer chains

Using quantum Monte Carlo, exact diagonalization and perturbation theory, we study the spectrum of the $S=1/2$ antiferromagnetic Heisenberg trimer chain by varying the ratio $g=J_2/J_1$ of the intertrimer and intratrimer coupling strengths. The doublet ground states of trimers form effective interacting $S=1/2$ degrees of freedom described by a Heisenberg chain. Therefore, the conventional two-spinon continuum of width $\propto J_1$ when $g=1$ evolves into to a similar continuum of width $\propto J_2$ when $g\to 0$. The intermediate-energy and high-energy modes are termed \emph{doublons} and \emph{quartons} which fractionalize with increasing $g$ to form the conventional spinon continuum. In particular, at $g \approx 0.716$, the gap between the low-energy spinon branch and the high-energy band with mixed doublons, quartons, and spinons closes. These features should be observable in inelastic neutron scattering experiments if a quasi-one-dimensional quantum magnet with the linear trimer structure and $J_2<J_1$ can be identified. Our results may open a window for exploring the high-energy fractional excitations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Criticality and Spin Liquid Phase in the Shastry-Sutherland model

Using the density-matrix renormalization group method for the ground state and excitations of the Shastry- Sutherland spin model, we demonstrate the existence of a narrow quantum spin liquid phase between the previously known plaquette-singlet and antiferromagnetic states. Our conclusions are based on finite-size scaling of excited level crossings and order parameters. Together with previous results on candidate models for deconfined quantum criticality and spin liquid phases, our results point to a unified quantum phase diagram where the deconfined quantum-critical point separates a line of first-order transitions and a gapless spin liquid phase. The frustrated Shastry-Sutherland model is close to the critical point but slightly inside the spin liquid phase, while previously studied unfrustrated models cross the first-order line. We also argue that recent heat capacity measurements in SrCu2(BO3)2 show evidence of the proposed spin liquid at pressures between 2.6 and 3 GPa.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum spin liquid phase in the Shastry-Sutherland model detected by an improved level spectroscopic method

We study the spin-$1/2$ two-dimensional Shastry-Sutherland spin model by exact diagonalization of clusters with periodic boundary conditions. We develop an improved level spectroscopic technique using energy gaps between states with different quantum numbers. The crossing points of some of the relative (composite) gaps have much weaker finite-size drifts than the normally used gaps defined only with respect to the ground state, thus allowing precise determination of quantum critical points even with small clusters. Our results support the picture of a spin liquid phase intervening between the well known plaquette-singlet and antiferromagnetic ground states, with phase boundaries in almost perfect agreement with a recent density matrix renormalization group study, where much larger cylindrical lattices were used [J. Yang et al., Phys. Rev. B {\bf 105}, L060409 (2022)]. The method of using composite low-energy gaps to reduce scaling corrections has potentially broad applications in numerical studies of quantum critical phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extreme Suppression of Antiferromagnetic Order and Critical Scaling in a Two-Dimensional Random Quantum Magnet

Sr$_2$CuTeO$_6$ is a square-lattice Néel antiferromagnet with superexchange between first-neighbor $S=1/2$ Cu spins mediated by plaquette centered Te ions. Substituting Te by W, the affected impurity plaquettes have predominantly second-neighbor interactions, thus causing local magnetic frustration. Here we report a study of Sr$_2$CuTe$_{1-x}$W$_x$O$_6$ using neutron diffraction and $μ$SR techniques, showing that the Néel order vanishes already at $x = 0.025 \pm 0.005$. We explain this extreme order suppression using a two-dimensional Heisenberg spin model, demonstrating that a W-type impurity induces a deformation of the order parameter that decays with distance as $1/r^2$ at temperature $T=0$. The associated logarithmic singularity leads to loss of order for any $x>0$. Order for small $x>0$ and $T>0$ is induced by weak interplane couplings. In the nonmagnetic phase of Sr$_2$CuTe$_{1-x}$W$_x$O$_6$, the $μ$SR relaxation rate exhibits quantum critical scaling with a large dynamic exponent, $z \approx 3$, consistent with a random-singlet state.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Consistent scaling exponents at the deconfined quantum-critical point

We report a quantum Monte Carlo study of the phase transition between antiferromagnetic and valence-bond solid ground states in the square-lattice $S=1/2$ $J$-$Q$ model. The critical correlation function of the $Q$ terms gives a scaling dimension corresponding to the value $ν= 0.455 \pm 0.002$ of the correlation-length exponent. This value agrees with previous (less precise) results from conventional methods, e.g., finite-size scaling of the near-critical order parameters. We also study the $Q$-derivatives of the Binder cumulants of the order parameters for $L^2$ lattices with $L$ up to $448$. The slope grows as $L^{1/ν}$ with a value of $ν$ consistent with the scaling dimension of the $Q$ term. There are no indications of runaway flow to a first-order phase transition. The mutually consistent estimates of $ν$ provide compelling support for a continuous deconfined quantum-critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deconfined quantum criticality in spin-1/2 chains with long-range interactions

We study spin-$1/2$ chains with long-range power-law decaying unfrustrated (bipartite) Heisenberg exchange $J_r \propto r^{-α}$ and multi-spin interactions $Q$ favoring a valence-bond solid (VBS) ground state. Employing quantum Monte Carlo techniques and Lanczos diagonalization, we analyze order parameters and excited-state level crossings to characterize quantum states and phase transitions in the $(α,Q)$ plane. For weak $Q$ and sufficiently slowly decaying Heisenberg interactions (small $α$), the system has a long-range-ordered antiferromagnetic (AFM) ground state, and upon increasing $α$ there is a continuous transition into a quasi long-range ordered (QLRO) critical state of the type in the standard Heisenberg chain. For rapidly decaying long-range interactions, there is transition between QLRO and VBS ground states of the same kind as in the frustrated $J_1$-$J_2$ Heisenberg chain. Our most important finding is a direct continuous quantum phase transition between the AFM and VBS states - a close analogy to the 2D deconfined quantum-critical point. In previous 1D analogies the ordered phases both have gapped fractional excitations, and the critical point is a conventional Luttinger Liquid. In our model the excitations fractionalize upon transitioning from the AFM state, changing from spin waves to deconfined spinons. We extract critical exponents at the AFM-VBS transition and use order-parameter distributions to study emergent symmetries. We find emergent O($4$) symmetry of the O($3$) AFM and scalar VBS order parameters. Thus, the order parameter fluctuations exhibit the covariance of a uniaxially deformed O($4$) sphere (an "elliptical" symmetry). This unusual quantum phase transition does not yet have any known field theory description, and our detailed results can serve to guide its construction. We discuss possible experimental realizations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Existence of a Spectral Gap in the Affleck-Kennedy-Lieb-Tasaki Model on the Hexagonal Lattice

The $S=1$ Affleck-Kennedy-Lieb-Tasaki (AKLT) quantum spin chain was the first rigorous example of an isotropic spin system in the Haldane phase. The conjecture that the $S=3/2$ AKLT model on the hexagonal lattice is also in a gapped phase has remained open, despite being a fundamental problem of ongoing relevance to condensed-matter physics and quantum information theory. Here we confirm this conjecture by demonstrating the size-independent lower bound $Δ>0.006$ on the spectral gap of the hexagonal model with periodic boundary conditions in the thermodynamic limit. Our approach consists of two steps combining mathematical physics and high-precision computational physics. We first prove a mathematical finite-size criterion which gives an analytical, size-independent bound on the spectral gap if the gap of a particular cut-out subsystem of 36 spins exceeds a certain threshold value. Then we verify the finite-size criterion numerically by performing state-of-the-art DMRG calculations on the subsystem.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hilbert Space Fragmentation and Ashkin-Teller Criticality in Fluctuation Coupled Ising Models

We discuss the effects of exponential fragmentation of the Hilbert space on phase transitions in the context of coupled ferromagnetic Ising models in arbitrary dimension with special emphasis on the one dimensional case. We show that the dynamics generated by quantum fluctuations is bounded within spatial partitions of the system and weak mixing of these partitions caused by global transverse fields leads to a zero temperature phase with ordering in the local product of both Ising copies but no long range order in either species. This leads to a natural connection with the Ashkin-Teller universality class for general lattices. We confirm this for the periodic chain using quantum Monte Carlo simulations. We also point out that our treatment provides an explanation for pseudo-first order behavior seen in the Binder cumulants of the classical frustrated $J_1-J_2$ Ising model and the $q=4$ Potts model in 2D.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Monte Carlo Renormalization Flows in the Space of Relevant and Irrelevant Operators: Application to Three-Dimensional Clock Models

We present a way to visualize and quantify renormalization group flows in a space of observables computed using Monte Carlo simulations. We apply the method to classical three-dimensional clock models, i.e., the planar (XY) spin model perturbed by a $Z_q$ symmetric anisotropy field. The method performs significantly better than standard techniques for determining the scaling dimension $y_q$ of the $Z_q$ field at the critical point if it is irrelevant ($q\ge4$). Furthermore, we analyze all stages of the complex renormalization flow, including the cross-over from the U(1) Nambu-Goldstone fixed point to the ultimate $Z_q$ symmetry-breaking fixed point due to the relevance of the $Z_q$ field inside the ordered phase. We expect our method to be particularly useful in the context of quantum-critical points with inherent dangerously irrelevant operators that cannot be tuned away microscopically but whose renormalization flows can be analyzed exactly as we do here for the clock models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum phases of SrCu2(BO3)2 from high-pressure thermodynamics

We report heat capacity measurements of SrCu$_2$(BO$_3$)$_2$ under high pressure along with simulations of relevant quantum spin models and map out the $(P,T)$ phase diagram of the material. We find a first-order quantum phase transition between the low-pressure quantum dimer paramagnet and a phase with signatures of a plaquette-singlet state below T = $2$ K. At higher pressures, we observe a transition into a previously unknown antiferromagnetic state below $4$ K. Our findings can be explained within the two-dimensional Shastry-Sutherland quantum spin model supplemented by weak inter-layer couplings. The possibility to tune SrCu$_2$(BO$_3$)$_2$ between the plaquette-singlet and antiferromagnetic states opens opportunities for experimental tests of quantum field theories and lattice models involving fractionalized excitations, emergent symmetries, and gauge fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum-critical scaling properties of the two-dimensional random-singlet state

We use QMC simulations to study effects of disorder on the $S=1/2$ Heisenberg model with exchange constant $J$ on the square lattice supplemented by multispin interactions $Q$. It was found recently [L. Lu et al., Phys. Rev. X 8, 041040 (2018)] that the ground state of this $J$-$Q$ model with random couplings undergoes a quantum phase transition from the Néel state into a randomness-induced spin-liquid-like state that is a close analogue to the well known random-singlet (RS) state of the random Heisenberg chain. The 2D RS state arises from spinons localized at topological defects. The interacting spinons form a critical state with mean spin-spin correlations decaying with distance $r$ as $r^{-2}$, as in the 1D RS state. The dynamic exponent $z \ge 2$, varying continuously with the model parameters. We here further investigate the properties of the RS state in the $J$-$Q$ model with random $Q$ couplings. We study the temperature dependence of the specific heat and various susceptibilities for large enough systems to reach the thermodynamic limit and also analyze the size dependence of the critical magnetic order parameter and its susceptibility in the ground state. For all these quantities, we find consistency with conventional quantum-critical scaling when the condition implied by the $r^{-2}$ form of the spin correlations is imposed. All quantities can be explained by the same value of the dynamic exponent $z$ at fixed model parameters. We argue that the RS state identified in the $J$-$Q$ model corresponds to a generic renormalization group fixed point that can be reached in many quantum magnets with random couplings, and may already have been observed experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Scaling and diabatic effects in quantum annealing with a D-Wave device

We discuss quantum annealing of the two-dimensional transverse-field Ising model on a D-Wave device, encoded on $L \times L$ lattices with $L \le 32$. Analyzing the residual energy and deviation from maximal magnetization in the final classical state, we find an optimal $L$ dependent annealing rate $v$ for which the two quantities are minimized. The results are well described by a phenomenological model with two powers of $v$ and $L$-dependent prefactors to describe the competing effects of reduced quantum fluctuations (for which we see evidence of the Kibble-Zurek mechanism) and increasing noise impact when $v$ is lowered. The same scaling form also describes results of numerical solutions of a transverse-field Ising model with the spins coupled to noise sources. We explain why the optimal annealing time is much longer than the coherence time of the individual qubits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tunable deconfined quantum criticality and interplay of different valence-bond solid phases

We use quantum Monte Carlo simulations to study a quantum $S=1/2$ spin model with competing multi-spin interactions. We find a quantum phase transition between a columnar valence-bond solid (cVBS) and a Néel antiferromagnet (AFM), as in the scenario of deconfined quantum-critical points, as well as a transition between the AFM and a staggered valence-bond solid (sVBS). By continuously varying a parameter, the sVBS--AFM and AFM--cVBS boundaries merge into a direct sVBS--cVBS transition. Unlike previous models with putative deconfined AFM--cVBS transitions, e.g., the standard $J$-$Q$ model, in our extended $J$-$Q$ model with competing cVBS and sVBS inducing terms the transition can be tuned from continuous to first-order. We find the expected emergent U(1) symmetry of the microscopically $Z_4$ symmetric cVBS order parameter when the transition is continuous. In contrast, when the transition changes to first-order the clock-like $Z_4$ fluctuations are absent and there is no emergent higher symmetry. We argue that the confined spinons in the sVBS phase are fracton-like. We also present results for an SU(3) symmetric model with a similar phase diagram. The new family of models can serve as a useful tool for further investigating open questions related to deconfined quantum criticality and its associated emergent symmetries.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Comment on "Gapless spin liquid ground state of the spin-$\frac{1}{2}$ $J_1$-$J_2$ Heisenberg model on square lattices"

Liu et al. [Phys.Rev.B 98, 241109 (2018)] used Monte Carlo sampling of the physical degrees of freedom of a Projected Entangled Pair State (PEPS) type wave function for the $S=1/2$ frustrated $J_1$-$J_2$ Heisenberg model on the square lattice and found a non-magnetic state argued to be a gapless spin liquid when the coupling ratio $g=J_2/J_1$ is in the range $g \in [0.42,0.6]$. Here we show that their definition of the order parameter for another candidate ground state within this coupling window---a spontaneously dimerized state---is problematic. The order parameter as defined will not detect dimer order when lattice symmeties are broken due to open boundaries or asymmetries originating from the calculation itself. Thus, a dimerized phase for some range of $g$ cannot be excluded (and is likely based on several other recent works).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Level spectroscopy in a two-dimensional quantum magnet: Linearly dispersing spinons at the deconfined quantum critical point

We study the level structure of excitations at the "deconfined" critical point separating antiferromagnetic and valence-bond-solid phases in two-dimensional quantum spin systems using the $J$-$Q$ model as an example. Energy gaps in different spin ($S$) and momentum (${\bf k}$) sectors are extracted from imaginary-time correlation functions obtained in quantum Monte Carlo simulations. We find strong quantitative evidence for deconfined linearly dispersing spinons with gapless points at ${\bf k}=(0,0)$, $(π,0)$, $(0,π)$, and $(π,π)$, as inferred from two-spinon excitations ($S=0$ and $S=1$ states) around these points. We also observe a duality between singlet and triplet excitations at the critical point and inside the ordered phases, in support of an enhanced symmetry, possibly SO(5).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Properties of the random-singlet phase: from the disordered Heisenberg chain to an amorphous valence-bond solid

We use a strong-disorder renormalization group (SDRG) method and ground-state quantum Monte Carlo (QMC) simulations to study S=1/2 spin chains with random couplings, calculating disorder-averaged spin and dimer correlations. The QMC simulations demonstrate logarithmic corrections to the power-law decaying correlations obtained with the SDRG scheme. The same asymptotic forms apply both for systems with standard Heisenberg exchange and for certain multi-spin couplings leading to spontaneous dimerization in the clean system. We show that the logarithmic corrections arise in the valence-bond (singlet pair) basis from a contribution that can not be generated by the SDRG scheme. In the model with multi-spin couplings, where the clean system dimerizes spontaneously, random singlets form between spinons localized at domain walls in the presence of disorder. This amorphous valence-bond solid is asymptotically a random-singlet state and only differs from the random-exchange Heisenberg chain in its short-distance properties.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum criticality with two length scales

The theory of deconfined quantum critical points describes phase transitions at temperature T = 0 outside the standard paradigm, predicting continuous transformations between certain ordered states where conventional theory requires discontinuities. Numerous computer simulations have offered no proof of such transitions, however, instead finding deviations from expected scaling relations that were neither predicted by the DQC theory nor conform to standard scenarios. Here we show that this enigma can be resolved by introducing a critical scaling form with two divergent length scales. Simulations of a quantum magnet with antiferromagnetic and dimerized ground states confirm the form, proving a continuous transition with deconfined excitations and also explaining anomalous scaling at T > 0. Our findings revise prevailing paradigms for quantum criticality, with potentially far-reaching implications for many strongly-correlated materials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum Monte Carlo studies of spinons in one-dimensional spin systems

Observing constituent particles with fractional quantum numbers in confined and deconfined states is an interesting and challenging problem in quantum many-body physics. Here we further explore a computational scheme [Y. Tang and A. W. Sandvik, Phys. Rev. Lett. {\bf 107}, 157201 (2011)] based on valence-bond quantum Monte Carlo simulations of quantum spin systems. Using several different one-dimensional models, we characterize $S=1/2$ spinon excitations using the spinon size and confinement length (the size of a bound state). The spinons have finite size in valence-bond-solid states, infinite size in the critical region, and become ill-defined in the Néel state. We also verify that pairs of spinons are deconfined in these uniform spin chains but become confined upon introducing a pattern of alternating coupling strengths (dimerization) or coupling two chains (forming a ladder). In the dimerized system an individual spinon can be small when the confinement length is large---this is the case when the imposed dimerization is weak but the ground state of the corresponding uniform chain is a spontaneously formed valence-bond-solid (where the spinons are deconfined). Based on our numerical results, we argue that the situation $λ\ll Λ$ is associated with weak repulsive short-range spinon-spinon interactions. In principle both the length-scales can be individually tuned from small to infinite (with $λ\le Λ$) by varying model parameters. In the ladder system the two lengths are always similar, and this is the case also in the dimerized systems when the corresponding uniform chain is in the critical phase. In these systems the effective spinon-spinon interactions are purely attractive and there is only a single large length scale close to criticality, which is reflected in the standard spin correlations as well as in the spinon characteristics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

1D valence bond solids in a magnetic field

A Valence bond solid (VBS) is a nonmagnetic, long-range ordered state of a quantum spin system where local spin singlets are formed in some regular pattern. We here study the competition between VBS order and a fully polarized ferromagnetic state as function of an external magnetic field in a one-dimensional extended Heisenberg model---the J-Q$_2$ model---using stochastic series expansion quantum Monte Carlo simulations with directed loop updates. We discuss the ground state phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anomalous quantum glass of bosons in a random potential in two dimensions

We present a quantum Monte Carlo study of the "quantum glass" phase of the 2D Bose-Hubbard model with random potentials at filling $ρ=1$. In the narrow region between the Mott and superfluid phases the compressibility has the form $κ\sim {\rm exp}(-b/T^α)+c$ with $α<1$ and $c$ vanishing or very small. Thus, at $T=0$ the system is either incompressible (a Mott glass) or nearly incompressible (a Mott-glass-like anomalous Bose glass). At stronger disorder, where a glass reappears from the superfluid, we find a conventional highly compressible Bose glass. On a path connecting these states, away from the superfluid at larger Hubbard repulsion, a change of the disorder strength by only $10\%$ changes the low-temperature compressibility by more than four orders of magnitude, lending support to two types of glass states separated by a phase transition or a sharp cross-over.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Emergent topological excitations in a two-dimensional quantum spin system

We study the mechanism of decay of a topological (winding-number) excitation due to finite-size effects in a two-dimensional valence-bond solid state, realized in an $S=1/2$ spin model ($J$-$Q$ model) and studied using projector Monte Carlo simulations in the valence bond basis. A topological excitation with winding number $|W|>0$ contains domain walls, which are unstable due to the emergence of long valence bonds in the wave function, unlike in effective descriptions with the quantum dimer model. We find that the life time of the winding number in imaginary time diverges as a power of the system length $L$. The energy can be computed within this time (i.e., it converges toward a "quasi-eigenvalue" before the winding number decays) and agrees for large $L$ with the domain-wall energy computed in an open lattice with boundary modifications enforcing a domain wall. Constructing a simplified two-state model and using the imaginary-time behavior from the simulations as input, we find that the real-time decay rate out of the initial winding sector is exponentially small in $L$. Thus, the winding number rapidly becomes a well-defined conserved quantum number for large systems, supporting the conclusions reached by computing the energy quasi-eigenvalues. Including Heisenberg exchange interactions which brings the system to a quantum-critical point separating the valence-bond solid from an antiferromagnetic ground state (the putative "deconfined" quantum-critical point), we can also converge the domain wall energy here and find that it decays as a power-law of the system size. Thus, the winding number is an emergent quantum number also at the critical point, with all winding number sectors becoming degenerate in the thermodynamic limit. This supports the description of the critical point in terms of a U(1) gauge-field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mott glass phase in a diluted bilayer Heisenberg quantum antiferromagnet

We use quantum Monte Carlo simulations to study a dimer-diluted $S=1/2$ Heisenberg model on a bilayer square lattice with intralayer interaction $J_{1}$ and interlayer interaction $J_{2}$. Below the classical percolation threshold $p_c$, the system has three phases reachable by tuning the interaction ratio $g=J_{2}/J_{1}$: a N$\acute{e}$el ordered phase, a gapless quantum glass phase, and a gapped quantum paramagnetic phase. We present the ground-state phase diagram in the plane of dilution $p$ and interaction ratio $g$. The quantum glass phase is certified to be of the gapless Mott glass type, having a uniform susceptibility vanishing at zero temperature $T$ and following a stretched exponential form at $T>0$; $χ_u \sim \exp(-b/T^α)$ with $α< 1$. At the phase transition point from N$\acute{e}$el ordered to Mott glass, we find that the critical exponents are different from those of the clean system described by the standard $O(3)$ universality class in 2+1 dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiplicative logarithmic corrections to quantum criticality in three-dimensional dimerized antiferromagnets

We investigate the quantum phase transition in an $S = 1/2$ dimerized Heisenberg antiferromagnet in three spatial dimensions. By performing large-scale quantum Monte Carlo simulations and detailed finite-size scaling analyses, we obtain high-precision results for the quantum critical properties at the transition from the magnetically disordered dimer-singlet phase to the antiferromagnetically ordered Néel phase. This transition breaks O($N$) symmetry with $N = 3$ in $D = 3 + 1$ dimensions. This is the upper critical dimension, where multiplicative logarithmic corrections to the leading mean-field critical properties are expected; we extract these corrections, establishing their precise forms for both the zero-temperature staggered magnetization, $m_s$, and the Néel temperature, $T_N$. We present a scaling Ansatz for $T_N$, including logarithmic corrections, which agrees with our data and indicates exact linearity with $m_s$, implying a complete decoupling of quantum and thermal fluctuation effects even arbitrarily close to the quantum critical point. We also demonstrate the predicted $N$-independent leading and subleading logarithmic corrections in the size-dependence of the staggered magnetic susceptibility. These logarithmic scaling forms have not previously been identified or verified by unbiased numerical methods and we discuss their relevance to experimental studies of dimerized quantum antiferromagnets such as TlCuCl$_3$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum versus classical annealing: insights from scaling theory and results for spin glasses on 3-regular graphs

We discuss an Ising spin glass where each $S=1/2$ spin is coupled antiferromagnetically to three other spins (3-regular graphs). Inducing quantum fluctuations by a time-dependent transverse field, we use out-of-equilibrium quantum Monte Carlo simulations to study dynamic scaling at the quantum glass transition. Comparing the dynamic exponent and other critical exponents with those of the classical (temperature-driven) transition, we conclude that quantum annealing is less efficient than classical simulated annealing in bringing the system into the glass phase. Quantum computing based on the quantum annealing paradigm is therefore inferior to classical simulated annealing for this class of problems. We also comment on previous simulations where a parameter is changed with the simulation time, which is very different from the true Hamiltonian dynamics simulated here.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Topological properties of a Valence-Bond-Solid

We present a projector quantum Monte Carlo study of the topological properties of the valence-bond-solid ground state in the $J$-$Q_3$ spin model on the square lattice. The winding number is a topological number counting the number of domain walls in the system and is a good quantum number in the thermodynamic limit. We study the finite-size behaviour and obtain the domain wall energy density for a topological nontrivial valence-bond-solid state.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Velocity of excitations in ordered, disordered and critical antiferromagnets

We test three different approaches, based on quantum Monte Carlo simulations, for computing the velocity $c$ of triplet excitations in antiferromagnets. We consider the standard $S=1/2$ one- and two-dimensional Heisenberg models, as well as a bilayer Heisenberg model at its critical point. Computing correlation functions in imaginary time and using their long-time behavior, we extract the lowest excitation energy versus momentum using improved fitting procedures and a generalized moment method. The velocity is then obtained from the dispersion relation. We also exploit winding numbers to define a cubic space-time geometry, where the velocity is obtained as the ratio of the spatial and temporal lengths of the system when all winding number fluctuations are equal. The two methods give consistent results for both ordered and critical systems, but the winding-number estimator is more precise. For the Heisenberg chain, we accurately reproduce the exactly known velocity. For the two-dimensional Heisenberg model, our results are consistent with other recent calculations, but with an improved statistical precision; $c=1.65847(4)$. We also use the hydrodynamic relation $c^2=ρ_s/χ_\perp(q\to 0)$ between $c$, the spin stiffness $ρ_s$, and the transversal susceptibility $χ_\perp$, using the smallest non-zero momentum $q=2π/L$. This method also is well controlled in two dimensions, but the cubic criterion for winding numbers delivers better numerical precision. In one dimension the hydrodynamic relation is affected by logarithmic corrections which make accurate extrapolations difficult. As an application of the winding-number method, for the quantum-critical bilayer model our high-precision determination of the velocity enables us to quantitatively test, at an unprecedented level, field-theoretic predictions for low-temperature scaling forms where $c$ enters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamic scaling at classical phase transitions approached through non-equilibrium quenching

We use Monte Carlo simulations to demonstrate generic scaling aspects of classical phase transitions approached through a quench (or annealing) protocol where the temperature changes as a function of time with velocity $v$. Using a generalized Kibble-Zurek ansatz, we demonstrate dynamic scaling for different types of stochastic dynamics (Metropolis, Swendsen-Wang, and Wolff) on Ising models in two and higher dimensions. We show that there are dual scaling functions governing the dynamic scaling, which together describe the scaling behavior in the entire velocity range $v \in [0, \infty)$. These functions have asymptotics corresponding to the adiabatic and diabatic limit, and close to these limits they are perturbative in $v$ and $1/v$, respectively. Away from their perturbative domains, both functions cross over into the same universal power-law scaling form governed by the static and dynamic critical exponents (as well as an exponent characterizing the quench protocol). As a by-product of the scaling studies, we obtain high-precision estimates of the dynamic exponent $z$ for the two-dimensional Ising model subject to the three variants of Monte Carlo dynamics; for single-spin Metropolis updates $z_M=2.1767(5)$, for Swendsen-Wang multi-cluster updates $z_{\rm SW}=0.297(3)$, and for Wolff single-cluster updates $z_{W}=0.30(2)$. For Wolff dynamics, we find an interesting behavior with a non-analytic break-down of the quasi-adiabatic and diabatic scaling, instead of the generic smooth cross-over described by a power law. We interpret this disconnect between the two scaling regimes as a dynamic phase transition of the Wolff algorithm, caused by an effective sudden loss of ergodicity at high velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Confinement and Deconfinement of Spinons in Two Dimensions

We use Monte Carlo methods to study spinons in two-dimensional quantum spin systems, characterizing their intrinsic size $λ$ and confinement length $Λ$. We confirm that spinons are deconfined, $Λ\to \infty$ and $λ$ finite, in a resonating valence-bond spin-liquid state. In a valence-bond solid, we find finite $λ$ and $Λ$, with $λ$ of a single spinon significantly larger than the bound-state---the spinon is soft and shrinks as the bound state is formed. Both $λ$ and $Λ$ diverge upon approaching the critical point separating valence-bond solid and Néel ground states. We conclude that the spinon deconfinement is marginal in the lowest-energy state in the spin-1 sector, due to weak attractive spinon interactions. Deconfinement in the vicinity of the critical point should occur at higher energies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quasi-adiabatic quantum Monte Carlo algorithm for quantum evolution in imaginary time

We propose a quantum Monte Carlo (QMC) algorithm for non-equilibrium dynamics in a system with a parameter varying as a function of time. The method is based on successive applications of an evolving Hamiltonian to an initial state and delivers results for a whole range of the tuning parameter in a single run, allowing for access to both static and dynamic properties of the system. This approach reduces to the standard Schrödinger dynamics in imaginary time for quasi-adiabatic evolutions, i.e., including the leading non-adiabatic correction to the adiabatic limit. We here demonstrate this quasi-adiabatic QMC (QAQMC) method for linear ramps of the transverse-field Ising model across its quantum-critical point in one and two dimensions. The critical behavior can be described by generalized dynamic scaling. For the two-dimensional square-lattice system we use the method to obtain a high-precision estimate of the quantum-critical point $(h/J)_c=3.04463(12)$, where $h$ is the transverse magnetic field and $J$ the nearest-neighbor Ising coupling. The QAQMC method can also be used to extract the Berry curvature and the metric tensor.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermal valence-bond-solid transition of quantum spins in two dimensions

We study the S=1/2 Heisenberg (J) model on the two-dimensional square lattice in the presence of additional higher-order spin interactions (Q) which lead to a valence-bond-solid (VBS) ground state. Using quantum Monte Carlo simulations, we analyze the thermal VBS transition. We find continuously varying exponents, with the correlation-length exponent "nu" close to the Ising value for large Q/J and diverging when Q/J approaches the quantum-critical point (the critical temperature Tc -> 0). This is in accord with the theory of deconfined quantum-critical points, which predicts that the transition should approach a Kosterlitz-Thouless (KT) fixed point when Tc -> 0+ (while the transition versus Q/J for T=0 is in a different class). We find explicit evidence for KT physics by studying the emergence of U(1) symmetry of the order parameter at T=T_c when Tc -> 0.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Antiferromagnetic order in systems with doublet $S_{\rm tot}=1/2$ ground states

We use projector Quantum Monte-Carlo methods to study the $S_{\rm tot}=1/2$ doublet ground states of two dimensional $S=1/2$ antiferromagnets on a $L \times L$ square lattice with an odd number of sites $N_{\rm tot}=L^2$. We compute the ground state spin texture $Φ^z(\vec{r}) = <S^z(\vec{r})>_{\uparrow}$ in $|G>_{\uparrow}$, the $S^z_{\rm tot}=1/2$ component of this doublet, and investigate the relationship between $n^z$, the thermodynamic limit of the staggered component of this ground state spin texture, and $m$, the thermodynamic limit of the magnitude of the staggered magnetization vector of the same system in the singlet ground state that obtains for even $N_{\rm tot}$. We find a univeral relationship between the two, that is independent of the microscopic details of the lattice level Hamiltonian and can be well approximated by a polynomial interpolation formula: $n^z \approx (1/3 - \frac{a}{2} -\frac{b}{4}) m + am^2+bm^3$, with $a \approx 0.288$ and $b\approx -0.306$. We also find that the full spin texture $Φ^z(\vec{r})$ is itself dominated by Fourier modes near the antiferromagnetic wavevector in a universal way. On the analytical side, we explore this question using spin-wave theory, a simple mean field model written in terms of the total spin of each sublattice, and a rotor model for the dynamics of $\vec{n}$. We find that spin-wave theory reproduces this universality of $Φ^z(\vec{r})$ and gives $n^z = (1-α-β/S)m + (α/S)m^2 +{\mathcal O}(S^{-2})$ with $α\approx 0.013$ and $β\approx 1.003$ for spin-$S$ antiferromagnets, while the sublattice-spin mean field theory and the rotor model both give $n^z = 1/3 m$ for $S=1/2$ antiferromagnets. We argue that this latter relationship becomes asymptotically exact in the limit of infinitely long-range {\em unfrustrated} exchange interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ashkin-Teller criticality and pseudo first-order behavior in a frustrated Ising model on the square lattice

We study the challenging thermal phase transition to stripe order in the frustrated square-lattice Ising model with couplings J1<0 (nearest-neighbor, ferromagnetic) and J2>0 (second-neighbor, antiferromagnetic) for g=J2/|J1|>1/2. Using Monte Carlo simulations and known analytical results, we demonstrate Ashkin-Teller criticality for g>= g*, i.e., the critical exponents vary continuously between those of the 4-state Potts model at g=g* and the Ising model for g -> infinity. Thus, stripe transitions offer a route to realizing a related class of conformal field theories with conformal charge c=1 and varying exponents. The transition is first-order for g<g*= 0.67(1), much lower than previously believed, and exhibits pseudo first-order behavior for g* < g < 1.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bridging lattice-scale physics and continuum field theory with quantum Monte Carlo simulations

We discuss designer Hamiltonians---lattice models tailored to be free from sign problems ("de-signed") when simulated with quantum Monte Carlo methods but which still host complex many-body states and quantum phase transitions of interest in condensed matter physics. We focus on quantum spin systems in which competing interactions lead to non-magnetic ground states. These states and the associated quantum phase transitions can be studied in great detail, enabling direct access to universal properties and connections with low-energy effective quantum field theories. As specific examples, we discuss the transition from a Neel antiferromagnet to either a uniform quantum paramagnet or a spontaneously symmetry-broken valence-bond solid in SU(2) and SU(N) invariant spin models. We also discuss anisotropic (XXZ) systems harboring topological Z2 spin liquids and the XY* transition. We briefly review recent progress on quantum Monte Carlo algorithms, including ground state projection in the valence-bond basis and direct computation of the Renyi variants of the entanglement entropy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Correlated valence-bond states

We study generalizations of the singlet-sector amplitude-product (AP) states in the valence-bond basis of S=1/2 quantum spin systems. In the standard AP states, the weight of a tiling of the system into valence bonds (singlets of two spins) is a product of amplitudes depending on the length of the bonds. We here introduce correlated AP (CAP) states, in which the amplitude product is further multiplied by factors depending on two bonds connected to a pair of sites (here nearest neighbors). While the standard AP states can describe a phase transition between an antiferromagnetic (Neel) state and a valence-bond solid (VBS) in one dimension (which we also study here), in two dimensions it cannot describe VBS order. With the CAP states, Neel-VBS transitions are realized as a function of some parameter describing the bond correlations. We here study such phase transitions of CAP wave-functions on the square lattice. We find examples of direct first-order Neel-VBS transitions, as well as cases where there is an extended U(1) spin liquid phase intervening between the Neel and VBS states. In the latter case the transitions are continuous and we extract critical exponents and address the issue of a possible emergent U(1) symmetry in the near-critical VBS. We also consider variationally optimized CAP states for the standard Heisenberg model in one and two dimensions and the J-Q model in two dimensions, with the latter including four-spin interactions (Q) in addition to the Heisenberg exchange (J) and harboring VBS order for large Q/J. The optimized CAP states lead to significantly lower variational energies than the simple AP states for these models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite-size scaling and boundary effects in two-dimensional valence-bond-solids

Various lattice geometries and boundaries are used to investigate valence-bond-solid (VBS) ordering in the ground state of an S=1/2 square-lattice quantum spin model---the J-Q model, in which 4- or 6-spin interactions Q are added to the Heisenberg exchange J. Ground state results for finite systems (with up to thousands of spins) are obtained using a projector QMC method. Great care has to be taken when extrapolating the order parameter to infinite size, in particular in cylinder geometry. Even though strong 2D VBS order exists and is established clearly with increasing system size on L*L lattices (or Lx* Ly lattices with a fixed Lx/Ly), only short-range VBS correlations are observed on long cylinders (when Lx -> infinity at fixed Ly). The correlation length increases with Ly, until long-range order sets in at a "critical" Ly. This width is large even when the 2D order is strong, e.g, for a system where the order parameter is 70% of the largest possible value, Ly=8 is required for ordering. Correlation functions for small L*L lattices can also be misleading. For a 20%-ordered system results for L up to 20 appear to extrapolate to a vanishing order parameter, while for larger L the behavior crosses over and extrapolates to a non-zero value (with exponentially small finite size corrections). The VBS order also exhibits interesting edge effects related to emergent U(1) symmetry, which, if not considered properly, can lead to wrong conclusions for the thermodynamic limit. The finite-size behavior for small L*L lattices and long cylinders is similar to that predicted for a Z2 spin liquid. The results raise concerns about recent works claiming Z2 spin liquid ground states in frustrated 2D systems, in particular, the Heisenberg model with nearest and next-nearest-neighbor couplings. VBS state in this system cannot be ruled out.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase transitions in the frustrated Ising model on the square lattice

We consider the thermal phase transition from a paramagnetic to stripe-antiferromagnetic phase in the frustrated two-dimensional square-lattice Ising model with competing interactions J1<0 (nearest neighbor, ferromagnetic) and J2 >0 (second neighbor, antiferromagnetic). The striped phase breaks a Z4 symmetry and is stabilized at low temperatures for g=J2/|J1|>1/2. Despite the simplicity of the model, it has proved difficult to precisely determine the order and the universality class of the phase transitions. This was done convincingly only recently by Jin et al. [PRL 108, 045702 (2012)]. Here, we further elucidate the nature of these transitions and their anomalies by employing a combination of cluster mean-field theory, Monte Carlo simulations, and transfer-matrix calculations. The J1-J2 model has a line of very weak first-order phase transitions in the whole region 1/2<g<g*, where g* = 0.67(1). Thereafter, the transitions from g above g* are continuous and can be fully mapped, using universality arguments, to the critical line of the well known Ashkin-Teller model from its 4-state Potts point to the decoupled Ising limit. We also comment on the pseudo-first-order behavior at the Potts point and its neighborhood in the Ashkin-Teller model on finite lattices, which in turn leads to the appearance of similar effects in the vicinity of the multicritical point g* in the J1-J2 model. The continuous transitions near g* can therefore be mistaken to be first-order transitions, and this realization was the key to understanding the paramagnetic-striped transition for the full range of g>1/2. Most of our results are based on Monte Carlo calculations, while the cluster mean-field and transfer-matrix results provide useful methodological bench-marks for weakly first-order behaviors and Ashkin-Teller criticality.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Universal Neel Temperature in Three-Dimensional Quantum Antiferromagnets

We study three-dimensional dimerized S=1/2 Heisenberg antiferromagnets, using quantum Monte Carlo simulations of systems with three different dimerization patterns. We propose a way to relate the Néel temperature T_N to the staggered moment m_s of the ground state. Mean-field arguments suggest that T_N is proportional to m_s close to a quantum-critical point. We find an almost perfect universality (including the prefactor) if T_N is normalized by a proper lattice-scale energy. We show that the temperature T* at which the magnetic susceptibility has a maximum is a good choise, i.e., T_N/T* versus m_s is a universal function (also beyond the linear regime). These results are useful for analyzing experiments on systems where the spin couplings are not known precisely, e.g., TlCuCl3.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A lattice model for the SU(N) Neel-VBS quantum phase transition at large N

We generalize the SU(N=2) $S=1/2$ square-lattice quantum magnet with nearest-neighbor antiferromagnetic coupling ($J_1$) and next-nearest-neighbor ferromagnetic coupling ($J_2$) to arbitrary $N$. For all $N>4$, the ground state has valence-bond-solid (VBS) order for $J_2=0$ and Néel order for $J_2/J_1\gg 1$, allowing us access to the transition between these types of states for large $N$. Using quantum Monte Carlo simulations, we show that both order parameters vanish at a single quantum-critical point, whose universal exponents for large enough $N$ (here up to N=12) approach the values obtained in a 1/N expansion of the non-compact CP$^{N-1}$ field theory. These results lend strong support to the deconfined quantum-criticality theory of the Néel--VBS transition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Computational Studies of Quantum Spin Systems

These lecture notes introduce quantum spin systems and several computational methods for studying their ground-state and finite-temperature properties. Symmetry-breaking and critical phenomena are first discussed in the simpler setting of Monte Carlo studies of classical spin systems, to illustrate finite-size scaling at continuous and first-order phase transitions. Exact diagonalization and quantum Monte Carlo (stochastic series expansion) algorithms and their computer implementations are then discussed in detail. Applications of the methods are illustrated by results for some of the most essential models in quantum magnetism, such as the S=1/2 Heisenberg antiferromagnet in one and two dimensions, as well as extended models useful for studying quantum phase transitions between antiferromagnetic and magnetically disordered states.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Plaquette Renormalization Scheme for Tensor Network States

We present a method for contracting a square-lattice tensor network in two dimensions, based on auxiliary tensors accomplishing successive truncations (renormalization) of 8-index tensors for 2 by 2 plaquettes into 4-index tensors. The scheme is variational, and thus the tensors can be optimized by minimizing the energy. Test results for the quantum phase transition of the transverse-field Ising model confirm that even the smallest possible tensors (two values for each tensor index at each renormalization level) produce much better results than the simple product (mean-field) state.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Properties of Resonating-Valence-Bond Spin Liquids and Critical Dimer Models

We use Monte Carlo simulations to study properties of Anderson's resonating-valence-bond (RVB) spin-liquid state on the square lattice (i.e., the equal superposition of all pairing of spins into nearest-neighbor singlet pairs) and compare with the classical dimer model (CDM). The latter system also corresponds to the ground state of the Rokhsar-Kivelson quantum dimer model at its critical point. We find that although spin-spin correlations decay exponentially in the RVB, four-spin valence-bond-solid (VBS) correlations are critical, qualitatively like the well-known dimer-dimer correlations of the CDM, but decaying more slowly (as $1/r^a$ with $a \approx 1.20$, compared with $a=2$ for the CDM). We also compute the distribution of monomer (defect) pair separations, which decay by a larger exponent in the RVB than in the CDM. We further study both models in their different winding number sectors and evaluate the relative weights of different sectors. Like the CDM, all the observed RVB behaviors can be understood in the framework of a mapping to a "height" model characterized by a gradient-squared stiffness constant $K$. Four independent measurements consistently show a value $K_{RVB} \approx 1.6 K_{CDM}$, with the same kinds of numerical evaluations of $K_{CDM}$ give results in agreement with the rigorously known value $K_{CDM}=π/16$. The background of a nonzero winding number gradient $W/L$ introduces spatial anisotropies and an increase in the effective K, both of which can be understood as a consequence of anharmonic terms in the height-model free energy, which are of relevance to the recently proposed scenario of "Cantor deconfinement" in extended quantum dimer models. We also study ensembles in which fourth-neighbor (bipartite) bonds are allowed, at a density controlled by a tunable fugacity, resulting (as expected) in a smooth reduction of K.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Continuous quantum phase transition between an antiferromagnet and a valence-bond-solid in two dimensions; evidence for logarithmic corrections to scaling

The antiferromagnetic to valence-bond-solid phase transition in the two-dimensional J-Q model (an S=1/2 Heisenberg model with four-spin interactions) is studied using large-scale quantum Monte Carlo simulations. The results support a continuous transition of the ground state, in agreement with the theory of "deconfined" quantum criticality. There are, however, large corrections to scaling, of logarithmic or very slowly decaying power-law form, which had not been anticipated. This suggests that either the SU($N$) symmetric noncompact CP^(N-1) field theory for deconfined quantum criticality has to be revised, or that the theory for N=2 (as in the system studied here) differs significantly from N -> infinity (where the field theory is analytically tractable).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Definitions of entanglement entropy of spin systems in the valence-bond basis

The valence-bond structure of spin-1/2 Heisenberg antiferromagnets is closely related to quantum entanglement. We investigate measures of entanglement entropy based on transition graphs, which characterize state overlaps in the overcomplete valence-bond basis. The transition graphs can be generated using projector Monte Carlo simulations of ground states of specific hamiltonians or using importance-sampling of valence-bond configurations of amplitude-product states. We consider definitions of entanglement entropy based on the bonds or loops shared by two subsystems (bipartite entanglement). Results for the bond-based definition agrees with a previously studied definition using valence-bond wave functions (instead of the transition graphs, which involve two states). For the one dimensional Heisenberg chain, with uniform or random coupling constants, the prefactor of the logarithmic divergence with the size of the smaller subsystem agrees with exact results. For the ground state of the two-dimensional Heisenberg model (and also Neel-ordered amplitude-product states), there is a similar multiplicative violation of the area law. In contrast, the loop-based entropy obeys the area law in two dimensions, while still violating it in one dimension - both behaviors in accord with expectations for proper measures of entanglement entropy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effects of edge disorder in nano-scale antiferromagnetic clusters

We study the distribution of local magnetic susceptibilities in the two-dimensional antiferromagnetic S=1/2 Heisenberg model on various random clusters, in order to determine whether effects of edge disorder could be detected in NMR experiments (through the line shape, as given by the distribution of local Knight shifts). Although the effects depend strongly on the nature of the edge and the cluster size, our results indicate that line widths broader than the average shift should be expected even in clusters as large as $\approx 1000$ lattice spacing in diameter. Experimental investigations of the NMR line width should give insights into the magnetic structure of the edges.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Example of a first-order Néel to Valence-Bond-Solid transition in two-dimensions

We consider the $S=1/2$ Heisenberg model with nearest-neighbor interaction $J$ and an additional multi-spin interaction $Q_3$ on the square lattice. The $Q_3$ term consists of three bond-singlet projectors and is chosen to favor the formation of a valence-bond solid (VBS) where the valence bonds (singlet pairs) form a staggered pattern. The model exhibits a quantum phase transition from the Néel state to the VBS as a function of $Q_3/J$. We study the model using quantum Monte Carlo (stochastic series expansion) simulations. The Néel-VBS transition in this case is strongly first-order in nature, in contrast to similar previously studied models with continuous transitions into columnar VBS states. The qualitatively different transitions illustrate the important role of an emerging U(1) symmetry in the latter case, which is not possible in the present model due to the staggered VBS pattern (which does not allow local fluctuations necessary to rotate the local VBS order parameter).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ground states of a frustrated quantum spin chain with long-range interactions

The ground state of a spin-1/2 Heisenberg chain with both frustration and long-range interactions is studied using Lanczos exact diagonalization. The evolution of the well known dimerization transition of the system with short-range frustrated interactions (the J1-J2 chain) is investigated in the presence of additional unfrustrated interactions decaying with distance as 1/r^a. It is shown that the continuous (infinite-order) dimerization transition develops into a first-order transition between a long-range ordered antiferromagnetic state and a state with coexisting dimerization and critical spin correlations at wave-number k=π/2. The relevance of the model to real systems is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum phase transitions in disordered dimerized quantum spin models and the Harris criterion

We use quantum Monte Carlo simulations to study effects of disorder on the quantum phase transition occurring versus the ratio g=J/J' in square-lattice dimerized S=1/2 Heisenberg antiferromagnets with intra- and inter-dimer couplings J and J'. The dimers are either randomly distributed (as in the classical dimer model), or come in parallel pairs with horizontal or vertical orientation. In both cases the transition violates the Harris criterion, according to which the correlation-length exponent should satisfy nu >= 1. We do not detect any deviations from the three-dimensional O(3) universality class obtaining in the absence of disorder (where nu = 0.71). We discuss special circumstances which allow nu<1 for the type of disorder considered here.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Symmetry breaking and criticality in tensor-product states

We discuss variationally optimized matrix-product states for the transverse-field Ising chain, using D*D matrices with small D=2-10. For finite system size N there are energy minimums for symmetric as well as symmetry-broken states, which cross each other at a field value hc(N,D); thus the transition is first-order. A continuous transition develops as N->infinity. The asymptotic critical behavior is then always of mean-field type (the magnetization exponent beta=1/2), but a window of field strengths where true Ising scaling holds (beta=1/8) emerges with increasing D. We also demonstrate asymptotic mean-field behavior for infinite-size two-dimensional tensor-product (iPEPS) states with small tensors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermodynamics of a gas of deconfined bosonic spinons in two dimensions

We consider the quantum phase transition between a Neel antiferromagnet and a valence-bond solid (VBS) in a two-dimensional system of S=1/2 spins. Assuming that the excitations of the critical ground state are linearly dispersing deconfined spinons obeying Bose statistics, we derive expressions for the specific heat and the magnetic susceptibility at low temperature T. Comparing with quantum Monte Carlo results for the J-Q model, which is a candidate for a deconfined Neel-VBS transition, we find excellent agreement, including a previously noted logarithmic correction in the susceptibility. In our treatment, this is a direct consequence of a confinement length scale Lambda which is proportional to the correlation length xi raised to a non-trivial power; Lambda ~ xi^(1+a) ~1/T^(1+a), with a>0 (with a approximately 0.2 in the model).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nature of the low-energy excitations of two-dimensional diluted Heisenberg quantum antiferromagnets

We present several different calculations pertaining to the nature of the low-energy excitations of the site-diluted S=1/2 Heisenberg antiferromagnet, in particular at the percolation point. We present a picture of excitations originating from an effective low-energy subsystem consisting of localized magnetic moments. At the percolation point, these moments lead to an anomalously large dynamic exponent; z=3.6(1). The magnetic moments are shown to originate from local sublattice imbalance, which we study quantitatively using a classical dimer-monomer model. The localization properties of triplet excitations of clusters with singlet ground state are examined using simulations in the valence bond basis. The triplet is shown to affect predominantly sites on the classical monomer regions. The number of sites affected grows as a non-trivial power of the cluster size. We also study a bilayer model, where there is no sublattice imbalance. Accordingly, we find a much smaller dynamic exponent, consistent with quantum rotor (or possibly fracton) excitations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Edge effects in the two-dimensional spin-1/2 Heisenberg antiferromagnet

We use quantum Monte Carlo simulations to study effects of free edges in the two-dimensional spin-1/2 Heisenberg antiferromagnet. We find that the magnetic response of an edge is smaller than the bulk susceptibility. This counter-intuitive quantum effect can be traced to enhanced antiferromagnetic nearest-neighbor spin correlations, i.e., tendency to local singlet formation, at and close to the edge. These correlations form a comb-like pattern, which can be reproduced with a simple variational valence-bond state. We also study rough edges, and find that these instead significantly enhance the susceptibility, due to local sublattice imbalance impeding singlet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Z4-U(1) crossover of the order parameter symmetry in a two-dimensional valence-bond-solid

We discuss ground-state projector simulations of a modified two-dimensional S=1/2 Heisenberg model in the valence bonds basis. Tuning matrix elements corresponding to the diagonal and off-diagonal terms in the quantum dimer model, we show that there is a quantum phase transition from the antiferromagnet into a columnar valence-bond-solid (VBS). There are no signs of discontinuities, suggesting a continuous or very weakly first-order transition. The Z4-symmetric VBS order parameter exhibits an emergent U(1) symmetry as the phase transition is approached. We extract the associated length-scale governing the U(1)-Z4 cross-over inside the VBS phase.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Emergence of U(1) symmetry in the 3D XY model with Zq anisotropy

We study the three-dimensional XY model with a Z_q anisotropic term. At temperatures T < Tc this dangerously irrelevant perturbation is relevant only above a length scale Lambda, which diverges as a power of the correlation length; Lambda ~ xi^a_q. Below Lambda the order parameter is U(1) symmetric. We derive the full scaling function controlling the emergence of U(1) symmetry and use Monte Carlo results to extract the exponent a_q for q=4,...,8. We find that a_q = a_4 (q/4)^2, with a_4 only marginally larger than 1. We discuss these results in the context of U(1) symmetry at "deconfined" quantum critical points separating antiferromagnetic and valence-bond-solid states in quantum spin systems.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Scale-renormalized matrix-product states for correlated quantum systems

A generalization of matrix product states (MPS) is introduced which is suitable for describing interacting quantum systems in two and three dimensions. These scale-renormalized matrix-product states (SR-MPS) are based on a course-graining of the lattice in which the blocks at each level are associated with matrix products that are further transformed (scale renormalized) with other matrices before they are assembled to form blocks at the next level. Using variational Monte Carlo simulations of the two-dimensional transverse-field Ising model as a test, it is shown that the SR-MPS converge much more rapidly with the matrix size than a standard MPS. It is also shown that the use of lattice-symmetries speeds up the convergence very significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Low-energy dynamics of the two-dimensional S=1/2 Heisenberg antiferromagnet on percolating clusters

We investigate the quantum dynamics of site diluted S=1/2 Heisenberg antiferromagnetic clusters at the percolation threshold. We use Lanczos diagonalization to calculate the lowest excitation gap Delta and, to reach larger sizes, study an upper bound for Delta obtained from sum rules involving the staggered structure factor and susceptibility, which we evaluate by quantum Monte Carlo simulations. Scaling the gap distribution with the cluster length L, Delta sim 1/L^z, we obtain a dynamic exponent z approximate 2D_f, where D_f=91/48 is the fractal dimensionality of the percolating cluster. This is in contrast to previous expectations of z=D_f. We argue that the low-energy excitations are due to weakly coupled effective moments formed due to local imbalance in sublattice occupation.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Ground state projection of quantum spin systems in the valence bond basis

A Monte Carlo method for quantum spin systems is formulated in the basis of valence bond (singlet pair) states. The non-orthogonality of this basis allows for an efficient importance-sampled projection of the ground state out of an arbitrary state. The method provides access to resonating valence-bond physics, enables a direct improved estimator for the singlet-triplet gap, and extends the class of models that can be studied without negative-sign problems. As a demonstration, the valence bond distribution in the ground state of the 2D Heisenberg antiferromagnet is calculated. Generalizations of the method to fermion systems are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

High-precision finite-size scaling analysis of the quantum-critical point of S=1/2 Heisenberg antiferromagnetic bilayers

We use quantum Monte Carlo (stochastic series expansion) and finite-size scaling to study the quantum critical points of two S=1/2 Heisenberg antiferromagnets in two dimensions: a bilayer and a Kondo-lattice-like system (incomplete bilayer), each with intra- and inter-plane couplings J and J_perp. We discuss the ground-state finite-size scaling properties of three different quantities--the Binder moment ratio, the spin stiffness, and the long-wavelength magnetic susceptibility--which we use to extract the critical value of the coupling ratio g=J_perp/J. The individual estimates of g_c are consistent provided that subleading finite-size corrections are properly taken into account. In the case of the complete bilayer, the Binder ratio leads to the most precise estimate of the critical coupling, although the subleading finite-size corrections to the stiffness are considerably smaller. For the incomplete bilayer, the subleading corrections to the stiffness are extremely small, and this quantity then gives the best estimate of the critical point. Contrary to predictions, we do not find a universal prefactor of the 1/L spin stiffness scaling at the critical point, whereas the Binder ratio is consistent with a universal value. Our results for the critical coupling ratios are g_c=2.52181(3) (full bilayer) and g_c=1.38882(2) (incomplete bilayer), which represent improvements of two orders of magnitude relative to the previous best estimates. For the correlation length exponent we obtain nu = 0.7106(9), consistent with the expected 3D Heisenberg universality class.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Quantum-Critical Behavior in a Two-Layer Antiferromagnet

We analyze quantum Monte Carlo data in the vicinity of the quantum transition between a Neel state and a quantum paramagnet in a two-layer, square lattice spin 1/2 Heisenberg antiferromagnet. The real-space correlation function and the universal amplitude ratio of the structure factor and the dynamic susceptibility show clear evidence of quantum critical behavior at low temperatures. The numerical results are in good quantitative agreement with $1/N$ calculations for the $O(N)$ non-linear sigma model. A discrepancy, reported earlier, between the critical properties of the antiferromagnet and the sigma model is resolved. We also discuss the values of prefactors of the dynamic susceptibility and the structure factor in a single layer antiferromagnet at low $T$.

0 citations0 reviewsNo code linkedOpen access