Fast and Fourier: extreme mass ratio inspiral waveforms in the frequency domain

• Published in: Frontiers in applied mathematics and statistics Vol. 9
• Main Authors: Speri, Lorenzo, Katz, Michael L., Chua, Alvin J. K., Hughes, Scott A., Warburton, Niels, Thompson, Jonathan E., Chapman-Bird, Christian E. A., Gair, Jonathan R.
• Format: Journal Article
• Language: English
• Published: Frontiers Media S.A 12-01-2024
• Subjects: data analysis; EMRI; frequency domain; gravitational wave; LISA
• ISSN: 2297-4687; 2297-4687
• DOI: 10.3389/fams.2023.1266739

Extreme Mass Ratio Inspirals (EMRIs) are one of the key sources for future space-based gravitational wave interferometers. Measurements of EMRI gravitational waves are expected to determine the characteristics of their sources with sub-percent precision. However, their waveform generation is challenging due to the long duration of the signal and the high harmonic content. Here, we present the first ready-to-use Schwarzschild eccentric EMRI waveform implementation in the frequency domain for use with either graphics processing units (GPUs) or central processing units (CPUs). We present the overall waveform implementation and test the accuracy and performance of the frequency domain waveforms against the time domain implementation. On GPUs, the frequency domain waveform takes in median 0.044 s to generate and is twice as fast to compute as its time domain counterpart when considering massive black hole masses
≥
2
×
1
0
6
M
⊙
and initial eccentricities
e
0
> 0.2. On CPUs, the median waveform evaluation time is 5 s, and it is five times faster in the frequency domain than in the time domain. Using a sparser frequency array can further speed up the waveform generation, reaching up to 0.3 s. This enables us to perform, for the first time, EMRI parameter inference with fully relativistic waveforms on CPUs. Future EMRI models, which encompass wider source characteristics (particularly black hole spin and generic orbit geometries), will require significantly more harmonics. Frequency domain models will be essential analysis tools for these astrophysically realistic and important signals.