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 review our contribution to infrared Renormalization Group (RG) effects to General Relativity in the context of galaxies. Considering the effective action approach to Quantum Field Theory in curved background, we argued that the proper RG energy scale, in the weak field limit, should be related to the Newtonian potential. In the galaxy context, even without dark matter, this led to a remarkably small gravitational coupling G variation (about or less than 10^{-12} of its value per light-year), while also capable of generating galaxy rotation curves about as good as the best phenomenological dark matter profiles (considering both the rotation curve shape and the expected mass-to-light ratios). Here we also comment on related developments, open issues and perspectives.

A method is presented to bend a thin massive line when the curvature is small. The procedure is applied to a homogeneous thin bar with two types of curvatures. One of them mimics a galactic bar with two spiral arms at its tips. It is showed that if the bending function is a linear combination of Legendre polynomials, then the bent potential is an exact solution of the Laplace equation. A transformation is applied on the thin bent bars to generate three-dimensional potential-density pairs without singularities. The potentials of the thin bent bars are also used to generate non-axisymmetric planar distributions of matter.

An identity that relates multipolar solutions of the Einstein equations to Newtonian potentials of bars with linear densities proportional to Legendre polynomials is used to construct analytical potential-density pairs of infinitesimally thin bars with a given linear density profile. By means of a suitable transformation, softened bars that are free of singularities are also obtained. As an application we study the equilibrium points and stability for the motion of test particles in the gravitational field for three models of rotating bars.

We consider the application of quantum corrections computed using renormalization group arguments in the astrophysical domain and show that, for the most natural interpretation of the renormalization group scale parameter, a gravitational coupling parameter $G$ varying $10^{-7}$ of its value across a galaxy (which is roughly a variation of $10^{-12}$ per light-year) is sufficient to generate galaxy rotation curves in agreement with the observations. The quality of the resulting fit is similar to the Isothermal profile quality once both the shape of the rotation curve and the mass-to-light ratios are considered for evaluation. In order to perform the analysis, we use recent high quality data from nine regular disk galaxies. For the sake of comparison, the same set of data is modeled also for the Modified Newtonian Dynamics (MOND) and for the recently proposed Scalar Tensor Vector Gravity (STVG). At face value, the model based on quantum corrections clearly leads to better fits than these two alternative theories.

The complex-shift method is applied to the Kuzmin-Toomre family of discs to generate a family of non-axisymmetric flat distributions of matter. These are then superposed to construct non-axisymmetric flat rings. We also consider triaxial potential-density pairs obtained from these non-axisymmetric flat systems by means of suitable transformations. The use of the imaginary part of complex-shifted potential-density pairs is also discussed.

The Kuzmin-Toomre family of discs is used to construct potential-density pairs that represent flat ring structures in terms of elementary functions. Systems composed of two concentric flat rings, a central disc surrounded by one ring and a ring with a centre of attraction are also presented. The circular velocity of test particles and the epicyclic frequency of small oscillations about circular orbits are calculated for these structures. A few examples of three-dimensional potential-density pairs of "inflated" flat rings (toroidal mass distributions) are presented.

We study the effect of noise on the evolution of the growth factor of density perturbations in the context of the LCDM model. Stochasticity is introduced as a Wiener process amplified by an intensity parameter alpha. By comparing the evolution of deterministic and stochastic cases for different values of alpha we estimate the intensity level necessary to make noise relevant for cosmological tests based on large-scale structure data. Our results indicate that the presence of random forces underlying the fluid description can lead to significant deviations from the nonstochastic solution at late times for alpha>0.001.

The advance of perihelion for geodesic motion on the galactic plane of some exact general relativistic disc solutions is calculated. Approximate analytical and numerical results are presented for the static Chazy-Curzon and the Schwarzschild discs in Weyl coordinates, the Schwarzschid disc in isotropic coordinantes and the stationary Kerr disc in the Weyl-Lewis-Papapetrou metrics. It is found that for these disc models the advance of perihelion may be an increasing or decreasing function of the orbital excentricity. The precession due to Newtonian gravity for these disc models is also calculated.

An analytical solution of Einstein-Maxwell equations with a static fluid as a source is presented. The spacetime is represented by the axially symmetric Weyl metric and the energy-momentum tensor describes a coupling of a fluid with an electromagnetic field. When appropriate limits are performed we recover the well-known solutions of Gutsunaev-Manko and Schwarzschild. Also, using Eckart's thermodynamics, we calculated the temperature, the mechanical pressure, the charge density and the energy density of the system. The analysis of thermodynamic quantities suggests that the solution can be used to represent a magnetized compact stellar object surrounded by a charged fluid.

A geometric flow based in the Riemann-Christoffel curvature tensor that in two dimensions has some common features with the usual Ricci flow is presented. For $n$ dimensional spaces this new flow takes into account all the components of the intrinsic curvature. For four dimensional Lorentzian manifolds it is found that the solutions of the Einstein equations associated to a "detonant" sphere of matter, as well, as a Friedman-Roberson-Walker cosmological model are examples of Riemann-Christoffel flows. Possible generalizations are mentioned.

The stability of a recently proposed general relativistic model of galaxies is studied in some detail. This model is a general relativistic version of the well known Miyamoto-Nagai model that represents well a thick galactic disk. The stability of the disk is investigated under a general first order perturbation keeping the spacetime metric frozen (no gravitational radiation is taken into account). We find that the stability is associated with the thickness of the disk. We have that flat galaxies have more not-stable modes than the thick ones i.e., flat galaxies have a tendency to form more complex structures like rings, bars and spiral arms.

We find a family of exact solutions of Einstein equations describing the field of a static axisymmetric thin disk living in six-dimensional space-time. In particular, we study the disks constructed (by cutting out the central part of the spacetime) from the conventional Schwarzschild and Chazy-Curzon solutions with simple extensions in the extra 2 dimensions. The disks are interpreted in terms of two counter-rotating streams of particles on free circular orbits. Two extra parameters -- the constants of motion resulting from projections of the particle-velocities' extra components -- are constrained by the requirement that the orbits within the disk be stable. The requirement is met just in cases when the radial profile of the disk orbital speed fits the rotation curves of many spiral galaxies. It thus turns out that the effective modification of the gravity law by extra dimensions could explain the observed flatness of these curves equally well as the usually given dark-matter interpretation. In the second part of the paper we show that the inclusion of extra dimensions also leads to better fits of the gravitational-lensing data for galaxy clusters, without changing results obtained in solar-system scales. Finally, we discuss whether the effect of extra dimensions could also be translated as following from the occurrence of extra matter ("Kaluza-Klein particles"). A comment on possible properties of such particles and a chance to detect them at LHC is given in the appendix.

Pairs potential-density in terms of elementary functions that represents flat rings structures are presented. We study structures representing one or several concentric flat rings. Also disks surrounded by concentric flat rings are exhibited. The stability of concentrically circular orbits of particles moving on a flat ring structure is analyzed for radial perturbations.

We present two triaxial analytical potential-density pairs that can be viewed as generalized versions of the axisymmetric Miyamoto and Nagai and Satoh galactic models. These potential-density pairs may be useful models for galaxies with box-shaped bulges. The resulting mass density distributions are everywhere non-negative and free from singularities. Also, a few numerically calculated orbits for the Miyamoto and Nagai-like triaxial potential are presented.

The time evolution of a set of 22 Mo unstable charged stars that collapse is computed integrating the Einstein-Maxwell equations. The model simulate the collapse of an spherical star that had exhausted its nuclear fuel and have or acquires a net electric charge in its core while collapsing. When the charge to mass ratio is Q/M >= 1 the star do not collapse and spreads. On the other hand, it is observed a different physical behavior with a charge to mass ratio 1 > Q/ M > 0.1. In this case, the collapsing matter forms a bubble enclosing a lower density core. We discuss an immediate astrophysical consequence of these results that is a more efficient neutrino trapping during the stellar collapse and an alternative mechanism for powerful supernova explosions. The outer space-time of the star is the Reissner-Nordstrom solution that match smoothly with our interior numerical solution, thus the collapsing models forms Reissner-Nordstrom black holes.

A class of exact general relativistic rotating thick disks surrounded by rotating dust is constructed by applying the ``displace, cut, fill and reflect'' method on the van Stockum class of metrics. Two particular solutions are considered in detail. We find that the disks have annular regions with negative energy density and heat flow is present in most of the radial extension of the disk. It is shown that thin disks are composed of exotic matter with an equation of state of cosmic strings or struts. We also comment the recent relativistic galactic model proposed by Cooperstock and Tieu and the presence of an additional disk of exotic matter in their model.

We study the stability of general relativistic static thick disks, as an application we consider the thick disk generated by applying the ``displace, cut, fill and reflect" method, usually known as the image method, to the Schwarzschild metric in isotropic coordinates. The isotropic Schwarzschild thick disk obtained from this method is the simplest model to describe, in the context of General Relativity, real thick galaxies. The stability under a general first order perturbation of the disk is investigated. The first order perturbation, when applying to the conservation equations, leads to a set of differential equations that are, in general, not self-consistent. This opens the possibility of performing various kinds of perturbations to transform the resulting system of equations into a self-consistent system. We perform a complete classification of the perturbations as well as the stability analysis for all the relevant physical perturbations. We found that, in general, the isotropic Schwarzschild thick disk is stable under these kinds of perturbations.

A pseudo-Newtonian Hill problem based on the Paczynski-Wiita pseudo-Newtonian potential that reproduces general relativistic effects is presented and compared with the usual Newtonian Hill problem. Poincare maps, Lyapunov exponents and fractal escape techniques are employed to study bounded and unbounded orbits. In particular we consider the systems composed by Sun, Earth and Moon and composed by the Milky Way, the M2 cluster and a star. We find that some pseudo-Newtonian systems - including the M2 system - are more stable than their Newtonian equivalent.

We present analytical potential-density pairs in three dimensions for the gravitational field of galaxies, obtained by thickening the multipolar expansion up to the quadrupole term. These may be interpreted as generalizations of the Miyamoto and Nagai potential-density pairs. With a suitable restriction on the possible values of the multipole moments, the density distributions are positive and monotone decreasing functions of the radial and axial coordinates.

In this comment we question some arguments presented in astro-ph/0512048 to refuse the presence of an singular mass surface layer. In particular, incorrect expressions are used for the disk's surface mass density. We also point out that the procedure of removing the descontinuity on the $z=0$ plane with a region of continuous density gradient generates other two regions of descontinuities with singular mass surface layers making the model unrealistic.

We describe an implicit procedure for solving linear equation systems resulting from the discretization of the three dimensional (seven variables) linear Fokker-Planck equation. The discretization of the Fokker-Planck equation is performed using a twenty-five point molecule that leads to a coefficient matrix with equal number of diagonals. The method is an extension of Stone's implicit procedure, includes a vast class of collision terms and can be applied to stationary or non stationary problems with different discretizations in time. Test calculations and comparisons with other methods are presented in two stationary examples, including an astrophysical application for the Miyamoto-Nagai disk potential for a typical galaxy.

We analyze the presence of an additional singular thin disk in the recent General Relativistic model of galactic gravitational field proposed by Cooperstock and Tieu. The physical variables of the disk's energy-momentum tensor are calculated. We show that the disk is made of exotic matter, either cosmic strings or struts with negative energy density.

A special form of the isotropic metric in cylindrical coordinates is used to construct what may be interpreted as the General Relativistic versions of some wellknown potential-density pairs used in Newtonian gravity to model three-dimensional distributions of matter in galaxies. The components of the energy-momentum tensor are calculated for the first two Miyamoto-Nagai potentials and a particular potential due to Satoh. The three potentials yield distributions of matter in which all tensions are pressures and all energy conditions are satisfied for certain ranges of the free parameters. A few non-planar geodesic orbits are computed for one of the potentials and compared with the Newtonian case. Rotation is also incorporated to the models and the effects of the source rotation on the rotation profile are calculated as first order corrections by using an approximate form of the Kerr metric in isotropic coordinates.

We find a numerical self-consistent stellar model by finding the distribution function of a thin disk that satisfies simultaneously the Fokker-Planck and Poisson equations. The solution of the Fokker-Planck equation is found by a direct numerical solver using finite differences and a variation of Stone's method. The collision term in the Fokker-Planck equation is found using the local approximation and the Rosenbluth potentials. The resulting diffusion coefficients are explicitly evaluated using a Maxwellian distribution for the field stars. As a paradigmatic example, we apply the numerical formalism to find the distribution function of a Kuzmin-Toomre thin disk. This example is studied in some detail showing that the method applies to a large family of actual galaxies.

Contrary to a recent claim by Anderson ["The Mathematical Theory of Cosmic Strings", I.O.P. Publishing, Bristol 2003], we show that the Letelier-Gal'tsov metric does represent a system of crossed straight infinite cosmic strings moving with arbitrary constant velocities.

New families of exact general relativistic thick disks are constructed using the ``displace, cut, fill and reflect'' method. A class of functions used to ``fill'' the disks is derived imposing conditions on the first and second derivatives to generate physically acceptable disks. The analysis of the function's curvature further restrict the ranges of the free parameters that allow phisically acceptable disks. Then this class of functions together with the Schwarzschild metric is employed to construct thick disks in isotropic, Weyl and Schwarzschild canonical coordinates. In these last coordinates an additional function must be added to one of the metric coefficients to generate exact disks. Disks in isotropic and Weyl coordinates satisfy all energy conditions, but those in Schwarzschild canonical coordinates do not satisfy the dominant energy condition.

An exact but simple general relativistic model for the gravitational field of active galactic nuclei is constructed, based on the superposition in Weyl coordinates of a black hole, a Chazy-Curzon disk and two rods, which represent matter jets. The influence of the rods on the matter properties of the disk and on its stability is examined. We find that in general they contribute to destabilize the disk. Also the oscillation frequencies for perturbed circular geodesics on the disk are computed, and some geodesic orbits for the superposed metric are numerically calculated.

The stability of general relativistic thin disks is investigated under a general first order perturbation of the energy momentum tensor. In particular, we consider temporal, radial and azimuthal "test matter" perturbations of the quantities involved on the plane $z=0$. We study the thin disks generated by applying the "displace, cut and reflect" method, usually known as the image method, to the Schwarzschild metric in isotropic coordinates and to the Chazy-Curzon metric and the Zipoy-Voorhees metric ($\gamma$-metric) in Weyl coordinates. In the case of the isotropic Schwarzschild thin disk, where a radial pressure is present to support the gravitational attraction, the disk is stable and the perturbation favors the formation of rings. Also, we found the expected result that the thin disk models generated by the Chazy-Curzon and Zipoy-Voorhees metric with only azimuthal pressure are not stable under a general first order perturbation

Using the well-known ``displace, cut and reflect'' method used to generate disks from given solutions of Einstein field equations, we construct static charged disks made of perfect fluid based on the Reissner-Nordstr\"{o}m solution in isotropic coordinates. We also derive a simple stability condition for charged and non charged perfect fluid disks. As expected, we find that the presence of charge increases the regions of instability of the disks.

A numerical solution of Einstein field equations for a spherical symmetric and stationary system of identical and auto-gravitating particles in phase transition is presented. The fluid possess a perfect fluid energy momentum tensor, and the internal interactions of the system are represented by a van der Walls like equation of state able to describe a first order phase transition of the type gas-liquid. We find that the space-time curvature, the radial component of the metric, and the pressure and density show discontinuities in their radial derivatives in the phase coexistence region. This region is found to be a spherical surface concentric with the star and the system can be thought as a foliation of acronal, concentric and isobaric surfaces in which the coexistence of phases occurs in only one of these surfaces. This kind of system can be used to represent a star with a high energy density core and low energy density mantle in hydrodynamic equilibrium.

A method to construct exact general relativistic thick disks that is a simple generalization of the ``displace, cut and reflect'' method commonly used in Newtonian, as well as, in Einstein theory of gravitation is presented. This generalization consists in the addition of a new step in the above mentioned method. The new method can be pictured as a ``displace, cut, {\it fill} and reflect'' method. In the Newtonian case, the method is illustrated in some detail with the Kuzmin-Toomre disk. We obtain a thick disk with acceptable physical properties. In the relativistic case two solutions of the Weyl equations, the Weyl gamma metric (also known as Zipoy-Voorhees metric) and the Chazy-Curzon metric are used to construct thick disks. Also the Schwarzschild metric in isotropic coordinates is employed to construct another family of thick disks. In all the considered cases we have non trivial ranges of the involved parameter that yield thick disks in which all the energy conditions are satisfied.

The Rayleigh criterion is used to study the stability of circular orbits of particles moving around static black holes surrounded by different axially symmetric structures with reflection symmetry, like disks, rings and halos. We consider three models of disks one of infinite extension and two finite, and one model of rings. The halos are represented by external quadrupole moments (either oblate or prolate). Internal quadrupole perturbation (oblate and prolate) are also considered. For this class of disks the counter-rotation hypothesis implies that the stability of the disks is equivalent to the stability of test particles. The stability of Newtonian systems is also considered and compared with the equivalent relativistic situation. We find that the general relativistic dynamics favors the formation of rings.

Using the well-known ``displace, cut and reflect'' method used to generate disks from given solutions of Einstein field equations, we construct static disks made of perfect fluid based on vacuum Schwarzschild's solution in isotropic coordinates. The same method is applied to different exactsolutions to the Einstein'sequations that represent static spheres of perfect fluids. We construct several models of disks with axially symmetric perfect fluid halos. All disks have some common features: surface energy density and pressures decrease monotonically and rapidly with radius. As the ``cut'' parameter $a$ decreases, the disks become more relativistic, with surface energy density and pressure more concentrated near the center. Also regions of unstable circular orbits are more likely to appear for high relativistic disks. Parameters can be chosen so that the sound velocity in the fluid and the tangential velocity of test particles in circular motion are less then the velocity of light. This tangential velocity first increases with radius and reaches a maximum.

We present a novel method to study interacting orbits in a fixed mean gravitational field associated with a solution of the Einstein field equations. The idea is to consider the Newton gravity among the orbiting particles in a geometry given by the main source. We apply the technique in the of study two and three self-gravitating particles moving around a black hole, i.e., in a Schwarzschild geometry. We also compare with the equivalent Newtonian problem and noted differences in the structural stability, e.g., binary systems were found only in the general relativistic approach.

The Melnikov method is applied to periodically perturbed open systems modeled by an inverse--square--law attraction center plus a quadrupolelike term. A compactification approach that regularizes periodic orbits at infinity is introduced. The (modified) Smale-Birkhoff homoclinic theorem is used to study transversal homoclinic intersections. A larger class of open systems with degenerated (nonhyperbolic) unstable periodic orbits after regularization is also briefly considered.

The exact solution to the Einstein equations that represents a static axially symmetric source deformed by an internal quadrupole is considered. By using the Poincare section method we numerically study the geodesic motion of test particles. For the prolate quadrupolar deformations we found chaotic motions contrary to the oblate case where only regular motion is found. We also study the metric that represents a rotating black hole deformed by a quadrupolar term. This metric is obtained as a two soliton solution in the context of Belinsky--Zakharov inverse scattering method. The stability of geodesics depends strongly on the relative direction of the spin of the center of attraction and the test particle angular momentum. The rotation does not alter the regularity of geodesic motions in the oblate case, i.e., the orbits in this case remain regular. We also employ the method of Lyapounov characteristic numbers to examine the stability of orbits evolving around deformed nonrotating centers of attraction. The typical time to observe instability of orbits is analyzed.

The motion of particles in the field of forces associated to an axially symmetric attraction center modeled by a monopolar term plus a prolate quadrupole deformation is studied using Poincare surface of sections and Lyapunov characteristic numbers. We find chaotic motion for certain values of the parameters, and that the instability of the orbits increases when the quadrupole parameter increases. A general relativistic analogue is briefly discussed.

The Newtonian as well as the special relativistic dynamics are used to study the stability of orbits of a test particle moving around a black hole plus a dipolar halo. The black hole is modeled by either the usual monopole potential or the Paczynki-Wiita pseudo-Newtonian potential. The full general relativistic similar case is also considered. The Poincare section method and the Lyapunov characteristic exponents show that the orbits for the pseudo-Newtonian potential models are more unstable than the corresponding general relativistic geodesics.

By considering full-field string solutions of the Abelian--Higgs model, we modify the model of a fluid of strings (which is composed of Nambu strings) to obtain a model for a ``fluid of vortices.'' With this model, and following closely Soleng's proposal of a fluid of strings as the source of a Milgrom-type correction to the Newton dynamics, we determine quantitatively the modified dynamics generated by a static, spherical fluid of vortices.

We raise the hypothesis that the density fluctuations field which originates the growth of large scale structures is a combination of two or more distributions. By applying the statistical analysis of finite mixture distributions to a specific combination of Gaussian plus non-Gaussian random fields, we studied the case where just a small departure from Gaussianity is allowed. Our results suggest that even a very small level of non-Gaussianity may introduce significant changes in the cluster abundance evolution rate.

The well-known ``displace, cut, and reflect'' method used to generate cold disks from given solutions of Einstein equations is extended to solutions of Einstein-Maxwell equations. Four exact solutions of the these last equations are used to construct models of hot disks with surface density, azimuthal pressure, and azimuthal current. The solutions are closely related to Kerr, Taub-NUT, Lynden-Bell-Pinault and to a one-soliton solution. We find that the presence of the magnetic field can change in a nontrivial way the different properties of the disks. In particular, the pure general relativistic instability studied by Bicak, Lynden-Bell and Katz [Phys. Rev. D47, 4334, 1993] can be enhanced or cured by different distributions of currents inside the disk. These currents, outside the disk, generate a variety of axial symmetric magnetic fields. As far as we know these are the first models of hot disks studied in the context of general relativity.

Using a Newtonian approximation, we developed a quantitative criterion for the collapse of a spherical distribution of matter under an isolated texture field. In particular, we found that the evolution of an overdense region is strongly determined by two parameters: the energy scale of symmetry breaking and the initial radius of the system.

We make a detailed analysis of Newtonian as well as relativistic core-shell models recently proposed to describe a black hole or neutron star surrounded by shells of matter, and in a seminal sense also galaxies, supernovae and star remnants since there are massive shell-like structures surrounding many of them and also evidences for many galactic nuclei hiding black holes. We discuss the unicity of the models in relation to their analyticity at the black hole horizon and also to the full elimination of conical singularities. Secondly, we study the role played by the presence/lack of discrete reflection symmetries about equatorial planes in the chaotic behavior of the orbits, which is to be contrasted with the almost universal acceptance of reflection symmetries as default assumptions in galactic modeling. We also compare the related effects if we change a true central black hole by a Newtonian central mass. The numerical findings are: 1- The breakdown of the reflection symmetry about the equatorial plane in both Newtonian and relativistic core-shell models does i) enhance in a significant way the chaoticity of orbits in reflection symmetric oblate shell models and ii) inhibit significantly also the occurrence of chaos in reflection symmetric prolate shell models. In particular, in the prolate case the lack of the reflection symmetry provides the phase space with a robust family of regular orbits that is otherwise not found at higher energies. 2- The relative extents of the chaotic regions in the relativistic cases (i. e. with a true central black hole) are significantly larger than in the corresponding Newtonian ones (which have just a $-1/r$ central potential).