{"id":"89ef2c98-4e5d-44fe-be90-fdf6d5e10853","arxiv_id":"2502.02037","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A perturbative scheme connects any static spherical density profile to gravitational deflection angles and lensed image positions through the Tolman-Oppenheimer-Volkoff equations.","lead":"This paper derives a perturbative method to compute how light or massive particles are deflected by a static, spherically symmetric matter distribution, starting directly from the density profile. It gives explicit series formulas for the deflection angle and lensed image positions for galaxy halo models including NFW, Hernquist, and power-law profiles.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Term-by-term integration of the small-u expansion is uncontrolled at b≈rm; the density-to-deflection claim needs an error-bound or convergence check.","rationale":"The reader's weakest assumption points to the lack of convergence proofs for the series ansatze and term-by-term integration. I agree that this is the soft spot, but the more specific and load-bearing issue is that the term-by-term integration is performed up to u=1, where the small-u expansion parameter is not small when b is only a few times rm. The paper's own figures use exactly this regime, so the agreement with numerics at fixed truncation order does not settle whether the series is convergent or merely asymptotic with an uncontrolled remainder. The concrete test I propose directly measures that remainder by increasing the truncation order at fixed parameters. Credit is due for the self-consistent derivation, the correct point-mass limit, and the numerical checks across several models; nevertheless, because the central claim is stated broadly as 'establishing a direct connection,' the absence of any error control in the regime where the formulas are showcased makes acceptance conditional on adding this convergence check or explicitly narrowing the claimed validity range.","tokens_in":30529,"tokens_out":9859,"duration_ms":117507,"concrete_test":"For the gNFW case of Fig. 3 (γ=0.5, ρc=4×10^8 M⊙/kpc3, rm=18 kpc, rs=rd=10^7 rm), compute the deflection angle from Eq. (70), truncating at orders N=3, 5, 7, 9 in rm/b, and compare with high-precision numerical integration of the geodesic equation (32) at b/rm=1, 2, 3, 5, 10. If |Δφ_N − Δφ_num| does not decrease systematically with increasing N at b/rm=1–3, the term-by-term integration is not justified there and the central claim should be restricted to b≫rm; if it does decrease, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The master formula (43) is obtained by inserting the small-u expansion (40) of y(u/b) into the deflection integral (37) and integrating term-by-term up to u=1. The expansion point is u=0, i.e., r→∞, but for a ray with impact parameter b the endpoint u=1 corresponds to r=b. For b only a few times rm—exactly the range displayed in Figs. 3, 4 and 6—the asymptotic metric expansions (12) are not uniformly small over the integration domain, so the remainder of the term-by-term series is uncontrolled. The paper itself states that the series are asymptotic and that the method is valid for b≳O(rm), but no explicit error estimate is given, and the numerical comparisons are made at fixed truncation order (typically third order in rm/b) without a convergence check. Thus the central claim that deflection angles and apparent image angles are determined directly by the density coefficients is a validated conjecture for the tested profiles rather than an established procedure in the regime where it is applied.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a perturbative pipeline that starts from a static, spherically symmetric perfect-fluid density profile ρ(r), solves the TOV equations by series methods, and then computes the weak-field deflection angle of null or timelike signals as a quasi-series in the impact parameter b and ln b. Two styles of solution are treated: asymptotic series for densities extending to infinity (integer and rational fractional powers) and Taylor expansions around a finite boundary for truncated models. The resulting master formula, Eq. (43), is applied to the gNFW/NFW, Hernquist, power-law, SIS, PIS, and uniform-density profiles. The deflection formulas are then fed into the lens equation to obtain apparent image angles, and the analytical results are compared with numerical integration of the TOV equations and geodesic equations for a range of model parameters.","tokens_in":30640,"tokens_out":8785,"duration_ms":89028,"significance":"If the method is fully justified, it provides a useful and efficient map from density-profile coefficients to observable lensing quantities, unifying null and timelike signals and including finite source/detector distance effects. The derivations are explicit and the series coefficients are shown to reduce correctly in the point-mass limit. The numerical comparisons in Figs. 2-8 cover several standard astrophysical profiles and show good agreement over the displayed ranges. At the same time, the central claim is stronger than what is rigorously established: the term-by-term integration that produces the master deflection formula is not supplied with an error estimate, and the pressure ans\"atze are assumptions rather than proven consequences of the TOV equations. The paper therefore currently establishes a well-tested numerical method for a class of profiles, with the 'direct connection' between density and lensing holding conditionally on those unproved steps.","major_comments":[{"comment":"The master formula (43) is obtained by expanding y(u/b) about u=0, i.e., about r→∞, and integrating term by term up to u=1, i.e., down to r=b. For b≳rm, which is exactly the range displayed in Figs. 3, 4(c), and 6, the asymptotic metric expansions (12) are not uniformly small over the entire integration interval, so the remainder of the truncated series is uncontrolled. The paper states that the series are asymptotic and that the method is valid for b≳O(rm), but no error bound or convergence check at increasing truncation order is provided. Because Eq. (43) is the bridge from density coefficients to deflection and lensing, this is a load-bearing gap. I request an explicit remainder estimate, or at minimum a systematic truncation-order convergence study for b/rm near unity, and a corresponding softening of the claim that the procedure 'establishes' the density-to-deflection connection in that regime.","section":"III.A, Eq. (43)"},{"comment":"The perturbative solution of the TOV equation relies on the assumption that the pressure admits the double-logarithmic ansatz (8) in the integer-power case and the fractional-power ansatz (16) in the rational-δ case. No existence or uniqueness statement is given for solutions of Eq. (3a) in these classes, nor is a characterization provided of the density coefficients for which the ansatz is consistent. The abstract's phrase 'fairly arbitrary density distributions' is therefore broader than what is demonstrated. Please state these ans\"atze explicitly as assumptions on ρ(r), and verify that each model in Section IV satisfies the required structure; otherwise the central formulas should be presented as conditional on the ansatz rather than as a general result.","section":"II.A, Eqs. (8) and (16)"},{"comment":"The numerical validation is carried out at very low compactness: for the gNFW and Hernquist plots M0/rm ≈ 6×10^{-6}, and the power-law cases are comparable. In this regime the pressure and nonlinear terms in the TOV solution are extremely small, so the comparison does not stress the pressure ansatz or the nonlinear structure of the equations. I ask the authors to test at least one case with substantially larger M0/rm (while remaining in the weak-deflection regime) or to explicitly restrict the validated claim to the low-compactness regime. Without such a test, the demonstrated agreement is evidence for the series method only in a narrow part of the parameter space relevant to the claimed general connection.","section":"IV.A, Figs. 3, 4(c), and 6"}],"minor_comments":[{"comment":"The notation 'ln b 2' should read '(ln b)/2', and 'O(ε)3' should be 'O(ε^3)'; as typeset these expressions are confusing and could be misread as powers of the logarithm.","section":"Eqs. (44) and (46)"},{"comment":"The captions of Figs. 7 and 8 label the perturbative curves as dashed and the numerical curves as solid, while the main text says the analytical results are solid and the numerical results are dashed; these should be made consistent.","section":"Figs. 7 and 8"},{"comment":"The symbol ε is used collectively for M/b, b/rs,d, and logarithmic higher-order terms without a precise definition; please define the ordering explicitly so that the claimed truncation order is unambiguous.","section":"Section III.A, Eq. (44)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the reader's positive assessment is understandable: the derivations are explicit, the formulas are concrete, and the numerical tests are relevant. My concern is not with the correctness of the algebraic manipulations but with the rigor of the asymptotic integration step, which is central to the paper's main claim. This is fixable by adding a remainder estimate or a truncation-order convergence check and by carefully restating the scope of the claim. I do not see a need for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does what it sets out to do: it gives a perturbative pipeline that goes directly from a mass density profile to deflection angles and lensing image angles for static, spherically symmetric perfect fluids. It extends the authors' earlier perturbative deflection framework to TOV-solved spacetimes, handling both logarithmic and fractional-power series, and works out explicit formulas for gNFW, NFW, Hernquist, power-law, SIS, and PIS profiles, for null and timelike signals, with finite-distance corrections included. The numerical comparisons against direct integration (Figs. 2–8) are genuine and match well in the claimed regime b ≳ rm. That is real work and a useful toolbox for lensing theory.\n\nThe main soft spot is the one the stress-test note flags: the deflection master formula (43) is obtained by expanding the integrand around u = 0 and integrating term-by-term to u = 1, but for b only a few times rm the expansion is asymptotic and the remainder is uncontrolled. The paper states the series are asymptotic and valid for b ≳ O(rm), and it checks against numerics at a fixed truncation order (third order in rm/b), but it does not establish an error bound or check convergence by varying the truncation order. So the central claim that density parameters directly determine deflection and image angles is a validated conjecture for the tested profiles, not an established general procedure. If a future density produced metric coefficients outside the assumed ansatz, the connection could break. This is serious but not fatal: the numerics do support the claim in the tested cases, and the paper is honest about the asymptotic nature.\n\nMinor issues: no code or data are provided, which makes independent reproduction slower; the pressure ansatz (8) and (23) is assumed rather than derived from general principles; and the finite-boundary method is only demonstrated for b ≳ R/2. A few equations (e.g., 77a) have typesetting slips that make them hard to parse, but they do not affect the argument.\n\nFor a reader working on gravitational lensing models of dark matter halos, this is a useful reference and a reasonable starting point. It deserves a serious referee. I would ask the authors to add a convergence check (e.g., comparing truncation orders) and to state clearly that the method is validated numerically rather than proven general.","headline":"A useful perturbative pipeline from density to deflection for SSS perfect fluids, with honest numerics, but the central claim is a validated conjecture rather than a proven procedure because the series remainder is uncontrolled.","tokens_in":31222,"tokens_out":2321,"would_cite":false,"duration_ms":22961,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.62.Sb"],"model":"deepseek-v4-flash","headline":"A perturbative method derives gravitational lensing straight from a galaxy's mass-density profile.","keywords":["gravitational lensing","deflection angle","Tolman-Oppenheimer-Volkoff equations","perturbative methods","Navarro-Frenk-White profile","Hernquist profile","perfect fluid","timelike geodesics"],"falsifier":"Take a truncated model with an exact solution, such as the uniform-density sphere or the singular isothermal sphere already treated in the paper, and compare the series deflection from Eq. (103) to the exact numerical deflection at an impact parameter just above $R/2$, increasing the truncation order. If the series does not converge to the numerical value as the order grows, the central claim that density coefficients determine the deflection fails in that regime.","tokens_in":30226,"feed_emoji":"🔭","tokens_out":5404,"duration_ms":51747,"temperature":0.7,"pith_summary":"This paper tries to establish a direct, step-by-step link between the density distribution of a spherical galaxy or cluster and the images it produces as a gravitational lens. Starting from a density profile, the authors solve the Tolman-Oppenheimer-Volkoff equations as series, use those series to write the bending angle of null or timelike rays as a power series in the impact parameter, and then solve the lens equation to get the apparent positions of the two images. The point is that for a wide family of realistic profiles—gNFW, NFW, Hernquist, power-law, pseudo-isothermal, and uniform—the whole chain can be carried out analytically, where previously one would need a numerical spacetime solution. The paper checks the series against numerical integration for several choices of the profile parameters and finds agreement where the impact parameter is at or beyond the halo scale.","feed_headline":"Density profiles predict lensed image angles directly","feed_subtitle":"A perturbative method turns spherical mass-density profiles into deflection angles and image positions for galaxies and clusters.","key_machinery":"The load-bearing object is a double-series ansatz for the pressure, $P(r)=\\sum_{n,m} P_{n,m}(\\ln r)^m/r^n$ for asymptotically flat profiles and $P(r)=\\sum_n P_n(r-R)^n$ for profiles with a finite boundary, together with the matching series for the metric functions $A(r)$ and $B(r)$. Inserting these into the TOV equations fixes the coefficients by undetermined coefficients, and a change of variables reduces the geodesic bending integral to integrals $I_{n,m}$ that obey a recurrence; the master formulas (43), (52), and (61) then deliver the deflection as a series whose terms are the density coefficients.","core_discovery":"In the paper's own terms, the central discovery is that the deflection angle in a static, spherically symmetric perfect-fluid spacetime is a quasi-series in $1/b$ and $\\ln b$ whose coefficients are fixed by the expansion coefficients of the density. A nonzero $\\rho_3$ term in the density's asymptotic expansion produces logarithms in the metric and hence in the deflection; the same perturbative machinery also handles non-integer power-law tails, whose leading deflection behaves as $b^{-\\delta}$, and finite matter distributions, for which the deflection splits into a Schwarzschild part outside the boundary and a series inside. Feeding the deflection into the lens equation gives the impact parameters $b_\\pm$ and the apparent angles $\\beta_\\pm$ of the two images, with finite source and detector distances included.","pith_inferences":["Because the series are asymptotic rather than proven convergent, a natural testable extension is to push the truncation order up for impact parameters just inside the halo scale and see whether the predicted deflection keeps tracking numerical integration or begins to diverge.","The paper's treatment of the logarithmic $\\rho_3$ term reroutes the total mass through a length scale $l$; a reader could infer that strong-lensing observables in gNFW halos carry a weak sensitivity to the mass-assembly history or truncation of the profile, which might be probed with stacked lensing data.","The same coefficient-matching machinery should apply to other observables the paper names as future work, notably time delays in strong lensing, since the metric series is already in hand.","One could test the boundary-expansion method's reliability by using the exact uniform-density and SIS solutions already cited in the paper as benchmarks, comparing the series deflection at $b\\approx R/2$ to the exact result."],"forward_implications":["For any density profile that fits the three-parameter family of Eq. (62), lensing observables—deflection, image positions, their dependence on halo parameters—can be written down analytically to arbitrary order in $1/b$ instead of being extracted from a numerical metric.","The appearance of $\\ln b$ terms tied to $\\rho_3$ means that the logarithmic mass divergence of profiles like gNFW is encoded directly in the deflection, and the paper's formulas show exactly which density coefficients drive the logarithms.","Timelike signals (with velocity $v<1$) are treated on the same footing as light, so the same series gives lensing of relativistic particles as well as photons.","Finite source and detector distances appear naturally in the deflection, which turns the lens equation into a solvable polynomial or quasi-polynomial for $b_\\pm$ and yields the apparent angles $\\beta_\\pm$.","For truncated profiles, the method automatically combines the exterior Schwarzschild bending with the interior matter contribution, so it applies to realistic halos cut off at a virial or truncation radius."],"supporting_citations":[{"why":"supplies the TOV equations and the Schwarzschild exterior that the perturbative solution must match.","marker":"[31]"},{"why":"supplies the perturbative change of variables and integrand expansion that the paper extends to logarithmic and fractional-power metrics.","marker":"[32]"},{"why":"supplies the deflection-angle expansion with finite source and detector distances used in the master formula.","marker":"[33]"},{"why":"supplies the Schwarzschild apparent-angle formula and lensing setup that the paper generalizes to extended density profiles.","marker":"[34]"},{"why":"introduces the NFW density model, the key astrophysical target of the method and a subclass of gNFW.","marker":"[29]"},{"why":"introduces the Hernquist profile, one of the models the paper treats explicitly.","marker":"[30]"}],"fun_headline_variants":["Density profile predicts lensing angles directly","Perturbative bridge: mass density to deflection","One method: density to deflection to images","Weak-field lensing from arbitrary density profiles","A perturbative path from density to Einstein angles"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method stands on the assumption that the pressure and metric admit the postulated series ansatze and that these series can be integrated term-by-term to give a deflection angle accurate for impact parameters at least of order the halo scale; if a realistic density profile produced pressure outside these ansatze, or the series stopped converging before the integration range, the density-to-deflection link would break.","fun_headline_variants_meta":{"raw":{"variants":["Density profile predicts lensing angles directly","Perturbative bridge: mass density to deflection","One method: density to deflection to images","Weak-field lensing from arbitrary density profiles","A perturbative path from density to Einstein angles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000157,"raw_usage":{"total_tokens":1198,"prompt_tokens":900,"completion_tokens":298,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":516,"completion_tokens_details":{"reasoning_tokens":229}},"tokens_in":516,"tokens_out":298,"duration_ms":3698,"temperature":1.0,"reasoning_tokens":229,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T13:36:46.817903+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a truncated model with an exact solution, such as the uniform-density sphere or the singular isothermal sphere already treated in the paper, and compare the series deflection from Eq. (103) to the exact numerical deflection at an impact parameter just above $R/2$, increasing the truncation order. If the series does not converge to the numerical value as the order grows, the central claim that density coefficients determine the deflection fails in that regime.","supporting_citations":[{"cited_title":"Joyce, L","cited_arxiv_id":null,"evidence_quote":"supplies the TOV equations and the Schwarzschild exterior that the perturbative solution must match."},{"cited_title":"Perlick, Living reviews in relativity 7, 1 (2004)","cited_arxiv_id":null,"evidence_quote":"supplies the perturbative change of variables and integrand expansion that the paper extends to logarithmic and fractional-power metrics."},{"cited_title":"Bozza, F","cited_arxiv_id":null,"evidence_quote":"supplies the deflection-angle expansion with finite source and detector distances used in the master formula."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduces the NFW density model, the key astrophysical target of the method and a subclass of gNFW."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduces the Hernquist profile, one of the models the paper treats explicitly."}],"review_version":1}