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

Niels Warburton

Niels Warburton 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.HEhep-th

Catalog footprint

What is connected

20works
3topics
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

20 published item(s)

preprint2022arXiv

Eccentric self-forced inspirals into a rotating black hole

We develop the first model for extreme mass-ratio inspirals (EMRIs) into a rotating massive black hole driven by the gravitational self-force. Our model is based on an action angle formulation of the method of osculating geodesics for eccentric, equatorial (i.e., spin-aligned) motion in Kerr spacetime. The forcing terms are provided by an efficient spectral interpolation of the first-order gravitational self-force in the outgoing radiation gauge. We apply a near-identity (averaging) transformation to eliminate all dependence of the orbital phases from the equations of motion, while maintaining all secular effects of the first-order gravitational self-force at post-adiabatic order. This implies that the model can be evolved without having to resolve all $\mathcal{O}(10^5)$ orbit cycles of an EMRI, yielding an inspiral model that can be evaluated in less than a second for any mass-ratio. In the case of a non-rotating central black hole, we compare inspirals evolved using self-force data computed in the Lorenz and radiation gauges. We find that the two gauges generally produce differing inspirals with a deviation of comparable magnitude to the conservative self-force correction. This emphasizes the need for including the (currently unknown) dissipative second order self-force to obtain gauge independent, post-adiabatic waveforms.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hyperboloidal method for frequency-domain self-force calculations

Gravitational self-force theory is the leading approach for modeling gravitational wave emission from small mass-ratio compact binaries. This method perturbatively expands the metric of the binary in powers of the mass ratio. The source for the perturbations depends on the orbital configuration, calculational approach, and the order of the perturbative expansion. These sources fall into three broad classes: (i) distributional, (ii) worldtube, and (iii) unbounded support. The latter, in particular, is important for emerging second-order (in the mass ratio) calculations. Traditional frequency domain approaches employ the variation of parameters method and compute the perturbation on standard time slices with numerical boundary conditions supplied at finite radius from series expansions of the asymptotic behavior. This approach has been very successful, but the boundary conditions calculations are tedious, and the approach is not well suited to unbounded sources where homogeneous solutions must be computed at all radii. This work develops an alternative approach where hyperboloidal slices foliate the spacetime, and compactifying coordinates simplify the boundary treatment. We implement this approach with a multi-domain spectral solver with analytic mesh refinement and use the scalar-field self-force on circular orbits around a Schwarzschild black hole as an example problem. The method works efficiently for all three source classes encountered in self-force calculations and has distinct advantages over the traditional approach. For example, our code efficiently computes the perturbation for orbits with extremely large orbital radii ($r_{p}>10^5M$) or modes with very high spherical harmonic mode index ($\ell \ge 100$). Our results indicate that hyperboloidal methods can play an essential role in self-force calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

FastEMRIWaveforms: New tools for millihertz gravitational-wave data analysis

We present the FastEMRIWaveforms (FEW) package, a collection of tools to build and analyze extreme mass ratio inspiral (EMRI) waveforms. Here, we expand on the Physical Review Letter that introduced the first fast and accurate fully-relativistic EMRI waveform template model. We discuss the construction of the overall framework; constituent modules; and the general methods used to accelerate EMRI waveforms. Because the fully relativistic FEW model waveforms are for now limited to eccentric orbits in the Schwarzschild spacetime, we also introduce an improved Augmented Analytic Kludge (AAK) model that describes generic Kerr inspirals. Both waveform models can be accelerated using graphics processing unit (GPU) hardware. With the GPU-accelerated waveforms in hand, a variety of studies are performed including an analysis of EMRI mode content, template mismatch, and fully Bayesian Markov Chain Monte Carlo-based EMRI parameter estimation. We find relativistic EMRI waveform templates can be generated with fewer harmonic modes ($\sim10-100$) without biasing signal extraction. However, we show for the first time that extraction of a relativistic injection with semi-relativistic amplitudes can lead to strong bias and anomalous structure in the posterior distribution for certain regions of parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dissipation in extreme-mass ratio binaries with a spinning secondary

We present the gravitational-wave flux balance law in an extreme mass-ratio binary with a spinning secondary. This law relates the flux of energy (angular momentum) radiated to null infinity and through the event horizon to the local change in the secondary's orbital energy (angular momentum) for generic (non-resonant) bound orbits in Kerr spacetime. As an explicit example we compute these quantities for a spin-aligned body moving on a circular orbit around a Schwarzschild black hole. We perform this calculation both analytically, via a high-order post-Newtonian expansion, and numerically in two different gauges. Using these results we demonstrate explicitly that our new balance law holds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rapid generation of fully relativistic extreme-mass-ratio-inspiral waveform templates for LISA data analysis

The future space mission LISA will observe a wealth of gravitational-wave sources at millihertz frequencies. Of these, the extreme-mass-ratio inspirals of compact objects into massive black holes are the only sources that combine the challenges of strong-field complexity with that of long-lived signals. Such signals are found and characterized by comparing them against a large number of accurate waveform templates during data analysis, but the rapid generation of such templates is hindered by computing the $\sim10^3$-$10^5$ harmonic modes in a fully relativistic waveform. We use order-reduction and deep-learning techniques to derive a global fit for these modes, and implement it in a complete waveform framework with hardware acceleration. Our high-fidelity waveforms can be generated in under $1\,\mathrm{s}$, and achieve a mismatch of $\lesssim 5\times 10^{-4}$ against reference waveforms that take $\gtrsim 10^4$ times longer. This marks the first time that analysis-length waveforms with full harmonic content can be produced on timescales useful for direct implementation in LISA analysis algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The location of the last stable orbit in Kerr spacetime

Black hole spacetimes, like the Kerr spacetime, admit both stable and plunging orbits, separated in parameter space by the separatrix. Determining the location of the separatrix is of fundamental interest in understanding black holes, and is of crucial importance for modeling extreme mass-ratio inspirals. Previous numerical approaches to locating the Kerr separatrix were not always efficient or stable across all of parameter space. In this paper we show that the Kerr separatrix is the zero set of a single polynomial in parameter space. This gives two main results. First, we thoroughly analyze special cases (extreme Kerr, polar orbits, etc.), finding strict bounds on the limits of roots, and unifying a number of results in the literature. Second, we pose a stable numerical method which is guaranteed to quickly and robustly converge to the separatrix. This new approach is implemented in the Black Hole Perturbation Toolkit, and results in a ~45x speedup over the prior robust approach.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Second-order self-force calculation of the gravitational binding energy in compact binaries

Self-force theory is the leading method of modeling extreme-mass-ratio inspirals (EMRIs), key sources for the gravitational-wave detector LISA. It is well known that for an accurate EMRI model, second-order self-force effects are critical, but calculations of these effects have been beset by obstacles. In this letter we present the first implementation of a complete scheme for second-order self-force computations, specialized to the case of quasicircular orbits about a Schwarzschild black hole. As a demonstration, we calculate the gravitational binding energy of these binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Highly eccentric inspirals into a black hole

We model the inspiral of a compact stellar-mass object into a massive nonrotating black hole including all dissipative and conservative first-order-in-the-mass-ratio effects on the orbital motion. The techniques we develop allow inspirals with initial eccentricities as high as $e\sim0.8$ and initial separations as large as $p\sim 50$ to be evolved through many thousands of orbits up to the onset of the plunge into the black hole. The inspiral is computed using an osculating elements scheme driven by a hybridized self-force model, which combines Lorenz-gauge self-force results with highly accurate flux data from a Regge-Wheeler-Zerilli code. The high accuracy of our hybrid self-force model allows the orbital phase of the inspirals to be tracked to within $\sim0.1$ radians or better. The difference between self-force models and inspirals computed in the radiative approximation is quantified.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inspiral into Gargantua

We model the inspiral of a compact object into a more massive black hole rotating very near the theoretical maximum. We find that once the body enters the near-horizon regime the gravitational radiation is characterized by a constant frequency, equal to (twice) the horizon frequency, with an exponentially damped profile. This contrasts with the usual "chirping" behavior and, if detected, would constitute a "smoking gun" for a near-extremal black hole in nature.

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

Particle on the Innermost Stable Circular Orbit of a Rapidly Spinning Black Hole

We compute the radiation emitted by a particle on the innermost stable circular orbit of a rapidly spinning black hole both (a) analytically, working to leading order in the deviation from extremality and (b) numerically, with a new high-precision Teukolsky code. We find excellent agreement between the two methods. We confirm previous estimates of the overall scaling of the power radiated, but show that there are also small oscillations all the way to extremality. Furthermore, we reveal an intricate mode-by-mode structure in the flux to infinity, with only certain modes having the dominant scaling. The scaling of each mode is controlled by its conformal weight, a quantity that arises naturally in the representation theory of the enhanced near-horizon symmetry group. We find relationships to previous work on particles orbiting in precisely extreme Kerr, including detailed agreement of quantities computed here with conformal field theory calculations performed in the context of the Kerr/CFT correspondence.

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 on a scalar charge in Kerr spacetime: inclined circular orbits

Accurately modeling astrophysical extreme-mass-ratio-insprials requires calculating the gravitational self-force for orbits in Kerr spacetime. The necessary calculation techniques are typically very complex and, consequently, toy scalar-field models are often developed in order to establish a particular calculational approach. To that end, I present a calculation of the scalar-field self-force for a particle moving on a (fixed) inclined circular geodesic of a background Kerr black hole. I make the calculation in the frequency-domain and demonstrate how to apply the mode-sum regularization procedure to all four components of the self-force. I present results for a number of strong-field orbits which can be used as benchmarks for emerging self-force calculation techniques in Kerr spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tidal invariants for compact binaries on quasi-circular orbits

We extend the gravitational self-force approach to encompass `self-interaction' tidal effects for a compact body of mass $μ$ on a quasi-circular orbit around a black hole of mass $M \gg μ$. Specifically, we define and calculate at $O(μ)$ (conservative) shifts in the eigenvalues of the electric- and magnetic-type tidal tensors, and a (dissipative) shift in a scalar product between their eigenbases. This approach yields four gauge-invariant functions, from which one may construct other tidal quantities such as the curvature scalars and the speciality index. First, we analyze the general case of a geodesic in a regular perturbed vacuum spacetime admitting a helical Killing vector and a reflection symmetry. Next, we specialize to focus on circular orbits in the equatorial plane of Kerr spacetime at $O(μ)$. We present accurate numerical results for the Schwarzschild case for orbital radii up to the light-ring, calculated via independent implementations in Lorenz and Regge-Wheeler gauges. We show that our results are consistent with leading-order post-Newtonian expansions, and demonstrate the existence of additional structure in the strong-field regime. We anticipate that our strong-field results will inform (e.g.) effective one-body models for the gravitational two-body problem that are invaluable in the ongoing search for gravitational waves.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Applying the effective-source approach to frequency-domain self-force calculations

The equations of motion of a point particle interacting with its own field are defined in terms of a certain regularized self-field. Two of the leading methods for computing this regularized field are the mode-sum and effective-source approaches. In this work we unite these two distinct regularization schemes by generalizing traditional frequency-domain mode-sum calculations to incorporate effective-source techniques. For a toy scalar-field model we analytically compute an appropriate puncture field from which the regularized residual field can be calculated. To demonstrate the method, we compute the self-force for a scalar particle on a circular orbit in Schwarzschild spacetime. We also demonstrate the relation between the worldtube and window function approaches to localizing the puncture field to the neighborhood of the worldline and show how the method reduces to the well-known mode-sum regularization scheme in a certain limit. This new computational scheme can be applied to cases where traditional mode-sum regularization is inadequate, such as in calculations at second perturbative order.

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

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