Dynamical chaos in the integrable Toda chain induced by time discretization.

We use the Toda chain model to demonstrate that numerical simulation of integrable Hamiltonian dynamics using time discretization destroys integrability and induces dynamical chaos. Specifically, we integrate this model with various symplectic integrators parametrized by the time step τ and measure the Lyapunov time TΛ (inverse of the largest Lyapunov exponent Λ). A key observation is that TΛ is finite whenever τ is finite but diverges when τ→0. We compare the Toda chain results with the nonintegrable Fermi-Pasta-Ulam-Tsingou chain dynamics. In addition, we observe a breakdown of the simulations at times TB≫TΛ due to certain positions and momenta becoming extremely large ("Not a Number"). This phenomenon originates from the periodic driving introduced by symplectic integrators and we also identify the concrete mechanism of the breakdown in the case of the Toda chain.

Fragile many-body ergodicity from action diffusion.

Weakly nonintegrable many-body systems can restore ergodicity in distinctive ways depending on the range of the interaction network in action space. Action resonances seed chaotic dynamics into the networks. Long-range networks provide well connected resonances with ergodization controlled by the individual resonance chaos time scales. Short-range networks instead yield a dramatic slowing down of ergodization in action space, and lead to rare resonance diffusion. We use Josephson junction chains as a paradigmatic study case. We exploit finite time average distributions to characterize the thermalizing dynamics of actions. We identify an action resonance diffusion regime responsible for the slowing down. We extract the diffusion coefficient of that slow process and measure its dependence on the proximity to the integrable limit. Independent measures of correlation functions confirm our findings. The observed fragile diffusion is relying on weakly chaotic dynamics in spatially isolated action resonances. It can be suppressed, and ergodization delayed, by adding weak action noise, as a proof of concept.

Dynamical glass in weakly nonintegrable Klein-Gordon chains.

Integrable many-body systems are characterized by a complete set of preserved actions. Close to an integrable limit, a nonintegrable perturbation creates a coupling network in action space which can be short or long ranged. We analyze the dynamics of observables which become the conserved actions in the integrable limit. We compute distributions of their finite time averages and obtain the ergodization time scale T_{E} on which these distributions converge to δ distributions. We relate T_{E} to the statistics of fluctuation times of the observables, which acquire fat-tailed distributions with standard deviations σ_{τ}^{+} dominating the means µ_{τ}^{+} and establish that T_{E}∼(σ_{τ}^{+})^{2}/µ_{τ}^{+}. The Lyapunov time T_{Λ} (the inverse of the largest Lyapunov exponent) is then compared to the above time scales. We use a simple Klein-Gordon chain to emulate long- and short-range coupling networks by tuning its energy density. For long-range coupling networks T_{Λ}≈σ_{τ}^{+}, which indicates that the Lyapunov time sets the ergodization time, with chaos quickly diffusing through the coupling network. For short-range coupling networks we observe a dynamical glass, where T_{E} grows dramatically by many orders of magnitude and greatly exceeds the Lyapunov time, which satisfies T_{Λ}≲µ_{τ}^{+}. This effect arises from the formation of highly fragmented inhomogeneous distributions of chaotic groups of actions, separated by growing volumes of nonchaotic regions. These structures persist up to the ergodization time T_{E}.

Dynamical Glass and Ergodization Times in Classical Josephson Junction Chains.

Models of classical Josephson junction chains turn integrable in the limit of large energy densities or small Josephson energies. Close to these limits the Josephson coupling between the superconducting grains induces a short-range nonintegrable network. We compute distributions of finite-time averages of grain charges and extract the ergodization time T_{E} which controls their convergence to ergodic δ distributions. We relate T_{E} to the statistics of fluctuation times of the charges, which are dominated by fat tails. T_{E} is growing anomalously fast upon approaching the integrable limit, as compared to the Lyapunov time T_{Λ}-the inverse of the largest Lyapunov exponent-reaching astonishing ratios T_{E}/T_{Λ}≥10^{8}. The microscopic reason for the observed dynamical glass is rooted in a growing number of grains evolving over long times in a regular almost integrable fashion due to the low probability of resonant interactions with the nearest neighbors. We conjecture that the observed dynamical glass is a generic property of Josephson junction networks irrespective of their space dimensionality.

Weakly Nonergodic Dynamics in the Gross-Pitaevskii Lattice.

The microcanonical Gross-Pitaevskii (also known as the semiclassical Bose-Hubbard) lattice model dynamics is characterized by a pair of energy and norm densities. The grand canonical Gibbs distribution fails to describe a part of the density space, due to the boundedness of its kinetic energy spectrum. We define Poincaré equilibrium manifolds and compute the statistics of microcanonical excursion times off them. The tails of the distribution functions quantify the proximity of the many-body dynamics to a weakly nonergodic phase, which occurs when the average excursion time is infinite. We find that a crossover to weakly nonergodic dynamics takes place inside the non-Gibbs phase, being unnoticed by the largest Lyapunov exponent. In the ergodic part of the non-Gibbs phase, the Gibbs distribution should be replaced by an unknown modified one. We relate our findings to the corresponding integrable limit, close to which the actions are interacting through a short range coupling network.

Flatbands under correlated perturbations.

Flatband networks are characterized by the coexistence of dispersive and flatbands. Flatbands (FBs) are generated by compact localized eigenstates (CLSs) with local network symmetries, based on destructive interference. Correlated disorder and quasiperiodic potentials hybridize CLSs without additional renormalization, yet with surprising consequences: (i) states are expelled from the FB energy E_{FB}, (ii) the localization length of eigenstates vanishes as ξ∼1/ln(E-E_{FB}), (iii) the density of states diverges logarithmically (particle-hole symmetry) and algebraically (no particle-hole symmetry), and (iv) mobility edge curves show algebraic singularities at E_{FB}. Our analytical results are based on perturbative expansions of the CLSs and supported by numerical data in one and two lattice dimensions.