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&

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

Simple tools to study global dynamics in non-axisymmetric galactic potentials - I

In a first part we discuss the well-known problem of the motion of a star in a general non-axisymmetric 2D galactic potential by means of a very simple but almost universal system: the pendulum model. It is shown that both loop and box families of orbits arise as a natural consequence of the dynamics of the pendulum. An approximate invariant of motion is derived. A critical value of the latter sharply separates the domains of loops and boxes and a very simple computation allows to get a clear picture of the distribution of orbits on a given energy surface. Besides, a geometrical representation of the global phase space using the natural surface of section for the problem, the 2D sphere, is presented. This provides a better visualization of the dynamics. In a second part we introduce a new indicator of the basic dynamics, the Mean Exponential Growth factor of Nearby Orbits (MEGNO), that is suitable to investigate the phase space structure associated to a general Hamiltonian. When applied to the 2D logarithmic potential it is shown to be effective to obtain a picture of the global dynamics and, also, to derive good estimates of the largest Lyapunov characteristic number in realistic physical times. Comparisons with other techniques reveal that the MEGNO provides more information about the dynamics in the phase space than other wide used tools. Finally, we discuss the structure of the phase space associated to the 2D logarithmic potential for several values of the semiaxis ratio and energy. We focus our attention on the stability analysis of the principal periodic orbits and on the chaotic component. We obtain critical energy values for which connections between the main stochastic zones take place. In any case, the whole chaotic domain appears to be always confined to narrow filaments, with a Lyapunov time about three characteristic periods.Facultad de Ciencias Astronómicas y Geofísica

A randomized, double-blind, placebo-controlled trial to assess safety and tolerability during treatment of type 2 diabetes with usual diabetes therapy and either Cycloset™ or placebo

Background: Cycloset™ is a quick-release formulation of bromocriptine mesylate, a dopamine agonist, which in animal models of insulin resistance and type 2 diabetes acts centrally to reduce resistance to insulin- mediated suppression of hepatic glucose output and tissue glucose disposal. In such animals, bromocriptine also reduces hepatic triglyceride synthesis and free fatty acid mobilization, manifesting decreases in both plasma triglycerides and free fatty acids. In clinical trials, morning administration of Cycloset™ either as monotherapy or adjunctive therapy to sulfonylurea or insulin reduces HbA1c levels relative to placebo by 0.55–1.2. Cycloset™ therapy also reduces plasma triglycerides and free fatty acid by approximately 25% and 20%, respectively, among those also receiving sulfonylurea therapies. The effects of once-daily morning Cycloset™ therapy on glycemic control and plasma lipids are demonstrable throughout the diurnal portion of the day (7 a.m. to 7 p.m.) across postprandial time points. Methods/Design: 3,095 individuals were randomized in a 2:1 ratio into a one year trial aimed to assess the safety and efficacy of Cycloset™ compared to placebo among individuals receiving a variety of treatments for type 2 diabetes. Eligibility criteria for this randomized placebo controlled trial included: age 30–80, HbA1c ≤ 10%, diabetes therapeutic regimen consisting of diet or no more than two hypoglycemic agents or insulin with or without one additional oral agent (usual diabetes therapy; UDT). The primary safety endpoint will test the hypothesis that the rate of all-cause serious adverse events after one year of usual diabetes therapy (UDT) plus Cycloset™ is not greater than that for UDT plus placebo by more than an acceptable margin defined as a hazard ratio of 1.5 with a secondary endpoint analysis of the difference in the rate of serious cardiovascular events, (myocardial infarction, stroke, coronary revascularization or hospitalization for or angina or congestive heart failure). Efficacy analyses will evaluate effects of Cycloset™ versus placebo on change from baseline in HbA1c, fasting glucose, body weight, waist circumference, blood pressure and plasma lipids. Discussion: This study will extend the current data on Cycloset™ safety, tolerability and efficacy in individuals with type 2 diabetes to include its effects in combination with thiazolodinediones, insulin secretagogues, metformin, alpha-glucosidase inhibitors and exogenous insulin regimens. Trial registration: clinical trials.gov NCT0037767

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

Chirikov and Nekhoroshev diffusion estimates : Bridging the two sides of the river

We present theoretical and numerical results pointing towards a strong connection between the estimates for the diffusion rate along simple resonances in multidimensional nonlinear Hamiltonian systems that can be obtained using the heuristic theory of Chirikov and a more formal one due to Nekhoroshev. We show that, despite a wide-spread impression, the two theories are complementary rather than antagonist. Indeed, although Chirikov’s 1979 review has thousands of citations, almost all of them refer to topics such as the resonance overlap criterion, fast diffusion, the Standard or Whisker Map, and not to the constructive theory providing a formula to measure diffusion along a single resonance. However, as will be demonstrated explicitly below, Chirikov’s formula provides values of the diffusion coefficient which are quite well comparable to the numerically computed ones, provided that it is implemented on the so-called optimal normal form derived as in the analytic part of Nekhoroshev’s theorem. On the other hand, Chirikov’s formula yields unrealistic values of the diffusion coefficient, in particular for very small values of the perturbation, when used in the original Hamiltonian instead of the optimal normal form. In the present paper, we take advantage of this complementarity in order to obtain accurate theoretical predictions for the local value of the diffusion coefficient along a resonance in a specific 3DoF nearly integrable Hamiltonian system. Besides, we compute numerically the diffusion coefficient and a full comparison of all estimates is made for ten values of the perturbation parameter, showing a very satisfactory agreement.Instituto de Astrofísica de La Plat

Chirikov and Nekhoroshev diffusion estimates : Bridging the two sides of the river

We present theoretical and numerical results pointing towards a strong connection between the estimates for the diffusion rate along simple resonances in multidimensional nonlinear Hamiltonian systems that can be obtained using the heuristic theory of Chirikov and a more formal one due to Nekhoroshev. We show that, despite a wide-spread impression, the two theories are complementary rather than antagonist. Indeed, although Chirikov’s 1979 review has thousands of citations, almost all of them refer to topics such as the resonance overlap criterion, fast diffusion, the Standard or Whisker Map, and not to the constructive theory providing a formula to measure diffusion along a single resonance. However, as will be demonstrated explicitly below, Chirikov’s formula provides values of the diffusion coefficient which are quite well comparable to the numerically computed ones, provided that it is implemented on the so-called optimal normal form derived as in the analytic part of Nekhoroshev’s theorem. On the other hand, Chirikov’s formula yields unrealistic values of the diffusion coefficient, in particular for very small values of the perturbation, when used in the original Hamiltonian instead of the optimal normal form. In the present paper, we take advantage of this complementarity in order to obtain accurate theoretical predictions for the local value of the diffusion coefficient along a resonance in a specific 3DoF nearly integrable Hamiltonian system. Besides, we compute numerically the diffusion coefficient and a full comparison of all estimates is made for ten values of the perturbation parameter, showing a very satisfactory agreement.Instituto de Astrofísica de La Plat

Phase space structures and ionization dynamics of hydrogen atom in elliptically polarized microwaves

The multiphoton ionization of hydrogen atoms in a strong elliptically polarized microwave field exhibits complex features that are not observed for ionization in circular and linear polarized fields. Experimental data reveal high sensitivity of ionization dynamics to the small changes of the field polarization. The multidimensional nature of the problem makes widely used diagnostics of dynamics, such as Poincar\'{e} surfaces of section, impractical. We analyze the phase space dynamics using finite time stability analysis rendered by the fast Lyapunov Indicators technique. The concept of zero--velocity surface is used to initialize the calculations and visualize the dynamics. Our analysis provides stability maps calculated for the initial energy at the maximum and below the saddle of the zero-velocity surface. We estimate qualitatively the dependence of ionization thresholds on the parameters of the applied field, such as polarization and scaled amplitude

Application of the MEGNO technique to the dynamics of Jovian irregular satellites

We apply the MEGNO (Mean Exponential Growth of Nearby Orbits) technique to the dynamics of Jovian irregular satellites. We demonstrate the efficiency of applying the MEGNO indicator to generate a mapping of relevant phase-space regions occupied by observed jovian irregular satellites. The construction of MEGNO maps of the Jovian phase-space region within its Hill-sphere is addressed and the obtained results are compared with previous studies regarding the dynamical stability of irregular satellites. Since this is the first time the MEGNO technique is applied to study the dynamics of irregular satellites we provide a review of the MEGNO theory. We consider the elliptic restricted three-body problem in which Jupiter is orbited by a massless test satellite subject to solar gravitational perturbations. The equations of motion of the system are integrated numerically and the MEGNO indicator computed from the systems variational equations. An unprecedented large set of initial conditions are studied to generate the MEGNO maps. The chaotic nature of initial conditions are demonstrated by studying a quasi-periodic orbit and a chaotic orbit. As a result we establish the existence of several high-order mean-motion resonances detected for retrograde orbits along with other interesting dynamical features. The computed MEGNO maps allows to qualitatively differentiate between chaotic and quasi-periodic regions of the irregular satellite phase-space given only a relatively short integration time. By comparing with previous published results we can establish a correlation between chaotic regions and corresponding regions of orbital instability.Comment: 15 pages, 13 figures, 2 tables, submitted to MNRA