We introduce an algorithm aimed to identify large-scale filaments founded on the conception that these structures are bridges of matter that connect high density peaks. Our method is based on two standard tools, the Minimal Spanning Tree (MST) and the Friends of Friends (FoF) algorithm. Briefly, the complete process consists of five stages. Initially, we use the FoF algorithm to select intermediate density regions to stave off underdense zones (voids). Next, we build the MST by restricting it only to the regions defined in the previous step. Then the tree is pruned according to the length of its branches keeping the most dominant and discarding the more tenuous ones. Finally, the filaments are individualised according to the mass of the ends and smoothed using a B-spline fitting routine. To assess the results of the filament finder we apply it to a cosmological simulation. By focusing our analysis on those filaments whose halos at the ends have large masses, we found that the radial density profile, at scales around $1\, h^{-1} \mathrm{Mpc}$, follow a power law function with index -2. Even though the method only relies on halo positions, it is capable to recover the expected velocity field in filamentary structures. Large infall velocities coming from low density environment approach perpendicularly to the filaments and diverges toward to the ends. By studying the transverse velocity dispersion, we estimate the dynamical linear density following Eisenstein et al. (1997), finding a good correspondence between that and the actual mass per unit length of the filaments.

The tidal torque theory (TTT) relates the origin and evolution of angular momentum with the environment in which dark matter (DM) haloes form. The deviations introduced by late non-linearities are commonly thought as noise in the model. In this work, we analyze a cosmological simulation looking for systematics on these deviations, finding that the classification of DM haloes according to their angular momentum growth results in samples with different internal alignment, spin parameter distribution and assembly history. Based on this classification, we obtain that low mass haloes are embedded in denser environments if they have acquired angular momentum below the TTT expectations (L haloes), whereas at high masses enhanced clustering is typically associated with higher angular momentum growths (W haloes). Additionally, we find that the low mass signal has a weak dependence on the direction, whereas the high mass signal is entirely due to the structure perpendicular to the angular momentum. Finally, we study the anisotropy of the matter distribution around haloes as a function of their mass. We find that the angular momentum direction of W (L) haloes remains statistically perpendicular (parallel) to the surrounding structure across the mass range $11<\mathrm{log}(M/h^{-1}\mathrm{M}_{\odot})<14$, whereas haloes following TTT show a "spin flip" mass consistent with previously reported values ($\sim 5 \times 10^{12}$ $h^{-1}\mathrm{M}_\odot$). Hence, whether the spin flip mass of the deviated samples is highly shifted or straightly undefined, our results indicate that is remarkably connected to the haloes angular momentum growth.

We present a new cosmological test using the distribution of galaxies around cosmic voids without assuming a fiducial cosmology. The test is based on a physical model for the void-galaxy cross-correlation function projected along and perpendicular to the line of sight. We treat correlations in terms of void-centric angular distances and redshift differences between void-galaxy pairs, hence it is not necessary to assume a fiducial cosmology. This model reproduces the coupled dynamical (Kaiser effect, RSD) and geometrical (Alcock-Paczynski effect, GD) distortions that affect the correlation measurements. It also takes into account the scale mixing due to the projection ranges in both directions. The model is general, so it can be applied to an arbitrary cylindrical binning scheme, not only in the case of the projected correlations. It primarily depends on two cosmological parameters: $\Omega_m$, the matter fraction of the Universe today (sensitive to GD), and $\beta$, the ratio between the growth rate factor of density perturbations and the tracer bias (sensitive to RSD). In the context of the new generation of galaxy spectroscopic surveys, we calibrated the test using the Millennium XXL simulation for different redshifts. The method successfully recovers the cosmological parameters. We studied the effect of measuring with different projection ranges, finding robust results up to wide ranges. The resulting data covariance matrices are relatively small, which reduces the noise in the Gaussian likelihood analysis and will allow the usage of a smaller number of mock catalogues. The performance evaluated in this work indicates that the developed method is a promising test to be applied on real data.

The expanding complex pattern of filaments, walls and voids build the evolving cosmic web with material flowing from underdense onto high density regions. Here we explore the dynamical behaviour of voids and galaxies in void shells relative to neighboring overdense superstructures, using the Millenium Simulation and the main galaxy catalogue in Sloan Digital Sky Survey data. We define a correlation measure to estimate the tendency of voids to be located at a given distance from a superstructure. We find voids-in-clouds (S-types) preferentially located closer to superstructures than voids-in-voids (R-types) although we obtain that voids within $\sim40~\mathrm{Mpc}\,\mathrm{h}^{-1}$ of superstructures are infalling in a similar fashion independently of void type. Galaxies residing in void shells show infall towards the closest superstructure, along with the void global motion, with a differential velocity component depending on their relative position in the shell with respect to the direction to the superstructure. This effect is produced by void expansion and therefore is stronger for R-types. We also find that galaxies in void shells facing the superstrucure flow towards the overdensities faster than galaxies elsewere at the same relative distance to the superstructure. The results obtained for the simulation are also reproduced for the SDSS data with a linearized velocity field implementation.

We analyse the clustering of cosmic voids using a numerical simulation and the main galaxy sample from the Sloan Digital Sky Survey. We take into account the classification of voids into two types that resemble different evolutionary modes: those with a rising integrated density profile (void-in-void mode, or R-type) and voids with shells (void-in-cloud mode, or S-type). The results show that voids of the same type have stronger clustering than the full sample. We use the correlation analysis to define void clumps, associations with at least two voids separated by a distance of at most the mean void separation. In order to study the spatial configuration of void clumps, we compute the minimal spanning tree and analyse their multiplicity, maximum length and elongation parameter. We further study the dynamics of the smaller sphere that encloses all the voids in each clump. Although the global densities of void clumps are different according to their member-void types, the bulk motions of these spheres are remarkably lower than those of randomly placed spheres with the same radii distribution. In addition, the coherence of pairwise void motions does not strongly depend on whether voids belong to the same clump. Void clumps are useful to analyse the large-scale flows around voids, since voids embedded in large underdense regions are mostly in the void-in-void regime, were the expansion of the larger region produces the separation of voids. Similarly, voids around overdense regions form clumps that are in collapse, as reflected in the relative velocities of voids that are mostly approaching.

We perform a statistical study of the global motion of cosmic voids using both a numerical simulation and observational data. We analyse their relation to large-scale mass flows and the physical effects that drive those motions. We analyse the bulk motions of voids, defined by the mean velocity of haloes in the surrounding shells in the numerical simulation, and by galaxies in the Sloan Digital Sky Survey Data Release 7. We find void mean bulk velocities close to 400 km/s, comparable to those of haloes (~ 500-600 km/s), depending on void size and the large-scale environment. Statistically, small voids move faster than large ones, and voids in relatively higher density environments have higher bulk velocities than those placed in large underdense regions. Also, we analyze the mean mass density around voids finding, as expected, large-scale overdensities (underdensities) along (opposite to) the void motion direction, suggesting that void motions respond to a pull-push mechanism. This contrasts with massive cluster motions who are mainly governed by the pull of the large-scale overdense regions. Our analysis of void pairwise velocities shows how their relative motions are generated by large-scale density fluctuations. In agreement with linear theory, voids embedded in low (high) density regions mutually recede (attract) each other, providing the general mechanism to understand the bimodal behavior of void motions. In order to compare the theoretical results and the observations we have inferred void motions in the SDSS using linear theory, finding that the estimated observational void motions are in statisticalagreement with the results of the simulation. Regarding large-scale flows, our results suggest a scenario of galaxies and galaxy systems flowing away from void centers with the additional, and morerelevant, contribution of the void bulk motion to the total velocity.

Merging galaxy systems provides observational evidence of the existence of dark matter and constraints on its properties. Therefore, statistical uniform samples of merging systems would be a powerful tool for several studies. In this work we presents a new methodology for merging systems identification and the results of its application to galaxy redshift surveys. We use as starting point a mock catalogue of galaxy systems, identified using traditional FoF algorithms, which experienced a major merger as indicated by its merger tree. Applying machine learning techniques in this training sample, and using several features computed from the observable properties of galaxy members, it is possible to select galaxy groups with a high probability of have been experienced a major merger. Next we apply clustering techniques on galaxy members in order to reconstruct the properties of the haloes involved in such merger. This methodology provides a highly reliable sample of merging systems with low contamination and precise recovered properties. We apply our techniques in samples of galaxy systems obtained from SDSS-DR7, WINGS and HeCS. Our results recover previously known merging systems and provide several new candidates. We present its measured properties and discuss future analysis on current and forthcoming samples.

We compute the bulk motions of cosmic voids, using a $\Lambda$CDM numerical simulation considering the mean velocities of the dark matter inside the void itself and that of the haloes in the surrounding shell. We find coincident values of these two measures in the range $\sim$ 300-400 km/s, not far from the expected mean peculiar velocities of groups and galaxy clusters. When analysing the distribution of the pairwise relative velocities of voids, we find a remarkable bimodal behaviour consistent with an excess of both systematically approaching and receding voids. We determine that the origin of this bimodality resides in the void large scale environment, since once voids are classified into void-in-void (R-type) or void-in-cloud (S-type), R-types are found mutually receding away, while S-types approach each other. The magnitude of these systematic relative velocities account for more than 100 km/s, reaching large coherence lengths of up to 200 h$^{-1}$ Mpc . We have used samples of voids from the Sloan Digital Sky Survey Data Release 7 (SDSS-DR7) and the peculiar velocity field inferred from linear theory, finding fully consistent results with the simulation predictions. Thus, their relative motion suggests a scenario of a sparkling Universe, with approaching and receding voids according to their local environment.

The estimation of cosmological constraints from observations of the large scale structure of the Universe, such as the power spectrum or the correlation function, requires the knowledge of the inverse of the associated covariance matrix, namely the precision matrix, $\mathbf{\Psi}$. In most analyses, $\mathbf{\Psi}$ is estimated from a limited set of mock catalogues. Depending on how many mocks are used, this estimation has an associated error which must be propagated into the final cosmological constraints. For future surveys such as Euclid and DESI, the control of this additional uncertainty requires a prohibitively large number of mock catalogues. In this work we test a novel technique for the estimation of the precision matrix, the covariance tapering method, in the context of baryon acoustic oscillation measurements. Even though this technique was originally devised as a way to speed up maximum likelihood estimations, our results show that it also reduces the impact of noisy precision matrix estimates on the derived confidence intervals, without introducing biases on the target parameters. The application of this technique can help future surveys to reach their true constraining power using a significantly smaller number of mock catalogues.

We study the dependence of the properties of group galaxies on the surrounding large-scale environment, using SDSS-DR7 data. Galaxies are ranked according to their luminosity within each group and classified morphologically by the S\'ersic index. We have considered samples of the host groups in superstructures of galaxies, and elsewhere. We find a significant dependence of the properties of late-type brightest group galaxies on the large-scale environment: they show statistically significant higher luminosities and stellar masses, redder u-r colours, lower star formation activity and longer star-formation time-scale when embedded in superstructures. By contrast, the properties of the early-type brightest group galaxies are remarkably similar regardless of the group global environment. The other group member galaxies exhibit only the local influence of the group they inhabit. Our analysis comprises tests against the dependence on the host group luminosity and we argue that group brightest member properties are not only determined by the host halo, but also by the large-scale structure which can influence the accretion process onto their late-type brightest galaxies.

Inspired on the well known dynamical dichotomy predicted in voids, where some underdense regions expand whereas others collapse due to overdense surrounding regions, we explored the interplay between the void inner dynamics and its large scale environment. The environment is classified depending on its density as in previous works. We analyse the dynamical properties of void-centered spherical shells at different void-centric distances depending on this classification. The above dynamical properties are given by the angular distribution of the radial velocity field, its smoothness, the field dependence on the tracer density and shape, and the field departures from linear theory. We found that the velocity field in expanding voids follows more closely the linear prediction, with a more smooth velocity field. However when using velocity tracers with large densities such deviations increase. Voids with sizes around $18\,h^{-1}\,Mpc$ are in a transition regime between regions with expansion overpredicted and underpredicted from linear theory. We also found that velocity smoothness increases as the void radius, indicating the laminar flow dominates the expansion of larger voids (more than $18\,h^{-1}\,Mpc$). The correlations observed suggest that nonlinear dynamics of the inner regions of voids could be dependent on the evolution of the surrounding structures. These also indicate possible scale couplings between the void inner expansion and the large scale regions where voids are embedded. These results shed some light to the origin of nonlinearities in voids, going beyond the fact that voids just quickly becomes nonlinear as they become emptier.

Cosmic voids are becoming key players in testing the physics of our Universe. Here we concentrate on the abundances and the dynamics of voids as these are among the best candidates to provide information on cosmological parameters. Cai, Padilla \& Li (2014) use the abundance of voids to tell apart Hu \& Sawicki $f(R)$ models from General Relativity. An interesting result is that even though, as expected, voids in the dark matter field are emptier in $f(R)$ gravity due to the fifth force expelling away from the void centres, this result is reversed when haloes are used to find voids. The abundance of voids in this case becomes even lower in $f(R)$ compared to GR for large voids. Still, the differences are significant and this provides a way to tell apart these models. The velocity field differences between $f(R)$ and GR, on the other hand, are the same for halo voids and for dark matter voids. Paz et al. (2013), concentrate on the velocity profiles around voids. First they show the necessity of four parameters to describe the density profiles around voids given two distinct void populations, voids-in-voids and voids-in-clouds. This profile is used to predict peculiar velocities around voids, and the combination of the latter with void density profiles allows the construction of model void-galaxy cross-correlation functions with redshift space distortions. When these models are tuned to fit the measured correlation functions for voids and galaxies in the Sloan Digital Sky Survey, small voids are found to be of the void-in-cloud type, whereas larger ones are consistent with being void-in-void. This is a novel result that is obtained directly from redshift space data around voids. These profiles can be used to remove systematics on void-galaxy Alcock-Pacinsky tests coming from redshift-space distortions.

Using the redshift-space distortions of void-galaxy cross-correlation function we analyse the dynamics of voids embedded in different environments. We compute the void-galaxy cross-correlation function in the Sloan Digital Sky Survey (SDSS) in terms of distances taken along the line of sight and projected into the sky. We analyse the distortions on the cross-correlation isodensity levels and we find anisotropic isocontours consistent with expansion for large voids with smoothly rising density profiles and collapse for small voids with overdense shells surrounding them. Based on the linear approach of gravitational collapse theory we developed a parametric model of the void-galaxy redshift space cross-correlation function. We show that this model can be used to successfully recover the underlying velocity and density profiles of voids from redshift space samples. By applying this technique to real data, we confirm the twofold nature of void dynamics: large voids typically are in an expansion phase whereas small voids tend to be surrounded by overdense and collapsing regions. These results are obtained from the SDSS spectroscopic galaxy catalogue and also from semi-analytic mock galaxy catalogues, thus supporting the viability of the standard LCDM model to reproduce large scale structure and dynamics.

We perform a statistical study focused on void environments. We examine galaxy density profiles around voids in the SDSS, finding a correlation between void--centric distance to the shell of maximum density and void radius when a maximum in overdensity exists. We analyze voids with and without a surrounding over-dense shell in the SDSS. We find that small voids are more frequently surrounded by over-dense shells whereas the radial galaxy density profile of large voids tends to rise smoothly towards the mean galaxy density. We analyse the fraction of voids surrounded by overdense shells finding a continuous trend with void radius. The differences between voids with and without an overdense shell around them can be understood in terms of whether the voids are, on average, in the process of collapsing or continuing their expansion, respectively, in agreement with previous theoretical expectations. We use numerical simulations coupled to semi-analytic models of galaxy formation in order to test and interpret our results. The very good agreement between the mock catalog results and the observations provides additional support to the viability of a LCDM model to reproduce the large scale structure of the universe as defined by the void network, in a way which has not been analysed previously.

In the present work, we extend the classic halo model for the large-scale matter distribution including a triaxial model for the halo profiles and their alignments. In particular, we derive general expressions for the halo-matter cross correlation function. In addition, by numerical integration, we obtain instances of the cross-correlation function depending on the directions given by halo shape axes. These functions are called anisotropic cross-correlations. With the aim of comparing our theoretical results with the simulations, we compute averaged anisotropic correlations in cones with their symmetry axis along each shape direction of the centre halo. From these comparisons we characterise and quantify the alignment of dark matter haloes on the $\Lambda$CDM context by means of the presented anisotropic halo model. As our model requires multidimensional integral computation we implement a Monte Carlo method on GPU hardware which allows us to increase the precision of the results whereas it improves the performance of the computation.

We analyse properties of galaxy groups and their dependence on the large scale environment as defined by superstructures. We find that group galaxy cross correlations depend only on group properties regardless the groups reside in superstructures. This indicates that the total galaxy density profile around groups is independent of the global environment. At a given global luminosity, a proxy to group total mass, groups have a larger stellar mass content by a factor 1.3, a relative excess independent of the group luminosity. Groups in superstructures have 40 per cent higher velocity dispersions and systematically larger minimal enclosing radii. We also find that the stellar population of galaxies in groups in superstructures is systematically older as infered from the galaxy spectra Dn4000 parameter. Although the galaxy number density profile of groups is independent of environment, the star formation rate and stellar mass profile of the groups residing in superstructures differs from groups elsewhere. For groups residing in superstructures, the combination of a larger stellar mass content and star formation rate produces a larger time scale for star formation regardless the distance to the group center. Our results provide evidence that groups in superstructures formed earlier than elsewhere, as expected in the assembly bias scenario.

We analyse the dependence of clustering properties of galaxies as a function of their large-scale environment. In order to characterize the environment on large scales, we use the catalogue of future virialized superstructures (FVS) by Luparello et al. and separate samples of luminous galaxies according to whether or not they belong to FVS. In order to avoid biases in the selection of galaxies, we have constructed different subsamples so that the distributions of luminosities and masses are comparable outside and within FVS. As expected, at large scales, there is a strong difference between the clustering of galaxies inside and outside FVS. However, this behaviour changes at scales r < 1 Mpc/h, where the correlations have similar amplitudes. The amplitude of the two-halo term of the correlation function for objects inside FVS does not depend on their mass, but rather on that of the FVS. This is confirmed by comparing this amplitude with that expected from extended Press-Schechter fits. In order to compare these observational results with current models for structure formation, we have performed a similar analysis using a semi-analytic implementation in a LCDM cosmological model. We find that the cross-correlation functions from the mock catalogue depend on the large-scale structures in a similar way to the observations. From our analysis, we conclude that the clustering of galaxies within the typical virialized regions of groups, mainly depends on the halo mass, irrespective of the large-scale environment.

In this paper we analyse the alignment of galaxy groups with the surrounding large scale structure traced by spectroscopic galaxies from the Sloan Digital Sky Survey Data Release 7. We characterise these alignments by means of an extension of the classical two-point cross-correlation function, developed by Paz et al. 2008 (arXiv:0804.4477, MNRAS 389 1127). We find a strong alignment signal between the projected major axis of group shapes and the surrounding galaxy distribution up to scales of 30 Mpc/h. This observed anisotropy signal becomes larger as the galaxy group mass increases, in excellent agreement with the corresponding predicted alignment obtained from mock catalogues and LCDM cosmological simulations. These measurements provide new direct evidence of the adequacy of the gravitational instability picture to describe the large-scale structure formation of our Universe.

We perform a statistical analysis of the peculiar velocity field around dark matter haloes in numerical simulations. We examine different properties of the infall of material onto haloes and its relation to central halo shapes and the shape of the large scale surrounding regions (LSSR). We find that the amplitude of the infall velocity field along the halo shape minor axis is larger than that along the major axis. This is consistent for general triaxial haloes, and for both prolate and oblate systems. We also report a strong anisotropy of the velocity field along the principal axes of the LSSR. The infall velocity field around dark matter haloes reaches a maximum value along the direction of the minor axis of the LSSR, whereas along the direction of its major axis, it exhibits the smallest velocities. We also analyse the dependence of the matter velocity field on the local environment. The amplitude of the infall velocity at high local density regions is larger than at low local density regions. The velocity field tends to be more laminar along the direction towards the minor axis of the LSSR, where the mean ratio between flow velocity and velocity dispersion is of order unity and nearly constant up to scales of 15 Mpc/h. We also detect anisotropies in the outflowing component of the velocity field, showing a maximum amplitude along the surrounding LSSR major axis.

We study the alignments between the angular momentum of individual objects and the large-scale structure in cosmological numerical simulations and real data from the Sloan Digital Sky Survey, Data Release 6. To this end we measure anisotropies in the two point cross-correlation function around simulated halos and observed galaxies, studying separately the 1- and 2-halo regimes. The alignment of the angular momentum of dark-matter haloes in LCDM simulations is found to be dependent on scale and halo mass. At large distances (2-halo regime), the spins of high mass haloes are preferentially oriented in the direction perpendicular to the distribution of matter; lower mass systems show a weaker trend that may even reverse to show an angular momentum in the plane of the matter distribution. In the 1-halo term regime, the angular momentum is aligned in the direction perpendicular to the matter distribution; the effect is stronger than for the 1-halo term and increases for higher mass systems. On the observational side, we focus our study on galaxies in the Sloan Digital Sky Survey, Data Release 6 (SDSS-DR6) with elongated apparent shapes, and study alignments with respect to the major semi-axis. We find an excess of structure in the direction of the major semi-axis for all samples; the red sample shows the highest alignment (2.7+-0.08%) and indicates that the angular momentum of flattened spheroidals tends to be perpendicular to the large-scale structure. (Abridged)

We study the properties of the 3-dimensional and projected shapes of haloes using high resolution numerical simulations and observational data where the latter comes from the 2PIGG (Eke et al. 2004) and SDSS-DR3GC group catalogues (Merchan & Zandivarez 2005). We investigate the dependence of halo shape on characteristics such as mass and number of members. In the 3-dimensional case, we find a significant correlation between the mass and halo shape; massive systems are more prolate than small haloes. We detect a source of strong systematics in estimates of the triaxiality of a halo, which is found to be a strong function of the number of members; LCDM haloes usually characterised by triaxial shapes, slightly bent toward prolate forms, appear more oblate when taking only a small subset of the halo particles. The ellipticities of observed 2PIGG and SDSS-DR3GC groups are found to be strongly dependent on the number of group members, so that poor groups appear more elongated than rich ones. However, this is again an artifact caused by poor statistics and not an intrinsic property of the galaxy groups, nor an effect from observational biases. We interpret these results with the aid of a GALFORM mock 2PIGG catalogue. When comparing the group ellipticities in mock and real catalogues, we find an excellent agreement between the trends of shapes with number of group members. When carefully taking into account the effects of low number statistics, we find that more massive groups are consistent with more elongated shapes. Finally, our studies find no significant correlations between the shape of observed 2PIGG or SDSS-DR3GC groups with the properties of galaxy members such as colour or spectral type index.