REVIEW 2 major objections 3 minor 2 cited by
Galaxy number counts at second order: an independent approach
T0 review · 2 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper derives the full second-order relativistic galaxy number-count overdensity in a flat, matter-dominated FLRW universe, assembling the redshift, angular diameter distance, and survey volume corrections needed for percent-level…
desk verdict The derivation is serious and the linear-order check is good, but Eq. (4.23) drops the background expansion term in the linearized Sachs shear equation, so the advertised second-order number counts in Eq. (5.17) are not reliable as written. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing machinery is the Sachs propagation equations for a bundle of null geodesics, which relate the expansion of the bundle to the Ricci and Weyl terms and give the angular diameter distance through $d^2 d_A/d\lambda^2 = -\frac12(R_{\mu\nu}k^\mu k^\nu + \Sigma_{\mu\nu}\Sigma^{\mu\nu})d_A$. The paper supplements these with the Kristian\,–\,Sachs series expansion of $d_A^2$, which supplies the boundary conditions $d_A(\lambda_o)=0$ and $d_A'(\lambda_o)=-E_o$. To simplify integrations, a conformal transformation maps the perturbed FLRW metric to a perturbed Minkowski metric, where the background distance is simply $\lambda_o-\lambda_s$ and the area distance maps back as $\hat{d}_A=a\,d_A$. The volume element $dV=-E\,d_A^2\,d\lambda\,d\Omega$ then combines energy and distance, and perturbative inversion of the redshift\,–\,affine-parameter relation converts the final result into redshift space.
What would settle it
Evaluate the paper's Eq. (5.17) side by side with the independent second-order number-count results of Refs. [13], [14], [15], and [17] in an Einstein\,–\,de Sitter universe, translating all expressions to a common gauge and notation; any residual second-order difference would show that the claimed full expression is not gauge-complete. Alternatively, recompute the angular diameter distance in a model with a deliberately non-zero first-order shear at the observer and check whether $\Delta_g^{(2)}$ shifts.
Extended reading notes
Core claim
The central claim is that Eq. (5.17) is the full second-order galaxy number-count overdensity $\Delta_g^{(2)}$ in a flat FLRW universe with pressureless matter and scalar perturbations, allowing non-zero anisotropic stress. The derivation proceeds along an independent route: instead of computing luminosity distance, the authors use the volume measure $dV = -E\,d_A^2\,d\lambda\,d\Omega$, with the angular diameter distance $d_A$ obtained from the Sachs propagation equations. After a conformal transformation to a Minkowski background, they integrate the null geodesic equations, obtain the second-order redshift, distance, and volume, and then invert the affine parameter to observed redshift so that the overdensity is expressed in terms of observable quantities. The paper verifies that the linear-order limit of its expression agrees exactly with the established results of Refs. [13], [14], [15], [17], and the linear-order literature, while the full second-order comparison is left to a follow-up paper.
Load-bearing premise
The derivation assumes that setting the first-order shear perturbation of the light bundle to zero at the observer's position is a free choice, so the integrated shear in the Sachs equation carries no observer-side boundary term; if that choice is not genuinely general, the final number-count expression changes.
Editorial extensions
If this is right
- The expression can be evaluated directly with metric potentials from field equations or simulations, giving predicted number counts as a function of line-of-sight position for future surveys.
- Surveys reaching percent-level precision can test the relativistic corrections beyond linear order, since lensing, Doppler, and integrated Sachs\,–\,Wolfe contributions now appear consistently at second order.
- Because the derivation is independent of the luminosity-distance route used by other groups, it provides a cross-check that can isolate where earlier second-order disagreements originate.
- The result is expressed directly in observed redshift and comoving distance, so it can be implemented in survey analysis pipelines without further gauge translation.
Reading between the lines
- Inference: the validity of setting the first-order shear at the observer to zero is testable; a non-zero boundary value would alter the integrated shear entering the angular diameter distance and hence Eq. (5.17), so the paper's claim of general validity rests on this choice being truly free.
- Inference: the deferred comparison against the independent second-order results of Refs. [13], [14], [15], and [17] is the decisive checkpoint, because if those expressions do not match after a common gauge and variable translation, the difference will point to whether the redshift, distance, volume, or $\lambda$-to-$z$ conversion step contains the discrepancy.
- Inference: the affine-parameter formulation suggests a numerical evaluation along simulated light cones without first inverting $\lambda(z)$, which would yield a direct observational test of the second-order corrections.
- Inference: one could construct a toy model with a controlled non-zero observer shear to see how much $\Delta_g^{(2)}$ shifts, providing a quantitative measure of the boundary condition's impact on next-generation survey predictions.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents an independent derivation of the relativistic galaxy number-count overdensity to second order in cosmological perturbation theory. Working in a flat FLRW universe with pressureless matter and scalar perturbations in the longitudinal gauge, the authors compute the perturbed null geodesics, the observed redshift, the angular diameter distance (via a conformal transformation to a perturbed Minkowski spacetime), and the survey volume element, and combine these into the second-order galaxy overdensity Delta_g^(2) in Eq. (5.17). The linear-order result is compared with the works of Di Dio et al., Bertacca et al., and Yoo and Zaldarriaga, and found to agree. A comparison at second order is explicitly deferred to a companion paper, Ref. [37]. The tensor algebra was performed with xAct and xPand.
Significance. If the second-order expression were correct, the paper would provide a useful independent derivation of a quantity needed for next-generation surveys, and the use of the affine parameter as the integration variable is a genuine methodological difference from most earlier work. The linear-order agreement with three independent groups is a strong consistency check, and the derivation is self-contained rather than fitted to any target. However, the central second-order result is not verified against any existing calculation, and the derivation contains a serious error in the linearized null-shear equation that propagates directly into Eq. (5.17). The paper's contribution therefore remains potentially valuable but is not reliable as it stands.
major comments (2)
- [Section IV, Eqs. (4.22)-(4.24)] The linearized Sachs shear equation is missing the background expansion term. In the conformal Minkowski background, a bundle focused at the observer has theta_bar = -2/(lambda_o - lambda), so linearizing Eq. (4.6) gives d delta^(1)Sigma_ij/dlambda = -theta_bar delta^(1)Sigma_ij + S_ij = [2/(lambda_o - lambda)] delta^(1)Sigma_ij + S_ij, not d delta^(1)Sigma_ij/dlambda = S_ij as written in Eq. (4.22). The regular solution with delta^(1)Sigma_ij(lambda_o)=0 is delta^(1)Sigma_ij(lambda_s) = -(lambda_o - lambda_s)^{-2} integral_{lambda_s}^{lambda_o} (lambda_o - lambda)^2 S_ij(lambda) dlambda, which for constant trace-free S_ij equals -S_ij (lambda_o - lambda_s)/3 rather than the S_ij (lambda_o - lambda_s) of Eq. (4.23). The squared contraction in Eq. (4.24) is therefore wrong by a factor of 9. Since delta^(1)Sigma delta^(1)Sigma enters the second-order area-distance equation (4.21), the error propagates through Eq. (C9), the volume element (C10), and the main result Eq. (5.17). All downstream expressions must be recomputed with the corrected shear.
- [Section VII and Eq. (5.17)] The central second-order result is not compared with any existing second-order calculation; the authors state in Section VII that a full comparison is left to the companion paper, Ref. [37]. Given that the paper itself notes that previous second-order results are in disagreement, and given the shear error identified above, Eq. (5.17) is currently an unverified expression at exactly the order at which the paper claims novelty. The authors should provide the second-order comparison, or at minimum a check in a simplified limit (for example Einstein-de Sitter with no anisotropic stress), before the main claim can be reliably assessed.
minor comments (3)
- [Section IV, near Eq. (4.23)] The phrase 'Without loss of generality, we set the perturbation of the shear at the observer delta^(1)Sigma_mu_nu|_o = 0' is misleading: for a bundle focused at the observer, this is the regular caustic boundary condition, not a free gauge choice. With the corrected linearized equation, this condition selects the unique regular solution.
- [Eq. (5.17)] After the change of variable to the comoving distance chi, several integrals still have upper limit lambda_s (for example the terms written as integral_0^{lambda_s} dchi [Phi_1 (dPhi_1/dzeta - 2 dPsi_1/dzeta)]); the upper limit should be chi_s for dimensional consistency.
- [Appendix C, Eq. (C9)] The nested integrals in Eq. (C9) are extremely difficult to verify by eye; providing the xAct/xPand notebook or an intermediate derivation step for the shear contribution would substantially help referees and readers.
Circularity Check
No significant circularity: the second-order galaxy overdensity is derived directly from standard geodesic, Sachs, and volume relations, with no fitted parameters and no load-bearing self-citation.
full rationale
The derivation is self-contained. The central result, Eq. (5.17), is obtained by substituting previously derived pieces, Eq. (C9) for the second-order angular diameter distance, Eq. (C10) for the second-order volume perturbation, Eq. (C11) for the redshift-space density perturbation, and the second-order redshift (C8), into the exact identity Δ_g = δ_z + δV/V (5.7)–(5.9). Each of these pieces is computed from the perturbed metric (2.1), the null geodesic equations (3.4)–(3.5), (C1)–(C4), and the Sachs equations (4.5)–(4.7). No parameter is fitted to data, and no target expression is used as an input. The only self-referential element is Ref. [37], the planned companion paper, cited twice (Section VI and Section VII) for the deferred second-order comparison; this comparison is not used anywhere in the derivation, so it is not load-bearing. The first-order result is checked against the independent external results of Di Dio et al., Bertacca et al., and Yoo & Zaldarriaga (Section VI), which confirms that the pipeline reproduces a known result rather than being defined to match it. The boundary condition δ(1)Σ_μν|_o = 0 after Eq. (4.22) is an assumption about the observer shear, and the skeptic's complaint that the integrated shear solution in Eq. (4.23) omits the background expansion term in the linearized Sachs equation is a physical-correctness concern, not a circularity: the final expression is not wired to equal that assumption by construction. Likewise, the acknowledged absence of a full second-order comparison to other groups (explicitly deferred to the companion paper) is a completeness limitation and not an instance of the derivation reducing to its own inputs.
Assumptions & free parameters
assumptions (6)
- domain assumption FLRW background with flat spatial curvature and a single pressureless matter component (rho proportional to (1+z)^3)
- domain assumption Only scalar metric perturbations in longitudinal gauge; vector and tensor modes are neglected
- standard math Photon propagation follows null geodesics and the Sachs optical equations
- standard math The conformal transformation maps null geodesics and sets hat_d_A = a d_A and d_hat_lambda = a^2 d_lambda
- ad hoc to paper The shear perturbation at the observer is set to zero without loss of generality (delta^(1)Sigma_mu_nu|o = 0)
- domain assumption Galaxy bias and evolution bias are omitted; the computed quantity is the matter density perturbation, not the biased galaxy density
Cite this review
Pith. "Pith review of Galaxy number counts at second order: an independent approach." pith.science (2026). https://pith.science/paper/L7YECPTK
@misc{pith2026190808400,
author = {Pith},
title = {Pith review of: Galaxy number counts at second order: an independent approach},
year = {2026},
howpublished = {\url{https://pith.science/paper/L7YECPTK}},
note = {Machine review of arXiv:1908.08400}
}
read the original abstract
Next generation surveys will be capable of determining cosmological parameters beyond percent level. To match this precision, theoretical descriptions should look beyond the linear perturbations to approximate the observables in large scale structure. A quantity of interest is the Number density of galaxies detected by our instruments. This has been focus of interest recently, and several efforts have been made to explain relativistic effects theoretically, thereby testing the full theory. However, the results at nonlinear level from previous works are in disagreement. We present a new and independent approach to computing the relativistic galaxy number counts to second order in cosmological perturbation theory. We derive analytical expressions for the full second order relativistic observed redshift, for the angular diameter distance and for the volume spanned by a survey. Finally, we compare our results with previous works which compute the general distance-redshift relation, finding that our result is in agreement at linear order.
Figures
Forward citations
Cited by 2 Pith papers
-
Observable Gravitational Wave Strain at Second Order
At second order, the gravitational-wave strain measured by geodesic observers exchanging light pulses is the transverse-traceless metric perturbation in the Newton gauge (h_N^(2)).
-
Tetrad formalism for exact cosmological observables
The paper derives exact, coordinate-independent equations for cosmological observables on a new 'observer space-time' manifold, using tetrads to track the observer frame.
Reference graph
Works this paper leans on
-
[37]
we perform a comparison of the leading terms of the second order expansion of the galaxy number counts. 20 Our result, as given in Eq. (5.7) is ∆(1) g (ni,z ) = [δ(1)ρ(ni,z ) ¯ρ(z) + 3Φ1 ] − 2 (Φ1 + Ψ1) + 1 H ( Ψ1 ′−∂χΦ1 + d ( v1ini) dς ) + (H′ H2 + 2 Hχ )[ Φ1− ( v1ini) + ∫ χ 0 d˜χ ( Φ1 ′ + Ψ1 ′)] − 4 χ ∫ χ 0 d˜χΨ1 − 1 χ ∫ χ 0 d˜χ (˜χ−χ) ˜χ [ ∇2 (Φ1 + Ψ1)...
work page 2016
-
[1]
Geodesic Equation Solving Eq. (2.8) at second order gives dδ(2)ν dλ =−1 2 [dΦ2 dλ + dΨ2 dλ ] + 2dδ(1)n i dλ δ(1)ni− 2dδ(1)ν dλ δ(1)ν− dδ(2)n i dλ ni (C1) − 4dδ(1)ν dλ [ Φ1 + Ψ1 ] − 4δ(1)ν [dΦ1 dλ + dΨ1 dλ ] − 4Φ1 [dΦ1 dλ + dΨ1 dλ ] − 4dΦ1 dλ [ Φ1 + Ψ1 ] , dδ(2)n i dλ = 2dΨ2 dλ ni− [ Φ2 i , + Ψ2 i , ] − 4Φ1Ψ1 ′ni + 4δ(1)n i Ψ1 ′ + 4Φ1Φ1 i , (C2) − 2δ(1)n j...
-
[2]
Energy The perturbed energy in terms of the metric potentials is given by δ(1)E = ( −2Φ1 ⏐⏐⏐ s o + ∫ λs λo ( Φ1 ′ + Ψ1 ′) dλ ) + Φ1−v1ini, (C5) δ(2)E = 1 2 [ 2 {( −2Φ1 ⏐⏐⏐ s o + ∫ λs λo ( Φ1 ′ + Ψ1 ′) dλ ) + Φ1−v1ini } Φ1 (C6) − Φ1 2− 5Ψ1 2− 6Φ1Ψ1 + 1 2 (Φ2− Ψ2)−niv2i − 4 (Φ1 + Ψ1) Φ1 ⏐⏐⏐ s o + 2 (Φ1 + Ψ1) ∫ λs λo ( Φ1 ′ + Ψ1 ′) dλ − 4niv1iΨ1 ⏐⏐⏐ s o + 2v...
-
[3]
Observed Redshift At second order the redshift is δ(2)z = [ Φ1 2 + ( v1ini)2 − 2Φ1 ( v1ini)] o − Φ1|oδ(1)ν|s +δ(1)ν|s ( v1ini) o− Φ1|sΦ1|o (C7) + Φ1|s ( v1ini) o + ( v1ini) s Φ1|o− ( v1ini) s ( v1ini) o + 1 2 [ 2δ(1)EΦ1− Φ1 2− 5Ψ1 2 − 6Φ1Ψ1 +δ(2)ν + 1 2 (Φ2− Ψ2) + 2δ(1)ν (Φ1 + Ψ1)− 2δ(1)n i v1i−niv2i ] s − 1 2 [ Φ1 2− 5Ψ1 2− 6Φ1Ψ1 + 1 2 (Φ2− Ψ2)− (v1ini)Φ...
-
[4]
Angular Diameter Distance Using Eqs. (3.4), (3.5), (C3), (C4), (C6), (4.18), (B1), (B2), (B3) and (4.24) we find that the second order pertur- bation to the angular diameter distance becomes δ(2)dA(λs) ¯dA(λs) = 1 2 [ Φ1 2− 5Ψ1 2− 6Φ1Ψ1 + 1 2 (Φ2− Ψ2)− (v1ini)Φ1−niv2i ] o (C9) − 1 λo−λs ∫ λs λo dλ ∫ λs λs d˜λ { 4 ( −2Φ1 ⏐⏐⏐ s s + ∫ λs λo ( Φ1 ′ + Ψ1 ′) dλ ...
-
[5]
Physical Volume The second order perturbation to the physical volume is dδ(2)V =− ¯E ¯d2 A [(δ(1)dA ¯dA )2 + (δ(1)E ¯E )(δ(1)dA ¯dA ) + δ(2)dA ¯dA + 1 2 δ(2)E ¯E ] dλdΩ (C10) =a2 (λs) (λs−λo)2 {( Φ1|o− (v1ini)o− Φ1|s o− 3 2 (Ψ1− Φ1) ⏐⏐⏐ s o + ∫ λs λo dλ ∫ λs λs d˜λ { (Φ1 + Ψ1)′′−∇ 2 (Φ1 + Ψ1) } − 1 λo−λs [ 2 ∫ λs λo dλ ∫ λs λs d˜λ (Φ1 + Ψ1)′− ∫ λs λo dλ ∫...
-
[6]
Redshift Density Using Eqs. (3.14), (C7) and (5.15) in Eq. (5.10) we find that the redshift density perturbation at second order is given by δ(2) z (ni,z ) = 1 2 δ(2)ρ(ni,z ) ¯ρ(z) + 3 2(1 + ¯z)δ(2)z(ni,z ) + 3 (1 +z)2 [ δ(1)z(ni,z ) ]2 (C11) = 1 2 δ(2)ρ(ni,z ) ¯ρ(z) + 3 2(1 + ¯z) ([ Φ1 2 + ( v1ini)2 − 2Φ1 ( v1ini)] o − Φ1|o [ −2Φ1 + ∫ λs λo ( Φ1 ′ + Ψ1 ′)...
-
[7]
Green et al., ArXiv e-prints (2012), 1208.4012
J. Green et al., ArXiv e-prints (2012), 1208.4012
arXiv 2012
Show all 45 references
-
[8]
The BOSS Collaboration, The Astronomical Journal 145, 10 (2013), 1208.0022. 35
2013 arXiv
-
[9]
Zhao et al., Mon
G.-B. Zhao et al., Mon. Not. Roy. Astron. Soc. 457, 2377 (2016), 1510.08216
2016 arXiv
- [10]
-
[11]
M. G. Santos, M. Cluver, M. Hilton, M. Jarvis, G. I. G. Jozsa, L. Leeuw, O. Smirnov, R. Taylor, F. Abdalla, J. Afonso, et al., Meerklass: Meerkat large area synoptic survey (2017), 1709.06099
2017 arXiv
-
[12]
Maartens, F
R. Maartens, F. B. Abdalla, M. Jarvis, and M. G. Santos, Cosmology with the ska – overview (2015), 1501.04076
2015 arXiv
-
[13]
L. S. Collaboration, P. A. Abell, J. Allison, S. F. Anderson, J. R. Andrew, J. R. P. Angel, L. Armus, D. Arnett, S. J. Asztalos, T. S. Axelrod, et al., Lsst science book, version 2.0 (2009), 0912.0201
2009 arXiv
-
[14]
Bertacca, R
D. Bertacca, R. Maartens, and C. Clarkson, Journal of Cosmology and Astroparticle Physics 2014, 013 (2014), URL http://stacks.iop.org/1475-7516/2014/i=11/a=013
2014
-
[15]
Jeong, F
D. Jeong, F. Schmidt, and C. M. Hirata, Phys. Rev. D85, 023504 (2012), URL https://journals.aps.org/prd/abstract/ 10.1103/PhysRevD.85.023504
2012 doi
-
[16]
Bonvin and R
C. Bonvin and R. Durrer, Phys. Rev. D 84, 063505 (2011), URL https://journals.aps.org/prd/abstract/10.1103/ PhysRevD.84.063505
2011
-
[17]
Challinor and A
A. Challinor and A. Lewis, Phys. Rev. D 84, 043516 (2011), URL https://journals.aps.org/prd/abstract/10.1103/ PhysRevD.84.043516
2011
- [18]
-
[19]
J. Yoo, A. L. Fitzpatrick, and M. Zaldarriaga, Phys. Rev. D 80, 083514 (2009), 0907.0707
2009 arXiv
-
[20]
Bertacca, R
D. Bertacca, R. Maartens, and C. Clarkson, Journal of Cosmology and Astroparticle Physics 2014, 037 (2014), URL http://stacks.iop.org/1475-7516/2014/i=09/a=037
2014
-
[21]
Marozzi, Class
G. Marozzi, Class. Quant. Grav. 32, 045004 (2015), [erratum: Class. Quant. Grav.32,179501(2015)], 1406.1135
2015 arXiv
-
[22]
Di Dio, R
E. Di Dio, R. Durrer, G. Marozzi, and F. Montanari, JCAP 1412, 017 (2014), [Erratum: JCAP1506,no.06,E01(2015)], 1407.0376
2014 arXiv
-
[23]
Yoo, Phys
J. Yoo, Phys. Rev. D 82, 083508 (2010), URL http://link.aps.org/doi/10.1103/PhysRevD.82.083508
2010 doi
- [24]
- [25]
-
[26]
Ben-Dayan, G
I. Ben-Dayan, G. Marozzi, F. Nugier, and G. Veneziano, JCAP 1211, 045 (2012), 1209.4326
2012 arXiv
-
[27]
Fanizza, M
G. Fanizza, M. Gasperini, G. Marozzi, and G. Veneziano, JCAP 1311, 019 (2013), 1308.4935
2013 arXiv
-
[28]
Durrer, The Cosmic Microwave Background (Cambridge University Press, 2008), ISBN 9780521847049, URL https: //books.google.co.uk/books?id=WZkN1QDOD1sC
R. Durrer, The Cosmic Microwave Background (Cambridge University Press, 2008), ISBN 9780521847049, URL https: //books.google.co.uk/books?id=WZkN1QDOD1sC
2008
-
[29]
J. T. Nielsen and R. Durrer, JCAP 1703, 010 (2017), 1606.02113
2017 arXiv
- [30]
- [31]
-
[32]
Ellis, R
G. Ellis, R. Maartens, and M. MacCallum, Relativistic Cosmology, Relativistic Cosmology (Cambridge University Press, 2012), ISBN 9780521381154, URL https://books.google.co.uk/books?id=FPRFi3L0h9kC
2012
-
[33]
K. A. Malik and D. Wands, Physics Reports 475, 1 (2009), ISSN 0370-1573, URL http://www.sciencedirect.com/ science/article/pii/S0370157309000635
2009
-
[34]
O. Umeh, C. Clarkson, and R. Maartens, Classical and Quantum Gravity 31, 205001 (2014), URL http://stacks.iop. org/0264-9381/31/i=20/a=205001
2014
-
[35]
Clarkson, G
C. Clarkson, G. F. R. Ellis, A. Faltenbacher, R. Maartens, O. Umeh, and J.-P. Uzan, Mon. Not. Roy. Astron. Soc. 426, 1121 (2012), 1109.2484
2012 arXiv
-
[36]
R. K. Sachs, Proc. Roy. Soc. Lond. A264, 309 (1961)
1961
-
[39]
Kristian and R
J. Kristian and R. K. Sachs, The Astrophysical Journal 143, 379 (1966)
1966
- [40]
-
[41]
Aretakis, Lecture notes on general relativity
S. Aretakis, Lecture notes on general relativity. Columbia University. (2013), https://web.math.princeton.edu/ ~aretakis/columbiaGR.pdf, URL https://web.math.princeton.edu/~aretakis/columbiaGR.pdf
2013
-
[42]
O. Umeh, C. Clarkson, and R. Maartens, Classical and Quantum Gravity 31, 202001 (2014), URL http://stacks.iop. org/0264-9381/31/i=20/a=202001
2014
- [43]
-
[44]
J. L. Fuentes, J. C. Hidalgo, and K. A. Malik, In preparation (2019)
2019
-
[45]
J. M. Mart´ ın-Garc´ ıa,xAct: Efficient tensor computer algebra for the Wolfram Language (2016), http://www.xact.es/, URL http://www.xact.es/
2016
-
[46]
Pitrou, X
C. Pitrou, X. Roy, and O. Umeh, Class. Quant. Grav. 30, 165002 (2013), 1302.6174
2013 arXiv
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.