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.

We compare the divergence of orbits and the reversibility error for discrete time dynamical systems. These two quantities are used to explore the behavior of the global error induced by round off in the computation of orbits. The similarity of results found for any system we have analysed suggests the use of the reversibility error, whose computation is straightforward since it does not require the knowledge of the exact orbit, as a dynamical indicator.The statistics of fluctuations induced by round off for an ensemble of initial conditions has been compared with the results obtained in the case of random perturbations. Significant differences are observed in the case of regular orbits due to the correlations of round off error, whereas the results obtained for the chaotic case are nearly the same.Both the reversibility error and the orbit divergence computed for the same number of iterations on the whole phase space provide an insight on the local dynamical properties with a detail comparable with other dynamical indicators based on variational methods such as the finite time maximum Lyapunov characteristic exponent, the mean exponential growth factor of nearby orbits and the smaller alignment index. For 2D symplectic maps the differentiation between regular and chaotic regions is well full-filled. For 4D symplectic maps the structure of the resonance web as well as the nearby weakly chaotic regions are accurately described.

It has been recently demonstrated that both, a classical Schwarzschild black hole (BH), and a dense concentration of self-gravitating fermionic dark matter (DM) placed at the Galaxy centre, can explain the precise astrometric data (positions and radial velocities) of the S-stars orbiting Sgr A*. This result encompasses the 17 best resolved S-stars, and includes the test of general relativistic effects such as the gravitational redshift in the S2-star. In addition, the DM model features another remarkable result: the dense core of fermions is the central region of a continuous density distribution of DM whose diluted halo explains the Galactic rotation curve. In this Letter, we complement the above findings by analyzing in both models the relativistic periapsis precession of the S2-star orbit. While the Schwarzschild BH scenario predicts a unique prograde precession for S2, in the DM scenario it can be either retrograde or prograde, depending on the amount of DM mass enclosed within the S2 orbit, which in turn is a function of the DM fermion mass. We show that all the current and publicly available data of S2 can not discriminate between the two models, but upcoming S2 astrometry close to next apocentre passage could potentially establish if Sgr A* is governed by a classical BH or by a quantum DM system.