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.

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.

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.

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.

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

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.

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.

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.