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.

The long-term dynamics of planetesimals in debris discs in models with parameters of binary star systems Kepler-16, Kepler-34 and Kepler-35 with planets is investigated. Our calculations have shown the formation of a stable coorbital with the planet ring is possible for Kepler-16 and Kepler-35 systems. In Kepler-34 system, significant eccentricities of the orbits of the binary and planets can to prevent the formation of such a structure. Detection circumbinary annular structures in observations of systems binary stars can be evidence of the existence of planets, retaining coorbital rings from dust and planetesimals.

It is shown that circumbinary planetary systems are subject to universal tidal decay (shrinkage of orbits), caused by the forced orbital eccentricity inherent to them. Circumbinary planets (CBP) are liberated from parent systems, when, owing to the shrinkage, they enter the circumbinary chaotic zone. On shorter timescales (less than the current age of the Universe), the effect may explain, at least partially, the observed lack of CBP of close-enough (with periods < 5 days) stellar binaries; on longer timescales (greater than the age of the Universe but well within stellar lifetimes), it may provide massive liberation of chemically evolved CBP. Observational signatures of the effect may comprise (1) a prevalence of large rocky planets (super-Earths) in the whole population of rogue planets (if this mechanism were the only source of rogue planets); (2) a mass-dependent paucity of CBP in systems of low-mass binaries: the lower the stellar mass, the greater the paucity.

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.

It is shown that several habitability conditions (in fact, at least seven such conditions) appear to be fulfilled automatically by circumbinary planets of main-sequence stars (CBP-MS), whereas on Earth these conditions are fulfilled only by chance. Therefore, it looks natural that most of the production of replicating biopolymers in the Galaxy is concentrated on particular classes of CBP-MS, and life on Earth is an outlier, in this sense. In this scenario, Lathe's mechanism for the tidal "chain reaction" abiogenesis on Earth is favored as generic for CBP-MS, due to photo-tidal synchronization inherent to them. Problems with this scenario are discussed in detail.

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.

In massive numerical experiments we show that a planet embedded in a planetesimal disk induces a characteristic multi-lane "planetosignature" representing a pattern of several stellar-centric rings. If the planet's mass is large enough, the multi-lane signature degenerates to a three-lane one: then it consists of three rings, one bright coorbital with the planet, and two dark gaps in the radial distribution of the particles. The gaps correspond to orbital resonances 2:1 and 1:2 with the planet. This theoretical prediction may explain recent ALMA observations of the disk of HL Tau.

A condition upon which sporadic bursts (intermittent behaviour) of the relative energy become possible is derived for the motion in the chaotic layer around the separatrix of non-linear resonance. This is a condition for the existence of a marginal resonance, i.e. a resonance located at the border of the layer. A separatrix map in Chirikov's form [Chirikov, B. V., Phys. Reports 52, 263 (1979)] is used to describe the motion. In order to provide a straightforward comparison with numeric integrations, the separatrix map is synchronized to the surface of the section farthest from the saddle point. The condition of intermittency is applied to clear out the nature of the phenomenon of bursts of the eccentricity of chaotic asteroidal trajectories in the 3/1 mean motion commensurability with Jupiter. On the basis of the condition, a new intermittent regime of resonant asteroidal motion is predicted and then identified in numeric simulations.

Conditions for the emergence of a statistical relationship between $T_r$, the chaotic transport (recurrence) time, and $T_L$, the local Lyapunov time (the inverse of the numerically measured largest Lyapunov characteristic exponent), are considered for the motion inside the chaotic layer around the separatrix of a nonlinear resonance. When numerical values of the Lyapunov exponents are measured on a time interval not greater than $T_r$, the relationship is shown to resemble the quadratic one. This tentatively explains numerical results presented in the literature.

Planet formation scenarios and the observed planetary dynamics in binaries pose a number of theoretical challenges, especially in what concerns circumbinary planetary systems. We explore the dynamical stirring of a planetesimal circumbinary disk in the epoch when the gas component disappears. For this purpose, following theoretical approaches by Heppenheimer (1978) and Moriwaki and Nakagawa (2004), we develop a secular theory for the dynamics of planetesimals in circumbinary disks. If the binary is eccentric and its components have unequal masses, a spiral density wave is generated, engulfing the disk on the secular timescale, which may exceed 10^7 yr, depending on the problem parameters. The spiral pattern is transient; thus, its observed presence may betray system's young age. We explore the pattern both analytically and in numerical experiments. The derived analytical spiral is a modified lituus; it matches the numerical density wave in the gas-free case perfectly. Using the SPH scheme, we explore the effect of residual gas on the wave propagation.

The extent of the continuous zone of chaotic orbits of a small-mass tertiary around a system of two gravitationally bound primaries (a double star, a double black hole, a binary asteroid, etc.) is estimated analytically, in function of the tertiary's orbital eccentricity. The separatrix map theory is used to demonstrate that the central continuous chaos zone emerges (above a threshold in the primaries mass ratio) due to overlapping of the orbital resonances corresponding to the integer ratios p:1 between the tertiary and the central binary periods. In this zone, the unlimited chaotic orbital diffusion of the tertiary takes place, up to its ejection from the system. The primaries mass ratio, above which such a chaotic zone is universally present at all initial eccentricities of the tertiary, is estimated. The diversity of the observed orbital configurations of biplanetary and circumbinary exosystems is shown to be in accord with the existence of the primaries mass parameter threshold.

The problem of stability of the triangular libration points in the planar circular restricted three-body problem is considered. A software package, intended for normalization of autonomous Hamiltonian systems by means of computer algebra, is designed so that normalization problems of high analytical complexity could be solved. It is used to obtain the Birkhoff normal form of the Hamiltonian in the given problem. The normalization is carried out up to the 6th order of expansion of the Hamiltonian in the coordinates and momenta. Analytical expressions for the coefficients of the normal form of the 6th order are derived. Though intermediary expressions occupy gigabytes of the computer memory, the obtained coefficients of the normal form are compact enough for presentation in typographic format. The analogue of the Deprit formula for the stability criterion is derived in the 6th order of normalization. The obtained floating-point numerical values for the normal form coefficients and the stability criterion confirm the results by Markeev (1969) and Coppola and Rand (1989), while the obtained analytical and exact numeric expressions confirm the results by Meyer and Schmidt (1986) and Schmidt (1989). The given computational problem is solved without constructing a specialized algebraic processor, i.e., the designed computer algebra package has a broad field of applicability.

On the basis of tidal despinning timescale arguments, Peale showed in 1977 that the majority of irregular satellites (with unknown rotation states) are expected to reside close to their initial (fast) rotation states. Here we investigate the problem of the current typical rotation states among all known satellites from a viewpoint of dynamical stability. We explore location of the known planetary satellites on the ($\omega_0$, $e$) stability diagram, where $\omega_0$ is an inertial parameter of a satellite and $e$ is its orbital eccentricity. We show that most of the satellites with unknown rotation states cannot rotate synchronously, because no stable synchronous 1:1 spin-orbit state exists for them. They rotate either much faster than synchronously (those tidally unevolved) or, what is much less probable, chaotically (tidally evolved objects or captured slow rotators).

The planet Kepler-16b is known to follow a circumbinary orbit around a system of two main-sequence stars. We construct stability diagrams in the "pericentric distance - eccentricity" plane, which show that Kepler-16b is in a hazardous vicinity to the chaos domain - just between the instability "teeth" in the space of orbital parameters. Kepler-16b survives, because it is close to the stable half-integer 11/2 orbital resonance with the central binary, safe inside a resonance cell bounded by the unstable 5/1 and 6/1 resonances. The neighboring resonance cells are vacant, because they are "purged" by Kepler-16b, due to overlap of first-order resonances with the planet. The newly discovered planets Kepler-34b and Kepler-35b are also safe inside resonance cells at the chaos border.

Possible rotation states of two satellites of Saturn, Prometheus (S16) and Pandora (S17), are studied by means of numerical experiments. The attitude stability of all possible modes of synchronous rotation and the motion close to these modes is analyzed by means of computation of the Lyapunov spectra of the motion. The stability analysis confirms that the rotation of Prometheus and Pandora might be chaotic, though the possibility of regular behaviour is not excluded. For the both satellites, the attitude instability zones form series of concentric belts enclosing the main synchronous resonance center in the phase space sections. A hypothesis is put forward that these belts might form "barriers" for capturing the satellites in synchronous rotation. The satellites in chaotic rotation can mimic ordinary regular synchronous behaviour: they preserve preferred orientation for long periods of time, the largest axis of satellite's figure being directed approximately towards Saturn.

An essential role in the asteroidal dynamics is played by the mean motion resonances. Two-body planet-asteroid resonances are widely known, due to the Kirkwood gaps. Besides, so-called three-body mean motion resonances exist, in which an asteroid and two planets participate. Identification of asteroids in three-body (namely, Jupiter-Saturn-asteroid) resonances was initially accomplished by D.Nesvorny and A.Morbidelli (1998), who, by means of visual analysis of the time behaviour of resonant arguments, found 255 asteroids to reside in such resonances. We develop specialized algorithms and software for massive automatic identification of asteroids in the three-body, as well as two-body, resonances of arbitrary order, by means of automatic analysis of the time behaviour of resonant arguments. In the computation of orbits, all essential perturbations are taken into account. We integrate the asteroidal orbits on the time interval of 100000 yr and identify main-belt asteroids in the three-body Jupiter-Saturn-asteroid resonances up to the 6th order inclusive, and in the two-body Jupiter-asteroid resonances up to the 9th order inclusive, in the set of ~250000 objects from the "Asteroids - Dynamic Site" (AstDyS) database. The percentages of resonant objects, including extrapolations for higher-order resonances, are determined. In particular, the observed fraction of pure-resonant asteroids (those exhibiting resonant libration on the whole interval of integration) in the three-body resonances up to the 6th order inclusive is approximately 0.9% of the whole set; and, using a higher-order extrapolation, the actual total fraction of pure-resonant asteroids in the three-body resonances of all orders is estimated as approximately 1.1% of the whole set.

We estimate the Lyapunov times (characteristic times of predictability of motion) in Quillen's models for the dynamics in the solar neighborhood. These models take into account perturbations due to the Galactic bar and spiral arms. For estimating the Lyapunov times, an approach based on the separatrix map theory is used. The Lyapunov times turn out to be typically of the order of 10 Galactic years. We show that only in a narrow range of possible values of the problem parameters the Galactic chaos is adiabatic; usually it is not slow. We also estimate the characteristic diffusion times in the chaotic domain. In a number of models, the diffusion times turn out to be small enough to permit migration of the Sun from the inner regions of the Milky Way to its current location. Moreover, due to the possibility of ballistic flights inside the chaotic layer, the chaotic mixing might be even far more effective and quicker than in the case of normal diffusion. This confirms the dynamical possibility of Minchev and Famaey's migration concept.

We consider statistics of the disruption and Lyapunov times in an hierarchical restricted three-body problem. We show that at the edge of disruption the orbital periods and the size of the orbit of the escaping body exhibit L\'evy flights. Due to them, the time decay of the survival probability is heavy-tailed with the power-law index equal to -2/3, while the relation between the Lyapunov and disruption times is quasilinear. Applicability of these results in an "hierarchical resonant scattering" setting for a three-body interaction is discussed.

We investigate the problem of the typical rotation states of the small planetary satellites from the viewpoint of the dynamical stability of their rotation. We show that the majority of the discovered satellites with unknown rotation periods cannot rotate synchronously, because no stable synchronous 1:1 spin-orbit state exists for them. They rotate either much faster than synchronously (those tidally unevolved) or, what is much less probable, chaotically (tidally evolved objects or captured slow rotators).

We investigate statistics of the decay process in the equal-mass three-body problem with randomized initial conditions. Contrary to earlier expectations of similarity with "radioactive decay", the lifetime distributions obtained in our numerical experiments turn out to be heavy-tailed, i.e. the tails are not exponential, but algebraic. The computed power-law index for the differential distribution is within the narrow range, approximately from -1.7 to -1.4, depending on the virial coefficient. Possible applications of our results to studies of the dynamics of triple stars known to be at the edge of disruption are considered.

It is known that the dependence of the emission-line luminosity of a typical cloud in the active galactic nuclei broad-line regions (BLRs) upon the incident flux of ionizing continuum can be nonlinear. We study how this nonlinearity can be taken into account in estimating the size of the BLR by means of the "reverberation" methods. We show that the BLR size estimates obtained by cross-correlation of emission-line and continuum light curves can be much (up to an order of magnitude) less than the values obtained by reverberation modelling. This is demonstrated by means of numerical cross-correlation and reverberation experiments with model continuum flares and emission-line transfer functions and by means of practical reverberation modelling of the observed optical spectral variability of NGC 4151. The time behaviour of NGC 4151 in the H_alpha and H_beta lines is modelled on the basis of the observational data by Kaspi et al. (1996, ApJ, 470, 336) and the theoretical BLR model by Shevchenko (1984, Sov. Astron. Lett., 10, 377; 1985, Sov. Astron. Lett., 11, 35). The values of the BLR parameters are estimated that allow to judge on the size and physical characteristics of the BLR. The small size of the BLR, as determined by the cross-correlation method from the data of Kaspi et al. (1996, ApJ, 470, 336), is shown to be an artifact of this method. So, the hypothesis that the BLR size varies in time is not necessitated by the observational data.

This paper has been withdrawn by the author, due to necessity of revision.