---
title: FastEMRIWaveforms (FEW) Overview
url: https://www.emergentmind.com/topics/fastemriwaveforms-few
type: topic
---

# FastEMRIWaveforms (FEW) Overview

FastEMRIWaveforms (FEW) is a modular, open-source software package for generating and analyzing gravitational waveforms from Extreme Mass Ratio Inspiral (EMRI) systems, with an initial focus on LISA-band sources [2104.04582]. It was introduced as a collection of tools to build and analyze extreme mass ratio inspiral waveforms, and it was expanded from the Physical Review Letter that introduced the first fast and accurate fully-relativistic EMRI waveform template model [2104.04582]. In its original formulation, the fully relativistic FEW model waveforms are for now limited to eccentric orbits in the Schwarzschild spacetime, while later work extends FEW to model eccentric equatorial inspirals into black holes with spin magnitudes $|a| \leq 0.999$ [2104.04582] [2506.09470]. Across these implementations, the framework is characterized by separating computationally expensive offline precomputations from rapid online waveform generation via interpolation, mode selection, and hardware-accelerated summation [2506.09470].

## 1. Origins, scope, and constituent modules

FEW is described as a highly modular design in which components can be used standalone or together, and the code structure is layered for flexibility and easy updating as waveform models improve [2104.04582]. The main workflow modules are the Trajectory Module, the Amplitude Module, the Waveform Summation Module, and Utility Modules [2104.04582]. The Trajectory Module computes the orbit-averaged motion using suitable ODE integration; the Amplitude Module provides for rapid prediction of all necessary harmonic mode amplitudes via reduced-order/deep-learning surrogates; the Waveform Summation Module efficiently sums all relevant harmonic modes across the EMRI trajectory; and the Utility Modules include tools for mode selection, separatrix finding, harmonic generation, frequency calculation, and mismatch analysis [2104.04582].

The package includes a Python API that allows users to easily switch between CPU and GPU computation by changing a single argument [2104.04582]. For all computationally intensive routines, FEW uses C++/CUDA backends with smart device-agnostic wrappers for NumPy/CuPy in Python [2104.04582]. This modular structure was later complemented by a frequency-domain implementation integrated into the FastEMRIWaveforms package for use with either graphics processing units (GPUs) or central processing units (CPUs) [2307.12585].

A plausible implication is that FEW was designed not as a single waveform approximant, but as an extensible analysis environment in which trajectory evolution, amplitude generation, and waveform synthesis can be updated independently as new self-force, Teukolsky, or surrogate ingredients become available.

## 2. Computational methodology and acceleration strategy

The core methodology relies on separating computationally expensive “offline” precomputations of physical data from the rapid “online” generation of specific waveforms via interpolation and integration [2506.09470]. In the offline phase, self-force data and spherical decomposition of mode amplitudes in the relevant parameter space are computed in advance; these involve expensive black hole perturbation theory or Teukolsky equation computations, but only need to be done once and can be parallelized [2506.09470]. These data are tabulated as dense multi-dimensional grids over parameters $(a, p, e)$, carefully parameterized to maximize interpolation accuracy, with specialized tapers and grids for efficiency [2506.09470].

In the online phase, the system’s adiabatic trajectory is computed by interpolating the stored fluxes to drive coupled ODEs for orbital elements $(p, e, x_I)$ and phases, while mode amplitudes are rapidly interpolated [2506.09470]. Dense and smooth variations allow the use of spline interpolators, with tricubic for fluxes and bicubic+linear for amplitudes, leading to sub-millisecond evaluations per grid point [2506.09470]. Trajectory evolution, mode spectra, and waveform summation are then handled online, exploiting hardware acceleration for both interpolation and summation across many harmonics [2506.09470].

Several acceleration mechanisms recur across FEW implementations. Sparse trajectory sampling exploits the slow evolution of inspirals, so that the ODE integrating the trajectory outputs only $\sim 100$ samples, which are efficiently upsampled by cubic spline during waveform synthesis [2104.04582]. Efficient mode selection picks only modes contributing above a user threshold to total power, and symmetry relations imply that only $m \geq 0$ amplitudes need be directly computed [2104.04582]. For the GPU implementation, spline interpolation and the large sum to construct the waveform are implemented as batched, parallelized GPU kernels, and the GPU summation is $\sim$2500x–3000x faster than CPU for waveform summation [2104.04582]. In the spinning equatorial extension, GPU support on an NVIDIA A100 or similar gives more than 1000x speedup vs. CPU for large waveforms, and typical one-year FEW waveforms with thousands of harmonics can be generated in $\sim 100$ ms [2506.09470].

This suggests that FEW’s speed derives less from a single surrogate component than from a pipeline architecture in which precomputation, interpolation, sparse trajectory sampling, mode truncation, and parallel reduction are all treated as first-class algorithmic objects.

## 3. Waveform representations and model families

The fully relativistic model for Schwarzschild eccentric orbits uses first-order gravitational self-force results (adiabatic inspiral) for motion and full Teukolsky-based amplitude computations for wave generation, restricted to Schwarzschild spacetime and thus to equatorial, eccentric orbits [2104.04582]. In that case the waveform may be written as
\[
h = \frac{\mu}{d_L}\sum_{lmn} A_{lmn}(t)\,_{-2}Y_{lm}(\theta, \phi)e^{-i\Phi_{mn}(t)},
\]
with phase accumulation
\[
\Phi_{\varphi,\theta,r}(t) = \int_0^t \Omega_{\varphi, \theta, r}\big(p(t'), e(t'), x_I(t')\big) dt'
\]
[2104.04582].

The improved Augmented Analytic Kludge (AAK) model is a semi-relativistic kludge model designed for speed and flexibility, applicable to generic Kerr and generic inspirals [2104.04582]. Its trajectory module uses robust 5PN-order inspiral ODEs for $(p, e, Y=\cos\iota)$, the frequency evolution uses actual Kerr fundamental frequencies $(\Omega_\varphi, \Omega_\theta, \Omega_r)$, and the output is built in the LISA detector frame [2104.04582]. Within FEW, the relativistic Schwarzschild model and the improved AAK model coexist as distinct waveform families with different domains of validity [2104.04582].

A separate development introduced 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) [2307.12585]. In that implementation, the waveform is decomposed voice-by-voice,
\[
h(t) = \sum_{V} H_V(t) e^{-i\Phi_V(t)},
\]
and the Fourier transform is generated using the stationary phase approximation [2307.12585]. On GPUs, the frequency domain waveform takes in median $0.044$ seconds to generate and is twice as fast to compute as its time domain counterpart when considering massive black hole masses $\geq 2 \times 10^6 \,{\rm M_\odot}$ and initial eccentricities $e_0 > 0.2$; on CPUs, the median waveform evaluation time is $5$ seconds, and it is five times faster in the frequency domain than in the time domain [2307.12585].

A plausible implication is that FEW should be understood as a family of interoperable waveform-generation strategies rather than a single domain-specific approximant: time-domain and frequency-domain synthesis, fully relativistic and kludge amplitudes, and surrogate-augmented modules all appear within the same package architecture.

## 4. Extension to spinning eccentric equatorial inspirals

The 2025 extension addresses eccentric equatorial inspirals into black holes with spin magnitudes $|a| \leq 0.999$ [2506.09470]. The model supports eccentricities $e < 0.9$ and semi-latus recta $p < 200$, and includes trajectories very close to the separatrix with $p - p_\text{sep} \geq 10^{-3}$ [2506.09470]. It is designed for $q = m_1/m_2 \gtrsim 10^2$ but can be used for moderately lower $q$, and it covers both prograde and retrograde spins, even beyond the so-called Thorne limit [2506.09470]. The type of orbits is eccentric equatorial inspirals, with the framework extendable to fully generic orbits once self-force data is available [2506.09470].

For equatorial inspirals $(k=0)$, the model is built around the multivoice sum
\[
h(t) = \frac{\mu}{d_\mathrm{L} \sum_{\ell, m, n} H_{\ell m n}(t, \theta, \phi) e^{-i \Phi_{mn}(t)},
\]
with
\[
\Phi_{mn}(t) = m\Phi_\phi(t) + n\Phi_r(t),
\]
\[
H_{\ell m n}(t) = A_{\ell m n}(t) {}_{-2}S_{\ell m n}(\theta; \hat\omega_{mn}) e^{i m \phi},
\]
\[
A_{\ell m n}(t) = -2 \frac{\hat Z_{\ell m n}^\infty(a, p, e, x_I)}{[\hat\omega_{mn}(a, p, e, x_I)]^2},
\]
and
\[
\hat\omega_{mn} = m\hat\Omega_\phi + n\hat\Omega_r
\]
[2506.09470]. The adiabatic evolution of orbital elements is written as
\[
\frac{d\alpha}{dt} = \frac{\nu}{M} \left[ \hat f^{(0)}_\alpha(a,p,e,x_I) + O(\nu) \right], \quad \alpha\in\{p, e\},
\]
\[
\frac{d\Phi_A}{dt} = \frac{1}{M}\hat\Omega_A(a,p,e,x_I), \quad A\in\{\phi, r\}
\]
[2506.09470].

Among the significant advancements listed for this version are extension to high spin and broad parameter space; fully relativistic amplitudes and harmonics incorporating thousands of $(\ell, m, n)$ modes up to $\ell=10, n=55$; advanced data gridding and interpolation; a Python-based, hardware-accelerated, modular codebase; robust validation; practical data analysis demonstrators; and immediate extensibility to inclined or generic orbits, merger-ringdown, post-adiabatic corrections, and environmental effects [2506.09470].

## 5. Accuracy, mode content, and biases from reduced models

A persistent result across FEW studies is that only a modest number of harmonics is often required if those harmonics are chosen by cumulative power rather than by a fixed low-multipole truncation. In the 2021 package study, relativistic EMRI waveform templates can be generated with fewer harmonic modes $(\sim10$–$100)$ without biasing signal extraction [2104.04582]. The on-the-fly mode selector keeps all modes whose cumulative power exceeds $1-\epsilon$, and only $\sim$10–100 harmonic modes are needed to capture $>99.9\%$ of emitted power for the purpose of data analysis [2104.04582]. In the spinning equatorial extension, mode selection is based on cumulative amplitude-power up to a tolerance $\kappa$ [2506.09470].

The same studies emphasize that reduced or kludge amplitude content is not interchangeable with full relativistic harmonic content. In the 2021 analysis, extraction of a relativistic injection with semi-relativistic amplitudes can lead to strong bias and anomalous structure in the posterior distribution for certain regions of parameter space [2104.04582]. Relativistic, quadrupole-only templates have mismatch of order $0.05$–$0.1$, while semi-relativistic amplitudes can have mismatch up to $0.6$ at high eccentricity [2104.04582]. The role of phase vs. amplitude is also explicitly differentiated: matching the phase is far more important than matching the amplitude for parameter estimation, while amplitude errors mostly manifest as mild biases unless gross [2104.04582].

In the spinning equatorial FEW model, the adiabatic error over most of parameter space is characterized by mismatches of $\sim 10^{-5}$ with respect to error-free adiabatic models, while errors grow in regions of very high eccentricity $(e \gtrsim 0.85)$ or high spin near the separatrix, but remain $\lesssim 10^{-3}$ for most cases [2506.09470]. By contrast, kludge models can overestimate the SNR by $+60\%$ or underestimate by $-40\%$, and induce marginal biases of up to $\sim 1\sigma$ in parameter recovery [2506.09470]. Using only the dominant $\ell=2$ multipoles induces similar biases, with inaccuracies scaling with spin and eccentricity [2506.09470].

A common misconception is that a few harmonics are sufficient because the signal is effectively quadrupolar. The FEW results support a more specific statement: a few harmonics are often sufficient only when they are selected to cover the significant power in the signal; quadrupole-only or semi-relativistic truncations can be insufficient for precise source characterization [2104.04582].

## 6. Data-analysis applications, astrophysical reach, and related extensions

FEW was developed for millihertz gravitational-wave data analysis, and with the GPU-accelerated waveforms in hand, a variety of studies were performed including an analysis of EMRI mode content, template mismatch, and fully Bayesian Markov Chain Monte Carlo-based EMRI parameter estimation [2104.04582]. The 2021 work states that fully Bayesian posteriors can be computed for two-year duration EMRI signals with fully relativistic templates, thanks to the speed of the FEW GPU implementation [2104.04582]. The 2023 frequency-domain implementation further states that this enables, for the first time, EMRI parameter inference with fully relativistic waveforms on CPUs [2307.12585].

The spinning equatorial extension quantifies the astrophysical reach of FEW-based I/EMRI templates for LISA. For a detection threshold SNR$=20$ and for high spin, $a=0.999$, the accurate waveform model gives EMRIs $(q \sim 10^5)$ to redshift $z \sim 3$ and IMRIs $(q \sim 10^2)$ to redshift up to $z \sim 15$, with only mildly dependent $(\sim 10\%)$ variation on eccentricity [2506.09470]. For SNR $\sim 50$, $\delta a \sim 10^{-7}$ and $\delta e_\mathrm{f} \sim 10^{-5}$; source-frame masses are measured to $\sim 10^{-3}$, luminosity distance to $\sim 1\%$, and 99% sky localization area down to $\sim 1\,\mathrm{deg}^2$ for EMRIs [2506.09470]. Full 12D parameter estimation confirms near-Gaussian posteriors and minimal parameter degeneracies for EMRIs, whereas IMRIs require including many harmonics because biases become significant if too few modes are kept [2506.09470].

A related extension, "A Deep Learning Framework for Amplitude Generation of Generic EMRIs" [2603.08635], addresses the major challenge for the fast generation of such waveforms: the generation of the Teukolsky amplitudes for generic Kerr orbits. It introduces a convolutional encoder-decoder architecture for a fast and end-to-end global fitting of the Teukolsky amplitudes, with transfer learning from Schwarzschild circular orbits to generic Kerr orbits [2603.08635]. Within this framework, the full harmonic amplitudes can be generated within milliseconds, while the median mode-distribution error for generic orbits is approximately $\sim10^{-3}$, and the surrogate can be used within adiabatic waveform frameworks such as FEW [2603.08635]. This suggests a pathway from the eccentric equatorial Kerr extension of FEW toward generic $(e,\iota)$ Kerr waveform generation once suitable amplitude datasets are available.

In aggregate, FEW occupies a specific place in EMRI and IMRI waveform modeling: it combines offline black hole perturbation theory or surrogate construction with online interpolated adiabatic evolution, mode selection, and hardware-accelerated summation, and it is used to test when reduced models remain adequate and when they induce detectable bias [2104.04582] [2506.09470].

Source: https://www.emergentmind.com/topics/fastemriwaveforms-few