Weakly Nonergodic Dynamics in the Gross-Pitaevskii Lattice

How close are integrable and nonintegrable models: A parametric case study based on the Salerno model

In the present work we revisit the Salerno model as a prototypical system that interpolates between a well-known integrable system (the Ablowitz-Ladik lattice) and an experimentally tractable, nonintegrable one (the discrete nonlinear Schrödinger model). The question we ask is, for "generic" initial data, how close are the integrable to the nonintegrable models? Our more precise formulation of this question is, How well is the constancy of formerly conserved quantities preserved in the nonintegrable case? Upon examining this, we find that even slight deviations from integrability can be sensitively felt by measuring these formerly conserved quantities in the case of the Salerno model. However, given that the knowledge of these quantities requires a deep physical and mathematical analysis of the system, we seek a more "generic" diagnostic towards a manifestation of integrability breaking. We argue, based on our Salerno model computations, that the full spectrum of Lyapunov exponents could be a sensitive diagnostic to that effect.

Condensation induced by coupled transport processes

Several lattice models display a condensation transition in real space when the density of a suitable order parameter exceeds a critical value. We consider one of such models with two conservation laws, in a onedimensional open setup where the system is attached to two external reservoirs. Both reservoirs impose subcritical boundary conditions at the chain ends. When such boundary conditions are equal, the system is in equilibrium below the condensation threshold and no condensate can appear. Instead, when the system is kept out of equilibrium, localization may arise in an internal portion of the lattice. We discuss the origin of this phenomenon, the relevance of the number of conservation laws, and the effect of the pinning of the condensate on the dynamics of the out-of-equilibrium state.

Chaos and thermalization in a classical chain of dipoles

We explore the connection between chaos, thermalization, and ergodicity in a linear chain of N interacting dipoles. Starting from the ground state, and considering chains of different numbers of dipoles, we introduce single site excitations with excess energy ΔK. The time evolution of the chaoticity of the system and the energy localization along the chain is analyzed by computing, up to a very long time, the statistical average of the finite-time Lyapunov exponent λ(t) and the participation ratio Π(t). For small ΔK, the evolution of λ(t) and Π(t) indicates that the system becomes chaotic at approximately the same time as Π(t) reaches a steady state. For the largest considered values of ΔK the system becomes chaotic at an extremely early stage in comparison with the energy relaxation times. We find that this fact is due to the presence of chaotic breathers that keep the system far from equipartition and ergodicity. Finally, we show numerically and analytically that the asymptotic values attained by the participation ratio Π(t) fairly correspond to thermal equilibrium.

Heat transport in long-ranged Fermi-Pasta-Ulam-Tsingou-type models

We perform a detailed study of heat transport in one-dimensional long-ranged anharmonic oscillator systems, such as the long-ranged Fermi-Pasta-Ulam-Tsingou model. For these systems, the long-ranged anharmonic potential decays with distance as a power law, controlled by an exponent δ≥0. For such a nonintegrable model, one of the recent results that has captured quite some attention is the puzzling ballisticlike transport observed for δ=2, reminiscent of integrable systems. Here, we first employ the reverse nonequilibrium molecular dynamics simulations to look closely at the δ=2 transport in three long-ranged models and point out a few problematic issues with this simulation method. Next, we examine the process of energy relaxation, and find that relaxation can be appreciably slow for δ=2 in some situations. We invoke the concept of nonlinear localized modes of excitation, also known as discrete breathers, and demonstrate that the slow relaxation and the ballisticlike transport properties can be consistently explained in terms of a novel depinning of the discrete breathers that makes them highly mobile at δ=2. Finally, in the presence of quartic pinning potentials we find that the long-ranged model exhibits Fourier (diffusive) transport at δ=2, as one would expect from short-ranged interacting systems with broken momentum conservation. Such a diffusive regime is not observed for harmonic pinning.

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.

Four-wave mixing and coherently coupled Schrödinger equations: Cascading processes and Fermi-Pasta-Ulam-Tsingou recurrence

Modulation instability, breather formation, and the Fermi-Pasta-Ulam-Tsingou recurrence (FPUT) phenomena are studied in this article. Physically, such nonlinear systems arise when the medium is slightly anisotropic, e.g., optical fibers with weak birefringence where the slowly varying pulse envelopes are governed by these coherently coupled Schrödinger equations. The Darboux transformation is used to calculate a class of breathers where the carrier envelope depends on the transverse coordinate of the Schrödinger equations. A "cascading mechanism" is utilized to elucidate the initial stages of FPUT. More precisely, higher order nonlinear terms that are exponentially small initially can grow rapidly. A breather is formed when the linear mode and higher order ones attain roughly the same magnitude. The conditions for generating various breathers and connections with modulation instability are elucidated. The growth phase then subsides and the cycle is repeated, leading to FPUT. Unequal initial conditions for the two waveguides produce symmetry breaking, with "eye-shaped" breathers in one waveguide and "four-petal" modes in the other. An analytical formula for the time or distance of breather formation for a two-waveguide system is proposed, based on the disturbance amplitude and instability growth rate. Excellent agreement with numerical simulations is achieved. Furthermore, the roles of modulation instability for FPUT are elucidated with illustrative case studies. In particular, depending on whether the second harmonic falls within the unstable band, FPUT patterns with one single or two distinct wavelength(s) are observed. For applications to temporal optical waveguides, the present formulation can predict the distance along a weakly birefringent fiber needed to observe FPUT.

Thermalization in the one-dimensional Salerno model lattice

The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (nonintegrable) discrete nonlinear Schrödinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS toward AL) is varied, the region in the space of initial energy and norm densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes.

Condensation transition and ensemble inequivalence in the discrete nonlinear Schrödinger equation

The thermodynamics of the discrete nonlinear Schrödinger equation in the vicinity of infinite temperature is explicitly solved in the microcanonical ensemble by means of large-deviation techniques. A first-order phase transition between a thermalized phase and a condensed (localized) one occurs at the infinite-temperature line. Inequivalence between statistical ensembles characterizes the condensed phase, where the grand-canonical representation does not apply. The control over finite-size corrections of the microcanonical partition function allows to design an experimental test of delocalized negative-temperature states in lattices of cold atoms.

Probing Band Topology Using Modulational Instability

We analyze the modulational instability of nonlinear Bloch waves in topological photonic lattices. In the initial phase of the instability development captured by the linear stability analysis, long wavelength instabilities and bifurcations of the nonlinear Bloch waves are sensitive to topological band inversions. At longer timescales, nonlinear wave mixing induces spreading of energy through the entire band and spontaneous creation of wave polarization singularities determined by the band Chern number. Our analytical and numerical results establish modulational instability as a tool to probe bulk topological invariants and create topologically nontrivial wave fields.

Chaotic wave-packet spreading in two-dimensional disordered nonlinear lattices

We reveal the generic characteristics of wave-packet delocalization in two-dimensional nonlinear disordered lattices by performing extensive numerical simulations in two basic disordered models: the Klein-Gordon system and the discrete nonlinear Schrödinger equation. We find that in both models (a) the wave packet's second moment asymptotically evolves as t^{a_{m}} with a_{m}≈1/5 (1/3) for the weak (strong) chaos dynamical regime, in agreement with previous theoretical predictions [S. Flach, Chem. Phys. 375, 548 (2010)CMPHC20301-010410.1016/j.chemphys.2010.02.022]; (b) chaos persists, but its strength decreases in time t since the finite-time maximum Lyapunov exponent Λ decays as Λ∝t^{α_{Λ}}, with α_{Λ}≈-0.37 (-0.46) for the weak (strong) chaos case; and (c) the deviation vector distributions show the wandering of localized chaotic seeds in the lattice's excited part, which induces the wave packet's thermalization. We also propose a dimension-independent scaling between the wave packet's spreading and chaoticity, which allows the prediction of the obtained α_{Λ} values.

Transient diffusion and two-regime localization of discrete breatherlike excitations in nonlinear Schrödinger lattice with disorder

We systematically simulate and analyze the motion of discrete breatherlike excitations (DBEs) in the nonlinear Schrödinger lattice with random potentials. A universal transient diffusion of the DBEs is observed for short timescales (t≲10^{3}). For longer timescales (t up to 10^{5}), the DBEs become localized. Such localization, depending on the DBE powers, has two different regimes: Regime I is the Anderson-like localization induced by the disorder, while Regime II is an enhanced localization attributed to both the disorder and discreteness. Our study is expected to shed light on understanding the interplay between disorder and strong nonlinearity, from the diffusive transport and localization properties of nonlinear localized excitations in random media.

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}.

On the Wave Turbulence Theory for the Nonlinear Schrödinger Equation with Random Potentials

We derive new kinetic and a porous medium equations from the nonlinear Schrödinger equation with random potentials. The kinetic equation has a very similar form compared to the four-wave turbulence kinetic equation in the wave turbulence theory. Moreover, we construct a class of self-similar solutions for the porous medium equation. These solutions spread with time, and this fact answers the “weak turbulence” question for the nonlinear Schrödinger equation with random potentials. We also derive Ohm’s law for the porous medium equation.

Multiple-time-scale analysis for pinned breathers in Bose-Hubbard chains

Localized and pinned discrete breathers in Bose-Einstein condensates in optical lattices or in arrays of optical waveguides oscillate with frequencies which are much higher than those present in the spectrum of the background. Hence, the interaction between localized breathers and their surroundings is extremely weak leading to a multiple-time-scale perturbation expansion. We identify the leading order in the asymptotic expansion of the breather amplitude which does not average to zero after one full oscillation. The reduced model predicts a lower bound of the breather drift times and explains the topological differences between breathers in dimers, trimers, and in spatially extended one-dimensional lattices even in the presence of transport from boundary heat-baths. These analytical boundaries hold true for lattices of any length, due to the highly localized nature of breathers.

Dynamical Freezing of Relaxation to Equilibrium

We provide evidence of an extremely slow thermalization occurring in the discrete nonlinear Schrödinger model. At variance with many similar processes encountered in statistical mechanics-typically ascribed to the presence of (free) energy barriers-here the slowness has a purely dynamical origin: it is due to the presence of an adiabatic invariant, which freezes the dynamics of a tall breather. Consequently, relaxation proceeds via rare events, where energy is suddenly released towards the background. We conjecture that this exponentially slow relaxation is a key ingredient contributing to the nonergodic behavior recently observed in the negative-temperature region of the discrete nonlinear Schrödinger equation.

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.

From Non-Normalizable Boltzmann-Gibbs Statistics to Infinite-Ergodic Theory

We study a particle immersed in a heat bath, in the presence of an external force which decays at least as rapidly as 1/x, e.g., a particle interacting with a surface through a Lennard-Jones or a logarithmic potential. As time increases, our system approaches a non-normalizable Boltzmann state. We study observables, such as the energy, which are integrable with respect to this asymptotic thermal state, calculating both time and ensemble averages. We derive a useful canonical-like ensemble which is defined out of equilibrium, using a maximum entropy principle, where the constraints are normalization, finite averaged energy, and a mean-squared displacement which increases linearly with time. Our work merges infinite-ergodic theory with Boltzmann-Gibbs statistics, thus extending the scope of the latter while shedding new light on the concept of ergodicity.