go to next pagego upgo to previous page

2.2 Chandrasekhar-Sasaki-Nakamura transformation

As seen from the asymptotic behaviors of the radial functions given in Equations (25View Equation) and (26View Equation), the Teukolsky equation is not in the form of a canonical wave equation near the horizon and infinity. Therefore, it is desirable to find a transformation that brings the radial Teukolsky equation into the form of a standard wave equation.

In the Schwarzschild case, Chandrasekhar found that the Teukolsky equation can be transformed to the Regge-Wheeler equation, which has the standard form of a wave equation with solutions having regular asymptotic behaviors at horizon and infinity [11Jump To The Next Citation Point]. The Regge-Wheeler equation was originally derived as an equation governing the odd parity metric perturbation [47Jump To The Next Citation Point]. The existence of this transformation implies that the Regge-Wheeler equation can describe the even parity metric perturbation simultaneously, though the explicit relation of the Regge-Wheeler function obtained by the Chandrasekhar transformation with the actual metric perturbation variables has not been given in the literature yet.

Later, Sasaki and Nakamura succeeded in generalizing the Chandrasekhar transformation to the Kerr case [50Jump To The Next Citation Point51Jump To The Next Citation Point]. The Chandrasekhar-Sasaki-Nakamura transformation was originally introduced to make the potential in the radial equation short-ranged, and to make the source term well-behaved at the horizon and at infinity. Since we are interested only in bound orbits, it is not necessary to perform this transformation. Nevertheless, because its flat-space limit reduces to the standard radial wave equation in the Minkowski spacetime, it is convenient to apply the transformation when dealing with the post-Minkowski or post-Newtonian expansion, at least at low orders of expansion.

We transform the homogeneous Teukolsky equation to the Sasaki-Nakamura equation [5051], which is given by

( 2 ) -d--- F(r) -d-- U(r) Xlmw = 0. (50) dr*2 dr*
The function F(r) is given by
j D F(r) = --,r-------, (51) j r2 + a2
where
c1 c2 c3 c4 j = c0 + -- + -2 + --3 + -4, (52) r r r r
with
c0 = -12iwM + c(c + 2) - 12aw(aw - m), c1 = 8ia[3aw - c(aw - m)], c = -24iaM (aw - m) + 12a2[1 - 2(aw - m)2], 2 (53) c3 = 24ia3(aw - m) - 24M a2, c4 = 12a4.
The function U(r) is given by
--DU1----- 2 -DG,r-- U (r) = (r2 + a2)2 + G + r2 + a2- FG, (54)
where
2(r---M--) ---rD----- G = - r2 + a2 + (r2 + a2)2, (55) [( ) ( )] D2- b,r j,r b,r U1 = V + b 2a + D - j a + D , (56) ,r Kb-- 6D- a = - i D2 + 3iK,r + c + r2 , (57) ( 2D ) b = 2D - iK + r- M - --- . (58) r

The relation between Rlmw and Xlmw is

[( ) ] 1 b,r b Rlmw = j- a + D-- xlmw - D-xlmw,r , (59)
where xlmw = XlmwD/(r2 + a2)1/2. Conversely, we can express Xlmw in terms of Rlmw as
( ) Xlmw = (r2 + a2)1/2 r2J- J- 1-Rlmw , (60) r2
where J- = (d/dr) - i(K/D).

If we set a = 0, this transformation reduces to the Chandrasekhar transformation for the Schwarzschild black hole [11]. The explicit form of the transformation is

D ( d ) r2 ( d ) Rlmw = -- ----+ iw -- ----+ iw rXlmw, (61) c0 ( dr* ) D ( dr* ) r5 d r2 d Rlmw Xlmw = -- --*-- iw -- --*-+ iw --2--, (62) D dr D dr r
where c0, defined in Equation (53View Equation), reduces to c0 = (l - 1)l(l + 1)(l + 2)- 12iM w. In this case, the Sasaki-Nakamura equation (50View Equation) reduces to the Regge-Wheeler equation [47], which is given by
( d2 ) --*2 + w2 - V (r) Xlw(r) = 0, (63) dr
where
( 2M ) (l(l + 1) 6M ) V (r) = 1 - ---- ----2---- --3- . (64) r r r
As clear from the above form of the equation, the lowest order solutions are given by the spherical Bessel functions. Hence it is intuitively straightforward to apply the post-Newtonian expansion to it. Some useful techniques for the post-Newtonian expansion were developed for the Schwarzschild case by Poisson [43] and Sasaki [49Jump To The Next Citation Point].

The asymptotic behavior of the ingoing wave solution Xin which corresponds to Equation (20View Equation) is

{ Aref eiwr* + Ainc e- iwr* for r* --> oo , Xinlmw --> lmw * lmw (65) Atrlamnwse-ikr for r* --> - oo .
The coefficients Ainc, Aref, and Atrans are related to Binc, Bref, and Btrans, defined in Equation (20View Equation), by
1 Binlmcw = - ---2Ailnmcw, (66) 4w 2 Bref = - 4w--Aref , (67) lmw c0 lmw 1 Btlrmanws = -----Atrlmawns, (68) dlmw
where
------ d = V~ 2M r [(8- 24iM w - 16M 2w2)r2 lmw + + 2 +(12iam - 16M + 16amM w + 24iM w)r+ - 4a2m2 - 12iamM + 8M 2]. (69)

In the following sections, we present a method of post-Newtonian expansion based on the above formalism in the case of the Schwarzschild background. In the Kerr case, although a post-Newtonian expansion method developed in previous work [52Jump To The Next Citation Point58Jump To The Next Citation Point] was based on the Sasaki-Nakamura equation, we will not present it in this paper. Instead, we present a different formalism, namely the one developed by Mano, Suzuki, and Takasugi which allows us to solve the Teukolsky equation in a more systematic manner, albeit very mathematical [33Jump To The Next Citation Point]. The reason is that the equations in the Kerr case are already complicated enough even if one uses the Sasaki-Nakamura equation, so that there is not much advantage in using it. In contrast, in the Schwarzschild case, it is much easier to obtain physical insight into the role of relativistic corrections if we deal with the Regge-Wheeler equation.



go to next pagego upgo to previous page