REVIEW
Non-Gaussian timing residual skewness in pulsar timing arrays
General Relativity and Quantum Cosmology · Submitted 2 Oct 2026 · 124 references
Pith's verdict: The polarization-frame covariance critique is clean and structural; the nonzero skewness benchmark is an explicitly acknowledged channel-level construction whose connection to a physical SGWB remains open.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
The paper shows that the simple CMB-inspired quadratic non-Gaussian ansatz for gravitational-wave amplitudes fails a polarization-frame test, and that rebuilding it at detector level gives a nonzero but extremely small timing-residual…
The reading
What carries the argument
The central object is the detector-projected gravitational-wave amplitude $\tilde{h}_A(\vec{k},\hat{e}) = h_A(\vec{k})\,d_{ij}(\hat{e})\,e^A_{ij}(\hat{k})$, where $d_{ij}(\hat{e}) = \hat{e}_i\hat{e}_j$ is the pulsar detector tensor. Under a local rotation of the polarization basis, the helicity amplitude and the polarization tensor acquire opposite spin-dependent phases, so their product is invariant; imposing the quadratic non-Gaussian ansatz on this projected amplitude cancels the phases explicitly in the three-point correlation. The quantitative evaluation uses the simplified response $F^A \sim d_{ij} e^A_{ij}/2$, retains only the Earth term, and reduces the angular integrals to the analytic kernel $\Upsilon(f,f',f'')$ built from $J_{\ell\ell 0}(f,f',f'')$ with triangle support $β(\Delta)$. This machinery is what turns the polarization-frame problem into a concrete, evaluable single-pulsar skewness.
Extended reading notes
Core claim
The central claim is twofold. First, the helicity-space ansatz (33–34), in which the nonlinear part is a quadratic convolution of Gaussian helicity amplitudes, does not produce a manifestly polarization-frame-invariant timing-residual bispectrum for tensor and vector modes: under a local polarization-basis rotation the quadratic convolution lacks the spin-dependent compensating phase required by the response coefficients, so the bispectrum (53) depends on the basis choice. Second, defining the same quadratic structure on the detector-projected amplitude $\tilde{h}_A(\vec{k},\hat{e}) = h_A(\vec{k})\, d_{ij}(\hat{e})\, e^A_{ij}(\hat{k})$ makes the phase cancellation explicit, and the resulting one-pulsar timing-residual third moment is nonzero and polarization-frame independent. With an isotropic, unpolarized tensor background, equal helicity powers, and the paper's simplified Earth-term response, the benchmark gives a third moment of order $10^{-36}\,\mathrm{s}^3$ and a Gaussian-rms-normalized magnitude $|K|$ of order $10^{-15}$ for $A_{\rm gw}=10^{-15}$ and $g_{\rm NL}=1$, with $K\propto g_{\rm NL}A_{\rm gw}$. The paper stresses that this is an approximate, detector-channel-level phenomenological benchmark, not a generic prediction for an isotropic tensor stochastic background.
Load-bearing premise
The nonzero benchmark skewness stands on treating the detector-projected nonlinear channels (61) as a physical description of non-Gaussian gravitational-wave statistics, which the paper explicitly says it has not shown; the magnitude also rests on the Earth-term-only, simplified-antenna-response approximations, whose corrections to the cubic statistic are not calculated.
Editorial extensions
If this is right
- The simple quadratic helicity ansatz (33–34) is not usable as written for PTA three-point statistics; consistent tensor and vector non-Gaussian models must carry compensating spin phases.
- Non-Gaussian timing-residual statistics can be formulated at detector level, where polarization-frame invariance is automatic, and that formulation predicts a nonzero single-pulsar skewness.
- At unit coupling and $A_{\rm gw}=10^{-15}$, the third moment is about $10^{-36}\,\mathrm{s}^3$ and $|K|\sim10^{-15}$, scaling as $K\propto g_{\rm NL}A_{\rm gw}$.
- Because the two-point function is unchanged at $\mathcal{O}(g_{\rm NL})$, a measured third moment would probe phase-coupled non-Gaussian information absent from the power spectrum.
- Within the illustrated band, redder spectral slopes increase the magnitude of the third moment.
Reading between the lines
- Editorial inference: Restoring the full antenna response and pulsar terms may change the third moment by more than an order of magnitude, so the reported $10^{-36}\,\mathrm{s}^3$ and $|K|\sim10^{-15}$ should be read as an order-of-magnitude benchmark rather than a sharp prediction.
- Editorial inference: The frame-invariance failure likely attaches to any scalar-inspired quadratic convolution of helicity amplitudes, not just the one written here, so spin-covariant structure is a general requirement for tensor and vector non-Gaussian SGWB models.
- Editorial inference: A single-pulsar skewness this small will be hard to detect; higher connected moments such as the four-point function, or multi-pulsar bispectrum combinations, may be better observational targets even when the bispectrum is zero.
- Editorial inference: The detector-projected construction could be adapted to vector and scalar polarizations in alternative theories of gravity, where the same polarization-phase issue appears.
Original abstract
Motivated by searches for weak primordial non-Gaussianity in the cosmic microwave background, we investigate perturbative departures from Gaussian stochastic gravitational wave background statistics using pulsar timing arrays. We consider a quadratic ansatz that shows three-point correlations at first order in the nonlinear coupling and covariance corrections at second order. We first examine the ansatz for gravitational wave helicity amplitudes. However, the timing residual bispectrum that this produces is not manifestly invariant under local polarization basis rotations for tensor and vector modes. Then, we formulate a polarization-frame-invariant phenomenological model for detector-projected amplitudes and evaluate its one-pulsar timing residual third moment with a simplified response. Our numerical benchmarks have shown that this gives a very small Gaussian-rms-normalized skewness.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
Motivated by searches for weak primordial non-Gaussianity in the cosmic microwave background, we investigate perturbative departures from Gaussian stochastic gravitational wave background statistics using pulsar timing arrays. We consider a quadratic ansatz that shows three-point correlations at first order in the nonlinear coupling and covariance corrections at second order. We first examine the ansatz for gravitational wave helicity amplitudes. However, the timing residual bispectrum that this produces is not manifestly invariant under local polarization basis rotations for tensor and vector modes. Then, we formulate a polarization-frame-invariant phenomenological model for detector-projected amplitudes and evaluate its one-pulsar timing residual third moment with a simplified response. Our numerical benchmarks have shown that this gives a very small Gaussian-rms-normalized skewness.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | Motivated by searches for weak primordial non-Gaussianity in the cosmic microwave background, we investigate perturbative departures from Gaussian stochastic gravitational wave background statistics using pulsar timing arrays. | contribution 0.51 | 0.62 | 0.05 | asserted 0.61 |
| s1 | We consider a quadratic ansatz that shows three-point correlations at first order in the nonlinear coupling and covariance corrections at second order. | method 0.90 | 0.83 | 0.06 | asserted 1.00 |
| s2 | We first examine the ansatz for gravitational wave helicity amplitudes. | method 0.70 | 0.82 | 0.02 | asserted 0.50 |
| s3 | However, the timing residual bispectrum that this produces is not manifestly invariant under local polarization basis rotations for tensor and vector modes. | limitation 0.86 | 0.78 | 0.07 | asserted 1.00 |
| s4 | Then, we formulate a polarization-frame-invariant phenomenological model for detector-projected amplitudes and evaluate its one-pulsar timing residual third moment with a simplified response. | method 0.86 | 0.95 | 0.08 | asserted 0.99 |
| s5 | Our numerical benchmarks have shown that this gives a very small Gaussian-rms-normalized skewness. | result 1.00 | 0.95 | 0.41 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:2610.02735 by the authors.
References
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.