5-Vector Resampling in GW Searches

• 5-vector resampling is a technique that demodulates binary orbital effects by resampling the time series, recovering the coherent signal-to-noise ratio lost due to Doppler and relativistic modulations.
• It utilizes a five-vector formalism to represent the residual sidereal amplitude modulation as five Fourier components, thereby concentrating the signal power for improved detection.
• The method underpins directed searches in advanced gravitational-wave detectors like LIGO, balancing sensitivity recovery and computational cost in targeted studies such as that of Scorpius X-1.

Searching arXiv for the cited five-vector resampling papers and closely related work. {"query":"five-vector resampling continuous wave Scorpius X-1 arXiv", "max_results": 10} {"query":"A directed continuous-wave search from Scorpius X-1 with the five-vector resampling technique", "max_results": 5} {"query":"(Amicucci et al., 10 Mar 2025)", "max_results": 5} 5-vector resampling is a data-analysis technique for directed and narrowband searches for continuous gravitational waves from neutron stars in binary systems. It combines time-domain resampling, which removes Doppler and relativistic phase modulations, with the five-vector formalism, which represents the residual sidereal amplitude modulation in five Fourier components. In the formulation applied to Scorpius X-1 and to directed searches from neutron stars in binary systems, the method is designed to recover the coherent signal-to-noise ratio otherwise lost when orbital motion spreads signal power across many frequency bins, while maintaining an affordable computational cost for searches over uncertain orbital parameters (Amicucci et al., 10 Mar 2025, Amicucci et al., 25 May 2025).

1. Physical setting and search motivation

Continuous gravitational-wave signals are treated as nearly monochromatic in the neutron star rest frame, with phase

$$\Phi^{\rm NS}(\tau)=\phi_0+2\pi f_{\rm GW}(\tau-\tau_0).$$

For sources in binaries, the phase observed at an Earth-based detector is modulated both by the Earth’s orbital and rotational motion and by the neutron star’s orbital motion around the binary barycentre. The latter is typically described by five Keplerian parameters, and inaccurate treatment of those parameters causes the signal to be smeared over multiple FFT bins, reducing the signal-to-noise ratio and potentially hindering detection (Amicucci et al., 25 May 2025).

In the detector frame, the phase can be written as

$$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$

where $\Delta\tau(t)$ includes Earth-barycentric corrections and the binary delay. In the low-eccentricity description emphasized in the method summary, the binary Rømer delay depends on the five Keplerian parameters $\{P,a_{\rm p},e,\omega,t_{\rm p}\}$; equivalently, for $e\ll 1$, one may use Laplace–Lagrange parameters $\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$ with $\Omega=2\pi/P$, $\kappa=e\cos\omega$, $\eta=e\sin\omega$, and $t_{\rm asc}\equiv t_{\rm p}+\omega/\Omega$ (Amicucci et al., 25 May 2025).

The search motivation is astrophysical as well as algorithmic. The Advanced LIGO-Virgo-KAGRA detectors operate in bands containing more than half of the known pulsars in the Galaxy existing in binary systems, and the method was developed specifically to conduct thorough directed and narrowband continuous-wave searches over such targets at affordable computational cost (Amicucci et al., 10 Mar 2025).

2. Phase demodulation through resampled time

The central step of 5-vector resampling is the construction of a new time coordinate that unwinds the detector-frame phase modulation. In the Scorpius X-1 implementation, the detector phase is written

$$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$0

with $$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$1 the projected radial binary delay and $$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$2 the geometric barycentric term (Amicucci et al., 10 Mar 2025).

The resampled time variable is then defined as

$$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$3

or, equivalently in the compact notation of the broader method description,

$$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$4

By construction, the phase in the new coordinate is strictly linear,

$$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$5

so the astrophysical signal becomes nearly strictly monochromatic in the resampled time series (Amicucci et al., 10 Mar 2025, Amicucci et al., 25 May 2025).

Operationally, the original uniformly sampled detector data $$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$6 are resampled onto the irregular grid $$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$7 and represented as a new time series $$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$8. The method summary explicitly notes that no interpolation of frequency is needed because the resampling is performed in the time domain (Amicucci et al., 10 Mar 2025). This is the defining distinction of the technique: the orbital and barycentric corrections are absorbed into the time coordinate rather than re-applied separately to every trial frequency.

3. Five-vector formalism and matched filtering

After phase demodulation, the signal retains the Earth’s sidereal amplitude modulation through the detector antenna patterns. In the Scorpius X-1 formulation,

$$\Phi^{\rm det}(t_{\rm arr})=\phi_0+2\pi f_{\rm GW}\bigl[t_{\rm arr}+\Delta\tau(t_{\rm arr})-\tau_0\bigr],$$9

where $\Delta\tau(t)$0 are known sidereal antenna patterns (Amicucci et al., 10 Mar 2025).

The five-vector is constructed by taking Fourier components at the five sidereal sidebands. For a real time series $\Delta\tau(t)$1,

$\Delta\tau(t)$2

and

$\Delta\tau(t)$3

The Earth sidereal frequency is given as $\Delta\tau(t)$4, so the signal is represented by five components at $\Delta\tau(t)$5 with $\Delta\tau(t)$6 (Amicucci et al., 10 Mar 2025, Amicucci et al., 25 May 2025).

Under the signal model, the data five-vector is

$\Delta\tau(t)$7

with $\Delta\tau(t)$8 the five-vectors of the sidereal templates and $\Delta\tau(t)$9 the noise vector (Amicucci et al., 10 Mar 2025). The broader formulation writes this as a $\{P,a_{\rm p},e,\omega,t_{\rm p}\}$0 linear model,

$\{P,a_{\rm p},e,\omega,t_{\rm p}\}$1

from which a generalized-likelihood-ratio statistic can be formed (Amicucci et al., 25 May 2025).

In the Scorpius X-1 pipeline, polarization estimators are computed as

$\{P,a_{\rm p},e,\omega,t_{\rm p}\}$2

and the detection statistic is

$\{P,a_{\rm p},e,\omega,t_{\rm p}\}$3

The companion methodological summary presents the corresponding optimal statistic as $\{P,a_{\rm p},e,\omega,t_{\rm p}\}$4, where $\{P,a_{\rm p},e,\omega,t_{\rm p}\}$5 is the inverse noise-covariance matrix in five-vector space (Amicucci et al., 25 May 2025). Taken together, these descriptions show that the method pairs explicit binary demodulation with a coherent sidereal matched filter.

4. Pipeline realization in directed and narrowband searches

The high-level implementation described for the Scorpius X-1 search begins with construction of the SFDB (Short-FFT Database). The full band is divided into 1-Hz sub-bands, short FFTs are formed, and calibrated strain is recorded. For each point in the search grid over orbital-parameter offsets $\{P,a_{\rm p},e,\omega,t_{\rm p}\}$6, the mapping $\{P,a_{\rm p},e,\omega,t_{\rm p}\}$7 is computed and the original time series is interpolated to obtain $\{P,a_{\rm p},e,\omega,t_{\rm p}\}$8 on a uniform grid. This time-domain interpolation is the most expensive step in CPU terms (Amicucci et al., 10 Mar 2025).

The broader methodological description expresses the same logic in a band-limited form. Starting from calibrated evenly sampled data $\{P,a_{\rm p},e,\omega,t_{\rm p}\}$9, one computes $e\ll 1$0, forms the resampled timestamps $e\ll 1$1, chooses a uniform grid in $e\ll 1$2 with $e\ll 1$3 to satisfy Nyquist sampling to $e\ll 1$4, and interpolates with a high-order kernel. The resampled series may then be heterodyned by $e\ll 1$5, low-pass filtered over a narrow bandwidth $e\ll 1$6, and downsampled before an FFT is taken and the five sideband bins are extracted (Amicucci et al., 25 May 2025).

Candidate selection is not based solely on the raw detection statistic. In the directed Scorpius X-1 pipeline, the noise-only distribution $e\ll 1$7 is estimated, a threshold $e\ll 1$8 is set by fixing the false-alarm rate through the look–elsewhere factor, and the statistic is normalized to

$e\ll 1$9

to provide robustness against band-by-band variations and instrumental lines. Candidates with $\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$0 are retained. They are then subjected to an internal same-detector coincidence veto requiring at least $\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$1 out of the 9 possible sidereal-sideband peaks to exceed threshold, followed by an inter-detector coincidence condition requiring appearance in both Hanford and Livingston within $\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$2 of the expected sidereal sidebands (Amicucci et al., 10 Mar 2025).

5. Sensitivity recovery and computational scaling

The method is motivated by a specific sensitivity loss mechanism. Without demodulation, orbital motion spreads the signal over $\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$3 bins, where the orbital-Doppler bandwidth is given as $\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$4. The coherent SNR is therefore reduced by approximately $\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$5. After exact resampling, the signal power is refocused into the five sidereal bins, and the SNR gain is summarized as

$\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$6

The Scorpius X-1 summary states the same point qualitatively: longer coherence time raises SNR as $\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$7, and the five-vector filter concentrates all signal power into nine peaks, yielding a matched-filter boost (Amicucci et al., 25 May 2025, Amicucci et al., 10 Mar 2025).

The computational scaling is likewise explicit. For one 1-Hz band and one Keplerian template in one detector, the total CPU cost is approximated by

$\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$8

with $\{\Omega,a_{\rm p},\kappa,\eta,t_{\rm asc}\}$9 for band extraction, $\Omega=2\pi/P$0 for resampling, $\Omega=2\pi/P$1 for the five-vector filter, and $\Omega=2\pi/P$2 for coincidence checking. For $\Omega=2\pi/P$3 1-Hz bands, $\Omega=2\pi/P$4 orbital templates, and $\Omega=2\pi/P$5 detectors,

$\Omega=2\pi/P$6

The dominant cost is therefore the time-domain resampling step (Amicucci et al., 10 Mar 2025).

A crucial optimization is that the time-resampling step is independent of frequency. A single resampling covers the entire 1 Hz band, and the same resampled time series can be reused for all sub-bins. The method summary identifies this as a major saving relative to recalculating Doppler corrections for each frequency and notes that this feature makes the approach ideally suited to narrowband searches of unknown spin frequency (Amicucci et al., 10 Mar 2025). A plausible implication is that the method’s practical niche is precisely the regime in which orbital-parameter uncertainty is significant but the source sky position is already well constrained.

6. Application to Scorpius X-1 and methodological limits

The first application described in the supplied material is a directed search for continuous waves from Scorpius X-1 using publicly available data from the third observing run of Advanced LIGO-Virgo-KAGRA. The search used the full O3 run from Hanford and Livingston over a frequency range of 10–1000 Hz and tested three ephemeris configurations: the nominal Scorpius X-1 parameters and two edge-of-uncertainty offsets in $\Omega=2\pi/P$7 (Amicucci et al., 10 Mar 2025, Amicucci et al., 25 May 2025).

No statistically significant continuous-wave signal was claimed, and no candidates survived the follow-up chain. The search therefore set 95% confidence-level upper limits in selected frequency bands and orbital-parameter ranges while also evaluating overall sensitivity (Amicucci et al., 10 Mar 2025). The method summary reports that continuous-wave signals down to strain amplitudes $\Omega=2\pi/P$8 can be recovered at 95% confidence in O3 data for Scorpius X-1, while the broader summary gives a best 95% confidence-level upper limit of $\Omega=2\pi/P$9 at $\kappa=e\cos\omega$0 and, more specifically, a best value $\kappa=e\cos\omega$1 at $\kappa=e\cos\omega$2 in LIGO Livingston (Amicucci et al., 10 Mar 2025, Amicucci et al., 25 May 2025).

Several practical limitations are emphasized. The robustness of the method depends on covering uncertain orbital parameters within maximum offsets that still produce a coherent peak. For Scorpius X-1, the electromagnetic uncertainties exceed those offsets, so three searches—null, “+”, and “−” offsets—were used to cover the ephemeris. Known spectral lines and their Doppler-broadened replicas must be zeroed before fitting the noise tail; otherwise the threshold can be overestimated by up to an order of magnitude. The method is therefore not a generic black-box demodulator: its performance depends on orbital ephemeris control, line handling, and a coincidence framework tailored to the sidereal-sideband structure (Amicucci et al., 10 Mar 2025).

A recurring misconception is that resampling renders all detector modulation irrelevant. The formalism shows otherwise. Resampling removes the phase modulation associated with barycentric and binary delays, but the Earth’s sidereal amplitude modulation remains and is intentionally retained as the structure exploited by the five-vector matched filter (Amicucci et al., 10 Mar 2025). That separation between phase demodulation and sidereal filtering is the defining feature of 5-vector resampling.