We use the Fokker-Planck equation and model the dispersive dynamics of solid particles in annular protoplanetary disks whose gas component is more massive than the particle phase. We model particle--gas interactions as hard sphere collisions, determine the functional form of diffusion coefficients, and show the existence of two global unstable modes in the particle phase. These modes have spiral patterns with the azimuthal wavenumber $m=1$ and rotate slowly. We show that in ring-shaped disks, the phase space density of solid particles increases linearly in time towards an accumulation point near the location of pressure maximum, while instabilities grow exponentially. Therefore, planetesimals and planetary cores can be efficiently produced near the peaks of unstable density waves. In this mechanism, particles migrating towards the accumulation point will not participate in the formation of planets, and should eventually form a debris ring like the main asteroid belt or classical Kuiper belt objects. We present the implications of global instabilities to the formation of ice giants and terrestrial planets in the solar system.

We study the linear perturbations of collisionless near-Keplerian discs. Such systems are models for debris discs around stars and the stellar discs surrounding supermassive black holes at the centres of galaxies. Using a finite-element method, we solve the linearized collisionless Boltzmann equation and Poisson's equation for a wide range of disc masses and rms orbital eccentricities to obtain the eigenfrequencies and shapes of normal modes. We find that these discs can support large-scale `slow' modes, in which the frequency is proportional to the disc mass. Slow modes are present for arbitrarily small disc mass so long as the self-gravity of the disc is the dominant source of apsidal precession. We find that slow modes are of two general types: parent modes and hybrid child modes, the latter arising from resonant interactions between parent modes and singular van Kampen modes. The most prominent slow modes have azimuthal wavenumbers $m=1$ and $m=2$. We illustrate how slow modes in debris discs are excited during a fly-by of a neighbouring star. Many of the non-axisymmetric features seen in debris discs (clumps, eccentricity, spiral waves) that are commonly attributed to planets could instead arise from slow modes; the two hypotheses can be distinguished by long-term measurements of the pattern speed of the features.

We describe a new finite element method (FEM) to construct continuous equilibrium distribution functions of stellar systems. The method is a generalization of Schwarzschild's orbit superposition method from the space of discrete functions to continuous ones. In contrast to Schwarzschild's method, FEM produces a continuous distribution function (DF) and satisfies the intra element continuity and Jeans equations. The method employs two finite-element meshes, one in configuration space and one in action space. The DF is represented by its values at the nodes of the action-space mesh and by interpolating functions inside the elements. The Galerkin projection of all equations that involve the DF leads to a linear system of equations, which can be solved for the nodal values of the DF using linear or quadratic programming, or other optimization methods. We illustrate the superior performance of FEM by constructing ergodic and anisotropic equilibrium DFs for spherical stellar systems (Hernquist models). We also show that explicitly constraining the DF by the Jeans equations leads to smoother and/or more accurate solutions with both Schwarzschild's method and FEM.

I formulate a general finite element method (FEM) for self-gravitating stellar systems. I split the configuration space to finite elements, and express the potential and density functions over each element in terms of their nodal values and suitable interpolating functions. General expressions are then introduced for the Hamiltonian and phase space distribution functions of the stars that visit a given element. Using the weighted residual form of Poisson's equation, I derive the Galerkin projection of the perturbed collisionless Boltzmann equation, and assemble the global evolutionary equations of nodal distribution functions. The FEM is highly adaptable to all kinds of potential and density profiles, and it can deal with density clumps and initially non-axisymmetric systems. I use ring elements of non-uniform widths, choose linear and quadratic interpolation functions in the radial direction, and apply the FEM to the stability analysis of the cutout Mestel disc. I also integrate the forced evolutionary equations and investigate the disturbances of a stable stellar disc due to the gravitational field of a distant satellite galaxy. The performance of the FEM and its prospects are discussed.

We use the weighted integral form of spherical Bessel functions, and introduce a new analytical set of complete and biorthogonal potential--density basis functions. The potential and density functions of the new set have finite central values and they fall off, respectively, similar to $r^{-(1+l)}$ and $r^{-(4+l)}$ at large radii where $l$ is the latitudinal quantum number of spherical harmonics. The lowest order term associated with $l=0$ is the perfect sphere of de Zeeuw. Our basis functions are intrinsically suitable for the modeling of three dimensional, soft-centred stellar systems and they complement the basis sets of Clutton-Brock, Hernquist & Ostriker and Zhao. We test the performance of our functions by expanding the density and potential profiles of some spherical and oblate galaxy models.

I decompose the unstable growing modes of stellar disks to their Fourier components and present the physical mechanism of instabilities in the context of resonances. When the equilibrium distribution function is a non-uniform function of the orbital angular momentum, the capture of stars into the corotation resonance imbalances the disk angular momentum and triggers growing bar and spiral modes. The stellar disk can then recover its angular momentum balance through the response of non-resonant stars. I carry out a complete analysis of orbital structure corresponding to each Fourier component in the radial angle, and present a mathematical condition for the occurrence of van Kampen modes, which constitute a continuous family. I discuss on the discreteness and allowable pattern speeds of unstable modes and argue that the mode growth is saturated due to the resonance overlapping mechanism. An individually growing mode can also be suppressed if the corotation and inner Lindblad resonances coexist and compete to capture a group of stars. Based on this mechanism, I show that self-consistent scale-free disks with a sufficient distribution of non-circular orbits should be stable under perturbations of angular wavenumber $m>1$. I also derive a criterion for the stability of stellar disks against non-axisymmetric excitations.

A set of bi-orthogonal potential-density basis functions is introduced to model the density and its associated gravitational field of three dimensional stellar systems. Radial components of our basis functions are weighted integral forms of spherical Bessel functions. We discuss the properties of our basis functions and demonstrate their shapes for the latitudinal Fourier number $l=2$.

We present N-body simulations of unstable spiral modes in a dynamically cool collisionless disc. We show that spiral modes grow in a thin collisionless disk in accordance with the analytical perturbation theory. We use the particle-mesh code SUPERBOX with nested grids to follow the evolution of unstable spirals that emerge from an unstable equilibrium state. We use a large number of particles (up to 40 million particles) and high-resolution spatial grids in our simulations (128^3 cells). These allow us to trace the dynamics of the unstable spiral modes until their wave amplitudes are saturated due to nonlinear effects. In general, the results of our simulations are in agreement with the analytical predictions. The growth rate and the pattern speed of the most unstable bar-mode measured in N-body simulations agree with the linear analysis. However the parameters of secondary unstable modes are in lesser agreement because of the still limited resolution of our simulations.

I utilize the Petrov-Galerkin formulation and develop a new method for solving the unsteady collisionless Boltzmann equation in both the linear and nonlinear regimes. In the first order approximation, the method reduces to a linear eigenvalue problem which is solved using standard numerical methods. I apply the method to the dynamics of a model stellar disk which is embedded in the field of a soft-centered logarithmic potential. The outcome is the full spectrum of eigenfrequencies and their conjugate normal modes for prescribed azimuthal wavenumbers. The results show that the fundamental bar mode is isolated in the frequency space while spiral modes belong to discrete families that bifurcate from the continuous family of van Kampen modes. The population of spiral modes in the bifurcating family increases by cooling the disk and declines by increasing the fraction of dark to luminous matter. It is shown that the variety of unstable modes is controlled by the shape of the dark matter density profile.

We present a perturbation theory for studying the instabilities of non-axisymmetric gaseous discs. We perturb the dynamical equations of self-gravitating fluids in the vicinity of a non-axisymmetric equilibrium, and expand the perturbed physical quantities in terms of a complete basis set and a small non-axisymmetry parameter $\epsilon$. We then derive a linear eigenvalue problem in matrix form, and determine the pattern speed, growth rate and mode shapes of the first three unstable modes. In non-axisymmetric discs the amplitude and phase angle of travelling waves are functions of both the radius $R$ and the azimuthal angle $\phi$. This is due to the interaction of different wave components in the response spectrum. We demonstrate that wave interaction in unstable discs with small initial asymmetries, can develop dense clumps during the phase of exponential growth. Local clumps, which occur on the major spiral arms, can constitute seeds of gas giant planets in accretion discs.

We study bisymmetric modes of angular wavenumber 2 for flat stellar disks in potentials with smooth cores. Stars either all circulate in the same direction or a small fraction may counter-rotate. The bisymmetric modes are unstable unless there is a sufficiently large external halo or bulge. We find two modes for each disk: a more central fundamental mode and a more extensive and more spiral (trailing) secondary mode. The fundamental mode is particularly sensitive to the orbital population in the central part of the disk. Depending on that population, it varies from a small compact bar to a trailing spiral that is almost as wound as the secondary mode. All modes rotate too rapidly for there to be an inner Lindblad resonance. They transfer angular momentum from the central to the outer regions of the disk. Most of them release gravitational energy and convert it to kinetic energy, which also flows outwards through the disk.

We introduce exact analytical solutions of the steady-state hydrodynamic equations of scale-free, self-gravitating gaseous disks with flat rotation curves. We express the velocity field in terms of a stream function and obtain a third-order ordinary differential equation (ODE) for the angular part of the stream function. We present the closed-form solutions of the obtained ODE and construct hydrodynamical counterparts of the power-law and elliptic disks, for which self-consistent stellar dynamical models are known. We show that the kinematics of the Large Magellanic Cloud can well be explained by our findings for scale-free elliptic disks.

We modify the curvature condition for the existence of self-consistent scale-free discs, introduced by Zhao, Carollo & de Zeeuw. We survey the parameter space of the power-law discs, and show that the modified curvature condition is in harmony with the results of Schwarzschild's numerical orbit superposition method. We study the orbital structure of the power-law discs, and find a correlation between the population of centrophobic banana orbits and the non-self-consistency index. We apply the curvature condition to other families of scale-free elongated discs and find that it rules out a large range of power-law slopes and axis ratios. We generalize the condition, and apply it, to three-dimensional scale-free axisymmetric galaxy models.

We construct phase-space distribution functions for the oblate, cuspy mass models of Sridhar & Touma, which may contain a central point mass (black hole) and have potentials of St\"ackel form in parabolic coordinates. The density in the ST models is proportional to a power $r^{-\gamma}$ of the radius, with $0<\gamma<1$. We derive distribution functions $f(E, L_z)$ for the scale-free ST models (no black hole) using a power series of the energy $E$ and the component $L_z$ of the angular momentum parallel to the symmetry axis. We use the contour integral method of Hunter & Qian to construct $f(E, L_z)$ for ST models with central black holes, and employ the scheme introduced by Dejonghe & de Zeeuw to derive more general distribution functions which depend on $E$, $L_z$ and the exact third integral $I_3$. We find that self-consistent two- and three-integral distribution functions exist for all values $0 < \gamma < 1$.

We introduce a new class of 2-D mass models, whose potentials are of St\"ackel form in elliptic coordinates. Our model galaxies have two separate strong cusps that form double nuclei. The potential and surface density distributions are locally axisymmetric near the nuclei and become {\it highly} non-axisymmetric outside the nucleus. The surface density diverges toward the cuspy nuclei with the law $\Sigma \propto r^{-2}$. Our model is sustained by four general types of regular orbits: {\it butterfly}, {\it nucleuphilic banana}, {\it horseshoe} and {\it aligned loop} orbits. Horseshoes and nucleuphilic bananas support the existence of cuspy regions. Butterflies and aligned loops control the non-axisymmetric shape of outer regions. Without any need for central black holes, our distributed mass models resemble the nuclei of M31 and NGC4486B. It is also shown that the self-gravity of the stellar disc can prevent the double nucleus to collapse.

We introduce a class of eccentric discs with "strong" density cusps whose potentials are of St\"ackel form in elliptic coordinates. Our models exhibit some striking features: sufficiently close to the location of the cusp, the potential and surface density distribution diverge as $\Phi \propto r^{-1}$ and $\Sigma \propto r^{-2}$, respectively. As we move outward from the centre, the model takes a non-axisymmetric, lopsided structure. In the limit, when $r$ tends to infinity, the isocontours of $\Phi$ and $\Sigma$ become spherically symmetric. It is shown that the configuration space is occupied by three families of regular orbits: {\it eccentric butterfly}, {\it aligned loop} and {\it horseshoe} orbits. These orbits are properly aligned with the surface density distribution and can be used to construct self-consistent equilibrium states.