From a Parametrized Post-Newtonian (PPN) [46-51] perspective, we address the question of whether or not the new degrees of freedom, represented by the PPN potentials, can lead to significant modifications in the dynamics of galaxies in the direction of rendering dark matter obsolete. Here, we focus on the study of rotation curves associated with spherically symmetric configurations. The values for the post-Newtonian parameters, which help us to classify the different metric theories of gravity, are tightly constrained mainly by solar system experiments [47,48]. However, in recent years, metric theories characterized by screening mechanisms [4,10,22,25] have become popular, due to the fact that they lead to the possibility of modifications in larger scales than the solar system while retaining the success of GR on it, allowing for violations of the constraints of the Post-Newtonian parameters. In such a context, we consider here two kinds of solutions for field equations: (i) Vacuum solutions and (ii) in the presence of a politropic distribution of matter. For the case (i) we find that the post-Newtonian corrections do not lead to modifications significant enough to be considered an alternative to the dark matter hypothesis. In the case (ii) we find that for a wide range of values for the PPN parameters {\gamma}, {\beta}, {\xi}, {\alpha}_3 , {\zeta}_1 and {\zeta}_2 , the need for dark matter is unavoidable, in order to find flat rotation curves. It is only for theories in which {\zeta}_3>0 that some resemblance of flat rotation curves is found. The latter suggests, at least for the models considered, that these are the only theories capable of replacing dark matter as a possible explanation for the dynamics of galaxies.

We recently found that regular orbits in axially symmetric galactic disks have their envelopes $Z(R)$ accurately described by the relation $Z(R)\propto[\Sigma_I(R)]^{-1/3}$, if their amplitudes are comparable to the disk thickness, where $\Sigma_I$ is the surface density of the disk (integrated over its whole vertical range). Moreover, the usual adiabatic approximation gives a good description of the orbits' envelopes for low vertical amplitudes. However, these two approaches are apparently disconnected, since their expressions differ qualitatively. Our purpose in this paper is to fill this gap by extending these previous formulae to regular orbits with arbitrary vertical amplitudes inside the disk. We compare existing $Z(R)$ estimates: the razor-thin disk case, the adiabatic approximation (low-amplitude orbits in three-dimensional disks), and the integrated surface-density estimate (high-amplitude orbits in three-dimensional disks) in order to establish a connection between them. The formula presented here links the aforementioned results in an elegant and continuous way, being valid for vertical amplitudes throughout the whole vertical extension of the disk and with an expression which has the same form for all regimes. The advantage of the present formalism is the dependence of $Z(R)$ only on observable quantities, namely the disk's vertically integrated surface density, without the need to obtain the gravitational potential for the system.

By considering the Einstein-Vlasov system for static spherically symmetric distributions of matter, we show that configurations with constant anisotropy parameter $\beta$ have, necessarily, a distribution function (DF) of the form $F=l^{-2\beta}\xi(\varepsilon)$, where $\varepsilon=E/m$ and $l=L/m$ are the relativistic energy and angular momentum per unit rest mass, respectively. We exploit this result to obtain DFs for the general relativistic extension of the Hypervirial family introduced by Nguyen and Lingam (2013), which Newtonian potential is given by $\phi(r)=-\phi_o /[1+(r/a)^{n}]^{1/n}$ ($a$ and $\phi_o$ are positive free parameters, $n=1,2,...$). Such DFs can be written in the form $F_{n}=l^{n-2}\xi_{n}(\varepsilon)$. For odd $n$, we find that $\xi_n$ is a polynomial of order $2n+1$ in $\varepsilon$, as in the case of the Hernquist model ($n=1$), for which $F_1\propto l^{-1}\left(2\varepsilon-1\right)\left(\varepsilon-1\right)^2$. For even $n$, we can write $\xi_n$ in terms of incomplete beta functions (Plummer model, $n=2$, is an example). Since we demand that $F\geq 0$ throughout the phase space, the particular form of each $\xi_n$ leads to restrictions for the values of $\phi_o$. For example, for the Hernquist model we find that $0\leq \phi_o \leq2/3$, i.e. an upper bounding value less than the one obtained for Nguyen and Lingam ($0\leq \phi_o \leq1$), based on energy conditions.

During the last few decades, there has been a growing interest in exact solutions of Einstein equations describing razor-thin disks. Despite the progress in the area, the analytical study of geodesic motion crossing the disk plane in these systems is not yet so developed. In the present work, we propose a definite vertical stability criterion for circular equatorial timelike geodesics in static, axially symmetric thin disks, possibly surrounded by other structures preserving axial symmetry. It turns out that the strong energy condition for the disk stress-energy content is sufficient for vertical stability of these orbits. Moreover, adiabatic invariance of the vertical action variable gives us an approximate third integral of motion for oblique orbits which deviate slightly from the equatorial plane. Such new approximate third integral certainly points to a better understanding of the analytical properties of these orbits. The results presented here, derived for static spacetimes, may be a starting point to study the motion around rotating, stationary razor-thin disks. Our results also allow us to conjecture that the strong energy condition should be sufficient to assure transversal stability of periodic orbits for any singular timelike hypersurface, provided it is invariant under the geodesic flow.

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.

We present an analytical simple formula for an approximated third integral of motion associated with nearly equatorial orbits in the Galaxy: $I_{3}=Z\Sigma_{I}^{1/3}$, where $Z(R)$ is the vertical amplitude of the orbit at galactocentric distance $R$ and $\Sigma_{I}(R)$ is the integrated dynamical surface mass density of the disk, a quantity which has recently become measurable. We also suggest that this relation is valid for disk-crossing orbits in a wide variety of axially symmetric galactic models, which range from razor-thin disks to disks with non-negligible thickness, whether or not the system includes bulges and halos. We apply our formalism to a Miyamoto-Nagai model and to a realistic model for the Milky Way. In both cases, the results provide fits for the shape of nearly equatorial orbits which are better than the corresponding curves obtained by the usual adiabatic approximation when the orbits have vertical amplitudes comparable to the disk's scale height. We also discuss the role of this approximate third integral of motion in modified theories of gravity.

We consider the three-dimensional bounded motion of a test particle around razor-thin disk configurations, by focusing on the adiabatic invariance of the vertical action associated with disk-crossing orbits. We find that it leads to an approximate third integral of motion predicting envelopes of the form $Z(R)\propto[\Sigma(R)]^{-1/3}$, where $R$ is the radial galactocentric coordinate, $Z$ is the z-amplitude (vertical amplitude) of the orbit and $\Sigma$ represents the surface mass density of the thin disk. This third integral, which was previously formulated for the case of flattened 3D configurations, is tested for a variety of trajectories in different thin-disk models.

The study of off-equatorial orbits in razor-thin disks is still in its beginnings. Contrary to what was presented in the literature in recent publications, the vertical stability criterion for equatorial circular orbits cannot be based on the vertical epicyclic frequency, because of the discontinuity in the gravitational field on the equatorial plane. We present a rigorous criterion for the vertical stability of circular orbits in systems composed by a razor-thin disk surrounded by a smooth axially symmetric distribution of matter, the latter representing additional structures such as thick disk, bulge and (dark matter) halo. This criterion is satisfied once the mass surface density of the thin disk is positive. Qualitative and quantitative analyses of nearly equatorial orbits are presented. In particular, the analysis of nearly equatorial orbits allows us to construct an approximate analytical third integral of motion in this region of phase-space, which describes the shape of these orbits in the meridional plane.

We make a revision of the stability criteria for equatorial circular orbits, obtained from the epicyclic approximation, which is widely used in Newtonian models for axisymmetric galaxies. We find that, for the case of thin disk models, the criteria of vertical stability must be reformulated, due to the discontinuity in the gravitational field. We show that, for a model characterized by a surface mass density $\Sigma$, the necessary and sufficient condition to have vertically stable orbits is that $\Sigma>0$. On the other hand, the criteria for radial stability is the same as in thick diks, i.e. that the epicyclic frequency is positive.

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.

Finite thin disc models of four galaxies in the Ursa Major cluster are presented. The models are obtained by means of the Hunter method and the particular solutions are choosen in such a way that the circular velocities are adjusted very accurately to the observed rotation curves of some specific spiral galaxies. We present particular models for the four galaxies NGC3877, NGC3917, NGC3949 and NGC4010 with data taken from the recent paper by Verheijen & Sancici (2001). By integrating the corresponding surface mass densities, we obtain the total mass M of these four galaxies, all of them being of the order of 10^10 solar masses. These obtained values for M may be taken as a quite accurately estimative of the mass upper bound of these galaxies, since in the model was considered that all their mass was concentrated at the galactic disc. The models can be consider as a first approximation to the obtaining of quite realistic models of spiral galaxies

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).

The motion of test particles in the gravitational fields generated by the first four members of the infinite family of generalized Kalnajs discs, is studied. In first instance, we analyze the stability of circular orbits under radial and vertical perturbations and describe the behavior of general equatorial orbits and so we find that radial stability and vertical instability dominate such disc models. Then we study bounded axially symmetric orbits by using the Poincare surfaces of section and Lyapunov characteristic numbers and find chaos in the case of disc-crossing orbits and completely regular motion in other cases.