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:1707.00393

Finite-time Lyapunov exponents in chaotic Bose-Hubbard chains

R. A. Kidd, M. K. Olsen, Joel Frederick Corney

Published Jul 3, 2017• 0 citations
arXivHTMLPDF

Abstract

Many-site Bose-Hubbard lattices display complex semiclassical dynamics, with both chaotic and regular features. We have characterised chaos in the semiclassical dynamics of short Bose-Hubbard chains using both stroboscopic phase space projections and finite-time Lyapunov exponents. We found that chaos was present for intermediate collisional nonlinearity in the open trimer and quatramer systems, with soft chaos and Kolmogoroff-Arnold-Moser islands evident. We have found that the finite-time Lyapunov exponents are consistent with stroboscopic maps for the prediction of chaos in these small systems. This gives us confidence that the finite-time Lyapunov exponents will be a useful tool for the characterisation of chaos in larger systems, where meaningful phase-space projections are not possible and the dimensionality of the problem can make the standard methods intractable.

DOI
10.48550/arxiv.1707.00393

Related papers

Hyperchaos in a Bose-Hubbard Chain with Rydberg-Dressed Interactions

Gary McCormack, Rejish Nath, Weibin Li

We study the chaos and hyperchaos of Rydberg-dressed Bose–Einstein condensates (BECs) in a one-dimensional optical lattice. Due to the long-range, soft-core interaction between the dressed atoms, the dynamics of the BECs are described by the extended Bose-Hubbard model. In the mean-field regime, we analyze the dynamical stability of the BEC by focusing on the ground state and localized state configurations. Lyapunov exponents of the two configurations are calculated by varying the soft-core interaction strength, potential bias, and length of the lattice. Both configurations can have multiple positive Lyapunov exponents, exhibiting hyperchaotic dynamics. We show the dependence of the number of the positive Lyapunov exponents and the largest Lyapunov exponent on the length of the optical lattice. The largest Lyapunov exponent is directly proportional to areas of phase space encompassed by the associated Poincaré sections. We demonstrate that linear and hysteresis quenches of the lattice potential and the dressed interaction lead to distinct dynamics due to the chaos and hyperchaos. Our work is relevant to current research on chaos as well as collective and emergent nonlinear dynamics of BECs with long-range interactions.

arXiv:2108.09683•Photonics•Dec 5, 2021•7 citations

Chaos in the three-site Bose-Hubbard model: Classical versus quantum

Goran Nakerst, Masudul Haque

We consider a quantum many-body system-the Bose-Hubbard system on three sites-which has a classical limit, and which is neither strongly chaotic nor integrable but rather shows a mixture of the two types of behavior. We compare quantum measures of chaos (eigenvalue statistics and eigenvector structure) in the quantum system, with classical measures of chaos (Lyapunov exponents) in the corresponding classical system. As a function of energy and interaction strength, we demonstrate a strong overall correspondence between the two cases. In contrast to both strongly chaotic and integrable systems, the largest Lyapunov exponent is shown to be a multivalued function of energy.

arXiv:2203.09953•Physical review. E•Feb 21, 2023•25 citations

Chaos and anomalous transport in a semiclassical Bose-Hubbard chain

Dragan Marković, Mihailo Čubrović

We study chaotic dynamics and anomalous transport in a Bose-Hubbard chain in the semiclassical regime (the limit when the number of particles goes to infinity). We find that the system has mixed phase space with both regular and chaotic dynamics, even for long chains with up to 100 wells. The consequence of the mixed phase space is strongly anomalous diffusion in the space of occupation numbers, with a discrete set of transport exponents. After very long times the system crosses over to the hydrodynamic regime with normal diffusion. Anomalous transport is quite universal and almost completely independent of the parameters of the model (Coulomb interaction and chemical potential): It is mainly determined by the initial distribution of particles along the chain. We corroborate our findings by analytical arguments: scaling analysis for the anomalous regime and the Langevin equation for the normal diffusion regime.

arXiv:2308.14720•Physical review. E•Mar 28, 2024•2 citations

Metastability, chaos and spectrum tomography for Bose-Hubbard rings and chains

Rajat, Doron Cohen

We analyze the metastability of Bose-Hubbard condensates for finite-size one-dimensional ring lattices and open chains, using a semiclassical tomographic perspective that emphasizes the relation of the many-body spectrum to the underlying classical phase-space structures. In order to address quantum ergodicity in far-from-equilibrium scenarios of experimental interest, we inspect both local aspects (via Bogoliubov analysis) and global aspects (mixed regular-chaotic dynamics). In particular, we highlight the roles of the two parameters that control metastability, clarify the essential differences between low- and high-dimensional chaos, and show how the dynamical instabilities diminish in the limit of the Gross-Pitaevskii equation. It is somewhat frustrating that with more degrees of freedom, the dynamically metastable islands become better distinct from the ergodic sea, while their borders become ill-defined topologically. This stands in opposition to the very structured phase-space of two-degree-of-freedom systems, as reflected in the tomographic quantum spectrum.

quant-pharXiv:2603.23109•Physical Review Research•Mar 25, 2026•0 citations

Classical out-of-time-order correlators for Bose-Hubbard systems and their relation to the finite-time Lyapunov exponent

Frank Großmann, Marcus W. Beims

The out-of-time-order correlator (OTOC) is studied for a bosonic quantum lattice model. We gain its classical analog through the replacement of both commutators appearing in the quantum correlator by a corresponding Poisson bracket. The evaluation of the Poisson bracket is then performed in a complex-valued description of the Hamiltonian dynamics and, for the initial choice of quantum operators to be site-specific annihilators, turns out to be given by the expectation value of the absolute square of a specific element of the complex-valued monodromy matrix. The growth rate of this expectation value is compared to a typical chaos indicator, the mean finite-time Lyapunov exponent (FTLE). In both cases the numerical phase-space average is weighted by a Wigner function corresponding to a multimode coherent state. For a three-well Bose-Hubbard model in the Mott insulator regime, it is found that although, on the level of single trajectories, FTLE and classical OTOC show similar long-time behavior, after averaging, they exhibit a marked difference [which for purely chaotic initial conditions is close to ln(sqrt[2])], rooting in the different order the logarithm and the average are taken. This observation is an example of the relevance of the fluctuations of the FTLE to correctly explain the quantitative difference between OTOC growth rate and FTLE in a prototypical many-body system.

Physical review. E•Sep 26, 2025•0 citations

Signatures of the interplay between chaos and local criticality on the dynamics of scrambling in many-body systems

Felix Meier, Mathias Steinhuber, Juan Diego Urbina +2 more

Fast scrambling, quantified by the exponential initial growth of out-of-time-ordered correlators (OTOCs), is the ability to efficiently spread quantum correlations among the degrees of freedom of interacting systems and constitutes a characteristic signature of local unstable dynamics. As such, it may equally manifest both in systems displaying chaos or in integrable systems around criticality. Here we go beyond these extreme regimes with an exhaustive study of the interplay between local criticality and chaos right at the intricate phase-space region where the integrability-chaos transition first appears. We address systems with a well-defined classical (mean-field) limit, as coupled large spins and Bose-Hubbard chains, thus allowing for semiclassical analysis. Our aim is to investigate the dependence of the exponential growth of the OTOCs, defining the quantum Lyapunov exponent λ_{q} on quantities derived from the classical system with mixed phase space, specifically the local stability exponent of a fixed point λ_{loc} as well as the maximal Lyapunov exponent λ_{L} of the chaotic region around it. By extensive numerical simulations covering a wide range of parameters we give support to a conjectured linear dependence 2λ_{q}=aλ_{L}+bλ_{loc}, providing a simple route to characterize scrambling at the border between chaos and integrability.

arXiv:2211.12147•Physical review. E•May 2, 2023•5 citations
Many-site Bose-Hubbard lattices display complex semiclassical dynamics, with both chaotic and regular features. We have characterised chaos in the semiclassical dynamics of short Bose-Hubbard chains using both stroboscopic phase space projections and finite-time Lyapunov exponents. We found that chaos was present for intermediate collisional nonlinearity in the open trimer and quatramer systems, with soft chaos and Kolmogoroff-Arnold-Moser islands evident. We have found that the finite-time Lyapunov exponents are consistent with stroboscopic maps for the prediction of chaos in these small systems. This gives us confidence that the finite-time Lyapunov exponents will be a useful tool for the characterisation of chaos in larger systems, where meaningful phase-space projections are not possible and the dimensionality of the problem can make the standard methods intractable.
We study the chaos and hyperchaos of Rydberg-dressed Bose–Einstein condensates (BECs) in a one-dimensional optical lattice. Due to the long-range, soft-core interaction between the dressed atoms, the dynamics of the BECs are described by the extended Bose-Hubbard model. In the mean-field regime, we analyze the dynamical stability of the BEC by focusing on the ground state and localized state configurations. Lyapunov exponents of the two configurations are calculated by varying the soft-core interaction strength, potential bias, and length of the lattice. Both configurations can have multiple positive Lyapunov exponents, exhibiting hyperchaotic dynamics. We show the dependence of the number of the positive Lyapunov exponents and the largest Lyapunov exponent on the length of the optical lattice. The largest Lyapunov exponent is directly proportional to areas of phase space encompassed by the associated Poincaré sections. We demonstrate that linear and hysteresis quenches of the lattice potential and the dressed interaction lead to distinct dynamics due to the chaos and hyperchaos. Our work is relevant to current research on chaos as well as collective and emergent nonlinear dynamics of BECs with long-range interactions.
We consider a quantum many-body system-the Bose-Hubbard system on three sites-which has a classical limit, and which is neither strongly chaotic nor integrable but rather shows a mixture of the two types of behavior. We compare quantum measures of chaos (eigenvalue statistics and eigenvector structure) in the quantum system, with classical measures of chaos (Lyapunov exponents) in the corresponding classical system. As a function of energy and interaction strength, we demonstrate a strong overall correspondence between the two cases. In contrast to both strongly chaotic and integrable systems, the largest Lyapunov exponent is shown to be a multivalued function of energy.
We study chaotic dynamics and anomalous transport in a Bose-Hubbard chain in the semiclassical regime (the limit when the number of particles goes to infinity). We find that the system has mixed phase space with both regular and chaotic dynamics, even for long chains with up to 100 wells. The consequence of the mixed phase space is strongly anomalous diffusion in the space of occupation numbers, with a discrete set of transport exponents. After very long times the system crosses over to the hydrodynamic regime with normal diffusion. Anomalous transport is quite universal and almost completely independent of the parameters of the model (Coulomb interaction and chemical potential): It is mainly determined by the initial distribution of particles along the chain. We corroborate our findings by analytical arguments: scaling analysis for the anomalous regime and the Langevin equation for the normal diffusion regime.
We analyze the metastability of Bose-Hubbard condensates for finite-size one-dimensional ring lattices and open chains, using a semiclassical tomographic perspective that emphasizes the relation of the many-body spectrum to the underlying classical phase-space structures. In order to address quantum ergodicity in far-from-equilibrium scenarios of experimental interest, we inspect both local aspects (via Bogoliubov analysis) and global aspects (mixed regular-chaotic dynamics). In particular, we highlight the roles of the two parameters that control metastability, clarify the essential differences between low- and high-dimensional chaos, and show how the dynamical instabilities diminish in the limit of the Gross-Pitaevskii equation. It is somewhat frustrating that with more degrees of freedom, the dynamically metastable islands become better distinct from the ergodic sea, while their borders become ill-defined topologically. This stands in opposition to the very structured phase-space of two-degree-of-freedom systems, as reflected in the tomographic quantum spectrum.
The out-of-time-order correlator (OTOC) is studied for a bosonic quantum lattice model. We gain its classical analog through the replacement of both commutators appearing in the quantum correlator by a corresponding Poisson bracket. The evaluation of the Poisson bracket is then performed in a complex-valued description of the Hamiltonian dynamics and, for the initial choice of quantum operators to be site-specific annihilators, turns out to be given by the expectation value of the absolute square of a specific element of the complex-valued monodromy matrix. The growth rate of this expectation value is compared to a typical chaos indicator, the mean finite-time Lyapunov exponent (FTLE). In both cases the numerical phase-space average is weighted by a Wigner function corresponding to a multimode coherent state. For a three-well Bose-Hubbard model in the Mott insulator regime, it is found that although, on the level of single trajectories, FTLE and classical OTOC show similar long-time behavior, after averaging, they exhibit a marked difference [which for purely chaotic initial conditions is close to ln(sqrt[2])], rooting in the different order the logarithm and the average are taken. This observation is an example of the relevance of the fluctuations of the FTLE to correctly explain the quantitative difference between OTOC growth rate and FTLE in a prototypical many-body system.
Fast scrambling, quantified by the exponential initial growth of out-of-time-ordered correlators (OTOCs), is the ability to efficiently spread quantum correlations among the degrees of freedom of interacting systems and constitutes a characteristic signature of local unstable dynamics. As such, it may equally manifest both in systems displaying chaos or in integrable systems around criticality. Here we go beyond these extreme regimes with an exhaustive study of the interplay between local criticality and chaos right at the intricate phase-space region where the integrability-chaos transition first appears. We address systems with a well-defined classical (mean-field) limit, as coupled large spins and Bose-Hubbard chains, thus allowing for semiclassical analysis. Our aim is to investigate the dependence of the exponential growth of the OTOCs, defining the quantum Lyapunov exponent λ_{q} on quantities derived from the classical system with mixed phase space, specifically the local stability exponent of a fixed point λ_{loc} as well as the maximal Lyapunov exponent λ_{L} of the chaotic region around it. By extensive numerical simulations covering a wide range of parameters we give support to a conjectured linear dependence 2λ_{q}=aλ_{L}+bλ_{loc}, providing a simple route to characterize scrambling at the border between chaos and integrability.