pith. sign in

arxiv: 2606.19472 · v2 · pith:GNNOXQSCnew · submitted 2026-06-17 · 🌌 astro-ph.CO · gr-qc

Statistical Field Theory for Weak Gravitational Lensing

Pith reviewed 2026-06-26 19:32 UTC · model grok-4.3

classification 🌌 astro-ph.CO gr-qc
keywords weak gravitational lensingSachs optical scalarsstochastic field theoryRicci focusingWeyl shearingdiagrammatic expansioncumulant hierarchyselection rule
0
0 comments X

The pith

Weak lensing is formulated as a stochastic field theory for Sachs optical scalars driven by random Ricci-focusing and Weyl-shearing fields.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper replaces the usual linear remapping picture of weak lensing with a path-integral formulation in which the Sachs optical scalars evolve as a stochastic field driven by independent random Ricci and Weyl fields. The integral produces a systematic diagrammatic series for any n-point lensing statistic, with the conventional Born-approximation result recovered as the leading linear-propagation term. Higher-order terms arise when nonlinear Sachs evolution couples to the non-Gaussian cumulants of the driving fields according to a selection rule that links an n-point observable directly to the n-point driving cumulant and its first correction to the (n+1)-point cumulant. A reader would care because the same framework automatically incorporates path corrections, small-scale mode coupling, and arbitrary matter statistics without having to add them by hand.

Core claim

We formulate lensing as a stochastic field theory for the Sachs optical scalars, driven by random Ricci-focusing and Weyl-shearing fields. The resulting path integral generates a diagrammatic expansion for arbitrary n-point correlation functions of lensing observables, organised into linear response, nonlinear propagation, and driving-field cumulants. The conventional calculation emerges as the lowest-order, linear-propagation limit. Beyond it, nonlinear Sachs evolution couples to driving-field non-Gaussianity, mixing the matter cumulant hierarchy into the lensing hierarchy. A selection rule governs the couplings: an n-point observable receives a direct contribution from the n-point driving-

What carries the argument

The path integral over Sachs optical scalars whose evolution is sourced by independent stochastic Ricci-focusing and Weyl-shearing fields, together with the selection rule that organises how driving-field cumulants enter the lensing hierarchy at successive orders.

If this is right

  • The two-point lensing function receives corrections from squeezed three-point cumulants of Ricci focusing and Weyl shearing, allowing small-scale modes to source larger-scale power.
  • These corrections feed the lensing E- and B-modes equally.
  • Higher-order lensing statistics automatically mix contributions from higher driving-field cumulants without separate modeling.
  • Path corrections and stochasticity appear at the same perturbative order as nonlinear propagation effects.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The selection rule suggests that measurements of lensing bispectra could be used to constrain the three-point statistics of the driving fields rather than treating them as nuisance parameters.
  • Because the formalism keeps the driving fields explicit, it offers a route to include baryonic or other small-scale physics by modifying only the cumulants of Ricci and Weyl fields.
  • The same diagrammatic structure could be applied to other line-of-sight integrals in cosmology, such as the integrated Sachs-Wolfe effect or cosmic magnification, by identifying the analogous optical scalars.

Load-bearing premise

The Sachs optical scalars can be treated as a stochastic field whose statistics are fully determined by the cumulants of independent random Ricci and Weyl driving fields, with those cumulants entering the lensing observables according to the stated selection rule.

What would settle it

A direct numerical extraction, from ray-tracing through N-body simulations, of the squeezed three-point cumulants of the Ricci and Weyl fields along lines of sight and a test of whether they produce the predicted correction to the lensing two-point function at the amplitude and scale dependence given by the selection rule.

Figures

Figures reproduced from arXiv: 2606.19472 by Andrina Nicola, Chris Clarkson, Philip Bull, Zheng Zhang.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic of a null geodesic congruence evolved along [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Second-order diagrams for the Sachs two-point correlation function in the truncation with the local cubic vertex [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. First-order diagrams for the Sachs three-point correlation function generated by either one local cubic vertex [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The selection rule of Eq. (48) at a glance: which driving-field cumulant [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The two scalar branches of the one-to-one light-cone coordinate map. [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Pair-aligned component dictionary and field pattern. The inset fixes the shared screen basis. The left dictionary [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Map-level relation between matter fluctuations and the Sachs driving fields in the scalar sector. Panel A shows a [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. FLRW response propagator [Eq. (78)] for [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Slices of different components of the correlation propagator: [PITH_FULL_IMAGE:figures/full_fig_p016_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Free-theory (Order-0) prediction of our path-integral pipeline against an independent [PITH_FULL_IMAGE:figures/full_fig_p019_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Order-0 and Order-2 decomposition of the four 2PCFs at [PITH_FULL_IMAGE:figures/full_fig_p020_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Angular power spectra of the four observables of Fig. 11: the full prediction (Order-0+FF+FK, grey), Order-0+FF [PITH_FULL_IMAGE:figures/full_fig_p021_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Polarization decomposition of the two Order-2 corrections at [PITH_FULL_IMAGE:figures/full_fig_p023_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. Source-redshift dependence of the convergence two-point statistics. Top: the convergence 2PCF [PITH_FULL_IMAGE:figures/full_fig_p023_14.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16 [PITH_FULL_IMAGE:figures/full_fig_p028_16.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. The equal-shell driving-field three-point cumulant on the squeezed/collapsed configuration (1 [PITH_FULL_IMAGE:figures/full_fig_p029_15.png] view at source ↗
read the original abstract

Standard weak-lensing calculations treat lensing as a linear remapping of the matter field along the line of sight. We instead formulate lensing as a stochastic field theory for the Sachs optical scalars, driven by random Ricci-focusing and Weyl-shearing fields. The resulting path integral generates a diagrammatic expansion for arbitrary $n$-point correlation functions of lensing observables, organised into linear response, nonlinear propagation, and driving-field cumulants. The conventional calculation emerges as the lowest-order, linear-propagation limit. Beyond it, nonlinear Sachs evolution couples to driving-field non-Gaussianity, mixing the matter cumulant hierarchy into the lensing hierarchy. A selection rule governs the couplings: an $n$-point observable receives a direct contribution from the $n$-point driving-field cumulant, and its leading hierarchy-mixing correction from the $(n+1)$-point cumulant via one nonlinear Sachs interaction, with higher cumulants entering only at higher order. The two-point function, for instance, is corrected by squeezed three-point cumulants of Ricci focusing and Weyl shearing, letting small-scale modes source larger scales and feeding the lensing $E$- and $B$-modes equally. Rather than a restrictive approximation scheme, the formalism is a paradigm shift: a unified framework naturally accommodating path corrections, higher-order matter statistics, stochasticity, and small-scale effects.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper formulates weak gravitational lensing as a stochastic field theory for the Sachs optical scalars, driven by random Ricci-focusing and Weyl-shearing fields. The resulting path integral yields a diagrammatic expansion for n-point lensing correlation functions organized by linear response, nonlinear propagation, and driving-field cumulants. The standard linear remapping result is recovered as the lowest-order limit. A selection rule is stated whereby an n-point lensing observable receives a direct drive from the n-point driving cumulant and its leading correction from the (n+1)-point cumulant via a single nonlinear Sachs interaction; this is illustrated for the two-point function by squeezed three-point cumulants of Ricci and Weyl fields that source E- and B-modes equally. The framework is presented as a paradigm shift that unifies path corrections, higher-order matter statistics, stochasticity, and small-scale effects.

Significance. If the selection rule and the underlying stochastic action can be derived explicitly from the Sachs equations, the approach would supply a systematic, non-perturbative organizing principle for incorporating non-Gaussian matter statistics and nonlinear propagation into weak-lensing observables. This could affect the interpretation of small-scale signals in upcoming surveys and the separation of E/B modes without ad-hoc cutoffs.

major comments (2)
  1. [Abstract] Abstract, final paragraph: the selection rule (n-point observable driven directly by n-point cumulant, leading correction from (n+1)-point cumulant via exactly one nonlinear Sachs interaction) is asserted without an explicit derivation from the nonlinear terms in the Sachs optical-scalar evolution or from the stochastic measure. Because this rule is the sole justification for the claim that the formalism 'automatically and exactly accommodates arbitrary higher-order matter statistics,' its status as a theorem rather than an assumption must be demonstrated before the paradigm-shift assertion can be evaluated.
  2. [Abstract] Abstract: the statement that 'the two-point function... is corrected by squeezed three-point cumulants... feeding the lensing E- and B-modes equally' is presented as a direct consequence of the selection rule, yet no explicit vertex or propagator is shown that would enforce equal E/B feeding at this order. Without the interaction Lagrangian or the resulting Feynman rules, it is impossible to confirm that lower-order mixings are forbidden at the claimed perturbative order.
minor comments (1)
  1. The abstract refers to 'the stated selection rule' but supplies neither the explicit Sachs equations nor the form of the stochastic action; a dedicated section deriving the vertices from the optical-scalar evolution would clarify the construction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments. We address each major comment below.

read point-by-point responses
  1. Referee: [Abstract] Abstract, final paragraph: the selection rule (n-point observable driven directly by n-point cumulant, leading correction from (n+1)-point cumulant via exactly one nonlinear Sachs interaction) is asserted without an explicit derivation from the nonlinear terms in the Sachs optical-scalar evolution or from the stochastic measure. Because this rule is the sole justification for the claim that the formalism 'automatically and exactly accommodates arbitrary higher-order matter statistics,' its status as a theorem rather than an assumption must be demonstrated before the paradigm-shift assertion can be evaluated.

    Authors: The selection rule is derived explicitly from the nonlinear Sachs equations and the stochastic path-integral measure in Sections 3–4 of the manuscript. The abstract summarizes this derived result. To make its status as a theorem unambiguous, we will revise the abstract to reference the relevant sections and add a concise statement that the rule follows from the explicit construction. revision: yes

  2. Referee: [Abstract] Abstract: the statement that 'the two-point function... is corrected by squeezed three-point cumulants... feeding the lensing E- and B-modes equally' is presented as a direct consequence of the selection rule, yet no explicit vertex or propagator is shown that would enforce equal E/B feeding at this order. Without the interaction Lagrangian or the resulting Feynman rules, it is impossible to confirm that lower-order mixings are forbidden at the claimed perturbative order.

    Authors: The interaction Lagrangian, vertices, and resulting Feynman rules are derived in Section 5; the equal E/B feeding at this order follows directly from the symmetry properties of the Weyl and Ricci terms under the selection rule, with no lower-order mixings permitted. We will revise the manuscript to include an explicit illustrative diagram and short calculation of this term, making the connection self-contained. revision: yes

Circularity Check

0 steps flagged

Derivation chain self-contained; no circular reductions identified

full rationale

The abstract presents the stochastic field theory and selection rule as generated by the path integral over Sachs optical scalars driven by Ricci/Weyl fields, with the linear result recovered as the explicit lowest-order limit. No equations are shown reducing a claimed prediction to a fitted input, self-citation, or definitional tautology. The selection rule is asserted to follow from the diagrammatic expansion rather than being presupposed as an input; absent any quoted reduction of the form 'n-point observable equals n-point cumulant by construction,' the framework remains non-circular. This is the expected outcome for a paper whose central construction is an independent reformulation rather than a re-labeling of prior results.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review provides no explicit free parameters, axioms, or invented entities beyond the modeling choice itself; the central claim rests on the unverified assumption that the path-integral construction correctly captures the Sachs evolution.

axioms (1)
  • domain assumption Lensing observables are generated by a path integral over stochastic Sachs optical scalars driven by random Ricci and Weyl fields.
    Stated in the abstract as the starting point of the formulation.

pith-pipeline@v0.9.1-grok · 5772 in / 1360 out tokens · 22986 ms · 2026-06-26T19:32:18.925393+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

45 extracted references · 10 linked inside Pith

  1. [1]

    Asgari, C.-A

    M. Asgari, C.-A. Lin, B. Joachimi, B. Giblin, C. Hey- mans, H. Hildebrandt, A. Kannawadi, B. St¨ olzner, T. Tr¨ oster, J. L. van den Busch,et al., Kids-1000 cos- mology: Cosmic shear constraints and comparison be- tween two point statistics, Astronomy & Astrophysics 645, A104 (2021). 29 1 3 10 30 100 300 [arcmin] 397 1823 2261 2314 2327affine shell [Mpc] ...

  2. [2]

    A. Amon, D. Gruen, M. A. Troxel, N. MacCrann, S. Do- delson, A. Choi, C. Doux, L. F. Secco, S. Samuroff, E. Krause,et al., Dark energy survey year 3 results: Cos- mology from cosmic shear and robustness to data calibra- tion, Physical Review D105, 023514 (2022)

  3. [3]

    Ivezi´ c, S

    ˇZ. Ivezi´ c, S. M. Kahn, J. A. Tyson, B. Abel, E. Acosta, R. Allsman, D. Alonso, Y. AlSayyad, S. F. Anderson, J. Andrew,et al., Lsst: from science drivers to reference design and anticipated data products, The Astrophysical Journal873, 111 (2019)

  4. [4]

    Laureijs, J

    R. Laureijs, J. Amiaux, S. Arduini, J.-L. Augueres, J. Brinchmann, R. Cole, M. Cropper, C. Dabin, L. Du- vet, A. Ealet,et al., Euclid definition study report, arXiv preprint arXiv:1110.3193 (2011)

  5. [5]

    Spergel, N

    D. Spergel, N. Gehrels, C. Baltay, D. Bennett, J. Breckin- ridge, M. Donahue, A. Dressler, B. S. Gaudi, T. Greene, O. Guyon,et al., Wide-field infrarred survey telescope- astrophysics focused telescope assets wfirst-afta 2015 re- port, ArXiv e-prints , arXiv (2015)

  6. [6]

    Bartelmann and P

    M. Bartelmann and P. Schneider, Weak gravitational lensing, Physics Reports340, 291 (2001)

  7. [7]

    Kilbinger, Cosmology with cosmic shear observations: a review, Reports on Progress in Physics78, 086901 (2015)

    M. Kilbinger, Cosmology with cosmic shear observations: a review, Reports on Progress in Physics78, 086901 (2015)

  8. [8]

    Mandelbaum, Weak lensing for precision cosmology, Annual Review of Astronomy and Astrophysics56, 393 (2018)

    R. Mandelbaum, Weak lensing for precision cosmology, Annual Review of Astronomy and Astrophysics56, 393 (2018)

  9. [9]

    Cooray and W

    A. Cooray and W. Hu, Second-order corrections to weak lensing by large-scale structure, The Astrophysical Jour- nal574, 19 (2002)

  10. [10]

    Dodelson, C

    S. Dodelson, C. Shapiro, and M. White, Reduced shear power spectrum, Physical Review D—Particles, Fields, Gravitation, and Cosmology73, 023009 (2006)

  11. [11]

    Bernardeau, C

    F. Bernardeau, C. Bonvin, and F. Vernizzi, Full-sky lens- ing shear at second order, Physical Review D—Particles, Fields, Gravitation, and Cosmology81, 083002 (2010)

  12. [12]

    Krause and C

    E. Krause and C. M. Hirata, Weak lensing power spec- tra for precision cosmology-multiple-deflection, reduced shear, and lensing bias corrections, Astronomy & Astro- physics523, A28 (2010)

  13. [13]

    Pratten and A

    G. Pratten and A. Lewis, Impact of post-born lensing on the CMB, Journal of Cosmology and Astroparticle Physics2016(08), 047

  14. [14]

    Fabbian, M

    G. Fabbian, M. Calabrese, and C. Carbone, Cmb weak- lensing beyond the born approximation: a numerical ap- proach, Journal of Cosmology and Astroparticle Physics 2018(02), 050

  15. [15]

    Marozzi, G

    G. Marozzi, G. Fanizza, E. Di Dio, and R. Durrer, Cmb lensing beyond the leading order: Temperature and po- larization anisotropies, Physical Review D98, 023535 (2018)

  16. [16]

    Fry, Gravity, bias, and the galaxy three-point correla- tion function, Physical Review Letters73, 215 (1994)

    J. Fry, Gravity, bias, and the galaxy three-point correla- tion function, Physical Review Letters73, 215 (1994)

  17. [17]

    Scoccimarro, H

    R. Scoccimarro, H. A. Feldman, J. N. Fry, and J. A. Frieman, The bispectrum of iras redshift catalogs, The Astrophysical Journal546, 652 (2001)

  18. [18]

    Jeong and E

    D. Jeong and E. Komatsu, Primordial non-gaussianity, scale-dependent bias, and the bispectrum of galaxies, The Astrophysical Journal703, 1230 (2009)

  19. [19]

    Halder, O

    A. Halder, O. Friedrich, S. Seitz, and T. N. Varga, The in- tegrated three-point correlation function of cosmic shear, Monthly Notices of the Royal Astronomical Society506, 2780 (2021)

  20. [20]

    Lemos, A

    P. Lemos, A. Challinor, and G. Efstathiou, The effect of Limber and flat-sky approximations on galaxy weak lens- ing, JCAP2017, 014 (2017), arXiv:1704.01054 [astro- ph.CO]

  21. [21]

    X. Fang, E. Krause, T. Eifler, and N. MacCrann, Beyond Limber: efficient computation of angular power spectra for galaxy clustering and weak lensing, JCAP2020, 010 (2020), arXiv:1911.11947 [astro-ph.CO]

  22. [22]

    T. D. Kitching, J. Alsing, A. F. Heavens, R. Jimenez, 30 J. D. McEwen, and L. Verde, The limits of cosmic shear, MNRAS469, 2737 (2017), arXiv:1611.04954 [astro- ph.CO]

  23. [23]

    Cuesta-Lazaro, A

    C. Cuesta-Lazaro, A. Quera-Bofarull, R. Reischke, and B. M. Sch¨ afer, Gravitational corrections to light propa- gation in a perturbed FLRW universe and correspond- ing weak-lensing spectra, MNRAS477, 741 (2018), arXiv:1801.03325 [astro-ph.CO]

  24. [24]

    Hilbert, S., Hartlap, J., White, S. D. M., and Schneider, P., Ray-tracing through the millennium simulation: Born corrections and lens-lens coupling in cosmic shear and galaxy-galaxy lensing, A&A499, 31 (2009)

  25. [25]

    Petri, Z

    A. Petri, Z. Haiman, L. Hui, M. May, and J. M. Kra- tochvil, Cosmology with Minkowski functionals and mo- ments of the weak lensing convergence field, Phys. Rev. D88, 123002 (2013), arXiv:1309.4460 [astro-ph.CO]

  26. [26]

    Takahashi, T

    R. Takahashi, T. Hamana, M. Shirasaki, T. Namikawa, T. Nishimichi, K. Osato, and K. Shiroyama, Full-sky Gravitational Lensing Simulation for Large-area Galaxy Surveys and Cosmic Microwave Background Experi- ments, Astrophys. J.850, 24 (2017), arXiv:1706.01472 [astro-ph.CO]

  27. [27]

    Sachs, Gravitational waves in general relativity

    R. Sachs, Gravitational waves in general relativity. vi. the outgoing radiation condition, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences264, 309 (1961)

  28. [28]

    P. C. Martin, E. D. Siggia, and H. A. Rose, Statisti- cal Dynamics of Classical Systems, Phys. Rev. A8, 423 (1973)

  29. [29]

    H.-K. Janssen, On a Lagrangean for classical field dynam- ics and renormalization group calculations of dynamical critical properties, Zeitschrift fur Physik B Condensed Matter23, 377 (1976)

  30. [30]

    J. A. Krommes, Fundamental statistical descriptions of plasma turbulence in magnetic fields, Phys. Rep.360, 1 (2002)

  31. [31]

    Zhang, sft-wick: A formalism and package for feynman-diagram expansion and evaluation in stochas- tic field theories (2026), arXiv:2606.19480 [physics.comp- ph]

    Z. Zhang, sft-wick: A formalism and package for feynman-diagram expansion and evaluation in stochas- tic field theories (2026), arXiv:2606.19480 [physics.comp- ph]

  32. [32]

    Bernardeau, S

    F. Bernardeau, S. Colombi, E. Gazta˜ naga, and R. Scoc- cimarro, Large-scale structure of the universe and cos- mological perturbation theory, Physics reports367, 1 (2002)

  33. [33]

    M. P. van Daalen, J. Schaye, C. M. Booth, and C. Dalla Vecchia, The effects of galaxy formation on the mat- ter power spectrum: a challenge for precision cosmol- ogy, MNRAS415, 3649 (2011), arXiv:1104.1174 [astro- ph.CO]

  34. [34]

    J. M. Diego, T. Broadhurst, and G. F. Smoot, Evidence for lensing of gravitational waves from LIGO-Virgo data, Phys. Rev. D104, 103529 (2021), arXiv:2106.06545 [gr- qc]

  35. [35]

    G. F. Ellis, R. Maartens, and M. A. MacCallum,Rela- tivistic cosmology(Cambridge University Press, 2012)

  36. [36]

    Maartens, Covariant cosmology: galaxy worldlines and lightrays (2026), working draft

    R. Maartens, Covariant cosmology: galaxy worldlines and lightrays (2026), working draft

  37. [37]

    Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Bacci- galupi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, N. Bartolo, E. Battaner, R. Battye, K. Benabed, A. Benoˆ ıt, A. Benoit-L´ evy, J.-P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger...

  38. [38]

    Zhang, P

    Z. Zhang, P. Bull, and A. Nicola,CANOES: Cosmologi- cal ANgular Observables and Estimator Suite, (compan- ion paper, in preparation) (2026)

  39. [39]

    Clarkson, G

    C. Clarkson, G. F. Ellis, A. Faltenbacher, R. Maartens, O. Umeh, and J.-P. Uzan, (mis) interpreting supernovae observations in a lumpy universe, Monthly Notices of the Royal Astronomical Society426, 1121 (2012)

  40. [40]

    N. E. Chisari, D. Alonso, E. Krause, C. D. Leonard, P. Bull, J. Neveu, A. Villarreal, S. Singh, T. McClintock, J. Ellison,et al., Core cosmology library: Precision cos- mological predictions for lsst, The Astrophysical Journal Supplement Series242, 2 (2019)

  41. [41]

    Fleury, J

    P. Fleury, J. Larena, and J.-P. Uzan, The theory of stochastic cosmological lensing, Journal of Cosmology and Astroparticle Physics2015(11), 022

  42. [42]

    Fleury, J

    P. Fleury, J. Larena, and J.-P. Uzan, Weak gravitational lensing of finite beams, Physical Review Letters119, 191101 (2017)

  43. [43]

    J. M. Martin-Garcia, A. Garc´ ıa-Parrado, A. Stecchina, 31 B. Wardell, C. Pitrou, D. Brizuela,et al., xact: Efficient tensor computer algebra for the wolfram language (2020)

  44. [44]

    J. M. Mart´ ın-Garc´ ıa, xperm: fast index canonicalization for tensor computer algebra, Computer physics commu- nications179, 597 (2008)

  45. [45]

    D. W. Hogg, Distance measures in cosmology, arXiv preprint astro-ph/9905116 (1999)