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:2306.14398 · Physical review. E

Thermalization slowing down in multidimensional Josephson junction networks

Gabriel M. Lando, Sergej Flach A

Published Dec 28, 2023• 15 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreve.108.l062301

Related papers

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

Merab Malishava, Sergej Flach A

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.

arXiv:2109.01361•Physical Review Letters•Mar 31, 2022•23 citations

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

Understanding thermalization of the magnetic gluons in a non-Abelian gauge theory

Sayak Guin, Harshit Pandey, Sayantan Sharma

We explore chaotic dynamics in SU(2) gauge theory by extracting the Lyapunov exponent under both in and out-of-equilibrium conditions to understand how chaos develops and influences the approach to thermalization. In thermal equilibrium at very high temperatures, the hard, electric, and magnetic scales are well separated, and one can write an effective theory for the soft gauge fields (magnetic modes) which are classical due to their large occupation number. In the vicinity of the critical temperature $T_c$, $SU(2)$ gauge theory falls within the same universality class as $Z_2$ scalar field theory. Exploiting this correspondence, we determine the Lyapunov exponent in the scalar theory through out-of-time-ordered correlators (OTOC). To probe non-equilibrium dynamics, we prepare the gauge fields in a highly occupied, non-thermal initial state. Under classical Yang–Mills evolution, the system evolves toward a self-similar scaling regime characterized by a scale separation similar to that of the thermal case. Within this regime, we again compute the Lyapunov exponent and compare it to its equilibrium counterpart at similar energy densities. This comparison allows us to quantify the timescale over which the system approaches thermal equilibrium.

International Symposium on Lattice Field Theory (LATTICE)•Apr 14, 2026•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 universality-class transition induced by Anderson localization

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

We study the disorder-induced crossover between the two recently discovered thermalization slowing-down universality classes, characterized by long- and short-range coupling, in classical unitary-circuit maps close to integrability. We compute Lyapunov spectra, which display qualitatively distinct features depending on whether the proximity to the integrable limit is short or long range. For sufficiently small nonlinearity, translationally invariant systems fall into the long-range class. Adding disorder to such a system triggers a transition to the short-range class, implying a breaking of this invariance, and in the very limit of vanishing nonlinearity Anderson localization emerges. The crossover from the long- to the short-range class is attained by tuning the localization length ξ from ξ≈N to ξ≪N , where N is the system size. As a consequence, the Lyapunov spectrum becomes exponentially suppressed, depending on the extent to which its translational invariance is destroyed. We expect that this disorder-induced crossover will lead to prethermalized phases and, following quantization, to many-body localization. Published by the American Physical Society 2024

arXiv:2308.08921•Physical Review Research•Mar 25, 2024•9 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
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 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 explore chaotic dynamics in SU(2) gauge theory by extracting the Lyapunov exponent under both in and out-of-equilibrium conditions to understand how chaos develops and influences the approach to thermalization. In thermal equilibrium at very high temperatures, the hard, electric, and magnetic scales are well separated, and one can write an effective theory for the soft gauge fields (magnetic modes) which are classical due to their large occupation number. In the vicinity of the critical temperature TcT_cTc​, SU(2)SU(2)SU(2) gauge theory falls within the same universality class as Z2Z_2Z2​ scalar field theory. Exploiting this correspondence, we determine the Lyapunov exponent in the scalar theory through out-of-time-ordered correlators (OTOC). To probe non-equilibrium dynamics, we prepare the gauge fields in a highly occupied, non-thermal initial state. Under classical Yang–Mills evolution, the system evolves toward a self-similar scaling regime characterized by a scale separation similar to that of the thermal case. Within this regime, we again compute the Lyapunov exponent and compare it to its equilibrium counterpart at similar energy densities. This comparison allows us to quantify the timescale over which the system approaches thermal equilibrium.
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 study the disorder-induced crossover between the two recently discovered thermalization slowing-down universality classes, characterized by long- and short-range coupling, in classical unitary-circuit maps close to integrability. We compute Lyapunov spectra, which display qualitatively distinct features depending on whether the proximity to the integrable limit is short or long range. For sufficiently small nonlinearity, translationally invariant systems fall into the long-range class. Adding disorder to such a system triggers a transition to the short-range class, implying a breaking of this invariance, and in the very limit of vanishing nonlinearity Anderson localization emerges. The crossover from the long- to the short-range class is attained by tuning the localization length ξ from ξ≈N to ξ≪N , where N is the system size. As a consequence, the Lyapunov spectrum becomes exponentially suppressed, depending on the extent to which its translational invariance is destroyed. We expect that this disorder-induced crossover will lead to prethermalized phases and, following quantization, to many-body localization. Published by the American Physical Society 2024
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.