Galaxy clusters are connected at their peripheries to the large scale structures by cosmic filaments that funnel accreting material. These filamentary structures are studied to investigate both environment-driven galaxy evolution and structure formation and evolution. In the present work, we probe in a statistical manner the azimuthal distribution of galaxies around clusters as a function of the cluster-centric distance, the cluster richness, and the galaxy activity (star-forming or passive).We perform a harmonic decomposition in large photometric galaxy catalogue around 6400 SDSS clusters with masses M >1e14 solar masses, in the redshift range of 0.1< z <0.3. The same analysis is performed on the mock galaxy catalogue from the light-cone of Magneticum hydrodynamical simulation. We use the multipole analysis to quantify asymmetries in the 2-D galaxy distribution. In the inner cluster regions at R <2 R500, we confirm that the galaxy distribution traces an ellipsoidal shape, which is more pronounced for richest clusters. In the clusters' outskirts (R= [2-8]R500), filamentary patterns are detected in harmonic space with a mean angular scale m_mean= 4.2+/-0.1. Massive clusters seem to have a larger number of connected filaments than low massive ones. We also find that passive galaxies appear to better trace the filamentary structures around clusters, even if the contribution of SF ones tend to increase with the cluster-centric distance, suggesting a gradient of galaxy activity in filaments around clusters.

We present a statistical study of the filamentary structures of the Cosmic Web in the large hydro-dynamical simulations TNG300-1 and MAGNETICUM at redshift z=0. We focus on the radial distribution of the galaxy density around filaments detected using the Discrete Persistent Structure Extractor (DisPerSE). We show that the average profile of filaments presents an excess of galaxy density (> 5 sigma) up to radial distances of 27 Mpc from the core. The relation between galaxy density and the length of filaments is further investigated showing that short (L_f < 9 Mpc) and long (L_f > 20 Mpc) filaments are two statistically different populations. Short filaments are puffier, denser and more connected to massive objects, whereas long filaments are thinner, less dense and more connected to less massive structures. These two populations trace different environments and may correspond to bridges of matter between over-dense structures (short filaments) and to the cosmic filaments shaping the skeleton of the Cosmic Web (long filaments). Through MCMC explorations, we find that the density profiles of both short and long filaments can be described by the same empirical models (GNFW, beta-model, a single and a double power law) with different and distinct sets of parameters.

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.

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.