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
- Abstract
- Authors
- Article Text
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.
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)
- Research Areas
Gravitational waves
- Physical Systems
Classical black holes
- 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
Access Options
- Buy Article »
- Log in with individual APS Journal Account »
- Log in with a username/password provided by your institution »
- Get access through a U.S. public or high school library »
Article part of CHORUS
Accepted manuscript will be available starting 2 November 2024.Images
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 (upper panel) and frequency of the (2, 2) waveform mode (lower panel). and are interpolants through evaluated at the pericenters (blue circles) and apocenters (pink squares), respectively. Equation(8) is used to compute given and . Right: time evolution of the mean anomaly (upper panel) and (lower panel). The vertical dashed gray lines denote the pericenter times. grows linearly in time from 0 to between successive pericenters [Eq.(10)].
Figure 2
Different methods to construct a monotonically increasing frequency to replace , in order to set the reference frequency for eccentric binaries. We consider two different approaches: (i), the mean of and , and (ii), an interpolant through the orbit-averaged [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 as it captures the correct frequency scale in an orbit-averaged sense (Sec.2e).
Figure 3
How to truncate time-domain eccentric waveforms while retaining all frequencies above . The orange, blue, and pink curves show different sections of 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 (pink dashed curve) crosses , 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 (blue dashed curve) to pick this time ensures that the discarded region (orange) contains no frequencies above 20Hz.
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.
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 occurs at the peak of . The quasicircular counterpart captures the secular growth in the amplitude and frequency of the eccentric waveform.
Figure 6
Illustration of the frequency fits method. Left: the blue circles indicate the extrema through which the fitting function Eq.(30) passes. The lower panel shows the envelope-subtracted data from which the extrema 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.
Figure 7
Comparison of and to the geodesic eccentricity in the limit, as a function of the orbit-averaged frequency . In the left panel, the colors show the absolute difference between and measured using Eq.(8) with the amplitude method. The right panel shows the same for . is closer to than by about 2 orders of magnitude.
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 , 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.
Figure 9
vs at the initial time, for SEOBNRv4EHM waveforms with varying , but keeping the other binary parameters fixed (given in figure title). is the model’s internal eccentricity, specified at . is evaluated at its first available time . We consider three different methods for locating pericenters and apocenters: amplitude, residual amplitude, and amplitude fits. The amplitude method breaks down for small eccentricities (), while the residual amplitude and amplitude fits method follow the expected trend down to .
Figure 10
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 after the start of the simulation, while for the rest we use before peak waveform amplitude. In both cases, is the first available time for . The inset shows the same but on a linear scale and focuses on the region.
Figure 11
for SEOBNRv4EHM waveforms with varying 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 , defined at . The amplitude method breaks down for small eccentricities , especially as one approaches the merger. The residual amplitude and amplitude fits methods continue to compute the eccentricity until . The features at arise from the waveform model itself (see Fig.12).
Figure 12
Tracing the noisy features in Fig.11 to the behavior of the SEOBNRv4EHM model at small eccentricities. The top panel shows for the case with at , from the middle panel in Fig.11. The bottom panel shows the corresponding [Eq.(26)], which helps highlight the modulations due to eccentricity. The drop in occurs at the same time as an abrupt drop in the eccentricity modulations in that arises from a transition function applied to the dynamical variables entering the NQC corrections in SEOBNRv4EHM [46].
Figure 13
Differences in due to different methods used to locate pericenters and apocenters, for the same SEOBNRv4EHM waveforms as Fig.11. Top left: the curves show obtained using the residual amplitude method with the quasicircular counterpart also obtained from SEOBNRv4EHM. The colors represent the absolute difference with respect to the obtained using the amplitude fits method, and the gray region shows the parts where the second method fails to compute . Top right: the same, but now the colors show the difference with respect to the 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 , with the only notable differences occurring for (i)small eccentricities () and near the merger, where the SEOBNRv4EHM model also has known issues (see Fig.12), and (ii)large eccentricities (), 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 decrease by an order of magnitude.
Figure 14
A smoothness test for at very high eccentricities. The curves show the time evolution of for SEOBNRv4EHM waveforms with initial eccentricities at . The colors represent at . 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.