Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userInteractions of renormalized waves in thermalized Fermi-Pasta-Ulam chains
72
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The dispersive interacting waves in Fermi-Pasta-Ulam (FPU) chains of particles in textit{thermal equilibrium} are studied from both statistical and wave resonance perspectives. It is shown that, even in a strongly nonlinear regime, the chain in thermal equilibrium can be effectively described by a system of weakly interacting textit{renormalized} nonlinear waves that possess (i) the Rayleigh-Jeans distribution and (ii) zero correlations between waves, just as noninteracting free waves would. This renormalization is achieved through a set of canonical transformations. The renormalized linear dispersion of these renormalized waves is obtained and shown to be in excellent agreement with numerical experiments. Moreover, a dynamical interpretation of the renormalization of the dispersion relation is provided via a self-consistency, mean-field argument. It turns out that this renormalization arises mainly from the trivial resonant wave interactions, i.e., interactions with no momentum exchange. Furthermore, using a multiple time-scale, statistical averaging method, we show that the interactions of near-resonant waves give rise to the broadening of the resonance peaks in the frequency spectrum of renormalized modes. The theoretical prediction for the resonance width for the thermalized $beta$-FPU chain is found to be in very good agreement with its numerically measured value.
rate research
Read More
In systems of N coupled anharmonic oscillators, exact resonant interactions play an important role in the energy exchange between normal modes. In the weakly nonlinear regime, those interactions may facilitate energy equipartition in Fourier space. We consider analytically resonant wave-wave interactions for the celebrated Fermi-Pasta-Ulam-Tsingou (FPUT) system. Using a number-theoretical approach based on cyclotomic polynomials, we show that the problem of finding exact resonances for a system of N particles is equivalent to a Diophantine equation whose solutions depend sensitively on the set of divisors of N. We provide an algorithm to construct all possible resonances, based on two methods: pairing-off and cyclotomic, which we introduce to build up explicit solutions to the 4-, 5- and 6-wave resonant conditions. Our results shed some light in the understanding of the long-standing FPUT paradox, regarding the sensitivity of the resonant manifolds with respect to the number of particles N and the corresponding time scale of the interactions leading to thermalisation. In this light we demonstrate that 6-wave resonances always exist for any N, while 5-wave resonances exist if N is divisible by 3 and N > 6. It is known (for finite N) that 4-wave resonances do not mix energy across the spectrum, so we investigate whether 5-wave resonances can produce energy mixing across a significant region of the Fourier spectrum by analysing the interconnected network of Fourier modes that can interact nonlinearly via resonances. The answer depends on the set of odd divisors of N that are not divisible by 3: the size of this set determines the number of dynamically independent components, corresponding to independent constants of motion (energies). We show that 6-wave resonances connect all these independent components, providing in principle a restoring mechanism for full-scale thermalisation.
We study the original $alpha$-Fermi-Pasta-Ulam (FPU) system with $N=16,32$ and $64$ masses connected by a nonlinear quadratic spring. Our approach is based on resonant wave-wave interaction theory, i.e. we assume that, in the weakly nonlinear regime (the one in which Fermi was originally interested), the large time dynamics is ruled by exact resonances. After a detailed analysis of the $alpha$-FPU equation of motion, we find that the first non trivial resonances correspond to six-wave interactions. Those are precisely the interactions responsible for the thermalization of the energy in the spectrum. We predict that for small amplitude random waves the time scale of such interactions is extremely large and it is of the order of $1/epsilon^8$, where $epsilon$ is the small parameter in the system. The wave-wave interaction theory is not based on any threshold: equipartition is predicted for arbitrary small nonlinearity. Our results are supported by extensive numerical simulations. A key role in our finding is played by the {it Umklapp} (flip over) resonant interactions, typical of discrete systems. The thermodynamic limit is also briefly discussed.
We perform a thorough investigation of the first FPUT recurrence in the $beta$-FPUT chain for both positive and negative $beta$. We show numerically that the rescaled FPUT recurrence time $T_{r}=t_{r}/(N+1)^{3}$ depends, for large $N$, only on the parameter $Sequiv Ebeta(N+1)$. Our numerics also reveal that for small $left|Sright|$, $T_{r}$ is linear in $S$ with positive slope for both positive and negative $beta$. For large $left|Sright|$, $T_{r}$ is proportional to $left|Sright|^{-1/2}$ for both positive and negative $beta$ but with different multiplicative constants. In the continuum limit, the $beta$-FPUT chain approaches the modified Korteweg-de Vries (mKdV) equation, which we investigate numerically to better understand the FPUT recurrences on the lattice. In the continuum, the recurrence time closely follows the $|S|^{-1/2}$ scaling and can be interpreted in terms of solitons, as in the case of the KdV equation for the $alpha$ chain. The difference in the multiplicative factors between positive and negative $beta$ arises from soliton-kink interactions which exist only in the negative $beta$ case. We complement our numerical results with analytical considerations in the nearly linear regime (small $left|Sright|$) and in the highly nonlinear regime (large $left|Sright|$). For the former, we extend previous results using a shifted-frequency perturbation theory and find a closed form for $T_{r}$ which depends only on $S$. In the latter regime, we show that $T_{r}proptoleft| Sright|^{-1/2}$ is predicted by the soliton theory in the continuum limit. We end by discussing the striking differences in the amount of energy mixing as well as the existence of the FPUT recurrences between positive and negative $beta$ and offer some remarks on the thermodynamic limit.
We consider a $pi$-mode solution of the Fermi-Pasta-Ulam $beta$ system. By perturbing it, we study the system as a function of the energy density from a regime where the solution is stable to a regime, where is unstable, first weakly and then strongly chaotic. We introduce, as indicator of stochasticity, the ratio $rho$ (when is defined) between the second and the first moment of a given probability distribution. We will show numerically that the transition between weak and strong chaos can be interpreted as the symmetry breaking of a set of suitable dynamical variables. Moreover, we show that in the region of weak chaos there is numerical evidence that the thermostatistic is governed by the Tsallis distribution.
The recently introduced GALI method is used for rapidly detecting chaos, determining the dimensionality of regular motion and predicting slow diffusion in multi--dimensional Hamiltonian systems. We propose an efficient computation of the GALI$_k$ indices, which represent volume elements of $k$ randomly chosen deviation vectors from a given orbit, based on the Singular Value Decomposition (SVD) algorithm. We obtain theoretically and verify numerically asymptotic estimates of GALIs long--time behavior in the case of regular orbits lying on low--dimensional tori. The GALI$_k$ indices are applied to rapidly detect chaotic oscillations, identify low--dimensional tori of Fermi--Pasta--Ulam (FPU) lattices at low energies and predict weak diffusion away from quasiperiodic motion, long before it is actually observed in the oscillations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
