BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Leor Barack

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-phmath.NANumerical Analysis

Catalog footprint

What is connected

26works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

26 published item(s)

preprint2022arXiv

Time-domain metric reconstruction for hyperbolic scattering

Self-force methods can be applied in calculations of the scatter angle in two-body hyperbolic encounters, working order by order in the mass ratio (assumed small) but with no recourse to a weak-field approximation. This, in turn, can inform ongoing efforts to construct an accurate model of the general-relativistic binary dynamics via an effective-one-body description and other semi-analytical approaches. Existing self-force methods are to a large extent specialised to bound, inspiral orbits. Here we develop a technique for (numerical) self-force calculations that can efficiently tackle scatter orbits. The method is based on a time-domain reconstruction of the metric perturbation from a scalar-like Hertz potential that satisfies the Teukolsky equation, an idea pursued so far only for bound orbits. The crucial ingredient in this formulation are certain jump conditions that (each multipole mode of) the Hertz potential must satisfy along the orbit, in a 1+1-dimensional multipole reduction of the problem. We obtain a closed-form expression for these jumps, for an arbitrary geodesic orbit in Schwarzschild spacetime, and present a full numerical implementation for a scatter orbit. In this paper we focus on method development, and go only as far as calculating the Hertz potential; a calculation of the self-force and its physical effects on the scatter orbit will be the subject of forthcoming work.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Discontinuous collocation methods and gravitational self-force applications

Numerical simulations of extereme mass ratio inspirals, the mostimportant sources for the LISA detector, face several computational challenges. We present a new approach to evolving partial differential equations occurring in black hole perturbation theory and calculations of the self-force acting on point particles orbiting supermassive black holes. Such equations are distributionally sourced, and standard numerical methods, such as finite-difference or spectral methods, face difficulties associated with approximating discontinuous functions. However, in the self-force problem we typically have access to full a-priori information about the local structure of the discontinuity at the particle. Using this information, we show that high-order accuracy can be recovered by adding to the Lagrange interpolation formula a linear combination of certain jump amplitudes. We construct discontinuous spatial and temporal discretizations by operating on the corrected Lagrange formula. In a method-of-lines framework, this provides a simple and efficient method of solving time-dependent partial differential equations, without loss of accuracy near moving singularities or discontinuities. This method is well-suited for the problem of time-domain reconstruction of the metric perturbation via the Teukolsky or Regge-Wheeler-Zerilli formalisms. Parallel implementations on modern CPU and GPU architectures are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Completion of metric reconstruction for a particle orbiting a Kerr black hole

Vacuum perturbations of the Kerr metric can be reconstructed from the corresponding perturbation in either of the two Weyl scalars $ψ_0$ or $ψ_4$, using a procedure described by Chrzanowski and others in the 1970s. More recent work, motivated within the context of self-force physics, extends the procedure to metric perturbations sourced by a particle in a bound geodesic orbit. However, the existing procedure leaves undetermined a certain stationary, axially-symmetric piece of the metric perturbation. In the vacuum region away from the particle, this "completion" piece corresponds simply to mass and angular-momentum perturbations of the Kerr background, with amplitudes that are, however, a priori unknown. Here we present and implement a rigorous method for finding the completion piece. The key idea is to impose continuity, off the particle, of certain gauge-invariant fields constructed from the full (completed) perturbation, in order to determine the unknown amplitude parameters of the completion piece. We implement this method in full for bound (eccentric) geodesic orbits in the equatorial plane of the Kerr black hole. Our results provide a rigorous underpinning of recent results by Friedman {\it et al.}\ for circular orbits, and extend them to non-circular orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comparison Between Self-Force and Post-Newtonian Dynamics: Beyond Circular Orbits

The gravitational self-force (GSF) and post-Newtonian (PN) schemes are complementary approximation methods for modelling the dynamics of compact binary systems. Comparison of their results in an overlapping domain of validity provides a crucial test for both methods, and can be used to enhance their accuracy, e.g. via the determination of previously unknown PN parameters. Here, for the first time, we extend such comparisons to noncircular orbits---specifically, to a system of two nonspinning objects in a bound (eccentric) orbit. To enable the comparison we use a certain orbital-averaged quantity $\langle U \rangle$ that generalizes Detweiler's redshift invariant. The functional relationship $\langle U \rangle(Ω_r,Ω_ϕ)$, where $Ω_r$ and $Ω_ϕ$ are the frequencies of the radial and azimuthal motions, is an invariant characteristic of the conservative dynamics. We compute $\langle U \rangle(Ω_r,Ω_ϕ)$ numerically through linear order in the mass ratio $q$, using a GSF code which is based on a frequency-domain treatment of the linearized Einstein equations in the Lorenz gauge. We also derive $\langle U \rangle(Ω_r,Ω_ϕ)$ analytically through 3PN order, for an arbitrary $q$, using the known near-zone 3PN metric and the generalized quasi-Keplerian representation of the motion. We demonstrate that the $\mathcal{O}(q)$ piece of the analytical PN prediction is perfectly consistent with the numerical GSF results, and we use the latter to estimate yet unknown pieces of the 4PN expression at $\mathcal{O}(q)$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gravitational self-force on a particle in eccentric orbit around a Schwarzschild black hole

We present a numerical code for calculating the local gravitational self-force acting on a pointlike particle in a generic (bound) geodesic orbit around a Schwarzschild black hole. The calculation is carried out in the Lorenz gauge: For a given geodesic orbit, we decompose the Lorenz-gauge metric perturbation equations (sourced by the delta-function particle) into tensorial harmonics, and solve for each harmonic using numerical evolution in the time domain (in 1+1 dimensions). The physical self-force along the orbit is then obtained via mode-sum regularization. The total self-force contains a dissipative piece as well as a conservative piece, and we describe a simple method for disentangling these two pieces in a time-domain framework. The dissipative component is responsible for the loss of orbital energy and angular momentum through gravitational radiation; as a test of our code we demonstrate that the work done by the dissipative component of the computed force is precisely balanced by the asymptotic fluxes of energy and angular momentum, which we extract independently from the wave-zone numerical solutions. The conservative piece of the self force does not affect the time-averaged rate of energy and angular-momentum loss, but it influences the evolution of the orbital phases; this piece is calculated here for the first time in eccentric strong-field orbits. As a first concrete application of our code we recently reported the value of the shift in the location and frequency of the innermost stable circular orbit due to the conservative self-force [Phys. Rev. Lett.\ {\bf 102}, 191101 (2009)]. Here we provide full details of this analysis, and discuss future applications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Overspinning a Kerr black hole: the effect of self-force

We study the scenario in which a massive particle is thrown into a rapidly rotating Kerr black hole in an attempt to spin it up beyond its extremal limit, challenging weak cosmic censorship. We work in black-hole perturbation theory, and focus on non-spinning, uncharged particles sent in on equatorial orbits. We first identify the complete parameter-space region in which overspinning occurs when back-reaction effects from the particle's self-gravity are ignored. We find, in particular, that overspinning can be achieved only with particles sent in from infinity. Gravitational self-force effects may prevent overspinning by radiating away a sufficient amount of the particle's angular momentum ("dissipative effect"), and/or by increasing the effective centrifugal repulsion, so that particles with suitable parameters never get captured ("conservative effect"). We analyze the full effect of the self-force, thereby completing previous studies by Jacobson and Sotiriou (who neglected the self-force) and by Barausse, Cardoso and Khanna (who considered the dissipative effect on a subset of orbits). Our main result is an inequality, involving certain self-force quantities, which describes a necessary and sufficient condition for the overspinning scenario to be overruled. This "censorship" condition is formulated on a certain one-parameter family of geodesics in an extremal Kerr geometry. We find that the censorship condition is insensitive to the dissipative effect (within the first-order self-force approximation used here), except for a subset of perfectly fine-tuned orbits, for which a separate censorship condition is derived. We do not obtain here the self-force input needed to evaluate either of our two conditions, but discuss the prospects for producing the necessary data using state-of-the-art numerical codes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self force on a scalar charge in Kerr spacetime: circular equatorial orbits

We present a calculation of the scalar field self-force (SSF) acting on a scalar-charge particle in a strong-field orbit around a Kerr black hole. Our calculation specializes to circular and equatorial geodesic orbits. The analysis is an implementation of the standard mode-sum regularization scheme: We first calculate the multipole modes of the scalar-field perturbation using numerical integration in the frequency domain, and then apply a certain regularization procedure to each of the modes. The dissipative piece of the SSF is found to be consistent with the flux of energy and angular momentum carried by the scalar waves through the event horizon and out to infinity. The conservative (radial) component of the SSF is found to be attractive (inward pointing) for $r_0>r_{\rm c}(a)$ and repulsive (outward pointing) for $r_0<r_{\rm c}(a)$, where $a$ is the Kerr spin parameter, $r_0$ is the Boyer-Lindquist orbital radius, and $r_{\rm c}$ is a critical $a$-dependent radius at which the conservative SSF vanishes. When the motion is retrograde the conservative SSF is repulsive for all $r_0$ (as in the Schwarzschild case). The dominant conservative effect of the SSF in Schwarzschild spacetime is known to be of 3rd post-Newtonian (PN) order (with a logarithmic running). Our numerical results suggest that the leading-order PN correction due to the black hole's spin arises from spin-orbit coupling at 3PN, which dominates the overall SSF effect at large $r_0$. In PN language, the change-of-sign of the radial SSF is attributed to an interplay between the spin-orbit term ($\propto -ar_0^{-4.5}$) and the "Schwarzschild" term ($\propto r_0^{-5}\log r_0$).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self force on a scalar charge in Kerr spacetime: eccentric equatorial orbits

We present a numerical code for calculating the self force on a scalar charge moving in a bound (eccentric) geodesic in the equatorial plane of a Kerr black hole. We work in the frequency domain and make use of the method of extended homogeneous solutions [Phys.\ Rev.\ D {\bf 78}, 084021 (2008)], in conjunction with mode-sum regularization. Our work is part of a program to develop a computational architecture for fast and efficient self-force calculations, alternative to time-domain methods. We find that our frequency-domain method outperforms existing time-domain schemes for small eccentricities, and, remarkably, remains competitive up to eccentricities as high as $\sim 0.7$. As an application of our code we (i) compute the conservative scalar-field self-force correction to the innermost stable circular equatorial orbit, as a function of the Kerr spin parameter; and (ii) calculate the variation in the rest mass of the scalar particle along the orbit, caused by the component of the self force tangent to the four-velocity.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational self-torque and spin precession in compact binaries

We calculate the effect of self-interaction on the "geodetic" spin precession of a compact body in a strong-field orbit around a black hole. Specifically, we consider the spin precession angle $ψ$ per radian of orbital revolution for a particle carrying mass $μ$ and spin $s \ll (G/c) μ^2$ in a circular orbit around a Schwarzschild black hole of mass $M \gg μ$. We compute $ψ$ through $O(μ/M)$ in perturbation theory, i.e, including the correction $δψ$ (obtained numerically) due to the torque exerted by the conservative piece of the gravitational self-field. Comparison with a post-Newtonian (PN) expression for $δψ$, derived here through 3PN order, shows good agreement but also reveals strong-field features which are not captured by the latter approximation. Our results can inform semi-analytical models of the strong-field dynamics in astrophysical binaries, important for ongoing and future gravitational-wave searches.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Frequency-domain algorithm for the Lorenz-gauge gravitational self-force

State-of-the-art computations of the gravitational self-force (GSF) on massive particles in black hole spacetimes involve numerical evolution of the metric perturbation equations in the time-domain, which is computationally very costly. We present here a new strategy, based on a frequency-domain treatment of the perturbation equations, which offers considerable computational saving. The essential ingredients of our method are (i) a Fourier-harmonic decomposition of the Lorenz-gauge metric perturbation equations and a numerical solution of the resulting coupled set of ordinary equations with suitable boundary conditions; (ii) a generalized version of the method of extended homogeneous solutions [Phys. Rev. D {\bf 78}, 084021 (2008)] used to circumvent the Gibbs phenomenon that would otherwise hamper the convergence of the Fourier mode-sum at the particle's location; and (iii) standard mode-sum regularization, which finally yields the physical GSF as a sum over regularized modal contributions. We present a working code that implements this strategy to calculate the Lorenz-gauge GSF along eccentric geodesic orbits around a Schwarzschild black hole. The code is far more efficient than existing time-domain methods; the gain in computation speed (at a given precision) is about an order of magnitude at an eccentricity of 0.2, and up to three orders of magnitude for circular or nearly circular orbits. This increased efficiency was crucial in enabling the recently reported calculation of the long-term orbital evolution of an extreme mass ratio inspiral [Phys. Rev. D {\bf 85}, 061501(R) (2012)]. Here we provide full technical details of our method to complement the above report.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitational self-force from radiation-gauge metric perturbations

Calculations of the gravitational self-force (GSF) in curved spacetime require as input the metric perturbation in a sufficiently regular gauge. A basic challenge in the program to compute the GSF for orbits around a Kerr black hole is that the standard procedure for reconstructing the perturbation is formulated in a class of radiation gauges, in which the particle singularity is non-isotropic and extends away from the particle's location. Here we present two practical schemes for calculating the GSF using a radiation-gauge reconstructed metric as input. The schemes are based on a detailed analysis of the local structure of the particle singularity in the radiation gauges. We identify 3 types of radiation gauges: two containing a radial string-like singularity emanating from the particle, either in one direction ("half-string" gauges) or both directions ("full-string" gauges); and a third type containing no strings but with a jump discontinuity across a surface intersecting the particle. Based on a flat-space example, we argue that the standard mode-by-mode reconstruction procedure yields the "regular half" of a half-string solution, or (equivalently) either of the regular halves of a no-string solution. For the half-string case, we formulate the GSF in a locally deformed radiation gauge that removes the string singularity near the particle. We derive a mode-sum formula for the GSF in this gauge, analogous to the standard Lorenz-gauge formula but with modified regularization parameters. For the no-string case, we formulate the GSF directly, without a local deformation, and we derive a mode-sum formula that requires no correction to the parameters but involves a certain averaging procedure. We explain the consistency of our results with Gralla's invariance theorem, and discuss the correspondence between our method and a related approach by Friedman et al.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Isofrequency pairing of geodesic orbits in Kerr geometry

Bound geodesic orbits around a Kerr black hole can be parametrized by three constants of the motion: the (specific) orbital energy, angular momentum and Carter constant. Generically, each orbit also has associated with it three frequencies, related to the radial, longitudinal and (mean) azimuthal motions. Here we note the curious fact that these two ways of characterizing bound geodesics are not in a one-to-one correspondence. While the former uniquely specifies an orbit up to initial conditions, the latter does not: there is a (strong-field) region of the parameter space in which pairs of physically distinct orbits can have the same three frequencies. In each such isofrequency pair the two orbits exhibit the same rate of periastron precession and the same rate of Lense-Thirring precession of the orbital plane, and (in a certain sense) they remain "synchronized" in phase.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Critical phenomena at the threshold of immediate merger in binary black hole systems: the extreme mass ratio case

In numerical simulations of black hole binaries, Pretorius and Khurana [Class. Quant. Grav. {\bf 24}, S83 (2007)] have observed critical behaviour at the threshold between scattering and immediate merger. The number of orbits scales as $n\simeq -γ\ln|p-p_*|$ along any one-parameter family of initial data such that the threshold is at $p=p_*$. Hence they conjecture that in ultrarelavistic collisions almost all the kinetic energy can be converted into gravitational waves if the impact parameter is fine-tuned to the threshold. As a toy model for the binary, they consider the geodesic motion of a test particle in a Kerr black hole spacetime, where the unstable circular geodesics play the role of critical solutions, and calculate the critical exponent $γ$. Here, we incorporate radiation reaction into this model using the self-force approximation. The critical solution now evolves adiabatically along a sequence of unstable circular geodesic orbits under the effect of the self-force. We confirm that almost all the initial energy and angular momentum are radiated on the critical solution. Our calculation suggests that, even for infinite initial energy, this happens over a finite number of orbits given by $n_\infty\simeq 0.41/η$, where $η$ is the (small) mass ratio. We derive expressions for the time spent on the critical solution, number of orbits and radiated energy as functions of the initial energy and impact parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational self-force and the effective-one-body formalism between the innermost stable circular orbit and the light ring

We compute the conservative piece of the gravitational self-force (GSF) acting on a particle of mass m_1 as it moves along an (unstable) circular geodesic orbit between the innermost stable circular orbit (ISCO) and the light ring of a Schwarzschild black hole of mass m_2>> m_1. More precisely, we construct the function h_{uu}(x) = h_{μν} u^μ u^ν (related to Detweiler's gauge-invariant "redshift" variable), where h_{μν} is the regularized metric perturbation in the Lorenz gauge, u^μ is the four-velocity of m_1, and x= [Gc^{-3}(m_1+m_2)Ω]^{2/3} is an invariant coordinate constructed from the orbital frequency Ω. In particular, we explore the behavior of h_{uu} just outside the "light ring" at x=1/3, where the circular orbit becomes null. Using the recently discovered link between h_{uu} and the piece a(u), linear in the symmetric mass ratio ν, of the main radial potential A(u,ν) of the Effective One Body (EOB) formalism, we compute a(u) over the entire domain 0<u<1/3. We find that a(u) diverges at the light-ring as ~0.25 (1-3u)^{-1/2}, explain the physical origin of this divergence, and discuss its consequences for the EOB formalism. We construct accurate global analytic fits for a(u), valid on the entire domain 0<u<1/3 (and possibly beyond), and give accurate numerical estimates of the values of a(u) and its first 3 derivatives at the ISCO, as well as the O(ν) shift in the ISCO frequency. In previous work we used GSF data on slightly eccentric orbits to compute a certain linear combination of a(u) and its first two derivatives, involving also the O(ν) piece \bar d(u) of a second EOB radial potential {\bar D}(u,ν). Combining these results with our present global analytic representation of a(u), we numerically compute {\bar d}(u)$ on the interval 0<u\leq 1/6.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Self-force via $m$-mode regularization and 2+1D evolution: III. Gravitational field on Schwarzschild spacetime

This is the third in a series of papers aimed at developing a practical time-domain method for self-force calculations in Kerr spacetime. The key elements of the method are (i) removal of a singular part of the perturbation field with a suitable analytic "puncture", (ii) decomposition of the perturbation equations in azimuthal ($m$-)modes, taking advantage of the axial symmetry of the Kerr background, (iii) numerical evolution of the individual $m$-modes in 2+1-dimensions with a finite difference scheme, and (iv) reconstruction of the local self-force from the mode sum. Here we report a first implementation of the method to compute the gravitational self-force. We work in the Lorenz gauge, solving directly for the metric perturbation in 2+1-dimensions. The modes $m=0,1$ contain nonradiative pieces, whose time-domain evolution is hampered by certain gauge instabilities. We study this problem in detail and propose ways around it. In the current work we use the Schwarzschild geometry as a platform for development; in a forthcoming paper---the fourth in the series---we apply our method to the gravitational self-force in Kerr geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Beyond the geodesic approximation: conservative effects of the gravitational self-force in eccentric orbits around a Schwarzschild black hole

We study conservative finite-mass corrections to the motion of a particle in a bound (eccentric) strong-field orbit around a Schwarzschild black hole. We assume the particle's mass $μ$ is much smaller than the black hole mass $M$, and explore post-geodesic corrections of $O(μ/M)$. Our analysis uses numerical data from a recently developed code that outputs the Lorenz-gauge gravitational self-force (GSF) acting on the particle along the eccentric geodesic. First, we calculate the $O(μ/M)$ conservative correction to the periastron advance of the orbit, as a function of the (gauge-dependent) semilatus rectum and eccentricity. A gauge-invariant description of the GSF precession effect is made possible in the circular-orbit limit, where we express the correction to the periastron advance as a function of the invariant azimuthal frequency. We compare this relation with results from fully nonlinear numerical-relativistic simulations. In order to obtain a gauge-invariant measure of the GSF effect for fully eccentric orbits, we introduce a suitable generalization of Detweiler's circular-orbit "redshift" invariant. We compute the $O(μ/M)$ conservative correction to this invariant, expressed as a function of the two invariant frequencies that parametrize the orbit. Our results are in good agreement with results from post-Newtonian calculations in the weak-field regime, as we shall report elsewhere. The results of our study can inform the development of analytical models for the dynamics of strongly gravitating binaries. They also provide an accurate benchmark for future numerical-relativistic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Periastron Advance in Black-Hole Binaries

The general relativistic (Mercury-type) periastron advance is calculated here for the first time with exquisite precision in full general relativity. We use accurate numerical relativity simulations of spinless black hole binaries with mass ratios 1/8 < m1/m2 < 1 and compare with the predictions of several analytic approximation schemes. We find the effective-one-body model to be remarkably accurate, and, surprisingly, so also the predictions of self-force theory [replacing m1/m2 --> m1m2/(m1+m2)^2]. Our results can inform a universal analytic model of the two-body dynamics, crucial for ongoing and future gravitational-wave searches.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self force via $m$-mode regularization and 2+1D evolution: II. Scalar-field implementation on Kerr spacetime

This is the second in a series of papers aimed at developing a practical time-domain method for self-force calculations in Kerr spacetime. The key elements of the method are (i) removal of a singular part of the perturbation field with a suitable analytic "puncture" based on the Detweiler--Whiting decomposition, (ii) decomposition of the perturbation equations in azimuthal ($m$-)modes, taking advantage of the axial symmetry of the Kerr background, (iii) numerical evolution of the individual $m$-modes in 2+1-dimensions with a finite difference scheme, and (iv) reconstruction of the physical self-force from the mode sum. Here we report an implementation of the method to compute the scalar-field self-force along circular equatorial geodesic orbits around a Kerr black hole. This constitutes a first time-domain computation of the self force in Kerr geometry. Our time-domain code reproduces the results of a recent frequency-domain calculation by Warburton and Barack, but has the added advantage of being readily adaptable to include the back-reaction from the self force in a self-consistent manner. In a forthcoming paper---the third in the series---we apply our method to the gravitational self-force (in the Lorenz gauge).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gravitational self force on a particle in circular orbit around a Schwarzschild black hole

We calculate the gravitational self force acting on a pointlike particle of mass $μ$, set in a circular geodesic orbit around a Schwarzschild black hole. Our calculation is done in the Lorenz gauge: For given orbital radius, we first solve directly for the Lorenz-gauge metric perturbation using numerical evolution in the time domain; We then compute the (finite) back-reaction force from each of the multipole modes of the perturbation; Finally, we apply the ``mode sum'' method to obtain the total, physical self force. The {\em temporal} component of the self force (which is gauge invariant) describes the dissipation of orbital energy through gravitational radiation. Our results for this component are consistent, to within the computational accuracy, with the total flux of gravitational-wave energy radiated to infinity and through the event horizon. The {\em radial} component of the self force (which is gauge dependent) is calculated here for the first time. It describes a conservative shift in the orbital parameters away from their geodesic values. We thus obtain the $O(μ)$ correction to the specific energy and angular momentum parameters (in the Lorenz gauge), as well as the $O(μ)$ shift in the orbital frequency (which is gauge invariant).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precession effect of the gravitational self-force in a Schwarzschild spacetime and the effective one-body formalism

Using a recently presented numerical code for calculating the Lorenz-gauge gravitational self-force (GSF), we compute the $O(m)$ conservative correction to the precession rate of the small-eccentricity orbits of a particle of mass $m$ moving around a Schwarzschild black hole of mass ${\mathsf M}\gg m$. Specifically, we study the gauge-invariant function $ρ(x)$, where $ρ$ is defined as the $O(m)$ part of the dimensionless ratio $(\hatΩ_r/\hatΩ_φ)^2$ between the squares of the radial and azimuthal frequencies of the orbit, and where $x=[Gc^{-3}({\mathsf M}+m)\hatΩ_φ]^{2/3}$ is a gauge-invariant measure of the dimensionless gravitational potential (mass over radius) associated with the mean circular orbit. Our GSF computation of the function $ρ(x)$ in the interval $0<x\leq 1/6$ determines, for the first time, the {\em strong-field behavior} of a combination of two of the basic functions entering the Effective One Body (EOB) description of the conservative dynamics of binary systems. We show that our results agree well in the weak-field regime (small $x$) with the 3rd post-Newtonian (PN) expansion of the EOB results, and that this agreement is improved when taking into account the analytic values of some of the logarithmic-running terms occurring at higher PN orders. Furthermore, we demonstrate that GSF data give access to higher-order PN terms of $ρ(x)$ and can be used to set useful new constraints on the values of yet-undetermined EOB parameters. Most significantly, we observe that an {\em excellent global representation} of $ρ(x)$ can be obtained using a simple `two-point' Padé approximant which combines 3PN knowledge at $x=0$ with GSF information at a single strong-field point (say, $x=1/6$).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Self force via m-mode regularization and 2+1D evolution: Foundations and a scalar-field implementation on Schwarzschild

To model the radiative evolution of extreme mass-ratio binary inspirals (a key target of the LISA mission), the community needs efficient methods for computation of the gravitational self-force (SF) on the Kerr spacetime. Here we further develop a practical `$m$-mode regularization' scheme for SF calculations, and give details of a first implementation. The key steps in the method are (i) removal of a singular part of the perturbation field with a suitable `puncture' to leave a sufficiently regular residual within a finite worldtube surrounding the particle's worldline, (ii) decomposition in azimuthal ($m$-)modes, (iii) numerical evolution of the $m$-modes in 2+1D with a finite difference scheme, and (iv) reconstruction of the SF from the mode sum. The method relies on a judicious choice of puncture, based on the Detweiler--Whiting decomposition. We give a working definition for the `order' of the puncture, and show how it determines the convergence rate of the $m$-mode sum. The dissipative piece of the SF displays an exponentially convergent mode sum, while the $m$-mode sum for the conservative piece converges with a power law. In the latter case the individual modal contributions fall off at large $m$ as $m^{-n}$ for even $n$ and as $m^{-n+1}$ for odd $n$, where $n$ is the puncture order. We describe an $m$-mode implementation with a 4th-order puncture to compute the scalar-field SF along circular geodesics on Schwarzschild. In a forthcoming companion paper we extend the calculation to the Kerr spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Gravitational self force in extreme mass-ratio inspirals

This review is concerned with the gravitational self-force acting on a mass particle in orbit around a large black hole. Renewed interest in this old problem is driven by the prospects of detecting gravitational waves from strongly gravitating binaries with extreme mass ratios. We begin here with a summary of recent advances in the theory of gravitational self-interaction in curved spacetime, and proceed to survey some of the ideas and computational strategies devised for implementing this theory in the case of a particle orbiting a Kerr black hole. We review in detail two of these methods: (i) the standard mode-sum method, in which the metric perturbation is regularized mode-by-mode in a multipole decomposition, and (ii) $m$-mode regularization, whereby individual azimuthal modes of the metric perturbation are regularized in 2+1 dimensions. We discuss several practical issues that arise, including the choice of gauge, the numerical representation of the particle singularity, and how high-frequency contributions near the particle are dealt with in frequency-domain calculations. As an example of a full end-to-end implementation of the mode-sum method, we discuss the computation of the gravitational self-force for eccentric geodesic orbits in Schwarzschild, using a direct integration of the Lorenz-gauge perturbation equations in the time domain. With the computational framework now in place, researchers have recently turned to explore the physical consequences of the gravitational self force; we will describe some preliminary results in this area. An appendix to this review presents, for the first time, a detailed derivation of the regularization parameters necessary for implementing the mode-sum method in Kerr spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Two approaches for the gravitational self force in black hole spacetime: Comparison of numerical results

Recently, two independent calculations have been presented of finite-mass ("self-force") effects on the orbit of a point mass around a Schwarzschild black hole. While both computations are based on the standard mode-sum method, they differ in several technical aspects, which makes comparison between their results difficult--but also interesting. Barack and Sago [Phys. Rev. D {\bf 75}, 064021 (2007)] invoke the notion of a self-accelerated motion in a background spacetime, and perform a direct calculation of the local self force in the Lorenz gauge (using numerical evolution of the perturbation equations in the time domain); Detweiler [Phys. Rev. D {\bf 77}, 124026 (2008)] describes the motion in terms a geodesic orbit of a (smooth) perturbed spacetime, and calculates the metric perturbation in the Regge--Wheeler gauge (using frequency-domain numerical analysis). Here we establish a formal correspondence between the two analyses, and demonstrate the consistency of their numerical results. Specifically, we compare the value of the conservative $O(μ)$ shift in $u^t$ (where $μ$ is the particle's mass and $u^t$ is the Schwarzschild $t$ component of the particle's four-velocity), suitably mapped between the two orbital descriptions and adjusted for gauge. We find that the two analyses yield the same value for this shift within mere fractional differences of $\sim 10^{-5}$--$10^{-7}$ (depending on the orbital radius)--comparable with the estimated numerical error.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Self-force on a scalar particle in spherically-symmetric spacetime via mode-sum regularization: radial trajectories

Recently, we proposed a method for calculating the ``radiation reaction'' self-force exerted on a charged particle moving in a strong field orbit in a black hole spacetime. In this approach, one first calculates the contribution to the ``tail'' part of the self force due to each multipole mode of the particle's self field. A certain analytic procedure is then applied to regularize the (otherwise divergent) sum over modes. This involves the derivation of certain regularization parameters using local analysis of the (retarded) Green's function. In the present paper we present a detailed formulation of this mode-sum regularization scheme for a scalar charge on a class of static spherically-symmetric backgrounds (including, e.g., the Schwarzschild, Reissner-Nordström, and Schwarzschild-de Sitter spacetimes). We fully implement the regularization scheme for an arbitrary radial trajectory (not necessarily geodesic) by explicitly calculating all necessary regularization parameters in this case.

0 citations0 reviewsNo code linkedOpen access