In this article, we explore both numerically and analytically how the dynamical environments of mildly relativistic binaries evolve with increasing the general relativity factor $\gamma$ (the normalized inverse of the binary size measured in the units of the gravitational radius corresponding to the total mass of the system). Analytically, we reveal a phenomenon of the relativistic shifting of mean-motion resonances: on increasing $\gamma$, the resonances between the test particle and the central binary shift, due to the relativistic variation of the mean motions of the primary and secondary binaries and the relativistic advance of the tertiary's pericenter. To exhibit the circumbinary dynamics globally, we numerically integrate equations of the circumbinary motion of a test particle, and construct relevant scans of the maximum Lyapunov exponents and stability diagrams in the "pericentric distance -- eccentricity" plane of initial conditions. In these scans and diagrams, regular and chaotic domains are identified straightforwardly. Our analytical and numerical estimates of the shift size are in a good agreement. Prospects for identification of the revealed effect in astronomical observations are discussed.

We present a novel numerical method aimed to characterize global behaviour, in particular chaotic diffusion, in dynamical systems. It is based on an analysis of the Poincar\'e recurrence statistics on massive grids of initial data or values of parameters. We concentrate on Hamiltonian systems, featuring the method separately for the cases of bounded and non-bounded phase spaces. The embodiments of the method in each of the cases are specific. We compare the performances of the proposed Poincar\'e recurrence method (PRM) and the custom Lyapunov exponent (LE) methods and show that they expose the global dynamics almost identically. However, a major advantage of the new method over the known global numerical tools, such as LE, FLI, MEGNO, and FA, is that it allows one to construct, in some approximation, charts of local diffusion timescales. Moreover, it is algorithmically simple and straightforward to apply.

We apply a generalized Kepler map theory to describe the qualitative chaotic dynamics around cometary nuclei, based on accessible observational data for five comets whose nuclei are well-documented to resemble dumb-bells. The sizes of chaotic zones around the nuclei and the Lyapunov times of the motion inside these zones are estimated. In the case of Comet 1P/Halley, the circumnuclear chaotic zone seems to engulf an essential part of the Hill sphere, at least for orbits of moderate to high eccentricity.

We study the dynamical chaos and integrable motion in the planar circular restricted three-body problem and determine the fractal dimension of the spiral strange repeller set of non-escaping orbits at different values of mass ratio of binary bodies and of Jacobi integral of motion. We find that the spiral fractal structure of the Poincar\'e section leads to a spiral density distribution of particles remaining in the system. We also show that the initial exponential drop of survival probability with time is followed by the algebraic decay related to the universal algebraic statistics of Poincar\'e recurrences in generic symplectic maps.

We study the capture of galactic dark matter particles (DMP) in two-body and few-body systems with a symplectic map description. This approach allows modeling the scattering of $10^{16}$ DMPs after following the time evolution of the captured particle on about $10^9$ orbital periods of the binary system. We obtain the DMP density distribution inside such systems and determine the enhancement factor of their density in a center vicinity compared to its galactic value as a function of the mass ratio of the bodies and the ratio of the body velocity to the velocity of the galactic DMP wind. We find that the enhancement factor can be on the order of tens of thousands.

The main features of 1P/Halley chaotic dynamics can be described by a two dimensional symplectic map. Using Mel'nikov integral we semi-analytically determine such a map for 1P/Halley taking into account gravitational interactions from the Sun and the eight planets. We determine the Solar system kick function ie the energy transfer to 1P/Halley along one passage through the Solar system. Our procedure allows to compute for each planet its contribution to the Solar system kick function which appears to be the sum of the Keplerian potential of the planet and of a rotating circular gravitational dipole potential due to the Sun movement around Solar system barycenter. We test the robustness of the symplectic Halley map by directly integrating Newton's equations over $\sim 2.4\cdot 10^4$ yr around Y2K and by reconstructing the Solar system kick function. Our results show that the Halley map with fixed parameters gives a reliable description of comet dynamics on time scales of $10^4$ yr while on a larger scales the parameters of the map are slowly changing due to slow oscillations of orbital momentum.