Researcher profile

Jonathan Thornburg

Jonathan Thornburg contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 published item(s)

preprint2019arXiv

Excitation of Kerr quasinormal modes in extreme--mass-ratio inspirals

If a small compact object orbits a black hole, it is known that it can excite the black hole's quasinormal modes (QNMs), leading to high-frequency oscillations ("wiggles") in the radiated field at $\mathcal{J}^+$, and in the radiation-reaction self-force acting on the object after its orbit passes through periapsis. Here we survey the phenomenology of these wiggles across a range of black hole spins and equatorial orbits. In both the scalar-field and gravitational cases we find that wiggles are a generic feature across a wide range of parameter space, and that they are observable in field perturbations at fixed spatial positions, in the self-force, and in radiated fields at $\mathcal{J}^+$. For a given charge or mass of the small body, the QNM excitations have the highest amplitudes for systems with a highly spinning central black hole, a prograde orbit with high eccentricity, and an orbital periapsis close to the light ring. The QNM amplitudes depend smoothly on the orbital parameters, with only very small amplitude changes when the orbit's (discrete) frequency spectrum is tuned to match QNM frequencies. The association of wiggles with QNM excitations suggest that they represent a situation where the \emph{nonlocal} nature of the self-force is particularly apparent, with the wiggles arising as result of QNM excitation by the compact object near periapsis, and then encountered later in the orbit. Astrophysically, the effects of wiggles at $\mathcal{J}^+$ might allow direct observation of Kerr QNMs in extreme-mass-ratio inspiral (EMRI) binary black hole systems, potentially enabling new tests of general relativity.

preprint2011arXiv

Gravitational Waves Notes, Issue #5 : "The Capra research programme for capture of small compact objects by massive black holes"

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - astrophysics, general relativity and data analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. In this 6th GW Note (since we started numbering with #0), we present the work of Jonathan Thornburg, who has been fully-refereed, on the Capra research programme for capture of small compact objects by massive black holes.

preprint2011arXiv

The Effect of Massive Perturbers on Extreme Mass-Ratio Inspiral Waveforms

Extreme mass ratio inspirals, in which a stellar-mass object merges with a supermassive black hole, are prime sources for space-based gravitational wave detectors because they will facilitate tests of strong gravity and probe the spacetime around rotating compact objects. In the last few years of such inspirals, the total phase is in the millions of radians and details of the waveforms are sensitive to small perturbations. We show that one potentially detectable perturbation is the presence of a second supermassive black hole within a few tenths of a parsec. The acceleration produced by the perturber on the extreme mass-ratio system produces a steady drift that causes the waveform to deviate systematically from that of an isolated system. If the perturber is a few tenths of a parsec from the extreme-mass ratio system (plausible in as many as a few percent of cases) higher derivatives of motion might also be detectable. In that case, the mass and distance of the perturber can be derived independently, which would allow a new probe of merger dynamics.

preprint2010arXiv

Highly accurate and efficient self-force computations using time-domain methods: Error estimates, validation, and optimization

If a small "particle" of mass $μM$ (with $μ\ll 1$) orbits a Schwarzschild or Kerr black hole of mass $M$, the particle is subject to an $Ø(μ)$ radiation-reaction "self-force". Here I argue that it's valuable to compute this self-force highly accurately (relative error of $\ltsim 10^{-6}$) and efficiently, and I describe techniques for doing this and for obtaining and validating error estimates for the computation. I use an adaptive-mesh-refinement (AMR) time-domain numerical integration of the perturbation equations in the Barack-Ori mode-sum regularization formalism; this is efficient, yet allows easy generalization to arbitrary particle orbits. I focus on the model problem of a scalar particle in a circular geodesic orbit in Schwarzschild spacetime. The mode-sum formalism gives the self-force as an infinite sum of regularized spherical-harmonic modes $\sum_{\ell=0}^\infty F_{\ell,\reg}$, with $F_{\ell,\reg}$ (and an "internal" error estimate) computed numerically for $\ell \ltsim 30$ and estimated for larger~$\ell$ by fitting an asymptotic "tail" series. Here I validate the internal error estimates for the individual $F_{\ell,\reg}$ using a large set of numerical self-force computations of widely-varying accuracies. I present numerical evidence that the actual numerical errors in $F_{\ell,\reg}$ for different~$\ell$ are at most weakly correlated, so the usual statistical error estimates are valid for computing the self-force. I show that the tail fit is numerically ill-conditioned, but this can be mostly alleviated by renormalizing the basis functions to have similar magnitudes. Using AMR, fixed mesh refinement, and extended-precision floating-point arithmetic, I obtain the (contravariant) radial component of the self-force for a particle in a circular geodesic orbit of areal radius $r = 10M$ to within $1$~ppm relative error.

preprint2007arXiv

Recoil Velocities from Equal-Mass Binary-Black-Hole Mergers

The final evolution of a binary black-hole system gives rise to a recoil velocity if an asymmetry is present in the emitted gravitational radiation. Measurements of this effect for non-spinning binaries with unequal masses have pointed out that kick velocities $\sim~175$ km/s can be reached for a mass ratio $\simeq 0.36$. However, a larger recoil can be obtained for equal-mass binaries if the asymmetry is provided by the spins. Using two independent methods we show that the merger of such binaries yields velocities as large as $\sim 440$ km/s for black holes having unequal spins that are antialigned and parallel to the orbital angular momentum.

preprint2006arXiv

Event and Apparent Horizon Finders for 3+1 Numerical Relativity

Event and apparent horizons are key diagnostics for the presence and properties of black holes. In this article I review numerical algorithms and codes for finding event and apparent horizons in numerically-computed spacetimes, focusing on calculations done using the 3+1 ADM formalism. There are 3 basic algorithms for finding event horizons, based respectively on integrating null geodesics \emph{forwards} in time, integrating null geodesics \emph{backwards} in time, and integrating null \emph{surfaces} backwards in time. The last of these is generally the most efficient and accurate. There are a large number of apparent-horizon finding algorithms, with differing trade-offs between speed, robustness, accuracy, and ease of programming. In axisymmetry, shooting algorithms work well and are fairly easy to program. In slices with no continuous symmetries, Nakamura et al.'s algorithm and elliptic-PDE algorithms are fast and accurate, but require good initial guesses to converge. In many cases Schnetter's "pretracking" algorithm can greatly improve an elliptic-PDE algorithm's robustness. Flow algorithms are generally quite slow, but can be very robust in their convergence.

preprint1998arXiv

A 3+1 Computational Scheme for Dynamic Spherically Symmetric Black Hole Spacetimes -- I: Initial Data

When using the black hole exclusion (horizon boundary condition) technique, $K$ is usually nonzero and spatially variable, so none of the special cases of York's conformal-decomposition algorithm apply, and the full 4-vector nonlinear York equations must be solved numerically. We discuss the construction of dynamic black hole initial data slices using this technique: We perturb a known black hole slice via some Ansatz, apply the York decomposition (using another Ansatz for the inner boundary conditions) to project the perturbed field variables back into the constraint hypersurface, and finally optionally apply a numerical 3-coordinate transformation to (eg) restore an areal radial coordinate. In comparison to other initial data algorithms, the key advantage of this algorithm is its flexibility: $K$ is unrestricted, allowing the use of whatever slicing is most suitable for (say) a time evolution. We have implemented this algorithm for the spherically symmetric scalar field system. We present numerical results for a number of Eddington- Finkelstein--like initial data slices containing black holes surrounded by scalar field shells. Using 4th order finite differencing with resolutions of $Δr/r \approx 0.02$ (0.01) near the (Gaussian) perturbations, the numerically computed energy and momentum constraints for the final slices are $\ltsim 10^{-8}$ ($10^{-9}$) and $\ltsim 10^{-9}$ ($10^{-10}$) in magnitude. Finally, we briefly discuss the errors incurred when interpolating data from one grid to another, as in numerical coordinate transformation or horizon finding. We show that for the usual moving-local-interpolation schemes, even for smooth functions the interpolation error is not smooth.