-
Gauge-invariant perturbations of Schwarzschild spacetime
Authors:
Abhay G. Shah,
Bernard F. Whiting,
Steffen Aksteiner,
Lars Andersson,
Thomas Bäckdahl
Abstract:
We study perturbations of Schwarzschild spacetime in a coordinate-free, covariant form. The GHP formulation, having the advantage of not only being covariant but also tetrad-rotation invariant, is used to write down the previously known odd- and even-parity gauge-invariants and the equations they satisfy. Additionally, in the even-parity sector, a new invariant and the second order hyperbolic equa…
▽ More
We study perturbations of Schwarzschild spacetime in a coordinate-free, covariant form. The GHP formulation, having the advantage of not only being covariant but also tetrad-rotation invariant, is used to write down the previously known odd- and even-parity gauge-invariants and the equations they satisfy. Additionally, in the even-parity sector, a new invariant and the second order hyperbolic equation it satisfies are presented. Chandrasekhar's work on transformations of solutions for perturbation equations on Schwarzschild spacetime is translated into the GHP form, i.e., solutions for the equations of the even- and odd-parity invariants are written in terms of one another, and the extreme Weyl scalars; and solutions for the equations of these latter invariants are also written in terms of one another. Recently, further gauge invariants previously used by Steven Detweiler have been described. His method is translated into GHP form and his basic invariants are presented here. We also show how these invariants can be written in terms of curvature invariants.
△ Less
Submitted 24 November, 2016;
originally announced November 2016.
-
iResum: a new paradigm for resumming gravitational wave amplitudes
Authors:
Alessandro Nagar,
Abhay G. Shah
Abstract:
We introduce a new, resummed, analytical form of the post-Newtonian (PN), factorized, multipolar amplitude corrections $f_{\ell m}$ of the effective-one-body (EOB) gravitational waveform of spinning, nonprecessing, circularized, coalescing black hole binaries (BBHs). This stems from the following two-step paradigm: (i) the factorization of the orbital (spin-independent) terms in $f_{\ell m}$; (ii)…
▽ More
We introduce a new, resummed, analytical form of the post-Newtonian (PN), factorized, multipolar amplitude corrections $f_{\ell m}$ of the effective-one-body (EOB) gravitational waveform of spinning, nonprecessing, circularized, coalescing black hole binaries (BBHs). This stems from the following two-step paradigm: (i) the factorization of the orbital (spin-independent) terms in $f_{\ell m}$; (ii) the resummation of the residual spin (or orbital) factors. We find that resumming the residual spin factor by taking its inverse resummed (iResum) is an efficient way to obtain amplitudes that are more accurate in the strong-field, fast-velocity regime. The performance of the method is illustrated on the $\ell=2$ and $m=(1,2)$ waveform multipoles, both for a test-mass orbiting around a Kerr black hole and for comparable-mass BBHs. In the first case, the iResum $f_{\ell m}$'s are much closer to the corresponding "exact" functions (obtained solving numerically the Teukolsky equation) up to the light-ring, than the nonresummed ones, especially when the black-hole spin is nearly extremal. The iResum paradigm is also more efficient than including higher post-Newtonian terms (up to 20PN order): the resummed 5PN information yields per se a rather good numerical/analytical agreement at the last-stable-orbit, and a well-controlled behavior up to the light-ring. For comparable mass binaries (including the highest PN-order information available, 3.5PN), comparing EOB with Numerical Relativity (NR) data shows that the analytical/numerical fractional disagreement at merger, without NR-calibration of the EOB waveform, is generically reduced by iResum, from a $40\%$ of the usual approach to just a few percents. This suggests that EOBNR waveform models for coalescing BBHs may be improved using iResum amplitudes.
△ Less
Submitted 28 September, 2016; v1 submitted 1 June, 2016;
originally announced June 2016.
-
Self-force as a cosmic censor in the Kerr overspinning problem
Authors:
Marta Colleoni,
Leor Barack,
Abhay G. Shah,
Maarten van de Meent
Abstract:
It is known that a near-extremal Kerr black hole can be spun up beyond its extremal limit by capturing a test particle. Here we show that overspinning is always averted once back-reaction from the particle's own gravity is properly taken into account. We focus on nonspinning, uncharged, massive particles thrown in along the equatorial plane, and work in the first-order self-force approximation (i.…
▽ More
It is known that a near-extremal Kerr black hole can be spun up beyond its extremal limit by capturing a test particle. Here we show that overspinning is always averted once back-reaction from the particle's own gravity is properly taken into account. We focus on nonspinning, uncharged, massive particles thrown in along the equatorial plane, and work in the first-order self-force approximation (i.e., we include all relevant corrections to the particle's acceleration through linear order in the ratio, assumed small, between the particle's energy and the black hole's mass). Our calculation is a numerical implementation of a recent analysis by two of us [Phys.\ Rev.\ D {\bf 91}, 104024 (2015)], in which a necessary and sufficient "censorship" condition was formulated for the capture scenario, involving certain self-force quantities calculated on the one-parameter family of unstable circular geodesics in the extremal limit. The self-force information accounts both for radiative losses and for the finite-mass correction to the critical value of the impact parameter. Here we obtain the required self-force data, and present strong evidence to suggest that captured particles never drive the black hole beyond its extremal limit. We show, however, that, within our first-order self-force approximation, it is possible to reach the extremal limit with a suitable choice of initial orbital parameters. To rule out such a possibility would require (currently unavailable) information about higher-order self-force corrections.
△ Less
Submitted 17 August, 2015;
originally announced August 2015.
-
Metric perturbations produced by eccentric equatorial orbits around a Kerr black hole
Authors:
Maarten van de Meent,
Abhay G. Shah
Abstract:
We present the first numerical calculation of the (local) metric perturbation produced by a small compact object moving on an eccentric equatorial geodesic around a Kerr black hole, accurate to first order in the mass ratio. The procedure starts by first solving the Teukolsky equation to obtain the Weyl scalar $ψ_4$ using semi-analytical methods. The metric perturbation is then reconstructed from…
▽ More
We present the first numerical calculation of the (local) metric perturbation produced by a small compact object moving on an eccentric equatorial geodesic around a Kerr black hole, accurate to first order in the mass ratio. The procedure starts by first solving the Teukolsky equation to obtain the Weyl scalar $ψ_4$ using semi-analytical methods. The metric perturbation is then reconstructed from $ψ_4$ in an (outgoing) radiation gauge, adding the appropriate non-radiative contributions arising from the shifts in mass and angular momentum of the spacetime.
As a demonstration we calculate the generalized redshift $U$ as a function of the orbital frequencies $Ω_r$ and $Ω_φ$ to linear order in the mass ratio, a gauge invariant measure of the conservative corrections to the orbit due to self-interactions. In Schwarzschild, the results surpass the existing result in the literature in accuracy, and we find new estimates for some of the unknown 4PN and 5PN terms in the post-Newtonian expansion of $U$. In Kerr, we provide completely novel values of $U$ for eccentric equatorial orbits. Calculation of the full self-force will appear in a forthcoming paper.
△ Less
Submitted 25 May, 2020; v1 submitted 15 June, 2015;
originally announced June 2015.
-
Experimental mathematics meets gravitational self-force
Authors:
Nathan K. Johnson-McDaniel,
Abhay G. Shah,
Bernard F. Whiting
Abstract:
It is now possible to compute linear in mass-ratio terms in the post-Newtonian (PN) expansion for compact binaries to very high orders using black hole perturbation theory applied to various invariants. For instance, a computation of the redshift invariant of a point particle in a circular orbit about a black hole in linear perturbation theory gives the linear-in-mass-ratio portion of the binding…
▽ More
It is now possible to compute linear in mass-ratio terms in the post-Newtonian (PN) expansion for compact binaries to very high orders using black hole perturbation theory applied to various invariants. For instance, a computation of the redshift invariant of a point particle in a circular orbit about a black hole in linear perturbation theory gives the linear-in-mass-ratio portion of the binding energy of a circular binary with arbitrary mass ratio. This binding energy, in turn, encodes the system's conservative dynamics. We give a method for extracting the analytic forms of these PN coefficients from high-accuracy numerical data using experimental mathematics techniques, notably an integer relation algorithm. Such methods should be particularly important when the calculations progress to the considerably more difficult case of perturbations of the Kerr metric. As an example, we apply this method to the redshift invariant in Schwarzschild. Here we obtain analytic coefficients to 12.5PN, and higher-order terms in mixed analytic-numerical form to 21.5PN, including analytic forms for the complete 13.5PN coefficient, and all the logarithmic terms at 13PN. At these high orders, an individual coefficient can have over 30 terms, including a wide variety of transcendental numbers, when written out in full. We are still able to obtain analytic forms for such coefficients from the numerical data through a careful study of the structure of the expansion. The structure we find also allows us to predict certain "leading logarithm"-type contributions to all orders. The additional terms in the expansion we obtain improve the accuracy of the PN series for the redshift observable, even in the very strong-field regime inside the innermost stable circular orbit, particularly when combined with exponential resummation.
△ Less
Submitted 27 July, 2015; v1 submitted 9 March, 2015;
originally announced March 2015.
-
Raising and Lowering operators of spin-weighted spheroidal harmonics
Authors:
Abhay G. Shah,
Bernard F. Whiting
Abstract:
In this paper we generalize the spin-raising and lowering operators of spin-weighted spherical harmonics to linear-in-$γ$ spin-weighted spheroidal harmonics where $γ$ is an additional parameter present in the second order ordinary differential equation governing these harmonics. One can then generalize these operators to higher powers in $γ$. Constructing these operators required calculating the…
▽ More
In this paper we generalize the spin-raising and lowering operators of spin-weighted spherical harmonics to linear-in-$γ$ spin-weighted spheroidal harmonics where $γ$ is an additional parameter present in the second order ordinary differential equation governing these harmonics. One can then generalize these operators to higher powers in $γ$. Constructing these operators required calculating the $\ell$-, $s$- and $m$-raising and lowering operators (and various combinations of them) of spin-weighted spherical harmonics which have been calculated and shown explicitly in this paper.
△ Less
Submitted 3 June, 2016; v1 submitted 9 March, 2015;
originally announced March 2015.
-
Linear-in-mass-ratio contribution to spin precession and tidal invariants in Schwarzschild spacetime at very high post-Newtonian order
Authors:
Abhay G. Shah,
Adam Pound
Abstract:
Using black hole perturbation theory and arbitrary-precision computer algebra, we obtain the post-Newtonian (pN) expansions of the linear-in-mass-ratio corrections to the spin-precession angle and tidal invariants for a particle in circular orbit around a Schwarzschild black hole. We extract coefficients up to 20pN order from numerical results that are calculated with an accuracy greater than 1 pa…
▽ More
Using black hole perturbation theory and arbitrary-precision computer algebra, we obtain the post-Newtonian (pN) expansions of the linear-in-mass-ratio corrections to the spin-precession angle and tidal invariants for a particle in circular orbit around a Schwarzschild black hole. We extract coefficients up to 20pN order from numerical results that are calculated with an accuracy greater than 1 part in $10^{500}$. These results can be used to calibrate parameters in effective-one-body models of compact binaries, specifically the spin-orbit part of the effective Hamiltonian and the dynamically significant tidal part of the main radial potential of the effective metric. Our calculations are performed in a radiation gauge, which is known to be singular away from the particle. To overcome this irregularity, we define suitable Detweiler-Whiting singular and regular fields in this gauge, and we devise a rigorous mode-sum regularization method to compute the invariants constructed from the regular field.
△ Less
Submitted 2 April, 2015; v1 submitted 9 March, 2015;
originally announced March 2015.
-
Self-force from reconstructed metric perturbations: numerical implementation in Schwarzschild spacetime
Authors:
Cesar Merlin,
Abhay G. Shah
Abstract:
We present a first numerical implementation of a new scheme by Pound et al. that enables the calculation of the gravitational self-force in Kerr spacetime from a reconstructed metric-perturbation in a radiation gauge. The numerical task of the metric reconstruction essentially reduces to solving the fully separable Teukolsky equation, rather than having to tackle the linearized Einstein's equation…
▽ More
We present a first numerical implementation of a new scheme by Pound et al. that enables the calculation of the gravitational self-force in Kerr spacetime from a reconstructed metric-perturbation in a radiation gauge. The numerical task of the metric reconstruction essentially reduces to solving the fully separable Teukolsky equation, rather than having to tackle the linearized Einstein's equations themselves. The method offers significant computational saving compared to existing methods in the Lorenz gauge, and we expect it to become a main workhorse for precision self-force calculations in the future. Here we implement the method for circular orbits on a Schwarzschild background, in order to illustrate its efficacy and accuracy. We use two independent methods for solving the Teukolsky equation, one based on a direct numerical integration, and the other on the analytical approach of Mano, Suzuki, and Takasugi. The relative accuracy of the output self-force is at least $10^{-7}$ using the first method, and at least $10^{-9}$ using the second; the two methods agree to within the error bars of the first. We comment on the relation to a related approach by Shah et al., and discuss foreseeable applications to more generic orbits in Kerr spacetime.
△ Less
Submitted 7 August, 2015; v1 submitted 11 October, 2014;
originally announced October 2014.
-
Gravitational Self-Force Correction to the Innermost Stable Circular Equatorial Orbit of a Kerr Black Hole
Authors:
Soichiro Isoyama,
Leor Barack,
Sam R. Dolan,
Alexandre Le Tiec,
Hiroyuki Nakano,
Abhay G. Shah,
Takahiro Tanaka,
Niels Warburton
Abstract:
For a self-gravitating particle of mass μin orbit around a Kerr black hole of mass M >> μ, we compute the O(μ/M) shift in the frequency of the innermost stable circular equatorial orbit (ISCEO) due to the conservative piece of the gravitational self-force acting on the particle. Our treatment is based on a Hamiltonian formulation of the dynamics in terms of geodesic motion in a certain locally-def…
▽ More
For a self-gravitating particle of mass μin orbit around a Kerr black hole of mass M >> μ, we compute the O(μ/M) shift in the frequency of the innermost stable circular equatorial orbit (ISCEO) due to the conservative piece of the gravitational self-force acting on the particle. Our treatment is based on a Hamiltonian formulation of the dynamics in terms of geodesic motion in a certain locally-defined effective smooth spacetime. We recover the same result using the so-called first law of binary black-hole mechanics. We give numerical results for the ISCEO frequency shift as a function of the black hole's spin amplitude, and compare with predictions based on the post-Newtonian approximation and the effective one-body model. Our results provide an accurate strong-field benchmark for spin effects in the general relativistic two-body problem.
△ Less
Submitted 16 October, 2014; v1 submitted 24 April, 2014;
originally announced April 2014.
-
Gravitational-wave flux for a particle orbiting a Kerr black hole to 20th post-Newtonian order: a numerical approach
Authors:
Abhay G Shah
Abstract:
In this article we present the post-Newtonian (pN) coefficients of the energy flux (and angular momentum flux) at infinity and event horizon for a particle in circular, equatorial orbits about a Kerr black hole (of mass $M$ and spin-parameter $a$) up to 20-pN order. When a pN term is not a polynomial in $a/M$ and includes irrational functions (like polygamma functions), it is written as a power se…
▽ More
In this article we present the post-Newtonian (pN) coefficients of the energy flux (and angular momentum flux) at infinity and event horizon for a particle in circular, equatorial orbits about a Kerr black hole (of mass $M$ and spin-parameter $a$) up to 20-pN order. When a pN term is not a polynomial in $a/M$ and includes irrational functions (like polygamma functions), it is written as a power series of $a/M$. This is achieved by calculating the fluxes numerically with an accuracy greater than 1 part in $10^{600}$. Such high accuracy allows us to extract analytical values of pN coefficients that are linear combinations of transcendentals like the Euler constant, logarithms of prime numbers and powers of $π$. We also present the 22-pN expansion (spin-independent pN expansion) of the ingoing energy flux at the event horizon for a particle in circular orbit about a Schwarzschild black hole.
△ Less
Submitted 24 May, 2014; v1 submitted 11 March, 2014;
originally announced March 2014.
-
Finding high-order analytic post-Newtonian parameters from a high-precision numerical self-force calculation
Authors:
Abhay G Shah,
John L Friedman,
Bernard F Whiting
Abstract:
We present a novel analytic extraction of high-order post-Newtonian (pN) parameters that govern quasi-circular binary systems. Coefficients in the pN expansion of the energy of a binary system can be found from corresponding coefficients in an extreme-mass-ratio inspiral (EMRI) computation of the change $ΔU$ in the redshift factor of a circular orbit at fixed angular velocity. Remarkably, by compu…
▽ More
We present a novel analytic extraction of high-order post-Newtonian (pN) parameters that govern quasi-circular binary systems. Coefficients in the pN expansion of the energy of a binary system can be found from corresponding coefficients in an extreme-mass-ratio inspiral (EMRI) computation of the change $ΔU$ in the redshift factor of a circular orbit at fixed angular velocity. Remarkably, by computing this essentially gauge-invariant quantity to accuracy greater than one part in $10^{225}$, and by assuming that a subset of pN coefficients are rational numbers or products of $π$ and a rational, we obtain the exact analytic coefficients. We find the previously unexpected result that the post-Newtonian expansion of $ΔU$ (and of the change $ΔΩ$ in the angular velocity at fixed redshift factor) have conservative terms at half-integral pN order beginning with a 5.5 pN term. This implies the existence of a corresponding 5.5 pN term in the expansion of the energy of a binary system.
Coefficients in the pN series that do not belong to the subset just described are obtained to accuracy better than 1 part in $10^{265-23n}$ at $n$th pN order. We work in a radiation gauge, finding the radiative part of the metric perturbation from the gauge-invariant Weyl scalar $ψ_0$ via a Hertz potential. We use mode-sum renormalization, and find high-order renormalization coefficients by matching a series in $L=\ell+1/2$ to the large-$L$ behavior of the expression for $ΔU$. The non-radiative parts of the perturbed metric associated with changes in mass and angular momentum are calculated in the Schwarzschild gauge.
△ Less
Submitted 11 March, 2014; v1 submitted 6 December, 2013;
originally announced December 2013.
-
EMRI corrections to the angular velocity and redshift factor of a mass in circular orbit about a Kerr black hole
Authors:
Abhay G. Shah,
John L. Friedman,
Tobias S. Keidl
Abstract:
This is the first of two papers on computing the self-force in a radiation gauge for a particle moving in circular, equatorial orbit about a Kerr black hole. In the EMRI (extreme-mass-ratio inspiral) framework, with mode-sum renormalization, we compute the renormalized value of the quantity $h_{αβ}u^αu^β$, gauge-invariant under gauge transformations generated by a helically symmetric gauge vector;…
▽ More
This is the first of two papers on computing the self-force in a radiation gauge for a particle moving in circular, equatorial orbit about a Kerr black hole. In the EMRI (extreme-mass-ratio inspiral) framework, with mode-sum renormalization, we compute the renormalized value of the quantity $h_{αβ}u^αu^β$, gauge-invariant under gauge transformations generated by a helically symmetric gauge vector; and we find the related order $\frak{m}$ correction to the particle's angular velocity at fixed renormalized redshift (and to its redshift at fixed angular velocity). The radiative part of the perturbed metric is constructed from the Hertz potential which is extracted from the Weyl scalar by an algebraic inversion\cite{sf2}. We then write the spin-weighted spheroidal harmonics as a sum over spin-weighted spherical harmonics and use mode-sum renormalization to find the renormalization coefficients by matching a series in $L=\ell+1/2$ to the large-$L$ behavior of the expression for $H := \frac12 h_{αβ}u^αu^β$. The non-radiative parts of the perturbed metric associated with changes in mass and angular momentum are calculated in the Kerr gauge.
△ Less
Submitted 26 September, 2012; v1 submitted 24 July, 2012;
originally announced July 2012.
-
Gravitational Self-force in a Radiation Gauge
Authors:
Tobias S. Keidl,
Abhay G. Shah,
John L. Friedman,
Dong-Hoon Kim,
Larry R. Price
Abstract:
In this, the first of two companion papers, we present a method for finding the gravitational self-force in a modified radiation gauge for a particle moving on a geodesic in a Schwarzschild or Kerr spacetime. An extension of an earlier result by Wald is used to show the spin-weight $\pm 2$ perturbed Weyl scalar ($ψ_0$ or $ψ_4$) determines the metric perturbation outside the particle up to a gauge…
▽ More
In this, the first of two companion papers, we present a method for finding the gravitational self-force in a modified radiation gauge for a particle moving on a geodesic in a Schwarzschild or Kerr spacetime. An extension of an earlier result by Wald is used to show the spin-weight $\pm 2$ perturbed Weyl scalar ($ψ_0$ or $ψ_4$) determines the metric perturbation outside the particle up to a gauge transformation and an infinitesimal change in mass and angular momentum. A Hertz potential is used to construct the part of the retarded metric perturbation that involves no change in mass or angular momentum from $ψ_0$ in a radiation gauge. The metric perturbation is completed by adding changes in the mass and angular momentum of the background spacetime outside the radial coordinate $r_0$ of the particle in any convenient gauge. The resulting metric perturbation is singular on the trajectory of the particle and discontinuous across the sphere $r=r_0$. A mode-sum method can be used to renormalize the self-force, but the justification given in the published version of this paper \cite{sf2} referred to work by Sam Gralla \cite{gralla10} to justify the use of the renormalized self-force, and the radiation gauge we use does not satisfy the regularity conditions required by Gralla. Instead we show that the renormalized self-force, computed either from the retarded field for $r>r_0$ or for $r<r_0$ gives the correct equations of motion in a gauge smoothly related to a Lorenz gauge; and Pound et al. \cite{pmb13} argue that the average of the self-force obtained in the way described in our paper for $r>r_0$ and for $r<r_0$ gives the correct equation of motion for our gauge (what Pound et al. call the no-string gauge).
△ Less
Submitted 22 October, 2014; v1 submitted 13 April, 2010;
originally announced April 2010.