Researcher profile

Leor Barack

Leor Barack contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
5topics
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

5 published item(s)

preprint2022arXiv

Time-domain metric reconstruction for hyperbolic scattering

Self-force methods can be applied in calculations of the scatter angle in two-body hyperbolic encounters, working order by order in the mass ratio (assumed small) but with no recourse to a weak-field approximation. This, in turn, can inform ongoing efforts to construct an accurate model of the general-relativistic binary dynamics via an effective-one-body description and other semi-analytical approaches. Existing self-force methods are to a large extent specialised to bound, inspiral orbits. Here we develop a technique for (numerical) self-force calculations that can efficiently tackle scatter orbits. The method is based on a time-domain reconstruction of the metric perturbation from a scalar-like Hertz potential that satisfies the Teukolsky equation, an idea pursued so far only for bound orbits. The crucial ingredient in this formulation are certain jump conditions that (each multipole mode of) the Hertz potential must satisfy along the orbit, in a 1+1-dimensional multipole reduction of the problem. We obtain a closed-form expression for these jumps, for an arbitrary geodesic orbit in Schwarzschild spacetime, and present a full numerical implementation for a scatter orbit. In this paper we focus on method development, and go only as far as calculating the Hertz potential; a calculation of the self-force and its physical effects on the scatter orbit will be the subject of forthcoming work.

preprint2021arXiv

Discontinuous collocation methods and gravitational self-force applications

Numerical simulations of extereme mass ratio inspirals, the mostimportant sources for the LISA detector, face several computational challenges. We present a new approach to evolving partial differential equations occurring in black hole perturbation theory and calculations of the self-force acting on point particles orbiting supermassive black holes. Such equations are distributionally sourced, and standard numerical methods, such as finite-difference or spectral methods, face difficulties associated with approximating discontinuous functions. However, in the self-force problem we typically have access to full a-priori information about the local structure of the discontinuity at the particle. Using this information, we show that high-order accuracy can be recovered by adding to the Lagrange interpolation formula a linear combination of certain jump amplitudes. We construct discontinuous spatial and temporal discretizations by operating on the corrected Lagrange formula. In a method-of-lines framework, this provides a simple and efficient method of solving time-dependent partial differential equations, without loss of accuracy near moving singularities or discontinuities. This method is well-suited for the problem of time-domain reconstruction of the metric perturbation via the Teukolsky or Regge-Wheeler-Zerilli formalisms. Parallel implementations on modern CPU and GPU architectures are discussed.

preprint2010arXiv

Gravitational self force on a particle in circular orbit around a Schwarzschild black hole

We calculate the gravitational self force acting on a pointlike particle of mass $μ$, set in a circular geodesic orbit around a Schwarzschild black hole. Our calculation is done in the Lorenz gauge: For given orbital radius, we first solve directly for the Lorenz-gauge metric perturbation using numerical evolution in the time domain; We then compute the (finite) back-reaction force from each of the multipole modes of the perturbation; Finally, we apply the ``mode sum'' method to obtain the total, physical self force. The {\em temporal} component of the self force (which is gauge invariant) describes the dissipation of orbital energy through gravitational radiation. Our results for this component are consistent, to within the computational accuracy, with the total flux of gravitational-wave energy radiated to infinity and through the event horizon. The {\em radial} component of the self force (which is gauge dependent) is calculated here for the first time. It describes a conservative shift in the orbital parameters away from their geodesic values. We thus obtain the $O(μ)$ correction to the specific energy and angular momentum parameters (in the Lorenz gauge), as well as the $O(μ)$ shift in the orbital frequency (which is gauge invariant).

preprint2009arXiv

Gravitational self force in extreme mass-ratio inspirals

This review is concerned with the gravitational self-force acting on a mass particle in orbit around a large black hole. Renewed interest in this old problem is driven by the prospects of detecting gravitational waves from strongly gravitating binaries with extreme mass ratios. We begin here with a summary of recent advances in the theory of gravitational self-interaction in curved spacetime, and proceed to survey some of the ideas and computational strategies devised for implementing this theory in the case of a particle orbiting a Kerr black hole. We review in detail two of these methods: (i) the standard mode-sum method, in which the metric perturbation is regularized mode-by-mode in a multipole decomposition, and (ii) $m$-mode regularization, whereby individual azimuthal modes of the metric perturbation are regularized in 2+1 dimensions. We discuss several practical issues that arise, including the choice of gauge, the numerical representation of the particle singularity, and how high-frequency contributions near the particle are dealt with in frequency-domain calculations. As an example of a full end-to-end implementation of the mode-sum method, we discuss the computation of the gravitational self-force for eccentric geodesic orbits in Schwarzschild, using a direct integration of the Lorenz-gauge perturbation equations in the time domain. With the computational framework now in place, researchers have recently turned to explore the physical consequences of the gravitational self force; we will describe some preliminary results in this area. An appendix to this review presents, for the first time, a detailed derivation of the regularization parameters necessary for implementing the mode-sum method in Kerr spacetime.

preprint2008arXiv

Two approaches for the gravitational self force in black hole spacetime: Comparison of numerical results

Recently, two independent calculations have been presented of finite-mass ("self-force") effects on the orbit of a point mass around a Schwarzschild black hole. While both computations are based on the standard mode-sum method, they differ in several technical aspects, which makes comparison between their results difficult--but also interesting. Barack and Sago [Phys. Rev. D {\bf 75}, 064021 (2007)] invoke the notion of a self-accelerated motion in a background spacetime, and perform a direct calculation of the local self force in the Lorenz gauge (using numerical evolution of the perturbation equations in the time domain); Detweiler [Phys. Rev. D {\bf 77}, 124026 (2008)] describes the motion in terms a geodesic orbit of a (smooth) perturbed spacetime, and calculates the metric perturbation in the Regge--Wheeler gauge (using frequency-domain numerical analysis). Here we establish a formal correspondence between the two analyses, and demonstrate the consistency of their numerical results. Specifically, we compare the value of the conservative $O(μ)$ shift in $u^t$ (where $μ$ is the particle's mass and $u^t$ is the Schwarzschild $t$ component of the particle's four-velocity), suitably mapped between the two orbital descriptions and adjusted for gauge. We find that the two analyses yield the same value for this shift within mere fractional differences of $\sim 10^{-5}$--$10^{-7}$ (depending on the orbital radius)--comparable with the estimated numerical error.