---
title: 'MaxWave: Rapid Gravitational-Wave Reconstruction'
url: https://www.emergentmind.com/topics/maxwave
type: topic
---

# MaxWave: Rapid Gravitational-Wave Reconstruction

Searching arXiv for MaxWave and closely related gravitational-wave reconstruction papers.
MaxWave is a gravitational-wave data-analysis framework for wavelet-based maximum-likelihood reconstruction of non-Gaussian structure in detector data. Introduced in "MaxWave: Rapid Maximum Likelihood Wavelet Reconstruction of Non-Gaussian features in Gravitational Wave Data" [2508.13377], it is described as an improvement of the existing BayesWave initialization algorithm and as a rapid, low-latency approximate maximum likelihood method for reconstructing non-Gaussian features. Its stated scope includes de-noising long-duration signals, reconstructing glitches in forms useful for machine-learning analysis, and, in a later extension, coherent multi-detector signal reconstruction and glitch rejection in real time [2606.23035].

## 1. Provenance, naming, and scope

MaxWave was published on 2025-08-18 as a gravitational-wave method centered on rapid maximum-likelihood wavelet reconstruction of non-Gaussian features [2508.13377]. A subsequent paper, "MaxWave Signal: Rapid, coherent maximum likelihood wavelet reconstruction of transient signals in gravitational wave data," published on 2026-06-22, extends the framework to coherent multi-detector signal reconstruction and glitch rejection [2606.23035].

The term should be distinguished from "WaveMax: Radar Waveform Design via Convex Maximization of FrFT Phase Retrieval," which is a radar waveform-design method based on ambiguity-function phase retrieval and convex optimization rather than a gravitational-wave reconstruction algorithm [2501.14164]. The similarity of names can obscure the fact that the two works address different inverse problems, different sensing modalities, and different mathematical objects.

Within gravitational-wave analysis, MaxWave is explicitly positioned relative to BayesWave: the 2025 paper states that it improves the existing BayesWave initialization algorithm [2508.13377]. The later coherent extension further states that it can complement existing burst search and reconstruction frameworks through a fundamentally distinct approach [2606.23035].

## 2. Problem setting in gravitational-wave data analysis

The motivating context is the increase in transient-event rates driven by improved detector sensitivity. The 2025 MaxWave paper states that advancements in the sensitivity of gravitational-wave detectors have increased the detection rate of transient astrophysical signals, while the 2026 extension states that faster automated analysis is therefore required [2508.13377].

MaxWave targets reconstruction of non-Gaussian features in detector data. In the gravitational-wave context, this category includes at least glitches and transient signal structure that departs from stationary Gaussian noise assumptions. The 2025 abstract emphasizes that the algorithm produces a rapid, low-latency approximate maximum likelihood solution for such features and can be used to de-noise long-duration signals, including the stochastic gravitational-wave background from numerous unresolved sources and continuous wave signals from isolated sources such as rotating neutron stars [2508.13377].

A further stated use-case is machine-learning support. MaxWave is described as providing non-whitened, time and frequency domain, Gaussian-noise-free reconstructions of any glitch, with the intended benefit of aiding applications that classify and understand glitches [2508.13377]. This suggests a dual role: direct reconstruction for data conditioning and feature extraction, and indirect support for downstream discriminative or representation-learning pipelines.

## 3. Core methodological modifications in the original MaxWave formulation

The original MaxWave abstract identifies three specific enhancements over the existing BayesWave initialization algorithm [2508.13377].

| Enhancement | Stated purpose |
|---|---|
| Modified wavelet basis | Eliminate redundant inner product calculations |
| Shift from TFQ to TF$\tau$ transforms | Optimize wavelet subtractions |
| Downsampled heterodyned wavelet transform | Accelerate initial calculations |

The first modification uses a modified wavelet basis to eliminate redundant inner product calculations [2508.13377]. In practical terms, this locates the efficiency gain in the linear-algebraic core of wavelet matching rather than only in search heuristics.

The second modification replaces traditional time-frequency-quality factor wavelet transforms with time-frequency-time extent transforms, written as TF$\tau$, in order to optimize wavelet subtractions [2508.13377]. Because the abstract pairs TFQ and TF$\tau$ directly, the intended emphasis is not merely a change of parameterization but a change aimed at more efficient subtraction behavior during reconstruction.

The third modification is a downsampled heterodyned wavelet transform used to accelerate initial calculations [2508.13377]. Together with the other two changes, this places MaxWave in the class of reconstruction methods that seek latency reduction by restructuring the transform and basis evaluation stages rather than by abandoning likelihood-based modeling.

## 4. Maximum-likelihood wavelet modeling in the MaxWave lineage

The 2025 abstract identifies MaxWave as an approximate maximum-likelihood reconstruction method but does not, in the material summarized here, provide the explicit objective function. A later extension, MaxWave Signal, makes the statistical structure explicit for the coherent multi-detector case by modeling each detector time series as
$$
D_i(t)=h_i(t)+n_i(t),
$$
with zero-mean, stationary Gaussian noise characterized by one-sided power spectral density $S_i(f)$, and likelihood
$$
\mathcal{L}(h\,|\,D)\propto \exp\!\left[-\frac12\sum_{i=1}^K \langle D_i-h_i \mid D_i-h_i\rangle_i\right].
$$
The corresponding log-likelihood is written as
$$
\ell(h)= -\frac12\sum_i \int_0^\infty \frac{|D_i(f)-h_i(f)|^2}{S_i(f)}\,df+\mathrm{const.}
$$
[2606.23035]

In that extension, the waveform is represented as a sparse superposition of complex wavelets,
$$
h(t)=\sum_{k=1}^{N_w} A_k\,\psi_k(t;\theta_k),
$$
with parameters
$$
\theta_k=\{t_0^k,f_0^k,Q^k,\phi_0^k\},
$$
and a frequency-domain sine-Gaussian form
$$
\Psi_k(f)=A_k \exp\!\left[-\frac{Q_k^2}{4}\left(\frac{f-f_0^k}{f_0^k}\right)^2\right] e^{-2\pi i f t_0^k+i\phi_0^k},
$$
together with $\tau_k=Q_k/(2\pi f_0^k)$ [2606.23035]. At fixed wavelet morphology, maximizing the likelihood over amplitudes yields the linear system
$$
\sum_{j=1}^{N_w} M_{kj}A_j=b_k,\qquad
M_{kj}=\sum_i \langle \psi_k\mid \psi_j\rangle_i,\qquad
b_k=\sum_i \langle \psi_k\mid D_i\rangle_i,
$$
so that
$$
\mathbf{A}=\mathbf{M}^{-1}\mathbf{b}.
$$

These equations come from the coherent extension rather than the 2025 abstract itself, but they clarify what the designation "maximum likelihood wavelet reconstruction" means within the MaxWave program: sparse wavelet morphology, noise-weighted inner products, and analytic or semi-analytic amplitude estimation embedded inside a nonlinear search over wavelet parameters [2606.23035].

## 5. Coherent multi-detector extension and glitch rejection

MaxWave Signal extends the framework from rapid single-stream reconstruction to coherent network analysis. Its procedure begins by running the single-detector MaxWave glitch model on each detector and selecting the detector with largest SNR as the dominant reconstruction; the other detectors are aligned to it using the FINDCHIRP $z$-statistic
$$
z_n(\Delta t)=4\int_0^\infty \frac{\tilde d_n(f)\,\tilde h^*(f)}{S_n(f)} e^{2\pi i f \Delta t}\,df,
$$
whose magnitude peak determines the time delay $\Delta t_n$, whose argument gives the phase shift $\Delta \phi_n$, and whose normalized magnitude gives the amplitude scale $a_n$ [2606.23035].

The aligned detector streams are then combined through adaptive noise weighting based on the geometric mean PSD
$$
S_{\rm geo}(f)=\left[\prod_{i=1}^K S_i(f)\right]^{1/K},
$$
with effective strain
$$
d_{\rm eff}(f)=\sum_{i=1}^K \frac{S_{\rm geo}(f)}{a_i S_i(f)}\,d'_i(f),
$$
so that the network inner product reduces to a single-detector form with effective PSD $S_{\rm eff}(f)=S_{\rm geo}(f)$ [2606.23035]. The paper states that this synthetic detector amplifies non-Gaussian features, down-weights noisy detectors, and preserves Gaussian noise statistics.

Glitch rejection is implemented in two stages. First, only detectors satisfying the inter-detector light-travel-time condition of approximately $20\,\mathrm{ms}$ and SNR $\ge 5$ enter the coherent combination. Second, after reconstructing the coherent signal on the synthetic detector, the method shifts and rescales the coherent reconstruction back to each detector, forms residuals, and re-runs MaxWave; if the residual reconstruction has more than one wavelet, the event is flagged as a coincident-event non-removal, i.e. glitch [2606.23035].

The extension reports quantitative performance. For GW150914, the coherent detector Q-scan shows clear amplification of signal power, and the time-domain match with the NR template improves from $93\%$ in H1 alone and $81\%$ in L1 alone to $94\%$ after coherence. In injections in real LIGO O3 noise, the coherent match improves by $6.3_{-1.7}^{+3.3}\%$ for network SNR $=10$ and by $3.2_{-0.7}^{+1.2}\%$ for SNR $=20$. The $z$-statistic test alone rejects $74.4\%$ of $2{,}960$ artificial simultaneous glitches; adding the residuals test raises rejection to $94.0\%$, while cross-type glitches reject at $95.9\%$. In a real $14\,\mathrm{h}$ O3 run with H1+L1, the method recovered GW190521_030229 at SNR $=14.8$ and GW190521_074359 at SNR $=25.1$ at the correct GPS times, with no false alarms in a $1\,\mathrm{s}$ time-shift analysis, corresponding to $\mathrm{FAR}\le 1.99\times 10^{-5}\,\mathrm{Hz}$, reported as less than $2$/day. Runtime scaling on a single CPU core for $T_{\rm obs}=4\,\mathrm{s}$ is given as $0.8$--$1.32\,\mathrm{s}$ for $2$ detectors and $6$ wavelet layers, $1.16$--$1.94\,\mathrm{s}$ for $3$ detectors, and $1.88$--$3.18\,\mathrm{s}$ for $5$ detectors; even with $8$ layers, runtime remains below $1.4\,\mathrm{s}$ for $2$ detectors and below $3.4\,\mathrm{s}$ for $5$ detectors [2606.23035].

## 6. Applications, significance, and relation to adjacent work

The original MaxWave abstract identifies three application domains. First, it can de-noise long-duration signals, specifically including the stochastic gravitational-wave background from numerous unresolved sources and continuous wave signals from isolated sources such as rotating neutron stars. Second, it reconstructs non-Gaussian features in a form suitable for low-latency analysis. Third, it can aid machine-learning applications by providing non-whitened, time and frequency domain, Gaussian-noise-free reconstructions of any glitch [2508.13377].

The coherent extension adds a fourth application domain: real-time, model-independent signal reconstruction that strengthens detection confidence and improves sensitivity to diverse signal morphologies while safely denoising gravitational-wave data without removing transient signals [2606.23035]. In that sense, MaxWave occupies a methodological position between fully Bayesian, potentially higher-latency reconstruction pipelines and purely discriminative trigger-generation systems.

A common misconception is to conflate MaxWave with WaveMax because of the near-anagrammatic titles. The distinction is categorical. MaxWave concerns wavelet reconstruction of non-Gaussian features in gravitational-wave detector data and coherent multi-detector signal extraction [2508.13377; 2606.23035]. WaveMax, by contrast, addresses radar waveform recovery from ambiguity-function magnitudes using the fractional Fourier transform and a convex low-rank formulation [2501.14164].

The present public summary of the original 2025 MaxWave paper is abstract-level: it specifies the algorithm’s stated enhancements, target use-cases, and low-latency objective, while the later 2026 extension provides the most explicit equations, workflow, and performance figures available in the material summarized here [2508.13377; 2606.23035]. This suggests that the MaxWave research program is best understood as a fast maximum-likelihood wavelet-reconstruction family whose defining themes are computationally efficient basis evaluation, low-latency subtraction and denoising, and increasingly coherent exploitation of multi-detector information.

Source: https://www.emergentmind.com/topics/maxwave