Emergent hydrodynamics (EHD) bridges short-time unitarity with late-time thermodynamics, universal transport phenomena characterize the manner and speed of transport and thermalization. Typical non-integrable systems with few conserved local quantities are expected to be diffusive. In contrast, strongly disordered systems which admit phases such as many-body localization, are predicted to inhibit thermalization and thus slow dynamical transport. Disordered systems represent a uniquely poised platform to probe the quantum-to-classical transition and the emergence of irreversible thermodynamics from the underlying unitary structure. Here, we study a strongly disordered nuclear spin ensemble, using local measurements enabled by the disordered-state technique. We observe an apparent phase transition into a sub-diffusive regime, which we model as a random walk on the emergent fractal structure of a percolating network in the dipolar spin ensemble. Our novel theoretical model provides a framework for characterizing the emergence of thermalization in closed quantum systems, even in the presence of strong disorder.
Emergent hydrodynamics (EHD) bridges short-time unitarity with late-time thermodynamics, universal transport phenomena characterize the manner and speed of transport and thermalization. Typical non-integrable systems with few conserved local quantities are expected to be diffusive. In contrast, strongly disordered systems which admit phases such as many-body localization, are predicted to inhibit thermalization and thus slow dynamical transport. Disordered systems represent a uniquely poised platform to probe the quantum-to-classical transition and the emergence of irreversible thermodynamics from the underlying unitary structure. Here, we study a strongly disordered nuclear spin ensemble, using local measurements enabled by the disordered-state technique. We observe an apparent phase transition into a sub-diffusive regime, which we model as a random walk on the emergent fractal structure of a percolating network in the dipolar spin ensemble. Our novel theoretical model provides a framework for characterizing the emergence of thermalization in closed quantum systems, even in the presence of strong disorder.
Understanding how closed quantum systems dynamically approach thermal equilibrium presents a major unresolved problem in statistical physics. Generically, non-integrable quantum systems are expected to thermalize as they comply with the Eigenstate Thermalization Hypothesis. However, in the presence of strong disorder, the dynamics can possibly slow down to a degree that systems fail to thermalize on experimentally accessible timescales, as in spin glasses or many-body localized systems. In general, particularly in long-range interacting quantum systems, the specific nature of the disorder necessary for the emergence of a prethermal, metastable state--distinctly separating the timescales of initial relaxation and subsequent slow thermalization--remains an open question. We study an ensemble of Heisenberg spins with a tunable distribution of random coupling strengths realized by a Rydberg quantum simulator. We observe a drastic change in the late-time magnetization when increasing disorder strength. The data is well described by models based on pairs of strongly interacting spins, which are treated as thermal for weak disorder and isolated for strong disorder. Our results indicate a crossover into a pair-localized prethermal regime in a closed quantum system of thousands of spins in the critical case where the exponent of the power law interaction matches the spatial dimension.
We examine the standard model of many-body localization (MBL), i.e., the disordered chain of interacting spinless fermions, by representing it as the network in the many-body (MB) basis of noninteracting localized Anderson states. By studying eigenstates of the full Hamiltonian, for strong disorders we find that the dynamics is confined up to very long times to disconnected MB clusters in the Fock space. By keeping only resonant contributions and simplifying the quantum problem to rate equations (REs) for MB states, in analogy with percolation problems, the MBL transition is located via the universal cluster distribution and the emergence of the macroscopic cluster. On the ergodic side, our approximate RE approach to the relaxation processes captures well the diffusion transport, as found for the full quantum model. In a broad transient regime, we find an anomalous, i.e., subdiffusivelike, transport, emerging from weak links between MB states.
Elucidating the emergence of irreversible macroscopic laws from reversible quantum many-body dynamics is a question of broad importance across all quantum science. Many-body decoherence plays a key role in this transition, yet connecting microscopic dynamics to emergent macroscopic behavior remains challenging. Here, in a doubly disordered electron-nuclear spin network, we uncover an emergent decoherence law for nuclear polarization, e−Rpte−Rdt, that is robust across broad parameter regimes. We trace its microscopic origins to two interdependent decoherence channels: long-range interactions mediated by the electron network and spin transport within the nuclear network exhibiting anomalous, sub-diffusive dynamics. We demonstrate the capacity to control--and even eliminate--either channel individually through a combination of Floquet engineering and (optical) environment modulation. We find that disorder, typically viewed as detrimental, here proves protective, generating isolated electron-free clusters that localize polarization and prolong coherence lifetimes. These findings establish a microscopic framework for manipulating decoherence pathways and suggests engineered disorder as a new design principle for realizing long-lived quantum memories and sensors.
We introduce and study a toy model for anomalous transport and Griffiths effects in one-dimensional (1D) quantum disordered isolated systems near the many-body localization (MBL) transition. The model is constituted by a collection of 1D tight-binding chains with on-site random energies, locally coupled to a weak perturbation, corresponding to a random matrix which belongs to the Gaussian orthogonal ensemble (GOE) and mimics the effect of thermal inclusions due to delocalizing interactions by providing a local broadening of the Poisson spectrum. While in the absence of such a coupling the model is localized as expected for the one-dimensional Anderson model, increasing the coupling with the GOE perturbation we find a transition to a conducting phase driven by the proliferation of quantum avalanches, which does not fit the standard paradigm of the Anderson transition. In particular an intermediate Griffiths region emerges, where exponentially distributed insulating segments coexist with a few rare resonances. Typical correlations decay exponentially fast, while average correlations decay as stretched exponentials and diverge with the length of the chain, indicating that the conducting inclusions have a fractal structure and that the localization length is broadly distributed at the critical point. This behavior is consistent with a Kosterlitz-Thouless-like criticality of the MBL transition. Transport and relaxation are dominated by rare resonances and rare strong insulating regions, and show anomalous behaviors strikingly similar to those observed in recent simulations and experiments in the bad metal delocalized phase preceding MBL. In particular, we find subdiffusive transport and power-law decay of the return probability at large times, with exponents that gradually change as one moves across the intermediate region. Concomitantly, the ac conductivity vanishes near zero frequency with an anomalous power law.
We study the dynamical melting of hot one-dimensional quantum glasses. We focus on quantum glasses that emerge at strong disorder: specifically, many-body localized systems and their critical variants. The latter category includes disordered Ising, Potts, and anyonic chains whose excited-state properties and dynamics can be accessed via real-space renormalization group methods. As disorder is weakened below a critical value these non-thermal quantum glasses melt via a continuous dynamical phase transition into a classical thermal liquid. By accounting for resonant tunneling of energy, we derive and solve an effective model for such quantum-to-classical transitions and compute their universal critical properties. Notably, the classical thermal liquid exhibits a broad regime of anomalously slow sub-diffusive equilibration dynamics and energy transport, crossing over to ordinary classical diffusion at asymptotically long length scales for critical glasses. The subdiffusive regime is characterized by a continuously evolving dynamical critical exponent that diverges with a universal power at the transition. For non-critical, many-body localized glasses, the crossover scale diverges and subdiffusion persists to the thermodynamic limit. Our approach provides a firm microscopic basis for understanding the generic structure of many-body delocalization transitions in one dimension, and also reveals a general scaling relation among the critical exponents of the transition.
We study the real-time dynamics of local occupation numbers in a one-dimensional model of spinless fermions with a random on-site potential for a certain class of initial states. The latter are thermal (mixed or pure) states of the model in the presence of an additional static force but become nonequilibrium states after a sudden removal of this static force. For this class and high temperatures, we show that the induced dynamics is given by a single correlation function at equilibrium, independent of the initial expectation values being prepared close to equilibrium (by a weak static force) or far away from equilibrium (by a strong static force). Remarkably, this type of universality holds true in both the ergodic phase and the many-body localized regime. Moreover, it does not depend on the specific choice of a unit cell for the local density. We particularly discuss two important consequences. First, the long-time expectation value of the local density is uniquely determined by the fluctuations of its diagonal matrix elements in the energy eigenbasis. Thus, the validity of the eigenstate thermalization hypothesis is not only a sufficient but also a necessary condition for thermalization. Second, the real-time broadening of density profiles is always given by the current autocorrelation function at equilibrium via a generalized Einstein relation. In the context of transport, we discuss the influence of disorder for large particle-particle interactions, where normal diffusion is known to occur in the disorder-free case. Our results suggest that normal diffusion is stable against weak disorder, while they are consistent with anomalous diffusion for stronger disorder below the localization transition. Particularly, for weak disorder, Gaussian density profiles can be observed for single disorder realizations, which we demonstrate for finite lattices up to 31 sites.