Pith. Sign in

REVIEW

Non-linear Pairwise Velocities as a Cosmological Probe

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Mariana Jaber (1), Antonela Taverna (1), Rasmus Strid (1), Wojciech A. Hellwing (2) ((1) National Centre for Nuclear Research, (2) Center for Theoretical Physics, Polish Academy of Sciences)

Cosmology and Nongalactic Astrophysics · Submitted 7 Sep 2026 · 38 references

Pith's verdict: Careful paper that quantifies the information gain from non-linear scales in pairwise velocities; the main limitation is that the model's key derivative (∂aξ̄) is validated only at the fiducial cosmology, but the analysis is otherwise thorough and honest.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 9, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2609.07976 (version 1)
Pith IDIEM4ESGA, Pith's permanent record of this paper and its review
Citation(2) Center for Theoretical Physics, Antonela Taverna (1), Mariana Jaber (1), Polish Academy of Sciences), Rasmus Strid (1), Wojciech A. Hellwing (2) ((1) National Centre for Nuclear Research. "Non-linear Pairwise Velocities as a Cosmological Probe." arXiv:2609.07976 (2026). https://pith.science/paper/IEM4ESGA
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Pairwise velocities on non-linear scales deliver percent-level cosmological constraints in simulation-based validation.

The reading

This paper establishes that the mean pairwise velocity of dark-matter particles, when modelled using the exact pair conservation equation and a stable non-linear power spectrum prescription, carries substantial cosmological information down to separations of 4–8 $h^{-1}\mathrm{Mpc}$. Using the Quijote simulation suite, the authors show that extending the fit from the linear regime ($r_{\rm min}=50\,h^{-1}\mathrm{Mpc}$) into the non-linear regime tightens the marginalized $1\sigma$ uncertainties on $\Omega_{\rm m}$ and $\sigma_8$ by up to 74% and 81% at $z=0.5$, and by 70% and 74% at $z=0$. The uncertainty on $f\sigma_8$ improves by up to 86%, and the joint $\Omega_{\rm m}{-}f\sigma_8$ Figure-of-Merit increases by up to $9.1\times$ relative to the linear baseline. These gains hold at fixed $h$; jointly sampling $h$ broadens the constraints markedly due to a near-degeneracy with $\sigma_8$. The work provides a validated route to harnessing non-linear velocity information from upcoming spectroscopic and kSZ surveys.

What carries the argument

The pair conservation equation $v_{12}(r,a) = -\frac{H(a) r a}{3[1+\xi(r,a)]}\, \partial_a \bar\xi(r,a)$, which links the mean pairwise velocity to the two-point correlation function $\xi$ and its volume-averaged time derivative $\partial_a \bar\xi$. The model is evaluated using HMcode-2020 (via pyCAMB) for the non-linear power spectrum, which provides stable derivatives; the resulting $v_{12}$ predictions are compared to simulation measurements within a multivariate $t$-distribution likelihood that marginalises over covariance sampling noise.

Extended reading notes

Core claim

The central claim is that non-linear pairwise velocities, modelled via the exact pair conservation equation with HMcode-2020 non-linear matter power spectrum, can be used to infer $\Omega_{\rm m}$, $\sigma_8$, and $f\sigma_8$ with much greater precision than previously possible with linear-scale analyses alone. The authors demonstrate this using 100 realisations of the Quijote simulation, showing that by including scales down to $4\,h^{-1}\mathrm{Mpc}$ (at $z=0.5$) the uncertainties shrink dramatically and the joint Figure-of-Merit improves by a factor of $\sim$6–9. They further validate the model against resolution, particle thinning, and box-size effects, and identify the sensitivity of the time derivative of the correlation function as the key numerical challenge, which HMcode-2020 resolves.

Load-bearing premise

The model's accuracy for the scale-factor derivative $\partial_a \bar\xi$ relies on HMcode-2020's prediction of the non-linear power spectrum; although the paper tests stability against the Takahashi prescription, it does not independently verify the derivative's accuracy against direct simulation measurements.

Editorial extensions

If this is right

  • Upcoming surveys (DESI, Euclid, 4MOST/4HS) that measure peculiar velocities or the pairwise kSZ effect can directly apply this non-linear likelihood to extract tighter constraints on growth parameters than linear-only fits allow.
  • The method demonstrates that the scale-dependent response of $v_{12}$ to $\Omega_{\rm m}$ and $\sigma_8$ is distinct from that of the two-point correlation function, offering complementary cosmological information in joint analyses.
  • Within the dark-matter-only idealisation, the constraints reach $\sigma(\Omega_{\rm m})/\Omega_{\rm m}\approx 1.6\%$ and $\sigma(\sigma_8)/\sigma_8\approx 0.5\%$ at $z=0.5$, setting a benchmark for what real-world analyses could achieve once tracer bias and systematics are controlled.
  • The strong degeneracy between $h$ and $\sigma_8$ in $v_{12}$ data implies that pairwise-velocity alone cannot constrain the Hubble constant; external priors or multi-probe combinations are required to break this degeneracy.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the HMcode-2020 fitting function has systematic errors on $k\sim10\,h/\mathrm{Mpc}$ scales at $z\sim0.5$ that are not captured by the comparison with Takahashi, the quoted constraint gains may be overestimated; a direct derivative measurement from high-resolution simulations would settle this.
  • The paper’s fixed-$h$ constraints are conditional on the fiducial cosmology; real analyses that marginalise over $h$ will see much wider posteriors unless combined with distance-scale or CMB data, reducing the effective gain from non-linear scales.
  • Extending the method to biased tracers (halos, galaxies) and redshift space requires modelling the pair-weighting and selection effects; the strong small-scale $v_{12}\to 0$ recovery failure in Quijote at $z=0$ hints that even dark-matter-only simulations may have resolution-dependent systematics that need careful scale cuts.
Original abstract

Peculiar velocities trace gravitational dynamics directly, complementing density-based clustering as wide-field spectroscopic surveys (DESI, Euclid, PFS) and kinetic Sunyaev-Zel'dovich (kSZ) measurements enter the precision era. Observational analyses of pairwise velocities and the pairwise kSZ signal remain largely restricted to linear and quasi-linear scales ($\gtrsim20-30$ $h^{-1}Mpc$), despite substantial cosmological information available at smaller separations. We present a simulation-calibrated likelihood for $v_{12}(r,a)$, built on the exact pair conservation equation and PyCAMB HMcode-2020 non-linear clustering, and validate it against the QUIJOTE and TNG300-3 suites across resolution, particle sampling, box size, and redshift. Using QUIJOTE, we show that extending the fit to non-linear scales tightens the $1\sigma$ uncertainty on $\Omega_\mathrm{m}$ and $\sigma_8$ by $74\%$ and $81\%$ at $z=0.5$ ($r_{\rm min}=4$ $h^{-1}Mpc$), and by $70\%$ and $74\%$ at $z=0$ ($r_{\rm min}=8$ $h^{-1}Mpc$), relative to the linear-regime baseline ($r_{\rm min}=50$ $h^{-1}Mpc$) (all at fixed $r_{max}=140$ $h^{-1}Mpc$); $\sigma(f\sigma_8)$ tightens by up to $86\%$, and the joint $\Omega_\mathrm{m}$-$f\sigma_8$ Figure-of-Merit improves by up to $9.1\times$ relative to the same baseline. These constraints are conditional on $h$ fixed at its QUIJOTE fiducial value, due to a near-degenerate response with $\sigma_8$; jointly sampling $h$ shifts and widens them substantially. This framework provides a validated route to exploiting that information in upcoming direct peculiar-velocity, redshift-space, and kSZ analyses.

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.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 5,614 input tokens.

What the abstract itself states, as probabilities: claims first 0.07 · states limitation 0.97 · human participants 0.01 · controlled comparison 0.87 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.95 · mentions code or data availability 0.04.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Peculiar velocities trace gravitational dynamics directly, complementing density-based clustering as wide-field spectroscopic surveys (DESI, Euclid, PFS) and kinetic Sunyaev-Zel'dovich (kSZ) measurements enter the precision era. Observational analyses of pairwise velocities and the pairwise kSZ signal remain largely restricted to linear and quasi-linear scales ($\gtrsim20-30$ $h^{-1}Mpc$), despite substantial cosmological information available at smaller separations. We present a simulation-calibrated likelihood for $v_{12}(r,a)$, built on the exact pair conservation equation and PyCAMB HMcode-2020 non-linear clustering, and validate it against the QUIJOTE and TNG300-3 suites across resolution, particle sampling, box size, and redshift. Using QUIJOTE, we show that extending the fit to non-linear scales tightens the $1\sigma$ uncertainty on $\Omega_\mathrm{m}$ and $\sigma_8$ by $74\%$ and $81\%$ at $z=0.5$ ($r_{\rm min}=4$ $h^{-1}Mpc$), and by $70\%$ and $74\%$ at $z=0$ ($r_{\rm min}=8$ $h^{-1}Mpc$), relative to the linear-regime baseline ($r_{\rm min}=50$ $h^{-1}Mpc$) (all at fixed $r_{max}=140$ $h^{-1}Mpc$); $\sigma(f\sigma_8)$ tightens by up to $86\%$, and the joint $\Omega_\mathrm{m}$-$f\sigma_8$ Figure-of-Merit improves by up to $9.1\times$ relative to the same baseline. These constraints are conditional on $h$ fixed at its QUIJOTE fiducial value, due to a near-degenerate response with $\sigma_8$; jointly sampling $h$ shifts and widens them substantially. This framework provides a validated route to exploiting that information in upcoming direct peculiar-velocity, redshift-space, and kSZ analyses.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Peculiar velocities trace gravitational dynamics directly, complementing density-based clustering as wide-field spectroscopic surveys (DESI, Euclid, PFS) and kinetic Sunyaev-Zel'dovich (kSZ) measurements enter the precision era. background 1.00 0.04 0.03 asserted 1.00
s1 Observational analyses of pairwise velocities and the pairwise kSZ signal remain largely restricted to linear and quasi-linear scales ($\gtrsim20-30$ $h^{-1}Mpc$), despite substantial cosmological information available at smaller separations. gap 0.93 0.06 0.03 asserted 1.00
s2 We present a simulation-calibrated likelihood for $v_{12}(r,a)$, built on the exact pair conservation equation and PyCAMB HMcode-2020 non-linear clustering, and validate it against the QUIJOTE and TNG300-3 suites across resolution, particle sampling, box size, and redshift. method 0.66 0.97 0.06 asserted 1.00
s3 Using QUIJOTE, we show that extending the fit to non-linear scales tightens the $1\sigma$ uncertainty on $\Omega_\mathrm{m}$ and $\sigma_8$ by $74\%$ and $81\%$ at $z=0.5$ ($r_{\rm min}=4$ $h^{-1}Mpc$), and by $70\%$ and $74\%$ at $z=0$ ($r_{\rm min}=8$ $h^{-1}Mpc$), relative to the linear-regime baseline ($r_{\rm min}=50$ $h^{-1}Mpc$) (all at fixed $r_{max}=140$ $h^{-1}Mpc$); $\sigma(f\sigma_8)$ tightens by up to $86\%$, and the joint $\Omega_\mathrm{m}$-$f\sigma_8$ Figure-of-Merit improves by up to $9.1\times$ relative to the same baseline. result 1.00 0.97 0.74 asserted 1.00
s4 These constraints are conditional on $h$ fixed at its QUIJOTE fiducial value, due to a near-degenerate response with $\sigma_8$; jointly sampling $h$ shifts and widens them substantially. limitation 1.00 0.87 0.06 asserted 1.00
s5 This framework provides a validated route to exploiting that information in upcoming direct peculiar-velocity, redshift-space, and kSZ analyses. future work 0.50 0.37 0.28 asserted 0.92

What kind of paper the classifier reports this is: method or tool 0.97 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Simulation validation of the pairwise velocity model. Top: Convergence of the Quijote v12 (𝑟, 𝑎) measure￾ment at 𝑧 = 0 across particle sampling and mass resolution. Mid￾Resolution measurements at 100% and 10% (production) thinning agree above their respective Nyquist scales, 𝜆 MR Nyq = 3.9 ℎ −1Mpc and 𝜆 MR Nyq,10% = 8.4 ℎ −1Mpc, as do the independent High-Resolution (10%) and TNG300-3 measurements, validating the use of the Mid￾Resolution 10% configuration as the production dataset. Middle: Mean pairwise velocity v12from Quijote Mid-Resolution (10% thinning, production…
Figure 2. Sensitivity of the v12 (𝑟, 𝑎) model to the non-linear power-spectrum prescription and to the discretization of the red￾shift array used in the model evaluation. Top: Ratio of the modelled non-linear matter power spectrum to the Quijote meas￾urement, 𝑃(𝑘) mod/𝑃(𝑘) data, at 𝑧 = 0 and 𝑧 = 0.5, computed with HMcode-2020 (solid) and the Takahashi [21] halofit prescription (dashed). The upper axis shows the corresponding scale 𝑟 = 2𝜋/𝑘. Bottom panels: Ratio 𝑣 mod 12 /𝑣 data 12 for the HMcode-2020 (middle) and Takahashi (lower) prescriptions. Each panel shows 𝑧 = 0 and 𝑧 = 0.5 models,…
Figure 4. shows the two-dimensional posterior distributions at 𝑧 = 0.5, comparing four choices of minimum fitting scale, 𝑟min = 50, 30, 8, and 4 ℎ −1Mpc, at fixed 𝑟max = 140 ℎ −1Mpc. In all four cases the recovered best-fit values remain consistent with the fiducial Quijote cosmology (dashed lines). We take 𝑟min = 50 ℎ −1Mpc as the baseline scale, since this marks the linear regime of streaming motions, where the full non-linear solution Eq. (1) agrees with the linear approximation Eq. (2) to within 5% (Section II A). Progressively extending the fit to smaller, increasingly non￾linear…
Figure 6. and Table III (Appendix C) quantify this mar￾ginal gain directly, taking 𝑟min = 50 ℎ −1Mpc as the baseline. At 𝑧 = 0.5, the 1𝜎 uncertainty on 𝑓 𝜎8 improves by 64% at 𝑟min = 30 ℎ −1Mpc and 88% at 𝑟min = 8 ℎ −1Mpc, settling at 86% at the adopted production scale 𝑟min = 4 ℎ −1Mpc. At 𝑧 = 0, the corresponding improvements are 67% at 𝑟min = 30 ℎ −1Mpc, reaching 85% at the adopted produc￾tion scale 𝑟min = 8 ℎ −1Mpc; as elsewhere in this section, 𝑟min = 4 ℎ −1Mpc is not adopted at 𝑧 = 0 and is therefore omitted from
Figure 7. Sensitivity of the mean pairwise velocity v
Figure 8. Two-dimensional and one-dimensional marginalized posterior distributions for 𝜎8, 𝛺m, and 𝑓 𝜎8 at 𝑧 = 0, comparing three choices of minimum fitting scale, 𝑟min ∈ {50, 30, 8} ℎ −1Mpc, at fixed 𝑟max = 140 ℎ −1Mpc. Shaded regions show the 68% and 95% credible levels; dashed lines mark the fiducial Quijote cosmology. Unlike at 𝑧 = 0.5 (

Showing 6 of the paper's 8 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2609.07976 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.