This Is How You Lose Me

I liked the intimate setting of the class at first. The silence before the professor walked in. The cramped room. It always smelled like citrus cleaning products. Some hair gel mixed in there, too. There was peanut butter stuck on the roof of my mouth — from my sandwich at lunch — when he walked in that day, throwing a stack of Junot Díaz’s short story, “Alma,” onto the center of the shared table. I liked Junot Díaz’s writing. Loved it, actually. The way he captures pain and molds stories by weaving together the language of diary entries and love letters and suicide notes. The way he uses magic and grit and blood. Darkness and love simultaneously. His work is tragically beautiful. At least, that’s what I think. [excerpt

Testing a Fast Dynamical Indicator: The MEGNO

To investigate non-linear dynamical systems, like for instance artificial satellites, Solar System, exoplanets or galactic models, it is necessary to have at hand several tools, such as a reliable dynamical indicator. The aim of the present work is to test a relatively new fast indicator, the Mean Exponential Growth factor of Nearby Orbits (MEGNO), since it is becoming a widespread technique for the study of Hamiltonian systems, particularly in the field of dynamical astronomy and astrodynamics, as well as molecular dynamics. In order to perform this test we make a detailed numerical and statistical study of a sample of orbits in a triaxial galactic system, whose dynamics was investigated by means of the computation of the Finite Time Lyapunov Characteristic Numbers (FT-LCNs) by other authors.Comment: 25 pages, 35 figure

Chaotic diffusion of orbits in systems with divided phase space

In this paper we discuss the relevance of diffusive processes in multidimensional Hamiltonian systems. By means of a rather simple model, we present evidence that for moderate-to-strong chaotic systems the stochastic motion remains confined to disjoint domains on the energy surface, at least for mild motion times. We show that only for extremely large timescales and for rather large perturbations, does the chaotic component appear almost fully connected through the relics of the resonance structure. The discussion whether diffusion over the energy surface could actually occur in asteroidal or galaxy dynamics is also included.Fil: Giordano, Claudia Marcela. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica La Plata; ArgentinaFil: Cincotta, Pablo Miguel. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas; Argentin

Global Dynamics in Galactic Triaxial Systems I

In this paper we present a theoretical analysis of the global dynamics in a triaxial galactic system using a 3D integrable Hamiltonian as a simple representation. We include a thorough discussion on the effect of adding a generic non--integrable perturbation to the global dynamics of the system. We adopt the triaxial Stackel Hamiltonian as the integrable model and compute its resonance structure in order to understand its global dynamics when a perturbation is introduced. Also do we take profit of this example in order to provide a theoretical discussion about diffussive processes taking place in phase space.Comment: Accepted A&

Stochastic approach to diffusion inside the chaotic layer of a resonance

We model chaotic diffusion, in a symplectic 4D map by using the result of a theorem that was developed for stochastically perturbed integrable Hamiltonian systems. We explicitly consider a map defined by a free rotator (FR) coupled to a standard map (SM). We focus in the diffusion process in the action, $I$, of the FR, obtaining a semi--numerical method to compute the diffusion coefficient. We study two cases corresponding to a thick and a thin chaotic layer in the SM phase space and we discuss a related conjecture stated in the past. In the first case the numerically computed probability density function for the action $I$ is well interpolated by the solution of a Fokker-Planck (F-P) equation, whereas it presents a non--constant time delay respect to the concomitant F-P solution in the second case suggesting the presence of an anomalous diffusion time scale. The explicit calculation of a diffusion coefficient for a 4D symplectic map can be useful to understand the slow diffusion observed in Celestial Mechanics and Accelerator Physics.Comment: This is the author's version of a work that was submitted to Physical Review E (http://pre.aps.org

A comparison of different indicators of chaos based on the deviation vectors. Application to symplectic mappings

The aim of this research work is to compare the reliability of several variational indicators of chaos on mappings. The Lyapunov Indicator (LI); the Mean Exponential Growth factor of Nearby Orbits (MEGNO); the Smaller Alignment Index (SALI); the Fast Lyapunov Indicator (FLI); the Dynamical Spectra of stretching numbers (SSN) and the corresponding Spectral Distance (D); and the Relative Lyapunov Indicator (RLI), which is based on the evolution of the difference between two close orbits, have been included. The experiments presented herein allow us to reliably suggest a group of chaos indicators to analyze a general mapping. We show that a package composed of the FLI and the RLI (to analyze the phase portrait globally) and the MEGNO and the SALI (to analyze orbits individually) is good enough to make a description of the systems' dynamics.Comment: 25 pages, 40 figures. Celestial Mechanics and Dynamical Astronomy, in pres

The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects

We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, we consider the relative Poincare variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator MEGNO (a measure of the maximal Lyapunov exponent) that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best fit solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best fit solutions are found in dynamically active region of the phase space. The long term stability of the system is determined by a net of low-order two-body and three-body mean motion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myrs of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behavior.Comment: 11 pages (total), 8 figures. Accepted for publication in MNRAS. The definitive version will be/is available at http://www.blackwellpublishing.com. The astro-ph version is prepared with low resolution figures. To obtain the manuscript with full-resolution figures, please visit http://www.astri.uni.torun.pl/~chris/mnrasIII.ps.g

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