---
title: TEOBResumSPA_NRPMw Waveform Model
url: https://www.emergentmind.com/topics/teobresumspa_nrpmw-waveform
type: topic
---

# TEOBResumSPA_NRPMw Waveform Model

Searching arXiv for TEOBResumSPA, NRPMw, and closely related TEOBResumS papers to ground the article in the relevant literature.
The **TEOBResumSPA\_NRPMw waveform** is a complete inspiral–merger–postmerger waveform model for binary neutron star systems that combines **TEOBResumSPA**—or its time-domain parent **TEOBResumD**—for the inspiral and merger up to the time-domain peak of the \((2,2)\) amplitude, with **NRPMw**, a numerical-relativity–informed **frequency-domain** model of the postmerger kilohertz signal from the remnant [2205.09112]. In this construction, the inspiral–merger sector inherits the effective-one-body dynamics, resummed multipolar amplitudes, NR-informed next-to-quasi-circular structure, and acceleration strategies developed in the TEOBResumS family [1805.03891], [2104.07533], while the postmerger sector is modeled analytically as a sum of complex Gaussian wavelets calibrated to numerical-relativity simulations through EOS-insensitive relations [2205.09112]. The result is a single waveform family intended to cover the signal from low frequencies through merger and postmerger out to several kilohertz [2205.09112].

## 1. Definition and signal domain

TEOBResumSPA\_NRPMw targets the gravitational-wave signal of binary neutron star remnants **after** the amplitude peak of merger, while retaining a consistent inspiral–merger description from the TEOBResumS/TEOBResumSPA framework [2205.09112]. The postmerger regime modeled by NRPMw is the dominant \((\ell,m)=(2,2)\) multipole emitted from the time when the two neutron-star cores have fused and a massive, hot, differentially rotating remnant has formed, up to the point where the remnant collapses to a black hole or relaxes to a long-lived neutron star on secular timescales [2205.09112].

In this regime, the signal is dominated by quasi-periodic oscillations at a nearly constant frequency \(f_2 \sim 2\!-\!4\) kHz, shows broad spectral peaks rather than discrete lines, with full-width-at-half-maximum typically \(300{-}600\) Hz, exhibits strong amplitude and frequency modulations associated with quasi-radial oscillations of the remnant at a characteristic frequency \(f_0 \sim \mathcal{O}(1\,{\rm kHz})\), and may end abruptly when the remnant collapses to a black hole [2205.09112]. TEOBResumSPA provides \(\tilde h_{\rm insp}(f)\) up to the merger frequency \(f_{\rm mrg}\), whereas NRPMw supplies \(\tilde h_{\rm PM}(f)\) in the kilohertz range, typically from \(\sim 1\) kHz up to \(\sim 8\) kHz, the band used in NR faithfulness tests with ET-D noise [2205.09112].

The inspiral–merger component is rooted in the TEOBResumS family of effective-one-body models. In that framework, the waveform is constructed mode by mode from resummed EOB dynamics, a multipolar flux, NR-informed NQC corrections, and an NR-informed merger–ringdown attachment [2104.07533]. The TEOBResumSPA variant is the stationary-phase–approximant representation of the same underlying EOB dynamics, and the post-adiabatic approximation provides an efficient route to compute inspiral waveforms up to a few orbits before merger with phase differences below \(10^{-3}\) rad in representative BBH and BNS cases [1805.03891].

## 2. Frequency-domain postmerger construction

NRPMw models the postmerger waveform as a sum of complex Gaussian wavelets with amplitude and frequency modulations [2205.09112]. Its basic time-domain building block is
\[
W(t;\alpha,\beta,\gamma,\tau) = \begin{cases} e^{\alpha t^2+\beta t +\gamma} & \text{if } t\in[0,\tau]\,,\\[2pt] 0 & \text{otherwise}\,, \end{cases}
\]
with complex coefficients \(\alpha,\beta,\gamma\in\mathbb{C}\) and support \([0,\tau]\) [2205.09112]. The Fourier transform can be written analytically as
\[
W(f) = \frac{e^{\gamma}}{2} \sqrt{\frac{\pi}{\alpha}}\, e^{-z^2}\, \left[{\rm erfi}{\left( z + \sqrt{\alpha} \tau \right)}- {\rm erfi}{\left( z \right)} \right],
\]
where
\[
z(f) = \frac{\beta-2\pi i f}{2\sqrt{\alpha}},
\]
and a global time shift is implemented through
\[
W(f;t_0) = W(f)\,e^{-2\pi i f t_0}.
\]
For \(\alpha=0\), the wavelet reduces to a damped sinusoid and \(W(f)\) becomes a Lorentzian [2205.09112].

Amplitude modulations are introduced through
\[
\tilde W_{\rm am}(t) = W(t)\left[ 1 + \Delta_{\rm am}  \sin\left( \Omega_{\rm am} t +\phi_{\rm am} \right) \right],
\]
which in frequency domain yields sidebands at \(f_2\pm f_0\), corresponding to the secondary spectral peaks \(f_{2\pm 0}=f_2\pm f_0\) [2205.09112]. Frequency modulations act on the instantaneous angular frequency,
\[
\omega_W(t) = -2 \,\Im(\alpha) t - \Im(\beta),
\]
through
\[
\omega_{\rm FM}(t) = \omega_W(t) - \Delta_{\rm fm} e^{-\Gamma_{\rm fm} t}\sin(\Omega_{\rm fm} t + \phi_{\rm fm}),
\]
with \(\Omega_{\rm fm}=2\pi f_0\), yielding an analytic frequency-domain representation with multiple sidebands at \(f_2 \pm n f_0\) [2205.09112].

The dominant \((2,2)\) postmerger mode is represented as
\[
h_{22}(t) \approx W_{\rm fus}(t) + W_{\rm bnc}(t) + W_{\rm pul}(t) + W_{\rm peak}(t),
\]
where the four components are physically motivated [2205.09112]. \(W_{\rm fus}\) models the early postmerger fusion stage from merger time \(t_{\rm mrg}\) to the first amplitude minimum at \(t_0\); \(W_{\rm bnc}\) models the first bounce of the merged core from \(t_0\) to the first postmerger amplitude maximum \(t_1\); \(W_{\rm pul}\) describes the strongly dynamical pulsating phase from \(t_1\) to \(t_3\) with both AM and FM active; and \(W_{\rm peak}\) models the late damped tail from \(t_3\) to the collapse time \(t_{\rm coll}\) or to a large cutoff for long-lived remnants [2205.09112]. A fifth piece \(W_{\rm coll}\) could represent black-hole ringdown after collapse, but it is set to zero in this work because its power is negligible in the 1–8 kHz band for current/ET sensitivities [2205.09112].

The corresponding frequency-domain waveform is
\[
\tilde h_{22}(f) = \tilde W_{\rm fus}(f) + \tilde W_{\rm bnc}(f) + \tilde W_{\rm pul}(f) + \tilde W_{\rm peak}(f),
\]
with each term an analytical function of the wavelet and modulation parameters [2205.09112]. The dominant peak at \(f_2\) corresponds to the main quadrupole oscillation of the remnant, while sidebands at \(f_{2-0}\) and \(f_{2+0}\) arise from AM/FM coupling with the quasi-radial mode at \(f_0\) [2205.09112].

## 3. Calibration dataset and EOS-insensitive mapping

NRPMw is calibrated on **618 numerical-relativity simulations** spanning total mass \(M\in[2.4,3.4]\,M_\odot\), mass ratio \(q\in[1,2]\), symmetric tidal coupling parameter \(\kappa_2^T\in[22,458]\), aligned spins \(|\chi_i|\lesssim 0.2\), and 21 EOSs, including 7 finite-temperature microphysical EOSs and 14 piecewise-polytropic EOSs [2205.09112]. Some EOSs include quark deconfinement or hyperons, and the dataset is constructed from three independent NR codes, {\sc BAM}, {\sc THC}, and the Kyoto/Tokyo code [2205.09112].

The binary parameter vector used to map inspiral information into postmerger quantities is
\[
\theta_{\rm bin} = \{m_1, m_2, \Lambda_1, \Lambda_2, \chi_1, \chi_2\},
\]
with tidal coupling constant
\[
\kappa_2^{\rm T} = 3\nu\,\left[\left(\frac{m_1}{M}\right)^3 \Lambda_1 + \left(\frac{m_2}{M}\right)^3 \Lambda_2\right],
\qquad
\nu = \frac{m_1 m_2}{M^2}
\]
[2205.09112]. For any postmerger quantity \(Q\), the fit ansatz is
\[
Q^{\rm fit} = a_0 \,Q^{\rm M}(X)\,Q^{\rm S}(\hat S,X)\,Q^{\rm T}(\kappa_2^{\rm T}, X),
\]
where
\[
X = 1-4\nu = \frac{(q-1)^2}{(q+1)^2},
\qquad
\hat S  = \left(\frac{m_1}{M}\right)^2 \chi_1 + \left(\frac{m_2}{M}\right)^2  \chi_2,
\]
and the tidal factor is the rational function
\[
Q^{\rm T}  =\frac{1+p_1^{\rm T} \,\kappa_2^{\rm T} + p_2^{\rm T}\, (\kappa_2^{\rm T})^2}{1+p_3^{\rm T}\, \kappa_2^{\rm T} + p_4^{\rm T} \,(\kappa_2^{\rm T})^2}
\]
[2205.09112].

The postmerger parameter vector is
\[
\theta_{\rm PM} = \{\phi_{\rm PM},\phi_{\rm fm},t_0, t_{\rm coll}, A_{\rm mrg}, A_0, A_1, A_2, A_3, f_{\rm mrg}, f_2, f_0, \Delta_{\rm fm}, \Gamma_{\rm fm}, \delta f_{\rm peak}, \delta f_{\rm fus} \},
\]
and the subset mapped from inspiral parameters is
\[
\theta_{\rm map} = \{A_{\rm mrg}, A_0, A_1, A_2, A_3, f_{\rm mrg}, f_2, f_0, t_0, \delta f_{\rm peak}, \delta f_{\rm fus}, \Delta_{\rm fm}, \Gamma_{\rm fm}\}
\]
[2205.09112]. The remaining parameters \(\{\phi_{\rm PM},t_{\rm coll},\phi_{\rm fm}\}\) are treated as free or fitted directly to data [2205.09112].

Several relations are comparatively robust. The merger amplitude \(A_{\rm mrg}/M \in [0.159,0.313]\) has a \(1\sigma\) error of about 1.8%, the merger frequency \(M f_{\rm mrg}/\nu\in[0.0554,0.141]\) about 2.6%, and the main postmerger peak \(M f_2 \in [0.0216,0.0512]\) about 3.9%, corresponding to roughly 100 Hz near 2.5 kHz [2205.09112]. By contrast, \(f_0\), the postmerger nodal amplitudes, and the FM parameters \(\Delta_{\rm fm},\Gamma_{\rm fm}\) are substantially less constrained, with the FM parameters exhibiting errors of \(\sim 70{-}100\%\) [2205.09112]. EOSs with strong phase transitions can mildly violate the \(f_2(\kappa_2^T)\) relation at the \(\gtrsim 2\sigma\) level in special cases, though the relation remains overall robust [2205.09112].

## 4. Coupling to TEOBResumSPA and EOB foundations

The inspiral–merger sector attached to NRPMw inherits the effective-one-body dynamics and waveform factorization of the TEOBResumS family [1805.03891], [2104.07533]. In the nonspinning multipolar construction, the EOB Hamiltonian is
\[
\hat{H}_{\rm EOB} = \frac{1}{\nu} \sqrt{1 + 2 \nu \left( \hat{H}_{\rm eff} - 1\right)},
\]
with
\[
\hat{H}_{\rm eff} = \sqrt{ p_{r_*}^2 + A(r,\nu) \left( 1 + \frac{p_\varphi^2}{r^2} + \mathcal{Q}_4 \right) },
\qquad
\mathcal{Q}_4 = 2 \nu (4 - 3 \nu)\, p_{r_*}^4
\]
[1904.09550]. The radial potential is written as a PN series augmented by an effective 5PN coefficient \(a_6^c(\nu)\), then Padé-resummed as
\[
A(u;\nu;a_6^c) = P^1_5\big[ A^{\rm PN}_{\rm orb}(u)\big]
\]
in the nonspinning multipolar model [1904.09550]. The calibrated NR-informed fit for the Padé-resummed, 6PN-hybrid amplitude model is
\[
a_6^c(\nu) = n_0\,\frac{1+n_1 \nu + n_2 \nu^2 + n_3 \nu^3}{1+d_1\nu}
\]
with numerical coefficients reported in the paper [1904.09550].

The inspiral waveform is decomposed into spherical-harmonic multipoles,
\[
h_+ - i h_\times = \frac{GM}{\mathcal{D}_L c^2} \sum_{\ell=2}^\infty\sum_{m=-\ell}^{\ell} h_{\ell m}(t)\,{}_{-2}Y_{\ell m},
\]
and each orbital mode is factorized as
\[
h_{\ell m} = h_{\ell m}^{(N,\epsilon)}\, \hat{S}_{\rm eff}^{(\epsilon)}\, \hat{h}_{\ell m}^{\rm tail}(y)\, \big[\rho_{\ell m}^{\rm orb}(x)\big]^\ell\, \hat{h}_{\ell m}^{\rm NQC}
\]
[1904.09550]. The improved nonspinning multipolar model upgrades the residual amplitude corrections \(\rho_{\ell m}^{\rm orb}\) by hybridizing full 3PN \(\nu\)-dependent information with test-mass results up to 6PN order and then Padé-resumming them [1904.09550]. The calibrated model includes NR-completed higher modes for \((2,1)\), \((3,3)\), \((3,2)\), \((3,1)\), \((4,4)\), \((4,3)\), \((4,2)\), \((4,1)\), and \((5,5)\), while other modes up to \(\ell=8\) remain purely analytical [1904.09550].

The TEOBResumS family also uses NR-informed NQC corrections and can be accelerated by a post-adiabatic treatment of the inspiral. In that approximation, the EOB momenta are solved algebraically on a radial grid and the time and phase are recovered by quadratures, yielding waveforms more than 100 times faster than a standard ODE solver in a nonoptimized Matlab implementation for a standard BNS system from 10 Hz [1805.03891]. A later assessment of consistent NQC corrections and PA approximation in multipolar BBH waveforms showed EOB/NR unfaithfulness well below \(0.01\) over 611 NR simulations, with 78.5% of cases below \(0.001\) for \(10M_\odot\leq M \leq 200M_\odot\) using Advanced LIGO noise [2104.07533]. This suggests that the TEOBResumSPA inspiral attached to NRPMw is designed to inherit an NR-calibrated EOB backbone rather than a purely PN one.

## 5. Joining prescription, recalibration, and parameter inference

TEOBResumSPA\_NRPMw is assembled by joining a frequency-domain inspiral–merger waveform \(\tilde h_{\rm IM}(f)\) to a frequency-domain remnant waveform \(\tilde h_{\rm PM}(f)\) near the merger frequency \(f_{\rm mrg}\) [2205.09112]. In practice, one computes \(f_{\rm mrg}\) and \(A_{\rm mrg}\) from the NRPMw quasiuniversal relations, chooses a joining frequency \(f_{\rm join}\) close to \(f_{\rm mrg}\), fixes \(A_{\rm mrg}\) and \(\phi_{\rm mrg}\) in the fusion wavelet to match TEOBResumSPA at merger, and sets the time origin at the TEOBResumSPA merger time so that postmerger nodal times are referenced consistently [2205.09112].

To avoid sharp spectral features, the transition is smoothed by tapering or blending:
\[
\tilde h(f) = [1-w(f)]\,\tilde h_{\rm IM}(f)+w(f)\,\tilde h_{\rm PM}(f),
\]
where \(w(f)\) is a Planck-like or raised-cosine window over a small interval around \(f_{\rm join}\) [2205.09112]. Because the fusion wavelet is parameterized by \(A_{\rm mrg}\), \(f_{\rm mrg}\), and \(\phi_{\rm mrg}\), the full waveform is phase-coherent across inspiral, merger, and postmerger [2205.09112].

A central feature of NRPMw is its use of **recalibration parameters** \(\delta_Q\) to represent the residual scatter of the EOS-insensitive relations [2205.09112]. For any calibrated quantity,
\[
Q = Q^{\rm fit}(1+\delta_Q),
\]
with \(\delta_Q\) treated as a dimensionless fractional correction [2205.09112]. In Bayesian analyses, \(\delta_Q\) is assigned a Gaussian prior with mean zero and variance equal to the variance of the NR residuals. This enables marginalization over theoretical modeling errors and allows the waveform to behave as an informed but flexible postmerger model [2205.09112]. The same framework also permits tests of quasiuniversality: if data require large \(\delta_{f_2}\) inconsistent with the priors, that may indicate an EOS beyond the calibration assumptions, such as a strong phase transition [2205.09112].

In a parameter-estimation pipeline such as the `bajes` code mentioned in the paper, one samples
\[
\theta_{\rm bin}=\{m_1,m_2,\chi_1,\chi_2,\Lambda_1,\Lambda_2\},
\]
imposes inspiral–postmerger consistency through \(\kappa_2^T\), chooses priors for the free postmerger parameters \(\phi_{\rm PM}\) and \(t_{\rm coll}\), includes \(\{\delta_Q\}\) with Gaussian priors, and evaluates the likelihood over the full frequency band, for example 10 Hz–8 kHz [2205.09112]. The output constrains masses, spins, tidal parameters, postmerger peak frequencies, collapse time, and possible deviations from EOS-insensitive relations [2205.09112].

## 6. Validation, scope, and limitations

NRPMw was validated on **102 independent NR simulations** not used in calibration, using ET-D noise in the frequency interval \(f\in[1,8]\) kHz [2205.09112]. The faithfulness is defined by maximizing the overlap over overall time and phase shifts,
\[
\mathcal{F} = \max_{t_{\rm mrg},\phi_{\rm mrg}} \frac{(h_1|h_2)}{\sqrt{(h_1|h_1)(h_2|h_2)}},
\qquad
\bar{\mathcal{F}} = 1-\mathcal{F}
\]
[2205.09112]. Without recalibration, and optimizing only over \(\{\phi_{\rm PM},t_{\rm coll},\delta f_{\rm peak}\}\), the median unfaithfulness is \(\sim 0.13\), with about 38% of cases below 0.1 [2205.09112]. When all recalibration parameters are allowed to vary within \(\pm 4\sigma\) of their NR-based priors, the median unfaithfulness improves to about \(2.5\times 10^{-2}\), about 94% of the validation sample lies below 0.1, and many long-lived and unequal-mass cases reach \(\sim 10^{-2}\) [2205.09112].

The validated regime is explicitly limited to
\(M\in[2.4,3.4]\,M_\odot\),
\(q\in[1,2]\),
\(\kappa_2^T\in[22,458]\),
\(|\chi_i|\lesssim 0.2\),
and EOSs similar to those used in calibration [2205.09112]. Extrapolation may be unreliable for very high total masses near prompt collapse, extreme mass ratios \(q\gg 2\), large spins \(|\chi|\gtrsim 0.2\), or EOSs with strong, sharp phase transitions not represented in the training set [2205.09112]. The FM parameters are particularly uncertain and may need to be treated as free or strongly recalibrated parameters, as they can encode effects such as turbulent viscosity [2205.09112].

The inspiral–merger foundation also has model-dependent systematics. In the TEOBResumS multipolar EOB literature, explicit mode mixing is neglected in the ringdown treatment of modes such as \((2,1)\), \((3,2)\), \((4,3)\), and \((4,2)\), and this simplification can mildly affect unfaithfulness while keeping it below the percent-to-few-percent level in the tested BBH ranges [1904.09550]. Later analyses of analytic systematics in TEOBResumS emphasized that correct control of the noncircular part of the dynamics during late plunge is essential, and that simply adding higher-PN information without retuning can worsen EOB/NR agreement [2304.09662]. This suggests that, in TEOBResumSPA\_NRPMw, the inspiral and postmerger parts should be understood as jointly calibrated components rather than separable approximations.

## 7. Scientific role and extensions within the TEOBResumS family

Within the TEOBResumS family, TEOBResumSPA\_NRPMw occupies the role of a frequency-domain, NR-informed waveform that links a calibrated EOB inspiral to a phenomenological but physics-motivated BNS postmerger [2205.09112], [2503.14580]. The broader TEOBResumS program has expanded from quasi-circular spin-aligned BBH and BNS models to generic-orbit frameworks such as TEOBResumS-Dalí, which incorporates tidal interactions, generic spins, multipolar radiation reaction and waveform, and NR information for arbitrary orbits [2503.14580]. That work explicitly states that the model inherits both the post-adiabatic approximation for fast quasi-circular dynamics and **the EOBSPA for the computation of frequency-domain EOB waveforms** [2503.14580]. It also embeds NRPM or NRPMw as the postmerger completion for BNS systems, so that the model delivers a representation of the complete BNS spectrum within a unified EOB framework [2503.14580].

For BNS inspirals, the tidal sector of TEOBResum has itself been refined through resummed gravitoelectric and gravitomagnetic tidal potentials and tidal multipolar waveform corrections [1812.02744]. That improved nonspinning tidal model validated EOB energetics and phasing against NR and identified the resummed gravitoelectric octupolar term as the most important new conservative contribution, capable of producing up to 1 rad of dephasing relative to its nonresummed version depending on the neutron-star model [1812.02744]. This provides additional context for TEOBResumSPA\_NRPMw: the inspiral tidal parameters that determine \(\kappa_2^T\) and enter the EOB Hamiltonian are also the parameters used by NRPMw to predict postmerger frequencies and morphology [2205.09112], [1812.02744].

A common misconception is that NRPMw is merely an agnostic spectral fit. The construction described in the literature is more constrained: its wavelet parameters are tied to binary masses, tidal deformabilities, and aligned spins through EOS-insensitive relations, and the recalibration parameters are explicitly introduced to encode the residuals of those relations rather than to replace them [2205.09112]. Conversely, it would also be misleading to treat TEOBResumSPA\_NRPMw as fully universal outside its calibration region. The stated parameter ranges, EOS coverage, and spin limits are integral to the model definition [2205.09112].

Taken together, these elements define TEOBResumSPA\_NRPMw as a waveform architecture in which the inspiral–merger signal is provided by a calibrated EOB stationary-phase model and the kilohertz remnant signal is supplied by an analytical frequency-domain postmerger representation informed by hundreds of numerical-relativity simulations [2205.09112]. This suggests a synthesis between low-frequency tidal inference and kilohertz spectroscopy: the same binary parameters that shape the inspiral phase also determine the postmerger peak frequency \(f_2\), sideband structure, and collapse behavior within a single coherent model [2205.09112].

Source: https://www.emergentmind.com/topics/teobresumspa_nrpmw-waveform