Researcher profile

Helen Christodoulidi

Helen Christodoulidi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
6works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2016arXiv

Dynamics and Statistics of the Fermi--Pasta--Ulam $β$--model with different ranges of particle interactions

In the present work we study the Fermi--Pasta--Ulam (FPU) $β$--model involving long--range interactions (LRI) in both the quadratic and quartic potentials, by introducing two independent exponents $α_1$ and $α_2$ respectively, which make the {forces decay} with distance $r$. Our results demonstrate that weak chaos, in the sense of decreasing Lyapunov exponents, and $q$--Gaussian probability density functions (pdfs) of sums of the momenta, occurs only when long--range interactions are included in the quartic part. More importantly, for $0\leq α_2<1$, we obtain extrapolated values for $q \equiv q_\infty >1$, as $N\rightarrow \infty$, suggesting that these pdfs persist in that limit. On the other hand, when long--range interactions are imposed only on the quadratic part, strong chaos and purely Gaussian pdfs are always obtained for the momenta. We have also focused on similar pdfs for the particle energies and have obtained $q_E$-exponentials (with $q_E>1$) when the quartic-term interactions are long--ranged, otherwise we get the standard Boltzmann-Gibbs weight, with $q=1$. The values of $q_E$ coincide, within small discrepancies, with the values of $q$ obtained by the momentum distributions.

preprint2015arXiv

Numerical integration of variational equations for Hamiltonian systems with long range interactions

We study numerically classical 1-dimensional Hamiltonian lattices involving inter-particle long range interactions that decay with distance like 1/r^alpha, for alpha>=0. We demonstrate that although such systems are generally characterized by strong chaos, they exhibit an unexpectedly organized behavior when the exponent alpha<1. This is shown by computing dynamical quantities such as the maximal Lyapunov exponent, which decreases as the number of degrees of freedom increases. We also discuss our numerical methods of symplectic integration implemented for the solution of the equations of motion together with their associated variational equations. The validity of our numerical simulations is estimated by showing that the total energy of the system is conserved within an accuracy of 4 digits (with integration step tau=0.02), even for as many as N=8000 particles and integration times as long as 10^6 units.

preprint2014arXiv

Fermi-Pasta-Ulam model with long-range interactions: Dynamics and thermostatistics

We introduce and numerically study a long-range-interaction generalization of the one-dimensional Fermi-Pasta-Ulam (FPU) $β-$ model. The standard quartic interaction is generalized through a coupling constant that decays as $1/r^α$ ($α\ge 0$)(with strength characterized by $b>0$). In the $α\to\infty$ limit we recover the original FPU model. Through classical molecular dynamics computations we show that (i) For $α\geq 1$ the maximal Lyapunov exponent remains finite and positive for increasing number of oscillators $N$ (thus yielding ergodicity), whereas, for $0 \le α<1$, it asymptotically decreases as $N^{- κ(α)}$ (consistent with violation of ergodicity); (ii) The distribution of time-averaged velocities is Maxwellian for $α$ large enough, whereas it is well approached by a $q$-Gaussian, with the index $q(α)$ monotonically decreasing from about 1.5 to 1 (Gaussian) when $α$ increases from zero to close to one. For $α$ small enough, the whole picture is consistent with a crossover at time $t_c$ from $q$-statistics to Boltzmann-Gibbs (BG) thermostatistics. More precisely, we construct a &#34;phase diagram&#34; for the system in which this crossover occurs through a frontier of the form $1/N \propto b^δ/t_c^γ$ with $γ>0$ and $δ>0$, in such a way that the $q=1$ ($q>1$) behavior dominates in the $\lim_{N \to\infty} \lim_{t \to\infty}$ ordering ($\lim_{t \to\infty} \lim_{N \to\infty}$ ordering).

preprint2013arXiv

Low-dimensional q-Tori in FPU Lattices: Dynamics and Localization Properties

This is a continuation of our study concerning q-tori, i.e. tori of low dimensionality in the phase space of nonlinear lattice models like the Fermi-Pasta-Ulam (FPU) model. In our previous work we focused on the beta FPU system, and we showed that the dynamical features of the q-tori serve as an interpretational tool to understand phenomena of energy localization in the FPU space of linear normal modes. In the present paper i) we employ the method of Poincare - Lindstedt series, for a fixed set of frequencies, in order to compute an explicit quasi-periodic representation of the trajectories lying on q-tori in the alpha model, and ii) we consider more general types of initial excitations in both the alpha and beta models. Furthermore we turn into questions of physical interest related to the dynamical features of the q-tori. We focus on particular q-tori solutions describing low-frequency `packets&#39; of modes, and excitations of a small set of modes with an arbitrary distribution in q-space. In the former case, we find formulae yielding an exponential profile of energy localization, following an analysis of the size of the leading order terms in the Poincare - Lindstedt series. In the latter case, we explain the observed localization patterns on the basis of a rigorous result concerning the propagation of non-zero terms in the Poincare - Lindstedt series from zeroth to subsequent orders. Finally, we discuss the extensive (i.e. independent of the number of degrees of freedom) properties of some q-tori solutions.

preprint2011arXiv

The two-stage dynamics in the Fermi-Pasta-Ulam problem: from regular to diffusive behavior

A numerical and analytical study of the relaxation to equilibrium of both the Fermi-Pasta-Ulam (FPU) alpha-model and the integrable Toda model, when the fundamental mode is initially excited, is reported. We show that the dynamics of both systems is almost identical on the short term, when the energies of the initially unexcited modes grow in geometric progression with time, through a secular avalanche process. At the end of this first stage of the dynamics the time-averaged modal energy spectrum of the Toda system stabilizes to its final profile, well described, at low energy, by the spectrum of a q-breather. The Toda equilibrium state is clearly shown to describe well the long-living quasi-state of the FPU system. On the long term, the modal energy spectrum of the FPU system slowly detaches from the Toda one by a diffusive-like rising of the tail modes, and eventually reaches the equilibrium flat shape. We find a simple law describing the growth of tail modes, which enables us to estimate the time-scale to equipartition of the FPU system, even when, at small energies, it becomes unobservable.

preprint2011arXiv

Weak chaos detection in the Fermi-Pasta-Ulam-$α$ system using $q$-Gaussian statistics

We study numerically statistical distributions of sums of orbit coordinates, viewed as independent random variables in the spirit of the Central Limit Theorem, in weakly chaotic regimes associated with the excitation of the first ($k=1$) and last ($k=N$) linear normal modes of the Fermi-Pasta-Ulam-$α$ system under fixed boundary conditions. We show that at low energies ($E=0.19$), when the $k=1$ linear mode is excited, chaotic diffusion occurs characterized by distributions that are well approximated for long times ($t>10^9$) by a $q$-Gaussian Quasi-Stationary State (QSS) with $q\approx1.4$. On the other hand, when the $k=N$ mode is excited at the same energy, diffusive phenomena are \textit{absent} and the motion is quasi-periodic. In fact, as the energy increases to $E=0.3$, the distributions in the former case pass through \textit{shorter} $q$-Gaussian states and tend rapidly to a Gaussian (i.e. $q\rightarrow 1$) where equipartition sets in, while in the latter we need to reach to E=4 to see a \textit{sudden transition} to Gaussian statistics, without any passage through an intermediate QSS. This may be explained by different energy localization properties and recurrence phenomena in the two cases, supporting the view that when the energy is placed in the first mode weak chaos and &#34;sticky&#34; dynamics lead to a more gradual process of energy sharing, while strong chaos and equipartition appear abruptly when only the last mode is initially excited.