<<

Global spectral representations of black hole in the complex

Maurice H. P. M. van Putten†

Massachusetts Institute of Technology–Laser Interferometer Gravitational-Wave Observatory, 175 Albary Street, NW17-161, Cambridge, MA 02139

Edited by Robert M. Wald, University of Chicago, Chicago, IL, and approved November 17, 2005 (received for review September 26, 2005) Binary black hole coalescence produces a finite burst of gravita- which shows that G is analytic in each coordinate away from the tional radiation, which propagates toward quiescent infinity. singularity at y ϭ p with domain of convergence r Ͼ p.Wenow These spacetimes are analytic about infinity and contain a dimen- complexify the real nonnegative radial coordinate r Ն 0 into a sionless coupling constant M͞s, where M denotes the total mass- complex variable z by energy and s is an imaginary . This work introduces a globally convergent Fourier series in a complex radial coordinate, z ϭ x ϩ is ͑s Ͼ p͒, [3] allowing spectral representations of black hole spacetimes in all three . This method circumvents singularities associ- where the real variable x assumes both positive and negative 1 ated with black holes and includes infinity, thereby avoiding the values. The analytic continuation of Eq. given by need to impose radiative boundary conditions. We illustrate this p lϩ1 ϭ Ϫ1 ͸ͩ ͪ ͑ ␪ representation theory on Boyer–Lindquist initial data of two black G p Pl cos ), [4] holes and a scalar wave with signal recovery by Cauchy’s z integral formula. is globally convergent on the straight Ϫϱ Ͻ x Ͻϱoutside ͉z͉ ϭ p and to the real line s ϭ 0. (The corresponding ͉ computation gravitation branch cut for the root is (Ϫϱ, 0].) To be precise, for each ␪ Eq. 4 defines a Taylor series in 1͞z, and, for each z, it defines apid progress in sensitivity in the broadband gravitational- a convergent expansion in Legendre functions. Thus, p͞s is a ASTRONOMY Rwave observatories (1–4) is creating new opportunities for running coupling constant, and it must be less than one for searches for the bursts of gravitational-wave emissions produced convergence. We shall refer to M͞s as the fixed coupling by coalescing black holes. Detection strategies will benefit constant, noting that p͞s is decreasing in binary coalescences. greatly from a priori understanding of their gravitational-wave Analyticity of the in Im(z) Ն ͉p͉ defines an optimal emissions, which includes matched filtering techniques against a distribution of points on z ϭ x ϩ is for numerical implementation catalog of signals precomputed by numerical relativity. It poses by the conformal transformation an interesting challenge of designing highly efficient and stable computational algorithms. M w ϭ . [5] Numerical relativity is particularly challenging in calculating z bursts and preceding chirps over many wave-periods in the Ϫ1 Ϫ1 presence of singularities associated with black hole spacetimes. Here, w describes a of (2s) about w0 ϭ (2is) , Spectral methods provide an attractive approach. They are whose interior corresponds to Im(z) Ն s. There exists a Taylor optimal in efficiency and accuracy, both in amplitude and phase, series in w about w0, which is convergent up to and including the Ϫ ␺ provided that the metric is smooth everywhere. Here, we focus boundary C: w(␺) ϭ (2is) 1(ei ϩ 1). The corresponding Fourier on bursts of radiation produced by black hole coalescence, which series in ␺ can be efficiently computed by using the fast Fourier Ϫmϩ1 propagate toward quiescent infinity. By causality, these - transform on a uniform distribution of points ␺k ϭ 2 ␲k(k ϭ times preserve quiescence at infinity for all finite time. 1,2,...,2m)onC. Their images on z ϭ x ϩ is produce a unique We can introduce globally smooth metrics for black hole nonuniform distribution of points spacetimes by using complex . To illustrate this ap- ϭ ͓ ␺ ͞ ͒ ϩ proach, we note that the asymptotic structure of black hole zk Ms tan( k 2 i], [6] spacetimes, which are asymptotically and quiescent at in- whose large-scale cut-off is by the proximity of the {␺k}to finity, shares the same asymptotic properties as the Green‘s Ϯ␲. function G of Minkowski . Consequently, these space- ͞ We apply these observations to Boyer–Lindquist initial data times carry a dimensionless coupling constant M s, where M for two black holes (5, 6), described by the potential denotes the total mass-energy of the space and s is a finite step

into the . The Green’s function of the Laplacian M1 M2 in flat spacetime, ⌽ ϭ 1 ϩ ϩ , ͱz2 ϩ p2 Ϫ 2pz cos ␪ ͱz2 ϩ q2 Ϫ 2qz cos ␪ 1 G͑y, p͒ ϭ , [1] [7] D͑y, p͒ for a pair of black holes with masses M1,2 at positions p and q is a function of the distance D between an observer at location along the polar axis of a spherical (r, ␪, ␾). y and a source at location p. We can expand G in Legendre The Boyer–Lindquist data represent an exact solution to polynomials Pl upon expressing y and p in spherical coordinates Laplace’s equation, representing the solution of the Hamilto- (r, ␪, ␾)

1 p lϩ1 Conflict of interest statement: No conflicts declared. ϭ ϭ Ϫ1 ͸ͩ ͪ ͑ ␪ G p P cos ), This paper was submitted directly (Track II) to the PNAS office. ͱr2 ϩ p2 Ϫ 2pr cos ␪ r l †E-mail: [email protected]. [2] © 2006 by The National Academy of Sciences of the USA

www.pnas.org͞cgi͞doi͞10.1073͞pnas.0508406103 PNAS ͉ January 17, 2006 ͉ vol. 103 ͉ no. 3 ͉ 519–522 Downloaded by guest on September 28, 2021 Fig. 1. The Boyer–Lindquist scalar field ⌽(z, ␪) possesses a Fourier–Legendre expansion in w ϭ M͞z, z ϭ x ϩ is (Upper) and cos␪ (Lower), here shown for a symmetric configuration consisting of two black holes of mass M1,2 ϭ 1 positioned at the north and south pole p ϭϮ1 of a spherical coordinate system (r, ␪, ␾). The Fourier coefficients show exponential decay by analyticity in Im(z) Ն s Ͼ s*(␪) ϭ psin␪ shown with cos␪ ϭ 0.25 and s* ϭ 0.9682. It gives rise to spectral accuracy in the first and second coordinate derivatives (Right Upper). The of the Legendre coefficients shows exponential decay when s is larger than p (Left). It gives rise to spectral accuracy in the first and second derivatives with respect to ␪ (Right Lower).

nian constraint for time-symmetric initial data. For this rea- binary black hole spacetime with total mass-energy M as mea- son, the Boyer–Lindquist data are representative in general of sured at infinity, defined by the residue the asymptotic properties of the Green’s function in curved ϱϩis spacetime. i 1/4 Ϯ ␪ ϭ E ͵ ͑ Ϫ ͒ [10] ϭ i M R ͭ hzz 1 dzͮ . It will be appreciated that the singularities z pe are ␲ essential in making ⌽ nontrivial: the alternative of an entire Ϫϱϩis function that is analytic at infinity would be constant by Li- For binary black hole coalescence, M is the sum M ϩ M of the ouville’s . For complex z ϭ x ϩ is with s larger than both 1 2 masses of the two black holes, plus the energy in gravitational p and q, ⌽(z, ␪) is finite and pointwise analytic everywhere on Ϫϱ ϩ Ͻ Ͻϱϩ radiation and minus the binding energy in the system. Thus, M is z is. Following Eqs. 4 and 6, we capture is a constant. In contrast, M is a nonincreasing function of time ϭ ͞ analyticity by the Fourier transform on the image w M z.It on any finite computational domain with outgoing radiation gives rise to an efficient and accurate functional representation boundary conditions. For large z, the line-element (Eq. 9) ␪ in (z, ) with spectral accuracy in derivatives by spectral differ- satisfies entiation. Fig. 1 shows the computational results on the alter- native representations M ds2 ϭ ϪN2dt2 ϩ ͩ 1 ϩ ͪ ͑dz2 ϩ d⌺2͒ ϩ ..., [11] 2 z ⌽͑ ␪͒ ϭ ͸N ͑ ␪ n ϭ ͸ L ͑ ͒ ͑ ␪ z, nϭ1an cos )w lϭ0bl z Pl cos ), [8] up to order (M͞z)2, where d⌺2 ϭ z2d␪2 ϩ z2 sin2␪d␾2 denotes the where the series are finite in our numerical implementations. line-element of the coordinate of radius z. Accordingly, The application to black hole spacetimes in the complex plane the three-metric hij is decomposed into diagonal and off- can be considered in the line-element diagonal elements given by the matrix factorization 1 Ϫ M͞2 z 4 2 ϭ Ϫ 2 2 ϩ i j ͑ ͒ ϭ A1 B1 B2 ds N dt hij dx dx , N z ϩ ͞ , [9] 4 1 M 2 z h ϭ ␥ ͩB1 A2 B3ͪ␥. [12] ij 4 B2 B3 A3 for a three-metric hij. Working in the complex plane, we are at liberty to choose a lapse function N and vanishing shift functions In accord with the spectral representation of Eq. 4, the potential in accord with the asymptotic Schwarzschild structure of space- Eqs. 7 and 8, the coefficients (Ai, Bi) satisfy the general expan- time at large distances. It applies without loss of generality to any sion in spherical harmonics

520 ͉ www.pnas.org͞cgi͞doi͞10.1073͞pnas.0508406103 van Putten Downloaded by guest on September 28, 2021 Fig. 3. Signal recovery of a wave propagating to quiescent infinity is shown Fig. 2. Normalized errors in the first () and second () derivatives for a quantity u(t, z)onz ϭ x ϩ is, here satisfying the scalar wave equation utt ϭ Ϫ2 2 Ϫ2x2 of the metric and the Ricci scalar tensor (asterisks) based on a Fourier– r (r ur)r with time-symmetric initial data at t ϭ 0 u(0, x) ϭ e cos(10x) and Legendre representation of Boyer–Lindquist initial data of Fig. 1 with M͞s ϭ ut(0, x) ϭ 0. The simulations use the Fourier representation on the complex 4 1͞2. The Fourier–Legendre results are computed with N ϭ 32 Fourier coeffi- radius under w ϭ M͞z and the wave-equation utt ϭ w uww using N ϭ 1,024 cients on w ϭ M͞z varying degrees L ϭ 4,...,32ofLegendre polynomials. The points. Plotted is the quantity ru on r Ն 0 for t0 ϭ 0 and t1 ϭ 4. The numerical results are compared with a mixed finite differencing-Legendre method with results (dots) accurately track the exact results (continuous line, by the method varying degree N ϭ 2m, m ϭ 5,..., 11ofthe number of points in the of characteristics). The errors are essentially independent of s, i.e., 4.5 ϫ 10Ϫ4 z-coordinate and fixed expansion in L ϭ 16 Legendre polynomials. This com- for s ϭ 0.5 and 4.7 ϫ 10Ϫ4 for s ϭ 1.0, as implied by analyticity in s Ͼ 0. ASTRONOMY parison highlights the accuracy and efficiency of the Fourier-spectral method, even at low resolutions. The results are practically unchanged for other configurations of Boyer–Lindquist data, and the errors decrease for larger s. ͑ ͒ ϭ 4 ͑ ϭ ͞ ͒ Ј utt t, z w uww w 1 z K : ͭ ͑ ͒ ϭ ͑ ͒ ͑ ͒ ϭ . [18] u 0, z u0 z , ut 0, z 0 ␥ ϭ ␥ ϭ ␥ ϭ ␪ 11 1, 22 z, 33 zsin . [13] [In general, the analytic continuation of the initial data are the extension of data on (r, ␪, ␾) and (r, ␲ Ϫ ␪, ␾ ϩ ␲).] We can ϭ ͸ ml͑ ͒ ͑␪ ␾͒ Ai M ͉ m ͉ Յ l ai t, w Yml , , [14] implement KЈ numerically using finite-differencing and leapfrog ϭ ͸ ml͑ ͒ ͑␪ ␾͒ time-stepping. Fig. 3 shows a representative numerical result, ͉ ͉ Յ [15] Bi M m l bi t, w Yml , . which is verified to satisfy independence of s. The scalar wave- Together with a Fourier expansion in the 2␲-periodic variable ␺ equation with time-symmetric initial data shows the appearance ␺ of both left and right movers, the first propagating toward z ϭ of w( ), Eqs. 14 and 15 define a spectral representation in all Ϫϱ ϩ ϭϱϩ three coordinates. Fig. 2 shows the results in case of the is and the second propagating toward z is. Only the axisymmetric Boyer–Lindquist initial data, including the com- right movers are projected onto the observer’s nonnegative radii z Ն 0. The left movers are ultimately dissipated ‘‘unseen’’ in their putational error behavior in the scalar Ricci tensor. Here, we use Ϫϱ ϩ the general expansions (Eq. 8) to probe independently the propagation to is, while maintaining a contribution to the ␺ total mass-energy M in Eq. 10. We emphasize that the projection dependencies on the degree N of the Fourier series of w( ) and ͞ the maximal degree L of the spherical harmonics. onto real radii is only required at the observer, that is, LIGO Time-evolution of h on z ϭ x ϩ is gives rise to burst of Virgo, and hence for large distances away from the source. ij (Projections are not needed to recover a real representation of radiation that propagates toward quiescent infinity, where the ͞ observer‘s signal is defined on nonnegative radii z Ն 0. This the source itself for the purpose of LIGO Virgo searches of signal can be calculated by projection according to Cauchy’s gravitational radiation.) integral formula To conclude, the complex plane offers a unique opportunity for circumventing spacetime singularities and enabling spectral ϱϩ representations in all three dimensions, consisting of spherical i is u͑z, ␪, ␾, t͒ u ͑r, ␪, ␾, t͒ ϭ REͭ ͵ dzͮ , [16] harmonics in the angular coordinates and Fourier series in the ␲ z Ϫ r Ϫϱϩis complex radial coordinate. This approach introduces a compu- tational duality between weakly nonlinear spacetimes in the where u denotes a relevant metric component or tensor quantity. complex plane and strongly nonlinear spacetimes on the real Time-evolution on z ϭ x ϩ is is stable with constant translation line, by virtue of Cauchy’s integral formula. It creates a most s, because it preserves reality of the dispersion relation. We can efficient representation of black hole spacetimes and promises a illustrate this on the linear wave-equation for a scalar u(t, reduction in computational effort by orders of magnitude com- r) in spherical coordinates (t, r, ␪, ␾) with time-symmetric initial pared with standard finite differencing approaches. It will be of data, given by interest to take advantage of this approach in time-dependent calculations described by the fully nonlinear for gen- ͑ ͒ ϭ Ϫ2͓ 2 ͑ ͔͒ utt t, r r r ur t, r r eral relativity in any of its hyperbolic forms (e.g., ref. 7). K: ͭ ͑ ͒ ϭ ͑ ͒ ͑ ͒ ϭ . [17] u 0, r u0 r , ut 0, r 0 I thank G. ’t Hooft for stimulating discussions and the referee for a By analytic continuation, the equivalent initial value problem KЈ detailed consideration of the manuscript. Part of this work was supported with complex radial coordinate on Ϫϱ ϩ is Ͻ z Ͻϱϩis is by the Laser Interferometer Gravitational-Wave Observatories, con-

van Putten PNAS ͉ January 17, 2006 ͉ vol. 103 ͉ no. 3 ͉ 521 Downloaded by guest on September 28, 2021 structed by Caltech and the Massachusetts Institute of Technology with PHY 9210038. The Laser Interferometer Gravitational-Wave Observa- funding from National Science Foundation Cooperative Agreement tory Laboratory operates under cooperative agreement PHY-0107417.

1. Abramovici, A., Althouse, W. E., Drever, R. W. P., Gursel, Y., Kawamura, S., 4. Bradaschia, C., del Fabbro, R., di Virgilio, A., Flaminio, R., Giazotto, A., Raab, F. J., Shoemaker, D., Sievers, L., Spero, R. E. & Thorne, K. S. (1992) Michelozzi, B., Montelatici, V., Passuello, D., Brillet, A., Man, N. C., et al. Science 292, 325–333. (1992) Phys. Lett. A 163, 15–20. 2. Barish, B. & Weiss, R. (1999) Phys. Today 52, 44–50. 5. Cook, G. B. (1991) Phys. Rev. D Part. Fields 44, 2983–3000. 3. Acernese, F., Amico, P., Arnaud, N., Arnault, C., Babusci, D., Ballardin, G., 6. Price, R. H. & Pullin, J. (1994) Phys. Rev. Lett. 72, 3297–3300. Barone, F., Barsuglia, M., Bellachia, F., Beney, J. L., et al. (2002) Classical 7. van Putten, M. H. P. M. & Eardley, D. M. (1996) Phys. Rev. D Part. Fields 53, Quant. Gravity 19, 1421–1428. 3056–3063.

522 ͉ www.pnas.org͞cgi͞doi͞10.1073͞pnas.0508406103 van Putten Downloaded by guest on September 28, 2021