Source author record

Jonathan Thornburg

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

ResearcherUnclaimed source record

gr-qc astro-ph.CO astro-ph.GA astro-ph.HE physics.comp-ph

Catalog footprint

What is connected

10works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

preprint2012arXiv

Generic effective source for scalar self-force calculations

A leading approach to the modelling of extreme mass ratio inspirals involves the treatment of the smaller mass as a point particle and the computation of a regularized self-force acting on that particle. In turn, this computation requires knowledge of the regularized retarded field generated by the particle. A direct calculation of this regularized field may be achieved by replacing the point particle with an effective source and solving directly a wave equation for the regularized field. This has the advantage that all quantities are finite and require no further regularization. In this work, we present a method for computing an effective source which is finite and continuous everywhere, and which is valid for a scalar point particle in arbitrary geodesic motion in an arbitrary background spacetime. We explain in detail various technical and practical considerations that underlie its use in several numerical self-force calculations. We consider as examples the cases of a particle in a circular orbit about Schwarzschild and Kerr black holes, and also the case of a particle following a generic time-like geodesic about a highly spinning Kerr black hole. We provide numerical C code for computing an effective source for various orbital configurations about Schwarzschild and Kerr black holes.

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

Adaptive Mesh Refinement for Characteristic Grids

I consider techniques for Berger-Oliger adaptive mesh refinement (AMR) when numerically solving partial differential equations with wave-like solutions, using characteristic (double-null) grids. Such AMR algorithms are naturally recursive, and the best-known past Berger-Oliger characteristic AMR algorithm, that of Pretorius & Lehner (J. Comp. Phys. 198 (2004), 10), recurses on individual "diamond" characteristic grid cells. This leads to the use of fine-grained memory management, with individual grid cells kept in 2-dimensional linked lists at each refinement level. This complicates the implementation and adds overhead in both space and time. Here I describe a Berger-Oliger characteristic AMR algorithm which instead recurses on null \emph{slices}. This algorithm is very similar to the usual Cauchy Berger-Oliger algorithm, and uses relatively coarse-grained memory management, allowing entire null slices to be stored in contiguous arrays in memory. The algorithm is very efficient in both space and time. I describe discretizations yielding both 2nd and 4th order global accuracy. My code implementing the algorithm described here is included in the electronic supplementary materials accompanying this paper, and is freely available to other researchers under the terms of the GNU general public license.

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.

preprint2000arXiv

Type II Critical Collapse of a Self-Gravitating Nonlinear $σ$-Model

We report on the existence and phenomenology of type II critical collapse within the one-parameter family of SU(2) $σ$-models coupled to gravity. Numerical investigations in spherical symmetry show discretely self-similar (DSS) behavior at the threshold of black hole formation for values of the dimensionless coupling constant $\ccbeta$ ranging from 0.2 to 100; at 0.18 we see small deviations from DSS. While the echoing period $Δ$ of the critical solution rises sharply towards the lower limit of this range, the characteristic mass scaling has a critical exponent $γ$ which is almost independent of $\ccbeta$, asymptoting to $0.1185 \pm 0.0005$ at large $\ccbeta$. We also find critical scaling of the scalar curvature for near-critical initial data. Our numerical results are based on an outgoing-null-cone formulation of the Einstein-matter equations, specialized to spherical symmetry. Our numerically computed initial-data critical parameters $p^*$ show 2nd order convergence with the grid resolution, and after compensating for this variation in $p^*$, our individual evolutions are uniformly 2nd order convergent even very close to criticality.

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.