Interacting many-body quantum systems and their dynamics, while fundamental to modern science and technology, are formidable to simulate and understand. However, by discovering their symmetries, conservation laws, and integrability, one can unravel their intricacies. Here, using up to 124 qubits of a fully programmable quantum computer, we uncover local conservation laws and integrability in one- and two-dimensional periodically-driven spin lattices in a regime previously inaccessible to such detailed analysis. We focus on the paradigmatic example of disorder-induced ergodicity breaking, where we first benchmark the system crossover into a localized regime through anomalies in the one-particle-density-matrix spectrum and other hallmark signatures. We then demonstrate that this regime stems from hidden local integrals of motion by faithfully reconstructing their quantum operators, thus providing a more detailed portrait of the system’s integrable dynamics. Our results demonstrate a versatile strategy for extracting the hidden dynamical structure from noisy experiments on large-scale quantum computers. Local integrals of motion are useful for understanding emergent integrability in many-body localized systems. Here the authors use a large-scale superconducting quantum processor with up to 124 qubits to simulate many-body dynamics in 1D and 2D systems and demonstrate extraction of local integrals of motion.
Interacting many-body quantum systems and their dynamics, while fundamental to modern science and technology, are formidable to simulate and understand. However, by discovering their symmetries, conservation laws, and integrability, one can unravel their intricacies. Here, using up to 124 qubits of a fully programmable quantum computer, we uncover local conservation laws and integrability in one- and two-dimensional periodically-driven spin lattices in a regime previously inaccessible to such detailed analysis. We focus on the paradigmatic example of disorder-induced ergodicity breaking, where we first benchmark the system crossover into a localized regime through anomalies in the one-particle-density-matrix spectrum and other hallmark signatures. We then demonstrate that this regime stems from hidden local integrals of motion by faithfully reconstructing their quantum operators, thus providing a more detailed portrait of the system’s integrable dynamics. Our results demonstrate a versatile strategy for extracting the hidden dynamical structure from noisy experiments on large-scale quantum computers. Local integrals of motion are useful for understanding emergent integrability in many-body localized systems. Here the authors use a large-scale superconducting quantum processor with up to 124 qubits to simulate many-body dynamics in 1D and 2D systems and demonstrate extraction of local integrals of motion.
Interacting many-body quantum systems and their dynamics, while fundamental to modern science and technology, are formidable to simulate and understand. However, by discovering their symmetries, conservation laws, and integrability one can unravel their intricacies. Here, using up to 124 qubits of a fully programmable quantum computer, we uncover local conservation laws and integrability in one- and two-dimensional periodically-driven spin lattices in a regime previously inaccessible to such detailed analysis. We focus on the paradigmatic example of disorder-induced ergodicity breaking, where we first benchmark the system crossover into a localized regime through anomalies in the one-particle-density-matrix spectrum and other hallmark signatures. We then demonstrate that this regime stems from hidden local integrals of motion by faithfully reconstructing their quantum operators, thus providing a detailed portrait of the system's integrable dynamics. Our results demonstrate a versatile strategy for extracting the hidden dynamical structure from noisy experiments on large-scale quantum computers.
Many-body localization provides a generic mechanism of ergodicity breaking in quantum systems. In contrast to conventional ergodic systems, many-body-localized (MBL) systems are characterized by extensively many local integrals of motion (LIOM), which underlie the absence of transport and thermalization in these systems. Here we report a physically motivated construction of local integrals of motion in the MBL phase. We show that any local operator (e.g., a local particle number or a spin-flip operator), evolved with the system's Hamiltonian and averaged over time, becomes a LIOM in the MBL phase. Such operators have a clear physical meaning, describing the response of the MBL system to a local perturbation. In particular, when a local operator represents a density of some globally conserved quantity, the corresponding LIOM describes how this conserved quantity propagates through the MBL phase. Being uniquely defined and experimentally measurable, these LIOMs provide a natural tool for characterizing the properties of the MBL phase, in both experiments and numerical simulations. We demonstrate the latter by numerically constructing an extensive set of LIOMs in the MBL phase of a disordered spin-chain model. We show that the resulting LIOMs are quasilocal and use their decay to extract the localization length and establish the location of the transition between the MBL and ergodic phases.
Understanding the effects of disorder in quantum many-body systems is one of the central challenges of contemporary statistical physics. Its interplay with interactions and quantum fluctuations profoundly modifies the nature of quantum states and their dynamics, giving rise to phenomena such as localization, ergodicity breaking, and the emergence of complex collective behavior. This thesis presents new results on two classes of disordered quantum systems. The first part is devoted to the study of the many-body localization transition, a paradigmatic example of disorder-induced ergodicity breaking. The critical properties of the transition are investigated using a numerical renormalization-group approach, establishing a connection between its scaling behavior and that of the Anderson localization transition on tree-like graphs. The second part focuses on the zero-temperature quantum spin-glass phase of the two-dimensional bond-disordered Heisenberg model. Combining large-scale numerical simulations with a semiclassical expansion, evidence is provided for the existence of the spin-glass phase. The effective model governing the semiclassical fluctuations is further analyzed to characterize the low-energy excitation spectrum and its localization properties. Finally, the weak-disorder regime is investigated and its implications for the ferromagnetic phase are discussed. Overall, the results presented in this thesis contribute to the understanding of the effects of disorder in interacting quantum many-body systems, encompassing both localization phenomena and the collective behavior of disordered quantum magnets.
A periodically driven quantum many-body system is believed to approach a locally infinite-temperature-like state as dictated by the ergodicity hypothesis, when all the couplings/parameters are of the same order. When the drive is stronger than a threshold value, a large class of systems is shown to break this mandate, exhibiting Dynamical Freezing (DF): ergodicity is broken due to the emergence of local operators with permanent memory of their initial values (PRX 11, 021008, 2021). However, an exception to DF has recently been noted in a particular setting (PRX 16, 011041, 2026). Here we identify the condition for ergodicity-breaking in a system under strong drive, and show how to stabilize DF by exact, local engineering. We further show that the Floquet Hamiltonian Heff=(i/T)ln[U(T)] is non-local when ordered by spin-flips, in a non-DF system, while it turns quasi-local under our simple, local DF engineering that keeps the Hamiltonian H(t) still non-integrable at each t. This perplexing non-universality is absent in equilibrium, where non-integrability is believed to imply ergodicity and thermalization universally. Our work paves a way to engineering novel ergodicity-broken interacting Floquet phases of matter in various quantum platforms, and also endows DF-assisted quantum error mitigation and state stabilization techniques (see Nature 653, 1027 (2026)) with unprecedented versatility and flexibility.
Isolated, many-body quantum systems, evolving under their intrinsic dynamics, exhibit a multitude of exotic phenomena and raise foundational questions about statistical mechanics. A flurry of theoretical work has been devoted to understanding how these systems reach thermal equilibrium in the absence of coupling to an external bath and, when thermalization does not occur, investigating the emergent non-equilibrium phases of matter. With the advent of synthetic quantum systems, such as ultra-cold atoms in optical lattices or trapped ions, these questions are no longer academic and can be directly studied in the laboratory. This dissertation explores the non-equilibrium phenomena that stem from the interplay between interactions, disorder, symmetry, topology, and external driving. First, we study how strong disorder, leading to many-body localization, can arrest the heating of a Floquet system and stabilize symmetry-protected topological order that does not have a static analogue. We analyze its dynamical and entanglement properties, highlight its duality to a discrete time crystal, and propose an experimental implementation in a cold-atom setting.Quenched disorder and the many-body localized state are crucial ingredients in protecting macroscopic quantum coherence. We explore the stability of many-body localization in two and higher dimensions and analyze its robustness to rare regions of weak disorder.We then study a second example of non-thermal behavior, namely integrability. We show that a class of random spin models, realizable in systems of atoms coupled to an optical cavity, gives rise to a rich dynamical phase diagram, which includes regions of integrability, classical chaos, and of a novel integrable structure whose conservation laws are reminiscent of the integrals of motion found in a many-body localized phase.The third group of disordered, non-ergodic systems we consider, spin glasses, have fascinating connections to complexity theory and the hardness of constraint satisfaction. We define a statistical ensemble that interpolates between the classical and quantum limits of such a problem and show that there exists a sharp boundary separating satisfiable and unsatisfiable phases.
Many-body quantum systems typically display fast dynamics and ballistic spreading of information. Here we address the open problem of how slow the dynamics can be after a generic breaking of integrability by local interactions. We develop a method based on degenerate perturbation theory that reveals slow dynamical regimes and delocalization processes in general translation invariant models, along with accurate estimates of their delocalization time scales. Our results shed light on the fundamental questions of the robustness of quantum integrable systems and the possibility of many-body localization without disorder. As an example, we construct a large class of one-dimensional lattice models where, despite the absence of asymptotic localization, the transient dynamics is exceptionally slow, i.e., the dynamics is indistinguishable from that of many-body localized systems for the system sizes and time scales accessible in experiments and numerical simulations.