In general relativity, the Weyl metrics (named after the German-American mathematician Hermann Weyl)^[1] are a class of static and axisymmetric solutions to Einstein's field equation. Three members in the renowned Kerr–Newman family solutions, namely the Schwarzschild, nonextremal Reissner–Nordström and extremal Reissner–Nordström metrics, can be identified as Weyl-type metrics.

Standard Weyl metrics

The Weyl class of solutions has the generic form^[2][3]

$${\displaystyle ds^{2}=-e^{2\psi (\rho ,z)}dt^{2}+e^{2\gamma (\rho ,z)-2\psi (\rho ,z)}(d\rho ^{2}+dz^{2})+e^{-2\psi (\rho ,z)}\rho ^{2}d\phi ^{2}\,,}$$

(1)

where ${\displaystyle \psi (\rho ,z)}$ and ${\displaystyle \gamma (\rho ,z)}$ are two metric potentials dependent on Weyl's canonical coordinates ${\displaystyle \{\rho \,,z\))$. The coordinate system ${\displaystyle \{t,\rho ,z,\phi \))$ serves best for symmetries of Weyl's spacetime (with two Killing vector fields being ${\displaystyle \xi ^{t}=\partial _{t))$ and ${\displaystyle \xi ^{\phi }=\partial _{\phi ))$) and often acts like cylindrical coordinates,^[2] but is incomplete when describing a black hole as ${\displaystyle \{\rho \,,z\))$ only cover the horizon and its exteriors.

Hence, to determine a static axisymmetric solution corresponding to a specific stress–energy tensor ${\displaystyle T_{ab))$, we just need to substitute the Weyl metric Eq(1) into Einstein's equation (with c=G=1):

$${\displaystyle R_{ab}-{\frac {1}{2))Rg_{ab}=8\pi T_{ab}\,,}$$

(2)

and work out the two functions ${\displaystyle \psi (\rho ,z)}$ and ${\displaystyle \gamma (\rho ,z)}$.

Reduced field equations for electrovac Weyl solutions

One of the best investigated and most useful Weyl solutions is the electrovac case, where ${\displaystyle T_{ab))$ comes from the existence of (Weyl-type) electromagnetic field (without matter and current flows). As we know, given the electromagnetic four-potential ${\displaystyle A_{a))$, the anti-symmetric electromagnetic field ${\displaystyle F_{ab))$ and the trace-free stress–energy tensor ${\displaystyle T_{ab))$ ${\displaystyle (T=g^{ab}T_{ab}=0)}$ will be respectively determined by

$${\displaystyle F_{ab}=A_{b\,;\,a}-A_{a\,;\,b}=A_{b\,,\,a}-A_{a\,,\,b))$$

(3)

$${\displaystyle T_{ab}={\frac {1}{4\pi ))\,\left(\,F_{ac}F_{b}^{\;c}-{\frac {1}{4))g_{ab}F_{cd}F^{cd}\right)\,,}$$

(4)

which respects the source-free covariant Maxwell equations:

$${\displaystyle {\big (}F^{ab}{\big )}_{;\,b}=0\,,\quad F_{[ab\,;\,c]}=0\,.}$$

(5.a)

Eq(5.a) can be simplified to:

$${\displaystyle \left({\sqrt {-g))\,F^{ab}\right)_{,\,b}=0\,,\quad F_{[ab\,,\,c]}=0}$$

(5.b)

in the calculations as ${\displaystyle \Gamma _{bc}^{a}=\Gamma _{cb}^{a))$. Also, since ${\displaystyle R=-8\pi T=0}$ for electrovacuum, Eq(2) reduces to

$${\displaystyle R_{ab}=8\pi T_{ab}\,.}$$

(6)

Now, suppose the Weyl-type axisymmetric electrostatic potential is ${\displaystyle A_{a}=\Phi (\rho ,z)[dt]_{a))$ (the component ${\displaystyle \Phi }$ is actually the electromagnetic scalar potential), and together with the Weyl metric Eq(1), Eqs(3)(4)(5)(6) imply that

$${\displaystyle \nabla ^{2}\psi =\,(\nabla \psi )^{2}+\gamma _{,\,\rho \rho }+\gamma _{,\,zz))$$

(7.a)

$${\displaystyle \nabla ^{2}\psi =\,e^{-2\psi }(\nabla \Phi )^{2))$$

(7.b)

$${\displaystyle {\frac {1}{\rho ))\,\gamma _{,\,\rho }=\,\psi _{,\,\rho }^{2}-\psi _{,\,z}^{2}-e^{-2\psi }{\big (}\Phi _{,\,\rho }^{2}-\Phi _{,\,z}^{2}{\big )))$$

(7.c)

$${\displaystyle {\frac {1}{\rho ))\,\gamma _{,\,z}=\,2\psi _{,\,\rho }\psi _{,\,z}-2e^{-2\psi }\Phi _{,\,\rho }\Phi _{,\,z))$$

(7.d)

$${\displaystyle \nabla ^{2}\Phi =\,2\nabla \psi \nabla \Phi \,,}$$

(7.e)

where ${\displaystyle R=0}$ yields Eq(7.a), ${\displaystyle R_{tt}=8\pi T_{tt))$ or ${\displaystyle R_{\varphi \varphi }=8\pi T_{\varphi \varphi ))$ yields Eq(7.b), ${\displaystyle R_{\rho \rho }=8\pi T_{\rho \rho ))$ or ${\displaystyle R_{zz}=8\pi T_{zz))$ yields Eq(7.c), ${\displaystyle R_{\rho z}=8\pi T_{\rho z))$ yields Eq(7.d), and Eq(5.b) yields Eq(7.e). Here ${\displaystyle \nabla ^{2}=\partial _{\rho \rho }+{\frac {1}{\rho ))\,\partial _{\rho }+\partial _{zz))$ and ${\displaystyle \nabla =\partial _{\rho }\,{\hat {e))_{\rho }+\partial _{z}\,{\hat {e))_{z))$ are respectively the Laplace and gradient operators. Moreover, if we suppose ${\displaystyle \psi =\psi (\Phi )}$ in the sense of matter-geometry interplay and assume asymptotic flatness, we will find that Eqs(7.a-e) implies a characteristic relation that

$${\displaystyle e^{\psi }=\,\Phi ^{2}-2C\Phi +1\,.}$$

(7.f)

Specifically in the simplest vacuum case with ${\displaystyle \Phi =0}$ and ${\displaystyle T_{ab}=0}$, Eqs(7.a-7.e) reduce to^[4]

$${\displaystyle \gamma _{,\,\rho \rho }+\gamma _{,\,zz}=-(\nabla \psi )^{2))$$

(8.a)

$${\displaystyle \nabla ^{2}\psi =0}$$

(8.b)

$${\displaystyle \gamma _{,\,\rho }=\rho \,{\Big (}\psi _{,\,\rho }^{2}-\psi _{,\,z}^{2}{\Big )))$$

(8.c)

$${\displaystyle \gamma _{,\,z}=2\,\rho \,\psi _{,\,\rho }\psi _{,\,z}\,.}$$

(8.d)

We can firstly obtain ${\displaystyle \psi (\rho ,z)}$ by solving Eq(8.b), and then integrate Eq(8.c) and Eq(8.d) for ${\displaystyle \gamma (\rho ,z)}$. Practically, Eq(8.a) arising from ${\displaystyle R=0}$ just works as a consistency relation or integrability condition.

Unlike the nonlinear Poisson's equation Eq(7.b), Eq(8.b) is the linear Laplace equation; that is to say, superposition of given vacuum solutions to Eq(8.b) is still a solution. This fact has a widely application, such as to analytically distort a Schwarzschild black hole.

We employed the axisymmetric Laplace and gradient operators to write Eqs(7.a-7.e) and Eqs(8.a-8.d) in a compact way, which is very useful in the derivation of the characteristic relation Eq(7.f). In the literature, Eqs(7.a-7.e) and Eqs(8.a-8.d) are often written in the following forms as well:

$${\displaystyle \psi _{,\,\rho \rho }+{\frac {1}{\rho ))\psi _{,\,\rho }+\psi _{,\,zz}=\,(\psi _{,\,\rho })^{2}+(\psi _{,\,z})^{2}+\gamma _{,\,\rho \rho }+\gamma _{,\,zz))$$

(A.1.a)

$${\displaystyle \psi _{,\,\rho \rho }+{\frac {1}{\rho ))\psi _{,\,\rho }+\psi _{,\,zz}=e^{-2\psi }{\big (}\Phi _{,\,\rho }^{2}+\Phi _{,\,z}^{2}{\big )))$$

(A.1.b)

$${\displaystyle {\frac {1}{\rho ))\,\gamma _{,\,\rho }=\,\psi _{,\,\rho }^{2}-\psi _{,\,z}^{2}-e^{-2\psi }{\big (}\Phi _{,\,\rho }^{2}-\Phi _{,\,z}^{2}{\big )))$$

(A.1.c)

$${\displaystyle {\frac {1}{\rho ))\,\gamma _{,\,z}=\,2\psi _{,\,\rho }\psi _{,\,z}-2e^{-2\psi }\Phi _{,\,\rho }\Phi _{,\,z))$$

(A.1.d)

$${\displaystyle \Phi _{,\,\rho \rho }+{\frac {1}{\rho ))\Phi _{,\,\rho }+\Phi _{,\,zz}=\,2\psi _{,\,\rho }\Phi _{,\,\rho }+2\psi _{,\,z}\Phi _{,\,z))$$

(A.1.e)

and

$${\displaystyle (\psi _{,\,\rho })^{2}+(\psi _{,\,z})^{2}=-\gamma _{,\,\rho \rho }-\gamma _{,\,zz))$$

(A.2.a)

$${\displaystyle \psi _{,\,\rho \rho }+{\frac {1}{\rho ))\psi _{,\,\rho }+\psi _{,\,zz}=0}$$

(A.2.b)

$${\displaystyle \gamma _{,\,\rho }=\rho \,{\Big (}\psi _{,\,\rho }^{2}-\psi _{,\,z}^{2}{\Big )))$$

(A.2.c)

$${\displaystyle \gamma _{,\,z}=2\,\rho \,\psi _{,\,\rho }\psi _{,\,z}\,.}$$

(A.2.d)

Considering the interplay between spacetime geometry and energy-matter distributions, it is natural to assume that in Eqs(7.a-7.e) the metric function ${\displaystyle \psi (\rho ,z)}$ relates with the electrostatic scalar potential ${\displaystyle \Phi (\rho ,z)}$ via a function ${\displaystyle \psi =\psi (\Phi )}$ (which means geometry depends on energy), and it follows that

$${\displaystyle \psi _{,\,i}=\psi _{,\,\Phi }\cdot \Phi _{,\,i}\quad ,\quad \nabla \psi =\psi _{,\,\Phi }\cdot \nabla \Phi \quad ,\quad \nabla ^{2}\psi =\psi _{,\,\Phi }\cdot \nabla ^{2}\Phi +\psi _{,\,\Phi \Phi }\cdot (\nabla \Phi )^{2},}$$

(B.1)

Eq(B.1) immediately turns Eq(7.b) and Eq(7.e) respectively into

$${\displaystyle \Psi _{,\,\Phi }\cdot \nabla ^{2}\Phi \,=\,{\big (}e^{-2\psi }-\psi _{,\,\Phi \Phi }{\big )}\cdot (\nabla \Phi )^{2},}$$

(B.2)

$${\displaystyle \nabla ^{2}\Phi \,=\,2\psi _{,\,\Phi }\cdot (\nabla \Phi )^{2},}$$

(B.3)

which give rise to

$${\displaystyle \psi _{,\,\Phi \Phi }+2\,{\big (}\psi _{,\,\Phi }{\big )}^{2}-e^{-2\psi }=0.}$$

(B.4)

Now replace the variable ${\displaystyle \psi }$ by ${\displaystyle \zeta :=e^{2\psi ))$, and Eq(B.4) is simplified to

$${\displaystyle \zeta _{,\,\Phi \Phi }-2=0.}$$

(B.5)

Direct quadrature of Eq(B.5) yields ${\displaystyle \zeta =e^{2\psi }=\Phi ^{2}+{\tilde {C))\Phi +B}$, with ${\displaystyle \{B,{\tilde {C))\))$ being integral constants. To resume asymptotic flatness at spatial infinity, we need ${\displaystyle \lim _{\rho ,z\to \infty }\Phi =0}$ and ${\displaystyle \lim _{\rho ,z\to \infty }e^{2\psi }=1}$, so there should be ${\displaystyle B=1}$. Also, rewrite the constant ${\displaystyle {\tilde {C))}$ as ${\displaystyle -2C}$ for mathematical convenience in subsequent calculations, and one finally obtains the characteristic relation implied by Eqs(7.a-7.e) that

$${\displaystyle e^{2\psi }=\Phi ^{2}-2C\Phi +1\,.}$$

(7.f)

This relation is important in linearize the Eqs(7.a-7.f) and superpose electrovac Weyl solutions.

Newtonian analogue of metric potential Ψ(ρ,z)

In Weyl's metric Eq(1), ${\textstyle e^{\pm 2\psi }=\sum _{n=0}^{\infty }{\frac {(\pm 2\psi )^{n)){n!))}$; thus in the approximation for weak field limit ${\displaystyle \psi \to 0}$, one has

$${\displaystyle g_{tt}=-(1+2\psi )-{\mathcal {O))(\psi ^{2})\,,\quad g_{\phi \phi }=1-2\psi +{\mathcal {O))(\psi ^{2})\,,}$$

(9)

and therefore

$${\displaystyle ds^{2}\approx -{\Big (}1+2\psi (\rho ,z){\Big )}\,dt^{2}+{\Big (}1-2\psi (\rho ,z){\Big )}\left[e^{2\gamma }(d\rho ^{2}+dz^{2})+\rho ^{2}d\phi ^{2}\right]\,.}$$

(10)

This is pretty analogous to the well-known approximate metric for static and weak gravitational fields generated by low-mass celestial bodies like the Sun and Earth,^[5]

$${\displaystyle ds^{2}=-{\Big (}1+2\Phi _{N}(\rho ,z){\Big )}\,dt^{2}+{\Big (}1-2\Phi _{N}(\rho ,z){\Big )}\,\left[d\rho ^{2}+dz^{2}+\rho ^{2}d\phi ^{2}\right]\,.}$$

(11)

where ${\displaystyle \Phi _{N}(\rho ,z)}$ is the usual Newtonian potential satisfying Poisson's equation ${\displaystyle \nabla _{L}^{2}\Phi _{N}=4\pi \varrho _{N))$, just like Eq(3.a) or Eq(4.a) for the Weyl metric potential ${\displaystyle \psi (\rho ,z)}$. The similarities between ${\displaystyle \psi (\rho ,z)}$ and ${\displaystyle \Phi _{N}(\rho ,z)}$ inspire people to find out the Newtonian analogue of ${\displaystyle \psi (\rho ,z)}$ when studying Weyl class of solutions; that is, to reproduce ${\displaystyle \psi (\rho ,z)}$ nonrelativistically by certain type of Newtonian sources. The Newtonian analogue of ${\displaystyle \psi (\rho ,z)}$ proves quite helpful in specifying particular Weyl-type solutions and extending existing Weyl-type solutions.^[2]

Schwarzschild solution

The Weyl potentials generating Schwarzschild's metric as solutions to the vacuum equations Eq(8) are given by^[2][3][4]

$${\displaystyle \psi _{SS}={\frac {1}{2))\ln {\frac {L-M}{L+M))\,,\quad \gamma _{SS}={\frac {1}{2))\ln {\frac {L^{2}-M^{2)){l_{+}l_{-))}\,,}$$

(12)

where

$${\displaystyle L={\frac {1}{2)){\big (}l_{+}+l_{-}{\big )}\,,\quad l_{+}={\sqrt {\rho ^{2}+(z+M)^{2))}\,,\quad l_{-}={\sqrt {\rho ^{2}+(z-M)^{2))}\,.}$$

(13)

From the perspective of Newtonian analogue, ${\displaystyle \psi _{SS))$ equals the gravitational potential produced by a rod of mass ${\displaystyle M}$ and length ${\displaystyle 2M}$ placed symmetrically on the ${\displaystyle z}$-axis; that is, by a line mass of uniform density ${\displaystyle \sigma =1/2}$ embedded the interval ${\displaystyle z\in [-M,M]}$. (Note: Based on this analogue, important extensions of the Schwarzschild metric have been developed, as discussed in ref.^[2])

Given ${\displaystyle \psi _{SS))$ and ${\displaystyle \gamma _{SS))$, Weyl's metric Eq(1) becomes

$${\displaystyle ds^{2}=-{\frac {L-M}{L+M))dt^{2}+{\frac {(L+M)^{2)){l_{+}l_{-))}(d\rho ^{2}+dz^{2})+{\frac {L+M}{L-M))\,\rho ^{2}d\phi ^{2}\,,}$$

(14)

and after substituting the following mutually consistent relations

$${\displaystyle {\begin{aligned}&L+M=r\,,\quad l_{+}-l_{-}=2M\cos \theta \,,\quad z=(r-M)\cos \theta \,,\\&\rho ={\sqrt {r^{2}-2Mr))\,\sin \theta \,,\quad l_{+}l_{-}=(r-M)^{2}-M^{2}\cos ^{2}\theta \,,\end{aligned))}$$

(15)

one can obtain the common form of Schwarzschild metric in the usual ${\displaystyle \{t,r,\theta ,\phi \))$ coordinates,

$${\displaystyle ds^{2}=-\left(1-{\frac {2M}{r))\right)\,dt^{2}+\left(1-{\frac {2M}{r))\right)^{-1}dr^{2}+r^{2}d\theta ^{2}+r^{2}\sin ^{2}\theta \,d\phi ^{2}\,.}$$

(16)

The metric Eq(14) cannot be directly transformed into Eq(16) by performing the standard cylindrical-spherical transformation ${\displaystyle (t,\rho ,z,\phi )=(t,r\sin \theta ,r\cos \theta ,\phi )}$, because ${\displaystyle \{t,r,\theta ,\phi \))$ is complete while ${\displaystyle (t,\rho ,z,\phi )}$ is incomplete. This is why we call ${\displaystyle \{t,\rho ,z,\phi \))$ in Eq(1) as Weyl's canonical coordinates rather than cylindrical coordinates, although they have a lot in common; for example, the Laplacian ${\displaystyle \nabla ^{2}:=\partial _{\rho \rho }+{\frac {1}{\rho ))\partial _{\rho }+\partial _{zz))$ in Eq(7) is exactly the two-dimensional geometric Laplacian in cylindrical coordinates.

Nonextremal Reissner–Nordström solution

The Weyl potentials generating the nonextremal Reissner–Nordström solution (${\displaystyle M>|Q|}$) as solutions to Eqs(7) are given by^[2][3][4]

$${\displaystyle \psi _{RN}={\frac {1}{2))\ln {\frac {L^{2}-\left(M^{2}-Q^{2}\right)}{\left(L+M\right)^{2))}\,,\quad \gamma _{RN}={\frac {1}{2))\ln {\frac {L^{2}-\left(M^{2}-Q^{2}\right)}{l_{+}l_{-))}\,,}$$

(17)

where

$${\displaystyle L={\frac {1}{2)){\big (}l_{+}+l_{-}{\big )}\,,\quad l_{+}={\sqrt {\rho ^{2}+\left(z+{\sqrt {M^{2}-Q^{2))}\right)^{2))}\,,\quad l_{-}={\sqrt {\rho ^{2}+\left(z-{\sqrt {M^{2}-Q^{2))}\right)^{2))}\,.}$$

(18)

Thus, given ${\displaystyle \psi _{RN))$ and ${\displaystyle \gamma _{RN))$, Weyl's metric becomes

$${\displaystyle ds^{2}=-{\frac {L^{2}-\left(M^{2}-Q^{2}\right)}{\left(L+M\right)^{2))}dt^{2}+{\frac {\left(L+M\right)^{2)){l_{+}l_{-))}(d\rho ^{2}+dz^{2})+{\frac {(L+M)^{2)){L^{2}-(M^{2}-Q^{2})))\rho ^{2}d\phi ^{2}\,,}$$

(19)

and employing the following transformations

$${\displaystyle {\begin{aligned}&L+M=r\,,\quad l_{+}-l_{-}=2{\sqrt {M^{2}-Q^{2))}\,\cos \theta \,,\quad z=(r-M)\cos \theta \,,\\&\rho ={\sqrt {r^{2}-2Mr+Q^{2))}\,\sin \theta \,,\quad l_{+}l_{-}=(r-M)^{2}-(M^{2}-Q^{2})\cos ^{2}\theta \,,\end{aligned))}$$

(20)

one can obtain the common form of non-extremal Reissner–Nordström metric in the usual ${\displaystyle \{t,r,\theta ,\phi \))$ coordinates,

$${\displaystyle ds^{2}=-\left(1-{\frac {2M}{r))+{\frac {Q^{2)){r^{2))}\right)dt^{2}+\left(1-{\frac {2M}{r))+{\frac {Q^{2)){r^{2))}\right)^{-1}dr^{2}+r^{2}d\theta ^{2}+r^{2}\sin ^{2}\theta \,d\phi ^{2}\,.}$$

(21)

Extremal Reissner–Nordström solution

The potentials generating the extremal Reissner–Nordström solution (${\displaystyle M=|Q|}$) as solutions to Eqs(7) are given by^[4] (Note: We treat the extremal solution separately because it is much more than the degenerate state of the nonextremal counterpart.)

$${\displaystyle \psi _{ERN}={\frac {1}{2))\ln {\frac {L^{2)){(L+M)^{2))}\,,\quad \gamma _{ERN}=0\,,\quad {\text{with))\quad L={\sqrt {\rho ^{2}+z^{2))}\,.}$$

(22)

Thus, the extremal Reissner–Nordström metric reads

$${\displaystyle ds^{2}=-{\frac {L^{2)){(L+M)^{2))}dt^{2}+{\frac {(L+M)^{2)){L^{2))}(d\rho ^{2}+dz^{2}+\rho ^{2}d\phi ^{2})\,,}$$

(23)

and by substituting

$${\displaystyle L+M=r\,,\quad z=L\cos \theta \,,\quad \rho =L\sin \theta \,,}$$

(24)

we obtain the extremal Reissner–Nordström metric in the usual ${\displaystyle \{t,r,\theta ,\phi \))$ coordinates,

$${\displaystyle ds^{2}=-\left(1-{\frac {M}{r))\right)^{2}dt^{2}+\left(1-{\frac {M}{r))\right)^{-2}dr^{2}+r^{2}d\theta ^{2}+r^{2}\sin ^{2}\theta \,d\phi ^{2}\,.}$$

(25)

Mathematically, the extremal Reissner–Nordström can be obtained by taking the limit ${\displaystyle Q\to M}$ of the corresponding nonextremal equation, and in the meantime we need to use the L'Hospital rule sometimes.

Remarks: Weyl's metrics Eq(1) with the vanishing potential ${\displaystyle \gamma (\rho ,z)}$ (like the extremal Reissner–Nordström metric) constitute a special subclass which have only one metric potential ${\displaystyle \psi (\rho ,z)}$ to be identified. Extending this subclass by canceling the restriction of axisymmetry, one obtains another useful class of solutions (still using Weyl's coordinates), namely the conformastatic metrics,^[6][7]

$${\displaystyle ds^{2}\,=-e^{2\lambda (\rho ,z,\phi )}dt^{2}+e^{-2\lambda (\rho ,z,\phi )}{\Big (}d\rho ^{2}+dz^{2}+\rho ^{2}d\phi ^{2}{\Big )}\,,}$$

(26)

where we use ${\displaystyle \lambda }$ in Eq(22) as the single metric function in place of ${\displaystyle \psi }$ in Eq(1) to emphasize that they are different by axial symmetry (${\displaystyle \phi }$-dependence).

Weyl vacuum solutions in spherical coordinates

Weyl's metric can also be expressed in spherical coordinates that

$${\displaystyle ds^{2}\,=-e^{2\psi (r,\theta )}dt^{2}+e^{2\gamma (r,\theta )-2\psi (r,\theta )}(dr^{2}+r^{2}d\theta ^{2})+e^{-2\psi (r,\theta )}\rho ^{2}d\phi ^{2}\,,}$$

(27)

which equals Eq(1) via the coordinate transformation ${\displaystyle (t,\rho ,z,\phi )\mapsto (t,r\sin \theta ,r\cos \theta ,\phi )}$ (Note: As shown by Eqs(15)(21)(24), this transformation is not always applicable.) In the vacuum case, Eq(8.b) for ${\displaystyle \psi (r,\theta )}$ becomes

$${\displaystyle r^{2}\psi _{,\,rr}+2r\,\psi _{,\,r}+\psi _{,\,\theta \theta }+\cot \theta \cdot \psi _{,\,\theta }\,=\,0\,.}$$

(28)

The asymptotically flat solutions to Eq(28) is^[2]

$${\displaystyle \psi (r,\theta )\,=-\sum _{n=0}^{\infty }a_{n}{\frac {P_{n}(\cos \theta )}{r^{n+1))}\,,}$$

(29)

where ${\displaystyle P_{n}(\cos \theta )}$ represent Legendre polynomials, and ${\displaystyle a_{n))$ are multipole coefficients. The other metric potential ${\displaystyle \gamma (r,\theta )}$is given by^[2]

$${\displaystyle \gamma (r,\theta )\,=-\sum _{l=0}^{\infty }\sum _{m=0}^{\infty }a_{l}a_{m}{\frac {(l+1)(m+1)}{l+m+2)){\frac {P_{l}P_{m}-P_{l+1}P_{m+1)){r^{l+m+2))}\,.}$$

(30)