Recent high-resolution N-body CDM simulations indicate that nonsingular three-parameter models such as the Einasto profile perform better than the singular two-parameter models, e.g. the Navarro, Frenk and White, in fitting a wide range of dark matter haloes. While many of the basic properties of the Einasto profile have been discussed in previous studies, a number of analytical properties are still not investigated. In particular, a general analytical formula for the surface density, an important quantity that defines the lensing properties of a dark matter halo, is still lacking to date. To this aim, we used a Mellin integral transform formalism to derive a closed expression for the Einasto surface density and related properties in terms of the Fox H and Meijer G functions, which can be written as series expansions. This enables arbitrary-precision calculations of the surface density and the lensing properties of realistic dark matter halo models. Furthermore, we compared the S\'ersic and Einasto surface mass densities and found differences between them, which implies that the lensing properties for both profiles differ.

Assuming the separable augmented density, it is always possible to construct a distribution function of a spherical population with any given density and anisotropy. We consider under what conditions the distribution constructed as such is in fact non-negative everywhere in the accessible phase-space. We first generalize known necessary conditions on the augmented density using fractional calculus. The condition on the radius part R(r^2) (whose logarithmic derivative is the anisotropy parameter) is equivalent to the complete monotonicity of R(1/w)/w. The condition on the potential part on the other hand is given by its derivative up to any order not greater than (3/2-beta) being non-negative where beta is the central anisotropy parameter. We also derive a specialized inversion formula for the distribution from the separable augmented density, which leads to sufficient conditions on separable augmented densities for the non-negativity of the distribution. The last generalizes the similar condition derived earlier for the generalized Cuddeford system to arbitrary separable systems.

The S\'ersic model is the de facto standard to describe the surface brightness distribution of hot stellar systems. An important inconvenience of this analytical model is that the corresponding luminosity density and associated properties cannot be expressed using elementary functions or even standard special functions. We present a set of compact and elegant analytical expressions for the luminosity density, cumulative luminosity and potential for the S\'ersic model in terms of the Fox H function for general values of the S\'ersic index. Furthermore, we present explicit series expansions of these quantities and discuss the asymptotic behaviour. Our analysis completes the work of Mazure & Capelato (2002) and Baes & Gentile (2011) and demonstrates the power of the underestimated Fox H function as a tool for analytical work.

Recently, some intriguing results have lead to speculations whether the central density slope -- velocity dispersion anisotropy inequality (An & Evans) actually holds at all radii for spherical dynamical systems. We extend these studies by providing a complete analysis of the global slope -- anisotropy inequality for all spherical systems in which the augmented density is a separable function of radius and potential. We prove that these systems indeed satisfy the global inequality if their central anisotropy is $\beta_0\leq 1/2$. Furthermore, we present several systems with $\beta_0 > 1/2$ for which the inequality does not hold, thus demonstrating that the global density slope -- anisotropy inequality is not a universal property. This analysis is a significant step towards an understanding of the relation for general spherical systems.

$N$-body simulations have unveiled several apparently universal properties of dark matter halos, including a cusped density profile, a power-law pseudo phase-space density $\rho/\sigma_r^3$, and a linear $\beta-\gamma$ relation between the density slope and the velocity anisotropy. We present a family of self-consistent phase-space distribution functions $F(E,L)$, based on the Dehnen-McLaughlin Jeans models, that incorporate these universal properties very accurately. These distribution functions, derived using a quadratic programming technique, are analytical, positive and smooth over the entire phase space and are able to generate four-parameter velocity anisotropy profiles $\beta(r)$ with arbitrary asymptotic values $\beta_0$ and $\beta_\infty$. We discuss the orbital structure of six radially anisotropic systems in detail and argue that, apart from its use for generating initial conditions for $N$-body studies, our dynamical modeling provides a valuable complementary approach to understand the processes involved in the formation of dark matter halos.

Both numerical simulations and observational evidence indicate that the outer regions of galaxies and dark matter haloes are typically mildly to significantly radially anisotropic. The inner regions can be significantly non-isotropic, depending on the dynamical formation and evolution processes. In an attempt to break the lack of simple dynamical models that can reproduce this behaviour, we explore a technique to construct dynamical models with an arbitrary density and an arbitrary anisotropy profile. We outline a general construction method and propose a more practical approach based on a parameterized anisotropy profile. This approach consists of fitting the density of the model with a set of dynamical components, each of which have the same anisotropy profile. Using this approach we avoid the delicate fine-tuning difficulties other fitting techniques typically encounter when constructing radially anisotropic models. We present a model anisotropy profile that generalizes the Osipkov-Merritt profile, and that can represent any smooth monotonic anisotropy profile. Based on this model anisotropy profile, we construct a very general seven-parameter set of dynamical components for which the most important dynamical properties can be calculated analytically. We use the results to look for simple one-component dynamical models that generate simple potential-density pairs while still supporting a flexible anisotropy profile. We present families of Plummer and Hernquist models in which the anisotropy at small and large radii can be chosen as free parameters. We also generalize these two families to a three-parameter family that self-consistently generates the set of Veltmann potential-density pairs. (Abridged...)

We present a method to investigate the radial stability of a spherical anisotropic system that hosts a central supermassive black hole (SBH). Such systems have never been tested before for stability, although high anisotropies have been considered in the dynamical models that were used to estimate the masses of the central putative supermassive black holes. A family of analytical anisotropic spherical Hernquist models with and without a black hole were investigated by means of N-body simulations. A clear trend emerges that the supermassive black hole has a significant effect on the overall stability of the system, i.e. an SBH with a mass of a few percent of the total mass of the galaxy can prevent or reduce the bar instabilities in anisotropic systems. Its mass not only determines the strength of the instability reduction, but also the time in which this occurs. These effects are most significant for models with strong radial anisotropies. Furthermore, our analysis shows that unstable systems with similar SBH but with different anisotropy radii evolve differently: highly radial systems become oblate, while more isotropic models tend to form into prolate structures. In addition to this study, we also present a Monte-Carlo algorithm to generate particles in spherical anisotropic systems.