91̽��

Abstract:

We analyze the effect of quenched disorder on spin-1/2 quantum magnets in which magnetic frustration promotes the formation of local singlets. Our results include a theory for 2D valence-bond solids subject to weak bond randomness, as well as extensions to stronger disorder regimes where we make connections with quantum spin liquids. We find, on various lattices, that the destruction of a valence-bond solid phase by weak quenched disorder leads inevitably to the nucleation of topological defects carrying spin-1/2 moments. This renormalizes the lattice into a strongly random spin network with interesting low-energy excitations. Similarly, when short-ranged valence bonds would be pinned by stronger disorder, we find that this putative glass is unstable to defects that carry spin-1/2 magnetic moments, and whose residual interactions decide the ultimate low-energy fate. Motivated by these results we conjecture Lieb-Schultz-Mattis-like restrictions on ground states for disordered magnets with spin 1/2 per statistical unit cell. These conjectures are 91̽��ed by an argument for 1D spin chains. We apply insights from this study to the phenomenology of YbMgGaO4, a recently discovered triangular lattice spin-1/2 insulator which was proposed to be a quantum spin liquid. We instead explore a description based on the present theory. Experimental signatures, including unusual specific heat, thermal conductivity, and dynamical structure factor, and their behavior in a magnetic field, are predicted from the theory, and compare favorably with existing measurements on YbMgGaO4 and related materials.

Abstract:

Quenched randomness can have a dramatic effect on the dynamics of isolated 1D quantum many-body systems, even for systems that thermalize. This is because transport, entanglement, and operator spreading can be hindered by “Griffiths” rare regions, which locally resemble the many-body-localized phase and thus act as weak links. We propose coarse-grained models for entanglement growth and for the spreading of quantum operators in the presence of such weak links. We also examine entanglement growth across a single weak link numerically. We show that these weak links have a stronger effect on entanglement growth than previously assumed: entanglement growth is subballistic whenever such weak links have a power-law probability distribution at low couplings, i.e., throughout the entire thermal Griffiths phase. We argue that the probability distribution of the entanglement entropy across a cut can be understood from a simple picture in terms of a classical surface growth model. We also discuss spreading of operators and conserved quantities. Surprisingly, the four length scales associated with (i) production of entanglement, (ii) spreading of conserved quantities, (iii) spreading of operators, and (iv) the width of the “front” of a spreading operator, are characterized by dynamical exponents that in general are all distinct. Our numerical analysis of entanglement growth between weakly coupled systems may be of independent interest.

Abstract:

Random quantum circuits yield minimally structured models for chaotic quantum dynamics, which are able to capture, for example, universal properties of entanglement growth. We provide exact results and coarse-grained models for the spreading of operators by quantum circuits made of Haar-random unitaries. We study both 1+1D and higher dimensions and argue that the coarse-grained pictures carry over to operator spreading in generic many-body systems. In 1+1D, we demonstrate that the out-of-time-order correlator (OTOC) satisfies a biased diffusion equation, which gives exact results for the spatial profile of the OTOC and determines the butterfly speed vB. We find that in 1+1D, the “front” of the OTOC broadens diffusively, with a width scaling in time as t1/2. We address fluctuations in the OTOC between different realizations of the random circuit, arguing that they are negligible in comparison to the broadening of the front within a realization. Turning to higher dimensions, we show that the averaged OTOC can be understood exactly via a remarkable correspondence with a purely classical droplet growth problem. This implies that the width of the front of the averaged OTOC scales as t1/3 in 2+1D and as t0.240 in 3+1D (exponents of the Kardar-Parisi-Zhang universality class). We 91̽�� our analytic argument with simulations in 2+1D. We point out that, in two or higher spatial dimensions, the shape of the spreading operator at late times is affected by underlying lattice symmetries and, in general, is not spherical. However, when full spatial rotational symmetry is present in 2+1D, our mapping implies an exact asymptotic form for the OTOC, in terms of the Tracy-Widom distribution. For an alternative perspective on the OTOC in 1+1D, we map it to the partition function of an Ising-like statistical mechanics model. As a result of special structure arising from unitarity, this partition function reduces to a random walk calculation which can be performed exactly. We also use this mapping to give exact results for entanglement growth in 1+1D circuits.