Using the canonical Hamilton-Jacobi approach we study the Lynden-Bell (1979) concept of bar formation based on the idea of orbital trapping parallel to the long or short axes of the oval potential distortion. The concept considered a single parameter -- a sign of the derivative of the precession rate over angular momentum, determining the orientation of the trapped orbits. We derived a perturbation Hamiltonian which includes two more parameters characterising the background disc and the perturbation, that just as important as the previously known. This allows us to link this old theory with the matrix approach in linear perturbation theory, the theory of weak bars, and explain some features of the nonlinear secular evolution observed in N-body simulations.

A new so-called `gravitational loss-cone instability' in stellar systems has recently been investigated theoretically in the framework of linear perturbation theory and proved to be potentially important in understanding the physical processes in centres of galaxies, star clusters, and the Oort comet cloud. Using N-body simulations of a toy model, we confirm previous findings for the harmonic dominating potential and go beyond the linear theory. Unlike the well-known instabilities, the new one shows no notable change in the spherical geometry of the cluster, but it significantly accelerates the speed of diffusion of particles in phase space leading to an early instability saturation.

A new so-called `gravitational loss-cone instability' in stellar systems has recently been investigated theoretically in the framework of linear perturbation theory and proved to be potentially important in understanding the physical processes in centres of galaxies, star clusters, and the Oort comet cloud. Using N-body simulations, we confirm previous findings and go beyond the linear theory. Unlike the well-known instabilities, the new one shows no notable change in spherical geometry of the cluster, but it significantly accelerates the speed of diffusion of particles in phase space leading to a repopulation of the loss cone and early instability saturation.

Interaction of spiral density waves with stars in the vicinity of the inner Lindblad resonance in galactic discs is investigated using the linear perturbation theory and the leading orders in the epicyclic and WKB approximations. In analogy with shear flows in hydrodynamics, we conjecture that a weak nonlinearity in a narrow resonance region modifies the standard (Landau-Lin) bypass rule of the singularity to the integration in the principal value sense. This indeed leads to the reflection of the spiral wave instead of absorption, but the detailed picture looks awkward: the intervals of the wave weakening alternate with the intervals of the wave growth, so that the net absorption is absent. Incidentally, we rectify the result concerning leading spiral waves obtained earlier for the standard bypass rule.

Spiral patterns in some disc galaxies have two arms in the centre, and three or more arms on the periphery. The same result is also obtained in numerical simulations of stellar and gaseous discs. We argue that such patterns may occur due to fast cooling of the gas, resulting in formation of giant molecular clouds. The timescale of this process is 50 Myr, the factor of 10 shorter than of ordinary secular instability. The giant molecular clouds give rise to multiarm spirals through the mechanism of swing amplification.

Gaseous rotating razor-thin discs are a testing ground for theories of spiral structure that try to explain appearance and diversity of disc galaxy patterns. These patterns are believed to arise spontaneously under the action of gravitational instability, but calculations of its characteristics in the gas are mostly obscured, presumably due to a difficult outer boundary condition. The paper suggests a new effective method for finding the spiral patterns based on an expansion of small amplitude perturbations over finite radial elements. The final matrix equation is extracted from the original hydrodynamical equations without the use of an approximate theory and has a form of the linear algebraic eigenvalue problem. The method is applied to an exactly solvable model with finite outer boundary and to a galactic disc model.

Stationary stellar systems with radially elongated orbits are subject to radial orbit instability -- an important phenomenon that structures galaxies. Antonov (1973) presented a formal proof of the instability for spherical systems in the limit of purely radial orbits. However, such spheres have highly inhomogeneous density distributions with singularity $\sim 1/r^2$, resulting in an inconsistency in the proof. The proof can be refined, if one considers an orbital distribution close to purely radial, but not entirely radial, which allows to avoid the central singularity. For this purpose we employ non-singular analogs of generalised polytropes elaborated recently in our work in order to derive and solve new integral equations adopted for calculation of unstable eigenmodes in systems with nearly radial orbits. In addition, we establish a link between our and Antonov's approaches and uncover the meaning of infinite entities in the purely radial case. Maximum growth rates tend to infinity as the system becomes more and more radially anisotropic. The instability takes place both for even and odd spherical harmonics, with all unstable modes developing rapidly, i.e. having eigenfrequencies comparable to or greater than typical orbital frequencies. This invalidates orbital approximation in the case of systems with all orbits very close to purely radial.

Aims. We aim to show how encounters with low-mass satellite galaxies may alter the bar formation in a Milky Way-like disc galaxy. Methods. We use high-resolution N-body simulations of a disc galaxy prone to mild bar instability. For realistic initial conditions of satellites, we take advantage of cosmological simulations of Milky Way-like dark matter haloes. Results. The satellites may have a significant impact on the time of bar formation. Some runs with satellites demonstrate a delay, while others show an advancement in bar formation compared to the isolated run, with such time differences reaching $\sim$ 1 Gyr. Meanwhile, the final bar configuration, including its very appearance and the bar characteristics such as the pattern speed and the exponential growth rate of its amplitude are independent of the number of encounters and their orbits. The contribution of satellites with masses below $10^9 M_{\odot}$ is insignificant, unless their pericentre distances are small. We suggest that the encounters act indirectly via inducing perturbations across the disc that evolve to delayed waves in the central part and interfere with an emerging seed bar. The predicted effect for the present-day host galaxy is expected to be even more significant at redshifts $z \gtrsim 0.5$.

Many barred galaxies, possibly including the Milky Way, have cusps in the centres. There is a widespread belief, however, that usual bar instability taking place in bulgeless galaxy models is impossible for the cuspy models, because of the presence of the inner Lindblad resonance for any pattern speed. At the same time there are numerical evidences that the bar instability can form a bar. We analyse this discrepancy, by accurate and diverse N-body simulations and using the calculation of normal modes. We show that bar formation in cuspy galaxies can be explained by taking into account the disc thickness. The exponential growth time is moderate for typical current disc masses (about 250 Myr), but considerably increases (factor 2 or more) upon substitution of the live halo and bulge with a rigid halo/bulge potential; meanwhile pattern speeds remain almost the same. Normal mode analysis with different disc mass favours a young bar hypothesis, according to which the bar instability saturated only recently.

Using 3D N-body simulations we analyse an onset of the bar in cuspy models, and argue that role of swing amplification is twofold. Amplified shot noise due to disc discreteness hampers bar formation, while induced resonance perturbations allow bar amplitude to overcome shots. A bar pattern speed and a growth rate obtained in N-body simulations agree well with global mode analysis.

We show by numerical simulations that a purely stellar dynamical model composed of an exponential disc, a cuspy bulge, and an NFW halo with parameters relevant to the Milky Way Galaxy is subject to bar formation. Taking into account the finite disc thickness, the bar formation can be explained by the usual bar instability, in spite of the presence of an inner Lindblad resonance, that is believed to damp any global modes. The effect of replacing the live halo and bulge by a fixed external axisymmetric potential (rigid models) is studied. It is shown that while the e-folding time of bar instability increases significantly (from 250 to 500 Myr), the bar pattern speed remains almost the same. For the latter, our average value of 55 km/s/kpc agrees with the assumption that the Hercules stream in the solar neighbourhood is an imprint of the bar--disc interaction at the outer Lindblad resonance of the bar. Vertical averaging of the radial force in the central disc region comparable to the characteristic scale length allows us to reproduce the bar pattern speed and the growth rate of the rigid models, using normal mode analysis of linear perturbation theory in a razor thin disc. The strong increase of the e-folding time with decreasing disc mass predicted by the mode analysis suggests that bars in galaxies similar to the Milky Way have formed only recently.

We consider a two-parametric family of radially anisotropic models with non-singular density distribution in the centre. If highly eccentric orbits are locked near the centre, the characteristic growth rate of the instability is much less than the Jeans and dynamic frequencies of the stars (slow modes). The instability occurs only for even spherical harmonics and the perturbations are purely growing (aperiodic). On the contrary, if all orbits nearly reach the outer radius of the sphere, both even and odd harmonics are unstable. Unstable odd modes oscillate having characteristic frequencies of the order of the dynamical frequencies (fast modes). Unstable even harmonics contain a single aperiodic mode and several oscillatory modes, the aperiodic mode being the most unstable. The question of the nature of the radial orbit instability (ROI) is revisited. Two main interpretations of ROI were suggested in the literature. The first one refers to the classical Jeans instability associated with the lack of velocity dispersion of stars in the transverse direction. The second one refers to Lynden-Bell's orbital approach to bar formation in disc galaxies, which implies slowness and bi-symmetry of the perturbation. Oscillatory modes, odd spherical harmonics modes, and non-slow modes found in one of the models show that the orbital interpretation is not the only possible.

We test methods for the determination of unstable modes in stellar discs: a point collocation scheme in the action sub-space, a scheme based on expansion of the density and potential on the biorthonormal basis, and a finite element method. Using models of galaxies with low and high mass concentration to the center, the existence of two different kinds of spectra of unstable modes is demonstrated. Characteristic features of methods and obtained spectra are discussed. Despite ignoring any constraint on the continuity or differentiability of the perturbed DF, the collocation scheme is reliable for obtaining spectra and patterns in both kinds of models. The method based on the expansion of the perturbed potential and surface density over a biorthogonal basis was not applicable to a model with high mass concentration. The finite element method successfully used in various fields of science and engineering is currently sensitive to the presence of resonant orbits due to the choice of interpolation functions for the distribution function.

We study a new class of equilibrium two-parametric distribution functions of spherical stellar systems with radially anisotropic velocity distribution of stars. The models are less singular counterparts of the so called generalized polytropes, widely used in works on equilibrium and stability of gravitating systems in the past. The offered models, unlike the generalized polytropes, have finite density and potential in the center. The absence of the singularity is necessary for proper consideration of the radial orbit instability, which is the most important instability in spherical stellar systems. Comparison of the main observed parameters (potential, density, anisotropy) predicted by the present models and other popular equilibrium models is provided.

Several mechanisms of bar formation in stellar galactic disks are considered, including Toomre swing amplification and normal mode approach. On example of the well-known model of Kuzmin--Toomre using N-body simulations it was shown that the stellar bar is developed as a result of the evolution of an unstable normal mode. The pattern speed and the growth rate found agree well with linear perturbation theory. Nonlinear evolution of the bar is followed. Role of the growing transient spirals in bar formation is discussed.

We discuss some contradictions found in the literature concerning the problem of stability of collisionless spherical stellar systems which are the simplest anisotropic generalization of the well-known polytrope models. Their distribution function $F(E,L)$ is a product of power-low functions of the energy $E$ and the angular momentum $L$, i.e. $F\propto L^{-s}(-E)^q$. On the one hand, calculation of the growth rates in the framework of linear stability theory and N-body simulations show that these systems become stable when the parameter $s$ characterizing the velocity anisotropy of the stellar distribution is lower than some finite threshold value, $s<s_\textrm{crit}$. On the other hand Palmer & Papaloizou (1987) showed that the instability remained up to the isotropic limit $s=0$. Using our method of determining the eigenmodes for stellar systems, we show that the growth rates in weakly radially-anisotropic systems are indeed positive, but decrease exponentially as the parameter $s$ approaches zero, i.e. $\gamma\propto \exp(-s_{\ast}/s)$. In fact, for the systems with finite lifetime this means stability.

Small perturbations in spherical and thin disk stellar clusters surrounding massive a black hole are studied. Due to the black hole, stars with sufficiently low angular momentum escape from the system through the loss cone. We show that stability properties of spherical clusters crucially depend on whether the distribution of stars is monotonic or non-monotonic in angular momentum. It turns out that only non-monotonic distributions can be unstable. At the same time the instability in disk clusters is possible for both types of distributions.

We study spherical and disk clusters in a near-Keplerian potential of galactic centers or massive black holes. In such a potential orbit precession is commonly retrograde, i.e. direction of the orbit precession is opposite to the orbital motion. It is assumed that stellar systems consist of nearly radial orbits. We show that if there is a loss cone at low angular momentum (e.g., due to consumption of stars by a black hole), an instability similar to loss-cone instability in plasma may occur. The gravitational loss-cone instability is expected to enhance black hole feeding rates. For spherical systems, the instability is possible for the number of spherical harmonics $l \ge 3$. If there is some amount of counter-rotating stars in flattened systems, they generally exhibit the instability independently of azimuthal number $m$. The results are compared with those obtained recently by Tremaine for distribution functions monotonically increasing with angular momentum. The analysis is based on simple characteristic equations describing small perturbations in a disk or a sphere of stellar orbits highly elongated in radius. These characteristic equations are derived from the linearized Vlasov equations (combining the collisionless Boltzmann kinetic equation and the Poisson equation), using the action-angle variables. We use two techniques for analyzing the characteristic equations: the first one is based on preliminary finding of neutral modes, and the second one employs a counterpart of the plasma Penrose-Nyquist criterion for disk and spherical gravitational systems.

This paper presents a new approach to studying galactic structures. They are considered as the low-frequency normal modes in a disc of orbits precessing at different angular speeds. Such a concept is an adequate alternative to the commonly used approach of treating the disc as a set of individual stars rotating at near-circular orbits around the centre. The problem of determining the normal modes is reduced to a simple integral equation in the form of the classical eigen-value problem, where the eigen-value is directly equal to the pattern speed of the mode, $\Omega_p$. An examination of the general properties of the basic integral equation shows that two types of solutions exist, bar-like and spiral. The numerical solutions of both types are obtained. The characteristic pattern speeds are of the order of the mean orbit precession speed, although for the bar-modes $\Omega_p$ can markedly exceed the maximum precessing speed of orbits. It is shown that the bar-mode grows due to the immediate action of its gravitational field on the stars at the resonance regions. As for the spiral mode, its excitation is probably due to the inner Lindblad resonance that can promote mode growth.

The arguments in favor of the unified formation mechanism for both slow (Lynden-Bell's) bars and common fast bars are given. This mechanism consists in a certain instability that is akin to the well-known radial orbit instability; it is caused by the mutual attraction and alignment of axes of precessing star orbits (up to now, such a way of formation was considered only for slow bars). The general theory of the low-frequency modes of a disk consisting of precessing orbits (at different angular velocities) is presented. The problem of determining these modes is reduced to the integral equations of a rather simple structure. The characteristic pattern speeds ($\Omega_p$) of the low-frequency modes are of order of the mean orbit precession speeds ($\bar\Omega_{pr}$). The bar-modes also belong to this type of modes. The slow bars have $\Omega_p \approx \bar\Omega_{pr}$; for the fast bars, $\Omega_p$ may far exceed even the maximum precessing speed of disk orbits (however, $\Omega_p$ remains to be of order of these precessing speeds). The possibility of such an excess of $\Omega_p$ over $\Omega_{pr}^{\max}$ is connected with the effect of ``repelling'' orbits that tend to move in the direction opposite to that they are being pushed. The preliminary analysis of the orbit precession patterns for a number of typical potentials is given. It is noted that the maximum radius of the ``attracting'' circular orbits ($r_c$) may be used as a reasonable estimate of a bar length.

We consider a scenario of formation of the spiral structure in barred galaxies. This scenario includes the new non-resonant mechanism of elongation of spirals, due to the characteristic behaviour of the gravitational potential beyond the principal spiral arms.

The velocity field of the nearly face-on galaxy NGC 3631, derived from observations in the H alpha line and the HI radio line, is analysed to study perturbations related to the spiral structure of the galaxy. We confirm our previous conclusion that the line-of-sight velocity field gives evidence of the wave nature of the observed two-armed spiral structure. Fourier analysis of the observed velocity field is used to determine the location of corotation of the spiral structure of this galaxy, and the radius of corotation R_c is found to be about 42 arcsec, or 3.2 kpc. The vector velocity field of the gas in the plane of the disc is restored, and taking into account that we previously investigated vertical motions, we now have a full 3D gaseous velocity field of the galaxy. We show clear evidence of the existence of two anticyclonic and four cyclonic vortices near corotation in a frame of reference rotating with the spiral pattern. The centres of the anticyclones lie between the observed spiral arms. The cyclones lie close to the observed spirals, but their centres are shifted from the maxima in brightness.

We report detection of giant cyclonic vortices in the gaseous disc of the spiral galaxy NGC 3631 in the reference frame rotating with the spiral pattern. A presence of such structures was predicted by the authors for galaxies, where the radial gradient of the perturbed velocity exceeds that of the rotational velocity. This situation really takes place in NGC 3631.