Gravitational field of slightly deformed naked singularities

We derive a particular approximate solution of Einstein equations, describing the gravitational field of a mass distribution that slightly deviates from spherical symmetry. The deviation is described by means of a quadrupole parameter that is responsible for the appearance of a curvature singularity, which is not covered by a horizon. We investigate the motion of test particles in the gravitational field of this naked singularity and show that the quadrupole parameter affects the properties of Schwarzschild trajectories. By investigating radial geodesics, we find that no effects of repulsive gravity are present. We interpreted this result as indicating that repulsive gravity is non-linear effect.


Introduction
Although the existence of naked singularities in Nature is the subject of intense debate nowadays, it has been well established that Einstein field equations for the gravitational field allow solutions that can be interpreted as describing naked singularities. In particular, black hole solutions are characterized by the existence of naked singularity counterparts [1]. However, it seems that the particular choice of the physical parameters, which is necessary for the formation of naked singularities, is difficult to be realizable in Nature. Indeed, a rotating naked singularity needs a specific angular momentum that must be greater that its mass (in geometric units), a condition that probably cannot be fulfilled in realistic configurations because it would imply such a high angular velocity that the object would destroy itself before reaching it [2]. a e-mail: saken.yan@yandex.com (corresponding author) b e-mail: quevedo@nucleares.unam.mx c e-mail: abishevme@gmail.com d e-mail: arai12@mail.ru There is, however, a simpler way to generate naked singularities, namely, by considering mass distributions with quadrupole moment [3]. Indeed, from the point of view of multipole moments, the uniqueness theorems prove that black holes can have only mass monopole and angular momentum dipole [4]. Consequently, the addition of a quadrupole to a mass distribution, even in the static case, would imply that the corresponding gravitational field describes a naked singularity. Consequently, a simple shape deviation from spherical symmetry in a mass distribution leads to the appearance of naked singularities. In previous works [3,[5][6][7], we used a particular static quadrupolar solution [8][9][10][11] to study the physical properties of naked singularities.
There are several solutions of Einstein field equations that can be used to describe the exterior gravitational field of a static mass distribution with quadrupole moment [12]. In the limiting case of vanishing quadrupole, they reduce to the spacetime of a Schwarzschild black hole. An interesting characteristic of all of them is that the outermost naked singularity is a sphere with radius r = 2m, where m is the mass of the gravitational source. This could be interpreted intuitively as if the presence of the quadrupole causes the destruction of the regular horizon turning it into a singular hypersurface.
In this work, we show that this is not always the case. Indeed, we will derive a new solution, whose singularity is located on a sphere of radius r = m. This means that in this case the quadrupole completely destroys the regular horizon at r = 2m, but generates a new special hypersurface at r = m which contains a singularity. To derive this new solution, we use the fact in realistic situations we expect that compact objects deviates only slightly from spherical symmetry. This implies that the quadrupole can be considered as a small quantity. With this in mind, we investigate a particular approximate line element, which is valid only up to the first order in the quadrupole parameter. Then, we find the general solution of the corresponding Einstein vacuum field equations and show that a particular solution is characterized by a curvature singularity located on a sphere of radius r = m.
This work is organized as follows. In Sect. 2, we consider a line element that is especially adapted to the study of interior and exterior solutions. We derive the general field equations for the case of vacuum gravitational fields. In Sect. 3, we find the most general solution that is linear in the quadrupole moment. We select a particular case that is characterized by the presence of naked singularity at a distance r = m from the origin of coordinates. We also calculate the Newtonian limit of the new approximate solution and show that it corresponds to a mass distribution with a small quadrupole. In Sect. 4, we investigate the motion of test particles in the spacetime described by the approximate solution. In general, we find that the quadrupole affects the behavior of Schwarzschild orbits. By analyzing the behavior of free falling particles we show that no effects associated with the presence of repulsive gravity can be detected in contrast to repulsive effects found previously in the case of exact solutions with quadrupole. We conclude that repulsive gravity in naked singularities is a non-linear phenomenon. Finally, Sect. 5 contains a summary of our results.

Line element and field equations
The search for and investigation of physically meaningful solutions of Einstein equations begins with the choice of an appropriate line element. In the case of static gravitational fields of deformed mass distributions, one can assume that the fields preserve axial symmetry. Moreover, if we are interested in describing the gravitational field outside as well as inside the mass distribution, it is convenient to choose a line element that can be used in both cases. In a previous work [13], we found interior perfect-fluid solutions that can be matched with exterior vacuum solutions under the assumption that the quadrupole moment is small. From a physical point of view, this implies that the mass distribution is only slightly deformed. We were able to find a particular line element that can be used to search for interior and exterior approximate solutions. It can be written as where the set of (t, r, θ, ϕ) can be interpreted as polar coordinates in the limiting case q → 0. Moreover, the functions ν = ν(r ), a = a(r ), c = c(r ),m =m(r ), and b = b(r, θ) are arbitrary.
In the particular case q = 0, the above line element can be used to describe the exterior Schwarzschild solutioñ and the interior perfect fluid Schwarzschild metric with with where ρ = const and p = p(r ) are the density and pressure of the fluid, respectively. Furthermore, for q = 0 the line element (1) contains the approximate (up to the first order in q) quadrupolar metric (q−metric), which has been interpreted as the simplest generalization of the Schwarzschild metric that includes a quadrupole moment [3]. In this limit, q has been interpreted as the quadrupole parameter.
Here, we will use the advantages of the line element (1) to search for more general approximate solutions. Then, the Einstein vacuum field equations, up to the first order in q, can be written as 2r where a comma represents partial differentiation with respect to the corresponding coordinate. For simplicity, we replaced the solution of the first equationm = m = const. in the remaining equations.

Approximate solutions
We now investigate the system of partial differential equation (6)- (11). Equations (6) and (7) can be integrated and yield where α 1 and α 2 are dimensionless integration constants. It turns out that the remaining system of partial differential equations can be integrated in general and yields where α 3 , α 4 , and α 5 are dimensionless integration constants.
We can see that the general approximate exterior solution with quadrupole moment is represented by the 5-parameter family of solutions (12)- (15). In this general solution, the additive constants α 4 and α 5 can be chosen such that at infinity the solution describes the Minkowski spacetime in spherical coordinates. This means that non asymptotically flat solutions are also contained in the above general solution.
To partially investigate the physical meaning of this solution, we calculate the Kretschmann scalar K = R abcd R abcd . We obtain where where the term proportional to q 2 has been neglected. We can see that this approximate spacetime is characterized by the presence of three different curvature singularities located at The geometric structure of the curvature singularities of the solutions (12), (13), (14) and (15) is illustrated in Fig. 1. Another interesting particular case corresponds to the choice which leads to the line element This is an asymptotically flat approximate solution with parameters m, q and α 2 .
Since α 2 appears always in combination with q, it can be absorbed in the definition of q. We, therefore, set α 2 = 1 without loss of generality. The singularity structure can be found by analyzing the Kretschmann invariant, which in this case reduces to We see that there are only two singularities, which are located at r = 0 and r = m. This is an interesting property of this solution because all the remaining solutions contained in (12)- (15) are singular on the hypersurface r = 2m. Moreover, all the known exact solutions with quadrupole turn out to be singular at r = 2m [12]. To our knowledge, the solution (22) is the only one in which the outermost singularity is located inside the sphere with r = 2m. This means that the spacetime is well defined behind in the interval r ∈ (m, 2m]. We are interested in investigating the properties of this spacetime near the singularity r = m. To further analyze the physical meaning of the solution (22), we calculate the corresponding Newtonian limit. To this end, we perform a coordinate transformation of the form (r, θ) → (R, ϑ) defined by the equations [13,14] and where the β 1 and β 2 are constants and we have neglected terms of the order higher that m 3 /R 3 . Inserting the above coordinates into the metric (22), we obtain the approximate line element where P 2 (cos ϑ) is the Legendre polynomial of degree 2, and we have chosen the free constants as α 2 = 2, β 1 = 1/3, and β 2 = 5/3. We recognize the metric (26) as the Newtonian limit of general relativity, where Φ represents the Newtonian potential. Moreover, the constants can be interpreted as the Newtonian mass and quadrupole moment of the corresponding mass distribution.

Motion of test particles
Consider the trajectory x a (τ ) of a test particle with 4-velocity u a = dx a /dτ =ẋ a . The 4-moment p a = μẋ a of the particle can be normalized so that the equations and constraint for geodesics are given as where = 0, 1, −1 for null, timelike, and spacelike curves, respectively [15]. For the approximate metric (22) we obtain from (30) that geodesics are determined, in general, by the following set of equations The 4-moment p a = μẋ a of the particle can be normalized so that where = 0, 1, −1 for null, timelike, and spacelike curves, respectively [16]. Then, for the approximate metric (22) we obtain from (36) that where is the effective potential and we have used the expression for the energy E = μẼ and the angular moment l = μl of the test particle which are constants of motion associated with the Killing vector fields ξ t = ∂ t and ξ ϕ = ∂ ϕ , respectively. For the sake of simplicity we set μ = 1 so that E = E andl = l. Figure 2 illustrates the behavior of the effective potential in terms of the parameter q for θ = π/2. The effective potential of the Schwarzschild spacetime is also shown for comparison. For positive (negative) values of q, the effective potential at a given point outside the outer singularity is always less (greater) than the Schwarzschild value. This indicates that the distribution of orbits on the equator of the metric (1) can depend drastically on the value of q.

Circular orbits
We will now investigate the properties of circular orbits on the equatorial plane, θ = π/2 of the gravitational described by the approximate metric (1). Circular orbits correspond to the limiting caseṙ = 0. Their stability properties are determined by the extrema of the effective potential. Following the conventional stability analysis of circular orbits involving a potential function, from Eq. (38), we obtain and where B 1 and B 2 are given by where we set = 1 and replaced the value of the angular momentum which can be derived from the condition V e f f,r = 0. The numerical analysis of the stability condition, V e f f,rr > 0, is depicted in Fig. 3. The green region contains only stable orbits whereas the yellow region corresponds to unstable orbits. For comparison, we include the limiting values of the Schwarzschild spacetime. We see that the quadrupole parameter q changes the value of the minimum allowed radius (3m) and of the inner most stable circular orbit radius (6m) of the Schwarzschild metric. In fact, the quadrupole leads to the appearance of a second stable region below and over the radius 3m, which is not present in the Schwarzschild limiting case. This region can reach the value of r/m ≈ 1, approaching the singularity which is located at r = m. Moreover, within the spacetime determined by the interval r ∈ (m, 2m], we notice that most of this region allows the existence of stable circular orbits with a disjoint region of instability for positive values of q.
Furthermore, the energy of test particles on circular orbits can be expressed as In Fig. 4, we plot the regions in which the energy E 2 and angular momentum l 2 are both positive or negative simultaneously. The red region denotes all the radii that are not allowed for circular orbits because either the squared of the energy or of the angular momentum are negative. A comparison with Fig. 3 shows that the region contained between the singularity r/m = 1 and around r/m = 3 is allowed for circular orbits by the stability condition but is excluded by the energy and angular momentum conditions. We conclude that the effect of the quadrupole on the properties of circular orbits is as follows. A positive quadrupole Fig. 4 Energy E 2 and angular momentum l 2 of test particles on circular orbits in terms of the radius r/m and the quadrupole parameter q leads to an increase of the minimum allowed radius whereas a negative quadrupole generates the opposite effect. This means that only in the case of an oblate object, test particles are allowed to exist on orbits closer to the singularity, which is situated at r = m.
Finally, we consider the angular velocity and the period of circular orbits. The behavior of the angular velocity and period are depicted in Fig. 5. We can see that the influence of the quadrupole on the value of the angular velocity increases as the radius of the orbit approaches the value of r ≈ 3m. This agrees with the behavior of the stability condition and the energy and angular momentum shown in Figs. 3 and 4, respectively. For a given orbit radius, the angular velocity increases (decreases) for positive (negative) values of the quadrupole. Notice that close to the singularity located at r = m, there is a region in which the angular velocity is a well behaved function of q and r . This region corresponds to the stability region that was also found in Fig. 3.

Bounded and unbounded orbits
We now study the influence of the quadrupole on the trajectories of massive test particles, moving along unbounded paths We consider first unbounded Schwarzschild trajectories with non-zero initial radial velocities under the influence of the quadrupole. In this case, we see that for the chosen initial angular and radial velocities all the particles escape from the gravitational field of central slightly deformed body. This is illustrated Fig. 6. The direction along which the particle escapes to infinity depends on the value of the quadrupole. It is worth noticing that, in principle, this effect could be used to measure the quadrupole of the central mass distribution. In Fig. 7, we consider a Schwarzschild bounded orbit with zero initial radial velocity and the same non-zero value for the initial angular velocity. The left panel shows the Schwarzschild geodesic. The central and right panels illustrate the influence of a negative and positive small quadrupole, respectively. We conclude that the small quadrupole does not affect the bounded character of the geodesic, but it does drastically modify the morphology of the trajectories.
We conclude that the quadrupole always affects the Schwarzschild trajectories. The explicit modifications depend on the properties of the original Schwarzschild trajectory and the value of the quadrupole parameter.

Radial geodesics and repulsive gravity
We now study the free fall of test particles. To this end, we consider the geodesics equations in their most general form as given in Eqs. (32)-(35). The important point is that in all the cases to be considered all the initial spatial velocities are assumed to vanish. The starting point can be chosen arbitrarily, but some special values of the angle θ are of interest, namely, the symmetry axis θ = 0, the equatorial plane θ = π/2, and some other different value that we can choose as θ = π/4. In fact, in [6], it was shown that by analyzing the behavior of radial geodesics one can detect the presence of repulsive gravity [17][18][19].
A free falling particle will continue its motion along the radial direction, unless a force acts on it and changes the original radial direction. This phenomenon has been reported in the case of an exact solution with quadrupole moment in [6].
We will now consider the same situation in the case of the approximate solution we are analyzing in this work. The result of the integration of the geodesic equations (32)-(35) is illustrated in Fig. 8. We see that, in fact, free falling particles move along the original directions (θ = 0, θ = π/4, and θ = π/2), independently of the value of quadrupole parameter q. This means that no repulsive gravity effects can be detected in the case of the approximate metric. Taking into account also the result of [6], we conclude that repulsive effects are non-linear, i.e., they appear only in the case of an exact quadrupolar metric.

Conclusions and remarks
In this work, we have derived a family of approximate vacuum solutions of Einstein equations with quadrupole. Among all the solutions contained in this family we choose a particular one, which presents a naked singularity at the hypersurface r = m, instead of r = 2m as in other quadrupolar metrics. To our knowledge this is the only metric with such a singularity.
By applying an appropriate coordinate transformation, we found the Newtonian limit of the approximate solution and showed that it corresponds to the gravitational potential of a mass with quadrupole.
Then, we investigated the motion of test particles along circular orbits. We established that a positive quadrupole leads to an increase of the minimum allowed radius whereas a nega- Fig. 8 Free fall of test particles with vanishing initial velocities. On the axes θ = 0 (top, left), θ = π/4 (top, right), and θ = π/2 (bottom). The sphere represents the singularity at the radius r = m (top). The trajectories of the free falling test particles do not depend on the value of the parameter q. For concreteness, we set the energy of the test particle as E = 1 tive quadrupole generates the opposite effect. This means that only in the case of an oblate object, test particles are allowed to exist on orbits closer to the singularity, which is situated at r = m. In the case of bounded and unbounded trajectories, we found that the quadrupole always affects the Schwarzschild trajectories. The explicit modifications depend on the properties of the original Schwarzschild trajectory and the value of the quadrupole parameter.
Finally, we analyzed radial geodesics that have been used previously to detect effects of repulsive gravity in an exact quadrupolar metric. However, in the case of the approximate solution no repulsive effects were found. We conclude that repulsive gravity in the quadrupolar naked singularities is a non-linear phenomenon.