The production of Tsallis entropy in the limit of weak chaos and a new indicator of chaoticity

We study the connection between the appearance of a `metastable' behavior of weakly chaotic orbits, characterized by a constant rate of increase of the Tsallis q-entropy (Tsallis 1988), and the solutions of the variational equations of motion for the same orbits. We demonstrate that the variational equations yield transient solutions, lasting for long time intervals, during which the length of deviation vectors of nearby orbits grows in time almost as a power-law. The associated power exponent can be simply related to the entropic exponent for which the q-entropy exhibits a constant rate of increase. This analysis leads to the definition of a new sensitive indicator distinguishing regular from weakly chaotic orbits, that we call `Average Power Law Exponent' (APLE). We compare the APLE with other established indicators of the literature. In particular, we give examples of application of the APLE in a) a thin separatrix layer of the standard map, b) the stickiness region around an island of stability in the same map, and c) the web of resonances of a 4D symplectic map. In all these cases we identify weakly chaotic orbits exhibiting the `metastable' behavior associated with the Tsallis q-entropy.Comment: 19 pages, 12 figures, accepted for publication by Physica

Stickiness in Chaos

We distinguish two types of stickiness in systems of two degrees of freedom (a) stickiness around an island of stability and (b) stickiness in chaos, along the unstable asymptotic curves of unstable periodic orbits. We studied these effects in the standard map with a rather large nonlinearity K=5, and we emphasized the role of the asymptotic curves U, S from the central orbit O and the asymptotic curves U+U-S+S- from the simplest unstable orbit around the island O1. We calculated the escape times (initial stickiness times) for many initial points outside but close to the island O1. The lines that separate the regions of the fast from the slow escape time follow the shape of the asymptotic curves S+,S-. We explained this phenomenon by noting that lines close to S+ on its inner side (closer to O1) approach a point of the orbit 4/9, say P1, and then follow the oscillations of the asymptotic curve U+, and escape after a rather long time, while the curves outside S+ after their approach to P1 follow the shape of the asymptotic curves U- and escape fast into the chaotic sea. All these curves return near the original arcs of U+,U- and contribute to the overall stickiness close to U+,U-. The isodensity curves follow the shape of the curves U+,U- and the maxima of density are along U+,U-. For a rather long time the stickiness effects along U+,U- are very pronounced. However after much longer times (about 1000 iterations) the overall stickiness effects are reduced and the distribution of points in the chaotic sea outside the islands tends to be uniform.Comment: 28 pages, 12 figure

Asymptotic Orbits in Barred Spiral Galaxies

We study the formation of the spiral structure of barred spiral galaxies, using an $N$-body model. The evolution of this $N$-body model in the adiabatic approximation maintains a strong spiral pattern for more than 10 bar rotations. We find that this longevity of the spiral arms is mainly due to the phenomenon of stickiness of chaotic orbits close to the unstable asymptotic manifolds originated from the main unstable periodic orbits, both inside and outside corotation. The stickiness along the manifolds corresponding to different energy levels supports parts of the spiral structure. The loci of the disc velocity minima (where the particles spend most of their time, in the configuration space) reveal the density maxima and therefore the main morphological structures of the system. We study the relation of these loci with those of the apocentres and pericentres at different energy levels. The diffusion of the sticky chaotic orbits outwards is slow and depends on the initial conditions and the corresponding Jacobi constant.Comment: 17 pages, 24 figure

Invariant manifolds and the response of spiral arms in barred galaxies

The unstable invariant manifolds of the short-period family of periodic orbits around the unstable Lagrangian points $L_1$ and $L_2$ of a barred galaxy define loci in the configuration space which take the form of a trailing spiral pattern. In the present paper we investigate this association in the case of the self-consistent models of Kaufmann & Contopoulos (1996) which provide an approximation of real barred-spiral galaxies. We also examine the relation of `response' models of barred-spiral galaxies with the theory of the invariant manifolds. Our main results are the following: The invariant manifolds yield the correct form of the imposed spiral pattern provided that their calculation is done with the spiral potential term turned on. We provide a theoretical model explaining the form of the invariant manifolds that supports the spiral structure. The azimuthal displacement of the Lagrangian points with respect to the bar's major axis is a crucial parameter in this modeling. When this is taken into account, the manifolds necessarily develop in a spiral-like domain of the configuration space, delimited from below by the boundary of a banana-like non-permitted domain, and from above either by rotational KAM tori or by cantori forming a stickiness zone. We construct `spiral response' models on the basis of the theory of the invariant manifolds and examine the connection of the latter to the `response' models (Patsis 2006) used to fit real barred-spiral galaxies, explaining how are the manifolds related to a number of morphological features seen in such models.Comment: 16 Page

Adaptive Memetic Particle Swarm Optimization with Variable Local Search Pool Size

We propose an adaptive Memetic Particle Swarm Optimization algorithm where local search is selected from a pool of different algorithms. The choice of local search is based on a probabilistic strategy that uses a simple metric to score the efficiency of local search. Our study investigates whether the pool size affects the memetic algorithm’s performance, as well as the possible benefit of using the adaptive strategy against a baseline static one. For this purpose, we employed the memetic algorithms framework provided in the recent MEMPSODE optimization software, and tested the proposed algorithms on the Benchmarking Black Box Optimization (BBOB 2012) test suite. The obtained results lead to a series of useful conclusions

Application of new dynamical spectra of orbits in Hamiltonian systems

In the present article, we investigate the properties of motion in Hamiltonian systems of two and three degrees of freedom, using the distribution of the values of two new dynamical parameters. The distribution functions of the new parameters, define the S(g) and the S(w) dynamical spectra. The first spectrum definition, that is the S(g) spectrum, will be applied in a Hamiltonian system of two degrees of freedom (2D), while the S(w) dynamical spectrum will be deployed in a Hamiltonian system of three degrees of freedom (3D). Both Hamiltonian systems, describe a very interesting dynamical system which displays a large variety of resonant orbits, different chaotic components and also several sticky regions. We test and prove the efficiency and the reliability of these new dynamical spectra, in detecting tiny ordered domains embedded in the chaotic sea, corresponding to complicated resonant orbits of higher multiplicity. The results of our extensive numerical calculations, suggest that both dynamical spectra are fast and reliable discriminants between different types of orbits in Hamiltonian systems, while requiring very short computation time in order to provide solid and conclusive evidence regarding the nature of an orbit. Furthermore, we establish numerical criteria in order to quantify the results obtained from our new dynamical spectra. A comparison to other previously used dynamical indicators, reveals the leading role of the new spectra.Comment: Published in Nonlinear Dynamics (NODY) journal. arXiv admin note: text overlap with arXiv:1009.1993 by other author

Acoustics of early universe. I. Flat versus open universe models

A simple perturbation description unique for all signs of curvature, and based on the gauge-invariant formalisms is proposed to demonstrate that: (1) The density perturbations propagate in the flat radiation-dominated universe in exactly the same way as electromagnetic or gravitational waves propagate in the epoch of the matter domination. (2) In the open universe, sounds are dispersed by curvature. The space curvature defines the minimal frequency $\omega_{\rm c}$ below which the propagation of perturbations is forbidden. Gaussian acoustic fields are considered and the curvature imprint in the perturbations spectrum is discussed.Comment: The new version extended by 2 sections. Changes in notation. Some important comments adde

The properties of the local spiral arms from RAVE data: two-dimensional density wave approach

Using the RAVE survey, we recently brought to light a gradient in the mean galactocentric radial velocity of stars in the extended solar neighbourhood. This gradient likely originates from non-axisymmetric perturbations of the potential, among which a perturbation by spiral arms is a possible explanation. Here, we apply the traditional density wave theory and analytically model the radial component of the two-dimensional velocity field. Provided that the radial velocity gradient is caused by relatively long-lived spiral arms that can affect stars substantially above the plane, this analytic model provides new independent estimates for the parameters of the Milky Way spiral structure. Our analysis favours a two-armed perturbation with the Sun close to the inner ultra-harmonic 4:1 resonance, with a pattern speed \Omega_p=18.6^{+0.3}_{-0.2} km/s/kpc and a small amplitude A=0.55 \pm 0.02% of the background potential (14% of the background density). This model can serve as a basis for numerical simulations in three dimensions, additionally including a possible influence of the galactic bar and/or other non-axisymmetric modes.Comment: 9 pages, 4 figures, accepted for publication in MNRA

A family of models of partially relaxed stellar systems I. Dynamical properties

[Abridged] Recently we have found that a family of models of partially relaxed, anisotropic stellar systems, inspired earlier by studies of incomplete violent relaxation, exhibits some interesting thermodynamic properties. Here we present a systematic investigation of its dynamical characteristics, in order to establish the basis for a detailed comparison with simulations of collisionless collapse, planned for a separate paper. For a full comparison with the observations of elliptical galaxies, the models should be extended to allow for the presence a sizable dark halo and of significant rotation. In the spherical limit, the family is characterized by two dimensionless parameters, i.e. \Psi, measuring the depth of the galaxy potential, and \nu, defining the form of a third global quantity Q, which is argued to be approximately conserved during collisionless collapse (in addition to the total energy and the total number of stars).Comment: 12 pages, 15 figures, accepted for publication in Astronomy & Astrophysic