## Aspects of electrostatics in BTZ geometries

Y. Herrera 1 ∗ V. Hurovich 1 † O. Santill´n a 2 ‡ and C. Simeone 1 § .

- 1

Departamento de F´ ısica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires and IFIBA, CONICET, Cuidad Universitaria, Buenos Aires 1428, Argentina.

2 CONICET-Instituto de Investigaciones Matem´ticas Luis Santal´, a o Ciudad Universitaria Pab. I, Buenos Aires 1428, Argentina.

## Abstract

In the present paper the electrostatic of charges in non rotating BTZ black hole and wormhole space times is studied. In particular, the self force of a point charge in the geometry is characterized analitically. The differences between the self force in both cases is a theoretical experiment for distinguishing both geometries, which otherwise are locally indistinguishable. This idea was applied before to four and higher dimensional black holes by the present and other authors. However, the particularities of the BTZ geometry makes the analysis considerable more complicated than usual electrostatic in a flat geometry, and its even harder than its four dimensional counterparts. First, the BTZ space times are not asymptotically flat but instead asymptotically AdS. In addition, the relative distance d r, r ( +1) between two particles located at a radius r and r +1 in the geometry tends to zero when r →∞ . This behavior, which is radically different in a flat geometry, changes the analysis of the asymptotic conditions for the electrostatic field. In addition, there are no summation formulas that allow a closed analytic expression for the self force. We find a method to calculate such force in series, and the resulting expansion is convergent to the real solution. However, we suspect that the convergence is not uniform. In other, for points that are far away from the black hole the calculation of the force requires higher order summation. These subtleties are carefully analyzed in the paper, and it is shown that they lead to severe problems when calculating the self force for asymptotic points in the geometry.

## 1. Introduction

Electrodynamics in General Relativity is described by the Maxwell equations in curved space-time [1]. A freely falling observer in such background would write the same equations valid for Minkowski space-time; however, these equations must have a different solution, because the curved geometry imposes a different asymptotic behavior than the flat one. In particular, the electric field around a static point charge in a curved background is not spherically symmetric in general, and this gives a non zero electrostatic self-force on the charge.

One of the earliest studies on the electrostatic self-force on static charges induced by a curved background was that on a Schwarzschild black hole geometry [3]. In that reference it was shown that

∗ e-mail: yagogaoc@hotmail.com

† e-mail: valehurovich@gmail.com

‡ e-mail: firenzecita@hotmail.com, osantil@dm.uba.ar

§ e-mail: csimeone@df.uba.ar

the self-force on a charge q is repulsive, i.e. it points outwards from the black hole, and the functional dependence on the position is given by

$$f \sim \frac { m q ^ { 2 } } { r ^ { 3 } },$$

where 2 m is the horizon radius of the black hole and r is the Schwarzschild radial coordinate of the charge. This result was first obtained within the framework of linearized general relativity [2], and was later recovered working within the full theory [3].

After the publication of these leading works the study of the self-interaction of a charge was extended to other geometries. A notable result was the self-force on a charge in the vicinity of a straight cosmic string arising from symmetry breaking in a system composed by a complex scalar field coupled to a gauge field [4]. The associated geometry is locally flat but includes a deficit angle determined by µ , the mass per unit length of the string [5]. The self-force in this case points outwards from the cosmic string and is proportional to µ/r 2 . This non null self-force in a locally flat background is of great interest because it shows how the global properties of a manifold (in this case, the existence of a deficit angle) are revealed by the electromagnetic field of the charge.

The results described above together with the calculation of the self-force on a point charge in a wormhole space-time [6], which turned out to be attractive, i.e. towards the wormhole throat, suggested the possibility of detecting thin-shell wormholes by means of electrostatics. Differing from well-known wormholes of the Morris-Thorne type [7] which are supported by non localized exotic matter, thin-shell wormhole geometries are supported by a shell of exotic matter located at the wormhole throat [8]. The throat connects two (equal or different) geometries which can be those of other astrophysical objects. For example, Schwarzschild thin-shell wormholes connect two exterior (that is, beyond the horizon) non charged black hole space-times; hence the geometry at each side of the throat is locally identical to the exterior of a black hole geometry. However, the topology of the wormhole geometry is non trivial, thus the global properties are essentially different in each case.

Inspired by the previous discussion, the philosophy of the present work is that global aspects, such as the existence of a throat or not, may be revealed by studying electrodynamics, in particular, by the electrostatic self-force on a point charge. In our recent article [9] this proposal was developed and applied it to the case of wormholes with a cylindrical throat which are mathematically constructed by removing the regions r &lt; a of two gauge cosmic string manifolds and pasting the two regions at r ≥ a . The self-force on a charge in the cylindrical wormhole geometry was calculated, and compared it with the self-force on a charge in the vicinity of a gauge cosmic string. The result is that the force in the wormhole case can be attractive or repulsive depending of the position of the charge; this result would then allow an observer to distinguish between two geometries which are locally equal. The same argument was applied to the Schwarzschild case by the authors in [10]. Related works are also [13]-[16].

It should be mentioned that there exist some works related to these ideas. For instance, in [11], the authors considered a minimally coupled scalar charge and an electromagnetic charge when a

Schwarzschild black hole interior is replaced by a material body and found that the leading term in a larger expansion of the force was independent of the central body type. Nevertheless, when the scalar charge is not minimally coupled, the self-force is dependent on the composition of the body. Another work in the same line is [12], where a spherical ball of perfect fluid in hydrostatic equilibrium with rest mass density and pressure related by some polytropic equations of state is considered. The authors found that the leading term of the force is universal and does not distinguish the internal body structure, but the next-to-leading order term is sensible to the equation of state. Thus the self-force distinguishes the body composition.

In the present work our studies about electrostatics in black hole geometries is extended to the three dimensional case, which is not a completely explored area. The natural candidates to consider are the BTZ black hole and wormhole. These geometries, although tridimensional and non realistic, have several features that makes them an interesting test laboratory. First of all, both have a negative cosmological constant Λ &lt; 0, which corresponds to an attraction instead of repulsion. On the other hand, their metric is not asymptotically flat, but asymptotically anti De Sitter. In addition there exist a radial coordinate r such that the circles of r constant have perimeter 2 πr , but the relative distance d r, r ( + 1) between points located on the same radial line at positions r and r + 1 goes to zero as r →∞ . This behavior is not characteristic for simple black holes in higher dimensions and is a consequence of the attractive cosmological constant term. This behavior has consequences on the boundary conditions of the electrostatic problem. These consequences will be elaborated in the paper.

It should be remarked that there are no summation formulas allowing to find a closed analytic for the self force, thus it is given in series expansion. But the series expansion of the singular part of the electrostatic field is rather complicated, and can be achieved by certain specific parameterization of the radial distance, which is explained in detail in section 7.1. These tricks result in a series expansion for the total force that is convergent to the real one. However, we argue that the convergence to the real solution is not uniform, in other words, as larger the coordinate r 0 of the charge becomes, the larger the quantity of terms that has to be summed in order to approximate the self force at that point. This results in a problem when truncating the series at r →∞ .

The present work is organized as follows. Section 2 contains a brief description of the BTZ black hole geometry. Section 3 contains a description of electrostatic equations in the geometry and the problems for fixing the boundary conditions for the physical solutions. Section 4 contains the expression for the electrostatic field for the BTZ black hole and wormhole. Section 5 contains a review of Synge calculus, which is a relevant tool for calculating the singular part of the electrostatic Green function of the geometry. In section 6 this singular part is calculated for the BTZ local geometry. In section 7 a series expansion for the electrostatic self force is calculated, and the problems mentioned above about the convergence at the asymptotic boundary is described. Section 8 contains and 9 contains interpretations of the results, which are rather non trivial. Section 10 is a summary of the obtained results. In the appendix there are collected some useful formulas which are applied along

the text.

## 2. The BTZ black hole

Since the appearance of the seminal works [17]-[18], General Relativity in (2+1) dimensions became a widely analyzed model for exploring classical and quantum gravity, since it is recognized as a useful laboratory for studying real system properties in (3+1) dimensions. In (2+1) dimensions GR there is no newtonian limit and there are no local degrees of freedom (that is, there are no gravitational waves in the classical theory or gravitons in the quantum theory). It came as a surprise for some then when the black hole BTZ solution was found [19]. This black hole has important differences with the Schwarzschild and Kerr black holes: it is asymptotically anti De Sitter and not asymptotically flat, and does not have any curvature singularity at the origin. Nevertheless, it is clearly a black hole: it has an even horizon and (in the rotating case) an internal horizon, and thermodynamical properties similar to black holes in (3+1) dimensions.

The BTZ solution is well known, but in order to fix the conventions a short description of the local and global properties of the geometry will be given. The discussion is not exhaustive, but focused in the aspects that are more important for the present work.

## 2.1 Parameters of the solution

The BTZ black hole is a solution of the Einstein field equations in (2+1) dimensions with cosmological constant Λ &lt; 0, which bears some similarities with black hole solutions in four dimensions [19]-[20]. These are for instance as the presence of a event horizon, an inner horizon and an ergosphere. Also, it has a non vanishing Hawking temperature and interesting thermodynamical properties [21]. Despite these similarities, there are several differences between BTZ black holes and Schwarzschild or Kerr ones. The later are asymptotically flat, the BTZ solution instead is asymptotically anti De Sitter. Furthermore, the BTZ solution does not have a singularity at the origin. But since the BTZ structure is simpler than its four dimensional counterparts, it may be a good testing laboratory for making exact calculations.

The local form of the BTZ solution is well known, but in order to fix the notation a brief review of these solutions will be given. Starting with the three dimensional action [20]

$$I = \frac { 1 } { 2 \pi } \int \sqrt { - g } [ R + 2 l ^ { - 2 } ] d ^ { 2 } x d t + B,$$

with B is a surface term and l is related to the cosmological constant by -Λ = l -2 , it follows that the extremal solutions corresponding to g µν ( x, t ) variations are given by the Einstein equations

$$R _ { \mu \nu } - \frac { 1 } { 2 } g _ { \mu \nu } ( R + 2 l ^ { - 2 } ) = 0,$$

which, in three dimensions only, completely determine the Riemann tensor as

$$R _ { \mu \nu \lambda \rho } = - l ^ { - 2 } ( g _ { \mu \lambda } g _ { \nu \rho } - g _ { \nu \lambda } g _ { \mu \rho } ).$$

This solution corresponds to a symmetric space with negative curvature. If one restricts the attention to solutions possessing a rotational Killing vector ∂/∂θ and a time like Killing vector ∂/∂t , then by an specific choice of the radial coordinate it follows that the line element is given by

$$d s ^ { 2 } = - N ^ { 2 } d t ^ { 2 } + N ^ { - 2 } d r ^ { 2 } + r ^ { 2 } ( N ^ { \theta } d t + d \theta ) ^ { 2 },$$

with N 2 ( r ) and N θ ( r ) the following radial functions

$$N ^ { 2 } ( r ) = - M + \frac { r ^ { 2 } } { l ^ { 2 } } + \frac { J ^ { 2 } } { 4 r ^ { 2 } },$$

$$N ^ { \theta } ( r ) = - \frac { J } { 2 r ^ { 2 } }.$$

The range of the coordinates is -∞ &lt; t &lt; ∞ , 0 &lt; r &lt; ∞ and 0 ≤ θ ≤ 2 π . The two integration constants in (2.5) and (2.6) are M and J and correspond to the mass and angular momentum of the solutions respectively [20].

The BTZ space time is not asymptotically flat. For large radial values r →∞ the metric becomes

$$d s \rightarrow - ( \frac { r } { l } ) ^ { 2 } d t ^ { 2 } + ( \frac { r } { l } ) ^ { - 2 } d r ^ { 2 } + r ^ { 2 } d \theta ^ { 2 },$$

which shows that this solution is asymptotically anti De Sitter.

The function N r ( ) vanish for the following two r values

$$r _ { \pm } = l \left [ \frac { M } { 2 } \left ( 1 \pm \sqrt { 1 - \left ( \frac { J } { M l } \right ) ^ { 2 } } \right ) \right ] ^ { 1 / 2 }.$$

The value r + corresponds to the horizon of the black hole. It exist if the following inequalities are satisfied

$$M > 0, \quad \ | J | \leq M l.$$

In the extreme case | J | = Ml , both roots of N 2 = 0 coalesce into one. The mass M and the angular momentum J can be expressed in terms of r ± as

$$M = \frac { r _ { + } ^ { 2 } + r _ { - } ^ { 2 } } { l ^ { 2 } }, \quad \ J = \frac { 2 r _ { + } r _ { - } } { l }.$$

For large l the exterior horizon tends to infinite and only the interior remains. The vacuum state is obtained when the black hole disappear, and this corresponds to take the horizon radius to zero. This is equivalent of taking M → 0, which implies J → 0 due to (2.9). In this case

$$d s _ { v a c } ^ { 2 } = - ( \frac { r } { l } ) ^ { 2 } d t ^ { 2 } + ( \frac { r } { l } ) ^ { - 2 } d r ^ { 2 } + r ^ { 2 } d \theta ^ { 2 }.$$

When M becomes negative, the solutions studied in [22] are found. The conical singularity that they posses is a naked one, such as the one as a black hole with negative mass in (3+1) dimension. Such value should be excluded from the spectrum. Nevertheless there exists an exceptional case. When M = -1 and J = 0 the naked singularity disappears. There is no horizon in this case, but also no singularity to hide. The solution corresponding to this regime is

$$d s ^ { 2 } = - ( 1 + ( r / l ) ^ { 2 } ) d t ^ { 2 } + ( 1 + ( r / l ) ^ { 2 } ) ^ { - 1 } d r ^ { 2 } + r ^ { 2 } d \theta ^ { 2 },$$

and is AdS as well.

## 2.2 Particular properties of the non rotating geometry

In this section some properties of the BTZ black hole will be pointed out, which will be relevant when analyzing the electrostatic properties of charges in the geometry. In the present work the non rotating case J = 0 will be only considered. The rotating case is leaved for a forthcoming publication.

An observation which will be of importance for interpreting the results of the present work is that, in the non rotating BTZ geometry, the distance d r, r ( +1) between two points with the same θ values and lying on the circles r and ( r +1) decreases when r increases. To see this, consider for simplicity the case M = l = 1. Then the distance from a point with coordinate r to the horizon r h is

$$d = \log ( r + \sqrt { r ^ { 2 } - 1 } ),$$

which can be inverted to give r = cosh( d ). When d &gt;&gt; 1 it follows that r ∼ e d . If two points lying on the same line θ = θ 0 are at positions r and r + δr , then the last formula gives that

$$\delta r = e ^ { d + \delta d } - e ^ { d } = e ^ { d } ( e ^ { \delta d } - 1 ),$$

which leads to

$$\delta d = \log \left ( 1 + \frac { 1 } { e ^ { d } } \right ) = \log \left ( 1 + \frac { 1 } { r } \right ).$$

From here it is seen that for r &gt;&gt; 1, which implies going far from the horizon d &gt;&gt; 1, the true distance δd between these points goes to zero δd → 0. This particularity holds for other values of M and will play a significant role in the interpretation of our results.

## 3. The equations of electrostatics in BTZ space times

In the present section the Maxwell equations corresponding to a static charge q located at r 0 and θ 0 = 0 in a BTZ black hole will be derived. The effect of the curved geometry is to deform the field lines and, as a consequence, the charge q experience a self-force due to its own electric field. As it will be shown below, the Maxwell equations are separable in this case. Nevertheless the analysis of the physical and unphysical solutions is more involved than in the flat case due to the particularities of the geometry mentioned in the previous section, in particular, the behavior (2.15). These aspects are carefully examined below and a criteria for discarding unphysical solutions is found.

## 3.1 Separation of variables

The Maxwell equations in three dimensional curved space times in natural units are given by [1]

$$\begin{array} { c c c c c c } & & \cdot & & \cdot & & \cdot & \cdot \\ \frac { 1 } { \sqrt { - g } } \partial _ { \alpha } ( \sqrt { - g } g ^ { \mu \alpha } g ^ { \nu \beta } F _ { \mu \nu } ) = 2 \pi j ^ { \alpha }, & & & \cdot & & \cdot & \cdot \\ \epsilon _ { \beta \gamma \delta } \partial _ { \beta } F _ { \gamma \delta } = 0. & & & & & & \cdot & & \\ \end{array}$$

Here F µν = ∂ A µ ν -∂ A ν µ is the field strength tensor, A µ is the vector potential and j α the three current in the geometry. For an static charge q in front of the non rotating geometry one has that

$$j ^ { t } = \frac { q } { r } \delta ( r - r _ { 0 } ) \delta ( \theta - \theta _ { 0 } ),$$

with ( r , θ 0 0 ) the coordinates of the position of the charge. The Maxwell equations (3.16) in this situation reduce to

$$- \, \partial _ { r } ( r F _ { t r } ) + \frac { r } { M r ^ { 2 } - \frac { r ^ { 4 } } { l ^ { 2 } } } \partial _ { \theta } F _ { t \theta } = 2 \pi r j ^ { t },$$

$$\partial _ { \theta } \left [ \frac { 1 } { r } ( M - \frac { r ^ { 2 } } { l ^ { 2 } } ) F _ { r \theta } \right ] = 0,$$

$$\partial _ { r } \left [ \frac { 1 } { r } ( M - \frac { r ^ { 2 } } { l ^ { 2 } } ) F _ { r \theta } \right ] = 0.$$

Assuming that the vector A µ is time independent it follows from these equations that there exist a gauge in which only the component A t is non zero, and the three (3.18) reduce to the following single equation

$$- \, \partial _ { r } ( r \partial _ { r } A _ { t } ) + \frac { r } { M r ^ { 2 } - \frac { r ^ { 4 } } { l ^ { 2 } } } \partial _ { \theta } ^ { 2 } A _ { t } = 2 \pi r \, j ^ { t }.$$

Outside the position of the charge this equation is homogeneous and can be solved by variable separation by postulating

$$A _ { t } ( r, \theta ) = R ( r ) \Theta ( \theta ).$$

When this is inserted into (3.20) it is obtained that

$$\Theta ( \theta ) = \exp [ i n ( \theta - \theta _ { 0 } ) ],$$

where n is an integer due to the periodicity on θ , and the following equation for R r ( )

$$r ^ { 2 } ( 1 - \frac { r ^ { 2 } } { M l ^ { 2 } } ) \partial _ { r } ^ { 2 } R ( r ) + r ( 1 - \frac { r ^ { 2 } } { M l ^ { 2 } } ) \partial _ { r } R ( r ) + \frac { n ^ { 2 } R ( r ) } { M } = 0.$$

By further defining the horizon radius r 2 h = Ml 2 and making the variable change r 2 → xr 2 h it is transformed into

$$\Pi \omega _ { x } ^ { 2 } ( x - 1 ) \partial _ { x } ^ { 2 } R ( x ) + x ( x - 1 ) \partial _ { x } R ( x ) - \frac { n ^ { 2 } R ( x ) } { 4 M } = 0.$$

This equation has two regular singular points, which corresponds to the horizon x = 1 and the infinite x →∞ . In order to analyze the behavior at the infinite it is customary to make the change of variables x → 1 u which transforms the last equation into

$$u ( 1 - u ) \partial _ { u } ^ { 2 } R + ( 1 - u ) \partial _ { u } R - \frac { n ^ { 2 } R } { 4 M } = 0.$$

The equation (3.25) is a particular case of the hypergeometric one

$$u ( 1 - u ) R ^ { \prime \prime } + [ \gamma - ( 1 + \alpha + \beta ) u ] R ^ { \prime } - \alpha \beta R = 0,$$

corresponding to the particular values

$$\gamma = 1, \quad \alpha + \beta = 0, \quad \beta ^ { 2 } = - n ^ { 2 } / 4 M.$$

It is important to remark that the change of variables u = 1 x just performed is regular in the exterior region r &gt; r h of the black hole, which is the region which we are interested in.

## 3.2 Solutions centered around the infinite

Having derived the equation (3.26) which characterizes the radial behavior of the electrostatic potential A t , the next task is to find their solutions. Since it is a linear equation of second order, it has two independent solutions. The most elementary one, which is centered around u = 0 ( r →∞ ), is given by the hypergeometric series [23]-[26]

$$f _ { n } = \, _ { 2 } F _ { 1 } ( \alpha _ { n }, - \alpha _ { n } ; 1 ; u ) = \sum _ { m = 0 } ^ { \infty } \frac { ( \alpha _ { n } ) _ { m } ( - \alpha _ { n } ) _ { m } } { ( m! ) ^ { 2 } } u ^ { m },$$

where

$$\alpha _ { n } = \frac { i n } { 2 \sqrt { M } },$$

and the Pochhammer symbols ( a ) n are defined by

$$( \alpha ) _ { m } = \alpha ( \alpha + 1 ) ( \alpha + 2 ) \dots ( \alpha + m - 1 ), \quad ( \alpha ) _ { 0 } = 1.$$

The elementary D' Alembert principle shows that this series is convergent for | u &lt; | 1. Besides, when

$$R e ( \gamma - \alpha - \beta ) > 0,$$

the series is also convergent in | u | = 1 [25]. This condition is satisfied in our situation since β = -α and γ = 1. The zone | u &gt; | 2 corresponds to the inner part of the black hole, which is of no interest to us.

The transformation R → u 1 -γ R applied to (3.26) transforms it into another hypergeometric equation but induces a parameter transformation ( α, β, γ ) → ( α -γ +1 , β -γ +1 2 , -γ ). Therefore, in general, the function

$$g _ { n } = u ^ { 1 - \gamma } F ( \alpha - \gamma + 1, \beta - \gamma + 1, 2 - \gamma, u ),$$

is also a solution of (3.26). Nevertheless, when γ = 1, as in our case, this solution is equivalent to f n , and gives no new information. In these particular cases, a new solution is obtained by postulating a series of the form

$$g _ { n } = f _ { n } \log u + \sum _ { n = 0 } ^ { \infty } c _ { r } u ^ { r },$$

with c r constant coefficients to be determined. By inserting this into (3.26) the following recurrence for c r is obtained

$$( r + 1 ) ^ { 2 } c _ { r } - r ( \alpha + \beta + 1 ) c _ { r + 1 } + \frac { ( \alpha \beta - \alpha - \beta ) ( \alpha ) _ { r } ( \beta ) _ { r } } { r! ( r + 1 )! } = 0,$$

which, when solved explicitly, gives the following solution

$$g _ { n } = \, _ { 2 } F _ { 1 } ( \alpha _ { n }, - \alpha _ { n } ; 1 ; u ) \log u + \sum _ { m = 0 } ^ { \infty } \frac { ( \alpha _ { n } ) _ { m } ( - \alpha _ { n } ) _ { m } } { ( m! ) ^ { 2 } } \, u ^ { m } \, S _ { n, m },$$

with

$$S _ { n, m } = \sum _ { k = 0 } ^ { m - 1 } \left ( \frac { 1 } { k + \alpha _ { n } } + \frac { 1 } { k - \alpha _ { n } } - \frac { 2 } { k + 1 } \right ).$$

Alternatively, this second solution may be though as the limit [25]

$$g _ { n } = \lim _ { \gamma \rightarrow 1 } \frac { u ^ { 1 - \gamma } \, _ { 2 } F _ { 1 } ( \alpha - \gamma + 1, \beta - \gamma + 1, 2 - \gamma, u ) - \, _ { 2 } F _ { 1 } ( \alpha, \beta, \gamma, u ) } { \gamma - 1 }.$$

An important property of hypergeometric functions is the following [23]

$$\ _ { 2 } F _ { 1 } ^ { \prime } ( \alpha, \beta, \gamma, u ) = \frac { \alpha \beta } { \gamma } \, _ { 2 } F _ { 1 } ( \alpha + 1, \beta + 1, \gamma + 1, u ),$$

which express it derivatives in terms of other hypergeometric functions. From this property and the definition (3.27)-(3.32) for f ν and g ν it follows that

$$\frac { \partial f _ { n } } { \partial r } = \frac { 2 \alpha _ { n } ^ { 2 } u ^ { 3 / 2 } } { r _ { h } } \, _ { 2 } F _ { 1 } ( \alpha _ { n } + 1, - \alpha _ { n } + 1, 2, u ),$$

$$\frac { \partial g _ { n } } { \partial r } = - \frac { 2 } { r _ { h } } \, u ^ { 3 / 2 } \left [ u ^ { - 1 } \, _ { 2 } F _ { 1 } ( \alpha _ { n }, - \alpha _ { n }, 1, u ) - \log ( u ) \, \alpha _ { n } ^ { 2 } \, _ { 2 } F _ { 1 } ( 1 + \alpha _ { n }, 1 - \alpha _ { n }, 2, u ) \, \quad \ \ ( 3. 3 7 )$$

$$+ \sum _ { m = 1 } ^ { \infty } \frac { ( \alpha _ { n } ) _ { m } ( - \alpha _ { n } ) _ { m } } { ( m! ) ^ { 2 } } m \, u ^ { m - 1 } S _ { n, m } \Big.$$

In deriving this formulas the definition u = r 2 h r -2 was taken into account. These formulas will be useful when evaluating the electrostatic field of the charge q as derivatives of the potential A t .

The behavior when r &gt;&gt; r h (which corresponds to u → 0) of the solutions is directly inferred from their definition, the result is

$$f _ { n } \to 1, \quad u \to 0,$$

$$g _ { n } \sim \log ( u ) \to - \infty, \quad u \to 0.$$

The behavior of their derivatives for large r &gt;&gt; r h is inferred by taking into account the following elementary limits

$$\lim _ { u \to 0 } F ( \alpha _ { n }, - \alpha _ { n }, 1, u ) = 1,$$

$$\lim _ { u \to 0 } \log ( u ) u ^ { n } = 0, \quad \forall \ n > 0,$$

$$\lim _ { u \to 0 } F ( \alpha _ { n } + 1, - \alpha _ { n } + 1, 2, u ) < \infty,$$

the last limit follows from the fact that any hypergeometric function is convergent at u = 0. These limits, together with (3.36)-(3.37) show that

$$\frac { d } { d r } f _ { n } \sim \frac { 1 } { r ^ { 3 } }, \quad r \to \infty$$

$$\frac { d } { d r } g _ { n } \sim \frac { 1 } { r }, \quad r \to \infty.$$

Thus none of the derivatives of the solutions is divergent at the asymptotic region. Note that this behavior is in contrast with ordinary electrodynamics in R 2 or R 3 , where there always exist a solution whose electrostatic field is divergent at infinite and is discarded in physical problems. This fact will play a crucial role in our analysis.

Consider now the behavior near the horizon r → r h or u → 1. Both solutions (3.27) y (3.32) are both finite for | u | = 1 since

$$f _ { n } \to \, _ { 2 } F _ { 1 } ( \alpha _ { n }, - \alpha _ { n }, 1, 1 ) = \frac { 1 } { \Gamma ( 1 - \alpha _ { n } ) \Gamma ( 1 + \alpha _ { n } ) },$$

$$g _ { n } \rightarrow \sum _ { m = 0 } ^ { \infty } \frac { ( \alpha _ { n } ) _ { m } ( - \alpha _ { n } ) _ { m } } { ( m! ) ^ { 2 } } S _ { n, m } < \infty,$$

the second inequality follows from the DAlembert criteria for series. More specifically, the function S n,m defined in (3.33) can be approximated by an integral whose result is

$$S _ { n, m } \sim \log \frac { m ^ { 2 } - \alpha _ { n } ^ { 2 } } { ( m + 1 ) ^ { 2 } } - \log ( - \alpha _ { n } ^ { 2 } ),$$

and remembering that α n is purely imaginary it follows that S n,m &lt; ∞ for all m . Therefore (3.45) is

$$| g _ { n } | \rightarrow | \sum _ { m = 0 } ^ { \infty } & \frac { ( \alpha _ { n } ) _ { m } ( - \alpha _ { n } ) _ { m } } { ( m! ) ^ { 2 } } S _ { n, m } | < | S _ { n, m } ^ { m a x } \sum _ { m = 0 } ^ { \infty } \frac { ( \alpha _ { n } ) _ { m } ( - \alpha _ { n } ) _ { m } } { ( m! ) ^ { 2 } } | \\ & = | S _ { n, m } ^ { m a x } \, _ { 2 } F _ { 1 } ( \alpha _ { n }, - \alpha _ { n }, 1, 1 ) | < \infty,$$

where in the last step (3.40) and the definition (3.27) has been taken into account. This shows that (3.45) holds. The derivatives f n with respect to r involve functions of the form 2 F 1 (1+ α , n 1 -α , n 2 , u ), which do not satisfy (3.28). This means that d dr f n is divergent in the horizon u = 1. The analysis for d dr g n is more involved. The first term (3.37) is convergent. The second is also convergent, but the third is divergent. Thus the final result is that

$$\frac { d } { d r } f _ { n } \to \infty, \quad r \to r _ { h }$$

$$\frac { d } { d r } g _ { n } \to \infty, \quad r \to r _ { h }.$$

It is not easy to work with functions with this divergent behavior. Fortunately, there exist a linear combination αf n + βg n of both solutions whose derivative is convergent at | u | = 1. This combination can be found by considering the set of solutions centered at the horizon u = 1.

## 3.3 Solutions centered around the horizon

As it was mentioned above, the equation (3.26) has three regular singular points. The solutions ( f n , g n ) found in the previous section are centered around the regular singular point u = 0, which corresponds to the asymptotic region, and are convergent in the interval 0 &lt; u &lt; | | 1. Consider now a set of solutions ( h , k n n ) centered around the regular singular point u = 1. These solutions will be convergent, as it will be shown below, in the interval 0 &lt; u &lt; 2, in particular for r h &lt; r &lt; ∞ . This means that in the overlapping region 0 &lt; u &lt; 1, which is r h &lt; r &lt; ∞ both sets ( f n , g n ) and ( h n , k n ) constitute a basis of solutions, therefore there should exist a relation of the form

$$h _ { n } = a \, f _ { n } + b \, g _ { n },$$

$$k _ { n } & = c \, f _ { n } + d \, g _ { n },$$

valid in the overlapping region, with a, b, c, d some constant coefficients. These coefficients can be found by evaluating these equalities and its first derivatives in an arbitrary point r 0 inside the overlapping region, the result is

$$\lambda \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \, \widehat { \lambda } \quad \, b = \frac { W ( h _ { n }, f _ { n } ) ( r _ { 0 } ) } { W ( f _ { n }, g _ { n } ) ( r _ { 0 } ) }.$$

$$c = \frac { W ( k _ { n }, g _ { n } ) ( r _ { 0 } ) } { W ( f _ { n }, g _ { n } ) ( r _ { 0 } ) }, \quad d = \frac { W ( k _ { n }, f _ { n } ) ( r _ { 0 } ) } { W ( f _ { n }, g _ { n } ) ( r _ { 0 } ) }.$$

Here W f,g ( ) = g r ∂ f ( ) r ( r ) -f ( r ∂ g r ) r ( ) is the wronskian of the two functions f n and g n . Naturally, the value of a, b, c, d does not depend on the choice of r 0 .

A method for finding the solutions ( h , k n n ) is the following. Consider the change of variables s = 1 -u . The equation (3.26) in this variable takes the form

$$s ( 1 - s ) \tilde { R } ^ { \prime \prime } ( s ) + s \tilde { R } ^ { \prime } ( s ) - \frac { n ^ { 2 } } { 4 M } \tilde { R } ( s ) = 0.$$

Clearly, the solutions of (3.26) around u = 1 correspond to solutions of (3.53) around s = 0. The equation (3.53) is an hypergeometric one with parameters

$$\gamma = 0, \quad \alpha _ { n } = \frac { i n } { 2 \sqrt { M } }, \quad \beta _ { n } = \frac { - i n } { 2 \sqrt { M } }.$$

Its solutions are given by

$$h _ { n } ( r ) = s _ { \, 2 } F _ { 1 } ( 1 + \alpha _ { n }, 1 - \alpha _ { n }, 2, s ),$$

$$k _ { n } ( r ) = s \ln ( s ) \, _ { 2 } F _ { 1 } ( 1 + \alpha _ { n }, 1 - \alpha _ { n }, 2, s ) + s \sum _ { m = 1 } ^ { \infty } \frac { ( 1 + \alpha _ { n } ) _ { m } ( 1 - \alpha _ { n } ) _ { m } } { ( 2 ) _ { m } ( m! ) } \, s ^ { m } \tilde { S } _ { n, m } - \frac { 1 } { \alpha _ { n } ^ { 2 } }, \quad \ ( 3. 5 5 )$$

where α n is the same as before and

$$s = 1 - u = 1 - \frac { r _ { h } ^ { 2 } } { r ^ { 2 } }.$$

Note that non of the solutions (3.54)-(3.55) are simply hypergeometric functions. This happens for some special choice of parameters, such as in our case. In addition

$$\tilde { S } _ { n, m } = \sum _ { k = 0 } ^ { m - 1 } \left ( \frac { 1 } { k + 1 + \alpha _ { n } } + \frac { 1 } { k + 1 - \alpha _ { n } } - \frac { 1 } { k + 1 } - \frac { 1 } { k + 2 } \right ) = \psi ( m + 1 + \alpha _ { n } ) - \psi ( 1 + \alpha _ { n } )$$

$$+ \psi ( m + 1 - \alpha _ { \nu } ) - \psi ( 1 - \alpha _ { \nu } ) - \psi ( m + 2 ) + \psi ( 2 ) - \psi ( m + 1 ) + \psi ( 1 ),$$

has been introduced, with ψ a ( )the digamma function. This new variable change maps the exterior region of the black hole to 0 ≤ s ≤ 1. The line s = 0 corresponding to the event horizon and s = 1 corresponds to the asymptotic region.

The derivatives with respect to r are given by

$$\frac { d } { d r } h _ { n } = \frac { 2 } { r _ { h } } \left ( 1 - s \right ) ^ { 3 / 2 } \left [ \tilde { F } _ { n } ( s ) + s \tilde { F } _ { n } ^ { \ \prime } ( s ) \right ]$$

$$\frac { d } { d r } k _ { n } = \frac { 2 } { r _ { h } } \left ( 1 - s \right ) ^ { 3 / 2 } \left [ \left ( 1 + \ln s \right ) \tilde { F } _ { n } ( s ) + s \ln s \tilde { F } _ { n } ^ { \ \prime } ( s )$$

$$+ \sum _ { m = 0 } ^ { \infty } \frac { ( 1 + \alpha _ { n } ) _ { m } ( 1 - \alpha _ { n } ) _ { m } } { ( 2 ) _ { m } m! } ( m + 1 ) \, s ^ { m } S _ { n, m } \Big ],$$

where the following notation has been introduced for simplicity

$$\tilde { F } _ { n } ( s ) \equiv \, _ { 2 } F _ { 1 } ( 1 + \alpha _ { n }, 1 - \alpha _ { n }, 2, s ),$$

$$\tilde { F } _ { n } ^ { \ \prime } ( s ) \equiv \frac { ( 1 - \alpha _ { n } ^ { 2 } ) } { 2 } \ _ { 2 } F _ { 1 } ( 2 + \alpha _ { n }, 2 - \alpha _ { n }, 3, s ).$$

The behavior far from the horizon r &gt;&gt; r h or s → 1 is not easily seen from (3.54)-(3.55) or (3.57)(3.58). For instance, factor (1 -s ) 3 2 / in (3.57)goes to zero at the infinite but the combination inside the parenthesis is divergent, so there is a 0 . ∞ ambiguity. But the result of this limit can be inferred from (3.49)-(3.50), since h n and k n are linear combinations of the functions f n or g n , which are centered at the infinite. Since g n ∼ -log( r ) is divergent at the infinite it follows that

$$h _ { n } \sim \log ( r ) \to \infty, \quad r \to \infty,$$

$$k _ { n } \sim \log ( r ) \to \infty, \quad r \to \infty.$$

The behavior from the derivatives follows from (3.43)-(3.44), since k n and h n are linear combinations of f n and g n and the last two have finite derivatives at the asymptotic region. From this simple fact it follows that

$$\frac { d } { d r } h _ { n } \sim \frac { 1 } { r } \to 0, \quad r \to \infty,$$

$$\frac { d } { d r } k _ { n } \sim \frac { 1 } { r } \to 0 \quad r \to \infty.$$

On the other hand, the behavior of the new solutions (3.54) and (3.55) at the horizon s → 0 is directly seen from its definitions, it is given by

$$h _ { n } \sim s F ( s ) \to 0, \quad s \to 0$$

$$k _ { n } \to s \log ( s ) F ( s ) - \frac { 1 } { \alpha _ { n } ^ { 2 } } \to - \frac { 1 } { \alpha _ { n } ^ { 2 } }, \quad s \to 0.$$

The behavior of their derivatives is

$$\frac { d } { d r } h _ { n } \sim \frac { 2 } { r _ { h } } ( F ( s ) + s F ^ { \prime } ( s ) ) \to \frac { 2 } { r _ { h } }, \quad s \to 0$$

$$\frac { d } { d r } k _ { n } \sim \frac { 2 } { r _ { h } } [ ( 1 + \log s ) F ( s ) + s \log s F ^ { \prime } ( s ) + s ] \to - \infty, \quad s \to 0.$$

Therefore we have reached our goal namely, to find a basis for which one of the eigenfunctions behave regularly at the horizon. It will be more convenient for our purposes to work with solutions satisfying (3.65)-(3.64) than using the ones satisfying (3.38)-(3.48). For this reason the following calculations will be referred to the set constituted by h n and k n .

## 3.4 The unphysical solutions

After elucidating the behavior of the solutions of the equation (3.20) for the potential A t , the next step is to discuss the boundary conditions for the electrostatic problem. The particularities of the BTZ geometry discussed in previous sections make the analysis different than in ordinary electrostatic in flat spaces, since the geometry is not asymptotically flat. In addition the behavior of the distance d r, r ( +1) given in (2.15) does not hold in a flat geometry. As a consequence of this behavior, when the usual boundary conditions of electrostatic in flat space are applied to this case, the electric field is not uniquely defined. This is an artifact which suggest that the new types of boundary conditions should be considered for the non asymptotically flat geometry.

The problems described above can be illustrated with an heuristic argument as follows. Consider a perfect dipole in flat space R 3 , constituted by two charges q and -q separated by a distance d . An elementary result in electrostatic states that the dipolar momenta p of such configuration is independent on the origin O of the coordinates. Thus this dipolar moment is the same near the origin or far away from it. This situation is radically different in the BTZ geometry. As it was discussed in (2.15), the distance d r ( ) between two points lying and on circle of radius r and r +1 and on the same line θ = θ 0 tends to zero when r → ∞ . Consider now two charges q and -q located at these points. If these charges are translated simultaneously in radial direction to r →∞ , since their mutual distance d r ( ) → 0, these charges become superposed one onto the other. It may then seem plausible that all the multipolar momenta tend to zero in this limit. The same reasoning holds for radially directed finite charged lines with total charge equal to zero.

The discussion given above suggest that one can send to the asymptotic region any finite number of radially directed neutral configurations, which will seem to disappear at the infinite. But if an infinite number of configurations is sent, then the result is ambiguous, since the resulting multipoles are an indetermination of the form 0 . ∞ .

These facts can be visualized by considering the multipole expansion of these radially directed neutral configurations. This expansion is expected to be of the form

$$A ( x, x ^ { \prime } ) = \sum _ { j = 1 } ^ { \infty } \frac { M _ { i _ { 1 } \dots i _ { j } } ^ { j } ( x ) \sigma ^ { i _ { 1 } } ( x, x ^ { \prime } ) \dots \sigma ^ { i _ { j } } ( x, x ^ { \prime } ) } { \sigma ^ { n _ { j } } ( x, x ^ { \prime } ) },$$

where σ i ( x, x ′ ) is a function with reduce to the usual difference σ i = ( x -x ′ ) i in a flat space. In addition M i 1 ...i j ( x ) are by definition multipolar momenta, and n j are positive numbers, whose explicit value is of no importance in this discussion. This is the essence of the Synge calculus [27], to be discussed in detail in the next sections. Here x is some characteristic point of the charged body and x ′ the observation point. The sum starts at j = 1 since the zero multipole, which is the total charge, is assumed to be zero.

Now, if the expansion (3.69) is applied to the BTZ case, one encounters the following ambiguity. When a neutral radially directed configuration whose center is at r 0 is sent to the asymptotic region, then r 0 take large values and the denominator tends to zero for in an small neighborhood. This follows

from the fact that σ a ( r , r 0 ) → 0 when r →∞ , r 0 →∞ such that r -r 0 &lt; ∞ . On the other hand, when r 0 →∞ the size of the system tends to zero, this follows from the behavior of d r ( ) in a BTZ geometry. It is plausible then that the multipoles M j i 1 ...i j are also zero in this limit, but this affirmation is to be taken with care. For instance, one can consider a dipole composed by two charges q and -q , which is sent to the asymptotic region while adding opposite charges at increasing positions, in such a way that when the dipole is centered at r the charges are q r ( ) and -q r ( ), with q r ( ) an arbitrary function of r . This function may be fixed to give a non zero value for the multipoles M j i 1 ...i j at infinite. In any case, if the multipoles are zero then there is an indetermination of the form 0 . ∞ for the potential (3.69). If instead the multipoles are finite or even infinite, then a divergence located at the infinite appears, which may give as a result a finite remanent electric field at finite r values. These arguments are of course heuristic, but suggest that the appropriate boundary conditions are not that straightforward as in the flat case.

The problems discussed in the previous paragraphs are reflected in the calculations as follows. The separation of variables for A t shows that the general electrostatic potential outside the source in a non rotating BTZ black hole admits an expansion of the form

$$A _ { t } ( r, \theta ) = \sum _ { n = 1 } ^ { \infty } \exp ( i n \theta ) ( A _ { n } h _ { n } ( r ) + B _ { n } k _ { n } ( r ) ) + \alpha \log ( r ) + \beta.$$

In general grounds one expect the electric field E i to vanish asymptotically and to be finite at the horizon. This field is obtained by taking derivatives of A t . More precisely, one expect the invariant

$$F _ { \mu \nu } F ^ { \mu \nu } = - ( \partial _ { r } A _ { t } ) ^ { 2 } - \frac { 1 } { r ^ { 2 } ( \frac { r ^ { 2 } } { l ^ { 2 } } - M ^ { 2 } ) } ( \partial _ { \theta } A _ { t } ) ^ { 2 },$$

to vanish at the infinite and to be finite at the horizon. Consider the simplest configuration first, namely, the one without charges. In this case B n should be zero since the derivatives of k n ( r ) are infinite at the horizon by (3.68). Thus, if B n were not zero, then the first term in (3.71) would be divergent. On the other hand the derivatives of h n ( r ) are well behaved at both the horizon and the asymptotic region. Nevertheless, its value is divergent at the infinite and since the derivative ∂ A θ t contains terms proportional to h n it follows that ∂ A θ t → ∞ at the asymptotic region. But this derivative is divided in (3.71) by a factor which diverges when r →∞ faster than h ′ n ( r ). In fact from (3.61) it follows that

$$\omega \varpi \ v a \varpi \quad & \frac { 1 } { r ^ { 2 } \left ( \frac { r ^ { 2 } } { l ^ { 2 } } - M \right ) } ( \partial _ { \theta } A _ { t } ) ^ { 2 } \sim \frac { 1 } { r ^ { 4 } } \log ^ { 2 } ( r ) \to 0, \quad r \to \infty.$$

So the invariant tends to zero at the infinite and so the electric field. The other term to be careful with is the denominator in the second term in (3.71), which gives a potentially divergence at the horizon. But taking into account (3.65) it follows that

$$\frac { 1 } { r ^ { 2 } ( \frac { r ^ { 2 } } { l ^ { 2 } } - M ) } ( \partial _ { \theta } A _ { t } ) ^ { 2 } \sim \left ( \frac { 1 } { r ^ { 4 } s } \right ) s ^ { 2 } F ^ { 2 } ( s ) \to 0, \quad r \to r _ { h }.$$

Thus, the presence of h n is not dangerous at the horizon either. The Gauss law fixes α = 0. Therefore

it is concluded that in absence of charges the most general potential is

$$A _ { t } ( r, \theta ) = \sum _ { n = 1 } ^ { \infty } A _ { n } \exp ( i n \theta ) h _ { n } ( r ),$$

with A n a arbitrary coefficients.

At first sight, this result may lead to the awkward conclusion that there exist an electric field, corresponding to (3.74), in absence of charge. The interpretation to be adopted in this work is that this conclusion is not true but instead, the solution (3.74) is unphysical , and corresponds to the electrostatic potential of 'configurations in the infinite' of the type mentioned above. These configurations are characteristic in a BTZ geometry due to the pathological behavior of the radial distance d r, r ( +1) explained in (2.15). Therefore, in an electrostatic problem in BTZ geometry our criteria for discarding solutions will not be the request that the radial solutions goes to zero at infinite, which is customary in ordinary electrodynamics in flat space. The task of determining appropriate boundary conditions is described in the next section.

## 4. Electrostatic field of BTZ black hole and wormhole

Having derived the eigenfunctions of the electrostatic problem in the BTZ geometry, the next step is the calculation of the electrostatic potential of a point charge q in front of the BTZ black hole and wormhole. The static charge is located at the position r 0 &gt; r h , θ = 0. Its electrostatic potential be expressed as

$$\text{seu as} A _ { t } ^ { I } ( r, \theta ) = \sum _ { n = 1 } ^ { \infty } \exp ( i n \theta ) ( A _ { n } h _ { n } ( r ) + B _ { n } k _ { n } ( r ) ) + \alpha \log ( r ) + \beta,$$

$$A _ { t } ^ { I I } ( r, \theta ) = \sum _ { n = 1 } ^ { \infty } \exp ( i n \theta ) ( C _ { n } h _ { n } ( r ) + D _ { n } k _ { n } ( r ) ) + \gamma \log ( r ) + \delta.$$

The potential A I t is the one in the region between the charge position r 0 and the horizon r h , and the A II t corresponds to the region between r 0 and the asymptotic boundary. In the first region B µ = 0 should be imposed, this is due to the fact that the derivatives of k µ ( r ) and the invariant (3.71) would not be bounded at the horizon. For the second region, the discussion below (3.74) suggest that if the coefficients C n multiplying h n ( r ) are non vanishing, then non trivial charge configurations at the infinite are turned on. This lead us to the following.

First type of boundary conditions: In order to avoid the residual electric field above one may impose that C n = 0 in the region between the charge and the asymptotic boundary. Note that if this boundary conditions is imposed, then (3.74) is automatically zero. This is expected by intuition namely, it implies the absence of electric field in absence of charges. The resulting potential for the charge has now the form

$$A _ { t } ^ { I } ( r, \theta ) = \sum _ { n = 1 } ^ { \infty } \exp ( i n \theta ) A _ { n } h _ { n } ( r ) + \alpha \log ( r ) + \beta,$$

$$A _ { t } ^ { I I } ( r, \theta ) = \sum _ { n = 1 } ^ { \infty } \exp ( i n \theta ) D _ { n } k _ { n } ( r ) + \gamma \log ( r ) + \delta.$$

However, there is some unpleasant detail concerning this choice of boundary conditions. First, it is not the only type of conditions that insures vanishing electric field in absence of charges. One may add to (4.77)-(4.78) a solution of the homogeneous Maxwell equation with the form (3.74), and with the coefficients A n proportional to the charge q , which will vanish when q → 0. This shows that there is an ambiguity for the choice of the boundary conditions.

Second type of boundary conditions: There exist an unique privileged type of boundary condition, which is seen as follows. The derivative of the functions k n ( r ), as shown in (3.64), decays as 1 /r at the asymptotic region. Instead, the derivatives of f n ( r ), as seen by (3.43), decay as 1 /r 3 . In view of this, it may look strange that the asymptotic behavior of (4.78) goes like 1 /r and not like 1 /r 3 . In other words, it may be reasonable to expect that, for a localized system of charges, the potential decays as fast as possible at the asymptotic region. This is of course what happens in an ordinary problem in electrostatics. And the boundary conditions imposing C n = 0 and leaving to the solution (4.77)-(4.78) do not respect this behavior. Therefore, one may consider the possibility of working with the base ( h , f n n ) instead of ( h , k n n ). In these terms the electrostatic potential in absence of charges is

$$A _ { t } ^ { I I } ( r, \theta ) = \sum _ { n = 1 } ^ { \infty } \exp ( i n \theta ) ( \alpha _ { n } h _ { n } ( r ) + \beta _ { n } f _ { n } ( r ) ) + \gamma \log ( r ) + \delta.$$

The requirement fast decay at the asymptotic implies α n = α = 0, otherwise the derivative of the potential would decay as 1 /r instead as 1 /r 3 . The regularity at the horizon requires β n = 0 since the derivatives of f n ( r ) are not bounded at the horizon. Thus these boundary conditions gives a constant potential δ and no electric field in absence of charges as well.

It should be remarked there are many boundary conditions giving no electric field when no charges are present, but only the second type one gives a decay of the form 1 /r 3 . This follows from the fact that the radial function f n ( r ) is the unique between the four f n ( r ), g n ( r ), h n ( r ) and k n ( r ) which this fast decay and therefore, the addition of other of these functions in the region II will spoil this behavior. In any case, it may be instructive to consider both types of conditions separately, and this will be done in the following subsection.

## 4.0.1 The black hole electrostatics for the first type of boundary conditions

Consider the conditions requiring C n = 0 first. We call this 'wrong' conditions since they are non unique, although this name will be justified better when studying the charge self-force. The matching conditions for the coefficients A n and D n are the request of continuity of the potential and the request of discontinuity of the electric field when crossing the surface r 0 along the radial line where the charge is located. These conditions are translated into the following linear equations for the unknown

coefficients

$$A _ { n } h _ { n } ( r _ { 0 } ) = D _ { n } k _ { n } ( r _ { 0 } ), \\ A _ { n } h _ { n } ^ { \prime } ( r _ { 0 } ) - D _ { n } k _ { n } ^ { \prime } ( r _ { 0 } ) = \frac { q } { r _ { 0 } }, \\ = - q, \quad \alpha = 0, \quad \beta = - q \log ( r _ { 0 } ) + \delta.$$

$$\gamma = - q, \quad \alpha = 0, \quad \beta = - q \log ( r _ { 0 } ) + \delta.$$

The constant δ can be fixed to zero without losing generality. The solution of this system is

$$A _ { n } = \frac { q k _ { n } ( r _ { 0 } ) } { r _ { 0 } W _ { n } ( r _ { 0 } ) }, \quad D _ { n } = \frac { q h _ { n } ( r _ { 0 } ) } { r _ { 0 } W _ { n } ( r _ { 0 } ) },$$

with

$$W ( r _ { 0 } ) = k _ { n } ( r _ { 0 } ) h _ { n } ^ { \prime } ( r _ { 0 } ) - k _ { n } ^ { \prime } ( r _ { 0 } ) h _ { n } ( r _ { 0 } ),$$

the Wronskian of the two solutions k n and h n at the charge radial position r 0 . Therefore the electrostatic potential is given by

$$A _ { t } ^ { I } ( r, \theta ) = - q \log ( r _ { 0 } ) + \sum _ { n = 0 } ^ { \infty } \frac { q k _ { n } ( r _ { 0 } ) h _ { n } ( r ) } { r _ { 0 } W _ { n } ( r _ { 0 } ) } \exp ( i n \theta ),$$

$$A _ { t } ^ { I I } ( r, \theta ) = - q \log ( r ) + \sum _ { n = 0 } ^ { \infty } \frac { q h _ { n } ( r _ { 0 } ) k _ { n } ( r ) } { r _ { 0 } W _ { n } ( r _ { 0 } ) } \exp ( i n \theta ).$$

The Wronskian (4.80) can be calculated as follows. Consider two arbitrary linearly independent solutions y 1 n and y 2 n of the hypergeometric equation

$$u ( 1 - u ) y _ { n } ^ { \prime \prime i } ( u ) + ( 1 - u ) y _ { n } ^ { \prime i } ( u ) - \frac { n ^ { 2 } y _ { n } ^ { i } ( u ) } { 4 M } = 0.$$

By multiplying the equation for y 1 n by y 2 n and by doing the opposite procedure for the equation for y 2 n , then after subtracting the resulting terms it is obtained the following equation

$$u \, \widetilde { W } ^ { \prime } ( u ) + \widetilde { W } ( u ) = 0,$$

for the wronskian

$$\widetilde { W } ( u ) = y _ { n } ^ { 1 } y _ { n } ^ { \prime 2 } - y _ { n } ^ { \prime 1 } y _ { n } ^ { 2 }.$$

If the wronskian ˜ W v ( ) at a point v is known, then the solution of the last equation is

$$\widetilde { W } ( u ) = \frac { v \widetilde { W } ( v ) } { u }.$$

The Wronskian just considered is referred to derivatives in u . The wronskian W r ( ) referred to derivatives of r is obtained by multiplying ˜ W by u ′ ( r ), the result is

$$W ( r ) = \frac { r ^ { \prime } W ( r ^ { \prime } ) } { r },$$

with r ′ an arbitrary fixed point. Therefore, once the wronksian at a given point r ′ is known, its values at a generic point r are determined by the last formula. For the case in consideration, it is convenient

to calculate W h ,k ( n n ) at r h , which corresponds to s = 0. The value follows directly from (3.65)-(3.68), the result is

$$W _ { n } ( r _ { h } ) = \frac { 2 F _ { n } ( 0 ) } { r _ { h } \alpha _ { n } ^ { 2 } }, \\ \omega \dots \quad \omega \dots \cdots \cdots \cdots$$

and taking into account the definition of the hypergeometric function

$$_ { 2 } F _ { 1 } ( \alpha, \beta ; \gamma ; s ) = \sum _ { n = 0 } ^ { \infty } \frac { ( \alpha ) _ { n } ( \beta ) _ { n } } { ( n ) ( \gamma ) _ { n } } s ^ { n },$$

it is concluded that

$$_ { 2 } F _ { 1 } ( \alpha, \beta ; \gamma ; 0 ) = 1.$$

By this and (3.59) the wronskian (4.84) takes the following form

$$W _ { n } ( r _ { h } ) = \frac { 8 M } { r _ { h } n ^ { 2 } }.$$

In these terms (4.81)-(4.82) become

$$A _ { t } ^ { I } ( r, \theta ) = - q \log ( r _ { 0 } ) + \sum _ { n = 1 } ^ { \infty } \frac { q n ^ { 2 } k _ { n } ( r _ { 0 } ) h _ { n } ( r ) } { 8 M$$

$$A _ { t } ^ { I I } ( r, \theta ) = - q \log ( r ) + \sum _ { n = 1 } ^ { \infty } \frac { q n ^ { 2 } h _ { n } ( r _ { 0 } ) k _ { n } ( r ) } { 8 M } \exp$$

This is the electrostatic potential corresponding for the first type of boundary conditions.

## 4.0.2 The black hole electrostatics for the second type of boundary conditions

The second boundary conditions, which will be called of 'right' type state that the region II should be described in terms of the fast decaying radial functions f n ( r ), while the region I should be described by h n ( r ), which are regular at the horizon. The electrostatic potential satisfying this requirement is generically

$$A _ { t } ^ { I } ( r, \theta ) = \sum _ { n = 1 } ^ { \infty } \exp ( i n \theta ) A _ { n } h _ { n } ( r ) + \alpha \log ( r ) + \beta,$$

$$A _ { t } ^ { I I } ( r, \theta ) = \sum _ { n = 1 } ^ { \dots } \exp ( i n \theta ) D _ { n } f _ { n } ( r ) + \gamma \log ( r ) + \delta. \quad \quad \quad$$

The analysis of the boundary conditions at the charge position r 0 is completely analogous to the one performed in the previous section. For the case M = 1 the Wronskian constructed in terms of ( h n ( r , f ) n ( r )) is given by r W

(

r

) =

2

0

n

0

Γ(1 +

in

.

2

)Γ(1

-

in

2

)

The electrostatic potential in this case is

$$A _ { t } ^ { I } ( r, \theta ) = - q \log ( r _ { 0 } ) + \frac { q } { 2 } \sum _ { n = 1 } ^ { \infty } \Gamma ( 1 + \frac { i n } { 2 } ) \Gamma ( 1 - \frac { i n } { 2 } ) f _ { n } ( r _ { 0 } ) h _ { n } ( r ) \exp ( i n \theta ),$$

$$A _ { t } ^ { I I } ( r, \theta ) = - q \log ( r ) + \frac { q } { 2 } \sum _ { n = 1 } ^ { \infty } \Gamma ( 1 + \frac { i n } { 2 } ) \Gamma ( 1 - \frac { i n } { 2 } ) h _ { n } ( r _ { 0 } ) f _ { n } ( r ) \exp ( i n \theta ). \quad \ \ ( 4. 9 1 )$$

This is the unique potential with the right discontinuity at the charge position, the right behavior at the horizon and and with the fastest decaying conditions.

## 4.1 The wormhole electrostatic field

In order to find the electrostatic potential of a charge in front of a BTZ wormhole it is convenient to divide the space time in the following three regions

$$\text{Region I} \colon \quad r _ { - } < r _ { g },$$

$$\text{Region II} \colon \quad r _ { g } < r _ { + } < r _ { 0 },$$

Region III: r + 0

$$& - c r _ { g }, \\ & < r _ { + } < r _ { 0 }, \\ & r _ { + } < r _ { 0 }.$$

Here r g indicates the throat position. The electrostatic solution in any of these regions is of the form

$$g \text{ indicates the throat position. }'I \text{ the electrostatic solution in any of these regions is of the form} \\ A _ { n } ^ { I I } & = - q \log ( r ) + \sum _ { n = 1 } ^ { \infty } ( A _ { n } h _ { n } ( r ) + B _ { n } k _ { n } ( r ) ) \exp [ n ( \theta - \theta _ { 0 } ) ], \\ A _ { n } ^ { I I } & = - q \log ( r _ { 0 } ) + \sum _ { n = 1 } ^ { \infty } ( C _ { n } h _ { n } ( r ) + D _ { n } k _ { n } ( r ) ) \exp [ n ( \theta - \theta _ { 0 } ) ], \\ A _ { n } ^ { I } & = - q \log ( r _ { 0 } ) + \sum _ { n = 1 } ^ { \infty } ( E _ { n } h _ { n } ( r ) + F _ { n } k _ { n } ( r ) ) \exp [ n ( \theta - \theta _ { 0 } ) ]. \\ \text{efficients } A _ { n },.., F _ { n } \text{ are given by the boundary conditions of the problem, which are the following.}$$

The coefficients A ,.., F n n are given by the boundary conditions of the problem, which are the following.

- 1. The potential is continuous in r + = r 0 ,

$$A ^ { I I I } ( r _ { + } \to r _ { 0 } ^ { + } ) = A ^ { I I } ( r _ { + } \to r _ { 0 } ^ { - } )$$

- 2. The potential is continuous in r -= r + = r g

$$A ^ { I I } ( r _ { + } \to r _ { g } ^ { + } ) = A ^ { I } ( r _ { - } \to r _ { g } ^ { - } )$$

- 3. Continuity of the field in r -= r + = r g

$$\partial _ { r _ { + } } A ^ { I I } ( r _ { + } \to r _ { g } ^ { + } ) = - \partial _ { r _ { - } } A ^ { I } ( r _ { - } \to r _ { g } ^ { - } )$$

- 4. Discontinuity of the electric field in r + = r 0

$$\partial _ { r _ { + } } A ^ { I I I } ( r _ { + } \to r _ { 0 } ^ { + } ) - \partial _ { r _ { + } } A ^ { I I } ( r _ { + } \to r _ { 0 } ^ { - } ) = - \frac { 2 \pi q } { r _ { 0 } } \delta ( \theta - \theta _ { 0 } )$$

- 5. The correct asymptotic behavior in M -and M + .

The last condition is related to the two boundary conditions discussed in the previous section for the black hole case. For the first type of boundary conditions

$$B _ { n } = F _ { n } = 0.$$

The other four conditions imply that

$$A _ { n } h _ { n } ( r _ { 0 } ) = C _ { n } h _ { n } ( r _ { 0 } ) + D _ { n } k _ { n } ( r _ { 0 } ),$$

$$q$$

$$A _ { n } h _ { n } ^ { \prime } ( r _ { 0 } ) = C _ { n } h _ { n } ^ { \prime } ( r _ { 0 } ) + D _ { n } k _ { n } ^ { \prime } ( r _ { 0 } ) + \frac { q } { r _ { 0 } }, \\ E _ { n } h _ { n } ( r _ { g } ) = C _ { n } h _ { n } ( r _ { g } ) + D _ { n } k _ { n } ( r _ { g } ), \\ - E _ { n } h _ { n } ^ { \prime } ( r _ { g } ) = C _ { n } h _ { n } ^ { \prime } ( r _ { g } ) + D _ { n } k _ { n } ^ { \prime } ( r _ { g } ).$$

This is system of four equations with four undetermined, whose solution gives

$$A _ { t } ^ { I I I } = - q \log ( r ) + \sum _ { n = 1 } ^ { \infty } \left [ \frac { q k _ { n } ( r _ { 0 } ) } { r _ { 0 } W _ { n } ( r _ { 0 } ) } - \frac { q h _ { n } ( r _ { 0 } ) } { 2 r _ { 0 } h _ { n } ( r _ { g } ) h _ { n } ^ { \prime } ( r _ { g } ) } \left [ \frac { W _ { n } ( r _ { g } ) + 2 h _ { n } ( r _ { g } ) k _ { n } ^ { \prime } ( r _ { g } ) } { W _ { n } ( r _ { 0 } ) } \right ] \right ] h _ { n } ( r ) \exp [ n ( \theta - \theta _ { 0 } ) ].$$

Attention will be paid only for the potential in the third region, since is the one to be used for calculating the charge self-force. It can be decomposed further as

$$A _ { t } ^ { I I I } = A _ { t } ^ { b h } - \frac { 1 } { 2 } \sum _ { n = 1 } ^ { \infty } \left [ \frac { q h _ { n } ( r _ { 0 } ) } { r _ { g } h _ { n } ( r _ { g } ) h _ { n } ^ { \prime } ( r _ { g } ) } + \frac { q n ^ { 2 } h _ { n } ( r _ { 0 } ) k _ { n } ^ { \prime } ( r _ { g } ) } { 8 M h _ { n } ^ { \prime } ( r _ { g } ) } \right ] h _ { n } ( r ) \exp [ n ( \theta - \theta _ { 0 } ) ],$$

with A bh t the potential corresponding to the black hole solution (4.86)-(4.87).

For the second type of boundary conditions the resulting potential is

$$A _ { t } ^ { I I I } = A _ { t } ^ { b h } - \frac { 1 } { 2 } \sum _ { n = 1 } ^ { \infty } \left [ \frac { q h _ { n } ( r _ { 0 } ) } { r _ { g } h _ { n } ( r _ { g } ) h _ { n } ^ { \prime } ( r _ { g } ) } + \Gamma ( 1 + \frac { i n } { 2 } ) \Gamma ( 1 - \frac { i n } { 2 } ) \frac { q h _ { n } ( r _ { 0 } ) k _ { n } ^ { \prime } ( r _ { g } ) } { 2 h _ { n } ^ { \prime } ( r _ { g } ) } \right ] h _ { n } ( r ) \exp [ n ( \theta - \theta _ { 0 } ) ], \ ( 4. 1 0 0 )$$

where now A bh t is the potential corresponding to the black hole solution (4.90)-(4.91). The remaining sum is due to the effect of the throat at r g , which deform the electric field lines. This shows that both geometries, which are locally the same, can be distinguished by electrostatic effects.

## 5. Coincident points limits and Taylor like expansions in curved space times

The electrostatic potential A t for the static charge q in any geometry is singular at the position where the charge is located. In a flat space, this charge does not experience any self-force, this is clear due to the rotational symmetry of the electrostatic field. In a curved space, this argument is not true, since the non trivial curvature of the geometry deforms the electric lines and gives a net force on the charge. A seminal work about electrostatic in curved space is the one of Haddamard [35], who started a program for calculating the singular part for A t in static geometries.

The electrostatic vector potential A x,x t ( ′ ) is an example of a bivector, since it depends on two arguments, the position of the charge x and the position of the observer x ′ . The self force on the charge is determined by the behavior of A x,x t ( ′ ) in an infinitesimal neighborhood of x , and the analogous of a Taylor expansion in a curved space plays an important role in determining this behavior. In the present section the main properties of these expansions are described, which requires Synge calculus [27]. The references [28]- [36] are more detailed and contains more information. Nevertheless a concise but self-contained description of the Synge calculus is given in the following subsections.

## 5.1 The Synge world function and its main properties

Our task is to calculate the self force of a static charge in front of a BTZ black hole. This requires to calculate its electrostatic field E and to substract the part that it is divergent at the position of the charge. There are several methods to extract this singular part. The one to be implemented here, which is better adapted to static geometries, has as a basic ingredients the parallel propagator bitensor g β α ( x, x ′ ) and the Synge world function σ x, x ( ′ ) [27]. To define them, consider a space time ( g , M ) and choose an einbein basis e a for the metric g such that

$$g _ { \mu \nu } ( x ) = \eta _ { a b } e _ { \mu } ^ { a } ( x ) e _ { \nu } ^ { b } ( x ).$$

There is an SO n ( -1 1) freedom for choosing this basis, since , SO n ( -1 1) rotations , R a b induce new one forms

$$e ^ { \prime a } ( x ) = R _ { b } ^ { a } e ^ { b } ( x )$$

which are still an einbein for the metric g µν . In particular, since always ∇ g ab = 0 one may choose an einbein e a ( x ) at x such that for any x ′ lying in the injectivity radius of x it is parallel transported along the unique geodesic γ joining the two points.

Now given a vector field A µ ( x ) defined at TM x one can express it in the basis e a ( x ) as

$$A _ { \mu } ( x ) = A _ { a } e _ { \mu } ^ { a } ( x ).$$

If this vector is parallel transported to x ′ along γ then its components at that point are

$$A _ { \mu } ( x ^ { \prime } ) = A _ { a } e _ { \mu } ^ { a } ( x ^ { \prime } ),$$

and it follows that with

$$g _ { \mu } ^ { \nu } ( x, x ^ { \prime } ) = e _ { \mu } ^ { a } ( x ) e _ { a } ^ { \nu } ( x ^ { \prime } ).$$

The object (5.101) then relates the components of the vector field A µ at x and x ′ . This object is by definition the parallel transport bitensor of the geometry.

Let us turn the attention to the Synge world function σ x, x ( ′ ). This function is defined as half of the square of the geodesic distance d x, x ( ′ ) between x and x ′

$$\sigma ( x, x ^ { \prime } ) = \frac { 1 } { 2 } d ( x, x ^ { \prime } ) ^ { 2 }.$$

This distance d x, x ( ′ ) can be represented in integral form as

$$d ( x, x ^ { \prime } ) = \int _ { 0 } ^ { 1 } \sqrt { g _ { a b } \dot { x } ^ { a } \dot { x } ^ { b } } d \lambda, \quad x ( 0 ) = x, \quad x ( 1 ) = x ^ { \prime },$$

with ˙ x a satisfying the geodesic equation

$$\dot { x } ^ { a } \nabla _ { a } \dot { x } ^ { b } = 0.$$

$$A _ { \mu } ( x ^ { \prime } ) = g _ { \mu } ^ { \nu } ( x, x ^ { \prime } ) A _ { \nu } ( x ),$$

The bivector

$$n _ { a } = \nabla _ { a } d ( x, x ^ { \prime } ), \quad n _ { a ^ { \prime } } = \nabla _ { a ^ { \prime } } d ( x, x ^ { \prime } ),$$

constructed by taking derivatives of the distance d x, x ( ′ ) with respect to the initial or final point x or x ′ , has unit length. This can be seen explicitly by calculating the variation of the distance

$$\delta d = d ( x + \delta x, x ^ { \prime } ) - d ( x, x ^ { \prime } ),$$

with δx subject to the boundary conditions

$$\delta x ( 0 ) = \delta x _ { 0 }, \quad \delta x ( 1 ) = 0.$$

Now the integral (5.103) represents the distance d x, x ( ′ ) as an action with lagrangian

$$\mathcal { L } & = \sqrt { g _ { a b } \dot { x } ^ { a } \dot { x } ^ { b } },$$

and the standard theory of Hamilton-Jacobi implies that the last variation is

$$\delta d & = n _ { a } \delta x ^ { a }, &$$

with n a the momentum corresponding to the coordinate x a calculated with the lagrangian (5.106), which is given by

$$n _ { a } = - \frac { \partial \mathcal { L } } { \partial \dot { x } _ { a } } = \frac { g _ { a b } ( x ) \dot { x } ^ { b } } { \sqrt { g _ { a b } \dot { x } ^ { a } \dot { x } ^ { b } } }.$$

Clearly, this is the bivector (5.105) and it follows from the last expression that

$$g ^ { a b } n _ { a } n _ { b } = 1,$$

which proves that n a has unit length, as stated. In addition

$$n _ { a ^ { \prime } } ( x, x ^ { \prime } ) = - g _ { a ^ { \prime } } ^ { a } ( x, x ^ { \prime } ) n _ { a } ( x, x ^ { \prime } ),$$

which follows from the definition of the parallel transport bitensor. Furthermore, the norm g ab ˙ x x a ˙ b is constant along a geodesic γ , this follows from the Levi-Civita condition ∇ c g ab = 0 together with the geodesic equation (5.104). This implies that the Synge function can be expressed in integral form as

$$\sigma ( x, x ^ { \prime } ) = \frac { 1 } { 2 } \int _ { 0 } ^ { 1 } g _ { a b } \dot { x } ^ { a } \dot { x } ^ { b } d \lambda, \quad x ( 0 ) = x, \quad x ( 1 ) = x ^ { \prime },$$

and their derivatives σ a = ∇ a σ also satisfy some useful identities analogous to (5.108). One of them is

$$g ^ { a b } \sigma _ { a } \sigma _ { b } = g ^ { a ^ { \prime } b ^ { \prime } } \sigma _ { a ^ { \prime } } \sigma _ { b ^ { \prime } } = \sigma,$$

its proof follows directly from the action representation (5.110) and the Hamilton-Jacobi theory. Note that for the flat metric σ = η ab ( x -x ′ ) a ( x -x ′ ) b and the identity (5.111) is immediate. In addition the following relation takes place

$$g _ { a } ^ { a ^ { \prime } } ( x, x ^ { \prime } ) \sigma ^ { a } ( x, x ^ { \prime } ) = - \sigma ^ { a ^ { \prime } } ( x, x ^ { \prime } ),$$

which is the analogous of (5.109). In the following the notation σ i 1 ..i n = ∇ 1 ... ∇ n σ will be employed. With this notation the formula (5.111) can be differentiated with respect to the coordinate x giving that

$$\sigma _ { a } ^ { a ^ { \prime } } \sigma ^ { a } = \sigma ^ { a ^ { \prime } }.$$

At this point it is convenient to introduce more formally the definition of a bitensor, since this is a notion to be used recurrently in the following.

## 5.2 Taylor expansions of bitensors

Consider an arbitrary manifold M and choose two of its points x y x ′ . A bitensor T β ..β 1 k β ..β ′ 1 ′ l α ..α 1 m α ...α ′ 1 ′ n is a linear application of the form

$$T \colon T M _ { x } \times \dots \times T M _ { x } \times T M _ { x ^ { \prime } } \times \dots \times T M _ { x ^ { \prime } } \times T M _ { x } ^ { * } \times \dots \times T M _ { x ^ { \prime } } ^ { * } \to C$$

with TM p the tangent space at the point p and TM ∗ p its dual. In particular n = l = 0 corresponds to a tensor m times covariant and k times contravariant. A bitensor field is a rule that assigns to the pair of points ( x, x ′ ) ∈ M × M the bitensor T β ..β 1 k β ..β ′ 1 ′ l α ..α 1 m α ...α ′ 1 ′ n ( x, x ′ ). Such objets can be differentiated with respect to x or x ′ . The derivative ∇ a T β ..β 1 k β ..β ′ 1 ′ l α ..α 1 m α ...α ′ 1 ′ n it is obtained by considering x ′ frozen and taking its covariant derivative by pretending that it is an m times covariant and k times contravariant tensor. As a simple example consider the object σ x, x ( ′ ), which is a biscalar. The covariant derivative σ a ( x, x ′ ) is a bivector and by taking successive derivatives one can construct a bitensor of arbitrary order.

When the points of a generic bitensor of rank two T α β ′ ′ ( x, x ′ ) are close enough, one may make a Taylor like expansion of the form

$$T _ { \alpha ^ { \prime } \beta ^ { \prime } } ( x, x ^ { \prime } ) = A _ { \alpha ^ { \prime } \beta ^ { \prime } } ( x ^ { \prime } ) + A _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } ( x ^ { \prime } ) \sigma ^ { \gamma ^ { \prime } } ( x, x ^ { \prime } ) + A _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } \beta ^ { \prime } } ( x ^ { \prime } ) \sigma ^ { \gamma ^ { \prime } } ( x, x ^ { \prime } ) \sigma ^ { \beta ^ { \prime } } ( x, x ^ { \prime } ) + O ( \epsilon ^ { 3 } ), \quad ( 5. 1 1 4 )$$

with glyph[epsilon1] a characteristic value of σ α ′ and A i 1 ..i n ( x ′ ) are ordinary tensors defined at x ′ . This is the analogous of an ordinary Taylor expansion in a flat space. Assuming that T α β ′ ′ ( x, x ′ ) is known the task is to calculate the coefficients A i 1 ..i n ( x ′ ) of (5.114) . This requires the analysis of the coincident points limit in (5.114) and all its derivatives [27], [33]. Given an arbitrary bitensor U x,x ( ′ ) this limit is defined by the formula

$$[ U ] ( x ^ { \prime } ) = \lim _ { x \rightarrow x ^ { \prime } } U ( x, x ^ { \prime } ).$$

In these terms it follows directly that the first coefficient of the expansion (5.114) is simply

$$A _ { \alpha ^ { \prime } \beta ^ { \prime } } ( x ^ { \prime } ) = [ T _ { \alpha ^ { \prime } \beta ^ { \prime } } ].$$

The calculation of the higher order terms A i 1 ...i n in (5.114) requires the knowledge of the coincident point limits of the covariant derivatives σ i 1 ..i n . These limits can be calculated as follows. First it is true that

$$[ \sigma ] = 0,$$

since the distance between to points when x → x ′ goes to zero. This condition together with (5.111) imply that

$$[ \sigma _ { \alpha } ] & = 0.$$

Besides one has

$$[ \sigma _ { \alpha \beta } ] = [ \sigma _ { \alpha ^ { \prime } \beta ^ { \prime } } ] = g _ { \alpha \beta } ( x ^ { \prime } ), \quad [ \sigma _ { \alpha ^ { \prime } \beta } ] = [ \sigma _ { \alpha \beta ^ { \prime } } ] = - g _ { \alpha \beta } ( x ^ { \prime } ).$$

The last relations are intuitive by considering the flat case and can be established by use of the last two formulas together with (5.113). Now, In order to calculate the other higher order coincident limits, it is convenient to take two covariant derivatives in (5.111) to obtain that

$$\sigma _ { \alpha \beta \gamma } = \sigma _ { \alpha \beta \gamma } ^ { \delta } \sigma _ { \delta } + \sigma _ { \alpha \beta } ^ { \delta } \sigma _ { \delta \gamma } + \sigma _ { \alpha \gamma } ^ { \delta } \sigma _ { \delta \beta } + \sigma _ { \alpha } ^ { \delta } \sigma _ { \delta \beta \gamma }.$$

and taking into account (5.116)-(5.118), it follows that

$$[ \sigma _ { \alpha \beta \gamma } ] = [ \sigma _ { \alpha \beta } ^ { \delta } ] g _ { \gamma ^ { \prime } \delta ^ { \prime } } ( x ^ { \prime } ) + [ \sigma _ { \alpha \gamma } ^ { \delta } ] g _ { \delta ^ { \prime } \beta ^ { \prime } } ( x ^ { \prime } ) + [ \sigma _ { \delta \beta \gamma } ] \delta _ { \alpha ^ { \prime } } ^ { \delta ^ { \prime } }.$$

from where it is obtained that

$$[ \sigma _ { \gamma \alpha \beta } ] + [ \sigma _ { \beta \alpha \gamma } ] = 0.$$

Also, since σ is a biscalar, it follows that σ αβ = σ βα . By use of this and the Ricci identity it follows that

$$2 [ \sigma _ { \alpha \beta \gamma } ] = R ^ { \delta } _ { \alpha \beta \gamma } ( x ^ { \prime } ) [ \sigma _ { \delta } ],$$

This, together with (5.117) shows that

$$[ \sigma _ { \alpha \beta \gamma } ] & = 0.$$

Analogously, it can be shown that

$$[ \sigma _ { \alpha \beta \gamma ^ { \prime } } ] = [ \sigma _ { \alpha \beta ^ { \prime } \gamma ^ { \prime } } ] = [ \sigma _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } ] = 0.$$

To proceed further requires to take the covariant derivative of (5.119) to obtain

$$\sigma _ { \alpha \beta \gamma \delta } = \sigma _ { \alpha \beta \gamma \delta } ^ { \epsilon } \sigma _ { \epsilon } + \sigma _ { \alpha \beta \gamma } ^ { \epsilon } \sigma _ { \epsilon \delta } + \sigma _ { \alpha \beta \delta } ^ { \epsilon } \sigma _ { \epsilon \gamma } + \sigma _ { \alpha \gamma \delta } ^ { \epsilon } \sigma _ { \epsilon \beta } + \sigma _ { \alpha \beta } ^ { \epsilon } \sigma _ { \epsilon \gamma \delta }$$

$$+ \sigma _ { \alpha \gamma } ^ { \epsilon } \sigma _ { \epsilon \beta \delta } + \sigma _ { \alpha \delta } ^ { \epsilon } \sigma _ { \epsilon \beta \gamma } + \sigma _ { \alpha } ^ { \epsilon } \sigma _ { \epsilon \beta \gamma \delta }.$$

The limit of coincident points in (5.124) shows that

$$[ \sigma _ { \alpha \beta \gamma \delta } ] + [ \sigma _ { \alpha \delta \beta \gamma } ] + [ \sigma _ { \alpha \gamma \delta \beta } ] = 0.$$

The last expression can be worked further by taking the derivative of the Ricci identity

$$\sigma _ { \alpha \beta \gamma } = \sigma _ { \alpha \gamma \beta } - R ^ { \delta } _ { \alpha \beta \gamma } \sigma _ { \delta },$$

with respect to x glyph[epsilon1] and taking the coincident point limit. The result is

$$[ \sigma _ { \alpha \beta \gamma \delta } ] = [ \sigma _ { \alpha \gamma \beta \delta } ] + R _ { \alpha ^ { \prime } \delta ^ { \prime } \gamma ^ { \prime } \beta ^ { \prime } }$$

Besides, the Ricci identity implies that

$$\sigma _ { \alpha \beta \gamma \delta } = \sigma _ { \alpha \beta \delta \gamma } - R ^ { \epsilon } _ { \alpha \gamma \delta } \sigma _ { \epsilon \beta } - R ^ { \epsilon } _ { \beta \gamma \delta } \sigma _ { \epsilon \alpha }$$

and this, together with the symmetry properties of the curvature tensor gives

$$[ \sigma _ { \alpha \beta \gamma \delta } ] = [ \sigma _ { \alpha \beta \delta \gamma } ].$$

In these terms it follows that (5.125) leads to

$$[ \sigma _ { \alpha \beta \gamma \delta } ] = - \frac { 1 } { 3 } ( R _ { \alpha ^ { \prime } \gamma ^ { \prime } \beta ^ { \prime } \delta ^ { \prime } } + R _ { \alpha ^ { \prime } \delta ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } ).$$

In analogous way the following identities

$$[ \sigma _ { ; \alpha \beta \gamma \delta ^ { \prime } } ] = \frac { 1 } { 3 } ( R _ { \alpha ^ { \prime } \gamma ^ { \prime } \beta ^ { \prime } \delta ^ { \prime } } + R _ { \alpha ^ { \prime } \delta ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } ),$$

$$[ \sigma _ { ; \alpha \beta \gamma ^ { \prime } \delta ^ { \prime } } ] = - \frac { 1 } { 3 } ( R _ { \alpha ^ { \prime } \gamma ^ { \prime } \beta ^ { \prime } \delta ^ { \prime } } + R _ { \alpha ^ { \prime } \delta ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } ),$$

$$[ \sigma _ { ; \alpha \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } } ] = - \frac { 1 } { 3 } ( R _ { \alpha ^ { \prime } \gamma ^ { \prime } \beta ^ { \prime } \delta ^ { \prime } } + R _ { \alpha ^ { \prime } \delta ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } ),$$

$$[ \sigma _ { ; \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } } ] = - \frac { 1 } { 3 } ( R _ { \alpha ^ { \prime } \gamma ^ { \prime } \beta ^ { \prime } \delta ^ { \prime } } + R _ { \alpha ^ { \prime } \delta ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } ),$$

can be proven.

Once the limits (5.116)-(5.130) are known the coefficients A i 1 ..i n ( x ′ ) of (5.114) can be calculated to the third order. The coincident point limit in (5.114) and (5.116)-(5.130) give the following recurrence formula

$$A _ { \alpha ^ { \prime } \beta ^ { \prime } } & = [ T _ { \alpha ^ { \prime } \beta ^ { \prime } } ],$$

$$A _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } = [ T _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } } ] - A _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } },$$

$$A _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } } = [ T _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } \delta ^ { \prime } } ] - A _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } \delta ^ { \prime } } - A _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } ; \delta ^ { \prime } } - A _ { \alpha ^ { \prime } \beta ^ { \prime } \delta ^ { \prime } ; \gamma ^ { \prime } }.$$

If T αβ is known, these formulas allow to determine the expansion coefficients up to order three.

The expansion (5.131)-(5.133) is valid for a bitensor with indices referred to the point x ′ . Consider now the expansion of a bitensor of the form T α β ′ ( x, x ′ ). In this case one can construct an associated tensor ˜ T α β ′ ′ ( x, x ′ ) given by

$$\widetilde { T } _ { \alpha ^ { \prime } \beta ^ { \prime } } ( x, x ^ { \prime } ) = g _ { \beta ^ { \prime } } ^ { \beta } ( x, x ^ { \prime } ) T _ { \alpha ^ { \prime } \beta } ( x, x ^ { \prime } ),$$

which can be expanded by use of (5.131)-(5.133) and (5.114) as

$$\widetilde { T } _ { \alpha ^ { \prime } \beta } = B _ { \alpha ^ { \prime } \beta ^ { \prime } } + B _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } \sigma ^ { \gamma ^ { \prime } } + B _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } } \sigma ^ { \gamma ^ { \prime } } \sigma ^ { \delta ^ { \prime } } + O ( \epsilon ^ { 3 } ).$$

The formula (5.134) can be inverted

$$T _ { \alpha ^ { \prime } \beta } ( x, x ^ { \prime } ) = g _ { \beta } ^ { \beta ^ { \prime } } ( x, x ^ { \prime } ) \widetilde { T } _ { \alpha ^ { \prime } \beta ^ { \prime } } ( x, x ^ { \prime } ),$$

and this together with (5.135) gives that

$$T _ { \alpha ^ { \prime } \beta } = g _ { \beta } ^ { \beta ^ { \prime } } ( B _ { \alpha ^ { \prime } \beta ^ { \prime } } + B _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } \sigma ^ { \gamma ^ { \prime } } + B _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } } \sigma ^ { \gamma ^ { \prime } } \sigma ^ { \delta ^ { \prime } } ) + O ( \epsilon ^ { 3 } ).$$

The evaluation of the coefficients B i 1 ..i n of this expansion requires the use of (5.116)-(5.130) and also the coincident point limit of the parallel propagator g β ′ α ( x, x ′ ) and its derivative. These can be calculated as follows. First of all, it is evident from the definition that

$$[ g _ { \beta } ^ { \alpha ^ { \prime } } ] = \delta _ { \beta } ^ { \alpha ^ { \prime } }.$$

Besides, the parallel transport propagator can be constructed as g α ′ β = e α ′ a e a β with e a β an einbein basis which is parallel transported along a geodesic, which means that

$$e _ { a ; \beta } ^ { \alpha } \sigma ^ { \beta } = 0.$$

This implies that and the use of the Ricci identity gives

$$g _ { \beta ; \gamma } ^ { \alpha ^ { \prime } } \sigma ^ { \gamma } = 0.$$

A differentiation of the last formula gives

$$g ^ { \alpha ^ { \prime } } _ { \beta ; \gamma \delta } \sigma ^ { \gamma } + g ^ { \alpha ^ { \prime } } _ { \beta ; \gamma } \sigma ^ { \gamma } _ { \delta } = 0,$$

and taking the coincidence limit and using (5.116)-(5.118) give that

$$[ g _ { \beta ; \gamma } ^ { \alpha ^ { \prime } } ] = [ g _ { \beta ; \gamma ^ { \prime } } ^ { \alpha ^ { \prime } } ] = 0.$$

Further differentiation gives

$$g _ { \beta ; \gamma \delta \epsilon } ^ { \alpha ^ { \prime } } \sigma ^ { \gamma } + g _ { \beta ; \gamma \delta } ^ { \alpha ^ { \prime } } \sigma _ { \epsilon } ^ { \gamma } + g _ { \beta ; \gamma \epsilon } ^ { \alpha ^ { \prime } } \sigma _ { \delta } ^ { \gamma } + g _ { \beta ; \gamma } ^ { \alpha ^ { \prime } } \sigma _ { \delta \epsilon } ^ { \gamma } = 0,$$

The coincident limit of this relation is

$$[ g ^ { \alpha ^ { \prime } } _ { \beta ; \gamma \delta } ] + [ g ^ { \alpha ^ { \prime } } _ { \beta ; \delta \gamma } ] = 0,$$

$$2 [ g _ { \beta ; \gamma \delta } ^ { \alpha ^ { \prime } } ] + R \alpha ^ { \prime } _ { \, \beta ^ { \prime } ; \delta ^ { \prime } \gamma ^ { \prime } } = 0,$$

from where it follows that

$$[ g ^ { \alpha } _ { \beta ^ { \prime } ; \gamma \delta } ] = - \frac { 1 } { 2 } R ^ { \alpha ^ { \prime } } _ { \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } }.$$

In analogous way it can be shown that

$$[ g ^ { \alpha } _ { \beta ^ { \prime } ; \gamma ^ { \prime } \delta } ] = - \frac { 1 } { 2 } R ^ { \alpha ^ { \prime } } _ { \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } },$$

$$[ g ^ { \alpha } _ { \beta ^ { \prime } ; \gamma \delta ^ { \prime } } ] = \frac { 1 } { 2 } R ^ { \alpha ^ { \prime } } _ { \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } },$$

and also that

$$g _ { b ^ { \prime } ; c } ^ { a } = \frac { 1 } { 2 } g _ { a ^ { \prime } } ^ { a } g _ { c } ^ { c ^ { \prime } } R _ { \, b ^ { \prime } c ^ { \prime } d ^ { \prime } } ^ { a ^ { \prime } } \sigma ^ { d ^ { \prime } } + O ( \epsilon ^ { 2 } ).$$

Finally, we quote without proof the expansion for a bivector K x,x a ( ′ )

$$K _ { a } ( x, x ^ { \prime } ) = g _ { a } ^ { a ^ { \prime } } \left [ K _ { a ^ { \prime } } - K _ { a ^ { \prime } ; c ^ { \prime } } \sigma ^ { c ^ { \prime } } + \frac { 1 } { 2 } K _ { a ^ { \prime } ; c ^ { \prime } d ^ { \prime } } \sigma ^ { c ^ { \prime } } \sigma ^ { d ^ { \prime } } + O ( \epsilon ^ { 3 } ) \right ].$$

We turn now our attention on the application of these formulas to the calculation of the singular part of the Green function.

$$[ g ^ { \alpha } _ { \beta ^ { \prime } ; \gamma ^ { \prime } \delta ^ { \prime } } ] = \frac { 1 } { 2 } R ^ { \alpha ^ { \prime } } _ { \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } },$$

With the help of (5.137),(5.142) together with (5.116) and (5.130) the limit of coincident (5.136) can be calculated in straightforward manner. The result is

$$B _ { \alpha ^ { \prime } \beta ^ { \prime } } & = [ T _ { \alpha ^ { \prime } \beta } ],$$

$$B _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } & = [ T _ { \alpha ^ { \prime } \beta ; \gamma ^ { \prime } } ] - B _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } },$$

$$B _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } } = [ T _ { \alpha ^ { \prime } \beta ; \gamma ^ { \prime } \delta ^ { \prime } } ] + \frac { 1 } { 2 } B _ { \alpha ^ { \prime } \epsilon ^ { \prime } } R _ { \beta ^ { \prime } ; \gamma ^ { \prime } \delta ^ { \prime } } ^ { \epsilon ^ { \prime } } - B _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } \delta ^ { \prime } } - B _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } ; \delta ^ { \prime } } - B _ { \alpha ^ { \prime } \beta ^ { \prime } \delta ^ { \prime } ; \gamma ^ { \prime } },$$

from where the coefficients B i 1 ..i n and consequently the expansion (5.136) are determined. Finally, in the case that T αβ ( x, x ′ ) is a tensor referred to x one may construct the auxiliary tensor

$$\widetilde { T } _ { \alpha ^ { \prime } \beta ^ { \prime } } ( x, x ^ { \prime } ) = g _ { \beta ^ { \prime } } ^ { \beta } g _ { \alpha ^ { \prime } } ^ { \alpha } T _ { \alpha \beta } ( x, x ^ { \prime } ),$$

and expand it using (5.131)-(5.133) together with (5.114) for ˜ T αβ ( x, x ′ )

$$T _ { \alpha \beta } = g _ { \alpha } ^ { \alpha ^ { \prime } } g _ { \beta } ^ { \beta ^ { \prime } } ( C _ { \alpha ^ { \prime } \beta ^ { \prime } } + C _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma } \sigma ^ { \gamma ^ { \prime } } + C _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } } \sigma ^ { \gamma ^ { \prime } } \sigma ^ { \delta ^ { \prime } } ) + O ( \epsilon ^ { 3 } ),$$

with

$$C _ { \alpha ^ { \prime } \beta ^ { \prime } } & = [ T _ { \alpha ^ { \prime } \beta } ],$$

$$C _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } } = [ T _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } } ] - C _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } },$$

$$C _ { \alpha ^ { \prime } \beta ^ { \prime } \gamma ^ { \prime } \delta ^ { \prime } } = [ T _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } \delta ^ { \prime } } ] + \frac { 1 } { 2 } C _ { \alpha ^ { \prime } \ell ^ { \prime } } R _ { \beta ^ { \prime } ; \gamma ^ { \prime } \delta ^ { \prime } } ^ { \ell ^ { \prime } } + \frac { 1 } { 2 } C _ { \ell ^ { \prime } \beta ^ { \prime } } R _ { \alpha ^ { \prime } ; \gamma ^ { \prime } \delta ^ { \prime } } ^ { \ell ^ { \prime } } - C _ { \alpha ^ { \prime } \beta ^ { \prime } ; \gamma ^ { \prime } \delta ^ { \prime } } - C _ { \alpha ^ { \prime } \beta ^ { \prime } \delta ^ { \prime } ; \gamma ^ { \prime } }.$$

these formulas determine the expansion to order two, and can be continued to arbitrary order.

The formulas described above are valid for an arbitrary bitensor T ab ( x, x ′ ). To give some concrete example consider for instance T ab = σ ab . The coincident point limits of this bitensor can be calculated directly from (5.116)-(5.130). By use of this and the recurrence (5.131)-(5.133) it follows that

$$\sigma _ { a b } = g _ { a } ^ { \, a ^ { \prime } } g _ { b } ^ { \, b ^ { \prime } } \left [ g _ { a ^ { \prime } b ^ { \prime } } - \frac { 1 } { 3 } R _ { a ^ { \prime } c ^ { \prime } b ^ { \prime } d ^ { \prime } } \sigma ^ { c ^ { \prime } } \sigma ^ { d ^ { \prime } } + \frac { 1 } { 4 } R _ { a ^ { \prime } c ^ { \prime } b ^ { \prime } d ^ { \prime } ; c ^ { \prime } } \sigma ^ { ; c ^ { \prime } } \sigma ^ { d ^ { \prime } } \sigma ^ { e ^ { \prime } } + O ( \epsilon ^ { 4 } ) \right ],$$

$$\sigma _ { a ^ { \prime } b } = - g _ { b } ^ { b ^ { \prime } } \left [ g _ { a ^ { \prime } b ^ { \prime } } + \frac { 1 } { 6 } R _ { a ^ { \prime } c ^ { \prime } b ^ { \prime } d ^ { \prime } } \sigma ^ { c ^ { \prime } } \sigma ^ { d ^ { \prime } } + O ( \epsilon ^ { 3 } ) \right ]$$

$$\sigma _ { a ^ { \prime } b ^ { \prime } } = g _ { a ^ { \prime } b ^ { \prime } } - \frac { 1 } { 3 } R _ { a ^ { \prime } c ^ { \prime } b ^ { \prime } d ^ { \prime } } \sigma ^ { c ^ { \prime } } \sigma ^ { d ^ { \prime } } + O ( \epsilon ^ { 3 } ),$$

## 6. Green function for static geometries in three dimensions

The divergences of the electrostatic potential A t at the position of the charge arise due to fact that the source in the Maxwell equations (3.16) have a Dirac delta type of singularity. Something analogous happens when a charge is in front of a perfect conductor, which deforms the field lines and give a net force on the charge. The net force is calculated by subtracting terms in the electrostatic field which are divergent at the charge position. The analogous procedure for curved geometries was started by Haddamard [35]. This technique is the one employed in our calculation of the charge self -force in a BTZ geometry. The Taylor like expansions described in the previous section are specially suited for this purpose.

## 6.1 The Haddamard anzatz

The Maxwell equations (3.16) for an static charge in a static geometry can be written in the following form

$$g _ { s } ^ { i j } \nabla _ { i } \nabla _ { j } A _ { t } - N ^ { i } \partial _ { i } A _ { t } = - 2 \pi g _ { t t } j ^ { t }.$$

Here g ij s ( x ) is the spatial part of the metric and ∇ i the spatial Levi-Civita connection, both evaluated at the observation position x . In addition

$$N _ { a } = \partial _ { a } \log \sqrt { - g _ { t t } }.$$

By expressing the potential as

$$A _ { t } & = - q \sqrt { - g _ { t t } ( x ^ { \prime } ) } G ( x, x ^ { \prime } ),$$

it follows that G x,x ( ′ ) satisfies the equation

$$g _ { s } ^ { i j } \nabla _ { i } \nabla _ { j } G _ { t } - N ^ { i } \partial _ { i } G _ { t } = - 2 \pi \delta ( x, x ^ { \prime } ).$$

Here δ x, x ( ′ ) represents the spatial Dirac delta in curved space. It is characterized by the property that

$$\int f ( x ^ { \prime } ) \delta ( x, x ^ { \prime } ) \sqrt { g _ { s } } d V _ { s } = f ( x ),$$

with f ( x ′ ) an arbitrary test function defined in a neighborhood of x . The Haddamard Green function G x,x ( ′ ) [35] is a solution of (6.158) that has the singularity structure enforced by the Dirac delta source, but does not necessarily respect the boundary conditions of the problem (such as the behavior in the asymptotic region). Our strategy will not consist in calculating the electrostatic potential for the BTZ case by the Haddamard method, since it has been already calculated it in (4.86)-(4.87). But the Haddamard calculation will be useful to identify the singular part of this potential. Once this singular part is removed, the resulting renormalized potential will give directly the net self force on the static charge q .

The space time dimension will be kept generically in the following, and we denote it as d = n +2. If this dimension n is even then the static solution can be expressed as

$$G ( x, x ^ { \prime } ) = \frac { 1 } { n - 1 } \frac { U ( x, x ^ { \prime } ) } { ( 2 \sigma ) ^ { ( n - 1 ) / 2 } },$$

while when n is odd one has [35]

$$G ( x, x ^ { \prime } ) = \frac { 1 } { n - 1 } \frac { U ( x, x ^ { \prime } ) } { ( 2 \sigma ) ^ { ( n - 1 ) / 2 } } + V ( x, x ^ { \prime } ) \log \frac { 2 \sigma } { \lambda } + W ( x, x ^ { \prime } ).$$

Here x ′ is the position of the charge singularity. In both cases one has that U x,x ( ) = 1. The functional form for U , V and W depends on the space time in consideration. For n even one can postulate an expansion of the form

$$U ( x, x ^ { \prime } ) = \sum _ { p = 0 } ^ { \infty } U _ { p } ( x, x ^ { \prime } ) ( 2 \sigma ) ^ { p },$$

which, when inserted into (6.158) gives the following recurrence formula

$$( 2 \sigma ^ { a } \nabla _ { a } - N ^ { a } \sigma _ { a } + \nabla ^ { 2 } \sigma + 2 p - n - 1 ) U _ { p } = - \frac { 2 p - n + 1 } { ( n - 1 ) ^ { 2 } } ( \nabla ^ { 2 } - N ^ { a } \nabla _ { a } ) U _ { p - 1 }, \quad \quad ( 6. 1 6 1 )$$

For n odd instead it is postulated that

$$U ( x, x ^ { \prime } ) = \sum _ { p = 0 } ^ { \frac { 1 } { 2 } ( n - 3 ) } U _ { p } ( x, x ^ { \prime } ) ( 2 \sigma ) ^ { p }$$

$$V ( x, x ^ { \prime } ) = \sum _ { p = 0 } ^ { \infty } V _ { p } ( x, x ^ { \prime } ) ( 2 \sigma ) ^ { p }, \quad W ( x, x ^ { \prime } ) = \sum _ { p = 0 } ^ { \infty } W _ { p } ( x, x ^ { \prime } ) ( 2 \sigma ) ^ { p },$$

and sustitution into (6.160) and (6.158) gives instead the following recurrence

$$( 2 \sigma ^ { a } \nabla _ { a } - N ^ { a } \sigma _ { a } + \nabla ^ { 2 } \sigma - 2 ) V _ { 0 } = - \frac { 1 } { 2 ( n - 1 ) } ( \nabla ^ { 2 } - N ^ { a } \nabla _ { a } ) U _ { \frac { 1 } { 2 } ( n - 3 ) },$$

$$( 2 \sigma ^ { a } \nabla _ { a } - N ^ { a } \sigma _ { a } + \nabla ^ { 2 } \sigma + 2 p - 2 ) V _ { p } = - \frac { 1 } { 2 p } ( \nabla ^ { 2 } - N ^ { a } \nabla _ { a } ) V _ { p - 1 },$$

together with

$$( 2 \sigma ^ { a } \nabla _ { a } - N ^ { a } \sigma _ { a } + \nabla ^ { 2 } \sigma + 2 p - 2 ) W _ { p } = - \frac { 1 } { n } ( 2 \sigma ^ { a } \nabla _ { a } - N ^ { a } \sigma _ { a } + \nabla ^ { 2 } \sigma + 4 p - 2 ) V _ { p }$$

$$- \frac { 1 } { 2 p } ( \nabla ^ { 2 } - N ^ { a } \nabla _ { a } ) W _ { p - 1 }.$$

$$( 2 \sigma ^ { a } \nabla _ { a } - N ^ { a } \sigma _ { a } + \nabla ^ { 2 } \sigma + 2 p - 2 ) W _ { p } = - \frac { 1 } { p } ( 2 \sigma ^ { a } \nabla _ { a } \\ - \, \frac { 1 } { 2 p } ( \nabla ^ { 2 } - N ^ { a } \nabla _ { a } ) W _ { p - }$$

These equations should be supplemented with (6.161), which also applies to the odd case. The freedom in choosing the parameter λ or W x,x 0 ( ′ ) corresponds to the gauge transformations in the model. In four dimensions, these type of recurrences where already considered in [34].

## 6.2 Singular terms of the Green function in 2+1 dimensions

The recurrence described in the previous subsection has been analyzed in several situations, for instance in the context of black holes in five dimensions [36]. Nevertheless, to the best of our knowledge, it was not applied to three dimensional cases. For this reason we made an independent analysis by use of the Synge calculus and the Hadamard anzatz described in previous section. Our analysis goes as follows. From (6.160) it is inferred that for 2 + 1 dimensions, which corresponds to n = 1, the biscalar

U x,x ( ′ ) can be set to zero redefining W x,x ( ′ ). Thus the recurrence (6.166) does not play any role in this case. The equation (6.164) for V 0 and (6.165) for V 1 are reduced in this case to

$$( 2 \sigma ^ { a } \nabla _ { a } - N ^ { a } \sigma _ { a } + \nabla ^ { 2 } \sigma - 2 ) V _ { 0 } = 0.$$

$$( 2 \sigma ^ { a } \nabla _ { a } - N ^ { a } \sigma _ { a } + \nabla ^ { 2 } \sigma ) V _ { 1 } = - \frac { 1 } { 2 } ( \nabla ^ { 2 } - N ^ { a } \nabla _ { a } ) V _ { 0 }.$$

The singular part of the Green function we are interested in is not just the collection of terms which are divergent when x → x ′ , but also those whose first derivatives are divergent in that limit, since they give an infinite force. To solve (6.167) one may postulate

$$V _ { 0 } = 1 + a _ { \alpha ^ { \prime } } \sigma ^ { a ^ { \prime } } + \frac { 1 } { 2 } a _ { a ^ { \prime } b ^ { \prime } } \sigma ^ { a ^ { \prime } } \sigma ^ { b ^ { \prime } } + \frac { 1 } { 2 } a _ { a ^ { \prime } b ^ { \prime } c ^ { \prime } } \sigma ^ { a ^ { \prime } } \sigma ^ { b ^ { \prime } } \sigma ^ { c ^ { \prime } } + O ( \epsilon ^ { 4 } ),$$

with a i 1 ..i n coefficients to be determined. In the following the attention will be restricted to terms of order three, since they will contain all the singular pieces. By introducing (6.169) into (6.167) and taking into account the identities deduced in (5.150)-(5.153) it follows that

$$a _ { a ^ { \prime } } & = - \frac { 1 } { 2 } N _ { a ^ { \prime } }, & ( 6. 1 7 0 )$$

$$a _ { a ^ { \prime } b ^ { \prime } } & = \frac { 1 } { 2 } N _ { a ^ { \prime } ; b ^ { \prime } } + \frac { 1 } { 4 } N _ { a ^ { \prime } } N _ { b ^ { \prime } } + \frac { 1 } { 6 } R _ { a ^ { \prime } b ^ { \prime } }, \\ _ { \cdot } & \quad.$$

$$a _ { a ^ { \prime } b ^ { \prime } e ^ { \prime } } = - \frac { 1 } { 2 } N _ { ( a ^ { \prime } ; b ^ { \prime } e ^ { \prime } ) } - \frac { 3 } { 4 } N _ { ( a ^ { \prime } } N _ { b ^ { \prime } ; e ^ { \prime } ) } - \frac { 1 } { 8 } N _ { a ^ { \prime } } N _ { b ^ { \prime } } N _ { e ^ { \prime } } - \frac { 1 } { 4 } N _ { ( a ^ { \prime } } R _ { b ^ { \prime } e ^ { \prime } ) } - \frac { 1 } { 4 } R _ { ( a ^ { \prime } b ^ { \prime } ; e ^ { \prime } ) }.$$

In deriving this result one has to take into account the following result

$$N _ { a } = g _ { a } ^ { a ^ { \prime } } \left [ N _ { a ^ { \prime } } - N _ { a ^ { \prime } ; c ^ { \prime } } \sigma ^ { c ^ { \prime } } + \frac { 1 } { 2 } N _ { a ^ { \prime } ; c ^ { \prime } d ^ { \prime } } \sigma ^ { c ^ { \prime } } \sigma ^ { d ^ { \prime } } + O ( \epsilon ^ { 3 } ) \right ],$$

which follows from (5.154) by identifying K a with N a . Consider now (6.168). To solve it, it is enough to consider the following terms

$$V _ { 1 } & = b + b _ { a ^ { \prime } } \sigma ^ { a ^ { \prime } } + O ( \epsilon ^ { 2 } ).$$

By introducing (6.174) into (6.168), and taking into account (5.150)-(5.153), (6.169) and (6.170)(6.172) the following result is obtained

$$b = \frac { 1 } { 4 } a _ { a ^ { \prime } } ^ { a ^ { \prime } } - \frac { 1 } { 4 } N ^ { a ^ { \prime } } a _ { a ^ { \prime } },$$

$$b _ { a ^ { \prime } } = - \frac { 1 } { 2 } b N _ { a ^ { \prime } } - \frac { 1 } { 8 } a _ { b ^ { \prime } } ^ { b ^ { \prime } } + \frac { 1 } { 1 2 } a _ { b ^ { \prime } } R _ { a ^ { \prime } } ^ { b ^ { \prime } } - \frac { 1 } { 8 } N ^ { b ^ { \prime } } a _ { b ^ { \prime } a ^ { \prime } } + \frac { 1 } { 8 } a _ { b ^ { \prime } } N _ { ; a ^ { \prime } } ^ { b ^ { \prime } },$$

which, by (6.170)-(6.172) can be finally expressed as

$$b = - \frac { 1 } { 8 } ( N _ { ; a ^ { \prime } } ^ { a ^ { \prime } } - \frac { 1 } { 2 } N ^ { a ^ { \prime } } N _ { a ^ { \prime } } + \frac { 1 } { 3 } R ^ { \prime } )$$

$$b _ { a ^ { \prime } } = \frac { 1 } { 1 6 } ( \nabla ^ { ^ { \prime } 2 } N _ { a ^ { \prime } } + N _ { ; b ^ { \prime } } ^ { b ^ { \prime } } N _ { a ^ { \prime } } - N ^ { b ^ { \prime } } N _ { a ^ { \prime } ; b ^ { \prime } } + \frac { 1 } { 2 } N ^ { b ^ { \prime } } N _ { b ^ { \prime } } N _ { a ^ { \prime } } + \frac { 1 } { 3 } R ^ { \prime } N _ { a ^ { \prime } } + \frac { 1 } { 3 } R _ { ; a ^ { \prime } } ^ { \prime } ),$$

The equations (6.166) for W i are not relevant for us, since they do not contain any singularity. Therefore from (6.167)-(6.168) it follows that the Green function is of the form

$$G ( x, x ^ { \prime } ) = \left [ 1 + 2 b \sigma + a _ { a ^ { \prime } } \sigma ^ { a ^ { \prime } } + 2 b _ { a ^ { \prime } } \sigma ^ { a ^ { \prime } } \sigma + \frac { 1 } { 2 } a _ { a ^ { \prime } B ^ { \prime } } \sigma ^ { a ^ { \prime } } \sigma ^ { b ^ { \prime } } + \frac { 1 } { 2 } a _ { a ^ { \prime } B ^ { \prime } c ^ { \prime } } \sigma ^ { a ^ { \prime } } \sigma ^ { b ^ { \prime } } \sigma ^ { c ^ { \prime } } \right ] \log \frac { \sigma ( x, x ^ { \prime } ) } { \lambda } + \dots \ \ ( 6. 1 7 7 )$$

with the coefficients given by (6.170)-(6.172) y (6.175)- (6.176). This expression is valid for any 2+1 dimensional static space time.

The terms in this expansion whose derivatives are divergent when x → x ′ are

$$G _ { s i n g } ( x, x ^ { \prime } ) = \left ( 1 - \frac { g _ { t t ; a ^ { \prime } } } { 4 g _ { t t } } \sigma ^ { a ^ { \prime } } \right ) \log \frac { \sigma ( x, x ^ { \prime } ) } { \lambda } + \dots$$

the other terms give no singularities. This formula combined with (6.157) gives the singular part of the electrostatic potential. This is the expression we were looking for. We turn our attention in the application of this formula for the electrostatic problem in the BTZ geometry.

## 7. Electrostatic self force in BTZ geometries

The singular terms of the Green function calculated in (7.183) are generic for any static three dimensional geometry. In the present section the Green function formalism will be specialized to the non rotating BTZ geometry, and applies to the calculation of the electrostatic self force of a static charge q in the outer region. The differences between the black hole and wormhole case are a theoretical experiment for distinguishing both cases without reaching the throat or the horizon.

In order to solve this task it is necessary to find singular part of the Green function (7.183) for the non rotating BTZ geometry. At first sight, the formula requires the calculation of σ x, x ( ′ ) corresponding to the BTZ metric. However symmetry arguments show that charge self force is radial in the BTZ geometry. Thus, it seems reasonable to limit oneself to the case in which x and x ′ lie on the same radial line, which can be chosen as θ = 0 without loosing generality. The distance corresponding to this situation was calculated already in (2.13) for the case M = 1 and it follows that

$$r = l \cosh d ( r, r _ { h } ),$$

with d r, r ( h ) the distance between the point r and the point of the horizon r h located at the same radial line θ = 0. In fact, by introducing the coordinate s defined by

$$r = l \cosh s, \quad s = \log \left ( \frac { r } { l } + \sqrt { \frac { r ^ { 2 } } { l ^ { 2 } } - 1 } \right )$$

the spatial BTZ metric with M = 1

$$g _ { 2 } = \frac { d r ^ { 2 } } { \frac { r ^ { 2 } } { l ^ { 2 } } - 1 } + r ^ { 2 } d \theta ^ { 2 },$$

is transformed into

$$g _ { 2 } = l ^ { 2 } ( d s ^ { 2 } + \cosh ^ { 2 } s d \theta ^ { 2 } ),$$

from where it follows that the distance between two points on the same radial line is simply

$$d ( r _ { 2 }, r _ { 1 } ) = l ( s _ { 2 } - s _ { 1 } ),$$

or in radial coordinates

$$d ( r, r ^ { \prime } ) = l \log \left ( \frac { r } { l } + \sqrt { \frac { r ^ { 2 } } { l ^ { 2 } } - 1 } \right ) - l \log \left ( \frac { r ^ { \prime } } { l } + \sqrt { \frac { r ^ { \prime 2 } } { l ^ { 2 } } - 1 } \right ),$$

In these terms the singular part of the Green function for two points on the same radial line θ = 0 is given by

$$\Lambda _ { \text{sink} } ^ { \quad \cdot } G _ { \text{sink} } ( r, r ^ { \prime } ) = & \left [ 1 - \frac { r } { 2 l ( \frac { r ^ { 2 } } { l ^ { 2 } } - 1 ) ^ { 3 / 2 } } \left ( \log \left ( \frac { r } { l } + \sqrt { \frac { r ^ { 2 } } { l ^ { 2 } } - 1 } \right ) - \log \left ( \frac { r ^ { \prime } } { l } + \sqrt { \frac { r ^ { \prime 2 } } { l ^ { 2 } } - 1 } \right ) \right ) \right ] \\ & \log \left ( \frac { 1 } { \lambda } \log \left ( \frac { r } { l } + \sqrt { \frac { r ^ { 2 } } { l ^ { 2 } } - 1 } \right ) - \frac { 1 } { \lambda } \log \left ( \frac { r ^ { \prime } } { l } + \sqrt { \frac { r ^ { \prime 2 } } { l ^ { 2 } } - 1 } \right ) \right ) + \dots \\ \text{The regular or normalized notential for the RTZ black hole is}$$

The regular or renormalized potential for the BTZ black hole is

$$A _ { r e n } ^ { a } ( r, r ^ { \prime } ) = A _ { t } ( r, r ^ { \prime } ) - A _ { s i n g } ^ { a } ( r, r ^ { \prime } ),$$

where the singular part

$$A _ { s i n g } ^ { a } ( r, r ^ { \prime } ) = - q \sqrt { \frac { r ^ { \prime 2 } } { l ^ { 2 } } - 1 } \, G _ { s i n g } ( r, r ^ { \prime } ),$$

follows from (6.157) and (7.183) and A r,r t ( ′ ) is the black hole solution (4.87) evaluated at θ = 0, namely

$$A _ { t } ^ { I } ( r, \theta ) = - q \log ( r _ { 0 } ) + \sum _ { n = 1 } ^ { \infty } \frac { q n ^ { 2 } k _ { n } ( r _ { 0 } ) h _ { n } ( r ) } { 8 M },$$

for r &lt; r 0 and

$$A _ { t } ^ { I I } ( r, \theta ) = - q \log ( r ) + \sum _ { n = 1 } ^ { \infty } \frac { q n ^ { 2 } h _ { n } ( r _ { 0 } ) k _ { n } ( r ) } { 8 M },$$

for r &gt; r 0 . The self force for the charge q is calculated by

$$F ^ { r ^ { \prime } } = - q F ^ { ( r ^ { \prime } ) } _ { ( t ) } u ^ { ( t ) } = - q ( e ^ { ( r ^ { \prime } ) } ) _ { r ^ { \prime } } F ^ { r ^ { \prime } } _ { t } u ^ { t } = - q \sqrt { g _ { r ^ { \prime } r ^ { \prime } } } g ^ { r ^ { \prime } r ^ { \prime } } F _ { r ^ { \prime } t } \frac { 1 } { \sqrt { - g _ { t t } } } = - q \partial _ { r ^ { \prime } } A _ { r e n } ( r ^ { \prime } ),$$

with

$$A _ { r e n } ( r ^ { \prime } ) = \lim _ { r \rightarrow r ^ { \prime } } A _ { r e n } ( r, r ^ { \prime } ).$$

The self force described above work fine in even space time dimensions, and was applied by one of the authors [10] to reproduce correctly the Will-Simth self force in a Schwarzschild black hole [ ? ]. However, in odd dimensions there exist an axiomatic approach to the self force which suggest that such force may pick additional terms [38]. This procedure is equivalent to the standard one in four dimensions, but not necessarily in even ones. We will denote the self-force (7.187) as the 'naive' self force and in the next sections the effect of the additional terms that arise by the methods proposed in [38] will be considered separately.

Although the expressions found above seem to be easy to deal with, there is a problem when limiting oneself to radial lines. The calculation for the self-force involves taking the coincident limit r → r ′ . The complication is that the electrostatic potential A a sing ( r, r ′ ) has a singularity in this limit, which

is explicitly seen by looking (7.183). The full potential (7.185)-(7.186) has the same singularity with opposite sign and the combination (7.184) is free of singularities by construction. But one expression is given as an infinite series and the other is a single term expression. And this is not suitable for solving the ∞-∞ indetermination. One approach may be to sum the series. This can be done in four dimensions, and the resulting potential is a simple algebraic expression [3]. Unfortunately, we ignore if there exist such summation formulas in three dimensions. Thus, a different approach should be employed. The one to be used below is based on the following observation: if one is able to expand the Green function as a Fourier series of the form

$$G ^ { s } ( x, x ^ { \prime } ) = \frac { G _ { 0 } ( r, r ^ { \prime } ) } { 2 } + \sum _ { n = 1 } ^ { \infty } G _ { n c } ^ { s } ( r, r ^ { \prime } ) \cos ( n \theta ) + \sum _ { n = 0 } ^ { \infty } G _ { n s } ^ { s } ( r, r ^ { \prime } ) \sin ( n \theta ),$$

then both (7.188) and (4.87) can be combined into a single non divergent series, which can be approximated to an arbitrary order to obtain the approximated self force. The symmetry argument requires to take the limit θ → 0 and r → r ′ in order to calculate it. However this Fourier expansion can be performed only if we know σ x, x ( ′ ) for points on the space time with arbitrary θ values. This function can be calculated explicitly for the BTZ geometry, and the calculation is performed in the appendix. It is given by

$$\sigma ( r, r ^ { \prime }, \theta, \theta ^ { \prime } ) = \frac { 1 } { 2 } d ^ { 2 } ( r, r ^ { \prime }, \theta, \theta ^ { \prime } )$$

with

$$d ( \mathbf x _ { 1 }, \mathbf x _ { 2 } ) = l \cosh ^ { - 1 } \left [ \frac { r _ { 1 } r _ { 2 } } { r _ { h } ^ { 2 } } \cosh ( \theta _ { 2 } - \theta _ { 1 } ) - \sqrt { \left ( \frac { r _ { 1 } } { r _ { h } } \right ) ^ { 2 } - 1 } \sqrt { \left ( \frac { r _ { 2 } } { r _ { h } } \right ) ^ { 2 } - 1 } \right ],$$

An important consistency test of (7.189) is to recover the radial distance (7.181) when θ 2 → θ 1 , that is, when both points lie on the same radial line. For this, it is convenient to express (7.189) in terms of the coordinates ( s, θ ) defined in (7.179). The result is

$$d ( x _ { 1 }, x _ { 2 } ) = l \cosh ^ { - 1 } \left [ \cosh s _ { 1 } \cosh s _ { 2 } \cosh ( \theta _ { 2 } - \theta _ { 1 } ) - \sinh s _ { 1 } \sinh s _ { 2 } \right ].$$

From the last expression it follows that when the points are on the same radial line then α = 0 and the distance becomes

$$d ( x _ { 1 }, x _ { 2 } ) = l \cosh ^ { - 1 } \left [ \cosh s _ { 1 } \cosh s _ { 2 } - \sinh s _ { 1 } \sinh s _ { 2 } \right ] = l \cosh ^ { - 1 } \left [ \cosh ( s _ { 2 } - s _ { 1 } ) \right ] = l ( s _ { 2 } - s _ { 1 } ).$$

Thus, the radial distance (7.181) has been recovered in this limit, as expected. Since the expression (7.189) for the geodesic distance between two arbitrary points in the BTZ geometry is explicit, the calculation of the singular Green function (7.183) for the electrostatic problem is immediate. Unfortunately, the expression that is obtained is very complicated and we were not able to find a closed expression for the Fourier coefficients of the expansion. But this is not to be discouraged, since there is a numerical trick that can be employed, which will allow us to approximate the real value of the self force to an arbitrary order. We turn our attention to this trick in the next section.

## 7.1 Fourier expansion of the singular Green function

In the previous section the singular part of the Green function has been found explicitly, but not its Fourier expansion. The analytic form of this expansion turns out to be out of our computational methods. In these situations, the following observation may be useful. In some mathematical applications it may be of interest to study a given function f ( x ) not in the full domain, but in an small interval ( x 0 -glyph[epsilon1], x 0 + ) around a point glyph[epsilon1] x 0 . For example, this may be because the points around x 0 are the ones that considerably influence the value of some integral. In these situations it may be valid to find a trial function g x ( ) such that g x ( 0 ) = f ( x 0 ) and such that they values are very close in the mentioned interval. If this new function g x ( ) has the additional property that its Fourier expansion is simpler, then it may be advantageous to work with it by considering a large number of terms of the Fourier expansion. For example, one can approximate the function f ( x ) around x 0 with a rectangle function with small width and with f ( x 0 ) as its height. Note that if all terms are considered, both the exact and the approximated Fourier expansions have the same value at x 0 , under suitable conditions for f ( x ) and g x ( ).

Based on this idea, our goal is to construct an approximated Green G a ( r, r , θ, θ ′ ′ ) function which coincides with the exact one when the two points x and x ′ lay on the same radial line θ = θ ′ . Since the symmetry of the problem implies that the force is obtained taking the limit r → r ′ for points lying on the same radial line, then the fact that this approximated Green function does not coincide with the exact one outside this line is not relevant for the calculation. This approximated Green function is found in terms of an approximated distance d a ( x, x ′ ) which coincide with (7.181) for two points located on the same radial line.

At first sight there are a large number of possible choices for the approximated Green function, and the idea is choose one which is simple to work with. A physically motivated one is given in the appendix for mass value M = 1. The approximated geodesic distance that is obtained is

$$d ^ { a } ( x, x _ { 0 } ) = l \sqrt { 1 + \frac { ( \theta - \theta _ { 0 } ) ^ { 2 } } { \left ( \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s _ { 0 } } { 2 } \right ) ^ { 2 } } } ( s - s _ { 0 } ), \quad 0 < | \theta - \theta _ { 0 } | < \pi,$$

$$d ^ { a } ( x, x _ { 0 } ) = l \sqrt { 1 + \frac { ( 2 \pi - \theta + \theta _ { 0 } ) ^ { 2 } } { \left ( \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s _ { 0 } } { 2 } \right ) ^ { 2 } } } ( s - s _ { 0 } ), \quad \pi < | \theta - \theta _ { 0 } | < 2 \pi. \quad ( 7. 1 9 2 )$$

$$d ^ { a } ( x, x _ { 0 } ) & = l \sqrt { 1 + \frac { ( \theta - \theta _ { 0 } ) ^ { 2 } } { \left ( \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s _ { 0 } } { 2 } \right ) } } ( s - s _ { 0 } ), \quad \ 0 < | \theta - \theta _ { 0 } | < \pi, \quad \ \ ( 7. 1 9 1 ) \\ d ^ { a } ( x, x _ { 0 } ) & = l \sqrt { 1 + \frac { ( 2 \pi - \theta + \theta _ { 0 } ) ^ { 2 } } { \left ( \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s _ { 0 } } { 2 } \right ) } } ( s - s _ { 0 } ), \quad \ \ \pi < | \theta - \theta _ { 0 } | < 2 \pi. \quad \ \ ( 7. 1 9 2 )$$

We will employ this expression in the following by choosing the θ values in the range -π &lt; θ &lt; π . By construction (7.192) is not strictly true, but tends to (7.181) when θ 1 = θ 2 , and it is even under the interchange θ ↔-θ , which is a property of the true distance. Note in addition, that the limit

$$\lim _ { s \to s _ { 0 } } \frac { ( s - s _ { 0 } ) ^ { 2 } } { \left ( \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s _ { 0 } } { 2 } \right ) ^ { 2 } } = \cosh ^ { 2 } ( s _ { 0 } ),$$

glyph[negationslash]

implies that the distance (7.191) is perfectly regular when θ = θ 0 . The same argument holds for (7.192).

The approximated world function σ a ( x, x ′ ) can be found directly found by taking the square of the distance element (7.191). It is convenient to locate the charge at θ 0 = 0 by simplicity. Then the singular Green function (7.183) takes the following form

$$G _ { s i n g } ^ { a } ( x, x ^ { \prime } ) = \left [ f ( r, r ^ { \prime } ) \theta ^ { 2 } + g ( r, r ^ { \prime } ) \right ] \log [ a ( r, r ^ { \prime } ) \theta ^ { 2 } + b ( r, r ^ { \prime } ) ],$$

with

$$a ( r, r ^ { \prime } ) = \frac { ( s - s ^ { \prime } ) ^ { 2 } } { 2 \left ( \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s ^ { \prime } } { 2 } \right ) ^ { 2 } },$$

$$& \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \colon \\ & b ( r, r ^ { \prime } ) = \frac { ( s - s ^ { \prime } ) ^ { 2 } } { 2 }, \\ f ( r, r ^ { \prime } ) = - \frac { r } { 2 } \frac { d } { d r ^ { \prime } } \left [ \frac { ( s - s ^ { \prime } ) ^ { 2 } } { ( \tan ^ { - 1 } \tanh \frac { s ^ { \prime } } { 2 } - \tan ^ { - 1 } \tanh \frac { s ^ { \prime } } { 2 } ) ^ { 2 } \right ]. \\ g ( r, r ^ { \prime } ) = 2 - \frac { r } { 2 } \frac { d } { d r ^ { \prime } } ( s - s ^ { \prime } ) ^ { 2 }, \\ \text{is given by } ( 7. 1 7 9 ). \text{ It is convenient to expand this function in the basis exp} ( i n \theta ) \text{ as}$$

Here s ′ ( r ′ ) is given by (7.179). It is convenient to expand this function in the basis exp( inθ ) as

$$G ^ { a } _ { s i n g } ( x, x ^ { \prime } ) = \frac { G _ { 0 } ( r, r ^ { \prime } ) } { 2 } + \sum _ { n = 1 } ^ { \infty } G ^ { a } _ { n c } ( r, r ^ { \prime } ) \cos ( n \theta ) + \sum _ { n = 0 } ^ { \infty } G ^ { a } _ { n s } ( r, r ^ { \prime } ) \sin ( n \theta ),$$

with the radial coefficients given by

$$G _ { n c } ^ { a } ( r, r ^ { \prime } ) = 2 \int _ { 0 } ^ { \pi } G ( x, x ^ { \prime } ) \cos ( n \theta ) d \theta,$$

and the analogous definition for G a ms ( x, x ′ ). Their explicit value of these integrals is

$$e \, \text{analogous definition for } G ^ { a } _ { m s } ( x, x ^ { \prime } ). \, \text{Their explicit value of these integrals is} \\ G ^ { a } _ { n c } ( r, r ^ { \prime } ) = \frac { 2 i } { a n ^ { 3 } } \Psi ^ { c } _ { n } \left [ C i \left ( n \pi - i n \sqrt { \frac { b } { a } } \right ) - C i \left ( n \pi + i n \sqrt { \frac { b } { a } } \right ) \\ - C i \left ( - i n \sqrt { \frac { b } { a } } \right ) + C i \left ( i n \sqrt { \frac { b } { a } } \right ) \right ] \\ + \, \frac { 4 } { a n ^ { 3 } } \Psi ^ { s } _ { n } \left [ S i \left ( n \pi + i n \sqrt { \frac { b } { a } } \right ) - S i \left ( - n \pi + i n \sqrt { \frac { b } { a } } \right ) \right ] \\ + \, \frac { 2 } { a n ^ { 3 } } \left [ 4 \, a \, f \, n \, \pi \, \cos ( n \pi ) \left ( \log ( a \pi ^ { 2 } + b ) + 1 \right ) \right ], \\ \, \text{we have introduced} \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,
 \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \  \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \cdot \, \, \,$$

where we have introduced

$$\Psi _ { n } ^ { c } = ( 2 a f - a g n ^ { 2 } + b f n ^ { 2 } ) \sinh \left ( n \sqrt { \frac { b } { a } } \right ) - 2 \sqrt { a b } f n \cosh \left ( n \sqrt { \frac { b } { a } } \right ),$$

$$\Psi _ { n } ^ { s } = ( 2 a f - a g n ^ { 2 } + b f n ^ { 2 } ) \cosh \left ( n \sqrt { \frac { b } { a } } \right ) - 2 \sqrt { a b } f n \sinh \left ( n \sqrt { \frac { b } { a } } \right ),$$

by simplicity. For n = 0 one has

$$\hat { G } _ { 0 } ^ { a } ( r, r ^ { \prime } ) = 2 \log ( a \pi ^ { 2 } + b ) \left ( g \pi + \frac { f \pi ^ { 3 } } { 3 } \right ) + 4 \left ( g - \frac { f b } { 3 a } \right ) \left [ \sqrt { \frac { b } { a } } \tan ^ { - 1 } \left ( \pi \sqrt { \frac { a } { b } } \right ) - \pi \right ] - 4 \frac { f \pi ^ { 3 } } { g }.$$

This is the part we will be interested in if the charge is at θ = 0 since sin( nθ ) vanish at this location. The notation Si x ( ) and Ci x ( ) denote the integral sine and cosine, whose definitions are given by

$$S i ( x ) = \int _ { 0 } ^ { x } \frac { \sin t } { t } d t,$$

$$C i ( x ) = - \int _ { x } ^ { \infty } \frac { \cos t } { t } d t.$$

Taking this and (6.157) into account it follows that the singular term of the potential A a sing are given by

$$A _ { s i n g } ^ { a } ( x, x ^ { \prime } ) = \frac { A _ { 0 } ( r, r ^ { \prime } ) } { 2 } + \sum _ { n = 1 } ^ { \infty } A _ { n c } ( r, r ^ { \prime } ) \cos ( n \theta ) + \sum _ { n = 0 } ^ { \infty } A _ { n s } ( r, r ^ { \prime } ) \sin ( n \theta ),$$

with

$$A _ { n s } ( r, r ^ { \prime } ) = - q \sqrt { - g _ { t ^ { \prime } t ^ { \prime } } } G _ { n s } ( r, r ^ { \prime } ) = - q \sqrt { \frac { r ^ { \prime 2 } } { l ^ { 2 } } - 1 } \, G _ { n s } ^ { a } ( r, r ^ { \prime } ),$$

$$A _ { n c } ( r, r ^ { \prime } ) = - q \sqrt { - g _ { t ^ { \prime } t ^ { \prime } } } G _ { n c } ( r, r ^ { \prime } ) = - q \sqrt { \frac { r ^ { \prime 2 } } { l ^ { 2 } } - 1 } \, G _ { n c } ^ { a } ( r, r ^ { \prime } ),$$

and it is completely determined by (7.197).

## 7.2 Analysis of the 'naive' self-force

After calculating the singular terms in the electrostatic potential A t , the renormalized potential A a ren ( r, r ′ ) is then find as

$$A _ { r e n } ^ { a } ( r, r ^ { \prime } ) = A _ { t } ( r, r ^ { \prime } ) - A _ { s i n g } ^ { a } ( r, r ^ { \prime } ),$$

with A r,r t ( ′ ) given by (4.86)-(4.87) for the BTZ black hole and by (4.99) for the BTZ wormhole, and A a sing ( r, r ′ ) is given by (7.198). As shown in (7.187) the self force is simply

$$F ^ { r ^ { \prime } } = - q \partial _ { r ^ { \prime } } A _ { r e n } ( r ^ { \prime } ),$$

with

$$A _ { r e n } ( r ^ { \prime } ) = \lim _ { r \rightarrow r ^ { \prime } } A _ { r e n } ( r, r ^ { \prime } ).$$

For the black hole it follows from (4.87) that

$$A _ { ( b k ) r e n } ^ { a } = - q \log ( r ^ { \prime } ) + q \sqrt { \frac { r ^ { \prime 2 } } { l ^ { 2 } } - 1 } \, G _ { 0 } ( r, r ^ { \prime } ) + q \sum _ { n = 1 } ^ { \infty } \left [ \frac { n ^ { 2 } h _ { n } ( r ^ { \prime } ) k _ { n } ( r ^ { \prime } ) } { 8 } + \sqrt { \frac { r ^ { \prime 2 } } { l ^ { 2 } } - 1 } \, G _ { n c } ( r ^ { \prime } ) \right ],$$

with

$$G _ { n c } ( r ^ { \prime } ) = \lim _ { r \rightarrow r ^ { \prime } } G _ { n c } ^ { a } ( r, r ^ { \prime } ).$$

The calculation of (7.200) requires the knowledge of the following limits

$$\lim _ { r \to r ^ { \prime } } a ( r, r ^ { \prime } ) = 2 \cosh ^ { 2 } s ^ { \prime } = \frac { 2 r ^ { \prime 2 } } { l ^ { 2 } },$$

$$\lim _ { r \to r ^ { \prime } } b ( r, r ^ { \prime } ) = 0,$$

$$\lim _ { r \to r ^ { \prime } } f ( r, r ^ { \prime } ) = - \cosh ^ { 2 } s ^ { \prime } = - \frac { r ^ { \prime 2 } } { l ^ { 2 } }, \\ \lim _ { r \to r ^ { \prime } } g ( r, r ^ { \prime } ) = 2,$$

$$r \rightarrow r ^ { \prime }$$

which follows directly from the definitions (7.195). From this it follows that (7.200) is

$$G _ { n c } ( r ^ { \prime } ) = ( - 1 ) ^ { n } \frac { 4 \pi } { n ^ { 2 } } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \left [ 2 \log \left ( \frac { r } { r _ { h } } \right ) + \log ( 2 \pi ^ { 2 } ) + 1 \right ] + \frac { 2 } { n } \left ( \frac { 2 } { n ^ { 2 } } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } - 1 \right ) S i ( n \pi ),$$

where in this case r h = . For l n = 0 one has

$$G _ { 0 } ( r ^ { \prime } ) = 2 \left [ 2 \log \left ( \frac { r } { r _ { h } } \right ) + \log ( 2 \pi ^ { 2 } ) \right ] \left ( \pi + \frac { \pi ^ { 3 } } { 3 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ) - 4 \pi - \frac { 4 \pi ^ { 3 } } { 9 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } }.$$

It is important to remark that from (7.201)-(7.202) it follows that the Green function at θ = 0 is given by

$$G ( r ) & = 2 \left [ 2 \log \left ( \frac { r } { r _ { h } } \right ) + \log ( 2 \pi ^ { 2 } ) \right ] \left ( \pi + \frac { \pi ^ { 3 } } { 3 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ) - 4 \pi - \frac { 4 \pi ^ { 3 } } { 9 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \\ & \quad + S _ { 1 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \left [ 2 \log \left ( \frac { r } { r _ { h } } \right ) + \log ( 2 \pi ) + 1 \right ] + S _ { 2 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } + S _ { 3 } \\ \text{series $S$} \text{, given } h v$$

with the series S i given by

$$S _ { 1 } = 4 \pi \sum _ { n = 1 } ^ { \infty } ( - 1 ) ^ { n } \frac { 1 } { n ^ { 2 } }, \quad S _ { 2 } = 4 \sum _ { n = 1 } ^ { \infty } \frac { 1 } { n ^ { 3 } }, \quad S _ { 3 } = 2 \sum _ { n = 1 } ^ { \infty } \frac { S i ( n \pi ) } { n }$$

The first two are convergent, the third is divergent. This is as expected, since this is the part of the Green function which is singular at coincident point limits. We turn now the attention for the self-force corresponding to the two boundary conditions discussed in section 4.

## 7.2.1 The self force first type of boundary conditions

The results discussed above together with the potential (4.86)-(4.87) gives the following renormalized potential for the black hole

$$A _ { b h } ^ { a } ( r ) = q \sum _ { n = 1 } ^ { \infty } \left \{ \frac { n ^ { 2 } h _ { n } ( r ) k _ { n } ( r ) } { 8 } + \sqrt { \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } } - 1 \left [ ( - 1 ) ^ { n } \frac { 4 \pi \ r ^ { 2 } } { n ^ { 2 } \ r _ { h } ^ { 2 } } \left ( \log \left ( 2 \pi ^ { 2 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ) + 1 \right ) + \frac { 2 } { n } \left ( \frac { 2 \ r ^ { 2 } } { n ^ { 2 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } } - 1 \right ) S i ( n \pi ) \right ] \right \}$$

$$- \, q \log \left ( \frac { r } { r _ { h } } \right ) + q \sqrt { \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } - 1 } \, \left [ 2 \log \left ( 2 \pi ^ { 2 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ) \left ( \pi + \frac { \pi ^ { 3 } } { 3 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ) - 4 \pi - \frac { 4 \pi ^ { 3 } } { 9 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ]$$

$$A _ { b h } ^ { a } ( r ) = q \sum _ { n = 1 } ^ { \infty } & \left \{ \frac { n ^ { 2 } h _ { n } ( r ) k _ { n } ( r ) } { 8 } + \sqrt { \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } - 1 } \left [ ( - 1 ) ^ { n } \frac { 4 \pi } { n ^ { 2 } } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \left ( \log \left ( 2 \pi ^ { 2 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ) + 1 \right ) + \frac { 2 } { n } \left ( \frac { 2 } { n ^ { 2 } } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } - 1 \right ) S i ( n \pi ) \right ] \right \} \\ & - q \log \left ( \frac { r } { r _ { h } } \right ) + q \sqrt { \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } - 1 } \ \left [ 2 \log \left ( 2 \pi ^ { 2 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ) \left ( \pi + \frac { \pi ^ { 3 } } { 3 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ) - 4 \pi - \frac { 4 \pi ^ { 3 } } { 9 } \frac { r ^ { 2 } } { r _ { h } ^ { 2 } } \right ]$$

The wormhole potential follows from (4.99) and the last expression (7.203), the result is

$$A _ { ( w h ) r e n } ^ { a } ( r ) = A _ { ( b k ) r e n } ^ { a } ( r ) - \frac { q } { 2 } \sum _ { n = 1 } ^ { \infty } \left [ \frac { 1 } { r _ { g } h _ { n } ( r _ { g } ) } + \frac { n ^ { 2 } k _ { n } ^ { \prime } ( r _ { g } ) } { 8 } \right ] \frac { h _ { n } ^ { 2 } ( r ^ { \prime } ) } { h _ { n } ^ { \prime } ( r _ { g } ) }.$$

These expressions is a combination of two series which are divergent at r → r ′ , but whose divergent behavior cancel exactly. Therefore when evaluated at the radial line θ = θ ′ , both (7.203)-(7.204) are finite and has finite derivatives, which give the self force.

The following figures show the approximate behavior for the force F t , for some values of n between n = 15 and n = 45 for the black hole case. This graphs were obtained by use of MATHEMATICA. The graphs show the behavior for large value of r , which corresponds to s ∼ 1. The functional form of the force is always the same in all the orders we considered namely, it starts to grow from the horizon s = 0 till it attach some maximum value, and then decays to zero at s = 1 (the asymptotic region). But as larger orders are considered, the maximum start to grow and to move closer to the asymptotic region s = 1. This suggest us that when all the series is summed up, this maximum becomes infinitely large, and the line s = 1 is an asymptote for the function. In other words, the decaying part is just an artifact of the truncation. If this is correct and is not related to a numerical problem in the program we are using, then it follows that the self force is repulsive and grows indefinitely as the charge goes to the asymptotic region. The repulsive nature is the same as in the Schwarzschild four dimensional counterpart [3]. But is mandatory to interpret the strange behavior at the infinite. The explanation for such non typical behavior may be that the boundary conditions considered are not correct, and include extra charges at the asymptotic region. Thus it is convenient to analyze the 'right' boundary conditions to see if this behavior is improved.

We have also considered the wormhole case. The effect of the throat is to give a contribution to the force of opposite sign for the black hole. This effect is also seen in the Schwarzschild case [10]. In that case, this contribution change the sign of the force in some region of the space time. However, with the numerical precision we are working with, we did not find such change of the behavior in the BTZ geometry, the net effect of the throat is just a change in the numerical value of the force. But this is not conclusive, and perhaps this effect may appear by considering higher orders in the series.

Figure 1: The black hole self force function of s = 1 -r -2 , from s = 0 to s = 0 9 to order . n = 25. The force seems to grow, but the next figure shows that it attach a maximum.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-68477cfa-d0d0-43f6-ab49-a34be7634b30.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-fd71359f-3b3a-46a0-88e2-3a3349a07dae.jpg)

Figure 2: The behavior of the self force near the boundary s = 1. It reaches a maximum and then goes to zero. We argue that this is just an effect of the truncation of the series and that the line s = 1 is in fact an asymptote.

Figure 3: The force for n = 45 in the interval from s = 0 to s = 0 9. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-d6ef0a46-e80e-45d5-9dc6-5d1649bf9dae.jpg)

Figure 4: The behavior near of the horizon s = 1 with n = 45. The qualitative form of the curve is similar to the case n = 25 but the maximum is again more pronounced and closer to s = 1. This suggest that when all the series is summed, the maximum will tend to infinite and the vertical line s = 1 will be the asymptote. See below for an explanation of this seemingly strange result.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-dba94797-2fd2-4ed3-98c7-680eccd7cb59.jpg)

## 7.2.2 The self force corresponding to the right type of boundary conditions

For this type of boundary conditions, the potential follows from (4.91) instead of (4.87). The result is

$$& \dashrightarrow \dashrightarrow r _ { h } \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashrightarrow \dashright \\ & A _ { b h } ^ { a } ( r ) = \frac { q } { 2 } \sum _ { n = 1 } ^ { \infty } \left \{ \Gamma ( 1 + \frac { i n } { 2 } ) \Gamma ( 1 - \frac { i n } { 2 } ) h _ { n } ( r ) f _ { n } ( r ) + \sqrt { \frac { r ^ { 2 } } { r ^ { 2 } _ { h } } - 1 } \left [ ( - 1 ) ^ { n } \frac { 4 \pi } { n ^ { 2 } } \frac { r ^ { 2 } } { r ^ { 2 } _ { h } } \left ( \log \left ( 2 \pi ^ { 2 } \frac { r ^ { 2 } } { r ^ { 2 } _ { h } } \right ) + 1 \right ) \\ & + \frac { 2 } { n } \left ( \frac { 2 } { n ^ { 2 } } \frac { r ^ { 2 } } { r ^ { 2 } _ { h } } - 1 \right ) S i ( n \pi ) \right ] \right \} - q \log \left ( \frac { r } { r _ { h } } \right ) + q \sqrt { \frac { r ^ { 2 } } { r ^ { 2 } _ { h } } - 1 } \, \left [ 2 \log \left ( 2 \pi ^ { 2 } \frac { r ^ { 2 } } { r ^ { 2 } _ { h } } \right ) \left ( \pi + \frac { \pi ^ { 3 } } { 3 } \frac { r ^ { 2 } } { r ^ { 2 } _ { h } } \right ) - 4 \pi - \frac { 4 \pi ^ { 3 } } { r ^ { 2 } _ { h } } \frac { r ^ { 2 } } { r ^ { 2 } _ { h } } \right ] \\ & \quad \end{array} \,.$$

Figure 5: The black hole self force function of s = 1 -r -2 , from s = 0 to s = 0 9 to order . n = 20. The force seems to be positive and divergent near the horizon and at the asymptotic boundary, and there is an intermediate values of s where it is negative.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-8acd5fea-f9d9-4532-ab19-93d487e793f4.jpg)

Figure 6: The behavior of the self force near the boundary s = 1, which shows that it is divergent.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-cbb675ea-9541-4336-a5d9-ddbc2855fdfc.jpg)

The wormhole potential follows from (4.100) and the last expression (7.203), the result is

$$A _ { ( w h ) r e n } ^ { a } ( r ) = A _ { ( b k ) r e n } ^ { a } ( r ) - \frac { q } { 2 } \sum _ { n = 1 } ^ { \infty } \left [ \frac { h _ { n } ( r ) } { r _ { g } h _ { n } ( r _ { g } ) h _ { n } ^ { \prime } ( r _ { g } ) } + \Gamma ( 1 + \frac { i n } { 2 } ) \Gamma ( 1 - \frac { i n } { 2 } ) \frac { h _ { n } ( r ) k _ { n } ^ { \prime } ( r _ { g } ) } { 2 h _ { n } ^ { \prime } ( r _ { g } ) } \right ] h _ { n } ( r ) \ \ ( 7. 2 0 6 )$$

The figures (5)-(8) show the self force behavior for this case. It is surprising to see that the self force is again divergent at the asymptotic boundary and also at the horizon. From a formal point of view,

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-f4eebfeb-f918-40a4-b644-022880f072eb.jpg)

Figure 7: The behavior of the self force near the horizon s = 0, which shows that it is also divergent.

Figure 8: The black hole self force function of s = 1 -r -2 , from s = 0 to s = 0 9 to order . n = 40. The result is qualitative the same as for n = 20.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-8fba6b03-7d6e-4a74-92b6-d0ef516b4c91.jpg)

Figure 9: The wormhole self force when the throat is r g = 3, with s varying from s = 0 to s = 0 9. . The effect of the throat is to increase the slope of the asymptote. It is not clearly seen in this graph if there is a region where the force is negative.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-df9ebd5d-fc9a-4a8f-823e-8f1d144f0ecc.jpg)

this divergent behavior can be seen directly by looking (7.205). The logarithmic terms are multiplied for functions of r and there are terms proportional to r 2 as well in this expression. The derivatives with respect to r of all these terms are not bounded for large r . This may seem the cause for the

Figure 10: The wormhole self force when the throat is r g = 3, with s varying from s = 0 4 to . s = 0 5. . Here it is seen that there is a negative force at some region. The difference with the black hole is that the point where the sign change is closer. So the location of that turning point may be a clue for distinguishing the geometries.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0575v3-figure-78128ea5-65b7-48c4-a2f5-99a3f0c274c3.jpg)

growth of the force at large values. On the other hand, th terms proportional h n ( r f ) n ( r ) go fast to zero since, as discussed in previous sections the derivative of f n ( r ) goes as 1 /r 3 and the one of h n ( r ) goes as 1 /r . Thus these terms does not influence what happens at large r . But they arguably the cause of the divergence near the horizon. As also discussed in previous section the function f n ( r ) possess the fastest decaying conditions. However, the behavior of f n ( r ) at the horizon is highly non regular and, even though h n ( r ) → 0, the full contribution is not bounded. In other words, the fastest decaying conditions also give the worse behavior at the horizon.

The other surprising behavior is that the self force is repulsive near the horizon and far away, but at intermediate radial positions, it is attractive. This is in contrast with the Schwarzschild case in four dimensions, where it is always repulsive [3]. However, without a proper understanding to the origin of the divergences, this analysis does not make any sense. We turn to this point below.

## 8. Is the force calculated the real one?

The procedure of previous section computes a regular Green function. However, there may exist other regular Green functions as well, since the addition to a solution of the homogeneous equation is still a full solution of the problem. Therefore the statement that the calculated self-force is the real one must be taken with care. To give an example, the half-advanced-half-retarded Green function, and the regularized ('R') force of the Detweiler-Whiting decomposition [29] are regular at the location of the particle. Nevertheless only the 'R' force corresponds to the correct expression for the self-force.

The method presented in the previous section works fine in even space time dimensions. In fact, this method was applied by one of the authors in [10] to reproduce correctly the Will-Smith self force for a Schwarzschild black hole in four dimensions [3]. However, the BTZ geometry is odd dimensional and this case is more tricky, as suggested by the axiomatic approach presented in [38] and applied recently for five dimensional space times in [36]. This approach introduce extra terms in our calculation. These

terms vanish in four dimensions and do not spoil the results obtained in [10]. But in odd dimensions, this may not be the case. A further motivation for studying these terms is that they may cure the divergent behavior at the asymptotic boundary found in the previous section. However, it will be shown below that this is not the case.

The self force calculated in the previous sections was given by

$$F _ { a } = \partial _ { a } A - \partial _ { a } A _ { s i n g },$$

with A the full potential satisfying the specific boundary conditions of the problem and A sing the part of the potential whose values and the values of its derivatives ∂ A a sing are divergent when the position x of the observation point tends to the charge position x ′ . The approach of [38] states instead that the self force should be calculated as

$$F _ { a } = \partial _ { a } A - \partial _ { a } A _ { s i n g } + < \partial _ { a } A _ { s i n g } >,$$

where &lt; ∂ A a sing &gt; is the average of the quantity inside over an small surface witth s = cte in the limit s → 0, from which all the contributions that are divergent in the limit s → 0 are removed. As it will be shown below, this average will pick some additional quantities in the calculation. The derivative is ∂ a is taken on a point of that surface. Now, since the singular part is given by

$$A _ { s i n g } = - q \sqrt { \frac { r ^ { \prime 2 } } { l ^ { 2 } } } - 1 G _ { s i n g } ( x, x ^ { \prime } ),$$

and the derivatives are taken at x , it follows that the relevant derivative to be taken is

$$\partial _ { a } G _ { s i n g } ( x, x ^ { \prime } ) = h _ { a } ^ { a ^ { \prime } } \left [ \left ( 1 - \frac { g _ { t t ; b } } { 4 g _ { t t } } \sigma ^ { b } \right ) \frac { \sigma _ { a ^ { \prime } } } { \sigma } - \frac { g _ { t t ; a ^ { \prime } } } { 4 g _ { t t } } \log \frac { \sigma ( x, x ^ { \prime } ) } { \lambda } + \dots \right ],$$

up to a factor which is evaluated at the charge position. In order to make the average, it is convenient to use Riemann normal coordinates x a . This procedure is extensively reviewed in [33] and we quote only the main ingredients. These coordinates are characterized in terms of the Synge world function as follows. The derivative σ a ′ = -x a and

$$2 \sigma = - 2 \delta _ { a b } x ^ { a } x ^ { b } = s ^ { 2 },$$

is the squared proper distance from the point x ′ . The parallel propagator in these coordinates is approximated by

$$h _ { b } ^ { a ^ { \prime } } = \delta _ { a } ^ { b } - \frac { 1 } { 6 } R _ { b ^ { \prime } c ^ { \prime } d ^ { \prime } } ^ { a ^ { \prime } } x ^ { b } x ^ { c } + O ( s ^ { 3 } ).$$

By denoting x a = Ω s a with Ω a some angular variables it follows that

$$\partial _ { a } G _ { s i n g } ( x, x ^ { \prime } ) = \frac { 2 \Omega _ { a } } { s } - \frac { g _ { t t ; b } } { 2 g _ { t t } } \Omega ^ { a } \Omega ^ { b } - \frac { g _ { t t ; a ^ { \prime } } } { 4 g _ { t t } } \log \frac { s ^ { 2 } } { \lambda } + \dots$$

In particular, one has that Ω a = δ ab Ω . The average along the surface b s = const is defined by

$$< \dots > = \frac { \int ( \dots ) d A } { \int d A },$$

with dA the area element of the surface. If this surface is parameterized by polar angles θ a then this and the relation x a = Ω s a implies that the metric is

$$d s ^ { 2 } = s ^ { 2 } h _ { a b } \Omega _ { i } ^ { a } \Omega _ { j } ^ { b } d \theta ^ { i } d \theta ^ { j },$$

where

$$\Omega _ { i } ^ { a } = \frac { \partial \Omega ^ { a } } { \partial \theta _ { i } }.$$

The area element is then calculated as

$$d A = s ^ { 2 } \left ( 1 - \frac { s ^ { 2 } } { 6 } R _ { a ^ { \prime } b ^ { \prime } } \Omega ^ { a } \Omega ^ { b } - \frac { s ^ { 3 } } { 1 2 } R _ { a ^ { \prime } b ^ { \prime } ; c ^ { \prime } } \Omega ^ { a } \Omega ^ { b } \Omega ^ { c } + O ( s ^ { 4 } ) \right ) d \Omega _ { n }.$$

Here d Ω n is the area element for the sphere S n , which satisfies

$$\Omega _ { n } = \int _ { S ^ { n } } d \Omega _ { n } = \frac { 2 \pi ^ { ( n + 2 ) / 2 } } { \Gamma ( \frac { n + 1 } { 2 } ) }.$$

The following identities take place [33]

$$\frac { 1 } { \Omega _ { n } } \int \Omega ^ { a } d \Omega _ { n } = 0,$$

$$\frac { 1 } { \Omega _ { n } } \int \Omega ^ { a } \Omega ^ { b } d \Omega _ { n } = \frac { \delta ^ { a b } } { n + 1 },$$

By use of these identities and the definition of the area element (8.208) it follows that

$$< \Omega ^ { a } > = - \frac { s ^ { 3 } } { 6 ( n + 1 ) ( n + 3 ) } \nabla ^ { a ^ { \prime } } R + O ( s ^ { 4 } )$$

$$< \Omega ^ { a } \Omega ^ { b } > = - \frac { 1 } { ( n + 1 ) } \left ( \delta ^ { a b } - \frac { s ^ { 2 } } { 6 ( n + 3 ) } R ^ { a ^ { \prime } b ^ { \prime } } + \frac { s ^ { 2 } } { 3 ( n + 1 ) ( n + 3 ) } R \delta ^ { a ^ { \prime } b ^ { \prime } } \right ) + O ( s ^ { 4 } )$$

$$\begin{array} { c } < \Omega ^ { a } > = - \frac { s ^ { 3 } } { 6 ( n + 1 ) ( n + 3 ) } \nabla ^ { a ^ { \prime } } R + O ( s ^ { 4 } ) \\ \\ < \Omega ^ { a } \Omega ^ { b } > = - \frac { 1 } { ( n + 1 ) } \left ( \delta ^ { a b } - \frac { s ^ { 2 } } { 6 ( n + 3 ) } R ^ { a ^ { \prime } b ^ { \prime } } + \frac { s ^ { 2 } } { 3 ( n + 1 ) ( n + 3 ) } R \delta ^ { a ^ { \prime } b ^ { \prime } } \right ) + O ( s ^ { 4 } ) \end{array}$$

Then the expression of the averaged singular Green function is

$$< \partial _ { a } G _ { s i n g } ( x, x ^ { \prime } ) > = \frac { 2 < \Omega _ { a } > } { s } - \frac { g _ { t t ; b } } { 2 g _ { t t } } < \Omega ^ { a } \Omega ^ { b } > - \frac { g _ { t t ; a ^ { \prime } } } { 4 g _ { t t } } \log \frac { s ^ { 2 } } { \lambda } + \dots$$

By taking into account the identities (8.209) it follows that in the limit s → 0 the surviving terms are

$$< \partial _ { a } G _ { s i n g } > ( x ^ { \prime } ) = - \frac { g _ { t t ; a ^ { \prime } } } { 6 g _ { t t } } - \frac { g _ { t t ; a ^ { \prime } } } { 4 g _ { t t } } \log \frac { s ^ { 2 } } { \lambda } + \dots$$

Thus, as in five dimensions [36], this average picks an extra term which is divergent as log s . This term can not be absorbed by a redefinition of the parameters of the problem. The nature of this divergence is different than the one of the previous section, it suggest that in odd dimensions the calculating is sensible to the size of the particle. The limit s → 0 gives a global divergence of the self force, while the divergence of the previous section is just asymptotic.

We would like to remark that we have considered the effect of the first term in (8.210) but it is not enough for avoiding the bad behavior at the infinite. Thus another interpretation is required. We turn to this point below.

## 9. Interpretation of the divergences

The most striking point of the results presented above is the behavior of the charge self force at the asymptotic region and at the horizon. For the black hole solution, it grows indefinitely when approaching both regions. For the wormhole, it diverges at the asymptotic boundary.

The physical origin of these divergencies is not easy to visualize. At first sight, it is not strange that the force does not vanish at the asymptotic region. In fact, when the geometry is not asymptotically flat, this situation may happen, as shown for four dimensions in [37]. When M = 1, on dimensional grounds, it is expected that the asymptotic self force should be proportional to

$$F _ { \infty } \sim \frac { e ^ { 2 } } { l }.$$

Now, the problem of the divergence at r →∞ may come from the expression (7.191)-(7.192) for the distance, namely

$$d ^ { a } ( x, x _ { 0 } ) = l \sqrt { 1 + \frac { ( \theta - \theta _ { 0 } ) ^ { 2 } } { \left ( \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s _ { 0 } } { 2 } \right ) ^ { 2 } } } ( s - s _ { 0 } ), \quad \ 0 < | \theta - \theta _ { 0 } | < \pi,$$

and the analogous formula for negative angle values. Recall that this distance is a fake one, but it reduces to the true radial distance when θ = θ 0 . The expressions that have been obtained for the self force are based on this fake distance. However, since the radial limit is correct, the resulting expressions coincide with the true self force when all the series is summed up. The choice of the fake distance was for simplicity, since the Fourier expansion for the potential is non tractable when the exact distance is considered (the true distance is derived in the Appendix). Nevertheless, an inspection of the formula d a ( x, x ′ ) shows the following potential numerical problem. This formula, as shown in the appendix, is similar to the real one when the points x and x 0 are almost on the same circle s ∼ s 0 . And coincides with real one when θ = θ 0 . The distance when the points are at different circles s &lt;&lt; s 0 or s &gt;&gt; s 0 is more complicated, and it turns out that the expressions differ considerably. Denote these distances as d c a ( x, x ′ ) and d r a ( x, x ′ ) respectively. Both distances are the same when evaluated on the line θ = θ 0 , but when θ = θ 0 + glyph[epsilon1] they may differ considerably. Moreover, their difference is a function of the radial coordinate | d r a ( x, x ′ ) -d c a ( x, x ′ ) | = f ( r, glyph[epsilon1] ). The series expansion presented above is convergent to the real solution but, due to the mentioned dependence in r , the convergence may be highly non uniform. In other words, one has that

$$| A _ { r e g } ^ { n } ( r ) - A _ { r e g } ( r ) | < \epsilon \quad \longleftrightarrow \quad n > n _ { 0 } ( r ),$$

with A n reg corresponding to a truncation of the series of A reg to order n . This means that, depending on the values of r , lower or higher orders may be required. We interpret that the divergence at r →∞ arises due to the fact that for large r values higher orders are required, and it is then an artifact of the truncation.

Another fact that suggest that this numerical problem is the cause of the divergence is the fact that when r and n takes large values the slope of the asymptote seem to grow indefinitely and it moves to the right of the graph ( s → 1). Arguably when n → ∞ this asymptote moves to r → ∞ ( s = 1) with infinite slope. Note that the point s = 1, since it corresponds to r → ∞ , is not a point in the manifold. Thus, we argue that when the series is summed up, the divergence is an artifact and the limit F r when r → ∞ is well defined. We were unable to overcome the problems described above, since the real distance leads to a Fourier expansion which is beyond our calculation technology.

It may be mentioned that the use of the approximate distance may be useful if a powerful summation formula were available allowing a closed analytical expression for the full potential A . This is the situation in four dimensions, as explained in [10] and references therein. But, if this formula do exist, we ignore it.

There is a further point to be discussed. When r &gt;&gt; lM , then one may consider M ∼ 0. For this situation, it follows by dimensional analysis that

$$F _ { r } = \frac { e ^ { 2 } } { l } f ( \frac { r } { l } ).$$

If the function f ( x ) is bounded, then when l → ∞ the self force goes to zero. This is the situation in four dimensional anti De Sitter spaces [37]. In this limit the geometry is flat and the self force vanishes. However, in three dimensions, such limit does not corresponds to a flat geometry, instead for M = 1 and J = 0 the BTZ metric (2.4) reduce to

$$g _ { 3 } = - d r ^ { 2 } + d t ^ { 2 } + r ^ { 2 } d \theta ^ { 2 },$$

in the limit l → ∞ . Note that r now is a 'time' coordinate and the compact radius is 'r-time' dependent. Thus the derivatives ∂ A r do not make sense as a self force, and perhaps this divergence is not signaling anything wrong. However, our opinion is that the behavior of the infinite is more likely a numerical artifact. This is suggested by the divergence at the horizon, for which we have no explanation. Thus, it is arguably a numerical problem what is causing the divergent behavior. Our results are of course non conclusive, and further research is mandatory to elucidate the real behavior at r &gt;&gt; l .

## 10. Summary

In the present work the particularities of the electrostatics of point charge in front of a BTZ black hole and wormhole were considered. This geometry is not flat neither asymptotically flat, and this results in several unexpected features. For instance, every radial functions for the electrostatic problems goes to zero at the asymptotic horizon. This complicates the boundary condition analysis, since there is no a clear criteria of which radial function discarding at the region between the charge position and the asymptotic infinite. We argue that, although both radial functions are well behaved, the right choice of the boundary condition is to choose the radial function which has the fastest decay when the radial

coordinate is large. This situation is the one which imitates better the electrostatic in flat space or in asymptotically flat spaces.

After discussing the boundary condition for the problem, we followed the program [9] -[16] in order to calculate the self force of a charge in front of the BTZ black hole horizon and compare it to the corresponding wormhole geometry. In addition we have considered the axiomatic approach of [38]. Unfortunately, both program seems to be harder to achieve. The regularized electrostatic potential is presented as a Fourier expansion, and we had no available summation formula for obtaining a closed analytical formula for such potential. As a consequence, we were forced to use an approximated geodesic distance which reduces to the real one when the two comparison points are on the same radial line. The use of such distance allows a series expansion which converges to the real vector potential, but arguably the convergence is highly non uniform. This is reflected in that the obtained diagrams for the self force, when truncated to a given order, have an asymptotic divergence. We interpret as fictitious, arising for the fact that a large number of orders are needed for large radial values. However, it is not clear which behavior should be expected in this geometry at asymptotic values, since it is not asymptotic flat neither it admits a flat limit.

There is a further possibility however. The boundary conditions that we have chosen are the ones that decay faster at the asymptotic region. This condition imitates ordinary electrostatics and is the one that is more comfortable to our intuition. However, there may be some symmetries of the BTZ space time similar to those in higher dimensional black holes [39]. We do not know how to implement this symmetry in the BTZ geometry and how is reflected in the Green function, but perhaps the result does not respect this boundary condition. After all, what the present work is showing is that flat geometry based intuitions may lead to wrong results. In any case our results are not conclusive about this issue. A further investigation about this point may be of special interest for the future.

Acknowledgements: The authors are supported by CONICET (Argentina).

## Appendix

## Exact BTZ geodesic distance

Our next task is to find the geodesic distance d r, r , θ, θ ( ′ ′ ) between two arbitrary points of the space time. This distance should reduce to (7.181) or (7.182) when θ ′ → θ . The geodesic curve joining two points x and x ′ on the BTZ geometry is describe by the Euler-Lagrangian equations derived with the following action

$$d ( { \mathbf x } _ { 1 }, { \mathbf x } _ { 2 } ) = \int _ { { \mathbf x } _ { 1 } } ^ { { \mathbf x } _ { 2 } } \sqrt { \frac { l ^ { 2 } } { r ^ { 2 } - r _ { h } } \left ( \frac { d r } { d \lambda } \right ) ^ { 2 } + r ^ { 2 } \left ( \frac { d \theta } { d \lambda } \right ) ^ { 2 } } d \lambda,$$

The geodesic parameter can be chosen as λ = . The corresponding lagrangian is θ

$$\mathcal { L } ( \dot { r }, r ) = \sqrt { \frac { l ^ { 2 } \dot { r } ^ { 2 } } { r ^ { 2 } - r _ { h } ^ { 2 } } + r ^ { 2 } }.$$

The value of (A.1) evaluated at the geodesic curve is by definition the geodesic distance. Since L does not depend explicitly on the 'time' parameter λ = θ there exists a conserved quantity, the hamiltonian

$$H = \frac { \partial \mathcal { L } } { \partial \dot { r } } \dot { r } - \mathcal { L } = \frac { - r ^ { 2 } } { \sqrt { \frac { l ^ { 2 } \dot { r } ^ { 2 } } { r ^ { 2 } - r _ { h } ^ { 2 } } + r ^ { 2 } } }.$$

In terms of this constant, the lagrangian evaluated in the physical trajectory is

$$\mathcal { L } ^ { 2 } = g _ { r r } \dot { r } ^ { 2 } + g _ { \theta \theta } ^ { - 1 } H ^ { 2 } \mathcal { L } ^ { 2 },$$

from where it follows that

$$\mathcal { L } ^ { 2 } = g _ { r r } g _ { \theta \theta } \dot { r } ^ { 2 } \left ( g _ { \theta \theta } - H ^ { 2 } \right ) ^ { - 1 }.$$

In these terms the geodesic distance is given as

$$d ( { \mathbf x } _ { 1 }, { \mathbf x } _ { 2 } ) = \int _ { \lambda _ { 1 } } ^ { \lambda _ { 2 } } \sqrt { \frac { g _ { r r } g _ { \theta \theta } } { ( g _ { \theta \theta } - H ^ { 2 } ) } } \, \dot { r } \, d \lambda$$

and since ˙ rdλ = dr , then by taking into account the explicit expression for the metric tensor, it follows that √

$$d ( { \mathbf x } _ { 1 }, { \mathbf x } _ { 2 } ) = l \int _ { r _ { 1 } } ^ { r _ { 2 } } \frac { r \, d r } { \sqrt { ( r ^ { 2 } - r _ { h } ^ { 2 } ) \, ( r ^ { 2 } - H ^ { 2 } ) } } = l \log \left [ \frac { \sqrt { r _ { 2 } ^ { 2 } - r _ { h } ^ { 2 } } + \sqrt { r _ { 2 } ^ { 2 } - H ^ { 2 } } } { \sqrt { r _ { 1 } ^ { 2 } - r _ { h } ^ { 2 } } + \sqrt { r _ { 1 } ^ { 2 } - H ^ { 2 } } } \right ].$$

The distance (A.4) is not completely determined, unless the value of H is given in terms of the two positions x and x ′ . This is achieved as follows. From (A.2) and (A.3) it is deduced that

$$\left ( \frac { d r } { d \theta } \right ) ^ { 2 } = \frac { \mathcal { L } ^ { 2 } } { g _ { r r } g _ { \theta \theta } } \left ( g _ { \theta \theta } - H ^ { 2 } \right ) \left ( \frac { g _ { \theta \theta } } { H \mathcal { L } } \right ) ^ { 2 } = \frac { r ^ { 2 } } { H ^ { 2 } l ^ { 2 } } \left ( r ^ { 2 } - H ^ { 2 } \right ) \left ( r ^ { 2 } - r _ { h } ^ { 2 } \right ),$$

and since ˙ r 2 &gt; 0 and r &gt; r h , then H &lt; r 2 2 . From the last equation it follows that

$$\theta _ { 2 } - \theta _ { 1 } = \pm l H \int _ { r _ { 1 } } ^ { r _ { 2 } } \frac { d r } { r \sqrt { r ^ { 2 } - H ^ { 2 } } \sqrt { r ^ { 2 } - r _ { h } ^ { 2 } } }.$$

The last integration is elementary, the result is

$$\mu _ { 2 } & = \log \left [ \frac { r _ { 1 } } { r _ { 2 } } \left ( \frac { r _ { h } \sqrt { r _ { 2 } ^ { 2 } - H ^ { 2 } } + L \sqrt { r _ { 2 } ^ { 2 } - r _ { h } ^ { 2 } } } { r _ { h } \sqrt { r _ { 1 } ^ { 2 } - H ^ { 2 } } + L \sqrt { r _ { 1 } ^ { 2 } - r _ { h } ^ { 2 } } } \right ) \right ].$$

The expression just obtained can be worked out further by defining

$$e ^ { \alpha } = e ^ { \theta _ { 2 } - \theta _ { 1 } } = \frac { \sqrt { H ^ { - 2 } - r _ { 2 } ^ { - 2 } } + \sqrt { r _ { h } ^ { - 2 } - r _ { 2 } ^ { - 2 } } } { \sqrt { H ^ { - 2 } - r _ { 1 } ^ { - 2 } } + \sqrt { r _ { h } ^ { - 2 } - r _ { 1 } ^ { - 2 } } }.$$

This definition is equivalent to

$$e ^ { \alpha / 2 } \sqrt { H ^ { - 2 } - r _ { 1 } ^ { - 2 } } - e ^ { - \alpha / 2 } \sqrt { H ^ { - 2 } - r _ { 2 } ^ { - 2 } } = e ^ { - \alpha / 2 } \sqrt { r _ { h } ^ { - 2 } - r _ { 2 } ^ { - 2 } } - e ^ { \alpha / 2 } \sqrt { r _ { h } ^ { - 2 } - r _ { 1 } ^ { - 2 } }.,$$

and by taking the square of both members it follows after some algebra that

$$\tilde { \Lambda } = \tilde { \eta } + ( H ^ { - 2 } - r _ { h } ^ { - 2 } ) \cosh \alpha,$$

where the following quantities

$$\tilde { \Lambda } = \sqrt { H ^ { - 2 } - r _ { 2 } ^ { - 2 } } \sqrt { H ^ { - 2 } - r _ { 1 } ^ { - 2 } }, \quad \tilde { \eta } = \sqrt { r _ { h } ^ { - 2 } - r _ { 1 } ^ { - 2 } } \sqrt { r _ { h } ^ { - 2 } - r _ { 2 } ^ { - 2 } },$$

has been introduced. On the other hand (A.4) implies that

$$e ^ { d / l } = \frac { \theta ^ { \quad \cdot } } { \sqrt { r _ { 1 } ^ { 2 } - r _ { h } ^ { 2 } } + \sqrt { r _ { 2 } ^ { 2 } - H ^ { 2 } } } { \sqrt { r _ { 1 } ^ { 2 } - r _ { h } ^ { 2 } } + \sqrt { r _ { 1 } ^ { 2 } - H ^ { 2 } } },$$

and procedure analogous to the one made above shows that

$$\Lambda & = \eta + ( r _ { h } ^ { \, 2 } - H ^ { 2 } ) \cosh \frac { d } { l }. & & ( \text{A.7} )$$

where in this case

$$\Lambda = \sqrt { r _ { 1 } ^ { 2 } - H ^ { 2 } } \sqrt { r _ { 2 } ^ { 2 } - H ^ { 2 } }, \quad \eta = \sqrt { r _ { 1 } ^ { 2 } - r _ { h } ^ { 2 } } \sqrt { r _ { 2 } ^ { 2 } - r _ { h } ^ { 2 } }.$$

By noticing the relations

$$\eta = r _ { 1 } r _ { 2 } r _ { h } ^ { 2 } \tilde { \eta }, \quad \Lambda = r _ { 1 } r _ { 2 } H ^ { 2 } \tilde { \Lambda },$$

then comparison between (A.6) and (A.7) shows that

$$r _ { h } ^ { 2 } \cosh \frac { d } { l } = r _ { 1 } r _ { 2 } \cosh \alpha - \eta.$$

From the last expression the explicit geodesic distance d can be obtained, the result is

$$d ( { \mathbf x } _ { 1 }, { \mathbf x } _ { 2 } ) = l \cosh ^ { - 1 } \left [ \frac { r _ { 1 } r _ { 2 } } { r _ { h } ^ { 2 } } \cosh ( \theta _ { 2 } - \theta _ { 1 } ) - \sqrt { \left ( \frac { r _ { 1 } } { r _ { h } } \right ) ^ { 2 } - 1 } \sqrt { \left ( \frac { r _ { 2 } } { r _ { h } } \right ) ^ { 2 } - 1 } \right ],$$

This is the formula (7.189) obtained in the text.

## Approximate geodesic distance

In the following the attention will be restricted to BTZ geometries with mass values M = 1. Based on the arguments given above and from the fact that the calculation of the self force requires to take the limit θ → 0, our next task is to find an approximated distance function d a ( x, x ′ ) which is equal to d x, x ( ′ ) when θ → 0 and satisfy the mandatory periodicity conditions. The upper index a enforce the fact that the distance is not exact and this notation will be used repeatedly in the following. At first sight there is a variety of candidates for d a ( x, x ′ ). But one that is physically motivated is the following. Consider the spatial BTZ metric (7.180) in the coordinates ( s, θ ). The geodesic distance is

$$d = l \, \int \sqrt { s ^ { \prime 2 } + \cosh ^ { 2 } s } \, d \theta,$$

with s ′ the derivate of s θ ( ) with respect to θ . The conserved quantity H related to the θ independence on the lagrangian is expressed in these coordinates as

$$H = \frac { l \cosh ^ { 2 } s } { \sqrt { s ^ { \prime 2 } + \cosh ^ { 2 } s } }.$$

$$l$$

It is important for the following to discuss the physical significance of this quantity. When the geodesic line is close to a circle, then s ′ ∼ 0 and H is essentially the radius of the circle r = l cosh s . Instead, when the geodesic is close to a radial line, then s ′ →∞ and H → 0 in this situation. The definition of H implies that

$$s ^ { \prime 2 } = \frac { l ^ { 2 } \cosh ^ { 4 } s } { H ^ { 2 } } - \cosh ^ { 2 } s,$$

and for nearly radial lines one can neglect the second term in (A.9) since H &lt;&lt; 1, thus concluding that

$$s ^ { \prime 2 } \sim \frac { l ^ { 2 } \cosh ^ { 4 } s } { H ^ { 2 } }.$$

For nearly circular geodesics one has that

$$H ^ { 2 } \sim ( l ^ { 2 } - \epsilon ^ { 2 } ) \cosh ^ { 2 } s,$$

with glyph[epsilon1] 2 &lt;&lt; l 2 , and within this approximation (A.9) becomes

$$s ^ { \prime 2 } & \sim \frac { \epsilon ^ { 2 } } { 2 l ^ { 2 } } \cosh ^ { 2 } s.$$

We have explicitly checked that the equation (A.10) leave to the simplest expression for d a ( x, x ′ ). The equation (A.10) implies that

$$\frac { 1 } { c } \frac { d s } { d \theta } = \pm \cosh ( s ),$$

c being a constant related to glyph[epsilon1] 2 whose value is fixed by the initial conditions of the problem. The equation (A.11) can be integrated in elementary form to give

1

[

]

θ

-

θ

s

s

0

0

=

±

-

1

-

1

c

tan tanh

2

-

tan tanh

2

.

This shows that the constant c is given in terms of the initial and final positions by

$$c = \frac { \left [ \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s _ { 0 } } { 2 } \right ] } { ( \theta - \theta _ { 0 } ) }.$$

Inserting (A.11) and (A.12) in (A.8) give a result the following expression for the approximated geodesic distance

$$d ^ { a } ( x, x _ { 0 } ) = l \sqrt { 1 + \frac { ( \theta - \theta _ { 0 } ) ^ { 2 } } { \left ( \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s _ { 0 } } { 2 } \right ) ^ { 2 } } } ( s - s _ { 0 } ), \quad \ 0 < | \theta - \theta _ { 0 } | < \pi,$$

$$d ^ { a } ( x, x _ { 0 } ) = l \sqrt { 1 + \frac { ( 2 \pi - \theta + \theta _ { 0 } ) ^ { 2 } } { \left ( \tan ^ { - 1 } \tanh \frac { s } { 2 } - \tan ^ { - 1 } \tanh \frac { s _ { 0 } } { 2 } \right ) ^ { 2 } } } ( s - s _ { 0 } ), \quad \pi < | \theta - \theta _ { 0 } | < 2 \pi.$$

This is the expression (7.191)-(7.192) found in the text.

## References

- [1] Wheeler, John Archibald (1963). Geometrodynamics. New York: Academic Press. LCCN 62013645.
- [2] A. Vilenkin, Phys. Rev. D 20 , 373 (1979).
- [3] A. G. Smith and C. M. Will, Phys. Rev. D 22 , 1276 (1980).
- [4] B. Linet, Phys. Rev. D 33 , 1833 (1986).
- [5] A. Vilenkin, Phys. Rev. D 23 , 852 (1981); J. R. Gott III, Astrophys. J. 288 , 422 (1985); W. A. Hiscock, Phys. Rev. D 31 , 3288 (1985).
- [6] N. R. Khusnutdinov and I. V. Bakhmatov Phys. Rev. D 76 , 124015 (2007).
- [7] M. S. Morris and K. S. Thorne, Am. J. Phys. 56 , 395 (1988).
- [8] M. Visser, Lorentzian Wormholes (AIP Press, New York, 1996).
- [9] E. Rub´ ın de Celis, O. P. Santill´n and C. Simeone, Phys. Rev. a D 86 , 124009 (2012).
- [10] E. Rubin de Celis, O. Santill´n and C. Simeone Phy. Rev. D a 86(12) (2013) 124009
- [11] T. Drivas and S. Gralla, Class. Quant. Grav. 28, 145025 (2011).
- [12] S. Isoyama and E. Poisson, Class. Quant. Grav. 29, 155012 (2012).
- [13] B. Boisseau, C. Charmousis and B. Linet Class. Quant. Grav. 13 (1996) 1797
- [14] B. Linet, Class.Quant.Grav. 16, 2947 (1999).
- [15] B. Linet, Class.Quant.Grav. 17, 4661 (2000).
- [16] N. R. Khusnutdinov, A. A. Popov and L. N. Lipatova, Class. Quant. Grav. 27, 215012 (2010).
- [17] S. Deser, R. Jackiw, and G. t Hooft, Ann. Phys. 152 (1984) 220; S. Deser and R. Jackiw, Commun. Math. Phys. 118 (1988) 495; G. t Hooft, Commun. Math. Phys. 117 (1988) 685.
- [18] E. Witten, Nucl. Phys. B311 (1988) 46; E. Witten, Nucl. Phys. B323 (1989) 113.
- [19] M. Banados, C. Teitelboim and J. Zanelli, Phys. Rev. Lett., 69, 1849 (1992)
- [20] M. Banados, M. Henneaux, C. Teitelboim, and J. Zanelli, Phys. Rev. D48, 1506 (1993).
- [21] S. Carlip Class. Quant. Grav. 22 (2005) R85.
- [22] S. Deser, R. Jackiw and G.t Hooft, Ann. Phys. (NY) 152, 220 (1984). S. Deser and R. Jackiw, Ann. Phys. (NY) 153, 405 (1984). A. Staruszkiewicz, Acta Phys. Pol. 24, 734 (1963). J. Gott and M. Alpert, Gen. Rel.and Grav. 16, 3 (1984).

- [23] E. T. Whittaker and G. N. Watson, A Course of Modern Analysis (Cambridge University Press; 4th edition, January 2, 1927).
- [24] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Dover, New York, 1972).
- [25] Ian N. Sneddon Special Functions of Mathematical Physics and Chemistry Oliver and Boyd 1961
- [26] Gradshteyn and Ryzhik's Table of Integrals, Series, and Products Alan Jeffrey and Daniel Zwillinger (eds.) Seventh edition (2007).
- [27] J. Synge Relativity: the general theory. Amsterdam: North-Holland 1960.
- [28] B. Allen and T. Jacobson Commun. Math. Phys. 103 (1986) 669.
- [29] S. Detweiler and B. F. Whiting, Phys. Rev. D 67 , 024025, (2003).
- [30] M. Casals, E. Poisson and I. Vega, Phys. Rev. D 86 , 064033, (2012).
- [31] B. S. DeWitt and R. W. Brehme, Radiation damping in a gravitational eld, Ann. Phys. (N.Y.) 9, 220 (1960).
- [32] F. G. Friedlander, The wave equation on a curved spacetime (Cambridge University Press, Cambridge, 1975).
- [33] Eric Poisson and Adam Pound and Ian Vega, 'The Motion of Point Particles in Curved Spacetime', Living Rev. Relativity 14, (2011), 7. URL: http://www.livingreviews.org/lrr-2011-7
- [34] B. DeWitt and R. Brehme Annals Phys. 9 (1960) 220.
- [35] J. Hadamard, Lectures on Cauchy's Problem in Linear Partial Differential Equations (Dover Publications, New York, 1952).
- [36] M.. Beach, E. Poisson and B. Nickel 'Self-force on a charge outside a five-dimensional black hole' arXiv:1404.1031.
- [37] J. Kuchar, E. Poisson and I. Vega 'Electromagnetic self-force on a static charge in Schwarzschildde Sitter spacetimes' arXiv:1307.8342.
- [38] T. C. Quinn and R. M. Wald Phys. Rev. D 56, 3381 (1997).
- [39] V. Frolov and A. Zelnikov 'Bi-conformal symmetry and static Green functions in the Schwarzschild-Tangherlini spacetimes' arXiv:1412.6719.
Y. Herrera V. Hurovich O. Santillán C. Simeone
Aspects of electrostatics in BTZ geometries
In the present paper the electrostatic of charges in non rotating BTZ black
hole and wormhole space times is studied. In particular, the self force of a
point charge in the geometry is characterized analitically. The differences
between the self force in both cases is a theoretical experiment for
distinguishing both geometries, which otherwise are locally indistinguishable.
This idea was applied before to four and higher dimensional black holes by the
present and other authors. However, the particularities of the BTZ geometry
makes the analysis considerable more complicated than usual electrostatic in a
flat geometry, and its even harder than its four dimensional counterparts.
First, the BTZ space times are not asymptotically flat but instead
asymptotically AdS. In addition, the relative distance $d(r,r+1)$ between two
particles located at a radius $r$ and $r+1$ in the geometry tends to zero when
$r\to\infty$. This behavior, which is radically different in a flat geometry,
changes the analysis of the asymptotic conditions for the electrostatic field.
In addition, there are no summation formulas that allow a closed analytic
expression for the self force. We find a method to calculate such force in
series, and the resulting expansion is convergent to the real solution.
However, we suspect that the convergence is not uniform. In other, for points
that are far away from the black hole the calculation of the force requires
higher order summation. These subtleties are carefully analyzed in the paper,
and it is shown that they lead to severe problems when calculating the self
force for asymptotic points in the geometry.