Source author record

Alessandro Nagar

Alessandro Nagar 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

40works
3topics
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

40 published item(s)

preprint2022arXiv

Assessment of Effective-One-Body Radiation Reactions for Generic Planar Orbits

In this paper we assess the performances of different analytical prescriptions for the effective-one-body (EOB) radiation reaction along generic planar orbits using exact numerical result in the test-mass limit. We consider three prescriptions put forward in the recent literature: (i) the quasicircular prescription (QC), (ii) the QC with second post-Newtonian (2PN) order noncircular corrections (QC2PN), and (iii) the QC corrected by the noncircular Newtonian prefactor (NCN). The analytical fluxes are then compared against the exact fluxes that are computed by solving the Teukolsky equation with a test-mass source in geodesic motion. We find that the NCN prescription is the most accurate for both eccentric and hyperbolic orbits and it is in robust agreement also for large values of the eccentricity. This result carries over to the comparable masses, as we discuss for a numerical-relativity (NR) case study. We also demonstrate that, while the EOB/NR waveform unfaithfulness is a necessary check for the precision of EOB models, the direct comparison of EOB/NR fluxes is a more stringent and informative test to select the best prescription. Finally, we propose an improved radiation reaction, NCN2PN, that includes noncircular 2PN corrections, in resummed form, as a further multiplicative contribution and that is valid for any mass ratio.

preprint2022arXiv

Effective one-body model for extreme-mass-ratio spinning binaries on eccentric equatorial orbits: testing radiation reaction and waveform

We provide a systematic analysis of the multipolar gravitational waveform, energy and angular momentum fluxes emitted by a nonspinning test particle orbiting a Kerr black hole along equatorial, eccentric orbits. These quantities are computed by numerically solving the Teukolsky equation in the time domain and are then used to test and improve the radiation reaction (and waveform) of an effective-one-body (EOB) model. Eccentricity is incorporated into EOB by replacing the quasi-circular Newtonian (or leading-order) prefactors in the EOB-factorized multipolar waveform (and fluxes) with their generic counterparts. The comparison between numerical and analytical quantities is carried out over a large portion of the parameter space, notably for orbits close to the separatrix and with high eccentricities. The analytical model agrees to $\sim 1\%$ with the numerical data for orbits with moderate eccentricities ($e\lesssim 0.3$) and moderate spins ($\hat{a}\lesssim 0.5$), although this increases up to $\sim 33\%$ for large, positive, black hole spins ($\sim 0.9$) and large eccentricities ($\sim 0.9$). For moderate eccentricities, the new EOB fluxes can be used to drive the dynamics through the nonadiabatic transition from eccentric inspiral to plunge, merger and ringdown, thus providing accurate an description of the merger. We also show that the same technique can be applied to hyperbolic encounters. In general, our approach to radiation reaction for eccentric inspirals should be seen as a first step toward EOB modelization of extreme-mass-ratio-inspirals waveforms for LISA.

preprint2022arXiv

Exploiting Newton-factorized, 2PN-accurate, waveform multipoles in effective-one-body models for spin-aligned noncircularized binaries

We present a new approach to factorize and resum the post-Newtonian (PN) waveform for generic equatorial motion to be used within effective-one-body (EOB) based waveform models. The new multipolar waveform factorization improves previous prescriptions in that: (i) the generic Newtonian contribution is factored out from each multipole; (ii) the circular part is factored out and resummed using standard EOB methods and (iii) the residual, 2PN-accurate, noncircular part, and in particular the tail contribution, is additionally resummed using Padé approximants. The resulting waveform is validated in the extreme-mass-ratio limit by comparisons with nine (mostly nonspinning) numerical waveforms either from eccentric inspirals, with eccentricities up to $e=0.9$, or dynamical captures . The resummation of the noncircular tail contribution is found essential to obtain excellent (${\lesssim}0.05$~rad at periastron for $e=0.9$) analytical/numerical agreement and to considerably improve the prescription with just the Newtonian prefactor. In the comparable mass case, the new 2PN waveform shows only a marginal improvement over the previous Newtonian factorization, though yielding maximal unfaithfulness $\simeq 10^{-3}$ with the 28 publicly available numerical relativity simulations with eccentricity up to $\sim 0.3$ (except for a single outlier that grazes $10^{-2}$). We finally use test-particle data to validate the waveform factorization proposed by Khalil et al.~[Phys.~Rev.~104 (2021) 2, 024046] and conclude that its amplitude can be considered reliable (though less accurate, $\sim 6\%$ fractional difference versus $1.5\%$ of our method) only up to eccentricities $\sim 0.3$.

preprint2022arXiv

New Avenue for Accurate Analytical Waveforms and Fluxes for Eccentric Compact Binaries

We introduce a new paradigm for constructing accurate analytic waveforms (and fluxes) for eccentric compact binaries. Our recipe builds on the standard Post-Newtonian (PN) approach but (i) retains implicit time-derivatives of the phase space variables in the instantaneous part of the noncircular waveform, and then (ii) suitably factorizes and resums this partly PN-implicit waveform using effective-one-body (EOB) procedures. We test our prescription against the exact results obtained by solving the Teukolsky equation with a test-mass source orbiting a Kerr black hole, and compare the use of the exact vs PN equations of motion for the time derivatives computation. Focusing only on the quadrupole contribution, we find that the use of the exact equations of motion yields an analytical/numerical agreement of the (averaged) angular momentum fluxes that is improved by $40\%$ with respect to previous work, with $4.5\%$ fractional difference for eccentricity $e=0.9$ and black hole dimensionless spin $-0.9\leq \hat{a}\leq +0.9$. We also find a remarkable convergence trend between Newtonian, 1PN and 2PN results. Our approach carries over to the comparable mass case using the resummed EOB equations of motion and paves the way to faithful EOB-based waveform model for long-inspiral eccentric binaries for current and future gravitational wave detectors.

preprint2022arXiv

Numerical-relativity validation of effective-one-body waveforms in the intermediate-mass-ratio regime

One of the open problems in developing binary black hole (BBH) waveforms for gravitational wave astronomy is to model the intermediate mass ratio regime and connect it to the extreme mass ratio regime. A natural approach is to employ the effective one body (EOB) approach to the two-body dynamics that, by design, can cover the entire mass ratio range and naturally incorporates the extreme mass ratio limit. Here we use recently obtained numerical relativity (NR) data with mass ratios $m_1/m_2=(7,15,\,32,\,64,\,128)$ to test the accuracy of the state-of-the-art EOB model TEOBResumS in the intermediate mass ratio regime. We generally find an excellent EOB/NR consistency around merger and ringdown for all mass ratios and for all available subdominant multipoles, except for the $\ell=m=5$ one. This mode can be crucially improved using the new large-mass ratio NR data of this paper. The EOB/NR inspirals are also consistent with the estimated NR uncertainties. We also use several NR datasets taken by different public catalogs to probe the universal behavior of the multipolar hierarchy of waveform amplitudes at merger, that smoothly connects the equal-mass BBH to the test-mass result. Interestingly, the universal behavior is strengthened if the nonoscillatory memory contribution is included in the NR waveform. Future NR simulations with improved accuracy will be necessary to further probe, and possibly quantitatively refine, the TEOBResumS transition from late inspiral to plunge in the intermediate mass ratio regime.

preprint2021arXiv

Waveforms and fluxes: Towards a self-consistent effective one body waveform model for nonprecessing, coalescing black-hole binaries for third generation detectors

We present a comprehensive comparison between numerical relativity (NR) angular momentum fluxes at infinity and the corresponding quantity entering the radiation reaction in TEOBResumS, an Effective-One-Body (EOB) waveform model for nonprecessing coalescing black hole binaries on quasi-circular orbits. This comparison prompted us to implement two changes in the model: (i) including Next-to-Quasi-Circular corrections in the $\ell=m$, $\ell\leq 5$ multipoles entering the radiation reaction and (ii) consequently updating the NR-informed spin-orbital sector of the model. This yields a new waveform model that presents a higher self-consistency between waveform and dynamics and an improved agreement with NR simulations. We test the model computing the EOB/NR unfaithfulness $\bar{F}_{\rm EOB/NR}$ over all 534 spin-aligned configurations available through the Simulating eXtreme Spacetime catalog, notably using the noise spectral density of Advanced LIGO, Einstein Telescope and Cosmic Explorer, for total mass up to $500M_\odot$. We find that the maximum unfaithfulness $\bar{F}^{\rm max}_{\rm EOB/NR}$ is mostly between $10^{-4}$ and $10^{-3}$, and the performance progressively worsens up to $\sim 5\times 10^{-3}$ as the effective spin of the system is increased. We perform similar analyses on the \SEOB{} model, that delivers $\bar{F}^{\rm max}_{\rm EOB/NR}$ values uniformly distributed versus effective spin and mostly between $10^{-3}$ and $10^{-2}$. We conclude that the improved TEOBResumS model already represents a reliable and robust first step towards the development of highly accurate waveform templates for third generation detectors.

preprint2020arXiv

A multipolar effective one body waveform model for spin-aligned black hole binaries

We introduce \TEOBiResumSM{}, an improved version of the effective-one-body (EOB) waveform model \TEOBResumS{} for spin-aligned, coalescing black hole binaries, that includes subdominant gravitational waveform modes completed through merger and ringdown. Beyond the dominant $(\ell,|m|)=(2,2)$ one, the more robust multipoles all over the parameter space are: $(2,1)$, $(3,3)$, $(3,2)$, $(4,4)$ and $(5,5)$. The multipolar ringdown EOB waveform stems from suitably fitting many numerical relativity (NR) waveform data from the Simulating eXtreme Spacetimes (SXS) collaboration together with test-mass waveform data. Mode-mixing effects are not incorporated. The orbital (nonspinning) part of the multipolar waveform amplitudes includes test-mass results up to (relative) 6PN order and, for most modes,is Padé resummed. The $m$=odd waveform multipoles (up to $\ell=5$) incorporate most of the currently available spin-dependent analytical information. Improving on previous work, we confirm that certain $m=\text{odd}$ modes, e.g. the $(2,1)$, and even the $(3,1)$, may develop a zero (or a minimum) in the amplitude for nearly equal-mass binaries and for several combinations of the individual spins. A remarkable EOB/NR agreement around such zero is found for these modes. The new waveform, and radiation reaction, prompts a new NR-calibration of the spinning sector of the model, done with only $32$ datasets.The maximum $(2,2)$ EOB/NR unfaithfulness $\bar{F}$ with Advanced LIGO noise against the SXS catalog ($\sim 595$ datasets) is always below $0.5\%$ for binaries with total mass $M$ as $10M_\odot\leq M \leq 200M_\odot$, except for a single outlier with $\max{(\bar{F})}\sim 0.85\%$. When $(2,1)$, $(3,3)$ and $(4,4)$ modes are included, one finds an excellent EOB/NR agreement up to $M\sim 120M_\odot$, above which the performance degrades slightly and moves above $3\%$

preprint2020arXiv

Gauge-invariant Non-spherical Metric Perturbations of Schwarzschild Black-Hole Spacetimes

The theory of gauge-invariant non-spherical metric perturbations of Schwarzschild black hole spacetimes is now well established. Yet, as different notations and conventions have been used throughout the years, the literature on the subject is often confusing and sometimes confused. The purpose of this paper is to review and collect the relevant expressions related to the Regge-Wheeler and Zerilli equations for the odd and even-parity perturbations of a Schwarzschild spacetime. Special attention is paid to the form they assume in the presence of matter-sources and, for the two most popular conventions in the literature, to the asymptotic expressions and gravitational-wave amplitudes. Besides pointing out some inconsistencies in the literature, the expressions collected here could serve as a quick reference for the calculation of the perturbations of Schwarzschild black hole spacetimes driven by generic sources and for those approaches in which gravitational waves are extracted from numerically generated spacetimes.

preprint2020arXiv

kiloHertz gravitational waves from binary neutron star remnants: time-domain model and constraints on extreme matter

The remnant star of a neutron star merger is an anticipated loud source of kiloHertz gravitational waves that conveys unique information on the equation of state of hot matter at extreme densities. Observations of such signals are hampered by the photon shot noise of ground-based interferometers and pose a challenge for gravitational-wave astronomy. We develop an analytical time-domain waveform model for postmerger signals informed by numerical relativity simulations. The model completes effective-one-body waveforms for quasi-circular nonspinning binaries in the kiloHertz regime. We show that a template-based analysis can detect postmerger signals with a minimal signal-to-noise ratios (SNR) of 8, corresponding to GW170817-like events for third-generation interferometers. Using Bayesian model selection and the complete inspiral-merger-postmerger waveform model it is possible to infer whether the merger outcome is a prompt collapse to a black hole or a remnant star. In the latter case, the radius of the maximum mass (most compact) nonrotating neutron star can be determined to kilometer precision. We demonstrate the feasibility of inferring the stiffness of the equation of state at extreme densities using the quasiuniversal relations deduced from numerical-relativity simulations.

preprint2020arXiv

Machine Learning Gravitational Waves from Binary Black Hole Mergers

We apply machine learning methods to build a time-domain model for gravitational waveforms from binary black hole mergers, called mlgw. The dimensionality of the problem is handled by representing the waveform's amplitude and phase using a principal component analysis. We train mlgw on about $\mathcal{O}(10^3)$ TEOBResumS and SEOBNRv4 effective-one-body waveforms with mass ratios $q\in[1,20]$ and aligned dimensionless spins $s\in[-0.80,0.95]$. The resulting models are faithful to the training sets at the ${\sim}10^{-3}$ level (averaged on the parameter space). The speed up for a single waveform generation is a factor 10 to 50 (depending on the binary mass and initial frequency) for TEOBResumS and approximately an order of magnitude more for SEOBNRv4. Furthermore, mlgw provides a closed form expression for the waveform and its gradient with respect to the orbital parameters; such an information might be useful for future improvements in GW data analysis. As demonstration of the capabilities of mlgw to perform a full parameter estimation, we re-analyze the public data from the first GW transient catalog (GWTC-1). We find broadly consistent results with previous analyses at a fraction of the cost, although the analysis with spin aligned waveforms gives systematic larger values of the effective spins with respect to previous analyses with precessing waveforms. Since the generation time does not depend on the length of the signal, our model is particularly suitable for the analysis of the long signals that are expected to be detected by third-generation detectors. Future applications include the analysis of waveform systematics and model selection in parameter estimation.

preprint2019arXiv

A Multipolar Effective One Body Model for Non-Spinning Black Hole Binaries

We introduce \TEOBiResumSM{}, a nonspinning inspiral-merger-ringdown waveform model built within the effective one body (EOB) framework that includes gravitational waveform modes beyond the dominant quadrupole $(\ell,|m|) = (2,2)$. The model incorporates: (i) an improved Padé resummation of the factorized waveform amplitudes $ρ_{\ell m}^{\rm orb}$ entering the EOB-resummed waveform where the 3PN, mass-ratio dependent, terms are hybridized with test-mass limit terms up to 6PN relative order for most of the multipoles up to $\ell=6$ included; (ii) an improved determination of the effective 5PN function $a_6^c(ν)$ entering the EOB interaction potential done using the most recent, error-controlled, nonspinning numerical relativity (NR) waveforms from the Simulating eXtreme Spacetimes (SXS) collaboration; and (iii) a NR-informed phenomenological description of the multipolar ringdown. Such representation stems from 19 NR waveforms with mass ratios up to $m_1/m_2=18$ as well as test-mass waveform data, although it does not incorporate mode-mixing effects. The NR-completed higher modes through merger and ringdown considered here are: $(\ell,|m|) = \lbrace (2,1), (3,3), (3,2),(3,1),(4,4), (4,3),(4,2), (4,1),(5,5)\rbrace$. For simplicity, the other subdominant modes, up to $\ell=8$, are approximated by the corresponding, purely analytical, factorized and resummed EOB waveform. To attempt an estimate of (some of) the underlying analytic uncertainties of the model, we also contrast the effect of the 6PN-hybrid Padé-resummed $ρ_{\ell m}$'s with the standard $3^{+2}$PN, Taylor-expanded, ones used in previous EOB works. The maximum unfaithfulness $\bar{F}$ against the SXS waveforms including all NR-completed modes up to $\ell=m=5$ is always $\lesssim 2\%$ for binaries with total mass $M$ as $50 M_{\odot} \leq M \lesssim 200 M_{\odot}$.

preprint2019arXiv

Comparing Effective One Body Hamiltonians for spin-aligned coalescing binaries

TEOBResumS and SEOBNRv4 are the two existing semi-analytical gravitational waveform models for spin-aligned coalescing black hole binaries based on the effective-one-body approach.They are informed by numerical relativity simulations and provide the relative dynamics and waveforms from early inspiral to plunge, merger and ringdown The central building block of each model is the EOB resummed Hamiltonian.The two models implement different Hamiltonians that are both deformations of the Hamiltonian of a test spinning black hole moving around a Kerr black hole.Here we analytically compare, element by element, the two Hamiltonians. In particular: we illustrate that one can introduce a centrifugal radius SEOBNRv4, so to rewrite the Hamiltonian in a more compact form that is analogous to the one of TEOBResumS.The latter centrifugal radius cannot, however, be identified with the one used in TEOBResumS because the two models differ in their ways of incorporating spin effects in their respective deformations of the background Kerr Hamiltonian. We performed extensive comparisons between the energetics corresponding to the two Hamiltonians using gauge-invariant quantities. Finally, as an exploratory investigation, we apply the post-adiabatic approximation to the newly rewritten SEOBNRv4 Hamiltonian, illustrating that it is possible to generate long-inspiral waveforms with negligible computational cost.

preprint2016arXiv

iResum: a new paradigm for resumming gravitational wave amplitudes

We introduce a new, resummed, analytical form of the post-Newtonian (PN), factorized, multipolar amplitude corrections $f_{\ell m}$ of the effective-one-body (EOB) gravitational waveform of spinning, nonprecessing, circularized, coalescing black hole binaries (BBHs). This stems from the following two-step paradigm: (i) the factorization of the orbital (spin-independent) terms in $f_{\ell m}$; (ii) the resummation of the residual spin (or orbital) factors. We find that resumming the residual spin factor by taking its inverse resummed (iResum) is an efficient way to obtain amplitudes that are more accurate in the strong-field, fast-velocity regime. The performance of the method is illustrated on the $\ell=2$ and $m=(1,2)$ waveform multipoles, both for a test-mass orbiting around a Kerr black hole and for comparable-mass BBHs. In the first case, the iResum $f_{\ell m}$'s are much closer to the corresponding "exact" functions (obtained solving numerically the Teukolsky equation) up to the light-ring, than the nonresummed ones, especially when the black-hole spin is nearly extremal. The iResum paradigm is also more efficient than including higher post-Newtonian terms (up to 20PN order): the resummed 5PN information yields per se a rather good numerical/analytical agreement at the last-stable-orbit, and a well-controlled behavior up to the light-ring. For comparable mass binaries (including the highest PN-order information available, 3.5PN), comparing EOB with Numerical Relativity (NR) data shows that the analytical/numerical fractional disagreement at merger, without NR-calibration of the EOB waveform, is generically reduced by iResum, from a $40\%$ of the usual approach to just a few percents. This suggests that EOBNR waveform models for coalescing BBHs may be improved using iResum amplitudes.

preprint2015arXiv

Energetics and phasing of nonprecessing spinning coalescing black hole binaries

We present an improved numerical relativity (NR) calibration of the new effective-one-body (EOB) model for coalescing non precessing spinning black hole binaries recently introduced by Damour and Nagar [Physical Review D 90, 044018 (2014)]. We do so by comparing the EOB predictions to both the phasing and the energetics provided by two independent sets of NR data covering mass ratios $1\leq q \leq 9.989$ and dimensionless spin range $-0.95\leq χ\leq +0.994$. One set of data is a subset of the Simulating eXtreme Spacetimes (SXS) catalog of public waveforms; the other set consists of new simulations obtained with the Llama code plus Cauchy Characteristic Evolution. We present the first systematic computation of the gauge-invariant relation between the binding energy and the total angular momentum, $E_{b}(j)$, for a large sample of, spin-aligned, SXS and Llama data. The dynamics of the EOB model presented here involves only two free functional parameters, one ($a_6^c(ν)$) entering the non spinning sector, as a 5PN effective correction to the interaction potential, and one ($c_3(\tilde{a}_1,\tilde{a}_2,ν))$ in the spinning sector, as an effective next-to-next-to-next-to-leading order correction to the spin-orbit coupling. These parameters are determined (together with a third functional parameter $Δt_{\rm NQC}(χ)$ entering the waveform) by comparing the EOB phasing with the SXS phasing, the consistency of the energetics being checked afterwards. The quality of the analytical model for gravitational wave data analysis purposes is assessed by computing the EOB/NR faithfulness. Over the NR data sample and when varying the total mass between 20 and 200~$M_\odot$ the EOB/NR unfaithfulness (integrated over the NR frequency range) is found to vary between $99.493\%$ and $99.984\%$ with a median value of $99.944\%$.

preprint2015arXiv

Modeling the complete gravitational wave spectrum of neutron star mergers

In the context of neutron star mergers, we study the gravitational wave spectrum of the merger remnant using numerical relativity simulations. Postmerger spectra are characterized by a main peak frequency $f_2$ related to the particular structure and dynamics of the remnant hot hypermassive neutron star. We show that $f_2$ is correlated with the tidal coupling constant $κ^T_2$ that characterizes the binary tidal interactions during the late-inspiral--merger. The relation $f_2(κ^T_2)$ depends very weakly on the binary total mass, mass-ratio, equation of state, and thermal effects. This observation opens up the possibility of developing a model of the gravitational spectrum of every merger unifying the late-inspiral and postmerger descriptions.

preprint2015arXiv

Modeling the dynamics of tidally-interacting binary neutron stars up to merger

The data analysis of the gravitational wave signals emitted by coalescing neutron star binaries requires the availability of an accurate analytical representation of the dynamics and waveforms of these systems. We propose an effective-one-body (EOB) model that describes the general relativistic dynamics of neutron star binaries from the early inspiral up to merger. Our EOB model incorporates an enhanced attractive tidal potential motivated by recent analytical advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. No fitting parameters are introduced for the description of tidal interaction in the late, strong-field dynamics. We compare the model energetics and the gravitational wave phasing with new high-resolution multi-orbit numerical relativity simulations of equal-mass configurations with different equations of state. We find agreement within the uncertainty of the numerical data for all configurations. Our model is the first semi-analytical model which captures the tidal amplification effects close to merger. It thereby provides the most accurate analytical representation of binary neutron star dynamics and waveforms currently available.

preprint2014arXiv

A new analytic representation of the ringdown waveform of coalescing spinning black hole binaries

We propose a new way of analyzing, and analytically representing, the ringdown part of the gravitational wave signal emitted by coalescing black hole binaries.By contrast with the usual {\it linear} decomposition of the multipolar complex waveform $h(t)$ in a sum of quasi-normal modes, our procedure relies on a {\it multiplicative} decomposition of $h(t)$ as the product of the fundamental quasi-normal mode with a remaining time-dependent complex factor whose amplitude and phase are separately fitted. As an illustrative example, we apply our analysis and fitting procedure to the ringdown part of a sample of sixteen $\ell=m=2$ equal-mass, spinning, nonprecessing, numerical waveforms computed with the SP$_{\rm E}$C code,now publicly available in the SXS catalogue. Our approach yields an efficient and accurate way to represent the ringdown waveform, thereby offering a new way to complete the analytical effective-one-body inspiral-plus-plunge waveform.

preprint2014arXiv

A new effective-one-body description of coalescing nonprecessing spinning black-hole binaries

We present a new, tunable effective-one-body (EOB) model of the motion and radiation of coalescing black hole binaries with arbitrary mass ratio and aligned spins. The most novel feature of our formalism is the introduction, and systematic use, of the (gauge-invariant) concept of {\it centrifugal radius} $r_{c}$. In the spinning small mass-ratio limit, the main radial potential expressed in terms of $r_{c}$ differs very little (and only multiplicatively so) from the usual Schwarzschild potential $1-2M/r_{c}$. This motivates a new, multiplicative way of blending finite-mass-ratio deformations with spin-deformations. In the present exploratory work we consider a minimal version of our spinning EOB model (containing essentially only two adjustable parameters: one in the Hamiltonian and one in the waveform) and calibrate its (dominant mode) waveform against a sample of fifteen equal-mass, equal-spin waveforms produced by the SXS collaboration, and covering the dimensionless spin range $-0.95\leq χ\leq +0.98$. The numerical relativity / EOB phasing disagreement remains remarkably small ($\lesssim \pm 0.15$ rad) over the entire spin range.

preprint2014arXiv

A new gravitational wave generation algorithm for particle perturbations of the Kerr spacetime

We present a new approach to solve the 2+1 Teukolsky equation for gravitational perturbations of a Kerr black hole. Our approach relies on a new horizon penetrating, hyperboloidal foliation of Kerr spacetime and spatial compactification. In particular, we present a framework for waveform generation from point-particle perturbations. Extensive tests of a time domain implementation in the code {\it Teukode} are presented. The code can efficiently deliver waveforms at future null infinity. As a first application of the method, we compute the gravitational waveforms from inspiraling and coalescing black-hole binaries in the large-mass-ratio limit. The smaller mass black hole is modeled as a point particle whose dynamics is driven by an effective-one-body-resummed analytical radiation reaction force. We compare the analytical angular momentum loss to the gravitational wave angular momentum flux. We find that higher-order post-Newtonian corrections are needed to improve the consistency for rapidly spinning binaries. Close to merger, the subdominant multipolar amplitudes (notably the $m=0$ ones) are enhanced for retrograde orbits with respect to prograde ones. We argue that this effect mirrors nonnegligible deviations from circularity of the dynamics during the late-plunge and merger phase. We compute the gravitational wave energy flux flowing into the black hole during the inspiral using a time-domain formalism proposed by Poisson. Finally, a self-consistent, iterative method to compute the gravitational wave fluxes at leading-order in the mass of the particle is presented. For a specific case study with $\hat{a}$=0.9, a simulation that uses the consistent flux differs from one that uses the analytical flux by $\sim35$ gravitational wave cycles over a total of about $250$ cycles. In this case the horizon absorption accounts for about $+5$ gravitational wave cycles.

preprint2014arXiv

Quasiuniversal properties of neutron star mergers

Binary neutron star mergers are studied using nonlinear 3+1 numerical relativity simulations and the analytical effective-one-body (EOB) model. The EOB model predicts quasiuniversal relations between the mass-rescaled gravitational wave frequency and the binding energy at the moment of merger, and certain dimensionless binary tidal coupling constants depending on the stars Love numbers, compactnesses and the binary mass ratio. These relations are quasiuniversal in the sense that, for a given value of the tidal coupling constant, they depend significantly neither on the equation of state nor on the mass ratio, though they do depend on stars spins. The spin dependence is approximately linear for small spins aligned with the orbital angular momentum. The quasiuniversality is a property of the conservative dynamics; nontrivial relations emerge as the binary interaction becomes tidally dominated. This analytical prediction is qualitatively consistent with new, multi-orbit numerical relativity results for the relevant case of equal-mass irrotational binaries. Universal relations are thus expected to characterize neutron star mergers dynamics. In the context of gravitational wave astronomy, these universal relations may be used to constrain the neutron star equation of state using waveforms that model the merger accurately.

preprint2014arXiv

Strong-Field Scattering of Two Black Holes: Numerics Versus Analytics

We probe the gravitational interaction of two black holes in the strong-field regime by computing the scattering angle $χ$ of hyperbolic-like, close binary-black-hole encounters as a function of the impact parameter. The fully general-relativistic result from numerical relativity is compared to two analytic approximations: post-Newtonian theory and the effective-one-body formalism. As the impact parameter decreases, so that black holes pass within a few times their Schwarzschild radii, we find that the post-Newtonian prediction becomes quite inaccurate, while the effective-one-body one keeps showing a good agreement with numerical results. Because we have explored a regime which is very different from the one considered so far with binaries in quasi-circular orbits, our results open a new avenue to improve analytic representations of the general-relativistic two-body Hamiltonian.

preprint2014arXiv

The antikick strikes back: recoil velocities for nearly-extremal binary black hole mergers in the test-mass limit

Gravitational waves emitted from a generic binary black-hole merger carry away linear momentum anisotropically, resulting in a gravitational recoil, or "kick", of the center of mass. For certain merger configurations the time evolution of the magnitude of the kick velocity has a local maximum followed by a sudden drop. Perturbative studies of this "antikick" in a limited range of black hole spins have found that the antikick decreases for retrograde orbits as a function of negative spin. We analyze this problem using a recently developed code to evolve gravitational perturbations from a point-particle in Kerr spacetime driven by an effective-one-body resummed radiation reaction force at linear order in the mass ratio $ν\ll 1$. Extending previous studies to nearly-extremal negative spins, we find that the well-known decrease of the antikick is overturned and, instead of approaching zero, the antikick increases again to reach $Δv/(cν^{2})=3.37\times10^{-3}$ for dimensionless spin $\hat{a}=-0.9999$. The corresponding final kick velocity is $v_{end}/(cν^{2})=0.076$. This result is connected to the nonadiabatic character of the emission of linear momentum during the plunge. We interpret it analytically by means of the quality factor of the flux to capture quantitatively the main properties of the kick velocity. The use of such quality factor of the flux does not require trajectories nor horizon curvature distributions and should therefore be useful both in perturbation theory and numerical relativity.

preprint2013arXiv

Effective one body Hamiltonian of two spinning black-holes with next-to-next-to-leading order spin-orbit coupling

Building on the recently computed next-to-next-to-leading order (NNLO) post-Newtonian (PN) spin-orbit Hamiltonian for spinning binaries \cite{Hartung:2011te} we extend the effective-one-body (EOB) description of the dynamics of two spinning black-holes to NNLO in the spin-orbit interaction. The calculation that is presented extends to NNLO the next-to-leading order (NLO) spin-orbit Hamiltonian computed in Ref. \cite{Damour:2008qf}. The present EOB Hamiltonian reproduces the spin-orbit coupling through NNLO in the test-particle limit case. In addition, in the case of spins parallel or antiparallel to the orbital angular momentum, when circular orbits exist, we find that the inclusion of NNLO spin-orbit terms moderates the effect of the NLO spin-orbit coupling.

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

Gravitational recoil in nonspinning black hole binaries: the span of test-mass results

We consider binary systems of coalescing, nonspinning, black holes of masses $m_{1}$ and $m_{2}$ and show that the gravitational recoil velocity for any mass ratio can be obtained accurately by extrapolating the waveform of the test-mass limit case. The waveform obtained in the limit $m_1/m_2\ll 1$ via a perturbative approach is extrapolated in $ν= m_{1} m_{2}/(m_{1}+m_{2})^{2} $ multipole by multipole using the corresponding, analytically known, leading-in-$ν$ behavior. The final kick velocity computed from this $ν$-flexed waveform is written as $v(ν)/c = 0.04457 ν^2\sqrt{1-4ν}\,(1-2.07106ν+ 3.93472ν^2 -4.78404ν^3+2.52040ν^{4})$ and is compatible with the outcome of numerical relativity simulations

preprint2013arXiv

Improved effective-one-body description of coalescing nonspinning black-hole binaries and its numerical-relativity completion

We improve the effective-one-body (EOB) description of nonspinning coalescing black hole binaries by incorporating several recent analytical advances, notably: (i) logarithmic contributions to the conservative dynamics; (ii) resummed horizon-absorption contribution to the orbital angular momentum loss; and (iii) a specific radial component of the radiation reaction force implied by consistency with the azimuthal one. We then complete this analytically improved EOB model by comparing it to accurate numerical relativity (NR) simulations performed by the Caltech-Cornell-CITA group for mass ratios $q=(1,2,3,4,6)$. In particular, the comparison to NR data allows us to determine with high-accuracy ($\sim 10^{-4}$) the value of the main EOB radial potential: $A(u;\,ν)$, where $u=GM/(R c^2)$ is the inter-body gravitational potential and $ν=q/(q+1)^2$ is the symmetric mass ratio. We introduce a new technique for extracting from NR data an intrinsic measure of the phase evolution, ($Q_ω(ω)$ diagnostics). Aligning the NR-completed EOB quadrupolar waveform and the NR one at low frequencies, we find that they keep agreeing (in phase and amplitude) within the NR uncertainties throughout the evolution for all mass ratios considered. We also find good agreement for several subdominant multipoles without having to introduce and tune any extra parameters.

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

Horizon-absorption effects in coalescing black-hole binaries: An effective-one-body study of the non-spinning case

We study the horizon absorption of gravitational waves in coalescing, circularized, nonspinning black hole binaries. The horizon absorbed fluxes of a binary with a large mass ratio (q=1000) obtained by numerical perturbative simulations are compared with an analytical, effective-one-body (EOB) resummed expression recently proposed. The perturbative method employs an analytical, linear in the mass ratio, effective-one-body (EOB) resummed radiation reaction, and the Regge-Wheeler-Zerilli (RWZ) formalism for wave extraction. Hyperboloidal (transmitting) layers are employed for the numerical solution of the RWZ equations to accurately compute horizon fluxes up to the late plunge phase. The horizon fluxes from perturbative simulations and the EOB-resummed expression agree at the level of a few percent down to the late plunge. An upgrade of the EOB model for nonspinning binaries that includes horizon absorption of angular momentum as an additional term in the resummed radiation reaction is then discussed. The effect of this term on the waveform phasing for binaries with mass ratios spanning 1 to 1000 is investigated. We confirm that for comparable and intermediate-mass-ratio binaries horizon absorbtion is practically negligible for detection with advanced LIGO and the Einstein Telescope (faithfulness greater than or equal to 0.997).

preprint2012arXiv

Measurability of the tidal polarizability of neutron stars in late-inspiral gravitational-wave signals

The gravitational wave signal from a binary neutron star inspiral contains information on the nuclear equation of state. This information is contained in a combination of the tidal polarizability parameters of the two neutron stars and is clearest in the late inspiral, just before merger. We use the recently defined tidal extension of the effective one-body formalism to construct a controlled analytical description of the frequency-domain phasing of neutron star inspirals up to merger. Exploiting this analytical description we find that the tidal polarizability parameters of neutron stars can be measured by the advanced LIGO-Virgo detector network from gravitational wave signals having a reasonable signal-to-noise ratio of $ρ=16$. This measurability result seems to hold for all the nuclear equations of state leading to a maximum mass larger than $1.97M_\odot$. We also propose a promising new way of extracting information on the nuclear equation of state from a coherent analysis of an ensemble of gravitational wave observations of separate binary merger events.

preprint2012arXiv

Tidal effects in binary neutron star coalescence

We compare dynamics and waveforms from binary neutron star coalescence as computed by new long-term ($\sim 10 $ orbits) numerical relativity simulations and by the tidal effective-one-body (EOB) model including analytical tidal corrections up to second post-Newtonian order (2PN). The current analytical knowledge encoded in the tidal EOB model is found to be sufficient to reproduce the numerical data up to contact and within their uncertainties. Remarkably, no calibration of any tidal EOB free parameters is required, beside those already fitted to binary black holes data. The inclusion of 2PN tidal corrections minimizes the differences with the numerical data, but it is not possible to significantly distinguish them from the leading-order tidal contribution. The presence of a relevant amplification of tidal effects is likely to be excluded, although it can appear as a consequence of numerical inaccuracies. We conclude that the tidally-completed effective-one-body model provides nowadays the most advanced and accurate tool for modelling gravitational waveforms from binary neutron star inspiral up to contact. This work also points out the importance of extensive tests to assess the uncertainties of the numerical data, and the potential need of new numerical strategies to perform accurate simulations.

preprint2011arXiv

Accurate numerical simulations of inspiralling binary neutron stars and their comparison with effective-one-body analytical models

Binary neutron-star systems represent one of the most promising sources of gravitational waves. In order to be able to extract important information, notably about the equation of state of matter at nuclear density, it is necessary to have in hands an accurate analytical model of the expected waveforms. Following our recent work, we here analyze more in detail two general-relativistic simulations spanning about 20 gravitational-wave cycles of the inspiral of equal-mass binary neutron stars with different compactnesses, and compare them with a tidal extension of the effective-one-body (EOB) analytical model. The latter tidally extended EOB model is analytically complete up to the 1.5 post-Newtonian level, and contains an analytically undetermined parameter representing a higher-order amplification of tidal effects. We find that, by calibrating this single parameter, the EOB model can reproduce, within the numerical error, the two numerical waveforms essentially up to the merger. By contrast, analytical models (either EOB, or Taylor-T4) that do not incorporate such a higher-order amplification of tidal effects, build a dephasing with respect to the numerical waveforms of several radians.

preprint2011arXiv

Binary black hole coalescence in the extreme-mass-ratio limit: testing and improving the effective-one-body multipolar waveform

We discuss the properties of the effective-one-body (EOB) multipolar gravitational waveform emitted by nonspinning black-hole binaries of masses $μ$ and $M$ in the extreme-mass-ratio limit, $μ/M=ν\ll 1$. We focus on the transition from quasicircular inspiral to plunge, merger and ringdown.We compare the EOB waveform to a Regge-Wheeler-Zerilli (RWZ) waveform computed using the hyperboloidal layer method and extracted at null infinity. Because the EOB waveform keeps track analytically of most phase differences in the early inspiral, we do not allow for any arbitrary time or phase shift between the waveforms. The dynamics of the particle, common to both wave-generation formalisms, is driven by leading-order ${\cal O}(ν)$ analytically--resummed radiation reaction. The EOB and the RWZ waveforms have an initial dephasing of about $5\times 10^{-4}$ rad and maintain then a remarkably accurate phase coherence during the long inspiral ($\sim 33$ orbits), accumulating only about $-2\times 10^{-3}$ rad until the last stable orbit, i.e. $Δϕ/ϕ\sim -5.95\times 10^{-6}$. We obtain such accuracy without calibrating the analytically-resummed EOB waveform to numerical data, which indicates the aptitude of the EOB waveform for LISA-oriented studies. We then improve the behavior of the EOB waveform around merger by introducing and tuning next-to-quasi-circular corrections both in the gravitational wave amplitude and phase. For each multipole we tune only four next-to-quasi-circular parameters by requiring compatibility between EOB and RWZ waveforms at the light-ring. The resulting phase difference around merger time is as small as $\pm 0.015$ rad, with a fractional amplitude agreement of 2.5%. This suggest that next-to-quasi-circular corrections to the phase can be a useful ingredient in comparisons between EOB and numerical relativity waveforms.

preprint2011arXiv

Binary black hole coalescence in the large-mass-ratio limit: the hyperboloidal layer method and waveforms at null infinity

We compute and analyze the gravitational waveform emitted to future null infinity by a system of two black holes in the large mass ratio limit. We consider the transition from the quasi-adiabatic inspiral to plunge, merger, and ringdown. The relative dynamics is driven by a leading order in the mass ratio, 5PN-resummed, effective-one-body (EOB), analytic radiation reaction. To compute the waveforms we solve the Regge-Wheeler-Zerilli equations in the time-domain on a spacelike foliation which coincides with the standard Schwarzschild foliation in the region including the motion of the small black hole, and is globally hyperboloidal, allowing us to include future null infinity in the computational domain by compactification. This method is called the hyperboloidal layer method, and is discussed here for the first time in a study of the gravitational radiation emitted by black hole binaries. We consider binaries characterized by five mass ratios, $ν=10^{-2,-3,-4,-5,-6}$, that are primary targets of space-based or third-generation gravitational wave detectors. We show significative phase differences between finite-radius and null-infinity waveforms. We test, in our context, the reliability of the extrapolation procedure routinely applied to numerical relativity waveforms. We present an updated calculation of the gravitational recoil imparted to the merger remnant by the gravitational wave emission. As a self consistency test of the method, we show an excellent fractional agreement (even during the plunge) between the 5PN EOB-resummed mechanical angular momentum loss and the gravitational wave angular momentum flux computed at null infinity. New results concerning the radiation emitted from unstable circular orbits are also presented.

preprint2011arXiv

Energy versus Angular Momentum in Black Hole Binaries

Using accurate numerical relativity simulations of (nonspinning) black-hole binaries with mass ratios 1:1, 2:1 and 3:1 we compute the gauge invariant relation between the (reduced) binding energy $E$ and the (reduced) angular momentum $j$ of the system. We show that the relation $E(j)$ is an accurate diagnostic of the dynamics of a black-hole binary in a highly relativistic regime. By comparing the numerical-relativity $E^{\rm NR} (j)$ curve with the predictions of several analytic approximation schemes, we find that, while the usual, non-resummed post-Newtonian-expanded $E^{\rm PN} (j)$ relation exhibits large and growing deviations from $E^{\rm NR} (j)$, the prediction of the effective one-body formalism, based purely on known analytical results (without any calibration to numerical relativity), agrees strikingly well with the numerical-relativity results.

preprint2011arXiv

Gravitational-wave extraction from neutron-star oscillations

We compare different gravitational-wave extraction methods used in three-dimensional nonlinear simulations against linear simulations of perturbations of spherical spacetimes with matter. We present results from fully general-relativistic simulations of a system composed by an oscillating and non-rotating star emitting gravitational radiation. Results about the onset of non-linear effects are also shown.

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 "effective source" 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

Analytic modelling of tidal effects in the relativistic inspiral of binary neutron stars

To detect the gravitational-wave (GW) signal from binary neutron stars and extract information about the equation of state of matter at nuclear density, it is necessary to match the signal with a bank of accurate templates. We present the two longest (to date) general-relativistic simulations of equal-mass binary neutron stars with different compactnesses, C=0.12 and C=0.14, and compare them with a tidal extension of the effective-one-body (EOB)model. The typical numerical phasing errors over the $\simeq 22$ GW cycles are $Δϕ\simeq \pm 0.24$ rad. By calibrating only one parameter (representing a higher-order amplification of tidal effects), the EOB model can reproduce, within the numerical error, the two numerical waveforms essentially up to the merger. By contrast, the third post-Newtonian Taylor-T4 approximant with leading-order tidal corrections dephases with respect to the numerical waveforms by several radians.

preprint2010arXiv

Binary black hole merger in the extreme-mass-ratio limit: a multipolar analysis

Building up on previous work, we present a new calculation of the gravitational wave (GW) emission generated during the transition from quasi-circular inspiral to plunge, merger and ringdown by a binary system of nonspinning black holes, of masses $m_1$ and $m_2$, in the extreme mass ratio limit, $m_1 m_2\ll(m_1+m_2)^2$. The relative dynamics of the system is computed {\it without making any adiabatic approximation} by using an effective one body (EOB) description, namely by representing the binary by an effective particle of mass $μ=m_1 m_2/(m_1+m_2)$ moving in a (quasi-)Schwarzschild background of mass $M=m_1+m_2$ and submitted to an $Ø(ν)$ 5PN-resummed analytical radiation reaction force, with $ν=μ/M$. The gravitational wave emission is calculated via a multipolar Regge-Wheeler-Zerilli type perturbative approach (valid in the limit $ν\ll 1$). We consider three mass ratios, $ν={10^{-2},10^{-3},10^{-4}}$,and we compute the multipolar waveform up to $\ell=8$. We estimate energy and angular momentum losses during the quasi-universal and quasi-geodesic part of the plunge phase and we analyze the structure of the ringdown. We calculate the gravitational recoil, or "kick", imparted to the merger remnant by the gravitational wave emission and we emphasize the importance of higher multipoles to get a final value of the recoil $v/(cν^2)=0.0446$. We finally show that there is an {\it excellent fractional agreement} ($\sim 10^{-3}$) (even during the plunge) between the 5PN EOB analytically-resummed radiation reaction flux and the numerically computed gravitational wave angular momentum flux. This is a further confirmation of the aptitude of the EOB formalism to accurately model extreme-mass-ratio inspirals, as needed for the future space-based LISA gravitational wave detector.

preprint2009arXiv

Effective One Body description of tidal effects in inspiralling compact binaries

The late part of the gravitational wave signal of binary neutron star inspirals can in principle yield crucial information on the nuclear equation of state via its dependence on relativistic tidal parameters. In the hope of analytically describing the gravitational wave phasing during the late inspiral (essentially up to contact) we propose an extension of the effective one body (EOB) formalism which includes tidal effects. We compare the prediction of this tidal-EOB formalism to recently computed nonconformally flat quasi-equilibrium circular sequences of binary neutron star systems. Our analysis suggests the importance of higher-order (post-Newtonian) corrections to tidal effects, even beyond the first post-Newtonian order, and their tendency to {\it significantly} increase the ``effective tidal polarizability'' of neutron stars. We compare the EOB predictions to some recently advocated, nonresummed, post-Newtonian based (``Taylor-T4'') description of the phasing of inspiralling systems. This comparison shows the strong sensitivity of the late-inspiral phasing to the choice of the analytical model, but raises the hope that a sufficiently accurate numerical--relativity--``calibrated'' EOB model might give us a reliable handle on the nuclear equation of state

preprint2009arXiv

Improved resummation of post-Newtonian multipolar waveforms from circularized compact binaries

We improve and generalize a resummation method of post-Newtonian multipolar waveforms from circular compact binaries introduced in Refs. \cite{Damour:2007xr,Damour:2007yf}. One of the characteristic features of this resummation method is to replace the usual {\it additive} decomposition of the standard post-Newtonian approach by a {\it multiplicative} decomposition of the complex multipolar waveform $h_{\lm}$ into several (physically motivated) factors: (i) the "Newtonian" waveform, (ii) a relativistic correction coming from an "effective source", (iii) leading-order tail effects linked to propagation on a Schwarzschild background, (iv) a residual tail dephasing, and (v) residual relativistic amplitude corrections $f_{\lm}$. We explore here a new route for resumming $f_{\lm}$ based on replacing it by its $\ell$-th root: $ρ_{\lm}=f_{\lm}^{1/\ell}$. In the extreme-mass-ratio case, this resummation procedure results in a much better agreement between analytical and numerical waveforms than when using standard post-Newtonian approximants. We then show that our best approximants behave in a robust and continuous manner as we "deform" them by increasing the symmetric mass ratio $ν\equiv m_1 m_2/(m_1+m_2)^2$ from 0 (extreme-mass-ratio case) to 1/4 (equal-mass case). The present paper also completes our knowledge of the first post-Newtonian corrections to multipole moments by computing ready-to-use explicit expressions for the first post-Newtonian contributions to the odd-parity (current) multipoles.