A generalization of q-metric

Shokoufe Faraji∗ University of Bremen, Center of Applied Space Technology and Microgravity (ZARM), 28359 Germany

In this paper, Weyl’s procedure is used to investigate an alternative class of the q-metric that describes a deformed compact object in the presence of an external distribution of matter via exercising quadrupole moments and taking advantage of combining different approaches.

I. INTRODUCTION space-time, the horizon remains regular [7]. However, a detailed analysis of the distorted Schwarzschild space- In finding the exact and approximate time’s global properties was introduced in 1982 by Ge- solutions describing the real source is always of high in- roch & Hartle [8]. Later the explicit form of this metric terest. There exists quite a large number of the solutions was presented in [9]. of the Einstein field equations in the literature. However, In this paper, the static q-metric has been general- when the solution supposed to describe a deviation from ized by considering an external matter distribution up the spherical symmetry, this is a challenging task indeed. to quadrupoles following the approach of Weyl and Eres- Furthermore, within the scope of general relativity, Rosen. Thus, the result is not asymptotically flat by its Birkoff’s theorem states that Schwarzschild space-time is construction. Mathematically, it is possible to find vari- the unique static solution of the asymptotically flat vac- ous vacuum solutions with mass and different quadrupole uum Einstein field equation with a regular . parameters. Nevertheless, differences appear only at the However, by considering a static mass distribution with higher multipoles [10]. This means that all the solutions a quadrupole moment, characterizing the deviation from can describe the exterior gravitational field of distorted spherical symmetry, this theorem will no longer be valid. mass distribution up to quadrupole moments, since they Therefore it is possible to find the various vacuum solu- satisfy all the necessary conditions to describe the exte- tions with mass and different quadrupole parameters. In rior gravitational field of realistic compact objects. In this respect, the first static and axially symmetric solu- this paper, we also focus on exterior solutions. The first tions with arbitrary quadrupole moment are described by reason for choosing q-metric to work with is its simplic- Weyl in 1917 [1]. Then Erez and Rosen discovered static ity. In fact, its mathematical structure is straightforward, solutions with arbitrary quadrupole in prolate spheroidal which facilitates the study and can be handled analyti- coordinates in 1959 [2]. Later Zipoy [3], and Voorhees [4] cally. From the physics perspective, this solution is sim- found a transformation that generates the simplest solu- ilar to the q-metric solution with an additional external tion with quadrupole. This metric is known as γ-metric gravitational field, similar to the addition of a magnetic or σ-metric, and possesses interesting physical aspects environment to the solution. with strange topologies. However, for exact spherical Furthermore, aside from the simplicity of q-metric, symmetry it reduces to the . Later there are several motives to consider this generalization. on, by representing this metric in terms of parameter q, For instance, in the relativistic astrophysical study, it it is known as q-metric [5]. is assumed that astrophysical compact objects are de- The above procedure, however, leads to having scribed by the Schwarzschild or Kerr space-times. How- Minkowski space as the limiting case. Indeed, this as- ever, besides these relevant solutions, others can imitate sumption is equivalent to restrict the solution to describ- a black hole’s properties, such as the electromagnetic sig- ing the isolated cases in the space. While this assumption nature. In this respect, this work aims at investigating is natural in the first place, the question of how the ex- how different backgrounds may describe an astrophysi- arXiv:2010.15723v4 [gr-qc] 19 Jul 2021 ternal distribution of mass may distort them might be of cal phenomenon with departure from spherical geometry. some interest. In this perspective, the second class of so- Besides, one can take quadrupole moments as the addi- lutions describing non-asymptotically flat space-time, is tional physical degrees of freedom to the central compact obtained by assuming the existence of a static and axially object and its surroundings. This can facilitate search- symmetric external distribution of matter in the vicinity ing for the link between observational phenomena and the of the central source. In fact, this space-time is a local compact object’s properties which is not isolated, since solution to the Einstein field equation by its construc- the identification of a compact object from initial obser- tion [6]. In 1965, Doroshkevich and his colleagues con- vations is a challenging task. For example, the geometric sidered an external gravitational field up to a quadrupole configuration of an accretion disk located in this space- in the Schwarzschild space-time. They also showed that time depends explicitly on the value of both quadrupole by adding quadrupole correction to the Schwarzschild parameters in a way that it is always possible to distin- guish between a distorted Schwarzschild black hole and a distorted, deformed compact object. For example, test- ing quasi-periodic oscillations in this background, which ∗ [email protected] is progressing, shows its ability to connect the data to the 2 model. Moreover, there is no doubt on the fundamental importance of gravitational waves in physics, where the  (1+α)  (1+α) experimental evidence finally supported the purely theo- x − 1 x + 1 ds2 = − dt2 + M 2(x2 − 1) retical research in this area. In fact, this metric can also x + 1 x − 1 apply to the study of gravitational waves generated by "  x2 − 1 α(2+α)  dx2 dy2  a slightly non-isolated and self-gravitating axisymmetric + distribution of mass. In fact, Ryan has shown [11] that x2 − y2 x2 − 1 1 − y2 we can extract the multipole moments of the central body +(1 − y2)dφ2 , (2) from the gravitational wave signal. Thus, one can par- ticularly apply this metric in the colliding gravitational with the following transformation law waves where asymptotic flatness is not a requirement. As another example, the gravitational field of a very large r plate of thin matter, where its gravity mainly focuses x = − 1 , y = cos θ . (3) at a finite distance, and therefore it stops behaving well M asymptotically [12]. This metric has the central curvature singularity at x = As the next step of this work, one can add rotation −1, where it occurs for all real values of α. Moreover, an into this setup to describe more realistic scenarios. How- additional singularity appears at x = 1, and the norm of ever, this is worth mentioning that it is feasible to explain the time-like Killing vector at this radius vanishes. Also, some of the observational data with a static setup. For outside this hypersurface, there exists no additional hori- example, it has been shown that the possible resonant zon. However, considering relatively small quadrupole oscillations of the thick accretion disk, can be observed moment, this hypersurface is located very close to the ori- even if the source of radiation is steady and perfectly gin of coordinates so that a physically reasonable interior axisymmetric. solution could be used to cover them [13]. Besides, out The paper’s organization is as follows: q-metric is of this region there is no more singularity, and the metric briefly explained in SectionII. Considering external is asymptotically flat. Although there are different defi- source is presented in SectionIII. The conclusion is sum- nitions of multipole moments, at the lower order all are marised in SectionV. equivalent and link to parameter α, where α ∈ (−1, ∞). In this paper, the geometrized units where c = 1 and Further, all odd multipole moments vanish because of the G = 1, is used. reflection symmetry with respect to the equatorial plane.

III. GENERALIZATION OF q-METRIC II. q-METRIC In this section, we explain briefly the Weyl procedure. The q-metric describes a static, axially symmetric, and The first static and axisymmetric solutions with arbitrary asymptotically flat solution to the Einstein equation with quadrupole moment are described via the Weyl metric arbitrary quadrupole moment. The metric represents the exterior gravitational field of an isolated static axisym- metric mass distribution. It can be used to investigate ds2 = −e2ψ0 dt2 + e2(γ0−ψ0)(dρ2 + dz2) + e−2ψ0 ρ2dφ2 , the exterior fields of slightly deformed astrophysical ob- (4) jects in the strong-field regime [13]. In fact, the presence of quadrupole, independent of its value, changes the ge- where ψ0 = ψ0(ρ, z) and γ0 = γ0(ρ, z) are the metric ometric properties of space-time drastically. The metric functions, and ψ0 plays the role of gravitational poten- is presented as follows tial. The relation between Schwarzschild coordinates and Weyl coordinates is given by

 2M 1+α  2M −α ds2 = 1 − dt2 − 1 − p r r ρ = r(r − 2M) sin θ, z = (r − M) cos θ, (5) " −α(2+α) !  M 2 sin2 θ  dr2 where M is a parameter that is considered as the mass 1 + + r2dθ2 r2 − 2Mr 1 − 2M of the compact object. If we consider three-dimension r manifold N, orthogonal to the static Killing vector field, 2 2 2 +r sin θdφ . (1) then this metric induces the flat metric on N. Besides, the metric function ψ0 with respect to this flat metric However, in 1959, Erez and Rosen pointed out that the obeys the Laplace equation multipole structure of a static axially symmetric solution has a simpler form in the prolate spheroidal coordinates [14] rather than the cylindrical coordinates of Weyl [2, 1 ψ0,ρρ + ψ0,ρ + ψ0,zz = 0, (6) 10, 15] ρ 3 which is the key factor in this approach. Also, the field To sum up, all the process reduces to choosing two fields γ0 is determined by the following relation [1] that provide a solution of the vacuum field equations and take them as the seed for introducing new fields that sat- isfy the above relations. In this approach, metric compo- ∇γ ∇ρ = ρ(∇ψ )2, (7) 0 0 nents are expressed through integrability conditions as where ∇ := ∂ρ + i∂z, or functions of a single harmonic function and its partial derivatives in a surprisingly simple mathematical proce-

2 2 dure. γ0,ρ =ρ(ψ0,ρ − ψ0,z), Following Weyl’s procedure by using the q-metric fields γ0,z =2ρψ0,ρψ0,z. (8) as the seed, we can obtain the solution of a deformed source surrounding with the external mass distribution, Now, Weyl’s technique takes advantage of the linearity which can be generalized to the whole class of q-metric. of Laplace’s equation for ψ0 and modifies our original Thus, we introduce auxiliary fields ψ and γ such that metric by replacing the field ψ0 by

 (1+α) ˆ x − 1 ˆ ψ = ψ0 + ψ. (9) e2ψ := e2ψ, (18) x + 1 Since equation (7) is not linear with respect to γ0, there-  x2 − 1 α(1+α) fore with the same modification as ψ, it is not possible to e2γ := e2ˆγ . (19) find a new fieldγ ˆ. In fact, it turns out that there will also x2 − y2 be a contribution of another term, say δγ, that satisfies Notice that the exponential factors, distorted terms, rep- ∇(δγ) = C(ρ)(∇ψ0), (10) resent the surrounding matter. As this process suggests, ˆ where C is a function of ρ [1,2]. This relation deter- we preserve the q-metric fields by takingγ ˆ = ψ = 0. It mines δγ up to some constant. However, the requirement can easily be checked that in the limits ψˆ = 0,γ ˆ = 0, of elementary flatness [16], in the neighbourhood of the and α = 0 we recover the Schwarzschild fields, and in the symmetry axis fixes the constant, and it should be set case α = 0, butγ ˆ 6= 0, ψˆ 6= 0 we obtain the distorted equal to zero. In what follows, we have chosen them via Schwarzschild fields written in the prolate spheroidal co- a standard procedure employing Legendre expansion. ordinates. However, following the Erez-Rosen approach, it is more In this way, field functions preserve their properties as appropriate to rewrite this metric in prolate spheroidal before by its construction. Ultimately, the metric is then coordinates (t, x, y, φ), via this transformation relation given by

1 p 2 2 p 2 2 (1+α) x = ( ρ + (z + M) + ρ + (z − M) ) , (11) x − 1 ˆ ˆ 2M ds2 = − e2ψdt2 + M 2(x2 − 1)e−2ψ x + 1 1 p p y = ( ρ2 + (z + M)2 − ρ2 + (z − M)2). (12) " 2M x + 1(1+α)  x2 − 1 α(2+α) e2ˆγ Equivalently x − 1 x2 − y2  dx2 dy2   p + + (1 − y2)dφ2 , (20) ρ = M (x2 − 1)(1 − y2), (13) x2 − 1 1 − y2 z = Mxy. (14) where t ∈ (−∞, +∞), x ∈ (1, +∞), y ∈ [−1, 1], and Also, the equations (6) and (8) in the prolate spheroidal φ ∈ [0, 2π). In this work for choosing the explicit form coordinates are written as of these distortion functions we followed the standard procedure via the separation of variables method [17]. In 2 y this approach ψˆ generally can be expressed in terms of (x − 1)ψ0,xx + 2xψ0,x + ψ0,yy + ψ0,y = 0, p1 − y2 Legendre polynomials of the first kind and second kind (15) [9, equation 6-7]. This method was introduced first by y2 − 1 Geroch and Hartle [8]. However, up to the quadrupole γ = [x ψ2 − (x2 − 1)ψ2  + β, these are expressed as follows [9][18], 0,x x2 − y2 0,y 0,x 2y(1 − x2) ψ0,xψ0,y], (16) β p1 − y2 ψˆ = − −3x2y2 + x2 + y2 − 1 , (21) 2 2 1 − y y 2 2 2  2 γ0,y = [ ψ − (x − 1)ψ + γˆ = −2xβ(1 − y ) x2 − y2 p 2 0,y 0,x 1 − y β2 + (x2 − 1)(1 − y2)(−9x2y2 + x2 + y2 − 1). (22) 2xψ0,yψ0,x]. (17) 4 4

We may refer to quadrupole α as a deformation param- initial values and parameters, particles can fall onto the eter and β as a distortion parameter to be compatible central object, escape from it, or trap in some region and with the proceeding related works. forming a toroidal shape around the central object. As Therefore, this metric contains three free parameters: one can see in the plots, only by a slight change in the the total mass, deformation parameter α, and distortion initial values one can obtain different trajectories for a β, which are taken to be relatively small and connected to chosen pair of parameters. the compact object deformation and external mass pres- ence distribution, respectively. The result is also locally valid by its construction [6,8], and may be considered as the local q-metric or distorted q-metric. From the physics perspective, this solution is similar to the q-metric solu- V. SUMMARY AND CONCLUSION tion with an additional external gravitational field, like adding a magnetic surrounding [19–21]. This paper has presented a class of q-metric for the rel- atively small quadrupole moment via Weyl’s procedure. IV. EFFECTIVE POTENTIAL This class explains the exterior of a deformed body lo- cally in the presence of an external distribution of matter Regarding the symmetries in the metric, there are two up to the quadrupole. It contains three free parame- constants of geodesic motion E and L, ters: the total mass, deformation parameter α, and the distortion parameter β, referring to the central object’s quadrupoles and its surrounding mass distribution.  (1+α) x − 1 2ψˆ E = −gttt˙ = e t,˙ Exercising the quadrupole moments may open the door x + 1 to a broader range of applications via considering them as  (1+α) new degrees of freedom to both central object and its sur- ˙ 2 2 x + 1 −2ψˆ ˙ L = gφφφ = M (x − 1) e φ. (23) rounding, to search for their observational fingerprints. x − 1 For instance, in the study of the gravitational wave, or By using these relations and the normalization condition quasi-periodic oscillations generated by a non-isolated, ρ ν gρν x˙ x˙ = −, where  can take values −1, 0 and 1, for self-gravitating axisymmetric distribution of mass. the space-like, light-like and for the time-like trajectories, We propose that this class has promising features, and respectively, then the geodesic equation is reduced to it may serve to link the metric to the physical nature of the compact object via its parameters. Although its simplicity still this generalization worth considering to  x2 − 1 α(2+α) 1 + M 2e2ˆγ x˙ 2 + V 2 = E2, (24) have more vacuum metric available, for example, to study x2 − y2 a vacuum solution that is accelerating or embedded in an external field. Therefore, it may be of some interest to where investigate how the external field affects the geometry and geodesics of the q-metric. Of course, real models x2 − 1  x2 − 1 α(2+α) will be inaccurate due to incomplete knowledge of the V 2 = M 2e2ˆγ y˙2 underlying physics, or perhaps a desire for simplicity. 1 − y2 x2 − y2 In the present contribution we further explored the ex- L2e4ψˆ x − 1(α+1) 2ψˆ istence of different trajectories for a particle in the equa- + 2 2 2 + e . (25) M (1 − y )(x + 1) x + 1 torial plane for correct choices of parameters. In general, One can interpret this equation as the motion along the test particles motion can be chaotic in this background x coordinate in terms of this potential, which is called for some combinations of parameters α and β. Also, colli- sion may occur in this vicinity of central object. Further, the effective potential Veff . However, this Furthermore, by analysing the effective potential in the the model can reproduce many of the relevant features equatorial plane, one can investigate general properties of the numerical results. At this stage, such qualitative of the particle dynamics in the equatorial plane. In fact, agreement is all that can be hoped for. this motion is determined by the energy’s boundaries The next step of this work could be a study on off- 2 given by E = Veff . In general, possible types of orbits, equatorial time-like and light-like geodesics. Studying dependent on the parameters , E, L, α and β. And the topological implication of this background may be in an analytical exploration depending on the number of the next stage of study in this area. Also, the con- positive real zeros and the sign of E2 − , one obtains struction of accretion disks, the study of gravitational different types of trajectories. In Figure1 the trajecto- waves and quasi-periodic oscillation in this background ries for some choices of the parameters have plotted. In are in progress. We expect future work in this field to fact, for each valid pair of parameters, one can distinguish be guided by questions of astrophysics, and other areas different trajectories in this background. Depending on where strong-field gravitational theory applies. 5

FIG. 1. Timelike geodesic for pairs of (α = 0.1, β = 0.000001). The trajectories in the (r, φ) section and in the complete 3D are plotted. Configurations in the first column correspond to a bounded trajectory, in the second column to a falling trajectory, and in the third column to an escaping trajectory.

VI. ACKNOWLEDGEMENTS mann for their valuable comments on this work. Also, thanks to Dr. Audry Trova for valuable discussions. This work is supported by the research training group GRK The author gratefully acknowledges Prof. Hernando 1620 ”Models of Gravity”, funded by the German Re- Quevedo, Prof. Domenico Giulini, and Dr. Eva Hack- search Foundation (DFG).

[1] . Zur gravitationstheorie. Annalen der [9] Nora Bret´on, Tatiana E. Denisova, and Vladimir S. Physik, 359(18):117–145, 1917. Manko. A Kerr black hole in the external gravitational [2] G. Erez and N. Rosen. The gravitational field of a particle field. Physics Letters A, 230:7–11, Feb 1997. possessing a multipole moment. Bull. Research Council [10] Hernando Quevedo. General static axisymmetric solution Israel, Vol: Sect. F.8, 9 1959. of einstein’s vacuum field equations in prolate spheroidal [3] David M. Zipoy. Topology of some spheroidal metrics. coordinates. Phys. Rev. D, 39:2904–2911, May 1989. Journal of Mathematical Physics, 7(6):1137–1143, 1966. [11] Fintan D. Ryan. Gravitational waves from the inspiral [4] B. H. Voorhees. Static axially symmetric gravitational of a compact object into a massive, axisymmetric body fields. Phys. Rev. D, 2:2119–2122, Nov 1970. with arbitrary multipole moments. PRD, 52(10):5707– [5] Hernando Quevedo. Mass Quadrupole as a Source of 5718, November 1995. Naked Singularities. International Journal of Modern [12] J. Horsk´y.A comment of thin shells in general relativity. Physics D, 20(10):1779–1787, January 2011. Physics Letters A, 28(9):599–599, February 1969. [6] S Chandrasekhar. The mathematical theory of black [13] Hernando Quevedo. Exterior and interior metrics with holes. Oxford classic texts in the physical sciences. Ox- quadrupole moment. Gen. Rel. Grav., 43:1141–1152, ford Univ. Press, Oxford, 2002. 2011. [7] A. G. Doroshkevich, Y. B. Zel’dovich, and I. D. Novikov. [14] Prolate spheroidal coordinates are three-dimensional or- Gravitational Collapse of Non-Symmetric and Rotating thogonal coordinates that result from rotating the two- Bodies. Zhurnal Eksperimentalnoi i Teoreticheskoi Fiziki, dimensional elliptic coordinates about the focal axis of 4:170, December 1965. the ellipse. [8] R. Geroch and J. B. Hartle. Distorted black holes. Jour- [15] Hernando Quevedo. Mass quadrupole as a source of nal of Mathematical Physics, 23:680–692, 1982. naked singularities. International Journal of Modern 6

Physics D, 20(10):1779–1787, 2011. [18] These functions also have been considered for the [16] In general, both functions γ and ψ should be regular at distorted Schwarzschild solution. In this case, the the symmetry axis y = ±1 [9]. Sometimes this condition Schwarzschild metric considers as a seed metric via this is referred to as the elementary flatness condition. approach [6,8,9]. [17] Milton Abramowitz. Handbook of Mathematical Func- [19] F. J. Ernst. Black holes in a magnetic universe. Journal tions, With Formulas, Graphs, and Mathematical Ta- of Mathematical Physics, 17(1):54–56, January 1976. bles,. Dover Publications, Inc., New York, NY, USA, [20] Frederick J. Ernst and Walter J. Wild. Kerr black holes 1974. in a magnetic universe. Journal of Mathematical Physics, 17(2):182–184, February 1976. [21] Robert M. Wald. Black hole in a uniform magnetic field. Phys. Rev. D, 10:1680–1685, Sep 1974.