Validity of the Born approximation for beyond Gaussian weak lensing observables

Accurate forward modeling of weak lensing (WL) observables from cosmological parameters is necessary for upcoming galaxy surveys. Because WL probes structures in the non-linear regime, analytical forward modeling is very challenging, if not impossible. Numerical simulations of WL features rely on ray-tracing through the outputs of N-body simulations, which requires knowledge of the gravitational potential and accurate solvers for light ray trajectories. A less accurate procedure, based on the Born approximation, only requires knowledge of the density field, and can be implemented more efficiently and at a lower computational cost. In this work, we use simulations to show that deviations of the Born-approximated convergence power spectrum, skewness and kurtosis from their fully ray--traced counterparts are consistent with the smallest non-trivial O(Φ3) post-Born corrections (so-called geodesic and lens-lens terms). Our results imply a cancellation among the larger O(Φ4) (and higher order) terms, consistent with previous analytic work. We also find that cosmological parameter bias induced by the Born approximated power spectrum is negligible even for an LSST-like survey, once galaxy shape noise is considered. When considering higher order statistics such as the κ skewness and kurtosis, however, we find significant bias of up to 2.5σ. Using the LensTools software suite, we show that the Born approximation saves a factor of 4 in computing time with respect to the full ray-tracing in reconstructing the convergence.