Towards a Fourier domain waveform for non-spinning binaries with arbitrary eccentricity

Abstract

Although the gravitational waves observed by advanced LIGO and Virgo are consistent with compact binaries in a quasi-circular inspiral prior to coalescence, eccentric inspirals are also expected to occur in Nature. Due to their complexity, we currently lack ready-to-use, analytic waveforms in the Fourier domain valid for sufficiently high eccentricity, and such models are crucial to coherently extract weak signals from the noise. We here take the first steps to derive and properly validate an analytic waveform model in the Fourier domain that is valid for inspirals of arbitrary orbital eccentricity. As a proof-of-concept, we build this model to leading post-Newtonian order by combining the stationary phase approximation, a truncated sum of harmonics, and an analytic representation of hypergeometric functions. Through comparisons with numerical post-Newtonian waveforms, we determine how many harmonics are required for a faithful (matches above 99%) representation of the signal up to orbital eccentricities as large as 0.9. As a first byproduct of this analysis, we present a novel technique to maximize the match of eccentric signals over time of coalescence and phase at coalescence. As a second byproduct, we determine which of the different approximations we employ leads to the largest loss in match, which could be used to systematically improve the model because of our analytic control. The future extension of this model to higher post-Newtonian order will allow for an accurate and fast phenomenological hybrid that can account for arbitrary eccentricity inspirals and mergers.PACS numbers: 04.30.-w,04.25.-g,04.25.Nx 1 The post-Newtonian approximation is one in which the field equations are solved assuming small velocities and weak gravitational fields in an expansion in powers of ( v c ), where v is the orbital velocity and c is the speed of light [20]. By nPN order we mean an expansion to order (v/c) 2n . arXiv:1807.07163v1 [gr-qc]

Pulsar timing array signals induced by black hole binaries in relativistic eccentric orbits

“…We begin by describing our computationally efficient way to obtain n(t) and e t (t), influenced by Refs. [42,44]. Our approach involves deriving certain analytic expressions for n(e t ) and t(e t ) and appropriately treating them numerically to obtain an accurate and efficient way to track the temporal evolution in n(t) and e t (t).…”

“…[42,43], while adapting recent results from Ref. [44]. This allows us to model PTA signals due to non-spinning SMBH binaries inspiraling along 3PN-accurate eccentric orbits in a semi-analytic manner.…”

Ready-to-use Fourier domain templates for compact binaries inspiraling along moderately eccentric orbits

“…It is customary to solve these two coupled differential equations numerically to obtain ω(t) and e t (t) and hence temporally evolving h(t). Interestingly, earlier efforts provided certain analytic way for obtaining temporal evolution for ω(t) and e t (t) that mainly involves the usage of hypergeometric functions [58][59][60][61] However, it is possible to obtain analytic frequency domain counterpart of the above h(t) as demonstrated in Ref. [44,48].…”

Improved constraints on modified gravity with eccentric gravitational waves

“…We truncate the waveforms with unit step functions Θ(x) to make sure each harmonic does not exceed its region of validity (see e.g. the discussion in Appendix A of [45]).…”