---
title: Frequentist Regularised GW Mapping Analysis
url: https://www.emergentmind.com/topics/frequentist-regularised-gravitational-wave-mapping-analysis
type: topic
---

# Frequentist Regularised GW Mapping Analysis

Frequentist regularised gravitational-wave mapping analysis constitutes a foundational methodology for reconstructing the angular power distribution of gravitational-wave (GW) signals on the sky using data from detector networks such as pulsar timing arrays (PTAs) or ground-based interferometers. Based on a maximum-likelihood formalism, it utilizes explicit sky decompositions, detector response modeling, and sophisticated regularisation techniques to stabilize inversion in the presence of ill-posedness. This framework is crucial for producing GW sky maps, quantifying anisotropy, and resolving point sources, as well as for understanding the effective angular resolution achievable given instrumental and geometric limitations [1406.4511, 2601.13957, 0708.2728].

## 1. Mathematical Foundations: Sky Decomposition and Linear Data Model

The analysis begins with a decomposition of the GW field $h_A(\hat\Omega)$ (for polarizations $A\in\{+,\times\}$) on the sphere into an orthonormal basis $\{\phi_i(\hat\Omega)\}$, satisfying
\[
\int_{S^2} d\Omega\, \phi_i(\hat\Omega)\, \phi_j(\hat\Omega) = \delta_{ij}.
\]
Both pixel-space and spherical-harmonic expansions can be employed, with $h_A(\hat\Omega) = \sum_{i=1}^N s_{A,i}\, \phi_i(\hat\Omega)$. All coefficients are collected into a $2N$-vector $s = (s_{+,1},...,s_{+,N}; s_{\times,1},...,s_{\times,N})^T$ [1406.4511].

For PTAs, each pulsar $p$ yields timing-residual data $\delta t_p$ modeled as
\[
\delta t_p = \sum_{A,i} F_{p,A,i} s_{A,i} + n_p,
\]
where $F_{p,A,i}$ is the detector response matrix projected onto the chosen basis, incorporating geometric quantities (antenna patterns, pulsar terms, light travel delays). The vector form is $\delta t = F s + n$, with noise $n$ having Gaussian covariance $N$ [1406.4511, 0708.2728].

Alternatively, one can construct cross-correlated "dirty maps" from pairwise detector data (e.g., for LIGO/Virgo or PTA pairs), integrating over frequency and time with optimal filters, producing a linear model $\boldsymbol{\rho} = R \boldsymbol{P} +$ noise, where $R$ is the response matrix linking spherical-harmonic powers $P_{\ell m}$ to observations [2601.13957, 0708.2728].

## 2. Likelihood Formulation and Maximum-Likelihood Estimation

The frequentist analysis is rooted in a composite Gaussian likelihood:
\[
\ln \mathcal{L}(s) = -\frac{1}{2} [\delta t - F s]^T N^{-1} [\delta t - F s],
\]
which, analogously, holds for the cross-correlation vector in radiometer or PTA map-making formalisms. Maximization yields normal equations:
\[
F^T N^{-1} F\, \hat{s} = F^T N^{-1} \delta t,
\]
\[
R^T \Sigma^{-1} R\, \hat{\boldsymbol{P}} = R^T \Sigma^{-1} \boldsymbol{\rho},
\]
with Fisher or normal matrices $M$ ($M = F^T N^{-1} F$ or $M = R^T \Sigma^{-1} R$) and "dirty map" vectors $b$ or $\boldsymbol{X}$ [1406.4511, 2601.13957].

The formal unregularised maximum-likelihood (ML) map solution is:
\[
\hat{s} = M^{-1} b, \quad \hat{\boldsymbol{P}} = M^{-1} \boldsymbol{X}.
\]
However, in practice, $M$ is typically ill-conditioned due to limited detector number, non-uniform sky sensitivity, and incomplete coverage, necessitating further regularisation [1406.4511, 2601.13957, 0708.2728].

## 3. Regularisation: Tikhonov and Truncated SVD

Direct inversion of $M$ amplifies noise in weakly constrained sky modes. Two primary regularisation strategies are employed:

- **Tikhonov (Ridge) Regularisation:** Addition of a penalty term $\lambda R$ (typically $R=I$) yields $M_\lambda = M + \lambda I$, with the regularised solution $\hat{s}(\lambda) = M_\lambda^{-1} b$. The parameter $\lambda$ is chosen via L-curve, cross-validation on noise, or empirical Bayes [1406.4511, 0708.2728].

- **Truncated Singular Value Decomposition (SVD):** $M$ is decomposed as $M = V \Xi V^T$ with eigenvalues $\{\xi_k\}$. Only the $r \leq N_\text{psr}$ largest singular values are retained, forming a pseudo-inverse $\tilde M^{-1} = V \Xi^+ V^T$ with $\Xi^+_{ii} = 1/\xi_i$ for $i \le r$, else $0$ [2601.13957]. This procedure discards noise-dominated sky modes, providing a stable "clean map" $\tilde{\boldsymbol{P}} = \tilde M^{-1} \boldsymbol{X}$ and corresponding uncertainties [1406.4511, 2601.13957, 0708.2728].

Table: Key Regularisation Approaches

| Regularisation Type      | Modification                | Typical Use Case                      |
|-------------------------|-----------------------------|---------------------------------------|
| Tikhonov                | $M_\lambda = M + \lambda I$ | General ill-posed linear inversion    |
| Truncated SVD           | Keep $r$ largest $\xi_k$    | Sky mapping with limited baselines    |

These regularisation approaches are robust, tractable, and adaptable to detector geometry and noise properties.

## 4. Sky-Map Statistics: Isotropy, Anisotropy, and Point-Source Recovery

### Isotropic Limit and Hellings–Downs Cross-Correlation

For a statistically isotropic GW background, the covariance structure simplifies:
\[
\langle h_A(\hat\Omega) h_B(\hat\Omega') \rangle = \frac{S_h}{2} \delta(\hat\Omega, \hat\Omega') \delta_{AB},
\]
leading to the classic Hellings–Downs overlap reduction function (for PTAs, $\Gamma_{pp'}$), and just the monopole sky mode survives regularisation. Marginalisation over sky coefficients analytically recovers the standard cross-correlation statistic [1406.4511].

### Anisotropic Backgrounds and Point Sources

No assumption of isotropy is necessary. Anisotropy is modeled by retaining higher spherical harmonic modes or pixel basis elements. For point-source recovery, localized delta-function templates $\phi^{(\mathrm{ps})}_k(\hat\Omega) = \delta(\hat\Omega, \hat\Omega_k)$ are included in the basis, with their amplitudes estimated via the same regularised normal equation machinery [1406.4511]. Broader prior covariances can be assigned for more flexibility in capturing extended anisotropy [1406.4511, 2601.13957].

## 5. Quantifying and Optimizing Angular Resolution

A defining innovation is the characterization and optimization of angular resolution using the point spread function (PSF) and distortion matrix formalism [2601.13957]. The PSF directly quantifies the sky-patch area to which a true point source's recovered flux is spread due to the measurement process and regularisation. The distortion matrix is constructed as:
\[
\Lambda^{\rm pix}_r = U (\tilde M^{-1} M) U^T,
\]
where $U$ relates the pixel and spherical-harmonic bases. The effective PSF area at each sky location $i$ is given by summing pixels above half-max:
\[
A_{\rm PSF}(i) = \sum_{\hat\Omega} \left[\Lambda^{\rm pix}_r \right]_{\hat\Omega,i} > \tfrac12 \max_{\hat\Omega} \left[\Lambda^{\rm pix}_r \right]_{\hat\Omega,i} \times \Delta\Omega_{\rm pix},
\]
yielding a PSF map $A_{\rm PSF}(\hat\Omega)$ across the sky. In regions with dense detector coverage, $A_{\rm PSF}$ approaches the geometric limit set by the nearest-neighbor spacing, expressed as
\[
A_{\rm nn}=2\pi [1 - \cos(\delta_{\rm nn}/2)].
\]
Simulations show that local resolution closely tracks the underlying detector geometry, and optimal $\ell_\text{max}$ is chosen such that the scale $\pi/(\ell_\text{max}+1) \simeq \bar\delta_{\rm nn}$, with $r \leq N_\text{psr}$ [2601.13957].

A significant improvement over earlier methods is the implementation of an adaptive, variable local-resolution scheme, allowing finer resolution where the detector network is more densely populated, as confirmed by MeerKAT PTA simulations [2601.13957].

## 6. Numerical Implementation and Computational Aspects

The map-making pipeline follows a sequence:

1. Choosing a sky basis (pixels, spherical harmonics).
2. Computing frequency-dependent detector response matrices.
3. Measuring data cross-spectra, assembling noise or covariance matrices.
4. Whitening data and response matrices ($Q$ such that $QQ^T = N^{-1}$).
5. SVD or eigendecomposition of the compactified response.
6. Regularisation via Tikhonov or truncation.
7. Solving for sky coefficients in the low-dimensional subspace.
8. Inverting transforms to reconstruct sky maps ($\hat h_A(\hat\Omega)$).
9. Quantifying angular resolution with PSF maps.
10. Examining higher-order or localized coefficients for anisotropy or individual source detection.

Table: Computational Steps in Frequentist Regularised GW Mapping

| Step | Major Operation                                     | Reference                |
|------|-----------------------------------------------------|--------------------------|
| 1    | Basis choice, response modeling                     | [1406.4511, 2601.13957]  |
| 2–3  | Data/covariance assembly                            | [1406.4511, 2601.13957]  |
| 4    | Whitening, low-rank projection                      | [1406.4511]              |
| 5–6  | SVD/truncation, regularisation                      | [2601.13957]             |
| 7–8  | Sky map solution and inverse transform              | [1406.4511]              |
| 9    | Local PSF/area calculation                         | [2601.13957]             |

The dominant computational costs arise in beam and response calculation, typically $O(N_\text{pix}^2 N_t N_f)$, which are handled by exploiting matrix symmetry, parallelisation, and fast transform techniques (e.g., FFT interpolation) [0708.2728]. The rank of the relevant matrices never exceeds $2N_\text{psr}$, allowing efficient low-rank representation and inversion [1406.4511].

## 7. Extensions, Applications, and Performance in Practice

The frequentist regularised framework is applicable across detector types, including GW radiometers for ground-based interferometers [0708.2728] and PTAs [1406.4511, 2601.13957]. It extends naturally to:

- Multiple-baseline, multi-detector networks by stacking response and noise matrices [0708.2728].
- Polarised backgrounds by including amplitude coefficients for each polarisation [0708.2728].
- Empirical-Bayes embedding for hyperparameter selection (e.g., regularisation strength) [1406.4511].
- Adaptive-resolution mapping to match local detector density, yielding significant gains (e.g., the clean-map S/N for a simulated continuous GW increased by ~factor 2 and hotspot area shrank by ~40% in MeerKAT simulations) [2601.13957].

In the isotropic limit, the method reduces to the standard Hellings–Downs cross-correlation statistic. For anisotropic or point-like signals, the framework provides a rigorous basis for statistical significance quantification and limits spurious anisotropy artifacts. Interpretation of map coefficients, PSF structure, and detection statistics are geometrically transparent within this formalism.

This approach provides a unified and robust methodology for extracting GW sky maps, probing anisotropy, and quantifying the true resolving power of current and near-future GW detector networks [1406.4511, 2601.13957, 0708.2728].

Source: https://www.emergentmind.com/topics/frequentist-regularised-gravitational-wave-mapping-analysis