---
title: Filtered-Squared Bispectrum
url: https://www.emergentmind.com/topics/filtered-squared-bispectrum
type: topic
---

# Filtered-Squared Bispectrum

The filtered-squared bispectrum (FSB) is a class of bispectrum estimators that compresses three-point information into a two-point statistic, typically the cross-power spectrum between a quadratic field built from a filtered map and a linear field. In three-dimensional large-scale structure this compression appears as density–quadratic-field cross-spectra, while in projected and full-sky analyses it is the cross-correlation between the square of a field filtered on a range of scales and the original field or a second field. In both settings, the estimator is a weighted integral of the underlying bispectrum, retains direct sensitivity to selected triangle families, and reuses much of the computational and covariance machinery developed for power-spectrum analysis [1411.6595][2409.07980].

## 1. Definition and compression principle

The defining operation is to filter a field, form a quadratic composite, and cross-correlate it with a linear field. On the flat sky, for a scalar field \(X(\theta)\) with filtered Fourier modes \(X_{f,\ell}=W(\ell)X_\ell\), one defines the squared map
\[
S(\theta)=[X_f(\theta)]^2,
\]
and the filtered-square cross-spectrum
\[
P^{SX}(L)\equiv V_L^{-1}\int_{|\ell|\in L} d^2\ell\,\langle S_\ell X_{-\ell}\rangle.
\]
With the bispectrum defined by \(\langle X_k X_q X_p\rangle=(2\pi)^2\delta_D(k+q+p)B(k,q,p)\), this becomes
\[
P^{SX}(L)=\int \frac{d^2\ell}{(2\pi)^2}W(\ell)W(|L-\ell|)B(\ell,|L-\ell|,L),
\]
so the FSB is an angle-averaged line integral through bispectrum configuration space [2409.07980].

The same compression structure appears in three-dimensional large-scale structure. If \(Q\) is a quadratic functional of the density field \(\delta\), the estimator is the cross-spectrum
\[
\hat P_{\delta,Q}(k)\equiv \langle \delta^R(k)Q_R(-k)\rangle_{\rm angles},
\]
where the superscript \(R\) denotes prior smoothing. Its expectation value is a bispectrum-weighted convolution over the kernel defining \(Q\), with one filtered linear leg and two filtered quadratic legs [1411.6595].

A central property of the FSB is that it turns a three-variable statistic into a one-variable data vector. This compression is not arbitrary: in the projected formulation, the estimator is most sensitive to triangles with two sides in the filter support and the third at the cross-spectrum multipole, so diagonal FSBs emphasize near-isosceles configurations, while generalized filtered-multiply constructions extend the coverage to neighboring or asymmetric bands [2409.07980]. In three-dimensional perturbation theory, the analogous compression can be chosen to isolate the Legendre components of the tree-level bispectrum kernel, making the estimator directly matched to the dominant gravitational and bias-induced bispectrum contributions [1411.6595].

## 2. Three-dimensional large-scale-structure formulation

For the matter or galaxy density contrast \(\delta(x)\), with Fourier transform \(\delta(k)\), the leading-order matter bispectrum in standard Eulerian perturbation theory is
\[
B_{mmm}(k_1,k_2,k_3)=2F_2(k_1,k_2)P_{\rm lin}(k_1)P_{\rm lin}(k_2)+2\,{\rm perms},
\]
with
\[
F_2(k_1,k_2)=\frac{5}{7}+\frac{1}{2}\left(\frac{k_1}{k_2}+\frac{k_2}{k_1}\right)\mu+\frac{2}{7}\mu^2,
\qquad \mu=\hat k_1\cdot \hat k_2.
\]
Using the Legendre basis \(P_0(\mu)=1\), \(P_1(\mu)=\mu\), and \(P_2(\mu)=\frac{3}{2}(\mu^2-1/3)\), this kernel decomposes into monopole, dipole, and quadrupole pieces with coefficients \(F_2^0=17/21\), \(F_2^1(k_1,k_2)=\frac{1}{2}(k_1/k_2+k_2/k_1)\), and \(F_2^2=4/21\) [1411.6595].

The FSB construction then introduces quadratic fields whose kernels match these angular structures. Writing
\[
Q(k)=\int \frac{d^3q}{(2\pi)^3}K_Q(q,k-q)\,\delta(q)\delta(k-q),
\]
the three basic quadratic fields are the squared density,
\[
Q_{\delta^2}(x)=\delta^2(x),\qquad K_{\delta^2}(q,p)=1,
\]
the shift term,
\[
Q_{\rm shift}(x)=-\Psi^i(x)\partial_i\delta(x),\qquad 
K_{\rm shift}(q,p)=-\frac{1}{2}\left(\frac{q}{p}+\frac{p}{q}\right)\mu,
\]
with displacement field \(\Psi_i(k)=-ik_i\delta(k)/k^2\), and the tidal term,
\[
Q_{s^2}(x)=s_{ij}(x)s_{ij}(x),\qquad 
K_{s^2}(q,p)=\frac{3}{2}(\mu^2-1/3),
\]
with
\[
s_{ij}(x)=\left(\partial_i\partial_j\nabla^{-2}-\delta_{ij}/3\right)\delta(x).
\]
These three operators form a complete basis for the angular dependence of the leading gravitational bispectrum and the leading nonlinear bias terms \(\delta^2\) and \(s^2\) [1411.6595].

To suppress nonlinear small-scale mode coupling, the fields are smoothed first. The standard choice is a Gaussian filter
\[
W_R(k)=\exp(-k^2R^2/2),
\qquad \delta^R(k)=W_R(k)\delta(k),
\]
with \(R\approx 20\,h^{-1}{\rm Mpc}\) in the baseline implementation. The quadratic fields are then built from \(\delta^R\), so the cross-spectrum \(P_{Q_R,\delta^R}(k)\) carries three factors of the smoothing window in ensemble averages [1411.6595].

This formulation has a direct physical interpretation. The shift cross-spectrum isolates the dipole component and, in the \(k\rightarrow 0\) limit and in the absence of velocity bias, has no independent nonlinear bias contribution. The squared-density and tidal cross-spectra isolate the monopole and quadrupole components associated with \(b_2\) and \(b_{s^2}\). In the large-scale limit the three cross-spectra become equal up to trivial window factors, but away from that limit they decorrelate and become complementary [1411.6595].

## 3. Near-optimality, bias estimation, and large-scale limits

The three-dimensional FSB was developed as a near-optimal estimator for bias determination. With Eulerian bias expanded to second order,
\[
\delta_h(x)=b_1\delta_m(x)+b_2[\delta_m^2(x)-\langle \delta_m^2\rangle]+\frac{2}{3}b_{s^2}[s_m^2(x)-\langle s_m^2\rangle],
\]
the matter–matter–halo cross-spectra satisfy
\[
P_{K_Q[\delta_m^R],\delta_h^R}(k)
=2b_1 I_{K_Q,F_2}^R+4b_1 I_{K_Q,F_2}^{\rm bare,R}
+2b_2 I_{K_Q,P_0}^R+\frac{4}{3}b_{s^2}I_{K_Q,P_2}^R.
\]
Subtracting \(b_1\) times the matter–matter result isolates the nonlinear bias terms,
\[
P_{K_Q[\delta_m^R],\delta_h^R}(k)-b_1P_{K_Q[\delta_m^R],\delta_m^R}(k)
=2b_2 I_{K_Q,P_0}^R+\frac{4}{3}b_{s^2}I_{K_Q,P_2}^R.
\]
For halo–halo–halo cross-spectra, the \(P_1\) shift piece is purely \(b_1\)-dependent, while the \(P_0\) and \(P_2\) pieces carry \(b_2\) and \(b_{s^2}\), so the three cross-spectra provide nearly orthogonal constraints on \(\{b_1,b_2,b_{s^2}\}\) [1411.6595].

The near-optimality claim follows from a maximum-likelihood derivation. For any separable theoretical bispectrum template \(B_{\rm theo}(k_1,k_2,k_3)=f(k_1)g(k_2)h(k_3)\), the weakly non-Gaussian maximum-likelihood estimator for its amplitude reduces to an inverse-variance weighted integral over a single cross-spectrum of a quadratic filtered density with a linearly filtered density,
\[
\hat f_{\rm NL}^{B_{\rm theo}}\propto \int dk\,\frac{k^2}{P(k)}\,
\hat P_{(f\delta/P)*(g\delta/P),\,h\delta}(k).
\]
Applied to the \(l=0,1,2\) pieces of the tree-level gravitational bispectrum, this yields the three quadratic fields \(\delta^2\), \(-\Psi\cdot\nabla\delta\), and \(s^2\) as matched filters. In practice, replacing exact inverse-variance weighting by Gaussian smoothing and using the theoretical covariance is sufficient to achieve near-optimality on large scales [1411.6595].

The same structure explains why joint fits to the power spectrum and FSB break the \(b_1\)–\(\sigma_8\) degeneracy. At fixed shape, the FSB amplitudes scale with the fluctuation amplitude through the linear power spectra inside the convolution integrals, while the three cross-spectra depend differently on \(b_1\), \(b_2\), and \(b_{s^2}\) than the power spectrum alone [1411.6595].

The large-scale limit makes this decomposition especially transparent. As \(k\rightarrow 0\), the halo–halo–halo cross-spectra approach a common expression involving \(b_1^3P_{\rm lin}(k)\), \(b_1^2b_2\), and \(b_1^2b_{s^2}\), but the shift cross-spectrum remains purely proportional to \(b_1\) at lowest order, whereas the \(\delta^2\) and \(s^2\) channels carry \(b_2\) and \(b_{s^2}\). A common misconception is therefore that the three estimators are redundant because they coincide on very large scales; in fact, their complementarity is realized on perturbative but nonzero \(k\), where they decorrelate [1411.6595].

## 4. Projected and full-sky generalizations

The projected FSB formalism recasts the same idea for maps on the plane or sphere. On the full sky, if \(a_{\ell m}\) are the spherical-harmonic coefficients of a scalar field and \((X_f)_{\ell m}=W_\ell a_{\ell m}\), then the full-sky FSB is the cross-power between the squared filtered map and the original map,
\[
\Phi_\ell^L\equiv \frac{1}{2\ell+1}\sum_m (X_f^2)_{\ell m}a_{\ell m}^*.
\]
Its expectation value is a linear functional of the reduced bispectrum,
\[
\Phi_b^L=\sum_{\ell_1\ell_2\ell_3} h_{\ell_1\ell_2\ell_3}^2
K_{\ell_1\ell_2\ell_3}^{LLb}\,b_{\ell_1\ell_2\ell_3},
\]
with
\[
h_{\ell_1\ell_2\ell_3}^2=\frac{(2\ell_1+1)(2\ell_2+1)(2\ell_3+1)}{4\pi}
\begin{pmatrix}
\ell_1&\ell_2&\ell_3\\
0&0&0
\end{pmatrix}^2,
\qquad
K_{\ell_1\ell_2\ell_3}^{LLb}=\frac{W_{\ell_1}W_{\ell_2}W_{\ell_3}^b}{N_b}.
\]
Thus the FSB is a bandpower-weighted compression of the full angular bispectrum [2409.07980].

Filter choice determines which triangles are emphasized. With a single filter used twice, the basic diagonal estimator weights triangles with two sides in the filter band and the third at the cross-spectrum multipole. Top-hat bands are simple and effective for wide large-scale-structure bands; smooth bands reduce ringing; matched filters can target specific shapes. A generalized filtered-multiply variant \(X_{f,{\rm low}}\times X_{f,{\rm high}}\) targets squeezed shapes more explicitly, and adding products of neighboring bands substantially increases triangle coverage [2409.07980].

Under masking, the estimator is naturally implemented in the pseudo-\(C_\ell\) framework. The practical pipeline is to filter the observed map in harmonic space, square it in real space, optionally re-mask the squared map, compute the pseudo-spectrum between the squared map and the observed field, and deconvolve the mask with the standard pseudo-\(C_\ell\) mixing matrix. Mean subtraction of the squared map removes the induced monopole, and a linear term built from Gaussian simulations can reduce variance in anisotropic conditions [2409.07980].

A mathematically sharper full-sky formulation appears in the optimal binned bispectrum literature. In the mask-free ideal limit, the binned estimator is
\[
\hat b_{\rm ideal}(\vec b)=F_{3,{\rm ideal}}^{-1}(\vec b)\,
\frac{1}{\Delta_3(\vec b)}\int d\Omega\,
H_{b_1}^{\rm ideal}[a]\,H_{b_2}^{\rm ideal}[a]\,H_{b_3}^{\rm ideal}[a],
\]
with filtered maps \(H_b^{\rm ideal}\) defined by harmonic binning and weighting. When two bins are equal, this becomes the filtered-squared statistic \(S=\int X\,Y^2\). In the presence of masking, beam convolution, or inpainting, the estimator remains unbiased after Fisher deconvolution, and the linear term is required for optimality at low \(\ell\) [2303.08828].

The projected formalism also extends directly to cross-bispectra of multiple fields. For maps \(\alpha\), \(\beta\), and \(\gamma\),
\[
\Phi_{LL\ell}^{\alpha\beta\gamma}\equiv \frac{1}{2\ell+1}\sum_m
(\alpha_L\beta_L)_{\ell m}\gamma_{\ell m}^*,
\]
compresses the bispectrum \(b_{\ell_1\ell_2\ell_3}^{\alpha\beta\gamma}\) with the same Wigner-\(3j\) kernel structure. This generalization underlies recent projected galaxy–galaxy–convergence analyses [2507.07968].

## 5. Covariance, implementation, and validation

A major attraction of the FSB is that its covariance is much simpler than the covariance of the full bispectrum. In three-dimensional large-scale structure, the leading Gaussian covariance of two cross-spectra at the same wavenumber is
\[
{\rm cov}\!\left[\hat P_{D[\delta_a^R],\delta_b^R}(k),
\hat P_{E[\delta_a^R],\delta_b^R}(k)\right]
=\frac{2}{N_{\rm modes}(k)}P_{bb}^R(k)\,
I_{DE}^{P_{aa}^RP_{aa}^R}(k),
\]
which implies the large-scale inverse-variance weight proportional to \(k^2/P_{bb}(k)\). For halo fields, Poisson stochasticity propagates into the cross-spectra, and phenomenological replacements \(1/\bar n\rightarrow 1/\bar n+\Delta_1\) and \(1/\bar n^2\rightarrow 1/\bar n^2+\Delta_2\) are used to model deviations from Poissonity [1411.6595].

For projected fields, the disconnected covariance can be treated as the covariance of a cross-power spectrum between a field and its squared filtered counterpart. In the full sky,
\[
{\rm Cov}_{\rm dc}[\Phi_\ell^L,\Phi_{\ell'}^{L'}]
=\frac{\delta_{\ell\ell'}}{2\ell+1}
\left[C_\ell^{aa}C_\ell^{a_L^2a_{L'}^2}
+C_\ell^{aa_L^2}C_\ell^{aa_{L'}^2}\right].
\]
With masks, the improved Narrow-Kernel Approximation (iNKA) reuses pseudo-\(C_\ell\) machinery with measured mode-coupled spectra. The dominant off-diagonal corrections are \(N_{222}\) for FSB auto-covariances and \(N_{32}\) for FSB–power cross-covariances, and these terms can be built from measured spectra and generalized FSBs in a model-independent manner [2409.07980].

The computational scaling follows the same logic. In the three-dimensional FFT implementation, one filters \(\delta(x)\), forms \(\delta_R^2(x)\), \(s_R^2(x)\), and \(-\Psi_i^R(x)\partial_i\delta^R(x)\), Fourier-transforms these quadratic fields, and measures three cross-spectra with \(\delta^R\). This reduces the task to three FFT-based power spectra rather than brute-force triangle counting with nominal cost \(O(N_{\rm grid}^6)\) [1411.6595]. On the sphere, each filter requires one forward spherical-harmonic transform, one inverse transform, one squaring step, and one pseudo-\(C_\ell\) cross-spectrum, so the total cost scales as \(O(B\times{\rm SHT})\) for \(B\) filters [2409.07980].

Validation results show that this compression is not merely formal. In the original three-dimensional study, with \(R\approx 20\,h^{-1}{\rm Mpc}\) at \(z\approx 0.55\), leading-order predictions agree with \(N\)-body results for matter–matter–matter and matter–matter–halo cross-spectra at the \(\lesssim 5\%\) level up to \(k\lesssim 0.09\,h\,{\rm Mpc}^{-1}\), while halo–halo–halo requires additional stochasticity corrections over the same range [1411.6595]. In projected analyses, 3D LPT full-sky simulations (100) and \(N\)-body GLAM lightcones (600) give FSB measurements unbiased within \(\pm 0.5\sigma\), and 2D LPT full-sky simulations (6000) show that the disconnected iNKA covariance underestimates the variance near filtered scales by \(\approx 20\%\), while adding \(N_{222}\) and \(N_{32}\) reproduces the full simulation covariance with \(\chi^2\) means differing from the simulation-based truth by \(<2.5\%\) [2409.07980].

The multi-field projected extension has also been validated on real survey geometry. In the DESI luminous red galaxy and Planck lensing analysis, the covariance combines Gaussian pseudo-\(C_\ell\) terms with multi-field generalizations of \(N_{222}\) and \(N_{32}\), rescaled by effective \(f_{\rm sky}\) factors calibrated with Gaussian simulations [2507.07968].

## 6. Relations, applications, and limitations

The FSB sits at the intersection of several pre-existing bispectrum compressions. In redshift-space-distortion studies, the skew-spectrum is the cross-power between a quadratic field and the density field; with the quadratic field taken as the square of a filtered density, it is precisely a filtered-squared bispectrum compression. The integrated bispectrum, by contrast, is a position-dependent power spectrum and correlates local small-scale power with the mean overdensity of a subvolume, so it targets the squeezed limit specifically. The line correlation function is different again: it uses phase information rather than amplitude-based quadratic filtering [1708.05303][1403.3411].

This relation is especially clear in the soft-limit literature. The position-dependent power spectrum implements a filtered–squared-field estimator of the squeezed-limit bispectrum, and in Fourier space the cross-spectrum of a squared small-scale filtered field with a long-wavelength field is
\[
P_{S,\delta}(q)=\int \frac{d^3k}{(2\pi)^3}F(k)F(|q-k|)\,B(q,k,-k-q),
\]
which is exactly the filtered-squared structure. In this sense, the integrated bispectrum, skew-spectrum, and FSB become equivalent weighted measures of squeezed configurations, differing mainly by the window and normalization conventions [1403.3411][1608.04345].

In CMB-style analyses, the FSB is closely analogous to KSW-like quadratic estimators: optimal or separable bispectrum estimators reduce to cross-correlating a linearly filtered field with a quadratic transform of the field. The full-sky optimal binned bispectrum estimator makes this correspondence exact, and the filtered-squared specialization \(S=\int X\,Y^2\) is simply the case with two equal bins [2303.08828].

Applications now span several observables. For the thermal Sunyaev–Zel'dovich effect, a filtered map \(T_f\) and its square \(S=T_f^2-\langle T_f^2\rangle\) yield either the scalar estimator \(E=\int d\hat n\,T_f S\) or the cross-spectrum \(C_\ell^{T_f\times S}\), whose expectation value is a Wigner-\(3j\)-weighted sum over the tSZ bispectrum. In that context the bispectrum amplitude scales as \(B_{\rm tSZ}\propto \sigma_8^{11\text{--}12}\), is principally sourced by massive clusters at redshifts around \(z\sim 0.4\), and is less sensitive to astrophysical uncertainties than the tSZ power spectrum at \(\ell\sim 3000\) [1203.6368].

For projected galaxy clustering and CMB lensing tomography, the multi-field FSB has enabled a combined power-spectrum and bispectrum analysis using DESI luminous red galaxies and Planck lensing. In that study the projected galaxy bispectrum is detected at very high significance, above \(30\sigma\) in all redshift bins, and the galaxy–galaxy–convergence bispectrum is detected above \(5\sigma\) in the three highest-redshift bins. The combination of \(\Phi^{ggg}\) with power spectra improves the \(\sigma_8\) precision by \(\sim 10\%\), while giving constraints consistent with the traditional \(2\times 2\)-point combination [2507.07968].

A more specialized extension appears in pulsar timing arrays. There the primary estimator is a filtered cubic statistic selecting a fixed frequency triangle \(f_1:f_2:f_3=a:b:c\), but in the squeezed limit \(c\ll a\approx b\) it approaches a filtered-squared-times-residual form, correlating a long-wavelength residual with a filtered square of short-wavelength residuals. The paper states that in the regime where the bispectrum is sharply peaked in the squeezed limit and the long mode is nearly constant across the observation, this filtered-squared estimator is equivalent up to normalization to the full bispectrum estimator with \((a,b,c)=(1,1,0)\) [1812.10615].

The main limitations are structural rather than conceptual. The basic diagonal FSB does not measure the full bispectrum shape space uniformly; it emphasizes near-isosceles configurations, and full shape recovery requires neighboring-band products or a filtered-multiply generalization. In three-dimensional perturbation theory, the choice of smoothing scale \(R\) trades signal against perturbative control: too small \(R\) increases signal but strains perturbation theory, while too large \(R\) suppresses signal-to-noise. The statement that the shift term has no independent nonlinear bias holds only in the \(k\rightarrow 0\) limit and in the absence of velocity bias. Real-data applications must additionally handle survey windows, redshift-space distortions, spin fields, or highly irregular sampling, and the current projected cosmology analyses therefore adopt conservative scale cuts and simple tree-level bias prescriptions [1411.6595][2409.07980][2507.07968].

Taken together, these formulations define the filtered-squared bispectrum as a general strategy for compressing bispectrum information into power-spectrum-like observables. Its distinguishing features are separable map-level construction, compatibility with FFT and pseudo-\(C_\ell\) pipelines, tractable covariance modeling, and the ability to target physically meaningful triangle families without reverting to full triangle-by-triangle estimation [1411.6595][2409.07980].

Source: https://www.emergentmind.com/topics/filtered-squared-bispectrum