Observations of galaxy clustering are made in redshift space, which results in distortions to the underlying isotropic distribution of galaxies. These redshift-space distortions (RSD) not only degrade important features of the matter density field, such as the baryonic acoustic oscillation (BAO) peaks, but also pose challenges for the theoretical modelling of observational probes. Here we introduce an iterative nonlinear reconstruction algorithm to remove RSD effects from galaxy clustering measurements, and assess its performance by using mock galaxy catalogues. The new method is found to be able to recover the real-space galaxy correlation function with an accuracy of $\sim1\%$, and restore the quadrupole accurately to $0$, on scales $s\gtrsim20\Mpch$. It also leads to an improvement in the reconstruction of the initial density field, which could help to accurately locate the BAO peaks. An `internal calibration' scheme is proposed to determine the values of cosmological parameters as a part of the reconstruction process, and possibilities to break parameter degeneracies are discussed. RSD reconstruction can offer a potential way to simultaneously extract the cosmological parameters, initial density field, real-space galaxy positions and large-scale peculiar velocity field (of the real Universe), making it an alternative to standard perturbative approaches in galaxy clustering analysis, bypassing the need for RSD modelling.

We perform an extensive review of the numerous studies and methods used to determine the total mass of the Milky Way. We group the various methods into seven broad classes, including: i) estimating Galactic escape velocity using high velocity objects; ii) measuring the rotation curve through terminal and circular velocities; iii) modeling halo stars, globular clusters and satellite galaxies with the Spherical Jeans equation and iv) with phase-space distribution functions; v) simulating and modeling the dynamics of stellar streams and their progenitors; vi) modeling the motion of the Milky Way, M31 and other distant satellites under the framework of Local Group timing argument; and vii) measurements made by linking the brightest Galactic satellites to their counterparts in simulations. For each class of methods, we introduce their theoretical and observational background, the method itself, the sample of available tracer objects, model assumptions, uncertainties, limits and the corresponding measurements that have been achieved in the past. Both the measured total masses within the radial range probed by tracer objects and the extrapolated virial masses are discussed and quoted. We also discuss the role of modern numerical simulations in terms of helping to validate model assumptions, understanding systematic uncertainties and calibrating the measurements. While measurements in the last two decades show a factor of two scatters, recent measurements using \textit{Gaia} DR2 data are approaching a higher precision. We end with a detailed discussion of future developments, especially as the size and quality of the observational data will increase tremendously with current and future surveys. In such cases, the systematic uncertainties will be dominant and thus will necessitate a much more rigorous testing and characterization of the various mass determination methods.

We analyse systems analogous to the Milky Way (MW) in the EAGLE cosmological hydrodynamics simulation in order to deduce the likely structure of the MW's dark matter halo. We identify MW-mass haloes in the simulation whose satellite galaxies have similar kinematics and spatial distribution to those of the bright satellites of the MW, specifically systems in which the majority of the satellites (8 out of 11) have nearly co-planar orbits that are also perpendicular to the central stellar disc. We find that the normal to the common orbital plane of the co-planar satellites is well aligned with the minor axis of the host dark matter halo, with a median misalignment angle of only $17.3^\circ$. Based on this result, we infer that the minor axis of the Galactic dark matter halo points towards $(l,b)=(182^\circ,-2^\circ)$, with an angular uncertainty at the 68 and 95 percentile confidence levels of 22$^\circ$ and 43$^\circ$ respectively. Thus, the inferred minor axis of the MW halo lies in the plane of the stellar disc. The halo, however, is not homologous and its flattening and orientation vary with radius. The inner parts of the halo are rounder than the outer parts and well-aligned with the stellar disc (that is the minor axis of the halo is perpendicular to the disc). Further out, the halo twists and the minor axis changes direction by $90^\circ$. This twist occurs over a very narrow radial range and reflects variations in the filamentary network along which mass was accreted into the MW.

Cosmic voids are a key component of the large-scale structure that contain a plethora of cosmological information. Typically, voids are identified from the underlying galaxy distribution, which is a biased tracer of the total matter field. Previous works have shown that 2D voids identified in weak lensing maps -- weak lensing voids -- correspond better to true underdense regions along the line of sight. In this work, we study how the properties of weak lensing voids depend on the choice of void finder, by adapting several popular void finders. We present and discuss the differences between identifying voids directly in the convergence maps, and in the distribution of weak lensing peaks. Particular effort has been made to test how these results are affected by galaxy shape noise, which is a dominant source of noise in weak lensing observations. By studying the signal-to-noise ratios (SNR) for the tangential shear profile of each void finder, we find that voids identified directly in the convergence maps have the highest SNR but are also the ones most affected by galaxy shape noise. Troughs are least affected by noise, but also have the lowest SNR. The tunnel algorithm, which identifies voids in the distribution of weak lensing peaks, represents a good compromise between finding a large tangential shear SNR and mitigating the effect of galaxy shape noise.

We determine the Milky Way (MW) mass profile inferred from fitting physically motivated models to the Gaia DR2 Galactic rotation curve and other data. Using various hydrodynamical simulations of MW-mass haloes, we show that the presence of baryons induces a contraction of the dark matter (DM) distribution in the inner regions, r<20 kpc. We provide an analytic expression that relates the baryonic distribution to the change in the DM halo profile. For our galaxy, the contraction increases the enclosed DM halo mass by factors of roughly 1.3, 2 and 4 at radial distances of 20, 8 and 1 kpc, respectively compared to an uncontracted halo. Ignoring this contraction results in systematic biases in the inferred halo mass and concentration. We provide a best-fitting contracted NFW halo model to the MW rotation curve that matches the data very well. The best-fit has a DM halo mass, $M_{200}^{\rm DM}=0.97_{-0.19}^{+0.24}\times10^{12} M_\odot$, and concentration before baryon contraction of $9.4_{-2.6}^{+1.9}$, which lie close to the median halo mass--concentration relation predicted in $\Lambda$CDM. The inferred total mass, $M_{200}^{\rm total}=1.08_{-0.14}^{+0.20} \times 10^{12} M_\odot$, is in good agreement with recent measurements. The model gives a MW stellar mass of $5.04_{-0.52}^{+0.43}\times10^{10} M_\odot$ and infers that the DM density at the Solar position is $\rho_{\odot}^{\rm DM}=8.8_{-0.5}^{+0.5}\times10^{-3} M_\odot \rm{pc}^{-3}\equiv0.33_{-0.02}^{+0.02}~\rm{GeV}~\rm{cm}^{-3}$. The rotation curve data can also be fitted with an uncontracted NFW halo model, but with very different DM and stellar parameters. The observations prefer the physically motivated contracted NFW halo, but the measurement uncertainties are too large to rule out the uncontracted NFW halo.

We study the orbital phase-space of dark matter (DM) halos in the AURIGA suite of cosmological hydrodynamics simulations of Milky Way analogues. We characterise halos by their spherical action distribution, a function of the specific angular momentum, and the radial action, of the DM particles. By comparing DM-only and hydrodynamical simulations of the same halos, we investigate the contraction of DM halos caused by the accumulation of baryons at the centre. We find a small systematic suppression of the radial action in the DM halos of the hydrodynamical simulations, suggesting that the commonly used adiabatic contraction approximation can result in an underestimate of the density by ~ 8%. We apply an iterative algorithm to contract the AURIGA DM halos given a baryon density profile and halo mass, recovering the true contracted DM profiles with an accuracy of ~15%, that reflects halo-to-halo variation. Using this algorithm, we infer the total mass profile of the Milky Way's contracted DM halo. We derive updated values for the key astrophysical inputs to DM direct detection experiments: the DM density and velocity distribution in the Solar neighbourhood.

We investigate if, for a fixed number density of targets and redshift, there is an optimal way to select a galaxy sample in order to measure the baryon acoustic oscillation (BAO) scale, which is used as a standard ruler to constrain the cosmic expansion. Using the mock galaxy catalogue built by Smith et al. in the Millennium-XXL N-body simulation with a technique to assign galaxies to dark matter haloes based on halo occupation distribution modelling, we consider the clustering of galaxies selected by luminosity, colour and local density. We assess how well the BAO scale can be extracted by fitting a template to the power spectrum measured for each sample. We find that the BAO peak position is recovered equally well for samples defined by luminosity or colour, while there is a bias in the BAO scale recovered for samples defined by density. The BAO position is contracted to smaller scales for the densest galaxy quartile and expanded to large scales for the two least dense galaxy quartiles. For fixed galaxy number density, density-selected samples have higher uncertainties in the recovered BAO scale than luminosity- or colour-selected samples.

Modifications to General Relativity (GR) often incorporate screening mechanisms in order to remain compatible with existing tests of gravity. The screening is less efficient in underdense regions, which suggests that cosmic voids can be a useful cosmological probe for constraining modified gravity models. In particular, weak lensing by voids has been proposed as a promising test of such theories. Usually, voids are identified from galaxy distributions, making them biased tracers of the underlying matter field. An alternative approach is to study voids identified in weak lensing maps -- weak lensing voids -- which have been shown to better correspond to true underdense regions. In this paper, we study the ability of weak lensing voids to detect the signatures of modified gravity. Focusing on the void abundance and weak lensing profiles, we find that both statistics are sensitive probes of gravity. These are quantified in terms of the signal-to-noise ratios (SNR) with which an LSST-like survey will be able to distinguish between different gravity models. We find that the tangential shear profiles of weak lensing voids are considerably better than galaxy voids at this, though voids have somewhat lower SNR than weak lensing peaks. The abundances of voids and peaks have respectively $\rm{SNR} = 50$ and $70$ for a popular class of modified gravity in an LSST-like survey.

We study the statistics of weak lensing convergence peaks, such as their abundance and two-point correlation function (2PCF), for a wide range of cosmological parameters $\Omega_m$ and $\sigma_8$ within the standard $\Lambda$CDM paradigm, focusing on intermediate-height peaks with signal-to-noise ratio (SNR) of $1.5$ to $3.5$. We find that the cosmology dependence of the peak abundance can be described by a one-parameter fitting formula that is accurate to within $\sim3\%$. The peak 2PCFs are shown to feature a self-similar behaviour: if the peak separation is rescaled by the mean inter-peak distance, catalogues with different minimum peak SNR values have identical clustering, which suggests that the peak abundance and clustering are closely interconnected. A simple fitting model for the rescaled 2PCF is given, which together with the peak abundance model above can predict peak 2PCFs with an accuracy better than $\sim5\%$. The abundance and 2PCFs for intermediate peaks have very different dependencies on $\Omega_m$ and $\sigma_8$, implying that their combination can be used to break the degeneracy between these two parameters.

We present the DTFE public software, a code for reconstructing fields from a discrete set of samples/measurements using the maximum of information contained in the point distribution. The code is written in C++ using the CGAL library and is parallelized using OpenMP. The software was designed for the analysis of cosmological data but can be used in other fields where one must interpolate quantities given at a discrete point set. The software comes with a wide suite of options to facilitate the analysis of 2- and 3-dimensional data and of both numerical simulations and galaxy redshift surveys. For comparison purposes, the code also implements the TSC and SPH grid interpolation methods. The code comes with an extensive user guide detailing the program options, examples and the inner workings of the code. The DTFE public software and further information can be found at http://www.astro.rug.nl/~voronoi/DTFE/dtfe.html .

We introduce the LOCal Universe Screening Test Suite (LOCUSTS) project, an effort to create 'screening maps' in the nearby Universe to identify regions in our neighbourhood which are screened, i.e., regions where deviations from General Relativity (GR) are suppressed, in various modified gravity (MG) models. In these models, deviations from the GR force law are often stronger for smaller astrophysical objects, making them ideal test beds of gravity in the local Universe. However, the actual behaviour of the modified gravity force also depends on the environment of the objects, and to make accurate predictions one has to take the latter into account. This can be done approximately using luminous objects in the local Universe as tracers of the underlying dark matter field. Here, we propose a new approach that takes advantage of state-of-the-art Bayesian reconstruction of the mass distribution in the Universe, which allows us to solve the modified gravity equations and predict the screening effect more accurately. This is the first of a series of works, in which we present our methodology and some qualitative results of screening for a specific MG model, $f(R)$ gravity. Applications to test models using observations and extensions to other classes of models will be studied in future works.

We study the formation of planes of dwarf galaxies around Milky Way (MW)-mass haloes in the EAGLE galaxy formation simulation. We focus on satellite systems similar to the one in the MW: spatially thin or with a large fraction of members orbiting in the same plane. To characterise the latter, we introduce a robust method to identify the subsets of satellites that have the most co-planar orbits. Out of the 11 MW classical dwarf satellites, 8 have highly clustered orbital planes whose poles are contained within a $22^\circ$ opening angle centred around $(l,b)=(182^\circ,-2^\circ)$. This configuration stands out when compared to both isotropic and typical $\Lambda$CDM satellite distributions. Purely flattened satellite systems are short-lived chance associations and persist for less than $1~\rm{Gyr}$. In contrast, satellite subsets that share roughly the same orbital plane are longer lived, with half of the MW-like systems being at least $4~\rm{Gyrs}$ old. On average, satellite systems were flatter in the past, with a minimum in their minor-to-major axes ratio about $9~\rm{Gyrs}$ ago, which is the typical infall time of the classical satellites. MW-like satellite distributions have on average always been flatter than the overall population of satellites in MW-mass haloes and, in particular, they correspond to systems with a high degree of anisotropic accretion of satellites. We also show that torques induced by the aspherical mass distribution of the host halo channel some satellite orbits into the host's equatorial plane, enhancing the fraction of satellites with co-planar orbits. In fact, the orbital poles of co-planar satellites are tightly aligned with the minor axis of the host halo.

We investigate the alignment of galaxies and haloes relative to cosmic web filaments using the EAGLE hydrodynamical simulation. We identify filaments by applying the NEXUS+ method to the mass distribution and the Bisous formalism to the galaxy distribution. Both web finders return similar filamentary structures that are well aligned and that contain comparable galaxy populations. EAGLE haloes have an identical spin alignment with filaments as their counterparts in dark matter only simulations: a complex mass dependent trend with low mass haloes spinning preferentially parallel to and high mass haloes spinning preferentially perpendicular to filaments. In contrast, galaxy spins do not show such a spin transition and have a propensity for perpendicular alignments at all masses, with the degree of alignment being largest for massive galaxies. This result is valid for both NEXUS+ and Bisous filaments. When splitting by morphology, we find that elliptical galaxies show a stronger orthogonal spin--filament alignment than spiral galaxies of similar mass. The same is true of their haloes, with the host haloes of elliptical galaxies having a larger degree of orthogonal alignment than the host haloes of spirals. Due to the misalignment between galaxy shape and spin, galaxy minor axes are oriented differently with filaments than galaxy spins. We find that the galaxies whose minor axis is perpendicular to a filament are much better aligned with their host haloes. This suggests that many of the same physical processes determine both the galaxy--filament and the galaxy--halo alignments.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

We present and apply a method to infer the mass of the Milky Way (MW) by comparing the dynamics of MW satellites to those of model satellites in the EAGLE cosmological hydrodynamics simulations. A distribution function (DF) for galactic satellites is constructed from EAGLE using specific angular momentum and specific energy, which are scaled so as to be independent of host halo mass. In this 2-dimensional space, the orbital properties of satellite galaxies vary according to the host halo mass. The halo mass can be inferred by calculating the likelihood that the observed satellite population is drawn from this DF. Our method is robustly calibrated on mock EAGLE systems. We validate it by applying it to the completely independent suite of 30 AURIGA high-resolution simulations of MW-like galaxies: the method accurately recovers their true mass and associated uncertainties. We then apply it to ten classical satellites of the MW with 6D phase-space measurements, including updated proper motions from the GAIA satellite. The mass of the MW is estimated to be $M_{200}^{\textnormal{MW}}=1.17_{-0.15}^{+0.21}\times10^{12}M_{\odot}$ (68\% confidence limits). We combine our total mass estimate with recent mass estimates in the inner regions of the Galaxy to infer an inner dark matter (DM) mass fraction $M^\textnormal{DM}(<20~\rm{kpc})/M^\textnormal{DM}_{200}=0.12$ which is typical of ${\sim}10^{12}M_{\odot}$ $\Lambda$CDM haloes in hydrodynamical galaxy formation simulations. Assuming an NFW profile, this is equivalent to a halo concentration of $c_{200}^{\textnormal{MW}}=10.9^{+2.6}_{-2.0}$.

The Milky Way (MW) offers a uniquely detailed view of galactic structure and is often regarded as a prototypical spiral galaxy. But recent observations indicate that the MW is atypical: it has an undersized supermassive black hole at its centre; it is surrounded by a very low mass, excessively metal-poor stellar halo; and it has an unusually large nearby satellite galaxy, the Large Magellanic Cloud (LMC). Here we show that the LMC is on a collision course with the MW with which it will merge in $2.4^{+1.2}_{-0.8}$ Gyrs (68% confidence level). This catastrophic and long-overdue event will restore the MW to normality. Using the EAGLE galaxy formation simulation, we show that, as a result of the merger, the central supermassive black hole will increase in mass by up to a factor of 8. The Galactic stellar halo will undergo an equally impressive transformation, becoming 5 times more massive. The additional stars will come predominantly from the disrupted LMC, but a sizeable number will be ejected onto the halo from the stellar disc. The post-merger stellar halo will have the median metallicity of the LMC, [Fe/H]=-0.5 dex, which is typical of other galaxies of similar mass to the MW. At the end of this exceptional event, the MW will become a true benchmark for spiral galaxies, at least temporarily.

We study cosmic voids in the normal-branch Dvali-Gabadadze-Porrati (nDGP) braneworld models, which are representative of a class of modified gravity theories where deviations from General Relativity are usually hidden by the Vainshtein screening in high-density environments. This screening is less efficient away from these environments, which makes voids ideally suited for testing this class of models. We use N-body simulations of $\Lambda$-cold dark matter ($\Lambda$CDM) and nGDP universes, where dark matter haloes are populated with mock galaxies that mimic the clustering and number densities of the BOSS CMASS galaxy sample. We measure the force, density and weak lensing profiles around voids identified with six different algorithms. Compared to $\Lambda$CDM, voids in nDGP are more under-dense due to the action of the fifth force that arises in these models, which leads to a faster evacuation of matter from voids. This leaves an imprint on the weak lensing tangential shear profile around nDGP voids, an effect that is particularly strong for 2D underdensities that are identified in the plane-of-the-sky. We make predictions for the feasibility of distinguishing between nDGP and $\Lambda$CDM using void lensing in upcoming large-scale surveys such as LSST and Euclid. We compare with the analysis of voids in chameleon gravity theories and find that the weak lensing signal for 3D voids is similar to nDGP, whereas for 2D voids the differences with $\Lambda$CDM are much stronger for the chameleon gravity case, a direct consequence of the different screening mechanisms operating in these theories.

The reconstruction of the initial conditions of the Universe is an important topic in cosmology, particularly in the context of sharpening the measurement of the baryon acoustic oscillation (BAO) peak. Nonlinear reconstruction algorithms developed in recent years, when applied to late-time matter fields, can recover to a substantial degree the initial density distribution, however, when applied to sparse tracers of the matter field, the performance is poorer. In this paper we apply the Shi et al. non-linear reconstruction method to biased tracers in order to establish what factors affect the reconstruction performance. We find that grid resolution, tracer number density and mass assignment scheme all have a significant impact on the performance of our reconstruction method, with triangular-shaped-cloud (TSC) mass assignment and a grid resolution of ${\sim}1{-}2h^{-1}$ Mpc being the optimal choice. We also show that our method can be easily adapted to include generic tracer biases up to quadratic order in the reconstruction formalism. Applying the reconstruction to halo and galaxy samples with a range of tracer number densities, we find that the linear bias is by far the most important bias term, while including nonlocal and nonlinear biases only leads to marginal improvements on the reconstruction performance. Overall, including bias in the reconstruction substantially improves the recovery of BAO wiggles, down to $k\sim0.25~h\text{Mpc}^{-1}$ for tracer number densities between $2\times10^{-4}$ and $2\times10^{-3}~(h^{-1}\text{Mpc})^{-3}$.

The Milky Way's (MW) satellite population is a powerful probe of warm dark matter (WDM) models as the abundance of small substructures is very sensitive to the properties of the WDM particle. However, only a partial census of the MW's complement of satellite galaxies exists because surveys of the MW's close environs are incomplete both in depth and in sky coverage. We present a new Bayesian analysis that combines the sample of satellites recently discovered by the Dark Energy Survey (DES) with those found in the Sloan Digital Sky Survey (SDSS) to estimate the total satellite galaxy luminosity function down to $M_{\rm V}=0$. We find that there should be at least $124^{+40}_{-27}$ ($68\%$ CL, statistical error) satellites as bright or brighter than $M_{\rm V}=0$ within 300 kpc of the Sun, with only a weak dependence on MW halo mass. When it comes online the Large Synoptic Survey Telescope should detect approximately half of this population. We also show that WDM models infer the same number of satellites as in $\Lambda$CDM, which will allow us to rule out those models that produce insufficient substructure to be viable.

We investigate the alignment of haloes with the filaments of the cosmic web using an unprecedently large sample of dark matter haloes taken from the P-Millennium $\Lambda$CDM cosmological N-body simulation. We use the state-of-the-art NEXUS morphological formalism which, due to its multiscale nature, simultaneously identifies structures at all scales. We find strong and highly significant alignments, with both the major axis of haloes and their peculiar velocity tending to orient along the filament. However, the spin - filament alignment displays a more complex trend changing from preferentially parallel at low masses to preferentially perpendicular at high masses. This "spin flip" occurs at an average mass of $5\times10^{11}~h^{-1}M_\odot$. This mass increases with increasing filament diameter, varying by more than an order of magnitude between the thinnest and thickest filament samples. We also find that the inner parts of haloes have a spin flip mass that is several times smaller than that of the halo as a whole. These results confirm that recent accretion is responsible for the complex behaviour of the halo spin - filament alignment. Low-mass haloes mainly accrete mass along directions perpendicular to their host filament and thus their spins tend to be oriented along the filaments. In contrast, high-mass haloes mainly accrete along their host filaments and have their spins preferentially perpendicular to them. Furthermore, haloes located in thinner filaments are more likely to accrete along their host filaments than haloes of the same mass located in thicker filaments.

We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $\Lambda$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at $|b|>20$ degrees. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. that preserves the phase-space distribution of the model stars. These publicly available catalogues contain 5-parameter astrometry, radial velocities, multi-band photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that: i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A- and B-dwarf stars up to radii of ~15 kpc; and ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to: i) develop and test analysis methods for the Gaia DR2 data; ii) gauge the limitations and biases of the data and iii) interpret the data in the light of theoretical predictions from realistic $ab$ $initio$ simulations of galaxy formation in the $\Lambda$CDM cosmological model.

Cosmic voids are an important probe of large-scale structure that can constrain cosmological parameters and test cosmological models. We present a new paradigm for void studies: void detection in weak lensing convergence maps. This approach identifies objects that relate directly to our theoretical understanding of voids as underdensities in the total matter field and presents several advantages compared to the customary method of finding voids in the galaxy distribution. We exemplify this approach by identifying voids using the weak lensing peaks as tracers of the large-scale structure. We find self-similarity in the void abundance across a range of peak signal-to-noise selection thresholds. The voids obtained via this approach give a tangential shear signal up to $\sim40$ times larger than voids identified in the galaxy distribution.

The total number and luminosity function of the population of dwarf galaxies of the Milky Way (MW) provide important constraints on the nature of the dark matter and on the astrophysics of galaxy formation at low masses. However, only a partial census of this population exists because of the flux limits and restricted sky coverage of existing Galactic surveys. We combine the sample of satellites recently discovered by the Dark Energy Survey (DES) with the satellites found in Sloan Digital Sky Survey (SDSS) Data Release 9 (together these surveys cover nearly half the sky) to estimate the total luminosity function of satellites down to $M_{\rm V}=0$. We apply a new Bayesian inference method in which we assume that the radial distribution of satellites independently of absolute magnitude follows that of subhaloes selected according to their peak maximum circular velocity. We find that there should be at least $124^{+40}_{-27}$ (68 per cent CL, statistical error) satellites brighter than $M_{\rm V}=0$ within $300$ kpc of the Sun. As a result of our use of new data and better simulations, and a more robust statistical method, we infer a much smaller population of satellites than reported in previous studies using earlier SDSS data only; we also address an underestimation of the uncertainties in earlier work by accounting for stochastic effects. We find that the inferred number of faint satellites depends only weakly on the assumed mass of the MW halo and we provide scaling relations to extend our results to different assumed halo masses and outer radii. We predict that half of our estimated total satellite population of the MW should be detected by the Large Synoptic Survey Telescope. The code implementing our estimation method is available online.

We investigate the population of dwarf galaxies with stellar masses similar to the Large Magellanic Cloud (LMC) and M33 in the EAGLE galaxy formation simulation. In the field, galaxies reside in haloes with stellar-to-halo mass ratios of $1.03^{+0.50}_{-0.31}\times10^{-2}$ (68% confidence level); systems like the LMC, which have an SMC-mass satellite, reside in haloes about 1.3 times more massive, which suggests an LMC halo mass at infall, $M_{200}=3.4^{+1.8}_{-1.2}\times10^{11}M_\odot$ (68% confidence level). The colour distribution of dwarfs is bimodal, with the red galaxies ($g-r>0.6$) being mostly satellites. The fraction of red LMC-mass dwarfs is 15% for centrals, and for satellites this fraction increases rapidly with host mass: from 10% for satellites of Milky Way (MW)-mass haloes to nearly 90% for satellites of groups and clusters. The quenching timescale, defined as the time after infall when half of the satellites have acquired red colours, decreases with host mass from ${>}5$ Gyrs for MW-mass hosts to $2.5$ Gyrs for cluster mass hosts. The satellites of MW-mass haloes have higher star formation rates and bluer colours than field galaxies. This is due to enhanced star formation triggered by gas compression shortly after accretion. Both the LMC and M33 have enhanced recent star formation that could be a manifestation of this process. After infall into their MW-mass hosts, the $g-r$ colours of LMC-mass dwarfs become bluer for the first 2 Gyrs, after which they rapidly redden. LMC-mass dwarfs fell into their MW-mass hosts only relatively recently, with more than half having an infall time of less than 3.5 Gyrs.

We study the incidence of group and filamentary dwarf galaxy accretion into Milky Way (MW) mass haloes using two types of hydrodynamical simulations: EAGLE, which resolves a large cosmological volume, and the AURIGA suite, which are very high resolution zoom-in simulations of individual MW-sized haloes. The present-day 11 most massive satellites are predominantly (75%) accreted in single events, 14% in pairs and 6% in triplets, with higher group multiplicities being unlikely. Group accretion becomes more common for fainter satellites, with 60% of the top 50 satellites accreted singly, 12% in pairs, and 28% in richer groups. A group similar in stellar mass to the Large Magellanic Cloud (LMC) would bring on average 15 members with stellar mass larger than $10^4{~\rm M_\odot}$. Half of the top 11 satellites are accreted along the two richest filaments. The accretion of dwarf galaxies is highly anisotropic, taking place preferentially perpendicular to the halo minor axis, and, within this plane, preferentially along the halo major axis. The satellite entry points tend to be aligned with the present-day central galaxy disc and satellite plane, but to a lesser extent than with the halo shape. Dwarfs accreted in groups or along the richest filament have entry points that show an even larger degree of alignment with the host halo than the full satellite population. We also find that having most satellites accreted as a single group or along a single filament is unlikely to explain the MW disc of satellites.

We present a systematic comparison of several existing and new void finding algorithms, focusing on their potential power to test a particular class of modified gravity models - chameleon $f(R)$ gravity. These models deviate from standard General Relativity (GR) more strongly in low-density regions and thus voids are a promising venue to test them. We use Halo Occupation Distribution (HOD) prescriptions to populate haloes with galaxies, and tune the HOD parameters such that the galaxy two-point correlation functions are the same in both f(R) and GR models. We identify both 3D voids as well as 2D underdensities in the plane-of-the-sky to find the same void abundance and void galaxy number density profiles across all models, which suggests that they do not contain much information beyond galaxy clustering. However, the underlying void dark matter density profiles are significantly different, with f(R) voids being more underdense than GR ones, which leads to f(R) voids having a larger tangential shear signal than their GR analogues. We investigate the potential of each void finder to test f(R) models with near-future lensing surveys such as EUCLID and LSST. The 2D voids have the largest power to probe f(R) gravity, with a LSST analysis of tunnel (which is a new type of 2D underdensity introduced here) lensing distinguishing at 80 and 11$\sigma$ (statistical error) f(R) models with $|f_{R0}|=10^{-5}$ and $10^{-6}$ from GR.

The Vainshtein mechanism, present in many models of gravity, is very effective at screening dark matter halos such that the fifth force is negligible and general relativity is recovered within their Vainshtein radii. Vainshtein screening is independent of halo mass and environment, in contrast to e.g. chameleon screening, making it difficult to test. However, our previous studies have found that the dark matter particles in filaments, walls, and voids are not screened by the Vainshtein mechanism. We therefore investigate whether cosmic voids, identified as local density minima using a watershed technique, can be used to test models of gravity that exhibit Vainshtein screening. We measure density, velocity, and screening profiles of stacked voids in cosmological $N$-body simulations using both dark matter particles and dark matter halos as tracers of the density field. We find that the voids are completely unscreened, and the tangential velocity and velocity dispersion profiles of stacked voids show a clear deviation from $\Lambda$CDM at all radii. Voids have the potential to provide a powerful test of gravity on cosmological scales.

A theoretically interesting and practically important question in cosmology is the reconstruction of the initial density distribution provided a late-time density field. This is a long-standing question with a revived interest recently, especially in the context of optimally extracting the baryonic acoustic oscillation (BAO) signals from observed galaxy distributions. We present a new efficient method to carry out this reconstruction, which is based on numerical solutions to the nonlinear partial differential equation that governs the mapping between the initial Lagrangian and final Eulerian coordinates of particles in evolved density fields. This is motivated by numerical simulations of the quartic Galileon gravity model, which has similar equations that can be solved effectively by multigrid Gauss-Seidel relaxation. The method is based on mass conservation, and does not assume any specific cosmological model. Our test shows that it has a performance comparable to that of state-of-the-art algorithms which were very recently put forward in the literature, with the reconstructed density field over $\sim80\%$ ($50\%$) correlated with the initial condition at $k\lesssim0.6h/{\rm Mpc}$ ($1.0h/{\rm Mpc}$). With an example, we demonstrate that this method can significantly improve the accuracy of BAO reconstruction.

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

We estimate the systemic orbital kinematics of the Milky Way classical satellites and compare them with predictions from the \Lambda{} cold dark matter (\Lambda{}CDM) model derived from a semi-analytical galaxy formation model applied to high resolution cosmological N-body simulations. We find that the Galactic satellite system is atypical of \Lambda{}CDM systems. The subset of 10 Galactic satellites with proper motion measurements has a velocity anisotropy, \beta=-2.2$\pm$0.4, that lies in the 2.9% tail of the \Lambda{}CDM distribution. Individually, the Milky Way satellites have radial velocities that are lower than expected for their proper motions, with 9 out of the 10 having at most 20% of their orbital kinetic energy invested in radial motion. Such extreme values are expected in only 1.5% of \Lambda{}CDM satellites systems. This tangential motion excess is unrelated to the existence of a Galactic "disc of satellites". We present theoretical predictions for larger satellite samples that may become available as more proper motion measurements are obtained.

We describe the Bayesian BARCODE formalism that has been designed towards the reconstruction of the Cosmic Web in a given volume on the basis of the sampled galaxy cluster distribution. Based on the realization that the massive compact clusters are responsible for the major share of the large scale tidal force field shaping the anisotropic and in particular filamentary features in the Cosmic Web. Given the nonlinearity of the constraints imposed by the cluster configurations, we resort to a state-of-the-art constrained reconstruction technique to find a proper statistically sampled realization of the original initial density and velocity field in the same cosmic region. Ultimately, the subsequent gravitational evolution of these initial conditions towards the implied Cosmic Web configuration can be followed on the basis of a proper analytical model or an N-body computer simulation. The BARCODE formalism includes an implicit treatment for redshift space distortions. This enables a direct reconstruction on the basis of observational data, without the need for a correction of redshift space artifacts. In this contribution we provide a general overview of the the Cosmic Web connection with clusters and a description of the Bayesian BARCODE formalism. We conclude with a presentation of its successful workings with respect to test runs based on a simulated large scale matter distribution, in physical space as well as in redshift space.

The spatial distribution of the satellite populations of the Milky Way and Andromeda are puzzling in that they are nearly perpendicular to the disks of their central galaxies. To understand the origin of such configurations we study the alignment of the central galaxy, satellite system and dark matter halo in the largest of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation. We find that centrals and their satellite systems tend to be well aligned with their haloes, with a median misalignment angle of $33^{\circ}$ in both cases. While the centrals are better aligned with the inner $10$ kpc halo, the satellite systems are better aligned with the entire halo indicating that satellites preferentially trace the outer halo. The central - satellite alignment is weak (median misalignment angle of $52^{\circ}$) and we find that around $20\%$ of systems have a misalignment angle larger than $78^{\circ}$, which is the value for the Milky Way. The central - satellite alignment is a consequence of the tendency of both components to align with the dark matter halo. As a consequence, when the central is parallel to the satellite system, it also tends to be parallel to the halo. In contrast, if the central is perpendicular to the satellite system, as in the case of the Milky Way and Andromeda, then the central - halo alignment is much weaker. Dispersion-dominated (spheroidal) centrals have a stronger alignment with both their halo and their satellites than rotation-dominated (disk) centrals. We also found that the halo, the central galaxy and the satellite system tend to be aligned with the surrounding large-scale distribution of matter, with the halo being the better aligned of the three.

We introduce a new method for stacking voids and deriving their profile that greatly increases the potential of voids as a tool for precision cosmology. Given that voids are highly non-spherical and have most of their mass at their edge, voids are better described relative to their boundary rather than relative to their centre, as in the conventional spherical stacking approach. The boundary profile is obtained by computing the distance of each volume element from the void boundary. Voids can then be stacked and their profiles computed as a function of this boundary distance. This approach enhances the weak lensing signal of voids, both shear and convergence, by a factor of two when compared to the spherical stacking method. It also results in steeper void density profiles that are characterised by a very slow rise inside the void and a pronounced density ridge at the void boundary. The resulting boundary density profile is self-similar when rescaled by the thickness of the density ridge, implying that the average rescaled profile is independent of void size. The boundary velocity profile is characterized by outflows in the inner regions whose amplitude scales with void size, and by a strong inflow into the filaments and walls delimiting the void. This new picture enables a straightforward discrimination between collapsing and expanding voids both for individual objects as well as for stacked samples.

We introduce Copernicus Complexio (COCO), a high-resolution cosmological N-body simulation of structure formation in the $\Lambda{\rm CDM}{}$ model. COCO follows an approximately spherical region of radius $\sim 17.4h^{-1}\,{\rm Mpc}$ embedded in a much larger periodic cube that is followed at lower resolution. The high resolution volume has a particle mass of $1.135\times10^5h^{-1}{\rm M}_{\odot}$ (60 times higher than the Millennium-II simulation). COCO gives the dark matter halo mass function over eight orders of magnitude in halo mass; it forms $\sim 60$ haloes of galactic size, each resolved with about 10 million particles. We confirm the power-law character of the subhalo mass function, $\bar{N}(>\mu)\propto\mu^{-s}$, down to a reduced subhalo mass $M_{sub}/M_{200}\equiv\mu=10^{-6}$, with a best-fit power-law index, $s=0.94$, for hosts of mass $\langle M_{200}\rangle=10^{12}h^{-1}{\rm M}_{\odot}$. The concentration-mass relation of COCO haloes deviates from a single power law for masses $M_{200}<\textrm{a few}\times 10^{8}h^{-1}{\rm M}_{\odot}$, where it flattens, in agreement with results by Sanchez-Conde et al. The host mass invariance of the reduced maximum circular velocity function of subhaloes, $\nu\equiv V_{max}/V_{200}$, hinted at in previous simulations, is clearly demonstrated over five orders of magnitude in host mass. Similarly, we find that the average, normalised radial distribution of subhaloes is approximately universal (i.e. independent of subhalo mass), as previously suggested by the Aquarius simulations of individual haloes. Finally, we find that at fixed physical subhalo size, subhaloes in lower mass hosts typically have lower central densities than those in higher mass hosts.

We study lensing by voids in Cubic Galileon and Nonlocal gravity cosmologies, which are examples of theories of gravity that modify the lensing potential. We find voids in the dark matter and halo density fields of N-body simulations and compute their lensing signal analytically from the void density profiles, which we show are well fit by a simple analytical formula. In the Cubic Galileon model, the modifications to gravity inside voids are not screened and they approximately double the size of the lensing effects compared to GR. The difference is largely determined by the direct effects of the fifth force on lensing and less so by the modified density profiles. For this model, we also discuss the subtle impact on the force and lensing calculations caused by the screening effects of haloes that exist in and around voids. In the Nonlocal model, the impact of the modified density profiles and the direct modifications to lensing are comparable, but they boost the lensing signal by only $\approx 10\%$, compared with that of GR. Overall, our results suggest that lensing by voids is a promising tool to test models of gravity that modify lensing.

The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($\Lambda$CDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $\Lambda$CDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high resolution cosmological simulations. We find that planar structures are very common, and that ~10% of $\Lambda$CDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and $\Lambda$CDM predictions. In fact, ~10% of $\Lambda$CDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively.

We investigate the angular and kinematic distributions of satellite galaxies around a large sample of bright isolated primaries in the spectroscopic and photometric catalogues of the Sloan Digital Sky Survey (SDSS). We detect significant anisotropy in the spatial distribution of satellites. To test whether this anisotropy could be related to the rotating disks of satellites recently found by Ibata et al. in a sample of SDSS galaxies, we repeat and extend their analysis. Ibata et al. found an excess of satellites on opposite sides of their primaries having anticorrelated radial velocities. We find that this excess is sensitive to small changes in the sample selection criteria which can greatly reduce its significance. In addition, we find no evidence for correspondingly correlated velocities for satellites observed on the same side of their primaries, which would be expected for rotating disks of satellites. We conclude that the detection of coherent rotation in the satellite population in current observational samples is not robust. We compare our data to the $\Lambda$CDM Millennium simulations populated with galaxies according to the semi-analytic model of Guo et al. We find excellent agreement with the spatial distribution of satellites in the SDSS data and the lack of a strong signal from coherent rotation.

We investigate the characteristics and the time evolution of the cosmic web from redshift, z=2, to present time, within the framework of the NEXUS+ algorithm. This necessitates the introduction of new analysis tools optimally suited to describe the very intricate and hierarchical pattern that is the cosmic web. In particular, we characterize filaments (walls) in terms of their linear (surface) mass density. This is very good in capturing the evolution of these structures. At early times the cosmos is dominated by tenuous filaments and sheets, which, during subsequent evolution, merge together, such that the present day web is dominated by fewer, but much more massive, structures. We also show that voids are more naturally described in terms of their boundaries and not their centres. We illustrate this for void density profiles, which, when expressed as a function of the distance from void boundary, show a universal profile in good qualitative agreement with the theoretical shell-crossing framework of expanding underdense regions.

We study the formation and evolution of the cosmic web, using the high-resolution CosmoGrid \$\Lambda\$CDM simulation. In particular, we investigate the evolution of the large-scale structure around void halo groups, and compare this to observations of the VGS-31 galaxy group, which consists of three interacting galaxies inside a large void. The structure around such haloes shows a great deal of tenuous structure, with most of such systems being embedded in intra-void filaments and walls. We use the Nexus+ algorithm to detect walls and filaments in CosmoGrid, and find them to be present and detectable at every scale. The void regions embed tenuous walls, which in turn embed tenuous filaments. We hypothesize that the void galaxy group of VGS-31 formed in such an environment.

We study the substructure population of Milky Way (MW)-mass halos in the $\Lambda$CDM cosmology using a novel procedure to extrapolate subhalo number statistics beyond the resolution limit of N-body simulations. The technique recovers the mean and the variance of the subhalo abundance, but not its spatial distribution. It extends the dynamic range over which precise statistical predictions can be made by the equivalent of performing a simulation with 50 times higher resolution, at no additional computational cost. We apply this technique to MW-mass halos, but it can easily be applied to halos of any mass. We find up to $20\%$ more substructures in MW-mass halos than found in previous studies. Our analysis lowers the mass of the MW halo required to accommodate the observation that the MW has only three satellites with a maximum circular velocity $V_{max}\ge30 km/s$ in the $\Lambda$CDM cosmology. The probability of having a subhalo population similar to that in the MW is $20\%$ for a virial mass, $M_{200}=1\times10^{12} M_\odot$ and practically zero for halos more massive than $M_{200}=2\times10^{12} M_\odot$.

We use the distribution of maximum circular velocities, $V_{max}$, of satellites in the Milky Way (MW) to constrain the virial mass, $M_{200}$, of the Galactic halo under an assumed prior of a $\Lambda$CDM universe. This is done by analysing the subhalo populations of a large sample of halos found in the Millennium II cosmological simulation. The observation that the MW has at most three subhalos with $V_{max}\ge30 km/s$ requires a halo mass $M_{200}\le1.4\times10^{12} M_\odot$, while the existence of the Magellanic Clouds (assumed to have $V_{max}\ge60 km/s$) requires $M_{200}\ge1.0\times10^{12} M_\odot$. The first of these conditions is necessary to avoid the "too-big-to-fail" problem highlighted by Boylan-Kolchin et al., while the second stems from the observation that massive satellites like the Magellanic Clouds are rare. When combining both requirements, we find that the MW halo mass must lie in the range $0.25 \le M_{200}/(10^{12} M_\odot) \le 1.4$ at $90\%$ confidence. The gap in the abundance of Galactic satellites between $30 km/s\le V_{max} \le 60 km/s$ places our galaxy in the tail of the expected satellite distribution.

The cosmic web is the largest scale manifestation of the anisotropic gravitational collapse of matter. It represents the transitional stage between linear and non-linear structures and contains easily accessible information about the early phases of structure formation processes. Here we investigate the characteristics and the time evolution of morphological components since. Our analysis involves the application of the NEXUS Multiscale Morphology Filter (MMF) technique, predominantly its NEXUS+ version, to high resolution and large volume cosmological simulations. We quantify the cosmic web components in terms of their mass and volume content, their density distribution and halo populations. We employ new analysis techniques to determine the spatial extent of filaments and sheets, like their total length and local width. This analysis identifies cluster and filaments as the most prominent components of the web. In contrast, while voids and sheets take most of the volume, they correspond to underdense environments and are devoid of group-sized and more massive haloes. At early times the cosmos is dominated by tenuous filaments and sheets, which, during subsequent evolution, merge together, such that the present day web is dominated by fewer, but much more massive, structures. The analysis of the mass transport between environments clearly shows how matter flows from voids into walls, and then via filaments into cluster regions, which form the nodes of the cosmic web. We also study the properties of individual filamentary branches, to find long, almost straight, filaments extending to distances larger than 100Mpc/h. These constitute the bridges between massive clusters, which seem to form along approximatively straight lines.

We investigate how long-range scalar interactions affect the properties of dark matter haloes. For doing so we employ the ReBEL model which implements an additional interaction between dark matter particles. On the phenomenological level this is equivalent to a modification of gravity. We analyse the differences between five ReBEL models and $\Lambda$CDM using a series of high resolution cosmological simulations. Emphasis is placed on investigating how halo properties change in the presence of a fifth force. We report that the density profile of ReBEL haloes is well described by the NFW profile but with mean concentrations from $5\%$ to a few times higher than the standard $\Lambda$CDM value. We also find a slight increase of the halo spin for haloes more massive than $5\times10^{11}\M_{\odot}$, reflecting a higher rotational support of those haloes due to scalar forces. In addition, the dark matter haloes in our models are more spherical than their counterparts in $\Lambda$CDM. The ReBEL haloes are also more virialised, with a large difference from $\Lambda$CDM for strong fifth forces and a much smaller change for weak scalar interactions.

We study the formation and evolution of filamentary configurations of dark matter haloes in voids. Our investigation uses the high-resolution LambdaCDM simulation CosmoGrid to look for void systems resembling the VGS_31 elongated system of three interacting galaxies that was recently discovered by the Void Galaxy Survey (VGS) inside a large void in the SDSS galaxy redshift survey. HI data revealed these galaxies to be embedded in a common elongated envelope, possibly embedded in intravoid filament. In the CosmoGrid simulation we look for systems similar to VGS_31 in mass, size and environment. We find a total of eight such systems. For these systems, we study the distribution of neighbour haloes, the assembly and evolution of the main haloes and the dynamical evolution of the haloes, as well as the evolution of the large-scale structure in which the systems are embedded. The spatial distribution of the haloes follows that of the dark matter environment. We find that VGS_31-like systems have a large variation in formation time, having formed between 10 Gyr ago and the present epoch. However, the environments in which the systems are embedded evolved resemble each other substantially. Each of the VGS_31-like systems is embedded in an intra-void wall, that no later than z = 0.5 became the only prominent feature in its environment. While part of the void walls retain a rather featureless character, we find that around half of them are marked by a pronounced and rapidly evolving substructure. Five haloes find themselves in a tenuous filament of a few Mpc/h long inside the intra-void wall. Finally, we compare the results to observed data from VGS_31. Our study implies that the VGS_31 galaxies formed in the same (proto)filament, and did not meet just recently. The diversity amongst the simulated halo systems indicates that VGS_31 may not be typical for groups of galaxies in voids.

One of the important unknowns of current cosmology concerns the effects of the large scale distribution of matter on the formation and evolution of dark matter haloes and galaxies. One main difficulty in answering this question lies in the absence of a robust and natural way of identifying the large scale environments and their characteristics. This work summarizes the NEXUS+ formalism which extends and improves our multiscale scale-space MMF method. The new algorithm is very successful in tracing the Cosmic Web components, mainly due to its novel filtering of the density in logarithmic space. The method, due to its multiscale and hierarchical character, has the advantage of detecting all the cosmic structures, either prominent or tenuous, without preference for a certain size or shape. The resulting filamentary and wall networks can easily be characterized by their direction, thickness, mass density and density profile. These additional environmental properties allows to us to investigate not only the effect of environment on haloes, but also how it correlates with the environment characteristics.

We introduce the NEXUS algorithm for the identification of Cosmic Web environments: clusters, filaments, walls and voids. This is a multiscale and automatic morphological analysis tool that identifies all the cosmic structures in a scale free way, without preference for a certain size or shape. We develop the NEXUS method to incorporate the density, tidal field, velocity divergence and velocity shear as tracers of the Cosmic Web. We also present the NEXUS+ procedure which, taking advantage of a novel filtering of the density in logarithmic space, is very successful at identifying the filament and wall environments in a robust and natural way. To asses the algorithms we apply them to an N-body simulation. We find that all methods correctly identify the most prominent filaments and walls, while there are differences in the detection of the more tenuous structures. In general, the structures traced by the density and tidal fields are clumpier and more rugged than those present in the velocity divergence and velocity shear fields. We find that the NEXUS+ method captures much better the filamentary and wall networks and is successful in detecting even the fainter structures. We also confirm the efficiency of our methods by examining the dark matter particle and halo distributions.