Source author record

Ch. Skokos

Ch. Skokos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

22 published item(s)

preprint2020arXiv

Properties of normal modes in a modified disordered Klein-Gordon lattice: From disorder to order

We introduce a modified version of the disordered Klein-Gordon lattice model, having two parameters for controlling the disorder strength: $D$, which determines the range of the coefficients of the on-site potentials, and $W$, which defines the strength of the nearest-neighbor interactions. We fix $W=4$ and investigate how the properties of the system's normal modes change as we approach its ordered version, i.e. $D\rightarrow 0$. We show that the probability density distribution of the normal modes' frequencies takes a `U'-shaped profile as $D$ decreases. Furthermore, we use two quantities for estimating the modes' spatial extent, the so-called localization volume $V$ (which is related to the mode's second moment) and the mode's participation number $P$. We show that both quantities scale as $\propto D^{-2}$ when $D$ approaches zero and we numerically verify a proportionality relation between them as $V/P \approx 2.6$.

preprint2020arXiv

Wave propagation in a strongly disordered 1D phononic lattice supporting rotational waves

We investigate the dynamical properties of a strongly disordered micropolar lattice made up of cubic block units. This phononic lattice model supports both transverse and rotational degrees of freedom hence its disordered variant posses an interesting problem as it can be used to model physically important systems like beam-like microstructures. Different kinds of single site excitations (momentum or displacement) on the two degrees of freedom are found to lead to different energy transport both superdiffusive and subdiffusive. We show that the energy spreading is facilitated both by the low frequency extended waves and a set of high frequency modes located at the edge of the upper branch of the periodic case for any initial condition. However, the second moment of the energy distribution strongly depends on the initial condition and it is slower than the underlying one dimensional harmonic lattice (with one degree of freedom). Finally, a limiting case of the micropolar lattice is studied where Anderson localization is found to persist and no energy spreading takes place.

preprint2015arXiv

Energy transport in one-dimensional disordered granular solids

We investigate the energy transport in one-dimensional disordered granular solids by extensive numerical simulations. In particular, we consider the case of a polydisperse granular chain composed of spherical beads of the same material and with radii taken from a random distribution. We start by examining the linear case, in which it is known that the energy transport strongly depends on the type of initial conditions. Thus, we consider two sets of initial conditions: i) an initial displacement and ii) an initial momentum excitation of a single bead. After establishing the regime of sufficiently strong disorder, we focus our studies on the role of nonlinearity for both sets of initial conditions. By increasing the initial excitation amplitudes we are able to identify three distinct dynamical regimes with different energy transport properties: a near linear, a weakly nonlinear and a highly nonlinear regime. Although energy spreading is found to be increasing for higher nonlinearities, in the weakly nonlinear regime no clear asymptotic behavior of the spreading is found. In this regime, we additionally find that energy, initially trapped in a localized region, can be eventually detrapped and this has a direct influence on the fluctuations of the energy spreading. We also demonstrate that in the highly nonlinear regime, the differences in energy transport between the two sets of initial conditions vanish. Actually, in this regime the energy is almost ballistically transported through shock-like excitations.

preprint2015arXiv

The Smaller (SALI) and the Generalized (GALI) Alignment Indices: Efficient Methods of Chaos Detection

We provide a concise presentation of the Smaller (SALI) and the Generalized Alignment Index (GALI) methods of chaos detection. These are efficient chaos indicators based on the evolution of two or more, initially distinct, deviation vectors from the studied orbit. After explaining the motivation behind the introduction of these indices, we sum up the behaviors they exhibit for regular and chaotic motion, as well as for stable and unstable periodic orbits, focusing mainly on finite-dimensional conservative systems: autonomous Hamiltonian models and symplectic maps. We emphasize the advantages of these methods in studying the global dynamics of a system, as well as their ability to identify regular motion on low dimensional tori. Finally we discuss several applications of these indices to problems originating from different scientific fields like celestial mechanics, galactic dynamics, accelerator physics and condensed matter physics.

preprint2014arXiv

Complex Statistics and Diffusion in Nonlinear Disordered Particle Chains

We investigate dynamically and statistically diffusive motion in a Klein-Gordon particle chain in the presence of disorder. In particular, we examine a low energy (subdiffusive) and a higher energy (self-trapping) case and verify that subdiffusive spreading is always observed. We then carry out a statistical analysis of the motion in both cases in the sense of the Central Limit Theorem and present evidence of different chaos behaviors, for various groups of particles. Integrating the equations of motion for times as long as $10^9$, our probability distribution functions always tend to Gaussians and show that the dynamics does not relax onto a quasi-periodic KAM torus and that diffusion continues to spread chaotically for arbitrarily long times.

preprint2014arXiv

Delocalization of wave packets in disordered nonlinear chains

We consider the spatiotemporal evolution of a wave packet in disordered nonlinear Schrödinger and anharmonic oscillator chains. In the absence of nonlinearity all eigenstates are spatially localized with an upper bound on the localization length (Anderson localization). Nonlinear terms in the equations of motion destroy Anderson localization due to nonintegrability and deterministic chaos. At least a finite part of an initially localized wave packet will subdiffusively spread without limits. We analyze the details of this spreading process. We compare the evolution of single site, single mode and general finite size excitations, and study the statistics of detrapping times. We investigate the properties of mode-mode resonances, which are responsible for the incoherent delocalization process.

preprint2014arXiv

High order three part split symplectic integrators: Efficient techniques for the long time simulation of the disordered discrete nonlinear Schroedinger equation

While symplectic integration methods based on operator splitting are well established in many branches of science, high order methods for Hamiltonian systems that split in more than two parts have not been studied in great detail. Here, we present several high order symplectic integrators for Hamiltonian systems that can be split in exactly three integrable parts. We apply these techniques, as a practical case, for the integration of the disordered, discrete nonlinear Schroedinger equation (DDNLS) and compare their efficiencies. Three part split algorithms provide effective means to numerically study the asymptotic behavior of wave packet spreading in the DDNLS - a hotly debated subject in current scientific literature.

preprint2013arXiv

Chaoticity without thermalisation in disordered lattices

We study chaoticity and thermalization in Bose-Einstein condensates in disordered lattices, described by the discrete nonlinear Schrödinger equation (DNLS). A symplectic integration method allows us to accurately obtain both the full phase space trajectories and their maximum Lyapunov exponents (mLEs), which characterize their chaoticity. We find that disorder destroys ergodicity by breaking up phase space into subsystems that are effectively disjoint on experimentally relevant timescales, even though energetically, classical localisation cannot occur. This leads us to conclude that the mLE is a very poor ergodicity indicator, since it is not sensitive to the trajectory being confined to a subregion of phase space. The eventual thermalization of a BEC in a disordered lattice cannot be predicted based only on the chaoticity of its phase space trajectory.

preprint2012arXiv

Interplay Between Chaotic and Regular Motion in a Time-Dependent Barred Galaxy Model

We study the distinction and quantification of chaotic and regular motion in a time-dependent Hamiltonian barred galaxy model. Recently, a strong correlation was found between the strength of the bar and the presence of chaotic motion in this system, as models with relatively strong bars were shown to exhibit stronger chaotic behavior compared to those having a weaker bar component. Here, we attempt to further explore this connection by studying the interplay between chaotic and regular behavior of star orbits when the parameters of the model evolve in time. This happens for example when one introduces linear time dependence in the mass parameters of the model to mimic, in some general sense, the effect of self-consistent interactions of the actual N-body problem. We thus observe, in this simple time-dependent model also, that the increase of the bar's mass leads to an increase of the system's chaoticity. We propose a new way of using the Generalized Alignment Index (GALI) method as a reliable criterion to estimate the relative fraction of chaotic vs. regular orbits in such time-dependent potentials, which proves to be much more efficient than the computation of Lyapunov exponents. In particular, GALI is able to capture subtle changes in the nature of an orbit (or ensemble of orbits) even for relatively small time intervals, which makes it ideal for detecting dynamical transitions in time-dependent systems.

preprint2011arXiv

Comparing the efficiency of numerical techniques for the integration of variational equations

We present a comparison of different numerical techniques for the integration of variational equations. The methods presented can be applied to any autonomous Hamiltonian system whose kinetic energy is quadratic in the generalized momenta, and whose potential is a function of the generalized positions. We apply the various techniques to the well-known Hénon-Heiles system, and use the Smaller Alignment Index (SALI) method of chaos detection to evaluate the percentage of its chaotic orbits. The accuracy and the speed of the integration schemes in evaluating this percentage are used to investigate the numerical efficiency of the various techniques.

preprint2011arXiv

Improving beam stability in particle accelerator models by using Hamiltonian control

We derive a Hamiltonian control theory which can be applied to a 4D symplectic map that models a ring particle accelerator composed of elements with sextupole nonlinearity. The controlled system is designed to exhibit a more regular orbital behavior than the uncontrolled one. Using the Smaller Alignement Index (SALI) chaos indicator, we are able to show that the controlled system has a dynamical aperture up to 1.7 times larger than the original mode

preprint2011arXiv

Nonlinear waves in disordered chains: probing the limits of chaos and spreading

We probe the limits of nonlinear wave spreading in disordered chains which are known to localize linear waves. We particularly extend recent studies on the regimes of strong and weak chaos during subdiffusive spreading of wave packets [EPL {\bf 91}, 30001 (2010)] and consider strong disorder, which favors Anderson localization. We probe the limit of infinite disorder strength and study Fröhlich-Spencer-Wayne models. We find that the assumption of chaotic wave packet dynamics and its impact on spreading is in accord with all studied cases. Spreading appears to be asymptotic, without any observable slowing down. We also consider chains with spatially inhomogeneous nonlinearity which give further support to our findings and conclusions.

preprint2011arXiv

Probing the local dynamics of periodic orbits by the generalized alignment index (GALI) method

As originally formulated, the Generalized Alignment Index (GALI) method of chaos detection has so far been applied to distinguish quasiperiodic from chaotic motion in conservative nonlinear dynamical systems. In this paper we extend its realm of applicability by using it to investigate the local dynamics of periodic orbits. We show theoretically and verify numerically that for stable periodic orbits the GALIs tend to zero following particular power laws for Hamiltonian flows, while they fluctuate around non-zero values for symplectic maps. By comparison, the GALIs of unstable periodic orbits tend exponentially to zero, both for flows and maps. We also apply the GALIs for investigating the dynamics in the neighborhood of periodic orbits, and show that for chaotic solutions influenced by the homoclinic tangle of unstable periodic orbits, the GALIs can exhibit a remarkable oscillatory behavior during which their amplitudes change by many orders of magnitude. Finally, we use the GALI method to elucidate further the connection between the dynamics of Hamiltonian flows and symplectic maps. In particular, we show that, using for the computation of GALIs the components of deviation vectors orthogonal to the direction of motion, the indices of stable periodic orbits behave for flows as they do for maps.

preprint2011arXiv

Wave interactions in localizing media - a coin with many faces

A variety of heterogeneous potentials are capable of localizing linear non-interacting waves. In this work, we review different examples of heterogeneous localizing potentials which were realized in experiments. We then discuss the impact of nonlinearity induced by wave interactions, in particular its destructive effect on the localizing properties of the heterogeneous potentials.

preprint2010arXiv

Numerical integration of variational equations

We present and compare different numerical schemes for the integration of the variational equations of autonomous Hamiltonian systems whose kinetic energy is quadratic in the generalized momenta and whose potential is a function of the generalized positions. We apply these techniques to Hamiltonian systems of various degrees of freedom, and investigate their efficiency in accurately reproducing well-known properties of chaos indicators like the Lyapunov Characteristic Exponents (LCEs) and the Generalized Alignment Indices (GALIs). We find that the best numerical performance is exhibited by the \textit{`tangent map (TM) method'}, a scheme based on symplectic integration techniques which proves to be optimal in speed and accuracy. According to this method, a symplectic integrator is used to approximate the solution of the Hamilton's equations of motion by the repeated action of a symplectic map $S$, while the corresponding tangent map $TS$, is used for the integration of the variational equations. A simple and systematic technique to construct $TS$ is also presented.

preprint2010arXiv

SALI: An Efficient Indicator of Chaos with Application to 2 and 3 Degrees of Freedom Hamiltonian Systems

The Smaller Alignment Index (SALI) is a very useful and efficient indicator that can distinguish rapidly and with certainty between ordered and chaotic motion in Hamiltonian systems. This is based on the different behavior of the SALI in the two cases, as it fluctuates around non-zero values for ordered orbits and converges exponentially to zero for chaotic orbits. In this paper we present a detailed numerical study, comparing the advantages of the SALI to those of the Maximal Lyapunov Exponent (MLE), which is the traditional indicator of chaos detection in dynamical systems. Exploiting these advantages, we demonstrate how the SALI can be used to identify even tiny regions of order and chaos in Hamiltonian systems of 2 and 3 degrees of freedom.

preprint2010arXiv

Spreading of wave packets in disordered systems with tunable nonlinearity

We study the spreading of single-site excitations in one-dimensional disordered Klein-Gordon chains with tunable nonlinearity $|u_{l}|^σ u_{l}$ for different values of $σ$. We perform extensive numerical simulations where wave packets are evolved a) without and, b) with dephasing in normal mode space. Subdiffusive spreading is observed with the second moment of wave packets growing as $t^α$. The dependence of the numerically computed exponent $α$ on $σ$ is in very good agreement with our theoretical predictions both for the evolution of the wave packet with and without dephasing (for $σ\geq 2$ in the latter case). We discuss evidence of the existence of a regime of strong chaos, and observe destruction of Anderson localization in the packet tails for small values of $σ$.

preprint2010arXiv

The crossover from strong to weak chaos for nonlinear waves in disordered systems

We observe a crossover from strong to weak chaos in the spatiotemporal evolution of multiple site excitations within disordered chains with cubic nonlinearity. Recent studies have shown that Anderson localization is destroyed, and the wave packet spreading is characterized by an asymptotic divergence of the second moment $m_2$ in time (as $t^{1/3}$), due to weak chaos. In the present paper, we observe the existence of a qualitatively new dynamical regime of strong chaos, in which the second moment spreads even faster (as $t^{1/2}$), with a crossover to the asymptotic law of weak chaos at larger times. We analyze the pecularities of these spreading regimes and perform extensive numerical simulations over large times with ensemble averaging. A technique of local derivatives on logarithmic scales is developed in order to quantitatively visualize the slow crossover processes.

preprint2008arXiv

Global dynamics of coupled standard maps

Understanding the dynamics of multi--dimensional conservative dynamical systems (Hamiltonian flows or symplectic maps) is a fundamental issue of non-linear science. The Generalized ALignment Index (GALI), which was recently introduced and applied successfully for the distinction between regular and chaotic motion in Hamiltonian systems \cite{sk:6}, is an ideal tool for this purpose. In the present paper we make a first step towards the dynamical study of multi--dimensional maps, by obtaining some interesting results for a 4--dimensional (4D) symplectic map consisting of N=2 coupled standard maps \cite{Kan:1}. In particular, using the new GALI$_3$ and GALI$_4$ indices, we compute the percentages of regular and chaotic motion of the map equally reliably but much faster than previously used indices, like GALI$_2$ (known in the literature as SALI).

preprint2007arXiv

Application of the Generalized Alignment Index (GALI) method to the dynamics of multi--dimensional symplectic maps

We study the phase space dynamics of multi--dimensional symplectic maps, using the method of the Generalized Alignment Index (GALI). In particular, we investigate the behavior of the GALI for a system of N=3 coupled standard maps and show that it provides an efficient criterion for rapidly distinguishing between regular and chaotic motion.

preprint2007arXiv

Studying the global dynamics of conservative dynamical systems using the SALI chaos detection method

We use the Smaller ALignment Index (SALI) method of chaos detection, to study the global dynamics of conservative dynamical systems described by differential or difference equations. In particular, we consider the well--known 2 and 4--dimensional symplectic standard map, as well as an autonomous Hamiltonian system of 2 and 3 degrees of freedom describing the motion of a star in models of barred galaxies. The application of SALI helped us to compute rapidly and accurately the percentage of regular and chaotic motion for particular values of the parameters of these systems. We were also able to perform a computationally efficient determination of the dependence of these percentages on the variation of several parameters of the studied models.

preprint2003arXiv

How does the Smaller Alignment Index (SALI) distinguish order from chaos?

The ability of the Smaller Alignment Index (SALI) to distinguish chaotic from ordered motion, has been demonstrated recently in several publications.\cite{Sk01,GRACM} Basically it is observed that in chaotic regions the SALI goes to zero very rapidly, while it fluctuates around a nonzero value in ordered regions. In this paper, we make a first step forward explaining these results by studying in detail the evolution of small deviations from regular orbits lying on the invariant tori of an {\bf integrable} 2D Hamiltonian system. We show that, in general, any two initial deviation vectors will eventually fall on the ``tangent space'' of the torus, pointing in different directions due to the different dynamics of the 2 integrals of motion, which means that the SALI (or the smaller angle between these vectors) will oscillate away from zero for all time.