Source author record

Barry Wardell

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

27 published item(s)

preprint2022arXiv

Gravitational perturbations of rotating black holes in Lorenz gauge

Perturbations of Kerr spacetime are typically studied with the Teukolsky formalism, in which a pair of invariant components of the perturbed Weyl tensor are expressed in terms of separable modes that satisfy ordinary differential equations. However, for certain applications it is desirable to construct the full metric perturbation in the Lorenz gauge, in which the linearized Einstein field equations take a manifestly hyperbolic form. Here we obtain a set of Lorenz-gauge solutions to the vacuum field equations in terms of homogeneous solutions to the spin-2, spin-1 and spin-0 Teukolsky equations; and completion pieces that represent perturbations to the mass and angular momentum of the spacetime. The solutions are valid in vacuum Petrov type-D spacetimes that admit a conformal Killing-Yano tensor.

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.

preprint2021arXiv

Self-Force Calculations with a Spinning Secondary

We compute the linear metric perturbation to a Schwarzschild black hole generated by a spinning compact object, specialising to circular equatorial orbits with an (anti-)aligned spin vector. We derive a two-timescale expansion of the field equations, with an attendant waveform-generation framework, that includes all effects through first post-adiabatic order, and we use the Regge-Wheeler-Zerilli formalism in the frequency domain to generate waveforms that include the complete effect of the spin on the waveform phase. We perform the calculations using expansions at fixed orbital frequency, increasing the computational efficiency and simplifying the procedure compared to previous approaches. Finally, we provide the first fully relativistic, first-principles regularisation procedure for gauge invariant self-force quantities to linear order in spin. We use this procedure to produce the first strong-field, conservative self-force calculation including the spin of the secondary -- computing Detweiler's redshift invariant.

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.

preprint2020arXiv

Transport Equation Approach to Calculations of Hadamard Green functions and non-coincident DeWitt coefficients

Building on an insight due to Avramidi, we provide a system of transport equations for determining key fundamental bi-tensors, including derivatives of the world-function, σ(x,x'), the square root of the Van Vleck determinant, Δ^{1/2}(x,x'), and the tail-term, V(x,x'), appearing in the Hadamard form of the Green function. These bi-tensors are central to a broad range of problems from radiation reaction to quantum field theory in curved spacetime and quantum gravity. Their transport equations may be used either in a semi-recursive approach to determining their covariant Taylor series expansions, or as the basis of numerical calculations. To illustrate the power of the semi-recursive approach, we present an implementation in \textsl{Mathematica} which computes very high order covariant series expansions of these objects. Using this code, a moderate laptop can, for example, calculate the coincidence limit a_7(x,x) and V(x,x') to order (σ^a)^{20} in a matter of minutes. Results may be output in either a compact notation or in xTensor form. In a second application of the approach, we present a scheme for numerically integrating the transport equations as a system of coupled ordinary differential equations. As an example application of the scheme, we integrate along null geodesics to solve for V(x,x') in Nariai and Schwarzschild spacetimes.

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.

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.

preprint2016arXiv

Science with the space-based interferometer eLISA. I: Supermassive black hole binaries

We compare the science capabilities of different eLISA mission designs, including four-link (two-arm) and six-link (three-arm) configurations with different arm lengths, low-frequency noise sensitivities and mission durations. For each of these configurations we consider a few representative massive black hole formation scenarios. These scenarios are chosen to explore two physical mechanisms that greatly affect eLISA rates, namely (i) black hole seeding, and (ii) the delays between the merger of two galaxies and the merger of the black holes hosted by those galaxies. We assess the eLISA parameter estimation accuracy using a Fisher matrix analysis with spin-precessing, inspiral-only waveforms. We quantify the information present in the merger and ringdown by rescaling the inspiral-only Fisher matrix estimates using the signal-to-noise ratio from non-precessing inspiral-merger-ringdown phenomenological waveforms, and from a reduced set of precessing numerical relativity/post-Newtonian hybrid waveforms. We find that all of the eLISA configurations considered in our study should detect some massive black hole binaries. However, configurations with six links and better low-frequency noise will provide much more information on the origin of black holes at high redshifts and on their accretion history, and they may allow the identification of electromagnetic counterparts to massive black hole mergers.

preprint2016arXiv

Second-order perturbation theory: the problem of infinite mode coupling

Second-order self-force computations, which will be essential in modeling extreme-mass-ratio inspirals, involve two major new difficulties that were not present at first order. One is the problem of large scales, discussed in [Phys. Rev. D 92, 104047 (2015)]. Here we discuss the second difficulty, which occurs instead on small scales: if we expand the field equations in spherical harmonics, then because the first-order field contains a singularity, we require an arbitrarily large number of first-order modes to accurately compute even a single second-order mode. This is a generic feature of nonlinear field equations containing singularities, allowing us to study it in the simple context of a scalar toy model in flat space. Using that model, we illustrate the problem and demonstrate a robust strategy for overcoming it.

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.

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.

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.

preprint2014arXiv

High-order expansions of the Detweiler-Whiting singular field in Kerr spacetime

In a previous paper, we computed expressions for the Detweiler-Whiting singular field of point scalar, electromagnetic and gravitational charges following a geodesic of the Schwarzschild spacetime. We now extend this to the case of equatorial orbits in Kerr spacetime, using coordinate and covariant approaches to compute expansions of the singular field in scalar, electromagnetic and gravitational cases. As an application, we give the calculation of previously unknown mode-sum regularization parameters. We also propose a new application of high-order approximations to the singular field, showing how they may be used to compute $m$-mode regularization parameters for use in the $m$-mode effective source approach to self-force calculations.

preprint2014arXiv

Self-force via Green functions and worldline integration

A compact object moving in curved spacetime interacts with its own gravitational field. This leads to both dissipative and conservative corrections to the motion, which can be interpreted as a self-force acting on the object. The original formalism describing this self-force relied heavily on the Green function of the linear differential operator that governs gravitational perturbations. However, because the global calculation of Green functions in non-trivial black hole spacetimes has been an open problem until recently, alternative methods were established to calculate self-force effects using sophisticated regularization techniques that avoid the computation of the global Green function. We present a method for calculating the self-force that employs the global Green function and is therefore closely modeled after the original self-force expressions. Our quantitative method involves two stages: (i) numerical approximation of the retarded Green function in the background spacetime; (ii) evaluation of convolution integrals along the worldline of the object. This novel approach can be used along arbitrary worldlines, including those currently inaccessible to more established computational techniques. Furthermore, it yields geometrical insight into the contributions to self-interaction from curved geometry (back-scattering) and trapping of null geodesics. We demonstrate the method on the motion of a scalar charge in Schwarzschild spacetime. This toy model retains the physical history-dependence of the self-force but avoids gauge issues and allows us to focus on basic principles. We compute the self-field and self-force for many worldlines including accelerated circular orbits, eccentric orbits at the separatrix, and radial infall. This method, closely modeled after the original formalism, provides a promising complementary approach to the self-force problem.

preprint2013arXiv

Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration

The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.

preprint2013arXiv

Falloff of the Weyl scalars in binary black hole spacetimes

The peeling theorem of general relativity predicts that the Weyl curvature scalars Psi_n (n=0...4), when constructed from a suitable null tetrad in an asymptotically flat spacetime, fall off asymptotically as r^(n-5) along outgoing radial null geodesics. This leads to the interpretation of Psi_4 as outgoing gravitational radiation at large distances from the source. We have performed numerical simulations in full general relativity of a binary black hole inspiral and merger, and have computed the Weyl scalars in the standard tetrad used in numerical relativity. In contrast with previous results, we observe that all the Weyl scalars fall off according to the predictions of the theorem.

preprint2013arXiv

High-order expansions of the Detweiler-Whiting singular field in Schwarzschild spacetime

The self field of a charged particle has a component that diverges at the particle. We use both coordinate and covariant approaches to compute an expansion of this singular field for generic geodesic orbits in Schwarzschild spacetime for scalar, electromagnetic and graviational cases. We check agreement of both approaches and give, as an application, the calculation of previously unknown regularisation parameters. In this so-called "mode-sum regularization" approach, each mode of the field is finite, while their sum diverges. The sum may be rendered finite and convergent by the subtraction of "regularization parameters". Higher order parameters lead to faster convergence in the mode-sum. As a second example application, we compute high order expressions for the effective source approach to self-force calculations.

preprint2013arXiv

Scalar self-force for eccentric orbits around a Schwarzschild black hole

We revisit the problem of computing the self-force on a scalar charge moving along an eccentric geodesic orbit around a Schwarzschild black hole. This work extends previous scalar self-force calculations for circular orbits, which were based on a regular "effective" point-particle source and a full 3D evolution code. We find good agreement between our results and previous calculations based on a (1+1) time-domain code. Finally, our data visualization is unconventional: we plot the self-force through full radial cycles to create "self-force loops", which reveal many interesting features that are less apparent in standard presentations of eccentric-orbit self-force data.

preprint2013arXiv

Self-Force and Green Function in Schwarzschild spacetime via Quasinormal Modes and Branch Cut

The motion of a small compact object in a curved background spacetime deviates from a geodesic due to the action of its own field, giving rise to a self-force. This self-force may be calculated by integrating the Green function for the wave equation over the past worldline of the small object. We compute the self-force in this way for the case of a scalar charge in Schwarzschild spacetime, making use of the semi-analytic method of matched Green function expansions. Inside a local neighbourhood of the compact object, this method uses the Hadamard form for the Green function in order to render regularization trivial. Outside this local neighbourhood, we calculate the Green function using a spectral decomposition into poles (quasinormal modes) and a branch cut integral in the complex-frequency plane. We show that both expansions overlap in a sufficiently large matching region for an accurate calculation of the self-force to be possible. The scalar case studied here is a useful and illustrative toy-model for the gravitational case, which serves to model astrophysical binary systems in the extreme mass-ratio limit.

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.

preprint2012arXiv

Green Functions and Radiation Reaction From a Spacetime Perspective

The basis of this work is the first full application of the Poisson-Wiseman-Anderson method of `matched expansions' to compute the self-force acting on a point particle moving in a curved spacetime. The method employs two expansions for the Green function which are respectively valid in the `quasilocal' and `distant past' regimes, and which are matched together within the normal neighborhood. Building on a fundamental insight due to Avramidi, we provide a system of transport equations for determining key fundamental bi-tensors, including the tail-term, V(x,x'), appearing in the Hadamard form of the Green function. These bitensors are central to a broad range of problems from radiation reaction and the self-force to quantum field theory in curved spacetime and quantum gravity. Using their transport equations, we show how the quasilocal Green function may be computed throughout the normal neighborhood both numerically and as a covariant Taylor series expansion. These calculations are carried out for several black hole spacetimes. Finally, we present a complete application of the method of matched expansions. The calculation is performed in a static region of the spherically symmetric Nariai spacetime (dS_2 X S^2), where the matched expansion method is applied to compute the scalar self-force acting on a static particle. We find that the matched expansion method provides insight into the non-local properties of the self-force. The Green function in Schwarzschild spacetime is expected to share certain key features with Nariai. In this way, the Nariai spacetime provides a fertile testing ground for developing insight into the non-local part of the self-force on black hole spacetimes.

preprint2012arXiv

Quasi-local contribution to the scalar self-force: Non-geodesic Motion

We extend our previous calculation of the quasi-local contribution to the self-force on a scalar particle to general (not necessarily geodesic) motion in a general spacetime. In addition to the general case and the case of a particle at rest in a stationary spacetime, we consider as examples a particle held at rest in Reissner-Nordstrom and Kerr-Newman space-times. This allows us to most easily analyse the effect of non-geodesic motion on our previous results and also allows for comparison to existing results for Schwarzschild spacetime.

preprint2012arXiv

Self-consistent orbital evolution of a particle around a Schwarzschild black hole

The motion of a charged particle is influenced by the self-force arising from the particle's interaction with its own field. In a curved spacetime, this self-force depends on the entire past history of the particle and is difficult to evaluate. As a result, all existing self-force evaluations in curved spacetime are for particles moving along a fixed trajectory. Here, for the first time, we overcome this long-standing limitation and present fully self-consistent orbits and waveforms of a scalar charged particle around a Schwarzschild black hole.

preprint2011arXiv

Effective source approach to self-force calculations

Numerical evaluation of the self-force on a point particle is made difficult by the use of delta functions as sources. Recent methods for self-force calculations avoid delta functions altogether, using instead a finite and extended "effective source" for a point particle. We provide a review of the general principles underlying this strategy, using the specific example of a scalar point charge moving in a black hole spacetime. We also report on two new developments: (i) the construction and evaluation of an effective source for a scalar charge moving along a generic orbit of an arbitrary spacetime, and (ii) the successful implementation of hyperboloidal slicing that significantly improves on previous treatments of boundary conditions used for effective-source-based self-force calculations. Finally, we identify some of the key issues related to the effective source approach that will need to be addressed by future work.

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).

preprint2010arXiv

A Transport Equation Approach to Green Functions and Self-force Calculations

In a recent work, we presented the first application of the Poisson-Wiseman-Anderson method of `matched expansions' to compute the self-force acting on a point particle moving in a curved spacetime. The method employs two expansions for the Green function which are respectively valid in the `quasilocal' and `distant past' regimes, and which may be matched together within the normal neighbourhood. In this article, we introduce the method of matched expansions and discuss transport equation methods for the calculation of the Green function in the quasilocal region. These methods allow the Green function to be evaluated throughout the normal neighborhood and are also relevant to a broad range of problems from radiation reaction to quantum field theory in curved spacetime and quantum gravity.

preprint2010arXiv

Method of Matched Expansions & the Singularity Structure of the Green Function

We present the first successful application of the method of Matched Expansions for the calculation of the self-force on a point particle in a curved spacetime. We investigate the case of a scalar charge in the Nariai spacetime, which serves as a toy model for a point mass moving in the Schwarzschild black hole background. We discuss the singularity structure of the Green function beyond the normal neighbourhood and the interesting effect of caustics on null wave propagation.