Defining eccentricity for gravitational wave astronomy (2024)

Physical Review D

covering particles, fields, gravitation, and cosmology
  • Highlights
  • Recent
  • Accepted
  • Collections
  • Authors
  • Referees
  • Search
  • Press
  • About
  • Editorial Team

Defining eccentricity for gravitational wave astronomy

Md Arif Shaikh, Vijay Varma, Harald P. Pfeiffer, Antoni Ramos-Buades, and Maarten van de Meent
Phys. Rev. D 108, 104007 – Published 3 November 2023
  • Article
  • References
  • Citing Articles (4)

PDFHTMLExport Citation

Defining eccentricity for gravitational wave astronomy (1)

Abstract
Authors
Article Text
  • INTRODUCTION
  • DEFINING ECCENTRICITY
  • METHODS TO LOCATE PERICENTERS AND…
  • ROBUSTNESS TESTS
  • CONCLUSION
  • ACKNOWLEDGMENTS
  • References

    Defining eccentricity for gravitational wave astronomy (2)

    Abstract

    Eccentric compact binary mergers are significant scientific targets for current and future gravitational wave observatories. To detect and analyze eccentric signals, there is an increasing effort to develop waveform models, numerical relativity simulations, and parameter estimation frameworks for eccentric binaries. Unfortunately, current models and simulations use different internal parametrizations of eccentricity in the absence of a unique natural definition of eccentricity in general relativity, which can result in incompatible eccentricity measurements. In this paper, we adopt a standardized definition of eccentricity and mean anomaly based solely on waveform quantities and make our implementation publicly available through an easy-to-use python package, gw_eccentricity. This definition is free of gauge ambiguities, has the correct Newtonian limit, and can be applied as a postprocessing step when comparing eccentricity measurements from different models. This standardization puts all models and simulations on the same footing and enables direct comparisons between eccentricity estimates from gravitational wave observations and astrophysical predictions. We demonstrate the applicability of this definition and the robustness of our implementation for waveforms of different origins, including post-Newtonian theory, effective-one-body, extreme mass ratio inspirals, and numerical relativity simulations. We focus on binaries without spin precession in this work, but possible generalizations to spin-precessing binaries are discussed.

    • Defining eccentricity for gravitational wave astronomy (3)
    • Defining eccentricity for gravitational wave astronomy (4)
    • Defining eccentricity for gravitational wave astronomy (5)
    • Defining eccentricity for gravitational wave astronomy (6)
    • Defining eccentricity for gravitational wave astronomy (7)
    • Defining eccentricity for gravitational wave astronomy (8)
    • Defining eccentricity for gravitational wave astronomy (9)

    7 More

    • Received 23 February 2023
    • Accepted 22 September 2023

    DOI:https://doi.org/10.1103/PhysRevD.108.104007

    © 2023 American Physical Society

    Physics Subject Headings (PhySH)

    1. Research Areas

    Gravitational waves

    1. Physical Systems

    Classical black holes

    1. Techniques

    Gravitational wave detectionNumerical relativity

    Gravitation, Cosmology & Astrophysics

    Authors & Affiliations

    Md Arif Shaikh1,2,*, Vijay Varma3,4,†, Harald P. Pfeiffer3, Antoni Ramos-Buades3, and Maarten van de Meent3,5

    • 1Department of Physics and Astronomy, Seoul National University, Seoul 08826, Korea
    • 2International Centre for Theoretical Sciences, Tata Institute of Fundamental Research, Bangalore 560089, India
    • 3Max Planck Institute for Gravitational Physics (Albert Einstein Institute), D-14476 Potsdam, Germany
    • 4Department of Mathematics, Center for Scientific Computing and Data Science Research, University of Massachusetts, Dartmouth, Massachusetts 02747, USA
    • 5Niels Bohr International Academy, Niels Bohr Institute, Blegdamsvej 17, 2100 Copenhagen, Denmark
    • *arifshaikh.astro@gmail.com
    • vijay.varma@aei.mpg.de

    Article Text (Subscription Required)

    Click to Expand

    References (Subscription Required)

    Click to Expand

    Issue

    Vol. 108, Iss. 10 — 15 November 2023

    Defining eccentricity for gravitational wave astronomy (10)
    Reuse & Permissions
    Access Options
    Defining eccentricity for gravitational wave astronomy (11)

    Article part of CHORUS

    Accepted manuscript will be available starting 2 November 2024.
    Defining eccentricity for gravitational wave astronomy (14)

    Authorization Required

    Other Options
    • Buy Article »
    • Find an Institution with the Article »

    ×

    Download & Share

    PDFExportReuse & PermissionsCiting Articles (4)

    ×

    Images

    • Defining eccentricity for gravitational wave astronomy (15)

      Figure 1

      Eccentricity and mean anomaly measured using the waveform from an equal-mass nonspinning eccentric NR simulation (SXS:BBH:2312 [48, 117]). Left: time evolution of the eccentricity egw (upper panel) and frequency of the (2, 2) waveform mode ω22 (lower panel). ω22p(t) and ω22a(t) are interpolants through ω22(t) evaluated at the pericenters (blue circles) and apocenters (pink squares), respectively. Equation(8) is used to compute egw(t) given ω22p(t) and ω22a(t). Right: time evolution of the mean anomaly lgw (upper panel) and ω22 (lower panel). The vertical dashed gray lines denote the pericenter times. lgw(t) grows linearly in time from 0 to 2π between successive pericenters [Eq.(10)].

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (16)

      Figure 2

      Different methods to construct a monotonically increasing frequency to replace ω22(t), in order to set the reference frequency fref for eccentric binaries. We consider two different approaches: (i)ω22mean(t), the mean of ω22p(t) and ω22a(t), and (ii)ω22(t), an interpolant through the orbit-averaged ω22 [Eq.(12)]. We show SEOBNRv4EHM waveforms with three different eccentricities; the binary parameters are given in the figure text. While the two approaches agree for small eccentricities, they deviate significantly at large eccentricities. We adopt ω22(t) as it captures the correct frequency scale in an orbit-averaged sense (Sec.2e).

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (17)

      Figure 3

      How to truncate time-domain eccentric waveforms while retaining all frequencies above flow=20Hz. The orange, blue, and pink curves show different sections of ω22(t) for an eccentric SEOBNRv4EHM waveform (with binary parameters shown in the title). If we discard all times below the point where the orbit-averaged frequency f22ω22/(2π) (pink dashed curve) crosses flow=20Hz, only the pink section is retained, and the blue section is discarded even though it contains some frequencies above 20Hz. On the other hand, using f22pω22p/(2π) (blue dashed curve) to pick this time ensures that the discarded region (orange) contains no frequencies above 20Hz.

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (18)

      Figure 4

      Limitations of the amplitude and frequency methods in identifying pericenters (blue circles) and apocenters (pink squares) for a low eccentricity waveform. These methods (top two rows) detect only the first few pericenters or apocenters and fail once sufficient eccentricity is radiated away. On the other hand, the residual amplitude and residual frequency methods (bottom two rows) can detect all of the pericenters or apocenters present. The waveform is generated using SEOBNRv4EHM, and the binary parameters are given in the title.

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (19)

      Figure 5

      Comparison of the amplitude (top) and the frequency (bottom) of an eccentric SEOBNRv4EHM waveform to those of its quasicircular counterpart. The binary parameters are shown in the figure text. Both waveforms are aligned so that t=0 occurs at the peak of A22. The quasicircular counterpart captures the secular growth in the amplitude and frequency of the eccentric waveform.

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (20)

      Figure 6

      Illustration of the frequency fits method. Left: the blue circles indicate the 2N+1=7 extrema through which the fitting function Eq.(30) passes. The lower panel shows the envelope-subtracted data from which the extrema Tα are determined. The solid blue circle indicates the central extremum, whose parameters are used for the eccentricity definition. The pink square and the pink dashed line show the analogous construction for the apocenter passages. Right: enlargement of the region around the solid markers in the upper panel on the left. The waveform is generated using SEOBNRv4EHM, and the binary parameters are given in the title.

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (21)

      Figure 7

      Comparison of egw and eω22 to the geodesic eccentricity egeo in the q limit, as a function of the orbit-averaged frequency ω22. In the left panel, the colors show the absolute difference between egeo and egw measured using Eq.(8) with the amplitude method. The right panel shows the same for eω22. egeo is closer to egw than eω22 by about 2 orders of magnitude.

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (22)

      Figure 8

      Demonstration of the measurement of eccentricity using the gw_eccentricity [93] package for waveforms of different origins: PN, EOB, NR, and EMRI. The binary parameters are indicated in the figure text. In each subplot, the lower panel shows the real part of h22, and the upper panel shows the measured eccentricity. We consider three different methods for identifying the pericenters and apocenters: amplitude, residual amplitude, and amplitude fits.

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (23)

      Figure 9

      egw vs eeob at the initial time, for SEOBNRv4EHM waveforms with varying eeob, but keeping the other binary parameters fixed (given in figure title). eeob is the model’s internal eccentricity, specified at t0=20000M. egw is evaluated at its first available time t^0. We consider three different methods for locating pericenters and apocenters: amplitude, residual amplitude, and amplitude fits. The amplitude method breaks down for small eccentricities (eeob103), while the residual amplitude and amplitude fits method follow the expected egw=eeob trend down to eeob=105.

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (24)

      Figure 10

      egw vs the internal definition of eccentricity, for waveforms of different origin, for equal-mass nonspinning binaries with varying eccentricity. For the NR waveforms (SpEC), we compute the internal eccentricity at t0=1500M after the start of the simulation, while for the rest we use t0=20000M before peak waveform amplitude. In both cases, t0^ is the first available time for egw(t). The inset shows the same but on a linear scale and focuses on the egw0.4 region.

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (25)

      Figure 11

      egw(t) for SEOBNRv4EHM waveforms with varying eeob but keeping the other binary parameters fixed (given in the figure title). The method used to locate pericenters and apocenters is indicated in the figure text. The colors indicate the value of eeob, defined at t0=20000M. The amplitude method breaks down for small eccentricities egw103102, especially as one approaches the merger. The residual amplitude and amplitude fits methods continue to compute the eccentricity until egw105. The features at egw105 arise from the waveform model itself (see Fig.12).

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (26)

      Figure 12

      Tracing the noisy features in Fig.11 to the behavior of the SEOBNRv4EHM model at small eccentricities. The top panel shows egw for the case with eeob=1.05×105 at t0=20000M, from the middle panel in Fig.11. The bottom panel shows the corresponding Δω22 [Eq.(26)], which helps highlight the modulations due to eccentricity. The drop in egw occurs at the same time as an abrupt drop in the eccentricity modulations in Δω22 that arises from a transition function applied to the dynamical variables entering the NQC corrections in SEOBNRv4EHM [46].

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (27)

      Figure 13

      Differences in egw(t) due to different methods used to locate pericenters and apocenters, for the same SEOBNRv4EHM waveforms as Fig.11. Top left: the curves show egw(t) obtained using the residual amplitude method with the quasicircular counterpart also obtained from SEOBNRv4EHM. The colors represent the absolute difference with respect to the egw(t) obtained using the amplitude fits method, and the gray region shows the parts where the second method fails to compute egw(t). Top right: the same, but now the colors show the difference with respect to the egw(t) obtained with residual amplitude method with the quasicircular counterpart obtained from the IMRPhenomT model. In both top panels, the different choices for locating pericenters and apocenters lead to broadly consistent results for egw(t), with the only notable differences occurring for (i)small eccentricities (egw5×103) and near the merger, where the SEOBNRv4EHM model also has known issues (see Fig.12), and (ii)large eccentricities (egw0.9), where locating apocenters is problematic. The bottom panels show the same as the top panels but when identifying the midpoints between pericenters as apocenters. This leads to more consistent results between different methods, and the largest differences in egw decrease by an order of magnitude.

      Reuse & Permissions

    • Defining eccentricity for gravitational wave astronomy (28)

      Figure 14

      A smoothness test for egw(t) at very high eccentricities. The curves show the time evolution of egw for SEOBNRv4EHM waveforms with initial eccentricities 0.9eeob0.999 at t0=5×106M. The colors represent eeob at t0. The binary parameters are shown in the figure title. We use the residual amplitude method to locate pericenters and identify the midpoints between pericenters as apocenters.

      Reuse & Permissions

    ×

    Defining eccentricity for gravitational wave astronomy (2024)
    Top Articles
    Latest Posts
    Article information

    Author: Patricia Veum II

    Last Updated:

    Views: 6224

    Rating: 4.3 / 5 (64 voted)

    Reviews: 87% of readers found this page helpful

    Author information

    Name: Patricia Veum II

    Birthday: 1994-12-16

    Address: 2064 Little Summit, Goldieton, MS 97651-0862

    Phone: +6873952696715

    Job: Principal Officer

    Hobby: Rafting, Cabaret, Candle making, Jigsaw puzzles, Inline skating, Magic, Graffiti

    Introduction: My name is Patricia Veum II, I am a vast, combative, smiling, famous, inexpensive, zealous, sparkling person who loves writing and wants to share my knowledge and understanding with you.