{"id":"e7744ae5-cfc9-45d5-b5a9-a5f1f961e9cb","arxiv_id":"1908.08400","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors present a new, independent derivation of the second-order relativistic galaxy number count in a flat expanding universe, but they only validate it at linear order against previous work.","lead":"This paper derives a new mathematical expression for how the number of galaxies counted in a survey is affected by the warping of space and time, up to second order in the standard cosmological perturbation expansion. It is an independent cross-check of existing calculations that currently disagree, and it matters for interpreting the next generation of galaxy surveys.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 4.23 drops the background expansion term in the linearized shear equation; the |σ|^2 term feeding Eq. 4.21 is therefore wrong, invalidating Eq. 5.17 as written.","rationale":"The reader's weakest_assumption focuses on the statement that δ(1)Σ|_o=0 is 'without loss of generality'. That boundary condition is actually the standard regular-caustic condition for a bundle focused at the observer, so it is unlikely to be the real weakness. The more serious problem is the integrated solution: Eq. 4.22 is presented as if the first-order shear obeys dδΣ/dλ = S, omitting the −θbar δΣ term from the linearized Sachs equation. In the conformal Minkowski background, θbar = −2/(λo−λ), so the solution necessarily contains a (λo−λ)^{-2}∫(λo−λ)^2 S weighting. Without that weighting, the δ(1)Σδ(1)Σ contribution to Eq. 4.21 is incorrect, and this error propagates into δ(2)d_A (Eq. C9) and then into the headline result Eq. 5.17. The missing second-order comparison with the literature is a separate validation gap, but the shear-integration issue is a concrete algebraic defect that can be checked independently. The use of xAct/xPand does not protect against a missing physical term in the equations fed into the algebra, and the paper provides no machine-checked proof. If the concrete test confirms the factor of 9 discrepancy, the central claim that Eq. 5.17 is the full second-order galaxy number count cannot stand as written.","tokens_in":34442,"tokens_out":24261,"duration_ms":268548,"concrete_test":"Take a constant trace-free tidal tensor S_ij, e.g. Φ1+Ψ1 = (1/2)A_ij x^i x^j with A_ii=0. Compute δΣ_ij from Eq. 4.23 and from the exact solution of the linearized Sachs equation including θbar = −2/(λo−λ), both with δΣ(λo)=0. The correct result is δΣ_ij(λs) = −A_ij(λo−λs)/3 (up to the sign convention), while Eq. 4.23 gives A_ij(λo−λs); the contraction δΣ_ij δΣ^ij entering Eq. 4.21 differs by a factor 9. Equivalently, derive the first-order Jacobi matrix J = (λo−λ)I + δJ from d^2δJ/dλ^2 = −K Jbar and extract the optical shear σ = (J' J^{-1})_{traceless}; compare that with Eq. 4.23 to settle whether the missing kernel is present.","verdict_should_be":"REJECT","load_bearing_attack":"In Section IV, the first-order null shear is obtained from dδ(1)Σ_ij/dλ = (1/2)δ_ij∇^2(Φ1+Ψ1) − (Φ1+Ψ1)_{,ij}, then integrated with δ(1)Σ|_o=0 to give δ(1)Σ_ij = ∫_{λo}^{λs} S_ij dλ (Eq. 4.23). This is not the linearization of the Sachs shear equation (4.6) in the conformal Minkowski background. For a bundle focused at the observer the background expansion is θbar = −2/(λo−λ), so the linearized equation is dδΣ_ij/dλ = −θbar δΣ_ij + S_ij = [2/(λo−λ)]δΣ_ij + S_ij. The regular solution with δΣ(λo)=0 is δΣ_ij(λs) = −(λo−λs)^{-2} ∫_{λs}^{λo} dλ (λo−λ)^2 S_ij(λ), not ∫ S_ij dλ. For a constant trace-free S, the correct δΣ is −S(λo−λs)/3, whereas Eq. 4.23 gives S(λo−λs); the squared contraction entering Eq. 4.21 differs by a factor 9. Since Eq. 4.21 is the origin of the δ(1)Σδ(1)Σ term in the second-order angular diameter distance, the error propagates through Eq. C9 into the final galaxy overdensity Eq. 5.17. The boundary condition δΣ|_o=0 is itself the regular-caustic condition, so the reader's concern is better directed at the missing θbar coupling in the integrated solution.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":34831,"tokens_out":6317,"duration_ms":61147,"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":[{"comment":"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":"Section IV, Eqs. (4.22)-(4.24)"},{"comment":"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.","section":"Section VII and Eq. (5.17)"}],"minor_comments":[{"comment":"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.","section":"Section IV, near Eq. (4.23)"},{"comment":"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.","section":"Eq. (5.17)"},{"comment":"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.","section":"Appendix C, Eq. (C9)"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern is correct and is a genuine technical error: Eq. (4.22) omits the background expansion coupling in the linearized Sachs shear equation, and the error affects the central result. The approach is salvageable and the linear-order comparison is useful, so I recommend major revision rather than rejection, but the authors must recompute the shear and all downstream expressions and should add an actual second-order comparison before the paper can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a careful, independent derivation of relativistic galaxy number counts up to second order, and the linear-order part checks out. But the central second-order result looks wrong as written, for a specific reason: the first-order shear used in the distance calculation is not the solution of the linearized Sachs equation.\n\nWhat the paper does well: it works with the affine parameter and a conformal map to Minkowski, computes the redshift, area distance, volume, and assembles the number overdensity Delta_g. The first-order expression is compared with Di Dio, Bertacca, and Yoo, and matches after translating notation. That is real work and a useful sanity check. The paper is also honest that it is not the first second-order calculation and that comparison at second order is deferred to a companion paper.\n\nThe soft spot: in Section IV, Eq. (4.22) gives a source S_ij for the first-order shear, and Eq. (4.23) integrates it as deltaSigma_ij = integral S_ij dlambda with deltaSigma|_o=0. That integration ignores the background expansion theta_bar = -2/(lambda_o-lambda) of a bundle focused at the observer. The linearized Sachs equation is d(deltaSigma_ij)/dlambda = -theta_bar deltaSigma_ij + S_ij, i.e. d(deltaSigma_ij)/dlambda = [2/(lambda_o-lambda)]deltaSigma_ij + S_ij. The regular solution with deltaSigma(lambda_o)=0 is deltaSigma_ij(lambda_s) = -(lambda_o-lambda_s)^(-2) integral_{lambda_s}^{lambda_o} dlambda (lambda_o-lambda)^2 S_ij(lambda), not the bare integral in Eq. (4.23). The squared shear entering Eq. (4.21), hence Eq. (C9) and finally Eq. (5.17), is therefore wrong. For a constant trace-free source the discrepancy is a factor of 9. The boundary condition itself is fine, since it is the regular caustic condition, so the reader's worry about the boundary condition is less sharp than the missing theta_bar coupling.\n\nBecause the shear-squared term is one of the second-order sources for the angular diameter distance, and because no second-order comparison is made with any existing result, the main claim of the paper is not supported as written. The paper is not a mess; the route is coherent and the linear-order parts are solid, but the advertised result needs revision.\n\nWho is this for? People working on relativistic corrections to large-scale structure, and would-be readers of the companion paper. It deserves a serious referee, since this is exactly the kind of calculation where referees earn their keep, but the referee should be asked to check the Sachs equation step carefully. I would not cite Eq. (5.17) until the shear issue is fixed and a second-order comparison appears.","headline":"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.","tokens_in":35347,"tokens_out":4590,"would_cite":false,"duration_ms":46068,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.Cq"],"model":"deepseek-v4-flash","headline":"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…","keywords":["galaxy number counts","cosmological perturbation theory","second order","angular diameter distance","Sachs equation","lensing","large-scale structure","relativistic redshift"],"falsifier":"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.","tokens_in":34240,"feed_emoji":"🌌","tokens_out":8068,"duration_ms":74785,"temperature":0.7,"pith_summary":"This paper computes the observed galaxy number-count overdensity up to second order in cosmological perturbation theory, an expression long and intricate enough that earlier groups disagree at this order. The authors build the observable from first principles in a flat universe filled with pressureless matter: the second-order redshift perturbations, the angular diameter distance from the Sachs optical equations, and the physical volume element of a survey. They carry out the distance calculation in a conformally mapped perturbed Minkowski background and convert final results from affine parameter to observed redshift. The main result is Eq. (5.17), $\\Delta_g^{(2)}$, which, if correct, gives the complete relativistic correction up to second order and provides an independent cross-check of previously published second-order expressions.","feed_headline":"Relativistic galaxy counts reach full second order","feed_subtitle":"Redshift, distance, and survey volume are derived together for percent-level cosmology.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the covariant volume element $dV = -E\\,d_A^2\\,d\\lambda\\,d\\Omega$ used to assemble the number counts.","marker":"[25]"},{"why":"Provides the Sachs propagation equations for null bundle expansion and shear that determine the angular diameter distance.","marker":"[29]"},{"why":"Gives the Kristian\\,–\\,Sachs series expansion for $d_A^2$ that fixes the boundary conditions at the observer.","marker":"[31]"},{"why":"The second-order number-count result that the paper uses as its benchmark for linear-order agreement and eventual second-order comparison.","marker":"[15]"},{"why":"The cosmic-rulers formulation whose linear-order number counts must be reproduced after gauge and sign translation.","marker":"[13]"},{"why":"The luminosity-distance route to galaxy overdensities with which the paper's linear-order expression is checked.","marker":"[17]"},{"why":"The planned follow-up in which the full second-order comparison with earlier results is to be performed.","marker":"[37]"}],"fun_headline_variants":["Second-order galaxy counts, derived independently","Full second-order galaxy counts via an independent route","Galaxy counts at second order, independently derived","Second-order galaxy number counts, freshly derived","Redshift, distance, volume: full second-order galaxy counts"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Second-order galaxy counts, derived independently","Full second-order galaxy counts via an independent route","Galaxy counts at second order, independently derived","Second-order galaxy number counts, freshly derived","Redshift, distance, volume: full second-order galaxy counts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000629,"raw_usage":{"total_tokens":2874,"prompt_tokens":879,"completion_tokens":1995,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":495,"completion_tokens_details":{"reasoning_tokens":1924}},"tokens_in":495,"tokens_out":1995,"duration_ms":15058,"temperature":1.0,"reasoning_tokens":1924,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:42:06.799419+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Higher order relativistic galaxy number counts: dominating terms","cited_arxiv_id":"1606.02113","evidence_quote":"Provides the Sachs propagation equations for null bundle expansion and shear that determine the angular diameter distance."},{"cited_title":"Challinor and A","cited_arxiv_id":null,"evidence_quote":"The luminosity-distance route to galaxy overdensities with which the paper's linear-order expression is checked."},{"cited_title":"cosmic rulers","cited_arxiv_id":null,"evidence_quote":"The planned follow-up in which the full second-order comparison with earlier results is to be performed."}],"review_version":1}