Source author record

David A. Huse

David A. Huse appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

78works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

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

78 published item(s)

preprint2026arXiv

Eigenstate thermalization in thermal first-order phase transitions

The eigenstate thermalization hypothesis (ETH) posits how isolated quantum many-body systems thermalize, assuming that individual eigenstates at the same energy density have identical expectation values of local observables in the limit of large systems. While the ETH apparently holds across a wide range of interacting quantum systems, in this work we show that it requires generalization in the presence of thermal first-order phase transitions. We introduce a class of all-to-all spin models, featuring first-order thermal phase transitions that stem from two distinct mean-field solutions (two ``branches'') that exchange dominance in the many-body density of states as the energy is varied. We argue that for energies in the vicinity of the thermal phase transition, eigenstate expectation values do not need to converge to the same thermal value. The system has a regime with coexistence of two classes of eigenstates corresponding to the two branches with distinct expectation values at the same energy density, and another regime with Schrodinger-cat-like eigenstates that are inter-branch superpositions; these two regimes are separated by an eigenstate phase transition. We support our results by semiclassical calculations and an exact diagonalization study of a microscopic spin model, and argue that the structure of eigenstates in the vicinity of thermal first-order phase transitions can be experimentally probed via non-equilibrium dynamics.

preprint2025arXiv

Crystalline Spectral Form Factors

We investigate crystalline-like behavior of the spectral form factor (SFF) in unitary quantum systems with extremely strong eigenvalue repulsion. Using a low-temperature Coulomb gas as a model of repulsive eigenvalues, we derive the Debye-Waller factor suppressing periodic oscillations of the SFF and estimate the order of its singularities at multiples of the Heisenberg time. We also reproduce this crystalline-like behavior using perturbed permutation circuits and random matrix ensembles associated with Lax matrices. Our results lay a foundation for future studies of quantum systems that exhibit intermediate level statistics between standard random matrix ensembles and permutation circuits.

preprint2023arXiv

Non-Abelian symmetry can increase entanglement entropy

The pillars of quantum theory include entanglement and operators' failure to commute. The Page curve quantifies the bipartite entanglement of a many-body system in a random pure state. This entanglement is known to decrease if one constrains extensive observables that commute with each other (Abelian ``charges''). Non-Abelian charges, which fail to commute with each other, are of current interest in quantum thermodynamics. For example, noncommuting charges were shown to reduce entropy-production rates and may enhance finite-size deviations from eigenstate thermalization. Bridging quantum thermodynamics to many-body physics, we quantify the effects of charges' noncommutation -- of a symmetry's non-Abelian nature -- on Page curves. First, we construct two models that are closely analogous but differ in whether their charges commute. We show analytically and numerically that the noncommuting-charge case has more entanglement. Hence charges' noncommutation can promote entanglement.

preprint2022arXiv

Avalanches and many-body resonances in many-body localized systems

We numerically study both the avalanche instability and many-body resonances in strongly-disordered spin chains exhibiting many-body localization (MBL). We distinguish between a finite-size/time MBL regime, and the asymptotic MBL phase, and identify some "landmarks" within the MBL regime. Our first landmark is an estimate of where the MBL phase becomes unstable to avalanches, obtained by measuring the slowest relaxation rate of a finite chain coupled to an infinite bath at one end. Our estimates indicate that the actual MBL-to-thermal phase transition, in infinite-length systems, occurs much deeper in the MBL regime than has been suggested by most previous studies. Our other landmarks involve system-wide resonances. We find that the effective matrix elements producing eigenstates with system-wide resonances are enormously broadly distributed. This means that the onset of such resonances in typical samples occurs quite deep in the MBL regime, and the first such resonances typically involve rare pairs of eigenstates that are farther apart in energy than the minimum gap. Thus we find that the resonance properties define two landmarks that divide the MBL regime in to three subregimes: (i) at strongest disorder, typical samples do not have any eigenstates that are involved in system-wide many-body resonances; (ii) there is a substantial intermediate regime where typical samples do have such resonances, but the pair of eigenstates with the minimum spectral gap does not; and (iii) in the weaker randomness regime, the minimum gap is involved in a many-body resonance and thus subject to level repulsion. Nevertheless, even in this third subregime, all but a vanishing fraction of eigenstates remain non-resonant and the system thus still appears MBL in many respects. Based on our estimates of the location of the avalanche instability, it might be that the MBL phase is only part of subregime (i).

preprint2022arXiv

Height-conserving quantum dimer models

We propose a height-conserving quantum dimer model (hQDM) such that the lattice sum of its associated height field is conserved, and that it admits a Rokhsar-Kivelson (RK) point. The hQDM with minimal interaction range on the square lattice exhibits Hilbert space fragmentation and maps exactly to the XXZ spin model on the square lattice in certain Krylov subspaces. We obtain the ground-state phase diagram of hQDM via quantum Monte Carlo simulations, and demonstrate that a large portion of it is within the Krylov subspaces which admit the exact mapping to the XXZ model, with dimer ordered phases corresponding to easy-axis and easy-plane spin orders. At the RK point, the apparent dynamical exponents obtained from the single mode approximation and the height correlation function show drastically different behavior across the Krylov subspaces, exemplifying Hilbert space fragmentation and emergent glassy phenomena.

preprint2022arXiv

Maximum entropy quantum state distributions

We propose an approach to the realization of many-body quantum state distributions inspired by combined principles of thermodynamics and mesoscopic physics. Its essence is a maximum entropy principle conditioned by conservation laws. We go beyond traditional thermodynamics and condition on the full distribution of the conserved quantities. The result are quantum state distributions whose deviations from `thermal states' get more pronounced in the limit of wide input distributions. We describe their properties in terms of entanglement measures and discuss strategies for state engineering by methods of current date experimentation.

preprint2022arXiv

Onset of many-body quantum chaos due to breaking integrability

Integrable quantum systems of finite size are generically robust against weak enough integrability-breaking perturbations, but become quantum chaotic and thermalizing if the integrability-breaking is strong enough. We argue that the onset of quantum chaos can be described as a Fock-space delocalization process, with the eigenstates of the integrable system being taken as the "Fock states". The integrability-breaking perturbation introduces hopping in this Fock space, and chaos sets in when this hopping delocalizes the many-body eigenstates in this space. Depending on the range of the dominant Fock-space hopping, delocalization can occur either through a crossover, or via a transition that becomes sharp in the appropriate large-system dynamic limit. In either case, the perturbation strength at the onset of chaos scales to zero in the usual thermodynamic limit, with a size-dependence that we estimate analytically and compute numerically for a few specific models. We also identify two intermediate finite-size-dependent regimes: There is generally an intermediate nonchaotic regime in which integrability is broken strongly enough to produce some system-wide many-body resonances but not enough to thermalize the system. In spatially extended systems (but not in quantum dots) there is also a crossover or transition between chaotic regimes where the ratio of the system size to the mean free path of the quasiparticles of the integrable system is small versus large compared to unity.

preprint2022arXiv

Operator scaling dimensions and multifractality at measurement-induced transitions

Repeated local measurements of quantum many body systems can induce a phase transition in their entanglement structure. These measurement-induced phase transitions (MIPTs) have been studied for various types of dynamics, yet most cases yield quantitatively similar values of the critical exponents, making it unclear if there is only one underlying universality class. Here, we directly probe the properties of the conformal field theories governing these MIPTs using a numerical transfer-matrix method, which allows us to extract the effective central charge, as well as the first few low-lying scaling dimensions of operators at these critical points. Our results provide convincing evidence that the generic and Clifford MIPTs for qubits lie in different universality classes and that both are distinct from the percolation transition for qudits in the limit of large onsite Hilbert space dimension. For the generic case, we find strong evidence of multifractal scaling of correlation functions at the critical point, reflected in a continuous spectrum of scaling dimensions.

preprint2021arXiv

Entanglement phase transitions in measurement-only dynamics

Unitary circuits subject to repeated projective measurements can undergo an entanglement phase transition (EPT) as a function of the measurement rate. This transition is generally understood in terms of a competition between the scrambling effects of unitary dynamics and the disentangling effects of measurements. We find that, surprisingly, EPTs are possible even in the absence of scrambling unitary dynamics, where they are best understood as arising from measurements alone. This motivates us to introduce \emph{measurement-only models}, in which the "scrambling" and "un-scrambling" effects driving the EPT are fundamentally intertwined and cannot be attributed to physically distinct processes. This represents a novel form of an EPT, conceptually distinct from that in hybrid unitary-projective circuits. We explore the entanglement phase diagrams, critical points, and quantum code properties of some of these measurement-only models. We find that the principle driving the EPTs in these models is \emph{frustration}, or mutual incompatibility, of the measurements. Suprisingly, an entangling (volume-law) phase is the generic outcome when measuring sufficiently long but still local ($\gtrsim 3$-body) operators. We identify a class of exceptions to this behavior ("bipartite ensembles") which cannot sustain an entangling phase, but display dual area-law phases, possibly with different kinds of quantum order, separated by self-dual critical points. Finally, we introduce a measure of information spreading in dynamics with measurements and use it to demonstrate the emergence of a statistical light-cone, despite the non-locality inherent to quantum measurements.

preprint2021arXiv

Field theory of charge sharpening in symmetric monitored quantum circuits

Monitored quantum circuits (MRCs) exhibit a measurement-induced phase transition between area-law and volume-law entanglement scaling. MRCs with a conserved charge additionally exhibit two distinct volume-law entangled phases that cannot be characterized by equilibrium notions of symmetry-breaking or topological order, but rather by the non-equilibrium dynamics and steady-state distribution of charge fluctuations. These include a charge-fuzzy phase in which charge information is rapidly scrambled leading to slowly decaying spatial fluctuations of charge in the steady state, and a charge-sharp phase in which measurements collapse quantum fluctuations of charge without destroying the volume-law entanglement of neutral degrees of freedom. By taking a continuous-time, weak-measurement limit, we construct a controlled replica field theory description of these phases and their intervening charge-sharpening transition in one spatial dimension. We find that the charge fuzzy phase is a critical phase with continuously evolving critical exponents that terminates in a modified Kosterlitz-Thouless transition to the short-range correlated charge-sharp phase. We numerically corroborate these scaling predictions also hold for discrete-time projective-measurement circuit models using large-scale matrix-product state simulations, and discuss generalizations to higher dimensions.

preprint2021arXiv

Observation of measurement-induced quantum phases in a trapped-ion quantum computer

Many-body open quantum systems balance internal dynamics against decoherence from interactions with an environment. Here, we explore this balance via random quantum circuits implemented on a trapped ion quantum computer, where the system evolution is represented by unitary gates with interspersed projective measurements. As the measurement rate is varied, a purification phase transition is predicted to emerge at a critical point akin to a fault-tolerent threshold. We probe the "pure" phase, where the system is rapidly projected to a deterministic state conditioned on the measurement outcomes, and the "mixed" or "coding" phase, where the initial state becomes partially encoded into a quantum error correcting codespace. We find convincing evidence of the two phases and show numerically that, with modest system scaling, critical properties of the transition clearly emerge.

preprint2020arXiv

Avoided quantum criticality in exact numerical simulations of a single disordered Weyl cone

Existing theoretical works differ on whether three-dimensional Dirac and Weyl semimetals are stable to a short-range-correlated random potential. Numerical evidence suggests the semimetal to be unstable, while some field-theoretic instanton calculations have found it to be stable. The differences go beyond method: the continuum field-theoretic works use a single, perfectly linear Weyl cone, while numerical works use tight-binding lattice models which inherently have band curvature and multiple Weyl cones. In this work, we bridge this gap by performing exact numerics on the same model used in analytic treatments, and we find that all phenomena associated with rare regions near the Weyl node energy found in lattice models persist in the continuum theory: The density of states is non-zero and exhibits an avoided transition. In addition to characterizing this transition, we find rare states and show that they have the expected behavior. The simulations utilize sparse matrix techniques with formally dense matrices; doing so allows us to reach Hilbert space sizes upwards of $10^7$ states, substantially larger than anything achieved before.

preprint2020arXiv

Critical properties of the measurement-induced transition in random quantum circuits

We numerically study the measurement-driven quantum phase transition of Haar-random quantum circuits in $1+1$ dimensions. By analyzing the tripartite mutual information we are able to make a precise estimate of the critical measurement rate $p_c = 0.17(1)$. We extract estimates for the associated bulk critical exponents that are consistent with the values for percolation, as well as those for stabilizer circuits, but differ from previous estimates for the Haar-random case. Our estimates of the surface order parameter exponent appear different from that for stabilizer circuits or percolation, but we are unable to definitively rule out the scenario where all exponents in the three cases match. Moreover, in the Haar case the prefactor for the entanglement entropies $S_n$ depends strongly on the Rényi index $n$; for stabilizer circuits and percolation this dependence is absent. Results on stabilizer circuits are used to guide our study and identify measures with weak finite-size effects. We discuss how our numerical estimates constrain theories of the transition.

preprint2020arXiv

How much entanglement is needed to reduce the energy variance?

We explore the relation between the entanglement of a pure state and its energy variance for a local one dimensional Hamiltonian, as the system size increases. In particular, we introduce a construction which creates a matrix product state of arbitrarily small energy variance $δ^2$ for $N$ spins, with bond dimension scaling as $\sqrt{N} D_0^{1/δ}$, where $D_0>1$ is a constant. This implies that a polynomially increasing bond dimension is enough to construct states with energy variance that vanishes with the inverse of the logarithm of the system size. We run numerical simulations to probe the construction on two different models, and compare the local reduced density matrices of the resulting states to the corresponding thermal equilibrium. Our results suggest that the spatially homogeneous states with logarithmically decreasing variance, which can be constructed efficiently, do converge to the thermal equilibrium in the thermodynamic limit, while the same is not true if the variance remains constant.

preprint2020arXiv

Kinetically constrained freezing transition in a dipole-conserving system

We study a stochastic lattice gas of particles in one dimension with strictly finite-range interactions that respect the fracton-like conservation laws of total charge and dipole moment. As the charge density is varied, the connectivity of the system's charge configurations under the dynamics changes qualitatively. We find two distinct phases: Near half filling the system thermalizes subdiffusively, with almost all configurations belonging to a single dynamically connected sector. As the charge density is tuned away from half filling there is a phase transition to a frozen phase where locally active finite bubbles cannot exchange particles and the system fails to thermalize. The two phases exemplify what has recently been referred to as weak and strong Hilbert space fragmentation, respectively. We study the static and dynamic scaling properties of this weak-to-strong fragmentation phase transition in a kinetically constrained classical Markov circuit model, obtaining some conjectured exact critical exponents.

preprint2020arXiv

Many-body localization near the critical point

We examine the many-body localization (MBL) phase transition in one-dimensional quantum systems with quenched randomness and short-range interactions. Following recent works, we use a strong-randomness renormalization group (RG) approach where the phase transition is due to the so-called avalanche instability of the MBL phase. We show that the critical behavior can be determined analytically within this RG. On a rough $\textit{qualitative}$ level the RG flow near the critical fixed point is similar to the Kosterlitz-Thouless (KT) flow as previously shown, but there are important differences in the critical behavior. Thus we show that this MBL transition is in a new universality class that is different from KT. The divergence of the correlation length corresponds to critical exponent $ν\rightarrow \infty$, but the divergence is weaker than for the KT transition.

preprint2020arXiv

Nonlinear Fluctuating Hydrodynamics for the Classical XXZ Spin Chain

Using the framework of nonlinear fluctuating hydrodynamics (NFH), we examine equilibrium spatio-temporal correlations in classical ferromagnetic spin chains with nearest neighbor interactions. In particular, we consider the classical XXZ-Heisenberg spin chain (also known as Lattice Landau Lifshitz or LLL model) evolving deterministically and chaotically via Hamiltonian dynamics, for which energy and $z$-magnetization are the only locally conserved fields. For the easy-plane case, this system has a low-temperature regime in which the difference between neighboring spin's angular orientations in the XY plane is an \textit{almost conserved} field. According to the predictions of NFH, the dynamic correlations in this regime exhibit a heat peak and propagating sound peaks, all with anomalous broadening. We present a detailed molecular dynamics test of these predictions and find a reasonably accurate verification. We find that, in a suitable intermediate temperature regime, the system shows two sound peaks with Kardar-Parisi-Zhang (KPZ) scaling and a heat peak where the expected anomalous broadening is less clear. In high temperature regimes of both easy plane and easy axis case of LLL, our numerics show clear diffusive spin and energy peaks and absence of any sound modes, as one would expect. We also simulate an integrable version of the XXZ-model, for which the ballistic component instead moves with a broad range of speeds rather than being concentrated in narrower peaks around the sound speed.

preprint2020arXiv

Perturbative instability towards delocalization at phase transitions between MBL phases

We examine the stability of marginally Anderson localized phase transitions between localized phases to the addition of many-body interactions, focusing in particular on the spin-glass to paramagnet transition in a disordered transverse field Ising model in one dimension. We find evidence for a perturbative instability of localization at finite energy densities once interactions are added, i.e. evidence for the relevance of interactions - in a renormalization group sense - to the non-interacting critical point governed by infinite randomness scaling. We introduce a novel diagnostic, the "susceptibility of entanglement", which allows us to perturbatively probe the effect of adding interactions on the entanglement properties of eigenstates, and helps us elucidate the resonant processes that can cause thermalization. The susceptibility serves as a much more sensitive probe, and its divergence can detect the perturbative beginnings of an incipient instability even in regimes and system sizes for which conventional diagnostics point towards localization. We expect this new measure to be of independent interest for analyzing the stability of localization in a variety of different settings.

preprint2020arXiv

Scalable probes of measurement-induced criticality

We uncover a local order parameter for measurement-induced phase transitions: the average entropy of a single reference qubit initially entangled with the system. Using this order parameter, we identify scalable probes of measurement-induced criticality (MIC) that are immediately applicable to advanced quantum computing platforms. We test our proposal on a 1+1 dimensional stabilizer circuit model that can be classically simulated in polynomial time. We introduce the concept of a "decoding light cone" to establish the local and efficiently measurable nature of this probe. We also estimate bulk and surface critical exponents for the transition. Developing scalable probes of MIC in more general models may be a useful application of noisy-intermediate scale quantum (NISQ) devices, as well as point to more efficient realizations of fault-tolerant quantum computation.

preprint2020arXiv

Spatio-temporal spread of perturbations in power-law models at low temperatures: Exact results for OTOC

We present exact results for the classical version of the Out-of-Time-Order Commutator (OTOC) for a family of power-law models consisting of $N$ particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form $\propto \sum_{\substack{i, j=1 (i\neq j)}}^N|x_i-x_j|^{-k}$ $\forall$ $k>1$ where $x_i$ is the position of the $i^\text{th}$ particle. We present numerical results for the OTOC for finite $N$ at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many body ground state. In the large-$N$ limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the "light cone". We find remarkably distinct features in OTOC below and above $k=3$ in terms of going from non-Airy behaviour ($1<k<3$) to an Airy universality class ($k>3$). We present certain additional rich features for the case $k=2$ that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatio-temporal spread of perturbations in long-range interacting systems.

preprint2020arXiv

Transport, correlations, and chaos in a classical disordered anharmonic chain

We explore transport properties in a disordered nonlinear chain of classical harmonic oscillators and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $N$, disorder strength $Δ$, and temperature $T$. The conductivity $κ_N$, obtained for finite length ($N$) systems, saturates to a value $κ_\infty >0$ in the large $N$ limit, for all values of disorder strength $Δ$ and temperature $T>0$. We show evidence that for any $Δ>0$ the conductivity goes to zero faster than any power of $T$ in the $(T/Δ) \to 0$ limit, and find that the form $κ_\infty \sim e^{-B |\ln(C Δ/T)|^3}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multi-oscillator chaotic islands. The finite-size effect can be $κ_N < κ_{\infty}$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_{\infty}$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the Out-of-Time-Ordered-Commutator.

preprint2018arXiv

Weyl semimetal to metal phase transitions driven by quasiperiodic potentials

We explore the stability of three-dimensional Weyl and Dirac semimetals subject to quasiperiodic potentials. We present numerical evidence that the semimetal is stable for weak quasiperiodic potentials, despite being unstable for weak random potentials. As the quasiperiodic potential strength increases, the semimetal transitions to a metal, then to an "inverted" semimetal, and then finally to a metal again. The semimetal and metal are distinguished by the density of states at the Weyl point, as well as by level statistics, transport, and the momentum-space structure of eigenstates near the Weyl point. The critical properties of the transitions in quasiperiodic systems differ from those in random systems: we do not find a clear critical scaling regime in energy; instead, at the quasiperiodic transitions, the density of states appears to jump abruptly (and discontinuously to within our resolution).

preprint2016arXiv

A Floquet Model for the Many-Body Localization Transition

The nature of the dynamical quantum phase transition between the many-body localized (MBL) phase and the thermal phase remains an open question, and one line of attack on this problem is to explore this transition numerically in finite-size systems. To maximize the contrast between the MBL phase and the thermal phase in such finite-size systems, we argue one should choose a Floquet model with no local conservation laws and rapid thermalization to "infinite temperature" in the thermal phase. Here we introduce and explore such a Floquet spin chain model, and show that standard diagnostics of the MBL-to-thermal transition behave well in this model even at modest sizes. We also introduce a physically motivated spacetime correlation function which peaks at the transition in the Floquet model, but is strongly affected by conservation laws in Hamiltonian models.

preprint2016arXiv

Exploring the many-body localization transition in two dimensions

One fundamental assumption in statistical physics is that generic closed quantum many-body systems thermalize under their own dynamics. Recently, the emergence of many-body localized systems has questioned this concept, challenging our understanding of the connection between statistical physics and quantum mechanics. Here we report on the observation of a many-body localization transition between thermal and localized phases for bosons in a two-dimensional disordered optical lattice. With our single site resolved measurements we track the relaxation dynamics of an initially prepared out-of-equilibrium density pattern and find strong evidence for a diverging length scale when approaching the localization transition. Our experiments mark the first demonstration and in-depth characterization of many-body localization in a regime not accessible with state-of-the-art simulations on classical computers.

preprint2016arXiv

Griffiths effects and slow dynamics in nearly many-body localized systems

The low-frequency response of systems near a many-body localization transition can be dominated by rare regions that are locally critical or "in the other phase". It is known that, in one dimension, these rare regions can cause the d.c. conductivity and diffusion constant to vanish even inside the delocalized thermal phase. Here, we present a general analysis of such Griffiths effects in the thermal phase near the many-body localization transition: we consider both one-dimensional and higher-dimensional systems, subject to quenched randomness, and discuss both linear response (including the frequency- and wavevector-dependent conductivity) and more general dynamics. In all the regimes we consider, we identify observables that are dominated by rare-region effects. In some cases (one-dimensional systems and Floquet systems with no extensive conserved quantities), essentially all long-time local observables are dominated by rare-region effects; in others, generic observables are instead dominated by hydrodynamic long-time tails throughout the thermal phase, and one must look at specific probes, such as spin echo, to see Griffiths behavior.

preprint2016arXiv

Heating and many-body resonances in a periodically driven two-band system

We study the dynamics and stability in a strongly interacting resonantly driven two-band model. Using exact numerical simulations, we find a stable regime at large driving frequencies where the time evolution is governed by a local Floquet Hamiltonian that is approximately conserved out to very long times. For slow driving, on the other hand, the system becomes unstable and heats up to infinite temperature. While thermalization is relatively fast in these two regimes (but to different "temperatures"), in the crossover between them we find slow nonthermalizing time evolution: temporal fluctuations become strong and temporal correlations long lived. Microscopically, we trace back the origin of this nonthermalizing time evolution to the properties of rare Floquet many-body resonances, whose proliferation at lower driving frequency removes the approximate energy conservation, and thus produces thermalization to infinite temperature.

preprint2016arXiv

Many-body localization phase transition: A simplified strong-randomness approximate renormalization group

We present a simplified strong-randomness renormalization group (RG) that captures some aspects of the many-body localization (MBL) phase transition in generic disordered one-dimensional systems. This RG can be formulated analytically, and the critical fixed point distribution and critical exponents (that satisfy the Chayes inequality) are obtained to numerical precision by solving integro-differential equations. This reproduces many, but not all, of the qualitative features of the MBL phase transition that are suggested by previous numerical work and approximate RG studies: our RG might serve as a "zeroth-order" approximation for future RG studies. One interesting feature that we highlight is that the rare Griffiths regions are fractal. For thermal Griffiths regions within the MBL phase, this feature might be qualitatively correctly captured by our RG. If this is correct beyond our approximations, then these Griffiths effects are stronger than has been previously assumed.

preprint2016arXiv

Rare region induced avoided quantum criticality in disordered three-dimensional Dirac and Weyl semimetals

We numerically study the effect of short ranged potential disorder on massless noninteracting three-dimensional Dirac and Weyl fermions, with a focus on the question of the proposed quantum critical point separating the semimetal and diffusive metal phases. We determine the properties of the eigenstates of the disordered Dirac Hamiltonian ($H$) and exactly calculate the density of states (DOS) near zero energy, using a combination of Lanczos on $H^2$ and the kernel polynomial method on $H$. We establish the existence of two distinct types of low energy eigenstates contributing to the disordered density of states in the weak disorder semimetal regime. These are (i) typical eigenstates that are well described by linearly dispersing perturbatively dressed Dirac states, and (ii) nonperturbative rare eigenstates that are weakly-dispersive and quasi-localized in the real space regions with the largest (and rarest) local random potential. Using twisted boundary conditions, we are able to systematically find and study these two types of eigenstates. We find that the Dirac states contribute low energy peaks in the finite-size DOS that arise from the clean eigenstates which shift and broaden in the presence of disorder. On the other hand, we establish that the rare quasi-localized eigenstates contribute a nonzero background DOS which is only weakly energy-dependent near zero energy and is exponentially small at weak disorder. We find that the expected semimetal to diffusive metal quantum critical point is converted to an {\it avoided} quantum criticality that is "rounded out" by nonperturbative effects, with no signs of any singular behavior in the DOS near the Dirac energy. We discuss the implications of our results for disordered Dirac and Weyl semimetals, and reconcile the large body of existing numerical work showing quantum criticality with the existence of the rare region effects.

preprint2016arXiv

Uncovering the hidden quantum critical point in disordered massless Dirac and Weyl semi-metals

We study the properties of the avoided or hidden quantum critical point (AQCP) in three dimensional Dirac and Weyl semi-metals in the presence of short range potential disorder. By computing the averaged density of states (along with its second and fourth derivative at zero energy) with the kernel polynomial method (KPM) we systematically tune the effective length scale that eventually rounds out the transition and leads to an AQCP. We show how to determine the strength of the avoidance, establishing that it is not controlled by the long wavelength component of the disorder. Instead, the amount of avoidance can be adjusted via the tails of the probability distribution of the local random potentials. A binary distribution with no tails produces much less avoidance than a Gaussian distribution. We introduce a double Gaussian distribution to interpolate between these two limits. As a result we are able to make the length scale of the avoidance sufficiently large so that we can accurately study the properties of the underlying transition (that is eventually rounded out), unambiguously identify its location, and provide accurate estimates of the critical exponents $ν=1.01\pm0.06$ and $z=1.50\pm0.04$. We also show that the KPM expansion order introduces an effective length scale that can also round out the transition in the scaling regime near the AQCP.

preprint2015arXiv

Fluctuating hydrodynamics for a discrete Gross-Pitaevskii equation: mapping to Kardar-Parisi-Zhang universality class

We show that several aspects of the low-temperature hydrodynamics of a discrete Gross-Pitaevskii equation (GPE) can be understood by mapping it to a nonlinear version of fluctuating hydrodynamics. This is achieved by first writing the GPE in a hydrodynamic form of a continuity and an Euler equation. Respecting conservation laws, dissipation and noise due to the system's chaos are added, thus giving us a nonlinear stochastic field theory in general and the Kardar-Parisi-Zhang (KPZ) equation in our particular case. This mapping to KPZ is benchmarked against exact Hamiltonian numerics on discrete GPE by investigating the non-zero temperature dynamical structure factor and its scaling form and exponent. Given the ubiquity of the Gross-Pitaevskii equation (a.k.a. nonlinear Schrodinger equation), ranging from nonlinear optics to cold gases, we expect this remarkable mapping to the KPZ equation to be of paramount importance and far reaching consequences.

preprint2015arXiv

Low-frequency conductivity in many-body localized systems

We argue that the a.c. conductivity $σ(ω)$ in the many-body localized phase is a power law of frequency $ω$ at low frequency: specifically, $σ(ω) \sim ω^α$ with the exponent $α$ approaching 1 at the phase transition to the thermal phase, and asymptoting to 2 deep in the localized phase. We identify two separate mechanisms giving rise to this power law: deep in the localized phase, the conductivity is dominated by rare resonant pairs of configurations; close to the transition, the dominant contributions are rare regions that are locally critical or in the thermal phase. We present numerical evidence supporting these claims, and discuss how these power laws can also be seen through polarization-decay measurements in ultracold atomic systems.

preprint2015arXiv

Many-body localization in a quantum simulator with programmable random disorder

When a system thermalizes it loses all local memory of its initial conditions. This is a general feature of open systems and is well described by equilibrium statistical mechanics. Even within a closed (or reversible) quantum system, where unitary time evolution retains all information about its initial state, subsystems can still thermalize using the rest of the system as an effective heat bath. Exceptions to quantum thermalization have been predicted and observed, but typically require inherent symmetries or noninteracting particles in the presence of static disorder. The prediction of many-body localization (MBL), in which disordered quantum systems can fail to thermalize in spite of strong interactions and high excitation energy, was therefore surprising and has attracted considerable theoretical attention. Here we experimentally generate MBL states by applying an Ising Hamiltonian with long-range interactions and programmably random disorder to ten spins initialized far from equilibrium. We observe the essential signatures of MBL: memory retention of the initial state, a Poissonian distribution of energy level spacings, and entanglement growth in the system at long times. Our platform can be scaled to higher numbers of spins, where detailed modeling of MBL becomes impossible due to the complexity of representing such entangled quantum states. Moreover, the high degree of control in our experiment may guide the use of MBL states as potential quantum memories in naturally disordered quantum systems.

preprint2015arXiv

Observation of antiferromagnetic correlations in the Hubbard model with ultracold atoms

Ultracold atoms in optical lattices have great potential to contribute to a better understanding of some of the most important issues in many-body physics, such as high-$T_c$ superconductivity. The Hubbard model describes many of the features shared by the copper oxides, including an interaction-driven Mott insulating state and an antiferromagnetic (AFM) state. Optical lattices filled with a two-spin-component Fermi gas of ultracold atoms can faithfully realise the Hubbard model with readily tunable parameters, and thus provide a platform for the systematic exploration of its phase diagram. Realisation of strongly correlated phases, however, has been hindered by the need to cool the atoms to temperatures as low as the magnetic exchange energy, and also by the lack of reliable thermometry. Here we demonstrate spin-sensitive Bragg scattering of light to measure AFM spin correlations in a realisation of the 3D Hubbard model at temperatures down to 1.4 times that of the AFM phase transition. This temperature regime is beyond the range of validity of a simple high-temperature series expansion, which brings our experiment close to the limit of the capabilities of current numerical techniques. We reach these low temperatures using a unique compensated optical lattice technique, in which the confinement of each lattice beam is compensated by a blue-detuned laser beam. The temperature of the atoms in the lattice is deduced by comparing the light scattering to determinantal quantum Monte Carlo and numerical linked-cluster expansion calculations. Further refinement of the compensated lattice may produce even lower temperatures which, along with light scattering thermometry, would open avenues for achieving and characterising other novel quantum states of matter, such as the pseudogap regime of the 2D Hubbard model.

preprint2015arXiv

Slowest local operators in quantum spin chains

We numerically construct slowly relaxing local operators in a nonintegrable spin-1/2 chain. Restricting the support of the operator to $M$ consecutive spins along the chain, we exhaustively search for the operator that minimizes the Frobenius norm of the commutator with the Hamiltonian. We first show that the Frobenius norm bounds the time scale of relaxation of the operator at high temperatures. We find operators with significantly slower relaxation than the slowest simple "hydrodynamic" mode due to energy diffusion. Then, we examine some properties of the nontrivial slow operators. Using both exhaustive search and tensor network techniques, we find similar slowly relaxing operators for a Floquet spin chain; this system is hydrodynamically "trivial", with no conservation laws restricting their dynamics. We argue that such slow relaxation may be a generic feature following from locality and unitarity.

preprint2015arXiv

Theory of the many-body localization transition in one dimensional systems

We formulate a theory of the many-body localization transition based on a novel real space renormalization group (RG) approach. The results of this theory are corroborated and intuitively explained with a phenomenological effective description of the critical point and of the "badly conducting" state found near the critical point on the delocalized side. The theory leads to the following sharp predictions: (i) The delocalized state established near the transition is a Griffiths phase, which exhibits sub-diffusive transport of conserved quantities and sub-ballistic spreading of entanglement. The anomalous diffusion exponent $α< 1/2$ vanishes continuously at the critical point. The system does thermalize in this Griffiths phase. (ii) The many-body localization transition is controlled by a new kind of infinite randomness RG fixed point, where the broadly distributed scaling variable is closely related to the eigenstate entanglement entropy. Dynamically, the entanglement grows as $\sim\log t$ at the critical point, as it also does in the localized phase. (iii) In the vicinity of the critical point the ratio of the entanglement entropy to the thermal entropy, and its variance (and in fact all moments) are scaling functions of $L/ξ$, where $L$ is the length of the system and $ξ$ is the correlation length, which has a power-law divergence at the critical point.

preprint2015arXiv

Thermal equilibrium of a macroscopic quantum system in a pure state

We consider the notion of thermal equilibrium for an individual closed macroscopic quantum system in a pure state, i.e., described by a wave function. The macroscopic properties in thermal equilibrium of such a system, determined by its wave function, must be the same as those obtained from thermodynamics, e.g., spatial uniformity of temperature and chemical potential. When this is true we say that the system is in macroscopic thermal equilibrium (MATE). Such a system may however not be in microscopic thermal equilibrium (MITE). The latter requires that the reduced density matrices of small subsystems be close to those obtained from the microcanonical, equivalently the canonical, ensemble for the whole system. The distinction between MITE and MATE is particularly relevant for systems with many-body localization (MBL) for which the energy eigenfunctions fail to be in MITE while necessarily most of them, but not all, are in MATE. We note however that for generic macroscopic systems, including those with MBL, most wave functions in an energy shell are in both MATE and MITE. For a classical macroscopic system, MATE holds for most phase points on the energy surface, but MITE fails to hold for any phase point.

preprint2015arXiv

Thermalization of entanglement

We explore the dynamics of the entanglement entropy near equilibrium in highly-entangled pure states of two quantum-chaotic spin chains undergoing unitary time evolution. We examine the relaxation to equilibrium from initial states with either less or more entanglement entropy than the equilibrium value, as well as the dynamics of the spontaneous fluctuations of the entanglement that occur in equilibrium. For the spin chain with a time-independent Hamiltonian and thus an extensive conserved energy, we find slow relaxation of the entanglement entropy near equilibration. Such slow relaxation is absent in a Floquet spin chain with a Hamiltonian that is periodic in time and thus has no local conservation law. Therefore, we argue that slow diffusive energy transport is responsible for the slow relaxation of the entanglement entropy in the Hamiltonian system.

preprint2014arXiv

Entanglement spreading in a many-body localized system

Motivated by the findings of logarithmic spreading of entanglement in a many-body localized system, we more closely examine the spreading of entanglement in the fully many-body localized phase, where all many-body eigenstates are localized. Performing full diagonalizations of an XXZ spin model with random longitudinal fields, we identify two factors contributing to the spreading rate: the localization length ($ξ$), which depends on the disorder strength, and the final value of entanglement per spin ($s_\infty$), which primarily depends on the initial state. We find that the entanglement entropy grows with time as $\sim ξ\times s_\infty \log t$, providing support for the phenomenology of many-body localized systems recently proposed by Huse and Oganesyan [arXiv:1305.4915v1].

preprint2014arXiv

Many body localization and thermalization in quantum statistical mechanics

We review some recent developments in the statistical mechanics of isolated quantum systems. We provide a brief introduction to quantum thermalization, paying particular attention to the `Eigenstate Thermalization Hypothesis' (ETH), and the resulting `single-eigenstate statistical mechanics'. We then focus on a class of systems which fail to quantum thermalize and whose eigenstates violate the ETH: These are the many-body Anderson localized systems; their long-time properties are not captured by the conventional ensembles of quantum statistical mechanics. These systems can locally remember forever information about their local initial conditions, and are thus of interest for possibilities of storing quantum information. We discuss key features of many-body localization (MBL), and review a phenomenology of the MBL phase. Single-eigenstate statistical mechanics within the MBL phase reveals dynamically-stable ordered phases, and phase transitions among them, that are invisible to equilibrium statistical mechanics and can occur at high energy and low spatial dimensionality where equilibrium ordering is forbidden.

preprint2014arXiv

Spectral features of a many-body localized system weakly coupled to a heat bath

We study many-body-localized (MBL) systems that are weakly coupled to thermalizing environments, focusing on the spectral functions of local operators. We argue that these spectral functions carry signatures of localization even away from the limit of perfectly isolated systems. We find that, in the limit of vanishing coupling to a bath, MBL systems come in two varieties, with either discrete or continuous local spectra. Both varieties of MBL systems exhibit a "soft gap" at zero frequency in the spatially-averaged spectral functions of local operators, which serves as a diagnostic for localization. We estimate the degree to which coupling to a bath broadens these spectral features, and find that characteristics of incipient localization survive as long as the system-bath coupling is much weaker than the characteristic energy scales of the system. Since perfect isolation is impossible, we expect the ideas discussed in this paper to be relevant for all experiments on many-body localization.

preprint2014arXiv

Testing whether all eigenstates obey the Eigenstate Thermalization Hypothesis

We ask whether the Eigenstate Thermalization Hypothesis (ETH) is valid in a strong sense: in the limit of an infinite system, {\it every} eigenstate is thermal. We examine expectation values of few-body operators in highly-excited many-body eigenstates and search for `outliers', the eigenstates that deviate the most from ETH. We use exact diagonalization of two one-dimensional nonintegrable models: a quantum Ising chain with transverse and longitudinal fields, and hard-core bosons at half-filling with nearest- and next-nearest-neighbor hopping and interaction. We show that even the most extreme outliers appear to obey ETH as the system size increases, and thus provide numerical evidences that support ETH in this strong sense. Finally, periodically driving the Ising Hamiltonian, we show that the eigenstates of the corresponding Floquet operator obey ETH even more closely. We attribute this better thermalization to removing the constraint of conservation of the total energy.

preprint2014arXiv

Three `species' of Schrödinger cat states in an infinite-range spin model

We explore a transverse-field Ising model that exhibits both spontaneous symmetry-breaking and eigenstate thermalization. Within its ferromagnetic phase, the exact eigenstates of the Hamiltonian of any large but finite-sized system are all Schrödinger cat states: superpositions of states with `up' and `down' spontaneous magnetization. This model exhibits two dynamical phase transitions {\it within} its ferromagnetic phase: In the lowest-temperature phase the magnetization can macroscopically oscillate between up and down. The relaxation of the magnetization is always overdamped in the remainder of the ferromagnetic phase, which is divided in to phases where the system thermally activates itself {\it over} the barrier between the up and down states, and where it quantum tunnels.

preprint2014arXiv

Unexpected z-Direction Ising Antiferromagnetic Order in a frustrated Spin-1/2 $J_1-J_2$ XY Model on the Honeycomb Lattice

Using the density matrix renormalization group (DMRG) on wide cylinders, we study the phase diagram of the spin-1/2 XY model on the honeycomb lattice, with first-neighbor ($J_1 = 1$) and frustrating second-neighbor ($J_2>0$) interactions. For the intermediate frustration regime $0.22\lesssim J_2\lesssim0.36$, we find a surprising antiferromagnetic Ising phase, with ordered moments pointing along the z axis, despite the absence of any S_z_z interactions in the Hamiltonian. Surrounding this phase as a function of $J_2$ are antiferromagnetic phases with the moments pointing in the $x-y$ plane for small $J_2$ and a close competition between an $x-y$ plane magnetic collinear phase and a dimer phase for large values of $J_2$. We do not find any spin liquid phases in this model.

preprint2013arXiv

A phenomenology of certain many-body-localized systems

We consider isolated quantum systems with all of their many-body eigenstates localized. We define a sense in which such systems are integrable, and discuss a method for finding their localized conserved quantum numbers ("constants of motion"). These localized operators are interacting pseudospins and are subject to dephasing but not to dissipation, so any quantum states of these pseudospins can in principle be recovered via (spin) echo procedures. We also discuss the spreading of entanglement in many-body localized systems, which is another aspect of the dephasing due to interactions between these localized conserved operators.

preprint2013arXiv

Anyonic quantum spin chains: Spin-1 generalizations and topological stability

There are many interesting parallels between systems of interacting non-Abelian anyons and quantum magnetism, occuring in ordinary SU(2) quantum magnets. Here we consider theories of so-called su(2)_k anyons, well-known deformations of SU(2), in which only the first k+1 angular momenta of SU(2) occur. In this manuscript, we discuss in particular anyonic generalizations of ordinary SU(2) spin chains with an emphasis on anyonic spin S=1 chains. We find that the overall phase diagrams for these anyonic spin-1 chains closely mirror the phase diagram of the ordinary bilinear-biquadratic spin-1 chain including anyonic generalizations of the Haldane phase, the AKLT construction, and supersymmetric quantum critical points. A novel feature of the anyonic spin-1 chains is an additional topological symmetry that protects the gapless phases. Distinctions further arise in the form of an even/odd effect in the deformation parameter k when considering su(2)_k anyonic theories with k>4, as well as for the special case of the su(2)_4 theory for which the spin-1 representation plays a special role. We also address anyonic generalizations of spin-1/2 chains with a focus on the topological protection provided for their gapless ground states. Finally, we put our results into context of earlier generalizations of SU(2) quantum spin chains, in particular so-called (fused) Temperley-Lieb chains.

preprint2013arXiv

Ballistic spreading of entanglement in a diffusive nonintegrable system

We study the time evolution of the entanglement entropy of a one-dimensional nonintegrable spin chain, starting from random nonentangled initial pure states. We use exact diagonalization of a nonintegrable quantum Ising chain with transverse and longitudinal fields to obtain the exact quantum dynamics. We show that the entanglement entropy increases linearly with time before finite-size saturation begins, demonstrating a ballistic spreading of the entanglement, while the energy transport in the same system is diffusive. Thus we explicitly demonstrate that the spreading of entanglement is much faster than the energy diffusion in this nonintegrable system.

preprint2013arXiv

Dirty Weyl fermions: rare region effects near 3D Dirac points

We study three-dimensional Dirac fermions with weak finite-range scalar potential disorder. We show that even though disorder is perturbatively irrelevant at 3D Dirac points, nonperturbative effects from rare regions give rise to a nonzero density of states and a finite mean free path, with the transport at the Dirac point being dominated by hopping between rare regions. As one moves in chemical potential away from the Dirac point, there are interesting intermediate-energy regimes where the rare regions produce scattering resonances that determine the DC conductivity. We also discuss the interplay of disorder with interactions at the Dirac point. Attractive interactions drive a transition into a granular superconductor, with a critical temperature that depends strongly on the disorder distribution. In the presence of Coulomb repulsion and weak retarded attraction, the system can be a Bose glass. Our results apply to all 3D systems with Dirac points, including Weyl semimetals, and overturn a thirty year old consensus regarding the irrelevance of weak disorder at 3D Dirac points.

preprint2013arXiv

Localization protected quantum order

Closed quantum systems with quenched randomness exhibit many-body localized regimes wherein they do not equilibrate even though prepared with macroscopic amounts of energy above their ground states. We show that such localized systems can order in that individual many-body eigenstates can break symmetries or display topological order in the infinite volume limit. Indeed, isolated localized quantum systems can order even at energy densities where the corresponding thermally equilibrated system is disordered, i.e.: localization protects order. In addition, localized systems can move between ordered and disordered localized phases via non-thermodynamic transitions in the properties of the many-body eigenstates. We give evidence that such transitions may proceed via localized critical points. We note that localization provides protection against decoherence that may allow experimental manipulation of macroscopic quantum states. We also identify a `spectral transition' involving a sharp change in the spectral statistics of the many-body Hamiltonian.

preprint2013arXiv

Many-Body Localization in a Quasiperiodic System

Recent theoretical and numerical evidence suggests that localization can survive in disordered many-body systems with very high energy density, provided that interactions are sufficiently weak. Stronger interactions can destroy localization, leading to a so-called many-body localization transition. This dynamical phase transition is relevant to questions of thermalization in extended quantum systems far from the zero-temperature limit. It separates a many-body localized phase, in which localization prevents transport and thermalization, from a conducting ("ergodic") phase in which the usual assumptions of quantum statistical mechanics hold. Here, we present numerical evidence that many-body localization also occurs in models without disorder but rather a quasiperiodic potential. In one dimension, these systems already have a single-particle localization transition, and we show that this transition becomes a many-body localization transition upon the introduction of interactions. We also comment on possible relevance of our results to experimental studies of many-body dynamics of cold atoms and non-linear light in quasiperiodic potentials.

preprint2013arXiv

Superconductivity of disordered Dirac fermions

We study the effect of disorder on massless, spinful Dirac fermions in two spatial dimensions with attractive interactions, and show that the combination of disorder and attractive interactions is deadly to the Dirac semimetal phase. First, we derive the zero temperature phase diagram of a clean Dirac fermion system with tunable doping level (μ) and attraction strength (g). We show that it contains two phases: a superconductor and a Dirac semimetal. Then, we show that arbitrarily weak disorder destroys the Dirac semimetal, turning it into a superconductor. We discuss the strength of the superconductivity for both long range and short range disorder. For long range disorder, the superconductivity is exponentially weak in the disorder strength. For short range disorder, a uniform mean field analysis predicts that superconductivity should be doubly exponentially weak in the disorder strength. However, a more careful treatment of mesoscopic fluctuations suggests that locally superconducting puddles should form at a much higher temperature, and should establish global phase coherence at a temperature that is only exponentially small in weak disorder. We also discuss the effect of disorder on the quantum critical point of the clean system, building in the effect of disorder through a replica field theory. We show that disorder is a relevant perturbation to the supersymmetric quantum critical point. We expect that in the presence of attractive interactions, the flow away from the critical point ends up in the superconducting phase, although firm conclusions cannot be drawn since the renormalization group analysis flows to strong coupling. We argue that although we expect the quantum critical point to get buried under a superconducting phase, signatures of the critical point may be visible in the finite temperature quantum critical regime.

preprint2013arXiv

Weak plaquette valence bond order in the $S=1/2$ honeycomb $J_1-J_2$ Heisenberg model

Using the density matrix renormalization group, we investigate the $S=1/2$ Heisenberg model on the honeycomb lattice with first- ($J_1$) and second-neighbor ($J_2$) interactions. We are able to study long open cylinders with widths up to 12 lattice spacings. For $J_2/J_1$ near 0.3, we find an apparently paramagnetic phase, bordered by an antiferromagnetic phase for $J_2\lesssim 0.26$ and by a valence bond crystal for $J_2\gtrsim 0.36$. The longest correlation length that we find in this intermediate phase is for plaquette valence bond (PVB) order. This correlation length grows strongly with cylinder circumference, indicating either quantum criticality or weak PVB order.

preprint2012arXiv

Enlarging and cooling the Néel state in an optical lattice

We propose an experimental scheme to favor both the realization and the detection of the Néel state in a two-component gas of ultracold fermions in a three-dimensional simple-cubic optical lattice. By adding three compensating Gaussian laser beams to the standard three pairs of retroreflected lattice beams, and adjusting the relative waists and intensities of the beams, one can significantly enhance the size of the Néel state in the trap, thus increasing the signal of optical Bragg scattering. Furthermore, the additional beams provide for adjustment of the local chemical potential and the possibility to evaporatively cool the gas while in the lattice. Our proposals are relevant to other attempts to realize many-body quantum phases in optical lattices.

preprint2012arXiv

Heat and spin transport in a cold atomic Fermi gas

Motivated by recent experiments measuring the spin transport in ultracold unitary atomic Fermi gases (Sommer et al., 2011; Sommer et al., 2011), we explore the theory of spin and heat transport in a three-dimensional spin-polarized atomic Fermi gas. We develop estimates of spin and thermal diffusivities and discuss magnetocaloric effects, namely the the spin Seebeck and spin Peltier effects. We estimate these transport coefficients using a Boltzmann kinetic equation in the classical regime and present experimentally accessible signatures of the spin Seebeck effect. We study an exactly solvable model that illustrates the role of momentum-dependent scattering in the magnetocaloric effects.

preprint2012arXiv

Non-equilibrium dynamic critical scaling of the quantum Ising chain

We solve for the time-dependent finite-size scaling functions of the 1D transverse-field Ising chain during a linear-in-time ramp of the field through the quantum critical point. We then simulate Mott-insulating bosons in a tilted potential, an experimentally-studied system in the same equilibrium universality class, and demonstrate that universality holds for the dynamics as well. We find qualitatively athermal features of the scaling functions, such as negative spin correlations, and show that they should be robustly observable within present cold atom experiments.

preprint2012arXiv

Superdiffusive nonequilibrium motion of an impurity in a Fermi sea

We treat the nonequilibrium motion of a single impurity atom in a low-temperature single-species Fermi sea, interacting via a contact interaction. In the nonequilibrium regime, the impurity does a superdiffusive geometric random walk where the typical distance traveled grows with time as $\sim t^{d/(d+1)}$ for the $d$-dimensional system with $d\geq 2$. For nonzero temperature $T$, this crosses over to diffusive motion at long times with diffusivity $D\sim T^{-(d-1)/2}$. These results apply also to a nonzero concentration of impurity atoms as long as they remain dilute and nondegenerate.

preprint2011arXiv

Disorder-induced Majorana metal in interacting non-Abelian anyon systems

We demonstrate that a thermal metal of Majorana fermions forms in a two-dimensional system of interacting non-Abelian (Ising) anyons in the presence of moderate disorder. This bulk metallic phase arises in the $ν=5/2$ quantum Hall state when disorder pins the anyonic quasiparticles. More generally, it naturally occurs for various proposed systems supporting Majorana fermion zero modes when disorder induces the random pinning of a finite density of vortices. This includes all two-dimensional topological superconductors in so-called symmetry class D. A distinct experimental signature of the thermal metal phase is the presence of bulk heat transport down to zero temperature.

preprint2010arXiv

Chaos and universality in two-dimensional Ising spin glasses

Recently extended precise numerical methods and droplet scaling arguments allow for a coherent picture of the glassy states of two-dimensional Ising spin glasses to be assembled. The length scale at which entropy becomes important and produces "chaos", the extreme sensitivity of the state to temperature, is found to depend on the type of randomness. For the $\pm J$ model this length scale dominates the low-temperature specific heat. Although there is a type of universality, some critical exponents do depend on the distribution of disorder.

preprint2010arXiv

Diagnosing Deconfinement and Topological Order

Topological or deconfined phases are characterized by emergent, weakly fluctuating, gauge fields. In condensed matter settings they inevitably come coupled to excitations that carry the corresponding gauge charges which invalidate the standard diagnostic of deconfinement---the Wilson loop. Inspired by a mapping between symmetric sponges and the deconfined phase of the $Z_2$ gauge theory, we construct a diagnostic for deconfinement that has the interpretation of a line tension. One operator version of this diagnostic turns out to be the Fredenhagen-Marcu order parameter known to lattice gauge theorists and we show that a different version is best suited to condensed matter systems. We discuss generalizations of the diagnostic, use it to establish the existence of finite temperature topological phases in $d \ge 3$ dimensions and show that multiplets of the diagnostic are useful in settings with multiple phases such as $U(1)$ gauge theories with charge $q$ matter. [Additionally we present an exact reduction of the partition function of the toric code in general dimensions to a well studied problem.]

preprint2010arXiv

Interfering directed paths and the sign phase transition

We revisit the question of the "sign phase transition" for interfering directed paths with real amplitudes in a random medium. The sign of the total amplitude of the paths to a given point may be viewed as an Ising order parameter, so we suggest that a coarse-grained theory for system is a dynamic Ising model coupled to a Kardar-Parisi-Zhang (KPZ) model. It appears that when the KPZ model is in its strong-coupling ("pinned") phase, the Ising model does not have a stable ferromagnetic phase, so there is no sign phase transition. We investigate this numerically for the case of ß1+1 dimensions, demonstrating the instability of the Ising ordered phase there.

preprint2010arXiv

Novel Ground-State Crystals with Controlled Vacancy Concentrations: From Kagomé to Honeycomb to Stripes

We introduce a one-parameter family, $0 \leq H \leq 1$, of pair potential functions with a single relative energy minimum that stabilize a range of vacancy-riddled crystals as ground states. The "quintic potential" is a short-ranged, nonnegative pair potential with a single local minimum of height $H$ at unit distance and vanishes cubically at a distance of $\rt$. We have developed this potential to produce ground states with the symmetry of the triangular lattice while favoring the presence of vacancies. After an exhaustive search using various optimization and simulation methods, we believe that we have determined the ground states for all pressures, densities, and $0 \leq H \leq 1$. For specific areas below $3\rt/2$, the ground states of the "quintic potential" include high-density and low-density triangular lattices, kagomé and honeycomb crystals, and stripes. We find that these ground states are mechanically stable but are difficult to self-assemble in computer simulations without defects. For specific areas above $3\rt/2$, these systems have a ground-state phase diagram that corresponds to hard disks with radius $\rt$. For the special case of H=0, a broad range of ground states is available. Analysis of this case suggests that among many ground states, a high-density triangular lattice, low-density triangular lattice, and striped phases have the highest entropy for certain densities. The simplicity of this potential makes it an attractive candidate for experimental realization with application to the development of novel colloidal crystals or photonic materials.

preprint2010arXiv

Optimized broad-histogram simulations for strong first-order phase transitions: Droplet transitions in the large-Q Potts model

The numerical simulation of strongly first-order phase transitions has remained a notoriously difficult problem even for classical systems due to the exponentially suppressed (thermal) equilibration in the vicinity of such a transition. In the absence of efficient update techniques, a common approach to improve equilibration in Monte Carlo simulations is to broaden the sampled statistical ensemble beyond the bimodal distribution of the canonical ensemble. Here we show how a recently developed feedback algorithm can systematically optimize such broad-histogram ensembles and significantly speed up equilibration in comparison with other extended ensemble techniques such as flat-histogram, multicanonical or Wang-Landau sampling. As a prototypical example of a strong first-order transition we simulate the two-dimensional Potts model with up to Q=250 different states on large systems. The optimized histogram develops a distinct multipeak structure, thereby resolving entropic barriers and their associated phase transitions in the phase coexistence region such as droplet nucleation and annihilation or droplet-strip transitions for systems with periodic boundary conditions. We characterize the efficiency of the optimized histogram sampling by measuring round-trip times tau(N,Q) across the phase transition for samples of size N spins. While we find power-law scaling of tau vs. N for small Q \lesssim 50 and N \lesssim 40^2, we observe a crossover to exponential scaling for larger Q. These results demonstrate that despite the ensemble optimization broad-histogram simulations cannot fully eliminate the supercritical slowing down at strongly first-order transitions.

preprint2010arXiv

Spin Liquid Ground State of the $S=1/2$ Kagome Heisenberg Model

Condensed matter physicists have long sought a realistic two-dimensional (2D) magnetic system whose ground state is a {\it spin liquid}---a zero temperature state in which quantum fluctuations have melted away any form of magnetic order. The nearest-neighbor $S=1/2$ Heisenberg model on the kagome lattice has seemed an ideal candidate, but in recent years some approximate numerical approaches to it have yielded instead a valence bond crystal. We have used the density matrix renormalization group to perform very accurate simulations on numerous cylinders with circumferences up to 12 lattice spacings, finding instead of the valence bond crystal a singlet-gapped spin liquid with substantially lower energy that appears to have $Z_2$ topological order. Our results, through a combination of very low energy, short correlation lengths and corresponding small finite size effects, a new rigorous energy bound, and consistent behavior on many cylinders, provide strong evidence that the 2D ground state of this model is a gapped spin liquid.

preprint2010arXiv

The many-body localization phase transition

We use exact diagonalization to explore the many-body localization transition in a random-field spin-1/2 chain. We examine the correlations within each many-body eigenstate, looking at all high-energy states and thus effectively working at infinite temperature. For weak random field the eigenstates are thermal, as expected in this nonlocalized, "ergodic" phase. For strong random field the eigenstates are localized, with only short-range entanglement. We roughly locate the localization transition and examine some of its finite-size scaling, finding that this quantum phase transition at nonzero temperature might be showing infinite-randomness scaling with a dynamic critical exponent $z\rightarrow\infty$.

preprint2010arXiv

The many-body localization transition

We use exact diagonalization to explore the many-body localization transition in a random-field spin-1/2 chain. We examine the correlations within each many-body eigenstate, looking at all states and thus effectively working at infinite temperature. For weak random field the eigenstates are thermal, as expected in this nonlocalized, ergodic phase. For strong random field the eigenstates are localized, with only short-range entanglement. We roughly locate the localization transition and examine the probability distributions of the correlations, finding that this quantum phase transition at nonzero temperature may be showing infinite-randomness scaling.

preprint2010arXiv

Trimers, molecules and polarons in imbalanced atomic Fermi gases

We consider the ground state of a single "spin-down" impurity atom interacting attractively with a "spin-up" atomic Fermi gas. By constructing variational wave functions for polarons, molecules and trimers, we perform a detailed study of the transitions between each of these dressed bound states as a function of mass ratio $r=m_\uparrow/m_\downarrow$ and interaction strength. We find that the presence of a Fermi sea enhances the stability of the $p$-wave trimer, which can be viewed as a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) molecule that has bound an additional majority atom. For sufficiently large $r$, we find that the transitions lie outside the region of phase separation in imbalanced Fermi gases and should thus be observable in experiment, unlike the well-studied equal-mass case.

preprint2010arXiv

Visualizing Critical Correlations near the Metal-Insulator Transition in Ga1-xMnxAs

Electronic states in disordered conductors on the verge of localization are predicted to exhibit critical spatial characteristics indicative of the proximity to a metal- insulator phase transition. We have used scanning tunneling microscopy to visualize electronic states in Ga1-xMnxAs samples close to this transition. Our measurements show that doping-induced disorder produces strong spatial variations in the local tunneling conductance across a wide range of energies. Near the Fermi energy, where spectroscopic signatures of electron-electron interaction are the most prominent, the electronic states exhibit a diverging spatial correlation length. Power-law decay of the spatial correlations is accompanied by log-normal distributions of the local density of states and multifractal spatial characteristics. Our method can be used to explore critical correlations in other materials close to a quantum critical point.

preprint2009arXiv

Energy transport in disordered classical spin chains

We present a numerical study of the diffusion of energy at high temperature in strongly disordered chains of interacting classical spins evolving deterministically. We find that quenched randomness strongly suppresses transport, with the diffusion constant becoming reduced by several orders of magnitude upon the introduction of moderate disorder. We have also looked for but not found signs of a classical many-body localization transition at any nonzero strength of the spin-spin interactions.

preprint2008arXiv

Maximizing the Neel temperature of fermions in a simple-cubic optical lattice

For a simple-cubic optical lattice with lattice spacing d, occupied by two species of fermionic atoms of mass m that interact repulsively, we ask what conditions maximize the Neel temperature in the Mott insulating phase at density one atom per site, with equal numbers of the two species. This maximum occurs near the edge of the regime where the system is well-approximated by the usual Hubbard model. The correction to the Hubbard-model approximation that produces a "direct" ferromagnetic interaction between atoms in nearest-neighbor Wannier orbitals is the leading term that limits how high the Neel temperature can be made.

preprint2007arXiv

Modulated phases in magnetic models frustrated by long-range interactions

We study an Ising model in one dimension with short range ferromagnetic and long range (power law) antiferromagnetic interactions. We show that the zero temperature phase diagram in a (longitudinal) field H involves a sequence of up and down domains whose size varies continuously with H, between -H_c and H_c which represent the edge of the ferromagnetic up and down phases. The implications of long range interaction in many body systems are discussed.

preprint2005arXiv

On the magnetization of two-dimensional superconductors

We calculate the magnetization of a two-dimensional superconductor in a perpendicular magnetic field near its Kosterlitz-Thouless transition and at lower temperatures. We find that the critical behavior is more complex than assumed in the literature and that, in particular, the critical magnetization is {\it not} field independent as naive scaling predicts. In the low temperature phase we find a substantial fluctuation renormalization of the mean-field result. We compare our analysis with the data on the cuprates.

preprint2003arXiv

Coulomb and Liquid Dimer Models in Three Dimensions

We study classical hard-core dimer models on three-dimensional lattices using analytical approaches and Monte Carlo simulations. On the bipartite cubic lattice, a local gauge field generalization of the height representation used on the square lattice predicts that the dimers are in a critical Coulomb phase with algebraic, dipolar, correlations, in excellent agreement with our large-scale Monte Carlo simulations. The non-bipartite FCC and Fisher lattices lack such a representation, and we find that these models have both confined and exponentially deconfined but no critical phases. We conjecture that extended critical phases are realized only on bipartite lattices, even in higher dimensions.

preprint2000arXiv

Dynamical response of a pinned two-dimensional Wigner crystal

We re-examine a long-standing problem of a finite-frequency conductivity of a weakly pinned two-dimensional classical Wigner crystal. In this system an inhomogeneously broadened absorption line (pinning mode) centered at disorder and magnetic field dependent frequency $ω_p$ is known to appear. We show that the relative linewidth $Δω_p / ω_p$ of the pinning mode is of the order of one in weak magnetic fields, exhibits a power-law decrease in intermediate fields, and eventually saturates at a small value in strong magnetic fields. The linewidth narrowing is due to a peculiar mechanism of mixing between the stiffer longitudinal and the softer transverse components of the collective excitations. The width of the high-field resonance proves to be related to the density of states in the low-frequency tail of the zero-field phonon spectrum. We find a qualitative agreement with recent experiments and point out differences from the previous theoretical work on the subject.

preprint2000arXiv

Superconducting Phase with Fractional Vortices in the Frustrated Kagome Wire Network at f=1/2

In classical XY kagome antiferromagnets, there can be a novel low temperature phase where $ψ^3=e^{i3θ}$ has quasi-long-range order but $ψ$ is disordered, as well as more conventional antiferromagnetic phases where $ψ$ is ordered in various possible patterns ($θ$ is the angle of orientation of the spin). To investigate when these phases exist in a physical system, we study superconducting kagome wire networks in a transverse magnetic field when the magnetic flux through an elementary triangle is a half of a flux quantum. Within Ginzburg-Landau theory, we calculate the helicity moduli of each phase to estimate the Kosterlitz-Thouless (KT) transition temperatures. Then at the KT temperatures, we estimate the barriers to move vortices and effects that lift the large degeneracy in the possible $ψ$ patterns. The effects we have considered are inductive couplings, non-zero wire width, and the order-by-disorder effect due to thermal fluctuations. The first two effects prefer $q=0$ patterns while the last one selects a $\sqrt{3}\times\sqrt{3}$ pattern of supercurrents. Using the parameters of recent experiments, we conclude that at the KT temperature, the non-zero wire width effect dominates, which stabilizes a conventional superconducting phase with a $q=0$ current pattern. However, by adjusting the experimental parameters, for example by bending the wires a little, it appears that the novel $ψ^3$ superconducting phase can instead be stabilized. The barriers to vortex motion are low enough that the system can equilibrate into this phase.

preprint1996arXiv

Nonlocal Conductivity in the Vortex-Liquid Regime of a Two-Dimensional Superconductor

We have simulated the time-dependent Ginzburg-Landau equation with thermal fluctuations, to study the nonlocal dc conductivity of a superconducting film. Having examined points in the phase diagram at a wide range of temperatures and fields below the mean-field upper critical field, we find a portion of the vortex-liquid regime in which the nonlocal ohmic conductivity in real space is negative over a distance several times the spacing between vortices. The effect is suppressed when driven beyond linear response. Earlier work had predicted the existence of such a regime, due to the high viscosity of a strongly-correlated vortex liquid. This behavior is clearly distinguishable from the monotonic spatial fall-off of the conductivity in the higher temperature or field regimes approaching the normal state. The possibilities for experimental study of the nonlocal transport properties are discussed.

preprint1996arXiv

Quantum Griffiths Singularities in the Transverse-Field Ising Spin Glass

We report a Monte Carlo study of the effects of {\it fluctuations} in the bond distribution of Ising spin glasses in a transverse magnetic field, in the {\it paramagnetic phase} in the $T\to 0$ limit. Rare, strong fluctuations give rise to Griffiths singularities, which can dominate the zero-temperature behavior of these quantum systems, as originally demonstrated by McCoy for one-dimensional ($d=1$) systems. Our simulations are done on a square lattice in $d=2$ and a cubic lattice in $d=3$, for a gaussian distribution of nearest neighbor (only) bonds. In $d=2$, where the {\it linear} susceptibility was found to diverge at the critical transverse field strength $Γ_c$ for the order-disorder phase transition at T=0, the average {\it nonlinear} susceptibility $χ_{nl}$ diverges in the paramagnetic phase for $Γ$ well above $Γ_c$, as is also demonstrated in the accompanying paper by Rieger and Young. In $d=3$, the linear susceptibility remains finite at $Γ_c$, and while Griffiths singularity effects are certainly observable in the paramagnetic phase, the nonlinear susceptibility appears to diverge only rather close to $Γ_c$. These results show that Griffiths singularities remain persistent in dimensions above one (where they are known to be strong), though their magnitude decreases monotonically with increasing dimensionality (there being no Griffiths singularities in the limit of infinite dimensionality).

preprint1995arXiv

Topological Defects in the Random-Field XY Model and Randomly Pinned Vortex Lattices

As a simplified model of randomly pinned vortex lattices or charge-density waves, we study the random-field XY model on square ($d=2$) and simple cubic ($d=3$) lattices. We argue, and confirm in simulations, that the spacing between topological defects (vortices) diverges more strongly than the Imry-Ma pinning length as the random field strength, $H$, is reduced. For $d=3$ the data are consistent with a topological phase transition at a nonzero $H_c$ to a vortex-free pinned phase.

preprint1994arXiv

Nonlocal Conductivity in Type-II Superconductors

Multiterminal transport measurements on YBCO crystals in the vortex liquid regime have shown nonlocal conductivity on length scales up to 50 microns. Motivated by these results we explore the wavevector ({\bf k}) dependence of the dc conductivity tensor, $σ_{μν} ({\bf k})$, in the Meissner, vortex lattice, and disordered phases of a type-II superconductor. Our results are based on time-dependent Ginzburg-Landau (TDGL) theory and on phenomenological arguments. We find four qualitatively different types of behavior. First, in the Meissner phase, the conductivity is infinite at $k=0$ and is a continuous function of $k$, monotonically decreasing with increasing $k$. Second, in the vortex lattice phase, in the absence of pinning, the conductivity is finite (due to flux flow) at $k=0$; it is discontinuous there and remains qualitatively like the Meissner phase for $k>0$. Third, in the vortex liquid regime in a magnetic field and at low temperature, the conductivity is finite, smooth and {\it non-monotonic}, first increasing with $k$ at small $k$ and then decreasing at larger $k$. This third behavior is expected to apply at temperatures just above the melting transition of the vortex lattice, where the vortex liquid shows strong short-range order and a large viscosity. Finally, at higher temperatures in the disordered phase, the conductivity is finite, smooth and again monotonically decreasing with $k$. This last, monotonic behavi or applies in zero magnetic field for the entire disordered phase, i.e. at all temperatures above $T_c$, while in

preprint1993arXiv

Statistical Mechanics of Vortices in Type-II Superconductors

Thermal fluctuations and disorder play an essential role in high-T$_c$ superconductors. After reviewing the mean-field phase diagram we describe significant modifications that result when the effects of finite temperature and disorder are incorporated. Thermal fluctuations cause the melting of the Abrikosov flux line lattice into a vortex liquid at high temperatures, and disorder produces novel glass phases at low temperatures. We analyze phase transitions between these new phases and describe transport properties in various regimes.