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

Relationship between nonlinearities and thermalization in classical open systems: The role of the interaction range

Roberto Onofrio, Bala Sundaram

Published May 13, 2022• 5 citations
arXivHTMLPDF

Abstract

We discuss results on the dynamics of thermalization for a model with Gaussian interactions between two classical many-body systems trapped in external harmonic potentials. Previous work showed an approximate power-law scaling of the interaction energy with the number of particles, with particular focus on the dependence of the anomalous exponent on the interaction strength. Here we explore the role of the interaction range in determining anomalous exponents, showing that it is a more relevant parameter to differentiate distinct regimes of responses of the system. More specifically, on varying the interaction range from its largest values while keeping the interaction strength constant, we observe a crossover from an integrable system, approximating the Caldeira-Leggett interaction term in the long-range limit, to an intermediate interaction range in which the system manifests anomalous scaling, and finally to a regime of local interactions in which anomalous scaling disappears. A Fourier analysis of the interaction energy shows that nonlinearities give rise to an effective bath with a broad band of frequencies, even when starting with only two distinct trapping frequencies, yielding efficient thermalization in the intermediate regime of interaction range. We provide qualitative arguments, based on an analogous Fourier analysis of the standard map, supporting the view that anomalous scaling and features of the Fourier spectrum may be used as proxies to identify the role of chaotic dynamics. Our work, that encompasses models developed in different contexts and unifies them in a common framework, may be relevant to the general understanding of the role of nonlinearities in a variety of many-body classical systems, ranging from plasmas to trapped atoms and ions.

DOI
10.1103/physreve.105.054122

Related papers

Universal and anomalous behavior in the thermalization of strongly interacting harmonically trapped gas mixtures

Francisco Jauffred, Roberto Onofrio, Bala Sundaram

Abstract We report on the dynamics of thermalization by extending a generalization of the Caldeira–Leggett model, developed in the context of cold atomic gases confined in a harmonic trap, to higher dimensions. Universal characteristics en route to thermalization which appear to be independent of dimensionality are highlighted, which additionally suggest a scaling analogous to turbulent mixing in fluid dynamics. We then focus on features dependent on dimensionality, with particular regard to the role of angular momentum of the two atomic clouds, having in mind the goal of efficient thermalization between the two species. Finally, by considering asymmetry in species numbers, we find that nonlinear interspecies interactions provide a mode locking mechanism between the majority and minority species, relevant to recent experiments involving Fermi–Bose mixtures in the normal phase.

arXiv:1704.06527•Journal of Physics B Atomic Molecular and Optical Physics•May 12, 2017•8 citations

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

Anomalous heat transport in harmonic long-range chains

Francesco Andreucci

In recent years there has been a remarkable rise of interest in the study of heat transport in long-range systems, i.e, systems in which the interaction between constituents scales as a power-law of their distance. The results present in the literature suggest that heat transport in long-range systems can be anomalous, but the understanding of transport properties in these systems remains a difficult and open problem, especially if one wants to obtain analytical results. We thus decided to analyze heat transport properties in linear long-range chains in one dimension, and this work contains the results of this analysis. We begin by considering the simplest long-range system: the fully-connected chain coupled to two external baths. In this case we are able to analytically extract the scaling of the heat-flux with the system's size: this system thus proved a reference result for more complex interactions. We then consider the case of a power-law long-range chain coupled to two baths: in this case we show numerically that the heat flux scales as non-trivial power of system's size and we provide an analytical estimate of the exponent. Finally, we consider a power-law long-range ring with stochastic collisions: in this case we are able to exactly compute the scaling exponent of the thermal conductivity using the Green-Kubo formula, and we find good agreement between our results and numerical simulations.

•Dec 11, 2023•0 citations

Effective microscopic models for sympathetic cooling of atomic gases

Roberto Onofrio, Bala Sundaram

Thermalization of a system in the presence of a heat bath has been the subject of many theoretical investigations especially in the framework of solid-state physics. In this setting, the presence of a large bandwidth for the frequency distribution of the harmonic oscillators schematizing the heat bath is crucial, as emphasized in the Caldeira-Leggett model. By contrast, ultracold gases in atomic traps oscillate at well-defined frequencies and therefore seem to lie outside the Caldeira-Leggett paradigm. We introduce interaction Hamiltonians which allow us to adapt the model to an atomic physics framework. The intrinsic nonlinearity of these models differentiates them from the original Caldeira-Leggett model and calls for a nontrivial stability analysis to determine effective ranges for the model parameters. These models allow for molecular-dynamics simulations of mixtures of ultracold gases, which is of current relevance for optimizing sympathetic cooling in degenerate Bose-Fermi mixtures.

arXiv:1510.02158•Physical Review A•Sep 28, 2015•9 citations

Scaled Caldeira-Leggett Dynamics: Deterministic Trajectories and the Classical Limit

S. V. Mousavi, S. Miret-Artés

The dynamics of an open quantum system in the high-temperature regime of the Caldeira-Leggett model using a Bohmian language is investigated. By expressing the density matrix in polar form, generalized continuity and Hamilton-Jacobi equations are obtained. Thermal effects appear explicitly in the former, while they influence the latter only indirectly through an effective potential. Remarkably, the effective potential does not vanish in the classical limit and cannot, in general, be decomposed into simple additive quantum and thermal contributions. Instead, it retains a nontrivial structure leading to a residual temperature dependence in the resulting dynamics. As a consequence, the resulting temperature-dependent-classical trajectories remain deterministic even in the presence of a thermal environment. This behavior contrasts with the stochastic dynamics observed in the Langevin description and reflects the ensemble-based nature of the present approach. To further explore the quantum-to-classical transition, a scaled version of the Caldeira-Leggett equation is introduced by scaling both the density matrix and Planck constant through a quantumness parameter. This parameter takes the value one in the fully quantum regime and smoothly approaches zero in the classical limit. Applications to Gaussian states demonstrate a gradual suppression of coherence and interference, governed jointly by thermal effects and the transition parameter.

quant-pharXiv:2608.02054•Physica Scripta•Aug 4, 2026•0 citations

Many-body localization in the age of classical computing *

Piotr Sierant, Maciej Lewenstein, Antonello Scardicchio +2 more

Abstract Statistical mechanics provides a framework for describing the physics of large, complex many-body systems using only a few macroscopic parameters to determine the state of the system. For isolated quantum many-body systems, such a description is achieved via the eigenstate thermalization hypothesis (ETH), which links thermalization, ergodicity and quantum chaotic behavior. However, tendency towards thermalization is not observed at finite system sizes and evolution times in a robust many-body localization (MBL) regime found numerically and experimentally in the dynamics of interacting many-body systems at strong disorder. Although the phenomenology of the MBL regime is well-established, the central question remains unanswered: under what conditions does the MBL regime give rise to an MBL phase , in which the thermalization does not occur even in the asymptotic limit of infinite system size and evolution time? This review focuses on recent numerical investigations aiming to clarify the status of the MBL phase, and it establishes the critical open questions about the dynamics of disordered many-body systems. The last decades of research have brought an unprecedented new variety of tools and indicators to study the breakdown of ergodicity, ranging from spectral and wave function measures, matrix elements of observables, through quantities probing unitary quantum dynamics, to transport and quantum information measures. We give a comprehensive overview of these approaches and attempt to provide a unified understanding of their main features. We emphasize general trends towards ergodicity with increasing length and time scales, which exclude naive single-parameter scaling hypothesis, necessitate the use of more refined scaling procedures, and prevent unambiguous extrapolations of numerical results to the asymptotic limit. Providing a concise description of numerical methods for studying ETH and MBL, we explore various approaches to tackle the question of the MBL phase. Persistent finite size drifts towards ergodicity consistently emerge in quantities derived from eigenvalues and eigenvectors of disordered many-body systems. The drifts are related to continuous inching towards ergodicity and non-vanishing transport observed in the dynamics of many-body systems, even at strong disorder. These phenomena impede the understanding of microscopic processes at the ETH-MBL crossover. Nevertheless, the abrupt slowdown of dynamics with increasing disorder strength provides premises suggesting the proximity of the MBL phase. This review concludes that the questions about thermalization and its failure in disordered many-body systems remain a captivating area open for further explorations.

arXiv:2403.07111•Reports on Progress in Physics•Nov 26, 2024•143 citations
We discuss results on the dynamics of thermalization for a model with Gaussian interactions between two classical many-body systems trapped in external harmonic potentials. Previous work showed an approximate power-law scaling of the interaction energy with the number of particles, with particular focus on the dependence of the anomalous exponent on the interaction strength. Here we explore the role of the interaction range in determining anomalous exponents, showing that it is a more relevant parameter to differentiate distinct regimes of responses of the system. More specifically, on varying the interaction range from its largest values while keeping the interaction strength constant, we observe a crossover from an integrable system, approximating the Caldeira-Leggett interaction term in the long-range limit, to an intermediate interaction range in which the system manifests anomalous scaling, and finally to a regime of local interactions in which anomalous scaling disappears. A Fourier analysis of the interaction energy shows that nonlinearities give rise to an effective bath with a broad band of frequencies, even when starting with only two distinct trapping frequencies, yielding efficient thermalization in the intermediate regime of interaction range. We provide qualitative arguments, based on an analogous Fourier analysis of the standard map, supporting the view that anomalous scaling and features of the Fourier spectrum may be used as proxies to identify the role of chaotic dynamics. Our work, that encompasses models developed in different contexts and unifies them in a common framework, may be relevant to the general understanding of the role of nonlinearities in a variety of many-body classical systems, ranging from plasmas to trapped atoms and ions.
Abstract We report on the dynamics of thermalization by extending a generalization of the Caldeira–Leggett model, developed in the context of cold atomic gases confined in a harmonic trap, to higher dimensions. Universal characteristics en route to thermalization which appear to be independent of dimensionality are highlighted, which additionally suggest a scaling analogous to turbulent mixing in fluid dynamics. We then focus on features dependent on dimensionality, with particular regard to the role of angular momentum of the two atomic clouds, having in mind the goal of efficient thermalization between the two species. Finally, by considering asymmetry in species numbers, we find that nonlinear interspecies interactions provide a mode locking mechanism between the majority and minority species, relevant to recent experiments involving Fermi–Bose mixtures in the normal phase.
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.
In recent years there has been a remarkable rise of interest in the study of heat transport in long-range systems, i.e, systems in which the interaction between constituents scales as a power-law of their distance. The results present in the literature suggest that heat transport in long-range systems can be anomalous, but the understanding of transport properties in these systems remains a difficult and open problem, especially if one wants to obtain analytical results. We thus decided to analyze heat transport properties in linear long-range chains in one dimension, and this work contains the results of this analysis. We begin by considering the simplest long-range system: the fully-connected chain coupled to two external baths. In this case we are able to analytically extract the scaling of the heat-flux with the system's size: this system thus proved a reference result for more complex interactions. We then consider the case of a power-law long-range chain coupled to two baths: in this case we show numerically that the heat flux scales as non-trivial power of system's size and we provide an analytical estimate of the exponent. Finally, we consider a power-law long-range ring with stochastic collisions: in this case we are able to exactly compute the scaling exponent of the thermal conductivity using the Green-Kubo formula, and we find good agreement between our results and numerical simulations.
Thermalization of a system in the presence of a heat bath has been the subject of many theoretical investigations especially in the framework of solid-state physics. In this setting, the presence of a large bandwidth for the frequency distribution of the harmonic oscillators schematizing the heat bath is crucial, as emphasized in the Caldeira-Leggett model. By contrast, ultracold gases in atomic traps oscillate at well-defined frequencies and therefore seem to lie outside the Caldeira-Leggett paradigm. We introduce interaction Hamiltonians which allow us to adapt the model to an atomic physics framework. The intrinsic nonlinearity of these models differentiates them from the original Caldeira-Leggett model and calls for a nontrivial stability analysis to determine effective ranges for the model parameters. These models allow for molecular-dynamics simulations of mixtures of ultracold gases, which is of current relevance for optimizing sympathetic cooling in degenerate Bose-Fermi mixtures.
The dynamics of an open quantum system in the high-temperature regime of the Caldeira-Leggett model using a Bohmian language is investigated. By expressing the density matrix in polar form, generalized continuity and Hamilton-Jacobi equations are obtained. Thermal effects appear explicitly in the former, while they influence the latter only indirectly through an effective potential. Remarkably, the effective potential does not vanish in the classical limit and cannot, in general, be decomposed into simple additive quantum and thermal contributions. Instead, it retains a nontrivial structure leading to a residual temperature dependence in the resulting dynamics. As a consequence, the resulting temperature-dependent-classical trajectories remain deterministic even in the presence of a thermal environment. This behavior contrasts with the stochastic dynamics observed in the Langevin description and reflects the ensemble-based nature of the present approach. To further explore the quantum-to-classical transition, a scaled version of the Caldeira-Leggett equation is introduced by scaling both the density matrix and Planck constant through a quantumness parameter. This parameter takes the value one in the fully quantum regime and smoothly approaches zero in the classical limit. Applications to Gaussian states demonstrate a gradual suppression of coherence and interference, governed jointly by thermal effects and the transition parameter.
Abstract Statistical mechanics provides a framework for describing the physics of large, complex many-body systems using only a few macroscopic parameters to determine the state of the system. For isolated quantum many-body systems, such a description is achieved via the eigenstate thermalization hypothesis (ETH), which links thermalization, ergodicity and quantum chaotic behavior. However, tendency towards thermalization is not observed at finite system sizes and evolution times in a robust many-body localization (MBL) regime found numerically and experimentally in the dynamics of interacting many-body systems at strong disorder. Although the phenomenology of the MBL regime is well-established, the central question remains unanswered: under what conditions does the MBL regime give rise to an MBL phase , in which the thermalization does not occur even in the asymptotic limit of infinite system size and evolution time? This review focuses on recent numerical investigations aiming to clarify the status of the MBL phase, and it establishes the critical open questions about the dynamics of disordered many-body systems. The last decades of research have brought an unprecedented new variety of tools and indicators to study the breakdown of ergodicity, ranging from spectral and wave function measures, matrix elements of observables, through quantities probing unitary quantum dynamics, to transport and quantum information measures. We give a comprehensive overview of these approaches and attempt to provide a unified understanding of their main features. We emphasize general trends towards ergodicity with increasing length and time scales, which exclude naive single-parameter scaling hypothesis, necessitate the use of more refined scaling procedures, and prevent unambiguous extrapolations of numerical results to the asymptotic limit. Providing a concise description of numerical methods for studying ETH and MBL, we explore various approaches to tackle the question of the MBL phase. Persistent finite size drifts towards ergodicity consistently emerge in quantities derived from eigenvalues and eigenvectors of disordered many-body systems. The drifts are related to continuous inching towards ergodicity and non-vanishing transport observed in the dynamics of many-body systems, even at strong disorder. These phenomena impede the understanding of microscopic processes at the ETH-MBL crossover. Nevertheless, the abrupt slowdown of dynamics with increasing disorder strength provides premises suggesting the proximity of the MBL phase. This review concludes that the questions about thermalization and its failure in disordered many-body systems remain a captivating area open for further explorations.