A numerical study of the stabilization effect of steepness

We present the results of numerical experiments about the influence of steepness on the resonant structure, stability and diffusion in a 4-dimensional symplectic map. The map is designed so that by changing a parameter, we smoothly switch steepness on and off by the change of the so called steepness coefficients. In both cases we measure the diffusion coefficients of the actions within a resonance. According to Nekhoroshev theorem we find that, in the steep case, the diffusion coefficients are definitely smaller than in the non steep one, thus confirming the threshold effect of the steepness coefficients which comes from the proof of Nekhoroshev theorem.     

An application of the Nekhoroshev theorem to the restricted three-body problem

We studied the stability of the restricted circular three-body problem. We introduced a model Hamiltonian in action-angle Delaunay variables. which is nearly-integrable with the perturbing parameter representing the mass ratio of the primaries. We performed a normal form reduction to remove the perturbation in the initial Hamiltonian to higher orders in the perturbing parameter. Next we applied a result on the Nekhoroshev theorem proved by Pöschel [13] to obtain the confinement in phase space of the action variables (related to the elliptic elements of the minor body) for an exponentially long time. As a concrete application. we selected the Sun-Ceres-Jupiter case, obtaining (after the proper normal form reduction) a stability result for a time comparable to the age of the solar system (i.e., 4.9 · 10⁹ years) and for a mass ratio of the primaries less or equal than 10^–6.     

Measure of the exponential splitting of the homoclinic tangle in four-dimensional symplectic mappings

Using four-dimensional symplectic maps as a model problem, we numerically compute the unstable manifolds of the hyperbolic manifolds of the phase space related to the single resonances. We measure an exponential dependence of the size of the lobes of these manifolds through many orders of magnitude of the perturbing parameter. This is an indirect numerical verification of the exponential decay of the normal form, as predicted by the Nekhoroshev theorem. The variation of the size of the lobes turns out to be correlated to the diffusion coefficient.     

Construction of a Nekhoroshev like result for the asteroid belt dynamical system

We investigate the possibility of obtaining a Nekhoroshev like result for the dynamical system describing the motion of an asteroid in the main belt, From the mathematical point of view this is a new result since the problem is degenerate and we want to control also the motion of degenerate actions, We find that there are regions, such as the resonances of low order among the fast angles (mean motion resonances), where a Nekhoroshev like result cannot be proved a priori, Conversely, we are able to confine the motions in the mean motion resonances of logarithmically large order in the perturbation parameters, as well as in the non-resonant region, We discuss also the connection with the existence of invariant tori.     

Kolmogorov and Nekhoroshev theory for the problem of three bodies

We investigate the long time stability in Nekhoroshev’s sense for the Sun– Jupiter–Saturn problem in the framework of the problem of three bodies. Using computer algebra in order to perform huge perturbation expansions we show that the stability for a time comparable with the age of the universe is actually reached, but with some strong truncations on the perturbation expansion of the Hamiltonian at some stage. An improvement of such results is currently under investigation.     

Probing the Nekhoroshev Stability of Asteroids

We apply the spectral formulation of the Nekhoroshev theorem to investigate the long-term stability of real main belt asteroids. We find numerical indication that some asteroids are in the so-called Nekhoroshev stability regime, that is they are on chaotic orbits but their motion is stable over very long times. We have analyzed the motion of bodies in different regions of the belt, to assess the sensitivity of our method. We found that it allows us to clearly discriminate between different dynamical regimes, such as the one described by the Nekhoroshev stability, the one well described by the KAM theory, and the unstable chaotic regime in which diffusion in phase space can be detected over time spans much shorter than the age of the solar system.     

Quantitative perturbation theory by successive elimination of harmonics

We revisit some results of perturbation theories by a method of successive elimination of harmonics inspired by some ideas of Delaunay. On the one hand, we give a connection between the KAM theorem and the Nekhoroshev theorem. On the other hand, we support in a quantitative fashion a semi-numerical method of analysis of a perturbed system recently introduced by one of the authors.     

A numerical study of the size of the homoclinic tangle of hyperbolic tori and its correlation with Arnold diffusion in Hamiltonian systems

Using a three degrees of freedom quasi-integrable Hamiltonian as a model problem, we numerically compute the unstable manifolds of the hyperbolic manifolds of the phase space related to single resonances. We measure an exponential dependence of the splitting of these manifolds through many orders of magnitude of the perturbing parameter. This is an indirect numerical verification of the exponential decay of the normal form, as predicted by the Nekhoroshev theorem. We also detect different transitions in the topology of these manifolds related to the local rational approximations of the frequencies. The variation of the size of the homoclinic tangle as well as the topological transitions turn out to be correlated to the speed of Arnold diffusion.     

On the connection between the Nekhoroshev theorem and Arnold diffusion

The analytical techniques of the Nekhoroshev theorem are used to provide estimates on the coefficient of Arnold diffusion along a particular resonance in the Hamiltonian model of Froeschlé et al. (Science 289:2108–2110, 2000). A resonant normal form is constructed by a computer program and the size of its remainder ||R _opt || at the optimal order of normalization is calculated as a function of the small parameter ${\epsilon}$ . We find that the diffusion coefficient scales as ${D \propto ||R_{opt}||^3}$ , while the size of the optimal remainder scales as ${||R_{opt}|| \propto {\rm exp}(1/\epsilon^{0.21})}$ in the range ${10^{-4} \leq \epsilon \leq 10^{-2}}$ . A comparison is made with the numerical results of Lega et al. (Physica D 182:179–187, 2003) in the same model.     

Local and global diffusion in the Arnold web of a priori unstable systems

Using the numerical techniques developed by Froeschlé et al. (Science 289 (5487): 2108–2110, 2000) and by Lega et al. (Physica D 182: 179–187, 2003) we have studied diffusion and stochastic properties of an a priori unstable 4D symplectic map. We have found two different kinds of diffusion that coexist for values of the perturbation below the critical value for the Chirikov overlapping of resonances. A fast diffusion along some resonant lines that exist already in the unperturbed case and a slow diffusion occurring in regions of the phase space far from such resonances. The latter one, although the system does not satisfy the Nekhoroshev hypothesis, decreases faster than a power law and possibly exponentially. We compare the diffusion coefficient to the indicators of stochasticity like the Lyapunov exponent and filling factor showing their behavior for chaotic orbits in regions of the Arnold web where the secondary resonances appear, or where they overlap.     

On the relationship between Lyapunov times and macroscopic instability times

On the basis of the general theory of Hamiltonian systems, we consider the relationship between Lyapunov times and macroscopic diffusion times. We find out that there are two regimes: the Nekhoroshev regime and the resonant overlapping regime. In the first case the diffusion time is exponentially long with respect to Lyapunov times. In the second case, the relationship is polynomial although we do not find any theoretical reason for the existence of a universal power law. We show numerical evidences which confirm our theoretical considerations.     

The Role of Hyperbolic Invariant Sets in Stickiness Effects

With the standard map model, we study the stickiness effect of invariant tori, particularly the role of hyperbolic sets in this effect. The diffusion of orbits originated from the neighborhoods of hyperbolic points, periodic islands and torus is studied. We find that they possess similar diffusion rules, but the diffusion of orbits originated from the neighborhood of a torus is faster than that originated near a hyperbolic set. The numerical results show that an orbit in the neighborhood of a torus spends most of time around hyperbolic invariant sets. We also calculate the areas of islands with different periods. The decay of areas with the periods obeys a power law, and the absolute values of the exponents increase monotonously with the perturbation parameter. According to the results obtained, we conclude that the stickiness effect of tori is caused mainly by the hyperbolic invariant sets near the tori, and the diffusion speed becomes larger when orbits diffuse away from the torus.     

The nekhoroshev theorem and the asteroid belt dynamical system

The present paper reviews the Nekhoroshev theorem from the point of view of physicists and astronomers. We point out that Nekhoroshev result is strictly connected with the existence of a specific structure of the phase space, the existence of which can be checked with several numerical tools. This is true also for a degenerate system such as the one describing the motion of an asteroid in the so called main belt. The main difference is that in some parts of the belt, the Nekhoroshev result cannot apply a priori. Mean motion resonances of order smaller than the logarithm of the mass of Jupiter and first order secular resonances must be excluded. In the remaining parts, conversely, the Nekhoroshev theorem can be proved, provided someparameters, such as the masses, the eccentricities and the inclinations of the planets are small enough. At the light of this result, a massive campaign of numerical integrations of real and fictitious asteroids should allow to understand which is the real dynamical structure of the asteroid belt.     

Stickiness in three-dimensional volume preserving mappings

We investigate the orbital diffusion and the stickiness effects in the phase space of a 3-dimensional volume preserving mapping. We first briefly review the main results about the stickiness effects in 2-dimensional mappings. Then we extend this study to the 3-dimensional case, studying for the first time the behavior of orbits wandering in the 3-dimensional phase space and analyzing the role played by the hyperbolic invariant sets during the diffusion process. Our numerical results show that an orbit initially close to a set of invariant tori stays for very long times around the hyperbolic invariant sets near the tori. Orbits starting from the vicinity of invariant tori or from hyperbolic invariant sets have the same diffusion rule. These results indicate that the hyperbolic invariant sets play an essential role in the stickiness effects. The volume of phase space surrounded by an invariant torus and its variation with respect to the perturbation parameter influences the stickiness effects as well as the development of the hyperbolic invariant sets. Our calculations show that this volume decreases exponentially with the perturbation parameter and that it shrinks down with the period very fast.     

Particle diffusion in degenerate,ionized gases

A perturbation treatment is used to find the distribution functions from the Boltzmann transport equation for a two-component ionized, degenerate gas. General expressions are found for the diffusion coefficient and relative diffusion velocity. Tables are given for specific values of spin and mass ratio.     

Negligibility of small divisor effects in the normal form theory for nearly-integrable Hamiltonians with decaying non-autonomous perturbations

The paper deals with the problem of the existence of a normal form for a nearly-integrable real-analytic Hamiltonian with aperiodically time-dependent perturbation decaying (slowly) in time. In particular, in the case of an isochronous integrable part, the system can be cast in an exact normal form, regardless of the properties of the frequency vector. The general case is treated by a suitable adaptation of the finite order normalization techniques usually used for Nekhoroshev arguments. The key point is that the so called “geometric part” is not necessary in this case. As a consequence, no hypotheses on the integrable part are required, apart from analyticity. The work, based on two different perturbative approaches developed by Giorgilli et al., is a generalisation of the techniques used by the same authors to treat more specific aperiodically time-dependent problems.     

Preface

We revisit the relegation algorithm by Deprit et al. (Celest. Mech. Dyn. Astron. 79:157–182, 2001) in the light of the rigorous Nekhoroshev’s like theory. This relatively recent algorithm is nowadays widely used for implementing closed form analytic perturbation theories, as it generalises the classical Birkhoff normalisation algorithm. The algorithm, here briefly explained by means of Lie transformations, has been so far introduced and used in a formal way, i.e. without providing any rigorous convergence or asymptotic estimates. The overall aim of this paper is to find such quantitative estimates and to show how the results about stability over exponentially long times can be recovered in a simple and effective way, at least in the non-resonant case.     

Strong bursts and diffusion in avalanching systems

The connection between avalanche dynamics and space physics has been studied for several years. In that context we recently suggested an avalanche model which explains the phenomena of reconnection. In this work the model is generalized to include the influence of an extremely strong perturbation, reflecting the effect of plasma storms originating from the sun. In addition, we allow for diffusion processes and show that the behavior changes with the onset of diffusion processes, rendering it quasi-periodic, along with the supression of small-size avalanches.     

Constraints on opacities from complex asteroseismology of B-type pulsators: the β Cephei star θ Ophiuchi

We present results of a comprehensive asteroseismic modelling of the β Cephei variable θ Ophiuchi. We call these studies complex asteroseismology because our goal is to reproduce both pulsational frequencies and corresponding values of a complex, non-adiabatic parameter, f , defined by the radiative flux perturbation. To this end, we apply the method of simultaneous determination of the spherical harmonic degree, ℓ, of excited pulsational mode and the corresponding non-adiabatic f parameter from combined multicolour photometry and radial velocity data. Using both the OP and OPAL opacity data, we find a family of seismic models which reproduce the radial and dipole centroid mode frequencies, as well as the f parameter associated with the radial mode. Adding the non-adiabatic parameter to seismic modelling of the B-type main-sequence pulsators yields very strong constraints on stellar opacities. In particular, only with one source of opacities it is possible to agree the empirical values of f with their theoretical counterparts. Our results for θ Oph point substantially to preference for the OPAL data.