Researcher profile

Anders W. Sandvik

Anders W. Sandvik contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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.

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.

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.

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.

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 &#34;elliptical&#34; 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.

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.

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.

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.

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.

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.

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.

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.

preprint2019arXiv

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

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