We present the Extreme-Horizon (EH) cosmological simulation: EH models galaxy formation with stellar and AGN feedback and uses a very high resolution in the intergalactic and circumgalactic medium. The high resolution in low-density regions results in smaller-size massive galaxies at redshift $z=2$, in better agreement with observations compared to other simulations. This results from the improved modeling of cold gas flows accreting onto galaxies. Besides, the EH simulation forms a population of particularly compact galaxies with stellar masses of $10^{10-11}$\,M$_\sun$ that are reminiscent of observed ultracompact galaxies at $z\simeq2$. These objects form mainly through repeated major mergers of low-mass progenitors, independently of baryonic feedback mechanisms. This formation process can be missed in simulations using a too low resolution in low-density intergalactic regions.

The merging rate of cosmic structures is computed, relying on the Ansatz that they can be predicted in the initial linear density field from the coalescence of critical points with increasing smoothing scale, used here as a proxy for cosmic time. Beyond the mergers of peaks with saddle points (a proxy for halo mergers), we consider the coalescence and nucleation of all sets of critical points, including wall-saddle to filament-saddle and wall-saddle to minima (a proxy for filament and void mergers respectively), as they impact the geometry of galactic infall, and in particular filament disconnection. Analytical predictions of the one-point statistics are validated against multiscale measurements in 2D and 3D realisations of Gaussian random fields (the corresponding code being available upon request) and compared qualitatively to cosmological $N$-body simulations at early times ($z\geq 10$) and large scales ($\geq 5\, \mathrm{Mpc}/h$). The rate of filament coalescence is compared to the merger rate of haloes and the two-point clustering of these events is computed, along with their cross-correlations with critical points. These correlations are qualitatively consistent with the preservation of the connectivity of dark matter haloes, and the impact of the large scale structures on assembly bias. The destruction rate of haloes and voids as a function of mass and redshift is quantified down to $z=0$ for a $\Lambda$CDM cosmology. The one-point statistics in higher dimensions are also presented, together with consistency relations between critical point and critical event counts.

Weak lensing of galaxies and CMB photons through the large-scale structure of the Universe is one of the most promising cosmological probes with upcoming experiments dedicated to its measurements such as Euclid/LSST and CMB Stage 4 experiments. With increasingly precise measurements, there is a dire need for accurate theoretical predictions. In this work, we focus on higher order statistics of the weak lensing convergence field, namely its cumulants such as skewness and kurtosis and its one-point probability distribution (PDF), and we quantify using perturbation theory the corrections coming from post-Born effects, meaning beyond the straight-line and independent lenses approximations. At first order, two such corrections arise: lens-lens couplings and geodesic deviation. Though the corrections are small for low source redshifts (below a few percents) and therefore for galaxy lensing, they become important at higher redshifts, notably in the context of CMB lensing, where the non-gaussianities computed from tree-order perturbation theory are found to be of the same order as the signal itself. We include these post-Born corrections on the skewness into a prediction for the one-point convergence PDF obtained with large deviation theory and successfully test these results against numerical simulations. The modelled PDF is indeed shown to perform better than the percent for apertures above ~ 10 arcminutes and typically in the three sigmas region around the mean.

The distribution of the cosmic convergence field is modeled using a large-deviation principle where all non-Gaussian contributions are computed from first principles. The geometry of the past light cone is accounted for by constructing the total weak-lensing signal from contributions of the matter density in thin disk slices. The prediction of this model is successfully tested against numerical simulation with ray tracing, and found to be accurate within at least 5 per cent in the tails at redshift 1 and opening angle of 10 arcmin and even more so with increasing source redshift and opening angle. An accurate analytical approximation to the theory is also provided for practical implementation. The lensing kernel that mixes physical scales along the line-of-sight tends to reduce the domain of validity of this theoretical approach compared to the three dimensional case of cosmic densities in spherical cells. This effect is shown to be avoidable if a nulling procedure is implemented in order to localise the lensing line-of-sight integrations in a tomographic analysis. Accuracy in the tails is thus achieved within a percent for source redshifts between 0.5 and 1.5 and an opening angle of 10 arcmin. Applications to future weak-lensing surveys like Euclid and the specific issue of shape noise are discussed.

The intrinsic correlations of galaxy shapes and orientations across the large-scale structure of the Universe are a known contaminant to weak gravitational lensing. They are known to be dependent on galaxy properties, such as their mass and morphologies. The complex interplay between alignments and the physical processes that drive galaxy evolution remains vastly unexplored. We assess the sensitivity of intrinsic alignments (shapes and angular momenta) to Active Galactic Nuclei -AGN- feedback by comparing galaxy alignment in twin runs of the cosmological hydrodynamical Horizon simulation, which do and do not include AGN feedback respectively. We measure intrinsic alignments in three dimensions and in projection at z=0 and z=1. We find that the projected alignment signal of all galaxies with resolved shapes with respect to the density field in the simulation is robust to AGN feedback, thus giving similar predictions for contamination to weak lensing. The relative alignment of galaxy shapes around galaxy positions is however significantly impacted, especially when considering high-mass ellipsoids. Using a sample of galaxy "twins" across simulations, we determine that AGN changes both the galaxy selection and their actual alignments. Finally, we measure the alignments of angular momenta of galaxies with their nearest filament. Overall, these are more significant in the presence of AGN as a result of the higher abundance of massive pressure-supported galaxies.

We comprehensively analyse the cosmology dependence of counts-in-cell statistics. We focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly nonlinear scales. Based on large-deviation statistics, we parametrise the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the nonlinear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass $M_\nu$ and can disentangle it from the clustering amplitude $\sigma_8$. Using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the $\nu\Lambda$CDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first $\nu\Lambda$CDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly nonlinear scales. The matter PDF is particularly powerful for constraining the matter density $\Omega_m$, clustering amplitude $\sigma_8$ and the total neutrino mass $M_\nu$. Adding the mildly nonlinear matter PDF to the mildly nonlinear matter power spectrum improves constraints on $\Omega_m$ by a factor of 5 and $\sigma_8$ by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc/h)$^3$. We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.

Elliptical galaxies today appear aligned with the large-scale structure of the Universe, but it is still an open question when they acquire this alignment. Observational data is currently insufficient to provide constraints on the time evolution of intrinsic alignments, and hence existing models range from assuming that galaxies gain some primordial alignment at formation, to suggesting that they react instantaneously to tidal interactions with the large-scale structure. Using the cosmological hydrodynamical simulation Horizon-AGN, we measure the relative alignments between the major axes of galaxies and eigenvectors of the tidal field as a function of redshift. We focus on constraining the time evolution of the alignment of the main progenitors of massive $z=0$ elliptical galaxies, the main weak lensing contaminant at low redshift. We show that this population, which at $z=0$ has a stellar mass above $10^{10.4}$ M$_\odot$, transitions from having no alignment with the tidal field at $z=3$, to a significant alignment by $z=1$. From $z=0.5$ they preserve their alignment at an approximately constant level until $z=0$. We find a mass-dependence of the alignment signal of elliptical progenitors, whereby ellipticals that are less massive today ($10^{10.4}<M/{\rm M}_\odot<10^{10.7}$) do not become aligned till later redshifts ($z<2$), compared to more massive counterparts. We also present an extended study of progenitor alignments in the parameter space of stellar mass and galaxy dynamics, the impact of shape definition and tidal field smoothing.

In the peaks approach, the formation sites of observable structures in the Universe are identified as peaks in the matter density field. The statistical properties of the clustering of peaks are particularly important in this respect. In this paper, we investigate the large-separation expansion of the correlation function of peaks in Gaussian random fields. The analytic formula up to third order is derived, and the resultant expression can be evaluated by a combination of one-dimensional fast Fourier transforms, which are evaluated very fast. The analytic formula obtained perturbatively in the large-separation limit is compared with a method of Monte-Carlo integrations, and a complementarity between the two methods is demonstrated.

We investigate the impact of the number of filaments connected to the nodes of the cosmic web on the physical properties of their galaxies using the Sloan Digital Sky Survey. We compare these measurements to the cosmological hydrodynamical simulations Horizon-(no)AGN and Simba. We find that more massive galaxies are more connected, in qualitative agreement with theoretical predictions and measurements in dark matter only simulation. The star formation activity and morphology of observed galaxies both display some dependence on the connectivity of the cosmic web at fixed stellar mass: less star forming and less rotation supported galaxies also tend to have higher connectivity. These results qualitatively hold both for observed and virtual galaxies, and can be understood given that the cosmic web is the main source of fuel for galaxy growth. The simulations show the same trends at fixed halo mass, suggesting that the geometry of filamentary infall impacts galaxy properties beyond the depth of the local potential well. Based on simulations, it is also found that AGN feedback is key in reversing the relationship between stellar mass and connectivity at fixed halo mass. Technically, connectivity is a practical observational proxy for past and present accretion (minor mergers or diffuse infall).

This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others.

Context. Accurate model predictions including the physics of baryons are required to make the most of the upcoming large cosmological surveys devoted to gravitational lensing. The advent of hydrodynamical cosmological simulations enables such predictions on sufficiently sizeable volumes. Aims. Lensing quantities (deflection, shear, convergence) and their statistics (convergence power spectrum, shear correlation functions, galaxy-galaxy lensing) are computed in the past lightcone built in the Horizon-AGN hydrodynamical cosmological simulation, which implements our best knowledge on baryonic physics at the galaxy scale in order to mimic galaxy populations over cosmic time. Methods. Lensing quantities are generated over a one square degree field of view by performing multiple-lens plane ray-tracing through the lightcone, taking full advantage of the 1 kpc resolution and splitting the line of sight over 500 planes all the way to redshift z~7. Two methods are explored (standard projection of particles with adaptive smoothing, and integration of the acceleration field) to assert a good implementation. The focus is on small scales where baryons matter most. Results. Standard cosmic shear statistics are impacted at the 10% level by the baryonic component for angular scales below a few arcmin. The galaxy-galaxy lensing signal, or galaxy-shear correlation function, is consistent with measurements for the redshift z~0.5 massive galaxy population. At higher redshift z>1, the impact of magnification bias on this correlation is relevant for separations greater than 1 Mpc. Conclusions. This work is pivotal for all current and upcoming weak lensing surveys and represents a first step towards building a full end-to-end generation of lensed mock images from large cosmological hydrodynamical simulations.

We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance.

Cosmic connectivity and multiplicity, i.e. the number of filaments globally or locally connected to a given cluster is a natural probe of the growth of structure and in particular of the nature of dark energy. It is also a critical ingredient driving the assembly history of galaxies as it controls mass and angular momentum accretion. The connectivity of the cosmic web is investigated here via the persistent skeleton. This tool identifies topologically the set of ridges of the cosmic landscape which allows us to investigate in details how the nodes of the cosmic web are connected together. When applied to Gaussian random fields, it is found that on average the nodes are connected to exactly $\kappa=4$ neighbours in two dimensions and ~6.1 in three dimensions. Investigating spatial dimensions up to d=11, typical departures from a cubic lattice $\kappa=2d$ are shown to scale like the power 7/4 of the dimension. These numbers strongly depend on the height of the peaks: the higher the peak the larger the connectivity. Locally, the number of filaments sticking out from a node is not equal to its number of connections. At small distance, the peak has an ellipsoidal shape so that only two filaments emerge. Further away, filaments are subdivided at bifurcation points. At the resolution of the maps, the local number of filaments is 3 on average in 2D and 4 in 3D, depending on the height of the peak. Predictions from first principles based on peak theory are shown to reproduce well the connectivity and multiplicity of Gaussian random fields and cosmological simulations. As an illustration, this connectivity is quantified in galaxy lensing convergence maps and large dark haloes catalogues. As a function of redshift and scale the mean connectivity decreases in a cosmology dependent way. As a function of halo mass it scales like 10/3 times the log of the mass.

The strikingly anisotropic large-scale distribution of matter made of an extended network of voids delimited by sheets, themselves segmented by filaments, within which matter flows towards compact nodes where they intersect, imprints its geometry on the dynamics of cosmic flows, ultimately shaping the distribution of galaxies and the redshift evolution of their properties. The (filament-type) saddle points of this cosmic web provide a local frame in which to quantify the induced physical and morphological evolution of galaxies on large scales. The properties of virtual galaxies within the Horizon-AGN simulation are stacked in such a frame. The iso-contours of the galactic number density, mass, specific star formation rate (sSFR), kinematics and age are clearly aligned with the filament axis with steep gradients perpendicular to the filaments. A comparison to a simulation without feedback from active galactic nuclei (AGN) illustrates its impact on quenching star formation of centrals away from the saddles. The redshift evolution of the properties of galaxies and their age distribution are consistent with the geometry of the bulk flow within that frame. They compare well with expectations from constrained Gaussian random fields and the scaling with the mass of non-linearity, modulo the redshift dependent impact of feedback processes. Physical properties such as sSFR and kinematics seem not to depend only on mean halo mass and density: the residuals trace the geometry of the saddle, which could point to other environment-sensitive physical processes, such as spin advection, and AGN feedback at high mass.

We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method.

This work investigates the alignment of galactic spins with the cosmic web across cosmic time using the cosmological hydrodynamical simulation Horizon-AGN. The cosmic web structure is extracted via the persistent skeleton as implemented in the DISPERSE algorithm. It is found that the spin of low-mass galaxies is more likely to be aligned with the filaments of the cosmic web and to lie within the plane of the walls while more massive galaxies tend to have a spin perpendicular to the axis of the filaments and to the walls. The mass transition is detected with a significance of 9 sigmas. This galactic alignment is consistent with the alignment of the spin of dark haloes found in pure dark matter simulations and with predictions from (anisotropic) tidal torque theory. However, unlike haloes, the alignment of low-mass galaxies is weak and disappears at low redshifts while the orthogonal spin orientation of massive galaxies is strong and increases with time, probably as a result of mergers. At fixed mass, alignments are correlated with galaxy morphology: the high-redshift alignment is dominated by spiral galaxies while elliptical centrals are mainly responsible for the perpendicular signal. These predictions for spin alignments with respect to cosmic filaments and unprecendently walls are successfully compared with existing observations. The alignment of the shape of galaxies with the different components of the cosmic web is also investigated. A coherent and stronger signal is found in terms of shape at high mass. The two regimes probed in this work induce competing galactic alignment signals for weak lensing, with opposite redshift and luminosity evolution. Understanding the details of these intrinsic alignments will be key to exploit future major cosmic shear surveys like Euclid or LSST.

Intensity mapping surveys will provide access to a coarse view of the cosmic large-scale structure in unprecedented large volumes at high redshifts. Given the large fractions of the sky that can be efficiently scanned using emission from cosmic neutral hydrogen (HI), intensity mapping is ideally suited to probe a wide range of density environments and hence to constrain cosmology and fundamental physics. To efficiently extract information from 21cm intensities beyond average, one needs non-Gaussian statistics that capture large deviations from mean HI density. Counts-in-cells statistics are ideally suited for this purpose, as the statistics of matter densities in spheres can be predicted accurately on scales where their variance is below unity. We use a large state-of-the-art magneto-hydrodynamic simulation from the IllustrisTNG project to determine the relation between neutral hydrogen and matter densities in cells. We demonstrate how our theoretical knowledge about the matter PDF for a given cosmology can be used to extract a parametrisation-independent HI bias function from a measured neutral hydrogen PDF. When combining the predicted matter PDFs with a simple bias fit to the simulation, we obtain a prediction for neutral hydrogen PDFs at a few percent accuracy at scale R=5 Mpc/h from redshift z=5 to z=1. Furthermore, we find a density-dependent HI clustering signal that is consistent with theoretical expectations and could allow for joint constraints of HI bias and the amplitude of matter fluctuations or the growth of structure.

The anisotropic distribution of satellites around the central galaxy of their host halo is well-documented. However the relative impact of baryons and dark matter in shaping this distribution is still debated. Using the simulation Horizon-AGN, the angular distribution of satellite galaxies with respect to their central counterpart and halo is quantified. Below one Rvir, satellites cluster more strongly in the plane of the central, rather than merely tracing the shape of their host halo. This is due to the increased isotropy of inner haloes acquired through their inside-out assembly in vorticity-rich flows along the cosmic web. While the effect of centrals decreases with distance, halos' triaxiality increases, impacting more and more the satellite's distribution. Effects become comparable just outside one virial radius. Above this scale, the filamentary infall also impacts the satellites distribution, dominating above two virial radii. The central's morphology plays a governing role: the alignment w.r.t. the central plane is four times stronger in haloes hosting stellar discs than in spheroids. But the impact of the galactic plane decreases for lower satellite-to-central mass ratios, suggesting this might not hold for dwarf satellites of the Local group. The orientation of the Milky-Way's satellites traces their cosmic filament, their level of coplanarity is consistent with systems of similar mass and cosmic location in Horizon-AGN. However, the strong impact of galactic planes in massive groups and clusters bounds the likelihood of finding a relaxed region where satellites can be used to infer halo shape. The minor-to-major axis ratios for haloes with log(M0/Msun)>13.5 is underestimated by 10%. This error soars quickly to 30-40% for individual halo measurements.

Large deviation statistics is implemented to predict the statistics of cosmic densities in cylinders applicable to photometric surveys. It yields few percent accurate analytical predictions for the one-point probability distribution function (PDF) of densities in concentric or compensated cylinders; and also captures the density-dependence of their angular clustering (cylinder bias). All predictions are found to be in excellent agreement with the cosmological simulation Horizon Run 4 in the quasi-linear regime where standard perturbation theory normally breaks down. These results are combined with a simple local bias model that relates dark matter and tracer densities in cylinders and validated on simulated halo catalogues. This formalism can be used to probe cosmology with existing and upcoming photometric surveys like DES, Euclid or WFIRST containing billions of galaxies.

The role of the cosmic web in shaping galaxy properties is investigated in the GAMA spectroscopic survey in the redshift range $0.03 \leq z \leq 0.25$. The stellar mass, $u - r$ dust corrected colour and specific star formation rate (sSFR) of galaxies are analysed as a function of their distances to the 3D cosmic web features, such as nodes, filaments and walls, as reconstructed by DisPerSE. Significant mass and type/colour gradients are found for the whole population, with more massive and/or passive galaxies being located closer to the filament and wall than their less massive and/or star-forming counterparts. Mass segregation persists among the star-forming population alone. The red fraction of galaxies increases when closing in on nodes, and on filaments regardless of the distance to nodes. Similarly, the star-forming population reddens (or lowers its sSFR) at fixed mass when closing in on filament, implying that some quenching takes place. Comparable trends are also found in the state-of-the-art hydrodynamical simulation Horizon-AGN. These results suggest that on top of stellar mass and large-scale density, the traceless component of the tides from the anisotropic large-scale environment also shapes galactic properties. An extension of excursion theory accounting for filamentary tides provides a qualitative explanation in terms of anisotropic assembly bias: at a given mass, the accretion rate varies with the orientation and distance to filaments. It also explains the absence of type/colour gradients in the data on smaller, non-linear scales.

The mass, accretion rate and formation time of dark matter haloes near proto-filaments (identified as saddle points of the potential) are analytically predicted using a conditional version of the excursion set approach in its so-called "upcrossing" approximation. The model predicts that at fixed mass, mass accretion rate and formation time vary with orientation and distance from the saddle, demonstrating that assembly bias is indeed influenced by the tides imposed by the cosmic web. Starved, early forming haloes of smaller mass lie preferentially along the main axis of filaments, while more massive and younger haloes are found closer to the nodes. Distinct gradients for distinct tracers such as typical mass and accretion rate occur because the saddle condition is anisotropic, and because the statistics of these observables depend on both the conditional means and their covariances. The theory is extended to other critical points of the potential field. The response of the mass function to variations of the matter density field (the so-called large scale bias) is computed, and its trend with accretion rate is shown to invert along the filament. The signature of this model should correspond at low redshift to an excess of reddened galactic hosts at fixed mass along preferred directions, as recently reported in spectroscopic and photometric surveys and in hydrodynamical simulations. The anisotropy of the cosmic web emerges therefore as a significant ingredient to describe jointly the dynamics and physics of galaxies, e.g. in the context of intrinsic alignments or morphological diversity.

Intrinsic alignments of galaxies are a significant astrophysical systematic affecting cosmological constraints from weak gravitational lensing. Obtaining numerical predictions from hydrodynamical simulations of expected survey volumes is expensive, and a cheaper alternative relies on populating large dark matter-only simulations with accurate models of alignments calibrated on smaller hydrodynamical runs. This requires connecting the shapes and orientations of galaxies to those of dark matter haloes and to the large-scale structure. In this paper, we characterise galaxy-halo alignments in the Horizon-AGN cosmological hydrodynamical simulation. We compare the shapes and orientations of galaxies in the redshift range $0<z<3$ to those of their embedding dark matter haloes, and to the matching haloes of a twin dark-matter only run with identical initial conditions. We find that galaxy ellipticities in general cannot be predicted directly from halo ellipticities. The mean misalignment angle between the minor axis of a galaxy and its embedding halo is a function of halo mass, with residuals arising from the dependence of alignment on galaxy type, but not on environment. Haloes are much more strongly aligned among themselves than galaxies, and they decrease their alignment towards low redshift. Galaxy alignments compete with this effect, as galaxies tend to increase their alignment with haloes towards low redshift. We discuss the implications of these results for current halo models of intrinsic alignments and suggest several avenues for improvement.

Non-Gaussianities of dynamical origin are disentangled from primordial ones using the formalism of large deviation statistics with spherical collapse dynamics. This is achieved by relying on accurate analytical predictions for the one-point probability distribution function (PDF) and the two-point clustering of spherically-averaged cosmic densities (sphere bias). Sphere bias extends the idea of halo bias to intermediate density environments and voids as underdense regions. In the presence of primordial non-Gaussianity, sphere bias displays a strong scale dependence relevant for both high and low density regions, which is predicted analytically. The statistics of densities in spheres are built to model primordial non-Gaussianity via an initial skewness with a scale-dependence that depends on the bispectrum of the underlying model. The analytical formulas with the measured nonlinear dark matter variance as input are successfully tested against numerical simulations. For local non-Gaussianity with a range from $f_{\rm NL}=-100$ to $+100$ they are found to agree within 2\% or better for densities $\rho\in[0.5,3]$ in spheres of radius 15 Mpc$/h$ down to $z=0.35$. The validity of the large deviation statistics formalism is thereby established for all observationally relevant local-type departures from perfectly Gaussian initial conditions. The corresponding estimators for the amplitude of the nonlinear variance $\sigma_8$ and primordial skewness $f_{\rm NL}$ are validated using a fiducial joint maximum likelihood experiment. The influence of observational effects and the prospects for a future detection of primordial non-Gaussianity from joint one- and two-point densities-in-spheres statistics are discussed.

Gravitational lensing allows to quantify the angular distribution of the convergence field around clusters of galaxies to constrain their connectivity to the cosmic web. We describe in this paper the corresponding theory in Lagrangian space where analytical results can be obtained by identifying clusters to peaks in the initial field. We derive the three-point Gaussian statistics of a two-dimensional field and its first and second derivatives. The formalism allows us to study the statistics of the field in a shell around a central peak, in particular its multipolar decomposition. The peak condition is shown to significantly remove power from the dipolar contribution and to modify the monopole and quadrupole. As expected, higher order multipoles are not significantly modified by the constraint. Analytical predictions are successfully checked against measurements in Gaussian random fields. The effect of substructures and radial weighting is shown to be small and does not change the qualitative picture. The non-linear evolution is shown to induce a non-linear bias of all multipoles proportional to the cluster mass.We predict the Gaussian and weakly non-Gaussian statistics of multipolar moments of a two-dimensional field around a peak as a proxy for the azimuthal distribution of the convergence field around a cluster of galaxies. A quantitative estimate of this multipolar decomposition of the convergence field around clusters in numerical simulations of structure formation and in observations will be presented in two forthcoming papers.

Context. Upcoming weak lensing surveys such as Euclid will provide an unprecedented opportunity to quantify the geometry and topology of the cosmic web, in particular in the vicinity of lensing clusters. Aims. Understanding the connectivity of the cosmic web with unbiased mass tracers, such as weak lensing, is of prime importance to probe the underlying cosmology, seek dynamical signatures of dark matter, and quantify environmental effects on galaxy formation. Methods. Mock catalogues of galaxy clusters are extracted from the N-body PLUS simulation. For each cluster, the aperture multipolar moments of the convergence are calculated in two annuli (inside and outside the virial radius). By stacking their modulus, a statistical estimator is built to characterise the angular mass distribution around clusters. The moments are compared to predictions from perturbation theory and spherical collapse. Results. The main weakly chromatic excess of multipolar power on large scales is understood as arising from the contraction of the primordial cosmic web driven by the growing potential well of the cluster. Besides this boost, the quadrupole prevails in the cluster (ellipsoidal) core, while at the outskirts, harmonic distortions are spread on small angular modes, and trace the non-linear sharpening of the filamentary structures. Predictions for the signal amplitude as a function of the cluster-centric distance, mass, and redshift are presented. The prospects of measuring this signal are estimated for current and future lensing data sets. Conclusions. The Euclid mission should provide all the necessary information for studying the cosmic evolution of the connectivity of the cosmic web around lensing clusters using multipolar moments and probing unique signatures of, for example, baryons and warm dark matter.

Starting from a very accurate model for density-in-cells statistics of dark matter based on large deviation theory, a bias model for the tracer density in spheres is formulated. It adopts a mean bias relation based on a quadratic bias model to relate the log-densities of dark matter to those of mass-weighted dark haloes in real and redshift space. The validity of the parametrised bias model is established using a parametrisation-independent extraction of the bias function. This average bias model is then combined with the dark matter PDF, neglecting any scatter around it: it nevertheless yields an excellent model for densities-in-cells statistics of mass tracers that is parametrised in terms of the underlying dark matter variance and three bias parameters. The procedure is validated on measurements of both the one and two point statistics of subhalo densities in the state-of-the-art Horizon Run 4 simulation showing excellent agreement for measured dark matter variance and bias parameters. Finally, it is demonstrated that this formalism allows for a joint estimation of the nonlinear dark matter variance and the bias parameters using solely the statistics of subhaloes. Having verified that galaxy counts in hydrodynamical simulations sampled on a scale of 10 Mpc/h closely resemble those of subhaloes, this work provides important steps towards making theoretical predictions for density-in-cells statistics applicable to upcoming galaxy surveys like Euclid or WFIRST.

Simple parameter-free analytic bias functions for the two-point correlation of densities in spheres at large separation are presented. These bias functions generalize the so-called Kaiser bias to the mildly non-linear regime for arbitrary density contrasts. The derivation is carried out in the context of large deviation statistics while relying on the spherical collapse model. A logarithmic transformation provides a saddle approximation which is valid for the whole range of densities and shown to be accurate against the 30 Gpc cube state-of-the-art Horizon Run 4 simulation. Special configurations of two concentric spheres that allow to identify peaks are employed to obtain the conditional bias and a proxy to BBKS extrema correlation functions. These analytic bias functions should be used jointly with extended perturbation theory to predict two-point clustering statistics as they capture the non-linear regime of structure formation at the percent level down to scales of about 10 Mpc/h at redshift 0. Conversely, the joint statistics also provide us with optimal dark matter two-point correlation estimates which can be applied either universally to all spheres or to a restricted set of biased (over- or underdense) pairs. Based on a simple fiducial survey, this estimator is shown to perform five times better than usual two-point function estimators. Extracting more information from correlations of different types of objects should prove essential in the context of upcoming surveys like Euclid, DESI, PFS or LSST.

The analytical formalism to obtain the probability distribution functions (PDFs) of spherically-averaged cosmic densities and velocity divergences in the mildly non-linear regime is presented. A large-deviation principle is applied to those cosmic fields assuming their most likely dynamics in spheres is set by the spherical collapse model. We validate our analytical results using state-of-the-art dark matter simulations with a phase-space resolved velocity field finding a 2% percent level agreement for a wide range of velocity divergences and densities in the mildly nonlinear regime (~10Mpc/h at redshift zero), usually inaccessible to perturbation theory. From the joint PDF of densities and velocity divergences measured in two concentric spheres, we extract with the same accuracy velocity profiles and conditional velocity PDF subject to a given over/under-density which are of interest to understand the non-linear evolution of velocity flows. Both PDFs are used to build a simple but accurate maximum likelihood estimators for the redshift evolution of the variance of both the density and velocity divergence fields, which have smaller relative errors than their sample variances when non-linearities appear. Given the dependence of the velocity divergence on the growth rate, there is a significant gain in using the full knowledge of both PDFs to derive constraints on the equation of state of dark energy. Thanks to the insensitivity of the velocity divergence to bias, its PDF can be used to obtain unbiased constraints on the growth of structures ($\sigma_8$,f) or it can be combined with the galaxy density PDF to extract bias parameters.

These notes are based on a lecture given at the 2016 Euclid Summer School in Narbonne. I will first give a quick overview of the concept of nuisance parameters in the context of large galaxy surveys. The second part will examine the case study of intrinsic alignments, a potential important contamination of weak lensing observables.

In the conference presentation we have reviewed the theory of non-Gaussian geometrical measures for the 3D Cosmic Web of the matter distribution in the Universe and 2D sky data, such as Cosmic Microwave Background (CMB) maps that was developed in a series of our papers. The theory leverages symmetry of isotropic statistics such as Minkowski functionals and extrema counts to develop post- Gaussian expansion of the statistics in orthogonal polynomials of invariant descriptors of the field, its first and second derivatives. The application of the approach to 2D fields defined on a spherical sky was suggested, but never rigorously developed. In this paper we present such development treating effects of the curvature and finiteness of the spherical space $S_2$ exactly, without relying on the flat-sky approximation. We present Minkowski functionals, including Euler characteristic and extrema counts to the first non-Gaussian correction, suitable for weakly non-Gaussian fields on a sphere, of which CMB is the prime example.

Intrinsic galaxy shape and angular momentum alignments can arise in cosmological large-scale structure due to tidal interactions or galaxy formation processes. Cosmological hydrodynamical simulations have recently come of age as a tool to study these alignments and their contamination to weak gravitational lensing. We probe the redshift and luminosity evolution of intrinsic alignments in Horizon-AGN between z=0 and z=3 for galaxies with an r-band absolute magnitude of <-20. Alignments transition from being radial at low redshifts and high luminosities, dominated by the contribution of ellipticals, to being tangential at high redshift and low luminosities, where discs dominate the signal. This cannot be explained by the evolution of the fraction of ellipticals and discs alone: intrinsic evolution in the amplitude of alignments is necessary. The alignment amplitude of elliptical galaxies alone is smaller in amplitude by a factor of ~2, but has similar luminosity and redshift evolution as in current observations and in the nonlinear tidal alignment model at projected separations of > 1 Mpc. Alignments of discs are null in projection and consistent with current low redshift observations. The combination of the two populations yields an overall amplitude a factor of ~4 lower than observed alignments of luminous red galaxies with a steeper luminosity dependence. The restriction on accurate galaxy shapes implies that the galaxy population in the simulation is complete only to an r-band absolute magnitude of <-20. Higher resolution simulations will be necessary to avoid extrapolation of the intrinsic alignment predictions to the range of luminosities probed by future surveys.

We present a first principle approach to obtain analytical predictions for spherically-averaged cosmic densities in the mildly non-linear regime that go well beyond what is usually achieved by standard perturbation theory. A large deviation principle allows us to compute the leading-order cumulants of average densities in concentric cells. In this symmetry, the spherical collapse model leads to cumulant generating functions that are robust for finite variances and free of critical points when logarithmic density transformations are implemented. They yield in turn accurate density probability distribution functions (PDFs) from a straightforward saddle-point approximation valid for all density values. Based on this easy-to-implement modification, explicit analytic formulas for the evaluation of the one- and two-cell PDF are provided. The theoretical predictions obtained for the PDFs are accurate to a few percent compared to the numerical integration, regardless of the density under consideration and in excellent agreement with N-body simulations for a wide range of densities. This formalism should prove valuable for accurately probing the quasi-linear scales of low redshift surveys for arbitrary primordial power spectra.

In order to quantify the error budget in the measured probability distribution functions of cell densities, the two-point statistics of cosmic densities in concentric spheres is investigated. Bias functions are introduced as the ratio of their two-point correlation function to the two-point correlation of the underlying dark matter distribution. They describe how cell densities are spatially correlated. They are computed here via the so-called large deviation principle in the quasi-linear regime. Their large-separation limit is presented and successfully compared to simulations for density and density slopes: this regime is shown to be rapidly reached allowing to get sub-percent precision for a wide range of densities and variances. The corresponding asymptotic limit provides an estimate of the cosmic variance of standard concentric cell statistics applied to finite surveys. More generally, no assumption on the separation is required for some specific moments of the two-point statistics, for instance when predicting the generating function of cumulants containing any powers of concentric densities in one location and one power of density at some arbitrary distance from the rest. This exact "one external leg" cumulant generating function is used in particular to probe the rate of convergence of the large-separation approximation.

The recently published analytic probability density function for the mildly non-linear cosmic density field within spherical cells is used to build a simple but accurate maximum likelihood estimate for the redshift evolution of the variance of the density, which, as expected, is shown to have smaller relative error than the sample variance. This estimator provides a competitive probe for the equation of state of dark energy, reaching a few percent accuracy on wp and wa for a Euclid-like survey. The corresponding likelihood function can take into account the configuration of the cells via their relative separations. A code to compute one-cell density probability density functions for arbitrary initial power spectrum, top-hat smoothing and various spherical collapse dynamics is made available online so as to provide straightforward means of testing the effect of alternative dark energy models and initial power-spectra on the low-redshift matter distribution.

The Lagrangian peaks of a 1D cosmological random field representing dark matter are used as a proxy for a catalogue of biased tracers in order to investigate the small-scale exclusion in the two-halo term. The two-point correlation function of peaks of a given height is numerically estimated and analytical approximations that are valid inside the exclusion zone are derived. The resulting power spectrum of these tracers is investigated and shows clear deviations from Poisson noise at low frequencies. On large scales, the convergence of a perturbative bias expansion is discussed. Finally, we go beyond Gaussian statistics for the initial conditions and investigate the subsequent evolution of the two-point clustering of peaks through their Zel'dovich ballistic displacement, to clarify how exclusion effects mix up with scale-dependencies induced by nonlinear gravitational evolution. While the expected large-scale separation limit is recovered, significant deviations are found in the exclusion zone that tends in particular to be reduced at later times. Even though these findings apply to the clustering of one-dimensional tracers, they provide useful insights into halo exclusion and its impact on the two-halo term.

The intrinsic alignments of galaxies are recognised as a contaminant to weak gravitational lensing measurements. In this work, we study the alignment of galaxy shapes and spins at low redshift ($z\sim 0.5$) in Horizon-AGN, an adaptive-mesh-refinement hydrodynamical cosmological simulation box of 100 Mpc/h a side with AGN feedback implementation. We find that spheroidal galaxies in the simulation show a tendency to be aligned radially towards over-densities in the dark matter density field and other spheroidals. This trend is in agreement with observations, but the amplitude of the signal depends strongly on how shapes are measured and how galaxies are selected in the simulation. Disc galaxies show a tendency to be oriented tangentially around spheroidals in three-dimensions. While this signal seems suppressed in projection, this does not guarantee that disc alignments can be safely ignored in future weak lensing surveys. The shape alignments of luminous galaxies in Horizon-AGN are in agreement with observations and other simulation works, but we find less alignment for lower luminosity populations. We also characterize the systematics of galaxy shapes in the simulation and show that they can be safely neglected when measuring the correlation of the density field and galaxy ellipticities.

The geometry of the cosmic web drives in part the spin acquisition of galaxies. This can be explained in a Lagrangian framework, by identifying the specific long-wavelength correlations within the primordial Gaussian random field which are relevant to spin acquisition. Tidal Torque Theory is revisited in the context of such anisotropic environments, biased by the presence of a filament within a wall. The point process of filament-type saddles represents it most efficiently. The constrained misalignment between the tidal and the inertia tensors in the vicinity of filament-type saddles simply explains the distribution of spin directions. This misalignment implies in particular an azimuthal orientation for the spins of more massive galaxies and a spin alignment with the filament for less massive galaxies. This prediction is found to be in qualitative agreement with measurements in Gaussian random fields and N-body simulations. It relates the transition mass to the geometry of the saddle, and accordingly predicts its measured scaling with the mass of non-linearity. Implications for galaxy formation and weak lensing are briefly discussed, as is the dual theory of spin alignments in walls.

The intrinsic alignment of galaxy shapes (by means of their angular momentum) and their cross-correlation with the surrounding dark matter tidal field are investigated using the 160 000, z=1.2 synthetic galaxies extracted from the high-resolution cosmological hydrodynamical simulation Horizon-AGN. One- and two-point statistics of the spin of the stellar component are measured as a function of mass and colour. For the low-mass galaxies, this spin is locally aligned with the tidal field `filamentary' direction while, for the high-mass galaxies, it is perpendicular to both filaments and walls. The bluest galaxies of our synthetic catalog are more strongly correlated with the surrounding tidal field than the reddest galaxies, and this correlation extends up to 10 Mpc/h comoving distance. We also report a correlation of the projected ellipticities of blue, intermediate mass galaxies on a similar scale at a level of 10^(-4) which could be a concern for cosmic shear measurements. We do not report any measurable intrinsic alignments of the reddest galaxies of our sample. This work is a first step toward the use of very realistic catalog of synthetic galaxies to evaluate the contamination of weak lensing measurement by the intrinsic galactic alignments.

In the context of count-in-cells statistics, the joint probability distribution of the density in two concentric spherical shells is predicted from first first principle for sigmas of the order of one. The agreement with simulation is found to be excellent. This statistics allows us to deduce the conditional one dimensional probability distribution function of the slope within under dense (resp. overdense) regions, or of the density for positive or negative slopes. The former conditional distribution is likely to be more robust in constraining the cosmological parameters as the underlying dynamics is less evolved in such regions. A fiducial dark energy experiment is implemented on such counts derived from Lambda-CDM simulations.

The kinematic analysis of dark matter and hydrodynamical simulations suggests that the vorticity in large-scale structure is mostly confined to, and predominantly aligned with their filaments, with an excess of probability of 20 per cent to have the angle between vorticity and filaments direction lower than 60 degrees relative to random orientations. The cross sections of these filaments are typically partitioned into four quadrants with opposite vorticity sign, arising from multiple flows, originating from neighbouring walls. The spins of halos embedded within these filaments are consistently aligned with this vorticity for any halo mass, with a stronger alignment for the most massive structures up to an excess of probability of 165 per cent. On large scales, adiabatic/cooling hydrodynamical simulations display the same vorticity in the gas as in the dark matter. The global geometry of the flow within the cosmic web is therefore qualitatively consistent with a spin acquisition for smaller halos induced by this large-scale coherence, as argued in Codis et al. (2012). In effect, secondary anisotropic infall (originating from the vortex-rich filament within which these lower-mass halos form) dominates the angular momentum budget of these halos. The transition mass from alignment to orthogonality is related to the size of a given multi-flow region with a given polarity. This transition may be reconciled with the standard tidal torque theory if the latter is augmented so as to account for the larger scale anisotropic environment of walls and filaments.

Filaments of the cosmic web drive spin acquisition of disc galaxies. The point process of filament-type saddle represent best this environment and can be used to revisit the Tidal Torque Theory in the context of an anisotropic peak (saddle) background split. The constrained misalignment between the tidal tensor and the Hessian of the density field generated in the vicinity of filament saddle points simply explains the corresponding transverse and longitudinal point-reflection symmetric geometry of spin distribution. It predicts in particular an azimuthal orientation of the spins of more massive galaxies and spin alignment with the filament for less massive galaxies. Its scale dependence also allows us to relate the transition mass corresponding to the alignment of dark matter halos spin relative to the direction of their neighboring filament to this geometry, and to predict accordingly it s scaling with the mass of non linearity, as was measured in simulations.

Intrinsic alignments are believed to be a major source of systematics for future generation of weak gravitational lensing surveys like Euclid or LSST. Direct measurements of the alignment of the projected light distribution of galaxies in wide field imaging data seem to agree on a contamination at a level of a few per cent of the shear correlation functions, although the amplitude of the effect depends on the population of galaxies considered. Given this dependency, it is difficult to use dark matter-only simulations as the sole resource to predict and control intrinsic alignments. We report here estimates on the level of intrinsic alignment in the cosmological hydrodynamical simulation Horizon-AGN that could be a major source of systematic errors in weak gravitational lensing measurements. In particular, assuming that the spin of galaxies is a good proxy for their ellipticity, we show how those spins are spatially correlated and how they couple to the tidal field in which they are embedded. We also present theoretical calculations that illustrate and qualitatively explain the observed signals.

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

In the context of upcoming large-scale structure surveys such as Euclid, it is of prime importance to quantify the effect of peculiar velocities on geometric probes. Hence the formalism to compute in redshift space the geometrical and topological one-point statistics of mildly non-Gaussian 2D and 3D cosmic fields is developed. Leveraging the partial isotropy of the target statistics, the Gram-Charlier expansion of the joint probability distribution of the field and its derivatives is reformulated in terms of the corresponding anisotropic variables. In particular, the cosmic non-linear evolution of the Minkowski functionals, together with the statistics of extrema are investigated in turn for 3D catalogues and 2D slabs. The amplitude of the non-Gaussian redshift distortion correction is estimated for these geometric probes. In 3D, gravitational perturbation theory is implemented in redshift space to predict the cosmic evolution of all relevant Gram-Charlier coefficients. Applications to the estimation of the cosmic parameters sigma(z) and beta=f/b1 from upcoming surveys is discussed. Such statistics are of interest for anisotropic fields beyond cosmology.

The joint probability distribution function (PDF) of the density within multiple concentric spherical cells is considered. It is shown how its cumulant generating function can be obtained at tree order in perturbation theory as the Legendre transform of a function directly built in terms of the initial moments. In the context of the upcoming generation of large-scale structure surveys, it is conjectured that this result correctly models such a function for finite values of the variance. Detailed consequences of this assumption are explored. In particular the corresponding one-cell density probability distribution at finite variance is computed for realistic power spectra, taking into account its scale variation. It is found to be in agreement with $\Lambda$-CDM simulations at the few percent level for a wide range of density values and parameters. Related explicit analytic expansions at the low and high density tails are given. The conditional (at fixed density) and marginal probability of the slope -- the density difference between adjacent cells -- and its fluctuations is also computed from the two-cells joint PDF; it also compares very well to simulations, in particular in under-dense regions, with a significant reduced cosmic scatter compared to over-dense regions. It is emphasized that this could prove useful when studying the statistical properties of voids as it can serve as a statistical indicator to test gravity models and/or probe key cosmological parameters.

We present a specific prescription for the calculation of cosmological power spectra, exploited here at two-loop order in perturbation theory (PT), based on the multi-point propagator expansion. In this approach power spectra are constructed from the regularized expressions of the propagators that reproduce both the resummed behavior in the high-k limit and the standard PT results at low-k. With the help of N-body simulations, we show that such a construction gives robust and accurate predictions for both the density power spectrum and the correlation function at percent-level in the weakly non-linear regime. We then present an algorithm that allows accelerated evaluations of all the required diagrams by reducing the computational tasks to one-dimensional integrals. This is achieved by means of pre-computed kernel sets defined for appropriately chosen fiducial models. The computational time for two-loop results is then reduced from a few minutes, with the direct method, to a few seconds with the fast one. The robustness and applicability of this method are tested against the power spectrum cosmic emulator from which a wide variety of cosmological models can be explored. The fortran program with which direct and fast calculations of power spectra can be done, RegPT, is publicly released as part of this paper.

We investigate the alignment of the spin of dark matter halos relative (i) to the surrounding large-scale filamentary structure, and (ii) to the tidal tensor eigenvectors using the Horizon 4pi dark matter simulation which resolves over 43 million dark matter halos at redshift zero. We detect a clear mass transition: the spin of dark matter halos above a critical mass tends to be perpendicular to the closest filament, and aligned with the intermediate axis of the tidal tensor, whereas the spin of low-mass halos is more likely to be aligned with the closest filament. Furthermore, this critical mass of 5 10^12 is redshift-dependent and scales as (1+z)^-2.5. We propose an interpretation of this signal in terms of large-scale cosmic flows. In this picture, most low-mass halos are formed through the winding of flows embedded in misaligned walls; hence they acquire a spin parallel to the axis of the resulting filaments forming at the intersection of these walls. On the other hand, more massive halos are typically the products of later mergers along such filaments, and thus they acquire a spin perpendicular to this direction when their orbital angular momentum is converted into spin. We show that this scenario is consistent with both the measured excess probabilities of alignment w.r.t. the eigen-directions of the tidal tensor, and halo merger histories. On a more qualitative level, it also seems compatible with 3D visualization of the structure of the cosmic web as traced by "smoothed" dark matter simulations or gas tracer particles. Finally, it provides extra support to the disc forming paradigm presented by Pichon et al (2011) as it extends it by characterizing the geometry of secondary infall at high redshift.