Dynamical energy equipartition of the Toda model with additional on-site potentials

We study the dynamical energy equipartition properties in the integrable Toda model with additional uniform or disordered on-site energies by extensive numerical simulations. The total energy is initially equidistributed among some of the lowest frequency linear modes. For the Toda model with uniform on-site potentials, the energy spectrum keeps its profile nearly unchanged in a relatively short time scale. On a much longer time scale, the energies of tail modes increase slowly with time. Energy equipartition is far away from being attached in our studied time scale. For the Toda model with disordered on-site potentials, the energy transfers continuously to the high frequency modes and eventually towards energy equipartition. We further perform a systematic study of the equipartition time teq depending on the energy density εand the nonlinear parameter α in the thermodynamic limit for the Toda model with disordered on-site potentials. We find teq∝(1/ε)~a(1/α)~b, where b ≈ 2 a. The values of a and b are increased when increasing the strengths of disordered on-site potentials or decreasing the number of initially excited modes.     

The Fermi-Pasta-Ulam Problem and Its Underlying Integrable Dynamics

This paper is devoted to a numerical study of the familiar α+β FPU model. Precisely, we here discuss, revisit and combine together two main ideas on the subject: (i) In the system, at small specific energy ε=E/N, two well separated time-scales are present: in the former one a kind of metastable state is produced, while in the second much larger one, such an intermediate state evolves and reaches statistical equilibrium. (ii) FPU should be interpreted as a perturbed Toda model, rather than (as is typical) as a linear model perturbed by nonlinear terms. In the view we here present and support, the former time scale is the one in which FPU is essentially integrable, its dynamics being almost indistinguishable from the Toda dynamics: the Toda actions stay constant for FPU too (while the usual linear normal modes do not), the angles fill their almost invariant torus, and nothing else happens. The second time scale is instead the one in which the Toda actions significantly evolve, and statistical equilibrium is possible. We study both FPU-like initial states, in which only a few degrees of freedom are excited, and generic initial states extracted randomly from an (approximated) microcanonical distribution. The study is based on a close comparison between the behavior of FPU and Toda in various situations. The main technical novelty is the study of the correlation functions of the Toda constants of motion in the FPU dynamics; such a study allows us to provide a good definition of the equilibrium time τ, i.e. of the second time scale, for generic initial data. Our investigation shows that τ is stable in the thermodynamic limit, i.e. the limit of large N at fixed ε, and that by reducing ε (ideally, the temperature), τ approximately grows following a power law τ～ε ^?a , with a=5/2.     

Nonlinear Dynamics and Fluctuations of Dissipative Toda Chains

The dynamics of a ring of masses including dissipative forces (passive or active friction) and Toda interactions between the masses is investigated. The characteristic attractor structure and the influence of noise by coupling to a heat bath are studied. The system may be driven from the thermodynamic equilibrium to far from equilibrium states by including negative friction. We show, that over-critical pumping with free energy may lead to a partition of the phase space into attractor regions corresponding to several types of collective motions including uniform rotations, one- and multiple soliton-like excitations and relative oscillations. The distribution functions in the phase space and the correlation functions of the forces and the spectra of nonlinear excitations are calculated. We show that a finite-size Toda ring with weak thermal coupling develops at intermediate temperatures a broadband colored noise spectrum with an 1/f tail at low frequencies.     

Modulational instability in isolated and driven Fermi–Pasta–Ulam lattices

We present a detailed analysis of the modulational instability of the zone-boundary mode for one and higher-dimensional Fermi–Pasta–Ulam (FPU) lattices. The growth of the instability is followed by a process of relaxation to equipartition, which we have called the Anti-FPU problem because the energy is initially fed into the highest frequency part of the spectrum, while in the original FPU problem low frequency excitations of the lattice were considered. This relaxation process leads to the formation of chaotic breathers in both one and two space dimensions. The system then relaxes to energy equipartition, on time scales that increase as the energy density is decreased. We supplement this study by considering the nonconservative case, where the FPU lattice is homogeneously driven at high frequencies. Standing and travelling nonlinear waves and solitonic patterns are detected in this case. Finally we investigate the dynamics of the FPU chain when one end is driven at a frequency located above the zone boundary. We show that this excitation stimulates nonlinear bandgap transmission effects.     

On Metastability in FPU

We present an analytical study of the Fermi–Pasta–Ulam (FPU) α–model with periodic boundary conditions. We analyze the dynamics corresponding to initial data with one low frequency Fourier mode excited. We show that, correspondingly, a pair of KdV equations constitute the resonant normal form of the system. We also use such a normal form in order to prove the existence of a metastability phenomenon. More precisely, we show that the time average of the modal energy spectrum rapidly attains a well defined distribution corresponding to a packet of low frequencies modes. Subsequently, the distribution remains unchanged up to the time scales of validity of our approximation. The phenomenon is controlled by the specific energy.     

The anti-FPU problem

We present a detailed analysis of the modulational instability of the zone-boundary mode for one and higher-dimensional Fermi-Pasta-Ulam (FPU) lattices. Following this instability, a process of relaxation to equipartition takes place, which we have called the Anti-FPU problem because the energy is initially fed into the highest frequency part of the spectrum, at variance with the original FPU problem (low frequency excitations of the lattice). This process leads to the formation of chaotic breathers in both one and two dimensions. Finally, the system relaxes to energy equipartition on time scales which increase as the energy density is decreased. We show that breathers formed when cooling the lattice at the edges, starting from a random initial state, bear strong qualitative similarities with chaotic breathers.     

Tail resonances of Fermi-Pasta-Ulam q-breathers and their impact on the pathway to equipartition

Upon initial excitation of a few normal modes the energy distribution among all modes of a nonlinear atomic chain (the Fermi-Pasta-Ulam model) exhibits exponential localization on large time scales. At the same time, resonant anomalies (peaks) are observed in its weakly excited tail for long times preceding equipartition. We observe a similar resonant tail structure also for exact time-periodic Lyapunov orbits, coined q-breathers due to their exponential localization in modal space. We give a simple explanation for this structure in terms of superharmonic resonances. The resonance analysis agrees very well with numerical results and has predictive power. We extend a previously developed perturbation method, based essentially on a Poincare-Lindstedt scheme, in order to account for these resonances, and in order to treat more general model cases, including truncated Toda potentials. Our results give a qualitative and semiquantitative account for the superharmonic resonances of q-breathers and natural packets.     

Classical dissipation and asymptotic equilibrium via interaction with chaotic systems

We study the energy flow between a one-dimensional oscillator and a chaotic system with two degrees of freedom in the weak coupling limit. The oscillator's observables are averaged over an initially microcanonical ensemble of trajectories of the chaotic system, which plays the role of an environment for the oscillator. We show numerically that the oscillator's average energy exhibits irreversible dynamics and ‘thermal’ equilibrium at long times. We use linear response theory to describe the dynamics at short times and we derive a condition for the absorption or dissipation of energy by the oscillator from the chaotic system. The equilibrium properties at long times, including the average equilibrium energies and the energy distributions, are explained with the help of statistical arguments. We also check that the concept of temperature defined in terms of the ‘volume entropy’ agrees very well with these energy distributions.     

Asymptotic Stability of Lattice Solitons in the Energy Space

Orbital and asymptotic stability for 1-soliton solutions of the Toda lattice equations as well as for small solitary waves of the FPU lattice equations are established in the energy space. Unlike analogous Hamiltonian PDEs, the lattice equations do not conserve the adjoint momentum. In fact, the Toda lattice equation is a bidirectional model that does not fit in with the existing theory for the Hamiltonian systems by Grillakis, Shatah and Strauss. To prove stability of 1-soliton solutions, we split a solution around a 1-soliton into a small solution that moves more slowly than the main solitary wave and an exponentially localized part. We apply a decay estimate for solutions to a linearized Toda equation which has been recently proved by Mizumachi and Pego to estimate the localized part. We improve the asymptotic stability results for FPU lattices in a weighted space obtained by Friesecke and Pego.     

A study of the Fermi-Pasta-Ulam problem in dimension two

Continuing the previous work on the same subject, we study here different two-dimensional Fermi-Pasta-Ulam (FPU)-like models, namely, planar models with a triangular cell, molecular-type potential and different boundary conditions, and perform on them both traditional FPU-like numerical experiments, i.e., experiments in which energy is initially concentrated on a small subset of normal modes, and other experiments, in which we test the time scale for the decay of a large fluctuation when all modes are excited almost to the same extent. For each experiment, we observe the behavior of the different two-dimensional systems and also make an accurate comparison with the behavior of a one-dimensional model with an identical potential. We assume the thermodynamic point of view and try to understand the behavior of the system for large n (the number of degrees of freedom) at fixed specific energy epsilon=En. As a result, it turns out that: (i) The difference between dimension one and two, if n is large, is substantial. In particular (making reference to FPU-like initial conditions) the "one-dimensional scenario," in which the dynamics is dominated for a long time scale by a weakly chaotic metastable situation, in dimension two is absent; moreover in dimension two, for large n, the time scale for energy sharing among normal modes is drastically shorter than in dimension one. (ii) The boundary conditions in dimension two play a relevant role. Indeed, models with fixed or open boundary conditions give fast equipartition, on a rather short time scale of order epsilon(-1), while a periodic model gives longer equilibrium times (although much shorter than in dimension one).     

Studies of thermal conductivity in Fermi-Pasta-Ulam-like lattices

The pioneering computer simulations of the energy relaxation mechanisms performed by Fermi, Pasta, and Ulam (FPU) can be considered as the first attempt of understanding energy relaxation and thus heat conduction in lattices of nonlinear oscillators. In this paper we describe the most recent achievements about the divergence of heat conductivity with the system size in one-dimensional (1D) and two-dimensional FPU-like lattices. The anomalous behavior is particularly evident at low energies, where it is enhanced by the quasiharmonic character of the lattice dynamics. Remarkably, anomalies persist also in the strongly chaotic region where long-time tails develop in the current autocorrelation function. A modal analysis of the 1D case is also presented in order to gain further insight about the role played by boundary conditions.     

Time-Scales to Equipartition in the Fermi–Pasta–Ulam Problem: Finite-Size Effects and Thermodynamic Limit

We investigate numerically the common α+β and the pure β FPU models, as well as some higher order generalizations. We consider initial conditions in which only low-frequency normal modes are excited, and perform a very accurate systematic study of the equilibrium time as a function of the number N of particles, the specific energy ε, and the parameters α and β. While at any fixed N the equilibrium time is found to be a stretched exponential in 1/ε, in the thermodynamic limit, i.e. for N→∞ at fixed ε, we observe a crossover to a power law. Concerning the (usually disregarded) dependence of T _eq on α and β, we find it is nontrivial, and propose and test a general law. A central role is played by the comparison of the FPU models with the Toda model.     

Strongly nonlinear beats in the dynamics of an elastic system with a strong local stiffness nonlinearity: Analysis and identification

We consider a linear cantilever beam attached to ground through a strongly nonlinear stiffness at its free boundary, and study its dynamics computationally by the assumed-modes method. The nonlinear stiffness of this system has no linear component, so it is essentially nonlinear and nonlinearizable. We find that the strong nonlinearity mostly affects the lower-frequency bending modes and gives rise to strongly nonlinear beat phenomena. Analysis of these beats proves that they are caused by internal resonance interactions of nonlinear normal modes (NNMs) of the system. These internal resonances are not of the classical type since they occur between bending modes whose linearized natural frequencies are not necessarily related by rational ratios; rather, they are due to the strong energy-dependence of the frequency of oscillation of the corresponding NNMs of the beam (arising from the strong local stiffness nonlinearity) and occur at energy ranges where the frequencies of these NNMs are rationally related. Nonlinear effects start at a different energy level for each mode. Lower modes are influenced at lower energies due to larger modal displacements than higher modes and thus, at certain energy levels, the NNMs become rationally related, which results in internal resonance. The internal resonances of NNMs are studied using a reduced order model of the beam system. Then, a nonlinear system identification method is developed, capable of identifying this type of strongly nonlinear modal interactions. It is based on an adaptive step-by-step application of empirical mode decomposition (EMD) to the measured time series, which makes it valid for multi-frequency beating signals. Our work extends an earlier nonlinear system identification approach developed for nearly mono-frequency (monochromatic) signals. The extended system identification method is applied to the identification of the strongly nonlinear dynamics of the considered cantilever beam with the local strong nonlinear stiffness at its free end.     

Phonon relaxation and heat conduction in one-dimensional Fermiben Pastaben Ulam β lattices by molecular dynamics simulations

The phonon relaxation and heat conduction in one-dimensional Fermi-Pasta-Ulam (FPU) β lattices are studied by using molecular dynamics simulations. The phonon relaxation rate, which dominates the length dependence of the FPU β lattice, is first calculated from the energy autocorrelation function for different modes at various temperatures through equilibrium molecular dynamics simulations. We find that the relaxation rate as a function of wave number k is proportional to k^1.688, which leads to a N^0.41 divergence of the thermal conductivity in the framework of Green-Kubo relation. This is also in good agreement with the data obtained by non-equilibrium molecular dynamics simulations which estimate the length dependence exponent of the thermal conductivity as 0.415. Our results confirm the N^2/5 divergence in one-dimensional FPU β lattices. The effects of the heat flux on the thermal conductivity are also studied by imposing different temperature differences on the two ends of the lattices. We find that the thermal conductivity is insensitive to the heat flux under our simulation conditions. It implies that the linear response theory is applicable towards the heat conduction in one-dimensional FPU β lattices.     

Time scale for energy equipartition in a two-dimensional FPU model

The FPU problem, i.e., the problem of energy equipartition among normal modes in a weakly nonlinear lattice, is here studied in dimension two, more precisely in a model with triangular cell and nearest-neighbors Lennard-Jones interaction. The number n of degrees of freedom ranges from 182 to 6338. Energy is initially equidistributed among a small number n(0) of low frequency modes, with n(0) proportional to n. We study numerically the time evolution of the so-called spectral entropy and the related "effective number" n(eff) of degrees of freedom involved in the dynamics; in this (rather typical) way we can estimate, for each n and each specific energy (energy per degree of freedom) epsilon, the time scale T(n)(epsilon) for energy equipartition. Numerical results indicate that in the thermodynamic limit the equipartition times are short: more precisely, for large n at fixed epsilon we find a limit curve T(infinity)(epsilon), and T(infinity) grows only as epsilon(-1) for small epsilon. Larger equipartition times are obtained by lowering epsilon, at fixed n, below a crossover value epsilon(c)(n). However, epsilon(c) appears to vanish by increasing n (faster than 1n), and the total energy E=nepsilon, rather than epsilon, appears to be the relevant variable when n is large and epsilon相似文献   

q-Breathers and the Fermi-Pasta-Ulam problem

The Fermi-Pasta-Ulam (FPU) paradox consists of the non-equipartition of energy among normal modes of a weakly anharmonic atomic chain model. In the harmonic limit each normal mode corresponds to a periodic orbit in phase space and is characterized by its wave number q. We continue normal modes from the harmonic limit into the FPU parameter regime and obtain persistence of these periodic orbits, termed here q-breathers (QB). They are characterized by time periodicity, exponential localization in the q-space of normal modes and linear stability up to a size-dependent threshold amplitude. Trajectories computed in the original FPU setting are perturbations around these exact QB solutions. The QB concept is applicable to other nonlinear lattices as well.     

Topological phase transition and electrically tunable diamagnetism in silicene

The movement and relaxation of the localized energy on FPU lattices have been studied by using Wavelet transforms methods. The energy relaxation mechanism for nonlinear chains involves the degradation of higher frequency excitations into lower frequencies. It is shown that low frequency modes decay more slowly in nonlinear chains. The wavelet spectrum exhibits a behavior involving the interplay of phonon modes and breather modes.     

Analysis of the modal energy distribution of an excited vibrating panel coupled with a heavy fluid cavity by a dual modal formulation

This paper describes the modal interaction between a panel and a heavy fluid cavity when the panel is excited by a broad band force in a given frequency band. The dual modal formulation (DMF) allows describing the fluid–structure coupling using the modes of each uncoupled subsystem. After having studied the convergence of the modal series on a test case, we estimate the modal energies and the total energy of each subsystem. An analysis of modal energy distribution is performed. It allows us to study the validity of SEA assumptions for this case. Added mass and added stiffness effects of the fluid are observed. These effects are related to the non-resonant acoustic modes below and above the frequency band of excitation. Moreover, the role of the spatial coupling of the resonant cavity modes with the non-resonant structure modes is also highlighted. As a result, the energy transmitted between the structure and the heavy fluid cavity generally cannot be deduced from the SEA relation established for a light fluid cavity.     

Full Fermi–Pasta–Ulam Recurrence Interaction with Fluctuations of the Medium

It is shown in a numerical study that the simultaneous influence of discrete and continuous random processes on the full Fermi–Pasta–Ulam (FPU) recurrence dynamics in two parametrically coupled identical chains of vibrators with open ends and under different initial conditions leads to stabilization of the FPU spectrum. A greater influence on chains of gaussian noise as compared to discrete noise is revealed. An increase in the amplitudes of both noises by one order of magnitude results in considerable parameter changes in their FPU spectrum. The full FPU recurrence interaction in coupled chains with discrete and continuous random noise of the medium manifests itself in extremely low-frequency periodic processes in stochastic dynamics of the medium (with frequency two orders of magnitude lower than the lowest frequency in the initial conditions of the chains).     

Double resonance in Raman scattering on single phonon modes of GaSe

Resonance in the Raman cross section for one-phonon modes of GaSe at 80°K has been observed for both incident and scattered photon energies equal to the direct exciton energy. The Raman efficiency spectrum within the region of resonance is shown to fit well with the shape computed from existing theories.