A Analytic approximation method for \omega_{*}\ll 1 and a_{*}\ll 1

Bulk scalar emission from a rotating black hole pierced by a tense brane


We study the emission of scalar fields into the bulk from a six-dimensional rotating black hole pierced by a 3-brane. We determine the angular eigenvalues in the presence of finite brane tension by using the continued fraction method. The radial equation is integrated numerically, giving the absorption probability (graybody factor) in a wider frequency range than in the preexisting literature. We then compute the power and angular momentum emission spectra for different values of the rotation parameter and brane tension, and compare their relative behavior in detail. As is expected from the earlier result for a nonrotating black hole, the finite brane tension suppresses the emission rates. As the rotation parameter increases, the power spectra are reduced at low frequencies due to the smaller Hawking temperature and are enhanced at high frequencies due to superradiance. The angular momentum spectra are enhanced over the whole frequency range as the rotation parameter increases. The spectra and the amounts of energy and angular momentum radiated away into the bulk are thus determined by the interplay of these effects.

04.50.-h, 04.70.Dy

I Introduction

Braneworld models with large extra dimensions (1); (2); (3) bring us an interesting possibility to address the hierarchy problem by lowering the fundamental scale of gravity down to order of TeV. It has been argued in the context of TeV scale gravity that mini black holes might be created through high-energy particle collision at future colliders (4). Much effort has been directed towards a theoretical understanding of the black hole formation at TeV energies (e.g., (5)). After their production, the black holes will decay via Hawking radiation (6). This process provides a window to probe high-energy physics, gravity at small distances, and properties of extra dimensions, which motivates recent extensive studies on this topic. A nonexhaustive sampling of the literature can be found in Refs. (7); (8); (9); (11); (10); (12); (13); (14); (15); (16); (17); (18); (19); (20); (21); (22); (23); (24); (25); (26); (27). For a review see Ref. (28). Most of the related work to date has ignored the effect of brane tension, treating black holes as “isolated” ones (see, however, Refs. (29); (30); (31); (32) for the effects of self-gravity of branes).

It is in general very difficult to obtain a black hole solution localized on a brane with finite tension because tension curves the brane as well as the bulk (cf. Refs. (33); (34); (35); (36)). However, codimension-2 branes exceptionally allow for a simple construction of localized black holes thanks to their special property; starting from the Myers-Perry solution (37) one rescales the polar angle around a symmetry axis as and then the brane tension is proportional to the deficit angle . In this way both nonrotating (38) and rotating (39) black holes on codimension-2 branes have been constructed. Following the work of (38), Hawking evaporation (40) and the quasi-normal modes for bulk scalars (41); (42) and fermions (43) have been investigated in the nonrotating background, showing that the finite brane tension modifies the standard result derived assuming negligible tension.

In this paper, we shall consider a six-dimensional rotating black hole pierced by a tense 3-brane and discuss the emission of massless scalar fields into the bulk. We intend to shed light on the spin-down phase in the life of a black hole, which is often neglected in the literature but could be of some significance. (In fact, a rotating black hole does not necessarily spin-down to zero, but evolves toward a nonzero angular momentum (44); (45).) Ignoring the brane tension, very recently Creek et al. studied the emission of scalars in the bulk in a higher-dimensional rotating black hole background (25). They employed matching techniques to obtain an analytic solution to the scalar field equation, which is a good approximation in the low-energy () and slow-rotation () regime, where is the energy of the emitted particle, is the black hole horizon radius, and is the rotation parameter. In the present paper, with the help of numerical computations we are able to handle the intermediate regime ( and ), and thus we not only include the effect of the finite tension but also extend the range of validity of (25).

This paper is organized as follows. In the next section we give a quick review of the rotating black hole solution on a codimension-2 brane. In Sec. III we present separated equations of motion for a massless scalar field and determine angular eigenvalues in the presence of the deficit angle. Then in Sec. IV the radial equation is solved numerically to give the power and angular momentum emission spectra. Finally we summarize our conclusions in Sec. V. Appendix contains the analytic calculation of the absorption probability, which complements the numerical results presented in the main text.

Ii A rotating black hole on a codimension two brane

We begin with a brief review of the rotating black hole solution on a codimension-2 brane. (For further detail see Refs. (38); (39).) The solution shares some properties with the Myers-Perry black hole (37). We are considering the models with five spatial dimensions, and so the rotation group is . The number of Casimirs (i.e., the number of mutually commuting elements of the group) is equal to rank. Hence, we have two axes of rotation associated with two angular momenta. However, in the present article we will be focusing on the special but simple case of a single rotation parameter with the angular momentum pointing along the brane. This is indeed an interesting case from the phenomenological point of view, because the black hole formed by the collision of two particles confined to the brane will have a single rotation parameter. The exact metric that describes such a rotating black hole is given by (39)




The coordinate ranges are , and . The parameter is related to the brane tension as


where is the six-dimensional fundamental scale. We assume that . When , the above metric reduces to the usual Myers-Perry solution with a single rotation parameter in six dimensions (37). When the solution is asymptotically conical. The parameters and denote the specific mass and angular momentum, respectively, related to the ADM mass and angular momentum of the black hole as


where is the area of a unit -sphere. Note here that the effect of the deficit angle is separated out from the definition of the area. The black hole horizon radius follows from . For later purpose it is convenient to define the dimensionless measure of the angular momentum . Since the sign flip simply changes the direction of rotation, in what follows we will assume without any loss of generality.

Note that has a root for arbitrary . This should be contrasted with the four-dimensional Kerr black hole, which has an upper bound on the rotation parameter. This point will become clearer by considering the black hole temperature. The generator of the horizon is parallel to the Killing vector , and the surface gravity of the black hole is given by . Here is the angular velocity of the horizon:


Then, the Hawking temperature is given by (6)


and one sees that the temperature never attains zero. This means the black hole is manifestly nondegenerate for any value of .

An interesting effect of the brane tension perhaps one is immediately aware of is the rescaling of the gravitational scale , as is observed in Eq. (4). This modifies the scenario of the black hole production and evaporation as follows (38); (39); (40). Suppose that a black hole is formed through a particle collision and its mass is fixed by the accelerator. The upper limit on the impact parameter is given by . Since , for fixed the horizon radius is larger than the naive estimate, , which is called the “lightning rod” effect (38); (39). In the present situation producing black holes will become easier due to this effect. Since the black hole formed in this way has the angular momentum , we have , and so the angular momentum can also be larger than naively expected, . Now the time scale of mass loss is estimated from (horizon area), yielding . Thus the lifetime of the black hole will become longer. Similarly, the time scale of angular momentum loss is estimated roughly from . This gives , and hence the time scale will also become larger.

The above mentioned modifications are caused by the rescaling of the gravitational scale. In the rest of the paper we separate this effect and discuss the impact of finite brane tension on Hawking emission with the horizon radius fixed.

Let us finally remark that in models with extra dimensions the bulk will be compact. Throughout the paper we assume that the horizon size of the black hole is smaller than the typical compactification scale. With this, the above solution provides a nice description of a six-dimensional black hole localized on a 3-brane.

Iii Scalar wave equation in a rotating black hole background

We now turn to discuss Hawking emission of massless scalar fields. To this end we solve the equation of motion for a massless scalar field in the rotating black hole background (1). The governing equation is given by


where . We can separate the above equation by assuming the ansatz


where The radial equation is




and the angular equations are given by


Here and are separation constants. The above equations of motion are the same as those derived in Ref. (16) except that in the last equation (13) there appears “.” We immediately see that


As is emphasized in Ref. (40), the eigenvalues are coupled and therefore all of the eigenvalues will be modified when . Since the eigenvalues and will be dependent on , we can restrict to without any loss of generality. Negative values of can be treated similarly.

Figure 1: Eigenvalue for selected values of and . Solid and dashed lines refer to the series expansion (28) truncated at 7th order, while points and crosses show the numerical results without relying on the small- expansion.

We proceed to determine the eigenvalue . Performing the change of the variable and function as , , we have


The boundary conditions are given by


As in (42), we first consider the special case in which with and hence is an integer. In this case it is easy to find Thus the eigenvalue can be written in terms of an integer as


where and the absolute value signs are inserted for clarity. Although Eq. (17) was derived for , this result can be generalized to arbitrary  (42). We have confirmed numerically that Eq. (17) indeed holds for noninteger values of .

When , is independent of (and as a result is also independent of ). For the tensional case, however, depends on (and hence does ). From Eq. (17) we see that the brane tension increases the eigenvalue relative to the tensionless case.

To determine the eigenvalue , we exploit the continued fraction method developed originally by Leaver (46). Following (47), we write with . Then, obeys

subject to the regularity boundary conditions


We seek for a series solution in the form of


Substituting this to Eq. (LABEL:calSeq) we obtain the three-term recursion relation




The continued fraction equation for the eigenvalue is then given by


The boundary conditions are automatically satisfied by the series solution (21).

We expand the eigenvalue in powers of around :


In order for the series to converge in the limit , we require that it has a finite number of terms. Imposing for some integer and identifying , we obtain


for . Here we have a restriction for the integer :


To determine the coefficients in (28) it is convenient to use the th inversion of (27):


With some manipulation we find and


One can confirm that Eqs. (32)–(34) correctly reproduce the known result in the case of  (47); (25).2 Although the expressions are too lengthy, higher order coefficients can be obtained easily.

One can instead integrate Eq. (11) numerically to determine . The numerical computation has an advantage that it does not rely on the small- expansion. In Fig. 1 we compare the series expansion truncated at 7th order with the numerical result that is free from any approximation. It can be seen from this that the analytic result can reproduce the numerical computation remarkably well even for . Therefore, in the next section we will safely use the analytic approximation for . We thus avoid numerical determination of the eigenvalues and so the rest of the problem simply reduces to solving the ordinary differential equation for the radial mode function.

We remark here that while the brane tension does not affect the mode, the eigenvalue (for ) increases as decreases (and so the tension increases), as can be seen from Fig. 1. This behavior has been observed for a nonrotating black hole in (40), and now it turns out that the similar thing generally holds for a rotating one.

Iv Power and angular momentum emission spectra

Figure 2: Absorption probabilities for , , and various .
Figure 3: Superradiance modes for and . The absorption probabilities are negative for , where for and for .

We are going to solve the radial equation (9) to compute the emission spectra. To do so we first specify the asymptotic form of the solution close to the horizon and far away from it. In terms of the new coordinate defined by and the function , Eq. (9) can be rewritten as


where and


It is easy to see that as and . Keeping this in mind the asymptotic solutions are found to be




and is the tortoise-like coordinate defined by


We choose the boundary condition , i.e., we impose that no outgoing modes exist near the horizon. The absorption probability is then given by


Note the explicit dependence on the angular eigenvalue of the absorption probability. This is due to the nonzero brane tension .

Using the series expansion of the eigenvalue (truncated at 7th order), we numerically integrate Eq. (9) and compute the absorption probability. We defer to Appendix the analytic calculation of the absorption probability in the limit and . The analytic and numerical results are found to be in good agreement in the regime and .

As was remarked above, both and become larger as the brane tension increases. This results in the enhancement of the “potential” (36), which will reduce the absorption probability. A typical example of the absorption probability (as a function of ) is plotted for various and in Fig. 2, showing that the absorption probability indeed decreases with decreasing . We can see in Fig. 3 that there appear superradiant modes with , for which the absorption probability is negative (48).

Figure 4: Power emission spectra for different values of and .
Figure 5: Angular momentum spectra for different values of and .

From the absorption probability we compute the energy and angular momentum emission rates. They are given by the formulas


where and were already defined in Eqs. (6) and (39), respectively. We summed up to modes in calculating the above quantities. Our numerical results are summarized in Figs. 4 and 5.

We find that the finite brane tension reduces the power and angular momentum emission spectra. This is anticipated from the behavior of the absorption probability stated above. The power emission rates are reduced in the low frequency regime as the rotation parameter increases. This is because the Hawking temperature becomes lower with increasing . However, the power emission rates are enhanced at high frequencies. This should be caused by superradiance. The angular momentum emission rates are enhanced over the whole frequency range with increasing because the effects of superradiance win out. (In the previous estimates (44); (45); (49) significant superradiance is observed as the angular momentum increases.) One can see from Fig. 6 that a large portion of the contribution to the power emission spectrum is coming from the modes, which show the superradiant behavior (except for ). We confirmed that the same is true for the angular momentum spectrum. As increases the amplitudes of these modes are enhanced, lifting up the total spectra. One can also see that the oscillatory behavior of the spectra is due to these modes, with each peak corresponding to each mode. However, the above things do not mean that only the modes are important in determining the behavior of the spectra; the eigenvalues for the modes are independent of the conical deficit and the other modes () determine the decrease in the spectra with decreasing .

Figure 6: Contributions from modes. Red (solid) line shows the total power emission spectrum depicted also in Fig. 4, while green (solid) line represents the sum of modes, each of which is shown by blue (dashed) line.

Since the brane induced metric is the same as the tensionless case, the Hawking spectrum of brane-localized fields does not depend on the tension (for fixed ).3 The emission of brane scalar fields has been studied in Refs. (21); (22); (10); (11); (13) for rotating black holes in various dimensions, and the brane-to-bulk ratio of the energy emission rates has been discussed in (25), implying that the dominant channels are the brane-localized modes. Since the finite brane tension further suppresses the bulk field contributions, it is likely that brane-localized scalar emission dominates bulk emission also in the present case. Although bulk modes will not be observed directly, bulk emission is still important because it indirectly determines the amount of energy and angular momentum left for brane-localized emission.

V Concluding remarks

The black hole production in TeV scale gravity offers us a possible window to explore the presence of extra dimensions. In this paper we have studied Hawking emission of scalar fields into the bulk from a rotating black hole localized on a codimension-2 brane. The exact solution we used is a rare example in which we can treat self-gravity of the brane consistently in the context of brane-localized black holes (39). This simple model enables us to elucidate how the brane tension modifies the Hawking spectra relative to the tensionless case.

Assuming the separable ansatz for the scalar field, we have determined the angular eigenvalue in the series expansion form: . This analytic approximation (truncated at 7th order) was turned out to be in excellent agreement with the numerical result in the regime . Using the analytic form of , we then integrated the radial equation numerically and computed the power and angular momentum emission spectra. Our finding is that the finite brane tension suppresses the power and angular momentum spectra. We also showed that the power emission rates are reduced at low frequencies and enhanced at high frequencies as the rotation parameter increases. The suppression at low frequencies is due to the smaller Hawking temperature and the enhancement at high frequencies is caused by superradiance. The angular momentum emission rates are enhanced over the whole frequency range as increases due to superradiance. The spectra and the amounts of energy and angular momentum radiated away into the bulk are thus determined by the interplay of these effects. To conclude, the brane tension plays an important role in the evaporation process in the life of a mini black hole.

It is possible to obtain an analytic but approximate solution to the radial equation, as has been done recently in (25) and is replicated in Appendix. The range of validity of this analytic approximation is restricted to and . Therefore, even in the tensionless case () our result is new, in that we have extended the range of validity of (25) by invoking the numerical approach. We however exploited the analytic expression for to simplify numerical calculations. The eigenvalues obtained here by the continued fraction method are quite accurate even for , which allows us to compute the spectra for and . In order to explore the full regime extending to , we need to solve both the radial and angular equations numerically, which is left to further investigation.

In this paper we considered only the emission of scalar fields. It would be interesting to study the bulk emission of higher spin fields and determine the brane-to-bulk ratio of the energy and angular momentum emission rates. Recently, it has been reported that for fermion fields the bulk emission dominates the brane-localized emission in a six- or higher dimensional Schwarzschild background (27). Therefore, investigating the effects of the finite brane tension (43) and the black hole rotation on the fermion emission would be of particular interest. Another open issue is to clarify the spin-down evolution of the rotating black hole and the effect of brane tension on it. This process can be studied along the line of (44); (45). We plan to return to this issue in the near future. Finally, it is fair to say that all of the results in this paper have been derived ignoring the compactification mechanism, which may affect both brane and bulk emissions. This point is worth exploring, though constructing brane-localized black hole solutions in compact space will be quite difficult.

We would like to thank Sam Dolan and Marc Casals for useful comments. TK, MN, and YT are supported by the JSPS under Contract Nos. 19-4199, 19-204, and 17-53192.

Appendix A Analytic approximation method for and

In this appendix we present an approximation method to obtain an analytic expression for the absorption probability. Our result here simply generalizes that of (25) to include the deficit angle (see also (13)). The procedure is as follows: first we obtain the asymptotic solutions in the far-field and near-horizon regions, and then match the two solutions in the intermediate zone. The approximation is valid in the low-energy () and slow-rotation () regime. We check that in this range the analytic result agrees with the numerical one displayed in the main text.

Let us first focus on the near-horizon zone (). It is convenient to work with a new radial variable defined by


for which we have


with . Now the horizon is located at , while the asymptotic infinity corresponds to . Near the horizon, the radial equation (9) reduces to



Performing the transformation with


we obtain a hypergeometric differential equation


where , , and . Thus, the near-horizon solution is given in terms of the hypergeometric function by


where are integration constants. In the limit () we have


where is defined earlier in the main text. Since we are imposing the boundary condition such that no outgoing wave is present at the horizon, we set . Thus we arrive at


One can check that the convergence condition for the hypergeometric function, Re, is indeed satisfied.

Figure 7: Absorption probabilities computed analytically and numerically. Solid lines follow from the analytic approximation, while dashed lines denote our numerical result.

Now we extend the solution (52) to go beyond the near-horizon zone. Using the formula

and taking the limit () we obtain



The expression (53) should be matched to the far-field solution which will be derived below.

Let us go on to the far-field solution. In the far-field zone () the radial equation (9) reduces to


The solution is given by


where () is the Bessel function of the first (second) kind, , and are the integration constants. This solution is in turn to be extended to the near-horizon zone. Taking the limit we get


The two solutions (53) and (57) appear to have different powers in , but by taking the low-energy () and slow-rotation () limit we are able to match them. Neglecting terms of order , , and , we have , , and , leading to and in Eq. (53), where we used . As for Eq. (57), we have and . Thus in this limit we achieve exact matching.

The absorption probability can be expressed in terms of the ratio of the coefficients . One finds


In the limit the solution (56) can be written as




The absorption probability is given by . In the low energy limit () we have . Hence,


Substituting (A) to (62) results in a lengthy expression, but when is not large we may expand the equation with respect to in the low-energy and slow-rotation regime. Then, at leading order in , one finds


In Fig. 7 we compare this analytic result with our numerical calculation employed in the main text. It can be seen that the numerical result indeed agrees with the analytic one in the low-energy and slow-rotation regime. Since all the terms except on the right hand side of Eq. (63) are positive, the sign of is controlled by . Therefore, the absorption probability is negative for superradiant modes ().

In the low energy regime the mode will be dominant, for which the absorption probability is given by


This does not depend on the deficit angle.

To obtain the expression for the absorption cross section we extract the ingoing s-wave from the plane wave (50); (51):


where the square root in the denominator is the normalization factor. In order to determine , one integrates both sides of Eq. (65) over the 4-sphere with a deficit angle, and then, looking at the far region , extracts only the ingoing modes. Thus, we arrive at


leading to the low energy absorption cross section


In the static limit we have (horizon area), which reproduces the general result for the spherically symmetric case (50).


  1. preprint: WU-AP/277/07
  2. The second order coefficient in (47) contains a sign error, which has been corrected in (25).
  3. Although this statement is true in the present model, it is not clear whether or not the induced metric is independent of the brane tension in more realistic situations with compact extra dimensions.


  1. N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998) [arXiv:hep-ph/9803315]; N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Rev. D 59, 086004 (1999) [arXiv:hep-ph/9807344]; I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436, 257 (1998) [arXiv:hep-ph/9804398].
  2. C. Kokorelis, Nucl. Phys. B 677, 115 (2004) [arXiv:hep-th/0207234].
  3. L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) [arXiv:hep-ph/9905221]; L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999) [arXiv:hep-th/9906064].
  4. S. B. Giddings and S. D. Thomas, Phys. Rev. D 65, 056010 (2002) [arXiv:hep-ph/0106219]; S. Dimopoulos and G. L. Landsberg, Phys. Rev. Lett. 87, 161602 (2001) [arXiv:hep-ph/0106295]; S. Hossenfelder, S. Hofmann, M. Bleicher and H. Stoecker, Phys. Rev. D 66, 101502 (2002) [arXiv:hep-ph/0109085].
  5. D. M. Eardley and S. B. Giddings, Phys. Rev. D 66, 044011 (2002) [arXiv:gr-qc/0201034]; H. Yoshino and Y. Nambu, Phys. Rev. D 66, 065004 (2002) [arXiv:gr-qc/0204060]; H. Yoshino and Y. Nambu, Phys. Rev. D 67, 024009 (2003) [arXiv:gr-qc/0209003]; E. Kohlprath and G. Veneziano, JHEP 0206, 057 (2002) [arXiv:gr-qc/0203093]; E. Berti, M. Cavaglia and L. Gualtieri, Phys. Rev. D 69, 124011 (2004) [arXiv:hep-th/0309203]; N. Kaloper and J. Terning, Gen. Rel. Grav. 39, 1525 (2007) [arXiv:0705.0408 [hep-th]].
  6. S. W. Hawking, Commun. Math. Phys. 43, 199 (1975) [Erratum-ibid. 46, 206 (1976)].
  7. R. Emparan, G. T. Horowitz and R. C. Myers, Phys. Rev. Lett. 85, 499 (2000) [arXiv:hep-th/0003118].
  8. P. Kanti and J. March-Russell, Phys. Rev. D 66, 024023 (2002) [arXiv:hep-ph/0203223]; P. Kanti and J. March-Russell, Phys. Rev. D 67, 104019 (2003) [arXiv:hep-ph/0212199].
  9. C. M. Harris and P. Kanti, JHEP 0310, 014 (2003) [arXiv:hep-ph/0309054]; P. Kanti, J. Grain and A. Barrau, Phys. Rev. D 71, 104002 (2005) [arXiv:hep-th/0501148]; S. Creek, O. Efthimiou, P. Kanti and K. Tamvakis, Phys. Lett. B 635, 39 (2006) [arXiv:hep-th/0601126].
  10. C. M. Harris and P. Kanti, Phys. Lett. B 633, 106 (2006) [arXiv:hep-th/0503010].
  11. G. Duffy, C. Harris, P. Kanti and E. Winstanley, JHEP 0509, 049 (2005) [arXiv:hep-th/0507274].
  12. M. Casals, P. Kanti and E. Winstanley, JHEP 0602, 051 (2006) [arXiv:hep-th/0511163]; M. Casals, S. R. Dolan, P. Kanti and E. Winstanley, JHEP 0703, 019 (2007) [arXiv:hep-th/0608193].
  13. S. Creek, O. Efthimiou, P. Kanti and K. Tamvakis, Phys. Rev. D 75, 084043 (2007) [arXiv:hep-th/0701288]
  14. S. Creek, O. Efthimiou, P. Kanti and K. Tamvakis, arXiv:0707.1768 [hep-th].
  15. E. Berti, K. D. Kokkotas and E. Papantonopoulos, Phys. Rev. D 68, 064020 (2003) [arXiv:gr-qc/0306106].
  16. D. Ida, Y. Uchida and Y. Morisawa, Phys. Rev. D 67, 084019 (2003) [arXiv:gr-qc/0212035].
  17. Y. Morisawa and D. Ida, Phys. Rev. D 71, 044022 (2005) [arXiv:gr-qc/0412070].
  18. V. Cardoso, G. Siopsis and S. Yoshida, Phys. Rev. D 71, 024019 (2005) [arXiv:hep-th/0412138].
  19. V. Cardoso, M. Cavaglia and L. Gualtieri, Phys. Rev. Lett. 96, 071301 (2006) [Erratum-ibid. 96, 219902 (2006)] [arXiv:hep-th/0512002]; V. Cardoso, M. Cavaglia and L. Gualtieri, JHEP 0602, 021 (2006) [arXiv:hep-th/0512116].
  20. P. Kanti and R. A. Konoplya, Phys. Rev. D 73, 044002 (2006) [arXiv:hep-th/0512257]; P. Kanti, R. A. Konoplya and A. Zhidenko, Phys. Rev. D 74, 064008 (2006) [arXiv:gr-qc/0607048].
  21. D. Ida, K. y. Oda and S. C. Park, Phys. Rev. D 67, 064025 (2003) [Erratum-ibid. D 69, 049901 (2004)] [arXiv:hep-th/0212108].
  22. D. Ida, K. y. Oda and S. C. Park, Phys. Rev. D 71, 124039 (2005) [arXiv:hep-th/0503052].
  23. D. Ida, K. y. Oda and S. C. Park, Phys. Rev. D 73, 124022 (2006) [arXiv:hep-th/0602188].
  24. V. P. Frolov and D. Stojkovic, Phys. Rev. D 67, 084004 (2003) [arXiv:gr-qc/0211055].
  25. S. Creek, O. Efthimiou, P. Kanti and K. Tamvakis, arXiv:0709.0241 [hep-th].
  26. A. S. Cornell, W. Naylor and M. Sasaki, JHEP 0602, 012 (2006) [arXiv:hep-th/0510009].
  27. H. T. Cho, A. S. Cornell, J. Doukas and W. Naylor, arXiv:0709.1661 [hep-th].
  28. P. Kanti, Int. J. Mod. Phys. A 19, 4899 (2004) [arXiv:hep-ph/0402168].
  29. V. P. Frolov and D. Stojkovic, Phys. Rev. D 66, 084002 (2002) [arXiv:hep-th/0206046]; V. P. Frolov, M. Snajdr and D. Stojkovic, Phys. Rev. D 68, 044002 (2003) [arXiv:gr-qc/0304083]; V. P. Frolov, D. V. Fursaev and D. Stojkovic, JHEP 0406, 057 (2004) [arXiv:gr-qc/0403002]; V. P. Frolov, D. V. Fursaev and D. Stojkovic, Class. Quant. Grav. 21, 3483 (2004) [arXiv:gr-qc/0403054].
  30. D. M. Gingrich, JHEP 0711, 064 (2007) [arXiv:0706.0623 [hep-ph]].
  31. D. C. Dai, G. Starkman, D. Stojkovic, C. Issever, E. Rizvi and J. Tseng, arXiv:0711.3012 [hep-ph].
  32. A. Flachi, O. Pujolas, M. Sasaki and T. Tanaka, Phys. Rev. D 74, 045013 (2006) [arXiv:hep-th/0604139]; A. Flachi and T. Tanaka, Phys. Rev. Lett. 95, 161302 (2005) [arXiv:hep-th/0506145]; A. Flachi and T. Tanaka, Phys. Rev. D 76, 025007 (2007) [arXiv:hep-th/0703019].
  33. R. Emparan, G. T. Horowitz and R. C. Myers, JHEP 0001, 007 (2000) [arXiv:hep-th/9911043].
  34. H. Kudoh, T. Tanaka and T. Nakamura, Phys. Rev. D 68, 024035 (2003) [arXiv:gr-qc/0301089].
  35. S. Creek, R. Gregory, P. Kanti and B. Mistry, Class. Quant. Grav. 23, 6633 (2006) [arXiv:hep-th/0606006].
  36. N. Tanahashi and T. Tanaka, arXiv:0712.3799 [gr-qc].
  37. F. R. Tangherlini, Nuovo Cim. 27, 636 (1963); R. C. Myers and M. J. Perry, Annals Phys. 172, 304 (1986).
  38. N. Kaloper and D. Kiley, JHEP 0603, 077 (2006) [arXiv:hep-th/0601110].
  39. D. Kiley, arXiv:0708.1016 [hep-th].
  40. D. C. Dai, N. Kaloper, G. D. Starkman and D. Stojkovic, Phys. Rev. D 75, 024043 (2007) [arXiv:hep-th/0611184].
  41. S. Chen, B. Wang and R. K. Su, Phys. Lett. B 647, 282 (2007) [arXiv:hep-th/0701209].
  42. U. A. al-Binni and G. Siopsis, arXiv:0708.3363 [hep-th].
  43. H. T. Cho, A. S. Cornell, J. Doukas and W. Naylor, arXiv:0710.5267 [hep-th].
  44. C. M. Chambers, W. A. Hiscock and B. Taylor, Phys. Rev. Lett. 78, 3249 (1997) [arXiv:gr-qc/9703018]; B. E. Taylor, C. M. Chambers and W. A. Hiscock, Phys. Rev. D 58, 044012 (1998) [arXiv:gr-qc/9801044].
  45. H. Nomura, S. Yoshida, M. Tanabe and K. i. Maeda, Prog. Theor. Phys. 114, 707 (2005) [arXiv:hep-th/0502179].
  46. E. W. Leaver, Proc. Roy. Soc. Lond. A 402 (1985) 285.
  47. E. Berti, V. Cardoso and M. Casals, Phys. Rev. D 73, 024013 (2006) [Erratum-ibid. D 73, 109902 (2006)] [arXiv:gr-qc/0511111].
  48. A. A. Starobinsky, JETP 37, 28 (1973); A. A. Starobinsky and S.M. Churilov, JETP 38, 1 (1973).
  49. M. Nozawa and K. i. Maeda, Phys. Rev. D 71, 084028 (2005) [arXiv:hep-th/0502166].
  50. S. R. Das, G. W. Gibbons and S. D. Mathur, Phys. Rev. Lett. 78, 417 (1997).
  51. S. R. Das and S. D. Mathur, Nucl. Phys. B 478, 561 (1996) [arXiv:hep-th/9606185].