Researcher profile

Sarp Akcay

Sarp Akcay contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
11works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

11 published item(s)

preprint2021arXiv

A hybrid post-Newtonian -- effective-one-body scheme for spin-precessing compact-binary waveforms

We introduce \texttt{TEOBResumSP}: an efficient, accurate hybrid scheme for generating gravitational waveforms from spin-precessing compact binaries. The precessing waveforms are generated via the established technique of Euler rotating the non-precessing \texttt{TEOBResumS} waveforms from a precessing frame to an inertial frame. We obtain the Euler angles by solving the post-Newtonian precession equations expanded to second post-Newtonian order. Current version of \texttt{TEOBResumSP} produces precessing waveforms through the inspiral phase up to the onset of the merger. We compare \texttt{TEOBResumSP} to current state-of-the-art precessing approximants \texttt{NRSur7dq4}, \texttt{SEOBNRv4PHM}, and \texttt{IMRPhenomPv3HM} for 200 cases of precessing compact binary inspirals with orbital inclinations up to 90 degrees, mass ratios up to four, and the effective precession parameter $χ_p$ up to 0.75. We further provide an extended comparison with \texttt{SEOBNRv4PHM} involving 1030 more inspirals with $χ_p\le 1$ and mass ratios up to 10. We find that 91\% of the \texttt{TEOBResumSP}-\texttt{NRSur7dq4} matches, 85\% of the \texttt{TEOBResumSP}-\texttt{SEOBNRv4PHM} matches, and 77\% of the \texttt{TEOBResumSP}-\texttt{IMRPhenomPv3HM} matches are greater than $0.965$. Most disagreements occur for large mass ratios and $χ_p \gtrsim 0.6$. We identify the mismatch of the \emph{non}-precessing $(2,1)$ mode as one of the leading causes of disagreements. We also introduce a new parameter, $χ_{\perp,\text{max}}$, to measure the strength of precession and hint that the mismatch between the above approximants shows an exponential dependence on $χ_{\perp,\text{max}}$ though this requires further study. Our results indicate that \texttt{TEOBResumSP} is on its way to becoming a robust precessing approximant to be employed in the parameter estimation of generic-spin compact binaries.

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.

preprint2016arXiv

Numerical computation of the EOB potential q using self-force results

The effective-one-body theory (EOB) describes the conservative dynamics of compact binary systems in terms of an effective Hamiltonian approach. The Hamiltonian for moderately eccentric motion of two non-spinning compact objects in the extreme mass-ratio limit is given in terms of three potentials: $a(v), \bar{d}(v), q(v)$. By generalizing the first law of mechanics for (non-spinning) black hole binaries to eccentric orbits, [\prd{\bf92}, 084021 (2015)] recently obtained new expressions for $\bar{d}(v)$ and $q(v)$ in terms of quantities that can be readily computed using the gravitational self-force approach. Using these expressions we present a new computation of the EOB potential $q(v)$ by combining results from two independent numerical self-force codes. We determine $q(v)$ for inverse binary separations in the range $1/1200 \le v \lesssim 1/6$. Our computation thus provides the first-ever strong-field results for $q(v)$. We also obtain $\bar{d}(v)$ in our entire domain to a fractional accuracy of $\gtrsim 10^{-8}$. We find to our results are compatible with the known post-Newtonian expansions for $\bar{d}(v)$ and $q(v)$ in the weak field, and agree with previous (less accurate) numerical results for $\bar{d}(v)$ in the strong field.

preprint2015arXiv

Comparison Between Self-Force and Post-Newtonian Dynamics: Beyond Circular Orbits

The gravitational self-force (GSF) and post-Newtonian (PN) schemes are complementary approximation methods for modelling the dynamics of compact binary systems. Comparison of their results in an overlapping domain of validity provides a crucial test for both methods, and can be used to enhance their accuracy, e.g. via the determination of previously unknown PN parameters. Here, for the first time, we extend such comparisons to noncircular orbits---specifically, to a system of two nonspinning objects in a bound (eccentric) orbit. To enable the comparison we use a certain orbital-averaged quantity $\langle U \rangle$ that generalizes Detweiler's redshift invariant. The functional relationship $\langle U \rangle(Ω_r,Ω_ϕ)$, where $Ω_r$ and $Ω_ϕ$ are the frequencies of the radial and azimuthal motions, is an invariant characteristic of the conservative dynamics. We compute $\langle U \rangle(Ω_r,Ω_ϕ)$ numerically through linear order in the mass ratio $q$, using a GSF code which is based on a frequency-domain treatment of the linearized Einstein equations in the Lorenz gauge. We also derive $\langle U \rangle(Ω_r,Ω_ϕ)$ analytically through 3PN order, for an arbitrary $q$, using the known near-zone 3PN metric and the generalized quasi-Keplerian representation of the motion. We demonstrate that the $\mathcal{O}(q)$ piece of the analytical PN prediction is perfectly consistent with the numerical GSF results, and we use the latter to estimate yet unknown pieces of the 4PN expression at $\mathcal{O}(q)$.

preprint2013arXiv

Frequency-domain algorithm for the Lorenz-gauge gravitational self-force

State-of-the-art computations of the gravitational self-force (GSF) on massive particles in black hole spacetimes involve numerical evolution of the metric perturbation equations in the time-domain, which is computationally very costly. We present here a new strategy, based on a frequency-domain treatment of the perturbation equations, which offers considerable computational saving. The essential ingredients of our method are (i) a Fourier-harmonic decomposition of the Lorenz-gauge metric perturbation equations and a numerical solution of the resulting coupled set of ordinary equations with suitable boundary conditions; (ii) a generalized version of the method of extended homogeneous solutions [Phys. Rev. D {\bf 78}, 084021 (2008)] used to circumvent the Gibbs phenomenon that would otherwise hamper the convergence of the Fourier mode-sum at the particle's location; and (iii) standard mode-sum regularization, which finally yields the physical GSF as a sum over regularized modal contributions. We present a working code that implements this strategy to calculate the Lorenz-gauge GSF along eccentric geodesic orbits around a Schwarzschild black hole. The code is far more efficient than existing time-domain methods; the gain in computation speed (at a given precision) is about an order of magnitude at an eccentricity of 0.2, and up to three orders of magnitude for circular or nearly circular orbits. This increased efficiency was crucial in enabling the recently reported calculation of the long-term orbital evolution of an extreme mass ratio inspiral [Phys. Rev. D {\bf 85}, 061501(R) (2012)]. Here we provide full technical details of our method to complement the above report.

preprint2012arXiv

A Fast Frequency-Domain Algorithm for Gravitational Self-Force: I, Circular Orbits in Schwarzschild Spacetime

Fast, reliable orbital evolutions of compact objects around massive black holes will be needed as input for gravitational wave search algorithms in the data stream generated by the planned Laser Interferometer Space Antenna (LISA). Currently, the state of the art is a time-domain code by [Phys. Rev. D{\bf 81}, 084021, (2010)] that computes the gravitational self-force on a point-particle in an eccentric orbit around a Schwarzschild black hole. Currently, time-domain codes take up to a few days to compute just one point in parameter space. In a series of articles, we advocate the use of a frequency-domain approach to the problem of gravitational self-force (GSF) with the ultimate goal of orbital evolution in mind. Here, we compute the GSF for a particle in a circular orbit in Schwarzschild spacetime. We solve the linearized Einstein equations for the metric perturbation in Lorenz gauge. Our frequency-domain code reproduces the time-domain results for the GSF up to $\sim 1000$ times faster for small orbital radii. In forthcoming companion papers, we will generalize our frequency-domain methods to include bound (eccentric) orbits in Schwarzschild and (eventually) Kerr spacetimes for computing the GSF, where we will employ the method of extended homogeneous solutions [Phys. Rev. D {\bf 78}, 084021 (2008)].

preprint2012arXiv

Critical phenomena at the threshold of immediate merger in binary black hole systems: the extreme mass ratio case

In numerical simulations of black hole binaries, Pretorius and Khurana [Class. Quant. Grav. {\bf 24}, S83 (2007)] have observed critical behaviour at the threshold between scattering and immediate merger. The number of orbits scales as $n\simeq -γ\ln|p-p_*|$ along any one-parameter family of initial data such that the threshold is at $p=p_*$. Hence they conjecture that in ultrarelavistic collisions almost all the kinetic energy can be converted into gravitational waves if the impact parameter is fine-tuned to the threshold. As a toy model for the binary, they consider the geodesic motion of a test particle in a Kerr black hole spacetime, where the unstable circular geodesics play the role of critical solutions, and calculate the critical exponent $γ$. Here, we incorporate radiation reaction into this model using the self-force approximation. The critical solution now evolves adiabatically along a sequence of unstable circular geodesic orbits under the effect of the self-force. We confirm that almost all the initial energy and angular momentum are radiated on the critical solution. Our calculation suggests that, even for infinite initial energy, this happens over a finite number of orbits given by $n_\infty\simeq 0.41/η$, where $η$ is the (small) mass ratio. We derive expressions for the time spent on the critical solution, number of orbits and radiated energy as functions of the initial energy and impact parameter.

preprint2012arXiv

Gravitational self-force and the effective-one-body formalism between the innermost stable circular orbit and the light ring

We compute the conservative piece of the gravitational self-force (GSF) acting on a particle of mass m_1 as it moves along an (unstable) circular geodesic orbit between the innermost stable circular orbit (ISCO) and the light ring of a Schwarzschild black hole of mass m_2>> m_1. More precisely, we construct the function h_{uu}(x) = h_{μν} u^μ u^ν (related to Detweiler&#39;s gauge-invariant &#34;redshift&#34; variable), where h_{μν} is the regularized metric perturbation in the Lorenz gauge, u^μ is the four-velocity of m_1, and x= [Gc^{-3}(m_1+m_2)Ω]^{2/3} is an invariant coordinate constructed from the orbital frequency Ω. In particular, we explore the behavior of h_{uu} just outside the &#34;light ring&#34; at x=1/3, where the circular orbit becomes null. Using the recently discovered link between h_{uu} and the piece a(u), linear in the symmetric mass ratio ν, of the main radial potential A(u,ν) of the Effective One Body (EOB) formalism, we compute a(u) over the entire domain 0<u<1/3. We find that a(u) diverges at the light-ring as ~0.25 (1-3u)^{-1/2}, explain the physical origin of this divergence, and discuss its consequences for the EOB formalism. We construct accurate global analytic fits for a(u), valid on the entire domain 0<u<1/3 (and possibly beyond), and give accurate numerical estimates of the values of a(u) and its first 3 derivatives at the ISCO, as well as the O(ν) shift in the ISCO frequency. In previous work we used GSF data on slightly eccentric orbits to compute a certain linear combination of a(u) and its first two derivatives, involving also the O(ν) piece \bar d(u) of a second EOB radial potential {\bar D}(u,ν). Combining these results with our present global analytic representation of a(u), we numerically compute {\bar d}(u)$ on the interval 0<u\leq 1/6.

preprint2011arXiv

Horizon-absorbed energy flux in circularized, nonspinning black-hole binaries and its effective-one-body representation

We propose, within the effective one body (EOB) approach, a new, resummed, analytical representation of the gravitational wave energy flux absorbed by a system of two circularized (nonspinning) black holes. This expression is such to be well-behaved in the strong-field, fast motion regime, notably up to the EOB-defined last unstable orbit. Building conceptually upon the procedure adopted to resum the multipolar asymptotic energy flux, we introduce a {\it multiplicative} decomposition of the multipolar absorbed flux made by three factors: (i) the leading-order contribution, (ii) an &#34;effective source&#34; and (iii) a new residual amplitude correction $(\tildeρ_\lm^H)^{2\ell}$. In the test-mass limit, we use a frequency-domain perturbative approach to accurately compute numerically the horizon-absorbed fluxes along a sequence of stable and unstable circular orbits and we extract from them the functions $\tildeρ_\lm^H$. These quantities are then fitted via rational functions. The resulting analytically represented test-mass knowledge is then suitably {\it hybridized} with lower-order analytical information that is valid for any mass ratio. This yields a resummed representation of the absorbed flux for a generic, circularized, nonspinning black-hole binary. Our result adds new information to the state-of-the-art calculation of the absorbed flux at fractional 5 post-Newtonian order [S. Taylor and E. Poisson, Phys. Rev. D {\bf 78} 084016 (2008)], that is recovered in the weak-field limit approximation by construction.

preprint2010arXiv

Kerr-de Sitter Universe

It is now widely accepted that the universe as we understand it is accelerating in expansion and fits the de Sitter model rather well. As such, a realistic assumption of black holes must place them on a de Sitter background and not Minkowski as is typically done in General Relativity. The most astrophysically relevant black hole is the uncharged, rotating Kerr solution, a member of the more general Kerr-Newman metrics. A generalization of the rotating Kerr black hole to a solution of the Einstein&#39;s equation with a cosmological constant $Λ$ was discovered by Carter \cite{DWDW}. It is typically referred to as the Kerr-de Sitter spacetime. Here, we discuss the horizon structure of this spacetime and its dependence on $Λ$. We recall that in a $\La>0$ universe, the term `extremal black hole&#39; refers to a black hole with angular momentum $J > M^2 $. We obtain explicit numerical results for the black hole&#39;s maximal spin value and get a distribution of admissible Kerr holes in the ($Λ$, spin) parameter space. We look at the conformal structure of the extended spacetime and the embedding of the 3-geometry of the spatial hypersurfaces. In analogy with Reissner-Nordström -de Sitter spacetime, in particular by considering the Kerr-de Sitter causal structure as a distortion of the Reissner-Nordström-de Sitter one, we show that spatial sections of the extended spacetime are 3-spheres containing 2-dimensional topologically spherical sections of the horizons of Kerr holes at the poles. Depending on how a $t=$ constant 3-space is defined these holes may be seen as black or white holes (four possible combinations).

preprint2009arXiv

Area Invariance of Apparent Horizons under Arbitrary Boosts

It is a well known analytic result in general relativity that the 2-dimensional area of the apparent horizon of a black hole remains invariant regardless of the motion of the observer, and in fact is independent of the $ t=constant $ slice, which can be quite arbitrary in general relativity. Nonetheless the explicit computation of horizon area is often substantially more difficult in some frames (complicated by the coordinate form of the metric), than in other frames. Here we give an explicit demonstration for very restricted metric forms of (Schwarzschild and Kerr) vacuum black holes. In the Kerr-Schild coordinate expression for these spacetimes they have an explicit Lorentz-invariant form. We consider {\it boosted} versions with the black hole moving through the coordinate system. Since these are stationary black hole spacetimes, the apparent horizons are two dimensional cross sections of their event horizons, so we compute the areas of apparent horizons in the boosted space with (boosted) $ t = constant $, and obtain the same result as in the unboosted case. Note that while the invariance of area is generic, we deal only with black holes in the Kerr-Schild form, and consider only one particularly simple change of slicing which amounts to a boost. Even with these restrictions we find that the results illuminate the physics of the horizon as a null surface and provide a useful pedagogical tool. As far as we can determine, this is the first explicit calculation of this type demonstrating the area invariance of horizons. Further, these calculations are directly relevant to transformations that arise in computational representation of moving black holes. We present an application of this result to initial data for boosted black holes.