Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2505.05882 · Physical review. B./Physical review. B

Finding local integrals of motion in quantum lattice models in the thermodynamic limit

Jakub Pawłowski, J. Herbrych, Marcin Mierzejewski

Published Sep 24, 2025• 4 citations
arXivHTMLPDF

Abstract

Local integrals of motion (LIOMs) play a key role in understanding the long-time properties of closed macroscopic systems. They were found for selected integrable systems via complex analytical calculations. The existence of LIOMs and their structure can also be studied via numerical methods, which, however, involve exact diagonalization of Hamiltonians, posing a bottleneck for such studies. We show that finding LIOMs in translationally invariant lattice models or unitary quantum circuits can be reduced to a problem for which one may numerically find an exact solution in the thermodynamic limit. We develop a simple algorithm and demonstrate its efficiency by calculating LIOMs and bounds on correlations (the Mazur bounds) for infinite integrable spin chains and unitary circuits. Finally, we demonstrate that this approach identifies slow modes in nearly integrable spin models and estimates their relaxation times.

DOI
10.1103/brst-wp36

Related papers

Local integrals of motion in many‐body localized systems

John Z. Imbrie, Valentina Ros, Antonello Scardicchio

We review the current (as of Fall 2016) status of the studies on the emergent integrability in many‐body localized models. We start by explaining how the phenomenology of fully many‐body localized systems can be recovered if one assumes the existence of a complete set of (quasi)local operators which commute with the Hamiltonian (local integrals of motions, or LIOMs). We describe the evolution of this idea from the initial conjecture, to the perturbative constructions, to the mathematical proof given for a disordered spin chain. We discuss the proposed numerical algorithms for the construction of LIOMs and the status of the debate on the existence and nature of such operators in systems with a many‐body mobility edge, and in dimensions larger than one. image

arXiv:1609.08076•Annalen der Physik•May 11, 2017•277 citations

Unified theory of local integrals of motion

Ben Craps, Oleg Evnin, Dmitry Kovrizhin +1 more

Conservation laws are of paramount importance in our understanding of classical and quantum dynamics. Here, we present a general framework for constructing exact quantum integrals of motion with the desired locality and quantum numbers, which will be illustrated for the case of many-body-localization (MBL). The latter has been understood theoretically in terms of the existence of an extensive number of local integrals of motion (LIOMs). Using our approach, we show that one can express the task of finding LIOMs as an optimization problem. For some specifications, this problem surprisingly connects to the question of finding classical ground states of spin-glass models. Our work unifies previous results obtained in the MBL context and reveals intriguing connections between many-body localization, spin-glass physics, and constrained optimization problems.

cond-mat.dis-nnCross-listedarXiv:2512.09595•Physical review. B./Physical review. B•Dec 11, 2025•0 citations

Constructing local integrals of motion in the many-body localized phase

Anushya Chandran, Isaac H. Kim, Guifré Vidal +1 more

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.

arXiv:1407.8480•Physical Review B•Feb 24, 2015•301 citations

Quantitative Impact of Integrals of Motion on the Eigenstate Thermalization Hypothesis

Marcin Mierzejewski, Lev Vidmar

Even though the eigenstate thermalization hypothesis (ETH) may be introduced as an extension of the random matrix theory, physical Hamiltonians and observables differ from random operators. One of the challenges is to embed local integrals of motion (LIOMs) within the ETH. Here we make steps towards a unified treatment of the ETH in integrable and nonintegrable models with translational invariance. Specifically, we focus on the impact of LIOMs on the fluctuations and structure of the diagonal matrix elements of local observables. We first show that nonvanishing fluctuations entail the presence of LIOMs. Then we introduce a generic protocol to construct observables, subtracted by their projections on LIOMs as well as products of LIOMs. The protocol systematically reduces fluctuations and/or the structure of the diagonal matrix elements. We verify our arguments by numerical results for integrable and nonintegrable models.

arXiv:1908.08569•Physical Review Letters•Jan 29, 2020•61 citations

Localization and slow-thermalization in a cluster spin model

Yoshihito Kuno, Takahiro Orito, Ikuo Ichinose

Abstract Novel cluster spin model with interactions and disorder is introduced and studied. In specific type of interactions, we find an extensive number of local integrals of motions (LIOMs), which are a modified version of the stabilizers in quantum information, i.e., mutually commuting operators specifying all quantum states in the system. These LIOMs can be defined for any strength of the interactions and disorder, and are of compact-support instead of exponentially-decaying tail. Hence, even under the presence of interactions, integrability is held, and all energy eigenstates are labeled by these LIOMs and can be explicitly obtained. Integrable dynamics is, then, expected to occur. The compact-support nature of the LIOMs crucially prevents the thermalization and entanglement spreading. We numerically investigate dynamics of the system governed by the existence of the compact-support LIOMs, and clarify the effects of additional interactions, which break the compact-support nature of the LIOMs. There, we find that the ordinary many-body localization behaviors emerge, such as the logarithmic growth of the entanglement entropy in the time evolution. Besides the ergodicity breaking dynamic, we find that symmetry protected topological order preserves for specific states even in the presence of the interactions.

arXiv:2110.05970•New Journal of Physics•Jun 29, 2022•2 citations

Local vs. translationally-invariant slowest operators in quantum Ising spin chains

Ekaterina Izotova

In this paper we study one-dimensional quantum Ising spin chains in external magnetic field close to an integrable point. We concentrate on the dynamics of the slowest operator, that plays a key role at the final period of thermalization. We introduce two independent definitions of the slowest operator: local and translationally-invariant ones. We construct both operators numerically using tensor networks and extensively compare their physical properties. We find that the local operator has a significant overlap with energy flux, it does not correspond to an integral of motion, and, as one goes away from the integrable point, its revivals get suppressed and the rate of delocalization changes from extremely slow to slower than diffusion. The translationally-invariant operator corresponds to an integral of motion; as the system becomes less integrable, at some point this operator changes its nature: from no overlap with any magnetization and fast rate of delocalization, to non-zero overlap with magnetizations $σ_{x}$ and $σ_{z}$ and slow rate of delocalization.

arXiv:2302.07995•Feb 16, 2023•0 citations
Local integrals of motion (LIOMs) play a key role in understanding the long-time properties of closed macroscopic systems. They were found for selected integrable systems via complex analytical calculations. The existence of LIOMs and their structure can also be studied via numerical methods, which, however, involve exact diagonalization of Hamiltonians, posing a bottleneck for such studies. We show that finding LIOMs in translationally invariant lattice models or unitary quantum circuits can be reduced to a problem for which one may numerically find an exact solution in the thermodynamic limit. We develop a simple algorithm and demonstrate its efficiency by calculating LIOMs and bounds on correlations (the Mazur bounds) for infinite integrable spin chains and unitary circuits. Finally, we demonstrate that this approach identifies slow modes in nearly integrable spin models and estimates their relaxation times.
We review the current (as of Fall 2016) status of the studies on the emergent integrability in many‐body localized models. We start by explaining how the phenomenology of fully many‐body localized systems can be recovered if one assumes the existence of a complete set of (quasi)local operators which commute with the Hamiltonian (local integrals of motions, or LIOMs). We describe the evolution of this idea from the initial conjecture, to the perturbative constructions, to the mathematical proof given for a disordered spin chain. We discuss the proposed numerical algorithms for the construction of LIOMs and the status of the debate on the existence and nature of such operators in systems with a many‐body mobility edge, and in dimensions larger than one. image
Conservation laws are of paramount importance in our understanding of classical and quantum dynamics. Here, we present a general framework for constructing exact quantum integrals of motion with the desired locality and quantum numbers, which will be illustrated for the case of many-body-localization (MBL). The latter has been understood theoretically in terms of the existence of an extensive number of local integrals of motion (LIOMs). Using our approach, we show that one can express the task of finding LIOMs as an optimization problem. For some specifications, this problem surprisingly connects to the question of finding classical ground states of spin-glass models. Our work unifies previous results obtained in the MBL context and reveals intriguing connections between many-body localization, spin-glass physics, and constrained optimization problems.
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.
Even though the eigenstate thermalization hypothesis (ETH) may be introduced as an extension of the random matrix theory, physical Hamiltonians and observables differ from random operators. One of the challenges is to embed local integrals of motion (LIOMs) within the ETH. Here we make steps towards a unified treatment of the ETH in integrable and nonintegrable models with translational invariance. Specifically, we focus on the impact of LIOMs on the fluctuations and structure of the diagonal matrix elements of local observables. We first show that nonvanishing fluctuations entail the presence of LIOMs. Then we introduce a generic protocol to construct observables, subtracted by their projections on LIOMs as well as products of LIOMs. The protocol systematically reduces fluctuations and/or the structure of the diagonal matrix elements. We verify our arguments by numerical results for integrable and nonintegrable models.
Abstract Novel cluster spin model with interactions and disorder is introduced and studied. In specific type of interactions, we find an extensive number of local integrals of motions (LIOMs), which are a modified version of the stabilizers in quantum information, i.e., mutually commuting operators specifying all quantum states in the system. These LIOMs can be defined for any strength of the interactions and disorder, and are of compact-support instead of exponentially-decaying tail. Hence, even under the presence of interactions, integrability is held, and all energy eigenstates are labeled by these LIOMs and can be explicitly obtained. Integrable dynamics is, then, expected to occur. The compact-support nature of the LIOMs crucially prevents the thermalization and entanglement spreading. We numerically investigate dynamics of the system governed by the existence of the compact-support LIOMs, and clarify the effects of additional interactions, which break the compact-support nature of the LIOMs. There, we find that the ordinary many-body localization behaviors emerge, such as the logarithmic growth of the entanglement entropy in the time evolution. Besides the ergodicity breaking dynamic, we find that symmetry protected topological order preserves for specific states even in the presence of the interactions.
In this paper we study one-dimensional quantum Ising spin chains in external magnetic field close to an integrable point. We concentrate on the dynamics of the slowest operator, that plays a key role at the final period of thermalization. We introduce two independent definitions of the slowest operator: local and translationally-invariant ones. We construct both operators numerically using tensor networks and extensively compare their physical properties. We find that the local operator has a significant overlap with energy flux, it does not correspond to an integral of motion, and, as one goes away from the integrable point, its revivals get suppressed and the rate of delocalization changes from extremely slow to slower than diffusion. The translationally-invariant operator corresponds to an integral of motion; as the system becomes less integrable, at some point this operator changes its nature: from no overlap with any magnetization and fast rate of delocalization, to non-zero overlap with magnetizations σxσ_{x}σx​ and σzσ_{z}σz​ and slow rate of delocalization.