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.

Small bodies of the Solar system, like asteroids, trans-Neptunian objects, cometary nuclei, planetary satellites, with diameters smaller than one thousand kilometers usually have irregular shapes, often resembling dumb-bells, or contact binaries. The spinning of such a gravitating dumb-bell creates around it a zone of chaotic orbits. We determine its extent analytically and numerically. We find that the chaotic zone swells significantly if the rotation rate is decreased, in particular, the zone swells more than twice if the rotation rate is decreased ten times with respect to the "centrifugal breakup" threshold. We illustrate the properties of the chaotic orbital zones in examples of the global orbital dynamics about asteroid 243 Ida (which has a moon, Dactyl, orbiting near the edge of the chaotic zone) and asteroid 25143 Itokawa.

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.

We study the capture of galactic dark matter particles in the Solar System produced by rotation of Jupiter. It is shown that the capture cross section is much larger than the area of Jupiter orbit being inversely diverging at small particle energy. We show that the dynamics of captured particles is chaotic and is well described by a simple symplectic dark map. This dark map description allows to simulate the scattering and dynamics of $10^{14}$ dark matter particles during the life time of the Solar System and to determine dark matter density profile as a function of distance from the Sun. The mass of captured dark matter in the radius of Neptune orbit is estimated to be $2 \cdot 10^{15} g$. The radial density of captured dark matter is found to be approximately constant behind Jupiter orbit being similar to the density profile found in galaxies.