Source author record

Abhay G. Shah

Abhay G. Shah appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2020arXiv

Metric perturbations produced by eccentric equatorial orbits around a Kerr black hole

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.

preprint2016arXiv

Gauge-invariant perturbations of Schwarzschild spacetime

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.

preprint2016arXiv

iResum: a new paradigm for resumming gravitational wave amplitudes

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.

preprint2016arXiv

Raising and Lowering operators of spin-weighted spheroidal harmonics

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.

preprint2015arXiv

Experimental mathematics meets gravitational self-force

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.

preprint2015arXiv

Linear-in-mass-ratio contribution to spin precession and tidal invariants in Schwarzschild spacetime at very high post-Newtonian order

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.

preprint2015arXiv

Self-force as a cosmic censor in the Kerr overspinning problem

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.

preprint2015arXiv

Self-force from reconstructed metric perturbations: numerical implementation in Schwarzschild spacetime

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.

preprint2014arXiv

Gravitational Self-Force Correction to the Innermost Stable Circular Equatorial Orbit of a Kerr Black Hole

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.

preprint2014arXiv

Gravitational Self-force in a Radiation Gauge

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).

preprint2012arXiv

EMRI corrections to the angular velocity and redshift factor of a mass in circular orbit about a Kerr black hole

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.