Local And Global Diffusion Along Resonant Lines in Discrete Quasi-integrable Dynamical Systems

We detect and measure diffusion along resonances in a quasi-integrable symplectic map for different values of the perturbation parameter. As in a previously studied Hamiltonian case (Lega et al., 2003) results agree with the prediction of the Nekhoroshev theorem. Moreover, for values of the perturbation parameter slightly below the critical value of the transition between Nekhoroshev and Chirikov regime we have also found a diffusion of some orbits along macroscopic portions of the phase space. Such a diffusion follows in a spectacular way the peculiar structure of resonant lines.     

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.     

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.     

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.     

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.     

Long-Term Stability Analysis of Quasi Integrable Degenerate Systems Through the Spectral Formulation of the Nekhoroshev Theorem

We describe a numerical application of the Nekhoroshev theorem to investigate the long-term stability of quasi-integrable systems. We extend the results of a previous paper to a class of degenerate systems, which are typical in celestial mechanics.     

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.     

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.     

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.     

A proof of Nekhoroshev's theorem for the stability times in nearly integrable Hamiltonian systems

In the present paper we give a proof of Nekhoroshev's theorem, which is concerned with an exponential estimate for the stability times in nearly integrable Hamiltonian systems. At variance with the already published proof, which refers to the case of an unperturbed Hamiltonian having the generic property of steepness, we consider here the particular case of a convex unperturbed Hamiltonian. The corresponding simplification in the proof might be convenient for an introduction to the subject.     

Directional diffusion coefficients of solar protons inside and outside the bow shock

The directional diffusion coefficients of low-energy (? 0.3 MeV) solar protons inside and outside the bow shock are examined during the solar flare event of 24 January 1969. The data are derived from simultaneous observations obtained by Explorer 33 inside the magnetosheath and by Explorer 35 in the interplanetary medium. Although the gross properties of the spin-averaged intensities on a diffusion-type plot appear to be the same in both media, the directional intensities show significant variations. It is shown that directional intensities of low-energy protons can be described reasonably well by anisotropic diffusion with an associated diffusion coefficient. Directional diffusion coefficients are found to differ by a factor as much as three among different directions in space, and from the spin-averaged diffusion coefficient. This suggests that anisotropic diffusion does indeed take place and that so called ‘isotropic’ diffusion coefficients derived in the past from spin-averaged intensities may actually be directional diffusion coefficients in cases where substantial anisotropies (> 50 per cent) exist. The typical postulated ratio of field aligned to cross-field diffusion coefficients is $\text{κ⊥}\text{κ∥}\text{< 0.1}$. The present data would indicate a ratio of ?0.3. This value of the anisotropy is to be taken only as an upper limit of the ratio because of the limitations introduced by the wide field of view of the detectors (~90°) and the lack of directional measurements over the entire sphere. Comparison between directional diffusion coefficients in the interplanetary medium and magnetosheath derived from identical directions in space implies changes in the parameters of the interplanetary magnetic field as it interacts with the bow shock.     

Nekhoroshev estimate for isochronous non resonant symplectic maps

We prove that non resonant isochronous symplectic maps in a neighborhood of an elliptic fixed point are stable for exponentially long times with the inverse of the distance from the fixed point. In the proof we make use of the majorant series method together with an idea for optimizing remainder estimates first applied to Hamiltonian problems by Nekhoroshev.     

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.     

Cosmic-ray diffusive acceleration at shock waves with finite upstream and downstream escape boundaries

In the present paper we discuss the modifications introduced into the first-order Fermi shock acceleration process due to a finite extent of diffusive regions near the shock or due to boundary conditions leading to an increased particle escape upstream and/or downstream of the shock. In the simple example of the planar shock wave considered we idealize the escape phenomenon by imposing a particle escape boundary at some distance from the shock. The presence of such a boundary (or boundaries) leads to coupled steepening of the accelerated particle spectrum and decreasing of the acceleration time scale. It allows for a semi-quantitative evaluation and, in some specific cases, also for modelling of the observed steep particle spectra as a result of the first-order Fermi shock acceleration. We also note that the particles close to the upper energy cut-off are younger than the estimate based on the respective acceleration time scale. In Appendix A we present a new time-dependent solution for infinite diffusive regions near the shock allowing for different constant diffusion coefficients upstream and downstream of the shock.     

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.     

The existence of a Smale horseshoe in a planar circular restricted four-body problem

In this paper we study the existence of a Smale horseshoe in a planar circular restricted four-body problem. For this planar four-body system there exists a transversal homoclinic orbit, but the fixed point is a degenerate saddle, so that the standard Smale–Birkhoff homoclinic theorem cannot be directly applied. We therefore apply the Conley–Moser conditions to prove the existence of a Smale horseshoe. Specifically, we first use the transversal structure of stable and unstable manifolds to make a linear transformation and then introduce a nonlinear Poincaré map $P$ by considering the truncated flow near the degenerate saddle; based on this Poincaré map $P$ , we define an invertible map $f$ , which is a composite function; by carefully checking the satisfiability of the Conley–Moser conditions for $f$ we finally prove that $f$ is a Smale horseshoe map, which implies that our restricted four-body problem has the chaotic dynamics of the Smale horseshoe type.     

Test Particle Simulations of Interaction Between Monochromatic Chorus Waves and Radiation Belt Relativistic Electrons

Chorus waves have been suggested to be effective in acceleration of radiation belt electrons. Here we perform gyro-averaged test-particle simulations to calculate the bounce-averaged pitch angle and energy diffusion coefficients for parallel-propagating monochromatic chorus waves, and perform a comparison of test-particle (TP) model with quasi-linear (QL) theory to evaluate the influence of nonlinear processes. For small amplitude chorus waves, the diffusion coefficients of TP and QL models are in good agreement. As the wave amplitude reaches a threshold value, two nonlinear processes (phase trapping and phase bunching) start to occur, especially at large equatorial pitch angles. Phase trapping yields rapid increases in pitch angle and kinetic energy. In contrast, phase bunching causes overall decreases in pitch angle and kinetic energy. For the waves with amplitudes slightly above the threshold value, the average behavior is dominated by the phase trapping, and TP diffusion coefficients are larger than QL ones. As wave amplitude increases, TP diffusion coefficients become smaller than QL ones, indicating that phase trapping gradually reduces the dominance over phase bunching.     

Steepness of Slopes at the <Emphasis Type="Italic">Luna-Glob</Emphasis> Landing Sites: Estimating by the Shaded Area Percentage in the LROC NAC Images

The paper presents estimates of the occurrence probability of slopes, whose steep surfaces could be dangerous for the landing of the Luna-Glob descent probe (Luna-25) given the baseline of the span between the landing pads (~3.5 m), for five potential landing ellipses. As a rule, digital terrain models built from stereo pairs of high-resolution images (here, the images taken by the Narrow Angle Camera onboard the Lunar Reconnaissance Orbiter (LROC NAC)) are used in such cases. However, the planned landing sites are at high latitudes (67°–74° S), which makes it impossible to build digital terrain models, since the difference in the observation angle of the overlapping images is insufficient at these latitudes. Because of this, to estimate the steepness of slopes, we considered the interrelation between the shaded area percentage in the image and the Sun angle over horizon at the moment of imaging. For five proposed landing ellipses, the LROC NAC images (175 images in total) with a resolution from 0.4 to 1.2 m/pixel were analyzed. From the results of the measurements in each of the ellipses, the dependence of the shaded area percentage on the solar angle were built, which was converted to the occurrence probability of slopes. For this, the data on the Apollo 16 landing region ware used, which is covered by both the LROC NAC images and the digital terrain model with high resolution. As a result, the occurrence probability of slopes with different steepness has been estimated on the baseline of 3.5 m for five landing ellipses according to the steepness categories of <7°, 7°–10°, 10°–15°, 15°–20°, and >20°.     

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.