We introduce a method to efficiently study the dynamical properties of many-body localized systems in the regime of strong disorder and weak interactions. Our method reproduces qualitatively and quantitatively the time evolution with a polynomial effort in system size and independent of the desired time scales. We use our method to study quantum information propagation, correlation functions, and temporal fluctuations in one- and two-dimensional many-body localization systems. Moreover, we outline strategies for a further systematic improvement of the accuracy and we point out relations of our method to recent attempts to simulate the time dynamics of quantum many-body systems in classical or artificial neural networks.
We introduce a method to efficiently study the dynamical properties of many-body localized systems in the regime of strong disorder and weak interactions. Our method reproduces qualitatively and quantitatively the time evolution with a polynomial effort in system size and independent of the desired time scales. We use our method to study quantum information propagation, correlation functions, and temporal fluctuations in one- and two-dimensional many-body localization systems. Moreover, we outline strategies for a further systematic improvement of the accuracy and we point out relations of our method to recent attempts to simulate the time dynamics of quantum many-body systems in classical or artificial neural networks.
Recent experiments in quantum simulators have provided evidence for the many-body localized (MBL) phase in one-dimensional (1D) and two-dimensional (2D) bosonic quantum matter. The theoretical study of such a bosonic MBL phase, however, is a daunting task due to the unbounded nature of its Hilbert space. In this work, we introduce a method to compute the long-time real-time evolution of 1D and 2D bosonic systems in an MBL phase at strong disorder and weak interactions. We focus on local dynamical indicators that are able to distinguish an MBL phase from an Anderson localized one. In particular, we consider the temporal fluctuations of local observables, the spatiotemporal behavior of two-time correlators and out-of-time-correlators. We show that these few-body observables can be computed with a computational effort that depends only polynomially on the system size but is independent of the target time, by extending a recently proposed numerical method [G. De Tomasi, F. Pollmann, and M. Heyl, Phys. Rev. B 99, 241114(R) (2019)] to mixed states and bosons. Our method also allows us to surrogate our numerical study with analytical considerations of the time-dependent behavior of the studied quantities.
In this work, we combine quantum renormalization group approaches with deep artificial neural networks for the description of the real-time evolution in strongly disordered quantum matter. We find that this allows us to accurately compute the long-time coherent dynamics of large many-body localized systems in nonperturbative regimes including the effects of many-body resonances. Concretely, we use this approach to describe the spatiotemporal buildup of many-body localized spin-glass order in random Ising chains. We observe a fundamental difference to a noninteracting Anderson insulating Ising chain, where the order only develops over a finite spatial range. We further apply the approach to strongly disordered two-dimensional Ising models, highlighting that our method can be used also for the description of the real-time dynamics of nonergodic quantum matter in a general context.
The efficient numerical simulation of nonequilibrium real-time evolution in isolated quantum matter constitutes a key challenge for current computational methods. This holds in particular in the regime of two spatial dimensions, whose experimental exploration is currently pursued with strong efforts in quantum simulators. In this work we present a versatile and efficient machine learning inspired approach based on a recently introduced artificial neural network encoding of quantum many-body wave functions. We identify and resolve key challenges for the simulation of time evolution, which previously imposed significant limitations on the accurate description of large systems and long-time dynamics. As a concrete example, we study the dynamics of the paradigmatic two-dimensional transverse-field Ising model, as recently also realized experimentally in systems of Rydberg atoms. Calculating the nonequilibrium real-time evolution across a broad range of parameters, we, for instance, observe collapse and revival oscillations of ferromagnetic order and demonstrate that the reached timescales are comparable to or exceed the capabilities of state-of-the-art tensor network methods.
Dynamics in correlated quantum matter is a hard problem, as its exact solution generally involves a computational effort that grows exponentially with the number of constituents. While remarkable progress has been witnessed in recent years for one-dimensional systems, much less has been achieved for interacting quantum models in higher dimensions, since they incorporate an additional layer of complexity. In this work, we employ a variational method that allows for an efficient and controlled computation of the dynamics of quantum many-body systems in one and higher dimensions. The approach presented here introduces a variational class of wave functions based on complex networks of classical spins akin to artificial neural networks, which can be constructed in a controlled fashion. We provide a detailed prescription for such constructions and illustrate their performance by studying quantum quenches in one- and two-dimensional models. In particular, we investigate the nonequilibrium dynamics of a genuinely interacting two-dimensional lattice gauge theory, the quantum link model, for which we have recently shown---employing the technique discussed thoroughly in this paper---that it features disorder-free localization dynamics [P. Karpov et al., Phys. Rev. Lett. 126, 130401 (2021)]. The present work not only supplies a framework to address purely theoretical questions but also could be used to provide a theoretical description of experiments in quantum simulators, which have recently seen an increased effort targeting two-dimensional geometries. Importantly, our method can be applied to any quantum many-body system with a well-defined classical limit.
Many-body localization is a striking mechanism that prevents interacting quantum systems from thermalizing. The absence of thermalization behavior manifests itself, for example, in a remanence of local particle number configurations, a quantity that is robust over a parameter range. Local particle numbers are directly accessible in programmable quantum simulators, in systems of cold atoms, even in two spatial dimensions. Yet, the classical simulation aimed at building trust in quantum simulations is highly challenging. In this work, we present a comprehensive tensor network simulation of a many-body localized systems in two spatial dimensions using a variant of an infinite projected entangled pair states algorithm. The required translational invariance can be restored by implementing the disorder into an auxiliary spin system, providing an exact disorder average under dynamics. We can quantitatively assess signatures of many-body localization for the infinite system: Our methods are powerful enough to provide crude dynamical estimates for the transition between localized and ergodic phases. Interestingly, in this setting of finitely many disorder values, which we also compare with simulations involving noninteracting fermions and for which we discuss the emergent physics, localization emerges in the interacting regime, for which we provide an intuitive argument, while Anderson localization is absent.
Abstract Conventional many-body quantum systems thermalize under their own dynamics, losing information about their initial configurations to the environment. However, it is known that a strong disorder results in many-body localization (MBL). A closed quantum systems with MBL retains local information even in the presence of interactions. Here, we numerically study the propagation and scrambling of quantum information of a closed system in the MBL phase from an information theoretic perspective. By simulating the dynamics and equilibration of the temporal mutual information for long times, we see that it can distinguish between MBL and ergodic phases.