Source author record

Thomas Osburn

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

ResearcherUnclaimed source record

gr-qc

Catalog footprint

What is connected

4works

1topics

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.

preprint2022arXiv

New self-force method via elliptic partial differential equations for Kerr inspiral models

We present a new method designed to avoid numerical challenges that have impeded calculation of the Lorenz gauge self-force acting on a compact object inspiraling into a Kerr black hole. This type of calculation is valuable in creating waveform templates for extreme mass-ratio inspirals, which are an important source of gravitational waves for the upcoming Laser Interferometer Space Antenna mission. Prior hyperbolic partial differential equation (PDE) formulations encountered numerical instabilities involving unchecked growth in time; our new method is based on elliptic PDEs, which do not exhibit instabilities of that kind. For proof of concept, we calculate the self-force acting on a scalar charge in a circular orbit around a Kerr black hole. We anticipate this method will subsequently facilitate calculation of first-order Lorenz gauge Kerr metric perturbations and self-force, which could serve as a foundation for second-order Kerr self-force investigations.

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.

preprint2015arXiv

Fast spectral source integration in black hole perturbation calculations

This paper presents a new technique for achieving spectral accuracy and fast computational performance in a class of black hole perturbation and gravitational self-force calculations involving extreme mass ratios and generic orbits. Called \emph{spectral source integration} (SSI), this method should see widespread future use in problems that entail (i) point-particle description of the small compact object, (ii) frequency domain decomposition, and (iii) use of the background eccentric geodesic motion. Frequency domain approaches are widely used in both perturbation theory flux-balance calculations and in local gravitational self-force calculations. Recent self-force calculations in Lorenz gauge, using the frequency domain and method of extended homogeneous solutions, have been able to accurately reach eccentricities as high as $e \simeq 0.7$. We show here SSI successfully applied to Lorenz gauge. In a double precision Lorenz gauge code, SSI enhances the accuracy of results and makes a factor of three improvement in the overall speed. The primary initial application of SSI--for us its \emph{raison d'être}--is in an arbitrary precision \emph{Mathematica} code that computes perturbations of eccentric orbits in the Regge-Wheeler gauge to extraordinarily high accuracy (e.g., 200 decimal places). These high accuracy eccentric orbit calculations would not be possible without the exponential convergence of SSI. We believe the method will extend to work for inspirals on Kerr, and will be the subject of a later publication. SSI borrows concepts from discrete-time signal processing and is used to calculate the mode normalization coefficients in perturbation theory via sums over modest numbers of points around an orbit. A variant of the idea is used to obtain spectral accuracy in solution of the geodesic orbital motion.

preprint2014arXiv

Lorenz gauge gravitational self-force calculations of eccentric binaries using a frequency domain procedure

We present an algorithm for calculating the metric perturbations and gravitational self-force for extreme-mass-ratio inspirals (EMRIs) with eccentric orbits. The massive black hole is taken to be Schwarzschild and metric perturbations are computed in Lorenz gauge. The perturbation equations are solved as coupled systems of ordinary differential equations in the frequency domain. Accurate local behavior of the metric is attained through use of the method of extended homogeneous solutions and mode-sum regularization is used to find the self-force. We focus on calculating the self-force with sufficient accuracy to ensure its error contributions to the phase in a long term orbital evolution will be $δΦ\lesssim 10^{-2}$ radians. This requires the orbit-averaged force to have fractional errors $\lesssim 10^{-8}$ and the oscillatory part of the self-force to have errors $\lesssim 10^{-3}$ (a level frequently easily exceeded). Our code meets this error requirement in the oscillatory part, extending the reach to EMRIs with eccentricities of $e \lesssim 0.8$, if augmented by use of fluxes for the orbit-averaged force, or to eccentricities of $e \lesssim 0.5$ when used as a stand-alone code. Further, we demonstrate accurate calculations up to orbital separations of $a \simeq 100 M$, beyond that required for EMRI models and useful for comparison with post-Newtonian theory. Our principal developments include (1) use of fully constrained field equations, (2) discovery of analytic solutions for even-parity static modes, (3) finding a pre-conditioning technique for outer homogeneous solutions, (4) adaptive use of quad-precision and (5) jump conditions to handle near-static modes, and (6) a hybrid scheme for high eccentricities.