We present the derivation of distribution functions for the first four members of a family of disks, previously obtained in (MNRAS, 371, 1873, 2006), which represent a family of axially symmetric galaxy models with finite radius and well behaved surface mass density. In order to do this we employ several approaches that have been developed starting from the potential-density pair and, essentially using the method introduced by Kalnajs (Ap. J., 205, 751, 1976) we obtain some distribution functions that depend on the Jacobi integral. Now, as this method demands that the mass density can be properly expressed as a function of the gravitational potential, we can do this only for the first four discs of the family. We also find another kind of distribution functions by starting with the even part of the previous distribution functions and using the maximum entropy principle in order to find the odd part and so a new distribution function, as it was pointed out by Dejonghe (Phys. Rep., 133, 217, 1986). The result is a wide variety of equilibrium states corresponding to several self-consistent finite flat galaxy models.

Spherical systems with polytropic equations of state are of great interest in astrophysics. They are widely used to describe neutron stars, red giants, white dwarfs, brown dwarfs, main sequence stars, galactic halos, and globular clusters of diverse sizes. In this paper we construct analytically a family of self-gravitating spherical models in the post-Newtonian approximation of general relativity. These models present interesting cusps in their density profiles which are appropriate for the modeling of galaxies and dark matter halos. The systems described here are anisotropic in the sense that their equiprobability surfaces in velocity space are nonspherical, leading to an overabundance of radial or circular orbits, depending on the parameters of the model under consideration. Among the family of models, we find the post-Newtonian generalization of the Plummer and Hernquist models. A close inspection of their equation of state reveals that these solutions interpolate smoothly between a polytropic sphere in the asymptotic region and an inner core that resembles an isothermal sphere. Finally, we study the thermodynamics of these models and argue for their stability.

In this paper we study the kinetic theory of many-particle astrophysical systems imposing axial symmetry and extending our previous analysis in Phys. Rev. D 83, 123007 (2011). Starting from a Newtonian model describing a collisionless self-gravitating gas, we develop a framework to include systematically the first general relativistic corrections to the matter distribution and gravitational potentials for general stationary systems. Then, we use our method to obtain particular solutions for the case of the Morgan & Morgan disks. The models obtained are fully analytical and correspond to the post-Newtonian generalizations of classical ones. We explore some properties of the models in order to estimate the importance of post-Newtonian corrections and we find that, contrary to the expectations, the main modifications appear far from the galaxy cores. As a by-product of this investigation we derive the corrected version of the tensor virial theorem. For stationary systems we recover the same result as in the Newtonian theory. However, for time dependent backgrounds we find that there is an extra piece that contributes to the variation of the inertia tensor.

In this paper we study the kinetic theory of many-particle astrophysical systems and we present a consistent version of the collisionless Boltzmann equation in the 1PN approximation. We argue that the equation presented by Rezania and Sobouti in A&A 354 1110 (2000) is not the correct expression to describe the evolution of a collisionless self-gravitating gas. One of the reasons that account for the previous statement is that the energy of a free-falling test particle, obeying the 1PN equations of motion for static gravitational fields, is not a static solution of the mentioned equation. The same statement holds for the angular momentum, in the case of spherical systems. We provide the necessary corrections and obtain an equation that is consistent with the corresponding equations of motion and the 1PN conserved quantities. We suggest some potential relevance for the study of high density astrophysical systems and as an application we construct the corrected version of the post-Newtonian polytropes.

We study the motion of test particles around a center of attraction represented by a monopole (with and without spheroidal deformation) surrounded by a ring, given as a superposition of Morgan & Morgan discs. We deal with two kinds of bounded orbits: (i) Equatorial circular orbits and (ii) general three-dimensional orbits. The first case provides a method to perform a linear stability analysis of these structures by studying the behavior of vertical and epicyclic frequencies as functions of the mass ratio, the size of the ring and/or the quadrupolar deformation. In the second case, we study the influence of these parameters in the regularity or chaoticity of motion. We find that there is a close connection between linear stability (or unstability) of equatorial circular orbits and regularity (or chaoticity) of the three-dimensional motion.

We present the first analytical superposition of a charged black hole with an annular disk of extremal dust. In order to obtain the solutions, we first solve the Einstein-Maxwell field equations for sources that represent disk-like configurations of matter in confomastatic spacetimes by assuming a functional dependence among the metric function, the electric potential and an auxiliary function,which is taken as a solution of the Laplace equation. We then employ the Lord Kelvin Inversion Method applied to models of finite extension in order to obtain annular disks. The structures obtained extend to infinity, but their total masses are finite and all the energy conditions are satisfied. Finally, we observe that the extremal Reissner-Nordstr\"{o}m black hole can be embedded into the center of the disks by adding a boundary term in the inversion.

A new formalism is presented for finding equilibrium distribution functions for axisymmetric systems. The formalism, obtainded by using the concept of fractional derivatives, generalizes the methods of Fricke (1952), Kalnajs (1972) and Jiang & Ossipkov (2007), and has the advantage that can be applied to a wider variety of models. We found that this approach can be applied both to tridimensional systems and to flat systems, without the necessity of dealing with a pseudo-volume mass density. As an application, we obtain the distribution functions of the Binney's logarithmic model and of the Mestel disc.

We present the formulation of a new infinite family of self-consistent stellar models, designed to describe axisymmetric flat galaxies. The corresponding density-potential pair is obtained as a superposition of members belonging to the generalized Kalnajs family, by imposing the condition that the density can be expressed as a regular function of the gravitational potential, in order to derive analytically the corresponding equilibrium distribution functions (DF). The resulting models are characterized by a well-behaved surface density, as in the case of generalized Kalnajs discs. Then, we present a study of the kinematical behavior which reveals, in some particular cases, a very satisfactory behavior of the rotational curves (without the assumption of a dark matter halo). We also analyze the equatorial orbit's stability and Poincare surfaces of section are performed for the 3-dimensional orbits. Finally, we obtain the corresponding equilibrium DFs, using the approaches introduced by Kalnajs (Ap. J., 205, 751, 1976) and Dejonghe (Phys. Rep., 133 (3-4), 217, 1986).