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:2109.01361 · Physical Review Letters

Lyapunov Spectrum Scaling for Classical Many-Body Dynamics Close to Integrability

Merab Malishava, Sergej Flach A

Published Mar 31, 2022• 23 citations
arXivHTMLPDF

Abstract

We propose a novel framework to characterize the thermalization of many-body dynamical systems close to integrable limits using the scaling properties of the full Lyapunov spectrum. We use a classical unitary map model to investigate macroscopic weakly nonintegrable dynamics beyond the limits set by the KAM regime. We perform our analysis in two fundamentally distinct long-range and short-range integrable limits which stem from the type of nonintegrable perturbations. Long-range limits result in a single parameter scaling of the Lyapunov spectrum, with the inverse largest Lyapunov exponent being the only diverging control parameter and the rescaled spectrum approaching an analytical function. Short-range limits result in a dramatic slowing down of thermalization which manifests through the rescaled Lyapunov spectrum approaching a non-analytic function. An additional diverging length scale controls the exponential suppression of all Lyapunov exponents relative to the largest one.

DOI
10.1103/physrevlett.128.134102

Related papers

Lyapunov spectrum scaling for many-body dynamics close to integrability

Merab Malishava, Sergej Flach A

We use unitary map evolution to address the thermalization of many-body dynamical systems close to an integrable limit. Integrable limits are classified as either long-range or short-range, depending on the connectivity between actions. The connectivity type is imposed by the nonintegrable perturbation. We compute the Lyapunov spectra and analyze their scaling properties. Long-range limits result in a single parameter scaling of the Lyapunov spectrum, with the largest Lyapunov exponent being the only control parameter. Short-range networks have an additional second scaling parameter which describes the exponential suppression of all relative to the largest Lyapunov exponent.

•Sep 3, 2021•0 citations

Thermalization of Classical Weakly Nonintegrable Many-Body Systems

Merab Malishava

We devote our studies to the subject of weakly nonintegrable dynamics of systems with a macroscopic number of degrees of freedom. Our main points of interest are the relations between the timescales of thermalization and the timescales of chaotization; the choice of appropriate observables and the structure of equations coupling them; identifying the classes of weakly nonintegrable dynamics and developing tools to diagnose the properties specific to such classes. We discuss the traditional in the field methods for thermalization timescale computation and employ them to study the scaling the timescale with the proximity to the integrable limit. We then elaborate on a novel framework based on the full Lyapunov spectra computation for large systems as a powerful tool for the characterization of weak nonintegrability. In particular, the Lyapunov spectrum scaling offers a quantitative description allowing us to infer the structure of the underlying network of observables. Proximity to integrable limit is associated with the rapid growth of thermalization timescales and, thus, potential numerical challenges. We solve these challenges by performing numerical tests using computationally efficient model - unitary maps. The great advantage of unitary maps for numerical applications is time-discrete error-free evolution. We use these advantages to perform large timescale and system size computations in extreme proximity to the integrable limit. To demonstrate the scope of obtained results we report on the application of the developed framework to Hamiltonian systems.

arXiv:2211.01887•Nov 3, 2022•0 citations

Thermalization slowing down in multidimensional Josephson junction networks

Gabriel M. Lando, Sergej Flach A

We characterize thermalization slowing down of Josephson junction networks in one, two, and three spatial dimensions for systems with hundreds of sites by computing their entire Lyapunov spectra. The ratio of Josephson coupling E_{J} to energy density h controls two different universality classes of thermalization slowing down, namely, the weak-coupling regime, E_{J}/h→0, and the strong-coupling regime, E_{J}/h→∞. We analyze the Lyapunov spectrum by measuring the largest Lyapunov exponent and by fitting the rescaled spectrum with a general ansatz. We then extract two scales: the Lyapunov time (inverse of the largest exponent) and the exponent for the decay of the rescaled spectrum. The two universality classes, which exist irrespective of network dimension, are characterized by different ways the extracted scales diverge. The universality class corresponding to the weak-coupling regime allows for the coexistence of chaos with a large number of near-conserved quantities and is shown to be characterized by universal critical exponents, in contrast with the strong-coupling regime. We expect our findings, which we explain using perturbation theory arguments, to be a general feature of diverse Hamiltonian systems.

arXiv:2306.14398•Physical review. E•Dec 28, 2023•15 citations

Observation of Prethermalization in Weakly Nonintegrable Unitary Maps

Xiaodong Zhang, Gabriel M. Lando, Barbara Dietz +1 more

We investigate prethermalization by studying the statistical properties of the time-dependent largest Lyapunov exponent \(\Lambda(t)\) for unitary-circuit maps upon approaching integrability. We follow the evolution of trajectories for different initial conditions and compute the mean \(\mu(t)\) and standard deviation \(\sigma(t)\) of \(\Lambda(t)\). Thermalization implies a temporal decay \(\sigma \sim t^{-1/2}\) at a converged finite value of \(\mu\). We report prethermalization plateaus that persist for long times where both \(\mu\) and \(\sigma\) appear to have converged to finite values, seemingly implying differing saturated Lyapunov exponent values for different trajectories. The lifetime of such plateaus furnishes a novel time scale characterizing the thermalization dynamics of many-body systems close to integrability. We also find that the plateaus converge to their respective thermal values for long enough times.

Qeios•Feb 7, 2025•0 citations

Observation of prethermalization in weakly nonintegrable unitary maps

Xiaodong Zhang, Gabriel M. Lando, Barbara Dietz +1 more

We investigate prethermalization by studying the statistical properties of the time-dependent largest Lyapunov exponent $Λ(t)$ for unitary-circuit maps upon approaching integrability. We follow the evolution of trajectories for different initial conditions and compute the mean $μ(t)$ and standard deviation $σ(t)$ of $Λ(t)$. Thermalization implies a temporal decay $σ\sim t^{-1/2}$ at a converged finite value of $μ$. We report prethermalization plateaus that persist for long times where both $μ$ and $σ$ appear to have converged to finite values, seemingly implying differing saturated Lyapunov exponent values for different trajectories. The lifetime of such plateaus furnishes a novel time scale characterizing the thermalization dynamics of many-body systems close to integrability. We also find that the plateaus converge to their respective thermal values for long enough times.

arXiv:2501.07956•Jan 14, 2025•0 citations

Observation of prethermalization in weakly nonintegrable unitary maps

Xiaodong Zhang, Gabriel M. Lando, Barbara Dietz +1 more

We investigate prethermalization by studying the statistical properties of the time-dependent largest Lyapunov exponent Λ(t) for unitary-circuit maps upon approaching integrability. We follow the evolution of trajectories for different initial conditions and compute the mean μ(t) and standard deviation σ(t) of Λ(t). Thermalization implies a temporal decay σ∼t−1/2 at a converged finite value of μ. We report prethermalization plateaus that persist for long times where both μ and σ appear to have converged to finite values, seemingly implying differing saturated Lyapunov exponent values for different trajectories. The lifetime of such plateaus furnishes a novel time scale characterizing the thermalization dynamics of many-body systems close to integrability. We also find that the plateaus converge to their respective thermal values for long enough times.

Low Temperature Physics•Jun 1, 2025•1 citation
We propose a novel framework to characterize the thermalization of many-body dynamical systems close to integrable limits using the scaling properties of the full Lyapunov spectrum. We use a classical unitary map model to investigate macroscopic weakly nonintegrable dynamics beyond the limits set by the KAM regime. We perform our analysis in two fundamentally distinct long-range and short-range integrable limits which stem from the type of nonintegrable perturbations. Long-range limits result in a single parameter scaling of the Lyapunov spectrum, with the inverse largest Lyapunov exponent being the only diverging control parameter and the rescaled spectrum approaching an analytical function. Short-range limits result in a dramatic slowing down of thermalization which manifests through the rescaled Lyapunov spectrum approaching a non-analytic function. An additional diverging length scale controls the exponential suppression of all Lyapunov exponents relative to the largest one.
We use unitary map evolution to address the thermalization of many-body dynamical systems close to an integrable limit. Integrable limits are classified as either long-range or short-range, depending on the connectivity between actions. The connectivity type is imposed by the nonintegrable perturbation. We compute the Lyapunov spectra and analyze their scaling properties. Long-range limits result in a single parameter scaling of the Lyapunov spectrum, with the largest Lyapunov exponent being the only control parameter. Short-range networks have an additional second scaling parameter which describes the exponential suppression of all relative to the largest Lyapunov exponent.
We devote our studies to the subject of weakly nonintegrable dynamics of systems with a macroscopic number of degrees of freedom. Our main points of interest are the relations between the timescales of thermalization and the timescales of chaotization; the choice of appropriate observables and the structure of equations coupling them; identifying the classes of weakly nonintegrable dynamics and developing tools to diagnose the properties specific to such classes. We discuss the traditional in the field methods for thermalization timescale computation and employ them to study the scaling the timescale with the proximity to the integrable limit. We then elaborate on a novel framework based on the full Lyapunov spectra computation for large systems as a powerful tool for the characterization of weak nonintegrability. In particular, the Lyapunov spectrum scaling offers a quantitative description allowing us to infer the structure of the underlying network of observables. Proximity to integrable limit is associated with the rapid growth of thermalization timescales and, thus, potential numerical challenges. We solve these challenges by performing numerical tests using computationally efficient model - unitary maps. The great advantage of unitary maps for numerical applications is time-discrete error-free evolution. We use these advantages to perform large timescale and system size computations in extreme proximity to the integrable limit. To demonstrate the scope of obtained results we report on the application of the developed framework to Hamiltonian systems.
We characterize thermalization slowing down of Josephson junction networks in one, two, and three spatial dimensions for systems with hundreds of sites by computing their entire Lyapunov spectra. The ratio of Josephson coupling E_{J} to energy density h controls two different universality classes of thermalization slowing down, namely, the weak-coupling regime, E_{J}/h→0, and the strong-coupling regime, E_{J}/h→∞. We analyze the Lyapunov spectrum by measuring the largest Lyapunov exponent and by fitting the rescaled spectrum with a general ansatz. We then extract two scales: the Lyapunov time (inverse of the largest exponent) and the exponent for the decay of the rescaled spectrum. The two universality classes, which exist irrespective of network dimension, are characterized by different ways the extracted scales diverge. The universality class corresponding to the weak-coupling regime allows for the coexistence of chaos with a large number of near-conserved quantities and is shown to be characterized by universal critical exponents, in contrast with the strong-coupling regime. We expect our findings, which we explain using perturbation theory arguments, to be a general feature of diverse Hamiltonian systems.
We investigate prethermalization by studying the statistical properties of the time-dependent largest Lyapunov exponent Λ(t)\Lambda(t)Λ(t) for unitary-circuit maps upon approaching integrability. We follow the evolution of trajectories for different initial conditions and compute the mean μ(t)\mu(t)μ(t) and standard deviation σ(t)\sigma(t)σ(t) of Λ(t)\Lambda(t)Λ(t). Thermalization implies a temporal decay σ∼t−1/2\sigma \sim t^{-1/2}σ∼t−1/2 at a converged finite value of μ\muμ. We report prethermalization plateaus that persist for long times where both μ\muμ and σ\sigmaσ appear to have converged to finite values, seemingly implying differing saturated Lyapunov exponent values for different trajectories. The lifetime of such plateaus furnishes a novel time scale characterizing the thermalization dynamics of many-body systems close to integrability. We also find that the plateaus converge to their respective thermal values for long enough times.
We investigate prethermalization by studying the statistical properties of the time-dependent largest Lyapunov exponent Λ(t)Λ(t)Λ(t) for unitary-circuit maps upon approaching integrability. We follow the evolution of trajectories for different initial conditions and compute the mean μ(t)μ(t)μ(t) and standard deviation σ(t)σ(t)σ(t) of Λ(t)Λ(t)Λ(t). Thermalization implies a temporal decay σ∼t−1/2σ\sim t^{-1/2}σ∼t−1/2 at a converged finite value of μμμ. We report prethermalization plateaus that persist for long times where both μμμ and σσσ appear to have converged to finite values, seemingly implying differing saturated Lyapunov exponent values for different trajectories. The lifetime of such plateaus furnishes a novel time scale characterizing the thermalization dynamics of many-body systems close to integrability. We also find that the plateaus converge to their respective thermal values for long enough times.
We investigate prethermalization by studying the statistical properties of the time-dependent largest Lyapunov exponent Λ(t) for unitary-circuit maps upon approaching integrability. We follow the evolution of trajectories for different initial conditions and compute the mean μ(t) and standard deviation σ(t) of Λ(t). Thermalization implies a temporal decay σ∼t−1/2 at a converged finite value of μ. We report prethermalization plateaus that persist for long times where both μ and σ appear to have converged to finite values, seemingly implying differing saturated Lyapunov exponent values for different trajectories. The lifetime of such plateaus furnishes a novel time scale characterizing the thermalization dynamics of many-body systems close to integrability. We also find that the plateaus converge to their respective thermal values for long enough times.