The cosmic web from perturbation theory

Context. Analysing the large-scale structure (LSS) in the Universe with galaxy surveys demands accurate structure formation models. Such models should ideally be fast and have a clear theoretical framework in order to rapidly scan a variety of cosmological parameter spaces without requiring large training data sets. Aims. This study aims to extend Lagrangian perturbation theory (LPT), including viscosity and vorticity, to reproduce the cosmic evolution from dark matter N-body calculations at the field level. Methods. We extend LPT to a Eulerian framework, which we dub eALPT. An ultraviolet regularisation through the spherical collapse model provided by Augmented LPT turns out to be crucial at low redshifts. This iterative method enables modelling of the stress tensor and introduces vorticity. The eALPT model has two free parameters apart from the choice of cosmology, redshift snapshots, cosmic volume, and the number of particles. Results. We find that compared to N-body solvers, the cross-correlation of the dark matter distribution increases at k = 1 h Mpc(-1) and z = 0 from similar to 55% with the Zel'dovich approximation (similar to 70% with ALPT), to similar to 95% with the three-timestep eALPT, and the power spectra show percentage accuracy up to k similar or equal to E; 0.3 h Mpc(-1).